Fossil Subduction Zone Research Articles

Seismic tomography reveals contrasting styles of subduction-channel and mantle-wedge exhumation controlled by upper plate divergent motion

In fossil subduction zones associated with massive exhumation of (ultra)high-pressure ((U)HP) rocks such as the Western Alps, the geometry and behavior of subduction-channel and mantle-wedge rocks during exhumation are still poorly constrained by independent geophysical observations. Here we use a new local earthquake tomography model of the entire fossil subduction zone of the Western Alps based on data collected during the CIFALPS and CIFALPS2 passive seismic experiments, and the first receiver-function profile across the Ligurian Alps, to investigate the styles of subduction-channel and mantle-wedge exhumation as a function of increasing upper-plate divergence motion. In the northern Western Alps (low divergence), a thickened subduction channel can be detected, but no exhumed mantle wedge is found beneath the Gran Paradiso (U)HP dome. In the southern Western Alps (intermediate divergence), an exhumed mantle wedge is detected beneath the Dora-Maira (U)HP dome above a serpentinized subduction channel. In the Ligurian Alps (high divergence), an exhumed mantle wedge and a former subduction channel are detected at much shallower levels beneath the Voltri-Valosio (U)HP dome, and above a shallow-dipping lower-plate Moho. In this latter case, the lower boundary of the exhumed subduction channel is the most evident seismic-velocity interface, which may be easily misinterpreted as a true Moho. Similar Moho-like interfaces are found beneath the exhumed (U)HP domes of eastern Papua New Guinea and the Dabie Shan, which suggests that the results of the CIFALPS experiments may be used as a reference case for the interpretation of other (U)HP terranes worldwide.

Tourmaline composition probes serpentinite-derived fluid mobility in subduction zones

Serpentinite dehydration in subduction zones plays a pivotal role in geochemical cycling on Earth. A number of geochemical studies on arc magmas have elucidated the contributions of serpentinite-derived fluids to mantle sources. However, due to complex geological overprints during subduction zone processes, discerning serpentinite signatures in exposed metamorphic rocks within fossil subduction zones remains challenging. In this study we address these difficulties through in-situ investigations of tourmaline, the geochemistry of which reflects the host environment as well as potential fluid-induced processes. The presence of zonations in tourmaline makes it an excellent recorder of consecutive geological events. Integrated major and trace elements along with in-situ boron isotopes of tourmaline from the high-pressure Sopron area (Hungary) in the Eastern Alps were used to unravel fluid action sourced from serpentinite. Despite the presence of color zoning, tourmaline in the orthogneiss (Tur-G) has low XMg [Mg/(Mg + Fe)] of ca. 0.3–0.6 and δ11B values of around −11 ‰, along with variable trace element compositions. Petrological observations and geochemical analyses suggest that the inner domains of Tur-G are of igneous origin, while the outer rims are likely affected by subsequent metamorphic events. Tourmaline in metasomatized kyanite-quartzite (Tur-K) veins exhibits distinct geochemical zoning, and preserves metamorphic cores and fluid-induced rims. The inner domains of Tur-K display low XMg (<0.6), relatively high trace element concentrations and δ11B values of less than −10 ‰, whereas the overgrowths exhibit extremely high XMg values (>0.99), low trace element concentrations and high δ11B values reaching up to +21 ‰, clearly indicating the incorporation of serpentinite-derived Mg-11B-rich fluids. Through comparison with other metamorphic and metasomatic tourmalines in (ultra)high-pressure rocks globally, we establish that tourmaline with high XMg > 0.85 and δ11B values >0 ‰ may serve as an effective proxy for detecting serpentinite-derived fluids in subduction zones.

Seismological evidence for a multifault network at the subduction interface

Subduction zones generate the largest earthquakes on Earth, yet their detailed structure, and its influence on seismic and aseismic slip, remains poorly understood. Geological studies of fossil subduction zones characterize the seismogenic interface as a 100 m–1 km thick zone1–3 in which deformation occurs mostly on metres-thick faults1,3–6. Conversely, seismological studies, with their larger spatial coverage and temporal resolution but lower spatial resolution, often image the seismogenic interface as a kilometres-wide band of seismicity7. Thus, how and when these metre-scale structures are active at the seismic-cycle timescale, and what influence they have on deformation is not known. Here we detect these metres-thick faults with seismicity and show their influence on afterslip propagation. Using a local three-dimensional velocity model and dense observations of more than 1,500 double-difference relocated earthquakes in Ecuador, we obtain an exceptionally detailed image of seismicity, showing that earthquakes occur sometimes on a single plane and sometimes on several metres-thick simultaneously active subparallel planes within the plate interface zone. This geometrical complexity affects afterslip propagation, demonstrating the influence of fault continuity and structure on slip at the seismogenic interface. Our findings can therefore help to create more realistic models of earthquake rupture, aseismic slip and earthquake hazard in subduction zones.

Subduction erosion revealed by exhumed lower arc crustal rocks in an accretionary complex, northeastern China

Abstract Subduction erosion at convergent margins is a leading mechanism for the destruction (recycling and reworking) of continental crust. But because of the lack of direct evidence, it is not straightforward to identify erosive events and their intensities in fossil subduction zones. The Heilongjiang accretionary complex in northeastern China was formed during the early Mesozoic subduction of the Paleo-Pacific Ocean. We investigated amphibolites from this accretionary complex, whose protoliths (based on whole-rock trace elements and Sr-Nd-Hf isotopes) were mafic continental arc magmatic rocks (255–249 Ma; zircon core U-Pb ages) from the upper plate. Phase equilibria modeling constrained by mineral geochemistry indicates that the amphibolites and their wall rocks were first heated to low granulite facies (750–800 °C, ~7 kbar) at 251–244 Ma (zircon rim U-Pb ages) and then cooled to ~700 °C with increasing pressure (8–9 kbar) before 213–187 Ma (titanite and apatite U-Pb ages). To explain the occurrence of the lower arc crustal lithologies in the accretionary complex and their metamorphic history, we propose that the subducting plate strongly eroded the forearc crust, allowing the plate interface to advance landward and scrape the amphibolites and wall rocks formed under the old arc, which finally were exhumed along the subduction channel and became components of the complex. The case study exemplifies direct petrological evidence of strong subduction erosion occurring in an ancient orogen, thus implying that consumption of the entire forearc crust could occur within only ~50 m.y.

Crustal‐Scale Seismic Reflection Profiling Constrains How the Paleo‐Asian Ocean Was Closed

AbstractThe Central Asian Orogenic Belt (CAOB) is the most significant accretionary orogenic belt since the Phanerozoic and the most ideal site for studying continental growth evolution processes. A 460‐km‐long high‐resolution crustal‐scale seismic reflection study was conducted across the eastern CAOB in North‐Central China to constrain the closure mode and location of the Paleo‐Asian Ocean, that is, the previous ocean of the CAOB. The resultant seismic reflection profile revealed opposite‐dipping reflectors in the northern and southern parts of the profile, which converge at the profile center to form an inverted U‐shaped reflector pattern near the crust–mantle transition zone beneath the Solonker Suture. The dipping reflectors represent bidirectional fossil subduction zones sloping to the north and south, and the convergence reflector pattern represents the ocean closure location. Integration of these results with available geological data facilitated model construction whereby Paleo‐Asian Ocean closure was accomplished by divergent subduction of the Paleo‐Asian oceanic plate, with northward subduction beneath the southern margin of the Mongolian Block and southward subduction beneath the northern margin of the North China Craton. The oceanic lithosphere contracted and deformed, yielding the observed inverted U‐shaped reflector pattern, representing Paleo‐Asian Ocean closure. This subsurface location lies beneath the Solonker Suture surface exposure, suggesting that this suture marks the ocean closure location, rather than the previously proposed Hegenshan–Heihe Suture to the north or Xar Moron Suture to the south. Our study suggests that divergently dipping subduction and associated accretion and magmatism may constitute the primary continental growth mode for accretionary‐type orogens.

Eastward extension of the Solonker Suture beneath the central Songliao Basin, NE China: Evidence from a deep seismic reflection profile

The Solonker Suture is one of the most important structures in central Asia and marks the location of closure of the Paleo-Asian Ocean and the assemblage of the North China and Siberia cratons during the Late Permian - Early Triassic. The suture is traceable for ∼ 700 km within the Central Asian Orogenic Belt, but its eastward extension beneath the Songliao Basin has not been defined. In this paper, we present a 127.3-km-long N–S deep seismic reflection profile across the central Songliao Basin. The data acquired from five large-size and four medium-size dynamite shots were used to generate a single-fold reflection profile imaging the Moho and the structures of the crust and upper mantle. The southern end of the profile reveals a northward-dipping reflector that extends 80 km into the mantle. This dipping reflector is interpreted as a fossil subduction zone. The new data from this study suggest that the Songliao Terrane may not be a single micro-continent block or orogenic belt but a composite terrane that comprises the North Accretionary Zone, South Accretionary Zone, and northern margin of the North China Craton from north to south. The fossil subduction zone is regarded as part of the NE-trending extension of the Solonker Suture beneath the central Songliao Basin.

Fossil divergent double-subduction zone in the Great Xing’an Range, NE China: Evidence from a deep seismic reflection profile

Divergent double-subduction is an important component of plate tectonics and continental growth at modern active convergent margins. However, such a process is rarely identified by deep seismic reflection (DSR) profiles, which casts doubt on its past significance. We report here a 416 km DSR profile of large dynamite shots across the Great Xing’an Range (NE China), which shows multiple arc-shaped reflections in the upper-middle crust of the Xing’an terrane, and densely layered reflection structures suggesting tectonic nappes in the lower crust. This new evidence, together with recent geological studies, indicates that the Xing’an terrane may not be a micro-continental block as previously thought, but a Paleozoic accretionary complex. In addition, there is a strong west-dipping reflection at the boundary between the Xing’an and Erguna terranes in the western DSR profile extending from the lower crust to the lithospheric mantle, suggesting that this fossil subduction zone is most likely the Toudaoqiao-Xinlin suture, which separates the Xing’an and Erguna terranes. In addition, an east-dipping strong reflection was also found in the eastern profile between the Xing’an and Songliao terranes, corresponding to the Hegenshan-Nenjiang suture at the surface. Given the newly identified crust-mantle structure, the opposite subduction zone, and other geological evidence, we propose a new model that the Xing’an terrane was developed by divergent double-subduction of the Paleo-Xing’an Ocean (part of the Paleo-Asian Ocean in the area) from ca. 500 Ma to ca. 300 Ma. Our study provides the first DSR evidence for a divergent double-subduction in the Paleozoic, which emphasizes its importance in past accretionary orogens.

Constraints on tectonic processes in subduction mélange: A review of insights from the Catalina Schist (CA, USA)

Subduction mélange, with its distinctive block-in-matrix structure, is documented in exhumed fossil subduction zones worldwide. Rocks from these terranes preserve features that record tectonic processes from the time that the rocks were at the subduction interface. Careful study of these features allows for connections to be made with tectonic processes occurring in active subduction zones. The Catalina Schist mélange has served as an exhumed analog in such studies as it records abundant evidence for tectonic processes occurring at the subduction interface. Focusing of fluids in mélange matrix at the subduction interface is documented, and this fluid-rich environment may have contributed to seismic activity. Deformation and tectonic mixing juxtaposed disparate materials (mafic, ultramafic, and sedimentary rocks) over length-scales of 10s of km along the interface, occurring in concert with metasomatism and mass transport by fluids to create mineralogically, chemically, and rheologically distinct compositions. These processes may have impacted seismic behavior and plate geodynamics along with influencing the chemistry of arc magmas that form above the subduction interface. Evidence suggests that the duration of tectonic formation of mélange may be variable from one locality to another, with relatively small differences in peak ages of blocks of <10 Myrs recorded in the amphibolite facies rocks of the Catalina Schist.

Effects of Tectonic Quiescence Between Orogeny and Rifting

AbstractTectonic quiescence is a period of nearly null relative plate motion. During the Atlantic Ocean opening, rifting cut across continental segments of supercontinent Pangea that had been through different periods of quiescence, which varied from tens to hundreds of Myr. Using numerical models, we explored the effects of quiescence between orogeny and rifting on the pre‐rift lithosphere and subsequent conjugate rifted margin configuration. Our model results showed that quiescence of ∼30–60 Myr caused wide orogenic wedges (&gt;300 km) to be heated beneath their centers, which led to the breakup of the lithospheric mantle before the breakup of the continental crust during rifting, resulting in hyperextended conjugate rifted margins. In these scenarios, continental lithosphere breakup away from the suture contributed to the preservation of previously subducted lower crust along the subduction zone. Long periods of tectonic quiescence (∼100 Myr) experienced by a wide orogen permitted efficient orogenic collapse and lateral spreading, with pronounced lithospheric thermal weakening. Consequently, post‐quiescence rifting reactivated the suture and effectively drove plate eduction, exhuming most of the subducted continental lower crust and mantle lithosphere. Ultimately, ultra‐wide domains developed into asymmetric conjugate margins. Contrastingly, in narrow orogenic wedges (&lt;200 km), different quiescence durations had minor effects on the margin architecture. In these cases, a longer quiescence time contributed to a larger preservation of subducted crust in the fossil subduction zone after rifting. Our models agree with conjugate rifted margins in the North and South Atlantic, which probably developed after the collapse of wide orogens and subsequent periods of tectonic quiescence.

Deformation Mechanisms of Blueschist Facies Continental Metasediments May Offer Insights Into Deep Episodic Tremor and Slow Slip Events

AbstractExhumed fossil subduction zones are archives of the deformation conditions and mechanisms from depths not directly accessible. Microstructural analysis of samples exhumed therefrom offers insights into the micromechanics and deformation processes associated with subduction such as earthquakes, slow earthquakes, and aseismic creep. Subducted and exhumed continental metasediments of the Tuscan Metamorphic Units of the Italian Northern Apennines contain a mylonitic foliation and quartz and carpholite dilational hydroshear veins with crack‐and‐seal textures, both developed at blueschist facies conditions (350–400°C, ∼1 GPa). As shown by overprinting relationships and mineral assemblages, these two structure types formed broadly coeval within the stability field of carpholite. Metaconglomerates and metaquartzarenites deformed mainly by dissolution‐precipitation creep and secondary by dislocation creep. Microstructural and electron backscatter diffraction analyses of the veins suggest only limited recrystallization of quartz fibers by subgrain rotation recrystallization, with adjacent metapelite bands acting as decollement horizons, likely by slip on the basal plane of phyllosilicates. We estimated differential stresses of 43–55 MPa and strain rates between 10−13 and 10−14 s−1 from the vein recrystallized quartz fibers. We propose these microstructures and deformation mechanisms to represent a geological evidence of deep episodic tremor and slow slip events in subducted continental metasediments. Pore pressure cyclically reached supralithostatic values triggering tremors causing fracturing of all involved lithotypes. Likely, slow slip was accommodated preferentially by slip on phyllosilicate bands. Aseismic creep occurred mainly by dislocation creep with subgrain rotation recrystallization in vein quartz, slip on the basal plane of phyllosilicates, and dissolution and precipitation creep in the host rock.

Along-strike variations in the fossil subduction zone of the Western Alps revealed by the CIFALPS seismic experiments and their implications for exhumation of (ultra-) high-pressure rocks

In complex plate-boundary settings, a reliable 3-D geophysical characterization of the deep tectonic structure is a fundamental starting point for a breakthrough in the analysis of processes controlling plate subduction and (U)HP rock exhumation. The Western Alps host one of the best-studied fossil subduction zones worldwide, with a well-defined deep structure in 2-D based on recent geophysical experiments. However, a full 3-D characterization of its deep tectonic structure is still lacking. Here we present a series of new receiver function cross-sections across the northern and southern Western Alps, validated and complemented by a S-wave velocity model from ambient-noise tomography that provides additional constraints between the profiles. We document a marked change in Moho attitude from the northern Western Alps, where the eastward-dipping European Moho reaches ∼45 km depth beneath the Gran Paradiso dome, to the southern Western Alps, where the European Moho reaches ∼70 km depth beneath the equivalent Dora-Maira dome. This change in Moho attitude takes place over a few tens of kilometers and was likely emphasized by deformation of the slab during subduction. The Western Alps subduction wedge is much thicker in the south than in the north, and the mantle-wedge rocks are deeply involved in orogeny exclusively in the south, where coesite is found in continental (U)HP rocks at several locations. Our detailed information on the 3-D structure of the subduction wedge provides first-order constraints for the next-generation of thermo-mechanical numerical models and may help explain the lateral variations in exhumation style revealed by the geologic record.

Early Late Triassic retro-foreland basin in response to flat subduction of the Paleo-Tethyan oceanic plate, SE Tibet

Syn-subduction basins bear significant implications to understand tectonic evolution of any fossil subduction zone. The late Paleozoic to early Mesozoic (Paleo-Tethyan) tectonics of the eastern and southeastern Tibetan Plateau (i.e., the Sanjiang Orogenic Belt) is featured by ocean-continent subduction systems. A huge pile of volcanic-absent sedimentary succession developed in the middle segment of the Sanjiang orogenic belt, its age and tectonic nature remain unclear. Detailed geological mapping and zircon U-Pb dating results demonstrate that the early Late Triassic volcanic-absent succession comprises the nonmarine Maichuqing Formation in the lower part and the shallow marine Sanhedong Formation in the upper part. The Maichuqing Formation consists of coarse to fine-grained sandstone, siltstone and mudstone with abundant basal erosional surfaces, trough and planar cross-beddings, ripples, mudcracks, and plant fragments. The Sanhedong Formation comprises predominantly bioclastic limestones interlayered with marl, calcareous-muddy siltstone, and calcareous sandstone with abundant bivalve fossils. Syn-sedimentation deformation structures, such as slump folds and associated normal faults are common, suggesting intense tectonism during deposition. Synthesizing sedimentary data, paleocurrent and provenance results, combined with other available data, demonstrate that the volcanic-absent succession deposited within a retro-foreland basin along the rear part of the Permian-Triassic Jomda-Weixi-Yunxian arc in response to flat-subduction of the Paleo-Tethyan Ocean during the early Late Triassic time.

Multiple volcanic episodes of the Kermanshah forearc basin, SW Iran: a record of the deactivation and re-initiation of Neotethyan subduction involving a mid-ocean ridge

Numerical simulations and theoretical analyses predict that forearc ophiolites probably record the end of subduction followed by re-initiation. We report here the results of a field study and new geochronological data that enable identification of the Kermanshah forearc basin, SW Iran. This basin is filled with two terrestrial volcaniclastic successions separated by an Upper Cretaceous reef facies limestone. Early Late Cretaceous ocean island basalts, early Eocene incompatible element-enriched mid-ocean ridge basalts/normal-type mid-ocean ridge basalts and late Eocene arc basalts of the basin indicate that Neotethyan subduction ended during the early Late Cretaceous and then re-initiated during the late Eocene. We suggest that the Neotethyan Ocean between the Iranian and Arabian continents consisted of leading and trailing oceanic plates separated by a mid-ocean ridge (MOR). Subduction of the leading plate generated the Jurassic–Early Cretaceous Sanandaj–Sirjan arc, which was ended by MOR–trench collision. The Zagros late Early Cretaceous forearc ophiolites represent the youngest/hottest segments of the leading plate emplaced during the MOR–trench collision. Subsequently, subduction of the MOR generated the forearc basin and basalts with diverse geochemistry. After c. 60 myr, the cooled trailing plate started to subduct, generating the late Eocene arc basalts. This model has significant implications for investigating and understanding other fossil subduction zones elsewhere in the world. Supplementary material: Supplementary material is available at https://doi.org/10.6084/m9.figshare.c.6086586

The effect of supercritical fluids on Nb-Ta fractionation in subduction zones: Geochemical insights from a coesite-bearing eclogite-vein system

Supercritical fluids are ideal media for mass transfer from the subducting slab into the mantle wedge. However, little is known about the role of natural supercritical fluids in subduction zones. A combined study of petrology, geochemistry and zirconology was carried out for a coesite-bearing eclogitic vein and its surrounding eclogites from the Dabie orogen, one of the typical continental subduction zones on Earth. The results are used to reveal the composition and source of supercritical fluids, the P-T conditions and timing of their formation and their geochemical effects in subduction zones. The eclogitic vein is composed of garnet, omphacite, quartz, amphibole, rutile and apatite, with minor amounts of epidote, plagioclase and zircon. Coesite was identified in omphacite in the vein, demonstrating its formation under ultrahigh-pressure (UHP) metamorphic conditions that correspond to subarc depths in oceanic subduction zones. Zircons from the vein give concordant U-Pb ages of 225 ± 5 Ma and flat HREE patterns for newly grown rims, consistent with their crystallization at the UHP eclogite facies. The vein shows similar Hf-O isotope compositions to the host UHP eclogites, indicating that the UHP vein-forming fluid is internally buffered within the UHP eclogites. Minerals in the vein contain not only higher contents of Cr, Ni, Sr, REE and HFSE, but also multiphase crystal inclusions (such as omphacite, quartz, epidote, apatite, amphibole, plagioclase, mica, rutile, calcite, and anhydrite) as well as liquid and gas phases of H2O. The calculated fluid composition of major elements is 31 wt% SiO2, 20 wt% CaO, 10 wt% (SO4)2−, 8 wt% (CO3)2−, 8 wt% FeO, 7 wt% Al2O3, and 9 wt% H2O, with traces of Na2O, K2O and TiO2. This composition is responsible for the UHP veining from supercritical silicate-rich fluids in equilibrium with peak UHP minerals. The P-T-t path of the UHP eclogite-vein system indicates that the second critical endpoint of the basalt-H2O system is located close to 3.4 GPa and 770°C. The high contents of sulfate and carbonate in the vein suggest that supercritical fluids were oxidized and likely a predominant agent to transfer sulfur and carbon at subarc depths. High Nb/Ta ratios of rutile crystals in the vein and quantitative modelling indicate that dehydration of subducting eclogites at subarc depths can produce supercritical fluids with suprachondritic Nb/Ta ratios. The UHP eclogites and mantle wedge peridotites metasomatized by supercritical fluids can acquire suprachondritic Nb/Ta ratios and thus provide a complementary reservoir to balance the subchondritic reservoirs on Earth. Based on a statistics of Nb-Ta data, two indices are proposed to identify the existence of supercritical fluids in fossil subduction zones: (1) rutile crystallized from supercritical fluids shows lower Nb contents than that from aqueous solutions and hydrous melts; (2) UHP eclogites associated with supercritical fluids exhibit a distance >0.1 in their Nb-Ta compositions to the line defined by basalts in the plot of log[Nb] vs log[Ta]. Therefore, the existence of supercritical fluids during UHP metamorphism at subarc depths can be recognized from their effect on the mobility of fluid-immobile trace elements in UHP metamorphic rocks.

Magnetic Structure and its Tectonic Implication Around Longmenshan Fault Zone Revealed by EMAG2v3

The formation of magnetic minerals is bound up with the tectonic evolution history, whereupon the distribution of magnetic anomalies has great meanings for regional tectonics. In this study, we use the latest global lithospheric magnetic field model EMAG2-v3, processed by various techniques including reduction to the pole (RTP), upward continuation, derivations, Euler deconvolution, estimation of total magnetization direction, and Curie point depth (CPD), to unveil the tectonics around Longmenshan fault zone (LFZ). LFZ is clearly displayed as a positive and negative anomaly transition zone in RTP anomalies and acts as a magnetic basement boundary. The Sichuan Basin (SB), located to the east of LFZ, is marked by strong magnetic basement and NE-strike banded induced positive anomalies which are associated with the Neoproterozoic magmatic activity. The banded shape, absence of radial pattern of anomalies, and existence of fossil subduction zone supports that the magnetic basement was formed in arc environment. The CPD in SB estimated by radial average power spectral is 30–51 km, which allows magnetic minerals in deep crust or even in lithospheric mantle to exhibit high magnetizations. The Songpan-Ganzi fold belt (SGFB), in contrast, is located to the west of LFZ and covered by thick and low-susceptibility Triassic deposits of flysch, manifests as weak negative anomalies caused by relatively shallow CPD and widespread remanent magnetization. Significant positive anomalies, appearing around the Manai and Rilonguan granitic massifs, indicate a strong magnetic basement beneath SGFB, which is conjectured as westward extension of the Yangtze Block at depth.

Subducted fragments of the Liguro-Piemont ocean, Western Alps: Spatial correlations and offscraping mechanisms during subduction

Fragments of subducted slow-spreading oceanic lithosphere are exposed continuously in the Liguro-Piemont domain of the Western Alps. By combining new and literature petrological data, interpolated maps of maximum temperatures, maximum Si contents of phengite as a proxy for peak pressure and thermodynamic modelling, we provide a detailed framework of the peak metamorphic conditions experienced by the distinct subduction slices. High-resolution mapping confirms the marked eastward increase in metamorphic grade throughout the domain, as well as within some slices. The compilation of lithostratigraphic, structural and radiochronological data and the estimation of sediment/mafic-ultramafic ratio for each slice allow refining the origin of these tectonometamorphic units within the former oceanic domain. The refined structural sketchmap allows to restore the geometries of the Alpine subduction at peak burial conditions. Results point to a trimodal distribution of the units with an increase in metamorphic conditions from the Upper (LPU; 320-400°C- 1.2-1.9 GPa) to the Middle (LPM; 415-475°C- 1.7-2.2 GPa) and to the Lower units (LPL; 500-580°C- 2.2-2.8 GPa). The blueschist-facies LPU and LPM units are dominated by sediments (>90%), whereas the eclogitic LPL units are far richer in mafic-ultramafic rocks (>40%). These characteristics, along with lithostratigraphic differences, reflect major differences in their initial paleogeography and/or in the mechanisms responsible for material offscraping from the downgoing slab. The peak burial depths of the LPU, LPM and LPL units are similar to those inferred for slicing and underplating in both modern and fossil subduction zones. Petrological and lithostratigraphic data suggest that the offscraping of the LPU and LPM units was mostly controlled by lithological contrasts, within pelagic shales or along contacts within the uppermost serpentinized mantle. In contrast, major dehydration reactions (such as lawsonite breakdown in sediments) likely controlled the offscraping of the LPL units at eclogite-facies conditions, possibly through high fluid pressure conditions and rocks embrittlement.

Regional and Local Patterns of Upper‐Plate Deformation in Cascadia: The Importance of the Down‐Dip Extent of Locking Relative to Upper‐Plate Strength Contrasts

AbstractThe details of subduction zone locking place constraints on the characteristics of megathrust events. Due to the lack of significant present‐day seismicity along the Cascadia subduction interface, geodetic data are used to assess subduction locking along the margin. We isolate the subduction signal from other tectonic signals within the Cascadia GPS field, to assess the details of plate‐interface locking. Apparent coupling determined by a simple homogenous elastic half‐space inversion cannot everywhere reproduce the subduction component of the GPS field. Consequently, we explore the relationships among upper‐plate strength, locking depth and the resulting surface velocity signal using 2D finite element models. When the upper plate is composed of a weak material, trenchward of a strong backstop, we find that the down‐dip limit of locking relative to the location of the weak‐to‐strong transition controls how upper‐plate deformation is spatially distributed. If locking extends into the stronger material, as observed in central Cascadia, the surface velocity field propagates farther inland than expected from a simple homogeneous elastic model. In contrast, in southern Cascadia, because locking terminates within the weak accretionary margin, upper‐plate shortening is localized within the weaker material, particularly in the region between the end of locking and the strong Klamath terrane. This behavior provides a possible mechanism for producing the high (geodetic and permanent) uplift rates, plate‐motion‐parallel shortening, and crustal exhumation observed in many active and fossil subduction zones.

Deep Sourced Fluids for Peridotite Carbonation in the Shallow Mantle Wedge of a Fossil Subduction Zone: Sr and C Isotope Profiles of OmanDP Hole BT1B

AbstractCompletely carbonated peridotites represent a window to study reactions of carbon‐rich fluids with mantle rocks. Here, we present details on the carbonation history of listvenites close to the basal thrust in the Samail ophiolite. We use samples from Oman Drilling Project Hole BT1B, which provides a continuous record of lithologic transitions, as well as outcrop samples from listvenites, metasediments, and metamafics below the basal thrust of the ophiolite. 87Sr/86Sr of listvenites and serpentinites, ranging from 0.7090 to 0.7145, are significantly more radiogenic than mantle values, Cretaceous seawater, and other peridotite hosted carbonates in Oman. The Hawasina sediments that underlie the ophiolite, on the other hand, show higher 87Sr/86Sr values of up to 0.7241. δ13C values of total carbon in the listvenites and serpentinites range from −10.6‰ to 1.92‰. We also identified a small organic carbon component with δ13C as low as −27‰. Based on these results, we propose that during subduction at temperatures above &gt;400°C, carbon‐rich fluids derived from decarbonation of the underlying sediments migrated updip and generated the radiogenic 87Sr/86Sr signature and the fractionated δ13C values of the serpentinites and listvenites in core BT1B.

Stress, microstructure, and deformation mechanisms during subduction underplating at the depth of tremor and slow slip, Franciscan Complex, northern California

Subduction zone shear stress is thought to influence the topography of accretionary wedges as well as earthquake magnitude and frequency, and it has traditionally been calculated indirectly using methods that include relating surface heat flow and wedge topography to shear stress on the megathrust. The Eastern Belt of the Franciscan Complex in northern California is part of an exhumed accretionary complex deformed and metamorphosed under high P/T ratio conditions. Its lack of a retrogressive overprint makes it an ideal location to examine the mechanics of accretion within a fossil subduction zone. Blueschist-facies metasediments were deformed by a combination of dislocation and pressure-solution creep. EBSD analysis of dynamically recrystallized quartz across both the upper and lower boundaries of an accreted packet of metasediments indicates shear stresses of 32.8–46.6 MPa along the bounding faults and ∼26 MPa within the main body of the accreted packet.A novel method of determining minimum strain due to pressure solution yields finite stretches of ∼1.6–1.9 and corresponds to strain rates generally on the order of 10−14 s−1 over the 2–6 m.y. period of deformation. A maximum strain rate on the order of 10−13 s−1 was calculated based on plate convergence rate and shear zone width. These independently determined strain rates are shown to be best predicted by a thin film pressure solution flow law. Deformation mechanism maps show that pressure solution and dislocation creep operated with similar strain rates, and that progressive differentiation then increased the favorability and domination of dislocation creep.Dilational microcracking associated with pressure solution produced opening increments in crenulation arcs of ∼ 1 mm. These events have the potential to produce microseismic events with an effective seismic moment comparable to that of low-frequency earthquakes. Concatenation of these events to form shear zones or folds on scales up to hundreds of meters could explain tremor bursts and slow slip events.

Abiotic methane generation through reduction of serpentinite-hosted dolomite: Implications for carbon mobility in subduction zones

Abiotic methane has been increasingly detected at the surface of Earth and other terrestrial planets, exerting a strong effect on the study of chemolithoautotrophic life and thus astrobiology. In contrast, abiotic methane generation in subduction zones, which is intimately linked to questions such as the mechanisms of deep carbon mobility, has received scarce attention. Experiments elucidated the significant production of abiotic methane through reduction of carbonate minerals under subduction zone conditions, whereas detailed geological conditions and processes for the reduction in natural rocks are hitherto poorly understood. Here, we report carbonate reduction and genesis of abiotic methane in dolomitized serpentinites (referred to as ophidolomites) from a fossil subduction zone (SW Tianshan, China). Detailed petrological, Raman spectroscopic, strontium and carbon isotopic, and thermodynamic results provide evidence for dolomite reduction into the phase assemblage of calcite + brucite + methane, likely associated with retrograde serpentinization starting at 7–9 kbar and 410–430 °C in the subduction zone. Microthermometric data for dolomite-hosted fluid inclusions are consistent with petrographic observations, indicative of fluid entrapment postdating the onset of dolomite reduction during exhumation. Model calculations suggest that water-rich fluids characterized by relatively high hydrogen fugacities can create favorable conditions for the reduction process, which, however, do not exclude the possibility of carbonate methanation by hydrogen-rich fluids as reported in previous studies. The widespread occurrence of methane in these rocks gives credence to the intricate redox transformations of subducted carbon, implying that the elevated hydrogen fugacities may facilitate abiotic synthesis of methane through dolomite reduction at convergent plate boundaries. Our work shows that alteration of dolomite-bearing lithologies represents a potential source for abiotic methane in subduction zones, which may have implications for the transfer of subducted carbon.