Mass loss from iceberg calving at marine-terminating glaciers is one of the largest and most poorly constrained contributors to sea-level rise. However, our understanding of the processes controlling ice fracturing and crevasse evolution is incomplete. Here, we use Gabor filter banks to automatically map crevasse density and orientation through time on a ~150 km terminus region of Narsap Sermia, an outlet glacier of the southwest Greenland ice sheet. We find that Narsap Sermia is dominated by transverse (flow-perpendicular) crevasses near the ice front and longitudinal (flow-aligned) crevasses across its central region. Measured crevasse orientation varies on sub-annual timescales by more than 45 in response to seasonal velocity changes, and also on multi-annual timescales in response to broader dynamic changes and glacier retreat. Our results show a gradual up-glacier propagation of the zone of flow-transverse crevassing coincident with frontal retreat and acceleration occurring in 2020/21, in addition to sub-annual crevasse changes primarily in transition zones between longitudinal to transverse crevasse orientation. This provides new insight into the dynamics of crevassing at large marine-terminating glaciers and a potential approach for the rapid identification of glacier dynamic change from a single pair of satellite images.

Changes in Greenland's marine-terminating outlet glaciers have led to changes in the flux of icebergs into Greenland's coastal waters, yet icebergs remain a relatively understudied component of the ice-ocean system. We developed a simple iceberg delineation algorithm for Landsat imagery. A machine learning-based cloud mask incorporated into the algorithm enables us to extract iceberg size distributions from open water even in partially cloudy scenes. We applied the algorithm to the Landsat archive covering Disko Bay, West Greenland, to derive a time series of iceberg size distributions from 2000-02 and 2013-15. The time series captures a change in iceberg size distributions, which we interpret as a result of changes in the calving regime of the parent glacier, Sermeq Kujalleq (Jakobshavn Isbræ). The change in calving style associated with the disintegration and disappearance of Sermeq Kujalleq's floating ice tongue resulted in the production of more small icebergs. The increased number of small icebergs resulted in increasingly negative power law slopes fit to iceberg size distributions in Disko Bay, suggesting that iceberg size distribution time series provide useful insights into changes in calving dynamics.

While numerous maps of Greenland ice flow velocity exist, most have gaps in coverage and/or accuracy is limited. We processed a large volume of synthetic aperture radar (SAR) and Landsat 8 imagery collected between 1995 and 2015 to produce a nearly complete map of ice flow velocity for Greenland at a far greater accuracy than most prior products. We evaluated the accuracy of this map by comparing it with a variety of measured and estimated velocities. For the slow-moving interior of the ice sheet, where estimates are determined from interferometric phase, the errors are ~2 m a or better. For coastal areas, where estimates are determined entirely from speckle- or feature-tracking methods, errors are 2-3 m a, which is in good agreement with the estimated formal errors. Especially for the slow-moving majority of the ice sheet, this map provides an important source of data for numerous types of glaciological studies.

Glaciers in Greenland are changing rapidly. To better understand these changes, we have produced a series of seven synthetic-aperture-radar (SAR) backscatter mosaics for seven winters during the period from 2000 to 2013. Six of the mosaics were created using RADARSAT Fine-Beam data and the seventh used ALOS PALSAR Fine-Beam Single-Polarization data. The RADARSAT mosaics are radiometrically calibrated and capture changes in the backscatter coefficient related to melt and other events, particularly the strong melting in the summer of 2012. Comparison of features in the ascending-orbit ALOS mosaic and the descending-orbit RADARSAT mosaics indicate that in areas of smooth to moderate topography their locations are consistent to within a few 10s of meters. The locations of features identifiable in the RADARAT mosaics, which were collected with the same imaging parameters, generally agree to within better than the 20-m posting of the data. With such geometric accuracy, these data establish a record of change in Greenland for the early part of the 21 Century, thus providing a baseline that can be compared with new radar and optical data sets.

Increasing ocean and air temperatures have contributed to the removal of floating ice shelves from several Greenland outlet glaciers; however, the specific contribution of these external forcings remains poorly understood. Here we use atmospheric, oceanographic, and glaciological time series data from the ice shelf of Petermann Gletscher, NW Greenland to quantify the forcing of the ocean and atmosphere on the ice shelf at a site ~16 km from the grounding line within a large sub-ice-shelf channel. Basal melt rates here indicate a strong seasonality, rising from a winter mean of 2 m yr to a maximum of 80 m yr during the summer melt season. This increase in basal melt rates confirms the direct link between summer atmospheric warming around Greenland and enhanced ocean-forced melting of its remaining ice shelves. We attribute this enhanced melting to increased discharge of subglacial runoff into the ocean at the grounding line, which strengthens under-ice currents and drives a greater ocean heat flux toward the ice base.