Seismic Anisotropy Research Articles

Effects of background elastic and permeability anisotropy on dynamic seismic signatures of a fluid-saturated porous rock with aligned fractures

Seismic dispersion, attenuation, and frequency-dependent anisotropy due to wave-induced fluid flow in fluid-saturated porous and fractured rocks have been widely studied. However, most of these studies do not consider the effects of background elastic and permeability anisotropy. In this work, we study such effects by considering a fluid-saturated porous rock with penny-shaped fractures embedded in the isotropic plane of a transversely isotropic (TI) background. Starting with P waves propagating perpendicularly to the fracture plane, we obtain the wavefields in the fluid-saturated porous TI background scattered from a single fracture. This allows for the calculation of the effective P wavenumber through the Foldy approximation and the derivation of a frequency-dependent fracture compliance matrix. The angle-dependent seismic dispersion and attenuation, as well as frequency-dependent anisotropy, can thus be calculated. Using this model, we study a fractured tight sandstone with a TI background. The results indicate that the background elastic anisotropy primarily affects the seismic frequency-dependent anisotropy, whereas the background permeability anisotropy not only affects seismic frequency-dependent anisotropy but also influences seismic dispersion and attenuation. In particular, the background permeability anisotropy alters the rate of change in seismic velocities and velocity anisotropy parameters with frequencies and the peaks of seismic attenuation and attenuation anisotropy parameters. By comparing the theoretical predictions of our model with other theories and with two sets of experimental data, our model is validated. Our model can be applied to characterize the fractured layered earth.

Modeling a pure qP wave in attenuative transverse isotropic media using a complex-valued fractional viscoacoustic wave equation

Seismic anisotropy, characterized by direction-dependent seismic velocity and attenuation, is a fundamental property of the earth, widely observed in field and laboratory measurements. This anisotropy significantly influences seismic wavefield attributes, such as amplitude, phase, and traveltime. Accurate quantification of anisotropy-related seismic wave propagation is essential for global- and regional-scale studies in seismology, such as imaging and inversion. Because the elastic assumption of the earth requires contending with S waves, which are often weak in our seismic recording, we derive a new general complex-valued spatial fractional Laplacian (SFL) viscoacoustic anisotropic pure quasi-P (qP) wave equation (WE) in attenuative transverse isotropic (A-TI) media under the acoustic assumption and use the fixed-order SFL operator to quantify seismic wave propagation in heterogeneous media. Our WE accurately captures anisotropic effects in seismic velocity and attenuation, even in strongly attenuative and anisotropic environments. It also describes the decoupled energy attenuation and velocity dispersion behaviors of the frequency-independent quality factor Q, beneficial for seismic imaging and inversion. Compared with existing viscoacoustic anisotropic pure qP WEs in A-TI media, our WE simplifies computational complexity and is easy to implement, as it decomposes into a complex-valued scalar operator and two differential operators. The evaluated dispersion curves further confirm the accuracy of our approach. Numerical examples in complex geologic settings demonstrate the new equation’s precision, stability, and efficiency in modeling seismic wave propagation.

Firn seismic anisotropy in the Northeast Greenland Ice Stream from ambient-noise surface waves

Abstract. We analyse ambient-noise seismic data from 23 three-component seismic nodes to study firn velocity structure and seismic anisotropy near the EastGRIP camp along the Northeast Greenland Ice Stream (NEGIS). Using nine-component correlation tensors, we derive dispersion curves of Rayleigh and Love wave group velocities from 3 to 40 Hz. These velocity distributions exhibit anisotropy along and across the flow. To assess these variations, we invert dispersion curves for shear wave velocities (Vsh and Vsv) in the top 150 m of the NEGIS using a Markov chain Monte Carlo approach. The reconstructed 1-D shear velocity model reveals radial anisotropy in the firn, with Vsh 12 %–15 % greater than Vsv, peaking at the critical density (550 kg m−3). We combine density data from firn cores drilled in 2016 and 2018 to create a new density parameterisation for the NEGIS, serving as a reference for our results. We link seismic anisotropy in the NEGIS to effective and intrinsic causes. Seasonal densification, wind crusts, and melt layers induce effective anisotropy, leading to faster Vsh waves. Changes in firn recrystallisation cause intrinsic anisotropy, altering the Vsv / Vsh ratio. We observe a shallower firn–ice transition across the flow (≈ 50 m) compared with along the flow (≈ 60 m), suggesting increased firn compaction due to the predominant wind direction and increased deformation towards the shear margin. We demonstrate that short-duration (9 d minimum), passive, seismic deployments and noise-based analysis can determine seismic anisotropy in firn, and we reveal 2-D firn structure and variability.

New constraints to subduction zone arc and backarc mantle temperatures: a test of the corner flow model

The heat source for the hot volcanic arc of subduction zones and the uniformly warm backarc is not well understood, particularly where backarc extension is absent. In the widely studied corner flow model, the traction of the underthrusting plate drags down the overlying viscous mantle, and the induced shallower seaward flow brings in heat from the landward and deeper asthenosphere. However, earlier comparisons with observations, especially backarc heat flow, gave poor agreement. There are several new constraints to the temperatures and structures above the underthrusting plate and in the backarc, especially those in western North America, for comparison with model predictions. Seismic receiver functions map horizontal boundaries, including the lithosphere-asthenosphere boundary (LAB). Seismic shear wave tomography gives mantle velocity-depth that provides an approximate proxy for temperature. Attenuation of seismic wave speed indicates high temperature or partial melting. The presence of recent backarc volcanics and their geochemistry give the depth and temperature of the LAB. From these constraints in western North America, the LAB reflects ponding of partial melt, is consistently at a depth of 65±3 km and temperature of 1350±25°C, and overlies nearly constant temperature (adiabatic) asthenosphere across the backarc from Mexico to Alaska. These depth and temperature estimates are greatly different from corner flow model predictions. An alternative hypothesis more consistent with these constraints is vigorous small-scale convection in the backarc asthenosphere. Whether or how this convection coexists with the large-scale flow patterns inferred from the interpretation of seismic anisotropy analyses, and to what degree it prevails in the arc-forearc region, deserve further research.

Modeling the effects of pore aspect ratio, porosity, and seismic anisotropy on wave velocity dispersion and attenuation patterns in oil- and brine-saturated carbonates using a dynamic self-consistent anisotropic approach

Modeling the effects of pore aspect ratio, porosity, and seismic anisotropy on wave velocity dispersion and attenuation patterns in oil- and brine-saturated carbonates using a dynamic self-consistent anisotropic approach

Prediction of shale bedding-parallel crack density based on the SH-SH wave inversion of the VTI media

Bedding-parallel cracks relate to hydrocarbon accumulation and can be beneficial for hydraulic fracturing. The prediction of bedding-parallel crack density (BCD) from seismic data is difficult because the development of bedding-parallel cracks has limited effects on the velocities and density of subsurface formation. Well logs of the Qingshankou (K2qn) shales in the Songliao Basin, northeast China, reveal that the target formations represent extremely strong seismic anisotropy (> 40%), and the magnitude of S-wave anisotropy highly relates to the density of bedding-parallel cracks (up to 2500 m−1). We present a novel BCD prediction method using the reflected seismic shear wave. The S-wave anisotropy parameter γ is first inverted from the SH-SH wave by using a stepwise inversion method based on a new shear wave reflection coefficient approximation. Furthermore, the highly correlated rock-physics relationship between BCD and the anisotropy parameter γ facilitates accurate BCD estimation. Field data application on a real seismic shear wave data set acquired from the study area demonstrates the viability of the proposed method. The investigation provides a useful tool for the characterization of ‘sweet-spots’ of shale-oil reservoir formation.

Martian seismic anisotropy underneath Elysium Planitia revealed by direct S wave splitting

Seismic anisotropy, arising from the crystallographic/lattice-preferred orientation of anisotropic minerals or the shape-preferred orientation of melts or cracks, can establish a critical link between Mars's past evolution and its current state. So far, seismic anisotropy in Mars has been proposed due to different velocities of vertically and horizontally polarized shear waves in the Martian crust, obtained from crustal converted waves, multiples, and surface waves recorded by the InSight seismometer. However, the shear wave splitting, which stands out as a straightforward indicator of seismic anisotropy, has not been reported using marsquake records. In this study, we employ low-frequency marsquakes detected by the InSight seismometer to reveal shear wave splitting in marsquake recordings. We find that the direct S waves of three marsquake recordings (S0173a, S0235b, and S1133c) with high signal-to-noise ratios exhibit the splitting phenomenon. We rule out the possibility of apparent anisotropy through synthetic tests, affirming the presence of seismic anisotropy in Mars. The delay time measured from the direct S wave splitting is too large to be solely attributed to the seismic anisotropy in the upper crust (0 – 10 km) beneath the InSight. Thus, seismic anisotropy in the deeper region of Mars is required. Combined with other geophysical evidence near the InSight landing site, the seismic anisotropy observed in this study implies the aligned cracks in the crust are greater than 10 km beneath the InSight and/or the existence of mantle flow underneath the Elysium Planitia of Mars.

Crystallographic Preferred Orientation of Phase D at High Pressure and Temperature: Implications for Seismic Anisotropy in the Mid‐Mantle

AbstractSeismic anisotropy has been widely observed in the lower mantle transition zone and the uppermost lower mantle near several subducting slabs, which is typically attributed to the crystallographic preferred orientation (CPO) of constituent minerals. In this study, well‐controlled uniaxial and simple shear deformation experiments were conducted on Mg‐endmember phase D and Al‐bearing phase D aggregates at 20 GPa and 500–1,000°C. Strong (0001) fabrics were observed in the uniaxially compressed samples. The CPO from all experiments, including extension and simple shear, indicated that the most active slip systems in phase D were <110>(0001) and <100>(0001), with a dominance of <110>(0001). The calculated seismic velocity suggests that the deformed phase D generates positive radial seismic anisotropy (VSH > VSV) in both subhorizontal shearing and subvertical compression, causing significant shear‐wave splitting time. Our results suggest that the observed seismic anisotropy in the mid‐mantle in several cold subduction zones can be well explained by the presence of deformed phase D.

Widespread D″ ${\mathbf{D}}^{\mathbf{{\prime\prime}}}$ Anisotropy Beneath North America and the Northeastern Pacific and Implications for Upper Mantle Anisotropy Measurements

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

The tectonic development of the Central African Plateau: evidence from shear-wave splitting

SUMMARY The Central African Plateau comprises a mosaic of numerous Archean terranes—the Congo, Bangweulu and Kalahari Cratons—sutured in a series of Proterozoic to early Cambrian orogenic events. Major upper-crustal deformation and complex craton margin fault zones reflect the region’s diverse tectonic history: rifting during the Neoproterozoic, collision during the Pan-African orogeny, and more recently, Permo-Triassic Karoo rifting and the Pliocene development of the Southwestern branch of the East African Rift. The tectonic evolution and extent to which the lithospheric mantle has been re-worked by each tectonic event is poorly understood. New seismograph networks across the Plateau provide fresh opportunity to place constraints on the plate-scale Precambrian-to-Phanerozoic processes that have acted across the region. Utilizing data from seismograph deployments across the Central African Plateau, including the new Copper Basin Exploration Science network—a NW–SE-trending, 750-km-long profile of 35 broad-band stations—we explore lithospheric deformation fabrics associated with past and present tectonic events via a shear-wave splitting study of mantle seismic anisotropy. Results reveal short length-scale variations in splitting parameters (fast direction: $\phi$, delay time: $\delta$t), suggestive of a fossil lithospheric fabric cause for the observed anisotropy. A lack of fault-parallel $\phi$ across the Mwembeshi shear zone, suggests it may be too narrow at mantle depths, a thin-skinned, crustal-scale feature, and/or did not experience sufficient fault parallel shear-strain during its last active phase to form a lithospheric deformation fabric discernible via teleseismic shear-wave splitting. In the heart of the Lufilian Arc, we observe abrupt changes in splitting parameters with NE–SW, N–S and NW–SE $\phi$ and 0.5 s $< $$\delta$t$< $ 1.2 s evident at short length-scales: no single, uniform, anisotropic lattice preferred orientation (LPO) fabric defines the entire region. This is consistent with the view that multiple episodes of deformation shaped the Lufilian Arc, or perhaps that pre-existing fabrics, relating to Neoproterozoic Katangan Basin development, have failed to be completely overprinted by the Pan-African orogeny. Near the Domes, where most intense crustal re-working is thought to have taken place during the Pan-African orogeny, there is a cluster of null and low $\delta$t splits which likely reflects the lack of organized LPO fabrics, perhaps due to the presence of depth-dependent anisotropy. The neighbouring Congo Craton margin is marked by consistently weak anisotropy ($\delta$t$\lt $ 0.7 s) indicating a weak horizontal alignment of olivine at mantle lithospheric depths, typical of several Archean terranes worldwide.

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.

Calculating traveltimes in 2D general tilted transversely isotropic media using fast sweeping method

Traveltime calculations play an important role in the field of exploration seismology, such as traveltime tomography and seismic imaging and so on. Seismic anisotropy poses a challenge for traveltime calculation, because anisotropic eikonal solvers are more complex than the isotropic counter part. To solve the eikonal equations in 2D tilted transversely isotropic (TTI) media, we have developed a fast algorithm combine with fast sweeping method to compute the first arrival traveltimes of quasi-P (qP)-, quasi-SV (qSV)-, and quasi-SH(qSH)-waves. For the qP- and qSV-waves, we analyzed the quartic coupled slowness surface equation derived from the Christoffel equation. Then, we constructed a local solver to relate traveltime and slowness. We found that in the local solver, one component of the slowness vector is known and the corresponding slowness equation is monotonic. This provides a strong basis for the fast iterative algorithm we proposed, where we use the Newton method to solve the qP- and qSV-wave slowness equation to determine the related traveltimes. For the qSH wave, the slowness equation is quadratic and simple to solve. Numerical experiments demonstrate that the proposed method can obtain accurate traveltimes for simple and complicated 2D TTI models.

Study and application of wide-azimuth seismic anisotropy analysis and correction in shale reservoir in Gulong Sag, Songliao Basin, China

Study and application of wide-azimuth seismic anisotropy analysis and correction in shale reservoir in Gulong Sag, Songliao Basin, China

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.

The sensitivity of lowermost mantle anisotropy to past mantle convection

It is widely believed that seismic anisotropy in the lowermost mantle is caused by the flow-induced alignment of anisotropic crystals such as post-perovskite. What is unclear, however, is whether the anisotropy observations in the lowermost mantle hold information about past mantle flow, or if they only inform us about the present-day flow field. To investigate this, we compare the general and seismic anisotropy calculated using Earth-like mantle convection models where one has a time-varying flow, and another where the present-day flow is constant throughout time. To do this, we track a post-perovskite polycrystal through the flow fields and calculate texture development using the sampled strain rate and the visco-plastic self-consistent approach. We assume dominant slip on (001) and test the effect of the relative importance of this glide plane over others by using three different plasticity models with different efficiencies at developing texture. We compare the radial anisotropy parameters and the anisotropic components of the elastic tensors produced by the flow field test cases at the same location. We find, under all ease-of-texturing cases, the radial anisotropy is very similar (difference <2%) in the majority of locations and in some regions, the difference can be very large (>10%). The same is true when comparing the elastic tensors directly. Varying the ease-of-texture development in the crystal aggregate suggests that easier-to-texture material may hold a stronger signal from past flow than harder-to-texture material. Our results imply that broad-scale observations of seismic anisotropy such as those from seismic tomography, 1-D estimates and normal mode observations, will be mainly sensitive to present-day flow. Shear-wave splitting measurements, however, could hold information about past mantle flow. In general, mantle memory expressed in anisotropy may be dependent on path length in the post-perovskite stability field. Our work implies that, as knowledge of the exact causative mechanism of lowermost mantle anisotropy develops, we may be able to constrain both present-day and past mantle convection.

Wave velocities and anisotropy of rocks: Implication for origin of low velocity zone of the Qinling Orogenic Belt, China

Structure and composition of Earth is of fundamental importance in exploring the dynamic evolution of the crust and mantle. The Qinling Orogenic Belt (QOB) is located between the North China plate and the South China Plate, and is one of the main orogenic belts in China. To explore the composition and origin of anisotropy and the low wave velocity zone of the QOB, ten rock samples (gneiss and schist) were collected from the five sites of the QOB and the P- and S-wave velocities of these samples were measured under 0.6 to 2.0 GPa and 100 to 550 °C. The wave velocities increase with increasing pressure and decreasing temperature. The VP and VS of the schist and gneiss match the velocity of the middle and lower crust of the QOB, indicating that schist and gneiss are important component of the QOB. All the schist and gneiss samples exhibit obvious seismic anisotropy with 1.64%–17.42% for VS and 2.93%–14.78% for VP under conditions of crust and upper mantle. The CPO/LPO and layering distribution of mica in rock samples are the main reasons for this anisotropy. The VS structures below the five sampled sites from seismic ambient noise tomography were built to explore the effect of schist and gneiss on the composition and structure of the QOB. The results indicate that orientation-arranged gneiss and schist driven by the tectonic stresses might be a new origin of the character of VP/VS, seismic anisotropy, and the low velocity zone in the QOB.

Spatial and temporal variations of seismic azimuthal anisotropy following the 2019 ridgecrest earthquake sequence in southern california

To investigate the spatial and temporal variations of seismic azimuthal anisotropy in the vicinity of the conjugate fault network activated by the 4 July 2019 (UTC) M6.4 foreshock and the 6 July 2019 M7.1 mainshock in Ridgecrest, California, a total of 1,470 shear wave splitting measurements were obtained from local events recorded by five seismic stations over the period of July 2019 to January 2020. The results suggest a strong asymmetry in anisotropy forming mechanisms across the NW-SE striking Eastern Little Lake Fault, which is the main fault of the area. In the area located to the northeast of the main fault, the observed fast orientations are dominantly N-S, which is parallel to the maximum horizontal compressive stress. In the area surrounding the main fault, the fast orientations are primarily parallel to the strikes of multiple fault zones including the main fault and the crossing faults, rather than solely consistent with that of the main fault, which may indicate along-strike variations in fault strength of the main fault. To the southwest of the main fault, the observed anisotropy is jointly controlled by faults or regional stress. The observed NE-SW-oriented anisotropy in the vicinity of two previously proposed blind faults confirms the existence of faults. The splitting times of anisotropic areas with different inducements are independent of the focal depths, suggesting that either the stress- or structure-induced anisotropy is mostly located in a shallow layer in the top several kilometers. A nearly 90-degree switch in the fast orientations and greatly reduced splitting times from a group of nearby earthquakes may indicate fault zone healing.

Seismic anisotropy of granitic rocks from a fracture stimulation well at Utah FORGE using ultrasonic measurements

For effective fracture stimulation in Enhanced Geothermal Systems (EGS), characterization of subsurface rock elastic properties is essential. Because seismic anisotropy significantly influences hydraulic fracturing, locating and characterizing microseismic events, and forecasting stress trajectories, our research undertakes laboratory experiments to study seismic anisotropy parameters of core samples from the Utah FORGE EGS site. We evaluated two rock samples from well 58–32, at two depths of 2074.16 m (6,805 feet) and 2268.01 m (7,441 feet), respectively. We analyze the mineralogy and classify the shallower rock as syenogranite, and the deeper rock as orthogneiss and monzonite. We cut cylindrical plugs from these two core samples, imaged them by computed tomography (CT) X-ray imaging, and measured their seismic anisotropy using ultrasonic waves at 1 MHz and under varying confining pressures (5–50 MPa) in ambient temperature conditions. With the assumption of the transversely isotropic (TI) medium, we measured three P-wave velocities and two S-wave velocities along different propagation directions to compute the Thomsen anisotropy parameters. Our results showed that the measured seismic velocities of each rock sample increase with increasing effective pressures, a behavior likely caused by closure of microcracks identified in our CT images. We also find that the measured anisotropies are higher at lower effective pressures. The maximum measured anisotropies for P-waves and S-waves are ∼23% and ∼20%, respectively. Therefore, we expect in the field operation that anisotropy should increase with increasing pore pressure caused by injection, which could be valuable in future EGS reservoir management and monitoring.

Upper-mantle seismic anisotropy in the southwestern North Island, New Zealand: Implications for regional upper-mantle and slab deformation

We employed shear-wave splitting analysis on both teleseismic SKS and S waves, and S waves from deep (150–250 km) local earthquakes collected from a dense array with 43 temporary broadband seismic stations and nine long-term seismic stations centered at Mount Taranaki to characterize the upper-mantle dynamics in the southwestern North Island of New Zealand, in areas previously unexamined for shear-wave splitting. We observed predominantly trench-parallel fast polarizations and strikingly large delay times over 3 s from teleseismic analysis. In contrast, local S analysis yielded a sharp transition of fast-polarization from trench-parallel in the northeast to trench-normal in the southwest. Trench-parallel fast-polarization from teleseismic analysis may be attributed to sub-slab trench-parallel flow or to trench-parallel fractures in the subducting slab. More importantly, we attribute large delay times to deep upper-mantle (200–400 km depth) deformation, possibly associated with the dynamic interaction between the downgoing slab and the 410-km discontinuity or with the lithosphere delamination near the Taranaki-Ruapehu line. In contrast, the trench-parallel anisotropy from the local S waves in the northeast could be caused by fluid-bearing cracks in the crust of the Taupō Volcanic Zone and/or by trench-parallel fractures in the subducting slab resulting from outer rise bending. The abrupt change to trench-normal may be related to stress variations in the downgoing slab at different depths.

Microstructures, hydrogen concentrations, and seismic properties of a tectonically exhumed sliver of oceanic mantle lithosphere, Moa Island, Timor-Tanimbar outer-arc, eastern Indonesia

We characterize and quantify the microstructure, hydrogen concentrations, and seismic properties of a tectonically exhumed sliver of oceanic lithospheric mantle outcropping in the Moa Island (Leti archipelago, Timor-Tanimbar outer-arc). The 18 spinel peridotites (lherzolites and harzburgites) have coarse-porphyroclastic microstructures and olivine crystal-preferred orientations (CPO) with axial-[010] (also known as AG-type) or [100](010) (A-type) patterns, similar to those observed in peridotitic xenoliths from oceanic mantle lithosphere. These coarse-porphyroclastic microstructures are variably overprinted by growth of strain-free olivine neoblasts and crystallization of secondary pyroxenes. Recrystallized fractions vary from 6.9 up to 31.3%. The interstitial (cuspate) shapes and CPOs of clinopyroxene, uncorrelated with the olivine CPOs, indicate that refertilization by a reactive melt percolation post-dated deformation. Seismic properties are calculated based on the modal compositions and CPOs of all samples. Increase in the recrystallized olivine fraction decreased the seismic anisotropy, since static recrystallization produced some dispersion of the CPO, but did not change drastically the texture acquired during deformation. Mean seismic velocities (mean Vp = 7.9 km.s−1; mean Vs = 4.5 km.s−1) and anisotropy (mean maximum S wave polarization anisotropy = 4.5%), estimated by considering coherent orientation of the foliation and lineation of all samples, are within the range of typical values for the uppermost mantle. The nominally anhydrous minerals contain small amounts of hydrogen (olivine: 13–18 ppm H2O by weight; orthopyroxene: 58–175 wt ppm H2O and clinopyroxenes: 244–288 wt ppm H2O). A bulk water content of 50 wt ppm H2O is estimated based on nomminally anhydrous minerals for the Moa peridotites, in agreement with previous estimates for the oceanic mantle lithosphere based on peridotitic xenoliths. This is the first direct measurement of hydrogen concentrations in peridotites from an oceanic mantle lithosphere which experienced melt extraction.