Teleseismic Shear Wave Research Articles

Tectonic Implications of Seismic Anisotropy Layering Beneath the Southern Tibetan Plateau Revealed by Integrated Shear Wave Splitting and Receiver Function Analyses

AbstractTo investigate continental dynamics underneath the south‐central Tibetan plateau, which composes the Himalayan, Lhasa, and Qiangtang blocks, we have conducted comprehensive examinations of seismic azimuthal anisotropy in the crust using receiver functions (RFs) and crustal and mantle anisotropy using teleseismic shear wave splitting (SWS) analysis. In the Qiangtang block, the observed predominantly E‐W fast orientations from RF and SWS analyses with similar magnitude are interpreted as resulting from eastward crustal flow with minor contributions from the mantle. In the Lhasa block, the crustal anisotropy is approximately N‐S oriented, which is parallel to the strike of rift basins and southward crustal flow. Anisotropy revealed by SWS demonstrates a rotation from E‐W in the north to NE‐SW in the south, which can be interpreted as reflecting mantle flow field induced by the northward movement of the subducting Indian plate. The addition of PKS and SKKS measurements and extension of epicentral distance range to 171.8° for SWS analysis revealed dominantly strong E‐W oriented anisotropy in most parts of the Himalayan block, where most previous studies reported pervasively null measurements. The absence of azimuthal anisotropy is observed in two regions in the Himalayan block which is attributable to mantle upwelling through a previously identified slab window. A two‐layered anisotropy structure with different fast orientations for the upper and lower layers can be constrained in the southern Qiangtang and the vicinity of the Main Boundary Thrust.

Insight into the Evolution of the Eastern Margin of the Wyoming Craton from Complex, Laterally Variable Shear Wave Splitting

Abstract Dense seismic arrays such as EarthScope’s Transportable Array (TA) have enabled high-resolution seismic observations that show the structure of cratonic lithosphere is more heterogeneous and complex than previously assumed. In this study, we pair TA data with data from the Bighorn Arch Seismic Experiment and the Crust and lithosphere Investigation of the Easternmost expression of the Laramide Orogeny (CIELO) to provide unprecedented detail on the seismic anisotropic structure of the eastern margin of the Wyoming Craton, where several orogens emerged from nominally strong cratonic lithosphere during the Laramide Orogeny. In this study, we use the splitting of teleseismic shear waves to characterize fabrics associated with deformation in the Earth’s crust and mantle. We constrain distinct anisotropic domains in the study area, and forward modeling shows that each of these domains can be explained by a single layer of anisotropy. Most significantly, we find a fast direction in the southern part of the Powder River Basin, which we refer to as the Thunder Basin Block (TBB), that deviates from absolute plate motion (APM). This change in splitting behavior coincides with changes in other modeled geophysical observations, such as active source P-wave velocity models, potential field modeling, and seismic attenuation analysis, which all show a significant change moving from the Bighorn Mountains to the TBB. We argue that these results correspond to structure predating the Laramide Orogeny, and most likely indicate a Neoarchean boundary preserved within the lithosphere.

Discriminating Lithospheric and Asthenospheric Anisotropy Beneath Northern Oman: Sharp Contrast Observed at the Semail Gap Fault Zone

AbstractTo gain a deeper understanding of the extensive and varied lithospheric deformations beneath northern Oman, we examine seismic anisotropy in this region using splitting analysis of teleseismic shear wave data. Our study utilizes data from a dense network consisting of 13 permanent and 45 temporary seismic stations, which were operational for approximately 2.5 years starting from October 2013. By examining the azimuthal distribution of shear wave splitting (SWS) parameters, we were able to divide the study area into three sub‐regions. The stations located to the west of the Hawasina window exhibit relatively azimuthally invariant SWS parameters suggesting a single anisotropic layer. On the other hand, most of the stations located in the central and eastern regions display variations versus back‐azimuth, indicating the potential presence of depth‐dependent anisotropy. The General NW‐SE trend of the Fast Polarization Directions (FPDs) of the one‐layer anisotropy in the west and FPDs of the upper layers in the east is concordant with the strike of the structures resulting from the collision between the continental and oceanic plates. A clear contrast in SWS parameters is observed in the Semail Gap Fault Zone (SGFZ), suggesting that the SGFZ can be a lithospheric‐scale structure that hampers the intrusion of mafic magma from the southeast. Furthermore, the FPDs of the lower layer in the east exhibit an NE‐SW trend, which may be indicative of the large‐scale mantle flow resulting from the present‐day plate motion.

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.

Crustal and Uppermost Mantle Azimuthal Seismic Anisotropy of Antarctica From Ambient Noise Tomography

AbstractSeismic anisotropy provides essential information for characterizing the orientation of deformation and flow in the crust and mantle. The isotropic structure of the Antarctic crust and upper mantle has been determined by previous studies, but the azimuthal anisotropy structure has only been constrained by mantle core phase (SKS) splitting observations. This study determines the azimuthal anisotropic structure of the crust and mantle beneath the central and West Antarctica based on 8—55 s Rayleigh wave phase velocities from ambient noise cross‐correlation. An anisotropic Rayleigh wave phase velocity map was created using a ray—based tomography method. These data are inverted using a Bayesian Monte Carlo method to obtain an azimuthal anisotropy model with uncertainties. The azimuthal anisotropy structure in most of the study region can be fit by a two‐layer structure, with one layer at depths of 0–15 km in the shallow crust and the other layer in the uppermost mantle. The azimuthal anisotropic layer in the shallow crust of West Antarctica, where it coincides with strong positive radial anisotropy quantified by the previous study, shows a fast direction that is subparallel to the inferred extension direction of the West Antarctic Rift System. Fast directions of upper mantle azimuthal anisotropy generally align with teleseismic shear wave splitting fast directions, suggesting a thin lithosphere or similar lithosphere‐asthenosphere deformation. However, inconsistencies in this exist in Marie Byrd Land, indicating differing ancient deformation patterns in the shallow mantle lithosphere sampled by the surface waves and deformation in the deeper mantle and asthenosphere sampled more strongly by splitting measurements.

SKS Polarization Anomalies Due to the Coriolis Force

ABSTRACT The Earth’s Coriolis force has been well-known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear-wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear-wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well-isolated phases. We recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

Imaging upper mantle anisotropy with traveltime and splitting intensity observations from teleseismic shear waves: insights from tomographic reconstructions of subduction simulations

SUMMARYTeleseismic traveltime tomography remains one of the most popular methods for obtaining images of Earth’s upper mantle. However, despite extensive evidence for an elastically anisotropic mantle, the isotropic assumption remains commonplace in such imaging studies. This can result in significant model artefacts which in turn may yield misguided inferences regarding mantle dynamics. The nature of anisotropy-induced apparent velocity anomalies has been well-documented in P-wave imaging and various strategies have been proposed to constrain both isotropic and anisotropic heterogeneity from these data. In contrast, few studies have explored the consequences for shear wave tomography and no practical framework for the anisotropic inversion of S-wave delays exists. Here, we propose a new method for constraining arbitrarily oriented hexagonal anisotropy using both traveltime and splitting intensity observations from direct S phases. Our approach accounts for polarization and finite-frequency effects and allows for isotropic starting models. The imaging method is validated through the tomographic analysis of a realistic synthetic dataset produced from waveform simulations through a geodynamic model of subduction. Results illustrate that neglecting anisotropy produces distortions in slab geometry and the appearance of sub- and supraslab low-velocity zones. Anisotropic inversions remove these artefacts while also constraining geodynamically relevant fabric properties including dip.

Teleseismic shear wave splitting and intracontinental collision deformation of the northern Tibetan Plateau and the eastern Tarim basin

Teleseismic shear wave splitting and intracontinental collision deformation of the northern Tibetan Plateau and the eastern Tarim basin

Mantle Anisotropy in NW Namibia From XKS Splitting: Effects of Asthenospheric Flow, Lithospheric Structures, and Magmatic Underplating

AbstractThe presence of the Etendeka flood basalts in northwestern Namibia is taken as evidence for the activity of the Tristan da Cunha mantle plume during the continental breakup between Africa and South America. We investigate seismic anisotropy beneath NW Namibia by splitting analysis of core‐refracted teleseismic shear waves (XKS phases) to probe mantle flow and lithospheric deformation related to the tectonic history of the region. We present the results of the joint splitting analysis of XKS data collected from 34 onshore stations and 12 ocean‐bottom seismometers. The fast polarization directions (FPDs) are consistent with a model that combines the effects of lithospheric deformation and large‐scale mantle flow due to the NE motion of the African plate. The dominantly NNW‐SSE‐oriented FPDs in the northern part are likely caused by shallow lithospheric structures. Our observations do not show any strong evidence of a pervasive effect of the Tristan da Cunha mantle plume.

Lithospheric S Wave Velocity Variations Beneath the Mackenzie Mountains and Northern Canadian Cordillera

AbstractThe Mackenzie Mountains (MMs) in the Yukon and Northwest Territories, Canada, are an enigmatic mountain range. They are currently uplifting (Leonard et al., 2008, https//doi.org/10.1029/2007JB005456), yet are about 700 km from the nearest plate boundary. Their arcuate shape is distinct and extends over 100 km eastward from the general trend of the Northern Canadian Cordillera. To better assess the cause and conditions of the current uplift, we processed ambient seismic noise data from a linear array of broadband seismographs crossing the mountains, along with other regional seismic stations, to estimate Rayleigh wave phase velocities between 6 and 40 s periods. From this, we estimated phase velocity dispersion and performed a tomographic inversion to estimate VS. Tomography reveals a low‐velocity structure that extends upward from the base of the ∼50–66 km thick lithosphere to the upper crust, and we hypothesize that inferred low density and low rigidity associated with the VS anomaly localizes the ongoing uplift and thrust‐dominated seismicity of the MMs. Additionally, we find relatively low crustal velocities that extend to the west of the MMs, suggesting that strain transfer from the Gulf of Alaska plate boundary plays a driving role as the crust translates to the northeast and buckles up against the craton consistent with the orogenic float hypothesis of Mazzotti and Hyndman (2002, https//doi.org/10.1130/0091-7613(2002)030〈0495:YCASTA〉2.0.CO;2). Finally, we observe lithospheric azimuthal anisotropy with an NW‐SE fast direction. This is nearly orthogonal to teleseismic shear wave splitting measurements in the central MMs, and suggests that asthenosphere flow and lithospheric strain are not aligned in this region.

Heterogeneous seismic anisotropy beneath Madeira and Canary archipelagos revealed by local and teleseismic shear wave splitting

SUMMARYMid-plate upward mantle flow is a key component of global mantle convection, but its patterns are poorly constrained. Seismic anisotropy is the most direct way to infer mantle flow as well as melt distribution, yet the convection patterns associated with plume-like mantle upwelling are understudied due to limited seismic data coverage. Here, we investigate seismic anisotropy beneath the Madeira and Canary hotspots using a dense set of shear wave splitting observations and combining teleseismic and local events recorded by three-component broad-band and short-period seismic stations. Using a total of 26 stations in the Madeira archipelago and 43 stations around the Canary Islands, we obtain 655 high-quality measurements that reveal heterogeneous flow patterns. Although local event results are sparse around most islands, we can observe a small average of S-wave splitting times of 0.16 ± 0.01 s, which significantly increase with source depth beneath El Hierro (>20 km) and Tenerife (>38 km) up to 0.58 ± 0.01 and 0.47 ± 0.05 s. This suggests an influence of melt pocket orientation in magma reservoirs developed at uppermost-mantle depths. Likewise, anisotropy increases significantly beneath the islands with shield stage volcanism (up to 9.81 ± 1.78 per cent at El Hierro, western Canaries, against values up to 1.76 ± 0.73 per cent at Lanzarote, eastern Canaries). On average, teleseismic SKS-wave splitting delay times are large (2.19 ± 0.05 s), indicating sublithospheric mantle flow as the primary source for anisotropy in the region. In the Canaries, the western islands show significantly smaller average SKS delay times (1.93 ± 0.07 s) than the eastern ones (2.25 ± 0.11 s), which could be explained by destructive interference above the mantle upwelling. Despite complex patterns of fast polarization directions throughout both regions, some azimuthal pattern across close stations can be observed and related to present-day mantle flow and anisotropy frozen in the lithosphere since before 60 Ma. Additionally, we infer that the current presence of a mantle plume beneath the archipelagos leads to the associated complex, small-scale heterogeneous anisotropy observations.

Varying Shear Wave Splitting Parameters Suggest Interaction Between Lithosphere and Asthenosphere in Arxan‐Chaihe Volcanic Field, NE China

Over 2 years of broadband seismic data from 27 temporary stations deployed in the Arxan-Chaihe Volcanic Field (ACVF) provides us a valued opportunity to study mantle deformation and related geodynamic processes in NE China intraplate region. We estimated 153 pairs of high-quality teleseismic shear wave splitting parameters at 22 stations and 98 non-splitting measurements at 19 stations. Most of the stations in ACVF are characterized by NNW–SSE trending fast directions with variable splitting delay times, suggesting the fossilized extensional deformation within the lithosphere is locally thermal erased by the hot mantle upwellings associated with the volcanism. In the region south of the 46.5°N, a sharp swing of fast polarization occurs, and both nearly NW-SE and W-E oriented fast directions are observed, which could be explained by small local mantle convection caused by the lithospheric delamination.

Classification of Teleseismic Shear Wave Splitting Measurements: A Convolutional Neural Network Approach

AbstractShear wave splitting (SWS) analysis is widely used to provide critical constraints on crustal and mantle structure and dynamic models. In order to obtain reliable splitting measurements, an essential step is to visually verify all the measurements to reject problematic measurements, a task that is increasingly time consuming due to the exponential increase in the amount of data. In this study, we utilized a convolutional neural network (CNN) based method to automatically select reliable SWS measurements. The CNN was trained by human‐verified teleseismic SWS measurements and tested using synthetic SWS measurements. Application of the trained CNN to broadband seismic data recorded in south central Alaska reveals that CNN classifies 97.0% of human selected measurements as acceptable, and revealed ∼30% additional measurements. To our knowledge, this is the first study to systematically explore the potential of a machine‐learning based technique to assist with SWS analysis.

Seismic Anisotropy and Mantle Flow Constrained by Shear Wave Splitting in Central Myanmar

AbstractThis study represents the first campaign‐style teleseismic shear wave splitting (SWS) investigation of central Myanmar, an area that is tectonically controlled by the oblique subduction of the Indian Plate underneath the Eurasian Plate. The resulting 678 well‐defined and 247 null SWS measurements obtained from recently deployed 71 broadband seismic stations show that the Indo‐Burma Ranges (IBR) possess mostly N‐S fast orientations that are parallel to the trend of the depth contours of the subducted slab. Relative to the global average of 1.0 s, extremely large splitting times with station‐averaged values ranging from 1.28 to 2.79 s and an area‐averaged value of 2.09 ± 0.55 s are observed in the IBR. In contrast, the Central Basin (CB) and the Shan Plateau (SP) are characterized by slightly larger than normal splitting times. The fast orientations observed in the CB are mostly NE‐SW in the northern part of the study area, N‐S in the central part, and NW‐SE in the southern part. The fast orientations change from nearly N‐S along the N‐S oriented Sagaing Fault, to NW‐SE in the central and eastern portions of the SP. These observations, together with SWS measurements using local S events, crustal anisotropy measurements using P‐to‐S receiver functions, and the estimated depth of the source of anisotropy using the spatial coherency of the splitting parameters, suggest the presence of a trench‐parallel sub‐slab flow system driven by slab rollback, a trench‐perpendicular corner flow, and a trench‐parallel flow possibly entering the mantle wedge through a slab window or gap.

Mantle Flow in the Vicinity of the Eastern Edge of the Pacific‐Yakutat Slab: Constraints From Shear Wave Splitting Analyses

AbstractTo investigate the effects of a slab edge and varying slab geometry on the mantle flow systems beneath south central Alaska, a total of 971 pairs of teleseismic shear wave (SKS, SKKS, and PKS) and 65 pairs of local S wave splitting parameters (fast orientations and splitting times) are measured using data from the USArray and other networks. The Pacific‐Yakutat slab edge separates two regions with different characteristics of the splitting measurements. The area to the west of the slab edge has greater splitting times and mostly trench parallel fast orientations, and the area to the east is dominated by smaller splitting times and spatially varying fast orientations. The spatial distribution of the splitting parameters and results of anisotropy layering and depth analyses can be explained by a model involving three flow systems. The sub‐slab flow initially entraining with the shallow‐dipping Yakutat slab deflects to a trench‐parallel direction due to slab retreat and an increase in slab dip, and flows northeastward toward the slab edge, where it splits into two branches. The first branch enters the mantle wedge as a toroidal flow and flows southwestward along the slab, and the second branch continues approximately eastward. The flowlines of the toroidal and continued flow systems are approximately orthogonal to each other in the vicinity of the slab edge, producing the observed small splitting times and spatially varying fast orientations.

Small Shear Wave Splitting Delays Suggest Weak Anisotropy in Cratonic Mantle Lithosphere

AbstractWe use splitting in core‐refracted teleseismic shear waves (SKS, PKS, and similar) to investigate anisotropic properties of the upper mantle beneath the Superior craton in eastern North America and the Yilgarn craton in Western Australia. At four sites in each craton, we assemble extensive data sets that emphasize directional coverage, and use three different measurement methods to develop mutually consistent constraints on the nature of splitting and on the likely anisotropic properties that cause it. In both cratons, we see evidence of clear directional variation in both delays and fast polarization directions, as well as lateral differences between sites. Relatively small (0.3–0.8 s) amounts of splitting imply weak anisotropy within 150–220 km thick mantle lithosphere. Anisotropy in the asthenosphere likely contributes to splitting in North America where fast directions align with absolute plate motion, but not in Western Australia where fast polarizations and plate motion are nearly orthogonal.

Crustal structure of southeast Australia from teleseismic receiver functions

Abstract. In an effort to improve our understanding of the seismic character of the crust beneath southeast Australia and how it relates to the tectonic evolution of the region, we analyse teleseismic earthquakes recorded by 24 temporary and 8 permanent broadband stations using the receiver function method. Due to the proximity of the temporary stations to Bass Strait, only 13 of these stations yielded usable receiver functions, whereas seven permanent stations produced receiver functions for subsequent analysis. Crustal thickness, bulk seismic velocity properties, and internal crustal structure of the southern Tasmanides – an assemblage of Palaeozoic accretionary orogens that occupy eastern Australia – are constrained by H–κ stacking and receiver function inversion, which point to the following: a ∼ 39.0 km thick crust; an intermediate–high Vp/Vs ratio (∼ 1.70–1.76), relative to ak135; and a broad (> 10 km) crust–mantle transition beneath the Lachlan Fold Belt. These results are interpreted to represent magmatic underplating of mafic materials at the base of the crust. a complex crustal structure beneath VanDieland, a putative Precambrian continental fragment embedded in the southernmost Tasmanides, that features strong variability in the crustal thickness (23–37 km) and Vp/Vs ratio (1.65–193), the latter of which likely represents compositional variability and the presence of melt. The complex origins of VanDieland, which comprises multiple continental ribbons, coupled with recent failed rifting and intraplate volcanism, likely contributes to these observations. stations located in the East Tasmania Terrane and eastern Bass Strait (ETT + EB) collectively indicate a crust of uniform thickness (31–32 km), which clearly distinguishes it from VanDieland to the west. Moho depths are also compared with the continent-wide AusMoho model in southeast Australia and are shown to be largely consistent, except in regions where AusMoho has few constraints (e.g. Flinders Island). A joint interpretation of the new results with ambient noise, teleseismic tomography, and teleseismic shear wave splitting anisotropy helps provide new insight into the way that the crust has been shaped by recent events, including failed rifting during the break-up of Australia and Antarctica and recent intraplate volcanism.

Seismic anisotropy in southern Costa Rica confirms upper mantle flow from the Pacific to the Caribbean

AbstractSurrounded by subducting slabs and continental keels, the upper mantle of the Pacific is largely prevented from mixing with surrounding areas. One possible outlet is beneath the southern part of the Central American isthmus, where regional observations of seismic anisotropy, temporal changes in isotopic composition of volcanic eruptions, and considerations of dynamic topography all suggest upper mantle flow from the Pacific to the Caribbean. We derive new constraints on the nature of seismic anisotropy in the upper mantle of southern Costa Rica from observations of birefringence in teleseismic shear waves. Fast and slow components separate by ∼1 s, with faster waves polarized along the 40°–50° (northeast) direction, near-orthogonally to the Central American convergent margin. Our results are consistent with upper mantle flow from the Pacific to the Caribbean and require an opening in the lithosphere subducting under the region.

Crustal seismic anisotropy of the Northeastern Tibetan Plateau and the adjacent areas from shear-wave splitting measurements

SUMMARYThe Northeastern Tibetan Plateau has thickened crust and is still undergoing strong active crustal shortening and deformation. Crustal anisotropy can provide clues to how the crust is currently deforming and evolving. We use an automatic method to analyse the upper-crustal anisotropy of the NE Tibetan Plateau and the adjacent region using local earthquakes recorded at 39 permanent seismic stations during the period 2009–2018. The majority of the dominant fast directions are consistent with the maximum horizontal stress orientation, suggesting that the upper-crustal anisotropy is mainly controlled by the regional or local stress field. Several fault-parallel measurements are observed for stations on or near to the main faults. These fault-parallel fast directions indicate that the main mechanism of upper-crustal anisotropy is associated with shear fabric caused by deformation. Fast directions neither fault-parallel nor stress-parallel are observed at stations lying several kilometres away from fault zones, likely reflecting the combined influence of stress-aligned microcracks and active faults. A comparison between our upper-crustal anisotropy parameters and those inferred from previous anisotropy studies that used receiver function and teleseismic shear wave splitting measurements suggests that the crust has the same deformation mechanisms as mantle anisotropy in the southern part of the Western Qinling Fault, whereas the upper-crustal anisotropic mechanism is different from those of lower crust and mantle anisotropy in the northern part of the Western Qinling Fault. These observations imply that the Western Qinling Fault may be an important boundary fault.

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.