Geology Reports⌕ Search

Geology topics

S. A. Leake

Publications and source records attributed to S. A. Leake.

At least 19 recordsLinked to original sources
USGS2013

Review: groundwater in Alaska (USA)

Groundwater in the US state of Alaska is critical to both humans and ecosystems. Interactions among physiography, ecology, geology, and current and past climate have largely determined the location and properties of aquifers as well as the timing and magnitude of fluxes to, from, and within the groundwater system. The climate ranges from maritime in the southern portion of the state to continental in the Interior, and arctic on the North Slope. During the Quaternary period, topography and rock type have combined with glacial and periglacial processes to develop the unconsolidated alluvial aquifers of Alaska and have resulted in highly heterogeneous hydrofacies. In addition, the long persistence of frozen ground, whether seasonal or permanent, greatly affects the distribution of aquifer recharge and discharge. Because of high runoff, a high proportion of groundwater use, and highly variable permeability controlled in part by permafrost and seasonally frozen ground, understanding groundwater/surface-water interactions and the effects of climate change is critical for understanding groundwater availability and the movement of natural and anthropogenic contaminants.

Alaska↗
USGS2010

A new capture fraction method to map how pumpage affects surface water flow

All groundwater pumped is balanced by removal of water somewhere, initially from storage in the aquifer and later from capture in the form of increase in recharge and decrease in discharge. Capture that results in a loss of water in streams, rivers, and wetlands now is a concern in many parts of the United States. Hydrologists commonly use analytical and numerical approaches to study temporal variations in sources of water to wells for select points of interest. Much can be learned about coupled surface/groundwater systems, however, by looking at the spatial distribution of theoretical capture for select times of interest. Development of maps of capture requires (1) a reasonably well-constructed transient or steady state model of an aquifer with head-dependent flow boundaries representing surface water features or evapotranspiration and (2) an automated procedure to run the model repeatedly and extract results, each time with a well in a different location. This paper presents new methods for simulating and mapping capture using three-dimensional groundwater flow models and presents examples from Arizona, Oregon, and Michigan. Journal compilation ?? 2010 National Ground Water Association. No claim to original US government works.

Ground Water↗
USGS2007

A new ghost-node method for linking different models and initial investigations of heterogeneity and nonmatching grids

A flexible, robust method for linking parent (regional-scale) and child (local-scale) grids of locally refined models that use different numerical methods is developed based on a new, iterative ghost-node method. Tests are presented for two-dimensional and three-dimensional pumped systems that are homogeneous or that have simple heterogeneity. The parent and child grids are simulated using the block-centered finite-difference MODFLOW and control-volume finite-element FEHM models, respectively. The models are solved iteratively through head-dependent (child model) and specified-flow (parent model) boundary conditions. Boundary conditions for models with nonmatching grids or zones of different hydraulic conductivity are derived and tested against heads and flows from analytical or globally-refined models. Results indicate that for homogeneous two- and three-dimensional models with matched grids (integer number of child cells per parent cell), the new method is nearly as accurate as the coupling of two MODFLOW models using the shared-node method and, surprisingly, errors are slightly lower for nonmatching grids (noninteger number of child cells per parent cell). For heterogeneous three-dimensional systems, this paper compares two methods for each of the two sets of boundary conditions: external heads at head-dependent boundary conditions for the child model are calculated using bilinear interpolation or a Darcy-weighted interpolation; specified-flow boundary conditions for the parent model are calculated using model-grid or hydrogeologic-unit hydraulic conductivities. Results suggest that significantly more accurate heads and flows are produced when both Darcy-weighted interpolation and hydrogeologic-unit hydraulic conductivities are used, while the other methods produce larger errors at the boundary between the regional and local models. The tests suggest that, if posed correctly, the ghost-node method performs well. Additional testing is needed for highly heterogeneous systems. ?? 2007 Elsevier Ltd. All rights reserved.

Advances in Water Resources↗
USGS2007

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↗
USGS2006

Comparison of local grid refinement methods for MODFLOW

Many ground water modeling efforts use a finite-difference method to solve the ground water flow equation, and many of these models require a relatively fine-grid discretization to accurately represent the selected process in limited areas of interest. Use of a fine grid over the entire domain can be computationally prohibitive; using a variably spaced grid can lead to cells with a large aspect ratio and refinement in areas where detail is not needed. One solution is to use local-grid refinement (LGR) whereby the grid is only refined in the area of interest. This work reviews some LGR methods and identifies advantages and drawbacks in test cases using MODFLOW-2000. The first test case is two dimensional and heterogeneous; the second is three dimensional and includes interaction with a meandering river. Results include simulations using a uniform fine grid, a variably spaced grid, a traditional method of LGR without feedback, and a new shared node method with feedback. Discrepancies from the solution obtained with the uniform fine grid are investigated. For the models tested, the traditional one-way coupled approaches produced discrepancies in head up to 6.8% and discrepancies in cell-to-cell fluxes up to 7.1%, while the new method has head and cell-to-cell flux discrepancies of 0.089% and 0.14%, respectively. Additional results highlight the accuracy, flexibility, and CPU time trade-off of these methods and demonstrate how the new method can be successfully implemented to model surface water-ground water interactions. Copyright ?? 2006 The Author(s).

Conference Paper↗
USGS2006

Ground-water surface-water interactions and long-term change in riverine riparian vegetation in the southwestern United States

Riverine riparian vegetation has changed throughout the southwestern United States, prompting concern about losses of habitat and biodiversity. Woody riparian vegetation grows in a variety of geomorphic settings ranging from bedrock-lined channels to perennial streams crossing deep alluvium and is dependent on interaction between ground-water and surface-water resources. Historically, few reaches in Arizona, southern Utah, or eastern California below 1530 m elevation had closed gallery forests of cottonwood and willow; instead, many alluvial reaches that now support riparian gallery forests once had marshy grasslands and most bedrock canyons were essentially barren. Repeat photography using more than 3000 historical images of rivers indicates that riparian vegetation has increased over much of the region. These increases appear to be related to several factors, notably the reduction in beaver populations by trappers in the 19th century, downcutting of arroyos that drained alluvial aquifers between 1880 and 1910, the frequent recurrence of winter floods during discrete periods of the 20th century, an increased growing season, and stable ground-water levels. Reductions in riparian vegetation result from agricultural clearing, excessive ground-water use, complete flow diversion, and impoundment of reservoirs. Elimination of riparian vegetation occurs either where high ground-water use lowers the water table below the rooting depth of riparian species, where base flow is completely diverted, or both. We illustrate regional changes using case histories of the San Pedro and Santa Cruz Rivers, which are adjacent watersheds in southern Arizona with long histories of water development and different trajectories of change in riparian vegetation.

Conference Paper↗
USGS2006

Hydrogeologic framework refinement, ground-water flow and storage, water-chemistry analyses, and water-budget components of the Yuma area, southwestern Arizona and southeastern California

The ground-water and surface-water system in the Yuma area in southwestern Arizona and southeastern California is managed intensely to meet water-delivery requirements of customers in the United States, to manage high ground-water levels in the valleys, and to maintain treaty-mandated water-quality and quantity requirements of Mexico. The following components in this report, which were identified to be useful in the development of a ground-water management model, are: (1) refinement of the hydrogeologic framework; (2) updated water-level maps, general ground-water flow patterns, and an estimate of the amount of ground water stored in the mound under Yuma Mesa; (3) review and documentation of the ground-water budget calculated by the Bureau of Reclamation, U.S. Department of the Interior (Reclamation); and (4) water-chemistry characterization to identify the spatial distribution of water quality, information on sources and ages of ground water, and information about the productive-interval depths of the aquifer. A refined three-dimensional digital hydrogeologic framework model includes the following hydrogeologic units from bottom to top: (1) the effective hydrologic basement of the basin aquifer, which includes the Pliocene Bouse Formation, Tertiary volcanic and sedimentary rocks, and pre-Tertiary metamorphic and plutonic rocks; (2) undifferentiated lower units to represent the Pliocene transition zone and wedge zone; (3) coarse-gravel unit; (4) lower, middle, and upper basin fill to represent the upper, fine-grained zone between the top of the coarse-gravel unit and the land surface; and (5) clay A and clay B. Data for the refined model includes digital elevation models, borehole lithology data, geophysical data, and structural data to represent the geometry of the hydrogeologic units. The top surface of the coarse-gravel unit, defined by using borehole and geophysical data, varies similarly to terraces resulting from the down cutting of the Colorado River. Clay A is nearly the same as the previous conceptual hydrogeologic model definition (Olmsted and others, 1973), except for a minor westward extension from the city of Yuma. Clay B is extended to the southerly international boundary and increased in areal extent by about two-thirds of the original extent (Olmsted and others, 1973). The other hydrogeologic units generally are the same as in the previous conceptual hydrogeologic model. Before development, the Colorado and Gila Rivers were the sources of nearly all the ground water in the Yuma area through direct infiltration of water from river channels and annual overbank flooding. After construction of upstream reservoirs and clearing and irrigation of the floodplains, the rivers now act as drains for the ground water. Ground-water levels in most of the Yuma area are higher now than they were in predevelopment time. A general gradient of ground-water flow toward the natural discharge area south of the Yuma area still exists, but many other changes in flow are evident. Ground water in Yuma Valley once flowed away from the Colorado River, but now has a component of flow towards the river and Mexicali Valley. A ground-water mound has formed under Yuma Mesa from long-term surface-water irrigation; about 600,000 to 800,000 acre-ft of water are stored in the mound. Ground-water withdrawals adjacent to the southerly international boundary have resulted in water-level declines in that area. The reviewed and documented water budget includes the following components: (1) recharge in irrigated areas, (2) evapotranspiration by irrigated crops and phreatophytes, (3) ground-water return flow to the Colorado River, and (4) ground-water withdrawals (including those in Mexicali Valley). Recharge components were calculated by subtracting the amount of water used by crops from the amount of water delivered. Evapotranspiration rates were calculated on the basis of established methods, thus were appropriate for input to the ground-water flow model developed by the Bureau of Reclamation (William Greer, hydrologist, Bureau of Reclamation, written commun., 2005). Evapotranspiration by crops and phreatophytes were calculated by using crop coefficient methods and meteorological data. Other methods of calculating evapotranspiration rates by using combinations of satellite imagery and ground-based data could be used for higher spatial and temporal resolution. Ground-water return flow during years of low flow on the Colorado River (1972–82, 1987–92, and 1994–96) averaged 79,000 acre-ft per year. Ground-water withdrawal data for 1970–99 were similar to other estimates made by the U.S. Geological Survey for the Yuma area. New water-chemistry data were collected in 12 wells and 8 canals/drains to characterize spatial patterns in chemical constituents, determine isotopic ages of water, infer possible sources of ground water, and locate the vertical intervals of the aquifer that contribute most water to wells. Depth-dependent samples were collected at one of the wells (YM-10). A large quantity of water-quality data were compiled from Bureau of Reclamation and U.S. Geological Survey records and merged into the U.S. Geological Survey National Water Information System database. New samples were analyzed for major ions, nutrients, stable isotopes of oxygen and hydrogen, tritium ( 3 H), and carbon-14 ( 14 C) (along with C 13 /C 12 ratios). Light values of oxygen-18 ( 18 O) and deuterium ( 2 H, D) in well 242-2 indicate recharge from the Colorado River. Heavy water samples from wells 242-22, CADC, and Mesa del Sol indicate local recharge sources. Tritium data indicate there is young water in wells in the valleys and near the edge of Yuma Mesa, while older water is found far from the Colorado River. 14 C data indicate that water from wells near the southerly international boundary is at least several thousand years old.

Arizona, California↗
USGS2006

User guide for the farm process (FMP1) for the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model, MODFLOW-2000

There is a need to estimate dynamically integrated supply-and-demand components of irrigated agriculture as part of the simulation of surface-water and ground-water flow. To meet this need, a computer program called the Farm Process (FMP1) was developed for the U.S. Geological Survey three-dimensional finite-difference modular ground-water flow model, MODFLOW- 2000 (MF2K). The FMP1 allows MF2K users to simulate conjunctive use of surface- and ground water for irrigated agriculture for historical and future simulations, water-rights issues and operational decisions, nondrought and drought scenarios. By dynamically integrating farm delivery requirement, surface- and ground-water delivery, as well as irrigation-return flow, the FMP1 allows for the estimation of supplemental well pumpage. While farm delivery requirement and irrigation return flow are simulated by the FMP1, the surface-water delivery to the farm can be simulated optionally by coupling the FMP1 with the Streamflow Routing Package (SFR1) and the farm well pumping can be simulated optionally by coupling the FMP1 to the Multi-Node Well (MNW) Package. In addition, semi-routed deliveries can be specified that are associated with points of diversion in the SFR1 stream network. Nonrouted surface-water deliveries can be specified independently of any stream network. The FMP1 maintains a dual mass balance of a farm budget and as part of the ground-water budget. Irrigation demand, supply, and return flow are in part subject to head-dependent sources and sinks such as evapotranspiration from ground water and leakage between the conveyance system and the aquifer. Farm well discharge and farm net recharge are source/sink terms in the FMP1, which depend on transpiration uptake from ground water and other head dependent consumptive use components. For heads rising above the bottom of the root zone, the actual transpiration is taken to vary proportionally with the depth of the active root zone, which can be restricted by anoxia or wilting. Depths corresponding to anoxia- or wilting-related pressure heads within the root zone are found using analytical solutions of a vertical pseudo steady-state pressure- head distribution over the depth of the total root zone (Consumptive Use Concept 1). Alternatively, a simpler, conceptual model is available, which defines how consumptive use (CU) components vary with changing head (CU Concept 2). Subtracting the ground water and precipitation transpiration components from the total transpiration yields a transpiratory irrigation requirement for each cell. The total farm delivery requirement (TFDR) then is determined as cumulative transpiratory and evaporative irrigation requirements of all farm cells and increased sufficiently to compensate for inefficient use from irrigation with respect to plant consumption. The TFDR subsequently is satisfied with surface- and ground-water delivery, respectively constrained by allotments, water rights, or maximum capacities. Five economic and noneconomic drought response policies can be applied optionally, if the potential supply of surface water and ground water is insufficient to meet the crop demand: acreage-optimization with or without a water conservation pool, deficit irrigation with or without water-stacking, and zero policy.

Techniques and Methods↗
USGS2006

New ghost-node method for linking different models with varied grid refinement

A flexible, robust method for linking grids of locally refined ground-water flow models constructed with different numerical methods is needed to address a variety of hydrologic problems. This work outlines and tests a new ghost-node model-linking method for a refined "child" model that is contained within a larger and coarser "parent" model that is based on the iterative method of Steffen W. Mehl and Mary C. Hill (2002, Advances in Water Res., 25, p. 497-511; 2004, Advances in Water Res., 27, p. 899-912). The method is applicable to steady-state solutions for ground-water flow. Tests are presented for a homogeneous two-dimensional system that has matching grids (parent cells border an integer number of child cells) or nonmatching grids. The coupled grids are simulated by using the finite-difference and finite-element models MODFLOW and FEHM, respectively. The simulations require no alteration of the MODFLOW or FEHM models and are executed using a batch file on Windows operating systems. Results indicate that when the grids are matched spatially so that nodes and child-cell boundaries are aligned, the new coupling technique has error nearly equal to that when coupling two MODFLOW models. When the grids are nonmatching, model accuracy is slightly increased compared to that for matching-grid cases. Overall, results indicate that the ghost-node technique is a viable means to couple distinct models because the overall head and flow errors relative to the analytical solution are less than if only the regional coarse-grid model was used to simulate flow in the child model's domain.

Conference Paper↗
USGS2005

Assignment of boundary conditions in embedded ground water flow models

Many small-scale ground water models are too small to incorporate distant aquifer boundaries. If a larger-scale model exists for the area of interest, flow and head values can be specified for boundaries in the smaller-scale model using values from the larger-scale model. Flow components along rows and columns of a large-scale block-centered finite-difference model can be interpolated to compute horizontal flow across any segment of a perimeter of a small-scale model. Head at cell centers of the larger-scale model can be interpolated to compute head at points on a model perimeter. Simple linear interpolation is proposed for horizontal interpolation of horizontal-flow components. Bilinear interpolation is proposed for horizontal interpolation of head values. The methods of interpolation provided satisfactory boundary conditions in tests using models of hypothetical aquifers.

Groundwater↗
USGS2005

Simulating reservoir leakage in ground-water models

Leakage to ground water resulting from the expansion and contraction of reservoirs cannot be easily simulated by most ground-water flow models. An algorithm, entitled the Reservoir Package, was developed for the United States Geological Survey (USGS) three-dimensional finite-difference modular ground-water flow model MODFLOW. The Reservoir Package automates the process of specifying head-dependent boundary cells, eliminating the need to divide a simulation into many stress periods while improving accuracy in simulating changes in ground-water levels resulting from transient reservoir stage. Leakage between the reservoir and the underlying aquifer is simulated for each model cell corrresponding to the inundated area by multiplying the head difference between the reservoir and the aquifer with the hydraulic conductance of the reservoir-bed sediments.

Groundwater↗
USGS2003

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↗
USGS1991

Documentation of a computer program to simulate aquifer-system compaction using the modular finite-difference ground-water flow model

Removal of ground water by pumping from aquifers may result in compaction of compressible fine-grained beds that are within or adjacent to the aquifers. Compaction of the sediments and resulting land subsidence may be permanent if the head declines result in vertical stresses beyond the previous maximum stress. The process of permanent compaction is not routinely included in simulations of ground-water flow. To simulate storage changes from both elastic and inelastic compaction, a computer program was written for use with the U.S. Geological Survey modular finite-difference ground- water flow model. The new program, the Interbed-Storage Package, is designed to be incorporated into this model. In the Interbed-Storage Package, elastic compaction or expansion is assumed to be proportional to change in head. The constant of proportionality is the product of the skeletal component of elastic specific storage and the thickness of the sediments. Similarly, inelastic compaction is assumed to be proportional to decline in head. The constant of proportionality is the product of the skeletal component of inelastic specific storage and the thickness of the sediments. Storage changes are incorporated into the ground-water flow model by adding an additional term to the right-hand side of the flow equation. Within a model time step, the package appropriately apportions storage changes between elastic and inelastic components on the basis of the relation of simulated head to the previous minimum (preconsolidation) head. Two tests were performed to verify that the package works correctly. The first test compared model-calculated storage and compaction changes to hand-calculated values for a three-dimensional simulation. Model and hand-calculated values were essentially equal. The second test was performed to compare the results of the Interbed-Storage Package with results of the one-dimensional Helm compaction model. This test problem simulated compaction in doubly draining confining beds stressed by head changes in adjacent aquifers. The Interbed-Storage Package and the Helm model computed essentially equal values of compaction. Documentation of the Interbed-Storage Package includes data input instructions, flow charts, narratives, and listings for each of the five modules included in the package. The documentation also includes an appendix describing input instructions and a listing of a computer program for time-variant specified-head boundaries. That package was developed to reduce the amount of data input and output associated with one of the Interbed-Storage Package test problems.

Techniques of Water-Resources Investigations↗
USGS1991

Simulation of vertical compaction in models of regional ground-water flow

A new computer program was developed to simulate vertical compaction in models of regional ground-water flow. The program accounts for ground-water storage changes and compaction in discontinuous interbeds or in extensive confining beds. The new program is a package for 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 load can be treated as a function of water-table elevation, and compaction is a function of computed changes in effective stress at the center of a model layer. Thickness of compressible sediments in an unconfined model layer can vary in proportion to saturated thickness. The new package was tested by comparison with an existing model of one-dimensional compaction.

Conference Paper↗
USGS1988

Documentation of a computer program to simulate aquifer-system compaction using the modular finite-difference ground-water flow model

The process of permanent compaction is not routinely included in simulations of groundwater flow. To simulate storage changes from both elastic and inelastic compaction, a computer program was written for use with the U. S. Geological Survey modular finite-difference groundwater flow model. The new program is called the Interbed-Storage Package. In the Interbed-Storage Package, elastic compaction or expansion is assumed to be proportional to change in head. The constant of proportionality is the product of skeletal component of elastic specific storage and thickness of the sediments. Similarly, inelastic compaction is assumed to be proportional to decline in head. The constant of proportionality is the product of the skeletal component of inelastic specific storage and the thickness of the sediments. Storage changes are incorporated into the groundwater flow model by adding an additional term to the flow equation. Within a model time step, the package appropriately apportions storage changes between elastic and inelastic components on the basis of the relation of simulated head to the previous minimum head. Another package that allows for a time-varying specified-head boundary is also documented. This package was written to reduce the data requirements for test simulations of the Interbed-Storage Package. (USGS)

Open-File Report↗