Search

The once-contiguous Ellesmere Ice Shelf, first reported in writing by European explorers in 1876, and now almost completely disintegrated, has rolling, wave-like surface topography, the origin of which we investigate using a viscous buckling instability analysis. We show that rolls can develop during a winter season (~ 100 d) if sea-ice pressure (depth-integrated horizontal stress applied to the seaward front of the Ellesmere Ice Shelf) is sufficiently large (1 MPa m) and ice thickness sufficiently low (1–10 m). Roll wavelength initially depends only on sea-ice pressure, but evolves over time depending on amplitude growth rate. This implies that a thinner ice shelf, with its faster amplitude growth rate, will have a shorter wavelength compared to a thicker ice shelf when sea-ice pressure is equal. A drawback of the viscous buckling mechanism is that roll amplitude decays once sea-ice pressure is removed. However, non-Newtonian ice rheology, where effective viscosity, and thus roll change rate, depends on total applied stress may constrain roll decay rate to be much slower than growth rate and allow roll persistence from year to year. Whether the viscous-buckling mechanism we explore here ultimately can be confirmed as the origin of the Ellesmere Ice Shelf rolls remains for future research.

Surface meltwater is becoming increasingly widespread on Antarctic ice shelves. It is stored within surface ponds and streams, or within firn pore spaces, which may saturate to form slush. Slush can reduce firn air content, increasing an ice-shelf's vulnerability to break-up. To date, no study has mapped the changing extent of slush across ice shelves. Here, we use Google Earth Engine and Landsat 8 images from six ice shelves to generate training classes using a k-means clustering algorithm, which are used to train a random forest classifier to identify both slush and ponded water. Validation using expert elicitation gives accuracies of 84% and 82% for the ponded water and slush classes, respectively. Errors result from subjectivity in identifying the ponded water/slush boundary, and from inclusion of cloud and shadows. We apply our classifier to the Roi Baudouin Ice Shelf for the entire 2013–20 Landsat 8 record. On average, 64% of all surface meltwater is classified as slush and 36% as ponded water. Total meltwater areal extent is greatest between late January and mid-February. This highlights the importance of mapping slush when studying surface meltwater on ice shelves. Future research will apply the classifier across all Antarctic ice shelves.

We develop a two-dimensional, plan-view formulation of ice-shelf flow and viscoelastic ice-shelf flexure. This formulation combines, for the first time, the shallow-shelf approximation for horizontal ice-shelf flow (and shallow-stream approximation for flow on lubricated beds such as where ice rises and rumples form), with the treatment of a thin-plate flexure. We demonstrate the treatment by performing two finite-element simulations: one of the relict pedestalled lake features that exist on some debris-covered ice shelves due to strong heterogeneity in surface ablation, and the other of ice rumpling in the grounding zone of an ice rise. The proposed treatment opens new venues to investigate physical processes that require coupling between the longitudinal deformation and vertical flexure, for instance, the effects of surface melting and supraglacial lakes on ice shelves, interactions with the sea swell, and many others.

Diurnal depth cycles of decimeter scale are observed in a supraglacial lake on the McMurdo Ice Shelf, Antarctica. We evaluate two possible causes: (1) tidal tilt of the ice shelf in response to the underlying ocean tide, and (2) meltwater input variation. We find the latter to be the most likely explanation of our observations. However, we do not rule out tidal tilt as a source of centimeter scale variations, and point to the possibility that other, larger supraglacial lake systems, particularly those on ice shelves that experience higher amplitude tidal tilts, such as in the Weddell Sea, may have depth cycles driven by ocean tide. The broader significance of diurnal cycles in meltwater depth is that, under circumstances where the ice shelf is thin, tidal-tilt amplitudes are high, and meltwater runoff rates are large, there may be associated flexure stresses that can contribute to ice-shelf fracture and destabilization. For the McMurdo Ice Shelf (~20–50 m thickness, ~ 1 m tidal amplitude and ~10 cm water-depth variations), these stresses amount to several 10's of kPa.

Seismograms acquired on the McMurdo Ice Shelf, Antarctica, during an Austral summer melt season (November 2016–January 2017) reveal a diurnal cycle of seismicity, consisting of hundreds of thousands of small ice quakes limited to a 6–12 hour period during the evening, in an area where there is substantial subsurface melting. This cycle is explained by thermally induced bending and fracture of a frozen surface superimposed on a subsurface slush/water layer that is supported by solar radiation penetration and absorption. A simple, one-dimensional model of heat transfer driven by observed surface air temperature and shortwave absorption reproduces the presence and absence (as daily weather dictated) of the observed diurnal seismicity cycle. Seismic event statistics comparing event occurrence with amplitude suggest that the events are generated in a fractured medium featuring relatively low stresses, as is consistent with a frozen surface superimposed on subsurface slush. Waveforms of the icequakes are consistent with hydroacoustic phases at frequency $ {\bf \gt} \bf 75\,{\bf Hz}$ and flexural-gravity waves at frequency $ \bf {\bf \lt}25\,{\bf Hz}$. Our results suggest that seismic observation may prove useful in monitoring subsurface melting in a manner that complements other ground-based methods as well as remote sensing.

Surface debris covers much of the western portion of the McMurdo Ice Shelf and has a strong influence on the local surface albedo and energy balance. Differential ablation between debris-covered and debris-free areas creates an unusual heterogeneous surface of topographically low, high-ablation, and topographically raised (‘pedestalled’), low-ablation areas. Analysis of Landsat and MODIS satellite imagery from 1999 to 2018, alongside field observations from the 2016/2017 austral summer, shows that pedestalled relict lakes (‘pedestals’) form when an active surface meltwater lake that develops in the summer, freezes-over in winter, resulting in the lake-bottom debris being masked by a high-albedo, superimposed, ice surface. If this ice surface fails to melt during a subsequent melt season, it experiences reduced surface ablation relative to the surrounding debris-covered areas of the ice shelf. We propose that this differential ablation, and resultant hydrostatic and flexural readjustments of the ice shelf, causes the former supraglacial lake surface to become increasingly pedestalled above the lower topography of the surrounding ice shelf. Consequently, meltwater streams cannot flow onto these pedestalled features, and instead divert around them. We suggest that the development of pedestals has a significant influence on the surface-energy balance, hydrology and flexure of the ice shelf.

Supraglacial lakes are known to trigger Antarctic ice-shelf instability and break-up. However, to date, no study has focused on lakes on Greenland's floating termini. Here, we apply lake boundary/area and depth algorithms to Landsat 8 imagery to analyse the inter- and intraseasonal evolution of supraglacial lakes across Petermann Glacier's (81°N) floating tongue from 2014 to 2016, while also comparing these lakes to those on the grounded ice. Lakes start to fill in June and quickly peak in total number, volume and area in late June/early July in response to increases in air temperatures. However, through July and August, total lake number, volume and area all decline, despite sustained high temperatures. These observations may be explained by the transportation of meltwater into the ocean by a river, and by lake drainage events on the floating tongue. Further, as mean lake depth remains relatively constant during this time, we suggest that a large proportion of the lakes that drain, do so completely, likely by rapid hydrofracture. The mean areas of lakes on the tongue are only ~20% of those on the grounded ice and exhibit lower variability in maximum and mean depth, differences likely attributable to the contrasting formation processes of lakes in each environment.

The controls on rapid surface lake drainage on the Greenland ice sheet (GrIS) remain uncertain, making it challenging to incorporate lake drainage into models of GrIS hydrology, and so to determine the ice-dynamic impact of meltwater reaching the ice-sheet bed. Here, we first use a lake area and volume tracking algorithm to identify rapidly draining lakes within West Greenland during summer 2014. Second, we derive hydrological, morphological, glaciological and surface-mass-balance data for various factors that may influence rapid lake drainage. Third, these factors are used within Exploratory Data Analysis to examine existing hypotheses for rapid lake drainage. This involves testing for statistical differences between the rapidly and non-rapidly draining lake types, as well as examining associations between lake size and the potential controlling factors. This study shows that the two lake types are statistically indistinguishable for almost all factors investigated, except lake area. Thus, we are unable to recommend an empirically supported, deterministic alternative to the fracture area threshold parameter for modelling rapid lake drainage within existing surface-hydrology models of the GrIS. However, if improved remotely sensed datasets (e.g. ice-velocity maps, climate model outputs) were included in future research, it may be possible to detect the causes of rapid drainage.

On 2 March 2016, several small en échelon tabular icebergs calved from the seaward front of the McMurdo Ice Shelf, and a previously inactive rift widened and propagated by ~3 km, ~25% of its previous length, setting the stage for the future calving of a ~14 km2 iceberg. Within 24 h of these events, all remaining land-fast sea ice that had been stabilizing the ice shelf broke-up. The events were witnessed by time-lapse cameras at nearby Scott Base, and put into context using nearby seismic and automatic weather station data, satellite imagery and subsequent ground observation. Although the exact trigger of calving and rifting cannot be identified definitively, seismic records reveal superimposed sets of both long-period (>10 s) sea swell propagating into McMurdo Sound from storm sources beyond Antarctica, and high-energy, locally-sourced, short-period (

Global Navigation Satellite System (GNSS) observations and ground-based timelapse photography obtained over the record-high 2019/2020 melt season are combined to characterise the flexure and fracture behaviour of a previously formed doline on George VI Ice Shelf, Antarctica. The GNSS timeseries shows a downward vertical displacement of the doline centre with respect to the doline rim of ~60 cm in response to loading from a central meltwater lake. The GNSS data also show a tens-of-days episode of rapid-onset, exponentially decaying horizontal displacement, where the horizontal distance between the doline rim and its centre increases by ~70 cm. We interpret this event as the initiation and/or widening of a fracture, aided by stress perturbations associated with meltwater loading in the doline basin. Viscous flexure modelling indicates that the meltwater loading generates tensile surface stresses exceeding 75 kPa. This, together with our timelapse photos of circular fractures around the doline, suggests the first such documentation of meltwater-loading-induced ‘ring fracture’ formation on an ice shelf, equivalent to the fracture type proposed as part of the chain-reaction lake drainage process involved in the 2002 breakup of the Larsen B Ice Shelf.