Antarctic ice shelves are losing mass at drastically different rates, primarily due to differing rates of oceanic heat supply to their bases. However, a generalized theory for the inflow of relatively warm water into ice shelf cavities is lacking. This study proposes such a theory based on a geostrophically constrained inflow, combined with a threshold bathymetric elevation, the Highest Unconnected isoBath (HUB), that obstructs warm water access to ice shelf grounding lines. This theory captures ∼ 90% of the variance in melt rates across a suite of idealized process‐oriented ocean/ice shelf simulations with quasi‐randomized geometries. Applied to observations of ice shelf geometries and offshore hydrography, the theory captures ∼80% of the variance in measured ice shelf melt rates. These findings provide a generalized theoretical framework for melt resulting from buoyancy‐driven warm water access to geometrically complex Antarctic ice shelf cavities.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Abstract -
Abstract Eastern boundary upwelling systems (EBUSs) host equatorward wind-driven near-surface currents overlying poleward subsurface undercurrents. Various previous theories for these undercurrents have emphasized the role of poleward alongshore pressure gradient forces (APFs). Energetic mesoscale variability may also serve to accelerate undercurrents via mesoscale stirring of the potential vorticity gradient imposed by the continental slope. However, it remains unclear whether this eddy rectification mechanism contributes substantially to driving poleward undercurrents in EBUS. This study isolates the influence of eddy rectification on undercurrents via a suite of idealized simulations forced either by alongshore winds, with or without an APF, or by randomly generated mesoscale eddies. It is found that the simulations develop undercurrents with strengths comparable to those found in nature in both wind-forced and randomly forced experiments. Analysis of the momentum budget reveals that the along-isobath undercurrent flow is accelerated by isopycnal advective eddy momentum fluxes and the APF and retarded by frictional drag. The undercurrent acceleration may manifest as eddy momentum fluxes or as topographic form stress depending on the coordinate system used to compute the momentum budget, which reconciles these findings with previous work that linked eddy acceleration of the undercurrent to topographic form stress. The leading-order momentum balance motivates a scaling for the strength of the undercurrent that explains most of the variance across the simulations. These findings indicate that eddy rectification is of comparable importance to the APF in driving poleward undercurrents in EBUSs and motivate further work to diagnose this effect in high-resolution models and observations and to parameterize it in coarse-resolution ocean/climate models.
-
Abstract Depth‐averaged eddy buoyancy diffusivities across continental shelves and slopes are investigated using a suite of eddy‐resolving, process‐oriented simulations of prograde frontal currents characterized by isopycnals tilted in the opposite direction to the seafloor, a flow regime commonly found along continental margins under downwelling‐favorable winds or occupied by buoyant boundary currents. The diagnosed cross‐slope eddy diffusivity varies by up to three orders of magnitude, decaying from
in the relatively flat‐bottomed region to over the steep continental slope, consistent with previously reported suppression effects of steep topography on baroclinic eddy fluxes. To theoretically constrain the simulated cross‐slope eddy fluxes, we examine extant scalings for eddy buoyancy diffusivities across prograde shelf/slope fronts and in flat‐bottomed oceans. Among all tested scalings, the GEOMETRIC framework developed by D. P. Marshall et al. (2012, https://doi.org/10.1175/JPO-D-11-048.1 ) and a parametrically similar Eady scale‐based scaling proposed by Jansen et al. (2015,https://doi.org/10.1016/j.ocemod.2015.05.007 ) most accurately reproduce the diagnosed eddy diffusivities across the entire shelf‐to‐open‐ocean analysis regions in our simulations. This result relies upon the incorporation of the topographic suppression effects on eddy fluxes, quantified via analytical functions of the slope Burger number, into the scaling prefactor coefficients. The predictive skills of the GEOMETRIC and Eady scale‐based scalings are shown to be insensitive to the presence of along‐slope topographic corrugations. This work lays a foundation for parameterizing eddy buoyancy fluxes across large‐scale prograde shelf/slope fronts in coarse‐resolution ocean models. -
Abstract A current along a sloping bottom gives rise to upwelling, or downwelling Ekman transport within the stratified bottom boundary layer (BBL), also known as the bottom Ekman layer. In 1D models of slope currents, geostrophic vertical shear resulting from horizontal buoyancy gradients within the BBL is predicted to eventually bring the bottom stress to zero, leading to a shutdown, or “arrest,” of the BBL. Using 3D ROMS simulations, we explore how the dynamics of buoyancy adjustment in a current‐ridge encounter problem differs from 1D and 2D temporal spin up problems. We show that in a downwelling BBL, the destruction of the ageostrophic BBL shear, and hence the bottom stress, is accomplished primarily by nonlinear straining effects during the initial topographic encounter. As the current advects along the ridge slopes, the BBL deepens and evolves toward thermal wind balance. The onset of negative potential vorticitymodes of instability and their subsequent dissipation partially offsets the reduction of the BBL dissipation during the ridge‐current interaction. On the upwelling side, although the bottom stress weakens substantially during the encounter, the BBL experiences a horizontal inflectional point instability and separates from the slopes before sustained along‐slope stress reduction can occur. In all our solutions, both the upwelling and downwelling BBLs are in a partially arrested state when the current separates from the ridge slope, characterized by a reduced, but non‐zero bottom stress on the slopes.
-
Abstract Fjord circulation modulates the connection between marine‐terminating glaciers and the ocean currents offshore. These fjords exhibit both overturning and horizontal recirculations, which are driven by water mass transformation at the head of the fjord via subglacial discharge plumes and distributed meltwater plumes. However, little is known about how various fjord characteristics influence the interaction between 3D fjord circulation and glacial melt. In this study, high‐resolution numerical simulations of idealized glacial fjords demonstrate that recirculation strength controls melt, which feeds back on overturning and recirculation. The relationships between overturning, recirculation, and melt rate are well predicted by vorticity balance, reduced‐order melt parameterizations, and empirical scaling arguments. These theories allow us to take into account the near‐glacier horizontal velocities, which yield improved predictions of fjord overturning, recirculation, and glacial melt.
-
Abstract The ongoing Arctic warming has been pronounced in winter and has been associated with an increase in downward longwave radiation. While previous studies have demonstrated that poleward moisture flux into the Arctic strengthens downward longwave radiation, less attention has been given to the impact of the accompanying increase in snowfall. Here, utilizing state-of-the-art sea ice models, we show that typical winter snowfall (snow water equivalent) anomalies of around 1.0 cm, accompanied by positive downward longwave radiation anomalies of ∼5 W m−2, can cause basinwide sea ice thinning by around 5 cm in the following spring over the Arctic seas in the Eurasian–Pacific seas. In extreme cases, this is followed by a shrinking of summer ice extent. In the winter of 2016/17, anomalously strong warm, moist air transport combined with ∼2.5-cm increase in snowfall (snow water equivalent) decreased spring ice thickness by ∼10 cm and decreased the following summer sea ice extent by 5%–30%. This study suggests that small changes in the pattern and volume of winter snowfall can strongly impact the sea ice thickness and extent in the following seasons.
-
Abstract 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.
-
Abstract 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. -
Free, publicly-accessible full text available December 1, 2025
-
Abstract 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.