Significance The ocean’s role in regulating atmospheric carbon dioxide on glacial–interglacial timescales remains an unresolved issue in paleoclimatology. Many apparently independent changes in ocean physics, chemistry, and biology need to be invoked to explain the full signal. Recent understanding of the deep ocean circulation and stratification is used to demonstrate that the major changes invoked in ocean physics are dynamically linked. In particular, the expansion of permanent sea ice in the Southern Hemisphere results in a volume increase of Antarctic-origin abyssal waters and a reduction in mixing between abyssal waters of Arctic and Antarctic origin.

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.

We review recent observations and models concerning the dynamics of Cryogenian global glaciation and their biological consequences.

No abstract

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.

Statistical analysis of observations (including atmospheric reanalysis and forced ocean model simulations) is used to address two questions: First, does an analogous mechanism to that of El Niñ o-Southern Oscillation (ENSO) exist in the equatorial Atlantic or Indian Ocean? Second, does the intrinsic variability in these basins matter for ENSO predictability? These questions are addressed by assessing the existence and strength of the Bjerknes and delayed negative feedbacks in each tropical basin, and by fitting conceptual recharge oscillator models, both with and without interactions among the basins.In the equatorial Atlantic the Bjerknes and delayed negative feedbacks exist, although weaker than in the Pacific. Equatorial Atlantic variability is well described by the recharge oscillator model, with an oscillatory mixed ocean dynamics-sea surface temperature (SST) mode present in boreal spring and summer. The dynamics of the tropical Indian Ocean, however, appear to be quite different: no recharge-discharge mechanism is found. Although a positive Bjerknes-like feedback from July to September is found, the role of heat content seems secondary.Results also show that Indian Ocean interaction with ENSO tends to damp the ENSO oscillation and is responsible for a frequency shift to shorter periods. However, the retrospective forecast skill of the conceptual model is hardly improved by explicitly including Indian Ocean SST. The interaction between ENSO and the equatorial Atlantic variability is weaker. However, a feedback from the Atlantic on ENSO appears to exist, which slightly improves the retrospective forecast skill of the conceptual model.

Recent studies suggest that the enhanced upper ocean mixing caused by tropical cyclones significantly contributes to the ocean heat transport. However, existing studies that try to quantify this contribution make the assumption that all heat pumped below the mixed layer by tropical cyclones is finally released in higher latitudes. Tropical cyclones occur primarily during summer and early fall, when the ocean mixed layer is generally shallow. As the mixed layer deepens in the following winter, any warm anomaly deposited within the seasonal thermocline will be reabsorbed by the mixed layer and lost to the atmosphere. Analysis of satellite sea surface temperature and sea surface height data, together with climatological subsurface ocean data, suggests that only about one quarter of the heat that is mixed downward by tropical cyclones eventually makes it into the permanent thermocline.

No abstract