Search USGSSearch

Geology topics

D. L. Galloway

Publications and source records attributed to D. L. Galloway.

16 recordsLinked to original sources

MODFLOW Ground-Water Model - User Guide to the Subsidence and Aquifer-System Compaction Package (SUB-WT) for Water-Table Aquifers

A new computer program was developed to simulate vertical compaction in models of regional ground-water flow. The program simulates ground-water storage changes and compaction in discontinuous interbeds or in extensive confining units, accounting for stress-dependent changes in storage properties. The new program is a package for MODFLOW, the U.S. Geological Survey modular finite-difference ground-water flow model. Several features of the program make it useful for application in shallow, unconfined flow systems. Geostatic stress can be treated as a function of water-table elevation, and compaction is a function of computed changes in effective stress at the bottom of a model layer. Thickness of compressible sediments in an unconfined model layer can vary in proportion to saturated thickness.

Techniques and Methods

The application of satellite differential SAR interferometry-derived ground displacements in hydrogeology

The application of satellite differential synthetic aperture radar (SAR) interferometry, principally coherent (InSAR) and to a lesser extent, persistent-scatterer (PSI) techniques to hydrogeologic studies has improved capabilities to map, monitor, analyze, and simulate groundwater flow, aquifer-system compaction and land subsidence. A number of investigations over the previous decade show how the spatially detailed images of ground displacements measured with InSAR have advanced hydrogeologic understanding, especially when a time series of images is used in conjunction with histories of changes in water levels and management practices. Important advances include: (1) identifying structural or lithostratigraphic boundaries (e.g. faults or transitional facies) of groundwater flow and deformation; (2) defining the material and hydraulic heterogeneity of deforming aquifer-systems; (3) estimating system properties (e.g. storage coefficients and hydraulic conductivities); and (4) constraining numerical models of groundwater flow, aquifer-system compaction, and land subsidence. As a component of an integrated approach to hydrogeologic monitoring and characterization of unconsolidated alluvial groundwater basins differential SAR interferometry contributes unique information that can facilitate improved management of groundwater resources. Future satellite SAR missions specifically designed for differential interferometry will enhance these contributions. ?? Springer-Verlag 2006.

Hydrogeology Journal

Variability in geyser eruptive timing and its causes: Yellowstone National Park

Field data from Upper Geyser Basin, Yellowstone, indicate that geyser frequency is less sensitive to elastic deformation than might be surmised from a review of the literature. Earth-tide influences are not identifiable in any of the geysers we monitored. Though atmospheric-pressure influences are observed, only long-period variations on the order of 5 mBars or greater seem to influence geyser frequency. Long-distance interconnections between geysers are common and add to the difficulty of identifying strain influences. Additional variations in geyser periodicity may be governed by the internal dynamics of the geysers rather than external influences.

Geophysical Research Letters

Inverse modeling of interbed storage parameters using land subsidence observations, Antelope Valley, California

We use land-subsidence observations from repeatedly surveyed benchmarks and interferometric synthetic aperture radar (InSAR) in Antelope Valley, California, to estimate spatially varying compaction time constants, ??, and inelastic specific skeletal storage coefficients, Skv*, in a previously calibrated regional groundwater flow and subsidence model. The observed subsidence patterns reflect both the spatial distribution of head declines and the spatially variable inelastic skeletal storage coefficient. Using the nonlinear parameter estimation program UCODE we estimate compaction time constants between 3.8 and 285 years. The Skv* values are estimated by linear estimation and range from 0 to almost 0.09. We find that subsidence observations over long time periods are necessary to constrain estimates of the large compaction time constants in Antelope Valley. The InSAR data used in this study cover only a three-year period, limiting their usefulness in constraining these time constants. This problem will be alleviated as more SAR data become available in the future or where time constants are small. By incorporating the resulting parameter estimates in the previously calibrated regional model of groundwater flow and land subsidence we can significantly improve the agreement between simulated and observed land subsidence both in terms of magnitude and spatial extent. The sum of weighted squared subsidence residuals, a common measure of model fit, was reduced by 73% with respect to the original model. However, the ability of the model to adequately reproduce the subsidence observed over only a few years is impaired by the fact that the simulated hydraulic heads over small time periods are often not representative of the actual aquifer hydraulic heads. Errors in the simulated hydraulic aquifer heads constitute the primary limitation of the approach presented here.

Water Resources Research

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

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. ?? 2003 Elsevier B.V. All rights reserved.

California

MODFLOW-2000 ground-water model-user guide to the Subsidence and Aquifer-System Compaction (SUB) Package

This report documents a computer program, the Subsidence and Aquifer-System Compaction (SUB) Package, to simulate aquifer-system compaction and land subsidence using the U.S. Geological Survey modular finite-difference ground-water flow model, MODFLOW-2000. The SUB Package simulates elastic (recoverable) compaction and expansion, and inelastic (permanent) compaction of compressible fine-grained beds (interbeds) within the aquifers. The deformation of the interbeds is caused by head or pore-pressure changes, and thus by changes in effective stress, within the interbeds. If the stress is less than the preconsolidation stress of the sediments, the deformation is elastic; if the stress is greater than the preconsolidation stress, the deformation is inelastic. The propagation of head changes within the interbeds is defined by a transient, one-dimensional (vertical) diffusion equation. This equation accounts for delayed release of water from storage or uptake of water into storage in the interbeds. Properties that control the timing of the storage changes are vertical hydraulic diffusivity and interbed thickness. The SUB Package supersedes the Interbed Storage Package (IBS1) for MODFLOW, which assumes that water is released from or taken into storage with changes in head in the aquifer within a single model time step and, therefore, can be reasonably used to simulate only thin interbeds. The SUB Package relaxes this assumption and can be used to simulate time-dependent drainage and compaction of thick interbeds and confining units. The time-dependent drainage can be turned off, in which case the SUB Package gives results identical to those from IBS1. Three sample problems illustrate the usefulness of the SUB Package. One sample problem verifies that the package works correctly. This sample problem simulates the drainage of a thick interbed in response to a step change in head in the adjacent aquifer and closely matches the analytical solution. A second sample problem illustrates the effects of seasonally varying discharge and recharge to an aquifer system with a thick interbed. A third sample problem simulates a multilayered regional ground-water basin. Model input files for the third sample problem are included in the appendix.

Open-File Report

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 intrussion 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 do not reflect actual geologic processes. ?? 2001 Elsevier Science B.V. All rights reserved.

Earth and Planetary Science Letters

Time-variation of hydrothermal discharge at selected sites in the Western United States: Implications for monitoring

We compiled time series of hydrothermal discharge consisting of 3593 chloride- or heat-flux measurements from 24 sites in the Yellowstone region, the northern Oregon Cascades, Lassen Volcanic National Park and vicinity, and Long Valley, California. At all of these sites the hydrothermal phenomena are believed to be as yet unaffected by human activity, though much of the data collection was driven by mandates to collect environmental-baseline data in acticipation of geothermal development. The time series average 19 years in length and some of the Yellowstone sites have been monitored intermittently for over 30 years. Many sites show strong seasonality but few show clear long-term trends, and at most sites statistically significant decadal-scale trends are absent. Thus, the data provide robust estimates of advective heat flow ranging from ~130 MW in the north-central Oregon Cascades to ~6100 MW in the Yellowstone region, and also document Yellowstone hydrothermal chloride and arsenic fluxes of 1740 and 15-20 g/s, respectively. The discharge time series show little sensitivity to regional tectonic events such as earthquakes or inflation/deflation cycles. Most long-term monitoring to date has focused on high-chloride springs and low-temperature fumaroles. The relative stability of these features suggests that discharge measurements done as part of volcano-monitoring programs should focus instead on high-temperature fumaroles, which may be more immediately linked to the magmatic heat source.

Journal of Volcanology and Geothermal Research

Detection and measurement of land subsidence using Global Positioning System and interferometric synthetic aperture radar, Coachella Valley, California, 1996-98

Land subsidence associated with ground-water-level declines has been recognized as a potential problem in Coachella Valley, California. Since the early 1920s, ground water has been a major source of agricultural, municipal, and domestic supply in the valley, resulting in water-level declines as large as 15 meters (50 feet) through the late 1940s. In 1949, the importation of Colorado River water to the lower Coachella Valley began, resulting in a reduction in ground-water pumping and a recovery of water levels from the 1950s through the 1970s. Since the late 1970s, the demand for water in the valley has exceeded the deliveries of imported surface water, again resulting in increased pumping and ground-water-level declines. The magnitude and temporal occurrence of land subsidence in the lower Coachella Valley are not well known; data are sparse and accuracy varies. Also, the area is tectonically active and has subsided during the past several million years, which further complicates interpretations of the data. Land-surface-elevation data have been collected by many agencies using various methods and different geographic scales; because of this, the -150 millimeters (-0.5 foot) of subsidence determined for the southern parts of the valley for 1930-96 may have a possible error of plus or minus (?)90 millimeters (?0.3 foot). The location, extent, and magnitude of vertical land-surface changes from 1996 to 1998 were determined using Global Positioning System (GPS) and interferometric synthetic aperture radar (InSAR) methods. GPS measurements for 14 monuments in the lower Coachella Valley indicate that the vertical land-surface changes from 1996 to 1998 ranged from -13 to -67 millimeters ? 40 millimeters (-0.04 to -0.22 foot ?0.13 foot). Changes at seven of the monuments exceeded the measurement error of ?40 millimeters (?0.13 foot), which indicates that small amounts of land subsidence occurred at these monuments between 1996 and 1998. Some of the water levels measured in wells near several of these monuments during 1996-98 were the lowest water levels in the recorded histories of the wells. The possible relation between the stresses caused by historically low water levels and the measured vertical changes in land surface suggests that the preconsolidation stress of the aquifer system may have been exceeded during this period and that subsidence may be permanent. Comparisons of several paired monuments and wells indicated that the relation between short-term ground-water-level changes and vertical changes in land surface in the lower Coachella Valley is not clearly defined. Results of InSAR measurements made between 1996 and 1998 indicate that vertical changes in land surface, ranging from about -20 to -70 millimeters ? 5-10 millimeters (-0.07 to -0.23 foot ? 0.02-0.03 foot), occurred in three areas of the Coachella Valley--near Palm Desert, Indian Wells, and Lake Cahuilla. The areas of subsidence near Palm Desert and Indian Wells coincide with areas of substantial ground-water production during 1996-98. The Coachella Valley Water District reported that they had no ground-water production wells in the Lake Cahuilla area but that there may be private production wells in the area. Production from these wells or possibly tectonic activity may be contributing to or causing the subsidence. The geodetic network used for the GPS measurements described in this report covers the area from the Salton Sea on the south to just northwest of Indio. The maps processed using InSAR overlap the part of the geodetic network west of Coachella and north of Lake Cahuilla, and include the Palm Desert area. Both methods of measuring vertical land-surface changes, GPS and InSAR, were used to characterize vertical land-surface changes from the Palm Desert area to the Salton Sea. Because InSAR produces more spatially detailed data over large areas, it generally was useful where vertical land-surface changes were previously unrecognized, such as the

California

Delta subsidence in California: The sinking heart of the state

The Sacramento-San Joaquin River Delta of California once was a great tidal freshwater marsh blanketed by peat and peaty alluvium. Beginning in the late 1800s, levees were built along the stream channels, and the land thus protected from flooding was drained, cleared, and planted. Although the Delta is now an exceptionally rich agricultural area (over a $500 million crop value in 1993), its unique value is as a source of freshwater for the rest of the State. It is the heart of a massive north-to-south waterdelivery system. Much of this water is pumped southward for use in the San Joaquin Valley and elsewhere in central and southern California. The leveed tracts and islands help to protect water-export facilities in the southern Delta from saltwater intrusion by displacing water and maintaining favorable freshwater gradients. However, ongoing subsidence behind the levees reduces levee stability and, thus, threatens to degrade water quality in the massive north-to-south water-transfer system.

California

Use of InSAR to identify land-surface displacements caused by aquifer-system compaction in the Paso Robles area, San Luis Obispo County, California, March to August 1997

The population in San Luis Obispo County has grown steadily during the 1990s, and some land use has been converted from dry farming to grazing to irrigated vineyards and urban areas. Because insufficient surface-water supplies are available to meet the growing demand, ground-water pumpage has increased and the resulting water-level declines have raised concern that this water resource may become overstressed. One particular concern is whether the larger ground-water basins within the county function as large individual basins or whether subsurface structures divide these large basins into smaller subbasins, as differences in ground-water-level data suggest. Interferometric Synthetic Aperture Radar (InSAR) is used in locating land-surface displacement, which may indicate subsurface structures in ground-water basins by determining seasonal and historical land-surface changes. Owing to the high spatial detail of InSAR imagery, the InSAR-derived displacement maps can be used with ground-water-level data to reveal differential aquifer-system compaction related to the presence of geological structures or the distribution of compressible sediments that may define subbasin boundaries. This report describes InSAR displacement maps of the Paso Robles area of San Luis Obispo County and compares them to maps of seasonal changes in ground-water levels to detect the presence of aquifer-system compaction.

Open-File Report

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 (∼10 3 km 2 ) 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.

California

Hydrologic effects associated with the January 17, 1994 Northridge, California, earthquake

This report compiles hydrologic observations in southern California associated with the 1994 Mw = 6.7 Northridge, California earthquake. In southern California, the largest ground water level change was a drop of 52 cm at Crystalaire. Most of the steplike water-level changes recorded following the Northridge earthquake agreed in direction with the sign of the calculated coseismic volume strain field. A spring discharge increase at Southern Pacific Springs was reported to have preceded the earthquake by several days. Outside of southern California, water-level changes were also observed, but are not consistent in sign or size with the static strain field of the earthquake sequence.

California, Nevada, Oregon