Nioghalvfjerdsfjorden is a major outlet glacier in Northeast-Greenland. Although earlier studies showed that the floating part near the grounding line thinned by 30% between 1999 and 2014, the temporal ice loss evolution, its relation to external forcing and the implications for the grounded ice sheet remain largely unclear. By combining observations of surface features, ice thickness and bedrock data, we find that the ice shelf mass balance has been out of equilibrium since 2001, with large variations of the thinning rates on annual/multiannual time scales. Changes in ice flux and surface ablation are too small to produce this variability. An increased ocean heat flux is the most plausible cause of the observed thinning. For sustained environmental conditions, the ice shelf will lose large parts of its area within a few decades and ice modeling shows a significant, but locally restricted thinning upstream of the grounding line in response.

Abstract. Ice shelf break-up and disintegration events over the past 5 decades have led to speed-up, thinning, and retreat of upstream tributary glaciers and increases to rates of global sea-level rise. The southward progression of these episodes indicates a climatic cause and in turn suggests that the larger Larsen C and George VI ice shelves may undergo a similar collapse in the future. However, the extent to which removal of the Larsen C and George VI ice shelves will affect upstream tributary glaciers and add to global sea levels is unknown. Here we apply numerical ice-sheet models of varying complexity to show that the centennial sea-level commitment of Larsen C embayment glaciers following immediate shelf collapse is low (<2.5 mm to 2100, <4.2 mm to 2300). Despite its large size, Larsen C does not provide strong buttressing forces to upstream basins and its collapse does not result in large additional discharge from its tributary glaciers in any of our model scenarios. In contrast, the response of inland glaciers to a collapse of the George VI Ice Shelf may add up to 8 mm to global sea levels by 2100 and 22 mm by 2300 due in part to the mechanism of marine ice sheet instability. Our results demonstrate the varying and relative importance to sea level of the large Antarctic Peninsula ice shelves considered to present a risk of collapse.

Abstract. The majority of Antarctic ice shelves are bounded by grounded ice rises. These ice rises exhibit local flow fields that partially oppose the flow of the surrounding ice shelves. Formation of ice rises is accompanied by a characteristic upward-arching internal stratigraphy (“Raymond arches”), whose geometry can be analysed to infer information about past ice-sheet changes in areas where other archives such as rock outcrops are missing. Here we present an improved modelling framework to study ice-rise evolution using a satellite-velocity calibrated, isothermal, and isotropic 3-D full-Stokes model including grounding-line dynamics at the required mesh resolution (<500 m). This overcomes limitations of previous studies where ice-rise modelling has been restricted to 2-D and excluded the coupling between the ice shelf and ice rise. We apply the model to the Ekström Ice Shelf, Antarctica, containing two ice rises. Our simulations investigate the effect of surface mass balance and ocean perturbations onto ice-rise divide position and interpret possible resulting unique Raymond arch geometries. Our results show that changes in the surface mass balance result in immediate and sustained divide migration (>2.0 m yr−1) of up to 3.5 km. In contrast, instantaneous ice-shelf disintegration causes a short-lived and delayed (by 60–100 years) response of smaller magnitude (<0.75 m yr−1). The model tracks migration of a triple junction and synchronous ice-divide migration in both ice rises with similar magnitude but differing rates. The model is suitable for glacial/interglacial simulations on the catchment scale, providing the next step forward to unravel the ice-dynamic history stored in ice rises all around Antarctica.

Ice shelves around Antarctica can provide back stress for outlet glaciers and control ice sheet mass loss. They often contain narrow bands of thin ice termed ice shelf channels. Ice shelf channel morphology can be interpreted through surface depressions and exhibits junctions and deflections from flowlines. Using ice flow modeling and radar, we investigate ice shelf channels in the Roi Baudouin Ice Shelf. These are aligned obliquely to the prevailing easterly winds. In the shallow radar stratigraphy, syncline and anticline stacks occur beneath the upwind and downwind side, respectively. The structures are horizontally and vertically coherent, except near an ice shelf channel junction where patterns change structurally with depth. Deeper layers truncate near basal incisions. Using ice flow modeling, we show that the stratigraphy is ∼9 times more sensitive to atmospheric variability than to oceanic variability. This is due to the continual adjustment toward flotation. We propose that syncline-anticline pairs in the shallow stratigraphy are caused by preferential snow deposition on the windward side and wind erosion at the downwind side. This drives downwind deflection of ice shelf channels of several meters per year. The depth variable structures indicate formation of an ice shelf channel junction by basal melting. We conclude that many ice shelf channels are seeded at the grounding line. Their morphology farther seaward is shaped on different length scales by ice dynamics, the ocean, and the atmosphere. These processes act on finer (subkilometer) scales than are captured by most ice, atmosphere, and ocean models, yet the dynamics of ice shelf channels may have broader implications for ice shelf stability.