Nearby Faults Research Articles

Dependence of horizontal PGA prediction on site conditions at critically short hypocentral distances: an analysis based on Japanese strong motion data

A fundamental element in seismic hazard assessment studies is the prediction of maximum values of ground accelerations. This is achieved through an attenuation relationship, also known as the ground motion prediction equation (GMPE). Seismic design of structures generally assumes the occurrence of a large earthquake on a nearby fault. The predictions from the GMPE are linked to the magnitude, distance, and local site conditions. Critical near‑fault ground motions have unique characteristics such as magnitude saturation and distance saturation. To investigate the impact of local site conditions in the near‑fault region, we analyzed the attenuation relations using the accelerations recorded in Japan by the Kiban Kyoshin network (KiK‑net) from 2000 to 2023. All events with a magnitude equal to or greater than 5.5, and with focal depth and epicentral distance less than 20 km, were considered. In this study, the hypocentral distance ranges from 5 to 23 km. Ordinary least squares regression analysis was conducted for the maximum horizontal vector of ground accelerations. The attenuation relation coefficients for horizontal vector peak ground acceleration were estimated using an effective distance and saturation effect. An additional independent term, the average shear wave velocity of the soil to a depth of 30 meters below the ground surface, and more cost‑effective shallower depths were included in the regression model. The analysis revealed that the site condition is statistically non‑significant at very short hypocentral distances of less than 23 kilometers.

Interplay of tectonic and dynamic processes shaping multilayer extensional system in southern-central Apennines

High-hazard seismic zones can remain silent over centuries with meager seismicity rates challenging our understanding of seismic processes. We focus on the comprehensive analysis of cascading episodes of swarms and seismic sequences following the 2009 L’Aquila mainshock (MW 6.3) in the southern-central Apennine that previously experienced ~ M7 earthquakes. We enhance the seismic catalog, unmasking low-magnitude seismicity down to completeness magnitude ML ~ 0, and we unveil that the microseismicity might be secondarily triggered by the L’Aquila mainshock, influencing the frictional properties in the nearby fault zones or opening fault valves generating the intense seismic activity detected from 2009 to 2013. The diffusivity, observed in the most seismic episodes, and the high Vp/Vs values (> 1.88) indicate fluid circulation promoting multilayered extensional seismicity within 11–15 km and 16–23 km depth ranges. Mapping the 3D distribution of seismicity alongside geological data reveals an evident tectonic influence, unveiling unknown geometric aspects and providing the first evidence of a NNE-dipping deformation zone bounding at depths of 11–15 km the overlying fault system. Deeper seismicity suggests a mantellic CO2 ascending shape. These findings enrich the literature on tectonic seismic swarms in extensional domains, providing essential constraints on fluid involvement in the seismic processes and contributing to forthcoming discussions on the seismotectonic setting in high-seismic-risk areas of the Apennines of Italy.

Triggered and recurrent slow slip in North Sulawesi, Indonesia

Nearby faults interact with each other through the exchange of stress. However, the extent of fault interaction is poorly understood. In particular, interactions may lead to slow-slip activity, resulting in episodes of transient surface motion. Our study concentrates on Northwest Sulawesi (Indonesia), which hosts two fault zones with potential for major earthquakes and tsunamis: the strike-slip Palu-Koro fault and the Minahassa subduction zone. Thanks to a 20-year-long effort in geodetic monitoring, we are able to identify multiple periods during which surface velocities deviate from their interseismic trend. We use a Bayesian methodology with forward predictions of slip on the two fault interfaces to match the observations following the 2018 Mw7.5 Palu earthquake, and infer that both deep afterslip on the Palu-Koro fault and slow slip on the Minahassa subduction interface have caused the observed transient surface motion. This finding represents the first recording of a slow slip event on the Minahassa subduction interface. We also infer that the subduction interface and the strike-slip fault are likely interacting on a regular basis.

Analysis of changes in shale mechanical properties and fault instability activation caused by drilling fluid invasion into formations

During the drilling process, the issue of drilling fluid loss can lead to changes in the mechanical properties of the formation, thereby altering the stress environment of nearby faults. In order to assess the risk of fault activation during drilling operations, the Ordos M area shale was selected as the research object. Mechanical experiments were conducted on rock samples immersed in water-based drilling fluid with a pressure differential of 2 MPa and a temperature of 50 °C. The changes in the mechanical properties of the shale before and after immersion in drilling fluid were determined. Based on the experimental results, combined with the spring combination model and fault activation theory, a quantitative evaluation of fault activation risk was conducted. The findings revealed that the shale in this region has a high clay content, demonstrating a certain level of water sensitivity. The presence of micro-pores and micro-fractures is well-developed, increasing the interaction probability between drilling fluids and clay minerals. After immersion in drilling fluid, there was a varied decline in all mechanical strengths of the shale. The elastic modulus is positively correlated with the shear strength and Coulomb stress of the fault plane. The Poisson’s ratio is positively correlated with the shear strength and negatively correlated with the Coulomb stress. The greater the internal friction and cohesion, the higher the shear strength of the fault plane, and the larger the friction coefficient, the smaller the Coulomb stress, resulting in a more stable fault.

Analysis of impulsive ground motions from February 2023 Kahramanmaraş earthquake sequence

AbstractOn the 6th of February 2023, a large magnitude earthquake (Pazarcık earthquake), $$M_{w}=$$ M w = 7.7, occurred in southeast Türkiye, which caused significant destruction in Türkiye and Syria. Relatively large magnitude aftershocks followed the main shock, and after 9 hours of the main event, another large magnitude earthquake (Elbistan earthquake) occurred, $$M_{w}=$$ M w = 7.6, on a nearby fault. This study analyzes the near-fault seismic signals from earthquakes larger than 5.5 recorded between the main shock and the 31st of March 2023. More than 60 impulsive motions are detected in 3 earthquakes, mostly concentrated in the Pazarcık and Elbistan earthquakes. In the Pazarcık earthquake, many impulsive motions are recorded in near-fault stations with periods of up to 14 s. In contrast, in the Elbistan earthquake, impulsive motions are spatially distributed, with pulse periods of up to 11 s and at distances greater than 150 km. Pulse periods mostly correlate with the magnitude of the earthquake, but pulse probability models do not predict impulsive motions over long distances. The presence of strong impulsive motions in vertical components is also observed. For both earthquakes, peak ground velocities (PGVs) are larger than predicted by ground motion prediction equations. The observation of long-period, large amplitude signals may indicate the presence of a directivity effect for both earthquakes. In some stations, spectral periods exceed the 2018 Turkish building design codes for long periods ($$\ge$$ ≥ 1 s).

Characteristic Analysis of the MS6.8 Luding Earthquake Sequence in Sichuan, China, on Sept. 5, 2022

The MS6.8 earthquake occurred in Luding County, Ganzi Prefecture, Sichuan Province, China at 12:52 p.m. on Sept. 5, 2022, is preliminarily analyzed in terms of regional tectonics, historical earthquakes, sequence characteristics and focal mechanism solutions. The results show that: 1) The magnitude difference ΔM between the main shock and the maximum aftershock, the energy ratio ER of the main shock to the aftershocks, the p-value and the estimated Mmax-value indicate that the MS6.8 Luding earthquake sequence can be divided into the main shock-aftershock type (MAT). 2) The spatial distribution of the aftershocks, their corresponding b-value and expected maximum magnitude are obviously segmented, which reflect the complexity of the seismogenic tectonics. 3) According to the focal mechanisms obtained from the HASH program, the geometry distribution of the sequence, and the relationship between the sequence and the nearby faults, it can be inferred that the Moxi fault segment of the Xianshuihe fault zone is the seismogenic tectonics of the MS6.8 Luding earthquake sequence.

Prediction of fractures by using the stress field in medium to deep shale gas reservoirs: A case study of Late Ordovician–early Silurian shale gas reservoirs in the Nanchuan region, South China

Containing representative shale gas reservoirs, the Wufeng–Longmaxi Formations in the Nanchuan region in South China have been modified by three compressional tectonic movements (in the Late Jurassic to Early Cretaceous, Late Cretaceous, and Cenozoic to present). These movements formed deeply buried reservoirs with complex in situ stresses and rapid directional changes, influencing the state of shale gas reservoirs fractures. Numerical modeling applied to model the stress field in the three periods showed that the maximum compressive horizontal principal stress (σHmax), strain rate, and fault slip rate regimes transformed from a thrust and strike-slip regime to a thrust and strike-slip regime to a thrust regime. The σHmax direction changed from NW–SE and WNW–ESE to WNW–ESE and nearly E–W and then to multiple directions. The magnitude of strain rate changed from -16 ∼ -15 to -17.7 ∼ -15.6 to -18.8 ∼ -15.8. The fault slip rate changed from 0 to 0.029 mm/a to 0–0.0026 mm/a to 0–0.0012 mm/a. These data strongly agree with data from drilled wells and outcrops and show that the present-day stress field is the superposed effect of multistage movements. Semiquantitative research was conducted on shale gas fracture openness and fracture development. Fracture openness is favorable when the angle between the σHmax direction and nearby faults is ≤45° and the distance from the fault is >2 km, and fracture development is favorable when the magnitude of the strain rate is <-17.6 and the fault slip rate is <0.00036 mm/a. Accordingly, we predicted fractured shale gas reservoirs.

Are You Ready to Shake, Dunedin?

In 1974, I was one and a half years old. We didn’t own a TV in my family’s home, but if we had, we would have perhaps been watching Coronation Street on the evening of the 9th of April (there was only one channel - black and white). My mother remembers the quake quite vividly - the sound like a train roaring outside confused her for a moment. Once she realised what was happening, Mum rushed her little family under the strong, oaken kitchen table as our little bungalow shook.&#x0D; I have been making a short film in which we find out about Dunedin’s moderate intensity 1974 earthquake through the firsthand memories of residents of the time and through our local earthquake experts. How did the earthquake manifest in our city? What damage was caused? How prepared were we?&#x0D; Taking this experience further: what is the potential for a larger quake to happen here - could nearby fault lines deliver an earthquake to Dunedin of catastrophic proportions? Have we learned from 1974 enough to be prepared for such an event?In my presentation, I’ll discuss my process for finding a narrative approach to this documentary in order to bring about audience engagement and even entertainment to this communication resource, thus amplifying retention of the important messaging of resilience and preparedness.&#x0D; Supervised by: Caroline Orchiston&#x0D; Scholarship Project Funded by: The Centre for Sustainability

On the role of poroelastic stressing and pore pressure diffusion in discrete fracture and fault system in triggering post-injection seismicity in enhanced geothermal systems

Injection-induced seismicity has become one of the most critical challenges for the widespread deployment of Enhanced Geothermal Systems (EGS). In particular, some EGS development projects have led to large, damaging earthquakes that unexpectedly occurred far off the stimulated reservoir region and, in particular, after stopping fluid injection. Yet, the causative mechanisms of these seismicity patterns remain highly elusive. Here, we identify a combination of mechanisms that could explain delayed seismicity in EGS sites by conducting fully-coupled hydromechanical simulations of the hydraulic stimulation of a naturally-fractured granitic reservoir. The model comprises a sparse network of long, variably-oriented fractures interacting with a nearby, critically-oriented fault. The results show that the presence of fractures introduces notable nonlinearities in the flow field and rock deformation and significantly expands the rock volume affected by fluid injection. First, the stimulated fracture network provides highly-permeable conduits for communicating elevated pore pressure over long distances. Second, the anisotropic expansion of fractures generates shear stress that is transmitted almost instantaneously across the reservoir. The pore pressure and stress perturbations can not only cause slip along fractures, inducing (micro)seismicity during injection, but also affect the stability of nearby faults, which may not necessarily be pressurized during injection. The transferred poroelastic stresses can increase or decrease the slip tendency along different fault segments. However, the fault may reactivate only after several months following injection when a progressive pore pressure diffusion modulated by the transient fault permeability evolution brings a critically-stressed fault segment to failure conditions. We also find that the spatiotemporal evolution of seismicity depends largely on the nearby fault orientation, hydromechanical properties, and hydraulic connection with the fracture network, as well as the initial state of stress. We conclude that accurate subsurface characterization and continuous monitoring during and after injection should allow for managing the risks posed by injection-induced seismicity and safely unlocking the immense potential for clean and sustainable geothermal energy.

Nearby fault interaction within the double-vergence suture in eastern Taiwan during the 2022 Chihshang earthquake sequence

Nearby faults interact through stress changes induced by fault slip and viscoelastic flow. The process is, however, often elusive and can be geometry-dependent and time-variant. Here, we combine geodetic and field observations to characterize the interaction of two head-to-head, conjugate faults in eastern Taiwan during the 2022 Chihshang earthquake sequence. We map the coseismic slip on the Central Range fault and dynamically-triggered shallow slip on the Longitudinal Valley fault, which has been creeping interseismically. Overlapping of seismic and aseismic slip suggests that the Longitudinal Valley fault is capable of hosting a variety of distinct slip behaviors. Moreover, substantial slip on the Central Range fault suppresses Coulomb stress on the Longitudinal Valley fault, and vice versa, resulting in seismic bursts in an out-of-phase pattern on the two faults as seen in the hundred-year historical records. Such fault interaction implies the need for time-dependent seismic hazard reassessment for the complex fault system.

Evidence of a Transient Aseismic Slip Driving the 2017 Valparaiso Earthquake Sequence, From Foreshocks to Aftershocks

AbstractFollowing laboratory experiments and friction theory, slow slip events and seismicity rate accelerations observed before mainshocks are sometimes interpreted as evidence of a nucleation phase. However, such precursory observations still remain scarce and are associated with different time and length scales, raising doubts about their actual preparatory nature. We study the 2017 Valparaiso Mw = 6.9 earthquake, which was preceded by aseismic slip accompanied by an intense seismicity, suspected to reflect its nucleation phase. We complement previous observations, which have focused only on precursory activity, with a continuous investigation of seismic and aseismic processes from the foreshock sequence to the post‐mainshock phase. By building a high‐resolution earthquake catalog and searching for anomalous seismicity rate increases compared to aftershock triggering models, we highlight an over‐productive seismicity starting within the foreshock sequence and persisting several days after the mainshock. Using repeating earthquakes and high‐rate GPS observations, we highlight a transient aseismic perturbation starting 1‐day before the first foreshock and continuing after the mainshock. The estimated slip rate over time is lightly impacted by large magnitude earthquakes and does not accelerate toward the mainshock. Therefore, the unusual seismic and aseismic activity observed during the 2017 Valparaiso sequence might be interpreted as the result of a slow slip event starting before the mainshock and continuing beyond it. Rather than pointing to a possible nucleation phase of the 2017 Valparaiso mainshock, the identified slow slip event acts as an aseismic loading of nearby faults, increasing the seismic activity, and thus the likelihood of a large rupture.

Holocene Earthquake Cycles of an Active Tectonic Block Boundary Fault Zone: A Case Study in the Qilian–Haiyuan Fault Zone, Northeastern Tibet Plateau

Abstract Fault zones along active tectonic block boundaries are a significant source of devastating continental earthquakes. Strong earthquakes produce disruptions of sediment and induce characteristic sediments near the fault, which serve as valuable sedimentary evidence for identifying and dating of paleoearthquakes. In this study, we aimed to reconstruct the earthquake history of the Qilian–Haiyuan fault zone in the northeastern Tibetan Plateau during the Holocene. We reanalyzed forty-four trenches and used the sedimentary sequences, event indicators, and age constraints to determine the earthquake history. Our analysis revealed the paleoearthquakes of 6 subsidiary faults of the Qilian–Haiyuan fault zone with accurate event ages and rupture extents. Based on the spatial and temporal distributions of strong earthquakes since 10 ka, we identified five earthquake clusters around the central-eastern Qilian–Haiyuan fault zone including seven rupture cascades where the earthquakes migrated gradually from east to west. The existing seismic gap reveals that the latest migration may not yet be complete and suggests a high probability of M ≥ 7 earthquakes occurring on the Jinqianghe fault, Maomaoshan fault, and the central part of the Lenglongling faults. We concluded that, in order to better understand earthquake cycles and seismic hazards, it is important to consider a fault zone as a whole, including multiple faults and their interaction on the earthquake triggering between nearby faults.

Characteristics and controlling factors of natural fractures in deep lacustrine shale oil reservoirs of the Permian Fengcheng Formation in the Mahu Sag, Junggar Basin, China

Maximum burial depth of the Permian Fengcheng Formation, which contains significant oil resources in the Mahu Sag, exceeds 4500 m constituting a deep reservoir. Fractures are widespread in the deep lacustrine shale, providing important storage space and fluid conduits. Based on cores, borehole image logs, and petrographic analysis, we identified and analyzed the types and characteristics of fractures and conducted X-ray diffraction (XRD) and Cathodoluminescence (CL) tests to investigate the factors influencing fracture development. We found faults and three types of opening-mode fractures: bed-bounded fractures, unbounded fractures and bed-parallel fractures. Mechanical stratigraphy (depositional rock types modified by diagenesis) and proximity to faults, are main controls of opening-mode fracture occurrence and attributes. Specifically, Dolomitic mudstone with a high carbonate content shows the highest fracture density of 5.24 m-1. Bed-bounded fractures are mainly developed within brittle layers, with fracture density lower in thicker beds. On the other hand, unbounded (tall) fractures primarily occur near reverse faults (0–1500 m), with orientations parallel to extension directions implied by the kinematics of nearby faults. Additionally, the effectiveness of early-formed fractures as fluid conduits is likely modified (reduced) by cementation and enhanced by dissolution. Fractures have a positive influence on oil production, corresponding to high fracture density and aperture in the vertical direction.

Study on Fault Activation Evaluation Mechanism and Sensitivity Analysis in Shale Gas Hydraulic Fracturing

Shale gas hydraulic fracturing usually activates nearby faults and makes them slip. In horizontal wells, fault slip can result in serious casing deformation. Casing deformation slows the fracturing process, lowers production, and raises the cost of a shale gas well. It is challenging to obtain underground data on fault activation because of the deep shale reservoirs. As a result, the current study needed to indicate how hydraulic fracturing affects fault activation length. This made it challenging to control casing deformation. The fluid-structure coupled finite element method was used in this study to create a coupled seepage-stress model for heterogeneous shale formation. With microseismic signs and hydraulic fracking for shale gas, this model examined the variation law of pore pressure and ground stress. The fault activation coefficient was created to assess the fault activation duration and the impact of hydraulic fracturing. The model was verified by the microseismic signal. The outcomes of the numerical simulation demonstrate how the rapid rise in formation pore pressure during hydraulic fracturing affected the ground stress at the fault interface. The influence of ground stress variation at the fault interface on fault activation could not be ignored. Increased fault elastic modulus, fracture pressure, fracture time, and the fault Poisson ratio result in longer fault activation lengths. The length of the fault’s activation was decreased by the increase in fracture stage, distance from the fault, friction angle within the fault, and fault angle. Finally, a shale gas horizontal well with casing deformation in block C was analyzed. The results showed that reducing the fracturing duration can reduce the activation length of the fault by 68.25%, resulting in a 9.1 mm fault slide and a 0.86 mm casing deformation, respectively. This study offers theoretical guidelines for preventing fault activation during hydraulic fracturing in horizontal shale gas wells.

Long silence on the East Anatolian Fault Zone (Southern Turkey) ends with devastating double earthquakes (6 February 2023) over a seismic gap: implications for the seismic potential in the Eastern Mediterranean region

On 6 February 2023, an earthquake with magnitude M w c . 7.0 on the Narlı Fault, a fault sub-parallel to the East Anatolian Fault Zone (EAFZ), initiated a chain of large earthquakes on the EAFZ. The earthquakes occurred in a seismic gap with low geodetic strain rates and low background seismicity, where deformation is distributed across a wide fault zone and a long recurrence time of historical earthquakes. The c. 50 km long rupture of the Narlı Fault towards Pazarcık led to an M w 7.8 left-lateral strike-slip earthquake breaking a c. 300 km section of the c. 600 km long EAFZ bilaterally with a total duration of more than 80 s. Toward the SW, the rupture propagated on the c. 100 km long Amanos segment with a peak surface offset of 5 m, before diminishing toward the Hatay graben. In the NE direction, the rupture reached a peak surface offset of 7 m before sharply declining at the termination of the 2020 M w 6.8 Sivrice earthquake rupture. A second large earthquake with M w 7.6 occurred 9 h later on the Çardak Fault, located at the western margin of (and sub-parallel to) the EAFZ and breaking the surface with almost 9 m of left-lateral slip (average of c. 4 m). Following these large earthquakes, the increase in the regional stress led to a rapid seismic activation in a broad region from central to eastern Anatolia, loading the faults at various scales and increasing seismic hazard. Two weeks after the initiation of the seismic crisis, a third earthquake with M w 6.4 occurred at the southern boundary of the Hatay graben, near the southwestern termination of the Amanos rupture. The earthquakes caused significant loss of human life, devastating 12 cities. We evaluate the observations prior to the ruptures, and present preliminary seismological results with surface displacements from sub-pixel correlation of optical satellite images and the stress perturbations computed on the nearby faults based on preliminary slip models. The re-evaluation of the seismic potential in light of the recent and historical earthquakes provides some new insight on seismic hazard assessment. The recent series of events on the EAFZ is an important reminder that large faults can generate very large earthquakes on multiple segments. The seismic potential of large earthquakes on these fault zones can be estimated only by considering multiple seismic cycles and moment deficits from very large earthquakes. Supplementary material : Supplementary figures and tables are available at https://doi.org/10.6084/m9.figshare.c.6567094

Rock avalanche caused by nearby fault creep: A case study of Muztag rock avalanche

Rock avalanche caused by nearby fault creep: A case study of Muztag rock avalanche

Rock avalanche caused by nearby fault creep: A case study of Muztag rock avalanche

Rock avalanche caused by nearby fault creep: A case study of Muztag rock avalanche

Evaluation of rockburst risk in deep tunnels considering structural planes based on energy dissipation rate criterion and numerical simulation

Risk evaluation of rockbursts is crucial for ensuring the safety and stability of deep-tunnel construction. Structural planes play an important role in the potential for rockbursts to occur, and more related studies are needed. In this paper, a new rockburst criterion based on energy dissipation rate (EDR) index was developed and the role of structural planes in rockburst potential was studied. Firstly, a method of generating structural planes and a viscoplastic constitutive model with internal variables were introduced for numerical simulation using FLAC3D. Then, tunnelling under different hydrostatic pressures was simulated to obtain the values of EDR index and stress-strength/strength-stress ratio indices, and EDR criterion was estimated based on stress-strength/strength-stress ratio criteria. The results show that EDR criterion is: 0∼1 Pa/s (no risk), 1∼3 Pa/s (low risk), 3∼6 Pa/s (medium risk), >6 Pa/s (high risk). Subsequently, the case study of rockbursts in the drainage tunnel of Jinping II hydropower station was conducted considering a nearby fault to validate EDR criterion. The simulation results using EDR are basically consistent with actual rockbursts in terms of the potential location, intensity and time. Finally, the numerical simulations of tunnel excavation considering a joint set with different dip angles were implemented, and the influences of joints on rockbursts were investigated using EDR. Analysis results show that rockbursts are prone to occur on tunnel walls parallel and perpendicular to joints; potential rockburst locations are controlled by local joints and are less affected by geostress distribution. Joints may increase the intensity, delay the occurrence time, and control the pit boundaries of rockbursts. Obtained EDR contours can be used to infer damaged zones of strainbursts and strain-structure slip rockbursts and understand inducing mechanism of rockbursts in Jinping II drainage tunnel. The findings of this study provide insight into the role of structural planes on rockbursts. The presented methodology can be applied to rockburst evaluation and prevention during the construction of deep tunnels.

P-wave velocity structure beneath reservoirs and surrounding areas in the lower Jinsha River

The lower reaches of the Jinsha River are rich in hydropower resources because of the high mountains, deep valleys, and swift currents in this area. This region also features complex tectonic structures and frequent earthquakes. After the impoundment of the reservoirs, seismic activity increased significantly. Therefore, it is necessary to study the P-wave velocity structure and earthquake locations in the lower reaches of the Jinsha River and surrounds, thus providing seismological support for subsequent earthquake prevention and disaster reduction work in reservoir areas. In this study, we selected the data of 7,670 seismic events recorded by the seismic networks in Sichuan, Yunnan, and Chongqing and the temporary seismic arrays deployed nearby. We then applied the double-difference tomography method to this data, to obtain the P-wave velocity structure and earthquake locations in the lower reaches of the Jinsha River and surrounds. The results showed that the Jinsha River basin has a complex lateral P-wave velocity structure. Seismic events are mainly distributed in the transition zones between high- and low-velocity anomalies, and seismic events are particularly intense in the Xiluodu and Baihetan reservoir areas. Vertical cross-sections through the Xiangjiaba and Xiluodu reservoir areas revealed an apparent high-velocity anomaly at approximately 6 km depth; this high-velocity anomaly plays a role in stress accumulation, with few earthquakes distributed inside the high-velocity body. After the impoundment of the Baihetan reservoir, the number of earthquakes in the reservoir area increased significantly. The seismic events in the reservoir area north of 27° N were related to the enhanced activity of nearby faults after impoundment; the earthquakes in the reservoir area south of 27° N were probably induced by additional loads (or regional stress changes), and the multiple microseismic events may have been caused by rock rupture near the main faults under high pore pressure.

Seismic Activity in the Celje Basin (Slovenia) in Roman Times—Archaeoseismological Evidence from Celeia

Searching for unknown earthquakes in Slovenia in the first millennium, we performed archaeoseismological analysis of Roman settlements. The Mesto pod mestom museum in Celje exhibits a paved Roman road, which suffered severe deformation. Built on fine gravel and sand from the Savinja River, the road displays a bulge and trench, pop-up structures, and pavement slabs tilted up to 40°. The city wall was built over the deformed road in Late Roman times, supported by a foundation containing recycled material (spolia) from public buildings, including an emperor’s statue. We hypothesize that a severe earthquake hit the town before 350 AD, causing widespread destruction. Seismic-induced liquefaction caused differential subsidence, deforming the road. One of the nearby faults from the strike-slip Periadriatic fault system was the seismic source of this event.