On the Green's function emergence from interferometry of seismic wave fields generated in high-melt glaciers: implications for passive imaging and monitoring

• Published in: The Cryosphere
• Main Authors: A. Sergeant, M. Chmiel, F. Lindner, F. Walter, P. Roux, J. Chaput, F. Gimbert, A. Mordret
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2020
• Subjects: geo; envir; Greenland; glacier; Ice Sheet; The Cryosphere
• Online Access: https://doi.org/10.5194/tc-14-1139-2020; https://www.the-cryosphere.net/14/1139/2020/tc-14-1139-2020.pdf; https://doaj.org/article/30342200b31445c2b41d1679d2adcbbf

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, ...