Ice-marginal glaciated landscapes demarcate former boundaries of the continental ice sheets. Throughout circumpolar regions, permafrost has preserved relict ground ice and glacigenic sediments, delaying the sequence of postglacial landscape change that transformed temperate environments millennia earlier. Here we show that within 7 × 10 6 km 2 of glaciated permafrost terrain, extensive landscapes remain poised for major climate-driven change. Across northwestern Canada, 60-100-km-wide concentric swaths of thaw slump-affected terrain delineate the maximum and recessional positions of the Laurentide Ice Sheet. These landscapes comprise ~17% of continuous permafrost terrain in a 1.27 × 10 6 km 2 study area, indicating widespread preservation of late Pleistocene ground ice. These thaw slump, relict ground ice, and glacigenic terrain associations are also evident at the circumpolar scale. Recent intensification of thaw slumping across northwestern Canada has mobilized primary glacial sediments, triggering a cascade of fluvial, lacustrine, and coastal effects. These geologically significant processes, highlighted by the spatial distribution of thaw slumps and patterns of fluvial sediment mobilization, signal the climate-driven renewal of deglaciation and postglacial permafrost landscape evolution.

No abstract

Spring melt pond fraction (f p) has been shown to influence September sea ice extent and, with a growing need to improve melt pond physics in climate and forecast models, observations at large spatial scales are needed. We present a novel technique for estimating f p on sea ice at high spatial resolution from the Sentinel-1 satellite during the winter period leading up to spring melt. A strong correlation (r = À0.85) is found between winter radar backscatter and f p from first-year and multiyear sea ice data collected in the Canadian Arctic Archipelago (CAA) in 2015. Observations made in the CAA in 2016 are used to validate a f p retrieval algorithm, and a f p prediction for the CAA in 2017 is made. The method is effective using the horizontal transmit and receive polarization channel only and shows promise for providing seasonal, pan-Arctic, f p maps for improved understanding of melt pond distributions and forecast model skill. Plain Language Summary Recent and well-documented changes in Arctic sea ice have introduced the need for timely and accurate seasonal forecasts of ice conditions. Seasonal forecasts of ice conditions will reduce the risks to humans and help preserve the fragile Arctic ecosystem by preventing accidents and spills. Recent studies have shown a link between the amount of surface meltwater flooding that occurs on sea ice in the spring, termed melt pond fraction, and the extent of sea ice that remains at the end of summer. This link is due to the ability of surface meltwater to absorb more sunlight compared to bare ice and snow. This study provides a new way to estimate the amount of surface meltwater flooding expected to occur on the sea ice in spring, using satellite data collected during the winter period. The results presented here provide a key link between winter and late summer sea ice conditions that will enhance the ability of forecasters to make accurate seasonal predictions several months in advance of the active summer period.

Abstract Two satellite datasets are used to characterize winter landfast first-year sea-ice (FYI), deformed FYI (DFYI) and multiyear sea-ice (MYI) roughness in the Canadian Arctic Archipelago (CAA): (1) optical Multi-angle Imaging SpectroRadiometer (MISR) and (2) synthetic aperture radar Sentinel-1. The Normalized Difference Angular Index (NDAI) roughness proxy derived from MISR, and backscatter from Sentinel-1 are intercompared. NDAI and backscatter are also compared to surface roughness derived from an airborne LiDAR track covering a subset of FYI and MYI (no DFYI). Overall, NDAI and backscatter are significantly positively correlated when all ice type samples are considered. When individual ice types are evaluated, NDAI and backscatter are only significantly correlated for DFYI. Both NDAI and backscatter are correlated with LiDAR-derived roughness (r = 0.71 and r = 0.74, respectively). The relationship between NDAI and roughness is greater for MYI than FYI, whereas for backscatter and ice roughness, the relationship is greater for FYI than MYI. Linear regression models are created for the estimation of FYI and MYI roughness from NDAI, and FYI roughness from backscatter. Results suggest that using a combination of Sentinel-1 backscatter for FYI and MISR NDAI for MYI may be optimal for mapping winter sea-ice roughness in the CAA.

Snow thickness measurements over relatively smooth Arctic first-year sea ice, obtained near Cambridge Bay in the Canadian Arctic (2014, 2016 and 2017) and near Elson Lagoon in the Alaskan Arctic (2003 and, are analyzed to quantify physical length-scales and their relevant scaling behaviors. We use the multi-fractal temporally weighted detrended fluctuation analysis method to detect two major physical length-scales from the two independent study locations. Our results suggest that physical processes underlying the formation of snow dunes are consistent and that the wind is the main process shaping the snow thickness variability and redistribution. One scale, around 10 m, appears to be related to the formation of the snow 'dunes', while the other scale, between 30 and 100 m, is likely associated with the various interactions of the snow dunes such as merging, calving and lateral linking. Results imply that snow on level sea ice shows self-organized characteristics.