Geology topics

W. E. Sanford

Publications and source records attributed to W. E. Sanford.

27 records · Page 2Linked to original sources

USGS2004

Numerical analysis of seawater circulation in carbonate platforms: II. The dynamic interaction between geothermal and brine reflux circulation

Density-driven seawater circulation may occur in carbonate platforms due to geothermal heating and / or reflux of water of elevated salinity. In geothermal circulation lateral contrasts in temperature between seawater and platform groundwaters warmed by the geothermal heat flux result in upward convective flow, with colder seawater drawn into the platform at depth. With reflux circulation, platform-top waters concentrated by evaporation flow downward, displacing less dense underlying groundwaters. We have used a variable density groundwater flow model to examine the pattern, magnitude and interaction of these two different circulation mechanisms, for mesosaline platform-top waters ( 50 ‰ ) and brines concentrated up to saturation with respect to gypsum ( 150 ‰ ) and halite ( 246 ‰ ). Geothermal circulation, most active around the platform margin, becomes restricted and eventually shut-off by reflux of brines from the platform interior towards the margin. The persistence of geothermal circulation is dependent on the rate of brine reflux, which is proportional to the concentration of platform-top brines and also critically dependent on the magnitude and distribution of permeability. Low permeability evaporites can severely restrict reflux whereas high permeability units in hydraulic continuity enhance brine transport. Reduction in permeability with depth and anisotropy of permeability (kv < < kh) focuses flow laterally in the shallow subsurface (<1 km), resulting in a horizontally elongated brine plume. Aquifer porosity and dispersivity are relatively minor controls on reflux. Platform brines can entrain surficial seawater when brine generating conditions cease but the platform-top remains submerged, a variant of reflux we term "latent reflux". Brines concentrated up to gypsum saturation have relatively long residence times of at least 100 times the duration of the reflux event. They thus represent a long-term control on post-reflux groundwater circulation, and consequently on the rates and spatial patterns of shallow burial diagenesis, such as dolomitization.

American Journal of Science↗

USGS2002

Source of solutes to the coastal sabkha of Abu Dhabi

An ascending-brine model is proposed to address the observed isotope geochemistry, solute composition, and solute and water fluxes in the coastal sabkha of the Emirate of Abu Dhabi. Mass-balance measurements document that >95% of the solutes are derived from ascending continental brines; minor amounts are derived from rainfall and from groundwater entering from upgradient areas. Nearly 100% of the annual water loss is from evaporation and not lateral discharge. Direct rainfall on the sabkha and subsequent recharge to the underlying aquifer account for ~ 90% of the annual water input to the system; the remaining 10% comes from both lateral and ascending groundwater flow. Thus, the water and solutes in this system are from different sources. Solute concentrations of conservative (i.e., nonreactive) elements in the coastal, sabkha-covered aquifer are consistent with the fluid pore volumes of ascending brine calculated from hydrologic properties. Calcium to sulfate ratios and sulfur isotopes are consistent with this source of solute from the underlying Tertiary formations. Recharging rainwater dissolves halite and other soluble minerals on the surface, causing the solution to become more dense and sink to the bottom of the aquifer where it vertically mixes with less dense ascending brines. Solutes are returned to the surface by capillary forces and recycled or lost from the system by eolian or fluvial processes. Thus, the system becomes vertically mixed, consistent with the presence of tritium throughout the aquifer; but there is essentially no horizontal mixing of seawater with groundwater. The observed seawater solutes in the supratidal zone come from interstitial seawater trapped by the rapid progradation of the sediments into the Arabian Gulf and are not refluxed or laterally mixed. The ascending-brine model contrasts significantly with both the seawater-flooding and evaporative-pumping models previously proposed as a source of solutes to the coastal sabkha of the Emirate of Abu Dhabi. Use of these earlier models leads to incorrect conclusions and raises serious questions about their applicability in the evaluation of sabkhat in the geologic record.

Abu Dhabi↗

USGS2002

Parameter estimation using carbon-14 ages: Lessons from the Danube-Tisza interfluvial region of Hungary

Parameter estimation was conducted on a groundwater model of the Danube-Tisza interfluvial region of Hungary. The model was calibrated using 300 water levels and 48 14C ages. The model provided a test of regression methods for a system with a large number of observations. Up to 103 parameters representing horizontal and vertical hydraulic conductivities and boundary conductances were assigned using point values and bilinear interpolation between points. The lowest errors were obtained using an iterative approach with groups of parameters, rather than estimating all of the parameters simultaneously. The model with 48 parameters yielded the lowest standard error of regression.

Acta Universitatis Carolinae, Geologica↗

USGS1997

Constant-concentration boundary condition: Lessons from the HYDROCOIN variable-density groundwater benchmark problem

In a solute-transport model, if a constant-concentration boundary condition is applied at a node in an active flow field, a solute flux can occur by both advective and dispersive processes. The potential for advective release is demonstrated by reexamining the Hydrologic Code Intercomparison (HYDROCOIN) project case 5 problem, which represents a salt dome overlain by a shallow groundwater system. The resulting flow field includes significant salinity and fluid density variations. Several independent teams simulated this problem using finite difference or finite element numerical models. We applied a method-of-characteristics model (MOCDENSE). The previous numerical implementations by HYDROCOIN teams of a constant-concentration boundary to represent salt release by lateral dispersion only (as stipulated in the original problem definition) was flawed because this boundary condition allows the release of salt into the flow field by both dispersion and advection. When the constant-concentration boundary is modified to allow salt release by dispersion only, significantly less salt is released into the flow field. The calculated brine distribution for case 5 depends very little on which numerical model is used, as long as the selected model is solving the proper equations. Instead, the accuracy of the solution depends strongly on the proper conceptualization of the problem, including the detailed design of the constant-concentration boundary condition. The importance and sensitivity to the manner of specification of this boundary does not appear to have been recognized previously in the analysis of this problem.

Water Resources Research↗

USGS1992

Large lake basins of the southern High Plains: Ground-water control of their origin?

The origin of the ∼40-50 topographically large lake basins on the southern High Plains of Texas and New Mexico has been an enigma. Previous workers have considered deflation or evaporite dissolution at depth and subsequent collapse as the most probable mechanisms. However, the eolian hypotheses have been unable to provide convincing arguments as to how the wind selectively erodes the thick, deflation-resistant calcrete "caprock" that is persistent over much of the southern High Plains. Furthermore, recent detailed studies on some of the basins show no significant evaporite dissolution at depth, and neither mechanism offers a satisfactory explanation as to why the basins are almost universally associated with subsurface topographic highs, or why they are absent where the High Plains aquifer thickness exceeds 60 m. We address these latter concerns and modify the deflation hypothesis by proposing that the calcrete caprock may never have been deposited in the areas now occupied by the basins. The absence of calcrete deposition is proposed to have resulted from high water tables caused by an increase in hydraulic gradient where aquifers thinned above bedrock highs. A high water table close to and/or intersecting the surface prevents deposition of calcrete, and, thus, the uncemented surface would be more susceptible to deflation than the surrounding calcrete-covered areas after decline of the water table. The rise in water table associated with bedrock highs is documented by numerical simulation using boundary conditions and hydrologic parameters representative of the southern High Plains.

Geology↗

USGS1991

Brine evolution and mineral deposition in hydrologically open evaporite basins

A lumped-parameter, solute mass-balance model is developed to define the role of water outflow from a well-mixed basin. A mass-balance model is analyzed with a geochemical model designed for waters with high ionic strengths. Two typical waters, seawater and a Na-HCO3 ground water, are analyzed to illustrate the control that the leakage ratio (or hydrologic openness of the basin) has on brine evolution and the suite and thicknesses of evaporite minerals deposited. The analysis suggests that brines evolve differently under different leakage conditions. -from Authors

American Journal of Science↗

USGS1990

Ground-water control of evaporite deposition

Topographically closed basins may be hydrologically open as a result of seepage losses to underlying or surrounding ground-water systems. In such cases, these losses can have a substantial control over the suite and the thicknesses of evaporite minerals formed in the basin. The ratio of ground-water outflow to inflow (flux ratio) in hydrologically open basins is as important in determining the mineralogy and thicknesses of evaporite deposits as the solute composition of the inflow water. Attainment of steady state flux ratios permits large thicknesses of two or three minerals to form rather than thin veneers of many minerals. The presence or absence of glauberite, mirabilite, halite, bloedite, polyhalite, and hexahydrite, caused by subtle changes in the ground-water seepage is illustrated using an example from the Southern High Plains of Texas and New Mexico. However, the model is general and is applicable with any solute composition including that of seawater and the use of surface rather than ground water. An analytical, lumped parameter, solute mass balance model is developed to define the concept of a ground-water flux ratio as it applies to topographically closed basins in which evaporation exceeds precipitation. Diffusion, advection, and density-driven flow are proposed as mechanisms by which solutes can escape to the ground water from these closed basins. The geochemical reaction computer program PHRQPITZ is used to document the effects of various flux ratios on the mineralogy and thickness of deposits. Solute analyses used in conjunction with the model can be used to screen prospective basins as well as to provide insights for exploratory drilling program.

Texas, New Mexico↗

USGS1989

Porosity development in coastal carbonate aquifers

Geochemical mixing theory suggests that the mixing of seawater and calcite-saturated fresh ground water can result in a solution that is undersaturated with respect to calcite. Previous studies of the mixing of such waters in carbonate rocks along certain coastlines have indicated that this mixing effect may be responsible for significant amounts of calcite dissolution and porosity development. In this study, potential rates of porosity development by calcite dissolution are assessed by combining geochemical mixing theory with the hydrodynamics of fresh-water-salt-water mixing zones in a coupled reaction- transport model. Results from the reaction-path model PHREEQE are used with a variable-density ground-water flow and solute-transport model to simulate an idealized cross section of a coastal carbonate aquifer. Results of the simulations indicate that the dissolution process is sensitive to fresh-water chemistry, ground-water velocities, and sea-level movement. Dissolution potential was evaluated at three field sites, and evidence from those sites is in general agreement with the simulation results. Dissolution rates indicated by the model show that under the proper conditions this dissolution mechanism can produce significant increases in porosity over relatively short spans of geologic time (tens of thousands of years).

Geology↗

USGS1985

A two-constituent solute-transport model for ground water having variable density

A numerical model has been developed to simulate solute transport and dispersion of either one or two constituents in groundwater where there is two-dimensional, density-dependent flow. The model is a modified version of the one documented by Konikow and Bredehoeft (1978), which uses finite-difference methods and the method of characteristics to solve the flow and transport equations. The model was tested on an idealized seawater intrusion problem for which an analytical solution has been developed. The results were nearly identical to those of other numerical models tested on the same problem. A description of the formats for the input data, a sample of input and output for a two-constituent example problem, and a listing of the Fortran program are presented. (Author 's abstract)

Water-Resources Investigations Report↗