USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computers & Geosciences”

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 109 records · Page 6Linked to original sources
USGS2012

Methods for simulating solute breakthrough curves in pumping groundwater wells

In modeling there is always a trade-off between execution time and accuracy. For gradient-based parameter estimation methods, where a simulation model is run repeatedly to populate a Jacobian (sensitivity) matrix, there exists a need for rapid simulation methods of known accuracy that can decrease execution time, and thus make the model more useful without sacrificing accuracy. Convolution-based methods can be executed rapidly for any desired input function once the residence-time distribution is known. The residence-time distribution can be calculated efficiently using particle tracking, but particle tracking can be ambiguous near a pumping well if the grid is too coarse. We present several embedded analytical expressions for improving particle tracking near a pumping well and compare them with a finely gridded finite-difference solution in terms of accuracy and CPU usage. Even though the embedded analytical approach can improve particle tracking near a well, particle methods reduce, but do not eliminate, reliance on a grid because velocity fields typically are calculated on a grid, and additional error is incurred using linear interpolation of velocity. A dilution rate can be calculated for a given grid and pumping well to determine if the grid is sufficiently refined. Embedded analytical expressions increase accuracy but add significantly to CPU usage. Structural error introduced by the numerical solution method may affect parameter estimates.

Computers & Geosciences↗
USGS2012

Relative azimuth inversion by way of damped maximum correlation estimates

Horizontal seismic data are utilized in a large number of Earth studies. Such work depends on the published orientations of the sensitive axes of seismic sensors relative to true North. These orientations can be estimated using a number of different techniques: SensOrLoc (Sensitivity, Orientation and Location), comparison to synthetics (Ekstrom and Busby, 2008), or by way of magnetic compass. Current methods for finding relative station azimuths are unable to do so with arbitrary precision quickly because of limitations in the algorithms (e.g. grid search methods). Furthermore, in order to determine instrument orientations during station visits, it is critical that any analysis software be easily run on a large number of different computer platforms and the results be obtained quickly while on site. We developed a new technique for estimating relative sensor azimuths by inverting for the orientation with the maximum correlation to a reference instrument, using a non-linear parameter estimation routine. By making use of overlapping windows, we are able to make multiple azimuth estimates, which helps to identify the confidence of our azimuth estimate, even when the signal-to-noise ratio (SNR) is low. Finally, our algorithm has been written as a stand-alone, platform independent, Java software package with a graphical user interface for reading and selecting data segments to be analyzed.

Computers & Geosciences↗
USGS2012

Social.Water - A crowdsourcing tool for environmental data acquisition

Remote telemetry has a long history of use for collection of environmental measurements. With the rise of mobile phones and SMS text-messaging capacity, many members of the general pubic carry communications equipment in their pockets at all times. Enabling the general public to provide environmental data through text messages has the potential both to provide additional data to scientific projects and also to raise awareness of the projects through participation. Hydrologic measurements – some of which can be made without training, involve a single measurement, and are often made in rural areas – are well-suited to text-message conveyance. Many other environmental measurements are similarly well-suited for this technology. Social.Water is a software package, written in Python, that collects, parses, and categorizes text messages sent to a dedicated phone number, updates a simple database, and posts both graphical results and the database on the Web. Social.Water was designed as the backend to the Crowdhydrology project and is written in an object-oriented design that makes customization and modification straightforward.

Computers & Geosciences↗
USGS2013

extrap: Software to assist the selection of extrapolation methods for moving-boat ADCP streamflow measurements

Selection of the appropriate extrapolation methods for computing the discharge in the unmeasured top and bottom parts of a moving-boat acoustic Doppler current profiler (ADCP) streamflow measurement is critical to the total discharge computation. The software tool, extrap, combines normalized velocity profiles from the entire cross section and multiple transects to determine a mean profile for the measurement. The use of an exponent derived from normalized data from the entire cross section is shown to be valid for application of the power velocity distribution law in the computation of the unmeasured discharge in a cross section. Selected statistics are combined with empirically derived criteria to automatically select the appropriate extrapolation methods. A graphical user interface (GUI) provides the user tools to visually evaluate the automatically selected extrapolation methods and manually change them, as necessary. The sensitivity of the total discharge to available extrapolation methods is presented in the GUI. Use of extrap by field hydrographers has demonstrated that extrap is a more accurate and efficient method of determining the appropriate extrapolation methods compared with tools currently (2012) provided in the ADCP manufacturers’ software.

Computers & Geosciences↗
USGS2013

A wetting and drying scheme for ROMS

The processes of wetting and drying have many important physical and biological impacts on shallow water systems. Inundation and dewatering effects on coastal mud flats and beaches occur on various time scales ranging from storm surge, periodic rise and fall of the tide, to infragravity wave motions. To correctly simulate these physical processes with a numerical model requires the capability of the computational cells to become inundated and dewatered. In this paper, we describe a method for wetting and drying based on an approach consistent with a cell-face blocking algorithm. The method allows water to always flow into any cell, but prevents outflow from a cell when the total depth in that cell is less than a user defined critical value. We describe the method, the implementation into the three-dimensional Regional Oceanographic Modeling System (ROMS), and exhibit the new capability under three scenarios: an analytical expression for shallow water flows, a dam break test case, and a realistic application to part of a wetland area along the Georgia Coast, USA.

Georgia↗
USGS2014

MTpy: A Python toolbox for magnetotellurics

We present the software package MTpy that allows handling, processing, and imaging of magnetotelluric (MT) data sets. Written in Python, the code is open source, containing sub-packages and modules for various tasks within the standard MT data processing and handling scheme. Besides the independent definition of classes and functions, MTpy provides wrappers and convenience scripts to call standard external data processing and modelling software. In its current state, modules and functions of MTpy work on raw and pre-processed MT data. However, opposite to providing a static compilation of software, we prefer to introduce MTpy as a flexible software toolbox, whose contents can be combined and utilised according to the respective needs of the user. Just as the overall functionality of a mechanical toolbox can be extended by adding new tools, MTpy is a flexible framework, which will be dynamically extended in the future. Furthermore, it can help to unify and extend existing codes and algorithms within the (academic) MT community. In this paper, we introduce the structure and concept of MTpy . Additionally, we show some examples from an everyday work-flow of MT data processing: the generation of standard EDI data files from raw electric ( E -) and magnetic flux density ( B -) field time series as input, the conversion into MiniSEED data format, as well as the generation of a graphical data representation in the form of a Phase Tensor pseudosection.

Computers & Geosciences↗
USGS2015

EverVIEW: a visualization platform for hydrologic and Earth science gridded data

The EverVIEW Data Viewer is a cross-platform desktop application that combines and builds upon multiple open source libraries to help users to explore spatially-explicit gridded data stored in Network Common Data Form (NetCDF). Datasets are displayed across multiple side-by-side geographic or tabular displays, showing colorized overlays on an Earth globe or grid cell values, respectively. Time-series datasets can be animated to see how water surface elevation changes through time or how habitat suitability for a particular species might change over time under a given scenario. Initially targeted toward Florida's Everglades restoration planning, EverVIEW has been flexible enough to address the varied needs of large-scale planning beyond Florida, and is currently being used in biological planning efforts nationally and internationally.

Computers & Geosciences↗
USGS2015

The data quality analyzer: a quality control program for seismic data

The U.S. Geological Survey's Albuquerque Seismological Laboratory (ASL) has several initiatives underway to enhance and track the quality of data produced from ASL seismic stations and to improve communication about data problems to the user community. The Data Quality Analyzer (DQA) is one such development and is designed to characterize seismic station data quality in a quantitative and automated manner. The DQA consists of a metric calculator, a PostgreSQL database, and a Web interface: The metric calculator, SEEDscan, is a Java application that reads and processes miniSEED data and generates metrics based on a configuration file. SEEDscan compares hashes of metadata and data to detect changes in either and performs subsequent recalculations as needed. This ensures that the metric values are up to date and accurate. SEEDscan can be run as a scheduled task or on demand. The PostgreSQL database acts as a central hub where metric values and limited station descriptions are stored at the channel level with one-day granularity. The Web interface dynamically loads station data from the database and allows the user to make requests for time periods of interest, review specific networks and stations, plot metrics as a function of time, and adjust the contribution of various metrics to the overall quality grade of the station. The quantification of data quality is based on the evaluation of various metrics (e.g., timing quality, daily noise levels relative to long-term noise models, and comparisons between broadband data and event synthetics). Users may select which metrics contribute to the assessment and those metrics are aggregated into a “grade” for each station. The DQA is being actively used for station diagnostics and evaluation based on the completed metrics (availability, gap count, timing quality, deviation from a global noise model, deviation from a station noise model, coherence between co-located sensors, and comparison between broadband data and synthetics for earthquakes) on stations in the Global Seismographic Network and Advanced National Seismic System.

Computers & Geosciences↗
USGS2015

Spatially explicit spectral analysis of point clouds and geospatial data

The increasing use of spatially explicit analyses of high-resolution spatially distributed data (imagery and point clouds) for the purposes of characterising spatial heterogeneity in geophysical phenomena necessitates the development of custom analytical and computational tools. In recent years, such analyses have become the basis of, for example, automated texture characterisation and segmentation, roughness and grain size calculation, and feature detection and classification, from a variety of data types. In this work, much use has been made of statistical descriptors of localised spatial variations in amplitude variance (roughness), however the horizontal scale (wavelength) and spacing of roughness elements is rarely considered. This is despite the fact that the ratio of characteristic vertical to horizontal scales is not constant and can yield important information about physical scaling relationships. Spectral analysis is a hitherto under-utilised but powerful means to acquire statistical information about relevant amplitude and wavelength scales, simultaneously and with computational efficiency. Further, quantifying spatially distributed data in the frequency domain lends itself to the development of stochastic models for probing the underlying mechanisms which govern the spatial distribution of geological and geophysical phenomena. The software package PySESA (Python program for Spatially Explicit Spectral Analysis) has been developed for generic analyses of spatially distributed data in both the spatial and frequency domains. Developed predominantly in Python , it accesses libraries written in Cython and C++ for efficiency. It is open source and modular, therefore readily incorporated into, and combined with, other data analysis tools and frameworks with particular utility for supporting research in the fields of geomorphology, geophysics, hydrography, photogrammetry and remote sensing. The analytical and computational structure of the toolbox is described, and its functionality illustrated with an example of a high-resolution bathymetric point cloud data collected with multibeam echosounder.

Computers & Geosciences↗
USGS2017

Computing spatial correlation of ground motion intensities for ShakeMap

Modeling the spatial correlation of ground motion residuals, caused by coherent contributions from source, path, and site, can provide valuable loss and hazard information, as well as a more realistic depiction of ground motion intensities. The U.S. Geological Survey (USGS) software package, ShakeMap, utilizes a deterministic empirical approach to estimate median ground shaking in conjunction with observed seismic data. ShakeMap-based shaking estimates are used in concert with loss estimation algorithms to estimate fatalities and economic losses after significant seismic events around the globe. Incorporating the spatial correlation of ground motion residuals has been shown to improve seismic loss estimates. In particular, Park, Bazzuro, and Baker (Applications of Statistics and Probability in Civil Engineering, 2007) investigated computing spatially correlated random fields of residuals. However, for large scale ShakeMap grids, computational requirements of the method are prohibitive. In this work, a memory efficient algorithm is developed to compute the random fields and implemented using the ShakeMap framework. This new, iterative parallel algorithm is based on decay properties of an associated ground motion correlation function and is shown to significantly reduce computational requirements associated with adding spatial variability to the ShakeMap g round motion estimates. Further, we demonstrate and quantify the impact of adding peak ground motion spatial variability on resulting earthquake loss estimates.

Computers & Geosciences↗
USGS2017

Development of a coupled wave-flow-vegetation interaction model

Emergent and submerged vegetation can significantly affect coastal hydrodynamics. However, most deterministic numerical models do not take into account their influence on currents, waves, and turbulence. In this paper, we describe the implementation of a wave-flow-vegetation module into a Coupled-Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system that includes a flow model (ROMS) and a wave model (SWAN), and illustrate various interacting processes using an idealized shallow basin application. The flow model has been modified to include plant posture-dependent three-dimensional drag, in-canopy wave-induced streaming, and production of turbulent kinetic energy and enstrophy to parameterize vertical mixing. The coupling framework has been updated to exchange vegetation-related variables between the flow model and the wave model to account for wave energy dissipation due to vegetation. This study i) demonstrates the validity of the plant posture-dependent drag parameterization against field measurements, ii) shows that the model is capable of reproducing the mean and turbulent flow field in the presence of vegetation as compared to various laboratory experiments, iii) provides insight into the flow-vegetation interaction through an analysis of the terms in the momentum balance, iv) describes the influence of a submerged vegetation patch on tidal currents and waves separately and combined, and v) proposes future directions for research and development.

Computers & Geosciences↗
USGS2017

Large crater clustering tool

In this paper we present the Large Crater Clustering (LCC) tool set, an ArcGIS plugin that supports the quantitative approximation of a primary impact location from user-identified locations of possible secondary impact craters or the long-axes of clustered secondary craters. The identification of primary impact craters directly supports planetary geologic mapping and topical science studies where the chronostratigraphic age of some geologic units may be known, but more distant features have questionable geologic ages. Previous works (e.g., McEwen et al., 2005; Dundas and McEwen, 2007) have shown that the source of secondary impact craters can be estimated from secondary impact craters. This work adapts those methods into a statistically robust tool set. We describe the four individual tools within the LCC tool set to support: (1) processing individually digitized point observations (craters), (2) estimating the directional distribution of a clustered set of craters, back projecting the potential flight paths (crater clusters or linearly approximated catenae or lineaments), (3) intersecting projected paths, and (4) intersecting back-projected trajectories to approximate the local of potential source primary craters. We present two case studies using secondary impact features mapped in two regions of Mars. We demonstrate that the tool is able to quantitatively identify primary impacts and supports the improved qualitative interpretation of potential secondary crater flight trajectories.

Computers & Geosciences↗
USGS2017

HERA: A dynamic web application for visualizing community exposure to flood hazards based on storm and sea level rise scenarios

The Hazard Exposure Reporting and Analytics (HERA) dynamic web application was created to provide a platform that makes research on community exposure to coastal-flooding hazards influenced by sea level rise accessible to planners, decision makers, and the public in a manner that is both easy to use and easily accessible. HERA allows users to (a) choose flood-hazard scenarios based on sea level rise and storm assumptions, (b) appreciate the modeling uncertainty behind a chosen hazard zone, (c) select one or several communities to examine exposure, (d) select the category of population or societal asset, and (e) choose how to look at results. The application is designed to highlight comparisons between (a) varying levels of sea level rise and coastal storms, (b) communities, (c) societal asset categories, and (d) spatial scales. Through a combination of spatial and graphical visualizations, HERA aims to help individuals and organizations to craft more informed mitigation and adaptation strategies for climate-driven coastal hazards. This paper summarizes the technologies used to maximize the user experience, in terms of interface design, visualization approaches, and data processing.

Computers & Geosciences↗
USGS2018

Determining on-fault earthquake magnitude distributions from integer programming

Earthquake magnitude distributions among faults within a fault system are determined from regional seismicity and fault slip rates using binary integer programming. A synthetic earthquake catalog (i.e., list of randomly sampled magnitudes) that spans millennia is first formed, assuming that regional seismicity follows a Gutenberg-Richter relation. Each earthquake in the synthetic catalog can occur on any fault and at any location. The objective is to minimize misfits in the target slip rate for each fault, where slip for each earthquake is scaled from its magnitude. The decision vector consists of binary variables indicating which locations are optimal among all possibilities. Uncertainty estimates in fault slip rates provide explicit upper and lower bounding constraints to the problem. An implicit constraint is that an earthquake can only be located on a fault if it is long enough to contain that earthquake. A general mixed-integer programming solver, consisting of a number of different algorithms, is used to determine the optimal decision vector. A case study is presented for the State of California, where a 4 kyr synthetic earthquake catalog is created and faults with slip ≥3 mm/yr are considered, resulting in >10 6 variables. The optimal magnitude distributions for each of the faults in the system span a rich diversity of shapes, ranging from characteristic to power-law distributions.

Computers & Geosciences↗
USGS2019

AVO-G2S: A modified, open-source Ground-to-Space atmospheric specification for infrasound modeling

To facilitate infrasound propagation studies, we present AVO-G2S, an open-source, Ground-to-Space model which provides temperature and wind specifications from the surface to an altitude of 225 km. This model provides a means of smoothly characterizing atmospheric conditions using multiple numerical weather prediction forecast and reanalysis products, along with upper-atmospheric empirical models. Regional atmospheric reconstructions only require a limited domain and can utilize high-resolution numerical weather prediction forecasts typically provided on a projected grid. The use of a projected grid allows for faster spectral transform libraries to be employed. The AVO-G2S software can also provide global reconstructions that rely on global numerical weather prediction products and spherical harmonic decompositions. AVO-G2S is inspired by a global Ground-to-Space model developed by the Naval Research Laboratory, and relies on their empirical descriptions of upper-atmospheric conditions. Alaska Volcano Observatory has implemented this model for near-real-time infrasound monitoring of volcanic eruptions and historical research projects.

Computers & Geosciences↗
USGS2019

HyCReWW: A hybrid coral reef wave and water level metamodel

Wave-induced flooding is a major coastal hazard on tropical islands fronted by coral reefs. The variability of shape, size, and physical characteristics of the reefs across the globe make it difficult to obtain a parameterization of wave run-up, which is needed for risk assessments. Therefore, we developed the HyCReWW metamodel to predict wave run-up under a wide range of reef morphometric and offshore forcing characteristics. Due to the complexity and high dimensionality of the problem, we assumed an idealized one-dimensional reef profile, characterized by seven primary parameters. XBeach Non-Hydrostatic was chosen to create the synthetic dataset, and Radial Basis Functions implemented in MATLAB® were chosen for interpolation. Results demonstrate the applicability of the metamodel to obtain fast and accurate results of wave run-up for a large range of intrinsic reef morphologic and extrinsic hydrodynamic forcing parameters, offering a useful tool for risk management and early warning systems.

Computers & Geosciences↗
USGS2020

Acoustic Sediment Estimation Toolbox (ASET): A software package for calibrating and processing TRDI ADCP data to compute suspended-sediment transport in sandy rivers

Quantifying suspended-sediment transport is critical for a variety of disciplines related to the management of water resources. However, the number of gauging stations and monitoring networks in most rivers around the world is insufficient to improve understanding of river dynamics and support water resource management decisions. This is mainly due to the high operational costs and intensive labor involved in traditional sediment measurement techniques, especially in sand bed rivers where coarse material varies spatially in the river cross section. Recently, the acoustic surrogate method has received attention as a potentially accurate surrogate technology for estimating suspended-sediment concentrations. In addition, the acoustic surrogate method, through use of acoustic Doppler current profilers (ADCPs), has the advantage of being able to simultaneously measure the flow velocity field and cross-sectional area when moving-boat measurements are performed. In spite of the important advances made in the implementation of this technique, there are no widely-available, free tools for processing the ADCP acoustic signal cross section measurements which include options to extrapolate velocity and sediment in unmeasured ADCP zones and develop calibrations with physical samples. This paper presents a new software called Acoustic Sediment Estimation Toolbox (ASET), which enables the user to develop a calibration between the acoustic signal collected with a down-looking Teledyne RD Instruments ADCP and sediment concentrations determined using traditional sediment sampling techniques. Moreover, ASET software uses dynamic ADCP measurements to estimate the total suspended-sediment transport through a river cross section. The theoretical framework and data processing routines applied by each module in ASET are presented. Finally, a comparison is made between the results obtained by ASET and by traditional methodologies for computing suspended-sediment transport in a large river system (Paraná River, Argentina).

Computers & Geosciences↗
USGS2021

Revisiting the declustering of spatial data with preferential sampling

Preferential sampling is a form of data collection that may significantly distort the histogram and the semivariogram of spatially correlated data . Typical situations are a higher sampling density at high-valued areas favorable for mining, and highly contaminated areas in need of environmental remediation. Multiple statistical procedures are devoted to obtaining representative statistics, whose magnitudes should be close to the respective population values. This paper proposes a resampling method that can compensate for preferential sampling of spatially correlated data without using declustering weights. The application of the method herein generates a dataset of median estimates of quantiles of multiple stratified resamples that is free of preferential sampling. The methodology is illustrated with two examples. The first one involves values actually measured in the field and has the advantage of representing a real scenario of spatial fluctuations and preferential sampling. A second dataset is synthetic and has the main benefit of a priori knowledge of the underlying spatial distribution, thus allowing a satisfactory evaluation of the results against the known baseline. Access to computer code is offered for practical application of the method.

Computers & Geosciences↗