Hanging glacier monitoring with icequake repeaters and seismic coda wave interferometry: a case study of the Eiger hanging glacier

Driven by the force of gravity, hanging glacier instabilities can lead to catastrophic rupture events. Reliable forecasting remains a challenge as englacial damage leading to large-scale failure is masked from modern sensing technology focusing on the ice surface. The Eiger hanging glacier, located in the Swiss Alps, was intensely monitored between April and August 2016 before a moderate 15,000 m3 break-off event from the ice cliff. Among different instruments, such as an automatic camera and interferometric radar, four 3-component seismometers were installed on the glacier. A single seismometer operated throughout the whole monitoring period. It recorded over 200,000 repeating icequakes showing strong englacial seismic coda waves. We propose a novel approach for hanging glacier monitoring by combining repeating icequake analysis, coda wave interferometry, and attenuation measurements. Our results show a seasonal 0.1 % decrease in relative englacial seismic velocity dv/v and an increase in coda wave attenuation Qc−1 (Qc decreases from ~50 to ~30). Comparison of dv/v and Qc with air temperature suggests that these changes are driven by a seasonal increase in the glacier’s ice and firn pack temperature that might affect the top 20 m of the glacier. Diurnal cycles of Qc−1, repeating icequake activity, and the velocity of the glacier front shift from cosinusoidal to sinusoidal variations under the presence of meltwater. The proposed approach extends the monitoring of the hanging glacier beyond the ice surface and allows for a better understanding of the glacier’s response to time-dependent external forcing, which is an important step towards improved break-off forecasting systems.

Triggering Mechanisms of Gayari Avalanche, Pakistan

A massive snow avalanche occurred on April, 2012 at Gayari, located in NE part of Pakistan, close to India and China Border. The catastrophic avalanche killed nearly 148 people, majority of which were Pakistan army personnel destroying army base camp. To mitigate its future hazard, different triggering mechanisms have been investigated in this study. We contemplate that the avalanche was triggered due to snow pack existence on favorable slope in combination with different meteorological conditions and anomalous ground vibration. The avalanche occurrence clock was advanced by two earthquakes: M4.1 at a distance ∼ 125 km that occurred about 21 hours before and another comparatively larger (M5.6) earthquake that occurred comparatively at larger distance (∼ 370 km) and longer time (∼ 25 days) before which have significantly changed the loading conditions. The latter event (M 5.6) has imparted maximum peak dynamic stress and cumulative seismic moment a month before the avalanche. Interestingly the avalanche occurred within the seismic coda of M2.8 earthquake from Hindu Kush region, located at 560 km distance. Although the size and its expected impact on avalanche might be minor but its role in instantaneous triggering cannot be ruled out. Even smaller events at larger distance have been reported to cause snow avalanches in same environments. The presence of cracks within the avalanche, were further weaken by persistence of extremely low temperature (lowest in the past decade), causing high precipitation rate along with altering the mechanical properties of the weak layer within the snow pack. Robust wind pressure pattern highest and lowest in March and April, 2012 respectively might be responsible for abrupt changes in loading conditions.

Source Separation and Medium Change of Contained Chemical Explosions from Coda Wave Interferometry

Differences in the seismic coda of neighboring events can be used to investigate source location offsets and medium change with coda wave interferometry (CWI). We employ CWI to infer the known relative location between two chemical explosions in Phase I of the Source Physics Experiment (SPE). The inferred displacement between the first, SPE-1, and second, SPE-2, chemical explosion is between 6 and 18 m, with an expectation of 9.2 m, where the known separation is close to 9.4 m. We also employ CWI to find any velocity perturbation due to damage from SPE-2, by comparing its coda with the collocated third SPE chemical explosion, SPE-3. We find that damage due to SPE-2 must be confined to a spherical region with radius less than 10 m and velocity perturbation less than 25%.

Shear wave velocity from inter-source interferometry

<p>Ambient noise correlation uses the recordings of multiple, statistically distributed seismic sources (the noise sources) at two seismometers. By cross-correlating this signal, one obtains a wave traveling between two seismometers. Due to the principle of reciprocity it is possible to interchange the role of sources and receivers. This cannot be done with ambient noise, but another stochastic signal, the seismic coda is used. Using a cross-correlation of the seismic coda of two earthquakes recorded at multiple seismometers, it is possible to construct a seismic wave traveling between the two earthquakes in depth (inter-source interferometry). Here, I use the the time lag of the maxima in the cross-correlation of the coda wave field to measure the shear wave velocity in the source volume of swarm earthquakes. This technique is different from previous studies analyzing the decorrelation of the coda wave field of nearby events or using the cross-correlation for relocation purposes.</p><p>The technique is applied to five event clusters of the 2018 West Bohemia earthquake swarm. With the help of a high quality earthquake catalog, I was able to determine the shear wave velocity in the region of the five clusters separately. The shear wave velocities range between 3.5 km/s and 4.2 km/s. The resolution of this novel method is given by the extent of the clusters and better than for a comparable classical tomography. The method can be incorporated into a tomographic inversion to map the shear wave velocity in the source region with unprecedented resolution. The influence of focal mechanisms and the attenuation properties on the polarity and location of the maxima in the cross-correlation functions is discussed.</p>

Sensitivity kernels for static and dynamic tomography of scattering and absorbing media with elastic waves: a probabilistic approach

SUMMARY Scattered seismic coda waves are frequently used to characterize small scale medium heterogeneities, intrinsic attenuation or temporal changes of wave velocity. Spatial variability of these properties raises questions about the spatial sensitivity of seismic coda waves. Especially the continuous monitoring of medium perturbations using ambient seismic noise led to a demand for approaches to image perturbations observed with coda waves. An efficient approach to localize spatial and temporal variations of medium properties is to invert the observations from different source–receiver combinations and different lapse times in the coda for the location of the perturbations. For such an inversion, it is key to calculate the coda-wave sensitivity kernels which describe the connection between observations and the perturbation. Most discussions of sensitivity kernels use the acoustic approximation in a spatially uniform medium and often assume wave propagation in the diffusion regime. We model 2-D multiple non-isotropic scattering in a random elastic medium with spatially variable heterogeneity and attenuation using the radiative transfer equations which we solve with the Monte Carlo method. Recording of the specific energy density of the wavefield that contains the complete information about the energy density at a given position, time and propagation direction allows us to calculate sensitivity kernels according to rigorous theoretical derivations. The practical calculation of the kernels involves the solution of the adjoint radiative transport equations. We investigate sensitivity kernels that describe the relationships between changes of the model in P- and S-wave velocity, P- and S-wave attenuation and the strength of fluctuation on the one hand and seismogram envelope, traveltime changes and waveform decorrelation as observables on the other hand. These sensitivity kernels reflect the effect of the spatial variations of medium properties on the wavefield and constitute the first step in the development of a tomographic inversion approach for the distribution of small-scale heterogeneity based on scattered waves.

Seismic Coda-Waves Imaging Based on Sensitivity Kernels Calculated Using an Heuristic Approach

In this paper we review and discuss the seismic method based on the analysis of seismic coda waves used in the last 10 years by the present authors and/or their co-workers, to produce separate images of intrinsic- and scattering attenuation in zones of peculiar geological interest (mainly volcanoes). Such separate attenuation images are considered by the scientific community as complementary to those from ordinary velocity-tomography and useful to improve the geological interpretation in volcanoes and in tectonically active zones. In this review we only list but do not discuss the most significative papers showing the images obtained, as we are focused to review the method and not the interpretation of data analysis. For sake of completeness, we anyway show also a new analysis applied to data from Stromboli volcano. We thus first introduce the physical model describing the seismogram Energy Envelope (derived from the solution of the Energy Transport integral Equation) and discuss its asymptotic approximations (Diffusion- and Single-scattering model). Then, we describe a numerical method to heuristically calculate the Sensitivity Kernels for the propagation of the scattered waves in the assumption of isotropic scattering. We attribute to these Sensitivity Kernels the physical meaning of probability that for a single source-receiver couple the measured attenuation parameters can be associated with the space coordinates. Based on this definition, the attenuation image can be obtained mapping the estimated attenuation parameters onto the zone under study weighting with the Sensitivity Kernels. We further discuss how to estimate the uncertainties associated with the results and report the list of the papers describing the (separated) scattering- and intrinsic-attenuation structures investigated using this approach.

Understanding seismic path biases and magmatic activity at Mount St Helens volcano before its 2004 eruption

SUMMARY In volcanoes, topography, shallow heterogeneity and even shallow morphology can substantially modify seismic coda signals. Coda waves are an essential tool to monitor eruption dynamics and model volcanic structures jointly and independently from velocity anomalies: it is thus fundamental to test their spatial sensitivity to seismic path effects. Here, we apply the Multiple Lapse Time Window Analysis (MLTWA) to measure the relative importance of scattering attenuation vs absorption at Mount St Helens volcano before its 2004 eruption. The results show the characteristic dominance of scattering attenuation in volcanoes at lower frequencies (3–6 Hz), while absorption is the primary attenuation mechanism at 12 and 18 Hz. Scattering attenuation is similar but seismic absorption is one order of magnitude lower than at open-conduit volcanoes, like Etna and Kilauea, a typical behaviour of a (relatively) cool magmatic plumbing system. Still, the seismic albedo (measuring the ratio between seismic energy emitted and received from the area) is anomalously high (0.95) at 3 Hz. A radiative-transfer forward model of far- and near-field envelopes confirms this is due to strong near-receiver scattering enhancing anomalous phases in the intermediate and late coda across the 1980 debris avalanche and central crater. Only above this frequency and in the far-field diffusion onsets at late lapse times. The scattering and absorption parameters derived from MLTWA are used as inputs to construct 2-D frequency-dependent bulk sensitivity kernels for the S-wave coda in the multiple-scattering (using the Energy Transport Equations—ETE) and diffusive (AD, independent of MLTWA results) regimes. At 12 Hz, high coda-attenuation anomalies characterize the eastern side of the volcano using both kernels, in spatial correlation with low-velocity anomalies from literature. At 3 Hz, the anomalous albedo, the forward modelling, and the results of the tomographic imaging confirm that shallow heterogeneity beneath the extended 1980 debris-avalanche and crater enhance anomalous intermediate and late coda phases, mapping shallow geological contrasts. We remark the effect this may have on coda-dependent source inversion and tomography, currently used across the world to image and monitor volcanoes. At Mount St Helens, higher frequencies and deep borehole data are necessary to reconstruct deep volcanic structures with coda waves.