scholarly journals Contact of the Samoan Plume with the Tonga Subduction from Intermediate and Deep-Focus Earthquakes

Author(s):  
Pavla Hrubcová ◽  
Václav Vavryčuk

AbstractThe Tonga subduction zone in the south-west Pacific is the fastest convergent plate boundary in the world with the most active mantle seismicity. This zone shows unique tectonic features including Samoan volcanic lineament of plume-driven origin near the northern rim of the Tonga subducting slab. The proximity of the Samoa hotspot to the slab is enigmatic and invokes debates on interactions between the Samoa plume and the Tonga subduction. Based on long-term observations of intermediate and deep-focus Tonga earthquakes reported in the Global Centroid Moment Tensor (CMT) catalog, we provide novel detailed imaging of this region. Accurate traveltime residua of the P- and S-waves recorded at two nearby seismic stations of the Global Seismographic Network are inverted for the P- and S-wave velocities and their ratio and reveal their pronounced lateral variations. In particular, they differ for the southern and northern parts of the Tonga subduction region. While no distinct anomalies are detected in the southern Tonga segment, striking low-velocity anomalies associated with a high Vp/Vs ratio are observed in the northern Tonga segment close to the Samoa plume. These anomalies spread through the whole upper mantle down to depths of ~ 600 km. Together with the fast extension of the northern back-arc Lau Basin, slab deformation and geochemical enrichment in the northern Tonga region, they trace deep-seated magmatic processes and evidence an interaction of the Tonga subduction with the Samoa plume.

2021 ◽  
Vol 9 ◽  
Author(s):  
Melissa O. Anderson ◽  
Chantal Norris-Julseth ◽  
Kenneth H. Rubin ◽  
Karsten Haase ◽  
Mark D. Hannington ◽  
...  

The transition from subduction to transform motion along horizontal terminations of trenches is associated with tearing of the subducting slab and strike-slip tectonics in the overriding plate. One prominent example is the northern Tonga subduction zone, where abundant strike-slip faulting in the NE Lau back-arc basin is associated with transform motion along the northern plate boundary and asymmetric slab rollback. Here, we address the fundamental question: how does this subduction-transform motion influence the structural and magmatic evolution of the back-arc region? To answer this, we undertake the first comprehensive study of the geology and geodynamics of this region through analyses of morphotectonics (remote-predictive geologic mapping) and fault kinematics interpreted from ship-based multibeam bathymetry and Centroid-Moment Tensor data. Our results highlight two notable features of the NE Lau Basin: 1) the occurrence of widely distributed off-axis volcanism, in contrast to typical ridge-centered back-arc volcanism, and 2) fault kinematics dominated by shallow-crustal strike slip-faulting (rather than normal faulting) extending over ∼120 km from the transform boundary. The orientations of these strike-slip faults are consistent with reactivation of earlier-formed normal faults in a sinistral megashear zone. Notably, two distinct sets of Riedel megashears are identified, indicating a recent counter-clockwise rotation of part of the stress field in the back-arc region closest to the arc. Importantly, the Riedel structures identified in this study directly control the development of complex volcanic-compositional provinces, which are characterized by variably-oriented spreading centers, off-axis volcanic ridges, extensive lava flows, and point-source rear-arc volcanoes. This study adds to our understanding of the geologic and structural evolution of modern backarc systems, including the association between subduction-transform motions and the siting and style of seafloor volcanism.


2021 ◽  
Author(s):  
JD Eccles ◽  
AK Gulley ◽  
PE Malin ◽  
CM Boese ◽  
John Townend ◽  
...  

© 2015. American Geophysical Union. All Rights Reserved. Fault Zone Guided Waves (FZGWs) have been observed for the first time within New Zealand's transpressional continental plate boundary, the Alpine Fault, which is late in its typical seismic cycle. Ongoing study of these phases provides the opportunity to monitor interseismic conditions in the fault zone. Distinctive dispersive seismic codas (~7-35Hz) have been recorded on shallow borehole seismometers installed within 20m of the principal slip zone. Near the central Alpine Fault, known for low background seismicity, FZGW-generating microseismic events are located beyond the catchment-scale partitioning of the fault indicating lateral connectivity of the low-velocity zone immediately below the near-surface segmentation. Initial modeling of the low-velocity zone indicates a waveguide width of 60-200m with a 10-40% reduction in S wave velocity, similar to that inferred for the fault core of other mature plate boundary faults such as the San Andreas and North Anatolian Faults.


2019 ◽  
Vol 219 (2) ◽  
pp. 1447-1462 ◽  
Author(s):  
Alexandre P Plourde ◽  
Michael G Bostock

SUMMARY We introduce a new relative moment tensor (MT) inversion method for clusters of nearby earthquakes. The method extends previous work by introducing constraints from S-waves that do not require modal decomposition and by employing principal component analysis to produce robust estimates of excitation. At each receiver, P and S waves from each event are independently aligned and decomposed into principal components. P-wave constraints on MTs are obtained from a ratio of coefficients corresponding to the first principal component, equivalent to a relative amplitude. For S waves we produce constraints on MTs involving three events, where one event is described as a linear combination of the other two, and coefficients are derived from the first two principal components. Nonlinear optimization is applied to efficiently find best-fitting tensile-earthquake and double-couple solutions for relative MT systems. Using synthetic data, we demonstrate the effectiveness of the P and S constraints both individually and in combination. We then apply the relative MT inversion to a set of 16 earthquakes from southern Alaska, at ∼125 km depth within the subducted Yakutat terrane. Most events are compatible with a stress tensor dominated by downdip tension, however, we observe several pairs of earthquakes with nearly antiparallel slip implying that the stress regime is heterogeneous and/or faults are extremely weak. The location of these events near the abrupt downdip termination of seismicity and the low-velocity zone suggest that they are caused by weakening via grain-size and volume reduction associated with eclogitization of the lower crustal gabbro layer.


2021 ◽  
Vol 873 (1) ◽  
pp. 012043
Author(s):  
Jaya Murjaya ◽  
Pepen Supendi ◽  
Dwikorita Karnawati ◽  
Subagyo Pramumijoyo

Abstract During the last one hundred years, there are no shallow seismicity in the north of Java. This area is dominated by intermediate and deep focus earthquakes due to the subducted Indo-Australian slab. An earthquake with magnitude ML 4.5 struck Indramayu, north of West Java on August 1, 2020. According to the Agency for Meteorology, Climatology, and Geophysics (BMKG), the earthquake was felt III MMI scale in Indramayu and its vicinity. We used waveform data from BMKG seismic station in West Java, then we picked P-and S-waves arrival times from each station and hypocenter location was determined by Geiger method. We have detected Pn before Pg phase on four BMKG seismic stations, indicating a shallow crustal earthquake. Our inversion show that the earthquake occurred in 6.1805° S, 108.2612° E with 5 km focus depth at 16:24:38 GMT+7. Our focal mechanism solution was determined by using moment tensor inversion shows a strike-slip faulting, which corresponds to the active fault in the north of Indramayu.


2020 ◽  
Author(s):  
Simona Gabrielli ◽  
Luca De Siena ◽  
Matteo Spagnolo

<p>In volcanoes, topography, shallow heterogeneity, and even shallow morphology can substantially modify seismic coda signals. Coda waves are an essential tool to monitor eruption dynamics and model volcanic structures jointly and independently from velocity anomalies: it is thus fundamental to test their spatial sensitivity to seismic path effects. Here, we apply the Multiple Lapse Time Window Analysis (MLTWA) to measure the relative importance of scattering attenuation vs absorption at Mount St. Helens volcano (MSH) before its 2004 eruption. The results show the typical dominance of scattering attenuation in volcanoes at lower frequencies (3 - 6 Hz), while absorption is the primary attenuation mechanism at 12 Hz and 18 Hz. Still, the seismic albedo (measuring the ratio between seismic energy emitted and received from the area) is anomalously-high (0.95) at 3 Hz.</p><p>A radiative-transfer forward model of far- and near-field envelopes confirms this is due to strong near-receiver scattering enhancing anomalous phases in the intermediate and late coda across the 1980 debris avalanche and central crater. Only above this frequency and in the far-field, diffusion onsets at late lapse times.  We also implemented a layered model with a shallower layer with increased scattering properties to model late coda envelopes. While the broadening of late coda phases improves, this model cannot explain the phases of the intermediate coda with higher amplitude than the direct waves.</p><p>The scattering and absorption parameters derived from MLTWA are used as inputs to construct 2D frequency-dependent bulk sensitivity kernels for the S-wave coda in the multiple-scattering (using the Energy Transport Equations - ETE) and diffusive (AD, independent of MLTWA results) regimes. At 12 Hz, high coda-attenuation anomalies characterise the eastern side of the volcano using both kernels, in spatial correlation with low-velocity anomalies from literature. At 3 Hz, the anomalous albedo, the forward modelling, and the results of the tomographic imaging confirm that shallow heterogeneity beneath the extended 1980 debris-avalanche and crater enhance anomalous intermediate and late coda phases, mapping shallow geological contrasts.</p><p>The geomorphological map of MSH highlights extremely rough landforms (hummocky structures) of the already complex morphology of the debris avalanche. The comparison with the attenuation tomography reveals several matches, not only with the debris avalanche itself but also with other areas in the south flank of MSH, like the volcanoclastic plane, affected by intense eruptions in the past (e.g. Cougar stage, 28-18 ka).</p><p>We remark the effect those may have on coda-dependent source inversion and tomography, currently used across the world to image and monitor volcanoes.</p>


2015 ◽  
Vol 7 (3) ◽  
pp. 1909-1939
Author(s):  
M. L. Kolstrup ◽  
V. Maupin

Abstract. We present a data processing routine to compute relative finite-frequency travel time residuals using a combination of the Iterative Cross-Correlation and Stack (ICCS) algorithm and the MultiChannel Cross-Correlation method (MCCC). The routine has been tailored for robust measurement of P and S wave travel times in several frequency bands and for avoiding cycle-skipping problems at the shortest periods. We also investigate the adequacy of ray theory to calculate crustal corrections for finite-frequency regional tomography in normal continental settings with non-thinned crust. We find that ray theory is valid for both P and S waves at all relevant frequencies as long as the crust does not contain low-velocity layers associated with sediments at the surface. Reverberations in the sediments perturb the arrival times of the S waves and the long-period P waves significantly, and need to be accounted for in crustal corrections. The data processing routine and crustal corrections are illustated using data from a network in southwestern Scandinavia.


1969 ◽  
Vol 59 (1) ◽  
pp. 385-398 ◽  
Author(s):  
Otto W. Nuttli

Abstract The underground Nevada explosions HALF-BEAK and GREELEY were unique in creating relatively large amplitude and long-period body S waves which could be detected at teleseismic distances. Observations of the travel times of these S waves provide a surface focus travel-time curve which in its major features is similar to a curve calculated from the upper mantle velocity model of Ibrahim and Nuttli (1967). This model includes a low-velocity channel at a depth of 150 to 200 km and regions of rapidly increasing velocity beginning at depths of 400 and 750 km. Observations of the S wave amplitudes suggest that a discontinuous increase in velocity occurs at 400 km, whereas at 750 km the velocity is continuous but the velocity gradient discontinuous. Body wave magnitudes calculated from S amplitudes are 5.3 ± 0.2 for GREELEY and 4.9 ± 0.2 for HALF-BEAK. These are about one unit less than body wave magnitudes from P amplitudes as reported by others. The shape and orientation of the radiation pattern of SH for both explosions are consistent with the Rayleigh and P-wave amplitude distribution of BILBY as given by Toksoz and Clermont (1967). This suggests that the regional stress field is the same at all three sites, and that the direction of cracking as well as the strain energy release in the elastic zone outside the cavity is determined by the regional stress field.


Solid Earth ◽  
2015 ◽  
Vol 6 (4) ◽  
pp. 1117-1130 ◽  
Author(s):  
M. L. Kolstrup ◽  
V. Maupin

Abstract. We present a data-processing routine to compute relative finite-frequency travel time residuals using a combination of the Iterative Cross-Correlation and Stack (ICCS) algorithm and the Multi-Channel Cross-Correlation method (MCCC). The routine has been tailored for robust measurement of P- and S-wave travel times in several frequency bands and for avoiding cycle-skipping problems at the shortest periods. We also investigate the adequacy of ray theory to calculate crustal corrections for finite-frequency regional tomography in normal continental settings with non-thinned crust. We find that ray theory is valid for both P and S waves at all relevant frequencies as long as the crust does not contain low-velocity layers associated with sediments at the surface. Reverberations in the sediments perturb the arrival times of the S waves and the long-period P waves significantly, and need to be accounted for in crustal corrections. The data-processing routine and crustal corrections are illustrated using data from a~network in southwestern Scandinavia.


2020 ◽  
Vol 222 (2) ◽  
pp. 1109-1125 ◽  
Author(s):  
Shunsuke Takemura ◽  
Ryo Okuwaki ◽  
Tatsuya Kubota ◽  
Katsuhiko Shiomi ◽  
Takeshi Kimura ◽  
...  

SUMMARY Due to complex 3-D heterogeneous structures, conventional 1-D analysis techniques using onshore seismograms can yield incorrect estimation of earthquake source parameters, especially dip angles and centroid depths of offshore earthquakes. Combining long-term onshore seismic observations and numerical simulations of seismic wave propagation in a 3-D model, we conducted centroid moment tensor (CMT) inversions of earthquakes along the Nankai Trough between April 2004 and August 2019 to evaluate decade-scale seismicity. Green's functions for CMT inversions of earthquakes with moment magnitudes of 4.3–6.5 were evaluated using finite-difference method simulations of seismic wave propagation in the regional 3-D velocity structure model. Significant differences of focal mechanisms and centroid depths between previous 1-D and our 3-D catalogues were found in the solutions of offshore earthquakes. By introducing the 3-D structures of the low-velocity accretionary prism and the Philippine Sea Plate, dip angles and centroid depths for offshore earthquakes were well-constrained. Teleseismic CMT also provides robust solutions, but our regional 3-D CMT could provide better constraints of dip angles. Our 3-D CMT catalogue and published slow earthquake catalogues depicted spatial distributions of slip behaviours on the plate boundary along the Nankai Trough. The regular and slow interplate earthquakes were separately distributed, with these distributions reflecting the heterogeneous distribution of effective strengths along the Nankai Trough plate boundary. By comparing the spatial distribution of seismic slip on the plate boundary with the slip-deficit rate distribution, regions with strong coupling were clearly identified.


Sign in / Sign up

Export Citation Format

Share Document