USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology”

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

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

At least 973 records · Page 54Linked to original sources

Summary of Hydrologic Data for the Tuscarawas River Basin, Ohio, with an Annotated Bibliography

The Tuscarawas River Basin drains approximately 2,600 square miles in eastern Ohio and is home to 600,000 residents that rely on the water resources of the basin. This report summarizes the hydrologic conditions in the basin, describes over 400 publications related to the many factors that affect the groundwater and surface-water resources, and presents new water-quality information and a new water-level map designed to provide decisionmakers with information to assist in future data-collection efforts and land-use decisions. The Tuscarawas River is 130 miles long, and the drainage basin includes four major tributary basins and seven man-made reservoirs designed primarily for flood control. The basin lies within two physiographic provinces-the Glaciated Appalachian Plateaus to the north and the unglaciated Allegheny Plateaus to the south. Topography, soil types, surficial geology, and the overall hydrology of the basin were strongly affected by glaciation, which covered the northern one-third of the basin over 10,000 years ago. Within the glaciated region, unconsolidated glacial deposits, which are predominantly clay-rich till, overlie gently sloping Pennsylvanian-age sandstone, limestone, coal, and shale bedrock. Stream valleys throughout the basin are filled with sands and gravels derived from glacial outwash and alluvial processes. The southern two-thirds of the basin is characterized by similar bedrock units; however, till is absent and topographic relief is greater. The primary aquifers are sand- and gravel-filled valleys and sandstone bedrock. These sands and gravels are part of a complex system of aquifers that may exceed 400 feet in thickness and fill glacially incised valleys. Sand and gravel aquifers in this basin are capable of supporting sustained well yields exceeding 1,000 gallons per minute. Underlying sandstones within 300 feet of the surface also provide substantial quantities of water, with typical well yields of up to 100 gallons per minute. Although hydraulic connection between the sandstone bedrock and the sands and gravels in valleys is likely, it has not been assessed in the Tuscarawas River Basin. In 2001, the major land uses in the basin were approximately 40 percent forested, 39 percent agricultural, and 17 percent urban/residential. Between 1992 and 2001, forested land use decreased by 2 percent with correspondingly small increases in agricultural and urban land uses, but from 1980 to 2005, the 13-county area that encompasses the basin experienced a 7.1-percent increase in population. Higher population density and percentages of urban land use were typical of the northern, headwaters parts of the basin in and around the cities of Akron, Canton, and New Philadelphia; the southern area was rural. The basin receives approximately 38 inches of precipitation per year that exits the basin through evapotranspiration, streamflow, and groundwater withdrawals. Recharge to groundwater is estimated to range from 6 to 10 inches per year across the basin. In 2000, approximately 89 percent of the 116 million gallons per day of water used in the basin came from groundwater sources, whereas 11 percent came from surface-water sources. To examine directions of groundwater flow in the basin, a new dataset of water-level contours was developed by the Ohio Department of Natural Resources. The contours were compiled on a map that shows that groundwater flows from the uplands towards the valleys and that the water-level surface mimics surface topography; however, there are areas where data were too sparse to adequately map the water-level surface. Additionally, little is known about deep groundwater that may be flowing into the basin from outside the basin and groundwater interactions with surface-water bodies. Many previous reports as well as new data collected as part of this study show that water quality in the streams and aquifers in the Tuscarawas River Basin has been degraded by urban, suburban, and rural

Scientific Investigations Report↗

Baseline hydrologic studies in the lower Elwha River prior to dam removal

After the removal of two large, long‑standing dams on the Elwha River, Washington, the additional load of sediment and wood is expected to affect the hydrology of the lower river, its estuary, and the alluvial aquifer underlying the surrounding flood plain. To better understand the surface-water and groundwater characteristics of the river and estuary before dam removal, several hydrologic data sets were collected and analyzed. An experiment using a dye tracer characterized transient storage, and it was determined that the low‑flow channel of the lower Elwha River was relatively simple; 1–6 percent of the median travel time of dye was attributed to transient‑storage processes. Water data from monitoring wells adjacent to the main‑stem river indicated a strong hydraulic connectivity between stage in the river and groundwater levels in the flood plain. Analysis of temperature data from the monitoring wells showed that changes in the groundwater temperature responded weeks or months after water temperature changed in the river. A seepage investigation indicated that water from the river was moving into the aquifer (losing reach) between 1.7 and 2.8 kilometers from the river mouth. Surface‑water measurements and temperature and salinity data collected throughout the estuary helped to characterize the magnitude and nature of water movement in and out of the estuary. Salinity and stage sensors positioned in the estuarine network showed a strong surface‑water connection between the river and estuary waters east of the river. In contrast, there was a weaker connection between the river and estuarine water bodies west of the river.

Washington↗

Estimated hydrologic budgets of kettle-hole ponds in coastal aquifers of southeastern Massachusetts

Kettle-hole ponds in southeastern Massachusetts are in good hydraulic connection to an extensive coastal aquifer system that includes the Plymouth-Carver aquifer system on the mainland and aquifers underlying Cape Cod. The ponds receive water from, and contribute water to, the underlying glacial aquifer; ponds also receive water from precipitation and lose water to evaporation from the pond surface. Some ponds are connected to surface-water drainage systems and receive water from or contribute water to streams or adjacent wetlands. The Massachusetts Department of Environmental Protection currently (2011) is developing Total Maximum Daily Loads of phosphorus for the freshwater ponds in the region to maintain the health of pond ecosystems; the amounts and sources of water fluxes into and out of the ponds are important factors in determining the amount of phosphorus that can be assimilated into a pond. To assist in this effort, the U.S. Geological Survey used groundwater-flow models of the coastal aquifer system to estimate hydrologic budgets-including inflows and outflows from the aquifer system and adjacent streams and wetlands, and recharge from precipitation-for 425 ponds in southeastern Massachusetts. Water fluxes through the ponds are a function of several factors, including the size, shape, and bathymetry of the pond, orientation of the pond relative to the regional hydraulic gradient, and hydrologic setting relative to the proximity of groundwater divides and discharge boundaries. Total steady-state fluxes through the ponds range from more than 3,300,000 to less than 2,000 cubic feet per day. For ponds without surface-water inlets or outlets, groundwater inflow accounts for 98 to 3 percent of total inflow; conversely, recharge onto the pond surface accounts for the remainder of inflow (between 2 and 97 percent). All natural flows from these ponds are through recharge from the pond into the aquifer. In one pond, about 94 percent of the total outflow is removed for water supply. For ponds that are connected to surface-water drainages, most inflow and outflow are through streams. Ponds that receive water from streams receive most (58 to 89 percent) of their water from those streams. Ponds that are drained by streams lose between 5 and 100 percent of their water to those streams.

Scientific Investigations Report↗

Groundwater hydrology and estimation of horizontal groundwater flux from the Rio Grande at selected locations in Albuquerque, New Mexico, 2003-9

The Albuquerque, New Mexico, area has two principal sources of water: groundwater from the Santa Fe Group aquifer system and surface water from the San Juan-Chama Diversion Project. From 1960 to 2002, groundwater withdrawals from the Santa Fe Group aquifer system have caused water levels to decline more than 120 feet in some places within the Albuquerque area, resulting in a great deal of interest in quantifying the river-aquifer interaction associated with the Rio Grande. In 2003, the U.S. Geological Survey in cooperation with the Bureau of Reclamation, the Middle Rio Grande Endangered Species Collaborative Program, and the U.S. Army Corps of Engineers began a detailed characterization of the hydrogeology of the Rio Grande riparian corridor in the Albuquerque, New Mexico, area to provide hydrologic data and enhance the understanding of rates of water leakage from the Rio Grande to the alluvial aquifer, groundwater flow through the aquifer, and discharge of water from the aquifer to the riverside drains. A simple conceptual model of flow indicates that the groundwater table gently slopes from the Rio Grande towards riverside drains and the outer boundaries of the inner valley. Water infiltrating from the Rio Grande initially moves vertically below the river, but, as flow spreads farther into the Rio Grande inner valley alluvial aquifer, flow becomes primarily horizontal. The slope of the water-table surface may be strongly controlled by the riverside drains and influenced by other more distal hydrologic boundary conditions, such as groundwater withdrawals by wells. Results from 35 slug tests performed in the Rio Grande inner valley alluvial aquifer during January and February 2009 indicate that hydraulic-conductivity values ranged from 5 feet per day to 160 feet per day with a median hydraulic-conductivity for all transects of 40 feet per day. Median annual horizontal hydraulic gradients in the Rio Grande inner valley alluvial aquifer ranged from 0.011 to 0.002. Groundwater fluxes through the alluvial aquifer calculated by using median slug-test results (qm slug ) and Darcy's law ranged from about 0.1 feet per day to about 0.7 feet per day. Groundwater fluxes calculated by using the Suzuki-Stallman method (qm heat ) ranged from 0.52 feet per day to 0.23 feet per day. Results from the Darcy's law and Suzuki-Stallman flux calculations were compared to discharge measured in riverside drains on both sides of the river north of the Montaño Bridge on February 26, 2009. Flow in the Corrales Riverside Drain increased by 1.4 cubic feet per second from mile 2 to mile 4, about 12 cubic feet per day per linear foot of drain. Flow in the Albuquerque Riverside Drain increased by 15 cubic feet per second between drain miles 0 and 3, about 82 cubic feet per day per linear foot of drain. The flux of water from the river to the aquifer was calculated to be 2.2 cubic feet per day per linear foot of river by using the median qm slug of 0.09 feet per day at Montaño transects west of the river. The total flux was calculated to be 6.0 cubic feet per day per linear foot of river by using the mean(qm heat of 0.24 feet per day for the Montaño transects west of the river. Assuming the Corrales Riverside Drain intercepted all of this flow, the qm slug or qm heat fluxes account for 18 to 50 percent, respectively, of the increase of flow in the drain. The flux of water from the river to the aquifer was calculated to be 15 cubic feet per day per linear foot of river by using the median qm slug of 0.30 feet per day at the Montaño transects east of the river. The flux of water from the river to the aquifer was calculated to be 17 cubic feet per day per linear foot of river by using the mean flux calculated from the Suzuki-Stallman method for the Montaño East transects of 0.34 feet per day. Assuming the Albuquerque Riverside Drain intercepted all this flow, the qm slug or (qm heat fluxes would only account for 18 to 21 percent, respectively, of the increase in flow in the drain. The comparison of these results with those of previous investigations suggests that calculated flux through the Rio Grande inner valley alluvial aquifer is strongly scale dependent and that the thickness of aquifer through which river water flows may be greater than indicated by the vertical temperature profiles.

Albuquerque, New Mexico↗

Preliminary analysis of the hydrologic effects of temporary shutdowns of the Rondout-West Branch Water Tunnel on the groundwater-flow system in Wawarsing, New York

Flooding of streets and residential basements, and bacterial contamination of private-supply wells with Escherichia coli (E. coli) are recurring problems in the Rondout Valley near the Town of Wawarsing, Ulster County, New York. Leakage from the Rondout-West Branch (RWB) Water Tunnel and above-normal precipitation have been suspected of causing elevated groundwater levels and basement flooding. The hydrology of a 12-square-mile study area within the Town of Wawarsing was studied during 2008-10. A network of 41 wells (23 unconsolidated-aquifer and 18 bedrock wells) and 2 surface-water sites was used to monitor the hydrologic effects of four RWB Water Tunnel shutdowns. The study area is underlain by a sequence of northeast-trending sedimentary rocks that include limestone, shale, and sandstone. The bedrock contains dissolution features, fractures, and faults. Inflows that ranged from less than 1 to more than 9,000 gallons per minute from the fractured bedrock were documented during construction of the 13.5-foot-diameter RWB Water Tunnel through the sedimentary-rock sequence 710 feet (ft) beneath the study-area valley. Glacial sediments infill the valley above the bedrock sequence and consist of clay, silt, sand, and gravel. The groundwater-flow system in the valley consists of both fractured-rock and unconsolidated aquifers. Water levels in both the bedrock and unconsolidated aquifers respond to variations in seasonal precipitation. During the past 9 years (2002-10), annual precipitation at Central Park, N.Y., has exceeded the 141-year mean. Potentiometric-surface maps indicate that groundwater in the bedrock flows from the surrounding hills on the east and west sides of the valley toward the center of the valley, and ultimately toward the northeast. On average, water levels in the bedrock aquifer had seasonal differences of 5.3 ft. Analysis of hydrographs of bedrock wells indicates that many of these wells are affected by the RWB Tunnel leakage. Tunnel-leakage influences (water level and temperature changes) in the bedrock aquifer were measured at distances up to 7,000 ft from the RWB Tunnel. Water levels in the bedrock changed as much as 12 ft within 0.5 hour during tunnel shutdowns. Nine of the 10 wells that responded to the shutdowns showed a water-level response of 5 ft or greater. Changes in water levels ranged from 1.5 to 12 ft, with tunnel-leakage influence delay times ranging from 0.5 to 60 hours. Many of the longest tunnel-influence delay times and smallest water-level changes were in wells located closest to the tunnel in shale. Tunnel-influence response of the bedrock aquifer is consistent with its preliminary characterization as an anisotropic aquifer with greater transmissivity along bedding strike than across bedding strike. This tunnel-influence response is also consistent with the likely presence of discrete high-transmissivity networks along fractured limestone beds that have undergone dissolution. A lack of bedrock observation wells in half of the study area hampered a more thorough analysis of the extent of leakage from the RWB Tunnel in the study area. On average, water levels in the unconsolidated aquifer had a seasonal difference of 5.0 ft. Some unconsolidated-aquifer wells indicated water-level changes due to tunnel leakage. The locations of unconsolidated-aquifer wells with measurable water-level changes due to tunnel leakage correlated with those in the bedrock. Water levels in the unconsolidated aquifer changed as much as 2.5 ft within 18 hours of tunnel shutdowns, but water-level changes in some unconsolidated-aquifer wells were smaller or nonexistent.

New York↗

Water quality, hydrology, and simulated response to changes in phosphorus loading of Mercer Lake, Iron County, Wisconsin, with special emphasis on the effects of wastewater discharges

Mercer Lake is a relatively shallow drainage lake in north-central Wisconsin. The area near the lake has gone through many changes over the past century, including urbanization and industrial development. To try to improve the water quality of the lake, actions have been taken, such as removal of the lumber mill and diversion of all effluent from the sewage treatment plant away from the lake; however, it is uncertain how these actions have affected water quality. Mercer Lake area residents and authorities would like to continue to try to improve the water quality of the lake; however, they would like to place their efforts in the actions that will have the most beneficial effects. To provide a better understanding of the factors affecting the water quality of Mercer Lake, a detailed study of the lake and its watershed was conducted by the U.S. Geological Survey in collaboration with the Mercer Lake Association. The purposes of the study were to describe the water quality of the lake and the composition of its sediments; quantify the sources of water and phosphorus loading to the lake, including sources associated with wastewater discharges; and evaluate the effects of past and future changes in phosphorus inputs on the water quality of the lake using eutrophication models (models that simulate changes in phosphorus and algae concentrations and water clarity in the lake). Based on analyses of sediment cores and monitoring data collected from the lake, the water quality of Mercer Lake appears to have degraded as a result of the activities in its watershed over the past 100 years. The water quality appears to have improved, however, since a sewage treatment plant was constructed in 1965 and its effluent was routed away from the lake in 1995. Since 2000, when a more consistent monitoring program began, the water quality of the lake appears to have changed very little. During the two monitoring years (MY) 2008-09, the average summer near-surface concentration of total phosphorus was 0.023 mg/L, indicating the lake is borderline mesotrophic-eutrophic, or has moderate to high concentrations of phosphorus, whereas the average summer chlorophyll a concentration was 3.3 mg/L and water clarity, as measured with a Secchi depth, was 10.4 ft, both indicating mesotrophic conditions or that the lake has a moderate amount of algae and water clarity. Although actions have been taken to eliminate the wastewater discharges, the bottom sediment still has slightly elevated concentrations of several pollutants from wastewater discharges, lumber operations, and roadway drainage, and a few naturally occurring metals (such as iron). None of the concentrations, however, were high enough above the defined thresholds to be of concern. Based on nitrogen to phosphorus ratios, the productivity (algal growth) in Mercer Lake should typically be limited by phosphorus; therefore, understanding the phosphorus input to the lake is important when management efforts to improve or prevent degradation of the lake water quality are considered. Total inputs of phosphorus to Mercer Lake were directly estimated for MY 2008-09 at about 340 lb/yr and for a recent year with more typical hydrology at about 475 lb/yr. During these years, the largest sources of phosphorus were from Little Turtle Inlet, which contributed about 45 percent, and the drainage area near the lake containing the adjacent urban and residential developments, which contributed about 24 percent. Prior to 1965, when there was no sewage treatment plant and septic systems and other untreated systems contributed nutrients to the watershed, phosphorus loadings were estimated to be about 71 percent higher than during around 2009. In 1965, a sewage treatment plant was built, but its effluent was released in the downstream end of the lake. Depending on various assumptions on how much effluent was retained in the lake, phosphorus inputs from wastewater may have ranged from 0 to 342 lb. Future highway and stormwater improvements have been identified in the Mercer Infrastructure Improvement Project, and if they are done with the proposed best management practices, then phosphorus inputs to the lake may decrease by about 40 lb. Eutrophication models [Canfield and Bachman model (1981) and Carlson Trophic State Index equations (1977)] were used to predict how the water quality of Mercer Lake should respond to changes in phosphorus loading. A relatively linear response was found between phosphorus loading and phosphorus and chlorophyll a concentrations in the lake, with changes in phosphorus concentrations being slightly less (about 80 percent) and changes in chlorophyll a concentrations being slightly more (about 120 percent) than the changes in phosphorus loadings to the lake. Water clarity, indicated by Secchi depths, responded more to decreases in phosphorus loading than to increases in loading. Results from the eutrophication models indicated that the lake should have been negatively affected by the wastewater discharges. Prior to 1965, when there was no sewage treatment plant effluent and inputs from the septic systems and other untreated systems were thought to be high, the lake should have been eutrophic; near the surface, average phosphorus concentrations were almost 0.035 mg/L, chlorophyll a concentrations were about 7 μg/L, and Secchi depths were about 6 ft, which agreed with the shallower Secchi depths during this time estimated from the sediment-core analysis. The models indicated that between 1965 and 1995, when the lake retained some of the effluent from the new sewage treatment plant, water quality should have been between the conditions estimated prior to 1965 and what was expected during typical hydrologic conditions around MY 2008-09. The models also indicated that if the future Mercer Infrastructure Improvement Project is conducted with the best management practices as proposed, the water quality in the lake could improve slightly from that measured during 2006-10. Because of the small amount of phosphorus that is presently input into Mercer Lake any additional phosphorus added to the lake could degrade water quality; therefore, management actions can usefully focus on minimizing future phosphorus inputs. Phosphorus released from the sediments of a degraded lake often delays its response to decreases in external phosphorus loading, especially in shallow, frequently mixed systems. Mercer Lake, however, remains stratified throughout most of the summer, and phosphorus released from the sediments represents only about 6 percent, or a small fraction, of the total phosphorus load to the lake. Therefore, the phosphorus trapped in the sediments should minimally affect the long-term water quality of the lake and should not delay the response in its productivity to future changes in nutrient loading from its watershed.

Wisconsin↗

Hydrogeology, groundwater seepage, nitrate distribution, and flux at the Raleigh hydrologic research station, Wake County, North Carolina, 2005-2007

rom 2005 to 2007, the U.S. Geological Survey and the North Carolina Department of Environment and Natural Resources, Division of Water Quality, conducted a study to describe the geologic framework, measure groundwater quality, characterize the groundwater-flow system, and describe the groundwater/surface-water interaction at the 60-acre Raleigh hydrogeologic research station (RHRS) located at the Neuse River Waste Water Treatment Plant in eastern Wake County, North Carolina. Previous studies have shown that the local groundwater quality of the surficial and bedrock aquifers at the RHRS had been affected by high levels of nutrients. Geologic, hydrologic, and water-quality data were collected from 3 coreholes, 12 wells, and 4 piezometers at 3 well clusters, as well as from 2 surface-water sites, 2 multiport piezometers, and 80 discrete locations in the streambed of the Neuse River. Data collected were used to evaluate the three primary zones of the Piedmont aquifer (regolith, transition zone, and fractured bedrock) and characterize the interaction of groundwater and surface water as a mechanism of nutrient transport to the Neuse River. A conceptual hydrogeologic cross section across the RHRS was constructed using new and existing data. Two previously unmapped north striking, nearly vertical diabase dikes intrude the granite beneath the site. Groundwater within the diabase dike appeared to be hydraulically isolated from the surrounding granite bedrock and regolith. A correlation exists between foliation and fracture orientation, with most fractures striking parallel to foliation. Flowmeter logging in two of the bedrock wells indicated that not all of the water-bearing fractures labeled as water bearing were hydraulically active, even when stressed by pumping. Groundwater levels measured in wells at the RHRS displayed climatic and seasonal trends, with elevated groundwater levels occurring during the late spring and declining to a low in the late fall. Vertical gradients in the groundwater discharge area near the Neuse River were complex and were affected by fluctuations in river stage, with the exception of a well completed in a diabase dike. Water-quality data from the wells and surface-water sites at the RHRS were collected continuously as well as during periodic sampling events. Surface-water samples collected from a tributary were most similar in chemical composition to groundwater found in the regolith and transition zone. Nitrate (measured as nitrite plus nitrate, as nitrogen) concentrations in the sampled wells and tributary ranged from about 5 to more than 120 milligrams per liter as nitrogen. Waterborne continuous resistivity profiling conducted on the Neuse River in the area of the RHRS measured areas of low apparent resistivity that likely represent groundwater contaminated by high concentrations of nitrate. These areas were located on either side of a diabase dike and at the outfall of two unnamed tributaries. The diabase dike preferentially directed the discharge of groundwater to the Neuse River and may isolate groundwater movement laterally. Discrete temperature measurements made within the pore water beneath the Neuse River revealed seeps of colder groundwater discharging into warmer surface water near a diabase dike. Water-quality samples collected from the pore water beneath the Neuse River indicated that nitrate was present at concentrations as high as 80 milligrams per liter as nitrogen on the RHRS side of the river. The highest concentrations of nitrate were located within pore water collected from an area near a diabase dike that was identified as a suspected seepage area. Hydraulic head was measured and pore water samples were collected from two 140-centimeter-deep (55.1-inch-deep) multiport piezometers that were installed in bed sediments on opposite sides of a diabase dike. The concentration of nitrate in pore water at a suspected seepage area ranged from 42 to 82 milligrams per liter as nitrogen with a median concentration of 79 milligrams per liter as nitrogen. On the opposite side of the dike, concentrations of nitrate in pore water samples ranged from 3 to 91 milligrams per liter as nitrogen with a median concentration of 52 milligrams per liter. At one of the multiport piezometers the vertical gradient of hydraulic head between the Neuse River and the groundwater was too small to measure. At the multiport piezometer located in the suspected seepage area, an upward gradient of about 0.1 was present and explains the occurrence of higher concentrations of nitrate near the sediment/water interface. Horizontal seepage flux from the surficial aquifer to the edge of the Neuse River was estimated for 2006. Along a 130-foot flow path, the estimated seepage flux ranged from –0.52 to 0.2 foot per day with a median of 0.09 foot per day. The estimated advective horizontal mass flux of nitrate along a 300-foot reach of the Neuse River ranged from –10.9 to 5 pounds per day with a median of 2.2 pounds per day. The total horizontal mass flux of nitrate from the surficial aquifer to the Neuse River along the 130-foot flow path was estimated to be about 750 pounds for all of 2006. Seepage meters were deployed on the bed of the Neuse River in the areas of the multiport piezometers on either side of the diabase dike to estimate rates of vertical groundwater discharge and flux of nitrate. The average estimated daily seepage flux differed by two orders of magnitude between seepage areas. The potential vertical flux of nitrate from groundwater to the Neuse River was estimated at an average of 2.5 grams per day near one of the multiport piezometers and an average of 784 grams per day at the other. These approximations suggest that under some hydrologic conditions there is the potential for substantial quantities of nitrate to discharge from the groundwater to the Neuse River.

North Carolina↗

Geology, water-quality, hydrology, and geomechanics of the Cuyama Valley groundwater basin, California, 2008--12

To assess the water resources of the Cuyama Valley groundwater basin in Santa Barbara County, California, a series of cooperative studies were undertaken by the U.S. Geological Survey and the Santa Barbara County Water Agency. Between 2008 and 2012, geologic, water-quality, hydrologic and geomechanical data were collected from selected sites throughout the Cuyama Valley groundwater basin. Geologic data were collected from three multiple-well groundwater monitoring sites and included lithologic descriptions of the drill cuttings, borehole geophysical logs, temperature logs, as well as bulk density and sonic velocity measurements of whole-core samples. Generalized lithologic characterization from the monitoring sites indicated the water-bearing units in the subsurface consist of unconsolidated to partly consolidated sand, gravel, silt, clay, and occasional cobbles within alluvial fan and stream deposits. Analysis of geophysical logs indicated alternating layers of finer- and coarser-grained material that range from less than 1 foot to more than 20 feet thick. On the basis of the geologic data collected, the principal water-bearing units beneath the monitoring-well sites were found to be composed of younger alluvium of Holocene age, older alluvium of Pleistocene age, and the Tertiary-Quaternary Morales Formation. At all three sites, the contact between the recent fill and younger alluvium is approximately 20 feet below land surface. Water-quality samples were collected from 12 monitoring wells, 27 domestic and supply wells, 2 springs, and 4 surface-water sites and were analyzed for a variety of constituents that differed by site, but, in general, included trace elements; nutrients; dissolved organic carbon; major and minor ions; silica; total dissolved solids; alkalinity; total arsenic and iron; arsenic, chromium, and iron species; and isotopic tracers, including the stable isotopes of hydrogen and oxygen, activities of tritium, and carbon-14 abundance. Of the 39 wells sampled, concentrations of total dissolved solids and sulfate from 38 and 37 well samples, respectively, were greater than the U.S. Environmental Protection Agency’s secondary maximum contaminant levels. Concentrations greater than the maximum contaminant levels for nitrate were observed in five wells and were observed for arsenic in four wells. Differences in the stable-isotopic values of hydrogen and oxygen among groundwater samples indicated that water does not move freely between different formations or between different zones within the Cuyama Valley. Variations in isotopic composition indicated that recharge is derived from several different sources. The age of the groundwater, expressed as time since recharge, was between 600 and 38,000 years before present. Detectable concentrations of tritium indicated that younger water, recharged since the early 1950s, is present in parts of the groundwater basin. Hydrologic data were collected from 12 monitoring wells, 56 domestic and supply wells, 3 surface-water sites, and 4 rainfall-gaging stations. Rainfall in the valley averaged about 8 inches annually, whereas the mountains to the south received between 12 and 19 inches. Stream discharge records showed seasonal variability in surface-water flows ranging from no-flow to over 1,500 cubic feet per second. During periods when inflow to the valley exceeds outflow, there is potential recharge from stream losses to the groundwater system Water-level records included manual quarterly depth-to-water measurements collected from 68 wells, time-series data collected from 20 of those wells, and historic water levels from 16 wells. Hydrographs of the manual measurements showed declining water levels in 16 wells, mostly in the South-Main zone, and rising water levels in 14 wells, mostly in the Southern Ventucopa Uplands. Time-series hydrographs showed daily, seasonal, and longer-term effects associated with local pumping. Water-level data from the multiple-well monitoring sites indicated seasonal fluctuations as great as 80 feet and water-level differences between aquifers as great as 40 feet during peak pumping season. Hydrographs from the multiple-well groundwater monitoring sites showed vertical hydraulic gradients were upward during the winter months and downward during the irrigation season. Historic hydrographs showed water-level declines in the Southern-Main, Western Basin, Caliente Northern-Main, and Southern Sierra Madre zone ranging from 1 to 7 feet per year. Hydrographs of wells in the Southern Ventucopa Uplands zone showed several years with marked increases in water levels that corresponded to increased precipitation in the Cuyama Valley. Investigation of hydraulic properties included hydraulic conductivity and transmissivity estimated from aquifer tests performed on 63 wells. Estimates of horizontal hydraulic conductivity ranged from about 1.5 to 28 feet per day and decreased with depth. The median estimated hydraulic conductivity for the older alluvium was about five times that estimated for the Morales Formation. Estimates of transmissivity ranged from 560 to 163,400 gallons per day per foot and decreased with depth. The median estimated transmissivity for the younger alluvium was about three times that estimated for the older alluvium. Geomechanical analysis included land-surface elevation changes at five continuously operating global positioning systems (GPS) and land-subsidence detection at five interferometric synthetic aperture radar (InSAR) reference points. Analysis of data collected from continuously operating GPS stations showed the mountains to the south and west moved upward about 1 millimeter (mm) annually, whereas the station in the center of the Southern-Main zone moved downward more than 7 mm annually, indicating subsidence. It is likely that this subsidence is inelastic (permanent) deformation and indicates reduced storage capacity in the aquifer sediments. Analysis of InSAR data showed local and regional changes that appeared to be dependent, in part, on the time span of the interferogram, seasonal variations in pumping, and tectonic uplift. Long-term InSAR time series showed a total maximum detected subsidence rate of approximately 12 mm per year at one location and approximately 8 mm per year at a second location, while short-term InSAR time series showed maximum subsidence of about 15 mm at one location and localized maximum uplift of about 10 mm at another location.

California↗

Historical and projected climate (1901–2050) and hydrologic response of karst aquifers, and species vulnerability in south-central Texas and western South Dakota

Two karst aquifers, the Edwards aquifer in the Balcones Escarpment region of south-central Texas and the Madison aquifer in the Black Hills of western South Dakota, were evaluated for hydrologic response to projected climate change through 2050. Edwards aquifer sites include Barton Springs, the Bexar County Index Well, and Comal Springs. Madison aquifer sites include Spearfish Creek and Rhoads Fork Spring. Climate projections at sites were based on output from the Community Climate System Model of global climate, linked to the Weather Research and Forecasting (WRF) model of regional climate. The WRF model output was bias adjusted to match means for 1981–2010 from records at weather stations near Madison and Edwards aquifer sites, including Boerne, Texas, and Custer and Lead, South Dakota. Hydrologic response at spring and well sites was based on the Rainfall-Response Aquifer and Watershed Flow (RRAWFLOW) model. The WRF model climate projections for 2011–50 indicate a significant upward trend in annual air temperature for all three weather stations and a significant downward trend in annual precipitation for the Boerne weather station. Annual springflow simulated by the RRAWFLOW model had a significant downward trend for Edwards aquifer sites and no trend for Madison aquifer sites. Flora and fauna that rely on springflow from Edwards and Madison aquifer sites were assessed for vulnerability to projected climate change on the basis of the Climate Change Vulnerability Index (CCVI). The CCVI is determined by the exposure of a species to climate, the sensitivity of the species, and the ability of the species to cope with climate change. Sixteen species associated with springs and groundwater were assessed in the Balcones Escarpment region. The Barton Springs salamander ( Eurycea sosorum ) was scored as highly vulnerable with moderate confidence. Nine species—three salamanders, a fountain darter ( Etheostoma fonticola ), three insects, and two amphipods—were scored as moderately vulnerable. The remaining six species—four vascular plants, the Barton cavesnail ( Stygopyrgus bartonensis ), and a cave shrimp—were scored as not vulnerable/presumed stable (not vulnerable and evidence does not support change in abundance or range of the species). Vulnerability of eight species associated with streams that receive springflow from the Madison aquifer in the Black Hills was assessed. Of these, the American dipper ( Cinclus mexicanus ) and the lesser yellow lady’s slipper ( Cypripedium parviflorum ) were scored as moderately vulernable with high confidence. The dwarf scouringrush ( Equisetum scirpoides ) and autumn willow ( Salix serissima ) were also scored as moderately vulnerable with moderate to low confidence, respectively. Other species were assessed as not vulnerable/presumed stable or not vulnerable/increase likely (not vulnerable and evidence supporting an increase in abundance or range of the species). Lower vulnerability scores for the Black Hills species in comparison to the Balcones Escarpment species reflect lower endemicity, higher projected springflow than in the historical period, and high thermal tolerance of many of the species for the Black Hills. Importantly, climate change vulnerability scores differed substantially for Edwards aquifer species when RRAWFLOW model projections were included, resulting in increased vulnerability scores for 12 of the 16 species.

South Dakota, Texas↗

Hydrogeology and hydrology of the Punta Cabullones wetland area, Ponce, southern Puerto Rico, 2007-08

The U.S. Geological Survey, in cooperation with the Municipio Autónomo de Ponce and the Puerto Rico Department of Natural and Environmental Resources, conducted a study of the hydrogeology and hydrology of the Punta Cabullones area in Ponce, southern Puerto Rico. (Punta Cabullones is also referred to as Punta Cabullón.) The Punta Cabullones area is about 9 square miles and is an ecological system made up of a wetland, tidal flats, saltflats, mangrove forests, and a small fringing reef located a short distance offshore. The swales or depressions between successive beach ridges became development avenues for saline to hypersaline wetlands. The Punta Cabullones area was designated by the U.S. Fish and Wildlife Service as a coastal barrier in the 1980s because of its capacity to act as a buffer zone to ameliorate the impacts of natural phenomenon such as storm surges. Since 2003, Punta Cabullones has been set aside for preservation as part of the mitigation effort mandated by Federal and State laws to compensate for the potential environmental effects that might be caused by the construction of the Las Américas Transshipment Port. Total rainfall measured during 2008 within the Punta Cabullones area was 36 inches, which is slightly greater than the long-term annual average of 32 inches for the coastal plain near Ponce. Two evapotranspiration estimates, 29 and 37 inches, were obtained for the subarea of the Punta Cabullones area that is underlain by fan-delta and alluvial deposits by using two variants of the Penman semi-empirical equation. The long-term water stage and chemical character of the wetland in Punta Cabullones are highly dependent on the seasonal and annual variations of both rainfall and sea-wave activity. Also, unseasonal short-term above-normal rainfall and sea-wave events resulting from passing storms may induce substantial changes in the water stage and the chemical character of the wetland. In general, tidal fluctuations exert a minor role in modifying the water quality and stage of the wetland in Punta Cabullones. The role of the tidal fluctuations becomes important during those times when the outlets/inlets to the sea are not blocked by a sand bar and is allowed to freely flow into the wetland interior. The salinity of the wetland varies from brackish to hypersaline. The hypersaline conditions, including the occurrence of saltflats, within the Punta Cabullones wetland area result from a high evapotranspiration rate. The hypersaline conditions are further enhanced by a sand bar that blocks the inlet/outlet of the wetland’s easternmost channel, particularly during the dry season. Groundwater in Punta Cabullones mostly is present within beds of silisiclastic sand and gravel. During the study period, the depth to groundwater did not exceed 4 feet below land surface. The movement and direction of the groundwater flow in Punta Cabullones are driven by density variations that in turn result from the wide range of salinities in the groundwater. The salinity of the groundwater decreases within the first 60 to 100 feet of depth and decreases outward from a mound of hypersaline groundwater centered on piezometer nest PN2. The main groundwater types within the Punta Cabullones area vary from calcium-bicarbonate type in the northernmost part of the study area to a predominantly sodium-potassium-chloride groundwater type southward. According to stable-isotope data, groundwater within the study area is both modern meteoric water and seawater highly affected by evaporation. The chemical and stable-isotopic character of local groundwater is highly influenced by evapotranspiration because of its shallow depth. Equivalent freshwater heads indicate groundwater moves away from a mound centered on piezometer nest PN2, in a pattern similar to the spatial distribution of groundwater salinity. Vertical groundwater flow occurs in Punta Cabullones due to local differences in density. In the wetland subarea of Punta Cabullones, groundwater and surface water are hydraulically coupled. Locally, surface-hypersaline water sinks into the aquifer, providing recharge and serving as a mechanism to redistribute salt throughout the study area. The evapotranspiration in the wetland subarea is estimated at about 11 million gallons per day (Mgal/d) that is equivalent to about 12,586 acre-feet per year. The balance of evapotranspiration, in excess of the about 0.5 Mgal/d of groundwater flow within the wetland, is supplied by saline to hypersaline surface water that may include seawater and meteoric water highly affected by evaporation with dissolved salts. In one of the extreme scenarios in which no groundwater is intercepted by pumpage at the Restaurada well field, the amount of saline to hypersaline water in the wetland consumed by evapotranspiration is about 10.5 Mgal/d. In the opposite extreme in which the entire regional groundwater flow is intercepted by pumpage in the Restaurada well field, the entire evapotranpiration requirement is met by saline to hypersaline water. Hydrologic, isotopic, and chemical data indicate that all of, or a large portion of, the historical groundwater flow to Punta Cabullones is being captured by the Puerto Rico Aqueducts and Sewer Authority pumpage at the Restaurada well field at a rate of about 2 Mgal/d. As a consequence, seawater intrusion into the aquifer at the Punta Cabullones area seems to be occurring, while the current pumpage at the Restaurada well field is sustained by storage depletion of the aquifer.

Puerto Rico↗

Hydrologic budget and conditions of Permian, Pennsylvanian, and Mississippian aquifers in the Appalachian Plateaus physiographic province

In response to challenges to groundwater availability posed by historic land-use practices, expanding development of hydrocarbon resources, and drought, the U.S. Geological Survey Groundwater Resources Program began a regional assessment of the Appalachian Plateaus aquifers in 2013 that incorporated a hydrologic landscape approach to estimate all components of the hydrologic system: surface runoff, base flow from groundwater, and interaction with atmospheric water (precipitation and evapotranspiration). This assessment was intended to complement other Federal and State investigations and provide foundational groundwater-related datasets in the Appalachian Plateaus. A regional Soil-Water-Balance model was constructed for a 160,000-square-mile study area that extended to the topographic divide of all streams originating outside but flowing into areas underlain by Appalachian Plateaus aquifers. The model incorporated soil, landscape, and climate variables to estimate an annual water budget for the 32-year period from 1980 to 2011 and was calibrated using base-flow data estimated by hydrograph separation techniques from 20 streamflow gaging stations across the study area. Over this period, an average of 47 inches per year (in/yr) of precipitation fell on Appalachian Plateaus aquifers. Simulations from the regional Soil-Water-Balance model indicate that only 19 percent of the precipitation or an average 9 in/yr recharged aquifers, and 19 percent resulted in surface runoff to streams. The remaining 62 percent, an average of 27 in/yr of water, was returned to the atmosphere via evapotranspiration. Because withdrawals from aquifers due to pumping equated to less than 1 percent of the water budget, differences in predevelopment and postdevelopment regional water budgets of the Appalachian Plateaus were minimal. Storage changes caused by filling of abandoned coal-mine aquifers and long-term differences in aquifer storage resulting from climate fluctuations constitute a small portion of the overall water budget. The percentage of precipitation that results in recharge, runoff, or evapotranspiration from the landscape varies annually by up to a factor of two depending on temporal changes in prevailing climate conditions and spatial changes in basin characteristics, precipitation patterns, and sources of atmospheric moisture over a large study area. A comparison of water-budget estimates from the regional Soil-Water-Balance model for a dry year (1988) and wet year (2004) showed that evapotranspiration accounts for most of the annual differences in precipitation. As a portion of annual precipitation, evapotranspiration ranged from 69 percent (dry year) to 52 percent (wet year), a range four times greater than the 15 percent (dry year) to 18 percent (wet year) range estimated for recharge. Evapotranspiration as a percentage of precipitation peaks during dry periods, whereas base flow and runoff tend to reach minimum values. During wet periods, this relationship is reversed and base flow and runoff as a percentage of precipitation generally peak while evapotranspiration percentages reach minimum values. Annual recharge in the Appalachian Plateaus reaches a maximum at near 20 percent of annual precipitation, regardless of the severity of wet conditions. Hydrograph separation data from 849 streamflow gaging stations in the study area were used to assess trends in streamflow, base flow, surface runoff, and base-flow index, or ratio of base flow to streamflow, in the Appalachian Plateaus for the period from 1930 to 2011. Annual data anomalies for each of the four variables were individually defined as the annual standard deviation from the mean at all 849 streamflow gaging stations. Annual data anomalies confirm the close relation of annual precipitation to both base flow and runoff components of streamflow, and both components increased during the period of analysis. Around 1970, conditions shifted streamflow from values generally below to above long-term means. At a regional scale, increases in base flow account for most of these observed increases in mean annual streamflow. The independence of the base-flow index to annual climate trends indicate that changes in the components of streamflow of the Appalachian Plateaus are probably in response to shifts in seasonal precipitation or widespread land-use practices. A subset of 77 index streamgages, defined as having 60 or more years of complete record between the years 1930 and 2011 with no more than 20 percent missing data, was selected to show spatial patterns of change in the water budget. Data from the index streamgages showed that the overall trends in base flow are dependent upon the period of evaluation. Long-term (1930–2011) increases in base flow were observed throughout the study area. For two shorter periods (1930–1969 and 1970–2011) trends in base flow were largely negative. In general, spatial patterns of change in streamflow, base flow, and runoff were mixed but generally consistent with prevailing climate patterns and land-use changes.

Alabama, Kentucky, Maryland, Ohio, Pennslyvania, V↗

Hydrology and water quality in 13 watersheds in Gwinnett County, Georgia, 2001–15

The U.S. Geological Survey (USGS), in cooperation with Gwinnett County Department of Water Resources, established a Long-Term Trend Monitoring (LTTM) program in 1996. The LTTM program is a comprehensive, long-term, water-quantity and water-quality monitoring program designed to document and analyze the hydrologic and water-quality conditions of selected watersheds in Gwinnett County, Georgia. Water-quality monitoring initially began in six watersheds and currently [2016] includes 13 watersheds. As part of the LTTM program, streamflow, precipitation, water temperature, specific conductance, and turbidity were measured every 15 minutes for water years 2001–15 at 12 of the 13 watershed monitoring stations and for water years 2010–15 at the other watershed. In addition, discrete water-quality samples were collected seasonally from May through October (summer) and November through April (winter), including one base-flow and three stormflow event composite samples, during the study period. Samples were analyzed for nutrients (nitrogen and phosphorus), total organic carbon, trace elements (total lead and total zinc), total dissolved solids, and total suspended sediment (total suspended solids and suspended-sediment concentrations). The sampling scheme was designed to identify variations in water quality both hydrologically and seasonally. The 13 watersheds were characterized for basin slope, population density, land use for 2012, and the percentage of impervious area from 2000 to 2014. Several droughts occurred during the study period—water years 2002, 2007–08, and 2011–12. Watersheds with the highest percentage of impervious areas had the highest runoff ratios, which is the portion of precipitation that occurs as runoff. Watershed base-flow indexes, the ratio of base-flow runoff to total runoff, were inversely correlated with watershed impervious area. Flood-frequency estimates were computed for 13 streamgages in the study area that have 10 or more years of annual peak flow data through water year 2015, using the expected moments algorithm to fit a Pearson Type III distribution to logarithms of annual peak flows. Kendall’s tau nonparametric test was used to determine the statistical significance of trends in the annual peak flows, with none of the 13 streamgages exhibiting significant trends. A comparison of base-flow and stormflow water-quality samples indicates that turbidity and concentrations of total ammonia plus organic nitrogen, total nitrogen, total phosphorus, total organic carbon, total lead, total zinc, total suspended solids, and suspended-sediment concentrations increased with increasing discharge at all watersheds. Specific conductance decreased during stormflow at all watersheds, and total dissolved solids concentrations decreased during stormflow at a few of the watersheds. Total suspended solids and suspended-sediment concentrations typically were two orders of magnitude higher in stormflow samples, turbidities were about 1.5 orders of magnitude higher, total phosphorus and total zinc were about one order of magnitude higher, and total ammonia plus organic nitrogen, total nitrogen, total organic carbon, and total lead were about twofold higher than in base-flow samples. Seasonality and long-term trends were identified for the period water years 2001–15 for 10 constituents—total nitrogen, total nitrate plus nitrite, total phosphorus, dissolved phosphorus, total organic carbon, total suspended solids, suspended-sediment concentration, total lead, total zinc, and total dissolved solids. Seasonal patterns were present in most watersheds for all constituents except total dissolved solids, and the watersheds had fairly similar patterns of higher concentrations in the summer and lower concentrations in the winter. A linear long-term trend analysis of residual concentrations from the flow-only load estimation model (without time-trend terms) identified significant trends in 67 of the 130 constituent-watershed combinations. Seventy percent of the significant trends were negative. Total organic carbon and total dissolved solids had predominantly positive trends. Total phosphorus, total suspended solids, suspended-sediment concentration, total lead, and total zinc had only negative trends. The other three constituents exhibited fewer trends, both positive and negative. Streamwater loads were estimated annually for the 13-year period water years 2003–15 for the same 10 constituents in the trend analysis. Loads were estimated using a regression-model-based approach developed by the USGS for the Gwinnett County LTTM program that accommodates the use of storm-event composited samples. Concentrations were modeled as a function of discharge, base flow, time, season, and turbidity to improve model predictions and reduce errors in load estimates. Total suspended solids annual loads have been identified in Gwinnett County’s Watershed Protection Plan for target performance criterion. Although the amount of annual runoff was the primary factor in variations in annual loads, climatic conditions (classified as dry, average, or wet) affected annual loads beyond what was attributed to climatic-related variations in annual runoff. Significant negative trends in loads were estimated for the combined area of the watersheds for all constituents except dissolved phosphorus, total organic carbon, and total dissolved solids. The trend analysis indicated that total suspended solids and suspended-sediment concentration loads in the study area were decreasing by 57,000 and 87,000 pounds per day per year, respectively. Variations in constituent yields between watersheds appeared to be related to various watershed characteristics. Suspended sediment (as either total suspended solids or suspended-sediment concentrations), along with constituents transported predominately in solid phase (total phosphorus, total organic carbon, total lead, and total zinc), and total dissolved solids typically had higher yields from watersheds that had high percentages of impervious areas or high basin slope. High total nitrogen yields were also associated with watersheds with high percentages of impervious areas. Low total nitrogen, total suspended solids, total lead, and total zinc yields appeared to be associated with watersheds that had a low percentage of high-density development.

Georgia↗

Hydrology of the Claiborne aquifer and interconnection with the Upper Floridan aquifer in southwest Georgia

The U.S. Geological Survey conducted a study, in cooperation with the Georgia Environmental Protection Division, to define the hydrologic properties of the Claiborne aquifer and evaluate its connection with the Upper Floridan aquifer in southwest Georgia. The effort involved collecting and compiling hydrologic data from the aquifer in subarea 4 of southwestern Georgia. Data collected for this study include borehole geophysical logs in 7 wells, and two 72-hour aquifer tests to determine aquifer properties. The top of the Claiborne aquifer extends from an altitude of about 200 feet above the North American Vertical Datum of 1988 (NAVD 88) in Terrell County to 402 feet below NAVD 88 in Decatur County, Georgia. The base of the aquifer extends from an altitude of about 60 feet above NAVD 88 in eastern Sumter County to about 750 feet below NAVD 88 in Decatur County. Aquifer thickness ranges from about 70 feet in eastern Early County to 400 feet in Decatur County. The transmissivity of the Claiborne aquifer, determined from two 72-hour aquifer tests, was estimated to be 1,500 and 700 feet squared per day in Mitchell and Early Counties, respectively. The storage coefficient was estimated to be 0.0006 and 0.0004 for the same sites, respectively. Aquifer test data from Mitchell County indicate a small amount of leakage occurred during the test. Groundwater-flow models suggest that the source of the leakage was the underlying Clayton aquifer, which produced about 2.5 feet of drawdown in response to pumping in the Claiborne aquifer. The vertical hydraulic conductivity of the confining unit between the Claiborne and Clayton aquifers was simulated to be about 0.02 foot per day. Results from the 72-hour aquifer tests run for this study indicated no interconnection between the Claiborne and overlying Upper Floridan aquifers at the two test sites. Additional data are needed to monitor the effects that increased withdrawals from the Claiborne aquifer may have on future water resources.

Georgia↗

Hydrologic assessment of the Edwin B. Forsythe National Wildlife Refuge

The Edwin B. Forsythe National Wildlife Refuge (hereafter Forsythe refuge or the refuge) is situated along the central New Jersey coast and provides a mixture of freshwater and saltwater habitats for numerous bird, wildlife, and plant species. Little data and information were previously available regarding the freshwater dynamics that support the refuge’s ecosystems. In cooperation with the U.S. Fish and Wildlife Service, the U.S. Geological Survey conducted an assessment of the hydrologic resources and processes in the refuge and surrounding areas to provide baseline information for evaluating restoration projects and future changes in the hydrologic system associated with climate change and other anthropogenic stressors. During spring 2015, water levels were measured at groundwater and surface-water sites in and near the Forsythe refuge. These water-level measurements, along with surface-water elevations obtained from digital elevation models, were used to construct water-table-elevation and depth-to-water maps of the refuge and surrounding areas. Water-table elevations in the refuge ranged from sea level to approximately 65 feet above sea level; in most of the refuge, the water-table elevation was within 3 feet of sea level. The water-table-elevation map indicates that the direction of shallow groundwater flow at the regional scale is generally from west to east (much of it from the northwest to the southeast), and groundwater moves downgradient from the uplands toward major groundwater discharge areas consisting of coastal streams and wetlands. The depth to water is estimated to be less than 2 feet for approximately 86 percent of the refuge, which coincides closely with the percentage of wetland area in the refuge. Depth to water in excess of 20 feet below land surface is limited to higher elevation areas of the refuge. Streamflow data collected at continuous-record streamgages and partial-record stations within the Mullica-Toms Basin were summarized. Hydrograph separation of streamflow data for eight streamgages (2004–13) reveals that base flow accounts for 68–94 percent of streamflow in basins upstream from the refuge. The high base-flow inputs underscore the importance of groundwater as a source of freshwater that supports both the streams that flow into the refuge and the hydroecology of the contributing basins. Mean annual flow typically ranged from 1.7 to 2.1 cubic feet per second per square mile at the streamgages (2004–13) and between 1.2 and 2.3 cubic feet per second per square mile at the partial-record stations (1965–2015) but was notably greater or lower than these ranges at several stations. Mean annual water budgets were estimated for multiple regions of the refuge for 2004–13 using data compiled from nearby meteorological stations and groundwater flows derived from previously calibrated groundwater-flow models. Precipitation, groundwater recharge, and evapotranspiration were estimated from available data; direct runoff was calculated as the residual component of the water balance. Groundwater recharge rates were greatest in the upland-dominated areas of the refuge with estimates of 14.4 to 18.9 inches per year, which are equivalent to 30 to 40 percent of precipitation. Groundwater recharge rates were nearly zero in the central coastal areas because these areas are major groundwater discharge zones, the water table is near land surface, the subsurface is close to saturation and cannot accept much recharge, and much of the area is underlain by thick marsh deposits likely with low permeability. Estimates of evapotranspiration varied from about 26 inches per year in the upland-dominated areas to more than 35 inches per year in the coastal wetlands, equivalent to 55–79 percent of mean annual precipitation, indicating that it is a major component of the hydrodynamics of the Forsythe refuge. On the basis of output from previously calibrated groundwater-flow models, nearly all of the groundwater exiting the surficial aquifer system in the central coastal areas of the refuge is discharged to wetlands, which highlights the importance of groundwater discharge in supporting the ecosystems of the Forsythe refuge. In the central coastal areas, horizontal flow contributes more than 90 percent of the groundwater flow to the surficial system, indicating that the upbasin areas are a substantial source of water that ultimately discharges to streams and wetlands in the refuge.

New Jersey↗

Hydrogeologic setting, conceptual groundwater flow system, and hydrologic conditions 1995–2010 in Florida and parts of Georgia, Alabama, and South Carolina

The hydrogeologic setting and groundwater flow system in Florida and parts of Georgia, Alabama, and South Carolina is dominated by the highly transmissive Floridan aquifer system. This principal aquifer is a vital source of freshwater for public and domestic supply, as well as for industrial and agricultural uses throughout the southeastern United States. Population growth, increased tourism, and increased agricultural production have led to increased demand on groundwater from the Floridan aquifer system, particularly since 1950. The response of the Floridan aquifer system to these stresses often poses regional challenges for water-resource management that commonly transcend political or jurisdictional boundaries. To help water-resource managers address these regional challenges, the U.S. Geological Survey (USGS) Water Availability and Use Science Program began assessing groundwater availability of the Floridan aquifer system in 2009. The current conceptual groundwater flow system was developed for the Floridan aquifer system and adjacent systems partly on the basis of previously published USGS Regional Aquifer-System Analysis (RASA) studies, specifically many of the potentiometric maps and the modeling efforts in these studies. The Floridan aquifer system extent was divided into eight hydrogeologically distinct subregional groundwater basins delineated on the basis of the estimated predevelopment (circa 1880s) potentiometric surface: (1) Panhandle, (2) Dougherty Plain-Apalachicola, (3) Thomasville-Tallahassee, (4) Southeast Georgia-Northeast Florida-South South Carolina, (5) Suwannee, (6) West-central Florida, (7) East-central Florida, and (8) South Florida. The use of these subregions allows for a more detailed analysis of the individual basins and the groundwater flow system as a whole. The hydrologic conditions and associated groundwater budget were updated relative to previous RASA studies to include additional data collected since the 1980s and to reflect the entire groundwater flow system, including the surficial, intermediate, and Floridan aquifer systems for a contemporary period (1995–2010). Inflow to the groundwater flow system of 33,700 million gallons per day (Mgal/d) was assumed to be exclusively from net recharge (precipitation minus evapotranspiration and surface runoff). Outflow from the groundwater flow system included spring discharge (7,700 Mgal/d) and groundwater withdrawals (5,200 Mgal/d). Estimates for all components of the groundwater system were not possible because of large uncertainties associated with internal leakage, coastal discharge, and discharge to streams and lakes. A numerical modeling analysis is required to improve this hydrologic budget calculation and to forecast future changes in groundwater levels and aquifer storage caused by groundwater withdrawals, land-use change, and the effects of climate variability and change.

Alabama, Florida, Georgia, South Carolina↗

Effects of surface-water and groundwater inflows and outflows on the hydrology of the Tsala Apopka Lake Basin in Citrus County, Florida

The U.S. Geological Survey, in cooperation with the Southwest Florida Water Management District, initiated a study to quantify the inflows and outflows in the Floral City, Inverness, and Hernando pools of the Tsala Apopka Lake Basin in Citrus County, Florida. This study assesses hydrologic changes in pool stages, groundwater levels, spring flows, and streamflows caused by the diversion of streamflow from the Withlacoochee River to the Tsala Apopka Lake Basin through water-control structures. A surface-water/groundwater flow model was developed using hydraulic parameters for lakes, streams, the unsaturated zone, and the underlying surficial and Upper Floridan aquifers estimated using an inverse modeling calibration technique. After calibration, the model was used to assess the relation between inflows and outflows in the Tsala Apopka Lake Basin and changes in pool stages. Simulation results using the calibrated surface-water/groundwater flow model showed that leakage rates from the pools to the Upper Floridan aquifer were largest at the deep lake cells and that these leakage rates to the Upper Floridan aquifer were the highest in the model area. Downward leakage to the Upper Floridan aquifer occurred beneath most of the extent of the Floral City, Inverness, and Hernando pools. These leakage rates depended on the lakebed leakance and the difference between lake stages and heads in the Upper Floridan aquifer. Leakage rates were higher for the Floral City pool than for the Inverness pool, and higher for the Inverness pool than for the Hernando pool. Lakebed leakance was higher for the Floral City pool than for the Hernando pool, and higher for the Hernando pool than for the Inverness pool. Simulation results showed that the average recharge rate to the surficial aquifer was 10.3 inches per year for the 2004 to 2012 simulation period. Areas that recharge the surficial aquifer covered about 86 percent of the model area. Simulations identified areas along segments of the Withlacoochee River and within land-surface depressions that receive water from the surficial aquifer. Recharge rates were largest in physiographic regions having a deep water table. Simulated heads in the Upper Floridan aquifer indicated the general flow directions in the active flow model area were from the northeast toward the southwest and then westward toward the coast, and from the southeast toward the northwest and then westward toward the coast, consistent with flow directions inferred from the estimated potentiometric surface map for May 2010. The largest inflow in the water budget of the Upper Floridan aquifer was downward leakage from the overlying hydrogeologic unit. The largest outflow in the water budget of the Upper Floridan aquifer was spring flow. The calibrated surface-water and groundwater flow model was used to simulate hydrologic scenarios that included changes in rainfall rates, projected increases in groundwater pumping rates for 2025 and 2035, no flow for the 2004–12 period through the eight water-control structures in the Tsala Apopka Lake Basin, and the removal of the Inglis Dam and the Inglis Bypass Spillway on Lake Rousseau. Scenario simulation results were compared to annual average calibrated water levels and flows from 2004 to 2012. Simulated declines in the Tsala Apopka Lake pool stages under the 10-percent lower rainfall scenario were about 0.8, 0.3, and 1.3 feet (ft) for the Floral City, Inverness, and Hernando pools, respectively. Simulated groundwater levels under the same scenario declined up to 5.4 ft in the surficial aquifer and up to 2.9 ft in the Upper Floridan aquifer. Under the projected increases in groundwater pumping rates for 2035 that represented an increase of 36 percent from average 2004 to 2012 pumping rates, the simulated declines in the Floral City, Inverness, and Hernando pool stages were, in downstream order, 0.02, 0.06, and 0.04 ft. The largest drawdown under the projected increases in groundwater pumping rates for 2035 was 2.1 ft in the surficial aquifer and about 1.8 ft in the Upper Floridan aquifer. A scenario of decreased rainfall by 10 percent caused greater declines in water levels and pool stages than projected increases in groundwater pumping rates. The simulation with no flow through the eight Tsala Apopka Lake water-control structures resulted in simulated declines in average pool stage of 1.8, 1.9, and 0.5 ft in the Floral City, Inverness, and Hernando pools, respectively. The simulated removal of the two water-control structures in Lake Rousseau caused flow to increase at Rainbow Springs by 28 cubic feet per second, an increase of 4.7 percent from the average calibrated flow for 2004 to 2012.

Florida↗

Hydrologic and hydraulic analyses of selected streams in Richland County, Ohio

Hydrologic and hydraulic analyses were done for selected reaches of Clear Fork Mohican River and Cedar Fork in Richland County, Ohio. To update and expand a portion of the Federal Emergency Management Agency detailed Flood Insurance Study, the U.S. Geological Survey (USGS) and the Muskingum Watershed Conservancy District initiated a cooperative study. The study comprised an 18.6-mile reach of the Clear Fork Mohican River and a 5.9-mile reach of Cedar Fork. Historical streamflow data from the streamgage Clear Fork Mohican River at Bellville, Ohio (USGS station number 03131982) and regional regression equations were used to estimate instantaneous peak streamflows for floods with 10-, 4-, 2-, 1-, and 0.2-percent and 1-percent plus annual exceedance probabilities. The 1-percent plus flood elevation is defined by the Federal Emergency Management Agency as a flood elevation derived by using streamflows that include the average predictive error for the regression equation streamflow calculation for the Flood Risk project. This error is then added to the 1-percent annual exceedance probability flood streamflow to calculate the new 1-percent plus streamflow. The annual exceedance probability streamflows were then used in a Hydrologic Engineering Center-River Analysis System step-backwater model to determine water-surface elevation profiles and flood-inundation boundaries for the 10-, 4-, 2-, 1-, and 0.2-percent and 1-percent plus annual exceedance probability floods and a regulatory floodway along a selected reach of each stream. The Clear Fork Mohican River model was calibrated to 16 flood events by using the current stage-streamflow relation at the streamgage Clear Fork Mohican River at Bellville, Ohio (USGS station number 03131982) and a submersible pressure transducer. Flood-inundation boundaries for the 1- and 0.2-percent annual exceedance probability floods and a regulatory floodway were mapped for each stream.

Ohio↗

Hydrologic study at Farm Creek Marsh, Dorchester County, Maryland, from April 2015 to April 2016

In 2015, the U.S. Geological Survey began a 1-year hydrologic study to investigate the extent and cause of inundation at Farm Creek Marsh, in Dorchester County, Maryland. In combination with a tide and precipitation gage, a representative section of the marsh was instrumented with surface-water monitors and shallow groundwater piezometers to capture the spatial and temporal extent of inundation. In addition, water-quality data (major ions and nutrients) were collected to help discern the cause of inundation. Results indicate that during the year-long study, all sites were periodically inundated, ranging from a total of 108 days to the entire study period of 353 days. The depth of inundation was typically between 0 and 0.2 feet (ft) (above land surface), with the exception of large storm events. Less than 0.5 ft of elevation was the difference between a site being inundated during the entire study period of 353 days and a site being inundated for 36 consecutive days out of 108 total days of inundation during the study period. Water-quality data showed a large difference in pH between marsh surface water (6.1 to 6.9 standard pH units) and shallow groundwater (3.0 to 3.6 standard pH units), with differences also observed in concentrations of silica, iron, manganese, and potassium. Collectively, the combination of water-quality, hydrologic, and soils data indicate that inundation is caused by tide and storm events rather than groundwater discharge.

Maryland↗