Fault Dip Research Articles

Contrasting faulting styles of salt domes and volcanoes: Can unsupervised learning techniques differentiate fault styles?

An attractive feature seen on seismic data, also known as funny-looking thing (FLT), has a wide range of interpretations, from noise patterns to amplitude anomalies. An example of an FLT is the similar faulting patterns between a volcanic intrusion and a salt intrusion from the point of view of a machine learning (ML) algorithm. Oftentimes, seismic interpreters do not have a complete data set or geologic background to determine the genesis of the observed features. This can be particularly perplexing when trying to determine if an intrusion is volcanic or halokinetic in origin because they exhibit similar geomorphologies. Examining the differences in these features in the Gulf of Mexico, a well-documented salt basin, and the Taranaki Basin in New Zealand, which is igneous prone, can help to understand the differences. The analysis aims to discern geologic features based on the geometries and attributes shared by seismic data and remote sensing tools. Seismic attributes and ML techniques highlight differences and similarities between the intrusions, including the discussion of using ML techniques, such as self-organized maps (SOMs), an unsupervised ML technique, and cluster fault systems, without regard to the geologic context. The attributes used in the SOM are fault probability, fault dip azimuth, fault dip magnitude, and thin-bed detector. Fault probability is performed through a combination of convolutional neural network fault prediction and a skeletonization process. Once the faults are clustered using SOM, the visualization of fault architecture due to the existing mount (either volcano or salt dome) is done considering high fault probabilities (>75%). The methodology consists of selecting the neurons from the SOM grid corresponding to the presence of faults and combining them with fault probability and a fault dip azimuth using a crossplot. The crossplot product assists in the automatic extraction of the fault planes using: (1) a voxel representation of the fault planes and (2) fault patches representing the fault planes. Moreover, the visualization technique defined demonstrates that the crossplot product yields better-defined fault planes. With the fault system characterized, compared horizon slices using coherence, fault dip magnitude, and azimuth against remote sensing images with similar attributes. In conclusion, our methodology combines technologies to differentiate the genesis of intrusion — salt or igneous — using the fault presence and could be helpful in frontier exploration or planetary exploration. Geologic feature: Comparison of extensional faults due to salt and igneous intrusions Seismic appearance: In the symmetric salt dome, crown faulting style with a set of synthetic and antithetic faults around the salt dome; for the igneous systems, symmetric radial faults Formation: Upper and Lower Wilcox (Gulf of Mexico); Ariki Formation — Mo hakatino Formation (Taranaki) Age: Upper Miocene and Lower Pliocene (Gulf of Mexico); Mid Pliocene (Taranaki) Location: Offshore Gulf of Mexico and Offshore Northeast New Zealand Seismic data: From the Gulf of Mexico: B-78-89-LA-3D, B-69-94-LA-3D, B-54d-94-LA-3D, and B-69a-94-LA-3D. From New Zealand: Nimitz and Kora 3D Analysis tools: Structural-oriented filter, Sobel similarity, fault dip magnitude and azimuth, CNN automatic fault prediction, and SOM classification

Coupled thermal-hydro-mechanical modeling on characteristics of excavation damage zone around deep tunnels crossing a major fault

An understanding of excavation damage zone (EDZ) is critical in safe construction of underground openings, but conventional continuum-based numerical methods such as finite element method or finite difference method have limitations in simulating crack development under coupled thermal-hydro-mechanical (THM) conditions. In this study, a numerical method of predicting EDZ subjected to coupled THM conditions is developed based on discrete element method (DEM), which can simulate propagation of newly formed cracks by adding fictitious joints. A physical model of a high temperature tunnel filled with pressurized water is first used to verify the accuracy and efficiency of the developed numerical method. Based on the case study of the Baoling Mountain tunnel of the Sichuan-Tibet railway, the effects of coupled THM conditions on characteristics of EDZ are investigated. The results show that the two EDZ indices of volume of plastic zones per length and area of cracks per length increase from M-only to THM scenarios. In addition, the shear displacements of the fault in the roof and floor of the tunnel also increase from M-only to THM scenarios. The increasing EDZ size and fault displacement in THM scenario is attributed to the higher ambient rock stress and water pressure in the fictitious joints and the fault. In high temperature regions, thermal stress increases the maximum principal stresses and induces rock damage. Effect of fault dip angle on characteristics of EDZ is also discussed.

Coseismic and early postseismic deformation of the 2020 Nima Mw 6.4 earthquake, central Tibet, from InSAR and GNSS observations

The 2020 Mw6.4 Nima earthquake is one of the largest normal-faulting earthquakes recently occurring north of the Banggong suture zone in remote central Tibet, where geologic investigation of active faults is extremely limited. We analyze jointly InSAR and GNSS observations over 9 months after the Nima earthquake and calculate the coseismic and postseismic displacement. The optimal coseismic slip model suggests this event is the result of moderate-angle down-dip slip on a complex reversed “S-shape” three-segmented structure at fault junctions of the West Yibu-Chaka fault, the Heishi fault, and an unmapped blind fault, with a small component of left-lateral slip. The superposition of seismic waves from faults with different strikes and dips accounted for a large non-double-couple component in the long-period point-source solutions. The geodetic moment released by the mainshock is 6.4 × 1018 N⋅m, equivalent to Mw 6.42. Coseismic rupture concentrated at a depth of 4–15 km, with a peak slip of 1.36 m at 8.5 km depth. The cumulative afterslip moment within 9 months after the mainshock is 1 × 1018 N m, about 15.6% of that released by the mainshock coseismic slip. The afterslips contributed largely to the release of additional strain energy. In addition, shallow creep on the northern part of the blind fault, and deep uplift on the east normal fault system are promoted by stress perturbations. A significant proportion of down-dip coseismic slip spreading to more than 20 km beneath the surface, and deep up-slip afterslip have implications for the rheology of down-dip extension of the dipping faults in northern Tibet. Two obvious stress loading zones of more than 1 bar highlight seismic hazards in the region, especially in the junction between normal faults and ends of the large-size sinistral Riganpei-Co and Jiangai-Zangbu faults. It is necessary to forecast accurately by longer-term afterslip observation over timescales of years for the faults. Compared with previous studies, our results suggest a more complex subsurface fault geometry linking the normal and strike-slip faults and dynamic stress adjustment in this poorly-known region of Tibet.

Mechanism Analysis of Floor Water Inrush Based on Criteria Importance though Intercrieria Correlation

Floor water inrush has always been an important problem affecting coal mine safety. The water inrush from the deep coal seam floor is the final manifestation of the water inrush channel formed by the expansion and penetration of rock mass cracks under the combined action of a high ground temperature, high ground stress, high karst water pressure and strong mining disturbance. The fault is a structurally weak plane developed in the complete rock layer, which can be regarded as a permeable channel with good permeability. In order to explore the location of confined water from the fault to rock stratum caused by the tensile fracture of a fault during floor water inrush, floor water inrush is regarded as the result of the joint action of water pressure and mining stress. The COMSOL Multiphysics numerical simulation software was used to simulate the advancing process of the working face of the coal seam with a fault structure. Different fault dip models were established for comparative analysis. The objective weight CRITIC based on data volatility was used to analyze the water pressure and principal stress near the fault. The results show that the uplift strength of highly confined water increases with the increase in the fault dip angle, and the weight of water pressure near the fault has a certain relationship with the location of the water body entering the rock stratum from the fault.

Induced seismicity in the Groningen gas field – arrest of ruptures by fault plane irregularities

AbstractFrom dynamic rupture simulations, we reveal under which conditions a rupture in the Groningen gas field stops along fault dip or along fault strike after it starts on a fault in the reservoir. The simulations focus on the capabilities of fault plane irregularities to arrest ruptures. Such irregularities can be recognised in sandstone outcrops. Fault planes in the Groningen field, extracted from the 3D seismic data by seismic attribute extraction methods, show similar irregularities. A detailed surface of a major fault plane in the field indicates that steps and jogs of tenths of metres are possible. Although these irregularities are close to seismic resolution and could be partially artificial, we investigated their effect on rupture arrest.For typical current stresses in the Groningen field, jogs and steps of this length scale are found to be remarkably effective to stop ruptures in the reservoir. Also, a significant increase in the fault dip along fault strike can stop these ruptures but a kink in the fault under a constant fault dip not.Including non-planar fault features and pressure diffusion in the Carboniferous, the simulations in this paper follow trends of previous simulations in the literature using 2D planar faults. In particular, the horizontal stress in this formation and the strength of the Carboniferous fault zone are important for rupture propagation. If the fault would have been reactivated in the Neogene or Quaternary and poorly healed in clay-rich parts, rupture propagation into the Carboniferous can only be prevented by jogs of sufficient size and lateral continuity under the present estimate of the horizontal field stress.

Fluid-rock interaction control on fault reactivation: A review of the Montmell-Vallès Fault System, central Catalan Coastal Ranges (NE Iberia)

Structural inheritance is a key factor controlling the tectonic evolution of the central Catalan Coastal Ranges. Up to two periods of tectonic inversion (one positive during the Paleogene and the other negative during the Neogene) affected a previously well-developed Mesozoic extensional basin system and characterized the Cenozoic evolution of the area. In this scenario, tectonic fault inversion is often observed along the Montmell-Vallès Fault System. Fault reactivation shows differences along strike from NE to SW and appears decoupled from surface to depth due to its kinked-planar geometry and the change of fault dip from >60° to 30° in depth. The ability of the Mesozoic faults to be reactivated appears also influenced by changes in the mechanical properties of the inherited fault zone. Whereas the deeper and less dipping panels of the major faults are reactivated in the entire zone (contractional during the Paleogene and extensional during the late Oligocene-Neogene), the upper and highly dipping parts of the faults only show local reactivations. The observations indicate that fault dip, the indirect role of hosting lithologies (granites and siliciclastic metasediments versus carbonate rocks) on fault rocks, the indirect role of mineral precipitation and cementation product of fluid circulation, and the direct role of the mechanical properties of the resulting fault rocks (gouge versus cemented breccias) significantly control the fault reactivation. Upper crust low angle fault segments are easily reactivated during contractional deformation but not during the extensional one. Conversely, the segments with a high-angle dip are more easily reactivated during the extensional deformation but not during the contractional one. In the study area, this resulted in the formation of footwall short-cuts developed during the Paleogene compression and extensional short-cuts occurred during the Neogene extension. On the other hand, reactivation is effective in areas where granites and siliciclastic metasediments characterize the host-rock, and non-cohesive fault gouge forms the pre-existent fault core. Instead, fault reactivation appears restricted or even prevented where the host-rock includes thick carbonate successions, and the pre-existent damage zone is formed by highly cemented and cohesive breccias.

A local earthquake tomography on the EAF shows dipping fault structure

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

Quantifying Dyke-Induced Graben and Dyke Structure Using 3D Seismic Reflection Data and The Role of Interpretation Bias

During dyke intrusion, tensile stresses concentrated within the overlying rock may lead to the formation of normal faults. These faults typically form graben-bounding pairs that are sub-parallel to, and dip toward, the upper tip of their underlying dyke. Many studies use geometric properties extracted from the surface expression of such dyke-induced faults to estimate the geometry of subsurface dykes. These methods assume dyke-induced faults are planar and nucleate at the surface. However, recent seismic reflection-based investigations of the 3D structure of dyke-induced faults confirm they can be non-planar and have complex growth histories. Here, we use 3D seismic reflection surveys from offshore NW Australia to: (1) examine how the surface expression of dyke-induced faults relates to subsurface dyke geometry and depth; and (2) test whether subjective bias may influence the quantitative analyses of dyke-induced faults using seismic reflection data. We show displacement and dip vary across dyke-induced faults, supporting previous suggestions that faults nucleate between dyke upper tips and the free surface. We also find that prediction of dyke upper tip depths using graben width and area of loss methods are sensitive to fault dip variations and interpretation biases, but often still produce similar results to measured dyke depths. Both measured and predicted dyke depths vary by several hundred metres along-strike, which we relate to the preservation of dyke heads, segmentation, and/or magma density changes. Overall, we show reflection seismology provides a better understanding of the 3D structure of dyke-induced faults and their relationship to the geometry and emplacement dynamics of their causal dykes.

Spatial coupling relationship between faults and gold deposits in the Jiaodong ore concentration area and the effect of thermal doming-extension on mineralisation

As the third largest gold ore-forming province in the world, the deep metallogenic potential of and exploration strategies for Jiaodong gold deposits have attracted substantial attention. By applying three-dimensional (3D) geological modelling to these super-giant and giant gold deposits, this study reveals the 3D morphological characteristics of the ore-bearing faults and the ore-related granites, analyses the relationships between ore body locations and slope dip angle changes in the faults, and discusses the coupling relationships among gold mineralisation, ore-controlling faults, and magmatic activities. In the 3D model, the dip angles of ore-controlling faults gradually decrease downwards and alternate between steep and gentle along the fault dip direction. The ore bodies mainly occur in areas with relatively high fault-dip-angle change rates. Hence, positions with fault dip angles of 20°–50° and change rates of ≥0.91 should be considered as target areas for prospecting deep gold ore bodies. In-depth analysis suggests that the ore-hosting faults, granites, and gold mineralisation resulted from regional thermal doming and extensional movements. Gold ore bodies in the Jiaodong Peninsula occur inside or on the periphery of Mesozoic granite. Strong magmatic activities during the transition from enriched to depleted mantle provided the necessary metallogenic and thermal dynamic environment for gold mineralisation. In particular, the magmatic evolution provided abundant ore-forming materials. The detachment faults with varied dip angles provided favourable space for gold mineralisation, as the fluctuations of ore-forming fluid pressure caused by sudden changes in these dip angles served as the primary mechanism underlying the step-like gold mineralisation pattern.

Dynamic simulations of coseismic slickenlines on non-planar and rough faults

SUMMARY Knowing the directions of rupture propagation of palaeo-earthquakes is a challenging, yet important task for our understanding of earthquake physics and seismic hazard, as the rupture direction significantly influences the distribution of strong ground motion. Recent studies proposed a relationship between the direction of rupture propagation and curvature of slickenlines formed during seismic slip. The relationship was established using a global catalogue of historic surface-rupturing earthquakes and dynamic models of idealized, planar faults. At the same time, some slickenlines previously documented on geometrically complex fault segments show their convexity opposite from the simple model prediction, which we refer to as ‘abnormal convexity’. To explain such observations, we perform simulations of spontaneous earthquake ruptures on non-planar and rough faults. We find that in the case of strike-slip faults, a non-planar fault model can lead to abnormal convexity of slickenlines in places where the fault dip angle changes abruptly from the average dip of the fault. Abnormal convexity of slickenlines is produced when the initial along-dip stresses are larger than, and are opposite in direction to, the dynamic stresses imparted by the mixed-mode rupture. Such results are also confirmed in our rupture simulations with a rough fault. Our results also show that the parameter space for which abnormal convexity of slickenlines occurs near Earth’s surface is narrow, especially when the fault strength and initial shear stresses increase with depth. Nevertheless, slickenlines on geometrically complex faults need to be carefully interpreted and investigation of rupture direction using curved slickenlines should focus on structurally simple parts of faults.

Thrust faulting on Venus: Tectonic modeling of the Vedma Dorsa Ridge Belt

Radar images from the Magellan mission reveal compressional ridge belts transecting several low-lying plains on Venus, including Atalanta, Vinmara, Lavinia, and Rusalka/Llorona/Vellamo Planitiae. Vedma Dorsa is an ∼1800 km long ridge belt with 0.5–1 km of relief in Llorona/Vellamo Planitiae. We analyze sections of the ridge belt that have an asymmetric topographic cross section, indicative of deep thrust faulting and shallow folding as the dominant formation mechanisms. We perform elastic dislocation and Monte Carlo tectonic modeling of thrust faulting at several locations along Vedma Dorsa to constrain the fault displacement (D), depth of faulting (z), and fault dip (θ). We find that the models are consistent with blind thrust faulting and folding as the dominant mechanism for producing Vedma, with preferred fault parameters of D ∼ 1–2 km, z ∼ 10–25 km, and θ ∼ 25–30°. The results support an upper fault tip >3 km below the surface (blind fault), and a change in dip at depth. Our preferred models and the additional fault complexities supported at each site suggest that this section of Vedma is an intermediate structure between lobate scraps and fold and thrust belts, incorporating morphometric aspects from both classifications. Observations of terrestrial fault depths and the best-fitting depth of faulting from our models lead to an estimated heat flux of 16–39 mW m−2 at the time of Vedma Dorsa's formation, consistent with formation over a cold, downwelling region of the mantle.

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

AbstractThe 2008 Mw 7.9 Wenchuan earthquake is the most disastrous seismic event in China since 1976. Both field and seismological investigations suggest a multi‐stage coseismic rupture with most damage associated with the Yingxiu‐Beichuan Fault (YBF) of spatially variable fault strike and dip angles. To investigate the effect of fault geometric complexity on coseismic rupture and paleoseismic pattern on the YBF, we perform earthquake sequence modeling on 3D fault geometry in the framework of rate‐and‐state friction law. Our model produces a long‐term earthquake sequence with quasi‐regular recurrences of whole‐fault ruptures and segmented ruptures. The along‐strike rupture segmentation, earthquake propagation speed, and slip rate are mainly controlled by along‐strike variations of fault dip and strike angles. Particularly, the YBF can be divided into two segments: the southern segment featured by large events (Mw > 8.0) and the northern segment with smaller events (Mw < 7.5), with recurrence intervals primarily determined by the tectonic loading rate. In a whole‐rupture event, the rupture speed mainly correlates with the fault dip angle; a smaller dip angle results in a wider seismogenic zone and higher rupture speeds, whereas the small‐scale variation in rupture speed is regulated by the fault strike angle. The effect of strike variation on the total coseismic slip amount is more pronounced due to the large strike angle gradient along the YBF. After varying the slip vector to reflect the northward transition from the dominant thrust‐slip to strike‐slip faulting during the Wenchuan coseismic rupture, we obtain reasonably good agreement of model simulated coseismic surface displacements with GPS observations.

Effect of Fault Slip on the Casing Passability of Horizontal Shale Oil and Gas Wells

Summary Fault slip, a secondary effect of the hydrofracture to extract shale oil and gas, deforms the casing locally and affects downhole tools going through it, hindering the exploitation plan. The casing passability is introduced to evaluate the casing passability for the downhole tool. A simulation model of the fault slip was established to analyze the effect of fault type, slip distance, and fault dip on the casing passability. Meanwhile, the simulation method adopted in this paper was verified by the large-scale physical experiments of the perforated casing subjected to uniform outer pressure. Based on the experimental results and the geometry of fault slip, an empirical formula is proposed for the relationship between casing passability and fault slip. The result shows that the sensitivity of casing passability affected by fault type is ranked from highest to lowest by strike-slip fault, normal fault, and reverse fault. This study provides a reference for determining the available diameter of downhole tools quickly and reduces the cost of exploring wells when downhole tools are obstructed.

Fault geometry and slip rates from the Nullarbor and Roe Plains of south‐central Australia: Insights into the spatial and temporal characteristics of intraplate seismicity

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

Research on local buckling failure range of X80 buried steel pipeline under oblique-reverse fault

Permanent ground deformation induced by the earthquake will cause serious damage to long-distance buried oil and gas pipelines, resulting in buckling failure or even cracking damage to pipelines. Based on the ABAQUS program, a 3-D numerical model of X80 buried steel pipeline crossing oblique-reverse fault was established, and the failure characteristics and failure range of local buckling of buried pipeline were systematically analyzed. According to the engineering parameters of the Second West-East Gas Pipeline Project (China), the influence of fault dip angle, displacement, diameter-thickness ratio, burial depth, and internal pressure on the local buckling failure range of the pipeline was discussed in detail. Subsequently, regression analysis was carried out on the numerical results of 1800 sets of the finite element cases, and a semi-empirical prediction equation (“Risk Segment” regression formula) of the local buckling failure range of X80 steel pipeline subjected to oblique-reverse fault displacement was derived. The accuracy and applicability of the formula were verified by the numerical results of 120 practical examples. The results show that the “Risk Segment” regression formula has higher accuracy and is less time-consuming, which is suitable for failure range prediction and local safety design of X80 buried steel pipeline under oblique-reverse fault.

Numerical Assessment of the Structural Damage of a Composite Lining Water Conveyance Tunnel Subjected to Reverse Fault Conditions

In this paper, the structural responses and failure characteristics of a new type of water conveyance tunnel lining structure subjected to reverse fault conditions were numerically investigated by considering multiple loads and interaction separation modes between different structural layers. This study proposes a new evaluation standard for the safety level of the damage state of the composite lining water conveyance tunnel. It also discusses the influences of fault dislocation displacement (Δf), dip angle (β), and the mechanical properties of the surrounding rock in the fault fracture zone on the water conveyance tunnel response and damage. The results indicate that the buckling failure of the steel tube under axial compression is the dominant failure mode of the composite lining structure. With increasing fault dislocation displacement, the axial compressive strain and circumferential shear strain of the composite lining are most severely damaged on the sliding plane. With decreasing fault dip angle, the axial compressive strain of the composite lining weakens, while the bending and shear strains increase. The increase in rock stiffness in the fault fracture zone reduces the damage scope but increases the composite lining structural damage severity. Overall, the numerical results of this study provide a better understanding of the failure mode and damage process of composite lining water conveyance tunnels under reverse fault conditions; therefore, this study can serve as a reference for composite lining structure disaster assessments.

Induced seismicity potential based on probabilistic geomechanics for the San Juan Basin CarbonSAFE project

Although geologic carbon sequestration projects have yet to induce — or may never induce — a damaging earthquake, experiences from other deep injection industries such as hydraulic fracturing, enhanced geothermal systems, and saltwater disposal suggest that effective quantitative seismic risk assessment is necessary for deep saline carbon capture and sequestration (CCS) projects. One such imminent CCS project is the San Juan Basin CarbonSAFE Phase III program. We use Monte Carlo probabilistic geomechanical analyses combined with observations of the geologic and operational parameters of the San Juan Basin site and suggest that this project is of low-induced seismic risk. The primary analysis is split into four sections. First, we assess the literature for faults and past seismicity, and at least five faulting scenarios are directly relevant. Second, we develop and calibrate an integrated earth model for the project site. Third, we perform Monte Carlo simulations that consider reasonable uncertainties of the geomechanical parameters. Of five tested faulting scenarios, only one presents high Coulomb failure functions, but fourth, we determine the risk to be low based on the combined lack of historical seismicity, the geologic framework of the faulting scenario, and the presence of saltwater injection at the same depth as the proposed supercritical carbon dioxide injection. The most sensitive parameters in the geomechanical calculations are the fault dip and the coefficient of friction. The least sensitive are the fault strike and the orientation of the maximum horizontal principal stress.

Structure, chronology, kinematics, and geodynamics of tectonic extension in the greater Catalina metamorphic core complex, southeastern Arizona, USA

Abstract Oligocene and early Miocene displacement on the Catalina–San Pedro detachment fault and its northern correlatives uncovered mylonitic fabrics that form the greater Catalina metamorphic core complex in southeastern Arizona, USA. Gently to moderately dipping mylonitic foliations in the complex are strongly lineated, with a lineation-azimuth average of 064–244° and dominantly top-southwest shear sense over the entire 115-km-long mylonite belt. Reconstruction of detachment fault displacement based on a variety of features indicates 40–60 km of displacement, with greater displacement in more southern areas. Widespread 26–28 Ma volcanism during early extensional basin genesis was followed by 24–26 Ma granitoid magmatism. Cooling of footwall mylonites continued until 22–24 Ma, as indicated by 40Ar/39Ar mica dates. Lower temperature thermochronometers suggest that footwall exhumation was still underway at ca. 20 Ma. Tectonic reconstruction places a variety of unmetamorphosed supracrustal units in the Tucson and Silver Bell Mountains above equivalent units that were metamorphosed and penetratively deformed in the Tortolita and Santa Catalina Mountains. This restored juxtaposition is interpreted as a consequence of older Laramide thrust burial of the deformed units, with northeast-directed thrusting occurring along the Wildhorse Mountain thrust in the Rincon Mountains and related but largely concealed thrusts to the northwest. Effective extensional exhumation of lower plate rocks resulted from a general lack of internal extension of the upper plate wedge. This is attributed to a stable sliding regime during the entire period of extension, with metamorphic core complex inflation by deep crustal flow leading to maintenance of wedge surface slope and detachment fault dip that favored stable sliding rather than internal wedge extension.

Crustal response to inflation imaged by mapping of subsurface faults, slip directions, and stress changes during the 1993–1998 unrest in Hengill, SW-Iceland

Between 1993 and 1998, the Hengill volcanic region in southwest Iceland was subject to a significant volcano-tectonic event which caused a total of 8 cm uplift and induced >90 thousand earthquakes, most of magnitude smaller than 2 ML and a small intrusion has been proposed as its cause. In this paper, we study fault activation around the uplift center to analyze whether the Coulomb stress change induced by such an intrusion would cause the observed crustal response. The locations of a subset of 12,000 micro-earthquakes are significantly improved by a double-difference relocation algorithm and the use of cross-correlation of waveforms to refine phase arrival times. The vast majority of seismic activity is scattered, but twenty-five (25) rectangular prisms, consisting of 11 to 363 earthquakes and 200 to 3000 m long, were identified and are interpreted as faults. Most of these faults are sub-vertical (within 5° from vertical) and only two faults dip at angles of 62° and 65°, respectively. Focal mechanism analysis of individual events was applied to estimate slip direction on these faults. Simulation of the Coulomb stress changes reveals that the stress changes would have been able to induce the observed fault movement on 22 of 25 faults. Discrepancies on the remaining faults are likely caused by uncertainties in the rake directions. Our results emphasize that the seismicity caused by a volume change in the crust is controlled by the interplay between the local stress field and the induced stresses.

Deformation and failure of overburden soil subjected to normal fault dislocation and its impact on tunnel

In this study, a large-scaled model test was conducted to explore the deformation and failure mechanism of tunnel with fault rupture. The strain evolution, ground fissure and displacement, earth pressure and tunnel deformation are collected in the test. Then, tunnel mechanical behavior under the fault dislocation is revealed based on the test result. The results show that a prismoid-like 3D shear zone appeared within fault fracture zone. The tunnel experiences an elongated ‘S’-shaped deformation along longitudinal axis of tunnel and the displacement changes of sectional linings within the shear zone are more significant. In transversely, the linings exhibit an oval horse-shoe deformation pattern. Three types of failure mode defined as Type Ⅰ, Type Ⅱ and Type III, corresponding to slight damage, moderate damage and severe damage, are summarized transversely in this paper. Under normal fault dislocation, the tunnel within the shear zone experienced a combination failure of bending moment and shear force. But compression failure of linings far from the fault would occur with longitudinal cracks on tunnel subjected to normal fault. Taking fault width (Wf) and tunnel span (St) into consideration, a reasonable tunnel fortification length was proposed with maximum of 2.25Wf and 3.5St under a specific fault dip of 60°. In the model test, when the fault dislocation reaches 10 mm, a slight crack begins to appear on arch foot of lining and then dislocation of joint occurs. Thus, corresponding to a prototype tunnel, a safe vertical displacement to ensure no damage to tunnel should be 34.6 cm when a segmentally designed tunnel is subjected to normal fault dislocation.