On the Green’s function emergence from interferometry of seismic wavefields generated in high-melt glaciers: implications for passive imaging and monitoring
Abstract. Ambient noise seismology has revolutionized seismic characterization of the Earth's crust from local to global scales. The elastic Green's function (GF) between two receivers can be reconstructed via cross-correlation of the ambient noise seismograms. An homogenized wavefield illuminating the propagation medium in all directions is a pre-requesite for obtaining accurate GF. For seismic data recorded on glaciers, this condition imposes strong limitations on GF convergence, because of minimal seismic scattering in homogeneous ice. We address this difficulty by investigating three patterns of seismic wavefields: a favourable 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. From the velocity measurements of reconstructed Rayleigh waves, we invert bed properties and depth profiles, and map seismic anisotropy, which is likely introduced by crevassing. In Greenland, we employ an advanced pre-processing scheme which include match-field processing and eigenspectral equalization of the cross-spectra to remove the moulin source signature and reduce the effect of inhomogeneous wavefields on the GF. At Gornergletscher, cross-correlations of icequake coda waves show evidence for homogenized wavefields. Optimization of coda correlation windows further promotes the GF convergence. This study presents new processing schemes on suitable array geometries for passive seismic imaging and monitoring of glaciers.