Publications

2021
Stewart AL, Chi X, Solodoch A, Hogg AM. High-frequency fluctuations in Antarctic Bottom Water transport driven by Southern Ocean winds. Geophysical Research Letters [Internet]. 2021;8 :e2021GL094569. Publisher's VersionAbstract
Northward flow of Antarctic Bottom Water (AABW) across the Southern Ocean comprises a key component of the global overturning circulation. Yet AABW transport remains poorly constrained by observations and state estimates, and there is presently no means of directly monitoring any component of the Southern Ocean overturning. However, AABW flow is dynamically linked to Southern Ocean surface circulation via the zonal momentum balance, offering potential routes to indirect monitoring of the transport. Exploiting this dynamical link, this study shows that wind stress (WS) fluctuations drive large AABW transport fluctuations on time scales shorter than grl62918-math-00012 years, which comprise almost all of the transport variance. This connection occurs due to differing time scales on which topographic and interfacial form stresses respond to wind variability, likely associated with differences in barotropic versus baroclinic Rossby wave propagation. These findings imply that AABW transport variability can largely be reconstructed from the surface WS alone.
Download PDF
Stewart AL. Mesoscale, tidal and seasonal/interannual drivers of the Weddell Sea overturning circulation. Journal of Physical Oceanography [Internet]. 2021;51 :3695–3722. Publisher's VersionAbstract
The Weddell Sea supplies 40–50% of the Antarctic BottomWaters that fill the global ocean abyss, and therefore exerts significant influence over global circulation and climate. Previous studies have identified a range of different processes that may contribute to dense shelf water (DSW) formation and export on the southern Weddell Sea continental shelf. However, the relative importance of these processes has not been quantified, which hampers prioritization of observational deployments and development of model parameterizations in this region. In this study a high-resolution (1/12°) regional model of the southern Weddell Sea is used to quantify the overturning circulation and decompose it into contributions due to multi-annual mean flows, seasonal/interannual variability, tides, and other sub-monthly variability. It is shown that tides primarily influence the overturning by changing the melt rate of the Filchner-Ronne Ice Shelf (FRIS). The resulting ~0.2 Sv decrease in DSW transport is comparable to the magnitude of the overturning in the FRIS cavity, but small compared to DSW export across the continental shelf break. Seasonal/interannual fluctuations exert a modest influence on the overturning circulation due to the relatively short (8-year) analysis period. Analysis of the transient energy budget indicates that the non-tidal, sub-monthly variability is primarily baroclinically-generated eddies associated with dense overflows. These eddies play a comparable role to the mean flow in exporting dense shelf waters across the continental shelf break, and account for 100% of the transfer of heat onto the continental shelf. The eddy component of the overturning is sensitive to model resolution, decreasing by a factor of ~2 as the horizontal grid spacing is refined from 1/3° to 1/12°.
Download PDF
Jagannathan A, Srinivasan K, McWilliams JC, Molemaker MJ, Stewart AL. Boundary layer-mediated vorticity generation in currents over sloping bathymetry. Journal of Physical Oceanography [Internet]. 2021;51 :1757-1778. Publisher's VersionAbstract
Current–topography interactions in the ocean give rise to eddies spanning a wide range of spatial and temporal scales. The latest modeling efforts indicate that coastal and underwater topography are important generation sites for submesoscale coherent vortices (SCVs), characterized by horizontal scales of O(0.1−10)km. Using idealized, submesoscale and bottom boundary layer (BBL)-resolving simulations and adopting an integrated vorticity balance formulation, we quantify precisely the role of BBLs in the vorticity generation process. In particular, we show that vorticity generation on topographic slopes is attributable primarily to the torque exerted by the vertical divergence of stress at the bottom. We refer to this as the bottom stress divergence torque (BSDT). BSDT is a fundamentally nonconservative torque that appears as a source term in the integrated vorticity budget and is to be distinguished from the more familiar bottom stress curl (BSC). It is closely connected to the bottom pressure torque (BPT) via the horizontal momentum balance at the bottom and is in fact shown to be the dominant component of BPT in solutions with a well-resolved BBL. This suggests an interpretation of BPT as the sum of a viscous, vorticity-generating component (BSDT) and an inviscid, “flow-turning” component. Companion simulations without bottom drag illustrate that although vorticity generation can still occur through the inviscid mechanisms of vortex stretching and tilting, the wake eddies tend to have weaker circulation, be substantially less energetic, and have smaller spatial scales.
Download PDF
Bai Y, Wang Y, Stewart AL. Does Topographic Form Stress Impede Prograde Ocean Currents?. Journal of Physical Oceanography [Internet]. 2021;51 :2617-2638. Publisher's VersionAbstract
Topographic form stress (TFS) plays a central role in constraining the transport of the Antarctic Circumpolar Current (ACC), and thus the rate of exchange between the major ocean basins. Topographic form stress generation in the ACC has been linked to the formation of standing Rossby waves, which occur because the current is retrograde (opposing the direction of Rossby wave propagation). However, it is unclear whether TFS similarly retards current systems that are prograde (in the direction of Rossby wave propagation), which cannot arrest Rossby waves. An isopycnal model is used to investigate the momentum balance of wind-driven prograde and retrograde flows in a zonal channel, with bathymetry consisting of either a single ridge or a continental shelf and slope with a meridional excursion. Consistent with previous studies, retrograde flows are almost entirely impeded by TFS, except in the limit of flat bathymetry, whereas prograde flows are typically impeded by a combination of TFS and bottom friction. A barotropic theory for standing waves shows that bottom friction serves to shift the phase of the standing wave’s pressure field from that of the bathymetry, which is necessary to produce TFS. The mechanism is the same in prograde and retrograde flows, but is most efficient when the mean flow arrests a Rossby wave with a wavelength comparable to that of the bathymetry. The asymmetry between prograde and retrograde momentum balances implies that prograde current systems may be more sensitive to changes in wind forcing, for example associated with climate shifts.
Download PDF
Cohanim K, Zhao KX, Stewart AL. Dynamics of Eddies Generated by Sea Ice Leads. Journal of Physical Oceanography [Internet]. 2021;51 :3071–3092. Publisher's VersionAbstract
Interaction between the atmosphere and ocean in sea ice–covered regions is largely concentrated in leads, which are long, narrow openings between sea ice floes. Refreezing and brine rejection in these leads inject salt that plays a key role in maintaining the polar halocline. The injected salt forms dense plumes that subsequently become baroclinically unstable, producing submesoscale eddies that facilitate horizontal spreading of the salt anomalies. However, it remains unclear which properties of the stratification and leads most strongly influence the vertical and horizontal spreading of lead-input salt anomalies. In this study, the spread of lead-injected buoyancy anomalies by mixed layer and eddy processes are investigated using a suite of idealized numerical simulations. The simulations are complemented by dynamical theories that predict the plume convection depth, horizontal eddy transfer coefficient, and eddy kinetic energy as functions of the ambient stratification and lead properties. It is shown that vertical penetration of buoyancy anomalies is accurately predicted by a mixed layer temperature and salinity budget until the onset of baroclinic instability (~3 days). Subsequently, these buoyancy anomalies are spread horizontally by eddies. The horizontal eddy diffusivity is accurately predicted by a mixing-length scaling, with a velocity scale set by the potential energy released by the sinking salt plume and a length scale set by the deformation radius of the ambient stratification. These findings indicate that the intermittent opening of leads can efficiently populate the polar halocline with submesoscale coherent vortices with diameters of ~10 km, and they provide a step toward parameterizing their effect on the horizontal redistribution of salinity anomalies.
Download PDF
Solodoch A, Stewart AL, McWilliams JC. Formation of anticyclones above topographic depressions. Journal of Physical Oceanography [Internet]. 2021;51 :207-228. Publisher's VersionAbstract
Long-lived anticyclonic eddies (ACs) have been repeatedly observed over several North Atlantic basins characterized by bowl-like topographic depressions. Motivated by these previous findings, the authors conduct numerical simulations of the spindown of eddies initialized in idealized topographic bowls. In experiments with one or two isopycnal layers, it is found that a bowl-trapped AC is an emergent circulation pattern under a wide range of parameters. The trapped AC, often formed by repeated mergers of ACs over the bowl interior, is characterized by anomalously low potential vorticity (PV). Several PV segregation mechanisms that can contribute to the AC formation are examined. In one-layer experiments, the dynamics of the AC are largely determined by a nonlinearity parameter ϵ that quantifies the vorticity of the AC relative to the bowl’s topographic PV gradient. The AC is trapped in the bowl for low ϵ≲1, but for moderate values (0.5≲ϵ≲1) partial PV segregation allows the AC to reside at finite distances from the center of the bowl. For higher ϵ≳1, eddies freely cross the topography and the AC is not confined to the bowl. These regimes are characterized across a suite of model experiments using ϵ and a PV homogenization parameter. Two-layer experiments show that the trapped AC can be top or bottom intensified, as determined by the domain-mean initial vertical energy distribution. These findings contrast with previous theories of mesoscale turbulence over topography that predict the formation of a prograde slope current, but do not predict a trapped AC.
Download PDF
Stewart AL, McWilliams JC, Solodoch A. On the role of bottom pressure torques in wind-driven gyres. Journal of Physical Oceanography [Internet]. 2021;51 :1441-1464. Publisher's VersionAbstract
Previous studies have concluded that the wind-input vorticity in ocean gyres is balanced by bottom pressure torques (BPT), when integrated over latitude bands. However, the BPT must vanish when integrated over any area enclosed by an isobath. This constraint raises ambiguities regarding the regions over which BPT should close the vorticity budget, and implies that BPT generated to balance a local wind stress curl necessitates the generation of a compensating, nonlocal BPT and thus nonlocal circulation. This study aims to clarify the role of BPT in wind-driven gyres using an idealized isopycnal model. Experiments performed with a single-signed wind stress curl in an enclosed, sloped basin reveal that BPT balances the winds only when integrated over latitude bands. Integrating over other, dynamically motivated definitions of the gyre, such as barotropic streamlines, yields a balance between wind stress curl and bottom frictional torques. This implies that bottom friction plays a nonnegligible role in structuring the gyre circulation. Nonlocal bottom pressure torques manifest in the form of along-slope pressure gradients associated with a weak basin-scale circulation, and are associated with a transition to a balance between wind stress and bottom friction around the coasts. Finally, a suite of perturbation experiments is used to investigate the dynamics of BPT. To predict the BPT, the authors extend a previous theory that describes propagation of surface pressure signals from the gyre interior toward the coast along planetary potential vorticity contours. This theory is shown to agree closely with the diagnosed contributions to the vorticity budget across the suite of model experiments.
Download PDF
Moscoso JE, Stewart AL, Bianchi D, McWilliams JC. The Meridionally Averaged Model of Eastern Boundary Upwelling Systems (MAMEBUSv1.0). Geoscientific Model Development [Internet]. 2021;14 :763-794. Publisher's VersionAbstract
Eastern boundary upwelling systems (EBUSs) are physically and biologically active regions of the ocean with substantial impacts on ocean biogeochemistry, ecology, and global fish catch. Previous studies have used models of varying complexity to study EBUS dynamics, ranging from minimal two-dimensional (2-D) models to comprehensive regional and global models. An advantage of 2-D models is that they are more computationally efficient and easier to interpret than comprehensive regional models, but their key drawback is the lack of explicit representations of important three-dimensional processes that control biology in upwelling systems. These processes include eddy quenching of nutrients and meridional transport of nutrients and heat. The authors present the Meridionally Averaged Model of Eastern Boundary Upwelling Systems (MAMEBUS) that aims at combining the benefits of 2-D and 3-D approaches to modeling EBUSs by parameterizing the key 3-D processes in a 2-D framework. MAMEBUS couples the primitive equations for the physical state of the ocean with a nutrient– phytoplankton–zooplankton–detritus model of the ecosystem, solved in terrain-following coordinates. This article defines the equations that describe the tracer, momentum, and biological evolution, along with physical parameterizations of eddy advection, isopycnal mixing, and boundary layer mixing. It describes the details of the numerical schemes and their implementation in the model code, and provides a reference solution validated against observations from the California Current. The goal of MAMEBUS is to facilitate future studies to efficiently explore the wide space of physical and biogeochemical parameters that control the zonal variations in EBUSs.
Download PDF
Zhao KX, Stewart AL, McWilliams JC. Geometric Constraints on Glacial Fjord-Shelf Exchange. Journal of Physical Oceanography [Internet]. 2021;51 :1223-1246. Publisher's VersionAbstract
The oceanic connections between tidewater glaciers and continental shelf waters are modulated and controlled by geometrically complex fjords. These fjords exhibit both overturning circulations and horizontal recirculations, driven by a combination of water mass transformation at the head of the fjord, variability on the continental shelf, and atmospheric forcing. However, it remains unclear which geometric and forcing parameters are the most important in exerting control on the overturning and horizontal recirculation. To address this, idealized numerical simulations are conducted using an isopycnal model of a fjord connected to a continental shelf, which is representative of regions in Greenland and the West Antarctic Peninsula. A range of sensitivity experiments demonstrate that sill height, wind direction/strength, subglacial discharge strength, and depth of offshore warm water are of first-order importance to the overturning circulation, while fjord width is also of leading importance to the horizontal recirculation. Dynamical predictions are developed and tested for the overturning circulation of the entire shelf-to-glacierface domain, subdivided into three regions: the continental shelf extending from the open ocean to the fjord mouth, the sill-overflow at the fjord mouth, and the plume-driven water mass transformation at the fjord head. A vorticity budget is also developed to predict the strength of the horizontal recirculation, which provides a scaling in terms of the overturning and bottom friction. Based on these theories, we may predict glacial melt rates that take into account overturning and recirculation, which may be used to refine estimates of ocean-driven melting of the Greenland and Antarctic ice sheets.
Download PDF
2020
McCoy DE, Bianchi D, Stewart AL. Global Observations of Submesoscale Coherent Vortices in the Ocean. Progress in Oceanography [Internet]. 2020;189 :102452. Publisher's VersionAbstract
Subsurface-intensified anticyclones are ubiquitous in the ocean, yet their impact on the large-scale transport of heat, salt and chemical tracers is poorly understood. These submesoscale coherent vortices (SCVs) can trap and advect waters thousands of kilometers away from the formation region, providing a transport pathway that is unresolved by low-resolution Earth System Models. However, knowledge of the importance of these eddies for the large scale circulation is hindered by the lack of systematic observations. Here, we take advantage of the global network of Argo floats to identify occurrences of these eddies, which appear as weakly stratified anomalous water masses with Gaussian-shaped vertical structures. We develop a general algorithm to detect subsurface eddies that have propagated away from their source region, and apply it to the database of Argo float profiles, resulting in roughly 4000 detections from more than 20 years of observations. We further group detections into regional populations to identify hot-spots of generation and mechanisms of formation. Analysis of regional SCV statistics reveals important sites of SCV generation in Eastern Boundary Upwelling Systems, marginal sea overflows, and mode water formation regions along major open-ocean fronts. Because of the heat and salt anomaly contained within their cores, SCV could leave a significant imprint on the hydrographic properties of water masses in regions of high SCV density.
Download PDF
Hazel JE, Stewart AL. Bi-stability of the Filchner-Ronne Ice Shelf Cavity Circulation and Basal Melt. Journal of Geophysical Research: Oceans [Internet]. 2020;125 :e2019JC015848. Publisher's VersionAbstract
Circulation and water mass transformation within the Filchner-Ronne Ice Shelf (FRIS) cavity create precursors to Antarctic bottom water, which closes the global overturning circulation. This water mass transformation is contingent upon a relative low rate of FRIS basal melt, currently around 100–200 Gt/yr. Previous studies have indicated that Antarctic climate changes may induce intrusions of warm modified Warm Deep Water (mWDW) and an order-of-magnitude increase in basal melt, and signatures of mWDW have recently been observed along the face of the FRIS. However, it remains unclear how changes in near-Antarctic climate translate mechanistically to changes in mWDW access to the FRIS cavity. In this study a regional model is developed to investigate FRIS circulation dependence on local atmospheric state. Experiments with modified initial cavity conditions but identical atmospheric states yield bistable “warm” and “cold” FRIS cavity states, with an order-of-magnitude difference in basal melt rates. Idealized atmospheric perturbation experiments reveal that relatively modest perturbations to the katabatic winds shift the FRIS cavity between “warm” and “cold” states, which occur when the FRIS cavity is filled by mWDW or High Salinity Shelf Water (HSSW), respectively. The authors present a conceptual model in which the FRIS cavity state is determined by whether mWDW or HSSW is denser and thus floods the cavity; these states are bistable because the basal melt rate feeds back on the salinity of HSSW. These findings highlight a key role for the katabatic winds in mediating the melt of the FRIS and other Antarctic ice shelves
Download PDF
Solodoch A, McWilliams JC, Stewart AL, Gula J, Renault L. Why Does the Deep Western Boundary Current "Leak" Around Flemish Cap?. Journal of Physical Oceanography [Internet]. 2020;50 :1989-2016. Publisher's VersionAbstract
The southward-flowing deep limb of the Atlantic meridional overturning circulation is composed of both the deep western boundary current (DWBC) and interior pathways. The latter are fed by “leakiness” from the DWBC in the Newfoundland Basin. However, the cause of this leakiness has not yet been explored mechanistically. Here the statistics and dynamics of the DWBC leakiness in the Newfoundland Basin are explored using two float datasets and a high-resolution numerical model. The float leakiness around Flemish Cap is found to be concentrated in several areas (hot spots) that are collocated with bathymetric curvature and steepening. Numerical particle advection experiments reveal that the Lagrangian mean velocity is offshore at these hot spots, while Lagrangian variability is minimal locally. Furthermore, model Eulerian mean streamlines separate from the DWBC to the interior at the leakiness hot spots. This suggests that the leakiness of Lagrangian particles is primarily accomplished by an Eulerian mean flow across isobaths, though eddies serve to transfer around 50% of the Lagrangian particles to the leakiness hot spots via chaotic advection, and rectified eddy transport accounts for around 50% of the offshore flow along the southern face of Flemish Cap. Analysis of the model’s energy and potential vorticity budgets suggests that the flow is baroclinically unstable after separation, but that the resulting eddies induce modest modifications of the mean potential vorticity along streamlines. These results suggest that mean uncompensated leakiness occurs mostly through inertial separation, for which a scaling analysis is presented. Implications for leakiness of other major boundary current systems are discussed.
Download PDF
Sun S, Eisenman I, Zanna L, Stewart AL. Surface constraints on the depth of the Atlantic Meridional Overturning Circulation: Southern Ocean vs North Atlantic. Journal of Climate [Internet]. 2020;33 :3125-3149. Publisher's VersionAbstract
Paleoclimate proxy evidence suggests that the Atlantic meridional overturning circulation (AMOC) was about 1000 m shallower at the Last Glacial Maximum (LGM) compared to the present. Yet it remains unresolved what caused this glacial shoaling of the AMOC, and many climate models instead simulate a deeper AMOC under LGM forcing. While some studies suggest that Southern Ocean surface buoyancy forcing controls the AMOC depth, others have suggested alternatively that North Atlantic surface forcing or interior diabatic mixing plays the dominant role. To investigate the key processes that set the AMOC depth, here we carry out a number of MITgcm ocean-only simulations with surface forcing fields specified from the simulation results of three coupled climate models that span much of the range of glacial AMOC depth changes in phase 3 of the Paleoclimate Model Intercomparison Project (PMIP3). We find that the MITgcm simulations successfully reproduce the changes in AMOC depth between glacial and modern conditions simulated in these three PMIP3 models. By varying the restoring time scale in the surface forcing, we show that the AMOC depth is more strongly constrained by the surface density field than the surface buoyancy flux field. Based on these results, we propose a mechanism by which the surface density fields in the high latitudes of both hemispheres are connected to the AMOC depth. We illustrate the mechanism using MITgcm simulations with idealized surface forcing perturbations as well as an idealized conceptual geometric model. These results suggest that the AMOC depth is largely determined by the surface density fields in both the North Atlantic and the Southern Ocean.
Download PDF
Wang Y, Stewart AL. Scalings for eddy buoyancy transfer across continental slopes under retrograde winds. Ocean Modelling [Internet]. 2020;147 :101579. Publisher's VersionAbstract

Baroclinic eddy restratification strongly influences the ocean’s general circulation and tracer budgets, and has been routinely parameterized via the Gent–McWilliams (GM) scheme in coarse-resolution ocean climate models. These parameterizations have been improved via refinements of the GM eddy transfer coefficient using eddy-resolving simulations and theoretical developments. However, previous efforts have focused primarily on the open ocean, and the applicability of existing GM parameterization approaches to continental slopes remains to be addressed. In this study, we use a suite of eddy-resolving, process-oriented simulations to test scaling relationships between eddy buoyancy diffusivity, mean flow properties, and topographic geometries in simulations of baroclinic turbulence over continental slopes. We focus on the case of retrograde (i.e., opposing the direction of topographic wave propagation) winds, a configuration that arises commonly around the margins of the subtropical gyres.

Three types of scalings are examined, namely, the GEOMETRIC framework developed by Marshall et al. (2012), a new ”Cross-Front” (CF) scaling derived via dimensional arguments, and the mixing length theory (MLT)-based scalings tested recently by Jansen et al. (2015) over a flat ocean bed. The present study emphasizes the crucial role of the local slope parameter, defined as the ratio between the topographic slope and the depth-averaged isopycnal slope, in controlling the nonlinear eddy buoyancy fluxes. Both the GEOMETRIC framework and the CF scaling can reproduce the depth-averaged eddy buoyancy transfer across alongshore-uniform continental slopes, for suitably chosen constant prefactors. Generalization of these scalings across both continental slope and open ocean environments requires the introduction of prefactors that depend on the local slope parameter via empirically derived analytical functions. In contrast, the MLT-based scalings fail to quantify the eddy buoyancy transfer across alongshore-uniform continental slopes when constant prefactors are adopted, but can reproduce the cross-slope eddy flux when the prefactors are adapted via empirical functions of the local slope parameter. Application of these scalings in prognostic ocean simulations also depends on an accurate representation of standing eddies associated with the topographic corrugations of the continental slope. These findings offer a basis for extending existing approaches to parameterizing transient eddies, and call for future efforts to parameterize standing eddies in coarse-resolution ocean climate models.

Download PDF
2019
Stewart AL. Approximating isoneutral ocean transport via the Temporal Residual Mean. Fluids [Internet]. 2019;4 :179. Publisher's VersionAbstract
Ocean volume and tracer transports are commonly computed on density surfaces because doing so approximates the semi-Lagrangian mean advective transport. The resulting density-averaged transport can be related approximately to Eulerian-averaged quantities via the Temporal Residual Mean (TRM), valid in the limit of small isopycnal height fluctuations. This article builds on a formulation of the TRM for volume fluxes within Neutral Density surfaces, (the “NDTRM”), selected because Neutral Density surfaces are constructed to be as neutral as possible while still forming well-defined surfaces. This article derives a TRM, referred to as the “Neutral TRM” (NTRM), that approximates volume fluxes within surfaces whose vertical fluctuations are defined directly by the neutral relation. The purpose of the NTRM is to more closely approximate the semi-Lagrangian mean transport than the NDTRM, because the latter introduces errors associated with differences between the instantaneous state of the modeled/observed ocean and the reference climatology used to assign the Neutral Density variable. It is shown that the NDTRM collapses to the NTRM in the limiting case of a Neutral Density variable defined with reference to the Eulerian-mean salinity, potential temperature and pressure, rather than an external reference climatology, and therefore that the NTRM approximately advects this density variable. This prediction is verified directly using output from an idealized eddy-resolving numerical model. The NTRM therefore offers an efficient and accurate estimate of modeled semi-Lagrangian mean transports without reference to an external reference climatology, but requires that a Neutral Density variable be computed once from the model’s time-mean state in order to estimate isopycnal and diapycnal components of the transport.
Download PDF
Park H-S, Kim S-J, Stewart AL, Son S-W, Seo K-H. Mid-Holocene Northern Hemisphere warming driven by Arctic amplification and sea ice loss. Science Advances [Internet]. 2019;5 :eaax8203. Publisher's VersionAbstract
The Holocene thermal maximum was characterized by strong summer solar heating that substantially increased the summertime temperature relative to preindustrial climate. However, the summer warming was compensated by weaker winter insolation, and the annual mean temperature of the Holocene thermal maximum remains ambiguous. Using multimodel mid-Holocene simulations, we show that the annual mean Northern Hemisphere temperature is strongly correlated with the degree of Arctic amplification and sea ice loss. Additional model experiments show that the summer Arctic sea ice loss persists into winter and increases the mid- and high-latitude temperatures. These results are evaluated against four proxy datasets to verify that the annual mean northern high-latitude temperature during the mid-Holocene was warmer than the preindustrial climate, because of the seasonally rectified temperature increase driven by the Arctic amplification. This study offers a resolution to the “Holocene temperature conundrum”, a well-known discrepancy between paleo-proxies and climate model simulations of Holocene thermal maximum.
Download PDF
Stewart AL, Klocker A, Menemenlis D. Acceleration and overturning of the Antarctic Slope Current by winds, eddies, and tides. Journal of Physical Oceanography [Internet]. 2019;49 :2043-2074. Publisher's VersionAbstract
All exchanges between the open ocean and the Antarctic continental shelf must cross the Antarctic Slope Current (ASC). Previous studies indicate that these exchanges are strongly influenced by mesoscale and tidal variability, yet the mechanisms responsible for setting the ASC’s transport and structure have received relatively little attention. In this study the roles of winds, eddies, and tides in accelerating the ASC are investigated using a global ocean–sea ice simulation with very high resolution (1/48° grid spacing). It is found that the circulation along the continental slope is accelerated both by surface stresses, ultimately sourced from the easterly winds, and by mesoscale eddy vorticity fluxes. At the continental shelf break, the ASC exhibits a narrow (~30–50 km), swift (>0.2 m s−1) jet, consistent with in situ observations. In this jet the surface stress is substantially reduced, and may even vanish or be directed eastward, because the ocean surface speed matches or exceeds that of the sea ice. The shelfbreak jet is shown to be accelerated by tidal momentum advection, consistent with the phenomenon of tidal rectification. Consequently, the shoreward Ekman transport vanishes and thus the mean overturning circulation that steepens the Antarctic Slope Front (ASF) is primarily due to tidal acceleration. These findings imply that the circulation and mean overturning of the ASC are not only determined by near-Antarctic winds, but also depend crucially on sea ice cover, regionally-dependent mesoscale eddy activity over the continental slope, and the amplitude of tidal flows across the continental shelf break.
Download PDF
Hazel JE, Stewart AL. Are the near-Antarctic easterly winds weakening in response to enhancement of the Southern Annular Mode?. Journal of Climate [Internet]. 2019;32 :1895-1918. Publisher's VersionAbstract
Previous studies have highlighted the sensitivity of the Southern Ocean circulation to the strengthening, poleward-shifting westerlies, associated with the increasingly positive southern annular mode (SAM). The impacts of the SAM have been hypothesized to weaken momentum input to the ocean from the easterly winds around the Antarctic margins. Using ERA-Interim data, the authors show that the circumpolar-averaged easterly wind stress has not weakened over the past 3–4 decades, and, if anything, has slightly strengthened by around 7%. However, there has been a substantial increase in the seasonality of the easterlies, with a weakening of the easterly winds during austral summer and a strengthening during winter. A similar trend in the seasonality of the easterlies is found in three other reanalysis products that compare favorably with Antarctic meteorological observations. The authors associate the strengthening of the easterly winds during winter with an increase in the pressure gradient between the coast and the pole. Although the trend in the overall easterly wind strength is small, the change in the seasonal cycle may be expected to reduce the shoreward Ekman transport of summer surface waters and also to admit more warm Circumpolar Deep Water to the continental shelf in summer. Changes in the seasonal cycle of the near-coastal winds may also project onto seasonal formation and export of sea ice, fluctuations in the strengths of the Weddell and Ross Gyres, and seasonal export of Antarctic Bottom Water from the continental shelf.
Download PDF
Zhao KX, Stewart AL, McWilliams JC. Sill-Influenced Exchange Flows in Ice Shelf Cavities. Journal of Physical Oceanography [Internet]. 2019;49 (1) :163-191. Publisher's VersionAbstract
Bathymetric sills are important features in the ocean-filled cavities beneath a few fast-retreating ice shelves in West Antarctica and northern Greenland. The sills can be high enough to obstruct the cavity circulation and thereby modulate glacial melt rates. This study focuses on the idealized problem of diabatically driven, sill-constrained overturning circulation in a cavity. The circulation beneath fast-melting ice shelves can generally be characterized by an inflow of relatively warm dense water (with temperatures of a few degrees Celsius above the local freezing point) at depth and cold, less-dense, outflowing water, which exhibits an approximately two-layer structure in observations. We use a two-layer isopycnal hydrostatic model to study the cross-sill exchange of these waters in ice shelf cavities wide enough to be rotationally dominated. A quasigeostrophic constraint is determined for the transport imposed by the stratification. Relative to this constraint, the key parameters controlling the transport and its variability are the sill height relative to the bottom layer thickness and the strength of the friction relative to the potential vorticity (PV) gradient imposed by the sill. By varying these two key parameters, we simulate a diversity of flow phenomena. For a given meridional pressure gradient, the cross-sill transport is controlled by sill height beyond a critical threshold in the eddy-permitting, low-friction regime, while it is insensitive to friction in both the low-friction and high-friction regimes. We present theoretical ideas to explain the flow characteristics: a Stommel boundary layer for the friction-dominated regime; mean–eddy PV balances and energy conversion in the low-friction, low-sill regime; and hydraulic control in the low-friction, high-sill regime, with various estimates for transport in each of these regimes.
Download PDF
2018
Thompson AF, Stewart AL, Spence P, Heywood KJ. The Antarctic Slope Front in a Changing Climate. Reviews of Geophysics [Internet]. 2018;56 :741-770. Publisher's VersionAbstract

The Antarctic Slope Current (ASC) is a coherent circulation feature that rings the Antarctic continental shelf and regulates the flow of water toward the Antarctic coastline. The structure and variability of the ASC influences key processes near the Antarctic coastline that have global implications, such as the melting of Antarctic ice shelves and water mass formation that determines the strength of the global overturning circulation. Recent theoretical, modeling, and observational advances have revealed new dynamical properties of the ASC, making it timely to review. Earlier reviews of the ASC focused largely on local classifications of water properties of the ASC’s primary front. Here we instead provide a classification of the current’s frontal structure based on the dynamical mechanisms that govern both the along-slope and cross-slope circulation; these two modes of circulation are strongly coupled, similar to the Antarctic Circumpolar Current. Highly variable motions, such as dense overflows, tides, and eddies are shown to be critical components of cross-slope and cross-shelf exchange, but understanding of how the distribution and intensity of these processes will evolve in a changing climate remains poor due to observational and modeling limitations. Results linking the ASC to larger modes of climate variability, such as El Niño, show that the ASC is an integral part of global climate. An improved dynamical understanding of the ASC is still needed to accurately model and predict future Antarctic sea ice extent, the stability of the Antarctic ice sheets, and the Southern Ocean’s contribution to the global carbon cycle.

Download PDF

Pages