Piotr Glowacki

High-resolution hydrothermal structure of Hansbreen, Spitsbergen, mapped by ground-penetrating radar

Detailed ground-penetrating radar (GPR) surveys at 50 and 200 MHz on Hansbreen, a polythermal glacier in southern Svalbard, are presented and interpreted. Comparison of the variations in character of the radar reflections with borehole thermometry and water levels in moulins suggests that GPR can be used to study the hydrothermal properties of the glacier. The high resolution of the GPR data shows that the hydrothermal structure of the glacier is highly variable both along the centre line and on transverse profiles. Water contents for many places and depths within the glacier were calculated by estimating radar-wave velocities to point reflectors. We find typical water contents of 1-2% for the temperate ice, but wetter ice associated with surface crevassing and moulins (typically 4% water content). There is evidence that wet ice sometimes overlays drier ice. The hydrothermal structure is thus shown to be very complex. Temperature gradients in the cold ice indicate freezing rates of temperate ice below cold ice of 0.1 0.5 m a, while isolated point reflectors within the cold ice indicate large water-filled bodies that are probably related to the regular drainage structure of the glacier.

The drainage pattern of Hansbreen and Werenskioldbreen, two polythermal glaciers in Svalbard

To improve our understanding of Svalbard-type polythermal glacier drainage, hydraulic geometry models of the subglacial hydrology of two contrasting glaciers in Svalbard have been constructed. The models are tested against a uniquely long and rich set of field observations spanning 45 years. Digital elevation models (DEMs) were constructed from bedrock data measured with ground penetrating radar and surface data of two medium-sized polythermal glaciers, Hansbreen and Werenskioldbreen, in south-west Spitsbergen. Hansbreen has a low angle bed with over-deepenings and a calving front, while Werenskioldbreen has steeper bed and terminates on land. Together they are representative of many Svalbard glaciers. The DEMs were used to derive maps of hydraulic potential and subglacial drainage networks. Validation of the models was done using field observations including location mapping and speleological exploration of active moulins, positions of main river outflows, dyetracing and water chemistry studies, and observations of water pressure inside moulins. Results suggest that the water pressure is generally close to ice overburden pressure but varies greatly depending on local conditions such as bed location, the thickness of cold ice layer, the thickness of the glacier and seasonal changes in meltwater input.

Assessment of interannual variations in the surface mass balance of 18 Svalbard glaciers from the Moderate Resolution Imaging Spectroradiometer/Terra albedo product

[1] We estimate annual anomalies of the surface mass balance of glaciers on Svalbard for the period 2000–2005 (six years), by calculating the so-called ‘‘satellite-derived mass balance’’ (Bsat) from time series of satellite-derived surface albedos. The method needs no other input variables. Surface albedos are extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra albedo product. We validate the MODIS albedos by comparing them with in situ measurements on Kongsvegen, and we find a low root-meansquare error of 0.04 for higher-quality MODIS data. Confidence in the MODIS product is also provided by realistic profiles of albedo along glacier centerlines. We apply the method to 18 glaciers that are evenly distributed over the archipelago. Correlation coefficients of time series of Bsat and direct measurements of the annual mass balance on Kongsvegen and Hansbreen are highly significant (0.94 and 0.82, respectively). Moreover, spatial distributions of the anomalies for individual years are coherent. Disadvantages of the method are that absolute values of the mass balance cannot be determined and that the interannual variability is underestimated. The latter might be corrected by equations to be established with mass balance models.

Glacier changes in southern Spitsbergen, Svalbard,1901–2000

Abstract High-resolution ground-penetrating radar surveys at 50 MHz on the polythermal glaciers Hornbreen, Hambergbreen and several surrounding glaciers in southern Spitsbergen, Svalbard, are presented and interpreted. Accurate positioning was obtained using differential global positioning system (DGPS). Digital elevation models (DEMs) of the bedrock and surface were constructed. Comparison of DGPS data and surface DEMs with data from the topographic mappings from 1936 oblique stereoscopic aerial photographs and from Mission Russe in 1899–1901 shows that the Hornbreen and Hambergbreen surfaces are about 60–100 m thinner today in the upper part than at the beginning of the 20th century. Hornbreen has retreated by 13.5 km from the central part of the front, and Hambergbreen by 16 km. All the fronts of the nearby east-coast glaciers in this area have retreated. The bedrock DEM shows that the Hornbreen and Hambergbreen beds lie at –25 to 25 m a.s.l. The combination of sub-sea-level fronts and increasing steepness of the glaciers suggests that the low-lying glaciated valley filled by Hornbreen and Hambergbreen may become a partially inundated ice-free isthmus within perhaps 100 years.

Temporal changes in the radiophysical properties of a polythermal glacier in Spitsbergen

Abstract In order to study the seasonal and inter-seasonal variations in radio-wave velocity (RWV), radiophysical investigations were made at Hansbreen, a polythermal glacier in Spitsbergen, in July– August 2003 and April 2004. These investigations included repeated radar profiling (20 and 25 MHz) along a transverse profile, repeated common-midpoint measurements, continuous radar measurements during 8 days at a fixed site, meteorological observations, and continuous ice surface velocity monitoring by differential GPS. Seasonal and inter-seasonal RWV changes in the temperate ice layer are attributed, respectively, to rapid water redistribution within it during the summer, and to variations in water content from 2.1% in summer to 0.4% in spring. The reflection properties of the temperate ice layer correlate well with the air temperature, with a nearly semi-diurnal time lag. The temporal variability of the reflection properties of the internal horizon suggests enlargement of water inclusions or water drainage from the horizon. Repeated profiling shows a stable spatial pattern in bed reflection power interpreted as changes in water content controlled by bedrock topography. The spatial variations of internal reflection energy along the repeated profile correlate with the thickness of the cold ice layer and the occurrence of drainage and crevasse systems.

Adjustment of regional climate model output for modeling the climatic mass balance of all glaciers on Svalbard

Abstract Large‐scale modeling of glacier mass balance relies often on the output from regional climate models (RCMs). However, the limited accuracy and spatial resolution of RCM output pose limitations on mass balance simulations at subregional or local scales. Moreover, RCM output is still rarely available over larger regions or for longer time periods. This study evaluates the extent to which it is possible to derive reliable region‐wide glacier mass balance estimates, using coarse resolution (10 km) RCM output for model forcing. Our data cover the entire Svalbard archipelago over one decade. To calculate mass balance, we use an index‐based model. Model parameters are not calibrated, but the RCM air temperature and precipitation fields are adjusted using in situ mass balance measurements as reference. We compare two different calibration methods: root mean square error minimization and regression optimization. The obtained air temperature shifts (+1.43°C versus +2.22°C) and precipitation scaling factors (1.23 versus 1.86) differ considerably between the two methods, which we attribute to inhomogeneities in the spatiotemporal distribution of the reference data. Our modeling suggests a mean annual climatic mass balance of −0.05 ± 0.40 m w.e. a−1 for Svalbard over 2000–2011 and a mean equilibrium line altitude of 452 ± 200 m  above sea level. We find that the limited spatial resolution of the RCM forcing with respect to real surface topography and the usage of spatially homogeneous RCM output adjustments and mass balance model parameters are responsible for much of the modeling uncertainty. Sensitivity of the results to model parameter uncertainty is comparably small and of minor importance.

Modeling and Hindcasting of the Mass Balance of Werenskioldbreen (Southern Svalbard)

Abstract The authors propose a model of glacial mass balance based on correlations with meteorological observations and data from climate re-analysis. The minimum input data required include the following: average monthly temperature on the glacier and in its vicinity during summertime for a reference time period, average monthly air temperature, and average precipitation total at the nearest weather station or from re-analysis. This model was used to hindcast the mass balance and its components at Werenskioldbreen (southern Svalbard) over the period 1912–2005. The hindcast specific mass balance was then used to estimate the change in the thickness of the snout of Werenskioldbreen over the period 1958–1990. These results were compared with results obtained using a cartographic method. Comparing the topographic maps, the glacier front lowered 28.7 m on average over 32 years. The average difference in the calculation of the change in glacier thickness between these two methods amounted to 3.7 m (based on meteorological data) and 0.2 m (using ERA-40). The discrepancy of less than 13% confirmed that the method is a reasonably accurate way of predicting past glacier mass balance. The proposed method can find a broad application in hindcasting the mass balances of small Svalbard glaciers where observation data are scarce or nonexistent.

Ice volume changes (1936-1990-2007) and ground-penetrating radar studies of Ariebreen, Hornsund, Spitsbergen

Ariebreen is a small (0.37 km2)-valley glacier located in southern Spitsbergen. Our ground-penetrating radar surveys of the glacier show that it is less than 30 m thick on average, with a maximum thickness of 82 m, and it appears to be entirely cold. By analysing digital terrain models of the ice surface from different dates, we determine the area and volume changes during two periods, 1936–1990 and 1990–2007. The total ice volume of the glacier has decreased by 73% during the entire period 1936–2007, which is equivalent to a mean mass balance rate of −0.61±0.17 m y−1 w.eq. The glacier thinning rate has increased markedly between the first and second periods, from −0.50±0.22 to −0.95±0.17 m y−1 w.eq. To access the supplementary material for this article, please see Supplementary Files under Article Tools online.

Coast formation in an Arctic area due to glacier surge and retreat: the Hornbreen - Hambergbreen case from Spistbergen

Glacierized coasts undergo faster geomorphic processes than unglaciated ones. We have studied changes of the coastal area in southern Svalbard with the glacier bridge between Torell Land and Sorkapp Land since the beginning of the 20th century. The existence of a continuous subglacial depression beneath the Hornbreen–Hambergbreen system has been debated since the 1960s, with inconclusive results. In this study we assess both the subglacial topography and the bathymetry of Hornsund Fjord and Hambergbukta bay. This included ~40 km of radar surveys over the glacial system and sea depth sounding. The extent of the glaciers from maps and satellite images together with digital terrain models and surface elevation data based on GPS profiling were used to analyse geometry changes of the glacier surfaces. The results confirm the existence a continuous subglacial depression below sea level (c. 40 m deep) between Hornsund and the Barents Sea. The Hornbreen-Hambergbreen system has changed in shape over the past century, reflecting its dynamic origin and activity, also exemplified by the sequential surges identified since 1899. There was a pre-surge build-up event of Flatbreen causing a surge and subsequent lowering of the Hornbreen-Hambergbreen terminus by the 1960s. After, the entire surface lowered, albeit with a delay in the Hornbreen terminal zone. Since the year 2000, Hornbreen terminus has retreated at an average rate of 106 m a-1; ~50% faster than that of Hambergbreen. If the retreat continues at the 2000 – 2015 average rate, the ice bridge between Hornsund and Hambergbukta will be broken sometime between 2055 and 2065 and the Hornsund strait will separate Sorkapp Land from the Spitsbergen island. The processes and events described in this study, particularly the effects of the glacier surge, may provide a model for changes likely to occur in other coastal glaciated regions experiencing rapid change.

A numerical algorithm for the assessment of the conjecture of a subglacial lake tested at Amundsenisen, Svalbard

The melting of glaciers coming with climate change threatens the heritage of the last glaciation of Europe likely contained in subglacial lakes in Greenland and Svalbard. This aspect urges specialists to focus their studies (theoretical, numerical, and on-field) on such fascinating objects. Along this line, we have approached the validation of the conjecture of the existence of a subglacial lake beneath the Amundsenisen Plateau at South-Spitzbergen, Svalbard, where ground penetrating radar measurements have revealed several flat signal spots, the sign of the presence of a body of water. The whole investigation aspects, mathematical modeling and numerical simulation procedure, and the numerical results are presented through a trilogy of papers of which the present one is the last. The time-dependent mathematical model in the background of the numerical algorithm includes the description of dynamics and thermodynamics of the icefield and of the subglacial lake, with heat exchange and liquid/solid phase-change mechanisms at the interface. Critical modeling choices and confidence in the algorithm are granted by the numerical results of the sensitivity analysis versus the contribution of ice water content, of firn and snow layers at top of the icefield and versus the approximation of ice sliding on bedrock. The two previous papers deal with these issues, show successful comparison with local measured quantities, and demonstrate numerically the likelihood of the subglacial lake. In this work, we aim at providing the studied case and the numerical algorithm with a possible paradigmatic value. At this aim, we introduce on-field measurement data related to the physical characteristics of the Amundsenisen Plateau that justify the adoption of significant modeling simplifications, here, focussed from physical viewpoint. Furthermore, we present the numerical algorithm and discuss several representative results from the numerical test to point out the type of results coming from the procedure. Such results might, eventually, provide a support to the decision to undertake drilling operations for tracing the subglacial water bio-chemicals generally present within the accreted ice above the presumed ice/water front.