Magmatic gas scrubbing: implications for volcano monitoring

Despite the abundance of SO_2(g) in magmatic gases, precursory increases in magmatic SO_2(g) are not always observed prior to volcanic eruption, probably because many terrestrial volcanoes contain abundant groundwater or surface water that scrubs magmatic gases until a dry pathway to the atmosphere is established. To better understand scrubbing and its implications for volcano monitoring, we model thermochemically the reaction of magmatic gases with water. First, we inject a 915°C magmatic gas from Merapi volcano into 25°C air-saturated water (ASW) over a wide range of gas/water mass ratios from 0.0002 to 100 and at a total pressure of 0.1 MPa. Then we model closed-system cooling of the magmatic gas, magmatic gas-ASW mixing at 5.0 MPa, runs with varied temperature and composition of the ASW, a case with a wide range of magmatic–gas compositions, and a reaction of a magmatic gas–ASW mixture with rock. The modeling predicts gas and water compositions, and, in one case, alteration assemblages for a wide range of scrubbing conditions; these results can be compared directly with samples from degassing volcanoes. The modeling suggests that CO_2(g) is the main species to monitor when scrubbing exists; another candidate is H₂S_(g), but it can be affected by reactions with aqueous ferrous iron. In contrast, scrubbing by water will prevent significant SO_2(g) and most HCl_(g) emissions until dry pathways are established, except for moderate HCl_(g) degassing from pH<0.5 hydrothermal waters. Furthermore, it appears that scrubbing will prevent much, if any, SO_2(g) degassing from long-resident boiling hydrothermal systems. Several processes can also decrease or increase H_2(g) emissions during scrubbing making H_2(g) a poor choice to detect changes in magma degassing.We applied the model results to interpret field observations and emission rate data from four eruptions: (1) Crater Peak on Mount Spurr (1992) where, except for a short post-eruptive period, scrubbing appears to have drastically diminished pre-, inter-, and post-eruptive SO_2(g) emissions, but had much less impact on CO_2(g) emissions. (2) Mount St. Helens where scrubbing of SO_2(g) was important prior to and three weeks after the 18 May 1980 eruption. Scrubbing was also active during a period of unrest in the summer of 1998. (3) Mount Pinatubo where early drying out prevented SO_2(g) scrubbing before the climactic 15 June 1991 eruption. (4) The ongoing eruption at Popocatépetl in an arid region of Mexico where there is little evidence of scrubbing.In most eruptive cycles, the impact of scrubbing will be greater during pre- and post-eruptive periods than during the main eruptive and intense passive degassing stages. Therefore, we recommend monitoring the following gases: CO_2(g) and H₂S_(g) in precursory stages; CO_2(g), H₂S_(g), SO_2(g), HCl_(g), and HF_(g) in eruptive and intense passive degassing stages; and CO_2(g) and H₂S_(g) again in the declining stages. CO_2(g) is clearly the main candidate for early emission rate monitoring, although significant early increases in the intensity and geographic distribution of H₂S_(g) emissions should be taken as an important sign of volcanic unrest and a potential precursor. Owing to the difficulty of extracting SO_2(g) from hydrothermal waters, the emergence of >100 t/d (tons per day) of SO_2(g) in addition to CO_2(g) and H₂S_(g) should be taken as a criterion of magma intrusion. Finally, the modeling suggests that the interpretation of gas-ratio data requires a case-by-case evaluation since ratio changes can often be produced by several mechanisms; nevertheless, several gas ratios may provide useful indices for monitoring the drying out of gas pathways.     

Thermal, seismic and infrasonic evidences of variable degassing rates at Stromboli volcano

During June 1999, we measured the amplitude and rate (number of events per second) at which gas exited the vent at Stromboli volcano as discrete gas bursts or puffs. This allowed us to identify two styles of gas burst (puffing) activity. The first is characterized by frequent, rapidly rising puffs, the second by less frequent, slowly rising puffs. Each style persisted over 5–40-min-long durations and was associated with a high and low number of strombolian explosions per hour, respectively. Each period was also associated with characteristic delay times between the arrival of the infrasonic and thermal signals during strombolian explosions; the delays were longer during vigorous puffing periods. To explain our observations, we propose a model in which the degassing process cycles between vigorous and weak degassing phases. During vigorous degassing phases, bubble layers ascend the conduit at a frequency of 0.5–1.0 s⁻¹. This high degassing level reflects a gas-rich magma column and leads to an increased rate in the formation of shallow foams and, hence, an increase in puffing and explosive activity, as well as a higher free surface level and/or gas jet velocity. During weak phases, bubble layers ascend the conduit at a reduced frequency of 0.2–0.3 s⁻¹. During such times the magma column is poor in gas. This leads to a decreased rate of foam layer formation and hence a reduction in puffing and explosive activity, as well as a lower free surface level and/or gas jet velocity. Variations in puffing activity can thus be used to track changes in the rate at which the shallow system is supplied by fresh, gas-rich magma. Our observations indicate that the two degassing styles last from 5 to 40 min and that the switch from one to the other occurs over a matter of minutes.     

Sulfur dioxide emissions and degassing behavior of Erebus volcano,Antarctica

Emission rates of sulfur dioxide (SO₂) were measured at Erebus volcano, Antarctica in December between 1992 and 2005. Since 1992 SO₂ emissions rates are normally distributed with a mean of 61 ± 27 Mg d^− 1 (0.7 ± 0.3 kg s^− 1) (n = 8064). The emission rates vary over minutes, hours, days and years. Hourly and daily variations often show systematic and cyclic trends. Long-wavelength, large amplitude trends appear related to lava lake area and both are likely controlled by processes occurring at depth. Time series analysis of continuous sequences of measurements obtained over periods of several hours reveals periodicity in SO₂ output ranging from 10 to 360 min, with a 10 min cycle being the most dominant. Closed and open-system degassing models are considered to explain observed variable degassing rates. Closed-system degassing is possible as rheological stiffening and stick/slip may occur within the system. However, the timescales represented in these models do not fit observations made on Erebus. Open-system degassing and convection fits the observations collected as the presented models were developed for a system similar to Erebus in terms of degassing, eruptive activity and process repose time. We show that with the observed emission rate (0.71 kg s^− 1) and a crystal content of 30%, magma will cool 65 °C to match observed heat fluxes; this cooling is sufficient enough to drive convection.     

Monitoring volcanic hazard using eddy covariance at Solfatara volcano, Naples, Italy

An eddy covariance (EC) station was deployed at Solfatara crater, Italy, June 8–25, 2001 to assess if EC could reliably monitor CO₂ fluxes continuously at this site. Deployment at six different locations within the crater allowed areas of focused gas venting to be variably included in the measured flux. Turbulent (EC) fluxes calculated in 30-min averages varied between 950 and 4460 g CO₂ m⁻² d⁻¹; the highest measurements were made downwind of degassing pools. Comparing turbulent fluxes with chamber measurements of surface fluxes using footprint models in diffuse degassing regions yielded an average difference of 0% (±4%), indicating that EC measurements are representative of surface fluxes at this volcanic site. Similar comparisons made downwind of degassing pools yielded emission rates from 12 to 27 t CO₂ d⁻¹ for these features. Reliable EC measurements (i.e. measurements with sufficient and stationary turbulence) were obtained primarily during daytime hours (08:00 and 20:00 local time) when the wind speed exceeded 2 m s⁻¹. Daily average EC fluxes varied by ±50% and variations were likely correlated to changes in atmospheric pressure. Variations in CO₂ emissions due to volcanic processes at depth would have to be on the same order of magnitude as the measured diurnal variability in order to be useful in predicting volcanic hazard. First-order models of magma emplacement suggest that emissions could exceed this rate for reasonable assumptions of magma movement. EC therefore provides a useful method of monitoring volcanic hazard at Solfatara. Further, EC can monitor significantly larger areas than can be monitored by previous methods.     

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.     

Leakage of Active Crater lake brine through the north flank at Rincón de la Vieja volcano, northwest Costa Rica, and implications for crater collapse

The Active Crater at Rincón de la Vieja volcano, Costa Rica, reaches an elevation of 1750 m and contains a warm, hyper-acidic crater lake that probably formed soon after the eruption of the Rio Blanco tephra deposit approximately 3500 years before present. The Active Crater is buttressed by volcanic ridges and older craters on all sides except the north, which dips steeply toward the Caribbean coastal plains. Acidic, above-ambient-temperature streams are found along the Active Crater's north flank at elevations between 800 and 1000 m. A geochemical survey of thermal and non-thermal waters at Rincón de la Vieja was done in 1989 to determine whether hyper-acidic fluids are leaking from the Active Crater through the north flank, affecting the composition of north-flank streams.Results of the water-chemistry survey reveal that three distinct thermal waters are found on the flanks of Rincón de la Vieja volcano: acid chloride–sulfate (ACS), acid sulfate (AS), and neutral chloride (NC) waters. The most extreme ACS water was collected from the crater lake that fills the Active Crater. Chemical analyses of the lake water reveal a hyper-acidic (pH0) chloride–sulfate brine with elevated concentrations of calcium, magnesium, aluminum, iron, manganese, copper, zinc, fluorine, and boron. The composition of the brine reflects the combined effects of magmatic degassing from a shallow magma body beneath the Active Crater, dissolution of andesitic volcanic rock, and evaporative concentration of dissolved constituents at above-ambient temperatures. Similar cation and anion enrichments are found in the above-ambient-temperature streams draining the north flank of the Active Crater. The pH of north-flank thermal waters range from 3.6 to 4.1 and chloride:sulfate ratios (1.2–1.4) that are a factor of two greater than that of the lake brine (0.60). The waters have an ACS composition that is quite different from the AS and NC thermal waters that occur along the southern flank of Rincón de la Vieja.The distribution of thermal water types at Rincón de la Vieja strongly indicates that formation of the north-flank ACS waters is not due to mixing of shallow, steam-heated AS water with deep-seated NC water. More likely, hyper-acidic brines formed in the Active Crater area are migrating through permeable zones in the volcanic strata that make up the Active Crater's north flank. Dissolution and shallow subsurface alteration of north-flank volcanoclastic material by interaction with acidic lake brine, particularly in the more permeable tephra units, could weaken the already oversteepened north flank of the Active Crater. Sector collapse of the Active Crater, with or without a volcanic eruption, represents a potential threat to human lives, property, and ecosystems at Rincón de la Vieja volcano.     

Magnetic monitoring at Merapi volcano, Indonesia

Merapi volcano, located 30 km north of the heavily populated city of Yogjakarta, Java, is one of the most active of the 129 volcanoes in Indonesia. About every 2 years a new phase of activity is observed. Depending on the past activity the unrest gives rise either to an endogenous dome which partly collapses in the southwest direction or to pyroclastic flows which travel as far as 15 km. The 1990–1997 period has involved a plume emission on 30 August 1990, an extrusion on 20 January 1992, and a pyroclastic eruption on 22 November 1994. The intensity of the Earth magnetic field has been measured simultaneously and digitally recorded at four stations since 1990. Two Overhauser magnetometers with resolution of 0.01 nT have been installed in the summit area to strengthen the volcano monitoring. Outstanding magnetic changes appear to correlate with volcanic activity. Three types of volcanomagnetic signals can be identified: long-term trends up to 15 nT with period >10 years; medium-term cyclic variations, at most 3 nT in amplitude and with 1–2 years period; and small events, reaching 1.5 nT, lasting a few months, and associated with any remarkable volcanic activity. Merapi volcano began a new cycle of activity in 1995 leading to a dome growth in July 1996, and accompanied by 27 nuées ardentes in August. The comparison between magnetic data, seismicity, and surface phenomena suggests that some long-term trends of decade periods could be of thermomagnetic origin, while mid-term volcanomagnetic variations associated with the cycles of Merapi activity could be of piezomagnetic origin. Short-term variations of a few weeks duration, less than 1.5 nT, are well correlated with the 1995–1996 seismic activity.     

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.     

Environmental influences on soil CO2 degassing at Furnas and Fogo volcanoes (São Miguel Island,Azores archipelago)

Since October 2001, four soil CO₂ flux stations were installed in the island of São Miguel (Azores archipelago), at Fogo and Furnas quiescent central volcanoes. These stations perform measurements by the accumulation chamber method and, as the gas flux may be influenced by external variables, the stations are equipped with several meteorological sensors. Multivariate regression analysis applied to the large datasets obtained allowed observing that the meteorological variables may influence the soil CO₂ flux oscillations from 18% to 50.5% at the different monitoring sites. Additionally, it was observed that meteorological variables (mainly soil water content, barometric pressure, wind speed and rainfall) play a different role in the control of the gas flux, depending on the selected monitoring site and may cause significant short-term (spike-like) fluctuations. These divergences may be potentially explained by the porosity and hydraulic conductivity of the soils, topographic effects, drainage area and different exposure of the monitoring sites to the weather conditions. Seasonal effects are responsible for long-term oscillations on the gas flux.     

Gravity changes and passive SO2 degassing at the Masaya caldera complex, Nicaragua

An understanding of the mechanisms responsible for persistent volcanism can be acquired through the integration of geophysical and geochemical data sets. By interpreting data on micro-gravity, ground deformation and SO₂ flux collected at Masaya Volcano since 1993, it is now clear that the characteristically cyclical nature of the activity is not driven by intrusion of additional magma into the system. Rather, it may be due in large part to the blocking and accumulation of gas by restrictions in the volcano substructure. The history of crater collapse and formation of caverns beneath the crater floor would greatly facilitate the trapping and storage of gas in a zone immediately beneath San Pedro and the other craters. Another mechanism that may explain the observed gravity and gas flux variations is the convective overturn of shallow, pre-existing, degassed, cooled, dense magma that is replaced periodically by lower density, hot, gas-rich magma from depth. Buoyant gas-rich magma rises from depth and is emplaced near the surface, resulting in the formation and fluctuation of a low-density gas-rich layer centred beneath Nindirí and Santiago craters. As this magma vigorously degasses, it must cool, increase in density and eventually sink. Five stages of activity have been identified at Masaya since 1853 and the most recent data suggest that the system may have been entering another period of reduced degassing in 2000. This type of analysis has important implications for hazard mitigation because periods of intense degassing are associated with poor agricultural yields and reduced quality of life. A better understanding of persistent cyclically active volcanoes will allow for more effective planning of urban development and agricultural land use.     

Coupled degassing and crystallization: experimental study at continuous pressure drop, with application to volcanic bombs

 Experiments on degassing of water-saturated granite melts with a pressure drop from 100 and 450 MPa to 40 and 120 MPa, respectively, at temperatures close to feldspar liquidus (750–700  °C), were carried out to determine the modality of water exsolution and vesicle formation at the liquidus temperature. Pressure-drop rates as small as approximately 100 bar/day were used. Uniform space distributions of bubbles of exsolved water were obtained with starting glass containing a small fraction (≈0.5 vol.%) of trapped air bubbles. Volume crystallization of feldspar was observed in degassed melts supplied with seeds. Bubble size distributions (BSD) measured in granite glasses after degassing are presented. Data on vesicle characteristics (number, radius, area, elongation) were acquired on images digitized with standard software, while the reconstruction of size distributions was performed with the Schwartz-Saltikov "unfolding" procedure. Bubble size distributions of size classes in the range 5–1000 μm were acquired with proper magnification and satisfactory statistical reliability of determined number densities. The BSDs of the experimental samples are compared with the results of measurements of rapidly degassed products of Mt. Etna and Vulcano Island. Many particular features of the bubble nucleation and growth can be distinguished in an individual BSD. However, the general BSD of the whole data set, including natural ones, can be relatively well described with linear regression in bilogarithmic coordinates. The slope of this regression is approximately 2.8±0.1. This dependence is in striking contrast with distributions theoretically predicted with classical nucleation models based on homogeneous nucleation of vesicles. The theoretical distribution requires the occurrence of strong maxima that are not observed in our experimental and natural samples, thus arguing for heterogeneous nucleation mechanisms. Received: 1 October 1998 / Accepted: 25 June 1999     

Diffuse CO<Subscript>2</Subscript> degassing at Vesuvio,Italy

At Vesuvio, a significant fraction of the rising hydrothermal–volcanic fluids is subjected to a condensation and separation process producing a CO₂–rich gas phase, mainly expulsed through soil diffuse degassing from well defined areas called diffuse degassing structures (DDS), and a liquid phase that flows towards the outer part of the volcanic cone. A large amount of thermal energy is associated with the steam condensation process and subsequent cooling of the liquid phase. The total amount of volcanic–hydrothermal CO₂ discharged through diffuse degassing has been computed through a sequential Gaussian simulation (sGs) approach based on several hundred accumulation chamber measurements and, at the time of the survey, amounted to 151 t d^–1. The steam associated with the CO₂ output, computed assuming that the original H₂O/CO₂ ratio of hydrothermal fluids is preserved in fumarolic effluents, is 553 t d^–1, and the energy produced by the steam condensation and cooling of the liquid phase is 1.47×10¹² J d^–1 (17 MW). The location of the CO₂ and temperature anomalies show that most of the gas is discharged from the inner part of the crater and suggests that crater morphology and local stratigraphy exert strong control on CO₂ degassing and subsurface steam condensation. The amounts of gas and energy released by Vesuvio are comparable to those released by other volcanic degassing areas of the world and their estimates, through periodic surveys of soil CO₂ flux, can constitute a useful tool to monitor volcanic activity.Editorial responsibility: H. Shinohara     

Anomalous Diffuse CO2 Emission prior to the January 2002 Short-term Unrest at San Miguel Volcano, El Salvador, Central America

On January 16, 2002, short-term unrest occurred at San Miguel volcano. A gas-and-steamash plume rose a few hundred meters above the summit crater. An anomalous microseismicity pattern, about 75 events between 7:30 and 10:30 hours, was also observed. Continuous monitoring of CO₂ efflux on the volcano started on November 24, 2001, in the attempt to provide a multidisciplinary approach for its volcanic surveillance. The background mean of the diffuse CO₂ emission is about 16 g m^-2 d^-1, but a 17- fold increase, up to 270 g m^-2 d^-1, was detected on January 7, nine days before the January 2002 short-term unrest at San Miguel volcano. These observed anomalous changes on diffuse CO₂ degassing could be related to either a sharp increase of CO₂ pressure within the volcanic-hydrothermal system or degassing from an uprising fresh gas-rich magma within the shallow plumbing system of the volcano since meteorological fluctuations cannot explain this observed increase of diffuse CO₂ emission.     

Major-ion bulk deposition around an active volcano (Mt. Etna, Italy)

Bulk atmospheric deposition of major cations (Na, K, Ca, Mg) and anions (Cl, F, SO₄) were measured at 15 sites around an active volcano, Mount Etna, from 2001 to 2003. Their composition indicates several natural sources, among which deposition of plume-derived volcanogenic gas compounds is prevalent for F, Cl and S. Plume-derived acidic compounds are also responsible for the prevailing acidic composition of the samples collected on the summit of the volcano (pH in the 2.45–5.57 range). Cation species have complex origin, including deposition of plume volcanogenic ash and aerosols and soil-dust wind re-suspension of either volcanic or carbonate sedimentary rocks. Variation of the deposition rates during the March 2001–March 2003 period, coupled with previous measurements from 1997 to 2000 (Appl Geochem 16:985–1000, 2001), were compared with the variation of SO₂ flux, volcanic activity and rainfall. The deposition rate was mainly controlled by rainfall. Commonly, about 0.1–0.9% of HF, HCl and SO₂ emitted by the summit crater's plume were deposited around the volcano. We estimate that ∼2 Gg of volcanogenic sulphur were deposited over the Etnean area during the 2002–2003 flank eruption, at an average rate of ∼24 Mg day⁻¹ which is two orders of magnitude higher than that typical of quiescent degassing phases.     

Catastrophic precipitation‐triggered lahar at Casita volcano,Nicaragua: occurrence,bulking and transformation

A catastrophic lahar began on 30 October 1998, as hurricane precipitation triggered a small ?ank collapse of Casita volcano, a complex and probably dormant stratovolcano. The initial rockslide‐debris avalanche evolved on the ?ank to yield a watery debris ?ood with a sediment concentration less than 60 per cent by volume at the base of the volcano. Within 2·5 km, however, the watery ?ow entrained (bulked) enough sediment to transform entirely to a debris ?ow. The debris ?ow, 6 km downstream and 1·2 km wide and 3 to 6 m deep, killed 2500 people, nearly the entire populations of the communities of El Porvenir and Rolando Rodriguez. These ‘new towns’ were developed in a prehistoric lahar pathway: at least three ?ows of similar size since 8330 ¹⁴C years bp are documented by stratigraphy in the same 30‐degree sector. Travel time between perception of the ?ow and destruction of the towns was only 2·5–3·0 minutes. The evolution of the ?ow wave occurred with hydraulic continuity and without pause or any extraordinary addition of water. The precipitation trigger of the Casita lahar emphasizes the need, in volcano hazard assessments, for including the potential for non‐eruption‐related collapse lahars with the more predictable potential of their syneruption analogues. The ?ow behaviour emphasizes that volcano collapses can yield not only volcanic debris avalanches with restricted runouts, but also mobile lahars that enlarge by bulking as they ?ow. Volumes and hence inundation areas of collapse‐runout lahars can increase greatly beyond their sources: the volume of the Casita lahar bulked to at least 2·6 times the contributing volume of the ?ank collapse and 4·2 times that of the debris ?ood. At least 78 per cent of the debris ?ow matrix (sediment < ?1·0Φ; 2 mm) was entrained during ?ow. Copyright © 2004 John Wiley & Sons, Ltd.     

Diffuse Emission of CO2 from Showa-Shinzan, Hokkaido, Japan: A Sign of Volcanic Dome Degassing

Two soil CO₂ efflux surveys were carried out in September 1999 and June 2002 to study the spatial distribution of diffuse CO₂ degassing and estimate the total CO₂ output from Showa-Shinzan volcanic dome, Japan. Seventy-six and 81 measurements of CO₂ efflux were performed in 1999 and 2002, respectively, covering most of Showa-Shinzan volcano. Soil CO₂ efflux data showed a wide range of values up to 552 g m^-2 d^-1. Carbon isotope signatures of the soil CO₂ ranged from -0.9‰ to -30.9‰, suggesting a mixing between different carbon reservoirs. Most of the study area showed CO₂ efflux background values during the 1999 and 2002 surveys (B = 8.2 and 4.4 g m^-2 d^-1, respectively). The spatial distribution of CO₂ efflux anomalies for both surveys showed a good correlation with the soil temperature, indicating a similar origin for the extensive soil degassing generated by condensation processes and fluids discharged by the fumarolic system of Showa-Shinzan. The total diffuse CO₂ output of Showa-Shinzan was estimated to be about 14.0–15.6 t d^-1 of CO₂ for an area of 0.53 km².     

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.     

Remote and in situ plume measurements of acid gas release from La Soufrière volcano, Guadeloupe

This paper presents the first remote measurements of La Soufrière gas emissions since the fumarolic and seismic reactivation in 1992. The chemical composition of the plumes has been measured from May 2003 to September 2004 using an Open Path Fourier Transform InfraRed (OP-FTIR) spectrometer, up to 15 m downwind the South Crater. HCl is clearly detected (concentration between 2.4 and 12 ppmv) whereas SO₂ and H₂S generally remain below the detection limit of the OP-FTIR. Direct measurements of SO₂ and H₂S near the South Crater with a Lancom III analyzer show a fast decrease of their concentrations with the distance. Calculated Cl / S mass ratios are high: from 9.4 ± 1.7 at 15 m from the vent to 2.8 ± 0.6 at 140 m. The enrichment in HCl of the gas emitted at La Soufrière, observed since 1998, corresponds to the degassing of a magma enriched in Cl and depleted in S. This result agrees with isotopic measurements which suggest a magmatic origin of the gases. Readjustments inside the volcanic system may have taken place during the seismic activity beginning in 1992 and enhance the transfer of magmatic gases to the summit.     

Linear spectral unmixing algorithm for modelling suspended sediment concentration of flash floods,upper Tekeze River,Ethiopia

Flash floods are the highest sediment transporting agent,but are inaccessible for in-situ sampling and have rarely been analyzed by remote sensing technology.Laboratory and field experiments were done to develop linear spectral unmixing(LSU) remote sensing model and evaluate its performance in simulating the suspended sediment concentration(SSC) in flash floods.The models were developed from continuous monitoring in the laboratory and the onsite spectral signature of river bed sediment deposits and flash floods in the Tekeze River and in its tributary,the Tsirare River.The Pearson correlation coefficient was used to determine the variability of correlations between reflectance and SSCs.The coefficient of determination(R2) and root mean square of error(RMSE) were used to evaluate the performance of the generated models.The results found that the Pearson correlation coefficient between SSCs and reflectance varied based on the level of the SSCs,geological colors,and grain sizes.The performance of the LSU model and empirical remote sensing approaches were computed to be R2=0.92,and RMSE=±0.76 g/1 in the Tsirare River and R2=0.91,and RMSE=±0.73 g/1 in the Tekeze River and R2=0.81,RMSE=±2.65 g/l in the Tsirare river and R2=0.76,RMSE=±10.87 g/l in the Tekeze River,respectively.Hence,the LSU approach of remote sensing was found to be relatively accurate in monitoring and modeling the variability of SSCs that could be applied to the upper Tekeze River basin.