Geology topics

R. J. Hunt

Publications and source records attributed to R. J. Hunt.

At least 19 recordsLinked to original sources

USGS2012

Improving wetland simulations by including heat transport in groundwater flow modeling

A procedure was developed to automatically calibrate a groundwater flow and heat transport model, resulting in the estimation of hydraulic conductivity and flux across the water table in wetland systems. This paper describes differences between previous approaches and this study, and summarizes some challenges in the method implementation. The procedure was validated in a sequence of hypothetical models with known structure and parameters, and applied to a wetlands system located near Wilton, Wisconsin. Copyright ASCE 2004.

Conference Paper↗

USGS2011

Application of the Local Grid Refinement package to an inset model simulating the interaction of lakes, wells, and shallow groundwater, northwestern Waukesha County, Wisconsin

Groundwater use from shallow, high-capacity wells is expected to increase across southeastern Wisconsin in the next decade (2010-2020), owing to residential and business growth and the need for shallow water to be blended with deeper water of lesser quality, containing, for example, excessive levels of radium. However, this increased pumping has the potential to affect surface-water features. A previously developed regional groundwater-flow model for southeastern Wisconsin was used as the starting point for a new model to characterize the hydrology of part of northwestern Waukesha County, with a particular focus on the relation between the shallow aquifer and several area lakes. An inset MODFLOW model was embedded in an updated version of the original regional model. Modifications made within the inset model domain include finer grid resolution; representation of Beaver, Pine, and North Lakes by use of the LAK3 package in MODFLOW; and representation of selected stream reaches with the SFR package. Additionally, the inset model is actively linked to the regional model by use of the recently released Local Grid Refinement package for MODFLOW-2005, which allows changes at the regional scale to propagate to the local scale and vice versa. The calibrated inset model was used to simulate the hydrologic system in the Chenequa area under various weather and pumping conditions. The simulated model results for base conditions show that groundwater is the largest inflow component for Beaver Lake (equal to 59 percent of total inflow). For Pine and North Lakes, it is still an important component (equal, respectively, to 16 and 5 percent of total inflow), but for both lakes it is less than the contribution from precipitation and surface water. Severe drought conditions (simulated in a rough way by reducing both precipitation and recharge rates for 5 years to two-thirds of base values) cause correspondingly severe reductions in lake stage and flows. The addition of a test well south of Chenequa at a pumping rate of 47 gal/min from a horizon approximately 200 feet below land surface has little effect on lake stages or flows even after 5 years of pumping. In these scenarios, the stage and the surface-water outflow from Pine Lake are simulated to decrease by only 0.03 feet and 3 percent, respectively, relative to base conditions. Likely explanations for these limited effects are the modest pumping rate simulated, the depth of the test well, and the large transmissivity of the unconsolidated aquifer, which allows the well to draw water from upstream along the bedrock valley and to capture inflow from the Bark River. However, if the pumping rate of the test well is assumed to increase to 200 gal/min, the decrease in simulated Pine Lake outflow is appreciably larger, dropping by 14 percent relative to base-flow conditions.

Wisconsin↗

USGS2010

Assessment of sewer source contamination of drinking water wells using tracers and human enteric viruses

This study investigated the source, transport, and occurrence of human enteric viruses in municipal well water, focusing on sanitary sewer sources. A total of 33 wells from 14 communities were sampled once for wastewater tracers and viruses. Wastewater tracers were detected in four of these wells, and five wells were virus- positive by qRT-PCR. These results, along with exclusion of wells with surface water sources, were used to select three wells for additional investigation. Viruses and wastewater tracers were found in the groundwater at all sites. Some wastewater tracers, such as ionic detergents, flame retardants, and cholesterol, were considered unambiguous evidence of wastewater. Sampling at any given time may not show concurrent virus and tracer presence; however, given sufficient sampling over time, a relation between wastewater tracers and virus occurrence was identified. Presence of infectious viruses at the wellhead demonstrates that high-capacity pumping induced sufficiently short travel times for the transport of infectious viruses. Therefore, drinking-water wells are vulnerable to contaminants that travel along fast groundwater flowpaths even if they contribute a small amount of virus-laden water to the well. These results suggest that vulnerability assessments require characterization of "low yield-fast transport" in addition to traditional "high yield-slow transport", pathways. ?? 2010 American Chemical Society.

Environmental Science & Technology↗

USGS2010

Regional groundwater-flow model of the Lake Michigan Basin in support of Great Lakes Basin water availability and use studies

A regional groundwater-flow model of the Lake Michigan Basin and surrounding areas has been developed in support of the Great Lakes Basin Pilot project under the U.S. Geological Survey's National Water Availability and Use Program. The transient 2-million-cell model incorporates multiple aquifers and pumping centers that create water-level drawdown that extends into deep saline waters. The 20-layer model simulates the exchange between a dense surface-water network and heterogeneous glacial deposits overlying stratified bedrock of the Wisconsin/Kankakee Arches and Michigan Basin in the Lower and Upper Peninsulas of Michigan; eastern Wisconsin; northern Indiana; and northeastern Illinois. The model is used to quantify changes in the groundwater system in response to pumping and variations in recharge from 1864 to 2005. Model results quantify the sources of water to major pumping centers, illustrate the dynamics of the groundwater system, and yield measures of water availability useful for water-resources management in the region. This report is a complete description of the methods and datasets used to develop the regional model, the underlying conceptual model, and model inputs, including specified values of material properties and the assignment of external and internal boundary conditions. The report also documents the application of the SEAWAT-2000 program for variable-density flow; it details the approach, advanced methods, and results associated with calibration through nonlinear regression using the PEST program; presents the water-level, drawdown, and groundwater flows for various geographic subregions and aquifer systems; and provides analyses of the effects of pumping from shallow and deep wells on sources of water to wells, the migration of groundwater divides, and direct and indirect groundwater discharge to Lake Michigan. The report considers the role of unconfined conditions at the regional scale as well as the influence of salinity on groundwater flow. Lastly, it describes several categories of limitations and discusses ways of extending the regional model to address issues at the local scale. Results of the simulations portray a regional groundwater-flow system that, over time, has largely maintained its natural predevelopment configuration but that locally has been strongly affected by well withdrawals. The quantity of rainfall in the Lake Michigan Basin and adjacent areas supports a dense surface-water network and recharge rates consistent with generally shallow water tables and predominantly shallow groundwater flow. At the regional scale, pumping has not caused major modifications of the shallow flow system, but it has resulted in decreases in base flow to streams and in direct discharge to Lake Michigan (about 2 percent of the groundwater discharged and about 0.5 cubic foot per second per mile of shoreline). On the other hand, well withdrawals have caused major reversals in regional flow patterns around pumping centers in deep, confined aquifers - most noticeably in the Cambrian-Ordovician aquifer system on the west side of Lake Michigan near the cities of Green Bay and Milwaukee in eastern Wisconsin, and around Chicago in northeastern Illinois, as well as in some shallow bedrock aquifers (for example, in the Marshall aquifer near Lansing, Mich.). The reversals in flow have been accompanied by large drawdowns with consequent local decrease in storage. On the west side of Lake Michigan, groundwater withdrawals have caused appreciable migration of the deep groundwater divides. Before the advent of pumping, the deep Lake Michigan groundwater-basin boundaries extended west of the Lake Michigan surface-water basin boundary, in some places by tens of miles. Over time, the pumping centers have replaced Lake Michigan as the regional sink for the deep flow system. The regional model is intended to support the framework pilot study of water availability and use for the Great Lakes Basin (Reeves, in press).

Illinois, Indiana, Michigan, Wisconsin↗

USGS2010

SWB: A modified Thornthwaite-Mather Soil-Water-Balance code for estimating groundwater recharge

A Soil-Water-Balance (SWB) computer code has been developed to calculate spatial and temporal variations in groundwater recharge. The SWB model calculates recharge by use of commonly available geographic information system (GIS) data layers in combination with tabular climatological data. The code is based on a modified Thornthwaite-Mather soil-water-balance approach, with components of the soil-water balance calculated at a daily timestep. Recharge calculations are made on a rectangular grid of computational elements that may be easily imported into a regional groundwater-flow model. Recharge estimates calculated by the code may be output as daily, monthly, or annual values.

Illinois, Michigan, Wisconsin↗

USGS2009

On constraining pilot point calibration with regularization in PEST

Ground water model calibration has made great advances in recent years with practical tools such as PEST being instrumental for making the latest techniques available to practitioners. As models and calibration tools get more sophisticated, however, the power of these tools can be misapplied, resulting in poor parameter estimates and/or nonoptimally calibrated models that do not suit their intended purpose. Here, we focus on an increasingly common technique for calibrating highly parameterized numerical models - pilot point parameterization with Tikhonov regularization. Pilot points are a popular method for spatially parameterizing complex hydrogeologic systems; however, additional flexibility offered by pilot points can become problematic if not constrained by Tikhonov regularization. The objective of this work is to explain and illustrate the specific roles played by control variables in the PEST software for Tikhonov regularization applied to pilot points. A recent study encountered difficulties implementing this approach, but through examination of that analysis, insight into underlying sources of potential misapplication can be gained and some guidelines for overcoming them developed. ?? 2009 National Ground Water Association.

Ground Water↗

USGS2008

Relating groundwater to seasonal wetlands in southeastern Wisconsin, USA

Historically, drier types of wetlands have been difficult to characterize and are not well researched. Nonetheless, they are considered to reflect the precipitation history with little, if any, regard for possible relation to groundwater. Two seasonal coastal wetland types (wet prairie, sedge meadow) were investigated during three growing seasons at three sites in the Lake Michigan Basin, Wisconsin, USA. The six seasonal wetlands were characterized using standard soil and vegetation techniques and groundwater measurements from the shallow and deep systems. They all met wetland hydrology criteria (e.g., water within 30 cm of land surface for 5% of the growing season) during the early portion of the growing season despite the lack of appreciable regional groundwater discharge into the wetland root zones. Although root-zone duration analyses did not fit a lognormal distribution previously noted in groundwater-dominated wetlands, they were able to discriminate between the plant communities and showed that wet prairie communities had shorter durations of continuous soil saturation than sedge meadow communities. These results demonstrate that the relative rates of groundwater outflows can be important for wetland hydrology and resulting wetland type. Thus, regional stresses to the shallow groundwater system such as pumping or low Great Lake levels can be expected to affect even drier wetland types. ?? Springer-Verlag 2008.

Hydrogeology Journal↗

USGS2008

Importance of unsaturated zone flow for simulating recharge in a humid climate

Transient recharge to the water table is often not well understood or quantified. Two approaches for simulating transient recharge in a ground water flow model were investigated using the Trout Lake watershed in north-central Wisconsin: (1) a traditional approach of adding recharge directly to the water table and (2) routing the same volume of water through an unsaturated zone column to the water table. Areas with thin (less than 1 m) unsaturated zones showed little difference in timing of recharge between the two approaches; when water was routed through the unsaturated zone, however, less recharge was delivered to the water table and more discharge occurred to the surface because recharge direction and magnitude changed when the water table rose to the land surface. Areas with a thick (15 to 26 m) unsaturated zone were characterized by multimonth lags between infiltration and recharge, and, in some cases, wetting fronts from precipitation events during the fall overtook and mixed with infiltration from the previous spring snowmelt. Thus, in thicker unsaturated zones, the volume of water infiltrated was properly simulated using the traditional approach, but the timing was different from simulations that included unsaturated zone flow. Routing of rejected recharge and ground water discharge at land surface to surface water features also provided a better simulation of the observed flow regime in a stream at the basin outlet. These results demonstrate that consideration of flow through the unsaturated zone may be important when simulating transient ground water flow in humid climates with shallow water tables.

Ground Water↗

USGS2008

Comparison of total mercury and methylmercury cycling at five sites using the small watershed approach

The small watershed approach is well-suited but underutilized in mercury research. We applied the small watershed approach to investigate total mercury (THg) and methylmercury (MeHg) dynamics in streamwater at the five diverse forested headwater catchments of the US Geological Survey Water, Energy, and Biogeochemical Budgets (WEBB) program. At all sites, baseflow THg was generally less than 1 ng L −1 and MeHg was less than 0.2 ng L −1 . THg and MeHg concentrations increased with streamflow, so export was primarily episodic. At three sites, THg and MeHg concentration and export were dominated by the particulate fraction in association with POC at high flows, with maximum THg (MeHg) concentrations of 94 (2.56) ng L −1 at Sleepers River, Vermont; 112 (0.75) ng L −1 at Rio Icacos, Puerto Rico; and 55 (0.80) ng L −1 at Panola Mt., Georgia. Filtered (<0.7 μm) THg increased more modestly with flow in association with the hydrophobic acid fraction (HPOA) of DOC, with maximum filtered THg concentrations near 5 ng L −1 at both Sleepers and Icacos. At Andrews Creek, Colorado, THg export was also episodic but was dominated by filtered THg, as POC concentrations were low. MeHg typically tracked THg so that each site had a fairly constant MeHg/THg ratio, which ranged from near zero at Andrews to 15% at the low-relief, groundwater-dominated Allequash Creek, Wisconsin. Allequash was the only site with filtered MeHg consistently above detection, and the filtered fraction dominated both THg and MeHg. Relative to inputs in wet deposition, watershed retention of THg (minus any subsequent volatilization) was 96.6% at Allequash, 60% at Sleepers, and 83% at Andrews. Icacos had a net export of THg, possibly due to historic gold mining or frequent disturbance from landslides. Quantification and interpretation of Hg dynamics was facilitated by the small watershed approach with emphasis on event sampling.

Environmental Pollution↗

USGS2008

Effects of climate and land management change on streamflow in the driftless area of Wisconsin

Baseflow and precipitation in the Kickapoo River Watershed, located in the Driftless Area of Wisconsin, exhibit a step increase around 1970, similar to minimum and median flows in many other central and eastern USA streams. Potential effects on streamflow due to climatic and land management changes were evaluated by comparing volumetric changes in the hydrologic budget before and after 1970. Increases in precipitation do not fully account for the increase in baseflow, which appears to be offset by a volumetric decrease in stormflow. This suggests that factors that influence the partitioning of precipitation into overland runoff or infiltration have changed. A transition from relatively more intensive to relatively less intensive agricultural land use is generally associated with higher infiltration rates, and likely influences partitioning of flow. Changes in agricultural land management practices in the Driftless Area, which began in the mid-1930s, do not coincide with the abrupt increase in baseflow around 1970. Instead, the timing of hydrologic change appears to coincide with changes in precipitation, whereas the magnitude of the change in baseflow and stormflow was likely amplified by changes in agricultural land management. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

USGS2006

Measuring groundwater-surface water interaction and its effect on wetland stream benthic productivity, Trout Lake watershed, northern Wisconsin, USA

Measurements of groundwater-surface water exchange at three wetland stream sites were related to patterns in benthic productivity as part of the US Geological Survey's Northern Temperate Lakes-Water, Energy and Biogeochemical Budgets (NTL-WEBB) project. The three sites included one high groundwater discharge (HGD) site, one weak groundwater discharge (WGD) site, and one groundwater recharge (GR) site. Large upward vertical gradients at the HGD site were associated with smallest variation in head below the stream and fewest gradient reversals between the stream and the groundwater beneath the stream, and the stream and the adjacent streambank. The WGD site had the highest number of gradient reversals reflecting the average condition being closest to zero vertical gradient. The duration of groundwater discharge events was related to the amount of discharge, where the HGD site had the longest strong-gradient durations for both horizontal and vertical groundwater flow. Strong groundwater discharge also controlled transient temperature and chemical hyporheic conditions by limiting the infiltration of surface water. Groundwater-surface water interactions were related to highly significant patterns in benthic invertebrate abundance, taxonomic richness, and periphyton respiration. The HGD site abundance was 35% greater than in the WGD site and 53% greater than the GR site; richness and periphyton respiration were also significantly greater (p???0.001, 31 and 44%, respectively) in the HGD site than in the GR site. The WGD site had greater abundance (27%), richness (19%) and periphyton respiration (39%) than the GR site. This work suggests groundwater-surface water interactions can strongly influence benthic productivity, thus emphasizing the importance of quantitative hydrology for management of wetland-stream ecosystems in the northern temperate regions. ?? 2005 Elsevier B.V. All rights reserved.

Conference Paper↗

USGS2006

The importance of diverse data types to calibrate a watershed model of the Trout Lake Basin, Northern Wisconsin, USA

As part of the USGS Water, Energy, and Biogeochemical Budgets project and the NSF Long-Term Ecological Research work, a parameter estimation code was used to calibrate a deterministic groundwater flow model of the Trout Lake Basin in northern Wisconsin. Observations included traditional calibration targets (head, lake stage, and baseflow observations) as well as unconventional targets such as groundwater flows to and from lakes, depth of a lake water plume, and time of travel. The unconventional data types were important for parameter estimation convergence and allowed the development of a more detailed parameterization capable of resolving model objectives with well-constrained parameter values. Independent estimates of groundwater inflow to lakes were most important for constraining lakebed leakance and the depth of the lake water plume was important for determining hydraulic conductivity and conceptual aquifer layering. The most important target overall, however, was a conventional regional baseflow target that led to correct distribution of flow between sub-basins and the regional system during model calibration. The use of an automated parameter estimation code: (1) facilitated the calibration process by providing a quantitative assessment of the model's ability to match disparate observed data types; and (2) allowed assessment of the influence of observed targets on the calibration process. The model calibration required the use of a 'universal' parameter estimation code in order to include all types of observations in the objective function. The methods described in this paper help address issues of watershed complexity and non-uniqueness common to deterministic watershed models. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

USGS2005

Using stable isotopes of water and strontium to investigate the hydrology of a natural and a constructed wetland

Wetlands cannot exist without water, but wetland hydrology is difficult to characterize. As a result, compensatory wetland mitigation often only assumes the proper hydrology has been created. In this study, water sources and mass transfer processes in a natural and constructed wetland complex were investigated using isotopes of water and strontium. Water isotope profiles in the saturated zone revealed that the natural wetland and one site in the constructed wetland were primarily fed by ground water; profiles in another constructed wetland site showed recent rain was the predominant source of water in the root zone. Water isotopes in the capillary fringe indicated that the residence time for rain is less in the natural wetland than in the constructed wetland, thus transpiration (an important water sink) was greater in the natural wetland. Strontium isotopes showed a systematic difference between the natural and constructed wetlands that we attribute to the presence or absence of peat. In the peat‐rich natural wetland, δ 87 Sr in the pore water increased along the flowline due to preferential weathering of minerals containing radiogenic Sr in response to elevated Fe concentrations in the water. In the constructed wetland, where peat thickness was thin and Fe concentrations in water were negligible, δ 87 Sr did not increase along the flowline. The source of the peat (on‐site or off‐site derived) applied in the constructed wetland controlled the δ 87 Sr at the top of the profile, but the effects were restricted by strong cation exchange in the underlying fluvial sediments. Based on the results of this study, neither constructed wetland site duplicated the water source and weathering environment of the adjoining natural wetland. Moreover, stable isotopes were shown to be effective tools for investigating wetlands and gaining insight not easily obtained using non‐isotopic techniques. These tools have potential widespread application to wetlands that have distinct isotopic endmember sources.

Ground Water↗

USGS2005

Delineating a recharge area for a spring using numerical modeling, Monte Carlo techniques, and geochemical investigation

Recharge areas of spring systems can be hard to identify, but they can be critically important for protection of a spring resource. A recharge area for a spring complex in southern Wisconsin was delineated using a variety of complementary techniques. A telescopic mesh refinement (TMR) model was constructed from an existing regional-scale ground water flow model. This TMR model was formally optimized using parameter estimation techniques; the optimized "best fit" to measured heads and fluxes was obtained by using a horizontal hydraulic conductivity 200% larger than the original regional model for the upper bedrock aquifer and 80% smaller for the lower bedrock aquifer. The uncertainty in hydraulic conductivity was formally considered using a stochastic Monte Carlo approach. Two-hundred model runs used uniformly distributed, randomly sampled, horizontal hydraulic conductivity values within the range given by the TMR optimized values and the previously constructed regional model. A probability distribution of particles captured by the spring, or a "probabilistic capture zone," was calculated from the realistic Monte Carlo results (136 runs of 200). In addition to portions of the local surface watershed, the capture zone encompassed areas outside of the watershed - demonstrating that the ground watershed and surface watershed do not coincide. Analysis of water collected from the site identified relatively large contrasts in chemistry, even for springs within 15 m of one another. The differences showed a distinct gradation from Ordovician-carbonate-dominated water in western spring vents to Cambrian-sandstone-influenced water in eastern spring vents. The difference in chemistry was attributed to distinctive bedrock geology as demonstrated by overlaying the capture zone derived from numerical modeling over a bedrock geology map for the area. This finding gives additional confidence to the capture zone calculated by modeling.

Ground Water↗

USGS2005

Stepwise use of GFLOW and MODFLOW to determine relative importance of shallow and deep receptors

A stepwise modeling approach is implemented in which a regional one-layer analytic element model is used to simulate the flow system and to furnish boundary conditions for an extracted local three-dimensional model. In this case study the stepwise approach is used to evaluate the fate of recharge in the Menomonee Valley adjacent to Lake Michigan. Two major receptors exist for recharge that flows through contaminated valley fill: the surface water estuary and a tunnel system constructed ???75 to 110 m below land surface to store storm runoff. The primary objective of the modeling is to delineate the contributing areas of recharge to each receptor. Of interest is the ability of the one-layer regional model to furnish flux boundary conditions to the local three-dimensional model despite the presence of vertical flow conditions at the boundaries of the local model. Sensitivity analysis suggests that the local model was insensitive to the vertical distribution of the flux. Each step of the modeling approach demonstrates that both receptors play an important role in capturing valley recharge. The pattern of capture of the one-layer model differed in shape from that delineated by the multi-layer local model in the presence of a flow system with pronounced vertical anisotropy and with sinks drawing water from different elevations.

Ground Water↗

USGS2005

Simulating ground water-lake interactions: Approaches and insights

Approaches for modeling lake-ground water interactions have evolved significantly from early simulations that used fixed lake stages specified as constant head to sophisticated LAK packages for MODFLOW. Although model input can be complex, the LAK package capabilities and output are superior to methods that rely on a fixed lake stage and compare well to other simple methods where lake stage can be calculated. Regardless of the approach, guidelines presented here for model grid size, location of three-dimensional flow, and extent of vertical capture can facilitate the construction of appropriately detailed models that simulate important lake-ground water interactions without adding unnecessary complexity. In addition to MODFLOW approaches, lake simulation has been formulated in terms of analytic elements. The analytic element lake package had acceptable agreement with a published LAK1 problem, even though there were differences in the total lake conductance and number of layers used in the two models. The grid size used in the original LAK1 problem, however, violated a grid size guideline presented in this paper. Grid sensitivity analyses demonstrated that an appreciable discrepancy in the distribution of stream and lake flux was related to the large grid size used in the original LAK1 problem. This artifact is expected regardless of MODFLOW LAK package used. When the grid size was reduced, a finite-difference formulation approached the analytic element results. These insights and guidelines can help ensure that the proper lake simulation tool is being selected and applied.

Ground Water↗

USGS2005

Scale effects of hydrostratigraphy and recharge zonation on base flow

Uncertainty regarding spatial variations of model parameters often results in the simplifying assumption that parameters are spatially uniform. However, spatial variability may be important in resource assessment and model calibration. In this paper, a methodology is presented for estimating a critical basin size, above which base flows appear to be relatively less sensitive to the spatial distribution of recharge and hydraulic conductivity, and below which base flows are relatively more sensitive to this spatial variability. Application of the method is illustrated for a watershed that exhibits distinct infiltration patterns and hydrostratigraphic layering. A ground water flow model (MODFLOW) and a parameter estimation code (UCODE) were used to evaluate the influence of recharge zonation and hydrostratigraphic layering on base flow distribution. Optimization after removing spatial recharge variability from the calibrated model altered base flow simulations up to 53% in watersheds smaller than 40 km2. Merging six hydrostratigraphic units into one unit with average properties increased base flow residuals up to 83% in basins smaller than 50 km2. Base flow residuals changed <5% in watersheds larger than 40 and 50 km2 when recharge and hydrostratigraphy were simplified, respectively; thus, the critical basin size for the example area is ???40 to 50 km2. Once identified for an area, a critical basin size could be used to guide the scale of future investigations. By ensuring that parameter discretization needed to capture base flow distribution is commensurate with the scope of the investigation, uncertainty caused by overextending uniform parameterization or by estimating extra parameter values is reduced. ?? 2006 National Ground Water Association.

Ground Water↗

USGS2005

Estimating recharge rates with analytic element models and parameter estimation

Quantifying the spatial and temporal distribution of recharge is usually a prerequisite for effective ground water flow modeling. In this study, an analytic element (AE) code (GFLOW) was used with a nonlinear parameter estimation code (UCODE) to quantify the spatial and temporal distribution of recharge using measured base flows as calibration targets. The ease and flexibility of AE model construction and evaluation make this approach well suited for recharge estimation. An AE flow model of an undeveloped watershed in northern Wisconsin was optimized to match median annual base flows at four stream gages for 1996 to 2000 to demonstrate the approach. Initial optimizations that assumed a constant distributed recharge rate provided good matches (within 5%) to most of the annual base flow estimates, but discrepancies of >12% at certain gages suggested that a single value of recharge for the entire watershed is inappropriate. Subsequent optimizations that allowed for spatially distributed recharge zones based on the distribution of vegetation types improved the fit and confirmed that vegetation can influence spatial recharge variability in this watershed. Temporally, the annual recharge values varied >2.5-fold between 1996 and 2000 during which there was an observed 1.7-fold difference in annual precipitation, underscoring the influence of nonclimatic factors on interannual recharge variability for regional flow modeling. The final recharge values compared favorably with more labor-intensive field measurements of recharge and results from studies, supporting the utility of using linked AE-parameter estimation codes for recharge estimation. Copyright ?? 2005 The Author(s).

Ground Water↗