Search USGS⌕ Search

Geology topics

Randall T. Hanson

Publications and source records attributed to Randall T. Hanson.

At least 19 recordsLinked to original sources
USGS2025

Application of Hydrologic Simulation Program—FORTRAN (HSPF) as part of an integrated hydrologic model for the Salinas Valley, California

The U.S. Geological Survey (USGS), in cooperation with the Monterey County Water Resources Agency, conducted studies to help evaluate the surface-water and groundwater resources of the Salinas Valley study area, consisting of the entire Salinas River watershed and several smaller, adjacent coastal watersheds draining into Monterey Bay. The Salinas Valley study area is a highly productive agricultural region that depends on the coordinated use of surface water and groundwater to meet demand for irrigation and public water supply. To continue to meet these demands, a better understanding of the historical water balance and the effects of water-resource development on the long-term sustainability of water resources in the Salinas Valley study area is needed.

California↗
USGS2021

Assessment of water availability in the Osage Nation using an integrated hydrologic-flow model

The Osage Nation of northeastern Oklahoma, conterminous with Osage County, covers about 2,900 square miles. The area is primarily rural with 62 percent of the land being native prairie grass, and much of the area is used for cattle ranching and extraction of petroleum and natural gas. Protection of water rights are important to the Osage Nation because of its reliance on cattle ranching and the potential for impairment of water quality by petroleum extraction. Additionally, the potential for future population increases, demands for water from neighboring areas such as the Tulsa metropolitan area, and expansion of petroleum and natural-gas extraction on water resources of this area further the need for the Osage Nation to better understand its water availability. Therefore, the U.S. Geological Survey, in cooperation with the Osage Nation, completed a hydrologic investigation to assess the status and availability of surface-water and groundwater resources in the Osage Nation. A transient integrated hydrologic-flow model was constructed using the U.S. Geological Survey fully integrated hydrologic-flow model called the MODFLOW One-Water Hydrologic Model. The integrated hydrologic-flow model, called the Osage Nation Integrated Hydrologic Model (ONIHM), was constructed and uses an orthogonal grid of 276 rows and 289 columns, and each grid cell measures 1,312.34 feet (ft; 400 meters) per side, with eight variably thick vertical layers that represented the alluvial and bedrock aquifers within the study area, including the alluvial aquifer, the Vamoosa-Ada aquifer, and the minor Pennsylvanian bedrock aquifers, and the confining units. Landscape and groundwater-flow processes were simulated for two periods: (1) the 1950–2014 period from January 1950 through September 2014 and (2) the forecast period from October 2014 through December 2099. The 1950–2014 period ONIHM simulated past conditions using measured or estimated inputs, and the forecast-period ONIHM simulated three separate potential forecast conditions under constant dry, average, or wet climate conditions using calibrated input values from the 1950–2014 period ONIHM. The 1950–2014 period ONIHM was calibrated by linking the Parameter Estimation software (PEST) with the MODFLOW One-Water Hydrologic Model. PEST uses statistical parameter estimation techniques to identify the best set of parameter values to minimize the difference between measured or estimated calibration targets and their simulated equivalent values (residuals). Tikhonov regularization and singular-value decomposition-assist features of PEST were used during the calibration process. The 1950–2014 period ONIHM was calibrated to 713 measured groundwater levels at 195 wells; 95,636 estimated monthly mean groundwater levels at 124 wells; 5,307 measured streamflows at 13 streamgages; and 8,679 simulated mean monthly streamflows at 10 streamgages extracted from a surface-water model by adjusting 231 parameters. The estimated groundwater-level observations and streamflows were included as observations to improve the spatial and temporal density of observation targets during calibration. The best set of parameter values obtained during the calibration process of the 1950–2014 model was then used as the input parameter values for the forecast model simulations. A comparison of the calibration targets to their corresponding simulated values indicated that the model adequately reproduced streamflows and groundwater levels for some streamgages and wells and underestimated streamflows and groundwater levels at other locations. Measured and simulated streamflows correlated adequately with a coefficient of determination of 0.938, as did water levels with a coefficient of determination of 0.795. The 1950–2014 period ONIHM underestimated certain groundwater levels and streamflows, but generally measured or estimated calibration targets correlated well with simulated equivalents, which indicated that the model can adequately simulate the response of the hydrologic system to stresses in the 1950–2014 and forecast periods. In the 1950–2014 period ONIHM, the calibrated mean horizontal hydraulic conductivity for layer 1 alluvial aquifer was 30.7 feet per day, and the seven lower layers had a calibrated mean horizontal hydraulic conductivity of less than 3.3 feet per day. The mean calibrated groundwater-level residual was 16.6 ft, and the mean calibrated streamflow residual of the Arkansas River at Ralston, Oklahoma, streamgage (U.S. Geological Survey station 07152500) was within 6 percent (373 cubic feet per second) of mean measured streamflow for the 1950–2014 period ONIHM. The ONIHM simulated landscape fluxes of precipitation; groundwater applied by irrigation wells; evapotranspiration from precipitation, groundwater, and irrigation; runoff from precipitation; and deep percolation from precipitation. The largest loss of water from the landscape was evapotranspiration from precipitation with a calibrated mean annual outflow of 32 inches (in.): mean annual precipitation was about 36 in. Calibrated mean annual runoff and deep percolation (recharge to the water table) rates were 4.7 inches per year (in/yr) and 0.70 in/yr, respectively, for the 1950–2014 period ONIHM. The calibrated 1950–2014 period ONIHM groundwater fluxes included net farm net recharge (calculated as the difference between the inflow of recharge to the water table and the outflow of evapotranspiration from the water table such that negative values indicate that evapotranspiration from the water table was greater than deep percolation [recharge to the water table] and vice versa). Net farm net recharge was the largest flux from the groundwater system with a mean annual net outflow of 153.4 cubic feet per second. Stream leakage was the largest flux to the groundwater system with a mean annual net inflow of 152.5 cubic feet per second, indicating that, on average, the groundwater/surface-water interaction was a “losing” system where stream water leaked into the subsurface and recharged the water table. Simulated monthly trends demonstrated that net stream leakage was the largest inflow to the groundwater-flow system for 10 of the 12 months; for the other 2 months (January and March), farm net recharge (January) and net storage (March) were the largest inflow to the groundwater-flow system. A saline groundwater interface map was created for the study and compared to the water levels from the final stress period of the 1950–2014 model to identify the presence of fresh/marginal groundwater throughout the study area. Fresh/marginal groundwater was characterized as groundwater with less than 1,500 milligrams per liter of total dissolved solids. Fresh/marginal groundwater thickness ranged from 0 to 438.2 ft within the study area. The thickest regions of fresh/marginal groundwater were in the eastern part of the study area near Sand Creek, Bird Creek, and Hominy Creek and in the Arkansas River alluvial aquifer in the region downstream from the Arkansas River at Ralston, Okla. Like the 1950–2014 model, forecast model results for the landscape indicated that transpiration from precipitation was the largest flux out of the landscape for all three forecasts, constituting 77, 73, and 58 percent of precipitation for the dry, average, and wet forecasts, respectively. The dry and average forecast landscape fluxes demonstrated similar trends and magnitudes, whereas the wet forecast landscape fluxes indicated the largest changes compared to the average forecast fluxes. Most notably, runoff increased from a mean of 1.1 and 1.6 in/yr for the dry and average forecasts, respectively, to 10 in/yr for the wet forecast. Similar changes occurred for the other wet forecast landscape fluxes. The calibrated 1950–2014 period ONIHM simulated three forecasts to assess the effects of potential climatic changes on the hydrologic system from October 2014 to December 2099. The three forecasts simulated theoretical dry, average, and wet conditions using precipitation and potential evapotranspiration datasets from selected years in the calibrated 1950–2014 period ONIHM. Annual precipitation amounts were 26.89, 35.47, and 50.73 in. for the dry, average, and wet forecasts, respectively. Groundwater-flow component forecast results indicated that stream leakage is always a net inflow to the groundwater-flow system for dry, average, and wet conditions, meaning the study area stream network is always predominantly a “losing” regime where stream water infiltrates into the underlying aquifer. Storage was only a net outflow from the groundwater-flow system and indicated a replenishment to groundwater storage that resulted in an increase in groundwater levels only during the wet forecast. Further, these gains in groundwater storage for the wet forecast occurred only during February through June. Mean fresh/marginal groundwater saturated thicknesses were 125 and 126 ft for the dry and average forecast conditions, respectively, and wet forecast average thickness was 145 ft and ranged from 0 to 443 ft. The spatial extents of fresh/marginal groundwater at the end of the dry, average, and wet forecast model periods (December 2099) did not change substantially from the end of the 1950–2014 model period (September 2014).

Kansas, Oklahoma↗
USGS2020

Precipitation runoff modeling system (PRMS) as part of an integrated hydrologic model for the Osage Nation, northeastern Oklahoma, 1915–2014

Executive Summary The Osage Nation lacks a comprehensive tribal water plan to describe the quality and quantity of water resources in the Osage Nation, a 2,304-square-mile (mi 2 ) area of rolling pastures, tallgrass prairie, and mixed woodlands in northeastern Oklahoma. A tribal water plan can be used to help manage the sustainable development of surface and groundwater resources, thereby helping to provide a better future for the Osage Nation and their neighbors, while preserving water resources for the benefit of the surrounding environment and future generations. To help meet these goals and contribute to increased knowledge of the quantity and quality of water resources and the hydrologic processes and factors affecting those resources, the U.S. Geological Survey (USGS) in cooperation with the Osage Nation began studies to evaluate the surface-water and groundwater resources of the Osage Nation. An important component of these studies is the development and application of numerical models to improve quantification and understanding of the hydrologic system. These models are needed to estimate and quantify the effects of historical and potential future water resource development for the Osage Nation. This report describes the development and application of a precipitation-runoff model, the Osage Nation watershed model (ONWM). The ONWM is needed as a component of the Osage Nation integrated hydrologic model (ONIHM). At the time of this study, the ONIHM was being developed using the USGS computer software MODFLOW-One Water Hydrologic Flow Model (MODFLOW-OWHM). The intended use of the ONIHM is to simulate all surface-water and groundwater components of the hydrologic system for a 2,905-mi 2 study area centered on the Osage Nation. The ONWM was developed using the USGS Precipitation Runoff Modeling System, version 4 (PRMS-IV) computer software, also referred to as PRMS in this report, for an 8,343-mi 2 study area in northeastern Oklahoma and southeastern Kansas, centered on and including the areas of the Osage Nation and the ONIHM. The ONWM is to be used as part of the ONIHM to provide a direct coupling with spatially and temporally varying daily climate conditions affecting the ONIHM study area. As an integral part of the ONIHM, the ONWM (1) simulates the inflow boundary conditions from tributary basins in the region outside and surrounding the ONIHM area; (2) provides estimates of spatially and temporally distributed precipitation, air temperature, potential evapotranspiration (PET), actual evapotranspiration (ET), soil moisture, recharge, and streamflow in the ONIHM area; and (3) provides a preliminary water budget for the ONIHM area and the surrounding region, including tributary drainage basins outside of and next to the ONIHM. The specific objectives of this study were to use the ONWM to (1) provide a systematic inventory of the historical distribution of water inflows from precipitation (rain or snow) falling on the land surface and flowing through the surface-water network, (2) provide a historical context of the variability and spatial and temporal distribution of these waters, and (3) provide estimates of water inflows and potential observations to the ONIHM. The application of the ONWM as a component of the ONIHM is needed for planned simulations using the ONIHM to improve the understanding of the hydrologic system and to develop a fully comprehensive water budget, including the use and movement of water across the landscape, in the surface-water network, and in groundwater aquifers under historical and potential future conditions.

Oklahoma↗
USGS2020

Rio Grande transboundary integrated hydrologic model and water-availability analysis, New Mexico and Texas, United States, and northern Chihuahua, Mexico

Changes in population, agricultural development and practices (including shifts to more water-intensive crops), and climate variability are increasing demands on available water resources, particularly groundwater, in one of the most productive agricultural regions in the Southwest—the Rincon and Mesilla Valley parts of Rio Grande Valley, Doña Ana and Sierra Counties, New Mexico, and El Paso County, Texas. The goal of this study was to produce an integrated hydrological simulation model to help evaluate water-management strategies, including conjunctive use of surface water and groundwater for historical conditions, and to support long-term planning for the Rio Grande Project. This report describes model construction and applications by the U.S. Geological Survey, working in cooperation and collaboration with the Bureau of Reclamation. This model, the Rio Grande Transboundary Integrated Hydrologic Model, simulates the most important natural and human components of the hydrologic system, including selected components related to variations in climate, thereby providing a reliable assessment of surface-water and groundwater conditions and processes that can inform water users and help improve planning for future conditions and sustained operations of the Rio Grande Project (RGP) by the Bureau of Reclamation. Model development included a revision of the conceptual model of the flow system, construction of a Transboundary Rio Grande Watershed Model (TRGWM) water-balance model using the Basin Characterization Model, and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model version 2 (referred to as MF-OWHM2). The hydrologic models were developed for and calibrated to historical conditions of water and land use, and parameters were adjusted so that simulated values closely matched available measurements (calibration). The calibrated model was then used to assess the use and movement of water in the Rincon Valley, Mesilla Basin, and northern part of the Conejos-Médanos Basin, with the entire region referred to as the “Transboundary Rio Grande” or TRG. These tools provide a means to understand hydrologic system response to the evolution of water use in the region, its availability, and potential operational constraints of the RGP. The conceptual model identified surface-water and groundwater inflows and outflows that included the movement and use of water both in natural and in anthropogenic systems. The groundwater-flow system is characterized by a layered geologic sedimentary sequence combined with the effects of groundwater pumping, operation of the RGP, natural runoff and recharge, and the application of irrigation water at the land surface that is captured and reused in an extensive network of canals and drains as part of the conjunctive use of water in the region. Historical groundwater-level fluctuations followed a cyclic pattern that were aligned with climate cycles, which collectively resulted in alternating periods of wet or dry years. Periods of drought that persisted for one or more years are associated with low surface-water availability that resulted in higher rates of groundwater-level decline. Rates of groundwater-level decline also increased during periods of agricultural intensification, which necessitated increasing use of groundwater as a source of irrigation water. Agriculture in the area was initially dominated by alfalfa and cotton, but since 1970 more water-intensive pecan orchards and vegetable production have become more common. Groundwater levels substantially declined in subregions where drier climate combined with increased demand, resulting in periods of reduced streamflows. Most of the groundwater was recharged in the Rio Grande Valley floor, and most of the pumpage and aquifer storage depletion was in Mesilla Basin agricultural subregions. A cyclic imbalance between inflows and outflows resulted in the modeled cyclic depletion (groundwater withdrawals in excess of natural recharge) of the groundwater basin during the 75-year simulation period of 1940–2014. Changes in groundwater storage can vary considerably from year to year, depending on land use, pumpage, and climate conditions. Climatic drivers of wet and dry years can greatly affect all inflows, outflows, and water use. Although streamflow and, to a minor extent, precipitation during inter-decadal wet-year periods replenished the groundwater historically, contemporary water use and storage depletion could have reduced the effects of these major recharge events. The average net groundwater flow-rate deficit for 1953–2014 was estimated to be about 1,090 acre-feet per year.

Chihuahua, New Mexico, Texas↗
USGS2020

One-Water Hydrologic Flow Model: A MODFLOW based conjunctive-use simulation software

The U.S. Geological Survey’s (USGS) Modular Ground-Water Flow Model (MODFLOW-2005) is a computer program that simulates groundwater flow by using finite differences. The MODFLOW-2005 framework uses a modular design that allows for the easy development and incorporation of new features called processes and packages that work with or modify inputs to the groundwater-flow equation. A process solves a flow equation or set of equations. For example, the central part of MODFLOW is the groundwater-flow process that solves the groundwater-flow equation; the surface-water routing process is an additional process that solves the surface-water flow equation. Packages are code related to the groundwater-flow process. For example, the subsidence package modifies the groundwater-flow process by including aquifer compaction effects on flow. With the development of new packages and processes, the MODFLOW-2005 base framework diverged into multiple independent versions designed for specific simulation needs. This divergence limited each independent MODFLOW release to its specific purpose, so that there was no longer a single, comprehensive, general-purpose hydraulic-simulation framework. The MODFLOW One-Water Hydrologic Flow Model (MF-OWHM, also informally known as OneWater) is an integrated hydrologic flow model that combines multiple MODFLOW-2005 variants in one cohesive simulation software; changes were made to enable multiple capabilities in one code. This fusion of the MODFLOW-2005 versions resulted in a simulation software that can be used to address and analyze a wide class of conjunctive-use, water-management, water-food-security, and climate-crop-water scenarios. As a second core version of MODFLOW-2005, MF-OWHM maintains backward compatibility with existing MODFLOW-2005 versions, with features that include the following: Process-based simulation. Saturated groundwater flow (three-dimensional). Surface-water flow (one- and two-dimensional). Stream and river flow. Lake and reservoir storage. Landscape simulation and irrigated agriculture. Land-use and crop simulation. Root uptake of groundwater. Actual evapotranspiration. Estimated irrigation demand. Reservoir operations. Aquifer compaction and subsidence by vertical model-grid deformation. Seawater intrusion by a sharp-interface assumption. Karst-aquifer and fractured-bedrock flow. Turbulent and laminar-pipe network flow. Unsaturated groundwater flow (one-dimensional). Internal linkages among the processes that couple hydraulic head, flow, and deformation. Redesigned code for faster simulation, increased user-input options, easier model updates, and more robust error reporting than in previous models. MF-OWHM is a MODFLOW-2005 based integrated hydrologic model that can simulate and analyze varying environmental conditions to allow for the evaluation of management options from many components of human and natural water movement through a physically based, supply and demand framework. The term “integrated,” in the context of this report, refers to the tight coupling of groundwater flow, surface-water flow, landscape processes, aquifer compaction and subsidence, reservoir operations, and conduit (karst) flow. Another benefit of this integrated hydrologic model is that models developed to run by MODFLOW-2005, MODFLOW-NWT, MODFLOW-CFP, or MODFLOW-FMP can also be simulated with MF-OWHM. At the time of this report’s publication, MF-OWHM version 2 (MF-OWHM2) does not include a direct internal simulation of snowmelt, advanced mountainous watershed rainfall-runoff simulation, detailed shallow soil-moisture accounting, or atmospheric moisture content. Atmospheric moisture may be accounted for indirectly by, optionally, specifying a pan-evaporation rate, reference evapotranspiration, and precipitation. These features are not included to ensure that simulation runtime remains short enough to enable the use of automated methods of calibrating model parameters to field observations, which typically require many simulation model runs. The MF-OWHM approach is to include as much detail as possible to simulate hydrological processes, providing the simulation runtimes remain reasonable enough to allow for robust parameter estimation and model calibration. To represent both natural and human-influenced flow, MF-OWHM integrates physically based flow processes derived from MODFLOW-2005 in a supply and demand framework. From this integration, the physically based movement of groundwater, surface water, imported water, and precipitation serve as supply to meet consumptive demands associated with irrigated and non-irrigated agriculture, natural vegetation, and urban water uses. Water consumption is determined by balancing the available water supply with water demand, leading to the concept of a demand-driven, supply-constrained simulation. The MF-OWHM Supply-and-Demand Framework is especially useful for the analysis of agricultural water use, where there are often few data available to describe changes in land-use through time, such as crop type and distribution, and the associated changes in groundwater pumpage. This framework attempts to satisfy each land-use water demand with available water supplies—that is, groundwater uptake, precipitation, and irrigation. An option provided in MF-OWHM2 is to automatically increase groundwater pumping for irrigation, which often is unknown, by the calculated residual between demand and the other available sources of supply. From large- to small-scale applications, the physically based supply and demand framework provides key capabilities for simulating and analyzing historical, current, and future conjunctive-use of surface water and groundwater. To achieve the physically based supply and demand framework, the MODFLOW-2005 standard of no inter-package and -process communication was relaxed for MF-OWHM2. Traditional MODFLOW simulation models required that all packages and processes interact through the groundwater-flow equation or by removing the water flow from the simulation domain. For example, the MODFLOW-2005 representation of a groundwater well extracts water from the groundwater-flow equation (by subtraction) and removes it from the simulation domain. This feature is available in the MF-OWHM framework, but options have been added to allow the specification of a use or destination of pumped groundwater within the model domain, for example, it can be used for irrigation, managed aquifer recharge, or return-flow to streams.

Techniques and Methods↗
USGS2019

Integrated hydrologic modeling of the Salinas River, California, for sustainable water management

The Salinas River is the largest river in California’s Central Coast region. Groundwater resources of the Salinas River basin are used to meet water supply needs, including crop irrigation and municipal water supply. Two large multipurpose reservoirs also supply irrigation and municipal water uses. Historical imbalances between supply and demand have resulted in sinking groundwater levels, seawater intrusion, regulatory actions on pumping, adjudication, and requirements for minimum in-stream fish flows. Present needs include finding replacement water supplies and improving watershed management to comply with legal mandates, adapt to future climate variability and landuse conversions, and improve environmental conditions. The Salinas Valley Integrated Hydrologic Model (SVIHM) was developed to help water managers evaluate and adjust to projected impacts on water supplies and demands in the Salinas Valley watershed caused by changes in land use, population, and climate. The SVIHM includes four modeling components: (1) the Basin Characterization Model (BCM), (2) the Hydrologic Simulation Program – FORTRAN (HSPF), (3) MODFLOW - One Water Hydrologic Model (MF-OWHM), and (4) the Surface Water Operations (SWO) package. The BCM and HSPF components compose the Salinas Valley Watershed Model (SVWM). The 4,530 square-mile (mi2) SVWM domain encompasses the entire Salinas River watershed, as well as coastal drainages adjacent to the Salinas River outflow, and includes two separate and connected HSPF model domains, the 2,540 mi2 upper Salinas River and the 1,990 mi2 lower Salinas River models. SVWM (1) simulates the water budget for the entire Salinas River basin containing both the SVIHM domain as well as the mountainous terrain of the tributary headwater areas not included in the SVIHM; and (2) was used to develop the 148 boundary inflows for the SVIHM. Simulated evapotranspiration (ET) is the largest component of the water budget after precipitation, with a 71-year average basin-wide ET of 13.9 in/yr, compared to the basin-wide average precipitation of 18.4 in/yr. Simulated ET ranges from 15 to 29 in/yr along the western side of the SVWM to less than 10 in/yr throughout the valley floor and in the southeast part of the Salinas River watershed. The simulated total 71-year average inflow to the SVIHM was 890 ft3/sec (about 640,000 acre-feet per year), with the highest average inflow of 270 ft3/sec simulated for the Nacimiento River; whereas, the simulated 71-year average streamflow at the mouth of the Salinas River was only about 190 ft3/sec, indicating that most of the streamflow generated in the Salinas River basin is lost to channel seepage. The lack of sustained baseflow causes streamflow to be highly sensitive to the temporal variability in precipitation, especially during the drier periods, and this increases the importance of developing adequate reservoir management, flow augmentation, and conjunctive water use scenarios for potential future drought periods and potentially increased temporal variability in precipitation.

California↗
USGS2018

Rio Grande transboundary integrated hydrologic model and water-availability analysis, New Mexico and Texas, United States, and Northern Chihuahua, Mexico

Errata **September 28, 2018: The purpose of a USGS Open-file report (OFR) is dissemination of information that must be released immediately to fill a public need or information that is not sufficiently refined to warrant publication in one of the other USGS series. As part of that refinement process, an error was discovered in one of the input data sets of the Rio Grande Transboundary Integrated Hydrologic Model (RGTIHM) that this OFR was based upon. The error involved the assignment of storage properties to “phantom cells.” Phantom cells are required for most variants of MODFLOW that use a structured finite-difference grid when individual stratigraphic layers are represented as separate layers. Using phantom cells is a common practice that allows separate model layers to be maintained without having to combine stratigraphic layers into equivalent model layers or to use an unstructured grid. Typically, phantom cell horizontal hydraulic conductivities and storage properties are set to a small number and vertical hydraulic conductivities are set to a number large enough to allow vertical flow between the vertically adjacent layers. In the RGTIHM, the specific storage properties of the phantom cells for the upper (RGTIHM layers 3 and 4), middle (RGTIHM layers 5 and 6), and lower (RGTIHM layers 7 and 8) members of the Santa Fe Group were inadvertently assigned a value of 1 feet -1 . The revision of these specific storage values to a small number (1.0 x 10 -09 feet -1 ) required additional trial-and-error model calibration and a new sensitivity analysis. After calibration, the overall model fit remained similar to the fit described in the OFR, but the fit for many individual features such as project water available for diversions at the American Canal and Acequia Madre improved due to the reduction in flow coming from lower layers. Overall, there is still an average net depletion of groundwater flow, and the conclusions of the report are not changed. The revised average annual groundwater flow depletion simulated for the period 1953-2014 is -1,480 acre-feet/year for the entire model region, and -3,660 acre-feet/year for the portion of the model in the United States. The final version of the model will be the basis of the USGS Scientific Investigations Report that will supersede this OFR. An updated Model Archive of RGTIHM is available upon request to the USGS California Water Science Center. The corrected version of the model WAS the basis for the USGS Scientific Investigations Report that SUPERSEDED this Open-File Report.** Abstract Changes in population, agricultural development and practices (including shifts to more water-intensive crops), and climate variability are increasing demands on available water resources, particularly groundwater, in one of the most productive agricultural regions in the Southwest—the Rincon and Mesilla Valley parts of Rio Grande Valley, Doña Ana and Sierra Counties, New Mexico, and El Paso County, Texas. The goal of this study was to produce an integrated hydrological simulation model to help evaluate water-management strategies, including conjunctive use of surface water and groundwater for historical conditions, and to support long-term planning for the Rio Grande Project. This report describes model construction and applications by the U.S. Geological Survey, working in cooperation and collaboration with the Bureau of Reclamation. This model, the Rio Grande Transboundary Integrated Hydrologic Model, simulates the most important natural and human components of the hydrologic system, including selected components related to variations in climate, thereby providing a reliable assessment of surface-water and groundwater conditions and processes that can inform water users and help improve planning for future conditions and sustained operations of the Rio Grande Project (RGP) by the Bureau of Reclamation. Model development included a revision of the conceptual model of the flow system, construction of a Transboundary Rio Grande Watershed Model (TRGWM) water-balance model using the Basin Characterization Model (BCM), and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model (referred to as One Water). The hydrologic models were developed for and calibrated to historical conditions of water and land use, and parameters were adjusted so that simulated values closely matched available measurements (calibration). The calibrated model was then used to assess the use and movement of water in the Rincon Valley, Mesilla Basin, and northern part of the Conejos-Médanos Basin, with the entire region referred to as the “Transboundary Rio Grande” or TRG. These tools provide a means to understand hydrologic system response to the evolution of water use in the region, its availability, and potential operational constraints of the RGP. The conceptual model identified surface-water and groundwater inflows and outflows that included the movement and use of water both in natural and in anthropogenic systems. The groundwater-flow system is characterized by a layered geologic sedimentary sequence combined with the effects of groundwater pumping, operation of the RGP, natural runoff and recharge, and the application of irrigation water at the land surface that is captured and reused in an extensive network of canals and drains as part of the conjunctive use of water in the region. Historical groundwater-level fluctuations followed a cyclic pattern that were aligned with climate cycles, which collectively resulted in alternating periods of wet or dry years. Periods of drought that persisted for one or more years are associated with low surface-water availability that resulted in higher rates of groundwater-level decline. Rates of groundwater-level decline also increased during periods of agricultural intensification, which necessitated increasing use of groundwater as a source of irrigation water. Agriculture in the area was initially dominated by alfalfa and cotton, but since 1970 more water-intensive pecan orchards and vegetable production have become more common. Groundwater levels substantially declined in subregions where drier climate combined with increased demand, resulting in periods of reduced streamflows. Most of the groundwater was recharged in the Rio Grande Valley floor, and most of the pumpage and aquifer storage depletion was in Mesilla Basin agricultural subregions. A cyclic imbalance between inflows and outflows resulted in the modeled cyclic depletion (groundwater withdrawals in excess of natural recharge) of the groundwater basin during the 75-year simulation period of 1940–2014. Changes in groundwater storage can vary considerably from year to year, depending on land use, pumpage, and climate conditions. Climatic drivers of wet and dry years can greatly affect all inflows, outflows, and water use. Although streamflow and, to a minor extent, precipitation during inter-decadal wet-year periods replenished the groundwater historically, contemporary water use and storage depletion could have reduced the effects of these major recharge events. The average net groundwater flow-rate deficit for 1953–2014 was estimated to be about 8,990 acre-feet per year.

New Mexico, Northern Chihuahua, Texas↗
USGS2016

Development and use of mathematical models and software frameworks for integrated analysis of agricultural systems and associated water use impacts

The development of appropriate water management strategies requires, in part, a methodology for quantifying and evaluating the impact of water policy decisions on regional stakeholders. In this work, we describe the framework we are developing to enhance the body of resources available to policy makers, farmers, and other community members in their e orts to understand, quantify, and assess the often competing objectives water consumers have with respect to usage. The foundation for the framework is the construction of a simulation-based optimization software tool using two existing software packages. In particular, we couple a robust optimization software suite (DAKOTA) with the USGS MF-OWHM water management simulation tool to provide a flexible software environment that will enable the evaluation of one or multiple (possibly competing) user-defined (or stakeholder) objectives. We introduce the individual software components and outline the communication strategy we defined for the coupled development. We present numerical results for case studies related to crop portfolio management with several defined objectives. The objectives are not optimally satisfied for any single user class, demonstrating the capability of the software tool to aid in the evaluation of a variety of competing interests.

AIMS Agriculture and Food↗
USGS2016

Integrated water flow model and modflow-farm process: A comparison of theory, approaches, and features of two integrated hydrologic models

Effective modeling of conjunctive use of surface and subsurface water resources requires simulation of land use-based root zone and surface flow processes as well as groundwater flows, streamflows, and their interactions. Recently, two computer models developed for this purpose, the Integrated Water Flow Model (IWFM) from the California Department of Water Resources and the MODFLOW with Farm Process (MF-FMP) from the US Geological Survey, have been applied to complex basins such as the Central Valley of California. As both IWFM and MFFMP are publicly available for download and can be applied to other basins, there is a need to objectively compare the main approaches and features used in both models. This paper compares the concepts, as well as the method and simulation features of each hydrologic model pertaining to groundwater, surface water, and landscape processes. The comparison is focused on the integrated simulation of water demand and supply, water use, and the flow between coupled hydrologic processes. The differences in the capabilities and features of these two models could affect the outcome and types of water resource problems that can be simulated.

California Natural Resources Agency↗
USGS2015

Examples of deformation-dependent flow simulations of conjunctive use with MF-OWHM

The dependency of surface- and groundwater flows and aquifer hydraulic properties on deformation induced by changes in aquifer head is not accounted for in the standard version of MODFLOW. A new USGS integrated hydrologic model, MODFLOW-OWHM, incorporates this dependency by linking subsidence and mesh deformation with changes in aquifer transmissivity and storage coefficient, and with flows that also depend on aquifer characteristics and land-surface geometry. This new deformation-dependent approach is being used for the further development of the integrated Central Valley hydrologic model (CVHM) in California. Preliminary results from this application and from hypothetical test cases of similar systems show that changes in canal flows, stream seepage, and evapotranspiration from groundwater (ET gw ) are sensitive to deformation. Deformation feedback has been shown to also have an indirect effect on conjunctive surface- and groundwater use components with increased stream seepage and streamflows influencing surface-water deliveries and return flows. In the Central Valley model, land subsidence may significantly degrade the ability of the major canals to deliver surface water from the Delta to the San Joaquin and Tulare basins. Subsidence can also affect irrigation demand and ET gw , which, along with altered surface-water supplies, causes a feedback response resulting in changed estimates of groundwater pumping for irrigation. This modeling feature also may improve the impact assessment of dewatering-induced land subsidence/uplift (following irrigation pumping or coal-seam gas extraction) on surface receptors, inter-basin transfers, and surface infrastructure integrity.

Proceedings of the International Association of Hy↗
USGS2015

Hydrogeologic framework of the Santa Clara Valley, California

The hydrologic framework of the Santa Clara Valley in northern California was redefined on the basis of new data and a new hydrologic model. The regional groundwater flow systems can be subdivided into upper-aquifer and lower-aquifer systems that form a convergent flow system within a basin bounded by mountains and hills on three sides and discharge to pumping wells and the southern San Francisco Bay. Faults also control the flow of groundwater within the Santa Clara Valley and subdivide the aquifer system into three subregions. After decades of development and groundwater depletion that resulted in substantial land subsidence, Santa Clara Valley Water District (SCVWD) and the local water purveyors have refilled the basin through conservation and importation of water for direct use and artificial recharge. The natural flow system has been altered by extensive development with flow paths toward major well fields. Climate has not only affected the cycles of sedimentation during the glacial periods over the past million years, but interannual to interdecadal climate cycles also have affected the supply and demand components of the natural and anthropogenic inflows and outflows of water in the valley. Streamflow has been affected by development of the aquifer system and regulated flow from reservoirs, as well as conjunctive use of groundwater and surface water. Interaquifer flow through water-supply wells screened across multiple aquifers is an important component to the flow of groundwater and recapture of artificial recharge in the Santa Clara Valley. Wellbore flow and depth-dependent chemical and isotopic data indicate that flow into wells from multiple aquifers, as well as capture of artificial recharge by pumping of water-supply wells, predominantly is occurring in the upper 500 ft (152 m) of the aquifer system. Artificial recharge represents about one-half of the inflow of water into the valley for the period 1970–1999. Most subsidence is occurring below 250 ft (76 m), and most pumpage occurs within the upper-aquifer system between 300 and 650 ft (between 91 and 198 m) below land surface. Overall, the natural quality of most groundwater in the Santa Clara Valley is good. Isotopic data indicate that artificial recharge is occurring throughout the shallower parts of the upper-aquifer system and that recent recharge (less than 50 yr old) occurs throughout most of the basin in the upper-aquifer system, but many of the wells in the center of the basin with deeper well screens do not contain tritium and recent recharge. Age dates indicate that the groundwater in the upper-aquifer system generally is less than 2000 yr old, and groundwater in the lower-aquifer system generally ranges from 16,700 to 39,900 yr old. Depth-dependent sampling indicates that wellbores are the main path for vertical flow between aquifer layers. Isotopic data indicate as much as 60% of water pumped from production wells originated as artificial recharge. Shallow aquifers not only contain more recent recharge but may be more susceptible to anthropogenic and natural contamination, as evidenced by trace occurrences of iron, nitrate, and volatile organic compounds (VOCs) in selected water-supply wells. Water-resource management issues are centered on sustaining a reliable and good-quality source of water to the residents and industries of the valley. While the basin has been refilled, increased demand owing to growth and droughts could result in renewed storage depletion and the related potential adverse effects of land subsidence and seawater intrusion. The new hydrologic model demonstrates the importance of the aquifer layering, faults, and stream channels in relation to groundwater flow and infiltration of recharge. This model provides a means to analyze water resource issues because it separates the supply and demand components of the inflows and outflows.

California↗
USGS2015

Physical subdivision and description of the water-bearing sediments of the Santa Clara Valley, California

A thick Quaternary alluvial section fills a sedimentary basin beneath the Santa Clara Valley, California, located within the San Andreas Fault system at the south end of San Francisco Bay. This section consists of an upper sequence about 1,000 feet thick containing eight sedimentary cycles and a lower fine-grained unit as thick as several hundred feet. Together these constitute the Quaternary Santa Clara Basin. The section overlies an irregular unconformity with more than 1,200 feet of relief cut into the underlying bedrock. This stratigraphy is determined through study of new wells and seismic reflection profiles, together with a sample of the many thousands of water wells in the valley. It represents a major change and improvement in understanding of the basin, particularly with regard to the upper cyclic sequence, which forms a large groundwater system that is an important resource in the San Francisco Bay region. Each of the eight sedimentary cycles consists of a coarse-grained bottom interval overlain by a fine-grained top, with the coarse bottom forming a permeable sheet that is more or less continuous around the basin and the fine top forming a similarly extensive, relatively impermeable confining layer. This stratigraphic organization contrasts with most previous views, which have considered the coarse sediment in the basin to occur as scattered, discrete lenses and (or) sinuous channel sands, all embedded in a predominantly fine-grained section. Temperature logs in several wells demonstrate that the fine cycle tops do limit vertical movement of groundwater, although this may not be the case where those tops are thin to perhaps locally absent around parts of the basin margin. Age control has been obtained from previous work, in which the sedimentary cycles were correlated with the marine oxygen isotope record and the ages of two deeper Quaternary unconformities were estimated, and from detailed paleomagnetic study of cores from the new wells by E.A. Mankinen. Despite careful search of the cores, very few fossils were found, and none that are helpful in subdividing the section. No tephra (volcanic ash) was recovered, and the few carbon samples found and dated radiometrically are limited to the upper 120 feet of the section. The upper cyclic section ranges in age from 0 to somewhat older than 718 thousand years (ka), and the lower fine-grained section lies between unconformities with estimated ages of 950 and 1500 ka. Reflections in the seismic profiles indicate that layering in the basin is subparallel to the ground surface, and this fact, together with the continuous stratigraphic detail provided by geophysical logs of the new wells, allows the confident interwell correlation required to delineate the sedimentary cycles. The sequence of layers within any one cycle tends to persist laterally between the wells in the dataset, which are spaced 1 to 3 km apart, with most changes occurring gradually. The eight cycles, in contrast, tend to differ from each other in the details of their internal organization. Maps and cross sections show the elevations of cycle boundaries and the underlying bedrock surface, the varying thicknesses of the cycles and of their fine tops and coarse bottoms, and the aggregate thickness of coarse layers in those bottom intervals. Coarse sediment is more abundant toward some parts of the basin margin and in the southern part of the basin. Cycle boundary surfaces are relatively smooth, and their shapes are consistent with having been intercycle topographic surfaces. The underlying bedrock surface has a relief of more than 1,200 feet and deepens toward the center of the basin and the west edge of the fault-bounded Evergreen Basin, which is concealed beneath the east side of the Quaternary basin. The absence of consistent abrupt changes in thicknesses or boundary elevations across the basin or in cross section indicates that the interior of the basin is largely unfaulted, with the Silver Creek strand of the San Andreas system at the west edge of the Evergreen Basin being the sole exception. The east and west margins of the Santa Clara Basin, in contrast, are marked by reverse and thrust fault systems.

California↗
USGS2015

Hindcast of water availability in regional aquifer systems using MODFLOW Farm Process

Coupled groundwater and surface-water components of the hydrologic cycle can be simulated by the Farm Process for MODFLOW (MF-FMP) in both irrigated and non-irrigated areas and aquifer-storage and recovery systems. MF-FMP is being applied to three productive agricultural regions of different scale in the State of California, USA, to assess the availability of water and the impacts of alternative management decisions. Hindcast simulations are conducted for similar periods from the 1960s to near recent times. Historical groundwater pumpage is mostly unknown in one region (Central Valley) and is estimated by MF-FMP. In another region (Pajaro Valley), recorded pumpage is used to calibrate model-estimated pumpage. Multiple types of observations are used to estimate uncertain parameters, such as hydraulic, land-use, and farm properties. MF-FMP simulates how climate variability and water-import availability affect water demand and supply. MF-FMP can be used to predict water availability based on anticipated changes in anthropogenic or natural water demands. Keywords groundwater; surface-water; irrigation; water availability; response to climate variability/change

Conference Paper↗
USGS2015

An integrated approach to conjunctive-use analysis with the one-water hydrologic flow model, MODFLOW-OWHM

The MODFLOW-2005 (MF) family of hydrologic simulators has diverged into multiple versions designed for specific needs, thus limiting their use to their respective designs. The One-Water Hydrologic Flow Model (MF-OWHM v1.0) is an integrated hydrologic flow model that is an enhanced fusion of multiple MF versions. While maintaining compatibility with existing MF versions, MF-OWHM includes: linkages for coupled heads, flows, and deformation; facilitation of self-updating models, additional observation and parameter options for higher-order calibrations; and redesigned code for faster simulations. This first release of MF-OWHM incorporates MODFLOW-2005 and the Farm Process (MF-FMP2), with new features (FMP3), combined with Local Grid Refinement (MF-LGR), Streamflow Routing (SFR), Surfacewater Routing Process (SWR), Seawater Intrusion (SWI), Riparian Evapotranspiration (RIP-ET), the Newton Formulation (MF-NWT), and more. MF-OWHM represents a complete integrated hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for consumption by agriculture and natural vegetation on the landscape, and for potable and other uses. By retaining and keeping track of the water during simulation of the hydrosphere, MF-OWHM accounts for “all of the water everywhere and all of the time.” This provides the foundation needed to address integrated hydrologic problems such as evaluation of conjunctive-use alternatives and sustainability analysis, including potential adaptation and mitigation strategies, and best management practices.

Conference Paper↗
USGS2014

Feedback of land subsidence on the movement and conjunctive use of water resources

The dependency of surface- or groundwater flows and aquifer hydraulic properties on dewatering-induced layer deformation is not available in the USGS's groundwater model MODFLOW. A new integrated hydrologic model, MODFLOW-OWHM, formulates this dependency by coupling mesh deformation with aquifer transmissivity and storage and by linking land subsidence/uplift with deformation-dependent flows that also depend on aquifer head and other flow terms. In a test example, flows most affected were stream seepage and evapotranspiration from groundwater (ET gw ). Deformation feedback also had an indirect effect on conjunctive surface- and groundwater use components: Changed stream seepage and streamflows influenced surface-water deliveries and returnflows. Changed ET gw affected irrigation demand, which jointly with altered surface-water supplies resulted in changed supplemental groundwater requirements and pumping and changed return runoff. This modeling feature will improve the impact assessment of dewatering-induced land subsidence/uplift (following irrigation pumping or coal-seam gas extraction) on surface receptors, inter-basin transfers, and surface-infrastructure integrity.

Environmental Modelling and Software↗
USGS2014

Analysis of projected water availability with current basin management plan, Pajaro Valley, California

The projection and analysis of the Pajaro Valley Hydrologic Model (PVHM) 34 years into the future using MODFLOW with the Farm Process (MF-FMP) facilitates assessment of potential future water availability. The projection is facilitated by the integrated hydrologic model, MF-FMP that fully couples the simulation of the use and movement of water from precipitation, streamflow, runoff, groundwater flow, and consumption by natural and agricultural vegetation throughout the hydrologic system at all times. MF-FMP 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” (WBS) comprised of one or more model cells that can represent a single farm, a group of farms, watersheds, or other hydrologic or geopolitical entities. Analysis of conjunctive use would be difficult without embedding the fully coupled supply-and-demand into a fully coupled simulation, and are difficult to estimate a priori. The analysis of projected supply and demand for the Pajaro Valley indicate that the current water supply facilities constructed to provide alternative local sources of supplemental water to replace coastal groundwater pumpage, but may not completely eliminate additional overdraft. The simulation of the coastal distribution system (CDS) replicates: 20 miles of conveyance pipeline, managed aquifer recharge and recovery (MARR) system that captures local runoff, and recycled-water treatment facility (RWF) from urban wastewater, along with the use of other blend water supplies, provide partial relief and substitution for coastal pumpage (aka in-lieu recharge). The effects of these Basin Management Plan (BMP) projects were analyzed subject to historical climate variations and assumptions of 2009 urban water demand and land use. Water supplied directly from precipitation, and indirectly from reuse, captured local runoff, and groundwater is necessary but inadequate to satisfy agricultural demand without coastal and regional storage depletion that facilitates seawater intrusion. These facilities reduce potential seawater intrusion by about 45% with groundwater levels in the four regions served by the CDS projected to recover to levels a few feet above sea level. The projected recoveries are not high enough to prevent additional seawater intrusion during dry-year periods or in the deeper aquifers where pumpage is greater. While these facilities could reduce coastal pumpage by about 55% of the historical 2000–2009 pumpage for these regions, and some of the water is delivered in excess of demand, other coastal regions continue to create demands on coastal pumpage that will need to be replaced to reduce seawater intrusion. In addition, inland urban and agricultural demands continue to sustain water levels below sea level causing regional landward gradients that also drive seawater intrusion. Seawater intrusion is reduced by about 45% but it supplies about 55% of the recovery of groundwater levels in the coastal regions served by the CDS. If economically feasible, water from summer agricultural runoff and tile-drain returnflows could be another potential local source of water that, if captured and reused, could offset the imbalance between supply and demand as well as reducing discharge of agricultural runoff into the National Marine Sanctuary of Monterey Bay. A BMP update (2012) identifies projects and programs that will fund a conservation program and will provide additional, alternative water sources to reduce or replace coastal and inland pumpage, and to replenish the aquifers with managed aquifer recharge in an inland portion of the Pajaro Valley.

California↗
USGS2014

Cuyama Valley, California hydrologic study: An assessment of water availability

Water resources are under pressure throughout California, particularly in agriculturally dominated valleys. Since 1949, the Cuyama Valley’s irrigated acreage has increased from 13 to 35 percent of the valley. Increased agriculture has contributed to the demand for water beyond natural recharge. The tools and information developed for this study can be used to help understand the Cuyama Valley aquifer system, an important resource of Santa Barbara County.

California↗