# Impacts of Ocean Acidification on Sediment Processes in Shallow Waters of the Arctic Ocean

Frédéric Gazeau, Pieter van Rijswijk, Lara Pozzato, Jack J. Middelburg

► To cite this version:

Frédéric Gazeau, Pieter van Rijswijk, Lara Pozzato, Jack J. Middelburg. Impacts of Ocean Acidification on Sediment Processes in Shallow Waters of the Arctic Ocean. PLoS ONE, 2014, 9 (4), pp.e94068. 10.1371/journal.pone.0094068 . hal-01358325

HAL Id: hal-01358325

https://hal.sorbonne-universite.fr/hal-01358325v1

Submitted on 31 Aug 2016

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 CC BY 4.0 - Attribution - International License
# Impacts of Ocean Acidification on Sediment Processes in Shallow Waters of the Arctic Ocean

Frédéric Gazeau¹,²*, Pieter van Rijswijk³, Lara Pozzato³, Jack J. Middelburg³,⁴

¹ Sorbonne Universités, Université Pierre et Marie Curie-Paris 06, Laboratoire d'Océanographie de Villefranche, Observatoire océanologique, Villefranche/mer, France,
² Centre National de la Recherche Scientifique-Institut National des Sciences de l'Univers, Laboratoire d'Océanographie de Villefranche, Observatoire océanologique,
Villefranche/mer, France,
³ Department of Ecosystem Studies, Royal Netherlands Institute for Sea Research, Yerseke, The Netherlands,
⁴ Faculty of Geosciences, Utrecht
University, Utrecht, The Netherlands

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

## Abstract

Despite the important roles of shallow-water sediments in global biogeochemical cycling, the effects of ocean acidification on sedimentary processes have received relatively little attention. As high-latitude cold waters can absorb more CO₂ and usually have a lower buffering capacity than warmer waters, acidification rates in these areas are faster than those in sub-tropical regions. The present study investigates the effects of ocean acidification on sediment composition, processes and sediment-water fluxes in an Arctic coastal system. Undisturbed sediment cores, exempt of large dwelling organisms, were collected, incubated for a period of 14 days, and subject to a gradient of pCO₂ covering the range of values projected for the end of the century. On five occasions during the experimental period, the sediment cores were isolated for flux measurements (oxygen, alkalinity, dissolved inorganic carbon, ammonium, nitrate, nitrite, phosphate and silicate). At the end of the experimental period, denitrification rates were measured and sediment samples were taken at several depth intervals for solid-phase analyses. Most of the parameters and processes (i.e. mineralization, denitrification) investigated showed no relationship with the overlying seawater pH, suggesting that ocean acidification will have limited impacts on the microbial activity and associated sediment-water fluxes on Arctic shelves, in the absence of active bio-irrigating organisms. Only following a pH decrease of 1 pH unit, not foreseen in the coming 300 years, significant enhancements of calcium carbonate dissolution and anammox rates were observed. Longer-term experiments on different sediment types are still required to confirm the limited impact of ocean acidification on shallow Arctic sediment processes as observed in this study.

**Citation:** Gazeau F, van Rijswijk P, Pozzato L, Middelburg JJ (2014) Impacts of Ocean Acidification on Sediment Processes in Shallow Waters of the Arctic Ocean. PLoS ONE 9(4): e94068. doi:10.1371/journal.pone.0094068

**Editor:** Howard I. Browman, Institute of Marine Research, Norway

**Received** December 2, 2013; **Accepted** March 10, 2014; **Published** April 9, 2014

**Copyright:** © 2014 Gazeau et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This research has received support from the Netherlands Organization of Scientific Research, the Darwin Centre for Biogeosciences and is a contribution to the "European Project on Ocean Acidification" (EPOCA) which received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under grant agreement n° 211384. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: f.gazeau@obs-vlfr.fr

## Introduction

Because about one third of anthropogenic CO₂ emissions (from fossil fuel, cement production and land-use changes) has been stored in the oceans since the industrial revolution [1], surface seawater pH has already declined by ~0.1 unit compared with pre-industrial values [2,3] and it is projected to decrease by 0.36 unit by the end of the century [3]. Decreasing pH levels are expected to have profound impacts on the physiology and metabolism of marine organisms through a disruption of intercellular transport mechanisms [4]. Moreover, this pH decrease will modify the equilibrium between the different forms of dissolved inorganic carbon (Cᵢ) in seawater with an increase of the proportion of both CO₂ and bicarbonate (HCO₃⁻) at the expense of carbonate ions (CO₃^{2-}).

While it has been shown that some photosynthetic organisms will benefit from elevated CO₂ conditions [5], most calcareous organisms have revealed a high sensitivity to the decreasing availability of CO₃^{2-} [6]. A decrease in CO₃^{2-} concentration leads to a reduction in the level of calcium carbonate saturation (Ω) of seawater. While it appears that most of the ocean surface will

stay supersaturated with respect to all CaCO₃ forms, high-latitude cold waters that are naturally more corrosive to CaCO₃ than warm waters, will absorb more CO₂ and become undersaturated with respect to high-magnesian calcite and aragonite, the most soluble forms of CaCO₃, in few decades [7]. Acidification rates in the Arctic Ocean have shown to be faster than in sub-tropical regions [8] with some areas, as in the Canada basin and the Chukchi Sea, already experiencing undersaturated conditions [9,10] due to a combination of sea-ice melting and anthropogenic CO₂ penetration. These effects of freshening and increased carbon uptake in response to sea-ice retreat due to global change make Arctic Ocean surface waters the area in the world that will experience the largest pH declines in the coming decades [11]. Moreover, river runoff to Arctic shelves causes a significant decrease of alkalinity and calcium ions, and therefore lowers saturation states with respect to aragonite and calcite [12].

Shallow-water sediments play important roles in the global carbonate cycle as they represent a large reservoir of CaCO₃ that can react to the decreasing saturation state of seawater, releasing alkalinity to the overlying water column. This dissolution of sedimentary carbonates has certainly played a significant role in
past variations of atmospheric partial pressure of CO₂ (pCO₂) for instance during glacial/interglacial transitions [13]. As the rate of this reaction is both kinetically and physically limited and as the amount of CaCO₃ is not large enough [14], it cannot compensate for the actual very fast increase of atmospheric CO₂ [15,16]. However, it has been suggested that dissolution of metastable carbonate phases could locally buffer changes in carbonate saturation in coral reef areas [17] for both muddy/fine grained [18,19] and permeable sediments [20–22]. As high-latitudes sediments represent a much smaller CaCO₃ reservoir than tropical areas, it is very likely that they will not have the capacity for locally buffering anthropogenic increases in CO₂ [23].

Shallow-water sediments are also significant contributors to many biogeochemical processes in the ocean, for instance providing a large proportion of the nutrients required for primary production via organic matter degradation. Indeed, about a quarter of the organic matter that is exported from the surface of the ocean sinks onto the sediment where 90% of this organic matter is remineralized via diverse oxic or anoxic pathways [24]. Sediments are also important contributors of oceanic organic matter production via photo- or chemo-synthesis [25,26]. Recently, several studies have discussed the role of shallow sediments in releasing alkalinity and therefore acting as a negative feedback to atmospheric CO₂ increase, at least on the regional scale [27–29]. Alkalinity generating anaerobic processes (e.g. denitrification, sulfate reduction etc...) appear as important as calcium carbonate dissolution in buffering anthropogenic CO₂ on the global scale [29]. However, despite these important roles of shallow water sediments in global biogeochemical cycling, the effects of ocean acidification on sedimentary processes have received relatively little attention. For instance, while water-column nitrification (the microbial process of ammonia oxidation to nitrite and nitrate) has been shown, through several studies, to be inhibited at low pH [30–32], only one study so far focused on ocean acidification effects on nitrification within sediments, actually showing a clear resilience of nitrifying microbial communities to this perturbation [32]. In contrast, several recent studies have demonstrated some indirect modifications of shallow sediment-water nutrient fluxes due to ocean acidification impacts on macro-organisms responsible for bioturbation and bioirrigation activities in the sediment [33–37]. Indeed, these organisms are known to enhance nutrient fluxes by increasing the penetration of oxygen (ventilation activity) and nutrients within the sediment and increasing the total sediment-water interface [38]. For instance, Dashfield et al. [36] have shown that within-sediment pH profiles and nematode community structure were significantly affected by a decrease of ~0.5 pH unit through a modification of the activity of burrowing sea urchins.

The present study contributes to the assessment of the future role of shallow-water sediments to global carbon and nutrient cycling by investigating the effects of ocean acidification on sediment composition, processes and sediment-water fluxes in a coastal Arctic fjord. Undisturbed sediment cores were collected, incubated for a period of 14 days, and subject to a gradient of pCO₂ covering the range of values projected for the end of the century (IPCC, 2007). While the present manuscript reports on sediment solid-phase composition and sediment-water fluxes, a companion paper provides information on the evolution of sediment microbial community composition [39].

## Materials and Methods

Sediment cores (n = 20, h = 48 cm, d = 10 cm) were collected by scuba divers at 4 m depth from an area of moderately sorted sandy

mud, close to the harbour in Kongsfjorden, Svalbard (Fig. 1) on May, 10th 2009 and immediately transported to the Kings Bay Marine Laboratory in Ny-Ålesund. No specific permissions were required for sampling at this location. Sediment was ~70% silt and ~20% very fine sand, and contained ~2.5% and ~2.6% of organic and inorganic carbon, respectively. In the laboratory, during 14 days, the sediment cores were kept in darkness and continuously supplied with temperature-regulated *in situ* seawater (~0°C), pumped at a rate of 5 ml min⁻¹ from header tanks (n = 5, 200 l, Fig. 2) in which pCO₂ was controlled and maintained through pure CO₂ bubbling using flat surface, combination pH electrodes (Walchem S650CD) and an automated feedback relay system. Five different pCO₂ levels were considered: *in situ* conditions (no CO₂ bubbling; ~317 µatm), 540, 760, 1120 and 3000 µatm. Four sediment cores were used per treatment.

This perturbation design has been considered in order to encompass the different scenarios that are projected by IPCC climate models for the end of the century (540 to over 1000 ppm depending on the considered emission scenario [40]). The 3000 µatm treatment has been used as an extreme case study.

Corresponding pH levels used for the regulation system were calculated every other day using the R package seacarb (Lavigne and Gattuso, 2010) based on total alkalinity ($A_T$, see below) measurements and desired pCO₂ levels. To ensure an optimal aeration and a better dissolution of CO₂ in seawater, header tanks were continuously and actively bubbled with ambient air. Four cores were supplied with seawater from each header tank. pH was measured every second day in the 5 header tanks using a pH meter (Metrohm, 826 pH mobile) with a glass electrode (Metrohm, electrode plus) calibrated on the total scale using Tris/HCl and 2-aminopyridine/HCl buffer solutions with a salinity of 35.0 at the temperature of the experimental setup (Dickson et al., 2007). Simultaneously to the pH recording, 100 ml of water from the control tank were carefully sampled and filtered on GF/F for $A_T$ measurements that were performed within two days using a potentiometric titration and a Metrohm titrator (Titrando 80). Measurements were carried out in triplicate on 25 ml samples at 21°C and $A_T$ was calculated with a precision of 1–2 µmol kg⁻¹ using a Gran function applied to the pH_T values ranging from 3.5 to 3.0 as described by Dickson et al. [41]. Titrations of standard seawater provided by A. G. Dickson (batch 95) yielded $A_T$ values within 1.7 µmol kg⁻¹ of the nominal value (2216.45 µmol kg⁻¹; standard deviation = 1.85 µmol kg⁻¹; n = 12). All the parameters of the carbonate chemistry, including dissolved inorganic carbon ($C_T$) concentrations were determined from pH_T, $A_T$, temperature and salinity using the R package seacarb (Lavigne and Gattuso, 2012).

On five occasions during the experimental period (14 days), the sediment cores were isolated during 24 h for flux measurements. Before and after incubation, seawater samples were taken for pH_T (see above), $A_T$ (see above), dissolved oxygen (DO) and nutrient analyses: ammonium ($NH_4^+$), nitrate ($NO_3^-$), nitrite ($NO_2^-$), dissolved inorganic phosphorus (DIP) and dissolved silicate (Si). Nutrients were measured with an automated colorimetric method [42] with a precision better than 4% and the accuracy was verified through use of certified standards and participation in QUASI-MEME robin round tests with a Z-score <0.5. DO concentrations were measured using an automated Winkler titration technique with a potentiometric end-point detection. Analyses were performed with a Metrohm redox electrode and a custom built titrator. Reagents and standardizations were similar to those described by Knap et al. [43].

At the end of the experimental period, denitrification rates were measured in each core over a 24 h incubation period using ¹⁵N
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Location of the sampling area and of the experimental set-up. Sediment cores were collected near the harbour of Kongsfjorden and experiments were carried out at the marine laboratory in Ny-Ålesund (78°55′N, 11°56′E). doi:10.1371/journal.pone.0094068.g001</figcaption></figure>

tracers. Enriched sodium nitrate (¹⁵N 98%, Cambridge Isotope Laboratories Inc) was added to core overlying water to reach a final enrichment of 50 at% ¹⁵N for nitrate. Samples of overlying water (20 ml) were taken at 0, 4, 12 and 24 h for ¹⁵N-N₂ measurements [44], which were performed using a Thermo Flash EA 1112 elemental analyzer coupled via a Conflo III interface to an Isotope Ratio Mass Spectrometer (Thermo Delta V Advantage). Denitrification rates were estimated using the following equation:

$$D = \frac{[(^{15}\text{N} - \text{N}_2) \times \text{N}_2 \times \text{h}]}{^{15}\text{N} - \text{NO}_3^-} \quad (1)$$

Where D is denitrification (in µmol m⁻² d⁻¹), ¹⁵N-N₂ is enrichment rate of N₂ in the overlying water with time (in atom % d⁻¹), N₂ is the N₂ concentration in the overlying water (in µmol

m⁻³) which was assumed constant, at saturation and calculated based on [45], h is the water height in the sediment cores (in m) and ¹⁵N-NO₃⁻ is the atom % ¹⁵N-NO₃⁻ in the overlying water. During these incubations, bio-irrigation rates were quantified using the bromide technique [46]. Sodium bromide was added in the overlying water at the start of the incubations (concentration of ~0.9 mmol L⁻¹). Bio-irrigation was estimated as the decrease of bromide concentration with time by sampling 20 ml of the overlying water at 0, 4, 12 and 24 h. Bromide was measured by reverse-phase high-performance liquid chromatography (HPLC).

After the end of the last incubations for denitrification and bio-irrigation measurements, sediment samples were taken at depth intervals of 0–0.5, 0.5–1, 1–2, 2–4 and >4 cm and freeze-dried for subsequent analyses of grain size distribution, pigment concentrations and total/organic carbon and total nitrogen content as well as isotopic compositions (¹³C) of total and organic carbon and isotopic compositions (¹⁵N) of organic nitrogen. Grain-size distribution was assessed with a laser diffraction technique
<figure><img src="image_4.png" /><figcaption>Figure 2. Experimental set-up used in this study. See text for details.</figcaption></figure>

doi:10.1371/journal.pone.0094068.g002

(Malvern Mastersizer 2000). Pigments (chlorophyll a, b and c) were extracted from freeze-dried samples with 10 ml 90/10 acetone/water buffered with ammonium acetate (5%), and analyzed by HPLC. Measurements of total and organic carbon ($C_{tot}$ and $C_{org}$, respectively) and total nitrogen ($N_{tot}$) were made using a Fisons NA-2500 elemental analyzer [47]. The carbon (total and organic) and organic nitrogen isotopic composition of the sediment was determined using a Thermo Flash EA 1112 elemental analyzer coupled via a Conflo III interface to an Isotope Ratio Mass Spectrometer (Thermo Delta V Advantage). The carbon and nitrogen isotope ratios are expressed in the delta notation $δ^{13}C$ and $δ^{15}N$, respectively, where: $δ = (R_{sample}/R_{standard} - 1)×1000$

Results are referred to Vienna PDB for C and to atmospheric nitrogen for N and expressed in units of ‰. Reproducibility of the measurements is better than 0.2 ‰ [48].

## Statistics

Principal components analyses (PCA) were used to test 1) whether the sediment cores used in this study were similar in terms of grain size distribution (5 variables, 96 samples), 2) the effects of $\rho CO_2$ and depth of sampling in the sediment on the composition of the sediment with respect to the different analyzed variables (pigments, $C_{tot}$, $C_{org}$, $N_{tot}$, $δ^{13}C_{tot}$, $δ^{13}C_{org}$ and $δ^{15}N_{org}$; 9 variables, 87 samples) and 3) the effects of the environmental conditions to which the sediment cores were exposed during the incubations ($\rho CO_2$ and temperature) and the effects of time on sediment-water fluxes (10 variables, 90 samples). A one-way ANOVA was used to test the effects of seawater $\rho CO_2$ on denitrification rates. Data were first checked to ensure they conformed to the assumptions of ANOVA (normality: Kolmogorov-Smirnov test and homoscedasticity: Bartlett test). All statistical analyses were conducted using the R software and PCA analyses were performed using the ADE4 package.

## Results

Environmental conditions in the header tanks, from which seawater was pumped and distributed to the sediment cores, are

presented in Table 1. Seawater temperature averaged 0.18±0.53°C during the experimental period. Seawater $pH_T$ was maintained at 8.13±0.08, 7.84±0.11, 7.71±0.21, 7.57±0.08 and 7.17±0.06 in the control, 540, 760, 1120 and 3000 µatm header tanks respectively. Variations were the highest in the 760 µatm tank due to a technical problem that lasted for few hours and during which seawater $pH_T$ dropped to 7.07. $C_T$ concentration increased from 2126±28 to 2454±31 µmol kg⁻¹ between the lowest and highest $\rho CO_2$ treatment. Seawater was, on average, only supersaturated with respect to aragonite in the control treatment while in the other tanks $Ω_{aragonite}$ decreased from 0.9±0.2 to 0.2±0.0 between the 540 and the 3000 µatm treatments. Seawater saturation with respect to calcite was below 1 for the two highest $\rho CO_2$ treatments, with the lowest value for the 3000 µatm tank at 0.33±0.04, and above 1 for the 3 other treatments.

Environmental conditions during the 24 h incubations (n = 5) are presented in Fig. 3. Data represent average values between initial and final sampling for the four cores used per treatment (± SD). For each $\rho CO_2$ treatment, $pH_T$ was relatively constant during the incubations and throughout the experimental period with the highest variation observed for the 540 µatm condition. Higher $A_T$ levels were observed for the 3000 µatm treatment due to the significant release of $A_T$ during the incubations. During the incubations, seawater was undersaturated with respect to aragonite in all cores except the ones from the control treatment. Sediment cores from the control and 540 µatm treatments were exposed to supersaturated conditions with respect to calcite (mean: 2.6±0.2 and 1.3±0.3 respectively) while those from the 1120 and 3000 µatm treatments were maintained below saturation (mean: 0.7±0.1 and 0.3±0.0 respectively) and those from the 760 µatm treatment were close to saturation (mean: 1.0±0.1).

The PCA on the sediment grain size distribution (data not shown) revealed some discrepancies between the different sediment cores used in this study (n = 20). Cores 13 and 14 that were exposed to 1120 µatm had a much larger fraction of medium and fine sand grains than the others. Also, the deepest depth-layer of core 20 exposed to 3000 µatm contained a much larger fraction
<table><caption>**Table 1.** Seawater temperature and carbonate chemistry controlled in the header tanks during the course of the experiment (~14 days; n = 7).</caption><thead><tr><th rowspan="2">Target pCO<sub>2</sub> (µatm)</th><th rowspan="2"></th><th>Tank 1</th><th>Tank 2</th><th>Tank 3</th><th>Tank 4</th><th>Tank 5</th></tr><tr><th>Control</th><th>540</th><th>760</th><th>1120</th><th>3000</th></tr></thead><tbody><tr><td rowspan="4">Water temperature (°C)</td><td>Mean</td><td>0.44</td><td>0.24</td><td>0.13</td><td>0.03</td><td>0.08</td></tr><tr><td>SD</td><td>0.70</td><td>0.78</td><td>0.31</td><td>0.31</td><td>0.27</td></tr><tr><td>Min</td><td>-0.40</td><td>-0.50</td><td>-0.50</td><td>-0.50</td><td>-0.40</td></tr><tr><td>Max</td><td>2.00</td><td>2.10</td><td>0.50</td><td>0.50</td><td>0.50</td></tr><tr><td rowspan="4">pH<sub>T</sub></td><td>Mean</td><td>8.13</td><td>7.84</td><td>7.71</td><td>7.57</td><td>7.17</td></tr><tr><td>SD</td><td>0.08</td><td>0.11</td><td>0.21</td><td>0.08</td><td>0.06</td></tr><tr><td>Min</td><td>8.01</td><td>7.65</td><td>7.07</td><td>7.44</td><td>7.07</td></tr><tr><td>Max</td><td>8.30</td><td>8.14</td><td>7.94</td><td>7.71</td><td>7.26</td></tr><tr><td rowspan="4">C<sub>T</sub> (µmol kg<sup>-1</sup>)</td><td>Mean</td><td>2126.1</td><td>2223.3</td><td>2267.7</td><td>2304.3</td><td>2454.3</td></tr><tr><td>SD</td><td>27.9</td><td>32.9</td><td>74.1</td><td>24.8</td><td>30.7</td></tr><tr><td>Min</td><td>2062.4</td><td>2134.2</td><td>2199.8</td><td>2263.6</td><td>2405.3</td></tr><tr><td>Max</td><td>2165.3</td><td>2288.3</td><td>2515.4</td><td>2345.3</td><td>2510.3</td></tr><tr><td rowspan="4">pCO<sub>2</sub> (µatm)</td><td>Mean</td><td>317.2</td><td>658.5</td><td>1031.2</td><td>1255.8</td><td>3169.4</td></tr><tr><td>SD</td><td>60.3</td><td>162.3</td><td>815.3</td><td>235.3</td><td>437.7</td></tr><tr><td>Min</td><td>200.2</td><td>306.0</td><td>507.2</td><td>874.4</td><td>2539.7</td></tr><tr><td>Max</td><td>422.3</td><td>1024.3</td><td>3959.6</td><td>1674.2</td><td>3981.6</td></tr><tr><td rowspan="4">Ω<sub>aragonite</sub></td><td>Mean</td><td>1.68</td><td>0.93</td><td>0.72</td><td>0.50</td><td>0.21</td></tr><tr><td>SD</td><td>0.24</td><td>0.22</td><td>0.23</td><td>0.09</td><td>0.03</td></tr><tr><td>Min</td><td>1.39</td><td>0.59</td><td>0.16</td><td>0.38</td><td>0.16</td></tr><tr><td>Max</td><td>2.25</td><td>1.66</td><td>1.12</td><td>0.69</td><td>0.25</td></tr><tr><td rowspan="4">Ω<sub>calcite</sub></td><td>Mean</td><td>2.67</td><td>1.48</td><td>1.15</td><td>0.80</td><td>0.33</td></tr><tr><td>SD</td><td>0.38</td><td>0.36</td><td>0.37</td><td>0.15</td><td>0.04</td></tr><tr><td>Min</td><td>2.22</td><td>0.94</td><td>0.25</td><td>0.60</td><td>0.25</td></tr><tr><td>Max</td><td>3.58</td><td>2.65</td><td>1.78</td><td>1.10</td><td>0.39</td></tr></tbody></table>

Salinity was constant throughout the experiment at 34.5. Measured pH in the total scale (pH<sub>T</sub>), dissolved inorganic carbon (C<sub>T</sub>) as well as calculated pCO<sub>2</sub> and saturation states with respect to aragonite (Ω<sub>aragonite</sub>) and calcite (Ω<sub>calcite</sub>), are presented as mean, standard deviation (SD), minimal (Min) and maximal (Max) values.
doi:10.1371/journal.pone.0094068.t001

of coarse sand grains. As all the other cores appeared to share the same grain size characteristics, it has been decided to remove from the subsequent analysis cores 13 and 14 as well as the deepest depth-layer of core 20. Averaged (±SD) sediment grain size distributions in the remaining cores for the different depth layers are shown in Table 2. The sediment cores had a thick coal layer at the bottom.

The correlation circle of the PCA performed on the composition of the sediment with respect to carbon and nitrogen (C<sub>tot</sub> and N<sub>tot</sub>), organic carbon (C<sub>org</sub>), the isotopic composition of C<sub>tot</sub> and C<sub>org</sub> (δ<sup>13</sup>C<sub>tot</sub> and δ<sup>13</sup>C<sub>org</sub>) and organic nitrogen (δ<sup>15</sup>N<sub>org</sub>), as well as the pigments concentrations (chlorophyll a, b and c) is presented in Fig. 4 for the two axis explaining most of the variance in the data (axis 1: 45% and axis 2: 27%). Seawater targeted pCO<sub>2</sub> (µatm) and the depth of sampling (cm) in the sediment have been added as supplementary variables. Data used for this analysis and the complete correlation matrix are presented in Tables S1 and S2, respectively. This analysis showed no significant effect of seawater pCO<sub>2</sub> on any of these parameters and pCO<sub>2</sub> was not correlated to the composite axes 1 and 2. Most of the variance in the dataset was explained by the depth at which parameters were measured in the sediment with significant decreases of N<sub>tot</sub> as well as of Chl a, b and c with depth, because of active diagenesis.

Results from the flux measurements during the 24 h incubations are shown in Fig. 5. At all pCO<sub>2</sub> levels and for all incubations, the sediment was a source of A<sub>T</sub>, C<sub>T</sub> and NO<sub>2</sub><sup>-</sup> and a sink of DO for the overlying water. The sediment was also a source of NH<sub>4</sub><sup>+</sup> and NO<sub>3</sub><sup>-</sup> for the overlying water, except for one incubation at 1120 µatm and at 3000 µatm for NH<sub>4</sub><sup>+</sup> and NO<sub>3</sub><sup>-</sup> respectively. The sediment was most of the time a source of Si except at the start of the experiment at 540 and 760 µatm. DIP fluxes were quite variable with a maximal efflux rate of 114.5 and a maximal influx rate of -14.6 µmol m<sup>-2</sup> h<sup>-1</sup>. The correlation circle of the PCA ran on these fluxes is presented in Fig. 6 for the two axis explaining most of the variance in the data (axis 1: 32% and axis 2: 21%). Seawater temperature and pCO<sub>2</sub> at which sediment cores were exposed during the 24 h incubations (average between initial and final sampling) as well as the time elapsed since the start of the experiment (in days) have been added as supplementary variables. The complete correlation matrix is presented in Table S3. This analysis shows that A<sub>T</sub>, C<sub>T</sub>, NH<sub>4</sub><sup>+</sup> and NO<sub>2</sub><sup>-</sup> fluxes were correlated with pCO<sub>2</sub> (p<0.01). A<sub>T</sub> and C<sub>T</sub> fluxes were positively correlated while NH<sub>4</sub><sup>+</sup> and NO<sub>2</sub><sup>-</sup> fluxes were negatively correlated. The second axis is mainly explained by the variance in Si fluxes that appeared to be significantly correlated with time (p<0.001). Indeed, as shown in Fig. 4, Si fluxes generally increased during the experimental period. Temperature did not have
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Experimental conditions (average between initial and final sampling) at which the sediment cores were exposed during the 24 h incubations (n = 5). :pHr: pH on the total scale, A<sub>T</sub>: total alkalinity, Ω<sub>aragonite</sub> and Ω<sub>calcite</sub>: saturation state of the overlying seawater with respect to, respectively, aragonite and calcite. Error bars represent the standard deviation.
doi:10.1371/journal.pone.0094068.g003</figcaption></figure>

significant effects on any of the observed fluxes. Visual inspections of the evolution of A<sub>T</sub>, C<sub>T</sub>, NH<sub>4</sub><sup>+</sup> and NO<sub>2</sub><sup>-</sup> presented in Fig. 5 suggest that the 3000 µatm treatment is dominant in the observed significant correlations between these fluxes and pCO<sub>2</sub> levels with larger C<sub>T</sub> and A<sub>T</sub> fluxes and relatively lower NH<sub>4</sub><sup>+</sup> and NO<sub>2</sub><sup>-</sup> fluxes compared to the ones measured at lower pCO<sub>2</sub> levels.

Denitrification rates estimated based on <sup>15</sup>N-NO<sub>3</sub><sup>-</sup> labeling are shown in Fig. 7. Increases in <sup>15</sup>N-N<sub>2</sub> were significant and linear over 24 h except for 2 cores from the 540 µatm treatment which, most certainly due to technical problems, did not show any denitrification activities and were therefore excluded from the analysis. Denitrification rates ranged between 11 and 35 µmol N

<table><caption>**Table 2.** Sediment size distribution (%) in the different depth layers, expressed as an average (SD) value between the different cores used in the experiment.</caption><thead><tr><th rowspan="2">Depth layer (cm)</th><th>Silt</th><th>Very fine sand</th><th>Fine sand</th><th>Medium sand</th><th>Coarse sand</th></tr><tr><th>&lt;62.5 µm</th><th>62.5-125 µm</th><th>125-250 µm</th><th>250-500 µm</th><th>500-1000 µm</th></tr></thead><tbody><tr><td>0-0.5</td><td>65.7 (3.0)</td><td>21.3 (1.6)</td><td>8.0 (0.9)</td><td>3.3 (0.5)</td><td>2.0 (0.6)</td></tr><tr><td>0.5-1</td><td>64.4 (2.0)</td><td>23.0 (1.2)</td><td>8.3 (0.9)</td><td>2.9 (0.5)</td><td>1.7 (0.7)</td></tr><tr><td>1-2</td><td>66.5 (2.9)</td><td>22.2 (2.0)</td><td>7.6 (1.1)</td><td>2.5 (0.4)</td><td>1.5 (0.8)</td></tr><tr><td>2-4</td><td>72.9 (3.3)</td><td>19.2 (2.6)</td><td>5.8 (0.8)</td><td>1.6 (0.2)</td><td>0.8 (0.3)</td></tr><tr><td>&gt;4</td><td>77.1 (3.0)</td><td>14.3 (1.5)</td><td>4.2 (0.5)</td><td>1.6 (0.4)</td><td>3.0 (2.0)</td></tr></tbody></table>

Cores 13 and 14 (1120 µatm treatment) as well as the layer deeper than 4 cm of core 20 (3000 µatm treatment) have not been considered in this analysis and in the subsequent ones as they have been showed to greatly differ from the other samples (See text for details).
doi:10.1371/journal.pone.0094068.t002
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Correlation circle from the principal component analysis (PCA) applied on the composition of the sediment with respect to carbon and nitrogen ($C_{tot}$ and $N_{tot}$), to organic carbon ($C_{org}$), to the isotopic composition of $C_{tot}$ and $C_{org}$ ($\delta^{13}C_{tot}$ and $\delta^{13}C_{org}$) and organic nitrogen ($\delta^{15}N_{org}$) as well as to the pigments concentrations (chlorophyll a, b and d). Seawater targeted $pCO_2$ (µatm) and the depth of sampling (cm) in the sediment have been added as supplementary variables. Axis 1 and axis 2 explained, respectively, 45 and 27% of the overall variance in the data. doi:10.1371/journal.pone.0094068.g004</figcaption></figure>

$m^{-2} d^{-1}$ and did not significantly differ (one-way ANOVA, $p=0.1$) between the different $pCO_2$ treatments. During these last incubations, there was no bio-irrigation activity as no significant decrease of bromide has been detected after the 24 h incubation in any of the treatments (data not shown).

## Discussion

Most of the parameters and processes investigated within this study showed no relationship with the overlying seawater pH. Solid-phase concentrations did not change over the 14 day-time period of the experiment. This does not appear as a surprise as changes in sediment solid-phases are usually difficult to detect over short time-scales and effects are stronger reflected in pore-water concentrations and sediment-water exchange rates [49]. The 14-day duration of our experiments may also have limited the impact on slowly reacting pools at depth in the sediments. The increase of silicate fluxes with time (Fig. 5) likely reflects a transient response of biogenic silica dissolution or acidification induced enhanced chemical weathering.

$C_T$ fluxes were the lowest at ~800 µatm and increased below and above this level. This pattern is due to the decrease in the $C_T$ gradient between water and sediment up to 800 µatm and potentially an increase in carbonate dissolution rates above this level. Interestingly, enhanced sediment carbonate dissolution (resulting in an increased $A_T$ release from the sediment) only started at the lowest pH conditions ($pCO_2$ ~3000 µatm) when strong undersaturations with respect to both aragonite and calcite occurred ($\Omega_{aragonite}$ = 0.21, $\Omega_{calcite}$ = 0.33). At 1120 µatm, the overlying water was also corrosive to calcite, but did not seem to induce enhanced dissolution rates. Unfortunately, as no pH, $C_T$ and/or $A_T$ sediment profiles were acquired during this study, it is impossible to investigate the capacity of the sediment to buffer pH changes in pore-waters.

It has to be noted that variations in the overlying water $A_T$ during the 24 h incubations are not only due to carbonate dissolution. Remineralization of particulate organic matter and subsequent nutrient release affect $A_T$ depending on the form of reactive nitrogen produced [50,51]. A release of 1 mole of ammonia or nitrate leads to an increase or decrease, respectively, of $A_T$ by 1 mole. Nitrification is also known to significantly impact $A_T$ (decrease of 2 moles per moles of $NO_3^-$ produced) while nitrification/denitrification coupling has no effect on $A_T$ fluxes. Release of phosphate during remineralization leads to a corresponding decrease of $A_T$. During the present study, $A_T$ release from the sediment was most likely due to oxic mineralization induced carbonate dissolution, as $C_T$ and $DO$ fluxes were by 1-2 orders of magnitude higher than rates of nutrient release/consumption. Since oxygen uptake, and by inference oxic mineralization was not impacted by decreases in pH, the enhanced dissolution rates as observed at the lowest pH levels were most likely directly linked to pH of the overlying water that penetrated to a certain extent in the sediment. It must be stressed that sulfate reduction rates have not been measured during the present study. Sulfate reduction may indeed increase alkalinity in the sediments, but most of the reduced products generated by anaerobic mineralization are re-oxidized leading to a decrease in $A_T$ via proton production, hence balancing the initial $A_T$ production [27,51]. It is very unlikely that an alteration of this process may have been responsible for the increased $A_T$ fluxes observed at low pH. $A_T$ fluxes measured during this study are significantly lower that rates measured for tropical carbonate muds [19] but in the range on what has been measured in carbonate sands [21] although the carbonate content of the sediment in the present study were much lower (21% $CaCO_3$). Values found in this study are in the same range that what has been reported for the Arctic in general (14±8%$CaCO_3$ [52]).

This is, to the best of our knowledge, the first report on the effect of ocean acidification on sediment oxygen fluxes in the Arctic. It must be stressed that, for logistic reasons, our experiment was conducted in darkness, and that the potential effects of increased $CO_2$ on autotrophic processes were not taken into account. During this study, no significant alteration of the oxygen consumption in the overlying water following pH decreases has been shown. Oxygen is used during both organic matter mineralization and re-oxidation processes such as nitrification, the biological oxidation of ammonia into nitrite followed by the oxidation of this nitrite into nitrate [53]. Significant effects of seawater acidification on pelagic nitrification rates have been shown by Huesemann et al. [31], Beman et al. [30] and Kitidis et al. [32], based on both laboratory and field studies. As the $NH_3/NH_4^+$ equilibrium is pH sensitive and in favor of $NH_4^+$ at lower pH, this has been formulated as an explanation for observed decreasing nitrification rates in these studies. Only one study so far focused on the effect of ocean acidification on sediment nitrification rates. In contrast to what has been found for pelagic nitrification, Kitidis et al. [32] showed no effect of a decrease in pH (up to -0.6) on benthic nitrification rates. Based on ammonium, nitrite and oxygen sediment-water fluxes as well as on denitrification rates, it seems unlikely that, in our study, nitrification rates have been impacted by lowered pH as most of the denitrified nitrate came from nitrification in the sediment. A modification of the nitrification activity would have undoubtedly had effects on coupled denitrification rates, which was clearly not the case in our study. The decrease in ammonium and nitrite release, observed only at the lowest experimental pH (i.e. 3000 µatm), could be attributed to an enhanced anammox activity, the anaerobic biological process converting nitrite and
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Sediment-water fluxes of total alkalinity (ΔA<sub>T</sub>), dissolved inorganic carbon (ΔC<sub>T</sub>), ammonium (ΔNH<sub>4</sub><sup>+</sup>), nitrate (ΔNO<sub>3</sub><sup>-</sup>), nitrite (ΔNO<sub>2</sub><sup>-</sup>), dissolved inorganic phosphorus (ΔDIP) and silicate (ΔSi) observed during the 24 h incubations (n = 5). Positive rates represent fluxes from the sediment to the overlying water.
doi:10.1371/journal.pone.0094068.g005</figcaption></figure>

ammonium directly into dinitrogen gas [54]. In a companion
paper [39], an increase in the abundance of Planctonomycete-
specific 16S rRNA, indicative of anammox bacteria, has been
shown and several studies already suggested that this process will
be favored in a higher CO₂ ocean [34]. Interestingly, Tait et al.
[39] also showed that effect of increasing pCO₂ was significant in
terms of bacterial community composition at the intermediate
levels of 540 and 760 µatm, with reduced abundance of bacterial
amoA transcript abundance (aerobic oxidizers) at 760 µatm and
above. However, the abundance of archaeal amoA transcripts did
not change within these treatments and actually increased at
3000 µatm, suggesting that bacterial ammonia oxidisers may be
more vulnerable to increased CO₂ than archaeal ammonia
oxidisers. The lack of response of nitrogen flux measurements at

the intermediate CO₂ levels may have been due to continued
archaeal amoA activity.

Previous studies focusing on the effects of ocean acidification on sediment-water fluxes reported stronger impacts than the ones presented here [33–37]. First of all, these experiments considered usually much lower pH levels than the ones projected for the end of the century. Moreover, the observed impacts were, in most cases, indirect and attributed to significant effects of ocean acidification on bioturbating and bio-irrigating organisms [33–37]. As mentioned previously, the sediment used in the present study was exempt of these organisms and transport processes were dominated by molecular diffusion. Although enhanced dissolution rates at the lowest pH levels showed that pH of overlying water had an impact on the sediment, the absence of sediment pH microprofiles during this study does not allow assessing the
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Correlation circle from the principal component analysis (PCA) applied on sediment-water fluxes of total alkalinity (ΔT), dissolved inorganic carbon (Ct), ammonium (NH4+), nitrate (NO3-), nitrite (NO2-), dissolved inorganic phosphorus (DIP) and silicate (Si). Seawater temperature (°C) and pCO2 (µatm) at which sediment cores were exposed during the incubations (average between initial and final sampling) as well as the time elapsed since the start of the experiment (in days) have been added as supplementary variables. Axis 1 and axis 2 explained, respectively, 32 and 21% of the overall variance in the data.
doi:10.1371/journal.pone.0094068.g006</figcaption></figure>

penetration of ocean acidification perturbation, after 14 days of exposure, into the sediment.

The absence of bioturbating organisms is a common feature of the sediments in the inner fjord where samples were taken from [55], but large abundances of bioturbating organisms can be found in the central and outer fjord. It was initially planned to sample the sediment cores in the central fjord but a thick ice cover blocked access to the majority of the fjord adjacent to Ny-Ålesund and made collection of samples and sediment cores extremely difficult. Net sediment-water fluxes of nitrogen and oxygen have been measured from sediments sampled in the same fjord but in deeper areas [56]. While oxygen uptakes were 6 times lower than the ones measured in the present study (4 vs. ~24 mmol m⁻² d⁻¹), their sediment acted as a sink of NH₄⁺ for the overlying water, as opposed to the constant net efflux observed in our study. In contrast, NO₃⁻ efflux appeared very similar between the two studies with values around 300 µmol m⁻² d⁻¹. As already mentioned, the sediments used in the present study were much shallower. As such, extrapolations based on our results should be done with care.

Based on the present results, we can conclude that ocean acidification, under the conditions of our experiment, will have limited impacts on sediment carbonate dissolution, on the microbial activity (e.g. oxic and anoxic mineralization, nitrification) and associated sediment-water fluxes on Arctic shelves. The potential effects that have been shown regarding carbonate dissolution and the anammox activity occurred at the lowest experimental pH (~3000 µatm), a level that is not foreseen within the next 100 years. However, as discussed previously, the sediment used in this study was exempt of bio-irrigating organisms and the penetration of low pH water inside the sediment might have been limited following an exposure time of 14 days. Longer-term experiments on different sediment types and under different light conditions (different seasons) are still required to confirm the

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Denitrification rates measured over the 24 h incubations performed at the end of the experimental period for the different pCO₂ levels considered in the present study.
doi:10.1371/journal.pone.0094068.g007</figcaption></figure>

limited impact of ocean acidification on shallow Arctic sediment processes as observed in this study.

## Supporting Information

<caption>**Table S1** Composition of all sediment cores (at several depths, in cm) after a 14 days exposure at various partial pressure of CO₂ levels (pCO₂ in µatm, 4 cores per treatment). The content in carbon and nitrogen (C<sub>tot</sub> and N<sub>tot</sub> in %), organic carbon (C<sub>org</sub> in %), the isotopic composition of C<sub>tot</sub> and C<sub>org</sub> (δ¹³C<sub>tot</sub> and δ¹³C<sub>org</sub> in ‰) and organic nitrogen (δ¹⁵N<sub>org</sub> in ‰), as well as the pigment concentrations (chlorophyll *a*, *b* and *c* in µg g⁻¹) are presented.
(DOCX)</caption>

<caption>**Table S2** Correlation matrix associated with the PCA performed on the composition of the sediment at the end of the experiment (14 days) with respect to carbon and nitrogen (C<sub>tot</sub> and N<sub>tot</sub>), organic carbon (C<sub>org</sub>), the isotopic composition of C<sub>tot</sub> and C<sub>org</sub> (δ¹³C<sub>tot</sub> and δ¹³C<sub>org</sub>) and organic nitrogen (δ¹⁵N<sub>org</sub>), as well as the pigments concentrations (chlorophyll *a*, *b* and *c*). Seawater targeted partial pressure of CO₂ (pCO₂ in µatm) and the depth of sampling (cm) in the sediment have been added as supplementary variables. The first two axes of the PCA represented 45 and 27% of the variance, respectively. All correlation values were tested for significance (Pearson correlation test performed on R, package Psych). Values in bold are associated with a *p* value below 0.01.
(DOCX)</caption>

<caption>**Table S3** Correlation matrix associated with the PCA performed on the sediment-water fluxes of total alkalinity (ΔΔT), dissolved inorganic carbon (ΔC<sub>T</sub>), dissolved oxygen (ΔDO), ammonium (ΔNH₄⁺), nitrite (ΔNO₂⁻), nitrate (ΔNO₃⁻), dissolved inorganic phosphate (ΔDIP) and silicate (ΔSi). Seawater temperature (Temp) and partial pressure of CO₂ (pCO₂) at which sediment cores were exposed during the 24 h incubations (average between initial and final sampling) as well as the time elapsed since the start of the experiment (in days) have been added as supplementary variables. The first two axes of the PCA represented 32 and 21% of the variance, respectively. All correlation values were tested for significance (Pearson correlation test performed using the R software, package Psych). Values in bold are associated with a *p* value below 0.01.
(DOCX)</caption>
Acknowledgments

We thank H. Findlay and H. Wood for setting up the mesocosm system, M. -D. Pizay and J. -P. Gattuso for providing carbonate data, M. Schwanitz, S. Treydte and J. Mönnich for collecting the sediment cores as well as three reviewers for their constructive comments.

References

1. Sabine CL, Feely RA, Gruber N, Key RM, Lee K, et al. (2004) The oceanic sink for anthropogenic CO₂. Science 305: 367–371.

2. Gattuso J-P, Hansson L (2011) Ocean acidification: background and history. In: J.-P. Gattuso and L. Hansson, editors. Ocean acidification. Oxford: Oxford University Press. pp. 1–20.

3. Orr JC (2011) Recent and future changes in ocean carbonate chemistry. In: J.-P. Gattuso and L. Hansson, editors. Oxford: Oxford University Press. pp. 41–66.

4. Pörtner HO, Langenbuch M, Reipschlager A (2004) Biological impact of elevated ocean CO₂ concentrations: Lessons from animal physiology and earth history. J Ocean 60: 705–718.

5. Riebesell U, Schulz KG, Bellerby RGJ, Botros M, Fritsche P, et al. (2007) Enhanced biological carbon consumption in a high CO₂ ocean. Nature 450: 545–U510.

6. Doney SC, Fabry VJ, Feely RA, Kleypas JA (2009) Ocean acidification: the other CO₂ problem. Ann Rev Mar Sci 1: 169–192.

7. Orr JC, Fabry VJ, Aumont O, Bopp L, Doney SC, et al. (2005) Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms. Nature 437: 681–686.

8. Olafsson J, Olafsdottir SR, Benoit-Cattin A, Danielsen M, Arnarson TS, et al. (2009) Rate of Iceland Sea acidification from time series measurements. Biogeosciences 6: 2661–2668.

9. Yamamoto-Kawai M, McLaughlin FA, Carmack EC, Nishino S, Shimada K (2009) Aragonite undersaturation in the Arctic Ocean: Effects of ocean acidification and sea ice melt. Science 326: 1098–1100.

10. Bates NR, Mathis JT, Cooper LW (2009) Ocean acidification and biologically induced seasonality of carbonate mineral saturation states in the western Arctic Ocean. J Geophys Res-Oceans 114.

11. Steinacher M, Joos F, Frölicher TL, Plattner GK, Doney SC (2009) Imminent ocean acidification in the Arctic projected with the NCAR global coupled carbon cycle-climate model. Biogeosciences 6: 515–533.

12. Chierici M, Fransson A (2009) Calcium carbonate saturation in the surface water of the Arctic Ocean: undersaturation in freshwater influenced shelves. Biogeosciences 6: 2421–2431.

13. Sigman DM, Boyle EA (2000) Glacial/interglacial variations in atmospheric carbon dioxide. Nature 407: 859–869.

14. Andersson A, Mackenzie FT, Ver LM (2003) Solution of shallow-water carbonates: an insignificant buffer against rising atmospheric CO₂. Geology 31: 513–516.

15. Feely RA, Sabine CL, Lee K, Berelson W, Kleypas J, et al. (2004) Impact of anthropogenic CO₂ on the CaCO₃ system in the oceans. Science 305: 362–366.

16. Caldeira K, Archer D, Barry JP, Bellerby RGJ, Brewer PG, et al. (2007) Comment on “Modern-age buildup of CO₂ and its effects on seawater acidity and salinity” by Hugo A. Loaiciga. Geophys Res Lett 34: 3.

17. Andersson AJ, Yeakel KL, Bates NR, de Patron SJ (2014) Partial offsets in ocean acidification from changing coral reef biogeochemistry. Nat Clim Chang 4: 56–61.

18. Alongi DM, Pfitzner J, Trott LA (2006) Deposition and cycling of carbon and nitrogen in carbonate mud of the lagoons of Arlington and Sudbury Reefs, Great Barrier Reef. Coral Reefs 25: 123–143.

19. Yates KK, Halley RB (2006) Diurnal variation in rates of calcification and carbonate sediment dissolution in Florida Bay. Estuaries Coast 29: 24–39.

20. Barnes DJ, Cuff C (2000) Solution of reef rock buffers seawater against rising atmospheric CO₂. Proceedings of the Ninth International Coral Reef Symposium, Indonesia: State Ministry for the Environment.pp. 248p.

21. Cyronak T, Santos IR, McMahon A, Eyre BD (2013) Carbon cycling hysteresis in permeable carbonate sands over a diel cycle: Implications for ocean acidification. Limnol Oceanogr 58: 131–143.

22. Rao AMF, Polerecky L, Ionescu D, Meysman FJR, de Beer D (2012) The influence of pore-water advection, benthic photosynthesis, and respiration on calcium carbonate dynamics in reef sands. Limnol Oceanogr 57: 809–825.

23. Hauck J, Arrigo KR, Hoppema M, van Dijken GL, Volker C, et al. (2013) Insignificant buffering capacity of Antarctic shelf carbonates. Global Biogeochem Cycles 27: 11–20.

24. Burdige DJ (2006) Geochemistry of marine sediments. 209 p.

25. Gattuso J-P, Gentili B, Kleypas JA, Middelburg JJ, Antoine D (2006) Light availability in the coastal ocean: impact on the distribution of benthic photosynthetic organisms and their contribution to primary production. Biogeosciences 3: 489–513.

26. Middelburg JJ (2011) Chemoautotrophy in the ocean. Geophys Res Lett 38.

27. Hu XP, Cai WJ (2011) An assessment of ocean margin anaerobic processes on oceanic alkalinity budget. Global Biogeochem Cycles 25.

28. Thomas H, Schietteccate LS, Suygens K, Kone YJM, Shadwick EH, et al. (2009) Enhanced ocean carbon storage from anaerobic alkalinity generation in coastal sediments. Biogeosciences 6: 267–274.

Author Contributions

Conceived and designed the experiments: FG JJM. Performed the experiments: FG PvR LP. Analyzed the data: FG LP JJM. Contributed reagents/materials/analysis tools: PvR JJM. Wrote the paper: FG LP JJM.

29. Krumins V, Gehlen M, Arndt S, Van Cappellen P, Regnier P (2013) Dissolved inorganic carbon and alkalinity fluxes from coastal marine sediments: model estimates for different shelf environments and sensitivity to global change. Biogeosciences 10: 371–398.

30. Beman JM, Chow CE, King AL, Feng YY, Fuhrman JA, et al. (2011) Global declines in oceanic nitrification rates as a consequence of ocean acidification. Proc Natl Acad Sci U S A 108: 208–213.

31. Huesemann MH, Skillman AD, Crecelius EA (2002) The inhibition of marine nitrification by ocean disposal of carbon dioxide. Mar Pollut Bull 44: 142–148.

32. Kitidis V, Laverock B, McNeill LC, Beesley A, Cummings D, et al. (2011) Impact of ocean acidification on benthic and water column ammonia oxidation. Geophys Res Lett 38.

33. Wood HL, Widdicombe S, Spicer JI (2009) The influence of hypercapnia and the infenal brittlestar Amphiprura filiformis on sediment nutrient flux - will ocean acidification affect nutrient exchange? Biogeosciences 6: 2015–2024.

34. Widdicombe S, Needham HR (2007) Impact of CO₂ induced seawater acidification on the burrowing activity of Nereis virens and sediment nutrient flux. Mar Ecol Prog Ser 341: 111–122.

35. Widdicombe S, Dashfield SL, McNeill CL, Needham HR, Beesley A, et al. (2009) Effects of CO₂ induced seawater acidification on infenal diversity and sediment nutrient fluxes. Mar Ecol Prog Ser 379: 59–75.

36. Dashfield SL, Sommerfield PJ, Widdicombe S, Austen MC, Nimmo M (2008) Impacts of ocean acidification and burrowing urchins on within-sediment pH profiles and subtidal nematode communities. J Exp Mar Bio Ecol 365: 46–52.

37. Donohue PJC, Calosi P, Bates AH, Laverock B, Rastrick S, et al. (2012) Impact of exposure to elevated pCO₂ on the physiology and behaviour of an important ecosystem engineer, the burrowing shrimp *Upogebia deltairea*. Aquat Biol 15: 73–86.

38. Meysman FJR, Galaktionov OS, Gribsholt B, Middelburg JJ (2006) Bio-irrigation in permeable sediments: Advective pore water transport induced by burrow ventilation Limnol Oceanogr 51: 142–156.

39. Tait K, Laverock B, Widdicombe S (2013) Response of an Arctic sediment nitrogen cycling community to increased CO₂. Estuaries Coast: 1–12.

40. Barry JP, Tyrrell T, Hansson L, Gattuso J-P, Plattner G-K (2010) Section 2.1: Atmospheric CO₂ Targets for Ocean Acidification Perturbation Experiments. In: U. Riebesell, V. J. Fabry, L. Hansson and J.-P. Gattuso, editors. Guide for Best Practices in Ocean Acidification Research and Data Reporting.

41. Dickson AG, Sabine CL, Christian JR (2007) Guide to best practices for ocean CO₂ measurements. PICES Special Publication 3, 191 pp.

42. Van den Meersche K, Middelburg JJ, Soetaert K, Rijswijk Pv, Boschker HTS, et al. (2004) Carbon-Nitrogen coupling and algal bacterial interactions during an experimental bloom: modelling a ¹³C tracer experiment. Limnol Oceanogr 49: 862–878.

43. Knap AH, Michaels AE, Close A, Ducklow HW, Dickson AG (1996) Protocols for the Joint Global Ocean Flux Study (JGOFS) core measurements. JGOFS Report No 19. Bergen, Norway: Unesco. pp. vi+170.

44. Gribsholt B, Boschker HTS, Struyf E, Andersson M, Tramper A, et al. (2005) Nitrogen processing in a tidal freshwater marsh: A whole-ecosystem ¹⁵N labeling study. Limnol Oceanogr 50: 1945–1959.

45. Weiss RF, Price BA (1980) Nitrous oxide solubility in water and seawater. Mar Chem 8: 347–359.

46. Forster S, Khalili A, Kitlar J (2003) Variation of nonlocal irrigation in a subtidal benthic community. J Mar Res 61: 335–357.

47. Nieuwenhuize J, Maas YEM, Middelburg JJ (1994) Rapid analysis of organic carbon and nitrogen in particulate materials. Mar Chem 45: 217–224.

48. Herman PMJ, Middelburg JJ, Widdows J, Lucas CH, Heip CHR (2000) Stable isotopes' as trophic tracers: combining field sampling and manipulative labelling of food resources for macrobenthos. Mar Ecol Prog Ser 204: 79–92.

49. Berner RA (1980) Early diagenesis: a theoretical approach. Princeton, NJ: Princeton University Press. 241 p.

50. Wolf-Gladrow DA, Zeebe RE, Klaas C, Körtzinger A, Dickson AG (2007) Total alkalinity: the explicit conservative expression and its application to biogeochemical processes. Mar Chem 106: 287–300.

51. Soetaert K, Hofmann AF, Middelburg JJ, Meysman FJR, Greenwood J (2007) The effect of biogeochemical processes on pH. Mar Chem 105: 30–51.

52. Darby DA, Naidu AS, Mowatt TC, Jones G (1990) Sediment composition and sedimentary processes in the Arctic Ocean. In: Y. Herman, editor. The Arctic Seas. Springer US. pp. 657–720.

53. Middelburg JJ, Levin LA (2009) Coastal hypoxia and sediment biogeochemistry. Biogeosciences 6: 1273–1293.

54. den Camp H, Kartal B, Guven D, van Niftrik L, Haaijer SCM, et al. (2006) Global impact and application of the anaerobic ammonium-oxidizing (anammox) bacteria. Biochem Soc Trans 34: 174–178.
55. Zaborska A, Pempkowiak J, Papucci C (2006) Some sediment characteristics and sedimentation rates in an Arctic Fjord (Kongsfjorden, Svalbard). *Annual Environmental Protection* 8: 79–96.

56. Gihring TM, Lavik G, Kuypers MMM, Kostka JE (2010) Direct determination of nitrogen cycling rates and pathways in Arctic fjord sediments (Svalbard, Norway). *Limnol Oceanogr* 55: 740–752.