Bridging the gap between marine biogeochemical and fisheries sciences; configuring the zooplankton link

Exploring climate and anthropogenic impacts on marine ecosystems requires an understanding of how trophic components interact. However, integrative end-to-end ecosystem studies (experimental and/or modelling) are rare. Experimental investigations often concentrate on a particular group or individual species within a trophic level, while tropho-dynamic field studies typically employ either a bottom-up approach concentrating on the phytoplankton community or a top-down approach concentrating on the fish community. Likewise the emphasis within modelling studies is usually placed upon phytoplankton-dominated biogeochemistry or on aspects of fisheries regulation. In consequence the roles of zooplankton communities (protists and metazoans) linking phytoplankton and fish communities are typically under-represented if not (especially in fisheries models) ignored. Where represented in ecosystem models, zooplankton are usually incorporated in an extremely simplistic fashion, using empirical descriptions merging various interacting physiological functions governing zooplankton growth and development, and thence ignoring physiological feedback mechanisms. Here we demonstrate, within a modelled plankton food-web system, how trophic dynamics are sensitive to small changes in parameter values describing zooplankton vital rates and thus the importance of using appropriate zooplankton descriptors. Through a comprehensive review, we reveal the mismatch between empirical understanding and modelling activities identifying important issues that warrant further experimental and modelling investigation. These include: food selectivity, kinetics of prey consumption and interactions with assimilation and growth, form of voided material, mortality rates at different age-stages relative to prior nutrient history. In particular there is a need for dynamic data series in which predator and prey of known nutrient history are studied interacting under varied pH and temperature regimes.

Have we been underestimating the effects of ocean acidification in zooplankton?

Understanding how copepods may respond to ocean acidification (OA) is critical for risk assessments of ocean ecology and biogeochemistry. The perception that copepods are insensitive to OA is largely based on experiments with adult females. Their apparent resilience to increased carbon dioxide (pCO2 ) concentrations has supported the view that copepods are 'winners' under OA. Here, we show that this conclusion is not robust, that sensitivity across different life stages is significantly misrepresented by studies solely using adult females. Stage-specific responses to pCO2 (385-6000 μatm) were studied across different life stages of a calanoid copepod, monitoring for lethal and sublethal responses. Mortality rates varied significantly across the different life stages, with nauplii showing the highest lethal effects; nauplii mortality rates increased threefold when pCO2 concentrations reached 1000 μatm (year 2100 scenario) with LC50 at 1084 μatm pCO2 . In comparison, eggs, early copepodite stages, and adult males and females were not affected lethally until pCO2 concentrations ≥3000 μatm. Adverse effects on reproduction were found, with >35% decline in nauplii recruitment at 1000 μatm pCO2 . This suppression of reproductive scope, coupled with the decreased survival of early stage progeny at this pCO2 concentration, has clear potential to damage population growth dynamics in this species. The disparity in responses seen across the different developmental stages emphasizes the need for a holistic life-cycle approach to make species-level projections to climate change. Significant misrepresentation and error propagation can develop from studies which attempt to project outcomes to future OA conditions solely based on single life history stage exposures.

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

The assimilation and regeneration of dissolved inorganic nitrogen, and the concentration of N2O, 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 NH4+, NO2− and NO3−. NH4+ was assimilated at 1.82–49.12 nmol N L−1 h−1 and regenerated at 3.46–14.60 nmol N L−1 h−1; NO2− was assimilated at 0–2.08 nmol N L−1 h−1 and regenerated at 0.01–1.85 nmol N L−1 h−1; NO3− was assimilated at 0.67–18.75 nmol N L−1 h−1 and regenerated at 0.05–28.97 nmol N L−1 h−1. Observations implied that these processes were closely coupled at the regional scale and nitrogen recycling played an important role in sustaining phytoplankton growth during the summer. The [N2O], measured in water column profiles, was 10.13 ± 1.11 nmol L−1 and did not strongly diverge from atmospheric equilibrium indicating that sampled marine regions where neither a strong source nor sink of N2O 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 5 further stations, Ocean Acidification (OA) bioassay experiments were conducted to investigate the response of NH4+ oxidising and regenerating organisms to simulated OA conditions, including the implications for [N2O]. Multivariate analysis was undertaken which considered the complete bioassay dataset of measured variables describing changes in N-regeneration rate, [N2O] 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 and that a mechanistic understanding of how NH4+ oxidation, NH4+ regeneration and N2O production responded to OA could be developed. 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.

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.

Challenges in integrative approaches to modelling the marine ecosystems of the North Atlantic: Physics to fish and coasts to ocean

It has long been recognised that there are strong interactions and feedbacks between climate, upper ocean biogeochemistry and marine food webs, and also that food web structure and phytoplankton community distribution are important determinants of variability in carbon production and export from the euphotic zone. Numerical models provide a vital tool to explore these interactions, given their capability to investigate multiple connected components of the system and the sensitivity to multiple drivers, including potential future conditions. A major driver for ecosystem model development is the demand for quantitative tools to support ecosystem-based management initiatives. The purpose of this paper is to review approaches to the modelling of marine ecosystems with a focus on the North Atlantic Ocean and its adjacent shelf seas, and to highlight the challenges they face and suggest ways forward. We consider the state of the art in simulating oceans and shelf sea physics, planktonic and higher trophic level ecosystems, and look towards building an integrative approach with these existing tools. We note how the different approaches have evolved historically and that many of the previous obstacles to harmonisation may no longer be present. We illustrate this with examples from the on-going and planned modelling effort in the Integrative Modelling Work Package of the EURO-BASIN programme.

A general framework for aquatic biogeochemical models

One of the most of challenging steps in the development of coupled hydrodynamic-biogeochemical models is the combination of multiple, often incompatible computer codes that describe individual physical, chemical, biological and geological processes. This “coupling” is time-consuming, error-prone, and demanding in terms of scientific and programming expertise. The open source, Fortran-based Framework for Aquatic Biogeochemical Models addresses these problems by providing a consistent set of programming interfaces through which hydrodynamic and biogeochemical models communicate. Models are coded once to connect to FABM, after which arbitrary combinations of hydrodynamic and biogeochemical models can be made. Thus, a biogeochemical model code works unmodified within models of a chemostat, a vertically structured water column, and a three-dimensional basin. Moreover, complex biogeochemistry can be distributed over many compact, self-contained modules, coupled at run-time. By enabling distributed development and user-controlled coupling of biogeochemical models, FABM enables optimal use of the expertise of scientists, programmers and end-users.

Temporal changes in total and size-fractioned chlorophyll-a in surface waters of three provinces in the Atlantic Ocean (September to November) between 2003 and 2010

Phytoplankton total chlorophyll concentration (TCHLa) and phytoplankton size structure are two important ecological indicators in biological oceanography. Using high performance liquid chromatography (HPLC) pigment data, collected from surface waters along the Atlantic Meridional Transect (AMT), we examine temporal changes in TCHLa and phytoplankton size class (PSC: micro-, nano- and pico-phytoplankton) between 2003 and 2010 (September to November cruises only), in three ecological provinces of the Atlantic Ocean. The HPLC data indicate no significant change in TCHLa in northern and equatorial provinces, and an increase in the southern province. These changes were not significantly different to changes in TCHLa derived using satellite ocean-colour data over the same study period. Despite no change in AMT TCHLa in northern and equatorial provinces, significant differences in PSC were observed, related to changes in key diagnostic pigments (fucoxanthin, peridinin, 19′-hexanoyloxyfucoxanthin and zeaxanthin), with an increase in small cells (nano- and pico-phytoplankton) and a decrease in larger cells (micro-phytoplankton). When fitting a three-component model of phytoplankton size structure — designed to quantify the relationship between PSC and TCHLa to each AMT cruise, model parameters varied over the study period. Changes in the relationship between PSC and TCHLa have wide implications in ecology and marine biogeochemistry, and provide key information for the development and use of empirical ocean-colour algorithms. Results illustrate the importance of maintaining a time-series of in-situ observations in remote regions of the ocean, such as that acquired in the AMT programme.

RRS Discovery DY021 Cruise Report, March 2015

The Shelf Sea Biogeochemistry research programme directly relates to the delivery of the NERC Earth system science theme and aims to provide evidence that supports a number of marine policy areas and statutory requirements, such as the Marine Strategy Framework Directive and Marine and Climate Acts. The shelf seas are highly productive compared to the open ocean, a productivity that underpins more than 90 per cent of global fisheries. Their importance to society extends beyond food production to include issues of biodiversity, carbon cycling and storage, waste disposal, nutrient cycling, recreation and renewable energy resources. The shelf seas have been estimated to be the most valuable biome on Earth, but they are under considerable stress, as a result of anthropogenic nutrient loading, overfishing, habitat disturbance, climate change and other impacts. However, even within the relatively well-studied European shelf seas, fundamental biogeochemical processes are poorly understood. For example: the role of shelf seas in carbon storage; in the global cycles of key nutrients (nitrogen, phosphorus, silicon and iron); and in determining primary and secondary production, and thereby underpinning the future delivery of many other ecosystem services. Improved knowledge of such factors is not only required by marine policymakers; it also has the potential to increase the quality and cost-effectiveness of management decisions at the local, national and international levels under conditions of climate change. The Shelf Sea Biogeochemistry research programme will take a holistic approach to the cycling of nutrients and carbon and the controls on primary and secondary production in UK and European shelf seas, to increase understanding of these processes and their role in wider biogeochemical cycles. It will thereby significantly improve predictive marine biogeochemical and ecosystem models over a range of scales. The scope of the programme includes exchanges with the open ocean (transport on and off the shelf to a depth of around 500m), together with cycling, storage and release processes on the shelf slope, and air-sea exchange of greenhouse gases (carbon dioxide and nitrous oxide). The DY021 cruise is the first of the 2015 Benthic SSB cruises to investigate the 4 main ‘representative’ sites in the Celtic Sea that will represent all the various sediment types found in the whole area, these being Mud, San, Sandy-Mud and Muddy-Sand. The cruise will also carry out complimentary sampling at the Pelagic SSB programme main site called CANDYFLOSS in the central Shelf area in order to better link the Benthic and Pelagic programmes.

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.

Marine regime shifts in ocean biogeochemical models: a case study in the Gulf of Alaska

Regime shifts have been reported in many marine ecosystems, and are often expressed as an abrupt change occurring in multiple physical and biological components of the system. In the Gulf of Alaska, a regime shift in the late 1970s was observed, indicated by an abrupt increase in sea surface temperature and major shifts in the catch of many fish species. This late 1970s regime shift in the Gulf of Alaska was followed by another shift in the late 1980s, not as pervasive as the 1977 shift, but which nevertheless did not return to the prior state. A thorough understanding of the extent and mechanisms leading to such regime shifts is challenged by data paucity in time and space. We investigate the ability of a suite of ocean biogeochemistry models of varying complexity to simulate regime shifts in the Gulf of Alaska by examining the presence of abrupt changes in time series of physical variables (sea surface temperature and mixed layer depth), nutrients and biological variables (chlorophyll, primary productivity and plankton biomass) using change-point analysis. Our study demonstrates that ocean biogeochemical models are capable of simulating the late 1970s shift, indicating an abrupt increase in sea surface temperature forcing followed by an abrupt decrease in nutrients and biological productivity. This predicted shift is consistent among all the models, although some of them exhibit an abrupt transition (i.e. a significant shift from one year to the next), whereas others simulate a smoother transition. Some models further suggest that the late 1980s shift was constrained by changes in mixed layer depth. Our study demonstrates that ocean biogeochemical can successfully simulate regime shifts in the Gulf of Alaska region, thereby providing better understanding of how changes in physical conditions are propagated from lower to upper trophic levels through bottom-up controls.

Decadal reanalysis of biogeochemical indicators and fluxes in the North West European shelf-sea ecosystem

In this paper we present the first decadal reanalysis simulation of the biogeochemistry of the North West European shelf, along with a full evaluation of its skill and value. An error-characterized satellite product for chlorophyll was assimilated into a physical-biogeochemical model of the North East Atlantic, applying a localized Ensemble Kalman filter. The results showed that the reanalysis improved the model predictions of assimilated chlorophyll in 60% of the study region. Model validation metrics showed that the reanalysis had skill in matching a large dataset of in situ observations for ten ecosystem variables. Spearman rank correlations were significant and higher than 0.7 for physical-chemical variables (temperature, salinity, oxygen), ∼0.6 for chlorophyll and nutrients (phosphate, nitrate, silicate), and significant, though lower in value, for partial pressure of dissolved carbon dioxide (∼0.4). The reanalysis captured the magnitude of pH and ammonia observations, but not their variability. The value of the reanalysis for assessing environmental status and variability has been exemplified in two case studies. The first shows that between 340,000-380,000 km2 of shelf bottom waters were oxygen deficient potentially threatening bottom fishes and benthos. The second application confirmed that the shelf is a net sink of atmospheric carbon dioxide, but the total amount of uptake varies between 36-46 Tg C yr−1 at a 90% confidence level. These results indicate that the reanalysis output dataset can inform the management of the North West European shelf ecosystem, in relation to eutrophication, fishery, and variability of the carbon cycle.

Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles

Global ocean biogeochemistry models currently employed in climate change projections use highly simplified representations of pelagic food webs. These food webs do not necessarily include critical pathways by which ecosystems interact with ocean biogeochemistry and climate. Here we present a global biogeochemical model which incorporates ecosystem dynamics based on the representation of ten plankton functional types (PFTs); six types of phytoplankton, three types of zooplankton, and heterotrophic bacteria. We improved the representation of zooplankton dynamics in our model through (a) the explicit inclusion of large, slow-growing zooplankton, and (b) the introduction of trophic cascades among the three zooplankton types. We use the model to quantitatively assess the relative roles of iron vs. grazing in determining phytoplankton biomass in the Southern Ocean High Nutrient Low Chlorophyll (HNLC) region during summer. When model simulations do not represent crustacean macrozooplankton grazing, they systematically overestimate Southern Ocean chlorophyll biomass during the summer, even when there was no iron deposition from dust. When model simulations included the developments of the zooplankton component, the simulation of phytoplankton biomass improved and the high chlorophyll summer bias in the Southern Ocean HNLC region largely disappeared. Our model results suggest that the observed low phytoplankton biomass in the Southern Ocean during summer is primarily explained by the dynamics of the Southern Ocean zooplankton community rather than iron limitation. This result has implications for the representation of global biogeochemical cycles in models as zooplankton faecal pellets sink rapidly and partly control the carbon export to the intermediate and deep ocean.