4. Hydraulics

Abstract

A solid state event recorder and reader suitable for use with a tipping bucket raingauge are described. The recorder is a small battery operated unit with very low power consumption, capable of recording 64, 8 bit, words. The portable reader unit is used to interrogate it.     

Hydraulics of runoff and loess loam deposition

On many more or less loamy soils, rill erosion is reported to start on slopes that are equal to or steeper than 2–3°; critical Froude numbers for the start of rill wash on these slopes vary between 2·0 and 3·0. This explains why colluvial deposition often occurs on slopes below 2–3°, when water spreads out at the downslope extremities of the rills. The critical hydraulic conditions for loess loam deposition were tested in the laboratory for slopes of 0·5° and 2°, applying unit-discharges (q) up to 10 cm²/s. It appeared from these experiments that for afterflow, without raindrop impact, deposition starts for critical load concentrations (c_cr) varying between several g/1 and about 60 g/l. Under rain c_cr amounts to a minimum value of 100–125 g/l and it increases when the runoff film becomes thinner. Nevertheless, deposition in pluvial runoff is also possible, as was the case during the Weichselian, according to data from quarries in Belgium and in The Netherlands. A modified Kalinske equation is proposed for c_cr prediction, with the introduction of a typical empirical coefficient C_r and considering such factors as shear stress and mean particle size. Massive sedimentation may occur when it stops raining and afterflow starts, since c_cr values are then much lower. It is shown from the Shields' diagram that loamy suspensions are more sensitive to sedimentation than sands in clear water.     

Hydraulics of recirculating well pairs for ground water remediation

Recirculating well pairs are a proven means of implementing bioremediation and may also be useful for applying other in situ ground water remediation technologies. A bromide tracer test was performed to characterize the hydraulic performance of a recirculating well pair installed at Moffett Field, California. In particular, we estimate two important properties of the recirculating well pair: (1) the fraction of captured water that is recycled between the wells, and (2) the travel-time distribution of ground water in the induced zone of recirculation. We also develop theoretical estimates of these two properties and demonstrate they depend upon a dimensionless pumping rate, denoted xi. The bromide breakthrough curve predicted from theory agrees well with that determined experimentally at Moffett Field. The minimum travel time between the wells is denoted t(min). In theory, t(min) depends inversely on Q, the pumping rate in the recirculating wells, and is proportional to a2, the square of the distance between the wells. Both the experimental and theoretical travel-time distributions indicate that at least half the recirculating water travels between the wells along fast flowpaths (travel time < 2*t(min)). Therefore, when designing recirculating well pairs, engineers should ensure that t(min) will be sufficiently high to allow biologically mediated reactions (or other in situ remediation processes) sufficient time to proceed.     

Hydraulics of sediment-laden sheetflow and the influence of simulated rainfall

Laminar sheetflows, transporting sediment at their capacity rates, both with and without rainfall disturbance, were investigated. Values of flow depth and relative submergence were very small. In the flows without rainfall, measured velocities exceeded the predictions of the smooth-surface, clear-water laminar model by an average of 12 per cent. Reduced flow resistance due to high sediment concentrations may explain this result. Velocities in the rainfall-disturbed flows were not significantly different from the predictions of the smooth-surface, clear-water model, and the velocity reduction due to rainfall was about 12 per cent. Although the uniformity of rainfall intensity under the single-nozzle rainfall simulator is high, variation of momentum and kinetic energy fluxes along the 1-5 m long flume was significant. The rainfall angle of incidence was highly correlated with deviations from expected flow velocities in the upper and lower sections of the flume.     

Estimating Flow Using Tracers and Hydraulics in Synthetic Heterogeneous Aquifers

Regional ground water flow is most usually estimated using Darcy's law, with hydraulic conductivities estimated from pumping tests, but can also be estimated using ground water residence times derived from radioactive tracers. The two methods agree reasonably well in relatively homogeneous aquifers but it is not clear which is likely to produce more reliable estimates of ground water flow rates in heterogeneous systems. The aim of this paper is to compare bias and uncertainty of tracer and hydraulic approaches to assess ground water flow in heterogeneous aquifers. Synthetic two-dimensional aquifers with different levels of heterogeneity (correlation lengths, variances) are used to simulate ground water flow, pumping tests, and transport of radioactive tracers. Results show that bias and uncertainty of flow rates increase with the variance of the hydraulic conductivity for both methods. The bias resulting from the nonlinearity of the concentration–time relationship can be reduced by choosing a tracer with a decay rate similar to the mean ground water residence time. The bias on flow rates estimated from pumping tests is reduced when performing long duration tests. The uncertainty on ground water flow is minimized when the sampling volume is large compared to the correlation length. For tracers, the uncertainty is related to the ratio of correlation length to the distance between sampling wells. For pumping tests, it is related to the ratio of correlation length to the pumping test's radius of influence. In regional systems, it may be easier to minimize this ratio for tracers than for pumping tests.     

In-Well Hydraulics of the Electromagnetic Borehole Flowmeter: Further Studies

The electromagnetic borehole flowmeter (EBF) is finding increasing application as a method for measuring hydraulic conductivity (K) distributions. A recent paper details an experimental/theoretical study of the effect of in-well hydraulics on calculated K distributions based on EBF measurements (Dinwiddie et al. 1999). Results showed that minimizing head loss associated with flow through the meter was the key to producing accurate K values. Using the same experimental procedures, the previous study has been extended to develop data from a larger diameter (1 inch) EBF, and to determine if an EBF can be calibrated effectively without using an inflatable packer to force all flow through the meter annulus. Both experiments were aimed at producing low head loss conditions. Results show that overall calibration can be accomplished in the absence of a packer, which reduces head losses to nonmeasurable levels, and use of the 1-inch EBF with a packer reduces head losses by a factor of 16 when compared with the 0.5-inch EBF studied by Dinwiddie et al. (1999).     

Hydraulics of overland flow influenced by litter incorporation under extreme rainfall

Water erosion on hillslopes is a worldwide environmental problem, which is a rainfall‐induced process, especially extreme rainfall. The great intensity of extreme rainfall strongly enhances the power of overland flow to detach soil and transport sediment. Plant litter is one of the most important constituents of ecosystems that often covers the soil surface and can be incorporated into topsoil. However, little attention has been paid to its effect on flow hydraulics owing to the veiled nature. This study aimed to examine the effects of incorporated litter on the hydraulic properties under extreme rainfall condition. To reach this goal, six litter rates of 0, 0.05, 0.10, 0.20, 0.35, and 0.50 kg m^?2 and four litter types collected from deciduous trees, coniferous trees, shrubs, and herbs were incorporated into topsoil. Then, simulated rainfall experiments were performed on five slope gradients (5°, 10°, 15°, 20°, and 25°) with an extreme rainfall intensity of 80 mm h^?1. The results showed that Froude number and flow velocity of the overland flow decreased, whereas flow resistance increased exponentially with litter incorporation rate. Litter type had an influence on flow hydraulics, which can mainly be attributed to the variations in surface coverage of the exposed litter and the litter morphology. Flow velocity and Darcy–Weisbach coefficient increased markedly with slope gradient. However, the variation of slope gradient did not modify the relationships between flow hydraulics and incorporated litter rate. The random roughness, resulting from heterogeneous erosion due to the uneven protection of surface exposed litter, increased linearly with litter incorporated rate. As rainfall proceeded, flow hydraulics varied with incorporated litter rate and slope gradient complicatedly due to the increases in flow rate and coverage of the exposed litter and the modification of soil surface roughness.     

Hydraulics of interrill overland flow on rough,bare soil surfaces

A set of laboratory experiments on bare, rough soil surfaces was carried out to study the relationship between soil surface roughness and its hydraulic resistance. Existing models relating roughness coefficients to a measure of surface roughness did not predict the hydraulic resistance well for these surfaces. Therefore, a new model is developed to predict the hydraulic resistance of the surface, based on detailed surface roughness data. Roughness profiles perpendicular to the flow are used to calculate the wet cross‐sectional area and hydraulic radius given a certain water level. The algorithm of Savat is then applied to calculate the hydraulic resistance. The value for the equivalent roughness, which is used in the algorithm of Savat, could be predicted from the roughness profiles. Here, the tortuosity of the submerged part of the surface was used, which means that the calculated roughness depends on flow depth. The roughness increased with discharge, due to the fact that rougher parts of the surface became submerged at higher discharges. Therefore, a single measure of surface roughness (e.g. random roughness) is not sufficient to predict the hydraulic resistance. The proposed model allows the extension of the flow over the surface with increasing discharge to be taken into account, as well as the roughness within the submerged part of the surface. Therefore, the model is able to predict flow velocities reasonably well from discharge and roughness data only. Copyright © 2000 John Wiley & Sons, Ltd.     

Numerical Modeling of the Near‐Field Hydraulics of Water Wells

Numerical flow models can be a useful tool for dimensioning water wells and to investigate the hydraulics in their near‐field. Fully laminar flow can be assumed for all models calculated up to the screen. Therefore models can be used to predict—at least qualitatively, neglecting turbulent losses inside the well—the spatial distribution of inflow into the well and the overall hydraulic performance of different combinations of aquifer parameters and technical installations. Models for both horizontal (plan view) and vertical flow (cross section) to wells were calculated for a variety of setups. For the latter, this included variations of hydraulic conductivity of the screen, pump position, and aquifer heterogeneity. Models of suction flow control devices showed that they indeed can homogenize inflow, albeit at the cost of elevated entrance losses.     

Simulation of Groundwater Level in Ephemeral Streams with an Improved Groundwater Hydraulics Model

As a component of arid ecosystems, groundwater plays an important role in plant growth; therefore, it is essential to use deterministic models to reconstruct the process of groundwater level change. Typically, the linearized solution of the one-dimensional (1-D) Boussinesq equation yields acceptable performance in simulating transient conditions over short recharge periods in ephemeral stream systems, but the ability of this solution to simulate multiyear changes in groundwater levels is limited. In this study, an improved groundwater hydraulics (GH-D2) model is built based on the groundwater hydraulics (GH) solution of the 1-D Boussinesq equation to simulate multiyear changes in the groundwater level in ephemeral stream systems. The model is validated in the lower reaches of the Tarim River to simulate groundwater level fluctuations within the scope of influence of the river (300, 500, 750, 1050 m) over a 16-year period (2000 to 2015). To evaluate the performance of the models, the bias, mean absolute error, root mean squared error, Nash-Sutcliffe efficiency (NSE), and coefficient of determination (R²) are calculated. The results show that the improved GH-D2 model, which considers ephemeral streamflow, unsteady flow theory and the delayed response effect of groundwater level changes, performs well in simulating multiyear changes in the groundwater level in the ephemeral stream system. The observed and simulated values of the groundwater level at different river distances are consistent, and the model provides a new basis for multiyear simulations of groundwater level fluctuations in ephemeral stream systems.     

Hydraulics of a partially penetrating well with skin zone in a confined aquifer

A steady/quasi-steady model is developed for predicting flow into a partially penetrating well with skin zone in a confined aquifer overlying an impervious layer. The model takes into account flow through the bottom of the wellbore, finite skin thickness and finite horizontal and vertical extent of the aquifer. Moreover, the solution can be easily extended to include the mixed-type boundary condition at the well face, where a Dirichlet in the form of a specified hydraulic head and a Neumann in the form of zero flux coexist at the same time at different portions of the well face. The validity of the proposed solution is tested by comparing a few results obtained from the developed model with corresponding results obtained by analytical and numerical means. The study shows that, among other factors remaining constant, both the horizontal and vertical extent of an artesian aquifer, thickness of the skin zone, bottom flow and conductivity contrast of the skin and formation zones, play an important part in deciding flow to a well dug in the aquifer, and hence these factors must be considered while analyzing the problem. The model proposed here can be used to estimate skin thickness as well as hydraulic conductivities of the skin and formation zones of a well with skin zone in an artesian aquifer underlain by an impervious layer by utilizing pumping test data falling in the steady or quasi-steady state of a typical pumping test. As the proposed solution is of a general nature in the sense that it can handle, apart from partial penetration and bottom flow, the finite size skin zone and finite horizontal and vertical extent of an artesian aquifer together with the mixed-type boundary condition at the well face, it is hoped that the predictions coming out of the model will be more realistic than those obtained using solutions developed with more stringent assumptions.     

Assessment of Bed Hydraulics and Metal Loadings in a Passive Vertical Flow Bioreactor in Commerce,Oklahoma

A treatment pond, with an engineered bed that served as a passive vertical flow bioreactor (VFBR), was operated as part of a passive sequenced treatment system for the removal of metals from groundwater at the Mayer Ranch in Commerce, Oklahoma. The groundwater was contaminated by mining activities in west Commerce and discharges at this location occurred as artesian springs through improperly abandoned, over-drilled, and cased legacy boreholes. The VFBR operated by establishing reducing conditions in the organic bed of a pond to promote metal sorption and precipitation as sulfides. In order to verify that operations were unhindered by nonuniform flow in the VFBR, an assessment of the flow uniformity in the pond was undertaken using the streambed point velocity probe (SBPVP). The velocity data were independently validated with a water balance. The outflow calculated from the SBPVP data came within 30% of the value suggested by measured inflow rates to the pond, supporting the conclusion that the SBPVP measurements were representative of flow in the VFBR, and that flow through the bed was occurring with a satisfactory level of uniformity. Water flow rates through the reactive bed were found to be up to an order of magnitude greater than those employed in the prior column testing, contributing to metal loading rates (of Cd, Pb, Zn, Ni) estimated to be two orders of magnitude greater than those tested in the columns (4.2 × 10⁴ and 3.2 × 10² mg/m³/d, respectively). However, apparently rapid chemical reactions that likely occurred close to the pond water-sediment interface contributed to the treatment system achieving its design objectives.