USGS Science⌕ Search

Geology topics

Kenneth L. Kipp

Publications and source records attributed to Kenneth L. Kipp.

11 recordsLinked to original sources

A geochemical transport model for redox-controlled movement of mineral fronts in groundwater flow systems: A case of nitrate removal by oxidation of pyrite

A one-dimensional prototype geochemical transport model was developed in order to handle simultaneous precipitation-dissolution and oxidation-reduction reactions governed by chemical equilibria. Total aqueous component concentrations are the primary dependent variables, and a sequential iterative approach is used for the calculation. The model was verified by analytical and numerical comparisons and is able to simulate sharp mineral fronts. At a site in Denmark, denitrification has been observed by oxidation of pyrite. Simulation of nitrate movement at this site showed a redox front movement rate of 0.58 m yr −1 , which agreed with calculations of others. It appears that the sequential iterative approach is the most practical for extension to multidimensional simulation and for handling large numbers of components and reactions. However, slow convergence may limit the size of redox systems that can be handled.

Water Resources Research↗

Adaptation of the Carter-Tracy water influx calculation to groundwater flow simulation

The Carter-Tracy calculation for water influx is adapted to groundwater flow simulation with additional clarifying explanation not present in the original papers. The Van Everdingen and Hurst aquifer-influence functions for radial flow from an outer aquifer region are employed. This technique, based on convolution of unit-step response functions, offers a simple but approximate method for embedding an inner region of groundwater flow simulation within a much larger aquifer region where flow can be treated in an approximate fashion. The use of aquifer-influence functions in groundwater flow modeling reduces the size of the computational grid with a corresponding reduction in computer storage and execution time. The Carter-Tracy approximation to the convolution integral enables the aquifer influence function calculation to be made with an additional storage requirement of only two times the number of boundary nodes more than that required for the inner region simulation. It is a good approximation for constant flow rates but is poor for time-varying flow rates where the variation is large relative to the mean. A variety of outer aquifer region geometries, exterior boundary conditions, and flow rate versus potentiometric head relations can be used. The radial, transient-flow case presented is representative. An analytical approximation to the functions of Van Everdingen and Hurst for the dimensionless potentiometric head versus dimensionless time is given.

Water Resources Research↗

Groundwater transport of strontium 90 in a glacial outwash environment

As part of the investigation of groundwater contamination at a uranium-scrap recovery plant at Wood River Junction, Rhode Island, laboratory experiments led to the development of a model for predicting the transport of strontium 90 in glacial outwash sediments based on an approximate mechanism for ion exchange. The multicomponent system was simplified to two components by regarding all exchangeable cations other than strontium 90 as a single component. The binary ion-exchange parameter was a function of the variable, total ion concentration. A one-dimensional solute transport model was formulated to evaluate the time necessary for natural groundwater flow to remove the strontium 90 contamination plume from the groundwater system to the Pawcatuck River. The finite difference transport equations were solved sequentially for total ion concentrations, then strontium 90 concentrations. Clay-free quartz and feldspar sands at the study site have little potential for strontium 90 sorption, and high calcium, magnesium, and sodium concentrations compete for the few ion exchange sites. As the total ion concentration plume moves out of the system, ion exchange of strontium 90 increases, reducing the strontium 90 concentration in the groundwater. Cleanout times predicted using the binary ion exchange mechanism were about two thirds of those predicted using a constant distribution coefficient. It is suggested that this type of model can simulate solute transport more realistically in many groundwater systems where the total ion concentration is not constant.

Rhode Island↗

PHAST version 2-A program for simulating groundwater flow, solute transport, and multicomponent geochemical reactions

The computer program PHAST (PHREEQC And HST3D) simulates multicomponent, reactive solute transport in three-dimensional saturated groundwater flow systems. PHAST is a versatile groundwater flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. Major enhancements in PHAST Version 2 allow spatial data to be defined in a combination of map and grid coordinate systems, independent of a specific model grid (without node-by-node input). At run time, aquifer properties are interpolated from the spatial data to the model grid; regridding requires only redefinition of the grid without modification of the spatial data. PHAST is applicable to the study of natural and contaminated groundwater systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock/water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, or density-dependent flow. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux (specified-flux), and leaky (head-dependent) conditions, as well as the special cases of rivers, drains, and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association or Pitzer specific interaction thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, ion exchange sites, surface complexation sites, solid solutions, and gases; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, exchangers, surfaces, gases, kinetic reactants, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and postprocessing programs; and in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST.

Techniques and Methods↗

Type curve analysis of inertial effects in the response of a well to a slug test.

The water level response to a slug or bailer test in a well completed in a confined aquifer has been evaluated taking into account well-bore storage and inertial effects of the water column in the well. The response range, from overdamped with negligible inertial effects to damped oscillation, was covered employing numerical inversions of the Laplace-transform solution. By scaling the time with respect to the undamped natural period of the well-aquifer system and by using the damping parameter for a second-order damped, inertial-elastic system, a set of type curves was constructed that enables water level response data from a slug or bailer test to be analyzed under conditions where the inertial parameter is large. Values of transmissivity and effective static water column length can be determined when an estimate of storage coefficient is available. The numerical solution and resulting type curves cover the transition range between the limiting cases of negligible inertial effects and of damped oscillation that have been treated by others. Two examples of slug test analysis show that precise results depend on accurate measurements of water level displacement (±1% of initial value).

Water Resources Research↗

Guide to the Revised Ground-Water Flow and Heat Transport Simulator: HYDROTHERM - Version 3

The HYDROTHERM computer program simulates multi-phase ground-water flow and associated thermal energy transport in three dimensions. It can handle high fluid pressures, up to 1 ? 109 pascals (104 atmospheres), and high temperatures, up to 1,200 degrees Celsius. This report documents the release of Version 3, which includes various additions, modifications, and corrections that have been made to the original simulator. Primary changes to the simulator include: (1) the ability to simulate unconfined ground-water flow, (2) a precipitation-recharge boundary condition, (3) a seepage-surface boundary condition at the land surface, (4) the removal of the limitation that a specified-pressure boundary also have a specified temperature, (5) a new iterative solver for the linear equations based on a generalized minimum-residual method, (6) the ability to use time- or depth-dependent functions for permeability, (7) the conversion of the program code to Fortran 90 to employ dynamic allocation of arrays, and (8) the incorporation of a graphical user interface (GUI) for input and output. The graphical user interface has been developed for defining a simulation, running the HYDROTHERM simulator interactively, and displaying the results. The combination of the graphical user interface and the HYDROTHERM simulator forms the HYDROTHERM INTERACTIVE (HTI) program. HTI can be used for two-dimensional simulations only. New features in Version 3 of the HYDROTHERM simulator have been verified using four test problems. Three problems come from the published literature and one problem was simulated by another partially saturated flow and thermal transport simulator. The test problems include: transient partially saturated vertical infiltration, transient one-dimensional horizontal infiltration, two-dimensional steady-state drainage with a seepage surface, and two-dimensional drainage with coupled heat transport. An example application to a hypothetical stratovolcano system with unconfined ground-water flow is presented in detail. It illustrates the use of HTI with the combination precipitation-recharge and seepage-surface boundary condition, and functions as a tutorial example problem for the new user.

Techniques and Methods↗

PHAST--a program for simulating ground-water flow, solute transport, and multicomponent geochemical reactions

The computer program PHAST simulates multi-component, reactive solute transport in three-dimensional saturated ground-water flow systems. PHAST is a versatile ground-water flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. PHAST is applicable to the study of natural and contaminated ground-water systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock-water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, density-dependent flow, or waters with high ionic strengths. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux, and leaky conditions, as well as the special cases of rivers and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, gases, surface complexation sites, ion exchange sites, and solid solutions; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, gases, exchangers, surfaces, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and post-processing programs; or in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST. Operator splitting of the flow, transport, and geochemical equations is used to separate the three processes into three sequential calculations. No iterations between transport and reaction calculations are implemented. A three-dimensional Cartesian coordinate system and finite-difference techniques are used for the spatial and temporal discretization of the flow and transport equations. The non-linear chemical equilibrium equations are solved by a Newton-Raphson method, and the kinetic reaction equations are solved by a Runge-Kutta or an implicit method for integrating ordinary differential equations. The PHAST simulator may require large amounts of memory and long Central Processing Unit (CPU) times. To reduce the long CPU times, a parallel version of PHAST has been developed that runs on a multiprocessor computer or on a collection of computers that are networked. The parallel version requires Message Passing Interface, which is currently (2004) freely available. The parallel version is effective in reducing simulation times. This report documents the use of the PHAST simulator, including running the simulator, preparing the input files, selecting the output files, and visualizing the results. It also presents four examples that verify the numerical method and demonstrate the capabilities of the simulator. PHAST requires three input files. Only the flow and transport file is described in detail in this report. The other two files, the chemistry data file and the database file, are identical to PHREEQC files and the detailed description of these files is found in the PHREEQC documentation.

Techniques and Methods↗

Hydrology of the Bonneville Salt Flats, northwestern Utah, and simulation of ground-water flow and solute transport in the shallow-brine aquifer

This report describes the hydrologic system of the Bonneville Salt Flats with emphasis on the mechanisms of solute transport. Variable-density, three-dimensional computer simulations of the near-surface part of the ground-water system were done to quantify both the transport of salt dissolved in subsurface brine that leaves the salt-crust area and the salt dissolved and precipitated on the land surface. The study was designed to define the hydrology of the brine ground-water system and the natural and anthropogenic processes causing salt loss, and where feasible, to quantify these processes. Specific areas of study include the transport of salt in solution by ground-water flow and the transport of salt in solution by wind-driven ponds and the subsequent salt precipitation on the surface of the playa upon evaporation or seepage into the subsurface. In addition, hydraulic and chemical changes in the hydrologic system since previous studies were documented.

Utah↗

Investigation of salt loss from the Bonneville Salt Flats, northwestern Utah

The Bonneville Salt Flats study area is located in the western part of the Great Salt Lake Desert in northwestern Utah, about 110 miles west of Salt Lake City. The salt crust covers about 50 square miles, but the extent varies yearly as a result of salt being dissolved by the formation and movement of surface ponds during the winter and redeposited with the evaporation of these ponds during the summer. A decrease in thickness and extent of the salt crust on the Bonneville Salt Flats has been documented during 1960-88 (S. Brooks, Bureau of Land Management, written commun., 1989). Maximum salt-crust thickness was 7 feet in 1960 and 5.5 feet in 1988. No definitive data are available to identify and quantify the processes that cause salt loss. More than 55 million tons of salt are estimated to have been lost from the salt crust during the 28-year period. The Bureau of Land Management needs to know the causes of salt loss to make appropriate management decisions.

Utah↗

Ground-water flow and contaminant transport at a radioactive-materials processing site, Wood River Junction, Rhode Island

Liquid wastes from an enriched-uranium cold-scrap recovery plant at Wood River Junction, Rhode Island, were discharged to the environment through evaporation ponds and trenches from 1966 through 1980. Leakage from the ponds and trenches resulted in a plume of contaminated ground water extending northwestward to the Pawcatuck River through a highly permeable sand and gravel aquifer of glacial origin.

Rhode Island↗