Michael L. Sorey

Carbon dioxide and helium emissions from a reservoir of magmatic gas beneath Mammoth Mountain, California

Carbon dioxide and helium with isotopic compositions indicative of a magmatic source ( δ13C = −4.5 to −5‰, 3He/ 4He = 4.5 to 6.7 RA) are discharging at anomalous rates from Mammoth Mountain, on the southwestern rim of the Long Valley caldera in eastern California. The gas is released mainly as diffuse emissions from normal-temperature soils, but some gas issues from steam vents or leaves the mountain dissolved in cold groundwater. The rate of gas discharge increased significantly in 1989 following a 6-month period of persistent earthquake swarms and associated strain and ground deformation that has been attributed to dike emplacement beneath the mountain. An increase in the magmatic component of helium discharging in a steam vent on the north side of Mammoth Mountain, which also began in 1989, has persisted until the present time. Anomalous CO2 discharge from soils first occurred during the winter of 1990 and was followed by observations of several areas of tree kill and/or heavier than normal needlecast the following summer. Subsequent measurements have confirmed that the tree kills are associated with CO2 concentrations of 30–90% in soil gas and gas flow rates of up to 31,000 g m−2 d−1 at the soil surface. Each of the tree-kill areas and one area of CO2 discharge above tree line occurs in close proximity to one or more normal faults, which may provide conduits for gas flow from depth. We estimate that the total diffuse CO2 flux from the mountain is approximately 520 t/d, and that 30–50 t/d of CO2 are dissolved in cold groundwater flowing off the flanks of the mountain. Isotopic and chemical analyses of soil and fumarolic gas demonstrate a remarkable homogeneity in composition, suggesting that the CO2 and associated helium and excess nitrogen may be derived from a common gas reservoir whose source is associated with some combination of magmatic degassing and thermal metamorphism of metasedimentary rocks. Furthermore, N2/Ar ratios and nitrogen isotopic values indicate that the Mammoth Mountain gases are derived from sources separate from those that supply gas to the hydrothermal system within the Long Valley caldera. Various data suggest that the Mammoth Mountain gas reservoir is a large, low-temperature cap over an isolated hydrothermal system, that it predates the 1989 intrusion, and that it could remain a source of gas discharge for some time.

Dynamics of carbon dioxide emission at Mammoth Mountain, California

Mammoth Mountain, a dormant volcano in the eastern Sierra Nevada, California, has been passively degassing large quantities of cold magmatic CO2 since 1990 following a 6-month-long earthquake swarm associated with a shallow magmatic intrusion in 1989. A search for any link between gas discharge and volcanic hazard at this popular recreation area led us to initiate a detailed study of the degassing process in 1997. Our continuous monitoring results elucidate some of the physical controls that influence dynamics in flank CO2 degassing at this volcano. High coherence between variations in CO2 efflux and variations in atmospheric pressure and wind speed imply that meteorological parameters account for much, if not all of the variability in CO2 efflux rates. Our results help explain differences among previously published estimates of CO2 efflux at Mammoth Mountain and indicate that the long-term (annual) CO2 degassing rate has in fact remained constant since ∼1997. Discounting the possibility of large meteorologically driven temporal variations in gas efflux at other volcanoes may result in spurious interpretations of transients that do not reflect actual geologic processes.

High CO2 emissions through porous media: transport mechanisms and implications for flux measurement and fractionation

Abstract Diffuse emissions of CO2 are known to be large around some volcanoes and hydrothermal areas. Accumulation-chamber measurements of CO2 flux are increasingly used to estimate the total magmatic or metamorphic CO2 released from such areas. To assess the performance of accumulation chamber systems at fluxes one to three orders of magnitude higher than normally encountered in soil respiration studies, a test system was constructed in the laboratory where known fluxes could be maintained through dry sand. Steady-state gas concentration profiles and fractionation effects observed in the 30-cm sand column nearly match those predicted by the Stefan-Maxwell equations, indicating that the test system was functioning successfully as a uniform porous medium. Eight groups of investigators tested their accumulation chamber equipment, all configured with continuous infrared gas analyzers (IRGA), in this system. Over a flux range of ∼200–12,000 g m−2 day−1, 90% of their 203 flux measurements were 0–25% lower than the imposed flux with a mean difference of −12.5%. Although this difference would seem to be within the range of acceptability for many geologic investigations, some potential sources for larger errors were discovered. A steady-state pressure gradient of −20 Pa/m was measured in the sand column at a flux of 11,200 g m−2 day−1. The derived permeability (50 darcies) was used in the dusty-gas model (DGM) of transport to quantify various diffusive and viscous flux components. These calculations were used to demonstrate that accumulation chambers, in addition to reducing the underlying diffusive gradient, severely disrupt the steady-state pressure gradient. The resultant diversion of the net gas flow is probably responsible for the systematically low flux measurements. It was also shown that the fractionating effects of a viscous CO2 efflux against a diffusive influx of air will have a major impact on some important geochemical indicators, such as N2/Ar, δ15N–N2, and 4He/22Ne.

Water-level changes induced by local and distant earthquakes at Long Valley caldera, California

Distant as well as local earthquakes have induced groundwater-level changes persisting for days to weeks at Long Valley caldera, California. Four wells open to formations as deep as 300 m have responded to 16 earthquakes, and responses to two earthquakes in the 3-km-deep Long Valley Exploratory Well (LVEW) show that these changes are not limited to weathered or unconsolidated near-surface rocks. All five wells exhibit water-level variations in response to earth tides, indicating they can be used as low-resolution strainmeters. Earthquakes induce gradual water-level changes that increase in amplitude for as long as 30 days, then return more slowly to pre-earthquake levels. The gradual water-level changes are always drops at wells LKT, LVEW, and CH-10B, and always rises at well CW-3. At a dilatometer just outside the caldera, earthquake-induced strain responses consist of either a step followed by a contractional strain-rate increase, or a transient contractional signal that reaches a maximum in about seven days and then returns toward the pre-earthquake value. The sizes of the gradual water-level changes generally increase with earthquake magnitude and decrease with hypocentral distance. Local earthquakes in Long Valley produce coseismic water-level steps; otherwise the responses to local earthquakes and distant earthquakes are indistinguishable. In particular, water-level and strain changes in Long Valley following the 1992 M7.3 Landers earthquake, 450 km distant, closely resemble those initiated by a M4.9 local earthquake on November 22, 1997, during a seismic swarm with features indicative of fluid involvement. At the LKT well, many of the response time histories are identical for 20 days after each earthquake, and can be matched by a theoretical solution giving the pore pressure as a function of time due to diffusion of a nearby, instantaneous, pressure drop. Such pressure drops could be produced by accelerated inflation of the resurgent dome by amounts too small to be detected by the two-color electronic distance-measuring network. Opening-mode displacement in the south moat, inferred to have followed a M4.9 earthquake on November 22, 1997, could also create extensional strain on the dome and lead to water-level changes similar to those following dome inflation. Contractional strain that could account for earthquake-induced water-level rises at the CW-3 well is inconsistent with geodetic observations. We instead attribute these water-level rises to diffusion of elevated fluid pressure localized in the south moat thermal aquifer. For hydraulic diffusivities appropriate to the upper few hundred meters at Long Valley, an influx of material at temperatures of 300°C can thermally generate pressure of 6 m of water or more, an order of magnitude larger than needed to account for the CW-3 water-level rises. If magma or hot aqueous fluid rises to within 1 km of the surface in the eastern part of the south moat, then hydraulic diffusivities are high enough to allow fluid pressure to propagate to CW-3 on the time scale observed. The data indicate that seismic waves from large distant earthquakes can stimulate upward movement of fluid in the hydrothermal system at Long Valley.

Tracing and quantifying magmatic carbon discharge in cold groundwaters: Lessons learned from Mammoth Mountain, USA

Abstract A major campaign to quantify the magmatic carbon discharge in cold groundwaters around Mammoth Mountain volcano in eastern California was carried out from 1996 to 1999. The total water flow from all sampled cold springs was ≥1.8×10 7 m 3 /yr draining an area that receives an estimated 2.5×10 7 m 3 /yr of recharge, suggesting that sample coverage of the groundwater system was essentially complete. Some of the waters contain magmatic helium with 3 He/ 4 He ratios as high as 4.5 times the atmospheric ratio, and a magmatic component in the dissolved inorganic carbon (DIC) can be identified in virtually every feature sampled. Many waters have a 14 C of 0–5 pmC, a δ 13 C near −5‰, and contain high concentrations (20–50 mmol/l) of CO 2(aq) ; but are otherwise dilute (specific conductance=100–300 μS/cm) with low pH values between 5 and 6. Such waters have previously escaped notice at Mammoth Mountain, and possibly at many other volcanoes, because CO 2 is rapidly lost to the air as the water flows away from the springs, leaving neutral pH waters containing only 1–3 mmol/l HCO 3 − . The total discharge of magmatic carbon in the cold groundwater system at Mammoth Mountain is ∼20 000 t/yr (as CO 2 ), ranging seasonally from about 30 to 90 t/day. Several types of evidence show that this high discharge of magmatic DIC arose in part because of shallow dike intrusion in 1989, but also demonstrate that a long-term discharge possibly half this magnitude (∼10 000 t/yr) predated that intrusion. To sustain a 10 000 t/yr DIC discharge would require a magma intrusion rate of 0.057 km 3 per century, assuming complete degassing of magma with 0.65 wt% CO 2 and a density of 2.7 t/m 3 . The geochemical data also identify a small (

New evidence on the hydrothermal system in Long Valley caldera, California, from wells, fluid sampling, electrical geophysics, and age determinations of hot-spring deposits

Abstract Data collected since 1985 from test drilling, fluid sampling, and geologic and geophysical investigations provide a clearer definition of the hydrothermal system in Long Valley caldera than was previously available. This information confirms the existence of high-temperature (> 200°C) reservoirs within the volcanic fill in parts of the west moat. These reservoirs contain fluids which are chemically similar to thermal fluids encountered in the central and eastern parts of the caldera. The roots of the present-day hydrothermal system (the source reservoir, principal zones of upflow, and the magmatic heat source) most likely occur within metamorphic basement rocks beneath the western part of the caldera. Geothermometer-temperature estimates for the source reservoir range from 214 to 248°C. Zones of upflow of hot water could exist beneath the plateau of moat rhyolite located west of the resurgent dome or beneath Mammoth Mountain. Lateral flow of thermal water away from such upflow zones through reservoirs in the Bishop Tuff and early rhyolite accounts for temperature reversals encountered in most existing wells. Dating of hot-spring deposits from active and inactive thermal areas confirms previous interpretations of the evolution of hydrothermal activity that suggest two periods of extensive hot-spring discharge, one peaking about 300 ka and another extending from about 40 ka to the present. The onset of hydrothermal activity around 40 ka coincides with the initiation of rhyolitic volcanism along the Mono-Inyo Craters volcanic chain that extends beneath the calderas west moat.

Response plan for volcano hazards in the Long Valley Caldera and Mono Craters region, California

Persistent unrest in Long Valley Caldera-characterized by recurring earthquake swarms, inflation of the resurgent dome in the central sections of the caldera, and emissions of magmatic carbon dioxide around Mammoth Mountain-during the last two decades and continuing into the 21st century emphasize that this geologically youthful volcanic system is capable of further volcanic activity. This document describes the U.S. Geological Surveys (USGS) response plan for future episodes of unrest that might augur the onset of renewed volcanism in the caldera or along the Inyo-Mono Craters chain to the north. Central to this response plan is a four-level color code with successive conditions, GREEN (no immediate risk) through RED (eruption under way), reflecting progressively more intense activity levels as summarized in table 1 and 2 and described in detail in section II.

Radiocarbon studies of plant leaves and tree rings from Mammoth Mountain, CA: a long-term record of magmatic CO2 release

Evaluation of 14C in tree rings provides a measure of the flux of magmatic CO2 from Mammoth Mountain both before and after 1994 when copious diffuse emissions were first discovered and linked to tree kill. We analyzed the annual rings of trees with two main purposes: (1) to track changes in the magnitude of magmatic CO2 emission over time, and (2) to determine the onset of magmatic CO2 emission at numerous sites on Mammoth Mountain. The onset of CO2 emission at different areas of tree kill was determined to be in 1990, closely following the seismic events of 1989. At Horseshoe Lake (HSL), CO2 emission was found to have peaked in 1991 and to have subsequently declined by a factor of two through 1998. The tree-ring data also show that emissions of magmatic carbon from cold springs below the tree-kill areas occurred well before 1989. Trees located on the margins of the kill areas or otherwise away from zones of maximum gas discharge were found to be better integrators of magmatic CO2 emission than those located in the center of tree kills. Although quantitative extrapolations from our data to a flux history will require that a relationship be established between 14C depletion in tree rings and average annual magmatic CO2 flux, the pattern of 14C depletion in tree rings is likely to be the most reliable indicator of the long-term changes in the magnitude of CO2 release from Mammoth Mountain.

Chemical and isotopic prediction of aquifer temperatures in the geothermal system at Long Valley, California

Abstract Temperatures of aquifers feeding thermal springs and wells in Long Valley, California, estimated using silica and Na-K-Ca geothermometers and warm spring mixing models, range from 160/dg to about 220°C. This information was used to construct a diagram showing enthalpy-chloride relations for the various thermal waters in the Long Valley region. The enthalpy-chloride information suggests that a 282 ± 10° C aquifer with water containing about 375 mg chloride per kilogram of water is present somewhere deep in the system. That deep water would be related to ∼ 220°C Casa Diablo water by mixing with cold water, and to Hot Creek water by first boiling with steam loss and then mixing with cold water. Oxygen and deuterium isotopic data are consistent with that interpretation. An aquifer at 282°C with 375 mg/kg chloride implies a convective heat flow in Long Valley of 6.6 × 10 7 cal/s.

Episodic thermal perturbations associated with groundwater flow: An example from Kilauea Volcano, Hawaii

present a new temperature-depth profile from a deep well on the summit of Kilauea Volcano, Hawaii, and analyze it in conjunction with a temperature profile measured 26 years earlier. We propose two groundwater flow models to interpret the complex temperature profiles. The first is a modified confined lateral flow model (CLFM) with a continuous flux of hydrothermal fluid. In the second, transient flow model (TFM), slow conductive cooling follows a brief, advective heating event. We carry out numerical simulations to examine the timescales associated with each of the models. Results for both models are sensitive to the initial conditions, and with realistic initial conditions it takes between 750 and 1000 simulation years for either model to match the measured temperature profiles. With somewhat hotter initial conditions, results are consistent with onset of a hydrothermal plume � 550 years ago, coincident with initiation of caldera subsidence. We show that the TFM is consistent with other data from hydrothermal systems and laboratory experiments and perhaps is more appropriate for this highly dynamic environment. The TFM implies that volcano-hydrothermal systems may be dominated by episodic events and that thermal perturbations may persist for several thousand years after hydrothermal flow has ceased. INDEX TERMS: 1878 Hydrology: Water/ energy interactions; 3210 Mathematical Geophysics: Modeling; 8424 Volcanology: Hydrothermal systems (8135); 3230 Mathematical Geophysics: Numerical solutions; 1829 Hydrology: Groundwater hydrology;