Two-dimensional modelling of subduction zone anisotropy with application to southwestern Japan

We present a series of 2-D numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate. We use a kinematically defined slab geometry approximating the subduction of the Philippine Sea plate beneath Eurasia. Through finite element modelling we explore the effects of different rheological and thermal constraints (e.g. a low-viscosity region in the wedge corner, power law versus Newtonian rheology, the inclusion of thermal buoyancy forces and a temperature-dependent viscosity law) on the velocity and finite strain field in the mantle wedge. From the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material. We then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory, and compare the predictions to ∼500 previously published shear wave splitting measurements from seventeen stations of the broad-band F-net array located in southwestern Japan. Although the use of different model parameters can have a substantial effect on the character of the finite strain field, the effect on the average predicted splitting parameters is small. However, the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different, and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge. The splitting observed in southwestern Japan agrees well with the predictions of trench-perpendicular flow in the mantle wedge along with B-type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench.     

Distribution of seismic anisotropy in the subduction zone beneath the Wellington region, New Zealand

Shear wave splitting measurements from S arrivals of local earthquakes recorded at the Incorporated Research Institutions for Seismology (IRIS) broadband sensor SNZO are used to determine a basic anisotropic structure for the subduction zone in the Wellington region. With the use of high-frequency filters, fast anisotropic polarization ( φ ) and splitting time ( δt ) measurements typical of crustal anisotropy are evident, but the larger splitting expected from the mantle is often not resolved. The small splitting seen agrees well with the results of previous studies concerning shallow crustal anisotropy. With the use of lower-frequency filters, measurements more consistent with mantle anisotropy are made. Anisotropy of 4.4 ± 0.9 per cent with a fast polarization of 29° ± 38° is calculated for the subducting slab, from 20 to 70  km depth. Using this result in addition to the results of previous studies, a model is proposed. The model requires a frequency-dependent anisotropy of less than 1.4 per cent when measured with a period of ~2  s to be present in the sub-slab mantle.
Separate from this population, a band of events in northern Cook Strait with an 86° ± 10° fast polarization is seen. This is at about 40° from the strike of the Hikurangi margin, and suggests a source of shear strain 40° removed from that found in the majority of the region. The cause of this is probably a deformation in the subducting slab in this region, as it moves towards a greater incline to the south.     

Seismic anisotropy beneath Dronning Maud Land, Antarctica, revealed by shear wave splitting

Shear wave splitting analyses have been carried out using teleseismic data from broad-band seismograph stations deployed at temporary and permanent locations in Dronning Maud Land (DML), Antarctica. In most cases, the observed anisotropy can be related to major tectonic events that formed the present-day Antarctic continent. We rule out an anisotropic contribution from recent asthenospheric flow. At the Russian base Novolazarevskaya near the coast in central DML, waveform inversion suggests a two-layer model where the fast direction of the upper layer is oriented parallel to Archean fabrics in the lithosphere, whereas the anisotropy of the lower layer is interpreted to have been created during the Jurassic Gondwana break-up. Recordings at the South African base Sanae IV, however, show enigmatic results. For narrow backazimuthal segments, splitting parameters show strong variations together with a multitude of isotropic measurements, indicative of complex scattering that cannot be explained by simple one- or two-layer anisotropic models. In the interior of the continent, the data are consistent with single-layer anisotropy, but show significant spatial variations in splitting parameters. A set of temporary stations across the Heimefront shear zone in western DML yield splitting directions that we interpret as frozen anisotropy from Proterozoic assembly of the craton. An abrupt change in fast axis direction appears to mark a suture between the Grunehogna craton, a fragment of the Kalahari–Kaapvaal craton in southern Africa and the Mesoproterozoic Maudheim Province.     

Complex upper-mantle deformation beneath the North China Craton: implications for lithospheric thinning

The mechanism of lithospheric thinning of the North China Craton (NCC) remains controversial. To constrain the mechanism, this study investigated the upper-mantle deformation pattern of the craton by measuring shear wave splitting at the cratonic edge. The results, derived from data recorded at 47 stations, reveal a complex pattern of mantle deformation. Inside the eastern craton, the majority of fast direction trends SE–NW parallel to the tectonic extension direction accompanying with the lithospheric thinning. At the cratonic edge, 15 stations with only null splitting results indicate undetectable anisotropy beneath the stations. This may be due to upwelling or chaotic ascension of mantle flow. To the north, off the craton, large delay times and variation of splitting parameter with backazimuth are generated by the combination of lithospheric and asthenospheric anisotropy. Based on comparison of the splitting results and the predicted ones by the compelling models, it is likely that lithospheric delamination dominated the lithospheric thinning at the north edge of the NCC during the Mesozoic to Cenozoic.     

Upper-mantle flow beneath French Polynesia from shear wave splitting

Upper-mantle flow beneath the South Pacific is investigated by analysing shear wave splitting parameters at eight permanent long-period and broad-band seismic stations and 10 broad-band stations deployed in French Polynesia from 2001 to 2005 in the framework of the Polynesian Lithosphere and Upper Mantle Experiment (PLUME). Despite the small number of events and the rather poor backazimuthal coverage due to the geographical distribution of the natural seismicity, upper-mantle seismic anisotropy has been detected at all stations except at Tahiti where two permanent stations with 15 yr of data show an apparent isotropy. The median value of fast polarization azimuths (N67.5°W) is parallel to the present Pacific absolute plate motion direction in French Polynesia (APM: N67°W). This suggests that the observed SKS fast polarization directions result mainly from olivine crystal preferred orientations produced by deformation in the sublithospheric mantle due to viscous entrainment by the moving Pacific Plate and preserved in the lithosphere as the plate cools. However, analysis of individual measurements highlights variations of splitting parameters with event backazimuth that imply an actual upper-mantle structure more complex than a single anisotropic layer with horizontal fast axis. A forward approach shows that a two-layer structure of anisotropy beneath French Polynesia better explains the splitting observations than a single anisotropic layer. Second-order variations in the measurements may also indicate the presence of small-scale lateral heterogeneities. The influence of plumes or fracture zones within the studied area does not appear to dominate the large-scale anisotropy pattern but may explain these second-order splitting variations across the network.