Geology Reports⌕ Search

SEARCH · Geology Reports

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.

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

At least 937 records · Page 52Linked to original sources
USGS1999

Use of environmental tracers to evaluate ground-water age and water-quality trends in a buried-valley aquifer, Dayton area, southwestern Ohio

Chlorofluorocarbons (CFC method) and tritium and helium isotopes (3H-3He method) were used as environmental tracers to estimate ground-water age in conjunction with efforts to develop a regional ground-water flow model of the buried-valley aquifer in the Dayton area, southwestern Ohio. This report describes results of CFC and water-quality sampling, summarizes relevant aspects of previously published work, and describes the use of 3H-3He ages to characterize temporal trends in ground-water quality of the buried-valley aquifer near Dayton, Ohio. Results of CFC sampling indicate that approximately 25 percent of the 137 sampled wells were contaminated with excess CFC's that rendered the ground water unsuitable for age dating. Evaluation of CFC ages obtained for the remaining samples indicated that the CFC compounds used for dating were being affected by microbial degradation. The degradation occurred under anoxic conditions that are found in most parts of the buried-valley aquifer. As a result, ground-water ages derived by the CFC method were too old and were inconsistent with measured tritium concentrations and independently derived 3H-3He ages. Limited data indicate that dissolved methane may play an important role in the degradation of the CFC's. In contrast, the 3H-3He technique was found to yield ground-water ages that were chemically and hydrologically reasonable. Ground-water ages derived by the 3H-3He technique were compared to values for selected water- quality characteristics to evaluate temporal trends in ground-water quality in the buried- valley aquifer. Distinct temporal trends were not identified for pH, alkalinity, or calcium and magnesium because of rapid equilibration of ground-water with calcite and dolomite in aquifer sediments. Temporal trends in which the amount of scatter and the number of outlier concentrations increased as ground-water age decreased were noted for sodium, potassium, boron, bromide, chloride, ammonia, nitrate, phosphate, sulfate, and organic carbon. Elevated concentrations of these constituents in shallow ground water are probably related to human activities. Temporal trends in which concentrations declined as ground-water age increased may reflect natural processes that reduce constituent concentrations to low levels. For example, the absence of nitrate detections in ground water recharged before 1980 may indicate natural removal of nitrate by bacterially mediated denitrification. Temporal trends observed for dissolved oxygen, iron, nitrate and silica indicate that these constituents may help identify recently (post-1990) recharged ground water.

Ohio↗
USGS2007

Simulation of streamflow and estimation of ground-water recharge in the Upper Cibolo Creek Watershed, south-central Texas, 1992-2004

A watershed model (Hydrological Simulation Program?FORTRAN) was developed, calibrated, and tested by the U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, San Antonio River Authority, San Antonio Water System, and Guadalupe-Blanco River Authority, to simulate streamflow and estimate ground-water recharge in the upper Cibolo Creek watershed in south-central Texas. Rainfall, evapotranspiration, and streamflow data were collected during 1992?2004 for model calibrations and simulations. Estimates of average ground-water recharge during 1992?2004 from simulation were 79,800 acre-feet (5.47 inches) per year or about 15 percent of rainfall. Most of the recharge (about 74 percent) occurred as infiltration of streamflow in Cibolo Creek. The remaining recharge occurred as diffuse infiltration of rainfall through the soil and rock layers and karst features. Most recharge (about 77 percent) occurred in the Trinity aquifer outcrop. The remaining 23 percent occurred in the downstream part of the watershed that includes the Edwards aquifer recharge zone (outcrop). Streamflow and recharge in the study area are greatly influenced by large storms. Storms during June 1997, October 1998, and July 2002 accounted for about 11 percent of study-area rainfall, 61 percent of streamflow, and 16 percent of the total ground-water recharge during 1992?2004. Annual streamflow and recharge also were highly variable. During 1999, a dry year with about 16 inches of rain and no measurable runoff at the watershed outlet, recharge in the watershed amounted to only 0.99 inch compared with 13.43 inches during 1992, a relatively wet year with about 54 inches of rainfall. Simulation of flood-control/recharge-enhancement structures showed that certain structures might reduce flood peaks and increase recharge. Simulation of individual structures on tributaries showed relatively little effect. Larger structures on the main stem of Cibolo Creek were more effective than structures on tributaries, both in terms of flood-peak reduction and recharge enhancement. One simulated scenario that incorporated two main-stem structures resulted in a 37-percent reduction of peak flow at the watershed outlet and increases in stream-channel recharge of 6.6 percent in the Trinity aquifer outcrop and 12.6 percent in the Edwards aquifer (recharge zone) outcrop.

Texas↗
USGS1996

Geochemical and isotopic composition of ground water with emphasis on sources of sulfate in the upper Floridan Aquifer in parts of Marion, Sumter, and Citrus counties, Florida

In inland areas of northwest central Florida, sulfate concentrations in the Upper Floridan aquifer are extremely variable and sometimes exceed drinking water standards (250 milligrams per liter). This is unusual because the aquifer is unconfined and near the surface, allowing for active recharge. The sources of sulfate and geochemical processes controlling ground-water composition were evaluated in this area. Water was sampled from thirty-three wells in parts of Marion, Sumter, and Citrus Counties, within the Southwest Florida Water Management District; these included at least a shallow and a deep well at fifteen separate locations. Ground water was analyzed for major ions, selected trace constituents, dissolved organic carbon, and stable isotopes (sulfur-34 of sulfate and sulfide, carbon-13 of inorganic carbon, deuterium, and oxygen-18). Sulfate concentrations ranged from less than 0.2 to 1,400 milligrams per liter, with higher sulfate concentrations usually in water from deeper wells. The samples can be categorized into a low sulfate group (less than 30 milligrams per liter) and a high sulfate group (greater than 30 milligrams per liter). For the high sulfate water, concentrations of calcium and magnesium increased concurrently with sulfate. Chemical and isotopic data and mass-balance modeling indicate that the composition of high sulfate waters is controlled by dedolomitization reactions (dolomite dissolution and calcite precipitation, driven by dissolution of gypsum). Gypsum occurs deeper in the aquifer than open intervals of sampled wells. Upward flow has been documented in deeper parts of the aquifer in the study area, which may be driven by localized discharge areas or rapid flow in shallow parts of the aquifer. Mixing between shallow ground water and sulfate-rich water that dissolved gypsum at the base of the aquifer is probably responsible for the range of concentrations observed in the study area. Other solutes that increased with sulfate apparently originate from the gypsum itself, from other mineral assemblages found deeper in the aquifer in association with gypsum, and from residual seawater from less- flushed, deeper parts of the aquifer. These ions are subsequently transported with sulfate to shallower parts of the aquifer where gypsum is not present. The composition of low sulfate ground water is controlled by differences in the extent of microbially mediated reactions, which produce carbon dioxide. This, in turn, influences the extent of calcite dissolution. Ground waters which underwent limited microbial reactions contained dissolved oxygen and were usually in ridge areas where recharge typically is rapid. Anaerobic waters were in lower lying areas of Sumter County, where soils are poorly drained and aquifer recharge is slow. Anaerobic waters had higher concentrations of calcium, bicarbonate, sulfide, dissolved organic carbon, iron, manganese, and silica, and had lower concentrations of nitrate than aerobic ground waters. For low sulfate waters, sulfate generally originates from meteoric sources (atmospheric precipitation), with variable amounts of oxidation of reduced sulfur and sulfate reduction. Sulfide is sometimes removed from solution, probably by precipitation of a sulfide minerals such as pyrite. In areas where deep ground water has low sulfate concentrations, the shallow flow system is apparently deeper than where high sulfate concentrations occur, and upwelling sulfate-rich water is negligible. The range of sulfate concentrations observed in the study areas and differences in sulfate concentrations with depth indicate a complex interaction between shallow and deep ground-water flow systems.

Water-Resources Investigations Report↗
USGS1996

Assessment of intrinsic bioremediation of gasoline contamination in the shallow aquifer, Laurel Bay Exchange, Marine Corps Air Station Beaufort, South Carolina

Laboratory, field, and digital solute-transport- modeling studies demonstrate that microorganisms indigenous to the shallow ground-water system at Laurel Bay Exchange, Marine Corps Air Station Beaufort, South Carolina, can degrade petroleum hydrocarbons in gasoline released at the site. Microorganisms in aquifer sediments incubated in the laboratory under aerobic and anaerobic conditions mineralized radiolabeled carbon 14-toluene to 14 C-carbon dioxide with first-order rate constants of K bio = -0.640 per day and K bio = -0.003 per day, respectively. Digital solute- transport modeling using the numerical code SUTRA revealed that anaerobic biodegradation of benzene occurs with a first-order rate constant near K bio = -0.00025 per day. Sandy aquifer material beneath Laurel Bay Exchange is characterized by relatively high hydraulic conductivities (K aq = 8.9 to 17.3 feet per day), average ground-water flow rate of about 60 feet per year, and a relatively uniform hydraulic gradient of 0.004 feet per foot. The sandy aquifer material also has low adsorptive potentials for toluene and benzene (both about K ad = 2.0 x 10 -9 cubic feet per milligram), because of the lack of natural organic matter in the aquifer. The combination of this ground-water-flow rate and absence of significant adsorptive capacity in the aquifer permits toluene and benzene concentrations to be detected downgradient from the source area in monitoring wells, even though biodegradation of these compounds has been demonstrated. Solute-transport simulations, however, indicate that toluene and benzene will not reach the Broad River, the nearest point of contact with wildlife or human populations, about 3,600 feet west of the site boundary. These simulations also show that contamination will not be transported to the nearest Marine Corps property line about 2,400 feet south of the site. This is primarily because the source of contaminants has essentially been removed, and the low adsorptive capacity of the aquifer sediments has prevented the occurrence of an adsorbed, continuous source of petroleum hydrocarbons. Therefore, digital simulations of toluene and benzene transport at Laurel Bay Exchange indicate that intrinsic bioremediation could be a successful remediation alternative for prohibiting transport of dissolved toluene and benzene to the Broad River.

South Carolina↗
USGS2002

Simulation of Transient Ground-Water Flow in the Valley-Fill Aquifers of the Upper Rockaway River Basin, Morris County, New Jersey

More than 90 percent of the public water supply in the upper Rockaway River Valley in Morris County, New Jersey, is obtained from ground-water withdrawals from the valley-fill aquifers. During 1997, an average of 9.6 million gallons per day of ground water was withdrawn from these aquifers. The aquifer system consists of an unconfined aquifer (upper aquifer) and a locally confined aquifer (lower aquifer), which are composed of sands and gravels. These aquifers are separated by a discontinuous confining unit that consists mostly of silt and clay. Increases in ground-water withdrawals can induce movement of water from streams to wells, increase flow from the upper aquifer to the lower aquifer, and reduce base flow in the Rockaway River downstream. A ground-water-flow model was used to simulate and quantify the effects of current withdrawals on the valley-fill aquifer system under transient monthly conditions. Recharge over the model area varies both spatially and temporally. Part of model calibration consisted of adjusting percentages of monthly precipitation that recharges the valley-fill aquifer system. More recharge occurs during winter and spring than during summer and fall. This seasonal variation affects ground-water discharge to the Rockaway River. Ground-water withdrawals from the valleyfill aquifers also affect ground-water discharge to the Rockaway River. Three scenarios were simulated to observe the effects of ground-water withdrawals on ground-water discharge to the Rockaway River and to determine the extent to which variations in rates of withdrawals correspond to variations in rates of streamflow depletion. Streamflow depletion was estimated by comparing model-computed ground-water discharge for the three scenarios with the modelcomputed ground-water discharge under transient conditions. In scenario 1, all pumpage was removed from the model. In scenarios 2 and 3, 1 million gallons per day of ground-water withdrawals in excess of the current pumpage was withdrawn from the valley-fill aquifers. In scenario 2, the additional 1 million gallons per day of withdrawals were made from a hypothetical well located in the upper aquifer about 250 feet from the river. In scenario 3, the additional withdrawals were made from a hypothetical well located in the lower aquifer about 1,750 feet from the river. Results of scenario 1 indicated that the difference between the streamflow depletion and withdrawals is small; increases in ground-water withdrawals from the valley-fill aquifers correspond to decreases in ground-water discharge to the Rockaway River of approximately the same amount. Results of scenario 2 and 3 indicated that a lag time could occur between the introduction of withdrawals and the full magnitude of the effects of the withdrawals on streamflow depletion. A lag time of about seven months occurred for scenario 2 with the well placed in the upper aquifer. A longer lag time of more than 1.5 years occurred with the well placed in the lower aquifer and separated from the upper aquifer by a confining unit (scenario 3). Extreme low flow in the Rockaway River is mostly base flow. A flow-duration analysis of the Rockaway River at the surface-water gaging station upstream from the Boonton Reservoir during the drought of 1961-66 indicated that streamflow from the upper Rockaway River Basin alone might not be sufficient to meet the minimum passing flow of 7 million gallons per day during a drought. Under similar drought conditions today, during 3.2 percent of the drought time, streamflow at this station upstream from the reservoir would be less than the minimum passing flow requirement downstream from the reservoir.

Water-Resources Investigations Report↗
USGS2000

Investigation of the geology and hydrology of the Mogollon Highlands of central Arizona: a project of the Arizona Rural Watershed Initiative

The Mogollon Highlands of east central Arizona is a region of forested plateau and mountains, deep, sheerwalled canyons, and desert valleys. Known for its scenic beauty and characterized by a generally mild climate, the area, though still sparsely populated, attracts an increasing number of tourists and summer residents. Furthermore, the permanent population is expected to nearly double over the next 50 years. Consequently, there is increased pressure on the water resources of this area for several sometimes conflicting uses. Rational management of water resources is necessary to meet increased domestic requirements while ensuring an adequate supply of water for commercial and agricultural use, for Indian lands, and for preservation of valued environmental elements, including surface waters, riparian woodlands, forest and grassland areas, and wildlife and aquatic habitat. Such management requires an understanding of the relations among different components of the hydrologic system—recharge areas, surface flows, shallow aquifers, deep aquifers, discharge areas, and the regional ground-water flow system—and how each is affected by geology, climate, topography, and human use. The U.S. Geological Survey (USGS) is conducting an assessment of the hydrogeology of the Mogollon Highlands in cooperation with the Arizona Department of Water Resources. The study, funded through the State’s Rural Watershed Initiative program, is one of three assessments being conducted by the USGS. Assessments also are underway in the Upper-Middle Verde River watershed and on the Coconino Plateau. Each study has as its objectives: (1) the collection, compilation, and evaluation of all existing geologic, hydrologic, and related data pertaining to the study area and the creation of a data base that is readily accessible to the public and (2) developing an understanding of the hydrogeologic framework, which is the relation between geologic and hydrologic properties, that can be used for water-- resources management purposes and that will support the development of an interpretive and predictive model to estimate the effects of climate and water use on the sustainability of regional water resources. Although the three contiguous areas in north-central Arizona are being studied separately, a single data base is being constructed from which data on each area can be extracted separately.

Arizona↗
USGS2000

Computer-model analysis of ground-water flow and simulated effects of contaminant remediation at Naval Weapons Industrial Reserve Plant, Dallas, Texas

In June 1993, the Department of the Navy, Southern Division Naval Facilities Engineering Command (SOUTHDIV), began a Resource Conservation and Recovery Act (RCRA) Facility Investigation (RFI) of the Naval Weapons Industrial Reserve Plant (NWIRP) in north-central Texas. The RFI has found trichloroethene, dichloroethene, vinyl chloride, as well as chromium, lead, and other metallic residuum in the shallow alluvial aquifer underlying NWIRP. These findings and the possibility of on-site or off-site migration of contaminants prompted the need for a ground-water-flow model of the NWIRP area. The resulting U.S. Geological Survey (USGS) model: (1) defines aquifer properties, (2) computes water budgets, (3) delineates major flowpaths, and (4) simulates hydrologic effects of remediation activity. In addition to assisting with particle-tracking analyses, the calibrated model could support solute-transport modeling as well as help evaluate the effects of potential corrective action. The USGS model simulates steadystate and transient conditions of ground-water flow within a single model layer. The alluvial aquifer is within fluvial terrace deposits of Pleistocene age, which unconformably overlie the relatively impermeable Eagle Ford Shale of Late Cretaceous age. Over small distances and short periods, finer grained parts of the aquifer are separated hydraulically; however, most of the aquifer is connected circuitously through randomly distributed coarser grained sediments. The top of the underlying Eagle Ford Shale, a regional confining unit, is assumed to be the effective lower limit of ground-water circulation and chemical contamination. The calibrated steady-state model reproduces long-term average water levels within +5.1 or –3.5 feet of those observed; the standard error of the estimate is 1.07 feet with a mean residual of 0.02 foot. Hydraulic conductivity values range from 0.75 to 7.5 feet per day, and average about 4 feet per day. Specific yield values range from 0.005 to 0.15 and average about 0.08. Simulated infiltration rates range from 0 to 2.5 inches per year, depending mostly on local patterns of ground cover. Computer simulation indicates that, as of December 31, 1998, remediation systems at NWIRP were removing 7,375 cubic feet of water per day from the alluvial aquifer, with 3,050 cubic feet per day coming from aquifer storage. The resulting drawdown prevented 1,800 cubic feet per day of ground water from discharging into Cottonwood Bay, as well as inducing another 1,325 cubic feet per day into the aquifer from the bay. An additional 1,200 cubic feet of water per day (compared to pre-remediation conditions) was prevented from discharging into the west lagoon, east lagoon, Mountain Creek Lake, and Mountain Creek swale. Particle-tracking simulations, assuming an aquifer porosity of 0.15, were made to delineate flowpath patterns, or contaminant “capture zones,” resulting from 2.5- and 5-year periods of remediation activity at NWIRP. The resulting flowlines indicate three such zones, or areas from which ground water is simulated to have been removed during July 1996–December 1998, as well as extended areas from which ground water would be removed during the next 2.5 years (January 1999– June 2001). Simulation indicates that, as of December 31, 1998, the recovery trench was intercepting about 827 cubic feet per day of ground water that—without the trench—would have discharged into Cottonwood Bay. During this time, the trench is simulated to have removed about 3,221 cubic feet per day of water from the aquifer, with about 934 cubic feet per day (29 percent) coming from the south (Cottonwood Bay) side of the trench.

Texas↗
USGS2005

Questa baseline and premining ground-water quality investigation. 8. Lake-sediment geochemical record from 1960 to 2002, Eagle Rock and Fawn Lakes, Taos County, New Mexico

Geochemical studies of lake sediment from Eagle Rock Lake and upper Fawn Lake were conducted to evaluate the effect of mining at the Molycorp Questa porphyry molybdenum deposit located immediately north of the Red River. Two cores were taken, one from each lake near the outlet where the sediment was thinnest, and they were sampled at 1-cm intervals to provide geochemical data at less than 1-year resolution. Samples from the core intervals were digested and analyzed for 34 elements using ICP–AES (inductively coupled plasma–atomic emission spectrometry). The activity of 137 Cs has been used to establish the beginning of sedimentation in the two lakes. Correlation of the geochemistry of heavy-mineral suites in the cores from both Fawn and Eagle Rock Lakes has been used to develop a sedimentation model to date the intervals sampled. The core from upper Fawn Lake, located upstream of the deposit, provided an annual sedimentary record of the geochemical baseline for material being transported in the Red River, whereas the core from Eagle Rock Lake, located downstream of the deposit, provided an annual record of the effect of mining at the Questa mine on the sediment in the Red River. Abrupt changes in the concentrations of many lithophile and deposit-related metals occur in the middle of the Eagle Rock Lake core, which we correlate with the major flood-of-record recorded at the Questa gage at Eagle Rock Lake in 1979. Sediment from the Red River collected at low flow in 2002 is a poor match for the geochemical data from the sediment core in Eagle Rock Lake. The change in sediment geochemistry in Eagle Rock Lake in the post-1979 interval is dramatic and requires that a new source of sediment be identified that has substantially different geochemistry from that in the pre-1979 core interval. Loss of mill tailings from pipeline breaks are most likely responsible for some of the spikes in trace-element concentrations in the Eagle Rock Lake core. Enrichment of Al 2 O 3 , Cu, and Zn occurred as a result of chemical precipitation of these metals from ground water upstream in the Red River. Comparisons of the geochemistry of the post-1979 sediment core with both mine wastes and with premining sediment from the vicinity of the Questa mine indicate that both are possible sources for this new component of sediment. Existing data have not resolved this enigma.

New Mexico↗
USGS1993

Geology and ground-water resources in the Zebulon area, Georgia

The current (1991) surface-water source of drinking-water supply for the city of Zebulon, Pike County, Georgia, no longer provides an adequate water supply and periodically does not meet water-quality standards. The hydrogeology of crystalline rocks in the Zebulon area was evaluated to assess the potential of ground-water resources as a supplemental or alternative source of water to present surface-water supplies. As part of the ground-water resource evaluation, well location and construction data were compiled, a geologic map was constructed, and ground water was sampled and analyzed. Three mappable geologic units delineated during this study provide a basic understanding of hydrogeologic settings in the Zebulon area. Rock types include a variety of aluminosilicate schists, granitic rocks, amphibolites/honblende gneisses, and gondites. Several geologic features that may enhance ground-water availability were identified in the study area. These features include contacts between contrasting rock types, where a high degree of differential weathering has occurred, and well-developed structural features, such as foliation and jointing are present. High-yielding wells (greater than 25 gallons per minute) and low-yielding wells (less than one gallon per minute) were located in all three geologic units in a variety of topographic settings. Well yields range from less than one gallon per minute to 250 gallons per minute. The variable total depths and wide ranges of casing depths of the high-yielding wells are indicative of variations in depths to water-bearing zones and regolith thicknesses, respectively. The depth of water-bearing zones is highly variable, even on a local scale. Analyses of ground-water samples indicate that the distribution of iron concentration is as variable as well yield in the study area and does not seem to be related to a particular rock type. Iron concentrations in ground-water samples ranged from 0.02 to 5.3 milligrams per liter. Both iron concentration and well yield vary substantially over a relatively small area. Implementation and Verification of a One-Dimensional, Unsteady-Flow Model for Spring Brook near Warrenville, Illinois By Mary J. Turner, Anthony P. Pulokas, and Audrey L. Ishii Abstract A one-dimensional, unsteady-flow model, Full EQuations (FEQ) model, based on de Saint-Venant equations for dynamic flow in open channels, was calibrated and verified for a 0.75-mile reach of Spring Brook, a tributary to the West Branch Du Page River, near Warrenville in northeastern Illinois. The model was used to simulate streamflow in a small urban stream reach with two short culverts, one with overbank flow around the culvert during high flows. Streamflow data were collected on the reach during three high-flow periods. Data from one period were used to calibrate the model, and data from the other two periods were used to verify the model. Stages and discharges over the periods were simulated, and the results were compared graphically with stage and discharge data collected at 10 sites in the study reach. Errors in simulated stage and discharge were small except when debris, not represented in the model, clogged the culvert. The effects of changes in physical and computational model parameters also were studied. The model was insensit'lve to replacement of measured cross sections with interpolated cross sections, especially if the measured thalweg elevation was preserved. Variation of the roughness, slope, and length of the culvert over-bank section, as well as the chosen representative measured cross section, caused only slight changes in the simulated peak stage and discharge. Changes in the modeled culvert area caused large differences in the simulated highflows in the vicinity of the culvert, whereas simulated low flows were unaffected. At all flows, the misrepresentation of the culvert area caused the simulated water-surface elevations to deviate from the measured elevations, especially on the falling

Georgia↗
USGS2007

Estimated infiltration, percolation, and recharge rates at the Rillito Creek focused recharge investigation site, Pima County, Arizona

A large fraction of ground water stored in the alluvial aquifers in the Southwest is recharged by water that percolates through ephemeral stream-channel deposits. The amount of water currently recharging many of these aquifers is insufficient to meet current and future demands. Improving the understanding of streambed infiltration and the subsequent redistribution of water within the unsaturated zone is fundamental to quantifying and forming an accurate description of streambed recharge. In addition, improved estimates of recharge from ephemeral-stream channels will reduce uncertainties in water-budget components used in current ground-water models. This chapter presents a summary of findings related to a focused recharge investigation along Rillito Creek in Tucson, Arizona. A variety of approaches used to estimate infiltration, percolation, and recharge fluxes are presented that provide a wide range of temporal- and spatial-scale measurements of recharge beneath Rillito Creek. The approaches discussed include analyses of (1) cores and cuttings for hydraulic and textural properties, (2) environmental tracers from the water extracted from the cores and cuttings, (3) seepage measurements made during sustained streamflow, (4) heat as a tracer and numerical simulations of the movement of heat through the streambed sediments, (5) water-content variations, (6) water-level responses to streamflow in piezometers within the stream channel, and (7) gravity changes in response to recharge events. Hydraulic properties of the materials underlying Rillito Creek were used to estimate long-term potential recharge rates. Seepage measurements and analyses of temperature and water content were used to estimate infiltration rates, and environmental tracers were used to estimate percolation rates through the thick unsaturated zone. The presence or lack of tritium in the water was used to determine whether or not water in the unsaturated zone infiltrated within the past 40 years. Analysis of water-level and temporal-gravity data were used to estimate recharge volumes. Data presented in this chapter were collected from 1999 though 2002. Precipitation and streamflow during this period were less than the long-term average; however, two periods of significant streamflow resulted in recharge—one in the summer of 1999 and the other in the fall/winter of 2000. Flux estimates of infiltration and recharge vary from less than 0.1 to 1.0 cubic meter per second per kilometer of streamflow. Recharge-flux estimates are larger than infiltration estimates. Larger recharge fluxes than infiltration fluxes are explained by the scale of measurements. Methods used to estimate recharge rates incorporate the largest volumetric and temporal scales and are likely to have fluxes from other nearby sources, such as unmeasured tributaries, whereas the methods used to estimate infiltration incorporate the smallest scales, reflecting infiltration rates at individual measurement sites.

Arizona↗
USGS2006

Simulation of proposed increases in ground-water withdrawals on the Atlantic City 800-foot sand, New Jersey Coastal Plain

The confined Atlantic City 800-foot sand and the unconfined Kirkwood-Cohansey aquifer system (surficial aquifer) are major sources of water for southeastern New Jersey. Because of recent concerns about streamflow depletion resulting from ground-water withdrawals and the potential ecological effects on stream habitat in the area, the focus on future withdrawals has been shifted away from the surficial aquifer to the confined Atlantic City 800-foot sand until the effects of increased withdrawals from the surficial aquifer can be investigated. A study was conducted to evaluate the effects of seven proposed increases in ground-water withdrawals from the Atlantic City 800-foot sand and the Kirkwood-Cohansey aquifer system on the Atlantic City 800-foot sand. The proposed withdrawals are increases above the 2004 allocated rates (full allocation). The effects of full-allocation ground-water withdrawals and the cumulative effect of withdrawals for each of seven proposed increases in withdrawals were simulated using three previously published ground-water flow models: the New Jersey Coastal Plain Regional Aquifer System Analysis model, the Coastal Plain Optimization model, and a model of the Atlantic City 800-foot sand in Atlantic County, New Jersey. These models were used to simulate changes in water levels, the source supplying the increased ground-water flow, and the effects on saltwater movement towards production wells in Cape May County as a result of the proposed increased withdrawals at proposed or existing wells. The results of the simulations represent the effects of the proposed increase from full-allocation withdrawals to an additional 1,825 Mgal/yr (million gallons per year) from the Atlantic City 800-foot sand and an additional 1,045 Mgal/yr from the deep part of the Kirkwood-Cohansey aquifer system near the updip limit of the Atlantic City 800-foot sand. Most of the simulated decline in water levels in Atlantic County occurred as the result of the proposed increased withdrawals simulated for the New Jersey American Water Company wells. Simulated declines in water levels in Cape May were caused mainly by the simulated increased withdrawals for the Cape May City Desalination Plant wells. The additional water to supply the proposed increases in the scenarios was primarily horizontal flow from the unconfined updip part of the Kirkwood-Cohansey aquifer system, which accounted for 63 percent of the inflow, and flow from the overlying Kirkwood-Cohansey aquifer system into the Atlantic City 800-foot sand, which supplied 27 percent of the additional water. Because the withdrawals were made from the confined aquifer and the deeper part of the unconfined aquifer, the effect on streamflow was substantially less than would have occurred had the withdrawals been made directly from the shallower parts of the unconfined aquifer. The travel times from the 250-mg/L isochlor to production wells in Stone Harbor were longer as a result of all the additional withdrawals. For some scenarios, withdrawals in Atlantic County caused the saltwater to move slightly faster towards the production wells. These effects were offset by the increase in travel time caused by the potential increased withdrawals simulated for the Cape May City desalination wells, which either diverted water towards the desalination wells or increased the travel time towards production wells.

New Jersey↗
USGS1985

The ground-water system in the LaGrange Aquifer near LaGrange, southeastern Wyoming

Groundwater is being developed from the La Grange aquifer in southeastern Wyoming. It consists of saturated permeable alluvium that is hydraulically connected with most of the underlying White River Group. In the area of principal interest east of Horse Creek, Hawk Springs Reservoir and 14 adjacent wells used to supplement surface-water supply in the reservoir are in a natural discharge area. Upgradient of the reservoir there are 28 irrigation wells in about a 6-square-mile area. In this area, water levels declined between 3 and 12 feet from 1973 to 1978 causing concern about the effects of well pumpage on the hydrologic system. A digital model was developed and used to simulate the two-dimensional groundwater flow system in the unconfined La Grange aquifer. Transient simulations were made for 1973-78 using 12 time periods and for 1978-80 using 25 time periods. The calibrated digital model was used to simulate four 6-month pumping alternatives including three hypothetical alternatives for the area of principal interest east of Horse Creek. The reservoir altitude was held constant, approximating a reservoir volume 7,000 acre-feet. Pumping alternative 1 simulated historic conditions for 1973-78 including monthly recharge from precipitation which was included in the next three pumping alternatives. For pumping alternative 1, the calculated rate of discharge at the end of the 6-month simulation was 5.2 cubic feet per second from the aquifer to the reservoir. At the end of the 6-month simulations for pumping alternatives 2 and 3, the calculated rate of discharge to the reservoir was decreased to 0.4 cubic feet per second by pumpage from the 14 wells and to 3.8 cubic feet per second by pumpage from the 28 irrigation wells. For pumping alternative 4, pumpage from the total 42 wells resulted in a 1.0 cubic feet per second loss from the reservoir to the aquifer. (USGS)

Water-Resources Investigations Report↗
USGS2004

Simulated effects of ground-water augmentation on the hydrology of Round and Halfmoon Lakes in northwestern Hillsborough County, Florida

Pumpage from the Upper Floridan aquifer in northwest Hillsborough County near Tampa, Florida, has induced downward leakage from the overlying surficial aquifer and lowered the water table in many areas. Leakage is highest where the confining layer separating the aquifers is breached, which is common beneath many of the lakes in the study area. Leakage of water to the Upper Floridan aquifer has lowered the water level in many lakes and drained many wetlands. Ground water from the Upper Floridan aquifer has been added (augmented) to some lakes in an effort to maintain lake levels, but the resulting lake-water chemistry and lake leakage patterns are substantially different from those of natural lakes. Changes in lake-water chemistry can cause changes in lake flora, fauna, and lake sediment composition, and large volumes of lake leakage are suspected to enhance the formation of sinkholes near the shoreline of augmented lakes. The leakage rate of lake water through the surficial aquifer to the Upper Floridan aquifer was estimated in this study using ground-water-flow models developed for an augmented lake (Round Lake) and non-augmented lake (Halfmoon Lake). Flow models developed with MODFLOW were calibrated through nonlinear regression with UCODE to measured water levels and monthly net ground-water-flow rates from the lakes estimated from lake-water budgets. Monthly estimates of ground-water recharge were computed using an unsaturated flow model (LEACHM) that simulated daily changes in storage of water in the soil profile, thus estimating recharge as drainage to the water table. Aquifer properties in the Round Lake model were estimated through transient-state simulations using two sets of monthly recharge rates computed during July 1996 to February 1999, which spanned both average conditions (July 1996 through October 1997), and an El Niño event (November 1997 through September 1998) when the recharge rate doubled. Aquifer properties in the Halfmoon Lake model were estimated through steady-state simulations of average conditions in July 1996. Simulated hydrographs computed by the Round and Halfmoon Lake models closely matched measured water-level fluctuations, except during El Niño, when the Halfmoon Lake model was unable to accurately reproduce water levels. Possibly, potential recharge during El Niño was diverted through ground-water-flow outlets that were not represented in the Halfmoon Lake model, or a large part of the rainfall was diverted into runoff before it could become recharge. Solute transport simulations with MT3D indicate that leakage of lake water extended 250 to 400 feet into the surficial aquifer around Round Lake, and from 75 to 150 feet around Halfmoon Lake before flowing to the underlying Upper Floridan aquifer. These results are in agreement with concentrations of stable isotopes of oxygen-18 (d18O) and deuterium (dD) in the surficial aquifer. Schedules of monthly augmentation rates to maintain constant stages in Round and Halfmoon Lakes were computed using an equation that accounted for changes in the Upper Floridan aquifer head and the deviation from the mean recharge rate. Resulting lake stages were nearly constant during the first half of the study, but increased above target lake stages during El Niño; modifying the computation of augmentation rates to account for the higher recharge rate during El Niño resulted in lake stages that were closer to the target lake stage. Substantially more lake leakage flows to the Upper Floridan aquifer from Round Lake than from Halfmoon Lake, because the estimated vertical hydraulic conductivities of lake and confining layer sediments and breaches in the confining layer beneath Round Lake are much greater. Augmentation rates required to maintain the low guidance stages in Round Lake (53 feet) and Halfmoon Lake (42 feet) under average Upper Floridan aquifer heads are estimated as 33,850 cubic feet per day and 1,330 to 10,000 cubic feet per day, respectively. These rates equate to 26 inches per month of water applied to the entire surface of Round Lake and 0.34 to 2.5 inches per month of water applied to Halfmoon Lake.

Florida↗
USGS2017

A review of surface energy balance models for estimating actual evapotranspiration with remote sensing at high spatiotemporal resolution over large extents

Many approaches have been developed for measuring or estimating actual evapotranspiration ( ET a ), and research over many years has led to the development of remote sensing methods that are reliably reproducible and effective in estimating ET a . Several remote sensing methods can be used to estimate ET a at the high spatial resolution of agricultural fields and the large extent of river basins. More complex remote sensing methods apply an analytical approach to ET a estimation using physically based models of varied complexity that require a combination of ground-based and remote sensing data, and are grounded in the theory behind the surface energy balance model. This report, funded through cooperation with the International Joint Commission, provides an overview of selected remote sensing methods used for estimating water consumed through ET a and focuses on Mapping Evapotranspiration at High Resolution with Internalized Calibration (METRIC) and Operational Simplified Surface Energy Balance (SSEBop), two energy balance models for estimating ET a that are currently applied successfully in the United States. The METRIC model can produce maps of ET a at high spatial resolution (30 meters using Landsat data) for specific areas smaller than several hundred square kilometers in extent, an improvement in practice over methods used more generally at larger scales. Many studies validating METRIC estimates of ET a against measurements from lysimeters have shown model accuracies on daily to seasonal time scales ranging from 85 to 95 percent. The METRIC model is accurate, but the greater complexity of METRIC results in greater data requirements, and the internalized calibration of METRIC leads to greater skill required for implementation. In contrast, SSEBop is a simpler model, having reduced data requirements and greater ease of implementation without a substantial loss of accuracy in estimating ET a . The SSEBop model has been used to produce maps of ET a over very large extents (the conterminous United States) using lower spatial resolution (1 kilometer) Moderate Resolution Imaging Spectroradiometer (MODIS) data. Model accuracies ranging from 80 to 95 percent on daily to annual time scales have been shown in numerous studies that validated ET a estimates from SSEBop against eddy covariance measurements. The METRIC and SSEBop models can incorporate low and high spatial resolution data from MODIS and Landsat, but the high spatiotemporal resolution of ET a estimates using Landsat data over large extents takes immense computing power. Cloud computing is providing an opportunity for processing an increasing amount of geospatial “big data” in a decreasing period of time. For example, Google Earth Engine TM has been used to implement METRIC with automated calibration for regional-scale estimates of ET a using Landsat data. The U.S. Geological Survey also is using Google Earth Engine TM to implement SSEBop for estimating ET a in the United States at a continental scale using Landsat data.

Scientific Investigations Report↗
USGS1985

Hydrogeology, water quality, and ground-water development alternatives in the Beaver-Pasquiset ground-water reservoir, Rhode Island

In a 23 sq mi study area, the Beaver-Pasquiset groundwater reservoir within the Pawcatuck River basin in southern Rhode Island, stratified drift is the only principal geologic unit capable of producing yields > 350 gal/min. Transmissivity of the aquifer ranges from 7,200 to 24,300 sq ft/day. Water table conditions prevail in the aquifer, which is in good hydraulic connection with perennial streams and ponds. A digital model of two-dimensional groundwater flow was used to simulate the interaction between surface water and groundwater, and to evaluate the impact of alternative schemes of groundwater development on groundwater levels, pond levels, and streamflow in the Beaver-Pasquiset groundwater reservoir. Transient simulations of theoretical pumpage were made for a drought period (1963-66) and a wet period (1976-78). The areas most favorable for development of high-capacity wells (350 gal/min or more) are along the Beaver River and near Pasquiset Pond. The water is soft and generally contains < 100 mg/L dissolved solids. Locally, groundwater contains elevated concentrations of iron and manganese (7.5 and 3.7 mg/L, respectively), southeast of Pasquiset Pond, and will require treatment if used for public supply. The groundwater reservoir was simulated with a two-dimensional finite-difference model using a block-centered grid consisting of 33 rows and 75 columns. Differences between measured and simulated water table altitudes for the final steady state run for 21 selected observation wells averaged +0.07 ft. Combined pumping rates for simulation of groundwater development alternatives at eight sites ranged from 3.25 to 7.00 Mgal/d. Pumping rates for individual wells ranged from 0.25 to 1.50 Mgal/d. Transient simulations suggest that the Beaver-Pasquiset groundwater reservoir is capable of sustaining a pumping rate of 4.25 Mgal/d during years of average groundwater recharge with minimal impact on groundwater levels, pond levels, and streamflow. During extreme drought periods (1965 and 1966) it would be necessary to reduce pumpage below 3.25 Mgal/d to maintain flow in both the Beaver River and Pasquiset Brook. (Author 's abstract)

Water-Resources Investigations Report↗
USGS1988

Simulation of ground-water flow in the lower sand unit of the Potomac-Raritan-Magothy aquifer system, Philadelphia, Pennsylvania

Ground-water flow in the lower sand unit of the Potomac-Raritan-Magothy aquifer system in Philadelphia was simulated with a two-dimensional finite- difference ground-water model. The modeled 133-square-mile area also included parts of Delaware County, Pennsylvania, and Camden and Gloucester Counties, New Jersey. The lower sand unit is Cretaceous in age and consists of well- sorted coarse sand and fine gravel that grades upward into medium to fine sand containing a few thin beds of clay. The modeled aquifer consists of the lower sand unit in Philadelphia and the lowermost sand unit of the Potomac-Raritan-Magothy aquifer system in New Jersey. Throughout most of the area, the lower sand unit is overlain by a clay confining unit. Where the clay is absent, the lower sand unit is unconfined. A hydraulic conductivity of 1.6 x 10-3 foot per second and a storage coefficient of 3.0 x 10-4 was assigned to the lower sand unit based on 15 aquifer tests, and a hydraulic conductivity of 4.0 x 10-8 foot per second was assigned to the upper confining unit based on transient-flow sensitivity analysis. Water levels were not sensitive to changes in the value for specific storage of the upper confining unit, indicating that most vertical leakage occurs as steady leakage. Changes in the potentiometric surface of the lower sand unit for 1904-78 simulated. Differences between simulated and observed head generally were less than 10 feet. Simulations were made to determine the effects on hydraulic head of increases in industrial pumpage of 5 and 10 Mgal/d (million gallons per day) and of an emergency 60 Mgal/d municipal water supply in Philadelphia. A 5- and 10-Mgal/d increase in industrial pumpage would lower heads in the lower sand unit by as much as 33 and 66 feet, respectively. Pumping 60 Mgal/d for 30 days for an emergency municipal supply would lower heads in the lower sand unit by as much as 121 feet.

Water-Resources Investigations Report↗
USGS2005

Water budgets for selected watersheds in the Delaware River basin, eastern Pennsylvania and western New Jersey

This pilot study, done by the U.S. Geological Survey in cooperation with the Delaware River Basin Commission, developed annual water budgets using available data for five watersheds in the Delaware River Basin with different degrees of urbanization and different geological settings. A basin water budget and a water-use budget were developed for each watershed. The basin water budget describes inputs to the watershed (precipitation and imported water), outputs of water from the watershed (streamflow, exported water, leakage, consumed water, and evapotranspiration), and changes in ground-water and surface-water storage. The water-use budget describes water withdrawals in the watershed (ground-water and surface-water withdrawals), discharges of water in the watershed (discharge to surface water and ground water), and movement of water of water into and out of the watershed (imports, exports, and consumed water). The water-budget equations developed for this study can be applied to any watershed in the Delaware River Basin. Data used to develop the water budgets were obtained from available long-term meteorological and hydrological data-collection stations and from water-use data collected by regulatory agencies. In the Coastal Plain watersheds, net ground-water loss from unconfined to confined aquifers was determined by using ground-water-flow-model simulations. Error in the water-budget terms is caused by missing data, poor or incomplete measurements, overestimated or underestimated quantities, measurement or reporting errors, and the use of point measurements, such as precipitation and water levels, to estimate an areal quantity, particularly if the watershed is hydrologically or geologically complex or the data-collection station is outside the watershed. The complexity of the water budgets increases with increasing watershed urbanization and interbasin transfer of water. In the Wissahickon Creek watershed, for example, some ground water is discharged to streams in the watershed, some is exported as wastewater, and some is exported for public supply. In addition, ground water withdrawn outside the watershed is imported for public supply or imported as wastewater for treatment and discharge in the watershed. A GIS analysis was necessary to quantify many of the water-budget components. The 89.9-square mile East Branch Brandywine Creek watershed in Pennsylvania is a rural watershed with reservoir storage that is underlain by fractured rock. Water budgets were developed for 1977-2001. Average annual precipitation, streamflow, and evapotranspiration were 46.89, 21.58, and 25.88 inches, respectively. Some water was imported (average of 0.68 inches) into the watershed for public-water supply and as wastewater for treatment and discharge; these imports resulted in a net gain of water to the watershed. More water was discharged to East Branch Brandywine Creek than was withdrawn from it; the net discharge resulted in an increase in streamflow. Most ground water was withdrawn (average of 0.25 inches) for public-water supply. Surface water was withdrawn (average of 0.58 inches) for public-water and industrial supply. Discharge of water by sewage-treatment plants and industries (average of 1.22 inches) and regulation by Marsh Creek Reservoir caused base flow to appear an average of 7.2 percent higher than it would have been without these additional sources. On average, 67 percent of the difference was caused by sewage-treatment-plant and industrial discharges, and 33 percent was caused by regulation of the Marsh Creek Reservoir. Water imports, withdrawals, and discharges have been increasing as the watershed becomes increasingly urbanized. The 64-square mile Wissahickon Creek watershed in Pennsylvania is an urban watershed underlain by fractured rock. Water budgets were developed for 1987-98. Average annual precipitation, streamflow, and evapotranspiration were 47.23, 22.24, and 23.12 inches, respectively. The watershed is highly u

Pennsylvania↗
USGS1993

A modular finite-element model (MODFE) for areal and axisymmetric ground-water-flow problems, Part 3: Design philosophy and programming details

A MODular Finite-Element, digital-computer program (MODFE) was developed to simulate steady or unsteady-state, two-dimensional or axisymmetric ground-water-flow. The modular structure of MODFE places the computationally independent tasks that are performed routinely by digital-computer programs simulating ground-water flow into separate subroutines, which are executed from the main program by control statements. Each subroutine consists of complete sets of computations, or modules, which are identified by comment statements, and can be modified by the user without affecting unrelated computations elsewhere in the program. Simulation capabilities can be added or modified by either adding or modifying subroutines that perform specific computational tasks, and the modular-program structure allows the user to create versions of MODFE that contain only the simulation capabilities that pertain to the ground-water problem of interest. MODFE is written in a Fortran programming language that makes it virtually device independent and compatible with desk-top personal computers and large mainframes. MODFE uses computer storage and execution time efficiently by taking advantage of symmetry and sparseness within the coefficient matrices of the finite-element equations. Parts of the matrix coefficients are computed and stored as single-subscripted variables, which are assembled into a complete coefficient just prior to solution. Computer storage is reused during simulation to decrease storage requirements. Descriptions of subroutines that execute the computational steps of the modular-program structure are given in tables that cross reference the subroutines with particular versions of MODFE. Programming details of linear and nonlinear hydrologic terms are provided. Structure diagrams for the main programs show the order in which subroutines are executed for each version and illustrate some of the linear and nonlinear versions of MODFE that are possible. Computational aspects of changing stresses and boundary conditions with time and of mass-balance and error terms are given for each hydrologic feature. Program variables are listed and defined according to their occurrence in the main programs and in subroutines. Listings of the main programs and subroutines are given.

Techniques of Water-Resources Investigations↗