SEARCH · USGS Science

Results for “Ground Water Modeling Investigation Report”

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

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25Linked to original sources

USGS1989

Geohydrology, simulation of ground-water flow, and ground-water quality at two landfills, Marion County, Indiana

Geologic, hydrologic, and water-quality data were collected at the Julietta and Tibbs-Banta landfills in Marion County. Both landfills were closed in the mid-1970's, and sewage sludge mixed with dirt was spread on the landfills in the mid-1980's as part of a revegetation project. The landfills were constructed in unconsolidated glacial sediments that consist of sand, gravel, silt, and clay. The maximum thickness of the sediments is 180 feet at Julietta and 100 feet at Tibbs-Banta. Both landfills are underlain by sand and gravel aquifers and are adjacent to gaining streams. Ground water flows toward and into the streams at each study area. Two sand and gravel aquifers were mapped at Julietta and four were mapped at TibbsBanta. The aquifers are separated in places by discontinuous clay layers. Ground-water-flow models, calibrated to simulate steady-state low-flow conditions, indicate that about 19,000 gallons of water per day move through the refuse at Julietta and 42,000 gallons per day move through the refuse at Tibbs-Banta. The Julietta model also indicates that recharge through the surface of the landfill is less than in the surrounding natural areas, probably because of the addition and compaction of the sludge/soil mixture. Concentrations of dissolved inorganic substances in ground-water samples indicate that leachate from both landfills is reaching the shallow aquifers. The effect on deeper aquifers is small because of the predominance of horizontal ground-water flow and discharge to the streams. Increases in almost all dissolved constituents were observed in shallow wells that are screened beneath and downgradient from the landfills. Several analyses, especially those for bromide, dissolved solids, and ammonia, were useful in delineating the plume of leachate at both landfills.

Indiana↗

USGS1995

Geohydrology and water quality of stratified-drift aquifers in the middle Merrimack River basin, south-central New Hampshire

The U.S. Geological Survey, in cooperation with the State of New Hampshire, Department of Environmental Services, Water Resources Division has assessed the geohydrology and water quality of stratified-drift aquifers in the middle Merrimack River basin in south-central New Hampshire. The middle Merrimack River basin drains 469 square miles; 98 square miles is underlain by stratified-drift aquifers. Saturated thickness of stratified drift within the study area is generally less than 40 feet but locally greater than 100 feet. Transmissivity of stratified-drift aquifers is generally less than 2,000 feet squared per day but locally exceeds 6, 000 feet squared per day. At present (1990), ground-water withdrawals from stratified drift for public supply are about 0.4 million gallons per day within the basin. Many of the stratified-drift aquifers within the study area are not developed to their fullest potential. The geohydrology of stratified-drift aquifers was investigated by focusing on basic aquifer properties, including aquifer boundaries; recharge, discharge, and direction of ground-water flow; saturated thickness and storage; and transmissivity. Surficial geologic mapping assisted in the determination of aquifer boundaries. Data from 757 wells and test borings were used to produce maps of water-table altitude, saturated thickness, and transmissivity of stratified drift. More than 10 miles of seismic-refraction profiling and 14 miles of seismic-reflection profiling were also used to construct the water table and saturated-thickness maps. Stratified-drift aquifers in the southern, western, and central parts of the study area are typically small and discontinuous, whereas aquifers in the eastern part along the Merrimack River valley are continuous. The Merrimack River valley aquifers formed in glacial Lakes Merrimack and Hooksett. Many other smaller discontinuous aquifers formed in small temporary ponds during deglaciation. A stratified-drift aquifer in Goffstown was analyzed for aquifer yield by use of a two-dimensional, finite-difference ground-water-flow model. Yield of the Goffstown aquifer was estimated to be 2.5 million gallons per day. Sensitivity analysis showed that the estimate of aquifer yield was most sensitive to changes in hydraulic conductivity. The amount of water induced into the aquifer from the Piscataquog River was most affected by changes in estimates of streambed conductance. Results of analysis of water samples from 10 test wells indicate that, with some exceptions, water in the stratified-drift aquifers generally meets U.S. Environmental Protection Agency primary and secondary drinking-water regulations. Water from two wells had elevated sodium concentrations, waterfront two wells had elevated concentrations of dissolved iron, and waterfront seven wells had elevated concentrations of manganese. Known areas of contamination were avoided during water-quality sampling.

New Hampshire↗

USGS2006

Hydraulic and solute-transport properties and simulated advective transport of contaminated ground water in a fractured-rock aquifer at the Naval Air Warfare Center, West Trenton, New Jersey, 2003

Volatile organic compounds, predominantly trichloroethylene and its degradation products, have been detected in ground water at the Naval Air Warfare Center (NAWC), West Trenton, New Jersey. An air-stripping pump-and-treat system has been in operation at the NAWC since 1998. An existing ground-water-flow model was used to evaluate the effect of a change in the configuration of the network of recovery wells in the pump-and-treat system on flow paths of contaminated ground water. The NAWC is underlain by a fractured-rock aquifer composed of dipping layers of sedimentary rocks of the Lockatong and Stockton Formations. Hydraulic and solute-transport properties of the part of the aquifer composed of the Lockatong Formation were measured using aquifer tests and tracer tests. The heterogeneity of the rocks causes a wide range of values of each parameter measured. Transmissivity ranges from 95 to 1,300 feet squared per day; the storage coefficient ranges from 9 x 10-5 to 5 x 10-3; and the effective porosity ranges from 0.0003 to 0.002. The average linear velocity of contaminated ground water was determined for ambient conditions (when no wells at the site are pumped) using an existing ground-water-flow model, particle-tracking techniques, and the porosity values determined in this study. The average linear velocity of flow paths beginning at each contaminated well and ending at the streams where the flow paths terminate ranges from 0.08 to 130 feet per day. As a result of a change in the pump-and-treat system (adding a 165-foot-deep well pumped at 5 gallons per minute and reducing the pumping rate at a nearby 41-foot-deep well by the same amount), water in the vicinity of three 100- to 165-foot-deep wells flows to the deep well rather than the shallower well.

New Jersey↗

USGS2001

A field and statistical modeling study to estimate irrigation water use at Benchmark Farms study sites in southwestern Georgia, 1995-96

A benchmark irrigation monitoring network of farms located in a 32-county area in southwestern Georgia was established in 1995 to improve estimates of irrigation water use. A stratified random sample of 500 permitted irrigators was selected from a data base--maintained by the Georgia Department of Natural Resources, Georgia Environmental Protection Division, Water Resources Management Branch--to obtain 180 voluntary participants in the study area. Site-specific irrigation data were collected at each farm using running-time totalizers and noninvasive flowmeters. Data were collected and compiled for 50 farms for 1995 and 130 additional farms for the 1996 growing season--a total of 180 farms. Irrigation data collected during the 1996 growing season were compiled for 180 benchmark farms and used to develop a statistical model to estimate irrigation water use in 32 counties in southwestern Georgia. The estimates derived were developed from using a statistical approach know as "bootstrap analysis" that allows for the estimation of precision. Five model components--whether-to-irrigate, acres irrigated, crop selected, seasonal-irrigation scheduling, and the amount of irrigation applied--compose the irrigation model and were developed to reflect patterns in the data collected at Benchmark Farms Study area sites. The model estimated that peak irrigation for all counties in the study area occurred during July with significant irrigation also occurring during May, June, and August. Irwin and Tift were the most irrigated and Schley and Houston were the least irrigated counties in the study area. High irrigation intensity primarily was located along the eastern border of the study area; whereas, low irrigation intensity was located in the southwestern quadrant where ground water was the dominant irrigation source. Crop-level estimates showed sizable variations across crops and considerable uncertainty for all crops other than peanuts and pecans. Counties having the most irrigated acres showed higher variations in annual irrigation than counties having the least irrigated acres. The Benchmark Farms Study model estimates were higher than previous irrigation estimates, with 20 percent of the bias a result of underestimating irrigation acreage in earlier studies. Model estimates showed evidence of an upward bias of about 15 percent with the likely cause being a misrepresented inches-applied model. A better understanding of the causes of bias in the model could be determined with a larger irrigation sample size and increased substantially by automating the reporting of monthly totalizer amounts.

Georgia↗

USGS1982

Method for estimating historical irrigation requirements from ground water in the High Plains in parts of Colorado, Kansas, Nebraska, New Mexico, Oklahoma, South Dakota, Texas, and Wyoming

Historical information on ground-water pumpage for irrigation is required for use with a computer model of the High Plains aquifer. Available historical data on ground-water pumpage in the High Plains were inadequate for use with the computer model. Consequently, a method was developed to estimate historical pumpage data. Two principal components were used to develop estimates of ground-water pumpage: (1) Acreage irrigated; and (2) irrigation demand. The Blaney-Criddle formula was used to calculate consumptive-use requirements for irrigated crops grown on the High Plains. The irrigation demand for each crop was estimated by subtracting precipitation available to the crop from the consumptive-use requirement of the crop. Irrigation demands were combined with irrigated acreages to estimate the volume of irrigation water required for areas in the High Plains. Estimates of irrigation water requirements were compiled for 1949, 1954, 1959, 1964, 1969, 1974, and 1978. Irrigation pumpage for use with the computer model can be estimated from the irrigation water requirements by applying a factor to account for irrigation system efficiency.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗

USGS1976

Mathematical model of the West Bolsa Ground-water Basin, San Benito County, California

Simulation of the West Bolsa ground-water basin hydrology in California had provided values of basin recharge and discharge and nodally distributed values of transmissivity and storage coefficient. Average net recharge from April 1945 to March 1969 was 6.2 cubic feet per second and occurred as subsurace recharge and infiltration of rain and minor streamflow. Discharge from the basin during the same period was 8.1 cubic feet per second and occurred as pumping and leakage from confined parts of the basin. Values of transmissivity used in the model generally range from 3,300 to 20,000 feet squared per day. Values of storage coefficient used in the model range from 0.0005 to 0.10. (Woodard-USGS)

California↗

USGS1984

Ground water in the Fresno area, California; preliminary report

The U.S. Geological Survey, in cooperation with Fresno County, is developing a computer model of the unconfined aquifer system in the Fresno area. In the area, consolidated rocks are overlain by unconsolidated deposits. The unconsolidated deposits are divided into a lower fine-grained unit and an overlying coarse-grained unit. The system is assumed to be a single unconfined unit. The top of the system is the water table and the bottom is the base of the coarse-grained deposits. The northeastern boundary is the only natural lateral hydrologic boundary. The other lateral boundaries are placed so that errors along them will not affect heads in the area of interest. Initial heads for 1943, hydraulic conductivity, and bottom elevations were used in a steady-state model that generated net flux values. It is anticipated that most of the remaining effort will center on model calibration. (USGS)

Water-Resources Investigations Report↗

USGS1994

The computer model Hydrotherm, a three-dimensional finite-difference model to simulate ground-water flow and heat transport in the temperature range of 0 to 1,200 degrees C

Quantitative modeling of the deep parts of magmatic- hydrothermal systems has been limited by the lack of publicly available, documented computer models for multiphase, high-temperature flow. This report documents HYDROTHERM, a finite-difference model for three-dimensional, multiphase flow of pure water and heat over a temperature range of 0 to 1,200 degrees Celsius and a pressure range of 0.5 *0 10,000 bars.HYDROTHEP*M is an extension of multiphase geothermal models developed by the U.S. Geological Survey in the 1970's. It solves numerical approximations to mass- and energy-balance equations that are posed in terms of pressure and enthalpy. Major assumptions are that the rock matrix can be treated as a porous medium; that water and rock are in thermal equilibrium; and that capillary pressure is negligible.

Water-Resources Investigations Report↗

USGS2001

Approach for delineation of contributing areas and zones of transport to selected public-supply wells using a regional ground-water flow model, Palm Beach County, Florida

Rapid urban development and population growth in Palm Beach County, Florida, have been accompanied with the need for additional freshwater withdrawals from the surficial aquifer system. To maintain water quality, County officials protect capture areas and determine zones of transport of municipal supply wells. A multistep process was used to help automate the delineation of wellhead protection areas. A modular ground-water flow model (MODFLOW) Telescopic Mesh Refinement program (MODTMR) was used to construct an embedded flow model and combined with particle tracking to delineate zones of transport to supply wells; model output was coupled with a geographic information system. An embedded flow MODFLOW model was constructed using input and output file data from a preexisting three-dimensional, calibrated model of the surficial aquifer system. Three graphical user interfaces for use with the geographic information software, ArcView, were developed to enhance the telescopic mesh refinement process. These interfaces include AvMODTMR for use with MODTMR; AvHDRD to build MODFLOW river and drain input files from dynamically segmented linear (canals) data sets; and AvWELL Refiner, an interface designed to examine and convert well coverage spatial data layers to a MODFLOW Well package input file. MODPATH (the U.S. Geological Survey particle-tracking postprocessing program) and MODTOOLS (the set of U.S. Geological Survey computer programs to translate MODFLOW and MODPATH output to a geographic information system) were used to map zones of transport. A steady-state, five-layer model of the Boca Raton area was created using the telescopic mesh refinement process and calibrated to average conditions during January 1989 to June 1990. A sensitivity analysis of various model parameters indicates that the model is most sensitive to changes in recharge rates, hydraulic conductivity for layer 1, and leakance for layers 3 and 4 (Biscayne aquifer). Recharge (58 percent); river (canal) leakance (29 percent); and inflow through the northern, western, and southern prescribed flux model boundaries (10 percent) represent the major inflow components. Principal outflow components in the Boca Raton well field area include well discharge (56 percent), river (canal) leakance (27 percent), and water that discharges along the coast (10 percent). A particle-tracking analysis using MODPATH was conducted to better understand well-field ground-water flow patterns and time of travel. MODTOOLS was used to construct zones-of-transport spatial data for municipal supply wells. Porosity estimates were uniformly increased to study the effect of porosity on zones of transport. Where porosity was increased, the size of the zones of transport were shown to decrease.

Water-Resources Investigations Report↗

USGS2004

Documentation of revisions to the regional aquifer system analysis model of the New Jersey coastal plain

The model, which simulates flow in the New Jersey Coastal Plain sediments, developed for the U.S. Geological Survey Regional Aquifer System Analysis (RASA) program was revised. The RASA model was revised with (1) a rediscretization of the model parameters with a finer cell size, (2) a spatially variable recharge rate that is based on rates determined by recent studies and, (3) ground-water withdrawal data from 1981 to 1998. The RASA model framework, which subdivided the Coastal Plain sediments into 10 aquifers and 9 confining units, was preserved in the revised model. A transient model that simulates flow conditions from January 1, 1968 to December 31, 1998, was constructed using 21 stress periods. The model was calibrated by attempting to match the simulated results with (1) estimated base flow for five river basins, (2) measured water levels in long-term hydrographs for 28 selected observation wells, and (3) potentiometric surfaces in the model area for 1978, 1983, 1998, 1993, and 1998 conditions. The estimated and simulated base flow in the five river basins compare well. In general, the simulated water levels matched the interpreted potentiometric surfaces and the measured water levels of the hydrographs within 25 feet.

Water-Resources Investigations Report↗

USGS1995

Assessment of intrinsic bioremediation of jet fuel contamination in a shallow aquifer, Beaufort, South Carolina

Field and laboratory studies show that microorganisms indigenous to the ground-water system underlying Tank Farm C, Marine Corps Air Station Beaufort, S.C., degrade petroleum hydrocarbons under aerobic and anaerobic conditions. Under aerobic conditions, sediments from the shallow aquifer underlying the site mineralized radiolabeled (14C) toluene to 14CO2 with first-order rate constants of about -0.29 per day. Sediments incubated under anaerobic conditions mineralized radiolabeled toluene more slowly, with first-order rate constants of -0.001 per day. Although anaerobic rates of biodegradation are low, they are significant in the hydrologic and geochemical context of the site. Because of low hydraulic conductivities (1.9-9.1 feet per day) and low hydraulic gradients (about 0.004 feet per feet), ground water flows slowly (approximately 20 feet per year) at this site. Furthermore, aquifer sediments contain organic-rich peat that has a high sorptive capacity. Under these conditions, hydrocarbon contaminants have moved no further than 10 feet downgradient of the jet fuel free product. Digital solute-transport simulations, using the range of model parameters measured at the site, show that dissolved contaminants will be completely degraded before they are discharged from the aquifer into adjacent surface-water bodies. These results show that natural attenuation processes are containing the migration of soluble hydrocarbons, and that intrinsic bioremediation is a potentially effective remedial strategy at this site.

South Carolina↗

USGS2002

Preliminary estimates of spatially distributed net infiltration and recharge for the Death Valley region, Nevada-California

A three-dimensional ground-water flow model has been developed to evaluate the Death Valley regional flow system, which includes ground water beneath the Nevada Test Site. Estimates of spatially distributed net infiltration and recharge are needed to define upper boundary conditions. This study presents a preliminary application of a conceptual and numerical model of net infiltration. The model was developed in studies at Yucca Mountain, Nevada, which is located in the approximate center of the Death Valley ground-water flow system. The conceptual model describes the effects of precipitation, runoff, evapotranspiration, and redistribution of water in the shallow unsaturated zone on predicted rates of net infiltration; precipitation and soil depth are the two most significant variables. The conceptual model was tested using a preliminary numerical model based on energy- and water-balance calculations. Daily precipitation for 1980 through 1995, averaging 202 millimeters per year over the 39,556 square kilometers area of the ground-water flow model, was input to the numerical model to simulate net infiltration ranging from zero for a soil thickness greater than 6 meters to over 350 millimeters per year for thin soils at high elevations in the Spring Mountains overlying permeable bedrock. Estimated average net infiltration over the entire ground-water flow model domain is 7.8 millimeters per year.To evaluate the application of the net-infiltration model developed on a local scale at Yucca Mountain, to net-infiltration estimates representing the magnitude and distribution of recharge on a regional scale, the net-infiltration results were compared with recharge estimates obtained using empirical methods. Comparison of model results with previous estimates of basinwide recharge suggests that the net-infiltration estimates obtained using this model may overestimate recharge because of uncertainty in modeled precipitation, bedrock permeability, and soil properties for locations such as the Spring Mountains. Although this model is preliminary and uncalibrated, it provides a first approximation of the spatial distribution of net infiltration for the Death Valley region under current climatic conditions.

Water-Resources Investigations Report↗

USGS1999

Water resources of the Batavia Kill basin at Windham, Greene County, New York

The water resources of a 27.6-square-mile section of the Batavia Kill Basin near the village of Windham, N.Y., which has undergone substantial development, were evaluated. The evaluation entailed (1) estimation of the magnitude and distribution of several hydrologic components, including recharge, (2) measurement of discharge and chemical quality of the Batavia Kill and selected tributaries, (3) analysis of ground-water flow and chemistry, and (4) a conceptualization of the ground-water flow system. The region consists of deeply dissected, relatively flat-lying, clastic sedimentary sequences variably overlain by as much as 120 feet of glacial deposits. The types of bedrock fractures and their distribution in the Batavia Kill valley are consistent with valley stress-relief characteristics. Till predominates in the uplands, and stratified drift typically dominates within the valley of the Batavia Kill and the lower section of its largest tributary valley (Mitchell Hollow). Fractured bedrock is the most commonly used water source within the study area. The areas of highest yielding bedrock generally are with valleys, where the shallow fractures are saturated. Stratified-drift aquifers are also limited to the largest valleys; the greatest saturated thicknesses are in the Batavia Kill valley at Windham. A conceptual model of ground-water flow within the study areas suggests that the zones of most active flow are shallow fractured bedrock in upland areas and the shallow stratified drift in the largest valleys. The hydrogeologic system has been altered by development; major effects include (1) chemical alteration of natural ground-water and surface-water quality by point- and nonpoint-source contaminants, (2) hydraulic interconnection of other-wise isolated bedrock fractures by wellbores, and (3) drawdowns in wells within the Batavia Kill valley by pumping from the bedrock aquifer. Water resource development of the most promising unconsolidated aquifer beneath Windham may be precluded by the potential for contamination by leachate from an abandoned landfill, road-salt stockpiles, and domestic septic systems in the area.

New York↗

USGS2000

Hydrogeology, hydrologic budget, and water chemistry of the Medina Lake area, Texas

A three-phase study of the Medina Lake area in Texas was done to assess the hydrogeology and hydrology of Medina and Diversion Lakes combined (the lake system) and to determine what fraction of seepage losses from the lake system might enter the regional ground-water-flow system of the Edwards and (or) Trinity aquifers. Phase 1 consisted of revising the geologic framework for the Medina Lake area. Results of field mapping show that the upper member of the Glen Rose Limestone underlies Medina Lake and the intervening stream channel from the outflow of Medina Lake to the midpoint of Diversion Lake, where the Diversion Lake fault intersects Diversion Lake. A thin sequence of strata consisting primarily of the basal nodular and dolomitic members of the Kainer Formation of the Edwards Group, is present in the southern part of the study area. On the southern side of Medina Lake, the contact between the upper member of the Glen Rose Limestone and the basal nodular member is approximately 1,000 feet above mean sea level, and the contact between the basal nodular member and the dolomitic member is approximately 1,050 feet above mean sea level. The most porous and permeable part of the basal nodular member is about 1,045 feet above mean sea level. At these altitudes, Medina Lake is in hydrologic connection with rocks in the Edwards aquifer recharge zone, and Medina Lake appears to lose more water to the ground-water system along this bedding plane contact. Hydrologic budgets calculated during phase 2 for Medina Lake, Diversion Lake, and Medina/Diversion Lakes combined indicate that: (1) losses from Medina and Diversion Lakes can be quantified; (2) a portion of those losses are entering the Edwards aquifer; and (3) losses to the Trinity aquifer in the Medina Lake area are minimal and within the error of the hydrologic budgets. Hydrologic budgets based on streamflow, precipitation, evaporation, and change in lake storage were used to quantify losses (recharge) to the ground-water system from Medina Lake, Diversion Lake, and Medina/Diversion Lakes combined during October 1995–September 1996. Water losses from Medina Lake to the Edwards/Trinity aquifers ranged from -14.0 to 135 acre-feet per day; Diversion Lake ranged from -1.2 to 93.1 acre-feet per day; and Medina/Diversion Lakes combined ranged from 36.1 to 119 acre-feet per day. Monthly average recharge during December 1995–July 1996 was estimated using an alternative method developed during this study (current study method) and compared to monthly average recharge during December 1995–July 1996 estimated using the existing USGS method and the Trans-Texas method. Recharge to the Edwards aquifer estimated using the current study method was about 69 and 73 percent of the recharge estimated using the USGS and Trans-Texas methods, respectively. The USGS and Trans-Texas methods overestimated recharge from Medina Lake compared to the recharge estimated with the current study method when Medina Lake stage was between about 1,027 and 1,032 feet above mean sea level and underestimated recharge from Medina Lake when lake stage was between about 1,036 and 1,045 feet above mean sea level. The USGS and Trans-Texas methods underestimated recharge from Diversion Lake compared to the recharge estimated with the current study method when Diversion Lake stage was greater than 913 feet above mean sea level and overestimated recharge from Diversion Lake when lake stage was less than 913 feet above mean sea level. The water quality of Medina Lake and Medina River and in selected wells and springs in the Edwards and Trinity aquifers was characterized during phase 3 of the study. Environmental isotope analyses and geochemical modeling also were used to determine where water losses from the lake system might be entering the ground-water-flow system. Isotopic ratios of deuterium, oxygen, and strontium were analyzed in selected surface-water, lake-water, and ground-water samples to trace the isotopic “signature” of the lake water as it mixes with the ground water and to determine the fraction of lake water and ground water in selected Edwards aquifer wells. Isotopic data and geochemical modeling were used to show that lake water is moving into the Edwards aquifer in two fault blocks in the eastern Medina storage unit. One fault block is bounded on the north by the Vandenburg School fault and on the south by the Haby Crossing fault, and the second fault block is bounded on the north by the Diversion Lake fault and on the south by the Haby Crossing fault. In selected Edwards aquifer wells located southwest of Medina Lake and west of Diversion Lake, the proportion of lake water ranged from about 10 to 45 percent. Geochemical modeling using NETPATH confirms the degree of mixing between lake water and aquifer water shown by the isotopes.

Texas↗

USGS1985

A two-constituent solute-transport model for ground water having variable density

A numerical model has been developed to simulate solute transport and dispersion of either one or two constituents in groundwater where there is two-dimensional, density-dependent flow. The model is a modified version of the one documented by Konikow and Bredehoeft (1978), which uses finite-difference methods and the method of characteristics to solve the flow and transport equations. The model was tested on an idealized seawater intrusion problem for which an analytical solution has been developed. The results were nearly identical to those of other numerical models tested on the same problem. A description of the formats for the input data, a sample of input and output for a two-constituent example problem, and a listing of the Fortran program are presented. (Author 's abstract)

Water-Resources Investigations Report↗

USGS1987

Evaluation of availability of water from drift aquifers near the Pomme de Terre and Chippewa rivers, western Minnesota

Ground-water flow in the confined- and unconfined-drift aquifers near Appleton and Benson, Minnesota, was simulated with a three-dimensional finite-difference ground-water-flow model. Model results indicate that 98 percent of the total inflow to the modeled area is from precipitation. Of the total outflow, 38 percent is ground-water discharge to the Pom me de Terre and Chippewa Rivers, 36 percent is evapotranspiration, 17 percent is ground-water pumpage, and 8 percent is ground-water discharge to the Minnesota River. The model was used to simulate the effects of below-normal precipitation (drought) and hypothetical increases in ground-water development. Model results indicate that reduced recharge and increased pumping during a three-year extended drought probably would lower water levels 2 to 6 feet regionally in the surficial aquifer and in the Appleton and Benson-middle aquifers and as much as 11 feet near aquifer boundaries. Ground-water discharge to the Pomme de Terre and Chippewa Rivers in the modeled area probably would be reduced during the simulated drought by 15.2 and 7.4 cubic feet per second, respectively, compared to 1982 conditions. The addition of 30 hypothetical wells in the Benson-middle aquifer near Benson, pumping a total of 810 million gallons per year, resulted in water-level declines of as much as 1.3 and 2.7 feet in the surficial and Benson-middle aquifers, respectively. The addition of 28 hypothetical wells in the Appleton aquifer east and southeast of Appleton, pumping a total of 756 million gallons per year, lowered water levels in the surficial and Appleton confined aquifers as much as 5 feet.

Minnesota↗

USGS1994

Hydrogeology and simulation of ground-water flow at U.S. Marine Corps Air Station, Cherry Point, North Carolina, 1987-90

Geophysical and lithologic well-log data from 30 wells and chloride data, and water-level data from oil-test wells, supply wells, and observation wells were evaluated to define the hydrogeologic framework at the U.S. Marine Corps Air Station, Cherry Point, North Carolina. Elements of the hydrogeologic framework important to this study include six aquifers and their respective confining units. In descending order, these aquifers are the surficial, Yorktown, Pungo River, upper and lower Castle Hayne, and Beaufort. The upper and lower Castle Hayne and Beaufort aquifers and related confining units are relatively continuous throughout the study area. The surficial, Yorktown, Pungo River, and upper and lower Castle Hayne aquifers contain freshwater. The upper and lower Castle Hayne aquifers serve as the Air Station?s principal supply of freshwater. However, the lower Castle Hayne aquifer contains brackish water near its base and there is potential for upward movement of this water to supply wells completed in this aquifer. The potential for brackish-water encroachment is greatest if wells are screened too deep in the lower Castle Hayne aquifer or if pumping rates are too high. Lateral movement of brackish water into aquifers incised by estuarine streams is also possible if ground-water flow gradients toward these bodies are reversed by pumping. The potential for the reversed movement of water from the surficial aquifer downward to the water-supply aquifer is greatest in areas where clay confining units are missing. These missing clay units could indicate the presence of a paleochannel of the Neuse River. A quasi three-dimensional finite-difference ground-water flow model was constructed and calibrated to simulate conditions at and in the vicinity of the Air Station for the period of 1987-90. Comparisons of 94 observed and computed heads were made, and the average difference between them is -0.2 feet with a root mean square error of 5.7 feet. An analysis was made to evaluate the sensitivity of the model to the absence of the Yorktown and Pungo River confining units in a 1-square-mile area in the southern part of the Air Station. This analysis resulted in a maximum simulated head increase of 2 feet in one 0.11-square-mile model cell in the Pungo River aquifer.

North Carolina↗

USGS2008

Simulated effects of ground-water withdrawals and artificial recharge on discharge to streams, springs, and riparian vegetation in the Sierra Vista Subwatershed of the Upper San Pedro Basin, southeastern Arizona

In the context of ground-water resources, “capture” or “streamflow depletion” refers to withdrawal-induced changes in inflow to or outflow from an aquifer. These concepts are helpful in understanding the effects of long-term development of ground-water resources. For the Upper San Pedro Basin in Arizona, USA and Sonora, Mexico, a recently developed ground-water flow model is available to help quantify capture of water from the river and riparian system. A common method of analysis is to compute curves of capture and aquifer-storage change for a range of time at select points of interest. This study, however, presents results of a method to show spatial distributions of total change in inflow and outflow from withdrawal or injection for select times of interest. The mapped areal distributions show the effect of a single well in terms of the ratio of the change in boundary flow rate to rate of withdrawal or injection by the well. To the extent that the system responds linearly to ground-water withdrawal or injection, fractional responses in the mapped distributions can be used to quantify response for any withdrawal or injection rate. Capture distributions calculated using the Upper San Pedro model include response to (1) withdrawal in the lower basin-fill aquifer for times of 10 and 50 years following the initiation of pumping from predevelopment conditions and (2) artificial recharge to the water table in the area underlain by the lower basin-fill aquifer after 10 and 50 years. The mapped distributions show that response to withdrawals and injections is greatest near the river/riparian system. Presence of clay layers in the vertical interval between withdrawal locations and the river/riparian system, however, can delay the response.

Arizona;Sonora↗