Predicting future thermal habitat suitability of competing native and invasive fish species: from metabolic scope to oceanographic modelling

Global increase in sea temperatures has been suggested to facilitate the incoming and spread of tropical invaders. The increasing success of these species may be related to their higher physiological performance compared with indigenous ones. Here, we determined the effect of temperature on the aerobic metabolic scope (MS) of two herbivorous fish species that occupy a similar ecological niche in the Mediterranean Sea: the native salema (Sarpa salpa) and the invasive marbled spinefoot (Siganus rivulatus). Our results demonstrate a large difference in the optimal temperature for aerobic scope between the salema (21.8°C) and the marbled spinefoot (29.1°C), highlighting the importance of temperature in determining the energy availability and, potentially, the distribution patterns of the two species. A modelling approach based on a present-day projection and a future scenario for oceanographic conditions was used to make predictions about the thermal habitat suitability (THS, an index based on the relationship between MS and temperature) of the two species, both at the basin level (the whole Mediterranean Sea) and at the regional level (the Sicilian Channel, a key area for the inflow of invasive species from the Eastern to the Western Mediterranean Sea). For the present-day projection, our basin-scale model shows higher THS of the marbled spinefoot than the salema in the Eastern compared with the Western Mediterranean Sea. However, by 2050, the THS of the marbled spinefoot is predicted to increase throughout the whole Mediterranean Sea, causing its westward expansion. Nevertheless, the regional-scale model suggests that the future thermal conditions of Western Sicily will remain relatively unsuitable for the invasive species and could act as a barrier for its spread westward. We suggest that metabolic scope can be used as a tool to evaluate the potential invasiveness of alien species and the resilience to global warming of native species.

Modelling the effects of climate change on the distribution and production of marine fishes: accounting for trophic interactions in a dynamic bioclimate envelope model

Climate change has already altered the distribution of marine fishes. Future predictions of fish distributions and catches based on bioclimate envelope models are available, but to date they have not considered interspecific interactions. We address this by combining the species-based Dynamic Bioclimate Envelope Model (DBEM) with a size-based trophic model. The new approach provides spatially and temporally resolved predictions of changes in species' size, abundance and catch potential that account for the effects of ecological interactions. Predicted latitudinal shifts are, on average, reduced by 20% when species interactions are incorporated, compared to DBEM predictions, with pelagic species showing the greatest reductions. Goodness-of-fit of biomass data from fish stock assessments in the North Atlantic between 1991 and 2003 is improved slightly by including species interactions. The differences between predictions from the two models may be relatively modest because, at the North Atlantic basin scale, (i) predators and competitors may respond to climate change together; (ii) existing parameterization of the DBEM might implicitly incorporate trophic interactions; and/or (iii) trophic interactions might not be the main driver of responses to climate. Future analyses using ecologically explicit models and data will improve understanding of the effects of inter-specific interactions on responses to climate change, and better inform managers about plausible ecological and fishery consequences of a changing environment.

The meso-scale response of subarctic North Pacific seabird community structure to lower trophic level abundance and diversity

Understanding the mechanisms that structure communities and influence biodiversity are fundamental goals of ecology. To test the hypothesis that the abundance and diversity of upper-trophic level predators (seabirds) is related to the underlying abundance and diversity of their prey (zooplankton) and ecosystem-wide energy availability (primary production), we initiated a monitoring program in 2002 that jointly and repeatedly surveys seabird and zooplankton populations across a 7,500 km British Columbia-Bering Sea-Japan transect. Seabird distributions were recorded by a single observer (MH) using a strip-width technique, mesozooplankton samples were collected with a Continuous Plankton Recorder, and primary production levels were derived using the appropriate satellite parameters and the Vertically Generalized Production Model (Behrenfeld and Falkowski 1997). Each trophic level showed clear spatio-temporal patterns over the course of the study. The strongest relationship between seabird abundance and diversity and the lower trophic levels was observed in March/April ('spring') and significant relationships were also found through June/July ('summer'). No discernable relationships were observed during the September/October ('fall') months. Overall, mesozooplankton abundance and biomass explained the dominant portion of seabird abundance and diversity indices (richness, Simpson's Index, and evenness), while primary production was only related to seabird richness. These findings underscore the notion that perturbations of ocean productivity and lower trophic level ecosystem constituents influenced by climate change, such as shifts in timing (phenology) and synchronicity (match-mismatch), could impart far-reaching consequences throughout the marine food web.

Scale-dependent foraging ecology of a marine top predator modelled using passive acoustic data

1.Understanding which environmental factors drive foraging preferences is critical for the development of effective management measures, but resource use patterns may emerge from processes that occur at different spatial and temporal scales. Direct observations of foraging are also especially challenging in marine predators, but passive acoustic techniques provide opportunities to study the behaviour of echolocating species over a range of scales. 2.We used an extensive passive acoustic data set to investigate the distribution and temporal dynamics of foraging in bottlenose dolphins using the Moray Firth (Scotland, UK). Echolocation buzzes were identified with a mixture model of detected echolocation inter-click intervals and used as a proxy of foraging activity. A robust modelling approach accounting for autocorrelation in the data was then used to evaluate which environmental factors were associated with the observed dynamics at two different spatial and temporal scales. 3.At a broad scale, foraging varied seasonally and was also affected by seabed slope and shelf-sea fronts. At a finer scale, we identified variation in seasonal use and local interactions with tidal processes. Foraging was best predicted at a daily scale, accounting for site specificity in the shape of the estimated relationships. 4.This study demonstrates how passive acoustic data can be used to understand foraging ecology in echolocating species and provides a robust analytical procedure for describing spatio-temporal patterns. Associations between foraging and environmental characteristics varied according to spatial and temporal scale, highlighting the need for a multi-scale approach. Our results indicate that dolphins respond to coarser scale temporal dynamics, but have a detailed understanding of finer-scale spatial distribution of resources.

Modelling large-scale CO2 leakages in the North Sea

A three dimensional hydrodynamic model with a coupled carbonate speciation sub-model is used to simulate large additions of CO2into the North Sea, representing leakages at potential carbon sequestration sites. A range of leakage scenarios are conducted at two distinct release sites, allowing an analysis of the seasonal, inter-annual and spatial variability of impacts to the marine ecosystem. Seasonally stratified regions are shown to be more vulnerable to CO2release during the summer as the added CO2remains trapped beneath the thermocline, preventing outgasing to the atmosphere. On average, CO2 injected into the northern North Sea is shown to reside within the water column twice as long as an equivalent addition in the southern North Sea before reaching the atmosphere. Short-term leakages of 5000 tonnes CO2over a single day result in substantial acidification at the release sites (up to -1.92 pH units), with significant perturbations (greater than 0.1 pH units) generally confined to a 10 km radius. Long-term CO2leakages sustained for a year may result in extensive plumes of acidified seawater, carried by major advective pathways. Whilst such scenarios could be harmful to marine biota over confined spatial scales, continued unmitigated CO2emissions from fossil fuels are predicted to result in greater and more long-lived perturbations to the carbonate system over the next few decades.

Are Fast Responses More Random? Testing the Effect of Response Time on Scale in an Online Choice Experiment

Scepticism over stated preference surveys conducted online revolves around the concerns over “professional respondents” who might rush through the questionnaire without sufficiently considering the information provided. To gain insight on the validity of this phenomenon and test the effect of response time on choice randomness, this study makes use of a recently conducted choice experiment survey on ecological and amenity effects of an offshore windfarm in the UK. The positive relationship between self-rated and inferred attribute attendance and response time is taken as evidence for a link between response time and cognitive effort. Subsequently, the generalised multinomial logit model is employed to test the effect of response time on scale, which indicates the weight of the deterministic relative to the error component in the random utility model. Results show that longer response time increases scale, i.e. decreases choice randomness. This positive scale effect of response time is further found to be non-linear and wear off at some point beyond which extreme response time decreases scale. While response time does not systematically affect welfare estimates, higher response time increases the precision of such estimates. These effects persist when self-reported choice certainty is controlled for. Implications of the results for online stated preference surveys and further research are discussed.

Influence of light in the mixed-layer on the parameters of a three-component model of phytoplankton size class

Phytoplankton size structure is an important indicator of the state of the pelagic ecosystem. Stimulated by the paucity of in situ observations on size structure, and by the sampling advantages of autonomous remote platforms, new efforts are being made to infer the size-structure of the phytoplankton from oceanographic variables that may be measured at high temporal and spatial resolution, such as total chlorophyll concentration. Large-scale analysis of in situ data has revealed coherent relationships between size-fractionated chlorophyll and total chlorophyll that can be quantified using the three-component model of Brewin et al. (2010). However, there are variations surrounding these general relationships. In this paper, we first revise the three-component model using a global dataset of surface phytoplankton pigment measurements. Then, using estimates of the average irradiance in the mixed-layer, we investigate the influence of ambient light on the parameters of the three-component model. We observe significant relationships between model parameters and the average irradiance in the mixed-layer, consistent with ecological knowledge. These relationships are incorporated explicitly into the three-component model to illustrate variations in the relationship between size-structure and total chlorophyll, ensuing from variations in light availability. The new model may be used as a tool to investigate modifications in size-structure in the context of a changing climate.

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.

Quantifying pursuit diving seabirds' associations with fine-scale physical features in a high tidal energy environment

The rapid increase in the number of tidal stream turbine arrays will create novel and unprecedented levels of anthropogenic activity within habitats characterized by horizontal current speeds exceeding 2 ms−1. However, the potential impacts on pursuit‐diving seabirds exploiting these tidal stream environments remain largely unknown. Identifying similarities between the fine‐scale physical features (100s of metres) suitable for array installations, and those associated with foraging pursuit‐diving seabirds, could identify which species are most vulnerable to either collisions with moving components, or displacement from these installations. A combination of vessel‐based observational surveys, Finite Volume Community Ocean Model outputs and hydroacoustic seabed surveys provided concurrent measures of foraging distributions and physical characteristics at a fine temporal (15 min) and spatial (500 m) resolution across a tidal stream environment suitable for array installations, during both breeding and non‐breeding seasons. These data sets were then used to test for associations between foraging pursuit‐diving seabirds (Atlantic puffins Fratercula arctica, black guillemots Cepphus grylle, common guillemots Uria aalge, European shags Phalacrocorax aristotelis) and physical features. These species were associated with areas of fast horizontal currents, slow horizontal currents, high turbulence, downward vertical currents and also hard–rough seabeds. The identity and strength of associations differed among species, and also within species between seasons, indicative of interspecific and intraspecific variations in habitat use. However, Atlantic puffins were associated particularly strongly with areas of fast horizontal currents during breeding seasons, and European shags with areas of rough–hard seabeds and downward vertical currents during non‐breeding seasons. Synthesis and applications. Atlantic puffins’ strong association with fast horizontal current speeds indicates that they are particularly likely to interact with installations during breeding seasons. Any post‐installation monitoring and mitigation measures should therefore focus on this species and season. The multi‐species associations with high turbulence and downward vertical currents, which often coincide with fast horizontal current speeds, also highlight useful pre‐installation mitigation measures via the omission of devices from these areas, reducing the overall likelihood of interactions. Environmental impact assessments (EIA) generally involve once‐a‐month surveys across 2‐year periods. However, the approaches used in this study show that more focussed surveys can greatly benefit management strategies aiming to reduce the likelihood of negative impacts by facilitating the development of targeted mitigation measures. It is therefore recommended that these approaches contribute towards EIA within development sites.

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.