Geology topics

Jennifer L. Lewicki

Publications and source records attributed to Jennifer L. Lewicki.

At least 19 recordsLinked to original sources

USGS2025

Gas chemistry and isotope data for volcano monitoring at the Lassen Volcanic Center, Lassen Volcanic National Park

This report presents chemical and isotopic compositions of volcanic gases collected from thermal areas within Lassen Volcanic National Park in northern California from 1974 through 2019. As the southernmost volcano in the Cascade Range and designated a very-high-threat volcano by the U.S. Geological Survey, the Lassen Volcanic Center (LVC) requires consistent monitoring to assess potential volcanic hazards. In 2014, the California Volcano Observatory established a gas geochemical monitoring program at LVC to provide baseline data to evaluate future changes.</p> Results demonstrate consistent spatial patterns in bulk gas chemistry that support a two-circulation-cell hydrothermal model previously established for LVC. Gas samples from circulation cell 1 thermal areas have higher helium isotope ratios (6.59–7.50 times the air value) than those from circulation cell 2 (5.86–6.52 times the air value), indicating a stronger magmatic signature. The Sulphur Works and Pilot Pinnacle thermal areas within circulation cell 1 consistently emit gases with the highest magmatic helium contents, suggesting gas at these areas best represents conditions in the underlying volcanic system. A slight decrease in helium isotope values since 1974 may indicate progressive dilution of magmatic helium-3 (3He) by radiogenic helium-4 (4He) in the absence of recent magma intrusion. Carbon isotope compositions of carbon dioxide across all thermal areas are relatively uniform (−9.7–−7.3 per mil), falling within the range observed at other Cascade Range volcanoes. Based on gas geochemical characteristics and site accessibility, the Sulphur Works and Pilot Pinnacle thermal areas represent optimal targets for continued monitoring of the LVC magmatic-hydrothermal system. This study includes the most comprehensive helium isotope dataset collected at LVC currently available and establishes critical baseline data for future volcanic monitoring efforts.</p>

Open-File Report↗

USGS2025

Gas emissions from the Sulphur Bank Mercury Mine hydrothermal system, Clear Lake volcanic field, California

The Sulphur Bank Mercury Mine (SBMM) hydrothermal system offers insights into active degassing processes in the Clear Lake volcanic field (CLVF), a high-threat region based on its record of Holocene eruptions and proximity to populated areas. Here we present chemical and isotopic analyses of gas samples collected between 2015 and 2023, along with the first comprehensive CO 2 flux survey of the SBMM area conducted in 2023. Sampled gases are CO 2 - and CH 4 -rich (≥84 and 6 mol% in dry gas, respectively) with high mantle-derived helium contributions ( 3 He/ 4 He = 6.54–7.86 R C /R A ). Carbon isotopic compositions of CO 2 (δ 13 C = −10.0 to −9.5 ‰) and CH 4 (δ 13 C = −35.8 ‰) indicate mixed sources, with significant contributions from metamorphism of organic-rich Franciscan Complex rocks hosting the hydrothermal system. Modeling of gas compositions shows that scrubbing by interaction with air-saturated groundwater strongly influences observed compositional variability. From our CO₂ flux measurements, we estimate the deeply derived CO 2 emission rate from the SBMM hydrothermal area (0.2 km 2 ) at 240 t d −1 , comparable to many quiescently degassing volcanoes worldwide. We also provide a first-order estimate of CH 4 emissions at approximately 0.5 t d −1 . Our findings establish crucial baseline data for future volcanic monitoring efforts, enhancing detection capabilities for potential changes in this active hydrothermal system. This work contributes to the broader understanding of volatile contributions from volcanic and metamorphic sources to the global carbon budget, while highlighting the strong influence of bedrock geology on gas compositions in the CLVF.

California↗

USGS2024

Volcanic gas monitoring

Introduction As magma rises through the crust, decreasing pressure conditions allow volatiles to exsolve from the magma. These volatiles then migrate upward through the crust, where they can be stored at shallower levels or escape to the atmosphere. Rising magma also heats rock masses beneath volcanic centers, causing water in shallow aquifers and hydrothermal systems to boil and release additional gases and steam (see chapter F , this volume; Ingebritsen and Hurwitz, 2024). The chemistry and quantity of gases that reach the surface during periods of quiescence or volcanic unrest can reveal that gas-rich magma is ascending, crystallizing, or alternatively stalling, with important implications for volcanic hazard (for example, Sutton and others, 1992; Aiuppa and others, 2007, 2021; Werner and others, 2009, 2011, 2012; Moretti and others, 2013; de Moor and others, 2016; Lewicki and others, 2019; Edmonds and others, 2022; Kern and others, 2022; Kunrat and others, 2022). Most volcanoes in Alaska and the western United States are characterized by weak degassing, with one or more low-temperature fumaroles (typically near the local boiling temperature of water) and connect to a deeper and sometimes extensive hydrothermal system (for example, McGee and others, 2001; Symonds and others, 2003a, b). Hydrothermal systems will affect the chemistry of rising gases exsolved from deeper magma (Symonds and others, 2001), including sulfur dioxide (SO 2 ), hydrogen chloride (HCl), and water vapor (for example, Doukas and Gerlach, 1995; Gerlach and others, 1998, 2008; Symonds and others, 2001; Werner and others, 2013). As an example, depending on factors such as temperature, pressure, and oxidation state, rising SO 2 will react with groundwater to form hydrogen sulfide (H 2 S) gas, dissolved sulfate (SO 4 2− ), or elemental sulfur (Christenson, 2000; Symonds and others, 2001; Werner and others, 2008). The reaction and dissolution of SO 2 into shallow groundwater is commonly referred to as scrubbing, and can reduce the likelihood that ascending, degassing magma can be detected. Carbon dioxide, however, in addition to exsolving from magma early in the ascent process, is not easily removed by hydrothermal fluids (Lowenstern, 2001). As scrubbing and other processes take place, the SO 2 /H 2 S, CO 2 /SO 2 , and CO 2 /H 2 S ratios may change. High rates of SO 2 emission indicate that magma has moved to relatively shallow levels in the volcano and that the system has heated up enough to establish dry pathways from depth to the surface. Monitoring multiple gas species and the total output of those species is thereby useful for volcano monitoring during both periods of quiescence, to establish background degassing conditions, and during unrest, when gas geochemistry and emission rates can provide information on changing conditions, such as magma ascent. To provide context for multidisciplinary volcano forecasts, we focus on the following two key required capabilities: (1) characterizing baseline geochemistry and gas discharge from volcanoes and volcanic regions and (2) monitoring changes in gas geochemistry and discharge to inform forecasts of volcanic eruptions and their effects. Sufficient baseline data must be collected to identify and interpret anomalous degassing associated with volcanic unrest (for example, Sorey and others, 1998; Rouwet and others, 2014). Differences in volcano type, baseline degassing rates, local hydrology, and geography (for example, high versus low latitude) will result in a different baseline for each volcano. Volcanoes of any threat level that exhibit one or more degassing phenomena would ideally be monitored by techniques needed to establish baseline degassing data, with the sampling frequency of baseline data dictated by the threat level (table E1). Additional monitoring techniques become necessary during periods of unrest. In general, three of the most important techniques for gas monitoring are (1) direct sampling of fumarole, spring, and soil gases for laboratory geochemical measurements, (2) measurements of the chemical composition of the volcanic plume and emission rates of major gas species (for example, H 2 O, CO 2 , SO 2 , and H 2 S) by satellite, airborne, or ground-based techniques, and (3) measurements of diffuse emissions of CO 2 and other gases through soils. Various methods and instruments may be useful both for baseline studies and during unrest.

Scientific Investigations Report↗

USGS2023

Out of the frying pan and into the fire: Effects of volcanic heat and other stressors on the conservation of a critically endangered plant in Hawaiʻi

Loss of local biodiversity resulting from abrupt environmental change is a significant environmental problem throughout the world. Extinctions of plants are particularly important yet are often overlooked. Drawing from a case in Hawai‘i, a global hotspot for plant and other extinctions, we demonstrate an effort to better understand and determine priorities for the management of an endangered plant (‘Ihi makole or Portulaca sclerocarpa ) in the face of rapid and extreme environmental change. Volcanic heat emissions and biological invasions have anecdotally been suggested as possible threats to the species. We integrated P. sclerocarpa outplanting with efforts to collect geological and ecological data to gauge the role of elevated soil temperatures and invasive grasses in driving P. sclerocarpa mortality and population decline. We measured soil temperature, soil depth, surrounding cover and P. sclerocarpa survivorship over three decades. The abundance of wild P. sclerocarpa decreased by 99.7% from the 1990s to 2021. Only 51% of outplantings persisted through 3–4 years. Binomial regression and structural equation modelling revealed that, among the variables we analysed, high soil temperatures were most strongly associated with population decline. Finding the niche where soil temperatures are low enough to allow P. sclerocarpa survival but high enough to limit other agents of P. sclerocarpa mortality may be necessary to increase population growth of this species.

Hawaii↗

USGS2023

Geochemistry and fluxes of gases from hydrothermal features at Newberry Volcano, Oregon, USA

We present the chemical and isotopic compositions of gases and fluxes of CO 2 from the hydrothermal features of Newberry Volcano, a large composite volcano located in Oregon's Cascade Range with a summit caldera that hosts two lakes, Paulina and East Lakes. Gas samples were collected from 1982 to 2021 from Paulina Hot Springs (PHS) on the shore of Paulina Lake, East Lake Hot Springs (ELHS) on the shore of East Lake, and Obsidian Flow Gas Seep (OFGS), an area of diffuse gas emissions. Surveys of CO 2 flux were conducted in 2020 at OFGS (1400 m 2 ) and East Lake (4.1 km 2 ). Gases from all three sites were CO 2 -rich (≥79 mol% in dry gas) but showed considerable compositional variability over time due to interaction with ground and surface water. An increase in H 2 S concentrations and decline in CO 2 /H 2 S ratios in ELHS gases coincided with a drop in East Lake water level from 1999 to 2021. ELHS and OFGS gases were high in CH 4 relative to PHS and the δ 13 C of CH 4 values for ELHS gases (−72.2 and − 63.6 ‰) reflected a predominantly biogenic origin. The dominant source of N 2 and Ar in PHS, ELHS, and OFGS samples was likely groundwater. Helium isotopic ratios (6.47 to 8.02 R c /R a ) support a persistent source of magmatic He beneath Newberry caldera and consistently high values measured at OFGS and PHS relative to ELHS suggest distinct fluid flow paths from depth to the surface features. The δ 13 C of CO 2 and CO 2 / 3 He values (−8.9 to −5.35 ‰ and 1.3 × 10 9 to 4.6 × 10 10 , respectively) measured in gases reflect contributions of CO 2 from both mantle and crustal sources. Measured CO 2 fluxes at OFGS and East Lake ranged from 1 to 8808 and < 1 to 364 g m −2 d −1 , respectively. A CO 2 emission rate of 0.5 t d −1 was calculated for OFGS. The CO 2 emission rate estimated for East Lake was 30 t d −1 and when compared to prior estimates, reflects steady-state lake degassing. An enhanced geochemical monitoring plan, including annual sampling of gases at ELHS, OFGS, and PHS for geochemical analysis, installation of a continuous lake-level monitoring station at East Lake, and annual CO 2 flux surveys at OFGS, would provide valuable background data and insights into any precursor volcanic activity. Integrating geochemical data with data from the real-time seismic and GPS network at Newberry Volcano could better resolve and interpret potential changes in its magma-hydrothermal system.

Oregon↗

USGS2022

Seasonal and multi-year changes in CO2 degassing at Mammoth Mountain explained by solid-earth-driven fault valving

Changes in CO 2 emissions from volcanoes may evidence volcanic unrest. We use a multiyear time series of CO 2 flux collected at the Horseshoe Lake Tree Kill area on Mammoth Mountain, CA, to understand processes that cause variations in flux from this system. Seasonal variations are systematically lowest during the winter months and reach maximum values during the summer season. A persistent ∼20% reduction in CO 2 flux occurred during the Spring of 2017, coincident with the emergence of the area from drought and earthquake swarms in Long Valley Caldera. We used continuous GNSS measurements to calculate seasonal strains and stresses across the Mammoth Mountain area, and resolved resultant stresses onto the Mammoth Mountain Fault, which appears to facilitate gas transport to the surface. The normal stress changes are consistent with seasonal and multiyear changes in CO 2 flux, suggesting that fault valving by solid earth processes can alter surface gas fluxes.

California↗

USGS2021

Long-term year-round observations of magmatic CO2 emissions on Mammoth Mountain, California, USA

Diffuse emission of magmatic CO 2 is one of the main indicators of volcanic unrest at Mammoth Mountain, but the presence of deep seasonal snowpack at the site has hindered year-round CO 2 flux observations. A permanent eddy covariance station was established at the largest area of diffuse CO 2 degassing on Mammoth Mountain (Horseshoe Lake tree kill) that measured CO 2 fluxes ( F c ) and meteorological parameters on a half-hourly basis. From July 22, 2014 to May 24, 2020, F c ranged from −35 to 10,546 g m −2 d −1 . F c decreased on average by 53% over the study period, tracking the long-term decline in CO 2 emissions following the last major increase that occurred at the Horseshoe Lake tree kill area from 2009 to 2011. Statistical and spectral analyses were applied to the F c and ancillary meteorological parameter time series to understand (1) relationships between these parameters, (2) their dominant periodicities, and (3) changes in F c that may be unexplained by meteorological forcing. Variations in detrended F c ( F cdt ) were most strongly correlated with wind direction and atmospheric temperature, followed by atmospheric pressure on diurnal to annual time scales, but wind direction likely exerted the most direct control on F cdt . Comparison of the smoothed (180-d span) F cdt time series to the time series of average-daily snow water equivalent measured ~1 km away suggested that snowpack may have suppressed CO 2 emissions. No evidence of a change in CO 2 emissions related to the last major seismic swarm beneath Mammoth Mountain on February 2–18, 2014 was observed.

California↗

USGS2020

High-resolution imaging of hydrothermal heat flux using optical and thermal Structure-from-Motion photogrammetry

Quantifying hydrothermal heat flux at meter-scale resolution over N0.25 km2 is required to bridge in-situ heat flux and satellite-based measurements. We advance a methodology that blends ground-based daytime optical and nighttime thermal infrared (TIR) imagery using Structure-from-Motion photogrammetry to map radiant hydrothermal heat flux over these scales at sites with low signal-to-noise ratios that would otherwise be difficult to characterize using, for example, unmanned aerial systems. The improved method uses a computerized telescopic mount to relocate and align daytime optical acquisitions with nighttime TIR imagery, thereby enabling TIR acquisition from multiple camera orientations positioned throughout a study area. This facilitates mapping of thermal features at sites of varying size and complexity and helps to ameliorate topographic occlusion effects and geometric distortions that can bias radiant hydrothermal heat flux estimates derived from the resulting orthorectified thermal maps. We assessed detection thresholds of this method at three sites across central California, which range in size, topography, and heat flux conditions. We found that blending of optical and thermal acquisitions successfully detected anomalous heat flux, even in cases where temperatures were only slightly greater than background. This approach might be applied to a variety of volcanic and hydrothermal systems to quantify the spatial distribution of heat flux, and how this may relate to factors such as the distribution of ground fractures and lava flow rheology.

Journal of Volcanology and Geothermal Research↗

USGS2019

Rate of magma supply beneath Mammoth Mountain, California based on helium isotopes and CO2 emissions

Mammoth Mountain, California, has exhibited unrest over the past ~30 years, characterized by seismicity over a broad range of depths, elevated 3 He/ 4 He ratios in fumarolic gas, and large-scale diffuse CO 2 emissions. This activity has been attributed to magmatic intrusion, but minimal ground deformation and the presence of a shallow crustal gas reservoir beneath Mammoth Mountain pose a challenge for estimating magma supply rate. Here, we use the record of fumarolic 3 He/ 4 He ratios and CO 2 emissions to estimate that of the ~5.2 Mt of CO 2 released from Mammoth Mountain between 1989 and 2016, 1.6 Mt was associated with active intrusion and degassing of ~0.05–0.07 km 3 of basaltic magma. Intrusion at an average rate of ~0.002–0.003 km 3 /year into a postulated zone of partial melt at ~15-km depth could occur without detection by local Global Navigation Satellite System stations.

California↗

USGS2019

Heat and mass transport in a vapor-dominated hydrothermal area in Yellowstone National Park, USA: Inferences from magnetic, electrical, electromagnetic, subsurface temperature and diffuse CO2 flux measurements

Vapor‐dominated hydrothermal systems are characterized by localized and elevated heat and gas flux. In these systems, steam and gas ascend from a boiling water reservoir, steam condenses beneath a low‐permeability cap layer, and liquid water descends, driven by gravity (“heat pipe” model). We combine magnetic, electromagnetic, and geoelectrical methods and CO 2 flux and subsurface temperature measurements in the Solfatara Plateau Thermal Area in the Yellowstone Caldera to address several fundamental questions: (1) What are the structural and/or lithological controls on heat and mass transport in vapor‐dominated areas? (2) What is the geometry and size of convecting multiphase thermal plumes? (3) Are thermal plumes associated with subsurface rock alteration and demagnetization? Magnetic and electromagnetic data inversions suggest an asymmetric 50‐ to 100‐m thick basin of glacial deposits with the thickest part adjacent to the margin of a rhyolite flow. The 3‐D electrical conductivity model in the glacial basin reveals a narrow vertical conductor interpreted as a focused multiphase plume, which coincides at the ground surface with the heat and CO 2 flux maxima. The magnetic data suggest that destruction of magnetic minerals due to rock alteration associated with the hydrothermal plume occurs mainly near the ground surface. We propose a model where the buoyant multiphase plume forms in response to decompression, boiling, and phase separation of pressurized thermal groundwater that discharges from the brecciated base of a rhyolite flow into the basin of glacial deposits. Results from multiphase groundwater flow and heat transport numerical simulations corroborate the first‐order characteristics of this model.

Wyoming↗

USGS2018

Ecosystem responses to elevated CO2 using airborne remote sensing at Mammoth Mountain, California

We present an exploratory study examining the use of airborne remote-sensing observations to detect ecological responses to elevated CO 2 emissions from active volcanic systems. To evaluate these ecosystem responses, existing spectroscopic, thermal, and lidar data acquired over forest ecosystems on Mammoth Mountain volcano, California, were exploited, along with in situ measurements of persistent volcanic soil CO 2 fluxes. The elevated CO 2 response was used to statistically model ecosystem structure, composition, and function, evaluated via data products including biomass, plant foliar traits and vegetation indices, and evapotranspiration (ET). Using regression ensemble models, we found that soil CO 2 flux was a significant predictor for ecological variables, including canopy greenness (normalized vegetation difference index, NDVI), canopy nitrogen, ET, and biomass. With increasing CO 2 , we found a decrease in ET and an increase in canopy nitrogen, both consistent with theory, suggesting more water- and nutrient-use-efficient canopies. However, we also observed a decrease in NDVI with increasing CO 2 (a mean NDVI of 0.27 at 200 g m −2 d −1 CO 2 reduced to a mean NDVI of 0.10 at 800 g m −2 d −1 CO 2 ). This is inconsistent with theory though consistent with increased efficiency of fewer leaves. We found a decrease in above-ground biomass with increasing CO 2 , also inconsistent with theory, but we did also find a decrease in biomass variance, pointing to a long-term homogenization of structure with elevated CO 2 . Additionally, the relationships between ecological variables changed with elevated CO 2 , suggesting a shift in coupling/decoupling among ecosystem structure, composition, and function synergies. For example, ET and biomass were significantly correlated for areas without elevated CO 2 flux but decoupled with elevated CO 2 flux. This study demonstrates that (a) volcanic systems show great potential as a means to study the properties of ecosystems and their responses to elevated CO 2 emissions and (b) these ecosystem responses are measurable using a suite of airborne remotely sensed data.

California↗

USGS2018

Unraveling the dynamics of magmatic CO2 degassing at Mammoth Mountain, California

The accumulation of magmatic CO 2 beneath low-permeability barriers may lead to the formation of CO 2 -rich gas reservoirs within volcanic systems. Such accumulation is often evidenced by high surface CO 2 emissions that fluctuate over time. The temporal variability in surface degassing is believed in part to reflect a complex interplay between deep magmatic degassing and the permeability of degassing pathways. A better understanding of the dynamics of CO 2 degassing is required to improve monitoring and hazards mitigation in these systems. Owing to the availability of long-term records of CO 2 emissions rates and seismicity, Mammoth Mountain in California constitutes an ideal site towards such predictive understanding. Mammoth Mountain is characterized by intense soil CO 2 degassing (up to ∼1000 t d −1 ) and tree kill areas that resulted from leakage of CO 2 from a CO 2 -rich gas reservoir located in the upper ∼4 km. The release of CO 2 -rich fluids from deeper basaltic intrusions towards the reservoir induces seismicity and potentially reactivates faults connecting the reservoir to the surface. While this conceptual model is well-accepted, there is still a debate whether temporally variable surface CO 2 fluxes directly reflect degassing of intrusions or variations in fault permeability. Here, we report the first large-scale numerical model of fluid and heat transport for Mammoth Mountain. We discuss processes (i) leading to the initial formation of the CO 2 -rich gas reservoir prior to the occurrence of high surface CO 2 degassing rates and (ii) controlling current CO 2 degassing at the surface. Although the modeling settings are site-specific, the key mechanisms discussed in this study are likely at play at other volcanic systems hosting CO 2 -rich gas reservoirs. In particular, our model results illustrate the role of convection in stripping a CO 2 -rich gas phase from a rising hydrothermal fluid and leading to an accumulation of a large mass of CO 2 (∼10 7 –10 8 t) in a shallow gas reservoir. Moreover, we show that both, short-lived (months to years) and long-lived (hundreds of years) events of magmatic fluid injection can lead to critical pressures within the reservoir and potentially trigger fault reactivation. Our sensitivity analysis suggests that observed temporal fluctuations in surface degassing are only indirectly controlled by variations in magmatic degassing and are mainly the result of temporally variable fault permeability. Finally, we suggest that long-term CO 2 emission monitoring, seismic tomography and coupled thermal–hydraulic–mechanical modeling are important for CO 2 -related hazard mitigation.

California↗

USGS2017

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↗

USGS2016

High spatio-temporal resolution observations of crater-lake temperatures at Kawah Ijen volcano, East Java, Indonesia

The crater lake of Kawah Ijen volcano, East Java, Indonesia, has displayed large and rapid changes in temperature at point locations during periods of unrest, but measurement techniques employed to-date have not resolved how the lake’s thermal regime has evolved over both space and time. We applied a novel approach for mapping and monitoring variations in crater-lake apparent surface (“skin”) temperatures at high spatial (~32 cm) and temporal (every two minutes) resolution at Kawah Ijen on 18 September 2014. We used a ground-based FLIR T650sc camera with digital and thermal infrared (TIR) sensors from the crater rim to collect (1) a set of visible imagery around the crater during the daytime and (2) a time series of co-located visible and TIR imagery at one location from pre-dawn to daytime. We processed daytime visible imagery with the Structure-from-Motion photogrammetric method to create a digital elevation model onto which the time series of TIR imagery was orthorectified and georeferenced. Lake apparent skin temperatures typically ranged from ~21 to 33oC. At two locations, apparent skin temperatures were ~ 4 and 7 oC less than in-situ lake temperature measurements at 1.5 and 5 m depth, respectively. These differences, as well as the large spatio-temporal variations observed in skin temperatures, were likely largely associated with atmospheric effects such as evaporative cooling of the lake surface and infrared absorption by water vapor and SO2. Calculations based on orthorectified TIR imagery thus yielded underestimates of volcanic heat fluxes into the lake, whereas volcanic heat fluxes estimated based on in-situ temperature measurements (68 to 111 MW) were likely more representative of Kawah Ijen in a quiescent state. The ground-based imaging technique should provide a valuable tool to continuously monitor crater-lake temperatures and contribute insight into the spatio-temporal evolution of these temperatures associated with volcanic activity.

Kawah Ijen volcano↗

USGS2015

Integrated thermal infrared imaging and Structure-from-Motion photogrametry to map apparent temperature and radiant hydrothermal heat flux at Mammoth Mountain, CA USA

This work presents a method to create high-resolution (cm-scale) orthorectified and georeferenced maps of apparent surface temperature and radiant hydrothermal heat flux and estimate the radiant hydrothermal heat emission rate from a study area. A ground-based thermal infrared (TIR) camera was used to collect (1) a set of overlapping and offset visible imagery around the study area during the daytime and (2) time series of co-located visible and TIR imagery at one or more sites within the study area from pre-dawn to daytime. Daytime visible imagery was processed using the Structure-from-Motion photogrammetric method to create a digital elevation model onto which pre-dawn TIR imagery was orthorectified and georeferenced. Three-dimensional maps of apparent surface temperature and radiant hydrothermal heat flux were then visualized and analyzed from various computer platforms (e.g., Google Earth, ArcGIS). We demonstrate this method at the Mammoth Mountain fumarole area on Mammoth Mountain, CA. Time-averaged apparent surface temperatures and radiant hydrothermal heat fluxes were observed up to 73.7 oC and 450 W m-2, respectively, while the estimated radiant hydrothermal heat emission rate from the area was 1.54 kW. Results should provide a basis for monitoring potential volcanic unrest and mitigating hydrothermal heat-related hazards on the volcano.

California↗

USGS2014

Crustal migration of CO2-rich magmatic fluids recorded by tree-ring radiocarbon and seismicity at Mammoth Mountain, CA, USA

Unrest at Mammoth Mountain over the past several decades, manifest by seismicity, ground deformation, diffuse CO 2 emissions, and elevated 3 He/ 4 He ratios in fumarolic gases has been driven by the release of CO 2 -rich fluids from basaltic intrusions in the middle to lower crust. Recent unrest included the occurrence of three lower-crustal (32–19 km depth) seismic swarms beneath Mammoth Mountain in 2006, 2008 and 2009 that were consistently followed by peaks in the occurrence rate of shallow (≤10 km depth) earthquakes. We measured 14 C in the growth rings (1998–2012) of a tree growing in the largest (∼0.3 km 2 ) area of diffuse CO 2 emissions on Mammoth Mountain (the Horseshoe Lake tree kill; HLTK) and applied atmospheric CO 2 concentration source area modeling to confirm that the tree was a reliable integrator of magmatic CO 2 emissions over most of this area. The tree-ring 14 C record implied that magmatic CO 2 emissions from the HLTK were relatively stable from 1998 to 2009, nearly doubled from 2009 to 2011, and then declined by the 2012 growing season. The initial increase in CO 2 emissions was detected during the growing season that immediately followed the largest (February 2010) peak in the occurrence rate of shallow earthquakes. Migration of CO 2 -rich magmatic fluids may have driven observed patterns of elevated deep, then shallow seismicity, while the relationship between pore fluid pressures within a shallow (upper 3 km of crust) fluid reservoir and permeability structure of the reservoir cap rock may have controlled the temporal pattern of surface CO 2 emissions.

California↗

USGS2014

Multi-scale observations of the variability of magmatic CO2 emissions, Mammoth Mountain, CA, USA

One of the primary indicators of volcanic unrest at Mammoth Mountain is diffuse emission of magmatic CO 2 , which can effectively track this unrest if its variability in space and time and relationship to near-surface meteorological and hydrologic phenomena versus those occurring at depth beneath the mountain are understood. In June–October 2013, we conducted accumulation chamber soil CO 2 flux surveys and made half-hourly CO 2 flux measurements with automated eddy covariance and accumulation chamber (auto-chamber) instrumentation at the largest area of diffuse CO 2 degassing on Mammoth Mountain (Horseshoe Lake tree kill; HLTK). Estimated CO 2 emission rates for HLTK based on 20 June, 30 July, and 24–25 October soil CO 2 flux surveys were 165, 172, and 231 t d − 1 , respectively. The average (June–October) CO 2 emission rate estimated for this area was 123 t d − 1 based on an inversion of 4527 eddy covariance CO 2 flux measurements and corresponding modeled source weight functions. Average daily eddy covariance and auto-chamber CO 2 fluxes consistently declined over the four-month observation time. Wavelet analysis of auto-chamber CO 2 flux and environmental parameter time series was used to evaluate the periodicity of, and local correlation between these variables in time–frequency space. Overall, CO 2 emissions at HLTK were highly dynamic, displaying short-term (hourly to weekly) temporal variability related to meteorological and hydrologic changes, as well as long-term (monthly to multi-year) variations related to migration of CO 2 -rich magmatic fluids beneath the volcano. Accumulation chamber soil CO 2 flux surveys were also conducted in the four additional areas of diffuse CO 2 degassing on Mammoth Mountain in July–August 2013. Summing CO 2 emission rates for all five areas yielded a total for the mountain of 311 t d − 1 , which may suggest that emissions returned to 1998–2009 levels, following an increase from 2009 to 2011.

California↗

USGS2014

Geochemical investigation of the hydrothermal system on Akutan Island, Alaska, July 2012

We have studied the geochemistry of the hot springs on Akutan Island in detail for the first time since the early 1980s. Springs in four discrete groups (A-D) along Hot Springs Creek showed generally higher temperatures and substantially higher Na, Ca, and Cl concentrations than previously reported, and total hot-spring discharge has also increased markedly. The springs now account for a heat output of ~29 MW, about an order of magnitude more than in 1981. Gas samples from the hot springs and from a fumarolic area on the flank of Akutan Volcano show high 3 He/ 4 He ratios (>6.4 RA) after correction for air contamination and reveal a common magmatic heat source. Hot-spring gases are unusually rich in N 2 , Ar, and CH 4 , suggesting that the water has boiled and lost CO 2 during upflow beneath the flank fumarole field. Gas geothermometry calculations applied to the flank fumarole field implies temperatures of 200–240 °C for the reservoir, and Na-K-Ca geothermometry implies temperatures near 180 °C for the outflow waters that feed the hot springs. The results of our study confirm the existence of a substantial geothermal resource on the island.

Alaska↗