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 973 records · Page 54Linked to original sources
USGS1985

Seismic-refraction study of suspected drift-filled bedrock valleys in Ramsey County, Minnesota

Seismic-refraction surveys were made across suspected buried, drift-filled bedrock valleys believed to underlie two sites of known ground-water contamination the Twin Cities Army Ammunition Plant (TCAAP) near New Brighton and the former Koppers Coke Plant in St. Paul, Ramsey County, Minnesota. Refraction data were collected along two lines at each site; each line traversed the axis of a suspected valley. Drift-filled bedrock valleys were thought to incise the Prairie du Chien-Jordan aquifer to an altitude between 500 and 560 feet above sea level at the TCAAP site. The interpretation of one seismic profile indicates that a valley probably does not exist under the area surveyed; the modeled depth profile shows the bedrock surface ranging between altitudes of about 700 to 780 feet. Interpretation of a second seismic profile just south of the TCAAP indicates that a shallow valley may exist under the northern half of the line; the modeled depth profile shows that the bedrock surface declines from an altitude of about 780 feet to about 690 feet over a horizontal distance of 400 feet. A drift-filled bedrock valley was thought to incise the St. Peter aquifer to an altitude between 770 and 800 feet above sea level at the Koppers site. The interpretation of a seismic profile just east of the Koppers site is not conclusive, but suggests that a bedrock valley may be present near the middle of the line. The interpretation of a second seismic profile across the westward extension of the same suspected valley also is not conclusive, but suggests that a bedrock valley may be present at the north end of the line. The optimal field layout for each line at the site (longer shot offsets) could not be obtained because of limited space available in the densely developed residential neighborhoods.

Minnesota↗
USGS1981

Effect of snowmelt on the quality of Filson Creek and Omaday Lake, northeastern Minnesota

Concentrations of major constituents were determined in the surface water, ground water, and precipitation in Filson Creek watershed to evaluate the effects of acid precipitation on surface-water quality during snowmelt. Concentrations of sulfate increased in Filson Creek and Omaday Lake during snowmelt from less than 2 to 12 milligrams per liter in 1977 and from less than 2 to 4 milligrams per liter in 1979. Hydrogen-ion activity increased from 10 -6.5 to 10 -5.0 during snowmelt in 1977 and from 10 -6.5 to 10 -6.0 in 1979. The magnitude of the increases of sulfate and hydrogen-ion concentrations depend on the (1) timing of the snowmelt, (2) concurrent precipitation, and (3) antecedent soil-moisture conditions. Alkalinity and concentrations of total calcium, magnesium, and sodium in the creek during snowmelt reflect the simple dilution of streamflow with more dilute precipitation. The mixing process is defined by hyperbolic curves of concentration versus discharge and can be modeled by use of a mass-balance equation. Concentrations of silica and total potassium in the creek during snowmelt are less than that at base flow because of dilution and uptake by diatoms and plants, respectively. A quality-assurance study was made to determine the accuracy and precision of the turbidimetric, complexometric (methylthymol blue), and ion-chromatographic methods used to determine concentrations of sulfate in water. Results suggest that the ion-chromatographic method is more accurate and more precise in the detection of sulfate at concentrations of 5 milligrams per liter or less. To obtain accuracy with either the turbidimetric or the complexometric methylthymol blue method, color blanks and turbidity should be determined if concentrations of sulfate are expected to be low, such as in precipitation, or where water is colored by organic material. Interference by color and turbidity is not a factor in the ion-chromatographic method.

Minnesota↗
USGS2001

Hydrogeology and ground-water flow in the Memphis and Fort Pillow aquifers in the Memphis area, Tennessee

On the basis of known hydrogeology of the Memphis and Fort Pillow aquifers in the Memphis area, a three-layer, finite-difference numerical model was constructed and calibrated as the primary tool to refine understanding of flow in the aquifers. The model was calibrated and tested for accuracy in simulating measured heads for nine periods of transient flow from 1886-1985. Testing and sensitivity analyses indicated that the model accurately simulated observed heads areally as well as through time. The study indicates that the flow system is currently dominated by the distribution of pumping in relation to the distribution of areally variable confining units. Current withdrawal of about 200 million gallons per day has altered the prepumping flow paths, and effectively captured most of the water flowing through the aquifers. Ground-water flow is controlled by the altitude and location of sources of recharge and discharge, and by the hydraulic characteristics of the hydrogeologic units. Leakage between the Fort Pillow aquifer and Memphis aquifer, and between the Memphis aquifer and the water-table aquifers (alluvium and fluvial deposits) is a major component of the hydrologic budget. The study indicates that more than 50 percent of the water withdrawn from the Memphis aquifer in 1980 is derived from vertical leakage across confining units, and the leakage from the shallow aquifer (potential source of contamination) is not uniformly distributed. Simulated leakage was concentrated along the upper reaches of the Wolf and Loosahatchie Rivers, along the upper reaches of Nonconnah Creek, and the surficial aquifer of the Mississippi River alluvial plain. These simulations are supported by the geologic and geophysical evidence suggesting relatively thin or sandy confining units in these general locations. Because water from surficial aquifers is inferior in quality and more susceptible to contamination than water in the deeper aquifers, high rates of leakage to the Memphis aquifer may be cause for concern. A significant component of flow (12 percent) discharging from the Fort Pillow aquifer was calculated as upward leakage to the Memphis aquifer. This upward leakage was generally limited to areas near major pumping centers in the Memphis aquifer, where heads in the Memphis aquifer have been drawn significantly below heads in the Fort Pillow aquifer. Although the Fort Pillow aquifer is not capable of producing as much water as the Memphis aquifer for similar conditions, it is nonetheless a valuable resource throughout the area.

Water-Resources Investigations Report↗
USGS2006

Physical and hydrochemical evidence of lake leakage near Jim Woodruff Lock and Dam and of ground-water inflow to Lake Seminole, and an assessment of karst features in and near the lake, southwestern Georgia and northwestern Florida

Hydrogeologic data and water-chemistry analyses indicate that Lake Seminole leaks into the Upper Floridan aquifer near Jim Woodruff Lock and Dam, southwestern Georgia and northwestern Florida, and that ground water enters Lake Seminole along upstream reaches of the lake’s four impoundment arms (Chattahoochee and Flint Rivers, Spring Creek, and Fishpond Drain). Written accounts by U.S. Army Corps of Engineers geologists during dam construction in the late 1940s and early 1950s, and construction-era photographs, document karst-solution features in the limestone that comprise the lake bottom and foundation rock to the dam, and confirm the hydraulic connection of the lake and aquifer. More than 250 karst features having the potential to connect the lake and aquifer were identified from preimpoundment aerial photographs taken during construction. An interactive map containing a photomosaic of 53 photographic negatives was orthorectfied to digital images of 1:24,000-scale topographic maps to aid in identifying karst features that function or have the potential to function as locations of water exchange between Lake Seminole and the Upper Floridan aquifer. Some identified karst features coincide with locations of mapped springs, spring runs, and depressions that are consistent with sinkholes and sinkhole ponds. Hydrographic surveys using a multibeam echosounder (sonar) with sidescan sonar identified sinkholes in the lake bottom along the western lakeshore and in front of the dam. Dye-tracing experiments indicate that lake water enters these sinkholes and is transported through the Upper Floridan aquifer around the west side of the dam at velocities of about 500 feet per hour to locations where water "boils up" on land (at Polk Lake Spring) and in the channel bottom of the Apalachicola River (at the "River Boil"). Water discharging from Polk Lake Spring joins flow from a spring-fed ground-water discharge zone located downstream of the dam; the combined flow disappears into a sinkhole located on the western floodplain of the river and is transmitted through the Upper Floridan aquifer, eventually discharging to the Apalachicola River at the River Boil. Acoustic Doppler current profiling yielded flow estimates from the River Boil in the range from about 140 to 220 cubic feet per second, which represents from about 1 to 3 percent of the average daily flow in the river. Binary mixing-model analysis using naturally occurring isotopes of oxygen and hydrogen (oxygen-18 and deuterium) indicates that discharge from the River Boil consists of a 13-to-1 ratio of lake water to ground water and that other sources of lake leakage and discharge to the boil probably exist. Analyses of major ions, nutrients, radon-222, and stable isotopes of hydrogen and oxygen contained in water samples collected from 29 wells, 7 lake locations, and 5 springs in the Lake Seminole area during 2000 indicate distinct chemical signatures for ground water and surface water. Ground-water samples contained higher concentrations of calcium and magnesium, and higher alkalinity and specific conductance than surface-water samples, which contained relatively high concentrations of total organic carbon and sulfate. Solute and isotopic tracers indicate that, from May to October 2000, springflow exhibited more ground-water qualities (high specific conductance, low dissolved oxygen, and low temperature) than surface water; however, the ratio of ground water to surface water of the springs was difficult to quantify from November to April because of reduced springflow and rapid mixing of springflow and lake water during sampling. The saturation index of calcite in surface-water samples indicates that while surface water is predominately undersaturated with regard to calcite year-round, a higher potential for dissolution of the limestone matrix exists from late fall through early spring than during summer. The relatively short residence time (5–7 hours) and rapid flow velocity (nearly 500 feet per hour) of lake water leaking into the Upper Floridan aquifer and exiting at the River Boil in the Apalachicola River implies that calcite-undersaturated water is in constant contact with the limestone, increasing the potential for limestone dissolution and enlargement of flow pathways by erosion. A relatively low potential exists, however, for limestone dissolution to cause sudden sinkhole collapse followed by catastrophic lake drainage because ground-water levels close to the lake, except near the dam, are nearly the same as lake stage, resulting in low vertical and lateral hydraulic gradients and low flow between the lake and aquifer. An increased potential for lake leakage and sinkhole formation and collapse exists near some in-lake springs during colder months of the year, as density differences and the hydraulic potential between lake water and ground water establish the conditions for calcite-undersaturated lake water to enter nonflowing springs and contact limestone.

Florida, Georgia↗
USGS2013

Method to support Total Maximum Daily Load development using hydrologic alteration as a surrogate to address aquatic life impairment in New Jersey streams

More than 300 ambient monitoring sites in New Jersey have been identified by the New Jersey Department of Environmental Protection (NJDEP) in its integrated water-quality monitoring and assessment report (that is, the 305(b) Report on general water quality and 303(d) List of waters that do not support their designated uses) as being impaired with respect to aquatic life; however, no unambiguous stressors (for example, nutrients or bacteria) have been identified. Because of the indeterminate nature of the broad range of possible impairments, surrogate measures that more holistically encapsulate the full suite of potential environmental stressors need to be developed. Streamflow alteration resulting from anthropogenic changes in the landscape is one such surrogate. For example, increases in impervious surface cover (ISC) commonly cause increases in surface runoff, which can result in “flashy” hydrology and other changes in the stream corridor that are associated with streamflow alteration. The NJDEP has indicated that methodologies to support a hydrologically based Total Maximum Daily Load (hydro-TMDL) need to be developed in order to identify hydrologic targets that represent a minimal percent deviation from a baseline condition (“minimally altered”) as a surrogate measure to meet criteria in support of designated uses. The primary objective of this study was to develop an applicable hydro-TMDL approach to address aquatic-life impairments associated with hydrologic alteration for New Jersey streams. The U.S. Geological Survey, in cooperation with the NJDEP, identified 51 non- to moderately impaired gaged streamflow sites in the Raritan River Basin for evaluation. Quantile regression (QR) analysis was used to compare flow and precipitation records and identify baseline hydrographs at 37 of these sites. At sites without an appropriately long period of record (POR) or where a baseline hydrograph could not be identified with QR, a rainfall-runoff model was used to develop simulated baseline hydrographs. The hydro-TMDL approach provided an opportunity to evaluate proportional differences in flow attributes between observed and baseline hydrographs and to develop complementary flow-ecology response relations at a subset of Raritan River Basin sites where available flow and ecological information overlapped. The New Jersey Stream Classification Tool (NJSCT) was used to determine the stream class of all 51 study sites by using either an observed or a simulated baseline hydrograph. Two New Jersey stream classes (A and C) were evaluated to help characterize the unique hydrology of the Raritan River Basin. In general, class C streams (1.99–40.7 square miles) had smaller drainage areas than class A streams (0.7–785 square miles). Many of the non-impaired and moderately impaired class A and C streams in the Raritan River Basin were found to have significant hydrologic alteration as indicated by numerous flow values that fell outside the established 25th-to-75th- and the more conservative 40th-to-60th-percentile boundaries. However, percent deviations for the class C streams (defined as moderately stable streams with moderately high base-flow contributions) were, in general, much larger than those for the class A streams (defined as semiflashy streams characterized by moderately low base flow). The greater deviations for class C streams in the hydro-TMDL assessments likely resulted from comparisons that were based solely on simulated baseline hydrographs, which were developed without considering any anthropogenic influences in the basin. In contrast, comparisons for many of the class A streams were made by using an observed baseline, which already includes an implicit level of ISC and other human influences on the landscape. By using the hydro-TMDL approach, numerous flow deviations were identified that were indicative of streams that are highly regulated by reservoirs or dams, streams that are affected by increasing amounts of surface runoff resulting from ISC, and streams that are affected by water abstraction (that is, groundwater or surface-water withdrawals used for agricultural and human supply). Eight of the reservoir- and (or) dam-affected sites showed flow deviations that are indicative of flow-managed systems. For example, indices that account for the timing and magnitude of high and low flows were often found to fall outside the 25th-to-75th-percentile range. In general, at regulated class C streams, annual summer low flows are arriving later and tend to be lower, and high flows are arriving earlier with higher magnitudes of longer duration. At class A streams, high and low flows are arriving later with an overall increase in discharge with respect to the prereservoir baseline conditions. The drainage basins of eight of the study sites had large values of ISC (>10 percent), most likely as a result of expanding urban development. In general, the magnitude and frequency of high flows at class A and C sites with high ISC are increasing and were commonly found to fall outside the 25th-to-75th-percentile range. Additionally, magnitudes of low flows are becoming lower and, although the timing of high flows was highly variable, low-flow events appeared to be arriving earlier than would be expected under normal low-flow conditions. Three of the study sites appeared to be affected by hydrologic changes associated with water abstraction. At these sites, the timing of flows appeared to be altered. For example, low flows tended to arrive earlier and high flows arrived later at two of the three sites. Additionally, the magnitude and duration of low flows were commonly less than the 25th-percentile value and the duration of high flows appeared to increase. A reduced set of hydrologic and ecological variables was used to develop univariate and multivariate flow-ecology response models for the aquatic-invertebrate assemblage. Many hydrologic variables accounting for the duration, magnitude, frequency, and timing of flows were significantly correlated with ecological response. Multiple linear regression (MLR) models were developed to provide a more holistic evaluation of the combined effects of hydrologic alteration and to identify models with two or three hydrologic variables that account for a significant proportion of the variability in invertebrate-assemblage condition as represented by assemblage metric scores. MLR models, derived on the basis of hydrologic attributes, accounted for 35 to 75 percent of the variability in assemblage condition. The hydro-TMDL method developed herein for non- to moderately impaired Raritan River Basin streams utilizes a “surrogate” approach in place of the traditional “pollutant of concern” approach commonly used for TMDL development. Managers can use the results obtained by using the hydro-TMDL method to offset the effects of impervious-surface runoff and altered streamflow and to implement measures designed to achieve the necessary load reductions for the “pollutant of concern” (that is, percentage deviations of stream-class-specific flow-index values outside the established 25th-to-75th-percentile range). In this case, such deviations could represent all or a subset of the altered flow indices that prevent the stream from meeting designated aquatic-life criteria. This hydro-TMDL uses a reference, or attainment stream approach for developing the TMDL endpoint. That is, either observed or simulated baseline hydrographs were selected as appropriate reference conditions on the basis of results of QR analysis and watershed modeling procedures, respectively. For any stream in the Raritan River Basin evaluated as part of this study, the hydro-TMDL can be expressed as the greatest amount of deviation in flow a stream can exhibit without violating the stream’s designated aquatic-life criteria. Use of this surrogate approach is appropriate because flows that fall outside the established percentile ranges are ultimately a function of many anthropogenic modifications of the landscape, including the amount of stormwater runoff generated from impervious surfaces within a given basin, the presence of manmade structures designed to retain or divert water, the magnitude of ground- and surface-water abstraction, and the presence of water-supply processes implemented to support human needs. In addition, the stream-type-specific flow indices used as the basis for the hydro-TMDL approach are useful for representing the hydrologic conditions of class A and C streams/basins because they incorporate the full spectrum of flow conditions (very low to very high) that occur in the stream system over a long period of time, as well as those flow properties that change as a result of seasonal variation. Ultimately, an estimate of the maximum percentage flow reduction that could be allowed will be needed to address the aquatic-life impairments in many of the study streams in the Raritan River Basin and will be necessary for identifying appropriate target flow conditions for hydro-TMDL implementation. As described in this report, a target flow value equal to the 25th- or 75th-percentile flow rate could be selected as the point useful for setting specific hydrologic targets. This selection, however, is a management decision that could vary depending on the designated use of the stream or other regulatory factors (for example, water-supply protection, trout production, antidegradation policies, or special protection designations). In New Jersey streams where no unambiguous stressors can be identified, State monitoring agencies, such as the NJDEP, could choose to require the implementation of a flow-based TMDL that not only supports designated uses, but meets the regulatory requirements under the Clean Water Act, and represents a balance between water supply intended to meet human needs and the conservation of ecosystem integrity.

New Jersey↗
USGS2005

Simulated changes in water levels caused by potential changes in pumping from shallow aquifers of Virginia Beach, Virginia

A steady-state ground-water flow model of the southern watersheds of Virginia Beach, Virginia, was refined and used to simulate changes in aquifer water levels caused by potential changes in pumping in the Transition Area of Virginia Beach, Va., a 20-square mile planning zone that runs through the middle of the city. Cessation of dewatering at borrow pits, pumping to irrigate a golf course, pumping to irrigate lawns of a hypothetical neighborhood, and pumping to irrigate both the golf course and lawns of the hypothetical neighborhood were simulated. Simulated recoveries from cessation of dewatering of borrow pits were generally restricted to the immediate area of the pits. The simulated recoveries averaged about 20 feet (ft) near the center of the cells representing the active areas of the pits and 2 ft at the cells representing the extent of the pits. At a golf course, 4 hypothetical wells pumping 300,000 gallons per day (gal/d) from the Yorktown sand aquifer resulted in drawdowns averaging 10 ft in the pumping cells and 1 ft at a distance of 1.2 miles (mi) from the center of the pumping cells. The extent of the 1-ft drawdown was virtually the same as that simulated previously and reported in a permit application for the golf course. Simulated pumping of 150,000 gal/d from 4 cells in the confined sand aquifer representing a 40-acre neighborhood resulted in drawdowns averaging 7 ft in the pumping cells and 1 ft at a distance of 0.8 mi from the center of the cells. Simulated pumping of 300,000 gal/d from the same 4 cells resulted in drawdowns averaging 15 ft in the pumping cells and 1 ft at a distance of 1.4 mi from the center of the cells. Simulated pumping of 150,000 gal/d at the golf course and another 150,000 gal/d in the hypothetical neighborhood resulted in drawdowns that averaged 5 ft around the cells representing the golf course wells spaced 1,300 ft apart and 7 ft around the contiguous cells representing the 40-acre neighborhood. A drawdown of 1 ft encompassed most of the eastern half of the Transition Area.

Virginia↗
USGS2004

PHAST--a program for simulating ground-water flow, solute transport, and multicomponent geochemical reactions

The computer program PHAST simulates multi-component, reactive solute transport in three-dimensional saturated ground-water flow systems. PHAST is a versatile ground-water flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. PHAST is applicable to the study of natural and contaminated ground-water systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock-water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, density-dependent flow, or waters with high ionic strengths. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux, and leaky conditions, as well as the special cases of rivers and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, gases, surface complexation sites, ion exchange sites, and solid solutions; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, gases, exchangers, surfaces, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and post-processing programs; or in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST. Operator splitting of the flow, transport, and geochemical equations is used to separate the three processes into three sequential calculations. No iterations between transport and reaction calculations are implemented. A three-dimensional Cartesian coordinate system and finite-difference techniques are used for the spatial and temporal discretization of the flow and transport equations. The non-linear chemical equilibrium equations are solved by a Newton-Raphson method, and the kinetic reaction equations are solved by a Runge-Kutta or an implicit method for integrating ordinary differential equations. The PHAST simulator may require large amounts of memory and long Central Processing Unit (CPU) times. To reduce the long CPU times, a parallel version of PHAST has been developed that runs on a multiprocessor computer or on a collection of computers that are networked. The parallel version requires Message Passing Interface, which is currently (2004) freely available. The parallel version is effective in reducing simulation times. This report documents the use of the PHAST simulator, including running the simulator, preparing the input files, selecting the output files, and visualizing the results. It also presents four examples that verify the numerical method and demonstrate the capabilities of the simulator. PHAST requires three input files. Only the flow and transport file is described in detail in this report. The other two files, the chemistry data file and the database file, are identical to PHREEQC files and the detailed description of these files is found in the PHREEQC documentation.

Techniques and Methods↗
USGS2000

Factors affecting nutrient trends in major rivers of the Chesapeake Bay Watershed

Trends in nutrient loads and flow-adjusted concentrations in the major rivers entering Chesapeake Bay were computed on the basis of water-quality data collected between 1985 and 1998 at 29 monitoring stations in the Susquehanna, Potomac, James, Rappahannock, York, Patuxent, and Choptank River Basins. Two computer models-the Chesapeake Bay Watershed Model (WSM) and the U.S. Geological Survey's 'Spatially Referenced Regressions on Watershed attributes' (SPARROW) Model-were used to help explain the major factors affecting the trends. Results from WSM simulations provided information on temporal changes in contributions from major nutrient sources, and results from SPARROW model simulations provided spatial detail on the distribution of nutrient yields in these basins. Additional data on nutrient sources, basin characteristics, implementation of management practices, and ground-water inputs to surface water were analyzed to help explain the trends. The major factors affecting the trends were changes in nutrient sources and natural variations in streamflow. The dominant source of nitrogen and phosphorus from 1985 to 1998 in six of the seven tributary basins to Chesapeake Bay was determined to be agriculture. Because of the predominance of agricultural inputs, changes in agricultural nutrient sources such as manure and fertilizer, combined with decreases in agricultural acreage and implementation of best management practices (BMPs), had the greatest impact on the trends in flow-adjusted nutrient concentrations. Urban acreage and population, however, were noted to be increasing throughout the Chesapeake Bay Watershed, and as a result, delivered loads of nutrients from urban areas increased during the study period. Overall, agricultural nutrient management, in combination with load decreases from point sources due to facility upgrades and the phosphate detergent ban, led to downward trends in flow-adjusted nutrient concentrations atmany of the monitoring stations in the watershed. The loads of nutrients, however, were not reduced significantly at most of the monitoring stations. This is due primarily to higher streamflow in the latter years of the monitoring period, which led to higher loading in those years. Results of this study indicate a need for more detailed information on BMP effectiveness under a full range of hydrologic conditions and in different areas of the watershed; an internally consistent fertilizer data set; greater consideration of the effects of watershed processes on nutrient transport; a refinement of current modeling efforts; and an expansion of the non-tidal monitoring network in the Chesapeake Bay Watershed.

Delaware, Maryland, New York, Pennsylvania, Virgin↗
USGS1995

Preliminary assessment of injection, storage, and recovery of freshwater in the lower Hawthorn aquifer, Cape Coral, Florida

A preliminary assessment of subsurface injection, storage and recovery of fresh canal water was made in the naturally brackish lower Hawthorn aquifer in Cape Coral, southwestern Florida. A digital modeling approach was used for this preliminary assessment, incorporating available data on hydrologic conditions, aquifer properties, and water quality to simulate density-dependent ground-water flow and advective-dispersive transport of a conservative ground-water solute (chloride ion). A baseline simulation was used as reference to compare the effects of changing various operational factors on the recovery efficiency. A recovery efficiency of 64 percent was estimated for the baseline simulation. Based on the model, the recovery efficiency increases if the injection rate and recovery rates are increased and if the ratio of recovery rate to injection rate is increased. Recovery efficiency decreases if the amount of water injected is increased; slightly decreases if the storage time is increased; is not changed significantly if the water is injected to a specific flow zone; increases with successive cycles of injection, storage, and recovery; and decreases if the chloride concentrations in either the injection water or native aquifer water are increased. In everal hypothetical tests, the recovery efficiency fluctuated between 22 and about 100 percent. Two successive cycles could bring the recovery efficiency from 60 to about 80 percent. Interlayer solute mass movement across the upper and lower boundaries seems to be the most important factor affecting the recovery efficiency. A sensitivity analysis was performed applying a technique in which the change in the various factors and the corresponding model responses are normalized so that meaningful comparisons among the responses could be made. The general results from the sensitivity analysis indicated that the permeabilities of the upper and lower flow zones were the most important factors that produced the greatest changes in the relative sensitivity of the recovery efficiency. Almost equally significant changes occurred in the relative sensitivity of the recovery efficiency when all porosity values of the upper and lower flow zones and the leaky confining units and the vertical anisotropy ratio were changed. The advective factors are the most important in the Cape Coral area according to the sensitivity analysis. However, the dispersivity values used in the model were extrapolated from studies conducted at the nearby Lee County Water Treatment Plant, and these values might not be representative of the actual dispersive characteristics of the lower Hawthorn aquifer in the Cape Coral area.

Florida↗
USGS1987

HST3D; a computer code for simulation of heat and solute transport in three-dimensional ground-water flow systems

The Heat- and Soil-Transport Program (HST3D) simulates groundwater flow and associated heat and solute transport in three dimensions. The three governing equations are coupled through the interstitial pore velocity, the dependence of the fluid density on pressure, temperature, the solute-mass fraction , and the dependence of the fluid viscosity on temperature and solute-mass fraction. The solute transport equation is for only a single, solute species with possible linear equilibrium sorption and linear decay. Finite difference techniques are used to discretize the governing equations using a point-distributed grid. The flow-, heat- and solute-transport equations are solved , in turn, after a particle Gauss-reduction scheme is used to modify them. The modified equations are more tightly coupled and have better stability for the numerical solutions. The basic source-sink term represents wells. A complex well flow model may be used to simulate specified flow rate and pressure conditions at the land surface or within the aquifer, with or without pressure and flow rate constraints. Boundary condition types offered include specified value, specified flux, leakage, heat conduction, and approximate free surface, and two types of aquifer influence functions. All boundary conditions can be functions of time. Two techniques are available for solution of the finite difference matrix equations. One technique is a direct-elimination solver, using equations reordered by alternating diagonal planes. The other technique is an iterative solver, using two-line successive over-relaxation. A restart option is available for storing intermediate results and restarting the simulation at an intermediate time with modified boundary conditions. This feature also can be used as protection against computer system failure. Data input and output may be in metric (SI) units or inch-pound units. Output may include tables of dependent variables and parameters, zoned-contour maps, and plots of the dependent variables versus time. (Lantz-PTT)

Water-Resources Investigations Report↗
USGS2005

Water-quality assessment of part of the Upper Mississippi River Basin, Minnesota and Wisconsin — Ground-water quality along a flow system in the Twin Cities metropolitan area, Minnesota, 1997-98

As part of a national analysis of the effects of land use on ground-water quality for the National Water-Quality Assessment Program, the U.S. Geological Survey sampled wells along a flow system in surficial glacial aquifers in the northwestern part of the Twin Cities metropolitan area during 1997 and 1998. In addition, a reconnaissance steady-state ground-water model was developed to estimate flowpaths and dates of ground-water recharge using a particle-tracking routine. Sediment samples collected during drilling had high horizontal hydraulic conductivities (ranging from about 3.1 to about 190 feet per day, based on sediment-size analysis of well cuttings), small organic carbon concentrations (ranging from less than 0.2 to 160 grams per kilogram), and pH values that were mostly alkaline (ranging from 4.9 to 8.2). Water samples were analyzed for physical properties, major ions, iron, manganese, nutrients, organic carbon, radon, pesticides, volatile organic compounds (VOCs), chlorofluorocarbons, tritium, and isotopes of nitrogen, hydrogen, and oxygen. Most of the water samples had small dissolved-oxygen concentrations (less than 1 milligram per liter). Calcium, magnesium, sodium, bicarbonate, sulfate, and chloride were the primary dissolved constituents in water samples. Nitrite plus nitrate as nitrogen (nitrate) concentrations were less than the U.S. Environmental Protection Agency (USEPA) Maximum Contaminant Level of 10 mg/L. Nitrogen isotope ratios indicated that the sources of nitrate primarily were soils, animal waste, or denitrification that increased nitrogen isotope ratios in nitrate. Small concentrations of pesticides were detected in the shallow parts of the aquifer. The herbicide prometon was the most frequently detected pesticide. Herbicides applied to control grasses and weeds in corn (atrazine, simazine, and metolachlor) also were frequently detected in water samples. All pesticide and VOCs detected were below USEPA Maximum Contaminant Levels or Health Advisory Limits. Chlorofluorocarbon compounds and tritium concentrations were used to estimate dates of recharge of ground-water samples. In general, shallower ground-water samples were more recently recharged although most water sampled from the aquifer was recharged after 1955. Although land use had substantial effects on ground-water quality, the distribution of contaminants in the aquifer also is affected by complex combinations of factors and processes that include sources of natural and anthropogenic contaminants, three-dimensional advective flow, physical and hydrologic settings, age and evolution of ground water, and transformation of chemical compounds along the flow system. Compounds such as nitrate and dissolved oxygen were greatest in water samples from the upgradient end of the flow system and near the water table. Specific conductance and dissolved solids increased along the flow system and with depth due to increase in residence time in the flow system and dissolution of aquifer materials.

Minnesota↗
USGS1990

Effects of storm-water runoff on local ground-water quality, Clarksville, Tennessee

Storm-related water-quality data were collected at a drainage-well site and at a spring site in Clarksville, Tennessee, to define the effects of storm-water runoff on the quality of ground water in the area. A dye-trace test verified the direct hydraulic connection between the drainage well and Mobley Spring. Samples of storm run off and spring flow were collected at these sites for nine storms during the period February to October 1988. Water samples were collected also from Mobley Spring and two other springs and two observation wells in the area during dry-weather conditions to assess the general quality of ground water in an urban karst terrain. Evaluation of the effect of storm-water runoff on the quality of local ground water is complicated by the presence of other sources of contaminants in the area Concentrations and load for most major constituents were much smaller in storm-water runoff at the drainage well than in the discharge of Mobley Spring, indicating that much of the chemical constituent load discharged from the spring comes from sources other than the drainage well. However, for some of the minor constituents associated with roadway runoff (arsenic, copper, lead, organic carbon, and oil and grease), the drainage well contributed relatively large amounts of these constituents to local ground water during storms. The close correlation between concentrations of total organic carbon and concentrations of most trace metals at the drainage-well and Mobley Spring sites indicates that these constituents are transported together. Many trace metals were flushed early during each runoff event. Mean storm loads for copper, lead, zinc, and four nutrient species (total nitrogen, ammonia nitrogen, total phosphorus, and orthophosphorus) in storm-water runoff at the drainage-well site were lower than mean storm load predicted from an existing regression model. The overprediction by the model may be a result of the small size of the drainage area relative to the range of drainage areas used in the development of the models, or to the below-normal amounts of rainfall during the period of sampling for this investigation. Loads& in storm-water runoff for 22 constituents were extrapolated from sampled storms to total loads for the period February to October 1988. Calculated loads for trace metals for the period ranged from 0.030pound.s for cadmium to 12pound.s for strontium. Loads of the primary nutrients ranged from 0.97pounds for nitrite as nitrogen to 34pounds of organic nitrogen. Storm-water quality at the drainage-well and Mobley Spring sites was compared to background water quality of the local aquifer; as characterized by dry-weather samples from three springs and two observation wells in the Clarksville area. Concentrations of total-recoverable cadmium, chromium, copper, lead, and nickel were higher in many stormwater samples from both the drainage-well and Mobley Spring sites than in samples from any other site. In addition, concentrations of total organic carbon, methylene blue active substances, and total-recoverable oil and grease were generally higher in storm-water samples from the drainage-well site than in any ground-water sample. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total recoverable iron, manganese, and methylene blue active substances in storm samples from the drainage-well site exceeded the maximum contaminant levels listed in Tennessee’s drinking-water standards (1988) by as much as 2,500 and 5,500 colonies per 100 milliliters, and 2.7, 0.29, and 0.05 milligrams per liter, respectively. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total-recoverable iron, manganese, and lead in storm samples from Mobley Spring exceeded the maximum contaminant levels by as much as 500 and 4,500 colonies per 100 milliliters, and 18.7,0.65, and 0.02 milligrams per liter, respectively. For iron, manganese, and bacteria, these undesirable levels are not necessarily attributable to storm-water recharge, because concentrations of these constituents also exceeded drinking-water standards in one or more of the dry-weather samples from selected springs and observation wells in the area.

Tennessee↗
USGS2009

Streamflow Simulations and Percolation Estimates Using the Soil and Water Assessment Tool for Selected Basins in North-Central Nebraska, 1940-2005

The U.S. Geological Survey, in cooperation with the Upper Elkhorn, Lower Elkhorn, Upper Loup, Lower Loup, Middle Niobrara, Lower Niobrara, Lewis and Clark, and Lower Platte North Natural Resources Districts, used the Soil and Water Assessment Tool to simulate streamflow and estimate percolation in north-central Nebraska to aid development of long-term strategies for management of hydrologically connected ground and surface water. Although groundwater models adequately simulate subsurface hydrologic processes, they often are not designed to simulate the hydrologically complex processes occurring at or near the land surface. The use of watershed models such as the Soil and Water Assessment Tool, which are designed specifically to simulate surface and near-subsurface processes, can provide helpful insight into the effects of surface-water hydrology on the groundwater system. The Soil and Water Assessment Tool was calibrated for five stream basins in the Elkhorn-Loup Groundwater Model study area in north-central Nebraska to obtain spatially variable estimates of percolation. Six watershed models were calibrated to recorded streamflow in each subbasin by modifying the adjustment parameters. The calibrated parameter sets were then used to simulate a validation period; the validation period was half of the total streamflow period of record with a minimum requirement of 10 years. If the statistical and water-balance results for the validation period were similar to those for the calibration period, a model was considered satisfactory. Statistical measures of each watershed model's performance were variable. These objective measures included the Nash-Sutcliffe measure of efficiency, the ratio of the root-mean-square error to the standard deviation of the measured data, and an estimate of bias. The model met performance criteria for the bias statistic, but failed to meet statistical adequacy criteria for the other two performance measures when evaluated at a monthly time step. A primary cause of the poor model validation results was the inability of the model to reproduce the sustained base flow and streamflow response to precipitation that was observed in the Sand Hills region. The watershed models also were evaluated based on how well they conformed to the annual mass balance (precipitation equals the sum of evapotranspiration, streamflow/runoff, and deep percolation). The model was able to adequately simulate annual values of evapotranspiration, runoff, and precipitation in comparison to reported values, which indicates the model may provide reasonable estimates of annual percolation. Mean annual percolation estimated by the model as basin averages varied within the study area from a maximum of 12.9 inches in the Loup River Basin to a minimum of 1.5 inches in the Shell Creek Basin. Percolation also varied within the studied basins; basin headwaters tended to have greater percolation rates than downstream areas. This variance in percolation rates was mainly was because of the predominance of sandy, highly permeable soils in the upstream areas of the modeled basins.

Scientific Investigations Report↗
USGS2005

Use of a ground-water flow model to delineate contributing areas to the Puchack Well Field, Pennsauken township and vicinity, Camden county, New Jersey

The New Jersey Department of Environmental Protection (NJDEP) Well Head Protection Program, developed in response to the 1986 Federal Safe Drinking Water Act Amendments, requires delineation of Well Head Protection Areas (WHPA's), commonly called contributing areas, for all public and non-community water-supply wells in New Jersey. Typically, WHPA's for public community water-supply wells in New Jersey are delineated using a two-dimensional ground-water flow model incorporating the regional hydraulic gradient; however, NJDEP guidelines allow for the use of a three-dimensional flow model to delineate contributing areas to wells in complex hydrogeologic settings. The Puchack well field in Pennsauken Township, Camden County, N.J., is an area of strong hydraulic connection between the Lower aquifer of the Potomac-Raritan-Magothy aquifer system and the Delaware River. Interactions among and within the public-supply well fields in the area are complex. To delineate the contributing area to the Puchack well field, the U.S. Geological Survey, in cooperation with the NJDEP, developed an 11-layer ground-water flow model of the Potomac-Raritan-Magothy aquifer system in the Pennsauken Township area to simulate flow in the vicinity of the well field. The model incorporates the interaction between the aquifer system and the Delaware River, and includes boundary flows from an existing regional model of the Camden area. Recharge used in the model ranged from 4.5 to 14 inches per year, and horizontal hydraulic conductivity ranged from 50 to 250 feet per day. Values of vertical hydraulic conductivity ranging from 0.001 to 0.5 feet per day were assigned to zones created on the basis of variations in hydrogeologic conditions observed in geophysical logs from wells. A steady-state simulation was used to calibrate the model to synoptic water-level data collected in March 1998. Near the Puchack well field, simulated heads generally were within 1 foot of the measured heads in both the Middle and Lower aquifers. Simulated water-level differences across the confining units at most of the nested wells were within ? 0.5 feet of the differences calculated from measured water levels. The existing flow model was modified to meet NJDEP guidelines for delineating contributing areas in complex hydrogeologic settings. These modifications included rediscretizing the model grid to a finer grid and preparing the water-use data set for use in the rediscretized model. The contributing area to the Puchack well field was delineated by means of particle tracking. An uncertainty analysis was conducted in which 36 model-input parameters were both increased and decreased until the resulting change in simulated heads exceeded the model-calibration criterion of ? 5 feet at any model cell. Porosity most affected the size and shape of the contributing area. The distribution of withdrawals at the Morris/Delair well field and variations in recharge affected both the size and shape of contributing area to the Puchack well field and the source of water to the Puchack wells. The results of the uncertainty analysis were combined to determine the 'aggregate' contributing area to the Puchack well field--a composite of areas on the land surface that contributed flow to the Puchack well field in less than 12 years in any uncertainty simulation. The shape of the aggregate contributing area was most similar to that associated with a reduction in porosity, which indirectly affected the size and shape of the contributing areas by changing travel time.

Scientific Investigations Report↗
USGS1980

Hydrology and model study of the proposed Prosperity Reservoir, Center Creek Basin, southwestern Missouri

A dam and reservoir have been proposed for construction on Center Creek, Jasper County, in southwestern Missouri. Ground-water levels in the hills adjacent to the reservoir will rise when the impoundment is completed. One of the problems is that the proposed site of Prosperity Reservoir is a few miles upstream from the lead-zinc mining area known as the Oronogo-Duenweg belt. In this belt transmissivities are variable but appear to be higher than they are in the immediate area of the reservoir. Grove Creek lies down-gradient from the reservoir area and separates it from the mining belt. A model study indicates that inflow from the proposed reservoir to the water table could cause water level rises varying from about 20 feet near the reservoir to 0.5 to 1.0 foot in the southern part of Grove Creek drainage basin. These rises will cause significant changes to the natural ground-water flow system. Increased ground-water elevations in the reservoir area could result in increased ground-water gradients and discharge to Grove and Center Creeks. The increase in ground-water discharge to Grove Creek, and in turn Center Creek, will have the beneficial effect of diluting mine-water discharge from the Oronogo-Duenweg belt during periods of low flow. However, if Grove Creek does not act as an effective drain and if conduits extend beneath Grove Creek to transfer the increased water available to the Oronogo-Duenweg belt, the flow regimen could change in the mining belt west of Grove Creek increasing mine-water discharge to Center Creek downstream from the reservoir. Bedrock in the area is Mississippian limestone, the deeply solutioned formation that contained the ore deposits. The limestone in the mining district was greatly altered by solution prior to ore deposition while the limestone in the area of the reservoir was altered less. The extent of the alteration is related to the aquifer characteristics in that high and low values of transmissivity and storage coefficient correspond to greatly altered brecciated rocks in the mining district and less altered, less brecciated rocks in the reservoir area, respectively. The authors suggest that an ancestral east-flowing White River drained the area about Joplin in Late Mississippian time. This is based on the configuration of the contact between Meramecian and Osagean rocks of Mississippian age. A high topographic area existed in the region about Joplin in which the water table stood 200 feet below the land surface when sinkholes and caverns of that depth were formed. The large number of Pennsylvanian-filled sinkholes in the Joplin area and the smaller number to the east suggest a higher land surface to the west than that to the east. The distribution of paleokarst sinkholes supports the conclusion based on the configuration of the Meramecian-Osagean contact.

Missouri↗
USGS2015

Geologic map of the Vashon 7.5' quadrangle and selected areas, King County, Washington

This map is an interpretation of a 6-ft-resolution lidar-derived digital elevation model combined with geology by Derek B. Booth and Kathy Goetz Troost. Field work by Booth and Troost was located on the 1:24,000-scale topographic map of the Vashon and Des Moines 7.5' quadrangles that were published in 1997 and 1995, respectively. Much of the geology was interpreted from landforms portrayed on the topographic maps, supplemented by field exposures, where available. In 2001, the Puget Sound Lidar Consortium (see http://pugetsoundlidar.org/) obtained a lidar-derived digital elevation model (DEM) for Vashon Island and the Des Moines quadrangle. For a brief description of lidar and this data acquisition program, see Haugerud and others (2003). This new DEM has a horizontal resolution of 6 ft (1.83 m) and mean vertical accuracy of about 1 ft (about 0.3 m). The greater resolution and accuracy of the lidar DEM facilitated a much-improved interpretation of many aspects of the surficial geology, especially the distribution and relative age of landforms and the materials inferred to comprise them. Booth and Troost were joined by Tabor to interpret the new lidar DEM but have done no futher field work for this map. This map, the Vashon quadrangle and selected adjacent areas, encompasses most of Vashon Island, Maury Island, and Three Tree Point in the south-central Puget Sound. One small area in the Vashon quadrangle on the east side of Puget Sound is excluded from this map but included on the adjacent Seattle quadrangle (Booth and others, 2005). The map displays a wide variety of surficial geologic deposits, which reflect many geologic environments and processes. Multiple ice-sheet glaciations and intervening nonglacial intervals have constructed a complexly layered sequence of deposits that underlie both islands to a depth of more than 300 m below sea level. These deposits not only record glacial and nonglacial history but also control the flow and availability of ground water, determine the susceptibility of the slopes to landslides, and provide economic reserves of sand and gravel. The islands are surrounded by channels of Puget Sound, some as deep as the islands are high (>600 ft (~200 m)). The shorelines provide many kilometers of well-exposed coastal outcrops that reveal abundant lithologic and stratigraphic details not ordinarily displayed in the heavily vegetated Puget Lowland.

Washington↗
USGS1987

Effects of proposed highway embankment modifications on water-surface elevations in the lower Pearl River flood plain near Slidell, Louisiana

Major flooding in the lower Pearl River basin in recent years has caused extensive damage to homes and highways in the area. In 1980 and 1983, Interstate Highway 10 and U.S. Highway 190 were overtopped. In 1983, the Interstate Highway 10 crossing was seriously damaged by the flood. The U.S. Geological Survey, in cooperation with the Louisiana Department of Transportation and Development, Office of Highways, used a two-dimensional finite-element surface-water flow model to evaluate the effects the proposed embankment modifications at Interstate Highway 10 and U.S. Highway 90 on the water-surface elevations in the lower Pearl River flood plain near Slidell, Louisiana. The proposed modifications that were considered for the 1983 flood are: (1) Removal of all highway embankments, the natural condition, (2) extension of the West Pearl River bridge by 1,000 feet at U.S. Highway 90, (3) construction of a new 250-foot bridge opening in the U.S. Highways 190 and 90, west of the intersection of the highways. The proposed highway bridge modifications also incorporated lowering of ground-surface elevations under the new bridges to sea level. The modification that provided the largest reduction in backwater, about 35 percent, was a new bridge in Interstate Highway 10. The modification of the West Pearl River bridge at U.S. Highway 90 and replacement of the bridge in U.S. Highway 190 provide about a 25% reduction in backwater each. For the other modification conditions that required structural modifications, maximum backwater computed on the west side of the flood plain ranges from 0.0 to 0.8 foot and on the east side from 0.0 to 0.6 foot. Results show that although backwater is greater on the west side of the flood plain than on the east side, upstream of highway embankments, backwater decreases more rapidly in the upstream direction on the west side of the flood plain than on the east side. Analysis of the proposed modifications indicates that backwater would still occur on the east and west sides of the flood plain, but values would be less than those computed with highway embankments in place. (Author 's abstract)

Water-Resources Investigations Report↗
USGS2003

Simulated effects of ground-water management scenarios on the Santa Fe group aquifer system, Middle Rio Grande Basin, New Mexico, 2001-40

Future conditions in the Santa Fe Group aquifer system through 2040 were simulated using the most recent revision of the U.S. Geological Survey groundwater- flow model for the Middle Rio Grande Basin. Three simulations were performed to investigate the likely effects of different scenarios of future groundwater pumping by the City of Albuquerque on the ground-water system. For simulation I, pumping was held constant at known year-2000 rates. For simulation II, pumping was increased to simulate the use of pumping to meet all projected city water demand through 2040. For simulation III, pumpingwas reduced in accordance with a plan by the City of Albuquerque to use surfacewater to meet most of the projectedwater demand. The simulations indicate that for each of the three pumping scenarios, substantial additional watertable declines would occur in some areas of the basin through 2040. However, the reduced pumping scenario of simulation III also results in water-table rise over a broad area of the city. All three scenarios indicate that the contributions of aquifer storage and river leakage to the ground-water system would change between 2000 and 2040. Comparisons among the results for simulations I, II, and III indicate that the various pumping scenarios have substantially different effects on water-level declines in the Albuquerque area and on the contribution of each water-budget component to the total budget for the ground-water system. Between 2000 and 2040, water-level declines for continued pumping at year-2000 rates are as much as 120 feet greater than for reduced pumping; water-level declines for increased pumping to meet all projected city demand are as much as 160 feet greater. Over the same time period, reduced pumping results in retention in aquifer storage of about 1,536,000 acre-feet of ground water as compared with continued pumping at year- 2000 rates and of about 2,257,000 acre-feet as compared with increased pumping. The quantity of water retained in the Rio Grande as a result of reduced pumping and the associated decrease in induced recharge from the river is about 731,000 acre-feet as compared with continued pumping at year-2000 rates and about 872,000 acre-feet as compared with increased pumping. Reduced pumping results in slight increases in the quantity of water lost from the groundwater system to evapotranspiration and agriculturaldrain flow compared with the other pumping scenarios.

Water-Resources Investigations Report↗