Coupling planktonic ecosystem and fisheries food web models for a pelagic ecosystem: Description and validation for the subarctic Pacific

Exploring multiple stressor effects with Ecopath, Ecosim, and Ecospace: Research designs, modeling techniques, and future directions

Understanding the cumulative effects of multiple stressors is a research priority in environmental science. Ecological models are a key component of tackling this challenge because they can simulate interactions between the components of an ecosystem. Here, we ask, how has the popular modeling platform Ecopath with Ecosim (EwE) been used to model human impacts related to climate change, land and sea use, pollution, and invasive species? We conducted a literature review encompassing 166 studies covering stressors other than fishing mostly in aquatic ecosystems. The most modeled stressors were physical climate change (60 studies), species introductions (22), habitat loss (21), and eutrophication (20), using a range of modeling techniques. Despite this comprehensive coverage, we identified four gaps that must be filled to harness the potential of EwE for studying multiple stressor effects. First, only 12% of studies investigated three or more stressors, with most studies focusing on single stressors. Furthermore, many studies modeled only one of many pathways through which each stressor is known to affect ecosystems. Second, various methods have been applied to define environmental response functions representing the effects of single stressors on species groups. These functions can have a large effect on the simulated ecological changes, but best practices for deriving them are yet to emerge. Third, human dimensions of environmental change - except for fisheries - were rarely considered. Fourth, only 3% of studies used statistical research designs that allow attribution of simulated ecosystem changes to stressors' direct effects and interactions, such as factorial (computational) experiments. None made full use of the statistical possibilities that arise when simulations can be repeated many times with controlled changes to the inputs. We argue that all four gaps are feasibly filled by integrating ecological modeling with advances in other subfields of environmental science and in computational statistics.

Do the Brazilian sardine commercial landings respond to local ocean circulation?

It has been reported that sea surface temperature (SST) anomalies, flow intensity and mesoscale ocean processes, all affect sardine production, both in eastern and western boundary current systems. Here we tested the hypothesis whether extreme high and low commercial landings of the Brazilian sardine fisheries in the South Brazil Bight (SBB) are sensitive to different oceanic conditions. An ocean model (ROMS) and an individual based model (Ichthyop) were used to assess the relationship between oceanic conditions during the spawning season and commercial landings of the Brazilian sardine one year later. Model output was compared with remote sensing and analysis data showing good consistency. Simulations indicate that mortality of eggs and larvae by low temperature prior to maximum and minimum landings are significantly higher than mortality caused by offshore advection. However, when periods of maximum and minimum sardine landings are compared with respect to these causes of mortality no significant differences were detected. Results indicate that mortality caused by prevailing oceanic conditions at early life stages alone can not be invoked to explain the observed extreme commercial landings of the Brazilian sardine. Likely influencing factors include starvation and predation interacting with the strategy of spawning “at the right place and at the right time”.

Formulation, General Features and Global Calibration of a Bioenergetically-Constrained Fishery Model

Human exploitation of marine resources is profoundly altering marine ecosystems, while climate change is expected to further impact commercially-harvested fish and other species. Although the global fishery is a highly complex system with many unpredictable aspects, the bioenergetic limits on fish production and the response of fishing effort to profit are both relatively tractable, and are sure to play important roles. Here we describe a generalized, coupled biological-economic model of the global marine fishery that represents both of these aspects in a unified framework, the BiOeconomic mArine Trophic Size-spectrum (BOATS) model. BOATS predicts fish production according to size spectra as a function of net primary production and temperature, and dynamically determines harvest spectra from the biomass density and interactive, prognostic fishing effort. Within this framework, the equilibrium fish biomass is determined by the economic forcings of catchability, ex-vessel price and cost per unit effort, while the peak harvest depends on the ecosystem parameters. Comparison of a large ensemble of idealized simulations with observational databases, focusing on historical biomass and peak harvests, allows us to narrow the range of several uncertain ecosystem parameters, rule out most parameter combinations, and select an optimal ensemble of model variants. Compared to the prior distributions, model variants with lower values of the mortality rate, trophic efficiency, and allometric constant agree better with observations. For most acceptable parameter combinations, natural mortality rates are more strongly affected by temperature than growth rates, suggesting different sensitivities of these processes to climate change. These results highlight the utility of adopting large-scale, aggregated data constraints to reduce model parameter uncertainties and to better predict the response of fisheries to human behaviour and climate change.

Southern Ocean biological iron cycling in the pre-whaling and present ecosystems

This study aimed to create the first model of biological iron (Fe) cycling in the Southern Ocean food web. Two biomass mass-balanced Ecopath models were built to represent pre- and post-whaling ecosystem states (1900 and 2008). Functional group biomasses (tonnes wet weight km^-2) were converted to biogenic Fe pools (kg Fe km^-2) using published Fe content ranges. In both models, biogenic Fe pools and consumption in the pelagic Southern Ocean were highest for plankton and small nektonic groups. The production of plankton biomass, particularly unicellular groups, accounted for the highest annual Fe demand. Microzooplankton contributed most to biological Fe recycling, followed by carnivorous zooplankton and krill. Biological Fe recycling matched previous estimates, and, under most conditions, could entirely meet the Fe demand of bacterioplankton and phytoplankton. Iron recycling by large baleen whales was reduced 10-fold by whaling between 1900 and 2008. However, even under the 1900 scenario, the contribution of whales to biological Fe recycling was negligible compared with that of planktonic consumers. These models are a first step in examining oceanic-scale biological Fe cycling, highlighting gaps in our present knowledge and key questions for future research on the role of marine food webs in the cycling of trace elements in the sea.This article is part of the themed issue 'Biological and climatic impacts of ocean trace element chemistry'.

Non-deterministic modelling of food-web dynamics

A novel approach to model food-web dynamics, based on a combination of chance (randomness) and necessity (system constraints), was presented by Mullon et al. in 2009. Based on simulations for the Benguela ecosystem, they concluded that observed patterns of ecosystem variability may simply result from basic structural constraints within which the ecosystem functions. To date, and despite the importance of these conclusions, this work has received little attention. The objective of the present paper is to replicate this original model and evaluate the conclusions that were derived from its simulations. For this purpose, we revisit the equations and input parameters that form the structure of the original model and implement a comparable simulation model. We restate the model principles and provide a detailed account of the model structure, equations, and parameters. Our model can reproduce several ecosystem dynamic patterns: pseudo-cycles, variation and volatility, diet, stock-recruitment relationships, and correlations between species biomass series. The original conclusions are supported to a large extent by the current replication of the model. Model parameterisation and computational aspects remain difficult and these need to be investigated further. Hopefully, the present contribution will make this approach available to a larger research community and will promote the use of non-deterministic-network-dynamics models as ‘null models of food-webs’ as originally advocated.