Abstract Many low-angle normal faults (dip ≤30°) accommodate tens of kilometers of crustal extension, but their mechanics remain contentious. Most models for low-angle normal fault slip assume vertical maximum principal stress σ1, leading many authors to conclude that low-angle normal faults are poorly oriented in the stress field (≥60° from σ1) and weak (low friction). In contrast, models for low-angle normal fault formation in isotropic rocks typically assume Coulomb failure and require inclined σ1 (no misorientation). Here, a data-based, mechanical-tectonic model is presented for formation of the Whipple detachment fault, southeastern California. The model honors local and regional geologic and tectonic history and laboratory friction measurements. The Whipple detachment fault formed progressively in the brittle-plastic transition by linking of “minidetachments,” which are small-scale analogs (meters to kilometers in length) in the upper footwall. Minidetachments followed mylonitic anisotropy along planes of maximum shear stress (45° from the maximum principal stress), not Coulomb fractures. They evolved from mylonitic flow to cataclasis and frictional slip at 300–400 °C and ∼9.5 km depth, while fluid pressure fell from lithostatic to hydrostatic levels. Minidetachment friction was presumably high (0.6–0.85), based upon formation of quartzofeldspathic cataclasite and pseudotachylyte. Similar mechanics are inferred for both the minidetachments and the Whipple detachment fault, driven by high differential stress (∼150–160 MPa). A Mohr construction is presented with the fault dip as the main free parameter. Using “Byerlee friction” (0.6–0.85) on the minidetachments and the Whipple detachment fault, and internal friction (1.0–1.7) on newly formed Reidel shears, the initial fault dips are calculated at 16°–26°, with σ1 plunging ∼61°–71° northeast. Linked minidetachments probably were not well aligned, and slip on the evolving Whipple detachment fault probably contributed to fault smoothing, by off-fault fracturing and cataclasis, and to formation of the fault core and fractured damage zone. Stress rotation may have occurred only within the mylonitic shear zone, but asymmetric tectonic forces applied to the brittle crust probably caused gradual rotation of σ1 above it as a result of: (1) the upward force applied to the base of marginal North America by buoyant asthenosphere upwelling into an opening slab-free window and/or (2) basal, top-to-the-NE shear traction due to midcrustal mylonitic flow during tectonic exhumation of the Orocopia Schist. The mechanical-tectonic model probably applies directly to low-angle normal faults of the lower Colorado River extensional corridor, and aspects of the model (e.g., significance of anisotropy, stress rotation) likely apply to formation of other strong low-angle normal faults.
more »
« less
Mechanical Implications of Creep and Partial Coupling on the World's Fastest Slipping Low‐Angle Normal Fault in Southeastern Papua New Guinea
We use densely spaced campaign GPS observations and laboratory friction experiments on fault rocks from one of the world's most rapidly slipping low‐angle normal faults, the Mai'iu fault in Papua New Guinea, to investigate the nature of interseismic deformation on active low‐angle normal faults. GPS velocities reveal 8.3 ± 1.2 mm/year of horizontal extension across the Mai'iu fault, and are fit well by dislocation models with shallow fault locking (above 2 km depth), or by deeper locking (from ~5–16 km depth) together with shallower creep. Laboratory friction experiments show that gouges from the shallowest portion of the fault zone are predominantly weak and velocity‐strengthening, while fault rocks deformed at greater depths are stronger and velocity‐weakening. Evaluating the geodetic and friction results together with geophysical and microstructural evidence for mixed‐mode seismic and aseismic slip at depth, we find that the Mai'iu fault is most likely strongly locked at depths of ~5–16 km and creeping updip and downdip of this region. Our results suggest that the Mai'iu fault and other active low‐angle normal faults can slip in large (Mw > 7) earthquakes despite near‐surface interseismic creep on frictionally stable clay‐rich gouges.
more » « less- NSF-PAR ID:
- 10381062
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 125
- Issue:
- 10
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Trench‐parallel translation of the Central American Forearc (CAFA) is the result of strain partitioning along the Cocos and Caribbean (CA) convergent margin. Unlike the tectonics of northwestern Costa Rica and El Salvador, CAFA‐CA relative motion in Nicaragua is not accommodated on margin‐parallel fault systems. Rather, the northwest‐trending dextral shear is accommodated on margin‐normal sinistral strike‐slip faults that approximate the motion of a margin‐parallel fault (i.e., bookshelf faulting). We compare a new Global Positioning System interseismic horizontal velocity field to analytical and numerical models to show that the bookshelf faulting model can produce the observed velocity field and provide insight into the kinematics and configuration of the margin‐normal fault system. We find that a fault system with 20 km‐long parallel to sub‐parallel margin‐normal sinistral faults, spaced ∼5 km apart, locked from the surface to 5 km depth, and with interseismic slip deficits of 4 mm yr−1, can replicate the observed velocity field. These findings have implications for the region's seismic hazard where shallow moderate‐magnitude earthquakes will have reoccurrence intervals of ∼50 years. These findings are also important for volcanic hazard estimation and unrest forecasting because the margin‐normal faults are in the volcanic arc and magma‐tectonic interactions have been documented along the CAFA.
-
more » « less
Abstract The frictional properties of faults control the initiation and propagation of earthquakes and the associated hazards. Although the ambient temperature and instantaneous slip velocity controls on friction in isobaric conditions are increasingly well understood, the role of normal stress on steady‐state and transient frictional behaviors remains elusive. The friction coefficient of rocks exhibits a strong dependence on normal stress at typical crustal depths. Furthermore, rapid changes in normal stress cause a direct effect on friction followed by an evolutionary response. Here, we derive a constitutive friction law that consistently explains the yield strength of rocks from atmospheric pressure to gigapascals while capturing the transient behavior following perturbations in normal stress. The model explains the frictional strength of a variety of sedimentary, metamorphic, and igneous rocks and the slip‐dependent response upon normal stress steps of Westerly granite bare contact and synthetic gouges made of quartz and a mixture of quartz and smectite. The nonlinear normal stress dependence of the frictional resistance may originate from the distribution of asperities that control the real area of contact. The direct and transient effects may be important for induced seismicity by hydraulic fracturing or for naturally occurring normal stress perturbations within fault zones in the brittle crust.
-
more » « less
Abstract Establishing a constitutive law for fault friction is a crucial objective of earthquake science. However, the complex frictional behavior of natural and synthetic gouges in laboratory experiments eludes explanations. Here, we present a constitutive framework that elucidates the rate, state, and temperature dependence of fault friction under the relevant sliding velocities and temperatures of the brittle lithosphere during seismic cycles. The competition between healing mechanisms, such as viscoelastic collapse, pressure‐solution creep, and crack sealing, explains the low‐temperature stability transition from steady‐state velocity‐strengthening to velocity‐weakening as a function of slip‐rate and temperature. In addition, capturing the transition from cataclastic flow to semi‐brittle creep accounts for the stabilization of fault slip at elevated temperatures. We calibrate the model using extensive laboratory data on synthetic albite and granite gouge, and on natural samples from the Alpine Fault and the Mugi Mélange in the Shimanto accretionary complex in Japan. The constitutive model consistently explains the evolving frictional response of fault gouge from room temperature to 600°C for sliding velocities ranging from nanometers to millimeters per second. The frictional response of faults can be uniquely determined by the in situ lithology and the prevailing hydrothermal conditions.
-
more » « less
Abstract Localized frictional sliding on faults in the continental crust transitions at depth to distributed deformation in viscous shear zones. This brittle‐ductile transition (BDT), and/or the transition from velocity‐weakening (VW) to velocity‐strengthening (VS) friction, are controlled by the lithospheric thermal structure and composition. Here, we investigate these transitions, and their effect on the depth extent of earthquakes, using 2D antiplane shear simulations of a strike‐slip fault with rate‐and‐state friction. The off‐fault material is viscoelastic, with temperature‐dependent dislocation creep. We solve the heat equation for temperature, accounting for frictional and viscous shear heating that creates a thermal anomaly relative to the ambient geotherm which reduces viscosity and facilitates viscous flow. We explore several geotherms and effective normal stress distributions (by changing pore pressure), quantifying the thermal anomaly, seismic and aseismic slip, and the transition from frictional sliding to viscous flow. The thermal anomaly can reach several hundred degrees below the seismogenic zone in models with hydrostatic pressure but is smaller for higher pressure (and these high‐pressure models are most consistent with San Andreas Fault heat flow constraints). Shear heating raises the BDT, sometimes to where it limits rupture depth rather than the frictional VW‐to‐VS transition. Our thermomechanical modeling framework can be used to evaluate lithospheric rheology and thermal models through predictions of earthquake ruptures, postseismic and interseismic crustal deformation, heat flow, and the geological structures that reflect the complex deformation beneath faults.
