Consistent increase in dimethyl sulfide (DMS) in response to high CO2 in five shipboard bioassays from contrasting NW European waters

The ubiquitous marine trace gas dimethyl sulfide (DMS) comprises the greatest natural source of sulfur to the atmosphere and is a key player in atmospheric chemistry and climate. We explore the short-term response of DMS production and cycling and that of its algal precursor dimethyl sulfoniopropionate (DMSP) to elevated carbon dioxide (CO₂) and ocean acidification (OA) in five 96 h shipboard bioassay experiments. Experiments were performed in June and July 2011, using water collected from contrasting sites in NW European waters (Outer Hebrides, Irish Sea, Bay of Biscay, North Sea). Concentrations of DMS and DMSP, alongside rates of DMSP synthesis and DMS production and consumption, were determined during all experiments for ambient CO₂ and three high-CO₂ treatments (550, 750, 1000 μatm). In general, the response to OA throughout this region showed little variation, despite encompassing a range of biological and biogeochemical conditions. We observed consistent and marked increases in DMS concentrations relative to ambient controls (110% (28–223%) at 550 μatm, 153% (56–295%) at 750 μatm and 225% (79–413%) at 1000 μatm), and decreases in DMSP concentrations (28% (18–40%) at 550 μatm, 44% (18–64%) at 750 μatm and 52% (24–72%) at 1000 μatm). Significant decreases in DMSP synthesis rate constants (μDMSP, d⁻¹) and DMSP production rates (nmol d⁻¹) were observed in two experiments (7–90% decrease), whilst the response under high CO₂ from the remaining experiments was generally indistinguishable from ambient controls. Rates of bacterial DMS gross consumption and production gave weak and inconsistent responses to high CO₂. The variables and rates we report increase our understanding of the processes behind the response to OA. This could provide the opportunity to improve upon mesocosm-derived empirical modelling relationships and to move towards a mechanistic approach for predicting future DMS concentrations.

The influence of ocean acidification on nitrogen regeneration and nitrous oxide production in the northwest European shelf sea

The assimilation and regeneration of dissolved inorganic nitrogen, and the concentration of N₂O, was investigated at stations located in the NW European shelf sea during June/July 2011. These observational measurements within the photic zone demonstrated the simultaneous regeneration and assimilation of NH₄⁺, NO₂⁻ and NO₃⁻. NH₄⁺ was assimilated at 1.82–49.12 nmol N L⁻¹ h⁻¹ and regenerated at 3.46–14.60 nmol N L⁻¹ h⁻¹; NO₂^- was assimilated at 0–2.08 nmol N L⁻¹ h⁻¹ and regenerated at 0.01–1.85 nmol N L⁻¹ h⁻¹; NO₃⁻ was assimilated at 0.67–18.75 nmol N L⁻¹ h⁻¹ and regenerated at 0.05–28.97 nmol N L⁻¹ h⁻¹. Observations implied that these processes were closely coupled at the regional scale and that nitrogen recycling played an important role in sustaining phytoplankton growth during the summer. The [N₂O], measured in water column profiles, was 10.13 ± 1.11 nmol L⁻¹ and did not strongly diverge from atmospheric equilibrium indicating that sampled marine regions were neither a strong source nor sink of N₂O to the atmosphere. Multivariate analysis of data describing water column biogeochemistry and its links to N-cycling activity failed to explain the observed variance in rates of N-regeneration and N-assimilation, possibly due to the limited number of process rate observations. In the surface waters of five further stations, ocean acidification (OA) bioassay experiments were conducted to investigate the response of NH₄⁺ oxidising and regenerating organisms to simulated OA conditions, including the implications for [N₂O]. Multivariate analysis was undertaken which considered the complete bioassay data set of measured variables describing changes in N-regeneration rate, [N₂O] and the biogeochemical composition of seawater. While anticipating biogeochemical differences between locations, we aimed to test the hypothesis that the underlying mechanism through which pelagic N-regeneration responded to simulated OA conditions was independent of location. Our objective was to develop a mechanistic understanding of how NH₄⁺ regeneration, NH₄⁺ oxidation and N₂O production responded to OA. Results indicated that N-regeneration process responses to OA treatments were location specific; no mechanistic understanding of how N-regeneration processes respond to OA in the surface ocean of the NW European shelf sea could be developed.

Free-ocean CO2 enrichment (FOCE) systems: present status and future developments

Free-ocean CO2 enrichment (FOCE) systems are designed to assess the impact of ocean acidification on biological communities in situ for extended periods of time (weeks to months). They overcome some of the drawbacks of laboratory experiments and field observations by enabling (1) precise control of CO2 enrichment by monitoring pH as an offset of ambient pH, (2) consideration of indirect effects such as those mediated through interspecific relationships and food webs, and (3) relatively long experiments with intact communities. Bringing perturbation experiments from the laboratory to the field is, however, extremely challenging. The main goal of this paper is to provide guidelines on the general design, engineering, and sensor options required to conduct FOCE experiments. Another goal is to introduce xFOCE, a community-led initiative to promote awareness, provide resources for in situ perturbation experiments, and build a user community. Present and existing FOCE systems are briefly described and examples of data collected presented. Future developments are also addressed as it is anticipated that the next generation of FOCE systems will include, in addition to pH, options for oxygen and/or temperature control. FOCE systems should become an important experimental approach for projecting the future response of marine ecosystems to environmental change.

Flexible C : N ratio enhances metabolism of large phytoplankton when resource supply is intermittent

Phytoplankton cell size influences particle sinking rate, food web interactions and biogeographical distributions. We present a model in which the uptake, storage and assimilation of nitrogen and carbon are explicitly resolved in different-sized phytoplankton cells. In the model, metabolism and cellular C :N ratio are influenced by the accumulation of carbon polymers such as carbohydrate and lipid, which is greatest when cells are nutrient starved, or exposed to high light. Allometric relations and empirical data sets are used to constrain the range of possible C : N, and indicate that larger cells can accumulate significantly more carbon storage compounds than smaller cells. When forced with extended periods of darkness combined with brief exposure to saturating irradiance, the model predicts organisms large enough to accumulate significant carbon reserves may on average synthesize protein and other functional apparatus up to five times faster than smaller organisms. The advantage of storage in terms of average daily protein synthesis rate is greatest when modeled organisms were previously nutrient starved, and carbon storage reservoirs saturated. Small organisms may therefore be at a disadvantage in terms of average daily growth rate in environments that involve prolonged periods of darkness and intermittent nutrient limitation. We suggest this mechanism is a significant constraint on phytoplankton C :N variability and cell size distribution in different oceanic regimes.

Dynamics of air–sea CO2 fluxes in the northwestern European shelf based on voluntary observing ship and satellite observations

From January 2011 to December 2013, we constructed a comprehensive pCO2 data set based on voluntary observing ship (VOS) measurements in the western English Channel (WEC). We subsequently estimated surface pCO2 and air–sea CO2 fluxes in northwestern European continental shelf waters using multiple linear regressions (MLRs) from remotely sensed sea surface temperature (SST), chlorophyll a concentration (Chl a), wind speed (WND), photosynthetically active radiation (PAR) and modeled mixed layer depth (MLD). We developed specific MLRs for the seasonally stratified northern WEC (nWEC) and the permanently well-mixed southern WEC (sWEC) and calculated surface pCO2 with uncertainties of 17 and 16 μatm, respectively. We extrapolated the relationships obtained for the WEC based on the 2011–2013 data set (1) temporally over a decade and (2) spatially in the adjacent Celtic and Irish seas (CS and IS), two regions which exhibit hydrographical and biogeochemical characteristics similar to those of WEC waters. We validated these extrapolations with pCO2 data from the SOCAT and LDEO databases and obtained good agreement between modeled and observed data. On an annual scale, seasonally stratified systems acted as a sink of CO2 from the atmosphere of −0.6 ± 0.3, −0.9 ± 0.3 and −0.5 ± 0.3 mol C m−2 yr−1 in the northern Celtic Sea, southern Celtic sea and nWEC, respectively, whereas permanently well-mixed systems acted as source of CO2 to the atmosphere of 0.2 ± 0.2 and 0.3 ± 0.2 mol C m−2 yr−1 in the sWEC and IS, respectively. Air–sea CO2 fluxes showed important inter-annual variability resulting in significant differences in the intensity and/or direction of annual fluxes. We scaled the mean annual fluxes over these provinces for the last decade and obtained the first annual average uptake of −1.11 ± 0.32 Tg C yr−1 for this part of the northwestern European continental shelf. Our study showed that combining VOS data with satellite observations can be a powerful tool to estimate and extrapolate air–sea CO2 fluxes in sparsely sampled area.

The biogeochemical structuring role of horizontal stirring: Lagrangian perspectives on iron delivery downstream of the Kerguelen Plateau

Field campaigns are instrumental in providing ground truth for understanding and modeling global ocean biogeochemical budgets. A survey however can only inspect a fraction of the global oceans, typically a region hundreds of kilometers wide for a temporal window of the order of (at most) several weeks. This spatiotemporal domain is also the one in which the mesoscale activity induces through horizontal stirring a strong variability in the biogeochemical tracers, with ephemeral, local contrasts which can easily mask the regional and seasonal gradients. Therefore, whenever local in situ measures are used to infer larger-scale budgets, one faces the challenge of identifying the mesoscale structuring effect, if not simply to filter it out. In the case of the KEOPS2 investigation of biogeochemical responses to natural iron fertilization, this problem was tackled by designing an adaptive sampling strategy based on regionally optimized multisatellite products analyzed in real time by specifically designed Lagrangian diagnostics. This strategy identified the different mesoscale and stirring structures present in the region and tracked the dynamical frontiers among them. It also enabled back trajectories for the ship-sampled stations to be estimated, providing important insights into the timing and pathways of iron supply, which were explored further using a model based on first-order iron removal. This context was essential for the interpretation of the field results. The mesoscale circulation-based strategy was also validated post-cruise by comparing the Lagrangian maps derived from satellites with the patterns of more than one hundred drifters, including some adaptively released during KEOPS2 and a subsequent research voyage. The KEOPS2 strategy was adapted to the specific biogeochemical characteristics of the region, but its principles are general and will be useful for future in situ biogeochemical surveys.

Effect of ocean acidification and elevated fCO2 on trace gas production by a Baltic Sea summer phytoplankton community

The Baltic Sea is a unique environment as the largest body of brackish water in the world. Acidification of the surface oceans due to absorption of anthropogenic CO2 emissions is an additional stressor facing the pelagic community of the already challenging Baltic Sea. To investigate its impact on trace gas biogeochemistry, a large-scale mesocosm experiment was performed off Tvärminne Research Station, Finland in summer 2012. During the second half of the experiment, dimethylsulphide (DMS) concentrations in the highest fCO2 mesocosms (1075–1333 μatm) were 34 % lower than at ambient CO2 (350 μatm). However the net production (as measured by concentration change) of seven halocarbons analysed was not significantly affected by even the highest CO2 levels after 5 weeks exposure. Methyl iodide (CH3I) and diiodomethane (CH2I2) showed 15 % and 57 % increases in mean mesocosm concentration (3.8 ± 0.6 pmol L−1 increasing to 4.3 ± 0.4 pmol L−1 and 87.4 ± 14.9 pmol L−1 increasing to 134.4 ± 24.1 pmol L−1 respectively) during Phase II of the experiment, which were unrelated to CO2 and corresponded to 30 % lower Chl-ɑ concentrations compared to Phase I. No other iodocarbons increased or showed a peak, with mean chloroiodomethane (CH2ClI) concentrations measured at 5.3 (± 0.9) pmol L−1 and iodoethane (C2H5I) at 0.5 (± 0.1) pmol L−1. Of the concentrations of bromoform (CHBr3; mean 88.1 ± 13.2 pmol L−1), dibromomethane (CH2Br2; mean 5.3 ± 0.8 pmol L−1) and dibromochloromethane (CHBr2Cl, mean 3.0 ± 0.5 pmol L−1), only CH2Br2 showed a decrease of 17 % between Phases I and II, with CHBr3 and CHBr2Cl showing similar mean concentrations in both Phases. Outside the mesocosms, an upwelling event was responsible for bringing colder, high CO2, low pH water to the surface starting on day t16 of the experiment; this variable CO2 system with frequent upwelling events implies the community of the Baltic Sea is acclimated to regular significant declines in pH caused by up to 800 μatm fCO2. After this upwelling, DMS concentrations declined, but halocarbon concentrations remained similar or increased compared to measurements prior to the change in conditions. Based on our findings, with future acidification of Baltic Sea waters, biogenic halocarbon emissions are likely to remain at similar values to today, however emissions of biogenic sulphur could significantly decrease from this region.

iMarNet: an ocean biogeochemistry model intercomparison project within a common physical ocean modelling framework

Ocean biogeochemistry (OBGC) models span a wide variety of complexities, including highly simplified nutrient-restoring schemes, nutrient–phytoplankton–zooplankton–detritus (NPZD) models that crudely represent the marine biota, models that represent a broader trophic structure by grouping organisms as plankton functional types (PFTs) based on their biogeochemical role (dynamic green ocean models) and ecosystem models that group organisms by ecological function and trait. OBGC models are now integral components of Earth system models (ESMs), but they compete for computing resources with higher resolution dynamical setups and with other components such as atmospheric chemistry and terrestrial vegetation schemes. As such, the choice of OBGC in ESMs needs to balance model complexity and realism alongside relative computing cost. Here we present an intercomparison of six OBGC models that were candidates for implementation within the next UK Earth system model (UKESM1). The models cover a large range of biological complexity (from 7 to 57 tracers) but all include representations of at least the nitrogen, carbon, alkalinity and oxygen cycles. Each OBGC model was coupled to the ocean general circulation model Nucleus for European Modelling of the Ocean (NEMO) and results from physically identical hindcast simulations were compared. Model skill was evaluated for biogeochemical metrics of global-scale bulk properties using conventional statistical techniques. The computing cost of each model was also measured in standardised tests run at two resource levels. No model is shown to consistently outperform all other models across all metrics. Nonetheless, the simpler models are broadly closer to observations across a number of fields and thus offer a high-efficiency option for ESMs that prioritise high-resolution climate dynamics. However, simpler models provide limited insight into more complex marine biogeochemical processes and ecosystem pathways, and a parallel approach of low-resolution climate dynamics and high-complexity biogeochemistry is desirable in order to provide additional insights into biogeochemistry–climate interactions.