# IDENTIFICATION IN PORE WATERS OF RECYCLED SEDIMENT ORGANIC MATTER USING THE DUAL ISOTOPIC COMPOSITION OF CARBON ($\delta$ 13 C AND $\Delta$ 14 C): NEW DATA FROM THE CONTINENTAL SHELF INFLUENCED BY THE RHÔNE RIVER

J-P Dumoulin, C Rabouille, S Pourtout, B Bombled, B Lansard, I Caffy, S

Hain, M Perron, M Sieudat, B Thellier, et al.

► To cite this version:

J-P Dumoulin, C Rabouille, S Pourtout, B Bombled, B Lansard, et al. IDENTIFICATION IN PORE WATERS OF RECYCLED SEDIMENT ORGANIC MATTER USING THE DUAL ISOTOPIC COMPOSITION OF CARBON ($\delta$ 13 C AND $\Delta$ 14 C): NEW DATA FROM THE CONTINENTAL SHELF INFLUENCED BY THE RHÔNE RIVER. Radiocarbon, 2022, 64 (6), pp.1617-1627. 10.1017/RDC.2022.71. hal-03971089

HAL Id: hal-03971089

https://hal.science/hal-03971089

Submitted on 6 Feb 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
Selected Papers from the 3rd Radiocarbon in the Environment Conference, Gliwice, Poland, 5–9 July 2021
© The Author(s), 2022. Published by Cambridge University Press for the Arizona Board of Regents on behalf of the University of Arizona. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.

# IDENTIFICATION IN PORE WATERS OF RECYCLED SEDIMENT ORGANIC MATTER USING THE DUAL ISOTOPIC COMPOSITION OF CARBON (δ¹³C AND Δ¹⁴C): NEW DATA FROM THE CONTINENTAL SHELF INFLUENCED BY THE RHÔNE RIVER

J-P Dumoulin<fnref n="1" />* • C Rabouille<fnref n="2" /> • S Pourtout<fnref n="2" /> • B Bombled<fnref n="2" /> • B Lansard<fnref n="2" /> • I Caffy<fnref n="1" /> • S Hain<fnref n="1" /> • M Perron<fnref n="1" /> • M Sieudat<fnref n="1" /> • B Thellier<fnref n="1" /> • E Delqué-Količ<fnref n="1" /> • C Moreau<fnref n="1" /> • L Beck<fnref n="1" />

**ABSTRACT.** Estuaries and deltas are crucial zones to better understand the interactions between continents and oceans, and to characterize the mineralization and burial of different sources of organic matter (OM) and their effect on the carbon cycle. In the present study, we focus on the continental shelf of the northwest Mediterranean Sea near the Rhône river delta. Sediment cores were collected and pore waters were sampled at different depths at one station (Station E) located on this shelf. For each layer, measurements of dissolved inorganic carbon concentration (DIC) and its isotopic composition (δ¹³C and Δ¹⁴C) were conducted and a mixing model was applied to target the original signature of the mineralized OM. The calculated δ¹³C signature of the mineralized organic matter is in accordance with previous results with a δ¹³COM of marine origin that is not significantly impacted by the terrestrial particulate inputs from the river. The evolution with depth of Δ¹⁴C shows two different trends indicating two different Δ¹⁴C signatures for the mineralised OM. In the first 15 cm, the mineralized OM is modern with a Δ¹⁴COM = 100 ± 17‰ and corresponds to the OM produced during the nuclear period of the last 50 years. Deeper in the sediment, the result is very different with a depleted value Δ¹⁴COM = -172 ± 60‰ which corresponds to the pre-nuclear period. In these two cases, the marine substrate was under the influence of the local marine reservoir effect with more extreme Δ¹⁴C results. These differences can be largely explained by the influence of the river plume on the local marine DIC during these two periods.

**KEYWORDS:** pore waters, radiocarbon AMS dating, rivers, sediments, stable isotopes.

## INTRODUCTION

A better understanding of the carbon cycle at the land-ocean interface is crucial to assess the influence of terrestrial inputs and the anthropic impact on the coastal seas (Bauer et al. 2013; Bianchi and Allison 2009). Deltas are particularly active in biogeochemical cycles because, despite their small surface areas compared to the oceans (0.3%), they represent nearly half of all the carbon buried in oceanic sediments (Berner 1982). These large accumulations of organic matter (OM) are accompanied by high mineralization rates and are a major source of CO₂ production, representing one fourth of the CO₂ absorbed by the whole ocean (Cai 2011). It is thus essential to properly constrain and understand the characteristics of OM deposited in deltaic sediments: its origin (land, river, sea), its nature (soils, vegetation, algae, phytoplankton) and even its age and reactivity (old and refractory or labile and energetic). The study of combined isotopic signatures (δ¹³C and Δ¹⁴C) of DIC in pore waters is a powerful tool to assess the sources and the mineralization mechanisms of OM deposited in deltas (Bauer et al. 1995; Aller et al. 2008). The δ¹³C signatures of the dissolved organic carbon (DIC) of sediment pore waters provide information about the origin of the OM mineralized in the sediment, with different signatures for terrestrial sources (δ¹³C around -25‰) and marine sources (δ¹³C around -20‰) (Goñi et al. 1998;

<footnote n="1">¹Laboratoire de Mesure du Carbone 14 (LMC14), LSCE/IPSL, CEA-CNRS-UVSQ, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France</footnote>

<footnote n="2">²Laboratoire des Sciences du Climat et de l'Environnement, UMR 8212 CEA-CNRS-UVSQ et IPSL, Université Paris-Saclay, 91190 Gif sur Yvette, France</footnote>

<footnote>*Corresponding author. Email: Jean-Pascal.Dumoulin@lsce.ipsl.fr</footnote>
Burdige 2006). The Δ¹⁴C signatures provide information about the age of the different OM but can also be used as a tracer for the freshly produced OM of the river (Aller et al. 2008; Aller and Blair 2004). The anthropogenic activities and the presence of nuclear plants on the Rhône river can lead to enriched Δ¹⁴C-OM signatures (Eyrolle et al. 2015) which can be used to follow the fate of the OM in the different pools of the system originating from the river, delta, or ocean (Pozzato et al. 2018; Dumoulin et al. 2018).

In these previous studies we showed that sedimentation and mineralization of terrestrial OM is very high in the proximal zone near the river mouth, while further offshore on the continental shelf, the organic matter sedimentation is dominated by the mineralization of marine OM. However, the data for the continental shelf were limited in these first studies, and a new sampling cruise, “MissRhoDia II”, was carried out on the RV Tethys II in May 2018. Here we report on these new results from the continental shelf sediments and investigate OM mineralization over a longer period of time (one century) in this zone under the influence of the river plume dominated by autochthonous production.

# MATERIALS AND METHODS

## Study Site

The Rhône river is one of the largest French rivers and provides the major input of freshwater and OM to the western Mediterranean Sea. The sources of riverine inputs discharged to the delta depend on the hydrological regime which can go from 700 m³ s⁻¹ during low flow periods to over 3000 m³ s⁻¹ during flood periods. Soils of various ages, vegetation debris and river phyto-planktonic production are mixed and carried by the river to the prodelta (Harmelin-Vivien et al. 2010; Cathalot et al. 2013; Copard et al. 2018). Large amounts of sediment accumulate in the river mouth at a very fast sedimentation rate of around 30 to 40 cm y⁻¹ (Charmasson et al. 1998) whereas further from the coast, on the continental shelf, the sedimentation rate is less than <0.3 cm y⁻¹ (Miralles et al. 2005).

The ¹⁴C activities in the Rhône river have varied over the last century with two main periods (Dumoulin et al. 2018). During the pre-nuclear period, before 1960, the Rhône river had a freshwater reservoir effect (FRE) for radiocarbon in dissolved inorganic carbon (DIC) because of its multiple tributaries and the dissolution of old carbonates from its drainage basin. The FRE was evaluated at Δ¹⁴C = -114‰, close to the value found for its main tributary, the Durance river which is known to be depleted in Δ¹⁴C with a Δ¹⁴C = -138‰ (Jean-Baptiste et al. 2018). The river OM signature was also impacted by this FRE. On the marine side, the DIC and consequently the marine OM, such as phytoplankton, showed a value of Δ¹⁴C = -50‰ to Δ¹⁴C = -60‰ due to the marine reservoir effect in this zone (Siani et al. 2000; Tisnérat-Laborde et al. 2013). Then, during the last 50 years, after the atmospheric nuclear bomb tests and the penetration of the bomb radiocarbon into the river and the ocean, enriched values were measured in the Mediterranean Sea surface waters with Δ¹⁴C = +90‰ around 1970 and a decrease to Δ¹⁴C = +60‰ in 2010 (Tisnérat-Laborde et al. 2013; Ayache et al. 2017). After 1980, nuclear energy developed in France. Today, the Rhône River is one of the most nuclearized rivers in Europe. Four nuclear plants are located along the Rhône and legally discharge cooling water in the river (Eyrolle et al. 2015; Jean-Baptiste et al. 2018). It has already been reported that DIC from the river can contain enriched ¹⁴C values which can reach a maximum of Δ¹⁴C = +393‰ with an average value of Δ¹⁴C = 60‰ (Jean-Baptiste et al. 2018). As a consequence, ¹⁴C can be used as a tracer of the fresh water DIC and the freshly produced OM such as
<figure><img src="image_2.png" /><figcaption>Figure 1 Map of the Rhône prodelta with the sampling station E.</figcaption></figure>

phytoplankton in the ecosystem of the prodelta and the continental shelf (Dumoulin et al. 2018).

In this study, we focused on station E (Figure 1), located in the distal zone of the continental shelf (43°13′12N 4°41′54E, 75 m depth at 17 Km from the River mouth, Cathalot et al. 2013). This station is further from the coast than station C (8.5 km) or D (13 km) studied in our previous work (Pozzato et al. 2018; Dumoulin et al. 2018). Station C provided interesting but scarce data on pore water isotopic signatures. It was important to better constrain the distal zone to refine our conclusions for the continental shelf. At this station, we chose one sediment core and collected pore waters at different depths to analyze the evolution of DIC, $\delta^{13}$C-DIC and $\Delta^{14}$C-DIC.

## Sampling

Sampling was carried out in May 2018. Interface sediment cores of nearly 40 cm were collected with a gravity corer (Uwitec) using PVC transparent tubes. At station E, with the sedimentation rate of the distal zone around $0.3 \text{ cm y}^{-1}$ a core represents between 100 or 150 years of deposition (Miralles et al. 2005). Pore waters were extracted using 0.2 µm rhizons (Rhizosphere Research Products, Seeberg-Elverfeldt et al. 2005) with a 20 mL rubber free syringe. The pore water samples were then sealed in a 15 mL Pyrex glass
ampoules and stored frozen at –20°C until ¹⁴C analysis. A small aliquot (2mL) was placed in an individual vial with mercuric chloride (HgCl₂) for ¹³C analysis. DIC was measured by the calibrated gas pressure gauge at the end of the CO₂ extraction-purification line. A second core was sampled to collect pore waters for more precise DIC measurements. These were also poisoned using mercuric chloride.

Bottom water samples were collected at around one meter from the bottom using a Niskin
bottle, poisoned with mercuric chloride (HgCl₂) and stored in a fridge.

**Measurements**

DIC Analysis

Inorganic carbon measured in pore waters and bottom waters refers to dissolved inorganic carbon (DIC) because during extraction from the sediment, the water passed through the 0.2 µm ceramic filter of the rhizon and the solid particles were removed. The DIC analysis was carried out with an Apollo Scitech Dissolved Inorganic Carbon Analyzer with a LI-COR CO₂ detector (Rassmann et al. 2016). The uncertainty of this analysis is ±10 µmol/L which is around 0.5%.

¹³C Analysis

The ¹³C aliquot was measured at GEOTOP-UQAM using an Isoprime 100-DI with a microgas system in continuous flow for both pore waters and bottom waters. Samples were acidified and heated at 60°C for 1 hour to ensure complete separation of CO₂, which was then introduced in the IRMS for isotope ratio measurements. The precision of the isotopic measurements is ± 0.1‰.

¹⁴C Analysis

The CO₂ extraction from the pore water and bottom water DIC was performed with the CO₂ extraction line installed at the LMC14 as described in Dumoulin et al. (2013). The water samples are introduced through an air-tight septum in the line and 2 mL of 85% phosphoric acid (H₃PO₄) is added to react with the water DIC and produce CO₂. A helium flow is used to push the CO₂ through the line and two water traps at –78°C remove the water from the gas. A liquid nitrogen trap at –190°C is used to collect the sample in a sealed tube before the graphitization step.

The ¹⁴C activities were calculated with the standard of “oxalic acid II”. Radiocarbon calculations were performed using the Mook and van der Plicht method (1999). The background correction was made with C1 AIEA samples and the Δ¹⁴C uncertainty values are ± 3‰ at a confidence interval of 1sigma. The ¹⁴C results are given in Δ¹⁴C and recalculated to 1950 with the correction of the delay between sampling year and measurement year (Mook and van der Plicht 1999).

**Graphitization and AMS Measurements**

The CO₂ was reduced to graphite at the LMC14 laboratory using hydrogen and iron powder at 600°C (Vogel et al. 1984; Dumoulin et al. 2017), then pressed into an aluminum cathode and finally loaded into the ion source of the ARTEMIS facility for the AMS measurement (Moreau et al. 2013, 2020).
<table><caption>Table 1 DIC, δ¹³C-DIC and Δ¹⁴C-DIC of pore waters at different depths. BW represents the bottom water sample. NA indicates samples which were not measurable. Uncertainties are 0.1‰ for δ¹³C and ± 3‰ for Δ¹⁴C.</caption><thead><tr><th>SacA</th><th>Sample</th><th>Depth (cm)</th><th>DIC (mM)</th><th>δ¹³C (‰)</th><th>pMC</th><th>Δ¹⁴C (‰)</th></tr></thead><tbody><tr><td>56788</td><td>BW</td><td>—</td><td>2.38</td><td>-0.30</td><td>103.10 ± 0.36</td><td>31.04</td></tr><tr><td>56789</td><td>2 cm</td><td>2</td><td>2.90</td><td>-1.30</td><td>103.33 ± 0.34</td><td>33.26</td></tr><tr><td>56790</td><td>4 cm</td><td>4</td><td>2.63</td><td>-0.72</td><td>103.99 ± 0.34</td><td>39.88</td></tr><tr><td>56791</td><td>8 cm</td><td>8</td><td>3.05</td><td>-2.01</td><td>104.62 ± 0.36</td><td>46.21</td></tr><tr><td>56792</td><td>10 cm</td><td>10</td><td>3.32</td><td>-6.00</td><td>105.50 ± 0.34</td><td>54.98</td></tr><tr><td>56793</td><td>14 cm</td><td>14</td><td>3.60</td><td>-1.48</td><td>105.19 ± 0.36</td><td>51.94</td></tr><tr><td>56794</td><td>18 cm</td><td>18</td><td>3.88</td><td>-5.73</td><td>102.32 ± 0.34</td><td>23.23</td></tr><tr><td>56795</td><td>22 cm</td><td>22</td><td>4.67</td><td>NA</td><td>101.42 ± 0.38</td><td>14.18</td></tr><tr><td>56796</td><td>26 cm</td><td>26</td><td>4.96</td><td>NA</td><td>99.63 ± 0.34</td><td>-3.72</td></tr><tr><td>56797</td><td>34 cm</td><td>34</td><td>4.96</td><td>-8.52</td><td>96.94 ± 0.33</td><td>-30.59</td></tr></tbody></table>

## Mixing Model

At each depth in the sediment, pore water DIC (DIC<sub>pore</sub>) is a mix of bottom water DIC (DIC<sub>BW</sub>, the original water trapped in the sediment) and DIC originating from the mineralized OM. To calculate the original isotopic signature of the organic matter mineralized in the sediment (δ¹³C<sub>OM</sub>, Δ¹⁴C<sub>OM</sub>) a mixing model similar to that of Bauer et al. (1995) was applied. The calculations were detailed in our previous studies (Pozzato et al. 2018; Dumoulin et al. 2018).

$$ \begin{gathered} \delta^{13}C_{\text{pore}} = x * \delta^{13}C_{\text{BW}} + (1-x) * \delta^{13}C_{\text{OM}} \quad \text{with} \quad x = \frac{\text{DIC}_{\text{BW}}}{\text{DIC}_{\text{pore}}} \\ \& \quad \downarrow \\ \delta^{13}C_{\text{pore}} * \text{DIC}_{\text{pore}} = \delta^{13}C_{\text{OM}} * \text{DIC}_{\text{pore}} + \text{DIC}_{\text{BW}} * \delta^{13}C_{\text{BW}} - \text{DIC}_{\text{BW}} * \delta^{13}C_{\text{OM}} \end{gathered} $$

And similarly:

$$ \Delta^{14}C_{\text{pore}} * \text{DIC}_{\text{pore}} = \Delta^{14}C_{\text{OM}} * \text{DIC}_{\text{pore}} + \text{DIC}_{\text{BW}} * \Delta^{14}C_{\text{BW}} - \text{DIC}_{\text{BW}} * \Delta^{14}C_{\text{OM}} $$

With this equation, the slope of Δ¹⁴C<sub>pore</sub> * DIC<sub>pore</sub> versus DIC<sub>pore</sub> provides the isotopic signature of the OM mineralized in the sediment Δ¹⁴C<sub>OM</sub> and similarly for δ¹³C<sub>OM</sub> with the δ¹³C*DIC versus DIC<sub>pore</sub>). The reported uncertainty on the OM signature was calculated as the uncertainty on the slope due to linear regression.

# RESULTS

## Sediment Pore Waters

The DIC concentration increases with depth (Table 1) with nearly 100 µmol/L per cm, reaching twice the DIC concentration of the bottom water at 25 cm (from 2.38 mmol/L to 4.96 mmol/L).

The δ¹³C-DIC values follow the opposite trend with a large decrease with depth from -0.3‰ for the bottom water to -8.52‰ at 34 cm. Unfortunately, two results are missing because two vials were broken before the δ¹³C analysis.
<figure><img src="image_3.png" /><figcaption>Figure 2 Evolution of the pore water DIC, its δ¹³C and Δ¹⁴C signature with depth in the sediment at station E. The bottom water value is reported above the sediment-water interface (horizontal line at 0 cm).</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>Figure 3 Mixing model for Δ¹⁴C and δ¹³C. The slope of the mixing curve provides the isotopic signature of the mineralized OM. For the Δ¹⁴C mixing model, the blue part (closed symbols) corresponds to the first 15 cm (nuclear period, see text) and the brown part (open symbols) corresponds to depths below 15 cm in the pre-nuclear period. (Please see online version for color figures.)</figcaption></figure>

The Δ¹⁴C-DIC in the pore waters follows two trends (Figure 2). In the first 15 cm, the activity increases from 30‰ to 55‰ but deeper down, it decreases from 55‰ to -30‰.

**Original isotopic signature of the organic matter mineralized in sediment pore waters:**

Figure 3 shows the mixing model and the isotopic signature of the OM mineralized in the sediment.

With the growth of DIC concentration in the pore water, the δ¹³Cₚore\*DICₚore decreases. The slope of the mixing line for δ¹³C gives a δ¹³Cₚor signature around -16.2 ± 3‰ for the OM mineralized in the sediment. For the Δ¹⁴C mixing model, the plot shows two different trends corresponding to different gradients observed on the Δ¹⁴C-DIC distribution. For the
first part of the graph, the slope of the mixing model indicates $\Delta^{14}C_{OM} = +99.7 \pm 17\%$ and for the second part of the graph, the slope of indicates $\Delta^{14}C_{OM} = -171.7 \pm 60\%$.

# DISCUSSION

## Mineralization of the Organic Matter and DIC Increase

The measured DIC concentrations in the pore waters of continental shelf sediments at station E double from the surface sediment to the bottom of the core at 34 cm. This DIC concentration increase with depth is the result of oxic and anoxic OM mineralization in the sediment conducted by microorganisms (Rassmann et al. 2016, 2020). On the continental shelf, the low concentration of nitrate and the low production of reduced metals in the pore waters highlight that oxic mineralization and sulfate reduction are the major organic carbon mineralization pathways (Pastor et al. 2011; Ait Ballagh et al. 2021). It follows the equation (Rassmann et al. 2016, 2020):

$$CH_2O + O_2 \rightarrow HCO_3^- + H^+ \text{ and } 2CH_2O + SO_4^{2-} \rightarrow 2HCO_3^- + H_2S$$

For the anoxic part of the sediments (directly after the first centimeters), carbonate dissolution can be excluded as a source of DIC since in this zone, carbonate minerals (calcite and aragonite) are supersaturated and cannot therefore dissolve (Rassmann et al. 2016).

## Isotopic signature of the mineralized OM

The slope of Figure 3 provides a $\delta^{13}C_{OM} = -16 \pm 3\%$ for the signature of the OM mineralized at station E, providing information about its origin. This result is in the lower range of the $\delta^{13}C$ signature corresponding to the mineralization of a marine substrate such as marine phytoplankton which commonly have a $\delta^{13}C$ signature between -19 to -21‰ (Harmelin-Vivien et al. 2010). This information confirms our previous results (Lansard et al. 2009; Cathalot et al. 2013; Pozzato et al. 2018; Dumoulin et al. 2018) indicating the limited amount of terrestrial organic matter in continental shelf sediments. The majority of OM delivered by the river settles near the river mouth, in the proximal zone. Further out at sea on the continental shelf, local marine OM is predominant. It is noteworthy that the input of C4 plants which have a heavier $\delta^{13}C$ signature of around -15‰ can be ruled out based on the $\delta^{13}C$ data of lignin phenols previously obtained (Cathalot et al. 2013). The absence of a large delivery of river particles to the continental shelf also explains why the sedimentation rates drastically decrease in this zone compared to those in the prodelta which receives large amounts of terrestrial particles and OM from the continent.

The $\Delta^{14}C$ isotopic signatures of the OM mineralized in the sediment (Figure 3) show two different trends. Line #1, in the first 15 cm, has a slope which gives a $\Delta^{14}C_{OM} = 100 \pm 17\%$. This value of $\Delta^{14}C_{OM}$ indicates modern organic matter with a slight $\Delta^{14}C$ enrichment and fits well with the post bomb period expected for the top of the core. A sedimentation rate of less than 0.3 cm/yr (Miralles et al. 2005) represents around 50-60 years for the first 15 cm and encompasses the $^{14}$C bomb peak period and the initiation of nuclear activity in the Rhône area which started in the 1980s. This $\Delta^{14}C_{OM}$ result is probably linked to the deposition of enriched phytoplankton grown in marine surface waters during the nuclear period which is mineralized in the sediment and which transfers its isotopic signature to the pore water DIC. Evidence for this is the fact that the $\Delta^{14}C$ values obtained in the NW Mediterranean Sea surface were around +90‰ in the 1970s and decreased to 60‰ in 2010
(Tisnérat-Laborde et al. 2013; Ayache et al. 2017), which is consistent with the signature found in the OM for the nuclear period. As our results are slightly higher than the values of the surface ocean, another possible source of enriched DIC originating from the river plume is the impact of the nuclear power plants along the Rhône during the last 40 years. DIC of the river plume with peaks of $\Delta^{14}C = +393\text{‰}$ and an average of $60\text{‰}$ (Jean-Baptiste et al. 2018) can be mixed with DIC of the surface waters and can be used by marine phytoplankton together with plume dissolved nutrients, creating an $\Delta^{14}C$ enrichment of marine OM.

A different trend is shown by line #2 for the signature of the organic matter mineralized under 15 cm. The slope of the mixing model gives a strongly depleted result with $\Delta^{14}C_{OM} = 172 \pm 60\text{‰}$. These sediment layers under 15 cm correspond to a period of time spanning from before 1900 to 1960 as estimated from the sedimentation rate (Miralles et al. 2005), i.e., the pre-nuclear period. Before the anthropogenic nuclear activities of the last 60 years, the $\Delta^{14}C$ signature of marine OM was directly linked to the local marine reservoir effect of the Mediterranean Sea evaluated around -50 to -60‰ (Siani et al. 2000; Tisnérat-Laborde et al. 2013). The value found for marine plankton degrading in these sediments is lower than these DIC signatures of the Mediterranean Sea, so it certainly indicates a larger influence of $^{14}$C depleted waters. As indicated above, the river plume brings a large amount of fresh water and dissolved nutrients to the Mediterranean Sea, modifying the DIC of the surface water. During the pre-nuclear period, the freshwater reservoir effect (FRE) for radiocarbon in DIC was strong in the Rhône River basin. It was generated by the dissolution of old carbonates from the drainage basin which created a freshwater reservoir effect and an apparent $^{14}$C aging of DIC (see for the Loire River, Coularis et al. 2016). The FRE measured on the Rhône river DIC was evaluated at $\Delta^{14}C = -114\text{‰}$, a value indicating the strong influence of one tributary, the Durance river which is known to be depleted in $\Delta^{14}C$ as indicated by a value obtained in May 2014 ($\Delta^{14}C = -138\text{‰}$; Jean-Baptiste et al. 2018). As the DIC signature of the river is partly transferred to the river plume, the local marine OM could have been more depleted than the average marine surface waters. Given the large uncertainty of the result found for $\Delta^{14}C_{OM}$ (-171.7 $\pm$ 60‰), there is an overall compatibility of the mineralized OM signature with the possible signature of river-influenced marine surface waters.

A final possibility to explain such a depleted result for the OM mineralized in the sediments below 15 cm is a possible “priming effect” (Bianchi 2011 and references therein). This phenomenon occurs when substrates with different labilities co-exist in the sediment. The labile OM is used by the micro-, meio- and macro-organisms to co-metabolize more refractory substrates. At station E in the distal zone, this “priming effect” is possible with the mixture of different sources of substrates. If the sources of energetic OM are limited, the fauna and the different organisms have to mineralize more refractory OM. Cathalot et al. (2013) described a dominant aged OM fraction in shelf sediments including station E, with a high content of black carbon (up to 50%) and a limited contribution of fresh marine OM, providing an old $\Delta^{14}C$ signature of particulate OM ($\Delta^{14}C = -312\text{‰}$). Therefore, the mineralization of a small part of this refractory material could explain part of the $\Delta^{14}C_{OM}$ depleted results in the deeper layers of pore waters.

## CONCLUSION

This new study of the continental shelf of the Rhône River prodelta has shed light on the origin and reactivity of organic matter present in the sediments and the coastal carbon cycle in river-
dominated ocean margins. The large increase in pore water DIC concentration indicates active OM mineralization by micro and macro fauna over a period of 100–150 yr. This process transfers the OM signature to the pore water DIC. The $\delta^{13}C$ signature of the pore water DIC indicates that mostly marine organic matter is mineralized in these sediments. It also confirms that a limited particulate substrate from the continent is mineralized in continental shelf sediments, confirming its efficient trapping in the proximal zone. The $\Delta^{14}C$ signature indicates two behaviors that coincide with the post-bomb/pre-bomb periods. Mineralized OM in the first half of the 20th century (below 15 cm) shows a depleted $^{14}C$ composition in line with low $\Delta^{14}C$ marine DIC values corresponding to the local marine reservoir effect and additions of river plume waters with lower $\Delta^{14}C$-DIC values impacted by a freshwater reservoir effect. Contrarily, mineralized OM of the top of the core is enriched and $\Delta^{14}C$ results fit well with the expected values for the nuclear period impacted first by the penetration of bomb radiocarbon in the ocean and later by the anthropic nuclear activities present along the Rhône River. The larger contrast between the $\Delta^{14}C$ signature of mineralized OM between the two periods (–172 ± 60‰ to +100 ± 17‰) and that of the marine surface waters (–60‰ to 90‰; Tisnérat-Laborde et al. 2013) could be explained by a significant influence of the river plume waters and its DIC on marine phytoplankton production. This could have contributed a more depleted DIC during the pre-nuclear period and a more enriched DIC during the nuclear period.

Overall, the $\delta^{13}C$ and $\Delta^{14}C$ signature of pore water DIC confirms our previous findings that reactive organic matter mineralized in sediments imprints its signature on the pore waters. We show here that this signature can be maintained on timescales of the order of a century.

## ACKNOWLEDGMENTS

The authors gratefully acknowledge the captains and crews of the RV Thetys II for their help during the MissRhoDia II sampling campaign. The authors thank A. Villedieu and L. Brethous who provided technical help during the expedition and in the laboratory. This work was supported by the INSU/EC2CO-MissRhoDia project, the Sulfat_IsoMic project funded by INSU/TELLUS Interrvie and IPSL, the BioGeoMethane project funded by Emergence-SPU-Paris Saclay and the French State program “Investissement d’avenir” run by the National Research Agency (AMORAD project ANR-11-RSNR-0002). The LMC14 is funded by five French organizations: CEA, CNRS, IRD, IRSN, and MCC. This is LSCE contribution number 7813.

## REFERENCES

Ait Ballagh FE, Rabouille C, Andrieux-Loyer F, Soetaert K, Lansard B, Bombled B, Monvoisin G, Elkalay K, Khalil K. 2021. Spatial variability of organic matter and phosphorus cycling in Rhône River prodelta sediments (NW Mediterranean Sea, France): a data-model approach, Estuar. Coast. doi: 10.1007/s12237-020-00889-9.

Aller RC, Blair NE. 2004. Early diagenetic remineralization of sedimentary organic C in the Gulf of Papua deltaic complex (Papua New Guinea): net loss of terrestrial C and diagenetic fractionation of C isotopes, Geochim.

Cosmochim. Acta 68:1815-1825. doi: 10.1016/j.gca.2003.10.028.

Aller RC, Blair NE, Brunskill GJ. 2008. Early diagenetic cycling, incineration, and burial of sedimentary organic carbon in the central Gulf of Papua (Papua New Guinea). Journal of Geophysical Research-Earth Surface 113.

Ayache M, Dutay J-C, Mouchet A, Tisnérat-Laborde N, Montagna P, Tanhua T, Siani G, Jean-Baptiste P. 2017. High resolution regional modeling of natural and anthropogenic radiocarbon in the Mediterranean Sea. Biogeosciences. doi: 10.5194/bg-2016-444.
Bauer JE, Cai WJ, Raymond PA, Bianchi ST, Hopkinson CS, Regnier PAG. 2013. The changing carbon cycle of the coastal ocean. Nature 504:61–70.

Bauer JE, Reimers CE, Druffel ERM, Willimas PM. 1995. Isotopic constraints on carbon exchange between deep ocean sediments and sea water. Nature 373:686–689.

Berner RA. 1982. A rate model for organic decomposition during bacterial sulfate reduction in marine sediments, Biogéochimie de la matière organique à l'interface eau-sédiment marin. Colloques internationaux du CNRS, Paris. p. 35-44.

Bianchi ST, Allison MA. 2009. Large-river delta-front estuaries as natural “recorders” of global environmental change. Proc. Nat. Acad. Sci. 106:8085–8092.

Bianchi ST. 2011. The role of terrestrially derived organic carbon in the coastal ocean: a changing paradigm and the priming effect. Proc. Nat. Acad. Sci. 108:19473–19481.

Burdige DJ. 2006. Geochemistry of marine sediments. Princeton (USA): Princeton University Press. 609 p.

Cai WJ. 2011. Estuarine and coastal ocean carbon paradox: CO₂ sinks or sites of terrestrial carbon incineration? Ann. Rev. Mar. Sci. 3:123–145.

Cathalot C, Rabouille C, Tisnerat-Laborde N, Toussaint F, Kerhervé P, Buscail R, Loftis K, Sun MY, Tronczynski J, Azoury S, et al. 2013. The fate of river organic carbon in coastal areas: A study in the Rhone River delta using multiple isotopic (delta C-13, Delta C-14) and organic tracers. Geochim. Cosmochim. Acta 118:33–55.

Charmasson S, Bouisset P, Radakovitch O, Pruchon AS, Arnaud M. 1998. Longcore profiles of Cs-137, Cs-134, Co-60 and Pb-210 in sediment near the Rhone River (Northwestern Mediterranean Sea). Estuaries 21:367–378.

Copard Y, Eyrolle F, Radakovitch O, Poirel A, Raimbault P, Gairoard S, Di-Giovanni C. 2018. Badlands as a hot spot of petrogenic contribution to riverine particulate organic carbon to the Gulf of Lion (NW Mediterranean Sea). Earth Surf. Process. Landform(bu). doi: 10.1002/esp.4409.

Coularis C, Tisnérat-Laborde N, Pastor L, Siclet F, Fontugne M. 2016. Temporal and spatial variations of freshwater reservoir ages in the Loire River watershed. Radiocarbon 58: 549–563.

Dumoulin JP, Caffy I, Comby-Zerbino C, Delqué-Količ E, Hain S, Massault M, Moreau C, Quiles A, Setti, V, Souprayen C, Tannau J-F, Thellier B, Vincent J. 2013. Development of a line for dissolved inorganic carbon extraction at LMC14 Artemis Laboratory in Saclay, France. Radiocarbon 55:1043–1049.

Dumoulin J-P, Comby-Zerbino C, Delqué-Količ E, Moreau C, Caffy I, Hain S, Perron M, Thellier B, Setti V, Berthier B, Beck L. 2017. Status report on sample preparation protocols developed at the LMC14 Laboratory, Saclay, France: from sample collection to ¹⁴C AMS measurement. Radiocarbon 59:713–726.

Dumoulin J-P, Pozzato L, Rassman J, Toussaint F, Fontugne M, Tisnérat-Laborde N, Beck L, Caffy I, Delqué-Količ E, Moreau C, Rabouille C. 2018. Isotopic dignature (δ¹³C, Δ¹⁴C) of DIC in sediment pore waters: an example from the Rhone River Delta. Radiocarbon 60(5):1465–1481. doi: 10.1017/RDC.2018.111.

Eyrolle F, Antonelli C, Renaud P, Tournieux D. 2015. Origins and trend of radionuclides within the lower Rhône River over the last decades. Radioprotection 50:27–34.

Goñi MA, Ruttenberg KC, Eglinton TI. 1998. A reassessment of the sources and importance of land-derived organic matter in surface sediments from the Gulf of Mexico. Geochimica Cosmochimica Acta 62:3055–3075.

Harmelin-Vivien M, Dierking J, Banaru D, Fontaine MF, Arlhac D. 2010. Seasonal variation in stable C and N isotope ratios of the Rhone River inputs to the Mediterranean Sea (2004–2005). Biogeochemistry 100:139–150.

Jean-Baptiste P, Fontugne M, Fourré E, Marang L, Antonelli C, Charmasson S, Siclet F. 2018. Tritium and radiocarbon levels in the Rhône river delta and along the French Mediterranean coastline. Journal of Environmental Radioactivity 187:53–64.

Lansard B, Rabouille C, Denis L, Grenz C. 2009. Benthic remineralization at the land-ocean interface: a case study of the Rhone River (NW Mediterranean Sea). Estuarine Coastal and Shelf Science 81(4):544–554. doi: 10.1016/j.ecss. 2008.11.025.

Miralles J, Radakovitch O, Aloisi JC. 2005. Pb-210 sedimentation rates from the northwestern Mediterranean margin. Mar. Geol. 216:155–167.

Mook WG, van der Plicht J. 1999. Reporting C-14 activities and concentrations. Radiocarbon 41:227–239.

Moreau C, Caffy I, Comby C, Delqué-Količ E, Dumoulin J-P, Hain S, Quiles A, Setti V, Souprayen C, Thellier B, et al. 2013. Research and development of the Artemis ¹⁴C AMS facility: status report. Radiocarbon 55(2–3): 331–337.

Moreau C et al. 2020. ARTEMIS, the ¹⁴C AMS facility of the LMC14 National Laboratory: a status report on quality control and microsample procedures. Radiocarbon 62(6):1755–1770. doi: 10.1017/RDC.2020.7.

Pastor LC, Cathalot B, Deflandre E, Viollier K, Soetaert FJR, Meysman C, Ulses E, Metzger E, Rabouille C. 2011. Modeling biogeochemical
processes in sediments from the Rhone River prodelta area (NW Mediterranean Sea). Biogeosciences 8(5):1351–1366. doi: 10.5194/bg-8-1351-2011.

Pozzato L, Rassmann J, Lansard B, Dumoulin J-P, van Brugel P, Rabouille C. 2018. Origin of remineralized organic matter in sediments from the Rhone River prodelta (NW Mediterranean) traced by Δ¹⁴C and δ¹³C signatures of pore water DIC. Progress in Oceanography 163:112–122.

Rassmann J, Eitel E, Cathalot C, Brandily C, Lansard B, Taillefert M, Rabouille C. 2020. Benthic alkalinity and DIC fluxes in the Rhône River prodelta generated by decoupled aerobic and anaerobic processes. Biogeosciences 17:13–33. doi: 10.5194/bg-17-13-2020.

Rassmann J, Lansard B, Pozzato L, Rabouille C. 2016. Carbonate chemistry in sediment pore waters of the Rhône River delta driven by early

diagenesis (NW Mediterranean). Biogeosc. 13:5379–5394. doi: 10.5194/bg-13-5379-2016.

Seeberg-Elverfeldt J, Schlüter M, Feseker T, Kölling, M. 2005. Rhizon sampling of pore waters near the sediment-water interface of aquatic systems. Limnol. Oceanogr.-Methods 3:361–371.

Siani G, Paterne M, Arnold M, Bard E, Métivier B, Tisnérat N, Bassinot F. 2000. Radiocarbon reservoir ages in the Mediterranean Sea and Black Sea. Radiocarbon 42(2):271–280.

Tisnérat-Laborde N, Montagna P, Frank N, Siani G, Silenzi S, Paterne M. 2013. A high-resolution coral-based D¹⁴C record of surface water processes in the western Mediterranean Sea. Radiocarbon 55(2–3):1917–1630.

Vogel JS, Southon JR, Nelson DE, Brown TA. 1984. Performance of catalytically condensed carbon for use in accelerator mass spectrometry. Nuclear Instruments and Methods in Physics Research B 5(2):289–93.