Effect of seepage on flow and bedforms dynamics

This study, using an experimental approach, focuses on the effect of downward seepage on a threshold alluvial channel morphology and corresponding turbulent flow characteristics. In all the experiments, we observed that the streamwise time‐averaged velocities and Reynolds shear stresses were increased under the influence of downward seepage. Scales of eddy length and eddy turnover time were significantly increased with the application of downward seepage, leading to sediment transport and initiation of bedforms along the channel length. As the amount of seepage discharge increased, eddy length and turnover time were further increased, causing the development of larger bedforms. It was revealed that the geometry of bedforms was linked with the size of eddies. In this work, statistics of bedform dynamics are presented in terms of multi‐scalar bedforms in the presence of seepage. These multi‐scalar ubiquitous bedforms cast a potential impact on flow turbulence as well as stream bed morphology in channels. We used wavelet to analyse temporally lagged spatial bed elevation profiles that were obtained from a set of laboratory experiments and synchronized the wavelet coefficients with bed elevation fluctuations at different length scales. A spatial cross‐correlation analysis, based on the wavelet coefficients, was performed on these bed elevation datasets to observe the effect of downward seepage on the dynamic behaviour of bedforms at different length scales. It was found that celerity of bedforms reduced with increase in seepage percentage. Bedform celerity was best approximated by a probability density function such as Rayleigh distribution under varying downward seepage. Further, statistical analysis of physical parameters of bedforms ascertained that the reduction in bedform celerity was a result of increased bedform size. Copyright © 2017 John Wiley & Sons, Ltd.     

Bedload transport processes in a braided gravel-bed river model

A 1:50 scale hydraulic model was designed, based on Froude number similarity and using hydrological and sediment data from a small braided gravel-bed river (the North Branch of the Ashburton River, Canterbury, New Zealand). Eighteen experiments were conducted; seven using steady flows, and eleven using unsteady flows. The experiments were carried out in a 20 m × 3 m tilting flume equipped with a continuous sediment feed and an automated data acquisition and control system. In all experiments water at 30°C was used to reduce viscosity-related scale effects. Analyses of the experimental data revealed that bedload transport rates in braided channels are highly variable, with relative variability being inversely related to mean bedload transport rate. Variability was also found to be cyclic with short-term variations being caused by the migration of bedforms. Bedload transport was found to be more efficient under steady flow than under unsteady flow, and it was postulated that this is caused by a tendency for channel form to evolve towards a condition which maximizes bedload transport for the occurring flow. Average bedload transport rate was found to vary with channel form, although insufficient measurements were made to define a relationship.     

Laboratory and Numerical Investigations of Air Sparging Using MTBE as a Tracer

Air sparging experiments were conducted in a laboratory column to investigate air flow and mass transfer behavior in different types of sand at different air injection rates. Methyl tertiary butyl ether (MTBE) was applied as a tracer, and by measuring the volatilization and the mean air content during the experiments, the air flow pattern and its influence on mass transfer were assessed. The experimental results showed large differences among the sand types. In fine sand, the mean air content was high and the volatilization of MTBE was rapid with total recovery after a few hours. In coarse sand, the mean air content was low and the volatilization of MTBE was limited. The results indicate two different air flow distributions. In fine-grained materials, a uniform air distribution can be expected compared to coarse-grained materials where isolated air channels will limit the mass transfer. Afterwards, the experiments were simulated using the numerical multiphase flow code T2VOC, and the results compared to those obtained in the laboratory. The experiments with fine sand were simulated well, while for coarser sand types the volatilization was highly overestimated. The differences between model and laboratory results were mainly attributed to the nonuniformity of the air saturation and the neglection of kinetics in the mass transfer formulation.     

Sheet flow hydrodynamics over a non-uniform sand bed channel

The current study experimentally investigates the flow characteristics and temporal variations in the sheet flow profile of a non-uniform sand bed channel. Experiments were done to explore turbulent structures in the presence of a sheet flow layer with and without seepage. The turbulent events, such as stream wise velocity, Reynolds shear stresses, and turbulence intensities were found to be increasing and vertical velocity was found decreasing with a sheet layer. The presence of a sheet layer also effects the turbulent energy production and energy dissipation. All the turbulence parameters with and without a sheet layer have also been influenced by the presence of downward seepage. The rate of sheet flow movement is increased with seepage, owing to increased turbulence with seepage. The current study used wavelet analysis on temporally lagged spatial bed elevation profiles obtained from a set of laboratory experiments and synchronized the wavelet coefficients with bed elevation fluctuation at different spatial scales. A spatial cross correlation analysis at multiple scales, based on the wavelet coefficients, has been done on these bed elevation datasets to observe the effect of downward seepage on the dynamic behavior of sheet flow at different length scales. It is found that seepage increases average bed celerity and also increases the celerity of sheet flow of similar length scales. This increase in the celerity has been hypothesized as the increase of sheet flow movement as well as the increase in turbulent parameters with seepage, which destabilizes the bed particles resulting in a disruption in the continuous propagation pattern of the sheet flow. The increase of sheet flow celerity with seepage is confirmed from the saturation level of the wavelet power spectra of the bed elevation series. The presence of seepage also affects the non-uniformity of collective sheet material.     

Natural degradation rates of BTEX compounds and naphthalene in a sulphate reducing groundwater environment

Abstract

Field and laboratory evidence show natural degradation of toluene, ethylbenzene, m-, p- and o-xylenes, 1,3,5-trimethylbenzene and naphthalene in sulphate reducing groundwater conditions of the Bassendean Sands in the Perth basin, Western Australia. Natural degradation rates were obtained from a groundwater tracer test with deuterated organic compounds injected into a dissolved hydrocarbon plume, down-gradient of a leaking underground storage tank at an urban service station. These were compared with similar data obtained from modelling of the whole contaminant plume itself and also with data obtained from large-scale laboratory column experiments with groundwater spiked with BTEX compounds. Toluene degradation rate was 200 to 500 times higher in the anaerobic laboratory columns than in the field. Degradation rates in the tracer test compared well with model-derived field estimates.     

Studies on emergent flow over vegetative channel bed with downward seepage

Experimental observations in a tilting flume having a bed covered with rice plants (Oryza sativa) are used to analyse the flow characteristics of flexible emergent vegetation with downward seepage. The flow velocity for no-seepage and with seepage is reduced by, on average, 52% and 33%, respectively, as the flow reaches the downstream end with vegetation. Higher Reynolds stress occurs at the start of the vegetation zone; hence, bed material transport occurs in this region. The results indicate that the bed is no longer the primary source of turbulence generation in vegetated flow; rather it is dominated by turbulence generated by the vegetation stems. The local effect of the presence of vegetation causes variations in the hydrodynamic characteristics along the vegetated portion of the channel, which leads to erosion and deposition in the vegetation zone. The experiments show that vegetation can provide considerable stability to channels by reducing channel erosion even with downward seepage.     

Bank undercutting and tension failure by groundwater seepage: predicting failure mechanisms

Groundwater seepage can lead to the erosion and failure of streambanks and hillslopes. Two groundwater instability mechanisms include (i) tension failure due to the seepage force exceeding the soil shear strength or (ii) undercutting by seepage erosion and eventual mass failure. Previous research on these mechanisms has been limited to non‐cohesive and low cohesion soils. This study utilized a constant‐head, seepage soil box packed with more cohesive (6% and 15% clay) sandy loam soils at prescribed bulk densities (1.30 to 1.70 Mg m^?3) and with a bank angle of 90° to investigate the controls on failure mechanisms due to seepage forces. A dimensionless seepage mechanism (SM) number was derived and evaluated based on the ratio of resistive cohesion forces to the driving forces leading to instability including seepage gradients with an assumed steady‐state seepage angle. Tension failures and undercutting were both observed dependent primarily on the saturated hydraulic conductivity, effective cohesion, and seepage gradient. Also, shapes of seepage undercuts for these more cohesive soils were wider and less deep compared to undercuts in sand and loamy sand soils. Direct shear tests were used to quantify the geotechnical properties of the soils packed at the various bulk densities. The SM number reasonably predicted the seepage failure mechanism (tension failure versus undercutting) based on the geotechnical properties and assumed steady‐state seepage gradients of the physical‐scale laboratory experiments, with some uncertainty due to measurement of geotechnical parameters, assumed seepage gradient direction, and the expected width of the failure block. It is hypothesized that the SM number can be used to evaluate seepage failure mechanisms when a streambank or hillslope experiences steady‐state seepage forces. When prevalent, seepage gradient forces should be considered when analyzing bank stability, and therefore should be incorporated into commonly used stability models. Copyright © 2013 John Wiley & Sons, Ltd.     

Laboratory and numerical investigation of flow and transport near a seepage-face boundary

Laboratory and numerical modeling investigations were completed to study the unconfined ground water flow and transport processes near a seepage-face boundary. The laboratory observations were made in a radial sand tank and included measurements of the height of the seepage face, flow velocity near the seepage face, travel time distribution of multiple tracer slugs, and streamlines. All the observations were reliably reproduced with a three-dimensional, axi-symmetric, variably saturated ground water flow model. Physical data presented in this work demonstrate and quantify the importance of three-dimensional transport patterns within a seepage-face zone. The results imply that vertically averaged flow models that employ Dupuit approximations might introduce error in the analysis of localized solute transport near a seepage-face boundary. The experimental dataset reported in this work will also be of interest for those who are attempting to validate a numerical algorithm for solving ground water and contaminant discharge patterns near a surface-water boundary.     

The emplacement of pahoehoe toes: field observations and comparison to laboratory simulations

We observed active pahoehoe lobes erupted on Kilauea during May-June 1996, and found a range of emplacement styles associated with variations in local effusion rate, flow velocity, and strain rate. These emplacement styles were documented and quantified for comparison with earlier laboratory experiments.At the lowest effusion rates, velocities, and strain rates, smooth-surfaced lobes were emplaced via swelling, where new crust formed along an incandescent lip at the front of the lobe and the rest of the lobe was covered with a dark crust. At higher effusion rates, strain rates and velocities, lobes were emplaced through tearing or cracking. Tearing was characterized by ripping of the ductile crust near the initial breakout point, and most of the lobe surface was incandescent during its emplacement. This mechanism was observed to generate both smooth-surfaced lobes, and, when the lava encountered an obstacle, folded lobes. Cracking lobes were similar to those emplaced via tearing, but involved breaking of a thicker, brittle crust at the initial breakout of the lobe and therefore required somewhat higher flow rates than did tearing. Cracking lobes typically formed ropy folds in the center of the lobe, and smooth margins. At the highest effusion rates, strain rates, and flow velocities, the lava formed open channels with distinct levees.The final lobe morphologies were compared to results from laboratory simulations, which were designed to infer effusion rate from final flow morphology, to quantitatively test the laboratory results on the scale of individual natural pahoehoe lobes. There is general agreement between results from laboratory simulations and natural lavas on the scale of individual pahoehoe lobes, but there are disparities between laboratory flows and lava flows on the scale of an entire pahoehoe lava flow field.Editorial responsibility: A. Woods     

An analytical solution for predicting the transient seepage from a subsurface drainage system

Subsurface drainage systems have been widely used to deal with soil salinization and waterlogging problems around the world. In this paper, a mathematical model was introduced to quantify the transient behavior of the groundwater table and the seepage from a subsurface drainage system. Based on the assumption of a hydrostatic pressure distribution, the model considered the pore-water flow in both the phreatic and vadose soil zones. An approximate analytical solution for the model was derived to quantify the drainage of soils which were initially water-saturated. The analytical solution was validated against laboratory experiments and a 2-D Richards equation-based model, and found to predict well the transient water seepage from the subsurface drainage system. A saturated flow-based model was also tested and found to over-predict the time required for drainage and the total water seepage by nearly one order of magnitude, in comparison with the experimental results and the present analytical solution. During drainage, a vadose zone with a significant water storage capacity developed above the phreatic surface. A considerable amount of water still remained in the vadose zone at the steady state with the water table situated at the drain bottom. Sensitivity analyses demonstrated that effects of the vadose zone were intensified with an increased thickness of capillary fringe, capillary rise and/or burying depth of drains, in terms of the required drainage time and total water seepage. The analytical solution provides guidance for assessing the capillary effects on the effectiveness and efficiency of subsurface drainage systems for combating soil salinization and waterlogging problems.     

Interferometric observations of the temperature structure in water cooled or heated from above

Controlled laboratory experiments were performed in a water tank to study energy transfer during cooling and radiant reheating of water from above. Accurate instantaneous temperature measurements were obtained using a Mach-Zehnder interferometer which did not disturb the temperature, radiation, and flow fields. The buoyancy induced flow field was visualized by employing an electro-chemical dye production technique. Cooling of initially uniform layer of water and stratification-cooling-restratification were studied.The experimental observation during cooling of an initially uniform temperature column of water indicated intermittent free convection near the surface. The phenomenon was characterized by randomly descending cooler plumes penetrating 10 to 25 cm into the warmer underlying region. The cooling rate had a decisive influence on the frequency and intensity of the descending parcels of water. For lower cooling rates the colder water usually descended in a form of sheets.Observations during cooling of initially stratified water showed that buoyancy driven convection occurred near the surface. The motion had an initial regular roll pattern, but as the cooling continued the roll pattern in the convective (mixed) layer deteriorated into motion similar to that observed during cooling of an initially uniform temperature layer of water. When the water was restratified by radiant heating, the circulation in the convective layer was suppressed, and the temperature profiles obtained were similar to those observed in a fluid during free convection between two solid parallel walls.     

Effects of vegetation on channel morphodynamics: results and insights from laboratory experiments

A series of laboratory experiments demonstrates that riparian vegetation can cause a braided channel to self‐organize to, and maintain, a dynamic, single‐thread channel. The initial condition for the experiments was steady‐state braiding in non‐cohesive sand under uniform discharge. From here, an experiment consisted of repeated cycles alternating a short duration high flow with a long duration low flow, and uniform dispersal of alfalfa seeds over the bed at the end of each high flow. Plants established on freshly deposited bars and areas of braidplain that were unoccupied during low flow. The presence of the plants had the effect of progressively focusing the high flow so that a single dominant channel developed. The single‐thread channel self‐adjusted to carry the high flow. Vegetation also slowed the rate of bank erosion. Matching of deposition along the point bar with erosion along the outer bend enabled the channel to develop sinuosity and migrate laterally while suppressing channel splitting and the creation of new channel width. The experimental channels spontaneously reproduced many of the mechanisms by which natural meandering channels migrate and maintain a single dominant channel, in particular bend growth and channel cutoff. In contrast with the braided system, where channel switching is a nearly continuous process, vegetation maintained a coherent channel until wholesale diversion of flow via cutoff and/or avulsion occurred, by which point the previous channel tended to be highly unfavorable for flow. Thus vegetation discouraged the coexistence of multiple channels. Varying discharge was key to allowing expression of feedbacks between the plants and the flow and promoting the transition from braiding to a single‐thread channel that was then dynamically maintained. Copyright © 2010 John Wiley & Sons, Ltd.     

CONTROL OF GULLY EROSION USING STIFF GRASSES

Concentrated flow can cause gully formation on sloping lands and in riparian zones. Current practice for riparian gully erosion control involves blocking the gully with a structure comprised of an earthen embankment and a metal or plastic pipe. Measures involving native vegetation would be more attractive for habitat recovery and economic reasons. To test the hypothesis that switchgrass （Panicum virgatum L.） hedges planted at 0.5-m vertical intervals within a gully would control erosion, a series of hedges was established in four concentrated flow channels. Two of the channels were previously eroded trapezoidal channels cut into compacted fill in an outdoor laboratory. The other two channels were natural gullies located at the edge of floodplain fields adjacent to an incised stream. While vegetation was dormant, artificial runoff events were created in the two laboratory gullies and one of the natural gullies using synthetic trapezoidal-shaped hydrographs with peak discharge rates of approximately 0.03, 0.07, and 0.16 m3/s. During these tests flow depth, velocity, turbidity, and soil pore water pressures were monitored. The fourth gully was subjected to a series of natural runoff events over a five-month period with peaks up to 0.09 m3/s. Flow depths in all tests were generally 〈 0.3 m, and flow velocities varied spatially and exceeded 2.0 m/s at the steepest points of the gullies. Erosion rates were negligible for controlled flow experiments, but natural flows in the fourth gully resulted in 1 m ofthalweg degradation, destroying the central portions of the grass hedges, most likely due to the highly erodible nature of the soils at this site. Geotechnical modeling of soil steps reinforced with switchgrass roots showed factors of safety 〉 1 for step heights 〈 0.5 m, but instability was indicated for step heights 〉1 m, consistent with the experimental observations.     

Investigation of Seepage Meter Measurements in Steady Flow and Wave Conditions

Water exchange between surface water and groundwater can modulate or generate ecologically important fluxes of solutes across the sediment‐water interface. Seepage meters can directly measure fluid flux, but mechanical resistance and surface water dynamics may lead to inaccurate measurements. Tank experiments were conducted to determine effects of mechanical resistance on measurement efficiency and occurrence of directional asymmetry that could lead to erroneous net flux measurements. Seepage meter efficiency was high (average of 93%) and consistent for inflow and outflow under steady flow conditions. Wave effects on seepage meter measurements were investigated in a wave flume. Seepage meter net flux measurements averaged 0.08 cm/h—greater than the expected net‐zero flux, but significantly less than theoretical wave‐driven unidirectional discharge or recharge. Calculations of unidirectional flux from pressure measurements (Darcy flux) and theory matched well for a ratio of wave length to water depth less than 5, but not when this ratio was greater. Both were higher than seepage meter measurements of unidirectional flux made with one‐way valves. Discharge averaged 23% greater than recharge in both seepage meter measurements and Darcy calculations of unidirectional flux. Removal of the collection bag reduced this net discharge. The presence of a seepage meter reduced the amplitude of pressure signals at the bed and resulted in a nearly uniform pressure distribution beneath the seepage meter. These results show that seepage meters may provide accurate measurements of both discharge and recharge under steady flow conditions and illustrate the potential measurement errors associated with dynamic wave environments.     

Three‐dimensional numerical modeling of the Bulle effect: the nonlinear distribution of near‐bed sediment at fluvial diversions

The Bulle effect is a phenomenon in which a disproportionately higher amount of near‐bed sediment load at a fluvial diversion moves into the diverted channel, even for cases in which the proportion of water (with respect to the main flow) entering the diversion channel is relatively small. This phenomenon has wide‐ranging implications for both engineered and natural systems: from efficient design of channels to redirect water and sediment for reclaiming sinking deltas, designing navigational channels that do not need frequent dredging, to morphological evolution of river bifurcations. The first ever, and one of the most extensive set of experiments conducted to explore this phenomenon, were conducted by Bulle in 1926 . In the current study the experiments conducted by Bulle have been simulated using an open‐source, free‐surface finite‐element‐based hydrodynamic solver. The main objectives were to explore to what extent the complex phenomenon of the Bulle effect at the scale of a laboratory experiment can be simulated accurately using Reynolds‐averaged Navier–Stokes (RANS)‐based hydrodynamic solver, and to understand the details of the hydrodynamics that Bulle could not analyze through his experiments. The hydrodynamics captured by the simulations were found to match the observations made by Bulle through his experiments, and the distributions of sediment at the diversion predicted by the numerical simulations were found to match the general trend observed in the laboratory experiments. The results from the numerical simulations were also compared with existing one‐dimensional models for sediment distribution at bifurcations, and the three‐dimensional numerical model was found to perform appreciably better. This is expected due to the complex flow features at the diversion, which can only be captured satisfactorily using a three‐dimensional hydrodynamic model. Copyright © 2017 John Wiley & Sons, Ltd.     

A laboratory simulation of pyroclastic flows down slopes

Laboratory experiments are described which explore the dynamical consequences of buoyant convective upflow observed above hot pyroclastic flows. In nature, the convection is produced by the hot ash particles exchanging heat with air mixed into the front and top of the pyroclastic flow. This effect on the buoyancy due to the mixing of air and ash has been modelled in the laboratory using mixtures of methanol and ethylene glycol (MEG), which have a nonlinear density behaviour when mixed with water. Intermediate mixtures of these fluids can be denser than either initial component, and so the laboratory experiments were inverted models of the natural situation. We studied MEG flowing up under a sloping roof in a tank filled with water. The experiments were performed both in a narrow channel and on a laterally unconfined slope. The flow patterns were also compared with those of conventional gravity currents formed using fresh and salt water. The presence of the region of reversed buoyancy outside the layer flowing along the slope had two significant effects. First, it periodically protected the flow from direct mixing with the environment, resulting in pulses of relatively undiluted fluid moving out intermittently ahead of the main flow. Second, it produced a lateral inflow towards the axis of the current which kept the current confined to a narrow tongue, even on a wide slope.In pyroclastic flows the basal avalanche portion has a much larger density contrast with its surroundings than the laboratory flows. Calculations show that mixing of air into the dense part of a pyroclastic flow cannot generate a mixture that is buoyant in the atmosphere. However, the overlying dilute ash cloud can behave as a gravity current comparable in density contrast to the laboratory flows and can become buoyant, depending on the temperature and ash content. In the August 7th pyroclastic flow of Mount St. Helens, Hoblitt (1986) describes pulsations in the flow front, which are reminiscent of those observed in the experiments. As proposed by Hoblitt, the pulsations are caused by the ash cloud accelerating away from the front of the dense avalanche as a density current. The ash cloud then mixes with more air, becomes buoyant and lifts off the ground, allowing the avalanche to catch up with and move ahead of the cloud. The pulsing behaviour at the fronts of pyroclastic flows could account for the occurrence of cross-bedded layer 1 deposits which occur beneath layer 2 deposits in many sequences.     

Thermal Imagery of Groundwater Seeps: Possibilities and Limitations

Quantifying groundwater flow at seepage faces is crucial because seepage faces influence the hydroecology and water budgets of watersheds, lakes, rivers and oceans, and because measuring groundwater fluxes directly in aquifers is extremely difficult. Seepage faces provide a direct and measurable groundwater flux but there is no existing method to quantitatively image groundwater processes at this boundary. Our objective is to determine the possibilities and limitations of thermal imagery in quantifying groundwater discharge from discrete seeps. We developed a conceptual model of temperature below discrete seeps, observed 20 seeps spectacularly exposed in three dimensions at an unused limestone quarry and conducted field experiments to examine the role of diurnal changes and rock face heterogeneity on thermal imagery. The conceptual model suggests that convective air‐water heat exchange driven by temperature differences is the dominant heat transfer mechanism. Thermal imagery is effective at locating and characterizing the flux of groundwater seeps. Areas of active groundwater flow and ice growth can be identified from thermal images in the winter, and seepage rates can be differentiated in the summer. However, the application of thermal imagery is limited by diverse factors including technical issues of image acquisition, diurnal changes in radiation and temperature, and rock face heterogeneity. Groundwater discharge rates could not be directly quantified from thermal imagery using our observations but our conceptual model and experiments suggest that thermal imagery could quantify groundwater discharge when there are large temperature differences, simple cliff faces, non‐freezing conditions, and no solar radiation.     

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework.The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy–sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements.The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux.As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers.     

Seepage face height, water table position, and well efficiency at steady state

When a fully penetrating well pumps an ideal unconfined aquifer at steady state, the water table usually does not join the water level in the well. There is a seepage face inside the well, which is a key element in evaluating the well performance. This problem is analyzed using the finite-element method, solving the complete equations for saturated and unsaturated flow. The seepage face position is found to be almost independent of the unsaturated zone properties. The numerical results are used to test the validity of several analytic approximations. Equations are proposed to predict the seepage face position at the pumping well for any well drawdown, and the water table position at any distance from the pumping well for any in-well drawdown. Practical hints are provided for installing monitoring wells and evaluating well efficiency.     

Hillside seepage and the steady water table. II: Applications

Solutions for steady hillside seepage depend on an accurate estimate of the water table location, before any of the flow parameters can be quantified. In this paper, solutions are obtained for hillslopes subjected to low recharge, using the series methods described in a companion paper (Adv. Water Resour., 19(2) (1996) 63–73). The recharge rate is bounded above, when the upstream boundary is not a vertical, impermeable dyke. The practical consequences of exceeding the maximum recharge rate are examined, for realistic hillslope geometry and recharge rates, and some modifications to the theory are suggested.