Antarctic Bottom Water (AABW) formation and transport constitute a key component of the global ocean circulation. Direct observations suggest that AABW volumes and transport rates may be decreasing, but these observations are too temporally or spatially sparse to determine the cause. To address this problem, we develop a new method to reconstruct AABW transport variability using data from the GRACE (Gravity Recovery and Climate Experiment) satellite mission. We use an ocean general circulation model to investigate the relationship between ocean bottom pressure and AABW: we calculate both of these quantities in the model, and link them using a regularised linear regression. Our reconstruction from modelled ocean bottom pressure can capture 65-90% of modelled AABW transport variability, depending on the ocean basin. When realistic observational uncertainty values are added to the modelled ocean bottom pressure, the reconstruction can still capture 30-80% of AABW transport variability. Using the same regression values, the reconstruction skill is within the same range in a second, independent, general circulation model. We conclude that our reconstruction method is not unique to the model in which it was developed and can be applied to GRACE satellite observations of ocean bottom pressure. These advances allow us to create the first global reconstruction of AABW transport variability over the satellite era. Our reconstruction provides information on the interannual variability of AABW transport, but more accurate observations are needed to discern AABW transport trends.

Recent ice loss on the western Antarctic Peninsula has been driven by warming ocean waters on the continental shelf. However, due to the short observational record, our understanding of the dynamics and variability in this region remains poor. High-resolution ocean model simulations show that the temperature variability along the western Antarctic Peninsula is controlled by the rate of dense water formation in the Weddell Sea. Passive tracer advection reveals connectivity between the Weddell Sea and the coastline of the western Antarctic Peninsula and Bellingshausen Sea. During multi-year periods of weak Weddell dense water formation, dense overflow transport in the Weddell Sea decreases, while the transport of cold water around the tip of the Antarctic Peninsula strengthens, driving a temperature decrease of 0.4oC along the western Antarctic Peninsula. This mechanism implies that western Antarctic Peninsula coastal ocean temperature may cool in the future if Weddell Dense Shelf Water production slows down.

The North Atlantic Oscillation (NAO) is a leading mode of atmospheric variability, affecting the North Atlantic Ocean on sub-seasonal to multi-decadal timescales. The NAO changes the atmospheric forcing at the ocean’s surface, including winds and surface buoyancy fluxes, both of which are known to impact large-scale gyre circulation. However, the relative role of other physical processes (such as mesoscale eddies and topography) in influencing gyre circulation under NAO variability is not fully understood. Here, we analyze a series of ocean–sea ice simulations using a barotropic vorticity budget to understand long-term response of the North Atlantic gyre circulation to NAO forcing. We find that for each standard deviation increase in the NAO index, the subtropical and subpolar gyres intensify by 0.90 Sv and 3.41 Sv (1 Sv = 10⁶ m³ s⁻¹) respectively. The NAO-induced wind stress anomalies drive approximately 90\% of the change in the subtropical gyre’s interior flow. However, in the subpolar gyre’s interior, in addition to wind stress, flow-topography interactions, stratification (influenced by surface heat fluxes), and non-linear advection significantly influence the circulation. Along the western boundary the bottom pressure torque plays a key role in steering the flow, and the vorticity input by the bottom pressure torque is partly redistributed by non-linear advection. Our study highlights the importance of both atmospheric forcing and oceanic dynamical processes in driving long-term gyre circulation responses to the NAO.

The Weddell Gyre’s variability on seasonal and interannual timescales is investigated using an ocean-sea ice model at three different horizontal resolutions. The model is evaluated against available observations to demonstrate that the highest resolution configuration (0.1$^{\circ}$ in the horizontal) best reproduces observed features of the region. The simulations suggest that the gyre is subject to large variability in its circulation that is not captured by summer-biased or short-term observations. The Weddell Gyre’s seasonal cycle consists of a summer minimum and a winter maximum and accounts for changes that are between one third and a half of its mean transport . On interannual time scales we find that the gyre’s strength is correlated with the local Antarctic easterlies and that extreme events of gyre circulation are associated with changes in sea ice concentration and the characteristics of warm inflow at the eastern boundary.

The response of the ventilation of mode and intermediate waters to abrupt changes in the Southern Annular Mode (SAM) is examined by analyzing the ideal age in a global ocean-sea ice model. The age response is shown to differ between the central Pacific Ocean and other basins. In the central Pacific there are large decreases in the age of subtropical mode and intermediate waters associated with a more positive SAM, contrasting only small age changes in the Atlantic and Indian Oceans, except near where intermediate water density surfaces outcrop. These interbasin differences hold for simulations at different horizontal resolutions, and can be explained by the combination of zonal variations in wind stress changes associated with the SAM, and differences in the age response to an increase or shift in the wind stress. These results suggest that the carbon and heat uptake associated with the SAM will likely vary between ocean basins.

Midlatitude gyres in the ocean are large scale horizontal circulations that are intensified on the western boundary of the ocean, giving rise to currents such as the Gulf Stream. The physical mechanism underlying gyres is widely recognised to involve the curl of the wind stress, which injects potential vorticity into the upper ocean. However, model results have highlighted the role of surface buoyancy fluxes (principally heating and cooling of the ocean surface) in driving circulation and enhancing gyre variability. Here we present results from numerical simulations in the fully turbulent regime which show that gyre-like circulation can be driven by surface buoyancy fluxes alone. We explore this phenomenon through a combination of modelling and linear theory to highlight that the transport of ocean gyres depends upon surface buoyancy fluxes as well as wind stress. Thus, the strength of gyres may be influenced by surface warming in response to climate change.

An important characteristic of geophysically turbulent flows is the transfer of energy between scales. It is expected that balanced flows pass energy from smaller to larger scales as part of the well-known upscale cascade while submesoscale and smaller scale flows can transfer energy eventually to smaller, dissipative scales. Much effort has been put into quantifying these transfers, but a complicating factor in realistic settings is that the underlying flows are often strongly spatially heterogeneous and anisotropic. Furthermore, the flows may be embedded in irregularly shaped domains that can be multiply connected. As a result, straightforward approaches like computing Fourier spatial spectra of nonlinear terms suffer from a number of conceptual issues. In this paper, we endeavor to compute cross-scale energy transfers in general settings, allowing for arbitrary flow structure, anisotropy and inhomogeneity. We employ a Green's function approach to the kinetic energy equation to relate kinetic energy at a point to its Lagrangian history. A spatial filtering of the resulting equation naturally decomposes kinetic energy into length scale dependent contributions and describes how the transfer of energy between those scales takes place. The method is applied to a numerical simulation of vortex merger, resulting in the demonstration of the expected upscale energy cascade. Somewhat novel results are that the energy transfers are dominated by pressure work, rather than kinetic energy exchange, and dissipation is a noticeable influence on the larger scale energy budgets.

Numerical mixing, defined here as the physically spurious diffusion of tracers due to the numerical discretization of advection, is known to contribute to biases in ocean circulation models. However, quantifying numerical mixing is non-trivial, with most studies utilizing specifically targeted experiments in idealized settings. Here, we present a precise, online water-mass transformation-based method for quantifying numerical mixing that can be applied to any conserved variable in global general circulation models. Furthermore, the method can be applied within individual fluid columns to provide a spatially-resolved metric. We apply the method to a suite of global ocean-sea ice model simulations with differing grid spacings and sub-grid scale parameterizations. In all configurations numerical mixing drives across-isotherm heat transport of comparable magnitude to that associated with explicitly-parameterized mixing. Numerical mixing is prominent at warm temperatures in the tropical thermocline, where it is sensitive to the vertical diffusivity and resolution. At colder temperatures, numerical mixing is sensitive to the presence of explicit neutral diffusion, suggesting that much of the numerical mixing in these regions acts as a proxy for neutral diffusion when it is explicitly absent. Comparison of equivalent (with respect to vertical resolution and explicit mixing parameters) $1/4^\circ$ and $1/10^\circ$ horizontal resolution configurations shows only a modest enhancement in numerical mixing at $1/4^\circ$. Our results provide a detailed view of numerical mixing in ocean models and pave the way for future improvements in numerical methods.

The positive trend of the Southern Annular Mode (SAM) will impact the Southern Ocean’s role in Earth’s climate, however the details of the Southern Ocean’s response remain uncertain. We introduce a methodology to examine the influence of SAM on the Southern Ocean, and apply this method to a global ocean–sea-ice model run at three resolutions (1$^{\circ}$, 1/4$^{\circ}$ and 1/10$^{\circ}$). Our methodology drives perturbation simulations with realistic atmospheric forcing of extreme SAM conditions. The thermal response agrees with previous studies; positive SAM perturbations warm the upper ocean north of the windspeed maximum and cool it to the south, with the opposite response for negative SAM. The overturning circulation exhibits a rapid response that increases/decreases for positive/negative SAM perturbations and is insensitive to model resolution. The longer term adjustment of the overturning circulation, however, depends on the representation of eddies, and is faster at higher resolutions.

Wind is an important driver of large-scale ocean currents, imparting momentum into the ocean at the sea surface. In particular, strong westerly winds help to drive the Antarctic Circumpolar Current, which of key importance for the global climate system. Over the past decades observations established that the strength of the westerlies over the Southern Ocean has increased as a result of climate change forcing. This increase is consistent with global climate model simulations. The future climate state depends strongly on how will the Antarctic Circumpolar Current respond to this strengthening. Eddy saturation is a theoretical regime where the transport of the current remains insensitive to the strengthening of the westerlies. Instead, the strengthening of the westerlies energizes transient eddies. Both satellite observations and numerical simulations suggest that the Antarctic Circumpolar Current is close to the eddy saturated limit. Traditionally eddy saturation has been attributed to baroclinic processes, but recent work suggests that barotropic processes that involve, e.g., standing meanders of the Antarctic Circumpolar Current, can also be responsible for producing eddy-saturated states. Here, we discus the different physical entities of the“usual” baroclinic eddy saturation as well as the recent notion of barotropic eddy saturation. We assess the relative importance of barotropic and baroclinic processes in producing eddy-saturated states using numerical simulations of primitive equations in an idealized setup. Lastly, we discuss potential implications these processes have on global ocean modeling.

We introduce a novel wind-derived metric to quantify variability in the Southern Ocean overturning circulation. This metric, which we call the Ekman streamfunction, integrates the Ekman pumping vertical velocity zonally and northwards from the Antarctic coastline to a given latitude. To evaluate the relationship between the Ekman streamfunction and Southern Ocean overturning circulation, we use a global 0.1 ocean–sea-ice model driven with interannual forcing (1958-2018). Throughout much of the Southern Ocean, strong correlations (r>0.9) exist between the Ekman streamfunction and the Southern Ocean overturning circulation on monthly and annual timescales. A regression analysis identifies regions where Ekman streamfunction variability coincides with >4Sv changes in the overturning; one such location is where the wind stress curl changes sign and the Ekman pumping velocity is highly variable. This study offers a new approach to infer recent changes in the Southern Ocean overturning circulation from existing datasets of wind stress.

In the Southern Ocean, mesoscale eddies contribute to the upwelling of deep waters along sloping isopycnals, helping to close the upper branch of the meridional overturning circulation. Eddy energy is not uniformly distributed along the Antarctic Circumpolar Current (ACC). Instead, ‘hotspots’ of eddy energy that are associated with enhanced eddy-induced upwelling exist downstream of topographic features. This study shows that, in idealized eddy-resolved simulations, a topographic feature in the ACC path can enhance and localize eddy-induced upwelling. However, the upwelling systematically occurs in regions where eddies grow through baroclinic instability, rather than in regions where eddy energy is large. Across a range of parameters, along-stream eddy growth rate is a more reliable indicator of eddy upwelling than traditional parameterizations such as eddy kinetic energy, eddy potential energy or isopycnal slope. Ocean eddy parameterizations should consider metrics specific to the growth of baroclinic instability to accurately model eddy upwelling near topography.