S-wave splitting in the transpressional zone of New Zealand plate-boundary: implications for deformation and dynamics

<p>Seismic anisotropy in the transpressional plate-boundary zone in New Zealand is investigated with shear-wave splitting to gain insights into lithospheric deformation and mantle flow. Constraints on plate-boundary deformation in the lithosphere of the oblique-collision and subduction regimes in South Island have been estimated from the local and regional shear-wave splitting parameters that are made at both inland and offshore seismographs. Mantle and lithospheric anisotropy of the southernmost Hikurangi subduction zone in the southern North Island is examined from SKS, ScS and teleseismic S-phases. The splitting of these phases measured on a recent transect crossing the Wellington region is analyzed to understand the lateral anisotropic structure of the fore-arc Hikurangi subduction zone. Local and regional splitting reveal both laterally and depth varying anisotropy in South Island. The scatter in splitting parameters at individual stations suggests the splitting of high-frequency S-phases is mainly controlled by heterogeneous anisotropic structure and S-wave propagation direction within those heterogeneities. When the average results are examined as a whole through 2-D delay time tomographic inversion and spatial averaging, consistent patterns in delay times and fast azimuths exist. Spatially averaged fast azimuths indicate a localized high strain zone in the southern central region of the South Island. Based on fast azimuths observed above 100 km depth, we suggest that the plate-boundary sub-parallel anisotropy that is produced by pervasive shear is mainly distributed within a zone extending ~130 km SE of the Alpine fault in the southern South Island and is widely distributed (at least 200 km wide) in the northern South Island. Average station delay times (δt) of ~0.1 - 0.4 s compared to 1.7 s SKS δt from previous studies in South Island further suggest a deep seated anisotropic zone or sensitivity of S-wave splitting to the layered and/or heterogeneous anisotropic structure of the plate-boundary zone in the inland South Island. The heterogeneous anisotropic structure further suggests that the lithosphere is not only characterized by the plate-boundary parallel shear related to Cenozoic deformation, but is also affected by anisotropic imprints from the other tectonic episodes and anisotropy that is governed by the contemporary stress. A shear-wave splitting anisotropy investigation in the offshore South Island regions is an extended study of the inland experiment and aims to provide a broad-scale understanding of the plate-boundary deformation. Individual splitting of local and regional S-phases yield a range of δt that varies between very small δt (~0.02 s), which may represent a nearly isotropic medium, and large δt (~0.6 s), which corresponds to lithospheric anisotropy. The average station δt of ~0.25 s and variable delays of the individual splitting measurements imply that the observed splitting is most likely controlled by the geometry of the ray paths. Long ray paths that are detected at the stations further away from the plate-boundary appear to penetrate to deeper lithosphere and capture a significant portion of the upper-mantle anisotropy to produce fast azimuths parallel to the plate-boundary shear (NE-SW). Thus, the long and deep ray paths respond to the deeper structure, but may not be re-split by the upper-most crustal structures. However, the observed variable delays suggest that changes in ray propagation direction with respect to the orientation of symmetry axes of the anisotropic media may have an effect on the measured anisotropy. Offshore measurements that are close to the land are consistent with the inland measurements and appear to be controlled by the regional stress field. This implies that short and shallow ray paths are mostly sensitive to the crustal anisotropy. The uneven distribution of ray paths from the shallow and deep events, therefore, plays a dominant role in controlling the observed splitting depending on their depth sensitivity and/or extent of anisotropy. Consequently, when fast directions are spatially averaged along with the inland measurements consistent patterns appear to correlate with the possible depth contribution of anisotropy in the region. We are unable to provide accurate constraints on the offshore extent of plate-boundary parallel shear because of the shallow stress-controlled anisotropy that likely overlies the mantle-shear zone. However, the splitting parameters from long and deep ray paths suggest a deep-seated, plate-boundary sub-parallel shear in a broad zone at least in the northern and upper-central South Island. Mantle anisotropy detected from teleseismic earthquakes recorded across the southern North Island displays NE-SW fast axis alignment, consistent with the strike of the Hikurangi trench and the predominant upper-plate faulting trends, with a range of δt (~0.5 - 3.0 s) and small-scale variation in NE-SW fast azimuths. When combined with the previous measurements in the western side of the array, δt from long period (>7 s) S-phases indicate an abrupt lateral variation across the fore-arc Hikurangi subduction zone. This lateral variation together with frequency dependence suggest that the shear wave splitting in the fore-arc of the Hikurangi subduction zone in the southern North Island is governed in part by the laterally varying crustal contribution of anisotropy or isotropic velocity variations within the shallow crust. Frequency dependent splitting also suggests that the anisotropic structure is governed by either multilayer or more complex anisotropy perhaps due to the combined effects of laterally varying multilayer structure. If the variations are due to lateral changes in crustal anisotropy, then mantle and crustal deformation are most likely coupled in the east of the Wairarapa fault where there is a possibility of strong crustal contribution.</p>

S-wave splitting in the transpressional zone of New Zealand plate-boundary: implications for deformation and dynamics

<p>Seismic anisotropy in the transpressional plate-boundary zone in New Zealand is investigated with shear-wave splitting to gain insights into lithospheric deformation and mantle flow. Constraints on plate-boundary deformation in the lithosphere of the oblique-collision and subduction regimes in South Island have been estimated from the local and regional shear-wave splitting parameters that are made at both inland and offshore seismographs. Mantle and lithospheric anisotropy of the southernmost Hikurangi subduction zone in the southern North Island is examined from SKS, ScS and teleseismic S-phases. The splitting of these phases measured on a recent transect crossing the Wellington region is analyzed to understand the lateral anisotropic structure of the fore-arc Hikurangi subduction zone. Local and regional splitting reveal both laterally and depth varying anisotropy in South Island. The scatter in splitting parameters at individual stations suggests the splitting of high-frequency S-phases is mainly controlled by heterogeneous anisotropic structure and S-wave propagation direction within those heterogeneities. When the average results are examined as a whole through 2-D delay time tomographic inversion and spatial averaging, consistent patterns in delay times and fast azimuths exist. Spatially averaged fast azimuths indicate a localized high strain zone in the southern central region of the South Island. Based on fast azimuths observed above 100 km depth, we suggest that the plate-boundary sub-parallel anisotropy that is produced by pervasive shear is mainly distributed within a zone extending ~130 km SE of the Alpine fault in the southern South Island and is widely distributed (at least 200 km wide) in the northern South Island. Average station delay times (δt) of ~0.1 - 0.4 s compared to 1.7 s SKS δt from previous studies in South Island further suggest a deep seated anisotropic zone or sensitivity of S-wave splitting to the layered and/or heterogeneous anisotropic structure of the plate-boundary zone in the inland South Island. The heterogeneous anisotropic structure further suggests that the lithosphere is not only characterized by the plate-boundary parallel shear related to Cenozoic deformation, but is also affected by anisotropic imprints from the other tectonic episodes and anisotropy that is governed by the contemporary stress. A shear-wave splitting anisotropy investigation in the offshore South Island regions is an extended study of the inland experiment and aims to provide a broad-scale understanding of the plate-boundary deformation. Individual splitting of local and regional S-phases yield a range of δt that varies between very small δt (~0.02 s), which may represent a nearly isotropic medium, and large δt (~0.6 s), which corresponds to lithospheric anisotropy. The average station δt of ~0.25 s and variable delays of the individual splitting measurements imply that the observed splitting is most likely controlled by the geometry of the ray paths. Long ray paths that are detected at the stations further away from the plate-boundary appear to penetrate to deeper lithosphere and capture a significant portion of the upper-mantle anisotropy to produce fast azimuths parallel to the plate-boundary shear (NE-SW). Thus, the long and deep ray paths respond to the deeper structure, but may not be re-split by the upper-most crustal structures. However, the observed variable delays suggest that changes in ray propagation direction with respect to the orientation of symmetry axes of the anisotropic media may have an effect on the measured anisotropy. Offshore measurements that are close to the land are consistent with the inland measurements and appear to be controlled by the regional stress field. This implies that short and shallow ray paths are mostly sensitive to the crustal anisotropy. The uneven distribution of ray paths from the shallow and deep events, therefore, plays a dominant role in controlling the observed splitting depending on their depth sensitivity and/or extent of anisotropy. Consequently, when fast directions are spatially averaged along with the inland measurements consistent patterns appear to correlate with the possible depth contribution of anisotropy in the region. We are unable to provide accurate constraints on the offshore extent of plate-boundary parallel shear because of the shallow stress-controlled anisotropy that likely overlies the mantle-shear zone. However, the splitting parameters from long and deep ray paths suggest a deep-seated, plate-boundary sub-parallel shear in a broad zone at least in the northern and upper-central South Island. Mantle anisotropy detected from teleseismic earthquakes recorded across the southern North Island displays NE-SW fast axis alignment, consistent with the strike of the Hikurangi trench and the predominant upper-plate faulting trends, with a range of δt (~0.5 - 3.0 s) and small-scale variation in NE-SW fast azimuths. When combined with the previous measurements in the western side of the array, δt from long period (>7 s) S-phases indicate an abrupt lateral variation across the fore-arc Hikurangi subduction zone. This lateral variation together with frequency dependence suggest that the shear wave splitting in the fore-arc of the Hikurangi subduction zone in the southern North Island is governed in part by the laterally varying crustal contribution of anisotropy or isotropic velocity variations within the shallow crust. Frequency dependent splitting also suggests that the anisotropic structure is governed by either multilayer or more complex anisotropy perhaps due to the combined effects of laterally varying multilayer structure. If the variations are due to lateral changes in crustal anisotropy, then mantle and crustal deformation are most likely coupled in the east of the Wairarapa fault where there is a possibility of strong crustal contribution.</p>

Ray paths of VLF/LF transmitter radio signals in the seismic Adriatic regions

&lt;p&gt;We analyze the radio wave propagations of VLF/LF transmitter signals along subionospheric paths using two different reception systems localized in the Graz seismo-electromagnetic facility (15.43E,47.06N). Those systems allow the simultaneous detection of more than fifteen transmitter signals emitting in the northern (i.e. France, Germany and United Kingdom) and southern (i.e. Italy and Turkey) parts of Europe. In this work, we investigate the transmitter radio wave propagations associated with two earthquakes (EQs) which occurred, at two occasions, in nearly the same Croatian regions (Geo. Long.=16&amp;#176;E; Geo. Lat.=45&amp;#176;N). The first and second EQs happened, respectively, on March 22 and December 29, 2020, with magnitudes M&lt;sub&gt;w&lt;/sub&gt; equal to 5.4 and 6.4. The use of two complementary reception systems, i.e. INFREP (Biagi et al., Open Journal of Earthquake Research, 8, 2019) and UltraMSK (Schwingenschuh et al., Nat. Hazards Earth Syst. Sci., 11, 2011), and the proximity to the epicenters lead us to characterize the behavior of the transmitter signal amplitudes particularly above the Croatian seismic regions. We analyze the amplitude variation for a given transmitter frequency starting few weeks before the earthquakes occurrences. We discuss the observed anomalies in the transmitter signals which may be considered as precursors due to the ionospheric disturbances of the transmitter ray paths above the earthquakes preparation zones.&amp;#160;&lt;/p&gt;

Uppermost Mantle Velocity beneath the Mid-Atlantic Ridge and Transform Faults in the Equatorial Atlantic Ocean

ABSTRACT Seismic rays traveling just below the Moho provide insights into the thermal and compositional properties of the upper mantle and can be detected as Pn phases from regional earthquakes. Such phases are routinely identified in the continents, but in the oceans, detection of Pn phases is limited by a lack of long-term instrument deployments. We present estimates of upper-mantle velocity in the equatorial Atlantic Ocean from Pn arrivals beneath, and flanking, the Mid-Atlantic Ridge and across several transform faults. We analyzed waveforms from 50 earthquakes with magnitude Mw&gt;3.5, recorded over 12 months in 2012–2013 by five autonomous hydrophones and a broadband seismograph located on the St. Peter and St. Paul archipelago. The resulting catalog of 152 ray paths allows us to resolve spatial variations in upper-mantle velocities, which are consistent with estimates from nearby wide-angle seismic experiments. We find relatively high velocities near the St. Paul transform system (∼8.4 km s−1), compared with lower ridge-parallel velocities (∼7.7 km s−1). Hence, this method is able to resolve ridge-transform scale velocity variations. Ray paths in the lithosphere younger than 10 Ma have mean velocities of 7.9±0.5 km s−1, which is slightly lower than those sampled in the lithosphere older than 20 Ma (8.1 km±0.3 s−1). There is no apparent systematic relationship between velocity and ray azimuth, which could be due to a thickened lithosphere or complex mantle upwelling, although uncertainties in our velocity estimates may obscure such patterns. We also do not find any correlation between Pn velocity and shear-wave speeds from the global SL2013sv model at depths &lt;150 km. Our results demonstrate that data from long-term deployments of autonomous hydrophones can be used to obtain rare and insightful estimates of uppermost mantle velocities over hundreds of kilometers in otherwise inaccessible parts of the deep oceans.

Supplemental Material: Seismic velocity structure of the Queen Charlotte terrace off western Canada in the region of the 2012 Haida Gwaii Mw 7.8 thrust earthquake

Modeled ray paths and comparisons to picked arrivals are shown here for the refractions in sediment layers 1 (Fig. S1) and 2 (Fig.S2), the reflections from the bottoms of sediment layers 1 (Fig. S3) and 2 (Fig. S4), the refraction from the bottom of the crustal velocity layer (Fig. S5), and the head-wave from the bottom of the crustal-velocity layer (Fig. S6).

Supplemental Material: Seismic velocity structure of the Queen Charlotte terrace off western Canada in the region of the 2012 Haida Gwaii Mw 7.8 thrust earthquake

Modeled ray paths and comparisons to picked arrivals are shown here for the refractions in sediment layers 1 (Fig. S1) and 2 (Fig.S2), the reflections from the bottoms of sediment layers 1 (Fig. S3) and 2 (Fig. S4), the refraction from the bottom of the crustal velocity layer (Fig. S5), and the head-wave from the bottom of the crustal-velocity layer (Fig. S6).

Supplemental Material: Hotspot signatures at the North American passive margin

Methods and additional figures, including ray paths, phase velocities, model resolution tests, and cross plots of shear wave velocity versus radial anisotropy.<br>

Supplemental Material: Hotspot signatures at the North American passive margin

Methods and additional figures, including ray paths, phase velocities, model resolution tests, and cross plots of shear wave velocity versus radial anisotropy.<br>

Supplemental Material: Seismic velocity structure of the Queen Charlotte terrace off western Canada in the region of the 2012 Haida Gwaii Mw 7.8 thrust earthquake

Modeled ray paths and comparisons to picked arrivals are shown here for the refractions in sediment layers 1 (Fig. S1) and 2 (Fig.S2), the reflections from the bottoms of sediment layers 1 (Fig. S3) and 2 (Fig. S4), the refraction from the bottom of the crustal velocity layer (Fig. S5), and the head-wave from the bottom of the crustal-velocity layer (Fig. S6).