An official website of the United States government
Abstract Sediment thermal history controls the progress of diagenetic reactions that can alter the mechanical behavior of material entering a subduction zone that then: accretes to the margin, hosts the plate boundary interface, or is carried deeper within the Earth. On the Cascadia margin offshore Oregon (USA), hydrothermal circulation in the oceanic crust affects thermally controlled processes, enhancing sediment alteration above the MARGIN seamount, which is buried by the Astoria Fan. Hydrothermal circulation increases temperatures at the summit of the seamount and in the overlying sediment by up to ∼100°C. We use sediment thermal history constrained by heat flux observations to model the expected progress of the smectite‐to‐illite reaction around the MARGIN seamount. Above the seamount, the smectite‐to‐illite reaction is expected to progress to completion ∼250 m below the seafloor; away from the seamount, smectite is likely unaltered to a burial depth of ∼800 m. The altered sediment above the seamount has higher rigidity and p‐wave velocity than the surrounding sediment. Spatial variability in sediment alteration may be present around other buried seamounts. We use vertical gravity gradient anomalies to estimate the locations and heights of additional seamounts. Each of these seamounts may have altered sediment around it, which could affect deformation and seismicity in the margin wedge. Because cemented sediment with greater elastic strength is better able to store elastic strain energy, enhanced sediment alteration and cementation above seamounts entering the subduction zone could facilitate earthquake nucleation for material in the margin wedge that was above a seamount prior to subduction.
more »
« less
Volcanic arc rigidity variations illuminated by coseismic deformation of the 2011 Tohoku-oki M9
Rock strength has long been linked to lithospheric deformation and seismicity. However, independent constraints on the related elastic heterogeneity are missing, yet could provide key information for solid Earth dynamics. Using coseismic Global Navigation Satellite Systems (GNSS) data for the 2011 M9 Tohoku-oki earthquake in Japan, we apply an inverse method to infer elastic structure and fault slip simultaneously. We find compliant material beneath the volcanic arc and in the mantle wedge within the partial melt generation zone inferred to lie above ~100 km slab depth. We also identify low-rigidity material closer to the trench matching seismicity patterns, likely associated with accretionary wedge structure. Along with traditional seismic and electromagnetic methods, our approach opens up avenues for multiphysics inversions. Those have the potential to advance earthquake and volcano science, and in particular once expanded to InSAR type constraints, may lead to a better understanding of transient lithospheric deformation across scales.
more »
« less
- Award ID(s):
- 2121666
- PAR ID:
- 10539388
- Publisher / Repository:
- Science
- Date Published:
- Journal Name:
- Science Advances
- Volume:
- 10
- Issue:
- 23
- ISSN:
- 2375-2548
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Lithospheric shortening can be described by one of two end-member modes: indentation of the lithosphere and subduction of the lithospheric mantle. Deciphering the difference between these modes is crucial in the interpretation of past and present orogens and in predicting their structural architecture at depth. It is therefore important to establish how observable upper crustal proxies reflect deep lithospheric kinematics and dynamics. Over the last few decades, geological and geophysical data have provided valuable constraints on the northern margin of the Tibetan Plateau. This margin is defined by the Qilian Shan thrust belt, which developed in response to the far-field convergence between the Indian and Eurasian plates. The primary mechanism for this development is the southward subduction of the Asian lithospheric mantle beneath the Tibetan Plateau. We conducted numerical modelling to simulate the kinematics and response of the upper crust to the southward subduction of the lithospheric mantle. Our results show that subduction of the lithospheric mantle can result in upper crustal deformation that matches the records in the Qilian Shan, where pure shear shortening alone does not generate similar upper crust proxies, including the broad width and architecture of the bivergent orogenic wedge, the timing of fault initiation and evolution, seismicity and fault activity, the topography and geomorphology. The geometry of the subducting lithosphere impacts the width and asymmetry of the bivergent orogenic wedge. Our results demonstrate how records of crustal strain can be used to better interpret the deep structural architecture of past and present orogenic wedges.more » « less
-
Significance The exothermic metamorphic reaction in orthopyroxene (Opx), a major component of oceanic lithospheric mantle, is shown to trigger brittle failure in laboratory deformation experiments under conditions where garnet exsolution takes place. The reaction product is an extremely fine-grained material, forming narrow reaction zones that are mechanically weak, thereby facilitating macroscopic faulting. Oceanic subduction zones are characterized by two separate bands of seismicity, known as the double seismic zone. The upper band of seismicity, located in the oceanic crust, is well explained by dehydration-induced mechanical instability. Our newly discovered metamorphism-induced mechanical instability provides an alternative physical mechanism for earthquakes in the lower band of seismicity (located in the oceanic lithospheric mantle), with no requirement of hydration/dehydration processes.more » « less
-
Abstract Subduction zone accretionary prisms are commonly modeled as elastic structures where permanent deformation is accommodated by faulting and folding of otherwise elastic materials, yet accretionary prisms may exhibit other deformation styles over relatively short time scales. In this study, we use 6.5‐year (2014–2021) Sentinel‐1 interferometric synthetic aperture radar (InSAR) time‐series of post‐seismic deformation in the Makran accretionary prism of southeast Pakistan to characterize non‐linear viscoelastic deformation within an active accretionary prism on short timescales (months to years). We constructed a series of 3‐D finite‐element models of the Makran subduction zone, including an accretionary prism, and constrained the elastic thickness of the upper wedge and the flow‐law parameters (power‐law exponent, activation enthalpy, and pre‐exponential constant) of the lower wedge through forward model fits to the InSAR time‐series. Our results show that the prism is elastically thin (8–12 km) and the non‐linear viscoelastic relaxation of the deep portions of the prism alone can sufficiently explain the post‐seismic surface deformation. Our best fitting flow‐law parameters (n = 3.76 ± 0.39,Q = 82.2 ± 37.73 kJ mol−1, andA = 10−3.36±4.69) are consistent with triggering of low temperature dislocation creep within fluid‐saturated siliciclastic rocks. We believe that the fluids necessary for this weakening originate from sedimentary underplating and/or the presence the hydrocarbons. The presence of power‐law rheology within the lower wedge impacts the estimated plate coupling and the stress state in the subduction system, with respect to the conventional elastic wedge model, and hence should to be considered in future earthquake cycle models.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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1126/sciadv.adl4264
