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

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↗
USGS2026

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↗
USGS2014

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↗
USGS2012

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↗
USGS2011

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↗
USGS1978

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↗
USGS1990

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↗
USGS2012

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↗
USGS2015

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↗
USGS2003

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↗
USGS2021

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↗
USGS2021

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↗
USGS2021

Heat flux from a vapor-dominated hydrothermal field beneath Yellowstone Lake

We report results from 149 heat flux measurements made over n ∼2-year interval at sites in and around a vapor-dominated geothermal field located at water depths of ∼100–120 m in Yellowstone Lake, Wyoming. Measurements of both in situ temperature and thermal conductivity as a function of depth were made with a 1 m probe via a remotely operated vehicle, and are combined to compute the vertical conductive heat flux. Inside the ∼55.5 × 10 3 m 2 bathymetric depression demarcating the vapor-dominated field, the median conductive flux is 13 W m −2 , with a conductive output of 0.72 MW. Outside the thermal field, the median conductive flux is 3.5 W m −2 . We observed 49 active vents inside the thermal field, with an estimated mass discharge rate of 56 kg s −1 , a median exit-fluid temperature of 132°C, and a total heat output of 29 MW. We find evidence for relatively weak secondary convection with a total output of 0.09 MW in thermal area lake floor sediments. Our data indicate that vapor beneath the thermal field is trapped by a low-permeability cap at a temperature of ∼189°C and a depth of ∼15 m below the lake floor. The thermal output of the Deep Hole is among the highest of any vapor-dominated field in Yellowstone, due in part to the high boiling temperatures associated with the elevated lake floor pressures.

Wyoming↗
USGS1971

Field and office instructions in stream gauging for the Hydrological Survey of Zambia

The importance of water to the basic needs of man is self- evident and needs no particular emphasis. The importance of water to a developing economy cannot be overemphasized. A few decades ago, hydrology was a division of hydraulic engineering and was a tool for project survey, plan, and design. Today hydrology still remains an important part of planning and management of water use projects, but it is imperative that surface and ground-water basic data networks be designed and operated from the standpoint of both present and future water needs. Water problems are ever increasing and ever changing and preparation for the future water demands of Zambia requires knowledge of the hydrology of the country instead of the examination of piecemeal samples for each water use project. The hydrologic survey of Zambia needs to be under the guidance of competent and imaginative hydrologists solidly trained in all elements of basic data collection and analysis and not in the hands of water project planners. Hydrology is a science which requires the highest order of teamwork and the hydrologist will need the help and advice of many employees within the organization to operate the network, provide adequate research, and examine the water needs of the country. It must be thoroughly understood that communication is essential between the hydrological survey and water project planners from both the government and private sectors. It is very important to define the aims and duties of the Hydrological Branch of the Water Affairs Department in a clear cut "Statement of Policy". Personal copies of the statement should be made available to all professional employees and technicians. The reasons for the existence of the Branch may be self-apparent to heads of the organization, but to all other employees the reasons may be vague and unknown. Every member of the technical and administrative staff would benefit by an understanding of the purpose of his work. Nebulous ideas of the function of a hydrological survey of Zambia serve little purpose to the employee making the plans or those who execute them. A dissemination and free flow of ideas in all directions is necessary in any scientific and creative organization. The primary functions of the Hydrological Branch of the Water Affairs Department of Zambia are twofold: Water resource appraisal which includes the inventory of all surface and ground waters; Research, both basic and applied, in the science of water. Possibly some function regarding water law is also conceivable, although attachment of a law enforcement function to water-resource investigations might weaken the position of the Branch as an unbiased scientific organization. As the collection of basic streamflow data is a primary function of the Branch the following instructions for the collection, analysis and computation of streamflow records have been written.

Open-File Report↗
USGS2019

Mapping cropland extent of Southeast and Northeast Asia using multi-year time-series Landsat 30-m data using Random Forest classifier on Google Earth Engine

Cropland extent maps are useful components for assessing food security. Ideally, such products are a useful addition to countrywide agricultural statistics since they are not politically biased and can be used to calculate cropland area for any spatial unit from an individual farm to various administrative unites (e.g., state, county, district) within and across nations, which in turn can be used to estimate agricultural productivity as well as degree of disturbance on food security from natural disasters and political conflict. However, existing cropland extent maps over large areas (e.g., Country, region, continent, world) are derived from coarse resolution imagery (250 m to 1 km pixels) and have many limitations such as missing fragmented and\or small farms with mixed signatures from different crop types and\or farming practices that can be, confused with other land cover. As a result, the coarse resolution maps have limited useflness in areas where fields are small (<1 ha), such as in Southeast Asia. Furthermore, coarse resolution cropland maps have known uncertainties in both geo-precision of cropland location as well as accuracies of the product. To overcome these limitations, this research was conducted using multi-date, multi-year 30-m Landsat time-series data for 3 years chosen from 2013 to 2016 for all Southeast and Northeast Asian Countries (SNACs), which included 7 refined agro-ecological zones (RAEZ) and 12 countries (Indonesia, Thailand, Myanmar, Vietnam, Malaysia, Philippines, Cambodia, Japan, North Korea, Laos, South Korea, and Brunei). The 30-m (1 pixel = 0.09 ha) data from Landsat 8 Operational Land Imager (OLI) and Landsat 7 Enhanced Thematic Mapper (ETM+) were used in the study. Ten Landsat bands were used in the analysis (blue, green, red, NIR, SWIR1, SWIR2, Thermal, NDVI, NDWI, LSWI) along with additional layers of standard deviation of these 10 bands across 1 year, and global digital elevation model (GDEM)-derived slope and elevation bands. To reduce the impact of clouds, the Landsat imagery was time-composited over four time-periods (Period 1: January- April, Period 2: May-August, and Period 3: September-December) over 3-years. Period 4 was the standard deviation of all 10 bands taken over all images acquired during the 2015 calendar year. These four period composites, totaling 42 band data-cube, were generated for each of the 7 RAEZs. The reference training data (N = 7849) generated for the 7 RAEZ using sub-meter to 5-m very high spatial resolution imagery (VHRI) helped generate the knowledge-base to separate croplands from non-croplands. This knowledge-base was used to code and run a pixel-based random forest (RF) supervised machine learning algorithm on the Google Earth Engine (GEE) cloud computing environment to separate croplands from non-croplands. The resulting cropland extent products were evaluated using an independent reference validation dataset (N = 1750) in each of the 7 RAEZs as well as for the entire SNAC area. For the entire SNAC area, the overall accuracy was 88.1% with a producer’s accuracy of 81.6% (errors of omissions = 18.4%) and user’s accuracy of 76.7% (errors of commissions = 23.3%). For each of the 7 RAEZs overall accuracies varied from 83.2 to 96.4%. Cropland areas calculated for the 12 countries were compared with country areas reported by the United Nations Food and Agriculture Organization and other national cropland statistics resulting in an R 2 value of 0.93. The cropland areas of provinces were compared with the province statistics that showed an R 2 = 0.95 for South Korea and R 2 = 0.94 for Thailand. The cropland products are made available on an interactive viewer at www.croplands.org and for download at National Aeronautics and Space Administration’s (NASA) Land Processes Distributed Active Archive Center (LP DAAC): https://lpdaac.usgs.gov/node/1281 .

International Journal of Applied Earth Observation↗
USGS2000

The 1999 Southern California Seismic Network bulletin

The Pasadena office of the U.S. Geological Survey (USGS), together with the Caltech Seismological Laboratory, operates a network of more than 350 remote seismometers in southern California called the S outhern C alifornia S eismic N etwork (SCSN). SCSN is part of TriNet, a cooperative project between the USGS, Caltech, and the California Division of Mines and Geology (CDMG). The TriNet project is halfway completed and is upgrading the existing network to digital, adding new stations, and developing real-time and earthquake-alert capabilities. Signals from the SCSN sites are telemetered to a central processing location at the Caltech Seismological Lab in Pasadena. Computers that detect and record thousands of earthquakes each year continuously monitor these signals. Phase arrival times for these events are picked by analysts and are archived along with digital seismograms. Data acquisition, processing, and archiving are achieved using the Caltech/USGS Seismic Processing (CUSP) system ( Dollar, 1989 ). These data have been compiled into the SCSN Catalog of Earthquakes, a list beginning in 1932 that currently contains more than 356,300 events. Waveform, phase, and catalog data are archived by the Southern California Earthquake Center Data Center (SCEC_DC). This data set is critical to the evaluation of earthquake hazards in California and to the advancement of geoscience as a whole.

California↗
USGS1998

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2024

Developing a photography-based harvest survey to estimate age and subspecies composition of midcontinent sandhill cranes

Midcontinent sandhill cranes ( Antigone canadensis ) are managed as a single population, but hunting regulations are structured so harvest is targeted towards the more numerous lesser sandhill cranes ( A. c. canadensis ). However, research indicates that greater sandhill cranes ( A. c. tabida ) have been disproportionally exposed to harvest at a rate exceeding their proportion within the midcontinent population. In addition, harvest has increased 22% per year in the U.S. Central Flyway states. The midcontinent population appears to be growing in recent years, but variability in annual abundance estimates has increased substantially. With limited resources and harvest management uncertainty increasing, we developed methods for a citizen science, photography-based harvest survey to estimate age and subspecies composition of harvested midcontinent sandhill cranes. To develop survey methods, we collected physical parts from 284 sandhill cranes in North Dakota in 2019 and 2020. We manually measured the culmen and tarsus using calipers, and digitally measured these parts using photographs and computer software. All digitally derived measurements were 2.5% to 5.9% larger than manual measurements; therefore, we developed linear models that adjusted digital measurements, facilitating subspecies prediction using an existing morphometric-based technique. In 2021, we requested an equal number of hunters to participate using 2 data collection methods to test if hunters could reliably take photographs suitable for digital measurement. Collection method 1 involved photographing the head and leg simultaneously, and Collection method 2 involved photographing the head only. Hunters submitted a total of 239 photographs. Only 80 of these photographs were submitted using Collection method 1, and 72% were suitable for digital measurement. Conversely, hunters submitted twice as many photographs using Collection method 2, and 88% of these photographs were deemed suitable. Although obtaining the tarsus measurement slightly improved subspecies predictability, Collection method 2 increased participation and usable data. We believe our results could be used to develop an operational survey of age and subspecies composition of midcontinent sandhill cranes, wherein a sample of crane hunters throughout the midcontinent population range would be asked to electronically submit a photograph of the head of each bird they harvest. A time series of age and subspecies composition of this population would provide managers with valuable information and improve harvest management at minimal additional cost and burden, compared to a traditional parts collection survey administered by mail.

Wildlife Society Bulletin↗