Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Nature, Scientific Data”

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

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

At least 1,153 records · Page 64Linked to original sources

Methods for estimating magnitude and frequency of floods in Arizona, developed with unregulated and rural peak-flow data through water year 2010

Flooding is among the worst natural disasters responsible for loss of life and property in Arizona, underscoring the importance of accurate estimation of flood magnitude for proper structural design and floodplain mapping. Twenty-four years of additional peak-flow data have been recorded since the last comprehensive regional flood frequency analysis conducted in Arizona. Periodically, flood frequency estimates and regional regression equations must be revised to maintain the accurate estimation of flood frequency and magnitude. Annual peak-flow data collected through water year 2010 were compiled from 448 unregulated streamflow-gaging stations, hereafter referred to as streamgages, in Arizona having a minimum of 10 years of record. Flood frequency estimates were first computed with station (or at-site) skew using the Expected Moments Algorithm with a multiple Grubbs-Beck test to identify multiple potentially influential low flows to fit a Pearson Type III distribution. Next, a multiple step Bayesian least-squares-regression approach was used to determine a new statewide regional skew of −0.09. No basin characteristics analyzed were statistically significant in explaining the variation in skew and as a result, the constant model was chosen as the best regional skew model for the Arizona study area. The mean square error used in Bulletin 17B (B17B) of the Interagency Advisory Committee on Water Data is used to describe the precision of the regional skew. The constant model had a mean square error equal to 0.08, which corresponds to an effective record length of 85 years. This is a marked improvement over a previous Arizona regional skew analysis, with a reported mean square error of 0.31, for a corresponding effective record length of around 17 years. Thus the new regional model had almost five times the information content (as measured by effective record length) of that calculated in USGS Water Supply Paper 2433, published in 1997, or the value of 0.302 reported in the B17B generalized skew map. The flood frequency estimates were recalculated using a weighted skew of the station and regional skew. Station flood frequency estimates for each streamgage are presented for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities. Geographical information systems were used to compute basin characteristic information for each streamgage for the purpose of developing regional equations to estimate flood statistics at ungaged basins. Five hydrologic flood regions in Arizona were defined in a multivariate regionalization process based on mean basin elevation, mean annual precipitation, and soil permeability. A regional generalized least-squares-regression analysis was used to develop five sets of equations from 344 nonredundant streamgages, corresponding to five regions, for estimating the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities at ungaged basins in Arizona. The regression equations developed for these five regions were based on one or more of the statistically significant explanatory variables: drainage area, mean basin elevation, and mean annual precipitation. Average standard errors of prediction for the regression regions for the five regions ranged from 27 to 122 percent and the pseudo-coefficients of determination (pseudo-R 2 ), a measure of the proportion of peak-flow variation that is explained by the basin characteristics, ranged from 68 to 98 percent. Regression equations for Central Highlands (region 4) had the lowest model error and the greatest pseudo-R 2 metrics. The equations for Colorado Plateau (region 2) regression equations generally had greater model error and lower pseudo-R 2 metrics. The improvement of regional regression equation model error and pseudo-R 2 metrics was related to higher numbers of streamgages, longer period of record, and even spatial coverage within a region. The regional regression equations were integrated into the U.S. Geological Survey’s StreamStats program. The StreamStats program is a national map-based web application that allows the public to easily access published flood frequency and basin characteristic statistics. The interactive web application allows a user to select a point within a watershed (gaged or ungaged) and retrieve flood-frequency estimates derived from the current regional regression equations and geographic information system data within the selected basin. StreamStats provides users with an efficient and accurate means for retrieving the most up to date flood frequency and basin characteristic data. StreamStats is intended to provide consistent statistics, minimize user error, and reduce the need for large datasets and costly geographic information system software.

Arizona↗

General water-quality conditions, long-term trends, and network analysis at selected sites within the Ambient Water-Quality Monitoring Network in Missouri, water years 1993–2017

The U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources, collects data pertaining to the surface-water resources of Missouri. Established in 1964, the Ambient Water-Quality Monitoring Network (AWQMN) consisted of 69 sites in 2017. Two additional sites from the National Water-Quality Program are included with the AWQMN sites for the analyses in this report. The sites are sampled typically from 2 to 12 times per year for physical properties, total suspended solids, nutrients, fecal indicator bacteria, and trace elements. The period of analysis for this study was from 1993 through 2017 and data analysis included 71 sites and 15 water-quality constituents plus discharge. Data analysis involved retrieving the data, conditioning the data for analysis, analyzing the data for trends, and analyzing the monitoring network to determine if potential data gaps or data redundancies exist in the network. Results from these analyses can be used to help manage the monitoring network into the future. Water-quality data were analyzed using several software packages to provide graphical and statistical information for interpretation of trends in the data at selected sites. Discharge data at selected sites were analyzed to determine the general trends during the analysis period and how the water-quality samples represented the range of daily mean discharges at each site. Water-quality data also were analyzed at selected sites to determine the relative sensitivity of selected sites and constituents to changes in data collection frequency. Trend analysis at selected sites using a simulated reduction in sampling frequency was completed to compare to trends obtained using monthly data to determine the potential degradation in the ability of determining trends from a reduced sampling frequency. The viability of using estimated discharge to evaluate long-term trends for sites with no continuous discharge was investigated. Data from sites were statistically compared in groups to determine the relative similarity (or difference) between sites for each water-quality constituent to identify potentially redundant sites in the monitoring network. Discharge-weighted long-term trends during 1993 through 2017 were analyzed for 15 water-quality constituents at 58 sites and results indicated there were significant single- or two-period trends in about 17 percent of the analyses. Some trends indicated improvement and some trends indicated deterioration of the general water quality at some sites in the AWQMN. No trend was indicated in about 31 percent of the analyses. The constituents pH, specific conductance, and total phosphorus showed the most frequent significant trends, and each of the 15 constituents examined had a significant trend at one or more sites. A total of 42 sites indicated at least 1 constituent with a significant single- or two-period trend, and 10 sites indicated 6 or more significant trends. Potential data gaps identified for computing discharge-weighted long-term trends in the monitoring network included the lack of collection of continuous discharge at 23 sites, insufficient sampling frequency for some constituents (dissolved chloride and total and dissolved lead and zinc) at some sites, insufficient temporal sample distribution (lack of at least one sample in each season per year) at some sites, and insufficient sampling frequency for some highly censored constituents (nutrients and total and dissolved lead and zinc) at some sites. Potential data gaps based on site spatial distribution were identified in 7 basins greater than 800 square miles. Potential site redundancies were identified in 4 basins that had an area greater than 500 square miles with a site density greater than 2 sites per 1,000 square miles. Potential site redundancies also were identified for nine site pairs by observing statistical similarities in the constituent data distributions. Sampling frequency was investigated to determine if reducing the sampling frequency of select constituents could provide a statistically similar data distribution. At 28 of 71 sites, 11 constituents had sufficient data collection frequency (approximately monthly) to allow for the creation of simulated datasets of various reduced data collection frequency. For the selected monitoring network sites analyzed, the data distribution of a simulated sampling frequency of four times per year or greater, roughly evenly distributed over the year, was not significantly different than the data distribution of the original monthly sampling frequency. Sites analyzed using varying simulated sampling frequencies tended to be more sensitive to sampling frequency changes if they were in basins classified as large or very large size and tended to be least sensitive in basins classified as small and medium size in the Ozark Plateaus Province. Simulated reduced frequency sampling analysis indicated that the constituents and measurements most sensitive to changes in sampling frequencies were water temperature, dissolved oxygen, discharge, and dissolved nitrate, and least sensitive were pH, total suspended solids, dissolved phosphorus, and total phosphorus. Discharge-weighted long-term trend analysis was repeated at 22 sites for 11 constituents using a simulated quarterly sampling frequency, and matched about 46 percent of the significant single-period trends identified using monthly data and about 65 percent of the analyses that indicated no trend using the monthly data.

Missouri↗

Benthic habitats and offshore geological resources of Kaloko-Honokohau National Historical Park, Hawai‘i

Kaloko-Honokōhau National Historical Park (KAHO) is one of three National Park lands located along the western coast of the Island of Hawai‘i and the only one to include submerged lands and marine resources within its official boundaries. The park was established in 1978 and is 1,160 acres in size, including 596 acres of marine area. The submerged lands are currently managed by the State of Hawaii, Department of Land and Natural Resources, Division of Aquatic Resources (DLNR-DAR). Marine resources located within KAHO include coral reef and habitat for many marine animals, including the green sea turtle and a variety of fish and invertebrates. In addition, many archeological, cultural, and recreational resources are located within the marine realm of the park. Potential threats and stressors to the modern marine environment include ground-water and surface-water contamination, invasive plants and algae, fishing pressure and use of monofilament gill nets (which can ensnare marine life or become tangled on reefs and be left behind as fishing debris), and visitor use impacts, such as scuba diving and snorkeling. Illegal dumping, oil releases, boat groundings, and other physical damage to reef resources are potential threats from users of the nearby harbor. Special issues of concern for the park include establishing baseline conditions of the offshore resources before the development of adjacent coastal lands. Until this study, only a general knowledge of the distribution of benthic habitats and the characteristics of the offshore region of Kaloko-Honokōhau National Historical Park was available. In 2003, a collaborative project between the U.S. Geological Survey (USGS) Coastal and Marine Geology Program and the National Park Service (NPS) was initiated to develop detailed benthic-habitat classification maps for the marine lands within and adjacent to the park. The intent of this project is to provide baseline maps and a Geographic Information System (GIS) database and description of the biological and geological resources of these marine lands in order to facilitate the management, interpretation, and understanding of park resources. A benthic-habitat classification map was created for the park using existing color aerial photography, Scanning Hydrographic Operational Airborne Lidar Survey (SHOALS) bathymetric data, georeferenced underwater video, and still photography. Individual habitat polygons were classified using five basic attributes: (1) major structure or substrate, (2) dominant structure, (3) major biologic cover on the substrate, (4) percentage of major biological cover, and (5) geographic zone. Additional information regarding geology, morphology, and coral species were also noted. Marine resources located within KAHO include coral reef and habitat for many marine animals, including the green sea turtle and a variety of fish and invertebrates. In addition, many archeological, cultural, and recreational resources are located within the marine realm of the park. Potential threats and stressors to the modern marine environment include ground-water and surface-water contamination, invasive plants and algae, fishing pressure and use of monofilament gill nets (which can ensnare marine life or become tangled on reefs and be left behind as fishing debris), and visitor use impacts, such as scuba diving and snorkeling. Illegal dumping, oil releases, boat groundings, and other physical damage to reef resources are potential threats from users of the nearby harbor. Special issues of concern for the park include establishing baseline conditions of the offshore resources before the development of adjacent coastal lands. Until this study, only a general knowledge of the distribution of benthic habitats and the characteristics of the offshore region of Kaloko-Honokōhau National Historical Park was available. In 2003, a collaborative project between the U.S. Geological Survey (USGS) Coastal and Marine Geology Program and the National Park Service (NPS) was initiated to develop detailed benthic-habitat classification maps for the marine lands within and adjacent to the park. The intent of this project is to provide baseline maps and a Geographic Information System (GIS) database and description of the biological and geological resources of these marine lands in order to facilitate the management, interpretation, and understanding of park resources. A benthic-habitat classification map was created for the park using existing color aerial photography, Scanning Hydrographic Operational Airborne Lidar Survey (SHOALS) bathymetric data, georeferenced underwater video, and still photography. Individual habitat polygons were classified using five basic attributes: (1) major structure or substrate, (2) dominant structure, (3) major biologic cover on the substrate, (4) percentage of major biological cover, and (5) geographic zone. Additional information regarding geology, morphology, and coral species were also noted.

Hawai'i↗

Evaluation of natural and anthropogenic (human-made) hexavalent chromium

Hexavalent chromium, Cr(VI), was released between 1952 and 1964 from the Pacific Gas and Electric Company (PG&E) Hinkley compressor station, in the Mojave Desert about 80 miles northeast of Los Angeles, California. Geologic, geochemical, and hydrologic data from more than 100 wells collected between March 2015 and November 2017 were interpreted using a summative-scale analysis to define the extent of anthropogenic (human-made) Cr(VI) in groundwater. The summative scale consisted of eight questions requiring binary (yes or no) answers for each sampled well. The questions were intended to (1) provide a transparent framework for data interpretation in which all stakeholders participated; (2) provide unbiased interpretation of data traceable to numerical measurements; (3) provide a framework that enabled geologic, geochemical, and hydrologic data to be considered collectively; and (4) consolidate different types of data into a simple, easy-to-understand interpretation. When data from each well are scored using questions and metrics within the summative scale, all stakeholders would score each well the same way and would draw the same summative-scale Cr(VI) plume extent. The areal extent of the summative-scale Cr(VI) plume was 5.5 square miles (mi 2 ); this is larger than the 2.2-mi 2 extent of the October–December 2015 (Q4 2015) regulatory Cr(VI) plume but smaller than the 8.3-mi2 maximum mapped extent of Cr(VI) greater than the interim regulatory Cr(VI) background value of 3.1 micrograms per liter (μg/L). The summative-scale Cr(VI) plume is within the area covered by the PG&E monitoring well network and lies within “Mojave-type” deposits composed of low-chromium stream and near-shore lake deposits sourced from the Mojave River. The summative-scale Cr(VI) plume included all shallow wells within the footprint of the Q4 2015 regulatory Cr(VI) plume, but summative-scale scores indicate that anthropogenic Cr(VI) was not present in several wells within the footprint of the regulatory Cr(VI) plume that were screened within the deep zone of the upper aquifer. The summative-scale Cr(VI) plume extent was consistent with mineralogic and geochemical data collected as part of this study that were not used within the summative-scale analysis. Data from wells outside the summative-scale Cr(VI) plume collected for regulatory purposes from April 2017 through March 2018 were used to estimate Cr(VI) background concentrations as the upper 95-percent tolerance limit (UTL 95 ) in different parts of Hinkley and Water Valleys. The UTL 95 values were calculated using the computer program ProUCL 5.1 and are suitable for use by regulatory agencies in support of (1) updating the regulatory Cr(VI) plume extent and management of Cr(VI) near the plume margins, (2) establishing cleanup goals for Cr(VI) within the updated regulatory Cr(VI) plume, and (3) identifying unusual Cr(VI) concentrations outside the regulatory Cr(VI) plume. The nonparametric UTL 95 values for wells screened in Mojave-type deposits in the eastern, western, and northern subareas of Hinkley Valley were 3.7, 3.9, and 4.0 μg/L, respectively. The normal UTL 95 values for wells screened in undifferentiated, unconsolidated deposits in the eastern and western subareas and the northern subarea upgradient from the Mount General fault were 2.8, 3.8, and 4.8 μg/L, respectively. An overall normal UTL 95 value of 3.8 μg/L was calculated for undifferentiated, unconsolidated deposits in these areas. This value is similar to the overall nonparametric UTL 95 value of 3.9 μg/L calculated for Mojave-type deposits and similar to the maximum Cr(VI) concentration of older groundwater in contact with Mojave-type deposits of 3.6 μg/L. The provenance of most PG&E monitoring wells is not precisely known, and the UTL 95 values for wells screened in undifferentiated, unconsolidated deposits in the different subareas may be more widely applicable for regulatory purposes than the UTL 95 values for Mojave-type deposits. The UTL 95 value of 2.8 μg/L for wells screened in undifferentiated, unconsolidated deposits in the eastern subarea is important for plume management because most of the summative-scale Cr(VI) plume is within the eastern subarea. A UTL 95 value of 5.8 μg/L was calculated for older (pre-1952) groundwater associated with mudflat/playa deposits in the eastern subarea near Mount General. A UTL 95 value of 2.3 μg/L was calculated for Mojave-type deposits within the Cr(VI) plume downgradient from the Hinkley compressor station after regulatory updates. This lower value is consistent with neutral to slightly alkaline, younger (post-1952) groundwater within coarse-textured, low-chromium Mojave-type deposits in this area and may be a suitable metric for Cr(VI) cleanup goals. The UTL 95 value of 4.8 μg/L for wells screened in undifferentiated, unconsolidated deposits in the northern subarea upgradient from the Mount General fault provides for possible increases in Cr(VI) concentrations if water levels continue to decline. Downgradient from the Q4 2015 regulatory Cr(VI) plume and the summative-scale Cr(VI) plume, UTL 95 values of 9.0 and 6.4 μg/L were calculated for wells screened in undifferentiated, unconsolidated deposits in the northern subarea downgradient from the Mount General fault and for Water Valley, respectively, consistent with different geologic and geochemical conditions in these areas. The UTL 95 values calculated as part of this study provide scientifically defensible estimates of background Cr(VI) concentrations that differ with local geologic, geochemical, and hydrologic conditions in Hinkley and Water Valleys. The regulatory Cr(VI) plume extent can be updated on the basis of these values. The summative-scale Cr(VI) plume extent may contain wells having anthropogenic Cr(VI) concentrations less than the UTL 95 values for their respective subareas that may not require regulatory attention, and an updated regulatory Cr(VI) plume extent may be less than the summative-scale Cr(VI) plume extent. The UTL 95 values are not background Cr(VI) concentrations for regulatory purposes, and the authority to establish regulatory values resides solely with the Lahontan Regional Water Quality Control Board.

California↗

Hydrology and geomorphology of the Taiya River near the West Creek Tributary, southeast Alaska

The Taiya River flows through the Chilkoot Trail Unit of Klondike Gold Rush National Historical Park in southeast Alaska, which was founded to preserve cultural and historical resources and further understanding of natural processes active in the surrounding coastal-to-subarctic basin. Riverine processes exert an important influence on ecologically important boreal toad ( Anaxryus boreas boreas ), salmon [chum salmon ( Oncorhynchus keta ), pink salmon ( O. gorbushca ), and coho salmon ( O. kisutch )], and eulachon ( Thaleichthys pacificus ) habitats, erosion of the historic ghost town of Dyea and other cultural and historical artifacts, and recreational opportunities in the lower 7.5 kilometers (km) of the Taiya River valley bottom. Recurrent consideration of hydroelectric development in West Creek upstream of the park since the 1980s has included proposals for damming and diverting West Creek, which could alter the delivery of water and sediment to this section of the Taiya River. To improve understanding of the hydrologic dependence of park resources for the purposes of guiding effective monitoring and conservation, this study, conducted by the U.S. Geological Survey in cooperation with the National Park Service, used a review of hydrologic data, collection of discrete suspended sediment data, geomorphic mapping, and analysis of historical aerial and ground photographs in a reconnaissance of formative geomorphic processes and hydrologic conditions in the lower 7.5 km of the Taiya River valley bottom. Streamflow and suspended sediment data collected at the U.S. Geological Survey streamgages on the Taiya River and West Creek, combined with historical data, document conditions consistent with streams draining strongly glacierized basins in Alaska. Suspended sediment concentrations from samples collected concurrently over six varying flow levels during 2017–18 ranged from 6 to 284 milligrams per liter (mg/L) for the Taiya River and 13 to 162 mg/L for West Creek, which are similar to or slightly higher than historical values. For the common period of record (1970–77), correlation of daily mean discharge between the two streams was strongest (Pearson’s r = 0.97) during the prolonged May–October high-flow season and weakest ( r = 0.90) during the November–April low-flow season, when West Creek daily mean discharge was proportionally higher. For the Taiya River, streamflow data compared between the available periods of record (1970–77 and 2004–17) showed no decadal-scale patterns in mean annual discharge but did show a shift toward an earlier spring snowmelt pulse. Notable flooding in the Taiya River Basin includes glacial lake outburst floods from the Nourse River valley prior to and during the 1897–98 Gold Rush, a 2002 glacial lake outburst flood from the West Creek valley, and a 1967 rainfall-generated flood. Geomorphic mapping identified four categories of surfaces in the valley bottom—active main stem, abandoned main stem, alluvial fans, and emergent tidal surfaces. Using the maps, main-stem surfaces were subdivided into age categories to identify channel migration patterns from prior to 1940s to 2018. The valley bottom is dominated by active or abandoned channels of the Taiya River except at the extensive low-angle West Creek fan. The active main stem presently supports a mostly single-thread channel with bars and a few sloughs, but the channel actively moved and sometimes was braided within multiple, wider unvegetated corridors in 1894 and earlier. An inventory of 29 off-main-stem channels identified for the study indicates that abandoned main stem channels provide local topographic lows that can intercept groundwater or sustain tributary flow, facilitating the formation of most nonestuarine wetlands in the valley and sustaining important boreal toad breeding habitat. Within the active main stem corridor, the channel has episodically built and reworked meanders and bars, eroding more than one-half of the historic Dyea townsite, in response to glacially controlled delivery of water and sediment, flooding, inputs from West Creek, local features including large woody debris and beaver dams, and rapid uplift from isostatic rebound. West Creek has constructed a large, persistent fan, provoked kilometer-scale Taiya River channel change near the confluence, constructively added to high-season streamflow that affects Taiya River channel migration capacity, disproportionately contributed early-season streamflow, and possibly contributed to groundwater levels in the valley bottom. The progressive narrowing and stability of the main stem corridor, possibly a result of reduction in the magnitude or frequency of glacial lake outburst floods or glacial sediment delivery to streams, indicates less active future reworking of abandoned main-stem surfaces or regeneration of wetland features. The fluvial history of the Taiya River valley bottom collectively indicates continued channel change within a limited corridor, relative stability in wetland locations but uncertainty in stability of groundwater supply to them, and channel incision and extension in response to uplift.

Alaska↗

Historical changes in bed elevation and water depth within the Nehalem Bay, Oregon, 1891–2019

Estuaries, at the nexus of rivers and the ocean, are depositional areas that respond to changes in streamflow, tides, sea level, and inputs of sediment from marine and watershed sources. Understanding changes in bed elevations, deposited and eroded sediment, and water depth throughout estuaries is relevant for understanding their present-day status and long-term evolution, identifying potential hazards to human communities, and informing estuarine conservation. In response to observations of sedimentation in the Nehalem Bay, northwestern Oregon, by the Port of Nehalem, the magnitudes and patterns of bathymetric change in the Bay were documented and described by two approaches. The first approach compared changes in bed elevation with estimated volumes of erosion and deposition from overlapping survey data acquired in 1957 and 2019 for the area of the Nehalem Bay from upstream of the Highway 101 bridge to downstream of Fishery Point. The second approach examined changes in water depth for seven zones from the confluence of the North Fork and Nehalem Rivers to the mouth of the Nehalem River using nautical charts (1891, 1947, 1970, 1990, and 2004). These two approaches were used because the bathymetric surveys from 1957 and 2019 could be tied to a common vertical datum, allowing for a direct comparison of changes in bed elevations, whereas the nautical charts could not be tied to a common vertical datum, which limited the analyses to a comparison of changes in water depths over a broader time frame. Bed elevation changes from 1957 to 2019 were assessed from upstream of the Highway 101 bridge to downstream of Fishery Point where the two surveys overlapped (2 square kilometers) using thalweg longitudinal profiles, channel cross sections, and digital elevation models (DEMs) showing the elevation differences between the two surveys (or DEMs of difference). The most prominent change between 1957 and 2019 was the migration of the thalweg (or deepest part of the channel) between the downstream end of Lazarus Island and downstream of Fishery Point; this migration resulted in sediment deposition in the former thalweg and sediment erosion in formerly shallow areas to form the new thalweg. Bed elevation changes in the thalweg also varied longitudinally between 1957 and 2019. The bed elevation of the thalweg in both surveys, however, was generally less than 1 meter (m). The thalweg in the area of overlapping surveys shortened from about 7.0 to 6.7 kilometers in length over that same period. The bed elevation changes between the DEMs showed that maximum erosion and deposition was 4.3 and 4.5 m, respectively. In this same time period, the net change in sediment volume was 230,000 cubic meters (m 3 ), indicating net deposition. However, the error estimated for the 95 percent confidence interval analyses is ±315,000 m 3 , and therefore does not preclude the possibility that net erosion may have occurred. Historical changes in water depth from soundings depicted on nautical charts from 1891, 1947, 1970, 1990, and 2004 were evaluated by assessing spatial and temporal changes for seven zones of the Nehalem Bay. Across all years and zones, water depths ranged from about 0.2 to 9.4 m, whereas median water depths ranged from 0.3 to 6.4 m. Median depths and the range of water depths did not systematically increase or decrease throughout all zones during the same periods. In all nautical charts, the zone at the mouth of the Nehalem River consistently had the deepest soundings (7.9 to 9.4 m) and the greatest range of water depths (7.3 to 8.8 m). Qualitative evaluation of the nautical charts showed minimal changes in the overall shape of the Nehalem Bay. The exception to this observation was at the mouth of the Bay, where two historical outlets to the Pacific Ocean depicted in the 1891 nautical chart were reduced to one outlet following the construction of jetties (1916 and 1918). The results of this study emphasize that bed elevations and water depths within the Nehalem Bay have varied between 1891 and 2019, as illustrated by the lateral and vertical changes in the thalweg and changes in water depths over time. Changes in thalweg position and related patterns of sediment erosion and deposition are expected in the future as the Nehalem Bay continues to respond to changes in tides, sea level, streamflow, and sediment inputs from watershed and marine sources. The results of this study and the surveys from 1957 and 2019 provide a foundation for documenting and evaluating future changes in the Nehalem Bay and prioritizing actions to manage and protect natural resources and recreational access to the Nehalem Bay.

Oregon↗

Integrated monitoring of hydrogeomorphic, vegetative, and edaphic conditions in riparian ecosystems of Great Basin National Park, Nevada

In semiarid regions such as the Great Basin, riparian areas function as oases of cooler and more stable microclimates, greater relative humidity, greater structural complexity, and a steady flow of water and nutrients relative to upland areas. These qualities make riparian areaʼs attractive not only to resident and migratory wildlife, but also to visitors in recreation areas such as Great Basin National Park in the Snake Range, east-central Nevada. To expand upon the system of ten permanent plots sampled in 1992 (Smith et al. 1994) and 2001 (Beever et al. in press ), we established a collection of 31 cross-sectional transects of 50-m width across the mainstems of Strawberry, Lehman, Baker, and Snake creeks. Our aims in this research were threefold: a) map riparian vegetative communities in greater detail than had been done by past efforts; b) provide a monitoring baseline of hydrogeomorphology; structure, composition, and function of upland- and riparianassociated vegetation; and edaphic properties potentially sensitive to management; and c) test whether instream conditions or physiographic variables predicted vegetation patterns across the four target streams. In each of the four watersheds, we performed walking transects from the lower-elevation boundary of the park along creek mainstems to a point well above the point at which vehicle access stopped. In these transects, we ranked, by cover, the riparian and upland woody species on each side of the creek, in 0.32-km segments. These walking transects also facilitated selection of a suite of cross-sectional transects that might serve as an early-warning signal of change for natural (e.g., aggradative) and anthropogenic changes (e.g., due to visitor impacts or climate change). At each cross-sectional transect, we used several methods: a) measurement of the number, approximate volume, and total length of instream logs greater than 10 cm in diameter that were within 5 m up- or downstream of the transect; b) counts of pebbles by size class, following Wolman (1954); c) line-point intercepts, which provided various measures of percent cover; d) gap-intercept transects, following Herrick et al. (in press), to measure susceptibility of uplands to erosion by wind or water; e) 1-m 2 quadrats, to obtain frequency of woody species; f) nested-frequency plots, to measure frequency of all plant species in quadrats of varying size; g) a field-based soil aggregate stability test following Herrick et al. (2001); and h) an impact penetrometer, to measure penetration resistance of soil horizons. We used species-accumulation curves to assess the ability of our methods to detect the majority of plant species at sites, using the most species-rich and species-poor sites as illustrations. We compared characteristics of hydrogeomorphic valley types (designated by Frissell and Liss 1993), vegetation types, and creeks individually and, using multivariate analyses for the first two ʻtypes,ʼ simultaneously. For the latter, using both the nested-frequency and 1-m 2 frequency data, we first used nonmetric multidimensional scaling (NMS) to assess relationships of plant communities among sites. Secondly, we used multi-response permutation procedures (MRPP) to test whether plant-community differences existed among either hydrogeomophic valley types or vegetation types. To increase the value of these comparisons for management, we used indicator species analyses to quantify the indicator value of each individual plant species for separating groups. In contrast to the more incised riparian channels of central Nevada, we observed knickzones, downcutting, and incision only rarely and usually with limited extent in the walking surveys. Downcutting occurred most frequently and extensively in Strawberry and Snake creeks, due in part to their more erodible soils. According to a hydrogeomorphologist with extensive experience in Great Basin riparian systems, the sediment-delivery and hydrologic systems appeared relatively undisturbed in most reaches, with respect to grazing animals and other types of anthropogenic alteration. Site elevation of the 31 transects ranged from 1,950-2,987 m, and stream slope (i.e., gradient) was relatively steep (mean = 9.3%, range 3-16%). Strawberry Creek averaged the lowest maximum water depth, and correspondingly had greatest width/depth ratios. Baker Creek sites averaged the smallest amount of tree-canopy gaps, whereas Snake Creek sites on average had the largest proportion of gaps in understory vegetation. Sites in terrace-bound valley types averaged the lowest slope in the channel as well as the least cover of trees, litter, and vegetation overall, whereas alluviated, boulder-bed canyon sites averaged the greatest widths of the active channel. Sites in Lehman Creek averaged nearly twice as much coarse woody debris as sites from any other creek, whereas Baker Creek sites averaged greatest tree cover (mean = 67%, range 40 – 96%) and species richness (mean = 17.3 species). Multivariate ordinations suggested that sites in leveed outwash valleys and alluvial-fan-influenced valleys had the greatest inter-site heterogeneity in plant composition, whereas sites in incised moraine-filled valleys appeared most homogeneous. Differences among homogeneity of sites within vegetation types were less pronounced, but sites dominated by either aspen and Woodsʼ rose or narrow-leaved cottonwood had the most similar plant communities among sites of the same vegetation type. A number of species were faithful indicators of various valley and vegetation types, using either set of plant-frequency data. We estimate that all 31 sites could be subsequently re-sampled in 14-18 field days by individuals possessing familiarity of the riparian flora of the southern Snake Range. As with any research, monitoring-focused investigations must balance the concerns for number of ecosystem attributes measured, extensiveness in time and space of sampling periods and locations, and the time and cost of sampling.

Nevada↗

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

In this report, SPAtially Referenced Regression On Watershed attributes (SPARROW) models developed to describe long-term (2000–14) mean-annual streamflow, total nitrogen (TN), total phosphorus (TP), and suspended-sediment (SS) transport in streams of the Midwestern part of the United States (the Mississippi River, Great Lakes, and Red River of the North Basins) are described. The nutrient and suspended-sediment models have a base year of 2012, which means they were developed based on source inputs and management practices similar to those existing during or near 2012 and average hydrological conditions detrended to 2012 (2000–14), whereas the streamflow model has base years of 2000–14, which means it was developed based on the average input precipitation minus actual evapotranspiration from 2000 to 2014. In developing the models, several updates and improvements were made to the data inputs and statistical approaches used to calibrate/develop the models from those used in the previous 2002 SPARROW models. The 2012 SPARROW models were constructed using a higher resolution stream network, which resulted in a mean catchment size of 2.7 square kilometers compared to 480 square kilometers in the 2002 models; more detailed and updated wastewater treatment plant contribution estimates; inputs from background phosphorus sources that were not included in the 2002 model; and more accurate loads for calibration that were computed using a modified Beale ratio-estimator technique whenever no trend in load was determined. Statistical approaches were added to compensate for the unequal effect of each monitoring site during the calibration process by adjusting for the fraction of the basin included in other upstream monitored sites (nested share) and thinning the calibration sites if a negative statistical correlation between nearby sites was determined. Results from 2012 SPARROW models describe how much of each water, TN, TP, and SS source was delivered to the stream network, and the major landscape factors that affected their delivery. Atmospheric deposition and natural (background) sources of TN and TP, respectively, were the dominant sources in anthropogenically unaffected areas (especially in the Rocky Mountains and north-central areas of the Midwest), whereas fertilizers, manure, and fixation were dominant sources in agricultural areas, especially in the Corn Belt and near the Mississippi River. Urban sources of TN and TP were typically localized, but they were still important for some large areas, especially the Lake Erie Basin. All of the land-to-water delivery variables in the nutrient and sediment SPARROW models, such as runoff, soil erodibility, basin slope, and the amount of tile drains, are commonly included in process-driven models. In the SPARROW TN and TP models, best management practices (BMPs) reduced the delivery of these nutrients to streams. Long-term mean-annual flows and nutrient and sediment loads were simulated in streams throughout the Midwest. The simulated flows from the SPARROW flow model were used in the SPARROW TN, TP, and SS models to help describe nutrient and sediment transport from the watershed and through the stream network. Outputs from the TN, TP, and SS models describe loads and yields of these constituents throughout the Midwest, and from major drainage basins throughout the Midwest. Highest TN, TP, and SS yields and delivered yields were from the Lake Erie, Ohio River, Upper Mississippi River, and Lower Mississippi River Basins, whereas lowest yields were spread over most other areas. Losses during downstream delivery resulted in part of the TN, TP, and SS that reach the stream network not reaching the downstream receiving bodies: 14, 15, and 28 percent of the TN, TP, and SS, respectively, are lost during delivery to the Great Lakes and 19, 23, and 52 percent of the TN, TP, and SS, respectively, are lost during delivery to the Gulf of Mexico. The largest losses of nutrients and sediments during transport were in the Missouri and Arkansas River Basins. Information from these SPARROW models can help guide nutrient and sediment reduction strategies throughout the Midwest. Model results provide information on what may be the most appropriate general type of actions to reduce total loading by describing the relative importance of each source, and where to most efficiently place the efforts to reduce loading by describing the distribution of nutrient and sediment loading. By implementing management efforts addressing the major sources of the loads in areas contributing the highest loads, it may be possible to reduce nutrient loading throughout the Mississippi River Basin and thus reduce the size of the hypoxic zone in the Gulf of Mexico; reduce nutrient loading into lakes, and thus reduce the occurrence of harmful algal blooms; and reduce sediment losses, and thus improve the benthic habitat in streams and rivers throughout the Midwest.

Scientific Investigations Report↗

Evaluation of the characteristics, discharge, and water quality of selected springs at Fort Irwin National Training Center, San Bernardino County, California

Eight springs and seeps at Fort Irwin National Training Center were described and categorized by their general characteristics, discharge, geophysical properties, and water quality between 2015 and 2017. The data collected establish a modern (2017) baseline of hydrologic conditions at the springs. Two types of springs were identified: (1) precipitation-fed upland springs (Cave, Desert King, Devouge, No Name, and Panther Springs) and (2) groundwater discharge-fed basin springs (Garlic, Bitter, and Jack Springs). Comparison of electrical resistivity tomography data collected at groundwater basin springs from 2015 to 2017 indicated that spring discharge and connection to the underlying groundwater system is highly focused, although the springs themselves appear diffuse and are spread out over a large area. Spring discharge was consistently less than reported by Thompson (1929), except at Garlic Spring where discharges and vegetation have increased in recent years. Multiple discrete flume and seepage meter measurements taken between October 2015 and April 2016 indicated that discharge changed predictably on diurnal and seasonal timescales in response to evapotranspiration. These preliminary results and the lush vegetation noted at some of the springs, particularly at Bitter, Garlic, and Jack Springs, indicated plant evapotranspiration accounts for a substantial part of the discharge from these springs. The quality of water ranges from fresh in precipitation-fed upland springs (Cave, Desert King, Devouge, and Panther Springs) to slightly saline (Garlic and Jack Springs) and moderately saline (Bitter Spring) in groundwater-fed discharge springs. Nitrate concentrations from water at most of the springs were less than 3 milligrams per liter, except for samples from Devouge and Desert King Springs and one sample from Jack Spring. An analysis of delta nitrogen-15 in nitrate (δ 15 N-NO 3 ) and delta oxygen-18 in nitrate (δ 18 O-NO 3 ) indicates high nitrate concentrations in excess of the U.S. Environmental Protection Agency maximum contaminant level at Jack Spring and Desert King Spring resulting from the dissolution of nitrate-bearing caliche deposits; nitrate concentrations at Devouge Spring are a result of algal growth within the spring, and the source of nitrate concentrations in Garlic Spring are consistent with a treated wastewater origin from Langford Valley-Irwin subbasin upgradient. The source of water in upland springs, indicated by values of delta oxygen-18 (δ 18 O) and delta deuterium (δD) are consistent with recharge from winter precipitation. In groundwater basin springs, values of δ 18 O and δD are consistent with groundwater sampled from nearby wells. Summer monsoonal precipitation appears to contribute little water to spring flow. Most springs contain low levels of tritium and appear to be primarily older (pre-1950s) groundwater. Groundwater basin springs with detectable tritium may result from occasional streamflow in nearby washes. These springs could be susceptible to decreases in flow during extended dry periods when the localized recharge may be reduced due to the loss of focused recharge through nearby washes. Groundwater samples from Garlic and Bitter Springs contained arsenic concentrations above the U.S. Environmental Protection Agency maximum contaminant level. Groundwater samples from all springs, except Cave, Desert King, and Devouge Springs, exceeded the State of California maximum contaminant level for fluoride. Garlic Spring was the only sampled spring that contained vanadium concentrations that exceeded the State of California notification level. Only a single water sample from Jack Spring contained uranium at a concentration that exceeded the U.S. Environmental Protection Agency maximum contaminant level. Many other constituents of concern were analyzed, including those from anthropogenic sources that may be a result of military activities. Most of these constituents were not detected above their respective reporting levels in spring water; only 15 were detected in spring waters. Diesel and gasoline degradants, many of which also occur naturally, were the most commonly detected compounds. Several other organic compounds, primarily solvents or their degradants, were detected in groundwater basin springs. These constituents, in order of decreasing detection frequency, were carbon disulfide; perchlorate; mercury; acetone; methylnaphthalene; toluene; methyl ethyl ketone; cyanide; and styrene; 4-iso-propyl-toluene; isopropylbenzene; methyl salicylate; and phenol. Except for Garlic Spring, which is affected by discharges of treated wastewater, the quality of water from most springs appears to be relatively unaffected by activities at the Fort Irwin National Training Center.

California↗

Hexavalent and total chromium at low reporting concentrations in source-water aquifers and surface waters used for public supply in Illinois, 2013

On the basis of their recent review of the human health effects of hexavalent chromium [Cr(VI)] in public drinking water, the U.S. Environmental Protection Agency is considering the need for Federal regulation of Cr(VI). Presently, only total chromium is regulated, at a Maximum Contaminant Level (MCL) of 100 micrograms per liter (µg/L). The occurrence of Cr(VI) in groundwater and surface waters generally is attributed to industrial sources, but can be of natural origin. California’s recently established MCL for Cr(VI) of 10 µg/L illustrates the drinking-water concerns associated with Cr(VI). To improve understanding of the possible impact of a Cr(VI)-specific standard that approximates the California level on the management of Illinois’ public drinking water, the U.S. Geological Survey, in cooperation with the Illinois Environmental Protection Agency, assessed the occurrence and distribution of Cr(VI) in the State’s public-water supplies. During 2013, untreated water samples were collected to be analyzed for Cr(VI) and total chromium [Cr(T)] at 119 water-supply wells and 32 surface-water intakes; also, 32 treated surface-water samples were collected near the point of treatment and 32 near the furthest point of distribution. Public-supply sample sites were selected by a stratified random method. Samples typically were analyzed within 24 hours of collection at reporting limits of 0.02 µg/L for Cr(VI) and 0.1 µg/L for Cr(T). The occurrence of Cr(VI) was compared with selected geophysical, physical, and sampling factors that might more fully explain its distribution and magnitude of concentrations. The maximum concentration of Cr(VI) in groundwater was 2.1 µg/L. Maximum concentrations in untreated and treated surface water were 0.29 µg/L and 2.4 µg/L, respectively. All sample concentrations were below the California MCL; only 35 percent were below that State’s non-enforceable public health goal of 0.02 of µg/L. Cr(VI) was undetected in 43 percent of untreated groundwater samples, with a median of 0.06 µg/L when detected. All but two (94 percent) of untreated surface-water samples had detections. In untreated surface water, the median concentration was 0.09 µg/L, whereas in treated (tap and distributed) water the median was 0.20 µg/L. Surface waters treated with lime for softening typically had the greatest Cr(VI) concentrations (maximum, 2.4 µg/L; median, 1.2 µg/L). The maximum concentration of Cr(T) in groundwater was 1.8 µg/L. Maximum concentrations in untreated and treated surface water were 1.8 µg/L and 2.5 µg/L, respectively. All sample concentrations were below the Federal MCL. Total chromium was detected in 65 percent of untreated groundwater samples, with a median of 0.40 µg/L, when detected. All but one (97 percent) of untreated surface-water samples had detections. In untreated surface water, the median concentration was 0.40 µg/L, whereas in treated (tap and distributed) water the median was 0.30 µg/L. As with Cr(VI), surface waters treated with lime typically had the greatest Cr(T) concentrations. Examination of factors that might account for or be associated with the occurrence of Cr(VI) in public-supply source waters found few clearly evident factors. Associations in frequencies of occurrence and range of concentrations indicate that surface waters and groundwaters of shallow, unconsolidated, unconfined aquifers, particularly alluvial aquifers, are possibly most commonly affected by anthropogenic sources of Cr(VI). Groundwaters of deep (greater than 500 feet) bedrock aquifers, particularly the Cambrian-Ordovician aquifer system, are possibly most commonly affected by geologic sources of Cr(VI). Additional study, with supporting geologic and geochemical data that were not collected in this study, would be necessary to verify these associations. There was a weak positive relation (ρ = 0.23) between concentrations of Cr(VI) and Cr(T) in untreated water samples, with a much stronger positive relation (ρ = 0.86 and ρ = 0.90, respectively) in samples collected soon after treatment and near the endpoint of distribution. The stronger relation and greater similarity between Cr(VI) and Cr(T) concentrations in treated water samples indicate that Cr(VI) represents a greater proportion of the measured concentrations of Cr(T) in treated waters than in untreated waters. The analysis of spikes and other quality-assurance samples indicate uncertainties associated with obtaining or confirming consistently accurate analytical results for Cr(VI) at near the applied reporting limit of 0.02 µg/L.

Illinois↗

Application of dimensionless sediment rating curves to predict suspended-sediment concentrations, bedload, and annual sediment loads for rivers in Minnesota

Consistent and reliable sediment data are needed by Federal, State, and local government agencies responsible for monitoring water quality, planning river restoration, quantifying sediment budgets, and evaluating the effectiveness of sediment reduction strategies. Heightened concerns about excessive sediment in rivers and the challenge to reduce costs and eliminate data gaps has guided Federal and State interests in pursuing alternative methods for measuring suspended and bedload sediment. Simple and dependable data collection and estimation techniques are needed to generate hydraulic and water-quality information for areas where data are unavailable or difficult to collect. The U.S. Geological Survey, in cooperation with the Minnesota Pollution Control Agency and the Minnesota Department of Natural Resources, completed a study to evaluate the use of dimensionless sediment rating curves (DSRCs) to accurately predict suspended-sediment concentrations (SSCs), bedload, and annual sediment loads for selected rivers and streams in Minnesota based on data collected during 2007 through 2013. This study included the application of DSRC models developed for a small group of streams located in the San Juan River Basin near Pagosa Springs in southwestern Colorado to rivers in Minnesota. Regionally based DSRC models for Minnesota also were developed and compared to DSRC models from Pagosa Springs, Colorado, to evaluate which model provided more accurate predictions of SSCs and bedload in Minnesota. Multiple measures of goodness-of-fit were developed to assess the effectiveness of DSRC models in predicting SSC and bedload for rivers in Minnesota. More than 600 dimensionless ratio values of SSC, bedload, and streamflow were evaluated and delineated according to Pfankuch stream stability categories of “good/fair” and “poor” to develop four Minnesota-based DSRC models. The basis for Pagosa Springs and Minnesota DSRC model effectiveness was founded on measures of goodness-of-fit that included proximity of the model(s) fitted line to the 95-percent confidence intervals of the site-specific model, Nash-Sutcliffe Efficiency values, model biases, and deviation of annual sediment loads from each model to the annual sediment loads calculated from measured data. Composite plots comparing Pagosa Springs DSRCs, Minnesota DSRCs, site-specific regression models, and measured data indicated that regionally developed DSRCs (Minnesota DSRC models) more closely approximated measured data for nearly every site. Pagosa Springs DSRC models had markedly larger exponents (slopes) when compared to the Minnesota DSRC models and the site-specific regression models, and over-represent SSC and bedload at streamflows exceeding bankfull. The Nash-Sutcliffe Efficiency values for the Minnesota DSRC model for suspended-sediment concentrations closely matched Nash-Sutcliffe Efficiency values of the site-specific regression models for 12 out of 16 sites. Nash-Sutcliffe Efficiency values associated with Minnesota DSRCs were greater than those associated with Pagosa Springs DSRCs for every site except the Whitewater River near Beaver, Minnesota site. Pagosa Springs DSRC models were less accurate than the mean of the measured data at predicting SSC values for one-half of the sites for good/fair stability sites and one-half of the sites for poor stability sites. Relative model biases were calculated and determined to be substantial (greater than 5 percent) for Pagosa Springs and Minnesota models, with Minnesota models having a lower mean model bias. For predicted annual suspended-sediment loads (SSL), the Minnesota DSRC models for good/fair and poor stream stability sites more closely approximated the annual SSLs calculated from the measured data as compared to the Pagosa Springs DSRC model. Results of data analyses indicate that DSRC models developed using data collected in Minnesota were more effective at compensating for differences in individual stream characteristics across a variety of basin sizes and flow regimes than DSRC models developed using data collected for Pagosa Springs, Colorado. Minnesota DSRC models retained a substantial portion of the unique sediment signatures for most rivers, although deviations were observed for streams with limited sediment supply and for rivers in southeastern Minnesota, which had markedly larger regression exponents. Compared to Pagosa Springs DSRC models, Minnesota DSRC models had regression slopes that more closely matched the slopes of site-specific regression models, had greater Nash-Sutcliffe Efficiency values, had lower model biases, and approximated measured annual sediment loads more closely. The results presented in this report indicate that regionally based DSRCs can be used to estimate reasonably accurate values of SSC and bedload. Practitioners are cautioned that DSRC reliability is dependent on representative measures of bankfull streamflow, SSC, and bedload. It is, therefore, important that samples of SSC and bedload, which will be used for estimating SSC and bedload at the bankfull streamflow, are collected over a range of conditions that includes the ascending and descending limbs of the event hydrograph. The use of DSRC models may have substantial limitations for certain conditions. For example, DSRC models should not be used to predict SSC and sediment loads for extreme streamflows, such as those that exceed twice the bankfull streamflow value because this constitutes conditions beyond the realm of current (2016) empirical modeling capability. Also, if relations between SSC and streamflow and between bedload and streamflow are not statistically significant, DSRC models should not be used to predict SSC or bedload, as this could result in large errors. For streams that do not violate these conditions, DSRC estimates of SSC and bedload can be used for stream restoration planning and design, and for estimating annual sediment loads for streams where little or no sediment data are available.

Minnesota↗

Ecological assessment of streams in the Powder River Structural Basin, Wyoming and Montana, 2005-06

Energy and mineral development, particularly coalbed natural gas development, is proceeding at a rapid pace in the Powder River Structural Basin (PRB) in northeastern Wyoming. Concerns about the potential effects of development led to formation of an interagency working group of primarily Federal and State agencies to address these issues in the PRB in Wyoming and in Montana where similar types of resources exist but are largely undeveloped. Under the direction of the interagency working group, an ecological assessment of streams in the PRB was initiated to determine the current status (2005–06) and to establish a baseline for future monitoring. The ecological assessment components include assessment of stream habitat and riparian zones as well as assessments of macroinvertebrate, algal, and fish communities. All of the components were sampled at 47 sites in the PRB during 2005. A reduced set of components, consisting primarily of macroinvertebrate and fish community assessments, was sampled in 2006. Related ecological data, such as habitat and fish community data collected from selected sites in 2004, also are included in this report. The stream habitat assessment included measurement of channel features, substrate size and embeddedness, riparian vegetation, and reachwide characteristics. The width-to-depth ratio (bankfull width/bankfull depth) tended to be higher at sites on the main-stem Powder River than at sites on the main-stem Tongue River and at sites on tributary streams. The streambed substrate particle size was largest at sites on the main-stem Tongue River and smallest at sites on small tributary streams such as Squirrel Creek and Otter Creek. Total vegetative cover at the ground level, understory, and canopy layers ranged from less than 40 percent at a few sites to more than 90 percent at many of the sites. A bank-stability index indicated that sites in the Tongue River drainage were less at risk of bank failure than sites on the main-stem Powder River. Macroinvertebrate communities showed similarity at the river-drainage scale. Macroinvertebrate communities at sites with mountainous headwaters and snowmelt-driven hydrology, such as Clear Creek, Crazy Woman Creek, and Goose Creek, showed similarity with communities from the main-stem Tongue River. The data also indicated similarity among sites on the main-stem Powder River and among small tributaries of the Tongue River. Data analyses using macroinvertebrate observed/expected models and multimetric indices developed by the States of Wyoming and Montana indicated a tendency toward declining biological condition in the downstream direction along the Tongue River. Biological condition for the main-stem Powder River generally improved downstream, from below Salt Creek to near the Wyoming/Montana border, followed by a general decline downstream from the border to the confluence with the Yellowstone River. The biological condition generally was not significantly different between 2005 and 2006, although streamflow was less in 2006 because of drought. Algal communities showed similarity at the river-drainage scale with slight differences from the pattern observed in the macroinvertebrate communities. Although the algal communities from Clear Creek and Goose Creek were similar to those from the main-stem Tongue River, as was true of the macroinvertebrate communities, the algal communities from Crazy Woman Creek had more similarity to those of main-stem Powder River sites than to the Tongue River sites, contrary to the macroinvertebrates. Ordination of algal communities, as well as diatom metrics including salinity and dominant taxa, indicated substantial variation at two sites along the main stem of the Powder River. Fish communities of the PRB were most diverse in the Tongue River drainage. In part due to the effects of Tongue River Reservoir, 15 species of fish were found in the Tongue River drainage that were not found in the Cheyenne, Belle Fourche, or Little Powder River drainages. The number of introduced species and relative abundance of introduced species of fish were higher in the Tongue River and other drainages than at sites on the main-stem Powder River. Although non-native species were identified in the Powder River, the native fish community is largely intact. Western silvery minnow and sturgeon chub—species of special concern—were identified only at sites on the main-stem Powder River and were most common in the Montana segment of the main stem. Fish and habitat sampling on the main-stem Powder River indicated affinity of some species for certain habitats such as pools, runs, riffles, backwaters, or shoals.

Montana, Wyoming↗

Status of groundwater quality in the Upper Santa Ana Watershed, November 2006--March 2007--California GAMA Priority Basin Project

Groundwater quality in the approximately 1,000-square-mile (2,590-square-kilometer) Upper Santa Ana Watershed (USAW) study unit was investigated as part of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is located in southern California in Riverside and San Bernardino Counties. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey and the Lawrence Livermore National Laboratory. The GAMA USAW study was designed to provide a spatially unbiased assessment of untreated groundwater quality within the primary aquifer systems in the study unit. The primary aquifer systems (hereinafter, primary aquifers) are defined as the perforation interval of wells listed in the California Department of Public Health (CDPH) database for the USAW study unit. The quality of groundwater in shallower or deeper water-bearing zones may differ from that in the primary aquifers; shallower groundwater may be more vulnerable to surficial contamination. The assessment is based on water-quality and ancillary data collected by the U.S. Geological Survey (USGS) from 90 wells during November 2006 through March 2007, and water-quality data from the CDPH database. The status of the current quality of the groundwater resource was assessed based on data from samples analyzed for volatile organic compounds (VOCs), pesticides, and naturally occurring inorganic constituents, such as major ions and trace elements. The status assessment is intended to characterize the quality of groundwater resources within the primary aquifers of the USAW study unit, not the treated drinking water delivered to consumers by water purveyors. Relative-concentrations (sample concentration divided by the health- or aesthetic-based benchmark concentration) were used for evaluating groundwater quality for those constituents that have Federal or California regulatory or non-regulatory benchmarks for drinking-water quality. A relative-concentration greater than (>) 1.0 indicates a concentration above a benchmark, and a relative-concentration less than or equal to (&le;) 1.0 indicates a concentration equal to or less than a benchmark. Organic and special-interest constituent relative-concentrations were classified as "high" (> 1.0), "moderate" (0.1 < relative-concentration &le; 1.0), or "low" (&le; 0.1). Inorganic constituent relative-concentrations were classified as "high" (> 1.0), "moderate" (0.5 < relative-concentration &le; 1.0), or "low" ( &le; 0.5). Aquifer-scale proportion was used as the primary metric in the status assessment for evaluating regional-scale groundwater quality. Aquifer-scale proportions are defined as the percentage of the area of the primary aquifer system with concentrations above or below specified thresholds relative to regulatory or aesthetic benchmarks. High aquifer-scale proportion is defined as the percentage of the area of the primary aquifers with a relative-concentration greater than 1.0 for a particular constituent or class of constituents; percentage is based on an areal, rather than a volumetric basis. Moderate and low aquifer-scale proportions were defined as the percentage of the primary aquifers with moderate and low relative-concentrations, respectively. Two statistical approaches&mdash;grid-based and spatially weighted&mdash;were used to evaluate aquifer-scale proportions for individual constituents and classes of constituents. Grid-based and spatially weighted estimates were comparable in the USAW study unit (within 90-percent confidence intervals). Inorganic constituents with human-health benchmarks had relative-concentrations that were high in 32.9 percent of the primary aquifers, moderate in 29.3 percent, and low in 37.8 percent. The high aquifer-scale proportion of these inorganic constituents primarily reflected high aquifer-scale proportions of nitrate (high relative-concentration in 25.3 percent of the aquifer), although seven other inorganic constituents with human-health benchmarks also were detected at high relative-concentrations in some percentage of the aquifer: arsenic, boron, fluoride, gross alpha activity, molybdenum, uranium, and vanadium. Perchlorate, as a constituent of special interest, was evaluated separately from other inorganic constituents, and had high relative-concentrations in 11.1 percent, moderate in 53.3 percent, and low or not detected in 35.6 percent of the primary aquifers. In contrast to the inorganic constituents, relative-concentrations of organic constituents (one or more) were high in 6.7 percent, moderate in 11.1 percent, and low or not detected in 82.2 percent of the primary aquifers. Of the 237 organic and special-interest constituents analyzed for, 39 constituents were detected (21 VOCs, 13 pesticides, 3 pharmaceuticals, and 2 constituents of special interest). All of the detected VOCs had health-based benchmarks, and five of these&mdash;1,1-dichloroethene, 1,2-dibromo-3-chloropropane (DBCP), tetrachloroethene (PCE), carbon tetrachloride, and trichloroethene (TCE)&mdash;were detected in at least one sample at a concentration above a benchmark (high relative-concentration). Seven of the 13 pesticides had health-based benchmarks, and none were detected above these benchmarks (no high relative-concentrations). Pharmaceuticals do not have health-based benchmarks. Thirteen organic constituents were frequently detected (detected in at least 10 percent of samples without regard to relative-concentrations): bromodichloromethane, chloroform, cis-1,2-dichloroethene, 1,1-dichloroethene, dichlorodifluoromethane (CFC-12), methyl tert-butyl ether (MTBE), PCE, TCE, trichlorofluoromethane (CFC-11), atrazine, bromacil, diuron, and simazine.

California↗

Quality of streams in Johnson County, Kansas, 2002--10

Stream quality in Johnson County, northeastern Kansas, was assessed on the basis of land use, hydrology, stream-water and streambed-sediment chemistry, riparian and in-stream habitat, and periphyton and macroinvertebrate community data collected from 22 sites during 2002 through 2010. Stream conditions at the end of the study period are evaluated and compared to previous years, stream biological communities and physical and chemical conditions are characterized, streams are described relative to Kansas Department of Health and Environment impairment categories and water-quality standards, and environmental factors that most strongly correlate with biological stream quality are evaluated. The information is useful for improving water-quality management programs, documenting changing conditions with time, and evaluating compliance with water-quality standards, total maximum daily loads (TMDLs), National Pollutant Discharge Elimination System (NPDES) permit conditions, and other established guidelines and goals. Constituent concentrations in water during base flow varied across the study area and 2010 conditions were not markedly different from those measured in 2003, 2004, and 2007. Generally the highest specific conductance and concentrations of dissolved solids and major ions in water occurred at urban sites except the upstream Cedar Creek site, which is rural and has a large area of commercial and industrial land less than 1 mile upstream on both sides of the creek. The highest base-flow nutrient concentrations in water occurred downstream from wastewater treatment facilities. Water chemistry data represent base-flow conditions only, and do not show the variability in concentrations that occurs during stormwater runoff. Constituent concentrations in streambed sediment also varied across the study area and some notable changes occurred from previously collected data. High organic carbon and nutrient concentrations at the rural Big Bull Creek site in 2003 decreased to at least one-fourth of those concentrations in 2007 and 2010 likely because of the reduction in upstream wastewater discharge contributions. The highest concentrations of trace metals in 2010 occurred at urban sites on Mill and Indian Creeks. Zinc was the only metal to exceed the probable effects concentration in 2010, which occurred at a site on Indian Creek. In 2007, chromium and nickel at the upstream urban Cedar Creek site exceeded the probable effects concentrations, and in 2003, no metals exceeded the probable effects concentrations. Of 72 organic compounds analyzed in streambed sediment, 26 were detected including pesticides, polycyclic aromatic hydrocarbons (PAHs), fuel products, fragrances, preservatives, plasticizers, manufacturing byproducts, flame retardants, and disinfectants. All 6 PAH compounds analyzed were detected, and the probable effects concentrations for 4 of the 6 PAH compounds analyzed were exceeded in 2010. Only five pesticide compounds were detected in streambed sediment, including carbazole and four pyrethroid compounds. Chronic toxicity guidelines for pyrethroid compounds were exceeded at five sites. Biological conditions reflected a gradient in urban land use, with the less disturbed streams located in rural areas of Johnson County. About 19 percent of sites in 2010 (four sites) were fully supporting of aquatic life on the basis of the four metrics used by Kansas Department of Health and Environment to categorize sites. This is a notable difference compared to previous years when no sites (in 2003 and 2004) or just one site (in 2007) was fully supporting of aquatic life. Multimetric macroinvertebrate scores improved at the Big Bull Creek site where wastewater discharges were reduced in 2007. Environmental variables that consistently were highly negatively correlated with biological conditions were percent impervious surface and percent urban land use. In addition, density of stormwater outfall points adjacent to streams was significantly negatively correlated with biological conditions. Specific conductance of water and sum of PAH concentrations in streambed sediment also were significantly negatively correlated with biological conditions. Total nitrogen in water and total phosphorus in streambed sediment were correlated with most of the invertebrate variables, which is a notable difference from previous analyses using smaller datasets, in which nutrient relations were weak or not detected. The most important habitat variables were sinuosity, length and continuity of natural buffers, riffle substrate embeddedness, and substrate cover diversity, each of which was correlated with all invertebrate metrics including a 10-metric combined score. Correlation analysis indicated that if riparian and in-stream habitat conditions improve then so might invertebrate communities and stream biological quality. Sixty-two percent of the variance in macroinvertebrate community metrics was explained by the single environmental factor, percent impervious surface. Invertebrate responses to urbanization in Johnson County indicated linearity rather than identifiable thresholds. Multiple linear regression models developed for each of the four macroinvertebrate metrics used to determine aquatic-life-support status indicated that percent impervious surface, as a measure of urban land use, explained 34 to 67 percent of the variability in biological communities. Results indicate that although multiple factors are correlated with stream quality degradation, general urbanization, as indicated by impervious surface area or urban land use, consistently is determined to be the fundamental factor causing change in stream quality. Effects of urbanization on Johnson County streams are similar to effects described in national studies that assess effects of urbanization on stream health. Individually important environmental factors such as specific conductance of water, PAHs in streambed sediment, and stream buffer conditions, are affected by urbanization and, collectively, all contribute to stream impairments. Policies and management practices that may be most important in protecting the health of streams in Johnson County are those minimizing the effects of impervious surface, protecting stream corridors, and decreasing the loads of sediment, nutrients, and toxic chemicals that directly enter streams through stormwater runoff and discharges.

Kansas↗

Aquatic biological communities and associated habitats at selected sites in the Big Wood River Watershed, south-central Idaho, 2014

Assessments of streamflow (discharge) parameters, water quality, physical habitat, and biological communities were completed between May and September 2014 as part of a monitoring program in the Big Wood River watershed of south-central Idaho. The sampling was conducted by the U.S. Geological Survey in cooperation with Blaine County, Trout Unlimited, the Nature Conservancy, and the Wood River Land Trust to help identify the status of aquatic resources at selected locations in the watershed. Information in this report provides a basis with which to evaluate and monitor the long-term health of the Big Wood River and its major tributaries. Sampling sites were co-located with existing U.S. Geological Survey streamgaging stations: three on the main stem Big Wood River and four on the North Fork Big Wood River (North Fork), Warm Springs Creek (Warm Sp), Trail Creek (Trail Ck), and East Fork Big Wood River (East Fork) tributaries. The analytical results and quality-assurance information for water quality, physical habitat, and biological community samples collected at study sites during 2 weeks in September 2014 are summarized. Water-quality data include concentrations of major nutrients, suspended sediment, dissolved oxygen, and fecal-coliform bacteria. To assess the potential effects of nutrient enrichment on algal growth, concentrations of periphyton biomass (chlorophyll- a and ash free dry weight) in riffle habitats were determined at each site. Physical habitat parameters include stream channel morphology, habitat volume, instream structure, substrate composition, and riparian vegetative cover. Biological data include taxa richness, abundance, and stream-health indicator metrics for macroinvertebrate and fish communities. Statistical summaries of the water-quality, habitat, and biological data are provided along with discussion of how these findings relate to the health of aquatic resources in the Big Wood River watershed. Seasonal discharge patterns using statistical summaries of daily discharge from selected sites are reported for water years 2012–15. Results showed that annual average daily mean discharge increased from the Big Wood River near Ketchum, ID (BW Ketchum) downstream to the Big Wood River at Hailey, ID (BW Hailey), but decreased by nearly 50 percent from BW Hailey downstream to Big Wood River at Stanton Crossing near Bellevue, ID (BW Stanton). Annual variability in daily mean discharge among main-stem sites was highest at BW Stanton, suggesting that this part of the river may be subject to some level of flow alteration. Hydrologic alterations resulting in flow reduction can contribute to higher water temperature, especially during the summer months when conditions are often most stressful to fish and other stream organisms. Daily water temperature and water temperature trends from June to September 2014 are reported for select tributary and main-stem sites on the Big Wood River and can be used to assess the potential for biological impairment based on aquatic life temperature criteria for cold-water streams. The State of Idaho maximum temperature criteria for protection of cold-water aquatic life of 22 °C was exceeded at Warm Sp and BW Stanton during summer 2014, but at none of the other main-stem or tributary sites. The 13 °C critical temperature criterion for salmonid spawning was exceeded in early July 2014 at BW Ketchum and BW Hailey near the end of the rainbow trout critical spawning and rearing period. Temperature exceedances were most frequent at BW Stanton, where exceedances for rainbow trout and brown trout occurred from May through early July 2014 during most of the critical spawning and rearing period. Water quality and habitat availability did not seem to be limiting for fish or other aquatic organisms at most sites in the Big Wood River watershed. Water quality assessments in September 2014 determined no exceedances of total maximum daily load target levels. The availability and quality of habitat was limited at BW Stanton, where shallow-water habitat conditions prevailed. Macroinvertebrate community diversity was high at all sites except for BW Stanton, where low community diversity was attributed to low species richness and high abundances of a few tolerant taxa. Presence of low species diversity and high macroinvertebrate tolerance values at BW Stanton indicates that benthic community condition and stream health were reduced at that location. Fish surveys done in September 2014 did not indicate any significant reductions in native fish communities in the Big Wood River or its tributaries. Native rainbow trout ( Oncorhynchus mykiss ) and Wood River sculpin ( Cottus leiopomus ) were the dominant fish species in the drainage and were found at all tributary and main-stem sites. Non-native brown ( Salmo trutta ) and brook trout ( Salvelinus fontinalis ) were limited to lower drainage sites on the Big Wood River (BW Hailey and BW Stanton), and occurred in relatively low numbers.

Idaho↗

Comparison of two regression-based approaches for determining nutrient and sediment fluxes and trends in the Chesapeake Bay watershed

Nutrient and sediment fluxes and changes in fluxes over time are key indicators that water resource managers can use to assess the progress being made in improving the structure and function of the Chesapeake Bay ecosystem. The U.S. Geological Survey collects annual nutrient (nitrogen and phosphorus) and sediment flux data and computes trends that describe the extent to which water-quality conditions are changing within the major Chesapeake Bay tributaries. Two regression-based approaches were compared for estimating annual nutrient and sediment fluxes and for characterizing how these annual fluxes are changing over time. The two regression models compared are the traditionally used ESTIMATOR and the newly developed Weighted Regression on Time, Discharge, and Season (WRTDS). The model comparison focused on answering three questions: (1) What are the differences between the functional form and construction of each model? (2) Which model produces estimates of flux with the greatest accuracy and least amount of bias? (3) How different would the historical estimates of annual flux be if WRTDS had been used instead of ESTIMATOR? One additional point of comparison between the two models is how each model determines trends in annual flux once the year-to-year variations in discharge have been determined. All comparisons were made using total nitrogen, nitrate, total phosphorus, orthophosphorus, and suspended-sediment concentration data collected at the nine U.S. Geological Survey River Input Monitoring stations located on the Susquehanna, Potomac, James, Rappahannock, Appomattox, Pamunkey, Mattaponi, Patuxent, and Choptank Rivers in the Chesapeake Bay watershed. Two model characteristics that uniquely distinguish ESTIMATOR and WRTDS are the fundamental model form and the determination of model coefficients. ESTIMATOR and WRTDS both predict water-quality constituent concentration by developing a linear relation between the natural logarithm of observed constituent concentration and three explanatory variables—the natural log of discharge, time, and season. ESTIMATOR uses two additional explanatory variables—the square of the log of discharge and time-squared. Both models determine coefficients for variables for a series of estimation windows. ESTIMATOR establishes variable coefficients for a series of 9-year moving windows; all observed constituent concentration data within the 9-year window are used to establish each coefficient. Conversely, WRTDS establishes variable coefficients for each combination of discharge and time using only observed concentration data that are similar in time, season, and discharge to the day being estimated. As a result of these distinguishing characteristics, ESTIMATOR reproduces concentration-discharge relations that are closely approximated by a quadratic or linear function with respect to both the log of discharge and time. Conversely, the linear model form of WRTDS coupled with extensive model windowing for each combination of discharge and time allows WRTDS to reproduce observed concentration-discharge relations that are more sinuous in form. Another distinction between ESTIMATOR and WRTDS is the reporting of uncertainty associated with the model estimates of flux and trend. ESTIMATOR quantifies the standard error of prediction associated with the determination of flux and trends. The standard error of prediction enables the determination of the 95-percent confidence intervals for flux and trend as well as the ability to test whether the reported trend is significantly different from zero (where zero equals no trend). Conversely, WRTDS is unable to propagate error through the many (over 5,000) models for unique combinations of flow and time to determine a total standard error. As a result, WRTDS flux estimates are not reported with confidence intervals and a level of significance is not determined for flow-normalized fluxes. The differences between ESTIMATOR and WRTDS, with regard to model form and determination of model coefficients, have an influence on the determination of nutrient and sediment fluxes and associated changes in flux over time as a result of management activities. The comparison between the model estimates of flux and trend was made for combinations of five water-quality constituents at nine River Input Monitoring stations. The major findings with regard to nutrient and sediment fluxes are as follows: (1)WRTDS produced estimates of flux for all combinations that were more accurate, based on reduction in root mean squared error, than flux estimates from ESTIMATOR; (2) for 67 percent of the combinations, WRTDS and ESTIMATOR both produced estimates of flux that were minimally biased compared to observed fluxes(flux bias = tendency to over or underpredict flux observations); however, for 33 percent of the combinations, WRTDS produced estimates of flux that were considerably less biased (by at least 10 percent) than flux estimates from ESTIMATOR; (3) the average percent difference in annual fluxes generated by ESTIMATOR and WRTDS was less than 10 percent at 80 percent of the combinations; and (4) the greatest differences related to flux bias and annual fluxes all occurred for combinations where the pattern in observed concentration-discharge relation was sinuous (two points of inflection) rather than linear or quadratic (zero or one point of inflection). The major findings with regard to trends are as follows: (1) both models produce water-quality trends that have factored in the year-to-year variations in flow; (2) trends in water-quality condition are represented by ESTIMATOR as a trend in flow-adjusted concentration and by WRTDS as a flow normalized flux; (3) for 67 percent of the combinations with trend estimates, the WRTDS trends in flow-normalized flux are in the same direction and magnitude to the ESTIMATOR trends in flow-adjusted concentration, and at the remaining 33 percent the differences in trend magnitude and direction are related to fundamental differences between concentration and flux; and (4) the majority (85 percent) of the total nitrogen, nitrate, and orthophosphorus combinations exhibited long-term (1985 to 2010) trends in WRTDS flow-normalized flux that indicate improvement or reduction in associated flux and the majority (83 percent) of the total phosphorus (from 1985 to 2010) and suspended sediment (from 2001 to 2010) combinations exhibited trends in WRTDS flow-normalized flux that indicate degradation or increases in the flux delivered.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

Simulation of groundwater flow in the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and public supply uses in the Pecos County region of western Texas. The U.S. Geological Survey completed a comprehensive, integrated analysis of available hydrogeologic data to develop a numerical groundwater-flow model of the Edwards-Trinity and related aquifers in the study area in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. The active model area covers about 3,400 square miles of the Pecos County region of Texas west of the Pecos River, and its boundaries were defined to include the saturated areas of the Edwards-Trinity aquifer. The model is a five-layer representation of the Pecos Valley, Edwards-Trinity, Dockum, and Rustler aquifers. The Pecos Valley aquifer is referred to as the alluvial layer, and the Edwards-Trinity aquifer is divided into layers representing the Edwards part of the Edwards-Trinity aquifer and the Trinity part of the Edwards-Trinity aquifer, respectively. The calibration period of the simulation extends from 1940 to 2010. Simulated hydraulic heads generally were in good agreement with observed values; 1,684 out of 2,860 (59 percent) of the simulated values were within 25 feet of the observed value. The average root mean square error value of hydraulic head for the Edwards-Trinity aquifer was 34.2 feet, which was approximately 4 percent of the average total observed change in groundwater-level altitude (groundwater level). Simulated spring flow representing Comanche Springs exhibits a pattern similar to observed spring flow. Independent geochemical modeling corroborates results of simulated groundwater flow that indicates groundwater in the Edwards-Trinity aquifer in the Leon-Belding and Fort Stockton areas is a mixture of recharge from the Barilla and Davis Mountains and groundwater that has upwelled from the Rustler aquifer. The model was used to simulate groundwater-level altitudes resulting from prolonged pumping to evaluate sustainability of current and projected water-use demands. Each of three scenarios utilized a continuation of the calibrated model. Scenario 1 extended recent (2008) irrigation and nonirrigation pumping values for a 30-year period from 2010 to 2040. Projected groundwater-level changes in and around the Fort Stockton area under scenario 1 change little from current conditions, indicating that the groundwater system is near equilibrium with respect to recent (2008) pumping stress. Projected groundwater-level declines in the eastern part of the model area ranging from 5.0 to 15.0 feet are likely the result of nonequilibrium conditions associated with recent increases in pumping after a prolonged water-level recovery period of little or no pumping. Projected groundwater-level declines (from 15.0 to 31.0 feet) occurred in localized areas by the end of scenario 1 in the Leon-Belding area. Scenario 2 evaluated the effects of extended recent (2008) pumping rates as assigned in scenario 1 with year-round maximum permitted pumping rates in the Belding area. Results of scenario 2 are similar in water-level decline and extent as those of scenario 1. The extent of the projected groundwater-level decline in the range from 5.0 to 15.0 feet in the Leon-Belding irrigation area expanded slightly (about a 2-percent increase) from that of scenario 1. Maximum projected groundwater-level declines in the Leon-Belding irrigation area were approximately 31.3 feet in small isolated areas. Scenario 3 evaluated the effects of periodic increases in pumping rates over the 30-year extended period. Results of scenario 3 are similar to those of scenario 2 in terms of the areas of groundwater-level decline; however, the maximum projected groundwater-level decline increased to approximately 34.5 feet in the Leon-Belding area, and the extent of the decline was larger in area (about a 17-percent increase) than that of scenario 2. Additionally, the area of projected groundwater-level declines in the eastern part of the model area increased from that of scenario 2&mdash;two individual areas of decline coalesced into one larger area. The localized nature of the projected groundwater-level declines is a reflection of the high degree of fractured control on storage and hydraulic conductivity in the Edwards-Trinity aquifer. Additionally, the finding that simulated spring flow is highly dependent on the transient nature of hydraulic heads in the underlying aquifer indicates the importance of adequately understanding and characterizing the entire groundwater system.

Texas↗

Monitoring the water-quality response of agricultural conservation practices in the Bucks Branch watershed, Sussex County, Delaware, 2014–16

The purpose of this study was to evaluate the effects of irrigation and cover crops as conservation practices on water quality in groundwater and streams. Bucks Branch, a stream in the Nanticoke River watershed in southwestern Delaware, was identified as having one of the highest concentrations of nitrate in all surface-water sites sampled by the Delaware Department of Natural Resources and Environmental Control (DNREC). The study site is on two adjacent fields bordering Bucks Branch, one that has used irrigation since 2000 and one with dryland farming; both under conservation tillage and long-term rotation of corn, soybean, and small grain crops. A streamgage was installed near the study site fields to measure streamflow and water quality. The study area is typical of farming practices and environmental conditions throughout much of the intensively farmed agricultural land of the Coastal Plain of Delaware and surrounding parts of Maryland. Monitoring was conducted from January 2014 through June 2016. Corn was grown on both fields during the two growing seasons of the study period, and cover crops were planted before or shortly after harvest on both fields. During the second year of data collection, the effects of radish and rye grass cover crops on nutrient transport were studied. The combined results from data collected for this study show that water and nitrate moved below the root zone year round when soil moisture was high, especially after significant rainfall and frequently after irrigation. Soil water sampled 2 to 3 weeks after nutrients were applied had nitrate concentrations greater than 50 milligrams per liter as nitrogen (mg/L as N) and may be a significant source of nitrate to groundwater. Whereas recharge containing elevated nitrate concentrations also occurred under the dryland field, it was less frequent and of lower concentration than recharge under the irrigated field. Nitrate was present in all groundwater samples from these sites. Groundwater estimated to have recharged within 10 years or less had higher median concentrations of nitrate than in older water samples. The oldest groundwater encountered was over 30 years old, and had traveled along the longest, deep flowpaths from upland fields to the stream. The median nitrate concentration was 18 mg/L as N in younger water (less than 10 years old) beneath the irrigated field, compared to about 10 mg/L as N in younger water beneath the dryland field. Samples from the shallow upland wells in both study fields showed little, if any, evidence of denitrification. Several samples from deeper wells and from wells near forested riparian zone wetlands that border both fields did show partial denitrification. A mixing model estimated that between 12 and 22 percent of the nitrate discharging to the stream was lost through uptake and denitrification upstream of the streamgage on Bucks Branch. Continuous data collected at this site and evidence of denitrification in the surface-water samples showed a greater potential for loss of nitrate during the warmer months than the colder months. This pattern was similar to that seen below the streamgage at the most downstream site in the watershed. A mixed cover crop of radishes and rye was planted prior to removal (radishes) and just after harvest (rye) of the corn crop on the irrigated field. Rye grass was planted shortly after crop harvest on the dryland field. Cover crop biomass samples collected while radishes were growing and after they were killed by freezing temperatures indicates that the early planted radish crop effectively scavenged available nitrogen from the soil. Whereas radish biomass initially held more nitrogen than rye, at 55 to 8 pounds per acre, respectively, leaching of inorganic nitrate following radish die-off was minimal. Soil-water nitrate concentrations during the cover-crop growing period were lower than during the growing season prior to planting of the cover crop. There also was an increase in soil fertility and dissolved organic nitrate in samples of soil water that was likely related to increased soil microbial metabolism. Results indicate that cover crops stored plant nutrients over the winter and did not increase shallow groundwater concentrations of nitrate. Although conservation practices such as cover crops and nutrient management have been applied to these fields, there was still significant leaching of nitrate to groundwater, especially under the irrigated field. This will likely continue to be a challenge in this area and other parts of the Coastal Plain where soil moisture capacity is relatively low and managing irrigation around rainfall is difficult. Cover crops, when planted in standing corn, are one practice that can effectively pull nitrate from below the root zone to the top layer of soil, thus limiting the amount of potential nitrate leaching to groundwater. Irrigation management that would lower average soil moisture conditions during the growing season also could potentially limit nitrogen transport.

Delaware↗