Search USGS⌕ Search

SEARCH · Search USGS

Results for “Mathematical and Computer Modelling”

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 127 records · Page 7Linked to original sources
USGS2011

A distribution-based parameterization for improved tomographic imaging of solute plumes

Difference geophysical tomography (e.g. radar, resistivity and seismic) is used increasingly for imaging fluid flow and mass transport associated with natural and engineered hydrologic phenomena, including tracer experiments, in situ remediation and aquifer storage and recovery. Tomographic data are collected over time, inverted and differenced against a background image to produce ‘snapshots’ revealing changes to the system; these snapshots readily provide qualitative information on the location and morphology of plumes of injected tracer, remedial amendment or stored water. In principle, geometric moments (i.e. total mass, centres of mass, spread, etc.) calculated from difference tomograms can provide further quantitative insight into the rates of advection, dispersion and mass transfer; however, recent work has shown that moments calculated from tomograms are commonly biased, as they are strongly affected by the subjective choice of regularization criteria. Conventional approaches to regularization (Tikhonov) and parametrization (image pixels) result in tomograms which are subject to artefacts such as smearing or pixel estimates taking on the sign opposite to that expected for the plume under study. Here, we demonstrate a novel parametrization for imaging plumes associated with hydrologic phenomena. Capitalizing on the mathematical analogy between moment-based descriptors of plumes and the moment-based parameters of probability distributions, we design an inverse problem that (1) is overdetermined and computationally efficient because the image is described by only a few parameters, (2) produces tomograms consistent with expected plume behaviour (e.g. changes of one sign relative to the background image), (3) yields parameter estimates that are readily interpreted for plume morphology and offer direct insight into hydrologic processes and (4) requires comparatively few data to achieve reasonable model estimates. We demonstrate the approach in a series of numerical examples based on straight-ray difference-attenuation radar monitoring of the transport of an ionic tracer, and show that the methodology outlined here is particularly effective when limited data are available.

Geophysical Journal International↗
USGS1977

Ground-water appraisal of the Pineland Sands area, central Minnesota

The Pineland Sands area consists of 770 square miles of surficial glacial outwash, which is undergoing increasing ground-water development for irrigation. The aquifer material is commonly very fine sand to fine gravel, and grain size generally increases from south to north. Thickness, transmissivity, and theoretical well yields are highest in the northern part. In places, saturated thickness exceeds 100 feet, transmissivity exceeds 20,000 feet squared per day, and expected well yields exceed 2,000 gallons per minute. Water in the aquifer is chemically suitable for irrigation. It is calcium bicarbonate water generally containing dissolvedsolids concentrations of 100-300 milligrams per liter. An estimated long-term water budget for the aquifer involves inflow and outflow of 295 cubic feet per second. Over 80 percent of recharge to the aquifer is directly from snowmelt and spring rain; the remainder is from underflow across the northern boundary of the aquifer and locally from streams and lakes. Over 90 percent of discharge from the aquifer is flow to streams and lakes; the remainder is through evapotranspiration. A mathematical model of the surficial aquifer was made to facilitate analysis of the hydrologic system. The model was calibrated through simulation of the natural ground-water flow system. Subsequent applications of the model to estimate effects of future ground-water development indicate that much of the system can safely support long-term, large-scale withdrawals. Results of model analysis show that present development (withdrawals totaling 3.3 cubic feet per second) has no significant effect on the aquifer system. Simulations of hypothetical withdrawals of 60 to 120 cubic feet per second resulted in computed water-table declines as great as 12 feet in places. Most pumpage is derived from intercepted base flow to streams, thus reducing streamflow. Similarly, some lake levels can be expected to decline in response to nearby intensive development.

Minnesota↗
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↗
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↗
USGS2023

Source contributions to suspended sediment and particulate selenium export from the Loutsenhizer Arroyo and Sunflower Drain watersheds in Colorado

Selenium in aquatic ecosystems of the lower Gunnison River Basin in Colorado is affecting the recovery of populations of endangered, native fish species. Dietary exposure is the primary pathway for bioaccumulation of selenium in fish, and particulate selenium can be consumed directly by fish or by the invertebrates on which fish feed. Although selenium can be incorporated into particulate matter via biogeochemical processes, particulate selenium can also enter aquatic ecosystems of the lower Gunnison River Basin from sediments derived from the selenium-rich Mancos Shale. The U.S. Geological Survey, in cooperation with the Colorado Water Conservation Board, conducted this study during 2018–19 to identify sources of selenium-rich suspended sediments from two watersheds underlain by Mancos Shale: Loutsenhizer Arroyo and Sunflower Drain, which is a locally known agricultural drainage near the municipality of Delta, Colorado. A multipronged approach (fieldwork, laboratory work, and computer modeling) referred to as “sediment fingerprinting” was used to evaluate sources of suspended sediments in the streams flowing out of the two studied watersheds. Four potential source types for suspended sediments were identified and sampled (using soil plugs) within the watersheds: rangelands, agricultural fields, arroyo walls, and streambanks. The sediment fingerprinting approach used elemental concentrations and naturally occurring fallout radionuclides as tracers to apportion percent contributions from the four source types of suspended sediments found in streamflow from both watersheds. To determine the dominant sources of suspended sediment in streamflow from both watersheds, a mathematical “unmixing” model was used. Unmixing models apportion source percentages to samples of material in which those sources are mixed. These models used elemental and isotopic data in the suspended sediments to unmix them into proportional contributions from source types. The results indicated that arroyo walls and streambanks generally dominated as sources of the suspended sediment. Arroyo walls and streambanks were channel-adjacent sources, with sediments mobilized by water flowing within the stream channel. These sources accounted for greater than 50 percent of suspended sediment in all but one sample and accounted for 100 percent of suspended sediment in 5 of the 11 samples collected. Rangeland and agricultural field sources were located in uplands outside of stream channels and were detected more often during the non-irrigation season. Rangeland and agricultural field sources each were found in 5 of the 11 samples collected. Concentrations of selenium in sediment-source samples were comparatively greater in streambanks and lower in rangelands, with agricultural fields and arroyo walls being intermediate. As a result, source apportionments for particulate selenium skewed towards sources adjacent to stream channels more than for suspended sediments. Water imports for irrigation have changed the hydrology of the watersheds, and a notable fraction of imported water passes through the watersheds rapidly. The rapid flowthrough water during the irrigation season likely contributes heavily to sediment erosion and transport in Loutsenhizer Arroyo and Sunflower Drain, particularly from channel-adjacent sources of sediment. Decreases in irrigation season streamflow, at least in Loutsenhizer Arroyo, may have decreased sediment erosion and transport during the 2018–20 irrigation seasons compared to the 2015–17 seasons.

Colorado↗
USGS1975

Optimal exploitation strategies for an animal population in a Markovian environment: A theory and an example

Optimal exploitation strategies were studied for an animal population in a Markovian (stochastic, serially correlated) environment. This is a general case and encompasses a number of important special cases as simplifications. Extensive empirical data on the Mallard (Anas platyrhynchos) were used as an example of general theory. The number of small ponds on the central breeding grounds was used as an index to the state of the environment. A general mathematical model was formulated to provide a synthesis of the existing literature, estimates of parameters developed from an analysis of data, and hypotheses regarding the specific effect of exploitation on total survival. The literature and analysis of data were inconclusive concerning the effect of exploitation on survival. Therefore, two hypotheses were explored: (1) exploitation mortality represents a largely additive form of mortality, and (2) exploitation mortality is compensatory with other forms of mortality, at least to some threshold level. Models incorporating these two hypotheses were formulated as stochastic dynamic programming models and optimal exploitation strategies were derived numerically on a digital computer. Optimal exploitation strategies were found to exist under the rather general conditions. Direct feedback control was an integral component in the optimal decision—making process. Optimal exploitation was found to be substantially different depending upon the hypothesis regarding the effect of exploitation on the population. If we assume that exploitation is largely an additive force of mortality in Mallards, then optimal exploitation decisions are a convex function of the size of the breeding population and a linear or slight concave function of the environmental conditions. Under the hypothesis of compensatory mortality forces, optimal exploitation decisions are approximately linearly related to the size of the Mallard breeding population. Dynamic programming is suggested as a very general formulation for realistic solutions to the general optimal exploitation problem. The concepts of state vectors and stage transformations are completely general. Populations can be modeled stochastically and the objective function can include extra—biological factors. The optimal level of exploitation in year t must be based on the observed size of the population and the state of the environment in year t unless the dynamics of the population, the state of the environment, and the result of the exploitation decisions are completely deterministic. Exploitation based on an average harvest, or harvest rate, or designed to maintain a constant breeding population size is inefficient.

Ecology↗
USGS2015

PRMS-IV, the precipitation-runoff modeling system, version 4

Computer models that simulate the hydrologic cycle at a watershed scale facilitate assessment of variability in climate, biota, geology, and human activities on water availability and flow. This report describes an updated version of the Precipitation-Runoff Modeling System. The Precipitation-Runoff Modeling System is a deterministic, distributed-parameter, physical-process-based modeling system developed to evaluate the response of various combinations of climate and land use on streamflow and general watershed hydrology. Several new model components were developed, and all existing components were updated, to enhance performance and supportability. This report describes the history, application, concepts, organization, and mathematical formulation of the Precipitation-Runoff Modeling System and its model components. This updated version provides improvements in (1) system flexibility for integrated science, (2) verification of conservation of water during simulation, (3) methods for spatial distribution of climate boundary conditions, and (4) methods for simulation of soil-water flow and storage.

Techniques and Methods↗
USGS1974

Optimal exploitation strategies for an animal population in a stochastic serially correlated environment

Optimal exploitation strategies were studied for an animal population in a stochastic, serially correlated environment. This is a general case and encompasses a number of important cases as simplifications. Data on the mallard (Anas platyrhynchos) were used to explore the exploitation strategies and test several hypotheses because relatively much is known concerning the life history and general ecology of this species and extensive empirical data are available for analysis. The number of small ponds on the central breeding grounds was used as an index to the state of the environment. Desirable properties of an optimal exploitation strategy were defined. A mathematical model was formulated to provide a synthesis of the existing literature, estimates of parameters developed from an analysis of data, and hypotheses regarding the specific effect of exploitation on total survival. Both the literature and the analysis of data were inconclusive concerning the effect of exploitation on survival. Therefore, alternative hypotheses were formulated: (1) exploitation mortality represents a largely additive form of mortality, or (2 ) exploitation mortality is compensatory with other forms of mortality, at least to some threshold level. Models incorporating these two hypotheses were formulated as stochastic dynamic programming models and optimal exploitation strategies were derived numerically on a digital computer. Optimal exploitation strategies were found to exist under rather general conditions. Direct feedback control was an integral component in the optimal decision-making process. Optimal exploitation was found to be substantially different depending upon the hypothesis regarding the effect of exploitation on the population. Assuming that exploitation is largely an additive force of mortality, optimal exploitation decisions are a convex function of the size of the breeding population and a linear or slightly concave function of the environmental conditions. Optimal exploitation under this hypothesis tends to reduce the variance of the size of the population. Under the hypothesis of compensatory mortality forces, optimal exploitation decisions are approximately linearly related to the size of the breeding population. Environmental variables may be somewhat more important than the size of the breeding population to the production of young mallards. In contrast, the size of the breeding population appears to be more important in the exploitation process than is the state of the environment. The form of the exploitation strategy appears to be relatively insensitive to small changes in the production rate. In general, the relative importance of the size of the breeding population may decrease as fecundity increases. The optimal level of exploitation in year t must be based on the observed size of the population and the state of the environment in year t unless the dynamics of the population, the state of the environment, and the result of the exploitation decisions are completely deterministic. Exploitation based on an average harvest, harvest rate, or designed to maintain a constant breeding population size is inefficient.

Book↗
USGS1974

Optimal exploitation strategies for an animal population in a stochastic serially correlated environment

Optimal exploitation strategies were studied for an animal population in a stochastic, serially correlated environment. This is a general case and encompasses a number of important cases as simplifications. Data on the mallard (Anas platyrhynchos) were used to explore the exploitation strategies and test several hypotheses because relatively much is known concerning the life history and general ecology of this species and extensive empirical data are available for analysis. The number of small ponds on the central breeding grounds was used as an index to the state of the environment. Desirable properties of an optimal exploitation strategy were defined. A mathematical model was formulated to provide a synthesis of the existing literature, estimates of parameters developed from an analysis of data, and hypotheses regarding the specific effect of exploitation on total survival. Both the literature and the analysis of data were inconclusive concerning the effect of exploitation on survival. Therefore, alternative hypotheses were formulated: (1) exploitation mortality represents a largely additive form of mortality, or (2 ) exploitation mortality is compensatory with other forms of mortality, at least to some threshold level. Models incorporating these two hypotheses were formulated as stochastic dynamic programming models and optimal exploitation strategies were derived numerically on a digital computer. Optimal exploitation strategies were found to exist under rather general conditions. Direct feedback control was an integral component in the optimal decision-making process. Optimal exploitation was found to be substantially different depending upon the hypothesis regarding the effect of exploitation on the population. Assuming that exploitation is largely an additive force of mortality, optimal exploitation decisions are a convex function of the size of the breeding population and a linear or slightly concave function of the environmental conditions. Optimal exploitation under this hypothesis tends to reduce the variance of the size of the population. Under the hypothesis of compensatory mortality forces, optimal exploitation decisions are approximately linearly related to the size of the breeding population. Environmental variables may be somewhat more important than the size of the breeding population to the production of young mallards. In contrast, the size of the breeding population appears to be more important in the exploitation process than is the state of the environment. The form of the exploitation strategy appears to be relatively insensitive to small changes in the production rate. In general, the relative importance of the size of the breeding population may decrease as fecundity increases. The optimal level of exploitation in year t must be based on the observed size of the population and the state of the environment in year t unless the dynamics of the population, the state of the environment, and the result of the exploitation decisions are completely deterministic. Exploitation based on an average harvest, harvest rate, or designed to maintain a constant breeding population size is inefficient.

Dissertation Abstracts International↗
USGS2024

Wave propagation in layered soil deposits

Recent advances in the general theory of viscoelastic waves and rays in layered media provide a rigorous mathematical framework for site-specific, soil-response models used for earthquake resistant design. The advances provide general closed-form anelastic solutions for the classic problems of the response of a stack of soil layers to S and P waves, ray theory for reflected and refracted waves, Rayleigh- and Love-Type surface waves, and head waves. These general solutions valid for anelastic media regardless of the amount of material damping yield new insights regarding the characteristics of seismic waves and their ray paths that are not provided by conventional models. They provide corresponding numerical ground-response models and ray-tracing computation algorithms that account for changes in velocity and attenuation of anelastic waves associated with changes in inhomogeneity of the waves induced by anelastic soil and soil-rock boundaries. These anelastic effects manifest themselves as variations in amplitude response, amplitude attenuation, raypath location, and travel time as observed at the Earth’s surface. Implications of these anelastic effects for soil-response models used for earthquake resistant design are provided herein.

Conference Paper↗
USGS2023

Calibration of imperfect geophysical models by multiple satellite interferograms with measurement bias

Model calibration consists of using experimental or field data to estimate the unknown parameters of a mathematical model. The presence of model discrepancy and measurement bias in the data complicates this task. Satellite interferograms, for instance, are widely used for calibrating geophysical models in geological hazard quantification. In this work, we used satellite interferograms to relate ground deformation observations to the properties of the magma chamber at K i¯lauea Volcano in Hawai‘i. We derived closed-form marginal likelihoods and implemented posterior sampling procedures that simultaneously estimate the model discrepancy of physical models, and the measurement bias from the atmospheric error in satellite interferograms. We found that model calibration by aggregating multiple interferograms and downsampling the pixels in the interferograms can reduce the computation complexity compared to calibration approaches based on multiple data sets. The conditions that lead to no loss of information from data aggregation and downsampling are studied. Simulation illustrates that both discrepancy and measurement bias can be estimated, and real applications demonstrate that modeling both effects helps obtain a reliable estimation of a physical model’s unobserved parameters and enhance its predictive accuracy. We implement the computational tools in the RobustCalibration package available on CRAN.

Technometrics↗
USGS2008

Steady and intermittent slipping in a model of landslide motion regulated by pore-pressure feedback

This paper studies a parsimonious model of landslide motion, which consists of the one-dimensional diffusion equation (for pore pressure) coupled through a boundary condition to a first-order ODE (Newton's second law). Velocity weakening of sliding friction gives rise to nonlinearity in the model. Analysis shows that solutions of the model equations exhibit a subcritical Hopf bifurcation in which stable, steady sliding can transition to cyclical, stick-slip motion. Numerical computations confirm the analytical predictions of the parameter values at which bifurcation occurs. The existence of stick-slip behavior in part of the parameter space is particularly noteworthy because, unlike stick-slip behavior in classical models , here it arises in the absence of a reversible (elastic) driving force. Instead, the driving force is static (gravitational), mediated by the effects of pore-pressure diffusion on frictional resistance.

SIAM Journal on Applied Mathematics↗
USGS2006

Forecasting runout of rock and debris avalanches

Physically based mathematical models and statistically based empirical equations each may provide useful means of forecasting runout of rock and debris avalanches. This paper compares the foundations, strengths, and limitations of a physically based model and a statistically based forecasting method, both of which were developed to predict runout across three-dimensional topography. The chief advantage of the physically based model results from its ties to physical conservation laws and well-tested axioms of soil and rock mechanics, such as the Coulomb friction rule and effective-stress principle. The output of this model provides detailed information about the dynamics of avalanche runout, at the expense of high demands for accurate input data, numerical computation, and experimental testing. In comparison, the statistical method requires relatively modest computation and no input data except identification of prospective avalanche source areas and a range of postulated avalanche volumes. Like the physically based model, the statistical method yields maps of predicted runout, but it provides no information on runout dynamics. Although the two methods differ significantly in their structure and objectives, insights gained from one method can aid refinement of the other.

Conference Paper↗
USGS2014

A suggestion for computing objective function in model calibration

A parameter-optimization process (model calibration) is usually required for numerical model applications, which involves the use of an objective function to determine the model cost (model-data errors). The sum of square errors (SSR) has been widely adopted as the objective function in various optimization procedures. However, ‘square error’ calculation was found to be more sensitive to extreme or high values. Thus, we proposed that the sum of absolute errors (SAR) may be a better option than SSR for model calibration. To test this hypothesis, we used two case studies—a hydrological model calibration and a biogeochemical model calibration—to investigate the behavior of a group of potential objective functions: SSR, SAR, sum of squared relative deviation (SSRD), and sum of absolute relative deviation (SARD). Mathematical evaluation of model performance demonstrates that ‘absolute error’ (SAR and SARD) are superior to ‘square error’ (SSR and SSRD) in calculating objective function for model calibration, and SAR behaved the best (with the least error and highest efficiency). This study suggests that SSR might be overly used in real applications, and SAR may be a reasonable choice in common optimization implementations without emphasizing either high or low values (e.g., modeling for supporting resources management).

Ecological Informatics↗
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↗
USGS1981

A transient laboratory method for determining the hydraulic properties of 'tight' rocks-I. Theory

Transient pulse testing has been employed increasingly in the laboratory to measure the hydraulic properties of rock samples with low permeability. Several investigators have proposed a mathematical model in terms of an initial-boundary value problem to describe fluid flow in a transient pulse test. However, the solution of this problem has not been available. In analyzing data from the transient pulse test, previous investigators have either employed analytical solutions that are derived with the use of additional, restrictive assumptions, or have resorted to numerical methods. In Part I of this paper, a general, analytical solution for the transient pulse test is presented. This solution is graphically illustrated by plots of dimensionless variables for several cases of interest. The solution is shown to contain, as limiting cases, the more restrictive analytical solutions that the previous investigators have derived. A method of computing both the permeability and specific storage of the test sample from experimental data will be presented in Part II.

International Journal of Rock Mechanics and Mining↗
USGS1972

Hydrograph simulation models of the Hillsborough and Alafia Rivers, Florida: a preliminary report

Mathematical (digital) models that simulate flood hydrographs from rainfall records have been developed for the following gaging stations in the Hillsborough and Alafia River basins of west-central Florida: Hillsborough River near Tampa, Alafia River at Lithia, and north Prong Alafia River near Keysville. These models, which were developed from historical streamflow and and rainfall records, are based on rainfall-runoff and unit-hydrograph procedures involving an arbitrary separation of the flood hydrograph. These models assume the flood hydrograph to be composed of only two flow components, direct (storm) runoff, and base flow. Expressions describing these two flow components are derived from streamflow and rainfall records and are combined analytically to form algorithms (models), which are programmed for processing on a digital computing system. Most Hillsborough and Alafia River flood discharges can be simulated with expected relative errors less than or equal to 30 percent and flood peaks can be simulated with average relative errors less than 15 percent. Because of the inadequate rainfall network that is used in obtaining input data for the North Prong Alafia River model, simulated peaks are frequently in error by more than 40 percent, particularly for storms having highly variable areal rainfall distribution. Simulation errors are the result of rainfall sample errors and, to a lesser extent, model inadequacy. Data errors associated with the determination of mean basin precipitation are the result of the small number and poor areal distribution of rainfall stations available for use in the study. Model inadequacy, however, is attributed to the basic underlying theory, particularly the rainfall-runoff relation. These models broaden and enhance existing water-management capabilities within these basins by allowing the establishment and implementation of programs providing for continued development in these areas. Specifically, the models serve not only as a basis for forecasting floods, but also for simulating hydrologic information needed in flood-plain mapping and delineating and evaluating alternative flood control and abatement plans.

Florida↗
USGS1987

Documentation of a deep percolation model for estimating ground-water recharge

A deep percolation model, which operates on a daily basis, was developed to estimate long-term average groundwater recharge from precipitation. It has been designed primarily to simulate recharge in large areas with variable weather, soils, and land uses, but it can also be used at any scale. The physical and mathematical concepts of the deep percolation model, its subroutines and data requirements, and input data sequence and formats are documented. The physical processes simulated are soil moisture accumulation, evaporation from bare soil, plant transpiration, surface water runoff, snow accumulation and melt, and accumulation and evaporation of intercepted precipitation. The minimum data sets for the operation of the model are daily values of precipitation and maximum and minimum air temperature, soil thickness and available water capacity, soil texture, and land use. Long-term average annual precipitation, actual daily stream discharge, monthly estimates of base flow, Soil Conservation Service surface runoff curve numbers, land surface altitude-slope-aspect, and temperature lapse rates are optional. The program is written in the FORTRAN 77 language with no enhancements and should run on most computer systems without modifications. Documentation has been prepared so that program modifications may be made for inclusions of additional physical processes or deletion of ones not considered important. (Author 's abstract)

Open-File Report↗