skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A pressure solution flow law for the seismogenic zone: Application to Cascadia
We develop a linear viscous constitutive relationship for pressure solution constrained by models of deformed metasedimentary rocks and observations of exposed rocks from ancient subduction zones. We include pressure and temperature dependence on the solubility of silica in fluid by parameterizing a practical van’t Hoff relationship. This general flow law is well suited for making predictions about interseismic behavior of subduction zones. We apply the flow law to Cascadia, where thermal structure, geometry, relative plate velocity, and Global Positioning System velocity field are well constrained. Results are consistent with the temperature conditions at which resolvable ductile strain is recorded in subducted mudstones (at depths near the updip limit of the seismogenic zone) and with relative plate motion accommodated completely by viscous deformation (at depths near the downdip limit of the seismogenic zone). The flow law also predicts the observed forearc tapering of slip rate deficit with depth.  more » « less
Award ID(s):
2214324 2054522
PAR ID:
10508663
Author(s) / Creator(s):
;
Editor(s):
Hodges, K
Publisher / Repository:
Science
Date Published:
Journal Name:
Science Advances
Volume:
10
Issue:
4
ISSN:
2375-2548
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract The plate interface undergoes two transitions between seismogenic depths and subarc depths. A brittle-ductile transition at 20–50 km depth is followed by a transition to full viscous coupling to the overlying mantle wedge at ∼80 km depth. We review evidence for both transitions, focusing on heat-flow and seismic-attenuation constraints on the deeper transition. The intervening ductile shear zone likely weakens considerably as temperature increases, such that its rheology exerts a stronger control on subduction-zone thermal structure than does frictional shear heating. We evaluate its role through analytic approximations and two-dimensional finite-element models for both idealized subduction geometries and those resembling real subduction zones. We show that a temperature-buffering process exists in the shear zone that results in temperatures being tightly controlled by the rheological strength of that shear zone’s material for a wide range of shear-heating behaviors of the shallower brittle region. Higher temperatures result in weaker shear zones and hence less heat generation, so temperatures stop increasing and shear zones stop weakening. The net result for many rheologies are temperatures limited to ≤350–420 °C along the plate interface below the cold forearc of most subduction zones until the hot coupled mantle is approached. Very young incoming plates are the exception. This rheological buffering desensitizes subduction-zone thermal structure to many parameters and may help explain the global constancy of the 80 km coupling limit. We recalculate water fluxes to the forearc wedge and deep mantle and find that shear heating has little effect on global water circulation. 
    more » « less
  2. null (Ed.)
    Drilling the input materials of the north Sumatran subduction zone, part of the 5000 km long Sunda subduction zone system and the origin of the Mw ~9.2 earthquake and tsunami that devastated coastal communities around the Indian Ocean in 2004, was designed to groundtruth the material properties causing unexpectedly shallow seismogenic slip and a distinctive forearc prism structure. The intriguing seismogenic behavior and forearc structure are not well explained by existing models or by relationships observed at margins where seismogenic slip typically occurs farther landward. The input materials of the north Sumatran subduction zone are a distinctively thick (as thick as 4–5 km) succession of primarily Bengal-Nicobar Fan–related sediments. The correspondence between the 2004 rupture location and the overlying prism plateau, as well as evidence for a strengthened input section, suggest the input materials are key to driving the distinctive slip behavior and long-term forearc structure. During Expedition 362, two sites on the Indian oceanic plate ~250 km southwest of the subduction zone, Sites U1480 and U1481, were drilled, cored, and logged to a maximum depth of 1500 meters below seafloor. The succession of sediment/rocks that will develop into the plate boundary detachment and will drive growth of the forearc were sampled, and their progressive mechanical, frictional, and hydrogeological property evolution will be analyzed through postcruise experimental and modeling studies. The large penetration depths with good core recovery and successful wireline logging in the challenging submarine fan materials will enable evaluation of the role of thick sedimentary subduction zone input sections in driving shallow slip and amplifying earthquake and tsunami magnitudes at the Sunda subduction zone and globally at other subduction zones where submarine fan–influenced sections are being subducted. 
    more » « less
  3. null (Ed.)
    Drilling the input materials of the north Sumatran subduction zone, part of the 5000 km long Sunda subduction zone system and the origin of the Mw ~9.2 earthquake and tsunami that devastated coastal communities around the Indian Ocean in 2004, was designed to groundtruth the material properties causing unexpectedly shallow seismogenic slip and a distinctive forearc prism structure. The intriguing seismogenic behavior and forearc structure are not well explained by existing models or by relationships observed at margins where seismogenic slip typically occurs farther landward. The input materials of the north Sumatran subduction zone are a distinctively thick (as thick as 4–5 km) succession of primarily Bengal-Nicobar Fan–related sediments. The correspondence between the 2004 rupture location and the overlying prism plateau, as well as evidence for a strengthened input section, suggest the input materials are key to driving the distinctive slip behavior and long-term forearc structure. During Expedition 362, two sites on the Indian oceanic plate ~250 km southwest of the subduction zone, Sites U1480 and U1481, were drilled, cored, and logged to a maximum depth of 1500 meters below seafloor. The succession of sediment/rocks that will develop into the plate boundary detachment and will drive growth of the forearc were sampled, and their progressive mechanical, frictional, and hydrogeological property evolution will be analyzed through postcruise experimental and modeling studies. Large penetration depths with good core recovery and successful wireline logging in the challenging submarine fan materials will enable evaluation of the role of thick sedimentary subduction zone input sections in driving shallow slip and amplifying earthquake and tsunami magnitudes, at the Sunda subduction zone and globally at other subduction zones where submarine fan–influenced sections are being subducted. 
    more » « less
  4. 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. 
    more » « less
  5. Mélange (or block-in-matrix structures) exerts a first-order control on both the mechanical and chemical evolution of subduction megathrusts. However, the timing and mechanisms that form mélanges are variable and debated. Field observations and (micro-) structural analyses from a metasedimentary mélange in the lawsonite blueschist unit of the Catalina Schist (Santa Catalina Island, California, USA) reveal that syn-subduction deformation and fluid-mediated processes led to mélange formation at the plate interface. Deposited as turbidites, early shear occurred parallel to bedding planes (S1 foliation). At near peak subduction conditions, at the base of the subduction seismogenic zone (∼1.0 GPa, 320 °C), the rocks were intensely deformed in recumbent open to tight folds (F2) with axial planar cleavages (S2). Fracturing, fluid flow, and quartz precipitation are preserved as extensional vein mesh networks in fold noses. Continued shearing led to boudinage of these strengthened noses and transformation into strong blocks within the weaker less-veined matrix composed of high-strain fold limbs (S1−2). Microstructures reveal viscous deformation in the high-strain fold limbs occurred by pressure-solution creep of fine-grained quartz ± albite. In contrast, the fold noses and/or blocks contain coarse-grained quartz veins with little evidence of deformation. These rocks record the development of syn-subduction block-in-matrix mélange structures through the interaction of deformation and mineral precipitation; pressure solution weakened fold limbs-turned-matrix and veining strengthened fold noses-turned-blocks. Although mélange structure is often invoked to explain tremor and slow slip, rheological analysis indicates that these metasedimentary rocks can host tectonic creeping but cannot accommodate slow-slip strain rates by the deformation mechanisms preserved in their microstructures. 
    more » « less