Search

Wind-generated near-inertial internal waves (NIWs) are triggered in the mixed layer and propagate down into the ocean interior. Observational and numerical studies have shown the effects of background vorticity and high shear on propagating NIWs. However, the impacts of the background mean flow on NIWs as a function of the waves’ horizontal wavelength have yet to be fully investigated. Here, two distinct cases are analyzed, namely, the propagation of wind-generated, large-scale NIWs in negative vorticity and the behavior of small-scale NIWs in high shear. The propagation and energetics of the respective NIWs are investigated using a realistic eddy-resolving numerical simulation of the Kuroshio region. The large-scale NIWs display a rapid vertical propagation to depth in negative vorticity areas, while the small-scale NIWs are confined to shallower depths in high-shear regions. Furthermore, the dominant energy sources and sinks of near-inertial energy are estimated as the respective NIWs propagate into the ocean’s interior. The qualitative analysis of NIW energetics reveals that the wind triggers the generation of both the large-scale and small-scale NIWs, but the waves experience further amplification as they draw energy from the background mean flow upon propagation in negative vorticity and high-shear regions, respectively. In addition, the study demonstrates that small-scale NIWs can be induced independently by wind fluctuations and do not necessarily rely on straining nor refraction of large-scale NIWs by mesoscale motions.

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 export of Antarctic Bottom Water (AABW) supplies the bottom cell of the global overturning circulation and plays a key role in regulating climate. This AABW outflow must cross, and is therefore mediated by, the Antarctic Circumpolar Current (ACC). Previous studies present widely-varying conceptions of the role of the ACC in directing AABW across the Southern Ocean, suggesting either that AABW may be zonally recirculated by the ACC, or that AABW may flow northward within deep western boundary currents (DWBC) against bathymetry. In this study the authors investigate how the forcing and geometry of the ACC influences the transport and transformation of AABW using a suite of process-oriented model simulations. The model exhibits a strong dependence on the elevation of bathymetry relative to AABW layer thickness: higher meridional ridges suppress zonal AABW exchange, increase the strength of flow in the DWBC, and reduce the meridional variation in AABW density across the ACC. Furthermore, the transport and transformation vary with density within the AABW layer, with denser varieties of AABW being less efficiently transported between basins. These findings indicate that changes in the thickness of the AABW layer, for example due to changes in Antarctic shelf processes, and tectonic changes in the sea floor shape may alter the pathways and transformation of AABW across the ACC.

The Antarctic Slope Current (ASC) circumnavigates the Antarctic continent following the continental slope and separating the waters on the continental shelf from the deeper offshore Southern Ocean. Water mass exchanges across the continental slope are critical for the global climate as they impact the global overturning circulation and the mass balance of the Antarctic ice sheet via basal melting. Despite the ASC’s global importance, little is known about its spatial and subannual variability, as direct measurements of the velocity field are sparse. Here, we describe the ASC in a global eddying ocean–sea ice model and reveal its large-scale spatial variability by characterizing the continental slope using three regimes: the surface-intensified ASC, the bottom-intensified ASC, and the reversed ASC. Each ASC regime corresponds to a distinct classification of the density field as previously introduced in the literature, suggesting that the velocity and density fields are governed by the same leading-order dynamics around the Antarctic continental slope. Only the surface-intensified ASC regime has a strong seasonality. However, large temporal variability at a range of other time scales occurs across all regimes, including frequent reversals of the current. We anticipate our description of the ASC’s spatial and subannual variability will be helpful to guide future studies of the ASC aiming to advance our understanding of the region’s response to a changing climate.

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.

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.4C 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 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 Weddell Gyre, located in the Weddell Sea is one of the southernmost open ocean reaches in the world and largest feature of the ocean circulation south of the Antarctic Circumpolar Current. It is adjacent to a major site of bottom water production in the southwestern Weddell Sea and participates in poleward heat transport via its cyclonic circulation that brings relatively warm waters south towards the Antarctic continent. The region is covered by sea ice through most of the year, which has historically prevented long, continuous observational efforts, both in situ and remote. As a result, ocean circulation models offer perhaps the best means of estimating the Weddell Gyre's variability. Using a coupled ocean/sea ice high-resolution global model, ACCESS-OM2 at 0.1∘ horizontal resolution, we assess the variability of the Weddell Gyre on seasonal - multi-decadal timescales and explore possible drivers of this variability. The simulations suggest that the gyre is exhibits large variability in its circulation that is not captured by summer-biased or short-term observations. Anomalous strong and weak periods of the gyre's circulation are linked to changes in sea ice concentration and other oceanic features of the region. We further explore the link to possible driving mechanisms, including surface stress forcing and surface buoyancy fluxes.

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.

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

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 ventilation is the transfer of tracers and young water from the surface down into the ocean interior. The tracers that can be transported to depth include anthropogenic heat and carbon, both of which are critical to understanding future climate trajectories. Ventilation occurs in both high- and midlatitude regions, but it is the southern midlatitudes that are responsible for the largest fraction of anthropogenic heat and carbon uptake; such Southern Ocean ventilation is the focus of this review. Southern Ocean ventilation occurs through a chain of interconnected mechanisms, including the zonally averaged meridional overturning circulation, localized subduction, eddy-driven mixing along isopycnals, and lateral transport by subtropical gyres. To unravel the complex pathways of ventilation and reconcile conflicting results, here we assess the relative contribution of each of thesemechanisms, emphasizing the three-dimensional and temporally varying nature of the ventilation of the Southern Ocean pycnocline. We conclude that Southern Ocean ventilation depends on multiple processes and that simplified frameworks that explain ventilation changes through a single process are insufficient.

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

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