SEARCH · USGS Science

Results for “Ground Water Modeling Investigation Report”

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

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48Linked to original sources

USGS1996

Geohydrology and potential water-supply development on Bumkin, Gallops, Georges, Grape, Lovell, and Peddocks Islands, eastern Massachusetts

An investigation of the geohydrology and of the potential for water-supply development on several of the Boston Harbor Islands, eastern Massachusetts, was conducted to evaluate the possibility of developing a permanent small-capacity water supply to support recreational activities, such as camping, hiking, and swimming. The Boston Harbor Islands, including Bumkin, Gallops, Georges, Grape, Lovell, and Peddocks Islands are part of a larger group of glacially deposited drumlins, which are composed of thick, dense, homogeneous till in their core that are overlain by a thin layer of stratified-beach deposits. The surficial materials over-lie a weathered zone of the metasedimentary Cambridge Argillite in the Boston Harbor area and were deposited by continental ice sheets that covered New England twice during the late Pleistocene Epoch, and by near-shore processes in the Holocene Epoch. The thickness of these materials range from less than 1 to about 300 feet where present. The till was deposited by glacial ice and is characterized as an unsorted matrix of sand, silt, and clay with variable amounts of stones and large boulders. The stratified deposits primarily consist of sorted and layered sand and gravel that accumulated and formed the beaches and tombolos of the harbor islands. These deposits overlie the till at altitudes generally less than 10 feet above sea level. A cross-sectional, ground-water-flow model was developed to estimate depth to the water table for a hypothetical drumlin-island flow system, which was assumed to be representative of the drumlin islands in Boston Harbor. Areas were identified in each island flow system with the greatest potential for small-capacity water-supply development based on the model-calculated depth to water and surficial geology of the islands. Model-calculated depth to water estimates were used because of the lack of available hydrologic data for the islands. Model results indicate that the simulated depth to water is less than 20 feet within 240 feet from the shore of the hypothetical drumlin-island flow system. This area on the topographic maps of the six Boston Harbor Islands roughly coincides with the high transmissivity zones of stratified-beach deposits and weathered till on the lower slopes of the drumlins where ground-water discharge and surface and subsurface runoff occurs.

Massachusetts↗

USGS1994

Simulation of ground-water flow and potential land subsidence, upper Santa Cruz Basin, Arizona

A numerical ground-water flow model of the upper Santa Cruz basin in Pinal, Pima, and Santa Cruz Counties was developed to evaluate predevelopment conditions in 1940, ground-water withdrawals for 1940-86, and potential water-level declines and land subsidence for 1987-2024. Simulations of steady-state ground-water conditions indicate 12,900 acre-feet of ground-water inflow, 15,260 acre-feet of outflow, 53,000 acre-feet of pre- development pumpage, 29,840 acre-feet of mountain- front recharge, and 34,020 acre-feet of streamflow infiltration in 1940. Simulations of transient ground-water conditions indicate a total of 6.6 million acre-feet of net pumpage and 3.4 million acre-feet of water removed from aquifer storage for 1941-86. A difference of 1.2 million acre-feet between estimated and net pumpage is attributed to increased recharge from irrigation return flow, mine return flow, and infiltration of sewage effluent. Estimated natural recharge represents 40 percent of pumpage for 1966-86 and averaged 63,860 acre-feet per year for 1940-57 and 76,250 acre-feet per year for 1958-86. The increase in recharge after 1958 was coincident with above- average winter streamflow in the Santa Cruz River for 1959-86. Increased recharge after 1958 and decreased pumpage after 1975 contributed to decreased water-level declines or to recoveries after 1977 in wells near the Santa Cruz River and its tributaries. The results of projection simu- lations indicate that a maximum potential subsi- dence for 1987-2024 ranges from 1.2 feet for an inelastic specific storage of .0001 ft to 12 feet for an inelastic specific storage of .0015 ft. The simulations were made on the basis of pumpage and recharge rates from 1986 and by using a preconso- lidation-stress threshold of 100 feet. A permanent reduction in acquitard storage can range from 1 to 12 percent of the potential loss of 3.9 million acre-feet in aquifer-system storage for 1987-2024.

Water-Resources Investigations Report↗

USGS1981

Hydrologic effects of proposed changes in management practices, Winnebago Pool, Wisconsin

Various changes in the management practices for the Winnebago pool have been proposed. The objectives of the proposed changes are protecting wetlands adjacent to the lakes in the pool and assuring adequate flow in the Fox River downstream from Lake Winnebago. Three proposed operation plans for the Winnebago pool were studied to determine the effects on lake stage and discharge. A digital reservoir model simulated daily lake stage and discharge from the Winnebago pool. Inflow for this simulation was computed from daily discharge, recorded at Rapide Croche Dam, and daily stage, recorded at Oshkosh. This computed inflow is net inflow, combining the effects of streamflow, precipitation, evaporation, and net ground-water inflow. The model was adapted to simulate operation of the Winnebago pool according to each of the proposed plans. Because the important features of each operation plan were expressed in terms of stages to be attained at certain dates, simulated discharge was adjusted daily to keep simulated lake stage close to that specified in the plan. The simulated lake stage for each of the proposed plans exceeded the allowable maximum in some years. However, the simulated stage exceeded the allowable maximum less frequently and for shorter periods than did the observed stage. For the 61 years simulated (water years 1918-78), the simulated 7-day, 10- year low flow ((57,10) at Rapide Croche Dam was 1,000 ft^/s compared with 900 ft^/s for the observed record. The 7-day, 2-year low flow (Qy 2) was 1,500 ft^/s (simulated) versus 1,600 ft3/s (observed). The model predicted that executing the proposed operation plan would result in monthly Q7 JQ for August and September of 1,100 ft-Vs. This is no change for August and a 12-percent increase for September. In addition, the study illustrated that the discharge from the Winnebago pool could be managed to guarantee that the discharge in the Fox River downstream would not fall below a specified discharge up to 1,500 ft^/s. In most years this would have no effect on the stage of the pool. In a few years there would be a conflict between releasing water from the pool to maintain the desired discharge and maintaining water levels above the required minimum stage. The frequency of this possible conflict varies with the level of discharge desired in the Fox River from 1 year in 60 for a discharge of 1,100 ft^/s to about 1 year in 15 for a discharge of 1,500 ft^/s.

Wisconsin↗

USGS1996

Numerical model analysis of the effects of ground-water withdrawals on discharge to streams and springs in small basins typical of the Puget Sound Lowland, Washington

A numerical ground-water flow model of a hypothetical basin was constructed and used to investigate the effects of ground-water withdrawals on rates of natural discharge to streams and springs in small basins of the Puget Sound Lowland. Definitions of the topography, geology, drainage, and climate of the hypothetical basin were based on the features of typical small basins in the Puget Sound Lowland. This information was used to construct a 13-layer numerical ground-water flow model capable of simulating water levels, hydraulic gradients, and discharge to streams and springs. Three sequences of glacial drift and interglacial deposits were simulated in the model; each sequence consisted of recessional outwash, till, advance, outwash, and fine-grained interglacial sediments. Alluvial sediments of the major stream valleys and undifferentiated glacial and interglacial deposits were also included in the model. The model was calibrated by comparing simulated hydrologic conditions with expected conditions and making adjustments to values of hydraulic characteristics as needed. The model was calibrated to predevelop- ment conditions (those prior to pumping), and then used to simulate the effects of pumping on natural discharge to streams and springs. Seven series of simulations were made to investigate the effects of (1) distance from the well to a stream, (2) the presence of confining layers, (3) pumping rate, (4) depth of the pumped aquifer, (5) distance from the well to a bluff, (6) well density, and (7) recharge rate. The discharge of wells pumping from unconfined outwash aquifers on the drift plains is derived almost entirely from capture of natural discharge to nearby stream reaches. Increasing the lateral distance between the well and stream caused more of the well discharge to be captured from other streams on the drift plain. Pumping from aquifers separated from the stream by one or more confining layers caused a reduction in the effects of pumping on discharge to nearby streams that was offset by an increase in the effects on discharge to more distant streams and springs. The percentage of well discharge captured from springs on the bluff was sensitive to the distance of wells from the bluff. Simulations also showed that increased well density caused greater water-level decline locally, but, at equilibrium, did not affect the extent of the area affected by reduction of natural discharge to streams and springs. Finally, decreased recharge in areas where development had created impervious surfaces had a direct effect on the natural discharge rates to streams and springs. Increased recharge, however, increased natural discharge and offset the effects of well withdrawals. Further analysis of the time-dependent effects of with- drawals would provide additional insights, but would require the development of a transient version of the model.

Open-File Report↗

USGS2007

Geochemical effects of induced stream-water and artificial recharge on the Equus Beds Aquifer, South-Central Kansas, 1995-2004

Artificial recharge of the Equus Beds aquifer is part of a strategy implemented by the city of Wichita, Kansas, to preserve future water supply and address declining water levels in the aquifer of as much as 30 feet caused by withdrawals for water supply and irrigation since the 1940s. Water-level declines represent a diminished water supply and also may accelerate migration of saltwater from the Burrton oil field to the northwest and the Arkansas River to the southwest into the freshwater of the Equus Beds aquifer. Artificial recharge, as a part of the Equus Beds Ground-Water Recharge Project, involves capturing flows larger than base flow from the Little Arkansas River and recharging the water to the Equus Beds aquifer by means of infiltration or injection. The geochemical effects on the Equus Beds aquifer of induced stream-water and artificial recharge at the Halstead and Sedgwick sites were determined through collection and analysis of hydrologic and water-quality data and the application of statistical, mixing, flow and solute-transport, and geochemical model simulations. Chloride and atrazine concentrations in the Little Arkansas River and arsenic concentrations in ground water at the Halstead recharge site frequently exceeded regulatory criteria. During 30 percent of the time from 1999 through 2004, continuous estimated chloride concentrations in the Little Arkansas River at Highway 50 near Halstead exceeded the Secondary Drinking-Water Regulation of 250 milligrams per liter established by the U.S. Environmental Protection Agency. Chloride concentrations in shallow monitoring wells located adjacent to the stream exceeded the drinking-water criterion five times from 1995 through 2004. Atrazine concentrations in water sampled from the Little Arkansas River had large variability and were at or near the drinking-water Maximum Contaminant Level of 3.0 micrograms per liter as an annual average established by the U.S. Environmental Protection Agency. Atrazine concentrations were much smaller than the drinking-water criterion and were detected at much smaller concentrations in shallow monitoring wells and diversion well water located adjacent to the stream probably because of sorption on aquifer sediment. Before and after artificial recharge, large, naturally occurring arsenic concentrations in the recharge water for the Halstead diversion well and recharge site exceeded the Maximum Contaminant Level of 10 micrograms per liter established by the U.S. Environmental Protection Agency for drinking water. Arsenic and iron concentrations decreased when water was recharged through recharge basins or a trench; however, chemical precipitation and potential biofouling eventually may decrease the artificial recharge efficiency through basins and trenches. At the Sedgwick site, chloride concentrations infrequently exceeded regulatory criteria. Large concentrations of atrazine were treated to decrease concentrations to less than regulatory criteria. Recharge of treated stream water through recharge basins avoids potentially large concentrations of arsenic and iron that exist at the Halstead diversion site. Results from a simple mixing model using chloride as a tracer indicated that the water chemistry in shallow monitoring well located adjacent to the Little Arkansas River was 80 percent of stream water, demonstrating effective recharge of the alluvial aquifer by the stream. Results also indicated that about 25 percent of the water chemistry of the diversion well water was from the shallow part of the aquifer. Additionally, diverting water through a diversion well located adjacent to the stream removed about 75 percent of the atrazine, probably through sorption to aquifer sediment, and decreased the need for additional water treatment to remove atrazine.

Kansas↗

USGS1998

Estimated predevelopment discharge to streams from the High Plains Aquifer in northwestern Oklahoma, southwestern Kansas, and northwestern Texas

A study of the High Plains aquifer in Okla homa was initiated in 1996 to: (1) provide the information needed by the Oklahoma Water Resources Board to manage the quantity of water produced from the aquifer; and (2) provide base line water-chemistry data. The approach used to meet the first objective is to develop a digital ground-water flow model. The model will be cali brated, in part, by comparing simulated and esti mated predevelopment discharge from the aquifer to streams and cross-boundary flow. This report presents the estimated predevelopment discharge to streams from the High Plains aquifer. Streamflow data were the primary source of information used to estimate predevelopment dis charge from the High Plains aquifer. Data from 30 streamflow stations between the Arkansas and Canadian Rivers were considered in the analysis, and winter low-flow frequencies for 7-, 14-, and 30-day periods were determined for 25 stations. The 14-day low flow with a recurrence interval of 2 years was the primary value used to estimate pre development discharge from the aquifer. The streams that drain the eastern part of the High Plains aquifer in Kansas (generally east of 99.5 longitude) are estimated to have had large predevelopment discharge from the aquifer, and most of them received discharge from near their headwaters. For streams with more than one streamflow gage, the upper perennial reaches appeared to have gained more discharge from the aquifer than the lower reaches. The total predevel opment discharge from the aquifer in this area to several streams is estimated to have been about 312 cubic feet per second, not including discharge that probably went directly to the Arkansas River. The Cimarron River and its tributaries are estimated to have gained about 78 cubic feet per second, but nearly one-half that amount was lost in the lower reaches of the river. The cause of the loss in the lower reaches is unknown. The Beaver River and its tributaries are estimated to have gained a net of about 10 cubic feet per second above Wolf Creek with the upper reaches gaining more than the lower reaches. Wolf Creek is estimated to have gained 30 cubic feet per second over its total length.

Water-Resources Investigations Report↗

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↗

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↗

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↗