Seismic Anisotropy Patterns Research Articles

Splitting intensity tomography to image depth-dependent seismic anisotropy patterns beneath the Italian Peninsula and surrounding regions

The region between central Europe and the centre of the Mediterranean is characterised by complex tectonics and kinematics. Here, the interaction between thickened crust, subducting lithosphere and surrounding asthenosphere produces strong and pervasive anisotropy in the upper mantle. Shear wave splitting measurements, the most adopted method to image seismic anisotropy so far, when interpreted in a ray-based framework result in little or no depth resolution, hampering a correct image of the anisotropy distribution with depth. In this study, we aim to better constrain the depth-dependent seismic anisotropy beneath Italy and surrounding regions, by isolating for the first time the source region of anisotropy at different depths. To do that, we perform an anisotropy tomography, adopting the splitting intensity inversion method. It is entirely based on the finite-frequency effect in the splitting of SKS waves. We first computed the splitting intensity using SKS waves recorded at all available permanent and temporary stations over the region, obtaining a huge dataset of measurements used as an input for the tomographic inversion. The large-scale 3D model of seismic anisotropy obtained with the inversion shows a clear change of anisotropy properties in terms of fast polarisation direction and intensity for different depths, thus improving the characterization of the main sources of anisotropy in the mantle as a function of depth. Shallower layers (70–100 km depth) are characterised by a complex and variable oriented pattern of anisotropy fast direction and intensity, which becomes progressively more organised with depth (100–300 km). This pattern suggests a strong control exerted by the geometry and motion of the different slab segments and the large-scale asthenospheric flow generated by subduction and roll-back processes. The strength of anisotropy increases with depth, with high values affecting the bulge of the Alps and Apennines chains and the southern Tyrrhenian subduction system. On the contrary, weaker anisotropy characterises the transition zone from the Apennines to Alpine domains beneath the Po plain, and both the Adriatic and European domains.

Anisotropy and XKS splitting from geodynamic models of double subduction: testing the limits of interpretation

SUMMARY The analysis of the splitting signature of XKS phases is crucial for constraining seismic anisotropy patterns, especially in complex subduction settings such as outward-dipping double subduction. A natural example of this is found in the Central Mediterranean, where the Apennine and the Dinaride slabs subduct in opposite directions, with the Adriatic plate separating them. To assess the capability of XKS-splitting analysis in revealing anisotropic seismic properties, such as fast polarization directions and shear wave anisotropy (in per cent), we use three-dimensional numerical geodynamic models combined with texture evolution simulations. In these models, two identical outward-dipping oceanic plates are separated by a continental plate. Using the full elastic tensors – directly derived from the texture evolution simulations – we compute anisotropic seismic properties and synthetic teleseismic waveforms. From these waveforms synthetic observables are determined, including apparent splitting parameters (fast polarization directions and delay times) and splitting intensities. Based on these observables, we (1) derive models for a single anisotropic layer (one-layer model), (2) identify regions with significant depth-dependent anisotropic seismic properties, and (3) perform inversions at selected locations in terms of two anisotropic layers (two-layer model). We consider two geodynamic models: one with a strong (M1) and one with a weak (M2) continental plate. Model M1 exhibits significant retreat of the subducting plates with no horizontal stretching of the continental plate, whereas Model M2 shows less retreat, substantial horizontal stretching, and detachment of the subducting plates. These different subduction styles result in distinct flow and deformation patterns in the upper mantle, which are reflected in the anisotropic seismic properties. In Model M1, the fast polarization directions below the continental plate are predominantly trench-parallel, whereas in Model M2, they are mostly trench-normal. In most regions of both models, the one-layer models are sufficient to resolve the anisotropic seismic properties, as these properties are nearly constant with depth. However, for both models, we identify some isolated regions – primarily near the tips of the subducting plates and beneath the continental plate – where fast polarization directions exhibit significant variations with depth. Inverting the apparent splitting parameters in these regions yields multiple two-layer models at each location that excellently fit the observables. However, their anisotropic seismic properties can vary significantly, and not all these two-layer models adequately approximate the true depth variations. This ambiguity can be partially reduced by selecting two-layer models in which the summed shear wave anisotropy closely matches that of one of the one-layer models, as these models better capture the true variations.

Corrigendum: Dynamic subsurface changes on El Hierro and La Palma during volcanic unrest revealed by temporal variations in seismic anisotropy patterns

Corrigendum: Dynamic subsurface changes on El Hierro and La Palma during volcanic unrest revealed by temporal variations in seismic anisotropy patterns

Advances in Mapping Lowermost Mantle Convective Flow With Seismic Anisotropy Observations

AbstractConvective 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.

Toroidal mantle flow beneath the NE termination of the Kuril–Kamchatka subduction zone from seismic anisotropy

SUMMARY This study presents the findings of a splitting analysis conducted on core-refracted teleseismic shear waves (SKS, SKKS and PKS, called together as XKS) and local shear waves, obtained from a dense seismological network spanning the Kamchatka Peninsula. The objective of the study is to examine the pattern of mantle flow beneath the study area through the investigation of seismic anisotropy. The peninsula is situated at the northeastern end of the Kuril–Kamchatka subduction zone, where the Kuril trench intersects with the western boundary of the Aleutian trench. The data set utilized in this study comprises waveform data from a dense network of seismic stations (99 broad-band and short-period stations for the local shear wave splitting analysis and 69 broad-band stations for the SKS splitting analysis). The seismograms were downloaded from publicly available data repositories including the IRIS Data Management Center and the GFZ Data Services (GEOFON program). The dense station coverage allows us to investigate the lateral variations in anisotropy, providing insights into the flow patterns within the mantle. The processing of the combined data sets of local shear wave and teleseismic XKS waves allowed us to partially decipher the source of anisotropy in the mantle. Small delay (splitting) times (∼0.35 s) observed from the local-S data suggest that anisotropy in the mantle wedge is relatively weak with lateral variations. Larger splitting times (∼1.1 s) observed for the XKS waves relative to local S suggest that the main part of splitting on the XKS waves occurs in the subslab mantle. On the other hand, the rotational pattern of seismic anisotropy observed by both the local S and XKS waves suggests the presence of a toroidal flow at the NE edge of the subducting slab, which affects both the mantle wedge and subslab mantle. For the regions away from the edge of the slab, the mantle flow seems to be governed mainly by the drag of the lithospheric plate over the underlying asthenosphere.

Dynamic subsurface changes on El Hierro and La Palma during volcanic unrest revealed by temporal variations in seismic anisotropy patterns

Active hotspot volcanism is the surface expression of ongoing dynamic subsurface changes, such as the generation, transport, and stalling of magmas within the upper mantle and crust. Magmatic influx and migration affects local stress patterns in the crust and lithospheric mantle, which influences seismic anisotropy. A better understanding of those patterns helps improve robustness of models forecasting the likelihood of an eruption and prolonged seismicity, with detailed studies being required to observe the significant variations that can occur on small spatial and temporal scales. Here, we investigate seismic anisotropy before, during and after volcanic eruptions. We use local seismicity around El Hierro and La Palma, the two westernmost islands in the Canaries and sites of the most recent volcanic eruptions in the archipelago. We obtained 215 results in El Hierro during and after the 2011/2012 eruption with five three-component broadband seismic stations and 908 results around the 2021 eruption in La Palma with two three-component broadband stations. On La Palma, the majority of seismicity and splitting results are recorded during the eruption and simultaneous deflation of the island. Seismicity locations do not change significantly and fast shear wave polarisation direction is mostly constant, but some variation can be attributed to changes in the magmatic plumbing system. On El Hierro, the general radial pattern reflects stresses induced by the overall uplift of the island during multiple magma intrusion events. Temporal subsets reveal significant variations in location and depth of the events, as well as significant variations in fast polarisation direction caused by ongoing dynamic changes of under- and overpressurisation. An increase of results starting in 2018 hints towards renewed subsurface activity within deeper parts of the plumbing system, affecting the rate of overall seismicity but not any vertical movement of the island.

Upper Mantle Anisotropy beneath Northern Algeria from Shear-Wave Splitting: Anisotropy beneath Northern Algeria

Northern Algeria is a large region in the north-western of Africa, lying on the collision boundary between the African and Eurasian plates. Few studies on Lithosphere and Mantle deformation have been conducted in this region. To better understand the seismic anisotropy patterns beneath this area, we used data from 17 broadband stations installed in Tellian and the Saharan atlases using the SKS shear wave splitting method by processing hundreds of teleseismic events. To estimate the seismic anisotropy, two parameters were calculated, the fast polarization direction and the delay time for each station-event pair. The results show that the seismic anisotropy can be described by two main orientations, ENE-WSW follows the general trend of the Saharan Atlas, particularly in the central and western parts, and ESE-WNW follows the Hodna Mountains in south-eastern Algeria. Our results show that the anisotropy can be explained by single and heterogeneity in the anisotropic structure, where the measurements are very scattered, and the delay time and fast direction changed with the events backazimuth. In the ABSD, CBBR and CDCN stations, which lie in the arc between the Sahara Atlas and the Aurès Mountains, the origin hypothesis of the seismic anisotropy can be linked to the existence of a detached slab. Moreover, in ABZH, OKGL and EARB stations located in the Kabylide and Western regions, it can be associated with the Gibraltar slab. Both slabs sinking in the African margin mantle were previously imaged by seismic tomography. For the remaining stations, the single-layer best explains the observed seismic anisotropy from their regular fast polarization direction. The comparison of the obtained fast directions with GPS measurements shows that anisotropy fast axes are nearly perpendicular to the convergence direction between the African and Eurasia plates.

Northeastward expansion of the Tibetan Plateau: Seismic anisotropy evidence from shear-wave splitting measurements

• Shear-wave splitting in and around the NE Tibetan Plateau mainly reflects anisotropy in the lithospheric mantle. • Dominant NW-SE anisotropy in the lithospheric mantle indicates NE-SW compression. • Fossil anisotropies exist in the interior Alxa Block and the western Ordos Block. The northeastern Tibetan Plateau (NETP) is a key region to clarify the nature and mechanism of the Tibetan Plateau expansion. This region includes the eastern Qilian Orogenic Belt, the Alxa Block, the Helan-Yinchuan Tectonic Belt and the Ordos Block. We make shear-wave splitting measurements using teleseismic waveforms recorded at 17 portable broadband stations that we newly deployed in the study region. Our measurements provide new information on seismic anisotropy and deformation patterns in the lithospheric mantle, which can be used to constrain the extent of the NETP northeastward expansion. The dominant fast polarization direction (FPD) is NW-SE in the eastern Qilian Orogenic Belt, the southern Alxa Block and the Helan-Yinchuan Tectonic Belt, subparallel to the prominent orientation of tectonics in the NETP. The same NW-SE FPD in the three tectonic units indicates NE-SW compression in the lithospheric mantle there caused by the NETP northeastward expansion. NE-SW FPDs are revealed in the interior Alxa Block, which may reflect fossil anisotropy probably formed by subduction and closure of the Paleo-Asian ocean. ENE-WSW FPDs exist in the western Ordos Block, reflecting fossil anisotropy probably caused by convergence of the South China and North China Cratons. Combined with the NW-SE FPDs in the southwestern margin of the Ordos Block revealed by previous studies, our results suggest that the NETP northeastward expansion has obviously affected the southern Alxa Block and the southwestern margin of the Ordos Block. However, its recent northeastward expansion mainly occurred along the Helan-Yinchuan Tectonic Belt.

Alpine Fault‐Related Microstructures and Anisotropy of the Mantle Beneath the Southern Alps, New Zealand

AbstractMantle xenoliths from the Southern Alps, New Zealand, provide insight into the origin of mantle seismic anisotropy related to the Australian‐Pacific plate boundary. Most xenoliths from within 100 km lateral distance of the Alpine Fault are coarse grained, but a small number are finer grained protomylonites. The protomylonites contain connected networks of fine grains with a different crystallographic preferred orientation (CPO) to coarse porphyroclasts in the same xenolith, suggesting that protomylonites and coarse‐grained samples record different deformation kinematics. The CPOs of fine grains in protomylonites have monoclinic symmetry, with the 2‐fold rotation axis normal to a plane that contains olivine [010] and orthopyroxene [100] maxima, suggesting that the protomylonite deformation involved significant simple shear. Some coarse‐grained samples contain unconnected lenses and layers of fine grains with the same CPO as the coarse grains. Microstructures suggest that these fine grains formed by subgrain rotation recrystallization and that protomylonites may represent an up‐strain progression of this microstructure, where the connectivity of fine grains has allowed them to localize shear and develop a new Alpine Fault CPO. The samples tell us about the state of the mantle at 25 Ma, in the early history of the plate boundary. If this suite of samples is representative of the mantle beneath the Alpine Fault in the present day, then we can interpret the complex seismic anisotropy patterns in the lithospheric mantle as representative of blocks containing variably rotated older CPOs juxtaposed by narrow shear zones associated with Alpine Fault deformation.

New SKS splitting anisotropy measurements point to tearing beneath Central Italy

In the middle of the Mediterranean, the partly still active Apennines subduction system has been usually defined using tomographic images and available shear wave splitting measurements. In this paper we describe the new seismic anisotropy dataset for Central Italy, the region where the transition between Northern and Southern Apennines occurs. The new measurements show NW-SE fast polarisation directions beneath the belt, due to the retreat of the slab, NNE-SSW orientations from proper Adriatic mantle sources, and E-W directed anisotropy, attributed to mantle convection flow at the Thyrranian side. Additionally, the new data suggest the presence of a toroidal mantle flow through a tear in the Apenninic slab, from the Adria to the Tyrrhenian side. However, mantle circulation and flows, identified by the pattern of shear wave splitting results, seem different from what was proposed in previous geodynamic models. Indeed, our results support the presence of a vertical slab tear with limited dimension. In the geodynamic model we propose, the tear acts to accommodate a differential slab retreat. The slab partitioning results in a different pattern and strength of seismic anisotropy traced from the Central Apennines with respect to the adjacent Northern and Southern Apennines.

Seismic anisotropy beneath Central Australia: A record of ancient lithospheric deformation

From 2008 to 2011 a broadband seismic array (BILBY) spanned the Australian continent from North to South, crossing the suture zone between the North and South Australian cratons, and traversing the Central Australian orogenic belt. Past tectonic events, including orogenies during the Proterozoic and Paleozoic, have left a long-lasting impression on the crustal and Moho structure of this inter-cratonic region. However, the impact of past tectonic activity on the lithosphere has been less clear. Here we present the first shear-wave splitting results for the BILBY array using a combination of SKS and PKS teleseismic phases to investigate patterns of deformation and seismic anisotropy within the upper mantle beneath central Australia. Null *KS observations are found to be abundant compared to observations of splitting, as has been widely reported by previous studies, but this appears to be largely due to a coincidental alignment of the inferred anisotropic fast direction with the back-azimuthal range at which most available events occur (140°-160°). Across the central Australian belt the station averaged fast directions tend to orientate ENE-WSW parallel with topographic, gravity, and magnetic trends. Northwards however, the fast directions switch orientation, instead following the NW-SE elongated geometry of the Tennant Creek Inlier, thus delineating a sharp lateral change in the underlying seismic anisotropy. Overall, evidence suggests that the splitting pattern likely reflects anisotropy inherent within the lithosphere generated by past deformational events over 300 million plus years ago, as opposed to the present-day mantle flow in the asthenosphere. While two distinct layers of anisotropy, present in both the asthenosphere and lithosphere, is supported by other evidence, it is not necessarily required by our current dataset. Instead, we can sufficiently model our results with only a single layer of anisotropy, consistent with the expected geometry of azimuthal anisotropy in the lithosphere.

Mantle dynamics in the SE Tibetan Plateau revealed by teleseismic shear-wave splitting analysis

The tectonic evolution of the southeastern Tibetan Plateau is key to understanding the mechanism of growth of the whole Tibetan Plateau. Its imprint on lithospheric and asthenospheric mantle rock fabrics is expected to have a strong signature on seismic anisotropy. In this study, we measure the splitting delays and fast polarization directions (FPD) with the minimum-transverse-energy and the splitting intensity methods using core-refracted phases (e.g., SKS, SKKS, PKS) recorded by the ChinArray experiment. The potential complex anisotropic structures, i.e., the frequency-dependent splitting and two-layer anisotropy, are investigated in particular. The FPDs are consistent in different frequency bands, but the delay times measured at low frequency (8–20 s) are often larger than those measured at high frequency (2–8 s). The measurements can be better explained by a two-layer anisotropy model at many stations. Splitting intensities are also inverted to constrain the depth distribution of seismic anisotropy in the upper mantle. Combined with the two-layer anisotropy models, these results suggest different patterns of seismic anisotropy in the lithosphere and asthenosphere. In the lithosphere, anisotropy shows distinct patterns in the different tectonic blocks, while the FPDs are oriented along the major tectonic boundaries. In the asthenosphere, continuous NW-SE FPDs spread from the Tibetan plateau to the Youjiang Orogen along the southwestern margin of the Yangtze craton, implying asthenospheric extrusion from the high plateau. In the northern Indochina block, the dominant E-W FPDs may indicate asthenospheric flow driven by the eastern subduction and subsequent rollback of the Indian plate.

Mantle-flow diversion beneath the Iranian plateau induced by Zagros\u2019 lithospheric keel

Previous investigation of seismic anisotropy indicates the presence of a simple mantle flow regime beneath the Turkish-Anatolian Plateau and Arabian Plate. Numerical modeling suggests that this simple flow is a component of a large-scale global mantle flow associated with the African superplume, which plays a key role in the geodynamic framework of the Arabia-Eurasia continental collision zone. However, the extent and impact of the flow pattern farther east beneath the Iranian Plateau and Zagros remains unclear. While the relatively smoothly varying lithospheric thickness beneath the Anatolian Plateau and Arabian Plate allows progress of the simple mantle flow, the variable lithospheric thickness across the Iranian Plateau is expected to impose additional boundary conditions on the mantle flow field. In this study, for the first time, we use an unprecedented data set of seismic waveforms from a network of 245 seismic stations to examine the mantle flow pattern and lithospheric deformation over the entire region of the Iranian Plateau and Zagros by investigation of seismic anisotropy. We also examine the correlation between the pattern of seismic anisotropy, plate motion using GPS velocities and surface strain fields. Our study reveals a complex pattern of seismic anisotropy that implies a similarly complex mantle flow field. The pattern of seismic anisotropy suggests that the regional simple mantle flow beneath the Arabian Platform and eastern Turkey deflects as a circular flow around the thick Zagros lithosphere. This circular flow merges into a toroidal component beneath the NW Zagros that is likely an indicator of a lateral discontinuity in the lithosphere. Our examination also suggests that the main lithospheric deformation in the Zagros occurs as an axial shortening across the belt, whereas in the eastern Alborz and Kopeh-Dagh a belt-parallel horizontal lithospheric deformation plays a major role.

Seismic Properties of a Unique Olivine-Rich Eclogite in the Western Gneiss Region, Norway

Investigating the seismic properties of natural eclogite is crucial for identifying the composition, density, and mechanical structure of the Earth’s deep crust and mantle. For this purpose, numerous studies have addressed the seismic properties of various types of eclogite, except for a rare eclogite type that contains abundant olivine and orthopyroxene. In this contribution, we calculated the ambient-condition seismic velocities and seismic anisotropies of this eclogite type using an olivine-rich eclogite from northwestern Flemsøya in the Nordøyane ultrahigh-pressure (UHP) domain of the Western Gneiss Region in Norway. Detailed analyses of the seismic properties data suggest that patterns of seismic anisotropy of the Flem eclogite were largely controlled by the strength of the crystal-preferred orientation (CPO) and characterized by significant destructive effects of the CPO interactions, which together, resulted in very weak bulk rock seismic anisotropies (AVp = 1.0–2.5%, max. AVs = 0.6–2.0%). The magnitudes of the seismic anisotropies of the Flem eclogite were similar to those of dry eclogite but much lower than those of gabbro, peridotite, hydrous-phase-bearing eclogite, and blueschist. Furthermore, we found that amphibole CPOs were the main contributors to the higher seismic anisotropies in some amphibole-rich samples. The average seismic velocities of Flem eclogite were greatly affected by the relative volume proportions of omphacite and amphibole. The Vp (8.00–8.33 km/s) and Vs (4.55–4.72 km/s) were remarkably larger than the hydrous-phase-bearing eclogite, blueschist, and gabbro, but lower than dry eclogite and peridotite. The Vp/Vs ratio was almost constant (avg. ≈ 1.765) among Flem eclogite, slightly larger than olivine-free dry eclogite, but similar to peridotite, indicating that an abundance of olivine is the source of their high Vp/Vs ratios. The Vp/Vs ratios of Flem eclogite were also higher than other (non-)retrograded eclogite and significantly lower than those of gabbro. The seismic features derived from the Flem eclogite can thus be used to distinguish olivine-rich eclogite from other common rock types (especially gabbro) in the deep continental crust or subduction channel when high-resolution seismic wave data are available.

Textural and Compositional Changes in the Lithospheric Mantle Atop the Hawaiian Plume: Consequences for Seismic Properties

AbstractWe characterized the texture, composition, and seismic properties of the lithospheric mantle atop the Hawaiian plume by petrostructural analysis of 48 spinel peridotite xenoliths from four localities in three Hawaiian islands. Coarse‐porphyroclastic peridotites with variable degrees of recrystallization, recorded by growth of strain‐free neoblasts onto the deformed microstructure, predominate. Full evolution of this process produced equigranular microstructures. Some peridotites have coarse‐granular microstructures. Coarse‐granular and coarse‐porphyroclastic peridotites have strong orthorhombic or axial‐[100] olivine crystal‐preferred orientations (CPOs). Recrystallization produced some dispersion and, locally, changed the olivine CPO towards axial‐[010]. Enrichment in pyroxenes relative to model melting trends and pyroxenes with interstitial shapes and CPO uncorrelated with the olivine CPO imply refertilization by reactive melt percolation. The unusual spatial distribution of the recrystallized fraction, Ti enrichment, and Rare Earth Element fractionation in recrystallized, equigranular, and coarse‐granular peridotites support that these microstructures are produced by static recrystallization triggered by melt percolation. However, there is no simple relation between microstructure and chemical or modal composition. This, together with marked variations in mineral chemistry among samples, implies multiple spatially heterogeneous melt‐rock reaction events. We interpret the coarse‐porphyroclastic microstructures and CPO as representative of the original oceanic lithosphere fabric. Annealing changed the microstructure to coarse‐granular, but did not modify significantly the olivine CPO. Recrystallization produced moderate dispersion of the CPO. “Normal” oceanic lithosphere seismic anisotropy patterns are therefore preserved. Yet Fe enrichment, refertilization, and limited heating of the base of the lithosphere may reduce seismic velocities by up to 2%, partially explaining negative velocity anomalies imaged at lithospheric depths beneath Hawaii.

Complex upper mantle anisotropy in the Pacific Northwest: Evidence from SKS splitting

We use a dense network of observations and an automated method of analysis to investigate complex patterns of seismic anisotropy in the Pacific Northwest. We present SKS splitting results for approximately 220 broadband seismic stations, including 33 stations from the new Wallowa2 array deployed between 2016 and 2018 in northeast Oregon. Our data set contains approximately 3300 splitting measurements.Our anisotropy measurements indicate that anisotropy in the Pacific Northwest can be divided into two domains: near the subduction zone, anisotropy is perpendicular to the slab; and inboard of the subduction zone anisotropy tends to parallel the velocity structure and is roughly oriented with “absolute” plate motion. However, splitting analysis performs poorly in northeast Oregon, and results are not consistent with uniform or layered anisotropy. We argue that the complex splitting behavior in northern Oregon is a result of spatially variable upper mantle anisotropy, related to the complex upper mantle seismic structure. Interpreting the mantle anisotropy in this region requires methods not limited to standard ray theoretical shear wave splitting analysis.

SKS Splitting Measurements in NE China: New Insights Into the Wudalianchi Intraplate Volcanism and Mantle Dynamics

AbstractTo better understand the Wudalianchi intraplate volcanism and mantle dynamics in NE China, we collect abundant teleseismic waveform data during 2008 to 2018 recorded at 259 seismic stations around the Wudalianchi volcano. These stations include 148 WAVESArray portable stations, 66 NECESSArray portable stations and 45 Chinese provincial stations. Our newly deployed WAVESArray makes the station intervals of ~20–30 km around the Wudalianchi volcano. Using the minimum transverse‐energy grid‐search method, we make effective SKS splitting measurements at 192 stations. Our results show that the fast direction in the study region is mainly NW‐SE, but there are obvious contrasts within and between different tectonic units. To the east of the North‐South Gravity Lineament, around the Wudalianchi volcano, our SKS splitting measurements show a southeastward opening U‐shaped pattern in the fast direction and splitting time. Such a pattern of seismic anisotropy could be caused by horizontal and vertical flows in the big mantle wedge due to the deep subduction, long stagnancy, and deep dehydration of the Pacific slab in the mantle transition zone. The horizontal flow in the big mantle wedge could be obstructed by the thicker lithosphere under the Great Xing'an range and disturbed by hot and wet upwelling flow under the Wudalianchi volcano. Around the Nuominhe and Halaha volcanic groups, variations in the splitting parameters may also be associated with lithospheric deformation and delamination, though they are located at the edge of our study region. These results shed new light on the intraplate volcanism and mantle dynamics in NE China.

Relationships Between Olivine CPO and Deformation Parameters in Naturally Deformed Rocks and Implications for Mantle Seismic Anisotropy

AbstractWe analyze peridotites from a wide range of tectonic settings to investigate relationships between olivine crystallographic preferred orientation (CPO) and deformation conditions in naturally deformed rocks. These samples preserve the five olivine CPO types (A‐ through E‐type) that rock deformation experiments have suggested are controlled by water content, temperature, stress magnitude, and pressure. The naturally deformed specimens newly investigated here (65 samples) and compiled from an extensive literature review (445 samples) reveal that these factors may matter less than deformation history and/or geometry. Some trends support those predicted by experimentally determined parametric dependence, but several observations disagree—namely, that all CPO types are able to form at very low water contents and stresses and that there is no clear relationship between water content and CPO type. This implies that at the low stresses typical of deformation in the mantle, CPO type more commonly varies as a function of strain geometry. Because olivine CPO is primarily responsible for seismic anisotropy in the upper mantle, the results of this study have several implications. These include (1) the many olivine CPO types recorded in samples from individual localities may explain some of the complex seismic anisotropy patterns observed in the continental mantle, and (2) B‐type CPO—where olivine's “fast axes” align perpendicular to flow direction—occurs under many more conditions than traditionally thought. This study highlights the need for more experiments and the difficulty in using olivine CPO in naturally deformed peridotites to infer deformation conditions.

Insights into the structure and dynamics of the upper mantle beneath Bass Strait, southeast Australia, using shear wave splitting

We investigate the structure of the upper mantle using teleseismic shear wave splitting measurements obtained at 32 broadband seismic stations located in Bass Strait and the surrounding region of southeast Australia. Our dataset includes ∼366 individual splitting measurements from SKS and SKKS phases. The pattern of seismic anisotropy from shear wave splitting analysis beneath the study area is complex and does not always correlate with magnetic lineaments or current N-S absolute plate motion. In the eastern Lachlan Fold Belt, fast shear waves are polarized parallel to the structural trend (∼N25E). Further south, fast shear wave polarization directions trend on average N25–75E from the Western Tasmania Terrane through Bass Strait to southern Victoria, which is consistent with the presence of an exotic Precambrian microcontinent in this region as previously postulated. Stations located on and around the Neogene-Quaternary Newer Volcanics Province in southern Victoria display sizeable delay times (∼2.7 s). These values are among the largest in the world and hence require either an unusually large intrinsic anisotropy frozen within the lithosphere, or a contribution from both the lithospheric and asthenospheric mantle. In the Eastern Tasmania Terrane, nearly all observed fast directions are approximately NW-SE. Although part of our data set strongly favours anisotropy originating from “fabric” frozen in the lithospheric mantle, a contribution from the asthenospheric flow related to the present day plate motion is also required to explain the observed splitting parameters. We suggest that deviation of asthenospheric mantle flow around lithospheric roots could be occurring, and so variations in anisotropy related to mantle flow may be expected. Alternatively, the pattern of fast polarisation orientations observed around Bass Strait may be consistent with radial mantle flow associated with a plume linked to the recently discovered Cosgrove volcanic track. However, it is difficult to characterise the relative contributions to the observed splitting from the lithospheric vs. asthenospheric upper mantle due to poor backazimuthal coverage of the data.

Seismic anisotropy and mantle flow beneath East Africa and Arabia

Neoteric geodynamic models have revealed a variety of mantle processes, yet the predicted relationship between the mantle flow and seismic anisotropy remains enigmatic. This ambiguity is caused by the implementation of different tomographic models, viscosity calculations, and the lack of a homogeneous database of shear-wave splitting (SWS) parameters. In this study, I conducted continental scale SWS measurements using data recorded by 285 stations to constrain similarities and difference between the observed the SWS and the calculated ones by recent geodynamical studies in East Africa and Arabia. A comparison of the SWS measurements with previous results of geodynamic modeling suggests that a combination of NNE oriented flow field related to the motion and the subduction of the African/Arabian plates beneath Eurasia over the past 150 Ma are mostly responsible for the seismic anisotropy patterns. These new findings can help eliminate some of the discrepancies among geodynamical studies.