USGS Science⌕ Search

SEARCH · USGS Science

Results for “Agricultural Water Management”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

754 records · Page 31Linked to original sources
USGS2025

Using probability difference to compare streamflow information of alternatives for efficient operation of monitoring networks

Efficient operation of streamflow monitoring networks requires investments in technology and labor that provide the greatest benefits from available resources. Economic analyses comparing the costs and benefits from different types of alternatives for monitoring have not been practical to implement. Streamflow information provides a generic measure of benefits that can be incorporated into operational decisions as an objective for monitoring networks. A methodology for comparing how accuracy, monitoring period, and monitoring instead of modeling affects streamflow information is developed from information-theoretic approaches for network design but contributes three novel features: (1) a probability-difference model for conditional probability of monotonically paired variables, (2) explicit discounting of unverified information that may exceed the accuracy of streamflow records, and (3) run analysis to account for non-stationarity in streamflow probabilities. Application of the methodology to the U.S. Geological Survey streamflow monitoring network indicates the value of monitoring period to reduce the uncertainty of streamflow probabilities and, thus, increase streamflow information. The methodology has important limitations, particularly for sites with non-perennial streamflow, but demonstrates that probability difference could be used to evaluate operational alternatives to increase the efficiency of monitoring networks.

PLOS Water↗
USGS2025

Natural source zone depletion of crude oil in the subsurface: Processes controlling mass losses of individual compounds

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o -xylene, and n -alkanes are more depleted compared to benzene, ethylbenzene, and m - and p -xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Minnesota↗
USGS1986

Regionalization of low-flow characteristics of Tennessee streams

Procedures for estimating 3-day 2-year, 3-day 10-year, 3-day 20-year, and 7-day 10-year low flows at ungaged stream sites in Tennessee are based on surface geology and drainage area size. One set of equations applies to west Tennessee streams, and another set applies to central and east Tennessee streams. The equations do not apply to streams where flow is significantly altered by activities of man. Standard errors of estimate of equations for west Tennessee are 24 to 32% and for central and east Tennessee 31 to 35%. Streamflow recession indexes, in days/log cycle, are used to account for effects of geology of the drainage basin on low flow of streams. The indexes in Tennessee range from 32 days/log cycle for clay and shale to 350 days/log cycle for gravel and sand, indicating different aquifer characteristics of the geologic units that sustain streamflows during periods of no surface runoff. Streamflow recession rate depends primarily on transmissivity and storage characteristics of the aquifers, and the average distance from stream channels to basin divides. Geology and drainage basin size are the most significant variables affecting low flow in Tennessee streams according to regression analyses. (Author 's abstract)

Tennessee↗
USGS1998

Methods and guidelines for effective model calibration; with application to UCODE, a computer code for universal inverse modeling, and MODFLOWP, a computer code for inverse modeling with MODFLOW

This report documents methods and guidelines for model calibration using inverse modeling. The inverse modeling and statistical methods discussed are broadly applicable, but are presented as implemented in the computer programs UCODE, a universal inverse code that can be used with any application model, and MODFLOWP, an inverse code limited to one application model. UCODE and MODFLOWP perform inverse modeling, posed as a parameter-estimation problem, by calculating parameter values that minimize a weighted least-squares objective function using nonlinear regression. Minimization is accomplished using a modified Gauss-Newton method, and prior, or direct, information on estimated parameters can be included in the regression. Inverse modeling in many fields is plagued by problems of instability and nonuniqueness, and obtaining useful results depends on (1) defining a tractable inverse problem using simplifications appropriate to the system under investigation and (2) wise use of statistics generated using calculated sensitivities and the match between observed and simulated values, and associated graphical analyses. Fourteen guidelines presented in this work suggest ways of constructing and calibrating models of complex systems such that the resulting model is as accurate and useful as possible.

Water-Resources Investigations Report↗
USGS2024

Characterization of groundwater resources near the southeastern part of Puget Sound, Washington

More than 1 million people live within the active model area (AMA) in the southeastern part of the lowlands surrounding Puget Sound, or Puget Lowland, Washington, and groundwater is the source for approximately one-half of their public, domestic, and irrigation water demands. The 887-square-mile AMA, located in King and Pierce Counties, represents the area of analysis for the conceptual hydrogeologic framework and numerical groundwater-flow models within the study area and includes the Puyallup River and Chambers-Clover Creek watersheds. To assess the potential hydrologic and anthropogenic impacts to groundwater and the connected surface-water resources, conceptual and numerical groundwater-flow models of groundwater flow were developed by the U.S. Geological Survey Washington Water Science Center in close cooperation with 18 water-resource agencies and stakeholders. This multichapter volume documents the development of the conceptual and numerical groundwater-flow models of groundwater flow. Chapters A, B, and C provide an overall introduction to the multichapter volume (Chapter A), the conceptual hydrogeologic framework (Chapter B), and the groundwater budget (Chapter C). Chapters D and E describe numerical groundwater-flow model construction and calibration (Chapter D) and the numerical groundwater-flow model results (Chapter E). Collectively, these reports present a characterization and simulation tool for groundwater resources near the southeastern part of Puget Sound, Washington.

Scientific Investigations Report↗
USGS2000

Use of a ground-penetrating radar system to detect pre- and post-flood scour at selected bridge sites in New Hampshire, 1996-98

Ground-penetrating radar was used to measure the depth and extent of existing and infilled scour holes and previous scour surfaces at seven bridges in New Hampshire from April 1996 to November 1998. Ground-penetrating-radar survey techniques initially were used by the U.S. Geological Survey to study streambed scour at 30 bridges. Sixteen of the 30 bridges were re-surveyed where floods exceeded a 2-year recurrence interval. A 300-megahertz signal was used in the ground-penetrating radar system that penetrated through depths as great as 20 feet of water and as great as 32 feet of streambed materials. Existing scour-hole dimensions, infilled thickness, previous scour surfaces, and streambed materials were detected using ground-penetrating radar. Depths to riprap materials and pier footings were identified and verified with bridge plans. Post data-collection-processing techniques were applied to assist in the interpretation of the data, and the processed data were displayed and printed as line plots. Processing included distance normalization, migration, and filtering but processing was kept to a minimum and some interference from multiple reflections was left in the record. Of the 16 post-flood bridges, 22 ground-penetrating-radar cross sections at 7 bridges were compared and presented in this report. Existing scour holes were detected during 1996 (pre-flood) data collection in nine cross sections where scour depths ranged from 1 to 3 feet. New scour holes were detected during 1998 (post-flood) data collection in four cross sections where scour depths were as great as 4 feet deep. Infilled scour holes were detected in seven cross sections, where depths of infilling ranged from less than 1 to 4 feet. Depth of infilling by means of steel rod and hammer was difficult to verify in the field because of cobble and boulder streambeds or deep water. Previous scour surfaces in streambed materials were identified in 15 cross sections and the depths to these surfaces ranged from 1 to 10 feet below the streambed. Riprap materials or pier footings were identified in all cross sections. Calculated record depths generally agree with bridge plans. Pier footings were exposed at two bridges and steel pile was exposed at one bridge. Exposures were verified by field observations.

New Hampshire↗
USGS2025

Road salt collection and redistribution at an urban rain garden on sandy soil, Gary, Indiana

Rain gardens installed as green infrastructure to divert storm runoff from entering combined sewers also collect dissolved constituents and particulates. An urban rain garden in northwestern Indiana, USA, was continuously monitored from November 2019 to May 2021 to evaluate the fate of dissolved constituents entering the rain garden in runoff. Physical and chemical properties of soils in the rain garden were also monitored, along with underlying groundwater. Linear regression models relating specific conductance to chloride concentration indicated that the 0.0371-ha (3998 square feet) rain garden collected approximately 1490 kg (3285 pounds) of road salt from the surrounding 0.2228 ha (24,500 square feet) of impervious surfaces. Soils and groundwater were seasonally affected by road salt application but carryover from year to year was not indicated. Rain garden soil permeability (5.20 × 10 −5 to 9.72 × 10 −5 m/s) remained unchanged during the study period and soil organic carbon generally increased under native vegetation. The results suggest that a rain garden built on sandy soil can divert substantial quantities of runoff and dissolved constituents from combined sewers; however, chloride is transported to sub-infrastructure groundwater that eventually discharges to adjacent waterways with concentrations lower than those observed in runoff.

Indiana↗
USGS1996

Techniques to estimate generalized skew coefficients of annual peak streamflow for natural basins in Texas

This report presents two techniques to estimate generalized skew coefficients used for log-Pearson Type III peak-streamflow frequency analysis of natural basins in Texas. A natural basin has less than 10 percent impervious cover, and less than 10 percent of its drainage area is controlled by reservoirs. The estimation of generalized skew coefficients is based on annual peak and historical peak streamflow for all U.S. Geological Survey streamflow-gaging stations having at least 20 years of annual peak-streamflow record from natural basins in Texas. Station skew coefficients calculated for each of 255 Texas stations were used to estimate generalized skew coefficients for Texas. One technique to estimate generalized skew coefficients involved the use of regression equations developed for each of eight regions in Texas, and the other involved development of a statewide map of generalized skew coefficients. The weighted mean of the weighted mean standard errors of the regression equations for the eight regions is 0.36 log 10 skew units, and the weighted mean standard error of the map is 0.35 log 10 skew units. The technique based on the map is preferred for estimating generalized skew coefficients because of its smooth transition from one region of the State to another.

Texas↗
USGS2025

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report↗
USGS1998

Characterization of hydrogeologic units using matrix properties, Yucca Mountain, Nevada

Determination of the suitability of Yucca Mountain, in southern Nevada, as a geologic repository for high-level radioactive waste requires the use of numerical flow and transport models. Input for these models includes parameters that describe hydrologic properties and the initial and boundary conditions for all rock materials within the unsaturated zone, as well as some of the upper rocks in the saturated zone. There are 30 hydrogeologic units in the unsaturated zone, and each unit is defined by limited ranges where a discrete volume of rock contains similar hydrogeologic properties. These hydrogeologic units can be easily located in space by using three-dimensional lithostratigraphic models based on relation- ships of the properties with the lithostratigraphy. Physical properties of bulk density, porosity, and particle density; flow properties of saturated hydraulic conductivity and moisture-retention characteristics; and the state variables (variables describing the current state of field conditions) of saturation and water potential were determined for each unit. Units were defined using (1) a data base developed from 4,892 rock samples collected from the coring of 23 shallow and 8 deep boreholes, (2) described lithostratigraphic boundaries and corresponding relations to porosity, (3) recognition of transition zones with pronounced changes in properties over short vertical distances, (4) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (5) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. This study describes the correlation of hydrologic properties to porosity, a property that is well related to the lithostratigraphy and depositional and cooling history of the volcanic deposits and can, therefore, be modeled to be distributed laterally. Parameters of the hydrogeologic units developed in this study and the relation of flow properties to porosity that are described can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.

Nevada↗
USGS2025

Hydrogeologic investigation, framework, and conceptual flow model of the Antlers aquifer, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. Every 20 years, the Oklahoma Water Resources Board is required to update the hydrologic investigation on which the maximum annual yield determinations were based. The maximum annual yield allocated per acre of land is used to set the equal-proportionate share pumping rate. The maximum annual yield of 5,913,600 acre-feet per year and equal-proportionate-share of 2.1 acre-feet per acre per year currently (2025) in place for the Antlers aquifer were issued by the Oklahoma Water Resources Board on February 14, 1995. Because more than 20 years have elapsed since the 1995 final order for the Antlers aquifer was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, completed an in-depth hydrologic study that included a hydrogeologic framework and conceptual groundwater-flow model for the 1980–2022 study period. The results of an analysis of land use, long-term climate patterns, streamflow and base-flow patterns, historical groundwater use, as well as groundwater-level fluctuations across the Antlers aquifer are described. In addition, groundwater quality was analyzed for total dissolved solids concentrations and major ions for the Antlers aquifer. An updated hydrogeologic framework was developed that included refining the aquifer boundary in Oklahoma, the creation of new potentiometric surface and saturated thickness of fresh groundwater maps, one multiple-well aquifer test, slug tests, and an analysis of lithologic logs across the aquifer. A conceptual groundwater flow model and water budget were developed by incorporating estimates of recharge from precipitation, saturated-zone evapotranspiration, streambed seepage, lateral groundwater flows, vertical leakage, and withdrawals from groundwater wells.

Oklahoma, Texas↗
USGS2026

Shallow hydrogeologic framework of the Tully Valley mudboil area, Onondaga County, New York

Mudboils have been documented in the Tully Valley in southern Onondaga County, New York, since the late 1890s. Sediment-laden water from the mudboils flows into Onondaga Creek, which empties into Onondaga Lake at Syracuse 15 miles to the north. Turbidity from the mudboils has degraded the water quality of Onondaga Creek despite a series of mitigation efforts that began in the early 1990s. Turbidity mitigation actions presently (2025) being considered include creek relocation and offline sediment settling. In support of these proposed actions during 2021–23, the U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, U.S. Environmental Protection Agency, Onondaga Nation, Onondaga Environmental Institute, and Central New York Regional Planning and Development Board, collected and analyzed geologic, hydrologic, geophysical, and geotechnical data to characterize the shallow hydrogeology along four proposed creek-relocation paths and in the proposed offline settling basin area. The investigation indicated that the four proposed creek-relocation paths, two east of Onondaga Creek and two west of Onondaga Creek, are underlain by sediments including muck, alluvium, mudboil deposits, alluvial-fan sand and gravel, and lacustrine fines. The proposed excavations would penetrate partially to fully saturated conditions: generally, the water table is shallow near the creek and deep on the alluvial fans. The shallowest excavation, about 5 feet below land surface, would be near the creek and primarily in alluvium, and the deepest excavation, as much as 30 feet below land surface, would be in the alluvial-fan deposits. Brackish waters would be penetrated by proposed channel excavations on the eastern side of Onondaga Creek in an area downgradient from a potentially leaking historical salt-exploration borehole and near the main mudboil area. Excavation in these areas likely would provide a continuous source of brackish groundwater to the relocated creek. Proposed channel excavations of muck, soft to very soft lacustrine fines, and mudboil-type sediments in mudboil and suspected mudboil areas would pose an excavation and slope stability challenge and would have the greatest potential to create new mudboils. Proposed channel excavations below the water table on the Rattlesnake Gulf and Rainbow Creek alluvial fans would intercept groundwater and make the constructed streambank susceptible to seepage-induced slope instability. The substantial water-level fluctuation in the sediments of both alluvial fans would aggravate the stability condition. In addition, excavation on the Rattlesnake Gulf alluvial fan would have the potential to affect water-supply springs at the toe of the fan. The proposed offline settling basin area is in the northern part of the Rattlesnake Gulf alluvial fan. Natural and man-made diversions of Rattlesnake Gulf have resulted in saturated conditions in the general area of the proposed basin. The proposed offline settling basin would be excavated in, and berms would be constructed on, alluvial-fan deposits and lacustrine fines. In the proposed basin area, the alluvial deposits overlying the lacustrine fines are less than 10 feet thick. Excavation, berm construction, and loading of the saturated, soft to very soft lacustrine fines may be problematic and require soil strengthening.

New York↗
USGS1984

Techniques for estimating magnitude and frequency of floods on streams in Indiana

Equations are presented for estimating the magnitude and frequency of floods at ungaged sites on unregulated and nonurban streams in Indiana. The equations were developed by multiple-regression, analysis of basin characteristics and peak-flow statistical data from 242 gaged locations in Indiana, Ohio, and Illinois. The State of Indiana was divided into seven areas on the basis of the regression analysis. A set of equations for estimating peak discharges with recurrence intervals of 2, IO, 25, 50, and 100 years was developed for each area. Significant basin characteristics in the equations are drainage area, channel length, channel slope, mean annual precipitation, storage, precipitation intensity, and a runoff coefficient. Standard errors of estimate for the equations range from 24 to 45 percent. Methods are also presented for estimating flood magnitude and frequency at sites on gaged streams. Flood-frequency data based on observed peaks are given for 270 gaged locations. Twenty of these are on regulated streams, and six are on urban streams. Basin characteristics are also included car 245 of the gaged locations on unregulated and nonurban streams. No techniques are given for estimating flood magnitude and frequency at ungaged sites on regulated or urban streams. A rainfall-runoff model was tlsed to synthesize long-term peak data at 11 gaged locations on small streams. Flood-frequency curves developed from the long-term synthetic data were combined with curves based on short-term observed data to provide weighted estimates of flood magnitude and frequency at the rainfall-runoff stations.

Indiana↗
USGS1983

Use of digital land-cover data from the Landsat satellite in estimating streamflow characteristics in the Cumberland Plateau of Tennessee

Characteristics of Streamflow at ungaged sites in coal-mining areas of the Cumberland Plateau may be estimated with regression equations that relate streamflow characteristics to the physiographic and climatic characteristics of the corresponding drainage basins. An experiment was performed to compare the accuracy of equations using basin characteristics derived from maps and climatological records (control group equations) with the accuracy of equations using basin characteristics derived from digital processing of Landsat spectral data as well as maps and climatological records (experimental group equations). Results of this experiment show that (with the exception of low flows and four of six annual peak logarithms) drainage area can explain more than 90 percent of the variance in all streamflow characteristics in both groups of equations. Seventeen of 39 experimental group equations that have two basin characteristics each are different from the corresponding control group equations. Five of the 17 differing experimental group equations have no counterpart in the control group because of lack of significance in control group basin characteristics, seven are slightly more accurate, and five are not measurably different in accuracy. When the equations in both groups are arranged into six flow categories, there is no substantial difference in accuracy between equations using basin characteristics derived from maps and climatological records (control group) and equations using basin characteristics derived from Landsat tapes as well as maps and climatological records (experimental group) for this particular study area, the Cumberland Plateau of Tennessee.

Tennessee↗
USGS2025

Hydrogeology of unconsolidated and bedrock aquifers along the Salmon River, including Malone, Franklin County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, investigated the hydrogeology of the unconsolidated and bedrock aquifers along the Salmon River corridor in northern Franklin County, New York. The study area covers roughly 147 square miles and includes the village of Malone and parts of the Towns of Malone, Bellmont, Burke, Constable, Westville, Bangor, Duane, and Franklin, New York. Groundwater is the primary source for water supply within the study area. Eighty-three percent of all residents use public water supplied from two production wells that draw water from a thick, highly productive sand-and-gravel aquifer likely receiving induced infiltration from the Salmon River. Twenty-four percent of the 187 verified domestic wells in the study area outside of the production well service area boundary are screened in typically discontinuous deposits of stratified sand and sand and gravel. Seventy-six percent of the wells are completed in bedrock aquifers including the Potsdam Sandstone, metamorphic rocks, Theresa Formation, and unknown bedrock. Characterizing and understanding potential groundwater resources is critical for protecting the quality of the groundwater. The information in this report may be used to guide delineation of groundwater contributing areas, assess potential threats to aquifers from both point and nonpoint sources, respond to contamination from spills or leaks from underground storage facilities or other sources, and to support assessments for future development of municipal water supplies.

New York↗
USGS2024

Chloride concentrations in groundwater from the western part of the Southern Hills regional aquifer system, Louisiana, 2021–22

Groundwater is heavily used for public supply and industrial uses in the Baton Rouge, Louisiana, area. Lowered water levels resulting from groundwater withdrawals have induced the movement of saltwater towards wells in East Baton Rouge and West Baton Rouge Parishes. Saltwater intrusion has the potential to affect water supply infrastructure, reduce water availability for some uses, and increase treatment costs. To document current conditions, samples were collected from 161 wells screened in 10 aquifers of the Southern Hills regional aquifer system during November 2021 through February 2022. The results were compared with historical data to identify where chloride concentrations are increasing, which could indicate that saltwater intrusion is occurring. Saltwater intrusion, to varying degrees and areal extents, was observed in most of the 10 aquifers. The limited availability of monitoring wells near or within some of the known saltwater plume areas restricts tracking of the movement or delineation of the plumes’ current extents.

Louisiana↗
USGS2024

Conceptual hydrogeologic framework and groundwater budget near the southeastern part of Puget Sound, Washington

More than 1 million people live within the active model area (AMA) in the southeastern part of the lowlands surrounding Puget Sound, or Puget Lowland, Washington, and groundwater is the source for approximately one-half of their public, domestic, and irrigation water demands. The 887-square-mile AMA, located in King and Pierce Counties, represents the area of analysis for the conceptual hydrogeologic framework and numerical groundwater-flow models within the study area and includes the Puyallup River and Chambers-Clover Creek watersheds. To assess the potential hydrologic and anthropogenic impacts to groundwater and the connected surface-water resources, conceptual and numerical groundwater-flow models of groundwater flow were developed by the U.S. Geological Survey Washington Water Science Center in close cooperation with 18 water-resource agencies and stakeholders. This report presents information used to characterize the groundwater-flow system and the development of a numerical model in the AMA. Included are descriptions of the geology and conceptual hydrogeologic framework, groundwater levels and flow directions, groundwater recharge and discharge, numerical groundwater-flow model construction and results, and model limitations. The study area encompasses the western part of Pierce County and the southwestern part of King County, Washington. The study area extends south to the Nisqually River, southwest to Tanwax Creek, northeast to the Green River, and north through the valley near Auburn and adjacent uplands. It is bounded on the east by foothills of the Cascade Range, and on the northwest by Puget Sound.

Washington↗
USGS2025

A regional model comparison between MODPATH and MT3D of groundwater travel time distributions

Groundwater quality changes in wells and streams lag behind changes to land use due to groundwater travel times. Two contaminant transport methods were compared to assess differences in their simulated travel time distributions (TTDs) to streams and wells in the Wisconsin Central Sands. MODPATH simulates advective groundwater flow with particle tracking, while MT3D simulates age-mass using a finite difference solution without dispersion to allow for direct comparison of the two methods. MODPATH appropriately simulates groundwater TTDs from the water table to surface discharge but is subject to inaccuracies at weak-sink well cells due to the flow-model grid discretization and imprecise location of well discharge within well cells. MT3D better represents weak-sink well cells since it removes mass in proportion to the prescribed pumping rate, although travel time within well cells is neglected. Conversely, MT3D's treatment of surface water boundary cells is not as accurate as MODPATH because mass should be removed from the water table rather than the full cell volume. MT3D simulations of TTDs can also be confounded by the instantaneous vertical distribution of mass introduced throughout recharge cells instead of at the water table, which initiates mass along deeper flow paths. We evaluated 9 MODPATH and 13 MT3D implementations, generating differences in median travel times of up to 18 years. Both methods have strengths and weaknesses, with MT3D better representing weak-sink well cell behavior and MODPATH better representing surficial recharge and discharge. The effect of these characteristics on simulated TTDs, along with ideas for ameliorating method weaknesses, is discussed.

Wisconsin↗