Christopher D. Farrar

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.

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.

Eddy covariance measurement of CO2 flux to the atmosphere from an area of high volcanogenic emissions, Mammoth Mountain, California

Abstract Three pilot studies were performed to assess application of the eddy covariance micrometeorological method in the measurement of carbon dioxide (CO 2 ) flux of volcanic origin. The selected study area is one of high diffuse CO 2 emission on Mammoth Mountain, CA. Because terrain and source characteristics make this a complex setting for this type of measurement, added consideration was given to source area and upwind fetch. Footprint analysis suggests that the eddy covariance measurements were representative of an upwind elliptical source area (3.8×10 3 m 2 ) which can vary with mean wind direction, surface roughness, and atmospheric stability. CO 2 flux averaged 8–16 mg m −2 s −1 (0.7–1.4 kg m −2 day −1 ). Eddy covariance measurements of flux were compared with surface chamber measurements made in separate studies [Geophys. Res. Lett. 25 (1998a) 1947; EOS Trans. 79 (1998) F941.] and were found to be similar.

Deformation near the Casa Diablo geothermal well field and related processes Long Valley caldera, Eastern California, 1993–2000

Abstract Regional first-order leveling lines, which extend from Lee Vining, CA, to Tom’s Place, CA, have been surveyed periodically since 1957 by the U.S. Geological Survey (USGS), the National Geodetic Survey (NGS), and Caltrans. Two of the regional survey lines, or leveling networks, intersect at the Casa Diablo geothermal well field. These leveling networks, referenced to a distant bench mark (C916) near Lee Vining, provide time-series vertical control data of land-surface deformation that began around 1980. These data are also useful for delineating localized subsidence at Casa Diablo related to reservoir pressure and temperature changes owing to geothermal development that began in 1985. A comparison of differences in bench-mark elevations for five time periods between 1983 and 1997 shows the development and expansion of a subsidence bowl at Casa Diablo. The subsidence coincides spatially with the geothermal well field and temporally with the increased production rates and the deepening of injection wells in 1991, which resulted in an increase in the rate of pressure decline. The subsidence, superimposed on a broad area of uplift, totaled about 310 mm by 1997. The USGS established orthogonal tilt arrays in 1983 to better monitor deformation across the caldera. One tilt array (DBR) was established near what would later become the Casa Diablo geothermal well field. This array responded to magmatic intrusions prior to geothermal development, tilting away from the well field. With the start of geothermal fluid extraction in 1985, tilt at the DBR array reversed direction and began tilting into the well field. In 1991, geothermal power production was increased by a factor of four, and reservoir pressures began a period of steep decline. These changes caused a temporary three-fold increase in the tilt rate. The tilt rate became stable in 1993 and was about 40% lower than that measured in 1991–1992, but still greater than the rates measured during 1985–1990. Data from the local leveling networks spanning the well field and the bounding graben were analyzed for several 2-year periods (1993–1995, 1995–1997, and 1997–1999). Annual rates of change across the normal faults bounding the graben have steadily decreased for each 2-year period between 1993 and 1999, reflecting the slowing decline in geothermal reservoir pressure. Horizontal control data from a two-color electronic distance meter (EDM) defined the lateral extent of subsidence at Casa Diablo. The EDM and leveling data elucidate the localized effect of the shallow source of subsidence and the broader effect of the deeper magmatic inflation source. Data from bench marks common to both the vertical and the horizontal control networks were used to assess the effect of subsidence on the EDM base station (CASA). Modeling of geodetic data collected during periods of little or no magmatic inflation indicated that the CASA two-color EDM station is being drawn toward the well field at a rate of 3–5 mm/yr.

Inferences on the hydrothermal system beneath the resurgent dome in Long Valley Caldera, east-central California, USA, from recent pumping tests and geochemical sampling

Abstract Quaternary volcanic unrest has provided heat for episodic hydrothermal circulation in the Long Valley caldera, including the present-day hydrothermal system, which has been active over the past 40 kyr. The most recent period of crustal unrest in this region of east-central California began around 1980 and has included periods of intense seismicity and ground deformation. Uplift totaling more than 0.7 m has been centered on the caldera’s resurgent dome, and is best modeled by a near-vertical ellipsoidal source centered at depths of 6–7 km. Modeling of both deformation and microgravity data now suggests that (1) there are two inflation sources beneath the caldera, a shallower source 7–10 km beneath the resurgent dome and a deeper source ∼15 km beneath the caldera’s south moat and (2) the shallower source may contain components of magmatic brine and gas. The Long Valley Exploration Well (LVEW), completed in 1998 on the resurgent dome, penetrates to a depth of 3 km directly above this shallower source, but bottoms in a zone of 100°C fluid with zero vertical thermal gradient. Although these results preclude extrapolations of temperatures at depths below 3 km, other information obtained from flow tests and fluid sampling at this well indicates the presence of magmatic volatiles and fault-related permeability within the metamorphic basement rocks underlying the volcanic fill. In this paper, we present recently acquired data from LVEW and compare them with information from other drill holes and thermal springs in Long Valley to delineate the likely flow paths and fluid system properties under the resurgent dome. Additional information from mineralogical assemblages in core obtained from fracture zones in LVEW documents a previous period of more vigorous and energetic fluid circulation beneath the resurgent dome. Although this system apparently died off as a result of mineral deposition and cooling (and/or deepening) of magmatic heat sources, flow testing and tidal analyses of LVEW water level data show that relatively high permeability and strain sensitivity still exist in the steeply dipping principal fracture zone penetrated at a depth of 2.6 km. The hydraulic properties of this zone would allow a pressure change induced at distances of several kilometers below the well to be observable within a matter of days. This indicates that continuous fluid pressure monitoring in the well could provide direct evidence of future intrusions of magma or high-temperature fluids at depths of 5–7 km.