Hydrogeology and groundwater geochemistry of till confining units and confined aquifers in glacial deposits near Litchfield, Cromwell, Akeley, and Olivia, Minnesota, 2014–18

Jared J. Trost; Anna-Turi Maher; William W. Simpkins; Alyssa N. Witt; James R. Stark; Justin Blum; Andrew M. Berg

doi:10.3133/sir20205127

Home/Archive/USGS

USGS · sir20205127

Hydrogeology and groundwater geochemistry of till confining units and confined aquifers in glacial deposits near Litchfield, Cromwell, Akeley, and Olivia, Minnesota, 2014–18

Jared J. Trost·Anna-Turi Maher·William W. Simpkins·Alyssa N. Witt·James R. Stark·Justin Blum·Andrew M. Berg

View original source ↗Read source document ↗Export metadata ↓

Abstract

Confined (or buried) aquifers of glacial origin overlain by till confining units provide drinking water to hundreds of thousands of Minnesota residents. The sustainability of these groundwater resources is not well understood because hydraulic properties of till that control vertical groundwater fluxes (leakage) to underlying aquifers are largely unknown. The U.S. Geological Survey, Iowa State University, Minnesota Geological Survey, and Minnesota Department of Health investigated hydraulic properties and groundwater flow through till confining units using field studies and heuristic MODFLOW simulations. Till confining units in the following late-Wisconsinan stratigraphic units (with locations in parentheses) were characterized: Des Moines lobe till of the New Ulm Formation (Litchfield, Minnesota), Superior lobe till of the Cromwell and Aitkin Formations (Cromwell, Minn.), and Wadena lobe till of the Hewitt Formation (hydrogeology field camp [HFC] near Akeley, Minn.). Pre-Illinoian till of the Good Thunder formation (Olivia, Minn.) was also characterized. Hydraulic and geochemical field data were collected from sediment cores and a series of five piezometer nests. Each nest consisted of five to eight piezometers screened at short vertical intervals in hydrostratigraphic units including (if present) surficial aquifers, till confining units, confined/buried aquifers, and underlying bedrock. Till hydraulic conductivity was estimated from slug tests (horizontal [ K h ]) and constant-rate aquifer tests in the confined aquifer (vertical [ K v ]). Travel times through the till were evaluated with Darcy’s law and stable isotope concentrations. A series of heuristic MODFLOW simulations were used to evaluate groundwater fluxes through till across the range of till hydraulic properties and pumping rates observed at the field sites. The field data demonstrated variability in hydraulic properties between and within till stratigraphic units horizontally and vertically. The variability in hydraulic properties within and between sites resulted in substantial differences in groundwater flux through till. A conceptual understanding that emerges from the vertical till profiles is that they are not homogeneous hydrostratigraphic units with uniform properties; rather, each vertical sequence is a heterogeneous mixture of glacial sediment with differing abilities to transmit water. Till thicknesses varied from 60 to 166 feet, and till textures ranged from a sandy loam (Hewitt Formation, HFC site) to a silt loam/clay loam (Good Thunder formation, Olivia site). Till Kh varied by one to three orders of magnitude within each piezometer nest. Four piezometer nests had downward hydraulic gradients ranging from 0.04 to 0.56, and one nest had a slight upward hydraulic gradient of 0.02. The Cromwell, HFC, and Litchfield 1 sites were examples of “leaky” tills with high Kv (0.001 to 1.1 feet per day [ft/d]) and geometric mean Kh (0.03 to 0.07 ft/d) and extensive vertical hydraulic connectivity between the confined aquifer and the overlying till. Estimated groundwater travel times through these sites ranged from 1 to 81 years, and two of these sites had tritium throughout their till profiles. The tills at the other two sites, Olivia and Litchfield 2, were effective confining units that had low Kv (0.001 to 0.0005 ft/d) and geometric mean Kh (0.0002 to 0.004 ft/d). The till piezometers at these sites had no drawdown response to short-term (up to 10 hours for Olivia and up to 5 days for Litchfield) high-capacity pumping from the confined aquifer. Estimated groundwater travel times through the tills at these sites ranged from 165 to nearly 1,800 years, and tritium was only detected in the upper one-third of these till profiles. Across all sites, the till vertical anisotropy (ratio of K h to K v ) ranged by four orders of magnitude from 0.05 at the Cromwell nest to 70 at the Litchfield 1 nest. Stable isotopes of oxygen and hydrogen indicate that groundwater throughout all five till profiles is younger than the last glacial advance into Minnesota at about 11,000 years ago. The heuristic modeling demonstrated that, for understanding sustainability of groundwater pumping from confined aquifers, knowledge of till hydraulic properties is just as important as knowledge of aquifer hydraulic properties. Substantial differences in groundwater fluxes into and through till were observed across hydrogeologic settings representative of the field sites. Over long periods of time (hundreds of years), pumping-induced hydraulic gradients are established in confined aquifer systems and, even in low hydraulic conductivity tills, these pumping-induced hydraulic gradients increase leakage into and through till compared to ambient conditions. In conclusion, groundwater flowing vertically downward through till confining units (leakage) replenishes water pumped from confined aquifers. Till hydraulic properties, such as those presented in this report, provide important information that can be used to quantify leakage rates through till. Till hydraulic properties are variable over short distances and profoundly affect leakage rates, demonstrating the importance of site-specific till hydraulic data for evaluating the sustainability of groundwater withdrawals from confined aquifers.

Explore related subjects

Scientific Investigations Report Report

Source-reported bounding extent: 44.74185630294231° to 47.01678007415223° latitude; -95.02418518066406° to -92.84614562988281° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Jared J. Trost, Anna-Turi Maher, William W. Simpkins, Alyssa N. Witt, James R. Stark, Justin Blum, Andrew M. Berg. 2020-12-22. Hydrogeology and groundwater geochemistry of till confining units and confined aquifers in glacial deposits near Litchfield, Cromwell, Akeley, and Olivia, Minnesota, 2014–18. https://doi.org/10.3133/sir20205127

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Analysis of groundwater geochemistry in the Lansing, Michigan, area

Groundwater resources in the Lansing, Michigan, area are protected by a wellhead protection plan, but previous studies have indicated that groundwater in this area may be susceptible to surface contamination. The U.S. Geological Survey, in cooperation with the Lansing Board of Water & Light, analyzed geochemical and aquifer characteristics to assess the potential for contamination of groundwater resources. Recharge age categories were determined from tritium concentrations for some wells, and by using ternary diagrams and principal component analysis for other wells. These age classifications, combined with hydrogeologic and geochemical data, were analyzed by using statistical methods to identify trends that may indicate areas of increased groundwater susceptibility. Wells with modern recharge tended to have higher sodium, chloride, and chloride:bromide ratios than did wells with mixed recharge, suggesting greater surface connectivity and potential influence from anthropogenic sources such as road salt. In contrast, wells in the mixed age category exhibited more stable geochemistry and lower constituent concentrations. Classification and regression tree models identified aquifer geometry—particularly the thickness of the uppermost bedrock unit and presence of confining layers—as key predictors of recharge age and geochemical indicators. The trends identified from these analyses could help to reveal locations in the Lansing area that are among the most susceptible to contamination.

Scientific Investigations Report↗

USGS2026

Evaluating the number of realizations (train-test splits) necessary to ensure robust estimates of conventional hydrothermal resource favorability using Monte Carlo cross-validation, the Great Basin, Western USA

Supervised machine learning approaches used to produce favorability maps in resource assessments commonly partition an input dataset into training and testing subsets (that is, through train-test splits). However, if a single split is not representative of underlying resource conditions, the resulting favorability map may be biased, particularly with sparse or heterogeneous datasets. To mitigate the potential for unrepresentative favorability maps, recent works by U.S. Geological Survey researchers have used Monte Carlo cross-validation to implement a large number ( n) of randomized train-test splits to produce n realizations. The realizations are then aggregated using cell-wise median predictions to produce the resulting final favorability map, recognizing that the median will remove extreme values that might have been estimated from a comparatively few unrepresentative realizations. Theoretically, the median prediction at every cell converges to a single estimate as n approaches infinity, and the question addressed herein is: How many realizations are necessary to get a robust estimate (that is, an estimate unaffected by additional realizations)? Therefore, convergence metrics are described to identify the number of train-test splits and realizations needed to make a reliable map using the data and supervised machine learning workflow for conventional hydrothermal resource assessment of the Great Basin of the Western United States. For this study, robust prediction occurs with n approximately equal to 600 realizations, and smaller ensembles (for example, n is approximately equal to 200 or even 100) may be suitable during preliminary or exploratory analyses. Although the resulting threshold numbers will vary depending on the data and modeling approach used, the convergence metrics proposed and applied herein are applicable to other supervised machine learning studies that use Monte Carlo cross-validation.

California, Idaho, Nevada, Oregon, Utah↗

USGS2026

Groundwater budget for the Portneuf River Basin, southeastern Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources, developed the first groundwater budget for the upper Portneuf River Basin in southeastern Idaho, a tributary of the Snake River. A groundwater budget for the Portneuf River Basin supports water resource and administrative management decisions within the basin and for downstream users of the eastern Snake River Plain aquifer. The analysis quantifies inflows, outflows, and residuals for an aquifer system defined by the Idaho Department of Water Resources under three recent precipitation conditions: an average year (2016), a dry year (2021), and a wet year (2023). The budget was calculated for the basin and for two subbasins delineated by the U.S. Geological Survey streamgage Portneuf River at Topaz, Idaho (13073000): the above Topaz subbasin, which includes the Portneuf-Gem valley, and the below Topaz subbasin, which includes Marsh Valley. Aquifer inflows include mountain-front recharge; areal recharge; recharge from applied water, including recharge from irrigation, canal seepage, and effluent returns; groundwater flow into the basin; and reservoir seepage. Outflows include groundwater pumping for agricultural, municipal, and domestic uses; subsurface outflow through the Portneuf gap; and discharge to streams as base flow. Estimated annual inflows for the full basin ranged from 116,880 acre-feet in 2021 to 154,880 acre-feet in 2023, and outflows ranged from 132,190 acre-feet in 2016 to 158,490 acre-feet in 2023. Basinwide residuals, the change in storage plus uncertainty, were −2,410 acre-feet in 2016, −18,410 acre-feet in 2021, and −3,610 acre-feet in 2023, and each subbasin had a mixture of negative and positive residuals. These basinwide results, interpreted as change in storage plus water-budget uncertainty, indicate minor net losses that increase with dry conditions, consistent with well depth-to-water measurements indicating gradual multi-decadal declines. Patterns in subbasin residuals are driven by irrigation practices and surface-water management; for example, the above Topaz subbasin relies more on groundwater for irrigation, resulting in more negative residuals during dry years. In contrast, the below Topaz subbasin depends more on surface water, producing positive residuals during dry years because of irrigation losses and canal seepage. The analyses for both subbasins indicate decades-long groundwater declines, which agree with the overall negative residuals calculated. The largest sources of uncertainty include climate inputs (precipitation and evapotranspiration) and poorly constrained components such as groundwater inflow and tributary stream infiltration. These findings provide insight into how climatic variables and irrigation practices affect groundwater inflows and outflows within the Portneuf River Basin.

Idaho↗