Climate and Surging of Donjek Glacier, Yukon, Canada

Links between climate and glacier surges are poorly understood but are required to enable prediction of surges and mitigation of associated hazards. Here, we investigate the role of snow accumulation, rain, and temperature on surge periodicity, area changes, and timing of surge initiation since the 1930s at Donjek Glacier, Yukon, Canada. Snow accumulation measured in three ice cores collected at Eclipse Icefield indicates that a cumulative accumulation of 15.5 ± 1.46 or 16.6 ± 2.0 m w.e. occurred in the ten to twelve years between each of its last eight surges, depending on ice motion spatiotemporal offset corrections. Although we find consistent snow accumulation between surges, the transient snow line has risen 10.3 m decade−1 vertically since the 1950s, and Burwash Landing weather station records indicate a 0.5°C decade−1 increase in mean annual air temperature since the 1960s. Changes in surface mass balance are accompanied by a consistent surge interval but decreasing surge extent. The three recent surge events initiated in years with the rainiest summers on record. These findings highlight a complex interplay between external (i.e., climate) and internal glacier processes that control surging at Donjek Glacier, with climate having a more direct influence on surge extent than on recurrence interval

Dynamic Mass Loss from Greenland\u27s Marine-Terminating Peripheral Glaciers (1985–2018)

Global glacier mass balance decreased rapidly over the last two decades, exceeding mass loss from the Greenland and Antarctic Ice Sheets. In Greenland, peripheral glaciers and ice caps (GICs) cover only ~5% of Greenland\u27s area but contributed ~20% of the island\u27s ice mass loss between 2000 and 2018. Although Greenland GIC mass loss due to surface meltwater runoff has been estimated using atmospheric models, mass lost to changes in ice discharge into oceans (i.e., dynamic mass loss) remains unquantified. We use the flux gate method to estimate discharge from Greenland\u27s 585 marine-terminating peripheral glaciers between 1985 and 2018, and compute dynamic mass loss as the discharge anomaly relative to the 1985–98 period. Greenland GICs discharged between 2.94 ± 0.23 and 4.03 ± 0.23 Gt a−1 from 1985 to 1998, depending on the gap-filling method, and abruptly increased to 5.10 ± 0.21 Gt a−1 from 1999 to 2018. The resultant ~1–2 Gt a−1 dynamic mass loss was driven by synchronous widespread acceleration around Greenland. The mass loss came predominantly from the southeast region, which contains 39% of the glaciers. Although changes in discharge over time were small relative to surface mass-balance changes, our speed and discharge time series suggest these glaciers may quickly accelerate in response to changes in climate

Brief Communication: Is Vertical Shear in an Ice Shelf (Still) Negligible?

Vertical shear is recognized today as a key component of the stress balance of ice shelves. However, the first ice shelf models were built on the neglect of vertical shear. Partly due to its historical treatment, it remains common to discuss vertical shear as though it were still considered negligible in ice shelf models. Here, we offer a historical perspective on the changing treatment of vertical shear over time, and we emphasize the term\u27s non-negligibility in current ice shelf modeling. We illustrate our discussion in the simplest context of an analytic, isothermal, shallow-ice-shelf model

Automated Detection of Marine Glacier Calving Fronts Using the 2-D Wavelet Transform Modulus Maxima Segmentation Method

Changes in the calving front position of marine-terminating glaciers strongly influence the mass balance of glaciers, ice caps, and ice sheets. At present, quantification of frontal position change primarily relies on time-consuming and subjective manual mapping techniques, limiting our ability to understand changes to glacier calving fronts. Here we describe a newly developed automated method of mapping glacier calving fronts in satellite imagery using observations from a representative sample of Greenland’s peripheral marine-terminating glaciers. Our method is adapted from the 2-D wavelet transform modulus maxima (WTMM) segmentation method, which has been used previously for image segmentation in biomedical and other applied science fields. The gradient-based method places edge detection lines along regions with the greatest intensity gradient in the image, such as the contrast between glacier ice and water or glacier ice and sea ice. The lines corresponding to the calving front are identified using thresholds for length, average gradient value, and orientation that minimize the misfit with respect to a manual validation data set. We demonstrate that the method is capable of mapping glacier calving fronts over a wide range of image conditions (light to intermediate cloud cover, dim or bright, mélange presence, etc.). With these time series, we are able to resolve subseasonal to multiyear temporal patterns as well as regional patterns in glacier frontal position change

Iceberg properties and distributions in three Greenlandic fjords using satellite imagery

Icebergs calved from tidewater glaciers represent about one third to one half of the freshwater flux from the Greenland ice sheet to the surrounding ocean. Using multiple satellite datasets, we quantify the first fjord-wide distributions of iceberg sizes and characteristics for three fjords with distinct hydrography and geometry: Sermilik Fjord, Rink Isbræ Fjord and Kangerlussuup Sermia Fjord. We estimate average total iceberg volumes in summer in the three fjords to be 6.4 ± 1.5, 1.7 ± 0.40 and 0.16 ± 0.09 km3, respectively. Iceberg properties are influenced by glacier calving style and grounding line depth, with variations in size distribution represented by exponents of power law distributions that are −1.95 ± 0.06, −1.87 ± 0.05 and −1.62 ± 0.04, respectively. The underwater surface area of icebergs exceeds the subsurface area of glacial termini by at least one order of magnitude in all three fjords, underscoring the need to include iceberg melt in fjord freshwater budgets. Indeed, in Sermilik Fjord, we calculate summertime freshwater flux from iceberg melt of 620 m3 s−1 (±140 m3 s−1), similar in magnitude to subglacial discharge. The method developed here can be extended across Greenland to assess relationships between glacier calving, iceberg discharge and freshwater production.NNX12AP50G55223

Assessing Controls on Ice Dynamics at Crane Glacier, Antarctic Peninsula, Using a Numerical Ice Flow Model

The Antarctic Peninsula\u27s widespread glacier retreat and ice shelf collapse have been attributed to atmospheric and oceanic warming. Following the initial post-collapse period of retreat, several former tributary glaciers of the Larsen A and B ice shelves have been slowly re-advancing for more than a decade. Here, we use a flowline model of Crane Glacier to gauge the sensitivity of former tributary glaciers to future climate change following this period of long-term dynamic adjustment. The glacier\u27s long-term geometry and speed changes are similar to those of other former Larsen A and B tributaries, suggesting that Crane Glacier is a reasonable representation of regional dynamics. For the unperturbed climate simulations, discharge remains nearly unchanged in 2018–2100, indicating that dynamic readjustment to shelf collapse took ~15 years. Despite large uncertainties in Crane Glacier\u27s past and future climate forcing, a wide range of future climate scenarios leads to a relatively modest range in grounding line discharge (0.8–1.5 Gt a−1) by 2100. Based on the model results for Crane, we infer that although former ice shelf tributaries may readvance following collapse, similar to the tidewater glacier cycle, their dynamic response to future climate perturbations should be less than their response to ice shelf collapse

Uncertainty of ICESat-2 ATL06- and ATL08-Derived Snow Depths for Glacierized and Vegetated Mountain Regions

Seasonal snow melt dominates the hydrologic budget across a large portion of the globe. Snow accumulation and melt vary over a broad range of spatial scales, preventing accurate extrapolation of sparse in situ observations to watershed scales. The lidar onboard the Ice, Cloud, and land Elevation, Satellite (ICESat-2) was designed for precise mapping of ice sheets and sea ice, and here we assess the feasibility of snow depth-mapping using ICESat-2 data in more complex and rugged mountain landscapes. We explore the utility of ATL08 Land and Vegetation Height and ATL06 Land Ice Height differencing from reference elevation datasets in two end member study sites. We analyze ∼3 years of data for Reynolds Creek Experimental Watershed in Idaho\u27s Owyhee Mountains and Wolverine Glacier in southcentral Alaska\u27s Kenai Mountains. Our analysis reveals decimeter-scale uncertainties in derived snow depth and glacier mass balance at the watershed scale. Both accuracy and precision decrease as slope increases: the magnitudes of the median and median of the absolute deviation of elevation errors (MAD) vary from ∼0.2 m for slopes \u3c 5° to \u3e 1 m for slopes \u3e 20°. For glacierized regions, failure to account for intra- and inter-annual evolution of glacier surface elevations can strongly bias ATL06 elevations, resulting in under-estimation of the mass balance gradient with elevation. Based on these results, we conclude that ATL08 and ATL06 observations are best suited for characterization of watershed-scale snow depth and mass balance gradients over relatively shallow slopes with thick snowpacks. In these regions, ICESat-2 elevation residual-derived snow depth and mass balance transects can provide valuable watershed scale constraints on terrain parameter- and model-derived estimates of snow accumulation and melt

Ice mass change in Greenland and Antarctica between 1993 and 2013 from satellite gravity measurements

We construct long-term time series of Greenland and Antarctic ice sheet mass change from satellite gravity measurements. A statistical reconstruction approach is developed based on a Principal Component Analysis to combine high-resolution spatial modes from the Gravity Recovery and Climate Experiment (GRACE) mission with the gravity information from conventional satellite track-ing data. Uncertainties of this reconstruction are rigorously assessed; they include temporal limitations for short GRACE measurements, spatial limitations for the low-resolution conventional tracking data measurements, and limitations of the estimated statistical relationships between low and high degree potential coe�cients re ected in the PCA modes. Trends of mass variations in Greenland and Antarctica are assessed against a number of previous studies. The resulting time series for Greenland show a higher rate of mass loss than other methods before 2000, while the Antarctic ice sheet appears heavily in uenced by interannual variations

Sharp Contrasts in Observed and Modeled Crevasse Patterns at Greenland’s Marine Terminating Glaciers

Crevasses are affected by and affect both the stresses and the surface mass balance of glaciers. These effects are brought on through potentially important controls on meltwater routing, glacier viscosity, and iceberg calving, yet there are few direct observations of crevasse sizes and locations to inform our understanding of these interactions. Here we extract depth estimates for the visible portion of crevasses from high-resolution surface elevation observations for 52 644 crevasses from 19 Greenland glaciers. We then compare our observed depths with those calculated using two popular models that assume crevasse depths are functions of local stresses: the Nye and linear elastic fracture mechanics (LEFM) formulations. When informed by the observed crevasse depths, the LEFM formulation produces kilometer-scale variations in crevasse depth, in decent agreement with observations. However, neither formulation accurately captures smaller-scale variations in the observed crevasse depths. Critically, we find that along-flow patterns in crevasse depths are unrelated to along-flow patterns in strain rates (and therefore stresses). Cumulative strain rate is moderately more predictive of crevasse depths at the majority of glaciers. Our reliance on lidar limits the inference we can make regarding fracture depths. However, given the discordant patterns in observed and modeled crevasses, we recommend additional in situ and remote sensing analyses before Nye and LEFM models are considered predictive. Such analyses should span extensional and compressive regions to better understand the influence of advection on crevasse geometry. Ultimately, such additional study will enable more reliable projection of terminus position change and supraglacial meltwater routing that relies on accurate modeling of crevasse occurrence