Computer models are necessary for understanding and predicting marine ice sheet behaviour. However, there is uncertainty over implementation of physical processes at the ice base, both for grounded ...and floating glacial ice. Here we implement several sliding relations in a marine ice sheet flow-line model accounting for all stress components and demonstrate that model resolution requirements are strongly dependent on both the choice of basal sliding relation and the spatial distribution of ice shelf basal melting.Sliding relations that reduce the magnitude of the step change in basal drag from grounded ice to floating ice (where basal drag is set to zero) show reduced dependence on resolution compared to a commonly used relation, in which basal drag is purely a power law function of basal ice velocity. Sliding relations in which basal drag goes smoothly to zero as the grounding line is approached from inland (due to a physically motivated incorporation of effective pressure at the bed) provide further reduction in resolution dependence.A similar issue is found with the imposition of basal melt under the floating part of the ice shelf: melt parameterisations that reduce the abruptness of change in basal melting from grounded ice (where basal melt is set to zero) to floating ice provide improved convergence with resolution compared to parameterisations in which high melt occurs adjacent to the grounding line.Thus physical processes, such as sub-glacial outflow (which could cause high melt near the grounding line), impact on capability to simulate marine ice sheets. If there exists an abrupt change across the grounding line in either basal drag or basal melting, then high resolution will be required to solve the problem. However, the plausible combination of a physical dependency of basal drag on effective pressure, and the possibility of low ice shelf basal melt rates next to the grounding line, may mean that some marine ice sheet systems can be reliably simulated at a coarser resolution than currently thought necessary.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
The Antarctic Ice Sheet loses about half its mass through ocean-driven melting of its fringing ice shelves. However, the ocean processes governing ice shelf melting are not well understood, ...contributing to uncertainty in projections of Antarctica's contribution to global sea level. We use high-resolution large-eddy simulation to examine ocean-driven melt, in a geophysical-scale model of the turbulent ice shelf-ocean boundary layer, focusing on the ocean conditions observed beneath the Ross Ice Shelf. We quantify the role of double-diffusive convection in determining ice shelf melt rates and oceanic mixed layer properties in relatively warm and low-velocity cavity environments. We demonstrate that double-diffusive convection is the first-order process controlling the melt rate and mixed layer evolution at these flow conditions, even more important than vertical shear due to a mean flow, and is responsible for the step-like temperature and salinity structure, or thermohaline staircase, observed beneath the ice. A robust feature of the multiday simulations is a growing saline diffusive sublayer that drives a time-dependent melt rate. This melt rate is lower than current ice-ocean parameterizations, which consider only shear-controlled turbulent melting, would predict. Our main finding is that double-diffusive convection is an important process beneath ice shelves, yet is currently neglected in ocean-climate models.
Full text
Available for:
BFBNIB, NMLJ, NUK, PNG, SAZU, UL, UM, UPUK
Previous studies of Totten Ice Shelf have employed surface velocity
measurements to estimate its mass balance and understand its sensitivities to
interannual changes in climate forcing. However, ...displacement measurements
acquired over timescales of days to weeks may not accurately characterize
long-term flow rates wherein ice velocity fluctuates with the seasons.
Quantifying annual mass budgets or analyzing interannual changes in ice
velocity requires knowing when and where observations of glacier velocity
could be aliased by subannual variability. Here, we analyze 16 years of
velocity data for Totten Ice Shelf, which we generate at subannual resolution
by applying feature-tracking algorithms to several hundred satellite image
pairs. We identify a seasonal cycle characterized by a spring to autumn
speedup of more than 100 m yr−1 close to the ice front. The amplitude
of the seasonal cycle diminishes with distance from the open ocean,
suggesting the presence of a resistive back stress at the ice front that is
strongest in winter. Springtime acceleration precedes summer surface melt and
is not attributable to thinning from basal melt. We attribute the onset of
ice shelf acceleration each spring to the loss of buttressing from the
breakup of seasonal landfast sea ice.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Water mass transformation (WMT) around the Antarctic margin controls Antarctica Bottom Water formation and the abyssal limb of the global meridional overturning circulation, besides mediating ...ocean-ice shelf exchange, ice sheet stability and its contribution to sea level rise. However, the mechanisms controlling the rate of WMT in the Antarctic shelf are poorly understood due to the lack of observations and the inability of climate models to simulate those mechanisms, in particular beneath the floating ice shelves. We used a circum-Antarctic ocean-ice shelf model to assess the contribution of surface fluxes, mixing, and ocean-ice shelf interaction to the WMT on the continental shelf. The salt budget dominates the WMT rates, with only a secondary contribution from the heat budget. Basal melt of ice shelves drives buoyancy gain at lighter density classes (27.2<σ
θ
< 27.6 kg m
-3
), while salt input associated with sea-ice growth in coastal polynyas drives buoyancy loss at heavier densities (σ
θ
> 27.6). We found a large sensitivity of the WMT rates to model horizontal resolution, tides and topography within the Filchner-Ronne, East and West Antarctica ice shelf cavities. In the Filchner-Ronne Ice Shelf, an anticyclonic circulation in front of the Ronne Depression regulates the rates of basal melting/refreezing and WMT and is substantially affected by tides and model resolution. Model resolution is also found to affect the Antarctic Slope Current in both East and West Antarctica, impacting the on-shelf heat delivery, basal melt and WMT. Moreover, the representation of the ice shelf draft associated with model resolution impacts the freezing temperature and thus basal melt and WMT rates in the East Antarctica. These results highlight the importance of resolving small-scale features of the flow and topography, and to include the effects of tidal forcing, to adequately represent water mass transformations on the shelf that directly influence the abyssal global overturning circulation.
The Totten Ice Shelf (IS) has a large drainage basin,
much of which is grounded below sea level, leaving the glacier vulnerable to
retreat through the marine ice sheet instability mechanism. The ice ...shelf
has also been shown to be sensitive to changes in calving rate, as a very
small retreat of the calving front from its current position is predicted to
cause a change in flow at the grounding line. Therefore understanding the
processes behind calving on the Totten IS is key to predicting its future
sea level rise contribution. Here we use the Helsinki Discrete Element Model (HiDEM)
to show that not all of the fractures visible at the front of the
Totten IS are produced locally, but that the across-flow basal crevasses,
which are part of the distinctive cross-cutting fracture pattern, are
advected into the calving front area from upstream. A separate simulation of
the grounding line shows that re-grounding points may be key areas of basal
crevasse production, and can produce basal crevasses in both an along- and
across-flow orientation. The along-flow basal crevasses at the grounding
line may be a possible precursor to basal channels, while we suggest the
across-flow grounding-line fractures are the source of the across-flow
features observed at the calving front. We use two additional models to
simulate the evolution of basal fractures as they advect downstream,
demonstrating that both strain and ocean melt have the potential to deform
narrow fractures into the broad basal features observed near the calving
front. The wide range of factors which influence fracture patterns and
calving on this glacier will be a challenge for predicting its future mass loss.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Changes in ocean-driven basal melting have a key influence on the stability of ice shelves, the mass loss from the ice sheet, ocean circulation, and global sea level rise. Coupled ice sheet–ocean ...models play a critical role in understanding future ice sheet evolution and examining the processes governing ice sheet responses to basal melting. However, as a new approach, coupled ice sheet–ocean systems come with new challenges, and the impacts of solutions implemented to date have not been investigated. An emergent feature in several contributing coupled models to the 1st Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP1) was a time-varying oscillation in basal melt rates. Here, we use a recently developed coupling framework, FISOC (v1.1), to connect the modified ocean model ROMSIceShelf (v1.0) and ice sheet model Elmer/Ice (v9.0), to investigate the origin and implications of the feature and, more generally, the impact of coupled modeling strategies on the simulated basal melt in an idealized ice shelf cavity based on the MISOMIP setup. We found the spatial-averaged basal melt rates (3.56 m yr−1) oscillated with an amplitude ∼0.7 m yr−1 and approximate period of ∼6 years between year 30 and 100 depending on the experimental design. The melt oscillations emerged in the coupled system and the standalone ocean model using a prescribed change of cavity geometry. We found that the oscillation feature is closely related to the discretized ungrounding of the ice sheet, exposing new ocean, and is likely strengthened by a combination of positive buoyancy–melt feedback and/or melt–geometry feedback near the grounding line, and the frequent coupling of ice geometry and ocean evolution. Sensitivity tests demonstrate that the oscillation feature is always present, regardless of the choice of coupling interval, vertical resolution in the ocean model, tracer properties of cells ungrounded by the retreating ice sheet, or the dependency of friction velocities to the vertical resolution. However, the amplitude, phase, and sub-cycle variability of the oscillation varied significantly across the different configurations. We were unable to ultimately determine whether the feature arises purely due to numerical issues (related to discretization) or a compounding of multiple physical processes amplifying a numerical artifact. We suggest a pathway and choices of physical parameters to help other efforts understand the coupled ice sheet–ocean system using numerical models.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Changes in ocean-driven basal melting have a key influence on the stability of ice shelves, the mass loss from the ice sheet, ocean circulation, and global sea level rise. Coupled ice sheet-ocean ...models play a critical role in understanding future ice sheet evolution and examining the processes governing ice sheet responses to basal melting. However, as a new approach, coupled ice sheet-ocean systems come with new challenges, and the impacts of solutions implemented to date have not been investigated. An emergent feature in several contributing coupled models to the 1st Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP1) was a time-varying oscillation in basal melt rates. Here, we use a recently developed coupling framework, FISOC (v1.1), to connect the modified ocean model ROMSIceShelf (v1.0) and ice sheet model Elmer/Ice (v9.0), to investigate the origin and implications of the feature and, more generally, the impact of coupled modeling strategies on the simulated basal melt in an idealized ice shelf cavity based on the MISOMIP setup. We found the spatial-averaged basal melt rates (3.56 m yr.sup.-1) oscillated with an amplitude â¼0.7 m yr.sup.-1 and approximate period of â¼6 years between year 30 and 100 depending on the experimental design. The melt oscillations emerged in the coupled system and the standalone ocean model using a prescribed change of cavity geometry. We found that the oscillation feature is closely related to the discretized ungrounding of the ice sheet, exposing new ocean, and is likely strengthened by a combination of positive buoyancy-melt feedback and/or melt-geometry feedback near the grounding line, and the frequent coupling of ice geometry and ocean evolution. Sensitivity tests demonstrate that the oscillation feature is always present, regardless of the choice of coupling interval, vertical resolution in the ocean model, tracer properties of cells ungrounded by the retreating ice sheet, or the dependency of friction velocities to the vertical resolution. However, the amplitude, phase, and sub-cycle variability of the oscillation varied significantly across the different configurations. We were unable to ultimately determine whether the feature arises purely due to numerical issues (related to discretization) or a compounding of multiple physical processes amplifying a numerical artifact. We suggest a pathway and choices of physical parameters to help other efforts understand the coupled ice sheet-ocean system using numerical models.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Abstract
The Antarctic Slope Front (ASF) is located along much of the Antarctic continental shelf break and helps to maintain a barrier to the movement of Circumpolar Deep Water (CDW) onto the ...continental shelf. The stability of the ASF has a major control on cross-shelf heat transport and ocean-driven basal melting of Antarctic ice shelves. Here, the ASF dynamics are investigated for continental shelves with weak dense shelf water (DSW) formation, which are thought to have a stable ASF, common for regions in East Antarctica. Using an ocean process model, this study demonstrates how offshore bottom Ekman transport of shelf waters leads to the development of a deep bottom mixed layer at the lower continental slope, and subsequently determines an intrinsic variability of the ASF. The ASF variability is characterized by instability events that affect the entire water column and occur every 5–10 years and last for approximately half a year. During these instability events, the cross-shelf density gradient weakens and CDW moves closer to the continent. Stronger winds increase the formation rate of the bottom mixed layer, which causes a subsequent increase of instability events. If the observed freshening trend of continental shelf waters leads to weaker DSW formation, more regions might be vulnerable for the ASF variability to develop in the future.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Totten Glacier is a fast‐moving East Antarctic outlet with the potential for significant future sea‐level contributions. We deployed four autonomous phase‐sensitive radars on its ice shelf to monitor ...ice‐ocean interactions near its grounding zone and made active source seismic observations to constrain gravity‐derived bathymetry models. We observe an asymmetry in basal melting with mean melt rates along the grounding zone differing by up to 20 m/a. Our new bathymetry model reveals that this melt rate asymmetry coincides with an asymmetry in water column thickness and that the low‐melting ice‐shelf portion is shielded from the main cavity circulation. A 2‐year record yields year‐to‐year melt rate variability of 7–9 m/a with no seasonal cycle. Our results highlight the key role of bathymetry near grounding lines for accurate modeling of ice‐shelf melt, and the importance of sustained multi‐year monitoring, especially at ice‐shelf cavities where the dominant melt rate drivers vary primarily inter‐annually.
Plain Language Summary
The point were the Antarctic Ice Sheet goes afloat on the ocean represents a critical region, where minor variations in melt rates can impact glacier flow and influence the rate of sea‐level rise. East Antarctica's Totten Glacier holds the potential to raise global sea level by several meters. Therefore, to understand the conditions it is exposed to, we measured melt rates for 2 years in several key locations near the point where the ice first touches the ocean. Our new measurements of the shape of the Totten Ice Shelf cavity help explain an observed spatial pattern of basal melting and together with local melt rate data resolve a disagreement between existing melt rate estimates from remote‐sensing methods.
Key Points
Totten Glacier melt rates vary spatially between 0 and over 20 m/a; differences are explained by water column thickness variations from updated bathymetry
Temporal melt rate variability is primarily inter‐annual; melt rates differ by 7–9 m/a over two observed years and there is no clear seasonal cycle
Contrary to previous findings, we find no topographic barriers to the intrusion of warm water to the Totten Glacier grounding zone
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
A phase‐sensitive radar (ApRES) was deployed on Totten Ice Shelf to provide the first in situ basal melt estimate at this dynamic East Antarctic ice shelf. Observations of internal ice dynamics at ...tidal time scales showed that early arrivals from off‐nadir reflectors obscure the true depth of the ice shelf base. Using the observed tidal deformation, the true base was found to lie at 1,910–1,950‐m depth, at 350–400 m greater range than the first reflection from an ice‐ocean interface. The robustness of the basal melt rate estimate was increased by using multiple basal reflections over the radar footprint, yielding a melt rate of 22 ± 2.1 m a−1. The ApRES estimate is over 40% lower than the three existing satellite estimates covering Totten Ice Shelf. This difference in basal melt is dynamically significant and highlights the need for independent melt rate estimates using complementary instrumentation and techniques that rely on different sets of assumptions.
Plain Language Summary
Observations of the rate of melting at the base of ice shelves are needed to model accurately ice sheet evolution. Local measurements are scarce, yet necessary for validation of satellite products and ocean models. We deployed a phase‐sensitive radar in the proximity of grounding zone of Totten Ice Shelf in East Antarctica, to measure basal melt in this dynamic region where uncertainties on melt rate estimates are high. We developed a method that accounts for basal geometry complexities and derived a melt rate estimate of ∼22 m per year, which is lower than previous estimates, but it confirms that the basal melt rate Totten Ice Shelf experiences is unusually high for East Antarctica.
Key Points
First in situ radar‐derived basal melt estimate on Totten Ice Shelf yields 22 ± 2.1 m a−1, at least 40% lower than existing satellite estimates
Radar‐derived observation of tidal ice dynamics constrains estimate of ice thickness for a complex base
The use of multiple basal reflections in melt derivation increases robustness of the estimate when early off‐nadir returns are present
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK