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.

Spatio-temporal characteristics of Agulhas leakage: a model inter-comparison study

Investigating the variability of Agulhas leakage, the volume transport of water from the Indian Ocean to the South Atlantic Ocean, is highly relevant due to its potential contribution to the Atlantic Meridional Overturning Circulation as well as the global circulation of heat and salt and hence global climate. Quantifying Agulhas leakage is challenging due to the non-linear nature of this process; current observations are insufficient to estimate its variability and ocean models all have biases in this region, even at high resolution . An Eulerian threshold integration method is developed to examine the mechanisms of Agulhas leakage variability in six ocean model simulations of varying resolution. This intercomparison, based on the circulation and thermo- haline structure at the Good Hope line, a transect to the south west of the southern tip of Africa, is used to identify features that are robust regardless of the model used and takes into account the thermohaline biases of each model. When determined by a passive tracer method, 60 % of the magnitude of Agulhas leakage is captured and more than 80 % of its temporal fluctuations, suggesting that the method is appropriate for investigating the variability of Agulhas leakage. In all simulations but one, the major driver of variability is associated with mesoscale features passing through the section. High resolution (<1/10 deg.) hindcast models agree on the temporal (2–4 cycles per year) and spatial (300–500 km) scales of these features corresponding to observed Agulhas Rings. Coarser resolution models (<1/4 deg.) reproduce similar time scale of variability of Agulhas leakage in spite of their difficulties in representing the Agulhas rings properties. A coarser resolution climate model (2 deg.) does not resolve the spatio-temporal mechanism of variability of Agulhas leakage. Hence it is expected to underestimate the contribution of Agulhas Current System to climate variability.

Spatio-temporal characteristics of Agulhas leakage: a model inter-comparison study

Investigating the variability of Agulhas leakage, the volume transport of water from the Indian Ocean to the South Atlantic Ocean, is highly relevant due to its potential contribution to the Atlantic Meridional Overturning Circulation as well as the global circulation of heat and salt and hence global climate. Quantifying Agulhas leakage is challenging due to the non-linear nature of this process; current observations are insufficient to estimate its variability and ocean models all have biases in this region, even at high resolution . An Eulerian threshold integration method is developed to examine the mechanisms of Agulhas leakage variability in six ocean model simulations of varying resolution. This intercomparison, based on the circulation and thermo- haline structure at the Good Hope line, a transect to the south west of the southern tip of Africa, is used to identify features that are robust regardless of the model used and takes into account the thermohaline biases of each model. When determined by a passive tracer method, 60 % of the magnitude of Agulhas leakage is captured and more than 80 % of its temporal fluctuations, suggesting that the method is appropriate for investigating the variability of Agulhas leakage. In all simulations but one, the major driver of variability is associated with mesoscale features passing through the section. High resolution (<1/10 deg.) hindcast models agree on the temporal (2–4 cycles per year) and spatial (300–500 km) scales of these features corresponding to observed Agulhas Rings. Coarser resolution models (<1/4 deg.) reproduce similar time scale of variability of Agulhas leakage in spite of their difficulties in representing the Agulhas rings properties. A coarser resolution climate model (2 deg.) does not resolve the spatio-temporal mechanism of variability of Agulhas leakage. Hence it is expected to underestimate the contribution of Agulhas Current System to climate variability.

Climate change-related regime shifts have altered spatial synchrony of plankton dynamics in the North Sea

During the 1980s, the North Sea plankton community underwent a well-documented ecosystem regime shift, including both spatial changes (northward species range shifts) and temporal changes (increases in the total abundances of warmer water species). This regime shift has been attributed to climate change. Plankton provide a link between climate and higher trophic-level organisms, which can forage on large spatial and temporal scales. It is therefore important to understand not only whether climate change affects purely spatial or temporal aspects of plankton dynamics, but also whether it affects spatiotemporal aspects such as metapopulation synchrony. If plankton synchrony is altered, higher trophic-level feeding patterns may be modified. A second motivation for investigating changes in synchrony is that the possibility of such alterations has been examined for few organisms, in spite of the fact that synchrony is ubiquitous and of major importance in ecology. This study uses correlation coefficients and spectral analysis to investigate whether synchrony changed between the periods 1959–1980 and 1989–2010. Twenty-three plankton taxa, sea surface temperature (SST), and wind speed were examined. Results revealed that synchrony in SST and plankton was altered. Changes were idiosyncratic, and were not explained by changes in abundance. Changes in the synchrony of Calanus helgolandicus and Para-pseudocalanus spp appeared to be driven by changes in SST synchrony. This study is one of few to document alterations of synchrony and climate-change impacts on synchrony. We discuss why climate-change impacts on synchrony may well be more common and consequential than previously recognized.

Climate change-related regime shifts have altered spatial synchrony of plankton dynamics in the North Sea

During the 1980s, the North Sea plankton community underwent a well-documented ecosystem regime shift, including both spatial changes (northward species range shifts) and temporal changes (increases in the total abundances of warmer water species). This regime shift has been attributed to climate change. Plankton provide a link between climate and higher trophic-level organisms, which can forage on large spatial and temporal scales. It is therefore important to understand not only whether climate change affects purely spatial or temporal aspects of plankton dynamics, but also whether it affects spatiotemporal aspects such as metapopulation synchrony. If plankton synchrony is altered, higher trophic-level feeding patterns may be modified. A second motivation for investigating changes in synchrony is that the possibility of such alterations has been examined for few organisms, in spite of the fact that synchrony is ubiquitous and of major importance in ecology. This study uses correlation coefficients and spectral analysis to investigate whether synchrony changed between the periods 1959–1980 and 1989–2010. Twenty-three plankton taxa, sea surface temperature (SST), and wind speed were examined. Results revealed that synchrony in SST and plankton was altered. Changes were idiosyncratic, and were not explained by changes in abundance. Changes in the synchrony of Calanus helgolandicus and Para-pseudocalanus spp appeared to be driven by changes in SST synchrony. This study is one of few to document alterations of synchrony and climate-change impacts on synchrony. We discuss why climate-change impacts on synchrony may well be more common and consequential than previously recognized.

Dynamics of phytoplankton communities during late summer around the tip of the Antarctic Peninsula

The composition and distribution of phytoplankton assemblages around the tip of the Antarctic Peninsula were studied during two summer cruises (February/March 2008 and 2009). Water samples were collected for HPLC/CHEMTAX pigment and microscopic analysis. A great spatial variability in chlorophyll a (Chl a) was observed in the study area: highest levels in the vicinity of the James Ross Island (exceeding 7 mg m−3 in 2009), intermediate values (0.5 to 2 mg m−3) in the Bransfield Strait, and low concentrations in the Weddell Sea and Drake Passage (below 0.5 mg m−3). Phytoplankton assemblages were generally dominated by diatoms, especially at coastal stations with high Chl a concentration, where diatom contribution was above 90% of total Chl a. Nanoflagellates, such as cryptophytes and/or Phaeocystis antarctica, replaced diatoms in open-ocean areas (e.g., Weddell Sea). Many species of peridinin-lacking autotrophic dinoflagellates (e.g., Gymnodinium spp.) were also important to total Chl a biomass at well-stratified stations of Bransfield Strait. Generally, water column structure was the most important environmental factor determining phytoplankton communities’ biomass and distribution. The HPLC pigment data also allowed the assessment of different physiological responses of phytoplankton to ambient light variation. The present study provides new insights about the dynamics of phytoplankton in an undersampled region of the Southern Ocean highly susceptible to global climate change.

Potential consequences of climate change for primary production and fish production in large marine ecosystems

Existing methods to predict the effects of climate change on the biomass and production of marine communities are predicated on modelling the interactions and dynamics of individual species, a very challenging approach when interactions and distributions are changing and little is known about the ecological mechanisms driving the responses of many species. An informative parallel approach is to develop size-based methods. These capture the properties of food webs that describe energy flux and production at a particular size, independent of species' ecology. We couple a physical-biogeochemical model with a dynamic, size-based food web model to predict the future effects of climate change on fish biomass and production in 11 large regional shelf seas, with and without fishing effects. Changes in potential fish production are shown to most strongly mirror changes in phytoplankton production. We project declines of 30-60% in potential fish production across some important areas of tropical shelf and upwelling seas, most notably in the eastern Indo-Pacific, the northern Humboldt and the North Canary Current. Conversely, in some areas of the high latitude shelf seas, the production of pelagic predators was projected to increase by 28-89%.

Recent Arctic Climate Change and Its Remote Forcing of Northwest Atlantic Shelf Ecosystems

During recent decades, historically unprecedented changes have been observed in the Arctic as climate warming has increased precipitation, river discharge, and glacial as well as sea-ice melting. Additionally, shifts in the Arctic's atmospheric pressure field have altered surface winds, ocean circulation, and freshwater storage in the Beaufort Gyre. These processes have resulted in variable patterns of freshwater export from the Arctic Ocean, including the emergence of great salinity anomalies propagating throughout the North Atlantic. Here, we link these variable patterns of freshwater export from the Arctic Ocean to the regime shifts observed in Northwest Atlantic shelf ecosystems. Specifically, we hypothesize that the corresponding salinity anomalies, both negative and positive, alter the timing and extent of water-column stratification, thereby impacting the production and seasonal cycles of phytoplankton, zooplankton, and higher-trophic-level consumers. Should this hypothesis hold up to critical evaluation, it has the potential to fundamentally alter our current understanding of the processes forcing the dynamics of Northwest Atlantic shelf ecosystems.

Variation among northeast Atlantic regions in the responses of zooplankton to climate change: Not all areas follow the same path

Broad scale climate forcing can interact with local environmental processes to affect the observed ecological phenomena. This causes potential problems of over-extrapolation for results from a limited number of sites or the averaging out of region-specific responses if data from too wide an area are combined. In this study, an area similar in extent to the Celtic Biscay Large Marine Ecosystem, but including off-shelf areas, was partitioned using clustering of satellite chlorophyll (chl-a) measurements. The resulting clusters were used to define areas over which to combine copepod data from the Continuous Plankton Recorder. Following filtering due to data limitations, nine regions were defined with sufficient records for analysis. These regions were consistent with known oceanographic structure in the study area. Off-shelf regions showed a progressively later timing in the seasonal peak of chl-a measurements moving northwards. Generalised additive models were used to estimate seasonal and multiannual signals in the adult and juvenile stages of Calanus finmarchicus, C. helgolandicus and the Paracalanus–Pseudocalanus group. Associations between variables (sea surface temperature (SST), phenology and annual abundance) differed among taxonomic groups, but even within taxonomic groups, relationships were not consistent across regions. For example, in the deep waters off Spain and Portugal the annual abundance of Calanus finmarchicus has a weak positive association with SST, in contrast to the pattern in most other regions. The regions defined in this study provide an objective basis for investigations into the long term dynamics of plankton populations and suggest suitable sub regions for deriving pelagic system indicators.

Challenges for implementing the Marine Strategy Framework Directive in a climate of macroecological change

Unprecedented basin-scale ecological changes are occurring in our seas. As temperature and carbon dioxide concentrations increase, the extent of sea ice is decreasing, stratification and nutrient regimes are changing and pH is decreasing. These unparalleled changes present new challenges for managing our seas, as we are only just beginning to understand the ecological manifestations of these climate alterations. The Marine Strategy Framework Directive requires all European Member States to achieve good environmental status (GES) in their seas by 2020; this means management towards GES will take place against a background of climate-driven macroecological change. Each Member State must set environmental targets to achieve GES; however, in order to do so, an understanding of large-scale ecological change in the marine ecosystem is necessary. Much of our knowledge of macroecological change in the North Atlantic is a result of research using data gathered by the Continuous Plankton Recorder (CPR) survey, a near-surface plankton monitoring programme that has been sampling in the North Atlantic since 1931. CPR data indicate that North Atlantic and North Sea plankton dynamics are responding to both climate and human-induced changes, presenting challenges to the development of pelagic targets for achievement of GES in European Seas. Thus, the continuation of long-term ecological time series such as the CPR survey is crucial for informing and supporting the sustainable management of European seas through policy mechanisms.

Non-linearities, regime shifts and recovery: The recent influence of climate on Black Sea chlorophyll

The Black Sea ecosystem experienced severe eutrophication-related degradation during the 1970s and 1980s. However, in recent years the Black Sea has shown some signs of recovery which are often attributed to a reduction in nutrient loading. Here, SeaWiFS chlorophyll a (chl a), a proxy for phytoplankton biomass, is used to investigate spatio-temporal patterns in Black Sea phytoplankton dynamics and to explore the potential role of climate in the Black Sea's recovery. Maps of chl a anomalies, calculated relative to the 8 year mean, emphasize spatial and temporal variability of phytoplankton biomass in the Black Sea, particularly between the riverine-influenced Northwest Shelf and the open Black Sea. Evolution of phytoplankton biomass has shown significant spatial variability of persistence of optimal bloom conditions between three major regions of the Black Sea. With the exception of 2001, chl a has generally decreased during our 8 year time-series. However, the winter of 2000–2001 was anomalously warm with low wind stress, resulting in reduced vertical mixing of the water column and retention of nutrients in the photic zone. These conditions were associated with anomalously high levels of chl a throughout much of the open Black Sea during the following spring and summer. The unusual climatic conditions occurring in 2001 may have triggered a shift in the Black Sea's chl a regime. The long-term significance of this recent shift is still uncertain but illustrates a non-linear response to climate forcing that makes future ecosystem changes in the pelagic Black Sea ecosystem difficult to predict.

Spatio-temporal Bayesian network models with latent variables for revealing trophic dynamics and functional networks in fisheries ecology

Ecosystems consist of complex dynamic interactions among species and the environment, the understanding of which has implications for predicting the environmental response to changes in climate and biodiversity. However, with the recent adoption of more explorative tools, like Bayesian networks, in predictive ecology, few assumptions can be made about the data and complex, spatially varying interactions can be recovered from collected field data. In this study, we compare Bayesian network modelling approaches accounting for latent effects to reveal species dynamics for 7 geographically and temporally varied areas within the North Sea. We also apply structure learning techniques to identify functional relationships such as prey–predator between trophic groups of species that vary across space and time. We examine if the use of a general hidden variable can reflect overall changes in the trophic dynamics of each spatial system and whether the inclusion of a specific hidden variable can model unmeasured group of species. The general hidden variable appears to capture changes in the variance of different groups of species biomass. Models that include both general and specific hidden variables resulted in identifying similarity with the underlying food web dynamics and modelling spatial unmeasured effect. We predict the biomass of the trophic groups and find that predictive accuracy varies with the models' features and across the different spatial areas thus proposing a model that allows for spatial autocorrelation and two hidden variables. Our proposed model was able to produce novel insights on this ecosystem's dynamics and ecological interactions mainly because we account for the heterogeneous nature of the driving factors within each area and their changes over time. Our findings demonstrate that accounting for additional sources of variation, by combining structure learning from data and experts' knowledge in the model architecture, has the potential for gaining deeper insights into the structure and stability of ecosystems. Finally, we were able to discover meaningful functional networks that were spatially and temporally differentiated with the particular mechanisms varying from trophic associations through interactions with climate and commercial fisheries.

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.

Is the Russell Cycle a true cycle? Multidecadal zooplankton and climate trends in the western English Channel

The Russell Cycle is one of the classical examples of climate influence on biological oceanography, represented as shifts in the marine plankton over several decades with warm and cool conditions. While the time-series data associated with the phenomenon indicate cyclical patterns, the question remains whether or not the Russell Cycle should be considered a “true cycle”. Zooplankton time-series data from 1924 to 2011 from the western English Channel were analysed with principal component (PC), correlation and spectral analyses to determine the dominant trends, and cyclic frequencies of the Russell Cycle indicators in relation to long-term hydroclimatic indices. PC1 accounted for 37.4% of the variability in the zooplankton data with the main contributions from non-clupeid fish larvae, southwestern zooplankton, and overall zooplankton biovolume. For PC2 (14.6% of data variance), the dominant groups were northern fish larvae, non-sardine eggs, and southern fish larvae. Sardine eggs were the major contributors to PC3 (representing 12.1% of data variance). No significant correlations were observed between the above three components and climate indices: Atlantic Multidecadal Oscillation, North Atlantic Oscillation, and local seawater temperature. Significant 44- and 29-year frequencies were observed for PC3, but the physical mechanisms driving the cycles are unclear. Harmonic analysis did not reveal any significant frequencies in the physical variables or in PCs 1 and 2. To a large extent, this is due to the dominant cycles in all datasets generally being long term (>50 years or so) and not readily resolved in the examined time frame of 88 years, hence restricting the ability to draw firm conclusions on the multidecadal relationship between zooplankton community dynamics in the western English Channel and environmental indices. Thus, the zooplankton time-series often associated and represented as the Russell Cycle cannot be concluded as being truly cyclical.

Is the Russell Cycle a true cycle? Multidecadal zooplankton and climate trends in the western English Channel

The Russell Cycle is one of the classical examples of climate influence on biological oceanography, represented as shifts in the marine plankton over several decades with warm and cool conditions. While the time-series data associated with the phenomenon indicate cyclical patterns, the question remains whether or not the Russell Cycle should be considered a “true cycle”. Zooplankton time-series data from 1924 to 2011 from the western English Channel were analysed with principal component (PC), correlation and spectral analyses to determine the dominant trends, and cyclic frequencies of the Russell Cycle indicators in relation to long-term hydroclimatic indices. PC1 accounted for 37.4% of the variability in the zooplankton data with the main contributions from non-clupeid fish larvae, southwestern zooplankton, and overall zooplankton biovolume. For PC2 (14.6% of data variance), the dominant groups were northern fish larvae, non-sardine eggs, and southern fish larvae. Sardine eggs were the major contributors to PC3 (representing 12.1% of data variance). No significant correlations were observed between the above three components and climate indices: Atlantic Multidecadal Oscillation, North Atlantic Oscillation, and local seawater temperature. Significant 44- and 29-year frequencies were observed for PC3, but the physical mechanisms driving the cycles are unclear. Harmonic analysis did not reveal any significant frequencies in the physical variables or in PCs 1 and 2. To a large extent, this is due to the dominant cycles in all datasets generally being long term (>50 years or so) and not readily resolved in the examined time frame of 88 years, hence restricting the ability to draw firm conclusions on the multidecadal relationship between zooplankton community dynamics in the western English Channel and environmental indices. Thus, the zooplankton time-series often associated and represented as the Russell Cycle cannot be concluded as being truly cyclical.