Geology Reports⌕ Search

Geology topics

Howard W. Reeves

Publications and source records attributed to Howard W. Reeves.

At least 19 recordsLinked to original sources

Teach me how to pycap: A high-capacity well decision support tool using analytical solutions in Python

Regulatory agencies in humid temperate environments rely on timely evaluations of streamflow depletion and drawdown to protect aquatic ecosystems and existing water users. Numerical models offer detailed insights, but their complexity and time demands often preclude their practical use in rapid decision-making. We present pycap-dss, an open-source Python package that implements a suite of analytical solutions for estimating streamflow depletion and drawdown. The tool supports superposition of multiple wells and time-varying pumping, enabling cumulative impact assessments in situations with multiple wells and streams. The software is modular and extensible, allowing users to interchange solutions or add new analytical methods. A YAML-based configuration supports batch processing of multiple wells, and an optional AnalysisProject class facilitates integration with regulatory workflows. Rigorous unit and regression testing ensures computational reliability, and continuous integration supports ongoing development. We demonstrate deterministic examples of drawdown where multiple solutions are readily compared and streamflow depletion with multiple wells in the Central Sands region of Wisconsin. We also show the value of Monte Carlo analyses of streamflow depletion in the same Central Sands example, leveraging computational efficiency to evaluate the uncertainty of individual and cumulative streamflow depletion calculations from over 200 high-capacity wells.

Wisconsin↗

Evaluation of an impulse-response emulator for groundwater contaminant transport modeling

There is a significant need to develop decision support tools capable of delivering accurate representations of environmental conditions, such as ground and surface water solute concentrations, in a timely and computationally efficient manner. Such tools can be leveraged to assess a large number of potential management strategies for mitigating non-point source pollutants. Here, we assess the effectiveness of the impulse-response emulation approach to approximate process-based groundwater model estimates of solute transport from MODFLOW and MT3D over a wide range of model inputs and parameters, with the goal of assessing where in parameter space the assumptions underlying this emulation approach are valid. The impulse-response emulator was developed using the sensitivity analysis utilities in the PEST++ software suite and is capable of approximating MODFLOW/MT3D estimates of solute transport over a large portion of the parameter space tested, except in cases where the Courant number is above 0.5. Across all runs tested, the highest percent errors were at the plume fronts. These results suggest that the impulse-response approach may be suitable for emulation of solute transport models for a wide range of cases, except when high-resolution outputs are needed, or when very low concentrations at plume edges are of particular interest.

Groundwater↗

Simulation of groundwater flow at the former Badger Army Ammunition Plant, Sauk County, Wisconsin

To help support remedial efforts at the former Badger Army Ammunition Plant the U.S. Geological Survey built and calibrated a transient groundwater flow model using the Newton Raphson formulation (MODFLOW–NWT) of the U.S. Geological Survey’s modular three-dimensional finite-difference code. The model simulates the groundwater flow system at the site from 1984 to 2020. The former Badger Army Ammunition Plant is a 7,275-acre site in Sauk County, Wisconsin. The plant produced smokeless gunpower and solid rocket propellent as munitions components. Peak production periods were during World War II, the Korean War, and the Vietnam War. Subsequent groundwater contamination investigations have found four plumes at the site. A health risk assessment identified at least one contaminant of concern for human health risk present in three of the plumes: the propellant burning ground plume, the deterrent burning ground plume, and the central plume. A cooperative study began between the U.S. Army Environmental Command and U.S. Geological Survey to better understand the groundwater flow system at the former Badger Army Ammunition Plant. Field data, including aquifer tests, streamflow measurements, continuous groundwater elevations, and groundwater gradients with the Wisconsin River were collected and used to inform and calibrate the groundwater flow model. The model was used to assess the variability of the groundwater system over the study period, the components of the groundwater budget, and groundwater flow directions from identified source areas towards the Wisconsin River. Model performance assessment focused on using particle tracking to compare groundwater flowpaths that originate in the contaminant source areas to the observed plume footprints. This focus on plume behavior geometry should help constrain the advective component of a future groundwater transport model of the site.

Wisconsin↗

Can hydrological models benefit from using global soil moisture, evapotranspiration, and runoff products as calibration targets?

Hydrological models are usually calibrated to in-situ streamflow observations with reasonably long and uninterrupted records. This is challenging for poorly gage or ungaged basins where such information is not available. Even for gaged basins, the single-objective calibration to gaged streamflow cannot guarantee reliable forecasts because, as has been documented elsewhere, the inverse problem is mathematically ill-posed. Therefore, the inclusion of other observations, and the reproduction of other hydrological variables beyond streamflow, become critical components of accurate hydrological forecasting. In this study, six single- and multi-objective model calibration schemes based on different combinations of gaged streamflow, global-scale gridded soil moisture, actual evapotranspiration (ET), and runoff products are used for the calibration of a process-based hydrological model for 20 catchments located within the Lake Michigan watershed, of the Laurentian Great Lakes. Results show that the addition of gridded soil moisture to gaged streamflow in model calibration improves the ET simulation performance for most of the catchments, leading to the overall best-performing models. The monthly streamflow simulation performance for the experiments using gridded runoff products to inform the model is outperformed by those using the gaged streamflow, but the discrepancy is mitigated with increasing catchment scale. A new visualization method that effectively synthesizes model performance for the simulations of streamflow, soil moisture, and ET was also proposed. Based on the method, it is revealed that the streamflow simulation performance is relatively weak for baseflow-dominated catchments; overall, the 20 catchment models simulate streamflow and ET better than soil moisture.

Lake Michigan↗

Quantifying streamflow depletion from groundwater pumping: A practical review of past and emerging approaches for water management

Groundwater pumping can cause reductions in streamflow (“streamflow depletion”) that must be quantified for conjunctive management of groundwater and surface water resources. However, streamflow depletion cannot be measured directly and is challenging to estimate because pumping impacts are masked by streamflow variability due to other factors. Here, we conduct a management-focused review of analytical, numerical, and statistical models for estimating streamflow depletion and highlight promising emerging approaches. Analytical models are easy to implement, but include many assumptions about the stream and aquifer. Numerical models are widely used for streamflow depletion assessment and can represent many processes affecting streamflow, but have high data, expertise, and computational needs. Statistical approaches are a historically underutilized tool due to difficulty in attributing causality, but emerging causal inference techniques merit future research and development. We propose that streamflow depletion-related management questions can be divided into three broad categories (attribution, impacts, and mitigation) that influence which methodology is most appropriate. We then develop decision criteria for method selection based on suitability for local conditions and the management goal, actionability with current or obtainable data and resources, transparency with respect to process and uncertainties, and reproducibility.

Journal of the American Water Resources Associatio↗

SFRmaker and Linesink-Maker: Rapid construction of streamflow routing networks from hydrography data

Groundwater models have evolved to encompass more aspects of the water cycle, but the incorporation of realistic boundary conditions representing surface water remains time-consuming and error-prone. We present two Python packages that robustly automate this process using readily available hydrography data as the primary input. SFRmaker creates input for the MODFLOW SFR package, while Linesink-maker creates linesink string input for the GFLOW analytic element program. These programs can reduce weeks or even months of manual effort to a few minutes of execution time, and carry the added advantages of reduced potential for error, improved reproducibility and facilitation of step-wise modeling through reduced dependency on a particular conceptual model or discretization. Two real-world examples at the county to multi-state scales are presented.

Groundwater↗

National-scale reservoir thermal energy storage pre-assessment for the United States

The U.S. Geological Survey is performing a pre-assessment of the cooling potential for reservoir thermal energy storage (RTES) in five generalized geologic regions (Basin and Range, Coastal Plains, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Reservoir models are developed for the metropolitan areas of eight cities (Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon) so that computed metrics can be compared to evaluate RTES potential across diverse climates, geologic settings, and physiography. Permeable, semi-confined/confined units that underlie more-utilized aquifers and contain low-quality groundwater are selected for each city. Energy storage metrics are computed for the anticipated total thickness of stratigraphy for which RTES might be feasible, including estimated required well spacing, thermal storage capacity, and thermal recovery efficiency over time. Falta et al. (2016) showed that for a modern 25,000 square-foot (2,323 square-meter), two-story office building, cooling needs exceed heating demand for almost every region of the country. We therefore use Falta et al.’s cooling demand for each city as the representative RTES stress condition for metric computation, allowing comparisons across regions. Results indicate that favorable RTES conditions exist in each region, particularly in the Illinois Basin, Coastal Plains, and Basin and Range. Thermal recovery efficiencies are very high in all regions and increase over time. The thermal storage capacity metric is most informative in the pre-assessment and underscores the importance of mapping reservoir thicknesses and porosities to permit detailed mapping of thermal storage capacity per unit area as a key RTES resource classification standard. This assessment provides a basic understanding of the RTES potential in several metropolitan areas and geologic regions throughout the United States and will aid further evaluation of national RTES efficacy.

Conference Paper↗

Continuing progress toward a national assessment of water availability and use

Executive Summary The Omnibus Public Land Management Act of 2009 (Public Law 111—11) was passed into law on March 30, 2009. Subtitle F, also known as the SECURE Water Act, calls for the establishment of a “national water availability and use assessment program” within the U.S. Geological Survey (USGS). The USGS issued the first report on the program in 2013. Program progress over the period 2013–17 is reported herein to fulfill the requirement to inform Congress on implementation of the national water availability and use assessment program, also referred to as the USGS National Water Census (the Water Census). Much work has been accomplished during 2013–17 on producing water budgets for the nation, a goal USGS outlined in its first report on program progress to Congress. The USGS has completed three geographic focus area studies and has begun three others. Work has advanced on nationwide efforts in streamflow analysis, groundwater assessment and research, evapotranspiration studies, water use, environmental water science, and drought science. The USGS works with Federal and non-Federal agencies, universities, and other organizations to ensure that the information can be aggregated with other types of water-availability and socioeconomic information, such as data on food and energy production. The USGS has also made great strides in measures for delivering data and information on the Water Census to stakeholders and the public. Much work remains to be accomplished for the Nation to have a comprehensive, ongoing Water Census. In this report, the USGS lays out activities to be accomplished in the next 5 years (2017–22), based upon current funding levels. These include selecting new focus area studies, conducting hydrologic modeling to complete water budgets for the conterminous United States, expanding groundwater modeling efforts, mapping a national classification system for environmental water science, and developing an inventory of interbasin water transfers. All of these steps are necessary in order for the Water Census to achieve the goals outlined by Congress in the SECURE Water Act.

Circular↗

Simulation of potential groundwater recharge for the glacial aquifer system east of the Rocky Mountains, 1980–2011, using the Soil-Water-Balance Model

An understanding of the spatial and temporal extent of groundwater recharge is critical for many types of hydrologic assessments involving water quality, contaminant transport, ecosystem health, and sustainable use of groundwater. Annual potential groundwater recharge was simulated at a 1-kilometer resolution with the Soil-Water-Balance (SWB) model for the glacial aquifer system east of the Rocky Mountains, from central Montana east to Maine, for calendar years 1980–2011. The SWB model used high resolution meteorological, land cover, and soil hydrology datasets that are nationally consistent and publicly available. The SWB model computed daily potential groundwater recharge as precipitation in excess of interception, runoff, evapotranspiration, and soil-water storage capacity. Daily potential recharge values within each year of the simulation were summed to produce annual potential recharge rates. Potential recharge as described in this report is water that infiltrates vertically below the plant rooting zone and is assumed to reach the water table. The calibrated SWB model in this report is called the glacial SWB model. Model calibration assumed that the area contributing to groundwater discharge equaled the surface watershed. The model was calibrated to stream base flows from 39 watersheds throughout the model domain that had hydrologic conditions appropriate for hydrograph separation. Base flows were calculated from daily streamflow records with the HYSEP local minimum hydrograph separation method The glacial SWB model reproduced the mean annual base-flow calibration targets well; the Nash-Sutcliffe efficiency coefficient was 0.94, and the root mean squared error was 1.28 inches per year. The glacial SWB model provides insight into the spatial and temporal variability in potential annual recharge across the glacial aquifer system. About 20 percent of the active model area had an average potential recharge rate of less than 1 inch per year. Total precipitation, total recharge, and recharge as a percentage of precipitation increased from west to east. A substantial amount of the recharge water (39 percent) entering the glacial aquifer system travels through developed (urbanized) and agricultural landscapes, which are known to cause water-quality impairments. Regional climatic events, such as the 1988 to 1989 drought, are apparent in the potential recharge time series. Potential recharge generally increased across the glacial aquifer system between 2001 and 2011. A comparison of the potential recharge from the glacial SWB model to previous broad-scale recharge estimates reveals several important considerations for future SWB modeling applications. Shifts in the overall distribution of potential recharge between separate models can be explained by methods used to generate base-flow calibration target datasets. Spatial patterns in potential recharge simulated by SWB models are strongly dependent on the data and assumptions used to assign model cells to hydrologic soil groups. A review of several SWB models used to estimate groundwater recharge (and not surface runoff) revealed that model results are most sensitive to input climatic data, followed by surface runoff (curve number) and root-zone depth parameters.

Scientific Investigations Report↗

Measuring and evaluating ecological flows from streams to regions: Steps towards national coverage

Living aquatic communities are largely determined and maintained by the volume and quality of flowing waters, both within lotic systems and in receiving waters of coastal systems. However, flow is one of the most frequently and extensively altered features of rivers and streams; alteration effects are likely to be exacerbated by climate change. Lotic systems vary and different fish species need different environmental conditions, and distinct problems are evident at various spatial scales. New synoptic flow and biological information now make it possible to evaluate the effects of altered flows throughout the Great Lakes Region at scales from the stream reach to the Region. We used estimates of river and streamflow and observed fish abundances to develop tools that specify the response of fish to alterations in those flows. We fit the logistic model to a cumulative fish abundance curve as a function of yield providing an empirical means to develop models of the response of cumulative fish abundance to flows. Response zones of yield for each species in each system type (based on size and thermal class) illustrate how criteria may be developed that can be used in decision‐making for management of flows. In our example application, we evaluate both the general response of brook trout ( Salvelinus fontinalis ) abundances (and fish diversity) to changes in flows and assess the sensitivity of each stream fish community to flow alteration. Mapping stream sensitivity to flow alteration throughout the US Great Lakes Region with a multiscale spatial framework showed how regional variability in sensitivity for any fish species or assemblage may be evaluated and provides managers with information to help determine where the best opportunities for protection or restoration of streamflows and associated communities exist. These results provide valuable tools and critical information to managers responsible for balancing water uses and maintaining high quality lotic ecosystems. These methods may be applied to any geographic region and can be extended nationally or globally, where flow, temperature, fish and landscape data are available.

Illinois, Indiana, Michigan, Minnesota, New York, ↗

Generalized hydrogeologic framework and groundwater budget for a groundwater availability study for the glacial aquifer system of the United States

The glacial aquifer system groundwater availability study seeks to quantify (1) the status of groundwater resources in the glacial aquifer system, (2) how these resources have changed over time, and (3) likely system response to future changes in anthropogenic and environmental conditions. The glacial aquifer system extends from Maine to Alaska, although the focus of this report is the part of the system in the conterminous United States east of the Rocky Mountains. The glacial sand and gravel principal aquifer is the largest source of public and self-supplied industrial supply for any principal aquifer and also is an important source for irrigation supply. Despite its importance for water supply, water levels in the glacial aquifer system are generally stable varying with climate and only locally from pumping. The hydrogeologic framework developed for this study includes the information from waterwell records and classification of material types from surficial geologic maps into likely aquifers dominated by sand and gravel deposits. Generalized groundwater budgets across the study area highlight the variation in recharge and discharge primarily driven by climate.

Glacial Aquifer System↗

Maps and grids of hydrogeologic information created from standardized water-well drillers’ records of the glaciated United States

As part of the National Water Availability and Use Program established by the U.S. Geological Survey (USGS) in 2005, this study took advantage of about 14 million records from State-managed collections of water-well drillers’ records and created a database of hydrogeologic properties for the glaciated United States. The water-well drillers’ records were standardized to be relatively complete and error-free and to provide consistent variables and naming conventions that span all State boundaries. Maps and geospatial grids were developed for (1) total thickness of glacial deposits, (2) total thickness of coarse-grained deposits, (3) specific-capacity based transmissivity and hydraulic conductivity, and (4) texture-based estimated equivalent horizontal and vertical hydraulic conductivity and transmissivity. The information included in these maps and grids is required for most assessments of groundwater availability, in addition to having applications to studies of groundwater flow and transport. The texture-based estimated equivalent horizontal and vertical hydraulic conductivity and transmissivity were based on an assumed range of hydraulic conductivity values for coarse- and fine-grained deposits and should only be used with complete awareness of the methods used to create them. However, the maps and grids of texture-based estimated equivalent hydraulic conductivity and transmissivity may be useful for application to areas where a range of measured values is available for re-scaling. Maps of hydrogeologic information for some States are presented as examples in this report but maps and grids for all States are available electronically at the project Web site (USGS Glacial Aquifer System Groundwater Availability Study, http://mi.water.usgs.gov/projects/WaterSmart/Map-SIR2015-5105.html ) and the Science Base Web site, https://www.sciencebase.gov/catalog/item/58756c7ee4b0a829a3276352 .

Alaska, Connecticut, Idaho, Illinois, Indiana, Iow↗

A semi-structured MODFLOW-USG model to evaluate local water sources to wells for decision support

In order to better represent the configuration of the stream network and simulate local groundwater-surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater-flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface-water depletion by wells) without recourse to labor-intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW-USG) allows application of regional models to address local problems. A “semi-structured” approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20-layer model with uniform 5000-foot (1524-m) lateral spacing was converted to 4 layers with 500-foot (152-m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi-structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi-structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface-water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).

Illinois, Indiana, Michigan, Wisconsin↗

Estimation of monthly water yields and flows for 1951-2012 for the United States portion of the Great Lakes Basin with AFINCH

Monthly water yields from 105,829 catchments and corresponding flows in 107,691 stream segments were estimated for water years 1951–2012 in the Great Lakes Basin in the United States. Both sets of estimates were computed by using the Analysis of Flows In Networks of CHannels (AFINCH) application within the NHDPlus geospatial data framework. AFINCH provides an environment to develop constrained regression models to integrate monthly streamflow and water-use data with monthly climatic data and fixed basin characteristics data available within NHDPlus or supplied by the user. For this study, the U.S. Great Lakes Basin was partitioned into seven study areas by grouping selected hydrologic subregions and adjoining cataloguing units. This report documents the regression models and data used to estimate monthly water yields and flows in each study area. Estimates of monthly water yields and flows are presented in a Web-based mapper application. Monthly flow time series for individual stream segments can be retrieved from the Web application and used to approximate monthly flow-duration characteristics and to identify possible trends.

Great Lakes Basin↗

U.S. Geological Survey water science strategy—Observing, understanding, predicting, and delivering water science to the Nation

Executive Summary This report expands the Water Science Strategy that began with the USGS Science Strategy, “Facing Tomorrow’s Challenges—U.S. Geological Survey Science in the Decade 2007–2017” (U.S. Geological Survey, 2007). This report looks at the relevant issues facing society and develops a strategy built around observing, understanding, predicting, and delivering water science for the next 5 to 10 years by building new capabilities, tools, and delivery systems to meet the Nation’s water-resource needs. This report begins by presenting the vision of water science for the USGS and the societal issues that are influenced by, and in turn influence, the water resources of our Nation. The essence of the Water Science Strategy is built on the concept of “water availability,” defined as spatial and temporal distribution of water quantity and quality, as related to human and ecosystem needs, as affected by human and natural influences. The report also describes the core capabilities of the USGS in water science—the strengths, partnerships, and science integrity that the USGS has built over its 134-year history.

Circular↗

Groundwater availability as constrained by hydrogeology and environmental flows

Groundwater pumping from aquifers in hydraulic connection with nearby streams has the potential to cause adverse impacts by decreasing flows to levels below those necessary to maintain aquatic ecosystems. The recent passage of the Great Lakes-St. Lawrence River Basin Water Resources Compact has brought attention to this issue in the Great Lakes region. In particular, the legislation requires the Great Lakes states to enact measures for limiting water withdrawals that can cause adverse ecosystem impacts. This study explores how both hydrogeologic and environmental flow limitations may constrain groundwater availability in the Great Lakes Basin. A methodology for calculating maximum allowable pumping rates is presented. Groundwater availability across the basin may be constrained by a combination of hydrogeologic yield and environmental flow limitations varying over both local and regional scales. The results are sensitive to factors such as pumping time, regional and local hydrogeology, streambed conductance, and streamflow depletion limits. Understanding how these restrictions constrain groundwater usage and which hydrogeologic characteristics and spatial variables have the most influence on potential streamflow depletions has important water resources policy and management implications.

Ground Water↗

Management of surface water and groundwater withdrawals to maintain environmental stream flows in Michigan

In 2008, the State of Michigan enacted legislation requiring that new or increased high-capacity withdrawals (greater than 100,000 gallons per day) from either surface water or groundwater be reviewed to prevent Adverse Resource Impacts (ARI). Science- based guidance was sought in defining how groundwater or surface-water withdrawals affect streamflow and in quantifying the relation between reduced streamflow and changes in stream ecology. The implementation of the legislation led to a risk-based system based on a gradient of risk, ecological response curves, and estimation of groundwater-surface water interaction. All Michigan streams are included in the legislation, and, accordingly, all Michigan streams were classified into management types defined by size of watershed, stream-water temperature, and predicted fish assemblages. Different streamflow removal percentages define risk-based thresholds allowed for each type. These removal percentages were informed by ecological response curves of characteristic fish populations and finalized through a legislative workgroup process. The assessment process includes an on-line screening tool that may be used to evaluate new or increased withdrawals against the risk-based zones and allows withdrawals that are not likely to cause an ARI to proceed to water-use registration. The system is designed to consider cumulative impacts of high-capacity withdrawals and to promote user involvement in water resource management by the establishment of water-user committees as cumulative withdrawals indicate greater potential for ARI in the watershed.

Michigan↗