Seismogenic Zone Research Articles

Enigmatic Surface Ruptures at Cape Rytyi and Surroundings, Baikal Rift, Siberia: Seismic Hazard Implication

The geomorphic expression of active faulting and distinction of paleoseismic events in areas that are rapidly obscured by erosion/sedimentation still remains a considerable scientific problem. The present article discusses the revealing of surface faulting ruptures and their parameters to identify capable faults without trenching and to estimate the magnitude of earthquakes. The case study was at Cape Rytyi, located in Baikal-Lena Nature Reserve on the northwestern shore of Lake Baikal. Based on unmanned aerial photography, GPR, and structural observations, we mapped and investigated the relation between geomorphological forms and ruptures. The obtained results show that past landslides and paleoruptures at Cape Rytyi and its surroundings are associated with at least two earthquakes. The Mw of the earlier event was 7.3 (Ms = 7.4); the Mw of the later one was 7.1 (Ms = 7.3). The paleoruptures in the distal part of the delta of the Rita River and on the southeastern slope of the Baikal Ridge were included in the seismogenic rupture zone, which traces some 37 km along the Kocherikovsky fault. The approximate intervals in which earthquakes occurred are 12–5 ka and 4–0.3 ka, respectively. The applied analysis methods can be useful for paleoseismology and assessing seismic hazards in similar regions elsewhere.

Analysis of Seismicity in Elbasani Geothermal Zone during Five Last Decades

During the last five decades since 1968, the Elbasani zone in Albania is dominated by low seismicity. This fault zone has evidenced some seismological phenomena. It is a transversal active fault zone with normal and strike-slip faults characterized by low-velocity layers. From instrumental seismicity in this area have been located 1831 seismic events, with magnitude ML > 1.1 Richter and 3 of them with magnitude ML > 5.0. Recently on May 19, 2014, a moderate earthquake ML5.2 occurred in the Elbasani seismogenic zone. This earthquake was located 6 km southeast of Belshi town and 18 km south of Elbasani town and was felt over Albania. The analysis of seismicity during the period of time 968-2022 years indicates that most earthquakes were located in the upper and middle earth crust and very few in the uppermost mantle. The focal mechanisms solution of some earthquakes indicates an extensional regime in this fault zone. The parameters of focal mechanisms indicate the predominance of normal fault motion and strike-slip fault motion is compatible with the present-day extension regime. The Elbasan zone has a complex geological structure with a significant number of geothermal water resources and some features of seismicity. The analysis of the features of seismicity during the five last decades provides us with correct interpretations of seismotectonic, low velocity lowers and understanding of geodynamic phenomena of the Earth's crust. The area of active faults of Elbasan presents a threat to the city of Elbasan and to the curative and tourist area of Llinxha-Hidrat.

Open Source Data-Based Solutions for Identifying Patterns of Urban Earthquake Systemic Vulnerability in High-Seismicity Areas

Urban settlements located in high-seismicity areas should benefit from comprehensive vulnerability analyses, which are essential for the proper implementation of vulnerability modelling actions. Alas, many developing countries face a shortage of knowledge on seismic vulnerability, particularly concerning its systemic component, as a consequence of a combination of data scarcity and a lack of interest from authorities. This paper aims to identify primary time-independent spatial patterns of earthquake systemic vulnerability based on the accessibility of key emergency management facilities (e.g., medical units, fire stations), focusing on the urban settlements located in the high-seismicity area nearby the Vrancea Seismogenic Zone in Romania. The proposed methodological framework relies on open source data extracted from OpenStreetMap, which are processed via GIS techniques and tools (i.e., Network Analyst, Weighted Overlay Analysis), to compute the service areas of emergency management centres, and to map earthquake systemic vulnerability levels. The analysis shows that accessibility and systemic vulnerability patterns are significantly impacted by a synergy of factors deeply rooted in the urban spatial layout. Although the overall accessibility was estimated to be medium-high, and the overall systemic vulnerability to be low-medium, higher systemic vulnerability levels in certain cities (e.g., Bacău, Onești, Tecuci, Urziceni). The presented findings have multi-scalar utility: they aid in the development of improved, locally tailored seismic vulnerability reduction plans, as well as the allocation of financial and human resources required to manage earthquake-induced crises at regional scale. Further to that, the paper provides a transparent methodological framework that can be replicated to put cities in high-seismicity areas on the map of systemic vulnerability assessments, laying the groundwork for positive change in countries where the challenges associated with high-level seismic risk are often overlooked.

Apparent Low‐Velocity Belt in the Shallow Anninghe Fault Zone in SW China and Its Implications for Seismotectonics and Earthquake Hazard Assessment

AbstractThe Anninghe fault forms the eastern boundary of the Sichuan‐Yunnan block in Southwest China and has been identified as an earthquake gap zone. This study intends to construct the upper crustal shear wave velocity (Vs) structure beneath the Anninghe fault to understand its seismotectonics and potential large earthquake hazards. We deployed a dense seismic array along the southern central Anninghe fault valley. From the 3‐month continuous records, we calculated vertical‐component cross‐correlation functions (CCFs). However, the surface wave signals in the CCFs are intensely interfered by near zero‐time‐lag noise. We proposed a mode separation method based on the high‐resolution linear Radon transform, which suppressed the interfered noise and greatly enhanced the surface wave signals for ambient noise tomography of the Vs structure. The fine upper crustal structure reveals a distinct narrow low‐velocity belt within a depth of 3 km beneath the Anninghe fault zone. At deeper depths (4.5–8 km), the narrow low‐velocity belt shifts to the east and correlates with the distribution of local earthquakes. Combining previous results with our new findings, we presented a seismotectonic model of the southern central Anninghe fault, which interprets the narrow low‐velocity belt as a water‐contained fracture zone that forms a seismogenic zone at deeper depths under transpression. In addition, we demonstrated through scenario earthquake simulations that fine structures play a significant role in the assessment of earthquake hazards along the Anninghe fault. As such, this study provides a typical window into seismotectonics and large earthquake hazards in the active southeastern Tibetan Plateau.

The 2013 Slab‐Wide Kamchatka Earthquake Sequence

AbstractStudies of initiation of large earthquakes are usually focused on frictional instabilities occurring in the near vicinity of the future rupture. Possible contributions of long‐distance interactions with large‐scale tectonic instabilities remain unknown. Here we analyze seismic catalogs and geodetic time series during a few months preceding the 2013 M = 8.3 deep‐focus Okhotsk earthquake. This deep‐focus event is preceded by four intense seismic clusters in the seismogenic zone. GNSS time series in Kamchatka revealed a transient landward motion episode 1 month prior to the mainshock, consistent with an increase of seismogenic zone loading. This transient loading episode is accompanied by a doubling of the intermediate depth seismicity rate suggesting a transient slab pull as the origin. These observations question the constant subducting velocity hypotheses and may have implications in the understanding of the long‐distance along‐slab stress interactions and in their contribution to initiation of large deep‐focus earthquakes.

The 2013–2020 seismic activity at Sabancaya Volcano (Peru): Long lasting unrest and eruption

Sabancaya volcano is the youngest and second most active volcano in Peru. It is part of the Ampato-Sabancaya volcanic complex which sits to the south of the ancient Hualca Hualca volcano and several frequently active faults, thus resulting in complex volcano-tectonic interactions. After 15 years of repose, in 2013, a series of 4 earthquakes with magnitude >4.5 occurred within 24 h, marking the beginning of a new episode of unrest. Several additional swarms of earthquakes occurred in the following years until magmatic eruptive activity started on 6 November 2016. This activity is ongoing as of this writing, with an average of 50 explosions per day. In this study, we present results of multiparametric monitoring of Sabancaya's activity observed during 2013–2020. Seismic data are used to create a one-dimensional seismic velocity model, to catalog, locate, and characterize earthquakes, to detect repeating earthquake families, and to monitor seismic velocity variations by ambient noise cross-correlation. These analyses are complemented by visual and remote sensing observations and ground deformation measurements. All monitored parameters showed significant changes on 6 November 2016, the day of eruption onset, thus dividing the eruptive activity into pre-eruptive and eruptive stages.The unrest is characterized by high levels of seismic activity with hundreds of events detected per day. Volcano-tectonic (VT) earthquakes were dominant during the pre-eruptive period while long-period (LP) events and explosions have been most numerous since the eruption onset. Earthquake locations highlight long-lasting seismogenic zones along multiple previously active regional faults, as well as along newly identified faults. This VT seismicity is mainly distributed in a sector from the northwest to the east of the volcanic complex at distances of up to 30 km from the crater. We focus our analysis on two eruptive episodes: the eruption onset and subsequent crater migration from south to north, and the increase of lava dome extrusion rate in 2019. Both episodes are accompanied by seismic velocity decreases of up to 0.2% and are preceded by a few weeks by bursts of distal VT activity, including numerous repeating earthquakes. These repeated events were located on several remote tectonic faults (5–25 km from the vent). We suggest that these phenomena could be due to the injection of a batch of magma in the deep reservoir and/or conduit, which would generate 1) a pressure wave propagating in the hydrothermal system, triggering the bursts of seismic activity and 2) slow rising of magma by melting old material filling the conduit that eventually produced the eruptive and dome growth acceleration events.

A Microphysical Model of Rock Friction and the Brittle‐Ductile Transition Controlled by Dislocation Glide and Backstress Evolution

AbstractRate‐ and state‐friction (RSF) is an empirical framework that describes the complex velocity‐, time‐, and slip‐dependent phenomena observed during frictional sliding of rocks and gouge in the laboratory. Despite its widespread use in earthquake nucleation and recurrence models, our understanding of RSF, particularly its time‐ and/or slip‐dependence, is still largely empirical, limiting our confidence in extrapolating laboratory behavior to the seismogenic zone. While many microphysical models have been proposed over the past few decades, none have explicitly incorporated the effects of strain hardening, anelasticity, or transient elastoplastic rheology. Here we present a new model of rock friction that incorporates these phenomena directly from the microphysical behavior of lattice dislocations. This model of rock friction exhibits the same logarithmic dependence on sliding velocity (strain rate) as RSF and displays a dependence on the internal backstress caused by long‐range interactions among geometrically necessary dislocations (GNDs). Changes in the backstress (internal stress) evolve exponentially with plastic strain of asperities and are dependent on both the current backstress and previous deformation, which give rise to phenomena consistent with interpretations of the “critical slip distance,” “memory effect,” and “evolution effect” of RSF. The rate dependence of friction in this model is primarily controlled by the evolution of backstress and temperature. We provide several analytical predictions for RSF‐like behavior and the “brittle‐ductile” transition based on microphysical mechanisms and measurable parameters such as the GND density and strain‐dependent hardening modulus.

Strain accumulation and release associated with the occurrence of SSEs in the subduction zones of the Japanese Islands

We investigated the spatiotemporal changes in strain associated with the occurrence of slow slip events (SSEs) in the subduction zones of the Japanese Islands and compared the spatial distribution of both the amount of strain accumulated for the period before and during the SSEs release using time series data from the Global Navigation Satellite System (GNSS). In this study, four SSEs were analysed: the Tokai long-term SSE (2000–2005), the Boso-Oki short-term SSE (2007), and the Bungo Channel long-term SSEs (2009–2011 and 2018–2019). As a result, we found strong negative correlations for all four dilatations before and during SSE occurrence. For these dilatations, we estimated the amount of strain released at the time of occurrence of the SSE relative to that accumulated during the period prior to the SSE. The result indicates that not all the accumulated strain before the SSEs was released when the SSEs occurred. Moreover, it is highly likely that the strain released by SSE is not only due to the strain accumulation just below the SSE occurrence region, but also due to the strain accumulation on the shallower plate boundary, which is a seismogenic zone for a future megathrust earthquake.

In-situ physical and elastic properties of Archaean basement granitoids in the Koyna seismogenic zone, western India from 3 km downhole geophysical well logs: Implications for water percolation at depth

The Koyna region in western India is characterized by recurrent artificial water reservoir triggered seismicity since the impoundment of the Shivaji Sagar lake in 1962. Physical and mechanical properties of the crystalline basement rock beneath the Koyna seismogenic zone are determined from the downhole geophysical logs in a 3 km research borehole. The borehole, KFD1, is located in the proximity of the Donichawadi fissure zone (DFZ), the surface manifestation of the fault associated with the 1967 M 6.3 Koyna earthquake. KFD1 penetrated 1247 m thick Deccan traps and continued 1767 m into the underlying granite-gneiss basement. Geophysical logs provide a rare opportunity to characterize the basement granitoids, which host the recurrent seismicity. Several major zones of anomalous physical and mechanical properties were delineated below a depth of 2 km from the analysis of downhole geophysical logs. High Poisson's ratio and high neutron porosity across anomalous zones imply the presence of water saturated rock. The analyses of Nuclear Magnetic Resonance log reveals the presences of clay bound and free fluid across these zones. Sharp decrease in dynamic Young's modulus with increasing Poisson's ratio imply the presence of fractures saturated with water. These hydraulically conductive fractures are potential locations for fault reactivation. Additionally, strike of major and minor fractures associated with the anomalous zones are compared with the strike of Donichawadi fault. High angle between strike of the major shear fractures in the anomalous zone and strike of the Donichawadi fault imply that borehole KFD1 may have either intersected subsidiaries of the main fault at multiple depths or a new hidden fault system.

Enhanced 3D velocity structure, seismicity relocation and basement characterization of Changning shale gas and salt mining regions in southern Sichuan basin

Current extensive seismicity in southern Sichuan Basin is ascribed to the reactivation of pre-existing faults, as a result of prolonged fluid injection for salt mining and shale gas development, respectively. However, the structural framework of the region remains poorly understood. Here, we apply Vp/Vsconsistency-constrained double-difference seismic tomography to high quality phase data from 36,314 earthquakes jointly recorded by our local array and a regional seismic network to determine high-resolution velocity models. Earthquake relocations reveal shallow hypocenters for the Ms>5.0 earthquakes and two distinct seismogenic zones corresponding to the salt mine and shale gas regions, with most induced seismic events forming widespread lineaments some of which extend to the basement and are remarkably similar to the fault and fracture trends interpreted on reflection seismic and outcrops, respectively. Our 3-D crustal velocity analyses show that seismicity beneath the Changing salt mining area is associated with a combination of relatively low Vp/Vs(1.6–1.74) and high Vp/Vs(1.75–1.86) expressions, while most of small earthquakes within the Xingwen shale gas block are associated with relatively high Vp/Vsvalues (1.77–1.87), indicating the earthquakes in these two areas are caused by unique inducing mechanisms. The two moderately strong 2018 Xingwen Ms5.7 and 2019 Gongxian Ms5.3 earthquakes in the Xingwen shale gas block are located around low Vp/Vs. zones, suggesting they could be structurally controlled. In comparison, the 2019 Changning Ms6.0 earthquake in the Changning salt mining area is associated with high Vp/Vs. expression, suggesting its occurrence is related to fluid injections. In addition, top of the crystalline early Neoproterozoic (pre-Sinian) Sichuan basement is characterized by the 6.5 km/s Vpcontour, which is new for earthquake tomographic studies in the region. Combined with outcrop analysis, we are able to construct a structural framework for induced seismicity in southern Sichuan basin, which unravels the structural architecture of induced seismicity.

Integrated structural, geophysical and remote sensing data for characterizing extensional linked fault systems and related land deformation hazards at Cairo-Suez District, Egypt

A tiered framework has been developed to map and characterize extensional linked fault systems at Cairo-Suez District, Egypt to better understand seismic and land deformation hazards. As a test site, the present tiered approach is applied on the Cairo-Suez District (CSD), Egypt, which is considered as an active seismogenic zone that encompasses complex deformation settings. Tier 1, which involves outcrop observations and structural mapping, reveals four geometries belonging to soft- and hard-linkage transfer zones (TZs). Considering the mapped TZs, the near-surface heterogeneity is characterized using a 2D-electrical resistivity tomography (2D-ERT) technique at “tier 2” applying an advanced inversion scheme (AIS). Inversion results of the 2D-ERT profiles indicate a subsurface continuity of the mapped surface geological and structural features including shale smears of fault zones and shale layer distributions. At “tier 2”, we particularly identify a highly deformed shale bed of 15-m thick, that accommodate fault slips as indicated by internal gouge zones. We further represent all local hazard factors including slope, drainage density, land use/land cover and displacement maps using the “Fuzzy-Shannon Entropy Model”. Overall, our findings show that most of the delineated deformation hazards in the test site are related to hard-linked TZs. At “tier 3”, we ultimately provide our recommendation for regional and site-specific engineering mitigation and to avoid active faults by defining setback-zone buffer (i.e., no build zone) correlated with the local fault activity and near-surface geometry. For Egypt, an economically developing country, this approach represents a major advancement from conventional geologic and geotechnical investigations through introducing a tiered framework for risk-relevant characterization of the extensional linked fault systems to apply adequate procedures/precautions by engineers for urban planning.

Active Faults Revealed and New Constraints on Their Seismogenic Depth from a High-Resolution Regional Focal Mechanism Catalog in Myanmar (2016–2021)

ABSTRACTWe derive a new earthquake focal mechanism catalog for 86 Mw>4.0 earthquakes that occurred in the Myanmar region from 2016 to 2021. We apply the generalized Cut-and-Paste inversion method to a new set of regional broadband waveform data to obtain the earthquake focal mechanism and centroid depth with uncertainties estimated in a bootstrapping manner. Compared with global earthquake catalogs, our results are better aligned with mapped, active faults and reveal seismic activity along unmapped, blind faults. Our new catalog shows that the Sagaing Fault is more active in its northern segment with deeper seismogenic zone (∼27 km) compared to its southern segment that has a shallower seismogenic zone (∼10 km), sandwiching a seismic gap in its central segment. Earthquakes that occurred on the unmapped, blind faults beneath the Central Myanmar Basin at shallow depths (3–12 km) suggest a dominating northeast–southwest compressional stress field. Shallow earthquakes beneath the Indo-Myanmar Range (IMR) are rare, instead, north–south-oriented strike-slip faults are active within the deep accretionary wedge or lower crust of the Myanmar plate between depths of 20 and 40 km. At the eastern edge of the IMR, earthquakes with high-angle thrust mechanisms occurred between depths of 30 and 48 km, likely along steep faults separating the accretionary wedge from the Myanmar forearc crust. High-resolution intraslab focal mechanisms show that to the north of 22° N, slab deformation is dominated by strike-slip earthquakes with subvertical fault planes down to a depth of ∼25 km beneath the slab, suggesting lateral shear within the slab due to the northward motion of the Indian plate. To the south, more normal-faulting earthquakes suggest a stronger role of plate-bending processes in the slab deformation.

A Discrete Elements Study of the Frictional Behavior of Fault Gouges.

A series of discrete elements simulations is presented for the study of fault gouges' frictional response. The gouge is considered to have previously undergone ultra-cataclastic flow and long-time consolidation loading. We explore the effect of different particle characteristics such as size, polydispersity, and also shearing velocities on gouge's response under the conditions met in the seismogenic zone. Monte-Carlo analyses suggest that the local stick-slip events disappear when averaging over a large number of numerical samples. Moreover, the apparent material frictional response remains almost unaffected by the spatial randomness of particles' position and by the particle's size distribution. On the contrary, the mean particle size controls the formation and thickness of the observed shear bands, which appear after the peak friction is met. Furthermore, the apparent friction evolution fits well to an exponential decay law with slip, which involves a particle size dependent critical slip distance. For the studied conditions and depth, the shearing velocity is found to play a secondary role on the apparent frictional response of the gouge, which highlights the importance of analyses involving multiphysics for studying the rheology of fault gouges. Besides improving the understanding of the underlying physics of the problem, the above findings are also useful for deriving pertinent constitutive models in the case of modeling with continuum theories.

Non-stationary problem of high earth dam under seismic impact

The design, construction, and operation of high earth dams on the territory of seismogenic zones of Uzbekistan set the task for researchers to improve the calculation methods and algorithms for various kinds of loads to which the structure is subjected. At present, the application of the dynamic method (a wave theory) to the calculation of an earth structure, i.e., the dynamic stress-strain state of high earth dams under seismic impact within the framework of wave theory, is the most difficult problem in the mechanics of a deformable rigid body. In the example of the operating Charvak high earth dam, the numerical method of finite differences solved the problem of studying the stress-strain state with shear dynamic impact on the base (such as the records of a real seismogram of an earthquake). At the boundaries of the considered finite region of the base, the so-called radiation conditions are set, i.e., diffraction is not considered. The result of solving the non-stationary problem is presented as isolines of displacements and stresses along the dam's body, depending on time. At that, the most vulnerable zones of the considered earth dam were identified.

Hydrothermal Friction Experiments on Simulated Basaltic Fault Gouge and Implications for Megathrust Earthquakes

AbstractNucleation of earthquake slip at the plate boundary fault (décollement) in subduction zones has been widely linked to the frictional properties of subducting sedimentary facies. However, recent seismological and geological observations suggest that the décollement develops in the subducting oceanic crust in the depth range of the seismogenic zone, at least in some cases. To understand the frictional properties of oceanic crustal material and their influence on seismogenesis, we performed hydrothermal friction experiments on simulated fault gouges of altered basalt, at temperatures of 100–550°C. The friction coefficient (μ) lies around 0.6 at most temperature conditions but a low μ down to 0.3 was observed at the highest temperature and lowest velocity condition. The velocity dependence of μ, (a−b), changes with increasing temperature from positive to negative at ∼100°C and from negative to positive at ∼450°C. Compared to gouges derived from sedimentary facies, the altered basalt gouge showed potentially unstable velocity weakening over a wider temperature range. Microstructural observations and microphysical interpretation infer that competition between dilatant granular flow and viscous compaction through pressure‐solution creep of albite contributed to the observed transition in (a−b). Alteration of oceanic crust during subduction produces fine grains of albite and chlorite through interactions with interstitial water, leading to reduction in its frictional strength and an increase in its seismogenic potential. Therefore, shear deformation possibly localizes within the altered oceanic crust leading to a larger potential for the nucleation of a megathrust earthquake in the depth range of the seismogenic zone.

Experimental Simulation of Fluidity State Evolution from Slow Earthquake to Strong Earthquake

Slow earthquakes have been frequently detected in the neighboring regions of seismogenic zones. However, the physics underlying the connection of slow earthquakes to huge earthquakes remains one of the least understood aspects in seismic research today. By observing the strain time evolution of granular materials in laboratory experiments, we report that the fluidity parameter, which represents a time dependent rate for strain evolution, can be rejuvenated under progressively increasing dynamic loading and this physical quantity finally reaches a critical threshold state, and beyond that catastrophic and uncontained dynamic strain response begins. Moreover, memory effect and anti-fluidization effect in granular materials are also observed, and these effects underlay the long-term slip events preceded some strong earthquakes, including 2011 Mw 9.1 Tohoku earthquake and 2014 Mw 8.1 Iquique earthquake. The observations indicate the key role of development of fluidity state on connecting slow earthquakes and to full-scale earthquakes.

B-value variations in the Central Chile seismic gap assessed by a Bayesian transdimensional approach

The b-value can be used to characterize the seismic activity for a given earthquake catalog and provide information on the stress level accumulated at active faults. Here we develop an algorithm to objectively estimate variations of b-value along one arbitrary dimension. To this end, we employ a Bayesian transdimensional approach where the seismic domains will be self-defined according to information in the seismic catalog. This makes it unnecessary to prescribe the location and extent of domains, as it is commonly done. We first show the algorithm’s robustness by performing regressions from synthetic catalogs, recovering the target models with great accuracy. We also apply the algorithm to a microseismicity catalog for the Central Chile region. This segment is considered a seismic gap where the last major earthquake with shallow slip was in 1730. Our results illuminate the downdip limit of the seismogenic zone and the transition to intraslab seismicity. In the along-strike direction, low b-value coincides with the extent of locked asperities, suggesting a high-stress loading at the Central Chile seismic gap. Our results indicate the reliability of the Bayesian transdimensional method for capturing robust b-value variations, allowing us to characterize the mechanical behavior on the plate interface of subduction zones.

Depth‐Varying Friction on a Ramp‐Flat Fault Illuminated by ∼3‐Year InSAR Observations Following the 2017 Mw 7.3 Sarpol‐e Zahab Earthquake

AbstractWe use interferometric synthetic aperture radar observations to investigate the fault geometry and afterslip evolution within 3 years after a mainshock. The postseismic observations favor a ramp‐flat structure in which the flat angle should be lower than 10°. The postseismic deformation is dominated by afterslip, while the viscoelastic response is negligible. A multisegment, stress‐driven afterslip model (hereafter called the SA‐2 model) with depth‐varying frictional properties better explains the spatiotemporal evolution of the postseismic deformation than a two‐segment, stress‐driven afterslip model (hereafter called the SA‐1 model). Although the SA‐2 model does not improve the misfit significantly, this multisegment fault with depth‐varying friction is more physically plausible given the depth‐varying mechanical stratigraphy in the region. Compared to the kinematic afterslip model, the mechanical afterslip models with friction variation tend to underestimate early postseismic deformation to the west, which may indicate more complex fault friction than we expected. Both the kinematic and stress‐driven models can resolve downdip afterslip, although it could be affected by data noise and model resolution. The transition depth of the sedimentary cover basement interface inferred by afterslip models is ∼12 km in the seismogenic zone, which coincides with the regional stratigraphic profile. Because the coseismic rupture propagated along a basement‐involved fault while the postseismic slip may activate the frontal structures and/or shallower detachments in the sedimentary cover, the 2017 Sarpol‐e Zahab earthquake may have acted as a typical event that contributed to both thick‐ and thin‐skinned shortening of the Zagros in both seismic and aseismic ways.

Laboratory Study on Frictional Behavior of Rock Blocks of Meter Scale. Methods and Preliminary Results

—The paper introduces a new, unique for Russia, meter-scale laboratory setup created in the Institute of Geosphere Dynamics of the Russian Academy of Sciences (IDG RAS) to study the development of different sliding regimes on rock discontinuities. The experimental procedure is described and the results of the first series of tests aimed at studying the formation of different sliding regimes on rock faults are presented. The laboratory fault was a loaded contact of two 75-cm long blocks made of diabase. The fault was filled with granular material (a fault gouge). Normal stresses on the fault can reach 10 MPa. By varying fault gouge composition and loading rate, we reproduced a wide range of sliding regimes: sliding with constant velocity, regular stick-slip, and aperiodic slow slip episodes. It is shown that a variation in the loading rate can cause a significant change in the sliding regime. Intense frictional crushing of gouge grains is detected in the experiments at relatively low normal pressure of 2 MPa. In the case of high-amplitude stick-slip, besides crushing of the gouge material, also structural phase transformations of quartz grains corresponding to a local temperature increase up to 700°C are revealed. A possible set of the problems related to deformation processes in seismogenic fault zones that can be addressed by modeling on such setups—the pre-seismic stage of inelastic behavior of the main fault zone at critical stress—is outlined.

Geometric Control on Seismic Rupture and Earthquake Sequence Along the Yingxiu‐Beichuan Fault With Implications for the 2008 Wenchuan Earthquake

AbstractThe 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.