Search USGS⌕ Search

SEARCH · Search USGS

Results for “Mathematics of Computation”

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
USGS2007

OPR-PPR, a computer program for assessing data importance to model predictions using linear statistics

The OPR-PPR program calculates the Observation-Prediction (OPR) and Parameter-Prediction (PPR) statistics that can be used to evaluate the relative importance of various kinds of data to simulated predictions. The data considered fall into three categories: (1) existing observations, (2) potential observations, and (3) potential information about parameters. The first two are addressed by the OPR statistic; the third is addressed by the PPR statistic. The statistics are based on linear theory and measure the leverage of the data, which depends on the location, the type, and possibly the time of the data being considered. For example, in a ground-water system the type of data might be a head measurement at a particular location and time. As a measure of leverage, the statistics do not take into account the value of the measurement. As linear measures, the OPR and PPR statistics require minimal computational effort once sensitivities have been calculated. Sensitivities need to be calculated for only one set of parameter values; commonly these are the values estimated through model calibration. OPR-PPR can calculate the OPR and PPR statistics for any mathematical model that produces the necessary OPR-PPR input files. In this report, OPR-PPR capabilities are presented in the context of using the ground-water model MODFLOW-2000 and the universal inverse program UCODE_2005. The method used to calculate the OPR and PPR statistics is based on the linear equation for prediction standard deviation. Using sensitivities and other information, OPR-PPR calculates (a) the percent increase in the prediction standard deviation that results when one or more existing observations are omitted from the calibration data set; (b) the percent decrease in the prediction standard deviation that results when one or more potential observations are added to the calibration data set; or (c) the percent decrease in the prediction standard deviation that results when potential information on one or more parameters is added.

Techniques and Methods↗
USGS1977

State of stress in the lithosphere: Inferences from the flow laws of olivine

The experimental flow data for rocks and minerals are reviewed and found to fit a law of the form {Mathematical expression} where {Mathematical expression} This law reduces to the familiar power-law stress dependency at low stress and to an exponential stress dependency at high stress. Using the material flow law parameters for olivine, stress profiles with depth and strain rate are computed for a representative range of temperature distributions in the lithosphere. The results show that the upper 15 to 25 km of the oceanic lithosphere must behave elastically or fail by fracture and that the remainder deforms by exponential law flow at intermediate depths and by power-law flow in the rest. A model computation of the gravitational sliding of a lithospheric plate using olivine rheology exhibits a very sharp decoupling zone which is a consequence of the combined effects of increasing stress and temperature on the flow law, which is a very sensitive function of both. ?? 1977 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS2012

Slope-Area Computation Program Graphical User Interface 1.0—A Preprocessing and Postprocessing Tool for Estimating Peak Flood Discharge Using the Slope-Area Method

The slope-area method is a technique for estimating the peak discharge of a flood after the water has receded (Dalrymple and Benson, 1967). This type of discharge estimate is called an “indirect measurement” because it relies on evidence left behind by the flood, such as high-water marks (HWMs) on trees or buildings. These indicators of flood stage are combined with measurements of the cross-sectional geometry of the stream, estimates of channel roughness, and a mathematical model that balances the total energy of the flow between cross sections. This is in contrast to a “direct” measurement of discharge during the flood where cross-sectional area is measured and a current meter or acoustic equipment is used to measure the water velocity. When a direct discharge measurement cannot be made at a gage during high flows because of logistics or safety reasons, an indirect measurement of a peak discharge is useful for defining the high-flow section of the stage-discharge relation (rating curve) at the stream gage, resulting in more accurate computation of high flows. The Slope-Area Computation program (SAC; Fulford, 1994) is an implementation of the slope-area method that computes a peak-discharge estimate from inputs of water-surface slope (from surveyed HWMs), channel geometry, and estimated channel roughness. SAC is a command line program written in Fortran that reads input data from a formatted text file and prints results to another formatted text file. Preparing the input file can be time-consuming and prone to errors. This document describes the SAC graphical user interface (GUI), a crossplatform “wrapper” application that prepares the SAC input file, executes the program, and helps the user interpret the output. The SAC GUI is an update and enhancement of the slope-area method (SAM; Hortness, 2004; Berenbrock, 1996), an earlier spreadsheet tool used to aid field personnel in the completion of a slope-area measurement. The SAC GUI reads survey data, develops a plan-view plot, water-surface profile, cross-section plots, and develops the SAC input file. The SAC GUI also develops HEC-2 files that can be imported into HEC–RAS.

Fact Sheet↗
USGS2017

A general modeling framework for describing spatially structured population dynamics

Variation in movement across time and space fundamentally shapes the abundance and distribution of populations. Although a variety of approaches model structured population dynamics, they are limited to specific types of spatially structured populations and lack a unifying framework. Here, we propose a unified network-based framework sufficiently novel in its flexibility to capture a wide variety of spatiotemporal processes including metapopulations and a range of migratory patterns. It can accommodate different kinds of age structures, forms of population growth, dispersal, nomadism and migration, and alternative life-history strategies. Our objective was to link three general elements common to all spatially structured populations (space, time and movement) under a single mathematical framework. To do this, we adopt a network modeling approach. The spatial structure of a population is represented by a weighted and directed network. Each node and each edge has a set of attributes which vary through time. The dynamics of our network-based population is modeled with discrete time steps. Using both theoretical and real-world examples, we show how common elements recur across species with disparate movement strategies and how they can be combined under a unified mathematical framework. We illustrate how metapopulations, various migratory patterns, and nomadism can be represented with this modeling approach. We also apply our network-based framework to four organisms spanning a wide range of life histories, movement patterns, and carrying capacities. General computer code to implement our framework is provided, which can be applied to almost any spatially structured population. This framework contributes to our theoretical understanding of population dynamics and has practical management applications, including understanding the impact of perturbations on population size, distribution, and movement patterns. By working within a common framework, there is less chance that comparative analyses are colored by model details rather than general principles

Ecology and Evolution↗
USGS1980

Temperature and solute-transport simulation in streamflow using a Lagrangian reference frame

A computer program for simulating one-dimensional, unsteady temperature and solute transport in a river has been developed and documented for general use. The solution approach to the convective-diffusion equation uses a moving reference frame (Lagrangian) which greatly simplifies the mathematics of the solution procedure and dramatically reduces errors caused by numerical dispersion. The solution procedure has the further advantages, relative to conventional Eulerian solution schemes, of being easy to understand in the physical sense, of being extremely stable numerically, and of providing an accounting system which is very useful for model calibration. The model documentation is presented as a series of four programs of increasing complexity. The conservative transport model can be used to route a single conservative substance, such as dye, through a reach of a river. The simplified temperature model is used to predict water temperature in rivers, either with or without thermal loading, when few meteorological data are available. Only equilibrium temperature and windspeed are required. It is suggested that air temperature can be used to approximate equilibrium temperature. The complete temperature model is highly accurate but requires rather complete meteorological data. Finally, the 10-parameter model can be used to route as many as 10 interacting constituents through a river reach. The mathematical description of the interaction between the constituents, which does not need to be linear, is generally up to the user to supply. An example problem is solved for a three-parameter system involving temperature, dissolved oxygen, and biochemical oxygen demand. For simplicity, all models are developed and presented assuming steady non-uniform flow. Generalization of the models to allow unsteady flow is extremely simple, involving the addition of no more than 18 cards to the program deck. The report is concluded by describing this generalization for any of the models. Before using the models with unsteady flow, a flow model must be used to calculate and store the necessary flow data at each cross section and time step. Such a flow model is available and documented.

Water-Resources Investigations Report↗
USGS1989

A vertically averaged spectral model for tidal circulation in estuaries

A frequency dependent computer model based on the two-dimensional vertically averaged shallow-water equations is described for general purpose application in tidally dominated embayments. This model simulates the response of both tides and tidal currents to user-specified geometries and boundary conditions. The mathematical formulation and practical application of the model are discussed in detail. Salient features of the model include the ability to specify: (1) stage at the open boundaries as well as within the model grid, (2) velocities on open boundaries (river inflows and so forth), (3) spatially variable wind stress, and (4) spatially variable bottom friction. Using harmonically analyzed field data as boundary conditions, this model can be used to make real time predictions of tides and tidal currents. (USGS)

Water-Resources Investigations Report↗
USGS1971

Reduction of lithologic-log data to numbers for use in the digital computer

The development of a standardized system for conveniently coding lithologic-log data for use in the digital computer has long been needed. The technique suggested involves a reduction of the original written alphanumeric log to a numeric log by use of computer programs. This numeric log can then be retrieved as a written log, interrogated for pertinent information, or analyzed statistically.

Journal of the International Association for Mathe↗
USGS2011

A two-phase debris-flow model that includes coupled evolution of volume fractions, granular dilatancy, and pore-fluid pressure

Pore-fluid pressure plays a crucial role in debris flows because it counteracts normal stresses at grain contacts and thereby reduces intergranular friction. Pore-pressure feedback accompanying debris deformation is particularly important during the onset of debrisflow motion, when it can dramatically influence the balance of forces governing downslope acceleration. We consider further effects of this feedback by formulating a new, depth-averaged mathematical model that simulates coupled evolution of granular dilatancy, solid and fluid volume fractions, pore-fluid pressure, and flow depth and velocity during all stages of debris-flow motion. To illustrate implications of the model, we use a finite-volume method to compute one-dimensional motion of a debris flow descending a rigid, uniformly inclined slope, and we compare model predictions with data obtained in large-scale experiments at the USGS debris-flow flume. Predictions for the first 1 s of motion show that increasing pore pressures (due to debris contraction) cause liquefaction that enhances flow acceleration. As acceleration continues, however, debris dilation causes dissipation of pore pressures, and this dissipation helps stabilize debris-flow motion. Our numerical predictions of this process match experimental data reasonably well, but predictions might be improved by accounting for the effects of grain-size segregation.

Conference Paper↗
USGS2021

The application of ensemble wave forcing to quantify uncertainty of shoreline change predictions

Reliable predictions and accompanying uncertainty estimates of coastal evolution on decadal to centennial time scales are increasingly sought. So far, most coastal change projections rely on a single, deterministic realization of the unknown future wave climate, often derived from a global climate model. Yet, deterministic projections do not account for the stochastic nature of future wave conditions across a variety of temporal scales (e.g., daily, weekly, seasonally, and interannually). Here, we present an ensemble Kalman filter shoreline change model to predict coastal erosion and uncertainty due to waves at a variety of time scales. We compare shoreline change projections, simulated with and without ensemble wave forcing conditions by applying ensemble wave time series produced by a computationally efficient statistical downscaling method. We demonstrate a sizable (site-dependent) increase in model uncertainty compared with the unrealistic case of model projections based on a single, deterministic realization (e.g., a single time series) of the wave forcing. We support model-derived uncertainty estimates with a novel mathematical analysis of ensembles of idealized process models. Here, the developed ensemble modeling approach is applied to a well-monitored beach in Tairua, New Zealand. However, the model and uncertainty quantification techniques derived here are generally applicable to a variety of coastal settings around the world.

JGR Earth Surface↗
USGS2022

Absolute accuracy assessment of lidar point cloud using amorphous objects

The accuracy assessment of airborne lidar point cloud typically estimates vertical accuracy by computing RMSEz (root mean square error of the z coordinate) from ground check points (GCPs). Due to the low point density of the airborne lidar point cloud, there is often not enough accurate semantic context to find an accurate conjugate point. To advance the accuracy assessment in full three-dimensional (3D) context, geometric features, such as the three-plane intersection point or two-line intersection point, are often used. Although the point density is still low, geometric features are mathematically modeled from many points. Thus, geometric features provide a robust determination of the intersection point, and the point is considered as a GCP. When no regular built objects are available, we describe the process of utilizing features of irregular shape called amorphous natural objects, such as a tree or a rock. When scanned to a high-density point cloud, an amorphous natural object can be used as ground truth reference data to estimate 3D georeferencing errors of the airborne lidar point cloud. The algorithm to estimate 3D accuracy is the optimization that minimizes the sum of the distance between the airborne lidar points to the ground scanned data. The search volume partitioning was the most important procedure to improve the computational efficiency. We also performed an extensive study to address the external uncertainty associated with the amorphous object method. We describe an accuracy assessment using amorphous objects (108 trees) spread over the project area. The accuracy results for ∆ x , ∆ y , and ∆ z obtained using the amorphous object method were 3.1 cm, 3.6 cm, and 1.7 cm RMSE, along with a mean error of 0.1 cm, 0.1 cm, and 4.5 cm, respectively, satisfying the accuracy requirement of U.S. Geological Survey lidar base specification. This approach shows strong promise as an alternative to geometric feature methods when artificial targets are scarce. The relative convenience and advantages of using amorphous targets, along with its good performance shown here, make this amorphous object method a practical way to perform 3D accuracy assessment.

Remote Sensing↗
USGS2000

Uncertainty estimation for resource assessment-an application to coal

The U.S. Geological Survey is conducting a national assessment of coal resources. As part of that assessment, a geostatistical procedure has been developed to estimate the uncertainty of coal resources for the historical categories of geological assurance: measured, indicated, inferred, and hypothetical coal. Data consist of spatially clustered coal thickness measurements from coal beds and/or zones that cover, in some cases, several thousand square kilometers. Our procedure involved trend removal, an examination of spatial correlation, computation of a sample semivariogram, and fitting a semivariogram model. This model provided standard deviations for the uncertainty estimates. The number of sample points (drill holes) in each historical category also was estimated. Measurement error in the thickness of the coal bed/zone was obtained from the fitted model or supplied exogenously. From this information approximate estimates of uncertainty on the historical categories were computed. We illustrate the methodology using drill hole data from the Harmon coal bed located in southwestern North Dakota. The methodology will be applied to approximately 50 coal data sets.

Mathematical Geology↗
USGS2010

Bayesian Inference: with ecological applications

This text provides a mathematically rigorous yet accessible and engaging introduction to Bayesian inference with relevant examples that will be of interest to biologists working in the fields of ecology, wildlife management and environmental studies as well as students in advanced undergraduate statistics.. This text opens the door to Bayesian inference, taking advantage of modern computational efficiencies and easily accessible software to evaluate complex hierarchical models.

Book↗
USGS1993

Relationship of geological and geothermal field properties: Midcontinent area, USA, an example

Quantitative approaches to data analysis in the last decade have become important in basin modeling and mineral-resource estimation. The interrelation of geological, geophysical, geochemical, and geohydrological variables is important in adjusting a model to a real-world situation. Revealing the interdependences of variables can contribute in understanding the processes interacting in sedimentary basins. It is reasonably simple to compare spatial data of the same type but more difficult if different properties are involved. Statistical techniques, such as cluster analysis or principal components analysis, or some algebraic approaches can be used to ascertain the relations of standardized spatial data. In this example, structural configuration on five different stratigraphic horizons, one total sediment thickness map, and four maps of geothermal data were copared. As expected, the structural maps are highly related because all had undergone about the same deformation with differing degrees of intensity. The temperature gradients derived (1) from shallow borehole logging measurements under equilibrium conditions with the surrounding rock, and (2) from non-equilibrium bottom-hole temperatures (BHT) from deeper depths are mainly independent of each other. This was expected and confirmed also for the two temperature maps at 1000 ft which were constructed using both types of gradient values. Thus, it is evident that the use of a 2-point (BHT and surface temperature) straightline calculation of a mean temperature gradient gives different information about the geothermal regime than using gradients from temperatures logged under equilibrium conditions. Nevertheless, it is useful to determine to what a degree the larger dataset of nonequilibrium temperatures could reflect quantitative relationships to geologic conditions. Comparing all maps of geothermal information vs. the structural and the sediment thickness maps, it was determined that all correlations are moderately negative or slightly positive. These results are clearly shown by the cluster analysis and the principal components. Considering a close relationship between temperature and thermal conductivity of the sediments as observed for most of the Midcontinent area and relatively homogeneous heat-flow density conditions for the study area these results support the following assumptions: (1) undifferentiated geothermal gradients, computed from temperatures of different depth intervals and differing sediment properties, cannot contribute to an improved understanding of the temperature structure and its controls within the sedimentary cover, and (2) the quantitative approach of revealing such relations needs refined datasets of temperature information valid for the different depth levels or stratigraphic units. ?? 1993 International Association for Mathematical Geology.

Mathematical Geology↗
USGS1987

Estimating threshold values for the land disposal of organic solvent-contaminated wastes

This paper describes conceptual and mathematical modeling approaches to evaluate environmentally suitable waste management practices that prevent ground water pollution from landfill disposal of solid wastes contaminated with organic solvents. Threshold values for disposal of solid wastes contaminated with acetone, nitrobenzene, chlorobenzene, and carbon tetrachloride were estimated using a disposal site model using a two-dimensional advection-,dispersion solute transport computer code. For example, the threshold value for disposal of carbon tetrachloride at the site model was approximately 600 g/acre, while the value for acetone was limited only by the characteristics of the waste to prevent free drainage of liquid solvent. The study indicated that the toxicity, mass loading rate, and the mass of solvent initially disposed are critical parameters determining the relative success of a given site to attenuate contaminants to environmentally acceptable levels.

Journal of Hazardous Materials↗
USGS2011

Partnership of Environmental Education and Research-A compilation of student research, 1999-2008

The U.S. Geological Survey (USGS) Tennessee Water Science Center and the College of Engineering and Technology at Tennessee State University developed a Partnership in Environmental Education and Research (PEER) to support environmental research at TSU and to expand the environmental research capabilities of the USGS in Tennessee. The PEER program is driven by the research needs to better define the occurrence, fate, and transport of contaminants in groundwater and surface water. Research in the PEER program has primarily focused on the transport and remediation of organic contamination in karst settings. Research conducted through the program has also expanded to a variety of media and settings. Research areas include contaminant occurrence and transport, natural and enhanced bioremediation, geochemical conditions in karst aquifers, mathematical modeling for contaminant transport and degradation, new methods to evaluate groundwater contamination, the resuspension of bacteria from sediment in streams, the use of bioluminescence and chemiluminescence to identify the presence of contaminants, and contaminant remediation in wetlands. The PEER program has increased research and education opportunities for students in the College of Engineering, Technology, and Computer Science and has provided students with experience in presenting the results of their research. Students in the program have participated in state, regional, national and international conferences with more than 140 presentations since 1998 and more than 40 student awards. The PEER program also supports TSU outreach activities and efforts to increase minority participation in environmental and earth science programs at the undergraduate and graduate levels. TSU students and USGS staff participate in the TSU summer programs for elementary and high school students to promote earth sciences. The 2007 summer camps included more than 130 students from 20 different States and Washington DC.

Tennessee↗
USGS2011

Computational considerations for collecting and using data in the equidistant cylindrical map projection and the bounds of sampling geographic data at progressively higher resolution

The Equidistant Cylindrical Map projection is popular with digital modelers and others for storing and processing worldwide data sets because of the simple association of latitude and longitude to cell values or pixels in the resulting grid. This projection does not accurately display area, and the diminished geographic area represented by cells at high latitudes is not often carefully considered. A simple mathematical analysis quantifies the discrepancy in area sampled by cells at different latitudes. The presence of this discrepancy indicates that the use of this projection can induce bias in data sets when both sampling and reporting data. It is demonstrated that as the resolution requirements of input data for models increase, the necessity of providing data to accurately describe smaller cells, particularly at high latitude, will be a challenge.

Techniques and Methods↗
USGS1988

The geometric signature: Quantifying landslide-terrain types from digital elevation models

Topography of various types and scales can be fingerprinted by computer analysis of altitude matrices (digital elevation models, or DEMs). The critical analytic tool is the geometric signature, a set of measures that describes topographic form well enough to distinguish among geomorphically disparate landscapes. Different surficial processes create topography with diagnostic forms that are recognizable in the field. The geometric signature abstracts those forms from contour maps or their DEMs and expresses them numerically. This multivariate characterization enables once-in-tractable problems to be addressed. The measures that constitute a geometric signature express different but complementary attributes of topographic form. Most parameters used here are statistical estimates of central tendency and dispersion for five major categories of terrain geometry; altitude, altitude variance spectrum, slope between slope reversals, and slope and its curvature at fixed slope lengths. As an experimental application of geometric signatures, two mapped terrain types associated with different processes of shallow landsliding in Marin County, California, were distinguished consistently by a 17-variable description of topography from 21??21 DEMs (30-m grid spacing). The small matrix is a statistical window that can be used to scan large DEMs by computer, thus potentially automating the mapping of contrasting terrain types. The two types in Marin County host either (1) slow slides: earth flows and slump-earth flows, or (2) rapid flows: debris avalanches and debris flows. The signature approach should adapt to terrain taxonomy and mapping in other areas, where conditions differ from those in Central California. ?? 1988 International Association for Mathematical Geology.

Mathematical Geology↗
USGS1988

Use of the truncated shifted Pareto distribution in assessing size distribution of oil and gas fields

The truncated shifted Pareto (TSP) distribution, a variant of the two-parameter Pareto distribution, in which one parameter is added to shift the distribution right and left and the right-hand side is truncated, is used to model size distributions of oil and gas fields for resource assessment. Assumptions about limits to the left-hand and right-hand side reduce the number of parameters to two. The TSP distribution has advantages over the more customary lognormal distribution because it has a simple analytic expression, allowing exact computation of several statistics of interest, has a "J-shape," and has more flexibility in the thickness of the right-hand tail. Oil field sizes from the Minnelusa play in the Powder River Basin, Wyoming and Montana, are used as a case study. Probability plotting procedures allow easy visualization of the fit and help the assessment. ?? 1988 International Association for Mathematical Geology.

Mathematical Geology↗