Search

The ocean mixed layer is the interface between the ocean interior and the atmosphere or sea ice and plays a key role in climate variability. It is thus critical that numerical models used in climate studies are capable of a good representation of the mixed layer, especially its depth. Here we evaluate the mixed-layer depth (MLD) in six pairs of non-eddying (1∘ grid spacing) and eddy-rich (up to 1/16∘) models from the Ocean Model Intercomparison Project (OMIP), forced by a common atmospheric state. For model evaluation, we use an updated MLD dataset computed from observations using the OMIP protocol (a constant density threshold). In winter, low-resolution models exhibit large biases in the deep-water formation regions. These biases are reduced in eddy-rich models but not uniformly across models and regions. The improvement is most noticeable in the mode-water formation regions of the Northern Hemisphere. Results in the Southern Ocean are more contrasted, with biases of either sign remaining at high resolution. In eddy-rich models, mesoscale eddies control the spatial variability in MLD in winter. Contrary to a hypothesis that the deepening of the mixed layer in anticyclones would make the MLD larger globally, eddy-rich models tend to have a shallower mixed layer at most latitudes than coarser models do. In addition, our study highlights the sensitivity of the MLD computation to the choice of a reference level and the spatio-temporal sampling, which motivates new recommendations for MLD computation in future model intercomparison projects.

The ocean mixed layer is the interface between the ocean interior and the atmosphere or sea ice and plays a key role in climate variability. It is thus critical that numerical models used in climate studies are capable of a good representation of the mixed layer, especially its depth. Here we evaluate the mixed-layer depth (MLD) in six pairs of non-eddying (1∘ grid spacing) and eddy-rich (up to 1/16∘) models from the Ocean Model Intercomparison Project (OMIP), forced by a common atmospheric state. For model evaluation, we use an updated MLD dataset computed from observations using the OMIP protocol (a constant density threshold). In winter, low-resolution models exhibit large biases in the deep-water formation regions. These biases are reduced in eddy-rich models but not uniformly across models and regions. The improvement is most noticeable in the mode-water formation regions of the Northern Hemisphere. Results in the Southern Ocean are more contrasted, with biases of either sign remaining at high resolution. In eddy-rich models, mesoscale eddies control the spatial variability in MLD in winter. Contrary to a hypothesis that the deepening of the mixed layer in anticyclones would make the MLD larger globally, eddy-rich models tend to have a shallower mixed layer at most latitudes than coarser models do. In addition, our study highlights the sensitivity of the MLD computation to the choice of a reference level and the spatio-temporal sampling, which motivates new recommendations for MLD computation in future model intercomparison projects.

The Weddell Gyre is one of the largest features of the ocean circulation adjacent to the Antarctic margins. The gyre is a dynamically complex region and participates in several processes relevant to global climate. For example, the gyre’s circulation and its strength have been linked to changes in the properties and rates of export of Antarctic Bottom Water into the global abyssal ocean. However, the dynamic controls of the Weddell Gyre’s variability are largely unknown, possibly due to the complexities of the region: the interplay of the Weddell Gyre with an overturning circulation, strong buoyancy fluxes associated with sea ice formation and melt, and open and permeable boundaries which allow for significant inflows and outflows. In this work we analyse the mechanisms controlling the Weddell Gyre’s variability using a barotropic vorticity budget of a MOM6 simulation coupled with SIS2 and forced with a repeat year 1990-91 atmospheric state derived from JRA55-do. Unlike past studies that focus on the stationary state of a control simulation, we focus on the evolution of our simulation and the response to different wind and buoyancy perturbations. We find that a balance is achieved between surface stress and bottom pressure torque, bottom drag curl and the curl of horizontal viscosity.&nbsp, The 28th IUGG General Assembly (IUGG2023) (Berlin 2023)

The abyssal ocean circulation is a key component of the global meridional overturning circulation, cycling heat, carbon, oxygen and nutrients throughout the world ocean. The strongest historical trend observed in the abyssal ocean is warming at high southern latitudes, yet it is unclear what processes have driven this warming, and whether it is linked to a slowdown in the ocean's overturning circulation. Furthermore, future change in the abyssal overturning remains uncertain, with the latest CMIP6 projections not accounting for dynamic ice-sheet melt. In this talk I will present new transient forced high-resolution coupled ocean – sea-ice model simulations to show that under a high emissions scenario, abyssal warming is set to accelerate over the next 30 years. We find that meltwater input around Antarctica drives a contraction of Antarctic Bottom Water (AABW), opening a pathway that allows warm Circumpolar Deep Water greater access to the continental shelf. The reduction in AABW formation results in warming and ageing of the abyssal ocean, consistent with recent measurements. In contrast, projected wind and thermal forcing has little impact on the properties, age, and volume of AABW. These results highlight the critical importance of Antarctic meltwater in setting the abyssal ocean overturning, with implications for global ocean biogeochemistry and climate that could last for centuries., The 28th IUGG General Assembly (IUGG2023) (Berlin 2023)

This paper contains a description of recent changes to the formulation and numerical implementation of the Quasi-Geostrophic Coupled Model (Q-GCM), which constitute a major update of the previous version of the model (Hogg et al., 2014). The Q-GCM model has been designed to provide an efficient numerical tool to study the dynamics of multi-scale midlatitude air–sea interactions and their climatic impacts. The present additions/alterations were motivated by an inquiry into the dynamics of mesoscale ocean–atmosphere coupling and, in particular, by an apparent lack of the Q-GCM atmosphere's sensitivity to mesoscale sea-surface temperature (SST) anomalies, even at high (mesoscale) atmospheric resolutions, contrary to ample theoretical and observational evidence otherwise. Major modifications aimed at alleviating this problem include an improved radiative-convective scheme resulting in a more realistic model mean state and associated model parameters; a new formulation of entrainment in the atmosphere, which prompts more efficient communication between the atmospheric mixed layer and free troposphere; and an addition of a temperature-dependent wind component in the atmospheric mixed layer and the resulting mesoscale feedbacks. The most drastic change is, however, the inclusion of moist dynamics in the model, which may be key to midlatitude ocean–atmosphere coupling. Accordingly, this version of the model is to be referred to as the MQ-GCM model. Overall, the MQ-GCM model is shown to exhibit a rich spectrum of behaviors reminiscent of many of the observed properties of the Earth's climate system. It remains to be seen whether the added processes are able to affect in fundamental ways the simulated dynamics of the midlatitude ocean–atmosphere system's coupled decadal variability.

The Southern Ocean has undergone significant climate-related changes over recent decades, including intensified westerly winds and increased radiative heating. The interplay between wind-driven cooling and radiative warming of the ocean is complex and remains unresolved. In this study, idealized wind and thermal perturbations are analyzed in a global ocean–sea ice model at two horizontal resolutions: nominally, 1° and 0.1°. The sea surface temperature (SST) response shows a clear transition from a wind-driven cooling phase to a warming phase. This warming transition is largely attributed to meridional and vertical Ekman heat advection, which are both sensitive to model resolution due to the model-dependent components of temperature gradients. At higher model resolution, due to a more accurate representation of near-surface vertical temperature inversion and upward Ekman heat advection around Antarctica, the anomalous SST warming is stronger and develops earlier. The mixed layer depth at midlatitudes initially increases due to a wind-driven increase in Ekman transport of cold dense surface water northward, but then decreases when the thermal forcing drives enhanced surface stratification; both responses are more sensitive at lower model resolution. With the wind intensification, the residual overturning circulation increases less in the 0.1° case because of the adequately resolved eddy compensation. Ocean heat subduction penetrates along more tilted isopycnals in the 1° case, but it orients to follow isopycnal layers in the 0.1° case. These findings have implications for understanding the ocean response to the combined effects of Southern Hemisphere westerly wind changes and anthropogenic warming.

The upwelling of cold water from the depths of the Southern Ocean to its surface closes the global overturning circulation and facilitates uptake of anthropogenic heat and carbon. Upwelling is often conceptualised in a zonally averaged framework as the result of isopycnal flattening via baroclinic eddies. However, upwelling is zonally non-uniform and occurs in discrete hotspots near topographic features. The mechanisms that facilitate topographically confined eddy upwelling remain poorly understood and thus limit the accuracy of parameterisations in coarse-resolution climate models.Using a high-resolution global ocean sea-ice model, we calculate spatial distributions of upwelling transport and energy conversions associated with barotropic and baroclinic instability, derived from a thickness-weighted energetics framework. We find that five major topographic hotspots of upwelling, covering less than 30% of the circumpolar longitude range, account for up to 76% of the southward eddy upwelling transport. The conversion of energy into eddies via baroclinic instability is highly spatially correlated with upwelling transport, unlike the barotropic energy conversion, which is also an order of magnitude smaller than the baroclinic conversion. This result suggests that eddy parameterisations that quantify baroclinic energy conversions could be used to improve the simulation of upwelling hotspots in climate models. We also find that eddy kinetic energy maxima are found on average 110 km downstream of upwelling hotspots in accordance with sparse observations. Our findings demonstrate the importance of localised mechanisms to Southern Ocean dynamics.

Ocean–atmosphere coupling modifies the variability of Earth’s climate over a wide range of time scales. However, attribution of the processes that generate this variability remains an outstanding problem. In this article, air–sea coupling is investigated in an eddy-resolving, medium-complexity, idealized ocean–atmosphere model. The model is run in three configurations: fully coupled, partially coupled (where the effect of the ocean geostrophic velocity on the sea surface temperature field is minimal), and atmosphere-only. A surface boundary layer temperature variance budget analysis computed in the frequency domain is shown to be a powerful tool for studying air–sea interactions, as it differentiates the relative contributions to the variability in the temperature field from each process across a range of time scales (from daily to multidecadal). This method compares terms in the ocean and atmosphere across the different model configurations to infer the underlying mechanisms driving temperature variability. Horizontal advection plays a dominant role in driving temperature variance in both the ocean and the atmosphere, particularly at time scales shorter than annual. At longer time scales, the temperature variance is dominated by strong coupling between atmosphere and ocean. Furthermore, the Ekman transport contribution to the ocean’s horizontal advection is found to underlie the low-frequency behavior in the atmosphere. The ocean geostrophic eddy field is an important driver of ocean variability across all frequencies and is reflected in the atmospheric variability in the western boundary current separation region at longer time scales.

Oceanic eddies play a profound role in mixing tracers such as heat, carbon, and nutrients, thereby regulating regional and global climate. Yet, it remains unclear how global oceanic eddy kinetic energy has evolved over the past few decades. Furthermore, coupled climate model predictions generally fail to resolve oceanic mesoscale dynamics, which could limit their accuracy in simulating future climate change. Here we show a global statistically significant increase of the eddy activity using two independent observational datasets of mesoscale variability, one directly measuring currents and the other from sea surface temperature. Regions characterized by different dynamical processes show distinct evolution in the eddy field. For example, eddy-rich regions such as boundary current extensions and the Antarctic Circumpolar Current show a significant increase of 2% and 5% per decade in eddy activity, respectively. In contrast, most of the regions of observed decrease are found in the tropical oceans. Because eddies play a fundamental role in the ocean transport of heat, momentum, and carbon, our results have far-reaching implications for ocean circulation and climate, and the modelling platforms we use to study future climate change.

The interaction of a barotropic flow with topography generates baroclinic motion that exerts a stress on the barotropic flow. Here, explicit solutions are calculated for the spatial-mean flow (i.e., the barotropic tide) resulting from a spatially uniform but time-varying body force (i.e., astronomical forcing) acting over rough topography. This approach of prescribing the force contrasts with that of previous authors who have prescribed the barotropic flow. It is found that the topographic stress, and thus the impact on the spatial-mean flow, depend on the nature of the baroclinic motion that is generated. Two types of stress are identified: (i) a “wave drag” force associated with propagating wave motion, which extracts energy from the spatial-mean flow, and (ii) a topographic “spring” force associated with standing motion at the seafloor, including bottom-trapped internal tides and propagating low-mode internal tides, which significantly damps the time-mean kinetic energy of the spatial-mean flow but extracts no energy in the time-mean. The topographic spring force is shown to be analogous to the force exerted by a mechanical spring in a forced-dissipative harmonic oscillator. Expressions for the topographic stresses appropriate for implementation as baroclinic drag parameterizations in global models are presented.

Circulation in the Southern Ocean is unique. The strong wind stress forcing and buoyancy fluxes, in concert with the lack of continental boundaries, conspire to drive the Antarctic Circumpolar Current replete with an intense eddy field. The effect of Southern Ocean eddies on the ocean circulation is significant – they modulate the momentum balance of the zonal flow, and the meridional transport of tracers and mass. The strength of the eddy field is controlled by a combination of forcing (primarily thought to be wind stress) and intrinsic, chaotic, variability associated with the turbulent flow field itself. Here, we present results from an eddy-permitting ensemble of ocean model simulations to investigate the relative contribution of forced and intrinsic processes in governing the variability of Southern Ocean eddy kinetic energy. We find that variations of the eddy field are mostly random, even on longer (interannual) timescales. Where correlations between the wind stress forcing and the eddy field exist, these interactions are dominated by two distinct timescales – a fast baroclinic instability response; and a multi-year process owing to feedback between bathymetry and the mean flow. These results suggest that understanding Southern Ocean eddy dynamics and its larger-scale impacts requires an ensemble approach to eliminate intrinsic variability, and therefore may not yield robust conclusions from observations alone. Plain Language Summary The Southern Ocean is the most turbulent part of the world’s oceans. This turbulence, often referred to as eddies, is critical to the evolution of the Southern Ocean under climate change. But it’s hard to get information about these eddies, because they occur on small scales in a large ocean basin that is poorly observed. In addition, the observational record is quite short, which makes it more difficult to use these observations to study what controls variations of these eddies. For this reason, we take an eddy-permitting ocean model, and run it 50 times with the same forcing (but a slightly different initial state). The chaotic nature of the turbulent ocean means that these model runs exhibit different evolutions. Then we use these simulations to study which eddy processes occur as a consequence of the chaotic nature of turbulence and which are forced by the external factors that are common to all model runs (such as wind stress). We conclude that monthly-to-interannual fluctuations of the Southern Ocean eddy field are dominated by chaotic processes; but that the forced variability responds to wind on particular timescales that are controlled by the mechanisms that generate ocean turbulence.

Antarctic Bottom Water (AABW), which fills the global ocean abyss, is derived from dense water that forms in several distinct Antarctic shelf regions. Previous modeling studies have reached conflicting conclusions regarding export pathways of AABW across the Southern Ocean and the degree to which AABW originating from distinct source regions are blended during their export. This study addresses these questions using passive tracer deployments in a 61-year global high-resolution (0.1°) ocean/sea-ice simulation. Two distinct export “conduits” are identified: Weddell Sea- and Prydz Bay-sourced AABW are blended together and exported mainly to the Atlantic and Indian Oceans, while Ross Sea- and Adelie Land-sourced AABW are exported mainly to the Pacific Ocean. Northward transport of each tracer occurs almost exclusively (>90%) within a single conduit. These findings imply that regional changes in AABW production may impact the three-dimensional structure of the global overturning circulation.

The response of near-Antarctic waters to freshening by increased glacial melt is investigated using a high-resolution (0.1°) global ocean–sea ice model with realistic Antarctic water-mass properties. Two meltwater perturbation experiments are conducted where the ocean model is forced with constant elevated glacial melt rates of 1.5 and 2.8 times the control rate. Within 10 years of the onset of enhanced meltwater forcing, the generation of Antarctic Bottom Water from Dense Shelf Water ceases, as shelf waters become increasingly buoyant. Increased ocean stratification triggers subsurface warming in Dense Shelf Water source regions, suggesting a localized positive feedback to melt. In a parallel response, meltwater forcing enhances the subsurface lateral density gradients of the Antarctic Slope Front that modulate the transport of warm Circumpolar Deep Water across the continental slope toward ice shelf grounding lines. Consequently, coastal freshening acts to isolate the Antarctic Ice Sheet from open ocean heat, suggesting a cooling response to melt that counteracts warming associated with stratification. Further, these strengthening density gradients accelerate westward geostrophic currents along the coast and shelf break, homogenizing shelf waters and amplifying remote feedbacks. The net effect on the continental shelf is transient warming, followed by cooling in both experiments; however, this signal is the aggregate of a complex pattern of regional warming and cooling responses. These results suggest coastal freshening by meltwater may alter the thermal forcing of the Antarctic ice sheet in ways that both accelerate and inhibit ice shelf melt at different locations along the Antarctic coastline.

Open-ocean convection is a common phenomenon that regulates mixed layer depth and ocean ventilation in the high-latitude oceans. However, many climate model simulations overestimate mixed layer depth during open-ocean convection, resulting in excessive formation of dense water in some regions. The physical processes controlling transient mixed layer depth during open-ocean convection are examined using two different numerical models: a high-resolution, turbulence-resolving nonhydrostatic model and a large-scale hydrostatic ocean model. An isolated destabilizing buoyancy flux is imposed at the surface of both models and a quasi-equilibrium flow is allowed to develop. Mixed layer depth in the turbulence-resolving and large-scale models closely aligns with existing scaling theories. However, the large-scale model has an anomalously deep mixed layer prior to quasi-equilibrium. This transient mixed layer depth bias is a consequence of the lack of resolved turbulent convection in the model, which delays the onset of baroclinic instability. These findings suggest that in order to reduce mixed layer biases in ocean simulations, parameterizations of the connection between baroclinic instability and convection need to be addressed.

Climate variability is investigated by identifying the energy sources and sinks in an idealized, coupled, ocean–atmosphere model, tuned to mimic the North Atlantic region. The spectral energy budget is calculated in the frequency domain to determine the processes that either deposit energy into or extract energy from each fluid, over time scales from one day up to 100 years. Nonlinear advection of kinetic energy is found to be the dominant source of low-frequency variability in both the ocean and the atmosphere, albeit in differing layers in each fluid. To understand the spatial patterns of the spectral energy budget, spatial maps of certain terms in the spectral energy budget are plotted, averaged over various frequency bands. These maps reveal three dynamically distinct regions: along the western boundary, the western boundary current separation, and the remainder of the domain. The western boundary current separation is found to be a preferred region to energize oceanic variability across a broad range of time scales (from monthly to decadal), while the western boundary itself acts as the dominant sink of energy in the domain at time scales longer than 50 days. This study paves the way for future work, using the same spectral methods, to address the question of forced versus intrinsic variability in a coupled climate system.

Sustained extreme temperature events in the ocean, referred to as marine heatwaves (MHWs), generate substantial ecological, social, and economic impacts. Ocean models provide insights to the drivers, persistence, and dissipation of MHWs. However, the sensitivity of MHW metrics to ocean model resolution is unknown. Here, we analyze global MHW metrics in three configurations of a global ocean‐sea ice model at coarse (1°), eddy‐permitting (0.25°), and eddy‐rich (0.1°) resolutions. We show that all configurations qualitatively represent broad‐scale global patterns of MHWs. These simulated MHWs are, however, weaker, longer‐lasting, and less frequent than in observations. The 0.1° configuration, despite local biases, performs best both globally and regionally. Based on these results, model projections of future MHW metrics using coarse‐resolution models are expected to be biased toward weaker and less frequent MHWs, when compared with results using an eddy‐rich model.

Identifying internal waves in complex flow fields is a long‐standing problem in fluid dynamics, oceanography and atmospheric science, owing to the overlap of internal waves temporal and spatial scales with other flow regimes. Lagrangian filtering—that is, temporal filtering in a frame of reference moving with the flow—is one proposed methodology for performing this separation. Here we (a) describe an improved implementation of the Lagrangian filtering methodology and (b) introduce a new freely available, parallelized Python package that applies the method. We show that the package can be used to directly filter output from a variety of common ocean models including MITgcm, Regional Ocean Modeling System and MOM5 for both regional and global domains at high resolution. The Lagrangian filtering is shown to be a useful tool to both identify (and thereby quantify) internal waves, and to remove internal waves to isolate the non‐wave flow field.

Ocean circulation and mixing regulate Earth's climate by moving heat vertically within the ocean. We present a new formalism to diagnose the role of ocean circulation and diabatic processes in setting vertical heat transport in ocean models. In this formalism we use temperature tendencies, rather than explicit vertical velocities, to diagnose circulation. Using quasi-steady-state simulations from the Australian Community Climate and Earth-System Simulator Ocean Model (ACCESS-OM2), we diagnose a diathermal overturning circulation in temperature-depth space. Furthermore, projection of tendencies due to diabatic processes onto this coordinate permits us to represent these as apparent overturning circulations. Our framework permits us to extend the concept of "Super Residual Transport,'' which combines mean and eddy advection terms with subgridscale isopycnal mixing due to mesoscale eddies but excludes small-scale threedimensional turbulent mixing effect, to construct a new overturning circulation-the "Super Residual Circulation'' (SRC). We find that in the coarse-resolution version of ACCESS-OM2 (nominally 1 degrees horizontal resolution) the SRC is dominated by an similar to 11-Sv (1 Sv [10(6) m(3) s(-1)) circulation that transports heat upward. The SRC's upward heat transport is;2 times as large in a finer-horizontal-resolution (0.1 degrees) version of ACCESS, suggesting that a differing balance of super-residual and parameterized small-scale processes may emerge as eddies are resolved. Our analysis adds new insight into superresidual processes, because the SRC elucidates the pathways in temperature and depth space along which water mass transformation occurs.

Lee waves play an important role in transferring energy from the geostrophic eddy field to turbulent mixing in the Southern Ocean. As such, lee waves can impact the Southern Ocean circulation and modulate its response to changing climate through their regulation on the eddy field and turbulent mixing. The drag effect of lee waves on the eddy field and the mixing effect of lee waves on the tracer field have been studied separately to show their importance. However, it remains unclear how the drag and mixing effects act together to modify the Southern Ocean circulation. In this study, a lee wave parameterization that includes both lee wave drag and its associated lee-wave-driven mixing is developed and implemented in an eddy-resolving idealized model of the Southern Ocean to simulate and quantify the impacts of lee waves on the Southern Ocean circulation. The results show that lee waves enhance the baroclinic transport of the Antarctic Circumpolar Current (ACC) and strengthen the lower overturning circulation. The impact of lee waves on the large-scale circulation are explained by the control of lee wave drag on isopycnal slopes through their effect on eddies, and by the control of lee-wave-driven mixing on deep stratification and water mass transformation. The results also show that the drag and mixing effects are coupled such that they act to weaken one another. The implication is that the future parameterization of lee waves in global ocean and climate models should take both drag and mixing effects into consideration for a more accurate representation of their impact on the ocean circulation.

The mesoscale eddy field plays a key role in the mixing and transport of physical and biological properties and redistribute energy budgets in the ocean. Eddy kinetic energy is commonly defined as the kinetic energy of the time-varying component of the velocity field. However, this definition contains all processes that vary in time, including coherent mesoscale eddies, jets, waves, and large-scale motions. The focus of this paper is on the eddy kinetic energy contained in coherent mesoscale eddies. We present a new method to decompose eddy kinetic energy into oceanic processes. The proposed method uses a new eddy-identification algorithm (TrackEddy). This algorithm is based on the premise that the sea level signature of a coherent eddy can be approximated as a Gaussian feature. The eddy Gaussian signature then allows for the calculation of kinetic energy of the eddy field through the geostrophic approximation. TrackEddy has been validated using synthetic sea surface height data, and then used to investigate trends of eddy kinetic energy in the Southern Ocean using Satellite Sea Surface Height anomaly (AVISO+). We detect an increasing trend of eddy kinetic energy associated with mesoscale eddies in the Southern Ocean. This trend is correlated with an increase of the coherent eddy amplitude and the strengthening of wind stress over the last two decades., Comment: 20 pages, 8 figures; TrackEddy available at https://github.com/josuemtzmo/trackeddy

Changes in ventilation of the Southern Hemisphere oceans in response to changes in midlatitude westerly winds are examined by analyzing the ideal age tracer from global eddy-permitting ocean–ice model simulations in which there is an abrupt increase and/or a meridional shift in the winds. The age response in mode and intermediate waters is found to be close to linear; the response of a combined increase and shift of peak winds is similar to the sum of the individual responses to an increase and a shift. Further, a barotropic response, following Sverdrup balance, can explain much of the age response to the changes in wind stress. There are similar peak decreases (of around 50 years) in the ideal age for a 40% increase or 2.5° poleward shift in the wind stress. However, while the age decreases throughout the thermocline for an increase in the winds, for a poleward shift in the winds the age increases in the north part of the thermocline and there are decreases in age only south of 35°S. As a consequence, the change in the volume of young water differs, with a 15% increase in the volume of water with ages younger than 50 years for a 40% increase in the winds but essentially no change in this volume for a 2.5° shift. As ventilation plays a critical role in the uptake of carbon and heat, these results suggest that the storage of anthropogenic carbon and heat in mode and intermediate waters will likely increase with a strengthening of the winds, but will be much less sensitive to a meridional shift in the peak wind stress.

Volume transport in the presence of fully resolved convection and turbulence is investigated in a reentrant channel model of the Southern Ocean. The response of the meridional overturning and zonal transport to variations in wind and buoyancy forcing is quantified. Our simulations show two overturning cells—a buoyancy‐driven lower cell and a wind‐driven upper cell. The lower overturning cell is much larger in meridional extent and magnitude than the upper cell. The mean component of the overturning transport is smaller than the fluctuating component in the lower overturning cell, indicating that transport is dominated by eddies and/or turbulent convective flow. In contrast, the upper cell can be primarily described as a mean flow (indicating minimal eddy compensation). Both cells strengthen with increasing winds, with the upper cell being more sensitive to increasing wind intensity than the lower cell. Scaling for the mean upper overturning is also derived and matches previous theories which predict a linear sensitivity of the upper cell to wind. Zonal transport remains insensitive to increasing wind stress, suggesting that the system is eddy saturated when turbulence and eddies are resolved. These results suggest that fine‐scale flows control zonal transport and abyssal meridional overturning in the Southern Ocean, while wind stress drives the upper meridional overturning transport, highlighting the importance of accurately characterizing turbulence and convection in large‐scale models.

The action of the barotropic tide over seafloor topography is the major source of internal waves at the bottom of the ocean. This internal tide has long been recognized to play an important role in ocean mixing. Here it is shown that the internal tide is also associated with a net (domain integrated) momentum flux. The net flux occurs as a result of the Doppler shifting of the internal tide at the point of generation by near-bottom mean flows. Linear theory is presented that predicts the amplitude of the wave momentum flux. The net flux scales with the bottom flow speed and the topographic wavenumber to the fourth power and is directed opposite to the bottom flow. For realistic topography, the predicted peak momentum flux occurs at scales of order 10 km and smaller, with magnitudes of order 10−3–10−2 N m−2. The theory is verified by comparison with a suite of idealized internal wave-resolving simulations. The simulations show that, for the topography considered, the wave momentum flux radiates away from the bottom and enhances mean and eddying flow when the tidal waves dissipate in the upper ocean. Our results suggest that internal tides may play an important role in forcing the upper ocean.

Observations suggest that enhanced turbulent dissipation and mixing over rough topography are modulated by the transient eddy field through the generation and breaking of lee waves in the Southern Ocean. Idealized simulations also suggest that lee waves are important in the energy pathway from eddies to turbulence. However, the energy loss from eddies due to lee wave generation remains poorly estimated. This study quantifies the relative energy loss from the time-mean and transient eddy flow in the Southern Ocean due to lee wave generation using an eddy-resolving global ocean model and three independent topographic datasets. The authors find that the energy loss from the transient eddy flow (0.12 TW; 1 TW = 1012 W) is larger than that from the time-mean flow (0.04 TW) due to lee wave generation; lee wave generation makes a larger contribution (0.12 TW) to the energy loss from the transient eddy flow than the dissipation in turbulent bottom boundary layer (0.05 TW). This study also shows that the energy loss from the time-mean flow is regulated by the transient eddy flow, and energy loss from the transient eddy flow is sensitive to the representation of anisotropy in small-scale topography. It is implied that lee waves should be parameterized in eddy-resolving global ocean models to improve the energetics of resolved flow.

Internal waves generated at the seafloor propagate through the interior of the ocean, driving mixing where they break and dissipate. However, existing theories only describe these waves in two limiting cases. In one limit, the presence of an upper boundary permits bottom-generated waves to reflect from the ocean surface back to the seafloor, and all the energy flux is at discrete wavenumbers corresponding to resonant modes. In the other limit, waves are strongly dissipated such that they do not interact with the upper boundary and the energy flux is continuous over wavenumber. Here, a novel linear theory is developed for internal tides and lee waves that spans the parameter space in between these two limits. The linear theory is compared with a set of numerical simulations of internal tide and lee wave generation at realistic abyssal hill topography. The linear theory is able to replicate the spatially-averaged kinetic energy and dissipation of even highly non-linear wave fields in the numerical simulations via an appropriate choice of the linear dissipation operator, which represents turbulent wave breaking processes.

The flow of tides over rough bathymetry in the deep ocean generates baroclinic motion including internal waves and bottom‐trapped tides. The stresses generated by this motion feedback on the amplitude and phase of the large‐scale tide. Quantifying the stresses associated with tidal flow over abyssal hills is especially important, as this scale of bathymetry is often unresolved in global baroclinic tide models, and the stresses must therefore be parameterized. Here, we extend the previous theoretical work of the authors to determine the amplitude, phasing, and vertical location of the stresses exerted on a flow driven by a time‐periodic body force when it encounters rough bathymetry. The theory compares favorably with a suite of fully nonlinear numerical simulations. It is shown that all topographic stresses are applied directly above the bathymetry, leading to a two‐layer baroclinic flow, with the near‐bottom spatial‐mean flow (the benthic tide) strongly modified by topographic stresses, and the flow at height unperturbed by the presence of topography. Our results provide a framework to improve baroclinic tide models by (i) providing a simple parameterization for the subinertial stress which is currently not included in any models, (ii) establishing that parameterized stresses should be applied in the diffusive boundary layer directly above the topography, independent of where internal tides may dissipate, and (iii) identifying a minimum resolution of ∼10 km for baroclinic tidal models to adequately capture wave resonance effects that can significantly impact the magnitude of the benthic tide.

Ocean eddies influence regional and global climate through mixing and transport of heat and properties. One of the most recognizable and ubiquitous feature of oceanic eddies are coherent vortices w...

Laboratory experiments are an important and powerful method for investigating stratified fluid dynamics. Visualising the physics of interest, however, is challenging and often limits the effectiveness of the laboratory approach. Improved techniques and methodologies that simplify visualisation will broaden the application of laboratory experiments with stratified fluid dynamics. Here, we propose a simple variation to the “double-bucket” method, commonly used to generate stratified environments for laboratory experiments, that leads to the formation of robust, discrete, uniformly-spaced double-diffusive layers that are able to be visualised with a shadowgraph. This technique, which we refer to as the “hot double-bucket” method, provides a realtime global shadowgraph view of the movement of isopycnal surfaces throughout the stratified environment. To demonstrate the utility of the hot double-bucket method we visualise the propagation of small-amplitude internal waves and compare our observations with the theoretically derived dispersion relation.

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

Numerical mixing, defined here as the physically spurious tracer diffusion due to the numerical discretization of advection, is known to contribute to biases in ocean models. However, quantifying numerical mixing is nontrivial, with most studies utilizing targeted experiments in idealized settings. Here, we present a water mass transformation‐based method for quantifying numerical mixing that can be applied to any conserved variable in general circulation models. Furthermore, the method can be applied within individual fluid columns to provide spatial information. We apply the method to a suite of global ocean model simulations with differing grid spacings and subgrid‐scale parameterizations. In all configurations numerical mixing drives diathermal heat transport of comparable magnitude to that associated with explicit parameterizations. Numerical mixing is prominent 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 it may act as a proxy for neutral diffusion when it is explicitly absent. Comparison of otherwise equivalent 1/4° and 1/10° configurations with grid‐scale dependent horizontal viscosity shows only a modest enhancement in numerical mixing at 1/4°. However, if the lateral viscosity is maintained while resolution is increased then numerical mixing is reduced by almost 35%. This result suggests that the common practice of reducing viscosity in order to maximize permitted variability must be considered carefully. Our results provide a detailed view of numerical mixing in ocean models and pave the way for improvements in parameter choices and numerical methods.

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