USGS ScienceSearch

USGS · 70219526

Reconstructing the dynamics of the highly similar May 2016 and June 2019 Iliamna Volcano, Alaska ice–rock avalanches from seismoacoustic data

Abstract

Surficial mass wasting events are a hazard worldwide. Seismic and acoustic signals from these often remote processes, combined with other geophysical observations, can provide key information for monitoring and rapid response efforts and enhance our understanding of event dynamics. Here, we present seismoacoustic data and analyses for two very large ice–rock avalanches occurring on Iliamna Volcano, Alaska (USA), on 22 May 2016 and 21 June 2019. Iliamna is a glacier-mantled stratovolcano located in the Cook Inlet, ∼200 km from Anchorage, Alaska. The volcano experiences massive, quasi-annual slope failures due to glacial instabilities and hydrothermal alteration of volcanic rocks near its summit. The May 2016 and June 2019 avalanches were particularly large and generated energetic seismic and infrasound signals which were recorded at numerous stations at ranges from ∼9 to over 600 km. Both avalanches initiated in the same location near the head of Iliamna's east-facing Red Glacier, and their ∼8 km long runout shapes are nearly identical. This repeatability – which is rare for large and rapid mass movements – provides an excellent opportunity for comparison and validation of seismoacoustic source characteristics. For both events, we invert long-period (15–80 s) seismic signals to obtain a force-time representation of the source. We model the avalanche as a sliding block which exerts a spatially static point force on the Earth. We use this force-time function to derive constraints on avalanche acceleration, velocity, and directionality, which are compatible with satellite imagery and observed terrain features. Our inversion results suggest that the avalanches reached speeds exceeding 70 m s−1, consistent with numerical modeling from previous Iliamna studies. We lack sufficient local infrasound data to test an acoustic source model for these processes. However, the acoustic data suggest that infrasound from these avalanches is produced after the mass movement regime transitions from cohesive block-type failure to granular and turbulent flow – little to no infrasound is generated by the initial failure. At Iliamna, synthesis of advanced numerical flow models and more detailed ground observations combined with increased geophysical station coverage could yield significant gains in our understanding of these events.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 59.712097173322924° to 63.31268278043484° latitude; -156.26953125° to -144.8876953125° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Liam Toney, David Fee, Kate E. Allstadt, Matthew M. Haney, Robin S. Matoza. 2021-04-08. Reconstructing the dynamics of the highly similar May 2016 and June 2019 Iliamna Volcano, Alaska ice–rock avalanches from seismoacoustic data. https://doi.org/10.5194/esurf-9-271-2021

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

River suspended-sand flux computation with uncertainty estimation using water samples and high-resolution ADCP measurements

Measuring suspended-sand fluxes in rivers remains a scientific challenge due to their high spatial and temporal variability. To capture the vertical and lateral gradients of concentration in the cross-section, measurements with point samples are performed. However, the uncertainty related to these measurements is rarely evaluated, as few studies of the major sources of error exist. Therefore, the aim of this study is to develop a method to determine the cross-sectional sand flux and estimate its uncertainty. This SDC (for sand discharge computing) method combines suspended-sand concentrations from point samples with ADCP (acoustic Doppler current profiler) high-resolution depth and velocity measurements. The MAP (for multitransect averaged profile) method allows obtaining an average of several ADCP transects on a regular grid, including the unmeasured areas. The suspended-sand concentrations are integrated vertically by fitting a theoretical exponential suspended-sand profile to the data using Bayesian modeling. The lateral integration is based on the water depth as a proxy for the local bed shear stress to evaluate the bed concentration and sediment diffusion along the river cross-section. The estimation of uncertainty combines ISO standards and semi-empirical methods with a Bayesian approach to estimate the uncertainty due to the vertical integration. The new method is applied to data collected in four rivers under various hydro-sedimentary conditions: the Colorado, Rhône, Isère, and Amazon rivers, with computed flux uncertainties ranging between 18 % and 32 %. The relative difference between the suspended-sand flux in 21 cases calculated with the proposed SDC method compared to the International Organization for Standardization (ISO) 4363 standard method ranges between − 40 % and + 23 %. This method that comes with a flexible, open-source code is the first to propose an applicable uncertainty estimation that could be adapted to other flux computation methods.

Earth Surface Dynamics

Implications for the resilience of modern coastal systems derived from mesoscale barrier dynamics at Fire Island, New York

Understanding the response of coastal barriers to future changes in rates of sea level rise, sediment availability, and storm intensity/frequency is essential for coastal planning, including socioeconomic and ecological management. Identifying drivers of past changes in barrier morphology, as well as barrier sensitivity to these forces, is necessary to accomplish this. Using remote sensing, field, and laboratory analyses, we reconstruct the mesoscale (decades–centuries) evolution of central Fire Island, a portion of a 50 km barrier island fronting Long Island, New York, USA. We find that the configuration of the modern beach and foredune at Fire Island is radically different from the system's relict morphostratigraphy. Central Fire Island is comprised of at least three formerly inlet-divided rotational barriers with distinct subaerial beach and dune–ridge systems that were active prior to the mid-19th century. Varying morphologic states reflected in the relict barriers (e.g., progradational and transgressive) contrast with the modern barrier, which is dominated by a tall and nearly continuous foredune and is relatively static, except for erosion and drowning of its fringing marsh. We suggest that this state shift indicates a transition from a regime dominated by inlet-mediated gradients in alongshore sediment availability to one where human impacts exerted greater influence on island evolution from the late 19th century onward. The retention of some geomorphic capital in Fire Island's relict subaerial features combined with its static nature renders the barrier increasingly susceptible to narrowing and passive submergence. This may lead to an abrupt geomorphic state shift in the future, a veiled vulnerability that may also exist in other stabilized barriers.

New York

Steady-state forms of channel profiles shaped by debris flow and fluvial processes

Debris flows regularly traverse bedrock channels that dissect steep landscapes, but our understanding of bedrock erosion by debris flows and their impact on steepland morphology is still rudimentary. Quantitative models of steep bedrock channel networks are based on geomorphic transport laws designed to represent erosion by water-dominated flows. To quantify the impact of debris flow erosion on steep channel network form, it is first necessary to develop methods to estimate spatial variations in bulk debris flow properties (e.g., flow depth, velocity) throughout the channel network that can be integrated into landscape evolution models. Here, we propose and evaluate two methods to estimate spatial variations in bulk debris flow properties along the length of a channel profile. We incorporate both methods into a model designed to simulate the evolution of longitudinal channel profiles that evolve in response to debris flow and fluvial processes. To explore this model framework, we propose a general family of debris flow erosion laws where erosion rate is a function of debris flow depth and channel slope. Model results indicate that erosion by debris flows can explain the occurrence of a scaling break in the slope–area curve at low-drainage areas and that upper-network channel morphology may be useful for inferring catchment-averaged erosion rates in quasi-steady landscapes. Validating specific forms of a debris flow incision law, however, would require more detailed model–data comparisons in specific landscapes where input parameters and channel morphometry can be better constrained. Results improve our ability to interpret topographic signals within steep channel networks and identify observational targets critical for constraining a debris flow incision law.

Earth Surface Dynamics