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.

Phytoplankton phenology indices in coral reef ecosystems: Application to ocean-color observations in the Red Sea

Phytoplankton, at the base of the marine food web, represent a fundamental food source in coral reef ecosystems. The timing (phenology) and magnitude of the phytoplankton biomass are major determinants of trophic interactions. The Red Sea is one of the warmest and most saline basins in the world, characterized by an arid tropical climate regulated by the monsoon. These extreme conditions are particularly challenging for marine life. Phytoplankton phenological indices provide objective and quantitative metrics to characterize phytoplank- ton seasonality. The indices i.e. timings of initiation, peak, termination and duration are estimated here using 15 years (1997–2012) of remote sensing ocean-color data from the European Space Agency (ESA) Climate Change Initiative project (OC-CCI) in the entire Red Sea basin. The OC-CCI product, comprising merged and bias-corrected observations from three independent ocean-color sensors (SeaWiFS, MODIS and MERIS), and processed using the POLYMER algorithm (MERIS period), shows a significant increase in chlorophyll data cover- age, especially in the southern Red Sea during the months of summer NW monsoon. In open and reef-bound coastal waters, the performance of OC-CCI chlorophyll data is shown to be comparable with the performance of other standard chlorophyll products for the global oceans. These features have permitted us to investigate phytoplankton phenology in the entire Red Sea basin, and during both winter SE monsoon and summer NW monsoon periods. The phenological indices are estimated in the four open water provinces of the basin, and further examined at six coral reef complexes of particular socio-economic importance in the Red Sea, including Siyal Islands, Sharm El Sheikh, Al Wajh bank, Thuwal reefs, Al Lith reefs and Farasan Islands. Most of the open and deeper waters of the basin show an apparent higher chlorophyll concentration and longer duration of phyto- plankton growth during the winter period (relative to the summer phytoplankton growth period). In contrast, most of the reef-bound coastal waters display equal or higher peak chlorophyll concentrations and equal or lon- ger duration of phytoplankton growth during the summer period (relative to the winter phytoplankton growth period). The ecological and biological significance of the phytoplankton seasonal characteristics are discussed in context of ecosystem state assessment, and particularly to support further understanding of the structure and functioning of coral reef ecosystems in the Red Sea.

Arctic Oceanography - Oceanography: Atmosphere-Ocean Exchange, Biogeochemistry & Physics

The Arctic Ocean is, on average, the shallowest of Earth’s oceans. Its vast continental shelf areas, which account for approximately half of the Arctic Ocean’s total area, are heavily influenced by the surrounding land masses through river run-off and coastal erosion. As a main area of deep water formation, the Arctic is one of the main «engines» of global ocean circulation, due to large freshwater inputs, it is also strongly stratified. The Arctic Ocean’s complex oceanographic configuration is tightly linked to the atmosphere, the land, and the cryosphere. The physical dynamics not only drive important climate and global circulation patterns, but also control biogeochemical cycles and ecosystem dynamics. Current changes in Arctic sea-ice thickness and distribution, air and water temperatures, and water column stability are resulting in measurable shifts in the properties and functioning of the ocean and its ecosystems. The Arctic Ocean is forecast to shift to a seasonally ice-free ocean resulting in changes to physical, chemical, and biological processes. These include the exchange of gases across the atmosphere-ocean interface, the wind-driven ciruclation and mixing regimes, light and nutrient availability for primary production, food web dynamics, and export of material to the deep ocean. In anticipation of these changes, extending our knowledge of the present Arctic oceanography and these complex changes has never been more urgent.

Ecological controls on biogeochemical fluxes in the western Antarctic Peninsula studied with an inverse foodweb model

Sea ice in the western Antarctic Peninsula (WAP) region is both highly variable and rapidly changing. In the Palmer Station region, the ice season duration has decreased by 92 d since 1978. The sea-ice changes affect ocean stratification and freshwater balance and in turn impact every component of the polar marine ecosystem. Long-term observations from the WAP nearshore and offshore regions show a pattern of chlorophyll (Chl) variability with three to five years of negative Chl anomalies interrupted by one or two years of positive anomalies (high and low Chl regimes). Both field observations and results from an inverse food-web model show that these high and low Chl regimes differed significantly from each other, with high primary productivity and net community production (NCP) and other rates associated with the high Chl years and low rates with low Chl years. Gross primary production rates (GPP) averaged 30 mmolC.m-2.d-1 in the low Chl years and 100 mmolC.m-2.d-1 in the high Chl years. Both large and small phytoplankton were more abundant and more productive in high Chl years than in low Chl years. Similarly, krill were more important as grazers in high Chl years, but did not differ from microzooplankton in high or low Chl years. Microzooplankton did not differ between high and low Chl years. Net community production differed significantly between high and low Chl years, but mobilized a similar proportion of GPP in both high and low Chl years. The composition of the NCP was uniform in high and low Chl years. These results mphasize the importance of microbial components in the WAP plankton system and suggest that food webs dominated by small phytoplankton can have pathways that funnel production into NCP, and likely, export.

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.

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.

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.

Isotopic variation complicates analysis of trophic relations within the fish community of Plußsee: a small, deep, stratifying lake

Analysis of carbon and nitrogen stable isotopes has allowed freshwater ecologists to examine lake food webs in increasing detail. Many such studies have highlighted the existence of separate within-lake pelagic and benthic-littoral food webs but are typically conducted on large (> 10 km2) lakes, whereas the majority of lakes are actually relatively small. We used stable isotope analysis (δ13C & δ15N) to examine trophic interactions between fish and their prey in Plu�see, as an example of a small, stratifying lake, and to determine whether separate pelagic/benthic-littoral food webs could be distinguished in such systems. Our results indicate that the Plu�see food web was complicated, and due to extensive intra-annual isotopic variation in zooplankton (e.g. cladoceran δ13C annual range = 25.6�), it may be impossible to definitively assign consumers from small, eutrophic stratified lakes to pelagic or benthic-littoral food webs. We present evidence that some components of the Plu�see food web (large bream) may be subsidised by carbon of methanogenic origin.

Lipid extraction has little effect on the delta N-15 of aquatic consumers

Proper application of stable isotopes (e. g., delta N-15 and delta C-13) to food web analysis requires an understanding of all nondietary factors that contribute to isotopic variability. Lipid extraction is often used during stable isotope analysis (SIA), because synthesized lipids have a low delta C-13 and can mask the delta C-13 of a consumer's diet. Recent studies indicate that lipid extraction intended to adjust delta C-13 may also cause shifts in delta N-15, but the magnitude of and reasons for the shift are highly uncertain. We examined a large data set (n = 854) for effects of lipid extraction (using Bligh and dyer's [ 1959] chloroform-methanol solvent mixtures) on the delta N-15 of aquatic consumers. We found no effect of chemically extracting lipids on the delta N-15 of whole zooplankton, unionid mussels, and fish liver samples, and found a small increase in fish muscle delta N-15 of similar to 0.4%. We also detected a negative relationship between the shift in delta N-15 following extraction and the C:N ratio in muscle tissue, suggesting that effects of extraction were greater for tissue with lower lipid content. As long as appropriate techniques such as those from Bligh and dyer (1959) are used, effects of lipid extraction on delta N-15 of aquatic consumers need not be a major consideration in the SIA of food webs.

Trophic interactions and consequent impacts of the invasive fish Pseudorasbora parva in a native aquatic foodweb: a field investigation in the UK

Introduction of the invasive Asian cyprinid fish Pseudorasbora parva into a 0.3 ha pond in England with a fish assemblage that included Cyprinus carpio, Rutilus rutilus and Scardinius erythrophthalmus resulted in their establishment of a numerically dominant population in only 2 years; density estimates exceeded 60 ind. m(-2) and they comprised > 99% of fish present. Stable isotope analysis (SIA) revealed significant trophic overlap between P. parva, R. rutilus and C. carpio, a shift associated with significantly depressed somatic growth in R. rutilus. Despite these changes, fish community composition remained similar between the ponds. Comparison with SIA values collected from an adjacent pond free of P. parva revealed a simplified food web in P. parva presence, but with an apparent trophic position shift for several fishes, including S. erythrophthalmus which appeared to assimilate energy at a higher trophic level, probably through P. parva consumption. The marked isotopic shifts shown in all taxa in the P. parva invaded pond (C-13-enriched, N-15 depleted) were indicative of a shift to a cyanobacteria-dominated phytoplankton community. These findings provide an increased understanding of the ecological consequences of the ongoing P. parva invasion of European freshwater ecosystems.

Implications of climate change for the fishes of the British Isles

Recent climatic change has been recorded across the globe. Although environmental change is a characteristic feature of life on Earth and has played a major role in the evolution and global distribution of biodiversity, predicted future rates of climatic change, especially in temperature, are such that they will exceed any that has occurred over recent geological time. Climate change is considered as a key threat to biodiversity and to the structure and function of ecosystems that may already be subject to significant anthropogenic stress. The current understanding of climate change and its likely consequences for the fishes of Britain and Ireland and the surrounding seas are reviewed through a series of case studies detailing the likely response of several marine, diadromous and freshwater fishes to climate change. Changes in climate, and in particular, temperature have and will continue to affect fish at all levels of biological organization: cellular, individual, population, species, community and ecosystem, influencing physiological and ecological processes in a number of direct, indirect and complex ways. The response of fishes and of other aquatic taxa will vary according to their tolerances and life stage and are complex and difficult to predict. Fishes may respond directly to climate-change-related shifts in environmental processes or indirectly to other influences, such as community-level interactions with other taxa. However, the ability to adapt to the predicted changes in climate will vary between species and between habitats and there will be winners and losers. In marine habitats, recent changes in fish community structure will continue as fishes shift their distributions relative to their temperature preferences. This may lead to the loss of some economically important cold-adapted species such as Gadus morhua and Clupea harengus from some areas around Britain and Ireland, and the establishment of some new, warm-adapted species. Increased temperatures are likely to favour cool-adapted (e.g. Perca fluviatilis) and warm-adapted freshwater fishes (e.g. roach Rutilus rutilus and other cyprinids) whose distribution and reproductive success may currently be constrained by temperature rather than by cold-adapted species (e.g. salmonids). Species that occur in Britain and Ireland that are at the edge of their distribution will be most affected, both negatively and positively. Populations of conservation importance (e.g. Salvelinus alpinus and Coregonus spp.) may decline irreversibly. However, changes in food-web dynamics and physiological adaptation, for example because of climate change, may obscure or alter predicted responses. The residual inertia in climate systems is such that even a complete cessation in emissions would still leave fishes exposed to continued climate change for at least half a century. Hence, regardless of the success or failure of programmes aimed at curbing climate change, major changes in fish communities can be expected over the next 50 years with a concomitant need to adapt management strategies accordingly.

The top-down mechanism for body-mass-abundance scaling

Scaling relationships between mean body masses and abundances of species in multitrophic communities continue to be a subject of intense research and debate. The top-down mechanism explored in this paper explains the frequently observed inverse linear relationship between body mass and abundance (i.e., constant biomass) in terms of a balancing of resource biomasses by behaviorally and evolutionarily adapting foragers, and the evolutionary response of resources to this foraging pressure. The mechanism is tested using an allometric, multitrophic community model with a complex food web structure. It is a statistical model describing the evolutionary and population dynamics of tens to hundreds of species in a uniform way. Particularities of the model are the detailed representation of the evolution and interaction of trophic traits to reproduce topological food web patterns, prey switching behavior modeled after experimental observations, and the evolutionary adaptation of attack rates. Model structure and design are discussed. For model states comparable to natural communities, we find that (1) the body-mass-abundance scaling does not depend on the allometric scaling exponent of physiological rates in the form expected from the energetic equivalence rule or other bottom-up theories; (2) the scaling exponent of abundance as a function of body mass is approximately -1, independent of the allometric exponent for physiological rates assumed; (3) removal of top-down control destroys this pattern, and energetic equivalence is recovered. We conclude that the top-down mechanism is active in the model, and that it is a viable alternative to bottom-up mechanisms for controlling body-mass-abundance relations in natural communities.

Estimating trophic link density from quantitative but incomplete diet data

The trophic link density and the stability of food webs are thought to be related, but the nature of this relation is controversial. This article introduces a method for estimating the link density from diet tables which do not cover the complete food web and do not resolve all diet items to species level. A simple formula for the error of this estimate is derived. Link density is determined as a function of a threshold diet fraction below which diet items are ignored (

Universal power-law diet partitioning by marine fish and squid with surprising stability–diversity implications

A central question in community ecology is how the number of trophic links relates to community species richness. For simple dynamical food-web models, link density (the ratio of links to species) is bounded from above as the number of species increases; but empirical data suggest that it increases without bounds. We found a new empirical upper bound on link density in large marine communities with emphasis on fish and squid, using novel methods that avoid known sources of bias in traditional approaches. Bounds are expressed in terms of the diet-partitioning function (DPF): the average number of resources contributing more than a fraction f to a consumer's diet, as a function of f. All observed DPF follow a functional form closely related to a power law, with power-law exponents indepen- dent of species richness at the measurement accuracy. Results imply universal upper bounds on link density across the oceans. However, the inherently scale-free nature of power-law diet partitioning suggests that the DPF itself is a better defined characterization of network structure than link density.

Loss of functionally unique species may gradually undermine ecosystems

Functionally unique species contribute to the functional diversity of natural systems, often enhancing ecosystem functioning. An abundance of weakly interacting species increases stability in natural systems, suggesting that loss of weakly linked species may reduce stability. Any link between the functional uniqueness of a species and the strength of its interactions in a food web could therefore have simultaneous effects on ecosystem functioning and stability. Here, we analyse patterns in 213 real food webs and show that highly unique species consistently tend to have the weakest mean interaction strength per unit biomass in the system. This relationship is not a simple consequence of the interdependence of both measures on body size and appears to be driven by the empirical pattern of size structuring in aquatic systems and the trophic position of each species in the web. Food web resolution also has an important effect, with aggregation of species into higher taxonomic groups producing a much weaker relationship. Food webs with fewer unique and less weakly interacting species also show significantly greater variability in their levels of primary production. Thus, the loss of highly unique, weakly interacting species may eventually lead to dramatic state changes and unpredictable levels of ecosystem functioning.