Steven E. Ingebritsen

DETECTION OF AQUIFER SYSTEM COMPACTION AND LAND SUBSIDENCE USING INTERFEROMETRIC SYNTHETIC APERTURE RADAR, ANTELOPE VALLEY, MOJAVE DESERT, CALIFORNIA

Interferometric synthetic aperture radar (InSAR) has great potential to detect and quantify land subsidence caused by aquifer system compaction. InSAR maps with high spatial detail and resolution of range displacement (±10 mm in change of land surface elevation) were developed for a groundwater basin (∼103 km2) in Antelope Valley, California, using radar data collected from the ERS-1 satellite. These data allow comprehensive comparison between recent (1993–1995) subsidence patterns and those detected historically (1926–1992) by more traditional methods. The changed subsidence patterns are generally compatible with recent shifts in land and water use. The InSAR-detected patterns are generally consistent with predictions based on a coupled model of groundwater flow and aquifer system compaction. The minor inconsistencies may reflect our imperfect knowledge of the distribution and properties of compressible sediments. When used in conjunction with coincident measurements of groundwater levels and other geologic information, InSAR data may be useful for constraining parameter estimates in simulations of aquifer system compaction.

Geological implications of a permeability-depth curve for the continental crust

The decrease in permeability ( k ) of the continental crust with depth ( z ), as constrained by geothermal data and calculated fluid flux during metamorphism, is given by log k = −14 − 3.2 log z, where k is in meters squared and z is in kilometers. At moderate to great crustal depths (>∼5 km), this curve is defined mainly by data from prograde metamorphic systems, and is thus applicable to orogenic belts where the crust is being thickened and/or heated; lower permeabilities may occur in stable cratonic regions. This k-z relation implies that typical metamorphic fluid flux values of ∼10−11 m/s are consistent with fluid pressures significantly above hydrostatic values. The k-z curve also predicts that metamorphic CO2 flux from large orogens may be sufficient to cause significant climatic effects, if retrograde carbonation reactions are minimal, and suggests a significant capacity for diffuse degassing of Earth (1015–1016 g/yr) in tectonically active regions.

Multiphase groundwater flow near cooling plutons

We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200°C, thermal pressurization, and temperature-dependent rock properties. A series of experiments examines the effects of host-rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host-rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam-rich systems develop where permeability is of the order of 10−15 m2. Temperatures and life spans of systems decrease with increasing permeability. Conduction-dominated systems, in which permeabilities are ≤10−16 m2, persist longer but exhibit relatively modest increases in near-surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low-permeability caprock effectively decouples the topographically and density-driven systems and stabilizes the mixing interface between them thereby defining a likely ore-forming environment.

Use of precipitation and groundwater isotopes to interpret regional hydrology on a tropical volcanic island: Kilauea volcano area, Hawaii

Isotope tracer methods were used to determine flow paths, recharge areas, and relative age for groundwater in the Kilauea volcano area of the Island of Hawaii. A network of up to 66 precipitation collectors was emplaced in the study area and sampled twice yearly for a 3-year period. Stable isotopes in rainfall show three distinct isotopic gradients with elevation, which are correlated with trade wind, rain shadow, and high- elevation climatological patterns. Temporal variations in precipitation isotopes are controlled more by the frequency of storms than by seasonal temperature fluctuations. Results from this study suggest that (1) sampling network design must take into account areal variations in rainfall patterns on islands and in continental coastal areas and (2) isotope/elevation gradients on other tropical islands may be predictable on the basis of similar climatology. Groundwater was sampled yearly in coastal springs, wells, and a few high-elevation springs. Areal contrasts in groundwater stable isotopes and tritium indicate that the volcanic rift zones compartmentalize the regional groundwater system, isolating the groundwater south of Kilaueas summit and rift zones. Part of the Southwest Rift Zone appears to act as a conduit for water from higher elevation, but there is no evidence for downrift flow in the springs and shallow wells sampled in the lower East Rift Zone.

Heat flow and hydrothermal circulation in the cascade range, north-central Oregon.

In north-central Oregon a large area of near-zero near-surface conductive heat flow occurs in young volcanic rocks of the Cascade Range. Recent advective heat flux measurements and a heat-budget analysis suggest that ground-water circulation sweeps sufficient heat out of areas where rocks younger than 6 Ma (million years ago) are exposed to account for the anomalously high advective and conductive heat discharge measured in older rocks at lower elevations. Earlier workers have proposed that an extensive midcrustal magmatic heat source is responsible for this anomalously high heat flow. Instead, high heat flow in the older rocks may be a relatively shallow phenomenon caused by regional ground-water flow. Any deeper anomaly may be relatively narrow, spatially variable, and essentially confined to the Quaternary (less than 2 Ma) arc. Magmatic intrusion at a rate of 9 to 33 cubic kilometers per kilometer of arc length per million years can account for the total heat flow anomaly. Deep drilling in the areas of high heat flow in the older rocks could indicate which model is more appropriate for the near-surface heat flow data.

Controls on geyser periodicity

Geyser eruption frequency is not constant over time and has been shown to vary with small (≤10–6) strains induced by seismic events, atmospheric loading, and Earth tides. The geyser system is approximated as a permeable conduit of intensely fractured rock surrounded by a less permeable rock matrix. Numerical simulation of this conceptual model yields a set of parameters that controls geyser existence and periodicity. Much of the responsiveness to remote seismicity and other small strains in the Earth can be explained in terms of variations in permeability and lateral recharge rates.

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;

Magmatic intrusion west of Three Sisters, central Oregon, USA: The perspective from spring geochemistry

A geochemical investigation of springs near Three Sisters volcanoes was conducted in response to the detection of crustal uplift west of the peaks. Dilute, low-temperature springs near the center of uplift show 3 He/ 4 He ratios

Implications of crustal permeability for fluid movementbetween terrestrial fluid reservoirs

7RA (RA is the ratio in air), and transport in total ;16 MW of heat and ;180 g/s of magmatic carbon (as CO2). These anomalous conditions clearly reflect the influence of magma, but they seemingly predate the onset of the present uplift and derive from a previous event. Episodes of intrusion may thus be more common in this area than the age of eruptive vents would imply.

Modeling the Formation of Porphyry-Copper Ores

Abstract A classic paper by Rubey [Geol. Soc. Amer. Bull 62 (1951) 1111] examined various hypotheses regarding the origin of sea water and concluded that the most likely hypothesis was volcanic outgassing, a view that was generally accepted by Earth scientists for the next several decades. More recent work suggests that the rate of subduction of water is much larger than the volcanic outgassing rate, lending support to hypotheses that either ocean volume has decreased with time, or that the imbalance is offset by continuous replenishment of water by cometary impacts. These alternatives are required in the absence of additional mechanisms for the return of water from subducting lithosphere to the Earths surface. Our recent work on crustal permeability suggests a large capacity for water upflow through tectonically active continental crust, resulting in a heretofore-unrecognized degassing pathway that can accommodate the water-subduction rate. Escape of recycled water via delivery from the mantle through zones of active metamorphism eliminates the mass-balance argument for the loss of ocean volume or extraterrestrial sources.