Program disk to accompany MAGMIX (Computers & Geosciences, 1988, v. 14, p. 213-228)for IBM PC/AT (or compatible) computers
Explore the source record for details and available documents.
Open-File Report↗
SEARCH · USGS Science
Results for “Computers and 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.
Earthcasting: Geomorphic forecasts for society
Over the last several decades, the study of Earth surface processes has progressed from a descriptive science to an increasingly quantitative one due to advances in theoretical, experimental, and computational geosciences. The importance of geomorphic forecasts has never been greater, as technological development and global climate change threaten to reshape the landscapes that support human societies and natural ecosystems. Here we explore best practices for developing socially relevant forecasts of Earth surface change, a goal we are calling “earthcasting”. We suggest that earthcasts have the following features: they focus on temporal (∼1–∼100 years) and spatial (∼1 m–∼10 km) scales relevant to planning; they are designed with direct involvement of stakeholders and public beneficiaries through the evaluation of the socioeconomic impacts of geomorphic processes; and they generate forecasts that are clearly stated, testable, and include quantitative uncertainties. Earthcasts bridge the gap between Earth surface researchers and decision-makers, stakeholders, researchers from other disciplines, and the general public. We investigate the defining features of earthcasts and evaluate some specific examples. This paper builds on previous studies of prediction in geomorphology by recommending a roadmap for (a) generating earthcasts, especially those based on modeling; (b) transforming a subset of geomorphic research into earthcasts; and (c) communicating earthcasts beyond the geomorphology research community. Earthcasting exemplifies the social benefit of geomorphology research, and it calls for renewed research efforts toward further understanding the limits of predictability of Earth surface systems and processes, and the uncertainties associated with modeling geomorphic processes and their impacts.
Earth's Future↗
A BASIC program for locating references cited in geoscience manuscripts
[No abstract available]
Computers & Geosciences↗
Statistical analysis of water-quality data containing multiple detection limits II: S-language software for nonparametric distribution modeling and hypothesis testing
Analysis of low concentrations of trace contaminants in environmental media often results in left-censored data that are below some limit of analytical precision. Interpretation of values becomes complicated when there are multiple detection limits in the data-perhaps as a result of changing analytical precision over time. Parametric and semi-parametric methods, such as maximum likelihood estimation and robust regression on order statistics, can be employed to model distributions of multiply censored data and provide estimates of summary statistics. However, these methods are based on assumptions about the underlying distribution of data. Nonparametric methods provide an alternative that does not require such assumptions. A standard nonparametric method for estimating summary statistics of multiply-censored data is the Kaplan-Meier (K-M) method. This method has seen widespread usage in the medical sciences within a general framework termed "survival analysis" where it is employed with right-censored time-to-failure data. However, K-M methods are equally valid for the left-censored data common in the geosciences. Our S-language software provides an analytical framework based on K-M methods that is tailored to the needs of the earth and environmental sciences community. This includes routines for the generation of empirical cumulative distribution functions, prediction or exceedance probabilities, and related confidence limits computation. Additionally, our software contains K-M-based routines for nonparametric hypothesis testing among an unlimited number of grouping variables. A primary characteristic of K-M methods is that they do not perform extrapolation and interpolation. Thus, these routines cannot be used to model statistics beyond the observed data range or when linear interpolation is desired. For such applications, the aforementioned parametric and semi-parametric methods must be used.
Computers & Geosciences↗
3D semantic mapping of surface geological features
Semantic mapping in 3D is fundamental to a wide range of geoscientific studies and applications, including geomorphology, hazard assessment, and environmental monitoring. However, automatically segmenting geological features from large-scale photogrammetric datasets remains a significant challenge. We present a methodology to address this gap. Using overlapping images collected over environments of interest, Structure-from-Motion (SfM) produces georeferenced point clouds and estimates camera poses. Existing large vision models, such as Segment Anything Model, segment objects in the images, generating pixel-segmentation associations. To produce pixel-point associations, we project the points back onto the camera image planes. As objects are independently segmented across multiple images with different perspectives, we develop a segmentation mosaicking algorithm to build probabilistic point-segmentation associations that combines the pixel-segmentation associations and pixel-point associations. Our methodology is validated using both synthetic data generated by Kubric and real-world UAV-SfM data. The implementation is designed to be compatible with existing SfM software, including Agisoft and OpenDroneMap, for photogrammetry mapping in geoscience studies. As a case study, we apply our method to the semantic mapping of precariously balanced rocks (PBRs), which provide upper-bound constraints on historical ground motion shaking intensity. To support object-level identification of PBRs, we additionally integrated Grounding DINO, enabling text-prompted segmentation of features of interest within UAV imagery. This case study demonstrates the effectiveness of our method in generating a 3D semantic map of PBRs, enabling spatial distribution of PBR fragility for earthquake hazard analysis.
Computers & Geosciences↗
Linear complementarity formulation for 3D frictional sliding problems
Frictional sliding on quasi-statically deforming faults and fractures can be modeled efficiently using a linear complementarity formulation. We review the formulation in two dimensions and expand the formulation to three-dimensional problems including problems of orthotropic friction. This formulation accurately reproduces analytical solutions to static Coulomb friction sliding problems. The formulation accounts for opening displacements that can occur near regions of non-planarity even under large confining pressures. Such problems are difficult to solve owing to the coupling of relative displacements and tractions; thus, many geomechanical problems tend to neglect these effects. Simple test cases highlight the importance of including friction and allowing for opening when solving quasi-static fault mechanics models. These results also underscore the importance of considering the effects of non-planarity in modeling processes associated with crustal faulting.
Computational Geosciences↗
Bayesian analysis of the impact of rainfall data product on simulated slope failure for North Carolina locations
In the past decades, many different approaches have been developed in the literature to quantify the load-carrying capacity and geotechnical stability (or the Factor of Safety, F_s) of variably saturated hillslopes. Much of this work has focused on a deterministic characterization of hillslope stability. Yet, simulated F_s values are subject to considerable uncertainty due to our inability to characterize accurately the soil mantle’s properties (hydraulic, geotechnical and geomorphologic) and spatiotemporal variability of the moisture content of the hillslope interior. This is particularly true at larger spatial scales. Thus, uncertainty-incorporating analyses of physically based models of rain-induced landslides are rare in the literature. Such landslide modeling is typically conducted at the hillslope scale using gauge-based rainfall forcing data with rather poor spatiotemporal coverage. For regional landslide modeling, the specific advantages and/or disadvantages of gauge-only, radar-merged and satellite-based rainfall products are not clearly established. Here, we compare and evaluate the performance of the Transient Rainfall Infiltration and Grid-based Regional Slope-stability analysis (TRIGRS) model for three different rainfall products using 112 observed landslides in the period between 2004 and 2011 from the North Carolina Geological Survey database. Our study includes the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis Version 7 (TMPA V7), the North American Land Data Assimilation System Phase 2 (NLDAS-2) analysis, and the reference ‘truth’ Stage IV precipitation. TRIGRS model performance was rather inferior with the use of literature values of the geotechnical parameters and soil hydraulic properties from ROSETTA using soil textural and bulk density data from SSURGO (Soil Survey Geographic database). The performance of TRIGRS improved considerably after Bayesian estimation of the parameters with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm using Stage IV precipitation data. Hereto, we use a likelihood function that combines binary slope failure information from landslide event and ‘null’ periods using multivariate frequency distribution-based metrics such as the False Discovery and False Omission Rates. Our results demonstrate that the Stage IV-inferred TRIGRS parameter distributions generalize well to TMPA and NLDAS-2 precipitation data, particularly at sites with considerably larger TMPA and NLDAS-2 rainfall amounts during landslide events than null periods. TRIGRS model performance is then rather similar for all three rainfall products. At higher elevations, however, the TMPA and NLDAS-2 precipitation volumes are insufficient and their performance with the Stage IV-derived parameter distributions indicate their inability to accurately characterize hillslope stability.
North Carolina↗
Q-mode factor analysis of compositional data
An extended form of the method of Q-mode factor analysis may be used if the rows of the data matrix sum to a constant. The constant row-sum provides a means for scaling the factor model to conform to the units of the original data in such forms as proportions, percentages, or parts per million. The scale factors for an m-factor model are used to derive composition scores and loadings, which lead to the same reproduced data regardless of the choice of reference axes (end-members). Factor-variance diagrams may be constructed to show the proportion of the variance that will be accounted for in each compositional variable by models containing various numbers of end-members. The method may be used to develop a petrologic mixing model by searching for end-members that are geologically plausible. ?? 1976.
Computers & Geosciences↗
Extended cabfac and Qmodel computer programs for Q-mode factor analysis of compositional data
The computer program CABFAC for Q-mode factor analysis of geologic data has been extended for use with data having constant row-sums. Another program, QMODEL, reads an output file from CABFAC and can be used to develop a variety of Q-mode models. The models serve to reproduce estimates of the original data rather than of the data in normalized form. ?? 1976.
Computers & Geosciences↗
Resource data bases-Resource assessment
The U.S. Geological Survey's Office of Resource Analysis is developing computer methods for the handling of mineral-resources data in order to provide improved means for addressing and manipulating data. These methods include: computerized data files and predictive resource models. Data files contain the raw or disaggregated information on mineral deposits and commodities. One operational data file is CRIB (Computerized Resources Information Bank) which is a general purpose inventory and reference file on metallic and nonmetallic mineral deposits. A computer file on resources should contain detailed information on the following main categories: record identification, name and location, description of deposit, analytical data, and production/reserves. A resource model employs postulates and inferences in conjunction with the data to make predictions about resources-as key variables concerning a mineral commodity are changed. The objective is to estimate the availability of minerals including: geological availability (occurrence models), technological availability (exploration and beneficiation models), and economic availability (economics models). ?? 1976.
Computers & Geosciences↗
The national coal-resources data system of the U.S. geological survey
The National Coal Resources Data System (NCRDS) was designed by the U.S. Geological Survey (USGS) to meet the increasing demands for rapid retrieval of information on coal location, quantity, quality, and accessibility. An interactive conversational query system devised by the USGS retrieves information from the data bank through a standard computer terminal. The system is being developed in two phases. Phase I, which currently is available on a limited basis, contains published areal resource and chemical data. The primary objective of this phase is to retrieve, calculate, and tabulate coal-resource data by area on a local, regional, or national scale. Factors available for retrieval include: state, county, quadrangle, township, coal field, coal bed, formation, geologic age, source and reliability of data, and coal-bed rank, thickness, overburden, and tonnage, or any combinations of variables. In addition, the chemical data items include individual values for proximate and ultimate analyses, BTU value, and several other physical and chemical tests. Information will be validated and deleted or updated as needed. Phase II is being developed to store, retrieve, and manipulate basic point source coal data (e.g., field observations, drill-hole logs), including geodetic location; bed thickness; depth of burial; moisture; ash; sulfur; major-, minor-, and trace-element content; heat value; and characteristics of overburden, roof rocks, and floor rocks. The computer system may be used to generate interactively structure-contour or isoline maps of the physical and chemical characteristics of a coal bed or to calculate coal resources. ?? 1976.
Computers & Geosciences↗
RESIN, a FORTRAN IV program for determining the area of influence of samples or drill holes in resource target search
A FORTRAN IV program that calculates the area of influence of drill holes or samples with respect to the size and shape of elliptical or circular resource targets is presented. Program options include determination of the degree to which areas within a region have been explored and estimation of probabilities that points are centers of undiscovered deposits. Errors of recognition can be utilized in the program input. ?? 1976.
Computers & Geosciences↗
Binorm-a fortran subroutine to calculate the percentiles of a standardized binormal distribution
BINORM is a FORTRAN subroutine for calculating the percentiles of a standardized binormal distribution. By using a linear transformation, the percentiles of a binormal distribution can be obtained. The percentiles of a binormal distribution are useful for plotting purposes, for establishing confidence intervals, and for sampling from a mixed population that consists of two normal distributions. ?? 1977.
Computers & Geosciences↗
The U.S. geological survey rass-statpac system for management and statistical reduction of geochemical data
RASS is an acronym for Rock Analysis Storage System and STATPAC, for Statistical Package. The RASS and STATPAC computer programs are integrated into the RASS-STATPAC system for the management and statistical reduction of geochemical data. The system, in its present form, has been in use for more than 9 yr by scores of U.S. Geological Survey geologists, geochemists, and other scientists engaged in a broad range of geologic and geochemical investigations. The principal advantage of the system is the flexibility afforded the user both in data searches and retrievals and in the manner of statistical treatment of data. The statistical programs provide for most types of statistical reduction normally used in geochemistry and petrology, but also contain bridges to other program systems for statistical processing and automatic plotting. ?? 1977.
Computers & Geosciences↗
Graphic analysis of resources by numerical evaluation techniques (Garnet)
An interactive computer program for graphical analysis has been developed by the U.S. Geological Survey. The program embodies five goals, (1) economical use of computer resources, (2) simplicity for user applications, (3) interactive on-line use, (4) minimal core requirements, and (5) portability. It is designed to aid (1) the rapid analysis of point-located data, (2) structural mapping, and (3) estimation of area resources. ?? 1977.
Computers & Geosciences↗
A method and fortran program for the computerized reconstruction of three-dimensional objects from serial sections
Three-dimensional objects can be reconstructed automatically from sets of serial sections by using a suitable algorithm. For two-phase systems it is sufficient to match phase boundary loops between adjacent sections. When all possible neighboring pairs of sections have been compared in this manner, the bounding surfaces of the object are described by sets of triangular and nonplanar quadrilateral patches. ?? 1977.
Computers & Geosciences↗
An exploratory drilling exhaustion sequence plot program
The exhaustion sequence plot program computes the conditional area of influence for wells in a specified rectangular region with respect to a fixed-size deposit. The deposit is represented by an ellipse whose size is chosen by the user. The area of influence may be displayed on computer printer plots consisting of a maximum of 10,000 grid points. At each point, a symbol is presented that indicates the probability of that point being exhausted by nearby wells with respect to a fixed-size ellipse. This output gives a pictorial view of the manner in which oil fields are exhausted. In addition, the exhaustion data may be used to estimate the number of deposits remaining in a basin. ?? 1977.
Computers & Geosciences↗