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

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.

Supraglacial lakes on the Greenland Ice Sheet are expanding inland, but the impact on ice flow is equivocal because interior surface conditions may preclude the transfer of surface water to the bed. Here we use a well-constrained 3D model to demonstrate that supraglacial lakes in Greenland drain when tensile-stress perturbations propagate fractures in areas where fractures are normally absent or closed. These melt-induced perturbations escalate when lakes as far as 80 km apart form expansive networks and drain in rapid succession. The result is a tensile shock that establishes new surface-to-bed hydraulic pathways in areas where crevasses transiently open. We show evidence for open crevasses 135 km inland from the ice margin, which is much farther inland than previously considered possible. We hypothesise that inland expansion of lakes will deliver water and heat to isolated regions of the ice sheet’s interior where the impact on ice flow is potentially large.

Supraglacial meltwater lakes on the western Greenland Ice Sheet (GrIS) are critical components of its surface hydrology and surface mass balance, and they also affect its ice dynamics. Estimates of lake volume, however, are limited by the availability of in situ measurements of water depth, which in turn also limits the assessment of remotely sensed lake depths. Given the logistical difficulty of collecting physical bathymetric measurements, methods relying upon in situ data are generally restricted to small areas and thus their application to largescale studies is difficult to validate. Here, we produce and validate spaceborne estimates of supraglacial lake volumes across a relatively large area (1250 km 2 ) of west Greenland's ablation region using data acquired by the WorldView-2 (WV-2) sensor, making use of both its stereo-imaging capability and its meter-scale resolution. We employ spectrally-derived depth retrieval models, which are either based on absolute reflectance (singlechannel model) or a ratio of spectral reflectances in two bands (dual-channel model). These models are calibrated by using WV-2 multispectral imagery acquired early in the melt season and depth measurements from a high resolution WV-2 DEM over the same lake basins when devoid of water. The calibrated models are then validated with different lakes in the area, for which we determined depths. Lake depth estimates based on measurements recorded in WV-2's blue (450-510 nm), green (510-580 nm), and red (630-690 nm) bands and dual-channel modes (blue/green, blue/red, and green/red band combinations) had near-zero bias, an average root-meansquared deviation of 0.4 m (relative to post-drainage DEMs), and an average volumetric error of b1%. The approach outlined in this study -image-based calibration of depth-retrieval models -significantly improves spaceborne supraglacial bathymetry retrievals, which are completely independent from in situ measurements.