Helium-isotope systematics at Nevado del Ruiz volcano, Colombia: implications for the volcanic hydrothermal system

We have collected 14 water and gas samples from 9 thermal springs and gas vents near Nevado del Ruiz volcano, Colombia. The ³He/⁴He and ⁴He/²⁰Ne ratios vary significantly from 0.98 R_atm (where R_atm is the atmospheric ³He/⁴He ratio of 1.4 × 10⁻⁶) to 6.30 R_atm, and from 0.37 to 7.0, respectively. The ³He/⁴He ratio (corrected for air contamination) decreases with increasing distance from the central crater of the volcano to the sampling site. The trend is very similar to that observed at Ontake volcano, Japan. A hydrodynamic porous-media dispersion model can explain the ³He/⁴He trend. The temporal variations in the ³He/⁴He ratio at four sites provide useful information on the apparent velocity of the magmatic fluid flow brought on by a volcanic eruption. The estimated value of several tens m day⁻¹ agrees well with the inferred velocity of flow in Oshima volcano, Japan and is comparable to the largest rate of groundwater movement in a deep sedimentary basin.     

Pyroclastic deposits of the November 13, 1985 eruption of Nevado del Ruiz volcano, Colombia

The November 13, 1985 eruption of Nevado del Ruiz produced a series of pyroclastic flows and surges that eroded channels on the surface of the summit glacier and generated lahars which descended down most of the rivers that drain the volcano. The stratigraphy of the proximal pyroclastic deposits indicates that there were at least four episodes to the eruption. Episode I, deposited an unusual surge consisting of small pieces of ice mixed with ash and exhibiting planar stratification. Ballistically emplaced fragments are also intercalated with this unit. During Episode II, at least two pyroclastic flows were erupted. Their deposits contain the most evolved pumice of the entire eruption; SiO₂ content of matrix glass ranges between 74.5 and 74.9%. Episode III is marked by the emplacement of a welded tuff with an average SiO₂ content of about 66% in the matrix glass. The final Episode IV was characterized by the development of a high-altitude eruption column and the emplacement of several nonwelded pyroclastic flows. Banded pumice are common in the pyroclastic flow as well as in the pumice fall deposits. Co-existing dark and light pumice bands differ in SiO₂ content by 3.5% and in general are similar to the composition of the welded pumice from Episode III.The compositional zonation of the pyroclastic deposits from Episode I to IV suggests that a nearsurface compositionally-stratified portion of the magma body was tapped during Episode II. During Episodes III and IV the main body of magma was involved although the coexistence of the compositionally distinct pumice clasts at similar stratigraphic levels argues for mixing of magma from different levels in the chamber during the eruptive process.     

Eruption and mass-wasting-induced processes during the late holocene destructive phase of Nevado del Ruiz volcano, Colombia

The November 13, 1985, eruption was characteristic of the Arenas eruptive stage of Nevado del Ruiz, the most recent of a series of twelve eruptive stages that have occurred in the past 11,000 years. Eruptive sequences, deposits and processes similar to that of 1985 have characterized the behavior of Nevado del Ruiz during three major prehistorical and historical eruptive stages: the approximately 3300-3100 yr. B.P. Hedionda, the 16th century Azufrado, and the mid-1800's Lagunillas eruptive stages, that partly destroyed the present Ruiz summit.According to the interpretation of the stratigraphic record of prehistorical eruptions and historical accounts, almost every recent magmatic event was small or short-lived. Nevertheless, rockslide-debris avalanches and catastrophic debris flows were triggered in all the eruptions owing to slope failures related to specific tectonic features of Ruiz volcano and/or to significant interactions between pyroclastic debris and the ice cap. Evidence for headward retreat of avalanche scarps during multiple eruptions reinforce the case that large slope failures can occur repeatedly at a large-volume volcano like Ruiz without reconstruction of the edifice. The latest Ruiz eruptions that involved rockslide-debris avalanches resemble in part the Shiveluch 1964 event for which evidence of lateral blast deposits are lacking, but differ in part from this type because non-eruptive and mass-wasting processes also triggered rockslide-debris avalanches.Many factors render the cluster of domes of the Ruiz summit unstable, including: (1) deeply dissected troughs opened toward the north-northeast (Azufrado), east (Lagunillas), and south (Recio) flanks; (2) strongly hydrothermally altered north and east flanks of the summit; (3) currently glaciated or recently deglaciated, high cliffs; (4) steep unstable margins of the ice cap on the north and east.Thus, in light of its past behavior, a small eruption or an earthquake might trigger catastrophic rockslide-debris avalanches. Furthermore, such avalanches as well as glacial outburst-floods and ice avalanches could induce debris flows by mobilizing weathered, water-saturated, and unconsolidated rocks or deposits.     

The 1985 Nevado del Ruiz volcano catastrophe: anatomy and retrospection

This paper seeks to analyze in an objective way the circumstances and events that contributed to the 1985 Nevado del Ruiz catastrophe, in order to provide useful guidelines for future emergencies. The paper is organized into two principal parts. In the first part, an Anatomy of the catastrophe is developed as a step-by-step chronicle of events and actions taken by individuals and organizations during the period November 1984 through November 1985. This chronicle provides the essential background for the crucial events of November 13. This year-long period is broken down further to emphasize important chapters: the gradual awareness of the awakening of the volcano; a long period of institutional skepticism reflecting an absence of credibility; the closure of the credibility gap with the September 11 phreatic eruption, followed by an intensive effort to gird for the worst; and a detailed account of the day of reckoning. The second part of the paper, Retrospection, examines the numerous complicated factors that influenced the catastrophic outcome, and attempts to cull a few “lessons from Armero” in order to avoid similar occurrences in the future.In a nutshell, the government on the whole acted responsibly but was not willing to bear the economic or political costs of early evacuation or a false alarm. Science accurately foresaw the hazards but was insufficiently precise to render reliable warning of the crucial event at the last possible minute. Catastrophe was therefore a calculated risk, and this combination - the limitations of prediction/detection, the refusal to bear a false alarm and the lack of will to act on the uncertain information available - provided its immediate and most obvious causes. But because the crucial event occurred just two days before the Armero emergency management plan was to be critically examined and improved, the numerous circumstances which delayed progress of emergency management over the previous year also may be said to have contributed to the outcome. Thus the catastrophe was not caused by technological ineffectiveness or defectiveness, nor by an overwhelming eruption, or by an improbable run of bad luck, but rather by cumulative human error - by misjudgment, indecision and bureaucratic shortsightedness. Armero could have produced no victims, and therein dwells its immense tragedy.     

Impacts of acid emissions from Nevado del Ruiz volcano, Colombia, on selected terrestrial and aquatic ecosystems

Emissions of acidic gases and thermal waters from Nevado del Ruiz volcano have recently increased in concert with the November 13, 1985 eruption. This study examines the downwind and downstream effects of these emissions on alpine ecosystems high on the slopes of the volcano (4100 m) and on coffee plantations at lower elevations (< 2000 m) and greater distances from the active vent (> 30 km). Samples of bulk deposition, rain, soils, soil solutions, and streams were collected over a six-month period (January–July, 1987) to examine the impacts of this volcanogenic acidity.Bulk deposition falling on the higher slopes of the volcano is usually acidified; however, deposition reaching the distal coffee plantations seldom is acidic. The sources of the acids are hydrogen chloride and sulfur dioxide in the plume of the volcano. Although sulfur dioxide is by far the more abundant gas, hydrogen chloride is most responsible for acidification of rain falling on the slopes of the volcano. With distance from the vent, the chloride/sulfate ratio drops exponentially. The only major influence on regional precipitation chemistry in addition to the volcano appears to be land-use-related activities around the coffee plantations. Deposition on these areas is enriched by an order of magnitude in nitrate and base cations, compared to all other stations.Throughfall chemistry in the coffee plantations shows a dramatic response to occasional acid-rain events. A base-leaching process on coffee plant leaves is triggered by acid rain. For each equivalent of hydrogen ion in rain on the leaf surface, over 23 equivalents of potassium ion are leached from the leaf.In spite of this dramatic response by the vegetation, the plantation soils appear relatively unaffected by acidic deposition. In contrast, the alpine soils on the volcano exhibit low pHs, high sulfate and chloride concentrations in soil solutions, and high extractable sulfate concentrations. All of these factors indicate that these soils have undergone significant acid loading.While the deposition of the region is acidified by hydrogen chloride, the streams flowing off the volcano are apparently acidified by sulfuric acid in thermal waters discharging into the streams. The acidity of these streams decreases downstream, while the silica concentrations increase downstream. The composition of stream water is most influenced by thermal-water discharges as well as equilibrium dissolution of amorphous silica glass and non-equilibrium leaching of unweathered ash.The impacts of acid gases and thermal water released from the volcano appear to be restricted to ecosystems on the slope of the volcano. The only impact of Nevado del Ruiz on surrounding coffee plantations appears to be potassium leaching of coffee leaves from occasional acid-deposition events.     

Analysis of seismic patterns observed at Nevado del Ruiz volcano, Colombia during August – September 1985

The seismic activity of the Nevado del Ruiz volcano was monitored during August–September 1985 using a three-component portable seismograph station placed on the upper part of the volcano. The objective was to investigate the frequency content of the seismic signals and the possible sources of the volcanic tremor. The seismicity showed a wide spectrum of signals, especially at the beginning of September. Some relevant patterns from the collected records, which have been analyzed by spectrum analysis, are presented. For the purpose of analysis, the records have been divided into several categories such as long-period events, tremor, cyclic tremor episodes, and strong seismic activity on September 8, 1985.The origin of the seismic signals must be considered in relation to the dynamical and acoustical properties of fluids and the shape and dimensions of the volcano's conduits.The main results of the present experiment and analysis show that the sources of the seismic signals are within the volcanic edifice. The signal characteristics indicate that the sources lie in fluid-phase interactions rather than in brittle fracturing of solid components.     

Analysis of swarms of high-frequency seismic events at Nevado del Ruiz volcano, Colombia (January 1986–August 1987): Development of a procedure

The determination of signal properties like frequency, amplitude, “signature” of the sequence of arrivals of P- and S-phases and the subsequent computation of hypocenter parameters, are the steps of a procedure in the analysis of swarms of high-frequency seismic events from Nevado del Ruiz volcano, which could be useful in other volcanic areas. The use of this procedure on Ruiz data led to the determination of various seismic source volumes, which together with hypocenters of non-swarmed events and the identification of geological features such as faults which cross the area, alignment of domes and volcanic edifices, and the record of low-frequency seismic events not centered in the area of the active vent, suggest the existence of a magma chamber of a rather complex geometry. This hypothesis would suggest that magma, emplaced in zones of weakness (faults) that affect the area, may act as the source of seismic activity. This demands a revision of the methodology used in the deformation measurements, and in the analysis of released energy and event occurrence too. The relations between hypocenters (sources), location and relative energy of the swarms, as well as the general volcanic activity (tremor, ash and gas emission, deformation) could turn out to be a possible means to assist in the identification of premonitory activity and thus be relevant to the process of monitoring and forecasting.     

Evaluation of magmatic processes for the products of the Nevado del Ruiz Volcano, Colombia from geochemical and petrological data

Trace-element concentrations of the eruptive products from the Nevado del Ruiz volcano indicate that at least for the last 1 Ma, the major controlling factor in the evolution of the magma has been a simple fractional crystallization.The Zr/Hf ratios suggest that the source material is mantle-like in origin and whilst subduction-related contamination of the source material has taken place, there is no evidence that assimilation of continental crust has affected the series.No convincing geochemical evidence for magma mixing was found, and data on banded pumices of Pleistocene, Historic and the 1985 eruptions strongly support a crystal fractionation model to explain the compositional variation of their glasses. Therefore, if mixing has taken place it can only be a mixing of layers within a stratified or partially stratified magma chamber.Modeling of fractional crystallization using a stepwise program has been undertaken and least-squares approximations with small residuals are commensurate with analyzed samples. Crystal fractionation models with plagioclase dominating two-pyroxenes and iron oxides give good fits for all calculations. We conclude that simple crystal fractionation is the process most consistent with major- and trace-element variations in the Ruiz series. This closed-system model requires at least 77% crystal fractionation of a basaltic parent to generate the observed compositions of dacites.     

Hydrological study of the potential effects of the melting of Nevado del Ruiz glacier on urban growth zones in Manizales,Colombia

ABSTRACT

This paper presents the results of a hydrological study aimed at characterizing flood-prone areas in the urban growth zone in the city of Manizales based on the potential effects of melting of the Nevado del Ruiz glacier, in Caldas, Colombia. These results constitute a basis for decision making regarding possible urban growth zones in Andean areas that face risks from volcanic eruptions producing lahars and floods caused by glacier melt. Conservative estimates of extreme flows in the Chinchiná River in the urban growth area of El Rosario can be obtained by considering the effects of rain triggered by airborne particulate material following a volcanic explosion combined with the effects of glacier melt. The effects of global warming on tropical glaciers contribute to their retreat, leading to their disappearance. Therefore, the worst scenario would take place if these events occurred in the short term as glacier volume decreases with time.

Editor Z.W. Kundzewicz; Associate editor not assigned     

Water contents, temperatures and diversity of the magmas of the catastrophic eruption of Nevado del Ruiz, Colombia, November 13, 1985

The petrology of the highly phyric two-pyroxene andesitic to dacitic pyroclastic rocks of the November 13, 1985 eruption of Nevado del Ruiz, Colombia, reveals evidence of: (1) increasingly fractionated bulk compositions with time; (2) tapping of a small magma chamber marginally zoned in regard to H₂O contents (1 to 4%), temperature (960–1090°C), and amount of residual melt (35 to 65%); (3) partial melting and assimilation of degassed zones in the hotter less dense interior of the magma chamber; (4) probable heating, thermal disruption and mineralogic and compositional contamination of the magma body by basaltic magma “underplating”; and (5) crustal contamination of the magmas during ascent and within the magma chamber. Near-crater fall-back or “spill-over” emitted in the middle of the eruptive sequence produced a small pyroclastic flow that became welded in its central and basal portions because of ponding and thus heat conservation on the flat glaciated summit near the Arenas crater. The heterogeneity of Ruiz magmas may be related to the comparatively small volume (0.03 km³) of the eruption, nearly ten times less than the 0.2 km³ of the Plinian phase of Mount St. Helens, and probable steep thermal and P_H2O gradients of a small source magma chamber, estimated at 300 m long and 100 m wide for an assumed ellipsoidal shape.     

Chronology of the principal scientific and governmental actions leading up to the November 13, 1985 eruption of Nevado del Ruiz, Colombia

The virtue of studying the chronology of an event is the insight it affords as to how man views and responds to the forces of nature. This chronology reflects the fumbling nature of mankind (as only appreciated in hindsight) when faced by an impending event of uncertain nature, in this case in a developing nation, but it might have happened anywhere. It is hoped that we as volcanologists might improve our response to such crises knowing how they have evolved previously.     

Remanent magnetization of andesitic and dacitic pumice from the 1985 eruption of Nevado del Ruiz (Colombia) reversed due to self-reversal

The natural remanent magnetization of andesitic pumice emitted during the 1985 eruption of the Nevado del Ruiz volcano (Colombia) has a direction opposite to the present geomagnetic field. The self-reversing mechanism can be re-activated in the laboratory during cycles of heating and subsequent cooling in air and zero magnetic field. Laboratory-produced thermoremanent magnetization is dominated by the same self-reversal process in fields up to several mT. Microchemical, optical and Curie temperature analyses indicate that the ferromagnetic minerals are members of the magnetite-ulvöspinel and hematite-ilmenite series with average compositions of Fe_2.73Ti_0.27O₄ and Fe_1.38Ti_0.62O₃, respectively. In analogy with the magnetic behaviour of synthetically grown antiferromagnetic-ferromagnetic FeMn-FeNi films, the self-reversal can probably be interpreted in terms of an exchange field acting between a Ti-poor canted antiferromagnetic and a Ti-rich ferrimagnetic phase in the hemoilmenite grains.     

Anhydrite in Nevado del Ruiz November 1985 pumice: relevance to the sulfur problem

The purpose of this letter is to document the occurrence of anhydrite within a pumice phenocryst, and to determine the sulfur contents of the quenched melts in this pumice and to consider if the anhydrite could have crystallized from the magma. Additionally, the presence is noted in the Central Cordillera, potentially below Nevado del Ruiz, of sulfide- and sulfate-bearing deposits of possible relevance to this problem.     

TOMS measurement of the sulfur dioxide emitted during the 1985 Nevado del Ruiz eruptions

The eruptions of Nevado del Ruiz in 1985 were unusually rich in sulfur dioxide. These eruptions were observed with the Nimbus 7 Total Ozone Mapping Spectrometer (TOMS) which can quantitatively map volcanic sulfur dioxide plumes on a global scale. A small eruption, originally believed to be of phreatic origin, took place on September 11, 1985. However, substantial amounts of sulfur dioxide from this eruption were detected with TOMS on the following day. The total mass of SO₂, approximately 9 ± 3 × 10⁴ metric tons, was deposited in two clouds, one in the upper troposphere, the other possibly at 15 km near the stratosphere.The devastating November 13 eruptions were first observed with TOMS at 1150 EST on November 14. Large amounts of sulfur dioxide were found in an arc extending 1100 km from south of Ruiz northeastward to the Gulf of Venezuela and as an isolated cloud centered at 7°N on the Colombia-Venezuela border. On November 15 the plume extended over 2700 km from the Pacific Ocean off the Colombia coast to Barbados, while the isolated mass was located over the Brazil-Guyana border, approximately 1600 km due east of the volcano. Based on wind data from Panama, most of the sulfur dioxide was located at 10–16 km in the troposphere and a small amount was quite likely deposited in the stratosphere at an altitude above 24 km.The total mass of sulfur dioxide in the eruption clouds was approximately 6.6 ± 1.9 × 10⁵ metric tons on November 14. When combined with quiescent sulfur dioxide emissions during this period, the ratio of sulfur dioxide to erupted magma from Ruiz was an order of magnitude greater than in the 1982 eruption of El Chichon or the 1980 eruption of Mount St. Helens.     

Challenges of modeling current very large lahars at Nevado del Huila Volcano,Colombia

Nevado del Huila, a glacier-covered volcano in the South of Colombia’s Cordillera Central, had not experienced any historical eruptions before 2007. In 2007 and 2008, the volcano erupted with phreatic and phreatomagmatic events which produced lahars with flow volumes of up to about 300 million m³ causing severe damage to infrastructure and loss of lives. The magnitude of these lahars and the prevailing potential for similar or even larger events, poses significant hazards to local people and makes appropriate modeling a real challenge. In this study, we analyze the recent lahars to better understand the main processes and then model possible scenarios for future events. We used lahar inundation depths, travel duration, and flow deposits to constrain the dimensions of the 2007 event and applied LAHARZ and FLO-2D for lahar modeling. Measured hydrographs, geophone seismic sensor data and calculated peak discharges served as input data for the reconstruction of flow hydrographs and for calibration of the models. For model validation, results were compared with field data collected along the Páez and Simbola Rivers. Based on the results of the 2007 lahar simulation, we modeled lahar scenarios with volumes between 300 million and 1 billion m³. The approach presented here represents a feasible solution for modeling high-magnitude flows like lahars and allows an assessment of potential future events and related consequences for population centers downstream of Nevado del Huila.     

A warning model based on temporal changes of coda Q for volcanic activity at Nevado Del Ruiz Volcano, Colombia

The coda Q has been calculated for Nevado del Ruiz Volcano, Colombia (NRV) from 1985 to 1999 by using a single scattering model. During this period, the inverse of Q (Q^-1 proportional to attenuation) exhibited a long-term decrease with time, as well as shorter-term variations related to the volcanic activity. Q^-1 increased prior to volcanic crises and decreased afterward. Based on these observations, a seismic warning criterion has been developed. The parameters (frequency band, size of moving average window, and threshold levels) necessary to evidence clear and significant short-term changes in Q^-1 have been investigated and appropriated values are proposed. We suggest a phenomenological model with three stages for the short-term temporal changes in Q^-1 at NRV. Firstly, Q^-1 increases before a volcanic crises because of accumulation of gas and/or liquid, which decreases the aspect ratio of fluid pockets and increases the fractional volume of fluid in the rocks and the pore aspect ratio. Secondly, Q^-1 starts to decrease during the crises by the discharging of fluids such as gas, water, etc. from the volcano. Finally, Q^-1 becomes more stable after the crisis at a lower value because of the degassing and/or increasing of rigidity of the medium because of the long-term crystallization and cooling processes. Q^-1 seems to be a promising monitoring tool at NRV. It is possible that the observed temporal changes of Q^-1, combined with other parameters, may help to predict with greater accuracy a volcanic crisis at NRV.     

Electrical resistivity tomography study of Taal volcano hydrothermal system, Philippines

Taal volcano (311?m in altitude) is located in The Philippines (14°N, 121°E) and since 1572 has erupted 33 times, causing more than 2,000 casualties during the most violent eruptions. In March 2010, the shallow structures in areas where present-day surface activity takes place were investigated by DC resistivity surveys. Electrical resistivity tomography (ERT) lines were performed above the two identified hydrothermal areas located on the northern flank of the volcano and in the Main Crater, respectively. Due to rough topography, deep valleys, and dense vegetation, most measurements were collected using a remote method based on a laboratory-made equipment. This allowed retrieval of information down to a depth of 250?m. ERTs results detail the outlines of the two geothermal fields defined by previous self-potential, CO₂ soil degassing, ground temperature, and magnetic mapping (Harada et al. Japan Acad Sci 81:261–266, 2005; Zlotnicki et al. Bull Volcanol 71:29–49, 2009a, Phys Chem Earth 34:294–408, 2009b). Hydrothermal fluids originate mainly from inside the northern part of the Main Crater at a depth greater than the bottom of the Crater Lake, and flow upward to the ground surface. Furthermore, water from the Main Crater Lake infiltrates inside the surrounding geological formations. The hydrothermal fluids, outlined by gas releases and high temperatures, cross the crater rim and interact with the northern geothermal field located outside the Main Crater.     

Volcanic hazard zonation of the Nevado de Toluca volcano,México

The Nevado de Toluca is a quiescent volcano located 20 km southwest of the City of Toluca and 70 km west of Mexico City. It has been quiescent since its last eruptive activity, dated at ∼ 3.3 ka BP. During the Pleistocene and Holocene, it experienced several eruptive phases, including five dome collapses with the emplacement of block-and-ash flows and four Plinian eruptions, including the 10.5 ka BP Plinian eruption that deposited more than 10 cm of sand-sized pumice in the area occupied today by Mexico City. A detailed geological map coupled with computer simulations (FLOW3D, TITAN2D, LAHARZ and HAZMAP softwares) were used to produce the volcanic hazard assessment. Based on the final hazard zonation the northern and eastern sectors of Nevado de Toluca would be affected by a greater number of phenomena in case of reappraisal activity. Block-and-ash flows will affect deep ravines up to a distance of 15 km and associated ash clouds could blanket the Toluca basin, whereas ash falls from Plinian events will have catastrophic effects for populated areas within a radius of 70 km, including the Mexico City Metropolitan area, inhabited by more than 20 million people. Independently of the activity of the volcano, lahars occur every year, affecting small villages settled down flow from main ravines.