An official website of the United States government
Abstract The Earth's long‐ and intermediate‐wavelength geoid anomalies are surface expressions of mantle convection and are sensitive to mantle viscosity. While previous studies of the geoid provide important constraints on the mantle radial viscosity variations, the mantle buoyancy in these studies, as derived from either seismic tomography or slab density models, may suffer significant uncertainties. In this study, we formulate 3‐D spherical mantle convection models with plate motion history since the Cretaceous that generate dynamically self‐consistent mantle thermal and buoyancy structures, and for the first time, use the dynamically generated slab structures and the observed geoid to place important constraints on the mantle viscosity. We found that non‐uniform weak plate margins and strong plate interiors are critical in reproducing the observed geoid and surface plate motion, especially the net lithosphere rotation (i.e., degree‐1 toroidal plate motion). In the best‐fit model, which leads to correlation of 0.61 between the modeled and observed geoid at degrees 4–12, the lower mantle viscosity is ∼1.3–2.5 × 1022 Pa⋅s and is ∼30 and ∼600–1,000 times higher than that in the transition zone and asthenosphere, respectively. Slab structures and the geoid are also strongly affected by slab strength, and the observations prefer moderately strong slabs that are ∼10–100 times stronger than the ambient mantle. Finally, a thin weak layer below the 670‐km phase change on a regional scale only in subduction zones produces stagnant slabs in the mantle transition zone as effectively as a weak layer on a global scale.
more »
« less
Multiscale, radially anisotropic shear wave imaging of the mantle underneath the contiguous United States through joint inversion of USArray and global data sets
SUMMARY EarthScope's USArray seismic component provided unprecedented coverage of the contiguous United States and has therefore spurred significant advances in tomographic imaging and geodynamic modelling. Here, we present a new global, radially anisotropic shear wave velocity tomography model to investigate upper mantle structure and North American Plate dynamics, with a focus on the contiguous United States. The model uses a data-adaptive mesh and traveltimes of both surface waves and body waves to constrain structure in the crust and mantle in order to arrive at a more consistent representation of the subsurface compared to what is provided by existing models. The resulting model is broadly consistent with previous global models at the largest scales, but there are substantial differences under the contiguous United States where we can achieve higher resolution. On these regional scales, the new model contains short wavelength anomalies consistent with regional models derived from USArray data alone. We use the model to explore the geometry of the subducting Farallon Slab, the presence of upper mantle high velocity anomalies, low velocity zones in the central and eastern United States and evaluate models of dynamic topography in the Cordillera. Our models indicate a single, shallowly dipping, discontinuous slab associated with the Farallon Plate, but there are remaining imaging challenges. Inferring dynamic topography from the new model captures both the long-wavelength anomalies common in global models and the short-wavelength anomalies apparent in regional models. Our model thus bridges the gap between high-resolution regional models within the proper uppermost mantle context provided by global models, which is crucial for understanding many of the fundamental questions in continental dynamics.
more »
« less
- PAR ID:
- 10340954
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 226
- Issue:
- 3
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 1730 to 1746
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Topside reverberations off mantle discontinuities are commonly observed at long periods, but their interpretation is complicated because they include both near‐source and near‐receiver reflections. We have developed a method to isolate the stationside reflectors in large data sets with many sources and receivers. Analysis of USArray transverse‐component data from 3,200 earthquakes, using directSas a reference phase, shows clear reflections off the 410‐ and 660‐km discontinuities, which can be used to map the depth and brightness of these features. Because our results are sensitive to the impedance contrast (velocity and density), they provide a useful complement to receiver‐function studies, which are primarily sensitive to theSvelocity jump alone. In addition, reflectors in our images are more spread out in time than in receiver functions, providing good depth resolution. Our images show strong discontinuities near 410 and 660 km across the entire USArray footprint, with intriguing reflectors at shallower depths in many regions. Overall, the discontinuities in the east appear simpler and more monotonous with a uniform transition zone thickness of 250 km compared to the western United States. In the west, we observe more complex discontinuity topography and small‐scale changes below the Great Basin and the Rocky Mountains, and a decrease in transition‐zone thickness along the western coast. We also observe a dipping reflector in the west that aligns with the top of the high‐velocity Farallon slab anomaly seen in some tomography models, but which also may be an artifact caused by near‐surface scattering of incomingSwaves.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. Modelexhibits 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 modelsuggest 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 modelwith 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 modelmight 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.more » « less
-
Abstract EarthScope's USArray Transportable Array has shortcomings for the purpose of interpreting geologic features of wavelengths less than the Transportable Array station spacing, but these can be overcome by using higher spatial resolution gravity data. In this study, we exploit USArray receiver functions to reduce nonuniqueness in the interpretation of gravity anomalies. We model gravity anomalies from previously derived density variations of sedimentary basins, crustalVp/Vsvariation, Moho variation, and upper mantle density variation derived from body wave imaging informed by surface wave tomography to estimateVp/Vs. Although average densities and density contrasts for these seismic variations can be derived, the gravity anomalies modeled from them do not explain the entire observed gravity anomaly field in the United States. We use the unmodeled gravity anomalies (residuals) to reconstruct local variations in densities of the crust associated with geologic sources. The approach uses velocity‐density relationships and differs from density computations that assume isostatic compensation. These intracrustal densities identify geologic sources not sampled by and, in some cases, aliased by the USArray station spacing. We show an example of this improvement in the vicinity of the Bloomfield Pluton, north of the bootheel of Missouri, in the central United States.more » « less
-
Approximately two-thirds of Earth’s outermost shell is composed of oceanic plates that form at spreading ridges and recycle back to Earth’s interior in subduction zones. A series of physical and chemical changes occur in the subducting lithospheric slab as the temperature and pressure increase with depth. In particular, olivine, the most abundant mineral in the upper mantle, progressively transforms to its high-pressure polymorphs near the mantle transition zone, which is bounded by the 410 km and 660 km discontinuities. However, whether olivine still exists in the core of slabs once they penetrate the 660 km discontinuity remains debated. Based on SKS and SKKS shear-wave differential splitting times, we report new evidence that reveals the presence of metastable olivine in the uppermost lower mantle within the ancient Farallon plate beneath the eastern United States. We estimate that the low-density olivine layer in the subducted Farallon slab may compensate the high density of the rest of the slab associated with the low temperature, leading to neutral buoyancy and preventing further sinking of the slab into the deeper part of the lower mantle.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggab185
