Chemical characteristics of the crater lakes of Popocatetetl, El Chichon, and Nevado de Toluca volcanoes, Mexico

Three crater lakes from Mexican volcanoes were sampled and analyzed at various dates to determine their chemical characteristics. Strong differences were observed in the chemistry among the three lakes: Nevado de Toluca, considered as dormant, El Chichón at a post-eruptive stage, and Popocatépetl at a pre-eruptive stage. Not surprisingly, no influence of volcanic activity was found at the Nevado de Toluca volcano, while the other volcanoes showed a correlation between the changing level of activity and the evolution of chemical trends. Low pHs (<3.0) were measured in the water from the active volcanoes, while a pH of 5.6 was measured at the Nevado de Toluca Sun lake. Changes with time were observed at Popocatépetl and El Chichón. Concentrations of volcanic-gas derived species like Cl⁻, SO₄²⁻ and F⁻ decreased irregularly at El Chichón from 1983 until 1997. Major cations concentrations also diminished at El Chichón. A 100% increase in the SO₄²⁻ content was measured at Popocatépetl between 1985 and 1994. An increase in the Mg/Cl ratio between 1992 (Mg/Cl=0.085) and 1994 (Mg/Cl=0.177) was observed at Popocatépetl, before the disappearance of the crater lake in 1994. It is concluded that chemical analysis of crater lakes may provide a useful additional tool for active-volcano monitoring.     

Precursory seismicity of the 1994 eruption of Popocatépetl Volcano,Central Mexico

Popocatépetl Volcano is located in the central Mexican Volcanic Belt, within a densely populated region inhabited by over 20 million people. The eruptive history of this volcano indicates that it is capable of producing a wide range of eruptions, including Plinian events. After nearly 70 years of quiescence, Popocatépetl reawakened in December 21, 1994. The eruptive activity has continued up until the date of this submission and has been characterized by a succession of lava dome growth-and-destruction episodes, similar to events that have apparently been typical for Popocatépetl since the fourteenth century. In this regime, the episodes of effusive and moderately explosive activity alternate with long periods of almost total quiescence. In this paper we analyze five years of volcano-tectonic seismicity preceding the initial eruption of the current episode. The evolution of the V-T seismicity shows four distinct stages, which we interpret in terms of the internal processes which precede an eruption after a long period of quiescence. The thermal effects of a magma intrusion at depth, the fracturing related to the slow development of magma-related fluid pathways, the concentration of stress causing a protracted acceleration of this process, and a final relaxation or redistribution of the stress shortly before the initial eruption are reflected in the rates of V-T seismic energy release. A hindsight analysis of this activity shows that the acceleration of the seismicity in the third stage asymptotically forecast the time of the eruption. The total seismic energy release needed to produce an eruption after a long period of quiescence is related to the volume of rock that must be fractured so imposing a characteristic threshold limit for polygenetic volcanoes, limit that was reached by Popocatépetl before the eruption.     

Englacial tephrostratigraphy of Erebus volcano,Antarctica

A tephrostratigraphy for Erebus volcano is presented, including tephra composition, stratigraphy, and eruption mechanism. Tephra from Erebus were collected from glacial ice and firn. Scanning electron microscope images of the ash morphologies help determine their eruption mechanisms The tephra resulted mainly from phreatomagmatic eruptions with fewer from Strombolian eruptions. Tephra having mixed phreatomagmatic–Strombolian origins are common. Two tephra deposited on the East Antarctic ice sheet, ~ 200 km from Erebus, resulted from Plinian and phreatomagmatic eruptions. Glass droplets in some tephra indicate that these shards were produced in both phreatomagmatic and Strombolian eruptions. A budding ash morphology results from small spheres quenched during the process of hydrodynamically splitting off from a parent melt globule. Clustered and rare single xenocrystic analcime crystals, undifferentiated zeolites, and clay are likely accidental clasts entrained from a hydrothermal system present prior to eruption. The phonolite compositions of glass shards confirm Erebus volcano as the eruptive source. The glasses show subtle trends in composition, which correlate with stratigraphic position. Trace element analyses of bulk tephra samples show slight differences that reflect varying feldspar contents.     

Phreatomagmatic eruptions associated with the caldera collapse during the Miyakejima 2000 eruption,Japan

The 2000 AD eruption of Miyakejima was characterized by a series of phreatomagmatic eruptions from the subsiding caldera. Six major eruptive events occurred, and they can be divided into the first and second periods separated by a 25-day hiatus. The phreatomagmatic eruptions produced a total of ~ 2 × 10¹⁰ kg of tephra, which mainly comprised fine-grained volcanic ash. The tephra layers could be divided into six fall units corresponding to the six major eruptive events.     

Late Holocene phases of dome growth and Plinian activity at Guagua Pichincha volcano (Ecuador)

Since the eruption which affected Quito in AD 1660, Guagua Pichincha has been considered a hazardous volcano. Based on field studies and twenty ¹⁴C dates, this paper discusses the eruptive activity of this volcano, especially that of the last 2000 years. Three major Plinian eruptions with substantial pumice discharge occurred in the 1st century, the 10th century, and in AD 1660. The ages of organic paleosols and charcoal from block-and-ash flow and fallout deposits indicate that these eruptions occurred near the end of 100 to 200 year-long cycles of discontinuous activity which was comprised of dome growth episodes and minor pumice fallouts. The first cycle took place from ~ AD 1 to 140. The second one developed during the 9th and 10th centuries, lasted 150–180 yr, and included the largest Plinian event, with a VEI of 5. The third, historic cycle, about 200 yr in duration, includes pyroclastic episodes around AD 1450 and AD 1500, explosive activity between AD 1566 and AD 1582, possible precursors of the 1660 eruption in the early decades of the 17th century, and finally the 1660 eruption (VEI 4). A fourth event probably occurred around AD 500, but its authenticity requires confirmation. The Plinian events occurred at the end of these cycles which were separated by repose periods of at least 300 yr. Older volcanic activity of similar type occurred between ~ 4000 and ~ 3000 yr BP.     

Remarkable geochemical changes and degassing at Voui crater lake,Ambae volcano,Vanuatu

Ambae (also known as Aoba), is a 38 × 16 km² lozenge-shaped island volcano with a coastal population of around 10 000. At the summit of the volcano is lake Voui — one of the largest active crater lakes worldwide, with 40 × 10⁶ m³ of acidic water perched 1400 m a.s.l. After more than 300 years of dormancy, Ambae volcano reawakened with phreatic eruptions through Voui in 1995, and culminating in a series of surtseyan eruptions in 2005, followed by a rapid and spectacular colour change of the lake from light blue to red in 2006. Integrating lake water chemistry with new measurements of SO₂ emissions from the volcano during the 2005–2006 eruptive period helps to explain the unusual and spectacular volcanic activity of Ambae — initially, a degassed magma approached the lake bed and triggered the surtseyan eruption. Depressurization of the conduit facilitated ascent of volatile-rich magma from the deeper plumbing system. The construction of a cone during eruption and the high degassing destabilised the equilibrium of lake stratification leading to a limnic event and subsequently the spectacular colour change.     

Sulfur dioxide emissions from Popocatépetl volcano (Mexico): case study of a high-emission rate, passively degassing erupting volcano

Popocatépetl volcano in central Mexico has been erupting explosively and effusively for almost 4 years. SO₂ emission rates from this volcano have been the largest ever measured using a COSPEC. Pre-eruptive average SO₂ emission rates (2–3 kt/d) were similar to the emission rates measured during the first part of the eruption (up to August 1995) in contrast with the effusive–explosive periods (March 1996–January 1998) during which SO₂ emission rates were higher by a factor of four (9–13 kt/d). Based on a chronology of the eruption and the average SO₂ emission rates per period, the total SO₂ emissions (up to 1 January 1998) are estimated to be about 9 Mt, roughly half as much as the SO₂ emissions from Mount Pinatubo in a shorter period. Popocatépetl volcano is thus considered as a high-emission rate, passively degassing eruptive volcano. SO₂ emission rates and SO₂ emissions are used here to make a mass balance of the erupted magma and related gases. Identified excess SO₂ is explained in terms of continuous degassing of unerupted magma and magma mixing. Fluctuations in SO₂ emission rate may be a result of convection and crystallization in the chamber or the conduits, cleaning and sealing of the plumbing system, and/or SO₂ scrubbing by the hydrothermal system.     

Quantification of the 1998–1999 explosion sequence at Popocatépetl volcano,Mexico

The sequence of large Vulcanian explosions occurring at the andesitic Popocatépetl volcano, Mexico during November 1998 to April 1999 was studied. The size of 26 largest explosions was estimated from broadband seismic records at the distance of 4 km from the crater. The sequence began with the largest explosion (E = 2.6 × 10¹² J) occurring on 25 November at 08:05, and following largest daily explosions were characterized by gradual decrease in the energy. The energy of 20 large (E ≥ 10¹¹ J) explosions was distributed as Student's t-distribution with a geometrical mean Log E = 11.81 (J).     

Probabilistic hazard analysis of Citlaltépetl (Pico de Orizaba) Volcano, eastern Mexican Volcanic Belt

Citlaltépetl or Pico de Orizaba is the highest active volcano in the North American continent. Although Citlaltépetl is at present in repose, its eruptive history reveals repetitive explosive eruptions in the past. Its relatively low eruption rate has favored significant population growth in areas that may be affected by a potential eruptive activity. The need of some criteria for hazards assessment and land-use planning has motivated the use of statistical methods to estimate the time and space distribution of volcanic hazards around this volcano. The analysis of past activity, from late Pleistocene to historic times, and the extent of some well-identified deposits are used to calculate the recurrence probabilities of eruptions of various size during time periods useful for land-use planning.     

A statistical method linking geological and historical eruption time series for volcanic hazard estimations: Applications to active polygenetic volcanoes

The probabilistic analysis of volcanic eruption time series is an essential step for the assessment of volcanic hazard and risk. Such series describe complex processes involving different types of eruptions over different time scales. A statistical method linking geological and historical eruption time series is proposed for calculating the probabilities of future eruptions. The first step of the analysis is to characterize the eruptions by their magnitudes. As is the case in most natural phenomena, lower magnitude events are more frequent, and the behavior of the eruption series may be biased by such events. On the other hand, eruptive series are commonly studied using conventional statistics and treated as homogeneous Poisson processes. However, time-dependent series, or sequences including rare or extreme events, represented by very few data of large eruptions require special methods of analysis, such as the extreme-value theory applied to non-homogeneous Poisson processes. Here we propose a general methodology for analyzing such processes attempting to obtain better estimates of the volcanic hazard. This is done in three steps: Firstly, the historical eruptive series is complemented with the available geological eruption data. The linking of these series is done assuming an inverse relationship between the eruption magnitudes and the occurrence rate of each magnitude class. Secondly, we perform a Weibull analysis of the distribution of repose time between successive eruptions. Thirdly, the linked eruption series are analyzed as a non-homogeneous Poisson process with a generalized Pareto distribution as intensity function. As an application, the method is tested on the eruption series of five active polygenetic Mexican volcanoes: Colima, Citlaltépetl, Nevado de Toluca, Popocatépetl and El Chichón, to obtain hazard estimates.     

Seismic characterization of the fall 2007 eruptive sequence at Bezymianny Volcano,Russia

We examine an eruptive sequence in late 2007 at Bezymianny Volcano to characterize the magmatic plumbing system and eruption-related seismicity. Earthquake locations reveal seismicity below and offset to the north of the volcano along a tectonic fault. Based on historical seismicity, the magma chamber is postulated to have a top at about 6 km depth. Minor dome explosions, large sub-plinian eruptions and dome collapses are analyzed using an automated event classification scheme. Low-frequency tremor, interpreted as gas escape, and low-frequency earthquakes are a dominant proportion of the energy released. We also examine multiplet earthquakes whose behavior during the study period changed significantly and systematically before the largest eruption, demonstrating the potential of tracking multiplets to assess changing conditions with the conduit.     

Ground-based thermal imaging of lava lakes at Erebus volcano,Antarctica

Mount Erebus, a large intraplate stratovolcano dominating Ross Island, Antarctica, hosts the world's only active phonolite lava lakes. The main manifestation of activity at Erebus volcano in December 2004 was as the presence of two convecting lava lakes within an inner crater. The long-lived Ray Lake, ~ 1400 m² in area, was the site of up to 10 small Strombolian eruptions per day. A new but short-lived, ~ 1000–1200 m² lake formed at Werner vent in December 2004 sourced by lava flowing from a crater formed in 1993 by a phreatic eruption. We measured the radiative heat flux from the two lakes in December 2004 using a compact infrared (IR) imaging camera. Daily thermal IR surveys from the Main Crater rim provide images of the lava lake surface temperatures and identify sites of upwelling and downwelling. The radiative heat outputs calculated for the Ray and Werner Lakes are 30–35 MW and 20 MW, respectively. We estimate that the magma flux needed to sustain the combined heat loss is ~ 250–710 kg s^− 1, that the minimum volume of the magma reservoir is 2 km³, and that the radius of the conduit feeding the Ray lake is ~ 2 m.     

Major and trace element geochemistry of neutral and acidic thermal springs at El Chichón volcano,Mexico: Implications for monitoring of the volcanic activity

Four groups of thermal springs with temperatures from 50 to 80 °C are located on the S–SW–W slopes of El Chichón volcano, a composite dome-tephra edifice, which exploded in 1982 with a 1 km wide, 160 m deep crater left. Very dynamic thermal activity inside the crater (variations in chemistry and migration of pools and fumaroles, drastic changes in the crater lake volume and chemistry) contrasts with the stable behavior of the flank hot springs during the time of observations (1974–2005). All known groups of hot springs are located on the contact of the basement and volcanic edifice, and only on the W–SW–S slopes of the volcano at almost same elevations 600–650 m asl and less than 3 km of direct distance from the crater. Three groups of near-neutral (pH ≈ 6) springs at SW–S slopes have the total thermal water outflow rate higher than 300 l/s and are similar in composition. The fourth and farthest group on the western slope discharges acidic (pH ≈ 2) saline (10 g/kg of Cl) water with a much lower outflow rate (< 10 l/s).     

Geology of a submarine volcanic caldera in the Tonga Arc: Dive results

A submersible dive conducted on Volcano #1 located near 21° 09′S–175° 45′W on the Tonga Arc showed that the volcanic edifice with a caldera floor area of 30 km² located at and 450 m deep (b.s.l.=below sea level) was constructed recently during episodic volcanism. The sequential volcanic events are recorded along a faulted terrain formed in response to the collapse of the caldera wall. The post-caldera events are marked by occasional eruptions that have built scoriaceous cones associated with low-temperature hydrothermal venting and localized small-scale collapse features. The stratigraphy of the caldera wall indicates that the volcano was built by explosive volcanism alternating with quieter eruptive events. The repeated, violent explosive events formed ≤ 20 m thick sequences composed of alternating fine-grained ash beds and sand- to boulder-sized pyroclastic layers. During quieter volcanic events, dykes and massive flows intruded and/or accompanied the eruption of the volcaniclastic deposits throughout the sections of the wall explored. Massive columnar-jointed flows consist of viscous, silica-rich lavas forming tabular and giant radial-jointed (GRJ) flows formed in large (> 8 m in diameter) conduits and extruded onto the sea floor. In addition, massive lava flows forming sill-like complexes were observed underneath and near the giant radial-jointed columnar flows. Also, an intermittent quiet type of eruption produced vesicular lava flows, which are interbedded within the pyroclastic layered deposits. The massive and vesicular lavas consist of andesites and dacites with Ca-depleted (pigeonite) and Ca-enriched (salite) pyroxene, and intermediate (andesine-labradorite) to calcic (bytownite) plagioclase. They are depleted in total alkalis (Na₂O + K₂O < 3%), K₂O (< 1%), Zr/Y (< 1.8), Nb/Zr (< 0.01) and light Rare Earth Elements. We interpret that these andesite–dacite series were erupted after undergoing crystal-liquid fractionation in a magma chamber located underneath the caldera floor.     

Hazard map for volcanic ballistic impacts at Popocatépetl volcano (Mexico)

During volcanic explosions, volcanic ballistic projectiles (VBP) are frequently ejected. These projectiles represent a threat to people, infrastructure, vegetation, and aircraft due to their high temperatures and impact velocities. In order to protect people adequately, it is necessary to delimit the projectiles’ maximum range within well-defined explosion scenarios likely to occur in a particular volcano. In this study, a general methodology to delimit the hazard zones for VBP during volcanic eruptions is applied to Popocatépetl volcano. Three explosion scenarios with different intensities have been defined based on the past activity of the volcano and parameterized by considering the maximum kinetic energy associated with VBP ejected during previous eruptions. A ballistic model is used to reconstruct the “launching” kinetic energy of VBP observed in the field. In the case of Vulcanian eruptions, the most common type of activity at Popocatépetl, the ballistic model was used in concert with an eruptive model to correlate ballistic range with initial pressure and gas content, parameters that can be estimated by monitoring techniques. The results are validated with field data and video observations of different Vulcanian eruptions at Popocatépetl. For each scenario, the ballistic model is used to calculate the maximum range of VBP under optimum “launching” conditions: ballistic diameter, ejection angle, topography, and wind velocity. Our results are presented in the form of a VBP hazard map with topographic profiles that depict the likely maximum ranges of VBP under explosion scenarios defined specifically for Popocatépetl volcano. The hazard zones shown on the map allow the responsible authorities to plan the definition and mitigation of restricted areas during volcanic crises.     

Understanding the collapse–eruption link at Stromboli,Italy: A microanalytical study on the products of the recent Secche di Lazzaro phreatomagmatic activity

The Secche di Lazzaro (SDL) phreatomagmatic activity, with the associated Neostromboli sector collapse, represents the most powerful activity of the last 6 ky at Stromboli. As revealed by its present-day activity, Stromboli is one of the most eruptive volcanoes in Italy, and flank instability, along its NW flank, is a common process. Volcano instability is often dramatised by explosive eruptions, thus it is of crucial importance to understand the linking between volcano collapse and the plumbing system itself. The possible role of pre-eruptive magmatic processes as triggers of explosive eruptions can be mainly preserved by minerals and revealed by petrochemical studies. We studied the juvenile components (scoria and pumice) of the pyroclastic deposits from the SDL phreatomagmatic activity with the aim to understand the eruption–collapse link.     

Seismicity and crustal deformation preceding the January 1996 eruptions at Karymsky Volcanic Center,Kamchatka

Two explosive eruptions occurred on 2 January 1996 at Karymsky Volcanic Center (KVC) in Kamchatka, Russia: the first, dacitic, from the central vent of Karymsky volcano, and the second, several hours later, from Karymskoye lake in the caldera of Akademia Nauk volcano. The main significance of the 1996 volcanic events in KVC was the phreatomagmatic eruption in Karymskoye lake, which was the first eruption in this lake in historical time, and was a basaltic eruption at the acidic volcanic center. The volcanic events were associated with the 1 January Ms 6.7 (Mw 7.1) earthquake that occurred at a distance of about 9–17 km southeast from the volcanoes just before the eruptions. We study the long-term (1972–1995) and short-term (1–2 January 1996) characteristics of crustal deformations and seismicity before the double eruptive event in KVC. The 1972–1995 crustal deformation was homogeneous and characterized by a gradual extension with a steady velocity. The seismic activity in 1972–1995 developed at the depth interval from 0 to 20 km below the Akademia Nauk volcano and spread to the southeast along a regional fault. The seismic activity in January 1996 began with a short sequence of very shallow microearthquakes (M ~0) beneath Karymsky volcano. Then seismic events sharply increased in magnitude (up to mb 4.9) and moved along the regional fault to the southeast, culminating in the Ms 6.7 earthquake. Its aftershocks were located to the southeast and northwest from the main shock, filling the space between the two active volcanoes and the ancient basaltic volcano of Zhupanovsky Vostryaki. The eruption in Karymskoye lake began during the aftershock sequence. We consider that the Ms 6.7 earthquake opened the passageway for basic magma located below Zhupanovsky Vostryaki volcano that fed the eruption in Karymskoye lake.     

The February–July 2007 eruption of the summit crater of Klyuchevskoi volcano, Kamchatka

Observations of the summit eruption of Klyuchevskoi volcano in the period from February 15, 2007 to July 9, 2007 are considered. This typical (for this volcano) summit eruption was explosive-effusive in character. The ejectamenta volume is estimated at 0.025 km³. Calculation of active phases of the volcano was carried out in accordance with V.A. Shirokov’s technique. The identified active phases agree well with the eruptive periods. The 2007 summit eruption corresponds to an active phase (May 2006 to May 2009) favorable for the volcano’s eruption. Geodetic observations carried out since 1979 along a radial profile have revealed uplifts and subsidences of the northeastern slope of the volcano. The maximum displacement of 23 cm was recorded in 2007 on the site closest to the volcano crater at a distance of 11 km from the summit crater center. In the course of two previous summit eruptions (2003–2004 and 2005) insignificant uplifts and subsidences of the slope were also noted, although the general ascent of the slope remained. This indicated possible repeated eruptions in the nearest future. Changes in the seismicity before, during and after the eruption are also discussed.     

Seismicity at White Island volcano, New Zealand: a revised classification and inferences about source mechanism

The classification of earthquakes at White Island volcano, New Zealand, has been revised to address problems in existing classification schemes, to better reflect new data and to try to focus more on source processes. Seismicity generated by the direct involvement of magmatic or hydrothermal fluids are referred to as volcanic, and that generated by fault movement in response to stresses caused by those fluids, regional stresses, thermal effects and so on are referred to as volcano-tectonic. Spasmodic bursts form a separate category, as we have insufficient information to classify them as volcanic or volcano-tectonic. Volcanic seismicity is divided into short-duration, long-period volcanic earthquakes, long-duration volcanic earthquakes, and harmonic- and non-harmonic volcanic tremor, while volcano-tectonic seismicity is divided into shallow and deep volcano-tectonic earthquakes. Harmonic volcanic tremor is related to sub-surface intrusive processes, while non-harmonic volcanic tremor originates close to active craters at shallow depth, and usually occurs during eruptive activity. Short-duration, long-period volcanic earthquakes come from a single source close to the active craters, but originate deeper than non-harmonic volcanic tremor, and are not related to eruptive activity. Long-duration volcanic earthquakes often accompany larger discrete eruptions. The waveform of these events consists of an initial low-frequency part from a deep source, and a later cigar-shaped part of mixed frequencies from a shallow crater source.     

Strombolian explosive styles and source conditions: insights from thermal (FLIR) video

Forward Looking Infrared Radiometer (FLIR) cameras offer a unique view of explosive volcanism by providing an image of calibrated temperatures. In this study, 344 eruptive events at Stromboli volcano, Italy, were imaged in 2001–2004 with a FLIR camera operating at up to 30 Hz. The FLIR was effective at revealing both ash plumes and coarse ballistic scoria, and a wide range of eruption styles was recorded. Eruptions at Stromboli can generally be classified into two groups: Type 1 eruptions, which are dominated by coarse ballistic particles, and Type 2 eruptions, which consist of an optically-thick, ash-rich plume, with (Type 2a) or without (Type 2b) large numbers of ballistic particles. Furthermore, Type 2a plumes exhibited gas thrust velocities (>15 m s⁻¹) while Type 2b plumes were limited to buoyant velocities (<15 m s⁻¹) above the crater rim. A given vent would normally maintain a particular gross eruption style (Type 1 vs. 2) for days to weeks, indicating stability of the uppermost conduit on these timescales. Velocities at the crater rim had a range of 3–101 m s⁻¹, with an overall mean value of 24 m s⁻¹. Mean crater rim velocities by eruption style were: Type 1 = 34 m s⁻¹, Type 2a = 31 m s⁻¹, Type 2b = 7 m s⁻¹. Eruption durations had a range of 6–41 s, with a mean of 15 s, similar among eruption styles. The ash in Type 2 eruptions originates from either backfilled material (crater wall slumping or ejecta rollback) or rheological changes in the uppermost magma column. Type 2a and 2b behaviors are shown to be a function of the overpressure of the bursting slug. In general, our imaging data support a broadening of the current paradigm for strombolian behavior, incorporating an uppermost conduit that can be more variable than is commonly considered.