SEARCH · Search USGS

Results for “Computers and Geology”

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

SourceResearch connectionFrom yearThrough year

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

At least 541 records · Page 30Linked to original sources

USGS1985

Spirit Lake dam-failure flood routing assessment

Potential clear-water floodflows resulting from uncontrolled breaching by Spirit Lake of a debris dam deposited by the May 18 , 1980, eruption of Mount St. Helens, Washington, were evaluated. U.S. Geological Survey dam-break model K-634 was utilized, first to compute clear-water flood hydrographs for various hypothetical breach scenarios, and then to hydraulically route them downstream to the mouth of the Toutle River. Dam-break computations were obtained for 12 breach-development scenarios. Eight of the scenarios modeled breaches caused by overtopping of the June 1982 crest of the debris dam, while the other four modeled breaches resulting from overtopping of the debris dam at a possible future crest altitude of 3,490 feet. Equal numbers of scenarios, six each, were modeled for breach development durations of 0.25 hour and 1.0 hour. Peak discharges and times of arrival at selected locations for scenarios with breach development durations of 1.0 hour are presented in tables. Peak discharges computed for dam-break scenarios with breach crest widths of 400 feet are almost twice the magnitude of those computed for scenarios identical in all other respects, but with widths of 200 feet. Peak discharges computed for breaches caused by overtopping of the debris dam at an altitude of 3,531.8 feet were generally found to be roughly five times greater than those caused by overtopping at an altitude of 3,490 feet, if the breaches had identical crest elevations and crest widths. Decreases in peak discharge from Spirit Lake to the mouth of the Toutle River ranged from 16 to 26 percent for the scenarios modeled. This lack of significant attenuation is primarily due to the large volume and surface area of Spirit lake (in excess of 360,000 acre-feet and 3,000 acres, respectively, for all scenarios modeled), which result in the continued discharge of high flows long after breach development is complete. The degree of attenuation is also minimized by the generally steep narrow valleys through which the North Fork Toutle and Toutle Rivers flow. Elapsed time from beginning of breach development to arrival of peak discharges at the mouth of the Toutle River ranged from 4.3 to 7.4 hours for the scenarios modeled. (USGS)

Open-File Report↗

USGS1991

System requirements specification for the U.S. Geological Survey's National Water Information System II

The U.S. Geological Survey's Water Resources Division is designing and developing a new computer software system for processing and storing hydrologic data. This system, the National Water Information System (NWIS-II) will replace the current water-data and information systems: the National Water Data Storage and Retrieval System (WATSTORE), the National Water Data Exchange (NAWDEX), and the National Water-Use Information System, as well as the current National Water Information System (NWIS-I). The new system will be utilized across the Nation on the Distributed Information System (DIS-II), a network of 32-bit microcomputers on a wide-area network of desktop workstations and related hardware. The purpose of the NWIS-II System Requirements Specification is to integrate requirements described by eight User Groups representing the disciplines for Surface Water, Ground Water, Water Quality, Water Sediment, Water Use, Biology, Geographic Information, and National Water Data Exchange. The requirements specified in this document will serve as the framework for subsequent design and development activities throughout the life cycle of NWIS-II to ensure that the development products reflect the needs of the users. Specification details and associated data-base plans are contained in the description of the integrated functional requirements, the logical data model and the data dictionary, the plan for transferring data from the existing system to NWIS-II, and the descriptions of performance requirements and design constraints. The NWIS-II will be consistent across hydrologic disciplines and the other major Division software systems. The color-graphics desktop environment of NWIS-II will utilize multi-window and multi-tasking capabilities with multilevel of help for the users. Data will be distributed among nodes and be automatically accessible during input, editing, verification and analysis of data. Hydrologic data will be cataloged independent of the Environmental Protection Agency's storage and retrieval (STORET) codes. Character, digital, and graphical data will be entered into NWIS-II from digital and analog recorders, external files, or keyboard through standard and user-defined input forms. Data will be verified manually or automatically during input, retrieval, and editing of data. Numerous reference lists will be maintained so that entry of codes on forms will not be required. User-defined series of computer actions for standard programs will be supported for processing basic data and other frequently repeated processes. Data will be verified and exchanged through links to other DIS-II software, including a Geographic Information System, and selected external software. Quality assurance of data will include aging of data through change on status type and tracking access violations. Output of data will be to a variety of media in standard and user-defined formats. The data base will be preserved by operations such as backup, recovery, audit trails, archiving, and data histories. In addition, an index of water data from USGS and other agencies will be maintained.

Open-File Report↗

USGS1999

The Silent Canyon caldera complex — A three-dimensional model based on drill-hole stratigraphy and gravity inversion

The structural framework of Pahute Mesa, Nevada, is dominated by the Silent Canyon caldera complex, a buried, multiple collapse caldera complex. Using the boundary surface between low density Tertiary volcanogenic rocks and denser granitic and weakly metamorphosed sedimentary rocks (basement) as the outer fault surfaces for the modeled collapse caldera complex, it is postulated that the caldera complex collapsed on steeply- dipping arcuate faults two, possibly three, times following eruption of at least two major ash-flow tuffs. The caldera and most of its eruptive products are now deeply buried below the surface of Pahute Mesa. Relatively low-density rocks in the caldera complex produce one of the largest gravity lows in the western conterminous United States. Gravity modeling defines a steep sided, cup-shaped depression as much as 6,000 meters (19,800 feet) deep that is surrounded and floored by denser rocks. The steeply dipping surface located between the low-density basin fill and the higher density external rocks is considered to be the surface of the ring faults of the multiple calderas. Extrapolation of this surface upward to the outer, or topographic rim, of the Silent Canyon caldera complex defines the upper part of the caldera collapse structure. Rock units within and outside the Silent Canyon caldera complex are combined into seven hydrostratigraphic units based on their predominant hydrologic characteristics. The caldera structures and other faults on Pahute Mesa are used with the seven hydrostratigraphic units to make a three-dimensional geologic model of Pahute Mesa using the "EarthVision" (Dynamic Graphics, Inc.) modeling computer program. This method allows graphic representation of the geometry of the rocks and produces computer generated cross sections, isopach maps, and three-dimensional oriented diagrams. These products have been created to aid in visualizing and modeling the ground-water flow system beneath Pahute Mesa.

Nevada↗

USGS2016

Tropical storm Irene flood of August 2011 in northwestern Massachusetts

A Presidential disaster was declared in northwestern Massachusetts, following flooding from tropical storm Irene on August 28, 2011. During the storm, 3 to 10 inches of rain fell on soils that were susceptible to flash flooding because of wet antecedent conditions. The gage height at one U.S. Geological Survey streamgage rose nearly 20 feet in less than 4 hours because of the combination of saturated soils and intense rainfall. On August 28, 2011, in the Deerfield and Hoosic River Basins in northwestern Massachusetts, new peaks of record were set at six of eight U.S. Geological Survey long-term streamgages with 46 to 100 years of record. Additionally, high-water marks were surveyed and indirect measurements of peak discharge were calculated at two discontinued streamgages in the Deerfield and Hoosic River Basins with 24 and 61 years of record, respectively. This data resulted in new historic peaks of record at the two discontinued streamgages from tropical storm Irene. Peak flows that resulted from tropical storm Irene (August 28, 2011) were determined at the U.S. Geological Survey streamgages by using stage-discharge rating curves and indirect computation methods. For six streamgages, indirect computation methods were used to compute the peak flows. Peak flows from tropical storm Irene had annual exceedance probabilities (AEPs) that ranged from 5.4 percent to less than 0.2 percent at 10 streamgages in northwestern Massachusetts. Discharges calculated for select AEPs as a part of this study were compared with discharges published for the same AEPs in the effective Federal Emergency Management Agency flood insurance studies (FISs) for communities in the study area. Discharges estimated for the 10-, 2-, 1-, and 0.2-percent AEPs at two streamgages on the main stem of the Deerfield River ranged from about 3 percent lower to 14 percent higher than discharges in the FISs. AEP discharges calculated for two streamgages on tributaries to the Deerfield River were 27 to 89 percent higher than the FISs. For the four streamgages in the Hoosic River Basin, the 10-, 2-, 1-, and 0.2-percent AEP discharges calculated ranged from about 33 percent lower to 5 percent higher than the FISs. The simulated 1-percent AEP discharge water-surface elevations (nonregulatory) from recent (2015–16) hydraulic models for river reaches in the study area, which include the Deerfield, Green, and North Rivers in the Deerfield River Basin and the Hoosic River in the Hoosic River Basin, were compared with water-surface profiles in the FISs. The water-surface elevation comparisons were generally done downstream and upstream from bridges, dams, and major tributaries. The simulated 1-percent AEP discharge water-surface elevations of the recent hydraulic studies averaged 2.2, 2.3, 0.3, and 0.7 ft higher than water-surface elevations in the FISs for the Deerfield, Green, North, and Hoosic Rivers, respectively. The differences in water-surface elevations between the recent (2015–16) hydraulic studies and the FISs likely are because of (1) improved land elevation data from light detection and ranging (lidar) data collected in 2012, (2) detailed surveying of hydraulic structures and cross sections throughout the river reaches in 2012–13 (reflecting structure and cross section changes during the last 30–35 years), (3) updated hydrology analyses (30–35 water years of additional peak flow data at streamgages), and (4) high-water marks from the 2011 tropical storm Irene flood being used for model calibration.

Massachusetts↗

USGS2017

Methods for estimating selected low-flow frequency statistics and mean annual flow for ungaged locations on streams in North Georgia

The U.S. Geological Survey, in cooperation with the Georgia Department of Natural Resources, Environmental Protection Division, developed regional regression equations for estimating selected low-flow frequency and mean annual flow statistics for ungaged streams in north Georgia that are not substantially affected by regulation, diversions, or urbanization. Selected low-flow frequency statistics and basin characteristics for 56 streamgage locations within north Georgia and 75 miles beyond the State’s borders in Alabama, Tennessee, North Carolina, and South Carolina were combined to form the final dataset used in the regional regression analysis. Because some of the streamgages in the study recorded zero flow, the final regression equations were developed using weighted left-censored regression analysis to analyze the flow data in an unbiased manner, with weights based on the number of years of record. The set of equations includes the annual minimum 1- and 7-day average streamflow with the 10-year recurrence interval (referred to as 1Q10 and 7Q10), monthly 7Q10, and mean annual flow. The final regional regression equations are functions of drainage area, mean annual precipitation, and relief ratio for the selected low-flow frequency statistics and drainage area and mean annual precipitation for mean annual flow. The average standard error of estimate was 13.7 percent for the mean annual flow regression equation and ranged from 26.1 to 91.6 percent for the selected low-flow frequency equations. The equations, which are based on data from streams with little to no flow alterations, can be used to provide estimates of the natural flows for selected ungaged stream locations in the area of Georgia north of the Fall Line. The regression equations are not to be used to estimate flows for streams that have been altered by the effects of major dams, surface-water withdrawals, groundwater withdrawals (pumping wells), diversions, or wastewater discharges. The regression equations should be used only for ungaged sites with drainage areas between 1.67 and 576 square miles, mean annual precipitation between 47.6 and 81.6 inches, and relief ratios between 0.146 and 0.607; these are the ranges of the explanatory variables used to develop the equations. An attempt was made to develop regional regression equations for the area of Georgia south of the Fall Line by using the same approach used during this study for north Georgia; however, the equations resulted with high average standard errors of estimates and poorly predicted flows below 0.5 cubic foot per second, which may be attributed to the karst topography common in that area. The final regression equations developed from this study are planned to be incorporated into the U.S. Geological Survey StreamStats program. StreamStats is a Web-based geographic information system that provides users with access to an assortment of analytical tools useful for water-resources planning and management, and for engineering design applications, such as the design of bridges. The StreamStats program provides streamflow statistics and basin characteristics for U.S. Geological Survey streamgage locations and ungaged sites of interest. StreamStats also can compute basin characteristics and provide estimates of streamflow statistics for ungaged sites when users select the location of a site along any stream in Georgia.

Georgia↗

USGS2000

Factors affecting nutrient trends in major rivers of the Chesapeake Bay Watershed

Trends in nutrient loads and flow-adjusted concentrations in the major rivers entering Chesapeake Bay were computed on the basis of water-quality data collected between 1985 and 1998 at 29 monitoring stations in the Susquehanna, Potomac, James, Rappahannock, York, Patuxent, and Choptank River Basins. Two computer models-the Chesapeake Bay Watershed Model (WSM) and the U.S. Geological Survey's 'Spatially Referenced Regressions on Watershed attributes' (SPARROW) Model-were used to help explain the major factors affecting the trends. Results from WSM simulations provided information on temporal changes in contributions from major nutrient sources, and results from SPARROW model simulations provided spatial detail on the distribution of nutrient yields in these basins. Additional data on nutrient sources, basin characteristics, implementation of management practices, and ground-water inputs to surface water were analyzed to help explain the trends. The major factors affecting the trends were changes in nutrient sources and natural variations in streamflow. The dominant source of nitrogen and phosphorus from 1985 to 1998 in six of the seven tributary basins to Chesapeake Bay was determined to be agriculture. Because of the predominance of agricultural inputs, changes in agricultural nutrient sources such as manure and fertilizer, combined with decreases in agricultural acreage and implementation of best management practices (BMPs), had the greatest impact on the trends in flow-adjusted nutrient concentrations. Urban acreage and population, however, were noted to be increasing throughout the Chesapeake Bay Watershed, and as a result, delivered loads of nutrients from urban areas increased during the study period. Overall, agricultural nutrient management, in combination with load decreases from point sources due to facility upgrades and the phosphate detergent ban, led to downward trends in flow-adjusted nutrient concentrations atmany of the monitoring stations in the watershed. The loads of nutrients, however, were not reduced significantly at most of the monitoring stations. This is due primarily to higher streamflow in the latter years of the monitoring period, which led to higher loading in those years. Results of this study indicate a need for more detailed information on BMP effectiveness under a full range of hydrologic conditions and in different areas of the watershed; an internally consistent fertilizer data set; greater consideration of the effects of watershed processes on nutrient transport; a refinement of current modeling efforts; and an expansion of the non-tidal monitoring network in the Chesapeake Bay Watershed.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

USGS2023

Magnitude conversion and earthquake recurrence rate models for the central and eastern United States

Development of Seismic Source Characterization (SSC) models, which is an essential part of Probabilistic Seismic Hazard Analyses (PSHA), can help forecast the temporal and spatial distribution of future damaging earthquakes (𝑀 w ≥ 5) in seismically active regions. Because it is impossible to associate all earthquakes with known faults, seismic source models for PSHA often include sources of diffuse seismicity in which future earthquake scenarios are not localized on mapped faults. These sources of diffuse seismicity are referred to as area source zones, distributed seismicity zones, or just source zones. During the early years of PSHA studies, it was assumed that earthquakes in seismotectonic zones have (1) uniform spatial distribution, (2) Poisson temporal distribution, and (3) exponential magnitude distribution (NRC, 2012). In seismically active regions (e.g., the Western United States), where active faults are readily identified, models of the spatial distribution of earthquakes include both the fault source geometries and the distributed seismicity (background) source zones. Source characterization of active faults is complemented by paleoseismic studies with estimates of earthquake magnitudes, dates of occurrences, and slip rates, which provide important information for PSHA studies. In the Central and Eastern United States (CEUS) very few Quaternary-active faults have the requisite information for use in PSHA (i.e., fault geometry and dimensions, event rates or slip rates, etc.), and we lack knowledge about the causative faults for most observed seismicity in the region. As a result, area source zones are frequently used in site-specific PSHA in the CEUS to represent diffuse seismicity that cannot be associated with faults. However, there are examples of active fault sources in the CEUS, such as the Meers fault, the Cheraw fault, and New Madrid region, where individual faults can be characterized. The source characterization models for background seismicity are based, to a large extent, on an assumption that spatial distribution of historical and recorded seismicity will not change substantially for time periods of interest for PSHA (approximately the next 50-100 years for engineered structures). Furthermore, studies such as those by Kafka (2007, 2009) found a correlation between the locations of small- to moderate-magnitude earthquakes and the locations of large-magnitude earthquakes, indicating that we can, with some level of confidence, use the spatial pattern of smaller earthquakes to forecast the future pattern of damaging earthquakes. Within background seismicity zones, the earthquake rate forecast is developed using spatial smoothing of the small to moderate magnitude events in earthquake catalogs. Different methodologies are used for this purpose and can predict varying distributions of seismicity rates. This in turn affects the results of a seismic hazard analysis. The U.S. Geological Survey (USGS) and Nuclear Regulatory Commission (NRC) use different methods for computing spatially smoothed seismicity rates in the CEUS; the USGS uses kernel-based spatial smoothing methods in developing the National Seismic Hazard Model (NSHM), and the method adopted in the Central and Eastern United States Seismic Source Characterization (CEUS-SSC) project is used when evaluating seismic hazard for nuclear power plant siting. These methods are described and the impact on seismic hazard are evaluated in this Research Information Letter (RIL). Another important input to estimating the rate of distributed seismicity is event magnitudes listed in earthquake catalogs. A substantial source of uncertainty in catalogs is the magnitude assigned to a given earthquake. Numerous different magnitude types exist, with each magnitude type computed in a different way. Therefore, for the sake of consistency, both the CEUS-SSC and the USGS NSHM have attempted to assemble a complete catalog with a uniform magnitude determination. To this end, moment magnitude, 𝑀 w , which is a physics-based measurement, has been adopted as the standard. However, 𝑀 w was not computed routinely until the past few decades. To address this issue, the CEUS-SSC conducted extensive analyses to determine conversion equations from which to take a routinely computed network (e.g., 𝑀 L or 𝑚 bLg ) and convert it into 𝑀 w . Another issue with using 𝑀 w is that it becomes increasingly difficult to compute for earthquakes with 𝑀 less than ~4. This study investigates the effects of moment magnitude estimation and spatial smoothing methods on estimation of the earthquake rate forecast and on seismic hazard. We investigate the validity of the magnitude conversion equations and their associated uncertainties by applying them to a case study for induced earthquakes in southern Kansas and northern Oklahoma, and summarize the use of the decay of the seismic coda to estimate 𝑀 w for small earthquakes (𝑀 w < 4. Furthermore, the study documents a comparison and assessment of background seismicity smoothing methods implemented by the USGS for the NSHM and used by the CEUS-SSC for siting nuclear facilities based on probabilistic seismic hazard estimates from multiple source zones in the CEUS and for multiple sites.

central and eastern United States↗

USGS1990

Preconditioned Conjugate-Gradient 2 (PCG2), a computer program for solving ground-water flow equations

This report documents PCG2: a numerical code to be used with the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model. PCG2 uses the preconditioned conjugate-gradient method to solve the equations produced by the model for hydraulic head. Linear or nonlinear flow conditions may be simulated. PCG2 includes two reconditioning options: modified incomplete Cholesky preconditioning, which is efficient on scalar computers; and polynomial preconditioning, which requires less computer storage and, with modifications that depend on the computer used, is most efficient on vector computers. Convergence of the solver is determined using both head-change and residual criteria. Nonlinear problems are solved using Picard iterations. This documentation provides a description of the preconditioned conjugate gradient method and the two preconditioners, detailed instructions for linking PCG2 to the modular model, sample data inputs, a brief description of PCG2, and a FORTRAN listing.

Water-Resources Investigations Report↗

USGS1991

Techniques and strategies for data integration in mineral resource assessment

The Geologic and the National Mapping divisions of the U.S. Geological Survey have been involved formally in cooperative research and development of computer-based geographic information systems (GISs) applied to mineral-resource assessment objectives since 1982. Experience in the Conterminous United States Mineral Assessment Program (CUSMAP) projects including the Rolla, Missouri; Dillon, Montana; Butte, Montana; and Tonopah, Nevada 1?? ?? 2?? quadrangles, has resulted in the definition of processing requirements for geographically and mineral-resource data that are common to these studies. The diverse formats of data sets collected and compiled for regional mineral-resource assessments necessitate capabilities for digitally encoding and entering data into appropriate tabular, vector, and raster subsystems of the GIS. Although many of the required data sets are either available or can be provided in a digital format suitable for direct entry, their utility is largely dependent on the original intent and consequent preprocessing of the data. In this respect, special care must be taken to ensure the digital data type, encoding, and format will meet assessment objectives. Data processing within the GIS is directed primarily toward the development and application of models that can be used to describe spatially geological, geophysical, and geochemical environments either known or inferred to be associated with specific types of mineral deposits. Consequently, capabilities to analyze spatially, aggregate, and display relations between data sets are principal processing requirements. To facilitate the development of these models within the GIS, interfaces must be developed among vector-, raster-, and tabular-based processing subsystems to reformat resident data sets for comparative analyses and multivariate display of relations.

Conference Paper↗

USGS2012

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

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

Delaware, Maryland, New York, Pennsylvania, Virgin↗

USGS2004

SAM 2.1—A computer program for plotting and formatting surveying data for estimating peak discharges by the slope-area method

The U.S. Geological Survey (USGS) measures discharge in streams using several methods. However, measurement of peak discharges is often impossible or impractical due to difficult access, inherent danger of making measurements during flood events, and timing often associated with flood events. Thus, many peak discharge values often are calculated after the fact by use of indirect methods. The most common indirect method for estimating peak dis- charges in streams is the slope-area method. This, like other indirect methods, requires measuring the flood profile through detailed surveys. Processing the survey data for efficient entry into computer streamflow models can be time demanding; SAM 2.1 is a program designed to expedite that process. The SAM 2.1 computer program is designed to be run in the field on a portable computer. The program processes digital surveying data obtained from an electronic surveying instrument during slope- area measurements. After all measurements have been completed, the program generates files to be input into the SAC (Slope-Area Computation program; Fulford, 1994) or HEC-RAS (Hydrologic Engineering Center-River Analysis System; Brunner, 2001) computer streamflow models so that an estimate of the peak discharge can be calculated.

Fact Sheet↗

USGS1993

Plans for a sensitivity analysis of bridge-scour computations

Plans for an analysis of the sensitivity of Level 2 bridge-scour computations are described. Cross-section data from 15 bridge sites in Texas are modified to reflect four levels of field effort ranging from no field surveys to complete surveys. Data from United States Geological Survey (USGS) topographic maps will be used to supplement incomplete field surveys. The cross sections are used to compute the water-surface profile through each bridge for several T-year recurrence-interval design discharges. The effect of determining the downstream energy grade-line slope from topographic maps is investigated by systematically varying the starting slope of each profile. The water-surface profile analyses are then used to compute potential scour resulting from each of the design discharges. The planned results will be presented in the form of exceedance-probability versus scour-depth plots with the maximum and minimum scour depths at each T-year discharge presented as error bars.

Conference Paper↗

USGS2013

Geologic map of Chickasaw National Recreation Area, Murray County, Oklahoma

This 1:24,000-scale geologic map is a compilation of previous geologic maps and new geologic mapping of areas in and around Chickasaw National Recreation Area. The geologic map includes revisions of numerous unit contacts and faults and a number of previously “undifferentiated” rock units were subdivided in some areas. Numerous circular-shaped hills in and around Chickasaw National Recreation Area are probably the result of karst-related collapse and may represent the erosional remnants of large, exhumed sinkholes. Geospatial registration of existing, smaller scale (1:72,000- and 1:100,000-scale) geologic maps of the area and construction of an accurate Geographic Information System (GIS) database preceded 2 years of fieldwork wherein previously mapped geology (unit contacts and faults) was verified and new geologic mapping was carried out. The geologic map of Chickasaw National Recreation Area and this pamphlet include information pertaining to how the geologic units and structural features in the map area relate to the formation of the northern Arbuckle Mountains and its Arbuckle-Simpson aquifer. The development of an accurate geospatial GIS database and the use of a handheld computer in the field greatly increased both the accuracy and efficiency in producing the 1:24,000-scale geologic map.

Oklahoma↗

USGS1997

GeoChange Global Change Data

GeoChange is an online data system providing access to research results and data generated by the U.S. Geological Survey's Global Change Research Program. Researchers in this program study climate history and the causes of climatic variations, as well as providing baseline data sets on contemporary atmospheric chemistry, high-resolution meteorology, and dust deposition. Research results are packaged as data sets, groups of digital files containing scientific observations, documentation, and interpretation. The data sets are arranged in a consistent manner using standard file formats so that users of a variety of computer systems can access and use them.

Fact Sheet↗

USGS1984

Geohydrologic conditions in the Vekol Valley, Pinal and Maricopa Counties, Arizona-1983

Vekol Valley is a north-trending valley in south-central Arizona. The valley includes 270 mi 2 and is bounded by the Table Top and Vekol Mountains on the east, the Sand Tank and Maricopa Mountains and the White Hills on the west, and the Both and Haley Hills on the north. The climate is semiarid and precipitation occurs as local summer thunderstorms and regional winter storms. The average annual precipitation is 7 to 8 in., and the potential evaporation rate is about 10 times that of the average annual rainfall (Sells and Hill, 1974, p. 7, 44). Storm runoff occurs mainly as sheetflow and floods of short duration. Although some runoff is diverted into catchment tanks for use by livestock, runoff is not used for irrigation or public supply. Vekol Valley is underlain by a thick sequence of basin sediments, sedimentary rocks, and volcanic rocks. The surrounding mountains are composed of crystalline and volcanic rocks, which also underlie the basin sediments to form the sides and bottom of the basin. The valley is divided into a northern part and a southern part by a bedrock outcrop of crystalline rocks between the White Hills on the west and the Table Top Mountains on the east (R. P. Wilson, 1979, p. 8). The northern part is about 90 mi 2 (J. R. Marie and K. J. Hollett, U.S. Geological Survey, written commun., 1984) and the southern part is about 180 mi 2 . The hydrologic data on which this map is based are available, for the most part, in computer-printout form and may be consulted at the Arizona Department of Water Resources, 99 East Virginia, Phoenix, and the U.S. Geological Survey offices in: Federal Building, 301 West Congress Street, Tucson , and 3788 North 16th Street, Suite E, Phoenix. Material from which copies can be made at private expense is available at the Tuscon and Phoenix offices of the U.S. Geological Survey.

Arizona↗

USGS2022

Dynamic rating method for computing discharge from time-series stage data

Ratings are used for a variety of reasons in water-resources investigations. The simplest rating relates discharge to the stage of the river. From a pure hydrodynamics perspective, all rivers and streams have some form of hysteresis in the relation between stage and discharge because of unsteady flow as a flood wave passes. Simple ratings are unable to represent hysteresis in a stage/discharge relation. A dynamic rating method is capable of capturing hysteresis owing to the variable energy slope caused by unsteady momentum and pressure. A dynamic rating method developed to compute discharge from stage for compact channel geometry, referred to as DYNMOD, previously has been developed through a simplification of the one-dimensional Saint-Venant equations. A dynamic rating method, which accommodates compound and compact channel geometry, referred to as DYNPOUND, has been developed through a similar simplification as a part of this study. The DYNMOD and DYNPOUND methods were implemented in the Python programming language. Discharge time series computed with the dynamic rating method implementations were then compared to simulated discharge time series and discrete discharge measurements made at U.S. Geological Survey streamgage sites. Four sets of stage and discharge time series were created using one-dimensional unsteady simulation software with compound channel geometry to compare the results of both dynamic rating methods to results from the full one-dimensional shallow water equations. Discharge time series were computed from stage time series using DYNMOD and DYNPOUND. DYNPOUND outperformed DYNMOD in all four scenarios. The minimum and maximum mean squared logarithmic error (MSLE) for the DYNMOD results were 2.75×10 −2 and 3.40×10 −2 , respectively. The minimum and maximum MSLE for the DYNPOUND results were 2.51×10 −7 and 1.91×10 −4 , respectively. The dynamic rating methods were calibrated for six U.S. Geological Survey streamgage sites using observed discharge data collected at the sites. The calibration objective for each site was to minimize the MSLE of the discharge computed with the rating method with respect to observed discharge. For each site, the calibration included all field measurements within a selected water year. The DYNMOD method failed to compute discharge for the full calibration time series for three sites. A method fails to compute when the implementation returns a nonfinite value at a time step. Because the values computed for following time steps are dependent on the previous time step, a nonfinite value results in nonfinite values that follow. For the three sites for which DYNMOD computed the complete discharge time series, the minimum MSLE for calibration was 2.19×10 −3 and the maximum was 9.77×10 −3 . The MSLE of the DYNPOUND computed discharge calibration time series for the six sites ranged from 3.70×10 −3 to 1.25. For each site, an event-based time period was selected to compare the discharge time series computed with the dynamic rating methods to discrete discharge field measurements made at the streamgage sites. The DYNMOD-computed discharge time series for the three sites had an MSLE range of 2.76×10 −3 to 3.14×10 −2 . The range of MSLE for the six DYNPOUND sites was 3.64×10 −3 to 7.23×10 −2 . Although the DYNMOD method outperforms the DYNPOUND method when calibrated streamgage sites are under consideration, the DYNMOD method failed to compute a discharge time series at three of the six sites. The DYNPOUND method, therefore, was more robust than the DYNMOD method. Improvements to the implementation of the DYNPOUND method may improve the accuracy of the method.

Open-File Report↗

USGS2010

Low-flow frequency and flow duration of selected South Carolina streams in the Broad River basin through March 2008

In 2008, the U.S. Geological Survey, in cooperation with the South Carolina Department of Health and Environmental Control, initiated a study to update low-flow statistics at continuous-record streamgaging stations operated by the U.S. Geological Survey in South Carolina. This report presents the low-flow statistics for 23 selected streamgaging stations in the Broad River basin in South Carolina, and includes flow durations of 5-, 10-, 25-, 50-, 75-, 90-, and 95-percent probability of exceedance and the annual minimum 1-, 3-, 7-, 14-, 30-, 60-, and 90-day mean flows with recurrence intervals of 2, 5, 10, 20, 30, and 50 years, depending on the length of record available at the streamgaging station. The low-flow statistics were computed from records available through March 31, 2008. In addition, flow duration information is presented for one streamgaging station 021556525, Pacolet River below Lake Blalock near Cowpens, SC, where recurrence interval computations were not appropriate. Of the 23 streamgaging stations for which recurrence interval computations were made, 14 had low-flow statistics that were published in previous U.S. Geological Survey reports. A comparison of the low-flow statistics for the minimum mean flow for a 7-consecutive-day period with a 10-year recurrence interval (7Q10) from this study with the most recently published values indicated that 8 of the 14 streamgaging stations had values that were within plus or minus 25 percent of the previous value. Ten of the 14 streamgaging stations had negative percent differences indicating the low-flow statistic had decreased since the previous study, and 4 streamgaging stations had positive percent differences indicating that the low-flow statistic had increased since the previous study. The low-flow statistics are influenced by length of record, hydrologic regime under which the record was collected, techniques used to do the analysis, and other changes, such as urbanization, diversions, and so on, that may have occurred in the basin.

South Carolina↗

USGS2010

A geostatistical approach to mapping site response spectral amplifications

If quantitative estimates of the seismic properties do not exist at a location of interest then the site response spectral amplifications must be estimated from data collected at other locations. Currently, the most common approach employs correlations of site class with maps of surficial geology. Analogously, correlations of site class with topographic slope can be employed where the surficial geology is unknown. Our goal is to identify and validate a method to estimate site response with greater spatial resolution and accuracy for regions where additional effort is warranted. This method consists of three components: region-specific data collection, a spatial model for interpolating seismic properties, and a theoretical method for computing spectral amplifications from the interpolated seismic properties. We consider three spatial interpolation schemes: correlations with surficial geology, termed the geologic trend (GT), ordinary kriging (OK), and kriging with a trend (KT). We estimate the spectral amplifications from seismic properties using the square root of impedance method, thereby linking the frequency-dependent spectral amplifications to the depth-dependent seismic properties. Thus, the range of periods for which this method is applicable is limited by the depth of exploration. A dense survey of near-surface S-wave slowness (Ss) throughout Kobe, Japan shows that the geostatistical methods give more accurate estimates of Ss than the topographic slope and GT methods, and the OK and KT methods perform equally well. We prefer the KT model because it can be seamlessly integrated with geologic maps that cover larger regions. Empirical spectral amplifications show that the region-specific data achieve more accurate estimates of observed median short-period amplifications than the topographic slope method. ?? 2010 Elsevier B.V.

Engineering Geology↗