We review the findings of approximately 60 years of in situ and remote sensing studies of glacier crevasses, as well as the three broad classes of numerical models now employed to simulate crevasse fracture. The relatively new insight that mixed-mode fracture in local stress equilibrium, rather than downstream advection alone, can introduce nontrivial curvature to crevasse geometry may merit the reinterpretation of some key historical observation studies. In the past three decades, there have been tremendous advances in the spatial resolution of satellite imagery, as well as fully automated algorithms capable of tracking crevasse displacements between repeat images. Despite considerable advances in developing fully transient three-dimensional ice flow models over the past two decades, both the zero stress and linear elastic fracture mechanics crevasse models have remained fundamentally unchanged over this time. In the past decade, however, multidimensional and transient formulations of the continuum damage mechanics approach to simulating ice fracture have emerged. The combination of employing damage mechanics to represent slow upstream deterioration of ice strength and fracture mechanics to represent rapid failure at downstream termini holds promise for implementation in large-scale ice sheet models. Finally, given the broad interest in the sea level rise implications of recent and future cryospheric change, we provide a synthesis of 10 mechanisms by which crevasses can influence glacier mass balance. Field ObservationsBy virtue of providing relatively easy access to a glacier's interior, crevasses can be an asset to observational glaciologists, permitting, for example, the efficient examination of glacier stratigraphy. Annual accumulation

No abstract

Liquid water stored on the surface of ice sheets and glaciers impacts surface mass balance, ice dynamics, and heat transport. Multispectral remote sensing can be used to detect supraglacial lakes and estimate their depth and area. In this study, we use in situ spectral and bathymetric data to assess lake depth retrieval using the recently launched Landsat 8 Operational Land Imager (OLI). We also extend our analysis to other multispectral sensors to evaluate their performance with similar methods. Digital elevation models derived from WorldView stereo imagery (pre-lake filling and post-drainage) are used to validate spectrally derived depths, combined with a lake edge determination from imagery. The optimal supraglacial lake depth retrieval is a physically based single-band model applied to two OLI bands independently (red and panchromatic) that are then averaged together. When OLI-and WorldView-derived depths are differenced, they yield a mean and standard deviation of 0.0 ± 1.6 m. This method is then applied to OLI data for the Sermeq Kujalleq (Jakobshavn Isbrae) region of Greenland to study the spatial and intra-seasonal variability of supraglacial lakes during summer 2014. We also give coefficients for estimating supraglacial lake depth using a similar method with other multispectral sensors.

Melt and supraglacial lakes are precursors to ice shelf collapse and subsequent accelerated ice sheet mass loss. We used data from the Landsat 8 and Sentinel-2 satellites to develop a threshold-based method for detection of lakes found on the Antarctic ice shelves, calculate their depths and thus their volumes. To achieve this, we focus on four key areas: the Amery, Roi Baudouin, Nivlisen, and Riiser-Larsen ice shelves, which are all characterized by extensive surface meltwater features. To validate our products, we compare our results against those obtained by an independent method based on a supervised classification scheme (e.g., Random Forest algorithm). Additional verification is provided by manual inspection of results for nearly 1000 Landsat 8 and Sentinel-2 images. Our dual-sensor approach will enable constructing high-resolution time series of lake volumes. Therefore, to ensure interoperability between the two datasets, we evaluate depths from contemporaneous Landsat 8 and Sentinel-2 image pairs. Our assessments point to a high degree of correspondence, producing an average R2 value of 0.85, no bias, and an average RMSE of 0.2 m. We demonstrate our method’s ability to characterize lake evolution by presenting first evidence of drainage events outside of the Antarctic Peninsula on the Amery Ice shelf. The methods presented here pave the way to upscaling throughout the Landsat 8 and Sentinel-2 observational record across Antarctica to produce a first-ever continental dataset of supraglacial lake volumes. Such a dataset will improve our understanding of the influence of surface hydrology on ice shelf stability, and thus, future projections of Antarctica’s contribution to sea level rise.