Search

The M s 5.5 earthquake struck on 24 October 2023, in Subei County, Gansu Province, China, occurring along the eastern segment of the Altyn Tagh fault. It raises the question of whether this earthquake is linked to the ongoing shortening slip rate along this segment or triggered by other seismic events. Analyzing the fault geometry of the Subei earthquake and understanding the significance of the weakening activity rate for seismic hazards in neighboring regions is crucial. The surface deformation from small- and medium-sized earthquakes (magnitudes less than M w5.5) is often subtle, and the coseismic deformation detected by interferometric synthetic aperture radar (InSAR) is vulnerable to atmospheric disturbances, leading to significant measurement errors. Moreover, inaccuracies in the regional crustal velocity structure can cause errors in earthquake localization based on seismic data. These challenges complicate the establishment of a rupture model for seismogenic faults and hinder the inversion of fault slip models. To overcome these limitations, we employed the time-series InSAR stacking method and aftershock relocation to determine the fault geometry of the Subei earthquake. A two-step inversion method was utilized to ascertain both the fault geometry and slip distribution. Our modeling indicates that the 2023 Subei earthquake had a thrust mechanism with a component of strike-slip. The rupture did not reach the surface, with the maximum fault slip measuring 0.45 m at a depth of 2.5–3.5 km. The fault dips westward, and the moment magnitude is calculated at 5.4. This earthquake is associated with the ongoing weakening of the left-lateral strike-slip rupture along the Altyn Tagh fault in the Subei region. Furthermore, retrograde thrust tectonics significantly contribute to the absorption of accumulated stress during this process.Our findings highlight the potential of utilizing time-series InSAR images to enhance earthquake catalogs with geodetic observations, offering valuable data for further studies of the earthquake cycle and active tectonics. This approach is also applicable in other tectonically active regions, enhancing understanding of seismic hazards and risk assessment. [ABSTRACT FROM AUTHOR]

Fault location and geometry are the most fundamental input data in seismic hazard analysis, the ultimate aim of which is to mitigate damage from future large earthquakes. In regions prone to large earthquakes or where cumulative deformation by multiple earthquake events are well expressed in the landscape, fault models are constructed primarily by (1) identifying active fault traces, mapped mostly by the surface ruptures associated with large earthquakes; (2) simplifying fault traces while capturing their geometrical characteristics; and (3) segmenting the simplified geometry, given that a single earthquake does not always rupture the entire length of a fault system. In slowly deforming regions, however, the construction of fault models is challenging, even though geologic records of large earthquakes exist, because of the lack of clear active fault traces. Indeed, surface-rupturing earthquakes may not be part of the historical periods owing to their long recurrence time of thousands of years or more. Nevertheless, seismic hazard analysis is required for densely populated and industrial areas in slowly deforming regions, such as South Korea. On the basis of criteria established previously for determining segmentation geometry in fault models, here we propose a methodology for identifying the segmentation geometry of strike-slip fault systems in slowly deforming regions. In terms of the criteria used to identify segment boundaries, we examine along-fault variations not only in fault geometry but also in fault-surrounding lithology and fault-related geomorphic features. We test the methodology for assessing the fault segmentation geometry in a case study of the Yangsan Fault, which is one of the most active seismogenic strike-slip faults on the Korean Peninsula. Results show that the ∼200 km length of the Yangsan Fault on land consists of 12 to 15 distinct fault segments. We discuss how models of fault segmentation geometry are able to improve seismic hazard analysis in regions that have not experienced surface-faulting earthquakes in historical period. [ABSTRACT FROM AUTHOR]

Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably‐dipping planar and variably‐curved nonplanar megathrusts using the volumetric, high‐order accurate code tandem to account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly‐dipping megathrusts to host larger earthquakes than steeply‐dipping ones. Under elevated pore pressure and less strongly velocity‐weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly‐dipping and sharply‐curved megathrusts host multi‐period supercycles of slow‐to‐fast, small‐to‐large slip events under higher effective stresses and more strongly velocity‐weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second‐order effects from structurally‐derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co‐ and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults. Plain Language Summary: Subduction zones, where one tectonic plate dives beneath another, generate the largest earthquakes worldwide. Our study investigates how the shape and tilt of these large offshore underground faults, termed "megathrusts," may determine the size of large earthquakes, how often they happen, and how similar or different subsequent events are. By creating computer simulations of earthquakes in subduction zones, we found that the angles and dimensions of the megathrust may set a limit on how big an earthquake can get. We also find that the presence of bends or curves along these faults can make earthquakes more unpredictable, sometimes leading to more variable series of smaller quakes before the biggest one hits. Our findings may help explain why some areas near subduction zones are prone to larger or more frequent earthquakes than others. Understanding these patterns can improve our ability to prepare for these natural disasters, potentially reducing their damaging effects on nearby communities and infrastructure. Key Points: Systematic earthquake cycle simulations reveal how subduction zone geometry controls megathrust earthquake size and timingDip and curvature affect characteristic slip style: periodic uniform ruptures versus supercycles of variably‐sized slow‐to‐fast eventsSeismogenic zone dip and dimensions limit maximum earthquake size; curvature‐linked stress and strength heterogeneity modulates recurrence [ABSTRACT FROM AUTHOR]

The Meishan-Chiayi area of western Taiwan has a large probability of producing a major earthquake in the near future. Historically, one of the largest and most damaging of Taiwan’s earthquakes occurred there. It is, therefore, important to have a well-constrained upper crustal 3-D shear-wave velocity model that can be used to accurately determine ground motion predictions and fault geometry models used in seismic hazard and risk modelling. In this study, we carried out an ambient noise tomography experiment using 100 seismometers deployed with a ∼2 km spacing on a 20 by 20 km grid. The reliable periods of phase velocity from Rayleigh waves are 0.6 to 6.8 s, providing a well-resolved Vs structure from the surface to a depth of around 4 km. The velocity model displays a prominent, roughly northeast-striking change in Vs that follows the projected surface trace of the blind Chiayi thrust. The uplift of relatively higher Vs rocks in its hanging wall, together with a negative to positive change in dVs suggests that it dips gently eastward across the study area. A northward thickening of the lower Vs crust, together with a high negative dVs in the north of the study area is related to an increased thickness of foreland basin rocks across the Meishan fault. The Vs and dVs models provide reasonable evidence that the Meishan fault can be traced at a high angle from its surface rupture to the base of the model at 4 km depth. It cuts the Chiayi thrust.

The Ms 5.5 earthquake struck on 24 October 2023, in Subei County, Gansu Province, China, occurring along the eastern segment of the Altyn Tagh fault. It raises the question of whether this earthquake is linked to the ongoing shortening slip rate along this segment or triggered by other seismic events. Analyzing the fault geometry of the Subei earthquake and understanding the significance of the weakening activity rate for seismic hazards in neighboring regions is crucial. The surface deformation from small- and medium-sized earthquakes (magnitudes less than Mw5.5) is often subtle, and the coseismic deformation detected by interferometric synthetic aperture radar (InSAR) is vulnerable to atmospheric disturbances, leading to significant measurement errors. Moreover, inaccuracies in the regional crustal velocity structure can cause errors in earthquake localization based on seismic data. These challenges complicate the establishment of a rupture model for seismogenic faults and hinder the inversion of fault slip models. To overcome these limitations, we employed the time-series InSAR stacking method and aftershock relocation to determine the fault geometry of the Subei earthquake. A two-step inversion method was utilized to ascertain both the fault geometry and slip distribution. Our modeling indicates that the 2023 Subei earthquake had a thrust mechanism with a component of strike-slip. The rupture did not reach the surface, with the maximum fault slip measuring 0.45 m at a depth of 2.5–3.5 km. The fault dips westward, and the moment magnitude is calculated at 5.4. This earthquake is associated with the ongoing weakening of the left-lateral strike-slip rupture along the Altyn Tagh fault in the Subei region. Furthermore, retrograde thrust tectonics significantly contribute to the absorption of accumulated stress during this process.Our findings highlight the potential of utilizing time-series InSAR images to enhance earthquake catalogs with geodetic observations, offering valuable data for further studies of the earthquake cycle and active tectonics. This approach is also applicable in other tectonically active regions, enhancing understanding of seismic hazards and risk assessment.

One of the three narrow rift belts that mark the southern end of the East African Rift System (EARS) intersects the Makgadikgadi Basin of northeastern Botswana. Although tectonic activity in the region is known to have influenced the evolution of these pans, the interrelationship between shoreline geometry, fault strikes, and the intersection of the underlying tectonic terranes has yet to be fully realized. We analyzed faults and subsurface structures in the region of the pans using a field investigation in combination with satellite imagery and geophysical data, to constrain the influence that the regional tectonic regime has had on the formation of the present-day pan geometry. We find that pan shorelines are controlled by the intersection of three preferred fault orientations which can be understood in the context of the “older” terranes they overlie, namely the Magondi Belt and the Limpopo Belt. We propose that the pronounced curvature of the southern Magondi Belt has influenced the eastward curvature of the rift-related faults and was likely produced by the impingement of the developing fold belt on the Zimbabwe Craton. Furthermore, limited focal mechanism solutions data from earthquakes north and south of the pans suggests a change in regional extension direction from NW-SE to NE-SW. Determining the relationship between these fault orientations and the underlying tectonic terrains is an important step in understanding the formation of the Makgadikgadi Basin, and more broadly the current tectonic regime of Botswana. The evidence of fault-controlled shorelines within an evaporitic environment may also have implications for regional groundwater activity.

A devastating earthquake with moment magnitude 7.5 occurred in the Noto Peninsula in central Japan on 1 January 2024. We estimate the rupture evolution of this earthquake from teleseismic P‐wave data using the potency‐density tensor inversion method, which provides information on the spatiotemporal slip distribution including fault orientations. The results show a long and quiet initial rupture phase that overlaps with regions of preceding earthquake swarms and associated aseismic deformation. The following three major rupture episodes evolve on segmented, differently oriented faults bounded by the initial rupture region. The irregular initial rupture process followed by the multi‐scale rupture growth is considered to be controlled by the preceding seismic and aseismic processes and the geometric complexity of the fault system. Such a discrete rupture scenario, including the triggering of an isolated fault rupture, adds critical inputs on the assessment of strong ground motion and associated damages for future earthquakes. Plain Language Summary: On 1 January 2024, a moment magnitude 7.5 earthquake occurred in the northern Noto Peninsula, Japan. The strong ground motion and tsunami associated with the earthquake caused severe damage to buildings and infrastructure, resulting in at least 245 causalities in the affected areas. The Noto Peninsula is affected by northwest‐southeast compression, and active reverse faults are known along the northern coast of the peninsula and its offshore region. Before the 2024 earthquake, the source region experienced long‐lasting earthquake swarm activity, which is a set of seismic events without an obvious mainshock‐aftershock pattern. Our seismological analysis found that there was a 10‐s‐long initial rupture episode around the hypocenter that overlapped with the earthquake swarm region. The initial rupture was followed by a series of three different rupture episodes on differently oriented fault segments. This earthquake highlights a multi‐scale rupture growth across a segmented fault network after a very quiet initial rupture process that was controlled by the preceding earthquake swarms and associated aseismic deformation related to fluid injection from depth. The rupture process advances our understanding of earthquake source physics and can lead to a better assessment of future earthquake hazards. Key Points: The 2024 Mw 7.5 Noto Peninsula earthquake involves a multi‐segmented rupture sequence on differently oriented faultsThe long and quiet initial rupture domain coincides with the preceding earthquake swarm regionFluid‐induced earthquake swarms and a segmented fault network control the complex earthquake rupture growth [ABSTRACT FROM AUTHOR]

Utilizing an updated dataset of 145 GNSS surface velocities, this study examines the fault slip rate and fault geometry along the Main Himalayan Thrust (MHT) in the central Himalaya. Employing a Bayesian inversion model, the present analysis reveals that the upper portion of the MHT ramp exhibits full locking, while the lower flat displays creeping motion. The estimated locking depth and fault depth of MFT range from 4.3 ± 2.6 km to 9.7 ± 2.2 km and 13.5 ± 3.1 km to 15.8 ± 1.9 km, respectively, along the central Himalaya. Further, the slip rate along the transition zone lies in the range of 1.4 ± 0.8 mm/yr to 2.7 ± 0.5 mm/yr. Considering the amount of uncertainties as ~1–2 mm/yr in GNSS velocities, the study suggests that the transition zone along the middle flat of the MHT also exhibits locking behavior. Thus, the estimated locking depth extends to ~15.0 km down-dip and covers a horizontal distance of ~90 km (locking line) on the surface, reaching the foothills of the Higher Himalaya. Furthermore, along the deeper flat of the MHT, the slip rate ranges from 19.4 ± 2.5 mm/yr in the west to 12.8 ± 1.6 mm/yr in the east along Nepal Himalaya. The analysis also calculates the slip deficit rate along the MHT fault plane, revealing values of ~15.1 mm/yr in western Nepal, ~12.7 mm/yr in central Nepal, and ~10.6 mm/yr in eastern Nepal. These slip deficit rates across different segments of central Nepal indicate the potential for large earthquakes in the region. The results are further supported by a resolution test using a checkerboard synthetic model, demonstrating the capability of the GNSS network to capture the slip rate along the MHT. These findings inevitably contribute to a comprehensive assessment of the seismic hazard potential in the central Himalayan region. [ABSTRACT FROM AUTHOR]

Fault geometric heterogeneities such as roughness, stepovers, or other irregularities are known to affect the spectra of radiated waves during an earthquake. To investigate the effect of normal stress heterogeneity on radiated spectra, we utilized a poly(methyl methacrylate) (PMMA) laboratory fault with a single, localized bump. By varying the normal stress on the bump and the fault-average normal stress, we produced earthquake-like ruptures that ranged from smooth, continuous ruptures to complex ruptures with variable rupture propagation velocity, slip distribution, and stress drop. High prominence bumps produced complex events that radiated more high frequency energy, relative to low frequency energy, than continuous events without a bump. In complex ruptures, the high frequency energy showed significant spatial variation correlated with heterogeneous peak slip rate and maximum local stress drop caused by the bump. Continuous ruptures emitted spatially uniform bursts of high frequency energy. Near-field peak ground acceleration (PGA) measurements of complex ruptures show nearly an order-of-magnitude higher PGA near the bump than elsewhere. We propose that for natural faults, geometric heterogeneities may be a plausible explanation for commonly observed order-of-magnitude variations in near-fault PGA.

Compared with the landslides in the general gravity environment, the earthquake-induced landslides are significantly different in formation mechanisms and kinetic characteristics. Under the normal and rainfall conditions, the fault fracture zone, as the discontinuous structural plane of the slope, often adversely affect the stability of the slope. Under the earthquake condition, the fault fracture zone within the slope has a limited filtering effect, which could weaken the seismic response of the slope. To investigate the influence of the reverse fault’s geometry on the slope’s seismic response, we took the Niumiangou landslide, the Woqian landslide, the Xiejiadianzi landslide and the Donghekou landslide as reference objects and generalized the geological model of the fault-crossing landslide in this study. The seismic response of slopes with faults of different widths, dips and positions are simulated using the 3DEC discrete element software. The simulation results show that (1) as the fault dip angle increases, the peak value of the total displacement of the slope and the peak acceleration of the slope surface show an increasing trend, and the slope stability is worse. (2) The peak acceleration of the monitoring point at the top of the slope is generally greater than that at the bottom and waist of the slope. As the width of fault fracture zone increases, the effect on the seismic response of the slope becomes obvious.(3) The presence of faults facilitates the probability of slope instability. When the fault is located at the top of the slope, the variation of the seismic response with the dip angle and the fault width shows a more obvious regularity. This study can provide a theoretical basis for further revealing the impact of fault fracture zone on the stability of slopes under the earthquake condition.

The potential of high-temperature (>150 °C) geothermal systems in crustal fault zones (fault cores and hundreds of meters wide networks of interconnected fractures in the damage zone) is underestimated. Based on numerical models, we show that topography-driven, poroelasticity-driven as well as buoyancy-driven forces play a significant role in the establishment of shallow (1-4 km) thermal anomalies in fault zones. We investigate the role of permeability, topography, fault dip, tectonic regime and fault geometry on the amplitude of thermal anomalies. Favorable conditions include: (i) a damage zone thickness > 100m, (ii) a minimum cumulative displacement of 100-150m and (iii) fault zone lengths of at least one kilometer. Based on these parameters, we propose new potential targets for the geothermal exploration of fault zones inWestern Europe. [ABSTRACT FROM AUTHOR]

The Tieluzi Fault is the largest structure in the East Qinling Mountains, and is considered to be the easternmost continuation of the Altyn Tagh‐Haiyuan‐Qinling Fault System (AHQFS) that allows the eastward extrusion of the Tibetan Plateau and South China Block. We studied the fault geometry and kinematics of the Tieluzi Fault using field investigations, detailed interpretations of high‐resolution satellite imagery and digital elevation models, and late Quaternary dating methods. Paleoseismic investigations indicate that the most recent earthquake along the Tieluzi Fault occurred before 1,500–1,300 cal. BP. Geological and geomorphological observations show that segments west of Lushi County are more active than those to the east. The spatial variations in tectonic activity along the Tieluzi Fault are interpreted to be related to four possible mechanisms: strike change, discontinuity, intersection, and branch. The late Quaternary left‐lateral slip rate is determined to be 0.9 ± 0.1 mm/yr on the Tieluzi Fault. The prominent left‐lateral faulting along the Tieluzi Fault suggests that most of the left‐lateral displacement along the eastern AHQFS has been accommodated by the Tieluzi Fault, which forms the most frontier of the eastward expansion of the Tibetan Plateau. Furthermore, we suggest that the left‐lateral faulting in the East Qinling Mountains is a response to relative eastward motion of the South China block pushed by the Tibetan Plateau with respect to the North China Plain Block. Also, our results indicate that the Tibetan Plateau has undergone a stepwise eastward expansion. Key Points: Paleoseismic investigations indicate that the most recent earthquake along the Tieluzi Fault occurred before 1,500–1,300 cal. BPThe late Quaternary left‐lateral slip rate is determined to be 0.9 ± 0.1 mm/yr on the Tieluzi FaultThe prominent left‐lateral faulting along the Tieluzi Fault is a far‐field response for the eastward expansion of the Tibetan Plateau [ABSTRACT FROM AUTHOR]

We use a comprehensive dataset of field observations, high spatial resolution drone orthomosaics and digital terrain models (DTMs) to map, quantify and characterize the extensive ground fracturing related to the 2021 seismo-tectonic and volcanic activity in the Reykjanes Peninsula (Iceland). The dataset, spans an area of about 30 km 2 , where we map nearly 20 000 ground cracks with metric to decametric lengths and centimetric extensional offsets, revealing a dominant dextral shear, in agreement with published seismic data. Although striking in a direction similar to the volcanic systems in the Reykjanes Peninsula (N030–040), most fractures appear as en-échelon structures globally aligned along N-S-striking fault zones up to 3–4 km long. By examining the timing of ground fracturing through repeated field observations, seismic data and InSAR images, we associate a fracture zone with most earthquakes of M ω ≥ 5.0 that occurred in the month preceding the March 2021 Fagradalsfjall eruption. We describe three preexisting N-S fault zones, with fault segments that were reactivated up to three times during the pre-eruptive seismic activity, while the magma intrusion did not trigger graben-related ground fractures typically observed during magmatic injections. Our depiction of a system dominated by strike-slip tectonic features helps in understanding the geometry and bookshelf-mode of tectonic activity along a diffuse and highly oblique extensional plate boundary. Evidence of transient fracturing is typically quickly lost because of erosion or lava flow burial, highlighting a potential under-representation of diffuse fracturing when studying old tectonic and volcanic systems. [ABSTRACT FROM AUTHOR]

Physics‐based earthquake simulators have been developed to overcome the relatively short duration and incompleteness of historical earthquake and paleoseismic records, respectively. These simulators have the potential to be a useful addition to seismic hazard assessment as they produce millions of synthetic earthquakes over thousands to millions of years using predefined fault geometries and slip rates. Due to the sparsity of fault data and the computational expense of the modeling, it is common to simplify earthquake simulator input parameters. Fault surfaces are typically characterized by simplified mapped traces and constant dip with depth, with slip rates that are often assumed to be uniform over the fault plane. This study uses the well‐defined 3D geometry of an active normal fault in offshore Aotearoa New Zealand, the Cape Egmont Fault, to demonstrate the impact of using nonuniform slip rates and nonplanar fault geometries on the resulting synthetic earthquakes from the earthquake simulator RSQSim. Adopting variable slip rates that decrease to zero along the fault tipline (rather than uniform slip rates) reduces unrealistically high nucleation rates of seismicity along fault edges. Introduction of complex 3D fault geometries, including fault segmentation and bends on kilometer scales, and variable slip rates produces less characteristic earthquake populations, increasing the number of M6‐7 moderate‐large magnitude events. Incorporation of variable fault geometries and slip rates in physics‐based simulators may modify the seismic hazards estimated from synthetic earthquake catalogs. Plain Language Summary: Earthquake simulators are used to study and anticipate earthquakes. These simulators produce millions of synthetic earthquakes over thousands to millions of years using predefined fault properties. We use a well‐defined active normal fault in offshore Aotearoa New Zealand, to examine the impact of nonuniformly distributed fault slip rates and nonplanar fault geometries on the resulting synthetic earthquake catalogs. We show that adopting variable fault slip rates, with a gradual decrease in slip rate toward fault‐surface edges, reduces the unrealistic occurrence of high concentrations of earthquakes along the fault tip line. The study also shows that utilization of complex 3D fault geometries together with variable slip rates produces populations of earthquakes in a less predictable manner, increasing the number of moderate to large magnitude events. Therefore, details of the input fault geometries and slip rate distributions can significantly affect the resulting synthetic earthquake catalogs from physics‐based earthquake simulators with significant implications for seismic hazard applications. Key Points: Variable, rather than constant, slip rates reduce high seismicity rates along the fault edges and can improve the depth distributions of eventsIntroduction of complex 3D fault geometries produces less characteristic earthquake populationsVariable fault geometries and slip rates in physics‐based simulators can improve the synthetic earthquake catalogs [ABSTRACT FROM AUTHOR]

Copyright of Hydrogeology & Engineering Geology / Shuiwendizhi Gongchengdizhi is the property of Hydrogeology & Engineering Geology Editorial Office and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

The Xiluodu (XLD) reservoir is the second largest reservoir in China and the largest in the Jinsha River basin. The occurrence of two M > 5 earthquakes after reservoir impoundment has aroused great interest among seismologists and plant operators. We comprehensively analyzed the seismicity of the XLD reservoir area using precise earthquake relocation results and focal mechanism solutions and found that the seismicity of this area was weak before impoundment. Following impoundment, earthquake activity increased significantly. The occurrence of M ≥ 3.5 earthquakes within five years of impoundment also appear to be closely related to rapid rises and falls in water level, though this correlation weakened after five years because earthquake activity was far from the reservoir area. Earthquakes in the XLD reservoir area are clustered; near the dam (Area A), small faults are intermittently distributed along the river, while Area B is composed of multiple NW-trending left-lateral strike-slip faults and a thrust fault and Area C is composed of a NW-trending left-lateral strike-slip main fault and a nearly EW-trending right-lateral strike-slip minor fault. The geometries of the deep and the shallow parts of the NW-trending fault differ. Under the action of the NW-trending background stress field, a series of NW-trending left-lateral strike-slip faults and NE-trending thrust faults in critical stress states were dislocated due to the stress caused by reservoir impoundment. The two largest earthquakes in the XLD reservoir area were tectonic earthquakes that were directly triggered by impoundment.

Rupture dynamics along a heterogeneous fault is studied in the framework of the recently developed damage‐breakage rheological model. The model utilizes a fault structure accommodating shear and wear in a heterogeneous weak fault‐core a few centimeters thick, separating two blocks with damage level gradually decreasing toward the host rock. We first demonstrate the similarity of the static stress field around a fault, represented by a three‐body tribosystem, to the one developed around a rough frictional interface. We show that both models predict heterogeneous stress field pattern. We then apply simulations of rupture on heterogeneously damaged faults. We show that increasing initial heterogeneity amplitudes is associated with smaller events with lower slip rates. The simulations further allow us to quantify the amount of accumulated damage correlative with wearing. During the propagating rupture, the strength of the fault‐core evolves, leading to higher wear generation along relatively strong zones or barriers. The total wear production in a given event is strongly dominated by the initial damage heterogeneity. Processing of synthetic seismograms shows excess energy radiation of high‐frequency seismic waves comparing to the expected radiation from planar faults. This radiation is enhanced with the increase in the variability of fault strength heterogeneity. Further calculations of the scaled energy show good fit with previous seismological and laboratory observations and demonstrate the impact of fault heterogeneity on radiated energy during dynamic rupture. Therefore, initial fault heterogeneity, manifested by fault‐core strength or by geometrical irregularity controls many aspects of earthquake rupture, including slip displacement and velocity. Plain Language Summary: Earthquakes are triggered when the stress build‐up on rocks exceeds the resistance to fracture of fault material. When faulting occurs, energy is generated in the fault zone and radiated through the rock body to the surface. The ground motions and their durations strongly depend on the earthquake rupture characteristics such as the amount of displacement, velocity, and directivity. Many studies modeled natural faults as planar frictional interfaces that failed and were displaced during earthquakes. Yet, geological observations and seismological measurements demonstrated that this is probably a naive picture, and fault zones are very heterogeneous, including complicated geometry and variation of strength along and across them. Here we examine the effect of such initial heterogeneities on earthquake and show that the heterogeneity of the rock strength in the fault core itself strongly affects earthquake dynamics. In particular, larger magnitude earthquakes are predicted for the more homogenous fault. We also find that higher fault heterogeneity would result in excess high‐frequency radiation. Considering these effects on ground motion around the fault, we show that seismological records may capture these differences in fault characteristics during earthquakes. Key Points: Initial damage level and its heterogeneity along the fault core control dynamics and seismic moment of simulated rupture eventsWear production is affected by the initial damage distribution. During slip, higher wear accumulation is calculated for barriers along faultsHigh variability of fault damage heterogeneity significantly enhances high‐frequency radiation [ABSTRACT FROM AUTHOR]

Surface traces of earthquake faults are complex and segmented on multiple scales. At seismogenic depth the detailed geometry of faults and earthquake rupture is mainly constrained by earthquake locations. Standard earthquake locations are usually too diffuse to constrain multi-scale fault geometry, while differential-timing relocation mainly improves finest scale precision. NLL-SSST-coherence, an enhanced, absolute-timing earthquake location procedure, iteratively generates traveltime corrections to improve multi-scale precision and uses waveform similarity to improve fine-scale precision. Here we apply NLL-SSST-coherence to large-earthquake sequences and background seismicity along strike-slip faults in California. Our relocated seismicity at seismogenic depth along major fault segments and around large-earthquake ruptures often defines smooth, planar or arcuate, near-vertical surfaces across the sub-km to 10’s of km scales. These results show that multi-scale smooth fault segments are characteristic of major, strike-slip fault zones and may be essential to large earthquake rupture. Our results suggest that smoothness and curvature of faults influences earthquake initiation, rupture, rupture direction and arrest, and can define earthquake hazard. The results corroborate that surface traces of strike-slip fault zones reflect complex, shallow deformation and not directly simpler, main slip surfaces at depth, and support use of planar or smoothly curved faults for modeling primary earthquake rupture.

Geologic reconstructions of the Main Himalayan Thrust in Nepal show a laterally extensive midcrustal ramp, hypothesized to form the downdip boundary of interseismic locking. Using a recent compilation of interseismic GPS velocities and a simplified model of fault coupling, we estimate the width of coupling across Nepal using a series of two-dimensional transects. We find that the downdip width of fault coupling increases smoothly from 70 to 90 km in eastern Nepal to 100–110 km in central Nepal, then narrows again in western Nepal. The inferred coupling transition is closely aligned with geologic reconstructions of the base of the midcrustal ramp in central and eastern Nepal, but in western Nepal, the data suggest that the location is intermediate between two proposed ramp locations. The result for western Nepal implies either an anomalous coupling transition that occurs along a shallowly dipping portion of the fault or that both ramps may be partially coupled and that a proposed crustal-scale duplexing process may be active during the interseismic period. We also find that the models require a convergence rate of 15.5 ± 2 mm/year throughout Nepal, reducing the geodetic moment accumulation rate by up to 30% compared with earlier models, partially resolving an inferred discrepancy between geodetic and paleoseismic estimates of moment release across the Himalaya.

The East Anatolian Fault Zone (EAFZ) is a left-lateral transform fault zone located between the Anatolian and Arabian plates. In this study, in order to image the upper crustal structure beneath the eastern segments of EAFZ, 3D seismic velocity variations are computed using local earthquake tomography. The initial catalog for the tomography process consists of 2200 well-located earthquakes recorded at 49 seismic stations around the study region between 2007 and 2020. 1D initial velocity model is constructed based on previous studies in the region. The maximum number of iterations and the velocity perturbations which sustain the linearity of the inversion are determined based on the detailed tests. Reliable zones of the final model are decided based on the Derivative Weighted Sum and Hit Count distribution. The resulting velocity model displays a clear velocity contrast across the surface trace of the EAFZ down to a depth of 12 km. While the Anatolian side of the fault displays higher velocities associated with the ophiolitic units in the region, the south of the fault zone is represented by lower velocities due to sedimentary deposits. The vertical cross-sections of tomographic models show a north dipping fault between Palu and Çelikhan. The complete earthquake catalog is relocated using the 3D velocity model. Together with the obtained velocity model, the relocated hypocenters indicate that the dip of the EAFZ is not uniform, the Palu segment dips to the north with an angle of ~80°, while the Pütürge and Erkenek segments dip to the north with a lower angle of ~60-70°. [ABSTRACT FROM AUTHOR]

A precise finite‐fault model including the fault geometry and slip distribution is essential to understand the physics of an earthquake. However, the conventional linear inversion of geodetic data for a finite‐fault model cannot fully resolve the fault geometry. In this study, we developed a Bayesian inversion framework that can comprehensively solve a non‐planar fault geometry, the corresponding fault slip distribution with spatially variable directions, and objective weighting for multiple data types. In the proposed framework, the probability distributions of all the model parameters are sampled using the Monte Carlo method. The developed methodology removes the requirement for manual intervention for the fault geometry and data weighting and can provide an ensemble of plausible model parameters. The performance of the developed method is tested and demonstrated through inversions for synthetic oblique‐slip faulting models. The results show that the constant rake assumption can significantly bias the estimates of fault geometry and data weighting, whereas additional consideration of the variability of slip orientations can allow plausible estimates of a non‐planar fault geometry with objective data weighting. We applied the method to the 2013 Mw 6.5 Lushan earthquake in Sichuan province, China. The result reveals dominant thrust slips with left‐lateral components and a curved fault geometry, with the confidence interval of the dip angles being between 20°–25° and 56°–58°. The proposed method provides useful insights into the scope of imaging a non‐planar fault geometry, and could help to interpret more complex earthquake sources in the future. Plain Language Summary: An earthquake originates from a rapid relative movement between two blocks on both sides of a fault. The earthquake can produce ground deformation, which can be measured by geodetic instruments deployed on the ground. Geophysicists can thus utilize the geodetic measurements to invert the earthquake attributes, including the fault geometry and the associated slips (i.e., the zone and amount of the relative movement between two blocks). Simultaneous estimation of the fault geometry parameters and the slip parameters is a nonlinear problem. The conventional inversion of geodetic data uses a linear adjustment method that cannot fully resolve the fault geometry. In addition, geodetic measurements have different precisions in terms of the types of instruments. We developed a nonlinear inversion based Bayesian inference framework to solve both these problems. The proposed method can correctly interpret an earthquake source with a non‐planar fault geometry. We applied the method to the 2013 Mw 6.5 Lushan earthquake in southwest China and determined a statistical estimate of the fault geometry and slips. The results indicate that the proposed method could help to improve our knowledge of complicated earthquake sources in the future. Key Points: We develop a Bayesian finite‐fault inversion for non‐planar fault geometry and spatially variable slip amplitude and rakeAccounting for variable directions of fault slips allows to better constrain non‐planar fault geometry and tune objective data weightingAccounting for the non‐planar fault geometry, the 2013 Lushan earthquake modeling reveals a dominant thrust slip with sinistral component [ABSTRACT FROM AUTHOR]

Analysis of TanDEM‐X and Shuttle Radar Topography Mission (SRTM) data reveals geomorphic evidence for 292 fault‐propagation fold scarps across the Miocene Nullarbor and Pliocene Roe Plains in south‐central Australia. Vertical displacements (VD) are determined using topographic profiling of a subset (n = 48) of the fold traces. Fault dips (mean = 44 +16/−14° at 1σ) are estimated from seismic reflection data; the mean dip is assigned to faults with unknown dip and combined with VD to estimate net displacements (ND) and average net displacements (AD) for each fault. AD exceeds single‐event displacements estimated from fault‐length scaling regressions, indicating the identified faults have hosted multiple earthquakes. Combining AD with (i) faulted surface ages (Nullarbor ~10–5 Ma, Roe ~2.5 Ma), (ii) ages of faulted erosional–depositional features (e.g. relic Late Miocene dune fields and Pliocene paleochannels), and (iii) onset of the neotectonic regime in Australia at ~10 Ma yields average slip rates from <0.1 m Myr−1 to >17 m Myr−1 (mean = 1.1 m Myr−1). Summation of displacements across faults yields crustal horizontal shortening rates lower than geodetically detectable resolution (≤0.01 mm yr−1) since the Late Miocene. The ca. 10 Myr‐long record of neotectonic faulting on the Nullarbor Plain provides important insights into earthquake spatial–temporal behaviours in a slowly deforming intraplate continental region. [ABSTRACT FROM AUTHOR]

The 2008 Mw 7.9 Wenchuan earthquake is the most disastrous seismic event in China since 1976. Both field and seismological investigations suggest a multi‐stage coseismic rupture with most damage associated with the Yingxiu‐Beichuan Fault (YBF) of spatially variable fault strike and dip angles. To investigate the effect of fault geometric complexity on coseismic rupture and paleoseismic pattern on the YBF, we perform earthquake sequence modeling on 3D fault geometry in the framework of rate‐and‐state friction law. Our model produces a long‐term earthquake sequence with quasi‐regular recurrences of whole‐fault ruptures and segmented ruptures. The along‐strike rupture segmentation, earthquake propagation speed, and slip rate are mainly controlled by along‐strike variations of fault dip and strike angles. Particularly, the YBF can be divided into two segments: the southern segment featured by large events (Mw > 8.0) and the northern segment with smaller events (Mw < 7.5), with recurrence intervals primarily determined by the tectonic loading rate. In a whole‐rupture event, the rupture speed mainly correlates with the fault dip angle; a smaller dip angle results in a wider seismogenic zone and higher rupture speeds, whereas the small‐scale variation in rupture speed is regulated by the fault strike angle. The effect of strike variation on the total coseismic slip amount is more pronounced due to the large strike angle gradient along the YBF. After varying the slip vector to reflect the northward transition from the dominant thrust‐slip to strike‐slip faulting during the Wenchuan coseismic rupture, we obtain reasonably good agreement of model simulated coseismic surface displacements with GPS observations. Plain Language Summary: The Yingxiu‐Beichuan Fault (YBF) hosting the 2008 Mw 7.9 Wenchuan earthquake has a complex nonplanar interface. Seismic studies after the Wenchuan earthquake have unveiled a complex rupture process. Field trenching studies find a complex recurrence interval and rupture pattern of paleoseismic events. Here, we conduct a numerical simulation of earthquake sequences on the 3D YBF to study the effect of fault geometry on the coseismic rupture and paleoseismic behaviors. Our modeling results are quantitatively comparable with most field observations. Both whole‐fault ruptures and partial rupture events are produced. The larger events (magnitude >8.0) initiate near the southern end of the YBF. The smaller events (magnitude <7.5) initiates near the northern end of the YBF. Their repeating time is mainly controlled by the background tectonic loading rate. In the partial rupture events, coseismic rupture stops near Nanba area when encountering a geometrical barrier due to a sharp change of fault dip. While the lateral rupture propagation speed is controlled by variation of both fault dip and strike angles, coseismic slip distribution is more controlled by the variation of fault strike (lateral fault orientation) than dip. We find that the orientation of slip vector is important in coseismic rupture modeling of Wenchuan‐type event. Key Points: A long‐term earthquake sequence with quasi‐regular recurrences is simulated, including whole‐fault and segmented rupturesFault dip angle and hence the seismogenic zone width has first‐order control on the along‐strike rupture segmentationFault strike angle controls small‐scale variations in rupture propagation speed and slip rate, which affect the cumulative coseismic slip [ABSTRACT FROM AUTHOR]

Badlands are extensively eroded landscapes consisting of weakly consolidated deposits within highly dense drainage systems. Their controlling and shaping factors can differ in relation to various internal and external conditions and processes that are not always well understood. This study focuses on the development of a badland landscape affecting Miocene and Quaternary sand-clay sediments in the extensional tectonic regime of Western Turkey with a multidisciplinary approach. The area between Kula and Selendi towns exhibits a badland topography with extensively eroded surface features, deepened gullies within poorly consolidated, sand clay-sized sediments. The results of structural field mapping and morphometric analyses using a 5 m resolution DEM to study the role of structural control in the development of badlands are presented in this study. Field data analysis supported by the quantitative assessment of longitudinal gully profiles illustrates the role of pre-existing structures as faults, their orientation and geometry in net erosion-sedimentation and the development of deepened gully networks. Representative illustrations, field photographs and block diagrams are presented to show the relationship between the rock structure and badland landscape. The connection between the extensional tectonics, erosional dynamics and geomorphology point to a structurally-controlled landscape in the Kula badlands in Western Turkey. [ABSTRACT FROM AUTHOR]

Paleoseismic evidence from the late‐interseismic Alpine Fault suggests key section boundaries conditionally inhibit rupture. We utilize a year of data from a two‐part seismometer network (Dense Westland Arrays Researching Fault Segmentation) to characterize ∼7,500 earthquakes (−0.7 ≤ MLv ≤ 4.2) and ∼800 focal mechanisms, producing high‐resolution structural images of these boundaries to study effects of material and structural heterogeneities on mode‐switching rupture behavior. Lithologically‐controlled frictional behavior and crustal strength appear to influence lateral and vertical on‐fault seismicity distributions. Ultramafic hanging‐wall serpentinite and serpentinite‐related fault core minerals along the South Westland (SW) boundary, result in a locally shallow seismogenic cuttoff (∼8 km) and abundant on‐fault seismicity. Maximum horizontal compressive stress rotations (14° anti‐clockwise and 20° clockwise near the SW and North Westland (NW) boundaries, respectively, relative to the Central Section), coupled with spatially variable fault frictional properties, are more important than geometry alone in controlling Sections' relative frictional stability. Whereas the SW and Central Sections are well‐oriented for failure, the NW Section is severely misoriented compared with favorably oriented faults of the Marlborough Fault Zone, which possibly facilitate a preferred rupture route. Geometrically, a 40° dip change at the SW boundary may be accommodated either by a single through‐going fault plane ‐ a difficult geometry across which to obtain multi‐segment earthquakes when considering rupture dynamics ‐ or by a deeper vertical fault strand truncated by a shallower listric plane. Our new observations have implications for Alpine Fault rupture scenarios and highlight the need to consider a range of spatially heterogeneous, interdependent physical factors when evaluating controls on rupture segmentation. Plain Language Summary: New Zealand's Alpine Fault last ruptured in 1717 AD and the likelihood of it producing a major earthquake in the next 50 years is estimated to be 75%. The magnitude of the next earthquake will be strongly determined by the area of fault rupture. Two key localities along the fault, which correspond to the boundaries between sections of the fault with differing geometry, slip rate and frictional behavior, have been proposed to sometimes, but not always, act as barriers to fault rupture, limiting the sizes of the resulting earthquakes. However, it remains unclear what physical properties control this process. Using seismic data collected from dense seismometer networks, we construct detailed catalogs of small earthquakes occurring near these localities, at Martyr River and Inchbonnie, and precisely map the distribution of active seismicity near the fault. From these data we improve estimates of (a) the subsurface geometry of the fault and how it changes along length, (b) how seismicity is distributed with depth, to understand possible maximum depth of rupture in conjunction with observations of crustal temperatures and structure, and (c) the orientation of the stresses acting on the fault, to understand its stability to tectonic loading and its frictional behavior. Key Points: Dense microseismicity catalogs constrain fault structure, seismogenic depths and the stress field bookending the fault's Central SectionOn‐fault seismicity, a 40° dip change and a shallow, rheologically controlled cutoff depth are observed at the South Westland boundaryStress analyses indicate South and North Westland Sections are favorably and unfavorably oriented for frictional failure, respectively [ABSTRACT FROM AUTHOR]

The eastward expansion of the Tibetan Plateau has resulted in different earthquake types in the eastern Tibetan Plateau, but the mechanism remains unclear. Here, we construct a three-dimensional visco-elastoplastic finite element model considering the topography to investigate the influence of fault geometry and rheological heterogeneity on stress fields. In our best-fitting model, the minimum principal stress is nearly vertical around the southern Huya fault zone, which is adjacent to the Longmen Shan fault zone, due to the significant mid-lowerWE crust lateral rheological heterogeneity, and the thrust stress regime accounts for the reverse fault and thrust-dominated earthquakes. In this scenario, the eastward horizontal motion of the mid-lower crust is obstructed and facilitates thrust faulting, suggesting the limited eastward expansion of the Tibetan Plateau. In contrast, the northern Huya fault zone, one of the terminal branches of the East Kunlun fault, accommodates the continuous eastward extrusion of the East Kunlun fault, where the stress regime under a more homogenized crust favors the strike-slip faulting process, along with the dominant strike-slip earthquakes. Moreover, the best-fitting of stress regime explains the thrust-dominated 2008 Ms. 8.0 Wenchuan and 2013 Ms. 7.0 Lushan earthquakes on the Longmen Shan fault zone. Combining geophysical and geodetic observations and model analyses, we propose that the hybrid deformation mode in the eastern Tibetan Plateau is accommodated by upper crustal shear and thrusting deformation and mid-lower crustal thickening driven by the gravitational potential energy gradient. Our results elucidate the mechanism for differences in strong historical earthquakes and, more importantly, isolate the effect of fault geometry from those of heterogeneous viscosity on crustal deformation and stress heterogeneity in the eastern Tibetan Plateau. • We model the three-dimensional stress regime in the eastern Tibetan Plateau. • The stress regime well explains the spatial differences of earthquake types in Huya fault. • The best-fitting model isolates the effect of fault geometry from heterogeneous viscosity on crustal deformation. • The eastern Tibetan Plateau experiences the present-day limited eastward expansion. [ABSTRACT FROM AUTHOR]

Simple mechanical arguments suggest that slip along interlocked, rough faults, damages surrounding rocks. The same arguments require that the scale of secondary damage is proportional to the size of geometric irregularities along the main fault. This relationship could apply at all scales, but has, so far, been difficult to observe at the 10s to 100 s of km scales of large, natural faults, often because large-scale deformation is distributed across wide, complex plate-boundary fault systems, like the San Andreas Fault. The geometry and geology of another large-scale plate-boundary strike slip fault—the Queen Charlotte Fault (QCF)—is, in contrast, especially simple. Here, we show that observations of secondary deformation are well-aligned with predictions of stress variations caused by geometric irregularities along the QCF, suggesting a geometric relationship between primary fault geometry and secondary deformation. The analytic stress solution reveals that the highest stresses and highest likelihood of failure are confined to a zone of influence (ZOI) with a width quantified by ZOI = λ / 2 π , where λ is the wavelength of geometric variations along the main fault. This simple model is consistent with ∼100-km-scale observations along the QCF and can theoretically be used to predict the width of secondary deformation at all scales. • Fault geometry controls width of stress perturbations and zone of influence (ZOI). • ZOI scales linearly with wavelength of fault geometry. • Secondary deformation along Queen Charlotte Fault within predicted ZOI. [ABSTRACT FROM AUTHOR]

In this paper, faults, one of the most important causes of geohazards, were investigated from a kinematic and geometric viewpoint in the northern part of the Sistan suture zone (SSZ), which serves as the boundary between the Afghan and Lut blocks. Furthermore, field evidence was analyzed in order to assess the structural type and deformation mechanism of the research area. In the northern Birjand mountain range, several ~E–W striking faults cut through geological units; geometric and kinematic analyses of these faults indicate that almost all faults have main reverse components, which reveals the existing compressional stress in the study area. The northern Birjand mountain range is characterized by four main reverse faults with ~E–W striking: F1–F4. The F1 and F2 reverse faults have southward dips, while the F3 and F4 reverse faults have northward dips. Moreover, the lengths of the F1, F2, F3, and F4 faults are 31, 17, 8, and 38 km, respectively. These faults, with reverse components that have interactive relationships with each other, form high relief structures. The study area's main reverse faults, including F1 to F4, are extensions of the Nehbandan fault system, while their kinematics and geometry in the northern Birjand mountain range point to an N–S pop-up structure. [ABSTRACT FROM AUTHOR]

Abstract The geometries (i.e., dip angles) of active faults from the surface to the seismogenic zone are the most important factors used to evaluate earthquake ground motion, which is crucial for seismic hazard assessments in urban areas. In Osaka, a metropolitan city in Japan, there are several active faults (e.g., the Uemachi and Ikoma faults), which are inferred from the topography, the attitude of active faults in surface trenches, the seismic reflection profile at shallow depths (less than 2 km), and the three-dimensional distribution of the Quaternary sedimentary layers. The Uemachi and Ikoma faults are N–S-striking fault systems with total lengths of 42 km and 38 km, respectively, with the former being located ~ 12 km west of the latter; however, the geometries of each of the active faults within the seismogenic zone are not clear. In this study, to examine the geometries of the Uemachi and Ikoma faults from the surface to the seismogenic zone, we analyze the development of the geological structures of sedimentary layers based on numerical simulations of a two-dimensional visco-elasto-plastic body under a horizontal compressive stress field, including preexisting high-strained weak zones (i.e., faults) and surface sedimentation processes, and evaluate the relationship between the observed geological structures of the Quaternary sediments (i.e., the Osaka Group) in the Osaka Plain and the model results. As a result, we propose geometries of the Uemachi and Ikoma faults from the surface to the seismogenic zone. When the friction coefficient of the faults is ~ 0.5, the dip angles of the Uemachi and Ikoma faults near the surface are ~ 30°–40° and the Uemachi fault has a downward convex curve at the bottom of the seismogenic zone, but does not converge to the Ikoma fault. Based on the analysis in this study, the dip angle of the Uemachi fault zone is estimated to be approximately 30°–40°, which is lower than that estimated in the previous studies. If the active fault has a low angle, the width of the fault plane is long, and thus the estimated seismic moment will be large.

The 2021 Mw 7.3 Maduo earthquake revealed the significant seismic hazard of faults developed within the Bayan Har Block of eastern Tibet, China (e.g., the Kunlun Pass–Jiangcuo Fault). Relocated aftershock data are in good agreement with the Interferometric Synthetic Aperture Radar (InSAR) coseismic displacement field and field investigations. In this study, we used aftershock point cloud fitting to model the relocated aftershocks of the Maduo earthquake, and obtained the detailed geometry and characteristics of the seismogenic fault. Based on InSAR coseismic deformation, the geometric model of the seismogenic fault and its slip distribution were retrieved. The results show that this event was shallow (0–10 km) and characterized by sinistral strike-slip motion. We identified four asperities along the fault strike; the maximum slip of 4.84 m occurred on the eastern segment of the fault, in an area where the strike changed. The results suggest that the central segment of the main seismogenic fault is mature and smooth, while western and eastern segments are complex and immature.

The geometry and kinematics of active faults have a significant impact on their seismic potential. In this work, a structural characterization of the active Baza Fault (central Betic Cordillera, southern Spain) combining surface and subsurface data is presented. Two sectors are defined based on their surface geometry: a northern sector striking N–S to NNW–SSE with a narrow damage zone and a southern sector striking NW–SE with a wide damage zone. A kinematic analysis shows pure normal fault kinematics along most of the fault. Geometric differences between the northern and southern sectors are caused by i) a heterogeneous basement controlling the fault geometry at depth and in the cover; ii) different orientations of the Baza Fault in the basement with respect to the regional extension direction and iii) interaction with other active faults. We use this structural characterization to analyse the segmentation of the Baza Fault. According to segmentation criteria, the entire Baza Fault should be considered a single fault seismogenic segment. Consequently, the seismic potential of the fault is defined for a complete rupture. Magnitude for the Mmax event is calculated using several scale relationships, obtaining values ranging between Mw 6.6 and Mw 7.1. Recurrence times range between approximately 2,000 and 2,200 years for Mmax events and between 5,300 and 5,400 years for palaeo-events. A geodetic scenario modelled for an Mmax event of Mw 6.7 shows permanent vertical displacements of more than 0.40m and an overall WSW–ENE extension during entire ruptures of the Baza Fault.

Intersecting orthogonal strike‐slip faults with opposite senses of slip pose the question of what allows rupture to propagate through the junction and through both faults versus confining rupture to a single fault. I conduct dynamic rupture simulations on simplified orthogonal strike‐slip fault systems, to determine which conditions produce rupture on both component faults. In models with uniform initial tractions on both faults, slip on the first fault must reduce normal stress on the second fault for it to rupture. If the first fault ends at the cross fault, a stopping phase causes the cross fault to rupture. In models where I resolve a uniform regional stress field on the faults, only a narrow range of stress orientations allow multifault ruptures. These results will be helpful for evaluating hazard near orthogonal strike‐slip faults. Plain Language Summary: There are many examples around the world where two strike‐slip earthquake faults cross each other at nearly 90° angles. This is not remarkable when only one of the faults in a pair causes an earthquake, but it becomes notable when two or more crossing faults move at the same time. This raises the question of what causes the second fault to get involved, or not. To address this question, I use computer simulations of the physics of the earthquake process to test dozens of different fault configurations and earthquake starting points. I find that the location where the earthquake starts on the first fault controls whether the second fault is made stronger versus weaker, and therefore whether both faults can move together in one earthquake. These results can help us understand earthquake hazard around crossing faults. Key Points: Nucleation location effectively controls whether multifault rupture occurs on orthogonal strike‐slip fault systemsA stopping phase from rupture reaching the end of one fault is often required to initiate rupture on the cross faultOnly a narrow range of regional stress orientations allows both cross‐faults to rupture [ABSTRACT FROM AUTHOR]

Mechanical stratigraphy controls the growth patterns and dimensions of relatively small normal faults, yet how it influences the development of much larger structures remains unclear. Here, we use 3D seismic reflection data from the Outer Kwanza Basin, offshore Angola to constrain the geometry and kinematics of several normal faults formed in a deep‐water clastic succession. The faults are up to 6.3‐km long and 1.9‐km tall and have up to 44 m of throw. Aspect ratios and lower‐tip throw gradients are greater for faults that terminate downward at a c. 100 m thick, mass‐transport complex (MTC; up to 5.2 and 0.12) than for those that offset it (up to 2.7 and 0.01). Faults that offset the MTC invariably have >30 m of throw. Based on their geometric properties and throw patterns, we interpret that the faults nucleated above the MTC and propagated down toward it. Upon encountering this unit, which we infer behaved in a more ductile manner than encasing strata, tip propagation was halted until tip stresses were sufficiently high (corresponding to minimum throw of c. 30 m) to breach it. Faults with smaller throw were unable to breach the MTC. We argue that using only geometric criteria to determine fault growth patterns can mask the significant control mechanical stratigraphy has on fault kinematics. Mechanical stratigraphy is therefore a key control on the growth of large, seismic‐scale normal faults, in a similar way to that observed for far smaller structures. Key Points: Mechanical stratigraphy and fault interaction control the extent to which faults adhere to displacement–length scaling relationshipsFault growth and displacement–length relationships must be considered in context of local controls to be adequately describedDetailed kinematic analysis on faults developing within the same overall tectono‐stratigraphic setting can shed light on the kinematics of faults that do not reach the free surface [ABSTRACT FROM AUTHOR]

The 700‐km‐long Chaman fault (CF) marks the western edge of the plate boundary between India and Eurasia. Although global plate models predict 2.3–3.6 cm/yr left‐lateral motion between both plates, the fault is known to have hosted few earthquakes in historical times. Recent geodetic measurements attested the presence of aseismic slip locally. To detail the interplay between fast and slow slip along the CF, we build three Interferometric Synthetic‐Aperture Radar time series of ground deformation covering the whole fault length over 5 years (2014–2019). We find that most of the active fault trace slips aseismically and continuously. From south to north, we identify three creeping fault portions: the Nushki, Central, and Qalat segments of lengths between 80 and 130 km. The loading rate is 1.2 ± 0.3 cm/yr for the two southernmost portions, while it is about 0.7 ± 0.2 cm/yr for the Qalat segment. The Central segment and the nearby locked segments have hosted the largest known historical earthquakes on the CF, and three moderate magnitude earthquakes in our observation period. We image these earthquakes for which modeled slip at depth (Mw 5–5.6), time series of surface slip and deformation patterns argue toward large triggered aseismic slip. The June 2018 event displays postseismic moment 3–15 times greater than coseismic moment. Over the two decades covered by geodetic observations, continuous or triggered aseismic slip dominates along most of the fault and co‐locates with earthquakes. We observe that fault geometrical complexities delimit active segments and may be responsible for the kilometer‐scale intertwining between seismic and aseismic events. Plain Language Summary: The Chaman fault (CF) separates Indian and Eurasian tectonic plates moving at about 3 cm/yr with respect to each other. This fault is known to have hosted very few earthquakes (sudden and rapid slip) in historical times, and appears to slip slowly. We measure slip on the CF with radar satellite images covering 2014–2019, and assess the average amount of slip attributed to earthquakes over the past century. We identify three 80–130‐km‐long portions of the 700‐km‐long fault that slip silently, without radiating seismic waves, at rates reaching 1.2 ± 0.3 cm/year along the southern half, and 0.7 ± 0.2 cm/year in the north. The largest historical earthquakes and three earthquakes of moderate size in our observation period occurred on the central part of the CF. We image the surface deformation induced by these three earthquakes and find that they were followed by significant slow slip on the fault, especially the third event in June 2018. Since the 2000s space‐based measurements indicate that the slow mode of slip dominates on the CF plane and co‐locates with earthquakes. Changes of fault orientation may delimit portions of the fault which either slip slowly or break into fast earthquakes. Key Points: Aseismic slip along the Chaman fault (CF) reaches 12 mm/yr with three large distinct aseismic sections, despite marked elastic strain increaseThree earthquakes are imaged along the central CF and two exhibit large post‐seismic afterslipAlong‐strike distribution of slip is compared with historical seismicity and fault geometry [ABSTRACT FROM AUTHOR]

In the tectonic framework of the seismically active Tian Shan, the Jiashi‐Keping region, located at the mountain's southwestern front, is characterized by frequent moderate (M 5.5–7.0) earthquakes, rather than by strong (M ≥ 7.5) earthquakes as observed in other Tian Shan regions. What structures are responsible for these moderate earthquakes and whether or not more destructive earthquakes can occur remain largely unexplored. Based on an integrated analysis of published geologic and geomorphic mapping data, Google Earth satellite images, published and our collected 2D petroleum seismic‐reflection profiles, as well as available seismologic data and InSAR observations, in this study, we investigate active structural geometries and their correlation with the 2020 Mw 6.0, the 2003 Mw 6.2, and other recorded major earthquakes. Our study indicates that active structures in the region can be grouped into the Kepintagh thrust imbricates, the Bachu‐Kepintagh transpressional fault system, and the Xiasuhong transtensional fault system. The recorded major earthquakes are correlated with the Kepintagh thrust imbricates and the Xiasuhong fault system, both of which are highly segmented along strike and have limited downdip extensions. Although the recorded seismicity is characterized by moderate earthquakes, larger earthquakes may be generated by (i) the Bachu transpressional fault system (Mw ≥ 7.3) and the Kepintagh lower ramp (Mw > 7.5). Our study highlights the control of fault geometry on the seismic rupture process and regional seismic activities and enhances our understanding of seismic hazards in the Tian Shan. Key points: Regional active faults include Kepintagh thrusting, Xiasuhong transtensional, and probably Bachu‐Kepintagh transpressional fault systemsRecorded moderate earthquakes are sourced by along‐strike discontinuous and downdip‐extension‐limited Kepintagh and Xiasuhong fault systemsLarger earthquakes (Mw≥7.3) cannot be excluded for the Bachu fault system and Kepintagh lower ramp in terms of their structural geometries [ABSTRACT FROM AUTHOR]

The Xiadian orogenic deposit with ~ 100 t of gold resources, located in the Jiaodong Peninsula, Eastern China, shows an economically attractive gold mineralization that is spatially associated with the Zhaoping detachment fault. This study presents a combination of field investigation, 3D modeling, spatial analysis of models, and prospectivity modeling using a multi-layer perceptron deep neural network to understand the association between gold mineralization and structural deformation and to identify targets for deep exploration. The principal gold mineralization at Xiadian occurs as disseminated ore associated with strongly fractured rocks. Spatial analysis determined that gold mineralization is located primarily in convex segments (0–50 m) of the fault footwall within 150 m of the fault buffer where the dip is gentle (~ 35° to 55°). Sulfide vein ore with high Au grades is, by contrast, typically sited within steeply dipping fracture-fill distal to the Zhaoping detachment (150–250 m), consistently in areas with NE-strike and associated with steep segments of the detachment fault. The distribution of the orebodies is attributed to brittle–ductile deformation proximal to the detachment and distal brittle extension. Given the characteristics of syn-deformation mineralization and the stress regime (NW σ3, NE σ2, and vertical σ1), the NE-plunging ore shoots were interpreted as the result of auriferous fluid flow into oblique dilation zones formed by the coupled normal and dextral strike-slip movement of the Zhaoping fault. Six potential gold targets were identified sub-horizontal to, and NE-trending from, known orebodies. This study highlights the utility of 3D prospectivity modeling as a robust tool for understanding orogenic gold spatial distribution and related structural controls. [ABSTRACT FROM AUTHOR]

Abstract The 2018, Hokkaido Eastern Iburi, Japan, earthquake is an event characterized by complexity of the rupture process and slip pattern, which may involve both reverse and strike-slip motion depending on the locations on the fault surface. We perform dynamic rupture simulations based on simple physical laws, conditions for stressing and fault friction, and the non-planar fault geometry constrained by the aftershock observation. The complex fault geometry is numerically treated by the boundary integral equation method accelerated by the fast domain portioning method. The fault geometry is characterized primarily by the combination of six fault planes. As a result, we are able to explain several observed features of the event, including the spatial variation of the final fault slip and rupture velocity, which are inferred from the kinematic slip inversion. We also succeed in refining the constraint of the regional stress field in the focal area based on the simulation. Our results show that the overall patterns of the complex rupture event can be reproduced by a relatively simple model of the regional stress and the fault friction, if the geometrical complexity of the fault is properly taken into account.

The moderate-size 2010 Yushu earthquake (Mw 6.9) caused severe seismic damage around the source region which is situated on the sparsely populated hinterland of the Tibetan Plateau. But until now the mechanisms have not been well understood. To this end, we constructed the model to simulate the fault spontaneous rupture propagation, in which the realistic fault with a curved bend is imbedded. Our modeling results show that the special fault geometry controlled the rupture behavior, which encouraged rupture propagation speed transition from subshear to supershear with the speed of 5.17 km/s larger than the local shear wave velocity. Moreover, calculation results demonstrated that strong ground motion acceleration was greatly intensified by supershear ruptures, leading to widespread destruction. This may be the main reason why serious earthquake damage happened in the Yushu earthquake. In particular, we can see from numerous numerical experiments that rupture styles will be different if geometries of fault are varied. It is confirmed that it is the special geometry of the seismogenic fault of the Yushu mainshock produced grave seismic hazard. Thus, deeply investigating fault geometry will be helpful to better understand seismic source process and seismic hazard assessment.Key Points Curved bend in the seismogenic fault for the 2010 Yushu earthquake encouraged the supershear rupture transition. Strong ground motion acceleration was greatly intensified by supershear ruptures. Special geometry of the fault of the Yushu mainshock lead to seriouos seismic damage.

The 2021 Mw7.3 Maduo (Qinghai, China) earthquake ruptured ∼160 km along a poorly known secondary fault inside the Bayanhar block on the northern Tibetan plateau, which is generally parallel to the Kunlun fault. Here we integrate the interseismic (2015–2020) and coseismic geodetic observations to quantify the interseismic strain rate, fault geometry and coseismic slip distribution. Our results reveal that the seismogenic fault is featured by the low (<20 nanostrain/yr) and nearly undetectable interseismic strain rate. Three‐dimensional displacement fields and coseismic strain maps demonstrate the spatial variations of rupture kinematics due to the change of fault geometry. Our study reveals the steeply north‐dipping fault geometry of the seismogenic fault. The majority of coseismic slip occurred between 0 and 15 km with slight shallow slip deficit, which not penetrates through the inferred elastic upper crust (∼20–25 km). Our study highlights the tectonic and geometric control on fault kinematics of the Maduo earthquake. Plain Language Summary: The Bayanhar block on the north‐central Tibetan plateau experienced several major earthquakes with Mw > 6.5 on the boundary faults since the 1997 Mw7.6 Manyi earthquake, which is consistent with fast long‐term slip rates and large locking depth. However, the potential of generating large earthquakes on the faults inside the Bayanhar block is typically considered to be insignificant due to the estimated low fault slip rate (<5 mm/yr) and limited fault length. The 2021 Mw7.3 Maduo earthquake occurred on a poorly known subsidiary fault inside the Bayanhar block. Here, we use InSAR measurements prior to (2015–2020) and during the earthquake to investigate interseismic strain accumulation, to map coseismic deformation and to constrain the fault geometry as well as coseismic slip distribution of the Maduo earthquake. The three‐dimensional displacement fields and the relocated aftershocks support a slightly off vertical, dipping to the north, fault geometry. Our inversion results show that the coseismic slip is mainly distributed at a depth of 0–15 km. The kinematic analysis of the Maduo earthquake demonstrates that long‐term quantification of earthquake hazard of the subsidiary strike‐slip fault inside the Bayanhar block is challenging, because the interseismic strain rate is remarkably low and is largely undetectable by the geodetic observations. Key Points: Interseismic strain rate, coseismic three‐dimensional displacement field, finite strain and short‐term postseismic deformation of the Maduo earthquakeThe 2021 Maduo earthquake occurred on a subsidiary fault with a low (<20 nanostrain/yr) interseismic shear strain rateLong (∼160 km) rupture of the Maduo earthquake features shallow slip deficit with multi‐segment slip not penetrating through the upper crust [ABSTRACT FROM AUTHOR]

Eastern Bangladesh sits on the seismically active Chittagong‐Myanmar fold and thrust belt (CMFB), a north‐trending accretionary wedge on the eastern side of the India‐Eurasia collision. Earthquakes on the basal décollement and associated thrusts within the CMFB present a hazard to this densely populated region. In this study, we interpret 28 seismic reflection profiles from both published and unpublished sources to constrain the depth of the basal décollement. To convert profiles from the time domain to the depth domain, we integrate sonic log and seismic stacking velocity data to generate time‐velocity relationships for different parts of the CMFB. Our analysis reveals that the décollement is ∼9 km deep in northeast and southeast Bangladesh, but shallows to ∼5 km in east‐central Bangladesh. The décollement has an area of 7.25 × 104 km2 (∼150 × 450 km), making it capable of an Mw 8.5 earthquake. However, the warped geometry of this fault might act as a rupture barrier were a large earthquake to occur on the décollement. Our combined velocity and fault model lay the groundwork for future studies to address seismic segmentation, ground shaking, and rupture modeling in the CMFB. Finally, we use our compiled data set to analyze the evolution of fold kinematics in the CMFB. We observe that folding style and failure mode varies, from mainly ductile deformation in the foreland to mainly brittle in the hinterland. The dual‐failure modes within the CMFB support the hypothesis that a region with ductile deformation may still be capable of seismic behavior. Plain Language Summary: In Bangladesh, the interface between two tectonic plates has created a large, nearly flat, earthquake‐producing fault called a décollement. Above the décollement, sedimentary layers are compressed to create a series of north‐south trending folds and thrust faults that extend for hundreds of kilometers. This study uses 28 seismic reflection data sets (a technology similar to ultrasounds), originally collected by the oil and gas industry, and reinterprets the subsurface structure of the folds as a way to determine the décollement depth. We find that the décollement has a curved shape, with a depth of ∼9 km in northeast and southeast Bangladesh, and ∼5 km in east‐central Bangladesh. We hypothesize that the greater amount of sediment deposited in north and south Bangladesh has weighed down the surface of the earth and warped the décollement. This is the first study to constrain the geometry of this décollement, and we find that it has the potential to host a magnitude 8.5+ earthquake. The fault model presented here can be incorporated into studies of the rupture patterns of earthquakes on this décollement to better understand the earthquake hazard in a region inhabited by over 160 million people. Key Points: We present a data‐based seismic velocity model for Bangladesh to 10 km depthThe décollement underlying the Chittagong‐Myanmar fold and thrust belt forms a broad arch at depths of 5–9 kmThis morphology is likely caused by lithospheric flexure due to sediment deposition and may impact rupture propagation [ABSTRACT FROM AUTHOR]

Abstract Interferometric synthetic aperture radar and teleseismic P‐wave data were combined to investigate the source rupture characteristics of the 2020 Mw 6.3 Yutian, China, earthquake. We first utilized the near‐field displacements together with 29 broadband teleseismic P waveforms to investigate the uniform slip model by using a Bayesian bootstrap optimization nonlinear inversion method to resolve the nucleation point location, origin time and fault geometrical parameters. Based on these results, the kinematic rupture process of the earthquake was inverted. We conclude that the 2020 Yutian earthquake occurred on a west‐dipping blind normal fault with a dip of ∼60° at the junction of Bayan Har block (BHB) and western Kunlun block (WKB). The rupture nucleation point was located at longitude = 82.435°E, latitude = 35.619°N, and a depth of 6.377 km, making it shallower than estimated from other point source inversions. The earthquake started at 21:05:20 UTC on June 25, 2020 and was delayed by ∼2.07 s compared with GeoForschungsZentrum. The rupture lasted for ∼12 s, with a total seismic moment of ∼3.28 × 1018 Nm, corresponding to a moment magnitude of 6.3. The slip was mainly confined between ∼3.0 and 9.0 km in depth and the peak slip was ∼1.40 m, which occurred at a depth of ∼6.377 km. The slip was predominantly normal slip with slight right‐lateral strike‐slip components, which agrees with the southeastward movement of the BHB relative to the WKB.

The subsurface fault geometry is the base for understanding a process of crust deformation and mountain building. Based on kinematic models for fault-related folds, a geomorphic method is recently applied to estimate the subsurface fault geometry, while the validation on its reliability is lacking. In this study, we surveyed a suit of river terrace surfaces across an active fold at the north front of the Qilian Shan. According to the deformation geometry of the terraces, the fold deformation is interpreted by a listric fault fold model, and based on this kinematic model, the fault geometry underlying the fold is estimated. In comparison between the estimated fault geometry and a seismic reflection profile, we found that the decollement depth and the back thrust are highly consistent with each other. Although some small fault bends or internal shearing cannot be estimated solely by the terrace deformation, the overall fault geometry is successfully revealed by the terrace deformation. Using this fault geometry and the terrace dating results, the region deformation kinematics are re-evaluated, which suggest that the dip slip (in a rate of 1.8 ± 0.4 mm/a) along the decollement is mainly accommodated by two structures, one is the blind-back-thrust fault within the piggy basin in a dip-slip rate of 0.9 ± 0.3 mm/a and another is the thrust and fold at the west portion of the Yumu Shan range.

Interferometric synthetic aperture radar and teleseismic P‐wave data were combined to investigate the source rupture characteristics of the 2020 Mw 6.3 Yutian, China, earthquake. We first utilized the near‐field displacements together with 29 broadband teleseismic P waveforms to investigate the uniform slip model by using a Bayesian bootstrap optimization nonlinear inversion method to resolve the nucleation point location, origin time and fault geometrical parameters. Based on these results, the kinematic rupture process of the earthquake was inverted. We conclude that the 2020 Yutian earthquake occurred on a west‐dipping blind normal fault with a dip of ∼60° at the junction of Bayan Har block (BHB) and western Kunlun block (WKB). The rupture nucleation point was located at longitude = 82.435°E, latitude = 35.619°N, and a depth of 6.377 km, making it shallower than estimated from other point source inversions. The earthquake started at 21:05:20 UTC on June 25, 2020 and was delayed by ∼2.07 s compared with GeoForschungsZentrum. The rupture lasted for ∼12 s, with a total seismic moment of ∼3.28 × 1018 Nm, corresponding to a moment magnitude of 6.3. The slip was mainly confined between ∼3.0 and 9.0 km in depth and the peak slip was ∼1.40 m, which occurred at a depth of ∼6.377 km. The slip was predominantly normal slip with slight right‐lateral strike‐slip components, which agrees with the southeastward movement of the BHB relative to the WKB. Plain Language Summary: The June 25, 2020 Mw 6.3 Yutian earthquake occurred on a normal fault within the northwestern Tibetan Plateau, a region characterized by EW extension. Because the rough focal mechanism solutions and source locations based on seismic wave data reported by several institutes are different, they cannot be used to create a reliable fault model for inverting the kinematic rupture process of this moderate earthquake. The initial rupture point and origin time, which are the prior information needed for revealing the slip characteristics, are also unclear. To resolve these rupture parameters, both Sentinel‐1 interferometric synthetic aperture radar and teleseismic P‐wave data were collected. We found that the rupture mainly occurred in the shallow part of a ∼60° west‐dipping normal fault, the nucleation point was near the point of maximum slip, and the slip lasted ∼12 s and extended northward and southward. Our research has important significance for understanding earthquake disasters and deformation mechanisms of faults within the northwestern Tibetan Plateau. Key Points: Both interferometric synthetic aperture radar and teleseismic P wave are used to resolving the ambiguity of fault geometryThe rupture nucleation point and origin time were clarified by a Bayesian bootstrap optimization inversion methodRupture process of this earthquake is predominantly normal‐sip with slight right‐lateral strike‐slip components [ABSTRACT FROM AUTHOR]

The geometries (i.e., dip angles) of active faults from the surface to the seismogenic zone are the most important factors used to evaluate earthquake ground motion, which is crucial for seismic hazard assessments in urban areas. In Osaka, a metropolitan city in Japan, there are several active faults (e.g., the Uemachi and Ikoma faults), which are inferred from the topography, the attitude of active faults in surface trenches, the seismic reflection profile at shallow depths (less than 2 km), and the three-dimensional distribution of the Quaternary sedimentary layers. The Uemachi and Ikoma faults are N–S-striking fault systems with total lengths of 42 km and 38 km, respectively, with the former being located ~ 12 km west of the latter; however, the geometries of each of the active faults within the seismogenic zone are not clear. In this study, to examine the geometries of the Uemachi and Ikoma faults from the surface to the seismogenic zone, we analyze the development of the geological structures of sedimentary layers based on numerical simulations of a two-dimensional visco-elasto-plastic body under a horizontal compressive stress field, including preexisting high-strained weak zones (i.e., faults) and surface sedimentation processes, and evaluate the relationship between the observed geological structures of the Quaternary sediments (i.e., the Osaka Group) in the Osaka Plain and the model results. As a result, we propose geometries of the Uemachi and Ikoma faults from the surface to the seismogenic zone. When the friction coefficient of the faults is ~ 0.5, the dip angles of the Uemachi and Ikoma faults near the surface are ~ 30°–40° and the Uemachi fault has a downward convex curve at the bottom of the seismogenic zone, but does not converge to the Ikoma fault. Based on the analysis in this study, the dip angle of the Uemachi fault zone is estimated to be approximately 30°–40°, which is lower than that estimated in the previous studies. If the active fault has a low angle, the width of the fault plane is long, and thus the estimated seismic moment will be large. [ABSTRACT FROM AUTHOR]