Geomorphologically-controlled seismic signals at Mount St. Helens volcano

In volcanoes, topography and shallow morphology can substantially modify seismic signals, tracing anisotropic signatures in the crust's most surficial layers. To better understand the influence of key morphologies, forward modelling of the seismic waveforms is fundamental.  Here, we introduce a forward model of the seismic wave equation developed with finite-differences schemes in anisotropic viscoelastic media. The observation of geomorphological features and the surficial geology map of Mount St. Helens are used to reproduce the scattering and anisotropic effects caused by shallow heterogeneity on seismic signals. The main aim is to understand if and to which lengths lateral anisotropic variations in geomorphological features control the generation and propagation of low-frequency seismic signals, focusing especially on the timing of surface-wave enhancement.

The model shows how the geomorphology-derived anisotropy controls the travel times of the horizontally polarized S waves (SH), in particular along with two directions: WNW-ESE, following the trend of a buried fault, and NS, consistent with the main morphological difference between southern (mostly untouched by the 1980 eruption) and northern (collapsed in 1980’s blast) flanks of the volcano. An analysis of the waveforms of a shallow event of 2005 (during the last eruption of Mt. St. Helens), located in the crater, shows how an isotropic model can reproduce the arrival of the SH wave at high frequencies (10 Hz). The introduction of an effective anisotropic medium is necessary to explain the arrivals for stations deployed across the north-northwestern flank of the volcano at lower frequencies (1 Hz and 6 Hz). The heterogeneity in the crater (e.g., the glacier inside the crater covered by a rock-debris layer) can create interfaces made mostly of unconsolidated materials. As also demonstrated by radiative transfer simulation, the crater acts as a primary source of surface waves dominating the seismic signals.