Occurrence Of Large Earthquakes Research Articles

Clastic injectites and seismic-induced liquefaction in latest Quaternary travertine deposits (Serre di Rapolano, Italy)

ABSTRACT The liquefaction and consequent injection of clastic dykes (injectites) is often associated with the occurrence of large earthquakes. Injectites are most documented in thick unconsolidated sedimentary sequences and less commonly in rocks as a result of localized clastic mud overpressure due to earthquake-induced shaking. This paper documents near-surface clastic dykes injected into Middle-Late Pleistocene travertine deposits following repeated seismic events. The novelty is that there is a clear relationship between the development of the injectites and the activity of seismogenic faults that cut the entire body of the travertine, deposited at least at 45.5 ± 8.5 ka, and the topographic surface. The clastic slurry, derived from the fluidization of silty-sandy continental deposits at the base of, or intercalated with the travertine body, was in geometric connection with the faults. When this occurred, the over-pressured slurry was injected into the fault zone, filling the fractures and assisting their development by over-pressure hydrofracturing and counteracting the frictional forces on the fault. This evidence, in addition to documenting another case of liquefaction and associated emplacement of clastic dykes, has a twofold significance: i) it records the effects of clastic dykes intruded at very shallow depths, up to the surface, within a capable fault zone; ii) it documents a polyphase injection of clastic material into travertine deposits, never before described, induced by repeated liquefaction phenomena following major seismic events, thus highlighting and reinforcing the importance of travertine for the study of active tectonics and earthquake studies.

Stick-to-slip transition characterized by nucleation and emission of dislocations and the implications in earthquake nucleation

Stick-slip friction exists widely in our life especially the occurrence of large earthquakes, but people cannot predict and control by a limited understanding of the mechanisms involved. In the present work, the whole process of stick-to-slip transition has been investigated through digital image correlation and acoustic emission. Two phases, namely, the nucleation and abrupt rupture phases have been discovered during the transition that are characterized well by nucleation and transient emission of dislocations, which may support the combination of pre-slip and cascade-up models. Based on the findings simple yet analytical expressions then have been obtained to predict the earthquake cycles consistent with available simulations and practical observations.

Remagnetization of Pre‐Fan Sediments Offshore Sumatra: Alteration Associated With Seismogenic Diagenetic Strengthening

AbstractIncreases in temperature and pressure caused by rapid burial of sediments seaward of the Sumatra subduction zone have been hypothesized to trigger dehydration reactions that diagenetically strengthen sediments and contribute to the formation of an over‐pressured pre‐décollement, which together facilitate the occurrence of large shallow earthquakes. We present paleomagnetic, rock magnetic, and electron microscopic analyses from drill cores collected offshore Sumatra at Site U1480 during IODP Expedition 362 that support this hypothesis. The older pre‐fan units (Late Cretaceous to early Paleocene) were deposited when Site U1480 was moving rapidly northward with the Indian plate from a paleolatitude of 50° to 30°S, which would equate to expected absolute paleomagnetic inclinations of 70°–43°. Most of the older pre‐fan sediments, however, have shallow observed inclinations (shallower than ±20°), indicating that the sediments were overprinted when Site U1480 was located near the paleoequator, as it has been since the early Oligocene. Electron microscopic observations reveal that the pre‐existing detrital magnetite grains have undergone pervasive dissolution and alteration by hydrothermal fluids. The diagenesis observed is consistent with mineral dehydration, possibly driven by rapid burial of pelagic sediments by the ∼1250 m thick Nicobar Fan sequence. In addition, the elevated burial temperature also facilitated the smectite to illite conversion reaction. We hypothesize that chemical reactions resulted in the formation of fine‐grained magnetite that records a chemical remanent magnetization overprint. This overprint is consistent with the alteration occurring after burial by the thick Nicobar Fan sequence sometime in the past few million years.

Unraveling the roles of fault asperities over earthquake cycles

Fault asperities can produce concentrated slips during large earthquakes and intensify damage. However, how asperities control fault behavior during other phases of the earthquake cycle remains poorly known. Here, we conduct friction experiment on a laboratory fault featuring two prominent geometric asperities to directly image their influences on long-term and short-term seismicity, and the nucleation and propagation of the mainshock. The laboratory observations and supporting numerical simulations reveal that one asperity located in the fault interior behaves firstly as a mechanical attractor to long-term seismicity, then a barrier to preseismic slow slip, and ultimately a source of large coseismic slip. In contrast, another asperity located near the fault margin primarily undergoes persistent aseismic slip throughout most phases of the earthquake cycle. These results provide new insights into how asperities partition strain across a broad spatiotemporal domain, establishing a physical link between long-term and short-term fault behaviors and the occurrence of large earthquakes.

Dynamic Triggering of Earthquakes in Northeast Japan before and after the 2011 M 9.0 Tohoku-Oki Earthquake

ABSTRACT Previous studies documented a relative scarcity of remote dynamic triggering of earthquakes in Japan and suggested that it could be related to Japan’s predominantly compressive tectonic regime or the more frequent occurrence of large earthquakes in Japan. For example, remote triggering in California, characterized by extensional tectonics, occurs at levels of stress change as small as 0.1 kPa, whereas in Japan, transient stresses ≥30 kPa are required. However, the dynamic triggering threshold in Japan, following the 2016 Mw 7.0 Kumamoto earthquake, has been found to be of just a few kilopascals, significantly smaller than reported previously. It was proposed that a decrease in the triggering threshold may have taken place in Japan, in particular at volcanic and geothermal areas, after the 2011 Mw 9.0 Tohoku-Oki earthquake. In this study, we analyze the possible change in dynamic triggering conditions in five areas in northeast Japan, where swarm earthquakes have occurred immediately after the Tohoku-Oki earthquake. The triggering thresholds in these five areas have been estimated based on the analysis of waveform recordings of 49 teleseismic earthquakes that occurred between 2004 and 2020. A decrease of the triggering threshold (or triggering ability) is apparent in all but one region. However, a statistical significance Kolmogorov–Smirnov test could not reject, at a 5% level, the null hypothesis stating that “the distribution of dynamic stress changes for triggering earthquakes that occurred before and after the Tohoku-Oki event is the same.” We interpret the changes in the triggering threshold to be related to the pore pressure increase (and thus a fault strength decrease) in the crust following the Tohoku-Oki earthquake. Our results also indicate that dynamic triggering in Japan is more common than reported previously.

Machine learning technique in the north zagros earthquake prediction

Studying the changes in seismicity, and the potential of the occurrences of large earthquakes in a seismic zone is not only extremely important from the aspect of seismological research, but it is additionally significant in the decisions of crisis management. Since, nowadays Machine learning techniques have proven the high ability for analyzing information, and discovering the relations among the parameters, in this research were tested some of these techniques for the earthquake prediction. For analysis, the north Zagros seismic catalogue was selected. A region that is an active seismic zone, and large cities are located there. Moreover, nine seismic parameters were used to study the possibility of large earthquake prediction for 1 month using three different Machine Learning (ML) techniques (Artificial Neural Network (ANN), Random Forest, and Support Vector Machine (SVM)). The accuracy of prediction models was evaluated using four different statistical measures (recall, accuracy, precision, and F1-score). The results showed that the (ANN) method is more accurate than other methods. Based on three investigated methodologies, greater accuracy results have been produced to forecast the earthquakes with bigger scale earthquakes about the completeness of the seismic catalogue in large magnitude. These achievements promise the possibility of successful prediction in a short period, which is hopeful for better crisis management performance.

Stress disturbance around Xianshuihe fault zone in the eastern Qinghai–Tibet Plateau and implication for fault stability

The Xianshuihe fault zone in the eastern Qinghai–Tibetan Plateau is an important active tectonic boundary. Understanding its stress state is important for characterizing the dynamic evolution of the Qinghai–Tibet Plateau and the mechanism of the frequent occurrence of large earthquakes. Using 30 years of in-situ stress data from the Xianshuihe active fault zone, we statistically analyzed the spatial distribution characteristics of the stress in the region. The study area is generally characterized by a strike-slip stress field. Nevertheless, the stress state is vulnerable to topography and shows high spatial variation near the Earth’s surface at a depth of 0–400 m. The local stress near the fault zone varies from the far-field stress. The orientations of the maximum horizontal principal stress possess an elliptical shape around the fault zone, while its magnitudes become hump-like as the distance increases from the fault. The large difference in properties between the fault zone and its adjacent rocks contributes to the differentiation of the direction of the local stress field near the fault. The results allow us to formulate a preliminary hypothesis that a rigid lateral extrusion model may control the nonuniformity of the local stress field in the Xianshuihe fault zone and preferentially interpret the tectonic uplift of the southeastern margin of the Qinghai–Tibet Plateau. Further, the stress accumulation in the shallow crustal regions of the Xianshuihe fault zone is relatively high, indicating that some segments of the fault zone are critically unstable. Kangding area (the Zheduotang segment and the Yalahe segment) and Luhuo segment hold relatively high potential for large earthquakes. The results of this study are of great significance for revealing the mechanism of fault–stress field interactions and for understanding the dynamic evolution mechanism of the uplift of the Qinghai–Tibet Plateau.

Seasonal Modulation of Crustal Seismicity in Northeastern Japan Driven by Snow Load

AbstractNumerous studies have reported that surface hydrological loading can seasonally modulate seismicity rates at crustal depths. For example, substantial winter snow accumulation occurs across the Japanese Islands, and these snowy regions appear to have seasonally modulated the occurrence of previous large inland earthquakes. Therefore, it is important to investigate the impact of seasonal stress changes on crustal seismicity to deepen our understanding of earthquake generation. Here we constrain seasonal changes in the surface load across northeastern Japan using Global Navigation Satellite System surface displacements and evaluate the potential relationship between temporal trends in inland seismicity and estimated seasonal stress changes. The spatial distribution of the seasonal surface load is consistent with snow depth along the Sea of Japan. The inland seismicity beneath northeastern Japan is modestly modulated by the seasonal stress changes that are induced by the annual snow load. However, this seasonal response is weaker than that in other regions. This weak modulation may be due to the small surface‐load‐induced stress perturbation relative to the long‐term‐averaged stressing rate and/or the limited presence of crustal fluids to trigger seismicity in Japan.

Unraveling Earthquake Clusters Composing the 2014 Alto Tiberina Earthquake Swarm via Unsupervised Learning

AbstractEarthquake swarms represent a particular mode of seismicity, not directly related to the occurrence of large earthquakes (e.g., aftershocks) but rather driven by external forcing such as aseismic deformation or fluid migration in fault systems. Sometimes their occurrence overlaps with observable geodetic signals in space and time, indicating a direct link. However, the low resolution of geodetic observations tends to obscure the small scale spatial and temporal dynamics of swarms. In this work, we automatically extract clusters of seismicity related to the 2014 Alto Tiberina swarm sequence (Italy) using an unsupervised clustering approach that exploits space and time information of the seismicity. The quantitative characterization of each cluster indicates that the overall swarm is composed of spatially and temporally confined (sub) swarms each of which could potentially be driven by small‐scale aseismic deformation process. This observation aligns with similar findings during slow slip events in subduction zones.

The Occurrence of Large Holocene Earthquakes in the Kultor Fault Zone, Northern Issyk-Kul Area, Tien Shan Based on Radioisotope Dating

The Occurrence of Large Holocene Earthquakes in the Kultor Fault Zone, Northern Issyk-Kul Area, Tien Shan Based on Radioisotope Dating

Unmanned Aerial Vehicle Aeromagnetic Analyses Reveal the Causative Faults and Seismogenic Mechanism of the 2021 Ms 6.0 Luxian, China, Earthquake

Abstract On 16 September 2021, an Ms 6.0 earthquake occurred in Luxian, Sichuan, China, breaking the historical record of no earthquake with magnitude ≥ M 6 along the Huaying Mountain fault belt. The regional geological structure is primarily controlled by the northeast-striking fault belt, but the long axis of the isoseismic line, distribution of early aftershocks and coseismic rupture plane all strike northwest, posing challenges to the seismogenic mechanism. To investigate this, we conducted a 400 km2 unmanned aerial vehicle (UAV) aeromagnetic survey near the epicenter, with a line spacing of 1 km. Through aeromagnetic analyses, combined with the spatial distribution of relocated foreshocks and aftershocks, we outline potential basement causative faults and a change in the structural trend between the shallow and deep portions of the seismic zone. We conjecture that the Luxian earthquake was triggered by the hydrofracturing-driven reactivation of a pre-existing northwest-striking and southwest-dipping basement fault, of which the upward propagation induced extrusion and dislocation at the hypocenter within the sedimentary layer. The special structural configurations for the focal area could contribute to the stress concentration and occurrence of large earthquakes.

Revisiting Vrancea (Romania) Intermediate-Depth Seismicity: Some Statistical Characteristics and Seismic Quiescence Testing

Background: The intermediate-depth seismicity in the Vrancea region (Romania) is characterized by localized and persistent earthquake activity that culminates about two or three times in a century with the occurrence of a large event (M ≥ 6.5). Here we have revisited some important seismicity characteristics, using earthquake catalog data spanning two different time periods: 1960–1999 and 2005–2013. Methods: we have determined the b-value of the frequency-magnitude distribution of earthquakes, using a maximum likelihood procedure, and estimated the parameter β to quantify anomalous seismicity rate decreases and increases. Results: by using data from the first period, we have confirmed the existence of a decreased b-value in the deepest part of the seismogenic zone; by using data from the second period, we have statistically confirmed the seismic quiescence that preceded the occurrence of the 1977 M7.4 Vrancea earthquake. Conclusions: the decreased b-value has been interpreted either in terms of an increased lithostatic stress with depth or as an indicator of the depth range where the next major Vrancea earthquake may occur. The time variation of the seismicity parameter β may reveal anomalous seismic quiescence and increased earthquake rates that may precede the occurrence of large earthquakes.

Sediment buoyancy controls the effective slab pull force and deviatoric stress along trenches: Insights from 3D free-subduction model

In the complicated force equilibrium of subduction dynamics, the slab pull resulting from the negative buoyancy of the lithosphere is considered the primary driving force of plate tectonics. In contrast, the isostatic restoring force due to lighter materials covering the oceanic crust (e.g., sediment) has been recognized as a relatively minor contributor. However, as the restoring force counteracts the gravitational pull, the resultant interaction has the potential to control the force balance governing the subduction system. Here, we conducted a series of three-dimensional viscoelastic-free subduction simulations with varying sediment distributions, thicknesses, and densities. We found that a higher restoring force suppressed the trench retreat and decreased the trench velocity. The trench curvatures increased with a high differential restoring force along the trench strike over time. Our results consistently showed a negative correlation between the effective slab pull and the isostatic restoring forces. Moreover, the magnitude and distribution of the isostatic restoring force along the trench strike determined the deviatoric stress level and location of the stress concentration within the subducting slab. We also found that locally abundant sediment margins may lead to a high deviatoric stress in the trench, which is associated with the occurrence of large earthquakes. We argue that sediment buoyancy can significantly affect the subduction kinematics and dynamics to a greater degree than previously suggested. Our findings provide fresh insights into the development of a global predicted velocity model, that reflects the effect of the isostatic restoring force, and provides a better understanding of subduction earthquakes.

Kalman Filter, ANN-MLP, LSTM and ACO Methods Showing Anomalous GPS-TEC Variations Concerning Turkey’s Powerful Earthquake (6 February 2023)

On 6 February 2023, at 1:17:34 UTC, a powerful Mw = 7.8 earthquake shook parts of Turkey and Syria. Investigating the behavior of different earthquake precursors around the time and location of this earthquake can facilitate the creation of an earthquake early warning system in the future. Total electron content (TEC) obtained from the measurements of GPS satellites is one of the ionospheric precursors, which in many cases has shown prominent anomalies before the occurrence of strong earthquakes. In this study, five classical and intelligent anomaly detection algorithms, including median, Kalman filter, artificial neural network (ANN)-multilayer perceptron (MLP), long short-term memory (LSTM), and ant colony optimization (ACO), have been used to detect seismo-anomalies in the time series of TEC changes in a period of about 4 months, from 1 November 2022 to 17 February 2023. All these algorithms show outstanding anomalies in the period of 10 days before the earthquake. The median method shows clear TEC anomalies in 1, 2 and, 3 days before the event. Since the behavior of the time series of a TEC parameter is complex and nonlinear, by implementing the Kalman filter method, pre-seismic anomalies were observed in 1, 2, 3, 5, and 10 days prior to the main shock. ANN as an intelligent-method-based machine learning also emphasizes the abnormal behavior of the TEC parameter in 1, 2, 3, 6, and 10 days before the earthquake. As a deep-learning-based predictor, LSTM indicates that the TEC value in the 10 days prior to the event has crossed the defined permissible limits. As an optimization algorithm, the ACO method shows behavior similar to Kalman filter and MLP algorithms by detecting anomalies 3, 7, and 10 days before the earthquake. In a previous paper, the author showed the findings of implementing a fuzzy inference system (FIS), indicating that the magnitude of the mentioned powerful earthquake could be predicted during about 9 to 1 day prior to the event. The results of this study also confirm the findings of another study. Therefore, considering that different lithosphere–atmosphere–ionosphere (LAI) precursors and different predictors show abnormal behavior in the time period before the occurrence of large earthquakes, the necessity of creating an earthquake early warning system based on intelligent monitoring of different precursors in earthquake-prone areas is emphasized.

Large Earthquakes in Subduction Zones around the Polar Regions as a Possible Reason for Rapid Climate Warming in the Arctic and Glacier Collapse in West Antarctica

A correlation is observed between changes in the level of Earth’s seismic activity and increments of the atmospheric methane concentration over the past 40 years. Trigger mechanisms are proposed for methane emissions and glacier collapse in polar regions. These mechanisms are due to deformation waves caused by large earthquakes in subduction zones located near the polar regions: the Aleutian and Kuril–Kamchatka subduction zones, closest to the Arctic, and the Antarctica–Chilean and Tonga–Kermadec–Macquarie subduction zones. Disturbances of the lithosphere are transmitted over the distances of 3000–4000 km and more at a speed of about 100 km/year. Additional associated stresses come to the Arctic and Antarctica several decades after the occurrence of large earthquakes. In the Arctic zone, additional stresses affect the low-permeability structure of gas bearing sedimentary strata, causing increased methane emission and climate warming. In West Antarctica, deformation waves could trigger the acceleration and intensive collapse of West Antarctic glaciers, which has been observed since the 1970s. These waves are also capable of activating dormant volcanoes located under the sheet glaciers of West Antarctica, leading to an increase in heat flux, to the melting of ice at the glaciers’ base, and to their accelerated sliding towards the ocean, as is happening with the Thwaites Glacier.

Python code for "Deciphering Influential Factors for Large Earthquake Occurrence from the Seismic Catalogs using Explainable Machine Learning"

Python code for "Deciphering Influential Factors for Large Earthquake Occurrence from the Seismic Catalogs using Explainable Machine Learning"

Earthquake hazard characterization by using entropy: application to northern Chilean earthquakes

Abstract. The mechanical description of the seismic cycle has an energetic analogy in terms of statistical physics and the second law of thermodynamics. In this context, an earthquake can be considered a phase transition, where continuous reorganization of stresses and forces reflects an evolution from equilibrium to non-equilibrium states, and we can use this analogy to characterize the earthquake hazard of a region. In this study, we used 8 years (2007–2014) of high-quality Integrated Plate Boundary Observatory Chile (IPOC) seismic data for > 100 000 earthquakes in northern Chile to test the theory that Shannon entropy, H, is an indicator of the equilibrium state of a seismically active region. We confirmed increasing H reflects the irreversible transition of a system and is linked to the occurrence of large earthquakes. Using variation in H, we could detect major earthquakes and their foreshocks and aftershocks, including the 2007 Mw 7.8 Tocopilla earthquake, the 2014 Mw 8.1 Iquique earthquake, and the 2010 and 2011 Calama earthquakes (Mw 6.6 and 6.8, respectively). Moreover, we identified possible periodic seismic behaviour between 80 and 160 km depth.

The 2022 Mw 6.7 Menyuan Earthquake on the Northeastern Margin of the Tibetan Plateau, China: Complex Surface Ruptures and Large Slip

ABSTRACT On 8 January 2022, the Mw 6.7 Menyuan earthquake occurred near the stepover of the Lenglongling (LLLF) and Tuolaishan (TLSF) faults of the Qilian–Haiyuan fault zone in the middle of the northeastern Tibetan plateau. Field investigations and unmanned aerial vehicle-based photogrammetry revealed that the earthquake generated five surface rupture zones with different strikes and kinematic properties. Two large rupture zones, R1 (∼22.8 km long) and R2 (∼3.9 km long), occurred along the northern branch of the western LLLF and the eastern segment of TLSF, respectively, and featured left-lateral strike slips. Among the three small rupture zones, the left-lateral strike-slip-type R3 (0.6 km long) was located in the extension direction of R2, whereas the thrust-type R4 (∼3.3 km long) and R5 (∼1.1 km long) zones were located north of the central section of R1. These complex multifault ruptures were caused mainly by the rupture of strike-slip faults on both sides of the stepover structure. A small amount of compressive shortening strain was released during the earthquake due to regional oblique compression. The total length of the rupture zone was ∼31.7 km; the maximum left-lateral and vertical offsets were 3.5 ± 0.3 and 0.47 ± 0.04 m, respectively. Compared with the relationship observed between coseismic slips and magnitudes in historical and modern earthquakes in western China, the 2022 Menyuan earthquake produced a large coseismic slip in relation to its magnitude. The distribution characteristics of the aftershock belts and their relationship with rupture zones showed that the seismogenic fault of the earthquake was nearly east–west-striking TLSF, which may have triggered the rupture of the northern branch of the western LLLF. In addition, only a small segment of TLSF was ruptured, indicating that the accumulated strain could not be released completely during the earthquake and that this remains the most likely area for the occurrence of large earthquakes in the future.

Revealing the hidden signature of fault slip history in the morphology of degrading scarps

Active faults accommodate tectonic plate motion through different slip modes, some stable and aseismic, others characterized by the occurrence of large earthquakes after long periods of inactivity. Although the slip mode estimation is of primary importance to improve seismic hazard assessment, this parameter inferred today from geodetic observations needs to be better constrained over many seismic cycles. From an analytical formulation developed for analyzing fault scarp formation and degradation in loosely consolidated material, we show that the final topographic shape generated by one earthquake rupture or by creep (i.e., continuous slip) deviates by as much as 10–20%, despite a similar cumulated slip and a constant diffusion coefficient. This result opens up the theoretical possibility of inverting, not only the cumulated slip or averaged slip rate, but also the number of earthquakes and their sizes from scarp morphologies. This approach is all the more relevant as the number of rupture events is limited. Estimating the fault slip history beyond a dozen earthquakes becomes very difficult as the effect of erosion on scarp morphology prevails. Our modeling also highlights the importance of trade-offs between fault slip history and diffusive processes. An identical topographic profile can be obtained either with a stable fault creep associated with rapid erosion, or a single earthquake rupture followed by slow erosion. These inferences, derived from the simplest possible diffusion model, are likely to be even more pronounced in nature.

The role of the west-dipping collision boundary fault in the Taiwan 2022 Chihshang earthquake sequence

On 17–18 September 2022, an earthquake sequence with a moment magnitude of 6.6 foreshock and a 7.0 mainshock occurred in southeast Taiwan along the Longitudinal Valley. Several surface breaks and collapsed buildings were observed after the event and one person died. The focal mechanisms of the foreshock and mainshock both had a west-dipping fault plane, which is different from the known active east-dipping boundary fault between the Eurasian Plate and the Philippine Sea Plate. Joint source inversions were performed to better understand the rupture mechanism of this earthquake sequence. The results show that the ruptures mainly occurred on a west-dipping fault. In the mainshock, the slip originated from the hypocenter and propagated toward the north with a rupture velocity of approximately 2.5 km/s. The east-dipping Longitudinal Valley Fault also ruptured, which could be passive and dynamically triggered by the significant rupture on the west-dipping fault. Most importantly, this source rupture model together with the occurrence of large local earthquakes over the past decade strongly supports the existence of the Central Range Fault, which is a west-dipping boundary fault that lies along the north to south ends of the Longitudinal Valley suture.