A document by Tyler Eddy. Click on the document to view its contents.

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 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.

Eddies play a crucial role in shaping ocean dynamics by affecting material transport, and generating spatio-temporal heterogeneity. However, how eddies at different scales modulate biogeochemical transformation rates remains an open question. Applying a multi-scale decomposition to a numerical simulation, we investigate the respective impact of mesoscale and submesoscale eddies on nutrient transport and biogeochemical cycling in the California Current System. First, the non-linear nature of biological nutrient uptake results in a 50% reduction in primary production in the presence of eddies. Second, eddies shape the vertical transport of nutrients with a strong compensation between mesoscale and submesoscale. Third, the eddy effect on nutrient uptake is controlled by the covariance of temperature, nutrient and phytoplankton fluctuations caused by eddies. Our results shed new light on the tight interaction between non-linear fluid dynamics and ecosystem processes in realistic eddy regimes, highlighting the importance of both mesoscale and submesoscale variability.

The High Seas, lying beyond the boundaries of nations’ Exclusive Economic Zones, cover the majority of the ocean surface and host roughly two thirds of marine primary production. Yet, only a small fraction of global wild fish catch comes from the High Seas, despite intensifying industrial fishing efforts. The surprisingly small fish catch could reflect economic features of the High Seas - such as the difficulty and cost of fishing in remote parts of the ocean surface - or ecological features resulting in a small biomass of fish relative to primary production. We use the coupled biological-economic model BOATS to estimate contributing factors, comparing observed catches with simulations where: (i) fishing cost depends on distance from shore and seafloor depth; (ii) catchability depends on seafloor depth or vertical habitat extent; (iii) regions with micronutrient limitation have reduced biomass production; (iv) the trophic transfer of energy from primary production to demersal food webs depends on depth; and (v) High Seas biomass migrates to coastal regions. Our results suggest that the most important features are ecological: demersal fish communities receive a large proportion of primary production in shallow waters, but very little in deep waters due to respiration by small organisms throughout the water column. Other factors play a secondary role, with migrations having a potentially large but uncertain role, and economic factors having the smallest effects. Our results stress the importance of properly representing the High Seas biomass in future fisheries projections, and clarify their limited role in global food provision.

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.

Release of iron (Fe) from continental shelves is a major source of this limiting nutrient for phytoplankton in the open ocean, including productive Eastern Boundary Upwelling Systems. The mechanisms governing the transport and fate of Fe along continental margins remain poorly understood, reflecting interaction of physical and biogeochemical processes that are crudely represented by global ocean biogeochemical models. Here, we use a submesoscale-permitting physical-biogeochemical model to investigate processes governing the delivery of shelf-derived Fe to the open ocean along the northern U.S. West Coast. We find that a significant fraction (∼20%) of the Fe released by sediments on the shelf is transported offshore, fertilizing the broader Northeast Pacific Ocean. This transport is governed by two main pathways that reflect interaction between the wind-driven ocean circulation and Fe release by low-oxygen sediments: the first in the surface boundary layer during upwelling events; the second in the bottom boundary layer, associated with pervasive interactions of the poleward California Undercurrent with bottom topography. In the water column interior, transient and standing eddies strengthen offshore transport, counteracting the onshore pull of the mean upwelling circulation. Several hot-spots of intense Fe delivery to the open ocean are maintained by standing meanders in the mean current and enhanced by transient eddies and seasonal oxygen depletion. Our results highlight the importance of fine-scale dynamics for the transport of Fe and shelf-derived elements from continental margins to the open ocean, and the need to improve representation of these processes in biogeochemical models used for climate studies.

The coastal ocean contributes to regulating atmospheric greenhouse gas concentrations by taking up carbon dioxide (CO2) and releasing nitrous oxide (N2O) and methane (CH4). Major advances have improved our understanding of the coastal air-sea exchanges of these three gasses since the first phase of the Regional Carbon Cycle Assessment and Processes (RECCAP in 2013), but a comprehensive view that integrates the three gasses at the global scale is still lacking. In this second phase (RECCAP2), we quantify global coastal ocean fluxes of CO2, N2O and CH4 using an ensemble of global gap-filled observation-based products and ocean biogeochemical models. The global coastal ocean is a net sink of CO2 in both observational products and models, but the magnitude of the median net global coastal uptake is ~60% larger in models (-0.72 vs. -0.44 PgC/yr, 1998-2018, coastal ocean area of 77 million km2). We attribute most of this model-product difference to the seasonality in sea surface CO2 partial pressure at mid- and high-latitudes, where models simulate stronger winter CO2 uptake. The global coastal ocean is a major source of N2O (+0.70 PgCO2-e /yr in observational product and +0.54 PgCO2-e /yr in model median) and of CH4 (+0.21 PgCO2-e /yr in observational product), which offsets a substantial proportion of the net radiative effect of coastal \co uptake (35-58% in CO2-equivalents). Data products and models need improvement to better resolve the spatio-temporal variability and long term trends in CO2, N2O and CH4 in the global coastal ocean.

The Southern California Bight (SCB), an eastern boundary upwelling system, is impacted by global warming, acidification and deoxygetation, and receives anthropogenic nutrients from a coastal population of 20 million people. We describe the configuration, forcing, and validation of a realistic, submesoscale resolving ocean model as a tool to investigate coastal eutrophication. This modeling system represents an important achievement because it strikes a balance of capturing the forcing by U.S. Pacific Coast-wide phenomena, while representing the bathymetric features and submesoscale circulation that affect the vertical and horizontal transport of nutrients from natural and human sources. Moreover, the model allows to run simulations at timescales that approach the interannual frequencies of ocean variability, making the grand challenge of disentangling natural variability, climate change, and local anthropogenic forcing a tractable task in the near-term. The model simulation is evaluated against a broad suite of observational data throughout the SCB, showing realistic depiction of mean state and its variability with remote sensing and in situ physical-biogeochemical measurements of state variables and biogeochemical rates. The simulation reproduces the main structure of the seasonal upwelling front, the mean current patterns, the dispersion of plumes, as well as their seasonal variability. It reproduces the mean distributions of key biogeochemical and ecosystem properties. Biogeochemical rates reproduced by the model, such as primary productivity and nitrification, are also consistent with measured rates. Results of this validation exercise demonstrate the utility of fine-scale resolution modeling in support of management decisions on local anthropogenic nutrient discharges to coastal zones.

Oceanic emissions of nitrous oxide (N2O) account for roughly one-third of all natural sources to the atmosphere. Hot-spots of N2O outgassing occur over oxygen minimum zones (OMZs), where the presence of steep oxygen gradients surrounding anoxic waters leads to enhanced N2O production from both nitrification and denitrification. However, the relative contributions from these pathways to N2O production and outgassing in these regions remains poorly constrained, in part due to shared intermediary nitrogen tracers, and the tight coupling of denitrification sources and sinks. To shed light on this problem, we embed a new, mechanistic model of the OMZ nitrogen cycle within a three-dimensional eddy-resolving physical-biogeochemical model of the ETSP, tracking contributions from remote advection, atmospheric exchange, and local nitrification and denitrification. Our results indicate that net N2O production from denitrification is approximately one order of magnitude greater than nitrification within the ETSP OMZ. However, only ~30% of denitrification-derived N2O production ultimately outgasses to the atmosphere in this region (contributing ~34% of the air-sea N2O flux on an annual basis), while the remaining is exported out of the domain. Instead, remotely-produced N2O advected into the OMZ region accounts for roughly half (~56%) of the total N2O outgassing, with smaller contributions from nitrification (~7%). Our results suggests that, together with enhanced production by denitrification, upwelling of remotely-derived N2O (likely produced via nitrification in the oxygenated ocean) contributes the most to N2O outgassing over the ETSP OMZ.

Continental margins play an essential role in global ocean biogeochemistry and the carbon cycle; however, global assessments of this role remain highly uncertain. This uncertainty arises from large variability over a broad range of temporal and spatial scales of the processes that characterize these environments. High-resolution simulations with ocean biogeochemical models have emerged as essential tools to advance biogeochemical assessments at regional scales. Here, we examine the processes and balances for carbon, oxygen, and nitrogen cycles along the U.S. West Coast in an 11-year hindcast simulation with a submesoscale-permitting oceanic circulation-biogeochemical model. We highlight the importance of biogeochemical cycles on the continental shelf, and their connection to the broader regional context encompassing the California Current System. On the shelf, coastal and wind stress curl upwelling drive a vigorous overturning circulation that supports biogeochemical rates and fluxes that are approximately twice as large as offshore. Exchanges with the proximate sediments, submesoscale shelf currents, bottom boundary layer transport, and intensified cross-shelf export of shelf-produced materials impact coastal and open-ocean balances. While regional variability prevents extrapolation of our results to global margins, our approach provides a powerful tool to identify the dominant dynamics in different shelf setting and quantify their large-scale consequences.

Export of sinking particles from the surface ocean is critical for carbon sequestration and to provide energy to the deep biosphere. The magnitude and spatial patterns of this export have been estimated in the past by \emph{in situ} particle flux observations, satellite-based algorithms, and ocean biogeochemical models; however, these estimates remain uncertain. Here, we use a recent machine learning reconstruction of global ocean particle size distributions from Underwater Vision Profiler 5 (UVP5) measurements to estimate carbon fluxes by sinking particles (35 $\mu$m - 5 mm equivalent spherical diameter) from the surface ocean. We combine global maps of particle size distribution properties with empirical relationships constrained against \emph{in situ} flux observations to calculate particulate carbon export from the euphotic zone and wintertime mixed layer depths. The new flux reconstructions suggest a less variable seasonal cycle in the tropical ocean, and a more persistent export in the Southern Ocean than previously recognized. Smaller particles (less than 420 $\mu$m) contribute most of the flux globally, while larger particles become more important at high latitudes and in tropical upwelling regions. Export from the wintertime mixed layer globally exceeds that from the euphotic zone, suggesting shallow particle recycling and net heterotrophy in the deep euphotic zone. These estimates open the way to fully three-dimensional global reconstructions of particle fluxes in the ocean, supported by the growing database of \emph{in situ} optical observations.