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.     

Volcanic tsunamis and prehistoric cultural transitions in Cook Inlet,Alaska

The 1883 eruption of Augustine Volcano produced a tsunami when a debris avalanche traveled into the waters of Cook Inlet. Older debris avalanches and coeval paleotsunami deposits from sites around Cook Inlet record several older volcanic tsunamis. A debris avalanche into the sea on the west side of Augustine Island ca. 450 years ago produced a wave that affected areas 17 m above high tide on Augustine Island. A large volcanic tsunami was generated by a debris avalanche on the east side of Augustine Island ca. 1600 yr BP, and affected areas more than 7 m above high tide at distances of 80 km from the volcano on the Kenai Peninsula. A tsunami deposit dated to ca. 3600 yr BP is tentatively correlated with a southward directed collapse of the summit of Redoubt Volcano, although little is known about the magnitude of the tsunami. The 1600 yr BP tsunami from Augustine Volcano occurred about the same time as the collapse of the well-developed Kachemak culture in the southern Cook Inlet area, suggesting a link between volcanic tsunamis and prehistoric cultural changes in this region of Alaska.     

The volcanic acidification of glacial Lake Caviahue,Province of Neuquen,Argentina

Lake Caviahue (northern Patagonia, Argentina) is a large glacial lake acidified by volcanic fluids from Copahue volcano. The lake and the feeding rivers were sampled annually from 1997 till early 2006, including the eruptive period of 2000. Lake Caviahue waters evolved over time, with the most concentrated waters in 2000 during the eruptive period, followed by gradual dilution that was interrupted by renewed acidification in 2003–2004. Inversion of the lake water data and application of a dynamic non-steady state model for the lake provides our best quantitative estimates for the variation in element fluxes over the 9-year period. The model flux results agree well with most of the measured fluxes. The Copahue hydrothermal system had gently declining element fluxes between 1997 and mid-1999, although the lake was still becoming more concentrated. About 2–3 months before the 2000 eruption, element fluxes increased strongly, but the hydrothermal fluxes almost shutoff directly after the main eruptive events. The fluxes of several elements recovered post-2001, with an increase in element fluxes in 2003–2004; the lake became more dilute between 2004 and 2006. The intrusion of new magma into the hydrothermal system just prior to the 2000 eruption led to enhanced water rock interaction, with higher concentrations of the rock forming elements in the fluids, and the hot spring flow rate increased as a result of the higher pressure in the reservoir. The fluids became saturated in alunite and jarosite, and they were already saturated with anhydrite. Precipitation of these minerals possibly led to a decreased permeability of the hydrothermal reservoir, leading to the strongly reduced element fluxes just after the eruption. In addition, K, Al and S were retained in the newly precipitated minerals as well, further diminishing their export. The acidification in 2003–2004 may have resulted from a new small intrusion of magma or resulted from seismic activity that created new permeability and fresh rock surfaces for water rock interaction. The volcano is a significant source of toxic trace elements such as F, As, B and Li as well as a nutrient (P) for the local watershed. Monitoring of the hydrothermal fluids in the river that drains Copahue, especially the S/Cl, Mg/Cl and Mg/K values as well as the magnitude of the element fluxes would provide the best information for eruption forecasting for this volcano.     

Vent conditions for expected eruptions at Vesuvius

Determining consistent sets of vent conditions for next expected eruptions at Vesuvius is crucial for the simulation of the sub-aerial processes originating the volcanic hazard and the eruption impact. Here we refer to the expected eruptive scales and conditions defined in the frame of the EC Exploris project, and simulate the dynamics of magma ascent along the volcanic conduit for sub-steady phases of next eruptions characterized by intensities of the Violent Strombolian (VS), Sub-Plinian 2 (SP2), and Sub-Plinian 1 (SP1) scale. Sets of conditions for the simulations are determined on the basis of the bulk of knowledge on the past history of Vesuvius [Cioni, R., Bertagnini, A., Santacroce, R., Andronico, D., Explosive activity and eruption scenarios at Somma–Vesuvius (Italy): towards a new classification scheme. Journal of Volcanology and Geothermal Research, this issue.]. Volatile contents (H₂O and CO₂) are parameterized in order to account for the uncertainty in their expected amounts for a next eruption. In all cases the flow in the conduit is found to be choked, with velocities at the conduit exit or vent corresponding to the sonic velocity in the two-phase non-equilibrium magmatic mixture. Conduit diameters and vent mixture densities are found to display minimum overlapping between the different eruptive scales, while exit gas and particle velocities, as well as vent pressures, largely overlap. Vent diameters vary from as low as about 5 m for VS eruptions, to 35–55 m for the most violent SP1 eruption scale. Vent pressures can be as low as less than 1 MPa for the lowest volatile content employed of 2 wt.% H₂O and no CO₂, to 7–8 MPa for highest volatile contents of 5 wt.% H₂O and 2 wt.% CO₂ and large eruptive scales. Gas and particle velocities at the vent range from 100–250 m/s, with a tendency to decrease, and to increase the mechanical decoupling between the phases, with increasing eruptive scale. Except for velocities, all relevant vent quantities are more sensitive to the volatile content of the discharged magma for the highest eruptive scales considered.     

Volcanoes in the pre-Columbian life,legend, and archaeology of Costa Rica (Central America)

Costa Rica is located geographically in the southern part of the Central American Volcanic Front, a zone where interaction between the Mesoamerican and South American cultures occurred in pre-Columbian times. Several volcanoes violently erupted during the Holocene, when the first nomadic human hunters and later settlers were present. Volcanic rocks were the most important geo-resource in making artifacts and as construction materials for pre-Columbian inhabitants. Some pottery products are believed to resemble smoking volcanoes, and the settlements around volcanoes would seem to indicate their influence on daily life. Undoubtedly, volcanic eruptions disrupted the life of early settlers, particularly in the vicinity of Arenal and Irazú volcanoes, where archaeological remains show transient effects and displacement caused by periodical eruptions, but later resilient occupations around the volcanoes. Most native languages are extinct, with the exception of those presently spoken in areas far away from active volcanoes, where no words are related to volcanic phenomena or structures. The preserved legends are ambiguous, suggesting that they were either produced during the early Spanish conquest or were altered following the pre-Columbian period.     

喀左盆地下白垩统九佛堂组油页岩特征及控矿条件

野外地质调查和岩石学、地球化学研究表明：喀左盆地3个油页岩层群分别发育于下白垩统九佛堂组的一、二、三段,以位于三段的第三油页岩层群的工业意义最大;油页岩常与上覆沉凝灰岩或其它凝灰质岩石呈伴生关系。油页岩（矿石）类型有泥岩、具粉砂质条带的含钙-钙质泥岩、泥灰岩和页岩4种类型,以前两者为主;低品级油页岩体积分数为50％,中品级油页岩体积分数为42％,高品级油页岩体积分数为8％,均属高灰分油页岩。油页岩的有机质类型为腐植-腐泥型,源于低等植物;浅湖-半深湖形成的油页岩（矿石）类型为含粉砂质条带的含钙-钙质泥岩,深湖形成的油页岩（矿石）类型主要为泥岩;物源区母岩类型对油页岩的化学成分特征及灰分成分类型具有明显的控制关系;火山喷发的灾变事件终止了油页岩的成矿作用。     

Dike propagation within active central volcanic edifices: constraints from Somma-Vesuvius,Etna and analogue models

Dikes within stratovolcanoes are commonly expected to have radial patterns. However, other patterns may also be found, due to regional stresses, magmatic reservoirs and topographic variations. Here, we investigate dike patterns within volcanic edifices by studying dike and fissure complexes at Somma-Vesuvius and Etna (Italy) using analogue models. At the surface, the dikes and fissures show a radial configuration. At depths of tens to several hundreds of metres, in areas exposed by erosion, tangential and oblique dikes are also present. Analogue models indicate that dikes approaching the flanks of cones, regardless of their initial orientation, reorient to become radial (parallel to the maximum gravitational stress). This re-orientation is a significant process in shallow magma migration and may also control the emplacement of dike-fed fissures reaching the lower slopes of the volcano.     

Phreatomagmatic volcanic hazards where rift-systems meet the sea, a study from Ambae Island, Vanuatu

Ambae Island is a mafic stratovolcano located in the northern Vanuatu volcanic arc and has a NE–SW rift-controlled elongated shape. Several hundred scoria cones and fissure-fed lava fields occur along its long axis. After many decades of quiescence, Ambae Island erupted on the 28th of November 2005, disrupting the lives of its 10,000 inhabitants. Its activity remained focused at the central (crater-lake filled) vent and this is where hazard-assessments were focused. These assessments initially neglected that maars, tephra cones and rings occur at each tip of the island where the eruptive activity occurred < 500 and < 300 yr B.P. The products of this explosive phreatomagmatic activity are located where the rift axis meets the sea. At the NE edge of the island five tephra rings occur, each comparable in size to those on the summit of Ambae. Along the NE coastline, a near-continuous cliff section exposes an up to 25 m thick succession of near-vent phreatomagmatic tephra units derived from closely spaced vents. This can be subdivided into two major lithofacies associations. The first association represents when the locus of explosions was below sea level and comprises matrix-supported, massive to weakly stratified beds of coarse ash and lapilli. These are dominant in the lowermost part of the sequence and commonly contain coral fragments, indicating that the loci of explosion were located within a reef or coral sediment near the syn-eruptive shoreline. The second type indicate more stable vent conditions and rapidly repeating explosions of high intensity, producing fine-grained tephra with undulatory bedding and cross-lamination as well as megaripple bedforms. These surge and fall beds are more common in the uppermost part of the succession and form a few-m-thick pile. An older tephra succession of similar character occurs below, and buried trees in growth position, as well as those flattened within base surge beds. This implies that the centre of this eruption was very near the coastline. The processes implied by these deposits are amongst the most violent forms of volcanism on this island. In addition, the lowland and coastal areas affected by these events are the most heavily populated. This circumstance is mirrored on many similar volcanic islands, including the nearby SW Pacific examples of Taveuni (Fiji), Upolu and Savai'i (Samoa), and Ambrym (Vanuatu). These locations are paradoxically often considered safe areas during summit/central-vent eruptions, simply because they are farthest from the central sources of ash-fall and lahar hazard. The observations presented here necessitate a revision of this view.     

Sedimentological constraints on hydrometeor-enhanced particle deposition: 1992 Eruptions of Crater Peak, Alaska

Water is a dominant component of volcanic clouds and has fundamental control on very fine particle deposition. Particle size characteristics of distal tephra-fall (100s km from source volcano) have a higher proportion of very fine particles compared to predictions based on single particle settling rates. In this study, sedimentological analyses of fallout from for the 18 August and 16–17 September 1992 eruptions of Crater Peak, Alaska, are combined with satellite observations, and cloud trajectory and microphysics modeling to investigate meteorological influences on particle sedimentation. Total grain size distributions of tephra fallout were reconstructed for both Crater Peak eruptions and indicate a predominance of fine particles < 125 μm. Polymodal analysis of the deposits has identified a particle subpopulation with mode ~ 15–18 μm involved in particle aggregation. Accounting for the magmatic water source only, calculated ice water content of the 3.7 hour old September 1992 Spurr cloud was ~ 4.5 × 10^− 2 g m^− 3 (based on an estimated cloud thickness of ~ 1000 m from trajectory modeling). Hydrometeor formation on particles in the volcanic cloud and subsequent sublimation may induce a cloud base instability that leads to rapid bulk (en masse) sedimentation of very fine particles through a mammatus-like mechanism.     

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.