Records of Alpine microseismicity are a powerful tool to study landscape-shaping processes and warn against hazardous mass movements. Unfortunately, seismic sensor coverage in Alpine regions is typically insufficient. Here we show that distributed acoustic sensing (DAS) bridges critical observational gaps of seismogenic processes in Alpine terrain. Dynamic strain measurements in a 1 km long fiber optic cable on a glacier surface produce high-quality seismograms related to glacier flow and nearby rock falls. The nearly 500 cable channels precisely locate a series of glacier stick-slip events (within 20-40 m) and reveal seismic phases from which thickness and material properties of the glacier and its bed can be derived. As seismic measurements can be acquired with fiber optic cables that are easy to transport, install and couple to the ground, our study demonstrates the potential of DAS technology for seismic monitoring of glacier dynamics and natural hazards.

Recent work in the field of cryo-seismology demonstrates that high-frequency (>1  Hz) seismic waves provide key constraints on a wide range of glacier processes, such as basal friction, surface crevassing, or subglacial water flow. Establishing quantitative links between the seismic signal and the processes of interest, however, requires detailed characterization of the wavefield, which, at high frequencies, necessitates the deployment of large and dense seismic arrays. Although dense seismic array monitoring has recently become increasingly common in geophysics, its application to glaciated environments remains limited. Here, we present a dense seismic array experiment made of 98 three-component seismic stations continuously recording during 35 days in early spring 2018 on the Argentière Glacier, French Alps. The seismic dataset is supplemented with a wide range of complementary observations obtained from ground-penetrating radar, drone imagery, Global Navigation Satellite Systems positioning, and in situ measurements of basal glacier sliding velocities and subglacial water discharge. We present first results through conducting spectral analysis, template matching, matched-field processing, and eikonal-wave tomography. We report enhanced spatial resolution on basal stick slip and englacial fracturing sources as well as novel constraints on the heterogeneous nature of the noise field generated by subglacial water flow and on the link between crevasse properties and englacial seismic velocities. We outline in which ways further work using this dataset could help tackle key remaining questions in the field.

Crevasses and englacial fracture networks route meltwater from a glacier's surface to the subglacial drainage system and thus influence glacial hydraulics. However, rapid fracture growth may also lead to sudden and potentially hazardous structural failure of unstable glaciers and ice dams, rifting of ice shelves, or iceberg calving. Here, we use passive seismic recordings to investigate the englacial fracture network on Glacier de la Plaine Morte, Switzerland. Glacier dynamics and the drainage of an ice-marginal lake give rise to numerous icequakes, the majority of which generate dispersed, high-frequency Rayleigh waves. A wide distribution of events allows us to study azimuthal anisotropy between 10 and 30 Hz in order to extract englacial seismic velocities in regions of preferentially oriented crevasses. Beamforming applied to a 100-m-aperture array reveals azimuthal anisotropy of Rayleigh-wave phase velocities reaching a strength of 8% at high frequencies. In addition, we find that the fast direction of wave propagation coincides with the observed surface strike of the narrow crevasses. Forward modeling and inversion of dispersion curves suggest that the azimuthal anisotropy is induced by a 40-m-thick crevassed layer at the surface of the glacier with 8% anisotropy in shear-wave velocity.

Abstract. Hydraulic processes impact viscous and brittle ice deformation. Water-driven fracturing as well as turbulent water flow within and beneath glaciers radiate seismic waves which provide insights into otherwise hard-to-access englacial and subglacial environments. In this study, we analyze glaciohydraulic tremors recorded by four seismic arrays installed in different parts of Glacier de la Plaine Morte, Switzerland. Data were recorded during the 2016 melt season including the sudden subglacial drainage of an ice-marginal lake. Together with our seismic data, discharge, lake level, and ice flow measurements provide constraints on glacier hydraulics. We find that the tremors are generated by subglacial water flow, in moulins, and by icequake bursts. The dominating process can vary on sub-kilometer and sub-daily scales. Consistent with field observations, continuous source tracking via matched-field processing suggests a gradual up-glacier progression of an efficient drainage system as the melt season progresses. The ice-marginal lake likely connects to this drainage system via hydrofracturing, which is indicated by sustained icequake signals emitted from the proximity of the lake basin and starting roughly 24 h prior to the lake drainage. To estimate the hydraulics associated with the drainage, we use tremor–discharge scaling relationships. Our analysis suggests a pressurization of the subglacial environment at the drainage onset, followed by an increase in the hydraulic radii of the conduits and a subsequent decrease in the subglacial water pressure as the capacity of the drainage system increases. The pressurization is in phase with the drop in the lake level, and its retrieved maximum coincides with ice uplift measured via GPS. Our results highlight the use of cryo-seismology for monitoring glacier hydraulics.

Abstract. Hydraulic processes impact viscous and brittle ice deformation. Water-driven fracturing as well as turbulent water flow within and beneath glaciers radiate seismic waves which provide insights into otherwise hard-to-access englacial and subglacial environments. In this study, we analyze glaciohydraulic tremors recorded by four seismic arrays installed in different parts of Glacier de la Plaine Morte, Switzerland. Data were recorded during the 2016 melt season including the sudden subglacial drainage of an ice-marginal lake. Together with our seismic data, discharge, lake level, and ice flow measurements provide constraints on glacier hydraulics. We find that the tremors are generated by subglacial water flow, in moulins, and by icequake bursts. The dominating process can vary on sub-kilometer and sub-daily scales. Consistent with field observations, continuous source tracking via matched-field processing suggests a gradual upglacier progression of an efficient drainage system as the melt season progresses. The ice-marginal lake likely connects to this drainage system via hydrofracturing, which is indicated by sustained icequake signals emitted from the proximity of the lake basin and starting roughly 24 hours prior to the lake drainage. To estimate the hydraulics associated with the drainage, we use tremor-discharge scaling relationships. Our analysis suggests a pressurization of the subglacial environment at the drainage onset, followed by an increase in the hydraulic radii of the conduits and a subsequent decrease in the subglacial water pressure as the capacity of the drainage system increases. The pressurization is in phase with the drop in the lake level and its retrieved maximum coincides with ice-uplift measured via GPS. Our results highlight the use of cryo-seismology for monitoring glacier hydraulics.

Sliding of glacial ice over its base is typically the dominant flow mechanism of glaciers and ice sheets. The term "sliding" is often used in a loose sense to include the deformation of "soft" till beds and the differential motion between basal ice and underlying bedrock. The amount of sliding controls fast (Clarke, 1987) and slow (Maier et al., 2019) ice flow regimes. Sliding variations are powerful enough to halt, or even reverse, ice stream flow (Conway et al., 2002) or rupture ice tongues, leading to massive ice avalanches (Faillettaz et al., 2015; Kääb et al., 2018). Nevertheless, sliding mechanisms remain elusive, leading to uncertain predictions of ice sheet motion and stability, and ultimately sea level rise (

Abstract. Ambient noise seismology has revolutionized seismic characterization of the Earth's crust from local to global scales. The estimate of Green's function (GF) between two receivers, representing the impulse response of elastic media, can be reconstructed via cross-correlation of the ambient noise seismograms. A homogenized wave field illuminating the propagation medium in all directions is a prerequisite for obtaining an accurate GF. For seismic data recorded on glaciers, this condition imposes strong limitations on GF convergence because of minimal seismic scattering in homogeneous ice and limitations in network coverage. We address this difficulty by investigating three patterns of seismic wave fields: a favorable distribution of icequakes and noise sources recorded on a dense array of 98 sensors on Glacier d'Argentière (France), a dominant noise source constituted by a moulin within a smaller seismic array on the Greenland Ice Sheet, and crevasse-generated scattering at Gornergletscher (Switzerland). In Glacier d'Argentière, surface melt routing through englacial channels produces turbulent water flow, creating sustained ambient seismic sources and thus favorable conditions for GF estimates. Analysis of the cross-correlation functions reveals non-equally distributed noise sources outside and within the recording network. The dense sampling of sensors allows for spatial averaging and accurate GF estimates when stacked on lines of receivers. The averaged GFs contain high-frequency (>30 Hz) direct and refracted P waves in addition to the fundamental mode of dispersive Rayleigh waves above 1 Hz. From seismic velocity measurements, we invert bed properties and depth profiles and map seismic anisotropy, which is likely introduced by crevassing. In Greenland, we employ an advanced preprocessing scheme which includes match-field processing and eigenspectral equalization of the cross spectra to remove the moulin source signature and reduce the effect of inhomogeneous wave fields on the GFs. At Gornergletscher, cross-correlations of icequake coda waves show evidence for homogenized incident directions of the scattered wave field. Optimization of coda correlation windows via a Bayesian inversion based on the GF cross coherency and symmetry further promotes the GF estimate convergence. This study presents new processing schemes on suitable array geometries for passive seismic imaging and monitoring of glaciers and ice sheets.

No abstract