Search USGSSearch

SEARCH · Search USGS

Results for “Systems Analysis in Forest Resources.”

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.

70 records · Page 4Linked to original sources

Executive summary. In Second State of the Carbon Cycle Report (SOCCR2): A Sustained Assessment Report

Central to life on Earth, carbon is essential to the molecular makeup of all living things and plays a key role in regulating global climate. To understand carbon’s role in these processes, researchers measure and evaluate carbon stocks and fluxes. A stock is the quantity of carbon contained in a pool or reservoir in the Earth system (e.g., carbon in forest trees), and a flux is the direction and rate of carbon’s transfer between pools (e.g., the movement of carbon from the atmosphere into forest trees during photosynthesis). This document, the Second State of the Carbon Cycle Report (SOCCR2), examines the patterns of carbon stocks and fluxes—collectively called the “carbon cycle.” Emphasis is given to these patterns in specific sectors (e.g., agriculture and energy) and ecosystems (e.g., forests and coastal waters) and to the response of the carbon cycle to human activity. The purpose of SOCCR2 is to assess the current state of the North American carbon cycle and to present recent advances in understanding the factors that influence it. Concentrating on North America—Canada, the United States, and Mexico—the report describes carbon cycling for air, land, inland waters (streams, rivers, lakes, and reservoirs), and coastal waters (see Figure ES.1, p. 23). The questions framing the publication A U.S. Carbon Cycle Science Plan (Michalak et al., 2011) inspired development of three slightly modified questions that guide SOCCR2’s content and focus on North America in a global context: How have natural processes and human actions affected the global carbon cycle on land, in the atmosphere, in the ocean and other aquatic systems, and at ecosystem interfaces (e.g., coastal, wetland, and urban-rural)? How have socioeconomic trends affected atmospheric levels of the primary carbon-containing gases, carbon dioxide (CO2) and methane (CH4)? How have species, ecosystems, natural resources, and human systems been impacted by increasing greenhouse gas (GHG) concentrations, associated changes in climate, and carbon management decisions and practices? SOCCR2 synthesizes the most recent understanding of carbon cycling in North America, assessing new carbon cycle findings and information, the state of knowledge regarding core methods used to study the carbon cycle, and future research needed to best inform carbon management and policy options. Focusing on scientific developments in the decade since the First State of the Carbon Cycle Report (SOCCR1; CCSP 2007), SOCCR2 summarizes the past, current, and projected state of carbon sources, sinks, and natural processes, as well as contributions by human activities. In addition to CO2 and CH4, the report sometimes discusses nitrous oxide (N2O), a GHG associated with activities and processes that affect fluxes of carbon gases.1 SOCCR2 also describes improvements in analysis tools; developments in decision support; and new insights into ecosystem carbon cycling, human causes of changes in the carbon cycle, and social science perspectives on carbon. Since publication of SOCCR1, coordinated research from agencies in the three North American countries has enabled innovative observational, analytical, and modeling capabilities to further advance understanding of the North American carbon cycle (see Appendix D: Carbon Measurement Approaches and Accounting Frameworks, p. 834). Some of the report’s main conclusions, based on the Key Findings of each chapter, are highlighted in Box ES.1, Main Findings of SOCCR2, p. 24.

Report

Estimating water temperatures in small streams in western Oregon using neural network models

Artificial neural network models were developed to estimate water temperatures in small streams using data collected at 148 sites throughout western Oregon from June to September 1999. The sites were located on 1st-, 2nd-, or 3rd-order streams having undisturbed or minimally disturbed conditions. Data collected at each site for model development included continuous hourly water temperature and description of riparian habitat. Additional data pertaining to the landscape characteristics of the basins upstream of the sites were assembled using geographic information system (GIS) techniques. Hourly meteorological time series data collected at 25 locations within the study region also were assembled. Clustering analysis was used to partition 142 sites into 3 groups. Separate models were developed for each group. The riparian habitat, basin characteristic, and meteorological time series data were independent variables and water temperature time series were dependent variables to the models, respectively. Approximately one-third of the data vectors were used for model training, and the remaining two-thirds were used for model testing. Critical input variables included riparian shade, site elevation, and percentage of forested area of the basin. Coefficient of determination and root mean square error for the models ranged from 0.88 to 0.99 and 0.05 to 0.59 oC, respectively. The models also were tested and validated using temperature time series, habitat, and basin landscape data from 6 sites that were separate from the 142 sites that were used to develop the models. The models are capable of estimating water temperatures at locations along 1st-, 2nd-, and 3rd-order streams in western Oregon. The model user must assemble riparian habitat and basin landscape characteristics data for a site of interest. These data, in addition to meteorological data, are model inputs. Output from the models include simulated hourly water temperatures for the June to September period. Adjustments can be made to the shade input data to simulate the effects of minimum or maximum shade on water temperatures.

Water-Resources Investigations Report

Synthesis of nutrient and sediment data for watersheds within the Chesapeake Bay drainage basin

Nutrient and sediment data collected by Federal and state agencies from 1972 through 1992 at 1,058 surface-water sites in nontidal parts of the Chesapeake Bay Basin were compiled into a large database. Adequate nutrient, sediment, and streamflow data were not available to compute annual loads for all sites because water-quality monitoring at many of the sites was either short term or noncontinuous or because stream-flow was not measured. Annual nutrient and sediment loads were calculated at a total of 127 sites. Annual loads of dissolved nitrate were calculated for 108 sites, but total nitrogen loads could be calculated for only 48 of these sites because ammonia plus organic nitrogen data were not available for many of these 108 sites. Annual loads of total phosphorus were calculated for 99 sites, and annual loads of suspended sediment were calculated for 33 sites. Loads could be calculated for only a very few sites in the Juniata River Basin (a tributary to the Susquehanna River), the York River Basin, the middle and lower reaches of the James River, and the nontidal parts of the eastern shore of the Bay. Geographic Information System (GIS) spatial data sets of land use, physiographic province, rock type, and watershed delineation were compiled for the entire Chesapeake Bay Basin (approximately 64,000 square miles). The nutrient- and sediment-yield were evaluated with respect to land use, physiographic province, rock type, and hydrologic characteristics. During years that the mean streamflow was about equal to the long-term mean streamflow, the Susquehanna River contributed about 50 percent of the freshwater, 66 percent of the total nitrogen, and 40 percent of the total phosphorus transported by tributaries to the Bay. Nutrient and sediment data were available for less than 18 percent of the predominantly agricultural areas underlain by siliciclastic rock and for less than 35 percent of the predominantly agricultural areas underlain by either carbonate rock or unconsolidated rock. Nutrient and sediment data were available for about 91 percent of the predominantly forested areas underlain by siliciclastic rock. Spatial and temporal gaps in the water-quality data and GIs data sets limited some data analysis. Correlations of annual yields or nutrients and sediment with respect to land use, physiographic province, and rock type indicated (1) basins with larger percentages of agricultural land had larger nutrient and sediment yields, (2) basins with larger percentages of forest land had smaller nutrient and sediment loads, (3) the largest total nitrogen yields were from agricultural basins underlain by carbonate rock, (4) yields or nutrients from urban basins were substantially less than yields from agricultural basins, and (5) basins with small amounts of agricultural and urban land had relatively small nutrient and sediment yields.

Chesapeake Bay drainage basin

Precipitation isotopes link regional climate patterns to water supply in a tropical mountain forest, eastern Puerto Rico

Like many mountainous areas in the tropics, watersheds in the Luquillo Mountains of eastern Puerto Rico have abundant rainfall and stream discharge and provide much of the water supply for the densely populated metropolitan areas nearby. Projected changes in regional temperature and atmospheric dynamics as a result of global warming suggest that water availability will be affected by changes in rainfall patterns. It is essential to understand the relative importance of different weather systems to water supply to determine how changes in rainfall patterns, interacting with geology and vegetation, will affect the water balance. To help determine the links between climate and water availability, stable isotope signatures of precipitation from different weather systems were established to identify those that are most important in maintaining streamflow and groundwater recharge. Precipitation stable isotope values in the Luquillo Mountains had a large range, from fog/cloud water with δ 2 H, δ 18 O values as high as +12 ‰, −0.73 ‰ to tropical storm rain with values as low as −127 ‰, −16.8 ‰. Temporal isotope values exhibit a reverse seasonality from those observed in higher latitude continental watersheds, with higher isotopic values in the winter and lower values in the summer. Despite the higher volume of convective and low-pressure system rainfall, stable isotope analyses indicated that under the current rainfall regime, frequent trade -wind orographic showers contribute much of the groundwater recharge and stream base flow. Analysis of rain events using 20 years of 15 -minute resolution data at a mountain station (643 m) showed an increasing trend in rainfall amount, in agreement with increased precipitable water in the atmosphere, but differing from climate model projections of drying in the region. The mean intensity of rain events also showed an increasing trend. The determination of recharge sources from stable isotope tracers indicates that water supply will be affected if regional atmospheric dynamics change trade- wind orographic rainfall patterns in the Caribbean.

Luquillo Mountains

Water-quality assessment of the Lower Susquehanna River Basin, Pennsylvania and Maryland: sources, characteristics, analysis and limitations of nutrient and suspended-sediment data, 1975-90

This report describes analyses of available information on nutrients and suspended sediment collected in the Lower Susquehanna River Basin during water years 1975-90. Most of the analyses were applied to data collected during water years 1980-89. The report describes the spatial and temporal availability of nutrient and suspended-sediment data and presents a preliminary concept of the spatial and temporal patterns of concentrations and loads within the basin. Where data were available, total and dissolved forms of nitrogen and phosphorus species from precipitation, surface water, ground water, and springwater, and bottom material from streams and reservoirs were evaluated. Suspended-sediment data from streams also were evaluated. The U.S. Geological Survey National Water Information System (NWIS) database was selected as the primary database for the analyses. Precipitation-quality data from the National Atmospheric Deposition Program (NADP) and bottom-material-quality data from the National Uranium Resource Evaluation (NURE) were used to supplement the water-quality data from NWIS. Concentrations of nutrients were available from 3 precipitation sites established for longterm monitoring purposes, 883 wells (854 synoptic areal survey sites and 29 project and research sites), 23 springs (17 synoptic areal survey sites and 6 project and research sites), and 894 bottom-material sites (840 synoptic areal survey sites and 54 project and research sites). Concentrations of nutrients and (or) suspended sediment were available from 128 streams (36 long-term monitoring sites, 51 synoptic areal survey sites, and 41 project and research sites). Concentrations of nutrients and suspended sediment in streams varied temporally and spatially and were related to land use, agricultural practices, and streamflow. A general north-to-south pattern of increasing median nitrate concentrations, from 2 to 5 mg/L, was detected in samples collected in study unit streams. In streams that drain areas dominated by agriculture, concentrations of nutrients and suspended sediment tend to be elevated with respect to those found in areas of other land-use types and are related to the amount of commercial fertilizer and animal manure applied to the area drained by the streams. Animal manure is the dominant source of nitrogen for the streams in the lower, agricultural part of the basin. Concentrations of nutrients in samples from wells varied with season and well depth and were related to hydrogeologic setting. Median concentrations of nitrate were 2.5 and 3.5 mg/L for wells drawing water at depths of 0 to 100 ft and 101 to 200 ft, respectively. The lowest median concentrations for nitrate in ground water from wells were generally found in siliciclastic-bedrock, forested settings of the Ridge and Valley Physiographic Province, and the highest were found in carbonate-bedrock agricultural settings of the Piedmont Physiographic Province. Twenty-five percent of the measurements from wells in carbonate rocks in the Piedmont Physiographic Province exceeded the Pennsylvania drinking-water standard. An estimate of mass balance of nutrient loads within the Lower Susquehanna River Basin was produced by combining the available information on stream loads, atmosphericdeposition loads, commercial-fertilizer applications, animal-manure production, privateseptic-system nonpoint-source loads, and municipal and industrial point-source loads. The percentage of the average annual nitrate load carried in base flow of streams in the study unit ranged from 45 to 76 percent, and the average annual phosphorus load carried in base flow ranged from 20 to 33 percent. Average annual yields of nutrients and suspended sediment from tributary basins are directly related to percentage of drainage area in agriculture and inversely to drainage area. Information required to compute loads of nitrogen and phosphorus were available for all sources except atmospheric deposition, for which only nitrogen data were available. Atmospheric deposition is the dominant source of nitrogen for the mostly forested basins draining the upper half of the study unit. The estimate of total annual nitrogen load to the study unit from precipitation is 98.8 million pounds. Nonpoint and point sources of nutrients were estimated. Nonpoint and point sources combined, including atmospheric deposition, provide a potential annual load of 390 million pounds of nitrogen and 79.5 million pounds of phosphorus. The range of percentages of the estimated nonpoint and point sources that were measured in the stream was 20 to 47 percent for nitrogen and 6 to 14 percent for phosphorus. On the average, the Susquehanna River discharges 141,000 pounds of nitrogen and 7,920 pounds of phosphorus to the Lower Susquehanna River reservoir system each year. About 98 percent of the nitrogen and 60 percent of the phosphorus passes through the reservoir system. Interpretations of available water-quality data and conclusions about the water quality of the Lower Susquehanna River Basin were limited by the scarcity of certain types of water-quality data and current ancillary data. A more complete assessment of the water quality of the basin with respect to nutrients and suspended sediment would be enhanced by the availability of additional data for multiple samples over time from all water environments; samples from streams in the northern and western part of the basin; samples from streams and springs throughout the basin during high base-flow or stormflow conditions; and information on current land-use, and nutrient loading from all types of land-use settings.

Maryland, Pennsylvania

Multi-resource analysis: A proof of concept study of natural resource tradeoffs in the Piceance Basin, Colorado, using the net resources assessment (NetRA) decision support tool

Executive Summary The U.S. Geological Survey (USGS) is developing a multi-resource analysis (MRA) line of products to inform land-use decision makers. Specifically, MRA products will integrate scientific information, include considerations for natural resource interrelations, and quantify the effects of resource management decisions in biophysical, economic, and societal terms. As part of the establishment of the MRA, the USGS, in collaboration with the University of New Mexico, has developed the Net Resources Assessment (NetRA) decision support tool. As a proof of concept analysis, the NetRA was applied to the Piceance basin in Colorado in a hypothetical example to illustrate how resource managers could use the NetRA to consider tradeoffs of natural resources among alternative development plans and land cover patterns within a geographic region. The NetRA is a policy-relevant approach to assess the availability of multiple natural resources. It is an analytical toolset that may be used to examine the spatiotemporal relations between development of energy and mineral resources and delivery of biological natural resources. The NetRA operates at multiple map scales and contains a set of integrated, compatible submodels with specific data requirements for natural resource stocks, engineering economics, biophysical, and ecological data for ecosystem services stocks, market prices, regulations, and nonmarket values. The NetRA includes an explicit process to consider the interdependence between development and conservation, which is a crucial consideration in land-management and land-use decisions. The NetRA is used to estimate an expected net resource value (NRV). The NRV is the expected, present value, economic benefit from the extraction of a resource (for example, natural gas) minus the total cost of production, which is the aggregation of the development, production, and social costs. Social costs include private costs plus any external costs. There can be external social benefits associated with natural gas production, such as increased demand for locally produced goods and increased employment in the local area through backward and forward linkages of natural gas production. The NRV is used to compare development outcomes (scenarios) from a range of exploration and development plans for cumulative energy production. The Piceance basin application of the NetRA uses the NRV to assess the tradeoff between continuous natural gas extraction and the effects to the local populations of Odocoileus hemionus (mule deer) and aquatic species and to consumptive water uses for an area the size and resolution of a USGS energy resource assessment unit. In the proof of concept simulation, the 2.9-square-mile-area of USGS oil and gas assessment unit 50200263 (Piceance basin continuous gas unit of the Mesaverde Total Petroleum System) was gridded into 588 cells. From this area, seven clusters with potential for development and three that cannot be developed were identified; the three clusters that cannot be developed were identified as wilderness study areas, areas of critical environmental concern, and national forests. On the basis of these criteria, there are 118 cells unsuitable for development in the oil and gas assessment unit: 84 are in national forests, 23 are areas of critical environmental concern, and 11 are wilderness study areas. The remaining cells in the oil and gas assessment unit can be developed on both private and public lands. Two scenarios were considered that are distinguished as plan 1 and plan 2. Plan 1 keeps the amount of land disturbance unchanged and limits the number of development locations to 140 grid cells for the production period, which constrains the amount of the energy resources available for development; the plan requires the usage of the Bureau of Land Management (BLM) unsuitability criteria. Plan 2 also limits the number of development locations to 140 grid cells for the production period but provides a constant volume of energy production by increasing the density of well pads within the cells. The effects of plan 2 to the NRV when there are five wells per pad and five pads per square mile happen mostly in the first 5 years of development, even though the effects on the population of mule deer continue in later years. This outcome is the result of the upfront development and investment costs and the initial effect to the ecosystem services.

Colorado

Numerical simulation of streamflow distribution, sediment transport, and sediment deposition along Long Beach Creek in Northeast Missouri

This report presents the results of a study conducted by the U.S. Geological Survey in cooperation with the Missouri Department of Conservation to describe the hydrology, sediment transport, and sediment deposition along a selected reach of Long Branch Creek in Macon County, Missouri. The study was designed to investigate spatial and temporal characteristics of sediment deposition in a remnant forested riparian area and compare these factors by magnitude of discharge events both within and outside the measured range of flood magnitudes. The two-dimensional finite-element numerical models RMA2-WES and SED2D-WES were used in conjunction with measured data to simulate streamflow and sediment transport/deposition characteristics during 2-, 5-, 10-, and 25-year recurrence interval floods. Spatial analysis of simulated sediment deposition results indicated that mean deposition in oxbows and secondary channels exceeded that of the remaining floodplain areas during the 2-, 5-, 10-, and 25-year recurrence interval floods. The simulated mass deposition per area for oxbows and secondary channels was 1.1 to 1.4 centimeters per square meter compared with 0.1 to 0.60 centimeters per square meter for the remaining floodplain. The temporal variability of total incremental floodplain deposition during a flood was found to be strongly tied to sediment inflow concentrations. Most floodplain deposition, therefore, occurred at the beginning of the streamflow events and corresponded to peaks in sediment discharge. Simulated total sediment deposition in oxbows and secondary channels increased in the 2-year through 10-year floods and decreased in the 25- year flood while remaining floodplain deposition was highest for the 25-year flood. Despite increases in sediment inflows from the 2-year through 25-year floods, the retention ratio of sediments (the ratio of floodplain deposition to inflow load) was greatest for the 5-year flood and least for the 25-year flood. The decrease in retention ratio at greater flows is likely the result of higher velocities on the floodplain, resulting in higher bed shear stress, greater suspension time of deposited material, and greater sediment transport through the system. Simulated sediment deposition was most sensitive to sediment inflow concentrations and modification of floodplain roughness—factors that can be controlled through management practices. The increase in floodplain sediment deposition resulting from a simulated increase in vegetation density (increase in floodplain roughness from a Manning's n of 0.11 to 0.12) was 142,000 kilograms, or 6.5 percent for a 10-year recurrence interval flood. This increase was comparable to total oxbow and secondary channel deposition mass in the simulations, but would result in a mean increase in floodplain deposition thickness of only 0.025 centimeter. The hydrodynamic model results show the importance of the secondary channels and meander cutoff channels in this system because these areas quickly bring floodwaters and sediment to areas not close to the main channel. The meander cutoff channels in the simulation also effectively decrease flow and velocities in some main channel sections thereby affecting sediment deposition in the vicinity of these features.

Missouri

Water quality in the Little Sac River basin near Springfield, Missouri, 1999-2001

The Little Sac River, north of Springfield, Missouri, flows through mainly agricultural and forest land. However, the quality of the river water is a concern because the river flows into Stockton Lake, which is a supplemental drinking water source for Springfield. Large bacterial densities and nutrient concentrations are primary concerns to the water quality of the river. A 29-river mile reach of the Little Sac River is on the 1998 list of waters of Missouri designated under section 303(d) of the Federal Clean Water Act because of fecal coliform densities larger than the Missouri Department of Natural Resources standard (hereinafter referred to as Missouri standard) of 200 colonies per 100 milliliters for whole-body contact recreation. During an investigation of the water quality in the Little Sac River by the U.S. Geological Survey, in cooperation with the Watershed Committee of the Ozarks, fecal coliform bacteria densities exceeded the Missouri standard (the standard applies from April 1 through October 31) in one sample from a site near Walnut Grove. At other sites on the Little Sac River, the Missouri standard was exceeded in two samples and equalled in one sample upstream from the Northwest Wastewater Treatment Plant (NW WTP) and in one sample immediately downstream from the NW WTP. Effluent from the NW WTP flows into the Little Sac River. Annually from April 1 through October 31, the effluent is disinfected to meet the Missouri standard for whole-body contact recreation. Fecal coliform bacteria densities in samples collected during this period generally were less than 100 colonies per 100 milliliters. For the rest of the year when the effluent was not disinfected, the bacteria densities in samples ranged from 50 (sample collected on November 1, 2000) to 10,100 colonies per 100 milliliters (both counts were non-ideal). When the effluent was disinfected and the fecal coliform bacteria density was small, samples from sites upstream and downstream from the NW WTP had a bacteria density larger than the density in the effluent. Other sources of bacteria are likely to be present in the study area in addition to the NW WTP. These potential sources include effluent from domestic septic systems and animal wastes. Nutrient concentrations in the Little Sac River immediately downstream from the NW WTP were affected by effluent from the NW WTP and possibly other sources. At two sites upstream from the NW WTP, median nitrite plus nitrate concentrations were 1.1 and 1.4 milligrams per liter. The median nitrite plus nitrate concentration for the effluent from the NW WTP was 6.4 milligrams per liter, and the median concentration decreased downstream in the Little Sac River to 2.2, 1.2, and 0.56 milligrams per liter. The effects of the effluent from the NW WTP on the water quality of the Little Sac River downstream from the NW WTP were reflected in an increase in discharge (effluent from the NW WTP can be as much as 50 percent of the flow at the site about 1.5 river miles downstream from the NW WTP), an increase in specific conductance values, an increase in several inorganic constituent concentrations, including calcium, magnesium, and sulfate, and a large increase in sodium and chloride concentrations. The effluent from the NW WTP seemed to have no effect on the pH value, temperature, and dissolved oxygen concentrations in the Little Sac River. Results of repetitive element polymerase chain reaction (rep-PCR) pattern analysis indicated that most Escherichia coli ( E. coli ) bacteria in water samples probably were from nonhuman sources, such as horses and cattle. The rep-PCR pattern analysis indicated that horses were an important source of E. coli downstream from the NW WTP, which was consistent with horses pastured adjacent to the sampling site. Fecal coliform bacteria loads increased upstream from the NW WTP from the most upstream site to the site immediately upstream from the NW WTP. Loads in the effluent from the NW WTP and also those in the Little Sac River downstream from the NW WTP were dependent on the treatment of the effluent. When the effluent was not disinfected, the loads in the effluent increased from those upstream. Downstream in the Little Sac River, the loads decreased, but then increased at the most downstream site. The increase may be a result of increased loads from tributaries and other sources not sampled during the study. When the effluent was disinfected, fecal coliform bacteria loads were less than loads in samples from the Little Sac River downstream from the NW WTP.

Water-Resources Investigations Report

Residence times and nitrate transport in ground water discharging to streams in the Chesapeake Bay Watershed

One of the major water-quality problems in the Chesapeake Bay is an overabundance of nutrients from the streams and rivers that discharge to the Bay. Some of these nutrients are from nonpoint sources such as atmospheric deposition, agricultural manure and fertilizer, and septic systems. The effects of efforts to control nonpoint sources, however, can be difficult to quantify because of the lag time between changes at the land surface and the response in the base-flow (ground water) component of streams. To help resource managers understand the lag time between implementation of management practices and subsequent response in the nutrient concentrations in the base-flow component of streamflow, a study of ground-water discharge, residence time, and nitrate transport in springs throughout the Chesapeake Bay Watershed and in four smaller watersheds in selected hydrogeomorphic regions (HGMRs) was conducted. The four watersheds were in the Coastal Plain Uplands, Piedmont crystalline, Valley and Ridge carbonate, and Valley and Ridge siliciclastic HGMRs. A study of springs to estimate an apparent age of the ground water was based on analyses for concentrations of chlorofluorocarbons in water samples collected from 48 springs in the Chesapeake Bay Watershed. Results of the analysis indicate that median age for all the samples was 10 years, with the 25th percentile having an age of 7 years and the 75th percentile having an age of 13 years. Although the number of samples collected in each HGMR was limited, there did not appear to be distinct differences in the ages between the HGMRs. The ranges were similar between the major HGMRs above the Fall Line (modern to about 50 years), with only two HGMRs of small geographic extent (Piedmont carbonate and Mesozoic Lowland) having ranges of modern to about 10 years. The median values of all the HGMRs ranged from 7 to 11 years. Not enough samples were collected in the Coastal Plain for comparison. Spring samples showed slightly younger water under wet conditions than under dry conditions. The apparent age of water from wells, springs, and other ground-water discharge points in the four targeted watersheds was modern to 60 years, which was similar to the apparent ages from the spring study. In the Pocomoke River Watershed in the Coastal Plain Uplands HGMR, the apparent age of ground-water samples ranged from 0 to 60 years; the ages in the vicinity of the streams ranged from 0 to 23 years. The apparent ages of ground water in the Polecat Creek Watershed in the Piedmont crystalline HGMR ranged from 2 to 30 years. The apparent ages of water from wells in the Muddy Creek Watershed in the Valley and Ridge carbonate HGMR ranged from 10 to 20 years (except for a single sample that was 45 years). The ages in the East Mahantango Creek Watershed in the Valley and Ridge siliciclastic HGMR ranged from 0 to 50 years. The distribution in apparent age of water from wells in the targeted watersheds, however, generally is older than that for water from the springs. The median age of water from wells in the Muddy Creek Watershed, for example, was 15 years, compared to 11 years for the water from the springs in that watershed, and less than 10 years for water from all springs in the spring study. The similarity in the ranges in apparent age of water from the wells and from the springs shows that the samples from the targeted watersheds and springs have bracketed the range of apparent ages that would be expected in the shallow ground-water-flow systems throughout the Chesapeake Bay Watershed. The apparent age of water from individual wells does not necessarily represent the entire distribution of ages of the discharging ground water, and it is this distribution of ages that affects the response of nutrient concentrations in stream base flow. Nutrient-reduction scenarios were modeled for two watersheds for which the distribution of apparent ground-water ages was available, the East Mahantango Creek Watershed in the Valley and Ridge siliciclastic HGMR and the Locust Grove Watershed in the Coastal Plain Uplands HGMR. A nutrient-reduction scenario was created for East Mahantango Creek, where the average residence time was determined to be approximately 10 years on the basis of the output of particle tracking from a ground-water-flow model. This scenario showed decreases of nearly 50 percent in base-flow concentrations of nitrate in streams within the first year after the reduction in nitrogen input; smaller reductions in nitrate concentration occurred in each subsequent year. A second scenario for that same watershed, in which the same 10-year average residence time was assumed and an exponential model was used for analysis, showed that a 50-percent reduction in base-flow concentrations of nitrate could take up to 5 years. For the Locust Grove Watershed, in which an average residence time of 32 years was assumed, simulation with the exponential model showed that it may take more than 20 years to achieve a 50-percent reduction in base-flow concentra-tions of nitrate. Although it was not possible to construct such scenarios for all watersheds, these examples show the range of possible responses to changes in nutrient inputs in two very different types of watersheds. Findings from this study include information on factors that affect ground-water age, spatial distribution of ages, and nitrogen transport. In the East Mahantango Creek Watershed and the Polecat Creek Watershed, the residence time varied spatially depending on the position of the flow path, and temporally depending on the recharge conditions. Generally, ground water in areas near the stream had short residence times and the water in upland areas had longer residence times. Water traveling through deep layers had longer residence times than water traveling through shallow layers, and residence times were faster under high recharge conditions than low recharge conditions. Ground water in the Pocomoke Watershed exhibits a similar pattern: younger water discharges to small order streams in headwater basins and older water discharges to larger streams near the basin outlet. Factors affecting nitrogen transport in ground water include spatial and temporal variation in input sources, ground-water age, and aquifer processes that lead to denitrification. Spatial and temporal variations in nitrogen sources affect all the watersheds. Tributaries with higher inputs of nitrogen have higher concentrations in stream base flow. Areas where nitrogen application rates have increased over time show an age-nitrate relation in ground-water samples. The age-nitrate relation can be affected by denitrification, which occurs in Pocomoke and East Mahantango Creeks but is not evident in Polecat and Muddy Creeks. In East Mahantango Creek, the level of denitrification is significant in water with residence times greater than 20 years, but because this is a small component of overall ground-water discharge to a stream, it may not remove a significant quantity of nitrogen from the system. Denitrification in Pocomoke Creek is significant and appears to affect mostly older water discharging to streams. Therefore, if most of the nitrogen entering these two streams is associated with the discharge of younger ground water, denitrification may not greatly affect the overall nitrogen delivery to these streams. Other findings of this study show that nitrate in ground water discharging along preferential flow paths may not be affected by natural processes, such as denitrification or uptake by riparian vegetation. Seeps to swales and ditches beneath the north uplands at Polecat Creek indicate a shallow water table and discharge of young ground water whereas the absence of such seeps on the south side indicates a deep water table and a lack of young ground water. Similarly, discharge at the base of the slope and to the valley wetland south of the creek but not north of the creek indicates a different role for the riparian forest on the two sides of the creek. In many of the systems where water discharges at the base of slopes to wetlands, ditches have been dug to drain the valley. Such drainage circumvents possible removal of nitrate by riparian vegetation. Because ground-water residence times do not appear directly related to the HGMRs, the targeting of management practices will achieve the most rapid response in water quality if directed at 1) watersheds with large agricultural sources of nitrate, 2) areas with the shortest ground-water-flow paths and 3) areas not affected by significant denitrification. The fastest response in stream base-flow concentrations of nitrogen to implementation of management practices would be to implement practices in those areas with the highest loads rather than attempt to target practices on the basis of HGMR stratification. Overall findings of the study indicate that 1) ground-water contributions to nitrogen in streamflow are significant, 2) some response to management practices should be evident in base-flow concentrations of nitrogen and loads within 1 to 5 years in watersheds with the shortest average residence times, but response time may be closer to 20 years in watersheds with longer average ground-water residence times, 3) the majority of the response in ground-water discharge to any changes in management practices will be distributed over a 10-year time period even in the watersheds with the fastest response times, and 4) given that half the streamflow is from ground-water discharge and the other half is runoff or soil water, about 90 percent of total water being discharged to a stream will be less than about a decade old; therefore, full implementation of nutrient reductions may result in improved streamwater quality in about a decade. In the more-likely scenario of gradual source reduction, the reduction in concentrations of nitrate in streams and aquifers would take longer than the examples shown here.

Chesapeake Bay watershed

Data compilation and assessment for water resources in Pennsylvania state forest and park lands

As a result of a cooperative study between the U.S. Geological Survey and the Pennsylvania Department of Conservation and Natural Resources (PaDCNR), available electronic data were compiled for Pennsylvania state lands (state forests and parks) to allow PaDCNR to initially determine if data exist to make an objective evaluation of water resources for specific basins. The data compiled included water-quantity and water-quality data and sample locations for benthic macroinvertebrates within state-owned lands (including a 100-meter buffer around each land parcel) in Pennsylvania. In addition, internet links or contacts for geographic information system coverages pertinent to water-resources studies also were compiled. Water-quantity and water-quality data primarily available through January 2007 were compiled and summarized for site types that included streams, lakes, ground-water wells, springs, and precipitation. Data were categorized relative to 35 watershed boundaries defined by the Pennsylvania Department of Environmental Protection for resource-management purposes. The primary sources of continuous water-quantity data for Pennsylvania state lands were the U.S. Geological Survey (USGS) and the National Weather Service (NWS). The USGS has streamflow data for 93 surface-water sites located in state lands; 38 of these sites have continuous-recording data available. As of January 2007, 22 of these 38 streamflow-gaging stations were active; the majority of active gaging stations have over 40 years of continuous record. The USGS database also contains continuous ground-water elevation data for 32 wells in Pennsylvania state lands, 18 of which were active as of January 2007. Sixty-eight active precipitation stations (primarily from the NWS network) are located in state lands. The four sources of available water-quality data for Pennsylvania state lands were the USGS, U.S. Environmental Protection Agency, Pennsylvania Department of Environmental Protection (PaDEP), and the Susquehanna River Basin Commission. The water-quality data, which were primarily collected after 1970, were summarized by categorizing the analytical data for each site into major groups (for example, trace metals, pesticides, major ions, etc.) for each type (streams, lakes, ground-water wells, and springs) of data compiled. The number of samples and number of detections for each analyte within each group also were summarized. A total of 410 stream sites and 205 ground-water wells in state lands had water-quality data from the available data sets, and these sites were well-distributed across the state. A total of 107 lakes and 47 springs in state lands had water-quality data from the available data sets, but these data types were not well-distributed across the state; the majority of water-quality data for lakes was in the western or eastern sections of the state and water-quality data for springs was primarily located in the central part of the Lower Susquehanna River Valley. The most common types of water-quality data collected were major ions, trace elements, and nutrients. Physical parameters, such as water temperature, stream discharge, or water level, typically were collected for most water-quality samples. Given the large database available from PaDEP for benthic macroinvertebrates, along with some data from other agencies, there is very good distribution of benthic-macroinvertebrate data for state lands. Benthic macroinvertebrate samples were collected at 1,077 locations in state lands from 1973 to 2006. Most (980 samples) of the benthic-macroinvertebrate samples were collected by PaDEP as part of the state assessment of stream conditions required by the Clean Water Act. Data compiled in this report can be used for various water-resource issues, such as basin-wide water-budget analysis, studies of ecological or instream flow, or water-quality assessments. The determination of an annual water budget in selected basins is best supported by the availab

Pennsylvania

Drain tiles and groundwater resources: Understanding the relations

Executive Summary Drainage for agricultural production over the past 150 years has been an integral component of human-driven change to Minnesota’s rural landscapes. Benefits of drainage Historically, poorly drained soils across much of the State would often remain saturated or flooded after spring snowmelt, preventing timely farm operations such as tilling and planting crops (Arneman, 1963). Installation of agricultural drainage, both surface ditches and subsurface drainage, accelerated transport of water off farm fields and imparted producers higher crop yields (Beauchamp, 1987; Stoner and others, 1993). Agricultural drainage offered many other benefits such as preventing crop drown out, aerating the soil profile for improved plant growth, limiting surface runoff and soil erosion, and allowing farmers better access to croplands (Fausey and others, 1987). Without agricultural drainage on much of Minnesota’s croplands, it would have been difficult to realize high enough crop yields to remain economically viable. Environmental concerns While drainage of Minnesota croplands provided the benefits mentioned above, several environmental concerns result. These include wetland loss, degradation of downstream water quality, and reduced [potential for] recharge. Early agricultural drainage efforts (pre-20th century) led to the disappearance of much of Minnesota’s natural wetlands. Increased focus on preventing or mitigating wetland loss over the last 50 years has helped curtail further losses, even as agricultural drainage proceeds. Prior to establishment of Minnesota statehood, wetlands accounted for more than 10 million acres in Minnesota, including prairie wetlands, peatlands, and forest wetlands that comprised approximately 19 percent of the total land area (Palmer, 1915; King, 1980). In 2018, only half of Minnesota’s pre-settlement wetlands remain, mostly in parts of the State that have not experienced widespread drainage, such as northern Minnesota. Water-quality monitoring has shown that agricultural drainage, in particular the practice of subsurface drainage, provides a direct flow path for nutrient (nitrogen and soluble phosphorus) losses to surface water resources. The negative consequences of agricultural drainage on surface water quality are well documented (for example, Dinnes and others, 2002; Kladivko and others, 2004; Richards and others, 2008; Rozemeijer and others, 2010; Schottler and others, 2013). Agricultural basins with a high percentage of agricultural drainage have been implicated as part of the cause of the Gulf of Mexico hypoxia zone due to excessive nitrogen export (Goolsby and Battaglin, 2001; Randall and Mulla, 2001). The connection of hydrological effects of agricultural subsurface drainage on groundwater recharge and aquifers, on the other hand, has not been well-established. Agricultural subsurface drainage intercepts infiltrating water below croplands and directly discharges the water to nearby surface waters. However, the size of the water balance shift from drained water that would have evapotranspired or run off the land to drained water that would recharge underlying aquifers has been poorly characterized (Schuh, 2008). Drain Tiles and Groundwater Given the poor accounting of subsurface drainage effects on groundwater resources, the Minnesota Ground Water Association (MGWA) deemed it imperative that we document these effects so that groundwater resources in agricultural regions with substantial drainage can be effectively managed. This white paper documents the relations of drain tiles and groundwater resources and discusses the historical significance of agricultural drainage practices, the recognized positive benefits and potential negative consequences of agricultural drainage practices, and the gaps in understanding of the connections between agricultural drainage and groundwater resources. The major messages emerged from the findings of this white paper are: Complex history. Minnesota has a long history of agricultural drainage, spanning over 150 years. Agricultural drainage, and the eventual widespread usage of subsurface drainage, can be separated into at least four distinct periods of time: (a) early drainage to get water off the land, pre-20th century; (b) the boom and bust era (1900-1945); (c) postwar resurgence of subsurface drainage and early conservation efforts (1945-1960); and, (d) emergence of the environmental movement (1960 to present). The State’s regulatory framework that both allowed for drainage and controlled its usage during these periods is complicated and has been governed by a patchwork of both State and Federal statutes. Drainage Provinces. This white paper advances the concept of tile drainage provinces to aid in the discussion of regional differences in subsurface drainage and its overall effect on groundwater resources. Built upon the concept of groundwater provinces (Minnesota Department of Natural Resources, 2001), three distinct tile drainage provinces were conceptualized: (1) the Southeastern Province; (2) the SouthCentral Province; and, (3) the Western Province. The distinct geology and the soils that developed in these regions have implications for each region’s subsurface drainage density and the potential implications for groundwater. Knowledge gaps. Several critical knowledge gaps are identified in this paper, creating opportunities for further research to improve our understanding for better managed water resources: Extent of drainage is unknown. Direct estimates of the extent of subsurface drainage do not exist in Minnesota. However, several indirect methods have been utilized to estimate subsurface drainage, from the field-scale to countylevel through the use of geographic information system (GIS) analysis and aerial photography. Based on a 2012 U.S. Geological Survey estimate of subsurface drainage extent (Nakagaki and Wieczorek, 2016), about 21% of the land area in Minnesota has some density of subsurface drainage. Effect of drainage on underlying aquifers is unknown. A basic understanding of unconfined and confined aquifers and their recharge is necessary to connect any hydrological effects from agricultural drainage to groundwater. The basic goal of subsurface drainage to efficiently drain saturated soils clearly alters the water balance in croplands. However, its overall effect on groundwater resources has been poorly characterized, and is in large part determined by the geology below drained areas and the arrangement of underlying aquifers. Water balance shifts. An improved understanding of historical water balance shifts from pre- to post-drainage periods is necessary to understand long-term implications on net groundwater recharge. Also, more direct field-scale studies and indirect modeling studies are needed to characterize water budgets for fields with subsurface drainage.

Minnesota

Delineation of water sources for public-supply wells in three fractured-bedrock aquifer systems in Massachusetts

Fractured-bedrock aquifer systems in West Newbury, Maynard, and Paxton, Massachusetts, were studied to advance methods of data collection and analysis for delineating contributing areas to public-supply wells completed in fractured rock and for determining the effects of pumping on streams and wetlands. Contributing areas, as defined for this study, include all areas through which ground water flows from recharge areas to wells. In West Newbury, exploratory public-supply wells at two locations were completed in phyllite of the Eliot Formation. Aquifer testing indicated that subhorizontal and steeply dipping fractures that parallel two sets of foliation form elongated transmissive zones in the bedrock aquifer near the two well locations and also form a vertical hydraulic connection to surficial materials consisting of till at one location and marine clay at the other location. Recharge to bedrock is largely through a thin veneer of till over bedrock, but leakage through thick drumlin tills also recharges bedrock. Simulated contributing areas for the three supply wells pumped at a combined rate of 251 gallons per minute encompass about 1.3 square miles and extend to ground-water divides within most of a subbasin of the Artichoke River. Pumping likely would reduce streamflow in the Artichoke River subbasin by approximately the pumping rate. Pumping is likely to affect wetland areas underlain by till near the wells because of the vertical hydraulic connection to surficial materials. In Maynard, three exploratory public-supply wells were completed in coarse-grained schist of the Nashoba Formation. Aquifer testing indicated that a dense network of fractures in bedrock forms a laterally extensive transmissive zone in bedrock that is well connected vertically to surficial materials consisting of sandy till, lacustrine silts, sand and gravel, and wetland deposits. The simulated contributing area for the three supply wells pumped at a combined rate of 780 gallons per minute encompasses about 1.8 square miles of the Fort Pond Brook drainage area. Pumping likely would reduce streamflow in Fort Pond Brook by about the same amount as the pumping rate, and wetland-water levels within a 2,000-foot radius from the wells are likely to be lowered below the land surface by pumping. In Paxton, three existing supply wells are completed in granofels and schist of the Paxton and Littleton Formations. Aquifer testing demonstrated that a shallow bedrock well completed to a depth of 150 feet is closely connected hydraulically to overlying till. Two deep wells, however, receive much of their water from fractures at depths below 500 feet. Ground-water flow in bedrock appears to be mostly through parting fractures along a foliation set that dips gently (10 degrees) eastward. These parting fractures at depth are poorly connected vertically to shallow bedrock

Water-Resources Investigations Report

Hydrogeologic characterization, groundwater chemistry, and vulnerability assessment, Ute Mountain Ute Reservation, Colorado and Utah

The U.S. Geological Survey, in cooperation with the Ute Mountain Ute Tribe (UMUT), initiated a study in 2016 to increase understanding of the hydrogeology and chemistry of groundwater within select areas of the Ute Mountain Ute Reservation (UMUR) in Colorado and Utah, identify vulnerabilities to the system and other natural resources, and outline information needs to aid in the understanding and protection of groundwater resources. The results presented for this study can be used to support the UMUT’s goal of protecting their vital groundwater resources on the UMUR. Hydrogeologic conditions were characterized for the surficial aquifer contained in Quaternary-age unconsolidated surficial deposits and the Dakota aquifer contained in the Cretaceous-age Dakota Sandstone. In the surficial aquifer, median depth to water ranges from about 5.4 to 17.2 feet below land surface in the Farm and Ranch Enterprise area and 11 to 34 feet below land surface in the Towaoc area, and the water table slopes generally southwest or south. A map of depth to the top of the Dakota Sandstone was constructed from existing well data. Depths range from zero in outcrop areas to more than 3,000 feet below land surface on mesas in the southeastern part of the UMUR. Groundwater-chemistry data were collected by the UMUT from 13 springs and 31 wells from 1996 through 2017. Specific conductance was much lower for samples from springs than from wells; median values were 512 and 6,024 microsiemens per centimeter at 25 degrees Celsius, respectively. Spring samples were well oxygenated. A few well samples were anoxic (dissolved oxygen concentrations less than 0.5 milligrams per liter [mg/L]), indicating reducing conditions in the aquifer. About 75 percent of spring samples had fresh water (total dissolved solids concentrations less than 1,000 mg/L), and about 85 percent of well samples had brackish or highly saline water (total dissolved solids concentrations greater than 1,000 mg/L). Water type for springs on the Ute Mountains was calcium bicarbonate. Lower-altitude springs had a calcium-sulfate water type. Most well samples had sodium as the dominant cation, and sulfate, bicarbonate, and chloride as the dominant anions. Fluoride concentrations in about 45 percent of well samples were greater than an agricultural-use standard of 2 mg/L. Nitrate plus nitrite concentrations in most spring and well samples were less than about 1.6 mg/L per liter. Concentrations in samples from wells in the irrigated agricultural area were elevated; the maximum concentration was 78.5 mg/L. About one-half of the trace-element samples had concentrations that were less than laboratory reporting limits. Only aluminum, arsenic, and selenium in spring samples, and boron and selenium in well samples, were detected at concentrations greater than surface-water standards or water-quality standards for agricultural use of groundwater. Only three organic compounds, the pesticides alachlor and atrazine and the volatile organic compound di(2-ethylhexyl) phthalate, were detected in well samples. The Escherichia coli bacteria was detected in 47 and 23 percent of samples from wells and springs, respectively. The E. coli detections included samples from three culturally significant springs, which did not meet the UMUT cultural-use standard of total absence of E. coli. Tritium and carbon-14 were the primary environmental tracers used for interpreting groundwater ages for Lopez 2 Spring and five wells (AP–1, 5000 Block, Cottonwood Spring, Goodknight, and SE Toe). Water from the AP–1 well contained a mixture of pre- and post-1950s recharge. Tritium and carbon-14 recharge ages for Lopez 2 Spring (post-1950s in age), Goodknight and SE Toe wells (pre-1950s in age), and Cottonwood Spring well (primarily pre-1950s in age) are supported by helium-4 data. The helium-4 data for the 5000 Block well are inconsistent with the tritium and carbon-14 age of pre-1950s recharge because of interference caused by high methane concentrations in the water. Springs and surficial deposits are more vulnerable to contamination from anthropogenic chemicals than deeper bedrock wells. Bedrock aquifers are vulnerable in areas where the geologic formations containing the aquifers are exposed at the land surface. Groundwater in deep bedrock aquifers is likely thousands of years old and is not currently affected by present-day land uses. Both shallow and deep groundwater are vulnerable to naturally occurring salts and minerals, such as of total dissolved solids, major ions, nitrate, and trace elements. Effects of a changing climate on water resources and other ecological characteristics of the UMUR could include changes in evapotranspiration, a decrease in snowpack, decreased aquifer recharge and flow of springs, a decrease in soil moisture, and increased occurrence of wildfires and forest mortality. Of particular interest for the UMUT are possible effects of a changing climate on medicinal and culturally important plants and springs Several information needs were identified during this study that would aid in the understanding and protection of groundwater resources on the UMUR. These include well-completion information for bedrock wells, the collection of environmental tracer data at additional wells, the addition of methane and hydrocarbon analysis to well sampling plans, and the resampling of springs and wells that were last sampled in 2002 or earlier.

Colorado

Water quality and quantity and simulated surface-water and groundwater flow in the Laurel Hill Creek Basin, southwestern Pennsylvania, 1991–2007

Laurel Hill Creek is considered one of the most pristine waterways in southwestern Pennsylvania and has high recreational value as a high-quality cold-water fishery; however, the upper parts of the basin have documented water-quality impairments. Groundwater and surface water are withdrawn for public water supply and the basin has been identified as a Critical Water Planning Area (CWPA) under the State Water Plan. The U.S. Geological Survey, in cooperation with the Somerset County Conservation District, collected data and developed modeling tools to support the assessment of water-quality and water-quantity issues for a basin designated as a CWPA. Streams, springs, and groundwater wells were sampled for water quality in 2007. Streamflows were measured concurrent with water-quality sampling at main-stem sites on Laurel Hill Creek and tributaries in 2007. Stream temperatures were monitored continuously at five main-stem sites from 2007 to 2010. Water usage in the basin was summarized for 2003 and 2009 and a Water-Analysis Screening Tool (WAST) developed for the Pennsylvania State Water Plan was implemented to determine whether the water use in the basin exceeded the “safe yield” or “ the amount of water that can be withdrawn from a water resource over a period of time without impairing the long-term utility of a water resource .” A groundwater and surface-water flow (GSFLOW) model was developed for Laurel Hill Creek and calibrated to the measured daily streamflow from 1991 to 2007 for the streamflow-gaging station near the outlet of the basin at Ursina, Pa. The CWPA designation requires an assessment of current and future water use. The calibrated GSFLOW model can be used to assess the hydrologic effects of future changes in water use and land use in the basin. Analyses of samples collected for surface-water quality during base-flow conditions indicate that the highest nutrient concentrations in the main stem of Laurel Hill Creek were at sites in the northeastern part of the basin where agricultural activity is prominent. All of the total nitrogen (N) and a majority of the total phosphorus (P) concentrations in the main stem exceeded regional nutrient criteria levels of 0.31 and 0.01 milligrams per liter (mg/L), respectively. The highest total N and total P concentrations in the main stem were 1.42 and 0.06 mg/L, respectively. Tributary sites with the highest nutrient concentrations are in subbasins where treated wastewater is discharged, such as Kooser Run and Lost Creek. The highest total N and total P concentrations in subbasins were 3.45 and 0.11 mg/L, respectively. Dissolved chloride and sodium concentrations were highest in the upper part of the basin downstream from Interstate 76 because of road deicing salts. The mean base-flow concentrations of dissolved chloride and sodium were 117 and 77 mg/L, respectively, in samples from the main stem just below Interstate 76, and the mean concentrations in Clear Run were 210 and 118 mg/L, compared to concentrations less than 15 mg/L in tributaries that were not affected by highway runoff. Water quality in forested tributary subbasins underlain by the Allegheny and Pottsville Formations was influenced by acidic precipitation and, to a lesser extent, the underlying geology as indicated by pH values less than 5.0 and corresponding specific conductance ranging from 26 to 288 microsiemens per centimeter at 25 degrees Celsius for some samples; in contrast, pH values for main stem sites ranged from 6.6 to 8.5. Manganese (Mn) was the only dissolved constituent in the surface-water samples that exceeded the secondary maximum contaminant level (SMCL). More than one-half the samples from the main stem had Mn concentrations exceeding the SMCL level of 50 micrograms per liter (μg/L), whereas only 19 percent of samples from tributaries exceeded the SMCL for Mn. Stream temperatures along the main stem of Laurel Hill Creek became higher moving downstream. During the summer months of June through August, the daily mean temperatures at the five sites exceeded the limit of 18.9 degrees Celsius (°C) for a cold-water fishery. The maximum instantaneous values for each site ranged from 27.2 to 32.8 °C. Water-quality samples collected at groundwater sites (wells and springs) indicate that wells developed within the Mauch Chunk Formation had the best water quality, whereas wells developed within the Allegheny and Pottsville Formations yielded the poorest water quality. Waters from the Mauch Chunk Formation had the highest median pH (7.6) and alkalinity (80 mg/L calcium carbonate) values. The lowest pH and alkalinity median values were in waters from the Allegheny and Pottsville Formations. Groundwater samples collected from wells in the Allegheny and Pottsville Formations also had the highest concentrations of dissolved iron (Fe) and dissolved Mn. Seventy-eight percent of the groundwater samples collected from the Allegheny Formation exceeded the SMCL of 300 μg/L for Fe and 50 μg/L for Mn. Forty-three and 62 percent of the groundwater samples collected from the Pottsville Formation exceeded the SMCL for iron and Mn, respectively. The highest Fe and Mn concentrations for surface waters were measured for tributaries draining the Pottsville Formation. The highest median Fe concentration for tributaries was in samples from streams draining the Allegheny Formation. During base-flow conditions, the streamflow per unit area along the main stem of Laurel Hill Creek was lowest in the upper parts of the basin [farthest upstream site 0.07 cubic foot per second per square mile (ft 3 /s/mi 2 )] and highest (two sites averaging about 0.20 (ft 3 /s/mi 2 ) immediately downstream from Laurel Hill Lake in the center of the basin. Tributaries with the highest streamflow per unit area were those subbasins that drain the western ridge of the Laurel Hill Creek Basin. The mean streamflow per unit area for tributaries draining areas that extend into the western ridge and draining eastern or central sections was 0.24 and 0.05 ft 3 /s/mi 2 , respectively. In general, as the drainage area increased for tributary basins, the streamflow per unit area increased. Criteria established by the Pennsylvania Department of Environmental Protection indicate that the safe yield of water withdrawals from the Laurel Hill Creek Basin is 1.43 million gallons per day (Mgal/d). Water-use data for 2009 indicate that net (water withdrawals subtracted by water discharges) water withdrawals from groundwater and surface-water sources in the basin were approximately 1.93 Mgal/d. Water withdrawals were concentrated in the upper part of the basin with approximately 80 percent of the withdrawals occurring in the upper 36 mi 2 of the basin. Three subbasins—Allen Creek, Kooser Run, and Shafer Run— in the upper part were affected the most by water withdrawals such that safe yields were exceeded by more than 1,000 percent in the first two and more than 500 percent in the other. In the subbasin of Shafer Run, intermittent streamflow characterizes sections that historically have been perennial. The GSFLOW model of the Laurel Hill Creek Basin is a simple one-layer representation of the groundwater flow system. The GSFLOW model was primarily calibrated to reduce the error term associated with base-flow periods. The total amount of observed streamflow at the Laurel Hill Creek at Ursina, Pa. streamflow-gaging station and the simulated streamflow were within 0.1 percent over the entire modeled period; however, annual differences between simulated and observed streamflow showed a range of -27 to 24 percent from 1992 to 2007 with nine of the years having less than a 10-percent difference. The primary source of simulated streamflow in the GSFLOW model was the subsurface (interflow; 62 percent), followed by groundwater (25 percent) and surface runoff (13 percent). Most of the simulated subsurface flow that reached the stream was in the form of slow flow as opposed to preferential (fast) interflow.

Pennsylvania

Genetic integrity, population status, and long-term viability of isolated populations of shoal bass in the upper Chattahoochee River basin, Georgia

This report characterizes the status of multiple isolated Shoal Bass ( Micropterus cataractae ) populations in the upper Chattahoochee River basin (UCRB), Georgia. The Shoal Bass, a sport fish endemic to the Apalachicola-Chattahoochee-Flint River (ACF) basin, is a fluvial-specialist species considered vulnerable to local extirpations and extinction due to habitat fragmentation and introgression with non-native congeners. Perhaps one of the most isolated populations of Shoal Bass exists in a 2-km reach of Big Creek, a tributary of the Chattahoochee River located near Roswell, Georgia. Big Creek is partially contained within the Chattahoochee River National Recreation Area, although the Big Creek watershed is riddled with urban land cover. Roswell Mill Dam limits the upstream extent of the Shoal Bass population at Big Creek, and the downstream extent is presumably limited to the confluence of Big Creek and the Chattahoochee River. This reach of the Chattahoochee River is thermally depressed because of coldwater releases from Lake Lanier, and is considered unsuitable for Shoal Bass. Herein, we examine the genetic integrity, population status, and long-term viability of the Shoal Bass population in Big Creek. We also examine two additional Shoal Bass populations that occur in the UCRB, specifically the Chestatee River and the upper Chattahoochee River, both of which are impounded at Lake Lanier. Together, the Shoal Bass inhabiting these three stream systems comprise a distinct genetic stock of Shoal Bass (Taylor 2017), underscoring the importance of conserving these populations towards maintaining the overall diversity and adaptive potential of the species. We assessed genetic diversity and estimated effective population sizes within these three rivers by genotyping fish with 16 microsatellite DNA markers. Results demonstrated that the Shoal Bass population in Big Creek has experienced high rates of introgression with non-native Smallmouth Bass (M. dolomieu), purportedly introduced into the Chattahoochee River in the past 10-15 years. Alarmingly, only 24% (15 of 62) of putative Shoal Bass collected from Big Creek were genetically pure Shoal Bass, whereas the majority of fish were first-filial (F1) generation hybrids and unidirectional backcrosses towards Shoal Bass. Fleeting opportunity may remain to conserve the native genome of the Shoal Bass population in Big Creek. High hybridization rates prevented genetic diversity analysis for the Big Creek population. Shoal Bass populations in the Chestatee and Chattahoochee rivers displayed levels of genetic diversity similar to populations that persist in other rivers in the ACF basin, namely the Flint and Chipola rivers. Effective population sizes of 93.8– 197.4 for the Chestatee and Chattahoochee rivers (combined) suggest that the conservation status of these populations is stable for the short-term, but may be at risk of losing genetic diversity and adaptive potential in the long-term. To estimate age and mortality of the three populations, we used fish scales and capture-markrecapture (CMR) as complementary, non-lethal methods for age estimation. Estimated ages of phenotypic Shoal Bass ranged from 1-12 years in all three populations, demonstrating increased longevity compared to populations elsewhere within the native range. Catch-curve estimates of annual mortality ranged from 18.4-23.7%, which are markedly lower than those observed in other Shoal Bass populations in the ACF basin. These differences in life-history characteristics underscore the need for the development of population-specific management and conservation strategies for Shoal Bass in the UCRB. The lowest recruitment variability (i.e., the variation in year-class strength) was observed in the Chestatee River, a forested watershed, whereas the highest variability was observed in Big Creek, an urbanized watershed. Recruitment strength in Big Creek was negatively influenced by discharge variability in the summer months, suggesting that flashy, sediment-laden flows hinder survival of recently hatched young. Other statistically significant models from Big Creek and the Chattahoochee River indicated that over-winter survival could be an important pinch-point for recruitment in UCRB populations. A multi-agency sampling effort was conducted from May 2013-May 2016 to estimate the population size of Shoal Bass occupying the 1-km of wadeable shoal habitats in Big Creek. Using CMR models, we estimated that approximately 219-348 Shoal Bass (≥ 70 mm total length) occupied the area throughout the duration of our study. These estimates largely reflect abundance of individuals aged 0-2 years, as only 9% (36 of 408) tagged fish were aged ≥ 7 years. Local abundance appeared similar to that reported for a population that inhabited Little Uchee Creek, a similar-sized tributary of the Chattahoochee River, prior to its recent functional extirpation. The low abundance of large, adult Shoal Bass further suggests the long-term viability of the Big Creek population may be in jeopardy. Perhaps most importantly, CMR estimates reflect abundance of phenotypic Shoal Bass – genetic analyses suggest the abundance of pure Shoal Bass could be an order of magnitude smaller. To evaluate the potential for adult Shoal Bass to emigrate from Big Creek into the mainstem Chattahoochee River, we tagged eight adults with acoustic telemetry tags and assessed their seasonal residency at two stationary receiver locations located in increasing proximity to the confluence with the Chattahoochee River. Fish took up residency near the confluence during the fall and winter months, during which time water temperatures in Big Creek were periodically colder than the Chattahoochee River. Although we were unable to document emigration, we conclude that the potential for emigration is highest during the winter months when the Chattahoochee River may be warmer than Big Creek. Two of the tagged fish were caught by anglers near the confluence, suggesting that angling pressure at Big Creek may be higher than previously suspected. Overall, this study observed unique life-history characteristics and characterized the population status of multiple Shoal Bass populations in the UCRB. Populations in the Chestatee and Chattahoochee rivers appear stable at present and likely represent the last remaining strongholds for pure Shoal Bass in the UCRB. Efforts to preserve forested watershed conditions, natural hydrology, and shoal habitats would contribute to the long-term persistence of Shoal Bass populations in these two rivers. Additionally, the detection of non-native Alabama Bass and their associated hybrids in both rivers is cause for concern. Diligent monitoring of hybridization dynamics between Alabama Bass and Shoal Bass is warranted, along with an assessment of Alabama Bass invasion extent upstream of Lake Lanier. The Shoal Bass population in Big Creek is threatened by elevated levels of introgression with nonnative Smallmouth Bass, recruitment variability, low abundance of adults, and isolation from other populations. Conservation intervention is urgently needed to restore and preserve this genetically distinct population, which would contribute to preservation of range wide genetic diversity and adaptability of the species. Additionally, an urban sport fishery for Shoal Bass at Big Creek has the potential to serve as a tool for increasing public awareness, engagement, and support of Shoal Bass conservation efforts in the UCRB. We suggest strategies for conservation of the remnant shoal habitats and Shoal Bass population in Big Creek, including potential development of a supplemental stocking program, selective removal of non-native congeners, and delivery of environmental education programs that could bolster awareness and appreciation.

Georgia

Spatial and temporal variation of stream chemistry associated with contrasting geology and land-use patterns in the Chesapeake Bay watershed—Summary of results from Smith Creek, Virginia; Upper Chester River, Maryland; Conewago Creek, Pennsylvania; and Difficult Run, Virginia, 2010–2013

Despite widespread and ongoing implementation of conservation practices throughout the Chesapeake Bay watershed, water quality continues to be degraded by excess sediment and nutrient inputs. While the Chesapeake Bay Program has developed and maintains a large-scale and long-term monitoring network to detect improvements in water quality throughout the watershed, fewer resources have been allocated for monitoring smaller watersheds, even though water-quality improvements that may result from the implementation of conservation practices are likely to be first detected at smaller watershed scales. In 2010, the U.S. Geological Survey partnered with the U.S. Environmental Protection Agency and the U.S. Department of Agriculture to initiate water-quality monitoring in four selected small watersheds that were targeted for increased implementation of conservation practices. Smith Creek watershed is an agricultural watershed in the Shenandoah Valley of Virginia that is dominated by cattle and poultry production, and the Upper Chester River watershed is an agricultural watershed on the Eastern Shore of Maryland that is dominated by row-cropping activities. The Conewago Creek watershed is an agricultural watershed in southeastern Pennsylvania that is characterized by mixed agricultural activities. The fourth watershed, Difficult Run, is a suburban watershed in northern Virginia that is dominated by medium density residential development. The objective of this study was to investigate spatial and temporal variations in water chemistry and suspended sediment in these four relatively small watersheds that represent a range of land-use patterns and underlying geology to (1) characterize current water-quality conditions in these watersheds, and (2) identify the dominant sources, sinks, and transport processes in each watershed. The general study design involved two components. The first included intensive routine water-quality monitoring at an existing streamgage within each study area (including continuous water-quality monitoring as well as discrete water-quality sampling) to develop a detailed understanding of the temporal and hydrologic variability in stream chemistry and sediment transport in each watershed. The second component involved extensive water-quality monitoring at various sites throughout each watershed to develop a detailed understanding of spatial patterns. Both components were used to improve understanding of sources and transport processes affecting stream chemistry, including nutrients and suspended sediments, and their implications for detecting long-term trends related to best management practices. This report summarizes the results of monitoring that was performed from April 2010 through September 2013. Individual Small Watershed Summaries Summaries for each of the four small watersheds are presented below. Each watershed has a more descriptive and detailed section in the report, but these summaries may be particularly useful for some watershed managers and stakeholders desiring slightly less technical detail. Smith Creek Smith Creek is a 105.39-mi 2 watershed within the Shenandoah Valley that drains to the North Fork Shenandoah River. The long-term Smith Creek base-flow index is 72.3 percent, indicating that on average, approximately 72 percent of Smith Creek flow was base flow, which suggests that Smith Creek streamflow is dominated by groundwater discharge rather than stormwater runoff. A series of cluster and principal components analyses demonstrated that the majority of the variability in Smith Creek water quality could be attributed to hydrologic and seasonal variability. Statistically significant positive correlations with flow were observed for turbidity, suspended sediments, total nitrogen, ammonium, orthophosphate, iron, total phosphorus, and the ratio of calcium to magnesium. Statistically significant inverse correlations with flow were observed for specific conductance, magnesium, δ 15 N of nitrate, pH, bicarbonate, calcium, and δ 18 O of nitrate. Of particular note, flow and nitrate were not statistically significantly correlated, likely because of the relatively complex concentration-discharge relationship observed in continuous and discrete datasets. Statistically significant seasonal patterns were observed for numerous water-quality constituents: water temperature, turbidity, orthophosphate, total phosphorus, suspended-sediment concentration, and silica were higher during the warm season, but pH, dissolved oxygen, and sulfate were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient loads in Smith Creek using the continuous water-quality monitors. The mean Smith Creek in-stream sediment load was approximately 6,900 tons per year, with nearly 90 percent of the sediment load over the 3-year study period contributed during the eight largest storm events during that period. The Smith Creek total phosphorus load was approximately 21,000 pounds of phosphorus per year, with the majority of the load contributed during stormflow periods, although a substantial phosphorus load still occurs during base-flow conditions. The Smith Creek total nitrogen load was approximately 400,000 pounds per year, with total nitrogen accumulation less dominated by stormflow contributions (as was the case for sediment and total phosphorus) and strongly affected by base-flow export of nitrogen from the basin. Extensive water-quality monitoring throughout the Smith Creek watershed revealed how the complex geology and hydrology interacted to result in variable water chemistry. During relatively dry and low base-flow periods, much of the discharge in Smith Creek was contributed by a single dominant spring—Lacey Spring. During wetter base-flow periods, the flows in Smith Creek were largely generated by a mixture of headwater springs and forested mountain tributaries with very different geochemical composition. The headwater springs generally issued from limestone bedrock and were characterized as having relatively high nitrate, specific conductance, calcium, and magnesium, as well as relatively low concentrations of phosphorus, ammonium, iron, and manganese. The undeveloped, high-gradient, forested mountain sites were generally characterized by low ionic strength waters with low nutrient concentrations. Nitrate isotope data from the limestone springs generally were consistent with manure-derived nitrogen sources (such as cattle and poultry), although the possibility of other mixed sources cannot be excluded. Nitrate isotope data from the undeveloped, high-gradient forested mountain sites were more consistent with nitrogen from undisturbed soils, atmospheric deposition, or nitrogen fixation. Regardless of the nitrogen source, oxygen isotope data indicate that the nitrate was largely a result of nitrification. Land-use data indicate that manure sources of nitrogen dominated watershed nitrogen inputs. Phosphorus sources were less well studied. The presence of a single point-source discharge near the town of New Market contributed the majority of the phosphorus to Smith Creek under base-flow conditions, but nonpoint sources of phosphorus dominated the loading to Smith Creek during stormflow periods. Implementation of conservation practices increased in the Smith Creek watershed during the study period, and even though a broad range of practice types was implemented, the most common practices included stream fencing (for cattle exclusion), the development of nutrient management plans, conservation crop rotation, and the planting of cover crops. While the implementation of these conservation practices is encouraging, results indicate small increases in nitrate concentrations at the streamgage over the last 29 years, concurrent with small decreases in nitrate fluxes. It will likely be years before the cumulative effect of these practices can be detected in the Smith Creek water quality, and the magnitude of the effect of these conservation practices detected in Smith Creek will depend largely on whether nutrient loading (of manure and commercial fertilizer) is reduced over time. Upper Chester River The Upper Chester River watershed includes the 36-square-mile (mi 2 ) watershed area around several nontidal tributaries that drain into the tidal Chester River. The streamgage is on Chesterville Branch, the largest nontidal tributary (approximately 6.12 mi 2 ) and is the site for continuous water-quality monitoring for this project. The base-flow index at Chesterville Branch is about 72 percent and indicates that, as in most of the Coastal Plain, groundwater is the greatest contributor to streamflow. As such, more than 90 percent of the nitrogen in the stream is in the form of nitrate from groundwater. Continuous and discrete data collected at Chesterville Branch show the effects of streamflow and season on water quality. Significantly positive correlations with flow were observed for ammonium, dissolved and total phosphorus, sediment, and turbidity as runoff carried these constituents from the land surface into Chesterville Branch. Other constituents that increased significantly with flow include potassium, sulfate, iron, and manganese, which are likely contributed from near-stream areas and ponds with high organic-matter content. Total nitrogen, pH, and specific conductance, along with chemical constituents associated with groundwater inputs including nitrate, calcium, ratio of calcium to magnesium, silica, bicarbonate, and sodium, were negatively correlated with flow because concentrations of these constituents were diluted by runoff. Seasonal differences in water chemistry, which are most likely related to increased biologic effects on the uptake and release of chemicals in the stream and near-stream areas, also were observed. Water temperature, orthophosphate, δ 15 N of nitrate, bicarbonate, sodium, and the ratio of sodium to chloride were higher during the warm season, and dissolved oxygen, total nitrogen, nitrate, magnesium, sulfate, and manganese were higher during the cool season. Surrogate-regression models developed by using continuous water-quality data showed that the annual sediment load for the 2013 water year was about 2,600 tons, with more than 90 percent of this sediment contributed during two storms. The total phosphorus load in 2013 was about 13,000 pounds with more than 90 percent contributed during the same two storms as sediment. The load of total nitrogen, 140,000 pounds, accumulated steadily throughout the 2013 water year as nitrate in groundwater continuously discharged into the stream. The same two large storms that contributed 90 percent of the suspended-sediment and total phosphorus load only contributed about 20 percent of the annual total nitrogen load. Extensive water-quality monitoring of stream base flow throughout the Upper Chester River watershed identified how differences in land use and hydrogeology affected water chemistry. In parts of the watershed with well-drained soil and thick sandy aquifer sediments, concentrations of nitrate and other chemicals associated with fertilizer and lime application increased in streams as agricultural land use increased. More than 90 percent of the nitrogen in streams from these areas was in the form of nitrate, and concentrations ranged from about 5 milligrams per liter (mg/L) to 8 mg/L as nitrogen in the two largest tributaries. Stream nitrate concentrations were about 1 mg/L as nitrogen where soils were more poorly drained, the surficial aquifer sediments were thinner, and forests and wetlands were more widespread than agriculture. Nitrate isotope data were consistent with inorganic fertilizers ± atmospheric deposition and N 2 fixation as sources of nitrogen, and with nitrification as the dominant nitrate-forming process. Nitrate reduction was indicated by elevated δ 15 N and δ 18 O values in some samples from streams draining watersheds with poorly drained soils. An analysis of land-use data and SPARROW modeling input data attributed almost 90 percent of the nitrogen sources in the Upper Chester River watershed to inorganic fertilizer and fixation of atmospheric nitrogen by legumes, which is in agreement with the isotopic characteristics of nitrate in this watershed. Local sources of manure are limited in this area. Total phosphorus concentrations during base flow ranged from below detection to about 0.2 mg/L. Stream phosphorus concentrations during base flow were generally lower than those measured during storms because most phosphorus transport likely occurs as phosphorus attached to sediment particles during runoff. Because manure is not widely used in this area, the major source of phosphorus is likely fertilizer. The implementation of conservation practices in the Upper Chester River watershed increased substantially during the study period, with a total implementation of 1,194 U.S. Department of Agriculture-compliant practices. The most frequently used practices were oriented towards nutrient and sediment control, including cover crops, nutrient management planning, conservation crop rotation, conservation tillage, and irrigation management. The current Chesapeake Bay model for this area predicts that implementation of best management practices should result in a 13-percent decrease in overall delivery of nitrogen to the Upper Chester River. Because most nitrogen travels through the groundwater system for years to decades before being discharged to streams, the time period of monitoring was not sufficient to see the effects of these practices on water quality. The magnitude of the effect that may eventually be detected will depend on the degree to which nitrate leaching into the groundwater system is reduced over time. Loadings of phosphorus and sediment are primarily transported during large runoff events and are difficult to control and analyze for trends because of their timing and episodic nature. Conewago Creek Conewago Creek has two primary monitoring locations—one near the middle of the 47-mi 2 watershed and the other near the outlet just upstream of the Susquehanna River. The base-flow index was 47.3 percent for 2012–2013, indicating that on average, approximately 53 percent of the streamflow in Conewago Creek exited the watershed as surface flow, which suggests that the stormwater runoff was somewhat greater than groundwater discharge (base flow). A series of cluster and principal components analyses demonstrated that the majority of the variability in the Conewago Creek water quality could be attributed to hydrologic and seasonal variability. Statistically significant positive correlations with flow were observed at both monitoring sites for ammonium, total phosphorus, orthophosphate, iron, and manganese; additionally, at the upstream monitoring station, total nitrogen demonstrated a statistically significant positive correlation with flow. Statistically significant inverse correlations with flow were observed at both sites for water temperature, specific conductance (at the downstream site only), sulfate, chloride, calcium, and magnesium. Statistically significant seasonal patterns were observed for several water-quality constituents. Water temperature, phosphorus (upstream site only), and orthophosphate were higher during the warm season, and nitrate and total nitrogen (upstream site only) were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient load in Conewago Creek by using the continuous water-quality monitors and water-quality samples. Conewago Creek sediment load was approximately 9,900 tons in 2012 and approximately 18,900 tons in 2013, with nearly 80 percent of the sediment load in 2013 contributed by the three largest storm events. Annual total nitrogen loads could not be estimated due to poor model performance. The addition of continued monitoring or a continuously recording nitrate sensor could improve estimates of total nitrogen loads. During 2012 and 2013, phosphorus loads in Conewago Creek were approximately 50,000 pounds in each year. Combining data from one high-flow synoptic sampling with the data from routine sampling revealed how the geology and hydrology interact to result in variable water chemistry throughout the Conewago Creek watershed. The areas above the upstream gage in the headwaters are generally underlain by forested non-carbonate bedrock and are characterized by relatively low nitrate, specific conductance, calcium, and magnesium, as well as relatively low concentrations of phosphorus, ammonium, iron, and manganese. The more developed, agricultural areas below the upstream site were generally characterized by higher ionic strength waters with higher nutrient and metal concentrations. An analysis of land-use data and SPAtially Referenced Regressions On Watershed (SPARROW) modeling data indicates that manure sources of nitrogen dominate the input of nitrogen to the watershed. Implementation of conservation practices increased in the Conewago Creek watershed during the study period, and while a broad range of practice types were implemented, the most common practices included residue and tillage management, cover crops, nutrient management, terracing, and stream fencing (for animal exclusion or bank restoration). While the implementation of these conservation practices is encouraging, the cumulative effects of these practices probably will not be detected in Conewago Creek water quality for several years. The magnitude of the effects of these conservation practices on water quality in Conewago Creek will depend largely on the extent to which nutrient loading (septic, manure, and commercial fertilizer) and sediment-producing activities are reduced over time. Difficult Run The Difficult Run watershed is a 57.82-mi 2 watershed that drains to the Potomac River. The long-term Difficult Run base-flow index (from 1936 to 2010) was 57.9, indicating that approximately 58 percent of streamflow exited the watershed as base flow and 42 percent as stormflow; however, with continued development and urbanization of the watershed, the base-flow index has decreased to 50 percent during the last 20 years. This base-flow index was less than those of the other watersheds evaluated in this study, likely because the Difficult Run watershed largely is underlain by crystalline piedmont metamorphic rocks and has a greater proportion of impervious urban land cover. A series of cluster and principal components analyses indicated that most of the variability in Difficult Run water quality could be attributed to hydrologic variability and seasonality. Statistically significant positive correlations with flow were observed for turbidity, dissolved oxygen, suspended sediments, ammonium, orthophosphate, iron, and total phosphorus. Statistically significant inverse correlations with flow were observed for water temperature, pH, specific conductance, bicarbonate, calcium, magnesium, nitrate, δ 15 N of nitrate, and silica. Statistically significant seasonal patterns were observed for numerous water-quality constituents: water temperature, ammonium, orthophosphate, and δ 15 N of nitrate were higher during the warm season, and dissolved oxygen, nitrate, and manganese were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient loading rates. The Difficult Run sediment load was approximately 8,000 tons per year, with greater than 95 percent of the sediment load in the 2013 water year contributed by the seven largest storm events. The total phosphorus load in Difficult Run was approximately 14,000 pounds of phosphorus per year, with the majority of the load contributed during stormflow periods. The total nitrogen load in Difficult Run is estimated to have been approximately 140,000 pounds per year, with total nitrogen accumulation less dominated by stormflow contributions than that of phosphorus and strongly affected by base-flow export of nitrogen from the basin. Extensive water-quality monitoring throughout the Difficult Run watershed revealed relatively uniform generation of flow per unit of watershed area, as well as spatial variation in water quality that is strongly related to land-use activities. Elevated nitrate concentrations were observed in a subset of monitoring sites that are inversely correlated with population density and positively correlated to the septic system density within each subwatershed. The majority of the elevated nitrate concentrations for these sites are hypothesized to be caused by nitrate leaching from septic systems, more so than homeowner fertilizer usage among these subwatersheds that have lower population densities than other parts of the watershed. Nitrate isotope data, temporal patterns in the water-quality data, mass-balance computations, and a separate land-use analysis all generally indicate that leachate from septic systems was the likely source of the elevated nitrate. Another group of water-quality sites have relatively low nitrogen concentrations, are located in areas that are served by city sewer lines, and have experienced stream restoration activities. A final group of sites drained the areas with the highest imperviousness and had strongly elevated specific conductance, chloride, and sodium, which were likely caused by a combination of road salting and other anthropogenic sources draining these urbanized areas in the watershed. A fourth group of sites represents a mixture of water sources and had water quality similar to that at the Difficult Run streamgage. Analysis of the nitrate isotope data generally indicates a broad range of composition indicative of mixed natural and anthropogenic nitrogen sources. Implementation of conservation practices increased in the Difficult Run watershed during the study period, and while a broad range of practice types was implemented, the most common practices included stream restoration. While the implementation of these conservation practices is encouraging, the cumulative effect of these practices probably will not be detected in Difficult Run water quality for several years.

Maryland, Pennsylvania, Virginia