Propagating Fault Tip Research Articles

Rheology Dependent on the Distance to the Propagating Thrust Tip—(Ultra‐)Mylonites and Pseudotachylytes of the Silvretta Basal Thrust

AbstractTo evaluate how the presence of pseudotachylytes affects the strength of crustal rocks, deformed pseudotachylytes and their relationship with pristine pseudotachylytes at the base of the Silvretta nappe are analyzed. Pseudotachylytes formed associated with high‐stress crystal plasticity (σd > 400 MPa), as indicated by twinned amphiboles in gneisses. Mylonitic quartz clasts enclosed within deformed pseudotachylytes and mylonitic vein‐quartz, hosting folded pseudotachylyte injection veins, reflect creep at lower stresses (ca. 100 MPa) after seismic rupturing. Deformed pseudotachylytes can be crosscut by pristine pseudotachylytes, indicating a second, independent stage of coseismic rupturing after creep. The evidence of dynamic dislocation creep of quartz and the presence of stilpnomelane and epidote associated with all fault rocks indicate similar ambient greenschist facies conditions during all deformation stages. Whereas the intermediate stage of creep is interpreted to represent deformation at large distance to the propagating thrust tip, the pristine pseudotachylytes represent the last stage of rupturing eventually leading to nappe decoupling from its basement. Gneiss clasts in an ultramylonitic matrix (i.e., deformed pseudotachylyte) reveal that pseudotachylytes have a lower strength during creep in relation to the hosting gneisses. In contrast, during coseismic high‐stress crystal plasticity, the coarse gneisses accumulate a higher amount of strain. This strength‐relationship explains that only those rocks rupture, which have not been previously deformed before. The study demonstrates the importance of different strengths of crustal rocks at specific stress‐ and strain‐rate conditions in dependence on the distance to the propagating fault tip.

Tidally induced fault motion within Europa's ice shell and implications for subsurface communication development

The processes that maintain a subsurface ocean at Europa may also cause geophysical hazards within the overlying ice shell that are important to understand for any subsurface access mission. In this study, we applied 3D finite element models to better constrain the potential magnitude of tidally induced fault motion in the ice shell that a cryobot and communication hardware could experience during a subsurface mission. Models showed that displacement along a fault resulting from a static tidal stress at orbital perijove was <5 cm at the subjovian and ∼ 5–82 cm at Thera Macula, depending on fault orientation, changes in tidal stress with depth, and coefficient of friction. Time-dependent tidal forcing over one tidal cycle resulted in fault displacements of ∼0.50–2.9 m for end-member fault geometries at Thera Macula. The potential strain imposed on a communication tether from the fault movement will depend greatly on the currently unknown adhesion of the tether to the ice as well as potential fault tip propagation, yet these analyses provide first-order approximations for end-member strain conditions using measurements of the motion across a fault for a cryobot at 1 km depth. Results suggest it is important to consider motion throughout a tidal cycle and that certain locations and surface features on Europa may prove to be more challenging for a subsurface mission than others. With NASA's upcoming Europa Clipper and ESA's JUICE missions, fault displacements on Europa could be observed, providing further information on regions or features more/less hazardous to cryobot missions.

Stresses during pseudotachylyte formation - Evidence from deformed amphibole and quartz in fault rocks from the Silvretta basal thrust (Austria)

Pseudotachylyte-bearing, amphibole-rich gneisses with quartz-rich layers concordant to the foliation from the base of the Silvretta nappe were analyzed by polarized light microscopy and scanning electron microscopy to constrain the deformation history. Amphibole shows microfractures, kink bands and mechanical twins associated with pseudotachylytes. (−101)[101] twins document high differential stresses of ≥ 400 MPa. Quartz grains in contact with twinned amphibole in gneisses show subbasal deformation lamellae, short-wavelength undulatory extinction and new grains along cracks with diameters of < 15 μm and crystallographic preferred orientation (CPO) controlled by the orientation of host grains. This quartz microstructure indicates dislocation glide-controlled deformation coeval with high-stress amphibole deformation and pseudotachylyte formation.Quartz-rich layers concordant to the gneiss' foliation commonly show evidence of dislocation creep by homogenous recrystallized grain aggregates with CPO indicating <a> dislocation glide. Systematic differences in distribution of the recrystallized grains and grain sizes, indicate deformation at differential stresses from a few tens to hundred MPa. Yet, no evidence of increasing strain towards the contact to amphibole-rich gneisses is observed. In contrast, a decreasing number of recrystallized grains towards the contact occurs locally. In relation to pseudotachylytes, the recrystallized quartz is affected by shear fractures offsetting the recrystallized microstructure. This, and the missing or decreasing strain gradient of dislocation creep towards the contact to amphibole-rich gneisses, indicate that 1) dislocation creep in quartz-rich layers is preceding and 2) not directly related to the high-stress and high strain-rate deformation, ruling out a major influence of thermal runaway for the formation of the pseudotachylytes. The differences in the deformation behavior is interpreted to be controlled by the distance to the propagating fault tip. This study demonstrates the potential of microstructures to constrain the deformation history controlling seismic faulting and creep at hypocentral depth not directly accessible for in-situ stress measurements.

Origin and growth mechanisms of strike-slip faults in the central Tarim cratonic basin, NW China

Through fault structure analysis and chronology study, we discuss the origin and growth mechanisms of strike-slip faults in the Tarim Basin. (1) Multiple stages strike-slip faults with inherited growth were developed in the central Tarim cratonic basin. The faults initiation time is constrained at the end of Middle Ordovician of about 460 Ma according to U-Pb dating of the fault cements and seismic interpretation. (2) The formation of the strike-slip faults was controlled by the near N-S direction stress field caused by far-field compression of the closing of the Proto-Tethys Ocean. (3) The faults localization and characteristics were influenced by the pre-existing structures of the NE trending weakening zones in the basement and lithofacies change from south to north. (4) Following the fault initiation under the Andersonian mechanism, the strike-slip fault growth was dominantly fault linkage, associated with fault tip propagation and interaction of non-Andersonian mechanisms. (5) Sequential slip accommodated deformation in the conjugate strike-slip fault interaction zones, strong localization of the main displacement and deformation occurred in the overlap zones in the northern Tarim, while the fault tips, particularly of narrow-deep grabens, and strike-slip segments in thrust zones accumulated more deformation and strain in the Central uplift. In conclusion, non-Andersonian mechanisms, dominantly fault linkage and interaction, resulted in the small displacement but long intraplate strike-slip fault development in the central Tarim Basin. The regional and localized field stress, and pre-existing structures and lithofacies difference had strong impacts on the diversity of the strike-slip faults in the Tarim cratonic basin.

Extensional fault-related folding of the North West shelf, Western Australia

The occurrence of fault-related folds in extensional settings has been gaining increased recognition during the past 20 yr. Passive margins and rift basins often present intriguing structural complexities that are potential sites for the development of extensional fault–related folds. Fault-related folds are important anticlinal traps for hydrocarbon exploration, and understanding their evolution is critical for derisking the timing of trap formation and hydrocarbon generation. This study presents folds from three areas across the North West shelf (NWS) of Western Australia and aims to determine their relationship to associated extensional faults. A series of kinematic forward models are presented as well as detailed structural and stratigraphic interpretation of folding at the three locations. The models suggest each of the three examples have formed under extension that impacted the region during the Mesozoic. The modeling also highlights different characteristics of folding that are caused by fault geometry (ramp–flat) and a propagating fault tip (trishear) that can be used as diagnostic features to distinguish different fold mechanisms. Historical exploration across the NWS has targeted a relatively small number of these anticlinal traps, some of which have been successful (e.g., Marina-1 well, Petrel subbasin) and some have not (e.g., Candace-1, Sholl Island fault anticline). This paper suggests that the prospectivity of anticlinal structures should be more carefully assessed to determine whether or not they are related to extensional activity and hence likely to have formed prior to the bulk of the hydrocarbon charge. Similar structures undoubtedly exist in other basins that have experienced repeated extension.

Variable fault tip propagation rates affected by near-surface lithology and implications for fault displacement hazard assessment

The fabric of reverse fault zones close to the surface is usually partitioned in between a narrow discrete rupture zone and a more distributed one, where folding is predominant. This makes quite challenging the adoption of proper setbacks in surface rupture hazard studies for critical facilities or microzoning. Some of the parameters controlling fault zone fabric are related to mechanics of near-surface geology (lithology, overburden thickness, cohesion and water content) whose interaction is complex and only partially understood. Nevertheless, these can be seldom measured or derived. Kinematic models, conversely, express such an interaction of complex variables as simple synthetic parameters, such as the amount of upward propagation of the fault tip for unit of slip, usually referred to as the P/S ratio (Propagation on Slip). Here, we discuss results on the trishear kinematic inverse modeling of a contractional fault propagation fold at Monte Netto Hill (Capriano del Colle, N. Italy), observing a two-stage fault and fold growth evolution, marked by a significant shift in the P/S parameter. At this site, exceptional sequence of exposures due to ca. 10 years of quarry excavations allowed to obtain a series of cross-sections across the fault zone. We use this detailed, high-resolution, example as a natural “analogue” for more general, large-scale surface ruptures involving a thick alluvial cover, a very common setting for the siting of critical facilities.During the early stage of displacement, the fault cut through clast-supported fluvial gravels with a high propagation rate (P/S = 7) and a discrete rupture width. Then, during the latest movements of the thrust, fault tip propagation slowed down to P/S ≈ 2.9, as the fault started cutting through several stacked bodies of clays and silty clays, pedogenized aeolian silts and overbank deposits, causing a pronounced folding of the layers over a wider deformation zone. These results strongly suggest that lithological changes in the underlying shallow stratigraphy, common in an alluvial plain depositional setting, would significantly affect the potential for surface faulting across the same tectonic structure, with relevant implications in the fault displacement hazard assessment.

Structural evolution of the active Waimea-Flaxmore Fault System in the Nelson-Richmond urban area, South Island, New Zealand

ABSTRACTThe Waimea-Flaxmore Fault System (W-FFS) comprises NE-striking, SE-dipping reverse faults superposing Eastern Province terranes on Cenozoic sedimentary units of the Moutere Depression. We present a revised geological map and cross sections of the W-FFS in the Nelson-Richmond urban area, aimed at reconstructing the structural evolution during multiple phases of reactivation. Our data and structural interpretations depict the inversion of the W-FFS from an inherited system of Late Cretaceous-Eocene normal faults to multiple reverse fault splays that propagated through the cover sequence over the last 10 Ma. Fault tip propagation with folding of the cover sequence persisted during the Quaternary, with breaching of the Waimea and Flaxmore faults to the surface and propagation of a splay thrust (Tahunanui Fault) still concealed beneath Quaternary deposits. Ongoing crustal shortening associated with the Pacific-Australian plate boundary in the northern South Island creates favourable conditions for reactivation of the long-lived W-FFS in the contemporary stress field.

Seismic and Aseismic Fault Growth Lead to Different Fault Orientations

AbstractOrientations of natural fault systems are subject to large variations. They often contradict classical Coulomb failure theory as they are misoriented relative to the regional Andersonian stress field. This is ascribed to local effects of structural or stress heterogeneities and reorientations of structures or stresses on the long term. To better understand the relation between fault orientation and regional stresses, we simulate spontaneous fault growth and its effect on the stress field. Our approach incorporates earthquake rupture dynamics, viscoelastoplastic brittle deformation and a rate‐ and state‐dependent friction formulation in a continuum mechanics framework. We investigate how strike‐slip faults orient according to local and far‐field stresses during their growth. We identify two modes of fault growth, seismic and aseismic, distinguished by different fault angles and slip velocities. Seismic fault growth causes a significant elevation of dynamic stresses and friction values ahead of the propagating fault tip. These elevated quantities result in a greater strike angle relative to the maximum principal regional stress than that of a fault segment formed aseismically. When compared to the near‐tip time‐dependent stress field the fault orientations produced by both growth modes follow the classical failure theory. We demonstrate how the two types of fault growth may be distinguished in natural faults by comparing their angles relative to the original regional maximum principal stress. A stress field analysis of the Landers‐Kickapoo fault suggests that an angle greater than ∼25° between two faults indicates seismic fault growth.

The Zhayao tectonic window of the Jurassic Yuantai thrust system in Liaodong Peninsula, NE China: Geometry, kinematics and tectonic implications

Liaodong Peninsula on the northeastern North China Block that witnessed the spectacular Late Mesozoic tectonics is a key area to probe the transition from compression to extension. Borehole data clearly reveals that the Cambrian limestone overlay the Jurassic limestone, which identified the Yuantai thrust system in Liaodong Peninsula. Zircon grains from a granodiorite dike that intruded the thrust system yielded U-Pb age of 128.3 ± 1.3 Ma, further constraining the timing of the thrusting between Late Jurassic to early Early Cretaceous. A detailed field study of three-dimensional fault-zone exposures in the Zhayao tectonic window are conducted to decipher the thrust tectonics before the extension in Liaodong Peninsula. The key observations are: (1) meter-scale asymmetric folds in deformed mid-Cambrian rocks; (2) decameter-scale recumbent folds and secondary imbricate thrusts in deformed Middle Jurassic rocks; and (3) non-foliated breccia, fault gouge, and brittle porphyroclasts adjacent to fault planes. The combined brittle–ductile style of deformation indicates that the fault zone operated under near-surface conditions. In addition, we infer that the folds were formed first at the propagating fault tip and subsequently dismembered by normal faults. We interpret the Zhayao thrust, which is located around the margin of the tectonic window, as the frontal thrust of the Yuantai thrust system. Furthermore, by applying the theory of fault-related folds, we deduce that the Zhayao frontal thrust is a secondary duplex. Based on our field observations and borehole data, we propose an imbricate thrust model for the fold-and-thrust belt of the southern Liaodong Peninsula. The structural data indicate that the transport direction of the Yuantai thrust system was top-to-the-southeast, suggesting that the thrust system was resulted due to northwest-ward subduction of the paleo-Pacific plate during the Mesozoic.

Impact of fault evolution in Gilbert-type fan deltas in the Evrostini area, south-central Gulf of Corinth, Greece

Abstract In order to understand how syn-rift deposits are associated with border fault evolution, the outcrops around the Evrostini area in the southern Gulf of Corinth, Greece, are investigated. Two major border fault systems are identified in the study area, which are the Vela–Valimi–Xylokastro and Akrata–Xylokastro border fault complexes. The Vela–Valimi and Xylokastro faults were active in the Calabrian age when the WSW-trending linking faults were initiated at the Evrostini Mountain area and propagated eastwards. The Xylokastro Fault was reactivated in the Middle Pleistocene, when the Akrata Fault was created. It is proposed that WSW-trending faults were developed to link the Akrata and Xylokastro border faults. Longitudinal propagation of the rift margin in the area was segmented by three N–S striking transfer faults. These transfer faults are defined by N–S-aligned discontinuity of abrupt fault displacement termination, large horizontal steps of E–W-striking fault segments, and lateral offsets of syn-rift deposition with major changes in facies. The regional drainages were strongly influenced by border fault growth and along-strike segmentation. Near the tips of the border faults, the drainage may be forced to shift along with fault tip propagation. In contrast, in the boundaries of rift segmentation, the drainage maintains its course and proceeded parallel to the orientation of N–S striking transfer faults. Drainage patterns and their configurations with border fault systems from the field study provide an analog model to predict the sub-surface syn-rift deposits in buried rift basins, such as the North Sea.

Structural evolution of a normal fault array in the Gambier Embayment, offshore Otway Basin, South Australia: insights from 3D seismic data

ABSTRACTThis study was undertaken to determine the structural evolution of a normal fault array using detailed kinematic analysis of normal fault tip propagation and linkage, adding to the growing pool of research on normal fault growth. In addition, we aim to provide further insight into the evolution of the offshore Otway Basin, Australia. We use three-dimensional (3D) seismic reflection data to analyse the temporal and spatial evolution of a Late Cretaceous–Cenozoic age normal fault array located in the Gambier Embayment of the offshore Otway Basin, South Australia. The seismic reflection data cover a NW–SE-oriented normal fault array consisting of six faults, which have grown from the linkage of numerous, smaller segments. This fault array overlies and has partial dip-linkage to E–W-striking, basement-involved faults that formed during the initial Tithonian–Barremian rifting event in the Otway Basin. Fault displacement analysis suggests four key stages in the post-Cenomanian growth history of the upper array: (1) nucleation of the majority of faults resulting from resumed crustal extension during the early Late Cretaceous; (2) an intra-Late Cretaceous period of general fault dormancy, with the nucleation of only one newly formed fault; (3) latest Cretaceous nucleation of another newly formed fault and further growth of all other faults; and (4) continued growth of all faults, leading to the formation of the Cenozoic Gambier Sub-basin in the Otway Basin. Our analysis also demonstrates that Late Cretaceous faults, which are located above and dip-link to basement-involved faults, display earlier nucleation and greater overall throw and length, compared with those which do not link to basement-involved faults. This is likely attributed to increased rift-related stress concentrations in cover sediments above the upper tips of basement faults. This study improves our understanding of the geological evolution of the presently under-explored Gambier Embayment, offshore Otway Basin, South Australia by documenting the segmented growth style of a Late Cretaceous normal fault array that is located over, and interacts with, a reactivated basement framework.

Subsurface structural interpretation by applying trishear algorithm: An example from the Lenghu5 fold-and-thrust belt, Qaidam Basin, Northern Tibetan Plateau

The application of trishear algorithm, in which deformation occurs in a triangle zone in front of a propagating fault tip, is often used to understand fault related folding. In comparison to kink-band methods, a key characteristic of trishear algorithm is that non-uniform deformation within the triangle zone allows the layer thickness and horizon length to change during deformation, which is commonly observed in natural structures. An example from the Lenghu5 fold-and-thrust belt (Qaidam Basin, Northern Tibetan Plateau) is interpreted to help understand how to employ trishear forward modelling to improve the accuracy of seismic interpretation. High resolution fieldwork data, including high-angle dips, ‘dragging structures’, thinning hanging-wall and thickening footwall, are used to determined best-fit trishear model to explain the deformation happened to the Lenghu5 fold-and-thrust belt. We also consider the factors that increase the complexity of trishear models, including: (a) fault-dip changes and (b) pre-existing faults. We integrate fault dip change and pre-existing faults to predict subsurface structures that are apparently under seismic resolution. The analogue analysis by trishear models indicates that the Lenghu5 fold-and-thrust belt is controlled by an upward-steepening reverse fault above a pre-existing opposite-thrusting fault in deeper subsurface. The validity of the trishear model is confirmed by the high accordance between the model and the high-resolution fieldwork. The validated trishear forward model provides geometric constraints to the faults and horizons in the seismic section, e.g., fault cutoffs and fault tip position, faults’ intersecting relationship and horizon/fault cross-cutting relationship. The subsurface prediction using trishear algorithm can significantly increase the accuracy of seismic interpretation, particularly in seismic sections with low signal/noise ratio.

Main Frontal thrust deformation and topographic growth of the Mohand Range, northwest Himalaya

The Main Frontal thrust (MFT) uplifts the Himalayan topographic front. Deciphering MFT deformation kinematics is crucial for understanding how the orogen accommodates continuing continental collision and assessing associated hazards. Here, we (a) detail newly discovered fault-zone exposures along the MFT at the Mohand Range front in northwestern India and (b) apply contemporary fault zone theory to show that the MFT is an emergent fault with a well-developed fault zone overlain by uplifted Quaternary gravels over a horizontal length of ∼700 m. Northward from the front, the fault zone grades from a central, gouge-dominated core to a hanging-wall, rock-dominated damage zone. We observed incohesive, non-foliated breccia, fault gouge, and brittle deformation microstructures within the fractured country rocks (Middle Siwaliks) and outcrop scale, non-plunging folds in the proximal hanging wall. We interpret these observations to suggest that (1) elastico-frictional (brittle) deformation processes operated in the fault zone at near surface (∼1–5 km depth) conditions and (2) the folds formed first at the propagating MFT fault tip, then were subsequently dismembered by the fault itself. Thus, we interpret the Mohand Range as a fault-propagation fold driven by an emergent MFT in contrast to the consensus view that it is a fault-bend fold. A fault-propagation fold model is more consistent with these new observations, the modern range-scale topography, and existing erosion estimates. To further evaluate our proposed structural model, we used a Boundary Element Method-based dislocation model to simulate topographic growth from excess slip at a propagating fault tip. Results show that the frontal topography could have evolved by slip along a (a) near-surface fault plane consistent with the present-day MFT location, or (b) blind MFT at ∼3 km depth farther north near the drainage divide. Comparing modelled vs. measured high resolution (∼16 cm) topographic profiles for each case provides permissible end-member scenarios of an either dynamically-evolving, high erosion, northward-migrating frontal scarp or a static, low, and symmetric, MHT-related fold, respectively. Our integrated approach is expected to deliver an improved understanding of coupled fault-generated deformation and topographic growth that may be applied more broadly across the entire Himalayan front.

The unusual 3D interplay of basement fault reactivation and fault-propagation-fold development: A case study of the Laramide-age Stillwell anticline, west Texas (USA)

Subsurface fault geometries have a systematic influence on folds formed above those faults. We use the extraordinarily well-exposed fold geometries of the Laramide-age Stillwell anticline in west Texas (USA) to develop a strain-predictive model of fault-propagation fold formation. The anticline is a 10-km long, NW-trending, NE-vergent, asymmetric fold system with an axis that displays a map-view left-stepping, en echelon pattern. We integrated field observations, geologic and structural data, cross-sections, and 2D kinematic modeling to establish an unusual 3D two-stage model of contractional fold formation, including: 1) reverse reactivation of a pre-existing, NW-striking, SW-dipping, left-stepping, en echelon normal fault system in Paleozoic basement rocks to generate monoclinal flexures in overlying layered Cretaceous carbonate rocks; and 2) the formation of a subsequent flat-ramp fault system that propagated horizontally along a mechanically-weak, clay-rich Cretaceous unit before ramping up at the hinge of the pre-existing monocline system. Strain is focused within the forelimb of the system, in front of the propagating fault tip, and is accommodated by a combination of interlayer slip, flat-ramp faulting, and fracturing proximal to planes of slip. This strain predictive model can be applied to similar, less-well-exposed contractional systems worldwide and provides a new, unusual example of Laramide-age contractional deformation.

Off-fault tip splay networks: A genetic and generic property of faults indicative of their long-term propagation

We use fault maps and fault propagation evidences available in the literature to examine geometrical relations between parent faults and off-fault splays. The population includes 47 worldwide crustal faults with lengths from millimetres to thousands of kilometres and of different slip modes. We show that fault splays form adjacent to any propagating fault tip, whereas they are absent at non-propagating fault ends. Independent of fault length, slip mode, context, etc., tip splay networks have a similar fan shape widening in direction of long-term propagation, a similar relative length and width (∼30 and∼10% of parent fault length, respectively), and a similar range of mean angles to parent fault (10–20°). We infer that tip splay networks are a genetic and a generic property of faults indicative of their long-term propagation. Their generic geometrical properties suggest they result from generic off-fault stress distribution at propagating fault ends.

Early synrift reservoir development on the flanks of extensional forced folds: A seismic-scale outcrop analog from the Hadahid fault system, Suez rift, Egypt

Forced folds typically develop above the tips of propagating normal faults in rifts that contain thick, prerift salt or mudstone sequences. This structural style is associated with the deposition of wedge-shaped synrift deposits that thin and onlap toward monoclinal growth folds overlying the vertically restricted fault tips. Subtle stratigraphic traps may develop on the flanks of these folds although, because of limited seismic resolution and sparse well data, the architecture, thickness, and distribution of these early synrift reservoirs are difficult to predict. To improve our understanding of early synrift reservoir development on the flanks of forced folds, we focus on seismic-scale outcrop analogs along the Hadahid fault system, Suez rift, Egypt. Our data indicate that forced folding dominated during early rifting and that the onset of folding was diachronous along strike. Fluvial systems incised the rotating monocline limbs, leading to the formation of valley-like erosional relief along the base synrift unconformity. Reservoir-prone fluvial facies are only locally developed along the forced-fold flank, with their distribution related to the degree of sediment bypass downdip into the adjacent basin. Early synrift relief not filled by fluvial strata was backfilled by transgressive, tidally influenced, reservoir-prone facies, with carbonates being locally developed in areas of low clastic sediment supply. Further extension and fault-tip propagation led to amplification of the forced folds, and deposition of shallow marine-to-shelf parasequences that became thinner toward the growing folds. Although displaying greater strike continuity than the underlying fluvial or tidal reservoirs, shoreface sandstone reservoirs amalgamate onto the flanks of the forced folds and may be absent toward the fold crest. This seismic-scale outcrop analog helps us better understand the subseismic stratigraphic architecture and facies distributions of early synrift reservoirs on the flanks of extensional forced folds. Observations from this and other well-exposed outcrop analogs should help reduce subsurface uncertainty and risk when exploring for hitherto under-explored, subtle, early synrift stratigraphic traps.

Defining a 3-dimensional trishear parameter space to understand the temporal evolution of fault propagation folds

The application of trishear, in which deformation occurs in a triangular zone in front of a propagating fault tip, is often used to understand fault related folding. A key element of trishear, in comparison to kink-band methods, is that non-uniform deformation within the triangle zone allows the layer thickness and length to change during deformation. By varying three controlling parameters independently (trishear propagation/slip ratio, trishear apical angle and fault dip), we construct a three-dimensional parameter space to demonstrate the variability of resultant geometry feasible with trishear. We plot published natural examples in this parameter space and identify two clusters and show that the most applicable typical trishear propagation/slip ratio is 2 to 3, while the trishear apical angle varies from 30° to 100°. We propose that these findings can help estimate the best-fit parameters for natural structures. We then consider the temporal evolution of specific geometric examples and factors that increase the complexity of trishear including: (1) fault-dip changes and (2) pre-existing faults.To illustrate the applicability of the parameter space and complex trishear models to natural examples, we apply our results to a sub-surface example from the Qaidam basin in northern Tibetan Plateau.

Evolution of coseismic fault-related folds induced by the Chi–Chi earthquake: A case study of the Wufeng site, Central Taiwan by using 2D distinct element modeling

Coseismic fault-propagation folds are generally associated with blind faults, which are of considerable academic interest and have recently been recognized as critical for assessing seismic hazards. The Wufeng excavation site was characterized by a major east-dipping basal thrust exhibiting a dip angle of 34° and 2 opposing vergent thrusts generated a pop-up anticlinal fold in soil cover induced by the 1999 Chi–Chi earthquake (Mw=7.6). In order to represent the soft soil behavior, we set up a direct shear simulation test to obtain valid parameters for soil. A series of 2D distinct element models possessing different bonding types and strengths was conducted to determine the deformation pattern near the surface and the evolution of the fault tip propagation. The coseismic deformation features of soft and plastic soil cover like the mutative limb thickness, complex ruptures and the overturned forelimbs in small-scale caused by the propagation of the blind fault tip were accurately predicted by contact-bond model in which grains were allowed to rotate and slip without cements breaking. Our results show that the alternative thrusts and pop-up structure developed before the main basal thrust fault ruptured through the ground surface at the Wufeng excavation site. We evaluated the slipping distance of the main fault at approximately 6m in Wufeng during the 1999 Chi–Chi earthquake.

Critical state finite element models of contractional fault-related folding: Part 1. Structural analysis

Fourteen Lagrangian finite element models with a critical state mechanics constitutive model illustrate some of the primary controls on the formation of fault propagation, fault bend, and diffuse folds. The models demonstrate how variable mechanical stratigraphy, initial fault dip and inter-layer detachments affect the way faults propagate and thus exert a significant control on resultant fold layer geometry. For example, models of uniform sandstone properties exhibit efficient strain localization and clear patterns of fault tip propagation. Uniform shale properties tend to inhibit fault propagation due to distributed plastic deformation. Models with mixed inter-layered sandstone and shale deform in a disharmonic manner, resembling lobate–cuspate arrangements that are common to many folds observed in outcrop. Detachments accommodate shortening by bed-parallel slip, resulting in fault-bend fold kinematics and poorly expressed fault propagation across layers. Structural analysis of the numerical model results reveals that contractional deformation is a composite of lateral compaction, pure shear shortening, fault propagation along narrowly localized zones of reverse shear, and flexure of layers. The relative proportions of these shortening components vary in time and with mechanical properties (shale vs. sandstone). Depth-to-detachment calculations performed on selected numerical models suggest that reasonably accurate predictions can be made for detachment fold-thrust belts and toe-of-slope contractional systems. However, our study suggests that applications to mid-crustal level detachments within the basement beneath a sedimentary cover may be inaccurate.

Scaling of fault damage zones with displacement and the implications for fault growth processes

[1] Knowledge of the spatial extent of damage surrounding fault zones is important for understanding crustal fluid flow and also for understanding the physical processes and mechanics by which fault zones develop with slip. There are few data available on the scaling of the fault damage zone with fault displacement, and of those that exist, deriving scaling relationships is hampered by comparing faults that run through different lithologies, have formed at different crustal depths or tectonic regimes (e.g., normal versus strike‐slip movement). We describe new data on the microfracture damage zone width from small displacement fault zones within the Atacama fault zone in northern Chile that formed at ∼6 km depth within a dioritic protolith. The microfracture damage zone is shown by an alteration halo surrounding the faults in which the density of the microfractures is much greater than background levels in the undeformed protolith. The data show that damage zone width increases with fault displacement and there appears to be a zero intercept to this relationship, meaning that at zero displacement, there is no microfracture damage zone. This is supported by field observations at fault tips that show a tapering out of fault damage zones. These data, combined with data from the literature, indicate that this same relationship might hold for much larger displacement faults. There is also a distinct asymmetry to the fracture damage. Several processes for the development of the observed scaling are discussed. The widely accepted theory of a process zone predicts that fault damage zone width increases with fault length and thus should always be largest at a propagating fault tip where displacement is lowest. This prediction is opposite to that seen in the current data set, leading to suggestion that other processes, such as damage zone growth with increasing displacement due to geometric irregularities or coseismic damage formation might better explain the spatial extent of damage surrounding even low‐displacement faults.