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
Lateral Variations in Teleseismic Attenuation of the Conterminous U.S. and New Insights Derived From Its Relationship to Mantle Seismic Velocity
Much of our knowledge of the North American lithosphere comes from imaging seismic velocities. Additional constraints on the subsurface can be gained by studying seismic attenuation, which has different sensitivity to physical properties. We produce a model of lateral variations in attenuation across the conterminous U.S. by analyzing data recorded by the EarthScope Transportable Array. We divide the study area into 12 overlapping tiles and differential attenuation is measured in each tile independently; and twice for four of the tiles. Measurements are combined into a smooth map using a set of linear inversions. Comparing results for adjacent tiles and for repeated tiles shows that the imaged features are robust. The final map shows generally higher attenuation west of the Rocky Mountain Front than east of it, with significant small length scale variations superimposed on that broad pattern. In general, there is a strong anticorrelation between differential attenuation and shear wave velocities at depths of 80–250 km. However, a given change in velocity may correspond to a large or small change in attenuation, depending on the area; suggesting that different physical mechanisms are operating. In the western and south‐central U.S., as well as the Appalachians, velocity variations are large compared to attenuation changes, while the opposite is true in the north‐central and southeastern U.S. Calculations with the Very Broadband Rheology calculator show that these results are consistent with the main source of heterogeneity being temperature and melt fraction in the former regions and grain size variability in the latter ones.
more »
« less
- Award ID(s):
- 1944998
- PAR ID:
- 10504266
- Publisher / Repository:
- American Geophysical Union / Wiley
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 128
- Issue:
- 12
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY A better understanding of damage accumulation before dynamic failure events in geological material is essential to improve seismic hazard assessment. Previous research has demonstrated the sensitivity of seismic velocities to variations in crack geometry, with established evidence indicating that initial crack closure induces rapid changes in velocity. Our study extends these findings by investigating velocity changes by applying coda wave interferometry (CWI). We use an array of 16 piezoceramic transducers to send and record ultrasonic pulses and to determine changes in seismic velocity on intact and faulted Westerly granite samples. Velocity changes are determined from CWI and direct phase arrivals. This study consists of three sets of experiments designed to characterize variations in seismic velocity under various initial and boundary conditions. The first set of experiments tracks velocity changes during hydrostatic compression from 2 and 191 MPa in intact Westerly granite samples. The second set of experiments focuses on saw-cut samples with different roughness and examines the effects of confining pressure increase from 2 to 120 MPa. The dynamic formation of a fracture and the preceding damage accumulation is the focus of the third type of experiment, during which we fractured an initially intact rock sample by increasing the differential stress up to 780 MPa while keeping the sample confined at 75 MPa. The tests show that: (i) The velocity change for rough saw cut samples suggests that the changes in bulk material properties have a more pronounced influence than fault surface apertures or roughness. (ii) Seismic velocities demonstrate higher sensitivity to damage accumulation under increasing differential stress than macroscopic measurements. Axial stress measured by an external load cell deviates from linearity around two-third through the experiment at a stress level of 290 MPa higher than during the initial drop in seismic velocities. (iii) Direct waves exhibit strong anisotropy with increasing differential stress and accumulating damage before rock fracture. Coda waves, on the other hand, effectively average over elastic wave propagation for both fast and slow directions, and the resulting velocity estimates show little evidence for anisotropy. The results demonstrate the sensitivity of seismic velocity to damage evolution at various boundary conditions and progressive microcrack generation with long lead times before dynamic fracture.more » « less
-
Abstract The San Gabriel, Chino, and San Bernardino sedimentary basins in Southern California amplify earthquake ground motions and prolong the duration of shaking due to the basins' shape and low seismic velocities. In the event of a major earthquake rupture along the southern segment of the San Andreas fault, their connection and physical proximity to Los Angeles (LA) can produce a waveguide effect and amplify strong ground motions. Improved estimates of the shape and depth of the sediment‐basement interface are needed for more accurate ground‐shaking models. We obtain a three‐dimensional basement map of the basins by integrating gravity and seismic measurements. The travel time of the sediment‐basementP‐to‐Sconversion, and the Bouguer gravity along 10 seismic lines, are combined to produce a linear relationship that is used to extend the 2D profiles to a 3D basin map. Basement depth is calculated using the predicted travel time constrained by gravity with anS‐wave velocity model of the area. The model is further constrained by the basement depths from 17 boreholes. The basement map shows the south‐central part of the San Gabriel basin is the deepest part and a significant gravity signature is associated with our interpretation of the Raymond fault. The Chino basin deepens toward the south and shallows northeastward. The San Bernardino basin deepens eastward along the edge of the San Jacinto Fault Zone. In addition, we demonstrate the benefit of using gravity data to aid in the interpretation of the sediment‐basement interface in receiver functions.more » « less
-
SUMMARY Temporal changes in seismic velocities are an important tool for tracking structural changes within the crust during transient deformation. Although many geophysical processes span the crust, including volcanic unrest and large-magnitude earthquakes, existing methods for seismic monitoring are limited to the shallow subsurface. We present an approach for deep seismic monitoring based on teleseismic receiver functions, which illuminate the crustal velocity structure from the bottom-up. Using synthetic waveform modelling, we show that receiver functions are uniformly sensitive to velocity changes throughout the crust and can locate the depth of the perturbation. We introduce a novel method based on optimal transport for measuring the non-linear time–amplitude signal variations characteristic of receiver function monitoring. We show that optimal transport enables comparison of full waveform distributions rather than relying on representative stacked waveforms. We further study a linearized version of optimal transport that renders time-warping signal variations into simple Euclidean perturbations, and use this capability to perform blind source separation in the space of waveform variations. This disentangles the effects of changes in the source–receiver path from changes in subsurface velocities. Collectively, these methods extend the reach of seismic monitoring to deep geophysical processes, and provide a tool that can be used to study heterogeneous velocity changes with different spatial extents and temporal dynamics.more » « less
-
SUMMARY The British Isles' lithospheric structure, shaped by a dynamic geological history, remains incompletely understood, particularly regarding anelastic parameters, such as attenuation. In this study, we present a teleseismic attenuation model for the British Isles, using time-domain analysis of teleseismic P-wave data from 28 deep earthquakes. We constructed a 2-D differential attenuation map (Δt*) that reveals significant regional variations. Our findings show a weak anticorrelation between Δt* and shear wave velocity at upper mantle depths, suggesting that variations in the lithosphere–asthenosphere system influence this pattern. A high-attenuation zone extends from Scotland across the Irish Sea to southwest England, potentially linked to mantle upwelling associated with the Iceland plume. This model provides new insights into the mantle dynamics beneath the British Isles, offering a crucial reference for future geophysical studies in the region.more » « less
An official website of the United States government

