Search USGS⌕ Search

SEARCH · Search USGS

Results for “Sustainable Water Resources Management”

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.

408 records · Page 23Linked to original sources

Pilot framework for fish habitat assessments across tidal and non tidal waters in the Patuxent River Basin

As part of the 2014 Chesapeake Bay Watershed Agreement, all Bay States and the District of Columbia have committed to improving the condition of the Bay, which includes a goal to achieve sustainable fisheries. One outcome under that broad goal is improved effectiveness of fish habitat conservation and preservation efforts. In support of that outcome, the U.S. Geological Survey Eastern Ecological Science Center (USGS-EESC) and the National Oceanic and Atmospheric Association’s National Centers for Coastal Ocean Science (NOAA-NCCOS) are actively developing datasets, methods, and analyses to conduct fish habitat assessments in the Chesapeake Bay watershed, guided by recommendations from a regional stakeholder workshop held by the Chesapeake Bay Program’s (CBP) Fish Habitat Action Team (FHAT) in 2018. The joint USGS and NOAA team has been collaborating on methods for conducting inland and estuarine assessments and exploring whether a seamless headwater to estuary assessment could be developed. The goals of this assessment are to benefit both State and Federal fisheries managers, help advance fisheries science, and provide beneficial information for the public. While past national and regional assessments (e.g. the National Fish Habitat Partnership National Assessment) treated inland and estuarine fish habitat conditions separately due to differences in environments, GIS data representation, and data availability, a seamless habitat assessment could be of value for a broad range of stakeholders as many fish species, several of which are invasive or under federal jurisdiction, use habitats across both inland and estuarine waters. This project developed a pilot framework, explored and tested methods necessary for a finer scale, seamless assessment across both inland and estuarine waters, and demonstrated its use. Although there was interest by the CBP FHAT for the generation of a Baywide fish habitat assessment that spanned tidal salt, tidal fresh, warm non-tidal and cold non-tidal waters, there are a myriad of implementation details and considerations around conducting a Baywide assessment across all four of these general habitat areas. Therefore, the practical need to conduct a tributary-specific pilot assessment arose. At the beginning of this pilot process, members of the FHAT were presented with a decision matrix to choose a study basin using factors such as data availability and tributary size. FHAT members chose the Patuxent River basin, which has been relatively well sampled and studied. Several spatial frameworks were considered before selection of an inclusive gridded framework for summary and analysis that represented inland drainage networks and landscape influences as well as estuarine bathymetry. A suite of landscape and in-water stressor variables were summarized into the framework and were largely generalized over time. In order to assess the viability of the framework, we chose to use species distribution modeling for each of the species to test the framework’s ability to predict habitat use of non-tidal resident, estuarine resident, and migratory species. Tessellated darter (Etheostoma olmstedi), American eel (Anguilla rostrata), and white perch (Morone americana) were chosen as illustrative fish species based on data availability, and differences in life history and habitat use. A nested modeling approach, which involved successive model runs at multiple scales (1000m, 100m, and 10m raster grids) was developed to examine differences in variable importance at different spatial scales and to enhance modeling efficiency. For white perch, a complementary modeling analysis was performed for variables available only in estuarine waters. For all testing, an ensemble modeling approach was conducted, using a suite of potential statistical techniques driven by model strength and variable predictive power. The statistical testing that we conducted was intended only to test the framework and modeling approach, and not to definitively predict all habitats where specific fish species might be present. The modeling we conducted to test the framework did have some limitations. For example, the spatial distribution of favorable habitat areas for white perch was likely influenced by the predominance of fish survey locations near the center channel of the river and the use of generalized in-water conditions. For all species, the use of juvenile and adult fish survey data limits the estimation of habitat use to those life stages. Despite such limitations of the data inputs and modeling approach, we found the framework could seamlessly predict fish habitat distribution across freshwater and tidal environments and integrate the influence of landscape stressors with local in-water factors. The developed framework presented to the Sustainable Fisheries Goal Implementation Team (GIT) and FHAT is informative and could potentially be used for other modeling applications in the Chesapeake Bay watershed and elsewhere. In particular the framework and modeling approach lend themselves to evaluating living resource distributions and underlying habitat conditions in shallow tidal waters and beyond, as recommended by the recent Comprehensive Evaluation of System Response (CESR) report from the Chesapeake Bay Program.

Maryland↗

Groundwater simulation and management models for the upper Klamath Basin, Oregon and California

The upper Klamath Basin encompasses about 8,000 square miles, extending from the Cascade Range east to the Basin and Range geologic province in south-central Oregon and northern California. The geography of the basin is dominated by forested volcanic uplands separated by broad interior basins. Most of the interior basins once held broad shallow lakes and extensive wetlands, but most of these areas have been drained or otherwise modified and are now cultivated. Major parts of the interior basins are managed as wildlife refuges, primarily for migratory waterfowl. The permeable volcanic bedrock of the upper Klamath Basin hosts a substantial regional groundwater system that provides much of the flow to major streams and lakes that, in turn, provide water for wildlife habitat and are the principal source of irrigation water for the basin's agricultural economy. Increased allocation of surface water for endangered species in the past decade has resulted in increased groundwater pumping and growing interest in the use of groundwater for irrigation. The potential effects of increased groundwater pumping on groundwater levels and discharge to springs and streams has caused concern among groundwater users, wildlife and Tribal interests, and State and Federal resource managers. To provide information on the potential impacts of increased groundwater development and to aid in the development of a groundwater management strategy, the U.S. Geological Survey, in collaboration with the Oregon Water Resources Department and the Bureau of Reclamation, has developed a groundwater model that can simulate the response of the hydrologic system to these new stresses. The groundwater model was developed using the U.S. Geological Survey MODFLOW finite-difference modeling code and calibrated using inverse methods to transient conditions from 1989 through 2004 with quarterly stress periods. Groundwater recharge and agricultural and municipal pumping are specified for each stress period. All major streams and most major tributaries for which a substantial part of the flow comes from groundwater discharge are included in the model. Groundwater discharge to agricultural drains, evapotranspiration from aquifers in areas of shallow groundwater, and groundwater flow to and from adjacent basins also are simulated in key areas. The model has the capability to calculate the effects of pumping and other external stresses on groundwater levels, discharge to streams, and other boundary fluxes, such as discharge to drains. Historical data indicate that the groundwater system in the upper Klamath Basin fluctuates in response to decadal climate cycles, with groundwater levels and spring flows rising and declining in response to wet and dry periods. Data also show that groundwater levels fluctuate seasonally and interannually in response to groundwater pumping. The most prominent response is to the marked increase in groundwater pumping starting in 2001. The calibrated model is able to simulate observed decadal-scale climate-driven fluctuations in the groundwater system as well as observed shorter-term pumping-related fluctuations. Example model simulations show that the timing and location of the effects of groundwater pumping vary markedly depending on the pumping location. Pumping from wells close (within a few miles) to groundwater discharge features, such as springs, drains, and certain streams, can affect those features within weeks or months of the onset of pumping, and the impacts can be essentially fully manifested in several years. Simulations indicate that seasonal variations in pumping rates are buffered by the groundwater system, and peak impacts are closer to mean annual pumping rates than to instantaneous rates. Thus, pumping effects are, to a large degree, spread out over the entire year. When pumping locations are distant (more than several miles) from discharge features, the effects take many years or decades to fully impact those features, and much of the pumped water comes from groundwater storage over a broad geographic area even after two decades. Moreover, because the effects are spread out over a broad area, the impacts to individual features are much smaller than in the case of nearby pumping. Simulations show that the discharge features most affected by pumping in the area of the Bureau of Reclamation's Klamath Irrigation Project are agricultural drains, and impacts to other surface-water features are small in comparison. A groundwater management model was developed that uses techniques of constrained optimization along with the groundwater flow model to identify the optimal strategy to meet water user needs while not violating defined constraints on impacts to groundwater levels and streamflows. The coupled groundwater simulation-optimization models were formulated to help identify strategies to meet water demand in the upper Klamath Basin. The models maximize groundwater pumping while simultaneously keeping the detrimental impacts of pumping on groundwater levels and groundwater discharge within prescribed limits. Total groundwater withdrawals were calculated under alternative constraints for drawdown, reductions in groundwater discharge to surface water, and water demand to understand the potential benefits and limitations for groundwater development in the upper Klamath Basin. The simulation-optimization model for the upper Klamath Basin provides an improved understanding of how the groundwater and surface-water system responds to sustained groundwater pumping within the Bureau of Reclamation's Klamath Project. Optimization model results demonstrate that a certain amount of supplemental groundwater pumping can occur without exceeding defined limits on drawdown and stream capture. The results of the different applications of the model demonstrate the importance of identifying constraint limits in order to better define the amount and distribution of groundwater withdrawal that is sustainable.

Oregon;California↗

Vegetation changes over 12 years in ungrazed and grazed Conservation Reserve Program Grasslands in the central and southern plains

The Conservation Reserve Program (CRP) established under the 1985 Food Security Act has the fundamental objectives of jointly providing economic support to segments of the agricultural community and conservation of natural resources (Osborn, 1997; Heard and others, 2000). Although soil loss on highly erodable lands was the principal natural resource conservation issue addressed in the 1985 CRP, improving water quality and wildlife habitat both became important considerations as the program evolved (Farmer and others, 1988). For example, Best and others (1997) found that production of young birds on CRP fields in the Midwest was ≥15 times the production on row-crop fields because of improved habitat. The increasing importance of wildlife habitat is reflected in continuing refinement of the Environmental Benefits Index (EBI) used by the U.S. Department of Agriculture (USDA) to quantify the potential benefits of enrolling lands in CRP (Osborn, 1997; Ribaudo and others, 2001). The refinements reflect input furnished by federal, state, and non-government organizations seeking greater wildlife habitat quality on CRP lands (Roseberry and David, 1994; Hughes and others, 1995; Millenbah and others, 1996; Patterson and Best, 1996; Rodgers, 1999; Allen and others, 2001). Refinement in the EBI has changed the types of grasses planted on newly enrolled land. In early CRP signups (1 through 11), 71% of new grassland acres were planted to introduced grasses and legumes [Conservation Practice (CP) 1] while 29% of the acres were planted to native grasses (CP2) (Osborn and others, 1992). By the 27th signup in July 2004, over 34.8 million acres (14 million ha) were enrolled in the CRP. More than 73% of these lands were planted to various mixtures of introduced (CP1) or native (CP2) grasses for a minimal contract period of 10 years (USDA, 2004). Continuation of grass plantings under the 2002 Farm Bill may result in CRP lands furnishing grass dominated cover for 20 or more consecutive years. The species of grass established in seeded grasslands can have a major influence on the potential quality of wildlife habitat where vegetation is maintained over a multi-year period. Different species of grass may have comparable abilities to alleviate soil erosion but furnish dissimilar qualities of wildlife habitat (fig. 1). For example, smooth brome, an introduced cool-season grass (grass species and scientific names are presented in table 1), is highly valued for its erosion control and forage attributes (Casler and Carlson, 1995). Switchgrass, a native warm-season grass, also is valued for its soil and water conservation qualities (Moser and Vogel, 1995) but provides greater benefits for some species of wildlife (Clubine, 1995). The quality of nesting and winter cover for ring-necked pheasants (Phasianus colchicus) furnished by smooth brome on northeastern Colorado CRP lands is inferior to that provided by the taller, more robust switchgrass (Allen, 1994). Characteristics of the agricultural landscape surrounding individual CRP fields also play a role in the wildlife habitat potential of CRP plantings (Weber and others, 2002; Nusser and others, 2004). of native and seeded grasslands change in response to the presence (and absence) of physical disturbances such as fire, grazing, tillage, and haying (Hobbs and Huenneke, 1992; Millenbah and others, 1996; Allen and others, 2001; Renfrew and Ribic, 2001; Swengel and Swengel, 2001). The perpetuation of diversity in species composition and vegetation structure following disturbance sustains desirable habitat for a variety of grassland-dependent wildlife (Hall and Willig, 1994; Barnes and others, 1995; King and Savidge, 1995; Granfors and others, 1996; Herkert and others, 1996; Kurzejeski, 1996; Patterson and Best, 1996; Klute and others, 1997). Undisturbed grasslands have lower grass and forb species diversity, greater amounts of dead plant material, decreased as well as seasonally delayed productivity, and diminished structural diversity of vegetation (Peet and others, 1975; Rice and Parenti, 1978; Butler and Briske, 1988; Campa and Winterstein, 1992). Recommendations for the timing of disturbance to increase grass and forb species diversity range from 3 to 8 years following establishment of seeded grasslands in the northern Great Plains and Midwest (Duebbert and others, 1981; Higgens, 1987; Millenbah and others, 1996). The management interval, however, is affected by climatic conditions, soils, grass species, and management history of the individual stand. We quantified changes in vegetation structure and species composition across the typical 10-year contract period in undisturbed southern and central Great Plains CRP fields (fig. 2) planted to introduced and native grasses. In addition, we compared changes in vegetation in fields grazed during the emergency release of 1996 by comparing conditions prior to grazing and two and four years post grazing relative to changes in similar fields that were not grazed. Documentation of long-term changes in vegetation structure and composition for fields planted to common grass seed mixtures across a wide range of environmental conditions provides information to improve long-term wildlife habitat potential, guide program administration, and define management practices that yield economic benefits to operators while still meeting wildlife and conservation objectives. Emergency grazing provisions of the CRP are controversial. Although grazing can alter vegetative characteristics and reduce habitat quality in the short-term (Temple and others, 1999), periodic disturbance may be necessary to maintain habitat quality, and more information is needed assessing long-term effects of emergency grazing on vegetative structure and species composition.

Conference Paper↗

Sediment distribution and hydrologic conditions of the Potomac aquifer in Virginia and parts of Maryland and North Carolina

Sediments of the heavily used Potomac aquifer broadly contrast across major structural features of the Atlantic Coastal Plain Physiographic Province in eastern Virginia and adjacent parts of Maryland and North Carolina. Thicknesses and relative dominance of the highly interbedded fluvial sediments vary regionally. Vertical intervals in boreholes of coarse-grained sediment commonly targeted for completion of water-supply wells are thickest and most widespread across the central and southern parts of the Virginia Coastal Plain. Designated as the Norfolk arch depositional subarea, the entire sediment thickness here functions hydraulically as a single interconnected aquifer. By contrast, coarse-grained sediment intervals are thinner and less widespread across the northern part of the Virginia Coastal Plain and into southern Maryland, designated as the Salisbury embayment depositional subarea. Fine-grained intervals that are generally avoided for completion of water-supply wells are increasingly thick and widespread northward. Fine-grained intervals collectively as thick as several hundred feet comprise two continuous confining units that hydraulically separate three vertically spaced subaquifers. The subaquifers are continuous northward but merge southward into the single undivided Potomac aquifer. Lastly, far southeastern Virginia and northeastern North Carolina are designated as the Albemarle embayment depositional subarea, where both coarse- and fine-grained intervals are of only moderate thickness. The entire sediment thickness functions hydraulically as a single interconnected aquifer. A substantial hydrologic separation from overlying aquifers is imposed by the upper Cenomanian confining unit. Potomac aquifer sediments were deposited by a fluvial depositional complex spanning the Virginia Coastal Plain approximately 100 to 145 million years ago. Westward, persistently uplifted granite and gneiss source rocks sustained a supply of coarse-grained sand and gravel. Immature, high-gradient braided streams deposited longitudinal bars and channel fills across the Norfolk arch subarea. By contrast, across the Salisbury and Albemarle embayment subareas, mature, medium- to low-gradient meandering streams deposited medium- to coarse-grained channel fills and point bars segregated from fine-grained overbank deposits. The Virginia depositional complex merged northward across the Salisbury embayment subarea with another complex in Maryland. Here, additional sediments were received from schist source rocks that underwent three cycles of initial uplift and rapid erosion followed by crustal stability and erosional leveling. Because of the predominance of coarse-grained sediments, transmissivity, hydraulic conductivity, and regional velocities of lateral flow through the Potomac aquifer are greatest across the Norfolk arch depositional subarea, but decrease progressively northward with increasingly fine-grained sediments. Confining units hydraulically separate the Potomac aquifer from overlying aquifers, as indicated by large vertical hydraulic gradients. By contrast, most of the Potomac aquifer internally functions hydraulically as a single interconnected aquifer, as indicated by uniformly small vertical gradients. Most fine-grained sediments within the aquifer do not hydraulically separate overlying and underlying coarse-grained sediments. Across the Salisbury embayment depositional subarea, however, hydraulic separation among the vertically spaced subaquifers is imposed by the intervening confining units. The Potomac aquifer is the largest and most heavily used source of groundwater in the Virginia Coastal Plain. Water-level declines as great as 200 feet create the potential for saltwater intrusion. Conventional stratigraphic correlation has been generally ineffective at accurately characterizing complexly distributed fluvial sediments that compose the Potomac aquifer. Consequently, the aquifer’s internal hydraulic connectivity and overall hydrologic function have not been well understood. Water-supply planning and development efforts have been hampered, and interpretations of regulatory criteria for allowable water-level declines have been ambiguous. An investigation undertaken during 2010–11 by the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, provides a comprehensive regional description of the spatial distribution of Potomac aquifer sediments and their relation to hydrologic conditions. Altitudes and thicknesses of 2,725 vertical sediment intervals represent the spatial distribution of Potomac aquifer sediments in the Virginia Coastal Plain and adjacent parts of Maryland and North Carolina. Sediment intervals are designated as either dominantly coarse or fine grained and were determined by interpretation of geophysical logs and ancillary information from 456 boreholes. Sediment-interval and borehole summary statistical data indicate regional trends in sediment lithology and stratigraphic continuity, upon which three structurally based and hydrologically distinct sediment depositional subareas are designated. Broad patterns of sediment deposition over time are inferred from published sediment pollen-age data. Discrepancies in previously drawn hydrostratigraphic relations between southeastern Virginia and northeastern North Carolina are partly resolved based on borehole geophysical logs and a recently documented geologic map and corehole. A conceptual model theorizes the depositional history of the sediments and geologically accounts for their distribution. Documented pumping tests of the Potomac aquifer at 197 locations produced 336 values of transmissivity and 127 values of storativity. Based on effective aquifer thicknesses, 296 values of sediment hydraulic conductivity and 113 values of sediment specific storage are calculated. Vertical hydraulic gradients are calculated from 9,479 pairs of water levels measured between November 17, 1953, and October 4, 2011, in 129 closely spaced pairs of wells. Borehole sediment-interval and related data provide a means to achieve high yielding production wells in the Potomac aquifer by site-specific targeting of drilling operations toward water-bearing coarse-grained sand and gravel. Advance knowledge of the potential of different parts of the aquifer also aids in planning optimal groundwater-development areas. Depositional subareas further provide a possible context for resource management. Current (2013) regulatory limits on water-level declines are relative to top surfaces of subdivided upper, middle, and lower Potomac aquifers across the entire Virginia Coastal Plain, but have the potential to exceed the same limit relative to a single undivided Potomac aquifer. By contrast, designation of the sediments as a single aquifer in the Norfolk arch and Albemarle embayment subareas—and as a series of vertically spaced subaquifers and intervening confining units in the Salisbury embayment subarea—best reflects understanding of the Potomac aquifer and can avoid the potential for excessive water-level declines. Simulation modeling to evaluate effects of groundwater withdrawals could be designed similarly, including vertical discretization and (or) zonation of the Potomac aquifer based on depositional subareas and a geostatistical distribution of aquifer properties derived from borehole sediment-interval data. Further resource-management information needs extend beyond the developed part of the Potomac aquifer, particularly across the Northern Neck and Middle Peninsula where only the shallowest part of the aquifer is known, and include structural aspects such as faults, basement bedrock, and the Chesapeake Bay impact crater.

Maryland, North Carolina, Virginia↗

Hydrogeology, chemical characteristics, and water sources and pathways in the zone of contribution of a public-supply well in San Antonio, Texas

In 2001, the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey initiated a series of studies on the transport of anthropogenic and natural contaminants (TANC) to public-supply wells (PSWs). The main goal of the TANC project was to better understand the source, transport, and receptor factors that control contaminant movement to PSWs in representative aquifers of the United States. Regional- and local-scale study areas were selected from within existing NAWQA study units, including the south-central Texas Edwards aquifer. The local-scale TANC study area, nested within the regional-scale NAWQA study area, is representative of the regional Edwards aquifer. The PSW selected for study is within a well field of six production wells. Although a single PSW was initially selected, because of constraints of well-field operation, samples were collected from different wells within the well field for different components of the study. Data collected from all of the well-field wells were considered comparable because of similar well construction, hydrogeology, and geochemistry. An additional 38 PSWs (mostly completed in the confined part of the aquifer) were sampled throughout the regional aquifer to characterize water quality. Two monitoring well clusters, with wells completed at different depths, were installed to the east and west of the well field (the Zarzamora and Timberhill monitoring well clusters, respectively). One of the monitoring wells was completed in the overburden to evaluate potential hydrologic connectivity with the Edwards aquifer. Geophysical and flowmeter logs were collected from one of the well-field PSWs to determine zones of contribution to the wellbore. These contributing zones, associated with different hydrogeologic units, were used to select monitoring well completion depths and groundwater sample collection depths for depth-dependent sampling. Depth-dependent samples were collected from the PSW from three different depths and under three different pumping conditions. Additionally, selected monitoring wells and one of the well-field PSWs were sampled several times in response to a rainfall and recharge event to assess short-term (event-scale) temporal variations in water quality. For comparison purposes, groundwater samples were categorized as being from regional aquifer PSWs, from the well field (wellhead samples), from the monitoring wells (excluding the overburden well), from the overburden well, from the PSW depth-dependent sampling, and from temporal sampling. Groundwater samples were analyzed for inorganic, organic, isotopic, and age-dating tracers to characterize geochemical conditions in the aquifer and provide understanding of the mechanisms of mobilization and movement of selected constituents from source areas to a PSW. Sources, tracers, and conditions used to assess water quality and processes affecting the PSW and the aquifer system included (1) carbonate host rock composition; (2) physicochemical constituents; (3) major and trace element concentrations; (4) saturation indices with respect to minerals in aquifer rocks; (5) elemental ratios, such as magnesium to calcium ratios, that are indicative of water-rock interaction processes; (6) oxidation-reduction conditions; (7) nutrient concentrations, in particular nitrate concentrations; (8) the isotopic composition of nitrate, which can point to specific nitrate sources; (9) strontium isotopes; (10) stable isotopes of hydrogen and oxygen; (11) organic contaminant concentrations, including pesticides and volatile organic compounds; (12) age tracers, apparent-age distribution, and dissolved gas data used in age interpretations; (13) depth-dependent water chemistry collected from the PSW under different pumping conditions to assess zones of contribution; and (14) temporal variability in groundwater composition from the PSW and selected monitoring wells in response to an aquifer recharge event. Geochemical results indicate that the well-field and monitoring well samples were largely representative of groundwater in the regional confined aquifer. Constituents of concern in the Edwards aquifer for the long-term sustainability of the groundwater resource include the nutrient nitrate and anthropogenic organic contaminants. Nitrate concentrations (as nitrogen) for regional aquifer PSWs had a median value of 1.9 milligrams per liter, which is similar to previously reported values for the regional aquifer. Nitrate-isotope compositions for groundwater samples collected from the well-field PSWs and monitoring wells had a narrow range, with values indicative of natural soil organic values. A comparison with historical nitrate-isotope values, however, suggests that a component of nitrate in groundwater from biogenic sources might have increased over the last 30 years. Several organic contaminants (the pesticide atrazine, its degradate deethylatrazine, trichloromethane (chloroform; a drinking-water disinfection byproduct), and the solvent tetrachloroethene (PCE)) were widely distributed throughout the regional aquifer and in the local-scale TANC study area at low concentrations (less than 1 microgram per liter). Higher concentrations of PCE were detected in samples from the well-field PSWs and Zarzamora monitoring wells relative to the regional aquifer PSWs. The urban environment is a likely source of contaminants to the aquifer, and these results indicate that one or more local urban sources might be supplying PCE to the Zarzamora monitoring wells and the well-field wells. Samples from the well field also had high concentrations of chloroform relative to the monitoring wells and regional aquifer PSWs. For samples from the regional aquifer PSWs, the most frequently detected organic contaminants generally decreased in concentration with increasing well depth. Deeper wells might intercept longer regional flow paths with higher fractions of older water or water recharged in rural recharge areas in the western part of the aquifer that have been less affected by anthropogenic contaminants. A scenario of hypothetical contaminant loading was evaluated by using results from groundwater-flow-model particle tracking to assess the response of the aquifer to potential contamination. Results indicate that the aquifer responds quickly (less than 1 year to several years) to contaminant loading; however, it takes a relatively long time (decades) for concentrations to reach peak values. The aquifer also responds quickly (less than 1 year to several years) to the removal of contaminant loading; however, it also takes a relatively long time (decades) to reach near background concentrations. Interpretation of geochemical age tracers in this well-mixed karst system was complicated by contamination of a majority of measured tracers and complexities of extensive mixing. Age-tracer results generally indicated that groundwater samples were composed of young, recently recharged water with piston-flow model ages ranging from less than 1 to 41 years, with a median of 17 years. Although a piston-flow model is typically not valid for karst aquifers, the model ages provide a basis for comparing relative ages of different samples and a reference point for more complex hydrogeologic models for apparent-age interpretations. Young groundwater ages are consistent with particle-tracking results from hydrogeologic modeling for the local-scale TANC study area. Age-tracer results compared poorly with other geochemical indicators of groundwater residence time and anthropogenic effects on water quality, indicating that hydrogeologic conceptual models used in groundwater age interpretations might not adequately account for mixing in this karst system. Groundwater samples collected from the well field under a variety of pumping conditions were relatively homogeneous and well mixed for numerous geochemical constituents (with the notable exception of age tracers). Groundwater contributions to the PSW were dominated by well-mixed, relatively homogeneous groundwater, typical of the regional confined aquifer. Zones of preferential flow were determined for the PSW, but groundwater samples from different stratigraphic units were not geochemically distinct. Variations in chemical constituents in response to a rainfall and aquifer recharge event occurred but were relatively minor in the PSW and monitoring wells. This observation is consistent with the hypothesis that the response to individual recharge events in the confined aquifer, unless intersecting conduit flow paths, might be attenuated by mixing processes along regional flow paths. Results of this study are consistent with the existing conceptual understanding of aquifer processes in this karst system and are useful for water-resource development and management practices.

Texas↗

Environmental Setting and Effects on Water Quality in the Great and Little Miami River Basins, Ohio and Indiana

The Great and Little Miami River Basins drain approximately 7,354 square miles in southwestern Ohio and southeastern Indiana and are included in the more than 50 major river basins and aquifer systems selected for water-quality assessment as part of the U.S. Geological Survey's National Water-Quality Assessment Program. Principal streams include the Great and Little Miami Rivers in Ohio and the Whitewater River in Indiana. The Great and Little Miami River Basins are almost entirely within the Till Plains section of the Central Lowland physiographic province and have a humid continental climate, characterized by well-defined summer and winter seasons. With the exception of a few areas near the Ohio River, Pleistocene glacial deposits, which are predominantly till, overlie lower Paleozoic limestone, dolomite, and shale bedrock. The principal aquifer is a complex buried-valley system of sand and gravel aquifers capable of supporting sustained well yields exceeding 1,000 gallons per minute. Designated by the U.S. Environmental Protection Agency as a sole-source aquifer, the Buried-Valley Aquifer System is the principal source of drinking water for 1.6 million people in the basins and is the dominant source of water for southwestern Ohio. Water use in the Great and Little Miami River Basins averaged 745 million gallons per day in 1995. Of this amount, 48 percent was supplied by surface water (including the Ohio River) and 52 percent was supplied by ground water. Land-use and waste-management practices influence the quality of water found in streams and aquifers in the Great and Little Miami River Basins. Land use is approximately 79 percent agriculture, 13 percent urban (residential, industrial, and commercial), and 7 percent forest. An estimated 2.8 million people live in the Great and Little Miami River Basins; major urban areas include Cincinnati and Dayton, Ohio. Fertilizers and pesticides associated with agricultural activity, discharges from municipal and industrial wastewater-treatment and thermoelectric plants, urban runoff, and disposal of solid and hazardous wastes contribute contaminants to surface water and ground water throughout the study area. Surface water and ground water in the Great and Little Miami River Basins are classified as very hard, calcium-magnesiumbicarbonate waters. The major-ion composition and hardness of surface water and ground water reflect extensive contact with the carbonate-rich soils, glacial sediments, and limestone or dolomite bedrock. Dieldrin, endrin, endosulfan II, and lindane are the most commonly reported organochlorine pesticides in streams draining the Great and Little Miami River Basins. Peak concentrations of the herbicides atrazine and metolachlor in streams commonly are associated with post-application runoff events. Nitrate concentrations in surface water average 3 to 4 mg/L (milligrams per liter) in the larger streams and also show strong seasonal variations related to application periods and runoff events. Ambient iron concentrations in ground water pumped from aquifers in the Great and Little Miami River Basins often exceed the U.S. Environmental Protection Agency Secondary Maximum Contaminant Level (300 micrograms per liter). Chloride concentrations are below aesthetic drinking-water guidelines (250 mg/L), except in ground water pumped from low-yielding Ordovician shale; chloride concentrations in sodium-chloriderich ground water pumped from the shale bedrock can exceed 1,000 mg/L. Some of the highest average nitrate concentrations in ground water in Ohio and Indiana are found in wells completed in the buried-valley aquifer; these concentrations typically are found in those parts of the sand and gravel aquifer that are not overlain by clay-rich till. Atrazine was the most commonly detected herbicide in private wells. Concentrations of volatile organic compounds in ground water generally were below Federal drinking-water standards, except near areas of known or suspected contamination. Evaluation of fish and macroinvertebrate community performance in streams and rivers draining the Great and Little Miami River Basins indicates that most streams meet basic aquatic-life-use criteria set by the Ohio Environmental Protection Agency for warmwater habitat. Stream reaches whose biological community performance meet aquatic-lifeuse criteria defined for exceptional warmwater habitat are found in Twin Creek, the Upper Great Miami River, the Little Miami River, and the Whitewater River Basins. Other streams have exhibited significant improvements in biological community performance (and water quality)'that are attributed primarily to reduced pollutant loadings from wastewater-treatment plants upgraded since 1972. Four hydrogeomorphic regions were delineated in the Great and Little Miami River Basins based on distinct and relatively homogeneous natural characteristics. Primary features used to delineate the hydrogeomorphic regions include bedrock geology, surficial geology, physiography, hydrology, soil types, and vegetation. These four regions Till Plains, Drift Plains/Unglaciated, Interlobate, and Fluvial are used in the Great and Little Miami River Basins study to assess the influence of natural features of the environmental setting on surface- and ground-water quality.

Indiana, Ohio↗

Hydrogeomorphic segments and hydraulic microhabitats of the Niobrara River, Nebraska— With special emphasis on the Niobrara National Scenic River

The Niobrara River is an ecologically and economically important resource in Nebraska. The Nebraska Department of Natural Resources’ recent designation of the hydraulically connected surface- and groundwater resources of the Niobrara River Basin as “fully appropriated” has emphasized the importance of understanding linkages between the physical and ecological dynamics of the Niobrara River so it can be sustainably managed. In cooperation with the Nebraska Game and Parks Commission, the U.S. Geological Survey investigated the hydrogeomorphic and hydraulic attributes of the Niobrara River in northern Nebraska. This report presents the results of an analysis of hydrogeomorphic segments and hydraulic microhabitats of the Niobrara River and its valley for the approximately 330-mile reach from Dunlap Diversion Dam to its confluence with the Missouri River. Two spatial scales were used to examine and quantify the hydrogeomorphic segments and hydraulic microhabitats of the Niobrara River: a basin scale and a reach scale. At the basin scale, digital spatial data and hydrologic data were analyzed to (1) test for differences between 36 previously determined longitudinal hydrogeomorphic segments; (2) quantitatively describe the hydrogeomorphic characteristics of the river and its valley; and (3) evaluate differences in hydraulic microhabitat over a range of flow regimes among three fluvial geomorphic provinces. The statistical analysis of hydrogeomorphic segments resulted in reclassification rates of 3 to 28 percent of the segments for the four descriptive geomorphic elements. The reassignment of classes by discriminant analysis resulted in a reduction from 36 to 25 total hydrogeomorphic segments because several adjoining segments shared the same ultimate class assignments. Virtually all of the segment mergers were in the Canyons and Restricted Bottoms (CRB) fluvial geomorphic province. The most frequent classes among hydrogeomorphic segments, and the dominant classes per unit length of river, are: a width-restricted valley confinement condition, sinuous-planview pattern, irregular channel width, and an alternate bar configuration. The Niobrara River in the study area flows through a diversity of fluvial geomorphic settings in its traverse across northern Nebraska. In the Meandering Bottoms (MB) fluvial geomorphic province, river discharge magnitudes are low, and the valley exerts little control on the channel-planview pattern. Within the CRB province, the river flows over a diversity of geologic formations, and the valley and river narrow and expand in approximate synchronicity. In the Braided Bottoms (BB) fluvial geomorphic province, the river primarily flows over Cretaceous Pierre Shale, the valley and channel are persistently wide, and the channel slope is generally uniform. The existence of vegetated islands and consequent multithread channel environments, indicated by a higher braided index, mostly coincided with reaches having gentler slopes and less unit stream power. Longitudinal hydrology curves indicate that the flow of the Niobrara River likely is dominated by groundwater as far downstream as Norden. Unit stream power values in the study area vary between 0 and almost 2 pounds per foot per second. Within the MB province, unit stream power steadily increases as the Niobrara gains discharge from groundwater inflow, and the channel slope steepens. The combination of steep slopes, a constrained channel width, and persistent flow within the CRB province results in unit stream power values that are between three and five times greater than those in less confined segments with comparable or greater discharges. With the exception of hydrogeomorphic segment 3, which is affected by Spencer Dam, unit stream power values in the BB province are generally uniform. Channel sinuosity values in the study area varied generally between 1 and 2.5, but with locally higher values measured in the MB province and at the entrenched bedrock meanders of hydrogeomorphic segment 18 in the CRB province. The differences in channel morphology and hydraulic geometries between fluvial geomorphic provinces are evident in the types, relative abundance, and response of hydraulic microhabitats to changing discharges. The four gaging stations chosen for hydraulic microhabitat analysis are distributed among three different fluvial geomorphic provinces. In the MB province, the smaller channel and lower discharges resulted in the dominance of shallow and intermediate-depth hydraulic environments with the vast majority of hydraulic microhabitat restricted to shallow categories even during upper-decile discharges. In the CRB province, intermediate depth hydraulic conditions, particularly intermediate-swift, dominate over all ranges of discharge. Hydraulic microhabitat conditions were most diverse in the BB province, with most hydraulic microhabitat categories present over the entire range of discharges analyzed. The calculated differences in hydraulic microhabitat distributions, abundance, and adjustments between streamflow-gaging stations were the result of differences in physical structure of the channel and subsequent channel hydraulic geometry. At the reach scale, field measurements made in water years 2008 and 2009 in four study reaches within the Scenic Reach were used to (1) characterize the elevation and geomorphic processes associated with fluvial landforms, (2) build hydraulic geometry relations, (3) examine flow hydraulics over a range of discharges, and (4) examine the types and responses of hydraulic microhabitats to a range of flow magnitudes. Four landform groups were identified and named in order of increasing elevation: low flood plains, intermediate flood plains, low terraces, and high terraces. The terraces were poorly characterized because the surveys did not extend across the full width of the alluvial valley bottom. The two lowest fluvial landforms are likely active in the modern hydroclimatic regime. Sediment samples obtained in the study reaches indicate that the primary bed material in the active channel ranged in size from coarse silt to coarse sand. Grain-size distributions from samples also indicate that the bed of the Niobrara River among the study reaches coarsens and has increasing grainsize variability in the downstream direction. Values of at-a-station hydraulic geometry exponents indicate that the Niobrara River in the study reaches adjusts its geometry to changing discharges primarily through increases in flow depth and velocity. Relations at one cross section indicated that, at least locally, changes in width were also an important channel adjustment mechanism. Hydraulic behavior over the range of flows measured was not consistent among all study reaches, but two general modes of hydraulic behavior were observed in the reaches with substantial coverage of the bed by fine sediment. At the Sunny Brook and Muleshoe study reaches, average boundary-shear stress remained approximately constant, and hydraulic resistance decreased, for discharges below 900 cubic feet per second (ft 3 /s). Above 900 ft 3 /s, average boundary shear stress and hydraulic resistance both increased. The Rock Barn study reach did not exhibit the same two-mode hydraulic behavior observed at the Sunny Brook and Muleshoe reaches. The coincident increase in boundary shear stress above 900 ft 3 /s observed at the Sunny Brook and Muleshoe study reaches represents a potential hydraulic threshold above which bedload transport rates were likely to increase markedly. No consistent bed-adjustment pattern (scour or fill) was identified in the study reaches over the range of flows or over the measurement season. Analysis of hydraulic microhabitats over the range of discharges measured at the study reaches indicates that some percentage of most habitat niche categories was available for at least one discharge condition, but the majority of hydraulic habitat available was within the intermediate-swift and deepswift habitat niche categories. Deep-swift conditions dominated nearly all study reaches under all measured discharge conditions. Slight differences in habitat distributions were observed between the study reaches with substantial coverage of the bed by fine sediment—Sunny Brook, Muleshoe, and Rock Barn—and the bedrock-dominated reach, Crooked Creek. Although the four study reaches occupy three different hydrogeomorphic segments, the types, relative abundance, and response of hydraulic microhabitat niche distributions to changing discharge conditions generally were similar among all reaches.

Nebraska↗

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗

Toward improved decision-support tools for Delta Smelt management actions

The Collaborative Science and Adaptive Management Program (CSAMP) has endorsed a goal of reversing the recent downward trajectory of the Delta Smelt population within 5-10 generations, with the long-term aim of establishing a self-sustaining population. An ambitious agenda of management actions is planned, and more management actions are being considered. This White Paper furthers one of the recommendations in the 2019 Delta Smelt Science Plan – the need to predict the potential ecological effects of taking a management action. Existing statistical models can be highly informative in assessing the response of Delta Smelt to changing system conditions and management actions. However, management actions can shift or alter conditions in ways that models based on analysis of historical data may not be able to represent, and short-term or localized effects may be missed with models designed to assess effects at the population level. Decision support tools (DSTs) are computer-based tools developed to assist decision-making, often combining computationally intensive analysis and spatial mapping of environmental relationships. DSTs can be used in planning processes that evaluate an array of actions, such as in Structured Decision Making (SDM), where DSTs are needed to compare among alternatives. DSTs can also be used to explore the potential effects of different approaches to implementing management actions. The goal of this White Paper is to identify plausible options for DSTs that could be developed for future use to evaluate management actions that seek to either reverse the decline of Delta Smelt or minimize or mitigate the effects of other water management actions. Different types of management actions lead to different needs for DSTs. This White Paper was developed using three types of actions currently being considered to enhance the Delta Smelt population: Supplementation with Hatchery Fish, Summer-Fall Habitat, and Food Enhancement actions. These three management actions target different parts of the estuary and different processes, with a variety of possible metrics to gauge performance. Three DSTs are proposed that collectively address management questions related to the management actions considered, with each requiring a slightly different set of processes to be included and producing an array of outputs at varying spatial and temporal scales: DST 1. Modeling Fish Movement, Survival, and Reproduction Across Their Range. This DST can address management questions that require information about Delta Smelt spatial distribution and movement. • DST 1 could be used to compare conditions with and without management actions in place, how the management action performs among different types of water years (with varied flow and associated abiotic conditions), and to assess relative change with different variations and strategies of the management actions. • DST 2. Changes in Habitat Conditions and Delta Smelt Response. This DST is intended to evaluate combinations of conditions that are considered to provide suitable habitat for Delta Smelt, and Delta Smelt response. Delta Smelt habitat is generally described as open water with low salinity (0 to 6), turbidity of at least 12 NTU, suitable temperature conditions, and sufficient food availability to support growth. • DST 3. Regional Effects of Food Subsidy. This DSTs seeks to evaluate effectiveness of food enhancement actions by providing information on responses of the immediate targets of the action (i.e., phytoplankton or zooplankton) and tracing those to projected growth responses of Delta Smelt. There is not a single DST that adequately addresses management questions relevant to all management actions, although there is some overlap in the management questions each of the three DSTs can address. For each of the DSTs a substantial foundation of models and approaches already exists and modeling has already been applied to several of the management actions described. However, a number of outstanding issues remain for further development of the proposed DSTs. These are summarized in this White Paper together with potential approaches that could be applied or tested. Some components for the DSTs are already available and thus development could be relatively easy. However, for several of the topics identified there are gaps in knowledge that currently limit formulation of model structure and process representations. This presents challenges to readily incorporate some needed mechanisms into the models. Eleven next steps, aligned with relevant DSTs, are outlined. The next steps vary in their complexity or technical ‘lift’ required. Many build on existing work, or methods and approaches that have already been developed or are underway, while others require additional thinking to establish a viable approach. Some interim utility for decisions could be gained during initial development of the DSTs with further features added over time. Development of a DST requires engagement of both managers and scientists. Identifying the outputs and resolution needed for management purposes early in development of any DST is essential for effective pursuit of next steps and suitable approaches to address challenges. Dialog between managers and technical experts also informs what process-based simulation can do, and what tradeoffs are acceptable to meet a given purpose. To further develop the DSTs outlined here for application in the estuary requires: - Engagement of a committed group of technical experts with appropriate expertise. - The development of a coordinated workplan including appropriate project management and tracking. - Dialog between potential users (i.e., managers and policy makers) and technical experts. - Resources to pursue DST development including personnel and computational resources. This White Paper demonstrates the potential for moving toward DSTs for a variety of management actions in support of Delta Smelt that include mechanistic representations of physical and biological processes. Through focused effort from technical experts, managers and policy makers, DSTs can be developed to provide quantitative predictions of management effects on the ecosystem, targeting the changes the management actions seek to achieve, how these effects compare to ambient conditions, and how the effects vary among water year types or with timing and location of actions. Importantly, solid foundations exist which can be leveraged, refined, and built upon to specifically inform current and future management decisions.

White Paper↗

Mobile Bay

Mobile Bay is the largest bay found in Alabama’s coastal area (Handley et al., 2007). It was named an Estuary of National Significance in 1995 under the U.S. Environmental Protection Agency’s (EPA) National Estuary Program (NEP), and its Comprehensive Conservation Management Plan was completed in 2002. Mobile Bay is 1,070 km 2 (413 miles 2 ) in area and 51 km (32 miles) long, making it the sixth largest estuary in the continental United States (Mobile Bay NEP, 2008). Its ecosystem provides habitat for more than 300 species of birds, 310 species of fish, 68 species of reptiles, 57 species of mammals, 40 species of amphibians, and 15 species of shrimp (Mobile Bay NEP, 1997). Mobile Bay lies between the Mississippi and Atlantic Flyways (Mobile Bay NEP, 2003). Commercial and residential development and industrial use is heavy in the Mobile Bay area. Although local growth and industrial markets support the Mobile Bay area economy, the resulting environmental damage to the very ecosystem upon which they depend remains a threat to the environment, economy, and population. The Mobile Bay ecosystem boasts high biological diversity and productivity and supports many freshwater and saltwater species of recreational and commercial importance. The great diversity of Mobile Bay reflects the diversity of Alabama, which is home to the largest number of different plant and animal species of all states east of the Mississippi River (Stein, 2002), and is bolstered by the unique climate and geographic conditions surrounding the bay. Freshwater inflow from the Mobile-Tensaw River Delta, ranging from 60,000 to 3,700,000 gallons per second (Wallace, 1996), mixes with saltwater from the Gulf of Mexico, which enters Mobile Bay via wind and tides (Burgan and Engle, 2006). Because of the unique conditions surrounding Mobile Bay, including shallow waters, a dynamic climate, and artificial hydrologic modifications—such as the construction of the Mobile Bay Causeway in the 1920s, which serves as an unintentional barrier between Delta waters north of the Causeway and saline waters south of the Causeway, the salinity of Mobile Bay is highly variable. Mobile Bay receives an average of 165 cm (65 inches) of rain per year from tropical storms, summer thunderstorms, and winter cold fronts (Stout et al., 1998). The climate and geography that have made Mobile Bay so rich in resources have also contributed to the threats surrounding its ecosystem. The extensive amount of rain in Mobile Bay creates large amounts of runoff, polluting the waters with fertilizers, chemicals, sediment, oil, trash, and sewage (Mobile Bay NEP, 1997). Tourism, ecotourism, recreational and commercial fishing, recreational boating, shipping, and chemical, pulp, and paper production are significant industries in Mobile Bay and the surrounding areas. Despite the approximate \$3 billion and 55,000 jobs these industries bring into the community (Alabama Tourism Department, 2010), the growth, development, and environmental stress they create are major threats to the Mobile Bay ecosystem. Among the nation’s states, Alabama ranks fifth in number of different species (144 endemic species), second in number of extinctions that have already occurred (90 extinct species) and fourth in number of species at risk for extinction (14.8% at risk out of 4,533 total species; Stein, 2002). Twenty-one of these threatened and endangered species are found in Mobile Bay, whose brackish waters provide a nursery area for many species of vertebrates and invertebrates. Some of these species include the Alabama sturgeon, Gulf sturgeon, heavy pigtoe mussel, inflated heel-splitter mussel, West Indian manatee, Alabama beach mouse, Perdido beach mouse, Alabama red-bellied turtle, gopher tortoise, Kemp’s ridley sea turtle, green sea turtle, loggerhead sea turtle, eastern indigo snake, flatwoods salamander, piping plover, red-cockaded woodpecker, and wood stork. Habitat loss underlies the decline of some bird species in Mobile Bay, and large mammals such as the red wolf, Florida panther, and Florida black bear are no longer found in the area. However, some rare species, such as the swallow-tailed kite, sandhill crane, and gopher tortoise can still be found (Duke and Kruczynski, 1992). The value of wetlands in Mobile Bay and the rest of the Gulf of Mexico is still being investigated. Although various monetary valuations of wetlands exist, critics remark that undervaluation of wetlands is inevitable (Mobile Bay NEP, 2008) and that estimates often do not place appropriate value on ecological services (Mitsch and Gosselink, 2000). Additionally, many estimates account only for anthropogenic values. One estimate concludes that one acre of wetlands performs \$3,000 worth of water purification each year (Mobile Bay NEP, 1997). With more than 76,890 hectares (190,000 acres) of wetlands in the Mobile Bay area, that equates to a value exceeding one-half billion dollars every year. Tourism, fishing, boating, production, and shipping are significant industries in the Mobile Bay area. More than 90% of fish landed in recreational and commercial fishing in the bay depend on bay habitat, including wetlands, for life requirements (Mobile Bay NEP, 1997). The Port of Mobile is Alabama’s only ocean-ship port (Mobile Bay NEP, 2008). Baldwin County, on the eastern side of the bay, experienced a population increase of 75% from 1990 to 2007, with an 89% increase in housing units (Mobile Bay NEP, 2008). Development and industry support the Mobile Bay economy, but they depend on the continued health, sustainability, and production of the water and living resources of the Mobile Bay ecosystem. Wetland loss, along with other forms of environmental degradation, remains a threat to the Mobile Bay ecosystem and Mobile Bay’s socioeconomic foundation.

Alabama↗

Comprehensive water quality of the Boulder Creek Watershed, Colorado, during high-flow and low-flow conditions, 2000

Executive Summary The Boulder Creek Watershed, Colorado, is 1160 square kilometers in area and ranges in elevation from 1480 to 4120 meters above sea level. Streamflow originates primarily as snowmelt near the Continental Divide, and thus discharge varies seasonally and annually (Chapter 1). Most of the water in Boulder Creek is diverted for domestic, agricultural, and industrial use. Some diverted water is returned to the creek as wastewater effluent and by ditch returns, and additional water enters as groundwater and by transbasin diversions. These diversions and returns lead to complex temporal and spatial variations in discharge. The variations in discharge, along with natural factors such as geology and climate, and anthropogenic factors such as wastewater treatment, agriculture, mining, and urbanization, can affect water chemistry. As with many watersheds in the American West, dependable water quality and sufficient water supply are issues facing local water managers and users. Detailed water-quality and sediment sampling allows the identification of sources and sinks of chemical constituents and an understanding of the processes at work in a river system. This study, the most comprehensive water-quality analysis performed for Boulder Creek to date, was a cooperative effort of the U.S. Geological Survey (USGS) and the city of Boulder. Geographic information systems and modeling programs were used to delineate watershed boundaries, land cover, and geology (Chapter 2). During high-flow (June 2000) and low-flow (October 2000) conditions, researchers evaluated 226 water-quality variables, including basic water-quality indicators (Chapter 3), major ions and trace elements (Chapter 4), wastewater-derived organic compounds (Chapter 5), and pesticides (Chapter 6). Discharge (Chapter 1) and bed-sediment particle size and mineralogy (Chapter 7) were also evaluated. This cooperative study was facilitated by the Boulder Area Sustainability Information Network (BASIN), which provides public access to environmental information about the Boulder Creek Watershed on a website, www.basin.org. In addition to the USGS and city of Boulder data, researchers at the Institute of Arctic and Alpine Research at the University of Colorado provided water chemistry data for the headwaters of North Boulder Creek, upstream of the reach of the USGS/city of Boulder sampling sites (Chapter 8). Snowmelt produces high flows in Boulder Creek in late spring to early summer (Chapter 1). Because precipitation falling in the headwaters is very dilute (specific conductance about 5 microsiemens per centimeter), most chemical constituents are present in lower concentrations during high flows (Chapters 3, 4, 5, 6, and 8). However, concentrations of some constituents, such as total suspended solids (Chapter 3) and organic carbon (Chapter 5), increase during the spring snowmelt flush. The upper basin, which consists of alpine, subalpine, montane, and foothills regions west of the mouth of Boulder Canyon, is underlain by Precambrian igneous and metamorphic rocks (Chapter 1). Major dissolved inorganic constituents in headwater sites were found to be enriched by factors of 10 to 20 relative to precipitation; this is consistent with minor weathering of the local crystalline bedrock (Chapter 4). Some anthropogenic input is observed in the headwaters; precipitation introduces nitrogen derived from fossil fuel combustion and agricultural activities (Chapter 8). The lower basin, which consists of the plains region east of the mouth of Boulder Canyon, is underlain by Mesozoic sedimentary rock and Quaternary alluvium, and has substantially more anthropogenic sources. Concentrations of most dissolved inorganic constituents increased in the lower basin. Differentiation between natural and anthropogenic sources of some dissolved constituents is difficult because both sources contribute to the water composition in this region. The increase of most major constituents (bicarbonate, calcium, chloride, magnesium, sodium, and sulfate) is consistent with weathering of the underlying sedimentary bedrock (Chapter 4). It is likely that anthropogenic loading of constituents in this reach occurs during storm events. Fecal coliform concentrations were variable and in some cases exceeded state standards, primarily during low-flow conditions (Chapter 3). Effluent from Boulder’s 75th Street Wastewater Treatment Plant (WWTP) has a substantial impact on the water chemistry of lower Boulder Creek. The WWTP increases the concentrations of nutrients such as nitrogen and phosphorus (Chapter 3), major ions and trace metals (Chapter 4), and organic carbon (Chapter 5) in Boulder Creek. The effluent contained a spike in gadolinium, a rare earth element that is ingested for magnetic resonance imaging as a contrasting agent and then excreted to the urban wastewater system. The effluent also contained trace organic compounds such as surfactants, pharmaceuticals, hormones (Chapter 5), and pesticides (Chapter 6), which also were detected at downstream Boulder Creek sites. Water chemistry of Boulder Creek downstream of the WWTP is largely controlled by the degree of dilution of the wastewater effluent, which varies depending on the baseflow of Boulder Creek, the volume of wastewater effluent, and depletion by agricultural diversions. Coal Creek, a tributary of Boulder Creek, contains wastewater effluent from four additional WWTPs, and increases the load of many constituents in Boulder Creek. In addition to the impact from wastewater effluent, lower Boulder Creek is affected by agricultural land use. Eleven of 84 analyzed pesticides were detected in Boulder Creek or its inflows, primarily in the eastern section of the watershed (Chapter 6). This collaborative study provides an in-depth evaluation of the hydrology, water chemistry, and sediment mineralogy of North Boulder Creek, Middle Boulder Creek, Boulder Creek, and major inflows. The detailed sampling and analysis in this report provide a baseline for future reference, as well as information on the effect of land use and geology on water chemistry.

Colorado↗

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↗