Search for: All records

The dependence of seismic wave speeds on propagation or polarization direction, called seismic anisotropy, is a relatively direct indicator of mantle deformation and flow. Mantle seismic anisotropy is often inferred from measurements of shear-wave splitting. A number of standard techniques to measure shear-wave splitting have been applied globally; for example, *KS splitting is often used to measure upper mantle anisotropy. In order to obtain robust constraints on anisotropic geometry, it is necessary to sample seismic anisotropy from different directions, ideally using different seismic phases with different incidence angles. However, many standard analysis techniques can only be applied for certain epicentral distances and source–receiver geometries. To search for new ways to detect mantle anisotropy, instead of focusing on the sensitivity of individual phases, we investigate the wavefield as a whole: we apply a ‘wavefield differencing’ approach to (systematically) understand what parts of the seismic wavefield are most affected by splitting due to seismic anisotropy in the mantle. We analyze differences between synthetic global wavefields calculated for isotropic and anisotropic input models, incorporating seismic anisotropy at different depths. Our results confirm that the seismic phases that are commonly used in splitting techniques are indeed strongly influenced by mantle anisotropy. However, we also identify less commonly used phases whose waveforms reflect the effects of anisotropy. For example, PS is strongly affected by splitting due to seismic anisotropy in the upper mantle. We show that PS can be used to fill in gaps in global coverage in shear-wave splitting data sets (for example, beneath ocean basins). We find that PcS is also a promising phase, and present a proof-of-concept example of PcS splitting analysis across the contiguous United States using an array processing approach. Because PcS is recorded at much shorter distances than *KS phases, PcS splitting can therefore fill in gaps in backazimuthal coverage. Our wavefield differencing results further hint at additional potential novel methods to detect and characterize splitting due to mantle seismic anisotropy.

The Yellowstone region (western United States) is a commonly cited example of intraplate volcanism whose origin has been a topic of debate for several decades. Recent work has suggested that a deep mantle plume, rooted beneath southern California, is the source of Yellowstone volcanism. Seismic anisotropy, which typically results from deformation, can be used to identify and characterize mantle flow. Here, we show that the proposed plume root location at the base of the mantle is strongly seismically anisotropic. This finding is complemented by geodynamic modeling results showing upwelling flow and high strains in the lowermost mantle beneath the Yellowstone region. Our results support the idea that the Yellowstone volcanism is caused by a plume rooted in the deepest mantle beneath southern California, connecting dynamics in the deepest mantle with phenomena at Earth's surface.

Many regions of the Earth's mantle are seismically anisotropic, including portions of the lowermost mantle, which may indicate deformation due to convective flow. The splitting of ScS phases, which reflect once off the core-mantle boundary (CMB), is commonly measured to identify lowermost mantle anisotropy, although some challenges exist. Here, we use global wavefield simulations to evaluate commonly used approaches to inferring a lowermost mantle contribution to ScS splitting. We show that due to effects of the CMB reflection, only the epicentral distance range between 60° and 70° is appropriate for ScS splitting measurements. For this distance range, splitting is diagnostic of deep mantle anisotropy if no upper mantle anisotropy is present; however, if ScS is also split due to upper mantle anisotropy, the reliable diagnosis of deep mantle anisotropy is challenging. Moreover, even in the case of a homogeneously anisotropic deep mantle region sampled from a single azimuth by multiple ScS waves with different source polarizations (in absence of upper mantle anisotropy), different apparent fast directions are produced. We suggest that ScS splitting should only be measured at null stations and conduct such an analysis worldwide. Our results indicate that seismic anisotropy is globally widespread in the deep mantle.

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

Observations of seismic waves that have passed through the Earth's lowermost mantle provide insight into deep mantle structure and dynamics, often on relatively small spatial scales. Here we use SKS, S2KS, S3KS, and PKS signals recorded across a large region including the United States, Mexico, and Central America to study the deepest mantle beneath large swaths of North America and the northeastern Pacific Ocean. These phases are enhanced via beamforming and then used to investigate polarization‐ and propagation direction‐dependent shear wave speeds (seismic anisotropy). A differential splitting approach enables us to robustly identify contributions from anisotropy. Our results show strong seismic anisotropy in approximately half of our study region, indicating that anisotropy may be more prevalent than commonly thought. In some regions, the anisotropy may be induced by flow driven by sinking cold slabs, and in other, more compact regions, by upwelling flow. Measured splitting due to lowermost mantle anisotropy is sufficiently strong to be non‐negligible in interpretations of SKS splitting due to upper mantle anisotropy in certain regions, which may prompt future re‐evaluations of upper mantle anisotropy beneath North and Central America.

Convective flow in the deep mantle controls Earth's dynamic evolution, influences plate tectonics, and has shaped Earth's current surface features. Present and past convection‐induced deformation manifests itself in seismic anisotropy, which is particularly strong in the mantle's uppermost and lowermost portions. While the general patterns of seismic anisotropy have been mapped for the upper mantle, anisotropy in the lowermost mantle (called D′′) is at an earlier stage of exploration. Here we review recent progress in methods to measure and interpret D′′ anisotropy. Our understanding of the limitations of existing methods and the development of new measurement strategies have been aided enormously by the availability of high‐performance computing resources. We give an overview of how measurements of seismic anisotropy can help constrain the mineralogy and fabric of the deep mantle. Specifically, new and creative strategies that combine multiple types of observations provide much tighter constraints on the geometry of anisotropy than have previously been possible. We also discuss how deep mantle seismic anisotropy provides insights into lowermost mantle dynamics. We summarize what we have learned so far from measurements of D′′ anisotropy, how inferences of lowermost mantle flow from measurements of seismic anisotropy relate to geodynamic models of mantle flow, and what challenges we face going forward. Finally, we discuss some of the important unsolved problems related to the dynamics of the lowermost mantle that can be elucidated in the future by combining observations of seismic anisotropy with geodynamic predictions of lowermost mantle flow.

Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (Sdiff) are potentially strongly influenced by lowermost mantle anisotropy. Sdiff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. Sdiff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SHdiff–SVdiff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of Sdiff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which Sdiff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on Sdiff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SHdiff waves. In particular, we investigate how SHdiff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to Sdiff splitting that can be clearly distinguished from other waveform effects. These results enable us to lay out a strategy for the analysis of Sdiff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of Sdiff and the influence of source-side anisotropy. We demonstrate our Sdiff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained Sdiff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.

Ultralow velocity zones (ULVZs) and seismic anisotropy are both commonly detected in the lowermost mantle at the edges of the two antipodal large low velocity provinces (LLVPs). The preferential occurrences of both ULVZs and anisotropy at LLVP edges are potentially connected to deep mantle dynamics; however, the two phenomena are typically investigated separately. Here we use waveforms from three deep earthquakes to jointly investigate ULVZ structure and lowermost mantle anisotropy near an edge of the Pacific LLVP to the southeast of Hawaii. We model global wave propagation through candidate lowermost mantle structures using AxiSEM3D. Two structures that cause ULVZ‐characteristic postcursors in our data are identified and are modeled as cylindrical ULVZs with radii of ∼1° and ∼3° and velocity reductions of ∼36% and ∼20%. One of these features has not been detected before. The ULVZs are located to the south of Hawaii and are part of the previously detected complex low velocity structure at the base of the mantle in our study region. The waveforms also reveal that, to first order, the base of the mantle in our study region is a broad and thin region of modestly low velocities. Measurements of S_diffshear wave splitting reveal evidence for lowermost mantle anisotropy that is approximately co‐located with ULVZ material. Our measurements of co‐located anisotropy and ULVZ material suggest plausible geodynamic scenarios for flow in the deep mantle near the Pacific LLVP edge.

Determinations of seismic anisotropy, or the dependence of seismic wave velocities on the polarization or propagation direction of the wave, can allow for inferences on the style of deformation and the patterns of flow in the Earth’s interior. While it is relatively straightforward to resolve seismic anisotropy in the uppermost mantle directly beneath a seismic station, measurements of deep mantle anisotropy are more challenging. This is due in large part to the fact that measurements of anisotropy in the deep mantle are typically blurred by the potential influence of upper mantle and/or crustal anisotropy beneath a seismic station. Several shear wave splitting techniques are commonly used that attempt resolve seismic anisotropy in deep mantle by considering the presence of multiple anisotropic layers along a raypath. Examples include source-side S-wave splitting, which is used to characterize anisotropy in the deep upper mantle and mantle transition zone beneath subduction zones, and differential S-ScS and differential SKS-SKKS splitting, which are used to study anisotropy in the D″ layer at the base of the mantle. Each of these methods has a series of assumptions built into them that allow for the consideration of multiple regions of anisotropy. In this work, we systematically assess the accuracy of these assumptions. To do this, we conduct global wavefield modelling using the spectral element solver AxiSEM3D. We compute synthetic seismograms for earth models that include seismic anisotropy at the periods relevant for shear wave splitting measurements (down to 5 s). We apply shear wave splitting algorithms to our synthetic seismograms and analyse whether the assumptions that underpin common measurement techniques are adequate, and whether these techniques can correctly resolve the anisotropy incorporated in our models. Our simulations reveal some inaccuracies and limitations of reliability in various methods. Specifically, explicit corrections for upper mantle anisotropy, which are often used in source-side direct S splitting and S-ScS differential splitting, are typically reliable for the fast polarization direction ϕ but not always for the time lag δt, and their accuracy depends on the details of the upper mantle elastic tensor. We find that several of the assumptions that underpin the S-ScS differential splitting technique are inaccurate under certain conditions, and we suggest modifications to traditional S-ScS differential splitting approaches that lead to improved reliability. We investigate the reliability of differential SKS-SKKS splitting intensity measurements as an indicator for lowermost mantle anisotropy and find that the assumptions built into the splitting intensity formula can break down for strong splitting cases. We suggest some guidelines to ensure the accuracy of SKS-SKKS splitting intensity comparisons that are often used to infer lowermost mantle anisotropy. Finally, we suggest a new strategy to detect lowermost mantle anisotropy which does not rely on explicit upper mantle corrections and use this method to analyse the lowermost mantle beneath east Asia.