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 109 records · Page 6Linked to original sources
USGS1987

Hydrogeology and water quality of areas with persistent ground-water contamination near Blackfoot, Bingham County, Idaho

The Groveland-Collins area near Blackfoot, Idaho, has a history of either periodic or persistent localized groundwater contamination. Water users in the area report offensive smell, metallic taste, rust deposits, and bacteria in water supplies. During 1984 and 1985, data were collected to define regional and local geologic, hydrologic, and groundwater quality conditions, and to identify factors that may have affected local groundwater quality. Infiltration or leakage of irrigation water is the major source of groundwater recharge, and water levels may fluctuate 15 ft or more during the irrigation season. Groundwater movement is generally northwestward. Groundwater contains predominantly calcium, magnesium, and bicarbonate ions and characteristically has more than 200 mg/L hardness. Groundwater near the Groveland-Collins area may be contaminated from one or more sources, including infiltration of sewage effluent, gasoline or liquid fertilizer spillage, or land application of food processing wastewater. Subsurface basalt ridges impede lateral movement of water in localized areas. Groundwater pools temporarily behind these ridges and anomalously high water levels result. Maximum concentrations or values of constituents that indicate contamination were 1,450 microsiemens/cm specific conductance, 630 mg/L bicarbonate (as HCO 3 ), 11 mg/L nitrite plus nitrate (as nitrogen), 7.3 mg/L ammonia (as nitrogen), 5.9 mg/L organic nitrogen, 4.4 mg/L dissolved organic carbon, 7,000 micrograms/L dissolved iron, 5 ,100 microgram/L dissolved manganese, and 320 microgram/L dissolved zinc. Dissolved oxygen concentrations ranged from 8.9 mg/L in uncontaminated areas to 0 mg/L in areas where food processing wastewater is applied to the land surface. Stable-isotope may be useful in differentiating between contamination from potato-processing wastewater and whey in areas where both are applied to the land surface. Development of a ground-water model to evaluate effects of land applications of organic wastewater and organic solute loading rates on subsurface water quality is not feasible at this time.

Idaho↗
USGS2001

Documentation and verification of VST2D; a model for simulating transient, Variably Saturated, coupled water-heat-solute Transport in heterogeneous, anisotropic 2-Dimensional, ground-water systems with variable fluid density

This report describes a model for simulating transient, Variably Saturated, coupled water-heatsolute Transport in heterogeneous, anisotropic, 2-Dimensional, ground-water systems with variable fluid density (VST2D). VST2D was developed to help understand the effects of natural and anthropogenic factors on quantity and quality of variably saturated ground-water systems. The model solves simultaneously for one or more dependent variables (pressure, temperature, and concentration) at nodes in a horizontal or vertical mesh using a quasi-linearized general minimum residual method. This approach enhances computational speed beyond the speed of a sequential approach. Heterogeneous and anisotropic conditions are implemented locally using individual element property descriptions. This implementation allows local principal directions to differ among elements and from the global solution domain coordinates. Boundary conditions can include time-varying pressure head (or moisture content), heat, and/or concentration; fluxes distributed along domain boundaries and/or at internal node points; and/or convective moisture, heat, and solute fluxes along the domain boundaries; and/or unit hydraulic gradient along domain boundaries. Other model features include temperature and concentration dependent density (liquid and vapor) and viscosity, sorption and/or decay of a solute, and capability to determine moisture content beyond residual to zero. These features are described in the documentation together with development of the governing equations, application of the finite-element formulation (using the Galerkin approach), solution procedure, mass and energy balance considerations, input requirements, and output options. The VST2D model was verified, and results included solutions for problems of water transport under isohaline and isothermal conditions, heat transport under isobaric and isohaline conditions, solute transport under isobaric and isothermal conditions, and coupled water-heat-solute transport. The first three problems considered in model verification were compared to either analytical or numerical solutions, whereas the coupled problem was compared to measured laboratory results for which no known analytic solutions or numerical models are available. The test results indicate the model is accurate and applicable for a wide range of conditions, including when water (liquid and vapor), heat (sensible and latent), and solute are coupled in ground-water systems. The cumulative residual errors for the coupled problem tested was less than 10–8 cubic centimeter per cubic centimeter, 10-5 moles per kilogram, and 102 calories per cubic meter for liquid water content, solute concentration and heat content, respectively. This model should be useful to hydrologists, engineers, and researchers interested in studying coupled processes associated with variably saturated transport in ground-water systems.

Water-Resources Investigations Report↗
USGS2003

Development and calibration of a ground-water flow model for the Sparta Aquifer of southeastern Arkansas and north-central Louisiana and simulated response to withdrawals, 1998-2027

The Sparta aquifer, which consists of the Sparta Sand, in southeastern Arkansas and north-central Louisiana is a major water resource and provides water for municipal, industrial, and agricultural uses. In recent years, the demand in some areas has resulted in withdrawals from the Sparta aquifer that substantially exceed replenishment of the aquifer. Considerable drawdown has occurred in the potentiometric surface forming regional cones of depression as water is removed from storage by withdrawals. These cones of depression are centered beneath the Grand Prairie area and the cities of Pine Bluff and El Dorado in Arkansas, and Monroe in Louisiana. The rate of decline for hydraulic heads in the aquifer has been greater than 1 foot per year for more than a decade in much of southern Arkansas and northern Louisiana where hydraulic heads are now below the top of the Sparta Sand. Continued hydraulic-head declines have caused water users and managers alike to question the ability of the aquifer to supply water for the long term. Concern over protecting the Sparta aquifer as a sustainable resource has resulted in a continued, cooperative effort by the Arkansas Soil and Water Conservation Commission, U.S. Army Corps of Engineers, and the U.S. Geological Survey to develop, maintain, and utilize numerical ground-water flow models to manage and further analyze the ground-water system. The work presented in this report describes the development and calibration of a ground-water flow model representing the Sparta aquifer to simulate observed hydraulic heads, documents major differences in the current Sparta model compared to the previous Sparta model calibrated in the mid-1980's, and presents the results of three hypothetical future withdrawal scenarios. The current Sparta model-a regional scale, three-dimensional numerical ground-water flow model-was constructed and calibrated using available hydrogeologic, hydraulic, and water-use data from 1898 to 1997. Significant changes from the previous model include grid rediscretization of the aquifer, extension of the active model area northward beyond the Cane River Formation facies change, and representation of model boundaries. The current model was calibrated with the aid of parameter estimation, a nonlinear regression technique, combined with trial and error parameter adjustment using a total of 795 observations from 316 wells over 4 different years-1970, 1985, 1990, and 1997. The calibration data set provides broad spatial and temporal coverage of aquifer conditions. Analysis of the residual statistics, spatial distribution of residuals, simulated compared to observed hydrographs, and simulated compared to observed potentiometric surfaces were used to analyze the ability of the calibrated model to simulate aquifer conditions within acceptable error. The calibrated model has a root mean square error of 18 feet for all observations, an improvement of more than 12 feet from the previous model. The current Sparta model was used to predict the effects of three hypothetical withdrawal scenarios on hydraulic heads over the period 1998-2027 with one of those extended indefinitely until equilibrium conditions were attained, or steady state. In scenario 1a, withdrawals representing the time period from 1990 to 1997 was held constant for 30 years from 1998 to 2027. Hydraulic heads in the middle of the cone of depression centered on El Dorado decreased by 10 feet from the 1997 simulation to 222 feet below NGVD of 1929 in 2027. Hydraulic heads in the Pine Bluff cone of depression showed a greater decline from 61 feet below NGVD of 1929 to 78 feet below NGVD of 1929 in the center of the cone. With these same withdrawals extended to steady state (scenario 1b), hydraulic heads in the Pine Bluff cone of depression center declined an 2 Development and Calibration of a Ground-Water Flow Model for the Sparta Aquifer of Southeastern Arkansas and North-Central Louisiana and Simulated Response to Withdrawa

Water-Resources Investigations Report↗
USGS2005

Simulated effects of water-level changes in the Mississippi River and Pokegama Reservoir on ground-water levels, Grand Rapids area, Minnesota

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, used an existing, three-dimensional, numerical ground-water flow model (referred to as the calibrated model) to assess the effects of water-level changes in the Mississippi River and Pokegama Reservoir on ground-water levels in adjacent glaciofluvial aquifers in the Grand Rapids area of north-central Minnesota. Pokegama Reservoir consists of Pokegama Lake, Little Jay Gould Lake, Jay Gould Lake, Cut-off Lake, and Blackwater Lake. Water levels in the Pokegama Reservoir are regulated at Pokegama Dam on the Mississippi River west of Grand Rapids. A steady-state model was used, and simulations represent “worse-case” scenarios for the effects of lowering or raising the river and lake water levels. The simulated ground-water levels represent levels that would result if the river and lake stages permanently declined or rose by the specified amounts. Eight hypothetical scenarios were simulated by varying water levels in the Mississippi River and Pokegama Reservoir from values used in the calibrated model. In the simulations, water levels for the Mississippi River, riverine wetlands of the Mississippi River, and lakes of the Pokegama Reservoir were raised and lowered uniformly by 0.50, 1.00, 2.00, and 3.00 feet from calibrated water levels. The extent of aquifer water-level changes resulting from these river, wetland, and lake water-level changes varied because of the complex hydrogeology of the study area. A 1.00-foot decline in reservoir/river water levels caused a maximum simulated ground-water-level decline in the middle aquifer near Jay Gould and Little Jay Gould Lakes of 1.09 feet and a maximum simulated ground-water-level decline of 1.00 foot in the lower aquifer near Cut-off and Blackwater Lakes. The amount and extent of ground-water-level changes in the middle and lower aquifers can be explained by the thickness, extent, and connectivity of the aquifers. Surface-water/ground-water interactions near wetlands and lakes with water levels unchanged from the calibrated model resulted in small water-table altitude differences among the simulations. Results of the ground-water modeling indicate that lowering of the reservoir and river water levels by 1.00 foot likely will not substantially affect water levels in the middle and lower aquifers.

Minnesota↗
USGS2008

Availability, sustainability, and suitability of ground water, Rogers Mesa, Delta County, Colorado: Types of analyses and data for use in subdivision water-supply reports

The population of Delta County, Colorado, like that in much of the Western United States, is forecast to increase substantially in the next few decades. A substantial portion of the increased population likely will reside in rural subdivisions and use residential wells for domestic water supplies. In Colorado, a subdivision developer is required to submit a water-supply plan through the county for approval by the Colorado Division of Water Resources. If the water supply is to be provided by wells, the water-supply plan must include a water-supply report. The water-supply report demonstrates the availability, sustainability, and suitability of the water supply for the proposed subdivision. During 2006, the U.S. Geological Survey, in cooperation with Delta County, Colorado, began a study to develop criteria that the Delta County Land Use Department can use to evaluate water-supply reports for proposed subdivisions. A table was prepared that lists the types of analyses and data that may be needed in a water-supply report for a water-supply plan that proposes the use of ground water. A preliminary analysis of the availability, sustainability, and suitability of the ground-water resources of Rogers Mesa, Delta County, Colorado, was prepared for a hypothetical subdivision to demonstrate hydrologic analyses and data that may be needed for water-supply reports for proposed subdivisions. Rogers Mesa is a 12-square-mile upland mesa located along the north side of the North Fork Gunnison River about 15 miles east of Delta, Colorado. The principal land use on Rogers Mesa is irrigated agriculture, with about 5,651 acres of irrigated cropland, grass pasture, and orchards. The principal source of irrigation water is surface water diverted from the North Fork Gunnison River and Leroux Creek. The estimated area of platted subdivisions on or partially on Rogers Mesa in 2007 was about 4,792 acres of which about 2,756 acres was irrigated land in 2000. The principal aquifer on Rogers Mesa consists of alluvial-fan deposits that overlie shale and, locally, sandstone. Maps of the base of the aquifer, the water table, and the saturated thickness of the aquifer were prepared from data from the well files of the Colorado Division of Water Resources. The base of the aquifer generally is topographically higher than the valleys of the North Fork Gunnison River and Leroux Creek, and direct hydraulic connection of the aquifer to North Fork Gunnison River and Leroux Creek is limited. The aquifer is recharged primarily by infiltration of surface water diverted for irrigation. Ground water discharges to seeps and springs and through slope deposits at the boundaries of the aquifer. Data from the well files also were used to estimate the specific capacity of wells and to estimate the transmissivity and hydraulic conductivity of the aquifer. A water budget was used to estimate recharge to and discharge from the aquifer. Although storage within the aquifer likely varies seasonally and from year to year, it was assumed that there were no long-term changes in ground-water storage. Estimated average annual recharge to and discharge from the aquifer during November 1998 through October 2006 were about 30,767 acre-feet per year. Although sufficient ground water is available on Rogers Mesa for additional domestic water supplies, conversion of irrigated land to residential land use likely would reduce recharge to the aquifer, affecting the sustainability of ground-water supplies on Rogers Mesa. Stream-depletion analyses indicate that the ground water in the aquifer likely would be considered tributary ground water and additional uses of ground water to supply new subdivisions likely would require implementation of augmentation plans. Although sufficient ground water is available on Rogers Mesa for additional domestic water supplies, conversion of irrigated land to residential land use likely would reduce recharge to the aquifer, affecting the sustainability of ground-water supplies on Rogers Mesa. Stream-depletion analyses indicate that the ground water in the aquifer likely would be considered tributary ground water and additional uses of ground water to supply new subdivisions likely would require implementation of augmentation plans. Although the dissolved solids and dissolved sulfate concentrations in ground water from Rogers Mesa aquifer commonly exceeded the U.S. Environmental Protection Agency Secondary Maximum Contaminant Levels for drinking-water supplies, the quality of ground water from the aquifer generally is suitable for residential use. Concentrations of total nitrogen (nitrite plus nitrate, as nitrogen) in ground water ranged from 0.38 to 3.2 milligrams per liter and were less than the State of Colorado maximum contaminant level of 10 milligrams per liter. Concentrations of selenium from seeps and springs at the boundaries of the aquifer commonly exceeded 50 micrograms per liter, the State of Colorado maximum contaminant level for drinking-water supplies. This preliminary evaluation of ground-water supplies on Rogers Mesa could be improved with the collection of additional data including: additional mapping of hydrogeologic features; more accurate locations and altitudes of wells; accurate estimates of water-budget components; measurements of ground-water levels; and collection and analyses of ground-water samples. The use of numerical models of ground-water flow could improve evaluations of the potential effects of changes in land and water use on the water budget, aquifer storage, stream depletion, and well interference.

Colorado↗
USGS1981

Hydrogeologic setting of the Glacial Lake Agassiz Peatlands, northern Minnesota

Seven test holes drilled in the Glacial Lake Agassiz Peatlands indicate that the thickness of surficial materials along a north-south traverse parallel to Minnesota Highway 72 ranges from 163 feet near Blackduck, Minnesota to 57 feet about 3 miles south of Upper Red Lake. Lenses of sand and gravel occur immediately above bedrock on the Itasca moraine and are interbedded with lake clay and till under the peatlands. Vertical head gradients measured in a piezometer nest near Blackduck on the moraine are downward, indicative of recharge to the regional ground-water-flow system. Vertical head gradients are upward in a piezometer nest on a sand beach ridge in the peatlands 12 miles north of Upper Red Lake. Numerical sectional models indicate that this discharge probably comes from local flow systems recharged from ground-water mounds located under large raised bogs.

Minnesota↗
USGS2006

Factors affecting occurrence and distribution of selected contaminants in ground water from selected areas in the Piedmont Aquifer System, Eastern United States, 1993-2003

Results of ground-water sampling from 255 wells and 19 springs in 11 studies done by the U.S. Geological Survey National Water-Quality Assessment (NAWQA) Program within the Piedmont Aquifer System (PAS) were analyzed to determine the factors affecting occurrence and distribution of selected contaminants. The contaminants, which were selected on the basis of potential human-health effects, included nitrate, pesticides, volatile organic compounds (VOCs), and radon. The PAS was subdivided on the basis of the general rock type of the aquifers into three areas for the study—crystalline, carbonate, and siliciclastic. The 11 studies were designed to areally represent an individual aquifer rock type and overall are representative of the PAS in their distribution; 7 studies are in the crystalline-rock aquifers, 3 studies are in the siliciclasticrock aquifers, and 1 study is in the carbonate-rock aquifers. Four of the studies were focused on land use, 1 in an agricultural area and 3 in urban areas. The remaining studies had wells representing a range of land-use types. Analysis of results of nitrate sampling indicated that in 8 of the 10 areas where nitrate concentrations were measured, median concentrations of nitrate were below 3 mg/L (milligrams per liter); 2 of the 10 areas had statistically significant higher median concentrations when compared to the other 8 areas. The agricultural land-use study in the carbonate-rock aquifer in the Lower Susquehanna River Basin had the highest median nitrate concentration (11 mg/L), and 60 percent of the wells sampled exceeded the U.S. Environmental Protection Agency (USEPA) Maximum Contaminant Level (MCL) of 10 mg/L. The major aquifer study in the crystalline-rock aquifer of the Lower Susquehanna River Basin Study Unit had the second-highest median nitrate concentration. Nitrate concentrations were positively correlated to the percentage of agricultural land use around the well, the total input of nitrogen from all sources, dissolved oxygen concentration, lithology, depth to water, and soil-matrix characteristics. A linear regression model was used to determine that increases in the percentage of agricultural land use, the input of nitrogen from all sources, and dissolved oxygen were the most significant variables affecting increased concentration of nitrate. A logistic regression model was used to determine that those same factors were the most significant variables affecting whether or not the nitrate concentration would exceed 4 mg/L. Of the analysis of samples from 253 wells and 19 springs for 47 pesticides, no sample had a pesticide concentration that exceeded any USEPA MCL. The most frequently detected pesticide was desethyl atrazine, a degradation product of atrazine; the detection frequency was 47 percent. Other frequently detected pesticides included atrazine, metolachlor, simazine, alachlor, prometon, and dieldrin. Detection frequency was affected by the analytical reporting limits; the frequency of detection was somewhat lower when all pesticides were censored to the highest common detection limit. Source factors such as agricultural land use (for agricultural herbicides), urban land use (for insecticides), and the application rate were found to have positive statistical correlations with pesticide concentration. Transport factors such as depth to water and percentage of well-drained soils, sand, or silt typically were positively correlated with higher pesticide concentrations. Sampling for VOCs was conducted in 187 wells and 19 springs that were sampled for 59 VOCs. There were 137 detections of VOCs above the common censoring limit of 0.2 µg/L. The most frequently detected VOCs were chloroform, a trihalomethane, and methyl-tert butyl ether (MTBE), a fuel oxygenate. Seventy-nine wells had at least one VOC detected. The detections were related to land use and well depth. Kendall’s tau correlations indicated a significant positive correlation between chloroform concentration and urban land use, leaking underground storage tanks, population density, and well depth. MTBE concentrations also were positively correlated to urban land use, leaking underground storage tanks, population density, and well depth. Radon was sampled at 205 sites. The subdivisions used for analysis of other contaminants were not adequate for analysis of radon because radon varies on the basis of variations in mineralogy that are not reflected by the general lithologic categories used for the rest of the studies. Concentrations of radon were highest in areas where the crystalline-rock aquifers had felsic mineralogy, and the lowest concentrations of radon were in areas where the crystalline-rocks aquifer had mafic mineralogy. Water from wells in siliciclastic-rock aquifers had concentrations of radon lower than that in the felsic crystalline-rock aquifers. More than 90 percent of the wells sampled for radon exceeded the proposed MCL of 300 pCi/L (picoCuries per liter); however, only 13 percent of those wells had concentrations in water that exceeded the alternative maximum contaminant level (AMCL), a higher level that can be used by municipalities addressing other sources of radon exposure. Overall, concentrations of constituents were related to land-use factors for nitrate, pesticides, VOCs, and to aquifer lithology for radon. None of the 47 pesticides or 59 VOCs analyzed exceeded the MCLs where those constituents were sampled. Concentrations exceeded the MCL for nitrate in 11 percent of the wells sampled. Nearly 91 percent of the wells sampled exceeded the proposed MCL for radon. Additional sampling in selected areas would improve overall understanding of the PAS and increase the possibility of creating predictive models of ground-water quality in this area.

Alabama, Delaware, Georgia, Maryland, New Jersey, ↗
USGS2011

Hydrogeology and simulation of groundwater flow in fractured rock in the Newark basin, Rockland County, New York

Groundwater in the Newark basin aquifer flows primarily through discrete water-bearing zones parallel to the strike and dip of bedding, whereas flow perpendicular to the strike is restricted, thereby imparting anisotropy to the groundwater flow field. The finite-element model SUTRA was used to represent bedrock structure in the aquifer by spatially varying the orientation of the hydraulic conductivity tensor to reflect variations in the strike and dip of the bedding. Directions of maximum and medium hydraulic conductivity were oriented parallel to the bedding, and the direction of minimum hydraulic conductivity was oriented perpendicular to the bedding. Groundwater flow models were prepared to simulate local flow in the vicinity of the Spring Valley well field and regional flow through the Newark basin aquifer. The Newark basin contains sedimentary rocks deposited as alluvium during the Late Triassic and is one of a series of basins that developed when Mesozoic rifting of the super continent Pangea created the Atlantic Ocean. The westward-dipping basin is filled with interbedded facies of coarse-grained to fine-grained rocks that were intruded by diabase associated with Jurassic volcanism. The Newark basin aquifer is bounded to the north and east by the Palisades sill and to the west by the Ramapo Fault. Although the general dip of bedding is toward the fault, mapping of conglomerate beds indicates the rocks are folded into broad anticlines and synclines. An alternative, more uniform pattern of regional structure, based on interpolated strike and dip measurements from a number of sources, has also been proposed. Two groundwater flow models (A for the former type of bedrock structure and B for the latter type) were developed to represent these alternative depictions of bedrock structure. Transient simulations were calibrated to reproduce measured water-level recoveries in a 9.3 mi&sup2 area surrounding the Spring Valley well field during a 5-day aquifer test in 1992. The models represented a 330-ft thick rock mass divided vertically into 10 equally spaced layers and were calibrated through nonlinear regression. Results of model B best matched the observed water-level recoveries with an estimated hydraulic conductivity of 9.5 ft/day, specific storage of 7.6 x 10 -6 ft -1, and K max : K min anisotropy ratio (hydraulic conductivity parallel to bedding: perpendicular to bedding) of 72:1. Model error was 50 percent greater in model A because the assumed structure did not match the actual strike of bedding in this area. Steady-state simulations of regional flow through the 85.4-mi2 modeled extent of the Newark basin aquifer represented both the alluvial aquifer beneath the Mawah River and the fractured bedrock. The rock mass was divided into two aquifer units: an upper 500-ft thick unit divided into 10 equally spaced layers through which most groundwater is assumed to flow and a lower unit divided into 7 layers with increasing thickness. Models were calibrated through nonlinear regression to average water levels measured in 140 wells from August 2005 through April 2007. Water levels simulated using the two models were similar and generally matched those observed, and the average recharge rate estimated using both models was 19 inches/year for the simulated period. Estimated transmissivity parallel to the strike of bedding (1,100 ft&sup2/d) was uniform in two transmissivity (T) zones in model A, but in model B the transmissivity of a high T zone (1,600 ft&sup2/d), delineated on the basis of aquifer test data, was slightly greater than in a low T zone (1,300 ft&sup2/d). The K max : K min anisotropy was estimated to be 58:1 in model A and 410:1 in model B, so the proportion of flow perpendicular to bedding is less in model B than in model A. Distributions of groundwater age simulated with models A and B are similar and indicate that most shallow ground-water (225 ft below the bedrock surface) is 5 t

Scientific Investigations Report↗
USGS2004

Regional ground-water-flow models of surficial sand and gravel aquifers along the Mississippi River between Brainerd and St. Cloud, central Minnesota

This report documents regional ground-waterflow models constructed by the U.S. Geological Survey in cooperation with the Minnesota Department of Health (MDH) to satisfy the requirements of their Source Water Protection Plan (SWPP). Steady-state single-layer ground-water-flow models were constructed with the computer program MODFLOW to simulate flow in surficial sand and gravel aquifers along the Mississippi River between Brainerd and St. Cloud in central Minnesota. The hydrogeologic data that were used to construct the models were compiled from available sources. Calibrated values of horizontal hydraulic conductivity and areal recharge for the aquifer in a northern model area were 70 m/d and 3.0x10-4 m/d, respectively. This model was sensitive to net areal recharge, vertical hydraulic conductivity of perennial streambed sediments, and horizontal hydraulic conductivity. The major source of net inflow to the model was from edge boundary cells. The major source of net outflow was ground-water discharge to perennial and ephemeral streams. Calibrated values of horizontal hydraulic conductivity and areal recharge for the aquifer in a southern model area were 70 m/d and 6.0x10-4 m/d, respectively. This model was sensitive mostly to horizontal hydraulic conductivity. Net areal recharge and ground-water discharge to perennial streams were the major sources of net inflow and outflow, respectively.

Minnesota↗
USGS2004

Geohydrology of the French Creek Basin and simulated effects of drought and ground-water withdrawals, Chester County, Pennsylvania

This report describes the results of a study by the U.S. Geological Survey, in cooperation with the Delaware River Basin Commission, to develop a regional ground-water-flow model of the French Creek Basin in Chester County, Pa. The model was used to assist water-resource managers by illustrating the interconnection between ground-water and surface-water systems. The 70.7-square mile French Creek Basin is in the Piedmont Physiographic Province and is underlain by crystalline and sedimentary fractured-rock aquifers. Annual water budgets were calculated for 1969-2001 for the French Creek Basin upstream of streamflow-measurement station French Creek near Phoenixville (01472157). Average annual precipitation was 46.28 in. (inches), average annual streamflow was 20.29 in., average annual base flow determined by hydrograph separation was 12.42 in., and estimated average annual ET (evapotranspiration) was 26.10 in. Estimated average annual recharge was 14.32 in. and is equal to 31 percent of the average annual precipitation. Base flow made up an average of 61 percent of streamflow. Ground-water flow in the French Creek Basin was simulated using the finite-difference MODFLOW-96 computer program. The model structure is based on a simplified two-dimensional conceptualization of the ground-water-flow system. The modeled area was extended outside the French Creek Basin to natural hydrologic boundaries; the modeled area includes 40 square miles of adjacent areas outside the basin. The hydraulic conductivity for each geologic unit was calculated from reported specific-capacity data determined from aquifer tests and was adjusted during model calibration. The model was calibrated for above-average conditions by simulating base-flow and water-level measurements made on May 1, 2001, using a recharge rate of 20 in/yr (inches per year). The model was calibrated for below-average conditions by simulating base-flow and water-level measurements made on September 11 and 17, 2001, using a recharge rate of 6.2 in/yr. Average conditions were simulated by adjusting the recharge rate until simulated streamflow at streamflow-measurement station 01472157 matched the long-term (1968-2001) average base flow of 54.1 cubic feet per second. The recharge rate used for average conditions was 15.7 in/yr. The effect of drought in the French Creek Basin was simulated using a drought year recharge rate of 8 in/yr for 3 months. After 3 months of drought, the simulated streamflow of French Creek at streamflow-measurement station 01472157 decreased 34 percent. The simulations show that after 6 months of average recharge (15.7 in/yr) following drought, streamflow and water levels recovered almost to pre-drought conditions. The effect of increased ground-water withdrawals on stream base flow in the South Branch French Creek Subbasin was simulated under average and drought conditions with pumping rates equal to 50, 75, and 100 percent of the Delaware River Basin Commission Ground Water Protected Area (GWPA) withdrawal limit (1,393 million gallons per year) with all pumped water removed from the basin. For average recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 18, 28, and 37 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. After 3 months of drought recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 27, 40, and 52 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. The effect of well location on base flow, water levels, and the sources of water to the well was simulated by locating a hypothetical well pumping 200 gallons per minute in different places in the Beaver Run Subbasin with all pumped water removed from the basin. The smallest reduction in the base flow of Beaver Run was from a well on the drainage divide between the French Creek Basin and the Marsh Creek Basin to the south; the simulated base flow of Beaver Run at the mouth was reduced 1 percent. The greatest reduction in the base flow of Beaver Creek was from a well close to Beaver Run; the simulated base flow of Beaver Run at the mouth was reduced 8 percent. The simulations showed that (1) if the contributing area of a well is in a basin, pumping will affect stream base flow and water levels in that basin whether the well is inside or outside that basin; (2) wells in different areas of a basin away from a divide produce a similar reduction in base flow; (3) a well within a basin will derive more water from diverted base flow and less water from storage than a well on or near a basin divide; and (4) the reduction in base flow at the mouth of the stream is the same for a well in the headwaters and a well downstream near the confluence. Model simulations illustrate some of the typical analyses and results that can be produced. The model was calibrated using annual values for recharge and ground-water ET and then was run using the annual values in a seasonally independent transient mode to show changes with time. The timing and relative magnitude of some of the changes simulated with the model when viewed in terms of a normal climatic year may be subject to considerable uncertainty because of the variability in seasonal recharge and ground-water ET rates. Transient model simulations for short-term periods are indicative of possible hydrologic system response and are considered an approximation.

Pennsyvlania↗
USGS2006

Water Resources Investigations at Edwards Air Force Base since 1988

Edwards Air Force Base (EAFB) in southern California (fig. 1) has relied on ground water to meet its water-supply needs. The extraction of ground water has led to two major problems that can directly affect the mission of EAFB: declining water levels (more than 120 ft since the 1920s) and land subsidence, a gradual downward movement of the land surface (more than 4 ft since the late 1920s). As water levels decline, this valuable resource becomes depleted, thus requiring mitigating measures. Land subsidence has caused cracked (fissured) runways and accelerated erosion on Rogers lakebed. In 1988, the U.S. Geological Survey (USGS), in cooperation with the U.S. Air Force, began investigations of the effects of declining water levels and land subsidence at EAFB and possible mitigation measures, such as the injection of imported surface water into the ground-water system. The cooperative investigations included data collection and analyses, numerical simulations of ground-water flow and land subsidence, and development of a preliminary simulation-optimization model. The results of these investigations indicate that the injection of imported water may help to control land subsidence; however, the potential ground-water-quality impacts are unknown.

Fact Sheet↗
USGS2005

Simulation of ground-water flow in coastal Georgia and adjacent parts of South Carolina and Florida-predevelopment, 1980, and 2000

A digital model was developed to simulate steady-state ground-water flow in a 42,155-square-mile area of coastal Georgia and adjacent parts of South Carolina and Florida. The model was developed to (1) understand and refine the conceptual model of regional ground-water flow, (2) serve as a framework for the development of digital subregional ground-water flow and solute-transport models, and (3) serve as a tool for future evaluations of hypothetical pumping scenarios used to facilitate water management in the coastal area. Single-density ground-water flow was simulated using the U.S. Geological Survey finite-difference code MODFLOW-2000 for mean-annual conditions during predevelopment (pre?1900) and the years 1980 and 2000. The model comprises seven layers: the surficial aquifer system, the Brunswick aquifer system, the Upper Floridan aquifer, the Lower Floridan aquifer, and the intervening confining units. A combination of boundary conditions was applied, including a general-head boundary condition on the top active cells of the model and a time-variable fixed-head boundary condition along part of the southern lateral boundary. Simulated heads for 1980 and 2000 conditions indicate a good match to observed values, based on a plus-or-minus 10-foot (ft) calibration target and calibration statistics. The root-mean square of residual water levels for the Upper Floridan aquifer was 13.0 ft for the 1980 calibration and 9.94 ft for the 2000 calibration. Some spatial patterns of residuals were indicated for the 1980 and 2000 simulations, and are likely a result of model-grid cell size and insufficiently detailed hydraulic-property and pumpage data in some areas. Simulated potentiometric surfaces for predevelopment, 1980, and 2000 conditions all show major flow system features that are indicated by estimated peotentiometric maps. During 1980?2000, simulated water levels at the centers of pumping at Savannah and Brunswick rose more than 20 ft and 8 ft, respectively, in response to decreased pumping. Simulated drawdown exceeded 10 ft in the Upper Floridan aquifer across much of the western half of the model area, with drawdown exceeding 20 ft along parts of the western, northern, and southern boundaries where irrigation pumping increased during this period. From predevelopment to 2000 conditions, the simulated water budget showed an increase in inflow from, and decrease in outflow to, the general-head boundaries, and a reversal from net seaward flow to net landward flow across the coastline. Simulated changes in recharge and discharge distribution from predevelopment to 2000 conditions showed an increase in extent and magnitude of net recharge cells in the northern part of the model area, and a decrease in discharge or change to recharge in cells containing major streams and beneath major pumping centers. The model is relatively sensitive to pumping and the controlling head at the fixed-head boundary and less sensitive to the distribution of aquifer properties in general. Model limitations include: (1) its spatial scale and discretization, (2) the extent to which data are available to physically define the flow system, (3) the type of boundary conditions and controlling parameters used, (4) uncertainty in the distribution of pumping, and (5) uncertainty in field-scale hydraulic properties. The model could be improved with more accurate estimates of ground-water pumpage and better characterization of recharge and discharge.

Florida, Georgia, South Carolina↗
USGS2006

Quantification of the contribution of nitrogen from septic tanks to ground water in Spanish Springs Valley, Nevada

Analysis of total dissolved nitrogen concentrations from soil water samples collected within the soil zone under septic tank leach fields in Spanish Springs Valley, Nevada, shows a median concentration of approximately 44 milligrams per liter (mg/L) from more than 300 measurements taken from four septic tank systems. Using two simple mass balance calculations, the concentration of total dissolved nitrogen potentially reaching the ground-water table ranges from 25 to 29 mg/L. This indicates that approximately 29 to 32 metric tons of nitrogen enters the aquifer every year from natural recharge and from the 2,070 houses that use septic tanks in the densely populated portion of Spanish Springs Valley. Natural recharge contributes only 0.25 metric tons because the total dissolved nitrogen concentration of natural recharge was estimated to be low (0.8 mg/L). Although there are many uncertainties in this estimate, the sensitivity of these uncertainties to the calculated load is relatively small, indicating that these values likely are accurate to within an order of magnitude. The nitrogen load calculation will be used as an input function for a ground-water flow and transport model that will be used to test management options for controlling nitrogen contamination in the basin.

Scientific Investigations Report↗
USGS2004

Age and source of water in springs associated with the Jacksonville Thrust Fault Complex, Calhoun County, Alabama

Water from wells and springs accounts for more than 90 percent of the public water supply in Calhoun County, Alabama. Springs associated with the Jacksonville Thrust Fault Complex are used for public water supply for the cities of Anniston and Jacksonville. The largest ground-water supply is Coldwater Spring, the primary source of water for Anniston, Alabama. The average discharge of Coldwater Spring is about 32 million gallons per day, and the variability of discharge is about 75 percent. Water-quality samples were collected from 6 springs and 15 wells in Calhoun County from November 2001 to January 2003. The pH of the ground water typically was greater than 6.0, and specific conductance was less than 300 microsiemens per centimeter. The water chemistry was dominated by calcium, carbonate, and bicarbonate ions. The hydrogen and oxygen isotopic composition of the water samples indicates the occurrence of a low-temperature, water-rock weathering reaction known as silicate hydrolysis. The residence time of the ground water, or ground-water age, was estimated by using analysis of chlorofluorocarbon, sulfur hexafluoride, and regression modeling. Estimated ground-water ages ranged from less than 10 to approximately 40 years, with a median age of about 18 years. The Spearman rho test was used to identify statistically significant covariance among selected physical properties and constituents in the ground water. The alkalinity, specific conductance, and dissolved solids increased as age increased; these correlations reflect common changes in ground-water quality that occur with increasing residence time and support the accuracy of the age estimates. The concentration of sodium and chloride increased as age increased; the correlation of these constituents is interpreted to indicate natural sources for chloride and sodium. The concentration of silica increased as the concentration of potassium increased; this correlation, in addition to the isotopic data, is evidence that silicate hydrolysis of clay minerals occurred. The geochemical modeling program NETPATH was used to investigate possible mixing scenarios that could yield the chemical composition of water collected from springs associated with the Jacksonville Thrust Fault Complex. The results of NETPATH modeling suggest that the primary source of water in Coldwater Spring is a deep aquifer, and only small amounts of rainwater from nearby sources are discharged from the spring. Starting with Piedmont Sports Spring and moving southwest along a conceptual ground-water flow path that parallels the Jacksonville Thrust Fault Complex, NETPATH simulated the observed water quality of each spring, in succession, by mixing rainwater and water from the spring just to the northeast of the spring being modeled. The percentage of rainwater and ground water needed to simulate the quality of water flowing from the springs ranged from 1 to 25 percent rainwater and 75 to 99 percent ground water.

Scientific Investigations Report↗
USGS2005

Borehole-geophysical and hydraulic investigation of the fractured-rock aquifer near the University of Connecticut Landfill, Storrs, Connecticut, 2000 to 2001

An integrated borehole-geophysical and hydraulic investigation was conducted at the former landfill area near the University of Connecticut in Storrs, Connecticut, where solvents and landfill leachate have contaminated a fractured-bedrock aquifer. Borehole-geophysical techniques and hydraulic methods were used to characterize the site bedrock lithology and structure, fractures, and hydraulic properties. The geophysical and hydraulic methods included conventional logs, borehole imaging, borehole radar, flowmeter under ambient- and stressed hydraulic conditions, and discrete-zone hydraulic testing, sampling, and monitoring. The conventional geophysical-logging methods included caliper, deviation, electromagnetic induction, gamma, specific conductance, and fluid temperature. The advanced methods included optical and acoustic imaging of the borehole wall, heat-pulse flowmeter, and directional radar reflection. Borehole-geophysical methods were used to further define conductive features identified with surface-geophysical methods in the first phase of the investigation. The results of the surface- and borehole-geophysical logging were evaluated in an iterative and integrated manner to develop a conceptual model of ground-water flow at the site. The rock type, foliation, and fractures at the site were characterized from high-resolution optical televiewer (OTV) images of rocks penetrated by the boreholes and were compared to drilling logs and conventional geophysical logs. The rocks are interpreted as fine- to mediumgrained quartz-feldspar-biotite-garnet gneiss and schist with local intrusions of quartz diorite and pegmatite and minor concentrations of sulfide mineralization similar to rocks described as the Bigelow Brook Formation on regional geologic maps. Layers containing high concentrations of sulfide minerals appear as high electrical conductivity zones on electromagnetic-induction and borehole-radar logs. Foliation in the rocks generally strikes to the southwest and northeast, and dips to the northwest and southeast consistent with previous investigations in this area. The orientation of foliation, however, varies locally and with depth in some of the boreholes. These results are consistent with geologic mapping that has identified small-scale folding. The orientations of the transmissive fractures identified in the six boreholes logged for this investigation are similar to the fracture orientations mapped in a previous investigation. Many of these fractures are oriented with a north-northwest strike and have a shallow dip to the west. Other transmissive fractures have a southwest strike and dip at shallow angles to the northwest, and some strike roughly east-west and dip to the north and south. Flowmeter logging was used to identify transmissive fractures and to estimate the hydraulic properties in the boreholes. Ambient down flow was measured in one borehole, and ambient up flow and down flow were measured in another borehole. The other four bedrock boreholes did not have measurable vertical flow. Under low-rate pumping conditions (0.25 to 0.5 gallons per minute), one to three inflow zones were identified in each well. Commonly, fractures that are active under ambient conditions contribute to the well under pumping conditions. The ambient conditions were incorporated into the determination of the relative proportions of transmissivity. Specific capacity and transmissivity were determined for these open-hole low-rate pumping tests. Quasi-steady-state water levels were reached in four of the boreholes, including MW201R, MW204R, MW302R, and W202-NE. When pumped at low-rate conditions for 0.5 to 4 hours, the specific capacity ranged from 0.03 to 0.18 gallons per minute per foot. The open-hole transmissivity estimates ranged from 4.9 to 30 feet squared per day (ft2/d). Open-hole transmissivity was determined for boreholes that did not reach quasi-steady-state conditions under low-rate pumping conditions. Transmissivity was estimated for MW201R, MW202R, and MW203R using non-equilibrium methods, pumping rate, and the transient drawdown data to estimate the open-hole transmissivity. Transmissivity in these boreholes ranged from 0.98 to 3.2 ft 2 /d. The transmissivity and head of individual fractures or zones of fractures were estimated from heat-pulse flowmeter data acquired under ambient and stressed conditions. In the absence of ambient flow, data from two profiles of heat-pulse flowmeter data under two different stressed conditions were used to estimate the transmissivity and head of individual fracture zones. Only two boreholes, MW302R and W202-NE, had sufficient data for these analyses. The estimated transmissivity of individual transmissive zones ranged from 1.2 to 9.2 ft 2 /d. The transmissivity values determined by this numerical simulation method were less than the open-hole estimations, which were 15 and 30 ft 2 /d. Transmissivity also was measured directly over discrete intervals of the borehole using a straddle-packer apparatus and constant-rate pumping tests. Pumping rates were less than or equal to 0.25 gallons per minute. These discretezone single-hole pumping tests were conducted over a short period of time, usually about 30 minutes to 1 hour in duration. Pumping continued until the test zone reached a steady-state water level or until it was determined that the zone could not yield water at the pumped rate. The estimated transmissivity of individual transmissive zones ranged from about 0.21 to 11 ft 2 /d. The zone at a depth of 197 feet in W202-NE was the only zone that had discrete-interval testing with a straddle packer and sufficient heat-pulse flowmeter data for modeling the flow and estimating transmissivity and head. The two methods produced similar results. The straddle-packer method estimated a transmissivity of 4.7 ft 2 /d, and the heat-pulse flowmeter modeling results estimated a transmissivity of 6.9 ft 2 /d. A comparison of the transmissivity estimates indicate estimates typically are within an order of magnitude. The heat-pulse flowmeter methods used in this investigation to determine transmissivity of the boreholes and the individual fractures measure only the upper two or three orders of magnitude of transmissivity. Hence, other fractures in these boreholes permit the movement of water; their transmissivities, however, are lower than the detection limits of the methods that were used for this investigation and very small compared to the transmissive fractures that were studied. The data collected in this investigation were used to design discrete-zone monitoring systems for four of the boreholes used for monitoring. The results of the investigation are useful for refining the conceptual site model of ground-water flow, and for providing critical information for interpreting the results of water-quality sampling.

Connecticut↗