Search USGSSearch

SEARCH · Search USGS

Results for “Computational Statistics and Data Analysis”

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 253 records · Page 14Linked to original sources

Recharge rates and aquifer hydraulic characteristics for selected drainage basins in middle and east Tennessee

Quantitative information concerning aquifer hydrologic and hydraulic characteristics is needed to manage the development of ground-water resources. These characteristics are poorly defined for the bedrock aquifers in Middle and East Tennessee where demand for water is increasing. This report presents estimates of recharge rate, storage coefficient, diffusivity, and transmissivity for representative drainage basins in Middle and East Tennessee, as determined from analyses of stream-aquifer interactions. The drainage basins have been grouped according to the underlying major aquifer, then statistical descriptions applied to each group, in order to define area1 distribution of these characteristics. Aquifer recharge rates are estimated for representative low, average, and high flow years for 63 drainage basins using hydrograph analysis techniques. Net annual recharge during average flow years for all basins ranges from 4.1 to 16.8 in/yr (inches per year), with a mean value of 7.3 in. In general, recharge rates are highest for basins underlain by the Blue Ridge aquifer (mean value11.7 in/yr) and lowest for basins underlain by the Central Basin aquifer (mean value 5.6 in/yr). Mean recharge values for the Cumberland Plateau, Highland Rim, and Valley and Ridge aquifers are 6.5, 7.4, and 6.6 in/yr, respectively. Gravity drainage characterizes ground-water flow in most surficial bedrock aquifer in Tennessee. Accordingly, a gravity yield analysis, which compares concurrent water-level and streamflow hydrographs, was used to estimate aquifer storage coefficient for nine study basins. The basin estimates range from 0.002 to 0.140; however, most estimates are within a narrow range of values, from 0.01 to 0.025. Accordingly, storage coefficient is estimated to be 0.01 for all aquifers in Middle and East Tennessee, with the exception of the aquifer in the inner part of the Central Basin, for which storage coefficient is estimated to be 0.002. Estimates of aquifer hydraulic diffusivity are derived from estimates of the streamflow recession index and drainage density for 75 drainage basins; values range from 3,300 to 130,000 ft^2/d (feet squared per day). Basin-specific and site-specific estimates of transmissivity are computed from estimates of hydraulic diffusivity and specific-capacity test data, respectively. Basin-specific, or areal, estimates of transmissivity range from 22 to 1,300 ft^2/d, with a mean of 240 ft^2/d In general, areal transmissivity is highest for basins underlain by the Cumberland Plateau aquifer (mean value 480 ft^2/d) and lowest for basins underlain by the Central Basin aquifer (mean value 79 ft^2/d). Mean transmissivity values for the Highland Rim, Valley and Ridge, and Blue Ridge aquifer are 320,140, and 120 ft^2/d respectively. Site-specific estimates of transmissivity, computed from specific-capacity data from 118 test wells in Middle and East Tennessee range from 2 to 93,000 ft^2/d with a mean of 2,600 ft^2/d Mean transmissivity values for the Cumberland Plateau, Highland Rim, Central Basin, Valley and Ridge, and Blue Ridge aquifers are 2,800,1,200, 7,800, 390, and 65Oft Id, respectively.

Tennessee

Evaluation of Street Sweeping as a Stormwater-Quality-Management Tool in Three Residential Basins in Madison, Wisconsin

Recent technological improvements have increased the ability of street sweepers to remove sediment and other debris from street surfaces; the effect of these technological advancements on stormwater quality is largely unknown. The U.S. Geological Survey, in cooperation with the City of Madison and the Wisconsin Department of Natural Resources, evaluated three street-sweeper technologies from 2002 through 2006. Regenerative-air, vacuum-assist, and mechanical-broom street sweepers were operated on a frequency of once per week (high frequency) in separate residential basins in Madison, Wis., to measure each sweeper's ability to not only reduce street-dirt yield but also improve the quality of stormwater runoff. A second mechanical-broom sweeper operating on a frequency of once per month (low frequency) was also evaluated to measure reductions in street-dirt yield only. A paired-basin study design was used to compare street-dirt and stormwater-quality samples during a calibration (no sweeping) and a treatment period (weekly sweeping). The basis of this paired-basin approach is that the relation between paired street-dirt and stormwater-quality loads for the control and tests basins is constant until a major change is made at one of the basins. At that time, a new relation will develop. Changes in either street-dirt and/or stormwater quality as a result of street sweeping could then be quantified by use of statistical tests. Street-dirt samples collected weekly during the calibration period and twice per week during the treatment period, once before and once after sweeping, were dried and separated into seven particle-size fractions ranging from less than 63 micrometers to greater than 2 millimeters. Street-dirt yield evaluation was based on a computed mass per unit length of pounds per curb-mile. An analysis of covariance was used to measure the significance of the effect of street sweeping at the end of the treatment period and to quantify any reduction in street-dirt yield. Both the regenerative-air and vacuum-assist sweepers produced reductions in street-dirt yield at the 5-percent significance level. Street-dirt yield was reduced by an average of 76, 63, and 20 percent in the regenerative-air, vacuum-assist, and high-frequency broom basins, respectively. The low-frequency broom basin showed no significant reductions in street-dirt yield. Sand-size particles (greater than 63 micrometers) recorded the greatest overall reduction. Street-sweeper pickup efficiency was determined by computing the difference between weekly street-dirt yields before and after sweeping cleaning. The regenerative-air and vacuum-assist sweepers had similar pickup efficiencies of 25 and 30 percent, respectively. The mechanical broom sweeper operating at high frequency was considerably less efficient, removing an average of 5 percent of street-dirt yield. The effects of street sweeping on stormwater quality were evaluated by use of statistical tests to compare event mean concentrations and loads computed for individual storms at the control and test basins. Loads were computed by multiplying the event mean concentrations by storm-runoff volumes. Only ammonia-nitrogen for the test basin with the vacuum-assist sweeper showed significant load increases over the control basin, at the 10-percent significance level, of 63 percent. Difficulty in detecting significant changes in constituent stormwater-quality loads could be due, in part, to the large amount of variability in the data. Coefficients of variation for the majority of constituent loads were greater than 1, indicating substantial variability. The ability to detect changes in constituent stormwater-quality loads was likely hampered by an inadequate number of samples in the data set. However, sediment transport in the storm-sewer pipe, sediment washing onto the street from other source areas, winter sand application, and sampling challenges were additional sources of variability within each study ba

Scientific Investigations Report

Statistical analysis and mapping of water levels in the Biscayne aquifer, water conservation areas, and Everglades National Park, Miami-Dade County, Florida, 2000–2009

Statistical analyses and maps representing mean, high, and low water-level conditions in the surface water and groundwater of Miami-Dade County were made by the U.S. Geological Survey, in cooperation with the Miami-Dade County Department of Regulatory and Economic Resources, to help inform decisions necessary for urban planning and development. Sixteen maps were created that show contours of (1) the mean of daily water levels at each site during October and May for the 2000–2009 water years; (2) the 25th, 50th, and 75th percentiles of the daily water levels at each site during October and May and for all months during 2000–2009; and (3) the differences between mean October and May water levels, as well as the differences in the percentiles of water levels for all months, between 1990–1999 and 2000–2009. The 80th, 90th, and 96th percentiles of the annual maximums of daily groundwater levels during 1974–2009 (a 35-year period) were computed to provide an indication of unusually high groundwater-level conditions. These maps and statistics provide a generalized understanding of the variations of water levels in the aquifer, rather than a survey of concurrent water levels. Water-level measurements from 473 sites in Miami-Dade County and surrounding counties were analyzed to generate statistical analyses. The monitored water levels included surface-water levels in canals and wetland areas and groundwater levels in the Biscayne aquifer. Maps were created by importing site coordinates, summary water-level statistics, and completeness of record statistics into a geographic information system, and by interpolating between water levels at monitoring sites in the canals and water levels along the coastline. Raster surfaces were created from these data by using the triangular irregular network interpolation method. The raster surfaces were contoured by using geographic information system software. These contours were imprecise in some areas because the software could not fully evaluate the hydrology given available information; therefore, contours were manually modified where necessary. The ability to evaluate differences in water levels between 1990–1999 and 2000–2009 is limited in some areas because most of the monitoring sites did not have 80 percent complete records for one or both of these periods. The quality of the analyses was limited by (1) deficiencies in spatial coverage; (2) the combination of pre- and post-construction water levels in areas where canals, levees, retention basins, detention basins, or water-control structures were installed or removed; (3) an inability to address the potential effects of the vertical hydraulic head gradient on water levels in wells of different depths; and (4) an inability to correct for the differences between daily water-level statistics. Contours are dashed in areas where the locations of contours have been approximated because of the uncertainty caused by these limitations. Although the ability of the maps to depict differences in water levels between 1990–1999 and 2000–2009 was limited by missing data, results indicate that near the coast water levels were generally higher in May during 2000–2009 than during 1990–1999; and that inland water levels were generally lower during 2000–2009 than during 1990–1999. Generally, the 25th, 50th, and 75th percentiles of water levels from all months were also higher near the coast and lower inland during 2000–2009 than during 1990–1999. Mean October water levels during 2000–2009 were generally higher than during 1990–1999 in much of western Miami-Dade County, but were lower in a large part of eastern Miami-Dade County.

Florida

An analysis of potential water availability from the Atwood, Leesville, and Tappan Lakes in the Muskingum River Watershed, Ohio

This report presents the results of a study to assess potential water availability from the Atwood, Leesville, and Tappan Lakes, located within the Muskingum River Watershed, Ohio. The assessment was based on the criterion that water withdrawals should not appreciably affect maintenance of recreation-season pool levels in current use. To facilitate and simplify the assessment, it was assumed that historical lake operations were successful in maintaining seasonal pool levels, and that any discharges from lakes constituted either water that was discharged to prevent exceeding seasonal pool levels or discharges intended to meet minimum in-stream flow targets downstream from the lakes. It further was assumed that the volume of water discharged in excess of the minimum in-stream flow target is available for use without negatively impacting seasonal pool levels or downstream water uses and that all or part of it is subject to withdrawal. Historical daily outflow data for the lakes were used to determine the quantity of water that potentially could be withdrawn and the resulting quantity of water that would flow downstream (referred to as “flow-by”) on a daily basis as a function of all combinations of three hypothetical target minimum flow-by amounts (1, 2, and 3 times current minimum in-stream flow targets) and three pumping capacities (1, 2, and 3 million gallons per day). Using both U.S. Geological Survey streamgage data and lake-outflow data provided by the U.S. Army Corps of Engineers resulted in analytical periods ranging from 51 calendar years for the Atwood Lake to 73 calendar years for the Leesville and Tappan Lakes. The observed outflow time series and the computed time series of daily flow-by amounts and potential withdrawals were analyzed to compute and report order statistics (95th, 75th, 50th, 25th, 10th, and 5th percentiles) and means for the analytical period, in aggregate, and broken down by calendar month. In addition, surplus-water mass curve data were tabulated for each of the lakes. Monthly order statistics of computed withdrawals indicated that, for the three pumping capacities considered, increasing the target minimum flow-by amount tended to reduce the amount of water that can be withdrawn. The reduction was greatest in the lower percentiles of withdrawal; however, increasing the flow-by amount had no impact on potential withdrawals during high flow. In addition, for a given target minimum flow-by amount, increasing the pumping rate increased the total amount of water that could be withdrawn; however, that increase was less than a direct multiple of the increase in pumping rate for most flow statistics. Potential monthly withdrawals were observed to be more variable and more limited in some calendar months than others. Monthly order statistics and means of computed daily mean flow-by amounts indicated that flow-by amounts generally tended to be lowest during June–October and February. Increasing the target minimum flow-by amount for a given pumping rate resulted in some small increases in the magnitudes of the mean and 50th percentile and lower order statistics of computed mean flow-by, but had no effect on the magnitudes of the higher percentile statistics. Increasing the pumping rate for a given target minimum flow-by amount resulted in decreases in magnitudes of higher-percentile flow-by statistics by an amount equal to the flow equivalent of the increase in pumping rate; however, some lower percentile statistics remained unchanged.

Ohio

Hydrology and water quality in 13 watersheds in Gwinnett County, Georgia, 2001–15

The U.S. Geological Survey (USGS), in cooperation with Gwinnett County Department of Water Resources, established a Long-Term Trend Monitoring (LTTM) program in 1996. The LTTM program is a comprehensive, long-term, water-quantity and water-quality monitoring program designed to document and analyze the hydrologic and water-quality conditions of selected watersheds in Gwinnett County, Georgia. Water-quality monitoring initially began in six watersheds and currently [2016] includes 13 watersheds. As part of the LTTM program, streamflow, precipitation, water temperature, specific conductance, and turbidity were measured every 15 minutes for water years 2001–15 at 12 of the 13 watershed monitoring stations and for water years 2010–15 at the other watershed. In addition, discrete water-quality samples were collected seasonally from May through October (summer) and November through April (winter), including one base-flow and three stormflow event composite samples, during the study period. Samples were analyzed for nutrients (nitrogen and phosphorus), total organic carbon, trace elements (total lead and total zinc), total dissolved solids, and total suspended sediment (total suspended solids and suspended-sediment concentrations). The sampling scheme was designed to identify variations in water quality both hydrologically and seasonally. The 13 watersheds were characterized for basin slope, population density, land use for 2012, and the percentage of impervious area from 2000 to 2014. Several droughts occurred during the study period—water years 2002, 2007–08, and 2011–12. Watersheds with the highest percentage of impervious areas had the highest runoff ratios, which is the portion of precipitation that occurs as runoff. Watershed base-flow indexes, the ratio of base-flow runoff to total runoff, were inversely correlated with watershed impervious area. Flood-frequency estimates were computed for 13 streamgages in the study area that have 10 or more years of annual peak flow data through water year 2015, using the expected moments algorithm to fit a Pearson Type III distribution to logarithms of annual peak flows. Kendall’s tau nonparametric test was used to determine the statistical significance of trends in the annual peak flows, with none of the 13 streamgages exhibiting significant trends. A comparison of base-flow and stormflow water-quality samples indicates that turbidity and concentrations of total ammonia plus organic nitrogen, total nitrogen, total phosphorus, total organic carbon, total lead, total zinc, total suspended solids, and suspended-sediment concentrations increased with increasing discharge at all watersheds. Specific conductance decreased during stormflow at all watersheds, and total dissolved solids concentrations decreased during stormflow at a few of the watersheds. Total suspended solids and suspended-sediment concentrations typically were two orders of magnitude higher in stormflow samples, turbidities were about 1.5 orders of magnitude higher, total phosphorus and total zinc were about one order of magnitude higher, and total ammonia plus organic nitrogen, total nitrogen, total organic carbon, and total lead were about twofold higher than in base-flow samples. Seasonality and long-term trends were identified for the period water years 2001–15 for 10 constituents—total nitrogen, total nitrate plus nitrite, total phosphorus, dissolved phosphorus, total organic carbon, total suspended solids, suspended-sediment concentration, total lead, total zinc, and total dissolved solids. Seasonal patterns were present in most watersheds for all constituents except total dissolved solids, and the watersheds had fairly similar patterns of higher concentrations in the summer and lower concentrations in the winter. A linear long-term trend analysis of residual concentrations from the flow-only load estimation model (without time-trend terms) identified significant trends in 67 of the 130 constituent-watershed combinations. Seventy percent of the significant trends were negative. Total organic carbon and total dissolved solids had predominantly positive trends. Total phosphorus, total suspended solids, suspended-sediment concentration, total lead, and total zinc had only negative trends. The other three constituents exhibited fewer trends, both positive and negative. Streamwater loads were estimated annually for the 13-year period water years 2003–15 for the same 10 constituents in the trend analysis. Loads were estimated using a regression-model-based approach developed by the USGS for the Gwinnett County LTTM program that accommodates the use of storm-event composited samples. Concentrations were modeled as a function of discharge, base flow, time, season, and turbidity to improve model predictions and reduce errors in load estimates. Total suspended solids annual loads have been identified in Gwinnett County’s Watershed Protection Plan for target performance criterion. Although the amount of annual runoff was the primary factor in variations in annual loads, climatic conditions (classified as dry, average, or wet) affected annual loads beyond what was attributed to climatic-related variations in annual runoff. Significant negative trends in loads were estimated for the combined area of the watersheds for all constituents except dissolved phosphorus, total organic carbon, and total dissolved solids. The trend analysis indicated that total suspended solids and suspended-sediment concentration loads in the study area were decreasing by 57,000 and 87,000 pounds per day per year, respectively. Variations in constituent yields between watersheds appeared to be related to various watershed characteristics. Suspended sediment (as either total suspended solids or suspended-sediment concentrations), along with constituents transported predominately in solid phase (total phosphorus, total organic carbon, total lead, and total zinc), and total dissolved solids typically had higher yields from watersheds that had high percentages of impervious areas or high basin slope. High total nitrogen yields were also associated with watersheds with high percentages of impervious areas. Low total nitrogen, total suspended solids, total lead, and total zinc yields appeared to be associated with watersheds that had a low percentage of high-density development.

Georgia

Analysis of the magnitude and frequency of the 4-day annual low flow and regression equations for estimating the 4-day, 3-year low-flow frequency at ungaged sites on unregulated streams in New Mexico

Two regression equations were developed for estimating the 4-day, 3-year (4Q3) low-flow frequency at ungaged sites on unregulated streams in New Mexico. The first, a statewide equation for estimating the 4Q3 low-flow frequency from drainage area and average basin mean winter precipitation, was developed from the data for 50 streamflow-gaging stations that had non-zero 4Q3 low-flow frequency. The 4Q3 low-flow frequency for the 50 gaging stations ranged from 0.08 to 18.7 cubic feet per second. For this statewide equation, the average standard error of estimate was 126 percent and the coefficient of determination was 0.48. The second, an equation for estimating the 4Q3 low-flowfrequency in mountainous regions from drainage area, average basin mean winter precipitation, and average basin slope, was developed from the data for 40 gaging stations located above 7,500 feet in elevation. For this regression equation, the average standard error of estimate was 94 percent and the coefficient of determination was 0.66. A U.S. Geological Survey computer-program interface for a geographical information system (GIS), called the GISWeasel,was used to determine basin and climatic characteristics for 84 gaging stations that were not affected by regulation. Mean monthly precipitation estimates from 1961 to 1990 were used in the GIS Weasel to compute the climatic characteristics of average basin winter precipitation and annual mean precipitation. The U.S. Geological Survey National Elevation Dataset, which currently consists of the 7.5-minute, 30-meter digital elevation model for each State, was used in the GISWeasel to compute the basin characteristics of drainage area, average basin slope, average basin elevation, and average basin aspect. Basin and climatic characteristics that were statistically significant in the regression equation with the 4Q3 lowflow frequency were drainage area, which ranged from 1.62 to 5,900 square miles; average basin mean winter precipitation, which ranged from 3.89 to 19.42 inches; and average basin slope, which ranged from 0.166 to 0.517 percent.

Water-Resources Investigations Report

Depth-Duration Frequency of Precipitation for Oklahoma

A regional frequency analysis was conducted to estimate the depth-duration frequency of precipitation for 12 durations in Oklahoma (15, 30, and 60 minutes; 1, 2, 3, 6, 12, and 24 hours; and 1, 3, and 7 days). Seven selected frequencies, expressed as recurrence intervals, were investigated (2, 5, 10, 25, 50, 100, and 500 years). L-moment statistics were used to summarize depth-duration data and to determine the appropriate statistical distributions. Three different rain-gage networks provided the data (15minute, 1-hour, and 1-day). The 60-minute, and 1-hour; and the 24-hour, and 1-day durations were analyzed separately. Data were used from rain-gage stations with at least 10-years of record and within Oklahoma or about 50 kilometers into bordering states. Precipitation annual maxima (depths) were determined from the data for 110 15-minute, 141 hourly, and 413 daily stations. The L-moment statistics for depths for all durations were calculated for each station using unbiased L-mo-ment estimators for the mean, L-scale, L-coefficient of variation, L-skew, and L-kur-tosis. The relation between L-skew and L-kurtosis (L-moment ratio diagram) and goodness-of-fit measures were used to select the frequency distributions. The three-parameter generalized logistic distribution was selected to model the frequencies of 15-, 30-, and 60-minute annual maxima; and the three-parameter generalized extreme-value distribution was selected to model the frequencies of 1-hour to 7-day annual maxima. The mean for each station and duration was corrected for the bias associated with fixed interval recording of precipitation amounts. The L-scale and spatially averaged L-skew statistics were used to compute the location, scale, and shape parameters of the selected distribution for each station and duration. The three parameters were used to calculate the depth-duration-frequency relations for each station. The precipitation depths for selected frequencies were contoured from weighted depth surfaces to produce maps from which the precipitation depth-duration-frequency curve for selected storm durations can be determined for any site in Oklahoma.

Water-Resources Investigations Report

Effects on water quality due to flood-water detention by Barker and Addicks Reservoirs, Houston, Texas

The Barker and Addicks Reservoirs, located about 16 miles west of Houston, Texas, provide flood-detention storage for storm runoff. Of interest are the water-quality characteristics in the study area and changes in water quality during detention. Study area sampling sites were selected upstream along Buffalo Bayou for Barker Reservoir and on Bear Creek and Langham Creek for Addicks Reservoir, within the reservoirs, near the reservoir outflows, and below the confluence of each reservoir outflow at the streamflow station Buffalo Bayou near Addicks. Flow data were available at all sites except in the reservoirs. Analyses of samples collected during both low flow and storm runoff show that, in general, the waters of the study area were low in mineralization, but the aesthetics of the water was a problem. The inorganic constituents, trace metals, and pesticides rarely exceeded maximum contaminant levels recommended by the U.S. Environmental Protection Agency for public supply using 1976 and 1977 criteria for primary and secondary standards. All species of nutrients, except ammonia nitrogen and phosphorus, almost always were below the recommended maximum contaminant levels. Phosphorus almost always exceeded these levels. Aesthetic problems are evident. Large values of suspended solids, turbidity, and color were common. Small dissolved-oxygen values commonly occurred in the reservoirs. Possible bacterial problems are indicated because coliformbacteria densities exceeded recommended levels in about 25 percent of the samples. The effects of the reservoirs on the water-quality characteristics of storm runoff were analyzed using three approaches. The first approach was a comparison of the discharge-weighted average values of nine selected constituents at each streamflow-gaging station during four storms. Reservoir effects on the quality of runoff detained 1 to 4 days in the two reservoirs were inconsistent. However, the reservoirs consistently had an effect on the water quality of runoff that was detained the longest (more than 8 days). Biochemical oxygen demand, suspended solids, turbidity, color, total nitrogen, and total organic carbon discharge-weighted average values were consistently smaller after flowing through the reservoirs. Dissolved solids and total phosphorus values were consistently larger after flowing through the reservoirs. The second approach was an analysis of the means of the discharge-weighted average values computed for the four hydrologic events using the Student t-test. Statistical results Indicate that reservoir detention significantly reduced suspended solids (the mean decreased from 178 milligrams per liter at the Inflows to 105 milligrams per liter at the outflows) and turbidity (the mean decreased from 119 nephelometrlc turbidity units at the inflows to 66 nephelometric turbidity units at the outflows). The third approach was a comparison at each site of the mean, maximum, and minimum values computed for seven constituents that did not correlate with discharge. These constituents or properties of water were temperature, pH, dissolved oxygen, dissolved oxygen percent saturation, total-coliform bacteria, fecal-conform bacteria, and fecal-streptococci bacteria. The only consistent water-quality changes observed were with the three bacteria groups, which were decreased by flood-water detention.

Texas

Linkage of the Soil and Water Assessment Tool and the Texas Water Availability Model to simulate the effects of brush management on monthly storage of Canyon Lake, south-central Texas, 1995-2010

The U.S. Geological Survey (USGS), in cooperation with the Texas State Soil and Water Conservation Board, developed and applied an approach to create a linkage between the published upper Guadalupe River Soil Water Assessment Tool (SWAT) brush-management (ashe juniper [ Juniperus ashei ]) model and the full authorization version Guadalupe River Water Availability Model (WAM). The SWAT model was published by the USGS, and the Guadalupe River WAM is available from the Texas Commission on Environmental Quality. The upper Guadalupe River watershed is a substantial component of the Guadalupe River WAM. This report serves in part as documentation of a proof of concept on the feasibility of linking these two water-resources planning models for the purpose of simulating possible increases in water storage in Canyon Lake as a result of different brush-management scenarios. The SWAT-WAM linkage for the upper Guadalupe River is documented with a principal objective to evaluate the distributional characteristics of the monthly water storage of Canyon Lake during selected drought conditions. Focus is on the relative evaluation of select scenarios of large-scale or “extensive” brush management within the upper Guadalupe River watershed. There are six SWAT simulations for the upper Guadalupe River watershed that include a baseline (0-percent management of treatable ashe juniper, the baseline scenario from a previous study in which no percentage of ashe juniper is numerically replaced with grassland) along with five scenarios (extensions of SWAT simulations from a previous study) of 20-, 40-, 60-, 80-, and 100-percent random (numerical) replacement of treatable ashe juniper with grasslands throughout the upper Guadalupe River watershed in south-central Texas. SWAT is a process-based, semidistributed, water-balance model designed to predict the effects of landscape management decisions on water yields. A watershed is subdivided into subbasins, and each subbasin is associated with a single reach on the stream network. In general a WAM, such as the Guadalupe River WAM, provides analysis of generalized water rights in a river and reservoir framework. A WAM accommodates hydrology and water usage through several input files containing water rights, watershed parameters, and naturalized streamflow time series. A WAM is generalized for application to rivers and reservoir systems, and input datasets are uniquely developed for a river basin of concern. The extractions of SWAT output for the five extensive brush-management and baseline scenarios were offset by –21 years and, in general, the results were then mapped to the WAM input-flow file. The offset of –21 years was chosen arbitrarily for technical reasons and means that the period of monthly record 1995–2010 of the upper Guadalupe River SWAT became the synthetic period of monthly record 1974–89, hereinafter 1974–89 (synthetic) period, of the Guadalupe River WAM. The relative (between scenario to baseline) effects of extensive brush-management scenarios by using the SWAT-WAM linkage were evaluated, and two critical intermediate results were total inflow to Canyon Lake from 1995 to 2010 and the monthly storage of Canyon Lake from 1974 to 1989 (synthetic). The first quartile or lower 25th percentile of monthly storage of Canyon Lake for the baseline scenario is 381,000 acre-feet (acre-ft) for the hereinafter 1974–89 (synthetic) period. This lower quartile was chosen for analysis for two critical purposes. First, Canyon Lake is managed with a conservation pool of about 386,200 acre-ft capacity (as recognized by the WAM) and is at or near conservation capacity about 50 percent or more of the time; further, there is intrinsic data censoring that occurs for the monthly storage distribution because Canyon Lake is at or near conservation pool elevation the majority of the time. This intrinsic censoring has the effect of creating a bounded distribution with a left or low-volume tail. Statistical assessment of the brush-management scenarios beginning with the 381,000 acre-ft censoring threshold provides readily interpretable results. Second, the quantification of brush management during periods lacking abundant rainfall, which were defined in this study as months for which Canyon Lake storage was below the 25th percentile for the simulation period, are of substantial interest to water-resource managers and stakeholders in the context of water-supply enhancement. A statistical assessment of the SWAT-WAM linkage for the low-volume tail of the distribution of monthly storage of Canyon Lake is the focus of analysis and interpretation. Drought periods for the analysis are defined as the months (consecutive or not) during which Canyon Lake is below the 25th percentile of storage (381,000 acre-ft) for the baseline scenario. Such months are referred to as being within the “Drought Quartile.” The Drought Quartile is a conceptual and heuristically determined waypoint for the analysis and is not related to any administrative definition of drought by stakeholders or policy makers. The five scenarios and the baseline scenario simulated in the upper Guadalupe River SWAT were all passed through the Guadalupe River WAM by the SWAT-WAM linkage described in this report. A comparison of the mean increase per month in reservoir storage for Canyon Lake conditioned for the Drought Quartile was made. For each of the five brush-management and baseline scenarios, the months with storage below 381,000 acre-ft were extracted. The mean monthly storages during the Drought Quartile were computed for each of the five scenarios and the baseline scenario. The mean of the baseline scenario was 376,458 acre-ft and subsequently was subtracted from the mean monthly storage during the Drought Quartile for each of the five scenarios. The mean monthly offset storages of Canyon Lake during the Drought Quartile were 110 acre-ft (20 percent); 448 acre-ft (40 percent); 754 acre-ft (60 percent); 1,080 acre-ft (80 percent); and 1,090 acre-ft (100 percent). A particular mean was interpreted as follows: the value of 754 acre-ft for the 60-percent brush-management scenario implies that, on average, this scenario indicates an additional 754 acre-ft per month of storage in Canyon Lake relative to the baseline during the Drought Quartile. All of the five scenarios resulted in an increase on average to water supply relative to the baseline scenario during the Drought Quartile through the SWAT-WAM linkage.

Texas

Techniques for estimating the magnitude and frequency of floods in rural basins of South Carolina, 1999

Data from 167 streamflow-gaging stations in or near South Carolina with 10 or more years of record through September 30, 1999, were used to develop two methods for estimating the magnitude and frequency of floods in South Carolina for rural ungaged basins that are not significantly affected by regulation. Floodfrequency estimates for 54 gaged sites in South Carolina were computed by fitting the water-year peak flows for each site to a log-Pearson Type III distribution. As part of the computation of flood-frequency estimates for gaged sites, new values for generalized skew coefficients were developed. Flood-frequency analyses also were made for gaging stations that drain basins from more than one physiographic province. The U.S. Geological Survey, in cooperation with the South Carolina Department of Transportation, updated these data from previous flood-frequency reports to aid officials who are active in floodplain management as well as those who design bridges, culverts, and levees, or other structures near streams where flooding is likely to occur. Regional regression analysis, using generalized least squares regression, was used to develop a set of predictive equations that can be used to estimate the 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence- interval flows for rural ungaged basins in the Blue Ridge, Piedmont, upper Coastal Plain, and lower Coastal Plain physiographic provinces of South Carolina. The predictive equations are all functions of drainage area. Average errors of prediction for these regression equations ranged from -16 to 19 percent for the 2-year recurrence-interval flow in the upper Coastal Plain to -34 to 52 percent for the 500-year recurrenceinterval flow in the lower Coastal Plain. A region-of-influence method also was developed that interactively estimates recurrence- interval flows for rural ungaged basins in the Blue Ridge of South Carolina. The region-of-influence method uses regression techniques to develop a unique relation between flow and basin characteristics for an individual watershed. This, then, can be used to estimate flows at ungaged sites. Because the computations required for this method are somewhat complex, a computer application was developed that performs the computations and compares the predictive errors for this method. The computer application includes the option of using the region-of-influence method, or the generalized least squares regression equations from this report to compute estimated flows and errors of prediction specific to each ungaged site. From a comparison of predictive errors using the region-of-influence method with those computed using the regional regression method, the region-of-influence method performed systematically better only in the Blue Ridge and is, therefore, not recommended for use in the other physiographic provinces. Peak-flow data for the South Carolina stations used in the regionalization study are provided in appendix A, which contains gaging station information, log- Pearson Type III statistics, information on stage-flow relations, and water-year peak stages and flows. For informational purposes, water-year peak-flow data for stations on regulated streams in South Carolina also are provided in appendix D. Other information pertaining to the regulated streams is provided in the text of the report.

South Carolina

An analysis of potential water availability from the Charles Mill, Clendening, Piedmont, Pleasant Hill, Senecaville, and Wills Creek Lakes in the Muskingum River Watershed, Ohio

This report presents the results of a study to assess potential water availability from the Charles Mill, Clendening, Piedmont, Pleasant Hill, Senecaville, and Wills Creek Lakes, located within the Muskingum River Watershed, Ohio. The assessment was based on the criterion that water withdrawals should not appreciably affect maintenance of recreation-season pool levels in current use. To facilitate and simplify the assessment, it was assumed that historical lake operations were successful in maintaining seasonal pool levels, and that any discharges from lakes constituted either water that was discharged to prevent exceeding seasonal pool levels or discharges intended to meet minimum in-stream flow targets downstream from the lakes. It further was assumed that the volume of water discharged in excess of the minimum in-stream flow target is available for use without negatively impacting seasonal pool levels or downstream water uses and that all or part of it is subject to withdrawal. Historical daily outflow data for the lakes were used to determine the quantity of water that potentially could be withdrawn and the resulting quantity of water that would flow downstream (referred to as “flow-by”) on a daily basis as a function of all combinations of three hypothetical target minimum flow-by amounts (1, 2, and 3 times current minimum in-stream flow targets) and three pumping capacities (1, 2, and 3 million gallons per day). Using both U.S. Geological Survey streamgage data (where available) and lake-outflow data provided by the U.S. Army Corps of Engineers resulted in analytical periods ranging from 51 calendar years for Charles Mill, Clendening, and Piedmont Lakes to 74 calendar years for Pleasant Hill, Senecaville, and Wills Creek Lakes. The observed outflow time series and the computed time series of daily flow-by amounts and potential withdrawals were analyzed to compute and report order statistics (95th, 75th, 50th, 25th, 10th, and 5th percentiles) and means for the analytical period, in aggregate, and broken down by calendar month. In addition, surplus-water mass curve data were tabulated for each of the lakes. Monthly order statistics of computed withdrawals indicated that, for the three pumping capacities considered, increasing the target minimum flow-by amount tended to reduce the amount of water that can be withdrawn. The reduction was greatest in the lower percentiles of withdrawal; however, increasing the flow-by amount had no impact on potential withdrawals during high flow. In addition, for a given target minimum flow-by amount, increasing the pumping rate typically increased the total amount of water that could be withdrawn; however, that increase was less than a direct multiple of the increase in pumping rate for most flow statistics. Potential monthly withdrawals were observed to be more variable and more limited in some calendar months than others. Monthly order statistics and means of computed daily mean flow-by amounts indicated that flow-by amounts generally tended to be lowest during June–October. Increasing the target minimum flow-by amount for a given pumping rate resulted in some small increases in the magnitudes of the mean and 50th percentile and lower order statistics of computed mean flow-by, but had no effect on the magnitudes of the higher percentile statistics. Increasing the pumping rate for a given target minimum flow-by amount resulted in decreases in magnitudes of higher-percentile flow-by statistics by an amount equal to the flow equivalent of the increase in pumping rate; however, some lower percentile statistics remained unchanged.

Ohio

Two models for evaluating landslide hazards

Two alternative procedures for estimating landslide hazards were evaluated using data on topographic digital elevation models (DEMs) and bedrock lithologies in an area adjacent to the Missouri River in Atchison County, Kansas, USA. The two procedures are based on the likelihood ratio model but utilize different assumptions. The empirical likelihood ratio model is based on non-parametric empirical univariate frequency distribution functions under an assumption of conditional independence while the multivariate logistic discriminant model assumes that likelihood ratios can be expressed in terms of logistic functions. The relative hazards of occurrence of landslides were estimated by an empirical likelihood ratio model and by multivariate logistic discriminant analysis. Predictor variables consisted of grids containing topographic elevations, slope angles, and slope aspects calculated from a 30-m DEM. An integer grid of coded bedrock lithologies taken from digitized geologic maps was also used as a predictor variable. Both statistical models yield relative estimates in the form of the proportion of total map area predicted to already contain or to be the site of future landslides. The stabilities of estimates were checked by cross-validation of results from random subsamples, using each of the two procedures. Cell-by-cell comparisons of hazard maps made by the two models show that the two sets of estimates are virtually identical. This suggests that the empirical likelihood ratio and the logistic discriminant analysis models are robust with respect to the conditional independent assumption and the logistic function assumption, respectively, and that either model can be used successfully to evaluate landslide hazards. ?? 2006.

Computers & Geosciences

Seasonal and long-term clarity trend assessment of Lake Tahoe, California–Nevada

The clarity of Lake Tahoe, observed using a Secchi disk on a regular basis since the late 1960s, continues to be a sentinel metric of lake health. Water clarity is influenced by physical and biological processes and has declined in the five decades of monitoring, revealing differences between summer (June–September) and winter (December–March). This document summarizes key findings of a study of Lake Tahoe water clarity, including long-term variability and the relative importance of several influencing variables and processes. This study, prepared in cooperation with the Nevada Division of Environmental Protection, focused on (1) an apparent divergence in clarity trends between summer and winter periods, (2) observed changes in in-lake physical and ecological variables that may influence or control seasonal and annual clarity trends, and (3) five research hypotheses regarding lake clarity that were developed by Lake Tahoe management agencies. Previously collected data were used to complete this study. Trend analysis confirmed that winter clarity stabilized (that is, there is no longer a statistically significant trend up or down) during the last 20 years. Evaluation of clarity for selected months in the 50-year Secchi disk clarity dataset showed that only two summer months, July and August, had statistically significant decreases in clarity from 2000–19. Different subsets of available data were analyzed to reveal the presence or absences of trends for each season, decade, and month. Five hypotheses related to lake clarity were part of the study described by this report. Hypothesis 1 stated that clarity is controlled predominantly by the distribution and volumetric density of fine particles in suspension. This hypothesis was studied using available data describing in-lake fine (0–20 micrometers) particles from 2008–19. Water clarity was negatively correlated with in-lake particle abundance, with particles in the 1.0-4.6 μm range having the greatest effect, consistent with light-scattering theory. Estimated abundances of diatoms of the genus Cyclotella also were found to be negatively correlated with clarity. Data limitations precluded a complete investigation of hypothesis 2, which stated that the observed improvements in winter water clarity are a response to decreasing fine suspended-sediment concentrations in the lake resulting from load reductions from upland sources in and near urbanized areas. Data describing fine-sediment loading from urban areas to the lake were only available since 2014, and only once or twice per month. A slight, statistically significant, negative correlation was identified between urban fine-particle loading and monthly lake clarity with a 4-month lag. Particle abundance in monitored streams is highly correlated with simultaneous particle abundance in the lake. Hypothesis 3 stated that changing hydrodynamic conditions in the lake are increasing thermal stability and resistance to mixing. Trend analyses performed on stability index and buoyancy frequency time series computed from long-term observations of lake temperatures support the hypothesis that hydrodynamic conditions have evolved since 1969 to increase the lake’s resistance to mixing. The date of maximum mixing in winter has become progressively earlier in the year. Lake density stratification, defined using the stability index, is commencing earlier in the year and extending a month longer than in the early years of the monitoring program. Hypothesis 4 stated that the trend of decreasing summer clarity is a result of earlier, prolonged, and more intense stratification. Statistically significant correlations were found between summer clarity and (1) date of onset of stratification, (2) duration of stratification, and (3) buoyancy frequency. Hypothesis 5 stated that ecological (food web) interactions are causing changes in the trends of seasonal or annual clarity; data supporting hypothesis 5 were limited to examples from other systems and to intermittent monitoring of Lake Tahoe and Emerald Bay. The resulting narrative assessment was motivated by a 6-year study of Mysis shrimp disappearance and return in Emerald Bay. The available data and a large body of published literature are consistent with the inference that Mysis shrimp-induced food web changes are causing changes in the trends of seasonal or annual clarity. This food-web study focused on the relations between introduced Mysis shrimp, the native cladocerans ( Daphnia and Bosmina ) that were largely eliminated following Mysis introduction, and the effect on fine particles within the lake. The records of Mysis and other zooplankton data for Lake Tahoe are episodic and have large gaps. Consequently, statistical analyses could not be conducted to compare zooplankton data with other variables. The long-term record, however, indicates that the key effect was a change to the phytoplankton assemblage, where larger diatoms disappeared, likely due to Mysis grazing, only to be replaced by Cyclotella that are an order-of-magnitude smaller and have increased the abundance and volumetric density of total fine particles in suspension (biotic and abiotic).

California

Evaluating external nutrient and suspended-sediment loads to Upper Klamath Lake, Oregon, using surrogate regressions with real-time turbidity and acoustic backscatter data

Executive Summary Suspended-sediment and total phosphorus loads were computed for two sites in the Upper Klamath Basin on the Wood and Williamson Rivers, the two main tributaries to Upper Klamath Lake. High temporal resolution turbidity and acoustic backscatter data were used to develop surrogate regression models to compute instantaneous concentrations and loads on these rivers. Regression models for the Williamson River site showed strong correlations of turbidity with total phosphorus and suspended-sediment concentrations (adjusted coefficients of determination [Adj R 2 ]=0.73 and 0.95, respectively). Regression models for the Wood River site had relatively poor, although statistically significant, relations of turbidity with total phosphorus, and turbidity and acoustic backscatter with suspended sediment concentration, with high prediction uncertainty. Total phosphorus loads for the partial 2014 water year (excluding October and November 2013) were 39 and 28 metric tons for the Williamson and Wood Rivers, respectively. These values are within the low range of phosphorus loads computed for these rivers from prior studies using water-quality data collected by the Klamath Tribes. The 2014 partial year total phosphorus loads on the Williamson and Wood Rivers are assumed to be biased low because of the absence of data from the first 2 months of water year 2014, and the drought conditions that were prevalent during that water year. Therefore, total phosphorus and suspended-sediment loads in this report should be considered as representative of a low-water year for the two study sites. Comparing loads from the Williamson and Wood River monitoring sites for November 2013–September 2014 shows that the Williamson and Sprague Rivers combined, as measured at the Williamson River site, contributed substantially more suspended sediment to Upper Klamath Lake than the Wood River, with 4,360 and 1,450 metric tons measured, respectively. Surrogate techniques have proven useful at the two study sites, particularly in using turbidity to compute suspended-sediment concentrations in the Williamson River. This proof-of-concept effort for computing total phosphorus concentrations using turbidity at the Williamson and Wood River sites also has shown that with additional samples over a wide range of flow regimes, high-temporal-resolution total phosphorus loads can be estimated on a daily, monthly, and annual basis, along with uncertainties for total phosphorus and suspended-sediment concentrations computed using regression models. Sediment-corrected backscatter at the Wood River has potential for estimating suspended-sediment loads from the Wood River Valley as well, with additional analysis of the variable streamflow measured at that site. Suspended-sediment and total phosphorus loads with a high level of temporal resolution will be useful to water managers, restoration practitioners, and scientists in the Upper Klamath Basin working toward the common goal of decreasing nutrient and sediment loads in Upper Klamath Lake.

Oregon

Estimates of median flows for streams on the Kansas surface water register

The Kansas State Legislature, by enacting Kansas Statute KSA 82a-2001 et. seq., mandated the criteria for determining which Kansas stream segments would be subject to classification by the State. One criterion for the selection as a classified stream segment is based on the statistic of median flow being equal to or greater than 1 cubic foot per second. As specified by KSA 82a-2001 et. seq., median flows were determined from U.S. Geological Survey streamflow-gaging-station data by using the most-recent 10-years of gaged data (KSA) for each streamflow-gaging station. Median flows also were determined by using gaged data from the entire period of record (all-available hydrology, AAH). Least-squares multiple regression techniques were used, along with Tobit analyses, to develop equations for estimating median flows for uncontrolled stream segments. The drainage area of the uncontrolled gaging stations used in the regression analyses ranged from 2.06 to 12,004 square miles. A logarithmic transformation of the data was needed to develop the best linear relation for computing median flows. In the regression analyses, the significant climatic and basin characteristics, in order of importance, were drainage area, mean annual precipitation, mean basin permeability, and mean basin slope. Tobit analyses of KSA data yielded a root mean square error of 0.285 logarithmic units, and the best equations using Tobit analyses of AAH data had a root mean square error of 0.247 logarithmic units. These equations and an interpolation procedure were used to compute median flows for the uncontrolled stream segments on the Kansas Surface Water Register. Measured median flows from gaging stations were incorporated into the regression-estimated median flows along the stream segments where available. The segments that were uncontrolled were interpolated using gaged data weighted according to the drainage area and the bias between the regression-estimated and gaged flow information. On controlled reaches of Kansas streams, the median flow information was interpolated between gaging stations using only gaged data weighted by drainage area. Of the 2,232 total stream segments on the Kansas Surface Water Register, 30 percent of the segments had an estimated median streamflow of less than 1 cubic foot per second when the KSA analysis was used. When the AAH analysis was used, 40 percent of the segments had an estimated median streamflow of less than 1 cubic foot per second.

Kansas

Characterization and data-gap analysis of surface-water quality data in the Piceance study area, western Colorado, 1959–2009

The U.S. Geological Survey, in cooperation with Federal, State, county, and industry partners, developed a Web-accessible common data repository to provide access to historical and current (as of August 2009) water-quality information (available on the Internet at http://rmgsc.cr.usgs.gov/cwqdr/Piceance/index.shtml ). Surface-water-quality data from public and private sources were compiled for the period 1931 to 2009 and loaded into the common data repository for the Piceance Basin. A subset of surface-water-quality data for 1959 to 2009 from the repository were compiled, reviewed, and checked for quality assurance for this report. This report contains data summaries, comparisons to water-quality standards, trend analyses, a generalized spatial analysis, and a data-gap analysis for select water-quality properties and constituents. Summary statistics and a comparison to standards were provided for 347 sites for 33 constituents including field properties, nutrients, major ions, trace elements, suspended sediment, Escherichia coli, and BTEX (benzene, toluene, ethylbenzene, and xylene). When sufficient data were available, trends over time were analyzed and loads were calculated for those sites where there were also continuous streamflow data. The majority of sites had information on field properties. Water temperature data was available for 316 sites where data were collected between 1959 and 2009. The only trend that was detected in temperature was an upward trend at the Gunnison River near Grand Junction, Colorado. There were 326 values out of a total of 32,006 values in the study area that exceeded the aquatic-life standard for daily maximum water temperature. For the entire study area, 196 sites had dissolved-oxygen data collected between 1970 and 2009, and median dissolved-oxygen concentrations ranged from 6.8 to11.2 milligrams per liter (mg/L). There were 185 concentrations that exceeded the dissolved oxygen aquatic-life standard out of a total of 11,248 values. The pH data were available for 276 sites, and median pH values ranged from 7.5 to 9.0. There were 241 values that exceeded the high pH standard and 13 values that were less than the low pH standard of the 16,790 values in the study area. Nutrients within the study area were not well represented in each basin and were often not being sampled currently. For the entire study area, 62 sites had nitrate data collected between 1958 and 2009, and median nitrate concentrations ranged from less than detection to 3.72 mg/L as nitrogen. The maximum contaminant level for domestic water supply for nitrate is 10 mg/L and was exceeded once in 3,736 samples. Total phosphorus was collected at 113 sites between 1974 and 2009, and median total phosphorus concentrations ranged from less than detection to 5.04 mg/L. The U.S. Environmental Protection Agency recommendation for phosphorus is less than 0.1 mg/L, and 1,469 of 4,842 samples exceeded this recommended standard. An upward trend in both nitrate and total phosphorus was detected in the White River above Coal Creek near Meeker, Colo. Standards for major ions exist only for chloride and sulfate. For the entire study area, 118 sites had both chloride and sulfate concentration data collected between 1958 and 2009. Median chloride concentrations ranged from 0.085 mg/L to 280 mg/L. Median sulfate concentrations ranged from 4.57 mg/L to 15,000 mg/L. Both chloride and sulfate domestic water-supply standards are 250 mg/L. There were 120 chloride concentrations and 1,111 sulfate concentration samples that exceeded these standards. A downward trend in dissolved solids was detected at the Colorado River near the Colorado-Utah state border and could be a result of salinity control work near Grand Junction, Colo. Trace elements were relatively well represented both temporally and spatially in the study area though the number of trace element samples per site was not typically enough to compute trends or loads except for selenium. There were 127 sites that had dissolved iron concentration data collected between 1961 and 2009, and median iron concentrations ranged from less than detection to 1,100 micrograms per liter (µg/L). The 30-day drinking-water standard for iron is 300 µg/L, and 203 samples exceeded the standard. Selenium was the best represented trace element with selenium concentration data collected at 197 sites between 1973 and 2009, and median selenium concentrations range from less than detection to 181 µg/L. The chronic standard of 4.6 µg/L for selenium concentrations was exceeded in 899 samples, and the acute aquatic-life standard of 18.4 µg/ for selenium was exceeded in 629 samples. High concentrations of selenium are of concern in the Lower Gunnison River Basin because of the combination of geologic formations and land use. There were significant downward trends in selenium at both main-stem sites on the Gunnison River at Delta, Colo., and the Gunnison River near Grand Junction, Colo. High selenium concentrations correlate with high salinity concentrations; thus, when salinity control efforts are conducted in selenium-rich areas in the Lower Gunnison River Basin, both salinity and selenium have the potential to decrease. Spatial, temporal, and analytical data gaps were identified in the study area. The spatial coverage of sampling sites could be expanded in the White River Basin by adding more tributary sites. No water-quality data exist for tributary streams in the area north of Rangely, Colo., where extensive energy development has occurred in a complex geologic setting. Douglas Creek has a drainage area of 425 square miles and has limited historic water-quality and water-quantity data. Limited data were available for field properties, major ions, nutrients, and trace elements on the main stem of the Colorado River between Glenwood Springs and Cameo, Colo. Nutrient data were minimally collected upstream from Colorado River at the Colorado-Utah state border and on the Gunnison River (major tributary in the reach). Approximately 30 percent of the samples for total phosphorus in the Lower Gunnison River Basin exceeded the recommended standard, yet there were insufficient data to do trends analysis in the Lower Gunnison River Basin except at the Gunnison near Grand Junction site. There is limited trace element data except for selenium in the Lower Gunnison River Basin. Additional sampling is necessary to understand the occurrence, concentrations, and loads of these constituents.

Colorado

Evaluation of statistically downscaled GCM output as input for hydrological and stream temperature simulation in the Apalachicola–Chattahoochee–Flint River Basin (1961–99)

The accuracy of statistically downscaled general circulation model (GCM) simulations of daily surface climate for historical conditions (1961–99) and the implications when they are used to drive hydrologic and stream temperature models were assessed for the Apalachicola–Chattahoochee–Flint River basin (ACFB). The ACFB is a 50 000 km 2 basin located in the southeastern United States. Three GCMs were statistically downscaled, using an asynchronous regional regression model (ARRM), to ⅛° grids of daily precipitation and minimum and maximum air temperature. These ARRM-based climate datasets were used as input to the Precipitation-Runoff Modeling System (PRMS), a deterministic, distributed-parameter, physical-process watershed model used to simulate and evaluate the effects of various combinations of climate and land use on watershed response. The ACFB was divided into 258 hydrologic response units (HRUs) in which the components of flow (groundwater, subsurface, and surface) are computed in response to climate, land surface, and subsurface characteristics of the basin. Daily simulations of flow components from PRMS were used with the climate to simulate in-stream water temperatures using the Stream Network Temperature (SNTemp) model, a mechanistic, one-dimensional heat transport model for branched stream networks. The climate, hydrology, and stream temperature for historical conditions were evaluated by comparing model outputs produced from historical climate forcings developed from gridded station data (GSD) versus those produced from the three statistically downscaled GCMs using the ARRM methodology. The PRMS and SNTemp models were forced with the GSD and the outputs produced were treated as “truth.” This allowed for a spatial comparison by HRU of the GSD-based output with ARRM-based output. Distributional similarities between GSD- and ARRM-based model outputs were compared using the two-sample Kolmogorov–Smirnov (KS) test in combination with descriptive metrics such as the mean and variance and an evaluation of rare and sustained events. In general, precipitation and streamflow quantities were negatively biased in the downscaled GCM outputs, and results indicate that the downscaled GCM simulations consistently underestimate the largest precipitation events relative to the GSD. The KS test results indicate that ARRM-based air temperatures are similar to GSD at the daily time step for the majority of the ACFB, with perhaps subweekly averaging for stream temperature. Depending on GCM and spatial location, ARRM-based precipitation and streamflow requires averaging of up to 30 days to become similar to the GSD-based output. Evaluation of the model skill for historical conditions suggests some guidelines for use of future projections; while it seems correct to place greater confidence in evaluation metrics which perform well historically, this does not necessarily mean those metrics will accurately reflect model outputs for future climatic conditions. Results from this study indicate no “best” overall model, but the breadth of analysis can be used to give the product users an indication of the applicability of the results to address their particular problem. Since results for historical conditions indicate that model outputs can have significant biases associated with them, the range in future projections examined in terms of change relative to historical conditions for each individual GCM may be more appropriate.

Florida, Georgia

User's manual for the Graphical Constituent Loading Analysis System (GCLAS)

This manual describes the Graphical Constituent Loading Analysis System (GCLAS), an interactive cross-platform program for computing the mass (load) and average concentration of a constituent that is transported in stream water over a period of time. GCLAS computes loads as a function of an equal-interval streamflow time series and an equal- or unequal-interval time series of constituent concentrations. The constituent-concentration time series may be composed of measured concentrations or a combination of measured and estimated concentrations. GCLAS is not intended for use in situations where concentration data (or an appropriate surrogate) are collected infrequently or where an appreciable amount of the concentration values are censored. It is assumed that the constituent-concentration time series used by GCLAS adequately represents the true time-varying concentration. Commonly, measured constituent concentrations are collected at a frequency that is less than ideal (from a load-computation standpoint), so estimated concentrations must be inserted in the time series to better approximate the expected chemograph. GCLAS provides tools to facilitate estimation and entry of instantaneous concentrations for that purpose. Water-quality samples collected for load computation frequently are collected in a single vertical or at single point in a stream cross section. Several factors, some of which may vary as a function of time and (or) streamflow, can affect whether the sample concentrations are representative of the mean concentration in the cross section. GCLAS provides tools to aid the analyst in assessing whether concentrations in samples collected in a single vertical or at single point in a stream cross section exhibit systematic bias with respect to the mean concentrations. In cases where bias is evident, the analyst can construct coefficient relations in GCLAS to reduce or eliminate the observed bias. GCLAS can export load and concentration data in formats suitable for entry into the U.S. Geological Survey's National Water Information System. GCLAS can also import and export data in formats that are compatible with various commonly used spreadsheet and statistics programs.

Techniques and Methods