Publications

2021
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
Park H-S, Kim S-J, Seo K-H, Stewart AL, Kim S-Y. The impact of Arctic sea ice loss on mid-Holocene climate. Nature Communications [Internet]. 2018;9 (1) :4571. Publisher's VersionAbstract
Mid-Holocene climate was characterized by strong summer solar heating that decreased Arctic sea ice cover. Motivated by recent studies identifying Arctic sea ice loss as a key driver of future climate change, we separate the influences of Arctic sea ice loss on mid-Holocene climate. By performing idealized climate model perturbation experiments, we show that Arctic sea ice loss causes zonally asymmetric surface temperature responses especially in winter: sea ice loss warms North America and the North Pacific, which would otherwise be much colder due to weaker winter insolation. In contrast, over East Asia, sea ice loss slightly decreases the temperature in early winter. These temperature responses are associated with the weakening of mid-high latitude westerlies and polar stratospheric warming. Sea ice loss also weakens the Atlantic meridional overturning circulation, although this weakening signal diminishes after 150–200 years of model integration. These results suggest that mid-Holocene climate changes should be interpreted in terms of both Arctic sea ice cover and insolation forcing.
Download PDF
Stewart AL, Klocker A, Menemenlis D. Circum-Antarctic shoreward heat transport derived from an eddy- and tide-resolving simulation. Geophysical Research Letters [Internet]. 2018;45 :834-845. Publisher's VersionAbstract

Almost all heat reaching the bases of Antarctica's ice shelves originates from warm Circumpolar Deep Water in the open Southern Ocean. This study quantifies the roles of mean and transient flows in transporting heat across almost the entire Antarctic continental slope and shelf using an ocean/sea ice model run at eddy- and tide-resolving (1/48°) horizontal resolution. Heat transfer by transient flows is approximately attributed to eddies and tides via a decomposition into time scales shorter than and longer than 1 day, respectively. It is shown that eddies transfer heat across the continental slope (ocean depths greater than 1,500 m), but tides produce a stronger shoreward heat flux across the shelf break (ocean depths between 500 m and 1,000 m). However, the tidal heat fluxes are approximately compensated by mean flows, leaving the eddy heat flux to balance the net shoreward heat transport. The eddy-driven cross-slope overturning circulation is too weak to account for the eddy heat flux. This suggests that isopycnal eddy stirring is the principal mechanism of shoreward heat transport around Antarctica, though likely modulated by tides and surface forcing.

Download PDF
Wang Y, Stewart AL. Eddy dynamics over continental slopes under retrograde winds: Insights from a model inter-comparison. Ocean Modelling [Internet]. 2018;121 :1-18. Publisher's VersionAbstract

Mesoscale eddies are ubiquitous in the ocean and play a key role in exchanges across continental slopes. In this study the properties of wind-driven baroclinic turbulence are investigated using eddy-resolving process simulations, focusing on the case of retrograde winds that arises around the margins of the subtropical gyres. In contrast to a flat-bottomed ocean, over steep slopes eddies develop from baroclinic instabilities are confined to the top few hundred meters. Deeper in the water column baroclinic instability and vertical momentum transfer are suppressed, so wind-input momentum is exported toward the open ocean by eddies before traversing down to the ocean bed. Close to the sloping topography, eddy energy sourced from the upper ocean is converted to potential energy, steepening isopycnals and driving bottom-trapped prograde flows. This process is associated with upgradient lateral buoyancy fluxes and downgradient isopycnal potential vorticity fluxes, and cannot be reproduced via linear stability calculations.

These properties of wind-driven shelf/slope turbulence are contrasted against simulations with flat bathymetry. The key differences described above hinge on the flow close to the steep topographic slope, which may be sensitive to the model’s vertical coordinate system. The simulations are therefore replicated using models that employ geopotential coordinates, terrain-following coordinates, and isopycnal coordinates. Quantitative inter-model discrepancies in the momentum and energy budgets are much more pronounced in the presence of a steep bottom slope. However, the key findings of this study are consistent across the models, suggesting that they are robust and warrant incorporation into parameterizations of eddy transfer across continental slopes.

Download PDF
Sun S, Eisenman I, Stewart AL. Does Southern Ocean surface forcing shape the global ocean overturning circulation?. Geophysical Research Letters [Internet]. 2018;45 (5) :2413-2423. Publisher's VersionAbstract

Paleoclimate proxy data suggest that the Atlantic Meridional Overturning Circulation (AMOC) was shallower at the Last Glacial Maximum (LGM) than its preindustrial (PI) depth. Previous studies have suggested that this shoaling necessarily accompanies Antarctic sea ice expansion at the LGM. Here the influence of Southern Ocean surface forcing on the AMOC depth is investigated using ocean‐only simulations from a state‐of‐the‐art climate model with surface forcing specified from the output of previous coupled PI and LGM simulations. In contrast to previous expectations, we find that applying LGM surface forcing in the Southern Ocean and PI surface forcing elsewhere causes the AMOC to shoal only about half as much as when LGM surface forcing is applied globally. We show that this occurs because diapycnal mixing renders the Southern Ocean overturning circulation more diabatic than previously assumed, which diminishes the influence of Southern Ocean surface buoyancy forcing on the depth of the AMOC.

Download PDF
Park H-S, Stewart AL, Son J-H. Dynamic and thermodynamic effects of the winter Arctic Oscillation on summer sea ice extent. Journal of Climate [Internet]. 2018;31 :1483-1497. Publisher's VersionAbstract

Arctic summer sea ice extent exhibits substantial interannual variability, as is highlighted by the remarkable recovery in sea ice extent in 2013 following the record minimum in the summer of 2012. Here, the mechanism via which Arctic Oscillation (AO)-induced ice thickness changes impact summer sea ice is explored, using observations and reanalysis data. A positive AO weakens the basin-scale anticyclonic sea ice drift and decreases the winter ice thickness by 15 and 10 cm in the Eurasian and the Pacific sectors of the Arctic, respectively. Three reanalysis datasets show that the upward surface heat fluxes are reduced over wide areas of the Arctic, suppressing the ice growth during the positive AO winters. The winter dynamic and thermodynamic thinning preconditions the ice for enhanced radiative forcing via the ice–albedo feedback in late spring–summer, leading to an additional 10 cm of thinning over the Pacific sector of the Arctic. Because of these winter AO-induced dynamic and thermodynamics effects, the winter AO explains about 22% (r = −0.48) of the interannual variance of September sea ice extent from 1980 to 2015.

Download PDF

Pages