Geology Reports⌕ Search

SEARCH · Geology Reports

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 91 records · Page 5Linked to original sources

A method to establish seismic noise baselines for automated station assessment

We present a method for quantifying station noise baselines and characterizing the spectral shape of out-of-nominal noise sources. Our intent is to automate this method in order to ensure that only the highest-quality data are used in rapid earthquake products at NEIC. In addition, the station noise baselines provide a valuable tool to support the quality control of GSN and ANSS backbone data and metadata. The procedures addressed here are currently in development at the NEIC, and work is underway to understand how quickly changes from nominal can be observed and used within the NEIC processing framework. The spectral methods and software used to compute station baselines and described herein (PQLX) can be useful to both permanent and portable seismic stations operators. Applications include: general seismic station and data quality control (QC), evaluation of instrument responses, assessment of near real-time communication system performance, characterization of site cultural noise conditions, and evaluation of sensor vault design, as well as assessment of gross network capabilities (McNamara et al. 2005). Future PQLX development plans include incorporating station baselines for automated QC methods and automating station status report generation and notification based on user-defined QC parameters. The PQLX software is available through the USGS (http://earthquake. usgs.gov/research/software/pqlx.php) and IRIS (http://www.iris.edu/software/ pqlx/).

Seismological Research Letters↗

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↗

Origins, evolutions, and future directions of Landsat science products for advancing global inland water and coastal ocean observations

In April 2020, the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center introduced a Level 2 provisional Aquatic Reflectance (AR) product for the Landsat 8 Operational Land Imager (OLI), marking the initial phase in developing a standardized global product for Landsat-derived surface water measurements. The goal of USGS EROS aquatic product research and development is to prepare for an operational processing architecture for Landsat Collection 3 in the late 2020s that will enable use of quality-controlled data for emerging Landsat aquatic science applications. To achieve this, we released a subset of the Landsat 8/9 provisional AR products (Crawford et al., 2025, https://doi.org/10.5066/P14MBBRM) and examined its general performance through the Science Algorithms to Operations (SATO) framework alongside quantitative assessment using community made inland water data records (GLObal Reflectance community dataset for Imaging and optical sensing of Aquatic environments, GLORIA) and radiometric coastal validation platforms (NASA’s Ocean Color component of the Aerosol Robotic Network, AERONET-OC). Variability within the validation datasets indicate that the performance of the Landsat 8/9 provisional AR retrieval is highly context-dependent; errors are minimal in optically simple waters (e.g., clear to moderately turbid coastal waters) but increase considerably in optically complex waters where factors such as elevated levels of turbidity, chlorophyll (Chl a ) concentrations, or colored dissolved organic matter (CDOM) dominate the water column. Additionally, this paper examines key algorithmic considerations for atmospheric correction, highlighting factors that influence accuracy, scalability, and computational efficiency necessary for collection processing in the operational Landsat Product Generation System (LPGS). This paper is intended to communicate with aquatic scientists, satellite oceanographers, and the broader Earth observation community on the origins, requirements, challenges, successes, and future objectives for operationalizing global AR data products for Landsat satellite missions.

Earth System Science Data↗

IVyTools: Space-Time Image Velocimetry for streamflow

Accurate streamflow measurements are essential for hydrologic monitoring, flood forecasting, and water resource management. Traditional in situ methods, although reliable, are often impractical or unsafe during flood conditions or at inaccessible locations. Image velocimetry techniques offer a potential non-contact alternative, yet operational adoption has been hindered by the lack of standardized, open-source tools, particularly for Space–Time Image Velocimetry (STIV). This paper introduces IVyTools, an open-source software application for operational streamflow measurement using STIV. Developed by the US Geological Survey (USGS), IVyTools integrates standardized streamflow measurement protocols into a reproducible, quality-controlled workflow that includes video preprocessing, orthorectification, velocity estimation, and streamflow (discharge) computation. The software incorporates a dual-method uncertainty framework based on ISO 748:2007 and an adapted Interpolated Variance Estimator (IVE), enabling users to quantify and diagnose measurement uncertainty. Validation against an independently published benchmark dataset demonstrates that IVyTools achieves a mean absolute percentage error of 4.56% (MAE = 5.16 m 3 /s ) from field-measured reference discharges, with acceptable performance across diverse river conditions. This work supports the broader adoption of image-based streamflow measurement methods and provides a foundation for potential future automation and integration into real-time hydrologic networks.

River Research and Applications↗

Inference of strata separation and gas emission paths in longwall overburden using continuous wavelet transform of well logs and geostatistical simulation

Prediction of potential methane emission pathways from various sources into active mine workings or sealed gobs from longwall overburden is important for controlling methane and for improving mining safety. The aim of this paper is to infer strata separation intervals and thus gas emission pathways from standard well log data. The proposed technique was applied to well logs acquired through the Mary Lee/Blue Creek coal seam of the Upper Pottsville Formation in the Black Warrior Basin, Alabama, using well logs from a series of boreholes aligned along a nearly linear profile. For this purpose, continuous wavelet transform (CWT) of digitized gamma well logs was performed by using Mexican hat and Morlet, as the mother wavelets, to identify potential discontinuities in the signal. Pointwise Hölder exponents (PHE) of gamma logs were also computed using the generalized quadratic variations (GQV) method to identify the location and strength of singularities of well log signals as a complementary analysis. PHEs and wavelet coefficients were analyzed to find the locations of singularities along the logs. Using the well logs in this study, locations of predicted singularities were used as indicators in single normal equation simulation (SNESIM) to generate equi-probable realizations of potential strata separation intervals. Horizontal and vertical variograms of realizations were then analyzed and compared with those of indicator data and training image (TI) data using the Kruskal–Wallis test. A sum of squared differences was employed to select the most probable realization representing the locations of potential strata separations and methane flow paths. Results indicated that singularities located in well log signals reliably correlated with strata transitions or discontinuities within the strata. Geostatistical simulation of these discontinuities provided information about the location and extents of the continuous channels that may form during mining. If there is a gas source within their zone of influence, paths may develop and allow methane movement towards sealed or active gobs under pressure differentials. Knowledge gained from this research will better prepare mine operations for potential methane inflows, thus improving mine safety.

Journal of Applied Geophysics↗

Profile measurements and data from the 2011 Optics, Acoustics, and Stress In Situ (OASIS) project at the Martha's Vineyard Coastal Observatory

This report documents data collected by the U.S. Geological Survey (USGS) for the Coastal Model Applications and Field Measurements project under the auspices of the U.S. Navy Office of Naval Research Optics, Acoustics, and Stress In Situ (OASIS) Project. The objective of the measurements was to relate optical and acoustic properties of suspended particles to changes in particle size, concentration, and vertical distribution in the bottom boundary layer near the seafloor caused by wave- and current-induced stresses. This information on the physics of particle resuspension and aggregation and light penetration and water clarity will help improve models of sediment transport, benthic primary productivity, and underwater visibility. There is well-established technology for acoustic profiling, but optical profiles are more difficult to obtain because of the rapid attenuation of light in water. A specially modified tripod with a moving arm was designed to solve this problem by moving instruments vertically in the bottom boundary layer, between the bottom and about 2 meters above the seafloor. The profiling arm was designed, built, and tested during spring and summer 2011 by a team of USGS scientists, engineers, and technicians. To accommodate power requirements and the large data files recorded by some of the optical instruments, the tripod was connected via underwater cable to the Martha's Vineyard Coastal Observatory, operated by the Woods Hole Oceanographic Institution (WHOI). This afforded real-time Internet communication with the embedded computers aboard the tripod. Instruments were mounted on the profiling arm, and additional instruments were mounted elsewhere on the tripod and nearby on the seafloor. The tripod and a small mooring for a profiling current meter were deployed on September 17, 2011, at the Martha's Vineyard Coastal Observatory 12-meter-deep underwater node about 2 kilometers south of Martha's Vineyard, Massachusetts. Divers assisted in the deployment and cleaned the instrument surfaces on the tripod approximately once per week until the tripod and current meter were recovered on October 23, 2011. There was a range of wave and current conditions during the 36-day deployment, including the distant passage of Hurricane Ophelia, several moderate wave events, and a significant local gale that generated wave heights greater than 4 meters at the 12-meter site and knocked over the tripod 3 days before it was recovered. All but one of the instruments functioned well and provided complete datasets. The details of these data and the location of files containing the best basic version of the data are described in this report.

Massachusetts↗

PQLX: A seismic data quality control system description, applications, and users manual

We present a detailed description and users manual for a new tool to evaluate seismic station performance and characteristics by providing quick and easy transitions between visualizations of the frequency and time domains. The software is based on the probability density functions (PDF) of power spectral densities (PSD) (McNamara and Buland, 2004) and builds on the original development of the PDF stand-alone software system (McNamara and Boaz, 2005) and the seismological data viewer application PQL (IRIS-PASSCAL Quick Look) and PQLII (available through the IRIS PASSCAL program: http://www.passcal.nmt.edu/content/pql-ii-program-viewing-data). With PQLX (PQL eXtended), computed PSDs are stored in a MySQL database, allowing a user to access specific time periods of PSDs (PDF subsets) and time series segments through a GUI-driven interface. The power of the method and software lies in the fact that there is no need to screen the data for system transients, earthquakes, or general data artifacts, because they map into a background probability level. In fact, examination of artifacts related to station operation and episodic cultural noise allow us to estimate both the overall station quality and a baseline level of Earth noise at each site. The output of this analysis tool is useful for both operational and scientific applications. Operationally, it is useful for characterizing the current and past performance of existing broadband stations, for conducting tests on potential new seismic station locations, for evaluating station baseline noise levels (McNamara and others, 2009), for detecting problems with the recording system or sensors, and for evaluating the overall quality of data and metadata. Scientifically, the tool allows for mining of PSDs for investigations on the evolution of seismic noise (for example, Aster and others, 2008; and Aster and others, 2010) and other phenomena. Currently, PQLX is operational at several organizations including the USGS National Earthquake Information Center (NEIC), the USGS Albuquerque Seismological Laboratory (ASL), and the Incorporated Research Institutions in Seismology (IRIS) Data Management Center (DMC) for station monitoring and instrument response quality control. The PQLX system is available to the community at large through the U.S. Geological Survey (USGS) (http://ehpm-earthquake.wr.usgs.gov/research/software/pqlx.php) and IRIS (http://www.iris.edu/software/pqlx). Also provided is a fully searchable website for bug reporting and enhancement requests (http://wush.net/bugzilla/PQLX). The first part of this document aims to describe and illustrate some of the features and capabilities of the software. The second part of this document is a detailed users manual that covers installation procedures, system requirements, operations, bug reporting, and software components (Appendix).

Open-File Report↗

Low-cost computer classification of land cover in the Portland area, Oregon, by signature extension techniques

Computer-aided techniques for interpreting multispectral data acquired by Landsat offer economies in the mapping of land cover. Even so, the actual establishment of the statistical classes, or "signatures," is one of the relatively more costly operations involved. Analysts have, therefore, been seeking cost-saving signature extension techniques that would accept training data acquired at one time or place and apply them to another. Signatures may be extended in preprocessing steps or in the classification steps that follow. In the present example, land cover classes were derived by the simplest and most direct form of signature extension: Classes statistically derived from a Landsat scene of the Puget Sound area, Wash., were applied to the adjacent Landsat scene of the Portland area, Oreg., acquired during the next 25 seconds down orbit. Many features can be recognized on the reduced-scale version of the Portland land cover map shown in this report although no statistical assessment of its accuracy is available.

Oregon, Washington↗

Reference manual for generation and analysis of Habitat Time Series: version II

The selection of an instream flow requirement for water resource management often requires the review of how the physical habitat changes through time. This review is referred to as 'Time Series Analysis." The Tune Series Library (fSLIB) is a group of programs to enter, transform, analyze, and display time series data for use in stream habitat assessment. A time series may be defined as a sequence of data recorded or calculated over time. Examples might be historical monthly flow, predicted monthly weighted usable area, daily electrical power generation, annual irrigation diversion, and so forth. The time series can be analyzed, both descriptively and analytically, to understand the importance of the variation in the events over time. This is especially useful in the development of instream flow needs based on habitat availability. The TSLIB group of programs assumes that you have an adequate study plan to guide you in your analysis. You need to already have knowledge about such things as time period and time step, species and life stages to consider, and appropriate comparisons or statistics to be produced and displayed or tabulated. Knowing your destination, you must first evaluate whether TSLIB can get you there. Remember, data are not answers . This publication is a reference manual to TSLIB and is intended to be a guide to the process of using the various programs in TSLIB. This manual is essentially limited to the hands-on use of the various programs. a TSLIB use interface program (called RTSM) has been developed to provide an integrated working environment where the use has a brief on-line description of each TSLIB program with the capability to run the TSLIB program while in the user interface. For information on the RTSM program, refer to Appendix F. Before applying the computer models described herein, it is recommended that the user enroll in the short course "Problem Solving with the Instream Flow Incremental Methodology (IFIM)." This course is offered by the Aquatic Systems Branch of the National Ecology Research Center. For more information about the TSLIB software, refer to the Memorandum of Understanding. Chapter 1 provides a brief introduction to the Instream Flow Incremental Methodology and TSLIB. Other chapters in this manual provide information on the different aspects of using the models. The information contained in the other chapters includes (2) acquisition, entry, manipulation, and listing of streamflow data; (3) entry, manipulation, and listing of the habitat-versus-streamflow function; (4) transferring streamflow data; (5) water resources systems analysis; (6) generation and analysis of daily streamflow and habitat values; (7) generation of the time series of monthly habitats; (8) manipulation, analysis, and display of month time series data; and (9) generation, analysis, and display of annual time series data. Each section includes documentation for the programs therein with at least one page of information for each program, including a program description, instructions for running the program, and sample output. The Appendixes contain the following: (A) sample file formats; (B) descriptions of default filenames; (C) alphabetical summary of batch-procedure files; (D) installing and running TSLIB on a microcomputer; (E) running TSLIB on a CDC Cyber computer; (F) using the TSLIB user interface program (RTSM); and (G) running WATSTORE on the USGS Amdahl mainframe computer. The number for this version of TSLIB--Version II-- is somewhat arbitrary, as the TSLIB programs were collected into a library some time ago; but operators tended to use and manage them as individual programs. Therefore, we will consider the group of programs from the past that were only on the CDC Cyber computer as Version 0; the programs from the past that were on both the Cyber and the IBM-compatible microcomputer as Version I; and the programs contained in this reference manual as Version II.

Report↗

A history of telemetry in fishery research

Biotelemetry has been defined as “the instrumental technique for gaining and transmitting information from a living organism and its environment to a remote observer” (Slater 1965). Biotelemetry typically utilizes wireless transmission of either an audible signal or electronic data to determine location of a tagged animal. Fisheries researchers use location information to gain a variety of insights into migration, habitat use, behavior, productivity, or survival of fish. Biotelemetry can be divided into two basic categories, acoustic or radio, based on mode of transmission, mechanical or electromagnetic energy, and operating frequency. Most acoustic systems in use today transmit at low frequency, between 30 and 300 kHz, while most radio systems transmit at very high frequency, between 30 and 300 MHz (Sisak and Lotimer 1998). Acoustic telemetry is based on the principals of sonar (sound navigation and ranging), which was developed to detect submarines during World War I. The properties of acoustic systems favor their use in deep waters with high conductivity and low turbulence (Winter 1996). Radio telemetry is based on the principals of wireless radio communication, which were first demonstrated by Nikola Tesla in 1893. Radio systems are best suited in shallow waters with relatively low conductivity but have the added benefit of improved signal detection in turbulent conditions and with aerial antennas. Advances in both technologies have resulted in highly efficient transmitter and receiving systems. Advancements in products used for animal telemetry over the past 50 years have generally followed those in the electronics field (Figure 1). Bell Laboratories1 ushered in the age of digital electronics with the invention of the transistor in 1947 (Mann 2000). Today transistors are common in everyday items such as radios, televisions, hearing aids, computers, cell phones and even MP3 players. Consumer demand for inexpensive small electronic devices with increased functionality has continually driven advancements in the field of electronics. These advancements have subsequently led to improvements in biotelemetry transmitters and receivers such as miniaturization of components, increased battery performance, and more powerful micro-processing.

Book chapter↗

Comparison of fluvial suspended-sediment concentrations and particle-size distributions measured with in-stream laser diffraction and in physical samples

Laser-diffraction technology, recently adapted for in-stream measurement of fluvial suspended-sediment concentrations (SSCs) and particle-size distributions (PSDs), was tested with a streamlined (SL), isokinetic version of the Laser In-Situ Scattering and Transmissometry (LISST) for measuring volumetric SSCs and PSDs ranging from 1.8-415 µm in 32 log-spaced size classes. Measured SSCs and PSDs from the LISST-SL were compared to a suite of 22 datasets (262 samples in all) of concurrent suspended-sediment and streamflow measurements using a physical sampler and acoustic Doppler current profiler collected during 2010-12 at 16 U.S. Geological Survey streamflow-gaging stations in Illinois and Washington (basin areas: 38 – 69,264 km 2 ). An unrealistically low computed effective density (mass SSC / volumetric SSC) of 1.24 g/ml (95% confidence interval: 1.05-1.45 g/ml) provided the best-fit value (R 2 = 0.95; RMSE = 143 mg/L) for converting volumetric SSC to mass SSC for over 2 orders of magnitude of SSC (12-2,170 mg/L; covering a substantial range of SSC that can be measured by the LISST-SL) despite being substantially lower than the sediment particle density of 2.67 g/ml (range: 2.56-2.87 g/ml, 23 samples). The PSDs measured by the LISST-SL were in good agreement with those derived from physical samples over the LISST-SL's measureable size range. Technical and operational limitations of the LISST-SL are provided to facilitate the collection of more accurate data in the future. Additionally, the spatial and temporal variability of SSC and PSD measured by the LISST-SL is briefly described to motivate its potential for advancing our understanding of suspended-sediment transport by rivers.

Water Resources Research↗

Computation and analysis of the instantaneous-discharge record for the Colorado River at Lees Ferry, Arizona — May 8, 1921, through September 30, 2000

A gaging station has been operated by the U.S. Geological Survey at Lees Ferry, Arizona, since May 8, 1921. In March 1963, Glen Canyon Dam was closed 15.5 miles upstream, cutting off the upstream sediment supply and regulating the discharge of the Colorado River at Lees Ferry for the first time in history. To evaluate the pre-dam variability in the hydrology of the Colorado River, and to determine the effect of the operation of Glen Canyon Dam on the downstream hydrology of the river, a continuous record of the instantaneous discharge of the river at Lees Ferry was constructed and analyzed for the entire period of record between May 8, 1921, and September 30, 2000. This effort involved retrieval from the Federal Records Centers and then synthesis of all the raw historical data collected by the U.S. Geological Survey at Lees Ferry. As part of this process, the peak discharges of the two largest historical floods at Lees Ferry, the 1884 and 1921 floods, were reanalyzed and recomputed. This reanalysis indicates that the peak discharge of the 1884 flood was 210,000±30,000 cubic feet per second (ft 3 /s), and the peak discharge of the 1921 flood was 170,000±20,000 ft 3 /s. These values are indistinguishable from the peak discharges of these floods originally estimated or published by the U.S. Geological Survey, but are substantially less than the currently accepted peak discharges of these floods. The entire continuous record of instantaneous discharge of the Colorado River at Lees Ferry can now be requested from the U.S. Geological Survey Grand Canyon Monitoring and Research Center, Flagstaff, Arizona, and is also available electronically at http://www.gcmrc.gov. This record is perhaps the longest (almost 80 years) high-resolution (mostly 15- to 30-minute precision) times series of river discharge available. Analyses of these data, therefore, provide an unparalleled characterization of both the natural variability in the discharge of a river and the effects of dam operations on a river. Following the construction and quality-control checks of the continuous record of instantaneous discharge, analyses of flow duration, sub-daily flow variability, and flood frequency were conducted on the pre- and post-dam parts of the record. These analyses indicate that although the discharge of the Colorado River varied substantially prior to the closure of Glen Canyon Dam in 1963, operation of the dam has caused changes in discharge that are more extreme than the pre-dam natural variability. Operation of the dam has eliminated flood flows and base flows, and thereby has effectively "flattened" the annual hydrograph. Prior to closure of the dam, the discharge of the Colorado River at Lees Ferry was lower than 7,980 ft 3 /s half of the time. Discharges lower than about 9,000 ft 3 /s were important for the seasonal accumulation and storage of sand in the pre-dam river downstream from Lees Ferry. The current operating plan for Glen Canyon Dam no longer allows sustained discharges lower than 8,000 ft 3 /s to be released. Thus, closure of the dam has not only cut off the upstream supply of sediment, but operation of the dam has also largely eliminated discharges during which sand could be demonstrated to accumulate in the river. In addition to radically changing the hydrology of the river, operation of the dam for hydroelectric-power generation has introduced large daily fluctuations in discharge. During the pre-dam era, the median daily range in discharge was only 542 ft 3 /s, although daily ranges in discharge exceeding 20,000 ft 3 /s were observed during the summer thunderstorm season. Relative to the pre-dam period of record, dam operations have increased the daily range in discharge during all but 0.1 percent of all days. The post-dam median daily range in discharge, 8,580 ft 3 /s, exceeds the pre-dam median discharge of 7,980 ft 3 /s. Operation of the dam has also radically changed the frequency of floods on the Colorado River at Lees Ferry. The frequency of floods with peak discharges larger than about 29,000 ft 3 /s has greatly decreased, while the frequency of smaller floods, with peak discharges between 18,500 and 29,000 ft 3 /s, has increased substantially. Operation of the dam has greatly extended the duration of smaller floods; for example, each of the four longest periods of sustained flows in excess of 18,500 ft 3 /s occurred after closure of the dam.

Arizona↗

Advancing cave detection using terrain analysis and thermal imagery

Since the initial experiments nearly 50 years ago, techniques for detecting caves using airborne and spacecraft acquired thermal imagery have improved markedly. These advances are largely due to a combination of higher instrument sensitivity, modern computing systems, and processor-intensive analytical techniques. Through applying these advancements, our goals were to: (1) Determine the efficacy of methods designed for terrain analysis and applied to thermal imagery; (2) evaluate the usefulness of predawn and midday imagery for detecting caves; and (3) ascertain which imagery type (predawn, midday, or the difference between those two times) was most informative. Using forward stepwise logistic (FSL) and Least Absolute Shrinkage and Selection Operator (LASSO) regression analyses for model selection, and a thermal imagery dataset acquired from the Mojave Desert, California, we examined the efficacy of three well-known terrain descriptors (i.e., slope, topographic position index (TPI), and curvature) on thermal imagery for cave detection. We also included the actual, untransformed thermal DN values (hereafter “unenhanced thermal”) as a fourth dataset. Thereafter, we compared the thermal signatures of known cave entrances to all non-cave surface locations. We determined these terrain-based analytical methods, which described the “shape” of the thermal landscape, hold significant promise for cave detection. All imagery types produced similar results. Down-selected covariates per imagery type, based upon the FSL models, were: Predawn— slope, TPI, curvature at 0 m from cave entrance, as well as slope at 1 m from cave entrance; midday— slope, TPI, and unenhanced thermal at 0 m from cave entrance; and difference— TPI and slope at 0 m from cave entrance, as well as unenhanced thermal and TPI at 3.5 m from cave entrance. We provide recommendations for future research directions in terrestrial and planetary cave detection using thermal imagery.

Remote Sensing↗

Slip distribution and rupture history of the August 11, 2012, double earthquakes in Ahar – Varzaghan, Iran, using joint inversion of teleseismic broadband and local strong motion data

We use combined teleseismic and strong motion data sets to investigate finite-fault slip models for a double of earthquakes that occurred on August 11, 2012, in northwestern Iran near the cities of Ahar and Varzaghan. The data include teleseismic P-waveforms retrieved from broadband seismic stations located between 30°–94° from the earthquakes and local strong motion data recorded by the Iran Strong Motion Network, installed and operated by the Building and Housing Research Centre. We first invert teleseismic P-waveforms and local strong motion data separately. For the first event (12:23 UTC), the teleseismic broadband inversion yields a somewhat deeper and simpler distribution of slip than the local strong motion inversion. The strong motion inversion results in a more complex distribution because of higher frequency content but can also be influenced by complexities in the propagation path. For the second event (12:34 UTC), the slip distribution from strong motion data is more similar to the teleseismic result and shows a simple slip area with a small relative movement to the west. To resolve the differences between the results of these two data sets and obtain a better constrained slip model, we perform a joint inversion of teleseismic broadband and local strong motion data. The joint inversion for the first event shows two asperities with a maximum slip of 3.9 m up- dip from the hypocenter and extending to the west between depths of 1 and 5 km. A second narrower high-slip area is seen just above the hypocenter from 6 to 10 km depth. The total moment for this earthquake is calculated to be M o = 3.8 × 10 25 dyn-cm (3.8 × 10 18 N.m) (M w 6.4). For the second event, the results of the joint inversion show a simple slip distribution that is mainly confined in a single patch around the hypocenter with a depth range from about 10 to 13 km and maximum slip of 1.9 m. We compute a total seismic moment of M o = 1.6 × 10 25 dyn-cm (1.6 × 10 18 N.m) (M w 6.1) for the second event. The largest stress drops for the first event occur above the hypocenter with an average stress drop over the rupture area of 120 bar (12 Mpa). For the second event, the maximum stress drop occurs at the reported focal depth with an average stress drop over the rupture area of 80 bar (8 Mpa).

East Anatolian Fault↗