Inter-station coda wavefield studies using a novel icequake database on Erebus volcano

Recent theoretical advances pertaining to the properties of multiply scattered wavefields have yielded a plethora of numerical and controlled source studies aiming to better understand what information may be derived from these otherwise chaotic signals. Practically, multiply scattered wavefields are difficult to compare to numerically derived models due to a combination of source paucity/directionality and array density limitations, particularly in passive seismology scenarios. Furthermore, in situations where data quantities are abundant, such as for ambient noise correlations, it remains very difficult to recover pseudo-Green's function symmetry in the ballistic components of the wavefield, let alone in the coda of the correlations. In this study, we use a large network of short period and broadband instruments on Erebus volcano to show that actual Green's function recovery is indeed possible in some cases. We make use of a large database of small impulsive icequakes distributed randomly on the summit plateau and, using fundamental theoretical properties of equipartitioned wavefields and interstation icequake coda correlations, are able to directly derive notoriously difficult quantities such as the bulk elastic mean free path for the volcano, demonstrations of correlation coda symmetry and its dependence on the number of icequakes used, and a theoretically predicted coherent backscattering amplification factor associated with weak localization. We furthermore show that stable equipartition and H^2/V^2 ratios may be consistently observed for icequake coda, and we perform simple depth inversions of these frequency dependent quantities to compare with known structures.