Ocean‐driven melting of ice shelves is a primary mechanism for ice loss from Antarctica. However, due to the difficulty in accessing the sub‐ice‐shelf ocean cavity, the relationship between ice‐shelf melting and ocean conditions is poorly understood, particularly near the grounding zone, where the ice transitions from grounded to floating. We present the first borehole oceanographic observations from the grounding zone of the Ross Ice Shelf, Antarctica's largest ice shelf by area. Contrary to predictions that tidal currents near grounding zones mix the water column, we found that Ross Ice Shelf waters were vertically stratified. Current velocities at mid‐depth in the ocean cavity did not change significantly over measurement periods at two different parts of the tidal cycle. The observed stratification resulted in low melt rates near this portion of the grounding zone, inferred from phase‐sensitive radar observations. These melt rates were generally <10 cm yr^‐1, which is lower than average for the Ross Ice Shelf (~20 cm yr^‐1). Melt rates may be higher at portions of the grounding zone that experience higher subglacial discharge or stronger tidal mixing. Stratification in the cavity at the borehole site was prone to diffusive convection as a result of ice‐shelf melting. Since diffusive convection influences vertical heat and salt fluxes differently than shear‐driven turbulence, this process may affect ice‐shelf melting and merits further consideration in ocean models of sub‐ice‐shelf circulation.

The processes responsible for freshwater flux from the Antarctic Ice Sheet (AIS), ice‐shelf basal melting and iceberg calving, are generally poorly represented in current Earth System Models (ESMs). Here we document the cryosphere configuration of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) v1.2. This includes simulating Antarctic ice‐shelf basal melting, which has been implemented through simulating the ocean circulation within static Antarctic ice‐shelf cavities, allowing for the ability to calculate ice‐shelf basal melt rates from the associated heat and freshwater fluxes. In addition, we added the capability to prescribe forcing from iceberg melt, allowing for realistic representation of the other dominant mass loss process from the AIS. In standard resolution simulations (using a noneddying ocean) under preindustrial climate forcing, we find high sensitivity of modeled ocean/ice shelf interactions to the ocean state, which can result in a transition to a high basal melt regime under the Filchner‐Ronne Ice Shelf (FRIS), presenting a significant challenge to representing the ocean/ice shelf system in a coupled ESM. We show that inclusion of a spatially dependent parameterization of eddy‐induced transport reduces biases in water mass properties on the Antarctic continental shelf. With these improvements, E3SM produces realistic ice‐shelf basal melt rates across the continent that are generally within the range inferred from observations. The accurate representation of ice‐shelf basal melting within a coupled ESM is an important step toward reducing uncertainties in projections of the Antarctic response to climate change and Antarctica's contribution to global sea‐level rise.

Knowledge gaps about how the ocean melts Antarctica’s ice shelves, borne from a lack of observations, lead to large uncertainties in sea level predictions. Using high-resolution maps of the underside of Dotson Ice Shelf, West Antarctica, we reveal the imprint that ice shelf basal melting leaves on the ice. Convection and intermittent warm water intrusions form widespread terraced features through slow melting in quiescent areas, while shear-driven turbulence rapidly melts smooth, eroded topographies in outflow areas, as well as enigmatic teardrop-shaped indentations that result from boundary-layer flow rotation. Full-thickness ice fractures, with bases modified by basal melting and convective processes, are observed throughout the area. This new wealth of processes, all active under a single ice shelf, must be considered to accurately predict future Antarctic ice shelf melt.

Mass loss from the Antarctic ice sheet is sensitive to conditions in ice shelf grounding zones, the transition between grounded and floating ice. To observe tidal dynamics in the grounding zone, we moored an ocean pressure sensor to Ross Ice Shelf, recording data for 54 days. In this region the ice shelf is brought out of hydrostatic equilibrium by the flexural rigidity of ice, yet we found that tidal pressure variations at a constant geopotential surface were similar within and outside of the grounding zone. This implies that the grounding zone ocean cavity was overpressurized at high tide and underpressurized at low tide by up to 10 kPa with respect to glaciostatic pressure at the ice shelf base. Phase lags between ocean pressure and vertical ice shelf motion were tens of minutes for diurnal and semidiurnal tides, an effect that has not been incorporated into ocean models of tidal currents below ice shelves. These tidal pressure variations may affect the production and export of meltwater in the subglacial environment and may increase basal crevasse heights in the grounding zone by several meters, according to linear elastic fracture mechanics. We find anomalously high tidal energy loss at the K₁ constituent in the grounding zone and hypothesize that this could be explained by seawater injection into the subglacial environment at high tide or internal tide generation through interactions with topography. These observations lay the foundation for improved representation of the grounding zone and its tidal dynamics in ocean circulation models of sub–ice shelf cavities.

Geothermal heat flux (GHF) is an important part of the basal heat budget of continental ice sheets. The difficulty of measuring GHF below ice sheets has directly hindered progress in the understanding of ice sheet dynamics. We present a new GHF measurement from below the West Antarctic Ice Sheet, made in subglacial sediment near the grounding zone of the Whillans Ice Stream. The measured GHF is 88 ± 7 mW m⁻², a relatively high value compared to other continental settings and to other GHF measurements along the eastern Ross Sea of 55 mW m⁻² and 69 ± 21 mW m⁻² but within the range of regional values indicated by geophysical estimates. The new GHF measurement was made ~100 km from the only other direct GHF measurement below the ice sheet, which was considerably higher at 285 ± 80 mW m⁻², suggesting spatial variability that could be explained by shallow magmatic intrusions or the advection of heat by crustal fluids. Analytical calculations suggest that spatial variability in GHF exceeds spatial variability in the conductive heat flux through ice along the Siple Coast. Accurate GHF measurements and high‐resolution GHF models may be necessary to reliably predict ice sheet evolution, including responses to ongoing and future climate change.

The U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 2.1 builds on E3SMv2 with several changes, with the most notable being the addition of the Fox-Kemper et al. (2011) mixed-layer eddy parameterization. This parameterization captures the effect of finite-amplitude, mixed-layer eddies as an overturning streamfunction and has the primary function of restratification. Herein, we outline the changes to the mean climate state of E3SM that were introduced by the addition of this parameterization. Overall, the presence of the submesoscale parameterization improves the fidelity of the v2.1 simulation by reducing the ocean surface biases in the North Atlantic present in v2, as illustrated by changes in the climatological sea surface temperature and salinity and the Arctic sea-ice extent. Other impacts include a slight shoaling of the mixed-layer depths in the North Atlantic and a small improvement in the Atlantic Meridional Overturning Circulation (AMOC). We note that the expected shoaling due to the parameterization is regionally dependent in our coupled configuration. In addition, we investigate why the parameterization and its impacts on mixed-layer depth have little impact on the simulated AMOC: despite increased dense-water formation in the Norwegian Sea, only a small fraction of the water formed makes its way south into the North Atlantic basin. Version 2.1 also exhibits small improvements in the atmospheric climatology, with smaller biases in many notable quantities and modes of variability.