Search

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

The basal melt rate of Antarctica's ice shelves is largely controlled by heat delivered from the Southern Ocean to the Antarctic continental shelf. The Antarctic Slope Current (ASC) is an almost circumpolar feature that encircles Antarctica along the continental shelf break in an anti-clockwise direction. Because the circulation is to first order oriented along the topographic slope, it inhibits exchange of water masses between the Southern Ocean and the Antarctic continental shelf and thereby impacts cross-slope heat supply. Direct observations of the ASC system are sparse, but indicate a highly variable flow field both in time and space. Given the importance of the circulation near the shelf break for cross-shelf exchange of heat, it is timely to further improve our knowledge of the ASC system. This study makes use of the global ocean-sea ice model ACCESS-OM2-01 with a 1/10 degree horizontal resolution and describes the spatial and temporal variability of the velocity field. We categorise the modelled ASC into three different regimes, similar to previous works for the associated Antarctic Slope Front: (i) A surface-intensified current found predominantly in East Antarctica, (ii) a bottom-intensified current found downstream of the dense shelf water formation sites in the Ross, Weddell, and Prydz Bay Seas, and (iii) a reversed current found in West Antarctica where the eastward flowing Antarctic Circumpolar Current impinges onto the continental shelf break. We find that the temporal variability of the Antarctic Slope Current varies between the regimes. In the bottom-intensified regions, the variability is set by the timing of the dense shelf water overflows, whereas the surface-intensified flow responds to the sub-monthly variability in the wind field.

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.

In both the atmosphere and ocean, large-scale (mean) flows over topography generate internal waves. A longstanding question in both fields is what forces – often known as ‘wave drag’ – are exerted on the mean flow in this process, as such forces must be parameterized in non-wave-resolving numerical models. For a time-invariant mean flow, it is well known that lee waves are generated which extract momentum from the solid earth and deposit it where they break and dissipate at height. Here, I address the equivalent problem for a time-periodic mean flow (e.g. the ocean tide) using theory and numerical simulations. In this situation, the waves influence the amplitude and phase of the periodic mean flow near the topography regardless of where they dissipate. Dissipation plays a role in terms of controlling the magnitude of wave reflections from an upper boundary (e.g. the ocean surface) which modifies the forces acting near the topography. Our results form a framework for parameterizing tidal internal wave drag in global ocean models.

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

The meridional variability of the Subtropical Front (STF) and the drivers of variability on interannual time scales in the New Zealand region are analysed using a multi‐decadal eddy‐resolving ocean hindcast model, in comparison with Argo data. The STF marks the water mass boundary between subtropical waters and subantarctic waters, and is defined as the southern‐most location of the 11°C isotherm and 34.8 isohaline between 100 m and 500 m. The STF shifts up to 650 km (6°) meridionally on seasonal timescales. In addition to seasonal variability, shifts of around 200 km (2°) occur on interannual time scales. These shifts are connected to regional wind stress curl anomalies in the eastern Tasman Sea and east of New Zealand, which trigger Ekman convergence/divergence and Rossby waves and result in meridional transport of heat and salt into/out of the Tasman Sea. The net transports across the northern boundary of the Tasman Sea show the largest sensitivity to these wind stress curl anomalies. During periods of positive wind stress curl anomalies and Ekman convergence, the heat and salt content increases shifting the position of the STF southward. The opposite tendency occurs during periods of negative wind stress curl anomalies. The migration of the STF does not appear to be directly linked to regional climate oscillations.

On behalf of the journal, AGU, and the scientific community, the editors would like to sincerely thank those who reviewed the manuscripts for Geophysical Research Letters in 2020. The hours reading and commenting on manuscripts not only improve the manuscripts but also increase the scientific rigor of future research in the field. We particularly appreciate the timely reviews in light of the demands imposed by the rapid review process at Geophysical Research Letters. The COVID pandemic imposed additional stresses on the review process, as many reviewers had to juggle increased family commitments, hours of online meetings, remote work and instruction, lack of physical access to library resources, and other hardships to maintain the quality and timeliness of their reviews. Although we witnessed an increase in the number of submissions, the average number of days to complete a review increased by less than one day!! That says a lot about the diligence of our reviewers. We deeply appreciate their contributions in these challenging times. With the advent of AGU's data policy, many reviewers have also helped immensely to evaluate the accessibility and availability of data, and many have provided insightful comments that helped to improve the data presentation and quality. We greatly appreciate the assistance of the reviewers in advancing open science, which is a key objective of AGU's data policy. Individuals in italics provided three or more reviews for Geophysical Research Letters during the year. In total, 5,177 referees contributed to 8,786 individual reviews. Thank you again. We look forward to the coming year of exciting advances in the field and communicating those advances to our community and to the broader public.

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

International audience

On behalf of the journal, AGU, and the scientific community, the editors would like to sincerely thank those who reviewed the manuscripts for Geophysical Research Letters in 2019. The hours reading and commenting on manuscripts not only improve the manuscripts but also increase the scientific rigor of future research in the field. We particularly appreciate the timely reviews in light of the demands imposed by the rapid review process at Geophysical Research Letters. With the revival of the “major revisions” decisions, we appreciate the reviewers' efforts on multiple versions of some manuscripts. With the advent of AGU's data policy, many reviewers have helped immensely to evaluate the accessibility and availability of data associated with the papers they have reviewed, and many have provided insightful comments that helped to improve the data presentation and quality. We greatly appreciate the assistance of the reviewers in advancing open science, which is a key objective of AGU's data policy. Many of those listed below went beyond and reviewed three or more manuscripts for our journal, and those are indicated in italics.

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 recognized 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 two numerical simulations—one in the fully turbulent regime and the second an eddy‐permitting ocean model—which show that gyre‐like circulation can be driven by surface buoyancy fluxes alone. We explore this phenomenon through a combination of modeling and linear theory to highlight that the transport of ocean gyres depends upon surface buoyancy fluxes as well as wind stress.

Mid-depth heat transport toward Antarctica is focused in “cold” regions, mechanically forced by overflowing dense water., Poleward transport of warm Circumpolar Deep Water (CDW) has been linked to melting of Antarctic ice shelves. However, even the steady-state spatial distribution and mechanisms of CDW transport remain poorly understood. Using a global, eddying ocean model, we explore the relationship between the cross-slope transports of CDW and descending Dense Shelf Water (DSW). We find large spatial variability in CDW heat and volume transport around Antarctica, with substantially enhanced flow where DSW descends in canyons. The CDW and DSW transports are highly spatially correlated within ~20 km and temporally correlated on subdaily time scales. Focusing on the Ross Sea, we show that the relationship is driven by pulses of overflowing DSW lowering sea surface height, leading to net onshore CDW transport. The majority of simulated onshore CDW transport is concentrated in cold-water regions, rather than warm-water regions, with potential implications for ice-ocean interactions and global sea level rise.

Climate change in the polar regions exerts a profound influence both locally and over all of our planet. Physical and ecosystem changes influence societies and economies, via factors that include food provision, transport and access to non-renewable resources. Sea level, global climate and potentially mid-latitude weather are influenced by the changing polar regions, through coupled feedback processes, sea ice changes and the melting of snow and land-based ice sheets and glaciers.Reflecting this importance, the IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (SROCC) features a chapter highlighting past, ongoing and future change in the polar regions, the impacts of these changes, and the possible options for response. The role of the polar oceans, both in determining the changes and impacts in the polar regions and in structuring the global influence, is an important component of this chapter.With emphasis on the Southern Ocean and through comparison with the Arctic, this talk will outline key findings from the polar regions chapter of SROCC. It will synthesise the latest information on the rates, patterns and causes of changes in sea ice, ocean circulation and properties. It will assess cryospheric driving of ocean change from ice sheets, ice shelves and glaciers, and the role of the oceans in determining the past and future evolutions of polar land-based ice. The implications of these changes for climate, ecosystems, sea level and the global system will be outlined.

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

International audience; An important characteristic of geophysically turbulent flows is the transfer of energy between scales. 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 develop a method 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 kinwetic 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 doubly periodic 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. We also describe, but do not employ here, a technique for developing filters to use in complex domains.

Atmosphere and ocean are coupled via air–sea interactions. The atmospheric conditions fuel the ocean circulation and its variability, but the extent to which ocean processes can affect the atmosphere at decadal time scales remains unclear. In particular, such low-frequency variability is difficult to extract from the short observational record, meaning that climate models are the primary tools deployed to resolve this question. Here, we assess how the ocean’s intrinsic variability leads to patterns of upper-ocean heat content that vary at decadal time scales. These patterns have the potential to feed back on the atmosphere and thereby affect climate modes of variability, such as El Niño or the Interdecadal Pacific Oscillation. We use the output from a global ocean–sea ice circulation model at three different horizontal resolutions, each driven by the same atmospheric reanalysis. To disentangle the variability of the ocean’s direct response to atmospheric forcing from the variability due to intrinsic ocean dynamics, we compare model runs driven with inter-annually varying forcing (1958-2019) and model runs driven with repeat-year forcing. Models with coarse resolution that rely on eddy parameterizations, show (i) significantly reduced variance of the upper-ocean heat content at decadal time scales and (ii) differences in the spatial patterns of low-frequency variability compared with higher resolution simulations. Climate projections are typically done with general circulation models with coarse-resolution ocean components. Therefore, these biases affect our ability to predict decadal climate modes of variability and, in turn, hinder climate projections. Our results suggest that for improving climate projections, the community should move towards coupled climate models with higher oceanic resolution.

We conduct a model-data analysis of the ocean, atmosphere and terrestrial carbon system to understand their effects on atmospheric CO2 during the last glacial cycle. We use a carbon cycle box model SCP-M, combined with multiple proxy data for the atmosphere and ocean, to test for variations in ocean circulation and biological productivity across marine isotope stages spanning 130 thousand years ago to the present. The model is constrained by proxy data associated with a range of environmental conditions including sea surface temperature, salinity, ocean volume, sea ice cover and shallow water carbonate production. Model parameters for global ocean circulation, Atlantic meridional overturning circulation and Southern Ocean biological export productivity are optimised in each marine isotope stage, against proxy data for atmospheric CO2, δ13C and ∆14C and deep ocean δ13C, ∆14C and carbonate ion. Our model-data results suggest that global overturning circulation weakened at marine isotope stage 5d, coincident with a ∼ 25 ppm fall in atmospheric CO2 from the penultimate interglacial level. This change was followed by a further slowdown in Atlantic meridional overturning circulation and enhanced Southern Ocean biological export productivity at marine isotope stage 4 (∼−30 ppm). There was also a transient slowdown in Atlantic meridional overturning circulation at MIS 5b. In this model, the last glacial maximum was characterised by relatively weak global ocean and Atlantic meridional overturning circulation, and increased Southern Ocean biological export productivity (∼−15–20 ppm during MIS 2–4). Ocean circulation and Southern Ocean biology rebounded to modern values by the Holocene period. The terrestrial biosphere decreased by ∼ 500 Pg C in the lead up to the last glacial maximum, followed by a period of intense regrowth during the Holocene (∼ 750 Pg C). Slowing ocean circulation, a cooler ocean and, to a lesser extent, shallow carbonate dissolution, contributed ∼−75 ppm to atmospheric CO2 in the ∼ 100 thousand-year lead-up to the last glacial maximum, with a further ∼−10 ppm contributed during the glacial maximum. Our model results also suggest that an increase in Southern Ocean biological productivity was one of the ingredients required to achieve the last glacial maximum atmospheric CO2 level. The incorporation of longer-timescale data into quantitative ocean transport models, provides useful insights into the timing of changes in ocean processes, enhancing our understanding of the last glacial maximum and Holocene carbon cycle transition.

Antarctica has long been considered biologically isolated1. Global warming will make parts of Antarctica more habitable for invasive taxa, yet presumed barriers to dispersal—especially the Southern Ocean’s strong, circumpolar winds, ocean currents and fronts—have been thought to protect the region from non-anthropogenic colonizations from the north1,2. We combine molecular and oceanographic tools to directly test for biological dispersal across the Southern Ocean. Genomic analyses reveal that rafting keystone kelps recently travelled >20,000 km and crossed several ocean-front ‘barriers’ to reach Antarctica from mid-latitude source populations. High-resolution ocean circulation models, incorporating both mesoscale eddies and wave-driven Stokes drift, indicate that such Antarctic incursions are remarkably frequent and rapid. Our results demonstrate that storm-forced surface waves and ocean eddies can dramatically enhance oceanographic connectivity for drift particles in surface layers, and show that Antarctica is not biologically isolated. We infer that Antarctica’s long-standing ecological differences have been the result of environmental extremes that have precluded the establishment of temperate-adapted taxa, but that such taxa nonetheless frequently disperse to the region. Global warming thus has the potential to allow the establishment of diverse new species—including keystone kelps that would drastically alter ecosystem dynamics—even without anthropogenic introductions. Genomic tools and ocean circulation models show that organisms surface-drift across the Southern Ocean frequently. The extreme cold therefore keeps Antarctica biologically isolated, but as the climate warms new species may establish quickly.

Mixing efficiency is a measure of the energy lost to mixing compared to that lost to viscous dissipation. In a turbulent stratified fluid the mixing efficiency is often assumed constant at η = 0.2, whereas with convection it takes values closer to 1. The value of mixing efficiency when both stratified shear flow and buoyancy-driven convection are active remains uncertain. We use a series of numerical simulations to determine the mixing efficiency in an idealized Southern Ocean model. The model is energetically closed and fully resolves convection and turbulence such that mixing efficiency can be diagnosed. Mixing efficiency decreases with increasing wind stress but is enhanced by turbulent convection and by large thermal gradients in regions with a strongly stratified thermocline. Using scaling theory and the model results, we predict an overall mixing efficiency for the Southern Ocean that is significantly greater than 0.2 while emphasizing that mixing efficiency is not constant.

The current generation of climate models exhibit a large spread in the steady-state and projected Southern Ocean upper and lower overturning circulation, with mechanisms for deep ocean variability remaining less well understood. Here, common Southern Ocean metrics in twelve models from the Coordinated Ocean-ice Reference Experiment Phase II (CORE-II) are assessed over a 60 year period. Specifically, stratification, surface buoyancy fluxes, and eddies are linked to the magnitude of the strengthening trend in the upper overturning circulation, and a decreasing trend in the lower overturning circulation across the CORE-II models. The models evolve similarly in the upper 1 km and the deep ocean, with an almost equivalent poleward intensification trend in the Southern Hemisphere westerly winds. However, the models differ substantially in their eddy parameterisation and surface buoyancy fluxes. In general, models with a larger heat-driven water mass transformation where deep waters upwell at the surface ( ∼ 55°S) transport warmer waters into intermediate depths, thus weakening the stratification in the upper 2 km. Models with a weak eddy induced overturning and a warm bias in the intermediate waters are more likely to exhibit larger increases in the upper overturning circulation, and more significant weakening of the lower overturning circulation. We find the opposite holds for a cool model bias in intermediate depths, combined with a more complex 3D eddy parameterisation that acts to reduce isopycnal slope. In summary, the Southern Ocean overturning circulation decadal trends in the coarse resolution CORE-II models are governed by biases in surface buoyancy fluxes and the ocean density field, and the configuration of the eddy parameterisation.

Recent numerical modeling studies have suggested significant spontaneous internal wave generation near the ocean surface and energy transfers to and from these waves in the ocean interior. Spontaneous generation is the emission of waves by unbalanced, large Rossby number flows in the absence of direct forcing. Here, the authors’ previous work is extended to investigate where and how these waves exchange energy with the nonwave (mean) flow. A novel double-filtering technique is adopted to separate first the wave and nonwave fields, then the individual upward- and downward-propagating wave fields, and thereby identify the pathways of energy transfer. These energy transfers are dominated by the interaction of the waves with the vertical shear in the mean flow. Spontaneously generated waves are found to be oriented such that the downward-propagating wave is amplified by the mean shear. The internal waves propagate through the entire model depth while dissipating energy and reflect back upward. The now-upward-propagating waves have the opposite sign interaction with the mean shear and decay, losing most of their energy to the nonwave flow in the upper 500 m. Overall, in the simulations described here, approximately 30% of the wave energy is dissipated, and 70% is returned to the mean flow. The apparent preferential orientation of spontaneous generation suggests a potentially unique role for these waves in the ocean energy budget in uniformly drawing net energy from mean flow in the upper-ocean interior and transporting it to depth.

We examine the separate contributions to spurious mixing from horizontal and vertical processes in an ALE ocean model, MOM6, using reference potential energy (RPE). The RPE is a global diagnostic which changes only due to mixing between density classes. We extend this diagnostic to a sub-timestep timescale in order to individually separate contributions to spurious mixing through horizontal (tracer advection) and vertical (regridding/remapping) processes within the model. We both evaluate the overall spurious mixing in MOM6 against previously published output from other models (MOM5, MITGCM and MPAS-O), and investigate impacts on the components of spurious mixing in MOM6 across a suite of test cases: a lock exchange, internal wave propagation, and a baroclinically-unstable eddying channel. The split RPE diagnostic demonstrates that the spurious mixing in a lock exchange test case is dominated by horizontal tracer advection, due to the spatial variability in the velocity field. In contrast, the vertical component of spurious mixing dominates in an internal waves test case. MOM6 performs well in this test case owing to its quasi-Lagrangian implementation of ALE. Finally, the effects of model resolution are examined in a baroclinic eddies test case. In particular, the vertical component of spurious mixing dominates as horizontal resolution increases, an important consideration as global models evolve towards higher horizontal resolutions.

The highest rates of Antarctic glacial ice mass loss are occurring to the west of the Antarctica Peninsula in regions where warming of subsurface continental shelf waters is also largest. However, the physical mechanisms responsible for this warming remain unknown. Here we show how localized changes in coastal winds off East Antarctica can produce significant subsurface temperature anomalies (>2 °C) around much of the continent. We demonstrate how coastal-trapped barotropic Kelvin waves communicate the wind disturbance around the Antarctic coastline. The warming is focused on the western flank of the Antarctic Peninsula because the circulation induced by the coastal-trapped waves is intensified by the steep continental slope there, and because of the presence of pre-existing warm subsurface water offshore. The adjustment to the coastal-trapped waves shoals the subsurface isotherms and brings warm deep water upwards onto the continental shelf and closer to the coast. This result demonstrates the vulnerability of the West Antarctic region to a changing climate. The subsurface waters west of the Antarctic Peninsula are warming rapidly. This study shows that changes in coastal winds in East Antarctica are remotely impacting this region and drive the upwelling of warm deep water.

In the Southern Ocean, strong eastward ocean jets interact with large topographic features, generating eddies that feed back onto the mean flow. Deep-reaching eddies interact with topography, where turbulent dissipation and generation of internal lee waves play an important role in the ocean’s energy budget. However, eddy effects in the deep ocean are difficult to observe and poorly characterized. This study investigates the energy contained in eddies at depth, when an ocean jet encounters topography. This study uses a two-layer ocean model in which an imposed unstable jet encounters a topographic obstacle (either a seamount or a meridional ridge) in a configuration relevant to an Antarctic Circumpolar Current frontal jet. The authors find that the presence of topography increases the eddy kinetic energy (EKE) at depth but that the dominant processes generating this deep EKE depend on the shape and height of the obstacle as well as on the baroclinicity of the jet before it encounters topography. In cases with high topography, horizontal shear instability is the dominant source of deep EKE, while a flat bottom or a strongly sheared inflow leads to deep EKE being generated primarily through baroclinic instability. These results suggest that the deep EKE is set by an interplay between the inflowing jet properties and topography and imply that the response of deep EKE to changes in the Southern Ocean circulation is likely to vary across locations depending on the topography characteristics.

The surface mixed layer (ML) governs atmosphere–ocean fluxes, and thereby affects Earth’s climate. Accurate representation of ML processes in ocean models remains a challenge, however. The O(100) m deep ML exhibits substantial horizontal thermohaline gradients, despite being near-homogenous vertically, making it an ideal location for processes that result from the nonlinearity of the equation of state, such as cabbeling and thermobaricity. Traditional approaches to investigate these processes focus on their roles in interior water-mass transformation and are ill suited to examine their influence on the ML. However, given the climatic significance of the ML, quantifying the extent to which cabbeling and thermobaricity influence the ML density field offers insight into improving ML representations in ocean models. A recent simplified equation of state of seawater allows the local effects of cabbeling and thermobaric processes in the ML to be expressed analytically as functions of the local temperature gradient and ML depth. These simplified expressions are used to estimate the extent to which cabbeling and thermobaricity contribute to local ML density differences. These estimates compare well with values calculated directly using the complete nonlinear equation of state. Cabbeling and thermobaricity predominantly influence the ML density field poleward of 30°. Mixed layer thermobaricity is basin-scale and winter intensified, while ML cabbeling is perennial and localized to intense, zonally coherent regions associated with strong temperature fronts, such as the Antarctic Circumpolar Current and the Kuroshio and Gulf Stream Extensions. For latitudes between 40° and 50° in both hemispheres, the zonally averaged effects of ML cabbeling and ML thermobaricity can contribute on the order of 10% of the local ML density difference.

Improvements in the horizontal resolution of global ocean models, motivated by the horizontal resolution requirements for specific flow features, has advanced modelling capabilities into the dynamical regime dominated by mesoscale variability. In contrast, the choice of the vertical grid remains a subjective choice, and it is not clear that efforts to improve vertical resolution adequately support their horizontal counterparts. Indeed, considering that the bulk of the vertical ocean dynamics (including convection) are parameterized, it is not immediately obvious what the vertical grid is supposed to resolve. Here, we propose that the primary purpose of the vertical grid in a hydrostatic ocean model is to resolve the vertical structure of horizontal flows, rather than to resolve vertical motion. With this principle we construct vertical grids based on their abilities to represent baroclinic modal structures commensurate with the theoretical capabilities of a given horizontal grid. This approach is designed to ensure that the vertical grids of global ocean models complement (and, importantly, to not undermine) the resolution capabilities of the horizontal grid. We find that for z -coordinate global ocean models, at least 50 well-positioned vertical levels are required to resolve the first baroclinic mode, with an additional 25 levels per subsequent mode. High-resolution ocean-sea ice simulations are used to illustrate some of the dynamical enhancements gained by improving the vertical resolution of a 1/10° global ocean model. These enhancements include substantial increases in the sea surface height variance (∼30% increase south of 40°S), the barotropic and baroclinic eddy kinetic energies (up to 200% increase on and surrounding the Antarctic continental shelf and slopes), and the overturning streamfunction in potential density space (near-tripling of the Antarctic Bottom Water cell at 65°S).

Ocean circulation models employ ‘turbulent’ viscosity and diffusivity to represent unresolved sub-gridscale processes such as breaking internal waves. Computational power has now advanced sufficiently to permit regional ocean circulation models to be run at sufficiently high (100 m–1 km) horizontal resolution to resolve a significant part of the internal wave spectrum. Here we develop theory for boundary generated internal waves in such models, and in particular, where the waves dissipate their energy. We focus specifically on the steady lee wave problem where stationary waves are generated by a large-scale flow acting across ocean bottom topography. We generalise the energy flux expressions of [Bell, T., 1975. Topographically generated internal waves in the open ocean. J. Geophys. Res. 80, 320–327] to include the effect of arbitrary viscosity and diffusivity. Applying these results for realistic parameter choices we show that in the present generation of models with O (1) m2s − 1 horizontal viscosity/diffusivity boundary-generated waves will inevitably dissipate the majority of their energy within a few hundred metres of the boundary. This dissipation is a direct consequence of the artificially high viscosity/diffusivity, which is not always physically justified in numerical models. Hence, caution is necessary in comparing model results to ocean observations. Our theory further predicts that O ( 10 − 2 ) m2s−1 horizontal and O ( 10 − 4 ) m2s−1 vertical viscosity/diffusivity is required to achieve a qualitatively inviscid representation of internal wave dynamics in ocean models.

Recent theories, models, and observations have suggested the presence of significant spontaneous internal wave generation at density fronts near the ocean surface. Spontaneous generation is the emission of waves by unbalanced, large Rossby number flows in the absence of direct forcing. Here, spontaneous generation is investigated in a zonally reentrant channel model using parameter values typical of the Southern Ocean. The model is carefully equilibrated to obtain a steady-state wave field for which a closed energy budget is formulated. There are two main results: First, waves are spontaneously generated at sharp fronts in the top 50 m of the model. The magnitude of the energy flux to the wave field at these fronts is comparable to that from other mechanisms of wave generation. Second, the surface-generated wave field is amplified in the model interior through interaction with horizontal density gradients within the main zonal current. The magnitude of the mean-to-wave conversion in the model interior is comparable to recent observational estimates and is the dominant source of wave energy in the model, exceeding the initial spontaneous generation. This second result suggests that internal amplification of the wave field may contribute to the ocean’s internal wave energy budget at a rate commensurate with known generation mechanisms.

The ocean’s meridional overturning circulation is closed by the upwelling of dense, carbon-rich waters to the surface of the Southern Ocean. It has been proposed that upwelling in this region is driven by strong westerly winds, implying that the intensification of Southern Ocean winds in recent decades may have enhanced the rate of upwelling, potentially affecting the global overturning circulation. However, there is no consensus on the sensitivity of upwelling to winds or on the nature of the connection between Southern Ocean processes and the global overturning circulation. In this study, the sensitivity of the overturning circulation to changes in Southern Ocean westerly wind stress is investigated using an eddy-permitting ocean–sea ice model. In addition to a suite of standard circulation metrics, an energy analysis is used to aid dynamical interpretation of the model response. Increased Southern Ocean wind stress enhances the upper cell of the overturning circulation through creation of available potential energy in the Southern Hemisphere, associated with stronger upwelling of deep water. Poleward shifts in the Southern Ocean westerlies lead to a complicated transient response, with the formation of bottom water induced by increased polynya activity in the Weddell Sea and a weakening of the upper overturning cell in the Northern Hemisphere. The energetic consequences of the upper overturning cell response indicate an interhemispheric connection to the input of available potential energy in the Northern Hemisphere.

Author(s): Rajaram, H; Diffenbaugh, N; Camargo, S; Cardenas, MB; Carey, R; Cobb, K; Cory, R; Cronin, M; Dombard, A; Donohue, K; Flesch, L; Giannini, A; Hayes, G; Hogg, A; Ilyina, T; Ivanov, V; Jacobsen, S; Korte, M; Lu, G; Morlighem, M; Magnusdottir, G; Newman, A; Opher, M; Passalacqua, P; Patricola, C; Ritsema, J; Sprintall, J; Su, H; Thornton, J; Williams, P; Yau, A | Abstract: On behalf of the journal, AGU, and the scientific community, the Editors would like to sincerely thank those who reviewed manuscripts for Geophysical Research Letters in 2018. The hours reading and commenting on manuscripts not only improves the manuscripts but also increases the scientific rigor of future research in the field. We particularly appreciate the timely reviews, in light of the demands imposed by the rapid review process at Geophysical Research Letters. With the revival of the “major revisions” decisions, we appreciate the reviewers' efforts on multiple versions of some manuscripts. Many of those listed below went beyond and reviewed three or more manuscripts for our journal, and those are indicated in italics. In total, 4,484 referees contributed to 7,557 individual reviews in journal. Thank you again. We look forward to the coming year of exciting advances in the field and communicating those advances to our community and to the broader public.