Fault Rupture Research Articles

Afterslip on Conjugate Faults of the 2020 Mw 6.3 Nima Earthquake in the Central Tibetan Plateau: Evidence from InSAR Measurements

ABSTRACT Afterslip could help to reveal seismogenic fault structure. The 2020 Mw 6.3 Nima earthquake happened in a pull-apart basin within the Qiangtang block, central Tibetan plateau. Previous studies have explained the coseismic and early (<6 mo) postseismic deformation by rupture and afterslip on a normal fault bounding the western side of the basin. Here, we resolved the 19-month Interferometric Synthetic Aperture Radar-measured sequences of postseismic displacements that revealed a second postseismic displacement center ~12 km to the east of the main fault. Fitting the postseismic displacement requires afterslip on both the main fault and an antithetic fault that probably forms a y-shaped pair of conjugate faults in a negative flower structure. Stress-driven afterslip models suggest that the required afterslip on the antithetic fault could be triggered by coseismic rupture of the main fault or by a simultaneous rupture on the antithetic fault. The afterslip on both faults occurred mainly up-dip to the coseismic slip and has released moment ~15%–19% of that by the coseismic rupture. These results provide insights into active extension in the central Tibetan plateau and highlight the complex nature of fault rupture and afterslip.

Coseismic river avulsion on surface rupturing faults: Assessing earthquake-induced flood hazard.

Surface-rupturing earthquakes can produce fault displacements that abruptly alter the established course of rivers. Several notable examples of fault rupture-induced river avulsions (FIRAs) have been documented, yet the factors influencing these phenomena have not been examined in detail. Here, we use a recent case study from New Zealand's 2016 Kaikōura earthquake to model the coseismic avulsion of a major braided river subjected to ~7-m vertical and ~4-m horizontal offset. We demonstrate that the salient characteristics of the avulsion can be reproduced with high accuracy by running a simple two-dimensional hydrodynamic model on synthetic (pre-earthquake) and "real" (post-earthquake) deformed lidar datasets. With adequate hydraulic inputs, deterministic and probabilistic hazard models can be precompiled for fault-river intersections to improve multihazard planning. Flood hazard models that ignore present and potential future fault deformation may underestimate the extent, frequency, and severity of inundation following large earthquakes.

Do soft soil layers reduce the seismic kinematic distress of onshore high-pressure gas pipelines?

AbstractOnshore high - pressure gas pipelines constitute critical infrastructure that usually cross seismic - prone regions and are vulnerable to Permanent Ground Deformations (PGDs) due to active seismic faults. In design, it may not be feasible to avoid fault rupture areas due to various technical, economical and topographic reasons. Moreover, the presence of soil layers affects the PGDs resulting from a tectonic fault, which in turn may alter the seismic demand on the pipeline. The current study investigates numerically the impact of soft soil layers on the seismic kinematic distress of onshore gas pipelines. For this purpose, a decoupled numerical modeling approach is adopted, consisting of two separate finite - element models for the simulation of soil response and pipeline distress, respectively. Soil non - linearities are taken into account utilizing the Mohr - Coulomb constitutive model with isotropic strain softening. An extensive parametric analysis is performed considering different faulting mechanisms and fault dip angles, as well as soil geometry and mechanical properties. Consequently, the maximum absolute values of both tensile and compressive pipeline strains are correlated with the seismic intensity level (i.e., in terms of bedrock offset which is associated with earthquake magnitude via simple relationships). The paper concludes with a set of design charts and tables for the preliminary seismic design of onshore high - pressure gas pipelines. These charts and tables predict with reasonable accuracy pipeline deformations, in terms of strains, for different magnitude, fault type, dip angle, sand type, and varying overlying soil layer thickness.

Impacts of Water and Stress Transfers from Ground Surface on the Shallow Earthquake of 11 November 2019 at Le Teil (France)

The 4.9 Mw earthquake of 11 November 2019 at Le Teil (France) occurred at a very shallow depth (about 1 km), inducing the surface rupture of La Rouvière fault. The question was raised shortly after about the potential impact of a nearby surface quarry. Thanks to satellite differential interferometry, here, we revealed the existence of a secondary surface rupture of the quasi-parallel Bayne Rocherenard fault. A newly processed seismic cross-section allowed us to shape the three-dimensional geometry of the local three-fault system. Assuming that the earthquake was triggered by the impact of meteoric water recharge, our numerical simulations show that the hydraulic pressure gradient at depth was at a maximum during the period of 2010–2019, just before the seismic event. The estimated overpressure at the intersection of the two faults, which is the most probable place of the hypocenter, was close to 1 MPa. This hydraulic effect is about two and a half times larger than the cumulative effect of mechanical stress release due to the mass removal from the surface quarry over the two past centuries. This work suggests a rapid hydraulic triggering mechanism on a network of faults at a shallow depth after a heavy rainfall episode.

Dueling dynamics of low-angle normal fault rupture with splay faulting and off-fault damage

Despite a lack of modern large earthquakes on shallowly dipping normal faults, Holocene Mw > 7 low-angle normal fault (LANF; dip<30°) ruptures are preserved paleoseismically and inferred from historical earthquake and tsunami accounts. Even in well-recorded megathrust earthquakes, the effects of non-linear off-fault plasticity and dynamically reactivated splay faults on shallow deformation and surface displacements, and thus hazard, remain elusive. We develop data-constrained 3D dynamic rupture models of the active Mai’iu LANF that highlight how multiple dynamic shallow deformation mechanisms compete during large LANF earthquakes. We show that shallowly-dipping synthetic splays host more coseismic slip and limit shallow LANF rupture more than steeper antithetic splays. Inelastic hanging-wall yielding localizes into subplanar shear bands indicative of newly initiated splay faults, most prominently above LANFs with thick sedimentary basins. Dynamic splay faulting and sediment failure limit shallow LANF rupture, modulating coseismic subsidence patterns, near-shore slip velocities, and the seismic and tsunami hazards posed by LANF earthquakes.

Dynamic rupture simulations based on interseismic locking models—taking the Suoerkuli section of the Altyn Tagh fault as an example

SUMMARY For simulating the dynamic rupture process in earthquake scenarios, the stress distribution along the fault remains unclear owing to a lack of direct measurements. Regional stress fields are often resolved onto the fault plane to determine the stress distribution along it. To overcome this limitation, we considered different interseismic locking models to better constrain the actual stress distribution. Specifically, we took the Suoerkuli section in the middle of the Altyn Tagh fault, China, and conducted dynamic rupture simulations to obtain possible earthquake scenarios. The surface rupture length and moment magnitude obtained from the simulations were consistent with those of historical earthquakes. Compared with the traditional stress field resolution method, our approach led to better constrained fault rupture extent and distribution characteristics of regional intensity, thereby avoiding overestimations of earthquake damage. We conclude that examining regional seismic hazards and risks based on seismic dynamic rupture simulations that account for the locking ratio of the fault plane is advantageous, and should be encouraged.

Fast and slow intraplate ruptures during the 19 October 2020 magnitude 7.6 Shumagin earthquake

Strong tsunami excitation from slow rupture of shallow subduction zone faults is recognized as a key concern for tsunami hazard assessment. Three months after the 22 July 2020 magnitude 7.8 thrust earthquake struck the plate boundary below the Shumagin Islands, Alaska, a magnitude 7.6 aftershock ruptured with complex intraplate faulting. Despite the smaller size and predominantly strike-slip faulting mechanism inferred from seismic waves for the aftershock, it generated much larger tsunami waves than the mainshock. Here we show through detailed analysis of seismic, geodetic, and tsunami observations of the aftershock that the event implicated unprecedented source complexity, involving weakly tsunamigenic fast rupture of two intraplate faults located below and most likely above the plate boundary, along with induced strongly tsunamigenic slow thrust slip on a third fault near the shelf break likely striking nearly perpendicular to the trench. The thrust slip took over 5 min, giving no clear expression in seismic or geodetic observations while producing the sizeable far-field tsunami.

Aseismic creep and gravitational sliding on the lower eastern flank of Mt. Etna: Insights from the 2002 and 2022 fault rupture events between Santa Venerina and Santa Tecla

Fault creep along the lower eastern flank of Mt. Etna volcano has been documented since the end of the 19th century and significantly contributes to the surface faulting hazard in the area. On 29 October 2002, during a seismic swarm related to dyke intrusions, two earthquakes caused extensive damage and surface faulting in an area between the Santa Venerina and Santa Tecla villages. On the same day after the two earthquakes, an episodic aseismic creep occurred along the Scalo Pennisi Fault close to the Santa Tecla coastline. On 8 February 2022, during another aseismic creep event along the Scalo Pennisi Fault, we observed the reopening of the pre-existing 2002 ground ruptures mostly as pure dilational fractures. We mapped the 2002 and 2022 surface ruptures, and collected data on displacement, length, and pattern of ground breaks. Ground ruptures affected structures located along the activated fault segments, including roads, walls and buildings. The 2002 surface faulting propagation can be ascribed to a sliding of the Mt. Etna eastern flank toward the SE, as also suggested by the related shallow seismicity, and InSAR and geodetic data between 2002 and 2005. For the 2022 event, differential InSAR data, acquired in both descending and ascending views, allowed us to decompose Line of Sight (LOS) displacement into horizontal and vertical components. We detect a ∼ 700 m long and ∼ 500 m wide deformation zone with a downward and eastward motion (max displacement ∼1,5 cm) consistent with a normal fault. We inverted the InSAR–detected surface deformation using a uniform-slip fault model and obtained a shallow detachment for the causative fault, located at ∼300 m depth, within the volcanic pile. This is the first in-depth study along the Scalo Pennisi Fault to suggest a shallow faulting that accommodates Mt. Etna E flank gravitational sliding.

A Three-Dimensional Parametric Study of the Effects of a Fault Rupture on a Multi-Story Building

Earthquake events have shown that besides the earthquake forces, the interaction between the fault rupture and structures could cause a lot of damage. Field observations have revealed the need to design structures for fault-induced loading in regions with active faults. In this study, three-dimensional numerical simulation was carried out to evaluate the performance of a 20-story frame building subjected to the propagation of normal fault rupture. A parametric study was conducted with the propagation of the fault slip (h) with outcrop to the right (s/B=-1.0), at the mid (s/B=0.5), and left (s/B=2.0) of the raft. An advanced hypoplastic sand model (which can capture small-strain stiffness and stress-state dependent dilatancy of sand) was adopted. The Concrete Damaged Plasticity (CDP) model was used to capture the cracking behavior in the concrete beams, columns, and piles. The computed results revealed that the performance of the building due to the propagation of the normal fault rupture depends on the position of the outcrop of the rupture with respect to the building. Among the three simulated cases, the maximum differential settlement occurred when s/B= 0.5. In each case, the lateral movement of the building caused inter-story drift which may induce distress in the structural components of the building. The beams and the columns of the first story were fully damaged in tension when s/B was equal to -1.0 or 0.5.

Influence of Viscous Lubricant on Nucleation and Propagation of Frictional Ruptures

AbstractFluids are pervasive in the Earth's crust and saturate fractures and faults. The combination of fluids and gouge layers developing along faults can generate fluids of different viscosities. Such viscous fluids were found to influence the reactivation, frictional stability of faults, and eventually the dynamics of propagating earthquake ruptures. We reproduced laboratory earthquakes on analog material (PMMA) to study the influence of viscous lubricant on fault frictional stability, rupture nucleation, and propagation under mixed lubrication conditions. Experiments were conducted in dry conditions, and with fluids presenting a viscosity ranging from 1 to 1,000 mPa.s. Through photoelasticity, high‐frequency strain gauge sensors, and accelerometer measurements, we obtained new insights about the influence of lubricant on a characteristic nucleation length, rupture propagation velocity, and local slip and slip rate evolution during the reproduced frictional ruptures. Our experiments show that the presence of a lubricant generating mixed lubricated conditions along the fault induces, (a) a reduction of the frictional resistance, (b) an increase in nucleation length, (c) a decrease in the fracture energy. In addition, the larger the viscosity of the fluids, the larger the reduction of frictional strength and the increase in the nucleation length. Moreover, ruptures occurring under mixed lubricated conditions showed a pulse‐like rather than crack‐like behavior, suggesting that viscous lubrication can induce the transition from crack‐like to pulse‐like rupture along natural faults. We demonstrate, supported by existing theory, that this transition is mainly governed by the stress acting on the fault at the onset of nucleation, which is drastically reduced in presence of a lubricant.

Towards a detailed urban structural mapping using outcrop, geophysical and borehole data integration: Implications for sustainable development

Nowadays, geological conditions and geohazards assessment are essential requirements for sustainable urban planning and development. Hence, with lack of outcrop and near-surface geology, additional geohazards can be imposed. The present study is the first such challenge in one of highly deformed areas in Egypt and, accordingly, an integrated approach for detailed urban geological and structural mapping is given. Therefore, the surface and subsurface field measurements converge on the detailed geological conditions characterization including structural settings in a new urban city, as a case study. In this regard, the measured surface structural data are presented on a detailed geological map considering the surface rock units geometry and distributions. In order to unravel the near-surface geology, geophysical direct current resistivity (DCR) sounding survey is performed to represent subsurface layers distribution considering the surface geological mapping. To reduce the limitations of DCR sounding interpretation, conventional and non-conventional inversion schemes are applied with reference to available borehole data. Evidently, the sum of surface and near-surface mapped fault segments are plotted and analyzed to provide a detailed geometrical structural model. Because the fault rupture hazards and shale layer expansion issues are closely related to the urban sustainable development, a preliminary geohazards assessment is illustrated in the form of fault setbacks calculation and the shale free swelling evaluation. The present approach provides a holistic view of the urban geological mapping and introduces a means of recognizing preliminary potential geohazards problems and/or opportunities at an early stage in any proposed sustainable development, which can support site investigation designing strategies.

The June 2022 Khost earthquake in southeastern Afghanistan: A complicated shallow slip event revealed with InSAR

On June 22, 2022, the Mw6.2 earthquake in southeastern Afghanistan caused a severe disaster. We used the Sentinel-1A ascending and descending track images of the European Space Agency and interferometric synthetic aperture radar (InSAR) to obtain the coseismic surface deformation field of the earthquake, which showed that the earthquake caused complex ruptures of multiple faults and various types. Using the dislocation model of the elastic half-space, we determined the focal parameters and slip distribution on the fault plane of this event. The results reveal that: (1) the seismogenic fault of this event is an unknown fault on the northeastern edge of the Katawaz microblock; (2) The slip on the fault plane is mainly in the range of 0–8 km along the dip, with maximum slips about 2 m at a depth of 2 km, which projected on the surface is 69.44°E, 32.96°N. This event suggests that, similar to the Chaman, Ghazaband and other large faults, the faults inside the microblock also play an important role in adjusting for the collision stress between India and Europe.

The risk of fluid-injection-induced fault reactivation in carbonate reservoirs: an investigation of a geothermal system in the Ruhr region (Germany)

In this study, we investigate the potential for fluid-injection-induced fault reactivation and induced seismicity risk during simultaneous injection-extraction operation in a theoretical geothermal doublet system in a carbonate reservoir in the Ruhr region. Using a coupled three-dimensional thermo-hydro-mechanical approach, we investigate the probability of injection-induced fault rupture. We perform a sensitivity study assuming variability of the fault and matrix permeability, injection/production flow rates, well placement options, rock thermal properties, and evaluate the influence of thermally induced stresses. The ruptured fault areas were calculated based on a Coulomb friction law and a notion that the shear slip is controlled by the ratio of shear to effective normal stresses acting on a pre-existing plane of weakness in the in situ stress field configuration. Ruptured fault areas in the intrinsically not critically-stressed environment, using location-specific empirical correlations, were used to compute local moment magnitudes of potential earthquakes. Based on this study, we conclude that, in the long-term, thermally-induced stresses play a dominant role during fault reactivation and greatly increase the likelihood for induced seismicity. We, therefore, propose that a minimum safe distance between an injection well and a fault should be based primarily on the radius of a thermal plume generated during the expected lifetime of a geothermal system. Results from this study provide valuable insights for the development of future deep geothermal systems in the Ruhr region and other geothermal reservoirs worldwide.

The March 2020, Mw 5.7 Magna, Utah, earthquake—documentation of geologic effects and summary of new research

The March 18, 2020, Mw 5.7 Magna earthquake was the largest earthquake in Utah since the 1992 ML 5.8 St. George earthquake. The Magna earthquake occurred in the northwest corner of the Salt Lake Valley, home to 1.2 million people. Immediately following the earthquake, the Utah Geological Survey organized teams to collect perishable field data on the geologic effects of ground shaking near the epicenter, as well as establish a web-based digital clearinghouse to collect, distribute, and archive data related to the earthquake. This earthquake also coincided with the beginning of the COVID-19 global pandemic, which added extra challenges to our earthquake response. Teams used a small, unmanned aircraft system to obtain aerial photos and videos of geologic effects to supplement ground-based reconnaissance. The observed geologic effects of ground motions from the Magna earthquake include liquefaction in the form of sand boils, tension cracks, lateral spreading, and localized subsidence. No primary surface fault rupture was observed. The areas with the highest observed concentration of liquefaction features were close to the shore of Great Salt Lake and near the epicenter, northeast of the town of Magna. Photos and other documentation of the geologic effects associated with this earthquake are critical in helping to understand the hazards associated with moderate magnitude earthquakes in the Wasatch Front region. The earthquake sequence and associated geologic effects were well documented, due to the proximity to a major metropolitan area and the mainshock and aftershocks occurring within the densest part of the Utah Regional Seismic Network. In the two years since the earthquake, numerous studies have been published documenting and interpreting data to characterize the Magna event and discuss how new data add to what is known about seismic hazards along the Wasatch Front.

Measurements of Wave‐Induced Attenuation in Saturated Metapelite and the Band‐Limitation of Low‐Frequency Earthquakes

AbstractThe most common explanation for the depletion of high frequency waves that defines low‐frequency earthquakes (LFEs) and very low‐frequency earthquakes (VLFEs) is that fault rupture and slip are slower than typical earthquakes. However, it is difficult to rule out the possibility that the high frequency waves are produced during slip, but attenuated near the LFE source. One reason this hypothesis has been poorly tested is that there are no measurements of attenuation on the relevant rocks. We present the results of forced oscillation experiments that measure the frequency‐dependent attenuation of a chlorite‐rich metapelitic schist, a lithology found along the subduction plate boundary where LFEs and VLFEs have been documented. Experiments were run on dry and water‐saturated schist at effective pressures of 2–10 MPa and at frequencies of 2 × 10−5–30 Hz. We find that pore fluids and low effective pressure result in the attenuation of high frequencies. The frequency‐dependent attenuation is consistent with the concomitant operation of two wave‐induced fluid flow mechanisms, squirt flow, and patchy saturation. When the effects of these mechanisms are extrapolated to geologic conditions using rock physics models, our results predict that attenuation is capable of completely diminishing the frequencies depleted in LFEs and VLFEs. Therefore, LFEs and VLFEs may not necessarily record slow fault slip, but possibly the presence of high fluid pressure.

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.

Normal fault rupture propagation in overburdened rock and its impact on the ground deformation profile based on a model experimental approach

Rupture propagation and ground deformation are two major issues encountered during fault propagation. However, few studies have investigated the rupture propagation and ground deformation of a crossing fault field. In this study, three large-scale experimental tests for simulating normal faulting in crossing fault field were conducted to investigate the propagation characteristics and patterns of fault ruptures. In the conducted tests, different factors were considered, such as the fault dip Df, fault width Wf and amount of normal faulting h. The result revealed that the fault rupture pattern was mainly affected by Df, while Wf could affect the initial rupture length and ground tensioned cracks. Three fault rupture patterns can be distinguished under different Df values. In the former two patterns, the main fault rupture trace is in the form of a nonlinear logarithmic spiral line with an outcropping direction of 45°-ψp/2. Although the ground deformation profiles can be fitted roughly by empirical equation, the local deformations at fault fracture zone are not usually predicted well. The Gompertz function can describe the growth of ground settlement in a fault-crossing rock mass field. Therefore, Gompertz was introduced to characterize the growth and geometrical features of ground settlement due to normal faulting. Finally, bell-shaped ground slope curves were identified mathematically according to the first derivation of the Gompertz function. Based on the geometrical parameters, the ground deformation characteristics were furtherly revealed. This research would assist engineers in selecting sites and designing facilities in regions with weak fault fracture zones.

Fault geometry of M6-class normal-faulting earthquakes in the outer trench slope of Japan Trench from ocean bottom seismograph observations

Since the 2011 Mw 9.0 Tohoku-oki earthquake, intra-plate normal-faulting earthquakes, including several M7-class earthquakes, have occurred in the outer trench slope area from the trench to the outer rise along the Japan Trench. Concerns regarding large earthquakes and associated tsunamis have also arisen. Based on aftershock distributions, several outer trench slope normal-faulting earthquakes (hereinafter referred to as outer-rise earthquakes) are likely related to the rupture of multiple faults. However, few observations have clearly shown how multiple faults act during outer-rise earthquakes. During the ocean bottom seismograph (OBS) observations in the outer trench slope area of the central Japan Trench from September 2017 to July 2018, three M6-class normal-faulting earthquakes (Mw 6.2 on September 20, Mw 6.2 on October 06, and Mw 6.0 on November 12) occurred around the OBS network. The near-field OBS observations provided detailed information on hypocenter locations and focal mechanisms of the mainshocks and aftershocks, including immediately after the mainshocks. We investigated the fault configurations of normal-faulting earthquakes based on OBS observations. During the September 2017 earthquake, the mainshock ruptured high-angle normal faults with a dip angle of 65°. Off-fault aftershock activities that were not directly related to the mainshock rupture and could be explained by the stress changes caused by the mainshock were confirmed. However, hypocenter distributions and focal mechanisms of the main and aftershocks of the October and November 2017 earthquakes suggest that the mainshock ruptured multiple faults with various dipping directions, angles, and strike orientations. The complicated fault geometry should be considered a possible fault model for large outer-rise earthquakes and related tsunamis.

Nonsystematic Rupture Directivity of Geothermal Energy Induced Microseismicity in Helsinki, Finland

AbstractThe rupture behavior of microseismicity in fluid‐injection settings with low fault stresses is generally believed to be controlled by the pore pressure, including a tendency of the larger induced earthquakes to rupture into the perturbed volume toward the injection well, implying a degree of predictability. Here, we examine directivity patterns to identify fault planes and rupture directions of the 21 largest earthquakes (local magnitudes, 1.3–1.9) recorded during the 2018 St1 Deep Heat geothermal project near Helsinki, Finland. We use the Empirical Green's Function technique to retrieve per‐seismic‐station corner frequencies, earthquake durations, and directivity trends. After combining the directivity trends with focal mechanisms calculated using principle component analysis, we resolve rupture planes and rupture directions of 10 events. In contrast to studies of induced events at other sites, we find that one event rupture toward, two rupture away, and the remaining rupture parallel to the well. Furthermore, we find that the events prefer mode II failures rather than mode III. These observations provide new constraints for mechanical models of rupture growth in pore‐pressure dominated settings.

Geothermal Anomalies and Coupling with the Ionosphere before the 2020 Jiashi Ms6.4 Earthquake

Rock temperature reflects the adjustment in crustal stress, and the fluctuation of ionospheric electron concentration is closely related to short-term disturbances of the stress field. Their coupling may reveal short-term effects before strong earthquakes. This study explores the rock temperature changes and mechanical-electrical coupling in the lithosphere–ionosphere before the Jiashi Ms6.4 earthquake on 19 January 2020. The observed data were detrended by general polynomial piecewise fitting; three observation points within 150 km of the epicenter were found to show significant temperature fluctuations in the 15 days before the earthquake. The peak occurred synchronously five days before the earthquake, and the variation range was approximately 10−3 orders of magnitude. Five days before the earthquake, the electromagnetic satellite Zhangheng-1 synchronously observed an anomalous electron concentration in orbit near the epicenter, with a maximum value of 2.01 × 1010 m−3. The loading/unloading response ratio (LURR) was calculated using small earthquakes within 100 km of the epicenter; it showed that the large changes in rock temperature and the ionosphere occurred at high LURR, indicating high-stress accumulation in the region. Various anomalies appeared simultaneously and may indicate fault rupture, which may be caused by an acoustic-gravity wave, indicating a synchronous coupling between the lithosphere atmosphere and the ionosphere.