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.

No abstract

Oxygen minimum zones (OMZs) are major sites of fixed nitrogen removal from the open ocean. However, commonly used gridded data sets such as the World Ocean Atlas (WOA) tend to overestimate the concentration of O2 compared to measurements in grids where O2 falls in the suboxic range (O2 < 2–10 mmol m−3), thereby underestimating the extent of O2 depletion in OMZs. We evaluate the distribution of the OMZs by (1) mapping high‐quality oxygen measurements from the WOCE program, and (2) by applying an empirical correction to the gridded WOA based on in situ observations. The resulting suboxic volumes are a factor 3 larger than in the uncorrected gridded WOA. We combine the new oxygen data sets with estimates of global export and simple models of remineralization to estimate global denitrification and N2O production. We obtain a removal of fixed nitrogen of 70 ± 50 Tg year−1 in the open ocean and 198 ± 64 Tg year−1 in the sediments, and a global N2O production of 6.2 ± 3.2 Tg year−1. Our results (1) reconcile water column denitrification rates based on global oxygen distributions with previous estimates based on nitrogen isotopes, (2) revise existing estimates of sediment denitrification down by 1/3d through the use of spatially explicit fluxes, and (3) provide independent evidence supporting the idea of a historically balanced oceanic nitrogen cycle. These estimates are most sensitive to uncertainties in the global export production, the oxygen threshold for suboxic processes, and the efficiency of particle respiration under suboxic conditions. Ocean deoxygenation, an expected response to anthropogenic climate change, could increase denitrification by 14 Tg year−1 of nitrogen per 1 mmol m−3 of oxygen reduction if uniformly distributed, while leaving N2O production relatively unchanged.

Nitrous oxide (N2O) is a powerful greenhouse gas and a major cause of stratospheric ozone depletion, yet its sources and sinks remain poorly quantified in the oceans. We used isotope tracers to directly measure N2O reduction rates in the eastern tropical North Pacific. Because of incomplete denitrification, N2O cycling rates are an order of magnitude higher than predicted by current models in suboxic regions, and the spatial distribution suggests strong dependence on both organic carbon and dissolved oxygen concentrations. Furthermore, N2O turnover is 20 times higher than the net atmospheric efflux. The rapid rate of this cycling coupled to an expected expansion of suboxic ocean waters implies future increases in N2O emissions.

Abstract. We analyse simulations of the Pacific Ocean oxygen minimum zones (OMZs) from 11 Earth system model contributions to the Coupled Model Intercomparison Project Phase 5, focusing on the mean state and climate change projections. The simulations tend to overestimate the volume of the OMZs, especially in the tropics and Southern Hemisphere. Compared to observations, five models introduce incorrect meridional asymmetries in the distribution of oxygen including larger southern OMZ and weaker northern OMZ, due to interhemispheric biases in intermediate water mass ventilation. Seven models show too deep an extent of the tropical hypoxia compared to observations, stemming from a deficient equatorial ventilation in the upper ocean, combined with too large a biologically driven downward flux of particulate organic carbon at depth, caused by particle export from the euphotic layer that is too high and remineralization in the upper ocean that is too weak. At interannual timescales, the dynamics of oxygen in the eastern tropical Pacific OMZ is dominated by biological consumption and linked to natural variability in the Walker circulation. However, under the climate change scenario RCP8.5, all simulations yield small and discrepant changes in oxygen concentration at mid depths in the tropical Pacific by the end of the 21st century due to an almost perfect compensation between warming-related decrease in oxygen saturation and decrease in biological oxygen utilization. Climate change projections are at odds with recent observations that show decreasing oxygen levels at mid depths in the tropical Pacific. Out of the OMZs, all the CMIP5 models predict a decrease of oxygen over most of the surface and deep ocean at low latitudes and over all depths at high latitudes due to an overall slow-down of ventilation and increased temperature.

No abstract

[1] Diel vertical migration (DVM) of zooplankton and micronekton is widespread in the ocean and forms a fundamental component of the biological pump, but is generally overlooked in global models of the Earth system. We develop a parameterization of DVM in the ocean and integrate it with a size-structured NPZD model. We assess the model's ability to recreate ecosystem and DVM patterns at three well-observed Pacific sites, ALOHA, K2, and EQPAC, and use it to estimate the impact of DVM on marine ecosystems and biogeochemical dynamics. Our model includes the following: (1) a representation of migration dynamics in response to food availability and light intensity; (2) a representation of the digestive and metabolic processes that decouple zooplankton feeding from excretion, egestion, and respiration; and (3) a light-dependent parameterization of visual predation on zooplankton. The model captures the first-order patterns in plankton biomass and productivity across the biomes, including the biomass of migrating organisms. We estimate that realistic migratory populations sustain active fluxes to the mesopelagic zone equivalent to between 15% and 40% of the particle export and contribute up to half of the total respiration within the layers affected by migration. The localized active transport has important consequences for the cycling of oxygen, nutrients, and carbon. We highlight the importance of decoupling zooplankton feeding and respiration and excretion with depth for capturing the impact of migration on the redistribution of carbon and nutrients in the upper ocean.

The distribution of radiocarbon (14C) in the ocean and atmosphere has fluctuated on time scales ranging from seasons to millennia. It is thought that these fluctuations partly reflect variability in the climate system, offering a rich potential source of information to help understand mechanisms of past climate change. Here, a long simulation with a new, coupled model is used to explore the mechanisms that redistribute 14C within the earth system on interannual to centennial time scales. The model, the Geophysical Fluid Dynamics Laboratory Climate Model version 2 (GFDL CM2) with Modular Ocean Model version 4p1(MOM4p1) at coarse-resolution (CM2Mc), is a lower-resolution version of the Geophysical Fluid Dynamics Laboratory’s CM2M model, uses no flux adjustments, and is run here with a simple prognostic ocean biogeochemistry model including 14C. The atmospheric 14C and radiative boundary conditions are held constant so that the oceanic distribution of 14C is only a function of internal climate variability. The simulation displays previously described relationships between tropical sea surface 14C and the model equivalents of the El Niño–Southern Oscillation and Indonesian Throughflow. Sea surface 14C variability also arises from fluctuations in the circulations of the subarctic Pacific and Southern Ocean, including North Pacific decadal variability and episodic ventilation events in the Weddell Sea that are reminiscent of the Weddell Polynya of 1974–76. Interannual variability in the air–sea balance of 14C is dominated by exchange within the belt of intense “Southern Westerly” winds, rather than at the convective locations where the surface 14C is most variable. Despite significant interannual variability, the simulated impact on air–sea exchange is an order of magnitude smaller than the recorded atmospheric 14C variability of the past millennium. This result partly reflects the importance of variability in the production rate of 14C in determining atmospheric 14C but may also reflect an underestimate of natural climate variability, particularly in the Southern Westerly winds.