The effect of idealized wind-driven circulation changes in the Southern Ocean on atmospheric CO 2 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 2 : 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 2 by *20 latm; 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 2 . Increased residual overturning deepens the global pycnocline, warming the upper ocean and decreasing the saturated carbon reservoir.

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

Quantifying variability in the ocean carbon sink remains problematic due to sparse observations and spatiotemporal variability in surface ocean pCO 2. 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 2 and air-sea CO 2 flux reconstructions in most biomes, particularly in the subtropical and equatorial regions. The largest differences between CO 2 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 2 sink (2.47 ± 0.50 Pg C year −1) 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. Plain Language Summary Data-driven estimates of how much carbon dioxide the ocean is absorbing (the so-called "ocean carbon sink") have improved substantially in recent years. However, computational ocean models that include biogeochemistry continue to play a critical role as they allow us to isolate and understand the individual processes that control ocean carbon sequestration. The ideal scenario is a combination of the above two methods, where data are ingested and then used to improve a model's fit to the observed ocean, also known as, data assimilation. While the physical oceanographic community has made great progress in developing data assimilation systems, for example, the Estimating the Circulation and Climate of the Ocean (ECCO) consortium, the biogeochemical community has generally lagged behind. The ECCO-Darwin model presented in this paper represents an important technological step forward as it is the first global ocean biogeochemistry model that (1) ingests both physical and biogeochemical observations into the model in a realistic manner and (2) considers how the nature of the ocean carbon sink has changed over multiple decades. As the ECCO ocean circulation

a rate of typically 1-1.5 parts per million/10 6 m 3 s −1 by enhancing the communication between the atmosphere and deep ocean. This metric can be used to interpret the longterm effect of Southern Ocean dynamics on the natural carbon cycle and atmospheric CO 2 , alongside other metrics, such as involving the proportion of preformed nutrients and the extent of sea ice cover.

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.