The rate of fracture-induced ice instability is an important factor contributing to uncertainties in sea level projections used for global flood mitigation planning. While the occurrence of ice fracturing at critical stress thresholds is well-documented, the detailed mechanisms controlling fracture timing, rate, and orientation are not fully understood. This gap is particularly evident in differences in fracture behaviour across varying ice types, such as meteoric ice and ice mélange. Observations on the Brunt Ice Shelf reveal a unique behaviour, where rifts deviate from the pathway predicted by the principal stresses to avoid thick blocks of meteoric ice. Their growth rate is significantly reduced when required to cross through these blocks. This stands in contrast to observations on other ice shelves, such as Larsen C, where rift propagation is slower in marine ice bands.

Here we use co-located geophysical methods, seismic and ground-penetrating radar (GPR), to assess the fracture pattern and dynamics and the relationship to ice properties at the leading edge of two active rifts, Halloween Crack and Chasm 2, on the Brunt Ice Shelf.

By determining the depth of seismic events using P to Rayleigh wave amplitude ratios, we estimate a theoretical maximum dry crevasse depth—the depth at which fracturing can occur without the presence of englacial water. Additionally, GPR data are used to precisely locate rift terminations and identify refrozen layers associated with seawater intrusion into the firn layer. Combining these data, we provide new insight into the mechanisms controlling fracture propagation within the Brunt Ice shelf. The synthesis of observations from Chasm 2 and Halloween Crack contributes to a comprehensive understanding of fracture mechanics, enhancing our knowledge of regional-scale ice dynamics.

How to cite: Pearce, E., Marsh, O., Brisbourne, A., and Hudson, T.:  Assessing the rate of ice fracture using co-located geophysical surveys on the Brunt Ice Shelf, Antarctica, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13663, https://doi.org/10.5194/egusphere-egu24-13663, 2024.

Reconstructing past ice sheets is important for understanding the response of modern ice sheets to changes in climate. The evolution of the Weddell Sea Sector’s grounding line since the last glacial maximum (LGM) to its present position remains ambiguous; previous authors have proposed hypotheses both of monotonic grounding line retreat and of rapid grounding line retreat followed by readvance. However, distinguishing these scenarios with current observations remains difficult. To explore these scenarios, we report seismic measurements of basal properties at KIR, an ice rise in the Weddell Sea Sector, West Antarctica. A three-component seismic survey enabled detection of the compressional (P) wave reflection and the converted (PS) wave reflection (an incident P wave converted to a shear wave at the base-ice reflector) from the base of KIR. Amplitude-vs-angle (AVA) analysis aims to constrain the physical properties (namely density, seismic velocity, by measuring the variation of reflectivity with incidence angle at the reflector. By jointly inverting the AVA responses of the PP wave reflection and the PS reflection, we increase the confidence in the interpretation of the base-ice properties. 

Analysis of PP and PS AVA responses at KIR indicates that the reflection arises from a material with a P wave velocity of 4.03 ± 0.05 km/s, an S wave velocity of 2.16 ± 0.06 km/s and a density of 1.44 ± 0.06 g/cm³; these properties are consistent with a reflection from a layer of entrained basal debris, with 20-30% debris by volume. The observed properties are not indicative of interference at a thin layer, as observed beneath glaciers elsewhere. The absence of deeper subglacial reflections indicates a poorly-defined boundary between this basal debris layer and the underlying subglacial material, which we therefore propose consists of frozen sediments . If this interpretation is correct, the presence of a debris layer overlying basal frozen sediment indicates a potential retreat/readvance scenario for KIR. A possible scenario is a previous episode of flow during which KIR may have been weakly grounded as an ice rumple, followed by grounding on the lee side of the bathymetric high and subsequent freezing of subglacial sediments. However, the origin of such a homogeneous and debris-rich layer remains unclear. The indication of a reflection from a basal debris layer raises questions about whether conventionally interpreted basal reflections can truly be considered as such, and whether these interpretations may mask the true nature of the underlying subglacial material. This ambiguity may be most effectively reconciled by borehole sampling.

How to cite: Agnew, R., Brisbourne, A., Booth, A., and Clark, R.: Basal debris at an Antarctic ice rise revealed by seismic amplitude-vs-angle analysis. , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-1383, https://doi.org/10.5194/egusphere-egu24-1383, 2024.

Airborne survey is one of the most important observational techniques in environmental science. This is especially true in polar settings where access is challenging and observational requirements, such as ice sounding radar, in situ study of turbulent atmospheric processes, cloud cover, or requirements for high resolution potential field data, limit use of satellite data. Although critical, airborne survey using traditional platforms, such as the versatile twin otter aircraft operated by the British Antarctic Survey (BAS), come with a relatively high logistical, financial, and environmental (CO₂) footprint. Larger UAV’s offer an alternative, but as yet un-realised, lower impact platform to deliver the same, if not more scientific data.

Through the Innovate UK SWARM project BAS is collaborating with Windracers to trial their large (10 m wing span) Ultra UAV as a platform for environmental science. Making use of the large (700 L/max 100 kg), easily accessible payload bay and a series of interchangeable payload floors this trial will be carried out in February/March 2023. The science payloads will include: Atmospheric (turbulence probe), environmental (hyperspectral and visual cameras), cryosphere (600-900 MHz accumulation radar), and potential field geophysics (gravity/magnetic sensors). The missions, between 10 and 330 km long, will be flown beyond visual line of sight (BVLOS) of the operator using the Distributed Avionics autopilot, including take-off and landing, which will be overseen by an in-field safety pilot.

Here we present the first results of this trial, including our experience integrating BVLOS UAV operations with traditional aircraft in an Antarctic context and initial results and lessons learned from the four trailed instrument suites. Our demonstration will be an important milestone in the transition to widespread use of larger UAVs for environmental science. We will discuss how the reduced environmental and logistical impact can open up new opportunities in Antarctic and beyond.

How to cite: Jordan, T. and Robinson, C.: Demonstrating a large UAV for Antarctic environmental science, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5435, https://doi.org/10.5194/egusphere-egu24-5435, 2024.

Estimations of bedrock topography below glaciers and ice thickness are vital for quantifying freshwater availability for surrounding populations and understanding the contribution of melting glaciers to sea-level rise in the context of global warming. While active seismology is commonly used for ice thickness estimation, the utilization of passive methods remains relatively rare. Passive seismology solutions offer cost-effectiveness, non-invasiveness and continuous monitoring capabilities that present valuable benefits in glaciological research.

Over the past two decades, numerous seismic stations have been deployed on glaciers worldwide for various purposes. Through passive seismology approaches, these seismic stations could show their potential as new sources of ice thickness measurements and feed the related database. For this purpose, we analyzed data of 3-components seismic sensors from different deployments as well as data from open access databases, such as IRIS, employing the Horizontal-to-Vertical Spectral Ratio (HVSR) technique. HVSR has been predominantly used in microzonation studies to determine site effects and the thickness of sediments in sedimentary basins.  Even though the use of this technique in glacial seismology is quite new, HVSR has been already utilized to estimate in-situ ice thickness, to retrieve the basal properties or to detect cavities under the ice.

Our primary objective is to demonstrate the potential of the HVSR technique to retrieve in-situ ice thickness on different glaciers. Subsequently, we will compare the HVSR results with different data sources including models, ground-penetrating radar or active seismology. By performing this comparison we evaluate the limitations of the HVSR method in an icy environment. We investigate these limitations by studying the effect of other natural agents such as wind on the H/V amplitude and fundamental frequency retrieved from the HVSR curves. Having a global understanding of these influences will eventually help deciphering variations in continuous H/V monitoring.

How to cite: Govoorts, J., Van Noten, K., Caudron, C., Einarsson, B., Lecocq, T., Nowé, S., Pálsson, F., Pätzel, J., and Zekollari, H.: Analysis of H/V spectral ratio curves from passive seismic data acquired on glaciers worldwide, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10814, https://doi.org/10.5194/egusphere-egu24-10814, 2024.

How to cite: Smith, E. C., Agnew, R., Booth, A. D., Christoffersen, P., Dawson, E. J., Gonzalez, L., Karplus, M., May, D. F., Nakata, N., Pretorius, A., Summers, P., Tulaczyk, S., Vietch, S., Walter, J., and Young, T. J.: Icequakes beneath Thwaites Glacier eastern shear margin , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11300, https://doi.org/10.5194/egusphere-egu24-11300, 2024.

Subglacial cavities may trap a considerable quantity of liquid water, causing devastating outburst floods in densely populated mountain areas. Both active and passive geophysical methods are employed for the glacier-bedrock interface and intra-glacial characterization, including Ground Penetrating Radar (GPR), refraction seismic, borehole measurements, and surface nuclear magnetic resonance (SNMR). 

Ambient seismic noise can be collected by light and dense arrays at a relatively moderate cost, and allows to access some mechanical properties of the glacier, including the detection and localization of ice cavities and, possibly, basal detachment, taking advantage of spectral anomalies in the horizontal-to-vertical-spectral ratio (HVSR) and in the Vertical-to-Horizontal spectral ratio (VHSR). Specifically, a peak in the VHSR indicates a low impedance volume beneath the surface [1,2]. As a simple picture, we can refer to the “bridge” vibrating mode, where the vertical displacement in the middle of the bridge largely dominates other components of the movement.  Antunes et al. [2] furthermore noticed that the VHSR gives information about seismic energy anomalies generated by fluids in reservoirs since the wavefield is polarized mainly in the vertical direction.
In this work, we apply the HVSR and VHSR techniques to locate a subglacial water-filled cavity in the Tête Rousse glacier (Mont Blanc area, French Alps), using 15 days of data collected in may, 2022 [3]. The results also confirm the general basal conditions of the glacier suggested by other methods, locating temperate areas of the glacier where basal detachments are possible.

We evaluate the optimal seismic noise record duration to obtain a reliable and stable mapping of the VHSR over the glacier to properly locate the main cavity (or secondary cavities). In our case, results suggest that 6 days of record are enough to detect and locate a cavity

[1] Saenger, E-H. et al: A passive seismic survey over a gas field: Analysis of low-frequency anomalies, Geophysics, 74 (2), O29–O40 (2009).

[2] Antunes V. et al: Insights into the dynamics of the Nirano Mud Volcano through seismic characterization of drumbeat signals and V/H analysis. Journal of Volcanology and Geothermal Research, 431 (2022).

[3] A. Guillemot, N. Bontemps, E. Larose, D. Teodor, S. Faller, L. Baillet, S. Garambois, E. Thibert, O. Gagliardini, C. Vincent: Investigating Subglacial Water-filled Cavities by Spectral Analysis of Ambient Seismic Noise : Results on the Polythermal Tête-Rousse Glacier (Mont Blanc, France), Geophys. Res. Lett. accepted (2024). DOI:10.1029/2023GL105038

How to cite: Larose, E., Bontemps, N., Guillemot, A., and Baillet, L.: Locating subglacial cavities and investigating basal conditions on glaciers with ambient seismic noise: toward acquisition optimization., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12583, https://doi.org/10.5194/egusphere-egu24-12583, 2024.

The Totten Glacier is a fast moving glacier that serves as a major outlet of the East Antarctic Ice Sheet. During December 2018 and January 2019, we deployed a 12 station broadband seismic array near the grounding zone of the Totten Glacier. We observed a significant number ( > 10,000) of repeating basal stick-slip icequakes across the region. Much of this seismic activity was dominated by higher frequency events (20-75 Hz) similar in size and temporal character (“bursty”) to those found in previous studies, such as those on the Rutford Ice Stream and Greenland Ice Sheet. Additionally, we observe a large number of repeating events dominated by lower frequencies (< 10 Hz) that have larger magnitudes and longer inter-event time than the high-frequency seismic activity. We will provide an overview into both the temporal and spatial variability of this seismic activity and discuss implications for fast flow in the region.

How to cite: Winberry, P.: Repeating Glacier Seismicity Near the Totten Glacier Grounding Zone., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12865, https://doi.org/10.5194/egusphere-egu24-12865, 2024.

The state and evolution of subglacial channels strongly impact glacier motion and as a result the mass balance of flowing ice bodies. Yet, the subglacial environment is difficult to access and thus often poorly constrained over significant temporal and spatial scales. This limits our understanding of complex subglacial hydraulic processes and consequently ice dynamics.

Seismology can help overcome these observational constraints, providing new insights into fundamental processes in the cryosphere, such as frictional sliding and subglacial water flow. However, different seismogenic processes of the cryosphere often overlap in both time and space. Differentiating between them and interpreting associated seismic signals require appropriate methodological and instrumental approaches.

Here, we investigate subglacial channel dynamics at the Rhone glacier (Switzerland) over one month in the summer of 2020, focusing on periods coinciding with glacier sliding episodes. To this end, we leverage the sensitivity of near-bed borehole geophones combined with seismic interferometry and beamforming techniques.

We show that the hydraulic tremor, generated by turbulent water flow and resulting pressure variations acting against the subglacial channel bed and walls, acts as a dominant, stable, and coherent noise source. Beamforming analysis reveals the directional stability of the hydraulic tremor and points toward the junction of two subglacial hydraulic channels from which stick-slip asperities originate. The analysis also reveals instances of sudden hydraulic tremor quieting, in agreement with previous observations before and after seismogenic sliding episodes. We explain this quieting as sudden changes in frictional conditions within the subglacial channel corresponding to a rapid transition between a fully and partially filled channel. We discuss channel properties (geometry and bed conditions) that are needed to satisfy the physical conditions for the frictional quieting mechanism. Our analysis offers new insights into the complex mechanical interactions between ice, water, and bed properties and the hydraulic control of glacier sliding.

How to cite: Chmiel, M., Caldera, N., Walter, F., Olivier, G., Farinotti, D., Guadagnini, A., Gräff, D., Köpfli, M., and Gimbert, F.: Quieting of hydraulic tremor: sudden changes in frictional conditions in subglacial channels, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16320, https://doi.org/10.5194/egusphere-egu24-16320, 2024.

An ice core is being drilled near Little Dome C, a small promontory about 30km downstream from the summit of Dome C, to extract a continuous record of climate over the last 1.5 million years. Present and past ice flow conditions are important to interpret the ice core because the surface velocity at the drilling site is about 40 mm/yr and the oldest ice in the record was deposited in the surface about 10km upstream of the drilling site. Here we explore newly acquired and existing geophysical data to describe present ice flow and investigate signs of past changes. We present new GNSS data that describes the subtle but complex local surface velocity, and ApRES radar data that provides englacial strain-rates along the flow path from the summit of Dome C and bulk englacial crystal orientation fabric. Ice currently flows from Dome C summit along the ridge to Little Dome C, even though a subtle uphill slope, but basal conditions are variable along the path due to the strong basal topography. Of special interest is an ice unit in contact with the bedrock with variable thickness up to about 300m that is vertically stagnant and produce a strong radar reflection.  This basal unit is not present in an area of strong melting about 5km upstream from the drilling site. The crystal orientation fabric reflects the ice flow horizontal extension along the path and changes with depth on ice flow properties following climatic transitions and, more intriguing, indicate a possible change in ice flow extension at the beginning of the Holocene. We aim to facilitate detailed ice flow models to better interpret the ice core data.  

How to cite: Mulvaney, R., Martín, C., Ritz, C., Vittuari, L., Frezzotti, M., and Eisen, O.: Glaciological characterization of Little Dome C: Influence of ice flow on the future Beyond Epica – Oldest Ice Core drilling project, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17767, https://doi.org/10.5194/egusphere-egu24-17767, 2024.

Ground Penetrating Radar (GPR) is a major tool to investigate, map and monitor polar ice sheets and alpine glaciers. Alpine glaciers are often composed of temperate ice, which has significantly different backscatter properties from cold ice. Radar attenuation is much stronger in temperate ice than in cold ice, because the radar signal encounters strong scattering in temperate ice. 
A major candidate for this scattering is the presence of liquid water inclusions, which are much smaller than the radar wavelength. The large contrast between water and ice dielectric permittivity would explain the diffuse radar scattering in temperate ice. Indeed, recent numerical modelling of the radar signal in temperate ice confirmed the contribution of liquid water inclusions on the scattering of the radar signal (Ogier, 2023). 
Here, we investigate if the strong scattering caused by liquid water inclusions, which is usually treated as noise, can be in fact exploited to unravel dynamic processes inside the glacier. This strong scattering results in large radar phase variations in space, which remain constant over short timescales (hours - days), during which the glacial water content remains constant. During that timescale, however, the mountain glacier might experience sudden internal deformation due to intermittent sliding at the glacier base, also called glacier stick-slip.  This deformation might be resolved using difference imaging and the spatio-temporal properties of the radar phase.
We numerically model radar wave propagation throughout temperate ice (i.e. with the presence of liquid water inclusions) before and after an idealized glacier deformation and show, that through phase difference imaging the internal movement of the sub-wavelength scatterers can be mapped. 
Finally, we discuss how this novel type of monitoring could be applied in the field, which is planned for spring 2024.

Ogier, Christophe, Dirk-Jan van Manen, Hansruedi Maurer, Ludovic Räss, Marian Hertrich, Andreas Bauder, and Daniel Farinotti. 2023. “Ground Penetrating Radar in Temperate Ice: Englacial Water Inclusions as Limiting Factor for Data Interpretation.” Journal of Glaciology, September, 1–12. https://doi.org/10.1017/jog.2023.68.

How to cite: Aichele, J., Ogier, C., and Anhorn, B.: Stick-slip imaging through the GPR phase: Turning  temperate ice 'noise' into signal, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18639, https://doi.org/10.5194/egusphere-egu24-18639, 2024.

Seismic events (icequakes) associated with floating ice sheets on lakes are a frequently observed phenomenon. We find at our study site on the frozen Lake St. Moritz in the Swiss Alps typically a clear diurnal pattern with hundreds to thousands of icequake signals per hour during night time, while the rate of observed events during daytime is about two orders of magnitude smaller. The seismicity rate shows a significant correlation with temperature changes. It is therefore assumed that the generation of the ice quakes is related to melting and freezing processes as well as the extension and contraction of the ice. Potentially the seismicity rate is also moderated by loading and unloading due to human activities on the ice and/or lake level changes.

These ice quakes generate seismic waves that propagate through the thin ice sheet as plate waves modulated by the air and water half-spaces above and below the ice (quasi-guided waves). One member of this wave-type family, the quasi-Scholte waves, are characterised by distinct dispersion that can be observed with seismic sensors on the ice. Furthermore, the seismic waves traveling through the ice couple into the air leading to audible seismo-acoustic signals. One particularity of the ice-air coupling is a so-called coincidence phenomenon. The particular velocity-frequency combination where the seismic wavelength in the ice matches the apparent acoustic wavelength in the air leads to a resonance phenomenon. Observation of the related coincidence frequency allows us, for example, to infer on the ice thickness from the acoustic observations with a low cost microphone above the ice only. Recording the acoustic signals with small microphone arrays enables additionally, for example, locating the source of the seismo-acoustic signal.

Combined observations of the seismic and acoustic signals provide new insights into the seismicity of lake ice which has rarely been studied in the past. The seismo-acoustic signals have the potential to provide information about the ice properties such as thickness and ice quality as well as waxing and waning processes of ice sheets. These observations are relevant for safe operations on the ice but also to complement other remote-sensing observations with autonomous in situ seismo-acoustic measurements for climate studies.

How to cite: Schmelzbach, C., Wetter, C., Stähler, S., Clinton, J., Lyu, Z., Mesimeri, M., and Massin, F.: Lake ice seismicity: seismic and acoustic observations, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17421, https://doi.org/10.5194/egusphere-egu24-17421, 2024.

Studying Arctic sea ice is essential as it plays an important role in climate regulation, influencing weather
patterns, as well as impacting the local ecosystems and living conditions of people. Among different
methods for collecting data and studying sea ice, seismology has proved to be an efficient way to extract
the ice properties, from which the mechanical behavior of sea ice can be explored. Seismic data recorded on
sea ice, using 3-component geophones, are used as a starting point to derive useful information regarding
the ice. To devise an efficient inversion method for deriving ice properties, an effective tool would be
necessary to generate synthetic data in a way that encompasses the physics of sea ice. While there are
approximate solutions to wave propagation problem in a floating ice layer based on plate theory, which is
based on the assumptions of homogeneity of the ice layer and valid at low values of frequency×wavelength,
numerical counterparts such as wavenumber integration method and finite element method have been
also used to to create synthetic waveforms. The numerical methods have shown the limitations of these
approximate solutions in modeling wave propagation; nonetheless, the effects of these limitations on the
estimations of the location of icequakes and thickness of ice need to be investigated.

In this study, these limitations are explored. To do this, two possible scenarios that can happen in
practice are taken into account: (1) when there is high-frequency content in the source generating the
seismic data, and (2) when the physical model includes a snow layer overlying the ice layer. First, we
will show the limitations of the approximate solutions for these two cases by comparing the waveforms,
derived from these approximate solutions, with those of a numerical method at a given distance from
the source. The numerical used here is spectral element method. Then, the effects of these limitations
on the estimations of icequake location and ice thickness are explored in an inversion process, in which
synthetic data are created using the approximate solutions. Results indicate that when there are high-
frequency content in the data and a snow layer on top of the ice, the use of the approximate solutions
to generate synthetic data introduces bias in the estimation of ice thickness and source-receiver distance
in the inversion process. This bias is in the form of underestimations, smaller ice thicknesses and smaller
source-receiver distances. Furthermore, to tackle the biases associated with the inversion method based on
the approximate solutions, a novel strategy is adopted, where a database of simulations using the proposed
numerical method is built for various models of ice and snow. Here the inversion comprises of searching
in the database to find the best ice thickness and source-receiver distance for each icequake. In addition,
the database-based inversion reduces the computational cost. Thanks to this inversion strategy, and
using real data recorded on sea ice, the ice thicknesses along different source-receiver paths are estimated
efficiently, from which a 3D map of ice thickness is constructed.

How to cite: Zandi, H., Moreau, L., Métivier, L., and Brossier, R.: Sensitivity Study for Seismic Waves Guided in an Ice Pack: Influence of the Frequency Content and Snow Layer Thickness Covering the Ice, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8245, https://doi.org/10.5194/egusphere-egu24-8245, 2024.

The observation of near-ground cosmogenic neutrons enables the monitoring of various water storage variations at the land surface at the field scale including soil moisture and the water content of snow layers. The parabolic neutron-count versus soil moisture function is quite uniform among different locations, and soil types and requires typically a one-time-only in situ reference observation. For the detection of snowpack water equivalent (SWE) variations by cosmic-ray neutron sensing such a uniform approach has so far not been developed. Therefore, the establishment of new cosmic-ray snow monitoring sites requires substantial in situ measurements for obtaining the local relationship of SWE amounts and neutron count rates. Observations suggest that the relationship is quite uniform for grass-vegetated locations which is different to what is found for stony ground.

Within the framework of the research unit Cosmic Sense, we generated extensive in situ measurements of snow water equivalent and cosmogenic neutron count rates at various sites with differing elevations in the German and Austrian Alps. From these data, we investigate commonalities among the site conditions and if the varying patterns of the relationships can be reasonably explained by physical reasons and therefore be modeled with a unified approach.

How to cite: Fersch, B., Součková, M., Schattan, P., Krebs, N., Weimar, J., Jahn, C., Grosse, P. M., and Schrön, M.: Towards a unified description of the count rate – snow water equivalent relationship in cosmic-ray neutron sensing, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11386, https://doi.org/10.5194/egusphere-egu24-11386, 2024.

Above-ground cosmic ray neutron sensing (CRNS) is an emerging technique for the investigation of dynamics in soil moisture, snow water equivalent (SWE), and vegetation at a spatial scale of several hectares. The measurement principle is based on the moderation of natural secondary cosmogenic neutrons by hydrogen atoms. On the earth surface hydrogen atoms are mainly bound in water molecules. However, at complex research sites the signal distinction between various water sources remains challenging. Especially in alpine terrain and at elevated topography, hydrological features are linked in an intricate patchwork, hampering signal discrimination. Satellite observations offer valuable complementary surface information and are commonly provided at a spatial resolution that meets the integrated footprint area of the CRNS detector. In this study we investigate if the interpretation of the CRNS signal can be enhanced by the use of remote sensing products. We compare three readily available fractional snow cover (FSC) products based on Sentinel (1 and 2) and MODIS and one reference FSC Sentinel-2 scene-based machine learning product at the approximate footprint resolution of CRNS, comprising a circular area of 250 m radius. The performance of all four products is assessed at five CRNS sites in the Austrian and Italian Alps that represent a variety of environmental properties, ranging from flat to steep topography, from low to high elevation and from sparse to abundant vegetation cover. At three sites, the presence and absence of snow can be validated by local snow height measurements. The analysis shows that remote sensing snow cover information can be extracted on around 80% of the analyzed days, demonstrating the use of FSC products for the estimation of snow cover duration and timing. Comparing the four products shows overall agreements and allows to deduce product-specific thresholds for the distinction of snow-covered and snow-free situations. Further, pairing remote FSC observations with neutron count measurements provides a first indication on the complexity of local hydrogen pool dynamics and consequent requirements on the calibration routine for ambient water monitoring with CRNS. We conclude that satellite-based FSC products can be used to fortify the choice of CRNS observation location and period prior to the detector installation and for a robust and viable first-order assessment of expected CRNS site conditions. Remote sensing FSC products and CRNS measurements hold complementary data that can mutually benefit snow observations and should be explored further in the future.

How to cite: Krebs, N., Schattan, P., Premier, V., Mejia-Aguilar, A., and Rutzinger, M.: Testing four Sentinel (1 and 2) and MODIS Fractional Snow Cover products for the evaluation of five Alpine Cosmic Ray Neutron Sensing sites, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13683, https://doi.org/10.5194/egusphere-egu24-13683, 2024.

Cosmic-Ray Neutron Sensing (CRNS) constitutes an emerging method for monitoring soil moisture and snow dynamics at intermediate spatial scales of several hectares. In complex environments such as mountain regions, however, the presence of areas with a high contrast of hydrogen content was found to cause a hysteresis in the relationship between neutron counts and water equivalent. A simulation study using the newly developed hierarchical scenario tool YULIA (Your URANOS Layer Integration Assistant) for the Monte-Carlo neutron simulation model URANOS was conducted to quantify the effect of snow-free areas on above-ground neutron sensing of the snow water equivalent (SWE). It was found that the size and distance of the snow free patches have the largest impact on the neutron flux. The simulations also showed a sensitivity of the signal towards soil moisture and SWE. Correction functions were developed and validated with observed CRNS measurements and LiDAR based distributed SWE maps. The main aim of the correction procedure is to estimate SWE under partly snow-covered conditions. Furthermore, also the soil moisture of the snow-free areas can be inferred if the SWE distribution is known. The latter can be used for other high-contrast CRNS applications like monitoring soil moisture in the presence of ponding water.

How to cite: Schattan, P., Schmieder, J., Köhli, M., Fey, C., and Schrön, M.: Applying Cosmic-Ray Neutron Sensing in Highly Heterogeneous Conditions: Monitoring Snow Water Equivalent in Periods with Partial Snow Coverage, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18177, https://doi.org/10.5194/egusphere-egu24-18177, 2024.

In the context of climate warming it is a common scientific goal to study and monitor surface and volume changes of glaciers and melting dynamics of its snow and ice. Therefore several measurement techniques exist to track permanently ice melting e.g. DGPS stations on glaciers, Smart stake, and snow and ice depth measurements via e.g. ultrasonic depth sensors to create time series of snow and ice loss or gain. None of the existing methods measure if actually liquid water is present and melting occurs, this is later concluded by interpretation of the geometric data. The capability of the laser sensor to do so via the reflectance value, in fact the received signal intensity, we consider as a big advantage and worth investigating further as a direct measure of snow or ice melt that helps not only to analyze glacier dynamics but is also important e.g. for providing reliable ground truth data for satellite remote sensing. When melting of snow and ice occurs, water changes the reflectance properties as due to absorption of the laser in water, only a portion of the laser is reflected. This allows determining if liquid water is present at the surface measured. We present the data collected in the last 2 melting seasons of the Austre Lovénbreen glacier near Ny Alesund, Svalbard. We show how we classify wet snow and wet ice hours with confidence and are able to compute melting rates. The single point measurement is put into context to area wide LiDAR measurements and melting dynamics of the glacier are analyzed. The data was verified against visual inspections from automatic cameras, data from an automatic weather station both located in the glacier catchment and ice melt was measured in close proximity with a SmartStake station.

How to cite: Prokop, A., Tolle, F., Friedt, J.-M., and Bernard, E.: Using received laser signal intensity to measure snow and ice surface properties automatically , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17917, https://doi.org/10.5194/egusphere-egu24-17917, 2024.

Abstract: Accurate estimation of glacier volume-to-mass conversion relies on a thorough understanding of firn density, both in-depth and over time. Ground-penetrating radar (GPR) serves as a suitable geophysical tool to trace internal reflection horizons (IRHs) and estimate the physical properties of different layers. Our goal is to characterize the IRHs as annual layers and ascertain the spatial firn density-depth profile in the accumulation zone of the Aletsch glacier.

The process involves identifying IRHs from radargrams and iteratively selecting the annual layers by excluding unreasonable layer structures. For an accurate estimation of firn density distribution, it is necessary to derive the velocity-depth profile of electromagnetic waves within the firn zone. The common mid-point (CMP) method was applied to track the velocity distribution within the firn body. Additionally, a method was introduced to estimate the velocity-depth profile for longer GPR profiles by backtracking the calculated velocity from the CMP gather.

To validate IRHs as annual firn layers, we utilized annual accumulation measurements at a nearby stake for Snow Water Equivalent (SWE) estimation. The resulting firn density-depth profile was compared to different firn densification models, considering regional meteorological information. This approach enables us to determine a reliable density-depth function for bulk SWE computations. The study also addresses uncertainties associated with selecting IRHs as annual layers and enhances the application of local volume-to-mass estimates.

How to cite: Patil, A., Mayer, C., and Braun, M.: Firn Density Distribution and Annual Snow Water Equivalent Estimates from Ground Penetrating Radar, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5751, https://doi.org/10.5194/egusphere-egu24-5751, 2024.

Geothermal heat flow (GHF) is a key element of Solid Earth-cryosphere interactions. However, polar regions as Antarctica are only sparsely covered with heat flow determinations from boreholes, so one must rely on interpolation or regression models of GHF (e.g. machine learning) from other sources to derive a regional map. Interpolation/regression of GHF in this manner depends strongly on the available sparse boreholes, which can distort the resulting regional map.

Additional information can be gained from the SMOS (Soil Moisture and Ocean Salinity) satellite by inferring ice temperature profiles with a Bayesian inversion from remote sensing microwave radiometer data. This retrieval uses geothermal heat flow as a free parameter so that it provides a posterior distribution of the GHF needed to explain the ice temperature profiles.

We aim to reconcile geophysical geothermal heat flow models with the ice temperature profiles and improve the estimates of GHF with this coupling.

We use stationary thermal modelling where we force the ice temperature and lithospheric temperature model to converge at the base of the ice. Using stochastic inversion, we estimate the thermal parameters in the lithosphere. The posterior distribution of the retrieval as constraint for the GHF is included as prior distribution to the inversion to the stationary thermal modelling so that the GHF with the highest likelihood can be estimated.

With our approach, we can evaluate a GHF distribution that both explains the ice temperature and lithospheric temperature models and covers large parts of Antarctica.

How to cite: Freienstein, J., Szwillus, W., Leduc-Leballeur, M., Macelloni, G., and Ebbing, J.: Coupling Solid Earth and ice temperature models to estimate geothermal heat flow, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7654, https://doi.org/10.5194/egusphere-egu24-7654, 2024.

The sedimentary units beneath the Ross Ice Shelf in the vicinity of the Kamb Ice Stream grounding line on the Siple Coast of the eastern Ross Ice Shelf play an important role in evaluating past advances and retreats of grounded ice in West Antarctica through the Quaternary. This region is an ongoing focus for drilling efforts that involve melting through the ice shelf and recovering sediments from beneath the seafloor. Seismic (and to a lesser extent gravity) methods have played a critical role in establishing a stratigraphic framework for these sediment sampling endeavours. Approximately 73 km of seismic data have been collected in this region during three seasons since early 2015, complemented by finely sampled gravity transects and a coarser regional gravity grid. Data acquisition provides localised coverage of the sub-ice-shelf ocean and sediments in a region where ROSETTA-Ice airborne-gravity data identified a gravity low. Seismic acquisition parameters have varied from survey to survey, but all involve explosive charges frozen into a hot-water-drilled holes that are recorded by conventional geophones buried in the firn. Such an acquisition configuration provides imaging of the ice shelf and underlying geological units. Processed seismic data show a mostly flat layered seafloor lying beneath the ocean cavity with at least 200 m of sub-horizontally layered sedimentary strata containing several mappable unconformities that are identified as distinct reflective horizons in the seismic data as well as reflection terminations and pinchouts in overlying and underlying units. These unconformities could correspond to past glacial erosion episodes as the position of the grounding line in this region has migrated landward and oceanward. Gravity modelling suggests that the thickness of the sedimentary basins in the region are variable beyond what we see in the shallow (few hundred metre) penetration of the seafloor.

How to cite: Gorman, A., Wilson, G., Horgan, H., Dunbar, G., Hall, C., Black, J., Dagg, B., Tankersley, M., and van Haastrecht, L.: Using seismic and gravity data to constrain subglacial seafloor stratigraphy in the vicinity of the Kamb Ice Stream grounding line, Ross Ice Shelf, Antarctica, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14420, https://doi.org/10.5194/egusphere-egu24-14420, 2024.

Geophysical monitoring becomes more and more popular in permafrost environments due to its remarkable success to detect permafrost thawing and spatio-temporal changes in the ground ice content. Mostly geoelectric methods such as Electrical Resistivity Tomography (ERT) are applied due to the strong differences in the electrical properties between frozen and unfrozen state. However, seismic properties also change markedly upon freezing/thawing and time-lapse refraction seismic tomography (RST) has been shown to be applicable to permafrost over smaller time scales (e.g., Hilbich 2010). The reason why only few studies employ long-term seismic monitoring in permafrost is probably due to the higher logistical effort required.

At two Swiss permafrost monitoring sites (Schilthorn and Stockhorn) yearly RST surveys are conducted using the same setup for more than 15 years, in addition to standard borehole temperature, climatic and ERT measurements (www.permos.ch). The monitoring aim is to image the interannual changes of the thickness of the active layer as well as differences in ice content within the permafrost layer below.

Additional long-term observations are available from RST (and contemporary ERT) surveys from several mountain permafrost sites in Norway that were initially conducted to characterise permafrost conditions around boreholes drilled in 1999/2008 (Juvvasshoe/Jotunheimen), and 2007/2008 (Iskoras/Finnmark, Guolasjavri/Troms, and Tronfjell, cf. Isaksen et al. 2011, Farbrot et al. 2013). These surveys were repeated with the same geometry in 2019 after 11 years in northern Norway, and after 8 and 20 years in southern Norway. As for the Swiss sites, temperatures from all these boreholes show a clear warming trend over the last 1-2 decades (Etzelmüller et al, 2020, 2023).

We here present the observed long-term changes in electrical resistivity and seismic P-wave velocity based on a) annually repeated measurements in the Swiss Alps, and b) on long-term repetition in northern and southern Norway. The geophysical changes are related to the observed borehole temperature increase during the same period (Etzelmüller et al. 2023) and analysed with respect to climate-induced thawing. We evaluate the advantages and disadvantages of seismic monitoring compared to the more standard ERT monitoring. Finally, the results are also analysed with respect to their suitability for future ERT-seismic joint inversion approaches in a monitoring context.

References

Etzelmüller B, Guglielmin M, Hauck C, Hilbich C, Hoelzle M, Isaksen K, Noetzli J, Oliva M and Ramos M 2020. Twenty years of European mountain permafrost dynamics—the PACE legacy. Environ. Res. Lett. 15 104070 DOI 10.1088/1748-9326/abae9d

Etzelmüller B, Isaksen K, Czekirda J, Westermann S, Hilbich C, Hauck C 2023. Rapid warming and degradation of mountain permafrost in Norway and Iceland. The Cryosphere. 17.5477-5497.10.5194/tc-17-5477-2023.

Farbrot H, Isaksen K, Etzelmüller B, Gisnås K 2013. Ground Thermal Regime and Permafrost Distribution under a Changing Climate in Northern Norway. Permafrost Periglac.,24(1):20-38. https://doi.org/10.1002/ppp.1763

Isaksen K, Ødegård RS, Etzelmüller B, Hilbich C, Hauck C, Farbrot H, Eiken T, Hygen HO, Hipp T 2011. Degrading mountain permafrost in southern Norway - spatial and temporal variability of mean ground temperatures 1999-2009. Permafrost Periglac.,22(4):361-377, https://doi 10.1002/ppp.728.

Hilbich C 2010. Time-lapse refraction seismic tomography for the detection of ground ice degradation, The Cryosphere, 4, 243–259, https://doi.org/10.5194/tc-4-243-2010, 2010.

How to cite: Hilbich, C., Etzelmüller, B., Isaksen, K., Mollaret, C., Morard, S., Pellet, C., and Hauck, C.: Long-term refraction seismic monitoring: a reliable method to detect ground ice loss at mountain permafrost sites, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15907, https://doi.org/10.5194/egusphere-egu24-15907, 2024.

Macro-scale transport properties (e.g., electrical conductivity, effective excess charge density and hydraulic conductivity) can be conceptualized as capillary bundle models, in which the pore structure of porous medium is viewed as a bundle of capillary tubes of varying sizes. This approach can be used to understand and address the relationship between the petrophysical properties and the geometry of soil phases. When the temperature of porous medium decreases below the freezing temperature, the soil physical properties (transport properties) change drastically. This is attributed to the complexity of the heterogeneous formation of ice in the porous medium. Therefore, understanding better pore ice formation from microscale insights is crucial to describe the evolution of electrical conductivity with temperature in frozen porous medium. In this study, we consider that capillary radius and tortuous length follow fractal distributions, and that total conductance at the microscale scale is determined by the Gibbs-Thomson and Young-Laplace effects as well as by the surface complexation model. A new capillary bundle model is then proposed using an upscaling procedure, which considers the effects of both bulk and surface conductions. Based primarily on an electrical resistance apparatus and the NMR method, a series of laboratory experiments are carried out to study the influence of initial water saturation and salinity on electrical conductivity under unfrozen and frozen conditions. Additionally, the rationality and validity of the proposed model were successfully verified with published data in the literature and experimental data of this study. Our new physically-based model for electrical conductivity opens up new possibilities to interpret electrical and electromagnetic monitoring to easily infer changes in key variables such as liquid water content and moisture gradients.

How to cite: Luo, H., Jougnot, D., Jost, A., Mendieta, A., and Thanh, L. D.: A physically-based fractal model for predicting the electrical conductivity in partially saturated frozen porous media, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11729, https://doi.org/10.5194/egusphere-egu24-11729, 2024.

Central Asian Mountain regions (Tien Shan and Pamir) are expected to be significantly impacted by climate change, affecting water availability and natural hazards. The cryosphere plays a crucial role in many watersheds of the region by providing water for hydropower station, irrigation, and domestic use downstream. At the same time, retreating glaciers and thawing permafrost increase the risk of natural hazards. Therefore, cryosphere monitoring systems are necessary to provide baseline data for estimating future water availability and detecting dangerous hazard zones. Despite the large areas underlain by permafrost in the Tien Shan and Pamir Mountain ranges, data on permafrost distribution, characteristics and evolution are scarce. However, quantitative estimations of permafrost subsurface components, especially water and ice contents, are needed to evaluate the consequences of current climate change on mountain permafrost environments.

Recent field-based investigations have emphasised the coupled use of geophysical techniques, e.g., by employing the Petrophysical Joint Inversion scheme (PJI, Wagner et al., 2019) that combines electrical resistivity and seismic refraction p-wave velocity data to estimate the four phases present in the subsurface (volumetric contents of air, water, ice, and rock). The traditional PJI implementation relies on Archie’s law (Archie, 1942) as one of the primary petrophysical equation to link resistivity to porosity and water content. Archie's law is generally considered valid when electrolytic conduction dominates, a condition that is not universally justified for dry and coarse blocky substrates and landforms in mountainous terrain. Recognizing this limitation, Mollaret et al. (2020) introduced the electrical Geometric Mean Model as an alternative implementation in the PJI. The Geometric Mean Model  assumes random distributions of the four phases and offers the advantage of including the fractions of ice and air in the petrophysical equation for resistivity, which are not present in Archie’s law. In this study, we assess the feasibility and effectiveness of using the Geometric Mean Model within the PJI framework across an extensive geophysical dataset comprising 22 profiles in Central Asia (Kyrgyzstan and Tajikistan). Our research encompasses diverse landforms, including moraines, rock glaciers, talus slopes, and fine-grained sediments. Our goals are to (i) evaluate the performance of the Geometric Mean Model in comparison to Archies law across different landforms and (ii) address the existing data gap concerning mountain permafrost and ground ice contents in the Central Asian region.

References

Archie, G. E. (1942). The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics. Transactions of the AIME, 146(01), 54–62. https://doi.org/10.2118/942054-G

Mollaret, C., Wagner, F. M., Hilbich, C., Scapozza, C., & Hauck, C. (2020). Petrophysical Joint Inversion Applied to Alpine Permafrost Field Sites to Image Subsurface Ice, Water, Air, and Rock Contents. Frontiers in Earth Science, 8, 85. https://doi.org/10.3389/feart.2020.00085

Wagner, F. M., Mollaret, C., Günther, T., Kemna, A., & Hauck, C. (2019). Quantitative imaging of water, ice and air in permafrost systems through petrophysical joint inversion of seismic refraction and electrical resistivity data. Geophysical Journal International, 219(3), 1866–1875. https://doi.org/10.1093/gji/ggz402

How to cite: Mathys, T., Hilbich, C., Mollaret, C., Hauck, C., Saks, T., Usubaliev, R., Moldobekov, B., Bektursunov, Z., Azimshoev, M., Navruzshoev, H., and Hoelzle, M.: Quantifying Ground Ice in Tien Shan and Pamir Permafrost: A Comprehensive Petrophysical Joint Inversion Study Applying the electrical Geometric Mean Model , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-21232, https://doi.org/10.5194/egusphere-egu24-21232, 2024.

Electrical resistivity tomography (ERT) is one of the most accurate geophysical techniques to distinguish between frozen and unfrozen ground in permafrost areas. Performing the measurements, however, requires considerable logistics and time efforts. This is mainly due to the fact that optimal galvanic contact between the electrodes and the ground surface is necessary to collect reliable ERT datasets. Therefore, the traditional steel-spike electrodes must be steadily coupled between the boulders and wet with salt water on coarse blocky surfaces. To further decrease the contact resistances, sponges soaked in salt water can be inserted between the spike and the surface of rocks. Nevertheless, this traditional coupling system is particularly time-consuming, making it challenging to collect several ERT survey lines in a single workday in mountain environments. Recently developed conductive textile electrodes were applied to facilitate the deployment of ERT arrays in rock glacier environments. Instead of hammering the steel spikes, the conductive textile electrodes can be easily pushed between the boulders and wet with less water (compared to the sponges). Consequently, this new electrode approach decreases the time needed to prepare an ERT array. In this work, we evaluate the performance of the textile electrodes by comparing these with the traditional electrode approach, considering common investigation lines. This comparative test has been carried out in three test sites, which present different lithologies, surface characteristics and using different electrode spacing. The collected datasets were statistically analysed with robust regression analysis and Wilcoxon rank-sum test to examine the accuracy and significant differences between the two electrode systems regarding contact resistances, injected electrical current, measured apparent resistivities, reciprocal error, and inverted resistivity values. The obtained results demonstrate that conductive textile electrodes are suitable to collect reliable ERT datasets and, consequently, applying this approach in future ERT measurements performed in high mountain environments with coarse blocky surfaces (e.g. rockfall deposits, blocky slopes, or rock glaciers) would allow to acquire more survey lines (e.g. realisation of pseudo-3D geometries) extending the characterisation of the subsurface.

How to cite: Pavoni, M., Boaga, J., Bast, A., Lichtenegger, M., and Buckel, J.: Conductive textile electrodes for time-efficient ERT surveys performed in coarse-blocky mountain environments, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4532, https://doi.org/10.5194/egusphere-egu24-4532, 2024.

Permafrost peatlands, located in the Arctic and high mountain regions, are typically known to be ice-rich. This is primarily linked to significant water content and often oversaturation, a characteristic property of peat soil. The current understanding of the effects of human-induced climate warming suggests that these regions are approaching a climatic tipping point with substantial permafrost thaw expected in the coming decades. Ice content is an important parameter for modelling permafrost evolution and at present limited studies exist that determine its in-situ spatial distribution in such areas.

The geophysical method known as high-frequency induced polarisation (HFIP) is advantageous for cryohydrological research in these environments. This method can capture the frequency-dependent polarisation of ice (also termed dielectric relaxation peak), which occurs within the range of 100 Hz to 100 kHz and is expressed by complex resistivity. Therefore, by analysing the spectral behaviour of this complex resistivity within the target frequency range the distribution and quantity of ice can be estimated.

The results from the latest field campaign conducted at Storflaket mire and Stordalen mire in Abisko, Sweden, are presented. Two-dimensional HFIP profiles were measured to resolve the near-surface unfrozen layer (no-ice) and the underlying frozen layer (ice-bearing). The measurements were performed in late summer when the depth of the unfrozen layer was at its maximum. Field data are inverted as independent frequencies to obtain the spectral variation of complex resistivity. No-ice and ice-bearing regions are classified by the presence of the relaxation peak. Subsequently, a two-component mixture model, with one component as ice and the second as the surrounding matrix, is applied to determine ice content distribution. Boundary constraints and starting parameters are chosen using the spectral analysis of the inverted complex resistivity. The model accuracy is evaluated using unfrozen layer probing and a permafrost core extracted along the HFIP profile. The HFIP-derived ice content distribution is consistent with unfrozen layer probing, i.e., the classification of no-ice and ice-bearing regions is successful. The model tends to underestimate ice content percentages compared to permafrost core laboratory measurements. This discrepancy can be explained since laboratory measurements are based on gravimetric water content and assumes all pore-water is frozen. However, it is known that residual pore-water is present in these soils even below 0°C. Additionally, it is observed that the model performs well when the ice content percentage is 10% or greater and its applicability might be limited in scenarios where the ice content is less than 10%.

The latest results are discussed in comparison with previous findings from Heliport, a permafrost mire also located in Abisko. In the Heliport study, HFIP successfully resolved the complex resistivity and ice content distribution on a larger scale. Building on the field knowledge gained at Heliport, this study incorporates improvements in electrode configuration setup, data acquisition speed, and minimising cable-earth coupling effects. The findings contribute to the understanding of the induced polarisation of permafrost peatlands, which is an underexplored area from a geophysical perspective.

How to cite: Sugand, M. and Hördt, A.: Broadband Spectral Induced Polarization in Permafrost Peatlands of Northern Sweden , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9857, https://doi.org/10.5194/egusphere-egu24-9857, 2024.