USGS Science⌕ Search

SEARCH · USGS Science

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,045 records · Page 58Linked to original sources

Water-quality trends for selected sites and constituents in the international Red River of the North Basin, Minnesota and North Dakota, United States, and Manitoba, Canada, 1970–2017

A comprehensive study to evaluate water-quality trends, while considering natural hydroclimatic variability, in the Red River of the North Basin and assess water-quality conditions for the Red River of the North crossing the international boundary near Emerson, Manitoba, Canada (the binational site), was completed by the U.S. Geological Survey in cooperation with the International Joint Commission, North Dakota Department of Environmental Quality, and Minnesota Pollution Control Agency and in collaboration with Manitoba Sustainable Development and Environment and Climate Change Canada. The international Red River of the North Basin encompasses 3 U.S. States (South Dakota, North Dakota, and Minnesota) and 1 Canadian Province (Manitoba). Water quality in the Red River of the North Basin is of concern for both Federal governments and State and Provincial governments. Water-quality objectives have been previously established for selected dissolved ions and recently (2019) proposed for selected nutrients for the binational site. In the current (2020) study, water-quality data from State, Provincial, and Federal agencies in the United States and Canada for sites in the Red River of the North Basin from 1970 to 2017 were compiled and used for trend analysis. Trend analysis using a water-quality dataset from multiple agencies that collect water-quality data for various objectives presented multiple challenges. The trend-analysis approach was able to accommodate differences in water-quality data caused by field-collection and laboratory-analytical method differences, disparities in sampling frequencies, and spatial and temporal gaps in data. Most of these challenges were overcome by the statistical tool, R–QWTREND, which identifies trends in concentration unrelated to variability in streamflow. The integrated basin approach used in the current study, combined with comparing current data trends with historical trends, provided valuable insights into understanding how water quality is changing spatially (34 sites analyzed for a recent period, 2000–15) and temporally (5 sites analyzed for a 45-year historical period, 1970–2015) within the Red River of the North Basin. One of the most consistent spatial and temporal changes observed in the current study was increasing concentrations of sulfate among tributary and main-stem sites since 2000. For some sites, increases were detected starting as early as 1985. Total dissolved solids and chloride concentrations had spatial and temporal patterns like sulfate. Although R–QWTREND removes the variability in constituent concentration caused by natural streamflow variability, all variability in sulfate caused by hydroclimatic variability may not be captured because of changes in hydrologic pathways and changes in the contributions of sulfate from various natural sources. Nutrient concentrations demonstrated less consistent spatial and temporal changes than sulfate, and changes in nutrient concentrations were assumed to be more closely tied to human-induced rather than natural changes. Nitrate-plus-nitrite concentrations were mostly increasing in the upper Red River of the North subbasin, and for nitrate plus nitrite and total nitrogen, the Sheyenne River subbasin had consistent decreasing concentrations. Since 2000, total phosphorus has decreased in the upper Red River of the North subbasin, but total phosphorus concentration has increased for sites in the lower Red River of the North subbasin, and for some main-stem sites, concentrations have been increasing since 1985. Unlike sulfate, the pattern in historical trends for total phosphorus for the main-stem sites differed from tributary sites, indicating that human-induced changes affected tributaries and main-stem sites differently. The more detailed evaluation of flow-averaged water-quality conditions for the binational site provided an understanding of how loads have changed over time and what proportion of the year and season concentrations are expected to exceed water-quality objectives. In a basin with highly variable streamflow like the Red River of the North, the trend in flow-averaged load (assuming streamflow conditions are the same year after year) provided a robust measure of change over time. Increasing concentrations of sulfate, chloride, total dissolved solids, and total phosphorus since 1985 for the binational site resulted in longer periods of exceedance of water-quality objectives per year occurring over time. For total nitrogen, decreasing concentrations resulted in shorter periods of exceedance per year during 1980 to 2015, but concentrations were still expected to exceed the water-quality objective about half the year. Periods of when exceedances were likely to occur during the year were affected by the source and transport mechanisms of the constituent. Trend results from this effort identified how water quality has changed across the basin, and further investigation would help to identify causes for the trends observed here. Information from the current study provides a basis for future trend attribution studies, evaluation of water-quality objectives, and development of comprehensive strategies for reducing nutrients to desired targets and establishes a baseline for tracking future progress in the Red River of the North Basin.

Minnesota, North Dakota, South Dakota, Manitoba↗

Simulated Effects of Year 2030 Water-Use and Land-Use Changes on Streamflow near the Interstate-495 Corridor, Assabet and Upper Charles River Basins, Eastern Massachusetts

Continued population growth and land development for commercial, industrial, and residential uses have created concerns regarding the future supply of potable water and the quantity of ground water discharging to streams in the area of Interstate 495 in eastern Massachusetts. Two ground-water models developed in 2002-2004 for the Assabet and Upper Charles River Basins were used to simulate water supply and land-use scenarios relevant for the entire Interstate-495 corridor. Future population growth, water demands, and commercial and residential growth were projected for year 2030 by the Metropolitan Area Planning Council. To assess the effects of future development on subbasin streamflows, seven scenarios were simulated by using existing computer-based ground-water-flow models with the data projected for year 2030. The scenarios incorporate three categories of projected 2030 water- and land-use data: (1) 2030 water use, (2) 2030 land use, and (3) a combination of 2030 water use and 2030 land use. Hydrologic, land-use, and water-use data from 1997 through 2001 for the Assabet River Basin study and 1989 through 1998 for the Upper Charles River Basin study were used to represent current conditions - referred to as 'basecase' conditions - in each basin to which each 2030 scenario was compared. The effects of projected 2030 land- and water-use change on streamflows in the Assabet River Basin depended upon the time of year, the hydrologic position of the subbasin in the larger basin, and the relative areas of new commercial and residential development projected for a subbasin. Effects of water use and land use on streamflow were evaluated by comparing average monthly nonstorm streamflow (base flow) for March and September simulated by using the models. The greatest decreases in streamflow (up to 76 percent in one subbasin), compared to the basecase, occurred in September, when streamflows are naturally at their lowest level. By contrast, simulated March streamflows decreased less than 6.5 percent from basecase streamflows in all subbasins for all scenarios. The simulations showed similar effects in the Upper Charles River Basin, but increased water use contributed to decreased simulated streamflow in most subbasins. Simulated changes in March streamflows for 2030 in the Upper Charles River Basin were within +- 6 percent of the basecase for all scenarios and subbasins. Percentage decreases in simulated September streamflows for 2030 were greater than in March but less than the September decreases that resulted for some subbasins in the Assabet River Basin. Only two subbasins of the Upper Charles River Basin had projected decreases greater than 5 percent. In the Mill River subbasin, the decrease was 11 percent, and in the Mine Brook subbasin, 6.6 percent. Changes in water use and wastewater return flow generally were found to have the greatest effect in the summer months when streamflow and aquifer recharge rates are low and water use is high. September increases in main-stem streamflow of both basins were due mainly to increased discharge of treated effluent from wastewater-treatment facilities on the main-stem rivers. In the Assabet River Basin, wastewater-treatment-facility discharge became a smaller proportion of total streamflow with distance downstream. In contrast, wastewater-treatment facility discharge in the Upper Charles River Basin became a greater proportion of streamflow with distance downstream. The effects of sewer-line extension and low-impact development on streamflows in two different subbasins of the Assabet River Basin also were simulated. The result of extending sewer lines with a corresponding decrease in septic-system return flow caused September streamflows to decrease as much as 15 percent in the Fort Pond Brook subbasin. The effect of low-impact development was simulated in the Hop Brook subbasin in areas projected for commercial development. In this simulation, the greater the area where low-i

Scientific Investigations Report↗

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming↗

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗

Factors Affecting Spatial and Temporal Variability in Nutrient and Pesticide Concentrations in the Surficial Aquifer on the Delmarva Peninsula

Water quality in the unconfined, unconsolidated surficial aquifer on the Delmarva Peninsula is influenced by the availability of soluble ions from natural and human sources, and by geochemical factors that affect the mobility and fate of these ions within the aquifer. Ground-water samples were collected from 60 wells completed in the surficial aquifer of the peninsula in 2001 and analyzed for major ions, nutrients, and selected pesticides and degradation products. Analytical results were compared to similar data from a subset of sampled wells in 1988, as well as to land use, soils, geology, depth, and other potential explanatory variables to demonstrate the effects of natural and human factors on water quality in the unconfined surficial aquifer. This study was conducted as part of the National Water-Quality Assessment Program of the U.S. Geological Survey, which is designed (in part) to describe the status and trends in ground-water quality and to provide an understanding of natural and human factors that affect ground-water chemistry in different parts of the United States. Results of this study may be useful for water-resources managers tasked with addressing water-quality issues of local and regional importance because the surficial aquifer on the Delmarva Peninsula is a major source of water for domestic and public supply and provides the majority of flow in local streams. Human impacts are apparent in ground-water quality throughout the surficial aquifer. The surficial aquifer on the Delmarva Peninsula is generally sandy and very permeable with well-oxygenated ground water. Dissolved constituents found throughout various depths of the unconfined aquifer are likely derived from the predominantly agricultural practices on the peninsula, although effects of road salt, mineral dissolution, and other natural and human influences are also apparent in some areas. Nitrate occurred at concentrations exceeding natural levels in many areas, and commonly exceeded 10 milligrams per liter (as nitrogen). In addition to land use in the aquifer recharge area, concentrations of nitrate in ground water are related to regional patterns in soil drainage that affect underlying aquifer redox conditions. Over the peninsula, nitrate concentrations are not related to recharge date of the water, but are positively correlated with depth in shallow wells screened beneath agricultural areas. Nitrate concentrations increased in oxic areas (dissolved oxygen greater than 1 milligram per liter) of the deeper part of the surficial aquifer used for domestic supply by an average of about 2 milligrams per liter between 1988 and 2001, although no changes were apparent in shallower parts of the aquifer over that same period. Water in the surficial aquifer generally flows from land-surface recharge to surface-water discharge areas in less than 30 years. As a result, the entire flow system in the surficial aquifer has likely been affected by human activities on and near the land surface over the past several decades. Pesticide compounds occurred widely at low levels throughout the surficial aquifer. The most commonly used herbicides (metolachlor, alachlor, and atrazine) were the most commonly detected. These pesticides primarily occurred in ground water in the form of degradation products. The widespread occurrence of these and other pesticide compounds reflects their abundant use as well as chemical properties and aquifer characteristics that allow their movement into ground water. Mixtures of pesticides are common. Most samples contained at least 3 different compounds; several samples contained as many as 11. Pesticide concentrations in the surficial aquifer are relatively high beneath recharge areas with well-drained soils in the shallow part of the aquifer and in oxic environments throughout the surficial aquifer. Concentrations are generally below existing drinking-water standards, although standards are not available for all of the pesticide compound

Scientific Investigations Report↗

Ecological effects of pinyon-juniper removal in the Western United States—A synthesis of scientific research, January 2014–March 2021

Executive Summary Increasing density of pinyon ( Pinus spp. ) and juniper ( Juniperus spp. ) woodlands (hereinafter “pinyon-juniper”), as well as expansion of these woodlands into adjacent shrublands and grasslands, has altered ecosystem function and wildlife habitat across large areas of the interior western United States. Although there are many natural and human-caused drivers of woodland infilling and expansion, restoration of sagebrush ( Artemisia spp. ) habitat through removal of pinyon-juniper is considered an urgent management objective in many locations, particularly in support of sagebrush-dependent wildlife species of conservation concern. In December 2020, the Bureau of Land Management (BLM) established the Pinyon-Juniper Management Categorical Exclusion (PJCX) to expedite the regulatory process for pinyon-juniper removal projects on public lands, largely intended to benefit mule deer ( Odocoileus hemionus ) and greater sage-grouse ( Centrocercus urophasianus ) habitats. During final preparation of this report, the BLM discontinued use of the PJCX (as of November, 2022), but the pinyon-juniper tree removal techniques assessed in this report are commonly used and understanding their effects remains relevant to land use planning. To address areas of uncertainty relative to potential ecological effects of the PJCX, we conducted a review of the peer-reviewed science literature to better understand the likely responses of vegetation, environmental (for example, soils), and wildlife variables to specific tree removal techniques permitted by the PJCX. In brief, the PJCX permitted removal of trees by either manual cutting, mechanical cutting, or mastication; allowed certain methods to redistribute or remove resulting tree biomass after treatment; and prohibited broadcast burning, roadbuilding, removal of old-growth, and seeding of non-native species. Specifically, we conducted our review to address the following questions: How will PJCX removal techniques affect plant communities, soils, and abiotic resources? How do these pinyon-juniper removal techniques affect wildlife communities, including both woodland- and sagebrush-dependent species? What are the potential ecological implications of different pinyon-juniper removal treatment types and implementation strategies (for example, treatment sizes) over time? What are the most important gaps in our scientific understanding of how treatments might affect targeted ecosystems over space and time (for example, potential effects of climate change)? To answer these questions, we considered studies related to pinyon-juniper ecosystems, focusing on research that occurred over a large portion of the interior western United States that is the primary focus of the PJCX. We also focused on papers published from 2014 onward, to avoid excessive overlap with other recent reviews on pinyon-juniper management effects. Using strict criteria, including only considering research that tested responses for statistical significance, we identified 48 papers that primarily examined treatment effects on vegetation and other environmental variables (1,709 responses), and 11 papers that addressed effects on wildlife (132 responses). Responses to the PJCX-permitted treatments were summarized as either positive (that is, a significant increase), negative (that is, a significant decrease), or non-significant (that is, no significant difference). Responses were assigned to categories (for example, Native Annual Grass/Forb Abundance) and hierarchical treatment levels. We found that there were large proportions of non-significant responses among all categories combined, with roughly half or more of all responses non-significant (48 percent for wildlife, 60 percent for vegetation-environmental), comparable to other recent systematic reviews of pinyon-juniper treatment effects. However, we also found that when there were significant responses, some important trends potentially emerged. Important undesirable outcomes included far more positive than negative responses of exotic grass and forb abundance among nearly all treatment types. Cutting treatments were also more likely to decrease biocrust cover and microbial activity. Potentially beneficial outcomes included mostly positive responses among sagebrush obligate species, including more positive than negative responses for mule deer and sage-grouse. Some treatment types (for example, mastication) also resulted in more positive than negative responses for native grasses and forbs (although, non-significant responses were the majority). We also highlighted many limitations of this review, including how responses often come from few studies, and how some response-treatment category combinations lack adequate response data. Moreover, the existing research is often insufficient to address many key questions about treatment effects, largely owing to short time-scales and limited spatial extents of observations, which do not match the size of treatments being implemented by land managers, nor capture long-term, post-treatment ecological dynamics. We also identify a lack of research that addresses key interactions that could undermine restoration objectives, including potential effects of climate change and grazing on post-treatment environments. Thus, we emphasize the importance of integrating these factors into future pinyon-juniper treatment research, and we stress the need for use of monitoring programs and research studies that partake in data collection and analysis over long durations and broad spatial scales.

Open-File Report↗

Effects of recreational flow releases on natural resources of the Indian and Hudson Rivers in the Central Adirondack Mountains, New York, 2004-06

The U.S. Geological Survey (USGS), the New York State Department of Environmental Conservation (NYSDEC), and Cornell University carried out a cooperative 2-year study from the fall of 2004 through the fall of 2006 to characterize the potential effects of recreational-flow releases from Lake Abanakee on natural resources in the Indian and Hudson Rivers. Researchers gathered baseline information on hydrology, temperature, habitat, nearshore wetlands, and macroinvertebrate and fish communities and assessed the behavior and thermoregulation of stocked brown trout in study reaches from both rivers and from a control river. The effects of recreational-flow releases (releases) were assessed by comparing data from affected reaches with data from the same reaches during nonrelease days, control reaches in a nearby run-of-the-river system (the Cedar River), and one reach in the Hudson River upstream from the confluence with the Indian River. A streamgage downstream from Lake Abanakee transmitted data by satellite from November 2004 to November 2006; these data were used as the basis for developing a rating curve that was used to estimate discharges for the study period. River habitat at most study reaches was delineated by using Global Positioning System and ArcMap software on a handheld computer, and wetlands were mapped by ground-based measurements of length, width, and areal density. River temperature in the Indian and Hudson Rivers was monitored continuously at eight sites during June through September of 2005 and 2006; temperature was mapped in 2005 by remote imaging made possible through collaboration with the Rochester Institute of Technology. Fish communities at all study reaches were surveyed and characterized through quantitative, nearshore electrofishing surveys. Macroinvertebrate communities in all study reaches were sampled using the traveling-kick method and characterized using standard indices. Radio telemetry was used to track the movement and persistence of stocked brown trout (implanted with temperature-sensitive transmitters) in the Indian and Hudson Rivers during the summer of 2005 and in all three rivers during the summer of 2006. The releases had little effect on river temperatures, but increased discharges by about one order of magnitude. Regardless of the releases, river temperatures at all study sites commonly exceeded the threshold known to be stressful to brown trout. At most sites, mean and median water temperatures on release days were not significantly different, or slightly lower, than water temperatures on nonrelease days. Most differences were very small and, thus, were probably not biologically meaningful. The releases generally increased the total surface area of fast-water habitat (rapids, runs, and riffles) and decreased the total surface area of slow-water habitat (pools, glides, backwater areas, and side channels). The total surface areas of wetlands bordering the Indian River were substantially smaller than the surface areas bordering the Cedar River; however, no channel geomorphology or watershed soil and topographic data were assessed to determine whether the releases or other factors were mainly responsible for observed differences. Results from surveys of resident biota indicate that the releases generally had a limited effect on fish and macroinvertebrate communities in the Indian River and had no effect on communities in the Hudson River. Compared to fish data from Cedar River control sites, the impoundment appeared to reduce total density, biomass, and richness in the Indian River at the first site downstream from Lake Abanakee, moderately reduce the indexes at the other two sites on the Indian River, and slightly reduce the indexes at the first Hudson River site downstream from the confluence with the Indian River. The densities of individual fish populations at all Indian River sites were also reduced, but related effects on fish populations in the Hudson River were less evident. Altho

Scientific Investigations Report↗

Median and low-flow characteristics for streams under natural and diverted conditions, northeast Maui, Hawaii

Flow-duration statistics under natural (undiverted) and diverted flow conditions were estimated for gaged and ungaged sites on 21 streams in northeast Maui, Hawaii. The estimates were made using the optimal combination of continuous-record gaging-station data, low-flow measurements, and values determined from regression equations developed as part of this study. Estimated 50- and 95-percent flow duration statistics for streams are presented and the analyses done to develop and evaluate the methods used in estimating the statistics are described. Estimated streamflow statistics are presented for sites where various amounts of streamflow data are available as well as for locations where no data are available. Daily mean flows were used to determine flow-duration statistics for continuous-record stream-gaging stations in the study area following U.S. Geological Survey established standard methods. Duration discharges of 50- and 95-percent were determined from total flow and base flow for each continuous-record station. The index-station method was used to adjust all of the streamflow records to a common, long-term period. The gaging station on West Wailuaiki Stream (16518000) was chosen as the index station because of its record length (1914-2003) and favorable geographic location. Adjustments based on the index-station method resulted in decreases to the 50-percent duration total flow, 50-percent duration base flow, 95-percent duration total flow, and 95-percent duration base flow computed on the basis of short-term records that averaged 7, 3, 4, and 1 percent, respectively. For the drainage basin of each continuous-record gaged site and selected ungaged sites, morphometric, geologic, soil, and rainfall characteristics were quantified using Geographic Information System techniques. Regression equations relating the non-diverted streamflow statistics to basin characteristics of the gaged basins were developed using ordinary-least-squares regression analyses. Rainfall rate, maximum basin elevation, and the elongation ratio of the basin were the basin characteristics used in the final regression equations for 50-percent duration total flow and base flow. Rainfall rate and maximum basin elevation were used in the final regression equations for the 95-percent duration total flow and base flow. The relative errors between observed and estimated flows ranged from 10 to 20 percent for the 50-percent duration total flow and base flow, and from 29 to 56 percent for the 95-percent duration total flow and base flow. The regression equations developed for this study were used to determine the 50-percent duration total flow, 50-percent duration base flow, 95-percent duration total flow, and 95-percent duration base flow at selected ungaged diverted and undiverted sites. Estimated streamflow, prediction intervals, and standard errors were determined for 48 ungaged sites in the study area and for three gaged sites west of the study area. Relative errors were determined for sites where measured values of 95-percent duration discharge of total flow were available. East of Keanae Valley, the 95-percent duration discharge equation generally underestimated flow, and within and west of Keanae Valley, the equation generally overestimated flow. Reduction in 50- and 95-percent flow-duration values in stream reaches affected by diversions throughout the study area average 58 to 60 percent.

Hawaii↗

Application of a watershed model (HSPF) for evaluating sources and transport of pathogen indicators in the Chino Basin drainage area, San Bernardino County, California

A watershed model using Hydrologic Simulation Program-FORTRAN (HSPF) was developed for the urbanized Chino Basin in southern California to simulate the transport of pathogen indicator bacteria, evaluate the flow-component and land-use contributions to bacteria contamination and water-quality degradation throughout the basin, and develop a better understanding of the potential effects of climate and land-use change on water quality. The calibration of the model for indicator bacteria was supported by historical data collected before this study and by samples collected by the U.S. Geological Survey from targeted land-use areas during storms in water-year 2004. The model was successfully calibrated for streamflow at 5 gage locations representing the Chino Creek and Mill Creek drainages. Although representing pathogens as dissolved constituents limits the model's ability to simulate the transport of pathogen indicator bacteria, the bacteria concentrations measured over the period 1998-2004 were well represented by the simulated concentrations for most locations. Hourly concentrations were more difficult to predict because of high variability in measured bacteria concentrations. In general, model simulations indicated that the residential and commercial land uses were the dominant sources for most of the pathogen indicator bacteria during low streamflows. However, simulations indicated that land used for intensive livestock (dairies and feedlots) and mixed agriculture contributed the most bacteria during storms. The calibrated model was used to evaluate how various land use, air temperature, and precipitation scenarios would affect flow and transport of bacteria. Results indicated that snow pack formation and melt were sensitive to changes in air temperature in the northern, mountainous part of the Chino Basin, causing the timing and magnitude of streamflow to shift in the natural drainages and impact the urbanized areas of the central Chino Basin. The relation between bacteria concentrations and air temperature was more complicated, and did not substantially affect the quality of water discharging from the Chino Basin into the Santa Ana River. Changes in precipitation had a greater basin-wide affect on bacteria concentrations than changes in air temperature, and varied according to location. Drainages representing natural conditions had a decrease in bacteria concentrations in correlation with an increase in precipitation, whereas drainages in the central and southern part of the Chino Basin had an increase in bacteria concentrations. Drier climate conditions tended to result in higher sensitivity of simulated bacteria concentrations to changes in precipitation. Simulated bacteria concentrations in wetter climates were usually less sensitive to changes in precipitation because bacteria transport becomes more dependent on the land-use specified bacteria loading rates and the storage limits. Bacteria contamination from impervious-area runoff is affected to a greater degree by drier climates, whereas contamination from pervious-area runoff is affected to a greater degree by wetter climates. Model results indicated that the relation between precipitation, runoff, and bacteria contamination is complicated because after the initial bacteria washoff and transport from the land surfaces during the beginning of a storm period, subsequent runoff has fewer bacteria available for washoff, which then dilutes the concentrations of bacteria in the downstream reach. It was illustrated that pathogen indicator bacteria transport depends most significantly on the relation of imperviousness to runoff, which controls the frequency, and often the magnitude, of transport, and on the contribution of higher bacteria loading rates used for pervious land areas, especially intensive feedlots, to the infrequent, but very high, peaks of bacteria contamination. The indicator bacteria transport model for the Chino Basin was based on the assumption that no

California↗

Development of regional curves relating bankfull-channel geometry and discharge to drainage area for streams in Pennsylvania and selected areas of Maryland

Natural-stream designs are commonly based on the dimensions of the bankfull channel, which is capable of conveying discharges that transport sediment without excessive erosion or deposition. Regional curves relate bankfull-channel geometry and discharge to drainage area in watersheds with similar runoff characteristics and commonly are utilized by practitioners of natural-stream design to confirm or refute selection of the field-identified bankfull channel. Data collected from 66 streamflow-gaging stations and associated stream reaches between December 1999 and December 2003 were used in one-variable ordinary least-squares regression analyses to develop regional curves relating drainage area to cross-sectional area, discharge, width, and mean depth of the bankfull channel. Watersheds draining to these stations are predominantly within the Piedmont, Ridge and Valley, and Appalachian Plateaus Physiographic Provinces of Pennsylvania and northern Maryland. Statistical analyses of physiography, percentage of watershed area underlain by carbonate bedrock, and percentage of watershed area that is glaciated indicate that carbonate bedrock, not physiography or glaciation, has a controlling influence on the slope of regional curves. Regional curves developed from stations in watersheds underlain by 30 percent or less carbonate bedrock generally had steeper slopes than the corresponding relations developed from watersheds underlain by greater than 30 percent carbonate bedrock. In contrast, there is little evidence to suggest that regional curves developed from stations in the Piedmont or Ridge and Valley Physiographic Province are different from the corresponding relations developed from stations in the Appalachian Plateaus Physiographic Province. On the basis of these findings, regional curves are presented to represent two settings that are independent of physiography: (1) noncarbonate settings characterized by watersheds with carbonate bedrock underlying 30 percent or less of watershed area, and (2) carbonate settings characterized by watersheds with carbonate bedrock underlying greater than 30 percent of watershed area. All regional curves presented in this report have slopes that are significantly different from zero and normally distributed residuals that vary randomly with drainage area. Drainage area explains the most variability in bankfull cross-sectional area and bankfull discharge in the noncarbonate setting (R2 = 0.92 for both). Less variability is explained in bankfull width and depth (R2 = 0.81 and 0.72, respectively). Regional curves representing the carbonate setting are generally not as statistically robust as the corresponding noncarbonate relations because there were only 11 stations available to develop these curves and drainage area cannot explain variance resulting from karst features. The carbonate regional curves generally are characterized by less confidence, lower R2 values, and higher residual standard errors. Poor representation by watersheds less than 40 mi2 causes the carbonate regional curves for bankfull discharge, cross-sectional area, and mean depth to be disproportionately influenced by the smallest watershed (values of Cook's Distance range from 3.6 to 8.4). Additional bankfull discharge and channel-geometry data from small watersheds might reduce this influence, increase confidence, and generally improve regional curves representing the carbonate setting. Limitations associated with streamflow-gaging station selection and development of the curves result in some constraints for the application of regional curves presented in this report. These curves apply only to streams within the study area in watersheds having land use, streamflow regulation, and drainage areas that are consistent with the criteria used for station selection. Regardless of the setting, the regional curves presented here are not intended for use as the sole method for estimation of bankfull characteristics; however, th

Scientific Investigations Report↗

Simulated effects of existing and proposed surface-water impoundments and gas-well pads on streamflow and suspended sediment in the Cypress Creek watershed, Arkansas

Cypress Creek is located in central Arkansas and is the main tributary to Brewer Lake, which serves as the primary water supply for Conway, Arkansas, and the surrounding areas. A model of the Cypress Creek watershed was developed and calibrated in cooperation with Southwestern Energy Company using detailed precipitation, streamflow, and discrete suspended-sediment data collected from 2009 through 2012. These data were used with a Hydrologic Simulation Program—FORTRAN model to address different potential gas-extraction activities within the watershed. The calibrated Hydrologic Simulation Program—FORTRAN model was used to simulate four land-use scenarios and examine the potential effects of these land-use changes on the streamflow and water quality within the Cypress Creek watershed. These simulated scenarios included (1) the conversion of all nonforested land to forest, representing a time period before extensive grazing activities and no gas-extraction activities; (2) a land-use change to that of 1949, representing a time period with some grazing activities and no gas-extraction activities; (3) a time period with current land-use conditions, but without any gas-extraction activities, that is, the exclusion of gas-well pads/pipelines, associated gravel roads, and surface-water impoundments; and (4) a time period with current land-use conditions, but with increased gas-extraction activities (for example, increased gas-well pad and surface-water impoundment activities) to represent a possible future natural gas full-development condition for the area. A current-conditions simulation also was built and calibrated and represents the current conditions (2013) within the watershed. This simulation was used as the comparison basis for the four land-use scenarios described above. The current-conditions simulation used the 2006 land-use conditions, which consisted primarily of forest and pasture, as well as the current (2013) 35 gas-well pads and pipelines and 6 surface-water impoundments, which account for approximately 1.6 percent of the land use. Simulating a time period before extensive-grazing activities and no gas-extraction activities for scenario 1 resulted in a decrease in suspended-sediment loads and volume of streamflow within the Cypress Creek watershed compared to the current-conditions simulation. Simulating a time period before any gas-extraction activities but with some grazing activities for scenario 2 also resulted in a decrease in suspended-sediment loads and volume of streamflow within the Cypress Creek watershed. Simulating current conditions, but without any natural gas-pad land use or related activities (including pipelines and associated gravel roads), for scenario 3 resulted in mostly unchanged suspended-sediment loads and volume of streamflow within the Cypress Creek watershed, as compared to the current-conditions simulation. Finally, simulating potential future conditions of increased gas-well pad and surface-water impoundment activities for scenario 4 resulted in a decrease (compared to the current-conditions simulation) in suspended-sediment loads and a slight increase of volume of streamflow within the Cypress Creek watershed. The Arkansas Natural Resources Commission and the Arkansas Department of Environmental Quality list suspended sediment from “poor pastures” as a primary source of nonpoint-source pollution in north-central Arkansas, but unpaved (gravel) roads are another important source of suspended sediment. Because of the high sediment-loading rates associated with gravel roads and the large amount of pasture within the watershed, the factors most responsible for suspended sediment within the Cypress Creek watershed are likely associated more with the pastureland and gravel roads, than factors associated with gas-well pads/pipelines.

Arkansas↗

Occurrence of Selected Nutrients, Trace Elements, and Organic Compounds in Streambed Sediment in the Lower Chena River Watershed near Fairbanks, Alaska, 2002-03

In 2002-03, the U.S. Geological Survey collected samples of streambed sediment at 18 sites in the lower Chena River watershed for analysis of selected nutrients, traces elements, and organic compounds. The purpose of the project was to provide Federal, State, and local agencies as well as neighborhood committees, with information for consideration in plans to improve environmental conditions in the watershed. The exploratory sampling program included analysis of streambed sediment from the Chena River and Chena Slough, a tributary to the Chena River. Results were compared to streambed-sediment guidelines for the protection of aquatic life and to 2001-02 sediment data from Noyes Slough, a side channel of the lower Chena River. The median total phosphorus concentration in Chena Slough sediment samples, 680 milligrams per kilogram (mg/kg), was two orders of magnitude greater than median total phosphorus concentration in Chena River sediment samples of 5.2 mg/kg. Median concentrations of chloride and sulfate also were greater in Chena Slough samples. Low concentrations of nitrate were detected in most of the Chena Slough samples; nitrate concentrations were below method reporting limits or not detected in Chena River sediment samples. Streambed-sediment samples were analyzed for 24 trace elements. Arsenic, nickel, and zinc were the only trace elements detected in concentrations that exceeded probable-effect levels for the protection of aquatic life. Concentrations of arsenic in Chena Slough samples ranged from 11 to 70 mg/kg and concentrations in most of the samples exceeded the probable-effect guideline for arsenic of 17 mg/kg. Arsenic concentrations in samples from the Chena River ranged from 9 to 12 mg/kg. The background level for arsenic in the lower Chena River watershed is naturally elevated because of significant concentrations of arsenic in local bedrock and ground water. Sources of elevated concentrations of zinc in one sample, and of nickel in two samples, are unknown. With the exception of elevated arsenic levels in samples from Chena Slough, the occurrence and concentration of trace elements in the streambed sediments of Chena Slough and Chena River were similar to those in Noyes Slough sediment. Sediment samples were analyzed for 78 semivolatile organic compounds and 32 organochlorine pesticides and polychlorinated biphenyls (PCBs). Low concentrations of dimethylnaphthalene and p-Cresol were detected in most Chena Slough and Chena River sediment samples. The number of semivolatile organic compounds detected ranged from 5 to 21 in most Chena Slough sediment samples. In contrast, three or fewer semivolatile organic compounds were detected in Chena River sediment samples, most likely because chemical-matrix interference resulted in elevated reporting limits for organochlorine compounds in the Chena River samples. Low concentrations of fluoranthene, pyrene, and phenanthrene were detected in Chena Slough sediment. Relatively low concentrations of DDT or its degradation products, DDD and DDE, were detected in all Chena Slough samples. Concentrations of total DDT (DDT+DDD+DDE) in two Chena Slough sediment samples exceeded the effectsrange median aquatic-life criteria of 46.1 micrograms per kilogram (ug/kg). DDT concentrations in Chena River streambed-sediment samples were less than 20 ug/kg. Low concentrations of PCB were detected in two Chena Slough streambed-sediment samples. None of the concentrations of the polychlorinated biphenyls or semivolatile organic compounds for which the samples were analyzed exceeded available guidelines for the protection of aquatic life. With the exception of elevated total DDT in two Chena Slough samples, the occurrence and concentration of organochlorine compounds in Chena Slough and Chena River sediment were similar to those in samples collected from Noyes Slough in 2001-02.

Scientific Investigations Report↗

Hydrology-driven chemical loads transported by the Green River to the Lower Duwamish Waterway near Seattle, Washington, 2013–17

The sediments in the Lower Duwamish Waterway Superfund site in Seattle, Washington, are contaminated with chemicals including metals such as arsenic, polychlorinated biphenyls (PCBs), carcinogenic polycyclic aromatic hydrocarbons (cPAHs), and dioxins/furans from decades of intense anthropogenic activities. The U.S. Geological Survey, in cooperation with the Washington State Department of Ecology, collected new data from 2013 to 2017 to estimate sediment and chemical loads transported by the Green/Duwamish River to the Lower Duwamish Waterway Superfund site (the final 8-kilometer reach of the river) in support of sediment remediation within the site. Chemical loads were calculated as the product of river suspended-sediment loads and suspended sediment-bound chemical concentrations measured at river kilometer 16.7. Using four different approaches, annual suspended sediment-bound chemical load estimates transported by the river to the Lower Duwamish Waterway were in the range of 1,120–1,470 kilograms arsenic, 2,810–8,200 grams (g) toxic equivalent cPAHs, 205–407 milligrams toxic equivalent dioxins/furans, and 340–1,180 g PCBs. Storm events contributed a disproportionately large amount of the load of anthropogenic organic compounds such as cPAHs (54 percent), dioxins/furans (44 percent), and PCBs (52 percent) as compared to overall time (17 percent). Chemical concentrations and load estimates often were underestimated using results from unfiltered water analysis only, especially in samples with high suspended-sediment concentrations and for hydrophobic organic chemicals such as cPAHs that prefer to sorb to particulates and are at low concentrations near or below the analytical limits of water methods. For metals and PCBs, the dissolved concentration was relatively low and consistent between sampling events, whereas the suspended sediment-bound chemical concentrations contributed most of the chemical concentration in the water column during periods of high river suspended-sediment concentrations. However, the dissolved fraction, on average, contributed more than one-third of the estimated total chemical load in the river system for arsenic and PCBs, even given the hydrophobic nature of the chemicals. These results suggest that the sum of the chemical concentrations measured on two separate fractions—the particulate fraction and the dissolved fraction—more fully represents the total chemical concentration as compared to analysis of an unfiltered water sample, especially in samples with high suspended-sediment concentrations. Most of the suspended-sediment load (97 percent) and sediment-bound chemical load (92–94 percent) occurred during the wet winter half of the year from October 15 to April 14. However, the highest sediment-bound chemical concentrations often occurred during short intense storms or “first flush” autumn runoff events during the dry summer half of the year from April 15 to October 14. Because of the highly variable and dynamic river system characteristics (including precipitation, discharge, sediment concentration, and tidal fluctuations), it is critical to characterize the occurrence, frequency, concentrations, and loads during extreme conditions (for example, when the river is affected by storm-derived runoff) rather than time-averaged conditions. These short extreme events have a high potential for acute effects on ecological and human health, and may have a great influence on the effectiveness of the sediment remediation activities that are underway in the Lower Duwamish Waterway.

Washington↗

Predicting lake trophic state by relating Secchi-disk transparency measurements to Landsat-satellite imagery for Michigan inland lakes, 2003-05 and 2007-08

Inland lakes are an important economic and environmental resource for Michigan. The U.S. Geological Survey and the Michigan Department of Natural Resources and Environment have been cooperatively monitoring the quality of selected lakes in Michigan through the Lake Water Quality Assessment program. Sampling for this program began in 2001; by 2010, 730 of Michigan’s 11,000 inland lakes are expected to have been sampled once. Volunteers coordinated by the Michigan Department of Natural Resources and Environment began sampling lakes in 1974 and continue to sample (in 2010) approximately 250 inland lakes each year through the Michigan Cooperative Lakes Monitoring Program. Despite these sampling efforts, it still is impossible to physically collect measurements for all Michigan inland lakes; however, Landsat-satellite imagery has been used successfully in Minnesota, Wisconsin, Michigan, and elsewhere to predict the trophic state of unsampled inland lakes greater than 20 acres by producing regression equations relating in-place Secchi-disk measurements to Landsat bands. This study tested three alternatives to methods previously used in Michigan to improve results for predicted statewide Trophic State Index (TSI) computed from Secchi-disk transparency (TSI (SDT)). The alternative methods were used on 14 Landsat-satellite scenes with statewide TSI (SDT) for two time periods (2003– 05 and 2007–08). Specifically, the methods were (1) satellitedata processing techniques to remove areas affected by clouds, cloud shadows, haze, shoreline, and dense vegetation for inland lakes greater than 20 acres in Michigan; (2) comparison of the previous method for producing a single open-water predicted TSI (SDT) value (which was based on an area of interest (AOI) and lake-average approach) to an alternative Gethist method for identifying open-water areas in inland lakes (which follows the initial satellite-data processing and targets the darkest pixels, representing the deepest water, before regression equations are created); and (3) checking to see whether the predicted TSI (SDT) values compared well between two regression equations, one previously used in Michigan and an alternative equation from the hydrologic literature. The combination of improved satellite-data processing techniques and the Gethist method to identify open-water areas in inland lakes during 2003–05 and 2007–08 provided a stronger relation and statistical significance between predicted TSI (SDT) and measured TSI than did the AOI lake-average method; differences in results for the two methods were significant at the 99-percent confidence level. With regard to the comparison of the regression equations, there were no statistically significant differences at the 95-percent confidence level between results from the two equations. The previously used equation, in combination with the Gethist method, yielded coefficient of determination (R 2 ) values of 0.71 and 0.77 for the periods 2003–05 and 2007–08, respectively. The alternative equation, in combination with the Gethist method, yielded R 2 values of 0.74 and 0.75 for 2003–05 and 2007–08, respectively. Predicted TSI (SDT) and measured TSI (SDT) values for lakes used in the regression equations compared well, with R 2 values of 0.95 and 0.96 for predicted TSI (SDT) for 2003–05 and 2007–08, respectively. The R 2 values for statewide predicted TSI (SDT) for all inland lakes with available open-water areas for 2003–05 and 2007–08 were 0.91 and 0.93, respectively. Although the two equations predicted similar trophic-state classes, the alternative equation is planned to be used for future prediction of TSI (SDT) values for Michigan inland lakes, to promote consistency in comparing predicted values between States and for potential use in trend analysis.

Michigan↗

Questa baseline and pre-mining ground-water quality investigation. 19. Leaching characteristics of composited materials from mine waste-rock piles and naturally altered areas near Questa, New Mexico

The goal of this study is to compare and contrast the leachability of metals and the acidity from individual mine waste-rock piles and natural erosional scars in the study area near Questa, New Mexico. Surficial multi-increment (composite) samples less than 2 millimeters in diameter from five waste-rock piles, nine erosional-scar areas, a less-altered site, and a tailings slurry-pipe sample were analyzed for bulk chemistry and mineralogy and subjected to two back-to-back leaching procedures. The first leaching procedure, the U.S. Geological Survey Field Leach Test (FLT), is a short-duration leach (5-minute shaking and 10-minute settling) and is intended to leach readily soluble materials. The FLT was immediately followed by an 18-hour, end-over-end rotation leaching procedure. Comparison of results from the back-to-back leaching procedures can provide information about reactions that may take place upon migration of leachates through changing geochemical conditions (for example, pH changes), both within the waste-rock and scar materials and away from the source materials. For the scar leachates, the concentrations of leachable metals varied substantially between the scar areas sampled. The scar leachates have low pH (pH 3.2-4.1). Under these low-pH conditions, cationic metals are solubilized and mobile, but anionic species, such as molybdenum, are less soluble and less mobile. Generally, metal concentrations in the waste-rock leachates did not exceed the upper range of those metal concentrations in the erosional-scar leachates. One exception is molybdenum, which is notably higher in the waste-rock leachates compared with the scar leachates. Most of the waste-rock leachates were at least mildly acidic (pH 3.0-6.2). The pH values in the waste-rock leachates span a large pH range that includes some pH-dependent solubility and metal-attenuation reactions. An increase in pH with leaching time and agitation indicates that there is pH-buffering capacity in some of the waste-rock piles. As pH increased in the waste-pile leachates, concentrations of several metals decreased with increasing time and agitation. Similar pH-dependent reactions may take place upon migration of the leachates in the waste-rock piles. Bulk chemistry, mineralogy, and leachate sulfur-isotope data indicate that the Capulin and Sugar Shack West waste-rock piles are compositionally different from the younger Sugar Shack South, Sugar Shack Middle, and Old Sulphur Gulch piles. The Capulin and Sugar Shack West piles have the lowest-pH leachates (pH 3.0-4.1) of the waste-pile samples, and the source material for the Capulin and Sugar Shack West piles appears to be similar to the source material for the erosional-scar areas. Calcite dissolution, in addition to gypsum dissolution, appears to produce the calcium and sulfate concentrations in leachates from the Sugar Shack South, Sugar Shack Middle, and Old Sulphur Gulch piles.

Scientific Investigations Report↗

Stratiform chromite deposit model

A new descriptive stratiform chromite deposit model was prepared which will provide a framework for understanding the characteristics of stratiform chromite deposits worldwide. Previous stratiform chromite deposit models developed by the U.S. Geological Survey (USGS) have been referred to as Bushveld chromium, because the Bushveld Complex in South Africa is the only stratified, mafic-ultramafic intrusion presently mined for chromite and is the most intensely researched. As part of the on-going effort by the USGS Mineral Resources Program to update existing deposit models for the upcoming national mineral resource assessment, this revised stratiform chromite deposit model includes new data on the geological, mineralogical, geophysical, and geochemical attributes of stratiform chromite deposits worldwide. This model will be a valuable tool in future chromite resource and environmental assessments and supplement previously published models used for mineral resource evaluation. Stratiform chromite deposits are found throughout the world, but the chromitite seams of the Bushveld Complex, South Africa, are the largest and most intensely researched. The chromite ore is located primarily in massive chromitite seams and, less abundantly, in disseminated chromite-bearing layers, both of which occur in the ultramafic section of large, layered mafic-ultramafic stratiform complexes. These mafic-ultramafic intrusions mainly formed in stable cratonic settings or during rift-related events during the Archean or early Proterozoic, although exceptions exist. The chromitite seams are cyclic in nature as well as laterally contiguous throughout the entire intrusion. Gangue minerals include olivine, pyroxenes (orthopyroxene and clinopyroxene), plagioclase, sulfides (pyrite, chalcopyrite, pyrrhotite, pentlandite, bornite), platinum group metals (mainly laurite, cooperite, braggite), and alteration minerals. A few deposits also contain rutile and ilmenite. The alteration phases include serpentine, chlorite, talc, magnetite, kaemmererite, uvarovite, hornblende, and carbonate minerals, such as calcite and dolomite. Stratiform chromite deposits are primarily hosted by peridotites, harzburgites, dunites, pyroxenites, troctolites, and anorthosites. Although metamorphism may have altered the ultramafic regions of layered intrusions postdeposition, only igneous processes are responsible for formation. From a diagnostic standpoint and for assessment purposes, they have no temporal or spatial relation to sedimentary rocks. The exact mechanisms responsible for the development of stratiform chromite deposits and the large, layered mafic-ultramafic intrusions where they are found are highly debated. The leading argument postulates that a parent magma mixed with a more primitive magma during magma chamber recharge. The partially differentiated magma could then be forced into the chromite stability field, resulting in the massive chromitite layers found in stratiform complexes. Contamination of the parent magma by localized assimilation of felsic country rock at the roof of the magma chamber has also been proposed as a mechanism of formation. Others suggest that changes in pressure or oxygen fugacity may be responsible for the occurrence of massive chromitite seams in layered mafic, ultramafic intrusions. The massive chromitite layers contain high levels of chromium and strong associations with platinum group elements. Anomalously high magnesium concentrations as well as low sodium, potassium, and phosphorus concentrations are also important geochemical features of stratiform chromite deposits. The presence of orthopyroxenite in many of the deposits suggests high silica and high magnesium concentrations in the parent magma. Most environmental concerns associated with the mining and processing of chromite ore focus on the solubility of chromium and its oxidation state. Although trivalent chromium (Cr 3+ ) is an essential micronutrient for humans, hexavalent chromium (Cr 6+ ) is highly toxic. Chromium-bearing solid phases that occur in the chromite ore-processing residue, for example, can effect the geochemical behavior and oxidation state of chromium in the environment.

Scientific Investigations Report↗

Organic wastewater compounds, pharmaceuticals, and coliphage in ground water receiving discharge from onsite wastewater treatment systems near La Pine, Oregon: Occurrence and implications for transport

The occurrence of organic wastewater compounds (components of 'personal care products' and other common household chemicals), pharmaceuticals (human prescription and nonprescription medical drugs), and coliphage (viruses that infect coliform bacteria, and found in high concentrations in municipal wastewater) in onsite wastewater (septic tank effluent) and in a shallow, unconfined, sandy aquifer that serves as the primary source of drinking water for most residents near La Pine, Oregon, was documented. Samples from two types of observation networks provided basic occurrence data for onsite wastewater and downgradient ground water. One observation network was a group of 28 traditional and innovative (advanced treatment) onsite wastewater treatment systems and associated downgradient drainfield monitoring wells, referred to as the 'innovative systems network'. The drainfield monitoring wells were located adjacent to or under onsite wastewater treatment system drainfield lines. Another observation network, termed the 'transect network', consisted of 31 wells distributed among three transects of temporary, stainless-steel-screened, direct-push monitoring wells installed along three plumes of onsite wastewater. The transect network, by virtue of its design, also provided a basis for increased understanding of the transport of analytes in natural systems. Coliphage were frequently detected in onsite wastewater. Coliphage concentrations in onsite wastewater were highly variable, ranging from less than 1 to 3,000,000 plaque forming units per 100 milliliters. Coliphage were occasionally detected (eight occurrences) at low concentrations in samples from wells located downgradient from onsite wastewater treatment system drainfield lines. However, coliphage concentrations were below method detection limits in replicate or repeat samples collected from the eight sites. The consistent absence of coliphage detections in the replicate or repeat samples is interpreted to indicate that the detections reported for ground-water samples represented low-level field or laboratory contamination, and it would appear that coliphage were effectively attenuated to less than 1 PFU/100 mL over distances of several feet of transport in the La Pine aquifer and (or) overlying unsaturated zone. Organic wastewater compounds were frequently detected in onsite wastewater. Of the 63 organic wastewater compounds in the analytical schedule, 45 were detected in the 21 samples of onsite wastewater. Concentrations of organic wastewater compounds reached a maximum of 1,300 ug/L (p-cresol). Caffeine was detected at concentrations as high as 320 ug/L. Fourteen of the 45 compounds were detected in more than 90 percent of onsite wastewater samples. Fewer (nine) organic wastewater compounds were detected in ground water, despite the presence of nitrate and chloride likely from onsite wastewater sources. The nine organic wastewater compounds that were detected in ground-water samples were acetyl-hexamethyl-tetrahydro-naphthalene (AHTN), caffeine, cholesterol, hexahydrohexamethyl-cyclopentabenzopyran, N,N-diethyl-meta-toluamide (DEET), tetrachloroethene, tris (2-chloroethyl) phosphate, tris (dichloroisopropyl) phosphate, and tributyl phosphate. Frequent detection of household-chemical type organic wastewater compounds in onsite wastewater provides evidence that some of these organic wastewater compounds may be useful indicators of human waste effluent dispersal in some hydrologic environments. The occurrence of organic wastewater compounds in ground water downgradient from onsite wastewater treatment systems demonstrates that a subgroup of organic wastewater compounds is transported in the La Pine aquifer. The consistently low concentrations (generally less than 1 ug/L) of organic wastewater compounds in water samples collected from wells located no more than 19 feet from drainfield lines indicates that the reactivity (sorption, degradation) of this suite of organic waste.

Oregon↗

Characteristics and Classification of Least Altered Streamflows in Massachusetts

Streamflow records from 85 streamflow-gaging stations at which streamflows were considered to be least altered were used to characterize natural streamflows within southern New England. Period-of-record streamflow data were used to determine annual hydrographs of median monthly flows. The shapes and magnitudes of annual hydrographs of median monthly flows, normalized by drainage area, differed among stations in different geographic areas of southern New England. These differences were gradational across southern New England and were attributed to differences in basin and climate characteristics. Period-of-record streamflow data were also used to analyze the statistical properties of daily streamflows at 61 stations across southern New England by using L-moment ratios. An L-moment ratio diagram of L-skewness and L-kurtosis showed a continuous gradation in these properties between stations and indicated differences between base-flow dominated and runoff-dominated rivers. Streamflow records from a concurrent period (1960-2004) for 61 stations were used in a multivariate statistical analysis to develop a hydrologic classification of rivers in southern New England. Missing records from 46 of these stations were extended by using a Maintenance of Variation Extension technique. The concurrent-period streamflows were used in the Indicators of Hydrologic Alteration and Hydrologic Index Tool programs to determine 224 hydrologic indices for the 61 stations. Principal-components analysis (PCA) was used to reduce the number of hydrologic indices to 20 that provided nonredundant information. The PCA also indicated that the major patterns of variability in the dataset are related to differences in flow variability and low-flow magnitude among the stations. Hierarchical cluster analysis was used to classify stations into groups with similar hydrologic properties. The cluster analysis classified rivers in southern New England into two broad groups: (1) base-flow dominated rivers, whose statistical properties indicated less flow variability and high magnitudes of low flow, and (2) runoff-dominated rivers, whose statistical properties indicated greater flow variability and lower magnitudes of low flow. A four-cluster classification further classified the runoff-dominated streams into three groups that varied in gradient, elevation, and differences in winter streamflow conditions: high-gradient runoff-dominated rivers, northern runoff-dominated rivers, and southern runoff-dominated rivers. A nine-cluster division indicated that basin size also becomes a distinguishing factor among basins at finer levels of classification. Smaller basins (less than 10 square miles) were classified into different groups than larger basins. A comparison of station classifications indicated that a classification based on multiple hydrologic indices that represent different aspects of the flow regime did not result in the same classification of stations as a classification based on a single type of statistic such as a monthly median. River basins identified by the cluster analysis as having similar hydrologic properties tended to have similar basin and climate characteristics and to be in close proximity to one another. Stations were not classified in the same cluster on the basis of geographic location alone; as a result, boundaries cannot be drawn between geographic regions with similar streamflow characteristics. Rivers with different basin and climate characteristics were classified in different clusters, even if they were in adjacent basins or upstream and downstream within the same basin.

Scientific Investigations Report↗