Modelling the impact of a hypothetical sub-Plinian eruption at La Soufrière of Guadeloupe (Lesser Antilles)

This paper describes the development and application of an impact model for a future hypothetical sub-Plinian eruption of La Soufrière of Guadeloupe. The model was designed to assess the impact from either a single or multiple eruption scenarios, each defined in terms of a map of the intensity of three volcanic hazards; volcanogenic earthquake, tephra fallout and pyroclastic density currents. The impact from the three hazards can be assessed independently or alternatively the joint impact of the three hazards can be assessed. The outputs that are produced from the model are; the number of buildings with collapsed roofs, and the number of fatal and non-fatal casualties.     

Reconstruction and analysis of sub-plinian tephra dispersal during the 1530 A.D. Soufrière (Guadeloupe) eruption: Implications for scenario definition and hazards assessment

The last magmatic eruption of Soufrière of Guadeloupe dated at 1530 A.D. (Soufrière eruption) is characterized by an onset with a partial flank-collapse and emplacement of a debris-avalanche that was followed by a sub-plinian VEI 2–3 explosive short-lived eruption (Phase-1) with a column that reached a height between 9 and 12 km producing about 3.9 × 10⁶ m³ DRE (16.3 × 10⁶ m³ bulk) of juvenile products. The column recurrently collapsed generating scoriaceous pyroclastic flows in radiating valleys up to a distance of 5–6 km with a maximum interpolated bulk deposit volume of 11.7 × 10⁶ m³ (5 × 10⁶ m³ DRE). We have used HAZMAP, a numerical simple first-order model of tephra dispersal [Macedonio, G., Costa, A., Longo, A., 2005. A computer model for volcanic ash fallout and assessment of subsequent hazard. Comput. Geosci. 31, 837–845] to reconstruct to a first approximation the potential dispersal of tephra and associated tephra mass loadings generated by the sub-plinian Phase 1 of the 1530 A.D. eruption. We have tested our model on a deterministic average dry season wind profile that best-fits the available data as well as on a set of randomly selected wind profiles over a 5 year interval that allows the elaboration of probabilistic maps for the exceedance of specific tephra mass load thresholds. Results show that in the hypothesis of a future 1530 A.D. scenario, populated areas to a distance of 3–4 km west–southwest of the vent could be subjected to a static load pressure between 2 and 10 kPa in case of wet tephra, susceptible to cause variable degrees of roof damage. Our results provide volcanological input parameters for scenario and event-tree definition, for assessing volcanic risks and evaluating their impact in case of a future sub-plinian eruption which could affect up to 70 000 people in southern Basse-Terre island and the region. They also provide a framework to aid decision-making concerning land management and development. A sub-plinian eruption is the most likely magmatic scenario in case of a future eruption of this volcano which has shown, since 1992, increasing signs of low-energy seismic, thermal, and acid degassing unrest without significant deformation.     

Discriminating the long distance dispersal of fine ash from sustained columns or near ground ash clouds: The example of the Pomici di Avellino eruption (Somma-Vesuvius,Italy)

Ash samples from tephra layers correlated with the Pomici di Avellino (Avellino Pumice) eruption of Somma-Vesuvius were collected in distal archives and their composition and particle morphology investigated in order to infer their behaviour of transportation and deposition. Differences in composition and particle morphologies were recognised for ash particles belonging to the magmatic Plinian and final phreatomagmatic phases of the eruption. The ash particles were dispersed in opposite directions during the two different phases of the eruption, and these directions are also different from that of coarse-grained fallout deposits. In particular, ash generated during magmatic phase and injected in the atmosphere to form a sustained column shows a prevailing SE dispersion, while ash particles generated during the final phreatomagmatic phase and carried by pyroclastic density currents show a general NW dispersion. These opposite dispersions indicate an ash dispersal influenced by both high and low atmosphere dynamics. In particular, the magmatic ash dispersal was first driven by stratospheric wind towards NE and then the falling particles encountered a variable wind field during their settling, which produced the observed preferential SE dispersal. The wind field encountered by the rising ash clouds that accompanied the pyroclastic density currents of the final phreatomagmatic phase was different with respect to that encountered by the magmatic ash, and produced a NW dispersal. These data demonstrate how ash transportation and deposition are greatly influenced by both high and low atmosphere dynamics. In particular, fine-grained particles transported in ash clouds of small-scale pyroclastic density currents may be dispersed over distances and cover areas comparable with those injected into the stratosphere by Plinian, sustained columns. This is a point not completely addressed by present day mitigation plans in case of renewal of activity at Somma-Vesuvius, and can yield important information also for other volcanoes potentially characterised by explosive activity.     

Simulation of the 1980 eruption of Mount St. Helens using the ash-tracking model PUFF

The dispersal of volcanic ash from the May 18, 1980 eruption of Mount St. Helens (MSH) has been simulated using the Lagrangian ash-tracking model PUFF. Previous applications of the model were limited to smaller, short-lived eruptions with ash dispersal occurring mainly within the troposphere. Two high-resolution atmospheric reanalysis datasets (ERA-40 and NCEP/NCAR-40) allowed MSH ash cloud dispersal to be simulated up to 30 km elevation. The 1980 eruption was divided into two distinct eruptive phases, (1) an initial, relatively short-lived blast/surge phase that injected ash up to 30 km and (2) a subsequent nine-hour plinian phase that maintained an average eruption column height of 16 km. Using PUFF, the two phases of the MSH eruption were modeled separately based on a range of individual input parameters and then combined to produce an integrated simulation of the entire eruption. The trajectory and areal extent of the modeled atmospheric ash cloud best match the actual distribution of MSH ash when input parameters are set to values inferred from satellite and radar data collected on May 18, 1980. The prevailing wind field exerts the strongest control on the advection and ultimate position of the modeled ash cloud, making the maximum column height and the vertical distribution of ash the most sensitive of the PUFF input parameters for this event. The results indicate that the PUFF model works well at simulating the dispersal of ash injected well into the lower stratosphere from a moderate, relatively long-lived eruption, such as MSH. However, attempts to use PUFF to recreate some granulometric aspects of the MSH fallout deposit, such as the maximum particle size as a function of distance from source, were not successful. PUFF consistently predicts much greater fallout distances for small ash particles (< 500 µm) than actually observed in the MSH deposit. The effective settling velocities used by the PUFF model appear to be too slow to accurately predict fallout distances of small ash particles. As a consequence the PUFF model may overestimate the duration of ash loading in the atmosphere associated with the distal fine ash component of explosive eruptions.     

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.     

Stretching characteristics and its dynamic significance of the northern continental margin of South China Sea

To investigate the thinning of the whole crust, and the contribution of the upper versus lower crust to the stretching since Cenozoic, we calculated the stretching factor of the northern margin of South China Sea with data such as whole crustal thickness, depth of Cenozoic sedimentary basement and the horizontal displacement of faults. An isometric line drawing on whole crustal stretching factor is then obtained. Along the seismic Line 1530 in Baiyun sag, we also calculated the stretching factors of the upper and lower crust. The results suggest that the whole crustal thickness decreases seaward while the whole crustal stretching factor increases from shelf to slope. The Moho upwells highest beneath where the crust is thinned most. The value of the whole crustal stretching factor ranges from 1.5 to 6. Two areas were thinned intensely: the center of Yinggehai Basin, and the Baiyun sag in the Pearl River Mouth Basin. The calculation of the upper and lower crustal stretching factors from DSP1530 in Baiyun sag shows that the original crust of Baiyun sag should be thinned before deformation. Its pre-Cenozoic evolution as well as tectonic position during Cenozoic might be responsible for that.