Search USGS⌕ Search

SEARCH · Search USGS

Results for “Computational 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 91 records · Page 5Linked to original sources
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↗
USGS2022

gTOOLS, an open-source MATLAB program for processing high precision, relative gravity data for time-lapse gravity monitoring

gTOOLS is an open-source software for the processing of relative gravity data. gTOOLS is available in MATLAB and as a compiled executable to be run under the free MATLAB Runtime Compiler. The software has been designed for time-lapse (temporal) gravity monitoring. Although programmed to read the Scintrex CG-5 and CG-6 gravimeters output data files, it can be easily modified to read data files from other gravimeters. The software binds together single-task processing modules within a very simple user interface that is based on one text file. Gravity processing involves three modules: (a) gravimeter calibration; (b) automatic processing of gravity data to find adjusted gravity differences; and (c) post processing of results. Each module is optional and runs independently from the others. Data processing includes (a) averaging out the measurements noise, and correction for solid Earth tides, and ocean loading, and residual instrumental drift, and (b) calculate the residual instrumental drift and gravity differences between the base station and monitoring sites, and their uncertainties, by a weighted least square analysis of the gravity data. The software allows the automatic processing of a gravity campaign spanning multiple days in a single run. The software is tested on gravity data from 2015 eruption at Cotopaxi volcano, Ecuador.

Computers & Geosciences↗
USGS2022

MTH5: An archive and exchangeable data format for magnetotelluric time series data

Magnetotellurics (MT) is a passive electromagnetic geophysical method that measures variations in subsurface electrical resistivity. MT data are collected in the time domain and processed in the frequency domain to produce estimates of a transfer function representing the Earth’s electrical structure. Unfortunately, the MT community lacks metadata and data standards for time series data. As the community grows and findability, accessibility, interoperability, and reuse of digital assets (FAIR) data principles are enforced by government and funding agencies, a standard is needed for time series data. Presented here is a hierarchical data format (MTH5) that is logically formatted to how MT data are collected. Open-source Python packages are also described to read, write, and manipulate MTH5 files. These include a package to deal with metadata ( mt_metadata ) based on standards developed by the Working Group for Magnetotelluric Data Handling and Software assembled by the Incorporated Research Institutions for Seismology (IRIS), and mth5 : a package to interact with MTH5 files that uses mt_metadata . Example code and workflows are presented.

Computers & Geosciences↗
USGS2000

Lithofacies identification using multiple adaptive resonance theory neural networks and group decision expert system

Lithofacies identification supplies qualitative information about rocks. Lithofacies represent rock textures and are important components of hydrocarbon reservoir description. Traditional techniques of lithofacies identification from core data are costly and different geologists may provide different interpretations. In this paper, we present a low-cost intelligent system consisting of three adaptive resonance theory neural networks and a rule-based expert system to consistently and objectively identify lithofacies from well-log data. The input data are altered into different forms representing different perspectives of observation of lithofacies. Each form of input is processed by a different adaptive resonance theory neural network. Among these three adaptive resonance theory neural networks, one neural network processes the raw continuous data, another processes categorial data, and the third processes fuzzy-set data. Outputs from these three networks are then combined by the expert system using fuzzy inference to determine to which facies the input data should be assigned. Rules are prioritized to emphasize the importance of firing order. This new approach combines the learning ability of neural networks, the adaptability of fuzzy logic, and the expertise of geologists to infer facies of the rocks. This approach is applied to the Appleton Field, an oil field located in Escambia County, Alabama. The hybrid intelligence system predicts lithofacies identity from log data with 87.6% accuracy. This prediction is more accurate than those of single adaptive resonance theory networks, 79.3%, 68.0% and 66.0%, using raw, fuzzy-set, and categorical data, respectively, and by an error-backpropagation neural network, 57.3%. (C) 2000 Published by Elsevier Science Ltd. All rights reserved.

Computers & Geosciences↗
USGS2001

A two-dimensional, time-dependent model of suspended sediment transport and bed reworking for continental shelves

A two-dimensional, time-dependent solution to the transport equation is formulated to account for advection and diffusion of sediment suspended in the bottom boundary layer of continental shelves. This model utilizes a semi-implicit, upwind-differencing scheme to solve the advection-diffusion equation across a two-dimensional transect that is configured so that one dimension is the vertical, and the other is a horizontal dimension usually aligned perpendicular to shelf bathymetry. The model calculates suspended sediment concentration and flux; and requires as input wave properties, current velocities, sediment size distributions, and hydrodynamic sediment properties. From the calculated two-dimensional suspended sediment fluxes, we quantify the redistribution of shelf sediment, bed erosion, and deposition for several sediment sizes during resuspension events. The two-dimensional, time-dependent approach directly accounts for cross-shelf gradients in bed shear stress and sediment properties, as well as transport that occurs before steady-state suspended sediment concentrations have been attained. By including the vertical dimension in the calculations, we avoid depth-averaging suspended sediment concentrations and fluxes, and directly account for differences in transport rates and directions for fine and coarse sediment in the bottom boundary layer. A flux condition is used as the bottom boundary condition for the transport equation in order to capture time-dependence of the suspended sediment field. Model calculations demonstrate the significance of both time-dependent and spatial terms on transport and depositional patterns on continental shelves. ?? 2001 Elsevier Science Ltd. All rights reserved.

Computers & Geosciences↗
USGS2002

Description and comparison of geologic maps with FRAGSTATS - A spatial statistics program

FRAGSTATS is a public-domain GIS implementation of a set of spatial statistics that address a fundamental problem in GIS applications, description and comparison of maps. The spatial statistics from the 1:2,500,000-scale United States geologic map of Nevada, the central United States, and the northeastern United States quantify the differences in complexity and variability between these three geologic terranes. Nevada is defined by a large number of patches of small size and low size variability, whereas the Central area has a small number of patches with smaller relative size variability. All three areas have similar map-unit shape complexity with Nevada having the lowest. Based on the density of edges, the areas can be ranked from highest to lowest, as Nevada, Central, and Northeast. The Shannon diversity index ranks the areas from highest to lowest, as Northeast, Nevada, and Central, but the Shannon evenness index ranks them from highest to lowest, as Northeast, Central, and Nevada. These rankings may reflect the influence of folding in the Northeast and Central areas as opposed to basin and range extension in Nevada. The core areas statistic ranks the areas for spatial accuracy from highest to lowest, as Central, Northeast, and Nevada, with Northeast and Nevada being similar. For a scale comparison, the FRAGSTATS statistics quantify the increased complexity and spatial accuracy that is inherent in going from small- to larger-scale maps. For example for 1:2,500,000-1:500,000-scale maps of Nevada, respectively, the area weighted fractal dimension increase from 1.1 to 1.18, and the total core areas index almost doubles from 39.09 to 63.38. In addition, the fractal dimensions discriminate gross lithology and tectonic terranes. ?? 2002 Elsevier Science Ltd. All rights reserved.

Computers & Geosciences↗
USGS2001

1r2dinv: A finite-difference model for inverse analysis of two dimensional linear or radial groundwater flow

We have developed a program for inverse analysis of two-dimensional linear or radial groundwater flow problems. The program, 1r2dinv, uses standard finite difference techniques to solve the groundwater flow equation for a horizontal or vertical plane with heterogeneous properties. In radial mode, the program simulates flow to a well in a vertical plane, transforming the radial flow equation into an equivalent problem in Cartesian coordinates. The physical parameters in the model are horizontal or x-direction hydraulic conductivity, anisotropy ratio (vertical to horizontal conductivity in a vertical model, y-direction to x-direction in a horizontal model), and specific storage. The program allows the user to specify arbitrary and independent zonations of these three parameters and also to specify which zonal parameter values are known and which are unknown. The Levenberg-Marquardt algorithm is used to estimate parameters from observed head values. Particularly powerful features of the program are the ability to perform simultaneous analysis of heads from different tests and the inclusion of the wellbore in the radial mode. These capabilities allow the program to be used for analysis of suites of well tests, such as multilevel slug tests or pumping tests in a tomographic format. The combination of information from tests stressing different vertical levels in an aquifer provides the means for accurately estimating vertical variations in conductivity, a factor profoundly influencing contaminant transport in the subsurface. ?? 2001 Elsevier Science Ltd. All rights reserved.

Computers & Geosciences↗