Orientation Of Stress Field Research Articles

The Ciclón - Exploradora epithermal polymetallic ore deposit: An Eocene Cordilleran-type of mineralization in the Domeyko Cordillera, Northern Chile

The Ciclón-Exploradora Project is located in the Exploradora Porphyry Copper District, which corresponds to an old artisanal mining and exploration area 65 km northeast from El Salvador Copper Mine, between 3200 and 4000 m.a.s.l. and inserted within the Eocene-Oligocene metallogenic porphyry belt and Domeyko Cordillera. Exploitation has been discontinuous since the late mid-19th century. It first started with production at the Copper Exploradora Mine, and later for Au–Ag at the Ciclón and San Carlos mines. Nevertheless, after decades of work and exploration the district has great potential for new discoveries.Epithermal polymetallic (Cu–Ag–Zn–Pb–Au) mineralization took place related to the evolution of the Exploradora Plutonic and Porphyry Complex. This Eocene Complex is a composite of dioritic to monzodioritic magmatic intrusions including several intermediate composition Cu–Au and Cu–Mo porphyry centers, mostly emplaced in NNE to NE-oriented corridors. The hosted rocks consists of folded Triassic and Jurassic calcareous sedimentary units partially altered by a kilometric metasomatic envelope that contains skarn mineralization.The complex was mainly emplaced syntectonically along the Sierra Castillo Fault System, which is part of the margin-parallel Domeyko Fault System. These faults together with NE and NW striking transversal faults were active during the Eocene, playing a key role in the emplacement of magmatism and mineralization at the district.The main magmatic phases of the Exploradora Complex record different levels of alteration exposure, mineralization and geochronological ages suggest at least three discrete and continuous episodes of magmatism during ∼4 m.y. (1) An early composite of dioritic to monzodioritic intrusions (38-36 Ma), generally associated with coarse Ca (K) alteration and spatially close to Zn or Fe exoskarn mineralization. (2) An early to inter mineral hypabyssal rocks (36-35 Ma) composed of dacitic to rhyolitic porphyries and aplitic rocks, mostly affected by k-feldspar, propylitic, quartz and argillic alteration, locally with secondary biotite and Cu anomalies. Include quartz monzodiorite porphyries exposing Cu–Au mineralization associated with potassic alteration (biotite-magnetite plus albite) and related quartz-veinlets. (3) CODELCO's Exploradora Cu–Mo mineralization (35-34 Ma) associated with late dacitic and dioritic porphyries, HS-Au mineralization veinlets and a post-mineral diatreme.Polymetallic mineralization occurs at the western part of the district, hosted mainly in sedimentary and intrusive rocks and partially at major intrusion contacts. The Ciclón-Exploradora polymetallic deposit is a 5 km-long NS- and NNW-striking vein-breccia-fault, with mainly sinistral movements, it includes carbonate replacement deposits and dissolution breccia bodies nearby.The tectonic stress field orientation related to the vein-breccia-fault mineralization records subhorizontal NW-trending σ1, SW-trending σ3, and subvertical σ2. Polymetallic hypogene mineralization shows a consistent evolution from low sulfidation conditions to intermediate sulfidation assemblages, spatially related to an early skarn stage characterized by the presence of magnetite, pyrrhotite, Fe-rich sphalerite, developed in calcareous sedimentary rocks. Skarn zonation includes garnet-pyroxene-wollastonite-scapolite minerals. Epithermal polymetallic mineralization shows four main stages. The first and second stages consist of pervasive chlorite followed by chlorite-pyrite veinlets that may or may not contain quartz, and massive pyrite mineralization bodies generally forming euhedral crystals. The second phase ends with a Cu-dominant stage given by chalcopyrite-quartz (±chlorite - sericite – pyrite) and lower contents of Ag-bearing sulfosalts. The third stage is recognized as a polymetallic event characterized by multiple veinlets emplaced both in intrusive and sedimentary rocks where it also forms replacement ore bodies along lithological contacts. It consists of Cd-rich sphalerite, galena, Fe–Mn–Zn carbonates, and minor amounts of chalcopyrite, tennantite-tetrahedrite, and other Ag sulfosalt minerals. Gold is present as native and spectrum minerals with quartz. Main alteration minerals are quartz-kaolinite and minor amounts of sericite. The last stage consists of calcite ± arsenopyrite veinlets crosscutting the previous mineralization. This paragenetic evolution is similar to the one observed in several epithermal porphyry-related polymetallic “Cordilleran” deposits (e.g., Cerro de Pasco, Perú).Crosscutting relationship and exposure of epithermal mineralization, porphyry alteration and magmatism suggests that, in a short period (38–34 Ma), the whole district was exhumed probably several kilometers (3–6 km), related to a compressional to transpressional deformation regime.

S Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

SUMMARY The orientation of SHmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed SHmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

The February 2018 seismic swarm in São Miguel, Azores

The island of São Miguel is among the most seismically active areas of the Azores archipelago. This work focuses on the most significant recent swarm, which occurred on February 2018. We set up an automated procedure to process continuous full seismic waveform data from local stations to generate high-quality earthquake information on the volcano unrest episode. First, we applied an automated detector software, next we located the detected events and then classified the earthquakes based on their waveform similarity, identifying three families of seismic events. We then extended the catalog by template matching. Finally, we computed moment tensors to investigate the source mechanisms of the largest earthquakes. Our results image the ∼2-week swarm evolution. The activity started with a precursory phase with low rate and low magnitude (ML < 2.0) seismicity and the activation of a deeper structure (∼10-15 km). After ∼1 week, a new earthquake family emerged at shallower depths (∼8–12 km) reaching magnitudes up to ML 3.4. Finally, a third slightly shallower family was activated. Moment tensors show mostly normal faulting mechanisms, striking ∼NW-SE, compatible with the orientation of the regional stress field. A surface deformation transient was recorded by geodetic stations, starting with the swarm, and continuing over the following ∼17 months, corresponding to either inflation or extension around the swarm region. The prolonged surface deformation implies a process that was initiated during the swarm and subsequently accommodated mostly aseismically. We interpret the seismicity observed at the early stage of deformation as indicating episodic fluid injection through the crust, related to the local hydrothermal or magmatic systems. We conclude that the Fogo-Congro region continues to be seismo-volcanically active, with both seismic and aseismic deformation observed and requiring close multidisciplinary monitoring. The proposed methology based on the automated analysis of continuous waveform data provides high-quality imaging of the spatio-temporal evolution of seismicity, which can be used elsewhere in the operational monitoring of seismo-volcanic crises to gain insight into the ongoing deformation processes, improve hazard assessment and help in the development of effective mitigation strategies.

In Situ Stress Determination Based on Acoustic Image Logs and Borehole Measurements in the In-Adaoui and Bourarhat Hydrocarbon Fields, Eastern Algeria

The study of in situ stress from image logs is a key factor for understanding regional stresses and the exploitation of hydrocarbon resources. This work presents a comprehensive geomechanical analysis of two eastern Algerian hydrocarbon fields to infer the magnitudes of principal stress components and stress field orientation. Acoustic image logs and borehole measurements were used in this research to aid our understanding of regional stress and field development. The studied In-Adaoui and Bourarhat fields encompass a combined thickness of 3050 m of Paleozoic and Mesozoic stratigraphy, with the primary reservoir facies in the Ordovician interval. The Ordovician sandstone reservoir interval indicates an average Poisson’s ratio (v) of 0.3, 100–150 MPa UCS, and 27–52 GPa Young’s modulus (E). Direct formation pressure measurements indicate that the sandstone reservoir is in a hydrostatic pore pressure regime. Density-derived vertical stress had a 1.1 PSI/feet gradient. Minimum horizontal stress modeled from both Poisson’s ratio and an effective stress ratio-based approach yielded an average 0.82 PSI/feet gradient, as validated with the leak-off test data. Drilling-induced tensile fractures (DITF) and compressive failures, i.e., breakouts (BO), were identified from acoustic image logs. On the basis of the DITF criterion, the maximum horizontal stress gradient was found to be 1.57–1.71 PSI/feet, while the BO width-derived gradient was 1.27–1.37 PSI/feet. Relative stress magnitudes indicate a strike-slip stress regime. A mean SHMax orientation of N130°E (NW-SE) was interpreted from the wellbore failures, classified as B-quality stress indicators following the World Stress Map (WSM) ranking scheme. The inferred stress magnitude and orientation were in agreement with the regional trend of the western Mediterranean region and provide a basis for field development and hydraulic fracturing in the low-permeable reservoir. On the basis of the geomechanical assessments, drilling and reservoir development strategies are discussed, and optimization opportunities are identified.

Magnetization Changes Induced by Stress Noncoaxial with the Magnetic Field in a Low-Carbon Steel

Ferromagnetic materials are widely used in the manufacturing of key parts of energy equipment, due to their good mechanical properties, such as in nuclear power and pipes. Mechanical stress exists inside of these key parts during operation. Stress can be estimated indirectly by nondestructive testing methods that measure the magnetic flux leakage signals on the surface of the structure, which is of great importance for ensuring the safety of the equipment. However, the physical mechanism of the stress and magnetic field in the magnetization of ferromagnetic materials is still unclear, leading to limited applications of the technique in practice. In this paper, magnetization tests were carried out to investigate magnetization changes under the coupling effect of stress and a noncoaxial magnetic field. Two identical Q195 low-carbon steel specimens were tested. Specimen 1 was subjected to magnetic field values successively increasing from 0 A/m to 6000 A/m under constant uniaxial tension at different angles θ between the field and stress axis. Specimen 2 was subjected to the same magnetic field under different levels of stress at an angle of 0°. The surface magnetic induction B of the specimens was measured and analyzed at each angle of stress–field orientation and at different levels of stress. It was found that there was a difference in the direction between the B and the magnetic field H at different angles θ. The magnetization curves correlated to the angle θ and the stress levels. The behavior of the derived maximum differential permeability and maximum magnetic induction could be used for the nondestructive evaluation of stress magnitude and direction in materials already in service.

Influence of regional background stress fields on the spontaneous rupture of the major faults around Xiluodu dam, China

Simulations of the spontaneous rupture of potential earthquakes in the vicinity of reservoir dams can provide accurate parameters for seismic resilience assessment, which is essential for improving the seismic performance of reservoir dams. In simulations of potential spontaneous ruptures, fault geometry, regional stress fields, constitutive parameters of the fault friction law, and many other factors control the slip rate, morphology, and dislocation of the rupture, thereby affecting the simulated ground motion parameters. The focus of this study was to elucidate the effects of the background stress field on the nucleation and propagation of spontaneous ruptures based on the factors influencing potential M > 7 earthquake events on the Leibo Middle Fault (LBMF) and the Mabian-Yanjing Fault (MB-YJF) in the Xiluodu dam (XLD) region. Our simulation results show that the magnitude of the regional background stress field plays a decisive role in whether a destructive earthquake exceeding the critical magnitude will occur. We found that the direction and magnitude of the regional stress significantly affect the range of rupture propagation on the fault plane, and fault geometry affects the spatial distribution of the rupture range. Under the same regional stress field magnitude and orientation, a more destructive, high-magnitude earthquake is more likely to occur on the LBMF than on the MB-YJF.

Geologic Evidence of Lithostatic Pore Fluid Pressures at the Base of the Subduction Seismogenic Zone

AbstractGeophysical and rheological work below the subduction seismogenic zone suggest that elevated pore fluid pressures (Pfs) may facilitate slow slip and tremor. However, quantification of Pf and its variation in this environment remain indirect and incomplete. We examine brittle‐frictional deformation within metasedimentary rocks from an exhumed subduction zone. Here, penetrative foliation formed under viscous deformation conditions by pressure‐solution creep at 0.9 GPa, 300–350°C, and differential stresses <35 MPa. Three sets of crack‐seal veins, precipitated at peak subduction temperatures, cross‐cut this foliation. Mechanical inversions of these veins indicate they formed within a single stress field orientation with foliation, but at different magnitudes of differential stress. These deformation structures reveal near‐lithostatic Pf that fluctuated in time by perhaps only ∼4 MPa, resulting in cyclic mode I and shear fracturing and viscous flow, and support slow slip and tremor mechanisms that require or are enhanced by high Pf and low differential stress.

Earthquakes Indicated Stress Field Change During the 2006 Unrest of Augustine Volcano, Alaska

AbstractTo examine controls on the local stress field at Augustine Volcano, Alaska, before its 2006 eruption, we calculated fault plane solutions for volcano‐tectonic earthquakes from 2002 to 2006. The P‐axis orientation was first aligned to the regional maximum compression (NW) and then rotated by about 90° (perpendicular to the dike alignment) after the onset of surface deformation in mid‐August 2005. Using 3D finite element models, we systematically evaluated the effects of tectonic stresses, volcanic edifice densities, and dike overpressures on the local stress field orientation. Combining data and models to generate “phase diagrams” of different stress controls by these competing effects, we argue that moderate tectonic stress of 2–3 MPa at 600 m above sea level slightly exceeded the edifice loading before the precursory deformation and was then overprinted by a local stress field from dike opening with an overpressure of ∼15 MPa.

Spatial and temporal stress field changes in the focal area of the 2016 Kaikōura earthquake, New Zealand: A multi-fault process interpretation

To understand the stress controls on the occurrence of a multi-fault rupture, we estimated the crustal stress between April 2013 to December 2018, i.e., before and after the Mw7.8 Kaikōura earthquake that occurred in New Zealand on 13 November 2016. We used both the focal mechanism solutions from the temporary seismic networks and the GeoNet moment tensor solutions and selected the solutions that differed significantly from the mainshock fault planes and rakes. Then, we performed stress tensor inversions for the selected focal mechanism solutions. Using the stress tensor inversion results, we also calculated the slip tendency. Prior to the Kaikōura earthquake, the stress regime was the strike-slip type, and the maximum eigenvalue of the stress tensor (σ1) was oriented WNW–ESE. The stress field orientation did not change significantly after the earthquake. This suggests that the stress change during the Kaikōura earthquake was too small to alter the stress orientations, implying that there may have been large differential stress prior to the Kaikōura earthquake. However, the average stress ratio in different clusters changed in two different patterns after the earthquake, suggesting possible changes in the magnitude of different components of the stress tensor, or of pore pressure in different regions. A high slip tendency was observed at the hypocentre, while a low slip tendency was observed at the northern end of the Kaikōura earthquake faults. This may suggest that the stress orientation and the stress ratio controlled the initiation and the end of the multi-fault rupture. These results corroborate previous fault propagation models.

Seasonal terrestrial water load modulation of seismicity at the southeastern margin of the Tibetan Plateau constrained by GNSS and GRACE data

SUMMARY The seasonal terrestrial water load modulation of seismicity is an important phenomenon to understand the mechanism behind earthquake triggering and nucleation. The presence of high-level seismicity and large seasonal water load amplitudes at the southeastern margin of the Tibetan Plateau provides a natural experimental environment for studying the modulation mechanism. The spatiotemporal distribution of the water load was inverted using Global Navigation Satellite System (GNSS) and Gravity Recovery and Climate Experiment (GRACE) data, and an earthquake catalogue (M ≥ 2.5) was declustered to obtain the background seismicity using the Epidemic Type Aftershock Sequences (ETAS) model. The Multichannel Singular-Spectrum Analysis (M-SSA) is adopted to the time-series of monthly averaged terrestrial water load and background seismicity rates, and the results show 1- and 2-yr periodicities in the seismicity rates and water load. The 1-yr periodicity in the seismicity rate is correlated with the rate of change of the monthly averaged water load. To evaluate the seasonal principal stress perturbations on the tectonic background stress orientations, the stress changes caused by the seasonal water load changes are projected onto the tectonic background stress field orientations constraining by 8 yr of earthquake focal mechanism solutions. The results show that the largest change of the seasonal principal stress perturbations is about 16 kPa. The number of excess earthquakes is evaluated with the background seismicity rates for discrete stress intervals. The results indicate a ∼10 per cent increase in the seismicity rates that correlate with the rates of the minimum and mean principal stress perturbations. The results above can be explained by the model of harmonic stress perturbations on rate-and-state fault. Based on this model, the nucleation period of the seasonal seismicity should be less than 1 yr at the southeastern margin of the Tibetan Plateau.

Current stress field and its relationship to tectonism in a coal mining district, central China, for underground coal energy exploration

Adopting 47 sets of in situ stress measurements taken at depths ranging from 417 to 1123 m in a coal mining district in central China by the overcoring method, we characterize the magnitude and orientation of the current stress field and its relationship to tectonism in this district to determine the controls on in situ stress in underground coal energy exploration and seismogeology. The measurements reveal that this district is in a relatively high horizontal stress environment. The stress regimes are generally favorable for strike-slip faulting and reverse faulting and are characterized by σH >σv >σh and σH >σh >σv, revealing that the regional stress field is absolutely dominated by horizontal tectonic stress and reflects a typical tectonic stress field type. The maximum horizontal principal stress is predominantly oriented between the ENE–WSW and WNW–ESE directions and averages nearly E–W,​ which is closely related to the tectonic plate movement and crustal deformation and agrees well with the direction derived from focal mechanism solutions and other stress indicators in and around the region. The current stress field in this district generally inherited the characteristics of the tectonic stress field in the third stage and partially retained the traces generated by the tectonic stress field in the first and second stages. In addition, the frictional parameter μm, which is calculated using the derived stress data, is used as a strong predictor to reflect the capability of the level of crustal stress accumulation. This region is in a state of critical or subcritical equilibrium, and regional faults at shallow depths have a medium to slightly high frictional strength, indicating that some well-oriented faults may theoretically approach or reach frictional equilibrium and that frictional sliding may occur.

Micro-fractures in the Ngubevu Gold and Base Metal Mines (Natal Thrust Front, Tugela Terrane, South Africa): Are they consistent with macro-scale deformation, and what do they tell us?

AbstractIn the Natal Thrust Front (NTF) of the Tugela Terrane, there is a lack of multi-scale characterisation of structural controls on shear-zone hosted precious and base metal mineralisation. This paper contributes to the current understanding of the structural controls on shear-zone-hosted precious and base metal mineralisation in the NTF with insight derived from the analyses of micro-fractures at and around the known shear-zone-hosted lode gold and base metal mineralisation in the Ngubevu area of the NTF. Based on the notion that the spatial distribution of stresses that control geological structures is fractal (i.e., possessing scale-invariance or self-similarity properties) across a range of geographic scales, this study (a) determines whether the stresses inferred from micro-fractures in the gold mines in Ngubevu are consistent with the stresses associated with regional- to local-scale geological structures in the NTF, and (b) interprets the mineralisation controls implied by these inferred and known stresses. From rose diagrams showing trends of micro-fractures in oriented vein material samples from Ngubevu, the relative micro-scale stress field per sample based on the maximum horizontal stress (SHmax) and the minimum horizontal stress (SHmin) were inferred. These were then compared to the local- to regional-scale stress field orientations derived from existing studies in the NTF.The north-northeast–south-southwest orientations of SHmax of micro-fractures in samples from Ngubevu are consistent with the north-northwest–south-southeast orientations of SHmax of regional- to district-scale thrust faults in the NTF of the Tugela Terrane. Likewise, the west-northwest–east-southeast orientations of micro-fractures in samples from the eastern part of Ngubevu are consistent with the general east–west orientations of local-scale S1 foliation. These orientations are associated with a D1 thrusting event, which was due to north to north-northeast-directed regional thrust tectonics. The clockwise change from the north-northwest–south-southeast orientations of SHmax of regional- to district-scale thrust faults to the north-northeast–south-southwest orientations of SHmax of micro-fractures in vein material attests to an earlier posited notion of clockwise rotation of regional- to district-scale thrust vergence. These suggest that V1 veins in Ngubevu were likely formed late in the D1 event. The orientations of micro-fractures in samples, mostly from the western part of Ngubevu, depict a strike-slip regime, which is consistent with D2 sinistral transpression that affected the earlier-formed east–west-trending structures. This suggests that some V1 veins in Ngubevu, particularly those in its western part, were affected by the D2 sinistral transpression (e.g., by remobilisation of gold in some of the V1 veins). The novel findings in this study are crucial to geological mapping, recognition of mineralisation controls and mineral prospecting in poorly-explored regions where the geology is permissive (i.e., represents processes that are favourable) for the formation of certain type(s) of mineralisation.

Rotation of the Stress Tensor in a Westerly Granite Sample During the Triaxial Compression Test

We simulated the spatiotemporal modelling of 3D stress and strain distributions during the triaxial compression laboratory test on a westerly granite sample using finite-difference numerical modelling implemented with FLAC3D software. The modelling was performed using a ubiquitous joint constitutive law with strain softening. The applied procedure is capable of reproducing the macroscopic stress and strain evolution in the sample during triaxial deformation until a failure process occurs. In addition, we calculated focal mechanisms of acoustic emission (AE) events and resolved local stress field orientations. This detailed stress information was compared with that from numerical modelling. The comparison was made based on the 3D rotation angle between the cardinal axes of the two stress tensors. To infer the differences in rotation, we applied ANOVA. We identified the two time levels as the plastic deformation phase and the after-failure phase. Additionally, we introduced the bin factor, which describes the location of the rotation scores in the rock sample. The p values of the test statistics F for the bin and phase effects are statistically significant. However, the interaction between them is insignificant. We can, therefore, conclude that there was a significant difference in the time between the rotation means in the particular bins, and we ran post hoc tests to obtain more information where the differences between the groups lie. The largest rotation of the stress field provided by the focal mechanisms of AE events from the numerically calculated stress field is observed in the edge bins, which do not frame the damage zone of the sample.

Shear Capacity Model with Variable Orientation of Concrete Stress Field for RC Beams Strengthened by FRP with Different Inclinations

A design-oriented analytical model able to evaluate the shear capacity of reinforced concrete (RC) beams strengthened with fiber-reinforced polymer (FRP) sheets or strips oriented in any direction is proposed. The formulation of the model is based on the variable-inclination stress-field approach, aiming to extend the provisions of current European standards to beams strengthened in shear using FRP. The main novelty of the model lies in taking into account the possible different inclination of steel stirrup and FRP reinforcement in determining the orientation of a compressed concrete stress field, and in shear strength evaluation, overcoming the approximation of the known models with variable inclination of the concrete strut in the assessment of concrete strut capacity, in which the value that has to be assigned to the shear reinforcement direction is not defined, that is, either that of the steel stirrup or the external FRP reinforcement. The proposed model is able to take into account different steel stirrup and external FRP shear reinforcement orientation in assessing the reduction of the steel transverse reinforcement efficiency due to the brittle failure of the composite and also as a function of the effective composite to yielding steel strain ratio. Moreover, regarding the former aspect, a simplified approximate procedure is proposed for solving the drawbacks related to verifying compressed concrete strength in the suggested method of application of code models for RC beams strengthened by means of FRP reinforcement inclined with a different slope from the pre-existing steel stirrup. Complete and U-shaped schemes are considered. The effectiveness of the proposed model adopting different relations for assessment of the FRP effective strains proposed in the literature is investigated, differentiating them by shape of the cross section and by the possible presence of fiber-anchoring devices. The shear capacity predicted by the model and those obtained using international codes and literature models are compared against the experimental results, proving that the proposed model is the most effective in predicting the shear strength when considering specimens having steel stirrups and FRP shear reinforcement arranged with different inclinations.

Earthquake Source Mechanisms and Stress Field Variations Associated With Wastewater‐Induced Seismicity in Southern Kansas, USA

AbstractThe strong increase of seismicity rates in the contiguous USA over the last 10 years is linked to the injection of huge amounts of wastewater from oil and gas production in unconventional hydrocarbon reservoirs. We calculated 549 moment tensors of induced earthquakes (MW ≤ 4.9) in southern Kansas to study their source mechanisms and their relation to injection activity. Seventeen percent of the events analyzed contained significant volumetric (ISO%) components, and these events mostly occurred near the two largest local earthquakes during the 4 months of largest active wastewater disposal. Mapping the local stress field, we determined that most of the region lies within a transtensional stress regime, with a maximum horizontal stress σHmax trending N75°E. In the epicentral area of the MW 4.9 Milan earthquake, the σHmax trend is rotated to about S80°E. Locally, two areas display a change in the stress field orientation with depth, from transtensional above 5.5 km depth to strike slip deeper in the basement. Relating the resolved fault geometries to the obtained local stress field orientation, we find that most of the activated fault planes were optimally oriented to the current stress field and thus small stress perturbations caused by the water injection could lead to failure.

Influence of preexisting structures on the variation of neotectonic paleostress orientation in an area of the Southeastern Brazilian Margin

Abstract The Southeastern Brazilian Margin presents a NE-SW-striking structural framework, which is characteristic of the Neoproterozoic Ribeira Belt. Other important structural trends along the area are associated with E-W and NW-SE-striking structures, being most related to the Mesozoic-Cenozoic tectonic evolution. A succession of Cenozoic tectonic events has been described by many authors in the Southeastern Brazilian. Among those Cenozoic tectonic events, a Neogene to Quaternary E-W dextral strike-slip regime (EW-DT) shows a wide distribution, presenting some variations in the stress field orientation. This research investigates the influence of the preexisting structural framework on the variation of the paleostress field orientation associated with the EW-DT tectonic event in a selected area in the onshore continental margin of the Southeastern Brazil. The study was carried out in three main stages: lineament analysis, paleostress field analysis, and numerical mechanical modeling. The paleostress field related to the EW-DT tectonic event in the studied area presents: NNW-SSE-striking maximum horizontal stress, and counterclockwise rotation influenced by the structural framework. Regarding this, the regional variation of the paleostress field is influenced by a preferential reactivation of NE-SW-striking structures while sets of approximately E-W-striking structures are the main local control of counterclockwise rotation of the stress field.

Aeromagnetic interpretation of basement structure and architecture of the Dahomey Basin, Southwestern Nigeria

ABSTRACT Significant yet untapped resources still abound in the Nigerian sector of the Dahomey Basin. Although the presence of an extensive Cretaceous Petroleum System has been confirmed following recent discoveries in Aje and Ogo Fields offshore Lagos, exploration outputs in the Dahomey Basin has so far not been encouraging. Proper understanding of the basement architectural framework and controls on tectonic development remains key to unlock the unrealised potentials in the basin. Hence, a geophysical interpretation of the basement structure and architecture of the Dahomey Basin southwestern Nigeria has been carried out in this study. Various edge enhancement techniques were applied to the high-resolution residual magnetic intensity (HRRMI) grid of the area. This includes first vertical derivative (FVD), total horizontal derivative (THDR), tilt derivative (TDR) and total horizontal derivative of upward continuation (10 km). Determination of the depth to magnetic sources and sedimentary thicknesses in the study area were achieved using standard Euler deconvolution and source parameter imaging (SPI) techniques, with depth range of 4.5–6.3 km attained in the two identified sub-basins located offshore of the study area. Lineament analysis gave insights on the tectonic trends and stress-field orientation in the basin with major trends in the NNE-SSW, NE-SW, NW-SE, and WNW-ESE directions. 2D forward modelling of some selected profiles was employed to characterise the basement pattern and architecture, which depicted a horst-graben architecture. The basement structure and architecture have a major control on the distribution of sub-basins, petroleum systems elements and trap styles in the basin. The study demonstrates the robust application of high-resolution aeromagnetic data in basin-wide mapping of regional subsurface geological features, basement architecture and determination of sedimentary thickness in a frontier basin.

Correlation between crustal anisotropy and seismogenic stress field beneath Shillong–Mikir Plateau and its vicinity in North East India

A systematic study towards understanding the correlation between polarization direction of crustal anisotropy with seismogenic stress field at different locations of the Shillong-Mikir Plateau and its vicinity in North East India is attempted. We used data from a 17-station broadband seismic network. In our earlier work , crustal anisotropic parameters were determined using ANISOMAT + for the 17 seismic stations. In this study, we have estimated stress field around the stations using focal mechanism solutions (FMS). Some 215 FMS are obtained by waveform inversion. These solutions are used for stress tensor inversion to estimate stress field around each location. It is observed that polarization direction of crustal anisotropy is consistent with that of the maximum horizontal stress (Ơmax) as well as the minimum horizontal stress (Ơmin). In addition to this, two orthogonal fast polarizations in some locations are also noted. The bivariate nature of correlations helps us to understand that the major mechanisms of seismic crustal anisotropy are not only due to the regional stress, but active faults and other geological conditions play a significant role in contemporary orientation of seismic crustal anisotropy and seismogenic stress field.

Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions

AbstractUnderstanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo‐Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High‐resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first‐order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. The results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.

Stress fields of ancient seismicity recorded in the dynamic geometry of pseudotachylyte in the Outer Hebrides Fault Zone, UK

Heterogeneous sequences of exhumed fault rocks preserve a record of the long-term evolution of fault strength and deformation behaviour during prolonged tectonic activity. Along the Outer Hebrides Fault Zone (OHFZ) in NW Scotland, numerous pseudotachylytes record palaeoseismic slip events within sequences of mylonites, cataclasites and phyllonites. To date, the kinematics and controls on seismicity within the long active history of the OHFZ have been poorly constrained. Additional uncertainties over the relative location of a meteorite impact and possible pre-OHFZ brittle faulting also complicate interpretation of the diffuse seismic record. We present kinematic analyses of seismicity in the OHFZ, combining observations of offset markers, en echelon injection veins and injection vein geometry to reconstruct slip directions and stress fields. This new dataset indicates that a range of fault orientations, slip directions and slip senses hosted seismicity in the OHFZ. Such complexity requires several stress field orientations, in contrast with the NW–SE Caledonian compression traditionally attributed to frictional melting along the OHFZ, indicating that seismicity had a long-term presence across the fault zone. Persistence of strong frictional failure alongside the simultaneous development of weak fault rocks and phyllonitic shear zones in parts of the OHFZ has significant implications for understanding seismic hazard along mature continental faults. Supplementary material: Tables listing analysed orientation measurements plus further information and sensitivity testing of palaeostress analysis parameters are available at https://doi.org/10.6084/m9.figshare.c.5134797