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 37 records · Page 2Linked to original sources
USGS1992

Mass-conserving method of characteristics for streamflow modeling

A robust numerical model is presented for the computation of unsteady streamflow on steep river slopes. The one-dimensional model uses the method of characteristics on a specified space-time grid to solve the Saint-Venant equations. An additional continuity equation requirement on each space-time element provides greatly improved conservation of mass over traditional implementations of the method of characteristics on a fixed grid. The space-time geometry of the problem is described in a finite element setting. Hermite interpolation of channel parameters is used to avoid numerical difficulties that may occur with steep slopes due to discontinuities in the derivatives of data such as channel top width. Manning's equation for friction slope can be modified by a factor to make the slope more appropriate for steep rivers. The standard Manning's friction slope can also be used, if preferred. The computer model is not restricted to steep slopes, and applies as well to gently sloping streams. Two numerical examples support the mathematical approach and computational algorithm.

Water Supply Paper↗
USGS1977

Image analysis as a method for determining grain-size distributions in sandstones

A Quantimet 720 image analyzing computer is useful for A Quantimet 720 image analyzing computer is useful for determination of grain-size distributions of sand-size material. Results including histograms and statistical measures of size can be produced in approximately 5 to 10 minutes per sample. Comparison with sieve results shows only small differences in the mean grain-sizes obtained. Research into operating methods for the image analyzing computer and into mathematical methods for combining data is continuing so that material coarser and finer than sand can be measured. The instrument, which is attached to a microscope via a black and white television scanner, uses the gray shades of features in the field of view as the means of detecting those features. The image analyzing computer can then measure the features including area, perimeter, and several types of diameters. These data are then analyzed by a computer and the results are plotted. Sample preparation for this technique is conventional. Samples can be run as loose disaggregated grains, grain mounts, or thin sections.

Open-File Report↗
USGS2000

Geomorphometry-Diversity in quantitative surface analysis

A widening variety of applications is diversifying geomorphometry ( digital terrain modelling ), the quantitative study of topography. An amalgam of earth science, mathematics, engineering and computer science, the discipline has been revolutionized by the computer manipulation of gridded terrain heights, or digital elevation models (DEMs). Its rapid expansion continues. This article reviews the remarkable diversity of recent morphometric work in 15 selected topics and discusses their significance and prospects. The quantitative analysis of industrial microsurface topography is introduced to the earth science community. The 14 other topics are Internet access to geomorphometry; global DEMs; DEM modelling of channel networks; self-organized criticality; fractal and wavelet analysis; soil resources; landslide hazards; barchan dunes; harvesting wind energy; sea-ice surfaces; sea-floor abyssal hills; Japanese work in morphometry; and the emerging fields of landscape ecology and image understanding. Closing remarks note reasons for the diversity within geomorphometry, speculate on future trends and recommend creating a unified field of surface representation.

Progress in Physical Geography: Earth and Environm↗
USGS2023

BatTool: Projecting bat populations facing multiple stressors using a demographic model

Bats provide ecologically and agriculturally important ecosystem services but are currently experiencing population declines caused by multiple environmental stressors, including mortality from white-nose syndrome and wind energy development. Analyses of the current and future health and viability of these species may support conservation management decision making. Demographic modeling provides a quantitative tool for decision makers and conservation managers to make more informed decisions, but widespread adoption of these tools can be limited because of the complexity of the mathematical, statistical, and computational components involved in implementing these models. In this work, we provide an exposition of the BatTool R package, detailing the primary components of the matrix projection model, a publicly accessible graphical user interface ( https://rconnect.usgs.gov/battool ) facilitating user-defined scenario analyses, and its intended uses and limitations (Wiens et al., US Geol Surv Data Release 2022; Wiens et al., US Geol Surv Softw Release 2022). We present a case study involving wind energy permitting, weighing the effects of potential mortality caused by a hypothetical wind energy facility on the projected abundance of four imperiled bat species in the Midwestern United States.

Methods in Ecology and Evolution↗
USGS2002

A bibliography of terrain modeling (geomorphometry), the quantitative representation of topography: Supplement 4.0

Terrain modeling, the practice of ground-surface quantification, is an amalgam of Earth science, mathematics, engineering, and computer science. The discipline is known variously as geomorphometry (or simply morphometry), terrain analysis, and quantitative geomorphology. It continues to grow through myriad applications to hydrology, geohazards mapping, tectonics, sea-floor and planetary exploration, and other fields. Dating nominally to the co-founders of academic geography, Alexander von Humboldt (1808, 1817) and Carl Ritter (1826, 1828), the field was revolutionized late in the 20th Century by the computer manipulation of spatial arrays of terrain heights, or digital elevation models (DEMs), which can quantify and portray ground-surface form over large areas (Maune, 2001). Morphometric procedures are implemented routinely by commercial geographic information systems (GIS) as well as specialized software (Harvey and Eash, 1996; Köthe and others, 1996; ESRI, 1997; Drzewiecki et al., 1999; Dikau and Saurer, 1999; Djokic and Maidment, 2000; Wilson and Gallant, 2000; Breuer, 2001; Guth, 2001; Eastman, 2002). The new Earth Surface edition of the Journal of Geophysical Research, specializing in surficial processes, is the latest of many publication venues for terrain modeling. This is the fourth update of a bibliography and introduction to terrain modeling (Pike, 1993, 1995, 1996, 1999) designed to collect the diverse, scattered literature on surface measurement as a resource for the research community. The use of DEMs in science and technology continues to accelerate and diversify (Pike, 2000a). New work appears so frequently that a sampling must suffice to represent the vast literature. This report adds 1636 entries to the 4374 in the four earlier publications1. Forty-eight additional entries correct dead Internet links and other errors found in the prior listings. Chronicling the history of terrain modeling, many entries in this report predate the 1999 supplement. Coverage is representative from about 1800 through early–mid 2002. Papers increasingly are published exclusively or in duplicate on the Internet's World Wide Web; the dates given here for Web addresses (URLs) that lack a print publication indicate a Web site's last update or my last access of it. The bibliography is arranged alphabetically and thus is not readily summarized. This introduction cites about 500 entries, a third of them grouped under 24 morphometric topics, as a guide to the listing's contents. Continuing the practice of previous bibliographies in the series to provide more information on a few applications (see summary of past topics in Pike, 2000a), this report elaborates further on topographic data, putative new parameters, tectonic geomorphology/neo-orometry, biogeography, ice-cap morphometry, results from the Mars Global DEM, landslide-hazard mapping, terrain modeling as physics, Hack's law, and broad-scale computer visualization. The literature of some of these subjects is large, and none of the summaries is intended to more than introduce the topic and comment on some of the current contributions of terrain modeling. Closing the essay is a discussion of pre-1900 papers that trace the evolution of ridge-line and watercourse quantification by descriptive geometry, as well as comments on some new books and an on-line bulletin board.

Open-File Report↗
USGS1980

Planning and design of studies for river-quality assessment in the Truckee and Carson River basins, California and Nevada

The objectives of the Geological Survey 's river-quality assessment in the Truckee and Carson River basins in California and Nevada are to identify the significant resource management problems; to develop techniques to assess the problems; and to effectively communicate results to responsible managers. Six major elements of the assessment to be completed by October 1981 are (1) a detailing of the legal, institutional, and structural development of water resources in the basins and the current problems and conflicts; (2) a compilation and synthesis of the physical hydrology of the basins; (3) development of a special workshop approach to involve local management in the direction and results of the study; (4) development of a comprehensive streamflow model emcompassing both basins to provide a quantitative hydrologic framework for water-quality analysis; (5) development of a water-quality transport model for selected constituents and characteristics on selected reaches of the Truckee River; and (6) a detailed examination of selected fish habitats for specified reaches of the Truckee River. Progress will be periodically reported in reports, maps, computer data files, mathematical models, a bibliography, and public presentations. In building a basic framework to develop techniques, the basins were viewed as a single hydrologic unit because of interconnecting diversion structures. The framework comprises 13 hydrographic subunits to facilitate modeling and sampling. Several significant issues beyond the scope of the assessment were considered as supplementary proposals; water-quality loadings in Truckee and Carson Rivers, urban runoff in Reno and management alternatives, and a model of limnological processes in Lahontan Reservoir. (USGS)

California, Nevada↗
USGS1968

Water-discharge determinations for the tidal reach of the Willamette River from Ross Island Bridge to Mile 10.3, Portland, Oregon

Water-discharge, velocity, and slope variations for a 3.7-mile-Iong tidal reach of the Willamette River at Portland, Oreg., were defined from discharge measurements and river stage data collected between July 1962 and January 1965. Observed water discharge during tide-affected flows, during floods, and during backwater from the Columbia River and recorded stages at each end of the river reach were used to determine water discharge from two mathematical models. These models use a finite-difference method to solve the equations of moderately unsteady open-channel streamflow, and discharges are computed by an electronic digital computer. Discharges computed by using the mathematical models compare satisfactorily with observed discharges, except during the period of backwater from the annual flood of the Columbia River. The flow resistance coefficients used in the models vary with discharge; for one model, the coefficients for discharges above 30,000 cfs (cubic feet per second) are 12 and 24 percent less than the coefficient used for discharges below 30,000 cfs. Daily mean discharges were determined by use of one mathematical model for approximately two-thirds of the water year, October 1963 through September 1964. Agreement of computed with routed daily mean discharges is fair; above 30,000 cfs, average differences between the two discharges are about 10 percent, and below 30,000 cfs, computed daily discharges are consistently greater (by as much as 25 percent) than routed discharges. The other model was used to compute discharges for the unusually high flood flows of December 1964.

Water Supply Paper↗
USGS2013

Conditional spectrum computation incorporating multiple causal earthquakes and ground-motion prediction models

The conditional spectrum (CS) is a target spectrum (with conditional mean and conditional standard deviation) that links seismic hazard information with ground-motion selection for nonlinear dynamic analysis. Probabilistic seismic hazard analysis (PSHA) estimates the ground-motion hazard by incorporating the aleatory uncertainties in all earthquake scenarios and resulting ground motions, as well as the epistemic uncertainties in ground-motion prediction models (GMPMs) and seismic source models. Typical CS calculations to date are produced for a single earthquake scenario using a single GMPM, but more precise use requires consideration of at least multiple causal earthquakes and multiple GMPMs that are often considered in a PSHA computation. This paper presents the mathematics underlying these more precise CS calculations. Despite requiring more effort to compute than approximate calculations using a single causal earthquake and GMPM, the proposed approach produces an exact output that has a theoretical basis. To demonstrate the results of this approach and compare the exact and approximate calculations, several example calculations are performed for real sites in the western United States. The results also provide some insights regarding the circumstances under which approximate results are likely to closely match more exact results. To facilitate these more precise calculations for real applications, the exact CS calculations can now be performed for real sites in the United States using new deaggregation features in the U.S. Geological Survey hazard mapping tools. Details regarding this implementation are discussed in this paper.

Bulletin of the Seismological Society of America↗
USGS2000

Information technology developments within the national biological information infrastructure

Looking out an office window or exploring a community park, one can easily see the tremendous challenges that biological information presents the computer science community. Biological information varies in format and content depending whether or not it is information pertaining to a particular species (i.e. Brown Tree Snake), or a specific ecosystem, which often includes multiple species, land use characteristics, and geospatially referenced information. The complexity and uniqueness of each individual species or ecosystem do not easily lend themselves to today's computer science tools and applications. To address the challenges that the biological enterprise presents the National Biological Information Infrastructure (NBII) (http://www.nbii.gov) was established in 1993. The NBII is designed to address these issues on a National scale within the United States, and through international partnerships abroad. This paper discusses current computer science efforts within the National Biological Information Infrastructure Program and future computer science research endeavors that are needed to address the ever-growing issues related to our Nation's biological concerns.

Mathematics and Computers in Modern Science - Acou↗
USGS1987

AQMAN; linear and quadratic programming matrix generator using two-dimensional ground-water flow simulation for aquifer management modeling

A FORTRAN-77 computer program code that helps solve a variety of aquifer management problems involving the control of groundwater hydraulics. It is intended for use with any standard mathematical programming package that uses Mathematical Programming System input format. The computer program creates the input files to be used by the optimization program. These files contain all the hydrologic information and management objectives needed to solve the management problem. Used in conjunction with a mathematical programming code, the computer program identifies the pumping or recharge strategy that achieves a user 's management objective while maintaining groundwater hydraulic conditions within desired limits. The objective may be linear or quadratic, and may involve the minimization of pumping and recharge rates or of variable pumping costs. The problem may contain constraints on groundwater heads, gradients, and velocities for a complex, transient hydrologic system. Linear superposition of solutions to the transient, two-dimensional groundwater flow equation is used by the computer program in conjunction with the response matrix optimization method. A unit stress is applied at each decision well and transient responses at all control locations are computed using a modified version of the U.S. Geological Survey two dimensional aquifer simulation model. The program also computes discounted cost coefficients for the objective function and accounts for transient aquifer conditions. (Author 's abstract)

Water-Resources Investigations Report↗
USGS1988

A simple approach to nonlinear estimation of physical systems

Recursive algorithms for estimating the states of nonlinear physical systems are developed. This requires some key hypotheses regarding the structure of the underlying processes. Members of this class of random processes have several desirable properties for the nonlinear estimation of random signals. An assumption is made about the form of the estimator, which may then take account of a wide range of applications. Under the above assumption, the estimation algorithm is mathematically suboptimal but effective and computationally attractive. It may be compared favorably to Taylor series-type filters, nonlinear filters which approximate the probability density by Edgeworth or Gram-Charlier series, as well as to conventional statistical linearization-type estimators. To link theory with practice, some numerical results for a simulated system are presented, in which the responses from the proposed and the extended Kalman algorithms are compared.

Mathematical and Computer Modelling↗
USGS2020

Carrying capacity of a population diffusing in a heterogeneous environment

The carrying capacity of the environment for a population is one of the key concepts in ecology and it is incorporated in the growth term of reaction-diffusion equations describing populations in space. Analysis of reaction-diffusion models of populations in heterogeneous space have shown that, when the maximum growth rate and carrying capacity in a logistic growth function vary in space, conditions exist for which the total population size at equilibrium (i) exceeds the total population that which would occur in the absence of diffusion and (ii) exceeds that which would occur if the system were homogeneous and the total carrying capacity, computed as the integral over the local carrying capacities, was the same in the heterogeneous and homogeneous cases. We review here work over the past few years that has explained these apparently counter-intuitive results in terms of the way input of energy or another limiting resource (e.g., a nutrient) varies across the system. We report on both mathematical analysis and laboratory experiments confirming that total population size in a heterogeneous system with diffusion can exceed that in the system without diffusion. We further report, however, that when the resource of the population in question is explicitly modeled as a coupled variable, as in a reaction-diffusion chemostat model rather than a model with logistic growth, the total population in the heterogeneous system with diffusion cannot exceed the total population size in the corresponding homogeneous system in which the total carrying capacities are the same.

Mathematics and Computers in Modern Science - Acou↗
USGS2024

Bayesian approaches to proxy uncertainty quantification in paleoecology: A mathematical justification and practical integration

Paleoenvironmental data are essential for reconstructing environmental conditions in the distant past, and these reconstructions strongly depend on proxies and age–depth models. Proxies are indirect measurements that substitute for variables that cannot be directly measured, such as past precipitation. Conversely, an age–depth model is a tool that correlates the observed proxy with a specific moment in time. Bayesian age–depth modelling has proved to be a powerful method for estimating sediment ages and their associated uncertainties. However, there remains considerable potential for further integration into proxy analysis. In this paper, we explore a mathematical justification and a computational approach that integrates uncertainty at the age–depth level and propagates it to the proxy scale in the form of a posterior predictive distribution. This method mitigates potential biases and errors by removing the need to assign a single age to a given proxy measurement. It allows for quantifying the likelihood that proxy data values correspond to modelled ages, thus enabling the quantification of uncertainty in both the temporal and proxy value domains. The use of Bayesian statistics in proxy analysis represents a relatively recent advancement. We aim to mathematically justify incorporating the Markov chain Monte Carlo output from age–depth models into proxy analysis and to present a novel methodology for constructing environmental reconstructions using this approach.

Journal of Agricultural, Biological and Environmen↗
USGS2005

Parallel Fortran-MPI software for numerical inversion of the Laplace transform and its application to oscillatory water levels in groundwater environments

A parallel Fortran-MPI (Message Passing Interface) software for numerical inversion of the Laplace transform based on a Fourier series method is developed to meet the need of solving intensive computational problems involving oscillatory water level's response to hydraulic tests in a groundwater environment. The software is a parallel version of ACM (The Association for Computing Machinery) Transactions on Mathematical Software (TOMS) Algorithm 796. Running 38 test examples indicated that implementation of MPI techniques with distributed memory architecture speedups the processing and improves the efficiency. Applications to oscillatory water levels in a well during aquifer tests are presented to illustrate how this package can be applied to solve complicated environmental problems involved in differential and integral equations. The package is free and is easy to use for people with little or no previous experience in using MPI but who wish to get off to a quick start in parallel computing. ?? 2004 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software↗
USGS2012

Use of gene-expression programming to estimate Manning’s roughness coefficient for high gradient streams

Manning’s roughness coefficient ( n ) has been widely used in the estimation of flood discharges or depths of flow in natural channels. Therefore, the selection of appropriate Manning’s n values is of paramount importance for hydraulic engineers and hydrologists and requires considerable experience, although extensive guidelines are available. Generally, the largest source of error in post-flood estimates (termed indirect measurements) is due to estimates of Manning’s n values, particularly when there has been minimal field verification of flow resistance. This emphasizes the need to improve methods for estimating n values. The objective of this study was to develop a soft computing model in the estimation of the Manning’s n values using 75 discharge measurements on 21 high gradient streams in Colorado, USA. The data are from high gradient (S > 0.002 m/m), cobble- and boulder-bed streams for within bank flows. This study presents Gene-Expression Programming (GEP), an extension of Genetic Programming (GP), as an improved approach to estimate Manning’s roughness coefficient for high gradient streams. This study uses field data and assessed the potential of gene-expression programming (GEP) to estimate Manning’s n values. GEP is a search technique that automatically simplifies genetic programs during an evolutionary processes (or evolves) to obtain the most robust computer program (e.g., simplify mathematical expressions, decision trees, polynomial constructs, and logical expressions). Field measurements collected by Jarrett (J Hydraulic Eng ASCE 110: 1519–1539, 1984 ) were used to train the GEP network and evolve programs. The developed network and evolved programs were validated by using observations that were not involved in training. GEP and ANN-RBF (artificial neural network-radial basis function) models were found to be substantially more effective (e.g., R 2 for testing/validation of GEP and RBF-ANN is 0.745 and 0.65, respectively) than Jarrett’s (J Hydraulic Eng ASCE 110: 1519–1539, 1984 ) equation (R 2 for testing/validation equals 0.58) in predicting the Manning’s n .

Water Resources Management↗
USGS1982

An electric-analog simulation of elliptic partial differential equations using finite element theory

Elliptic partial differential equations can be solved using the Galerkin-finite element method to generate the approximating algebraic equations, and an electrical network to solve the resulting matrices. Some element configurations require the use of networks containing negative resistances which, while physically realizable, are more expensive and time-consuming to construct.

Mathematics and Computers in Simulation↗
USGS1987

The space transformation in the simulation of multidimensional random fields

Space transformations are proposed as a mathematically meaningful and practically comprehensive approach to simulate multidimensional random fields. Within this context the turning bands method of simulation is reconsidered and improved in both the space and frequency domains.

Mathematics and Computers in Simulation↗