USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computers & Operations Research”

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 73 records · Page 4Linked to original sources

An expert system for prediction of aquatic toxicity of contaminants

The National Fisheries Research Center-Great Lakes has developed an interactive computer program in muLISP that runs on an IBM-compatible microcomputer and uses a linear solvation energy relationship (LSER) to predict acute toxicity to four representative aquatic species from the detailed structure of an organic molecule. Using the SMILES formalism for a chemical structure, the expert system identifies all structural components and uses a knowledge base of rules based on an LSER to generate four structure-related parameter values. A separate module then relates these values to toxicity. The system is designed for rapid screening of potential chemical hazards before laboratory or field investigations are conducted and can be operated by users with little toxicological background. This is the first expert system based on LSER, relying on the first comprehensive compilation of rules and values for the estimation of LSER parameters.

Book chapter↗

Metrics for comparing dynamic earthquake rupture simulations

Earthquakes are complex events that involve a myriad of interactions among multiple geologic features and processes. One of the tools that is available to assist with their study is computer simulation, particularly dynamic rupture simulation. A dynamic rupture simulation is a numerical model of the physical processes that occur during an earthquake. Starting with the fault geometry, friction constitutive law, initial stress conditions, and assumptions about the condition and response of the near‐fault rocks, a dynamic earthquake rupture simulation calculates the evolution of fault slip and stress over time as part of the elastodynamic numerical solution (Ⓔ see the simulation description in the electronic supplement to this article). The complexity of the computations in a dynamic rupture simulation make it challenging to verify that the computer code is operating as intended, because there are no exact analytic solutions against which these codes’ results can be directly compared. One approach for checking if dynamic rupture computer codes are working satisfactorily is to compare each code’s results with the results of other dynamic rupture codes running the same earthquake simulation benchmark. To perform such a comparison consistently, it is necessary to have quantitative metrics. In this paper, we present a new method for quantitatively comparing the results of dynamic earthquake rupture computer simulation codes.

Seismological Research Letters↗

A new towed platform for the unobtrusive surveying of benthic habitats and organisms

Maps of coral ecosystems are needed to support many conservation and management objectives, as well as research activities. Examples include ground-truthing aerial and satellite imagery, characterizing essential habitat, assessing changes, and monitoring the progress of restoration efforts. To address some of these needs, the U.S. Geological Survey developed the Along-Track Reef-Imaging System (ATRIS), a boat-based sensor package for mapping shallow-water benthic environments. ATRIS consists of a digital still camera, a video camera, and an acoustic depth sounder affixed to a moveable pole. This design, however, restricts its deployment to clear waters less than 10 m deep. To overcome this limitation, a towed version has been developed, referred to as Deep ATRIS. The system is based on a light-weight, computer-controlled, towed vehicle that is capable of following a programmed diving profile. The vehicle is 1.3 m long with a 63-cm wing span and can carry a wide variety of research instruments, including CTDs, fluorometers, transmissometers, and cameras. Deep ATRIS is currently equipped with a high-speed (20 frames · s-1) digital camera, custom-built light-emitting-diode lights, a compass, a 3-axis orientation sensor, and a nadir-looking altimeter. The vehicle dynamically adjusts its altitude to maintain a fixed height above the seafloor. The camera has a 29° x 22° field-of-view and captures color images that are 1360 x 1024 pixels in size. GPS coordinates are recorded for each image. A gigabit ethernet connection enables the images to be displayed and archived in real time on the surface computer. Deep ATRIS has a maximum tow speed of 2.6 m · s-1and a theoretical operating tow-depth limit of 27 m. With an improved tow cable, the operating depth can be extended to 90 m. Here, we present results from the initial sea trials in the Gulf of Mexico and Biscayne National Park, Florida, USA, and discuss the utility of Deep ATRIS for map-ping coral reef habitats. Several example mosaics illustrate the high-quality imagery that can be obtained with this system. The images also reveal the potential for unobtrusive animal observations; fish and sea turtles are unperturbed by the presence of Deep ATRIS

Florida↗

Overview of the land analysis system (LAS)

The Land Analysis System (LAS) is a fully integrated digital analysis system designed to support remote sensing, image processing, and geographic information systems research. LAS is being developed through a cooperative effort between the National Aeronautics and Space Administration Goddard Space Flight Center and the U. S. Geological Survey Earth Resources Observation Systems (EROS) Data Center. LAS has over 275 analysis modules capable to performing input and output, radiometric correction, geometric registration, signal processing, logical operations, data transformation, classification, spatial analysis, nominal filtering, conversion between raster and vector data types, and display manipulation of image and ancillary data. LAS is currently implant using the Transportable Applications Executive (TAE). While TAE was designed primarily to be transportable, it still provides the necessary components for a standard user interface, terminal handling, input and output services, display management, and intersystem communications. With TAE the analyst uses the same interface to the processing modules regardless of the host computer or operating system. LAS was originally implemented at EROS on a Digital Equipment Corporation computer system under the Virtual Memorial System operating system with DeAnza displays and is presently being converted to run on a Gould Power Node and Sun workstation under the Berkeley System Distribution UNIX operating system.

Conference Paper↗

Near-field remote sensing of Alaskan Rivers

The U.S. Geological Survey (USGS) Geomorphology and Sediment Transport Laboratory (GSTL), in collaboration with the U.S. Army Corps of Engineers Cold Regions Research and Engineering Laboratory (CRREL), acquired remotely sensed data from several Alaskan rivers in 2017 and 2018 with the goal of developing a methodology for measuring streamflow from a helicopter. CRREL operates a custom airborne lidar system that can be deployed in a helicopter-based pod (HeliPod). Data were collected with the HeliPod near existing USGS streamflow information stations on the Knik, Matanuska, Chena, and Salcha Rivers in both 2017 and 2018. Sites on the Tanana and Snow Rivers were added in 2018. In 2018, the HeliPod was modified to accommodate both a thermal infrared and a visible camera. The cameras were integrated with the flight management software to simultaneously acquire imagery with lidar. The Global Navigation Satellite System (GNSS) and inertial measurement unit (IMU) in the HeliPod were used to compute trajectories with precise position and orientation information needed for image orthorectification. The HeliPod sensors provide data for measuring river channel characteristics. Lidar can map the elevation of the water surface and thus be used to measure water-surface slopes and return intensity can be used to delineate the extent of the wetted river channel. Various approaches are currently being evaluated to estimate surface flow velocity from visible and thermal image time series. In this paper, we examine and compare water-surface elevation returns and slopes derived from the HeliPod lidar and found good agreement with measurements made using conventional field-based techniques.

Alaska↗

Computer Program for Point Location And Calculation of ERror (PLACER)

A program designed for point location and calculation of error (PLACER) was developed as part of the Quality Assurance Program of the Federal Highway Administration/U.S. Geological Survey (USGS) National Data and Methodology Synthesis (NDAMS) review process. The program provides a standard method to derive study-site locations from site maps in highwayrunoff, urban-runoff, and other research reports. This report provides a guide for using PLACER, documents methods used to estimate study-site locations, documents the NDAMS Study-Site Locator Form, and documents the FORTRAN code used to implement the method. PLACER is a simple program that calculates the latitude and longitude coordinates of one or more study sites plotted on a published map and estimates the uncertainty of these calculated coordinates. PLACER calculates the latitude and longitude of each study site by interpolating between the coordinates of known features and the locations of study sites using any consistent, linear, user-defined coordinate system. This program will read data entered from the computer keyboard and(or) from a formatted text file, and will write the results to the computer screen and to a text file. PLACER is readily transferable to different computers and operating systems with few (if any) modifications because it is written in standard FORTRAN. PLACER can be used to calculate study site locations in latitude and longitude, using known map coordinates or features that are identifiable in geographic information data bases such as USGS Geographic Names Information System, which is available on the World Wide Web.

Open-File Report↗

Identifying and preserving high-water mark data

High-water marks provide valuable data for understanding recent and historical flood events. The proper collection and recording of high-water mark data from perishable and preserved evidence informs flood assessments, research, and water resource management. Given the high cost of flooding in developed areas, experienced hydrographers, using the best available techniques, can contribute high-quality data toward efforts such as public education of flood risk, flood inundation mapping, flood frequency computations, indirect streamflow measurement, and hazard assessments. This manual presents guidance for skilled high-water mark identification, including marks left behind in natural and man-made environments by tranquil and rapid flowing water. This manual also presents pitfalls and challenges associated with various types of flood evidence that help hydrographers identify the best high-water marks and assess the uncertainty associated with a given mark. Proficient high-water mark data collection contributes to better understanding of the flooding process and reduces risk through greater ability to estimate flood probability. The U.S. Geological Survey, operating the Nation’s premier water data collection network, encourages readers of this manual to familiarize themselves with the art and science of high-water mark collection. The U.S. Geological survey maintains a national database at http://water.usgs.gov/floods/FEV/ that includes high-water mark information for many flood events, and local U.S. Geological Survey Water Science Centers can provide information to interested readers about participation in data collection and flood documentation efforts as volunteers or observers.

Techniques and Methods↗

Comparison of methods for estimating ground-water recharge and base flow at a small watershed underlain by fractured bedrock in the Eastern United States

This study by the U.S. Geological Survey (USGS), in cooperation with the Agricultural Research Service (ARS), U.S. Department of Agriculture, compared multiple methods for estimating ground-water recharge and base flow (as a proxy for recharge) at sites in east-central Pennsylvania underlain by fractured bedrock and representative of a humid-continental climate. This study was one of several within the USGS Ground-Water Resources Program designed to provide an improved understanding of methods for estimating recharge in the eastern United States. Recharge was estimated on a monthly and annual basis using four methods?(1) unsaturated-zone drainage collected in gravity lysimeters, (2) daily water balance, (3) water-table fluctuations in wells, and (4) equations of Rorabaugh. Base flow was estimated by streamflow-hydrograph separation using the computer programs PART and HYSEP. Estimates of recharge and base flow were compared for an 8-year period (1994-2001) coinciding with operation of the gravity lysimeters at an experimental recharge site (Masser Recharge Site) and a longer 34-year period (1968-2001), for which climate and streamflow data were available on a 2.8-square-mile watershed (WE-38 watershed). Estimates of mean-annual recharge at the Masser Recharge Site and WE-38 watershed for 1994-2001 ranged from 9.9 to 14.0 inches (24 to 33 percent of precipitation). Recharge, in inches, from the various methods was: unsaturated-zone drainage, 12.2; daily water balance, 12.3; Rorabaugh equations with PULSE, 10.2, or RORA, 14.0; and water-table fluctuations, 9.9. Mean-annual base flow from streamflow-hydrograph separation ranged from 9.0 to 11.6 inches (21-28 percent of precipitation). Base flow, in inches, from the various methods was: PART, 10.7; HYSEP Local Minimum, 9.0; HYSEP Sliding Interval, 11.5; and HYSEP Fixed Interval, 11.6. Estimating recharge from multiple methods is useful, but the inherent differences of the methods must be considered when comparing results. For example, although unsaturated-zone drainage from the gravity lysimeters provided the most direct measure of potential recharge, it does not incorporate spatial variability that is contained in watershed-wide estimates of net recharge from the Rorabaugh equations or base flow from streamflow-hydrograph separation. This study showed that water-level fluctuations, in particular, should be used with caution to estimate recharge in low-storage fractured-rock aquifers because of the variability of water-level response among wells and sensitivity of recharge to small errors in estimating specific yield. To bracket the largest range of plausible recharge, results from this study indicate that recharge derived from RORA should be compared with base flow from the Local-Minimum version of HYSEP.

Scientific Investigations Report↗

Digital Twin Earth - Coasts: Developing a fast and physics-informed surrogate model for coastal floods via neural operators

Developing fast and accurate surrogates for physics-based coastal and ocean mod- els is an urgent need due to the coastal flood risk under accelerating sea level rise, and the computational expense of deterministic numerical models. For this purpose, we develop the first digital twin of Earth coastlines with new physics-informed machine learning techniques extending the state-of-art Neural Operator. As a proof-of-concept study, we built Fourier Neural Operator (FNO) surrogates on the simulations of an industry-standard coastal and ocean model – Nucleus for Euro- pean Modelling of the Ocean (NEMO). The resulting FNO surrogate accurately predicts the sea surface height in most regions while achieving upwards of 45x acceleration of NEMO. We delivered an open-source CoastalTwin platform in an end-to-end and modular way, to enable easy extensions to other simulations and ML-based surrogate methods. Our results and deliverable provide a promising approach to massively accelerate coastal dynamics simulators, which can enable scientists to efficiently execute many simulations for decision-making, uncertainty quantification, and other research activities.

Conference Paper↗

Evaluation of the fish passage effectiveness of the Bonneville I prototype surface collector using three-dimensional ultrasonic fish tracking - Final Report

This report describes tests conducted at Bonneville Dam on the Columbia River in the spring of 2000. The studies used three-dimensional (3D) acoustic telemetry and computational fluid dynamics (CFD) hydraulic modeling techniques to evaluate the response of outmigrating juvenile steelhead (Oncorhynchus mykiss) and yearling chinook (O. tshawytscha) to the Prototype Surface Collector (PSC) installed at Powerhouse I of Bonneville Dam in 1998 to test the concept of using a deep-slot surface bypass collector to divert downstream migrating salmon from turbines. The study was conducted by Pacific Northwest National Laboratory (PNNL), the Waterways Experiment Station of the U.S. Army Corp of Engineers (COE), Asci Corporation, and the U.S. Geological Survey (USGS), and was sponsored by COE’s Portland District. The goal of the study was to observe the three-dimensional behavior of tagged fish (fish bearing ultrasonic micro-transmitters) within 100 meters (m) of the surface flow bypass structure to test hypotheses about the response of migrants to flow stimuli generated by the presence of the surface flow bypass prototype and its operation. Research was done in parallel with radio telemetry studies conducted by USGS and hydroacoustic studies conducted by WES & Asci to evaluate the prototype surface collector.

Report↗

Software for analysis of chemical mixtures--composition, occurrence, distribution, and possible toxicity

The composition, occurrence, distribution, and possible toxicity of chemical mixtures in the environment are research concerns of the U.S. Geological Survey and others. The presence of specific chemical mixtures may serve as indicators of natural phenomena or human-caused events. Chemical mixtures may also have ecological, industrial, geochemical, or toxicological effects. Chemical-mixture occurrences vary by analyte composition and concentration. Four related computer programs have been developed by the National Water-Quality Assessment Program of the U.S. Geological Survey for research of chemical-mixture compositions, occurrences, distributions, and possible toxicities. The compositions and occurrences are identified for the user-supplied data, and therefore the resultant counts are constrained by the user’s choices for the selection of chemicals, reporting limits for the analytical methods, spatial coverage, and time span for the data supplied. The distribution of chemical mixtures may be spatial, temporal, and (or) related to some other variable, such as chemical usage. Possible toxicities optionally are estimated from user-supplied benchmark data. The software for the analysis of chemical mixtures described in this report is designed to work with chemical-analysis data files retrieved from the U.S. Geological Survey National Water Information System but can also be used with appropriately formatted data from other sources. Installation and usage of the mixture software are documented. This mixture software was designed to function with minimal changes on a variety of computer-operating systems. To obtain the software described herein and other U.S. Geological Survey software, visit http://water.usgs.gov/software/.

Scientific Investigations Report↗

The proposed EROSpace institute, a national center operated by space grant universities

The "EROSpace Institute" is a proposed visiting scientist program in associated with the U.S. Geological Survey's EROS Data Center (EDC). The Institute would be operated by a consortium of universities, possible drawn from NASA's Space Grant College and Fellowship Program consortia and the group of 17 capability-enhancement consortia, or perhaps from consortia though out the nation with a topical interest in remote sensing. The National Center for Atmospheric Research or the Goddard Institute for Space Studies provide models for the structure of such an institute. The objectives of the Institute are to provide ready access to the body of data housed at the EDC and to increase the cadre of knowledgeable and trained scientists able to deal with the increasing volume of remote sensing data to become available from the Earth Observing System. The Institute would have a staff of about 100 scientists at any one time, about half permanent staff, and half visiting scientists. The latter would include graduate and undergraduate students, as well as faculty on temporary visits, summer fellowships, or sabbatical leaves. The Institute would provide office and computing facilities, as well as Internet linkages to the home institutions so that scientists could continue to participate in the program from their home base.

Pecora 12 Symposium↗

Operational earthquake forecasting during the 2019 Ridgecrest, California, earthquake sequence with the UCERF3-ETAS model

The first Uniform California Earthquake Rupture Forecast, Version 3–epidemic‐type aftershock sequence (UCERF3‐ETAS) aftershock simulations were running on a high‐performance computing cluster within 33 min of the 4 July 2019 M 6.4 Searles Valley earthquake. UCERF3‐ETAS, an extension of the third Uniform California Earthquake Rupture Forecast (UCERF3), is the first comprehensive, fault‐based, epidemic‐type aftershock sequence (ETAS) model. It produces ensembles of synthetic aftershock sequences both on and off explicitly modeled UCERF3 faults to answer a key question repeatedly asked during the Ridgecrest sequence: What are the chances that the earthquake that just occurred will turn out to be the foreshock of an even bigger event?

California↗

IRIS/USGS plans for upgrading the Global Seismograph Network

This report has been prepared to provide information to organizations that may be asked to participate in a program to upgrade the global seismographic network. In most cases, the organizations that will be offered new instrumentation by the U.S. Geological Survey currently operate stations in the World-Wide Standardized Seismograph Network (WWSSN) or the Global Digital Seismograph Network (GDSN). The deployment of the WWSSN in the 1960's and the subsequent equipping of some WWSSN stations with digital equipment and borehole seismometers during the 1970's has been a remarkably successful program that generated the high- quality data needed to fuel an unprecedented period of progress in earthquake and tectonic research. The success of the WWSSN can be attributed to the importance of the data, to the strong commitment by participating organizations to international scientific cooperation, to the dedication and skill of the station operators, and to the resourcefulness of the staff supporting the network. Benefits have been widespread. The community of scientists world-wide has benefited from unrestricted access to a standardized base of calibrated data, and the participating stations have benefited from the donation of modern observatory instruments that have been useful for local earthquake studies and for the training of scientists and engineers. Now, an exciting opportunity has arisen to deploy a new generation of seismograph systems to replace the outdated equipment at many of the WWSSN and GDSN stations. The U.S. Geological Survey (USGS) is cooperating with the Incorporated Research Institutions for Seismology (IRIS) in a program to upgrade the global seismograph network. The equipment development phase is nearly complete with a prototype of the new broadband seismograph system currently undergoing final testing at the USGS Albuquerque Seismological Laboratory. Deployment of the new equipment is expected to begin in early 1990. As this report will demonstrate, the IRIS broadband seismograph system combines the very latest data acquisition and computer technology to produce seismic data with unprecedented bandwidth and dynamic range. Moreover, the system has been designed so that the high-quality digital data are accessible for local display and analysis. The functional design of the new system, which uses off-the-shelf modules and a standard computer bus, will make it much easier than it has been in the past to modify and upgrade the data acquisition system as improvements in technology become available. With adequate support for the program, the new IRIS seismograph system need never become obsolete. We want you to be aware of our plans and the possibility that you may be asked to participate in this ,program. The schedule for upgrading WWSSN and GDSN stations depends on the level of funding earmarked for the program by our National Science Foundation. We hope to deploy at least ten new GSN data systems each year. If you have any questions concerning this program, please contact the Albuquerque Seismological Laboratory, U.S. Geological Survey, Albuquerque, New Mexico 87115-5000. This report was revised in February 1992 in order to update information concerning the current program and instrumentation. The amp in Figure 1 was revised in June 1993, April 1994, December 1994, and September 1996 to reflect updated siting information. In September 1996 a composite photo of standard and optional components of the IRIS-2 GSN system hardware was added as a separate page between Figures 9 and 10.

Open-File Report↗

A General Lake Model (GLM 3.0) for linking with high-frequency sensor data from the Global Lake Ecological Observatory Network (GLEON)

The General Lake Model (GLM) is a one-dimensional open-source code designed to simulate the hydrodynamics of lakes, reservoirs, and wetlands. GLM was developed to support the science needs of the Global Lake Ecological Observatory Network (GLEON), a network of researchers using sensors to understand lake functioning and address questions about how lakes around the world respond to climate and land use change. The scale and diversity of lake types, locations, and sizes, and the expanding observational datasets created the need for a robust community model of lake dynamics with sufficient flexibility to accommodate a range of scientific and management questions relevant to the GLEON community. This paper summarizes the scientific basis and numerical implementation of the model algorithms, including details of sub-models that simulate surface heat exchange and ice cover dynamics, vertical mixing, and inflow–outflow dynamics. We demonstrate the suitability of the model for different lake types that vary substantially in their morphology, hydrology, and climatic conditions. GLM supports a dynamic coupling with biogeochemical and ecological modelling libraries for integrated simulations of water quality and ecosystem health, and options for integration with other environmental models are outlined. Finally, we discuss utilities for the analysis of model outputs and uncertainty assessments, model operation within a distributed cloud-computing environment, and as a tool to support the learning of network participants.

Geoscientific Model Development↗

Synoptic evaluation of scale-dependent metrics for hydrographic line feature geometry

Methods of acquisition and feature simplification for vector feature data impact cartographic representations and scientific investigations of these data, and are therefore important considerations for geographic information science (Haunert and Sester 2008). After initial collection, linear features may be simplified to reduce excessive detail or to furnish a reduced-scale version of the features through cartographic generalization (Regnauld and McMaster 2008, Stanislawski et al. 2014). A variety of algorithms exist to simplify linear cartographic features, and all of the methods affect the positional accuracy of the features (Shahriari and Tao 2002, Regnauld and McMaster 2008, Stanislawski et al. 2012). In general, simplification operations are controlled by one or more tolerance parameters that limit the amount of positional change the operation can make to features. Using a single tolerance value can have varying levels of positional change on features; depending on local shape, texture, or geometric characteristics of the original features (McMaster and Shea 1992, Shahriari and Tao 2002, Buttenfield et al. 2010). Consequently, numerous researchers have advocated calibration of simplification parameters to control quantifiable properties of resulting changes to the features (Li and Openshaw 1990, Raposo 2013, Tobler 1988, Veregin 2000, and Buttenfield, 1986, 1989). This research identifies relations between local topographic conditions and geometric characteristics of linear features that are available in the National Hydrography Dataset (NHD). The NHD is a comprehensive vector dataset of surface 18 th ICA Workshop on Generalisation and Multiple Representation, Rio de Janiero, Brazil 2015 2 water features within the United States that is maintained by the U.S. Geological Survey (USGS). In this paper, geometric characteristics of cartographic representations for natural stream and river features are summarized for subbasin watersheds within entire regions of the conterminous United States and compared to topographic metrics. A concurrent processing workflow is implemented using a Linux high-performance computing cluster to simultaneously process multiple subbasins, and thereby complete the work in a fraction of the time required for a single-process environment. In addition, similar metrics are generated for several levels of simplification of the hydrographic features to quantify the effects of simplification over the various landscape conditions. Objectives of this exploratory investigation are to quantify geometric characteristics of linear hydrographic features over the various terrain conditions within the conterminous United States and thereby illuminate relations between stream geomorphological conditions and cartographic representation. The synoptic view of these characteristics over regional watersheds that is afforded through concurrent processing, in conjunction with terrain conditions, may reveal patterns for classifying cartographic stream features into stream geomorphological classes. Furthermore, the synoptic measurement of the amount of change in geometric characteristics caused by the several levels of simplification can enable estimation of tolerance values that appropriately control simplification-induced geometric change of the cartographic features within the various geomorphological classes in the country. Hence, these empirically derived rules or relations could help generate multiscale-representations of features through automated generalization that adequately maintain surface drainage variations and patterns reflective of the natural stream geomorphological conditions across the country.

Conference Paper↗

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↗

Evaluating the applicability of the generalized power-law rating curve model: With applications to paired discharge-stage data from Iceland, Sweden, and the United States

Hydrologic research and operations make extensive use of streamflow time series. In most applications, these time series are estimated from rating curves, which relate flow to some easy-to-measure surrogate, typically stage. The conventional stage-discharge rating takes the form of a segmented power law, with one segment for each hydrologic control at the stream gauge. However, these ratings are notoriously difficult to estimate with numerical methods, so that most are still developed manually. A few automated algorithms have emerged, but their use is sporadic, and their relative merits have not been rigorously assessed. One recently developed approach, the generalized power-law, avoids the segmenting problem by representing the power-law exponent as a Gaussian process. On the one hand, this representation is more flexible and easier to fit, but its flexibility might allow unrealistic solutions, so it needs to be tested under a range of conditions to assess its operational viability. This study evaluates the generalized power-law rating curve model by applying it to observations from 180 streams in Iceland, Sweden, and the United States. Overall, the model proved flexible and computationally robust, generating convincing rating curves across a range of geographic settings and was comparable to curves generated by a segmented rating model. Lastly, we propose a model-selection algorithm based on information theory to help identify the best rating curve model for a particular stream gauge.

Journal of Hydrology↗