While the physical dimensions of climate change are now routinely assessed through multimodel intercomparisons, projected impacts on the global ocean ecosystem generally rely on individual models with a specific set of assumptions. To address these single-model limitations, we present standardized ensemble projections from six global marine ecosystem models forced with two Earth system models and four emis- sion scenarios with and without fishing. We derive average biomass trends and associated uncertainties across the marine food web. Without fishing, mean global animal biomass decreased by 5% (±4% SD) under low emissions and 17% (±11% SD) under high emissions by 2100, with an average 5% decline for every 1 °C of warming. Projected biomass declines were primarily driven by increasing temperature and decreasing primary production, and were more pronounced at higher trophic levels, a process known as trophic amplification. Fishing did not substantially alter the effects of climate change. Considerable regional variation featured strong biomass increases at high latitudes and de- creases at middle to low latitudes, with good model agreement on the direction of change but variable magnitude. Uncertainties due to vari- ations in marine ecosystem and Earth system models were similar. Ensemble projections performed well compared with empirical data, emphasizing the benefits of multimodel inference to project future outcomes. Our results indicate that global ocean animal biomass con- sistently declines with climate change, and that these impacts are am- plified at higher trophic levels. Next steps for model development include dynamic scenarios of fishing, cumulative human impacts, and the effects of management measures on future ocean biomass trends.

Aim: Climate change will reshape marine ecosystems over the 21st century through diverse and complex mechanisms that are difficult to assess quantitatively. Here, we characterize expectations for how marine community biomass will respond to the energetic consequences of changes in primary production and temperature‐depend‐ ent metabolic rates, under a range of fishing/conservation scenarios.

Location: Global ocean.

Time period: 1950–2100.

Major taxa studied: Commercially harvested marine ectotherms (‘fish’).

Methods: We use a size‐structured macroecological model of the marine ecosystem, coupled with a catch model that allows for calibration with global historical data and simulation of fishing. We examine the four energetic mechanisms that, within the model framework, determine the community response to climate change: net pri‐ mary production, phytoplankton cell size, and the temperature dependencies of growth and natural mortality.

Results: Climate change decreases the modelled global fish community biomass by as much as 30% by 2100. This results from a diminished energy supply to upper trophic levels as photosynthesis becomes more nutrient limited and phytoplankton cells shrink, and from a temperature‐driven increase of natural mortality that, together, overwhelm the effect of accelerated somatic growth rates. Ocean circulation changes drive regional variations of primary production, producing patterns of winners and los‐ ers that largely compensate each other when averaged globally, whereas decreasing phytoplankton size drives weaker but more uniformly negative changes. The climate impacts are similar across the range of conservation scenarios but are slightly amplified in the strong conservation scenarios owing to the greater role of natural mortality. Main conclusions: The spatial pattern of climate impacts is mostly determined by changes in primary production. The overall decline of community biomass is attributed to a temperature‐driven increase of natural mortality, alongside an overall decrease in phytoplankton size, despite faster somatic growth. Our results highlight the importance of the competition between accelerated growth and mortality in a warming ocean.

Horizontal transport at the boundaries of the subtropical gyres plays a crucial role in providing the nutrients that fuel gyre primary productivity, the heat that helps restratify the surface mixed layer, and the dissolved inorganic carbon that influences air‐sea carbon exchange. Mesoscale eddies may be an important component of these horizontal transports; however, previous studies have not quantified the horizontal tracer transport due to eddies across the subtropical gyre boundaries. Here we assess the physical mechanisms that control the horizontal transport of mass, heat, nutrients, and carbon across the North Pacific and North Atlantic subtropical gyre boundaries using the eddy‐rich ocean component of a climate model (Geophysical Fluid Dynamics Laboratory Climate Model version 2.6) coupled to a simple biogeochemical model (mini‐BLING). Our results suggest that horizontal transport across the gyre boundaries supplies a substantial amount of mass and tracers to the ventilated layer of both Northern Hemisphere subtropical gyres, with the Kuroshio and Gulf Stream acting as main exchange gateways. Mass, heat, and dissolved inorganic carbon supply is principally driven by the time‐mean circulation, whereas nutrient transport differs markedly from the other tracers, as nutrients are mainly supplied to both subtropical gyres by downgradient eddy mixing across gyre boundaries. A budget analysis further reveals that the horizontal nutrient transport, combining the roles of both mean and eddy components, is responsible for more than three quarters of the total nutrient supply into the subtropical gyres, surpassing a recent estimate based on a coarse‐resolution model and thus further highlighting the importance of horizontal nutrient transport.

In ocean waters, anaerobic microbial respiration should be confined to the anoxic waters found in coastal regions and tropical oxygen minimum zones, where it is energetically favourable. However, recent molecular and geochemical evidence has pointed to a much broader distribution of denitrifying and sulfate-reducing microbes. Anaerobic metabolisms are thought to thrive in microenvironments that develop inside sinking organic aggregates, but the global distribution and geochemical significance of these microenvironments is poorly understood. Here, we develop a new size-resolved particle model to predict anaerobic respiration from aggregate properties and seawater chemistry. Constrained by observations of the size spectrum of sinking particles, the model predicts that denitrification and sulfate reduction can be sustained throughout vast, hypoxic expanses of the ocean, and could explain the trace metal enrichment observed in particles due to sulfide precipitation. Globally, the expansion of the anaerobic niche due to particle microenvironments doubles the rate of water column denitrification compared with estimates based on anoxic zones alone, and changes the sensitivity of the marine nitrogen cycle to deoxygenation in a warming climate.

Subsurface coherent eddies are well-known features of ocean circulation, but the sparsity of observations prevents an assessment of their importance for biogeochemistry. Here, we use a global eddying (0.1° ) ocean-biogeochemical model to carry out a census of subsurface coherent eddies originating from eastern boundary upwelling systems (EBUS), and quantify their biogeochemical effects as they propagate westward into the subtropical gyres. While most eddies exist for a few months, moving over distances of 100s of km, a small fraction (< 5%) of long-lived eddies propagates over distances greater than 1000km, carrying the oxygen-poor and nutrient-rich signature of EBUS into the gyre interiors. In the Pacific, transport by subsurface coherent eddies accounts for roughly 10% of the offshore transport of oxygen and nutrients in pycnocline waters. This "leakage" of subsurface waters can be a significant fraction of the transport by nutrient-rich poleward undercurrents, and may contribute to the well-known reduction of productivity by eddies in EBUS. Furthermore, at the density layer of their cores, eddies decrease climatological oxygen locally by close to 10%, thereby expanding oxygen minimum zones. Finally, eddies represent low-oxygen extreme events in otherwise oxygenated waters, increasing the area of hypoxic waters by several percent and producing dramatic short-term changes that may play an important ecological role. Capturing these non-local effects in global climate models, which typically include non-eddying oceans, would require dedicated parameterizations.

The global wild marine fish harvest increased fourfold between 1950 and a peak value near the end of the 20th century, reflecting interactions between anthropogenic and ecological forces. Here, we examine these interactions in a bio-energetically constrained, spatially and temporally resolved model of global fisheries. We conduct historical hindcasts with the model, which suggest that technological progress can explain most of the 20th century increase of fish harvest. In contrast, projections extending this rate of technological progress into the future under open access suggest a long-term decrease in harvest due to over-fishing. Climate change is predicted to gradually decrease the global fish production capacity, though our model suggests that this is of secondary importance to social and economic factors. Our study represents a novel way to integrate human-ecological interactions within a single model framework for long-term simulations.

Abstract Hydrothermal vents along the ocean's tectonic ridge systems inject superheated water and large amounts of dissolved metals that impact the deep ocean circulation and the oceanic cycling of trace metals. The hydrothermal fluid contains dissolved mantle helium that is enriched in 3He relative to the atmosphere, providing an isotopic tracer of the ocean's deep circulation and a marker of hydrothermal sources. This work investigates the potential for the 3He/4He isotope ratio to constrain the ocean's mantle 3He source and to provide constraints on the ocean's deep circulation. We use an ensemble of 11 data-assimilated steady-state ocean circulation models and a mantle helium source based on geographically varying sea-floor spreading rates. The global source distribution is partitioned into 6 regions, and the vertical profile and source amplitude of each region are varied independently to determine the optimal 3He source distribution that minimizes the mismatch between modeled and observed δ3He. In this way, we are able to fit the observed δ3He distribution to within a relative error of ∼15%, with a global 3He source that ranges from 640 to 850 mol yr−1, depending on circulation. The fit captures the vertical and interbasin gradients of the δ3He distribution very well and reproduces its jet-sheared saddle point in the deep equatorial Pacific. This demonstrates that the data-assimilated models have much greater fidelity to the deep ocean circulation than other coarse-resolution ocean models. Nonetheless, the modelled δ3He distributions still display some systematic biases, especially in the deep North Pacific where δ3He is overpredicted by our models, and in the southeastern tropical Pacific, where observed westward-spreading δ3He plumes are not well captured. Sources inferred by the data-assimilated transport with and without isopycnally aligned eddy diffusivity differ widely in the Southern Ocean, in spite of the ability to match the observed distributions of \CFCs\ and radiocarbon for either eddy parameterization.

Environmental change and the exploitation of marine resources have had profound impacts on marine communities, with potential implications for ocean biogeochemistry and food security. In order to study such global-scale problems, it is helpful to have computationally efficient numerical models that predict the first-order features of fish biomass production as a function of the environment, based on empirical and mechanistic understandings of marine ecosystems. Here we describe the ecological module of the BiOeconomic mArine Trophic Size-spectrum (BOATS) model, which takes an Earth-system approach to modelling fish biomass at the global scale. The ecological model is designed to be used on an Earth-system model grid, and determines size spectra of fish biomass by explicitly resolving life history as a function of local temperature and net primary production. Biomass production is limited by the availability of photosynthetic energy to upper trophic levels, following empirical trophic efficiency scalings, and by well-established empirical temperature-dependent growth rates. Natural mortality is calculated using an empirical size-based relationship, while reproduction and recruitment depend on both the food availability to larvae from net primary production and the production of eggs by mature adult fish. We describe predicted biomass spectra and compare them to observations, and conduct a sensitivity study to determine how they change as a function of net primary production and temperature. The model relies on a limited number of parameters compared to similar modelling efforts, while retaining reasonably realistic representations of biological and ecological processes, and is computationally efficient, allowing extensive parameter-space analyses even when implemented globally. As such, it enables the exploration of the linkages between ocean biogeochemistry, climate, and upper trophic levels at the global scale, as well as a representation of fish biomass for idealized studies of fisheries.