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.

Abstract. The success of geological carbon storage depends on the assurance of permanent containment for injected carbon dioxide (CO2) in the storage formation at depth. One of the critical elements of the safekeeping of CO2 is the sealing capacity of the caprock overlying the storage formation despite faults and/or fractures, which may occur in it. In this work, we present an ongoing injection experiment performed in a fault hosted in clay at the Mont Terri underground rock laboratory (NW Switzerland). The experiment aims to improve our understanding of the main physical and chemical mechanisms controlling (i) the migration of CO2 through a fault damage zone, (ii) the interaction of the CO2 with the neighboring intact rock, and (iii) the impact of the injection on the transmissivity in the fault. To this end, we inject CO2-saturated saline water in the top of a 3 m thick fault in the Opalinus Clay, a clay formation that is a good analog of common caprock for CO2 storage at depth. The mobility of the CO2 within the fault is studied at the decameter scale by using a comprehensive monitoring system. Our experiment aims to close the knowledge gap between laboratory and reservoir scales. Therefore, an important aspect of the experiment is the decameter scale and the prolonged duration of observations over many months. We collect observations and data from a wide range of monitoring systems, such as a seismic network, pressure temperature and electrical conductivity sensors, fiber optics, extensometers, and an in situ mass spectrometer for dissolved gas monitoring. The observations are complemented by laboratory data on collected fluids and rock samples. Here we show the details of the experimental concept and installed instrumentation, as well as the first results of the preliminary characterization. An analysis of borehole logging allows for identifying potential hydraulic transmissive structures within the fault zone. A preliminary analysis of the injection tests helped estimate the transmissivity of such structures within the fault zone and the pressure required to mechanically open such features. The preliminary tests did not record any induced microseismic events. Active seismic tomography enabled sharp imaging the fault zone.

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. Englacial conduits act as water pathways to feed surface meltwater into the subglacial drainage system. A change of meltwater into the subglacial drainage system can alter the glacier's dynamics. Between 2012 and 2019, repeated 25 MHz ground-penetrating radar (GPR) surveys were carried out over an active englacial conduit network within the ablation area of the temperate Rhonegletscher, Switzerland. In 2012, 2016, and 2017 GPR measurements were carried out only once a year, and an englacial conduit was detected in 2017. In 2018 and 2019 the repetition survey rate was increased to monitor seasonal variations in the detected englacial conduit. The resulting GPR data were processed using an impedance inversion workflow to compute GPR reflection coefficients and layer impedances, which are indicative of the conduit's infill material. The spatial and temporal evolution of the reflection coefficients also provided insights into the morphology of the Rhonegletscher's englacial conduit network. During the summer melt seasons, we observed an active, water-filled, sediment-transporting englacial conduit network that yielded large negative GPR reflection coefficients (<-0.2). The GPR surveys conducted during the summer provided evidence that the englacial conduit was 15–20 m±6 m wide, ∼0.4m±0.35m thick, ∼250m±6m long with a shallow inclination (2∘), and having a sinusoidal shape from the GPR data. We speculate that extensional hydraulic fracturing is responsible for the formation of the conduit as a result of the conduit network geometry observed and from borehole observations. Synthetic GPR waveform modelling using a thin water-filled conduit showed that a conduit thickness larger than 0.4 m (0.3× minimum wavelength) thick can be correctly identified using 25 MHz GPR data. During the winter periods, the englacial conduit no longer transports water and either physically closed or became very thin (<0.1 m), thereby producing small negative reflection coefficients that are caused by either sediments lying within the closed conduit or water within the very thin conduit. Furthermore, the englacial conduit reactivated during the following melt season at an identical position as in the previous year.

Traditionally, helicopter-borne ground-penetrating radar (GPR) systems are operated with a single pair of bistatic dipole antennas to measure the thickness of glaciers. We demonstrate numerically that the directivity of the radiation pattern of single airborne dipoles do not correspond to an ideal fullspace solution if the antennas are employed at typical flight heights. These directionality effects can degrade the quality of the subsurface images significantly, when the GPR antennas are orientated unfavorably. Since an adjustment of the antenna orientation is impractical during flight, we have developed a novel dual-polarization helicopter-borne GPR system consisting of two orthogonal pairs of commercial antennas in broadside configuration. To overcome the image quality deficits of the individual channels, we apply a pseudo-scalar approach in which we combine the data of both polarizations. Results of helicopter-borne GPR surveys on two alpine glaciers in Switzerland reveal more coherent bedrock reflections in the summed data compared with single dipole pair profiles. Generally, the dual-polarization setup is more suitable than a single antenna systems, because it is more versatile and less prone to directional effects caused by the placement of the dipole antennas in relation to undulating subsurface reflectors.

Abstract. Seismic investigations of geothermal reservoirs over the last 20 years have sought to interpret the resulting tomograms and reflection images in terms of the degree of reservoir fracturing and fluid content. Since the former provides the pathways and the latter acts as the medium for transporting geothermal energy, such information is needed to evaluate the quality of the reservoir. In conventional rock physics-based interpretations, this hydro-mechanical information is approximated from seismic velocities computed at the low-frequency (field-based) and high-frequency (labbased) limits. In this paper, we demonstrate how seismic properties of fluid-filled, fractured reservoirs can be modeled over the full frequency spectrum using a numerical simulation technique which has become popular in recent years. This technique is based on Biot's theory of poroelasticity and enables the modeling of the seismic velocity dispersion and the frequency dependent seismic attenuation due to waveinduced fluid flow. These properties are sensitive to key parameters such as the hydraulic permeability of fractures as well as the compressibility and viscosity of the pore fluids. Applying the poroelastic modeling technique to the specific case of a magmatic geothermal system under stress due to the weight of the overlying rocks requires careful parameterization of the model. This includes consideration of the diversity of rock types occurring in the magmatic system and examination of the confining-pressure dependency of each input parameter. After the evaluation of all input parameters, we use our modeling technique to determine the seismic attenuation factors and phase velocities of a rock containing a complex interconnected fracture network, whose geometry is based on a fractured geothermal reservoir in Iceland. Our results indicate that in a magmatic geothermal reservoir the overall seismic velocity structure mainly reflects the lithological heterogeneity of the system, whereas indicators for reservoir permeability and fluid content are deducible from the magnitude of seismic attenuation and the critical frequency at which the peak of attenuation and maximum velocity dispersion occur. The study demonstrates how numerical modeling provides a valuable tool to overcome interpretation ambiguity and to gain a better understanding of the hydrology of geothermal systems, which are embedded in a highly heterogeneous host medium.

No abstract