Blast dynamics at Mount St Helens on 18 May 1980

Abstract

At 8.32 a.m. on 18 May 1980, failure of the upper part of the north slope of Mount St Helens triggered a lateral eruption (‘the blast’) that devastated the conifer forests in a sector covering ∼500 km2 north of the volcano. I present here a steady flow model for the blast dynamics and propose that through much of the devastated area the blast was a supersonic flow of a complex multiphase (solid, liquid, vapour) mixture. The shape of the blast zone; pressure, temperature, velocity (Mach number) and density distributions within the flow; positions of weak and strong internal shocks; and mass flux, energy flux, and total energy are calculated. The shape of blast zone was determined by the initial areal expansion from the reservoir, by internal expansion and compression waves (including shocks), and by the density of the expanding mixture. The pressure within the flow dropped rapidly away from the source of the blast until, at a distance of ∼11 km, the flow became underpressured relative to the surrounding atmosphere. Weak shocks within the flow subparallel to the east and west margins coalesced at about this distance into a strong Mach disk shock, across which the flow velocities would have dropped from supersonic to subsonic as the pressure rose back towards ambient. The positions of the shocks may be reflected in differences in the patterns of felled trees. At the limits of the devastated area, the temperature had dropped only 20% from the reservoir temperature because the entrained solids thermally buffered the flow (the dynamic and thermodynamic effects of the admixture of the surrounding atmosphere and the uprooted forest and soils into the flow are not considered). The density of the flow decreased with distance until, at the limits of the blast zone, 20–25 km from the volcano, the density became comparable with that of the surrounding (dirty) atmosphere and the flow became buoyant and ramped up into the atmosphere. According to the model, the mass flux per unit area at the source was 0.6 × 104 g s−1 cm−2 and the energy flux per unit area was 2.5 MW cm−2. From the measured total ejected mass, 0.25 × 1015 g, the total energy released during the eruption was 1024 erg or 24 megatons. The model, triggering of the eruption and the transition from unsteady to steady flow, and applications to eyewitness observations and atmospheric effects are discussed in ref. 1.