Debate continues on the amount and distribution of radioactive heat producing elements (i.e., U, Th, and K) in the Earth, with estimates for mantle heat production varying by an order of magnitude. Constraints on the bulk‐silicate Earth's (BSE) radiogenic power also places constraints on overall BSE composition. Geoneutrino detection is a direct measure of the Earth's decay rate of Th and U. The geoneutrino signal has contributions from the local (
The composition of the lower continental crust is well studied but poorly understood because of the difficulty of sampling large portions of it. Petrological and geochemical analyses of this deepest portion of the continental crust are limited to the study of high‐grade metamorphic lithologies, such as granulite. In situ lower crustal studies require geophysical experiments to determine regional‐scale phenomena. Since geophysical properties, such as shear wave velocity (Vs), are nonunique among different compositions and temperatures, the most informative lower crustal models combine both geochemical and geophysical knowledge. We explored a combined modeling technique by analyzing the Basin and Range and Colorado Plateau of the United States, a region for which plentiful geochemical and geophysical data are available. By comparing seismic velocity predictions based on composition and thermodynamic principles to ambient noise inversions, we identified three compositional trends in the southwestern United States that reflect three different geologic settings. The Colorado Plateau (thick crust), Northern Basin and Range (medium crust), and Southern Basin and Range (thin crust) have intermediate, intermediate‐mafic, and mafic deep crustal compositions. Identifying the composition of the lower crust depends heavily on its temperature because of the effect it has on rock mineralogy and physical properties. In this region, we see evidence for a lower crust that overall is intermediate‐mafic in composition (53.7
- Award ID(s):
- 1650365
- NSF-PAR ID:
- 10449818
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 125
- Issue:
- 3
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract 40%) and global ( 35%) continental lithosphere and the underlying inaccessible mantle ( 25%). Geophysical models are combined with geochemical data sets to predict the geoneutrino signal at current and future geoneutrino detectors. We propagated uncertainties, both chemical and physical, through Monte Carlo methods. Estimated total signal uncertainties are on the order of 20%, proportionally with geophysical and geochemical inputs contributing 30% and 70%, respectively. We find that estimated signals, calculated using CRUST2.0, CRUST1.0, and LITHO1.0, are within physical uncertainty of each other, suggesting that the choice of underlying geophysical model will not change results significantly, but will shift the central value by up to 15%. Similarly, we see no significant difference between calculated layer abundances and bulk crustal heat production when using these geophysical models. The bulk crustal heat production is calculated as 7 2 TW, which includes an increase of 1 TW in uncertainty relative to previous studies. Combination of our predicted lithospheric signal with measured signals yield an estimated BSE heat production of 21.5 10.4 TW. Future improvements, including uncertainty attribution and near‐field modeling, are discussed. -
more » « less
Abstract Using data from 186 stations belonging to the USArray Transportable Array, a three‐dimensional shear wave velocity model for the southeastern United States is constructed for the top 180 km by a joint inversion of receiver functions and Rayleigh wave phase velocity dispersion computed from ambient noise and teleseismic earthquake data. The resulting shear wave velocity model and the crustal thickness and Vp/Vs (
) measurements show a clear spatial correspondence with major surficial geological features. The distinct low velocities observed in the depth range of 0–25 km beneath the eastern Gulf Coastal Plain reflect the thick layer of unconsolidated or poorly consolidated sediments atop the crystalline crust. The low (1.70–1.74) and slow lowermost crustal velocities observed beneath the eastern Southern Appalachian Mountains (including the Carolina Terrane and Inner Piedmont) relative to the adjacent Blue Ridge Mountains and Valley and Ridge can be interpreted by lower crustal delamination followed by relamination. The Osceola intrusive complex in the central Suwannee Terrane has similar crustal characteristics as the eastern Southern Appalachian Mountains and thus can similarly be attributed to crustal delamination/relamination processes. The Grenville Province and adjacent areas possess relatively high values which can be attributed to mafic intrusion associated with crustal extension in a recently recognized segments of the eastern arm of the Proterozoic Midcontinent Rift. -
more » « less
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 Indo‐Burma subduction zone is a highly oblique subduction system where the Indian plate is converging with the Eurasian plate. How strain is partitioned between the Indo‐Burma interface and upper plate Kabaw Fault, and whether the megathrust is a locked and active zone of convergence that can generate great earthquakes are ongoing debates. Here, we use data from a total of 68 Global Navigation Satellite System (GNSS) stations, including newly installed stations across the Kabaw Fault and compute an updated horizontal and vertical GNSS velocity field. We correct vertical rates for fluctuating seasonal signals by accounting for the elastic response of monsoon water on the crust. We model the geodetic data by inverting for 11,000 planar and non‐planar megathrust fault geometries and two geologically viable structural interpretations of the Kabaw Fault that we construct from field geological data, considering a basin‐scale wedge‐fault and a crustal‐scale reverse fault. We demonstrate that the Indo‐Burma megathrust is locked, converging at a rate of
mm/yr, and capable of hosting >8.2 M w megathrust events. We also show that the Kabaw Fault is locked and accommodating strike‐slip motion at a rate ofmm/yr and converging at a rate of mm/yr. Our interpretation of the geological, geophysical, and geodetic datasets indicates the Kabaw Fault is a crustal‐scale structure that actively absorbs a portion of the convergence previously ascribed to the Indo‐Burma megathrust. This reveals a previously unrecognized seismic hazard associated with the Kabaw Fault and slightly reduces the estimated hazard posed by megathrust earthquakes in the region. -
more » « less
Abstract A new azimuthal anisotropy model for the North American and Caribbean Plates, namely,
, is constructed based on full waveform inversion and records from the USArray and other temporary/permanent networks deployed in the study region. A total of 180 earthquakes and 4,516 seismographic stations are employed in the inversion to simultaneously constrain radially and azimuthally anisotropic model parameters: , , , and , within the crust and mantle. Thirty‐two preconditioned conjugate gradient iterations have been utilized to minimize frequency‐dependent phase discrepancies between observed and predicted seismograms for three‐component short‐period (15–40 s) body waves and long‐period (25–100 s) surface waves. Model exhibits complicated variations in anisotropic fabrics underneath the western and eastern United States, especially at depths shallower than 100 km. For instance, the fast axis orientations in model suggest the presence of trench‐perpendicular mantle flows underneath the Cascadia Subduction Zone and also follow the strikes of the Snake River Plain, the Ouachita Orogenic Front, and the Grenville and Appalachian Orogenic Belts. The amplitudes of azimuthal anisotropy reduce to around 1% at depths greater than 200 km, and the orientations are subparallel to the global plate motion directions to the east of the Rocky Mountain, except for large discrepancies in central and eastern Canada. At a depth of 700 km, the fast axes change along the trajectory of the Farallon slab underneath the Great Lakes region and Gulf of Mexico, which might indicate the development of 2‐D poloidal‐mode mantle flows perpendicular to the strike of the sinking slab within the uppermost lower mantle. Comparisons between model with a western U.S. model from ambient noise tomography and SKS splitting measurements demonstrate a relatively good agreement for the fast axis orientations, considering the usage of different data sets and imaging techniques. However, the absolute magnitude of azimuthal anisotropy in model might be underestimated, especially at greater depths, given the poor agreement on the amplitudes of predicted and observed SKS splitting times. At the current stage, the agreement among different azimuthal anisotropy models at global and continental scales is still poor even for the United States with a dense station coverage.
