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 199 records · Page 11Linked to original sources

Assessment of long-term salinity changes in an irrigated stream-aquifer system

Changes in salinity in groundwater and surface water in the Arkansas River valley of southeastern Colorado are primarily related to irrigation practices. A solute transport model was applied to an 11-mile reach of the valley to compute salinity changes in response to spatially and temporally varying stresses. The model was calibrated in 1973 using detailed field measurements made during 1971 and 1972. In 1973 the calibrated model was used to predict that a gradual long-term increase in groundwater salinity of about 2–3% per year would occur if the observed irrigation practices continued. The study area was resampled during the winter of 1982 to help evaluate if any long-term changes in salinity are actually occurring. Nonparametric and parametric statistical tests were used to help assess the significance of observed changes in groundwater salinity. These tests indicate that a statistically significant increase in salinity occurred between the winters of 1971 and 1972 (the model calibration period). However, a comparison of the winter 1972 and winter 1982 data indicates that no significant net change in salinity has occurred during this 10-year period. An analysis of the few available historical data (1895, 1923, 1959–1961, and 1964) supports the hypothesis that groundwater salinity in this irrigated area has reached a long-term dynamic equilibrium in response to irrigation practices. The model predictions of long-term salinity increases were invalid probably because the calibration period occurred during a short-term annual trend of increasing salinity in the river (and hence in leaky irrigation canals and in applied irrigation water), which was not representative of the long-term trend.

Colorado

Statistical Approaches to Interpretation of Local, Regional, and National Highway-Runoff and Urban-Stormwater Data

Decision makers need viable methods for the interpretation of local, regional, and national-highway runoff and urban-stormwater data including flows, concentrations and loads of chemical constituents and sediment, potential effects on receiving waters, and the potential effectiveness of various best management practices (BMPs). Valid (useful for intended purposes), current, and technically defensible stormwater-runoff models are needed to interpret data collected in field studies, to support existing highway and urban-runoffplanning processes, to meet National Pollutant Discharge Elimination System (NPDES) requirements, and to provide methods for computation of Total Maximum Daily Loads (TMDLs) systematically and economically. Historically, conceptual, simulation, empirical, and statistical models of varying levels of detail, complexity, and uncertainty have been used to meet various data-quality objectives in the decision-making processes necessary for the planning, design, construction, and maintenance of highways and for other land-use applications. Water-quality simulation models attempt a detailed representation of the physical processes and mechanisms at a given site. Empirical and statistical regional water-quality assessment models provide a more general picture of water quality or changes in water quality over a region. All these modeling techniques share one common aspect-their predictive ability is poor without suitable site-specific data for calibration. To properly apply the correct model, one must understand the classification of variables, the unique characteristics of water-resources data, and the concept of population structure and analysis. Classifying variables being used to analyze data may determine which statistical methods are appropriate for data analysis. An understanding of the characteristics of water-resources data is necessary to evaluate the applicability of different statistical methods, to interpret the results of these techniques, and to use tools and techniques that account for the unique nature of water-resources data sets. Populations of data on stormwater-runoff quantity and quality are often best modeled as logarithmic transformations. Therefore, these factors need to be considered to form valid, current, and technically defensible stormwater-runoff models. Regression analysis is an accepted method for interpretation of water-resources data and for prediction of current or future conditions at sites that fit the input data model. Regression analysis is designed to provide an estimate of the average response of a system as it relates to variation in one or more known variables. To produce valid models, however, regression analysis should include visual analysis of scatterplots, an examination of the regression equation, evaluation of the method design assumptions, and regression diagnostics. A number of statistical techniques are described in the text and in the appendixes to provide information necessary to interpret data by use of appropriate methods. Uncertainty is an important part of any decisionmaking process. In order to deal with uncertainty problems, the analyst needs to know the severity of the statistical uncertainty of the methods used to predict water quality. Statistical models need to be based on information that is meaningful, representative, complete, precise, accurate, and comparable to be deemed valid, up to date, and technically supportable. To assess uncertainty in the analytical tools, the modeling methods, and the underlying data set, all of these components need be documented and communicated in an accessible format within project publications.

Open-File Report

Techniques for simulating flood hydrographs and estimating flood volumes for ungaged basins in east and west Tennessee

A dimensionless hydrograph developed for a variety of basin conditions in Georgia was tested for its applicability to streams in East and West Tennessee by comparing it to a similar dimensionless hydrograph developed for streams in East and West Tennessee. Hydrographs of observed discharge at 83 streams in East Tennessee and 38 in West Tennessee were used in the study. Statistical analyses were performed by comparing simulated (or computed) hydrographs, derived by application of the Georgia dimensionless hydrograph, and dimensionless hydrographs developed from Tennessee data, with the observed hydrographs at 50 and 75% of their peak-flow widths. Results of the tests indicate that the Georgia dimensionless hydrography is virtually the same as the one developed for streams in East Tennessee, but that it is different from the dimensionless hydrograph developed for streams in West Tennessee. Because of the extensive testing of the Georgia dimensionless hydrograph, it was determined to be applicable for East Tennessee, whereas the dimensionless hydrograph developed from data on streams in West Tennessee was determined to be applicable in West Tennessee. As part of the dimensionless hydrograph development, an average lagtime in hours for each study basin, and the volume in inches of flood runoff for each flood event were computed. By use of multiple-regression analysis, equations were developed that relate basin lagtime to drainage area size, basin length, and percent impervious area. Similarly, flood volumes were related to drainage area size, peak discharge, and basin lagtime. These equations, along with the appropriate dimensionless hydrograph, can be used to estimate a typical (average) flood hydrograph and volume for recurrence-intervals up to 100 years at any ungaged site draining less than 50 sq mi in East and West Tennessee.

Tennessee

Effective groundwater model calibration: With analysis of data, sensitivities, predictions, and uncertainty

Methods and guidelines for developing and using mathematical models Turn to Effective Groundwater Model Calibration for a set of methods and guidelines that can help produce more accurate and transparent mathematical models. The models can represent groundwater flow and transport and other natural and engineered systems. Use this book and its extensive exercises to learn methods to fully exploit the data on hand, maximize the model's potential, and troubleshoot any problems that arise. Use the methods to perform: Sensitivity analysis to evaluate the information content of data Data assessment to identify (a) existing measurements that dominate model development and predictions and (b) potential measurements likely to improve the reliability of predictions Calibration to develop models that are consistent with the data in an optimal manner Uncertainty evaluation to quantify and communicate errors in simulated results that are often used to make important societal decisions Most of the methods are based on linear and nonlinear regression theory. Fourteen guidelines show the reader how to use the methods advantageously in practical situations. Exercises focus on a groundwater flow system and management problem, enabling readers to apply all the methods presented in the text. The exercises can be completed using the material provided in the book, or as hands-on computer exercises using instructions and files available on the text's accompanying Web site. Throughout the book, the authors stress the need for valid statistical concepts and easily understood presentation methods required to achieve well-tested, transparent models. Most of the examples and all of the exercises focus on simulating groundwater systems; other examples come from surface-water hydrology and geophysics. The methods and guidelines in the text are broadly applicable and can be used by students, researchers, and engineers to simulate many kinds systems.

Book

Using sensitivity analysis in model calibration efforts

In models of natural and engineered systems, sensitivity analysis can be used to assess relations among system state observations, model parameters, and model predictions. The model itself links these three entities, and model sensitivities can be used to quantify the links. Sensitivities are defined as the derivatives of simulated quantities (such as simulated equivalents of observations, or model predictions) with respect to model parameters. We present four measures calculated from model sensitivities that quantify the observation-parameter-prediction links and that are especially useful during the calibration and prediction phases of modeling. These four measures are composite scaled sensitivities (CSS), prediction scaled sensitivities (PSS), the value of improved information (VOII) statistic, and the observation prediction (OPR) statistic. These measures can be used to help guide initial calibration of models, collection of field data beneficial to model predictions, and recalibration of models updated with new field information. Once model sensitivities have been calculated, each of the four measures requires minimal computational effort. We apply the four measures to a three-layer MODFLOW-2000 (Harbaugh et al., 2000; Hill et al., 2000) model of the Death Valley regional ground-water flow system (DVRFS), located in southern Nevada and California. D’Agnese et al. (1997, 1999) developed and calibrated the model using nonlinear regression methods. Figure 1 shows some of the observations, parameters, and predictions for the DVRFS model. Observed quantities include hydraulic heads and spring flows. The 23 defined model parameters include hydraulic conductivities, vertical anisotropies, recharge rates, evapotranspiration rates, and pumpage. Predictions of interest for this regional-scale model are advective transport paths from potential contamination sites underlying the Nevada Test Site and Yucca Mountain.

California, Nevada

Assessing the robustness of quantitative fatty acid signature analysis to assumption violations

Knowledge of animal diets can provide important insights into life history and ecology, relationships among species in a community and potential response to ecosystem change or perturbation. Quantitative fatty acid signature analysis (QFASA) is a method of estimating diets from data on the composition, or signature, of fatty acids stored in adipose tissue. Given data on signatures of potential prey, a predator diet is estimated by minimizing the distance between its signature and a mixture of prey signatures. Calibration coefficients, constants derived from feeding trials, are used to account for differential metabolism of individual fatty acids. QFASA has been widely applied since its introduction and several variants of the original estimator have appeared in the literature. However, work to compare the statistical properties of QFASA estimators has been limited. One important characteristic of an estimator is its robustness to violations of model assumptions. The primary assumptions of QFASA are that prey signature data contain representatives of all prey types consumed and the calibration coefficients are known without error. We investigated the robustness of two QFASA estimators to a range of violations of these assumptions using computer simulation and recorded the resulting bias in diet estimates. We found that the Aitchison distance measure was most robust to errors in the calibration coefficients. Conversely, the Kullback–Leibler distance measure was most robust to the consumption of prey without representation in the prey signature data. In most QFASA applications, investigators will generally have some knowledge of the prey available to predators and be able to assess the completeness of prey signature data and sample additional prey as necessary. Conversely, because calibration coefficients are derived from feeding trials with captive animals and their values may be sensitive to consumer physiology and nutritional status, their applicability to free-ranging animals is difficult to establish. We therefore recommend that investigators first make any improvements to the prey signature data that seem warranted and then base estimation on the Aitchison distance measure, as it appears to minimize risk from violations of the assumption that is most difficult to verify.

Methods in Ecology and Evolution

QADATA user's manual; an interactive computer program for the retrieval and analysis of the results from the external blind sample quality- assurance project of the U.S. Geological Survey

The U.S. Geological Survey conducts an external blind sample quality assurance project for its National Water Quality Laboratory in Denver, Colorado, based on the analysis of reference water samples. Reference samples containing selected inorganic and nutrient constituents are disguised as environmental samples at the Survey 's office in Ocala, Florida, and are sent periodically through other Survey offices to the laboratory. The results of this blind sample project indicate the quality of analytical data produced by the laboratory. This report provides instructions on the use of QADATA, an interactive, menu-driven program that allows users to retrieve the results of the blind sample quality- assurance project. The QADATA program, which is available on the U.S. Geological Survey 's national computer network, accesses a blind sample data base that contains more than 50,000 determinations from the last five water years for approximately 40 constituents at various concentrations. The data can be retrieved from the database for any user- defined time period and for any or all available constituents. After the user defines the retrieval, the program prepares statistical tables, control charts, and precision plots and generates a report which can be transferred to the user 's office through the computer network. A discussion of the interpretation of the program output is also included. This quality assurance information will permit users to document the quality of the analytical results received from the laboratory. The blind sample data is entered into the database within weeks after being produced by the laboratory and can be retrieved to meet the needs of specific projects or programs. (USGS)

Open-File Report

Methodology for Estimation of Flood Magnitude and Frequency for New Jersey Streams

Methodologies were developed for estimating flood magnitudes at the 2-, 5-, 10-, 25-, 50-, 100-, and 500-year recurrence intervals for unregulated or slightly regulated streams in New Jersey. Regression equations that incorporate basin characteristics were developed to estimate flood magnitude and frequency for streams throughout the State by use of a generalized least squares regression analysis. Relations between flood-frequency estimates based on streamflow-gaging-station discharge and basin characteristics were determined by multiple regression analysis, and weighted by effective years of record. The State was divided into five hydrologically similar regions to refine the regression equations. The regression analysis indicated that flood discharge, as determined by the streamflow-gaging-station annual peak flows, is related to the drainage area, main channel slope, percentage of lake and wetland areas in the basin, population density, and the flood-frequency region, at the 95-percent confidence level. The standard errors of estimate for the various recurrence-interval floods ranged from 48.1 to 62.7 percent. Annual-maximum peak flows observed at streamflow-gaging stations through water year 2007 and basin characteristics determined using geographic information system techniques for 254 streamflow-gaging stations were used for the regression analysis. Drainage areas of the streamflow-gaging stations range from 0.18 to 779 mi2. Peak-flow data and basin characteristics for 191 streamflow-gaging stations located in New Jersey were used, along with peak-flow data for stations located in adjoining States, including 25 stations in Pennsylvania, 17 stations in New York, 16 stations in Delaware, and 5 stations in Maryland. Streamflow records for selected stations outside of New Jersey were included in the present study because hydrologic, physiographic, and geologic boundaries commonly extend beyond political boundaries. The StreamStats web application was developed cooperatively by the U.S. Geological Survey and the Environmental Systems Research Institute, Inc., and was designed for national implementation. This web application has been recently implemented for use in New Jersey. This program used in conjunction with a geographic information system provides the computation of values for selected basin characteristics, estimates of flood magnitudes and frequencies, and statistics for stream locations in New Jersey chosen by the user, whether the site is gaged or ungaged.

Scientific Investigations Report

Peak-, mean-, and low-streamflow regional-regression equations for natural streamflow in central and western Colorado, 2019

The U.S. Geological Survey (USGS), in cooperation with the Colorado Department of Transportation, developed peak-, mean-, and low-streamflow regional-regression equations for estimating various statistics for natural streamflow in hydrologic regions of central and western Colorado. The peak-streamflow regression equations were developed using data from 418 streamgages, consisting of 15,202 years of record and a mean of approximately 36 years of record per streamgage. The mean- and low-streamflow regional-regression equations were developed using data from 323 streamgages where daily streamflow data were collected year-round. The annual exceedance-probability discharges for each streamgage were computed using the USGS software program PeakFQ. Mean monthly and 7-day minimum and maximum streamflows were computed using the USGS software program SWToolbox. Streamflow-duration values were computed using an R script. The regional-regression equations were determined using data for the period of record for a given streamgage through water year 2019. Geographic information systems datasets were used to develop 55 basin and 42 climatic characteristics, which were evaluated as candidate explanatory variables in the regression analysis. For the peak-streamflow regional-regression equations, the study area was divided into four hydrologic regions based on mean basin elevation, including the Plateau (less than 8,014 feet), Mid-Elevation (8,015 feet to 9,492 feet), Sub-Alpine (9,493 feet to 10,490 feet), and Alpine (greater than 10,490 feet) regions. For the peak-streamflow equations, the selection of basin and climatic characteristics was based on the 1-percent annual exceedance-probability discharge for each hydrologic region. For the mean streamflow, streamflow-duration values, and 7-day minimum and maximum streamflows, the study area was divided into four hydrologic regions based on river basin, including the (1) Colorado-East Slope Headwaters, (2) Green River, (3) Rio Grande, and (4) San Juan-Dolores. For mean streamflows, basin and climatic characteristics were evaluated separately for the annual period and each month for each hydrologic region. Regional regression equations published in this report are available for use in the USGS web-based program StreamStats.

Colorado

Effects of remedial grouting on the ground-water flow system at Red Rock Dam near Pella, Iowa

Previous studies have shown direct evidence of under-seepage at Red Rock Dam on the Des Moines River near Pella, Iowa. Underseepage is thought to occur primarily on the northeast side of the dam in the lower bedrock of the St. Louis Limestone, which consists of discontinuous basal evaporite beds and an overlying cavity zone. Because of concerns about the integrity of the dam, the U.S. Army Corps of Engineers initiated a remedial grouting program in September 1991. To assess the effectiveness of the remedial grouting program and to evaluate methods for future assessments, a study was conducted by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. Potentiometric surface maps of the overburden and bedrock indicate that the direction of ground-water flow on the northeast side of the dam has changed little from pre-grout to post-grout periods. A comparison of water levels, between a pre-grout date and a post-grout date, shows that water levels decreased but that the decrease may be more attributable to changes in dam operations than to remedial grouting. Waterlevel data for the same two dates indicate that a more gradual potentiometric surface exists on the northeast side of the dam than on the southwest side of the dam, which suggests that the hydraulic connection between Lake Red Rock and downgradient bedrock wells still is greater on the northeast side of the dam than on the southwest side. Hydrographs for some wells on the northeast side of the dam indicated a departure from pre-grout trends at approximately the same time grouting was initiated. To varying degrees, hydrographs for the same wells then appear to return to a trend similar to pre-grout years, possibly as a result of new flow paths developing over time after remedial grouting. Spearman correlation coefficients computed for water levels in wells, pool, and tailwater indicate that some areas on the northeast side of the dam appear to be less under the influence of changing pool elevations after grouting than before grouting. This suggests that the hydraulic connection between the Red Rock pool and some downgradient areas has decreased. Analysis of water samples collected from selected wells on the northeast side of the dam shows significant increases in sulfate concentrations beginning about the same time remedial grouting was done upgradient from the wells, possibly indicating that flow paths were cut off to these wells, thereby reducing the amount of mixing with fresh reservoir water. Observable changes in chloride concentrations or trends as a result of remedial grouting were not apparent. Analysis results for hydrogen and oxygen stable isotope samples collected since 1995 indicate large seasonal fluctuations of isotope ratios in the tailwater (assumed representative of the reservoir). Similar but more subdued fluctuations were observed at some wells, but other wells appeared to have little seasonal change. Stable sulfur isotope results indicate the presence of distinct water types between Lake Red Rock and in ground water from downgradient bedrock wells. Sulfur isotope values from samples from a bedrock well located upgradient from the grout curtain indicate a mixture of pool and ground water, whereas samples from downgradient overburden wells have values similar to the pool. Samples from the bedrock wells downgradient from the grout curtain have sulfur isotope values similar to a value obtained from analysis of a gypsum and anhydrite core sample. Hydrographs, statistical analysis of waterlevel data, and water-chemistry data suggest that underseepage on the northeast side of the dam has been reduced but not completely eliminated. Some areas appear to have been affected to a greater degree and for a longer period of time than other areas. Future monitoring of water levels, water chemistry, and stable isotopes can aid in the evaluation of the long-term effectiveness of remedial grouting.

Iowa

Magnitude, frequency, and trends of floods at gaged and ungaged sites in Washington, based on data through water year 2014

An investigation into the magnitude and frequency of floods in Washington State computed the annual exceedance probability (AEP) statistics for 648 U.S. Geological Survey unregulated streamgages in and near the borders of Washington using the recorded annual peak flows through water year 2014. This is an updated report from a previous report published in 1998 that used annual peak flows through the water year 1996. New in this report, a regional skew coefficient was developed for the Pacific Northwest region that includes areas in Oregon, Washington, Idaho and western Montana within the Columbia River drainage basin south of the United States-Canada border, the coastal areas of Oregon and western Washington, and watersheds draining into Puget Sound, Washington. The skew coefficient is an important term in the Log Pearson Type III equation used to define the distribution of the log-transformed annual peaks. The Expected Moments Algorithm was used to fit historical and censored peak-flow data to the log Pearson Type III distribution. A Multiple Grubb-Beck test was employed to censor low outliers of annual peak flows to improve on the frequency distribution. This investigation also includes a section on observed trends in annual peak flows that showed significant trends ( p -value < 0.05) in 21 of 83 long-term sites, but with small magnitude Kendall tau values suggesting a limited monotonic trend in the time series of annual peaks. Most of the sites with a significant trend in western Washington were positive and all the sites with significant trends (three sites) in eastern Washington were negative. Multivariate regression analysis with measured basin characteristics and the AEP statistics at long-term, unregulated, and un-urbanized (defined as drainage basins with less than 5 percent impervious land cover for this investigation) streamgages within Washington and some in Idaho and Oregon that are near the Washington border was used to develop equations to estimate AEP statistics at ungaged basins. Washington was divided into four regions to improve the accuracy of the regression equations; a set of equations for eight selected AEPs and for each region were constructed. Selected AEP statistics included the annual peak flows that equaled or exceeded 50, 20, 10, 4, 2, 1, 0.5 and 0.2 percent of the time equivalent to peak flows for peaks with a 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals, respectively. Annual precipitation and drainage area were the significant basin characteristics in the regression equations for all four regression regions in Washington and forest cover was significant for the two regression regions in eastern Washington. Average standard error of prediction for the regional regression equations ranged from 70.19 to 125.72 percent for Regression Regions 1 and 2 on the eastern side of the Cascade Mountains and from 43.22 to 58.04 percent for Regression Regions 3 and 4 on the western side of the Cascade Mountains. The pseudo coefficient of determination (where a value of 100 signifies a perfect regression model) ranged from 68.39 to 90.68 for Regression Regions 1 and 2, and 92.35 to 95.44 for Regions 3 and 4. The calculated AEP statistics for the streamgages and the regional regression equations are expected to be incorporated into StreamStats after the publication of this report. StreamStats is the interactive Web-based map tool created by the U.S. Geological Survey to allow the user to choose a streamgage and obtain published statistics or choose ungaged locations where the program automatically applies the regional regression equations and computes the estimates of the AEP statistics.

Washington

Concentrations of Escherichia coli in streams in the Kankakee and lower Wabash River watersheds in Indiana, June-September 1999

Water samples collected from 58 surface- water sites in the Kankakee and Lower Wabash River Watersheds from June through September 1999 were analyzed for concentrations of Escherichia coli bacteria. Each sitewas sampled five times in a 30-day period. Twentynine sites were sampled during June and July, and 29 different sites were sampled during August and September. A five-sample geometric mean of concentrations was computed for each site. Concentrations of Escherichia coli (E. coli) in 126 of the 289 samples exceeded the State of Indiana single-sample standard of 235 colonies per 100 milliliters for waters used for recreation. Concentrations in samples from 38 of the 58 sites exceeded the State of Indiana standard for a five-sample geometric mean of 125 colonies per 100 milliliters for waters used for recreation. Ten of the 58 sites were at or near U.S. Geological Survey streamflow-gaging stations. Based on records from the streamflowgaging stations, 18 percent of the samples collected at these sites were collected at streamflows above the median daily discharge for each station. E. coli concentrations and turbidity measurements collected during 1999 were analyzed in concert with similar concentration and turbidity data collected in 1998 at streams within the Upper Wabash River Watershed in Indiana to investigate the relation between concentrations of bacteria and turbidity. The analysis indicated a statistically significant correlation between concentrations of E. coli and turbidity. If the turbidity was greater than 83 nephelometric turbidity units, the E. coli concentration always exceeded the singlesample standard. If, however, the turbidity was less than 83 nephelometric turbidity units, concentrations of E. coli were not always below the single-sample standard.

Illinois, Indiana

U.S. Geological Survey national shoreline change— Summary statistics for updated vector shorelines (1800s–2010s) and associated shoreline change data for the Georgia and Florida coasts

Rates of shoreline change have been updated for the open-ocean sandy coastlines of Georgia and Florida as part of the U.S. Geological Survey’s Coastal Change Hazards programmatic focus. This work was formerly within the National Assessment of Shoreline Change project. Shorelines were compiled from the original report published in 2005, recent update reports, and additional light detection and ranging (lidar) shorelines which were extracted from lidar data collected prior to and following Hurricane Irma, which made landfall in September 2017. These shorelines were used to compute long- and short-term rates that incorporate the proxy-datum bias on a transect-by-transect basis. The proxy-datum bias accounts for the unidirectional onshore bias of proxy-based high water line shorelines relative to datum-based mean high water shorelines. In this study, the coast of Georgia exhibited the highest average rates of erosion and accretion in both the long term (approximately 150 years) and the short term (approximately 30 years). Shoreline positions from the mid-1800s through 2018 were used to update the shoreline change rates for Florida and Georgia using the Digital Shoreline Analysis System (DSAS) software.

Florida, Georgia

Flood-frequency comparison from 1995 to 2016 and trends in peak streamflow in Arkansas, water years 1930–2016

In 2016, the U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers and the Federal Emergency Management Agency, began a study in Arkansas to investigate possible increasing trends in annual peak streamflow data and the possible resulting increase in the annual exceedance probability flood (AEPF) predictions. Temporal trends of peak streamflow were investigated at 15 selected streamgages on unregulated streams in Arkansas having 30 or more years of peak streamflow data through the 2016 water year. For the period of record at each streamgage, the Mann-Kendall trend test indicated that 14 of the 15 streamgages had no statistically significant peak streamflow trends and 1 streamgage had a statistically significant decreasing peak streamflow trend. Visual examination of the locally estimated scatterplot smoothing technique trend lines of the peak streamflow data indicated a possible increasing peak streamflow trend at 8 of the 15 streamgages since the 1990s. A sequential series analysis of the 1-percent AEPF at each of the 15 selected streamgages was completed by selecting an initial subset of the oldest peak streamflow data from each site to estimate the initial 1-percent AEPF. This initial peak streamflow data subset was subsequently appended with 10-year increments of additional peak streamflow data until the full period of peak streamflow data was analyzed. The maximum increase in the 1-percent AEPF was 113 percent, and the maximum decrease was 31.9 percent. Percentage differences between the AEPFs derived from regional regression equations presented in the 1995 and 2016 Arkansas flood-frequency reports were compared. The average percentage differences for the 74 selected locations indicate that the 4-, 2-, 1-, and 0.2-percent AEPFs computed using the 2016 regional regression equations were higher by 3.52, 5.10, 8.59, and 13.31 percent, respectively (25-, 50-, 100-, and 500-year recurrence interval floods), than the same percentage AEPFs computed using the 1995 regional regression equations. The average percentage differences between the 1995 and 2016 AEPFs for the 10-percent AEPF (10-year recurrence interval flood) resulted in 2016 AEPF predictions being 0.41 percent higher. For the 50- and 20-percent AEPFs (2- and 5-year recurrence interval floods), the 2016 AEPFs were less than the 1995 AEPFs by 2.53 and 0.31 percent, respectively.

Arkansas

User's Guide, software for reduction and analysis of daily weather and surface-water data: Tools for time series analysis of precipitation, temperature, and streamflow data

The software described here is used to process and analyze daily weather and surface-water data. The programs are refinements of earlier versions that include minor corrections and routines to calculate frequencies above a threshold on an annual or seasonal basis. Earlier versions of this software were used successfully to analyze historical precipitation patterns of the Mojave Desert and the southern Colorado Plateau regions, ecosystem response to climate variation, and variation of sediment-runoff frequency related to climate (Hereford and others, 2003; 2004; in press; Griffiths and others, 2006). The main program described here (Day_Cli_Ann_v5.3) uses daily data to develop a time series of various statistics for a user specified accounting period such as a year or season. The statistics include averages and totals, but the emphasis is on the frequency of occurrence in days of relatively rare weather or runoff events. These statistics are indices of climate variation; for a discussion of climate indices, see the Climate Research Unit website of the University of East Anglia (http://www.cru.uea.ac.uk/projects/stardex/) and the Climate Change Indices web site (http://cccma.seos.uvic.ca/ETCCDMI/indices.html). Specifically, the indices computed with this software are the frequency of high intensity 24-hour rainfall, unusually warm temperature, and unusually high runoff. These rare, or extreme events, are those greater than the 90th percentile of precipitation, streamflow, or temperature computed for the period of record of weather or gaging stations. If they cluster in time over several decades, extreme events may produce detectable change in the physical landscape and ecosystem of a given region. Although the software has been tested on a variety of data, as with any software, the user should carefully evaluate the results with their data. The programs were designed for the range of precipitation, temperature, and streamflow measurements expected in the semiarid Southwest United States. The user is encouraged to review the examples provided with the software. The software is written in Fortran 90 with Fortran 95 extensions and was compiled with the Digital Visual Fortran compiler version 6.6. The executables run on Windows 2000 and XP, and they operate in a MS-DOS console window that has only very simple graphical options such as font size and color, background color, and size of the window. Error trapping was not written into the programs. Typically, when an error occurs, the console window closes without a message.

Open-File Report

Subsampling large-scale digital elevation models to expedite geospatial analyses in coastal regions

Large-area, high-resolution digital elevation models (DEMs) created from light detection and ranging (LIDAR) and/or multibeam echosounder data sets are commonly used in many scientific disciplines. These DEMs can span thousands of square kilometers, typically with a spatial resolution of 1 m or finer, and can be difficult to process and analyze without specialized computers and software. Such DEMs often can be subsampled to expedite analysis with negligible impact on results for large-scale geospatial analyses. Subsampling can be achieved by creating a grid of points that specify the locations from which to extract elevation values from the DEM. This paper presents a method that can be used to accurately perform subsampling of large-scale, high-resolution DEMs using GIS software. This subsampling method was applied to two LIDAR-derived DEMs encompassing 242 km 2 of the northern Florida Reef Tract as an example application and to test subsampling accuracy. Results indicate that subsampling 1-m-resolution DEMs using a 2-m-spaced grid results in no significant difference in mean elevation or other basic statistics for analyses performed over multiple spatial scales ranging from 1 km 2 to 242 km 2 .

Journal of Coastal Research

Trends in selected streamflow statistics at 19 long-term streamflow-gaging stations indicative of outflows from Texas to Arkansas, Louisiana, Galveston Bay, and the Gulf of Mexico, 1922-2009

Trends in selected streamflow statistics during 1922-2009 were evaluated at 19 long-term streamflow-gaging stations considered indicative of outflows from Texas to Arkansas, Louisiana, Galveston Bay, and the Gulf of Mexico. The U.S. Geological Survey, in cooperation with the Texas Water Development Board, evaluated streamflow data from streamflow-gaging stations with more than 50 years of record that were active as of 2009. The outflows into Arkansas and Louisiana were represented by 3 streamflow-gaging stations, and outflows into the Gulf of Mexico, including Galveston Bay, were represented by 16 streamflow-gaging stations. Monotonic trend analyses were done using the following three streamflow statistics generated from daily mean values of streamflow: (1) annual mean daily discharge, (2) annual maximum daily discharge, and (3) annual minimum daily discharge. The trend analyses were based on the nonparametric Kendall's Tau test, which is useful for the detection of monotonic upward or downward trends with time. A total of 69 trend analyses by Kendall's Tau were computed - 19 periods of streamflow multiplied by the 3 streamflow statistics plus 12 additional trend analyses because the periods of record for 2 streamflow-gaging stations were divided into periods representing pre- and post-reservoir impoundment. Unless otherwise described, each trend analysis used the entire period of record for each streamflow-gaging station. The monotonic trend analysis detected 11 statistically significant downward trends, 37 instances of no trend, and 21 statistically significant upward trends. One general region studied, which seemingly has relatively more upward trends for many of the streamflow statistics analyzed, includes the rivers and associated creeks and bayous to Galveston Bay in the Houston metropolitan area. Lastly, the most western river basins considered (the Nueces and Rio Grande) had statistically significant downward trends for many of the streamflow statistics analyzed.

Arkansas, Colorado, Louisiana, Oklahoma, New Mexic

Use of set blanks in reporting pesticide results at the U.S. Geological Survey National Water Quality Laboratory, 2001-15

Executive Summary Background .—Pesticide results from the U.S. Geological Survey (USGS) National Water Quality Laboratory (NWQL) are used for water-quality assessments by many agencies and organizations. The USGS is committed to providing data of the highest possible quality to the consumers of its data. A cooperator’s inquiries about specific pesticide detections in water revealed potential laboratory contamination issues for some results. Consequently, the USGS conducted an extensive evaluation of potential low-level contamination related to processing or analysis of water-quality samples at NWQL for 21 pesticide compounds of interest to the cooperator. This is the most comprehensive study of NWQL pesticide quality-control (QC) results to date. Purpose and scope .—The purpose of this study was to document protocols used by the NWQL to censor pesticide results and to determine the effects of laboratory contamination—as determined from detections in laboratory set blanks—on pesticide detections in groundwater and surface-water samples. More than 30,000 pesticide results from 113 selected batches of samples (2 percent or less of total batches) analyzed by the NWQL during the 15 years from 2001 to 2015 were reviewed. All laboratory results from the selected batches, including results from environmental (surface water and groundwater) and QC (set-blank, blind-blank, and blind-spike) samples, were evaluated. The study includes results for 21 pesticide compounds analyzed in groundwater and surface-water samples collected across the United States. Eleven pesticide compounds were analyzed by a gas chromatography/mass spectrometry method and 10 compounds by a liquid chromatography/mass spectrometry method. Objectives and methods .—The objectives of this study were to (1) determine the characteristics of laboratory contamination over time, (2) compare distributions of pesticide results in set blanks with distributions in environmental samples, (3) evaluate the potential for false-positive and false-negative reporting of results, and (4) evaluate the effects of reevaluating historical pesticide results using 2017 compound identification protocols on detections of pesticides in groundwater and surface-water samples. The 113 instrument batches selected for this study contained detections of one or more of the 21 pesticide compounds in set blanks or were among those batches with the highest pesticide detection frequencies in set blanks. As a result, the dataset for this study was targeted toward pesticides and batches with laboratory contamination. The objectives were addressed by statistically comparing environmental and set-blank results; computing moving averages of set-blank detection frequencies to identify periods of episodic contamination; and using summary statistics, tabular summaries, and graphical approaches, such as time-series plots and cumulative distribution functions. Results .— Objective 1 : Laboratory contamination, as determined by pesticide detections in set blanks, was found in 13 percent of set-blank results from the 113 targeted batches included in this study (as compared to 6 percent of set-blank results from all 7,620 batches analyzed during the study period). It is estimated that 92 percent of the laboratory contamination during the study period was episodic, meaning that it occurred during discrete periods of time. All 21 of the targeted pesticide compounds had periods of episodic contamination, with most episodes ranging in duration from about 1 to 8 months. The remaining 8 percent of laboratory contamination was random or from a known source (deterministic). Objective 2 : For some compounds, graphs of cumulative distribution functions of the entire distributions of set-blank and environmental samples overlap, suggesting that there is no difference in the distributions of the two types of samples. However, time-series graphs show that detections in set blanks often occur at different times (sometimes separated by years) than detections in environmental samples, indicating clear differences in those distributions, and indicating the importance of evaluating the timing of detections in all sample types. For most compounds detected in set-blank and environmental samples, detection frequencies were significantly greater in set blanks than in groundwater or surface-water samples ( p <0.05). There are several explanations for this finding, including that the 113 batches of samples chosen for this study targeted batches with detections in set blanks or that detections in set-blank samples were historically determined with less stringent identification criteria than for environmental samples (groundwater and surface-water samples). Objective 3 : The false-positive and false-negative rates from blind samples submitted during the study period by the USGS Quality Systems Branch generally were less than 1 and 5 percent, respectively, for the 21 pesticides. The only compound with a false-positive rate greater than 1 percent was flumetsulam (2.6 percent), indicating that there is a higher likelihood of flumetsulam being reported as a detection when it is not present in an environmental sample compared with the reporting of other compounds. Objective 4 : Altogether, for data in targeted batches, NWQL would have reported 0.1 percent of results from groundwater samples and 1.4 percent of results from surface-water samples differently if 2017 identification protocols were applied to historical pesticide results. In most of these cases, detections observed in historical results would change to nondetections. The small percentages of changes that would occur if historical data were reevaluated indicate that historical protocols used by the NWQL to identify detections in environmental samples were robust and produced results that are predominantly consistent with current [2017] practices. Conclusions .—The NWQL produces high-quality pesticide results at environmentally relevant concentrations. NWQL identification protocols and censoring practices are largely effective at minimizing the reporting of false-positive and false-negative results. Laboratory contamination, when it occurred, tended to occur in episodes; thus, evaluating the timing and magnitude of detections in set blanks relative to detections in environmental samples was determined to be an important consideration for analysis of environmental results. Because NWQL censoring practices do not address all types and occurrences of laboratory contamination, options for additional censoring practices are provided for data users with more specific or stringent data-quality objectives. The methods used to analyze the 21 compounds for this report can similarly be applied to all 173 pesticide compounds that were analyzed by the NWQL during the same time period. This study also has helped to identify potential improvements in reporting USGS data, such as conducting more frequent review of set-blank datasets.

Scientific Investigations Report