Numerical Modelling of Lithospheric Block-and-Fault Dynamics: What Did We Learn About Large Earthquake Occurrences and Their Frequency?

Dynamics of lithospheric plates resulting in localisation of tectonic stresses and their release in large earthquakes provides important information for seismic hazard assessments. Numerical modelling of the dynamics and earthquake simulations have been changing our view about occurrences of large earthquakes in a system of major regional faults and about the recurrence time of the earthquakes. Here, we overview quantitative models of tectonic stress generation and stress transfer, models of dynamic systems reproducing basic features of seismicity, and fault dynamics models. Then, we review the thirty-year efforts in the modelling of lithospheric block-and-fault dynamics, which allowed us to better understand how the blocks react to the plate motion, how stresses are localised and released in earthquakes, how rheological properties of fault zones exert influence on the earthquake dynamics, where large seismic events occur, and what is the recurrence time of these events. A few key factors influencing the earthquake sequences, clustering, and magnitude are identified including lithospheric plate driving forces, the geometry of fault zones, and their physical properties. We illustrate the effects of the key factors by analysing the block-and-fault dynamics models applied to several earthquake-prone regions, such as Carpathians, Caucasus, Tibet-Himalaya, and the Sunda arc, as well as to the global tectonic plate dynamics.

Mapping Main Structures and Related Mineralization of the Arabian Shield (Saudi Arabia) Using Sharp Edge Detector of Transformed Gravity Data

Saudi Arabia covers most of the Arabian Peninsula and is characterized by tectonic regimes ranging from Precambrian to Recent. Using gravity data to produce the lateral boundaries of subsurface density bodies, and edge detection of potential field data, a new subsurface structural map was created to decipher the structural framework controls on the distribution of gold deposits in Saudi Arabia. Moreover, we detected the relationships between major structures and mineral accumulations, thereby simultaneously solving the problem of edge detectors over complex tectonic patterns for both deeper and shallower origins. Analytic signal (ASg), theta map (TM), TDX, and softsign function (SF) filters were applied to gravity data of Saudi Arabia. The results unveil low connectivity along the Najd fault system (NFS) with depth, except perhaps for the central zones along each segment. The central zones are the location of significant gold mineralization, i.e., Fawarah, Gariat Avala, Hamdah, and Ghadarah. Moreover, major fault zones parallel to the Red Sea extend northward from the south, and their connectivity increases with depth and controls numerous gold mines, i.e., Jadmah, Wadi Bidah, Mamilah, and Wadi Leif. These fault zones intersect the NFS in the Midyan Terrane at the northern part of the AS, and their conjugation is suggested to be favorable for gold mineralization. The SF maps revealed the boundary between the Arabian Shield and Arabian Shelf, which comprises major shear zones, implying that most known mineralization sites are linked to post-accretionary structures and are not limited to the Najd fault system (NFS).

Emergence of Low-Frequency Aftershocks of the 2019 Ridgecrest Earthquake Sequence

ABSTRACT Low-frequency earthquakes (LFEs) generally have relatively stronger spectral components in the lower frequency range compared with what is expected for regular earthquakes based on their magnitude. LFEs generally occur in volcanic systems or deep (>∼15 km) in plate boundary fault zones; however, LFEs have also been observed in nonvolcanic, upper crustal settings. Because there are few studies that explore the spatiotemporal behaviors of LFEs in the shallow crust, it remains unclear whether the shallow-crustal LFEs reflect local attenuation in their immediate vicinity or differences in their source mechanism. Therefore, it is important to identify shallow-crustal LFEs and to characterize their spatiotemporal activity, which may also improve our understanding of LFEs. In this study, we focus on detecting shallow-crustal LFEs and explore the possible generation mechanisms. We analyze 29,646 aftershocks in the 2019 Ridgecrest, California, earthquake sequence, by measuring the frequency index (FI) to identify candidate low-frequency aftershocks (LFAs), while accounting for the magnitude dependency of the FI. Using small earthquakes (ML 1–3) recorded in the borehole stations to minimize the attenuation effects in near-surface layers, we identify 68 clear LFAs in total. Based on their distribution and comparisons with other seismic parameters measured by Trugman (2020), the LFAs possess distinct features from regular events in the same depths range, including low corner frequencies and low stress drops. Events in the close vicinity of LFAs exhibit lower average FI values than regular aftershocks, particularly if the hypocentral distance between an LFA and its neighbors is less than 1 km. Our results suggest that LFAs are related to local heterogeneity or a highly fractured fault zone correlated with an abundance of cross faults induced by the aftershock sequence at shallow depths. Zones of high pore-fluid pressure in intensely fractured fault zones could cause the bandlimited nature of LFAs and LFEs in general.

Tectonic control on the palaeogeographical evolution of the Miocene seaway along the Western Alpine foreland basin

The Miocene of the Western Alpine foreland basin were deposited in a North-South seaway along the active alpine orogenic front. In the subalpine massifs and the southern Jura mountains, the revised Miocene stratigraphy documents a detailed chronology of thrust propagation at the western alpine front, where tectonic activity had a primary influence on seaway palaeogeographical evolution. Here we propose nine palaeogeographical maps during the Miocene, the first of which depicts the initial Miocene transgression at ∼21.0 Ma. Between ∼18.05 Ma and ∼12.0 Ma, a westward retreat of the Miocene Sea occurred in response to the activation of the basal thrust of the Belledonne massif, which in turn triggered successive fault zones from east to west. At ∼10.0 Ma, a major uplift phase intervened and induced a rapid southward retreat of the Miocene Sea. The reconstructed palaeogeographical maps outline the main controls on the foreland basin seaway evolution: (i) the timing of the main thrusts, (ii) the inherited palaeotopography and (iii) eustatic sea-level changes during the Miocene. These reconstructions are integrated at the basin scale, highlighting the southward to westward-directed seaway migration in response to the Belledonne thrust activity that deeply shaped the palaeogeographical evolution during the early to middle Miocene.

Hydrogeochemistry and Precursory Anomalies in Thermal Springs of Fujian (Southeastern China) Associated with Earthquakes in the Taiwan Strait

Analyzing the hydrochemical composition in thermal springs is an advantageous method for studying the coupling mechanism of the deep and shallow fluids in active fault zones. Here we conducted sampling in 30 thermal springs near fault zones in Fujian Province, and the major elements, trace elements, silica, stable isotopes (δD and δ18O) and strontium isotopes were tested in the laboratory. The results show that (1) the thermal springs in the study area can be divided into six types according to the content of the major elements: HCO3-Na, HCO3·SO4-Na, Cl·HCO3-Na, Cl-Na, Cl-Na·Ca and HCO3·SO4-Ca; (2) hydrogen and oxygen isotopes indicate that precipitation is the main source of recharge for thermal springs in the study area, and the recharge height is between 258 m and 1859 m; (3) the content of SiO2 in the thermal spring varies from 18.1 mg/L to 59.3 mg/L. The geothermal reservoir temperature calculated is 90~226 °C, and the circulation depth is 2.9~5.4 km, except for the W10 thermal spring, whose circulation depth is 8.4 km; and (4) the 87Sr/86Sr of the thermal springs in southwestern Fujian and eastern Fujian has obviously different characteristics, indicating the influence of different rock formations on the groundwater cycle process. Additionally, a continuous measurement of the main anions and cations was performed in five thermal springs every three days since January 2020. There were obvious abnormal changes in the hydrochemical compositions, chlorine in four of the five springs, sodium at three springs, and four ions at one spring, which all showed abnormal high-value changes by 15% to 80%, and which occurred 85~168 days prior to the M6.1 earthquake in Hualien, Taiwan. An inspiration could be provided for obtaining effective earthquake precursor anomalies by monitoring the change in ion concentration in thermal springs.

Sigmoidal normal faults and evidence for vertical-axis block rotation in an oblique convergent margin: A 3D seismic example from offshore Colombia

Sigmoidal fold and fault geometries are typical kinematic indicators of strike-slip fault zones. We document kilometer-scale, normal faults with sigmoidal plan-view geometries within the dextral pull-apart Bahia Basin, at the rear of the obliquely convergent South Caribbean Deformed Belt, offshore Colombia. Using 3D seismic reflection data calibrated to wells, closely spaced, low-displacement, planar normal faults are mapped within the Miocene strata. A series of seismic horizontal (time) slices and computed seismic attributes are used to interpret the 3D configuration of these faults. The closely spaced faults display an east–west trend with a progressive rotation into a northwest–southeast trend. In map view, the fault traces curve toward their tips, describing a sigmoidal-Z geometry that terminates at discrete northeast–southwest-trending fault zones. The structures observed may correspond to either tension fractures, which form theoretically at 45°, or antithetic shear fractures with normal displacement formed at 50°–70° to the boundaries of a dextral shear zone. These scenarios lead to a clockwise block rotation of between 20° and 40° within the shear zone. This study shows the first example of vertical-axis block rotations observed offshore in the western end of the South Caribbean margin and is an important example of the use of 3D seismic data to identify rotations where paleomagnetic studies are not available.

Influence of source and site effects of 2003 Tokachi-oki earthquake on the generation of high PGA in the near-fault zones

Source and site effects of 2003 Tokachi-oki earthquake (Japan, Mw~8.3) and their influence on the distribution of peak ground accelerations (PGA) in the near-fault zones are studied. Based on records of KiK-net vertical arrays, models of soil behavior are constructed, i.e. vertical distributions of stresses and strains induced in soil layers by strong motion. The method is used suggested by Pavlenko and Irikura and previously applied for studying soil behavior during 1995 Kobe, 2000 Tottori, and 2011 Tohoku earthquakes. During the Tokachi-oki earthquake, we did not find a widespread nonlinearity of soft soil behavior. Manifestations of soil nonlinearity were observed at sites closest to the source; at remote sites where high PGA were recorded, soil behavior was virtually linear, and shear moduli in soils increased till the moments of the highest intensity of motion, then decreased. The shapes of acceleration time histories at remote sites indicate directivity effects: seismic waves radiated by the crack tip during its propagation along a section of the fault plane came to the stations simultaneously. Soil hardening occurred at these sites that increased amplification and PGA on the surface. Similar effects were observed during 2011 Tohoku earthquake; evidently, they can occur during future strong earthquakes.