Geology topics

Wesley R. Henson

Publications and source records attributed to Wesley R. Henson.

17 recordsLinked to original sources

USGS2026

Improving estimates of spatially distributed agricultural land use using pesticide data

Accurate, spatially resolved agricultural land-use data are essential for estimating water demand, evaluating agricultural production, and understanding long-term changes in managed landscapes. However, existing land-use datasets are often temporally discontinuous, particularly for historical periods, and remotely sensed classifications face persistent challenges distinguishing spectrally similar crop types, capturing sub-annual rotations, and separating active from abandoned plantings. We developed the California Pesticide Use Reporting Land-Use Estimator (CalPUR-LUE), which converts grower-reported pesticide-application records into annual estimates of land use by crop class and cropped acreage at approximately 1-square-mile resolution across California. CalPUR-LUE performance was evaluated for 2018 using complementary spatial, categorical, and quantitative comparisons with California Department of Water Resources (CADWR) remotely sensed land-use data and county-level agricultural economic reports. At the county-meridian-township-range-section (COMTRS) scale, CalPUR-LUE and CADWR acreage estimates showed substantial proportional agreement (log-transformed Lin’s concordance correlation coefficient = 0.77) with a mean proportional bias of −15%. Crop-class agreement was high when CalPUR-LUE identified agricultural activity, with user’s accuracy exceeding 83% for all evaluated crop classes. Overall categorical agreement was 67% (Cohen’s kappa = 0.56). Grain crop underestimation contributed substantially to disagreement; however, there is evidence of CADWR systematic overestimation of grain during 2016–2019. County-level agricultural reports further corroborated CalPUR-LUE estimates and, in several cases, aligned more closely with CalPUR-LUE than with CADWR. These results demonstrate that pesticide-use reporting provides an independent and complementary source of agricultural land-use information that can corroborate contemporary estimates and reconstruct annual spatially distributed land use where historical data are limited.

California↗

USGS2026

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes substantially to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides substantial volumes of agricultural products not just for California but for the United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface-water and groundwater resources in the valley. This situation has resulted in water management challenges in the Salinas Valley that generally relate to the distribution of the water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. Historically, to deal with the distribution issue, water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface-water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface-water supplies are insufficient. The availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand for groundwater and surface-water resources throughout the Salinas Valley has resulted in undesirable effects from unsustainable water use, such as surface-water depletion, groundwater-level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface-water and groundwater resources. To meet water demands and reduce the undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and climate variability. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resources Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley hydrogeologic system called the Salinas Valley System Model. The geologic framework is known as the Salinas Valley Geologic Framework and was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The surface-water model is called the Salinas Valley Watershed Model and simulates the Salinas River watershed. Monthly surface-water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

California↗

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↗

USGS2025

Estimating irrigation consumptive use for the conterminous United States: Coupling satellite-sourced estimates of actual evapotranspiration with a national hydrologic model

Irrigation consumptive use is crucial information for water resource management and this information is generally not available for most regional and national assessments. From 2000 to 2020, daily estimates of consumptive use were produced for the conterminous U.S. at the 12-digit hydrologic unit code watershed level. Using the Simplified Surface Energy Balance model within the OpenET tool on Google Earth Engine, actual evapotranspiration estimates were derived from Landsat land surface temperature data. These estimates were mapped to the National Hydrologic Model (NHM) spatial units for irrigated areas defined by the Landsat Irrigation Dataset (LANID). The NHM allows for simultaneous estimation of consumptive use and effective precipitation. Results were compared to the 2015 USGS water use report, revealing annual maximum difference in consumptive use of 37 cm, with an average absolute difference of 13 cm and a bias of 2.7 %. Comparisons to field resolution irrigation withdrawals reported by growers in 4 major irrigation regions across the nation were made using efficiencies reported in the literature. These comparisons resulted in root mean square errors that ranged between 10 and 18 cm per year and biases of −18.2 to 15.8 %. Significant differences in these different reported estimates represent challenges in accurate data availability and inconsistent methods and definitions used in assessments. Annual consumptive use in the U.S. can vary by 37 % due to droughts, particularly in regions more heavily dependent on precipitation. Significant increases in consumptive use were noted over the 20-year period in western states, especially the Colorado River basin from 2010 to 2020, exacerbating water supply declines.

conterminous United States↗

USGS2025

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes significantly to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides significant volumes of agricultural products not just for California but the entire United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface water and groundwater resources in the valley. This has resulted in water management challenges in the Salinas Valley that are predominantly related to distribution of water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. To deal with the distribution issue, historically water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface water supplies are insufficient. Availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand on groundwater and surface water resources throughout the Salinas Valley has resulted in undesirable effects of unsustainable water use, such as surface water depletion, groundwater level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface water and groundwater resources. To meet water demands and reduce undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and a variable climate. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resource Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley Hydrogeologic system called the Salinas Valley System Model. The Salinas Valley Geologic Framework was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The Salinas Valley Watershed Model simulates the entire Salinas River watershed. Monthly surface water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

EarthArXiv↗

USGS2024

Evaluating seawater intrusion forecast uncertainty under climate change in the Pajaro Valley, California

Climate change and climate variability impacts such as rising sea levels have the potential to exacerbate seawater intrusion and the strain on coastal freshwater resources in already stressed groundwater basins such as those in the Pajaro Valley groundwater basin, California. Regional hydrologic models are often coupled with climate projections to forecast future hydrologic conditions and inform adaptive resources management strategies. However, there is high uncertainty in the future projections of water resources due to uncertainties from downscaling global general circulation models (GCMs) to local scale climate change projections, future land use changes, and the inherent uncertainty of developed hydrologic models. Future climate projections and the magnitude of their influence on modeled hydrologic drivers are highly variable. Therefore, to develop a forecast model, an ensemble of different projections can be used to capture a wider range of basin responses and the associated uncertainties in the modeled forecasts. Understanding the reliability and uncertainty of forecasts is important for developing climate adaptation strategies such as developing protective thresholds, particularly at the basin scale where the impacts are felt, and adaptation is implemented. To demonstrate this, an uncertainty analysis of groundwater level and seawater intrusion forecasts for the Pajaro Valley groundwater basin was performed using an ensemble of three future climate projections with the Pajaro Valley Integrated Hydrologic Model (PVIHM) and the first-order second moment (FOSM) method. FOSM uncertainty analysis of hydrologic forecasts across a multi-GCM climate ensemble provides an upper and lower bound of potential impacts of climate change on sustainability targets related to mitigating seawater intrusion. The groundwater level forecasts’ narrow range of variability can help policymakers in adaptation planning by constraining possible outcomes to a focused range for risk-management decisions. However, less than one-third of groundwater level forecasts met the current protection thresholds to prevent chronic lowering of groundwater. Therefore, sustainability targets may need to be reassessed. Relative to groundwater level changes, the seawater intrusion forecasts had larger uncertainty due to the GCM climate projections and the simulated hydrologic response that were compounded by the propagation of scaling and bias from the GCMs and model simplifications in simulating the coastal boundary.

California↗

USGS2021

Detection and measurement of land-surface deformation, Pajaro Valley, Santa Cruz and Monterey counties, California, 2015–18

Land-surface deformation (subsidence) caused by groundwater withdrawal is identified as an undesirable result in the Pajaro Valley Water Management Agency’s Basin Management Plan and California’s Sustainable Groundwater Management Act. In Pajaro Valley, groundwater provides nearly 90 percent of the total water supply. To aid the development of sustainable groundwater management criteria, the U.S. Geological Survey, in cooperation with the Pajaro Valley Water Management Agency, performed an analysis of land-surface deformation (subsidence and uplift) in Pajaro Valley for 2015–18, using Interferometric Synthetic Aperture Radar and continuous Global Positioning System methods. Land-surface deformation results were then compared with subsurface geology and groundwater altitudes to better understand the hydromechanical response of the coastal aquifer system. The results indicate the land surface is generally stable with only small magnitudes (less than 1 inch) of seasonal land-surface deformation (subsidence in the summer and uplift in the winter) during 2015–18. During this time, the largest magnitude of land-surface deformation was less than 2 inches of subsidence and was localized in one area just north of the city limits of Watsonville, California. Groundwater altitudes during 2015–18 demonstrated seasonal variability and annual to multi-annual increases after reaching historical lows by the mid-1990s. The small magnitudes of land-surface deformation coupled with groundwater-altitude increases in most areas indicate that the subsidence likely is largely elastic and recoverable. The Corralitos-Pajaro Valley groundwater basin contains fine-grained (clay) sediments that have the potential for permanent aquifer-system compaction and resultant land subsidence. However, groundwater altitudes throughout the Pajaro Valley have increased above historical lows, and observed increases in groundwater altitudes coincided with changes in groundwater management activities. Observed relations between groundwater management activities and groundwater altitudes indicate that management of groundwater supplies could minimize the potential for permanent land-surface deformation in Pajaro Valley.

California↗

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

Spatially distributed denitrification in a karst springshed

Karst spring measurements assess biogeochemical processes occurring within groundwater contributing areas to springs (springsheds) but can only provide aggregated information. To better understand spatially distributed processes that comprise these aggregated measures, we investigated aquifer denitrification evidence in groundwater wells ( n = 16) distributed throughout a springshed in the Upper Floridan aquifer in northern Florida. Aquifer geochemistry, nitrate isotopes, and dissolved gases were compared against similar measurements at the spring outlet to evaluate spatial heterogeneity of denitrification evidence in relation to land surface–aquifer connectivity. Sample locations spanned spatial variation in recharge processes (i.e., diffuse vs. focused recharge) and proximity to sources of denitrification reactants (e.g., wetlands). Although no distinct spatial pattern in denitrification was uncovered, excess dissolved N 2 gas measurements were only above detection in the unconfined springshed, with some evidence of a wetland proximity effect. Measured oxidation–reduction potential and dissolved oxygen poorly predicted denitrification, indicating that measured denitrification may be occurring upgradient from sampled wells. Despite dramatic spatial chemical heterogeneity across wells, mean values for recharge nitrate concentrations (0.02 to 5.56 mg N L −1 ) and excess N 2 from aquifer denitrification (below detection to 1.37 mg N L −1 ) corresponded reasonably with mean spring outlet measurements for initial nitrate (0.78 to 1.36 mg N L −1 ) and excess N 2 (0.15 to 1.04 mg N L −1 ). Congruence between groundwater and spring measurements indicates that combining sampling at the spring outlet and across the springshed is useful for understanding spatial aquifer denitrification. However, this approach would be improved with a high‐density sampling network with transects of wells along distinct groundwater flow paths.

Florida↗

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

UZIG research: Measurement and characterization of unsaturated zone processes under wide-ranging climates and changing conditions

Unsaturated zone properties and processes are central to understanding the interacting effects of land-use change, contamination, and hydroclimate on our ability to grow food, sustain clean water supplies, and minimize loss of life and property. Advances in unsaturated zone science are being achieved through collaborations across traditional boundaries where information from biological, physical, and chemical disciplines is combined for new insights. The Unsaturated Zone Interest Group (UZIG) is an organization that exists principally to promote multidisciplinary collaborations and the sharing of ideas, expertise, and technical assets. Here we summarize key findings from 14 papers, several of which originated from a meeting convened by UZIG in 2017 at the University of Florida in Gainesville titled “Land-Use Change, Climate Change, and Hydrologic Extremes: Unsaturated Zone Responses and Feedbacks.” This special section of Vadose Zone Journal contains multidisciplinary research in three general categories relevant to measuring and understanding unsaturated zone responses to changing land uses and climate: (i) unsaturated zone properties and processes; (ii) soil–plant–atmosphere interactions; and (iii) novel field sampling devices. A strong cross-cutting theme in these papers is the value of continuous monitoring data and ways of utilizing them to discover novel hydrologic, biologic, and pedologic information. As climatic and land-use conditions change and demands for resources and stresses on ecosystems continue to intensify, it is vital to improve our fundamental understanding of the processes at work in the unsaturated zone. Toward that goal, we discuss the need for improved ground-based unsaturated zone monitoring networks.

Vadose Zone Journal↗

USGS2018

What makes a first‐magnitude spring?: Global sensitivity analysis of a speleogenesis model to gain insight into karst network and spring genesis

Often, karstic conduit network geometry is unknown. This lack of knowledge represents a significant limitation when modeling flow and solute transport in karst systems. In this study, we apply Morris Method Global Sensitivity Analysis to a speleogenesis model to identify model input parameters, and combinations thereof, that most significantly influence evolution of karst conduit networks, development of first‐magnitude springs, and resulting flow and solute transport pulse responses. Based on an idealized model of the Silver Springshed in Central Florida USA, results showed that porous matrix hydraulic conductivity and parameters that govern connectivity of vertical and horizontal preferential flow paths (proto‐conduits) are the most influential parameters. In particular, a lower porous matrix conductivity is more likely to produce a first‐order magnitude spring. For the boundary conditions assumed in this application, conduits tend to develop in low topographic regions that drained nearby high regions. Morris ensemble realizations that generated first‐magnitude springs exhibit similar flow and solute transport pulse responses at the spring vent, despite differences in network configuration. However, distributed head fields are highly spatially variable, implying substantial spatial variability among solute flow paths and travel times from the land surface to the spring across realizations.

Water Resources Research↗

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↗

USGS2018

Input data processing tools for the integrated hydrologic model GSFLOW

Integrated hydrologic modeling (IHM) encompasses a vast number of processes and specifications, variable in time and space, and development of models can be arduous. Model input construction techniques have not been formalized or made easily reproducible. Creating the input files for integrated hydrologic models requires complex GIS processing of raster and vector datasets from various sources. Developing stream network topology that is consistent with the model grid-scale digital elevation model (DEM) is important for robust simulation of surface water and groundwater exchanges. Distribution of meteorological data over the model domain is difficult in complex terrain at the model-grid scale, but is necessary for realistic simulations. As model development requires extensive GIS and computer programming expertise, the use of IHMs has mostly been limited to research groups with available financial, human, and technical resources. Here we present a series of open-source Python scripts that are combined with ESRI ArcGIS to provide a formalized technique for the parameterization and development of inputs for the readily available IHM called GSFLOW. This Python toolkit automates many of the necessary and laborious processes of parameterization, including stream network development, land coverages, and meteorological distribution over the model domain. The final products of the toolkit are PRMS ready Parameter Files, along with several input parameters for a MODFLOW model, including input for the Streamflow Routing Package. A demonstration of the toolkit is provided to illustrate its capabilities.

Environmental Modelling and Software↗

USGS2017

Nitrate reduction mechanisms and rates in an unconfined eogenetic karst aquifer in two sites with different redox potential

This study integrates push-pull tracer tests (PPTT) with microbial characterization of extracted water via quantitative polymerase chain reaction (qPCR) and reverse transcriptase qPCR (RT-qPCR) of selected functional N transformation genes to quantify nitrate reduction mechanisms and rates in sites with different redox potential in a karst aquifer. PPTT treatments with nitrate (AN) and nitrate-fumarate (ANC) were executed in two wells representing anoxic and oxic geochemical end-members. Oxic aquifer zero-order nitrate loss rates (mmol L −1 h −1 ) were similar for AN and ANC treatment, ranging from 0.03 ± 0.01 to 0.05 ± 0.01. Anoxic aquifer zero-order nitrate loss rates ranged from 0.03 ± 0.02 (AN) to 0.13 ± 0.02 (ANC). Microbial characterization indicates mechanisms influencing these rates were dissimilatory nitrate reduction to ammonium (DNRA) at the anoxic site with AN treatment, assimilatory reduction of nitrate to ammonium (ANRA) with ANC treatment in the water column at both sites, and additional documented nitrate reduction that occurred in unsampled biofilms. With carbon treatment, total numbers of microbes (16S rRNA genes) significantly increased (fourteenfold to thirtyfold), supporting stimulated growth with resulting ANRA. Decreased DNRA gene concentrations ( nrfA DNA) and increased DNRA activity ratio (nrfA-cDNA/DNA) supported the assertion that DNRA occurred in the anoxic zone with AN and ANC treatment. Furthermore, decreased DNRA gene copy numbers at the anoxic site with ANC treatment suggests that DNRA microbes in the anoxic site are chemolithoautotrophic. Increased RT-qPCR denitrification gene expression ( nirK and nirS ) was not observed in water samples, supporting that any observed NO 3 -N loss due to denitrification may be occurring in unsampled microbial biofilms.

Florida↗

USGS2013

CRT--Cascade Routing Tool to define and visualize flow paths for grid-based watershed models

The U.S. Geological Survey Cascade Routing Tool (CRT) is a computer application for watershed models that include the coupled Groundwater and Surface-water FLOW model, GSFLOW, and the Precipitation-Runoff Modeling System (PRMS). CRT generates output to define cascading surface and shallow subsurface flow paths for grid-based model domains. CRT requires a land-surface elevation for each hydrologic response unit (HRU) of the model grid; these elevations can be derived from a Digital Elevation Model raster data set of the area containing the model domain. Additionally, a list is required of the HRUs containing streams, swales, lakes, and other cascade termination features along with indices that uniquely define these features. Cascade flow paths are determined from the altitudes of each HRU. Cascade paths can cross any of the four faces of an HRU to a stream or to a lake within or adjacent to an HRU. Cascades can terminate at a stream, lake, or HRU that has been designated as a watershed outflow location.

Techniques and Methods↗