Geology Reports⌕ Search

Geology topics

Andrew T. Leaf

Publications and source records attributed to Andrew T. Leaf.

At least 19 recordsLinked to original sources
USGS2025

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi↗
USGS2025

PEST++IES how many iterations & realizations, finding the point of diminishing returns

PEST++IES (White 2018; White et al. 2020) is widely used in the groundwater modeling community for its ability to perform computationally efficient history matching and uncertainty analysis in a highly parameterized context. One primary advantage of using an iterative ensemble smoother is that the number of model runs required per iteration depends on the number of realizations in an ensemble, not the number of parameters in each realization. However, this raises the question: what is the optimal number of realizations and iterations to use for any one model before the point of diminishing returns? Using a modified version of the Freyberg model (Freyberg 1988; Hunt et al. 2020), different parameter and observation scenarios were evaluated for four iterations and ensembles of 10, 25, 50, 100, 250, 500, 1000, and 2000 realizations. To match observations, PEST++IES altered hydraulic conductivity ( k ), both globally across the model and locally at three different pilot point densities, as well as global recharge (via a single multiplier), global river conductance, and individual well flow rates. Risk-based well capture zone results (Fienen et al. 2022a) and estimated hydraulic conductivity fields from each scenario were quantitatively and qualitatively compared against the “truth” model and its outputs. Across the cases examined, ensemble sizes of 100 to 250 realizations and two PEST++IES iterations were generally sufficient to achieve good results.

Groundwater↗
USGS2024

Simulating present and future groundwater/surface-water interactions and stream temperatures in Beaver Creek, Kenai Peninsula, Alaska

In many places, coldwater ecosystems are facing increasing pressure from anthropogenic warming. This study examined stream temperatures and the water balance in the Beaver Creek watershed on the Kenai Peninsula in south-central Alaska—an area that is experiencing rapid warming. Low-gradient streams near the Kenai coast provide important spawning and rearing habitat for salmon but may be especially vulnerable to rising temperatures, because of long residence times, inflows from abundant riparian wetlands, and reliance on groundwater discharge that may also warm, or decrease in volume with rising evapotranspiration. In recent decades, observed maximum 7-day temperatures have consistently exceeded statistical (regression-based) projections. Here we simulate total streamflows and temperatures with a physics-based model that links the Soil Water Balance, MODFLOW 6 and SNTEMP simulation codes on a 7-day timestep. The model is based on existing data and groundwater levels, instream flows, and stream temperatures collected during 2019–23. Future climate scenarios were developed for 2023–50 from downscaled climate projections. Results indicate that groundwater discharge is about 64 percent of the total streamflow during the months of May through September. Total streamflow and groundwater discharge are expected to remain similar to current conditions through 2050. Stream temperatures are expected to rise; by midcentury, near the Beaver Creek mouth the model predicts 34 to 63 additional days per year with average weekly temperatures above 13 degrees Celsius, 14 to 81 additional days with average weekly temperatures above 15 degrees Celsius, and routine exceedances of 20 degrees Celsius during the warmest periods. Projected stream temperatures vary spatially. Areas of high groundwater inflows in the lower main stem and some tributaries may be most resilient to warming air temperatures during dry conditions. During storm events, groundwater-dominated tributaries may have the coolest stream temperatures.

Alaska↗
USGS2024

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas↗
USGS2023

Automated construction of Streamflow-Routing networks for MODFLOW—Application in the Mississippi Embayment region

In humid regions with dense stream networks, surface water exerts a fundamental control on the water levels and flow directions of shallow groundwater. Understanding interactions between groundwater and surface water is critical for managing groundwater resources and groundwater-dependent ecosystems. Representing streams in groundwater models has historically been arduous and error prone. In recent years, however, all the information needed to numerically describe stream boundary conditions for a model area has become readily available online, as have robust open-source software tools for translating that information to a model grid. The SFRmaker Python package leverages geospatial capabilities in the scientific Python ecosystem to robustly automate the production of input to the Streamflow-Routing (SFR) Package of MODFLOW from the National Hydrography Dataset Plus or other hydrography data. This report documents an application of SFRmaker to automate production of SFR Package input for groundwater models within the Mississippi Embayment Regional Aquifer Study area. SFR Package input was developed in three steps: (1) preprocessing to develop a single set of grid-independent flowlines from National Hydrography Dataset Plus version 2 data; (2) setting up the SFR package from the preprocessed flowlines, and (3) correcting streambed top elevations after an initial model run. Separating the hydrography preprocessing from the construction of SFR Package input was advantageous in that it minimized the need to repeat computationally expensive geoprocessing (thereby speeding model construction) and also allowed for the curation of a single set of grid-independent SFR input data that can be used for any MODFLOW model within the study area.

Mississippi Embayment Region↗
USGS2023

Simulating groundwater flow in the Mississippi Alluvial Plain with a focus on the Mississippi Delta

The Mississippi Alluvial Plain has become one of the most important agricultural regions in the United States but relies heavily on groundwater for irrigation. On average, more than 12 billion gallons are withdrawn daily from the Mississippi River Valley alluvial aquifer. Declining groundwater levels, especially in the Delta region of northwest Mississippi and the Cache and Grand Prairie regions of eastern Arkansas, have led to concerns about future sustainability. The U.S. Geological Survey Mississippi Alluvial Plain Project is focused on quantifying the groundwater system in the alluvial plain and the response of groundwater resources to future development. A key objective of the project is to provide updated groundwater flow models supported by extensive data collection and analyses. MODFLOW 6, PEST++, and several open-source python packages were used to develop a simplified, faster running version of the Mississippi Embayment Regional Aquifer Study model that can provide boundary conditions for local inset models, including the Mississippi Delta model described in this report. An automated workflow was used for model construction, history matching, and development of baseline future climate scenarios. The models incorporate information from a Soil-Water-Balance code simulation of the terrestrial water balance, metering-based estimates of water use from thousands of wells, measured and estimated streamflow and stages, and the largest airborne electromagnetic survey flown to date in the United States. Baseline scenarios for the Mississippi Delta under potential future climates were constructed using recharge, surface runoff and irrigation pumping forcings from a future version of the Soil-Water-Balance model, driven by downscaled temperature and precipitation output from 10 general circulation model simulations, including high and moderate carbon emissions pathways. Results indicate a complex water balance that varies in time and space in terms of the terrestrial recharge, stream leakage, and regional groundwater flow components, which are affected by seasonal forcings, human activity, and alluvial geomorphology. The general circulation model outputs indicate a continued rise in average temperatures but no clear precipitation trend. Increased crop water demand is anticipated from the higher temperatures, resulting in increased irrigation withdrawals to sustain current levels of irrigated agriculture. Simulated drawdowns in groundwater levels at the mid-21st century vary greatly. Under moderate or wet climate scenarios, and in parts of the aquifer that are well connected to surface water, little to no additional drawdown is anticipated. Under dry or warm scenarios, drawdowns of as much as 10 meters or more are possible in parts of the aquifer that are relatively disconnected from surface water. Under dry or warm scenarios, the portion of the Delta with greater than 60 feet of saturated thickness could be reduced from near 100 percent currently (2018) to 80–90 percent by mid-century. Future simulations with the model could include alternative management scenarios to identify options for improving groundwater sustainability. The automated model construction workflows are designed to facilitate regular updating, making this a “living” framework that the Mississippi Department of Environmental Quality and other stakeholders can use for adaptive management going forward.

Arkansas, Louisiana, Mississippi↗
USGS2023

System-scale airborne electromagnetic surveys in the lower Mississippi River Valley support multidisciplinary applications

The lower Mississippi River Valley spans over 200,000 square kilometres in parts of seven states, encompassing areas of critical groundwater supplies, natural hazards, infrastructure, and low-lying coastal regions. From 2018 - 2022, the U.S. Geological Survey acquired over 82,000 line-kilometres of airborne electromagnetic, radiometric, and magnetic data over this region to provide comprehensive and systematic information about subsurface geologic and hydrologic properties that support multiple scientific and societal interests. Most of the data were acquired on a regional grid of west-east flight lines separated by 3 - 6 kilometres; however, several high-resolution inset grids with line spacing as close as 200 m were acquired in targeted areas of interest. Approximately 8,000 line-kilometres were acquired along streams and rivers to characterise the potential for surface water-groundwater connection, and another 6,000 line-kilometres were acquired along the Mississippi and Arkansas River levees to characterise this critical infrastructure. Here, we present a summary of the data along with several examples of how they are being used to inform regional groundwater model development, inferences of groundwater salinity, identification of faults in the New Madrid seismic zone, and levee infrastructure.

lower Mississippi River Valley↗
USGS2022

Modflow-setup: Robust automation of groundwater model construction

In an age of both big data and increasing strain on water resources, sound management decisions often rely on numerical models. Numerical models provide a physics-based framework for assimilating and making sense of information that by itself only provides a limited description of the hydrologic system. Often, numerical models are the best option for quantifying even intuitively obvious connections between human activities and water resource impacts. However, despite many recent advances in model data assimilation and uncertainty quantification, the process of constructing numerical models remains laborious, expensive, and opaque, often precluding their use in decision making. Modflow-setup aims to provide rapid and consistent construction of MODFLOW groundwater models through robust and repeatable automation. Common model construction tasks are distilled in an open-source, online code base that is tested and extensible through collaborative version control. Input to Modflow-setup consists of a single configuration file that summarizes the workflow for building a model, including source data, construction options, and output packages. Source data providing model structure and parameter information including shapefiles, rasters, NetCDF files, tables, and other (geolocated) sources to MODFLOW models are read in and mapped to the model discretization, using Flopy and other general open-source scientific Python libraries. In a few minutes, an external array-based MODFLOW model amenable to parameter estimation and uncertainty quantification is produced. This paper describes the core functionality of Modflow-setup, including a worked example of a MODFLOW 6 model for evaluating pumping impacts to a lake in central Wisconsin, United States.

Wisconsin↗
USGS2022

Simulation of regional groundwater flow and groundwater/lake interactions in the Central Sands, Wisconsin

A multiscale, multiprocess modeling approach was applied to the Wisconsin Central Sands region in central Wisconsin to quantify the connections between the groundwater system, land use, and lake levels in three seepage lakes in Waushara County, Wisconsin: Long and Plainfield (The Plainfield Tunnel Channel Lakes), and Pleasant Lakes. A regional groundwater-flow model, the Newton Raphson formulation of the U.S. Geological Survey modular finite-difference flow model groundwater-modeling package (MODFLOW-NWT), centered on the lakes, was used to extend regional surface-water boundaries to provide boundary conditions for two focused inset models, in the hydrologic simulation modeling package (MODFLOW 6), at higher resolution around the lakes. Land use and groundwater use were simulated at a regional scale using the Soil Water Balance model, which provided recharge and water-use boundary conditions for the MODFLOW models. Agricultural irrigation is the primary groundwater use in the area. Land and groundwater-use scenarios representing no irrigation, current (2018) irrigation, and potential future irrigation were simulated with the groundwater-flow model and the lake levels over a 38-year representative climate period.

Wisconsin↗
USGS2022

Areas contributing recharge to priority wells in valley-fill aquifers in the Neversink River and Rondout Creek drainage basins, New York

In southeastern New York, the villages of Ellenville, Wurtsboro, Woodridge, the hamlet of Mountain Dale, and surrounding communities in the Neversink River and Rondout Creek drainage basins rely on wells that pump groundwater from valley-fill glacial aquifers for public water supply. Glacial aquifers are vulnerable to contamination because they are highly permeable and have a shallow depth to water table. To protect the quality of these water resources, water managers need accurate information about the areas that contribute recharge to production wells that pump from these aquifers. The New York State Department of Environmental Conservation and the New York State Department of Health designated eight priority wells in this region for which water supply protection is of primary concern. The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation and the New York State Department of Health, began an investigation in 2019 with the general objectives of (1) improving understanding of regional groundwater-flow system, (2) delineating areas contributing recharge to eight priority production wells, and (3) quantifying the uncertainty of these contributing areas in a probabilistic way that can be used to inform decision-making related to priority well source-water protection. To complete these objectives, a MODFLOW 6 groundwater model was created encompassing the eight priority wells and the surrounding flow system, which includes parts of the Neversink River and Rondout Creek Basins in Sullivan County and Ulster County, New York. The model was built using Python tools (such as flopy , modflow-setup , and sfrmaker ) that facilitate transparent and repeatable model development using existing datasets. The model parameters were estimated with a stepwise approach using an iterative ensemble smoother implementation of the Parameter ESTimation software PEST++ (version 5.0.0). We evaluated initial “best guess” parameter bounds with a prior Monte Carlo analysis. Results of the first prior Monte Carlo analysis were used to make informed adjustments to model parameter bounds (typically resulting in expanded bounds), and a second prior Monte Carlo analysis was run to identify improved ranges for model parameters during history matching. The history matching effort produced an ensemble of parameter values for the groundwater-flow model that spans the range of values within prior uncertainty bounds. The ensemble is informed by the historical observation data, within a reasonable range of uncertainty on those observations. This history-matched ensemble was used in a particle tracking Monte Carlo analysis to delineate the areas contributing recharge to priority wells. The groundwater-flow and particle tracking (MODPATH7) models were run once for each ensemble member. Deterministic contributing areas computed for each ensemble member were aggregated to produce maps showing the probability that a location contributes recharge to priority wells. Finally, the particle tracking Monte Carlo analysis was repeated for six pumping scenarios, representing a wide range of possible pumping levels, to incorporate uncertainty in future pumping rates related to population growth or other management decisions. Increasing pumping rates generally led to larger contributing recharge areas and larger areas of high probability that a location contributes recharge to priority wells. These maps show the overall uncertainty of the areas contributing recharge to priority wells in the study area and provide a tool for risk-based decision making for protection of well source water.

New York↗
USGS2021

Risk-based wellhead protection decision support: A repeatable workflow approach

Environmental water management often benefits from a risk-based approach where information on the area of interest is characterized, assembled, and incorporated into a decision model considering uncertainty. This includes prior information from literature, field measurements, professional interpretation, and data assimilation resulting in a decision tool with a posterior uncertainty assessment accounting for prior understanding and what is learned through model development and data assimilation. Model construction and data assimilation are time consuming and prone to errors, which motivates a repeatable workflow where revisions resulting from new interpretations or discovery of errors can be addressed and the analyses repeated efficiently and rigorously. In this work, motivated by the real world application of delineating risk-based (probabilistic) sources of water to supply wells in a humid temperate climate, a scripted workflow was generated for groundwater model construction, data assimilation, particle-tracking and post-processing. The workflow leverages existing datasets describing hydrogeology, hydrography, water use, recharge, and lateral boundaries. These specific data are available in the United States but the tools can be applied to similar datasets worldwide. The workflow builds the model, performs ensemble-based history matching, and uses a posterior Monte Carlo approach to provide probabilistic capture zones describing source water to wells in a risk-based framework. The water managers can then select areas of varying levels of protection based on their tolerance for risk of potential wrongness of the underlying models. All the tools in this workflow are open-source and free, which facilitates testing of this repeatable and transparent approach to other environmental problems.

Groundwater↗
USGS2021

Risk-based decision-support groundwater modeling for the lower San Antonio River Basin, Texas, USA

A numerical surface-water/groundwater model was developed for the lower San Antonio River Basin to evaluate the responses of low base flows and groundwater levels within the basin under conditions of reduced recharge and increased groundwater withdrawals. Batch data assimilation through history matching used a simulation of historical conditions (2006-2013); this process included history-matching to groundwater levels and base-flow estimates at several gages, and was completed in a high-dimensional (highly parameterized) framework. The model was developed in an uncertainty framework such that parameters, observations, and scenarios of interest are envisioned stochastically as distributions of potential values. Results indicate that groundwater contributions to surface water during periods of low flow may be reduced from 6% to 25% with a corresponding 25% reduction in recharge and a 25% increase in groundwater pumping over an 8-year planning period. Furthermore, results indicate groundwater-level reductions in some hydrostratigraphic units are more likely than in other hydrostratigraphic units over an 8-year period under drought conditions with the higher groundwater withdrawal scenario.

Texas↗
USGS2021

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↗
USGS2021

Appendix C: Central sands lakes study technical report: Modeling documentation

This report provides the necessary documentation of the numerical models developed for the Central Sands Lake study in central Wisconsin and will be included as a technical appendix in the report to the Wisconsin State Legislature by the Wisconsin Department of Natural Resources (WDNR) in response to 2017 Wisconsin Act 10. This legislation directed WDNR to determine whether existing and potential groundwater withdrawals are causing or are likely to cause significant reduction of mean seasonal water levels at Pleasant Lake, Long Lake, and Plainfield Lake (s. 281.34(7m)(2)(b), Wis. Stats.) in Waushara County, Wisconsin. To evaluate the potential hydrologic connection between groundwater withdrawals and the nearby study lakes, hydrologic models were created that focused on the lakes of interest and yet were large enough to cover a broad enough region to extend to the major hydrologic boundaries of the natural flow system. The areas near the lakes require finer-scale grid discretization (or spacing) to better represent the lakes and streams in the model, but also need to cover a large enough area to include the groundwater withdrawal locations that have the potential to cause reduction in water levels in the lakes. To accomplish these goals, three groundwater models were created: a regional model extending to major hydrologic boundaries; and two inset models, inheriting boundaries from the regional model but focused near the lakes. Each of the inset models, in turn, included a detailed area close to the lakes surrounded by an area at the same spatial scale as the regional model (Figure 1). To support WDNR in evaluating the connection between groundwater withdrawals and lake levels, a representative time period was required over which to compare land use with and without irrigated agriculture and for WDNR to evaluate potential lake stage and flux changes related to irrigated agriculture. WDNR chose the climate period of 1981-2018 to be representative of a typical period and provided two land use scenarios—one with no irrigated agriculture and one with assumed crop rotations similar to current conditions—to simulate with groundwater models to, then, compare lake responses with. As a result, simulations over this climate record are not intended to recreate the history of 1981-2018 because land use changed over that time. These runs are, instead, intended to provide a basis on which to compare land use with and without irrigation-related groundwater withdrawals based on the current arrangement of land use and a varied climatic record. Groundwater withdrawals focused on irrigated-agriculture-related water use because greater than 95% of groundwater withdrawal in the two inset models around the study lakes is for irrigated agriculture water use. The period of 2012-2018 was used for parameter estimation (synonymously referred to as “history matching”) for the groundwater models. This time period was chosen because it includes the most complete water use records to simulate groundwater withdrawals. History matching was performed using groundwater elevations, lake stages, and streamflow observations over the 2012-2018 time period and processed observations derived from those raw data. Climatic data were incorporated into the model using a soil-water balance approach. A soil water balance model was constructed at the scale of the regional groundwater model to both calculate recharge based on land use and climate, and in the long-term climate-period runs, to estimate water use required by irrigated agriculture to apply as well boundary conditions in the groundwater model in the absence of reported water use values over that period.

Wisconsin↗
USGS2021

Hydrogeology and simulation of groundwater flow in Columbia County, Wisconsin

This report describes the regional hydrogeology and groundwater resources of Columbia County, Wisconsin, and documents a regional groundwater flow model developed for the county. Regional hydrostratigraphic units include the unlithified aquifer, the upper bedrock aquifer, and the Elk Mound aquifer. The unlithified aquifer consists of deposits that range in composition from sand and gravel outwash and stream deposits to silty, sandy till. This aquifer is less than 25 ft thick in much of eastern Columbia County, but consists of permeable sand and gravel extending to over 250 ft in depth in the Wisconsin River valley bottom. The upper bedrock aquifer consists of Ordovician and upper Cambrian sedimentary formations, including sandstone, siltstone and dolomitic strata. The upper bedrock aquifer underlies the unlithified aquifer in eastern portions of the County, but is absent to the west, where these formations are largely eroded. The contact between the Tunnel City Group and Wonewoc Formation (Top of Elk Mound Group) forms the base of the upper bedrock aquifer. Bedding plane fractures are common to this aquifer, although only a portion of the observed fractures appear to be hydraulically active. The upper bedrock aquifer is a significant source of groundwater at a regional scale. Measurements of hydraulic head showed a difference of several feet across the bottom of this aquifer to the underlying Wonewoc sandstone, indicating that the basal facies of the Tunnel City Group functions as an aquitard separating the upper bedrock aquifer from the Elk Mound aquifer. Conditions vary considerably within this aquifer, depending on the local lithostratigraphy. For example, where present, the St. Lawrence Fm. and fine-grained intervals of the Tunnel City Group may be locally-extensive aquitards. The Elk Mound aquifer consists of Cambrian sandstone of the Wonewoc, Eau Claire, and Mount Simon Formations. It is thin to absent in several locations but ranges up to 600 ft in thickness over much of southern Columbia County. The variation in thickness is due in large part to the irregular topography of the underlying Precambrian crystalline rock, which generally serves as the base of the groundwater system. In neighboring counties, a fine-grained facies within the Eau Claire Fm. acts as a regionally extensive aquitard, referred to as the Eau Claire aquitard. Much of the data collected and compiled for this study suggest that shale or dolomite within the Eau Claire Fm., which is the equivalent of the Eau Claire aquitard, occurs only within southwestern Columbia County. There is little to no evidence of the Eau Claire aquitard over most of the county. Where the dolomite and shale are absent, the Elk Mound aquifer is relatively homogenous and does not include a mappable aquitard. A three-dimensional steady-state flow model presented here represents long-term, average conditions in the regional groundwater system since about 1970. The model was constructed with the U.S. Geological Survey’s MODFLOW-NWT code; it has six layers with a uniform grid of 300 ft x 300 ft cells. Layers 1 and 2 simulate the unlithified aquifer and layer 3 represents the upper bedrock aquifer. The Elk Mound aquifer is simulated by layers 4, 5 and 6, representing the Wonewoc, Eau Claire, and Mount Simon Formations, respectively. The model extends beyond the boundaries of Columbia County to ensure that hydrologic conditions simulated within the County are consistent with regional conditions. Recharge to the groundwater flow model is based on results from a GIS-based soil-water-balance model. Recharge was simulated with the unsaturated zone flow (UZF) package in MODFLOW. This approach is particularly useful for quantifying groundwater discharge to riparian wetlands because UZF tracks recharge that would lead to the simulated water table exceeding the land surface (represented by the top of model layer 1) and reroutes it to nearby stream segments. The model includes pumping from 256 wells, and 178 of these are located within Columbia County. Pumping totaled about 28 million gallons per day (mgd) on average since 1970, with 7.2 mgd of the withdrawal from within the County. Model calibration was performed with the PEST parameter estimation code. Calibration targets included approximately 3,900 head measurements and 91 stream flow measurements. Four vertical-head differences across hydrogeologic units, calculated from data collected during packer testing in wells in Columbia County, were also used in model calibration. Results from the calibrated model provide a groundwater balance for the region. About 83 percent of groundwater originates as recharge to the water table, 12 percent comes from leakage from streams, and about 5 percent of the groundwater flows into the model domain from surrounding areas. About 95 percent of the simulated groundwater discharges to steams and other surface water features, about 3 percent flows across model boundaries to surrounding areas of the groundwater system, and pumping accounts for 2 percent of discharge. Simulated flow paths are relatively local, from recharge in upland areas to discharge in nearby streams and wetlands. The model has many potential applications, including simulation of the effects of existing or proposed high-capacity wells, estimating the zone of contribution for these wells, and understanding relationships between surface water and groundwater. Future refinements to the model, such as incorporating new information about the extent and hydraulic characteristics of the Tunnel City Group, will improve its utility in understanding advective flow between the upper bedrock and Elk Mound aquifers. If seasonal or annual variations in the groundwater system are of interest, this steady-state model could be brought into a transient mode.

Wisconsin↗
USGS2020

A distributed temperature sensing investigation of groundwater discharge to Haskell Lake, Lac du Flambeau Reservation, Wisconsin, July 27–August 1, 2016

Haskell Lake is a shallow, 89-acre drainage lake in the headwaters of the Squirrel River, on the Lac du Flambeau Reservation in northern Wisconsin. Historically, this lake was an important producer of wild rice for the Lac du Flambeau Band of Lake Superior Chippewa Indians (LDF Tribe); but, beginning in the late 1970s, the rice began to diminish and by the late 1990s, the lake no longer had harvestable stands. Restoring wild rice to Haskell Lake is a long-term priority for the LDF Tribe. A first step towards that effort is the cleanup of a petroleum-contamination plume in the shallow aquifer upgradient of the northern end of the lake. Knowledge of the downgradient extent of the plume and the locations where contaminated water is discharging to the lake is needed to inform cleanup efforts. A cooperative study between the U.S. Geological Survey and the LDF Tribe was initiated to characterize the distribution of groundwater discharge to Haskell Lake in the areas downgradient of the contamination plume. A fiber optic distributed temperature sensing system was used to monitor temperatures at the sediment-water interface for a 7-day period in July and August 2016. Challenges during the investigation included data storage and power supply limitations, maintenance of calibration baths, accurate location of the cable in space, cable placement in weeds and soft sediment, the confounding effects of solar radiation, and contamination of the data by multiple sources of instrument noise. The problem of instrument noise was overcome by solving the fiber optic distributed temperature sensing calibration equation for two parameters that describe temporal variation in the source laser and the photon detectors that observe the backscatter. Early morning temperatures, when the influence of solar radiation via direct warming of the sediment-water interface is minimized, were used to evaluate groundwater discharge, similar to other studies. The results indicate a persistent, horizontal variation in temperature of as much as 5.5 degrees Celsius across the study area, with cooler temperatures interpreted to indicate spatially discrete preferential groundwater discharge. Results of the study can be used to determine locations for collecting lakebed pore water samples to better define the extent of contamination discharging to the lake.

Wisconsin↗
USGS2020

Hydrology of Haskell Lake and investigation of a groundwater contamination plume, Lac du Flambeau Reservation, Wisconsin

Haskell Lake is a shallow, 89-acre drainage lake in the headwaters of the Squirrel River, on the Lac du Flambeau Reservation in northern Wisconsin. The lake has long been valued by the Lac du Flambeau Band of Lake Superior Chippewa Indians (LDF Tribe) for abundant wild rice and game fish. In recent decades, however, wild rice has mostly disappeared from the lake and the fishery has declined. A petroleum contamination plume discovered in the 1990s in the shallow aquifer upgradient from the northern end of the lake poses a threat to the ecological health of the lake and the aquifer, which is the sole drinking water source for nearby residents and businesses. Understanding of the lake’s hydrology is important to the LDF Tribe as they seek to restore wild rice and maintain the ecological health of the Haskell Lake/Tower Creek watershed. An improved understanding of lithology in the area of the contamination plume, documentation of a contamination pathway from groundwater in the plume source area to Haskell Lake, and an understanding of the plume extent beneath the lake are needed to advance remediation efforts. Evaluation of the fraction of groundwater discharge that is contaminated relative to the overall lake water budget is desired as a first step towards determining the extent of ecological effects from the plume. A cooperative study between the U.S. Geological Survey and the LDF Tribe was initiated to quantify the lake water budget and the sources of water to the lake, to provide a rough estimate of the maximum quantity of groundwater discharge to the lake that may be contaminated, and to improve the conceptual understanding of the plume extent and subsurface materials in the area of contamination. The results of this study can help inform natural resource management of the Haskell Lake/Tower Creek watershed, including planned wild rice restoration and cleanup of the contaminant plume. During 2016–17, field data on lake and groundwater levels, gradients, fluxes, and subsurface lithology were collected using a variety of techniques that ranged from basic measurement of water levels and streamflows to distributed temperature sensing, vertical temperature profiling, and several shallow geophysical methods. The data were used to inform a MODFLOW–NWT model that simulated the contributing groundwatershed, including the water budget for Haskell Lake and Tower Creek using the Lake, Streamflow-Routing, and Unsaturated Zone-Flow Packages. Particle tracking with the MODFLOW solution (using MODPATH 6) was used to improve understanding of the downgradient extent of the contamination plume, estimate groundwater flux through the plume area, and delineate the groundwater contributing area (groundwatershed) for the lake/creek system. Linear uncertainty estimates for model results were computed during model parameter estimation using the software package PEST++. Results indicate groundwater discharge along the perimeter of Haskell Lake, with groundwater accounting for about 22 (± 11.5) percent of the lake water budget. Field data and particle tracking results indicate discharge of the entire contamination plume to Haskell Lake. Although the exact locations where contaminated groundwater enters the lake are unknown, the downgradient extent of the plume beneath Haskell Lake is likely limited to within about 700 feet from the shore. Groundwater flux through the plume accounts for at most about 1.4 percent of total groundwater discharge to Haskell Lake, or about 0.3 percent of the lake water budget. Most groundwater discharging to Haskell Lake and Tower Creek originates as terrestrial recharge. A lesser amount originates in or passes through neighboring lakes, including Buckskin, Crawling Stone, Broken Bow, Tippecanoe, and Jerms Lakes, as well as several unnamed kettles. The average age of simulated groundwater discharge to the lake is about 20 years.

Wisconsin↗