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

Characteristics of Southern California coastal aquifer systems

Most groundwater produced within coastal Southern California occurs within three main types of siliciclastic basins: (1) deep (>600 m), elongate basins of the Transverse Ranges Physiographic Province, where basin axes and related fluvial systems strike parallel to tectonic structure, (2) deep (>6000 m), broad basins of the Los Angeles and Orange County coastal plains in the northern part of the Peninsular Ranges Physiographic Province, where fluvial systems cut across tectonic structure at high angles, and (3) shallow (75-350 m), relatively narrow fluvial valleys of the generally mountainous southern part of the Peninsular Ranges Physiographic Province in San Diego County. Groundwater pumped for agricultural, industrial, municipal, and private use from coastal aquifers within these basins increased with population growth since the mid-1850s. Despite a significant influx of imported water into the region in recent times, groundwater, although reduced as a component of total consumption, still constitutes a significant component of water supply. Historically, overdraft from the aquifers has caused land surface subsidence, flow between water basins with related migration of groundwater contaminants, as well as seawater intrusion into many shallow coastal aquifers. Although these effects have impacted water quality, most basins, particularly those with deeper aquifer systems, meet or exceed state and national primary and secondary drinking water standards. Municipalities, academicians, and local water and governmental agencies have studied the stratigraphy of these basins intensely since the early 1900s with the goals of understanding and better managing the important groundwater resource. Lack of a coordinated effort, due in part to jurisdictional issues, combined with the application of lithostratigraphic correlation techniques (based primarily on well cuttings coupled with limited borehole geophysics) have produced an often confusing, and occasionally conflicting, litany of names for the various formations, lithofacies, and aquifer systems identified within these basins. Despite these nomenclatural problems, available data show that most basins contain similar sequences of deposits and share similar geologic histories dominated by glacio-eustatic sea-level fluctuations, and overprinted by syndepositional and postdepositional tectonic deformation. Impermeable, indurated mid-Tertiary units typically form the base of each siliciclastic groundwater basin. These units are overlain by stacked sequences of Pliocene to Holocene interbedded marine, paralic, fluvial, and alluvial sediment (weakly indurated, folded, and fractured) that commonly contain the historically named "80-foot sand," "200-foot sand," and "400-foot gravel" in the upper part of the section. An unconformity, cut during the latest Pleistocene lowstand (?? 18 O stage 2; ca. 18 ka), forms a major sequence boundary that separates these units from the overlying Holocene fluvial sands and gravels. Unconfined aquifers occur in amalgamated coarse facies near the bounding mountains (forebay area). These units are inferred to become lithologically more complex toward the center of the basins and coast line, where interbedded permeable and low-permeability alluvial, fluvial, paralic, and marine facies contain confined aquifers (pressure area). Coastal bounding faults limit intrabasin and/or interbasin flow in parts of many basins. ?? 2009 Geological Society of America.

Special Paper of the Geological Society of America↗

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

Simulation of Intra- or transboundary surface-water-rights hierarchies using the farm process for MODFLOW-2000

Water-rights driven surface-water allocations for irrigated agriculture can be simulated using the farm process for MODFLOW-2000. This paper describes and develops a model, which simulates routed surface-water deliveries to farms limited by streamflow, equal-appropriation allotments, or a ranked prior-appropriation system. Simulated diversions account for deliveries to all farms along a canal according to their water-rights ranking and for conveyance losses and gains. Simulated minimum streamflow requirements on diversions help guarantee supplies to senior farms located on downstream diverting canals. Prior appropriation can be applied to individual farms or to groups of farms modeled as “virtual farms” representing irrigation districts, irrigated regions in transboundary settings, or natural vegetation habitats. The integrated approach of jointly simulating canal diversions, surface-water deliveries subject to water-rights constraints, and groundwater allocations is verified on numerical experiments based on a realistic, but hypothetical, system of ranked virtual farms. Results are discussed in light of transboundary water appropriation and demonstrate the approach’s suitability for simulating effects of water-rights hierarchies represented by international treaties, interstate stream compacts, intrastate water rights, or ecological requirements.

Journal of Water Resources Planning and Management↗

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

Relations between climatic variability and hydrologic time series from four alluvial basins across the southwestern United States

Hydrologic time series of groundwater levels, streamflow, precipitation, and tree-ring indices from four alluvial basins in the southwestern United States were spectrally analyzed, and then frequency components were reconstructed to isolate variability due to climatic variations on four time scales. Reconstructed components (RCs), from each time series, were compared to climatic indices like the Pacific Decadal Oscillation (PDO), North American Monsoon (NAM), and El Nin??o-Southern Oscillation (ENSO), to reveal that as much as 80% of RC variation can be correlated with climate variations on corresponding time scales. In most cases, the hydrologic RCs lag behind the climate indices by 1-36 months. In all four basins, PDO-like components were the largest contributors to cyclic hydrologic variability. Generally, California time series have more variation associated with PDO and ENSO than the Arizona series, and Arizona basins have more variation associated with NAM. ENSO cycles were present in all four basins but were the largest relative contributors in southeastern Arizona. Groundwater levels show a wide range of climate responses that can be correlated from well to well in the various basins, with climate responses found in unconfined and confined aquifers from pumping centers to mountain fronts. ?? Springer-Verlag 2006.

Hydrogeology 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

Inferring time‐varying recharge from inverse analysis of long‐term water levels

Water levels in aquifers typically vary in response to time‐varying rates of recharge, suggesting the possibility of inferring time‐varying recharge rates on the basis of long‐term water level records. Presumably, in the southwestern United States (Arizona, Nevada, New Mexico, southern California, and southern Utah), rates of mountain front recharge to alluvial aquifers depend on variations in precipitation rates due to known climate cycles such as the El Niño‐Southern Oscillation index and the Pacific Decadal Oscillation. This investigation examined the inverse application of a one‐dimensional analytical model for periodic flow described by Lloyd R. Townley in 1995 to estimate periodic recharge variations on the basis of variations in long‐term water level records using southwest aquifers as the case study. Time‐varying water level records at various locations along the flow line were obtained by simulation of forward models of synthetic basins with applied sinusoidal recharge of either a single period or composite of multiple periods of length similar to known climate cycles. Periodic water level components, reconstructed using singular spectrum analysis (SSA), were used to calibrate the analytical model to estimate each recharge component. The results demonstrated that periodic recharge estimates were most accurate in basins with nearly uniform transmissivity and the accuracy of the recharge estimates depends on monitoring well location. A case study of the San Pedro Basin, Arizona, is presented as an example of calibrating the analytical model to real data.

Water Resources Research↗

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↗