Geology Reports⌕ Search

Geology topics

R. T. Hanson

Publications and source records attributed to R. T. Hanson.

At least 19 recordsLinked to original sources
USGS2014

Hydrologic models and analysis of water availability in Cuyama Valley, California

Changes in population, agricultural development practices (including shifts to more water-intensive crops), and climate variability are placing increasingly larger demands on available water resources, particularly groundwater, in the Cuyama Valley, one of the most productive agricultural regions in Santa Barbara County. The goal of this study was to produce a model capable of being accurate at scales relevant to water management decisions that could be considered in the evaluation of the sustainable water supply. The Cuyama Valley Hydrologic Model (CUVHM) was designed to simulate the most important natural and human components of the hydrologic system, including components dependent on variations in climate, thereby providing a reliable assessment of groundwater conditions and processes that can inform water users and help to improve planning for future conditions. Model development included a revision of the conceptual model of the flow system, construction of a precipitation-runoff model using the Basin Characterization Model (BCM), and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model (MF-OWHM). The hydrologic models were calibrated to historical conditions of water and land use and, then, used to assess the use and movement of water throughout the Valley. These tools provide a means to understand the evolution of water use in the Valley, its availability, and the limits of sustainability. The conceptual model identified inflows and outflows that include the movement and use of water in both natural and anthropogenic systems. The groundwater flow system is characterized by a layered geologic sedimentary sequence that—in combination with the effects of groundwater pumping, natural recharge, and the application of irrigation water at the land surface—displays vertical hydraulic-head gradients. Overall, most of the agricultural demand for water in the Cuyama Valley in the initial part of the growing season is supplied by groundwater, which is augmented by precipitation during wet winter and spring seasons. In addition, the amount of groundwater used for irrigation varies from year to year in response to climate variation and can increase dramatically in dry years. Model simulation results, however, also indicated that irrigation may have been less efficient during wet years. Agricultural pumpage is a major component to simulated outflow that is often poorly recorded. Therefore, an integrated, coupled farm-process model is used to estimate historical pumpage for water-balance subregions that evolved with the development of groundwater in the Valley from 1949 through 2010. The integrated hydrologic model includes these water-balance subregions and delineates natural, municipal, and agricultural land use; streamflow networks; and groundwater flow systems. The redefinition of the geohydrologic framework (including the internal architecture of the sedimentary units) and incorporation of these units into the simulation of the regional groundwater flow system indicated that faults have compartmentalized the alluvial deposits into subregions, which have responded differently to regional groundwater flow, locations of recharge, and the effects of development. The Cuyama Valley comprises nine subregions grouped into three regional zones, the Main, Ventucopa Uplands, and Sierra Madre Foothills, which are fault bounded, represent different proportions of the three alluvial aquifers, and have different water quality. The CUVHM uses MF-OWHM to simulate and assess the use and movement of water, including the evolution of land use and related water-balance regions. The model is capable of being accurate at annual to interannual time frames and at subregional to valley-wide spatial scales, which allows for analysis of the groundwater hydrologic budget for the water years 1950–2010, as well as potential assessment of the sustainable use of groundwater. Simulated changes in storage over time showed that significant withdrawals from storage generally occurred not only during drought years (1976–77 and 1988–92) but also during the early stages of industrial agriculture, which was initially dominated by alfalfa production. Since the 1990s, agriculture has shifted to more water-intensive crops. Measured and simulated groundwater levels indicated substantial declines in selected subregions, mining of groundwater that is thousands to tens of thousands of years old, increased groundwater storage depletion, and land subsidence. Most of the recharge occurs in the upland regions of Ventucopa and Sierra Madre Foothills, and the largest fractions of pumpage and storage depletion occur in the Main subregion. The long-term imbalance between inflows and outflows resulted in simulated overdraft (groundwater withdrawals in excess of natural recharge) of the groundwater basin over the 61-year period of 1949–2010. Changes in storage varied considerably from year to year, depending on land use, pumpage, and climate conditions. Climatically driven factors can greatly affect inflows, outflows, and water use by more than a factor of two between wet and dry years. Although precipitation during inter-decadal wet years previously replenished the basin, the water use and storage depletion have lessened the effects of these major recharge events. Simulated and measured water-level altitudes indicated the presence of large areas where depressed water levels have resulted in large desaturated zones in the younger and Older Alluvium layers in the Main-zone subregions. The results of modeled projection of the base-case scenario 61 years into the future indicated that current supply-and-demand are unsustainable and will result in additional groundwater-level declines and related storage depletion and land subsidence. The reduced-supply and reduced-demand projections reduced groundwater storage depletion but may not allow for sustainable agriculture under current demands, agricultural practices, and land use.

California↗
USGS2013

Advective transport observations with MODPATH-OBS--documentation of the MODPATH observation process

The MODPATH-OBS computer program described in this report is designed to calculate simulated equivalents for observations related to advective groundwater transport that can be represented in a quantitative way by using simulated particle-tracking data. The simulated equivalents supported by MODPATH-OBS are (1) distance from a source location at a defined time, or proximity to an observed location; (2) time of travel from an initial location to defined locations, areas, or volumes of the simulated system; (3) concentrations used to simulate groundwater age; and (4) percentages of water derived from contributing source areas. Although particle tracking only simulates the advective component of conservative transport, effects of non-conservative processes such as retardation can be approximated through manipulation of the effective-porosity value used to calculate velocity based on the properties of selected conservative tracers. This program can also account for simple decay or production, but it cannot account for diffusion. Dispersion can be represented through direct simulation of subsurface heterogeneity and the use of many particles. MODPATH-OBS acts as a postprocessor to MODPATH, so that the sequence of model runs generally required is MODFLOW, MODPATH, and MODPATH-OBS. The version of MODFLOW and MODPATH that support the version of MODPATH-OBS presented in this report are MODFLOW-2005 or MODFLOW-LGR, and MODPATH-LGR. MODFLOW-LGR is derived from MODFLOW-2005, MODPATH 5, and MODPATH 6 and supports local grid refinement. MODPATH-LGR is derived from MODPATH 5. It supports the forward and backward tracking of particles through locally refined grids and provides the output needed for MODPATH_OBS. For a single grid and no observations, MODPATH-LGR results are equivalent to MODPATH 5. MODPATH-LGR and MODPATH-OBS simulations can use nearly all of the capabilities of MODFLOW-2005 and MODFLOW-LGR; for example, simulations may be steady-state, transient, or a combination. Though the program name MODPATH-OBS specifically refers to observations, the program also can be used to calculate model prediction of observations. MODPATH-OBS is primarily intended for use with separate programs that conduct sensitivity analysis, data needs assessment, parameter estimation, and uncertainty analysis, such as UCODE_2005, and PEST. In many circumstances, refined grids in selected parts of a model are important to simulated hydraulics, detailed inflows and outflows, or other system characteristics. MODFLOW-LGR and MODPATH-LGR support accurate local grid refinement in which both mass (flows) and energy (head) are conserved across the local grid boundary. MODPATH-OBS is designed to take advantage of these capabilities. For example, particles tracked between a pumping well and a nearby stream, which are simulated poorly if a river and well are located in a single large grid cell, can be simulated with improved accuracy using a locally refined grid in MODFLOW-LGR, MODPATH-LGR, and MODPATH-OBS. The locally-refined-grid approach can provide more accurate simulated equivalents to observed transport between the well and the river. The documentation presented here includes a brief discussion of previous work, description of the methods, and detailed descriptions of the required input files and how the output files are typically used.

Techniques and Methods↗
USGS2012

A method for physically based model analysis of conjunctive use in response to potential climate changes

Potential climate change effects on aspects of conjunctive management of water resources can be evaluated by linking climate models with fully integrated groundwater-surface water models. The objective of this study is to develop a modeling system that links global climate models with regional hydrologic models, using the California Central Valley as a case study. The new method is a supply and demand modeling framework that can be used to simulate and analyze potential climate change and conjunctive use. Supply-constrained and demand-driven linkages in the water system in the Central Valley are represented with the linked climate models, precipitation-runoff models, agricultural and native vegetation water use, and hydrologic flow models to demonstrate the feasibility of this method. Simulated precipitation and temperature were used from the GFDL-A2 climate change scenario through the 21st century to drive a regional water balance mountain hydrologic watershed model (MHWM) for the surrounding watersheds in combination with a regional integrated hydrologic model of the Central Valley (CVHM). Application of this method demonstrates the potential transition from predominantly surface water to groundwater supply for agriculture with secondary effects that may limit this transition of conjunctive use. The particular scenario considered includes intermittent climatic droughts in the first half of the 21st century followed by severe persistent droughts in the second half of the 21st century. These climatic droughts do not yield a valley-wide operational drought but do cause reduced surface water deliveries and increased groundwater abstractions that may cause additional land subsidence, reduced water for riparian habitat, or changes in flows at the Sacramento-San Joaquin River Delta. The method developed here can be used to explore conjunctive use adaptation options and hydrologic risk assessments in regional hydrologic systems throughout the world.

Water Resources Research↗
USGS2012

RIP-ET: A riparian evapotranspiration package for MODFLOW-2005

A new evapotranspiration package for the U.S. Geological Survey's groundwater-flow model, MODFLOW, is documented. The Riparian Evapotranspiration Package (RIP-ET) provides flexibility in simulating riparian and wetland transpiration not provided by the Evapotranspiration (EVT) or Segmented Function Evapotranspiration (ETS1) Packages for MODFLOW 2005. This report describes how the RIP-ET package was conceptualized and provides input instructions, listings and explanations of the source code, and an example. Traditional approaches to modeling evapotranspiration (ET) processes assume a piecewise linear relationship between ET flux and hydraulic head. The RIP-ET replaces this traditional relationship with a segmented, nonlinear dimensionless curve that reflects the eco-physiology of riparian and wetland ecosystems. Evapotranspiration losses from these ecosystems are dependent not only on hydraulic head, but on the plant types present. User-defined plant functional groups (PFGs) are used to elucidate the interaction between plant transpiration and groundwater conditions. Five generalized plant functional groups based on transpiration rates, plant rooting depth, and water tolerance ranges are presented: obligate wetland, shallow-rooted riparian, deep-rooted riparian, transitional riparian and bare ground/open water. Plant functional groups can be further divided into subgroups (PFSGs) based on plant size, density or other characteristics. The RIP-ET allows for partial habitat coverage and mixtures of plant functional subgroups to be present in a single model cell. RIP-ET also distinguishes between plant transpiration and bare-ground evaporation. Habitat areas are designated by polygons; each polygon can contain a mixture of PFSGs and bare ground, and is assigned a surface elevation. This process requires a determination of fractional coverage for each of the plant functional subgroups present in a polygon to account for the mixture of coverage types and resulting transpiration. The fractional cover within a cell has two components: (1) the polygonal fraction of active habitat (excluding area of bare ground, dead trees, or brush) in a cell, and (2) fraction of plant type area or bare ground area in a polygon. RIP-ET determines the transpiration rate for each plant functional group and evaporation from bare ground/open water in a cell, the total ET in the cell, and the total ET rate over the region of simulation.

Techniques and Methods↗
USGS2011

Kirschenmann Road multi-well monitoring site, Cuyama Valley, Santa Barbara County, California

The U.S. Geological Survey (USGS), in cooperation with the Water Agency Division of the Santa Barbara County Department of Public Works, is evaluating the geohydrology and water availability of the Cuyama Valley, California (fig. 1). As part of this evaluation, the USGS installed the Cuyama Valley Kirschenmann Road multiple-well monitoring site (CVKR) in the South-Main subregion of the Cuyama Valley (fig. 1). The CVKR well site is designed to allow for the collection of depth-specific water-level and water-quality data. Data collected at this site provides information about the geology, hydrology, geophysics, and geochemistry of the local aquifer system, thus, enhancing the understanding of the geohydrologic framework of the Cuyama Valley. This report presents the construction information and initial geohydrologic data collected from the CVKR monitoring site, along with a brief comparison to selected supply and irrigation wells from the major subregions of the Cuyama Valley (fig. 1).

California↗
USGS2011

MODPATH-LGR; documentation of a computer program for particle tracking in shared-node locally refined grids by using MODFLOW-LGR

The computer program described in this report, MODPATH-LGR, is designed to allow simulation of particle tracking in locally refined grids. The locally refined grids are simulated by using MODFLOW-LGR, which is based on MODFLOW-2005, the three-dimensional groundwater-flow model published by the U.S. Geological Survey. The documentation includes brief descriptions of the methods used and detailed descriptions of the required input files and how the output files are typically used. The code for this model is available for downloading from the World Wide Web from a U.S. Geological Survey software repository. The repository is accessible from the U.S. Geological Survey Water Resources Information Web page at http://water.usgs.gov/software/ground_water.html. The performance of the MODPATH-LGR program has been tested in a variety of applications. Future applications, however, might reveal errors that were not detected in the test simulations. Users are requested to notify the U.S. Geological Survey of any errors found in this document or the computer program by using the email address available on the Web site. Updates might occasionally be made to this document and to the MODPATH-LGR program, and users should check the Web site periodically.

Techniques and Methods↗
USGS2010

Simulation and analysis of conjunctive use with MODFLOW's farm process

The extension of MODFLOW onto the landscape with the Farm Process (MF-FMP) facilitates fully coupled simulation of the use and movement of water from precipitation, streamflow and runoff, groundwater flow, and consumption by natural and agricultural vegetation throughout the hydrologic system at all times. This allows for more complete analysis of conjunctive use water-resource systems than previously possible with MODFLOW by combining relevant aspects of the landscape with the groundwater and surface water components. This analysis is accomplished using distributed cell-by-cell supply-constrained and demand-driven components across the landscape within " water-balance subregions" comprised of one or more model cells that can represent a single farm, a group of farms, or other hydrologic or geopolitical entities. Simulation of micro-agriculture in the Pajaro Valley and macro-agriculture in the Central Valley are used to demonstrate the utility of MF-FMP. For Pajaro Valley, the simulation of an aquifer storage and recovery system and related coastal water distribution system to supplant coastal pumpage was analyzed subject to climate variations and additional supplemental sources such as local runoff. For the Central Valley, analysis of conjunctive use from different hydrologic settings of northern and southern subregions shows how and when precipitation, surface water, and groundwater are important to conjunctive use. The examples show that through MF-FMP's ability to simulate natural and anthropogenic components of the hydrologic cycle, the distribution and dynamics of supply and demand can be analyzed, understood, and managed. This analysis of conjunctive use would be difficult without embedding them in the simulation and are difficult to estimate a priori. Journal compilation ?? 2010 National Ground Water Association. No claim to original US government works.

Ground Water↗
USGS2009

The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1

The ability to dynamically simulate the integrated supply-and-demand components of irrigated agricultural is needed to thoroughly understand the interrelation between surface water and groundwater flow in areas where the water-use by vegetation is an important component of the water budget. To meet this need, the computer program Farm Process (FMP1) was updated and refined for use with the U.S. Geological Survey's MODFLOW-2005 groundwater-flow model, and is referred to as MF2005-FMP2. The updated program allows the simulation, analysis, and management of nearly all components of human and natural water use. MF2005-FMP2 represents a complete hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for water consumption of irrigated agriculture, but also of urban use, and of natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, the MF2005-FMP2 has the unique set of capabilities to simulate and analyze historical, present, and future conditions. MF2005-FMP2 facilitates the analysis of agricultural water use where little data is available for pumpage, land use, or agricultural information. The features presented in this new version of FMP2 along with the linkages to the Streamflow Routing (SFR), Multi-Node Well (MNW), and Unsaturated Zone Flow (UZF) Packages prevents mass loss to an open system and helps to account for 'all of the water everywhere and all of the time'. The first version, FMP1 for MODFLOW-2000, is limited to (a) transpiration uptake from unsaturated root zones, (b) on-farm efficiency defined solely by farm and not by crop type, (c) a simulation of water use and returnflows related only to irrigated agriculture and not also to non-irrigated vegetation, (d) a definition of consumptive use as potential crop evapotranspiration, (e) percolation being instantly recharged to the uppermost active aquifer, (f) automatic routing of returnflow from runoff either to reaches of tributary stream segments adjacent to a farm or to one reach nearest to the farm's lowest elevation, (g) farm-well pumping from cell locations regardless of whether an irrigation requirement from these cells exists or not, and (h) specified non-routed water transfers from an undefined source outside the model domain. All of these limitations are overcome in MF2005-FMP2. The new features include (a) simulation of transpiration uptake from variably saturated, fully saturated, or ponded root zones (for example, for crops like rice or riparian vegetation), (b) definition of on-farm efficiency not only by farm but also by crop, (c) simulation of water use and returnflow from non-irrigated vegetation (for example, rain-fed agriculture or native vegetation), (d) use of crop coefficients and reference evapotranspiration, (e) simulation of the delay between percolation from farms through the unsaturated zone and recharge into the uppermost active aquifer by linking FMP2 to the UZF Package, (f) an option to manually control the routing of returnflow from farm runoff to streams, (g) an option to limit pumping to wells located only in cells where an irrigation requirement exists, and (h) simulation of water transfers to farms from a series of well fields (for example, recovery well field of an aquifer-storage-and-recovery system, ASR). In addition to the output of an economic budget for each farm between irrigation demand and supply ('Farm Demand and Supply Budget' in FMP1), a new output option called 'Farm Budget' was created for FMP2, which allows the user to track all physical flows into and out of a water accounting unit at all times. Such a unit can represent individual farms, farming districts, natural areas, or urban areas. The example model demonstrates the application of MF2005-FMP2 with delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in de

Techniques and Methods↗
USGS2007

Climate variability controls on unsaturated water and chemical movement, High Plains aquifer, USA

Responses in the vadose zone and groundwater to interannual, interdecadal, and multidecadal climate variability have important implications for groundwater resource sustainability, yet they are poorly documented and not well understood in most aquifers of the USA. This investigation systematically examines the role of interannual to multidecadal climate variability on groundwater levels, deep infiltration (3–23 m) events, and downward displacement (>1 m) of chloride and nitrate reservoirs in thick (15–50 m) vadose zones across the regionally extensive High Plains aquifer. Such vadose zone responses are unexpected across much of the aquifer given a priori that unsaturated total-potential profiles indicate upward water movement from the water table toward the root zone, mean annual potential evapotranspiration exceeds mean annual precipitation, and millennia-scale evapoconcentration results in substantial vadose zone chloride and nitrate reservoirs. Using singular spectrum analysis (SSA) to reconstruct precipitation and groundwater level time-series components, variability was identified in all time series as partially coincident with known climate cycles, such as the Pacific Decadal Oscillation (PDO) (10–25 yr) and the El Niño/Southern Oscillation (ENSO) (2–6 yr). Using these lag-correlated hydrologic time series, a new method is demonstrated to estimate climate-varying unsaturated water flux. The results suggest the importance of interannual to interdecadal climate variability on water-flux estimation in thick vadose zones and provide better understanding of the climate-induced transients responsible for the observed deep infiltration and chemical-mobilization events. Based on these results, we discuss implications for climate-related sustainability of the High Plains aquifer.

Vadose Zone Journal↗
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↗
USGS2005

Combined use of flowmeter and time-drawdown data to estimate hydraulic conductivities in layered aquifer systems

The vertical distribution of hydraulic conductivity in layered aquifer systems commonly is needed for model simulations of ground-water flow and transport. In previous studies, time-drawdown data or flowmeter data were used individually, but not in combination, to estimate hydraulic conductivity. In this study, flowmeter data and time-drawdown data collected from a long-screened production well and nearby monitoring wells are combined to estimate the vertical distribution of hydraulic conductivity in a complex multilayer coastal aquifer system. Flowmeter measurements recorded as a function of depth delineate nonuniform inflow to the wellbore, and this information is used to better discretize the vertical distribution of hydraulic conductivity using analytical and numerical methods. The time-drawdown data complement the flowmeter data by giving insight into the hydraulic response of aquitards when flow rates within the wellbore are below the detection limit of the flowmeter. The combination of these field data allows for the testing of alternative conceptual models of radial flow to the wellbore.

Ground Water↗
USGS2004

A methodology to asess relations between climatic variability and variations in hydrologic time series in the southwestern United States

A new method for frequency analysis of hydrologic time series was developed to facilitate the estimation and reconstruction of individual or groups of frequencies from hydrologic time-series and facilitate the comparison of these isolated time-series components across data types, between different hydrologic settings within a watershed, between watersheds, and across frequencies. While climate-related variations in inflow to and outflow from aquifers have often been neglected, the development and management of ground-water and surface-water resources has required the inclusion of the assessment of the effects of climatic variability on the supply and demand and sustainability of use. The regional assessment of climatic variability of surface-water and ground-water flow throughout the southwestern United States required this new systematic method of hydrologic time-series analysis. To demonstrate the application of this new method, six hydrologic time-series from the Mojave River Basin, California were analyzed. The results indicate that climatic variability exists in all the data types and are partially coincident with known climate cycles such as the Pacific Decadal Oscillation and the El Nino-Southern Oscillation. The time-series also indicate lagged correlations between tree-ring indices, streamflow, stream base flow, and ground-water levels. These correlations and reconstructed time-series can be used to better understand the relation of hydrologic response to climatic forcings and to facilitate the simulation of streamflow and ground-water recharge for a more realistic approach to water-resource management. Published by Elsevier B.V.

Journal of Hydrology↗
USGS2004

Documentation of the Santa Clara Valley regional ground-water/surface-water flow model, Santa Clara Valley, California

The Santa Clara Valley is a long, narrow trough extending about 35 miles southeast from the southern end of San Francisco Bay where the regional alluvial-aquifer system has been a major source of water. Intensive agricultural and urban development throughout the 20th century and related ground-water development resulted in ground-water-level declines of more than 200 feet and land subsidence of as much as 12.7 feet between the early 1900s and the mid-1960s. Since the 1960s, Santa Clara Valley Water District has imported surface water to meet growing demands and reduce dependence on ground-water supplies. This importation of water has resulted in a sustained recovery of the ground-water flow system. To help support effective management of the ground-water resources, a regional ground-water/surface-water flow model was developed. This model simulates the flow of ground water and surface water, changes in ground-water storage, and related effects such as land subsidence. A numerical ground-water/surface-water flow model of the Santa Clara Valley subbasin of the Santa Clara Valley was developed as part of a cooperative investigation with the Santa Clara Valley Water District. The model better defines the geohydrologic framework of the regional flow system and better delineates the supply and demand components that affect the inflows to and outflows from the regional ground-water flow system. Development of the model includes revisions to the previous ground-water flow model that upgraded the temporal and spatial discretization, added source-specific inflows and outflows, simulated additional flow features such as land subsidence and multi-aquifer wellbore flow, and extended the period of simulation through September 1999. The transient-state model was calibrated to historical surface-water and ground-water data for the period 1970–99 and to historical subsidence for the period 1983–99. The regional ground-water flow system consists of multiple aquifers that are grouped into upper- and lower-aquifer systems. Ground-water inflow occurs as natural recharge in the form of streamflow infiltration and areal infiltration of precipitation along stream channels, artificial recharge from infiltration of imported water at recharge ponds and along selected stream channels, and leakage along selected transmission pipelines. Ground-water outflow occurs as evapotranspiration, stream base flow, discharge through pumpage from wells, and subsurface flow to the San Francisco Bay. The geohydrologic framework of the regional ground-water flow system was represented as six model layers. The hydraulic properties were redefined on the basis of cell-based lithologic properties that were delineated in terms of aggregate thicknesses of coarse-grained, fine-grained, and mixed textural categories. The regional aquifer systems also are dissected by several laterally extensive faults that may form at least partial barriers to the lateral flow of ground water. The spatial extent of the ground-water flow model was extended and refined to cover the entire Santa Clara Valley, including the Evergreen subregion. The temporal discretization was refined and the period of simulation was extended to 1970–99. The model was upgraded to MODFLOW-2000 (MF2K) and was calibrated to fit historical ground-water levels, streamflow, and land subsidence for the period 1970–99. The revised model slightly overestimates measured water levels with an root-mean-square error of -7.34 feet. The streamflow generally shows a good match on gaged creeks and rivers for flows greater than 1.2 cubic feet per second. The revised model also fits the measured deformation at the borehole extensometer site located near San Jose within 16 to 27 percent and the extensometer site near Sunnyvale within 3 percent of the maximum measured seasonal deformation for the deepest extensometers. The total ground-water inflow and outflow of about 225,500 acre-feet per year (acre-ft/yr) for the period 1970–89 and of about 205,300 acre-feet per year for the period for the period 1970–99 is comparable with that of the previous model, 207,200 acre-ft/yr for the period 1970–89. Overall the simulated net change in storage increased by about 189,500 acre-ft/yr for the entire period of simulation, which represents about one and a half years of the 1970–99 average pumping. The changes in ground-water flow and storage generally reflect the major climate cycles and the additional importation of water by Santa Clara Valley Water District, with the basin in recovery since the drought of the late 1980s and early 1990s. The average total recharge rate, from natural and artificial recharge and from streamflow infiltration for the revised model for the entire simulation period 1970–99, was about 157,100 acre-ft/yr, which represents about 59 percent of the inflow to the ground-water flow system. The average rate of artificial recharge of about 77,600 acre-ft/yr represents about 30 percent of the inflow to the ground-water flow system. The average pumpage for the entire 29.75-year simulation period is about 133,400 acre-ft/yr and represents about 69 percent of the outflow from the ground-water flow system. Most of the simulated recharge infiltrates and flows through the uppermost layers (i.e. model layers 1 and 3) of the aquifer system. Most of the water that flows to the deeper model layers is occurring through wellbores, with wellbore flow representing 19 percent of the total ground-water inflow between model layers.

California↗
USGS2004

Geologic, water-chemistry, and hydrologic data from multiple-well monitoring sites and selected water-supply wells in the Santa Clara Valley, California, 1999-2003

To better identify the three-dimensional geohydrologic framework of the Santa Clara Valley, lithologic, geologic, geophysical, geomechanical, hydraulic, and water-chemistry data were collected from eight ground-water multiple-well monitoring sites constructed in Santa Clara County, California, as part of a series of cooperative studies between the U.S. Geological Survey and the Santa Clara Valley Water District. The data are being used to update and improve the three-dimensional geohydrologic framework of the basin and to address issues related to water supply, water chemistry, sequence stratigraphy, geology, and geological hazards. This report represents a compilation of data collected from 1999 to 2003, including location and design of the monitoring sites, cone penetrometer borings, geologic logs, lithologic logs, geophysical logs, core analysis, water-chemistry analysis, ground-water-level measurements, and hydraulic and geomechanical properties from wells and core samples. Exploratory cone penetrometer borings taken in the upper 17 to 130 feet at six of the monitoring sites identified the base of Holocene as no deeper than 75 feet in the central confined area and no deeper than 35 feet in the southern unconfined areas of the valley. Generalized lithologic characterization from the monitoring sites indicates about four to six different aquifer units separated by relatively fine-grained units occur within the alluvial deposits shallower than 860 feet deep. Analysis of geophysical logs indicates that coarse-grained units varied in thickness between 10 and 25 feet in the southeastern unconfined area of the valley and between 50 and 200 feet in the south-central and southwestern areas of the valley. Deviations from temperature-gradient logs indicate that the majority of horizontal ground-water flow occurs above a depth of 775 feet in the south central and above 510 feet in the southeastern areas of the valley. Bulk physical properties from more than 1,150 feet of core samples indicate an average primary-wave velocity of about 5,515 feet per second, a bimodal distribution of density between 2.19 and 2.32 grams per cubic centimeter with an average of 2.16 grams per cubic centimeter, and a magnetic susceptibility that generally ranged between 9 and 40 with an average of 22. Water-chemistry data indicate that the ground water in the alluvial aquifers generally is low in total dissolved solids and chloride and of good quality. Isotopic data indicate that water from artificial recharge is present throughout the shallower parts of the aquifer system but may not be present toward the center of the valley. The percentage of water from artificial recharge present in ground water ranges from 0 to 61 percent for water-supply wells. The age of most shallow ground water is less than 2,000 years before present, and the age of deeper ground water is as much as 39,900 years before present, as determined from carbon age dates. Initial water-level data from the multiple-well monitoring sites indicate seasonal water-level fluctuations as great as 60 feet and water-level differences between aquifers as great as 10 feet. The water-level hydrographs indicate different water-level changes and relations between aquifers in different parts of the basin. However, most of these hydrographs indicate the potential for downward water-level gradients, with lower hydraulic heads in the deeper monitoring wells. Hydraulic properties of selected new monitoring wells indicate that horizontal hydraulic conductivities range from 0.1 to 583 feet per day. Hydraulic testing of selected core samples yielded vertical hydraulic conductivity values ranging from 8 x 10-4 to 0.3 feet per day, and effective porosity values ranging from 0.21 to 0.4. Geomechanical properties estimated from one-dimensional consolidation tests of selected core samples resulted in geometric mean inelastic and elastic specific storage values of 1.5 x 10-

Scientific Investigations Report↗
USGS2002

Santa Clara Valley water district multi-aquifer monitoring-well site, Coyote Creek Outdoor Classroom, San Jose, California

The U.S. Geological Survey (USGS), in cooperation with the Santa Clara Valley Water District (SCVWD), has completed the first of several multiple-aquifer monitoring-well sites in the Santa Clara Valley. This site monitors ground-water levels and chemistry in the one of the major historic subsidence regions south of San Jose, California, at the Coyote Creek Outdoor Classroom (CCOC) (fig. 1) and provides additional basic information about the geology, hydrology, geochemistry, and subsidence potential of the upper- and lower-aquifer systems that is a major source of public water supply in the Santa Clara Valley. The site also serves as a science education exhibit at the outdoor classroom operated by SCVWD.

Open-File Report↗
USGS2002

Three-dimensional measurements of flow in uncased wells completed in basalt, Mountain Home Air Force Base, Idaho, March 2000

Several ground-water monitoring wells on the Mountain Home Air Force Base in southwestern Idaho were constructed in February 2000 to replace existing monitoring wells that became ineffective as a result of declining water levels. Upon completion of the replacement wells, borehole geophysical logs were collected, including natural gamma radiation, electromagnetic induction, caliper, fluid temperature, and resistivity. A prototype borehole acoustic doppler velocimeter (B-ADV) was used to make experimental three-dimensional measurements of lateral and vertical flow in two of the replacement wells, MW11–2 and MW3–2, each 450 feet deep, to better understand ground-water flow in the basalt underlying this area. Measurements indicated two independent flow zones in each well: unit B, from about 380 to about 415 feet below land surface, and unit C, from about 415 to about 430 feet below land surface. In each well, direction of flow in unit B was north-northwest toward Canyon Creek and, in unit C, south-southwest toward the Snake River. Measurements also indicated downward intraborehole flow in both wells. Unit B appeared to represent a local-scale flow regime; unit C appeared to represent a regional flow regime. This information suggests the existence of a complex three-dimensional hydrogeologic setting that cannot be discriminated easily on conventional waterlevel maps. Although data from only two wells are insufficient to construct a conceptual model of the ground-water flow regime, these experimental results demonstrate the capability of the B-ADV to obtain detailed flow measurements that, combined with data from other types of geophysical logs, discrete measurements of hydraulic head, and water chemistry, would aid in future studies and management of the ground-water resources, including contaminant transport and remediation.

Idaho↗
USGS1996

Postaudit of head and transmissivity estimates and ground-water flow models of Avra Valley, Arizona

Ground water from regional alluvial-aquifer systems is the main source of water in the alluvial basins of Arizona, such as Avra Valley. Ground-water flow models are used to assess ground-water availability and the effects of development on the regional ground-water resources. A postaudit of regional-head and transmissivity estimates and the ground-water flow models of Avra Valley was used to evaluate potential errors in the distribution of aquifer properties and recharge that can cause predictive errors in ground-water models. Simulations of predevelopment conditions in 1940 and historical development conditions for 1960-79 provided the basis of comparison for assessing predictive errors of historical conditions for two regional ground-water flow models. Potential errors in the estimation of the regional-head and transmissivity and alternate conceptual models were compared with an existing calibrated two-layer flow model for predevelopment (1940) and developed conditions (1940-85). Measured heads can be subdivided into a north-central region and a region south of the basin constriction. A more variable regional-head surface typical of developed aquifer systems was indicated by kriged developed heads (1985) with about 50 percent more uncertainty than predevelopment heads (1940). Incorporating heads from adjacent basins at the ground-water inflow and outflow regions reduced uncertainty in kriged heads for these boundary areas. Universal cokriging of heads with the strongly correlated land-surface altitudes may improve regional-head estimates and model comparisons where head data are sparse. Local transmissivity estimates can be subdivided into northern and south-central regions that are distributed along the valley axis and the Santa Cruz River. Regional geostatistical estimates of transmissivity, which are based solely on local estimates, are low in the northern part of Avra Valley and are high in the south-central part when compared with the head-conditioned model-derived estimates. These differences may be related to a systematic bias between aquifer-test conditions and methods of aquifer-test analysis. Cokriging transmissivities with specific capacity and silt-and-clay content provided the least uncertainty of all kriged estimates. Predictive errors for the Avra Valley model are the result of a different combination of factors that become significant in the simulation of ground-water flow for the periods representing predevelopment, historical development, and future development conditions. Predictive errors for simulation of predevelopment conditions are caused by potential systematic errors in estimates of local transmissivity, uncertainty in long-term mountain-front recharge, and uncertainty in predevelopment heads along the margins of the basin where recharge and transmissivity estimates are constrained by heads during model calibration. Analog-model historical predictions of future development indicate changes to 1985 were as much as 50 to 100 feet different from actual declines that were caused by errors in the spatial distribution and not the total amount of estimated future pumpage. Predictive errors for simulation of historical development (1960-79) appear to be caused to a greater extent by combined errors in estimates of transmissivity and storage properties and to a lesser extent by estimates of net withdrawal and subsidence. Comparison of the two digital models resulted in differences in transmissivity of as much as 30,000 feet squared per day and differences in specific yield of as much as 0.1. In combination with some differences in net withdrawal, these model-parameter differences resulted in local differences in change in storage of as much as 4,000 acre-feet per square mile and are equivalent to historical predictive errors in water levels of as much as 40 feet. Areas with no differences in model parameters yield comparable simulated water-level declines that are similar to measured declines. The pattern of differences in transmissivity and storage parameters are similar to differences between model-derived estimates conditioned on heads and related geostatistical estimates derived from aquifer-test estimates. A postaudit analysis of alternate conceptual models was explored on the basis of well-by-well comparisons of reductions in mean error and variance, and through the use of standardized calibration-error maps for predevelopment heads (1940) and developed heads (1985). Calibration-error maps provide a useful tool for exploring the spatial structure of model errors and the relative adequacy of model fit that is not available from traditional methods of model comparison. Calibration-error maps indicate estimated heads were too high in the southern part of Avra Valley, and estimated heads were too low in the northern part. Increased transmissivity in the southern part of the lower model layer; decreased hydraulic conductivity in the southwestern part of the upper layer; reduced ground-water inflow from Altar Valley; and increased recharge along the Tortolita Mountains, Tucson Mountains, and Brawley Wash yielded a significantly better model for predevelopment but not for developed conditions (1940-85). This may indicate that alternate conceptual models are different for different time periods or require analysis of time-varying model parameters for developed conditions, such as climatically variable recharge. Predictive errors for future simulations (1986-2025) also could potentially include errors of more than 40 ft from omission of subsidence from the simulation of regional ground-water flow in Avra Valley. Further refinement of the changing conceptual model of an aquifer system under continuing development and variable climate, such as Avra Valley, will require a variety of additional geophysical, geochemical, and hydraulic field data.

Arizona↗