SEARCH · Search USGS

Results for “Land and Reservoir Management”

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

SourceResearch connectionFrom yearThrough year

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

212 records · Page 12Linked to original sources

USGS2014

Hydrogeologic framework and groundwater/surface-water interactions of the upper Yakima River Basin, Kittitas County, central Washington

The hydrogeology, hydrology, and geochemistry of groundwater and surface water in the upper (western) 860 square miles of the Yakima River Basin in Kittitas County, Washington, were studied to evaluate the groundwater-flow system, occurrence and availability of groundwater, and the extent of groundwater/surface-water interactions. The study area ranged in altitude from 7,960 feet in its headwaters in the Cascade Range to 1,730 feet at the confluence of the Yakima River with Swauk Creek. A west-to-east precipitation gradient exists in the basin with the western, high-altitude headwaters of the basin receiving more than 100 inches of precipitation per year and the eastern, low-altitude part of the basin receiving about 20 inches of precipitation per year. From the early 20th century onward, reservoirs in the upper part of the basin (for example, Keechelus, Kachess, and Cle Elum Lakes) have been managed to store snowmelt for irrigation in the greater Yakima River Basin. Canals transport water from these reservoirs for irrigation in the study area; additional water use is met through groundwater withdrawals from wells and surface-water withdrawals from streams and rivers. Estimated groundwater use for domestic, commercial, and irrigation purposes is reported for the study area. A complex assemblage of sedimentary, metamorphic, and igneous bedrock underlies the study area. In a structural basin in the southeastern part of the study area, the bedrock is overlain by unconsolidated sediments of glacial and alluvial origin. Rocks and sediments were grouped into six hydrogeologic units based on their lithologic and hydraulic characteristics. A map of their extent was developed from previous geologic mapping and lithostratigraphic information from drillers’ logs. Water flows through interstitial space in unconsolidated sediments, but largely flows through fractures and other sources of secondary porosity in bedrock. Generalized groundwater-flow directions within the unconfined part of the aquifers in unconsolidated sediments indicate generalized groundwater movement toward the Yakima River and its tributaries and the outlet of the study area. Groundwater movement through fractures within the bedrock aquifers is complex and varies over spatial scales depending on the architecture of the fracture-flow system and its hydraulic properties. The complexity of the fracturedbedrock groundwater-flow system is supported by a wide range of groundwater ages determined from geochemical analyses of carbon-14, sulfur hexafluoride, and tritium in groundwater. These geochemical data also indicate that the shallow groundwater system is actively flushing with young, isotopically heavy groundwater, but isotopicallylight, Pleistocene-age groundwater with a geochemicallyevolved composition occurs at depth within the fracturedbedrock aquifers of upper Kittitas County. An eastward depletion of stable isotopes in groundwater is consistent with hydrologically separate subbasins. This suggests that groundwater that recharges in one subbasin is not generally available for withdrawal or discharge into surface-water features within other subbasins. Water budget components were calculated for 11 subbasins using a watershed model and varied based on the climate, land uses, and geology of the subbasin. Synoptic streamflow measurements made in August 2011 indicate that groundwater discharges into several tributaries of the Yakima River with several losses of streamflow measured where the streams exit bedrock uplands and flow over unconsolidated sediments. Profiles of stream temperature during late summer suggest cool groundwater inflow over discrete sections of streams. This groundwater/surfacewater connection is further supported by the stable-isotope composition of stream water, which reflects the local stableisotope composition of groundwater measured at some wells and springs. Collectively, these hydrogeologic, hydrologic, and geochemical data support a framework for evaluating the potential effects of future groundwater appropriations on senior surface-water and groundwater rights and streamflows. Although total pumping rates in upper Kittitas County of about 3.5 cubic feet per second are small relative to other components of the water budget, the magnitude, timing, and location of withdrawals may have important effects on the hydrologic system. The heterogeneous and variably fractured bedrock in the study area precluded a detailed evaluation of localized effects of pumping, but several generalizations about the groundwater and surface-water systems can be made. These generalizations include evidence for the continuity between the groundwater and surface-water system apparent from synoptic streamflow measurements, stream-temperature profiles, and stable-isotope data of groundwater and surface waters.

Washington↗

USGS2001

A spatially referenced regression model (SPARROW) for suspended sediment in streams of the Conterminous U.S.

Suspended sediment has long been recognized as an important contaminant affecting water resources. Besides its direct role in determining water clarity, bridge scour and reservoir storage, sediment serves as a vehicle for the transport of many binding contaminants, including nutrients, trace metals, semi-volatile organic compounds, a nd numerous pesticides (U.S. Environmental Protection Agency, 2000a). Recent efforts to addr ess water-quality concerns through the Total Maximum Daily Load (TMDL) process have iden tified sediment as the single most prevalent cause of impairment in the Nation’s streams a nd rivers (U.S. Environmental Protection Agency, 2000b). Moreover, sediment has been identified as a medium for the tran sport and sequestration of organic carbon, playing a potentia lly important role in understa nding sources and sinks in the global carbon budget (Stallard, 1998). A comprehensive understanding of sediment fate a nd transport is considered essential to the design and implementation of effective plans for sediment management (Osterkamp and others, 1998, U.S. General Accounting Office, 1990). An exte nsive literature addr essing the problem of quantifying sediment transport has produced a nu mber of methods for estimating its flux (see Cohn, 1995, and Robertson and Roerish, 1999, for us eful surveys). The accuracy of these methods is compromised by uncertainty in the concentration measurements and by the highly episodic nature of sediment movement, particul arly when the methods are applied to smaller basins. However, for annual or decadal flux es timates, the methods are generally reliable if calibrated with extended periods of data (Robertson and Roerish, 1999). A substantial literature also supports the Universal Soil Loss Equation (U SLE) (Soil Conservation Service, 1983), an engineering method for estimating sheet and rill erosion, although the empirical credentials of the USLE have recently been questioned (Tri mble and Crosson, 2000). Conversely, relatively little direct evidence is available concerning the fate of sediment. The common practice of quantifying sediment fate with a sediment deliv ery ratio, estimated from a simple empirical relation with upstream basin area, does not artic ulate the relative importance of individual storage sites within a basin (Wolman, 1977). Rates of sediment deposition in reservoirs and flood plains can be determined from empirical measurement s , but only a limited number of sites have been monitored, and net rates of deposition or loss from other potential sinks and sources is largely unknown (Stallard, 1998). In particular, little is known about how much sediment loss from fields ultimately makes its way to stream channels, and how much sediment is subsequently stored in or lost from th e streambed (Meade and Parker, 1985, Trimble and Crosson, 2000). This paper reports on recent progress made to a ddress empirically the question of sediment fate and transport on a national scale. The model pres ented here is based on the SPAtially Referenced Regression On Watershed attr ibutes (SPARROW) methodology, fi rst used to estimate the distribution of nutrients in str eams and rivers of the United Stat es, and subsequently shown to describe land and stream processes affecting the delivery of nutrients (Smith and others, 1997, Alexander and others, 2000, Preston and Brakeb ill, 1999). The model makes use of numerous spatial datasets, available at the national level, to explain long-term sediment water-quality conditions in major streams and rivers throughou t the United States. Sediment sources are identified using sediment erosion rates from the National Resources I nventory (NRI) (Natural Resources Conservation Service, 2000) and apportioned over the landscape according to 30- meter resolution land-use information from th e National Land Cover Data set (NLCD) (U.S. Geological Survey, 2000a). More than 76,000 reservoirs from the National Inventory of Dams (NID) (U.S. Army Corps of Engin eers, 1996) are identified as pot ential sediment sinks. Other, non-anthropogenic sources and sinks are identified using soil in formation from the State Soil Survey Geographic (STATSGO) data base (Schwarz and Alexander, 1995) and spatial coverages representing surficial rock t ype and vegetative cover. The SPA RROW model empirically relates these diverse spatial datasets to estimates of long-term, mean annual sediment flux computed from concentration and flow measurements co llected over the period 1985 -95 from more than 400 monitoring stations maintained by the Na tional Stream Quality Accounting Network (Alexander and others, 1998), the National Wa ter Quality Assessment Program, and U.S. Geological Survey District offices (Turcios and Gray, in press). Th e calibrated model is used to estimate sediment flux for over 60,000 stream segments included in the River Reach File 1 (RF1) stream network (Alexander and others, 1999). SPARROW uses statis tical methods to calibrate a simple, structural model of riverine water quality, one that imposes mass ba lance in accounting for changes in contaminant flux. As applied here, the mass-balance approach facilitates the interpretation of model results in terms of physical processes affecting sediment transport, and makes possible the estimation of various rates of sediment generation and loss associated with stream channels and features of the landscape. The statistical approach provides a basi s for assessing the error of these inferred rates and of the error in extrapolated estimates of sediment flux made for streams in the RF1 network. An important implication of the holistic modeling approach adopted in this analysis is that estimates of sediment production and loss ar e based on, and therefore consistent with, measurements of in-stream flux. Other ancillary information, such as direct measurements of long-term sediment storage and release from rese rvoirs (Steffen, 1996), is incorporated into the analysis by specifying additional equations expl aining these ancillary variables. The imposition of cross-equation constraints affords this info rmation a statistically consistent weight in explaining in-stream sediment flux. Thus, the me thodology described here represents a general framework for synthesizing a wide spectrum of available information relevant to the understanding of sediment fate and transport.

Conterminous United States↗

USGS1973

Water-supply development and management alternatives for Clinton, Eaton, and Ingham Counties, Michigan

The Tri-County region, consisting of Clinton, Eaton, and Ingham Counties, is an area of 1,697 square miles in Michigan's Lower Peninsula and has as its hub the Lansing metropolitan area. The land surface ranges in altitude from about 700 to about 1,000 feet. The region receives an average of about 31 inches of precipitation each year. The population is nearing 400,000 and by 1990 will be near 600,000. Average daily water use is slightly more than 30 million gallons today; by 1980 it will be about 50 million gallons, and by 1990 it will probably be about 70 million gallons. The Tri-County region is drained by seven river systems. The median annual 7 -day mean low flows of the principal streams in these systems were measured at the point farthest downstream within the region. These values, in cubic feet per second, are as follows: Grand River, 180; Maple River, 34; Looking Glass River, 28; Red Cedar River, 30; Portage Creek, 15; Battle Creek, 20; and Thornapple River, 24-a total of 331 cubic feet per second or about 220 million gallons per day. The areal variance in 7-day low-flow runoff ranges from 0 to 0.15 cubic foot per second per square mile. The principal source of ground water in the Tri-County region is a complex aquifer system composed of the Saginaw and Grand River Formations and some of the overlying glacial sediments. This aquifer yields between 300 and 700 gallons per minute to individual wells in much of the western half of Ingham County, in the eastern half of Clinton County, in a small area in southeastern Clinton County, and in northeastern Eaton County. In some parts of the region, the glacial sediments are favorable for development of moderate to large supplies of water. Minor aquifers in the region are the Bayport, Michigan, and Marshall Formations. Providing water supplies in the future requires complete and comprehensive water-management programs. Such management programs involve determining which of several alternative water-development systems is the best. Some of the chief factors and methods that must be considered when planning these systems are combined use of ground and surface water, artificial recharge, treatment of wastes, use of storage reservoirs, and importation of water from the Great Lakes.

Michigan↗

USGS2006

Hydrogeologic framework refinement, ground-water flow and storage, water-chemistry analyses, and water-budget components of the Yuma area, southwestern Arizona and southeastern California

The ground-water and surface-water system in the Yuma area in southwestern Arizona and southeastern California is managed intensely to meet water-delivery requirements of customers in the United States, to manage high ground-water levels in the valleys, and to maintain treaty-mandated water-quality and quantity requirements of Mexico. The following components in this report, which were identified to be useful in the development of a ground-water management model, are: (1) refinement of the hydrogeologic framework; (2) updated water-level maps, general ground-water flow patterns, and an estimate of the amount of ground water stored in the mound under Yuma Mesa; (3) review and documentation of the ground-water budget calculated by the Bureau of Reclamation, U.S. Department of the Interior (Reclamation); and (4) water-chemistry characterization to identify the spatial distribution of water quality, information on sources and ages of ground water, and information about the productive-interval depths of the aquifer. A refined three-dimensional digital hydrogeologic framework model includes the following hydrogeologic units from bottom to top: (1) the effective hydrologic basement of the basin aquifer, which includes the Pliocene Bouse Formation, Tertiary volcanic and sedimentary rocks, and pre-Tertiary metamorphic and plutonic rocks; (2) undifferentiated lower units to represent the Pliocene transition zone and wedge zone; (3) coarse-gravel unit; (4) lower, middle, and upper basin fill to represent the upper, fine-grained zone between the top of the coarse-gravel unit and the land surface; and (5) clay A and clay B. Data for the refined model includes digital elevation models, borehole lithology data, geophysical data, and structural data to represent the geometry of the hydrogeologic units. The top surface of the coarse-gravel unit, defined by using borehole and geophysical data, varies similarly to terraces resulting from the down cutting of the Colorado River. Clay A is nearly the same as the previous conceptual hydrogeologic model definition (Olmsted and others, 1973), except for a minor westward extension from the city of Yuma. Clay B is extended to the southerly international boundary and increased in areal extent by about two-thirds of the original extent (Olmsted and others, 1973). The other hydrogeologic units generally are the same as in the previous conceptual hydrogeologic model. Before development, the Colorado and Gila Rivers were the sources of nearly all the ground water in the Yuma area through direct infiltration of water from river channels and annual overbank flooding. After construction of upstream reservoirs and clearing and irrigation of the floodplains, the rivers now act as drains for the ground water. Ground-water levels in most of the Yuma area are higher now than they were in predevelopment time. A general gradient of ground-water flow toward the natural discharge area south of the Yuma area still exists, but many other changes in flow are evident. Ground water in Yuma Valley once flowed away from the Colorado River, but now has a component of flow towards the river and Mexicali Valley. A ground-water mound has formed under Yuma Mesa from long-term surface-water irrigation; about 600,000 to 800,000 acre-ft of water are stored in the mound. Ground-water withdrawals adjacent to the southerly international boundary have resulted in water-level declines in that area. The reviewed and documented water budget includes the following components: (1) recharge in irrigated areas, (2) evapotranspiration by irrigated crops and phreatophytes, (3) ground-water return flow to the Colorado River, and (4) ground-water withdrawals (including those in Mexicali Valley). Recharge components were calculated by subtracting the amount of water used by crops from the amount of water delivered. Evapotranspiration rates were calculated on the basis of established methods, thus were appropriate for input to the ground-water flow model developed by the Bureau of Reclamation (William Greer, hydrologist, Bureau of Reclamation, written commun., 2005). Evapotranspiration by crops and phreatophytes were calculated by using crop coefficient methods and meteorological data. Other methods of calculating evapotranspiration rates by using combinations of satellite imagery and ground-based data could be used for higher spatial and temporal resolution. Ground-water return flow during years of low flow on the Colorado River (1972–82, 1987–92, and 1994–96) averaged 79,000 acre-ft per year. Ground-water withdrawal data for 1970–99 were similar to other estimates made by the U.S. Geological Survey for the Yuma area. New water-chemistry data were collected in 12 wells and 8 canals/drains to characterize spatial patterns in chemical constituents, determine isotopic ages of water, infer possible sources of ground water, and locate the vertical intervals of the aquifer that contribute most water to wells. Depth-dependent samples were collected at one of the wells (YM-10). A large quantity of water-quality data were compiled from Bureau of Reclamation and U.S. Geological Survey records and merged into the U.S. Geological Survey National Water Information System database. New samples were analyzed for major ions, nutrients, stable isotopes of oxygen and hydrogen, tritium ( 3 H), and carbon-14 ( 14 C) (along with C 13 /C 12 ratios). Light values of oxygen-18 ( 18 O) and deuterium ( 2 H, D) in well 242-2 indicate recharge from the Colorado River. Heavy water samples from wells 242-22, CADC, and Mesa del Sol indicate local recharge sources. Tritium data indicate there is young water in wells in the valleys and near the edge of Yuma Mesa, while older water is found far from the Colorado River. 14 C data indicate that water from wells near the southerly international boundary is at least several thousand years old.

Arizona, California↗

USGS2006

Forecasting selenium discharges to the San Francisco Bay-Delta Estuary: Ecological effects of a proposed San Luis Drain extension

Selenium discharges to the San Francisco Bay-Delta Estuary (Bay-Delta) could change significantly if federal and state agencies (1) approve an extension of the San Luis Drain to convey agricultural drainage from the western San Joaquin Valley to the North Bay (Suisun Bay, Carquinez Strait, and San Pablo Bay); (2) allow changes in flow patterns of the lower San Joaquin River and Bay-Delta while using an existing portion of the San Luis Drain to convey agricultural drainage to a tributary of the San Joaquin River; or (3) revise selenium criteria for the protection of aquatic life or issue criteria for the protection of wildlife. Understanding the biotransfer of selenium is essential to evaluating effects of selenium on Bay-Delta ecosystems. Confusion about selenium threats to fish and wildlife stem from (1) monitoring programs that do not address specific protocols necessary for an element that bioaccumulates; and (2) failure to consider the full complexity of the processes that result in selenium toxicity. Past studies show that predators are more at risk from selenium contamination than their prey, making it difficult to use traditional methods to predict risk from environmental concentrations alone. This report presents an approach to conceptualize and model the fate and effects of selenium under various load scenarios from the San Joaquin Valley. For each potential load, progressive forecasts show resulting (1) water-column concentration; (2) speciation; (3) transformation to particulate form; (4) particulate concentration; (5) bioaccumulation by invertebrates; (6) trophic transfer to predators; and (7) effects on those predators. Enough is known to establish a first-order understanding of relevant conditions, biological response, and ecological risks should selenium be discharged directly into the North Bay through a conveyance such as a proposed extension of the San Luis Drain. The approach presented here, the Bay-Delta selenium model, determines the mass, fate, and effects of selenium released to the Bay-Delta through use of (1) historical land-use, drainage, alluvial-fill, and runoff databases; (2) existing knowledge concerning biogeochemical reactions and physiological parameters of selenium (e.g., speciation, partitioning between dissolved and particulate forms, and bivalve assimilation efficiency); and (3) site-specific data mainly from 1986 to 1996 for clams and bottom-feeding fish and birds. Selenium load scenarios consider effluents from North Bay oil refineries and discharges of agricultural drainage from the San Joaquin Valley to enable calculation of (a) a composite freshwater endmember selenium concentration at the head of the estuary; and (b) a selenium concentration at a selected seawater location (Carquinez Strait) as a foundation for modeling. Analysis of selenium effects also takes into account the mode of conveyance for agricultural drainage (i.e., the San Luis Drain or San Joaquin River); and flows of the Sacramento River and San Joaquin River on a seasonal or monthly basis. Load scenarios for San Joaquin Valley mirror predictions made since 1955 of a worsening salt (and by inference, selenium) build-up exacerbated by an arid climate and massive irrigation. The reservoir of selenium in the San Joaquin Valley is sufficient to provide loading at an annual rate of approximately 42,500 pounds of selenium to a Bay-Delta disposal point for 63 to 304 years at the lower range of projections presented here, even if influx of selenium from the California Coast Ranges could be curtailed. Disposal of wastewaters on an annual basis outside of the San Joaquin Valley may slow the degradation of valley resources, but drainage alone cannot alleviate the salt and selenium build-up in the San Joaquin Valley, at least within a century. Load scenarios also show the different proportions of selenium loading to the Bay-Delta. Oil refinery loads from 1986 to 1992 ranged from 8.5 to 20 pounds of selenium per day; with treatment and cleanup, loads decreased to 3.0 pounds of selenium per day in 1999. In contrast, San Joaquin Valley agricultural drainage loads disposed of in a San Luis Drain extension could range from 45 to 117 pounds of selenium per day across a set of historical and future conditions. Components of this valley-wide load include five source subareas (i.e., Grassland, Westlands, Tulare, Kern, and Northern) defined by water and drainage management. Loads vary per subarea mainly because of proximity of the subarea to geologic sources of selenium and irrigation history. Loads from the Sacramento River, depending on flow conditions, range from 0.8 to 10 pounds of selenium per day. Loads from the San Joaquin River vary depending on restoration and flow conditions, which are considered. A consistent picture of ecological risk emerges under modeled selenium discharges from a proposed San Luis Drain extension. The threat to the estuary is greatest during low flow seasons and critically dry years. Where selenium undergoes reactions typical of low flow or longer residence time, highly problematic bioaccumulation in prey (food) is forecast. Surf scoter, greater and lesser scaup, and white sturgeon appear to be most at risk because these Bay-Delta predators feed on deposit and filter-feeding bivalves. Recent findings add Sacramento splittail and Dungeness crab to that list. During the low flow season of critically dry years, forecasted selenium concentrations in water, particulate matter, prey (diet), and predator tissue exceed guidelines with a high certainty of producing adverse effects under the most likely load scenario from a proposed San Luis Drain extension. High flows afford some protection under certain conditions in modeled San Joaquin River scenarios. However, meeting a combined goal of releasing a specific load during maximum flows and keeping selenium concentrations in the river below a certain objective to protect against bioaccumulation may not always be attainable. Management of the San Joaquin River on a constant concentration basis also could create problematic bioaccumulation during a wet year, especially during the low flow season, because high flows translate to high loads that are not always offset by seasonal river inflows. Prior to refinery cleanup, selenium contamination was sufficient to threaten reproduction in key species within the Bay-Delta ecosystems and human health advisories were posted based on selenium concentrations in tissues of diving ducks. During this time, selenium concentrations in the Bay-Delta were well below the most stringent recommended water quality criterion [1 microgram per liter (1 µg/L)]. Enhanced biogeochemical transformations to bioavailable particulate selenium and efficient bioaccumulation by bivalves characterized the system. If these biogeochemical conditions continue to prevail and agricultural selenium sources replace or exceed refinery sources, ecological forecasts suggest the risk of adverse effects will be difficult to eliminate under an out-of-valley resolution to the selenium problem. The Bay-Delta selenium model presented here is a systematic approach for conducting forecasts of the ecological effects from selenium on aquatic food webs. It is a new tool that links and models the major processes leading from loads through consumer organisms to predators. It also is a feasible approach for site-specific analysis and could provide a framework for developing new protective selenium foodweb guidelines and predator criteria. Model components that help ensure understanding ecosystems and the basis of environmental protection are (1) contaminant concentrations and speciation in sources, such as particulate material, that most influence bioavailability; (2) bioaccumulation models that calculate concentrations in diet, specifically in bivalves of the Bay-Delta that act as sensitive indicators of selenium contamination; (3) food-web type that determines what animals are threatened and when; and (4) multiple media concentrations (water, particulate material, and tissue of prey and predators) that, in-combination, determine risk or hazard.

California↗

USGS2009

Simulation of the Groundwater-Flow System in Pierce, Polk, and St. Croix Counties, Wisconsin

Groundwater is the sole source of residential water supply in Pierce, Polk, and St. Croix Counties, Wisconsin. A regional three-dimensional groundwater-flow model and three associated demonstration inset models were developed to simulate the groundwater-flow systems in the three-county area. The models were developed by the U.S. Geological Survey in cooperation with the three county governments. The objectives of the regional model of Pierce, Polk, and St. Croix Counties were to improve understanding of the groundwaterflow system and to develop a tool suitable for evaluating the effects of potential water-management programs. The regional groundwater-flow model described in this report simulates the major hydrogeologic features of the modeled area, including bedrock and surficial aquifers, groundwater/surface-water interactions, and groundwater withdrawals from high-capacity wells. Results from the regional model indicate that about 82 percent of groundwater in the three counties is from recharge within the counties; 15 percent is from surface-water sources, consisting primarily of recirculated groundwater seepage in areas with abrupt surface-water-level changes, such as near waterfalls, dams, and the downgradient side of reservoirs and lakes; and 4 percent is from inflow across the county boundaries. Groundwater flow out of the counties is to streams (85 percent), outflow across county boundaries (14 percent), and pumping wells (1 percent). These results demonstrate that the primary source of groundwater withdrawn by pumping wells is water that recharges within the counties and would otherwise discharge to local streams and lakes. Under current conditions, the St. Croix and Mississippi Rivers are groundwater discharge locations (gaining reaches) and appear to function as 'fully penetrating' hydraulic boundaries such that groundwater does not cross between Wisconsin and Minnesota beneath them. Being hydraulic boundaries, however, they can change in response to water withdrawals. Tributary rivers act as 'partially penetrating' hydraulic boundaries such that groundwater can flow underneath them through the deep sandstone aquifers. The model also demonstrates the effects of development on groundwater in the study area. Water-level declines since predevelopment (no withdrawal wells) are most pronounced where pumping is greatest and flow between layered aquifers is impeded by confining units or faults. The maximum simulated water-level decline is about 40 feet in the deep Mount Simon aquifer below the city of Hudson, Wisconsin. Three inset models were extracted from the regional model to demonstrate the process and additional capabilities of the U.S. Geological Survey MODFLOW code. Although the inset models were designed to provide information about the groundwater-flow system, results from the inset models are presented for demonstration purposes only and are not sufficiently detailed or calibrated to be used for decisionmaking purposes without refinement. Simulation of groundwater/lake-water interaction around Twin Lakes near Roberts, in St. Croix County, Wisconsin, showed that groundwater represents approximately 5 to 20 percent of the overall lake-water budget. Groundwater-contributing areas to streams in western Pierce County are generally similar in size to the surface-water-contributing areas but do not necessarily correspond to the same land area. Transient streamflow simulations of Osceola Creek in Polk County demonstrate how stream base flow can be influenced not only by seasonal precipitation and recharge variability but also by systematic changes to the system, such as groundwater withdrawal from wells.

Scientific Investigations Report↗

USGS2020

Spatially referenced models of streamflow and nitrogen, phosphorus, and suspended-sediment loads in the southeastern United States

Spatially Referenced Regression On Watershed attributes (SPARROW) models were applied to describe and estimate mean-annual streamflow and transport of total nitrogen (TN), total phosphorus (TP), and suspended sediment (SS) in streams and delivered to coastal waters of the southeastern United States on the basis of inputs and management practices centered near 2012, the base year of the model. Previously published TN and TP models for 2002 served as a starting point and reference for comparison. The datasets developed for the 2012 models not only represent updates of previous conditions but also incorporate new approaches for characterizing sources and transport processes that were not available for previous models. Variability in streamflow across the southeastern United States was explained as a function of precipitation adjusted for evapotranspiration, spring discharge, and municipal and domestic wastewater discharges to streams. Results from the streamflow model were used as input to the water-quality SPARROW models, and areas with large streamflow prediction errors—urban areas and karst areas—were used to provide guidance on where additional data are needed to improve routing of flow. Variability in TN transport in Southeast streams was explained by the following five sources in order of decreasing mass contribution to streams: atmospheric deposition, agricultural fertilizer, municipal wastewater, manure from livestock, and urban land. Variable rates of TN delivery from source to stream were attributed to variation among catchments in climate, soil texture, and vegetative cover, including the extent of cover crops in the watershed. Variability in TP transport in Southeast streams was explained by the following six sources in order of decreasing mass contribution to streams: parent-rock minerals, urban land, manure from livestock, municipal wastewater, agricultural fertilizer, and phosphate mining. Varying rates of TP delivery were attributed to variation in climate, soil erodibility, depth to water table, and the extent of conservation tillage practices in the watershed. Variability in SS transport in Southeast streams was explained by variable sediment export rates for different combinations of land cover and geologic setting (for upland sources of sediment) and by gains in stream power caused by longitudinal changes in channel hydraulics (for channel sources of sediment). Sediment yields for the transitional land cover (shrub, scrub, herbaceous, and barren) varied widely depending on geologic setting and on agricultural land cover. Varying rates of SS delivery, like those for TP, were attributed to variation in climate, soil erodibility, and the extent of conservation tillage practices in the watershed, as well as to areal extent of canopy land cover in the 100-meter buffer along the channel. Relatively large uncertainty, compared to the other three models, for almost all the SS source coefficients indicates the need for caution when interpreting the results from the sediment model. TN, TP, and SS inputs to streams from sources were balanced in the models with losses from physical processes in streams and reservoirs and with water withdrawals. The losses in streams and reservoirs along with withdrawals removed 35, 44, and 65 percent of the TN, TP, and SS load, respectively, that entered streams before reaching coastal waters.

Scientific Investigations Report↗

USGS2011

Sources, fate, and transport of nitrogen and phosphorus in the Chesapeake Bay watershed: An empirical model

Spatially Referenced Regression on Watershed Attributes (SPARROW) was used to provide empirical estimates of the sources, fate, and transport of total nitrogen (TN) and total phosphorus (TP) in the Chesapeake Bay watershed, and the mean annual TN and TP flux to the bay and in each of 80,579 nontidal tributary stream reaches. Restoration efforts in recent decades have been insufficient to meet established standards for water quality and ecological conditions in Chesapeake Bay. The bay watershed includes 166,000 square kilometers of mixed land uses, multiple nutrient sources, and variable hydrogeologic, soil, and weather conditions, and bay restoration is complicated by the multitude of nutrient sources and complex interacting factors affecting the occurrence, fate, and transport of nitrogen and phosphorus from source areas to streams and the estuary. Effective and efficient nutrient management at the regional scale in support of Chesapeake Bay restoration requires a comprehensive understanding of the sources, fate, and transport of nitrogen and phosphorus in the watershed, which is only available through regional models. The current models, Chesapeake Bay nutrient SPARROW models, version 4 (CBTN_v4 and CBTP_v4), were constructed at a finer spatial resolution than previous SPARROW models for the Chesapeake Bay watershed (versions 1, 2, and 3), and include an updated timeframe and modified sources and other explantory terms. Chesapeake Bay receives an estimated 1.32 x 10 8 kilograms (132,000 metric tons) of nitrogen and 9.74 x 10 6 kilograms (9,740 metric tons) of phosphorus annually from its watershed, mainly through its two largest tributaries, the Susquehanna and Potomac Rivers. Significant (α=0.10) sources of nutrients to streams in the watershed include fertilizer and manure applications in agricultural areas, undifferentiated urban activities, point sources, atmospheric deposition and direct fixation by crops (for nitrogen), and mineral sources (for phosphorus). Agriculture (primarily fertilizer applications and crop fixation) contributes more than half of the nitrogen delivered from the watershed to the bay; phosphorus contributions are more mixed, and fairly evenly distributed among agricultural (fertilizer and manure applications) and urban (including point) sources. Natural mineral dissolution contributes approximately 14 percent of the phosphorus flux from the watershed to the Chesapeake Bay. Empirical estimates of average yields from different source areas and of the portion of selected applications delivered to streams agree closely with previously reported values in the literature. Nutrient fate and transport through the Chesapeake Bay watershed to the bay reflect the diferent physical and chemical properties of nitrogen and phosphorus compounds. Groundwater is an important pathway for nitrogen transport (as nitrate), and TN flux is greatest in areas with greater groundwater flow and in areas of the Piedmont underlain by carbonate rocks. TN flux decreases with increasing vegetative growth (likely indicative of plant uptake) and soil available water capacity (likely indicative of reducing conditions). Phosphorus transport to streams, conversely, is greatest in areas most likely to generate overland runoff and related erosion, including those with less permeable and more erodible soils and greater precipitation. Phosphorus transport also is greater in the Coastal Plain than in other areas, possibly due to saturation of soils with historical phosphorus applications. Both nitrogen and phosphorus are lost within watershed impoundments (lakes, ponds, or reservoirs), and nitrogen is also lost significantly along flowing reaches, particularly in small streams and in larger streams in warmer areas.

Chesapeake Bay↗

USGS2022

Sources and characteristics of dissolved organic carbon in the McKenzie River, Oregon, related to the formation of disinfection by-products in treated drinking water

Executive Summary This study characterized the concentration and quality of dissolved organic carbon (DOC) in the McKenzie River, a relatively undeveloped watershed in western Oregon, and its link to forming disinfection by-products (DBPs) in treated drinking water. The study aimed to identify the primary source(s) of DOC in source water for the Eugene Water & Electric Board’s (EWEB) conventional treatment plant on the McKenzie River near river mile 11, upstream of Hayden Bridge. The two classes of regulated compounds examined—trihalomethanes (THMs) and haloacetic acids (HAAs)—form when organic carbon in raw source water reacts with chlorine and (or) bromine during water treatment. The objectives of the study were to: characterize the amount and quality of DOC in the McKenzie River and select tributaries during storms; identify the most common types of carbon using UV-vis spectroscopy and other methods; evaluate optical properties for predicting DBP precursors in surface water; and identify land cover classes or vegetation types that may be important sources of organic carbon and DBP precursors in EWEB’s source water. Eleven storms were sampled synoptically in upstream-to-downstream fashion to provide a “snapshot” of water quality conditions at four sites on the McKenzie River from Frissell Bridge (6 miles downstream from Trail Bridge Reservoir) to the EWEB water treatment plant at Hayden Bridge and nine contributing tributaries. Storms included late summer and early autumn “first flush” events and late autumn, winter, and spring storms spanning a range in streamflows from 3,000 to 26,000 cubic feet per second as measured in the main stem McKenzie River at the EWEB water intake. Water samples were analyzed for DOC concentrations and optical properties (fluorescence and ultraviolet absorbance [UVA]) across a range of wavelengths to characterize the quantity and quality of dissolved organic matter (DOM) in the McKenzie River at the drinking water intake and upstream locations. Paired sets of source and finished water samples were collected at the EWEB treatment plant to identify DOC quality parameters in raw source water that might predict DBP concentrations in finished drinking water. DOC concentrations were relatively low in the McKenzie River (0.4–3 milligrams per liter [mg/L]; average 1.5 mg/L) but much higher in the tributaries. The highest DOC concentrations occurred during “first flush” storms in October 2012 and September 2013; the highest value (16 mg/L) was measured at the 52nd Street stormwater outfall. The average DOC concentration in the lower basin-tributaries was 3.8 mg/L; three middle basin tributaries—Quartz, Gate, and Haagen Creeks, which drain private forestland with less coniferous forest compared with other higher elevation tributaries— had slightly lower average DOC concentrations (2.8 mg/L). These middle-basin watersheds may be important sources of DOC and DBP precursors to the McKenzie River, even more so than the lower basin tributaries, depending on their flows (and loads). This is particularly true after the September 2020 Holiday Farm fire, which burned much of this area. DOC concentrations increased 68 percent in the McKenzie River between the uppermost reference site at Frissell Bridge and Vida; this includes drainage from Quartz Creek, Blue River Lake and Cougar Reservoir, which all contributed DOC to the main stem. In contrast, the lowermost tributaries draining most of the agricultural and urban land did not have a large effect on DOC in the McKenzie River despite their higher DOC concentrations because of their presumed relatively low streamflows and, consequently, DOC loads. Apart from the continuous flow monitors in the McKenzie River and some tributaries (Blue River and South Fork McKenzie River, and streamflow at Hayden Bridge and Vida, Camp Creek and some other locations), streamflow was not assessed during sample collection for this study. This lack of streamflow data precludes a detailed analysis of loads, which is discussed in the future studies section. All DBP concentrations in finished drinking water were less than EPA maximum contaminant levels (MCLs) of 0.080 mg/L for the four trihalomethanes (THM4) and 0.060 mg/L for five haloacetic acids (HAA5). During the 11 storm sampling events the maximum summed concentrations were about 0.040 mg/L for both THM4 and HAA5. Compliance monitoring samples, collected separately by EWEB, yielded some higher concentrations—0.046 mg/L THM4 and 0.047 HAA5—during the December 2012 storm. The corresponding benchmark quotient (BQ) values, which indicate how close a measured DBP concentration is to the MCL, were 0.58 and 0.78, respectively, for THM4 and HAA5. Compared with a similar 2007–08 McKenzie River study that did not target storm events, concentrations of THM4 and HAA5 in finished water were 68 percent and 33 percent higher, respectively, during the current study. Due to the high dilution rates in the McKenzie River main stem, many of the individual fluorescence excitation-emission measurements were low (<0.1 Raman units) and approached analytical detection limits. Parallel factor analysis (PARAFAC) resulted in a five-component model (C1–C5) that represents five unique organic fluorophores. Components C1, C2, and C3 represent DOM associated with soil-derived, humic-like, more degraded organic matter. In contrast, components C4 and C5 represent “fresher” DOM, derived from terrestrial and aquatic plants, including algae and cyanobacteria that are common in the McKenzie River and its tributaries and reservoirs. The fluorescence data and PARAFAC modeling suggest that most of the DOC in the McKenzie River originated from terrestrial sources (primarily components C1 and C2). The largest increases in DOC in the main stem occurred in the reach upstream of Vida, from inflows by Quartz Creek, Blue River, South Fork McKenzie River, and other tributaries. Concentrations of DBPs in EWEB’s finished drinking water were positively correlated with DOC concentrations in raw source water (THM4, p <0.05; HAA5, p <0.01) for paired samples collected 12−24 hours apart. DOC concentrations were significantly positively correlated ( p <0.001) with laboratory-based fluorescent dissolved organic matter (fDOM) measurements, suggesting fDOM as a useful parameter for monitoring and predicting DOC concentration in surface water and DBP concentrations in finished water. Of all the PARAFAC components in surface water, C5 had the highest correlations with DBPs in finished water (rho = 0.77–0.84, p <0.01), followed by components C1 and C2 (rho = 0.75 and 0.71, respectively, p <0.01). This C5 carbon is associated with recently produced DOM, possibly from decomposed terrestrial and aquatic vegetation. Model loadings of these three components were considerably higher in the sampled tributaries relative to the main stem McKenzie River, with most of the observed increases in the main stem apparent at Vida. This points to Quartz Creek or other tributaries in the reach between Frissell Bridge and the sampling site near Vida (South Fork McKenzie and Blue Rivers) as potentially key contributors of DOM source material that leads to the production of DBPs in treated drinking water. A limited load analysis showed that the reservoirs contributed 8–37 percent of the instantaneous DOC loads observed at Vida at the time of sampling, which suggests other sources such as Quartz Creek and other streams in the reach between Frissell Bridge and Vida are more important. Random forest analyses identified PARAFAC components C1 and C5 and fluorescence peaks A, C, M, T and N as the best predictors for HAA5 concentrations in finished drinking water, explaining 62.5 percent of the variation. The best predictors for THM4 were C1, C4 + C5, and peaks T, A, and N, which explained 33 percent of the variation. Several land cover and vegetation classes were correlated with DOC concentration and other optical measurements. The percentage of evergreen forest in each of the subwatersheds sampled was negatively correlated ( p <0.001) with DOC concentration and many optical indicators of DOM quantity: UVA 254 , fDOM, and all of the fluorescence peaks. In contrast, mixed (deciduous) forest was positively correlated ( p <0.001) with DOC, fDOM, UVA 254 , and several fluorescence peaks, demonstrating the importance of deciduous leaf fall in generating DOC and DBP precursors. The high level of human activities in the middle and lower portion of the basin—including timber harvesting and road construction on private forestland, agricultural, rural, industrial, and urban development—have resulted in the greatest loss in native coniferous and mixed deciduous forests in the basin. DOC loading from these tributaries and reservoir releases, which contain DOC from terrestrial and aquatic productivity, both enrich the McKenzie River. Concentrations of DOC increased an average of 71 percent (range 30–120 percent) in the McKenzie River between Frissell Bridge, the upstream reference site, and Vida. PARAFAC components C1, C2, and C5—which were correlated with DBPs in finished water—increased, on average, 109–136 percent (range 20–250 percent) in this same Frissell-to-Vida reach. These increases occur from input of tributaries in the middle basin such as Quartz Creek and others, as noted above. Future monitoring, field, and lab studies can improve our understanding of seasonal and spatial sources of organic carbon contributing DBP precursors to the McKenzie River and allow detection of long-term trends resulting from the recent Holiday Farm Fire, which burned 173,393 acres of forestland, including riparian areas along the main stem, and numerous structures, homes, and outbuildings in September 2020. Future studies could examine DOC fluxes and flushing of carbon from the watershed, investigate the role of precipitation amount and intensity in mobilizing carbon and sediment, and evaluate impacts to aquatic communities and human health as part of a post-fire assessment. Other areas ripe for study include evaluating the impacts of potential temperature increases on carbon sequestration and decomposition in the burned and unburned forests and identifying practices that foster sequestration of carbon in forest soils. The use of fluorescence sensors such as fDOM to monitor the concentration and composition of raw water supplies may be improved for detection of specific DBP precursors, to provide continuous and real-time information to treatment plant operators. Future studies that monitor DOM amount and quality, and DBP Formation Potential (FP), particularly during storm events, paired with streamflow measurements, as suggested above, could help identify areas that contribute high DOC loads and thus help managers identify the key areas to focus restoration activities. Other studies could examine treatment options for currently regulated DBPs and potentially unregulated compounds, including advanced biological treatments for their removal. This study was a collaboration between the U.S. Geological Survey (USGS) and EWEB in Eugene, Oregon, with additional funding provided from USGS Cooperative Matching Funds Program.

Oregon↗

USGS2010

Geohydrology of the stratified-drift aquifer system in the lower Sixmile Creek and Willseyville Creek trough, Tompkins County, New York

In 2002, the U.S. Geological Survey, in cooperation with the Tompkins County Planning Department began a series of studies of the stratified-drift aquifers in Tompkins County to provide geohydrologic data for planners to develop a strategy to manage and protect their water resources. This aquifer study in lower Sixmile Creek and Willseyville Creek trough is the second in a series of aquifer studies in Tompkins County. The study area is within the northern area of the Appalachian Plateau and extends about 9 miles from the boundary between Tompkins County and Tioga County in the south to just south of the City of Ithaca in the north. In lower Sixmile Creek and Willseyville Creek trough, confined sand and gravel aquifers comprise the major water-bearing units while less extensive unconfined units form minor aquifers. About 600 people who live in lower Sixmile Creek and Willseyville Creek trough rely on groundwater from the stratified-drift aquifer system. In addition, water is used by non-permanent residents such as staff at commercial facilities. The estimated total groundwater withdrawn for domestic use is about 45,000 gallons per day (gal/d) or 0.07 cubic foot per second (ft3/s) based on an average water use of 75 gal/d per person for self-supplied water systems in New York. Scouring of bedrock in the preglacial lower Sixmile Creek and Willseyville Creek valleys by glaciers and subglacial meltwaters truncated hillside spurs, formed U-shaped, transverse valley profiles, smoothed valley walls, and deepened the valleys by as much as 300 feet (ft), forming a continuous trough. The unconsolidated deposits in the study area consist mostly of glacial drift, both unstratified drift (till) and stratified drift (laminated lake, deltaic, and glaciofluvial sediments), as well as some post-glacial stratified sediments (lake-bottom sediments that were deposited in reservoirs, peat and muck that were deposited in wetlands, and alluvium deposited by streams). Multiple advances and retreats of the ice in the study area resulted in several sequences of various types of glacial deposits. A large moraine (Valley Heads Moraine) dominates the southern part of the study area, a large delta dominates the central part, and ground moraine (mostly till) dominates the northern part. Glacial sediments in the center of the lower Sixmile Creek and Willseyville Creek trough typically range from 150 to 200 ft but can be greater than 300 ft in some places. Where the sediments are composed of sand and gravel they form aquifers. In most parts of the lower Sixmile Creek and Willseyville Creek trough, there is an upper and a basal confined aquifer. However, underlying the central parts of the Brooktondale delta, there are as many as four confined aquifers, whereas in the northern part of the study area, only one extensive confined aquifer is present. The major sources of recharge to these confined aquifers are (1) direct infiltration of precipitation where confined aquifers crop out at land surface (mostly along the western trough wall in the southern and central parts of the study area and, to a lesser degree, along the eastern trough wall); (2) unchanneled surface and subsurface runoff from adjacent upland areas that seeps into the aquifer along the western trough walls; (3) subsurface flow from underlying till or bedrock at the lateral contacts at trough walls; (4) adjacent fine-grained stratified drift, especially when the aquifer is pumped; and (5) discharge from bedrock at the bottom and sides of the trough. In the central part of the study area, the surficial coarse-grained sediments (sand and gravel) comprise a delta near Brooktondale and form a small unconfined aquifer (0.3 square mile). Although much of the upper part of the delta has been removed by several aggregate mining operations, sufficient amounts of sand and gravel remain in most places to form a thin unconfined aquifer. The major sources of recharge to the unconfined aquifer are (1)

Scientific Investigations Report↗

USGS2018

Science support for evaluating natural recovery of polychlorinated biphenyl concentrations in fish from Crab Orchard Lake, Crab Orchard National Wildlife Refuge, Illinois

Introduction Crab Orchard Lake in southern Illinois is one of the largest and most popular recreational lakes in the state. Construction of the nearly 7,000-acre reservoir in the late 1930s created employment opportunities through the Works Progress Administration, and the lake itself was intended to supply water, control flooding, and provide recreational opportunities for local communities (Stall, 1954). In 1942, the Department of War appropriated or purchased more than 20,000 acres of land around Crab Orchard Lake and constructed the Illinois Ordnance Plant, which manufactured bombs and anti-tank mines during World War II. After the war, an Act of Congress transferred the property to the U.S. Department of the Interior. Crab Orchard National Wildlife Refuge was established on August 5, 1947, for the joint purposes of wildlife conservation, agriculture, recreation, and industry. Production of explosives continued, but new industries also moved onsite. More than 200 tenants have held leases with Crab Orchard National Wildlife Refuge and have operated a variety of manufacturing plants (electrical components, plated metal parts, ink, machined parts, painted products, and boats) on-site. Soils, water, and sediments in several areas of the refuge were contaminated with hazardous substances from handling and disposal methods that are no longer acceptable environmental practice (for example, direct discharge to surface water, use of unlined landfills). Polychlorinated biphenyl (PCB) contamination at the refuge was identified in the 1970s, and a PCB-based fish-consumption advisory has been in effect since 1988 for Crab Orchard Lake. The present advisory covers common carp ( Cyprinus carpio ) and channel catfish ( Ictalurus punctatus ); see Illinois Department of Public Health (2017). Some of the most contaminated areas of the refuge were actively remediated, and natural ecosystem recovery processes are expected to further reduce residual PCB concentrations in the lake. The U.S. Fish and Wildlife Service sought technical support to understand environmental drivers of current (2017) PCB residues in fish tissue and patterns in PCB residues through time to inform the fish-consumption advisory for Crab Orchard Lake. This project is planned in two phases (Tasks 1 and 2); the first phase is included in this report. Task 1, reported here, includes a review of existing literature and a brief overview focused on environmental and biochemical/physiological processes that drive PCB residues in fish tissue. This review specifically targets processes that are relevant for freshwater lacustrine environments such as those at Crab Orchard Lake. In addition to discussions of environmental fate, metabolism, and accumulation of PCBs, this review includes a brief scientifically based explanation of approaches used to establish fish-consumption advisories. A planned second task (Task 2) will include a compilation and summary of existing data on PCB residues in fish tissue samples from Crab Orchard Lake. This summary will also place Crab Orchard Lake data in a broader geographic context through a comparison with fish data from other Midwestern lakes. When Task 1 and Task 2 are complete, resource managers will have (a) a synthesis of existing literature that characterizes the processes influencing the fate of residual PCBs remaining in systems such as Crab Orchard Lake, (b) a summary of natural PCB attenuation processes for use in risk communication with the public, and (c) a summary of existing data on PCBs in fish tissues from Crab Orchard Lake, including exploratory plots of tissue residues through time. Overall, this project will provide data to help resource managers better understand the ecological and public health consequences of residual PCBs in Crab Orchard Lake.

Illinois↗

USGS1979

Summary appraisals of the nation's ground-water resources – Lower Mississippi region

The Lower Mississippi Region comprises an area of 102,400 square miles (265,200 square kilometers). Almost all this area is in the physiographic province known as the Gulf Coastal Plain. Three small areas on the northwest boundary of the region are in the Interior Highlands. The Lower Mississippi Region has an abundance of ground water. The geologic structure in that part of the region within the Coastal Plain is an elongated trough which has been filled with permeable materials, resulting in vast subsurface reservoirs. Except in local areas where continued large withdrawals have caused significant water-level declines, these reservoirs are full. Recharge to the region's aquifers is primarily from rainfall. Annual rainfall in most of the region is well distributed throughout the year and is sufficient to satisfy evapotranspiration requirements and still provide recharge to the aquifers. An estimated 844 billion cubic feet (24 billion cubic meters) of fresh ground water is available for withdrawal annually in the region. Only about one-third of this quantity is being utilized. Therefore, on this basis alone, the region still has much potential for ground-water development. The Coastal Plain aquifers within the Lower Mississippi Region contain large reserves of saltwater in the downdip limits of the aquifers. The quantity of saltwater in the region is several times that of freshwater. As desalinization techniques are developed and as more uses are found for saltwater, this reserve could become an important source of water for the region. At present (1976), the most productive and potentially productive aquifers or aquifer systems in the region are the Mississippi River valley alluvial aquifer of Quaternary age and the Sparta Sand and the Memphis aquifer (Memphis Sand in Tennessee) of Tertiary age. The Sparta Sand and the Memphis aquifer are heavily utilized and have shown significant water-level declines. However, selected well hydrographs indicate that water levels may be stabilizing under present pumping conditions. The Mississippi River valley alluvial aquifer is the most extensive high-yielding aquifer in the region; yields of several thousand gallons per minute may be obtained at depths of less than 200 feet (61 meters). To obtain maximum benefit from the vast quantities of ground water in the region, adequate attention must be given to the effects of proposed development upon the ground-water regime. Knowledge of the geologic structure and hydraulic properties of the aquifer systems is essential to an evaluation of the effects of such development. Some studies have been made in sufficient detail to provide this knowledge, but additional studies are needed. Activities that could cause significant changes in the groundwater regime should be undertaken only after all available information has been considered. Failure to seek out and use such information may result in inefficient development of the groundwater resource and, in some instances, degradation of the quality of the resource. Some changes always result from ground-water development. The possible changes can be grouped into three categories: hydraulic, water quality, and those affecting the physical framework of the aquifers. Generally, they are small in magnitude and areal extent. Because these changes occur below the ground surface, they are unknown to the ground-water user unless they noticeably affect the quantity or quality of water produced or cause obvious physical effects, such as land subsidence. Great advances have been made in hydrologic technology in recent years. Predictive models have been developed that make it possible for the hydroiogist to simulate aquifer responses to proposed development or other stresses. These models would be invaluable tools in progressive water-resources planning and management.

Arkansas, Kentucky, Louisiana, Mississippi, Missou↗

USGS2012

A brief history and summary of the effects of river engineering and dams on the Mississippi River system and delta

The U.S. Geological Survey Forecast Mekong project is providing technical assistance and information to aid management decisions and build science capacity of institutions in the Mekong River Basin. A component of this effort is to produce a synthesis of the effects of dams and other engineering structures on large-river hydrology, sediment transport, geomorphology, ecology, water quality, and deltaic systems. The Mississippi River Basin (MRB) of the United States was used as the backdrop and context for this synthesis because it is a continental scale river system with a total annual water discharge proportional to the Mekong River, has been highly engineered over the past two centuries, and the effects of engineering have been widely studied and documented by scientists and engineers. The MRB is controlled and regulated by dams and river-engineering structures. These modifications have resulted in multiple benefits including navigation, flood control, hydropower, bank stabilization, and recreation. Dams and other river-engineering structures in the MRB have afforded the United States substantial socioeconomic benefits; however, these benefits also have transformed the hydrologic, sediment transport, geomorphic, water-quality, and ecologic characteristics of the river and its delta. Large dams on the middle Missouri River have substantially reduced the magnitude of peak floods, increased base discharges, and reduced the overall variability of intraannual discharges. The extensive system of levees and wing dikes throughout the MRB, although providing protection from intermediate magnitude floods, have reduced overall channel capacity and increased flood stage by up to 4 meters for higher magnitude floods. Prior to major river engineering, the estimated average annual sediment yield of the Mississippi River Basin was approximately 400 million metric tons. The construction of large main-channel reservoirs on the Missouri and Arkansas Rivers, sedimentation in dike fields, and protection of channel banks by revetments throughout the basin, have reduced the overall sediment yield of the MRB by more than 60 percent. The primary alterations to channel morphology by dams and other engineering projects have been (1) channel simplification and reduced dynamism; (2) lowering of channel-bed elevation; and (3) disconnection of the river channel from the flood plain, except during extreme flood events. Freshwater discharge from the Mississippi River and its associated sediment and nutrient loads strongly influence the physical and biological components in the northern Gulf of Mexico. Ninety percent of the nitrogen load reaching the Gulf of Mexico is from nonpoint sources with about 60 percent coming from fertilizer and mineralized soil nitrogen. Much of the phosphorus is from animal manure from pasture and rangelands followed by fertilizer applied to corn and soybeans. Increased nutrient enrichment in the northern Gulf of Mexico has resulted in the degradation of water quality as more phytoplankton grow, which increases turbidity and depletes oxygen in the lower depths creating what is known as the "dead zone." In 2002, the dead zone was 22,000 square kilometers (km2), an area similar to the size of the State of Massachusetts. Changes in the flow regime from engineered structures have had direct and indirect effects on the fish communities. The navigation pools in the upper Mississippi River have aged, and these overwintering habitats, which were created when the pools filled, have declined as sedimentation reduces water depth. Reproduction of paddlefish may have been adversely affected by dams, which impede access to suitable spawning habitats. Fishes that inhabit swift-current habitats in the unimpounded lower Mississippi River have not declined as much as in the upper Mississippi River. The decline of the pallid sturgeon may be attributable to channelization of the Missouri River above St. Louis, Missouri. The Missouri River supports a rich fish community and remains relatively intact. Nevertheless, the widespread and long history of human intervention in river discharge has contributed to the declines of about 25 percent of the species. The Mississippi River Delta Plain is built from six delta complexes composed of a massive area of coastal wetlands that support the largest commercial fishery in the conterminous United States. Since the early 20th century, approximately 4,900 km2 of coastal lands have been lost in Louisiana. One of the primary mechanisms of wetland loss on the Plaquemines-Balize complex is believed to be the disconnection of the river distributary network from the delta plain by the massive system of levees on the delta top, which prevent overbank flooding and replenishment of the delta top by sediment and nutrient deliveries. Efforts by Federal and State agencies to conserve and restore the Mississippi River Delta Plain began over three decades ago and have accelerated over the past decade. Regardless of these efforts, however, land losses are expected to continue because the reduced upstream sediment supplies are not sufficient to keep up with the projected depositional space being created by the combined forces of delta plain subsidence and global sea-level rise.

Alabama;Alberta;Arkansas;Colorado;Georgia;Illinois↗

USGS1967

Swatara Creek basin of southeastern Pennsylvania: An evaluation of its hydrologic system

Local concentrations of population in the Swatara Creek basin of Pennsylvania find it necessary to store, transport, and treat water because local supplies are either deficient or have been contaminated by disposal of wastes in upstream areas. Water in the basin is available for the deficient areas and for dilution of the coal-mine drainage in the northern parts and the sewage wastes in the southern parts. Swatara Creek drains 576 square miles just east of Harrisburg, Pa., and is the largest tributary to the Susquehanna River from the north side below Harrisburg. It rises in the southern Pocono Mountains and flows southwestward across the Lebanon Plateau. On an average day Swatara Creek discharges more than 630 million gallons into the Susquehanna River at Middletown, Pa. In a year this amounts to about 23 inches of water over the entire basin and is the residual from an average annual precipitation of 45.5 inches. During an average year the flow in Swatara Creek from the upper third of the basin above Harper Tavern is always greater than 1,300 mgd (million gallons per day) for at least 15 days and is always greater than 25 mgd for at least 350 days. The daily streamflow from the basin averages 1.1 mgd per sq mi, but yields from different areas range from 0.97 to 1.22 mgd per sq mi. These variations are caused chiefly by differences in precipitation and land cover. The area of lowest yield is in the valleys west of Tremont, and the highest yields are in the Upper and Lower Little Swatara Creek subbasins. At high and medium stages the chemical character of the water in the streams is suitable for public and private supplies. At lower stages, defending on the areas and the amounts of contamination by coal-mine drainage and sewage pollution, the natural flow may require some treatment. At low stages the chemical characteristics of the natural flow not affected by man is almost identical with that of the ground water in the area drained by the stream. In general, the total dissolved solids range from about 25 to 400 parts per million and the hardness is as much as about 300 parts per million. The ground-water increment to the base flow of Swatara Creek averages about 240 mgd, or about 8.8 inches annually, for the basin. Generally, ground-water supplies in amounts of less than 0.5 mgd can be developed south of Blue Mountain. Supplies of several million gallons per day have been developed for industrial use from the permeable limestones in the south-central part of the basin. More intensive investigation in other parts of the basin would indicate areas where supplies of more than 0.5 mgd could be developed from properly spaced wells. The chemical character of water from wells depends largely on the host rock. In highly soluble rocks water contains large amount of dissolved solids; in more resistant rocks concentrations are lower. The chemical character of unpolluted ground water generally reflects the composition of the more readily soluble minerals in the local geologic environment. Areas contaminated by septic- tank effluent may have above normal amounts of nitrate and detergent products. Except where polluted, most ground water is suitable for public and industrial uses without extensive treatment. Sites for storage of surface water exist in the part of the basin lying in the valley and ridge area. As much as 30 to 40 percent of the annual flow could be impounded for release as low-flow augmentation for dilution of mine drainage and other wastes in the basin. Low sediment yields of supplying drainage areas would ensure a long life expectancy of reservoirs at these sites. Overbank flooding of the main stem of the Swatara Creek and its tributaries has occurred many times in the past. However, it has not been a hazard because urban development has not encroached on the flood plain. An inundation map of the August 1933 flood provides a basis that urban planners may use to avoid future damage. As water in the Swatara Creek moves downstream to the Susquehanna River, the flow is influenced consecutively by a large annual rainfall on the northern valley and ridge area, the wastes of surface and subsurface coal-mining activities, and less annual rainfall on the part of the basin lying in the Lebanon Plateau area; the flow is supplemented and further influenced by many tributaries and by the industrial and domestic wastes that are carried by these secondary streams. The annual precipitation ranges from 52 inches at the east edge and 49 inches at the west edge of the mountainous part of the basin to about 41 inches at the southwestern part at Middletown. The rainfall generally is adequate during the growing season to mature the crops. The mean annual temperature at Lebanon is about 52°F, and the growing season is about 180 days. In this report the basin has been divided into eight hydrologic zones, leased on runoff, natural use of water, and chemical character of water. Four zones lie in the valley and ridge area, three lie in the Lebanon Plateau area, and one lies in the highland along the southeastern basin boundary. In each of the zones the hydrologic characteristics are virtually the same, but they may be completely different from those in adjacent zones. The boundaries of the zones generally coincide with boundaries between geologic formations, and the areas in each zone include rocks of similar influence on water. Streams in zone 4 at the northeast edge of the plateau have the highest average surface runoff from 1.2 to 1.1 mgd per sq mi whereas those in zone 2 at the northwest edge of the valley and ridge area have the lowest, about 1.0 mgd. Streams in zone 8, along the southeast edge of the basin, have the largest sustained low-flow yield, about 0.26 to 0.19 mgd per sq mi; those in zone 5 overlying the Martinsburg Shale east of Harrisburg have the smallest sustained low-flow yields, 0.03 to 0.01 mgd. Streams in the limestone area of zone 7 have the greatest range in low-flow yields in any one zone from 0.60 to 0 mgd per sq mi. Low-flow yields in zones 1 through 4 range from 0.13 to 0.03 mgd per sq mi. Surface flows from zones 1 and 2 are generally acidic and contain high concentrations of sulfate, iron, and total dissolved solids especially where contaminated with mine wastes. Surface flows from zones 3 and 4 are dilute, slightly alkaline, and suitable for public water supplies. Surface flows from zones 5, 6, and 7 are alkaline and contain moderate concentrations of dissolved solids with waters of highest hardness occurring in zone 7. Surface flows from zone 8 are dilute to moderately mineralized and are relatively high in silica concentration. Nitrate concentrations are high in surf Fee flows below sewage outfalls and in ground water contaminated by septic tank effluent and industrial wastes. Average annual sediment yields of 550 to 650 tons per square mile are characteristic of zones 1 and 2 where strip mining has destroyed the forest cover and coal culm is carried into the streams. From agricultural lands on the Martinsburg Shale in zones 5 and 6, annual sediment yields range from 300 to 350 tons per square mile; but from agricultural lands on the siliceous rocks in zone 8 and zones 3 and 4 in the valley and ridge area, the sediment yield ranges from 200 to 250 tons annually per square mile. Lowest annual sediment yields in the basin are in the forested areas of siliceous rocks in zones 2, 3, 4, and 5, and in the sinkhole topography of the limestones in zone 7 where the yield ranges from 30 to 35 tons and 50 to 60 tons per square mile, respectively. The amount of ground water that can be developed in the basin is dependent on the ability of the underlying rocks to yield water to wells. More than 300 gpm (gallons per minute) can be obtained from wells in alluvial materials in the valley bottoms and in some of the limestones where large solution channels and fractures are penetrated by the wells. From 50 to 300 gpm can be obtained from wells in loosely cemented sandstones and in fractured limestones. From 10 to 50 gpm can be developed from wells in the shales and harder sandstones. The most dense rocks will yield from 1 to 10 gpm from fractures and crevices. Most wells yield water from the upper 350 feet of the formation, for this part contains the most fractures or solution channels. Studies show that the velocity at which a contaminant will move downstream in the basin is related to the discharge of the stream at the time. At a stream discharge of about 400 mgd at Pine Grove, a contaminant in Swatara Creek would require about 40 hours to move from Pine Grove to Middletown. As a result of dispersion and dilution, the maximum concentration of the contaminant at Middletown would be less than 20 percent the concentration at Pine Grove under these conditions. An evaluation of the availability of water in the basin indicates that about I,239 mgd enters as precipitation, 630 mgd leaves as streamflow, 580 mgd is evaporated and transpired, and 56 mgd is diverted for use by man. Not all the diversions for man's use are lost to the basin, as about 27 mgd is returned as sewage for reuse. About one-fourth of the waste water is returned to the ground and the remainder to stream drainageways. Of that diverted by man, 11.6 mgd is used for public supply and 44.4 mgd for industrial and private supplies. Diversions of streamflow furnish 86 percent of the public supply and 27 percent of the industrial supply, and ground-water sources yield the remainder. Municipal and private sewage treatment plants are upgrading the waste water in many places, but no provisions are being made for treatment other than natural dilution and assimilation for the 15 mgd of coal-nine drainage in the northern part of the basin. Technology for economic treatment of mine water is not available at this time, although research in this field is being done. Urbanization eastward from Harrisburg and around Lebanon has increased the population density of the basin. Densities of 500 people per square mile and water use exceeding 2.0 mgd per sq mi can be expected in the future. By the year 2000 the population of the basin may increase 60 percent; and if the per capita rate of use increases 0.5 percent per year the domestic requirements for water will be about two times the present use, or 23 mgd. Similarly, if the present 1:4 ratio of domestic use to industrial use of water continues, at least 89 mgd will be needed for industry in the future. Although an increase to twice the present use of water can be foreseen, or 112 mgd, water for the dilution and assimilation of wastes from treatment systems are not included. Providing water for dilution of wastes from treatment plants has not been a problem, but in the future the amounts needed for this purpose will be greater as the population increases. As water becomes more valuable, treatment of sewage wastes to reduce the biochemical-oxygen-demand load by at least 80 to 90 percent will be necessary to conserve water for more productive uses. As much as 100 mgd may be needed for waste dilution in the basin by year 2000. The present trends in suburban and light industrial development will probably persist in the basin. Problems arising through changes in economic value of water, conflicts in use, and alternatives in development are typical of those confronting the manager of a water-resource system.

Pennsylvania↗