Search

Kilauea Volcano's 2018 lower East Rift Zone (LERZ) eruption produced exceptionally high lava effusion rates and record-setting SO2 emissions. The eruption involved a diverse range of magmas, including primitive basalts sourced from Kilauea's summit reservoirs. We analyzed LERZ matrix glasses, melt inclusions, and host minerals to identify melt volatile contents and magma storage depths. The LERZ glasses and melt inclusions span nearly the entire compositional range previously recognized at Kilauea. Melt inclusions in Fo(86-89) olivine from the main eruptive vent (fissure 8) underwent 70-170 degrees C cooling during transport in LERZ carrier melts, causing extensive post-entrapment crystallization and sulfide precipitation. Many of these melt inclusions have low sulfur (400-900 ppm) even after correction for sulfide formation. CO2 and H2O vapor saturation pressures indicate shallow melt inclusion trapping depths (1-5 km), consistent with formation within Kilauea's Halema'uma'u and South Caldera reservoirs. Many of these inclusions also have degassed delta S-34 values (1.5 to -0.5%). Collectively, these results indicate that some primitive melts experienced near-surface degassing before being trapped into melt inclusions. We propose that decades-to-centuries of repeated lava lake activity and lava drain-back during eruptions (e.g., 1959 Kilauea Iki) recycled substantial volumes of degassed magma into Kilauea's shallow reservoir system. Degassing and magma recycling from the 2008-2018 Halema'uma'u lava lake likely reduced the volatile contents of LERZ fissure 8 magmas, resulting in lower fountain heights compared to many prior Kilauea eruptions. The eruption's extreme SO2 emissions were due to high lava effusion rates rather than particularly volatile-rich melts. U.S. Geological Survey (USGS) Volcano Science Center; Hawaiian Volcano Observatory (HVO); University of Hawaii-Hilo; HVO volunteer program; Department of Earth Sciences at the University of Oregon; Mineralogical Society of America; Geological Society of America; Mazamas student research grant program; National Science Foundation (NSF) Graduate Research Fellowship ProgramNational Science Foundation (NSF); NSF Graduate Research Internship Program (GRIP) Published version The authors would like to thank Matthew Loewen, Nicole Metrich, and Matt Patrick for constructive input that significantly improved this manuscript. The authors also thank the U.S. Geological Survey (USGS) Volcano Science Center, Hawaiian Volcano Observatory (HVO), University of Hawaii-Hilo, partner agencies, and the residents of Hawaii for support, field access, data sharing, and for their great care in documenting and responding to the 2018 LERZ eruption crisis. AHL thanks Tina Neal and the HVO volunteer program for support, Mike Zoeller for map assistance, and Carolyn Parcheta for collecting and sharing samples. Geochemical analyses were conducted with the help of John Donovan and Julie Chouinard (EPMA), and Brian Monteleone and Glenn Gaetani (SIMS and MI rehomogenization). AHL thanks Michelle Muth and Madison Myers for discussions on methodology and melt inclusion interpretation. AHL acknowledges funding support from Department of Earth Sciences at the University of Oregon, the Mineralogical Society of America, the Geological Society of America, the Mazamas student research grant program, the National Science Foundation (NSF) Graduate Research Fellowship Program, and the NSF Graduate Research Internship Program (GRIP). Coordination of GRIP at the USGS is through the Youth and Education in Science programs within the Office of Science Quality and Integrity. Public domain – authored by a U.S. government employee

We characterized volatile emissions based upon diffuse soil degassing measurements and fumarolic gas chemistry at Crater Hills, a thermally-altered area adjoining the Sour Creek resurgent dome that is located within the Yellowstone Caldera. The objective of this study was to investigate the source and flux of CO2 and H2S gases to improve our understanding of both the total emissions and origin of the spatial distribution. The total emission of CO2 estimated using the sequential Gaussian simulation method (sGs) was 66 to 109 t day−1 with 95% confidence, which is an underestimation due to the: (1) inability to measure a high flux area on a steep slope, and (2) absence of measurements from fumarole and hot pool emissions. Based on gas chemistry data obtained for a fumarole at Crater Hills in 2007, the proportion of CO2 calculated to be derived from magma would be at least 38%, but could be as high as 50%. The spatial distribution of prominent geothermal features with the highest gas flux are broadly consistent with the regional fault pattern and, therefore, likely reflect the pattern of blind faults and/or fractures covered by overlying alluvium. The estimated emission of H2S was 0.39 t day−1, based on the linear correlation between H2S and CO2. The heat output was also estimated to be ~35 MW with an average heat flux of ~100 W m−2 based upon CO2-H2O-heat relations.

Starting on 3 May 2018, a series of eruptive fissures opened in Kīlauea Volcano’s lower East Rift Zone (LERZ). Over the course of the next 3 months, intense degassing accompanied lava effusion from these fissures. Here, we report on ground-based observations of the gas emissions associated with Kīlauea’s 2018 eruption. Visual observations combined with radiative transfer modeling show that ultraviolet light could not efficiently penetrate the gas and aerosol plume in the LERZ, complicating SO2 measurements by differential optical absorption spectroscopy (DOAS). By applying a statistical method that integrates a radiative transfer model with the DOAS retrievals, we were able to calculate sulfur dioxide (SO2) emission rates along with estimates of their uncertainty. We find that sustained SO2 emissions were highest in June and early July, when approximately 200 kt SO2 were emitted daily. At the 68% confidence interval, we estimate that 7.1–13.6 Mt SO2 were released from the LERZ during the entire May to September eruptive episode. Scaling our results with in situ measurements of plume composition, we calculate that 11–21 Mt H2O and 1.5–2.8 Mt CO2 were also emitted. The gas and aerosol emissions caused hazardous conditions in areas proximal to the active vents, but plume dispersion modeling shows that the eruption also significantly impacted air quality hundreds of kilometers downwind. Combined with petrologic studies of the erupted lavas, our measurements indicate that 1.1–2.3 km3 dense-rock equivalent of lava were erupted from the LERZ, which is approximately twice the concomitant collapse volume of the volcano’s summit.

Carbon dioxide emissions from active subaerial volcanoes represent 20–50% of the annual global volcanic CO2 flux (Barry et al., 2014). Passive degassing of carbon from the flanks of volcanoes, and the associated accumulation of dissolved inorganic carbon (DIC) within nearby groundwater, also represents a potentially important, yet poorly constrained flux of carbon to the surface (Werner et al., 2019). Here we investigate sources and sinks of DIC in groundwaters in the Lassen Peak region of California. Specifically, we report and interpret the relative abundance and isotopic composition of helium (3He, 4He) and carbon (12C, 13C, 14C) in 37 groundwater samples, from 24 distinct wells, collected between 20 and 60 km from Lassen Peak. Measured groundwater samples have air-corrected 3He/4He values between 0.19 and 7.44 RA (where RA = air 3He/4He = 1.39 × 10−6), all in excess of the radiogenic production value (~0.05 RA), indicating pervasive mantle-derived helium additions to the groundwater system in the Lassen Peak region. Stable carbon isotope ratios of DIC (δ13C) vary between −12.6 and − 27.7‰ (vs. VPDB). Measured groundwater DIC/3He values fall in the range of 2.2 × 1010 to 1.1 × 1012. Using helium and carbon isotope data, we explore several conceptual models to estimate surface carbon contributions and to differentiate between DIC derived from soil CO2 versus DIC derived from external (slab and mantle) carbon sources. Specifically, if we use 14C to identify soil-derived DIC (assuming decadal-to-centennial groundwater ages and a soil CO2 14C activity equal to that of the atmosphere), we calculate that a hypothetical external carbon source would have an apparent δ13C signature between −10.3 and − 59.3‰ (vs. Vienna Pee Dee Belemnite (VPDB)) and an apparent C/3He between 7.0 × 109 and 1.0 × 1012. These apparent δ13C and C/3He values are substantially isotopically lighter than and greater than canonical MORB values, respectively. We suggest that >95% of any external (non-soil-derived) DIC in groundwater must thus be non-mantle in origin (i.e., slab derived or assimilated organic carbon). We further investigate possible sources of external DIC to groundwater using two idealized conceptual approaches: a pure (unfractionated) source mixing model (after Sano and Marty, 1995) and a scenario that invokes fractionation due to calcite precipitation. Because the former model requires carbon contributions from an organic source component with unrealistically low δ13C (~ − 60‰), we suggest that the second scenario is more plausible. Importantly, however, we caution that all conceptual models are dependent on assumptions about initial 14C activity. Thus, we cannot rule out the possibility that the true fraction of non-surface-derived DIC in these samples is lower or negligible, despite the pervasive mantle-derived He isotope signatures throughout the region. Following the 14C approach to deconvolving sources of DIC, we determine that the maximum passive carbon flux could be up to ~2.2 × 106 kg/yr, which is lower than previous magmatic carbon flux estimates from the Lassen region ( Rose and Davisson, 1996 ). We find that the passive dissolved carbon flux could represent a maximum of ~4–18% of the total Lassen geothermal CO2 degassing flux (estimated to be ~3.5 × 107 kg/yr Rose and Davisson, 1996 ; Gerlach et al., 2008 ), which is still more than an order of magnitude smaller than soil gas CO2 flux estimates (7.3–11 × 107 kg/yr) for nearby volcanoes ( Sorey et al., 1998 ; Gerlach et al., 1999 ; Evans et al., 2002 ; Werner et al., 2014 ). We conclude that passive dissolved carbon fluxes should be combined with geothermal fluxes and soil gas fluxes to obtain a complete picture of volcanic carbon emissions globally. Our approach highlights the utility of measuring helium isotopes in concert with the full suite of noble gas abundances, tritium, δ13C and 14C, which when interpreted together can be used to better elucidate the various sources of DIC in groundwater.

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. SO2 emission rates ranged from 400 to 860 t d− 1 and CO2/SO2 ratios were

Cleveland and Tana are remote volcanoes located in the central Aleutian volcanic arc on the eastern end of the Islands of Four Mountains (IFM). The persistently active Mount Cleveland volcano, on the western side of Chuginadak Island, is surrounded by several closely spaced Quaternary volcanic centers including Carlisle, Herbert, Kagamil, Tana, and Uliaga, and numerous small satellite vents on Chiginadak between Cleveland and Tana. The Alaska Volcano Observatory (AVO) installed two permanent broadband seismometers on Chuginadak Island in 2014, and we operated a temporary broadband network focused on the western side of the island in 2015–2016. Collectively, these stations provided the first seismic observations of this frequently active volcano and the surrounding Holocene-aged volcanic vents. During the study period (July 2014–January 2019), eruptive activity at Cleveland was characterized by small explosions separated by periods of lava effusion that formed small domes in the volcano's summit crater. We characterize seismicity beneath Chuginadak Island through automated analysis of event waveform frequency content, development of a one-dimensional P-wave velocity model, calculation of earthquake hypocenters, magnitudes, focal mechanisms, and identification of earthquake families. This analysis reveals the full range of seismic event types expected in a highly active volcanic environment and includes Volcano-Tectonic (VT) earthquakes, Long-Period (LP) events, and explosion signals. LP events appear to cluster at shallow depth beneath the active crater of Mount Cleveland and almost all of the explosions occur without identifiable short-term (hours to days) seismic precursors. VT earthquakes beneath Mount Cleveland occur at depths of 2 to 8 km below sea level (BSL) and range in magnitude from −0.2 to 1.8. VT focal mechanisms have horizontal P-axes that align with the regional axis of maximum stress. These observations, and a relatively slow one-dimensional seismic velocity model, are consistent with a shallow body of magma that is fed through a deeper conduit system. The time-history of VT earthquakes and shallow LP events suggest their occurrence may track the transfer of magma and fluids from the mid-crust to the shallow portions of the conduit system and may provide a means to anticipate future explosions and periods of dome growth. VT hypocenters also extend ~7 km northeast of Cleveland's summit at depths of 5 to 10 km BSL, under a group of Holocene-aged vents between Mount Cleveland and Tana. These earthquakes have vertically-oriented P-axes and a greater percentage occur in families. These observations, combined with observations of vent orientation and morphology and gas flux, suggest the area between Cleveland and Tana represents a zone of complicated volcano-tectonic interaction, similar to calderas elsewhere in the Aleutian arc. The presence of a larger volcanic system in the eastern IFM could influence magmatism and account for the multiple closely spaced volcanic centers in this region.

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 m below and 180 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 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.

Magmatic hydrothermal systems are of increasing interest as a renewable energy source. Surface heat flux indicates system resource potential, and can be inferred from soil CO 2 flux measurements and fumarole gas chemistry. Here we compile and reanalyze results from previous CO 2 flux surveys worldwide to compare heat flux from a variety of magma-hydrothermal areas. We infer that availability of water to recharge magmatic hydrothermal systems is correlated with heat flux. Recharge availability is in turn governed by permeability, structure, lithology, rainfall, topography, and perhaps unsurprisingly, proximity to a large supply of water such as the ocean. The relationship between recharge and heat flux interpreted by this study is consistent with recent numerical modeling that relates hydrothermal system heat output to rainfall catchment area. This result highlights the importance of recharge as a consideration when evaluating hydrothermal systems for electricity generation, and the utility of CO 2 flux as a resource evaluation tool.

SO 2 camera measurements at Kīlauea Volcano, Hawaii, in May of 2010 captured two occurrences of lava lake rise and fall within the Halema'um'au Crater summit vent. During high lava stands we observed diminished SO 2 emission rates and decreased seismic tremor. Similar events at Kīlauea have been described as the result of sporadic degassing following gas accumulation beneath a mostly impermeable lava lake surface. Incorporation of SO 2 camera data into a multi-parameter dataset gives credence to the likelihood of shallow gas accumulation as the cause of these high stand events, with accumulated gas release upon lake-level drop compensating for the gas deficit reached during accumulation.

SO 2 camera systems allow rapid two-dimensional imaging of sulfur dioxide (SO 2 ) emitted from volcanic vents. Here, we describe the development of an SO 2 camera system specifically designed for semi-permanent field installation and continuous use. The integration of innovative but largely “off-the-shelf” components allowed us to assemble a robust and highly customizable instrument capable of continuous, long-term deployment at Kīlauea Volcano's summit Overlook Crater. Recorded imagery is telemetered to the USGS Hawaiian Volcano Observatory (HVO) where a novel automatic retrieval algorithm derives SO 2 column densities and emission rates in real-time. Imagery and corresponding emission rates displayed in the HVO operations center and on the internal observatory website provide HVO staff with useful information for assessing the volcano's current activity. The ever-growing archive of continuous imagery and high-resolution emission rates in combination with continuous data from other monitoring techniques provides insight into shallow volcanic processes occurring at the Overlook Crater. An exemplary dataset from September 2013 is discussed in which a variation in the efficiency of shallow circulation and convection, the processes that transport volatile-rich magma to the surface of the summit lava lake, appears to have caused two distinctly different phases of lake activity and degassing. This first successful deployment of an SO 2 camera for continuous, real-time volcano monitoring shows how this versatile technique might soon be adapted and applied to monitor SO 2 degassing at other volcanoes around the world.

Mammoth Mountain, California, is a dacitic volcano that has experienced several periods of unrest since 1989. The onset of diffuse soil CO2 emissions at numerous locations on the flanks of the volcano began in 1989–1990 following an 11-month period of heightened seismicity. CO2 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 CO2 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). CO2 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 CO2-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 CO2 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 CO2 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 CO2 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).

The quantification of heat and mass flow between deep reservoirs and the surface is important for understanding magmatic and hydrothermal systems. Here, we use high-resolution measurement of carbon dioxide flux (φCO2) and heat flow at the surface to characterize the mass (CO2 and steam) and heat released to the atmosphere from two magma-hydrothermal systems. Our soil gas and heat flow surveys at Rotokawa and White Island in the Taupō Volcanic Zone, New Zealand, include over 3000 direct measurements of φCO2 and soil temperature and 60 carbon isotopic values on soil gases. Carbon dioxide flux was separated into background and magmatic/hydrothermal populations based on the measured values and isotopic characterization. Total CO2 emission rates (ΣCO2) of 441 ± 84 t d−1 and 124 ± 18 t d−1 were calculated for Rotokawa (2.9 km2) and for the crater floor at White Island (0.3 km2), respectively. The total CO2 emissions differ from previously published values by +386 t d−1 at Rotokawa and +25 t d−1 at White Island, demonstrating that earlier research underestimated emissions by 700% (Rotokawa) and 25% (White Island). These differences suggest that soil CO2 emissions facilitate more robust estimates of the thermal energy and mass flux in geothermal systems than traditional approaches. Combining the magmatic/hydrothermal-sourced CO2 emission (constrained using stable isotopes) with reservoir H2O:CO2 mass ratios and the enthalpy of evaporation, the surface expression of thermal energy release for the Rotokawa hydrothermal system (226 MWt) is 10 times greater than the White Island crater floor (22.5 MWt).

Ultraviolet (UV) camera systems represent an exciting new technology for measuring two dimensional sulfur dioxide (SO 2 ) distributions in volcanic plumes. The high frame rate of the cameras allows the retrieval of SO 2 emission rates at time scales of 1 Hz or higher, thus allowing the investigation of high-frequency signals and making integrated and comparative studies with other high-data-rate volcano monitoring techniques possible. One drawback of the technique, however, is the limited spectral information recorded by the imaging systems. Here, a framework for simulating the sensitivity of UV cameras to various SO 2 distributions is introduced. Both the wavelength-dependent transmittance of the optical imaging system and the radiative transfer in the atmosphere are modeled. The framework is then applied to study the behavior of different optical setups and used to simulate the response of these instruments to volcanic plumes containing varying SO 2 and aerosol abundances located at various distances from the sensor. Results show that UV radiative transfer in and around distant and/or optically thick plumes typically leads to a lower sensitivity to SO 2 than expected when assuming a standard Beer–Lambert absorption model. Furthermore, camera response is often non-linear in SO 2 and dependent on distance to the plume and plume aerosol optical thickness and single scatter albedo. The model results are compared with camera measurements made at Kilauea Volcano (Hawaii) and a method for integrating moderate resolution differential optical absorption spectroscopy data with UV imagery to retrieve improved SO 2 column densities is discussed.

Filter pack data from six airborne campaigns at Redoubt Volcano, Alaska are reported here. These measurements provide a rare constraint on Cl output from an andesitic eruption at high emission rate (> 104 t d− 1 SO2). Four S/Cl ratios measured during a period of lava dome growth indicate a depth of last magma equilibration of 2–5 km. The S/Cl ratios in combination with COSPEC SO2 emission rate measurements indicate HCl emission rates of 1500–3600 t d− 1 during dome growth. SO2 and HCl emission rates at Redoubt Volcano correlate with each other and were low prior to the eruption, high during the eruption, and low after the eruption. S/Cl ratios measured by filter pack at andesitic volcanoes have a small range of variance, with no clear trends seen for eruptive versus passive activity. The very few S/Cl ratio measurements by filter pack at andesitic volcanoes are not as predictive of future volcanic activity as has been demonstrated for basaltic volcanoes. This may be because there are so few of these measurements. We have demonstrated it is possible to collect these samples by air between explosions during lava dome-building eruptions. We recommend more filter pack sampling be performed at andesitic volcanoes to determine the technique's utility for volcano monitoring. Filter pack data has been demonstrated to be useful for calculating the depth of magma equilibration at volcanoes including Redoubt Volcano.

The 2009 eruption of Redoubt Volcano, Alaska, provided a rare opportunity to compare satellite measurements of sulfur dioxide (SO2) by the Ozone Monitoring Instrument (OMI) with airborne SO2 measurements by the Alaska Volcano Observatory (AVO). Herein we: (1) compare OMI and airborne SO2 column density values for Redoubt's tropospheric plume, (2) calculate daily SO2 masses from Mount Redoubt for the first three months of the eruption, (3) develop simple methods to convert daily measured SO2 masses into emission rates to allow satellite data to be directly integrated with the airborne SO2 emissions dataset, (4) calculate cumulative SO2 emissions from the eruption, and (5) evaluate OMI as a monitoring tool for high-latitude degassing volcanoes. A linear correlation (R2 ~ 0.75) is observed between OMI and airborne SO2 column densities. OMI daily SO2 masses for the sample period ranged from ~ 60.1 kt on 24 March to below detection limit, with an average daily SO2 mass of ~ 6.7 kt. The highest SO2 emissions were observed during the initial part of the explosive phase and the emissions exhibited an overall decreasing trend with time. OMI SO2 emission rates were derived using three methods and compared to airborne measurements. This comparison yields a linear correlation (R2 ~ 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 SO2 masses calculated from daily OMI data for the sample period are estimated to range from 542 to 615 kt, with approximately half of this SO2 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 SO2 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 SO2 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 SO2 emissions, and provide high temporal resolution SO2 emissions data that can be used with interdisciplinary datasets to illuminate volcanic processes.

Carbon dioxide emissions and heat flow have been determined from the Ohaaki hydrothermal field, Taupo Volcanic Zone (TVZ), New Zealand following 20 a of production (116 MW e ). Soil CO 2 degassing was quantified with 2663 CO 2 flux measurements using the accumulation chamber method, and 2563 soil temperatures were measured and converted to equivalent heat flow (W m −2 ) using published soil temperature heat flow functions. Both CO 2 flux and heat flow were analysed statistically and then modelled using 500 sequential Gaussian simulations. Forty subsoil CO 2 gas samples were also analysed for stable C isotopes. Following 20 a of production, current CO 2 emissions equated to 111 ± 6.7 T/d. Observed heat flow was 70 ± 6.4 MW, compared with a pre-production value of 122 MW. This 52 MW reduction in surface heat flow is due to production-induced drying up of all alkali–Cl outflows (61.5 MW) and steam-heated pools (8.6 MW) within the Ohaaki West thermal area (OHW). The drying up of all alkali–Cl outflows at Ohaaki means that the soil zone is now the major natural pathway of heat release from the high-temperature reservoir. On the other hand, a net gain in thermal ground heat flow of 18 MW (from 25 MW to 43.3 ± 5 MW) at OHW is associated with permeability increases resulting from surface unit fracturing by production-induced ground subsidence. The Ohaaki East (OHE) thermal area showed no change in distribution of shallow and deep soil temperature contours despite 20 a of production, with an observed heat flow of 26.7 ± 3 MW and a CO 2 emission rate of 39 ± 3 T/d. The negligible change in the thermal status of the OHE thermal area is attributed to the low permeability of the reservoir beneath this area, which has limited production (mass extraction) and sheltered the area from the pressure decline within the main reservoir. Chemistry suggests that although alkali–Cl outflows once contributed significantly to the natural surface heat flow (∼50%) they contributed little ( 2 emissions due to the loss of >99% of the original CO 2 content due to depressurisation and boiling as the fluids ascended to the surface. Consequently, the soil has persisted as the major (99%) pathway of CO 2 release to the atmosphere from the high temperature reservoir at Ohaaki. The CO 2 flux and heat flow surveys indicate that despite 20 a of production the variability in location, spatial extent and magnitude of CO 2 flux remains consistent with established geochemical and geophysical models of the Ohaaki Field. At both OHW and OHE carbon isotopic analyses of soil gas indicate a two-stage fractionation process for moderate-flux (>60 g m −2 d −1 ) sites; boiling during fluid ascent within the underlying reservoir and isotopic enrichment as CO 2 diffuses through porous media of the soil zone. For high-flux sites (>300 g m −2 d −1 ), the δ 13 CO 2 signature (−7.4 ± 0.3‰ OHW and −6.5 ± 0.6‰ OHE) is unaffected by near-surface (soil zone) fractionation processes and reflects the composition of the boiled magmatic CO 2 source for each respective upflow. Flux thresholds of −2 d −1 for purely diffusive gas transport, between 30 and 300 g m −2 d −1 for combined diffusive–advective transport, and ⩾300 g m −2 d −1 for purely advective gas transport at Ohaaki were assigned. δ 13 CO 2 values and cumulative probability plots of CO 2 flux data both identified a threshold of ∼15 g m −2 d −1 by which background (atmospheric and soil respired) CO 2 may be differentiated from hydrothermal CO 2 .

article i nfo On September 17th, 2006, Fourpeaked volcano had a widely-observed phreatic eruption. At the time, Fourpeaked was an unmonitored volcano with no known Holocene activity, based on limited field work. Airborne gas sampling began within days of the eruption and a modest seismic network was installed in stages. Vigorous steaming continued for months; however, there were no further eruptions similar in scale to the September 17 event. This eruption was followed by several months of sustained seismicity punctuated by vigorous swarms, and SO2 emissions exceeding a thousand tons/day. Based on observations during and after the phreatic eruption, and assuming no recent pre-historical eruptive activity at Fourpeaked, we propose that the activity was caused by a minor injection of new magma at or near 5 km depth beneath Fourpeaked, which remained active over several months as this magma equilibrated into the crust. By early 2007 declining seismicity and SO2 emission signaled the end of unrest. Because the Fourpeaked seismic network was installed in stages and the seismicity was punctuated by discrete swarms, we use Fourpeaked to illustrate quantitatively the efficacy and shortcomings of rapid response seismic networks for tracking volcanic earthquakes.

Mount Baker volcano displayed a short interval of seismically-quiescent thermal unrest in 1975, with high emissions of magmatic gas that slowly waned during the following three decades. The area of snow-free ground in the active crater has not returned to pre-unrest levels, and fumarole gas geochemistry shows a decreasing magmatic signature over that same interval. A relative microgravity survey revealed a substantial gravity increase in the ~30 years since the unrest, while deformation measurements suggest slight deflation of the edifice between 1981–83 and 2006–07. The volcano remains seismically quiet with regard to impulsive volcano-tectonic events, but experiences shallow ( 10 km) long-period earthquakes. Reviewing the observations from the 1975 unrest in combination with geophysical and geochemical data collected in the decades that followed, we infer that elevated gas and thermal emissions at Mount Baker in 1975 resulted from magmatic activity beneath the volcano: either the emplacement of magma at mid-crustal levels, or opening of a conduit to a deep existing source of magmatic volatiles. Decadal-timescale, multi-parameter observations were essential to this assessment of magmatic activity.

Cook Inlet volcanoes that experienced an eruption between 1989 and 2006 had mean gas emission rates that were roughly an order of magnitude higher than at volcanoes where unrest stalled. For the six events studied, mean emission rates for eruptions were ∼13,000 t/d CO2 and 5200 t/d SO2, but only ∼1200 t/d CO2 and 500 t/d SO2 for non-eruptive events (‘failed eruptions’). Statistical analysis suggests degassing thresholds for eruption on the order of 1500 and 1000 t/d for CO2 and SO2, respectively. Emission rates greater than 4000 and 2000 t/d for CO2 and SO2, respectively, almost exclusively resulted during eruptive events (the only exception being two measurements at Fourpeaked). While this analysis could suggest that unerupted magmas have lower pre-eruptive volatile contents, we favor the explanations that either the amount of magma feeding actual eruptions is larger than that driving failed eruptions, or that magmas from failed eruptions experience less decompression such that the majority of H2O remains dissolved and thus insufficient permeability is produced to release the trapped volatile phase (or both). In the majority of unrest and eruption sequences, increases in CO2 emission relative to SO2 emission were observed early in the sequence. With time, all events converged to a common molar value of C/S between 0.5 and 2. These geochemical trends argue for roughly similar decompression histories until shallow levels are reached beneath the edifice (i.e., from 20–35 to ∼4–6 km) and perhaps roughly similar initial volatile contents in all cases. Early elevated CO2 levels that we find at these high-latitude, andesitic arc volcanoes have also been observed at mid-latitude, relatively snow-free, basaltic volcanoes such as Stromboli and Etna. Typically such patterns are attributed to injection and decompression of deep (CO2-rich) magma into a shallower chamber and open system degassing prior to eruption. Here we argue that the C/S trends probably represent tapping of vapor-saturated regions with high C/S, and then gradual degassing of remaining dissolved volatiles as the magma progresses toward the surface. At these volcanoes, however, C/S is often accentuated due to early preferential scrubbing of sulfur gases. The range of equilibrium degassing is consistent with the bulk degassing of a magma with initial CO2 and S of 0.6 and 0.2 wt.%, respectively, similar to what has been suggested for primitive Redoubt magmas.

Digital cameras, sensitive to specific regions of the ultra-violet (UV) spectrum, have been employed for quantifying sulfur dioxide (SO2) emissions in recent years. The instruments make use of the selective absorption of UV light by SO2 molecules to determine pathlength concentration. Many monitoring advantages are gained by using this technique, but the accuracy and limitations have not been thoroughly investigated. The effect of some user-controlled parameters, including image exposure duration, the diameter of the lens aperture, the frequency of calibration cell imaging, and the use of the single or paired bandpass filters, have not yet been addressed. In order to clarify methodological consequences and quantify accuracy, laboratory and field experiments were conducted. Images were collected of calibration cells under varying observational conditions, and our conclusions provide guidance for enhanced image collection. Results indicate that the calibration cell response is reliably linear below 1500 ppm m, but that the response is significantly affected by changing light conditions. Exposure durations that produced maximum image digital numbers above 32 500 counts can reduce noise in plume images. Sulfur dioxide retrieval results from a coal-fired power plant plume were compared to direct sampling measurements and the results indicate that the accuracy of the UV camera retrieval method is within the range of current spectrometric methods.

Long-term changes have occurred in the chemistry, isotopic ratios, and emission rates of gas at Mount Baker volcano following a major thermal perturbation in 1975. In mid-1975 a large pulse in sulfur and carbon dioxide output was observed both in emission rates and in fumarole samples. Emission rates of CO 2 and H 2 S were ∼ 950 and 112 t/d, respectively, in 1975; these decreased to ∼ 150 and 2 /CH 4 ratio since 1975 suggest a long steady trend back toward a more hydrothermal gas signature. The helium isotope ratio is very high (> 7 R c / R A ), but has declined slightly since the mid-1970s, and δ 13 C–CO 2 has decreased by ≥ 1‰ over time. Both trends are expected from a gradually crystallizing magma. While other scenarios are investigated, we conclude that magma intruded the mid- to shallow-crust beneath Mount Baker during the thermal awakening of 1975. Since that time, evidence for fresh magma has waned, but the continued emission of CO 2 and the presence of a long-term hydrothermal system leads us to suspect some continuing connection between the surface and deep convecting magma.

We characterize and quantify volatile emissions at Hot Spring Basin (HSB), a large acid-sulfate region that lies just outside the northeastern edge of the 640 ka Yellowstone Caldera. Relative to other thermal areas in Yellowstone, HSB gases are rich in He and H2, and mildly enriched in CH4 and H2S. Gas compositions are consistent with boiling directly off a deep geothermal liquid at depth as it migrates toward the surface. This fluid, and the gases evolved from it, carries geochemical signatures of magmatic volatiles and water–rock reactions with multiple crustal sources, including limestones or quartz-rich sediments with low K/U (or 40⁎Ar/4⁎He). Variations in gas chemistry across the region reflect reservoir heterogeneity and variable degrees of boiling. Gas-geothermometer temperatures approach 300 °C and suggest that the reservoir feeding HSB is one of the hottest at Yellowstone. Diffuse CO2 flux in the western basin of HSB, as measured by accumulation-chamber methods, is similar in magnitude to other acid-sulfate areas of Yellowstone and is well correlated to shallow soil temperatures. The extrapolation of diffuse CO2 fluxes across all the thermal/altered area suggests that 410 ± 140 t d− 1 CO2 are emitted at HSB (vent emissions not included). Diffuse fluxes of H2S were measured in Yellowstone for the first time and likely exceed 2.4 t d− 1 at HSB. Comparing estimates of the total estimated diffuse H2S emission to the amount of sulfur as SO42− in streams indicates ~ 50% of the original H2S in the gas emission is lost into shallow groundwater, precipitated as native sulfur, or vented through fumaroles. We estimate the heat output of HSB as ~ 140–370 MW using CO2 as a tracer for steam condensate, but not including the contribution from fumaroles and hydrothermal vents. Overall, the diffuse heat and volatile fluxes of HSB are as great as some active volcanoes, but they are a small fraction (1–3% for CO2, 2–8% for heat) of that estimated for the entire Yellowstone system.

The 2004–2008 dome-building eruption at Mount St. Helens ended during winter 2007–2008 at a time when field observations were hampered by persistent bad weather. As a result, recognizing the end of the eruption was challenging—but important for scientists trying to understand how and why long-lived eruptions end and for public officials and land managers responsible for hazards mitigation and access restrictions. In hindsight, the end of the eruption was presaged by a slight increase in seismicity in December 2007 that culminated on January 12–13, 2008, with a burst of more than 500 events, most of which occurred in association with several tremor-like signals and a spasmodic burst of long-period earthquakes. At about the same time, a series of regular, localized, small-amplitude tilt events—thousands of which had been recorded during earlier phases of the eruption—came to an end. Thereafter, seismicity declined to 10–20 events per day until January 27–28, when a spasmodic burst of about 50 volcano-tectonic earthquakes occurred over a span of 3 h. This was followed by a brief return of repetitive “drumbeat” earthquakes that characterized much of the eruption. By January 31, however, seismicity had declined to 1–2 earthquakes per day, a rate similar to pre-eruption levels. We attribute the tilt and seismic observations to convulsive stagnation of a semisolid magma plug in the upper part of the conduit. The upward movement of the plug ceased when the excess driving pressure, which had gradually decreased throughout the eruption as a result of reservoir deflation and increasing overburden from the growing dome, was overcome by increasing friction as a result of cooling and crystallization of the plug.

We present the first routine measurements of emissions (CO 2 , SO 2 , and H 2 S) from Ruapehu volcano, New Zealand during a crater lake heating cycle. Emissions were generally at a low level consistent with quiescent degassing of the volcano and the presence of a crater lake. Maximum concentrations (∼ 3–4 km downwind) during the highest measurement in May, 2004 were 14 ppm, 0.37 ppm, and 7 ppb for CO 2 , SO 2 and H 2 S, respectively. CO 2 emissions varied between not detectable to 900 t d − 1 over periods of months. SO 2 was not detectable until February, but gradually increased to 35 t d − 1 . Emissions of H 2 S were detected in April and May, although were − 1 . The CO 2 /SO 2 ratio in the plume was ∼ 37 by weight for each measurement during the peak of the cycle suggesting that significant scrubbing of SO 2 occurs through the crater lake, and that a common source exists for both gases. Magma volumes estimated from CO 2 emissions (0.001–0.004 km 3 ) are consistent with eruptive volumes given repose periods of 20–30 years. Delays between peaks in crater-lake heating and degassing suggests that volcanic emissions do not primarily reside in the shallow hydrothermal system directly beneath the crater. Further data is needed to adequately model the system, but first indications of the travel time associated with the degassing cycle suggests that the gas resides at minimally 300 m depth beneath the crater lake. Depths greater than 300 m are consistent with the top of the single-phase steam zone in heat pipe models (200–700 m depth) or perhaps even below the plastic–brittle transition zone above the cooling magmas (> 1 km deep).

The first measurements of volcanic/hydrothermal water vapor and heat flux using eddy covariance (EC) were made at Solfatara crater, Italy, June 8–25, 2001. 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 of water vapor varied between 680 and 11 200 g H2O m− 2 d− 1. Heat fluxes varied diurnally with the solar input, and the volcanic component of sensible heat ranged from ∼25 to 238 W m− 2. The highest measurements of both sensible and latent heat flux were made downwind of hot soil regions and degassing pools and during mid-day. The ratio of average volcanic heat (both latent and sensible) to CO2 flux resulted in an equivalent H2O/CO2 flux ratio of 2.2 by weight, which reflects the deep source H2O/CO2 gas ratio. The amount latent heat flux/evaporation was determined to be consistent both with what would be expected from the magnitude of CO2 fluxes and the fumarolic H2O/CO2 ratio, as well as with observed surface temperatures and wind speeds given a moist soil. This suggests that the water vapor that condenses in the shallow subsurface is remobilized at the soil–atmosphere interface through variable evaporation dependent on the deep heat flux and surface temperature. The results suggest that EC provides a quick and easy method to monitor average H2O/CO2 ratios continuously in volcanic regions, providing another important tool for volcanic hazards monitoring.

An eddy covariance (EC) station was deployed at Solfatara crater, Italy, June 8–25, 2001 to assess if EC could reliably monitor CO2 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 CO2 m−2 d−1; 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 CO2 d−1 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−1. Daily average EC fluxes varied by ±50% and variations were likely correlated to changes in atmospheric pressure. Variations in CO2 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.

The 2008-current summit eruption at Kīlauea Volcano, Hawai'i offers a unique opportunity to test models of degassing and magma plumbing and to improve our understanding of the volatile budget. The aim of this work was to test the hypothesis that gases emitted from a summit lava lake will be rich in carbon dioxide (CO) and similar to those measured during the persistent lava lake activity in the early 20th century at Kīlauea Volcano (Gerlach and Graeber, 1985). We measured the sulfur dioxide (SO) and CO concentrations in the gas plume from Halema'uma'u using electrochemical and non-dispersive infrared sensors during April 2009. We also analysed olivine-hosted melt inclusions from tephra erupted in 2008 and 2010 for major, trace and volatile elements. The gas and melt data are both consistent with the equilibration of a relatively evolved magma batch at depths of 1.2-2.0 km beneath Halema'uma'u prior to the current degassing activity. The differences in the volatile concentrations between the melt inclusions and matrix glasses are consistent with the observed gas composition. The degassing of sulfur and halogen gases from the melt requires low pressures and hence we invoke convection to bring the magma close to the surface to degas, before sinking back into the conduit. The fluxes of gases (900 and 80 t/d SO and CO respectively) are used to estimate magma fluxes (1.2-3.4 m/s) to the surface for April 2009. The observation of minimal loss of hydrogen from the melt inclusions implies a rapid rise rate (less than a few hours), which constrains the conduit radius to 1-2 m. The inferred conduit radius is much narrower than the lava lake at the surface, implying a flared geometry. The melt inclusion data suggest that there is a progressive decrease in melt volatile concentrations with time during 2008-2010, consistent with convection, degassing and mixing in a closed, or semi-closed magma system. The degassing regime of the current summit lava lake activity is not similar to that observed in the early 20th century; instead the gases are extensively depleted in CO. © 2013 The Authors.

Statistical sampling theory offers several sampling designs that are well-suited to geophysical research. Simple random sampling is the most commonly used formal methodology, in which the probability that any unit is sampled is known. Generally, units to be sampled are selected using a random number generator so that each unit has an equal chance of selection. In stratified sampling, the sampling frame is partitioned so that similar observed variable measurements are in each stratum. Systematic sampling attempts to capture the full variability within an area by obtaining samples at some predetermined interval. Adaptive cluster sampling is a relatively new methodology that allows additional units to be sampled, based on preceding sampling observations. The systematic adaptive sampling method applies adaptive cluster sampling to a systematic design. Selecting an appropriate design requires an understanding of the variability in the parameter of interest as well as the spatial covariance structure. This paper addresses such critical issues as sampling density, number of samples to collect, and assessment of the dimensions of a subregion to be sampled. While many researchers use the sampling designs presented herein, it is possible that inappropriate estimating equations were used. Formulae required to calculate unbiased, precise estimators for each of the designs are presented. The collection of CO2 flux degassing measurements from the Mud Volcano area, Yellowstone during the 1997 field season is used as a case study to outline the logical steps in selecting a sampling design. Although the case study used a stratified adaptive cluster sampling design, simpler designs can be developed based on the theory presented herein.

A stratified adaptive sampling plan was designed to estimate CO2 degassing in Yellowstone National Park and was applied in the Mud Volcano thermal area. The stratified component of the sampling design focused effort in thermal areas and the adaptive component in high-flux regions, yet neither sampling technique biased the estimate of total degassing. Both diffuse soil fluxes (up to ∼30,000 g m−2d−1) and emission rates from thermal vents (up to 1.7×108 mol yr−1) were measured in thermal areas. Soil fluxes observed in most nonthermal regions were similar to values reported for conifer forests (≤ 15 g m−2 d−1). However, through adaptive sampling, high-flux vegetated sites were identified in Mud Volcano that likely would not have been found if sampling was focused in obvious thermal or altered regions. A simple model applied to flux measurements suggests that ∼40% of the analyzed measurements were dominated by possible advective transport and ∼30% by diffusive transport. Isotopic signatures of soil CO2 generally suggest a deep origin (δ13C = −2.3 to 0.0) in thermal areas and biogenic origin (δ13C = −20.5) in nonthermal, low-flux areas. Vent emissions accounted for ∼32–63% of the total degassing observed at Mud Volcano (2.4 to 4.0×109 mol yr−1). The largest source of error in the estimation of total degassing (factor of ∼2) resulted because the population distribution of thermal feature emissions was indeterminate. Total CO2 emissions at Mud Volcano are comparable to other hydrothermal regions worldwide, suggesting that the Yellowstone volcanic system is likely a large contributor to global volcanic/metamoohic/hydrothermal (VMH) emissions.

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.

[1] Sulfur dioxide (SO2) 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 Kīlauea 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 SO2 and aerosol loads present in the volcanic plume. The internal consistency of the method was validated by measuring and modeling SO2absorption 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 SO2 plumes.

the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 � 0.4 MW and 7.5 � 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 � 0.3 MW and 1.2 � 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 � 1.4 MW. These estimates result in a total heat output of 11.8 � 1.4 MW and 8.8 � 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. ��

[1] We report CO2, SO2, and H2S emission rates and C/S ratios during the five months leading up to the 2009 eruption of Redoubt Volcano, Alaska. CO2emission 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 SO2, not the >2100 t/d SO2expected 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 CO2-rich magmatic gas for at least 5 months prior to eruption, and minimal scrubbing of SO2by near-surface groundwater. The high C/S ratios observed could reflect bulk degassing of mid-crustal magma followed by nearly complete loss of SO2in 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.

[1] We report on 4 years of airborne measurements of CO2, SO2, and H2S emission rates during a quiescent period at White Island volcano, New Zealand, beginning in 2003. During this time a significant crater lake emerged, allowing scrubbing processes to be investigated. CO2 emissions varied from a baseline of 250 to >2000 t d−1 and demonstrated clear annual cycling that was consistent with numbers of earthquake detections and annual changes in sea level. The annual variability was found to be most likely related to increases in the strain on the volcano during sea level highs, temporarily causing fractures to reduce in size in the upper conduit. SO2 emissions varied from 0 to >400 t d−1 and were clearly affected by scrubbing processes within the first year of lake development. Scrubbing caused increases of SO42− and Cl− in lake waters, and the ratio of carbon to total sulphur suggested that elemental sulphur deposition was also significant in the lake during the first year. Careful measurements of the lake level and chemistry allowed estimates of the rate of H2O(g) and HCl(g) input into the lake and suggested that the molar abundances of major gas species (H2O, CO2, SO2, and HCl) during this quiescent phase were similar to fumarolic ratios observed between earlier eruptive periods. The volume of magma estimated from CO2 emissions (0.015–0.04 km3) was validated by Cl− increases in the lake, suggesting that the gas and magma are transported from deep to shallow depths as a closed system and likely become open in the upper conduit region. The absence of surface deformation further leads to a necessity of magma convection to supply and remove magma from the degassing depths. Two models of convection configurations are discussed.

The diurnal and vertical variability of heat and carbon dioxide (CO2) in the atmospheric surface layer are studied by analyzing measurements from a 213 m tower in Cabauw (Netherlands). Observations of thermodynamic variables and CO2 mixing ratio as well as vertical profiles of the turbulent fluxes are used to retrieve the contribution of the budget terms in the scalar conservation equation. On the basis of the daytime evolution of turbulent fluxes, we calculate the budget terms by assuming that turbulent fluxes follow a linear profile with height. This assumption is carefully tested and the deviation from linearity is quantified. The budget calculation allows us to assess the importance of advection of heat and CO2 during day hours for three selected days. It is found that, under nonadvective conditions, the diurnal variability of temperature and CO2 is well reproduced from the flux divergence measurements. Consequently, the vertical transport due to the turbulent flux plays a major role in the daytime evolution of both scalars and the advection is a relatively small contribution. During the analyzed days with a strong contribution of advection of either heat or carbon dioxide, the flux divergence is still an important contribution to the budget. For heat, the quantification of the advection contribution is in close agreement with results from a numerical model. For carbon dioxide, we qualitatively corroborate the results with a Lagrangian transport model. Our estimation of advection is compared with traditional estimations based on the Net Ecosystem-atmosphere Exchange (NEE).

[1] We present a reliable methodology to estimate the energy associated with the subaerial diffuse degassing of volcanic-hydrothermal fluids. The fumaroles of 15 diffuse degassing structures (DDSs) located in eight volcanic systems in the world were sampled and analyzed. Furthermore, each area was measured for soil temperature gradients and for soil CO2 fluxes. The results show that each hydrothermal or volcanic system is characterized by a typical source fluid which feeds both the fumaroles and diffuse degassing through the soil. Experimental data and the results of physical numerical modeling of the process demonstrate that the heat released by condensation of steam at depth is almost totally transferred by conduction in the uppermost part of the soil. A linear relationship is observed between the log of the steam/gas ratio measured in the fumaroles and the log of the ratio between soil thermal gradient and soil-gas flux. The main parameter controlling this relation is the thermal conductivity of the soil (Kc). For each area, we computed the values of Kc which range from 0.4 to 2.3 W m−1 °C−1. Using the CO2 soil fluxes as a tracer of the deep fluids, we estimated that the total heat released by steam condensation in the systems considered varies from 1 to 100 MW.

We present a comparative study of soil CO2 flux (\(F_{{\rm CO}_2 }\)) measured by five groups (Groups 1–5) at the IAVCEI-CCVG Eighth Workshop on Volcanic Gases on Masaya volcano, Nicaragua. Groups 1–5 measured \(F_{{\rm CO}_2 }\) using the accumulation chamber method at 5-m spacing within a 900 m2 grid during a morning (AM) period. These measurements were repeated by Groups 1–3 during an afternoon (PM) period. Measured \(F_{{\rm CO}_2 }\) ranged from 218 to 14,719 g m−2 day−1. The variability of the five measurements made at each grid point ranged from ±5 to 167%. However, the arithmetic means of fluxes measured over the entire grid and associated total CO2 emission rate estimates varied between groups by only ±22%. All three groups that made PM measurements reported an 8–19% increase in total emissions over the AM results. Based on a comparison of measurements made during AM and PM times, we argue that this change is due in large part to natural temporal variability of gas flow, rather than to measurement error. In order to estimate the mean and associated CO2 emission rate of one data set and to map the spatial \(F_{{\rm CO}_2 }\) distribution, we compared six geostatistical methods: arithmetic and minimum variance unbiased estimator means of uninterpolated data, and arithmetic means of data interpolated by the multiquadric radial basis function, ordinary kriging, multi-Gaussian kriging, and sequential Gaussian simulation methods. While the total CO2 emission rates estimated using the different techniques only varied by ±4.4%, the \(F_{{\rm CO}_2 }\) maps showed important differences. We suggest that the sequential Gaussian simulation method yields the most realistic representation of the spatial distribution of \(F_{{\rm CO}_2 }\), but a variety of geostatistical methods are appropriate to estimate the total CO2 emission rate from a study area, which is a primary goal in volcano monitoring research.

[1] Two methods are used to estimate CO2 degassing from the Yellowstone magmatic-hydrothermal system. The amount of magmatic CO2 released as basaltic magma emplaces from the mantle into the crust beneath the Yellowstone caldera is calculated and compared to CO2 fluxes measured in three different types of hydrothermal regions within Yellowstone. Comparison of modeled estimates with surface measurements suggests that 3.7 ± 1.3 × 1011 mol y−1 (45 ± 16 kt d−1) of CO2 are released from Yellowstone due to diffuse degassing. Flux measurements suggest that the diffuse flux in acid-sulfate regions is significant in total calculations (>96% of the total), whereas the diffuse flux in neutral-chloride and travertine-precipitating areas is not significant. Analyses of carbon and helium isotopes suggest that ∼50% of the CO2 emitted is derived from sedimentary sources at locations outside the caldera, whereas locations inside the caldera likely have sedimentary contributions