It is widely recognized that the stoichiometry of nutrient elements in phytoplankton varies within the ocean. However, there are many conflicting mechanistic explanations for this variability, and it is often ignored in global biogeochemical models and carbon cycle simulations. Here we show that globally distributed particulate P:C varies as a linear function of ambient phosphate concentrations, whereas the N:C varies with ambient nitrate concentrations, but only when nitrate is most scarce. This observation is consistent with the adjustment of the phytoplankton community to local nutrient availability, with greater flexibility of phytoplankton P:C because P is a less abundant cellular component than N. This simple relationship is shown to predict the large-scale, long-term average composition of surface particles throughout large parts of the ocean remarkably well. The relationship implies that most of the observed variation in N:P actually arises from a greater plasticity in the cellular P:C content, relative to N:C, such that as overall macronutrient concentrations decrease, N:P rises. Although other mechanisms are certainly also relevant, this simple relationship can be applied as a first-order basis for predicting organic matter stoichiometry in large-scale biogeochemical models, as illustrated using a simple box model. The results show that including variable P:C makes atmospheric CO2 more sensitive to changes in low latitude export and ocean circulation than a fixed-stoichiometry model. In addition, variable P:C weakens the relationship between preformed phosphate and atmospheric CO2 while implying a more important role for the nitrogen cycle.

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.

The burial of organic carbon in marine sediments removes carbon dioxide from the ocean-atmosphere pool, provides energy to the deep biosphere, and on geological timescales drives the oxygenation of the atmosphere. Here we quantify natural variations in the burial of organic carbon in deep-sea sediments over the last glacial cycle. Using a new data compilation of hundreds of sediment cores, we show that the accumulation rate of organic carbon in the deep sea was consistently higher (50%) during glacial maxima than during interglacials. The spatial pattern and temporal progression of the changes suggest that enhanced nutrient supply to parts of the surface ocean contributed to the glacial burial pulses, with likely additional contributions from more efficient transfer of organic matter to the deep sea and better preservation of organic matter due to reduced oxygen exposure. These results demonstrate a pronounced climate sensitivity for this global carbon cycle sink.

The biomass and biogeochemical roles of fish in the ocean are ecologically important but poorly known. Here, we use a data-constrained marine ecosystem model to provide a first-order estimate of the historical reduction of fish biomass due to fishing and the associated change in biogeochemical cycling rates. The pre-exploitation global biomass of exploited fish (10 g to 100 kg) was 3.3 ± 0.5 Gt, cycling roughly 2% of global primary production (9.4 ± 1.6 Gt year⁻¹) and producing 10% of surface biological export. Particulate organic matter produced by exploited fish drove roughly 10% of the oxygen consumption and biological carbon storage at depth. By the 1990s, biomass and cycling rates had been reduced by nearly half, suggesting that the biogeochemical impact of fisheries has been comparable to that of anthropogenic climate change. Our results highlight the importance of developing a better mechanistic understanding of how fish alter ocean biogeochemistry.

Global observations show that the ocean lost approximately 2% of its oxygen inventory over the last five decades 1-3, with important implications for marine ecosystems 4, 5. The rate of change varies with northwest Atlantic coastal waters showing a long-term drop 6, 7 that vastly outpaces the global and North Atlantic basin mean deoxygenation rates 5, 8. However, past work has been unable to resolve mechanisms of large-scale climate forcing from local processes. Here, we use hydrographic evidence to show a Labrador Current retreat is playing a key role in the deoxygenation on the northwest Atlantic shelf. A high-resolution global coupled climate-biogeochemistry model 9 reproduces the observed decline of saturation oxygen concentrations in the region, driven by a retreat of the equatorward-flowing Labrador Current and an associated shift toward more oxygen-poor subtropical waters on the shelf. The dynamical changes underlying the shift in shelf water properties are correlated with a slowdown in the simulated Atlantic Meridional Overturning Circulation 10. Our results provide strong evidence that a major, centennial-scale change of the Labrador Current is underway, and highlight the potential for ocean dynamics to impact coastal deoxygenation over the coming century.

The carbon isotope composition (δ13C) of seawater provides valuable insight on ocean circulation, air-sea exchange, the biological pump, and the global carbon cycle and is reflected by the δ13C of foraminifera tests. Here more than 1700 δ13C observations of the benthic foraminifera genus Cibicides from late Holocene sediments (δ13CCibnat) are compiled and compared with newly updated estimates of the natural (preindustrial) water column δ13C of dissolved inorganic carbon (δ13CDICnat) as part of the international Ocean Circulation and Carbon Cycling (OC3) project. Using selection criteria based on the spatial distance between samples, we find high correlation between δ13CCibnat and δ13CDICnat, confirming earlier work. Regression analyses indicate significant carbonate ion (−2.6 ± 0.4) × 10−3‰/(μmol kg−1) [CO32−] and pressure (−4.9 ± 1.7) × 10−5‰ m−1 (depth) effects, which we use to propose a new global calibration for predicting δ13CDICnat from δ13CCibnat. This calibration is shown to remove some systematic regional biases and decrease errors compared with the one-to-one relationship (δ13CDICnat = δ13CCibnat). However, these effects and the error reductions are relatively small, which suggests that most conclusions from previous studies using a one-to-one relationship remain robust. The remaining standard error of the regression is generally σ ≅ 0.25‰, with larger values found in the southeast Atlantic and Antarctic (σ ≅ 0.4‰) and for species other than Cibicides wuellerstorfi. Discussion of species effects and possible sources of the remaining errors may aid future attempts to improve the use of the benthic δ13C record.

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 emission 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 decreases at middle to low latitudes, with good model agreement on the direction of change but variable magnitude. Uncertainties due to variations 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 consistently declines with climate change, and that these impacts are amplified 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.

Projections of climate change impacts on marine ecosystems have revealed long-term declines in global marine animal biomass and unevenly distributed impacts on fisheries. Here we apply an enhanced suite of global marine ecosystem models from the Fisheries and Marine Ecosystem Model Intercomparison Project (Fish-MIP), forced by new-generation Earth system model outputs from Phase 6 of the Coupled Model Intercomparison Project (CMIP6), to provide insights into how projected climate change will affect future ocean ecosystems. Compared with the previous generation CMIP5-forced Fish-MIP ensemble, the new ensemble ecosystem simulations show a greater decline in mean global ocean animal biomass under both strong-mitigation and high-emissions scenarios due to elevated warming, despite greater uncertainty in net primary production in the high-emissions scenario. Regional shifts in the direction of biomass changes highlight the continued and urgent need to reduce uncertainty in the projected responses of marine ecosystems to climate change to help support adaptation planning.