Basking sharks and oceanographic fronts: quantifying associations in the north-east Atlantic

Understanding the mechanisms linking oceanographic processes and marine vertebrate habitat use is critical to effective management of populations of conservation concern. The basking shark Cetorhinus maximus has been shown to associate with oceanographic fronts – physical interfaces at the transitions between water masses – to exploit foraging opportunities resulting from aggregation of zooplankton. However, the scale, significance and variability of these observed associations have not yet been established. Here, we quantify the influence of mesoscale (10s – 100s km) frontal activity on habitat use over timescales of weeks to months. We use animal-mounted archival tracking with composite front mapping via Earth Observation (EO) remote sensing to provide an oceanographic context to individual shark movements. We investigate levels of association with fronts occurring over two spatio-temporal scales, (i) broad-scale seasonally persistent frontal zones and (ii) contemporaneous mesoscale thermal and chl-a fronts. Using random walk simulations and logistic regression within an iterative generalised linear mixed modelling (GLMM) framework, we find that seasonal front frequency is a significant predictor of shark presence. Temporally-matched oceanographic metrics also indicate that sharks demonstrate a preference for productive regions, and associate with contemporaneous thermal and chl-a fronts more frequently than could be expected at random. Moreover, we highlight the importance of cross-frontal temperature change and persistence, which appear to interact to affect the degree of prey aggregation along thermal fronts. These insights have clear implications for understanding the preferred habitats of basking sharks in the context of anthropogenic threat management and marine spatial planning in the northeast Atlantic.

The assessment of a global marine ecosystem model on the basis of emergent properties and ecosystem function: a case study with ERSEM

Ecosystem models are often assessed using quantitative metrics of absolute ecosystem state, but these model-data comparisons are disproportionately vulnerable to discrepancies in the location of important circulation features. An alternative method is to demonstrate the models capacity to represent ecosystem function; the emergence of a coherent natural relationship in a simulation indicates that the model may have an appropriate representation of the ecosystem functions that lead to the emergent relationship. Furthermore, as emergent properties are large-scale properties of the system, model validation with emergent properties is possible even when there is very little or no appropriate data for the region under study, or when the hydrodynamic component of the model differs significantly from that observed in nature at the same location and time. A selection of published meta-analyses are used to establish the validity of a complex marine ecosystem model and to demonstrate the power of validation with emergent properties. These relationships include the phytoplankton community structure, the ratio of carbon to chlorophyll in phytoplankton and particulate organic matter, the ratio of particulate organic carbon to particulate organic nitrogen and the stoichiometric balance of the ecosystem. These metrics can also inform aspects of the marine ecosystem model not available from traditional quantitative and qualitative methods. For instance, these emergent properties can be used to validate the design decisions of the model, such as the range of phytoplankton functional types and their behaviour, the stoichiometric flexibility with regards to each nutrient, and the choice of fixed or variable carbon to nitrogen ratios.

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.

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.