Plate Boundary Research Articles

Low‐Frequency Earthquakes Downdip of Deep Slow Slip Beneath the North Island of New Zealand

AbstractWe report the first catalog of low‐frequency earthquakes in the Hikurangi subduction zone, located beneath the Kaimanawa Range of the North Island at 50 km depth, downdip of regularly recurring (every 4–5 years) deep M7 slow slip events. To systematically detect low‐frequency earthquakes within the regional continuous seismic data, we utilized a matched‐filter approach with template waveforms derived from previous observations of tectonic tremor. We built our catalog of 36 low‐frequency earthquake sources, that produced almost 21,000 events over more than a decade, with two matched‐filter search iterations. In each iteration, the detections were gathered into families and their coherent waveforms processed and stacked to extract high‐quality waveforms, allowing us to pick seismic phase arrivals to locate the low‐frequency earthquakes. We highlight three characteristic features to validate that our detected events are indeed low‐frequency earthquakes: the eponymous deficit of high frequencies in their seismic waveforms, the episodic swarms of activity that define their activity through time, and their location at the plate boundary with a double‐couple source mechanism and geometry consistent with the subduction interface. Considering the observed low‐frequency earthquakes' relationship to neighboring slow slip, we observe the event swarms to occur much more frequently than the M7 slow slip events located just updip. Similar to other deep low‐frequency earthquakes in other subduction zones, we suggest that this characteristic clustering in time is driven by more frequent, smaller slow slip events that are not clearly observable at the surface.

Incipient Subduction and Slip Partitioning at High Obliquity: The Haida Gwaii Plate Boundary

AbstractPlate motion obliquity along the dominantly transform Queen Charlotte plate boundary (QCPB) peaks offshore Haida Gwaii. To investigate the effects of obliquity on plate boundary deformation, we analyze continuous seismic waveforms from temporary and permanent stations from 1998 to 2020 to generate a catalog of ∼50,000 earthquakes across Haida Gwaii. We use an automated technique based on auto‐regressive phase detection and onset estimation to obtain the initial seismic catalog, integrate existing catalogs, invert for 3D velocity structure using data from the best constrained period, and relocate the entire catalog using the new 3D velocity model. We investigate the seismically active sections of the transcurrent Queen Charlotte fault (QCF), noting that little seismicity locates directly along its bathymetrically defined trace. Instead, seismicity illuminates a complex system of segmented structures with variable geometries along strike. Other clusters highlight active shallow faults within the highly deformed Queen Charlotte terrace. Few aftershocks appear on the thrust plane of the 2012 Mw 7.8 Haida Gwaii earthquake except near its inferred intersection with the QCF at 15–20 km depths, suggesting elevated residual stress at the juncture of slip‐partitioning. Deep crustal seismicity (up to ∼20 km depths) beneath central Haida Gwaii aligned parallel to the strike of the thrust plane may represent landward underthrusting of the Pacific plate. Our results suggest possible coseismic strike‐slip rupture on the QCF during the 2012 earthquake and add support to the thesis that highly oblique transform boundaries are viable settings for subduction initiation.

Mantle Flow Induced by the Interplay of Downgoing Slabs Revealed by Seismic Anisotropy Beneath the Sula Block in Eastern Indonesia

AbstractThe North Sulawesi subduction zone is characterized by southward subduction of the Celebes Sea slab to a depth of ∼250 km, mainly overlying the Sangihe slab that subducts west from the Molucca Sea and penetrates the mantle transition zone. The palaeo‐subducted Sula slab dips northward and partially underlies both the Sangihe and Celebes Sea slabs. Adjacent subduction zones with horizontal overlapping subducting slabs in the upper mantle have unclear dynamic interactions. An extensive strike‐slip fault forms the western boundary of the active North Sulawesi subduction zone, providing an ideal setting to study mantle flow between overlapping slabs. We use local S‐wave and teleseismic S and SK(K)S waveform splitting analysis to measure seismic anisotropy in the northern Sulawesi region. Our observations reveal typical mantle wedge corner flow within the Sangihe subduction system. In the Gulf of Tomini, the observed trench‐oblique fast‐axis orientations above the Celebes Sea slab are likely a consequence of the interaction between two subducting slabs. The southernmost measurement with an E–W‐trending fast direction in the mantle wedge might be related to the subduction of the Sula slab. Furthermore, fault‐parallel fast‐axis orientations of anisotropy near the southern segment of the Palu‐Koro fault are attributed to large‐scale shearing across this lithospheric‐scale strike‐slip fault system. Overall, our observations suggest that the strain caused by lithospheric and asthenospheric deformation is mainly confined within the microplate, displaying a restricted flow pattern and localized effects due to the size of the plate boundaries, such as the Palu‐Koro fault.

Testing Megathrust Rupture Models Using Tsunami Deposits

AbstractThe 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D‐FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

Plate interface geometry complexity and persistent heterogenous coupling revealed by a high-resolution earthquake focal mechanism catalog in Mentawai, Sumatra

Geometry and frictional states of the megathrust in subduction zones play critical roles on the rupture extent, hence the magnitude, of great earthquakes. However, their details are often difficult to obtain, primarily due to limited offshore geophysical observations. Here we overcome this challenge by precisely determining source parameters (i.e., location and focal mechanism) of medium-sized earthquakes using global broadband seismic waveform data. The dip angles of 108 Mw 5.0+ plate interface events from 1990 to 2017 in the central Sumatran subduction zone are tightly constrained by waveform inversion, and their locations are refined with surface-wave relocation and depth phases modeling. Our results reveal compact trench parallel earthquake belts at the up-dip and down-dip boundaries of the ruptures of the Mentawai 2007 Mw 8.4 and Mw 7.9 earthquakes. The down-dip seismicity belt marks ∼8º bending of the plate interface, suggesting that the change of fault geometry may have played a role in halting the megathrust rupture along with thermal effects. The up-dip seismicity belt shows a clear spatial complementarity to the largest events’ slip distribution. This belt is quite narrow (∼20 km) for the Mw 8.4 rupture but much broader (∼50 km) for the Mw 7.9 event. The combination of the up-dip seismicity belt and slip distributions matches well with the positive residual bathymetry and gravity anomalies, suggesting that these events resulted from long-term coupling on the plate boundary. Our findings suggest the importance of considering both up-dip and down-dip seismicity belts, along with residual bathymetry and gravity data, in understanding megathrust rupture behaviors and assessing the associated hazards.

Stress Variations in Southern Tonga Slab Derived From Deep‐Focus Earthquakes

AbstractTonga is a convergent plate boundary between the Pacific and Australian plates and is the fastest and the most seismically active deep subduction system in the world. We focused on southern Tonga (south of latitude 22°S) and the mantle transition zone (depths of 410–670 km), where seismic activity forms two subparallel bands of events in the east and west. We performed stress analysis by inverting focal mechanisms of earthquakes available in the Global Centroid Moment Tensor catalog and revealed two distinct stress regimes in the slab. While the stress orientation in the eastern slab segment conforms to the down‐dip compressional stress along the entire slab, the stress orientation in the western slab segment is different, having the maximum compression in the vertical direction. This suggests that the western segment can represent a stagnant slab with flattening and bending, as proposed by modeling studies. Its connection with the younger actively subducting slab is supported by the horizontal westward shift at 520 km depth. The stress analysis also indicates substantially different fault orientations in both segments. In the actively dipping slab, the majority of activated faults are predominantly sub‐horizontal. However, they are significantly inclined from vertical in the stagnant slab segment. A higher scatter in fault orientations in the stagnant slab suggests deformation, fragmentation and rheological complexity resulting from bending and flattening.

Two-stage exhumation of high-pressure rocks by corner flow and mud volcanism within an active subduction zone – A case study from serpentinite mud volcanoes along the Mariana convergent margin

During International Ocean Discovery Program (IODP) Expedition 366 cores were recovered from three serpentinite mud volcanoes containing clasts that originate from the subduction-channel along the Philippine Sea Plate – Pacific Plate boundary. The drilled and sampled mud volcanoes (Yinazao, Fantangisña, and Asùt Tesoru) are located at distances of 55 to 72 km from the Mariana Trench. Mafic rock clasts, embedded within a serpentinite mud matrix, from the flanks and summits of both Asùt Tesoru and Fantangisña Seamounts were analyzed for reconstruction of their metamorphic and deformational overprint in order to reveal the tectono-metamorphic conditions at the metamorphic peak within the subduction channel and the subsequent low-grade overprint during exhumation. Several seamounts comprise clasts of lower plate metabasites with different metamorphic overprint (from low-grade sub-greenschist facies to lower blueschist facies). The metabasites are also associated with clasts of fossiliferous carbonates and cherts with different degrees of metamorphic and deformational overprint, that also originated from the Pacific lower Plate. This implies that these rocks were exhumed from different depths within the subduction channel before being regurgitated within a serpentinite mud matrix. The blueschist facies metamorphic rocks, being affected by metamorphic pressures in the range of 11 to 13.8 kbar at minimum, were very likely exhumed from greater depth within the subduction channel before being captured by uprising, localized serpentinite mudflows, indicating evidence that corner flow is actually taking place along the Mariana convergent margin, and that the high-pressure rocks were exhumed by corner flow in an active subduction zone. Final exhumation, however, is related to the embedding of the rocks within a serpentinite mud matrix and the buoyant ascent of serpentinite mudflows along forearc fracture zones extending from the plate boundary to the upper plate sea floor.

Lower secondary school students’ conversations when constructing representations of plate boundaries in Minecraft

This article explores students’ (13 y) conversations when they construct representations of plate boundaries in the digital gaming environment of Minecraft Education. It also investigates how the spatial affordance of Minecraft Education may support students’ meaning making of the theory of plate tectonics. Video data of classroom activities with Minecraft Education was collected with head cameras. Analysis of the data shows that students overall used everyday language during the activity of constructing plate boundaries. The study suggests i.) a need for more teacher-centeredness and structured scaffolding of scientific language during small group work in Minecraft, and ii.) a need to increase the spatial extent in Minecraft Education to include both the subsurface and more regional scales of resolution.

Site-Specific Ground-Motion Waveform Generation Using a Conditional Generative Adversarial Network and Generalized Inversion Technique

ABSTRACT Accurate ground-motion simulations are essential for seismic hazard assessments and engineering practices. Herein, we propose a novel method combining conditional generative adversarial networks (cGANs) and the generalized inversion technique (GIT) to generate site-specific and variability-controlled strong-motion seismograms. The cGANs calculate synthetic seismogram without amplitude scales. The GIT is to separate the source, path, and site characteristics from the Fourier amplitude spectrum (FAS) of the observed seismograms. This method is applied to plate boundary earthquakes off the Pacific coast of Tohoku, Japan. It successfully generates a set of strong-motion seismograms at a given magnitude, distance, and observation station. The output waveforms reproduce the P and S waves as well as coda waves. We validate the method through a quantitative comparison with observed seismograms in terms of both time-domain duration and frequency-domain amplitude characteristics, using metrics of peak ground acceleration (PGA), peak ground velocity, FASs, response spectra, and waveform duration time. The validation results show that the variation in the PGA of the observed seismograms and the synthetic seismograms has a standard deviation of 0.643, and the duration of the seismograms has a standard deviation of 0.346, comparable to the standard deviations seen in the previous studies. Our approach offers high accuracy in stochastic finite-source modeling for a period of 1 s or shorter. The two features of the method, site-specificity and variability control, can contribute to further improvements in seismic hazard assessment by incorporating empirical information based on observed seismograms.

Impact of Backing Plate and Thermal Boundary Conditions for High-Speed Friction Stir Welding of 25-mm Thick Aluminum Alloy 7175-T79

Impact of Backing Plate and Thermal Boundary Conditions for High-Speed Friction Stir Welding of 25-mm Thick Aluminum Alloy 7175-T79

Along‐strike variation of fault slip rate of a transform plate boundary in Tierra del Fuego (South Patagonia)

AbstractThe Magallanes–Fagnano Fault is an active strike‐slip structure accommodating the relative displacement between South America and the Scotia Plate. The Chilean portion of the fault is poorly studied because most of it runs below the sea level in the Strait of Magellan. Our tectonic geomorphological study is focused on a rare onshore fault section, along which streams horizontally deflected by hundreds of metres since the last main deglaciation are compatible with a dominant left‐lateral fault kinematics and yield a slip rate of 15.7 ± 2.4 mm/year. This rate is between 2 and 3 times higher than the one estimated on the Argentinian portion of the fault over the same period. This spatial variation may be due to both glacial unloading on the fault zone and/or structural factors. These results point out the need to study strike‐slip faults on several portions to unravel behaviour changes related to internal or external forcing.

Rift obliquity in the Northern Volcanic Zone in Iceland using UAV-based structural data

Rift obliquity is known to be very frequent at divergent plate boundaries and shear components in the opening of fissures have been detected through seismicity in the monitoring of recent rifting events. However, the field structural evidence of obliquity in active rifts for Holocene events is still poorly studied, despite its relevance at both the lithospheric and upper-crust scales.We performed extensive high resolution UAV surveys on Holocene terrains in four areas along the active rift in the Northern Volcanic Zone of Iceland, to map structures and morphotectonic features in detail. We did a kinematic analysis and found an increase of the obliquity from N to S.We further discuss the origins of this obliquity and suggest that the combination of inherited structures reactivation, magmatic intrusions, and topography influence on the propagation of both fractures and magma intrusions can have an impact on the mismatch between rift opening direction and the far-field strain.We conclude that obliquity can be inferred for Holocene rifting events from surface structural data. It is thus fundamental to take into account the existing structural pattern and its (re)activation history in active rift zones to avoid missing out on possible preferential structure reactivation, magma propagation ways, and eruption sites during the monitoring of ongoing rifting events.

Robust sulfur and oxygen isotope evidence for a highly anoxic paleoenvironment in Late Devonian seawater: Insights from marine anhydrites in the Zaige Formation, South China

A global-scale mass extinction event occurred in the Late Devonian, and its triggering mechanisms have remained a subject of controversy. The key to resolving this controversy lies in elucidating paleoenvironmental characteristics and their coupling relationship with the mass extinction. Marine evaporites serve as sensitive indicators that can be used for reconstructing the geochemical history of ancient oceans. In this study, we present a dual sulfur and oxygen isotopic dataset of marine anhydrites in the Lunatian gypsum deposit (Yunnan, China), and investigate its paleoenvironmental significance within the framework of the Late Devonian mass extinction. These anhydrites originate in arid and hot continental margin seas near convergent plate boundaries, forming as a result of the dehydration of marine gypsums during burial processes. Anhydrites exhibit a conspicuous positive anomaly in δ34S values (23.6–25.1‰, averaging 24.3‰), primarily attributed to bacteria-mediated sulfate reduction. The δ18O values of the Lunatian anhydrites (12.2–15.2‰, averaging 14.0‰) indicate a slightly higher oxygen isotope composition for Late Devonian seawater (8.7–11.7‰, averaging 10.5‰) compared to those of contemporary seawater (9.5–10.1‰). Such an oxygen isotopic characteristic for the Lunatian anhydrites is intricately associated with the chemical processes that occur during the diffusion of sulfides into the upper seawater. The anhydrite-forming evaporating system has a noticeable brine layering where the surface seawater is less dense than the deeper brine, and bacterial processes gradually increase δ34S in the deep brine. When the deep brine and surface seawater rapidly mix, it leads to a significant positive shift in the sulfur isotope composition of the resulting evaporites. Our sulfur and oxygen isotope data from Lunatian anhydrites validate the existence of an exceedingly anoxic marine environment during the Late Devonian, potentially serving as a pivotal factor in the mass extinction event.

Double-Vergent Plate Boundary Faults and Triggered Coseismic Rupture of the 2022 Chihshang Doublet Earthquake Occurred in Eastern Taiwan

Abstract The earthquake doublet of Mw 6.5 and 6.9 occurred along the west-dipping Central Range fault (CRF) adjacent to the east-dipping creep segment of the Longitudinal Valley fault in eastern Taiwan on 17–18 September 2022. The faulting model derived from the Interferometric Synthetic Aperture Radar and Global Positioning System observation suggests that the west-dipping CRF rupture is responsible for the two mainshocks. Meanwhile, the two major earthquakes resulted in ∼100 km of fault slip along the double-convergent plate boundary faults near the Longitudinal Valley in eastern Taiwan. The 2022 Chihshang earthquake sequence filled the seismic gap of the CRF located between the 2006 ML 6.1 Taitung earthquake and the 2013 ML 6.1 Ruisui earthquake. Finally, the 2022 Chihshang earthquake sequence increased the Coulomb failure stress in the southernmost segment of the CRF, which may cause a higher seismic potential in the future.

Unfolding rotational tectonics and topographic evolution from localized verses diffuse plate boundary counterparts

We present a kinematic model developed from geodetic observations, topography analysis and analogue tectonic modelling results, which reveals a striking similarity between the rotational tectonic settings of the Gakkel Ridge-Chersky Range system in the Arctic, and the Central Indian Tectonic Zone within the Indian subcontinent. A crucial aspect of large-scale extensional rift systems is the gradual variation of extension along the rift axis, due to plate rotation about a Euler pole, which may lead to contraction on the opposite side of the Euler pole to form a rotational tectonic system. Our geodetic and topographic analysis, combined with the reanalysis of analogue tectonic modelling results demonstrates such rotational tectonic plate motion in both the Arctic and Indian case. However, the plate boundary between the North American and Eurasian Plates as represented by the Arctic Gakkel Ridge-Chersky Range system is strongly localized, whereas the Central Indian Tectonic Zone that separates the North and South India Plates involves diffuse deformation instead. Furthermore, in both the Arctic and Central Indian we find that the relative Euler rotation pole is located near an indenter-like feature, which possibly controls the present-day rotational tectonics and contrasting topography on opposite sides of the Euler pole.

Dynamical Modeling of Fault Slip Rates at the New Zealand Plate Boundary Indicates Fault Weakness

AbstractWe construct a thin‐sheet dynamical model of the New Zealand plate boundary that includes faults. Our model fits fault slip rates, style of distributed deformation, and is constrained by relative plate boundary motion. We assume a pseudo‐plastic rheology and achieve a best fit to slip rate observations with a deviatoric stress magnitude of 20 MPa. Modeled local forces are significant at Puysegur and Hikurangi subduction zones, and smaller forces are related to mantle downwelling beneath South Island and Havre Trough mantle upwelling. Modeled tractions on faults are mostly 5–20 MPa, similar to or slightly smaller than stress magnitudes adjacent to faults. Modeled shear tractions are generally 2–10 MPa, comparable to stress drops during earthquakes. Modeled stress orientations and fault dips suggest that many faults are not optimally oriented for their style of faulting. Notably small traction magnitudes of <5 MPa and shear tractions of <0.5 MPa are modeled for faults in the central North Island Dextral Fault Belt (NIDFB), which we infer to be very poorly oriented. Friction coefficients on faults (ratio of shear stress to effective normal stress) are in the range 0.1–0.3 for major crustal faults such as the Alpine Fault and Marlborough faults, but subduction zones and the NIDFB have values <0.1. We propose that low values of long‐term fault strength, shear stress resolved onto the fault, and overall magnitudes of deviatoric stress in the crust are a consequence of dynamic weakening of faults during fault slip.

Continued Convergence After the Occurrence of a Slab Break‐Off: The Case of the Cyprian Arc

AbstractThe detachment (i.e., break‐off) of down‐going subducting oceanic slabs is a major geodynamic event with far‐reaching consequences, one of which is the reduction of the slab pull force acting on the trailing plate. We investigate the motion of the Sinai Microplate where a recent (∼1 Myr ago) slab break‐off occurred along its sole converging plate boundary (Cyprian Arc) with the overriding Anatolia Microplate. Based on new bathymetric mapping, high‐resolution seismic reflection imaging, geodetic and earthquake data, we show that Sinai is actively moving in a northwest direction with respect to Nubia. Our results indicate that despite the recent slab break‐off, Sinai has and is still being pulled (or pushed) toward the overriding Anatolia Microplate. The continued convergence possibly occurs because of a persistent slab pull force, a suction force induced by the down‐going detached slab and/or by the upper mantle flow induced by the Afar Plume.

Effect of surface radiation on the stability of hypersonic flat plate boundary layer

Transition of hypersonic boundary layer flow is an important problem in aerodynamic research. To obtain the natural transition position using the flow stability analysis method, it is necessary to analysis the neutral curve and the perturbation growth rule. Accurate calculation of the basic flow is the basis and prerequisite for stability analysis. This paper proposes that surface radiation should be considered in the analysis of hypersonic boundary layer flow stability due to the high temperatures at the vehicle surface caused by aerodynamic heating. In order to obtain the influence of surface radiation on the flow stability, we assume that the surface satisfies the energy balance and the gas satisfies calorimetrically perfect gas properties. To simplify the analysis, gas radiation has been neglected. The boundary layer analysis of a hypersonic flat plate with and without radiation has been performed by a combination of flow stability theory and direct numerical simulation. It is found that high temperature radiation on the surface causes energy to transfer outwards. The results showed that surface temperature decreases along the streamwise direction and is no longer an adiabatic surface temperature. Instability interval decreases and the maximum growth rate increases. The range of the instability region expressed in dimensionless frequency values moves to lower frequency as the Mach number increases. Maximum growth rate shows a tendency to increase and then decrease as the Mach number increases.

Thermal Modeling of the Sanbagawa and Ryoke Belts

The Sanbagawa and Ryoke belts were formed in a convergent plate boundary along the eastern margin of Eurasia. Thermal modeling using the geological records of these belts as constraints allows quantitative estimates of both shear heating along the Wadati-Benioff zone and magma fluxes beneath the volcanic arc. In contrast to real-time observations of crustal movement and heat flow, rocks record changes in pressures and temperatures that occur over periods of several million years and can be used to examine conditions from the surface to the mantle. Thermal modeling combined with such geological records helps to bridge the gap in our knowledge between real-time observations of ongoing geological processes and the development of orogenies in convergent plate margins over geological time.

Path and Slip Dependent Behavior of Shallow Subduction Shear Zones During Fluid Overpressure

AbstractElevated pore fluid pressure is proposed to contribute to slow earthquakes along shallow subduction plate boundaries. However, the processes that create high fluid pressure, disequilibrium compaction and dehydration reactions, lead to different effective stress paths in fault rocks. These paths are predicted by granular mechanics frameworks to lead to different strengths and deformation modes, yet granular mechanics do not predict their effects on fault stability. To evaluate the role of fluid overpressure on shallow megathrust strength and slip behavior, we conducted triaxial shear experiments on chlorite and celadonite rich gouge layers. Experiments were conducted at constant temperature (130 and 100°C), confining pressure (130 and 140 MPa), and pore fluid pressures (between 10 and 120 MPa). Fluid overpressure due to disequilibrium compaction was simulated by increasing confining and pore fluid pressure synchronously without exceeding the target effective pressure, whereas overpressure due to dehydration reactions was simulated by first loading the sample to a target isotropic effective pressure and then increasing pore fluid pressure to reduce the effective pressure. We find that the effects of fluid pressure and stress path on the mechanical behavior of the chlorite and celadonite gouges can generally be described using the critical state soil mechanics (CSSM) framework. However, path effects are more pronounced and persist to greater displacements in chlorite because its microstructure is more influenced by stress path. Due to its effects on microstructure, the stress path also imparts greater control on the rate‐dependence of chlorite strength, which is not predicted by CSSM.