Search USGS⌕ Search

Geology topics

K.L. Kipp

Publications and source records attributed to K.L. Kipp.

10 recordsLinked to original sources
USGS2008

Source and transport controls on the movement of nitrate to public supply wells in selected principal aquifers of the United States

In 2003–2005, systematic studies in four contrasting hydrogeologic settings were undertaken to improve understanding of source and transport controls on nitrate movement to public supply wells (PSW) in principal aquifers of the United States. Chemical, isotopic, and age tracer data show that agricultural fertilizers and urban septic leachate were the primary sources of large nitrate concentrations in PSW capture zones at Modesto, California (Central Valley aquifer system) and York, Nebraska (High Plains aquifer). Urban septic leachate and fertilizer (possibly nonfarm) were the primary sources of large nitrate concentrations in PSW capture zones at Woodbury, Connecticut (glacial aquifer system), and Tampa, Florida (Floridan aquifer system), respectively. Nitrate fluxes to the water table were larger in agricultural settings than urban settings, indicating that it would be beneficial to reduce PSW capture zone areas in agricultural regions. Mixing calculations indicate that about 50 to 85% of the nitrate in water from the PSW could be from those modern anthropogenic sources, with the remainder coming from sources in old (>50 years) recharge or sources in young recharge in undisturbed settings such as forests. Excess N 2 concentrations and age tracers showed that denitrification at Modesto occurred gradually (first‐order rate constant of 0.02/a) in a thick reaction zone following a ∼30‐year lag time after recharge. Denitrification generally was not an important nitrate sink at Woodbury. At York and Tampa, denitrification occurred rapidly (0.5 to 6/a) in thin reaction zones in fine‐grained sediments that separated the anoxic PSW producing zones from overlying oxic, high‐nitrate ground water. Particle tracking showed that a major pathway by which anthropogenic nitrate reached the York and Tampa PSW was by movement through long well screens crossing multiple hydrogeologic units (York) and by movement through karst features (Tampa), processes which reduced ground water residence times in the denitrifying zones. These results illustrate how PSW vulnerability to nitrate contamination depends on complex variations and interactions between contaminant sources, reaction rates, transit times, mixing, and perturbation of ground water flow in contrasting hydrogeologic settings.

Water Resources Research↗
USGS2007

Peat porewater chloride concentration profiles in the Everglades during wet/dry cycles from January 1996 to June 1998: Field measurements and theoretical analysis

Water quality is a key aspect of the Everglades Restoration Project, the largest water reclamation and ecosystem management project proposed in the United States. Movement of nutrients and contaminants to and from Everglades peat porewater could have important consequences for Everglades water quality and ecosystem restoration activities. In a study of Everglades porewater, we observed complex, seasonally variable peat porewater chloride concentration profiles at several locations. Analyses and interpretation of these changing peat porewater chloride concentration profiles identifies processes controlling conservative solute movement at the peat–surface water interface, that is, solutes whose transport is minimally affected by chemical and biological reactions. We examine, with an advection–diffusion model, how alternating wet and dry climatic conditions in the Florida Everglades mediate movement of chloride between peat porewater and marsh surface water. Changing surface water–chloride concentrations alter gradients at the interface between peat and overlying water and hence alter chloride flux across that interface. Surface water chloride concentrations at two frequently monitored sites vary with marsh water depth, and a transfer function was developed to describe daily marsh surface water chloride concentration as a function of marsh water depth. Model results demonstrate that porewater chloride concentrations are driven by changing surface water chloride concentrations, and a sensitivity analysis suggests that inclusion of advective transport in the model improves the agreement between the calculated and the observed chloride concentration profiles.

Hydrological Processes↗
USGS2003

Groundwater flow, heat transport, and water table position within volcanic edifices: Implications for volcanic processes in the Cascade Range

The position of the water table within a volcanic edifice has significant implications for volcano hazards, geothermal energy, and epithermal mineralization. We have modified the HYDROTHERM numerical simulator to allow for a free-surface (water table) upper boundary condition and a wide range of recharge rates, heat input rates, and thermodynamic conditions representative of continental volcano-hydrothermal systems. An extensive set of simulations was performed on a hypothetical stratovolcano system with unconfined groundwater flow. Simulation results suggest that the permeability structure of the volcanic edifice and underlying material is the dominant control on water table elevation and the distribution of pressures, temperatures, and fluid phases at depth. When permeabilities are isotropic, water table elevation decreases with increasing heat flux and increases with increasing recharge, but when permeabilities are anisotropic, these effects can be much less pronounced. Several conditions facilitate the ascent of a hydrothermal plume into a volcanic edifice: a sufficient source of heat and magmatic volatiles at depth, strong buoyancy forces, and a relatively weak topography-driven flow system. Further, the plume must be connected to a deep heat source through a pathway with a time-averaged effective permeability ???1 ?? 10-16 m2, which may be maintained by frequent seismicity. Topography-driven flow may be retarded by low permeability in the edifice and/or the lack of precipitation recharge; in the latter case, the water table may be relatively deep. Simulation results were compared with observations from the Quaternary stratovolcanoes along the Cascade Range of the western United States to infer hydrothermal processes within the edifices. Extensive ice caps on many Cascade Range stratovolcanoes may restrict recharge on the summits and uppermost flanks. Both the simulation results and limited observational data allow for the possibility that the water table beneath the stratovolcanoes is relatively deep.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

MODFLOW-2000, the U.S. Geological Survey modular ground-water model -- Documentation of the Model-Layer Variable-Direction Horizontal Anisotropy (LVDA) capability of the Hydrogeologic-Unit Flow (HUF) package

This report documents the model-layer variable-direction horizontal anisotropy (LVDA) capability of the Hydrogeologic-Unit Flow (HUF) Package of MODFLOW-2000. The LVDA capability allows the principal directions of horizontal anisotropy to be different than the model-grid row and column directions, and for the directions to vary on a cell-by-cell basis within model layers. The HUF Package calculates effective hydraulic properties for model grid cells based on hydraulic properties of hydrogeologic units with thicknesses defined independently of the model layers. These hydraulic properties include, among other characteristics, hydraulic conductivity and a horizontal anisotropy ratio. Using the LVDA capability, horizontal anisotropy direction is defined for model grid cells within which one or more hydrogeologic units may occur. For each grid cell, the HUF Package calculates the effective horizontal hydraulic conductivity along the primary direction of anisotropy using the hydrogeologic-unit hydraulic conductivities, and calculates the effective horizontal hydraulic conductivity along the orthogonal anisotropy direction using the effective primary direction hydraulic conductivities and horizontal anisotropy ratios. The direction assigned to the model layer effective primary hydraulic conductivity is specified using a new data set defined by the LVDA capability, when active, to calculate coefficients needed to solve the ground-water flow equation. Use of the LVDA capability is illustrated in four simulation examples, which also serve to verify hydraulic heads, advective-travel paths, and sensitivities calculated using the LVDA capability. This version of the LVDA capability defines variable-direction horizontal anisotropy using model layers, not the hydrogeologic units defined by the HUF Package. This difference needs to be taken into account when designing model layers and hydrogeologic units to produce simulations that accurately represent a given field problem. This might be a reason, for example, to make model layer boundaries coincide with hydrogeologic-unit boundaries in all or part of a model grid.

Open-File Report↗
USGS1998

An implicit dispersive transport algorithm for the US Geological Survey MOC3D solute-transport model

This report documents an extension to the U.S. Geological Survey MOC3D transport model that incorporates an implicit-in-time difference approximation for the dispersive transport equation, including source/sink terms. The original MOC3D transport model (Version 1) uses the method of characteristics to solve the transport equation on the basis of the velocity field. The original MOC3D solution algorithm incorporates particle tracking to represent advective processes and an explicit finite-difference formulation to calculate dispersive fluxes. The new implicit procedure eliminates several stability criteria required for the previous explicit formulation. This allows much larger transport time increments to be used in dispersion-dominated problems. The decoupling of advective and dispersive transport in MOC3D, however, is unchanged. With the implicit extension, the MOC3D model is upgraded to Version 2. A description of the numerical method of the implicit dispersion calculation, the data-input requirements and output options, and the results of simulator testing and evaluation are presented. Version 2 of MOC3D was evaluated for the same set of problems used for verification of Version 1. These test results indicate that the implicit calculation of Version 2 matches the accuracy of Version 1, yet is more efficient than the explicit calculation for transport problems that are characterized by a grid Peclet number less than about 1.0.

Water-Resources Investigations Report↗
USGS1992

D4Z: A new renumbering for iterative solution of ground-water flow and solute-transport equations

D4 zig-zag (D4Z) is a new renumbering scheme for producing a reduced matrix to be solved by an incomplete LU preconditioned, restarted conjugate-gradient iterative solver. By renumbering alternate diagonals in a zig-zag fashion, a very low sensitivity of convergence rate to renumbering direction is obtained. For two demonstration problems involving groundwater flow and solute transport, iteration counts are related to condition numbers and spectra of the reduced matrices.

Conference Paper↗
USGS1987

HST3D; a computer code for simulation of heat and solute transport in three-dimensional ground-water flow systems

The Heat- and Soil-Transport Program (HST3D) simulates groundwater flow and associated heat and solute transport in three dimensions. The three governing equations are coupled through the interstitial pore velocity, the dependence of the fluid density on pressure, temperature, the solute-mass fraction , and the dependence of the fluid viscosity on temperature and solute-mass fraction. The solute transport equation is for only a single, solute species with possible linear equilibrium sorption and linear decay. Finite difference techniques are used to discretize the governing equations using a point-distributed grid. The flow-, heat- and solute-transport equations are solved , in turn, after a particle Gauss-reduction scheme is used to modify them. The modified equations are more tightly coupled and have better stability for the numerical solutions. The basic source-sink term represents wells. A complex well flow model may be used to simulate specified flow rate and pressure conditions at the land surface or within the aquifer, with or without pressure and flow rate constraints. Boundary condition types offered include specified value, specified flux, leakage, heat conduction, and approximate free surface, and two types of aquifer influence functions. All boundary conditions can be functions of time. Two techniques are available for solution of the finite difference matrix equations. One technique is a direct-elimination solver, using equations reordered by alternating diagonal planes. The other technique is an iterative solver, using two-line successive over-relaxation. A restart option is available for storing intermediate results and restarting the simulation at an intermediate time with modified boundary conditions. This feature also can be used as protection against computer system failure. Data input and output may be in metric (SI) units or inch-pound units. Output may include tables of dependent variables and parameters, zoned-contour maps, and plots of the dependent variables versus time. (Lantz-PTT)

Water-Resources Investigations Report↗
USGS1983

A two-dimensional graphing program for the Tektronix 4050-series graphics computers

A refined, two-dimensional graph-plotting program was developed for use on Tektronix 4050-series graphics computers. Important features of this program include: any combination of logarithmic and linear axes, optional automatic scaling and numbering of the axes, multiple-curve plots, character or drawn symbol-point plotting, optional cartridge-tape data input and plot-format storage, optional spline fitting for smooth curves, and built-in data-editing options. The program is run while the Tektronix is not connected to any large auxiliary computer, although data from files on an auxiliary computer easily can be transferred to data-cartridge for later plotting. The user is led through the plot-construction process by a series of questions and requests for data input. Five example plots are presented to illustrate program capability and the sequence of program operation. (USGS)

Water-Resources Investigations Report↗