scholarly journals Heterogeneous stresses and deformation mechanisms at shallow crustal conditions, Hungaroa Fault Zone, Hikurangi Subduction Margin, New Zealand

Author(s):  
Carolyn Boulton ◽  
Marcel Mizera ◽  
Maartje Hamers ◽  
Inigo Müller ◽  
Martin Ziegler ◽  
...  

<p>The Hungaroa Fault Zone (HFZ), an inactive thrust fault along the Hikurangi Subduction Margin, accommodated large displacements (~4–10 km) at the onset of subduction in the early Miocene. Within a 40 m-wide high-strain fault core, calcareous mudstones and marls display evidence for mixed-mode viscous flow and brittle fracture, including: discrete faults; extensional veins containing stretched calcite fibers; shear veins with calcite slickenfibers; calcite foliation-boudinage structures; calcite pressure fringes; dark dissolution seams; stylolites; embayed calcite grains; and an anastomosing phyllosilicate foliation.</p><p>Multiple observations indicate a heterogeneous stress state within the fault core. Detailed optical and electron backscatter diffraction-based texture analysis of syntectonic calcite veins and isoclinally folded limestone layers within the fault core reveal that calcite grains have experienced intracrystalline plasticity and interface mobility, and local subgrain development and dynamic recrystallisation. The recrystallized grain size in two calcite veins of 6.0±3.9 µm (n=1339; 1SD; HFZ-H4-5.2m_A;) and 7.2±4.2µm (n=406; 1SD; HFZ-H4-19.9m) indicate high differential stresses (~76–134 MPa). Hydrothermal friction experiments on a foliated, calcareous mudstone yield a friction coefficient of μ≈0.35. Using this friction coefficient in the Mohr-Coulomb failure criterion yields a maximum differential stress of 55 MPa at 4 km depth, assuming a minimum principal stress equal to the vertical stress, an average sediment density of 2350 kg/m<sup>3</sup>, and hydrostatic pore fluid pressure. Interestingly, calcareous microfossils within the foliated mudstone matrix are undeformed. Moreover, calcite veins are oriented both parallel to and highly oblique to the foliation, indicating spatial and/or temporal variations in the maximum principle stress azimuth.</p><p>To further constrain HFZ deformation conditions, clumped isotope geothermometry was performed on six syntectonic calcite veins, yielding formation temperatures of 79.3±19.9°C (95% confidence interval). These temperatures are well below those at which dynamic recrystallisation of calcite is anticipated and exclude shear heating and the migration of hotter fluids as an explanation for dynamic recrystallisation of calcite at shallow crustal levels (<5 km depth).</p><p>Our results indicate that: (1) stresses are spatiotemporally heterogeneous in crustal fault zones containing mixtures of competent and incompetent minerals; (2) heterogeneous deformation mechanisms, including frictional sliding, pressure solution, dynamic recrystallization, and mixed-mode fracturing accommodate slip in shallow crustal fault zones; and (3) brittle fractures play a pivotal role in fault zone deformation by providing fluid pathways that promote fluid-enhanced recovery and dynamic recrystallisation in the deforming calcite at remarkably low temperatures. Together, field geology, microscopy, and clumped isotope geothermometry provide a powerful method for constraining the multiscale slip behavior of large-displacement fault zones.</p>

2020 ◽  
Author(s):  
Javiera Ruz ◽  
Muriel Gerbault ◽  
José Cembrano ◽  
Pablo Iturrieta ◽  
Camila Novoa Lizama ◽  
...  

<p> The Chilean margin is amongst the most active seismic and volcanic areas on Earth. It hosts active and fossil geothermal and mineralized systems of economic interest documenting significant geofluid migration through the crust. By comparing numerical models with field and geophysical data, we aim at pinning when and where fluid migration occurs through porous domains, fault zone conduits, or remains stored at depth awaiting a more appropriate stress field. <span>Dyking and volcanic activity occur within fault zones</span> <span>along the S</span><span>A</span><span>VZ, linked with stress field variations</span> <span>in spatial and temporal association with</span> –<span>short therm-</span> <span>seismicity</span> <span>and -long term- oblique </span><span>plate </span><span>convergence.</span> <span>Volcanoes and geothermal domains are mostly located along or at the intersection of margin-oblique fault zones (Andean Transverse Faults), and along margin-parallel strike slip zones, some which may cut the entire lithosphere (Liquiñe-Ofqui fault system). Wh</span><span>ereas</span><span> the big picture displays</span> <span>fluid flow straight to the surface, at close look significant offsets between crustal structures occur. 3D numerical models using conventional elasto-plastic rheology provide insights on the interaction of (i) an inflating magmatic cavity, (ii) a slipping fault zone, and (iii) regional tectonic stresses. Applying either (i) a magmatic overpressure or (ii) a given fault slip can trigger failure of the intervening rock, and generate either i) fault motion or ii) magmatic reservoir failure, respectively, but only for distances less than the structures' breadth even at low rock</span> <span>strength. However, at greater inter-distances the bedrock domain in between the fault zone and the magmatic cavity undergoes dilatational strain of the order of 1-5x10-5. This dilation opens the bedrock’s pore space and forms «pocket domains» that may store up-flowing over-pressurized fluids, which may then further chemically</span> interact<span> with the bedrock, for the length of time</span> <span>that</span> <span>these pockets remain open. These porous pockets</span> <span>can reach kilometric size, questioning their parental link with outcropping plutons along the margin. Moreover, bedrock permeability may also increase as fluid flow diminishes effective bedrock friction and cohesion. Comparison with rock experiments indicates that such stress and fluid pressure changes may eventually trigger failure at the intermediate timescale (repeated slip or repeated inflation). Finally, incorporating far field compression (iii)</span> <span>loads the bedrock to</span> <span>a state of stress at the verge of failure. Then, failure around the magmatic </span><span>reservoir</span><span> or </span><span>at</span> <span>the fault zone occurs for lower load</span><span>ing</span><span>.</span> <span>Permanent tectonic loading on the one hand, far field episodic seismic inversion of the stress field on the other, and localized failure all together promote a transient stress field, thus explaining the occurrence of transient fluid pathways on seemingly independent timescales. These synthetic models are then discussed with regards to specific cases along the SVZ, particularly the Tatara-San Pedro area (~36°S), where magnetotelluric profiles </span><span>document</span><span> conductive volumes at different depths underneath active faults, volcanic edifices and geothermal vents. We discuss the mechanical link between these deep sources and surface structures</span>.</p>


2021 ◽  
Vol 8 (1) ◽  
Author(s):  
Dyuti Prakash Sarkar ◽  
Jun-ichi Ando ◽  
Akihiro Kano ◽  
Hirokazu Kato ◽  
Gautam Ghosh ◽  
...  

AbstractDisentangling the temperature and depth of formation of fault rocks is critical for understanding their rheology, exhumation, and the evolution of fault zones. Estimation of fault rock temperatures mostly relies on conventional geothermometers of metamorphic minerals and element partitioning analysis, which are largely inapplicable in shallow crustal fault rocks. Here, we demonstrate the applicability of the carbonate clumped isotope thermometer in low-grade carbonate-bearing fault rocks from the Himalayan frontal wedge (northwest India). Coalescing carbonate clumped isotope thermometry and calcite e-twin morphology allows us to constrain the temperature and depth of formation of the two main thrusts of the Himalayan frontal wedge, the Nahan thrust (170 ± 10 °C; 6–7 km depth), and the Main Boundary thrust (262 ± 30 °C; 10–11 km depth). The integration of the adopted analytical techniques can promote the application of calcite-based clumped isotope thermometry to the fault zone processes and refinement of shallow crustal fault zone models.


2020 ◽  
Author(s):  
Ulrich Kelka ◽  
Thomas Poulet ◽  
Luk Peeters

<p>Fault and fracture networks can govern fluid flow patterns in the subsurface and predicting fluid flow on a regional scale is of interest in a variety of fields like groundwater management, mining engineering, energy, and mineral resources. Especially the pore fluid pressure can have a strong impact on the strength of fault zones and might be one of the drivers for fault reactivation. Reliable simulations of the transient changes in fluid pressure need to account for the generic architecture of fault zones that comprises strong permeability contrast between the fault core and damage zone.</p><p>Particularly, the distribution and connectivity of large-scale fault zones can have a strong impact on the flow field. Yet, modelling numerically such features in their full complexity remains challenging. Often faults zones are conceptualized as forming exclusively either barriers or conduits to fluid flow. However, a generic architecture of fault zones often comprises a discrete fault core surrounded by a diffuse damage zone and conceptualizing large scale discontinuities simply as a barrier or conduit is unlikely to capture the regional scale fluid flow dynamics. It is known that if the fault zone is hosted in low-permeability strata, such as clays or crystalline rocks, a transversal flow barrier can form along the fault core whereas the fracture-rich fault damage zone represents a longitudinal conduit. In more permeable host-rocks (i.e. sandstones or carbonates) the reverse situation can occur, and the permeability distributions in the damage zones can be governed by the abundance of low-permeability deformation features. A reliable numerical model needs to account for the difference and strong contrasts in fluid flow properties of the core and the damage zone, both transversally and longitudinally, in order to make prediction about the regional fluid flow pattern.</p><p>Here, we present a numerical method that accounts for the generic fault zone architecture as lower dimensional interfaces in conforming meshes during fluid flow simulations in fault networks. With this method we aim to decipher the impact of fault zone architecture on subsurface flow pattern and fluid pressure evolution in fractures and faulted porous media. The method is implemented in a finite element framework for Multiphysics simulations. We demonstrate the impact of considering the more generic geological structure of individual faults on the flow field by conceptualizing discontinuities either as barriers, conduits or as a conduit-barrier system and show were these conceptualizations are applicable in natural systems. We further show that a reliable regional scale fluid flow simulation in faulted porous media needs to account for the generic fault zone architecture. The approach is finally used to evaluate the fluid flow response of statistically parameterised faulted media, in order to investigate the impact and sensitivity of each variable parameter.</p>


2014 ◽  
Vol 18 (sup4) ◽  
pp. S4-1095-S4-1101 ◽  
Author(s):  
B. F. Luan ◽  
X. Y. Li ◽  
M. Zhang ◽  
Adrien Chapuis ◽  
Q. Liu

2021 ◽  
Author(s):  
Ake Fagereng ◽  
Adam Beall

<p>Current conceptual fault models define a seismogenic zone, where earthquakes nucleate, characterised by velocity-weakening fault rocks in a dominantly frictional regime. The base of the seismogenic zone is commonly inferred to coincide with a thermally controlled onset of velocity-strengthening slip or distributed viscous deformation. The top of the seismogenic zone may be determined by low-temperature diagenetic processes and the state of consolidation and alteration. Overall, the seismogenic zone is therefore described as bounded by transitions in frictional and rheological properties. These properties are relatively well-determined for monomineralic systems and simple, planar geometries; but, many exceptions, including deep earthquakes, slow slip, and shallow creep, imply processes involving compositional, structural, or environmental heterogeneities. We explore how such heterogeneities may alter the extent of the seismogenic zone.</p><p> </p><p>We consider mixed viscous-frictional deformation and suggest a simple rule of thumb to estimate the role of heterogeneities by a combination of the viscosity contrast within the fault, and the ratio between the bulk shear stress and the yield strength of the strongest fault zone component. In this model, slip behaviour can change dynamically in response to stress and strength variations with depth and time. We quantify the model numerically, and illustrate the idea with a few field-based examples: 1) earthquakes within the viscous regime, deeper than the thermally-controlled seismogenic zone, can be triggered by an increase in the ratio of shear stress to yield strength, either by increased fluid pressure or increased local stress; 2) there is commonly a depth range of transitional behaviour at the base of the seismogenic zone – the thickness of this zone increases markedly with increased viscosity contrast within the fault zone; and 3) fault zone weakening by phyllosilicate growth and foliation development increases viscosity ratio and decreases bulk shear stress, leading to efficient, stable, fault zone creep. These examples are not new interpretations or observations, but given the substantial complexity of heterogeneous fault zones, we suggest that a simplified, conceptual model based on basic strength and stress parameters is useful in describing and assessing the effect of heterogeneities on fault slip behaviour.         </p>


Geophysics ◽  
1967 ◽  
Vol 32 (2) ◽  
pp. 297-301 ◽  
Author(s):  
S. N. Domenico

A gravity profile was obtained from closely spaced readings along a traverse approximately nine miles in length across the San Andreas fault zone immediately south of Palmdale, California in the western Mojave Desert. Corrected gravity values show a slight but distinctive minimum associated with the fault zone which may be attributed to the reduced density of the shattered rock masses in the fault zone. The existence of this minimum suggests that major fault zones may be traced across terrain, on which surface expression of the fault does not exist, by successive profiles across the suspected position of the fault zone.


2007 ◽  
Vol 40 (4) ◽  
pp. 1586 ◽  
Author(s):  
N. Palyvos ◽  
D. Pantosti ◽  
L. Stamatopoulos ◽  
P. M. De Martini

In this communication we discuss reconnaissance geomorphological observations along the active Psathopyrgos and Rion-Patras (NE part) fault zones. These fault zones correspond to more or less complex rangefronts, the geomorphic characteristics of which provide hints on the details of the fault zone geometries, adding to the existing geological data in the bibliography. Aiming at the identification of locations suitable or potentially suitable for geomorphological and geological studies for the determination of fault slip rates in the Holocene, we describe cases of faulted Holocene landforms and associated surficial deposits. We also discuss problems involved in finding locations suitable for geological (paleoseismological) studies for the determination of the timing of recent earthquake ruptures, problems due to both man-made and natural causes.


2019 ◽  
Vol 10 (2) ◽  
pp. 561-567
Author(s):  
N. F. Vasilenko ◽  
A. S. Prytkov

In the northern Sakhalin Island, the tectonic activity of the fault zones is a potential threat to the industrial infrastructure of the petroleum fields. Recently, the background seismicity has increased at the Hokkaido‐Sakhalin fault that consists of several segments, including the Garomai active fault. In the studies of the regional deformation processes, it is important not only to analyze the seismic activity, but also to quantitatively assess the dynamics of deformation accumulation in the fault zones. In order to study the contemporary geodynamics of the Garomai fault, a local GPS/GLONASS network has been established in the area wherein trunk oil and gas pipelines are installed across the fault zone. Based on the annual periodic measurements taken in 2006–2016, we study the features of surface deformation and calculate the rates of displacements caused by the tectonic activity in the fault zone. During the survey period, no significant displacement of the fault wings was revealed. In the immediate vicinity of the fault zone, multidirectional horizontal displacements occur at a rate up to 1.6 mm/yr, and uplifting of the ground surface takes place at a rate of 3.4 mm/yr. This pattern of displacements is a reflection of local deformation processes in the fault zone. At the western wing of the fault, a maximum deformation rate amounts to 1110–6 per year. The fault is a boundary mark of a transition from lower deformation rates at the eastern wing to higher ones at the west wing. In contrast to the general regional compression setting that is typical of the northern Sakhalin Island, extension is currently dominant in the Garomai fault zone. The estimated rates of relative deformation in the vicinity of the Garomai fault give grounds to classify it as ‘hazardous’.


2003 ◽  
Vol 78-79 ◽  
pp. 191-195
Author(s):  
Kagan Tuncay ◽  
Gonca Ozkan ◽  
Peter Ortoleva
Keyword(s):  

Sign in / Sign up

Export Citation Format

Share Document