
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39231975
52200
10.1038/s41467-024-52200-0
Article
Transient overturning changes cause an upper-ocean nutrient decline in a warming climate
http://orcid.org/0000-0002-6932-5589
Sun Shantong stsun@qnlm.ac

123
http://orcid.org/0000-0003-0322-4811
Thompson Andrew F. 3
Yu Jimin 145
http://orcid.org/0000-0002-4694-5531
Wu Lixin 16
1 https://ror.org/041w4c980 Laoshan Laboratory, Qingdao, China
2 https://ror.org/05g3dte14 grid.255986.5 0000 0004 0472 0419 Earth, Ocean and Atmospheric Science, Florida State University, Tallahassee, FL USA
3 https://ror.org/05dxps055 grid.20861.3d 0000 0001 0706 8890 Environmental Science and Engineering, California Institute of Technology, Pasadena, CA USA
4 grid.1001.0 0000 0001 2180 7477 Research School of Earth Sciences, Australian National University, Canberra, ACT Australia
5 grid.9227.e 0000000119573309 SKLLQG, Institute of Earth Environment, Chinese Academy of Sciences, Xi’an, China
6 https://ror.org/04rdtx186 grid.4422.0 0000 0001 2152 3263 Frontiers Science Center for Deep Ocean Multispheres and Earth System/Key Laboratory of Physical Oceanography, Ocean University of China, Qingdao, China
4 9 2024
4 9 2024
2024
15 772712 4 2024
29 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Models and proxy data suggest multi-centennial nutrient reorganization and biological productivity changes under sustained climate warming. These changes have traditionally been attributed to processes in the Southern Ocean. Here we instead show that transient overturning circulation adjustments, associated with changes in the Atlantic Meridional Overturning Circulation (AMOC), dominate the global nutrient reorganization on centennial timescales. Following an AMOC weakening, a typical feature of a warming climate, a transient overturning circulation develops in the Indo-Pacific basins, characterized by enhanced southward transport in the deep ocean. Coupled with the vertical nutrient structure, these transient overturning changes produce a net transport of nutrients from the Indo-Pacific into the Southern Ocean. Meanwhile, isopycnal surfaces deepen and bring nutrient-depleted waters to greater depths, causing nutrient concentrations to decline in much of the global upper ocean. Given the close link between nutrients and carbon, our findings suggest that transient overturning circulation changes across different basins can critically affect the marine carbon cycle.

Sun et al. showed that transient overturning circulation responses, involving interactions between the Atlantic and Indo-Pacific, dominate global nutrient redistribution and biological productivity changes in a warming climate on centennial time scales.

Subject terms

Carbon cycle
Physical oceanography
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Due to rising atmospheric CO2, the ocean is experiencing changes in physical, chemical, and biological environments with critical impacts on ocean ecosystems1–3. Increases in surface stratification due to global warming are projected to limit surface nutrient supply and decrease biological productivity in low-to-mid latitudes during the 21st century. A global decline in biological productivity is likely to continue beyond the 21st century, but attributing this long-term change to enhanced surface stratification alone is difficult, because upper ocean changes occur on relatively short timescales4,5. Instead, long-term productivity changes are more likely caused by nutrient redistribution between different ocean basins and/or different water depths2,6, with nutrient concentrations decreasing in much of the upper Atlantic and Indo-Pacific basins but increasing in the Southern Ocean. The associated mechanisms for this nutrient redistribution remain debated.

A recent study attributes global nutrient redistribution on centennial timescales to nutrient trapping in the Southern Ocean2, hinging on coupled changes in wind-driven deep water upwelling and biological activity in the Southern Ocean. As the climate warms, Southern Hemisphere westerly winds are predicted to shift poleward and intensify7,8, potentially increasing circumpolar upwelling of nutrient-rich deep waters in the Southern Ocean. The resulting biological productivity enhancement in the Southern Ocean, due to nutrient increases and favorable warm growth conditions, would increase the vertical export of organic matter that transfers nutrients from the surface to the deep ocean, where sinking particles are remineralized. This enhanced biological export reduces nutrient concentrations in the northward subducting waters that sustain low-latitude biological productivity today9. Together, these processes are thought to result in nutrient trapping in the Southern Ocean2. Critically, if these Southern Ocean processes dominate the centennial global nutrient redistribution, the upper-ocean nutrient reduction would be most prominent in the Southern Hemisphere, due to the northward subducting intermediate and mode waters. However, this cannot fully explain the nutrient structure changes simulated by earth system models2, where nutrient concentrations decline more significantly in the Northern Hemisphere. Thus, additional processes are required to explain Southern Ocean nutrient trapping and global nutrient redistribution.

Investigation of nutrient cycling in the ocean is also important for understanding past ocean and climate changes. Multi-centennial nutrient reorganization occurred during the last deglaciation. Paleoclimate records suggest enhanced diatom production around the Antarctic Polar Front during Heinrich events10,11, when the Atlantic Meridional Overturning Circulation (AMOC) was believed to be substantially weakened12,13. The enhanced diatom production was interpreted as evidence for increased Southern Ocean upwelling of nutrient-rich deep waters10,11,14, driven by an intensification and/or southward shift of Southern Hemisphere westerly winds15,16. However, more recent theoretical and modeling studies suggest a much weaker sensitivity of the Southern Ocean upwelling to westerly wind changes due to compensation effects from transient and standing eddies17–20. Thus, the dominant factors affecting the centennial global nutrient redistribution warrant further investigation.

In this study, we highlight the role of transient overturning circulation responses, involving interactions between the Atlantic and the Indo-Pacific basins, in the global redistribution of nutrients and biological productivity reduction in a warming climate on centennial time scales (Fig. 1). Following a weakening of the AMOC, a typical feature of a warming climate21,22, the Indo-Pacific ocean develops an overturning circulation anomaly that opposes the Atlantic changes (black dashed arrows in Fig. 1)23. This inter-ocean response is established via the propagation of planetary waves, associated with a global adiabatic deepening of density surfaces (blue dashed arrows in Fig. 1)23–25. This process is distinct from previous views of Indo-Pacific overturning changes that involve the initiation/enhancement of deep and/or intermediate water formation in the Pacific26–29.Fig. 1 Schematic showing tracer redistributions by transient ocean overturning changes.

Color shadings show climatological distributions of dissolved inorganic silicate (Si) from GLODAPv252 along a 30° W in the Atlantic and b 170° W in the Pacific. Following a slowdown of the Atlantic overturning (black dashed line with arrow in (a) shows the overturning anomaly), the Indo-Pacific develops an opposite overturning circulation anomaly (black dashed line with arrow in (b)). These transient overturning circulation changes redistribute Si between the Southern Ocean and the northern basins, which increases Si in the Southern Ocean (“+”) and reduces Si in the Indo-Pacific (“−”). Associated with these overturning adjustments, isopycnal surfaces deepen (blue dashed arrows) and bring nutrient-depleted surface waters to greater depths, decreasing nutrient concentrations in the upper Atlantic and Indo-Pacific. The red dotted lines represent the boundary between the Southern Ocean and the northern basins.

More specifically, changes in the AMOC modify the density structure in the North Atlantic. The associated stratification changes propagate into the South Atlantic and the Indo-Pacific, leading to geostrophically-balanced, meridional volume transport anomalies that have opposing signs between the Atlantic and Indo-Pacific basins (dashed arrows in Fig. 1). This process affects the global distribution of heat25 and biogeochemical tracers, but has been underappreciated and underexplored in previous studies. Coupled with the vertical structure of nutrient concentrations, the transient overturning response in the Indo-Pacific causes a global reorganization of nutrients between different ocean basins (green arrow in Fig. 1), affecting global biological productivity in a warming climate.

Results

Nutrient and productivity changes in a warming climate

We use outputs from the Community Earth System Model 2 (CESM2; see “CESM2 simulations” in “Methods”), as part of the Coupled Model Intercomparison Project Phase 630. In response to atmospheric CO2 quadrupling (“4×CO2”), the surface ocean stratification increases rapidly over the first 100 years, due to surface warming, and then levels off (Supplementary Fig. 1). This rapid stratification increase, which restricts nutrient entrainment into the euphotic layer (“CESM2 simulations” in “Methods”), suppresses diatom production and organic matter export in the Atlantic and Indo-Pacific basins (Supplementary Figs. 2 and 3).

Primary production continues to decrease in much of the Indo-Pacific and Atlantic basins over the next few centuries (Fig. 2a, b and Supplementary Fig. 2). However, these long-term biogeochemical responses can no longer be explained by surface stratification, which only increases minimally after  ~100 years. Instead, the multi-centennial decline in biological production is related to a global reorganization of nutrients (Fig. 2c–f and Supplementary Figs. 4 and 5). Here we focus on the distribution of dissolved inorganic silicate (Si), which is critical to diatom production and surface organic matter export31. Si decreases in much of the subsurface Indo-Pacific basin, with maximum reductions at around 1 km water depth in the North Pacific, and in the upper 1.5 km of the Atlantic, but Si increases in much of the Southern Ocean and in the deep Atlantic, exhibiting a nutrient redistribution between different ocean basins (Fig. 2).Fig. 2 Nutrients and primary production in CESM2.

a Primary production in the pre-industrial Control run (piControl), averaged over the last 50 years. b Primary production changes after CO2 quadrupling, averaged over years 150–200. c Si concentration along 170° W in the Pacific. d Si concentration changes (ΔSi) after CO2 quadrupling, averaged over years 150–200. e Same as (c) but along 30° W in the Atlantic. f Same as (d) but along 30° W in the Atlantic. The red dotted lines represent the boundary between the Southern Ocean and the northern basins. Source data are provided as a Source Data file.

This nutrient redistribution was previously attributed to wind-driven circumpolar enhancement of deep water upwelling and increases in biological production in the Southern Ocean2, but this mechanism cannot fully explain the structure of nutrient changes in the Pacific Ocean. Specifically, if the global nutrient redistribution is dominated by Southern Ocean processes, the enhanced biological production would lead to the largest nutrient reductions in the cores of intermediate/mode waters in the Southern Hemisphere near 1 km water depth, inconsistent with the simulated substantial nutrient reductions in the North Pacific (Fig. 2d). We note that a long-term decline of biological productivity in the northern basins may reduce subsurface nutrient levels in the North Pacific, but the associated processes would primarily redistribute nutrients vertically along the water column, not between ocean basins (Supplementary Fig. 5).

Transient overturning circulation adjustments and nutrient redistribution

The global nutrient redistribution, occurring on centennial timescales, is better explained by transient adjustments of the global ocean overturning circulation. Following atmospheric CO2 quadrupling, the ocean's overturning circulation adjusts to changes in surface forcing. A prominent feature is a slow-down of the AMOC (Fig. 3a and Supplementary Fig. 6)22. At similar water depths where the AMOC weakens, an overturning circulation anomaly with an opposing sign develops in the Indo-Pacific (Fig. 3a, b), characterized by a northward transport anomaly in the upper 1 km and a southward transport anomaly at greater depths (dashed arrow in Fig. 1b). During the overturning circulation adjustments, the Indo-Pacific overturning response, ΔψP, largely compensates for the AMOC weakening, ΔψA, at 30° S on centennial timescales, i.e., ΔψP ≈ −ψA (Supplementary Fig. 6b). However, nutrient distributions differ between the Atlantic and Indo-Pacific. Nutrients are more enriched in the deep Indo-Pacific and have stronger vertical gradients over the upper 3 km due to weaker ventilation, as compared to the Atlantic (Figs. 1 and 2c, e). As such, when the compensating overturning circulation changes in the Indo-Pacific and the Atlantic basins are coupled with the different vertical nutrient gradients, a net redistribution of nutrients into the Southern Ocean occurs (see “Nutrient transport” in “Methods”).Fig. 3 Overturning responses and nutrient redistribution in CESM2 4×CO2 experiment.

a Changes in the AMOC stream function after CO2 quadrupling (negative for anticlockwise), averaged over years 150–200. b Same as (a) but for the Indo-Pacific overturning. c Si transport anomalies (ΔTA, ΔTP) and their contributions due to overturning circulation changes (ΔTAψ, ΔTPψ) across the southern boundary (30° S) of the Atlantic (subscript “A”) and Indo-Pacific basins (subscript “P”). d Same as (c) but for the net Si transport anomaly between the Southern Ocean and the northern basins. The Si transport anomaly due to Southern Ocean overturning changes, ΔTSψ, is plotted as a pink curve. Negative values in (c and d) represent southward transport anomalies. Source data are provided as a Source Data file.

To demonstrate the above mechanism, we quantify the Si transport anomalies, ΔT, between the Southern Ocean and the northern basins across 30° S in the model (see “Nutrient transport” in “Methods”). We find a southward Si transport anomaly, ΔT<0, into the Southern Ocean (purple line in Fig. 3d), which is driven by circulation changes (Supplementary Fig. 7) and dominates over changes in diffusive nutrient transport (Supplementary Fig. 8). The southward Si redistribution causes nutrient trapping in the Southern Ocean with the Southern Ocean Si content increasing by 37% in 1,000 years (Supplementary Fig. 8). A further decomposition of the Si transport into different basins shows that the Southern Ocean Si increase arises from exchanges with the Indo-Pacific basin (red line in Fig. 3c), and thus, is not a circumpolar response as suggested in previous studies. The associated southward Si transport anomaly is due to transient overturning circulation responses, linked to the concomitant AMOC weakening, coupled to characteristic Si distributions with enriched Si and stronger vertical gradient in the Indo-Pacific (orange line in Fig. 3c). The Si transport is northward across 30° S in the Atlantic basin (blue line in Fig. 3c), largely due to the AMOC weakening (green line in Fig. 3c). Taken together, we show that transient overturning circulation changes can substantially impact large-scale nutrient transport and redistribution in the world oceans (Fig. 3).

We reach the same conclusion using previous model outputs2, i.e., the centennial nutrient redistribution between the Southern Ocean and the northern basins in a warming climate is predominantly due to transient overturning circulation adjustments between the Indo-Pacific and the Atlantic (Supplementary Fig. 9). While the previously proposed mechanism for Southern Ocean nutrient trapping may be useful to explain intermediate depth nutrient changes, we continue to find that overturning circulation changes and the associated isopycnal deepening are necessary to account for most of the nutrient changes in the subsurface Indo-Pacific and the upper Atlantic oceans (see “Nutrient changes due to adiabatic isopycnal movements” in “Methods”; Supplementary Fig. 10). Specifically, isopycnal surfaces deepen to balance the transient overturning adjustment23 and bring nutrient-depleted water to greater depths (blue arrows in Fig. 1), causing nutrient concentrations to decline in much of the global upper ocean.

Roles of wind and buoyancy

The poleward shift and strengthening of Southern Hemisphere westerly winds (Supplementary Fig. 11) has been thought to cause stronger upwelling of nutrient-rich deep waters15,32 and also contribute to nutrient increases in the Southern Ocean. However, in the CESM2 simulations, the Southern Ocean overturning circulation changes only briefly produce a small southward Si transport anomaly during the first 100 years, and then the sign reverses (pink line in Fig. 3d). Thus, the nutrient redistribution into the Southern Ocean is unlikely due to surface wind induced circumpolar upwelling changes2. To explore the underlying processes that affect the overturning circulation and nutrient redistribution, we assess the respective roles of surface winds and buoyancy changes, which have been extensively studied but heavily debated in the literature.

We use an ocean-biogeochemical model, based on MITgcm, to carry out a set of process-based numerical experiments, which allow us to separate the respective roles of wind and buoyancy in affecting the overturning circulation and nutrient redistribution (see “Ocean-biogeochemical model setup” in “Methods”). Under surface wind and buoyancy forcings from CESM2 4×CO2, our MITgcm runs largely reproduce the global ocean nutrient reorganization in CESM2 (Supplementary Fig. 12). Most of the global nutrient redistribution is attributable to changes in surface buoyancy forcing, rather than surface winds (Fig. 4 and Supplementary Fig. 13). The strengthening and poleward shift of Southern Hemisphere westerly winds increases Ekman pumping and surface nutrients only in regions close to Antarctica (Fig. 4a–c). In the broader Antarctic Circumpolar Current region, however, the southward shift of westerly winds weakens Ekman pumping and shifts the nutrient-rich deep waters to greater depths, leading to a decrease in surface nutrient concentrations (Supplementary Fig. 14; see “Nutrient changes due to adiabatic isopycnal movements” in “Methods”). Changes in surface buoyancy forcing, due to CO2 quadrupling, cause a reduction in nutrient concentrations in much of the global upper 2 km (Fig. 4d–f) that resembles nutrient changes in CESM2 4×CO2 (Fig. 2d, f), along with the southward nutrient redistribution driven by transient overturning changes (Supplementary Fig. 13c). Despite nutrient increases in the subsurface Southern Ocean, surface nutrient concentrations decline in the Southern Ocean. This likely occurs because surface warming and freshening increase surface stratification33 and reduce nutrient supply into the surface mixed layer from the deep ocean.Fig. 4 Si responses to surface forcing changes in the ocean-biogeochemical model, based on MITgcm.

a–c Si changes, ΔSi, when forced with surface wind from CESM2 4×CO2. d–f Si changes, ΔSi, when forced with surface buoyancy forcing from CESM2 4×CO2. ΔSi values are (a, d) vertically averaged in the upper 200 m and b–f zonally averaged in the b, e Atlantic and c, f Indo-Pacific basins. ΔSi is calculated as the difference from the MITgcm control run, which is forced with surface forcing taken from CESM2 piControl. For modeling details, see “Ocean-biogeochemical model setup” in “Methods”. Source data are provided as a Source Data file.

The MITgcm results show that changes in Southern Hemisphere westerly winds only marginally affect the global nutrient distribution. Therefore, we conclude that changes in surface buoyancy cause transient adjustments of the global overturning circulation that redistribute nutrients from the Indo-Pacific into the Southern Ocean.

Discussion

The deep Indo-Pacific contains a large amount of dissolved inorganic nutrients and carbon. Yet, previous studies on the ocean’s role in climate transitions have largely focused on physical processes that connect the overturning circulation between the Atlantic and the Southern Ocean. The impact of Indo-Pacific overturning changes remains critically under-explored. A change in the strength and structure of the AMOC is likely to be a key component of the climate’s response to anthropogenic forcing. Transient adjustments of the overturning circulation on centennial timescales associated with AMOC changes, as shown recently, involve interactions primarily between the Atlantic and the Indo-Pacific23 and could affect distributions of biogeochemical tracers, including nutrients and carbon. Here, we show that transient Indo-Pacific overturning circulation changes, associated with AMOC weakening, cause a southward nutrient redistribution from the Indo-Pacific into the Southern Ocean in a warming climate and lead to nutrient trapping in the Southern Ocean. This nutrient reorganization contributes to a long-term decline of oceanic biological productivity in the Indo-Pacific, as documented previously34–36.

The transient overturning circulation adjustments also provide mechanistic insights into past climate changes. For example, coral temperature records show significant subsurface warming in the eastern equatorial Pacific during Heinrich events and the Younger Dryas during the last deglaciation37. The subsurface warming on multi-centennial time scales is most likely due to isopycnal deepening in the broad low-to-mid latitudes that balances the Indo-Pacific overturning circulation anomalies associated with AMOC weakening25. The nutrient redistribution by transient overturning changes might have also supplied more Si to the Southern Ocean and contributed to increases in diatom production close to the Antarctic Polar Front during Heinrich events10,14. However, we note that diatom records during Heinrich stadials represent mean climate variations under different background climate conditions and on timescales (multi-centennial to millennial) longer than those discussed in this study. Further studies are required to understand processes controlling past ocean circulation, nutrient reorganizations, and the time scales over which these processes occur.

Nutrient cycling is intrinsically linked to the global ocean carbon cycle and atmospheric CO2 changes. Previous studies have heavily relied on Southern Ocean processes (e.g., changes to westerly winds) to explain changes in nutrient distribution and carbon cycle2,32,38. Our study suggests that processes in the North Atlantic, via transient overturning circulation changes between the Atlantic and the Indo-Pacific basins, play a critical role in regulating global nutrient distributions on centennial timescales. The global upper-ocean nutrient decline (Fig. 2) due to transient overturning changes suggests a more efficient biological carbon pump that may help increase oceanic CO2 uptake39. Dissolved inorganic carbon and Si have similar distributions in the ocean. Thus, transient overturning circulation adjustments, associated with AMOC changes, also inevitably redistribute carbon between different ocean basins and trap more carbon in the deep ocean as well as in the Southern Ocean (Supplementary Fig. 15), increasing deep ocean carbon storage.

Methods

CESM2 simulations

CESM2 is a fully coupled earth system model, consisting of ocean, atmosphere, land, sea ice, land ice, and wave components, which exchange states and fluxes through a coupler40. The ocean component is Parallel Ocean Program version 2 and has a nominal horizontal resolution of 1°. Ocean biogeochemistry in CESM2 is represented by the Marine Biogeochemistry Library (MARBL), which represents multiple nutrient colimitations (N, P, Si, and Fe). MARBL includes fully prognostic carbonate chemistry and simulates sinking particulate organic matter implicitly. There are three explicit phytoplankton functional groups (diatoms, diazotrophs, and pico/nano phytoplankton), one implicit group (calcifiers), and one zooplankton group. The piControl run approximately simulates the observed structure of nutrients in the upper 3 km (Fig. 2), but overestimates the abyssal nutrients level due to the too-weak abyssal overturning circulation41. The CESM2 CO2 quadrupling (4×CO2) experiment is initialized from piControl, but with the atmospheric CO2 instantaneously quadrupled. The 4×CO2 run is continued for another 1,000 years. Unless otherwise noted, we compare 4×CO2 averaged over years 150–200 with piControl.

During the first few decades, the reduced nutrient entrainment into the surface due to stronger stratification leads to a slight increase in nutrient concentrations below the surface mixed layer (Supplementary Fig. 5). The surface nutrient depletion favors picophytoplankton, which are better adapted to the low-nutrient environment, and their production increases due to the favorable warm growth conditions during the first few decades (green lines in Supplementary Fig. 2b, c). The picophytoplankton increase overwhelms the diatom decrease in primary production in the short term but does not contribute substantially to organic matter export to the subsurface ocean (Supplementary Fig. 3).

Regional changes in primary production are not exclusively controlled by nutrient concentrations (Supplementary Fig. 4). In low-to-mid latitudes, where diatom growth is somewhat limited by Si availability, diatom production tends to decrease in regions where Si concentration declines. However, primary production and organic matter export increase in the Southern Ocean (Supplementary Figs. 2 and 3), a typical high-nutrient low-chlorophyll region42, despite a decrease in surface Si concentration. This likely occurs due to a reduction in summer sea ice coverage close to Antarctica and surface warming in lower latitudes43, which relieves light limitation.

Nutrient transport

We calculate the nutrient transport across 30° S, T, as,1 T=∫zbot0∫xwxevNdxdz,

where N is nutrient concentration, v represents residual velocity, including the Eulerian-mean flow and eddies, xw and xe denote the western and eastern boundaries, respectively. Differences in nutrient transport between piControl and 4×CO2, ΔT, can be driven by changes in ocean circulation, v′, as well as changes in nutrient concentrations due to other processes, N′,2 ΔT=∫zbot0∫xwxev′N0+v0N′dxdz,

where N0 and v0 represents nutrient concentration and velocity, respectively, in piControl, and we neglected contributions due to v′N′. Our calculation suggests that circulation changes in 4×CO2 relative to piControl,ΔTv≡∫zbot0∫xwxev′N0dxdz,

dominates the nutrient transport anomaly, ΔT (Supplementary Fig. 7), i.e., ΔT≈ΔTv.

The nutrient transport changes due to circulation differences, ΔTv, can be further decomposed into contributions due to changes in ocean overturning circulation, ΔTψ, and horizontal circulation (e.g., gyres), ΔTG,3 ΔTv=∫zbot0∫xwxev′¯+v′^N0dxdz=∫zbot0Δψ∂N¯0∂zdz⏟ΔTψ+∫zbot0∫xwxev′^N^0dxdz⏟ΔTG,

where overline means zonal average between xw and xe, hat means deviations from zonal average, and Δψ represents overturning circulation differences between piControl and 4×CO2, defined as,4 Δψ(z)=∫zbotz∫xwxev′¯dxdz.

The nutrient exchange across 30° S, ΔT, occurs in both the Atlantic, ΔTA, and the Indo-Pacific sectors, TP,ΔT≡ΔTA+TP≈ΔTAv+ΔTPv.

Applying the decomposition in Eq. (3) for each basin, we have5 ΔT≈ΔTAψ+TPψ⏟ΔTiψ+ΔTAG+ΔTPG,

where the nutrient transport anomaly due to the compensating overturning transient between the Atlantic and Indo-Pacific basins,6 ΔTiψ=∫zbot0ΔψP∂N¯P∂z+ΔψA∂N¯A∂zdz≈∫zbot0ΔψP∂N¯P∂z−∂N¯A∂zdz<0

where N¯A and N¯P represent the zonal-mean nutrient concentration calculated in the Atlantic and the Indo-Pacific basin, respectively, in piControl. The Indo-Pacific overturning circulation anomaly at 30° S, ΔψP, is positive (Fig. 3b). The vertical nutrient gradient is much stronger in the Indo-Pacific basin, ∂N¯P∂z−∂N¯A∂z<0, due to the slower ventilation rate in the deep Pacific than in the deep Atlantic44. As such, the compensating overturning circulation transients drive a southward nutrient transport anomaly into the Southern Ocean.

Nutrient changes due to adiabatic isopycnal movements

The Indo-Pacific overturning responses in Fig. 3b are largely balanced by an adiabatic deepening of isopycnals. This process brings nutrient-depleted water to greater depth (Fig. 1) and contributes to the nutrient decrease in Fig. 2d. For conciseness, we consider a one-dimensional case. For piControl,7 N=N0(z),z=ζ0(b),b=B0(z),

and for 4×CO2,8 N=N1(z),z=ζ1(b),b=B1(z),

where N=N(z) represents nutrient distribution as a function of depth, z = ζ(b) denotes the depth of isopycnal of buoyancy b, and b = B(z) is the vertical buoyancy distribution. As such, the nutrient changes between piControl and 4×CO2 due to adiabatic isopycnal movements can be quantified as9 ΔNh(z)=N0ζ0(B1(z))−N0(z),

Changes in biological processes, horizontal advection, and diffusion may also contribute to the nutrient changes and are represented as the residual between the total nutrient changes, ΔN, and ΔNh.

Ocean-biogeochemical model setup

We carry out a suite of ocean-biogeochemical simulations to compare the role of wind and surface buoyancy forcing in global nutrient redistribution. The ocean component is the Massachusetts Institute of Technology General Circulation Model (MITgcm) with a nonlinear equation of state45. The model resolution is 1° in longitude and 0.5° in latitude. The model domain is similar to CESM2 but without the Arctic to avoid the singularity problem in the polar region. In the vertical, there are 30 layers ranging from 10 m at the top to 369 m at the bottom. Vertical mixing is handled by the K-Profile parameterization46. Unresolved eddies are represented using the skew-flux form of the Gent–McWilliams parameterization with a background diffusivity of 1,000 m2 s−1. Isopycnal diffusion due to eddies is parameterized using the Redi scheme with a diffusivity of 1,000 m2 s−1. The model is forced with surface forcing taken from CESM2. Temperature at the surface is restored to monthly sea surface temperature with a relaxation timescale of 14 days47 and a net surface heat flux. Surface salinity is forced by monthly surface freshwater flux, which includes precipitation, evaporation, ice melt, and river runoff. A weak relaxation of the surface salinity to CESM2 piControl climatology is also included to avoid possible bi-stability of the overturning circulation due to the mixed boundary conditions48.

The biogeochemistry component is a modified version of Biogeochemistry with Light, Iron, Nutrients, and Gas (BLING)49,50 to include the silicon cycle. The original BLING model includes an implicit physiological model of growth co-limitation by light, temperature, phosphate, and iron and an implicit ecological model of partitioning into remineralized and dissolved pools. It solves for dissolved inorganic phosphate, dissolved inorganic carbon, dissolved organic phosphate, alkalinity, oxygen, and dissolved iron. The modifications to BLING include an added phytoplankton group, diatoms, and silicon cycle, largely following a more sophisticated biogeochemical model, TOPAZv251. We force the biogeochemistry model with climatological monthly iron fluxes, surface wind speed, and sea ice coverage from CESM2 piControl, which are kept constant throughout the suite of different simulations.

A total of four experiments were carried out using the ocean-biogeochemical model: “Control”, “Wind”, “Buoyancy”, and “4×CO2”. In “Control”, both surface wind and buoyancy forcing (surface heat flux, surface temperature, surface freshwater flux, and surface salinity) are taken from CESM2 piControl. “Control” is run for over 4,000 years to achieve an approximately steady state. The other three experiments were initialized from the end of “Control” and continued for 300 years. In “Wind”, the surface wind is taken from CESM2 4×CO2 (monthly climatology over years 150–200) but surface buoyancy forcing is kept the same as in “Control”. In “Buoyancy”, the surface wind is kept the same as in “Control” but surface buoyancy forcing is taken from CESM2 4×CO2 (monthly climatology over years 150–200). In “4×CO2”, both surface wind and buoyancy forcing are taken from CESM2 4×CO2 (monthly climatology over years 150–200). “Control” largely reproduces the structure of the nutrients in CESM2 (Supplementary Fig. 12). With surface forcing from CESM2 4×CO2, “4×CO2” also approximately reproduces the global nutrients reorganization (Supplementary Fig. 12). Comparison of “Wind” and “Buoyancy” with “Control”, quantified as their differences averaged over years 250–300 after perturbations, shows the impacts of changes in surface wind and surface buoyancy forcing, respectively, on the overturning responses and nutrients reorganization.

Supplementary information

Supplementary Information

Peer Review File

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52200-0.

Acknowledgements

Without implying their endorsement, we are grateful for helpful discussions with Daniel Sigman, Jess Adkins, Keith Moore, Frankie Pavia, and Ariane Verdy. J.Y. acknowledges support from NSF China (42330403, 42076056). A.F.T. was supported by NSF OCE-2023259, the Resnick Sustainability Institute, and the David and Lucille Packard Foundation.

Author contributions

S.S. and A.F.T. conceived the study. S.S. carried out the model simulations and analyzed the model output. S.S. and A.F.T. drafted the manuscript. J.Y. and L.W. contributed to discussions and improving the manuscript.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

The CESM2 model output was downloaded from the Earth System Grid Federation node (https://esgf-node.llnl.gov/search/cmip6/). The climatological Si data used in creating Fig. 1 were downloaded from GLODAPv2 (https://glodap.info/). Source data for Figs. 2–4 are provided in this paper. Model outputs from the MITgcm experiments are made publicly available from the online open access repository, Figshare, with 10.6084/m9.figshare.26787751. Source data are provided in this paper.

Code availability

The code for reproducing the ocean-biogeochemical experiments in this study is made publicly available from the online open-access repository, Figshare, with 10.6084/m9.figshare.26129629.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Bopp L Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models Biogeosciences 2013 10 6225 6245 10.5194/bg-10-6225-2013
Bopp, L. et al. Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models. Biogeosciences 10, 6225–6245 (2013).10.5194/bg-10-6225-2013
2. Moore JK Sustained climate warming drives declining marine biological productivity Science 2018 359 1139 1143 10.1126/science.aao6379 29590043
Moore, J. K. et al. Sustained climate warming drives declining marine biological productivity. Science 359, 1139–1143 (2018).29590043 10.1126/science.aao6379
3. Kwiatkowski L Twenty-first century ocean warming, acidification, deoxygenation, and upper-ocean nutrient and primary production decline from CMIP6 model projections Biogeosciences 2020 17 3439 3470 10.5194/bg-17-3439-2020
Kwiatkowski, L. et al. Twenty-first century ocean warming, acidification, deoxygenation, and upper-ocean nutrient and primary production decline from CMIP6 model projections. Biogeosciences 17, 3439–3470 (2020).10.5194/bg-17-3439-2020
4. Held IM Probing the fast and slow components of global warming by returning abruptly to preindustrial forcing J. Clim. 2010 23 2418 2427 10.1175/2009JCLI3466.1
Held, I. M. et al. Probing the fast and slow components of global warming by returning abruptly to preindustrial forcing. J. Clim. 23, 2418–2427 (2010).10.1175/2009JCLI3466.1
5. Long S-M Xie S-P Zheng X-T Liu Q Fast and slow responses to global warming: sea surface temperature and precipitation patterns J. Clim. 2014 27 285 299 10.1175/JCLI-D-13-00297.1
Long, S.-M., Xie, S.-P., Zheng, X.-T. & Liu, Q. Fast and slow responses to global warming: sea surface temperature and precipitation patterns. J. Clim. 27, 285–299 (2014).10.1175/JCLI-D-13-00297.1
6. Liu, Y., Moore, J. K., Primeau, F. & Wang, W. L. Reduced CO2 uptake and growing nutrient sequestration from slowing overturning circulation. Nat. Clim. Change. 10.1038/s41558-022-01555-7 (2022).
7. Yin, J. H. A consistent poleward shift of the storm tracks in simulations of 21st-century climate. Geophys. Res. Lett. 32, L18701 (2005).
8. Goyal R Sen Gupta A Jucker M England MH Historical and projected changes in the southern hemisphere surface westerlies Geophys. Res. Lett. 2021 48 e2020GL090849 10.1029/2020GL090849
Goyal, R., Sen Gupta, A., Jucker, M. & England, M. H. Historical and projected changes in the southern hemisphere surface westerlies. Geophys. Res. Lett. 48, e2020GL090849 (2021).10.1029/2020GL090849
9. Sarmiento JL Gruber N Brzezinski M Dunne J High-latitude controls of thermocline nutrients and low-latitude biological productivity Nature 2004 427 56 60 10.1038/nature02127 14702082
Sarmiento, J. L., Gruber, N., Brzezinski, M. & Dunne, J. High-latitude controls of thermocline nutrients and low-latitude biological productivity. Nature 427, 56–60 (2004).14702082 10.1038/nature02127
10. Anderson RF Wind-driven upwelling in the southern ocean and the deglacial rise in atmospheric CO2 Science 2009 323 1443 1448 10.1126/science.1167441 19286547
Anderson, R. F. et al. Wind-driven upwelling in the southern ocean and the deglacial rise in atmospheric CO2. Science 323, 1443–1448 (2009).19286547 10.1126/science.1167441
11. Dumont M The nature of deep overturning and reconfigurations of the silicon cycle across the last deglaciation Nat. Commun. 2020 11 1 11 10.1038/s41467-020-15101-6 31911652
Dumont, M. et al. The nature of deep overturning and reconfigurations of the silicon cycle across the last deglaciation. Nat. Commun. 11, 1–11 (2020).31911652 10.1038/s41467-020-15101-6
12. Dansgaard W North Atlantic climatic oscillations revealed by deep Greenland ice cores Clim. Process. Clim. Sensit. 1984 29 288 298 10.1029/GM029p0288
Dansgaard, W. et al. North Atlantic climatic oscillations revealed by deep Greenland ice cores. Clim. Process. Clim. Sensit. 29, 288–298 (1984).10.1029/GM029p0288
13. Menviel L Skinner LC Tarasov L Tzedakis PC An ice–climate oscillatory framework for Dansgaard–Oeschger cycles Nat. Rev. Earth Environ. 2020 1 677 693 10.1038/s43017-020-00106-y
Menviel, L., Skinner, L. C., Tarasov, L. & Tzedakis, P. C. An ice–climate oscillatory framework for Dansgaard–Oeschger cycles. Nat. Rev. Earth Environ. 1, 677–693 (2020).10.1038/s43017-020-00106-y
14. Sigman DM Hain MP Haug GH The polar ocean and glacial cycles in atmospheric CO2 concentration Nature 2010 466 47 55 10.1038/nature09149 20596012
Sigman, D. M., Hain, M. P. & Haug, G. H. The polar ocean and glacial cycles in atmospheric CO2 concentration. Nature 466, 47–55 (2010).20596012 10.1038/nature09149
15. Toggweiler J Samuels B Effect of Drake passage on the global thermohaline circulation Deep Sea Res. I 1995 42 477 500 10.1016/0967-0637(95)00012-U
Toggweiler, J. & Samuels, B. Effect of Drake passage on the global thermohaline circulation. Deep Sea Res. I 42, 477–500 (1995).10.1016/0967-0637(95)00012-U
16. Toggweiler, J., Russell, J. L. & Carson, S. R. Midlatitude westerlies, atmospheric CO2, and climate change during the ice ages. Paleoceanogr. 10.1029/2005PA001154 (2006).
17. Abernathey R Marshall J Ferreira D The dependence of Southern Ocean meridional overturning on wind stress J. Phys. Oceanogr. 2011 41 2261 2278 10.1175/JPO-D-11-023.1
Abernathey, R., Marshall, J. & Ferreira, D. The dependence of Southern Ocean meridional overturning on wind stress. J. Phys. Oceanogr. 41, 2261–2278 (2011).10.1175/JPO-D-11-023.1
18. Dufour CO Standing and transient eddies in the response of the Southern Ocean meridional overturning to the southern annular mode J. Clim. 2012 25 6958 6974 10.1175/JCLI-D-11-00309.1
Dufour, C. O. et al. Standing and transient eddies in the response of the Southern Ocean meridional overturning to the southern annular mode. J. Clim. 25, 6958–6974 (2012).10.1175/JCLI-D-11-00309.1
19. Meredith MP Naveira Garabato AC Hogg AM Farneti R Sensitivity of the overturning circulation in the Southern Ocean to decadal changes in wind forcing J. Clim. 2012 25 99 110 10.1175/2011JCLI4204.1
Meredith, M. P., Naveira Garabato, A. C., Hogg, A. M. & Farneti, R. Sensitivity of the overturning circulation in the Southern Ocean to decadal changes in wind forcing. J. Clim. 25, 99–110 (2012).10.1175/2011JCLI4204.1
20. Bishop SP Southern ocean overturning compensation in an Eddy-resolving climate simulation J. Phys. Oceanogr. 2016 46 1575 1592 10.1175/JPO-D-15-0177.1
Bishop, S. P. et al. Southern ocean overturning compensation in an Eddy-resolving climate simulation. J. Phys. Oceanogr. 46, 1575–1592 (2016).10.1175/JPO-D-15-0177.1
21. Cheng W Chiang JC Zhang D Atlantic meridional overturning circulation (AMOC) in CMIP5 models: RCP and historical simulations J. Clim. 2013 26 7187 7197 10.1175/JCLI-D-12-00496.1
Cheng, W., Chiang, J. C. & Zhang, D. Atlantic meridional overturning circulation (AMOC) in CMIP5 models: RCP and historical simulations. J. Clim. 26, 7187–7197 (2013).10.1175/JCLI-D-12-00496.1
22. Weijer W Cheng W Garuba O Hu A Nadiga B CMIP6 models predict significant 21st century decline of the Atlantic Meridional Overturning Circulation Geophys. Res. Lett. 2020 47 e2019GL086075 10.1029/2019GL086075
Weijer, W., Cheng, W., Garuba, O., Hu, A. & Nadiga, B. CMIP6 models predict significant 21st century decline of the Atlantic Meridional Overturning Circulation. Geophys. Res. Lett. 47, e2019GL086075 (2020).10.1029/2019GL086075
23. Sun S Thompson AF Eisenman I Transient overturning compensation between Atlantic and Indo-Pacific basins J. Phys. Oceanogr. 2020 50 2151 2172 10.1175/JPO-D-20-0060.1
Sun, S., Thompson, A. F. & Eisenman, I. Transient overturning compensation between Atlantic and Indo-Pacific basins. J. Phys. Oceanogr. 50, 2151–2172 (2020).10.1175/JPO-D-20-0060.1
24. Sun S Thompson AF Centennial changes in the Indonesian throughflow connected to the Atlantic Meridional overturning circulation: the ocean’s transient conveyor belt Geophys. Res. Lett. 2020 47 e2020GL090615 10.1029/2020GL090615
Sun, S. & Thompson, A. F. Centennial changes in the Indonesian throughflow connected to the Atlantic Meridional overturning circulation: the ocean’s transient conveyor belt. Geophys. Res. Lett. 47, e2020GL090615 (2020).10.1029/2020GL090615
25. Sun S Thompson AF Xie S-P Long S-M Indo-Pacific warming induced by a weakening of the Atlantic meridional overturning circulation J. Clim. 2022 35 815 832 10.1175/JCLI-D-21-0346.1
Sun, S., Thompson, A. F., Xie, S.-P. & Long, S.-M. Indo-Pacific warming induced by a weakening of the Atlantic meridional overturning circulation. J. Clim. 35, 815–832 (2022).10.1175/JCLI-D-21-0346.1
26. Saenko OA Schmittner A Weaver AJ The Atlantic–Pacific seesaw J. Clim. 2004 17 2033 2038 10.1175/1520-0442(2004)017<2033:TAS>2.0.CO;2
Saenko, O. A., Schmittner, A. & Weaver, A. J. The Atlantic–Pacific seesaw. J. Clim. 17, 2033–2038 (2004).10.1175/1520-0442(2004)017<2033:TAS>2.0.CO;2
27. Okazaki Y Deepwater formation in the North Pacific during the last glacial termination Science 2010 329 200 204 10.1126/science.1190612 20616275
Okazaki, Y. et al. Deepwater formation in the North Pacific during the last glacial termination. Science 329, 200–204 (2010).20616275 10.1126/science.1190612
28. Chikamoto MO Variability in North Pacific intermediate and deep water ventilation during Heinrich events in two coupled climate models Deep Sea Res. II 2012 61 114 126 10.1016/j.dsr2.2011.12.002
Chikamoto, M. O. et al. Variability in North Pacific intermediate and deep water ventilation during Heinrich events in two coupled climate models. Deep Sea Res. II 61, 114–126 (2012).10.1016/j.dsr2.2011.12.002
29. Menviel L England MH Meissner K Mouchet A Yu J Atlantic–Pacific seesaw and its role in outgassing CO2 during Heinrich events Paleoceanogr. 2014 29 58 70 10.1002/2013PA002542
Menviel, L., England, M. H., Meissner, K., Mouchet, A. & Yu, J. Atlantic–Pacific seesaw and its role in outgassing CO2 during Heinrich events. Paleoceanogr. 29, 58–70 (2014).10.1002/2013PA002542
30. Eyring V Overview of the coupled model intercomparison project phase 6 (CMIP6) experimental design and organization Geosci. Model Dev. 2016 9 1937 1958 10.5194/gmd-9-1937-2016
Eyring, V. et al. Overview of the coupled model intercomparison project phase 6 (CMIP6) experimental design and organization. Geosci. Model Dev. 9, 1937–1958 (2016).10.5194/gmd-9-1937-2016
31. Sarmiento, J. L. & Gruber, N. Ocean Biogeochemical Dynamics (Princeton University Press, 2006).
32. Menviel, L. et al. Southern Hemisphere westerlies as a driver of the early deglacial atmospheric CO2 rise. Nat. Commun. 9, 2503 (2018).
33. Swart NC Gille ST Fyfe JC Gillett NP Recent southern ocean warming and freshening driven by greenhouse gas emissions and ozone depletion Nat. Geosci. 2018 11 836 10.1038/s41561-018-0226-1
Swart, N. C., Gille, S. T., Fyfe, J. C. & Gillett, N. P. Recent southern ocean warming and freshening driven by greenhouse gas emissions and ozone depletion. Nat. Geosci. 11, 836 (2018).10.1038/s41561-018-0226-1
34. Schmittner A Decline of the marine ecosystem caused by a reduction in the Atlantic overturning circulation Nature 2005 434 628 633 10.1038/nature03476 15800620
Schmittner, A. Decline of the marine ecosystem caused by a reduction in the Atlantic overturning circulation. Nature 434, 628–633 (2005).15800620 10.1038/nature03476
35. Schmittner A Galbraith ED Glacial greenhouse-gas fluctuations controlled by ocean circulation changes Nature 2008 456 373 376 10.1038/nature07531 19020618
Schmittner, A. & Galbraith, E. D. Glacial greenhouse-gas fluctuations controlled by ocean circulation changes. Nature 456, 373–376 (2008).19020618 10.1038/nature07531
36. Nielsen SB Jochum M Pedro JB Eden C Nuterman R Two-timescale carbon cycle response to an AMOC collapse Paleoceanogr. Paleoclimatol. 2019 34 511 523 10.1029/2018PA003481
Nielsen, S. B., Jochum, M., Pedro, J. B., Eden, C. & Nuterman, R. Two-timescale carbon cycle response to an AMOC collapse. Paleoceanogr. Paleoclimatol. 34, 511–523 (2019).10.1029/2018PA003481
37. Stewart JA Arctic and Antarctic forcing of ocean interior warming during the last deglaciation Sci. Rep. 2023 13 22410 10.1038/s41598-023-49435-0 38104174
Stewart, J. A. et al. Arctic and Antarctic forcing of ocean interior warming during the last deglaciation. Sci. Rep. 13, 22410 (2023).38104174 10.1038/s41598-023-49435-0
38. Sigman, D. M. et al. The Southern Ocean during the Ice Ages: a review of the Antarctic surface isolation hypothesis, with comparison to the North Pacific. Quat. Sci. Rev. 254, 106732 (2020).
39. Ito T Follows MJ Preformed phosphate, soft tissue pump and atmospheric CO2 J. Mar. Res. 2005 63 813 839 10.1357/0022240054663231
Ito, T. & Follows, M. J. Preformed phosphate, soft tissue pump and atmospheric CO2. J. Mar. Res. 63, 813–839 (2005).10.1357/0022240054663231
40. Danabasoglu G The community earth system model version 2 (CESM2) J. Adv. Model. Earth Syst. 2020 12 e2019MS001916 10.1029/2019MS001916
Danabasoglu, G. et al. The community earth system model version 2 (CESM2). J. Adv. Model. Earth Syst. 12, e2019MS001916 (2020).10.1029/2019MS001916
41. Heuzé C Heywood KJ Stevens DP Ridley JK Southern Ocean bottom water characteristics in CMIP5 models Geophys. Res. Lett. 2013 40 1409 1414 10.1002/grl.50287
Heuzé, C., Heywood, K. J., Stevens, D. P. & Ridley, J. K. Southern Ocean bottom water characteristics in CMIP5 models. Geophys. Res. Lett. 40, 1409–1414 (2013).10.1002/grl.50287
42. Martin JH Fitzwater SE Gordon RM Iron deficiency limits phytoplankton growth in Antarctic waters Glob. Biogeochem. Cycles 1990 4 5 12 10.1029/GB004i001p00005
Martin, J. H., Fitzwater, S. E. & Gordon, R. M. Iron deficiency limits phytoplankton growth in Antarctic waters. Glob. Biogeochem. Cycles 4, 5–12 (1990).10.1029/GB004i001p00005
43. Hauck J On the Southern Ocean CO2 uptake and the role of the biological carbon pump in the 21st century Global Biogeochem. Cycle 2015 29 1451 1470 10.1002/2015GB005140
Hauck, J. et al. On the Southern Ocean CO2 uptake and the role of the biological carbon pump in the 21st century. Global Biogeochem. Cycle 29, 1451–1470 (2015).10.1002/2015GB005140
44. Talley LD Closure of the global overturning circulation through the Indian, Pacific, and Southern Oceans: schematics and transports Oceanography 2013 26 80 97 10.5670/oceanog.2013.07
Talley, L. D. Closure of the global overturning circulation through the Indian, Pacific, and Southern Oceans: schematics and transports. Oceanography 26, 80–97 (2013).10.5670/oceanog.2013.07
45. Jackett DR Mcdougall TJ Minimal adjustment of hydrographic profiles to achieve static stability J. Atmos. Ocean. Technol. 1995 12 381 389 10.1175/1520-0426(1995)012<0381:MAOHPT>2.0.CO;2
Jackett, D. R. & Mcdougall, T. J. Minimal adjustment of hydrographic profiles to achieve static stability. J. Atmos. Ocean. Technol. 12, 381–389 (1995).10.1175/1520-0426(1995)012<0381:MAOHPT>2.0.CO;2
46. Large WG McWilliams JC Doney SC Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization Rev. Geophys. 1994 32 363 403 10.1029/94RG01872
Large, W. G., McWilliams, J. C. & Doney, S. C. Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization. Rev. Geophys. 32, 363–403 (1994).10.1029/94RG01872
47. Haney RL Surface thermal boundary condition for ocean circulation models J. Phys. Oceanogr. 1971 1 241 248 10.1175/1520-0485(1971)001<0241:STBCFO>2.0.CO;2
Haney, R. L. Surface thermal boundary condition for ocean circulation models. J. Phys. Oceanogr. 1, 241–248 (1971).10.1175/1520-0485(1971)001<0241:STBCFO>2.0.CO;2
48. Stommel H Thermohaline convection with two stable regimes of flow Tellus 1961 13 224 230 10.3402/tellusa.v13i2.9491
Stommel, H. Thermohaline convection with two stable regimes of flow. Tellus 13, 224–230 (1961).10.3402/tellusa.v13i2.9491
49. Dunne JP Simple global ocean biogeochemistry with light, iron, nutrients and gas version 2 (BLINGv2): model description and simulation characteristics in GFDL’s CM4. 0 J. Adv. Model. Earth Syst. 2020 12 e2019MS002008 10.1029/2019MS002008
Dunne, J. P. et al. Simple global ocean biogeochemistry with light, iron, nutrients and gas version 2 (BLINGv2): model description and simulation characteristics in GFDL’s CM4. 0. J. Adv. Model. Earth Syst. 12, e2019MS002008 (2020).10.1029/2019MS002008
50. Verdy A Mazloff MR A data assimilating model for estimating Southern Ocean biogeochemistry J. Geophys. Res. Oceans 2017 122 6968 6988 10.1002/2016JC012650
Verdy, A. & Mazloff, M. R. A data assimilating model for estimating Southern Ocean biogeochemistry. J. Geophys. Res. Oceans 122, 6968–6988 (2017).10.1002/2016JC012650
51. Dunne JP GFDL’s ESM2 global coupled climate-carbon earth system models. Part I: Physical formulation and baseline simulation characteristics J. Clim. 2012 25 6646 6665 10.1175/JCLI-D-11-00560.1
Dunne, J. P. et al. GFDL’s ESM2 global coupled climate-carbon earth system models. Part I: Physical formulation and baseline simulation characteristics. J. Clim. 25, 6646–6665 (2012).10.1175/JCLI-D-11-00560.1
52. Olsen A The Global Ocean Data Analysis Project version 2 (GLODAPv2)—an internally consistent data product for the world ocean Earth Syst. Sci. Data 2016 8 297 323 10.5194/essd-8-297-2016
Olsen, A. et al. The Global Ocean Data Analysis Project version 2 (GLODAPv2)—an internally consistent data product for the world ocean. Earth Syst. Sci. Data 8, 297–323 (2016).10.5194/essd-8-297-2016
