We investigate seismic anisotropy in the lowermost mantle in the vicinity of the African large low shear velocity province (LLSVP) using observations of differential
Iceland represents one of the most well-known examples of hotspot volcanism, but the details of how surface volcanism connects to geodynamic processes in the deep mantle remain poorly understood. Recent work has identified evidence for an ultra-low velocity zone in the lowermost mantle beneath Iceland and argued for a cylindrically symmetric upwelling at the base of a deep mantle plume. This scenario makes a specific prediction about flow and deformation in the lowermost mantle, which can potentially be tested with observations of seismic anisotropy. Here we present an investigation of seismic anisotropy in the lowermost mantle beneath Iceland, using differential shear wave splitting measurements of S–ScS and SKS–SKKS phases. We apply our techniques to waves propagating at multiple azimuths, with the goal of gaining good geographical and azimuthal coverage of the region. Practical limitations imposed by the suboptimal distribution of global seismicity at the relevant distance ranges resulted in a relatively small data set, particularly for S–ScS. Despite this, however, our measurements of ScS splitting due to lowermost mantle anisotropy clearly show a rotation of the fast splitting direction from nearly horizontal for two sets of paths that sample away from the low velocity region (implying VSH > VSV) to nearly vertical for a set of paths that sample directly beneath Iceland (implying VSV > VSH). We also find evidence for sporadic SKS–SKKS discrepancies beneath our study region; while the geographic distribution of discrepant pairs is scattered, those pairs that sample closest to the base of the Iceland plume tend to be discrepant. Our measurements do not uniquely constrain the pattern of mantle flow. However, we carried out simple ray-theoretical forward modelling for a suite of plausible anisotropy mechanisms, including those based on single-crystal elastic tensors, those obtained via effective medium modelling for partial melt scenarios, and those derived from global or regional models of flow and texture development in the deep mantle. These simplified models do not take into account details such as possible transitions in anisotropy mechanism or deformation regime, and test a simplified flow field (vertical flow beneath the plume and horizontal flow outside it) rather than more detailed flow scenarios. Nevertheless, our modelling results demonstrate that our ScS splitting observations are generally consistent with a flow scenario that invokes nearly vertical flow directly beneath the Iceland hotspot, with horizontal flow just outside this region.
more » « less- NSF-PAR ID:
- 10121500
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 219
- Issue:
- Supplement_1
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- p. S152-S166
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract SKS ‐SKKS shear‐wave splitting. We use data from 375 permanent and temporary stations in Africa which enable us to map the spatial distribution of the anisotropic regions of the lowermost mantle in unprecedented detail. Our results corroborate previous findings that anisotropy is most clearly observed at the margins of the LLSVP, indicating strong deformation at its border, and they are generally consistent with a mostly isotropic LLSVP interior. We find that most discrepantSKS ‐SKKS measurements sample the lowermost mantle close to what is inferred to be the root of the Afar plume. We also identify strongly discrepant splitting in the vicinity of a previously mapped ultralow velocity zone (ULVZ) at the base of the LLSVP, beneath Central Africa. This represents an unusual observation of lowermost mantle anisotropy that is spatially coincident with a ULVZ and may reflect a unique anisotropic mechanism such as alignment of partial melt or the presence of strongly anisotropic magnesiowüstite. We interpret discrepant measurements outside of the LLSVP as likely reflecting a change in flow direction from the horizontal plane to a more vertical direction, which may be caused by deflection at the steep LLSVP border. We propose that our observations of D″ anisotropy associated with the African LLSVP can be explained by a mantle flow regime that maintains passive thermochemical piles with slab‐driven flow and allows for the formation of upwellings at their edges. -
more » « less
SUMMARY 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.
-
more » « less
SUMMARY It is well known that regions of the lowermost mantle—D″—exhibit significant seismic anisotropy. Identifying a unique mechanism for seismic anisotropy in D″ and interpreting results in terms of mantle flow has proved challenging. In an attempt to address this, we outline a method for the direct inversion of shear wave waveform data for the orientation and strength of seismic anisotropy. We demonstrate our method by jointly inverting SKS, SKKS and ScS shear wave data for seismic anisotropy in a fast shear wave velocity anomaly beneath the Eastern Pacific Ocean. Using our inversion method we evaluate four candidate mechanisms for seismic anisotropy in D″: elliptical transverse isotropy (representing layering or inclusions), bridgmanite and post-perovskite (for fabrics dominated by either [100](001) or [100](010) slip). We find that all candidate mechanisms can reasonably explain our input data, with synthetic inversions demonstrating that improved backazimuthal coverage is required to identity a single best-fitting mechanism. By inverting for orientation and anisotropic strength parameters we are able to discount bridgmanite as a candidate mechanism as less plausible solution, as our inversion requires an unreasonable ca. 40 per cent of D″ to consist of aligned bridgmanite crystals. The inversion results for the 4 candidate mechanisms predict two different mantle flow regimes, near vertical upwelling (or downwelling) or predominantly horizontal Southwesterly (or Northwesterly) deformation, both of which are inconsistent with recent mantle flow models. These results show that our new inversion method gives seismologists a powerful new tool to constrain lowermost mantle anisotropy, allowing us to test predictions of lowermost mantle flow.
-
more » « less
Abstract Seismic azimuthal anisotropy beneath Australia is investigated using splitting of the teleseismic PKS, SKKS, and SKS phases to delineate asthenospheric flow and lithospheric deformation beneath one of the oldest and fast‐moving continents on Earth. In total 511 pairs of high‐quality splitting parameters were observed at 116 seismic stations. Unlike other stable continental areas in Africa, East Asia, and North America, where spatially consistent splitting parameters dominate, the fast orientations and splitting times observed in Australia show a complex pattern, with a slightly smaller than normal average splitting time of 0.85 ± 0.33 s. On the North Australian Craton, the fast orientations are mostly N‐S, which is parallel to the absolute plate motion (APM) direction in the hotspot frame. Those observed in the South Australian Craton are mostly NE‐SW and E‐W, which are perpendicular to the maximum lithospheric horizontal shortening direction. In east Australia, the observed azimuthal anisotropy can be attributed to either APM induced simple shear or lithospheric fabric parallel to the strike of the orogenic belts. The observed spatial variations of the seismic azimuthal anisotropy, when combined with results from depth estimation utilizing the spatial coherency of the splitting parameters and seismic tomography studies, suggest that the azimuthal anisotropy in Australia can mostly be related to simple shear in the rheologically transition layer between the lithosphere and asthenosphere. Non‐APM parallel anisotropy is attributable to modulations of the mantle flow system by undulations of the bottom of the lithosphere, with a spatially variable degree of contribution from lithospheric fabric.
-
more » « less
Abstract To systematically investigate seismic azimuthal anisotropy in the Sumatra subduction zone and probe mantle dynamics associated with the subduction of the Australian Plate beneath the Sunda Plate, a total of 169 pairs of teleseismic XKS (including PKS, SKKS, SKS) and 115 pairs of local
S splitting parameters are obtained using broadband seismic data recorded at ~70 stations. Additionally, crustal anisotropy in the overriding Sunda Plate is measured by analyzing the moveout ofP ‐to‐S conversions from the Moho using a sinusoidal function. Comparison between the three sets of anisotropy measurements obtained using shear waves with different depths of origin suggests that (1) the crust of the Sunda Plate is anisotropic with mostly trench‐parallel fast orientations and a mean splitting time of 0.28 ± 0.05 s; (2) the mantle wedge is azimuthally anisotropic with dominantly trench‐parallel fast orientations and splitting times ranging from 0.22 to 0.81 s, which generally increase with the focal depth; and (3) subslab anisotropy is mostly trench‐normal beneath the fore‐arc region with an averaged splitting time of 1.48 ± 0.06 s, and becomes trench‐parallel beneath the arc and back‐arc areas with a mean splitting time of 0.33 ± 0.04 s. The resulting lateral and vertical distributions of anisotropy obtained using splitting of three types of shear waves advocate the presence of an entrained subslab flow that is deflected by the mantle transition zone. The flow enters the mantle wedge through a slab window and flows horizontally parallel to the trench.
