ERSEM 15.06: a generic model for marine biogeochemistry and the ecosystem dynamics of the lower trophic levels

The ERSEM model is one of the most established ecosystem models for the lower trophic levels of the marine food-web in the scientific literature. Since its original development in the early nineties it has evolved significantly from a coastal ecosystem model for the North-Sea to a generic tool for ecosystem simulations from shelf seas to the global ocean. The current model release contains all essential elements for the pelagic and benthic part of the marine ecosystem, including the microbial food-web, the carbonate system and calcification. Its distribution is accompanied by a testing framework enabling the analysis of individual parts of the model. Here we provide a detailed mathematical description of all ERSEM components along with case-studies of mesocosm type simulations, water column implementations and a brief example of a full-scale application for the North-West European shelf. Validation against in situ data demonstrates the capability of the model to represent the marine ecosystem in contrasting environments.

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.

ERSEM 15.06: a generic model for marine biogeochemistry and the ecosystem dynamics of the lower trophic levels

The European Regional Seas Ecosystem Model (ERSEM) is one of the most established ecosystem models for the lower trophic levels of the marine food web in the scientific literature. Since its original development in the early nineties it has evolved significantly from a coastal ecosystem model for the North Sea to a generic tool for ecosystem simulations from shelf seas to the global ocean. The current model release contains all essential elements for the pelagic and benthic parts of the marine ecosystem, including the microbial food web, the carbonate system, and calcification. Its distribution is accompanied by a testing framework enabling the analysis of individual parts of the model. Here we provide a detailed mathematical description of all ERSEM components along with case studies of mesocosm-type simulations, water column implementations, and a brief example of a full-scale application for the north-western European shelf. Validation against in situ data demonstrates the capability of the model to represent the marine ecosystem in contrasting environments.

A genetic marker to separate Emiliania huxleyi (Prymnesiophyceae) morphotypes

Emiliania huxleyi (Lohm.) Hay and Mohler is a ubiquitous unicellular marine alga surrounded by an elaborate covering of calcite platelets called coccoliths. It is an important primary producer involved in oceanic biogeochemistry and climate regulation. Currently, E. huxleyi is separated into five morphotypes based on morphometric, physiological, biochemical, and immunological differences. However, a genetic marker has yet to be found to characterize these morphotypes. With the use of sequence analysis and denaturing gradient gel electrophoresis, we discovered a genetic marker that correlates significantly with the separation of the most widely recognized A and B morphotypes. Furthermore, we reveal that the A morphotype is composed of a number of distinct genotypes. This marker lies within the 3' untranslated region of a coccolith associated protein mRNA, which is implicated in regulating coccolith calcification. Consequently, we tentatively termed this marker the coccolith morphology motif.

Phaeocystis colony distribution in the North Atlantic Ocean since 1948, and interpretation of long-term changes in the Phaeocystis hotspot in the North Sea

Monitoring of Phaeocystis since 1948 during the Continuous Plankton Recorder survey indicates that over the last 5.5 decades the distribution of its colonies in the North Atlantic Ocean was not restricted to neritic waters: occurrence was also recorded in the open Atlantic regions sampled, most frequently in the spring. Apparently, environmental conditions in open ocean waters, also those far oVshore, are suitable for complete lifecycle development of colonies (the only stage recorded in the survey). In the North Sea the frequency of occurrence was also highest in spring. Its southeastern part was the Phaeocystis abundance hotspot of the whole area covered by the survey. Frequency was especially high before the 1960s and after the 1980s, i.e., in the periods when anthropogenic nutrient enrichment was relatively low. Changes in eutrophication have obviously not been a major cause of long-term Phaeocystis variation in the southeastern North Sea, where total phytoplankton biomass was related signiWcantly to river discharge. Evidence is presented for the suggestion that Phaeocystis abundance in the southern North Sea is to a large extent determined by the amount of Atlantic Ocean water Xushed in through the Dover Strait. Since Phaeocystis plays a key role in element Xuxes relevant to climate the results presented here have implications for biogeochemical models of cycling of carbon and sulphur. Sea-to-air exchange of CO2 and dimethyl sulphide (DMS) has been calculated on the basis of measurements during single-year cruises. The considerable annual variation in phytoplankton and in its Phaeocystis component reported here does not warrant extrapolation of such figures.

Nitrogen cycle feedbacks as a control on euxinia in the mid-Proterozoic ocean

Geochemical evidence invokes anoxic deep oceans until the terminal Neoproterozoic similar to 0.55 Ma, despite oxygenation of Earth's atmosphere nearly 2 Gyr earlier. Marine sediments from the intervening period suggest predominantly ferruginous (anoxic Fe(II)-rich) waters, interspersed with euxinia (anoxic H2S-rich conditions) along productive continental margins. Today, sustained biotic H2S production requires NO3- depletion because denitrifiers outcompete sulphate reducers. Thus, euxinia is rare, only occurring concurrently with (steady state) organic carbon availability when N-2-fixers dominate the production in the photic zone. Here we use a simple box model of a generic Proterozoic coastal upwelling zone to show how these feedbacks caused the mid-Proterozoic ocean to exhibit a spatial/temporal separation between two states: photic zone NO3- with denitrification in lower anoxic waters, and N-2-fixation- driven production overlying euxinia. Interchange between these states likely explains the varying H2S concentration implied by existing data, which persisted until the Neoproterozoic oxygenation event gave rise to modern marine biogeochemistry.

A mechanistic explanation of the Sargasso Sea DMS summer paradox

In the Sargasso Sea, maximum dimethylsulfide (DMS) accumulation occurs in summer, concomitant with the minimum of chlorophyll and 2 months later than its precursor, dimethylsulfoniopropionate (DMSP). This phenomenon is often referred to as the DMS "summer paradox". It has been previously suggested that the main agent triggering this pattern is increasing irradiance leading to light stress-induced DMS release from phytoplankton cells. We have developed a new model describing DMS(P) dynamics in the water column and used it to investigate how and to what extent processes other than light induced DMS exudation from phytoplankton, may contribute to the DMS summer paradox. To do this, we have conceptually divided the DMS "summer paradox" into two components: (1) the temporal decoupling between chlorophyll and DMSP and (2) the temporal decoupling between DMSP and DMS. Our results suggest that it is possible to explain the above cited patterns by means of two different dynamics, respectively: (1) a succession of phytoplankton types in the surface water and (2) the bacterially mediated DMSP(d) to DMS conversion, seasonally varying as a function of nutrient limitation. This work differs from previous modelling studies in that the presented model suggests that phytoplankton light-stress induced processes may only partially explain the summer paradox, not being able to explain the decoupling between DMSP and DMS, which is possibly the more challenging aspect of this phenomenon. Our study, therefore, provides an "alternative" explanation to the summer paradox further underlining the major role that bacteria potentially play in DMS production and fate.

Minor impact of ocean acidification to the composition of the active microbial community in an Arctic sediment

Effects of ocean acidification on the composition of the active bacterial and archaeal community within Arctic surface sediment was analysed in detail using 16S rRNA 454 pyrosequencing. Intact sediment cores were collected and exposed to one of five different pCO(2) concentrations [380 (present day), 540, 750, 1120 and 3000 atm] and RNA extracted after a period of 14 days exposure. Measurements of diversity and multivariate similarity indicated very little difference between pCO(2) treatments. Only when the highest and lowest pCO(2) treatments were compared were significant differences evident, namely increases in the abundance of operational taxonomic units most closely related to the Halobacteria and differences to the presence/absence structure of the Planctomycetes. The relative abundance of members of the classes Planctomycetacia and Nitrospira increased with increasing pCO(2) concentration, indicating that these groups may be able to take advantage of changing pH or pCO(2) conditions. The modest response of the active microbial communities associated with these sediments may be due to the low and fluctuating pore-water pH already experienced by sediment microbes, a result of the pH buffering capacity of marine sediments, or due to currently unknown factors. Further research is required to fully understand the impact of elevated CO2 on sediment physicochemical parameters, biogeochemistry and microbial community dynamics.

Temporal and spatial variation in the Nazaré Canyon (Western Iberian margin): Inter-annual and canyon heterogeneity effects on meiofauna biomass and diversity

The Nazaré Canyon on the Portuguese Margin (NE Atlantic) was sampled during spring-summer for three consecutive years (2005–2007), permitting the first inter-annual study of the meiofaunal communities at the Iberian Margin at two abyssal depths (~3500 m and ~4400 m). Using new and already published data, the meiofauna standing stocks (abundance and biomass) and nematode structural and functional diversity were investigated in relation to the sediment biogeochemistry (e.g. organic carbon, nitrogen, chlorophyll a, phaeopigments) and grain size. A conspicuous increase in sand content from 2005 to 2006 and decrease of phytodetritus at both sites, suggested the occurrence of one or more physical disturbance events. Nematode standing stocks and trophic diversity decreased after these events, seemingly followed by a recovery/recolonisation period in 2007, which was strongly correlated with an increase in the quantity and bioavailability of phytodetrital organic matter supplied. Changes in meiofauna assemblages, however, also differed between stations, likely because of the contrasting hydrodynamic and food supply conditions. Higher meiofauna and nematode abundances, biomass and trophic complexity were found at the shallowest canyon station, where the quantity, quality and bioavailability of food material were higher than at the deeper site. The present results suggest that even though inter-annual variations in the sedimentary environment can regulate the meiofauna in the abyssal Nazaré Canyon, heterogeneity between sampling locations in the canyon were more pronounced.

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.

Investigating the inter-relationships between water attenuated irradiance, primary production and DMS(P)

Both solar irradiance and primary production have been proposed as independent controls on seawater dimethyl sulphide (DMS) and dimethylsulphoniopropionate (DMSP) concentrations. However, irradiance also drives photosynthesis, and thus influences a complex set of inter-related processes that modulate marine DMS. We investigate the potential inter-relationships between the rate of primary production (carbon assimilation), water-attenuated irradiance and DMS/DMSP dynamics by applying correlation analysis to a high resolution, concurrently sampled in situ data set from a range of latitudes covering multiple biogeochemical provinces from 3 of the 4 Longhurst biogeochemical domains. The combination of primary production (PP) and underwater irradiance (Iz) within a multivariate regression model is able to explain 55% of the variance in DMS concentrations from all depths within the euphotic zone and 66% of the variance in surface DMS concentrations. Contrary to some previous studies we find a variable representing biological processes is necessary to better account for the variance in DMS. We find that the inclusion of Iz accounts for variance in DMS that is independent from the variance explained by PP. This suggests an important role for solar irradiance (beyond the influence of irradiance upon primary production) in mediating the relationship between the productivity of the ecosystem, DMS/DMSP production and ambient seawater DMS concentrations.