Quantifying variability in the ocean carbon sink remains problematic due to sparse observations and spatiotemporal variability in surface ocean pCO₂. To address this challenge, we have updated and improved ECCO‐Darwin, a global ocean biogeochemistry model that assimilates both physical and biogeochemical observations. The model consists of an adjoint‐based ocean circulation estimate from the Estimating the Circulation and Climate of the Ocean (ECCO) consortium and an ecosystem model developed by the Massachusetts Institute of Technology Darwin Project. In addition to the data‐constrained ECCO physics, a Green's function approach is used to optimize the biogeochemistry by adjusting initial conditions and six biogeochemical parameters. Over seasonal to multidecadal timescales (1995–2017), ECCO‐Darwin exhibits broad‐scale consistency with observed surface ocean pCO₂ and air‐sea CO₂ flux reconstructions in most biomes, particularly in the subtropical and equatorial regions. The largest differences between CO₂ uptake occur in subpolar seasonally stratified biomes, where ECCO‐Darwin results in stronger winter uptake. Compared to the Global Carbon Project OBMs, ECCO‐Darwin has a time‐mean global ocean CO₂ sink (2.47 ± 0.50 Pg C year⁻¹) and interannual variability that are more consistent with interpolation‐based products. Compared to interpolation‐based methods, ECCO‐Darwin is less sensitive to sparse and irregularly sampled observations. Thus, ECCO‐Darwin provides a basis for identifying and predicting the consequences of natural and anthropogenic perturbations to the ocean carbon cycle, as well as the climate‐related sensitivity of marine ecosystems. Our study further highlights the importance of physically consistent, property‐conserving reconstructions, as are provided by ECCO, for ocean biogeochemistry studies.

The inventory and variability of oceanic dissolved inorganic carbon (DIC) is driven by the interplay of physical, chemical, and biological processes. Quantifying the spatiotemporal variability of these drivers is crucial for a mechanistic understanding of the ocean carbon sink and its future trajectory. Here, we use the Estimating the Circulation and Climate of the Ocean‐Darwin ocean biogeochemistry state estimate to generate a global‐ocean, data‐constrained DIC budget and investigate how spatial and seasonal‐to‐interannual variability in three‐dimensional circulation, air‐sea CO ₂ flux, and biological processes have modulated the ocean sink for 1995–2018. Our results demonstrate substantial compensation between budget terms, resulting in distinct upper‐ocean carbon regimes. For example, boundary current regions have strong contributions from vertical diffusion while equatorial regions exhibit compensation between upwelling and biological processes. When integrated across the full ocean depth, the 24‐year DIC mass increase of 64 Pg C (2.7 Pg C year ⁻¹ ) primarily tracks the anthropogenic CO ₂ growth rate, with biological processes providing a small contribution of 2% (1.4 Pg C). In the upper 100 m, which stores roughly 13% (8.1 Pg C) of the global increase, we find that circulation provides the largest DIC gain (6.3 Pg C year ⁻¹ ) and biological processes are the largest loss (8.6 Pg C year ⁻¹ ). Interannual variability is dominated by vertical advection in equatorial regions, with the 1997–1998 El Niño‐Southern Oscillation causing the largest year‐to‐year change in upper‐ocean DIC (2.1 Pg C). Our results provide a novel, data‐constrained framework for an improved mechanistic understanding of natural and anthropogenic perturbations to the ocean sink.

The effect of idealized wind-driven circulation changes in the Southern Ocean on atmospheric CO₂ and the ocean carbon inventory is investigated using a suite of coarse-resolution, global coupled ocean circulation and biogeochemistry experiments with parameterized eddy activity and only modest changes in surface buoyancy forcing, each experiment integrated for 5,000 years. A positive correlation is obtained between the meridional overturning or residual circulation in the Southern Ocean and atmospheric CO₂: stronger or northward-shifted westerly winds in the Southern Hemisphere result in increased residual circulation, greater upwelling of carbon-rich deep waters and oceanic outgassing, which increases atmospheric pCO₂ by ∼20 μatm; weaker or southward-shifted winds lead to the opposing result. The ocean carbon inventory in our model varies through contrasting changes in the saturated, disequilibrium and biogenic (soft-tissue and carbonate) reservoirs, each varying by O(10–100) PgC, all of which contribute to the net anomaly in atmospheric CO₂. Increased residual overturning deepens the global pycnocline, warming the upper ocean and decreasing the saturated carbon reservoir. Increased upwelling of carbon- and nutrient-rich deep waters and inefficient biological activity results in subduction of unutilized nutrients into the ocean interior, decreasing the biogenic carbon reservoir of intermediate and mode waters ventilating the Northern Hemisphere, and making the disequilibrium carbon reservoir more positive in the mode waters due to the reduced residence time at the surface. Wind-induced changes in the model carbon inventory are dominated by the response of the global pycnocline, although there is an additional abyssal response when the peak westerly winds change their latitude, altering their proximity to Drake Passage and changing the depth extent of the southward return flow of the overturning: a northward shift of the westerly winds isolates dense isopycnals, allowing biogenic carbon to accumulate in the deep ocean of the Southern Hemisphere, while a southward shift shoals dense isopycnals that outcrop in the Southern Ocean and reduces the biogenic carbon store in the deep ocean.

The expansive gyres of the subtropical ocean account for a significant fraction of global organic carbon export from the upper ocean. In the gyre interior, vertical mixing and the heaving of nutrient-rich waters into the euphotic layer sustain local productivity, in turn depleting the layers below. However, the nutrient pathways by which these subeuphotic layers are themselves replenished remain unclear. Using a global, eddy-permitting simulation of ocean physics and biogeochemistry, we quantify nutrient resupply mechanisms along and across density surfaces, including the contribution of eddy-scale motions that are challenging to observe. We find that mesoscale eddies (10 to 100 km) flux nutrients from the shallow flanks of the gyre into the recirculating interior, through time-varying motions along density surfaces. The subeuphotic layers are ultimately replenished in approximately equal contributions by this mesoscale eddy transport and the remineralization of sinking particles. The mesoscale eddy resupply is most important in the lower thermocline for the whole subtropical region but is dominant at all depths within the gyre interior. Subtropical gyre productivity may therefore be sustained by a nutrient relay, where the lateral transport resupplies nutrients to the thermocline and allows vertical exchanges to maintain surface biological production and carbon export.

Mechanistic description of the transition from aerobic to anaerobic metabolism is necessary for diagnostic and predictive modeling of fixed nitrogen loss in anoxic marine zones (AMZs). In a metabolic model where diverse oxygen- and nitrogen-cycling microbial metabolisms are described by underlying redox chemical reactions, we predict a transition from strictly aerobic to predominantly anaerobic regimes as the outcome of ecological interactions along an oxygen gradient, obviating the need for prescribed critical oxygen concentrations. Competing aerobic and anaerobic metabolisms can coexist in anoxic conditions whether these metabolisms represent obligate or facultative populations. In the coexistence regime, relative rates of aerobic and anaerobic activity are determined by the ratio of oxygen to electron donor supply. The model simulates key characteristics of AMZs, such as the accumulation of nitrite and the sustainability of anammox at higher oxygen concentrations than denitrification, and articulates how microbial biomass concentrations relate to associated water column transformation rates as a function of redox stoichiometry and energetics. Incorporating the metabolic model into an idealized two-dimensional ocean circulation results in a simulated AMZ, in which a secondary chlorophyll maximum emerges from oxygen-limited grazing, and where vertical mixing and dispersal in the oxycline also contribute to metabolic co-occurrence. The modeling approach is mechanistic yet computationally economical and suitable for global change applications.

Iron is the limiting factor for biological production over a large fraction of the surface ocean because free iron is rapidly scavenged or precipitated under aerobic conditions. Standing stocks of dissolved iron are maintained by association with organic molecules (ligands) produced by biological processes. We hypothesize a positive feedback between iron cycling, microbial activity, and ligand abundance: External iron input fuels microbial production, creating organic ligands that support more iron in seawater, leading to further macronutrient consumption until other microbial requirements such as macronutrients or light become limiting, and additional iron no longer increases productivity. This feedback emerges in numerical simulations of the coupled marine cycles of macronutrients and iron that resolve the dynamic microbial production and loss of iron-chelating ligands. The model solutions resemble modern nutrient distributions only over a finite range of prescribed ligand source/sink ratios where the model ocean is driven to global-scale colimitation by micronutrients and macronutrients and global production is maximized. We hypothesize that a global-scale selection for microbial ligand cycling may have occurred to maintain “just enough” iron in the ocean.