Climate-driven ecosystem changes are increasingly affecting the world’s ocean ecosystems, necessitating urgent guidance on adaptation strategies to limit or prevent catastrophic impacts. The Fisheries and Marine Ecosystem Model Intercomparison Project (FishMIP) is a network and framework that provides standardised ensemble projections of the impacts of climate change and fisheries on ocean life and the benefits that it provides to people through fisheries. Since its official launch in 2013 as a small, self-organised project within the larger Inter-Sectoral Impact Model Intercomparison Project, the FishMIP community has grown substantially and contributed to key international policy processes, such as the IPCC AR5 and AR6, and the IPBES Global Biodiversity Assessment. While not without challenges, particularly around comparing heterogeneous ecosystem models, integrating fisheries scenarios, and standardising regional-scale ecosystem models, FishMIP outputs are now being used across a variety of applications (e.g., climate change targets, fisheries management, marine conservation, Sustainable Development Goals). Over the next decade, FishMIP will focus on improving ecosystem model ensembles to provide more robust and policy-relevant projections for different regions of the world under multiple climate and societal change scenarios, and continue to be open to a broad spectrum of marine ecosystem models and modellers. FishMIP also intends to enhance leadership diversity and capacity-building to improve representation of early- and mid-career researchers from under-represented countries and ocean regions. As we look ahead, FishMIP aims to continue enhancing our understanding of how marine life and its contributions to people may change over the coming century at both global and regional scales.

As the urgency to evaluate the impacts of climate change on marine ecosystems increases, there is a need to develop robust projections and improve the uptake of ecosystem model outputs in policy and planning. Standardising input and output data is a crucial step in evaluating and communicating results, but can be challenging when using models with diverse structures, assumptions, and outputs that address region-specific issues. We developed an implementation framework and workflow to standardise the climate and fishing forcings used by regional models contributing to the Fisheries and Marine Ecosystem Model Intercomparison Project (FishMIP) and to facilitate comparative analyses across models and a wide range of regions, in line with the FishMIP 3a protocol. We applied our workflow to three case study areas-models: the Baltic Sea Mizer, Hawai’i-based Longline fisheries therMizer, and the southern Benguela ecosystem Atlantis marine ecosystem models. We then selected the most challenging steps of the workflow and illustrated their implementation in different model types and regions. Our workflow is adaptable across a wide range of regional models, from non-spatially explicit to spatially explicit and fully-depth resolved models and models that include one or several fishing fleets. This workflow will facilitate the development of regional marine ecosystem model ensembles and enhance future research on marine ecosystem model development and applications, model evaluation and benchmarking, and global-to-regional model comparisons.

The Fisheries and Marine Ecosystems Model Intercomparison Project (FishMIP) has dedicated a decade to unravelling the future impacts of climate change on marine animal biomass. FishMIP is now preparing a new simulation protocol to assess the combined effects of both climate and socio-economic changes on marine fisheries and ecosystems. This protocol will be based on the Ocean System Pathways (OSPs), a new set of socio-economic scenarios derived from the Shared Socioeconomic Pathways (SSPs) widely used by the Intergovernmental Panel on Climate Change (IPCC). The OSPs extend the SSPs to the economic, governance, management and socio-cultural contexts of large pelagic, small pelagic, benthic-demersal and emerging fisheries, as well as mariculture. Comprising qualitative storylines, quantitative model driver pathways and a “plug-in-model” framework, the OSPs will enable a heterogeneous suite of ecosystem models to simulate fisheries dynamics in a standardised way. This paper introduces this OSP framework and the simulation protocol that FishMIP will implement to explore future ocean social-ecological systems holistically, with a focus on critical issues such as climate justice, global food security, equitable fisheries, aquaculture development, fisheries management, and biodiversity conservation. Ultimately, the OSP framework is tailored to contribute to the synthesis work of the IPCC. It also aims to inform ongoing policy processes within the United Nations Food and Agriculture Organisation (FAO). Finally, it seeks to support the synthesis work of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES), with a particular focus on studying pathways relevant for the United Nations Convention on Biological Diversity (CBD).

Climate change could irreversibly modify Southern Ocean ecosystems. Marine ecosystem model (MEM) ensembles can assist policy making by projecting future changes and allowing the evaluation and assessment of alternative management approaches. However, projected future changes in total consumer biomass from the Fisheries and Marine Ecosystem Model Intercomparison Project (FishMIP) global MEM ensemble highlight an uncertain future for the Southern Ocean, indicating the need for a region-specific ensemble. A large source of model uncertainty originates from the Earth system models (ESMs) used to force FishMIP models, particularly future changes to lower trophic level biomass and sea ice coverage. To build confidence in regional MEMs as ecosystem-based management tools in a changing climate that can better account for uncertainty, we propose the development of a Southern Ocean Marine Ecosystem Model Ensemble (SOMEME) contributing to the FishMIP 2.0 regional model intercomparison initiative. One of the challenges hampering progress of regional MEM ensembles is achieving the balance of global standardised inputs with regional relevance. As a first step, we design a SOMEME simulation protocol, that builds on and extends the existing FishMIP framework, in stages that include: detailed skill assessment of climate forcing variables for Southern Ocean regions, extension of fishing forcing data to include whaling, and new simulations that assess ecological links to sea-ice processes in an ensemble of candidate regional MEMs. These extensions will help advance assessments of urgently needed climate change impacts on Southern Ocean ecosystems.

Understanding climate change impacts on global marine ecosystems and fisheries requires complex marine ecosystem models, forced by global climate projections, that can robustly detect and project changes. The Fisheries and Marine Ecosystems Model Intercomparison Project (FishMIP) uses an ensemble modelling approach to fill this crucial gap. Yet FishMIP does not have a standardised skill assessment framework to quantify the ability of member models to reproduce past observations and to guide model improvement. In this study, we apply a comprehensive model skill assessment framework to a subset of global FishMIP models that produce historical fisheries catches. We consider a suite of metrics and assess their utility in illustrating the models’ ability to reproduce observed fisheries catches. Our findings reveal improvement in model performance at both global and regional (Large Marine Ecosystem) scales from the Coupled Model Intercomparison Project Phase 5 and 6 simulation rounds. Our analysis underscores the importance of employing easily interpretable, relative skill metrics to estimate the capability of models to capture temporal variations, alongside absolute error measures to characterise shifts in the magnitude of these variations between models and across simulation rounds. The skill assessment framework developed and tested here provides a first objective assessment and a baseline of the FishMIP ensemble’s skill in reproducing historical catch at the global and regional scale. This assessment can be further improved and systematically applied to test the reliability of FishMIP models across the whole model ensemble from future simulation rounds and include more variables like fish biomass or production.

Climate change is affecting ocean temperature, acidity, currents, and primary production, causing shifts in species distributions, marine ecosystems, and ultimately fisheries. Earth system models simulate climate change impacts on physical and biogeochemical properties of future oceans under varying emissions scenarios. Coupling these simulations with an ensemble of global marine ecosystem models indicates decreasing global fish biomass with warming. However, regional projections of these impacts remain much more uncertain. Here, we employ CMIP5 and CMIP6 climate change impact projections using two Earth system models coupled with four regional and nine global marine ecosystem models in ten ocean regions to evaluate model agreement at regional scales. We find that models developed at different scales can lead to stark differences in biomass projections. On average, global models projected greater biomass declines by the end of the 21st century than regional models. For both global and regional models, greater biomass declines were projected using CMIP6 than CMIP5 simulations. Global models projected biomass declines in 86% of CMIP5 simulations for ocean regions compared to 50% for regional models in the same ocean regions. In CMIP6 simulations, all global model simulations projected biomass declines in ocean regions by 2100, while regional models projected biomass declines in 67% of the ocean region simulations. Our analysis suggests that improved understanding of the causes of differences between global and regional marine ecosystem model climate change projections is needed, alongside observational evaluation of modelled responses.

Robust projections of future trends in global fish biomass, production and catches under different fishing scenarios are needed to inform fisheries policy in a changing climate. Trust in future projections, however, relies on establishing that the models used can accurately simulate past relationships between exploitation rates, catches and ecosystem states. Here we use fisheries catch and catch-only assessment models in combination with effort data to estimate regional fishing exploitation levels (defined as the fishing mortality relative to fishing mortality at maximum sustainable yield, F/FMSY). These estimates are given for large pelagic, forage and demersal fish types across all large marine ecosystems and the high seas between 1961 and 2004; and with a ‘ramp-up’ between 1841-1960. We find that global exploitation rates for both large pelagic and demersal fish are consistently higher than for forage fish and reached their peaks in the late 1980s. We use the exploitation rates to globally simulate historical fishing patterns in a mechanistic fish community model. We find a good match between model and reconstructed fisheries catch, both for total catch as well as catch distribution by functional type. Simulations show a clear deviation from an unfished model state, with a 25% reduction in fish biomass in large pelagic and demersal fish in shelf regions in the most recent years and a 50% increase in forage fish, primarily due to the release of predation pressure. These results can set a baseline from which the effect of climate change relative to fishing could be estimated.

There is an urgent need for models that can robustly detect past and project future ecosystem changes and risks to the services that they provide to people. The Fisheries and Marine Ecosystem Model Intercomparison Project (FishMIP) was established to develop model ensembles for projecting long-term impacts of climate change on fisheries and marine ecosystems while informing policy at spatio-temporal scales relevant to the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) framework. While contributing FishMIP models have improved over time, large uncertainties in projections remain, particularly in coastal and shelf seas where most of the world’s fisheries occur. Furthermore, previous FishMIP climate impact projections have mostly ignored fishing activity due to a lack of standardized historical and scenario-based human activity forcing and uneven capabilities to dynamically model fisheries across the FishMIP community. This, in addition to underrepresentation of coastal processes, has limited the ability to evaluate the FishMIP ensemble’s ability to adequately capture past states - a crucial step for building confidence in future projections. To address these issues, we have developed two parallel simulation experiments (FishMIP 2.0) on: 1) model evaluation and detection of past changes and 2) future scenarios and projections. Key advances include historical climate forcing, that captures oceanographic features not previously resolved, and standardized fishing forcing to systematically test fishing effects across models. FishMIP 2.0 is a key step towards a detection and attribution framework for marine ecosystem change at regional and global scales, and towards enhanced policy relevance through increased confidence in future ensemble projections.

Aquatic ectotherms often attain smaller body sizes at higher temperatures. By analysing ~15,000 coastal-reef fish surveys across a 15oC spatial sea surface temperature (SST) gradient, we found that the mean length of fish length in communities decreases by ~5% for each 1oC temperature increase across space. This equated to a 50% decrease in mean length from 14 to 29oC mean annual SST. We found that trophic guild composition shifts from domination by herbivores and planktivores in the tropics, to invertivores and piscivores in cooler waters. By investigating the contribution of trophic composition to community-level mean length, we show ~25% of temperature-related changes could be attributed to trophic composition at the warmest sites, but <1% at colder temperatures. Our findings suggest that small changes in temperature will lead to large changes in fish community body sizes, driven both by community trophic composition in warm waters and mean sizes within trophic guilds.

Aquatic ectotherms often reach smaller body sizes at higher temperatures. By analysing ~15,000 coastal-reef fish surveys across a 15oC spatial sea surface temperature gradient, we show mean community-level fish length decreased by ~5% for each 1oC warming, equating to a decrease in mean community-level body-length by 50%, or mean weight by 90% from 14-29oC. We further show dominant trophic guild composition shifts from invertivores and piscivores, to herbivores and planktivores, as water temperature increases. By investigating the contribution of trophic-composition to overall community-mean-length, we show ~25% of temperature-related changes could be attributed to trophic-composition at the warmest sites, but close to zero at colder temperatures. Our findings suggest that small changes in temperature will lead to large changes in fish community body sizes, driven both by changes in community trophic composition in warm waters, and by changes to the average body sizes of fish within trophic guilds across all temperatures.

Fed aquaculture is one of the fastest growing and most valuable food production industries. The efficiency with which farmed fish convert feed into biomass influences both environmental impact and economic revenue. Salmonid species, such as king salmon (Oncorhynchus tshawytscha), exhibit high levels of plasticity in vital rates such as feed intake and growth rates. Accurate estimations of individual variability in vital rates are important for production management. The use of mean trait values to evaluate feeding and growth performance can mask individual-level differences that potentially contribute to inefficiencies. Here, we apply an integral projection model (IPM) to investigate individual variation in growth performance of 1625 individually tagged king salmon fed one of three distinct rations and tracked over 276 days. To capture the observed sigmoidal growth, we compared a non-linear mixed-effects (logistic) model to a linear regression model used within the IPM framework. Ration significantly influenced several aspects of growth. Mean final body mass and mean growth rate increased with ration, however, variance in body mass and feed intake also increased significantly over time. Trends in body mass mean and variance were captured by both logistic and linear models, suggesting the linear model to be suitable for use in the IPM. Higher rations resulted in a decreasing proportion of individuals reaching the cohort’s mean size or larger by the end of the experiment. This suggests that, in our trial, feeding to satiation did not produce the desired effects of efficient and uniform growth in juvenile king salmon. While monitoring individuals through time is challenging in commercial aquaculture settings, recent technological advances combined with an IPM approach could provide new scope for tracking growth performance in experimental and farmed populations. The IPM framework also allows the exploration of other size-dependent processes affecting vital rate functions, such as competition and mortality.

Climate change and fisheries exploitation are dramatically changing the species composition, abundances, and size spectra of fish communities. We explore whether variation in abundance-size spectra, a widely studied ecosystem feature, is influenced by a critical parameter thought to govern the shape of size-structured ecosystems—the relationship between the sizes of predators and their prey (predator-prey mass ratios, or PPMRs). PPMR estimates are lacking for vast numbers of fish species, including at the broader trophic guild scale. Using measurements of 8,128 prey items in gut contents of 97 reef fish species, we established PPMRs for four major trophic guilds (piscivores, invertivores, planktivores and herbivores) using linear mixed effects models. To assess theoretical predictions that higher mean community-level PPMR leads to shallower size spectrum slopes, we compared observations of mean community-level PPMR with size spectrum slopes for coastal reef sites distributed around Australia. PPMRs of individual fishes were remarkably high (median ~71,000), with significant variation between different trophic guilds (~890 for piscivores; ~83,000 for planktivores), and ~8,700 for whole communities. Community-level PPMRs were positively related to size spectrum slopes, broadly consistent with theory, however, this pattern was also influenced by the latitudinal temperature gradient. Tropical reefs showed a stronger relationship between community-level PPMRs and community size spectrum slopes than temperate reefs. The extent that these patterns apply outside Australia, and consequences for community structure and dynamics, are key areas for future investigation.