An official website of the United States government
Abstract Although seismic velocity and electrical conductivity are both sensitive to temperature, thermal lithosphere properties are derived almost exclusively from seismic data because conductivity is often too strongly affected by minor highly conductive phases to be a reliable indicator of temperature. However, in certain circumstances, electrical observations can provide strong constraints on mantle temperatures. In the southeastern United States (SEUS), magnetotelluric (MT) data require high resistivity values (>300 Ωm) to at least 200‐km depth. As dry mantle mineral conduction laws provide an upper bound on temperature for an observed resistivity value, the only interpretation is that lithospheric temperatures (<1330 °C) persist to 200 km. However, seismic tomography shows that velocities in this region are generally slightly slow with respect to references models; this observation has led to a view of relatively thin (<150 km), eroded thermal lithosphere beneath the SEUS. We show that MT and seismic (tomography, attenuation, receiver function) results are consistent with thick (~200 km), coherent thermal lithosphere in this region. Reduced seismic velocities (relative to reference models) can be explained by considering the effect of finite grain size (anelasticity). Calculated velocity as a function of temperature is overall slower when including anelastic effects, even at reasonable grain sizes of 1 mm to 1 cm; this permits mantle temperatures that are colder than would typically be inferred. We argue for a geodynamic scenario in which the present thermal lithosphere in the SEUS formed in association with the Central Atlantic Magmatic Province and has subsequently survived intact for ~200 Ma.
more »
« less
This content will become publicly available on January 30, 2026
Global 3D model of mantle attenuation using seismic normal modes
Seismic tomographic models based only on wave velocities have limited ability to distinguish between a thermal or compositional origin for Earth’s 3D structure1 . Complementing wave velocities with attenuation observations can make that distinction, which is fundamental for understanding mantle convection evolution. However, global 3D attenuation models are only available for the upper mantle at present2–5. Here we present a 3D global model of attenuation for the whole mantle made using whole-Earth oscillations, constraining even spherical harmonics up to degree four. In the upper mantle, we find that high attenuation correlates with low velocity, indicating a thermal origin, in agreement with previous studies6,7. In the lower mantle, we find the opposite and observe the highest attenuation in the ‘ring around the Pacific’, which is seismically fast, and the lowest attenuation in the large low-seismic-velocity provinces (LLSVPs). Comparing our model with wave speeds and attenuation predicted by a laboratory-based viscoelastic model8 suggests that the circum-Pacific is a colder and small-grain-size region9, surrounding the warmer and large-grain-size LLSVPs. Viscosities calculated for the inferred variations in grain size and temperature confirm LLSVPs as long-lived, stable features10
more »
« less
- Award ID(s):
- 2125895
- PAR ID:
- 10655669
- Publisher / Repository:
- Nature
- Date Published:
- Journal Name:
- Nature
- Volume:
- 637
- Issue:
- 8048
- ISSN:
- 0028-0836
- Page Range / eLocation ID:
- 1131 to 1135
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Earth's lowermost mantle is characterized by two large low shear velocity provinces (LLSVPs). The regions outside the LLSVPs have been suggested to be strongly influenced by subducted slabs and, therefore, much colder than the LLSVPs. However, localized low‐velocity seismic anomalies have been detected in the subduction‐influenced regions, whose origin remains unclear. Here, three‐dimensional geodynamic calculations are performed, and they show that linear, ridge‐like hot thermal anomalies, or thermal ridges, form in the relatively cold, downwelling regions of the lowermost mantle. Like the formation of Richter rolls due to sublithosphere small‐scale convection (SSC), the thermal ridges form as a result of SSC from the basal thermal boundary layer and they extend in directions parallel to the surrounding mantle flow. The formation of thermal ridges in subduction regions of the lowermost mantle is very sensitive to the thermal structures of the subducted materials, and thermal heterogeneities brought to the bottom of the mantle by subducting slabs greatly promote the formation of thermal ridges. The formation of thermal ridges is also facilitated by the increase of core‐mantle boundary heat flux and vigor of lowermost mantle convection. The thermal ridges may explain the low‐velocity seismic anomalies outside of the LLSVPs in the lowermost mantle. The results suggest that the relatively cold, subduction‐influenced regions of the Earth's lowermost mantle may contain localized hot anomalies.more » « less
-
Abstract We present two new seismic velocity models for Alaska from joint inversions of body-wave and ambient-noise-derived surface-wave data, using two different methods. Our work takes advantage of data from many recent temporary seismic networks, including the Incorporated Research Institutions for Seismology Alaska Transportable Array, Southern Alaska Lithosphere and Mantle Observation Network, and onshore stations of the Alaska Amphibious Community Seismic Experiment. The first model primarily covers south-central Alaska and uses body-wave arrival times with Rayleigh-wave group-velocity maps accounting for their period-dependent lateral sensitivity. The second model results from direct inversion of body-wave arrival times and surface-wave phase travel times, and covers the entire state of Alaska. The two models provide 3D compressional- (VP) and shear-wave velocity (VS) information at depths ∼0–100 km. There are many similarities as well as differences between the two models. The first model provides a clear image of the high-velocity subducting plate and the low-velocity mantle wedge, in terms of the seismic velocities and the VP/VS ratio. The statewide model provides clearer images of many features such as sedimentary basins, a high-velocity anomaly in the mantle wedge under the Denali volcanic gap, low VP in the lower crust under Brooks Range, and low velocities at the eastern edge of Yakutat terrane under the Wrangell volcanic field. From simultaneously relocated earthquakes, we also find that the depth to the subducting Pacific plate beneath southern Alaska appears to be deeper than previous models.more » « less
-
SUMMARY Seismic tomography is a principal method for studying mantle structure, but imaging of Earth’s wave speed anomalies is conditioned by seismic wave sampling. Global models use misfit criteria that may strive for balance between portions of the data set but can leave important regional domains underserved. We evaluate two full-waveform global tomography wave speed models, GLAD-M25 (Global adjoint tomography model) and SEMUCB-WM1 (whole-mantle tomography model derived from fully numerical spectral element method forward modelling), in the mantle below the Pacific Ocean. The region of the South Pacific Superswell contains multiple hotspots which may be fed by plumes anchored in the Large Low Shear-Velocity Province at the base of the mantle. The uneven distribution of seismic receivers worldwide leaves several candidate plumes beneath various hotspots poorly resolved. We assess the regional quality of GLAD-M25 relative to its global performance using a partition of the seismic waveform data used in its construction. We evaluate synthetic waveforms computed using the spectral-element method to determine how well they fit the data according to a variety of criteria measured across multiple seismic phases and frequency bands. The distributions of traveltime anomalies that remain in GLAD-M25 are wider for trans-Pacific paths than globally, suggesting comparatively insufficiently resolved seismic velocity structure in the region of interest. Hence, Pacific-centred regional inversions, based on (augmented) subsets of the global data set have the potential to enhance the resolution of velocity structure. We compare GLAD-M25 and SEMUCB-WM1 by cross-validation with a new, independent, data set. Our results reveal that short- and long-wavelength structure is captured differently by the two models. Our findings lead us to recommend focusing future model iteration on and around the Pacific Superswell and adding data that sample new corridors, especially using ocean sensors, to better constrain seismic velocity structure in this area of significant geodynamic complexity.more » « less
-
Abstract Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.more » « less
