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 1,117 records · Page 62Linked to original sources
USGS1996

Geohydrology and conceptual model of a ground-water-flow system near a Superfund site in Cheshire, Connecticut

Degradation of ground-water quality has been identified in an area of the north-central part of the town of Cheshire, Connecticut. An investigation by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, was done during 1994-95 to characterize the unconsolidated glacial deposits and the sedimentary bedrock, integrate the local geohydrologic conditions with the regional geohydrologic system, and develop a conceptual understanding of ground-water flow in the study area. A regional ground-water-flow model developed for the region near the study area indicates that perennial streams, including Judd Brook and the Tenmile River, form hydrologic divides that separate the larger region into hydraulically independent flow systems. In the local study area, synoptic water-level measurements made in June 1995 indicate that ground water near the water table flows west and southwestward from the low hill on the eastern side of the area toward the pond and wetlands along Judd Brook. Water-level data indicate that there is good hydraulic connection between the unconsolidated materials and underlying fractured bedrock. Unconsolidated materials in the study area consist principally of glacial stratified deposits that are fine sand, silt, and clay of glaci- olacustrine origin; locally these overlie thin glacial till. The glacial sediments range in thickness from a few feet to about 25 ft in the eastern part of the study area and are as much as 100 ft thick in the western and southern part of the study area beneath the Judd Brook and Tenmile River valleys. Fluvial redbeds of the New Haven Arkose underlie the glacial deposits in the region; in the study area, the redbeds consist of (1) channel sandstone units, which are coarse sandstone to fine conglomerate, generally in 6- to 15-ft- thick sequences; and (2) overbank mudstone units, which are siltstone and silty sandstone with some fine sandstone, generally in 6- to 50-ft-thick sequences. Thin-bedded zones of siltstone that are particularly fissile are present locally within the mudstone units. Rock units strike northward and dip eastward at about 20. The eastward-dipping strata are cut by a consistent set of west to west-northwest dipping, high-angle fractures. These fractures are oriented perpendicular to bedding and are present mostly in the channel sandstone units, but locally extend into the mudstone units as well. Borehole-geophysical logging indicates that ground water flows along bedding planes in fissile zones and between fissile zones in high-angle fractures, which are perpendicular to bedding. The combined fracture types form an aquifer system in which ground water follows a stair-step flowpath, flowing horizontally through fissile zones and vertically through high-angle fractures. Heat-pulse flow meter measurements and borehole fluid-conductivity and temperature logs indicate that only a small subset of the fissile zones and some high-angle fractures are hydraulically significant. A generalized local-scale ground-water flow model based on a nonspecific, but realistic, rock and fracture geometry was developed for the study area. Simulations show that under nonpumping conditions at a hypothetical well located in the middle of the model, ground-water flow was separated into upper and lower zones in which flow paths differed but were generally from northeast to southwest. Several short-duration aquifer tests conducted in the study area indicate that there is good hydraulic connection in the fractures between the pumping well (CS-221) and two bedrock wells located approximately 100 ft to the north and south along bedding strike. During the short duration of the aquifer tests, there was no hydraulic connection in bedrock wells located to the east, perpendicular to the strike. A range of transmissivity of 27 to 46 ft2/d was calculated from the aquifer-test data for the fractured-bedrock aquifer at CS-221 and TH-2. Individual fracture zones identified by bo

Open-File Report↗
USGS2005

Simulated water sources and effects of pumping on surface and ground water, Sagamore and Monomoy flow lenses, Cape Cod, Massachusetts

The sandy sediments underlying Cape Cod, Massachusetts, compose an important aquifer that is the sole source of water for a region undergoing rapid development. Population increases and urbanization on Cape Cod lead to two primary environmental effects that relate directly to water supply: (1) adverse effects of land use on the quality of water in the aquifer and (2) increases in pumping that can adversely affect environmentally sensitive surface waters, such as ponds and streams. These considerations are particularly important on the Sagamore and Monomoy flow lenses, which underlie the largest and most populous areas on Cape Cod. Numerical models of the two flow lenses were developed to simulate ground-water-flow conditions in the aquifer and to (1) delineate areas at the water table contributing water to wells and (2) estimate the effects of pumping and natural changes in recharge on surface waters. About 350 million gallons per day (Mgal/d) of water recharges the aquifer at the water table in this area; most water (about 65 percent) discharges at the coast and most of the remaining water (about 28 percent) discharges into streams. A total of about 24.9 Mgal/d, or about 7 percent, of water in the aquifer is withdrawn for water supply; most pumped water is returned to the hydrologic system as return flow creating a state of near mass balance in the aquifer. Areas at the water table that contribute water directly to production wells total about 17 square miles; some water (about 10 percent) pumped from the wells flows through ponds prior to reaching the wells. Current (2003) steady-state pumping reduces simulated ground-water levels in some areas by more than 4 feet; projected (2020) pumping may reduce water levels by an additional 3 feet or more in these same areas. Current (2003) and future (2020) pumping reduces total streamflow by about 4 and 9 cubic feet per second (ft3/s), corresponding to about 5 percent and 9 percent, respectively, of total streamflow. Natural recharge varies with time, over both monthly and multiyear time scales. Monthly changes in recharge cause pond levels to vary between 1 and 2 feet in an average year; annual changes in recharge, which can be much larger than monthly variations, can cause pond levels to vary by more than 10 feet in some areas over a period of years. Streamflow, which also changes in response to changes in recharge, varies by a factor of two over an average year and can vary more over multiyear periods. On average, monthly pumping ranges from 15.8 Mgal/d in March to 45.3 Mgal/d in August. Pumping and the distribution of return flow can seasonally affect the hydrologic system by lowering ground-water and pond levels and by depleting streamflows, particularly in the summer months. Maximum drawdowns in March and August exceed 3 feet and 6 feet, respectively, for current (2003) pumping. Simulated drawdowns from projected (2020) pumping, relative to water levels representing 2003 pumping conditions, exceed 2 feet in March and 5 feet in August. Current (2003) and future (2020) pumping can decrease pond levels in some areas by more than 3 feet; drawdown generally is largest during the month of August of an average year. Over multiyear periods, seasonal pumping can lower pond levels in some areas by more than 4 feet; the effects of seasonal pumping are largest during periods of reduced recharge. Monthly streamflow depletion varies in individual streams but can exceed 2 ft3/s in some streams. The combined effects of seasonal pumping and drought can reduce pond levels by more than 10 feet below average levels. Water levels in Mary Dunn Pond, which is in an area of large current and projected pumping, are predicted (2020) to decline during drought conditions by about 10.6 feet: about 6.9 feet from lower recharge, about 2.3 feet from current (2003) pumping, and about 1.4 feet from additional future (2020) pumping. The results indicate that pumping generally does not cause substantial streamflow depletion and that the primary effect of pumping is on water levels in ponds. Natural changes in recharge account for most of the variation in pond levels; however, pumping can cause substantial declines in the levels of ponds near pumping wells. Also, the effects of pumping and recharge can combine to cause drawdowns of more than 10 feet in some areas.

Massachusetts↗
USGS1996

Hydraulic characteristics and nutrient transport and transformation beneath a rapid infiltration basin, Reedy Creek Improvement District, Orange County, Florida

The Reedy Creek Improvement District disposes of about 7.5 million gallons per day (1992) of reclaimed water through 85 1-acre rapid infiltration basins within a 1,000-acre area of sandy soils in Orange County, Florida. The U.S. Geological Survey conducted field experiments in 1992 at an individual basin to examine and better understand the hydraulic characteristics and nutrient transport and transformation of reclaimed water beneath a rapid infiltration basin. At the time, concentrations of total nitrogen and total phosphorus in reclaimed water were about 3 and 0.25 milligrams per liter, respectively. A two-dimensional, radial, unsaturated/saturated numerical flow model was applied to describe the flow system beneath a rapid infiltration basin under current and hypothetical basin loading scenarios and to estimate the hydraulic properties of the soil and sediment beneath a basin. The thicknesses of the unsaturated and saturated parts of the surficial aquifer system at the basin investigated were about 37 and 52 feet, respectively. The model successfully replicated the field-monitored infiltration rate (about 5.5 feet per day during the daily flooding periods of about 17 hours) and ground-water mounding response during basin operation. Horizontal and vertical hydraulic conductivity of the saturated part of the surficial aquifer system were estimated to be 150 and 45 feet per day, respectively. The field-saturated vertical hydraulic conductivity of the shallow soil, estimated to be about 5.1 feet per day, was considered to have been less than the full-saturation value because of the effects of air entrapment. Specific yield of the surficial aquifer was estimated to be 0.41. The upper 20 feet of the basin subsurface profile probably served as a system control on infiltration because of the relatively low field-saturated, vertical hydraulic conductivity of the sediments within this layer. The flow model indicates that, in the vicinity of the basin, flow in the deeper, saturated zone was relatively slow compared to the more vigorous flow in the shallow saturated zone. The large radial component of flow below the water table in the vicinity of the basin implies that reclaimed water moves preferentially in the shallow part of the saturated zone upon reaching the water table. Therefore, there may be some vertical stratification in the saturated zone, with recently infiltrated water overlying ambient water. The infiltration capacity at the basin would be unaffected by a small (less than 10 feet) increase in background water-table altitude, because the water table would remain below the system control on infiltration. However, water-table rises of 15 and 20 feet were estimated to reduce the infiltration capacity of the basin by 8 and 25 percent, respectively. Model simulations indicate that increasing ponded depth within the basin from 4 to 12 inches and from 4 to 24 inches would increase basin infiltration capacity by less than 6 and 11 percent, respectively. A loading strategy at the basin that relies on long, uninterrupted flooding was shown to offer the possibility of inducing a more anaerobic environment conducive to denitrification while maintaining reclaimed-water disposal capacity. Field measurements indicated that transient, elevated concentrations or "spikes" of nitrate (as high as 33 milligrams per liter as nitrogen) occurred at the leading edge of the infiltrating water and in the shallow saturated zone following a prolonged basin rest period. This phenomenon probably is the result of mineralization and nitrification of organic nitrogen retained within the subsurface during earlier basin loading events. The organic nitrogen was retained in the shallow soil (due to adsorption/straining) and the shallow saturated zone (due to deposition under slacking pore-water velocity). The magnitude of the nitrate spikes appears to be influenced by the scheduling of basin loading, with short flooding and resting periods being most favorable to minimization of nitrate spikes. Removal of nitrogen by denitrification from the percolating reclaimed water is minimal in the vicinity of the basin, probably because of the lack of reducing conditions and a relative paucity of organic carbon substrates. It is speculated that longer flooding periods could induce reducing conditions that would be favorable for nitrogen removal from the system, but probably would lead to more pronounced nitrate spiking. Phosphorus concentrations were decreased by about 90 percent from concentrations in reclaimed water after moving through the upper 15 feet of the soil profile. This most likely was a result of adsorption onto abundant iron and aluminum hydrous oxyhydroxide coatings on sand grains. However, some phosphorus (coarse fraction, organic) passes through the shallow soil and accumulates below the water table under slacking pore-water velocity. This phosphorus is immobilized by adsorption, fixation, or precipitation reactions during basin rest periods.

Florida↗
USGS2020

Storage capacity and sedimentation characteristics of the San Antonio Reservoir, California, 2018

The San Antonio Reservoir is a large water storage facility in Alameda County, California, and is a major component of the Hetch Hetchy Regional Water System (RWS). The RWS is a water-supply system owned and operated by the San Francisco Public Utilities Commission (SFPUC) and provides water for about 2.7 million people in the San Francisco, Santa Clara, Alameda, and San Mateo Counties. The San Antonio Reservoir is one of two RWS reservoirs in Alameda County and the third largest of the RWS reservoirs in the San Francisco Bay Area. The reservoir was formed by the James H. Turner Dam, which was completed in 1965. At the time of construction, the reservoir was estimated to have 50,500 acre-feet (acre-ft) of storage capacity. That early estimate was based on a 1963 pre-construction topographic map, which was drawn from aerial photographs. The capacity of the reservoir was later surveyed in 1994 and 2000. These two later surveys did not include the upper 18 feet (ft) of the reservoir, which represents roughly 30 percent of the overall storage volume. To determine the storage capacity and provide updated stage-capacity curves up to the spillway, the U.S. Geological Survey, in cooperation with the SFPUC, surveyed the bathymetry and shoreline of the reservoir in April 2018. The bathymetric survey was performed by making depth soundings using a boat-mounted, multibeam echosounder. At the time of the survey, the water level was between 13 and 14 ft below the spillway elevation. To measure capacity between the water line up to the spillway elevation, topography along most of the shoreline was surveyed from the boat using a terrestrial Light Detection and Ranging (LiDAR) scanner and in other areas by using ground-survey techniques. Location during bathymetric and topographic data collection was determined using a Global Navigation Satellite System-Real Time Network system. Vertical profiles of sound speed were collected periodically. The sound-speed profiles were used to spatially and temporally adjust the sound-speed calculations used to determine depth from the soundings. Approximately 125 kilometers (78 miles) of transects with a total of about 560 million depth soundings and topographic LiDAR points were collected (about 160 per square meter). In addition, approximately 500 topographic survey points were collected in shallow, wadable areas and on land near the upper reservoir area using a Global Navigation Satellite System receiver attached to a fixed length survey rod. Depth soundings, terrestrial LiDAR points, topographic survey points, and a digitized shoreline were merged and interpolated to generate a digital elevation model (DEM) of the reservoir. Gridded elevation data extracted from the DEM were then tabulated to determine total reservoir capacity and create reservoir stage-surface area and stage-storage capacity tables. Results of the reservoir capacity analysis indicated that the reservoir has 53,266 (plus or minus 140) acre-ft of storage capacity, which is an increase of 2,766 acre-ft (or 5.5 percent) greater than the original 1965 estimate; the increase is likely due to improved survey methods. Also, at the time of this 2018 survey, Intake #1 (the lowest intake) was not in operation. Intake #1 is estimated to be buried approximately 10 ft below the bed, whereas Intake #2 is about 20 ft above the bed. There are five intakes at different elevation levels; however, when consecutive lower intakes become inoperable due to sedimentation, the live storage capacity (capacity available for use) is reduced. At the time of this survey, the remaining live storage (above Intake #2) was approximately 52,363 acre-ft. The 2018 stage-capacity curve was compared to the original 1965 stage-capacity curve. Although overall, the changes indicate an increase in storage capacity, the change in volume at 372.7 ft North American Vertical Datum of 1988 (370 ft National Geodetic Vertical Datum of 1929, NGVD 29) shows a decrease of 733 acre-ft (the elevation of 370 ft NGVD 29 was used because it is the lowest elevation available for the 1965 stage-capacity curves). This finding agrees with the observed accumulation of sediment over Intake #1. That volume was converted to an annual sediment yield of 0.35 acre-ft per square mile (or 165 cubic meters per square kilometer), which is of the same order of magnitude as that found in other watersheds for the Coast Ranges in California. A decrease of 733 acre-ft between 1965 and 2018 thus represents a loss of 1.5 percent of the overall storage capacity in the reservoir. The updated stage-surface area and stage-capacity tables provided in this report and online ( https://doi.org/10.5066/P9KC9DU8 ) can be used by the SFPUC to improve reservoir operations and serve as an accurate baseline to monitor bathymetric changes in the future.

California↗
USGS2002

Ground-water, surface-water, and water-chemistry data, Black Mesa Area, northeastern Arizona — 2000–2001, and performance and sensitivity of the 1988 USGS numerical model of the N aquifer

The N aquifer is the major source of water in the 5,400-square-mile area of Black Mesa in northeastern Arizona. Availability of water is an important issue in this area because of continued industrial and municipal use, a growing population, and precipitation of about 6 to 14 inches per year. The monitoring program in Black Mesa has been operating since 1971 and is designed to determine the long-term effects of ground-water withdrawals from the N aquifer for industrial and municipal uses. The monitoring program includes measurements of (1) ground-water pumping, (2) ground-water levels, (3) spring discharge, (4) surface-water discharge, and (5) ground-water chemistry. In 2000, total ground-water withdrawals were 7,740 acre-feet, industrial use was 4,490 acre-feet, and municipal use was 3,250 acre-feet. From 1999 to 2000, total withdrawals increased by 9 percent, industrial use increased by 7 percent, and municipal use increased by 12 percent. From 1999 to 2001, water levels declined in 10 of 15 wells in the unconfined part of the aquifer, and the median change was -0.4 foot. Water levels declined in 8 of 16 wells in the confined part of the aquifer, and the median change was -0.2 foot. From the prestress period (prior to 1965) to 2001, the median water-level change for 33 wells was -17.2 feet. Median water-level changes were -1.2 feet for 15 wells in the unconfined part of the aquifer and -31.0 feet for 18 wells in the confined part. Discharges were measured once in 1999 and once in 2001 at four springs. Discharges decreased by 5 percent and 33 percent at two springs and increased by 3 percent and 81 percent at two springs. For about the past 10 years, discharges did not significantly change in Burro Spring, the unnamed spring near Dennehotso, and Moenkopi School Spring. The record of discharge from a consistent measuring point for Pasture Canyon Spring is too short for statistical analysis of trends. Continuous records of surface-water discharge have been collected from July 1976 to 2000 at Moenkopi Wash, July 1996 to 2000 at Laguna Creek, June 1993 to 2000 at Dinnebito Wash, and April 1994 to 2000 at Polacca Wash. Median flows for November, December, January, and February of each water year were used as an index of ground-water discharge to those streams. There is no significant trend in the median winter flows for Moenkopi Wash from 1977 to 2000. The records for the other three streams are too short for a statistical analysis of trends. The median winter flows for Dinnebito Wash and Polacca Wash, however, appear to have decreased during the last 6 years. There is no apparent trend in the median winter flows for Laguna Creek since 1997. In 2001, water samples were collected from 12 wells and 4 springs and analyzed for selected chemical constituents. Dissolved-solids concentrations ranged from 102 to 628 milligrams per liter. Water samples from 9 of the wells and from the 4 springs had less than 350 milligrams per liter of dissolved solids. Water-chemistry data with sufficient years of record for a statistical analysis of trends over time are available from 7 wells and 4 springs. From about the mid-1980s or early 1990s to 2001 there are no significant trends in the concentrations of dissolved solids, chloride, and sulfate in water samples from 6 of the 7 wells. The concentration of one tested constituent (dissolved solids) in samples from Rocky Ridge PM3 significantly increased from 1990 to 2001. From the late 1980s to 2001, there are no significant trends in the concentrations of dissolved solids, chloride, and sulfate in water samples from Burro Spring, the unnamed spring near Dennehotso, and Pasture Canyon Spring. From 1987 to 2001, concentrations of chloride and sulfate significantly increased in water samples from Moenkopi School Spring and concentrations of dissolved solids did not significantly change. The performance and sensitivity of the 1988 USGS numerical model of the N aquifer were analyzed. The overall performance of the model in steady-state conditions is reasonable for residuals of heads (difference between observed and simulated steady-state heads); 80 percent of the absolute values of residuals are less than 38 feet. Simulated flows are about 40 percent different than estimated flows at two of three discharge areas; however, this comparison is only a rough approximation of performance because the accuracy of the estimated steady-state flows is uncertain. The overall performance of the model for transient conditions is fair for residuals of changes in head (difference between observed and simulated changes in head from steady state to 1999); 80 percent of the absolute values of residuals are less than 31 feet. The model is biased in two areas. In the Tuba City area, simulated changes in head are more negative than observed changes in head; all six residuals are positive, and three residuals are between 75 and 155 feet. In the confined area of the aquifer, observed changes in head are more negative than simulated changes in head; 12 of the 17 residuals are negative, and 8 residuals are between -57 and -20 feet. Analysis of model sensitivity indicates that recharge, transmissivity, and storage coefficient are the most important parameters for estimating heads, changes in heads, and flows. A strong correlation between recharge and transmissivity and a lack of independent and reliable estimates of recharge, transmissivity, and discharge create a uniqueness problem in model calibration. Several models could be constructed and calibrated with different values of recharge or transmissivity and still have similar fits to the observed data. Information from recent data and studies and more advanced modeling techniques could be used to develop a more representative and less uncertain model. Future data collection and studies should focus on obtaining a better definition of recharge, discharge, transmissivity, and storage coefficient.

Arizona↗
USGS2001

Cross-hole radar scanning of two vertical, permeable, reactive-iron walls at the Massachusetts Military Reservation, Cape Cod, Massachusetts

A pilot-scale study was conducted by the U.S. Army National Guard (USANG) at the Massachusetts Military Reservation (MMR) on Cape Cod, Massachusetts, to assess the use of a hydraulic-fracturing method to create vertical, permeable walls of zero-valent iron to passively remediate ground water contaminated with chlorinated solvents. The study was conducted near the source area of the Chemical Spill-10 (CS-10) plume, a plume containing chlorinated solvents that underlies the MMR. Ground-water contamination near the source area extends from about 24 m (meters) to 35 m below land surface. The USANG designed two reactive-iron walls to be 12 m long and positioned 24 to 37 m below land surface to intersect and remediate part of the CS-10 plume.Because iron, as an electrical conductor, absorbs electromagnetic energy, the US Geological Survey used a cross-hole common-depth, radar scanning method to assess the continuity and to estimate the lateral and vertical extent of the two reactive-iron walls. The cross-hole radar surveys were conducted in boreholes on opposite sides of the iron injection zones using electric-dipole antennas with dominant center frequencies of 100 and 250 MHz. Significant decreases in the radar-pulse amplitudes observed in scans that traversed the injection zones were interpreted by comparing field data to results of two-dimensional finite-difference time-domain numerical models and laboratory-scale physical models.The numerical and physical models simulate a wall of perfectly conducting material embedded in saturated sand. Results from the numerical and physical models show that the amplitude of the radar pulse transmitted across the edge of a conductive wall is about 43 percent of the amplitude of a radar pulse transmitted across background material. The amplitude of a radar pulse transmitted through a hole in a conductive wall increases as the aperture of the hole increases. The modeling results indicate that holes with an aperture of less than 40 percent of the dominant wavelength of the radar pulse are not likely to be detected.Based on the results of the numerical and physical modeling, the decreases in radar-pulse amplitudes observed in scans traversing the injection zones are interpreted as electrically conductive zones that outline the distribution of iron. The area interpreted as iron in the northern A-wall contains two zones -- an upper zone about 10 m wide, extending from about 25 to 31 m below land surface, and a lower zone about 8 m wide, extending from 31.5 to 34.5 m below land surface. The area interpreted as iron in the southern B-wall is about 9 m wide, extending from about 27 to 34.5 m below land surface. No discrete holes were interpreted in either the A- or B-wall zones.The interpretation of the field data suggests that (1) the hydraulic-fracturing method introduced iron into the subsurface, but not in the dimensions originally proposed; (2) the iron within the treatment zones is distributed in a generally continuous manner; and (3) excluding the discontinuity in the A-wall, holes within the iron treatment zone, if any, exist at scales smaller than about 10 cm, the resolution limit of the radar antennas and acquisition geometry used for this study. The cross-hole radar method appears to have been an effective method for delineating the distribution of iron in the two walls; however, the veracity of the results cannot be ascertained without excavation or drilling into the treatment zone.

Massachusetts↗
USGS1984

Appraisal of the surficial aquifers in the Pomme de Terre and Chippewa River Valleys, western Minnesota

The surf icial sands in the Pomme de Terre and Chippewa River valleys in Grant, Pope, Stevens, and Swift Counties have been studied to determine the occurrence, availability, and quality of ground water in these aquifers. In the northern part of the Pomme de Terre and Chippewa River valleys, the aquifers consist of coarse sand and gravel ranging from 0 to 100 feet in thickness; transmissivities range from 0 to 35,000 feet squared per day in narrow, steep-sided erosional valleys. In the north, well yields commonly exceed 1,000 gallons per minute and may be as much as 4,000 gallons per minute locally. Farther south, the deposits are medium to fine grained, range from 0 to 90 feet thick, and reach a maximum width of 10 miles near Benson, Minnesota. Transmissivities range from 0 to 25,000 feet squared per day. Wells may yield as much as 1,500 gallons per minute locally. Southeast of Clontarf, well yields generally do not exceed 500 gallons per minute because the deposits are thinly saturated and fine grained. Ground water in the surf icial aquifer is a mixed calcium magnesium-sulf ate bicarbonate type that is chemically suitable for most uses. Concentrations of most constituents analyzed were below limits recommended by the Minnesota Pollution Control Agency for drinking water, but concentrations of manganese, iron, nitrite plus nitrate, and dissolved solids exceed recommended limits locally. Salinity, as indicated by the specific conductance (values ranged from 580 to 1,000 micromhos per centimeter) was in the medium to high range at several locations. An analytical model was used to estimate the effect on streamflow of pumpage from the surficial aquifer in the narrow, 50-mile reach of the Pomme de Terre River valley in Stevens and Grant Counties. The model indicates that the 43 existing wells pumping at maximum potential yields could reduce streamflow by 55 cubic feet per second. Addition of 23 wells also pumping at maximum potential yields could reduce streamflow by 77 cubic feet per second; this rate exceeds low base flow of the Pomme de Terre River. Finite-difference models were used to simulate flow in the surficial aquifer along the Pomme de Terre River near Appleton in Swift County and along the Chippewa River between Cyrus in Pope County and Danvers in Swift County. In the Appleton area/ model analyses indicate that pumping lowered water levels as much as 3 feet from 1973-80 and reduced streamflow by about 14 cubic feet per second. Additional regional water-level declines of 1 to 2 feet/ and up to 4 feet locally near aquifer-till boundaries, can be expected after about 4 years if pumping continues at the 1980 rate and area! recharge from precipitation is near normal. However/ simulation of increased pumping rates and decreased area! recharge during a 3-year drought indicates that water levels may decline as much as 9 feet near aquifer-till boundaries and streamflow may be reduced by about 41 cubic feet per second, which is about 95 percent of the available flow in the Pomme de Terre River at the 55-percent flow duration. Model results also suggest that/ during the first year of a drought/ the combined pumpage from wells operated during 1980 along the Pomme de Terre River in Stevens and Grant Counties and in the Appleton area could reduce streamflow to zero during base flow. Model-computed streamflow deficiencies are 48 and 60 cubic feet per second at the 55- and 70-percent flow duration/ respectively. Under such conditions/ pumping could not be sustained at the rates simulated unless there was sufficient water stored in the stream channel or streamflow was augmented. In the Cyrus-Benson area/ model results indicate that tinder 1980 development and average area! recharge/ dynamic equilibrium would be reached in less than 4 years and additional drawdown would be less than 2 feet. A 3-year drought coupled with increased pumping from irrigation wells operated during 1980 would lower water levels as much as 6 feet and reduce flow in the Chippewa River by about 26 cubic feet per second. At maximum hypothetical development in terms of the number of wells and normal area! recharge/ water levels would be lowered as much as 9 feet and streamflow would be reduced about 12 cubic feet per second. At maximum hypothetical development/ drought conditions and increased pumping would lower water levels as much as 12 feet and reduce flow in the Chippewa River by about 30 cubic feet per second/ which equals about 75 percent of available streamflow at the 70-percent flow duration.

Minnesota↗
USGS2001

Aquifer-characteristics data for West Virginia

Specific-capacity, storage-coefficient, and specific-yield data for wells in West Virginia were compiled to provide a data set from which transmissivity could be estimated. This data can be used for analytical and mathematical groundwater flow modeling. Analysis of available storage-coefficient and (or) specific-yield data indicates the Ohio River alluvial aquifer has a median specific yield of 0.20, which is characteristic of an unconfined aquifer. The Kanawha River alluvial aquifer has a median specific yield of 0.003, which is characteristic of a semi-confined aquifer. The median storage coefficient of fractured-bedrock aquifers is only 0.007, which is characteristic of confined aquifers. The highest median transmissivity of a specific aquifer in West Virginia occurs in Ohio River alluvium (4,800 ft2/d); the second highest occurs in Kanawha River alluvium (1,600 ft2/d). The lowest median transmissivity (23 ft2/d) is for the McKenzie-Rose Hill-Tuscarora aquifer. Rocks of Cambrian age within the Waynesboro-Tomstown-Harpers-Weverton-Loudon aquifer had a low median transmissivity of only 67 ft2/d. Other aquifers with low transmissivities include the Hampshire Formation, Brallier-Harrell Formations, Mahantango Formations, Oriskany Sandstone, and the Conococheague Formation with median transmissivities of 74, 72, 92, 82, and 92 ft2/d, respectively. All other aquifers within the State had intermediate values of transmissivity (130-920 ft2/d). The highest median transmissivities among bedrock aquifers were those for aquifers within the Pennsylvanian age Pocahontas Formation (1,200 ft2/d) and Pottsville Group (1,300 ft2/d), and the Mississippian age Mauch Chunk Group (1,300 ft2/d). These rocks crop out primarily in the southern part of the State and to a lesser extent within the Valley and Ridge Physiographic Province in West Virginia's Eastern Panhandle. The highest mean annual ground-water recharge rates within West Virginia (24.6 in.) occur within a band that extends through the central part of the State within the eastern part of the Kanawha River Basin. This area of relatively high relief has peaks higher than 4,000 ft and precipitation greater than 50 in./yr. The band of high recharge rates extends northward towards Pennsylvania and includes the Monongahela River Basin, which has a mean annual recharge of 21.4 inches. To the west of this central band lies a region of lower relief with much lower mean annual precipitation rates. Mean annual recharge for the Tug Fork, Twelvepole Creek, and Guyandotte River Basins is only 12.6 inches. For the western part of the Kanawha River Basin, mean recharge is 11.9 inches. The lowest mean annual recharge rates (8.4 in.) within the State occur in the Little Kanawha River Basin and the tributary streams in the region that discharge directly to the Ohio River. West Virginia's Eastern Panhandle is an area characterized by long linear northeast to southwest trending ridges and valleys. The mean annual ground-water recharge rate for this region, which is drained almost entirely by the Potomac River and its tributaries, is 9.4 inches. This area, which is located within a rain shadow resulting from orographic lifting in the higher altitude area to the west, receives less precipitation (approximately 30 in.) than the region to the west.

Water-Resources Investigations Report↗
USGS2005

Simulated water-level responses, ground-water fluxes, and storage changes for recharge scenarios along Rillito Creek, Tucson, Arizona

A local ground-water flow model is used to simulate four recharge scenarios along Rillito Creek in northern Tucson to evaluate mitigating effects on ground-water deficits and water-level declines in Tucson's Central Well Field. The local model, which derives boundary conditions from a basin-scale model, spans the 12-mile reach of Rillito Creek and extends 9 miles south into the Central Well Field. Recharge scenarios along Rillito Creek range from 5,000 to 60,000 acre-feet per year and are simulated to begin in 2005 and extend through 2225 to estimate long-term changes in ground-water level, ground-water storage, ground-water flux, and evapotranspiration. The base case for comparison of simulated water levels and flows, referred to as scenario A, uses a long-term recharge rate of 5,000 acre-feet per year to 2225. Scenario B, which increases the recharge along Rillito Creek by 9,500 acre-feet per year, has simulated water-level rises beneath Rillito Creek that range from about 53 feet to 86 feet. Water-level rises within the Central Well Field range from about 60 feet to 80 feet. More than half of these rises occur by 2050, and more than 95 percent occur by 2188. Scenario C, which increases the recharge along Rillito Creek by 16,700 acre-feet per year relative to scenario A, has simulated water-level rises beneath Rillito Creek that range from about 71 feet to 102 feet. Water-level rises within the Central Well Field range from about 80 feet to 95 feet. More than half of the rises occur by 2036, and more than 95 percent occur by 2100. Scenario D, which initially increases the recharge rate by about 55,000 acre-feet per year relative to scenario A, resulted in simulated water levels that rise to land surface along Rillito Creek. This rise in water level resulted in rejected recharge. As the water table continued to rise, the area of stream-channel surface intersected by the water table increased causing continual decline in the recharge rate until a long-term recharge rate of about 34,000 acre-feet per year was sustained. The long-term recharge rate for scenario D is about 29,000 acre-feet per year greater than the long-term recharge rate for scenario A. Simulated long-term water-level rises beneath Rillito Creek range from about 97 feet to 131 feet, resulting in water levels near or at the land surface. Shallow depths to water associated with this scenario have implications for contamination owing to the presence of landfills within or adjacent to Rillito Creek. Water-level rises for cells within the Central Well Field range from about 96 feet to 109 feet. More than half of the water-level rises occur by 2018 and more than 95 percent occur by 2041. Almost all the increased water added to the ground-water system in the recharge scenarios can be accounted for by a combination of increased storage near Rillito Creek, ground-water flux to the south, ground-water flux to the northwest, and increased discharge as evapotranspiration along Rillito Creek. The percentage of newly added water accounted for by storage changes is large relative to the percentage accounted for by changes in flux and evapotranspiration at the onset of each scenario; however, the changes in storage become smaller throughout the simulation, and the long-term component accounted for by storage is minimal. Long-term ground-water fluxes to the south increase by about 3,300, 4,840, and 7,500 acre-feet per year for scenarios B, C, and D, respectively. The percentage of increased recharge that flows south toward the Central Well Field, therefore, is 35, 29, and 26 percent for scenarios B, C, and D, respectively. Long-term ground-water fluxes to the northwest increase by about 3,100, 3,900, and 6,980 acre-feet per year for scenarios B, C, and D, respectively. The long-term percentage of increased recharge flowing northwestward is about 31, 25, and 21 percent for scenarios B, C, and D, respectively. Shallow ground-water evapotranspiration along Rillito Creek incr

Scientific Investigations Report↗
USGS2007

Assessment of Areal Recharge to the Spokane Valley-Rathdrum Prairie Aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

A numerical flow model of the Spokane Valley-Rathdrum Prairie aquifer currently (2007) being developed requires the input of values for areally-distributed recharge, a parameter that is often the most uncertain component of water budgets and ground-water flow models because it is virtually impossible to measure over large areas. Data from six active weather stations in and near the study area were used in four recharge-calculation techniques or approaches; the Langbein method, in which recharge is estimated on the basis of empirical data from other basins; a method developed by the U.S. Department of Agriculture (USDA), in which crop consumptive use and effective precipitation are first calculated and then subtracted from actual precipitation to yield an estimate of recharge; an approach developed as part of the Eastern Snake Plain Aquifer Model (ESPAM) Enhancement Project in which recharge is calculated on the basis of precipitation-recharge relations from other basins; and an approach in which reference evapotranspiration is calculated by the Food and Agriculture Organization (FAO) Penman-Monteith equation, crop consumptive use is determined (using a single or dual coefficient approach), and recharge is calculated. Annual recharge calculated by the Langbein method for the six weather stations was 4 percent of annual mean precipitation, yielding the lowest values of the methods discussed in this report, however, the Langbein method can be only applied to annual time periods. Mean monthly recharge calculated by the USDA method ranged from 53 to 73 percent of mean monthly precipitation. Mean annual recharge ranged from 64 to 69 percent of mean annual precipitation. Separate mean monthly recharge calculations were made with the ESPAM method using initial input parameters to represent thin-soil, thick-soil, and lava-rock conditions. The lava-rock parameters yielded the highest recharge values and the thick-soil parameters the lowest. For thin-soil parameters, calculated monthly recharge ranged from 10 to 29 percent of mean monthly precipitation and annual recharge ranged from 16 to 23 percent of mean annual precipitation. For thick-soil parameters, calculated monthly recharge ranged from 1 to 5 percent of mean monthly precipitation and mean annual recharge ranged from 2 to 4 percent of mean annual precipitation. For lava-rock parameters, calculated mean monthly recharge ranged from 37 to 57 percent of mean monthly precipitation and mean annual recharge ranged from 45 to 52 percent of mean annual precipitation. Single-coefficient (crop coefficient) FAO Penman-Monteith mean monthly recharge values were calculated for Spokane Weather Service Office (WSO) Airport, the only station for which the necessary meteorological data were available. Grass-referenced values of mean monthly recharge ranged from 0 to 81 percent of mean monthly precipitation and mean annual recharge was 21 percent of mean annual precipitation; alfalfa-referenced values of mean monthly recharge ranged from 0 to 85 percent of mean monthly precipitation and mean annual recharge was 24 percent of mean annual precipitation. Single-coefficient FAO Penman-Monteith calculations yielded a mean monthly recharge of zero during the eight warmest and driest months of the year (March-October). In order to refine the mean monthly recharge estimates, dual-coefficient (basal crop and soil evaporation coefficients) FAO Penman-Monteith dual-crop evapotranspiration and deep-percolation calculations were applied to daily values from the Spokane WSO Airport for January 1990 through December 2005. The resultant monthly totals display a temporal variability that is absent from the mean monthly values and demonstrate that the daily amount and timing of precipitation dramatically affect calculated recharge. The dual-coefficient FAO Penman-Monteith calculations were made for the remaining five stations using wind-speed values for Spokane WSO Airport and other assumptions regarding

Scientific Investigations Report↗
USGS2004

Simulation of solute transport of tetrachloroethylene in ground water of the glacial-drift aquifer at the Savage Municipal Well Superfund Site, Milford, New Hampshire, 1960-2000

The Savage Municipal Well Superfund site, named after the former municipal water-supply well for the town of Milford, is underlain by a 0.5-square mile plume of volatile organic compounds (VOCs), primarily tetrachloroethylene (PCE). The plume occurs mostly within a highly transmissive sand-and-gravel unit, but also extends to an underlying till and bedrock unit. The plume logistically is divided into two areas termed Operable Unit No. 1 (OU1), which contains the primary source area, and Operable Unit No. 2 (OU2), which is the extended plume area. PCE concentrations in excess of 100,000 parts per billion (ppb) had been detected in the OU1 area in 1995, indicating a likely Dense Non-Aqueous Phase Liquid (DNAPL) source. In the fall of 1998, the New Hampshire Department of Environmental Services (NHDES) and the U.S. Environmental Protection Agency (USEPA) installed a remedial system in OU1. The OU1 remedial system includes a low-permeability barrier that encircles the highest detected concentrations of PCE, and a series of injection and extraction wells. The barrier primarily sits atop bedrock and penetrates the full thickness of the sand and gravel; and in some places, the full thickness of the underlying basal till. The sand and gravel unit and the till comprise the aquifer termed the Milford-Souhegan glacial-drift aquifer (MSGD). Two-dimensional and three-dimensional finite-difference solute-transport models of the unconsolidated sediments (MSGD aquifer) were constructed to help evaluate solute-transport processes, assess the effectiveness of remedial activities in OU1, and to help design remedial strategies in OU2. The solute-transport models simulate PCE concentrations, and model results were compared to observed concentrations of PCE. Simulations were grouped into the following three time periods: an historical calibration of the distribution of PCE from the initial input (circa 1960) of PCE into the subsurface to the 1990s, a pre-remedial calibration from 1995 to 1998, and a remedial (post-barrier wall) calibration from 1998 to 1999. Model results also were checked against observed PCE concentrations from May and June 2000 as a post-audit of model performance. Results of the simulations of the two-dimensional model for the historical calibration indicate that the model-computed length of the plume is affected by the retardation factor (retardation). Values of retardation greater than 3 caused the longitudinal length of the computed plume to be too short compared to the observed plume. A retardation of 2-2.5 produced a reasonable comparison between computed and observed PCE concentrations. Testing of different starting times and rates of mass input of PCE indicated that the plume reaches a quasi steady-state distribution in about 20 years regardless of the rate of mass input or values of the solute-transport parameters (retardation, dispersion, and irreversible reaction) assigned the model. Results of the simulations of the three-dimensional model for the pre-remedial (1995-98) calibration of PCE for the OU2 area identified some spatial biases in computed concentrations that generally were unaffected by changes in retardation. The computed PCE concentrations exceeded observed concentrations along the northern part of the plume in OU2, where PCE increases were observed in a bedrock well. These results indicate that some PCE in this area may be entering the bedrock, which is not simulated in the model. Conversely, computed PCE concentrations were less than observed concentrations along the southern part of the plume in OU2. Because testing of high (above 4) values of retardation did little to reduce residuals, it is concluded that the low computed PCE concentrations along the southern flank are likely the result of an underestimation of the initial PCE mass in this area or an unaccounted source of PCE. Results of the simulations of the three-dimensional model for the remedial calibration period (1998-99) and po

Scientific Investigations Report↗
USGS2009

Occurrence and distribution of iron, manganese, and selected trace elements in ground water in the glacial aquifer system of the northern United States

Dissolved trace elements, including iron and manganese, are often an important factor in use of ground water for drinking-water supplies in the glacial aquifer system of the United States. The glacial aquifer system underlies most of New England, extends through the Midwest, and underlies portions of the Pacific Northwest and Alaska. Concentrations of dissolved trace elements in ground water can vary over several orders of magnitude across local well networks as well as across regions of the United States. Characterization of this variability is a step toward a regional screening-level assessment of potential human-health implications. Ground-water sampling, from 1991 through 2003, of the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey determined trace element concentrations in water from 847 wells in the glacial aquifer system. Dissolved iron and manganese concentrations were analyzed in those well samples and in water from an additional 743 NAWQA land-use and major-aquifer survey wells. The samples are from monitoring and water-supply wells. Concentrations of antimony, barium, beryllium, cadmium, chromium, cobalt, copper, iron, lead, manganese, molybdenum, nickel, selenium, strontium, thallium, uranium, and zinc vary as much within NAWQA study units (local scale; ranging in size from a few thousand to tens of thousands of square miles) as over the entire glacial aquifer system. Patterns of trace element concentrations in glacial aquifer system ground water were examined by using techniques suitable for a dataset with zero to 80 percent of analytical results reported as below detection. During the period of sampling, the analytical techniques changed, which generally improved the analytical sensitivity. Multiple reporting limits complicated the comparison of detections and concentrations. Regression on Order Statistics was used to model probability distributions and estimate the medians and other quantiles of the trace element concentrations. Strontium and barium were the most frequently detected and usually were present in the highest concentrations. Iron and manganese were the next most commonly detected and next highest in concentrations. Iron concentrations were the most variable with respect to the range of variations (both within local networks and aquifer-wide) and with respect to the disparity between magnitude of concentrations (detections) and the frequency of samples below reporting limits (nondetections). Antimony, beryllium, cadmium, silver, and thallium were detected too infrequently for substantial interpretation of their occurrence or distributions or potential human-health implications. For those elements that were more frequently detected, there are some geographic patterns in their occurrence that primarily reflect climate effects. The highest concentrations of several elements were found in the West-Central glacial framework area (High Plains and northern Plains areas). There are few important patterns for any element in relation to land use, well type, or network type. Shallow land-use (monitor) wells had iron concentrations generally lower than the glacial aquifer system wells overall and much lower than major-aquifer survey wells, which comprise mostly private- and public-supply wells. Unlike those for iron, concentration patterns for manganese were similar among shallow land-use wells and major-aquifer survey wells. An apparent relation between low pH and relatively low concentrations of many elements, except lead, may be more indicative of the relatively low dissolved-solids content in wells in the Northeastern United States that comprise the majority of low pH wells, than of a pH dependent pattern. Iron and manganese have higher concentrations and larger ranges of concentrations especially under more reducing conditions. Dissolved oxygen and well depth were related to iron and manganese concentrations. Redox conditions also affect several trace elements such as arsenic and copper; however, a comparison of redox categories, based in part on iron and manganese concentrations, indicated that the concentrations of many redox-sensitive elements were not significantly different among redox categories. Some of the redox-related patterns were not what would be expected on the basis of solubility constraints. Furthermore, barium is affected by redox conditions in at least one well network even though it is not a redox-sensitive element. Concentrations of barium in portions of the glacial aquifer system are limited by sulfate, which is strongly affected by redox conditions. Few samples had concentrations of any trace element that exceeded drinking-water standards (Maximum Contaminant Levels), for compounds regulated in drinking water or Health-Based Screening Levels for unregulated trace elements. More unregulated trace elements had concentrations greater than benchmarks than regulated trace elements. More samples had manganese concentrations greater its benchmark than any other element in the glacial aquifer system wells. Of the 1,590 wells sampled for manganese, only 556 are for private or public drinking-water supplies, and of those, 9.9 percent (55) exceeded the manganese Lifetime Health Advisory. Concentrations of arsenic, selenium, and uranium less frequently exceeded Maximum Contaminant Levels. There are 29 wells that had 2 element concentrations that exceeded their respective benchmarks. Most concentrations that exceeded a health-based benchmark were from wells in the West-Central area (Iowa, Minnesota, North and South Dakota, Nebraska, and Kansas); however, there is little geographical pattern to the wells with element concentrations of concern.

Scientific Investigations Report↗
USGS2011

Geohydrology, simulation of regional groundwater flow, and assessment of water-management strategies, Twentynine Palms area, California

The Marine Corps Air Ground Combat Center (MCAGCC) Twentynine Palms, California, overlies the Surprise Spring, Deadman, Mesquite, and Mainside subbasins of the Morongo groundwater basin in the southern Mojave Desert. Historically, the MCAGCC has relied on groundwater pumped from the Surprise Spring subbasin to provide all of its potable water supply. Groundwater pumpage in the Surprise Spring subbasin has caused groundwater levels in the subbasin to decline by as much as 190 feet (ft) from 1953 through 2007. Groundwater from the other subbasins contains relatively high concentrations of fluoride, arsenic, and (or) dissolved solids, making it unsuitable for potable uses without treatment. The potable groundwater supply in Surprise Spring subbasin is diminishing because of pumping-induced overdraft and because of more restrictive Federal drinking-water standards on arsenic concentrations. The U.S. Geological Survey, in cooperation with the MCAGCC, completed this study to better understand groundwater resources in the area and to help establish a long-term strategy for regional water-resource development. The Surprise Spring, Deadman, Mesquite, and Mainside subbasins are filled with sedimentary deposits of Tertiary age, alluvial fan deposits of Quaternary-Tertiary age, and younger alluvial and playa deposits of Quaternary age. Combined, this sedimentary sequence reaches a maximum thickness of more than 16,000 ft in the Deadman and Mesquite subbasins. The sedimentary deposits of Tertiary age yield a small amount of water to wells, and this water commonly contains high concentrations of fluoride, arsenic, and dissolved solids. The alluvial fan deposits form the principal water-bearing unit in the study area and have a combined thickness of 250 to more than 1,000 ft. The younger alluvial and playa deposits are unsaturated throughout most of the study area. Lithologic and downhole geophysical logs were used to divide the Quaternary/ Tertiary alluvial fan deposits into two aquifers (referred to as the upper and the middle aquifers) and the Tertiary sedimentary deposits into a single aquifer (referred to as the lower aquifer). In general, wells perforated in the upper aquifer yield more water than wells perforated in the middle and lower aquifers. The study area is dominated by extensive faulting and moderate to intense folding that has displaced or deformed the pre-Tertiary basement complex as well as the overlying Tertiary and Quaternary deposits. Many of these faults act as barriers to the lateral movement of groundwater flow and form many of the boundaries of the groundwater subbasins. The principal recharge to the study area is groundwater underflow across the western and southern boundaries that originates as runoff in the surrounding mountains. Groundwater discharges naturally from the study area as spring flow, as groundwater underflow to downstream basins, and as water vapor to the atmosphere by transpiration of phreatophytes and direct evaporation from moist soil. The annual volume of water that naturally recharged to or discharged from the groundwater flow system in the study area during predevelopment conditions was estimated to be 1,010 acre-feet per year (acre-ft/yr). About 90 percent of this recharge originated as runoff from the Little San Bernardino and the Pinto Mountains to the south, and the remainder originated as runoff from the San Bernardino Mountains to the west. Evapotranspiration by phreatophytes near Mesquite Lake (dry) was the primary form of predevelopment groundwater discharge. From 1953 through 2007, approximately 139,400 acre-feet (acre-ft) of groundwater was pumped by the MCAGCC from the Surprise Spring subbasin. A regional-scale numerical groundwater flow model was developed using MODFLOW-2000 for the Surprise Spring, Deadman, Mesquite, and Mainside subbasins. The aquifer system was simulated by using three model layers representing the upper, middle, and lower aquifers. Measured groundwater levels

California↗
USGS2014

Delineation of marsh types of the Texas coast from Corpus Christi Bay to the Sabine River in 2010

Coastal zone managers and researchers often require detailed information regarding emergent marsh vegetation types for modeling habitat capacities and needs of marsh-reliant wildlife (such as waterfowl and alligator). Detailed information on the extent and distribution of marsh vegetation zones throughout the Texas coast has been historically unavailable. In response, the U.S. Geological Survey, in cooperation and collaboration with the U.S. Fish and Wildlife Service via the Gulf Coast Joint Venture, Texas A&M University-Kingsville, the University of Louisiana-Lafayette, and Ducks Unlimited, Inc., has produced a classification of marsh vegetation types along the middle and upper Texas coast from Corpus Christi Bay to the Sabine River. This study incorporates approximately 1,000 ground reference locations collected via helicopter surveys in coastal marsh areas and about 2,000 supplemental locations from fresh marsh, water, and “other” (that is, nonmarsh) areas. About two-thirds of these data were used for training, and about one-third were used for assessing accuracy. Decision-tree analyses using Rulequest See5 were used to classify emergent marsh vegetation types by using these data, multitemporal satellite-based multispectral imagery from 2009 to 2011, a bare-earth digital elevation model (DEM) based on airborne light detection and ranging (lidar), alternative contemporary land cover classifications, and other spatially explicit variables believed to be important for delineating the extent and distribution of marsh vegetation communities. Image objects were generated from segmentation of high-resolution airborne imagery acquired in 2010 and were used to refine the classification. The classification is dated 2010 because the year is both the midpoint of the multitemporal satellite-based imagery (2009–11) classified and the date of the high-resolution airborne imagery that was used to develop image objects. Overall accuracy corrected for bias (accuracy estimate incorporates true marginal proportions) was 91 percent (95 percent confidence interval [CI]: 89.2–92.8), with a kappa statistic of 0.79 (95 percent CI: 0.77–0.81). The classification performed best for saline marsh (user’s accuracy 81.5 percent; producer’s accuracy corrected for bias 62.9 percent) but showed a lesser ability to discriminate intermediate marsh (user’s accuracy 47.7 percent; producer’s accuracy corrected for bias 49.5 percent). Because of confusion in intermediate and brackish marsh classes, an alternative classification containing only three marsh types was created in which intermediate and brackish marshes were combined into a single class. Image objects were reattributed by using this alternative three-marsh-type classification. Overall accuracy, corrected for bias, of this more general classification was 92.4 percent (95 percent CI: 90.7–94.2), and the kappa statistic was 0.83 (95 percent CI: 0.81–0.85). Mean user’s accuracy for marshes within the four-marsh-type and three-marsh-type classifications was 65.4 percent and 75.6 percent, respectively, whereas mean producer’s accuracy was 56.7 percent and 65.1 percent, respectively. This study provides a more objective and repeatable method for classifying marsh types of the middle and upper Texas coast at an extent and greater level of detail than previously available for the study area. The seamless classification produced through this work is now available to help State agencies (such as the Texas Parks and Wildlife Department) and landscape-scale conservation partnerships (such as the Gulf Coast Prairie Landscape Conservation Cooperative and the Gulf Coast Joint Venture) to develop and (or) refine conservation plans targeting priority natural resources. Moreover, these data may improve projections of landscape change and serve as a baseline for monitoring future changes resulting from chronic and episodic stressors.

Texas↗
USGS1999

Trends in base flows and extreme flows in the Beaver Kill Basin, Catskill Mountains, New York, 1915-94

Long-term records from five streamflow-gaging stations within and near the 300-square mile Beaver Kill Basin were analyzed to determine whether construction and presence of New York State Route 17 (NY 17), which was completed in the late 1960's, could have altered hydrologic processes in the basin and thereby adversely affected the basin's trout populations. The hypothesis investigated is that NY 17 has altered surface-water and shallow ground-water flowpaths where it parallels the stream and has increased runoff rates and thereby (1) increased the range in stream discharge (prolonged the base flows, decreased the low flows, and increased the high flows), and (2) altered stream-channel morphology through increased volume and velocity of stormflows. Analyses of base flows, discharge-duration curves, stage-to-discharge relations, peak and bankfull discharges, and flow extremes at a downstream (Beaver Kill at Cooks Falls) and a small tributary (Little Beaver Kill at Livingston Manor) site provide only limited evidence that NY 17 affected hydrologic processes within the basin. These effects are best indicated by significant increases in the magnitude and (or) the frequency of moderate to large discharges (exceedence probabilities) on an instantaneous basis at the Beaver Kill at Cooks Falls site after 1965. Increases in stormflows can not be attributed solely to NY 17, however, because the trend was evident long before NY 17 was constructed. Changes in land use in parts of the watershed may have contributed to gradual and continuous increases in stormflows throughout the entire 80 (plus) years of record. Changes in most base-flow and low-flow statistics for the downstream (Beaver Kill at Cooks Falls) site after 1965 are not statistically significant, but, changes in flow-duration curves and annual peak flows are evident. Flow-duration curves at this site indicate that there is a 16 percent increase in average daily flows after 1965. Annual peak flow data indicate that peak flows from storms recurring at 2-year (and longer) intervals after 1965 are significantly larger than those that recur at the same frequencies before 1965. The lack of comparable increases in peak flows from several nearby reference sites after 1965 indicate that the observed increases in peak flows may be unique to the Beaver Kill Basin. Flow-duration curves and many base-flow and high-flow statistics for the small tributary paralleled by NY 17 in the in the upper reaches of the basin (Little Beaver Kill at Livingston Manor) appear to be considerably altered since NY 17 was constructed. Flow-duration curves at this site indicate that there is about a 54 percent increase in average daily flows after 1965. Increases in the ratio of average annual base flow to average annual flow until 1965 then subsequent decreases suggest an extreme affect of NY 17 on hydrology of the subbasin. The effect of NY 17 on hydrology of the Little Beaver Kill subbasin cannot be defined with certainty, however, because the flow record after 1965 is too short; discharge monitoring was discontinued in 1981. The increases in peak stormflows in the lower Beaver Kill basin through the period of record may have increased the rates of bed-sediment erosion (degradation) and deposition and accelerated changes in stream-channel morphology, however, these possible effects were not examined. Suggestions for further investigation of the effects of NY 17 and of other factors on hydrology, channel morphology, fish habitat, and fish populations in the Beaver Kill Basin include (1) addition of streamflow gages or a creststage gage network at critical locations, (2) a review of engineering records and other aerial photographs for indications of changes in channel morphology, (3) compilation of temperature data and modeling spatial extent and magnitude of stressful summer temperatures (to selected trout species), and (4) confirming the extent and severity of toxic thermal episodes using in-situ fish toxicity tests.

Open-File Report↗
USGS2023

Integrated science strategy for assessing and monitoring water availability and migratory birds for terminal lakes across the Great Basin, United States

Executive Summary In 2022, the U.S. Geological Survey (USGS) established the Saline Lake Ecosystems Integrated Water Availability Assessment (IWAAs) to monitor and assess the hydrology of terminal lakes in the Great Basin and the migratory birds and other wildlife dependent on those habitats. Scientists from across the USGS (with specialties in water quantity, water quality, limnology, avian biology, data science, landscape ecology, and science communication) formed the Saline Lake Ecosystems IWAAs Team. The team has developed this regional strategic science plan to guide data collection and assessment activities at terminal lakes in the Great Basin. The U.S. Congress requested the USGS to establish the Saline Lake Ecosystems IWAAs in response to historically low water levels at terminal lakes and associated wetlands across the Great Basin. Not all Great Basin terminal lakes have high salinity; however, all terminal lakes occur in endorheic, closed, basins with no surface-water outflow. Low lake levels across the Great Basin are the result of increased water use for agriculture and municipalities, drought conditions, and a warming climate. Great Basin terminal lake water extents have decreased by as much as 90 percent over the last 150 years, and terminal lake wetlands have decreased in area by as much as 47 percent since 1984. Lake elevations and wetland areas are primarily supported by freshwater inputs from snowmelt feeding upgradient rivers, streams, and springs. These freshwater inputs have been severely reduced because of continued and increased surface-water diversions and surface-water capture through groundwater pumping for agriculture, mining, and public supply as well as unprecedented drought conditions and warming temperatures related to climate change. Water quality, specifically salinity, is highly variable for terminal lakes of the Great Basin, and this variability is a result of the balance between freshwater inflow and evaporation. Variability of salinity at each of the terminal lakes can be affected by lake morphology, hydrogeologic features of the basin, annual variability in weather patterns, and changes in upgradient water use. Hypersaline terminal lakes provide abundant food resources such as brine shrimp and brine flies that support nesting and migrating birds. The density and composition of invertebrates are closely tied to lake salinity. Increased salinity can exceed the tolerance of invertebrates, severely limiting their biomass. In contrast, decreased salinity can lead to altered invertebrate community composition, reducing the abundance of optimal avian prey resources. Great Basin terminal lake ecosystems, including open-water and adjacent aquatic and terrestrial environments, provide resources necessary to sustain many animal populations throughout the year. Although a variety of taxa use terminal lakes, these ecosystems are of acute importance for the millions of migratory waterbirds (for example, shorebirds, wading birds, and waterfowl) dependent on the network of terminal lakes and their associated wetlands. Migratory birds transiting the Pacific and Central Flyways use Great Basin terminal lake ecosystems throughout the year to feed, nest, and transit between wintering and breeding ranges. As such, successful conservation of birds and their habitats requires coordinated management of water and habitats across the Great Basin network of terminal lakes and wetlands. The linkages between water availability and ecosystem vulnerability of terminal lakes in the Great Basin are not well understood. The vulnerability of terminal lakes is related to the factors driving change and adaptive capacity of the lake ecosystem. Saline lake ecosystems are vulnerable when changes in water quantity affect ecosystem function. Water quantity affects salinity, which affects food webs and habitat; these linkages can be investigated with water-quality and food web monitoring. Water quantity also affects inundated habitat, which can be quantified through remote sensing. It is necessary to quantify hydroclimatic and water use controls on water availability to terminal lakes to assess the response of the ecosystems. Remotely sensed data can provide a broad-scale and long-term synoptic view of terminal lake hydrologic characteristics, but ground observations are required to interpret changes in water quality and ecological functions. Some terminal lake basins have ongoing monitoring and modeling efforts within the Great Basin (for example, Great Salt Lake, Carson River Basin), yet most monitoring locations are hydrologically upgradient and too far away from lake inflows to provide an accurate assessment of hydrological trends for the lake ecosystems. Other terminal lakes have no long-term hydrological monitoring in their respective watersheds (for example, Lake Abert). Ecological data collection in the Great Basin is also insufficient to understand how many birds exist on the landscape, how birds use the mosaic of terminal-lake habitats as an interconnected system, and how Great Basin terminal lakes are linked to the larger continental system of the Pacific and Central Flyways. Across agencies and organizations, tracking bird movement, abundance, and diversity is inconsistent, with some lakes having once- or twice-a-year bird survey efforts and a few locations having more intensive ecological data-gathering efforts (for example, Great Salt Lake, Lake Abert). Bridging hydrological and ecological information gaps will improve understanding of the trends in water supply and water quality, habitat availability and usage, and impacts on vulnerable waterbird species, all of which would be used by managers in coordinated conservation of this unique network of terminal-lake habitats. The terminal lakes of the Great Basin are part of the Basin and Range physiographic province that extends from the Colorado Plateau on the east to the Sierra Nevada on the west, and from the Snake River Plain on the north to the Garlock fault and the Mojave block on the south. The Great Basin is larger than 650,000 square kilometers and encompasses most of the State of Nevada but also extends to western Utah, eastern California, southeastern Idaho, southwestern Wyoming, and southeastern Oregon. The climate is arid to semiarid with a hydrologic regime that is snowmelt dominated, providing as much as 75 percent of total annual runoff for the region. Terminal lakes of the Great Basin occupy the lowest areas of closed (endorheic) drainage basins, such that lake levels and water quality respond rapidly to surface-water inflow. Terminal lakes provide local and regional economic value to the States in the Great Basin, including mineral extraction, aquaculture, public works, and recreational uses. As an example, assessments of Great Salt Lake’s ecological health and economic impact find hemispheric importance for the former and regional importance for the latter. Great Salt Lake creates about 7,000 jobs and $2 billion of economic output per year, most of which would be lost with further declines in lake level. The objectives of this Science Strategy are threefold: (1) to identify how changing water availability affects the quality, diversity, and abundance of habitats supporting continental waterbird populations; (2) to highlight the scientific monitoring and assessment needs of Great Basin terminal lakes; and (3) to support coordinated management and conservation actions to benefit those ecosystems, migratory birds, and other wildlife. There are long-term hydrological, ecological, and societal challenges associated with terminal lakes ecosystems in the Great Basin. This Science Strategy benefits partners by providing a conceptual model, nested at different spatial extents, that identifies key scientific information needs to inform coordinated implementation of management and conservation plans within and among hydrologic basins to address these complex challenges.

California, Idaho, Nevada, Oregon, Utah, Wyoming↗
USGS2010

Estimates of groundwater age from till and carbonate bedrock hydrogeologic units at Jefferson Proving Ground, Southeastern Indiana, 2007-08

During 2007-08, the U.S. Geological Survey, in cooperation with the U.S. Department of the Army, conducted a study to evaluate the relative age of groundwater in Pre-Wisconsinan till and underlying shallow and deep carbonate bedrock units in and near an area at Jefferson Proving Ground (JPG), southeastern Indiana, which was used during 1984-94 to test fire depleted uranium (DU) penetrators. The shallow carbonate unit includes about the upper 40 feet of bedrock below the bedrock-till surface; the deeper carbonate unit includes wells completed at greater depth. Samples collected during April 2008 from 15 wells were analyzed for field water-quality parameters, dissolved gases, tritium, and chlorofluorocarbon (CFC) compounds; samples from 14 additional wells were analyzed for tritium only. Water-level gradients in the Pre-Wisconsinan till and the shallow carbonate unit were from topographically higher areas toward Big Creek and Middle Fork Creek, and their tributaries. Vertical gradients were strongly downward from the shallow carbonate unit toward the deep carbonate unit at 3 of 4 paired wells where water levels recovered after development; indicating the general lack of flow between the two units. The lack of post development recovery of water levels at 4 other wells in the deep carbonate unit indicate that parts of that unit have no appreciable permeability. CFC and tritium-based age dates of Pre-Wisconsinan till groundwater are consistent with infiltration of younger (typically post-1960 age) recharge that 'mixes' with older recharge from less permeable or less interconnected strata. Part of the recharge to three till wells dated from the early to mid-1980s (JPG-DU-03O, JPG-DU-09O, and JPG-DU-10O). Age dates of young recharge in water from two till wells predated 1980 (JPG-DU-04O and JPG-DU-06O). Tritium-based age dates of water from seven other till wells indicated post-1972 age recharge. Most wells in the Pre-Wisconsinan till have the potential to produce groundwater that partially was recharged during or after DU penetrator testing; their water quality can indicate the presence of DU-related contaminants. The shallow carbonate unit near Big Creek is a karst flow system that may be recharged in part from areas with smaller thicknesses of overlying till or through more permeable parts of the till. This is indicated by CFC- and tritium-based piston-flow (non-mixing) model age dates of early-1980s for water from JPG-DU-02I, similar tritium-based ages of water produced from nearby wells MW-5 and MW-11, and cave development along the creek. The CFC and tritium-based age dates indicate that water samples from JPG-DU-01I and JPG-DU-03I were best described as mixtures of post-1984 modern recharge and submodern (1953 or older) recharge. These five wells produced groundwater that was recharged, at least partially, during or after DU-penetrator testing and are within or downgradient from the DU Impact Area with respect to groundwater flow directions inferred from water-level contours. Wells with groundwater age dates that are near to or after the onset (1984) of DU penetrator testing and that have a plausible connection to a contaminant source can be used to indicate the presence or absence of contaminants from DU penetrator or DU-related corrosion products in groundwater. Groundwater-age dates indicate that the ages of recharge sampled from shallow carbonate unit wells JPG-DU-04I, JPG-DU-05I, JPG-DU-06I, JPG-DU-09I, and JPG-DU-10D in easternmost (upgradient) and southernmost wells in the shallow carbonate unit are submodern (1953 or older) and predate the DU testing by at least 30 or more years. Water-quality data from these five wells are not likely to represent effects from DU-projectile testing or corrosion for years. Well JPG-DU-09D in the deep carbonate unit produced groundwater samples with a submodern (1953 or older) age date. The slow recovery of water levels in most wells in the deep carbonate unit is consis

Indiana↗
USGS1996

Hydrogeologic conditions and simulation of ground-water flow in the Greater Orlando Metropolitan Area, East-Central Florida

A finite-difference ground-water flow model was used to simulate the effects of both modern-day (1988) and projected 2010 ground-water withdrawals on the Floridan aquifer system in the greater Orlando metropolitan area. This area covers about 2,500 square miles and includes all of Orange and Seminole Counties and parts of Lake, Volusia, Brevard, Osceola, and Polk Counties. The hydrogeology of the area is characterized by a thin surficial aquifer underlain by the thick, highly productive rocks of the Floridan aquifer system. Water in the Upper Floridan aquifer is brackish (chloride concentrations greater than 1,000 milligrams per liter) in discharge areas beneath and near the St. Johns and Wekiva Rivers and is freshest (chloride concentrations less than 100 milligrams per liter) inrecharge areas. A slight trend toward increasing concentrations of dissolved solids, chloride, and sulfate has been observed at Upper Floridan aquifer springs. Chloride concentrations in the Upper Floridan aquifer measured between 1966 and 1993 at the Cocoa well field have increased from 50 milligrams per liter to 120 milligrams per liter; concentrations measured in the Lower Floridan aquifer between 1966 and 1993 have increasedfrom 600 milligrams per liter to 3,000 milligrams per liter. The flow model was calibrated by comparing (a) simulated and estimated Upper Floridan aquifer predevelopment (unstressed) potentiometric surfaces, (b) simulated and measured heads at 142 Upper Floridan aquifer monitoring wells in 1988 (averageabsolute error of 1.8 feet), (c) simulated and measured discharge rates at 15 Upper Floridan aquifer springs in 1988 (306 cubic feet per second), and (d) simulated and measured drawdowns at 134 Upper Floridan aquifer monitoring wells between 1988 and May 1990 (58 and 95 percent of simulated drawdowns were within plus or minus 25and 50 percent of measured drawdowns, respectively). Relative to predevelopment conditions, model simulations indicate that about half of the 305 million gallons per day of water pumped from the Floridan aquifer system in 1988 was accounted for by increased recharge from the surficial aquifer system. About 23 cubic feet persecond was derived from increased lateral inflow. A storage coefficient of 1x10-3 provided the best comparisons of measured-to-simulated data during the transient simulation from January to May 1990. This storativity probably is greater than the true storativity of the Upper Floridan aquifer because storage contributions from the intermediateconfining unit were not accounted for during model design and development. Calibrated transmissivity ranged from 10,000 to greater than 400,000 feet squared per day in the Upper Floridan aquifer, and from 5,000 to 600,000 feet squared per day in the Lower Floridan aquifer. Calibrated intermediate confining unit leakance ranged from 1x10-5 to 4x10-3 per day and was highest in areas where the unit is thin or has been breached by numerous sinkholes. In general,calibrated transmissivity and leakance values were higher than associated aquifer-test values. Simulated recharge rates to the Upper Floridan aquifer from the surficial aquifer system ranged from less than 3 to 21 inches per year. Recharge rates of greater than 10 inches per year were simulated in areas of west Seminole, west Orange, east Lake, and southwest Volusia Counties. Recharge rates of less than 3 inches per year were simulated in east Orange and northeast Osceola Counties. The calibrated model was used to simulate the effects of increased Floridan aquifer withdrawals in the year 2010 (542 million gallons per day) on water levels and spring flow. Projected effects were simulated for both "wet" conditions (using 1988 fixed-head arrays) and for "dry" conditions (using May 1990 fixed-head arrays), thus bracketing a potential range of effects. Relative to simulated 1988 conditions, simulated 2010 spring flow decreased by 43 cubic f

Water-Resources Investigations Report↗