Geology Reports⌕ Search

Geology topics

Peter J. Kelly

Publications and source records attributed to Peter J. Kelly.

33 records · Page 2Linked to original sources

Insights into the mechanisms of phreatic eruptions from continuous high frequency volcanic gas monitoring: Rincón de la Vieja volcano, Costa Rica

Understanding the trigger mechanisms of phreatic eruptions is key to mitigating the effects of these hazardous but poorly forecastable volcanic events. It has recently been established that high-rate volcanic gas observations are potentially very suitable to identifying the source processes driving phreatic eruptions, and to eventually detecting precursory changes prior to individual phreatic blasts. In February-May 2017, we deployed a Multi-GAS instrument to continuously monitor gas concentrations in the crater lake plume of Rincón de la Vieja, a remote and poorly monitored active volcano in Costa Rica, site of frequent phreatic/phreatomagmatic eruptions. Forty-two phreatic/phreatomagmatic eruptions were seismically recorded during our investigated period, 9 of which were also recorded for gas by the Multi-GAS. To the best of our knowledge, these represent the first instrumentally measured gas compositions during individual phreatic/phreatomagmatic explosions at an active volcano. Our results show that during background quiescent degassing the Rincón de la Vieja crater lake plume was characterized by high CO 2 /SO 2 ratios of 64 ± 59 and H 2 S/SO 2 ratios of 0.57 ± 0.20. This composition is interpreted as reflecting hydrothermal (re)processing of magma-sourced gas in the sub-limnic environment. Phreatic blasts were recorded by the Multi-GAS as brief (1–2 min long) pulses of elevated gas mixing ratios (up to ~52 ppmv SO 2 and >3,000 ppmv CO 2 ), or more than an order of magnitude higher than during background degassing (~1 ppmv SO 2 and ~450 ppmv CO 2 ). During the phreatic eruption(s), the H 2 S/SO 2 ratio was systematically lower (<0.18) than during background degassing, but the CO 2 /SO 2 ratio remained high (and variable), ranging from 37 to 390. These S-poor compositions for the eruptive gas imply extensive processing of the source magmatic gas during pre-eruptive hydrothermal storage, likely by deposition of native S and/or sulfate. Our gas results are thus overall consistent with a mechanism of phreatic eruptions triggered by accumulation of magmatic-hydrothermal gases beneath a hydrothermal seal. We claim that real-time Multi-GAS monitoring is urgently needed at other crater lake-hosting volcanoes (e.g., Ruapehu, Aso), where phreatic eruptions may similarly be preceded by phases of reduced S degassing at the surface.

Rincón de la Vieja Volcano↗

Construction of probabilistic event trees for eruption forecasting at Sinabung volcano, Indonesia 2013–14

Eruptions of Sinabung volcano , Indonesia have been ongoing since 2013. Since that time, the character of eruptions has changed, from phreatic to phreatomagmatic to magmatic explosive eruptions, and from production of a lava dome that collapsed to a subsequent thick lava flow that slowly ceased to be active, and later, to a new lava dome. As the eruption progressed, event trees were constructed to forecast eruptive behavior six times, with forecast windows that ranged from 2 weeks to 1 year: November 7–10, December 12–14, and December 27, 2013; and January 9–10, May 13, and October 7, 2014. These event trees were successful in helping to frame the forecast scenarios, to collate current monitoring information, and to document outstanding questions and unknowns. The highest probability forecasts closely matched outcomes of eruption size (including extrusion of the first dome), production of pyroclastic density currents, and pyroclastic density current runout distances. Events assigned low probabilities also occurred, including total collapse of the lava dome in January 2014 and production of a small blast pyroclastic density current in February 2014.

Sinabung volcano↗

The 2018 rift eruption and summit collapse of Kilauea Volcano

In 2018, Kīlauea Volcano experienced its largest lower East Rift Zone (LERZ) eruption and caldera collapse in at least 200 years. After collapse of the Pu'u 'Ō'ō vent on 30 April, magma propagated downrift. Eruptive fissures opened in the LERZ on 3 May, eventually extending ~6.8 km. A 4 May earthquake (M6.9) produced ~5 m of fault slip. Lava erupted at rates exceeding 100 m3/s, eventually covering 35.5 km2. The summit magma system partially drained, producing minor explosions and near-daily collapses releasing energy equivalent to M4.7-M5.4 earthquakes. Activity declined rapidly on 4 August. Summit collapse and lava flow volume estimates are roughly equivalent--about 0.8 km3. Careful historical observation and monitoring of Kīlauea enabled successful forecasting of hazardous events.

Hawaii↗

Evolution of the 2015 Cotopaxi eruption revealed by combined geochemical & seismic observations

Through integration of multiple data streams to monitor volcanic unrest scientists are able to make more robust eruption forecast and to obtain a more holistic interpretation of volcanic systems. We examined gas emission and gas geochemistry, seismic and petrologic data recorded during the 2015 unrest of Cotopaxi (Ecuador) in order to decipher the origin and temporal evolution of this eruption. Identification of families of similar seismic events and the use of seismic amplitude ratios reveals temporal changes in volcanic processes. SO2 (300 to 24000 t/d), BrO/SO2 (5-10 x10-5), SO2/HCl (5.8 ± 4.8 and 6.6 ± 3.0) and CO2/SO2 (0.6 to 2.1) measured throughout the eruption indicate a shallow magmatic source. Bulk ash and glass chemistry indicate a homogenous andesitic (SiO2 wt%=56.94 ± 0.25) magma having undergone extensive S-exsolution and degassing during ascent. These data lead us to interpret this eruption as a magma intrusion and ascend to shallow levels. The intrusion progressively interacted with the hydrothermal system, boiled off water, and produced hydromagmatic explosions. A small volume of this intrusion continued to fragment and produced episodic ash emissions until it was sufficiently degassed and rheologically stiff. Based on the 470 kt of measured SO2 we estimate that ~ 65.3 x106 m3 of magma were required to supply the emitted gases. This volume exceeds the volume of erupted juvenile material by a factor of 50. This result emphasizes the importance of careful monitoring of Cotopaxi to identify the intrusion of a new batch of magma, which could rejuvenate the non-erupted material.

Geochemistry, Geophysics, Geosystems↗

Multi-year high-frequency hydrothermal monitoring of selected high-threat Cascade Range volcanoes

From 2009 to 2015 the U.S. Geological Survey (USGS) systematically monitored hydrothermal behavior at selected Cascade Range volcanoes in order to define baseline hydrothermal and geochemical conditions. Gas and water data were collected regularly at 25 sites on 10 of the highest-risk volcanoes in the Cascade Range. These sites include near-summit fumarole groups and springs/streams that show clear evidence of magmatic influence (high 3 He/ 4 He ratios and/or large fluxes of magmatic CO 2 or heat). Site records consist mainly of hourly temperature and hydrothermal-flux data. Having established baseline conditions during a multiyear quiescent period, the USGS reduced monitoring frequency from 2015 to present. The archived monitoring data are housed at (doi:10.5066/F72N5088). These data (1) are suitable for retrospective comparison with other continuous geophysical monitoring data and (2) will provide context during future episodes of volcanic unrest, such that unrest-related variations at these thoroughly characterized sites will be more clearly recognizable. Relatively high-frequency year-round data are essential to achieve these objectives, because many of the time series reveal significant diurnal, seasonal, and inter-annual variability that would tend to mask unrest signals in the absence of baseline data. Here we characterize normal variability for each site, suggest strategies to detect future volcanic unrest, and explore deviations from background associated with recent unrest.

California, Oregon, Washington↗

Magmatic degassing, lava dome extrusion, and explosions from Mount Cleveland volcano, Alaska, 2011–2015: Insight into the continuous nature of volcanic activity over multi-year timescales

Mount Cleveland volcano (1730 m) is one of the most active volcanoes in the Aleutian arc, Alaska, but heightened activity is rarely accompanied by geophysical signals, which makes interpretation of the activity difficult. In this study, we combine volcanic gas emissions measured for the first time in August 2015 with longer-term measurements of thermal output and lava extrusion rates between 2011 and 2015 calculated from MODIS satellite data with the aim to develop a better understanding of the nature of volcanic activity at Mount Cleveland. Degassing measurements were made in the month following two explosive events (21 July and 7 August 2015) and during a period of new dome growth in the summit crater. SO 2 emission rates ranged from 400 to 860 t d − 1 and CO 2 /SO 2 ratios were < 3, consistent with the presence of shallow magma in the conduit and the observed growth of a new lava dome. Thermal anomalies derived from MODIS data from 2011 to 2015 had an average repose time of only 4 days, pointing to the continuous nature of volcanic activity at this volcano. Rapid increases in the cumulative thermal output were often coincident with visual confirmation of dome growth or accumulations of tephra in the crater. The average rate of lava extrusion calculated for 9 periods of rapid increase in thermal output was 0.28 m 3 s − 1 , and the total volume extruded from 2011 to 2015 was 1.9–5.8 Mm 3 . The thermal output from the lava extrusion events only accounts for roughly half of the thermal budget, suggesting a continued presence of shallow magma in the upper conduit, likely driven by convection. Axisymmetric dome morphology and occasional drain back of lava into the conduit suggests low-viscosity magmas drive volcanism at Mount Cleveland. It follows also that only small overpressures can be maintained given the small domes and fluid magmas, which is consistent with the low explosivity of most of Mount Cleveland's eruptions. Changes between phases of dome growth and explosive activity are somewhat unpredictable and likely result from plugs that are related to the dome obtaining a critical dimension, or from small variations in the magma ascent rate that lead to crystallization-induced blockages in the upper conduit, thereby reducing the ability of magma to degas. We suggest the small magma volumes, slow ascent rates, and low magma viscosity lead to the overall lack of anomalous geophysical signals prior to eruptions, and that more continuous volcanic degassing measurements might lead to more successful eruption forecasting at this continuously-active open-vent volcano.

Alaska↗

Evidence for degassing of fresh magma during the 2004-2008 eruption of Mount St. Helens: Subtle signals from the hydrothermal system

Results from chemical and isotopic analyses of water and gas collected between 2002 and 2016 from sites on and around Mount St. Helens are used to assess magmatic degassing related to the 2004-2008 eruption. During 2005 the chemistry of hot springs in The Breach of Mount St. Helens showed no obvious response to the eruption, and over the next few years, changes were subtle, giving only slight indications of perturbations in the system. By 2010 however, water chemistry, temperatures, and isotope compositions (δD and δ 18 O) clearly indicated some inputs of volatiles and heat associated with the eruption, but the changes were such that they could be attributed to a pre-existing, gas depleted magma. An increase of ~ 1.5‰ in the δ 13 C values of dissolved carbon in the springs was noted in 2006 and continued through 2009, a change that was mirrored by a similar shift in δ 13 C-CO 2 in bubble gas emissions. These changes require input of a new source of carbon to the hydrothermal system and provide clear evidence of CO 2 from an undegassed body of magma. Rising trends in 3 He/ 4 He ratios in gas also accompanied the increases in δ 13 C. Since 2011 maximum R C /R A values are ≥ 6.4 and are distinctly higher than 5 samples collected between 1986 and 2002, and provide additional evidence for some involvement of new magma as early as 2006, and possibly earlier, given the unknown time needed for CO 2 and He to traverse the system and arrive at the springs.

Washington↗

Monitoring gas and heat emissions at Norris Geyser Basin, Yellowstone National Park, USA based on a combined eddy covariance and Multi-GAS approach

We quantified gas and heat emissions in an acid-sulfate, vapor-dominated area (0.04-km 2 ) of Norris Geyser Basin, located just north of the 0.63 Ma Yellowstone Caldera and near an area of anomalous uplift. From 14 May to 3 October 2016, an eddy covariance system measured half-hourly CO 2 , H 2 O and sensible ( H ) and latent ( LE ) heat fluxes and a Multi-GAS instrument measured (1 Hz frequency) atmospheric H 2 O, CO 2 and H 2 S volumetric mixing ratios. We also measured soil CO 2 fluxes using the accumulation chamber method and temperature profiles on a grid and collected fumarole gas samples for geochemical analysis. Eddy covariance CO 2 fluxes ranged from − 56 to 885 g m − 2 d − 1 . Using wavelet analysis, average daily eddy covariance CO 2 fluxes were locally correlated with average daily environmental parameters on several-day to monthly time scales. Estimates of CO 2 emission rate from the study area ranged from 8.6 t d − 1 based on eddy covariance measurements to 9.8 t d − 1 based on accumulation chamber measurements. Eddy covariance water vapor fluxes ranged from 1178 to 24,600 g m − 2 d − 1 . Nighttime H and LE were considered representative of hydrothermal heat fluxes and ranged from 4 to 183 and 38 to 504 W m − 2 , respectively. The total hydrothermal heat emission rate ( H + LE + radiant) estimated for the study area was 11.6 MW and LE contributed 69% of the output. The mean ± standard deviation of H 2 O, CO 2 and H 2 S mixing ratios measured by the Multi-GAS system were 9.3 ± 3.1 parts per thousand, 467 ± 61 ppmv, and 0.5 ± 0.6 ppmv, respectively, and variations in the gas compositions were strongly correlated with diurnal variations in environmental parameters (wind speed and direction, atmospheric temperature). After removing ambient H 2 O and CO 2 , the observed variations in the Multi-GAS data could be explained by the mixing of relatively H 2 O-CO 2 -H 2 S-rich fumarole gases with CO 2 -rich and H 2 O-H 2 S-poor soil gases. The fumarole H 2 O/CO 2 and CO 2 /H 2 S end member ratios (101.7 and 27.1, respectively, on average) were invariant during the measurement period and fell within the range of values measured in direct fumarole gas samples. The soil gas H 2 O/CO 2 end member ratios (~ 15–30) were variable and low relative to the fumarole end member, likely resulting from water vapor loss during cooling and condensation in the shallow subsurface, whereas the CO 2 /H 2 S end member ratio was high (~ 160), presumably related to transport of CO 2 -dominated soil gas emissions mixed with trace fumarolic emissions to the Multi-GAS station. Nighttime eddy covariance ratios of H 2 O to CO 2 flux were typically between the soil gas and fumarole end member H 2 O/CO 2 ratios defined by Multi-GAS measurements. Overall, the combined eddy covariance and Multi-GAS approach provides a powerful tool for quasi-continuous measurements of gas and heat emissions for improved volcano-hydrothermal monitoring.

Wyoming↗

Long period seismicity and very long period infrasound driven by shallow magmatic degassing at Mount Pagan, Mariana Islands

Long period (LP) seismicity and very long period infrasound (iVLP) were recorded during continuous degassing from Mount Pagan, Mariana Islands, in July 2013 to January 2014. The frequency content of the LP and iVLP events and delay times between the two arrivals were remarkably stable and indicate nearly co-located sources. Using phase-weighted stacking over similar events to dampen noise, we find that the LP source centroid is located 60&thinsp;m below and 180&thinsp;m west of the summit vent. The moment tensor reveals a volumetric source modeled as resonance of a subhorizontal sill intersecting a dike. We model the seismoacoustic wavefields with a coupled earth-air 3-D finite difference code. The ratios of pressure to velocity measured at the infrasound arrays are an order of magnitude larger than the synthetic ratios, so the iVLP is not the result of LP energy transmitting into the atmosphere at its epicenter. Based on crater shape and dimensions determined by structure from motion, we model the iVLP as acoustic resonance of an exponential horn. The source of the continuous plume from gas analysis is shallow magmatic degassing, which repeatedly pressurized the dike-sill portion of the conduit over the 7&thinsp;months of observation. Periodic gas release caused the geologically controlled sill to partially collapse and resonate, while venting of gas at the surface triggered resonance in the crater. LP degassing only accounts for ~12% of total degassing, indicating that most degassing is relatively aseismic and that multiple active pathways exist beneath the vent.

Mount Pagan↗

Decadal-scale variability of diffuse CO2 emissions and seismicity revealed from long-term monitoring (1995–2013) at Mammoth Mountain, California, USA

Mammoth Mountain, California, is a dacitic volcano that has experienced several periods of unrest since 1989. The onset of diffuse soil CO 2 emissions at numerous locations on the flanks of the volcano began in 1989–1990 following an 11-month period of heightened seismicity. CO 2 emission rates were measured yearly from 1995 to 2013 at Horseshoe Lake (HSL), the largest tree kill area on Mammoth Mountain, and measured intermittently at four smaller degassing areas around Mammoth from 2006 to 2013. The long-term record at HSL shows decadal-scale variations in CO 2 emissions with two peaks in 2000–2001 and 2011–2012, both of which follow peaks in seismicity by 2–3 years. Between 2000 and 2004 emissions gradually declined during a seismically quiet period, and from 2004 to 2009 were steady at ~ 100 metric tonnes per day (t d − 1 ). CO 2 emissions at the four smaller tree-kill areas also increased by factors of 2–3 between 2006 and 2011–2012, demonstrating a mountain-wide increase in degassing. Delays between the peaks in seismicity and degassing have been observed at other volcanic and hydrothermal areas worldwide, and are thought to result from an injection of deep CO 2 -rich fluid into shallow subsurface reservoirs causing a pressurization event with a delayed transport to the surface. Such processes are consistent with previous studies at Mammoth, and here we highlight (1) the mountain-wide response, (2) the characteristic delay of 2–3 years, and (3) the roughly decadal reoccurrence interval for such behavior. Our best estimate of total CO 2 degassing from Mammoth Mountain was 416 t d − 1 in 2011 during the peak of emissions, over half of which was emitted from HSL. The cumulative release of CO 2 between 1995 and 2013 from diffuse emissions is estimated to be ~ 2–3 Mt, and extrapolation back to 1989 gives ~ 4.8 Mt. This amount of CO 2 release is similar to that produced by the mid-sized (VEI 3) 2009 eruption of Redoubt Volcano in Alaska (~ 2.3 Mt over 11 months), and significantly lower than long-term emissions from hydrothermal areas such as Solfatara in Campi Flegrei, Italy (16 Mt over 28 years).

California↗

Evaluation of Redoubt Volcano's sulfur dioxide emissions by the Ozone Monitoring Instrument

The 2009 eruption of Redoubt Volcano, Alaska, provided a rare opportunity to compare satellite measurements of sulfur dioxide (SO 2 ) by the Ozone Monitoring Instrument (OMI) with airborne SO 2 measurements by the Alaska Volcano Observatory (AVO). Herein we: (1) compare OMI and airborne SO 2 column density values for Redoubt's tropospheric plume, (2) calculate daily SO 2 masses from Mount Redoubt for the first three months of the eruption, (3) develop simple methods to convert daily measured SO 2 masses into emission rates to allow satellite data to be directly integrated with the airborne SO 2 emissions dataset, (4) calculate cumulative SO 2 emissions from the eruption, and (5) evaluate OMI as a monitoring tool for high-latitude degassing volcanoes. A linear correlation (R 2 ~ 0.75) is observed between OMI and airborne SO 2 column densities. OMI daily SO 2 masses for the sample period ranged from ~ 60.1 kt on 24 March to below detection limit, with an average daily SO 2 mass of ~ 6.7 kt. The highest SO 2 emissions were observed during the initial part of the explosive phase and the emissions exhibited an overall decreasing trend with time. OMI SO 2 emission rates were derived using three methods and compared to airborne measurements. This comparison yields a linear correlation (R 2 ~ 0.82) with OMI-derived emission rates consistently lower than airborne measurements. The comparison results suggest that OMI's detection limit for high latitude, springtime conditions varies from ~ 2000 to 4000 t/d. Cumulative SO 2 masses calculated from daily OMI data for the sample period are estimated to range from 542 to 615 kt, with approximately half of this SO 2 produced during the explosive phase of the eruption. These cumulative masses are similar in magnitude to those estimated for the 1989–90 Redoubt eruption. Strong correlations between daily OMI SO 2 mass and both tephra mass and acoustic energy during the explosive phase of the eruption suggest that OMI data may be used to infer relative eruption size and explosivity. Further, when used in conjunction with complementary datasets, OMI daily SO 2 masses may be used to help distinguish explosive from effusive activity and identify changes in lava extrusion rates. The results of this study suggest that OMI is a useful volcano monitoring tool to complement airborne measurements, capture explosive SO 2 emissions, and provide high temporal resolution SO 2 emissions data that can be used with interdisciplinary datasets to illuminate volcanic processes.

Alaska↗

Degassing of CO2, SO2, and H2S associated with the 2009 eruption of Redoubt Volcano, Alaska

The 2009 eruption of Redoubt Volcano, Alaska was particularly well monitored for volcanic gas emissions. We report 35 airborne measurements of CO 2 , SO 2 , and H 2 S emission rates that span from October 2008 to August 2010. The magmatic system degassed primarily as a closed system although minor amounts of open system degassing were observed in the 6 months prior to eruption on March 15, 2009 and over 1 year following cessation of dome extrusion. Only 14% of the total CO 2 was emitted prior to eruption even though high emissions rates (between 3630 and 9020 t/d) were observed in the final 6 weeks preceding the eruption. A minor amount of the total SO 2 was observed prior to eruption (4%), which was consistent with the low emission rates at that time (up to 180 t/d). The amount of the gas emitted during the explosive and dome growth period (March 15–July 1, 2009) was 59 and 66% of the total CO 2 and SO 2 , respectively. Maximum emission rates were 33,110 t/d CO 2 , 16,650 t/d SO 2 , and 1230 t/d H 2 S. Post-eruptive passive degassing was responsible for 27 and 30% of the total CO 2 and SO 2 , respectively. SO 2 made up on average 92% of the total sulfur degassing throughout the eruption. Magmas were vapor saturated with a C- and S-rich volatile phase, and regardless of composition, the magmas appear to be buffered by a volatile composition with a molar CO 2 /SO 2 ratio of ~ 2.4. Primary volatile contents calculated from degassing and erupted magma volumes range from 0.9 to 2.1 wt.% CO 2 and 0.27–0.56 wt.% S; whole-rock normalized values are slightly lower (0.8–1.7 wt.% CO 2 and 0.22–0.47 wt.% S) and are similar to what was calculated for the 1989–90 eruption of Redoubt. Such contents argue that primary arc magmas are rich in CO 2 and S. Similar trends between volumes of estimated degassed magma and observed erupted magma during the eruptive period point to primary volatile contents of 1.25 wt.% CO 2 and 0.35 wt.% S. Assuming these values, up to 30% additional unerupted magma degassed in the year following final dome emplacement.

Alaska↗

Rapid chemical evolution of tropospheric volcanic emissions from Redoubt Volcano, Alaska, based on observations of ozone and halogen-containing gases

We report results from an observational and modeling study of reactive chemistry in the tropospheric plume emitted by Redoubt Volcano, Alaska. Our measurements include the first observations of Br and I degassing from an Alaskan volcano, the first study of O 3 evolution in a volcanic plume, as well as the first detection of BrO in the plume of a passively degassing Alaskan volcano. This study also represents the first detailed spatially-resolved comparison of measured and modeled O 3 depletion in a volcanic plume. The composition of the plume was measured on June 20, 2010 using base-treated filter packs (for F, Cl, Br, I, and S) at the crater rim and by an instrumented fixed-wing aircraft on June 21 and August 19, 2010. The aircraft was used to track the chemical evolution of the plume up to ~ 30 km downwind (2 h plume travel time) from the volcano and was equipped to make in situ observations of O 3 , water vapor, CO 2 , SO 2 , and H 2 S during both flights plus remote spectroscopic observations of SO 2 and BrO on the August 19th flight. The airborne data from June 21 reveal rapid chemical O 3 destruction in the plume as well as the strong influence chemical heterogeneity in background air had on plume composition. Spectroscopic retrievals from airborne traverses made under the plume on August 19 show that BrO was present ~ 6 km downwind (20 min plume travel time) and in situ measurements revealed several ppbv of O 3 loss near the center of the plume at a similar location downwind. Simulations with the PlumeChem model reproduce the timing and magnitude of the observed O 3 deficits and suggest that autocatalytic release of reactive bromine and in-plume formation of BrO were primarily responsible for the observed O 3 destruction in the plume. The measurements are therefore in general agreement with recent model studies of reactive halogen formation in volcanic plumes, but also show that field studies must pay close attention to variations in the composition of ambient air entrained into volcanic plumes in order to unambiguously attribute observed O 3 anomalies to specific chemical or dynamic processes. Our results suggest that volcanic eruptions in Alaska are sources of reactive halogen species to the subarctic troposphere.

Alaska↗

Improving the accuracy of S0 2 column densities and emission rates obtained from upward-looking UV-spectroscopic measurements of volcanic plumes by taking realistic radiative transfer into account

Sulfur dioxide (SO 2 ) is monitored using ultraviolet (UV) absorption spectroscopy at numerous volcanoes around the world due to its importance as a measure of volcanic activity and a tracer for other gaseous species. Recent studies have shown that failure to take realistic radiative transfer into account during the spectral retrieval of the collected data often leads to large errors in the calculated emission rates. Here, the framework for a new evaluation method which couples a radiative transfer model to the spectral retrieval is described. In it, absorption spectra are simulated, and atmospheric parameters are iteratively updated in the model until a best match to the measurement data is achieved. The evaluation algorithm is applied to two example Differential Optical Absorption Spectroscopy (DOAS) measurements conducted at Kilauea volcano (Hawaii). The resulting emission rates were 20 and 90% higher than those obtained with a conventional DOAS retrieval performed between 305 and 315 nm, respectively, depending on the different SO 2 and aerosol loads present in the volcanic plume. The internal consistency of the method was validated by measuring and modeling SO2 absorption features in a separate wavelength region around 375 nm and comparing the results. Although additional information about the measurement geometry and atmospheric conditions is needed in addition to the acquired spectral data, this method for the first time provides a means of taking realistic three-dimensional radiative transfer into account when analyzing UV-spectral absorption measurements of volcanic SO 2 plumes.

Journal of Geophysical Research D: Atmospheres↗

Deep magmatic degassing versus scrubbing: Elevated CO2 emissions and C/S in the lead-up to the 2009 eruption of Redoubt Volcano, Alaska

We report CO 2, SO 2 , and H 2 S emission rates and C/S ratios during the five months leading up to the 2009 eruption of Redoubt Volcano, Alaska. CO 2 emission rates up to 9018 t/d and C/S ratios ≥30 measured in the months prior to the eruption were critical for fully informed forecasting efforts. Observations of ice-melt rates, meltwater discharge, and water chemistry suggest that surface waters represented drainage from surficial, perched reservoirs of condensed magmatic steam and glacial meltwater. These fluids scrubbed only a few hundred tonnes/day of SO 2 , not the >2100 t/d SO 2 expected from degassing of magma in the mid- to upper crust (3–6.5 km), where petrologic analysis shows the final magmatic equilibration occurred. All data are consistent with upflow of a CO 2 -rich magmatic gas for at least 5 months prior to eruption, and minimal scrubbing of SO 2 by near-surface groundwater. The high C/S ratios observed could reflect bulk degassing of mid-crustal magma followed by nearly complete loss of SO 2 in a deep magmatic-hydrothermal system. Alternatively, high C/S ratios could be attributed to decompressional degassing of low silica andesitic magma that intruded into the mid-crust in the 5 months prior to eruption, thereby mobilizing the pre-existing high silica andesite magma or mush in this region. The latter scenario is supported by several lines of evidence, including deep long-period earthquakes (−28 to −32 km) prior to and during the eruption, and far-field deformation following the onset of eruptive activity.

Alaska↗