Earth’s climate has undergone dramatic shifts between glacial and interglacial time periods, with high-latitude temperature changes on the order of 5–10 °C. These climatic shifts have been associated ...with major rearrangements in the deep ocean circulation and stratification, which have likely played an important role in the observed atmospheric carbon dioxide swings by affecting the partitioning of carbon between the atmosphere and the ocean. The mechanisms by which the deep ocean circulation changed, however, are still unclear and represent a major challenge to our understanding of glacial climates. This study shows that various inferred changes in the deep ocean circulation and stratification between glacial and interglacial climates can be interpreted as a direct consequence of atmospheric temperature differences. Colder atmospheric temperatures lead to increased sea ice cover and formation rate around Antarctica. The associated enhanced brine rejection leads to a strongly increased deep ocean stratification, consistent with high abyssal salinities inferred for the last glacial maximum. The increased stratification goes together with a weakening and shoaling of the interhemispheric overturning circulation, again consistent with proxy evidence for the last glacial. The shallower interhemispheric overturning circulation makes room for slowly moving water of Antarctic origin, which explains the observed middepth radiocarbon age maximum and may play an important role in ocean carbon storage.
Abstract
Antarctic sea‐ice formation plays a key role in shaping the abyssal overturning circulation and stratification in all ocean basins, by driving surface buoyancy loss through the associated ...brine rejection. Changes in Antarctic sea ice have therefore been suggested as drivers of major glacial‐interglacial ocean circulation rearrangements. Here, the relationship between Antarctic sea ice, buoyancy loss, deep‐ocean stratification, and overturning circulation is investigated in Last Glacial Maximum and preindustrial simulations from the Paleoclimate Modelling Intercomparison Project (PMIP). The simulations show substantial intermodel differences in their representation of the glacial deep‐ocean state and circulation, which is often at odds with the geological evidence. We argue that these apparent inconsistencies can largely be attributed to differing (and likely insufficient) Antarctic sea‐ice formation. Discrepancies can be further amplified by short integration times. Deep‐ocean equilibration and sea‐ice representation should, therefore, be carefully evaluated in the forthcoming PMIP4 simulations.
Key Points
Antarctic sea ice, buoyancy loss, abyssal stratification, and ocean overturning circulation are related across glacial and modern simulations
Insufficient sea ice formation and short integration times lead to discrepancies between different glacial climate simulations and paleodata
In the modern climate, the ocean below 2 km is mainly filled by waters sinking into the abyss around Antarctica and in the North Atlantic. Paleoproxies indicate that waters of North Atlantic origin ...were instead absent below 2 km at the Last Glacial Maximum, resulting in an expansion of the volume occupied by Antarctic origin waters. In this study we show that this rearrangement of deep water masses is dynamically linked to the expansion of summer sea ice around Antarctica. A simple theory further suggests that these deep waters only came to the surface under sea ice, which insulated them from atmospheric forcing, and were weakly mixed with overlying waters, thus being able to store carbon for long times. This unappreciated link between the expansion of sea ice and the appearance of a voluminous and insulated water mass may help quantify the ocean’s role in regulating atmospheric carbon dioxide on glacial–interglacial timescales. Previous studies pointed to many independent changes in ocean physics to account for the observed swings in atmospheric carbon dioxide. Here it is shown that many of these changes are dynamically linked and therefore must co-occur.
Abstract
Theoretical arguments are developed to derive general properties of the ocean circulation in a “snowball” world, and the predictions are confirmed in a series of idealized numerical ...simulations. As suggested previously, a turbulent flow is driven by geothermal heating at the seafloor, which is balanced by a similar heat loss through the ice sheet above. It is argued that the expected horizontal inhomogeneities in these heat fluxes are sufficient to generate baroclinic instability, which drives geostrophic turbulence. Turbulent eddies then transport heat upward and poleward along isolines of constant density, thereby maintaining a statically stable stratification, contrary to previous findings from numerical models that do not adequately resolve the geostrophic turbulence. The kinetic energy of the turbulent flow is expected to be controlled by a balance between the potential energy input by the diabatic forcing and frictional dissipation in the bottom boundary layer. The resulting characteristic flow speed is estimated to be on the order of 1 cm s
−1
, which is in agreement with previous numerical simulations. Eddy diffusivities are estimated to be on the order of 100 m
2
s
−1
, which is smaller than in the present-day ocean but probably within one order of magnitude. Because of the weak forcing, the resulting gradients of temperature and salinity would be much smaller than in the present-day ocean, with global-scale potential temperature variations on the order of 0.1 K, again in agreement with previous numerical simulations. The presented theoretical arguments may also be relevant to other planetary bodies with an ice-covered ocean.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Much of the existing theory for the ocean’s overturning circulation considers steady-state equilibrium solutions. However, Earth’s climate is not in a steady state, and a better understanding of the ...ocean’s non-equilibrium response to changes in the surface climate is urgently needed. Here, the time-dependent response of the deep-ocean overturning circulation to atmospheric warming is examined using a hierarchy of idealized ocean models. The transient response to surface warming is characterized by a shoaling and weakening of the Atlantic meridional overturning circulation (AMOC)—consistent with results from coupled climate simulations. The initial shoaling and weakening of the AMOC occurs on decadal time scales and is attributed to a rapid warming of northern-sourced deep water. The equilibrium response to warming, in contrast, is associated with a deepening and strengthening of the AMOC. The eventual deepening of the AMOC is argued to be associated with abyssal density changes and driven by modified surface fluxes in the Southern Ocean, following a reduction of the Antarctic sea ice cover. Full equilibration of the AMOC requires a diffusive adjustment of the abyss and takes many millennia. The equilibration time scale is much longer than most coupled climate model simulations, highlighting the importance of considering integration time and initial conditions when interpreting the deep-ocean circulation in climate models. The results also show that past climates are unlikely to be an adequate analog for changes in the overturning circulation during the coming decades or centuries.
Celotno besedilo
Dostopno za:
BFBNIB, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abstract
It remains uncertain how the Southern Ocean circulation responds to changes in surface wind stress, and whether coarse-resolution simulations, where mesoscale eddy fluxes are parameterized, ...can adequately capture the response. We address this problem using two idealized model setups mimicking the Southern Ocean: a flat-bottom channel and a channel with moderately complex topography. Under each topographic configuration and varying wind stress, we compare several coarse-resolution simulations, configured with different eddy parameterizations, against an eddy-resolving simulation. We find that 1) without topography, sensitivity of the Antarctic Circumpolar Current (ACC) to wind stress is overestimated by coarse-resolution simulations, due to an underestimate of the sensitivity of the eddy diffusivity; 2) in the presence of topography, stationary eddies dominate over transient eddies in counteracting the direct response of the ACC and overturning circulation to wind stress changes; and 3) coarse-resolution simulations with parameterized eddies capture this counteracting effect reasonably well, largely due to their ability to resolve stationary eddies. Our results highlight the importance of topography in modulating the response of the Southern Ocean circulation to changes in surface wind stress. The interaction between mesoscale eddies and stationary meanders induced by topography requires more attention in future development and testing of eddy parameterizations.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
A subgrid‐scale eddy parameterization is developed, which makes use of an explicit eddy kinetic energy budget and can be applied at both “non‐eddying” and “eddy‐permitting” resolutions. The ...subgrid‐scale eddies exchange energy with the resolved flow in both directions via a parameterization of baroclinic instability (based on the established formulation of Gent and McWilliams) and biharmonic and negative Laplacian viscosity terms. This formulation represents the turbulent cascades of energy and enstrophy consistent with our current understanding of the turbulent eddy energy cycle. At the same time, the approach is simple and general enough to be readily implemented in ocean climate models, without adding significant computational cost. The closure has been implemented in the Modular Ocean Model Version 6 and tested in the “Neverworld” configuration, which employs an idealized analytically defined topography designed as a testbed for mesoscale eddy parameterizations. The parameterization performs well over a range of resolutions, seamlessly connecting non‐eddying and eddy‐resolving regimes.
Key Points
A scale‐aware energy budget‐based eddy parameterization is introduced
Bidirectional energy transfer between resolved flow and subgrid scales can be represented
The parameterization allows for a smooth transition between non‐eddying and eddying resolution regimes
Abstract
Although the reconfiguration of the abyssal overturning circulation has been argued to be a salient feature of Earth’s past climate changes, our understanding of the physical mechanisms ...controlling its strength remains limited. In particular, existing scaling theories disagree on the relative importance of the dynamics in the Southern Ocean versus the dynamics in the basins to the north. In this study, we systematically investigate these theories and compare them with a set of numerical simulations generated from an ocean general circulation model with idealized geometry, designed to capture only the basic ingredients considered by the theories. It is shown that the disagreement between existing theories can be partially explained by the fact that the overturning strengths measured in the channel and in the basin scale distinctly with the external parameters, including surface buoyancy loss, diapycnal diffusivity, wind stress, and eddy diffusivity. The overturning in the reentrant channel, which represents the Southern Ocean, is found to be sensitive to all these parameters, in addition to a strong dependence on bottom topography. By contrast, the basin overturning varies with the integrated surface buoyancy loss rate and diapycnal diffusivity but is mostly unaffected by winds and channel topography. The simulated parameter dependence of the basin overturning can be described by a scaling theory that is based only on basin dynamics.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL ...climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in pure z* vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution.
Key Points
Documentation is provided for a new generation of NOAA‐GFDL CMIP6/OMIP ocean ice climate models
Dynamical core and physical parameterizations are described and key features of interannual CORE simulations are assessed
Using hybrid vertical coordinates reduces spurious ocean heat drift