An official website of the United States government
The composition of the crust is one of the most uncertain and controversial components of continental estimates due to (1) limited direct access and (2) inconsistent indirect assessments. Here we show that by combining high-resolution shear velocity (Vs) models with newly measured with newly measured ratio of compressional wave velocity (Vp) and Vs, or Vp/Vs ratio, for the crystalline crust, a 3-D composition (SiO2 wt%) model of the continental crust can be derived with uncertainty estimates. Comparing the model with local xenolith data shows consistency at mid and lower crustal depths. The spatial patterns in bulk SiO2 content correlate with major geological provinces, including the footprints of Cenozoic and Mesozoic mafic volcanism in the western U.S., and offer new insight into the composition and evolution of the continental U.S.
more »
« less
A crustal thermal model of the conterminous United States constrained by multiple data sets: a Monte–Carlo approach
SUMMARY The thermal structure of the continental crust plays a critical role in understanding its elastic and rheologic properties as well as its dynamic processes. Thermal parameter data sets on continental scales have been used to constrain the crustal thermal structure, including both the direct (e.g. temperature, heat flux and heat conductivity measured at the surface) and indirect (e.g. seismically derived Mohorovičić discontinuity (Moho) temperature, geomagnetically derived Curie depth) observations. In this study, we present a new continental scale crustal heat generation model with additional information from seismologically inferred crustal composition. Together with previous direct and indirect thermal parameter data sets in the conterminous United States, we use the new crustal heat generation model to construct a 3-D crustal temperature model under a newly developed Bayesian framework. Specifically, we first derive profiles of crustal heat generation based on an empirical geochemical relationship at 1683 locations where seismologically derived crustal composition information is available. Then for each of these locations, the average heat generation values in the upper, middle and lower crust are combined with other thermal parameters through a Markov Chain Monte-Carlo inversion for a conductive, vertically smooth temperature profile. The results, posterior distributions of temperature profiles, are used to generate a 3-D crustal thermal model with the uncertainties systematically assessed. The new temperature model overall exhibits similar patterns to that from the U.S. Geological Survey National Crustal Model, but also reduces possible biases and the model's dependence on a single thermal parameter.
more »
« less
- PAR ID:
- 10583493
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 241
- Issue:
- 3
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 1712-1725
- Size(s):
- p. 1712-1725
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract This study provides a global assessment of the abundance of the major oxides in the deep continental crust. The combination of geochemistry and seismology better constrains the composition of the middle and lower continental crust better than either discipline can achieve alone. The inaccessible nature of the deep crust (typically >15 km) forces reliance on analog samples and modeling results to interpret its bulk composition, evolution, and physical properties. A common practice relates major oxide compositions of small‐ to medium‐scale samples (e.g., medium to high metamorphic grade terrains and xenoliths) to large scale measurements of seismic velocities (Vp, Vs, Vp/Vs) to determine the composition of the deep crust. We provide a framework for building crustal models with multidisciplinary constraints on composition. We present a global deep crustal model that documents compositional changes with depth and accounts for uncertainties in Moho depth, temperature, and physical and chemical properties. Our 3D compositional model of the deep crust uses the USGS Global Seismic Structure Catalog (Mooney, 2015) and a compilation of geochemical analyses on amphibolite and granulite facies lithologies (Sammon & McDonough, 2021,https://doi.org/10.1029/2021JB022791). We find a SiO2gradient from 61.2 ± 7.3 to 53.3 ± 4.8 wt.% from the middle to the base of the crust, with the equivalent lithological gradient ranging from quartz monzonite to gabbronorite. In addition, we calculate trace element abundances as a function of depth from their correlations with major oxides. From here, other lithospheric properties, such as Moho heat flux ( mW/m2), are derived.more » « less
-
Abstract The chemical composition of the deep continental crust is key to understanding the formation and evolution of the continental crust. Constraining the chemical composition of present‐day deep continental crust is, however, limited by indirect accessibility. This paper presents a modeling method for constraining deep crustal chemical structures from observed crustal seismic structures. We compiled a set of published composition models for the continental crust to construct functional relationships between seismic wave speed and major oxide content in the crust. Phase equilibria and compressional wave speeds (VP) for each composition model were calculated over a range of depths and temperatures of the deep crust. For conditions within the alpha(α)‐quartz stability field, robust functional relationships were obtained betweenVPand major oxide contents of the crust. Based on these relationships, observedVPof the deep crust can be inverted to chemical compositions for regions with given geotherms. We provide a MATLAB code for this process (CalcCrustComp). We apply this method to constrain compositions from deep crustalVPof global typical tectonic settings and the North China Craton (NCC). Our modeling results suggest that the lower crust in subduction‐related and rifting‐related tectonic settings may be more mafic than platforms/shields and orogens. The lowVPsignature in the deep crust of the NCC can be explained by intermediate crustal compositions, higher water contents, and/or higher temperatures. The chemical structure obtained by this method can serve as a reference model to further identify deep crustal features.more » « less
-
The formation of continental crust in magmatic arcs involves cooling of hot magmas to a relatively colder crust enhanced by exhumation and hydrothermal circulation in the upper crust. To quantify the influence of these processes on the thermal and rheological states of the crust, we developed a one-dimensional thermal evolution model, which invokes conductive cooling, advection of crust by erosion-driven exhumation, and cooling by hydrothermal circulation. We parameterized hydrothermal cooling by adopting depth-dependent effective thermal conductivity, which is determined by the crustal permeability structure and the prescribed Nusselt number at the surface. Different combinations of erosion rate and Nusselt number were tested to study the evolution of crustal geotherms, surface heat flux, and cooling rate. Simulations and scaling analyses quantify how erosion and hydrothermal circulation promote cooling via increasing total surface heat flux compared to pure conductive cooling. Hydrothermal circulation imposes intense short-term and persistent long-term cooling effects. Thinner, warmer, fast exhuming crust, with higher permeability and more vigorous hydrothermal circulation, leads to higher steady-state total surface heat flux. Hydrothermal cooling at steady state is more effective when the Péclet number is small. Hydrothermal cooling also changes crustal rheological state and thickens the brittle crust. This in turn promotes the initiation of brittle deformation in the upper crust in magmatic arcs or in regions undergoing exhumation. Interpretation of low-temperature thermochronological data could overestimate average cooling rates if hydrothermal cooling is not considered.more » « less
-
Abstract Accurately determining the seismic structure of the continental deep crust is crucial for understanding its geological evolution and continental dynamics in general. However, traditional tools such as surface waves often face challenges in solving the trade‐offs between elastic parameters and discontinuities. In this work, we present a new approach that combines two established inversion techniques, receiver function H‐κstacking and joint inversion of surface wave dispersion and receiver function waveforms, within a Bayesian Monte Carlo (MC) framework to address these challenges. Demonstrated by synthetic tests, the new method greatly reduces trade‐offs between critical parameters, such as the deep crustal Vs, Moho depth, and crustal Vp/Vs ratio. This eliminates the need for assumptions regarding crustal Vp/Vs ratios in joint inversion, leading to a more accurate outcome. Furthermore, it improves the precision of the upper mantle velocity structure by reducing its trade‐off with Moho depth. Additional notes on the sources of bias in the results are also included. Application of the new approach to USArray stations in the Northwestern US reveals consistency with previous studies and identifies new features. Notably, we find elevated Vp/Vs ratios in the crystalline crust of regions such as coastal Oregon, suggesting potential mafic composition or fluid presence. Shallower Moho depth in the Basin and Range indicates reduced crustal support to the elevation. The uppermost mantle Vs, averaging 5 km below Moho, aligns well with the Pn‐derived Moho temperature variations, offering the potential of using Vs as an additional constraint to Moho temperature and crustal thermal properties.more » « less
