This is the first (Paper I) of three companion papers focused respectively, on the estimates of the errors in ice thickness retrieved from pulsed ground-penetrating radar (GPR) data, on how to estimate the errors at the grid points of an ice-thickness DEM, and on how the latter errors, plus the boundary delineation errors, affect the ice-volume estimates. We here present a comprehensive analysis of the various errors involved in the computation of ice thickness from pulsed GPR data, assuming they have been properly migrated. We split the ice-thickness error into independent components that can be estimated separately. We consider, among others, the effects of the errors in radio-wave velocity and timing. A novel aspect is the estimate of the error in thickness due to the uncertainty in horizontal positioning of the GPR measurements, based on the local thickness gradient. Another novel contribution is the estimate of the horizontal positioning error of the GPR measurements due to the velocity of the GPR system while profiling, and the periods of GPS refreshing and GPR triggering. Their effects are particularly important for airborne profiling. We illustrate our methodology through a case study of Werenskioldbreen, Svalbard.

Calving from tidewater glaciers and ice shelves accounts for around half the mass loss from both polar ice sheets, yet the process is not well represented in prognostic models of ice dynamics. Benn and others proposed a calving criterion appropriate for both grounded and floating glacier tongues or ice shelves, based on the penetration depth of transverse crevasses near the calving front, computed using Nye’s formula. The criterion is readily incorporated into glacier and ice-sheet models, but has not been fully validated with observations. We apply a three-dimensional extension of Benn and others’ criterion, incorporated into a full-Stokes model of glacier dynamics, to estimate the current position of the calving front of Johnsons Glacier, Antarctica. We find that two improvements to the original model are necessary to accurately reproduce the observed calving front: (1) computation of the tensile deviatoric stress opening the crevasse using the full-stress solution and (2) consideration of such a tensile stress as a function of depth. Our modelling results also suggest that Johnsons Glacier has a polythermal structure, rather than the temperate structure suggested by earlier studies.

This work combines very detailed measurements from terrestrial laser scanner (TLS), ground-based interferometry radar (GB-SAR) and ground-penetrating radar (GPR) to diagnose current conditions and to analyse the recent evolution of the Monte Perdido Glacier in the Spanish Pyrenees from 2011 to 2017. Thus, this is currently one of the best monitored small glacier (<0.5 km ² ) worldwide. The evolution of the glacier surface was surveyed with a TLS evidencing an important decline of 6.1 ± 0.3 m on average, with ice losses mainly concentrated over 3 years (2012, 2015 and 2017). Ice loss is unevenly distributed throughout the study period, with 10–15 m thinning in some areas while unchanged areas in others. GB-SAR revealed that areas with higher ice losses are those that are currently with no or very low ice motion. In contrast, sectors located beneath the areas with less ice loss are those that still exhibit noticeable ice movement (average 2–4.5 cm d ^─1 in summer, and annual movement of 9.98 ma ^─1 from ablation stakes data). GPR informed that ice thickness was generally <30 m, though locally 30–50 m. Glacier thinning is still accelerating and will lead to extinction of the glacier over the next 50 years.

We present ground-penetrating radar (GPR)—based volume calculations, with associated error estimates, for eight glaciers on Wedel Jarlsberg Land, southwestern Spitsbergen, Svalbard, and compare them with those obtained from volume-area scaling relationships. The volume estimates are based upon GPR ice-thickness data collected during the period 2004–2013. The total area and volume of the ensemble are 502.91 ± 18.60 km2 and 91.91 ± 3.12 km3, respectively. The individual areas, volumes, and average ice thickness lie within 0.37–140.99 km2, 0.01–31.98 km3, and 28–227 m, respectively, with a maximum recorded ice thickness of 619 ± 13 m on Austre Torellbreen. To estimate the ice volume of unsurveyed tributary glaciers, we combine polynomial cross-sections with a function providing the best fit to the measured ice thickness along the center line of a collection of 22 surveyed tributaries. For the time-to-depth conversion of GPR data, we test the use of a glacierwide constant radio-wave velocity chosen on the basis of local or regional common midpoint measurements, versus the use of distinct velocities for the firn, cold ice, and temperate ice layers, concluding that the corresponding volume calculations agree with each other within their error bounds

We use an automatic weather station and surface mass balance dataset spanning four melt seasons collected on Hurd Peninsula Glaciers, South Shetland Islands, to investigate the point surface energy balance, to determine the absolute and relative contribution of the various energy fluxes acting on the glacier surface and to estimate the sensitivity of melt to ambient temperature changes. Long-wave incoming radiation is the main energy source for melt, while short-wave radiation is the most important flux controlling the variation of both seasonal and daily mean surface energy balance. Short-wave and long-wave radiation fluxes do, in general, balance each other, resulting in a high correspondence between daily mean net radiation flux and available melt energy flux. We calibrate a distributed melt model driven by air temperature and an expression for the incoming short-wave radiation. The model is calibrated with the data from one of the melt seasons and validated with the data of the three remaining seasons. The model results deviate at most 140 mm w.e. from the corresponding observations using the glaciological method. The model is very sensitive to changes in ambient temperature: a 0.5 °C increase results in 56 % higher melt rates.

Calving is an important mass-loss process at ice sheet and marine-terminating glacier margins, but identifying and quantifying its principal driving mechanisms remains challenging. Hansbreen is a grounded tidewater glacier in southern Spitsbergen, Svalbard, with a rich history of field and remote sensing observations. The available data make this glacier suitable for evaluating mechanisms and controls on calving, some of which are considered in this paper. We use a full-Stokes thermomechanical 2D flow model (Elmer/Ice), paired with a crevasse-depth calving criterion, to estimate Hansbreen’s front position at a weekly time resolution. The basal sliding coefficient is re-calibrated every four weeks by solving an inverse model. We investigate the possible role of backpressure at the front (a function of ice mélange concentration) and the depth of water filling crevasses by examining the model’s ability to reproduce the observed seasonal cycles of terminus advance and retreat. Our results suggest that the ice-mélange pressure plays an important role in the seasonal advance and retreat of the ice front, and that the crevasse-depth calving criterion, when driven by modelled surface meltwater, closely replicates observed variations in terminus position. These results suggest that tidewater glacier behavior is influenced by both oceanic and atmospheric processes, and that neither of them should be ignored.