Search

The seismogenesis of rocks is mainly affected by their mineral composition and in situ conditions (temperature and state of stress). Diverse laboratory experiments have explored the frictional behaviour of the rocks and rock-forming minerals most common in the crust and uppermost mantle. However, it is debated how to "upscale" these results to the lithosphere. In particular, most earthquakes in the crust nucleate down to the crustal seismogenic depth (CSD), which is a proxy for the maximum depth of crustal earthquake ruptures in seismic hazard assessments. In this study we propose a workflow to upscale and validate those laboratory experiments to natural geological conditions relevant for crustal and upper-mantle rocks. We used the southern Caribbean and northwestern South America as a case study to explore the three-dimensional spatial variation of the CSD (mapped as D90, the 90 % percentile of hypocentral depths) and the temperatures at which crustal earthquakes likely occur. A 3D steady-state thermal field was computed for the region with a finite-element scheme using the software GOLEM, considering the uppermost 75 km of a previously published 3D data-integrative lithospheric configuration, lithology-constrained thermal parameters, and appropriate upper and lower boundary conditions. The model was validated using additional, independent measurements of downhole temperatures and heat flow. We found that the majority of crustal earthquakes nucleate at temperatures less than 350 ∘ C, in agreement with frictional experiments of typical crustal rocks. A few outliers with larger hypocentral temperatures evidence nucleation conditions consistent with the seismogenic window of olivine-rich rocks, and can be due to either uncertainties in the Moho depths and/or in the earthquake hypocentres or the presence of ultramafic rocks within different crustal blocks and allochthonous terranes accreted to this complex margin. Moreover, the spatial distribution of crustal seismicity in the region correlates with the geothermal gradient, with no crustal earthquakes occurring in domains with low thermal gradient. Finally, we find that the largest earthquake recorded in the region (Mw=7.1 , Murindó sequence, in 1992) nucleated close to the CSD, highlighting the importance of considering this lower-stability transition for seismogenesis when characterizing the depth of seismogenic sources in hazard assessments. The approach presented in this study goes beyond a statistical approach in that the local heterogeneity of physical properties is considered in our simulations and additionally validated by the observed depth distribution of earthquakes. The coherence of the calculated hypocentral temperatures with those expected from laboratory measurements provides additional support to our modelling workflow. This approach can be applied to other tectonic settings worldwide, and it could be further refined as new, high-quality hypocentral locations and heat flow and temperature observations become available. [ABSTRACT FROM AUTHOR]

Seismic hazard assessment in slow straining regions is challenging because earthquake catalogues only record events from approximately the last 100 years, whereas earthquake recurrence times on individual faults can exceed 1000 years. Systematic mapping of active faults allows fault sources to be used within probabilistic seismic hazard assessment, which overcomes the problems of short-term earthquake records. We use Shuttle Radar Topography Mission (SRTM) data to analyse surface deformation in the Luangwa Rift in Zambia and develop the Luangwa Rift Active Fault Database (LRAFD). The LRAFD is an open-source geospatial database containing active fault traces and their attributes and is freely available at 10.5281/zenodo.6513691. We identified 18 faults that display evidence for Quaternary activity, and empirical relationships suggest that these faults could cause earthquakes up to Mw 8.1, which would exceed the magnitude of historically recorded events in southern Africa. On the four most prominent faults, the median height of Quaternary fault scarps varies between 12.9 ± 0.4 and 19.2 ± 0.9 m, which suggests they were formed by multiple earthquakes. Deformation is focused on the edges of the Luangwa Rift: the most prominent Quaternary fault scarps occur along the 207 km long Chipola and 142 km long Molaza faults, which are the rift border faults and the longest faults in the region. We associate the scarp on the Molaza Fault with possible surface ruptures from two 20th century earthquakes. Thus, the LRAFD reveals new insights into active faulting in southern Africa and presents a framework for evaluating future seismic hazard. [ABSTRACT FROM AUTHOR]

Crustal shear zones generate earthquakes, which are at present unpredictable, but advances in mechanistic understanding of the earthquake cycle offer hope for future advances in earthquake forecasting. Studies of fault zone architecture have the potential to reveal the controls on fault rupture, locking, and reloading that control the temporal and spatial patterns of earthquakes. The Pofadder Shear Zone exposed in the Orange River in South Africa is an ancient, exhumed, paleoseismogenic continental transform which preserves the architecture of the earthquake source near the base of the seismogenic zone. To investigate the controls on earthquake rupture geometries in the seismogenic crust, we produced a high-resolution geologic map of the shear zone core mylonite zone. The core consists of ∼ 1–200 cm , pinch-and-swell layers of mylonites of variable mineralogic composition, reflecting the diversity of regional rock types which were dragged into the shear zone. Our map displays centimetric layers of a unique black ultramylonite along some mylonite interfaces, locally adding to thick composite layers suggesting reactivation or bifurcation. We present a set of criteria for identifying recrystallised pseudotachylytes (preserved earthquake frictional melts) and show that the black ultramylonite is a recrystallised pseudotachylyte, with its distribution representing a map of ancient earthquake rupture surfaces. Pseudotachylytes are most abundant on interfaces between the strongest wall rocks. We find that the geometry of lithologic interfaces which hosted earthquakes differs from interfaces lacking pseudotachylyte at wavelengths of ≳ 10 m. We argue that the pinch-and-swell structure of the mylonitic layering, enhanced by viscosity contrasts between layers of different mineralogy, is expected to generate spatially heterogeneous stress during viscous creep in the shear zone, which dictated the path of earthquake ruptures. The condition of rheologically layered materials causing heterogeneous stresses should be reasonably expected in any major shear zone, is enhanced by creep, and represents the pre-seismic background conditions through which earthquakes nucleate and propagate. This has implications for patterns of earthquake recurrence and explains why some potential geologic surfaces are favored for earthquake rupture over others. [ABSTRACT FROM AUTHOR]

The Dinaric Fault System in western Slovenia, consisting of NW–SE-trending, right-lateral strike-slip faults, accommodates the northward motion of Adria with respect to Eurasia. These active faults show a clear imprint in the morphology, and some of them hosted moderate instrumental earthquakes. However, it is largely unknown if the faults also had strong earthquakes in the late Quaternary. This hampers our understanding of the regional tectonics and the seismic hazard. Geological evidence of co-seismic surface ruptures only exists for one historical event, the 1511 Idrija earthquake with a magnitude of ∼ M 6.8, but the causative fault is still disputed. Here we use geomorphological data, near-surface geophysical surveys, and paleoseismological trenching to study two of these faults: the Predjama Fault and the Idrija Fault. In a paleoseismological trench across the Predjama Fault we found deformation features that may have been caused by an earthquake between 13–0.7 ka, very likely not earlier than 8.4 ka. At the Idrija Fault, a surface-rupturing earthquake happened around 2.5 ka. We show that instrumental and historical seismicity data do not capture the strongest events in this area. [ABSTRACT FROM AUTHOR]

Understanding diagenetic reactions in accreted sediments is critical for establishing the balance of fluid sources and sinks in accretionary prisms, which is in turn important for assessing the fluid pressure field and the ability for faults to host seismic slip. For this reason, we studied diagenetic reactions in deformation bands (shear zones and veins) within deep mud sediments from the Nankai accretionary prism (SW Japan) drilled at site C0001 during IODP Expedition 315, by means of microscopic observation, X-ray diffraction, and major- and trace-element analyses. Deformation bands are not only more compacted than the host sediment but are also enriched in framboidal pyrite, as observed under microscopy and confirmed by chalcophile-element enrichments (Fe, S, Cu, As, Sb, Pb). In tandem, one shear zone sample displays a destabilization of smectite or illite–smectite mixed layers and a slight crystallization of illite relative to its sediment matrix, and another sample shows correlated increases in B and Li in shear zones and veins compared to the host sediment, both effects suggesting a transformation of smectite into illite in deformation bands. The two diagenetic reactions of sulfide precipitation and smectite-to-illite transformation are explained by a combined action of sulfate-reducing and methanogen bacteria, which strongly suggests an increased activity of anaerobic microbial communities localized in deformation bands. This local bacterial proliferation was possibly enhanced by the liberation of hydrogen from strained phyllosilicates. We suggest that the proliferation of anoxic bacteria, boosted by deformation, may contribute to the pore water freshening observed at depth in accretionary prisms. Deformation-enhanced metabolic reactions may also explain the illitization observed in major faults of accretionary prisms. Care is therefore needed before interpreting illitization, and other diagenetic reactions as well, as evidence of shear heating, as these might be biogenic instead of thermogenic. [ABSTRACT FROM AUTHOR]

With the aim of unveiling evidence of Late Quaternary faulting, a series of ground-penetrating radar (GPR) profiles were acquired across the southern portion of the Fosso della Valle–Campotenese normal fault (VCT), located at the Campotenese continental basin (Mt. Pollino region) in the southern Apennines active extensional belt (Italy). A set of 49 GPR profiles, traced nearly perpendicular to this normal fault, was acquired using 300 and 500 MHz antennas and carefully processed through a customized workflow. The data interpretation allowed us to reconstruct a pseudo-3D model depicting the boundary between the Mesozoic bedrock and the sedimentary fill of the basin, which were in close proximity to the fault. Once the GPR signature of faulting was reviewed and defined, we interpret near-surface alluvial and colluvial sediments dislocated by a set of conjugate (W- and E-dipping) discontinuities that penetrate inside the underlying Triassic dolostones. Close to the contact between the continental deposits and the bedrock, some buried scarps which offset wedge-shaped deposits are interpreted as coseismic ruptures, subsequently sealed by later deposits. Our pseudo-3D GPR dataset represented a good trade-off between a dense 3D-GPR volume and conventional 2D data, which normally requires a higher degree of subjectivity during the interpretation. We have thus reconstructed a reliable subsurface fault pattern, discriminating master faults and a series of secondary splays. This contribution better characterizes active Quaternary faults in an area which falls within the Pollino seismic gap and is considered prone to severe surface faulting. Our results encourage further research at the study site, whilst we also recommend our workflow for similar regions characterized by high seismic hazard and scarcity of near-surface geophysical data. [ABSTRACT FROM AUTHOR]

The region around the town Albstadt, SW Germany, was struck by four damaging earthquakes with magnitudes greater than 5 during the last century. These earthquakes occurred along the Albstadt Shear Zone (ASZ), which is characterized by more or less continuous microseismicity. As there are no visible surface ruptures that may be connected to the fault zone, we study its characteristics by its seismicity distribution and faulting pattern. We use the earthquake data of the state earthquake service of Baden-Württemberg from 2011 to 2018 and complement it with additional phase picks beginning in 2016 at the AlpArray and StressTransfer seismic networks in the vicinity of the ASZ. This extended data set is used to determine new minimum 1-D seismic vp and vs velocity models and corresponding station delay times for earthquake relocation. Fault plane solutions are determined for selected events, and the principal stress directions are derived. The minimum 1-D seismic velocity models have a simple and stable layering with increasing velocity with depth in the upper crust. The corresponding station delay times can be explained well by the lateral depth variation of the crystalline basement. The relocated events align about north–south with most of the seismic activity between the towns of Tübingen and Albstadt, east of the 9 ∘ E meridian. The events can be separated into several subclusters that indicate a segmentation of the ASZ. The majority of the 25 determined fault plane solutions feature an NNE–SSW strike but NNW–SSE-striking fault planes are also observed. The main fault plane associated with the ASZ dips steeply, and the rake indicates mainly sinistral strike-slip, but we also find minor components of normal and reverse faulting. The determined direction of the maximum horizontal stress of 140–149 ∘ is in good agreement with prior studies. Down to ca. 7–8 km depth SHmax is bigger than SV ; below this depth, SV is the main stress component. The direction of SHmax indicates that the stress field in the area of the ASZ is mainly generated by the regional plate driving forces and the Alpine topography. [ABSTRACT FROM AUTHOR]

Significant uncertainties occur through varying methodologies when interpreting faults using seismic data. These uncertainties are carried through to the interpretation of how faults may act as baffles or barriers, or increase fluid flow. How fault segments are picked when interpreting structures, i.e. which seismic line orientation, bin spacing and line spacing are specified, as well as what surface generation algorithm is used, will dictate how rugose the surface is and hence will impact any further interpretation such as fault seal or fault growth models. We can observe that an optimum spacing for fault interpretation for this case study is set at approximately 100 m, both for accuracy of analysis but also for considering time invested. It appears that any additional detail through interpretation with a line spacing of ≤ 50 m adds complexity associated with sensitivities by the individual interpreter. Further, the locations of all seismic-scale fault segmentation identified on throw–distance plots using the finest line spacing are also observed when 100 m line spacing is used. Hence, interpreting at a finer scale may not necessarily improve the subsurface model and any related analysis but in fact lead to the production of very rough surfaces, which impacts any further fault analysis. Interpreting on spacing greater than 100 m often leads to overly smoothed fault surfaces that miss details that could be crucial, both for fault seal as well as for fault growth models. Uncertainty in seismic interpretation methodology will follow through to fault seal analysis, specifically for analysis of whether in situ stresses combined with increased pressure through CO 2 injection will act to reactivate the faults, leading to up-fault fluid flow. We have shown that changing picking strategies alter the interpreted stability of the fault, where picking with an increased line spacing has shown to increase the overall fault stability. Picking strategy has shown to have a minor, although potentially crucial, impact on the predicted shale gouge ratio. [ABSTRACT FROM AUTHOR]

Coseismic surface faulting is a significant source of hazard for critical plants and distributive infrastructure; it may occur either on the principal fault or as distributed rupture on nearby faults. Hazard assessment for distributed faulting is based on empirical relations which, in the case of normal faults, were derived almost 15 years ago using a dataset of US earthquakes. We collected additional case histories worldwide, for a total of 21 earthquakes, and calculated the conditional probability of distributed faulting as a function of distance from the principal fault. We found no clear dependency on the magnitude nor the time of occurrence of the earthquakes, but our data consistently show a higher probability of rupture when compared with the scaling relations currently adopted in engineering practice. We derive updated empirical regressions and show that the results are strongly conditioned by the averaging of earthquakes effectively generating distributed faulting at a given distance and those which did not generate faulting; thus, we introduce a more conservative scenario that can be included in a logic tree approach to consider the full spectrum of potential ruptures. Our results can be applied in the framework of probabilistic assessment of fault displacement hazard. [ABSTRACT FROM AUTHOR]

Seismic faulting causes wall rock damage, which is driven by both mechanical and thermal stress. In the lower crust, co-seismic damage increases wall rock permeability, permits fluid infiltration and triggers metamorphic reactions that transform rock rheology. Wall rock microstructures reveal high-stress conditions near earthquake faults; however, there is limited documentation on the effects of a thermal pulse coupled with fluid infiltration. Here, we present a transmission electron microscopy study of co-seismic microfractures in plagioclase feldspar from lower crustal granulites from the Bergen Arcs, Western Norway. Focused ion beam foils are collected 1.25 mm and 1.8 cm from a 1.3 mm thick eclogite facies pseudotachylyte vein. Dislocation-free plagioclase and K-feldspar aggregates in the microfractures record a history of fluid introduction and recovery from a short-lived high-stress state caused by slip along the nearby fault. The feldspar aggregates retain the crystallographic orientation of their host and are elongated subparallel to the pseudotachylyte. We propose that plagioclase partially amorphized along the microfractures at peak stress conditions followed by repolymerization to form dislocation-free grain aggregates. Repolymerization and recrystallization were enhanced by the infiltration of fluids that transported Ca and K into the microfractures. Subsequent cooling led to exsolution of intermediate plagioclase compositions and the formation of the Bøggild–Huttenlocher intergrowth in the grains from the fracture closest to the pseudotachylyte. Our findings provide unequivocal evidence that the introduction of fluids in the microfractures occurred within the timescale of the thermal perturbation, prompting rapid annealing of damaged wall rock soon after earthquake rupture. [ABSTRACT FROM AUTHOR]

Substantial insight into earthquake source processes has resulted from considering frictional ruptures analogous to cohesive-zone shear cracks from fracture mechanics. This analogy holds for slip-weakening representations of fault friction that encapsulate the resistance to rupture propagation in the form of breakdown energy, analogous to fracture energy, prescribed in advance as if it were a material property of the fault interface. Here, we use numerical models of earthquake sequences with enhanced weakening due to thermal pressurization of pore fluids to show how accounting for thermo-hydro-mechanical processes during dynamic shear ruptures makes breakdown energy rupture-dependent. We find that local breakdown energy is neither a constant material property nor uniquely defined by the amount of slip attained during rupture, but depends on how that slip is achieved through the history of slip rate and dynamic stress changes during the rupture process. As a consequence, the frictional breakdown energy of the same location along the fault can vary significantly in different earthquake ruptures that pass through. These results suggest the need to reexamine the assumption of predetermined frictional breakdown energy common in dynamic rupture modeling and to better understand the factors that control rupture dynamics in the presence of thermo-hydro-mechanical processes. [ABSTRACT FROM AUTHOR]

The strength properties of fault rocks at shearing rates spanning the transition from crystal–plastic flow to frictional slip play a central role in determining the distribution of crustal stress, strain, and seismicity in tectonically active regions. We review experimental and microphysical modelling work, which is aimed at elucidating the processes that control the transition from pervasive ductile flow of fault rock to rate-and-state-dependent frictional (RSF) slip and to runaway rupture, carried out at Utrecht University in the past 2 decades or so. We address shear experiments on simulated gouges composed of calcite, halite–phyllosilicate mixtures, and phyllosilicate–quartz mixtures performed under laboratory conditions spanning the brittle–ductile transition. With increasing shear rate (or decreasing temperature), the results consistently show transitions from (1) stable velocity-strengthening (v -strengthening) behaviour, to potentially unstable v -weakening behaviour, and (2) back to v strengthening. Sample microstructures show that the first transition seen at low shear rates and/or high temperatures represents a switch from pervasive, fully ductile deformation to frictional sliding involving dilatant granular flow in localized shear bands where intergranular slip is incompletely accommodated by creep of individual mineral grains. A recent microphysical model, which treats fault rock deformation as controlled by competition between rate-sensitive (diffusional or crystal–plastic) deformation of individual grains and rate-insensitive sliding interactions between grains (granular flow), predicts both transitions well. Unlike classical RSF approaches, this model quantitatively reproduces a wide range of (transient) frictional behaviours using input parameters with direct physical meaning, with the latest progress focusing on incorporation of dynamic weakening processes characterizing co-seismic fault rupture. When implemented in numerical codes for crustal fault slip, the model offers a single unified framework for understanding slip patch nucleation and growth to critical (seismogenic) dimensions, as well as for simulating the entire seismic cycle. [ABSTRACT FROM AUTHOR]

Faults and fractures represent unique features of the solid Earth and are especially pervasive in the shallow crust. Aside from directly relating to crustal dynamics and the systematic assessment of associated risk, fault and fracture networks enable the efficient migration of fluids and therefore have a direct impact on concrete topics relevant to society, including climate-change-mitigating measures like CO2 sequestration or geothermal exploration and production. Due to their small-scale complexity, fault zones and fracture networks are typically poorly resolved, and their presence can often only be inferred indirectly in seismic and ground-penetrating radar (GPR) subsurface reconstructions. We suggest a largely data-driven framework for the direct imaging of these features by making use of the faint and still often underexplored diffracted portion of the wave field. Finding inspiration in the fields of optics and visual perception, we introduce two different conceptual pathways for coherent diffraction imaging and discuss respective advantages and disadvantages in different contexts of application. At the heart of both of these strategies lies the assessment of data coherence, for which a range of quantitative measures is introduced. To illustrate the versatility and effectiveness of the approach for high-resolution geophysical imaging, several seismic and GPR field data examples are presented, in which the diffracted wave field sheds new light on crustal features like fluvial channels, erosional surfaces, and intricate fault and fracture networks on land and in the marine environment. [ABSTRACT FROM AUTHOR]

Travertine deposits present above the St. Johns Dome natural CO2 reservoir in Arizona, USA, document a long (>400 kyr) history of surface leakage of CO2 from a subsurface reservoir. These deposits are concentrated along surface traces of faults, implying that there has been a structural control on the migration pathway of CO2 -rich fluids. Here, we combine slip tendency and fracture stability to analyse the geomechanical stability of the reservoir-bounding Coyote Wash Fault for three different stress fields and two interpreted fault rock types to predict areas with high leakage risks. We find that these areas coincide with the travertine deposits on the surface, indicating that high-permeability pathways as a result of critically stressed fracture networks exist in both a fault damage zone and around a fault tip. We conclude that these structural features control leakage. Importantly, we find that even without in situ stress field data, the known leakage points can be predicted using geomechanical analyses, despite the unconstrained tectonic setting. Whilst acquiring high-quality stress field data for secure subsurface CO2 or energy storage remains critical, we shown that a first-order assessment of leakage risks during site selection can be made with limited stress field knowledge. [ABSTRACT FROM AUTHOR]

Natural fault networks are geometrically complex systems that evolve through time. The evolution of faults and their off-fault damage patterns are influenced by both dynamic earthquake ruptures and aseismic deformation in the interseismic period. To better understand each of their contributions to faulting we simulate both earthquake rupture dynamics and long-term deformation in a visco-elasto-plastic crust subjected to rate- and state-dependent friction. The continuum mechanics-based numerical model presented here includes three new features. First, a 2.5-D approximation is created to incorporate the effects of a viscoelastic lower crustal substrate below a finite depth. Second, we introduce a dynamically adaptive (slip-velocity-dependent) measure of fault width to ensure grid size convergence of fault angles for evolving faults. Third, fault localization is facilitated by plastic strain weakening of bulk rate and state friction parameters as inspired by laboratory experiments. This allows us to simulate sequences of episodic fault growth due to earthquakes and aseismic creep for the first time. Localized fault growth is simulated for four bulk rheologies ranging from persistent velocity weakening to velocity strengthening. Interestingly, in each of these bulk rheologies, faults predominantly localize and grow due to aseismic deformation. Yet, cyclic fault growth at more realistic growth rates is obtained for a bulk rheology that transitions from velocity-strengthening friction to velocity-weakening friction. Fault growth occurs under Riedel and conjugate angles and transitions towards wing cracks. Off-fault deformation, both distributed and localized, is typically formed during dynamic earthquake ruptures. Simulated off-fault deformation structures range from fan-shaped distributed deformation to localized splay faults. We observe that the fault-normal width of the outer damage zone saturates with increasing fault length due to the finite depth of the seismogenic zone. We also observe that dynamically and statically evolving stress fields from neighboring fault strands affect primary and secondary fault growth and thus that normal stress variations affect earthquake sequences. Finally, we find that the amount of off-fault deformation distinctly depends on the degree of optimality of a fault with respect to the prevailing but dynamically changing stress field. Typically, we simulate off-fault deformation on faults parallel to the loading direction. This produces a 6.5-fold higher off-fault energy dissipation than on an optimally oriented fault, which in turn has a 1.5-fold larger stress drop. The misalignment of the fault with respect to the static stress field thus facilitates off-fault deformation. These results imply that fault geometries bend, individual fault strands interact, and optimal orientations and off-fault deformation vary through space and time. With our work we establish the basis for simulations and analyses of complex evolving fault networks subject to both long-term and short-term dynamics. [ABSTRACT FROM AUTHOR]

Monitoring of CO2 degassing in seismoactive areas allows the study of correlations of gas release and seismic activity. Reliable continuous monitoring of the gas flow rate in rough field conditions requires robust methods capable of measuring gas flow at different types of gas outlets such as wet mofettes, mineral springs, and boreholes. In this paper we focus on the methods and results of the long-term monitoring of CO2 degassing in the West Bohemia/Vogtland region in central Europe, which is typified by the occurrence of earthquake swarms and discharge of carbon dioxide of magmatic origin. Besides direct flow measurement using flowmeters, we introduce a novel indirect technique based on quantifying the gas bubble contents in a water column, which is capable of functioning in severe environmental conditions. The method calculates the mean bubble fraction in a water–gas mixture from the pressure difference along a fixed depth interval in a water column. Laboratory tests indicate the nonlinear dependence of the bubble fraction on the flow rate, which is confirmed by empirical models found in the chemical and nuclear engineering literature. Application of the method in a pilot borehole shows a high correlation between the bubble fraction and measured gas flow rate. This was specifically the case for two coseismic anomalies in 2008 and 2014, when the flow rate rose during a seismic swarm to a multitude of the preseismic level for several months and was followed by a long-term flow rate decline. However, three more seismic swarms occurring in the same fault zone were not associated with any significant CO2 flow anomaly. We surmise that this could be related to the slightly farther distance of the hypocenters of these swarms compared to the two ones which caused the coseismic CO2 flow rise. Further long-term CO2 -flow monitoring is required to verify the mutual influence of CO2 degassing and seismic activity in the area. [ABSTRACT FROM AUTHOR]

Normal faulting drives extensional growth folding of the Earth's upper crust during continental extension, yet we know little of how fold geometry relates to the structural segmentation of the underlying fault. We use field data from the Hadahid Fault System, Suez Rift, Egypt, to investigate the geometry and kinematics of a large (30 km long, up to 2.5 km displacement), exceptionally well-exposed normal fault system and to test and develop models for extensional growth folding. The Hadahid Fault System comprises eight up to 5 km long segments that are defined by unbreached or breached monoclines. These segments are soft-linked, hard-linked, or defined by a more subtle along-strike transition in overall structural style. High overlap : separation (O:S) ratios between its segments suggest the Hadahid Fault System comprises a single, now hard-linked structure at depth. We demonstrate that a progressive loss of at-surface displacement along-strike of the Hadahid Fault System results in surface-breaking faults and breached monoclines being replaced by unbreached monoclines developed above blind faults. However, shorter along-strike length-scale variations in structural style also occur, with unbreached monoclines developed between breached monoclines. The origin of this variability is unclear, but it might reflect local variations in host rock material properties that drive short length-scale variations in fault propagation-to-slip ratio, and thus the timing and location of fold breaching. We show that folding is a key expression of the strain that accumulates in areas of continental extension, arguing that tectono-sedimentary models for rift development should capture the related structural complexity. [ABSTRACT FROM AUTHOR]

A shallow Mw 5.3 earthquake near Lake Muir in the stable continental region (SCR) crust of southwest Western Australia on the 16 September 2018 was followed on the 8 November by a proximal Mw 5.2 event. Focal mechanisms produced for the events suggest reverse and strike-slip rupture, respectively. Field mapping, guided by Sentinel-1 InSAR data, reveals that the first event produced an approximately 3 km long and up to 0.4–0.6 m high west-facing surface rupture, consistent with reverse slip on a moderately east-dipping fault. The InSAR data also show that the surface scarp relates to a subsurface rupture ∼ 5 km long, bound at its northern and southern extremities by bedrock structures. The November event produced a surface deformation envelope that is spatially coincident with that of the September event but did not result in discrete surface rupture. Almost 900 aftershocks were recorded by a temporary seismometer deployment. Hypocentre locations correlate poorly with the rupture plane of their respective mainshocks but correlate well with regions of increased Coulomb stress. The spatial and temporal relationships between the Mw>5.0 events and their aftershocks reveals dependencies with implications for how other less well-documented SCR earthquake sequences could be interpreted. Furthermore, the September Mw 5.3 Lake Muir earthquake was the ninth event documented to have produced surface rupture in Australia in historical times. These nine ruptures are located exclusively in the Precambrian non-extended SCR rocks of central and western Australia, and none could have been identified and mapped using topographic signature prior to the historical event. Consistent, though fragmentary, evidence exists from analogous regions worldwide. Our analysis of the Lake Muir earthquake sequence therefore provides constraint on models describing mechanisms for strain accumulation and localized release as earthquakes in non-extended SRC crust. [ABSTRACT FROM AUTHOR]

During September–October of 2013 an intense swarm of earthquakes occurred off the east coast of Spain associated with the injection of the base gas in an offshore underground gas storage. Two weeks after the end of the injection operations, three moderate-sized earthquakes (Mw 4.0–4.1) occurred near the storage. These events were widely felt by the nearby population, leading to the indefinite shut-down of the facility. Here we investigate the source parameters (focal depth and mechanism) of the largest earthquakes in the sequence in order to identify the faults reactivated by the gas injection and to help understand the processes that caused the earthquakes. Our waveform modeling results indicate that the largest earthquakes occurred at depths of 6–8 km beneath the sea floor, significantly deeper than the injection depth (∼1800 m). Although we cannot undoubtedly discriminate the fault plane from the two nodal planes of the mechanisms, most evidence seems to favor a NW–SE-striking fault plane. We propose that the gas injection reactivated faults in the Paleozoic basement, with regional orientation possibly inherited from the opening of the Valencia Trough. [ABSTRACT FROM AUTHOR]

Empirical scaling relationships between fault or slip dimensions and earthquake magnitudes are often used to assess the maximum possible earthquake magnitude of a territory. Upon the assumption of the reactivability of any fault, the earthquake magnitudes derived from the surface fault length (FLEM) are compared at the national scale in Italy against catalogued magnitudes. FLEMs are obtained by considering a comprehensive fault dataset regardless of fault age, stress field orientation, strain rate, etc. In particular, (1) a comprehensive catalogue of all known faults is compiled by merging the most complete databases available; (2) FLEM is then derived from fault length; and (3) the resulting FLEMs are compared (i.e. the mathematical difference) with catalogued earthquake magnitudes. Results show that the largest FLEMs as well as the largest differences between FLEMs and catalogued magnitudes are observed for poorly constrained faults, mainly inferred from subsurface data. It is suggested that these areas have to be further characterized to better estimate fault dimension and segmentation and hence properly assess the FLEM. Where, in contrast, the knowledge of faults is geologically well constrained, the calculated FLEM is often consistent with the catalogued seismicity, with the 2 σ value of the distribution of differences being 1.47 and reducing to 0.53 when considering only the Mw≥6.5 earthquakes. Our work highlights areas, in Italy, where further detailed studies on faults are required. [ABSTRACT FROM AUTHOR]

Laboratory experiments involving unconfined compressive failure of borosilicate glass cylinders quantified the elastic strain energy released at failure and the size distribution of the resulting fragments. The data were carefully assessed for potential inaccuracies in surface-area calculation, the contribution of energy from the compression machine relaxation during specimen failure, and possible variations in the specific fracture energy of the specimens. The data showed that more new surface area was created during the failures than would be possible if the long-standing assumption, which is that all the energy involved in creating new rock surface area in brittle material is taken up by the newly created surfaces as surface potential energy and is not available to do further work, were valid. We therefore conclude that the assumption is false. This conclusion is supported by independent data from a previous investigation whose authors did not pursue this particular application. Our result does not affect the validity of Griffith fracture mechanics and is significant only when large numbers of very fine fragments are created by brittle fracture, as in rock-avalanche motion and earthquake rupture, and are identified in particle-size distributions. In such situations our result is very significant to understanding fracture energetics. [ABSTRACT FROM AUTHOR]

Geologic carbon storage, as well as other geo-energy applications, such as geothermal energy, seasonal natural gas storage and subsurface energy storage imply fluid injection and/or extraction that causes changes in rock stress field and may induce (micro)seismicity. If felt, seismicity has a negative effect on public perception and may jeopardize wellbore stability and damage infrastructure. Thus, induced earthquakes should be minimized to successfully deploy geo-energies. However, numerous processes may trigger induced seismicity, which contribute to making it complex and translates into a limited forecast ability of current predictive models. We review the triggering mechanisms of induced seismicity. Specifically, we analyze (1) the impact of pore pressure evolution and the effect that properties of the injected fluid have on fracture and/or fault stability; (2) non-isothermal effects caused by the fact that the injected fluid usually reaches the injection formation at a lower temperature than that of the rock, inducing rock contraction, thermal stress reduction and stress redistribution around the cooled region; (3) local stress changes induced when low-permeability faults cross the injection formation, which may reduce their stability and eventually cause fault reactivation; (4) stress transfer caused by seismic or aseismic slip; and (5) geochemical effects, which may be especially relevant in carbonate-containing formations. We also review characterization techniques developed by the authors to reduce the uncertainty in rock properties and subsurface heterogeneity both for the screening of injection sites and for the operation of projects. Based on the review, we propose a methodology based on proper site characterization, monitoring and pressure management to minimize induced seismicity. [ABSTRACT FROM AUTHOR]

New field work, combined with analysis of high-resolution aerial photographs, digital elevation models, and satellite imagery, has identified an active fault that is traceable for ∼90 km across the Seymchan Basin and is part of the Ulakhan fault system, which is believed to form the Okhotsk–North America plate boundary. Age dating of alluvial fan sediments in a channel system that is disturbed by fault activity suggests the current scarp is a result of a series of large earthquakes (≥Mw7.5) that have occurred since 11.6±2.7 ka. A possible channel feature offset by 62±4 m associated with these sediments yields a slip rate of 5.3±1.3 mm yr -1 , in broad agreement with rates suggested from global plate tectonics. Our results clearly identify the Ulakhan fault as the Okhotsk–North America plate boundary and show that tectonic strain release is strongly concentrated on the boundaries of Okhotsk. In light of our results, the likelihood of recurrence of Mw7.5 earthquakes is high, suggesting a previously underestimated seismic hazard across the region. [ABSTRACT FROM AUTHOR]

Along-strike variation in scarp morphology reflects differences in a fault's geomorphic and structural development and can thus indicate fault rupture history and mechanical segmentation. Parameters that define scarp morphology (height, width, slope) are typically measured or calculated manually. The time-consuming manual approach reduces the density and objectivity of measurements and can lead to oversight of small-scale morphological variations that occur at a resolution impractical to capture. Furthermore, inconsistencies in the manual approach may also lead to unknown discrepancies and uncertainties between, and also within, individual fault scarp studies. Here, we aim to improve the efficiency, transparency and uniformity of calculating scarp morphological parameters by developing a semi-automated Scarp PARameTer Algorithm (SPARTA). We compare our findings against a traditional, manual analysis and assess the performance of the algorithm using a range of digital elevation model (DEM) resolutions. We then apply our new algorithm to a 12 m resolution TanDEM-X DEM for four southern Malawi fault scarps, located at the southern end of the East African Rift system: the Bilila–Mtakataka fault (BMF) and three previously unreported scarps – Thyolo, Muona and Malombe. All but Muona exhibit first-order structural segmentation at their surface. By using a 5 m resolution DEM derived from high-resolution (50 cm pixel -1) Pleiades stereo-satellite imagery for the Bilila–Mtakataka fault scarp, we quantify secondary structural segmentation. Our scarp height calculations from all four fault scarps suggest that if each scarp was formed by a single, complete rupture, the slip–length ratio for each earthquake exceeds the maximum typical value observed in historical normal faulting earthquakes around the world. The high slip–length ratios therefore imply that the Malawi fault scarps likely formed in multiple earthquakes. The scarp height distribution implies the structural segments of both the BMF and Thyolo fault have merged via rupture of discrete faults (hard links) through several earthquake cycles, and the segments of the Malombe fault have connected via distributed deformation zones (soft links). For all faults studied here, the length of earthquake ruptures may therefore exceed the length of each segment. Thus, our findings shed new light on the seismic hazard in southern Malawi, indicating evidence for a number of large (Mw 7–8) prehistoric earthquakes, as well as providing a new semi-automated methodology (SPARTA) for calculating scarp morphological parameters, which can be used on other fault scarps to infer structural development. [ABSTRACT FROM AUTHOR]

Contiguous continental high-pressure terranes in orogens offer insight into deep recycling and transformation processes that occur in subduction zones. These remain poorly understood, and currently debated ideas need testing. The approach we chose is to investigate, in detail, the record in suitable rock samples that preserve textures and robust mineral assemblages that withstood overprinting during exhumation. We document complex garnet zoning in eclogitic mica schists from the Sesia Zone (western Italian Alps). These retain evidence of two orogenic cycles and provide detailed insight into resorption, growth, and diffusion processes induced by fluid pulses in high-pressure conditions. We analysed local textures and garnet compositional patterns, which turned out remarkably complex. By combining these with thermodynamic modelling, we could unravel and quantify repeated fluid-rock interaction processes. Garnet shows low-Ca porphyroclastic cores that were stable under (Permian) granulite facies conditions. The series of rims that surround these cores provide insight into the subsequent evolution: the first garnet rim that surrounds the pre-Alpine granulite facies core in one sample indicates that pre-Alpine amphibolite facies metamorphism followed the granulite facies event. In all samples documented, cores show lobate edges and preserve inner fractures, which are sealed by high-Ca garnet that reflects high-pressure Alpine conditions. These observations suggest that during early stages of subduction, before hydration of the granulites, brittle failure of garnet occurred, indicating high strain rates that may be due to seismic failure. Several Alpine rims show conspicuous textures indicative of interaction with hydrous fluid: (a) resorption-dominated textures produced lobate edges, at the expense of the outer part of the granulite core; (b) peninsulas and atoll garnet are the result of replacement reactions; and (c) spatially limited resorption and enhanced transport of elements due to the fluid phase are evident along brittle fractures and in their immediate proximity. Thermodynamic modelling shows that all of these Alpine rims formed under eclogite facies conditions. Structurally controlled samples allow these fluid-garnet interaction phenomena to be traced across a portion of the Sesia Zone, with a general decrease in fluid-garnet interaction observed towards the external, structurally lower parts of the terrane. Replacement of the Permian HT assemblages by hydrate-rich Alpine assemblages can reach nearly 100% of the rock volume. Since we found no clear relationship between discrete deformation structures (e.g. shear zones) observed in the field and the fluid pulses that triggered the transformation to eclogite facies assemblages, we conclude that disperse fluid flow was responsible for the hydration. [ABSTRACT FROM AUTHOR]

Rhythmic seismicity associated with spine extrusion is a well-documented phenomenon at a number of dome-forming volcanic systems. At Unzen volcano, Japan, a 4-year dome-forming eruption concluded with the emplacement of a spine from October 1994 to February 1995, offering a valuable opportunity to further investigate seismogenic processes at dome-forming volcanoes. Using continuous data recorded at a seismic station located close to the dome, this study explores trends in the seismic activity during the extrusion of the spine. We identify a total of 12 208 volcano-seismic events in the period between October 1994 and February 1995. Hourly event counts indicate cyclic activity with periods of ~ 40 to ~ 100 h, attributed to pulsatory ascent defined by strain localisation and faulting at the conduit margins. Waveform correlation revealed two strong clusters (a.k.a. multiplets, families) which are attributed to fracturing along the margins of the shallow, ascending spine. Further analysis indicates variable seismic velocities during the spine extrusion as well as migration of the cluster sources along the spine margins. Our interpretation of the results from seismic data analyses is supported by previously published field and experimental observations, suggesting that the spine was extruded along an inclined conduit with brittle and ductile deformation occurring along the margins. We infer that changes in stress conditions acting on the upper and lower spine margins led to deepening and shallowing of the faulting sources, respectively. We demonstrate that the combination of geophysical, field and experimental evidence can help improve physical models of shallow conduit processes. [ABSTRACT FROM AUTHOR]

The article presents a study which investigates the impact of seamount subduction and aseismic ridge on the overriding plate deformation, and seismic coupling and earthquake rupture behavior in various localities. The study used a computational recommendation technology, an information filtering system technique, and the Internet. Results show that smaller volcanic chains and ridges observed to have fragile heterogeneous internal structure.

The article presents a study which examines the potential of temperature-sensitive analogue materials for three-dimensional (3-D) thermo-mechanical modeling of plate tectonic processes. It describes the advantages of the experimental set-up to investigate the 3-D solid-mechanics collisions of lithospheric plates in various settings. It demonstrates the experiment and suggests that the analogue model enables estimation of the stress conditions along the interplate zone and the plate boundary.

The article presents a study which examines the atmospheric effect of electromagnetic charges occurring in the surface layer through a model of self-generated seismo-electromagnetic oscillation of the local segment of the lithosphere-atmosphere-ionosphere (LAI) system. It estimates the frequency and characteristics of a seismic activity and the relation of the LAI system. It notes that the model explains the precursor effects of ruptures and suggests that it may likely to forecast earthquakes.

The article presents a study which examines the seismic activity provided by the reflection seismic data over the end-glacial Burträsk fault in Skellefteå, Sweden. It explores the correlation of the seismic activity between the Burträsk fault and the locations of the fault systems. It shows that no large displacements were identified along the fault zone as the glacial ice receded. It suggests that fault constraints will help determine the factors activating the fault.

The article presents a study which explores the potentials and limitations of time-reversal wave field imaging of finite seismic sources. It characterizes the seismic sources with time-reversed fields in the process of fault rupture localization, slip distribution estimation, and source complexity inversions. It shows that time-reversal imaging does not improve the frequency band modifications and suggests that such limitation starts at the point where the source stops being point-localised.

The article presents a study which investigates the potential of the Open Plot Project open-source software for three-dimensional (3-D) structural data analysis in geological modeling processes. It describes the funcationalities of the Open Plot and the internal data workflow in providing georeferenced images that represent surface faults and seismic zones. It suggests that the Open Plot addresses the limitations in data selection, management, and subsequent analysis integration.