Accurate knowledge of the ice thickness distribution and glacier bed topography is essential for predicting dynamic glacier changes and the future developments of downstream hydrology, which are impacting the energy sector, tourism industry and natural hazard management. Using AIR-ETH, a new helicopter-borne ground-penetrating radar (GPR) platform, we measured the ice thickness of all large and most medium-sized glaciers in the Swiss Alps during the years 2016–20. Most of these had either never or only partially been surveyed before. With this new dataset, 251 glaciers – making up 81% of the glacierized area – are now covered by GPR surveys. For obtaining a comprehensive estimate of the overall glacier ice volume, ice thickness distribution and glacier bed topography, we combined this large amount of data with two independent modeling algorithms. This resulted in new maps of the glacier bed topography with unprecedented accuracy. The total glacier volume in the Swiss Alps was determined to be 58.7 ± 2.5 km3 in the year 2016. By projecting these results based on mass-balance data, we estimated a total ice volume of 52.9 ± 2.7 km3 for the year 2020. Data and modeling results are accessible in the form of the SwissGlacierThickness-R2020 data package.

Knowing the ice thickness distribution of a glacier is of fundamental importance for a number of applications, ranging from the planning of glaciological fieldwork to the assessments of future sea-level change. Across spatial scales, however, this knowledge is limited by the paucity and discrete character of available thickness observations. To obtain a spatially coherent distribution of the glacier ice thickness, interpolation or numerical models have to be used. Whilst the first phase of the Ice Thickness Models Intercomparison eXperiment (ITMIX) focused on approaches that estimate such spatial information from characteristics of the glacier surface alone, ITMIX2 sought insights for the capability of the models to extract information from a limited number of thickness observations. The analyses were designed around 23 test cases comprising both real-world and synthetic glaciers, with each test case comprising a set of 16 different experiments mimicking possible scenarios of data availability. A total of 13 models participated in the experiments. The results show that the inter-model variability in the calculated local thickness is high, and that for unmeasured locations, deviations of 16% of the mean glacier thickness are typical (median estimate, three-quarters of the deviations within 37% of the mean glacier thickness). This notwithstanding, limited sets of ice thickness observations are shown to be effective in constraining the mean glacier thickness, demonstrating the value of even partial surveys. Whilst the results are only weakly affected by the spatial distribution of the observations, surveys that preferentially sample the lowest glacier elevations are found to cause a systematic underestimation of the thickness in several models. Conversely, a preferential sampling of the thickest glacier parts proves effective in reducing the deviations. The response to the availability of ice thickness observations is characteristic to each approach and varies across models. On average across models, the deviation between modeled and observed thickness increase by 8.5% of the mean ice thickness every time the distance to the closest observation increases by a factor of 10. No single best model emerges from the analyses, confirming the added value of using model ensembles.

Englacial hydrology plays an important role in routing surface water to the glacier's bed and it consequently affects the glacier's dynamics. However, it is often difficult to observe englacial conduit conditions on temperate glaciers because of their short-lived nature. We acquired repeated active surface seismic data over the Rhone Glacier, Switzerland to monitor and characterise englacial conduit conditions. Amplitude-versus-angle analysis suggested that the englacial conduit is water filled and between 0.5 and 4 m thick. A grid of GPR profiles, acquired during the 2018 melt season, showed the englacial conduit network persisting and covering ~ 14,000 m2. In late summer 2018, several boreholes were drilled into the conduit network. We observed generally stable water pressure, but there were also short sudden increases. A borehole camera provided images of a fast flowing englacial stream transporting sediment through the conduit. From these observations, we infer that the englacial conduit network is fed by surface meltwater and morainal streams. The surface and morainal streams merge together, enter the glacier subglacially and flow through subglacial channels along the flank. These subglacial channels flow into highly efficient englacial conduits traversing the up-glacier section of the overdeepening before connecting with the subglacial drainage system.

No abstract

Abstract. Ground-penetrating radar (GPR) and seismic imaging have proven to be important tools for the characterization of rock volumes. Both methods provide information about the physical rock mass properties and geological structures away from boreholes or tunnel walls. Here, we present the results from a geophysical characterization campaign that was conducted as part of a decametre-scale hydraulic stimulation experiment in the crystalline rock volume of the Grimsel Test Site (central Switzerland). For this characterization experiment, we used tunnel-based GPR reflection imaging as well as seismic travel-time tomography to investigate the volumes between several tunnels and boreholes. The interpretation of the GPR data with respect to geological structures is based on the unmigrated and migrated images. For the tomographic analysis of the seismic first-arrival travel-time data, we inverted for an anisotropic velocity model described by the Thomsen parameters v0, ϵ and δ to account for the rock mass foliation. Subsequently, the GPR and seismic images were interpreted in combination with the geological model of the test volume and the known in situ stress states. We found that the ductile shear zones are clearly imaged by GPR and show an increase in seismic anisotropy due to a stronger foliation, while they are not visible in the p-wave (v0) velocity model. Regions of decreased seismic p-wave velocity, however, correlate with regions of high fracture density. For geophysical characterization of potential deep geothermal reservoirs, our results imply that wireline-compatible borehole GPR should be considered for shear zone characterization, and that seismic anisotropy and velocity information are desirable to acquire in order to gain information about ductile shear zones and fracture density, respectively.