Fault Offsets Research Articles

A Novel Scalable and Reliable Control for DC Microgrids With Varying Number of Agents.

Existing scalable control methods mainly rely on a fixed block-diagonal structure for the Lyapunov matrix, potentially resulting in numerical infeasibility issues. To overcome this limitation, this article proposes a novel scalable and reliable control scheme for dc microgrids. Initially, a general model for dc microgrids is established to enhance reliability, considering scenarios involving loss of control effectiveness (LoCE) and offset faults. Subsequently, a structured free-weight matrix technique is introduced to mitigate negative coupling effects of power lines, and to address numerical infeasibility by avoiding the assumption about the Lyapunov matrix. Furthermore, the stability of the entire dc microgrid is guaranteed by checking local agent conditions, independently of power line couplings. Therefore, the proposed control scheme ensures plug-and-play scalability with varying number of agents. Finally, theoretical results are validated through numerical simulations using the MATLAB/SimPowerSystems toolbox.

Structural inheritance, morphotectonics and Holocene activity of the Cucalón-Pancrudo extensional fault (Iberian Chain, Spain)

A large active extensional fault is characterized, which will contribute to improve seismic hazard assessment in an intraplate region of the Mediterranean domain whose seismic potential until today has been underestimated. The Cucalón-Pancrudo fault (CPF) is a NNW-SSE striking extensional fault that represents the negative (extensional) inversion of a previous, Late Variscan and Paleogene, transpressional fault. It is part of the Río Grío-Pancrudo fault zone, the largest active structure in the intraplate central-eastern Iberian Chain. This study focuses on characterizing the CPF structure, morphotectonics and paleoseismology. The fault offsets a well-known regional planation surface (lower level of the Fundamental Erosion Surface, FES3, 3.5 Ma), providing a basis for calculating the maximum fault throw (c. 280–300 m, including the contribution of a minor synthetic fault and a gentle accommodation monocline), and hence estimating the net slip (c. 305–325 m) and the long-term slip rate (0.09 mm/a). A trench dug in a Holocene fluvial terrace reveals an ensemble of listric, domino-style ruptures synthetic to the main fault. They were activated during at least two paleoseismic events (X, Z), 14.9 ± 1.4 ka and 6.9 ± 0.4 ka in age. An intermediate event (Y), dated to 11.0 ± 1.0 ka, as an uncertainty because the attribution of their ruptures to event Z cannot be ruled out. The total net slip accumulated during these events (1.15–1.25 m) provides an apparent short-term slip rate of 0.07–0.09 mm/a, very similar to the long-term slip rate. Nevertheless, it is very likely that the faults exposed in the surveyed trench do not represent a complete paleoseismic record. Therefore, the actual slip rate during Late Pleistocene-Holocene times could be significantly higher, approaching that reported for most of the recent extensional faults in the central Iberian Chain. The inclusion of CPF in the map of Active Faults of Iberia (QAFI database) will result in improving seismic hazard assessment in Spain.

The performance of transmission pipelines on February 6th, 2023 Kahramanmaras earthquake: a series of case studies

AbstractThe Mw 7.8 earthquake that struck Kahramanmaras on February 6, 2023, caused significant damage to hydrocarbon and water transmission lines due to permanent ground deformations. Explosions occurred in the Kahramanmaras–Gaziantep Natural Gas transmission pipeline causing disruption of gas in the cities of Gaziantep and Hatay Fault Displacement Hazard. Water transmission pipelines were also damaged due to extreme fault offsets. This paper presents the findings from the earthquake affected zones. Particular attention was given to the performance of the 2600 mm diameter Duzbag–Gaziantep water transmission pipeline due to three reported challenges: (#1) The behavior of the pipe in the Duzbag tunnel, (#2) The pipe bend behavior at the fault crossing, and (#3) The pull-out of the pipe at valve room. Based on field observations, extreme damage forms at pipes were presented. A three-dimensional nonlinear numerical analysis was performed to simulate the observed damage of the reoriented water pipe at the fault crossing for Case #2. The nonlinear finite element model was able to capture the complex nature of the post-yield behavior of steel pipe as well as the damage locations along the pipe and their sequence of occurrences. As a mitigation measure for new pipes, the need for using flexible connections at critical crossings was highlighted. For existing pipes, the importance of developing secondary (by-pass) lines was emphasized.

A dynamic rifting model of the Caroline Ridge, West Pacific

The dynamic rifting model of the Caroline Ridge, an oceanic plateau in the West Pacific, remains unclear. Previous studies have revealed that the crustal width of the Caroline Ridge clearly varies from the northwest to the southeast. Here, we investigate Caroline Ridge rifting using numerical simulation and reveal the relationships between different plateau sizes and rifting. For small-scale oceanic plateaus, the initial rift migration forms a symmetric structure. The strain transfers toward one side of the margin. The rifting manifests as a process from landward migration to seaward migration, accompanied by a transition from a symmetric to an asymmetric structure. With larger initial lithospheric thinning, larger ranges of crustal deformation appear with less basement tectonic subsidence. For large-scale oceanic plateaus, basement tectonic subsidence occurs within narrower ranges. The proximal half-graben structures rapidly abandon with relatively horizontal basement and small fault offsets on the rift flank, leading to a symmetric structure during breakup. With larger initial lithospheric thinning, crustal deformation occurs within a narrower region during the early stage. Landward rift migration dominates the deformation, resulting in an asymmetric structure. On the basis of these results, for the large-scale section of the Caroline Ridge, NW-trending normal faults originated from the landward rift migration. At present, seaward rift migration dominates the Caroline Ridge deformation, forming a relatively symmetric structure. For the small-scale section of the Caroline Ridge, the initial landward rift migration occurred. Larger initial lithospheric thinning may exist beneath the Caroline Ridge during the initiation of Miocene rifting.

Land-use interactions, Oil-Field infrastructure, and natural processes control hydrocarbon and arsenic concentrations in groundwater, Poso Creek Oil Field, California, USA

Like many hydrocarbon production areas in the U.S., the Poso Creek Oil Field in California includes and is adjacent to other land uses (agricultural and other developed lands) that affect the hydrology and geochemistry of the aquifer overlying and adjacent to oil development. We hypothesize that the distributions of hydrocarbons and arsenic in groundwater in such areas will be controlled by complex interactions between mixed land uses, oil-field infrastructure, and natural processes. In 2020–2021, samples of groundwater and surface water were collected and analyzed for a large suite of inorganic and organic chemicals and isotope and gas tracers to test this hypothesis. Those data are supplemented with ancillary data on historical geochemistry, hydrology, geology, and oil-field infrastructure. Hydrocarbons in groundwater (e.g., methane through pentane gases and benzene) are associated with natural processes (e.g., fault offsets or transition in sediment depositional environment) and oil-field infrastructure (e.g., fluid-migration pathways associated with uncemented annulus in oil wells or unlined pits). Arsenic concentrations >10 μg per liter (μg/L; maximum concentration 12.9 μg/L) are associated with natural processes in old, high-pH groundwater, and more recent recharge of water from natural and/or engineered recharge processes. Along the southwest margin of the oil field, pumping for drinking-water and irrigation supplies in combination with engineered groundwater recharge produce a depression in groundwater elevations where groundwater with elevated sulfate concentrations from agricultural areas and groundwater with hydrocarbons from the oil field mix to produce a zone of sulfate reduction that removes hydrocarbons and arsenic from groundwater but produces elevated sulfide (S2-) concentrations (maximum concentration 29 mg per liter, mg/L). In this study, multiple approaches were required to resolve the overlapping effects of land uses, oil-field infrastructure, and natural processes on the distributions of hydrocarbons and arsenic in groundwater. The combined use of geographic, historical, physical, chemical, isotopic, and other information to constrain processes could be a useful approach for studies in other hydrocarbon-production areas. This is particularly important where land uses affect aquifer hydrology to an extent that causes mixing of groundwaters with different chemical compositions.

High-resolution co-seismic fault offsets of the 2023 Türkiye earthquake ruptures using satellite imagery

On February 6, 2023, southern Türkiye was struck by two large earthquakes at 01:17 UTC (Mw=7.8, Pazarcık, Kahramanmaraş) and 10:30 UTC (Mw = 7.6, Elbistan, Kahramanmaraş), causing severe damage at the complex junction of the Dead Sea Fault (DSF), the Cyprus Arc and the East Anatolian Fault Zone (EAFZ). The ruptures propagated along several known strands of the southwestern termination of the EAFZ, the main Pazarcık and Karasu valley faults, and the Çardak-Sürgü fault. Here we present the high-resolution mapping of the entire coseismic surface rupture and an estimate of the rupture width, total and on-fault offset, and diffuse deformation obtained a few days to three months after the two mainshocks. The mapping is derived from image correlation of Sentinel-2 optical satellite imagery and validated with offset measurements collected on the ground. We find that the ruptures extend over lengths of 310 km and 140 km for the Mw 7.8 and Mw 7.6 mainshocks, respectively. The maximum offsets reach 7.5 ± 0.8 m and 8.7 ± 0.8 m near the epicenters of the Mw 7.8 and Mw 7.6 events, respectively. We propose a segmentation of the two ruptures based on these observations, and further discuss the location of the potential supershear rupture. The use of optical image correlation, complemented by field investigations along earthquake faults, provides new insights into seismic hazard assessment.

Transient ground bulge derived from a dense broadband seismic array during an aquifer step-drawdown pumping test

We deployed nine broadband seismometers within 32 m of a well where a four-hour step-drawdown pumping test was conducted at a site near Shuangxi, Taiwan in August 2018. The whole seismic array simultaneously recorded hour-long seismic signals that we interpret as ground tilt evolution during and after the pumping. Here we document data from the six seismic stations within 10 m of the boreholes to show the ground deformation closest to the pump site. Instead of ground subsidence due to groundwater over-draw as documented in many long-term and even short-term GPS and tiltmeter studies, four out of the six seismometers show ground bulge at the pumping site with radial directional tilts up to 1.0 × 10-5 rad for a drawdown of 10 m. We interpret the bulge as the dominant transient crustal response to the unloading of discharged groundwater, outpacing the minor subsidence due to reduced effective stress from reduced pore fluid head. After the pumping stopped, the ground tilt signals started to level back, though not simultaneously: some signals recorded delays up to an hour, possibly related to replenishment of groundwater from far-field or different aquifers and shallow depth hydrological heterogeneity including minor faults and fractures. We conducted seismic reflection surveys and found mostly layered structures without large fault offsets. Using an analytical solution that calculates the ground tilt due to surface load changes, we calculated the effects of unloading on ground deformation. The predicted tilt signals are in the range of 10-5 rad, the same order as the seismically derived tilts from this study. For 10 m of drawdown we estimate 10-5 to10-4 m of vertical ground uplift for this particular site. Our study demonstrates that dense broadband seismic arrays may offer a relatively less expensive and non-invasive tool to study temporal and spatial evolution of ground motions, and thus the ground water flow patterns, during pumping tests.

Fault diagnosis for multiple current sensors in grid‐connected inverter based on average modulation voltage

AbstractGrid‐connected inverters are the core equipment in the renewable power system. There are multiple current sensors which may affect the driving module of the switch in inverter. In multiple current sensor fault diagnosis, the coupling between fault components makes fault diagnosis difficult. This paper presents an offset fault diagnosis method of multiple current sensors based on the average modulation voltage model. Based on the influence of current sensor offset fault in the system, the average modulation voltage model is established in three‐phase stationary coordinates. The difference between the measured value of the current and the actual value is estimated through the model when the offset fault occurs. And then the fault is located. Experimental results show that the fault can be located accurately and fault tolerant control can be performed by this method when there are offset faults in multiple sensors simultaneously.

Evidence for Cenozoic topographic rejuvenation associated with the Laurel Creek Lineament in the Spruce Pine 7.5-minute quadrangle, western North Carolina, USA

ABSTRACT Linear fracture systems that strike obliquely to rock units occur across western North Carolina and are associated with Cenozoic topographic rejuvenation of the region. A 2005 earthquake on fractures within one of these systems, the Laurel Creek Lineament, suggests it's an active system. Geologic mapping of the Spruce Pine 7.5-minute quadrangle was conducted to document fault kinematics on this lineament. Near-vertical joints and faults are exposed throughout the quadrangle and dominantly strike toward 080°/260°. Fault offsets were centimeter to meter in magnitude. Fault slip indicators suggest the north block moved up relative to the south block, supported by stream knickpoints that indicate disequilibrium on the northern block. Previous studies on the Boone Lineament, north of Spruce Pine, showed that fracture system moved with the south side up. These data indicate the block of crust between the Boone Lineament and the Laurel Creek Lineament has been uplifted by these fracture systems.

Structural inheritance in Carolinia, external to the excised Rheic Suture: Savannah River Site, South Carolina, USA.

Abstract Phanerozoic reactivations of basement fault zones are documented in 5000 m of basement core recovered from beneath the updip Atlantic Coastal Plain underlying the US Department of Energy Savannah River Site (SRS) in South Carolina. These basement fault zones are adjacent to the excised Rheic Ocean suture. Meta-intrusive rocks from c. 620 and 625 Ma contain a mylonitic fabric and intrude foliated mafic metavolcanic rocks. At c. 305 Ma, granulite facies orthogneisses were thrust over amphibolite facies meta-igneous rocks in the transpressive Tinker Creek Nappe. The overturned limb of the nappe localizes the Triassic Dunbarton Basin Border Fault. The border fault acted as a conduit for fluids in the Mesozoic and Cenozoic. At c. 220 ± 5 Ma, a potassium and silica metasomatic event affected the SRS basement. A propylitic event flushed reducing fluids through rocks as young as the Santonian. The remains of a Triassic sub-basin were identified in the northwesten part of the site. A Cretaceous and younger vein paragenesis overprints the previous events. More than 30 pseudotachylytes are found in the SRS basement and are preferentially localized on metasomatized Alleghanian chloritic fractures. Pseudotachylyte post-dates mineralized fractures. The Pen Branch Fault offsets the basement–Cretaceous unconformity and is present in c. 242 m of core between PBF-7-419 m and PBF-7-660.8 m. The Pen Branch Fault cross-cuts mineralized fractures and must post-date strike-normal zeolites.

Effects of Paleoregolith and Fault Offset on the Formation of Unconformity-Type Uranium Deposits

Regional paleoregolith is found to exist immediately below unconformities separating basin fills from basement rocks in sedimentary basins. However, the controlling role of paleoregolith on unconformity-type uranium mineralization has not been quantitatively addressed before. Coupled hydrothermal fluid flow and reactive mass transport modeling are therefore performed in this study by using the software TOUGHREACT. The modeling results reveal that preferential flow occurs in the regolith due to its relatively high permeability in comparison with that of the host rocks. The thicker the regolith is, the more concentrated the fluids in the footwall of a fault zone are, leading to more compact and higher-grade deposits therein, and vice versa. Also, displacement of the regolith caused by fault offset plays an important role, as it appears to control the shape of uranium deposits. When the displacement is less than 30 m, the deposits are characterized by a more compact shape. When the displacement is over 60 m, the deposits extend more laterally and even exhibit a ‘discrete’ shape due to the expelling effect of downslope flow that occurs at the fault offset site.

3D Dynamic Rupture Modeling of the 6 February 2023, Kahramanmaraş, Turkey Mw 7.8 and 7.7 Earthquake Doublet Using Early Observations

Abstract The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The Mw 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral Mw 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data-driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.

Crustal Block and Muted Ring Development During the Formation of Mercury's Caloris Megabasin

AbstractIn contrast to multiring basins, megabasins lack clear ring structures and display crustal profiles characterized by a thin layer of crust at their center that gradually increases in thickness to the surrounding region. The 1,550‐km‐diameter Caloris basin on Mercury is often speculated to be a multiring basin; however, its crustal structure is indicative of a megabasin, perhaps explaining the lack of discernible rings. Here, we model the formation of Caloris basin using iSALE‐2D under a range of preimpact thermal conditions to determine whether the processes responsible for its megabasin crustal structure could lead to the formation of basin rings. We find that thermal gradients between 22 and 30 K/km—hotter than previously inferred—facilitate the reproduction of Caloris basin's crustal structure through crustal flowback, though such preimpact thermal structures preclude its formation as a multiring basin. Instead, repeated thrusting and necking events during transient crater collapse induce instances of fault reactivation, which diminish initial fault offsets and produce a series of discrete crustal blocks. Models assuming cooler thermal gradients lead to overly thin or non‐existent crust at the basin center, in contrast to observations. Even if crustal deficits were made up via a differentiating melt pool, the crust elsewhere would be overly thick, supporting the idea that the crustal structure of megabasins is primarily associated with transient crater collapse. We conclude that the transition to a megabasin morphology on terrestrial planetary surfaces, like the multiring transition, is dependent upon the size of the basin compared to the lithospheric thickness.

Mitigation of reverse faulting in foundation soils using geosynthetic-encased granular columns

In this study, a series of reduced model tests was conducted on soil foundations reinforced by Geosynthetic-Encased Granular Columns (GECs) placed across a reverse fault. These tests aimed at evaluating the effectiveness, reinforcing mechanism and optimal GEC horizontal spacing to mitigate the ground surface deformation associated with reverse faulting. For comparison, reduced model tests were also performed on unreinforced and Geosynthetic-Reinforced Soil (GRS) foundations. The reduced model tests were conducted to simulate a prototype 3-m-thick foundation layer subjected to a reverse fault displacement of 0.9 m. Digital Image Analysis (DIA) techniques were adopted to determine the surface displacement profile, angular distortion and shear rupture propagation considering various reverse fault offsets. Test results revealed that the GEC foundation can considerably reduce the fault-induced angular distortion at the ground surface. A reduction of 23.3% on the maximum angular distortion at the ground surface was observed as the fault displacement reached 30% of the foundation height (S/H = 30%), indicating the GEC foundation can mitigate the risk of surface fault hazards associated with large reverse fault movement. Two mechanisms, shear rupture diffusion and diversion effects, were identified for the GEC foundation, depending on the magnitude of the fault displacement and the system stiffness of the GEC foundation. Because of complex mechanisms between diversion and diffusion of the shear rupture in the GEC foundations, an optimal GEC spacing of Sh/dc = 3.3 was observed, which exhibited the most significant reduction in βmax at large fault offsets.

Using shallow hydroacoustic data to image seafloor mass transport deposits on the North West Shelf of Australia: links to neotectonics

Abstract Mass transport deposits have long been known on the Exmouth Plateau, offshore NW Australia, identified in 2D and 3D industry seismic lines. The expedition SO257 in 2017 collected 30 high-resolution, shallow seismic lines along targeted transects on the northern Australian margin. Many of these imaged mass transport deposits, with the top 700–800 m of the section captured in detail not available with industry seismic data. We present nine new high-resolution seismic lines from three separate areas of the North West Shelf. Slides in the Roebuck Basin show complex anastomosing ductile extensional mechanisms, with multiple slip surfaces and no headscarps or adjacent faults. Slides on the Exmouth Plateau have fault control, with surface fault offsets of up to 300 m, indicating seismicity as a likely triggering mechanism. Slumps along the western margin of Western Australia are more limited in extent, associated with surface notches, with indications of previous activity at depth. All areas show a repeated history of mass transport deposits. The area of the active landslide province offshore of NW Australia is far larger than the individual slides recognized on the Exmouth Plateau.

Model-Based Detection and Estimation of DC Offset of Phase Current Sensors for Field Oriented PMSM Drives

Failure of a phase current sensor in Permanent Magnet Synchronous Machine (PMSM) drive can highly impact the performance of the drive. Especially, offset faults can cause torque oscillations, which can lead to an unacceptable deterioration of the drive's performance and mechanical damages. The diagnosis of this fault condition is of fundamental importance in many applications and it is often solved by means of highly computational state estimation techniques. This paper presents a new model-based approach to systematically and accurately detect and estimate phase current sensor offsets in PMSM drives. The steady-state analytical solution of the PMSM drive model, including the control loop of the Field Oriented Control (FOC) and the sensor offset disturbance, has been found, allowing for estimating the effect of the fault on the electric drive performance. Then, a simple and effective diagnostic algorithm has been developed to detect, isolate and accurately estimate the location and the magnitude of the fault eliminating any need for state estimation or observers. Numerical and experimental results have been reported to validate the proposed model and diagnostic algorithm even in case of motor parameters variation.

The Destructive Earthquake Doublet of 6 February 2023 in South-Central Türkiye and Northwestern Syria: Initial Observations and Analyses

AbstractOn 6 February 2023, two large earthquakes with magnitude 7.8 and 7.6 rocked south-central Türkiye and northwestern Syria. At the time of writing, the death toll exceeded 50,000 in Türkiye and 7200 in Syria. The epicenter of the first mainshock was located ∼15 km east of the east Anatolian fault (EAF), the second large earthquake (9 hr later) initiated ∼90 km to the north on the east–west-trending Sürgü fault. Aftershocks delineate fault lengths of ∼350 and ∼170 km, respectively. Using satellite and seismic data for first-order analyses of surface-fault offsets, space–time rupture evolution, and recorded ground motions, our study sheds light on the reasons for the extensive destruction. The first event ruptured the EAF bilaterally, lasted for ∼80 s, and created surface fault offsets of over 6 m. The second event also ruptured bilaterally with a duration of ∼35 s and more than 7 m surface offsets. Horizontal ground accelerations reached locally up to 2g in the first mainshock; severe and widespread shaking occurred in the Hatay-Antakia area with values near 0.5g. Both earthquakes are characterized by directivity effects and abrupt rupture cessation generating stopping phases that contributed to strong seismic radiation. Shaking was further aggravated locally by site-amplification effects.

The development of multiple phases of superposed rifting in the Turkana Depression, East Africa: Evidence from receiver functions

The Turkana Depression in Eastern Africa separates the elevated plateaus of East Africa to the south and Ethiopia-Yemen to the north. It remains unclear whether the Depression lacks dynamic mantle support, or if the entire East Africa region is dynamically supported and the Depression compensated isostatically by thinned crust. Also poorly understood is how Miocene-Recent extension has developed across the Depression, connecting spatially separated magmatic rift zones in Ethiopia and Kenya. Receiver function analysis is used to constrain Moho depth and bulk-crustal VP/VS ratio below new seismograph networks in the Depression, and on the northern Tanzania craton. Crustal thickness is ∼40km below northern Uganda and 30–35km below southern Ethiopia, but 20–30km below most of the Depression, where mass-balance calculations reveal low elevations can be explained adequately by crustal thinning alone. Despite the fact that magmatism has occurred for 45Ma across the Depression, more than 15Ma before East African Rift (EAR) extension initiated, bulk crustal VP/VS across southern Ethiopia and the Turkana Depression (∼1.74) is similar to that observed in areas unaffected by Cenozoic rifting and magmatism. Evidence for voluminous lower crustal intrusions and/or melt, widespread below the Ethiopian rift and Ethiopian plateau to the north, is therefore lacking. These observations, when reviewed in light of high stretching factors (β≤2.11), suggest Cenozoic extension has been dominated until recently by faulting and plate stretching, rather than magma intrusion, which is likely an incipient process, operating directly below seismically-active Lake Turkana. Early-stage EAR basins to the west of Lake Turkana, with associated stretching factors of β≈2, formed in crust only moderately thinned during earlier rifting episodes. Conversely, ∼23km-thick crust beneath the Kino Sogo Fault Belt (KSFB) has small offset faults and thin sedimentary strata, suggesting almost all of the observed stretching occurred in Mesozoic times. Despite the KSFB marking the shortest path between focused extensional zones to the north and south, seismicity and GPS data show that modern extension is localized below Lake Turkana to the west. Failed Mesozoic rift zones, now characterized by thinned crust and relatively refractory mantle lithosphere, are being circumnavigated, not exploited by EAR rifting.

3-D Simulations of earthquakes rupture jumps: 1. Homogeneous pre-stress conditions

SUMMARYObservational and modelling studies indicate that earthquake ruptures can jump between fault sections as large as ∼3 and ∼5 km for compressional and extensional offsets, respectively. Here, we compare characteristics of the rupture jump process on parallel but offset fault sections from traditional 3-D dynamic rupture simulations governed by slip weakening friction using the finite element code, FaultMod, to those from quasi-dynamic simulations governed by rate- and state-dependent friction (rate-state friction) using the code RSQSim. These simulations use spatially uniform initial stresses. For a variety of measures the rupture renucleation position on the offset fault, the rate-state friction and slip weakening friction models produce very similar results. The principal difference is the additional occurrence of delayed rupture jumps that arise from the time- and stress-dependent nucleation that is characteristic of rate-state friction. For immediate rupture jumps, models with slip weakening friction span greater offsets than those with rate-state friction. However the jump distances are nearly identical when delayed rupture jumps are included in the comparisons. We propose that delayed rupture jumps are the likely mechanism for adjacent large-earthquake pairs and clusters. Based on the similarity of renucleation positions with both dynamic and quasi-dynamic models, we conclude that the renucleation positions for rupture initiation on the receiver fault (separated by less than ∼3 km from the source fault) are primarily controlled by static stress changes induced by slip on the initiating fault. However, in light of the slightly greater maximum jump distances (>3 km) seen with the dynamic slip weakening friction model, dynamic stress changes from seismic waves play an increasingly important role as offset distances increase.

From pilot knowledge via integrated reservoir characterization to utilization perspectives of deep geothermal reservoirs: the 3D model of Groß Schönebeck (North German Basin)

The Groß Schönebeck site in the North German Basin serves as research platform to study the geothermal potential of deeply buried Permian reservoir rocks and the technical feasibility of heat extraction. The structural setting of the site was investigated in more detail by a newly acquired 3D-seismic survey to improve the former conceptual model that was based on several old 2D seismic lines. The new data allow a revision of the geological interpretation, enabling the setup of a new reservoir model and providing base information for a possible further site development of Permo-Carboniferous targets. The 3D seismic allows for the first time a consistent geological interpretation and model parameterization of the well-studied geothermal site. Main reflector horizons and the corresponding stratigraphic units were mapped and the structural pattern of the subsurface presented in the 8 km × 8 km × 4 km large seismic volume. Attribute analysis revealed some fracture and fault patterns in the upper Zechstein and post-Permian units, while formerly hypothesized large offset faults are not present in the Rotliegend reservoir. However, a well-established graben-like structure at the top of the Zechstein succession is most likely related to broken anhydritic brittle intra-salt layers of some meter of thickness. Most reflectors above the salt show a rather undisturbed pattern. The main reservoir sandstone of the Dethlingen Formation (Rotliegend) was mapped and characterized. The base of the underlying Permo-Carboniferous volcanic rock sequence and hence its thickness could not be depicted reliably from the geophysical data. Based on the seismic data and the available reconnaissance drilling, logging, and laboratory data of the Groß Schönebeck research site, the thickness and distribution of the sedimentary Rotliegend (with emphasis of the sandy reservoir section) and of the volcanic rock sequence was modelled and stochastically parameterized with petrophysical properties guided by seismic facies pattern correlation, providing a more realistic reservoir description. Properties include total and effective porosity, permeability, bulk density, thermal conductivity, thermal diffusivity, and specific heat capacity. The data and interpretation constitute the basis for a better understanding of the thermo and hydromechanical processes at the site and for future measures. Further site development could include a deepening of one well to provide evidence on the volcanic rock sequence and consider deviated wells into favourable zones and the design of a fracture-dominated utilization approach.