Search USGS⌕ Search

SEARCH · Search USGS

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 631 records · Page 35Linked to original sources
USGS1988

A modular three-dimensional finite-difference ground-water flow model

This report presents a finite-difference model and its associated modular computer program. The model simulates flow in three dimensions. The report includes detailed explanations of physical and mathematical concepts on which the model is based and an explanation of how those concepts are incorporated in the modular structure of the computer program. The modular structure consists of a Main Program and a series of highly independent subroutines called 'modules.' The modules are grouped into 'packages.' Each package deals with a specific feature of the hydrologic system which is to be simulated, such as flow from rivers or flow into drains, or with a specific method of solving linear equations which describe the flow system, such as the Strongly Implicit Procedure or Slice-Successive Overrelaxation. The division of the program into modules permits the user to examine specific hydrologic features of the model independently. This also facilita development of additional capabilities because new packages can be added to the program without modifying the existing packages. The input and output systems of the computer program are also designed to permit maximum flexibility. Ground-water flow within the aquifer is simulated using a block-centered finite-difference approach. Layers can be simulated as confined, unconfined, or a combination of confined and unconfined. Flow associated with external stresses, such as wells, areal recharge, evapotranspiration, drains, and streams, can also be simulated. The finite-difference equations can be solved using either the Strongly Implicit Procedure or Slice-Successive Overrelaxation. The program is written in FORTRAN 77 and will run without modification on most computers that have a FORTRAN 77 compiler. For each program ,module, this report includes a narrative description, a flow chart, a list of variables, and a module listing.

Techniques of Water-Resources Investigations↗
USGS1996

Hydrologic feasibility of water-supply-development alternatives in Cape May County, New Jersey

Increasing public-supply withdrawals in Cape May County, New Jersey associated with increasing residential and seasonal tourist populations have led to regionally lowered ground-water levels, a reversal of ground-water flow directions toward onshore, and landward encroachment of saltwater in the shallow aquifer system. The three aquifers composing the shallow system are, in order of increasing depth, the unconfined Holly Beach water-bearing zone and the confined estuarine sand and Cohansey aquifers. The changes to the ground-water system have been greatest in the confined aquifers near the three major well fields on the Cape May peninsula. Formerly productive water-supply wells have been abandoned because of saltwater contamination. Concern about anthropogenic contamination has prevented shifting of withdrawals to the unconfined aquifer. Surface- water sources have also been little used. Further development on the peninsula involving increased water demand will exacerbate the current saltwater-encroachment problems. The purpose of this study was to test the feasibility of possible water-supply-development alternatives by use of predictive ground-water flow simulations. The alternatives involve (1) injection of tertiary- treated wastewater to replenish aquifer storage and create a hydraulic barrier to saltwater encroachment, (2) withdrawal of brackish water in order to create a hydraulic barrier, (3) conjunctive use of ground water and surface water, enabling the reduction of ground-water withdrawals, and (4) redistribution of withdrawals inland to the unconfined aquifer. Results of these simulations can potentially be used in the design of a water-supply-development strategy that preserves supply and a monitoring program that ensures early warning of saltwater encroachment, thereby allowing sufficient time for development of an alternative supply. The water-supply- development alternatives were evaluated by comparison of results of predictive simulations made with a previously calibrated ground-water flow model of the shallow aquifer system. The quasi-three-dimensional sharp-interface model was calibrated to 1988 annual average hydrologic conditions. The planning period for the predictive simulations is 1989-2049. For the planning period, total public-supply withdrawals were increased 100 percent over average 1983-88 withdrawals. Results of a baseline simulation involving only the increased withdrawals were compared to each of the simulated alternatives, which also include the withdrawals. Hydraulic heads, saltwater- freshwater interface movement, and ground-water flows were compared. Simulation results indicate that the barrier-injection or barrier-withdrawal scheme could be useful in managing the water supply for a specific location. The conjunctive- use scheme would provide a marginal regional hydrologic benefit. Redistribution of withdrawals appears to be the only regional alternative that would result in recovery of ground-water levels and would substantially slow saltwater encroachment; however, anthropogenic land-surface contamination of the unconfined aquifer would have to be considered if the redistribution alternative is acted upon.

Water-Resources Investigations Report↗
USGS1995

Water budgets, water quality, and analysis of nutrient loading of the Winter Park Chain of Lakes, central Florida, 1989-92, with a section on littoral vegetation

The Winter Park chain of lakes (Lakes Maitland, Virginia, Osceola, and Mizell) has a combined area of about 900 acres, an immediate drainage area of about 3,100 acres, and mean depths ranging from 11 to 15 feet. The lakes are an important recreational resource for the surrounding communities, but there is concern about the possible effects of stormwater runoff and seepage of nutrient-enriched ground water on the quality of water in the lakes. The lakes receive water from several sources: rainfall on lake surfaces, inflow from other surface-water bodies, stormflow that enters the lakes through storm drains or by direct runoff from land adjacent to the lakes, and ground-water seepage. Water leaves the lakes by evaporation, surface outflow, and ground-water outflow. Of the three, only surface outflow can be measured directly. Rainfall, surface inflow and outflow, and lake-stage data were collected from October 1, 1989, to September 30, 1992. Stormflow, evaporation and ground-water inflow and outflow were estimated for the 3 years of the study. Ground-water outflow was calculated by evaluating the rate of lake-stage decline during dry periods. Estimated ground-water outflow was compared to downward leakage rates estimated by ground-water flow models. Lateral ground-water inflow from surficial sediments was calculated as the residual of the flow budget. Flow budgets were calculated for the 3 years of the study. In water year 1992 (a year with about average rainfall), inflow consisted of rainfall, 48 inches; stormflow, 15 inches; surface inflow, 67 inches; and ground water, 40 inches. The calculated outflows were evaporation, 47 inches; surface outflow, 90 inches; and ground water, 33 inches. Water-quality data also were used to calculate nutrient budgets for the lakes. Bimonthly water samples were collected from the lakes and at surface inflow and outflow sites, and were analyzed for physical characteristics, dissolved oxygen, pH, specific conductance, major ions, the nutrients nitrogen and phosphorus, and chlorophyll (collected at lake sites only). Specific conductance ranged from about 190 to 230 microsiemens per centimeter at 25 degrees Celsius in Lakes Maitland, Virginia and Osceola and from about 226 to 260 microsiemens per centimeter at 25 degrees Celsius in Lake Mizell. The median concentrations of total ammonia-plus-organic nitrogen in all the lakes ranged from 0.79 to 0.99 milligrams per liter. Median total phosphorus concentrations ranged from less than 0.02 to 0.20 milligrams per liter. Stormwater samples were collected for 17 storms at one storm-drain site and 16 storms at another storm-drain site on Lake Osceola. Median total nitrogen concentrations at the sites were 2.23 and 3.06 milligrams per liter and median total phosphorus concentrations were 0.34 and 0.40 milligrams per liter. The water quality in the Winter Park lakes generally is fair to good, based on a trophic-state index used by the Florida Department of Environmental Protection for assessing the tropic state of Florida lakes. This index was determined from median total nitrogen, total phosphorus, and chlorophyll-a concentrations, and median Secchi-disk transparency for all lakes for the period September 1989 to June 1992. Based on a one-time sampling of 20 sites around the lakes, surficial ground-water quality is highly variable. Nutrient concentrations were highly variable and could not be correlated to the proximity of septic tanks. Fertilizer probably is the primary source of nutrients in the surficial ground water. Nutrient budgets were calculated for the lakes for the 3 years of the study. The most variable source of nutrient loading to the lakes is stormwater. Nutrient-loading modeling indicates that reduction of nutrients in stormflow probably would improve lake-water quality. However, even with complete removal of nitrogen and phosphorus from stormwater, the lakes might still be mesotrophic with respect to both nutrients during periods of below average rainfall because of the input from the other sources of inflow to the lakes. Littoral vegetation in the lakes was surveyed in March 1992. The length of shoreline containing vegetation was 44 percent in Lake Maitland, 62 percent in Lake Virginia, 46 percent in Lake Osceola, and 76 percent in Lake Mizell. The types of vegetation present generally were similar for all four lakes.

Florida↗
USGS1998

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Naval Station Mayport, Florida

Ground-water flow through the surficial aquifer system at Naval Station Mayport near Jacksonville, Florida, was simulated with a two-layer finite-difference model as part of an investigation conducted by the U.S. Geological Survey. The model was calibrated to 229 water-level measurements from 181 wells during three synoptic surveys (July 17, 1995; July 31, 1996; and October 24, 1996). A quantifiable understanding of ground-water flow through the surficial aquifer was needed to evaluate remedial-action alternatives under consideration by the Naval Station Mayport to control the possible movement of contaminants from sites on the station. Multi-well aquifer tests, single-well tests, and slug tests were conducted to estimate the hydraulic properties of the surficial aquifer system, which was divided into three geohydrologic units�an S-zone and an I-zone separated by a marsh-muck confining unit. The recharge rate was estimated to range from 4 to 15 inches per year (95 percent confidence limits), based on a chloride-ratio method. Most of the simulations following model calibration were based on a recharge rate of 8 inches per year to unirrigated pervious areas. The advective displacement of saline pore water during the last 200 years was simulated using a particle-tracking routine, MODPATH, applied to calibrated steady-state and transient models of the Mayport peninsula. The surficial aquifer system at Naval Station Mayport has been modified greatly by natural and anthropogenic forces so that the freshwater flow system is expanding and saltwater is being flushed from the system. A new MODFLOW package (VAR1) was written to simulate the temporal variation of hydraulic properties caused by construction activities at Naval Station Mayport. The transiently simulated saltwater distribution after 200 years of displacement described the chloride distribution in the I-zone (determined from measurements made during 1993 and 1996) better than the steady-state simulation. The advective movement of contaminants from selected sites within the solid waste management units to discharge points was simulated using MODPATH. Most of the particles were discharged to the nearest surface-water feature after traveling less than 1,000 feet in the ground-water system. Most areas within 1,000 feet of a surface-water feature or storm sewer had traveltimes of less than 50 years, based on an effective porosity of 40 percent. Contributing areas, traveltimes, and pathlines were identified for 224 wells at Naval Station Mayport under steady-state and transient conditions by back-tracking a particle from the midpoint of the wetted screen of each well. Traveltimes to contributing areas that ranged between 15 and 50 years, estimated by the steady-state model, differed most from the transient traveltime estimates. Estimates of traveltimes and pathlines based on steady-state model results typically were 10 to 20 years more and about twice as long as corresponding estimates from the transient model. The models differed because the steady-state model simulated 1996 conditions when Naval Station Mayport had more impervious surfaces than at any earlier time. The expansion of the impervious surfaces increased the average distance between contributing areas and observation wells.

Florida↗
USGS1996

Geohydrology of stratified drift and streamflow in the Deerfield River basin, northwestern Massachusetts

This report presents the results of a study of the geohydrology of stratified drift and streamflow in the Deerfield River Basin, northwestern Massachusetts. Detailed hydrologic information is needed to plan for the optimal use of ground-water and surface-water resources and for development of new drinking-water supplies in the basin. Sources and percentage of water available for recharge on an annual basis from October 1993 to September 1994, to the fine-grained stratified-drift in a narrow valley bordered by upland till and bedrock were: (1) direct infiltration of precipitation on the valley (30 percent); (2) tributary loss from an upland brook as it crosses the valley (7 percent); and (3) ground- and surface-water runoff from the uplands (63 percent). Seventy percent of recharge was available from upland sources. Seasonal variation in recharge caused changes in ground-water levels and flow directions. In early spring, the direction of flow is toward the valley axis, but in late summer, the direction of flow is nearly parallel to the valley axis. Field observations and results of a ground-water flow simulation indicated that water available for recharge was greater than actual recharge during the spring snowmelt and during intense precipitation events. In 1994, estimates of water available for recharge were greater than actual recharge by 10 percent in March and by 60 percent in April; actual recharge to the valley on an annual basis from October 1993 to September 1994 was 20 percent less than original estimates. A map showing thickness of stratified drift in the Connecticut Valley Lowlands indicates a deep north-south trending buried valley. Maximum thickness of the stratified drift is 385 feet. Interpretation of a seismic-reflection survey indicates fine-grained stratified drift may be underlain by coarse-grained deposits ranging in thickness from 0 to 150 feet. Hydraulic properties of the stratified drift were calculated from ground-water-level fluctuations induced by river stage changes using a ground-water-flow model for a site adjacent to the Deerfield River. A comparison of measured and simulated heads resulted in a vertical riverbed hydraulic conductivity of 3 feet per day, anisotropic ratio of horizontal to vertical hydraulic conductivity of 40:1, and storage of 0.040 and 0.0002 for the unconfined and confined layers of the stratified drift. Hydraulic diffusivity (transmissivity divided by the unconfined storage) at the site is about 168,000 feet squared per day. Streamflows at times of low flow were determined for 27 sites that drain areas ranging from 0.57 to 15.8 percent stratified drift. Streamflows exceeded between 80 and 99 percent of the time were determined for sites on the unregulated tributaries to the Deerfield River. Streamflows per square mile of drainage area were greatest from sites at the downstream ends of the North River-Colrain and the Deerfield River-Charlemont stratified-drift valleys. Flow-duration curves for three continuous streamflow-gaging stations on the regulated Deerfield River were compared to flow-duration curves for three continuous streamflow-gaging stations on unregulated tributaries to show the effects of dam regulation on streamflow. Flow- duration curves constructed using instantaneous discharges for the three regulated gaging stations have flat sections that correspond to the predominant streamflows when water is being released from storage from the dams.

Water-Resources Investigations Report↗
USGS2012

Helicopter electromagnetic survey of the Model Land Area, Southeastern Miami-Dade County, Florida

This report describes a helicopter electromagnetic survey flown over the Model Land Area in southeastern Miami-Dade County, Florida, to map saltwater intrusion in the Biscayne aquifer. The survey, which is located south and east of Florida City, Florida, covers an area of 115 square kilometers with a flight-line spacing of 400 meters. A five-frequency, horizontal, coplanar bird with frequencies ranging from 400 to 100,000 Hertz was used. The data were interpreted using differential resistivity analysis and inversion to produce cross sections and resistivity depth-slice maps. The depth of investigation is as deep as 100 meters in freshwater-saturated portions of the Biscayne aquifer and the depth diminishes to about 50 meters in areas that are intruded by saltwater. The results compare favorably with ground-based, time-domain electromagnetic soundings and induction logs from observation wells in the area. The base of a high-resistivity, freshwater-saturated zone mapped in the northern 2 kilometers of the survey area corresponds quite well with the base of the surficial aquifer that has been determined by drilling. In general, saltwater in the survey area extends 9 to 12 kilometers inland from the coast; however, there is a long nose of saltwater centered along the Card Sound Road Canal that extends 15 kilometers inland. The cause of this preferential intrusion is likely due to uncontrolled surface flow along the canal and subsequent leakage of saltwater into the aquifer. Saltwater also extends farther inland in the area between U.S. Highway 1 and Card Sound Road than it does to the west of this area. Until 1944, a railroad grade occupied the current location of U.S. Highway 1. Borrow ditches associated with the railroad grade connected to Barnes Sound and allowed saltwater to flow during droughts and storm surges to within a few kilometers of Florida City. Relicts of this saltwater that settled to the bottom of the Biscayne aquifer can be seen in the helicopter electromagnetic data. The area to the west of U.S. Highway 1 is more resistive in the upper 10 meters than the area to the east of the road; this reflects the influence of surface-water flows that are blocked by U.S. Highway 1. Between Card Sound Road and U.S. Highway 1, resistivities are slightly lower compared to adjacent areas. In the southern portion of the survey area, the surficial aquifer underlying the Biscayne aquifer is more resistive; this indicates that it contains fresher water than that found at the base of the Biscayne aquifer.

Florida↗
USGS1991

Hydrogeology and simulation of ground-water flow in the Rochester area, southeastern Minnesota, 1987-88

Ground-water flow in the St. Peter-Prairie du Chien-Jordan aquifer was studied in a 700 square-mile area surrounding Rochester, Minnesota. The aquifer consisting of sandstone, limestone, and dolomite is locally confined by the Decorah-Platteville-Glenwood sequence of shales and limestones. Regional flow in the aquifer is from a ground-water divide on the western, southern, and eastern sides of the city toward various rivers. A 140-square-mile area of the aquifer is a source of water supply for the Rochester area. A cone of depression in the potentiometric surface of the aquifer throughout most of the year is centered around high-capacity (greater than about 200 gallons per minute) wells in downtown Rochester. The cone covered an area of about 2.3 square miles in August 1988. Most streams in the area gain water from the ground-water system. One reach of the South Fork Zumbro River, however, loses water to the system. This loss is probably caused by the pumping of nearby high-capacity wells. A ground-water-flow model was used to simulate the effects of an extended drought near Rochester. Conclusions based on the simulations are that (1) reduced recharge and increased pumping, conditions that could exist during a 3-year drought, would probably lower water levels 5 to 10 feet regionally and more than 30 feet in the city; (2) pumping of six additional municipal wells on the perimeter of the city would lower regional water levels about 1 to 5 feet; and (3) that water levels would recover 1 to 18 feet if pumping from six municipal wells in downtown Rochester were discontinued. The area encompasses five recharge zones that can be delineated on the basis of recharge rate. About 54 percent of recharge to the aquifer in the area contributing water to Rochester is from a zone along the edge of the Decorah-Platteville-Glenwood confining unit. About 10 percent of recharge in this contributing area is to the sewered area of Rochester.

Minnesota↗
USGS2007

Geostatistical Modeling of Sediment Abundance in a Heterogeneous Basalt Aquifer at the Idaho National Laboratory, Idaho

The spatial distribution of sediment in the eastern Snake River Plain aquifer was evaluated and modeled to improve the parameterization of hydraulic conductivity (K) for a subregional-scale ground-water flow model being developed by the U.S. Geological Survey. The aquifer is hosted within a layered series of permeable basalts within which intercalated beds of fine-grained sediment constitute local confining units. These sediments have K values as much as six orders of magnitude lower than the most permeable basalt, and previous flow-model calibrations have shown that hydraulic conductivity is sensitive to the proportion of intercalated sediment. Stratigraphic data in the form of sediment thicknesses from 333 boreholes in and around the Idaho National Laboratory were evaluated as grouped subsets of lithologic units (composite units) corresponding to their relative time-stratigraphic position. The results indicate that median sediment abundances of the stratigraphic units below the water table are statistically invariant (stationary) in a spatial sense and provide evidence of stationarity across geologic time, as well. Based on these results, the borehole data were kriged as two-dimensional spatial data sets representing the sediment content of the layers that discretize the ground-water flow model in the uppermost 300 feet of the aquifer. Multiple indicator kriging (mIK) was used to model the geographic distribution of median sediment abundance within each layer by defining the local cumulative frequency distribution (CFD) of sediment via indicator variograms defined at multiple thresholds. The mIK approach is superior to ordinary kriging because it provides a statistically best estimate of sediment abundance (the local median) drawn from the distribution of local borehole data, independent of any assumption of normality. A methodology is proposed for delineating and constraining the assignment of hydraulic conductivity zones for parameter estimation, based on the locally estimated CFDs and relative kriging uncertainty. A kriging-based methodology improves the spatial resolution of hydraulic property zones that can be considered during parameter estimation and should improve calibration performance and sensitivity by more accurately reflecting the nuances of sediment distribution within the aquifer.

Scientific Investigations Report↗
USGS2008

Digital surfaces and thicknesses of selected hydrogeologic units within the Mississippi Embayment Regional Aquifer Study (MERAS)

Digital surfaces of selected Tertiary and younger age hydrogeologic units within the Mississippi embayment aquifer system were created using more than 2,600 geophysical logs for an area that covers approximately 70,000 square miles and encompasses parts of eight states. The digital surfaces were developed to define and display the hydrogeologic framework for the Mississippi Embayment Regional Aquifer Study (MERAS). The digital surfaces also provide a foundation of the selected hydrogeologic units for development of a steady-state and transient regional ground-water flow model of the Mississippi embayment aquifer system from the top of the Midway confining unit upwards to land surface. The ground-water flow model is under development as part of the U.S. Geological Survey Ground-Water Resources Program. Using a Geographic Information System, nine digital surfaces of the tops of selected hydrogeologic units were created using the Australian National University Digital Elevation Model method as an interpolation scheme. Thickness maps also were constructed using the Geographic Information System by calculating the difference between the altitude of the interpreted base of an overlying unit and the altitude of the interpreted top of an underlying unit. In general, the highest hydrogeologic unit altitudes are located along the eastern edge of the study area in the outcrop, and the lowest altitudes, in general, are located along the southern edge of the study area along the axis of the embayment. The Mississippi River Valley alluvial aquifer and the lower Claiborne aquifer are the thinnest aquifers of importance in the study area; the thickest aquifer of importance is the middle Claiborne aquifer.

Mississippi Embayment Regional Aquifer↗
USGS1988

Simulation of ground-water flow in aquifers along the Susquehanna River in Columbia County, Pennsylvania

A numerical model of groundwater flow was developed for a 10.3 sq mi area along the Susquehanna River in Columbia County, east central Pennsylvania. Groundwater in the model area primarily is in secondary openings in the carbonate--and clastic-rock aquifers and primary openings in the glacial-outwash aquifer that discontinuously overlies bedrock. The groundwater flow model was calibrated under average steady-state conditions for 1981. The simulated 1981 water budget indicates an average inflow rate of 7.24 cu ft/sec. Of this, 93% is recharge from precipitation and 6.6% is boundary flow. 62% of the outflow is leakage to streams, 21% to pumpage, and 17% to evapotranspiration. The model was calibrated under transient conditions for December 22, 1980 through April 21, 1982. Water level fluctuations caused by natural stresses were more successfully simulated than those caused by pumping stresses. Three 10-year, hypothetical stress periods were simulated with the calibrated, transient model. The general impact of three pumping schemes under hypothetical drought and drought-recovery conditions were simulated. (USGS)

Pennyslvania↗
USGS1985

Potential incremental seepage losses in an alluvial channel in the Rio Grande Basin, New Mexico

A two-dimensional, digital, cross-sectional model was used to simulate seepage of water from an alluvial channel, which had the general characteristic of the Rio Grande channel, into the underlying alluvium within the reach from Cochiti Dam to Elephant Butte Reservoir. Seepage rates were determined for losing and gaining reaches, and reaches affected by pumping of ground water. The seepage rates were computed for stream surcharges (height of additional water applied on top of base flow) ranging from 0.5 foot to 3 feet and for application periods ranging from 1 to 100 days. The net seepage rates, which were nearly identical for each type of reach, ranged from 0.0 cubic foot per second per mile of channel length for a 0.5 foot surcharge applied for 1 day to 0.37 cubic foot per second per mile of channel length for a 3 feet surcharge applied for 100 days, followed by a 180 day seepage return flow from the aquifer. (USGS)

Water-Resources Investigations Report↗
USGS2007

A Science Plan for a Comprehensive Regional Assessment of the Atlantic Coastal Plain Aquifer System in Maryland

The Maryland Coastal Plain region is, at present, largely dependent upon ground water for its water supply. Decades of increasing pumpage have caused ground-water levels in parts of the Maryland Coastal Plain to decline by as much as 2 feet per year in some areas of southern Maryland. Continued declines at this rate could affect the long-term sustainability of ground-water resources in Maryland's heavily populated Coastal Plain communities and the agricultural industry of the Eastern Shore. In response to a recommendation in 2004 by the Advisory Committee on the Management and Protection of the State's Water Resources, the Maryland Geological Survey and the U.S. Geological Survey have developed a science plan for a comprehensive assessment that will provide new scientific information and new data management and analysis tools for the State to use in allocating ground water in the Coastal Plain. The comprehensive assessment has five goals aimed at improving the current information and tools used to understand the resource potential of the aquifer system: (1) document the geologic and hydrologic characteristics of the aquifer system in the Maryland Coastal Plain and appropriate areas of adjacent states; (2) conduct detailed studies of the regional ground-water-flow system and water budget for the aquifer system; (3) improve documentation of patterns of water quality in all Coastal Plain aquifers, including the distribution of saltwater; (4) enhance ground-water-level, streamflow, and water-quality-monitoring networks in the Maryland Coastal Plain; and (5) develop science-based tools to facilitate sound management of the ground-water resources in the Maryland Coastal Plain. The assessment, as designed, will be conducted in three phases and if fully implemented, is expected to take 7 to 8 years to complete. Phase I, which was initiated in January 2006, is an effort to assemble all the information and investigation tools needed to do a more comprehensive assessment of the aquifer system. The work will include updating the hydrogeologic framework, developing a Geographic Information System-based aquifer information system, refinement of water-use information, assessment of existing water-quality data, and development of detailed plans for ground-water-flow and management models. Phase II is an intensive study phase during which a regional ground-water-flow model will be developed and calibrated for the entire region of Maryland in the Atlantic Coastal Plain as well as appropriate areas of Delaware and Virginia. The model will be used to simulate flow and water levels in the aquifer system and to study the water budget of the system. The model analysis will be based on published information but will be supplemented with field investigations of recharge and leakage in the aquifer system. Localized and finely discretized ground-water-flow models that are embedded in the regional model will be developed for selected areas of heavy withdrawals. Other modeling studies will be conducted to better understand flow in the unconfined parts of the aquifer system and to support the recharge studies. Phase II will also include selected water-quality studies and a study to determine how hydrologic and water-quality-monitoring networks need to be enhanced to appropriately assess the sustainability of the Coastal Plain aquifer system. Phase III will be largely devoted to the development and application of a ground-water optimization model. This model will be linked to the ground-water-flow model to create a model package that can be used to test different water-management scenarios. The management criteria that will be used to develop these scenarios will be determined in consultation with a variety of state and local stakeholders and policy makers in Phases I and II of the assessment. The development of the aquifer information system is a key component of the assessment. The system will store all relevant aquifer data

Open-File Report↗
USGS2008

Recovery of Ground-Water Levels From 1988 to 2003 and Analysis of Potential Water-Supply Management Options in Critical Area 1, East-Central New Jersey

Water levels in four confined aquifers in the New Jersey Coastal Plain within Water Supply Critical Area 1 have recovered as a result of reductions in ground-water withdrawals initiated by the State in the late 1980s. The aquifers are the Wenonah-Mount Laurel, the Upper and Middle Potomac-Raritan-Magothy, and Englishtown aquifer system. Because of increased water demand due to increased development in Monmouth, Ocean, and Middlesex Counties, five base and nine alternate management models were designed for the four aquifers to evaluate the effects resulting from potential reallocation of part of the Critical Area 1 reductions in withdrawals. The change in withdrawals and associated water-level changes in the aquifers for 1988-2003 are discussed. Generally, withdrawals decreased 25 to 30 Mgal/d (million gallons per day), and water levels increased 0 to 80 ft (feet). The Regional Aquifer-System Analysis (RASA) ground-water-flow model of the New Jersey Coastal Plain developed by the U.S. Geological Survey was used to simulate ground-water flow and optimize withdrawals using the Ground-Water Management Process (GWM) for MODFLOW. Results of the model were used to evaluate the effects of several possible water-supply management options in order to provide the information to water managers. The optimization method, which provides a means to set constraints that support mandated hydrologic conditions, then determine the maximum withdrawals that meet the constraints, is a more cost-effective approach than simulating a range of withdrawals to determine the effects on the aquifer system. The optimization method is particularly beneficial for a regional-scale study of this kind because of the large number of wells to be evaluated. Before the model was run, a buffer analysis was done to define an area with no additional withdrawals that minimizes changes in simulated streamflow in aquifer outcrop areas and simulated movement of ground water toward the wells from areas of possible high chloride concentrations in the northern and southern parts of the Critical Area. Five base water-supply management models were developed. Each management model has an objective function, decision variables, and constraints. Two of the five management models were test cases: clean slate option and reallocation from the Wenonah-Mount Laurel aquifer and Englishtown aquifer system to small volume wells for potable water use. Nine other models also were developed as part of a trade-off analysis between withdrawal amounts and constraint values. The 14 management models included current (2003) or regularly spaced well locations with variations on the constraints of ground-water head, drawdown, velocity at the 250-mg/L (milligram per liter) isochlor, and withdrawal rate. Results of each management model were evaluated in terms of withdrawals, heads, saltwater intrusion, and source of water by aquifer. Each trade-off curve was defined by using six to nine separate management model runs. Results of the management models designed in this study indicate that a withdrawal reallocation of 5 to 20 Mgal/d within Critical Area 1 would increase the area of heads below -30 ft and the velocity at the 250-mg/L isochlor by up to 4 times that of the simulated 2003 results; the range of values are 0 to 521 square miles and 1 to 20 feet per year, respectively. The increase in area of heads below -30 ft was larger in the Middle Potomac-Raritan-Magothy aquifer than in other aquifers because that area was negligible in 2003. The range of modeled withdrawals is closely tied to management-model design. Interpretation of management model results is provided as well as a discussion of limitations.

Scientific Investigations Report↗
USGS1987

Simulation of the ground-water flow system and proposed withdrawals in the northern part of Vekol Valley, Arizona

Pursuant to the Ak-Chin Indian Community Water Rights Settlement Act (Public Law 95-328-enacted on July 28, 1978) a study was undertaken to assess the effect of proposed groundwater withdrawal from Federal lands near the reservation. The first area to be evaluated was the northern part of the Vekol Valley. The evaluation was made using a numerical model based on detailed geohydrologic concepts developed during the study. The numerical model, which was calibrated to steady-state and transient groundwater conditions in the northern part of Vekol Valley, adequately duplicated the conceptual model and was used to estimate the effect of withdrawing approximately 174,000 acre-ft from the system during a 25-yr period. At the end of the 25-yr period, the water level was drawn down an average of about 95 ft, and about 150,5000 acre-ft of water was removed from storage. The 150,500 acre-ft of water represents 43% of the estimated recoverable groundwater in storage. (Author 's abstract)

Water-Resources Investigations Report↗
USGS2003

Methodology for estimating times of remediation associated with monitored natural attenuation

Natural attenuation processes combine to disperse, immobilize, and biologically transform anthropogenic contaminants, such as petroleum hydrocarbons and chlorinated ethenes, in ground-water systems. The time required for these processes to lower contaminant concentrations to levels protective of human health and the environment, however, varies widely between different hydrologic systems, different chemical contaminants, and varying amounts of contaminants. This report outlines a method for estimating timeframes required for natural attenuation processes, such as dispersion, sorption, and biodegradation, to lower contaminant concentrations and mass to predetermined regulatory goals in groundwater systems. The time-of-remediation (TOR) problem described in this report is formulated as three interactive components: (1) estimating the length of a contaminant plume once it has achieved a steady-state configuration from a source area of constant contaminant concentration, (2) estimating the time required for a plume to shrink to a smaller, regulatoryacceptable configuration when source-area contaminant concentrations are lowered by engineered methods, and (3) estimating the time needed for nonaqueous phase liquid (NAPL) contaminants to dissolve, disperse, and biodegrade below predetermined levels in contaminant source areas. This conceptualization was used to develop Natural Attenuation Software (NAS), an interactive computer aquifers. NAS was designed as a screening tool and requires the input of detailed site information about hydrogeology, redox conditions, and the distribution of contaminants. Because NAS is based on numerous simplifications of hydrologic, microbial, and geochemical processes, the program may introduce unacceptable errors for highly heterogeneous hydrologic systems. In such cases, application of the TOR framework outlined in this report may require more detailed, site-specific digital modeling. The NAS software may be downloaded from the Web site http://www.cee.vt.edu/NAS/ Application of NAS illustrates several general characteristics shared by all TOR problems. First, the distance of stabilization of a contaminant plume is strongly dependent on the natural attenuation capacity of particular ground-water systems. The time that it takes a plume to reach a steady-state configuration, however, is independent of natural attenuation capacity. Rather, the time of stabilization is most strongly affected by the sorptive capacity of the aquifer, which is dependent on the organic matter content of the aquifer sediments, as well as the sorptive properties of individual contaminants. As a general rule, a high sorptive capacity retards a plume.s growth or shrinkage, and increases the time of stabilization. Finally, the time of NAPL dissolution depends largely on NAPL mass, composition, geometry, and hydrologic factors, such as ground-water flow rates. An example TOR analysis for petroleum hydrocarbon NAPL was performed for the Laurel Bay site in South Carolina. About 500 to 1,000 pounds of gasoline leaked into the aquifer at this site in 1991, and the NAS simulations suggested that TOR would be on the order of 10 years for soluble and poorly sorbed compounds, such as benzene and methyl tertiary-butyl ether (MTBE). Conversely, TOR would be on the order of 40 years for less soluble, more strongly sorbed compounds, such as toluene, ethylbenzene, and xylenes (TEX). These TOR estimates are roughly consistent with contaminant concentrations observed over 10 years of monitoring at this site where benzene and MTBE concentrations were observed to decrease rapidly and are approaching regulatory maximum concentration limits, whereas toluene, ethylbenzene, and xylene concentrations decreased at a slower rate and have remained relatively high. An example TOR analysis for petroleum hydrocarbon NAPL was performed for the Laurel Bay site in South Carolina.

South Carolina↗
USGS1980

Effect of pumpage on ground-water levels as modeled in Laramie County, Wyoming

Groundwater is being extensively developed for domestic, agricultural, and industrial use in a 2,320-square mile area in Laramie County, WY., bounded approximately by Horse Creek on the north, Nebraska on the east, Colorado on the south, and pre-Tertiary outcrops on the west. Currently (1977) about 47,300 acres of land are irrigated with groundwater. Groundwater levels are declining in some areas as much as 4 feet per year. The investigation was made to provide State water administrators with data on water-level changes resulting from present (1977) groundwater withdrawals and to provide a means of predicting the future effect of groundwater development. A digital model was developed of the hydrologic system in the post-Cretaceous rocks. The ability of the model to simulate the hydrologic system was determined by comparing the water-level changes measured at 37 observation wells located in areas of irrigation pumping with the water-level changes calculated by the model for 1971-77. Comparison of the measured and calculated changes showed agreement with a root-mean-square deviation of + or - 3.6 feet with 8 feet as the maximum deviation. It is concluded that the model adequately simulates present hydrologic conditions in the post-Cretaceous rocks and may be used to predict the effect of applied stress to the system. (USGS)

Open-File Report↗
USGS2006

Simulation and particle-tracking analysis of ground-water flow near the Savannah River site, Georgia and South Carolina, 2002, and for selected ground-water management scenarios, 2002 and 2020

Ground-water flow under 2002 hydrologic conditions was evaluated in an eight-county area in Georgia and South Carolina near the Savannah River Site (SRS), by updating boundary conditions and pumping rates in an existing U.S. Geological Survey (USGS) ground-water model. The original ground-water model, developed to simulate hydrologic conditions during 1987-92, used the quasi-three-dimensional approach by dividing the Floridan, Dublin, and Midville aquifer systems into seven aquifers. The hydrogeologic system was modeled using six active layers (A2-A7) that were separated by confining units with an overlying source-sink layer to simulate the unconfined Upper Three Runs aquifer (layer A1). Potentiometric- surface maps depicting September 2002 for major aquifers were used to update, evaluate, and modify boundary conditions used by the earlier ground-water flow model. The model was updated using the USGS finite-difference code MODFLOW-2000 for mean-annual conditions during 1987-92 and 2002. The specified heads in the source-sink layer A1 were lowered to reflect observed water-level declines during the 1998-2002 drought. These declines resulted in a decrease of 12.1 million gallons per day (Mgal/d) in simulated recharge or vertical inflow to the uppermost confined aquifer (Gordon, layer A2). Although ground-water pumpage in the study area has increased by 32 Mgal/d since 1995, most of this increase (17.5 Mgal/d) was from the unconfined Upper Three Runs aquifer (source-sink layer A1) with the remaining 14.5 Mgal/d assigned to the active layers within the model (A2-A7). The simulated water budget for 2002 shows a decrease from the 1987-92 model from 1,040 Mgal/d to 1,035 Mgal/d. The decreased ground-water inflows and increased ground-water withdrawal rates reduced the simulated ground-water outflow to river cells in the active layers of the model by 43 Mgal/d. The calibration statistics for all layers of the 2002 simulation resulted in a decrease in the root mean square (RMS) of the residuals from 10.6 to 8.0 feet (ft). The residuals indicate 83.3 percent of the values for the 2002 simulation met the calibration error criteria established in the original model, whereas 88.8 percent was within the specified range for the 1987-92 simulation. Simulated ground-water outflow to the Savannah River and its tributaries during water year 2002 was 560 cubic feet per second (ft3/s), or 86 percent of the observed gain in mean-annual streamflow between streamflow gaging stations at the Millhaven, Ga., and Augusta, Ga. At Upper Three Runs Creek, simulated ground-water discharge during 2002 was 110 ft3/s, or 83 percent of the observed streamflow at two streamflow gaging stations near the SRS. These results indicate that the constructed model calibrated to 1987-92 conditions and modified for 2002 dry conditions is still representative of the hydrologic system. The USGS particle-tracking code MODPATH was used to generate advective water-particle pathlines and their associated time-of-travel based on MODFLOW simulations for 1987-92, 2002, and each of four hypothetical ground-water management scenarios. The four hypothetical ground-water management scenarios represent hydrologic conditions for (1) reported pumping for 2002 and boundary conditions for an average year; (2) reported pumping for 2002 with SRS pumping discontinued and boundary conditions for an average year; (3) projected 2020 pumping and boundary conditions for an average year; and (4) projected 2020 pumping and boundary conditions for a dry year. The MODPATH code was used in forward-tracking mode to evaluate flowpaths from areas on the SRS and in backtracking mode to evaluate further areas of previously documented trans-river flow on the Georgia side of the Savannah River. Trans-river flow is a condition in which the local head gradients might allow migration of contaminants from the SRS into the underlying aquifers and beneath the Savannah River into Georgia. More...

Georgia, South Carolina↗
USGS2018

Evaluation of the Source and Transport of High Nitrate Concentrations in Ground Water, Warren Subbasin, California

Ground water historically has been the sole source of water supply for the Town of Yucca Valley in the Warren subbasin of the Morongo ground-water basin, California. An imbalance between ground-water recharge and pumpage caused ground-water levels in the subbasin to decline by as much as 300 feet from the late 1940s through 1994. In response, the local water district, Hi-Desert Water District, instituted an artificial recharge program in February 1995 using imported surface water to replenish the ground water. The artificial recharge program resulted in water-level recoveries of as much as 250 feet in the vicinity of the recharge ponds between February 1995 and December 2001; however, nitrate concentrations in some wells also increased from a background concentration of 10 milligrams per liter to more than the U.S. Environmental Protection Agency (USEPA) maximum contaminant level (MCL) of 44 milligrams per liter (10 milligrams per liter as nitrogen). The objectives of this study were to: (1) evaluate the sources of the high-nitrate concentrations that occurred after the start of the artificial-recharge program, (2) develop a ground-water flow and solute-transport model to better understand the source and transport of nitrates in the aquifer system, and (3) utilize the calibrated models to evaluate the possible effect of a proposed conjunctive-use project. These objectives were accomplished by collecting water-level and water-quality data for the subbasin and assessing changes that have occurred since artificial recharge began. Collected data were used to calibrate the ground-water flow and solute-transport models. Data collected for this study indicate that the areal extent of the water-bearing deposits is much smaller (about 5.5 square miles versus 19 square miles) than that of the subbasin. These water-bearing deposits are referred to in this report as the Warren ground-water basin. Faults separate the ground-water basin into five hydrogeologic units: the west, the midwest, the mideast, the east and the northeast hydrogeologic units. Water-quality analyses indicate that septage from septic tanks is the primary source of the high-nitrate concentrations measured in the Warren ground-water basin. Water-quality and stable-isotope data, collected after the start of the artificial recharge program, indicate that mixing occurs between imported water and native ground water, with the highest recorded nitrate concentrations in the midwest and the mideast hydrogeologic units. In general, the timing of the increase in measured nitrate concentrations in the midwest hydrogeologic unit is directly related to the distance of the monitoring well from a recharge site, indicating that the increase in nitrate concentrations is related to the artificial recharge program. Nitrate-to-chloride and nitrogen-isotope data indicate that septage is the source of the measured increase in nitrate concentrations in the midwest and the mideast hydrogeologic units. Samples from four wells in the Warren ground-water basin were analyzed for caffeine and selected human pharmaceutical products; these analyses suggest that septage is reaching the water table. There are two possible conceptual models that explain how high-nitrate septage reaches the water table: (1) the continued downward migration of septage through the unsaturated zone to the water table and (2) rising water levels, a result of the artificial recharge program, entraining septage in the unsaturated zone. The observations that nitrate concentrations increase in ground-water samples from wells soon after the start of the artificial recharge program in 1995 and that the largest increase in nitrate concentrations occur in the midwest and mideast hydrogeologic units where the largest increase in water levels occur indicate the validity of the second conceptual model (rising water levels). The potential nitrate concentration resulting from a water-level rise in the midwest and mideast hydrogeologic units was estimated using a simple mixing-cell model. The estimated value is within the range of concentrations measured in samples from wells, further indicating the validity of the second conceptual model. A ground-water flow model and a solute-transport model were developed for the Warren ground-water basin for the period 1956-2001. MODFLOW-96 was used for the ground-water flow model and MOC3D was used for the solute-transport model. The model cell size is about 500 feet by 500 feet and the models were discretized vertically into three layers. The models were calibrated using a trial-and-error approach using water-level and nitrate-concentration data collected between 1956-2001. In order to better match the measured data, low fault hydraulic characteristic values were required, thereby compartmentalizing the ground-water basin. In addition, it was necessary to parameterize the specific yield distribution for the top model layer where unconfined ground-water conditions occur into three homogeneous zones. Separate sets of specific- yield values were needed to simulate the drawdown and subsequent water-level recovery. In addition, the calibrated natural recharge was about 83 acre-feet per year. The entrainment of unsaturated-zone septage was simulated as recharge having an associated nitrate concentration. The volume of recharge was a function of the measured water-level rise between 1994-98 and the moisture content of the unsaturated zone. The nitrate concentration of the recharge water was a weighted function of the assumed nitrate concentration in the infiltrating water associated with the overlying land use. The simulated hydraulic head and nitrate concentration results were in good agreement with the measured data indicating that the mechanism for the increase in nitrate concentrations was rising water levels entraining high-nitrate septage in the unsaturated-zone. The calibrated models were used to simulate the possible effects of a planned conjunctive-use project in the western part of the ground-water basin. The simulated project included the addition of a new recharge pond and a new extraction well. In addition, recharge at two existing recharge ponds was increased and three existing production wells were pumped, treated in a nitrate-removal facility, and used for water supply. The simulated hydraulic heads increased in the west, the mideast, and parts of the east hydrogeologic units; however, the simulated hydraulic heads decreased in the midwest and northeast hydrogeologic units. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter (10 milligrams per liter as nitrogen) in parts of the west as a result of the increase in simulated hydraulic head. The simulated nitrate concentrations decreased in part of the midwest hydrogeologic unit as a result of the artificial recharge and pumping from the nitrate-removal wells. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter in part of the mideast and parts of the east hydrogeologic units beneath commercial land-use areas.

California↗