Journal of Geophysical Research: Solid Earth最新文献

Recent observations show that certain rupture phase can propagate backward relative to the earlier one during a single earthquake event. Such back-propagating rupture (BPR) was not well considered by the conventional earthquake source studies and remains a mystery to the seismological community. Here we present a comprehensive analysis of BPR, by combining theoretical considerations, numerical simulations, and observational evidences. First, we argue that BPR in terms of back-propagating stress wave is an intrinsic feature during dynamic ruptures; however, its signature can be easily masked by the destructive interference behind the primary rupture front. Then, we propose an idea that perturbation to an otherwise smooth rupture process may make some phases of BPR observable. We test and verify this idea by numerically simulating rupture propagation under a variety of perturbations, including a sudden change of stress, bulk or interfacial property and fault geometry along rupture propagation path. We further cross-validate the numerical results by available observations from laboratory and natural earthquakes, and confirm that rupture “reflection” at free surface, rupture coalescence and breakage of prominent asperity are very efficient for exciting observable BPR. Based on the simulated and observed results, we classify BPR into two general types: interface wave and high-order re-rupture, depending on the stress recovery and drop before and after the arrival of BPR, respectively. Our work clarifies the nature and excitation of BPR, and can help improve the understanding of earthquake physics, the inference of fault property distribution and evolution, and the assessment of earthquake hazard.

Observations 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 $��D^{�″�}�$ anisotropy. Our results show strong seismic anisotropy in approximately half of our study region, indicating that $��D^{�″�}�$ 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 elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

We present observations showing that during episodes of volcanic tremors, the phase of inter-station cross-correlations becomes stable. We propose a new quantity, the phase coherence, to identify the differential phase stability in recordings obtained from a single pair of stations, which is extrapolated to the seismic network. Then, we present a new approach based on the estimation of differential travel times through the differential phase measurements, to locate the sources of tremors occurring at the end of 2015 at the Klyuchevskoy Volcanic Group in Kamchatka, Russia. We present evidence supporting the existence of two types of activity happening simultaneously during the tremor episode: the main tremor source, originating from a region located between 7 and 9 km depth under the main volcanoes, and the widespread occurrence of weak low-frequency earthquakes occurring at random locations. We show how the phase coherence and the differential phases can be used to provide information on the stability of the tremor source position and to estimate its location.

Settling-driven gravitational instabilities (SDGIs) can form at the base of buoyant particle-laden suspensions, modulating particle sedimentation in various settings such as meteorological and volcanic clouds, fluvial plumes, magma chambers, submarine hydrothermal plumes, or industrial emissions. These instabilities result in the formation of rapidly descending currents called ‘fingers’ within which fine particles settle faster collectively than individually. This study investigates SDGI triggering conditions underneath volcanic ash clouds through analogue experiments considering sedimentation from aqueous particle suspensions. We confirm that the conditions for which SDGIs develop are controlled by two dimensionless numbers: B_f (ratio of the characteristic finger velocity to the individual particle settling velocity); and B_i (ratio of timescale for individual particle settling to that for collective settling controlled by inertial drag). SDGIs are triggered for values of B_f and B_i > 1 for which particles are fully coupled with the flow within fingers. Using these parameters, we produce a regime diagram for the 2010 eruption of Eyjafjallajökull (Iceland) that describes particle settling as a function of particle concentration and size. More studies are needed to produce a general regime diagram accounting for the evolution of SDGIs properties with eruption and atmospheric parameters. Nonetheless, our study confirms that fingers affect sedimentation from volcanic clouds with high ash volume fractions above 10⁻⁶ vol.%. Our validation of criteria predicting the onset of fingers due to SDGIs constitutes a step forward toward the incorporation of these collective settling processes in volcanic ash transport and dispersion models.

Earthquake interaction across multiple time scales can reveal complex stress evolution and rupture patterns. Here, we investigate the role of static stress change in the 2023 Mw 7.8 and 7.6 earthquake doublet along the East Anatolian Fault (EAF), using simulations of 19 historical earthquakes (M ≥ 6.1) and the 2023 earthquake doublet from 1822 to 2023. Focusing on six cascading sub-events during the 2023 Kahramanmaraş earthquake doublet, we reveal how one sub-event's stress alteration can impact the emergence and rupture of subsequent sub-events. Our analysis unveils that the 2023 Mw 7.8 earthquake was delayed due to stress shadow effects from historical events, while the 2023 Mw 7.6 earthquake was accelerated as a result of stress increases from historical events and ultimately triggered by the 2023 Mw 7.8 earthquake. This study underscores the importance of grasping earthquake preparation, rupture initiation, propagation, and termination in the context of intricate fault systems worldwide. Based on these results, we draw attention to increased seismic hazards in the Elazig-Bingol seismic gap of the EAF and the northern section of the Dead Sea Fault (DSF), necessitating increased monitoring and preparedness efforts.

Interseismic coupling maps and, especially, estimates of the location of the fully coupled (locked) zone relative to the trench, coastline, and slow slip events are crucial for determining megathrust earthquake hazard at subduction zones. We present an interseismic coupling inversion that estimates the locations of the upper and lower boundaries of the locked zone, the lower boundary of the deep transition zone, and downdip gradient of creep rate in the transition from locked to freely creeping in the downdip transition zone. We show that the locked zone at Cascadia is west of the coastline and 10 km updip of the slow slip zone along much of the margin, widest (25–125 km, extending to ∼19 km depth) in northern Cascadia, narrowest (0–70 km) in central Cascadia, with moment accumulation rate equivalent to a M_w 8.71 and M_w 8.85 earthquake for 300- and 500-year earthquake cycles. We find a steep gradient in creep immediately below the locked zone, indicative of propagating creep, along the entire margin. At Nankai, we find three distinct zones of locking (offshore Shikoku, offshore southeast Kii peninsula, and offshore Shima peninsula) with a total moment accumulation rate equivalent to a M_w 8.70 earthquake for a 150-year earthquake cycle. The bottom of the locked zone is nearly under the coastline for all three locked regions at Nankai and is positioned 0–5 km updip of the slow slip zone. In contrast with Cascadia, creep rate gradients below the locked zone at Nankai are generally gradual, consistent with stationary locking.