USGS Science⌕ Search

SEARCH · USGS Science

Results for “Final Project Report”

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

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

388 records · Page 22Linked to original sources
USGS2020

Simulation of discharge, water-surface elevations, and water temperatures for the St. Louis River estuary, Minnesota-Wisconsin, 2016–17

The St. Louis River estuary is a large freshwater estuary, next to Duluth, Minnesota, that encompasses the headwaters of Lake Superior. The St. Louis River estuary is one of the most complex and compromised near-shore systems in the upper Great Lakes with a long history of environmental contamination caused by logging, mining, paper mills, and other heavy industrial activities. Presently (2020), a widely available, science-based assessment tool capable of evaluating ecosystem-level responses to remediation and restoration projects has not existed for the estuary. To address this need, the U.S. Geological Survey (USGS) built a predictive, mechanistic, three-dimensional hydrodynamic model for the estuary using the Environmental Fluid Dynamics Code framework. In the current version, the model can simulate continuous discharge, water-surface elevations, water temperature, and flow velocity, although the modular framework allows for future additions of water-quality modeling. The model was calibrated using data collected from April 2016 through November 2016 and validated with data collected from April 2017 through November 2017. The four types of data used to evaluate model performance were water-surface elevations, discharge, water temperature, and flow velocities. Streamflow and temperature boundary condition data included a mixture of USGS streamgage data, Minnesota Department of Natural Resources gage data, and estimates derived from the gage data. The model was able to simulate the water-surface elevations with generally good agreement between the simulated and measured values for both years at the daily time step. Specifically, the model was able to demonstrate excellent agreement with the measured data with Nash-Sutcliffe efficiency coefficients greater than 0.8 for all three locations; however, the model was unable to produce hourly water-surface elevations with such accuracy for 2016–17. Discharge was more dynamic than the water-surface elevations, both for the measured and simulated data. Generally, most of the discharge ranged from −650 to 1,200 cubic meters per second, but the constantly changing flux exiting the estuary into Lake Superior (positive flows) and entering the estuary from Lake Superior (negative flows) occurred throughout the year. Even upstream at the St. Louis River at Oliver, Wisconsin, gage (USGS station 0402403250), the effect of flows into the estuary from Lake Superior did occur, demonstrating the strong effect of the Lake Superior seiche on flows for the estuary. From a performance standpoint, the model was able to simulate discharge with generally good agreement in both years, although the 2017 validation was better than the 2016 calibration period. For the daily Nash-Sutcliffe efficiency coefficients, the simulated values were 0.98, 0.62, 0.49, and 0.71 for the Oliver gage; the Superior Bay entry channel at Superior, Wisc., (USGS station 464226092005600); the Superior Bay Duluth Ship Canal at Duluth, Minn., (USGS station 464646092052900); and total entries (combination of the Superior entry and Duluth entry), respectively. For the hourly evaluation criteria, the model performed poorly, with Nash-Sutcliffe efficiency coefficients less than 0 for the two entries into Lake Superior; therefore, as a predictor of discharge at the hourly scale, the model performed worse than using the measured data average. Similar to discharge, the model was a good predictor of flow velocity at the daily time scale but had difficulty matching the measured data at the hourly scale. For discharge and flow velocity, matching at subdaily time steps for a system as complicated as the St. Louis River estuary is considered difficult because the match is highly sensitive to coordinating the exact measurement location to the simulated value. The final calibration target was water temperature, calibrated for the Oliver gage and the Duluth entry. For calibration purposes, the Duluth entry was the more important water temperature target because the Oliver gage was more of an internal check on the model. The Nash-Sutcliffe efficiency coefficients for the Duluth entry were high; hourly Nash-Sutcliffe efficiency coefficients at the Duluth entry were either at or greater than 0.7 for both years, and daily values were 0.84 and 0.82 for 2016 and 2017, respectively.

Minnesota, Wisconsin↗
USGS2005

Upper Auglaize watershed AGNPS modeling project

The Upper Auglaize Watershed agricultural non-point source modeling project was an interagency effort to use a Geographic Information System (GIS)-based modeling approach for assessing and reducing pollution from agricultural runoff and other non-point sources. This project applied the U.S. Department of Agriculture (USDA), Agricultural Research Service’s AGricultural Non-Point Source (AGNPS) suite of models to the Upper Auglaize River Watershed, a major watershed within the Maumee River Basin. This modeling project was conducted by an interagency team consisting of a partnership between the: (1) USDA, Agricultural Research Service (ARS); (2) USDA, Natural Resources Conservation Service (NRCS); (3) U.S. Army Corps of Engineers (USACE); (4) U.S. Geological Survey (USGS); (5) Ohio State University; (6) University of Toledo (UT); (7) Heidelberg College; (8) Ohio Department of Natural Resources (ODNR), Division of Soil and Water Conservation; (9) Ohio Environmental Protection Agency (OEPA); and (10) Allen, Auglaize, Van Wert, and Putnam Soil and Water Conservation Districts. The partnership was the first step in a process to eventually apply the model in a portioned subset of watersheds for the Maumee Basin, and then to link them to form a comprehensive basin-wide model This work was performed under the authority of Section 516(e) of the Water Resources Development Act (WRDA) of 1996, as amended, for the purpose of assisting State and local watershed managers with their evaluation, prioritization and implementation of alternatives for soil conservation, sediment trapping and non-point source pollution prevention in the Upper Auglaize River watershed. The project team, working in a cooperative effort, used the models to determine sediment sources, contributing locations, and the effect of application of best management practices (BMPs) on rates of sediment delivery to the mouth of the watershed. The results will be used to guide conservation incentive and land treatment programs. The team relied heavily on Geographic Information System (GIS)-based applications to expedite the application of the model. The results of the analysis demonstrated that the application of BMPs would have a positive effect on reducing the loadings of sediment leaving the mouth of the Upper Auglaize Watershed. An application of 17 percent new no-till acres and eight percent new grassland acres, when randomly applied to the watershed, reduced loadings at the mouth to 82 percent of the simulated existing condition loadings. No-till, conversion of cropland to grassland, other uses including grass buffers, and reforestation of parts of the watershed, were all shown by the model to have a measurable effect on reducing sediment loads. Conversion of all of the cropland in the watershed to no-till would reduce the average unit load (tons of sediment per acre) leaving the mouth of the watershed to a level that is 42 percent of the simulated existing condition load. Ephemeral gullies were found to be the primary source of erosion (72 percent), sediment yield (73 percent), and sediment loading (73 percent). Controlling sediment load means controlling gully erosion and possibly trapping sediment yield before it reaches the stream system. Most BMPs (e.g., no-till, conversion of cropland, etc.) that reduce sheet and rill erosion and its sediment yield will also reduce gully erosion and its sediment yield. However, grassed waterways, which have no effect on sheet and rill erosion, are frequently an effective BMP to prevent ephemeral gullies. And, of course, riparian vegetation and sediment traps would reduce the delivery ratios of all types of landscape erosion. New techniques were developed by the team to quantify the ephemeral gully erosion within the model. When calibrated to available stream gage data the model suggests that more (73% in the existing condition simulation) of the sediment load originates from ephemeral gully erosion than from traditional sheet and rill erosion. The model quantified the value of tile drainage in reducing the sediment load from the watershed. Loadings under drained conditions were always less than loadings under undrained conditions for otherwise identical land uses. The average sediment load of all alternatives for drained loadings was 89.2 percent of the load for the corresponding undrained loadings. The model established that while many conservation incentive programs treat tile drainage as a production practice, significant erosion and sediment control benefits are provided by the practice in comparison to cultivation in an undrained state.

Ohio↗
USGS2016

Accuracy assessment/validation methodology and results of 2010–11 land-cover/land-use data for Pools 13, 26, La Grange, and Open River South, Upper Mississippi River System

Introduction/Background The U.S. Geological Survey (USGS)-Upper Midwest Environmental Sciences Center (UMESC) was responsible for development of several land cover/land use (LCU) systemic datasets of the Upper Mississippi River System (UMRS). These efforts (1989 and 2000) were funded by the U.S. Army Corps of Engineers’ Upper Mississippi River Restoration Program (UMRR) Long Term Resource Monitoring (LTRM) element. Development of systemic datasets includes the acquisition, processing, and serving of high-resolution aerial photography and land cover/land use spatial datasets ( http://www.umesc.usgs.gov/data_library/land_cover_use/land_cover_use_data.html ). In 2008, the UMRR reached a collaborative agreement with the U.S. Fish and Wildlife Service-Region 3 to collect high-resolution digital imagery of the entire UMRS floodplain during 2010–11 for LTRM. The UMESC helped acquire, process, and serve this imagery, as well as produce and serve the 2010–11 LCU systemic dataset of the UMRS floodplain. Digital imagery for Pools 13, 26, La Grange, and Open River South was collected using an Applanix DSS 439 digital sensor system with a 40 millimeter lens and Color Infrared (CIR) filter. The imagery was collected at a resolution of 20 centimeters/pixel (8 inches/pixel) for Pool 13 and 40 centimeters/pixel (16 inches/pixel) for Pools 26, Open River South, and La Grange. All imagery was projected to Universal Transverse Mercator (UTM) Zone 15, North American Datum of 1983 (NAD 83). The General Wetland Vegetation Classification (GWVC) system used for mapping is hierarchical, and its 31 classes can be collapsed into broader categories using either a 15- or 7-class level. While the 1989 and 2000 LCU systemic datasets have not gone through a traditional thematic accuracy assessment (AA) in the past, nor have they undergone a validation analysis, the end products are of high quality. For each systemic dataset produced (1989, 2000, 2010–11), extensive field reconnaissance is performed before photointerpretation. The intent of this field reconnaissance is to learn, test, and verify image signatures as they relate to vegetation types. Questionable areas on the imagery are visited, and the plants or land features observed in the area are recorded for reference. This procedure verifies vegetation signatures on the imagery with those on the ground. In addition, once the photointerpretation is complete, the final LCU dataset undergoes extensive quality assurance/quality control to ensure the imagery is mapped correctly. Since the 2000 LCU systemic dataset was developed, there has been a growing interest in completing thematic AAs for the LTRM LCU spatial datasets. The objective of an AA is to measure the probability that a particular location has been assigned its correct vegetation class. An AA estimates thematic (map class) errors in the data, giving users information needed to determine data suitability for a particular application. At the same time, data producers are able to learn more about the nature of errors in the data. Thus, the two attributes of an AA are “producers’ accuracy,” which is the probability that an AA point has been mapped correctly (also referred to as an error of omission); and “users’ accuracy,” which is the probability that the map actually represents what was found on the ground (also referred to as error of commission). Producers’ and users’ accuracies can be obtained from the same set of data by using different analyses. Accuracy assessment is an extensive effort that requires seasonal field personnel and equipment, data entry, analyses, and post processing—tasks that are costly and time consuming. The geospatial team at the UMESC has suggested a validation process for understanding the accuracy of the spatial datasets, which will be tested on at least some areas of the UMRS. Validation is not a true verification of map-class type in the field; however, it can provide the user of the map with useful information that is similar to a field AA. Similar to an AA, validation involves generating random points based on the total area for each map class. However, instead of collecting field data, two or three individuals not involved with the photo-interpretative mapping separately review each of the points onscreen and record a best-fit vegetation type(s) for each site. Once the individual analyses are complete, results are joined together and a comparative analysis is performed. The objective of this initial analysis is to identify areas where the validation results were in agreement (matches) and areas where validation results were in disagreement (mismatches). The two or three individuals then perform an analysis, looking at each mismatched site, and agree upon a final validation class. (If two vegetation types at a specific site appear to be equally prevalent, the validation team is permitted to assign the site two best-fit vegetation types.) Following the validation team’s comparative analysis of vegetation assignments, the data are entered into a database and compared to the mappers’ vegetation assignments. Agreements and disagreements between the map and validation classes are identified, and a contingency table is produced. This document presents the AA processes/results for Pools 13 and La Grange, as well as the validation process/results for Pools 13 and 26 and Open River South.

Iowa, Illinois, Minnesota, Missouri, Wisconsin↗
USGS2012

The potential effects of sodium bicarbonate, a major constituent from coalbed natural gas production, on aquatic life

The production water from coalbed natural gas (CBNG) extraction contains many constituents. The U.S. Environmental Protection Agency has established aquatic life criteria for some of these constituents, and it is therefore possible to evaluate their risk to aquatic life. However, of the major ions associated with produced waters, chloride is the only one with an established aquatic life criterion (U.S. Environmental Protection Agency, 1988). The focus of this research was NaHCO 3 , a compound that is a major constituent of coalbed natural gas produced waters in the Tongue and Powder River Basins. This project included laboratory experiments, field in situ experiments, a field mixing zone study, and a fishery presence/absence assessment. Though this investigation focuses on the Tongue and Powder River Basins, the information is applicable to other watersheds where sodium bicarbonate is a principle component of product water either from CBNG or from traditional or unconventional oil and gas development. These data can also be used to separate effects of saline discharges from those potentially posed by other constituents. Finally, this research effort and the additional collaboration with USGS Water Resources and Mapping, Bureau of Land Management, US Environmental Protection Agency, State of Montana, State of Wyoming, Montana State University, University of Wyoming, and others as part of a Powder River Aquatic Task Group, can be used as a model for successful approaches to studying landscapes with energy development. The laboratory acute toxicity experiments were completed with a suite of organisms, including 7 species of fish, 5 species of invertebrates, and 1 amphibian species. Experiments performed on these multiple species resulted in LC50s that ranged from 1,120 to greater than (>) 8,000 milligrams sodium bicarbonate per liter (mg NaHCO 3 /L) (also defined as 769 to >8,000 milligrams bicarbonate per liter (mg HCO 3 -/L) or total alkalinity expressed as 608 to >4,181 milligrams calcium carbonate per liter (mg CaCO 3 /L)) that varied across species and lifestage within a species. The age at which fish were exposed to NaHCO 3 significantly affected the severity of toxic responses for some organisms. The chronic toxicity of NaHCO 3 was defined in experiments that lasted from 7—60 days post-hatch. For these experiments, sublethal effects such as growth and reproduction, in addition to significant reductions in survival were included in the final determination of effects. Chronic toxicity was observed at concentrations that ranged from 450 to 800mg NaHCO 3 /L (also defined as 430 to 657 mg HCO 3 -/L or total alkalinity expressed as 354 to 539 mg CaCO 3 /L) and the specific concentration depended on the sensitivity of the four species of invertebrates and fish exposed. Sublethal investigations during chronic studies revealed percent decrease in the activity of sodium-potassium adenosine triphosphatase (Na/K ATPase, an enzyme involved in ionoregulation) and the age of the fish at the onset of the decrease may affect the ability of fathead minnow to survive exposures to NaHCO 3 . A database of toxicity evaluations of NaHCO 3 on aquatic life has been constructed. Using these data, sample acute and chronic criteria of 459 and 381 mg NaHCO 3 /L, respectively, can be calculated for the protection of aquatic life. The final derivation and implementation of such criteria is, of course, left to the discretion of the concerned management agencies. A combination of in situ experiments, static-renewal experiments performed simultaneously with in situ experiments, and static renewal experiments performed with site water in the laboratory, demonstrated that untreated coalbed natural gas (CBNG) product water from the Tongue and Powder River Basins reduces survival of fathead minnow and pallid sturgeon. More precisely, the survival of early-lifestage fathead minnow, especially those less than 6-days post hatch (dph), likely is reduced significantly in the field when concentrations of NaHCO 3 rise above 1,500 mg/L. However, age was not a factor for pallid sturgeon and they were sensitive to product water regardless of age. Treatment with the Higgins Loop™ technology and dilution of untreated water increased survival in the laboratory. Both of these situations reduced ammonia in addition to the concentrations of NaHCO 3 . These experiments addressed the acute toxicity of effluent waters being added to the main stem rivers, but did not address issues related to the volumes of water that may be added to the watershed. Mixing zones of the three outfalls studied ranged from approximately 800—1,200 m below the confluence and the areas within these mixing zones with acutely lethal concentrations of NaHCO 3 (as defined by the presence of concentrated dye) are limited. The areas with concentrations of NaHCO 3 more than the concentrations likely to cause significant mortality, and more than the calculated sample water-quality criteria in the Tongue and Powder River Basins appear to be limited to tributaries and parts of mixing zones with considerable additions of untreated effluent.

Montana;Wyoming↗
USGS2012

Hydrological information products for the Off-Project Water Program of the Klamath Basin Restoration Agreement

The Klamath Basin Restoration Agreement (KBRA) was developed by a diverse group of stakeholders, Federal and State resource management agencies, Tribal representatives, and interest groups to provide a comprehensive solution to ecological and water-supply issues in the Klamath Basin. The Off-Project Water Program (OPWP), one component of the KBRA, has as one of its purposes to permanently provide an additional 30,000 acre-feet of water per year on an average annual basis to Upper Klamath Lake through "voluntary retirement of water rights or water uses or other means as agreed to by the Klamath Tribes, to improve fisheries habitat and also provide for stability of irrigation water deliveries." The geographic area where the water rights could be retired encompasses approximately 1,900 square miles. The OPWP area is defined as including the Sprague River drainage, the Sycan River drainage downstream of Sycan Marsh, the Wood River drainage, and the Williamson River drainage from Kirk Reef at the southern end of Klamath Marsh downstream to the confluence with the Sprague River. Extensive, broad, flat, poorly drained uplands, valleys, and wetlands characterize much of the study area. Irrigation is almost entirely used for pasture. To assist parties involved with decisionmaking and implementation of the OPWP, the U.S. Geological Survey (USGS), in cooperation with the Klamath Tribes and other stakeholders, created five hydrological information products. These products include GIS digital maps and datasets containing spatial information on evapotranspiration, subirrigation indicators, water rights, subbasin streamflow statistics, and return-flow indicators. The evapotranspiration (ET) datasets were created under contract for this study by Evapotranspiration, Plus, LLC, of Twin Falls, Idaho. A high-resolution remote sensing technique known as Mapping Evapotranspiration at High Resolution and Internalized Calibration (METRIC) was used to create estimates of the spatial distribution of ET. The METRIC technique uses thermal infrared Landsat imagery to quantify actual evapotranspiration at a 30-meter resolution that can be related to individual irrigated fields. Because evaporation uses heat energy, ground surfaces with large ET rates are left cooler as a result of ET than ground surfaces that have less ET. As a consequence, irrigated fields appear in the Landsat images as cooler than nonirrigated fields. Products produced from this study include total seasonal and total monthly (April-October) actual evapotranspiration maps for 2004 (a dry year) and 2006 (a wet year). Maps showing indicators of natural subirrigation were also provided by this study. "Subirrigation" as used here is the evapotranspiration of shallow groundwater by plants with roots that penetrate to or near the water table. Subirrigation often occurs at locations where the water table is at or above the plant rooting depth. Natural consumptive use by plants diminishes the benefit of retiring water rights in subirrigated areas. Some agricultural production may be possible, however, on subirrigated lands for which water rights are retired. Because of the difficulty in precisely mapping and quantifying subirrigation, this study presents several sources of spatially mapped data that can be used as indicators of higher subirrigation probability. These include the floodplain boundaries defined by stream geomorphology, water-table depth defined in Natural Resources Conservation Service (NRCS) soil surveys, and soil rooting depth defined in NRCS soil surveys. The two water-rights mapping products created in the study were "points of diversion" (POD) and "place of use" (POU) for surface-water irrigation rights. To create these maps, all surface-water rights data, decrees, certificates, permits, and unadjudicated claims within the entire 1,900 square mile study area were aggregated into a common GIS geodatabase. Surface-water irrigation rights within a 5-mile buffer of the study area were then selected and identified. The POU area was then totaled by water right for primary and supplemental water rights. The maximum annual volume (acre-feet) allowed under each water right also was calculated using the POU area and duty (allowable annual irrigation application in feet). In cases where a water right has more than one designated POD, the total volume for the water right was equally distributed to each POD listed for the water right. Because of this, mapped distribution of diversion rates for some rights may differ from actual practice. Water-right information in the map products was from digital datasets obtained from the Oregon Water Resources Department and was, at the time acquired, the best available compilation of water-right information available. Because the completeness and accuracy of the water-right data could not be verified, users are encouraged to check directly with the Oregon Water Resources Department where specific information on individual rights or locations is essential. A dataset containing streamflow statistics for 72 subbasins in the study area was created for the study area. The statistics include annual flow durations (5-, 10-, 25-, 50-, and 95-percent exceedances) and 7-day, 10-year (7Q10) and 7-day, 2-year (7Q2) low flows, and were computed using regional regression equations based on measured streamflow records in the region. Daily streamflow records used were adjusted as needed for crop consumptive use; therefore the statistics represent streamflow under more natural conditions as though irrigation diversions did not exist. Statistics are provided for flow rates resulting from streamflow originating from within the entire drainage area upstream of the subbasin pour point (referring to the outlet of the contributing drainage basin). The statistics were computed for the purpose of providing decision makers with the ability to estimate streamflow that would be expected after water conservation techniques have been implemented or a water right has been retired. A final product from the study are datasets of indicators of the potential for subsurface return flow of irrigation water from agricultural areas to nearby streams. The datasets contain information on factors such as proximity to surface-water features, geomorphic floodplain characteristics, and depth to water. The digital data, metadata, and example illustrations for the datasets described in this report are available on-line from the USGS Water Resources National Spatial Data Infrastructure (NSDI) Node Website http://water.usgs.gov/lookup/getgislist or from the U.S. Government website DATA.gov at http://www.data.gov with links provided in a Microsoft® Excel® workbook in appendix A.

California;Oregon↗
USGS1993

Flow recommendations for maintaining riparian vegetation along the Upper Missouri River, Montana

Montana Power Company, Inc. (MPC) submitted a final license application to the Federal Energy Regulatory Commission (FERC) on November 30, 1992. In this application, MPC proposed a plan for the protection of fish, wildlife, habitat, and water-quality resources. One concern was maintenance of woody riparian vegetation along the Missouri River, especially along the Wild and Scenic reach of the river, where the riparian forest occurs in relatively small discontinuous stands. The objectives of this project were 1) to recommend flows that would protect and enhance riparian forests along the Missouri River, and 2) to develop elements of an environmental monitoring program that could be used to assess the effectiveness of the recommended flows. Plains cottonwood ( Populus deltoides subsp. monilifera ) is the key structural component of riparian forests along the Missouri River. Therefore, we focused our analysis on factors affecting populations of this species. Previous work had demonstrated that the age structure of cottonwood populations is strongly influenced by aspects of flow that promote successfully establishment. In this study our approach was to determine the precise age of plains cottonwood trees growing along the Upper Missouri River and to relate years of establishment to the flow record. Our work was carried out between Coal Banks Landing and the Fred G. Robinson Bridge within the Wild and Scenic portion of the Missouri River. This segment of the river occupies a narrow valley and exhibits little channel migration. Maps and notes from the journals of Lewis and Clark (1804-1806) suggest that the present distribution and abundance of cottonwoods within the study reach is generally similar to presettlement conditions. Flows in the study reach are influenced by a number of dams and diversions, most importantly, Canyon Ferry and Tiber Dams. Although flow regulation has decreased peak flows and increased low flows, the gross seasonal pattern of flow has not been greatly altered. Most cottonwood establishment in our study reach occurred in years with a peak mean daily flow greater than 1,400 m 3 /s (49,434 cfs), or in the two years following such a flow. These years include 35 out of the 111 years of record, and account for establishment of 47 of 60 trees examined, a highly significant relationship. Infrequent establishment of cottonwood trees is not the result of scarcity of seed or seedlings. In the study reach seedlings become established most years on bare, relatively low surfaces deposited by the river. However, the high elevation of establishment of all trees dating to before 1978 indicates that only individuals established on high flood deposits are able to survive subsequent floods and ice jams. In order to maintain the present abundance of plains cottonwood in the study area we recommend flood flows in excess of 1,400 m 3 /s (49,434 cfs) measures as mean daily discharge at Fort Benton (U.S. Geological Survey gage 06090800) with a recurrence interval of approximately 9 years. Because cottonwood seeds remain viable for only a few weeks, and because seedling require a moist, bare surface, we further recommend maintenance of the historic timing of flooding with peak flood flows occurring between mid-May and late-June. Flow is not the only factor influencing cottonwood regeneration along this reach of the Missouri River. Land management, especially cattle grazing, is clearly having an impact, and changes in cottonwood populations could be expected if these practices were altered. However, the dependence of cottonwood establishment on high flow is clear in this reach in spite of the effects of other factors. Given the value of the resource, we strongly suggest establishment of a monitoring program to determine the effectiveness of the recommended flows and to provide the data necessary for refining them. We recommend a monitoring program that would include: 1) ten permanent, widely space channel cross sections for annual measurement of channel geometry and cottonwood establishment, growth, and survival; 2) five livestock enclosures to monitor the influence of grazing in the study area; and 3) low-elevation aerial photography of the reach every five years and after every flood to detect changes in channel geometry and forested area. Because cottonwood establishment is episodic, a long-term commitment to the monitoring effort is essential. In addition, cross sections and exclosures should be easy enough to access that measurements during flood years are possible.

Montana↗
USGS2012

Simulating potential structural and operational changes for Detroit Dam on the North Santiam River, Oregon, for downstream temperature management

Detroit Dam was constructed in 1953 on the North Santiam River in western Oregon and resulted in the formation of Detroit Lake. With a full-pool storage volume of 455,100 acre-feet and a dam height of 463 feet, Detroit Lake is one of the largest and most important reservoirs in the Willamette River basin in terms of power generation, recreation, and water storage and releases. The U.S. Army Corps of Engineers operates Detroit Dam as part of a system of 13 reservoirs in the Willamette Project to meet multiple goals, which include flood-damage protection, power generation, downstream navigation, recreation, and irrigation. A distinct cycle in water temperature occurs in Detroit Lake as spring and summer heating through solar radiation creates a warm layer of water near the surface and isolates cold water below. Controlling the temperature of releases from Detroit Dam, therefore, is highly dependent on the location, characteristics, and usage of the dam's outlet structures. Prior to operational changes in 2007, Detroit Dam had a well-documented effect on downstream water temperature that was problematic for endangered salmonid fish species, releasing water that was too cold in midsummer and too warm in autumn. This unnatural seasonal temperature pattern caused problems in the timing of fish migration, spawning, and emergence. In this study, an existing calibrated 2-dimensional hydrodynamic water-quality model [CE-QUAL-W2] of Detroit Lake was used to determine how changes in dam operation or changes to the structural release points of Detroit Dam might affect downstream water temperatures under a range of historical hydrologic and meteorological conditions. The results from a subset of the Detroit Lake model scenarios then were used as forcing conditions for downstream CE-QUAL-W2 models of Big Cliff Reservoir (the small reregulating reservoir just downstream of Detroit Dam) and the North Santiam and Santiam Rivers. Many combinations of environmental, operational, and structural options were explored with the model scenarios. Multiple downstream temperature targets were used along with three sets of environmental forcing conditions representing cool/wet, normal, and hot/dry conditions. Five structural options at Detroit Dam were modeled, including the use of existing outlets, one hypothetical variable-elevation outlet such as a sliding gate, a hypothetical combination of a floating outlet and a fixed-elevation outlet, and a hypothetical combination of a floating outlet and a sliding gate. Finally, 14 sets of operational guidelines for Detroit Dam were explored to gain an understanding of the effects of imposing different downstream minimum streamflows, imposing minimum outflow rules to specific outlets, and managing the level of the lake with different timelines through the year. Selected subsets of these combinations of operational and structural scenarios were run through the downstream models of Big Cliff Reservoir and the North Santiam and Santiam Rivers to explore how hypothetical changes at Detroit Dam might provide improved temperatures for endangered salmonids downstream of the Detroit-Big Cliff Dam complex. Conclusions that can be drawn from these model scenarios include: *The water-temperature targets set by the U.S. Army Corps of Engineers for releases from Detroit Dam can be met through a combination of new dam outlets or a delayed drawdown of the lake in autumn. *Spring and summer dam operations greatly affect the available release temperatures and operational flexibility later in the autumn. Releasing warm water during midsummer tends to keep more cool water available for release in autumn. *The ability to meet downstream temperature targets during spring depends on the characteristics of the available outlets. Under existing conditions, although warm water sometimes is present at the lake surface in spring and early summer, such water may not be available for release if the lake level is either well below or well above the spillway crest. *Managing lake releases to meet downstream temperature targets depends on having outlet structures that can access both (warm) lake surface water and (cold) deeper lake water throughout the year. The existing outlets at Detroit Dam do not allow near-surface waters to be released during times when the lake surface level is below the spillway (spring and autumn). *Using the existing outlets at Detroit Dam, lake level management is important to the water temperature of releases because it controls the availability and depth of water at the spillway. When lake level is lowered below the spillway crest in late summer, the loss of access to warm water at the lake surface can result in abrupt changes to release temperatures. *Because the power-generation intakes (penstocks) are 166 feet below the full-pool lake level, imposing minimum power production requirements at Detroit Dam limits the amount of warm surface water that can be expelled from the lake in midsummer, thereby postponing and amplifying warm outflows from Detroit Lake into the autumn spawning season. *Likewise, imposing minimum power production requirements at Detroit Dam in autumn can limit the amount of cool hypolimnetic water that is released from the lake, thereby limiting cool outflows from Detroit Lake during the autumn spawning season. *Model simulations indicate that a delayed drawdown of Detroit Lake in autumn would result in better control over release temperatures in the immediate downstream vicinity of Big Cliff Dam, but the reduced outflows necessary to retain more water in the lake in late summer are more susceptible to rapid heating downstream. *Compared to the existing outlets at Detroit Dam, floating or sliding-gate outlet structures can provide greater control over release temperatures because they provide better access to warm water at the lake surface and cooler water at depth. These conclusions can be grouped into several common themes. First, optimal and flexible management and achievement of downstream temperature goals requires that releases of warm water near the surface of the lake and cold water below the thermocline are both possible with the available dam outlets during spring, summer, and autumn. This constraint can be met to some extent with existing outlets, but only if access to the spillway is extended into autumn by keeping the lake level higher than called for by the current rule curve (the typical target water-surface elevation throughout the year). If new outlets are considered, a variable-elevation outlet such as a sliding gate structure, or a floating outlet in combination with a fixed-elevation outlet at sufficient depth to access cold water, is likely to work well in terms of accessing a range of water temperatures and achieving downstream temperature targets. Furthermore, model results indicate that it is important to release warm water from near the lake surface during midsummer. If not released downstream, the warm water will build up at the top of the lake as a result of solar energy inputs and the thermocline will deepen, potentially causing warm water to reach the depth of deeper fixed-elevation outlets in autumn, particularly when the lake level is drawn down to make room for flood storage. Delaying the drawdown in autumn can help to keep the thermocline above such outlets and preserve access to cold water. Although it is important to generate hydropower at Detroit Dam, minimum power-production requirements limit the ability of dam operators to meet downstream temperature targets with existing outlet structures. The location of the power penstocks below the thermocline in spring and most of summer causes the release of more cool water during summer than is optimal. Reducing the power-production constraint allows the temperature target to be met more frequently, but at the cost of less power generation. Finally, running the Detroit Dam, Big Cliff Dam, and North Santiam and Santiam River models in series allows dam operators to evaluate how different operational strategies or combinations of new dam outlets might affect downstream temperatures for many miles of critical endangered salmonid habitat. Temperatures can change quickly in these downstream reaches as the river exchanges heat with its surroundings, and heating or cooling of 6 degrees Celsius is not unusual in the 40–50 miles downstream of Big Cliff Dam. The results published in this report supersede preliminary results published in U.S. Geological Survey Open-File Report 2011-1268 (Buccola and Rounds, 2011). Those preliminary results are still valid, but the results in this report are more current and comprehensive.

Oregon↗
USGS1993

Field observations, preliminary model analysis, and aquifer thermal efficiency

In May 1980, the University of Minnesota began a project to evaluate the feasibility of storing heated (150 degrees Celsius (°C) water in the deep (180 to 240 meters (m)) Franconia-Ironton-Galesville aquifer and later recovering it for space heating. The Aquifer Thermal-Energy Storage (ATES) system was a doublet-well design in which the injection and withdrawal wells were spaced approximately 250 m apart. High-temperature water from the university's steam-generation facilities supplied heat for injection. Water was pumped from one of the wells through a heat exchanger, where heat was added or removed. Water then was injected back into the aquifer through the other well. The experimental plan for testing the ATES system consisted of a series of short-term hot-water injection, storage, and withdrawal cycles. Each cycle was 24 days long, and each injection, storage, and withdrawal step of the cycle was 8 days. The Franconia-Ironton-Galesville aquifer is a consolidated sandstone, approximately 60 m thick, the top of which is approximately 180 m below the land surface. It is confined above by the St. Lawrence Formation a dolomitic sandstone 8 m thick and below by the Eau Claire Formation a shale 30 m thick. Initial hydraulic testing with inflatable packers indicated that the aquifer has four hydraulic zones with distinctly different values of relative horizontal hydraulic conductivity. The thickness of each zone was determined by correlating data from geophysical logs, core samples, and the inflatable-packer tests. A comprehensive network for data collection, storage, and analysis was designed to monitor temperature and pressure changes during the ATES test cycles. A total of 22 pressure transducers and 56 thermocouples monitored pressures and temperatures in the aquifer and in the upper and lower confining units. Temperature and pressure measurements were collected in observation well nests at distances of approximately 7 and 14 m from the production wells. All pressure and temperature data were transmitted through buried cables to a central data logger, where the measurements were viewed independently or stored on computer magnetic tape for later analysis. Interactive computer programs were available to display data stored on magnetic tapes as individual measurements or as plots of pressure and temperature versus time. Analyses of step-drawdown and constant-discharge aquifer tests indicate that the Franconia-Ironton-Galesville aquifer is anisotropic in the horizontal plane. Major and minor transmissivities are 101.5 and 44.6 m 2 /d (square meters per day), respectively, for the Ironton and Galesville Sandstones and 40.0 and 24.0 m 2 /d, respectively, for the upper part of the Franconia Formation. The average transmissivity of the entire Franconia-Ironton-Galesville aquifer is about 98 m 2 /d. Effective porosity ranges from 0.25 to 0.31, and the average storage coefficient is 4.5x10 -5 . Two computer models were constructed to simulate the movement of ground water and heat. The first was a nonisothermal, isotropic, single-phase, radial, ground-water flow and thermal-energy-transport model that was constructed to examine the sensitivity of model results to various hydraulic and thermal properties. The model also was used to study the potential for buoyancy flow within the aquifer and the effect of various cyclic injection and withdrawal schemes on the relative thermal efficiency of the aquifer. The second model was a threedimensional ground-water flow and thermal-energy-transport model that was constructed to incorporate the anisotropy of the aquifer. In the first model, the sensitivity analysis assumed 8 days of injection of 150°C water at 18.9 liters per second (L/s), 8 days of storage, and 8 days of withdrawal of hot water at 18.9 L/s. The analysis indicates that, for practical ranges of hydraulic and thermal properties, the ratio of horizontal to vertical hydraulic conductivity is the least important property and thermal dispersivity is the most important property used to compute temperature and aquifer thermal efficiency. Buoyancy flow was examined for several values of hydraulic conductivity and ratios of horizontal to vertical hydraulic conductivities. For the assumed base values of hydraulic and thermal properties, buoyancy flow was negligible. The greatest simulated buoyancy flow resulted from simulations in which horizontal hydraulic conductivity was increased to ten times the base value and in which the vertical hydraulic conductivity was set equal to the horizontal hydraulic conductivity. The effects of various injection and withdrawal rates and durations on computed values of aquifer relative-thermal efficiency and final well-bore temperature were studied for five 1-year hypothetical test cycles of injection and withdrawal. The least efficient scheme was 8 months injection of 150°C water at 18.9 L/s and 4 months of withdrawal of hot water at 18.9 L/s. The most efficient scheme was obtained with 6 months of injection of 150°C water at 18.9 L/s and 6 months of withdrawal of hot water at 37.8 L/s. The hypothetical simulations indicate that the calibrated model of the doublet-well system would be a valuable tool for use by the university in selecting a highly efficient system operation. In the second model, analytical solutions of anisotropic hydraulic flow around the doublet-well system were obtained to provide fluid-flux boundary conditions around the heat-injection well in three dimensions. This information simplified simulation of the doublet-well system because only the heat injection well needed to be simulated. This second model was calibrated with data from an 8-day ambienttemperature injection test at 18.9 L/s. Boundary-flux conditions were examined for nonisothermal conditions by simulating 8 days of injection of 150°C water at 18.9 L/s. Results of simulations using both models indicate that the fluxboundary conditions are adequate for simulations of short-term heatinjection testing.

Professional Paper↗
USGS2016

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↗
USGS2012

Simulation of groundwater flow and hydrologic effects of groundwater withdrawals from the Kirkwood-Cohansey aquifer system in the Pinelands of southern New Jersey

The Kirkwood-Cohansey aquifer system is an important source of present and future water supply in southern New Jersey. Because this unconfined aquifer system also supports sensitive wetland and aquatic habitats within the New Jersey Pinelands (Pinelands), water managers and policy makers need up-to-date information, data, and projections that show the effects of potential increases in groundwater withdrawals on these habitats. Finite-difference groundwater flow models (MODFLOW) were constructed for three drainage basins (McDonalds Branch Basin, 14.3 square kilometers (km 2 ); Morses Mill Stream Basin, 21.63 km 2 ; and Albertson Brook Basin, 52.27 km 2 ) to estimate the effects of potential increases in groundwater withdrawals on water levels and the base-flow portion of streamflow, in wetland and aquatic habitats. Three models were constructed for each drainage basin: a transient model consisting of twenty-four 1-month stress periods (October 2004 through September 2006); a transient model to simulate the 5- to 10-day aquifer tests that were performed as part of the study; and a high-resolution, steady-state model used to assess long-term effects of increased groundwater withdrawals on water levels in wetlands and on base flow. All models were constructed with the same eight-layer structure. The smallest horizontal cell dimensions among the three model areas were 150 meters (m) for the 24-month transient models, 10 m for the steady-state models, and 3 m for the transient aquifer-test models. Boundary flows of particular interest to this study and represented separately are those for wetlands, streams, and evapotranspiration. The final variables calibrated from both transient models were then used in steady-state models to assess the long-term effects of increased groundwater withdrawals on water levels in wetlands and on base flow. Results of aquifer tests conducted in the three study areas illustrate the effects of withdrawals on water levels in wetlands and on base flow. Pumping stresses at aquifer-test sites resulted in measurable drawdown in each observation well installed for the tests. The magnitude of drawdown in shallow wetland observation wells at the end of pumping ranged from 5.5 to 16.7 centimeters (cm). The stresses induced by the respective tests reduced the flow of the smallest stream (McDonalds Branch) by 75 percent and slightly reduced flow in a side channel of Morses Mill Stream, but did not measurably affect the flow of Morses Mill Stream or Albertson Brook. Results of aquifer-test simulations were used to refine the estimates of hydraulic properties used in the models and to confirm the ability of the model to replicate observed hydrologic responses to pumping. Steady-state sensitivity simulation results for a variety of single well locations and depths were used to define overall “best-case” (smallest effect on wetland water levels and base flow) and “worst-case” (greatest effect on wetland water levels and base flow) groundwater withdrawal configurations. “Best-case” configurations are those for which the extent of the wetland areas within a 1-kilometer (km) radius of the withdrawal well is minimized, the well is located at least 100 m and as far from wetland boundaries as possible, and the withdrawal is from a deep well (50–90 m deep). “Worst-case” configurations are those for which the extent of wetlands within a 1-km radius of the withdrawal well is maximized, the well is located 100 m or less from a wetland boundary, and the withdrawal is from a relatively shallow well (30–67 m deep). “Best-” and “worst-case” simulations were applied by locating hypothetical wells across the study areas and assigning groundwater withdrawals so that the sum of the withdrawals for the basin is equal to 5, 10, 15, and 30 percent of overall recharge. The results were compared to the results of simulations of no groundwater withdrawals. Results for withdrawals of 5 percent of recharge show that the area of wetland water-level decline that exceeded 15 cm was as much as 1.5 percent of the total wetland area for the “best-case” simulations and as much as 9.7 percent of the total wetland area for the “worst-case” simulations. For the same withdrawals, base-flow reduction was as much as 5.1 percent for the “best-case” simulations and as much as 8.6 percent for the “worst-case” simulations. Results for withdrawals of 30 percent of recharge show that the area of wetland water-level decline that exceeded 15 cm was as much as 70 percent of the total wetland area for the “best-case” simulations and as much as 84 percent of the total wetland area for the “worst-case” simulations. For the same withdrawals, base-flow reduction was as much as 30 percent for the “best-case” simulations and as much as 51 percent for the “worst-case” simulations. Results for withdrawals of 10 and 15 percent of recharge show decreased water levels and base flow that are intermediate between those simulated for 5 and 30 percent of recharge. Several approaches for applying the results of this study to other parts of the Pinelands were explored. An analytical-modeling technique based on the Thiem equation and image-well theory was developed to estimate local drawdown distributions resulting from withdrawals in other areas within the Pinelands. Results of example applications of this technique were compared with those of the numerical simulations used in this study and were shown to be useful. Differences among the three basins in the simulated percentage of basin wetlands affected by drawdown were found to be related to the proximity of wetlands to streams, the proximity of wetlands to pumped wells, and the vertical conductance of the aquifer system. These factors formed the basis for an index of wetland vulnerability to drawdown. An empirically-derived model based on the Gompertz function and the wetland vulnerability index was developed, tested, and shown to be an effective means to evaluate potential drawdown in wetlands at a basin scale throughout the Pinelands. Base-flow reduction can be estimated from generalized results of the numerical models, estimates of evapotranspiration reduction, or available regional groundwater flow models. These approaches could be used to evaluate alternative water-supply strategies and, in conjunction with ecological-modeling results, to determine maximum basin withdrawal rates within the limits of acceptable ecological change.

New Jersey↗