High-pressure synthesis and storage of solid organic
compounds in active subduction zones

Baptiste Debret, Bénédicte Ménez, Bastien Walter, Hélène Bouquerel, Pierre Bouilhol, Nadine Mattielli, Céline Pisapia, Thomas Rigaudier, Helen Myfanwy Williams

► To cite this version:

Baptiste Debret, Bénédicte Ménez, Bastien Walter, Hélène Bouquerel, Pierre Bouilhol, et al.. High-pressure synthesis and storage of solid organic compounds in active subduction zones. Science Advances, 2022, 8 (37), 10.1126/sciadv.abo2397. hal-04053952

HAL Id: hal-04053952

https://hal.science/hal-04053952

Submitted on 1 Apr 2023

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# High-pressure synthesis and storage of solid organic compounds in active subduction zones

Baptiste Debret<fnref n="1" />*, Bénédicte Ménez<fnref n="1" />, Bastien Walter<fnref n="2" />, Hélène Bouquerel<fnref n="1" />, Pierre Bouilhol³, Nadine Mattielli⁴, Céline Pisapia<fnref n="1" />, Thomas Rigaudier³, Helen Myfanwy Williams⁵

Recent thermodynamic and experimental studies have suggested that volatile organic compounds (e.g., methane, formate, and acetate) can be produced and stabilized in subduction zones, potentially playing an important role in the deep carbon cycle. However, field evidence for the high-pressure production and storage of solid organic compounds is missing. Here, we examine forearc serpentinite clasts recovered by drilling mud volcanoes above the Mariana subduction zone. Notable correlations between carbon and iron stable-isotope signatures and fluid-mobile element (B, As and Sb) concentrations provide evidence for the percolation of slab-derived CO₂-rich aqueous fluids through the forearc mantle. The presence of carbonaceous matter rich in aliphatic moieties within high-temperature clasts (>350°C) demonstrates that molecular hydrogen production associated with forearc serpentinization is an efficient mechanism for the reduction and conversion of slab-derived CO₂-rich fluids into solid organic compounds. These findings emphasize the need to consider the forearc mantle as an important reservoir of organic carbon on Earth.

Copyright © 2022
The Authors, some
rights reserved;
exclusive licensee
American Association
for the Advancement
of Science. No claim to
original U.S. Government
Works. Distributed
under a Creative
Commons Attribution
License 4.0 (CC BY).

## INTRODUCTION

In modern Earth, the extraction of carbon from the mantle to the surface plays a major role regulating the temperature, composition, and oxygenation levels of the atmosphere and hydrosphere (1-3). Mantle carbon extraction mainly occurs at mid-ocean ridges and subduction zones through mantle partial melting and magmatic degassing processes (4, 5). Carbon is then returned to the mantle via subduction zones, although the nature and the fraction of carbon that is recycled into Earth's interior remains controversial (6-8). The sedimentary, crustal, and mantle parts of the subducting slab and the mantle wedge are all potential hosts for inorganic and organic carbon in subduction zones. However, while many studies have provided estimates of the budgets and distribution of inorganic and organic carbon in the oceanic crust and its sedimentary cover (9-13), relatively little is known about the role that mantle rocks may play in the subduction zone carbon cycle, especially regarding organic carbon.

At Earth's surface, the hydration of mantle-derived ultramafic rocks (a series of complex reactions encompassed under the term serpentinization) plays a major role in the conversion of inorganic carbon, such as CO₂ or CO, into organic compounds, made of carbon and hydrogen often associated with heteroatoms, such as oxygen or nitrogen, via abiotic redox reactions (14). In these reactions, molecular hydrogen (H₂) is produced from water reduction predominantly associated with the oxidation of primary Fe<fnref n="2" />⁺-bearing minerals to form Fe³⁺-bearing minerals, such as magnetite at high temperatures (>200°C) or serpentine at low temperatures (<200°C) (15, 16). The reducing power of H₂ is widely recognized as a source for the abiotic production of hydrocarbons (CₙHₘ), including methane (CH₄) and volatile organic acids (e.g., formic acid) detected in hydrothermal fluids [e.g., (17)].

Recent advances in in situ spectroscopic techniques have allowed the identification of serpentinization-derived solid organic compounds in abyssal serpentinites (18-20). These results provide fundamental evidence for the abiotic synthesis and storage of a diverse range of organic compounds in these rocks, including carbonaceous matter related to polycyclic aromatic hydrocarbons (PAHs). It is now recognized that there is a broad range of environments across Earth's surface and interior that may support conditions suitable for the abiotic synthesis of both solid and volatile organic compounds and that these processes are not limited to low-pressure (LP) hydrothermal environments. In particular, solid carbon compounds represent a host for carbon in metamorphic rocks and may represent an important, although largely unconstrained, fraction of carbon in the lithosphere, potentially contributing to the deep carbon cycle on Earth. Recent thermodynamic (21) and experimental (22, 23) studies show that the formation and stability of organic compounds in metamorphic fluids are predominantly controlled by variations in oxygen fugacity (fO₂) and pH, as opposed to pressure and temperature. Abiotic redox reactions leading to the synthesis of diverse organic components may therefore occur at depth in subduction zones, despite the extreme temperatures and pressures.

Although organic compounds (e.g., CH₄, formate, and acetate) are predicted to be present as dissolved species in slab-derived fluids (21), no consensus for the formation, role, and variability of solid organic compounds in subduction zones exists. At this time, few natural occurrences of organic compound formation in high-pressure (HP) metamorphic fluids have been documented, although exclusively limited to orogenic settings. For example, the observation of H₂- and/or CH₄-rich fluid inclusions in HP ultramafic rocks from the Western Alps (France, Italy) and Vermont (USA) mountain belts provides pieces of evidence for the circulation of reduced and organic-rich fluids in paleosubduction zones (24-26). Vitale-Brovarone et al. (24) proposed that the production of H₂-rich fluids during mantle peridotite serpentinization at HP and their interaction with carbonate minerals could generate abiotic CH₄. Under favorable conditions, these H₂-CH₄ fluids may also favor the abiotic precipitation of abiotic graphite-bearing lithologies as

<footnote n="1">¹Université Paris Cité, Institut de physique du globe de Paris, CNRS, Paris, France.</footnote>

<footnote n="2">²Université de Lorraine, CNRS, GeoRessources, Vandoeuvre-lès-Nancy, France. ³Université de Lorraine, CNRS, CRPG, Nancy, France. ⁴Laboratoire G-Time, DGES, Université Libre de Bruxelles (ULB), Brussels, Belgium. ⁵Department of Earth Sciences, University of Cambridge, Cambridge, UK.</footnote>

<footnote>*Corresponding author. Email: debret@ipgp.fr</footnote>
they migrate through the mantle (27, 28). In addition, Bouilhol et al. (29) observed solid organic compounds with variable relative content in aliphatic and aromatic moieties in HP-ultramafic rocks from eclogite-facies Alpine meta-ophiolites (550° to 600°C, 2 to 2.5 GPa), which demonstrates that organic compounds are not limited to fluids but can also exist and persist as rock-hosted organic matter in HP environments. These observations are consistent with thermodynamic (30) and experimental (31) studies that have shown a progressive increase in PAH stability in rocks relative to fluids with increasing pressure. Furthermore, at HP, the formation of PAHs requires lower activity of carbon and H₂ relative to LP environments (30). It is therefore conceivable that abiotic synthesis of solid organic compounds could be highly efficient in extreme environments, such as subduction zones, favoring the storage of organic carbon in rocks relative to fluids and potentially providing an important global reservoir for carbon. Critically, the formation and stabilization of these solid organic compounds will exert a major control on subduction zone redox reactions and fO₂, as for example, the conversion of carboxylic acids into methane requires an addition of C—H or H₂ (hydrogenation) and/or a loss of oxygen through redox processes. However, although solid organic compounds have been observed in active subduction settings, their presence has been ascribed to microbial activity associated with LP retrograde processes rather than abiotic reactions (32). Furthermore, there are almost no constraints on the nature and abundance of solid organic compounds generated at HP through abiotic redox processes in subduction zones and their potential role in the deep carbon cycle. Here, we address this knowledge gap and revisit the abiotic and biological formation of solid

organic compounds in the active Mariana subduction system through a study of serpentinized clasts recovered during the Integrated Ocean Discovery Program (IODP) Expedition 366 from mud volcanoes formed above the Mariana forearc. We demonstrate that the abiotic synthesis of solid organic compounds in the forearc mantle wedge results from the reduction of CO₂-rich slab-derived fluids at HP and high temperatures (>350°C). We show that although solid inorganic carbonates are abundant near the surface, they relate to shallow microbial activity while, with increasing depth, rock-hosted solid organic compounds become the main carrier of carbon, calling for a reappraisal of serpentinized forearc as a major reservoir for abiotic production and storage of organic carbon in subduction systems.

## RESULTS

### Geological setting and sample selection

The Mariana subduction zone belongs to a nonaccretionary subduction system, spreading from near Tokyo (Japan) to beyond Guam (USA) and involving the subduction of the Mesozoic Pacific plate below the Philippine Sea plate (Fig. 1A). It is the only place in the world where active metasomatized forearc can be sampled (33–35). Along the Mariana forearc, serpentinite mud volcanoes erupt as a consequence of the interaction of aqueous fluids released from the dehydrating subducting slab with forearc mantle wedge peridotites. These muds bring to the seafloor a cargo of serpentinized forearc-derived peridotite clasts that provide a direct window on chemical exchange between the slab and the forearc mantle at various depths. The clast samples studied here were recovered by drilling three mud

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Bathymetry map, pressure and temperature (P-T) path and conceptual cross section of the Mariana forearc system. Figure modified after (36, 37). (A) The map shows the locations of the three mud volcanoes, Yinzao, Fantangisña, and Asút Tesoru, drilled during IODP Expedition 366. Hole locations are indicated in red circles. (B) P-T path presenting the different conditions of serpentinization recorded by the sampled forearc ultramafic clasts (i.e., Liz-, Atg-, and blue-serpentinites, from 1 to 3, respectively). Life limit is after Kashefi and Lovley (46). (C) Conceptual model illustrating serpentinization processes in relation to fluid circulation and mantle flow within the Mariana forearc.</figcaption></figure>
volcanoes, namely, Yinazao, Fantangisña, and Asùt Tesoru, located at increasing distance from the trench (Fig. 1A). Three types of serpentinized forearc peridotite clasts, namely, lizardite, antigorite, and blue serpentine-bearing clasts, were identified on the basis of their petrographic and geochemical features (35–37). These clasts record various stages of serpentinization from deep forearc mantle hydration by slab-derived fluids to shallow fluid-rock interaction during clast exhumation nearby the seafloor (Fig. 1, B and C).

The lizardite- and antigorite-bearing serpentinite clasts preserved a record of deep forearc mantle serpentinization (Fig. 1, B and C) by slab-derived fluids and do not display evidence for notable retrogression (e.g., blue-serpentine or sulfide crystallization, late-stage carbonation) (36). The lizardite-bearing serpentinites (Liz-serpentinites) display variable degrees of serpentinization (30 to 100%, as defined by petrographic observations and bulk rock geochemical analyses) and are primarily composed of low-temperature brown serpentine (i.e., lizardite and/or chrysotile) and Fe-rich brucite assemblages forming mesh and bastite textures after olivine and orthopyroxene, respectively (fig. S1A). These rocks represent the first stage of forearc mantle wedge hydration at low temperatures (180° to 230°C, based on oxygen isotope estimates) and shallow depths (<13 km) (Fig. 1B) (36). The antigorite-bearing serpentinites (Atg/Liz- and Atg-serpentinites) are highly serpentinized (~100%) and display progressive replacement of lizardite-bearing textures (i.e., mesh and bastite) by antigorite, magnetite, and Fe-poor brucite ± andradite (fig. S1B). Oxygen isotope thermometry data suggest that the recrystallization of lizardite into antigorite occurred at temperatures ranging from 250° to 320°C (36), consistent with thermodynamic predictions (38) and natural temperature estimates (39) for antigorite crystallization in subduction zones (Fig. 1B). These observations provide evidence for the formation of these samples during the progressive burial and hydration of forearc mantle at depths ranging from 13 to 18 km (36).

The blue-serpentinites are considered to be former forearc peridotites or serpentinites (i.e., Liz- or Atg-serpentinites) recrystallized at low temperatures (<130°C) during the final stages of clast exhumation. These samples are highly serpentinized (serpentinization degree >80%) and are characterized by the crystallization of serpentine textures (i.e., mesh or bastite textures) with a strong blue color (fig. S1C). These blue-serpentine textures largely overprint any primary mantle mineral or former serpentine (i.e., lizardite and antigorite) textures. These textures can also be associated with framboidal sulfides (mainly pyrite), probably formed by tardive bacterially mediated sulfate reduction (36) and Fe hydroxides.

## Iron distribution in forearc ultramafic clasts

The Liz-serpentinite clasts display high $Fe^{3+}/\Sigma Fe$ [= $Fe^{3+}/(Fe^{2+} + Fe^{3+})$] ratios, ranging from 0.24 to 0.37 (table S1), relative to the primitive mantle ($Fe^{3+}/\Sigma Fe$ ~ 0.1) (40). In these clasts, opaque phases form large subautomorphic crystals that are interpreted as relict mantle spinels (fig. S1A). Magnetic susceptibility and Curie temperature measurements (figs. S2 and S3) reveal that magnetite is nearly absent, with modal amount estimates ranging from <0.1 to 0.7 weight % (wt %; table S2). In agreement with these observations, only rare and small (<50-µm width) magnetite-bearing veins were observed by x-ray computed tomography (CT) analyses and in thin section (fig. S1A). The near-total absence of magnetite despite high bulk rock $Fe^{3+}/\Sigma Fe$ suggests the preferential partitioning of $Fe^{3+}$ into serpentine minerals, in agreement with the inferred low temperatures

of serpentinization (<200°C) (15). The Atg/Liz- and Atg-serpentinites display even higher $Fe^{3+}/\Sigma Fe$ ratios (0.33 to 0.67; table S1) relative to Liz-serpentinites. In these clasts, the growth of antigorite is associated with the formation of a dense network of discrete and elongated particles, identified as magnetite, organized along a single preferred orientation (fig. S1B). The magnetic susceptibility and Curie temperature measurements (figs. S2 to S3) confirm that magnetite is the main magnetic phase crystallizing in Atg-serpentinites, with modal amount estimates ranging from 0.2 to 5.9 wt % (table S2). Magnetite crystallization is always associated with Fe-poor brucite and sometimes, to a lesser extent, andradite (36).

The blue-serpentinites are characterized by a large range in $Fe^{3+}/\Sigma Fe$ (0.33 to 0.77; table S1). Curiously, samples with high $Fe^{3+}/\Sigma Fe$ ratios display relatively low magnetite contents (<1 wt %; fig. S2 and table S2), and Curie temperature measurements reveal an irreversible behavior during successive heating and cooling (fig. S3). These features can be attributed to the presence of sulfate and/or Fe sulfides (e.g., greigite), which are weakly magnetic and are irreversibly converted to magnetite and/or pyrrhotite during Curie temperature measurements. This is supported by x-ray CT analyses that show the existence of isolated large pyramidal sulfides in those rocks (fig. S1C). The quasi-absence of magnetite in blue-serpentinites and the irreversible magnetic behavior observed during Curie temperature measurements suggest that serpentine is the main $Fe^{3+}$ carrier in those clasts. These observations are consistent with low temperatures of serpentinization favoring the incorporation of $Fe^{3+}$ in serpentine minerals (15) coupled with oxidizing conditions accompanying near-seafloor alteration (35).

## Carbon bulk rock isotope analyses

Seventeen representative samples were analyzed for bulk carbon elemental and isotope compositions. Carbon analyses were performed for both inorganic [TIC (total inorganic carbon)] and organic [TOC (total organic carbon)] carbon compounds in addition to total carbon (TC; table S3). The Mariana forearc serpentinite clasts have TC concentrations [(C<sub>TC</sub>) = 389 to 12,523 parts per million (ppm)] and isotope compositions (δ<sup>13</sup>C<sub>TC</sub> = -11.8 to -25.5‰) that overlap with those of abyssal and HP-serpentinites (fig. S4A). There are no systematic variations in (C<sub>TC</sub>) or δ<sup>13</sup>C<sub>TC</sub> from Liz- (408 to 637 ppm; -11.8 to -19.2‰) to Atg/Liz- (403 to 1960 ppm; -12.4 to -17.7‰), Atg- (389 to 1142 ppm; -14.0 to -21.3‰), or blue-serpentinites (652 to 12,523 ppm; -15.6 to -25.5‰), the latter group displaying the most extreme variations in both (C<sub>TC</sub>) and δ<sup>13</sup>C<sub>TC</sub>.

To put these results in context, the TC concentrations (C<sub>TC</sub>) of the Mariana forearc serpentinite clasts are significantly greater than that inferred for the depleted mantle [(C<sub>TC</sub>) = 50.3 ppm] (41). Moreover, the Mariana forearc serpentinites are characterized by similar total concentration of organic carbon (C<sub>TOC</sub>) (194 to 10,291 ppm) and δ<sup>13</sup>C<sub>TOC</sub> (-23.3 to -29.4‰) relative to those of abyssal (42, 43) or HP-serpentinites (44, 45) (fig. S4B). However, these samples also display highly variable concentration of TIC (C<sub>TIC</sub>) (67 to 2323 ppm) and δ<sup>13</sup>C<sub>TIC</sub> (-2.2 to -16.3‰), with a few samples having abnormally low δ<sup>13</sup>C<sub>TIC</sub> when compared with literature data for similar samples (fig. S4C).

## Carbon microimaging

The serpentinite clasts contain significant and highly variable amounts of carbon (bulk-rock concentrations of 400 ppm to 1.2 wt %; table S3). Although no carbon-bearing minerals (e.g., calcite) and graphite
are observed in thin sections, Fourier transform infrared (FTIR) microimaging shows the presence of large areas (several hundreds of micrometers wide) of carbonaceous matter displaying vibrational bands characteristic of organic compounds in the Atg- and blue-serpentinites (Figs. 2 and 3).

In the Atg-serpentinites, the FTIR spectra of the carbonaceous matter are characterized by intense methyl (CH₃—) and methylene (CH₂—) C—H stretching bands within the 3000 to 2800 cm⁻¹ and

1500 to 1200 cm⁻¹ ranges, corresponding to aliphatic chains and by a marked band characteristic of oxygen-bearing functional groups at ~1730 cm⁻¹ (C=O from aliphatic aldehyde) (Fig. 2B). The computed R<sub>CH2/CH3</sub> ratio does not exceed 1, giving a raw carbon number for the aliphatic chains below 8 to 10 when isomers are considered. The organic signal in the FTIR spectra is consistently associated with additional broad bands around 1640, 1460, and 1400 cm⁻¹, which correspond to antigorite. Scanning electron microscopy (SEM)

<figure><img src="image_3.png" /><figcaption>**Fig. 2. Characterization of carbonaceous matter wetting antigorite-bearing Mariana forearc serpentinites.** (A) FTIR map of the aliphatic CH₂-CH₃ C—H stretching band area between 2800 and 3000 cm⁻¹ showing the distribution of organic carbon encountered in an Atg-serpentinite. Areas depleted in organic compounds (blue) spatially correlate with the presence of Fe-poor brucite ±antigorite, whereas the organic signal (from green to red) is systematically associated with antigorite. The organic fraction appears chemically homogeneous throughout the analyzed areas (fig. S6). a.u., arbitrary unit. (B) Examples of raw FTIR spectra of brucite, antigorite, and carbonaceous matter associated with antigorite (in red), extracted from the map shown in (A). The organic-related spectrum is characterized by bands at 1: ~2960 cm⁻¹, CH₃ asymmetric C—H stretching, 2: ~2920 cm⁻¹, CH₂ asymmetric C—H stretching, 3: ~2870 cm⁻¹, CH₃ symmetric C—H stretching, 4: ~2850 cm⁻¹, CH₂ symmetric C—C—H stretching, 5: ~1730 cm⁻¹, C=O stretching of aliphatic aldehyde, 6: ~1470 cm⁻¹, CH₃ asymmetric C—H bending, CH₂ scissoring, and 7: ~1380 cm⁻¹, CH₃ symmetric C—H bending. Band assignment are from (71, 72). (C and D) SEM-SE (secondary electron) images of carbonaceous matter (CM) films embedding antigorite. (C) The discontinuous carbonaceous matter film covers the mineral surface outcropping at some places. (D) The film of carbonaceous matter completely embeds antigorite needles, locally connecting two crystals (white arrow with magnified view as inset). (E) SEM-BSE (backscatter electron) image with associated EDS element maps (C, Si, Ca, and Fe) of an Atg-serpentinite showing the relationship between CM and the HP paragenesis, made of antigorite, Fe-poor brucite (Brc), andradite (Adr), and magnetite (Mgt).</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Carbonaceous matter characterization in blue-serpentinites. (A) SEM-BSE image and associated EDS element maps for carbon and iron in a mesh core crystallizing in blue-serpentinites. The mesh core is composed of a crystal mush with a high porosity and containing brucite, lizardite, and Fe hydroxides and sulfides. The mush porosity is filled with carbonaceous matter. The red and black stars correspond to the localization of punctual EDS (shown in inset) and FTIR (B) analyses. The black inset locates (C). (B) Raw FTIR spectra of the carbonaceous matter (in red) associated with lizardite (in black). The organic-related spectrum is characterized by bands at 1 to 4: 3105, 3082, 3061, and 3028 cm⁻¹, aromatic C–H stretching, 5 and 6: ~2920 and ~2850 cm⁻¹, CH₂ asymmetric and symetric C–H stretching, 7 to 9: 1603, 1493, and 1450 cm⁻¹, C=C stretching ± CH₂ scissoring, 1493 cm⁻¹, and 10: 1370 cm⁻¹ aromatic C–N stretching. (C) Magnified SEM-SE image of carbonaceous matter filling the interstitial porosity of the crystal mush. Black arrows indicate micro–Fe hydroxides and sulfides associated with the crystal mush. (D) SEM-SE image of carbonaceous matter with a tubular texture.</figcaption></figure>

observations reveal that carbonaceous matter forms large patches up to 100 µm wide, either covering the surface of antigorite crystals (Fig. 2, C to E) or forming thin films embedding antigorite needles (Fig. 2D). These films have a gel-like appearance and occasionally connect bathed antigorite crystals (Fig. 2E). Carbonaceous matter was also observed at the interface between antigorite and brucite or andradite crystals (Fig. 2E), although this was rare.

In blue-serpentinites, the carbonaceous matter show a contrasted FTIR spectral signature compared to Atg-serpentinites, with methyl and methylene C–H stretching bands at 2850 and 2925 cm⁻¹, corresponding to aliphatic chains, associated with aromatic C–H stretching bands at 3025, 3060, and 3085 cm⁻¹ (Fig. 3B). Three additional marked bands at 1603, 1493, and 1450 cm⁻¹ can be assigned to aromatic C=C stretching (± CH₂ scissoring) and an additional band at 1370 cm⁻¹ can be attributed to C–N stretching. SEM observations show that the carbonaceous matter is located in a porous crystal mush of serpentine, brucite, andradite, and Fe-rich phases (including Fe sulfides and Fe hydroxides; Fig. 3A). The carbonaceous matter precipitates as gelatinous patches of about 50 µm in diameter, filling the textural porosity (Fig. 3, A to C). Magnified images reveal that the organic matter can also display a tubular filament-like texture

(Fig. 3D). Energy dispersive x-ray spectrometry (EDS) analyses show that this carbonaceous matter signal is associated with a marked nitrogen peak (Fig. 3A), suggesting that the solid organic compounds observed in the blue-serpentinites may potentially host nitrogen.

DISCUSSION

Deciphering the biologic versus abiotic origin of HP
carbonaceous matter

The crystallization of the blue-serpentinites is considered to have taken place during clast exhumation, in subsurface low-temperature environments (Fig. 1), as shown by δ¹⁸O temperature estimates (<180°C) (36) and the low magnetite modal amounts (figs. S2 and S3 and table S2) determined for these samples. These conditions are compatible with microbial activity (i.e., ≤121°C, the current established limit for life) (46) and are comparable to those estimated by Plümper et al. (32), who consider subduction zones as the largest terrestrial microbial habitat, using the example of the Mariana forearc system. The δ¹³C$_{TIC}$ compositions of the blue-serpentinites are highly heterogeneous, varying from –3.4 to –16.3‰. These values are abnormally “light” (depleted in heavy carbon isotopes) relative to the
$\delta^{13}C_{TIC}$ compositions of abyssal (from 6.2 to -9.9‰) or HP-serpentinites (from 1.4 to -8.4‰) (fig. S4C) (46). These $\delta^{13}C_{TIC}$ variations are correlated with a decrease in TIC (%TIC) in the bulk rock, with organic-rich samples displaying high $\delta^{13}C_{TIC}$ signatures and inorganic carbon-rich samples having low $\delta^{13}C_{TIC}$ signatures (Fig. 4A). This observation is unexpected as carbonates are known to preferentially uptake isotopically “heavy” carbon (enriched in heavy carbon isotopes), while organic carbon is expected to preferentially uptake isotopically light carbon (42, 43). Such low $\delta^{13}C_{TIC}$ values may therefore reflect the oxidation of isotopically light organic carbon into carbonate, such that the carbonate inherits an “organic carbon” isotope signature, e.g.

$$CH_4 + 2H_2O = CO_2 + 4H_2 \quad (1)$$

<figure><img src="image_5.png" /><figcaption>Fig. 4. Covariations of carbonate isotopic composition in serpentinite clasts with elemental tracers. (A) $\delta^{13}C_{TIC}$ versus carbonate abundance (%TIC). (B) $\delta^{13}C_{TIC}$ versus $\delta^{56}Fe$. (C) $\delta^{13}C_{TIC}$ versus As. (D) $\delta^{13}C_{TIC}$ versus (B). Plots (B to D) only show data for Atg-serpentinites. HP-serpentinite data are from (44, 45).</figcaption></figure>

This mechanism is commonly observed in passive margins where authigenic carbonates displaying abnormally low $\delta^{13}C_{TIC}$, down to -60‰, are formed by anaerobic microbial consortia through syntrophic methane (± other hydrocarbons) oxidation and sulfate reduction (47, 48)

$$CH_4 + SO_4^{2-} = HCO_3^- + HS^- + H_2O \quad (2)$$

The microbial oxidation of organic carbon compounds in blue-serpentinites is also supported by the predominance of authigenic aragonite and calcite in the top meters of the mud volcanoes (49-51). On the basis of Sr and C isotope analyses, their precipitation was interpreted as resulting from microbially mediated fluid interactions between shallow seawater and deep-sourced and slab-related fluids. Furthermore, in the blue-serpentinites, carbonaceous matter is found within serpentinite pore spaces, in crystal mushes composed of lizardite, brucite, andradite, Fe sulfide, and Fe hydroxides (Fig. 3, A to D). This close association of carbonaceous matter with sulfides, which can display framboidal textures (36), is compatible with Eq. 2, the later textures being commonly attributed to sulfate-reducing bacteria (52). Microbial activity is also supported by the tubular, i.e., filament-like, textures sometimes found in blue-serpentinites (Fig. 3D) and the spectral signature of the carbonaceous matter, which is composed of a mixture of nitrogen-bearing aliphatic and aromatic organic compounds (Fig. 3B), as described by Plümper et al. (32). Together, these geochemical and petrographic observations support the biologically mediated breakdown of organic compounds to form carbonate minerals during exhumation of the blue-serpentinites to near-seafloor depths.

Compared to the blue-serpentinites, the Atg-serpentinites display relatively homogeneous $\delta^{13}C_{TIC}$ values (-2.2 to -6.6‰) for a large range of %TIC (13.6 to 67.6%). No obvious correlation between $\delta^{13}C_{TIC}$ and %TIC is observed (Fig. 4A). The $\delta^{13}C_{TIC}$ values of these samples also fall within the range traditionally observed for HP-serpentinites (Fig. 4A and fig. S4C), as do the carbon isotopes values of TOC and TC. Cannao et al. (45) postulate that the devolatilization of carbon-rich sediments under oxidizing conditions during subduction is likely to produce fluids with relatively low $\delta^{13}C_{TIC}$ (<0‰) that could then react with peridotites and/or HP-serpentinites and affect low $\delta^{13}C_{TIC}$ signatures to these rocks. The similar $\delta^{13}C_{TIC}$ signatures of HP-serpentinites and the Atg-serpentinites suggest that fluids derived from the breakdown of carbon-rich slab sediments under oxidizing conditions may also be important in forearc peridotite serpentinization. Strong correlations between bulk rock concentrations of fluid mobile elements (FMEs) such as As and B, which are characteristic of sedimentary components in subduction zones (53-55), and $\delta^{13}C_{TIC}$ in the Atg-serpentinites (Fig. 4, C and D) further support the mobilization of carbonate derived from fluids released by sediment during subduction.

In the Atg-serpentinites, carbonaceous matter forms thin gel textures covering antigorite or patches filling the interstices between brucite, magnetite, and andradite. The close and systematic association of gel-like carbonaceous matter wetting antigorite suggests that these phases are syngenetic. In addition, the high crystallization temperatures of antigorite in the forearc (310° to 410°C; Fig. 1B) (36) relative to those considered compatible with life (≤121°C) (46), the textural association between carbonaceous matter and antigorite provides additional evidence for an abiotic origin of the carbonaceous matter. Consistent with this interpretation, the FTIR spectral
signature of the Atg-serpentinite carbonaceous matter differs to those observed in blue-serpentinites. In these samples, the FTIR spectral signature is characterized by a strong aliphatic organic signal, with carbon number <8 to 10, high content in O-bearing functional groups, and a low level of aromatic moieties (Fig. 2), suggestive of organic acids that are common by products of abiotic organic synthesis.

Overall, these observations suggest that biotic and abiotic processes lead to the formation of solid organic compounds with different petrological habits, as well as geochemical and FTIR spectral signatures in blue- and Atg-serpentinites. The Liz- and Atg/Liz-serpentinites display intermediate δ¹³CₗtIC (-4.7 to -14.6‰) and %TIC (24 to 83%) between blue- and Atg-serpentinites when plotted on Fig. 4A. Notably, these samples do not display any correlations between δ¹³CₗtIC and FME concentrations, suggesting that these samples and the blue-serpentinites do not record fluids released by sediment during subduction. This observation combined with low temperature of serpentinization (36) suggest that these samples were serpentinized at shallow depths by CO₂-free fluids.

## Linking redox reactions to HP abiotic synthesis and storage of solid carbon compounds

In the deep part of the forearc (i.e., > 10-km depth), the serpentinization of mantle peridotites through interaction with water-rich fluids derived from the slab is commonly associated with the oxidation of ferrous Fe leading to a progressive increase of Fe³⁺/ΣFe in the bulk rock and H₂ production in the fluid phase through water reduction

$$3 \text{FeO} + \text{H}_2 \text{O}_{\text{(aq)}} = \text{Fe}_3\text{O}_4 + \text{H}_{2\text{(aq)}} \quad (3)$$

In agreement with this, the Mariana serpentinized clasts display high Fe³⁺/ΣFe ratios (0.24 to 0.67) relative to mantle peridotites (~0.1) (40). However, for a given serpentinization degree, the Atg-serpentinites display higher Fe³⁺/ΣFe ratios than Liz- or Atg/Liz-serpentinites (fig. S5), which is consistent with the progressive recrystallization of serpentinite at high temperatures and an increase in the modal amount of magnetite in bulk rock (fig. S2). These results demonstrate the progressive oxidation of iron through the high-temperature recrystallization of lizardite and Fe-rich brucite into antigorite, Fe-poor brucite, and magnetite. Negative correlations have also been observed between Fe³⁺/ΣFe and serpentinite bulk rock stable iron isotopes (δ⁵⁶Fe) in the same sample suite (37). Debret et al. (37) interpreted these correlations and the high Fe³⁺/ΣFe ratios of Atg-serpentinites in terms of the transfer of slab-derived CO₂-bearing fluids to the forearc mantle wedge. In this scenario, these fluids preferentially complex isotopically light Fe and can modify the redox state of the forearc mantle at high temperature according to the simplified reaction

$$\text{FeCO}_{3\text{(aq)}} + 3\text{FeO} = 2\text{Fe}_2\text{O}_3 + \text{C} \quad (4)$$

A notable correlation between δ¹³CₗtIC values, which are believed to retain the isotopic signature of sediment-derived fluids, and δ⁵⁶Fe values is observed in the Atg-serpentinites (Fig. 4B). Such a correlation is not observed in blue-serpentinites, where the carbon isotope signature appears to have been reset by microbial activity associated with clast exhumation. Any late carbon recrystallization in an open system or redistribution in a closed system during clast exhumation should induce an isotope fractionation of δ¹³CₗtIC and would destroy any correlations between δ¹³CₗtIC and δ⁵⁶Fe, as for example,

seafloor-related serpentinization is not associated with Fe isotope fractionation (56). We therefore interpret the correlation between δ¹³CₗtIC and δ⁵⁶Fe in the Atg-serpentinites in terms of transfer of oxidizing fluids (e.g., CO₂) from slab to the forearc mantle at depth. Such a scenario is consistent with the large fluxes of subducting carbonate-rich sediments reported at ODP legs 129 and 189 nearby Mariana trench that carry an isotopically heavy signature in carbon (in average, 2.3‰) (57). Heavy carbon isotope signatures (-0.5‰) have also been reported by Hilton et al. (58) in geothermal fluids (from geothermal wells, fumaroles, and hot springs) collected from the Mariana Arc islands. These authors conclude that 87% of the emitted CO₂ in Mariana volcanic system is slab derived. Similarly, Resing et al. (59) analyzed He isotopes in hydrothermal fluids from active submarine volcanoes and proposed that >80% of the CO₂ venting from these volcanoes had a slab source. Further evidence in favor of a slab origin comes in the form of combined ¹⁴³Nd/¹⁴⁴Nd and ¹³⁸Ce/¹⁴²Ce isotope data for Mariana volcanic rocks, which reveals a strong subducted sedimentary component in the source of Mariana arc system (60). In contrast, the He and C isotope signatures of Mariana back arc basalts are indistinguishable from those of mid-ocean ridge basalt, suggesting a mantle source-like signature (between -6.1 and -5.5‰) and, hence, limited carbon transfer beyond subarc depth (61). Overall, these observations suggest that most of carbon transfer between the slab and the mantle wedge occurs under oxidized form (e.g., CO₂) during the first stages of subduction, placing the serpentinized forearc mantle in an ideal position to become a carbon sink in this system.

Despite recording the transfer of CO₂ in slab-derived fluids, the carbon budget of the Atg-serpentinites is, in fact, dominated by solid organic carbon species presumably formed by the reduction of the fluid-derived CO₂ component. The exact forms of organic carbon that can be formed through these redox processes (e.g., Eq. 4) operating in the forearc are numerous and depend on carbon sources, redox conditions associated with slab decarbonation, and local chemical and redox equilibria. For example, oxidized conditions may lead to a wide variety of organic compounds in addition to inorganic carbon species, as suggested by the O-bearing functional groups associated with the carbonaceous matter in Atg-serpentinite (Fig. 2B). The transfer of both H₂O (Eq. 3) and CO₃²⁻ (Eq. 4) to the forearc mantle may therefore also lead to the abiotic synthesis of organic acids, e.g.

$$2 \text{FeCO}_{3\text{(aq)}} + 6 \text{FeO} + 2 \text{H}_2 \text{O}_{\text{(aq)}} = 4 \text{Fe}_2 \text{O}_3 + \text{CH}_3 \text{COO}^- + \text{H}^+_{\text{(aq)}} \quad (5)$$

In agreement with such a scenario, high concentrations of dissolved H₂ (10 to 1000 µM), CH₄ (1 to 10,000 µM), and C₂H₆ (0.1 to 1000 µM), as well as formate (HCOO⁻, ~100 µM) and acetate (CH₃COO⁻, ~40 µM) organic acids were measured in the interstitial waters most distant to the trench mud volcanoes (i.e., Asùt Tesoru, South Chamorro, and Conical), which are considered to sample the deepest fluid/rock interactions between the slab and the forearc mantle wedge (35, 62, 63). In contrast, comparatively low levels (H₂ = 1 to 100 µM; CH₄ = 0.1 to 10 µM; C₂H₆ = not determined, HCOO⁻ ~2 µM; CH₃COO⁻ ~ 2 µM) were detected at the mud volcanoes located close to the trench (i.e., Yinazao) reflecting the limitation of carbonate reduction processes in the shallow parts of the forearc (35, 62). The link between serpentinization reactions and organic carbon production with depth is also supported by the
presence of solid organic compounds associated with antigorite and that have a strong aliphatic signature (Fig. 2). It shows that the rate of iron oxidation and, consequently, the H₂ production increases with depth where it enhances the reduction of CO₂-bearing slab fluids. This results in the formation of organic carbon compounds in the form of organic carbon-rich fluids discharged through mud volcanoes or carbonaceous matter with an aliphatic signature trapped within mantle rocks. Given that the solubility of alkanes decreases with carbon number (64) and that organic acids >C5 are insoluble in water owing to hydrophobicity of long aliphatic chains, we hypothesize that both the volatile hydrocarbons and organic acids detected at mud volcanoes and the clast-hosted carbonaceous matter may derive from the same redox reactions, with the fraction soluble in water, feeding hydrothermal activity toward the surface and the insoluble fraction being stored in HP ultramafic clasts.

Despite the fact that abundant carbonate minerals may be sampled in the shallow Mariana forearc, as late calcite and aragonite veins, chimneys (49–51), or stored within blue-serpentinite (this study), their carbon isotope composition suggest that they formed nearby the surface through microbially mediated fluid interaction between shallow seawater and deep-sourced fluids. It is therefore not clear how much of this carbon is derived from the subduction process and how much is derived from late seawater circulation. Hence, these authigenic carbonates provide few constraints on deep carbon transfers between the slab and the mantle wedge. In contrast, the Atg-serpentinites preserve a deep geochemical signature reflecting CO₂ transfer and reduction during the serpentinization of the forearc (Fig. 4). In these samples, the solid organic carbon fraction dominates over the inorganic carbon fraction, the former representing up to 80% of stored carbon.

Using our data in combination with recent estimates of total water flux in the Mariana subduction zone (65), we can provide a semiquantitative estimate of the maximum amount of solid organic carbon stored during forearc serpentinization (table S4). This calculation considers that between 20 and 74% of total water injected at the trench is transferred to the forearc mantle wedge by slab-derived fluids (66, 67), where it is stored in serpentine minerals during peridotite serpentinization. The amount of carbon stored as solid organic compounds is then considered constant and constrained by the molar abundance of C<sub>TOC</sub> relative to H₂O in Atg-serpentinites (table S4). Using this calculation, we can propose that between 0.05 and 0.39 megatons year⁻¹ of solid organic carbon is currently formed and stored in the Mariana forearc area. This represents between 5 and 54% of the carbon input of the Mariana subduction zone according to recent estimates of subducting sediments input by Clift (68) and Sadofsky and Bebout (57). Our estimate of the organic carbon input serves to underline the potential of the serpentinized forearc mantle as a major reservoir of solid organic carbon in subduction zones. The fate of this solid organic carbon reservoir needs further investigation. It may be transferred to the surface during clast and mud exhumation at the seafloor where it could serve as a carbon source for microbial ecosystems (32, 36). If so, the low δ¹³C<sub>TIC</sub> (<< 10‰; Fig. 4A) of the blue-serpentinites may reflect the microbial conversion of deep carbonaceous matter to authigenic carbonate. On the other hand, solid organic carbon may accumulate in the forearc along subduction until its transfer is operated to greater depths by mantle convection (69) or to the surface by plate tectonic (e.g., orogeny). The stability and devolatilization of such a solid organic carbon reservoir will then be intimately related

to sharp changes of fO₂ and pH over a large range of pressure and temperature (up to 6 GPa and 1000°C) (21, 31).

## MATERIALS AND METHODS

The IODP identification numbers, bulk and in situ major (including bulk Fe³⁺/ΣFe quantification), and trace element chemistry of the studied samples are detailed in (36, 37).

### Magnetic properties

Density and magnetic susceptibility measurements were performed on 21 clasts from Yinazao, 37 clasts from Fantangisña, and 31 clasts from Asùt Tesoru mud volcanoes (table S2). Analyzed sample volume was between 5 and 15 cm³. Measurements were performed either offshore during the IODP Expedition 366 (35) or onshore at the University of Lorraine (Nancy, France), using the same analytical devices and protocols. Density measurements were carried out by constant-volume gas pycnometry, using a Micromeritics AccuPyc II 1340 helium pycnometer at room temperature with a testing chamber of 3.5 cm³. Bulk magnetic susceptibility was measured at room temperature with an MFK1-A Kappabridge instrument (constant magnetizing field and operating frequency) based on the precise sample volume measurements from gas pycnometry analyses. A subset of 15 clasts with different mineralogy (i.e., blue-, Liz-, and Atg-serpentinites) was then selected for heat- and cool-induced magnetic transition measurements using Kappabridge devices CS-L Cryostat and CS-4 Furnace at the CRPG (Nancy, France). Measurements were performed under argon atmosphere at three temperature intervals: from -194° to 20°C using liquid nitrogen, from 20° to 700°C (heating curve), and from 700° to 20°C (cooling curve).

### X-ray CT

X-ray CT scanning of four samples with different mineralogy (i.e., blue-, Liz-, and Atg-serpentinites) was performed using a Phoenix Nanotom X-ray CT scanner at the University of Lorraine, GeoRessources. Scans were performed at 90 keV and 125 mA, with a voxel size of 5 to 6 µm³. Scans were reconstructed with Avizo software using a filtered rear projection method and a beam hardening correction. Reconstructed volumes were saved as tiff 16-bit uncompressed images.

### Sample preparation for carbon imaging

Sample preparation for organic carbon characterization was performed at the Institut de physique du globe de Paris (IPGP, France). Drilled rock samples were sawed with sterile ultrapure water to extract the inner core, free of possible postsampling contamination. The saw was previously rinsed twice with sterile ultrapure water. The inner core was then manipulated using clean pliers, thinned, and polished on both faces (down to a thickness of 10 µm) with pure ethanol using silicon carbide polishing disks without any use of resin or glue.

### FTIR microimaging

FTIR microspectroscopy was performed at IPGP on a Thermo Fisher Scientific iN10 MX microscope (Ever-Glo conventional infrared source) equipped with a 15× objective (numerical aperture = 0.7) and a liquid nitrogen-cooled MCT-A detector. The incident beam was collimated to a sample area of 20 µm by 20 µm. The thinned rock sample was deposited on a BaF₂ window without any
treatment. FTIR hyperspectral maps were acquired in transmission in the range 4000 to 675 cm⁻¹ using step size of 20 µm, and 64 accumulations per spectrum per pixel and the spectral resolution was set to 8 cm⁻¹. Spectrum analyses and reconstructions based on band areas were performed using OMNIC software (Thermo Fisher Scientific). A second derivative was applied to raw data to remove baseline variations using the Stavitsky-Golay method (polynomial order, 2; filter width, 15). Second derivative maps were obtained from unit vector-normalized maps using MATLAB and PLS_Toolbox software (Eigenvector Research Inc.). Curve fitting in the absorbance region of the aliphatic moieties (i.e., 2800 to 3000 cm⁻¹) was computed with Wire 3.2 software (Renishaw). Five mixed Lorentzian/Gaussian functions were used to account for all the symmetric and asymmetric C–H stretching bands. The methyl-to-methylene ratio (R<sub>CH2/CH3</sub>) was computed as the area ratio of the fitted asymmetric C–H stretching bands of CH<sub>2</sub> to CH<sub>3</sub> at ~2930/2960 cm⁻¹, respectively. R<sub>CH2/CH3</sub> allows rough assessment of the length and degree of branching of side groups associated with aliphatic chains (70).

## Scanning electron microscopy

SEM observations were performed at IPGP using a Zeiss Auriga FEG-FIB field emission scanning electron microscope. Samples were Au-coated. Images were collected using backscattered electron and secondary electron detector at high and low current, respectively, with accelerating voltage ranging from 10 to 15 kV. EDS measurements were performed at 15-kV accelerating voltage using a Bruker detector.

## Carbon bulk rock analyses

The TC concentrations (C<sub>TC</sub>) and isotope ratios (δ<sup>13</sup>C<sub>TC</sub>) analyses were performed using a Thermo Fisher Scientific EA IsoLink IRMS System at the CRPG (Nancy, France). Two to three aliquots of 30 mg of each powdered samples were wrapped in tin capsules and then burned at 1020°C by flash combustion in a quartz reactor filled with chromium oxide, pure copper, and silvered cobalt oxide. Produced gases (N<sub>2</sub> and CO<sub>2</sub>) were separated on a chromatographic column maintained at 70°C, and TC concentrations and isotope ratios were then measured using a Thermo Fisher Scientific Delta V Advantage continuous flow isotope ratio mass spectrometer (IRMS). TC isotopic composition was determined by comparison with CRPG and international standards that were analyzed along with samples: (i) δ<sup>13</sup>C<sub>TC</sub> (BFSd) = -21.48‰; (ii) δ<sup>13</sup>C<sub>TC</sub> (CRPG_M2) = -24.97‰; and (iii) δ<sup>13</sup>C<sub>TC</sub> (PSd) = -22.58‰. Values are quoted in the delta notation in ‰ relative to the V-PDB (Vienna PeeDee Belemnite) standard, and the reproducibility is better than 0.2‰. Internal standards were used to determine carbon concentrations: (i) [C<sub>TC</sub> (BFSd)] = 0.52 wt %; (ii) [C<sub>TC</sub> (CRPG_M2)] = 0.41 wt %. For each sample, the 2σ error is expected to be lower than 5% for (C<sub>TC</sub>) and 1‰ for δ<sup>13</sup>C<sub>TC</sub>. TOC concentrations (C<sub>TOC</sub>) and isotopic compositions (δ<sup>13</sup>C<sub>TOC</sub>) were determined for each whole-rock sample using the same instruments and calibration described for analyses of the TC. The inorganic carbon component was removed, before each analysis, by HCl fumigation during 5 days at 65°C.

TIC isotope ratios (δ<sup>13</sup>C<sub>TIC</sub>) were determined by using an auto sampler Gasbench coupled to a Thermo Fisher Scientific MAT 253 continuous flow IRMS at the CRPG (Nancy, France). For each sample, an aliquot of 1 to 100 mg of powder was reacted with 2 ml of supersaturated orthophosphoric acid at 70°C for 10 hours under a He atmosphere. Values are quoted in the delta notation in ‰ relative

to V-PDB. All sample measurements were adjusted to an internal reference (MCt: pure calcite; δ<sup>13</sup>C<sub>TIC</sub> = -8.64‰) calibrated on international standards (i) δ<sup>13</sup>C<sub>TIC</sub> (BR516) = -1.1‰ and (ii) δ<sup>13</sup>C<sub>TIC</sub> (BR8107) = -0.9‰. The reproducibility of the standards was better than 0.2‰. TIC contents (C<sub>TIC</sub>) of the samples were determined by comparison with two internal standards consisting in fine-grained marine sediments from the Bay of Bengal and routinely included during the analysis: (i) [CaCO<sub>3</sub> (BR516)] = 3.54 wt % and (ii) [CaCO<sub>3</sub> (BR8107)] = 6.16 wt %. The 2σ errors for carbonate contents are estimated to be lower than 5%.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at https://science.org/doi/10.1126/sciadv.abo2397

## REFERENCES AND NOTES

1. J. F. Kasting, D. H. Eggle, S. P. Raeburn, Mantle redox evolution and the oxidation state of the archean atmosphere. *J. Geol.* **101**, 245–257 (1993).

2. F. Gaillard, B. Scaillet, N. T. Arndt, Atmospheric oxygenation caused by a change in volcanic degassing pressure. *Nature* **478**, 229–232 (2011).

3. M. S. Duncan, R. Dasgupta, Rise of Earth's atmospheric oxygen controlled by efficient subduction of organic carbon. *Nat. Geosci.* **10**, 387–392 (2017).

4. V. Stagno, D. O. Ojwang, C. A. McCammon, D. J. Frost, The oxidation state of the mantle and the extraction of carbon from Earth's interior. *Nature* **493**, 84–88 (2013).

5. S. Poli, Carbon mobilized at shallow depths in subduction zones by carbonatitic liquids. *Nat. Geosci.* **8**, 633–636 (2015).

6. P. J. Gorman, D. M. Kerrick, J. A. D. Connolly, Modeling open system metamorphic decarbonation of subducting slabs. *Geochem. Geophys. Geosyst.* **7**, Q04007 (2006).

7. R. Dasgupta, M. M. Hirschmann, The deep carbon cycle and melting in Earth's interior. *Earth Planet. Sci. Lett.* **298**, 1–13 (2010).

8. P. B. Kelemen, C. E. Manning, Reevaluating carbon fluxes in subduction zones, what goes down, mostly comes up. *Proc. Natl. Acad. Sci. U.S.A.* **112**, E3997–E4006 (2015).

9. J. C. Alt, D. A. H. Teagle, The uptake of carbon during alteration of ocean crust. *Geochem. Cosmochim. Acta* **63**, 1527–1535 (1999).

10. T. Plank, The chemical composition of subducting sediments, in *Treatise on Geochemistry: Second Edition*, R. F. Keeling, Ed. (Elsevier, 2014), vol. 2, pp. 607–629.

11. J. C. Alt, E. M. Schwarzenbach, G. L. Früh-green, W. C. Shanks III, S. M. Bernasconi, C. J. Garrido, L. Crispini, L. Gaggero, J. A. Padrón-navarta, C. Marchesi, The role of serpentinites in cycling of carbon and sulfur: Seafloor serpentinization and subduction metamorphism. *Lithos* **178**, 40–54 (2013).

12. S. Shilobreeva, I. Martinez, V. Busigny, P. Agrinier, C. Laverne, Insights into C and H storage in the altered oceanic crust: Results from ODP/IODP Hole 1256D. *Geochim. Cosmochim. Acta* **75**, 2237–2255 (2011).

13. D. V. Bekaert, S. J. Turner, M. W. Broadley, J. D. Barnes, S. A. Halldórsson, J. Labidi, J. Wade, K. J. Walowski, P. H. Barry, Subduction-driven volatile recycling: A global mass balance. *Annu. Rev. Earth Planet. Sci.* **49**, 37–70 (2021).

14. T. M. McColлом, J. S. Seewald, Abiotic synthesis of organic compounds in deep-sea hydrothermal environments. *Chem. Rev.* **107**, 382–401 (2007).

15. F. Klein, W. Bach, S. E. Humphris, W. A. Kahl, N. Jöns, B. Moskowitz, T. S. Berquó, Magnetite in seafloor serpentinite-Some like it hot. *Geology* **42**, 135–138 (2014).

16. M. Andreani, M. Muñoz, C. Marcaillou, A. Delacour, μXANES study of iron redox state in serpentine during oceanic serpentinization. *Lithos* **178**, 70–83 (2013).

17. N. G. Holm, J. L. Charlou, Initial indications of abiotic formation of hydrocarbons in the Rainbow ultramafic hydrothermal system, Mid-Atlantic Ridge. *Earth Planet. Sci. Lett.* **191**, 1–8 (2001).

18. B. Ménez, C. Pisapia, M. Andreani, F. Jamme, Q. P. Vanbellingen, A. Brunelle, L. Richard, P. Dumas, M. Réfrégiers, Abiotic synthesis of amino acids in the recesses of the oceanic lithosphere. *Nature* **564**, 59–63 (2018).

19. M. Andreani, B. Ménez, New perspectives on abiotic organic synthesis and processing during hydrothermal alteration of the oceanic lithosphere, in *Deep Carbon: Past to Present*, B. N. Orcutt, I. Daniel, R. Dasgupta, Eds. (Cambridge Univ. Press, 2019), pp. 447–479.

20. M. C. Sforna, D. Brunelli, C. Pisapia, V. Pasini, D. Malferrari, B. Ménez, Abiotic formation of condensed carbonaceous matter in the hydrating oceanic crust. *Nat. Commun.* **9**, 5049 (2018).

21. D. A. Sverjensky, V. Stagno, F. Huang, Important role for organic carbon in subduction-zone fluids in the deep carbon cycle. *Nat. Geosci.* **7**, 909–913 (2014).

22. R. Tao, L. Zhang, M. Tian, J. Zhu, X. Liu, J. Liu, H. E. Höfer, V. Stagno, Y. Fei, Formation of abiotic hydrocarbon from reduction of carbonate in subduction zones: Constraints
from petrological observation and experimental simulation. *Geochem. Cosmochim. Acta* **239**, 390–408 (2018).

23. F. Huang, I. Daniel, H. Cardon, G. Montagnac, D. A. Sverjensky, Immiscible hydrocarbon fluids in the deep carbon cycle. *Nat. Commun.* **8**, 15798 (2017).

24. A. Vitale Brovarone, I. Martinez, A. Elmaleh, R. Compagnoni, C. Chaduteau, C. Ferraris, I. Esteve, Massive production of abiotic methane during subduction evidenced in metamorphosed ophiarcarbonates from the Italian Alps. *Nat. Commun.* **8**, 14134 (2017).

25. A. Vitale-Brovarone, D. A. Sverjensky, F. Piccoli, F. Ressico, D. Giovannelli, I. Daniel, Subduction hides high-pressure sources of energy that may feed the deep subsurface biosphere. *Nat. Commun.* **11**, 3880 (2020).

26. A. Boutier, A. Vitale Brovarone, I. Martinez, O. Sissmann, S. Mana, High-pressure serpentinization and abiotic methane formation in metaperidotite from the Appalachian subduction, northern Vermont. *Lithos* **396–397**, 106190 (2021).

27. J. Zhu, L. Zhang, R. Tao, Y. Fei, The formation of graphite-rich eclogite vein in S.W. Tianshan (China) and its implication for deep carbon cycling in subduction zone. *Chem. Geol.* **533**, 119430 (2020).

28. M. E. Galvez, O. Beyssac, I. Martinez, K. Benzerara, C. Chaduteau, B. Malvoisin, J. Malavieille, Graphite formation by carbonate reduction during subduction. *Nat. Geosci.* **6**, 473–477 (2013).

29. P. Bouilhol, B. Debret, E. C. Inglis, M. Warembourg, T. Grocolas, T. Rigaudier, J. Villeneuve, K. W. Burton, Decoupling of inorganic and organic carbon during slab mantle devolatilisation. *Nat. Commun.* **13**, 308 (2022).

30. M. Zolotov, E. Shock, Abiotic synthesis of polycyclic aromatic hydrocarbons on Mars. *J. Geophys. Res. E Planets* **104**, 14033–14049 (1999).

31. A. D. Chanyshev, K. D. Litasov, A. F. Shatskiy, I. S. Sharygin, Y. Higo, E. Ohtani, Transition from melting to carbonization of naphthalene, anthracene, pyrene and coronene at high pressure. *Phys. Earth Planet. Inter.* **270**, 29–39 (2017).

32. O. Plümper, H. E. King, T. Geisler, Y. Liu, S. Pabst, I. P. Savov, D. Rost, T. Zack, Subduction zone forearc serpentinites as incubators for deep microbial life. *Proc. Natl. Acad. Sci. U.S.A.* **114**, 4324–4329 (2017).

33. P. Fryer, J. P. Lockwood, N. Becker, S. Phipps, C. S. Todd, Significance of serpentine mud volcanism in convergent margins, in *Ophiolites and Oceanic Crust; New Insights from Field Studies and the Ocean Drilling Program*, Y. Dilek, E. M. Moores, D. Elthon, A. Nicolas, Eds. (Geological Society of America, 2000), pp. 35–51.

34. P. Fryer, Serpentinite mud volcanism: Observations, processes, and implications. *Ann. Rev. Mar. Sci.* **4**, 345–373 (2012).

35. P. Fryer, C. G. Wheat, T. Williams; the Expedition 366 Scientists, Mariana convergent margin and south chamorro seamount. *Proc. Int. Ocean Discov. Progr.* **366**, 1–12 (2018).

36. B. Debret, E. Albers, B. Walter, R. Price, J. D. Barnes, H. Beunon, S. Faqc, D. P. Gillikin, N. Mattielli, H. Williams, Shallow forearc mantle dynamics and geochemistry: New insights from the IODP expedition 366. *Lithos* **326–327**, 230–245 (2019).

37. B. Debret, C. D. J. Reekie, N. Mattielli, H. Beunon, B. Ménez, I. P. Savov, H. M. Williams, Redox transfer at subduction zones: Insights from Fe isotopes in the Mariana forearc. *Geochemical Perspective Lett.* **12**, 46–51 (2020).

38. B. W. Evans, The serpentinite multisystem revisited: Chrysotile is metastable. *Int. Rev. Geol.* **46**, 479–506 (2004).

39. S. Schwartz, S. Guillot, B. Reynard, R. Lafay, B. Debret, C. Nicollet, P. Lanari, A. L. Auzende, Pressure-temperature estimates of the lizardite/antigorite transition in high pressure serpentinites. *Lithos* **178**, 197–210 (2013).

40. D. Canil, H. S. C. O'Neill, D. G. Pearson, R. L. Rudnick, W. F. McDonough, D. A. Carswell, Ferric iron in peridotites and mantle oxidation states. *Earth Planet. Sci. Lett.* **123**, 205–220 (1994).

41. V. J. M. Salters, A. Stracke, Composition of the depleted mantle. *Geochem. Geophys. Geosyst.* **5**, Q05B07 (2004).

42. A. Delacour, G. L. Früh-Green, S. M. Bernasconi, P. Schaeffer, D. S. Kelley, Carbon geochemistry of serpentinites in the Lost City Hydrothermal System (30°N, MAR). *Geochem. Cosmochim. Acta* **72**, 3681–3702 (2008).

43. E. M. Schwarzenbach, G. L. Früh-Green, S. M. Bernasconi, J. C. Alt, A. Plas, Serpentinization and carbon sequestration: A study of two ancient peridotite-hosted hydrothermal systems. *Chem. Geol.* **351**, 115–133 (2013).

44. E. M. Schwarzenbach, M. J. Caddick, M. Petroff, B. C. Gill, E. H. G. Cooperdock, J. D. Barnes, Sulphur and carbon cycling in the subduction zone mélangé. *Sci. Rep.*, 15517 (2018).

45. E. Cannàò, M. Tiepolo, G. E. Debout, M. Scambelluri, Into the deep and beyond: Carbon and nitrogen subduction recycling in secondary peridotites. *Earth Planet. Sci. Lett.* **543**, 116328 (2020).

46. K. Kashefi, D. R. Lovley, Extending the upper temperature limit for life. *Science* **301**, 934 (2003).

47. A. Y. Lein, Authigenic carbonate formation in the ocean. *Lithol. Miner. Resour.* **39**, 1–30 (2004).

48. D. L. Valentine, W. S. Reeburgh, New perspectives on anaerobic methane oxidation. *Envir. Microbiol.* **2**, 477–484 (2000).

49. E. Frery, P. Fryer, W. Kurz, A. Nguyen, O. Sissmann, T. Uysal, J. Zhao, Episodicity of structural flow in an active subduction system, new insights from mud volcano's carbonate veins - Scientific Ocean drilling expedition IODP 366. *Mar. Geol.* **434**, 106431 (2021).

50. M. Freitas, V. H. Magalhães, M. R. Azevedo, L. Pinheiro, E. Salgueiro, F. F. Abrantes, Authigenic carbonate precipitation at the Yinazao serpentinite mud volcano. *Comun. Geol.* **107**, 61–64 (2019).

51. J. A. Haggerty, Petrology and geochemistry of neogene sedimentary rocks from Mariana Forearc seamounts: Implications for emplacement of the seamounts, in *Seamounts Islands and Atollas, Geophysical Monograph Series* **43**, B. Keating, P. Fryer, R. Batiza, G. Boehlert, Eds. (AGU, 1987), pp. 175–185.

52. R. Popa, B. K. Kinkle, A. Badescu, Pyrite framboids as biomarkers for iron-sulfur systems. *Geomicrobiol. J.* **21**, 193–206 (2004).

53. G. E. Debout, M. L. Fogel, P. Cartigny, Nitrogen: Highly volatile yet surprisingly compatible. *Elements* **9**, 333–338 (2013).

54. B. Debret, P. Bouilhol, M. L. Pons, H. Williams, Carbonate transfer during the onset of slab devolatilization: New insights from Fe and Zn stable isotopes. *J. Petrol.* **59**, 1145–1166 (2018).

55. I. P. Savov, J. G. Ryan, M. D'Antonio, P. Fryer, Shallow slab fluid release across and along the Mariana arc-basin system: Insights from geochemistry of serpentinized peridotites from the Mariana fore arc. *J. Geophys. Res. Solid Earth* **112**, B09205 (2007).

56. B. Debret, H. Beunon, N. Mattielli, M. Andreani, I. Ribeiro da Costa, J. Escartín, Ore component mobility, transport and mineralization at mid-oceanic ridges: A stable isotopes (Zn, Cu and Fe) study of the Rainbow massif (Mid-Atlantic Ridge 36°14′N). *Earth Planet. Sci. Lett.* **503**, 170–180 (2018).

57. S. J. Sadofsky, G. E. Debout, Nitrogen geochemistry of subducting sediments: New results from the Izu-Bonin-Mariana margin and insights regarding global nitrogen subduction. *Geochem. Geophys. Geosyst.* **5**, Q0315 (2004).

58. D. R. Hilton, T. P. Fischer, E. Hauri, A. M. Shaw, Controls on the He–C systematics of the Izu-Bonin-Marianas (IBM) subduction zone. *Geochem. Cosmochim. Acta* **70**, A252 (2006).

59. J. A. Resing, E. T. Baker, J. E. Lupton, S. L. Walker, D. A. Butterfield, G. J. Massoth, K.-I. Nakamura, Chemistry of hydrothermal plumes above submarine volcanoes of the mariana arc. *Geochem. Geophys. Geosyst.* **10**, Q02009 (2009).

60. N. Bellot, M. Boyet, R. Doucelance, P. Bonnand, I. P. Savov, T. Plank, T. Elliott, Origin of negative cerium anomalies in subduction-related volcanic samples: Constraints from Ce and Nd isotopes. *Chem. Geol.* **500**, 46–63 (2018).

61. C. G. Macpherson, D. R. Hilton, K. Hammerschmidt, No slab-derived CO₂ in Mariana Trough back-arc basals: Implications for carbon subduction and for temporary storage of CO₂ beneath slow spreading ridges. *Geochem. Geophys. Geosyst.* **11**, Q11007 (2010).

62. C. D. Menzies, R. E. Price, J. Ryan, O. Sissmann, K. Takai, C. G. Wheat, Spatial variation of subduction zone fluids during progressive subduction: Insights from Serpentinite mud volcanoes. *Geochem. Cosmochim. Acta* **319**, 118–134 (2022).

63. P. Eickenbusch, K. Takai, O. Sissman, S. Suzuki, C. Menzies, S. Sakai, P. Sansjofre, E. Tasumi, S. M. Bernasconi, C. Glombitza, B. B. Jørgensen, Y. Morono, M. A. Lever, Origin of short-chain organic acids in serpentinite mud volcanoes of the Mariana convergent margin. *Front. Microbiol.* **10**, 1–21 (2019).

64. C. Marche, C. Ferronato, J. Jose, Solubilities of *n*-Alkanes (C₆ to C₈) in water at 30 °C to 180 °C. *J. Chem. Eng. Data* **48**, 967–971 (2003).

65. C. Cai, D. A. Wiens, W. Shen, M. Eimer, Water input into the Mariana subduction zone estimated from ocean-bottom seismic data. *Nature* **563**, 389–392 (2018).

66. V. Magni, P. Bouilhol, J. van Hunen, Deep water recycling through time. *Geochem. Geophys. Geosyst.* **15**, 4203–4216 (2014).

67. J. M. Ribeiro, C. T. A. Lee, An imbalance in the deep water cycle at subduction zones: The potential importance of the fore-arc mantle. *Earth Planet. Sci. Lett.* **479**, 298–309 (2017).

68. P. D. Clift, A revised budget for Cenozoic sedimentary carbon subduction. *Rev. Geophys.* **55**, 97–125 (2017).

69. T. Nagaya, A. M. Walker, J. Wooley, S. R. Wallis, K. Ishii, J.-M. Kendall, Seismic evidence for flow in the hydrated mantle wedge of the Ryukyu subduction zone. *Sci. Rep.* **6**, 29981 (2016).

70. R. Lin, G. P. Ritz, Studying individual macerals using i.r. microspectrometry, and implications on oil versus gas/condensate prereness and “low-rank” generation. *Org. Geochem.* **20**, 695–706 (1993).

71. V. Pasini, D. Brunelli, P. Dumas, C. Sandt, J. Frederick, K. Benzerara, S. Bernard, B. Ménez, Low temperature hydrothermal oil and associated biological precursors in serpentinites from mid-ocean ridge. *Lithos* **178**, 84–95 (2013).

72. B. H. Stuart, *Infrared Spectroscopy: Fundamentals and Applications* (John Wiley & Sons, 2004).

73. D. Bonnemains, J. Carlut, J. Escartín, C. Mével, M. Andreani, B. Debret, Magnetic signatures of serpentinization at ophiolite complexes. *Geochem. Geophys. Geosyst.* **17**, 2969–2986 (2016).
74. L. B. Stokking, D. L. Merrill, R. B. Haston, J. R. Ali, K. L. Saboda, Rock magnetic studies of serpentinite seamounts in the Mariana and Izu-Bonin regions, in *Proceedings of the Ocean Drilling Program, 125 Scientific Results*, P. Fryer, J. A. Pearce, L. B. Stokking, Eds. (College Station, TX, 1992), pp. 561–579; doi:10.2973/odp.proc.sr.125.158.1992.

75. Shipboard Scientific Party, SITE 1200, in *Proceedings of the Ocean Drilling Program, Initial Reports*, M. H. Salisbury, M. Shinohara, C. Richter, Eds. (College Station, TX, 2002), vol. 195, pp. 1–173; http://www-odp.tamu.edu/publications/195_IR/chap_03/chap_03.htm.

76. A. P. Roberts, L. Chang, C. J. Rowan, C.-S. Horng, F. Florindo, Magnetic properties of sedimentary greigite (Fe3S4): An update. *Rev. Geophys.* **49**, RG1002 (2011).

multidisciplinary program PARI and by Région Ile-de-France SESAME grant no. 12015908.
**Author contributions:** B.D. and B.W. selected the samples during the IODP Expedition 366 and did the petrological characterization, the magnetic measurements, and the CT scanning. B.D., B.M., H.B., and C.P. did the sample preparation and the carbon characterization by FTIR and SEM. T.R. and P.B. did the carbon isotope measurements. All the authors contributed to the project conception and the manuscript writing. **Competing interests:** The authors declare that they have no competing interests. **Data and materials availability:** All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.

**Acknowledgments:** This research used samples and data provided by the International Ocean Discovery Program (IODP). **Funding:** Financial support was provided by LabEx UnivEarthS Frontier MYSTHIC (ANR-10-LABX-0023 and ANR-18-IDEX-0001), the INSU Syster program, a NERC Moratorium Award (no. NE/P020860/1), the F.N.R.S., the Fondation Wiener Anspach (RG97553), and the Deep Carbon Observatory. Part of this work was supported by the IPGP

Submitted 24 January 2022

Accepted 19 July 2022

Published 16 September 2022

10.1126/sciadv.abo2397
High-pressure synthesis and storage of solid organic compounds in active
subduction zones

Baptiste Debret, Bndicte Mnez, Bastien Walter, Hlne Bouquerel, Pierre Bouilhol, Nadine Mattielli, Cline Pisapia, Thomas Rigaudier, and Helen Myfanwy Williams

Sci. Adv., 8 (37), eabo2397.

DOI: 10.1126/sciadv.abo2397

View the article online

https://www.science.org/doi/10.1126/sciadv.abo2397

Permissions

https://www.science.org/help/reprints-and-permissions