Seismic activity related to the 2000 eruption of the Hekla volcano, Iceland

The 2000 Hekla eruption took place from February 26 to March 8. Its seismic expressions were a swarm of numerous small earthquakes related to its onset, and low-frequency volcanic tremor that continued throughout the eruption. A swarm of small earthquakes was observed some 80 min before the onset of the eruption, and the size of the events increased with time. Low-frequency volcanic tremor, with a characteristic frequency band of 0.5–1.5 Hz and dominant spectral peak(s) at 0.7–0.9 Hz, became visible at 18:19 GMT on February 26, marking the onset of the eruption. The tremor amplitude rose quickly and was very high in the beginning of the eruption. However, it soon began to decrease after about an hour. In general, the seismic activity related to the 2000 Hekla eruption was very similar to what was observed in the previous eruption in 1991. Based on knowledge gained from seismicity and strain observations from 1991, this was the first time that a Hekla eruption was predicted.Editorial responsibility: J Stix     

The millennium eruption of Hekla in February 2000

The 18th historic eruption of Hekla started on 26 February, 2000. It was a short-lived but intense event, emitting basaltic andesitic (55.5 wt% SiO₂) pyroclastic fragments and lava. During the course of the eruption, monitoring was done by both instruments and direct observations, together providing unique insight into the current activity of Hekla. During the 12-day eruption, a total of 0.189 km³ DRE of magma was emitted. The eruptive fissure split into five segments. The segments at the highest altitude were active during the first hours, while the segments at lower altitude continued throughout the eruption. The eruption started in a highly explosive manner giving rise to a Subplinian eruptive column and consequent basaltic pyroclastic flows fed by column collapses. After the explosive phase reached its maximum, the eruption went through three more phases, namely fire-fountaining, Strombolian bursts and lava effusion. In this paper, we describe the course of events of the eruption of Hekla and the origin of its magma, and then show that the discharge rate can be linked to different style of eruptive activity, which are controlled by fissure geometry. We also show that the eruption phases observed at Hekla can be linked with inferred magma chamber overpressure prior to the eruption.     

Seismicity around the Hekla and Torfajökull volcanoes, Iceland, during a volcanically quiet period, 1991–1995

The volcano Hekla in south Iceland had its latest eruption in January–March 1991. The eruption was accompanied and followed by considerable seismic activity. This study examines the seismicity in the Hekla region (63°42′–64°18′N, 18°30′–20°12′W) during a period when the high activity related to the eruption had ceased, from July 1991 to October 1995. The aim is to define the level of the normal background seismicity of the area that can be compared to the eruption-related activity. The Hekla Volcano proper was generally aseismic during the study period. The most prominent earthquake cluster is in the neighbouring Torfajökull Volcano. The epicentres are concentrated in the western part of the caldera and west of it. The hypocentres are located at all depths from the surface down to 14?km, with highest activity at 5–12?km. Inside this cluster, in the northwest part of the caldera, is a spherical volume void of earthquakes, approximately 4?km in diameter and centred at 8?km depth. This is interpreted as a cooling magma body. Small, low-frequency events of volcanic origin were occasionally recorded at Torfajökull. This activity has mainly occurred in swarms and was most abundant during the first year of the study period, presumably reflecting some kind of connection to the 1991 Hekla eruption. Our study area also includes the easternmost section of the South Iceland seismic zone, a transform zone characterized by bookshelf faulting on transverse faults. Two lineaments of epicentres were identified, roughly corresponding to mapped faults of the South Iceland seismic zone. The hypocentres are relatively deep, mainly at 6–12?km, matching the general trend of hypocentral depth increasing toward the east. The seismicity is highest in the area of the mapped faults. However, the epicentres extend beyond them and indicate greater width of the South Iceland seismic zone, or 20–30?km rather than approximately 10?km as indicated by the length of the surface faults. The seismicity in the volcanic systems of Hekla and Vatnafjöll shows some characteristics of the South Iceland seismic zone. Epicentres are concentrated into two N–S lineaments, one of which coincides with the location of the 1987 Vatnafjöll earthquake (M_w=5.9), a strike-slip event on a N- to S-trending fault. The hypocentres of the Hekla–Vatnafjöll events are mainly at 8–13?km depth, which indicates a continuation of the depth trend of the earthquakes of the South Iceland seismic zone. The events located at Hekla proper and immediately north of it are all of low-frequency character, which can be held as an indication of volcanic origin. On the other hand, they show clear S arrivals at observing stations like normal high-frequency tectonic earthquakes.     

Doppler radar sounding of volcanic eruption dynamics at Mount Etna

Besides their common use in atmospheric studies, Doppler radars are promising tools for the active remote sensing of volcanic eruptions but were little applied to this field. We present the observations made with a mid-power UHF Doppler radar (Voldorad) during a 7-h Strombolian eruption at the SE crater of Mount Etna on 11–12 October 1998. Main characteristics of radar echoes are retrieved from analysis of Doppler spectra recorded in the two range gates on either side of the jet axis. From the geometry of the sounding, the contribution of uprising and falling ejecta to each Doppler spectrum can be discriminated. The temporal evolution of total power backscattered by uprising targets is quite similar to the temporal evolution of the volcanic tremor and closely reproduces the overall evolution of the eruption before, during and after its paroxysm. Moreover, during the sharp decrease of eruptive activity following the paroxysm, detailed analysis of video (from camera recording), radar and seismic measurements reveals that radar and video signals start to decrease simultaneously, approximately 2.5 min after the tremor decline. This delay is interpreted as the ascent time through a magma conduit of large gas slugs from a shallow source roughly estimated at about 500 m beneath the SE crater. Detailed analysis of eruptive processes has been also made with Voldorad operating in a high sampling rate mode. Signature of individual outburst is clearly identified on the half part of Doppler spectra corresponding to rising ejecta: temporal variations of the backscattered power exhibit quasi periodic undulations, whereas the maximum velocity measured on each spectrum displays a sharp peak at the onset of each outburst followed by a slow decay with time. Periodicity of power variations (between 3.8 and 5.5 s) is in agreement with the occurrence of explosions visually observed at the SE vent. Maximum vertical velocities of over 160 m s^–1 were measured during the paraoxysmal stage and the renewed activity. Finally, by using a simplified model simulating the radar echoes characteristics, we show that when Voldorad is operating in high sampling rate mode, the power and maximum velocity variations are directly related to the difference in size and velocity of particles crossing the antenna beam.Editorial responsibility: A. Woods     

New physical characterization of the Fontana Lapilli basaltic Plinian eruption,Nicaragua

The Fontana Lapilli deposit was erupted in the late Pleistocene from a vent, or multiple vents, located near Masaya volcano (Nicaragua) and is the product of one of the largest basaltic Plinian eruptions studied so far. This eruption evolved from an initial sequence of fluctuating fountain-like events and moderately explosive pulses to a sustained Plinian episode depositing fall beds of highly vesicular basaltic-andesite scoria (SiO₂ > 53 wt%). Samples show unimodal grain size distribution and a moderate sorting that are uniform in time. The juvenile component predominates (> 96 wt%) and consists of vesicular clasts with both sub-angular and fluidal, elongated shapes. We obtain a maximum plume height of 32 km and an associated mass eruption rate of 1.4 × 10⁸ kg s⁻¹ for the Plinian phase. Estimates of erupted volume are strongly sensitive to the technique used for the calculation and to the distribution of field data. Our best estimate for the erupted volume of the majority of the climactic Plinian phase is between 2.9 and 3.8 km³ and was obtained by applying a power-law fitting technique with different integration limits. The estimated eruption duration varies between 4 and 6 h. Marine-core data confirm that the tephra thinning is better fitted by a power-law than by an exponential trend.     

Sulfur,chlorine, and fluorine degassing and atmospheric loading by the 1783–1784 AD Laki (Skaftár Fires) eruption in Iceland

 The 1783–1784 Laki tholeiitic basalt fissure eruption in Iceland was one of the greatest atmospheric pollution events of the past 250 years, with widespread effects in the northern hemisphere. The degassing history and volatile budget of this event are determined by measurements of pre-eruption and residual contents of sulfur, chlorine, and fluorine in the products of all phases of the eruption. In fissure eruptions such as Laki, degassing occurs in two stages: by explosive activity or lava fountaining at the vents, and from the lava as it flows away from the vents. Using the measured sulfur concentrations in glass inclusions in phenocrysts and in groundmass glasses of quenched eruption products, we calculate that the total accumulative atmospheric mass loading of sulfur dioxide was 122 Mt over a period of 8 months. This volatile release is sufficient to have generated ∼250 Mt of H₂SO₄ aerosols, an amount which agrees with an independent estimate of the Laki aerosol yield based on atmospheric turbidity measurements. Most of this volatile mass (∼60 wt.%) was released during the first 1.5 months of activity. The measured chlorine and fluorine concentrations in the samples indicate that the atmospheric loading of hydrochloric acid and hydrofluoric acid was ∼7.0 and 15.0 Mt, respectively. Furthermore, ∼75% of the volatile mass dissolved by the Laki magma was released at the vents and carried by eruption columns to altitudes between 6 and 13 km. The high degree of degassing at the vents is attributed to development of a separated two-phase flow in the upper magma conduit, and implies that high-discharge basaltic eruptions such as Laki are able to loft huge quantities of gas to altitudes where the resulting aerosols can reside for months or even 1–2 years. The atmospheric volatile contribution due to subsequent degassing of the Laki lava flow is only 18 wt.% of the total dissolved in the magma, and these emissions were confined to the lowest regions of the troposphere and therefore important only over Iceland. This study indicates that determination of the amount of sulfur degassed from the Laki magma batch by measurements of sulfur in the volcanic products (the petrologic method) yields a result which is sufficient to account for the mass of aerosols estimated by other methods. Received: 30 May 1995 / Accepted: 19 April 1996     

Earthquake swarms preceding the 2000 eruption of Miyakejima volcano, Japan

The first sign of magma accumulating beneath Miyakejima, an island volcano in the northern Izu islands, Japan, came at around 18:00 on 26 June 2000, when a swarm of earthquakes was detected by a volcano seismic network on the island. Earthquakes occurred initially beneath the southwest flank near the summit and gradually migrated west of the island, where a submarine eruption occurred the next morning. Earthquakes then migrated further to the northwest between Miyakejima and Kozushima, another volcanic island and developed to the most intense earthquake swarm ever observed in and around Japanese archipelago. To better image how the initial magma intrusion occurred, we relocated hypocenters by using a station-correction method and a double-difference method. The relocated epicenters are generally concentrated near the upper bound of dyke intrusions inferred from geodetic studies throughout the initial stages of the 2000 eruption at Miyakejima from 26 to 27 June 2000. As for seismic activity westward off Miyakejima in the morning on 27 June, hypocenters from both a nationwide seismic network that were relocated by the double-difference method, and those from the volcano seismic network relocated by the station-correction method, formed a very shallow cluster that ascended slowly with time as it propagated northwestward from Miyakejima. This suggests that the dykes have both a radial and upward component of movement.Editorial responsibility: S. Nakada, T. Druitt     

Observations of volcanic tremor during the January–February 2005 eruption of Mt. Veniaminof,Alaska

Mt. Veniaminof, Alaska Peninsula, is a stratovolcano with a summit ice-filled caldera containing a small intracaldera cone and active vent. From January 2 to February 21, 2005, Mt. Veniaminof erupted. The eruption was characterized by numerous small ash emissions (VEI 0 to 1) and accompanied by low-frequency earthquake activity and volcanic tremor. We have performed spectral analyses of the seismic signals in order to characterize them and to constrain their source. Continuous tremor has durations of minutes to hours with dominant energy in the band 0.5–4.0 Hz, and spectra characterized by narrow peaks either irregularly (non-harmonic tremor) or regularly spaced (harmonic tremor). The spectra of non-harmonic tremor resemble those of low-frequency events recorded simultaneously with surface ash explosions, suggesting that the source mechanisms might be similar or related. We propose that non-harmonic tremor at Mt. Veniaminof results from the coalescence of gas bubbles while low-frequency events are related to the disruption of large gas pockets within the conduit. Harmonic tremor, characterized by regular and quasi-sinusoidal waveforms, has duration of hours. Spectra containing up to five harmonics suggest the presence of a resonating source volume that vibrates in a longitudinal acoustic mode. An interesting feature of harmonic tremor is that frequency is observed to change over time; spectral lines move towards higher or lower values while the harmonic nature of the spectra is maintained. Factors controlling the variable characteristics of harmonic tremor include changes in acoustic velocity at the source and variations of the effective size of the resonator.     

Geology of the Heimaey volcanic centre, south Iceland: early evolution of a central volcano in a propagating rift?

Heimaey is the southernmost and also the youngest of nine volcanic centres in the southward-propagating Eastern Volcanic Zone, Iceland. The island of Heimaey belongs to the Vestmannaeyjar volcanic system (850 km²) and is situated 10 km off the south coast of Iceland. Although Heimaey probably started to form during the Upper Pleistocene all the exposed subaerial volcanics (10 monogenetic vents covering an area of 13.4 km²) are of Holocene age. Heimaey is composed of roughly equal amounts of tuff/tuff-breccias and lavas as most eruptions involve both a phreatomagmatic and an effusive phase. The compositions of the extrusives are predominantly alkali basalts belonging to the sodic series. Repeated eruptions on Heimaey, and the occurrence of slightly more evolved rocks (i.e. hawaiite approaching mugearite), might indicate that the island is in an early stage of forming a central volcano in the Vestmannaeyjar system. This is further substantiated by the development of a magma chamber at 10–20 km depth during the most recent eruption in 1973 and by the fact that the average volume of material produced in a single eruption on Heimaey is 0.32 km³ (dense rock equivalent), which is twice the value reported for the Vestmannaeyjar system as a whole. We find no support for the previously postulated episodic behaviour of the volcanism in the Vestmannaeyjar system. However, the oldest units exposed above sea level, i.e. the Norðurklettar ridge, probably formed over a 500-year interval during the deglaciation of southern Iceland. The absence of equilibrium phenocryst assemblages in the Heimaey lavas suggests that magma rose quickly from depth, without long-time ponding in shallow-seated crustal magma chambers. Eruptions on Heimaey have occurred along two main lineaments (N45°E and N65°E), which indicate that it is seismic events associated with the southward propagation of the Eastern Volcanic Zone that open pathways for the magma to reach the surface. Continuing southward propagation of the Eastern Volcanic Zone suggests that the frequency of volcanic eruptions in the Vestmannaeyjar system might increase with time, and that Heimaey may develop into a central volcano like the mature volcanic centres situated on the Icelandic mainland.     

Magma plumbing system of the 2000 eruption of Miyakejima Volcano, Japan

During the 2000 eruption at Miyakejima Volcano, two magmas with different compositions erupted successively from different craters. Magma erupted as spatter from the submarine craters on 27 June is aphyric basaltic andesite (<5 vol% phenocrysts, 51.4–52.2 wt% SiO₂), whereas magma issued as volcanic bombs from the summit caldera on 18 August is plagioclase-phyric basalt (20 vol% phenocrysts, 50.8–51.3 wt% SiO₂). The submarine spatter contains two types of crystal-clots, A-type and A-type (andesitic type). The phenocryst assemblages (plagioclase, pyroxenes and magnetite) and compositions of clinopyroxene in these clots are nearly the same, but only A-type clots contain Ca-poor plagioclase (An < 70). We consider that the A-type clots could have crystallized from a more differentiated andesitic magma than the A-type clots, because FeO*/MgO is not strongly influenced during shallow andesitic differentiation. The summit bombs contain only B-type (basaltic type) crystal-clots of Ca-rich plagioclase, olivine and clinopyroxene. The A-type and B-type clots have often coexisted in Miyakejima lavas of the period 1469–1983, suggesting that the magma storage system consists of independent batches of andesitic and basaltic magmas. According to the temporal variations of mineral compositions in crystal-clots, the andesitic magma became less evolved, and the basaltic magma more evolved, over the past 500 years. We conclude that gradually differentiating basaltic magma has been repeatedly injected into the shallower andesitic magma over this period, causing the andesitic magma to become less evolved with time. The mineral chemistries in crystal-clots of the submarine spatter and 18 August summit bombs of the 2000 eruption fall on the evolution trends of the A-type and B-type clots respectively, suggesting that the shallow andesitic and deeper basaltic magmas existing since 1469 had successively erupted from different craters. The 2000 summit collapse occurred due to drainage of the andesitic magma from the shallower chamber; as the collapse occurred, it may have caused disruption of crustal cumulates which then contaminated the ascending, deeper basalt. Thus, porphyritic basaltic magma could erupt alone without mixing with the andesitic magma from the summit caldera. The historical magma plumbing system of Miyakejima was probably destroyed during the 2000 eruption, and a new one may now form.Editorial responsibility: S Nakada, T Druitt     

Nature and significance of small volume fall deposits at composite volcanoes: Insights from the October 14, 1974 Fuego eruption,Guatemala

The first of four successive pulses of the 1974 explosive eruption of Fuego volcano, Guatemala, produced a small volume (∼0.02 km³ DRE) basaltic sub-plinian tephra fall and flow deposit. Samples collected within 48 h after deposition over much of the dispersal area (7–80 km from the volcano) have been size analyzed down to 8 φ (4 μm). Tephra along the dispersal axis were all well-sorted (σ _φ = 0.25–1.00), and sorting increased whereas thickness and median grain size decreased systematically downwind. Skewness varied from slightly positive near the vent to slightly negative in distal regions and is consistent with decoupling between coarse ejecta falling off the rising eruption column and fine ash falling off the windblown volcanic cloud advecting at the final level of rise. Less dense, vesicular coarse particles form a log normal sub-population when separated from the smaller (Md_φ < 3φ or < 0.125 mm), denser shard and crystal sub-population. A unimodal, relatively coarse (Md_φ = 0.58φ or 0.7 mm σ _φ = 1.2) initial grain size population is estimated for the whole (fall and flow) deposit. Only a small part of the fine-grained, thin 1974 Fuego tephra deposit has survived erosion to the present day. The initial October 14 pulse, with an estimated column height of 15 km above sea level, was a primary cause of a detectable perturbation in the northern hemisphere stratospheric aerosol layer in late 1974 to early 1975. Such small, sulfur-rich, explosive eruptions may substantially contribute to the overall stratospheric sulfur budget, yet leave only transient deposits, which have little chance of survival even in the recent geologic record. The fraction of finest particles (Md_φ = 4–8φ or 4–63 μm) in the Fuego tephra makes up a separate but minor size mode in the size distribution of samples around the margin of the deposit. A previously undocumented bimodal–unimodal–bimodal change in grain size distribution across the dispersal axis at 20 km downwind from the vent is best accounted for as the result of fallout dispersal of ash from a higher subplinian column and a lower “co-pf” cloud resulting from pyroclastic flows. In addition, there is a degree of asymmetry in the documented grain-size fallout pattern which is attributed to vertically veering wind direction and changing windspeeds, especially across the tropopause. The distribution of fine particles (<8 μm diameter) in the tephra deposit is asymmetrical, mainly along the N edge, with a small enrichment along the S edge. This pattern has hazard significance.     

Petrochemistry of the Gjálp-1996 subglacial eruption, Vatnajökull, SE Iceland

In October 1996 a subglacial fissure to the north of the Grimsvötn caldera in W-Vatnajökull produced about 0.4 km³ of Fe-rich basaltic andesite–icelandite—in an area characterized mostly by tholeiitic basalt. In this paper the chemical composition of volcanic systems in the region is discussed with the help of six new analyses and others from the literature, and a tentative model for their evolution is proposed, in which magma produced by the partial melting of a two-component mantle mixes with hydrous, silicic melt in the crust. The Vatnajökull 1996 magma belongs to a separate volcanic system, intermediate between Bardarbunga and Grimsvötn.     

Submarine flank eruption preceding caldera subsidence during the 2000 eruption of Miyakejima Volcano, Japan

During the early part of a seismic swarm preceding eruption and caldera formation at Miyakejima Volcano, discoloured sea surfaces were observed 1.5 km off the western coast of Miyakejima on 27 June 2000. A later survey of the area using a multi-beam side scan sonar and a remotely operated small submarine revealed four craters of 20–30 m diameter aligned east-west in a 100×10–30 m area on the seafloor, with hot water at 140°C being released from one of the centres. Each crater consists of submarine spatter overlain in part by scoria lapilli. Dredged spatter from the craters was fresh, and there was no evidence of activity of marine organisms on the spatter surface, indicating that the discoloured sea surface resulted from magmatic eruption on the seafloor. This eruption occurred when a westward-propagating seismic swarm, initiated beneath Miyakejimas summit, passed through the area. Finding new magma on the seafloor demonstrates that this seismic swarm was associated with intruding magma, moving outward from beneath Miyakejima. Submarine spatter shows flattened shapes with a brittle crust formed by cooling in water, and its composition is aphyric andesite of 54 wt% SiO₂. The spatter is similar in whole rock and mineral composition to spatter erupted in 1983. However, the wide range of Cl in melt inclusions in plagioclase of the 27 June submarine spatter shows that it is not simply a remnant of the 1983 magma, which has only high Cl melt inclusions in plagioclase. The mixed character of melt inclusions suggests involvement of a magma with low Cl melt inclusions. The magma erupted explosively on 18 August from Miyakejimas summit, considered as the second juvenile magma in this eruption, contains low Cl melt inclusions in plagioclase. Based on these observations and the eruption sequence, we present the following model: (1) A shallow magma chamber was filled with a remnant of 1983 magma that had evolved to a composition of 54–55 wt% SiO₂. (2) Injection of the 18 August magma into this chamber generated a mixed magma having a wide range of Cl in melt inclusions contained plagioclase. The magma mixing might have occurred shortly before the submarine eruption and could have been a trigger for the initiation of the removal of magma from the chamber as an extensive dyke, which eventually led to caldera subsidence.Editorial responsibility: S Nakada, T Druitt     

Tephra-fall deposits from the 1992 eruption of Crater Peak, Alaska: implications of clast textures for eruptive processes

 The 1992 eruption of Crater Peak, Mount Spurr, Alaska, involved three subplinian tephra-producing events of similar volume and duration. The tephra consists of two dense juvenile clast types that are identified by color, one tan and one gray, of similar chemistry, mineral assemblage, and glass composition. In two of the eruptive events, the clast types are strongly stratified with tan clasts dominating the basal two thirds of the deposits and gray clasts the upper one third. Tan clasts have average densities between 1.5 and 1.7 g/cc and vesicularities (phenocryst free) of approximately 42%. Gray clasts have average densities between 2.1 and 2.3 g/cc, and vesicularities of approximately 20%; both contain abundant microlites. Average maximum plagioclase microlite lengths (13–15 μm) in gray clasts in the upper layer are similar regardless of eruptive event (and therefore the repose time between them) and are larger than average maximum plagioclase microlite lengths (9–11 μm) in the tan clasts in the lower layer. This suggests that microlite growth is a response to eruptive processes and not to magma reservoir heterogeneity or dynamics. Furthermore, we suggest that the low vesicularities of the clasts are due to syneruptive magmatic degassing resulting in microlitic growth prior to fragmentation and not to quenching of clasts by external groundwater. Received: 5 September 1997 / Accepted: 1 February 1998     

Integrated satellite observations of the 2001 eruption of Mt. Cleveland, Alaska

Satellite data were the primary source of information for the eruption of Mt. Cleveland, Alaska on 19 February, and 11 and 19 March 2001. Multiple data sets were used pre-, syn- and post-eruption to mitigate the hazard and determine an eruption chronology. The 19 February eruption was the largest of the three, resulting in a volcanic cloud that formed an arc over 1000 km long, moved to the NE across Alaska and was tracked using satellite data over more than a 50-h period. The volcanic cloud was “concurrently” detected on the GOES, AVHRR and MODIS data at various times and their respective signals compared. All three sensors detected a cloud that had a very similar shape and position but there were differences in their areal extent and internal structural detail. GOES data showed the largest volcanic cloud in terms of area, probably due to its oblique geometry. MODIS bands 31 and 32, which are comparable to GOES and AVHRR thermal infrared wavelengths, were the least effective single channels at detecting the volcanic cloud of those investigated (MODIS bands 28, 29, 31 and 32). MODIS bands 28 and 29 detected the largest volcanic clouds that could easily be distinguished from weather clouds. Of the split-window data, MODIS bands 29 minus band 32 detected the largest cloud, but the band 31 minus band 32 data showed the volcanic cloud with the most internal structural detail. The Puff tracking model accurately tracked the movement, and predicted the extent and shape of this complex cloud even into areas beyond satellite detection. Numerous thermal anomalies were also observed during the eruption on the twice-daily AVHRR data and the high spatial-resolution Landsat data. The high-resolution Radarsat data showed that the AVHRR thermal anomalies were due to lava and debris flow features and a newly formed fan along the west coast of the island. Field observations and images from a hand-held Forward Looking Infrared Radiometer (FLIR) showed that the flow features were ′a′a lava, debris flows and a warm debris fan along the west coast. Real-time satellite data were the primary tool used to monitor the eruption, track changes and to mitigate hazards. High-resolution data, even though coverage is infrequent, were critical in helping to identify volcanic processes and to compile an eruption chronology.     

Implication of the temporal sulphur isotopic variation during the 2000 eruption of Miyakejima Volcano, Japan

Abstract   The variation of sulphur isotopic composition during the 2000 eruption of the Miyakejima Volcano was examined in order to monitor the temporal change of the volcanic activity. The δ³⁴S values of water-soluble sulphate leached from volcanic ash effused during intermittent eruptions from July to September 2005 range from +5 to +11‰ with a fluctuation of ca 3‰ within a single eruption. The δ³⁴S value of sulphuric acid mist collected with 'Cu-metal trap' placed on the flank of the volcano from December 2000 to January 2001 is +6.2‰. These sulphur isotopic compositions of sulphate, which were isotopically equilibrated in the subvolcanic hydrothermal system, indicate that the temperature of the hydrothermal system beneath the caldera increased after the period of intermittent phreatic and phreatomagmatic eruptions. Then, the δ³⁴S value of sulphuric acid trapped from January to March 2001 was +9.0‰ and the δ³⁴S value of water-soluble sulphate on volcanic ash emitted with minor eruption in May 2001 was +11.0‰, suggesting a decrease in temperature of the subvolcanic hydrothermal system.     

Under the surface: Pressure-induced planetary-scale waves,volcanic lightning,and gaseous clouds caused by the submarine eruption of Hunga Tonga-Hunga Ha'apai volcano

We present a narrative of the eruptive events culminating in the cataclysmic January 15, 2022 eruption of Hunga Tonga-Hunga Ha'apai Volcano by synthesizing diverse preliminary seismic, volcanological, sound wave, and lightning data available within the first few weeks after the eruption occurred. The first hour of eruptive activity produced fast-propagating tsunami waves, long-period seismic waves, loud audible sound waves, infrasonic waves, exceptionally intense volcanic lightning and an unsteady volcanic plume that transiently reached—at 58 ?km—the Earth's mesosphere. Energetic seismic signals were recorded worldwide and the globally stacked seismogram showed episodic seismic events within the most intense periods of phreatoplinian activity, and they correlated well with the infrasound pressure waveform recorded in Fiji. Gravity wave signals were strong enough to be observed over the entire planet in just the first few hours, with some circling the Earth multiple times subsequently. These large-amplitude, long-wavelength atmospheric disturbances come from the Earth's atmosphere being forced by the magmatic mixture of tephra, melt and gasses emitted by the unsteady but quasi-continuous eruption from 0402±1–1800 UTC on January 15, 2022. Atmospheric forcing lasted much longer than rupturing from large earthquakes recorded on modern instruments, producing a type of shock wave that originated from the interaction between compressed air and ambient (wavy) sea surface. This scenario differs from conventional ideas of earthquake slip, landslides, or caldera collapse-generated tsunami waves because of the enormous (～1000x) volumetric change due to the supercritical nature of volatiles associated with the hot, volatile-rich phreatoplinian plume. The time series of plume altitude can be translated to volumetric discharge and mass flow rate. For an eruption duration of ～12 ?h, the eruptive volume and mass are estimated at 1.9 ?km³ and ～2 900 ?Tg, respectively, corresponding to a VEI of 5–6 for this event. The high frequency and intensity of lightning was enhanced by the production of fine ash due to magma—seawater interaction with concomitant high charge per unit mass and the high pre-eruptive concentration of dissolved volatiles. Analysis of lightning flash frequencies provides a rapid metric for plume activity and eruption magnitude. Many aspects of this eruption await further investigation by multidisciplinary teams. It represents a unique opportunity for fundamental research regarding the complex, non-linear behavior of high energetic volcanic eruptions and attendant phenomena, with critical implications for hazard mitigation, volcano forecasting, and first-response efforts in future disasters.     

Characterization of shape and terminal velocity of tephra particles erupted during the 2002 eruption of Etna volcano,Italy

In this paper, we present a complete morphological characterization of the ash particles erupted on 18 December 2002 from Etna volcano, Italy. The work is based on the acquisition and processing of bidimensional digital images carried out by scanning electron microscopy (SEM) to obtain shape parameters by image analysis. We measure aspect ratio (AR), form factor (FF), compactness (CC), and rectangularity (RT) of 2065 ash particles with size between 0.026 and 1.122 mm. We evaluate the variation of these parameters as a function of the grain-size. Ash particles with a diameter of <0.125 mm vary from mostly equant to very equant, ash particles between 0.125 and 0.250 mm have an intermediate shape, and particles with diameters >0.250 mm are subelongate. We find that, on average, particles with a diameter of <0.250 mm are subrounded, particles between 0.250 and 0.50 mm are subangular, and particles >0.50 mm are angular. Using this morphological analysis and an empirical relation between the drag coefficient (C _D) and Reynolds number (R _e) of Wilson and Huang (Earth Planet Sci Lett 44:311–324, 1979), we calculate the terminal settling velocities (V _WH). The comparisons between these velocities and those calculated with the formula of Kunii and Levenspiel (Fluidization engineering. Wiley, New York, pp 97, 1969) (V _KL), which considers ash particles as spheres, show that V _KL are in average 1.28 greater than V _WH. Hence, we quantify the systematic error on the spatial distribution of the mass computed around the volcano carried out by tephra dispersal models when the particles are assumed to be spherical.     

Preliminary sensitivity study of eruption source parameters for operational volcanic ash cloud transport and dispersion models — A case study of the August 1992 eruption of the Crater Peak vent, Mount Spurr, Alaska

Ash clouds are one of the major hazards that result from volcanic eruptions. Once an eruption is reported, volcanic ash transport and dispersion (VATD) models are used to forecast the location of the ash cloud. These models require source parameters to describe the ash column for initialization. These parameters include: eruption cloud height and vertical distribution, particle size distribution, and start and end time of the eruption. Further, if downwind concentrations are needed, the eruption mass rate and/or volume of ash need to be known. Upon notification of an eruption, few constraints are typically available on many of these source parameters. Recently, scientists have defined classes of eruption types, each with a set of pre-defined eruption source parameters (ESP). We analyze the August 18, 1992 eruption of the Crater Peak vent at Mount Spurr, Alaska, which is the example case for the Medium Silicic eruption type. We have evaluated the sensitivity of two of the ESP – the grain size distribution (GSD) and the vertical distribution of ash – on the modeled ash cloud. HYSPLIT and Puff VATD models are used to simulate the ash clouds from the different sets of source parameters. We use satellite data, processed through the reverse absorption method, as reference for computing statistics that describe the modeled-to-observed comparison. With the grain size distribution, the three options chosen, (1) an estimated distribution based on past eruption studies, (2) a distribution with finer particles and (3) the National Oceanic and Atmospheric Administration HYSPLIT GSD, have little effect on the modeled ash cloud. For the initial vertical distribution, both linear (uniform concentration throughout the vertical column) and umbrella shapes were chosen. For HYSPLIT, the defined umbrella distribution (no ash below the umbrella), apparently underestimates the lower altitude portions of the ash cloud and as a result has a worse agreement with the satellite detected ash cloud compared to that with the linear vertical distribution for this particular eruption. The Puff model, with a Poisson function to represent the umbrella cloud, gave similar results as for a linear distribution, both having reasonable agreement with the satellite detected cloud. Further sensitivity studies of this eruption, as well as studies using the other source parameters, are needed.     

Mechanism of explosive eruption revealed by geophysical observations at the Sakurajima,Suwanosejima and Semeru volcanoes

A common sequence of phenomena associated with volcanic explosions is extracted based on seismic and ground deformation observations at 3 active volcanoes in Japan and Indonesia. Macroscopic inflation-related ground deformations are detected prior to individual explosions, while deflations are observed during eruptions. Precursory inflation occurs 5 min to several hours before eruption at the Sakurajima volcano, but just 1–2 min at Suwanosejima and 3–30 min at the Semeru volcano. The sequence commences with minor contraction, which is detected by extensometers 1.5 min before eruption at Sakurajima, as a dilatant first motion of the explosion earthquakes 0.2–0.3 s before surface explosions at Suwanosejima, and as downward tilt 4–5 s prior to eruption at the Semeru volcano. The sequence is detected for explosive eruptions with > 0.1 μrad tilt change at Sakurajima, 90% at Suwanosejima and 75% at Semeru volcanoes. It is inferred that the minor contraction is caused by a volume and pressure decrease due to the release of gas from a pocket at the top of the conduit as the gas pressure exceeds the strength of the confining plug. The subsequent violent expansion may be triggered by sudden outgassing of the water-saturated magma induced by the decrease in confining pressure.