Geology ReportsSearch

Geology topics

Gregory E. Granato

Publications and source records attributed to Gregory E. Granato.

At least 37 records · Page 2Linked to original sources

Compilation of streamflow statistics calculated from daily mean streamflow data collected during water years 1901–2015 for selected U.S. Geological Survey streamgages

Streamflow statistics are needed by decision makers for many planning, management, and design activities. The U.S. Geological Survey (USGS) StreamStats Web application provides convenient access to streamflow statistics for many streamgages by accessing the underlying StreamStatsDB database. In 2016, non-interpretive streamflow statistics were compiled for streamgages located throughout the Nation and stored in StreamStatsDB for use with StreamStats and other applications. Two previously published USGS computer programs that were designed to help calculate streamflow statistics were updated to better support StreamStats as part of this effort. These programs are named “GNWISQ” (Get National Water Information System Streamflow (Q) files), updated to version 1.1.1, and “QSTATS” (Streamflow (Q) Statistics), updated to version 1.1.2. Statistics for 20,438 streamgages that had 1 or more complete years of record during water years 1901 through 2015 were calculated from daily mean streamflow data; 19,415 of these streamgages were within the conterminous United States. About 89 percent of the 20,438 streamgages had 3 or more years of record, and about 65 percent had 10 or more years of record. Drainage areas of the 20,438 streamgages ranged from 0.01 to 1,144,500 square miles. The magnitude of annual average streamflow yields (streamflow per square mile) for these streamgages varied by almost six orders of magnitude, from 0.000029 to 34 cubic feet per second per square mile. About 64 percent of these streamgages did not have any zero-flow days during their available period of record. The 18,122 streamgages with 3 or more years of record were included in the StreamStatsDB compilation so they would be available via the StreamStats interface for user-selected streamgages. All the statistics are available in a USGS ScienceBase data release.

Open-File Report

Estimating risks for water-quality exceedances of total-copper from highway and urban runoff under predevelopment and current conditions with the Stochastic Empirical Loading and Dilution Model (SELDM)

The stochastic empirical loading and dilution model (SELDM) was used to demonstrate methods for estimating risks for water-quality exceedances of event-mean concentrations (EMCs) of total-copper. Monte Carlo methods were used to simulate stormflow, total-hardness, suspended-sediment, and total-copper EMCs as stochastic variables. These simulations were done for the Charles River Basin upstream of Interstate 495 in Bellingham, Massachusetts. The hydrology and water quality of this site were simulated with SELDM by using data from nearby, hydrologically similar sites. Three simulations were done to assess the potential effects of the highway on receiving-water quality with and without highway-runoff treatment by a structural best-management practice (BMP). In the low-development scenario, total copper in the receiving stream was simulated by using a sediment transport curve, sediment chemistry, and sediment-water partition coefficients. In this scenario, neither the highway runoff nor the BMP effluent caused concentration exceedances in the receiving stream that exceed the once in three-year threshold (about 0.54 percent). In the second scenario, without the highway, runoff from the large urban areas in the basin caused exceedances in the receiving stream in 2.24 percent of runoff events. In the third scenario, which included the effects of the urban runoff, neither the highway runoff nor the BMP effluent increased the percentage of exceedances in the receiving stream. Comparison of the simulated geometric mean EMCs with data collected at a downstream monitoring site indicates that these simulated values are within the 95-percent confidence interval of the geometric mean of the measured EMCs.

Conference Paper

Estimating total maximum daily loads with the Stochastic Empirical Loading and Dilution Model

The Massachusetts Department of Transportation (DOT) and the Rhode Island DOT are assessing and addressing roadway contributions to total maximum daily loads (TMDLs). Example analyses for total nitrogen, total phosphorus, suspended sediment, and total zinc in highway runoff were done by the U.S. Geological Survey in cooperation with FHWA to simulate long-term annual loads for TMDL analyses with the stochastic empirical loading and dilution model known as SELDM. Concentration statistics from 19 highway runoff monitoring sites in Massachusetts were used with precipitation statistics from 11 long-term monitoring sites to simulate long-term pavement yields (loads per unit area). Highway sites were stratified by traffic volume or surrounding land use to calculate concentration statistics for rural roads, low-volume highways, high-volume highways, and ultraurban highways. The median of the event mean concentration statistics in each traffic volume category was used to simulate annual yields from pavement for a 29- or 30-year period. Long-term average yields for total nitrogen, phosphorus, and zinc from rural roads are lower than yields from the other categories, but yields of sediment are higher than for the low-volume highways. The average yields of the selected water quality constituents from high-volume highways are 1.35 to 2.52 times the associated yields from low-volume highways. The average yields of the selected constituents from ultraurban highways are 1.52 to 3.46 times the associated yields from high-volume highways. Example simulations indicate that both concentration reduction and flow reduction by structural best management practices are crucial for reducing runoff yields.

Transportation Research Record

Methods for evaluating potential sources of chloride in surface waters and groundwaters of the conterminous United States

Chloride exists as a major ion in most natural waters, but many anthropogenic sources are increasing concentrations of chloride in many receiving waters. Although natural concentrations in continental waters can be as high as 200,000 milligrams per liter, chloride concentrations that are suitable for freshwater ecology, human consumption, and agricultural and industrial water uses commonly are on the order of 10 to 1,000 milligrams per liter. “Road salt” frequently is identified as the sole source of anthropogenic chloride, but only about 30 percent of the salt consumed and released to the environment is used for deicing. Furthermore, several studies in Southern States where the use of deicing salt is minimal also show anthropogenic chloride in rising concentrations and in strong correlation to imperviousness and road density. This is because imperviousness and road density also are strongly correlated to population density. The term “road salt” is a misnomer because deicers applied to parking lots, sidewalks, and driveways can be a substantial source of chloride in some catchments because these land covers are comparable to roadways as a percentage of the total impervious area and commonly receive higher salt application rates than some roadways. Other sources of anthropogenic chloride include wastewater, dust control on unpaved roads, fertilizer, animal waste, irrigation, aquaculture, energy production wastes, and landfill leachates. The assumption that rising chloride concentrations in surface water or groundwater is indicative of contamination by deicing chemicals rather than one or more other potential sources may preclude the identification of toxic, carcinogenic, mutagenic, or endocrine-disrupting contaminants that are associated with many sources of elevated chloride concentrations. Once the sources of anthropogenic chloride in an area of interest have been identified and measured, water and solute budgets can be estimated to guide decisionmakers to identify and apply potential mitigation measures that can reduce the problem. Scientists, engineers, regulators, and decisionmakers need information about potential sources of chloride, water and solute budgets, and methods for collecting water-quality data to help identify potential sources. This information is needed to evaluate potential sources of chloride in areas where chloride may have adverse ecological effects or may degrade water supplies used for drinking water, agriculture, or industry. Knowledge of potential sources will help decisionmakers identify the best mitigation measures to reduce the total background chloride load, thereby reducing the potential for water-quality exceedances that occur because of superposition on rising background concentrations. Also, knowledge of potential sources may help decisionmakers identify the potential for the presence of contaminants that have toxic, carcinogenic, mutagenic, or endocrine-disrupting effects at concentrations that are lower by orders of magnitude than the chloride concentrations in the source water. This report is a comprehensive synthesis of relevant information, but it is not the result of an exhaustive search for literature on each topic. The potential adverse effects of chloride on infrastructure and the environment are not discussed in this report because these issues have been extensively documented elsewhere.

Open-File Report

A case study demonstrating analysis of stormflows, concentrations, and loads of nutrients in highway runoff and swale discharge with the Stochastic Empirical Loading and Dilution Model (SELDM)

Decisionmakers need information about the quality and quantity of stormwater runoff, the risk for adverse effects of runoff on receiving waters, and the potential effectiveness of mitigation measures to reduce these risks. The Stochastic Empirical Loading and Dilution Model (SELDM) uses Monte Carlo methods to generate stormflows, concentrations, and loads from a highway site and an upstream basin to provide needed risk-based information. SELDM was designed to help inform water-management decisions for streams and lakes receiving runoff from a highway or other land-use site. The purpose of this paper is to provide a brief description of SELDM and a hypothetical case study demonstrating the type of risk-based information that SELDM can provide. Total nitrogen (TN) and total phosphorus (TP) were selected as example constituents because nutrients are a common concern throughout the Nation and data for receiving waters, highway runoff, and the performance of best management practices (BMPs) are readily available for these constituents. The case study is hypothetical, but was formulated by using actual data from selected monitoring sites in New England. Data representing streamflow and water-quality were collected at U.S. Geological Survey (USGS) streamgage 01208950 Sasco Brook near Southport, CT, which has a drainage area of 7.38 square miles. In this hypothetical case study a 4-lane highway would replace the current 2-lane road and would have a contributing area of 2.2 acres between the topographic basin divides. Concentrations of TN and TP in highway runoff were simulated with data from USGS highway-runoff monitoring station 423027071291301 along State Route 2 in Littleton Massachusetts. Results of a highway-runoff analysis are shown in relation to three hypothetical discharge criteria for TN and two hypothetical discharge criteria for TP. The risks for exceeding TN discharge criteria of 3, 5, and 8 mg/L for highway runoff are 7.4, 0.83, and 0.13 percent of 1,721 runoff events that may occur during a stochastic 30-year simulation. If a grassy swale is used to treat the runoff, the risks for TN exceedances are reduced to 3.2, 0.33 and 0.03 percent, respectively. The risks for exceeding TP discharge criteria of 0.1 and 0.5 mg/L for highway runoff are 49 and 1.2 percent, respectively. If a grassy swale is used to treat the runoff, the risks for TP exceedances are 57 and 0.8 percent, respectively. The risks for the 0.1 mg/L criterion increase because swales can be a source of TP if pavement concentrations are low. The risks for the 0.5 mg/L criterion decrease because the swale is effective for reducing high TP concentrations. Although the results are mixed for storm-event concentrations, the grassy swale effectively reduces annual loads. Annual loads from the swale are, on average, about 49 percent of highway loads for TN and 62 percent of highway loads of TP because the swale reduces high runoff concentrations and stormflow volumes. Analysis of upstream and downstream concentrations indicates that runoff from the site of interest does not have a substantial effect on instream stormflow concentrations in this example simulation.

Conference Paper

Estimating the risks for adverse effects of total phosphorus in receiving streams with the Stochastic Empirical Loading and Dilution Model (SELDM)

Studies from North Carolina (NC) indicate that increasing concentrations of total phosphorus (TP) and other constituents are correlated to adverse effects on stream ecosystems as evidenced by differences in benthic macroinvertebrate populations in streams across the state. As a result, stringent in-stream criteria based on the Water Quality Assessed by Benthic macroinvertebrate health ratings (WQABI) have been proposed for regulating TP concentrations in stormwater discharges and for selecting stormwater best management practices (BMPs). The WQABI criteria concentrations may not be suitable for evaluating stormwater discharges because they are based on baseflow concentration statistics, the criteria do not include a clearly defined allowable exceedance frequency, and there are substantial uncertainties in estimating the quality of runoff, BMP discharge, and receiving waters for sites without monitoring data. The Stochastic Empirical Loading and Dilution Model (SELDM), which was developed by the U.S. Geological Survey in cooperation with the Federal Highway Administration, was used to simulate the quality of runoff, BMP discharge, and receiving waters to evaluate risks for water-quality exceedances with different criteria concentrations, allowable exceedance frequencies, and selected water-quality statistics. Water-quality data from two neighboring basins in the Piedmont ecoregion in NC were used to simulate in-stream stormwater quality. Data collected at 15 sites in NC were used to simulate runoff quality. Statistics for stochastic modeling of volume reduction, hydrograph extension, and water-quality treatment by BMPs, were used to simulate potential effect of these treatments on discharge quality and downstream stormwater quality. Results of these long-term 30-year simulations were used to evaluate criteria concentrations, the potential frequency of water-quality exceedances, and the effect of data selection on risks for water-quality exceedances. The simulations indicate that the potential frequency for exceeding instream and stormwater discharge criteria depend on the detailed definition of the criteria and the data that are selected for simulating water quality. Data and simulation results indicate that the baseflow concentrations do not represent stormwater concentrations, even in predominantly forested basins. There is substantial uncertainty in applying stormwater statistics to unmonitored sites, even if these statistics are applied to neighboring basins such as in this example. Over a period of several years (or more) it would be impossible to meet many of the proposed instream and stormwater discharge quality criteria unless these criteria include an allowable exceedance frequency because stormwater concentrations commonly vary by orders of magnitude. Selection of BMPs by using concentration reduction as the sole criteria may underestimate potential benefits of BMPs that also provide volume reduction, which reduces discharge loads, and hydrograph extension, which increases the dilution of runoff into a larger proportion of the upstream stormflow. Results of this study indicate the potential benefits of the multi-decade simulations that SELDM provides because these simulations quantify risks and uncertainties that affect decisions made with available data and statistics. Results of the SELDM simulations indicate that the WQABI criteria concentrations may be too stringent for evaluating the stormwater quality in receiving streams, highway runoff, and BMP discharges; especially with the substantial uncertainties inherent in selecting representative data.

Conference Paper

Stochastic empirical loading and dilution model for analysis of flows, concentrations, and loads of highway runoff constituents

In cooperation with FHWA, the U.S. Geological Survey developed the stochastic empirical loading and dilution model (SELDM) to supersede the 1990 FHWA runoff quality model. The SELDM tool is designed to transform disparate and complex scientific data into meaningful information about the adverse risks of runoff on receiving waters, the potential need for mitigation measures, and the potential effectiveness of such measures for reducing such risks. The SELDM tool is easy to use because much of the information and data needed to run it are embedded in the model and obtained by defining the site location and five simple basin properties. Information and data from thousands of sites across the country were compiled to facilitate the use of the SELDM tool. A case study illustrates how to use the SELDM tool for conducting the types of sensitivity analyses needed to properly assess water quality risks. For example, the use of deterministic values to model upstream stormflows instead of representative variations in prestorm flow and runoff may substantially overestimate the proportion of highway runoff in downstream flows. Also, the risks for total phosphorus excursions are substantially affected by the selected criteria and the modeling methods used. For example, if a single deterministic concentration is used rather than a stochastic population of values to model upstream concentrations, then the percentage of water quality excursions in the downstream receiving waters may depend entirely on the selected upstream concentration.

Transportation Research Record

Hydrologic Drought Decision Support System (HyDroDSS)

The hydrologic drought decision support system (HyDroDSS) was developed by the U.S. Geological Survey (USGS) in cooperation with the Rhode Island Water Resources Board (RIWRB) for use in the analysis of hydrologic variables that may indicate the risk for streamflows to be below user-defined flow targets at a designated site of interest, which is defined herein as data-collection site on a stream that may be adversely affected by pumping. Hydrologic drought is defined for this study as a period of lower than normal streamflows caused by precipitation deficits and (or) water withdrawals. The HyDroDSS is designed to provide water managers with risk-based information for balancing water-supply needs and aquatic-habitat protection goals to mitigate potential effects of hydrologic drought. This report describes the theory and methods for retrospective streamflow-depletion analysis, rank correlation analysis, and drought-projection analysis. All three methods are designed to inform decisions made by drought steering committees and decisionmakers on the basis of quantitative risk assessment. All three methods use estimates of unaltered streamflow, which is the measured or modeled flow without major withdrawals or discharges, to approximate a natural low-flow regime. Retrospective streamflow-depletion analysis can be used by water-resource managers to evaluate relations between withdrawal plans and the potential effects of withdrawal plans on streams at one or more sites of interest in an area. Retrospective streamflow-depletion analysis indicates the historical risk of being below user-defined flow targets if different pumping plans were implemented for the period of record. Retrospective streamflow-depletion analysis also indicates the risk for creating hydrologic drought conditions caused by use of a pumping plan. Retrospective streamflow-depletion analysis is done by calculating the net streamflow depletions from withdrawals and discharges and applying these depletions to a simulated record of unaltered streamflow. Rank correlation analysis in the HyDroDSS indicates the persistence of hydrologic measurements from month to month for the prediction of developing hydrologic drought conditions and quantitatively indicates which hydrologic variables may be used to indicate the onset of hydrologic drought conditions. Rank correlation analysis also indicates the potential use of each variable for estimating the monthly minimum unaltered flow at a site of interest for use in the drought-projection analysis. Rank correlation analysis in the HyDroDSS is done by calculating Spearman’s rho for paired samples and the 95-percent confidence limits of this rho value. Rank correlation analysis can be done by using precipitation, groundwater levels, measured streamflows, and estimated unaltered streamflows. Serial correlation analysis, which indicates relations between current and future values, can be done for a single site. Cross correlation analysis, which indicates relations among current values at one site and current and future values at a second site, also can be done. Drought-projection analysis in the HyDroDSS indicates the risk for being in a hydrologic drought condition during the current month and the five following months with and without pumping. Drought-projection analysis also indicates the potential effectiveness of water-conservation methods for mitigating the effect of withdrawals in the coming months on the basis of the amount of depletion caused by different pumping plans and on the risk of unaltered flows being below streamflow targets. Drought-projection analysis in the HyDroDSS is done with Monte Carlo methods by using the position analysis method. In this method the initial value of estimated unaltered streamflows is calculated by correlation to a measured hydrologic variable (monthly precipitation, groundwater levels, or streamflows from an index station identified with the rank correlation analysis). Then a pseudorandom number generator is used to create 251 six-month-long flow traces by using a bootstrap method. Serial correlation of the estimated unaltered monthly minimum streamflows determined from the rank correlation analysis is preserved within each flow trace. The sample of unaltered streamflows indicates the risk of being below flow targets in the coming months under simulated natural conditions (without historic withdrawals). The streamflow-depletion algorithms are then used to estimate risks of flow being below targets if selected pumping plans are used. This report also describes the implementation of the HyDroDSS. The HyDroDSS was developed as a Microsoft Access® database application to facilitate storage, handling, and use of hydrologic datasets with a simple graphical user interface. The program is implemented in the database by using the Visual Basic for Applications® (VBA) programming language. Program source code for the analytical techniques is provided in the HyDroDSS and in electronic text files accompanying this report. Program source code for the graphical user interface and for data-handling code, which is specific to Microsoft Access® and the HyDroDSS, is provided in the database. An installation package with a run-time version of the software is available with this report for potential users who do not have a compatible copy of Microsoft Access®. Administrative rights are needed to install this version of the HyDroDSS. A case study, to demonstrate the use of HyDroDSS and interpretation of results for a site of interest, is detailed for the USGS streamgage on the Hunt River (station 01117000) near East Greenwich in central Rhode Island. The Hunt River streamgage was used because it has a long record of streamflow and is in a well-studied basin with a substantial amount of hydrologic and water-use data including groundwater pumping for municipal water supply.

Rhode Island

Statistics for stochastic modeling of volume reduction, hydrograph extension, and water-quality treatment by structural stormwater runoff best management practices (BMPs)

The U.S. Geological Survey (USGS) developed the Stochastic Empirical Loading and Dilution Model (SELDM) in cooperation with the Federal Highway Administration (FHWA) to indicate the risk for stormwater concentrations, flows, and loads to be above user-selected water-quality goals and the potential effectiveness of mitigation measures to reduce such risks. SELDM models the potential effect of mitigation measures by using Monte Carlo methods with statistics that approximate the net effects of structural and nonstructural best management practices (BMPs). In this report, structural BMPs are defined as the components of the drainage pathway between the source of runoff and a stormwater discharge location that affect the volume, timing, or quality of runoff. SELDM uses a simple stochastic statistical model of BMP performance to develop planning-level estimates of runoff-event characteristics. This statistical approach can be used to represent a single BMP or an assemblage of BMPs. The SELDM BMP-treatment module has provisions for stochastic modeling of three stormwater treatments: volume reduction, hydrograph extension, and water-quality treatment. In SELDM, these three treatment variables are modeled by using the trapezoidal distribution and the rank correlation with the associated highway-runoff variables. This report describes methods for calculating the trapezoidal-distribution statistics and rank correlation coefficients for stochastic modeling of volume reduction, hydrograph extension, and water-quality treatment by structural stormwater BMPs and provides the calculated values for these variables. This report also provides robust methods for estimating the minimum irreducible concentration (MIC), which is the lowest expected effluent concentration from a particular BMP site or a class of BMPs. These statistics are different from the statistics commonly used to characterize or compare BMPs. They are designed to provide a stochastic transfer function to approximate the quantity, duration, and quality of BMP effluent given the associated inflow values for a population of storm events. A database application and several spreadsheet tools are included in the digital media accompanying this report for further documentation of methods and for future use. In this study, analyses were done with data extracted from a modified copy of the January 2012 version of International Stormwater Best Management Practices Database, designated herein as the January 2012a version. Statistics for volume reduction, hydrograph extension, and water-quality treatment were developed with selected data. Sufficient data were available to estimate statistics for 5 to 10 BMP categories by using data from 40 to more than 165 monitoring sites. Water-quality treatment statistics were developed for 13 runoff-quality constituents commonly measured in highway and urban runoff studies including turbidity, sediment and solids; nutrients; total metals; organic carbon; and fecal coliforms. The medians of the best-fit statistics for each category were selected to construct generalized cumulative distribution functions for the three treatment variables. For volume reduction and hydrograph extension, interpretation of available data indicates that selection of a Spearman’s rho value that is the average of the median and maximum values for the BMP category may help generate realistic simulation results in SELDM. The median rho value may be selected to help generate realistic simulation results for water-quality treatment variables. MIC statistics were developed for 12 runoff-quality constituents commonly measured in highway and urban runoff studies by using data from 11 BMP categories and more than 167 monitoring sites. Four statistical techniques were applied for estimating MIC values with monitoring data from each site. These techniques produce a range of lower-bound estimates for each site. Four MIC estimators are proposed as alternatives for selecting a value from among the estimates from multiple sites. Correlation analysis indicates that the MIC estimates from multiple sites were weakly correlated with the geometric mean of inflow values, which indicates that there may be a qualitative or semiquantitative link between the inflow quality and the MIC. Correlations probably are weak because the MIC is influenced by the inflow water quality and the capability of each individual BMP site to reduce inflow concentrations.

Scientific Investigations Report

Assessing potential effects of highway runoff on receiving-water quality at selected sites in Oregon with the Stochastic Empirical Loading and Dilution Model (SELDM)

In 2012, the U.S. Geological Survey and the Oregon Department of Transportation began a cooperative study to demonstrate use of the Stochastic Empirical Loading and Dilution Model (SELDM) for runoff-quality analyses in Oregon. SELDM can be used to estimate stormflows, constituent concentrations, and loads from the area upstream of a stormflow discharge site, from the site of interest and in the receiving waters downstream of the discharge. SELDM also can be used to assess the potential effectiveness of best management practices (BMP) for mitigating potential effects of runoff in receiving waters. Nominally, SELDM is a highway-runoff model, but it is well suited for analysis of runoff from other land uses as well. This report provides case studies and examples to demonstrate stochastic-runoff modeling concepts and to demonstrate application of the model. Basin characteristics from six Oregon highway study sites were used to demonstrate various applications of the model. The highway catchment and upstream basin drainage areas of these study sites ranged from 3.85 to 11.83 acres and from 0.16 to 6.56 square miles, respectively. The upstream basins of two sites are urbanized, and the remaining four sites are less than 5 percent impervious. SELDM facilitates analysis by providing precipitation, pre-storm streamflow, and other variables by region or from hydrologically similar sites. In Oregon, there can be large variations in precipitation and streamflow among nearby sites. Therefore, spatially interpolated geographic information system data layers containing storm-event precipitation and pre-storm streamflow statistics specific to Oregon were created for the study using Kriging techniques. Concentrations and loads of cadmium, chloride, chromium, copper, iron, lead, nickel, phosphorus, and zinc were simulated at the six Oregon highway study sites by using statistics from sites in other areas of the country. Water‑quality datasets measured at hydrologically similar basins in the vicinity of the study sites in Oregon were selected and compiled to estimate stormflow-quality statistics for the upstream basins. The quality of highway runoff and some upstream stormflow constituents were simulated by using statistical moments (average, standard deviation, and skew) of the logarithms of data. Some upstream stormflow constituents were simulated by using transport curves, which are relations between stormflow and constituent concentrations. Stochastic analyses were done by using SELDM to demonstrate use of the model and to illustrate the types of information that stochastic analyses may provide: 1. An analysis was done to demonstrate use of dilution factors as an initial reconnaissance tool for comparing relative risk among sites. 2. An analysis of hardness-dependent, water-quality criteria was done to illustrate the effects of variations in hardness and flow on the application and interpretation of such criteria. This analysis shows that hardness-dependent criteria can vary by an order of magnitude among storm events because hardness is diluted by stormflows. 3. An analysis of uncertainties in input and output values was done to demonstrate that properly selected robust datasets are needed to represent conditions at a site of interest. This analysis shows that the rate of water-quality exceedances that are measured or simulated may depend on sample size and the luck of the draw. 4. An analysis was done to demonstrate that SELDM and other Monte Carlo models may generate extreme values from input statistics, which may or may not be feasible based on physicochemical or hydrological limits. 5. An analysis of BMP modeling methods was done to demonstrate use of the model for estimating treatment requirements for meeting water-quality objectives. 6. An analysis of the use of grab sampling and nonstochastic upstream modeling methods was done to evaluate the potential effects on modeling outcomes. Additional analyses using surrogate water-quality datasets for the upstream basin and highway catchment were provided for six Oregon study sites to illustrate the risk-based information that SELDM will produce. These analyses show that the potential effects of highway runoff on receiving-water quality downstream of the outfall depends on the ratio of drainage areas (dilution), the quality of the receiving water upstream of the highway, and the concentration of the criteria of the constituent of interest. These analyses also show that the probability of exceeding a water-quality criterion may depend on the input statistics used, thus careful selection of representative values is important.

Oregon

Numerical simulation of groundwater and surface-water interactions in the Big River Management Area, central Rhode Island

The Rhode Island Water Resources Board is considering use of groundwater resources from the Big River Management Area in central Rhode Island because increasing water demands in Rhode Island may exceed the capacity of current sources. Previous water-resources investigations in this glacially derived, valley-fill aquifer system have focused primarily on the effects of potential groundwater-pumping scenarios on streamflow depletion; however, the effects of groundwater withdrawals on wetlands have not been assessed, and such assessments are a requirement of the State’s permitting process to develop a water supply in this area. A need for an assessment of the potential effects of pumping on wetlands in the Big River Management Area led to a cooperative agreement in 2008 between the Rhode Island Water Resources Board, the U.S. Geological Survey, and the University of Rhode Island. This partnership was formed with the goal of developing methods for characterizing wetland vegetation, soil type, and hydrologic conditions, and monitoring and modeling water levels for pre- and post-water-supply development to assess potential effects of groundwater withdrawals on wetlands. This report describes the hydrogeology of the area and the numerical simulations that were used to analyze the interaction between groundwater and surface water in response to simulated groundwater withdrawals. The results of this analysis suggest that, given the hydrogeologic conditions in the Big River Management Area, a standard 5-day aquifer test may not be sufficient to determine the effects of pumping on water levels in nearby wetlands. Model simulations showed water levels beneath Reynolds Swamp declined by about 0.1 foot after 5 days of continuous pumping, but continued to decline by an additional 4 to 6 feet as pumping times were increased from a 5-day simulation period to a simulation period representative of long-term average monthly conditions. This continued decline in water levels with increased pumping time is related to the shift from the primary source of water to the pumped wells being derived from aquifer storage during the early-time (5 days) simulation to being derived more from induced infiltration from the flooded portion of the Big River (southernmost extent of the Flat River Reservoir) during the months of March through October or from captured groundwater discharge to this portion of the Big River when the downstream Flat River Reservoir is drained for weed control during the months of November through February, as was the case for the long-term monthly conditions.

Rhode Island

Stochastic empirical loading and dilution model (SELDM) version 1.0.0

The Stochastic Empirical Loading and Dilution Model (SELDM) is designed to transform complex scientific data into meaningful information about the risk of adverse effects of runoff on receiving waters, the potential need for mitigation measures, and the potential effectiveness of such management measures for reducing these risks. The U.S. Geological Survey developed SELDM in cooperation with the Federal Highway Administration to help develop planning-level estimates of event mean concentrations, flows, and loads in stormwater from a site of interest and from an upstream basin. Planning-level estimates are defined as the results of analyses used to evaluate alternative management measures; planning-level estimates are recognized to include substantial uncertainties (commonly orders of magnitude). SELDM uses information about a highway site, the associated receiving-water basin, precipitation events, stormflow, water quality, and the performance of mitigation measures to produce a stochastic population of runoff-quality variables. SELDM provides input statistics for precipitation, prestorm flow, runoff coefficients, and concentrations of selected water-quality constituents from National datasets. Input statistics may be selected on the basis of the latitude, longitude, and physical characteristics of the site of interest and the upstream basin. The user also may derive and input statistics for each variable that are specific to a given site of interest or a given area. SELDM is a stochastic model because it uses Monte Carlo methods to produce the random combinations of input variable values needed to generate the stochastic population of values for each component variable. SELDM calculates the dilution of runoff in the receiving waters and the resulting downstream event mean concentrations and annual average lake concentrations. Results are ranked, and plotting positions are calculated, to indicate the level of risk of adverse effects caused by runoff concentrations, flows, and loads on receiving waters by storm and by year. Unlike deterministic hydrologic models, SELDM is not calibrated by changing values of input variables to match a historical record of values. Instead, input values for SELDM are based on site characteristics and representative statistics for each hydrologic variable. Thus, SELDM is an empirical model based on data and statistics rather than theoretical physiochemical equations. SELDM is a lumped parameter model because the highway site, the upstream basin, and the lake basin each are represented as a single homogeneous unit. Each of these source areas is represented by average basin properties, and results from SELDM are calculated as point estimates for the site of interest. Use of the lumped parameter approach facilitates rapid specification of model parameters to develop planning-level estimates with available data. The approach allows for parsimony in the required inputs to and outputs from the model and flexibility in the use of the model. For example, SELDM can be used to model runoff from various land covers or land uses by using the highway-site definition as long as representative water quality and impervious-fraction data are available.

Techniques and Methods

Estimating basin lagtime and hydrograph-timing indexes used to characterize stormflows for runoff-quality analysis

A nationwide study to better define triangular-hydrograph statistics for use with runoff-quality and flood-flow studies was done by the U.S. Geological Survey (USGS) in cooperation with the Federal Highway Administration. Although the triangular hydrograph is a simple linear approximation, the cumulative distribution of stormflow with a triangular hydrograph is a curvilinear S-curve that closely approximates the cumulative distribution of stormflows from measured data. The temporal distribution of flow within a runoff event can be estimated using the basin lagtime, (which is the time from the centroid of rainfall excess to the centroid of the corresponding runoff hydrograph) and the hydrograph recession ratio (which is the ratio of the duration of the falling limb to the rising limb of the hydrograph). This report documents results of the study, methods used to estimate the variables, and electronic files that facilitate calculation of variables. Ten viable multiple-linear regression equations were developed to estimate basin lagtimes from readily determined drainage basin properties using data published in 37 stormflow studies. Regression equations using the basin lag factor (BLF, which is a variable calculated as the main-channel length, in miles, divided by the square root of the main-channel slope in feet per mile) and two variables describing development in the drainage basin were selected as the best candidates, because each equation explains about 70 percent of the variability in the data. The variables describing development are the USGS basin development factor (BDF, which is a function of the amount of channel modifications, storm sewers, and curb-and-gutter streets in a basin) and the total impervious area variable (IMPERV) in the basin. Two datasets were used to develop regression equations. The primary dataset included data from 493 sites that have values for the BLF, BDF, and IMPERV variables. This dataset was used to develop the best-fit regression equation using the BLF and BDF variables. The secondary dataset included data from 896 sites that have values for the BLF and IMPERV variables. This dataset was used to develop the best-fit regression equation using the BLF and IMPERV variables. Analysis of hydrograph recession ratios and basin characteristics for 41 sites indicated that recession ratios are random variables. Thus, recession ratios cannot be estimated quantitatively using multiple linear regression equations developed using the data available for these sites. The minimums of recession ratios for different streamgages are well characterized by a value of one. The most probable values and maximum values of recession ratios for different streamgages are, however, more variable than the minimums. The most probable values of recession ratios for the 41 streamgages analyzed ranged from 1.0 to 3.52 and had a median of 1.85. The maximum values ranged from 2.66 to 11.3 and had a median of 4.36.

Scientific Investigations Report

Quality of stormwater runoff discharged from Massachusetts highways, 2005-07

The U.S. Geological Survey (USGS), in cooperation with U.S. Department of Transportation Federal Highway Administration and the Massachusetts Department of Transportation, conducted a field study from September 2005 through September 2007 to characterize the quality of highway runoff for a wide range of constituents. The highways studied had annual average daily traffic (AADT) volumes from about 3,000 to more than 190,000 vehicles per day. Highway-monitoring stations were installed at 12 locations in Massachusetts on 8 highways. The 12 monitoring stations were subdivided into 4 primary, 4 secondary, and 4 test stations. Each site contained a 100-percent impervious drainage area that included two or more catch basins sharing a common outflow pipe. Paired primary and secondary stations were located within a few miles of each other on a limited-access section of the same highway. Most of the data were collected at the primary and secondary stations, which were located on four principal highways (Route 119, Route 2, Interstate 495, and Interstate 95). The secondary stations were operated simultaneously with the primary stations for at least a year. Data from the four test stations (Route 8, Interstate 195, Interstate 190, and Interstate 93) were used to determine the transferability of the data collected from the principal highways to other highways characterized by different construction techniques, land use, and geography. Automatic-monitoring techniques were used to collect composite samples of highway runoff and make continuous measurements of several physical characteristics. Flowweighted samples of highway runoff were collected automatically during approximately 140 rain and mixed rain, sleet, and snowstorms. These samples were analyzed for physical characteristics and concentrations of 6 dissolved major ions, total nutrients, 8 total-recoverable metals, suspended sediment, and 85 semivolatile organic compounds (SVOCs), which include priority polyaromatic hydrocarbons (PAHs), phthalate esters, and other anthropogenic or naturally occurring organic compounds. The distribution of particle size of suspended sediment also was determined for composite samples of highway runoff. Samples of highway runoff were collected year round and under various dry antecedent conditions throughout the 2-year sampling period. In addition to samples of highway runoff, supplemental samples also were collected of sediment in highway runoff, background soils, berm materials, maintenance sands, deicing compounds, and vegetation matter. These additional samples were collected near or on the highways to support data analysis. There were few statistically significant differences between populations of constituent concentrations in samples from the primary and secondary stations on the same principal highways (Mann-Whitney test, 95-percent confidence level). Similarly, there were few statistically significant differences between populations of constituent concentrations for the four principal highways (data from the paired primary and secondary stations for each principal highway) and populations for test stations with similar AADT volumes. Exceptions to this include several total-recoverable metals for stations on Route 2 and Interstate 195 (highways with moderate AADT volumes), and for stations on Interstate 95 and Interstate 93 (highways with high AADT volumes). Supplemental data collected during this study indicate that many of these differences may be explained by the quantity, as well as the quality, of the sediment in samples of highway runoff. Nonparametric statistical methods also were used to test for differences between populations of sample constituent concentrations among the four principal highways that differed mainly in traffic volume. These results indicate that there were few statistically significant differences (Mann-Whitney test, 95-percent confidence level) for populations of concentrations of most total-recoverable metals

Scientific Investigations Report

Methods for development of planning-level estimates of stormflow at unmonitored stream sites in the conterminous United States

This report documents methods for data compilation and analysis of statistics for stormflows that meet data-quality objectives for order-of-magnitude planning-level water-quality estimates at unmonitored sites in the conterminous United States. Statistics for prestorm streamflow, precipitation, and runoff coefficients are used to model stormflows for use with the Stochastic Empirical Loading and Dilution Model (SELDM), which is a highway-runoff model. SELDM is designed to better quantify the risk of exceeding water-quality criteria as precipitation, discharge, ambient water quality, and highway-runoff quality vary from storm to storm. Summary statistics also may be used to help estimate annual-average water-quality loads. Streamflow statistics are used to estimate prestorm flows. Streamflow statistics are estimated by analysis of data from 2,873 U.S. Geological Survey streamgages in the conterminous United States with drainage areas ranging from 10 to 500 square miles and at least 24 years of record during the period 1960−2004. Streamflow statistics are regionalized using U.S. Environmental Protection Agency Level III nutrient ecoregions. Storm-event precipitation statistics are estimated by analysis of data from 2,610 National Oceanic and Atmospheric Administration hourly-precipitation data stations in the conterminous United States with at least 25 years of data during the 1965−2006 period. Storm-event precipitation statistics are regionalized using U.S. Environmental Protection Agency rain zones. Statistics to characterize volumetric runoff coefficients are estimated using data from 6,142 storm events at 306 study sites. Runoff coefficient statistics are not regionalized, but are organized by total impervious area. All of the geographic information system files, computer programs, data files, and regression results developed for this study are included on the CD−ROM accompanying this report.

Report

Effects of Water-Management Strategies on Water Resources in the Pawcatuck River Basin, Southwestern Rhode Island and Southeastern Connecticut

The Pawcatuck River Basin in southwestern Rhode Island and southeastern Connecticut is an important high-quality water resource for domestic and public supplies, irrigation, recreation, and the aquatic ecosystem. Concerns about the effects of water withdrawals on aquatic habitat in the basin have prompted local, State, and Federal agencies to explore water-management strategies that minimize the effects of withdrawals on the aquatic habitat. As part of this process, the U.S. Geological Survey in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service and the Rhode Island Water Resources Board completed a study to assess the effects of current (2000-04) and potential water withdrawals on streamflows and groundwater levels using hydrologic simulation models developed for the basin. The major findings of the model simulations are: *Moving highly variable seasonal irrigation withdrawals from streams to groundwater wells away from streams reduces short-term fluctuations in streamflow and increases streamflow in the summer when flows are lowest. This occurs because of the inherent time lag between when water is withdrawn from the aquifer and when it affects streamflow. *A pumped well in the vicinity of small streams indicates that if withdrawals exceed available streamflow, groundwater levels drop substantially as a consequence of water lost from aquifer storage, which may reduce the time wetlands and vernal pools are saturated, affecting the animal and plant life that depend on these habitats. *The effects of pumping on water resources such as ponds, streams, and wetlands can be minimized by relocating pumping wells, implementing seasonal pumping schemes that utilize different wells and pumping rates, or both. *The effects of projected land-use change, mostly from forest to low- and medium density housing, indicate only minor changes in streamflow at the subbasin scale examined; however, at a local scale, high flows could increase, and low flows could decrease as a result of increased impervious area. In some instances, low flows could increase slightly as a result of decreased evapotranspiration from the loss of deeprooted vegetation (forest) associated with development. *In some subbasins where large areas of agricultural lands were converted to low- and medium-density housing, low flows increase because the consumptive domestic water use was projected to be less than consumptive agricultural water use. All agricultural water use was for irrigation purposes and was assumed to be lost from the basin through evapotranspiration.

Circular

Computer programs for obtaining and analyzing daily mean streamflow data from the U.S. Geological Survey National Water Information System web site

Five computer programs were developed for obtaining and analyzing streamflow from the National Water Information System (NWISWeb). The programs were developed as part of a study by the U.S. Geological Survey, in cooperation with the Federal Highway Administration, to develop a stochastic empirical loading and dilution model. The programs were developed because reliable, efficient, and repeatable methods are needed to access and process streamflow information and data. The first program is designed to facilitate the downloading and reformatting of NWISWeb streamflow data. The second program is designed to facilitate graphical analysis of streamflow data. The third program is designed to facilitate streamflow-record extension and augmentation to help develop long-term statistical estimates for sites with limited data. The fourth program is designed to facilitate statistical analysis of streamflow data. The fifth program is a preprocessor to create batch input files for the U.S. Environmental Protection Agency DFLOW3 program for calculating low-flow statistics. These computer programs were developed to facilitate the analysis of daily mean streamflow data for planning-level water-quality analyses but also are useful for many other applications pertaining to streamflow data and statistics. These programs and the associated documentation are included on the CD-ROM accompanying this report. This report and the appendixes on the CD-ROM describe the implementation and use of the programs and the interpretation of results from the programs. The body of this report provides an overview of the five programs included on this CD-ROM. The appendixes are the software manuals for each program. These manuals describe statistical and numerical methods used to implement each program, input-file formats, output-file formats, installation of the programs, and use of the programs. Each appendix is written as a self-contained manual because each program may have many uses alone or in tandem with other programs on the CD-ROM. Each of these programs uses graphical user interface that follows standard Microsoft Windows interface conventions.

Open-File Report