John L. Nieber

Regionalized drought flow hydrographs from a mature glaciated plateau

The drought or base flow characteristics of six basins in the Finger Lakes region are obtained by considering for each available record the lower envelope of |dQ/dt| as a function of Q, where Q is the flow rate. This procedure avoids the uncertainty regarding a proper time reference after each rainfall event, and it eliminates the effects of evapotranspiration. The results suggest that among several expressions, Boussinesqs nonlinear solution of free surface groundwater flow is best suited to parameterize the observed hydrographs. The obtained parameters can be related to the basin characteristics, viz., drainage area and the total stream length, in accordance with relationships derived on the basis of the Dupuit-Boussinesq aquifer model. This result allows the determination of drought flow parameters for ungaged sites within the region.

Modeling and field evidence of finger formation and finger recurrence in a water repellent sandy soil

With prolonged rainfall, infiltrating wetting fronts in water repellent soils may become unstable, leading to the formation of high-velocity flow paths, the so-called fingers. Finger formation is generally regarded as a potential cause for the rapid transport of water and contaminants through the unsaturated zone of soils. For the first time, field evidence of the process of finger formation and finger recurrence is given for a water repellent sandy soil. Theoretical analysis and model simulations indicate that finger formation results from hysteresis in the water retention function, and the character of the formation depends on the shape of the main wetting and main drainage branches of that function. Once fingers are established, hysteresis causes fingers to recur along the same pathways during following rain events. Leaching of hydrophobic substances from these fingered pathways makes the soil within the pathways more wettable than the surrounding soil. Thus, in the long-term, instability-driven fingers might become heterogeneity-driven fingers.

Physics of water repellent soils

Although it is generally well known that water repellent soils have distinct preferential flow patterns, the physics of this phenomenon is not well understood. In this paper, we show that water repellency affects the soil water contact angle and this, in turn, has a distinct effect on the constitutive relationships during imbibing. Using these constitutive relationships, unstable flow theory developed for coarse grained soils can be used to predict the shape and water content distribution for water repellent soils. A practical result of this paper is that with a basic experimental setup, we can characterize the imbibing front behavior by measuring the water entry pressure and the imbibing soil characteristic curve from the same heat treated soil. q 2000 Elsevier Science B.V. All rights reserved.

MODELING FINGER DEVELOPMENT AND PERSISTENCE IN INITIALLY DRY POROUS MEDIA

Abstract The mechanism for the growth and persistence of gravity-driven fingered flow of water in initially dry porous media is described. A Galerkin finite element solution of the two-dimensional Richards equation with the associated parameter equations for capillary hysteresis in the water retention function is presented. A scheme for upstream weighting of internodal unsaturated hydraulic conductivities is applied to limit smearing of steep wetting fronts. The growth and persistence of a single finger in an initially dry porous media is simulated using this numerical solution scheme. To adequately simulate fingered flow, it was found that the upstream weighting factor had to be negative, meaning that the internodal unsaturated hydraulic conductivities were weighted more by the downstream node. It is shown that the growth and persistence of a finger is sensitive to the character of the porous media water retention functions. For porous media where the water-entry capillary pressure on the main wetting function is less than the air-entry capillary pressure on the main drainage function, a small perturbation will grow into a finger, and during sequential drainage and wetting the finger will persist. In contrast, for porous media where the water-entry capillary pressure on the main wetting function is greater than the air-entry capillary pressure on the main drainage function, the same small perturbation will dissipate by capillary diffusion. The finger widths derived from the numerical simulation are similar to those predicted by analytical theory.

Stable or unstable wetting fronts in water repellent soils - effect of antecedent soil moisture content

Dry water repellent soils are known to inhibit water infiltration, ultimately forcing water to flow via preferential paths through the vadose zone. To study water flow and transport in a water repellent sandy soil, a bromide tracer experiment had been carried out, which started in the fall after winter wheat had been sown. Despite the uniform tracer application, soil core sampling indicated that bromide concentrations varied largely from place to place. Wetter sites in the experimental field received more bromide, due to lateral transport through a thin top layer. Wetting fronts infiltrated deeper here, leading to perturbed wetting fronts in the experimental field. In contrast to what was expected, the wetting front perturbations did not grow to fingers. Numerical results indicate that this was attributed to the relatively high soil water contents during the experiment, which caused the soil to be wettable instead of water repellent. The water-entry capillary pressure of the secondary wetting branch exceeds the air-entry capillary pressure of the primary drainage branch in this case. In the opposite situation, with the water-entry capillary pressure of the secondary wetting branch beneath the air-entry capillary pressure of the primary drainage branch, perturbations would have grown to fingers. Such a situation occurs during infiltration in initially dry, water repellent soil. The results presented illustrate the effect of antecedent moisture conditions on the formation of stable and unstable wetting fronts, and its relation to the moment of tracer application.

Numerical simulation of experimental gravity-driven unstable flow in water repellent sand

Abstract Laboratory experiments related to gravity-driven unstable flows in water repellent porous media contained in two-dimensional chambers have been reported [Bauters, T.W.J., DiCarlo, D.A., Steenhuis, T.S., Parlange, J.-Y., 1998. Preferential flow in water-repellent sands. Soil Sci. Soc. Am. J. 62, 1185–1190]. These experiments demonstrate that water repellency has a significant impact on the stability of flow. As a follow up to these experiments, numerical solutions of the Richards equation for a two-dimensional domain are derived to examine the effect of water repellency on flow characteristics. Of particular interest is the development of gravity-driven unstable flow conditions caused by water repellency. The degree of water repellency of the porous medium is manifested in the water saturation—capillary pressure and water saturation—hydraulic conductivity relationships for the porous medium. To derive the numerical solutions, parameters closely representing the flow domain boundary conditions and the porous medium properties in the experiments of Bauters et al., were employed. In this paper we present the results of simulations for two cases: a water wettable sand and an extremely water repellent sand. The numerical solution for the water wettable sand led to a stable flow condition, while for the water-repellent sand the flow was unstable as manifested by the development of a single finger of flow. A new feature of these modeling results, in comparison to previous modeling results for gravity-driven unstable flow, is that the water pressure inside the finger core is positive. In testing the numerical solutions we compared the solution results to the laboratory results in terms of flow patterns, water pressure at a single reference point, and wetting front velocity. The degree of agreement between the laboratory results and the numerical solutions in terms of these measures is quite good.

Modeling the effects of nonlinear equilibrium sorption on the transport of solute plumes in saturated heterogeneous porous media

Transport of sorbing solutes in 2D steady and heterogeneous flow fields is modeled using a particle tracking random walk technique. The solute is injected as an instantaneous pulse over a finite area. Cases of linear and Freundlich sorption isotherms are considered. Local pore velocity and mechanical dispersion are used to describe the solute transport mechanisms at the local scale. This paper addresses the impact of the degree of heterogeneity and correlation lengths of the log-hydraulic conductivity field as well as negative correlation between the log-hydraulic conductivity field and the log-sorption affinity field on the behavior of the plume of a sorbing chemical. Behavior of the plume is quantified in terms of longitudinal spatial moments: center-of-mass displacement, variance, 95% range, and skewness. The range appears to be a better measure of the spread in the plumes with Freundlich sorption because of plume asymmetry. It has been found that the range varied linearly with the travelled distance, regardless of the sorption isotherm. This linear relationship is important for extrapolation of results to predict behavior beyond simulated times and distances. It was observed that the flow domain heterogeneity slightly enhanced the spreading of nonlinearly sorbing solutes in comparison to that which occurred for the homogeneous flow domain, whereas the spreading enhancement in the case of linear sorption was much more pronounced. In the case of Freundlich sorption, this enhancement led to further deceleration of the solute plume movement as a result of increased retardation coefficients produced by smaller concentrations. It was also observed that, except for plumes with linear sorption, correlation between the hydraulic conductivity and the sorption affinity fields had minimal effect on the spatial moments of solute plumes with nonlinear sorption.

Modeling plume behavior for nonlinearly sorbing solutes in saturated homogeneous porous media

Abstract Transport of a sorbing solute in a two-dimensional steady and uniform flow field is modeled using a particle tracking random walk method. The solute is initially introduced from an instantaneous point source. Cases of linear and nonlinear sorption isotherms are considered. Local pore velocity and mechanical dispersion are used to describe the solute transport mechanisms at the local scale. The numerical simulation of solute particle transport yields the large scale behavior of the solute plume. Behavior of the plume is quantified in terms of the center-of-mass displacement distance, relative velocity of the center-of-mass, mass breakthrough curves, spread variance, and longitudinal skewness. The nonlinear sorption isotherm affects the plume behavior in the following way relative to the linear isotherm: (1) the plume velocity decreases exponentially with time; (2) the longitudinal variance increases nonlinearly with time; (3) the solute front is steepened and tailing is enhanced

Remediation to improve infiltration into compact soils

Urban development usually involves soil compaction through converting large pervious land into developed land. This change typically increases runoff during runoff events and consequently may add to flooding and additional volume of runoff. The wash off of pollutants may also create numerous water quality and environmental problems for receiving waters. To alleviate this problem many municipalities are considering low impact development. One technique to reduce runoff in an urban area is to improve the soil infiltration. This study is specifically undertaken to investigate tilling and compost addition to improve infiltration rate, and to investigate measurement tools to assess the effectiveness of remediated soil. Soil remediation was performed at three sites in an urban area metropolitan area. Each site was divided into three plots: tilled, tilled with compost addition, and a control plot with no treatment. The infiltration effectiveness within each plot was assessed by measuring saturated hydraulic conductivity (K(sat)) using the modified Philip Dunne (MPD) infiltrometer during pre- and post-treatment. In addition, the use of soil bulk density and soil strength as surrogate parameters for K(sat) was investigated. Results showed that deep tillage was effective at reducing the level of soil strength. Soil strength was approximately half that of the control plot in the first six inches of soil. At two of the sites, tilling was also ineffective at improving the infiltration capacity of the soil. The geometric mean of K(sat) was 0.5-2.3 times that of the control plot, indicating little overall improvement. Compost addition was more effective than tilling by reducing the soil strength and compaction and increasing soil infiltration. The geometric mean of K(sat) on the compost plots was 2.7-5.7 times that of the control plot. No strong correlations were observed before remediation between either soil bulk density or soil strength and K(sat). Simulation results showed that the amount of runoff generated from a selection of design storms for treated soil was less than for untreated soil. The results presented in this study may be used as guidance for urban hydromodification and stormwater management plans.

Internal Erosion during Soil Pipeflow: State of the Science for Experimental and Numerical Analysis

Keywords: Ephemeral gully erosion Erodibility Internal erosion Landslides Pipeflow Soil pipes. Abstract. Many field observations have led to speculation on the role of piping in embankment failures, landslides, and gully erosion. However, there has not been a consensus on the subsurface flow and erosion processes involved, and inconsistent use of terms have exacerbated the problem. One such piping process that has been the focus in numerous field observations, but with very limited mechanistic experimental work, is flow through a discrete macropore or soil pipe. Questions exist as to the conditions under which preferential flow through soil pipes results in internal erosion, stabilizes hillslopes by acting as drains, destabilizes hillslopes via pore-pressure buildups, and results in gully formation or reformation of filled-in ephemeral gullies. The objectives of this article are to review discrepancies in terminology in order to represent the piping processes better, to highlight past experimental work on the specific processes of soil pipeflow and internal erosion, and to assess the state-of-the-art modeling of pipeflow and internal erosion. The studies reviewed include those that examined the process of slope stability as affected by the clogging of soil pipes, the process of gullies forming due to mass failures caused by flow into discontinuous soil pipes, and the process of gully initiation by tunnel collapse due to pipes enlarging by internal erosion. In some of these studies, the soil pipes were simulated with perforated tubes placed in the soil, while in others the soil pipes were formed from the soil itself. Analytical solutions of the excess shear stress equation have been applied to experimental data of internal erosion of soil pipes to calculate critical shear stress and erodibility properties of soils. The most common numerical models for pipeflow have been based on Richards’ equation, with the soil pipe treated as a highly conductive porous medium instead of a void. Incorporating internal erosion into such models has proven complicated due to enlargement of the pipe with time, turbulent flow, and episodic clogging of soil pipes. These studies and modeling approaches are described, and gaps in our understanding of pipeflow and internal erosion processes and our ability to model these processes are identified, along with recommendations for future research.