USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Hydrology”

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 307 records · Page 17Linked to original sources
USGS1989

Seepage through a hazardous-waste trench cover

Water movement through a waste-trench cover under natural conditions at a low-level radioactive waste disposal site in northwestern Illinois was studied from July 1982 to June 1984, using tensiometers, a moisture probe, and meteorological instruments. Four methods were used to estimate seepage: the Darcy, zero-flux plane, surface-based water-budget, and groundwater-based water-budget methods. Annual seepage estimates ranged from 48 to 216 mm (5-23% of total precipitation), with most seepage occurring in spring. The Darcy method, although limited in accuracy by uncertainty in hydraulic conductivity, was capable of discretizing seepage in space and time and indicated that seepage varied by almost an order of magnitude across the width of the trench. Lowest seepage rates occurred near the center of the cover, where seepage was gradual. Highest rates occurred along the edge of the cover, where seepage was highly episodic, with 84% of the total there being traced to wetting fronts from 28 individual storms. Limitations of the zero-flux-plane method were severe enough for the method to be judged inappropriate for use in this study.Water movement through a waste-trench cover under natural conditions at a low-level radioactive waste disposal site in northwestern Illinois was studied from July 1982 to June 1984, using tensiometers, a moisture probe, and meteorological instruments. Four methods were used to estimate seepage: the Darcy, zero-flux plane, surface-based water-budget, and groundwater-based water-budget methods. Annual seepage estimates ranged from 48 to 216mm (5-23% of total precipitation), with most seepage occurring in spring. The Darcy method, although limited in accuracy by uncertainty in hydraulic conductivity, was capable of discretizing seepage in space and time and indicated that seepage varied by almost an order of magnitude across the width of the trench. Lowest seepage rates occurred near the center of the cover, where seepage was gradual. Highest rates occurred along the edge of the cover, where seepage was highly episodic, with 84% of the total there being traced to wetting fronts from 28 individual storms. Limitations of the zero-flux-plane method were severe enough for the method to be judged inappropriate for use in this study.

Illinois↗
USGS1990

Dolomite dissolution rates and possible Holocene dedolomitization of water-bearing units in the Edwards aquifer, south-central Texas

Rates of dolomite dissolution can be used to test the concept, based on geomorphologic evidence, that a major part of the Edwards aquifer could have formed within the Holocene, a timeframe of approximately 10,000 years. During formation of the aquifer in the Edwards limestone (Cretaceous, Albian) of the Balcones fault zone, dolomite dissolution and porosity development were synchronous and the result of mixing-zone dedolomitization. Initiation of the mixing zone in the early Holocene (∼11,000 years before present) is suggested by the maximum age of formation of major discharge sites that allowed the influx of meteoric water into brine-filled, dolomitic preaquifer units. Dedolomitization, the dissolution of dolomite and net precipitation of calcite, has left aquifer units that are calcitic, and 40 vol.% interconnected pore space. The mass of dolomite missing is obtained by comparison of stratigraphically equivalent altered and unaltered units. One dissolution rate (1.76 × 10 −4 mmol dolomite kg H2O −1 yr −1 ) is determined from this mass, 10 4 yr reaction time, and a log-linear function describing the increase in mass discharge (three orders of magnitude) during aquifer formation. The second estimated dissolution rate is obtained from the mass transfer of dolomite to solution calculated from the increase in magnesium in pore fluids selected from the modern aquifer to represent a typical flowpath during aquifer formation. A reaction time of 10 4 yr for this mass transfer yields a rate of 0.56 × 10 −4 mmol dolomite kg H2O −1 yr −1 Both of these rates are comparable to modern rates of dolomite dissolution (0.3 to 4.5 × 10 −4 mmol dolomite kg H2O −1 yr −1 ) calculated from measured reaction times in the Tertiary Floridan aquifer system in Florida and the Madison aquifer in the Mississippian Madison Limestone of the Northern Great Plains. Similarity of these rates to the estimated paleo-rates of dolomite dissolution supports a 10 4 yr reaction timeframe. Both of these rates are comparable to modern rates of dolomite dissolution (0.3 to 4.5 × 10 −4 mmol dolomite kg H2O −1 yr −1 ) calculated from measured reaction times in the Tertiary Floridan aquifer system in Florida and the Madison aquifer in the Mississippian Madison Limestone of the Northern Great Plains. Similarity of these rates to the estimated paleo-rates of dolomite dissolution supports a 10 4 yr reaction timeframe. The Holocene reaction time also can be compared to a series of reaction times calculated by assuming that the mass of dolomite missing from the Edwards was removed at rates observed in the Floridan and Madison aquifers. These reaction times (for complete removal of dolomite) range from 2700 to 58,500 yr and span the Pleistocene-Holocene boundary. Finally, an estimated dolomite reaction rate during dedolomitization of the Edwards aquifer based on surface area of exposed dolomite [mmol cm −2 s −1 (millimoles per square centimeter per second)] may be approximated from reaction times. This rate is directly a function of the mass of dolomite removed and the surface area exposed per pore volume passing through the rock. The surface area is available from the observed dolomite rhomb size in unaltered rock. The rate of pore fluid movement is obtained from the averaged annual discharge. Rates during formation of the Edwards aquifer calculated from all reaction times range from 10 −13 to 10 −14 mmol dolomite cm −2 s −1 . These rates are faster than rates (10 −18 mmol cm −2 s −1 ), measured in the pure laboratory system, CaMg(CO 3 ) 2 CO 2 H 2 O, but slower than rates determined in an alpine stream study (10 −10 to 10 −11 mmol cm −2 s −1 ) where cold glacial melt water flows over dolostone. Dolomite dissolution rates from both the Edwards and other aquifers support the concept that a major part of the Edwards aquifer could have formed within the Holocene.

Texas↗
USGS1990

Monitoring moisture storage in trees using time domain reflectometry

Laboratory and field tests were performed to examine the feasibility of using time domain reflectometry (TDR) to monitor changes in the moisture storage of the woody parts of trees. To serve as wave guides for the TDR signal, pairs of stainless steel rods (13 cm long, 0.32 cm in diameter, and 2.5 cm separation) were driven into parallel pilot holes drilled into the woody parts of trees, and a cable testing oscilloscope was used to determine the apparent dielectric constant. A laboratory calibration test was performed on two sapwood samples, so that the relation between the volumetric water content and the apparent dielectric constant of the sapwood could be determined over a range of water contents. The resulting calibration curve for these sapwood samples was significantly different than the general calibration curve used for soils, showing a smaller change in the apparent dielectric constant for a given change in the volumetric water content than is typical for soils. The calibration curve was used to estimate the average volumetric water content to a depth of 13 cm in living trees. One field experiment was conducted on an English walnut tree ( Juglans regia ) with a diameter of 40 cm, growing in a flood-irrigated orchard on a Hanford sandy loam near Modesto, California (U.S.A.). Rods were driven into the tree at about 50 cm above the soil surface and monitored hourly for the month of August, 1988. The moisture content determined by TDR showed a gradual decrease from 0.44 to 0.42 cm 3 cm −3 over a two week period prior to flood irrigation, followed by a rapid rise to 0.47 cm 3 cm −3 over a four day period after irrigation, then again a gradual decline approaching the next irrigation. A second field experiment was made on ten evergreen and deciduous trees with diameters ranging from 30 to 120 cm, growing in the foothills of the Coast Range of central California. Rods were driven into each tree at 50 to 100 cm above the soil surface and monitored on a biweekly to monthly basis for over a year. Most trees showed an early spring maximum in moisture content determined by TDR associated with leaf growth, and a late summer minimum in moisture content associated with the end of the dry season. Moisture contents ranged from 0.20 to 0.70 cm 3 cm −3 , with an annual percentage change in moisture of 15% to 70% depending on species and environmental conditions. A final field test was performed in northern New Mexico (U.S.A.) to examine the effect of trunk freezing on TDR measurements. This test confirmed that freezing conditions were recorded as a total loss of liquid water by the TDR method. These results suggest that further TDR calibration for wood, plus some understanding of the relation between tree moisture and physiological stress could be useful to several disciplines, ranging from irrigation scheduling to watershed management to forest ecology.

California↗
USGS1990

Simulation of dispersion in layered coastal aquifer systems

A density-dependent solute-transport formulation is used to examine ground-water flow in layered coastal aquifers. The numerical experiments indicate that although the transition zone may be thought of as an impermeable 'sharp' interface with freshwater flow parallel to the transition zone in homogeneous aquifers, this is not the case for layered systems. Freshwater can discharge through the transition zone in the confining units. Further, for the best simulation of layered coastal aquifer systems, either a flow-direction-dependent dispersion formulation is required, or the dispersivities must change spatially to reflect the tight thin confining unit.

Journal of Hydrology↗
USGS1990

Storm-runoff generation in the Permanente Creek drainage basin, west central California - An example of flood-wave effects on runoff composition

Variations in the isotopic and chemical composition of storm runoff in the 10.6-km 2 Permanente Creek basin, Santa Clara County, California, indicate that changes in water composition lag behind changes in streamflow. This lag occurs even though field observations and rainfall-runoff modeling indicate that much of the storm runoff must be composed of "new" water running off impervious surfaces. The apparent incompatibility posed by the presence of "old" water and the direct and indirect evidence that surface runoff contributes substantially to storm runoff can be explained if initial rises in streamflow result from effects of flood waves. Flood waves composed of old channel water reach downstream locations ahead of the new water derived from impervious areas. By this mechanism, streamflow can rise rapidly in response to surface runoff and still be composed of large amounts of old water. Data collected in Permanente Creek indicate that flood waves can occur even in small basins, at least when those basins contain impervious surfaces.

California↗
USGS1990

Effect of faults on fluid flow and chloride contamination in a carbonate aquifer system

A unified, multidiscipline hypothesis is proposed to explain the anomalous pattern by which chloride has been found in water of the Upper Floridan aquifer in Brunswick, Glynn County, Georgia. Analyses of geophysical, hydraulic, water chemistry, and aquifer test data using the equivalent porous medium (EPM) approach are used to support the hypothesis and to improve further the understanding of the fracture-flow system in this area. Using the data presented herein we show that: (1) four major northeast-southwest trending faults, capable of affecting the flow system of the Upper Floridan aquifer, can be inferred from structural analysis of geophysical data and from regional fault patterns; (2) the proposed faults account for the anomalous northeastward elongation of the potentiometric surface of the Upper Floridan aquifer; (3) the faults breach the nearly impermeable units that confine the Upper Floridan aquifer from below, allowing substantial quantities of water to leak vertically upward; as a result, aquifer transmissivity need not be excessively large (as previously reported) to sustain the heavy, long-term pumpage at Brunswick without developing a steep cone of depression in the potentiometric surface; (4) increased fracturing at the intersection of the faults enhances the development of conduits that allow the upward migration of high-chloride water in response to pumping from the Upper Floridan aquifer; and (5) the anomalous movement of the chloride plume is almost entirely controlled by the faults.

Georgia↗
USGS1990

Movement and fate of atrazine and bromide in central Kansas croplands

Two flooding experiments were conducted at two sites with different soils to study the transport and fate of the commonly used herbicide atrazine and inorganic chemicals in the Great Bend Prairie croplands of south-central Kansas. The instantaneous profile method supplemented by the use of an organic (atrazine) and an inorganic (bromide) tracer chemical was used to characterize in situ the hydraulic and chemical properties of the appropriately instrumented field sites. Atrazine readily degraded to hydroxyatrazine and biodegradation by-products and was not detected deeper in the soil profile and underlying shallow aquifer. The classical processes of chemical movement based on porous media-equilibrium-diffuse flow did not fit the data well at either site. Incompletely mixed, slug flow appeared to predominate at one of the sites and preferential flow at the other. The slug movement caused 'piston-type' displacement of more saline solutions in the soil profile to the shallow water table. Recommendations for conducting related field studies based on our sampling experience are given.

Kansas↗
USGS1990

Theory and application of an approximate model of saltwater upconing in aquifers

Motion and mixing of salt water and fresh water are vitally important for water-resource development throughout the world. An approximate model of saltwater upconing in aquifers is developed, which results in three non-linear coupled equations for the freshwater zone, the saltwater zone, and the transition zone. The description of the transition zone uses the concept of a boundary layer. This model invokes some assumptions to give a reasonably tractable model, considerably better than the sharp interface approximation but considerably simpler than a fully three-dimensional model with variable density. We assume the validity of the Dupuit-Forchheimer approximation of horizontal flow in each layer. Vertical hydrodynamic dispersion into the base of the transition zone is assumed and concentration of the saltwater zone is assumed constant. Solute in the transition zone is assumed to be moved by advection only. Velocity and concentration are allowed to vary vertically in the transition zone by using shape functions. Several numerical techniques can be used to solve the model equations, and simple analytical solutions can be useful in validating the numerical solution procedures. We find that the model equations can be solved with adequate accuracy using the procedures presented. The approximate model is applied to the Smoky Hill River valley in central Kansas. This model can reproduce earlier sharp interface results as well as evaluate the importance of hydrodynamic dispersion for feeding salt water to the river. We use a wide range of dispersivity values and find that unstable upconing always occurs. Therefore, in this case, hydrodynamic dispersion is not the only mechanism feeding salt water to the river. Calculations imply that unstable upconing and hydrodynamic dispersion could be equally important in transporting salt water. For example, if groundwater flux to the Smoky Hill River were only about 40% of its expected value, stable upconing could exist where hydrodynamic dispersion into a transition zone is the primary mechanism for moving salt water to the river. The current model could be useful in situations involving dense saltwater layers.

Journal of Hydrology↗
USGS1990

A method to extract soil water for stable isotope analysis

A method has been developed to extract soil water for determination of deuterium (D) and 18 O content. The principle of this method is based on the observation that water and toluene form an azeotropic mixture at 84.1°C, but are completely immiscible at ambient temperature. In a specially designed distillation apparatus, the soil water is distilled at 84.1°C with toluene and is separated quantitatively in the collecting funnel at ambient temperature. Traces of toluene are removed and the sample can be analyzed by mass spectrometry. Kerosene may be substituted for toluene. The accuracy of this technique is ± 2 and ± 0.2‰ "> ± 0.2‰ , respectively, for δD and δ 18 O. Reduced accuracy is obtained at low water contents.

Journal of Hydrology↗
USGS1991

Estimation of suspended-sediment rating curves and mean suspended-sediment loads

Suspended-sediment loads are often estimated from an empirical relation between suspended-sediment load ( L ) and streamflow ( S ). This relation is usually defined as a power function, L = aS h , and is referred to as a suspended-sediment rating curve. This function can be formulated as either a linear or non-linear model to find the solution of the rating-curve parameters ( a and b ). Formulation of the power function as a linear model requires a logarithmic transformation to linearize the function and a subsequent correction for transformation bias. Rating-curve parameter estimates for both the bias-corrected, transformed-linear or non-linear models can be obtained by the method of least squares. Each model has distinct advantages and disadvantages. A unique solution of the parameters of the transformed-linear model may be obtained algebraically. These parameter estimates have some optimal properties when certain attainable conditions are met. However, the parameter estimates must be corrected for transformation bias when obtained this way. Parameter estimates obtained for the non-linear model do not require a correction for transformation bias. However, these estimates must be obtained by iterative methods which do not always converge to a solution. In addition, the residual errors of the non-linear model typically are not identically distributed throughout the range of streamflow values. This problem adversely affects the precision of the parameter estimates. Weighted non-linear least squares can be used to improve the parameter estimates for the non-linear model, but the weights must be approximated and their appropriate form may be difficult to determine. A simulation study was done to evaluate: (1) the accuracy and precision of parameter estimates for the bias-corrected, transformed-linear and non-linear models obtained by the method of least squares; (2) the accuracy of mean suspended-sediment loads calculated by the flow-duration, rating-curve method using model parameters obtained by the alternative methods. Parameter estimates obtained by least squares for the bias-corrected, transformed-linear model were considerably more precise than those obtained for the non-linear or weighted non-linear model. The accuracy of parameter estimates obtained for the biascorrected, transformed-linear and weighted non-linear model was similar and was much greater than the accuracy obtained by non-linear least squares. The improved parameter estimates obtained by the biascorrected, transformed-linear or weighted non-linear model yield estimates of mean suspended-sediment load calculated by the flow-duration, rating-curve method that are more accurate and precise than those obtained for the non-linear model.

Journal of Hydrology↗
USGS1991

Fate of acetone in an outdoor model stream with a nitrate supplement, southern Mississippi, U.S.A.

The fate of acetone in an outdoor model stream to which nitrate was added as a nutrient supplement was determined. The stream, in southern Mississippi, U.S.A. was 234 m long. Water was supplied to the stream by an artesian well at about 1.21 s −1 , resulting in a mean water velocity of about 0.5 m min −1 . Acetone was injected continuously for 26 days resulting in concentrations of 20–40 mg l −1 . A nitrate solution was injected for 21 days resulting in an instream concentration of about 1.7 mg l −1 at the upstream end of the stream. Rhodamine-WT dye was used to determine the travel time and dispersion characteristics of the stream, and t-butyl alcohol was used to determine the volatilization characteristics. Volatilization controlled the fate of acetone in the model stream. The lack of substantial bacterial degradation of acetone was contrary to expectations based on the results of laboratory degradation studies using model stream water enriched with nitrate. A possible explanation for the lack of significant degradation in the model stream may be the limited 6-h residence time of the acetone in the stream.

Journal of Hydrology↗
USGS1991

Geochemistry of dissolved inorganic carbon in a Coastal Plain aquifer. 1. Sulfate from confining beds as an oxidant in microbial CO2 production

A primary source of dissolved inorganic carbon (DIC) in the Black Creek aquifer of South Carolina is carbon dioxide produced by microbially mediated oxidation of sedimentary organic matter. Groundwater chemistry data indicate, however, that the available mass of inorganic electron acceptors (oxygen, Fe(III), and sulfate) and observed methane production is inadequate to account for observed CO 2 production. Although sulfate concentrations are low (approximately 0.05–0.10 mM) in aquifer water throughout the flow system, sulfate concentrations are greater in confining-bed pore water (0.4–20 mM). The distribution of culturable sulfate-reducing bacteria in these sediments suggests that this concentration gradient is maintained by greater sulfate-reducing activity in sands than in clays. Calculations based on Fick's Law indicate that possible rates of sulfate diffusion to aquifer sediments are sufficient to explain observed rates of CO 2 production (about 10 −5 mmoll −1 year −1 ), thus eliminating the apparent electron-acceptor deficit. Furthermore, concentrations of dissolved hydrogen in aquifer water are in the range characteristic of sulfate reduction (2–6 nM), which provides independent evidence that sulfate reduction is the predominant terminal electron-accepting process in this system. The observed accumulation of pyrite- and calcite-cemented sandstones at sand-clay interfaces is direct physical evidence that these processes have been continuing over the history of these sediments.

Journal of Hydrology↗
USGS1991

Geohydrologic, geochemical, and geologic controls on the occurrence of radon in ground water near Conifer, Colorado, USA

Integrated studies of geohydrology, geochemistry, and geology of crystalline rocks in the vicinity of Conifer, Colorado, reveal that radon concentrations do not correlate with variations in concentrations of other dissolved species. Concentrations of major ions show systematic variations along selected groundwater flowpaths, whereas radon concentrations are dependent on local geochemical and geologic phenomena (such as localized uranium concentration in the rock or the presence of faults or folds). When radon enters the flow system, concentrations do not increase along flowpaths because its decay rate is fast relative to groundwater flow rates. Radon-222 is not in secular equilibrium with 238 U and 226 Ra in the water. Therefore, most of the 238 U and 226 Ra necessary to support the waterborne 222 Rn must be present locally in the rock. High concentrations of dissolved radon are not found in zones of high transmissivity, and transmissivity is not correlated with rock type in the study area. A higher transmissivity can be indicative of higher water-volume to rock-surface-area ratios, which could effectively dilute 222 Rn entering the water and/or may indicate that emanated radon is carried away more rapidly. Water samples collected from individual wells over periods of several months showed significant fluctuations in the dissolved 222 Rn content. This fluctuation may be controlled by changes in the contributions of water-producing zones within the well resulting from seasonal fluctuations of the water table and/or pumping stresses.

Colorado↗
USGS1991

Pumping tests in non-uniform aquifers - The linear strip case

Many pumping tests are performed in geologic settings that can be conceptualized as a linear infinite strip of one material embedded in a matrix of differing flow properties. A semi-analytical solution is presented to aid the analysis of drawdown data obtained from pumping tests performed in settings that can be represented by such a conceptual model. Integral transform techniques are employed to obtain a solution in transform space that can be numerically inverted to real space. Examination of the numerically transformed solution reveals several interesting features of flow in this configuration. If the transmissivity of the strip is much higher than that of the matrix, linear and bilinear flow are the primary flow regimes during a pumping test. If the contrast between matrix and strip properties is not as extreme, then radial flow should be the primary flow mechanism. Sensitivity analysis is employed to develop insight into the controls on drawdown in this conceptual model and to demonstrate the importance of temporal and spatial placement of observations. Changes in drawdown are sensitive to the transmissivity of the strip for a limited time duration. After that time, only the total drawdown remains a function of strip transmissivity. In the case of storativity, both the total drawdown and changes in drawdown are sensitive to the storativity of the strip for a time of quite limited duration. After that time, essentially no information can be gained about the storage properties of the strip from drawdown data. An example analysis is performed using data previously presented in the literature to demonstrate the viability of the semi-analytical solution and to illustrate a general procedure for analysis of drawdown data in complex geologic settings. This example reinforces the importance of observation well placement and the time of data collection in constraining parameter correlation, a major source of the uncertainty that arises in the parameter estimation procedure.

Journal of Hydrology↗
USGS1991

The effects of overburden stress on the specific storage and hydraulic conductivity of artesian aquifers

A mathematical algorithm is developed to determine the depth-dependent profiles of specific storage and hydraulic conductivity resulting from overburden stress in horizontally isotropic artesian aquifers. Vertical variations in the void ratio of the aquifer matrix brought about by overburden stresses determine the pore-volume compressibility and matrix permeability at specific depths within the aquifers which, in turn, determine the depth-dependent profiles of specific storage and hydraulic conductivity. Time-drawdown curves are obtained for two sands subjected to different overburden stresses. For shallow artesian aquifers with low overburden stress and high aquifer matrix compressibility, noticeable vertical gradients in specific storage occur. These vertical gradients cause deviations from the classical time-drawdown curves defined by the Theis solution. These deviations are negligible for deep artesian aquifers.

Journal of Hydrology↗
USGS1991

Effects of wetlands creation on groundwater flow

Changes in groundwater flow were observed near four Experimental Wetland Areas (EWAs) constructed along a reach of the Des Plaines River in northeastern Illinois. These changes were observed during monthly monitoring of groundwater elevation in nested piezometers and shallow observation wells before and after the wetlands were filled with water. A numerical model was calibrated with observed data and used to estimate seepage from the wetlands into the Des Plaines River. After the wetlands became operational, groundwater levels in adjacent wells increased by about 0.5m, while water levels in wells distant from the wetlands decreased. The increase in groundwater levels near the wetlands is a result of seepage from the wetlands. Numerical predictions of seepage from the wetlands are 60–150 m 3 day −1 for two wetlands situated over sand and gravel and less than 1 m 3 day −1 for two wetlands situated over clayey till. The difference in seepage rates is attributed to two factors. First, the hydraulic conductivity of the sand and gravel unit is greater than that of the till, and thus there is less mounding and a greater capacity for transmitting water beneath the wetlands overlying this deposit. Secondly, the wetlands located over till are groundwater flow-through ponds, whereas the wetlands over the sand and gravel are primarily groundwater recharge areas. The model was used to estimate that seepage from the wetlands will double groundwater discharge into the Des Plaines River and a tributary relative to pre-operational discharge from the study area. Overall, the wetlands have acted as a constant head boundary, stabilizing groundwater flow patterns.

Illinois↗
USGS1991

A comparison of short-term measurements of lake evaporation using eddy correlation and energy budget methods

Concurrent short-term measurements of evaporation from a shallow lake, using eddy correlation and energy budget methods, indicate that sensible and latent heat flux between lake and atmosphere, and energy storage in the lake, may vary considerably across the lake. Measuring net radiation with a net radiometer on the lake appeared to be more accurate than measuring incoming radiation nearby and modeling outgoing radiation. Short-term agreement between the two evaporation measurements was obtained by using an energy storage term that was weighted to account for the area-of-influence of the eddy correlation sensors. Relatively short bursts of evaporation were indicated by the eddy correlation sensors shortly after midnight on two of three occasions.

Oklahoma↗
USGS1991

The temperature dependence of ponded infiltration under isothermal conditions

A simple temperature-sensitive modification to the Green and Ampt infiltration equation is described; this assumes that the temperature dependence of the hydraulic conductivity is reciprocally equal to the temperature dependence of the viscosity of liquid water, and that both the transmission zone saturation and the wetting front matric potential gradient are independent of temperature. This modified Green and Ampt equation is compared with ponded, isothermal infiltration experiments run on repacked columns of Olympic Sand and Aiken Loam at 5, 25, and 60°C. Experimental results showed increases in infiltration rates of at least 300% between 5 and 60°C for both soil materials, with subsequent increases in cumulative infiltration of even greater magnitudes for the loam. There is good agreement between measured and predicted initial infiltration rates at 25°C for both soil materials, yet at 60°C, the predicted results overestimate initial infiltration rates for the sand and underestimate initial rates for the loam. Measurements of the wetting depth vs. cumulative infiltration indicate that the transmission zone saturation increased with increasing temperature for both soil materials. In spite of this increased saturation with temperature, the final infiltration rates at both 25 and 60°C were predicted accurately using the modified Green and Ampt equation. This suggests that increased saturation occurred primarily in dead-end pore spaces, so that transmission zone hydraulic conductivities were unaffected by these temperature-induced changes in saturation. In conclusion, except for initial infiltration rates at 60°C, the measured influence of temperature on infiltration rates was fully accounted for by the temperature dependence of the viscosity of liquid water.

Journal of Hydrology↗