Terrestrial laser scanning soil surfaces: a field methodology to examine soil surface roughness and overland flow hydraulics

Conventional roughness–resistance relationships developed for pipe and open‐channel flows cannot accurately describe shallow overland flows over natural rough surfaces. This paper develops a new field methodology combining terrestrial laser scanning (TLS) and overland flow simulation to provide a high‐resolution dataset of surface roughness and overland flow hydraulics as simulated on natural bare soil surfaces. This method permits a close examination of the factors controlling flow velocity and a re‐evaluation of the relationship between surface roughness and flow resistance. The aggregate effect of flow dynamics, infiltration and depression storage on retarding the passage of water over a surface is important where runoff‐generating areas are distant from well‐defined channels. Experiments to separate these effects show that this ‘effective resistance’ is dominated by surface roughness. Eight measurements of surface roughness are found to be related to flow resistance: standard deviation of elevations, inundation ratio, pit density (measured both perpendicular and parallel to the flow direction), slope, median depth, skewness of the depth distribution and frontal area. Hillslope position is found to affect the significant roughness measures. In contrast, infiltration rate has little effect on the velocity of water fronts advancing over the soil surfaces examined here and the effect of depression storage is limited. Overland flow resistance is depth dependent where complex microtopographic structures are progressively inundated. Copyright © 2010 John Wiley & Sons, Ltd.     

Preferential percolation quantified by large water content sensors with artifactual macroporous envelopes

For many scientific and practical tasks, it is important to estimate the soil–water percolation fluxes. This paper builds on measurements with large horizontal time‐domain reflectometry water content sensors in a loamy Mollisol. The sensors were installed into pre‐drilled holes and the gaps between them, and the soil was filled with a slurry of local soil with water. This gave rise to envelopes around them that contained artificial macropores. The sensors reacted to intensive rains by a rapid increase of their readings, often above the native soil's porosity, followed by an almost equally rapid decrease. The paper explores the feasibility of quantifying the rapid percolation, based on these anomalous water content peaks, and demonstrates that this is possible in principle, if the processes are simulated by a suitable model. A two‐dimensional dual porosity non‐equilibrium (mobile‐immobile) model was tried. The envelope around the sensor was modelled as an annulus with higher porosity and hydraulic conductivity, which attracts preferential flow and amplifies the percolation signal. With the model at hand, the flux hydrographs can be derived from model simulations and measured precipitation. For contrast, the Durner equilibrium dual porosity model was tried but was found little suitable. However, even the mobile‐immobile model did not perform perfectly. Simulated water contents were similar to the measured ones at some depths but not in the others, and the percolation fluxes were overestimated, compared to cumulative soil–water balance. Efforts to improve model performance were not successful. Hence, the model structure needs to be improved. Copyright © 2015 John Wiley & Sons, Ltd.     

Evolution of overland flow connectivity in bare agricultural plots

Soil surface roughness not only delays overland flow generation but also strongly affects the spatial distribution and concentration of overland flow. Previous studies generally aimed at predicting the delay in overland flow generation by means of a single parameter characterizing soil roughness. However, little work has been done to find a link between soil roughness and overland flow dynamics. This is made difficult because soil roughness and hence overland flow characteristics evolve differently depending on whether diffuse or concentrated erosion dominates. The present study examined whether the concept of connectivity can be used to link roughness characteristics to overland flow dynamics. For this purpose, soil roughness of three 30‐m² tilled plots exposed to natural rainfall was monitored for two years. Soil micro‐topography was characterized by means of photogrammetry on a monthly basis. Soil roughness was characterized by the variogram, the surface stream network was characterized by network‐based indices and overland flow connectivity was characterized by Relative Surface Connection function (RSCf) functional connectivity indicator. Overland flow hydrographs were generated by means of a physically‐based overland flow model based on 1‐cm resolution digital elevation models. The development of eroded flow paths at the soil surface not only reduced the delay in overland flow generation but also resulted in a higher continuity of high flow velocity paths, an increase in erosive energy and a higher rate of increase of the overland flow hydrograph. Overland flow dynamics were found to be highly correlated to the RSCf characteristic points. By providing information regarding overland flow dynamics, the RSCf may thus serve as a quantitative link between soil roughness and overland flow generation in order to improve the overland flow hydrograph prediction. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling investigation of water partitioning at a semiarid ponderosa pine hillslope

The effects of vegetation root distribution on near‐surface water partitioning can be two‐fold. On the one hand, the roots facilitate deep percolation by root‐induced macropore flow; on the other hand, they reduce the potential for deep percolation by root‐water‐uptake processes. Whether the roots impede or facilitate deep percolation depends on various conditions, including climate, soil, and vegetation characteristics. This paper examines the effects of root distribution on deep percolation into the underlying permeable bedrock for a given soil profile and climate condition using HYDRUS modelling. The simulations were based on previously field experiments on a semiarid ponderosa pine (Pinus ponderosa) hillslope. An equivalent single continuum model for simulating root macropore flow on hillslopes is presented, with root macropore hydraulic parameterization estimated based on observed root distribution. The sensitivity analysis results indicate that the root macropore effect dominates saturated soil water flow in low conductivity soils (K_matrix below 10^?7 m/s), while it is insignificant in soils with a K_matrix larger than 10^?5 m/s, consistent with observations in this and other studies. At the ponderosa pine site, the model with simple root‐macropore parameterization reasonably well reproduces soil moisture distribution and some major runoff events. The results indicate that the clay‐rich soil layer without root‐induced macropores acts as an impeding layer for potential groundwater recharge. This impeding layer results in a bedrock percolation of less than 1% of the annual precipitation. Without this impeding layer, percolation into the underlying permeable bedrock could be as much as 20% of the annual precipitation. This suggests that at a surface with low‐permeability soil overlying permeable bedrock, the root penetration depth in the soil is critical condition for whether or not significant percolation occurs. Copyright © 2010 John Wiley & Sons, Ltd.     

Infiltration,runoff and sediment production in blanket peat catchments: implications of field rainfall simulation experiments

Blanket peat covers the headwaters of many major European rivers. Runoff production in upland blanket peat catchments is flashy with large flood peaks and short lag times; there is minimal baseflow. Little is known about the exact processes of infiltration and runoff generation within these upland headwaters. This paper presents results from a set of rainfall simulation experiments performed on the blanket peat moorland of the North Pennines, UK. Rainfall was simulated at low intensities (3–12 mm h^?1), typical of natural rainfall, on bare and vegetated peat surfaces. Runoff response shows that infiltration rate increases with rainfall intensity; the use of low‐intensity rainfall therefore allows a more realistic evaluation of infiltration rates and flow processes than previous studies. Overland flow is shown to be common on both vegetated and bare peat surfaces although surface cover does exert some control. Most runoff is produced within the top few centimetres of the peat and runoff response decreases rapidly with depth. Little vertical percolation takes place to depths greater than 10 cm owing to the saturation of the peat mass. This study provides evidence that the quickflow response of upland blanket peat catchments is a result of saturation‐excess overland flow generation. Rainfall–runoff response from small plots varies with season. Following warm, dry weather, rainfall tends to infiltrate more readily into blanket peat, not just initially but to the extent that steady‐state surface runoff rates are reduced and more flow takes place within the peat, albeit at shallow depth. Sediment erosion from bare peat plots tends to be supply limited. Seasonal weather conditions may affect this in that after a warm, dry spell, surface desiccation allows sediment erosion to become transport limited. Copyright © 2002 John Wiley & Sons, Ltd.     

High spatial resolution mapping of surface velocities and depths for shallow overland flow

Point measurements of flow rate, depth or velocity are not sufficient to validate overland flow models, particularly when the interaction of the water with the soil surface creates a complex flow geometry. In this study, we present the coupling of two techniques obtaining spatial data of flow depths and surface velocity measurements for water depths as low as 1 mm. Overland flow experiments were performed in the laboratory at various flow rates and slopes on two surfaces. The first surface was 120 cm by 120 cm showing three undulations of sinusoidal shape with an amplitude of 1 cm and a wavelength of 20 cm, while the second was a 60 cm by 60 cm moulded reproduction of a seedbed with aggregates up to 2 cm in size. Large scale particle image velocimetry (LSPIV) was used for velocity measurements with a sub‐centimetre spatial resolution. An instantaneous‐profile laser scanner was used to map flow depths with a sub‐millimetre spatial resolution. A sensitivity analysis of the image processing of the LSPIV showed good robustness of the method. Comparison with measurements performed with hot film anemometer and salt velocity gauge showed that LSPIV surface velocities were representative of the flow. Water depths measured with the laser scanner were also in good agreement with single‐point measurements performed with a dial indicator. Spatially‐distributed flow rates could be computed by combining both presented techniques with a mean relative error less than 20%. Copyright © 2012 John Wiley & Sons, Ltd.     

Overland flow generation and continuity on stone-covered soil surfaces

The influence of stone cover on the generation and continuity of overland flow is a function of several variables, primarily stone size, distance (spacing) between stones, and stone position (on top of the soil surface or partially embedded). The initial hypothesis of the present study were that stone size affects overland flow generation by inducing ‘concentration overland flow’ and that the distance between stones affects overland flow continuity. With respect to stone size and distance between stones, opposing results have been found in the literature. Accordingly, the present paper reports on laboratory experiments which were conducted to investigate, systematically, the effect on overland flow of stone size, distance between stones, and stone position. The main conclusions were:

• 1 Stone cover tended to induce overland flow, relative to bare soil. Small stones, however, especially in a low-cover percentage setting, and when resting on top of the soil surface, produced less overland flow than bare soil.
• 2 Overland flow was positively related to stone size, but inversely related to distance between stones.
• 3 Overland flow yield was always greater when stones were embedded than when on top of the soil surface.
• 4 The effect of stone size on increasing the degree of overland flow was of greater significance than the effect of the distance between stones on reducing overland flow.

     

Hydrologic procedures of storm event watershed models: a comprehensive review and comparison

Currently, many watershed models are available that have various complexities, strengths, and weaknesses. The basic mathematical foundations of these mathematical models are often overlooked due to high demands on convenient applications with graphical user interfaces. Although this and other factors are important while selecting a model, the mathematical foundation should also be taken into account, as performance or efficiency and accuracy of a model depend on its simplicity or complexity. A comprehensive review of 14 storm event watershed models was conducted. Hydrologic procedures (rainfall excess, flow routing, and subsurface flow) of the models are presented and compiled. Among the procedures, flow routing has the most influence on model performances (speed and accuracy). Overland and channel flow routing procedures using different flow‐governing equations, having various approximations and solved by different methods, are compared based on their relative levels of physical bases, complexities, and expected accuracies in simulating the dynamics of water flow. Models using more mathematical terms in the flow‐governing equations are more physically based and expected to be more accurate than models using approximations, however, are more complex due to more intensive but approximate numerical schemes (inefficient). Models using approximate equations with analytical solutions may provide a balance between complexity and accuracy. The review and comparisons are useful to modellers, water resources managers, and researchers in understanding the basic foundations of the models and making informed selections for practical applications or further developments. Other factors such as data intensiveness, user friendliness, and resource requirements are also important considerations. Copyright © 2011 John Wiley & Sons, Ltd.     

The impact of ditch blocking on the hydrological functioning of blanket peatlands

Ditch blocking in blanket peatlands is common as part of peatland restoration. The effects of ditch blocking on flow regimes and nearby water tables were examined in a field trial. After an initial 6‐month monitoring period, eight ditches had peat dams installed 10 m apart along their entire length (dammed), four of these ditches were also partially infilled through bank reprofiling (reprofiled). Four ditches were left open with no dams or reprofiling (open). These 12 ditches and the surrounding peat were monitored for 4 more years. An initial five‐fold reduction in discharge occurred in the dammed and the reprofiled ditches with the displaced water being diverted to overland flow and pathways away from the ditches. However, there was a gradual change over time in ditch flow regime in subsequent years, with the overall volume of water leaving the dammed and the reprofiled ditches increasing per unit of rainfall to around twice that which occurred in the first year after blocking. Hence, monitoring for greater than one year is important for understanding hydrological impacts of peatland restoration. Overland flow and flow in the upper ~4 cm of peat was common and occurred in the inter‐ditch areas for over half of the time after ditch blocking. There was strong evidence that topographic boundaries of small ditch catchments, despite being defined using a high‐resolution Light Detection And Ranging‐based terrain model, were not always equivalent to actual catchment areas. Hence, caution is needed when upscaling area‐based fluxes, such as aquatic carbon fluxes, from smaller scale studies including those using ditches and small streams. The effect of ditch blocking on local water tables was spatially highly variable but small overall (time‐weighted mean effect <2 cm). Practitioners seeking to raise water tables through peatland restoration should first be informed either by prior measurement of water tables or by spatial modelling to show whether the peatland already has shallow water tables or whether there are locations that could potentially undergo large water‐table recoveries.     

Effect of snow structure on water flow and solute transport

Water flow through a melting snow pack modifies its structure and stability and affects the release of water and nutrients into soils and surface waters. Field and laboratory observations indicate a large spatial variability on various scales of the liquid water content and flow, a dominant system feature currently not included in numerical models. We investigated experimentally water and dye tracer movement through microstructurally different snow pack horizons and the persistence of preferential flow paths. Naturally rounded snow of varying grain size was artificially packed to obtain well known conditions by sieving it into rectangular bins. Surface melt was induced with infrared lamps. The flow paths were visualized with tracers and liquid water content was monitored with time domain reflectometry probes. Vertical cuts through the snow pack were imaged. The dye tracer patterns allowed the two flow regimes ‘matrix flow’ and ‘preferential flow’ to be distinguished. Matrix flow is apparently dominated by film and capillary flow in the unsaturated snow matrix. The capillary barrier effect at a boundary between a fine over a coarse textured layer on matrix flow in snow was confirmed. In contrast, preferential flow appears as well‐defined flow fingers that advance from 0·1 to 1 cm s^?1. During a melt phase, the advancing flow fingers enlarge and are only partially time invariant. It remains to be shown whether the continuum concept, including the Darcy–Buckingham law is apt to describe the extremely non‐linear nature of water flow and the travel time of solutes in snow under conditions of melt water percolation. Probably, snow packs that include faceted crystals and large variations in bulk density, feature more pronounced capillary barriers and preferential flow triggering, but also stronger impeding of fingers by lateral dispersion. Further, triggering and persistence of preferential flow is complicated by the usually transient infiltration rate. Copyright © 2004 John Wiley & Sons, Ltd.     

A new modeling approach for simulating microtopography-dominated,discontinuous overland flow on infiltrating surfaces

Realistic modeling of discontinuous overland flow on irregular topographic surfaces has been proven to be a challenge. This study is aimed to develop a new modeling framework to simulate the discontinuous puddle-to-puddle (P2P) overland flow dynamics for infiltrating surfaces with various microtopographic characteristics. In the P2P model, puddles were integrated in a well-delineated, cascaded drainage system to facilitate explicit simulation of their dynamic behaviors and interactions. Overland flow and infiltration were respectively simulated by using the diffusion wave model and a modified Green–Ampt model for the DEM-derived flow drainage network that consisted of a series of puddle-based units (PBUs). The P2P model was tested by using a series of data from laboratory overland flow experiments for various microtopography, soil, and rainfall conditions. The modeling results indicated that the hierarchical relationships and microtopographic properties of puddles significantly affected their connectivity, filling–spilling dynamics, and the associated threshold flow. Surface microtopography and rainfall characteristics also exhibited strong influences on the spatio-temporal distributions of infiltration rates, runoff fluxes, and unsaturated flow. The model tests demonstrated its applicability in simulating microtopography-dominated overland flow on infiltrating surfaces.     

Percolation losses in paddy fields with a dynamic soil structure: model development and applications

The hydraulic characteristics of the plough pan of paddy fields provide continuous ponding conditions during the growing season and control the water use efficiency in wet rice production. Its saturated hydraulic conductivity K_s, however, exhibits a large spatiotemporal variability as a consequence of a highly dynamic soil structure involving temporary shrinkage cracks. Water flow through the earthen bunds surrounding the fields further contributes to the uncertainty in water flux calculations. The objective of this study was to develop a simple deterministic model with stochastic elements (‘PADDY‐FLUX’) for depiction of deep percolation, and to assess the effect of different water management scenarios on percolation in two channel command areas. Darcy's law is used as the fundamental equation for water flow calculations with the ponding water depth h as a time‐dependent variable. Flux uncertainty is estimated by a Monte‐Carlo‐type implementation. K_s is treated as a random variable of a bimodal probability density function (PDF), which is the weighted sum of two Gaussian PDFs (accounting for a matrix and a preferential flow domain). The weighing factor α is a function of h, reflecting an increasing risk for preferential flow situations after desiccation and the development of shrinkage cracks. Under‐bund percolation is calculated using transfer functions. The results demonstrate that percolation losses increase in the following order: continuous soil saturation < continuous flooding (CF) < mid‐season drainage and intermittent irrigation (MD + II) < mid‐season drainage and continuous flooding. The bunds contribute up to 54 and 17% to total fluxes under CF and MD + II, respectively. Preferential water fluxes are responsible for the major part of water losses as soon as desiccation causes the formation of shrinkage cracks. As a conclusion, continuous soil saturation should be promoted as the least water‐intensive irrigation regime, while intermittent irrigation is recommended only in case that irreversible shrinkage cracks have already developed. Copyright © 2010 John Wiley & Sons, Ltd.     

Soil water dynamics in cropped and uncropped fields in northern Greece using a dual-permeability model

Abstract

This study focuses on the calibration and validation of a dual-permeability soil water flow model for simulating soil water dynamics during the growing period in an irrigated corn field and during the rainy winter period in an uncropped field in northern Greece. The 1D numerical transient dual-permeability model MACRO 5.0 was used to describe the soil water dynamics, the water balance and deep percolation considering both macropore (two-domain) flow and non-macropore (one-domain) flow. The simulated results were compared with measurements of total soil water content at different depths in the soils. The values of the statistical criteria RMSE, E and CRM were better when macroporosity flow was considered; the soil water content showed better redistribution in the soil profile. The limited irrigation of the corn field during the growing period and the irrigation rates did not create conditions for deep percolation of water. In the uncropped field (bare soil), the wet conditions and the high rainfall during the simulation period created conditions for significant deep percolation, whether macropore flow was included in the model or not. The two-domain approach significantly affects the actual evaporation and the deep percolation. The difference between these two approaches is in the amount of deep percolation and the flow path of drainage flow. In the two-domain approach, most deep percolation follows the macropore domain (79.8%). The errors due to macropore parameter uncertainty and to the difficulties of measuring the macropore water content and flow were estimated by a sensitivity analysis for the more important parameters of the model.

Editor Z.W. Kundzewicz

Citation Antonopoulos, V.Z., Georgiou, P.E., and Kolotouros, C.A., 2013. Soil water dynamics in cropped and uncropped fields in northern Greece using a dual-permeability model. Hydrological Sciences Journal, 58 (8), 1748–1759.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Do not only connect: a model of infiltration‐excess overland flow based on simulation

The paper focusses on connectivity in the context of infiltration‐excess overland flow and its integrated response as slope‐base overland flow hydrographs. Overland flow is simulated on a sloping surface with some minor topographic expression and spatially differing infiltration rates. In each cell of a 128 × 128 grid, water from upslope is combined with incident rainfall to generate local overland flow, which is stochastically routed downslope, partitioning the flow between downslope neighbours. Simulations show the evolution of connectivity during simple storms. As a first approximation, total storm runoff is similar everywhere, discharge increasing proportionally with drainage area. Moderate differences in plan topography appear to have only a second‐order impact on hydrograph form and runoff amount. Total storm response is expressed as total runoff, runoff coefficient or total volume infiltrated; each plotted against total storm rainfall, and allowing variations in average gradient, overland flow roughness, infiltration rate and storm duration. A one‐parameter algebraic expression is proposed that fits simulation results for total runoff, has appropriate asymptotic behaviour and responds rationally to the variables tested. Slope length is seen to influence connectivity, expressed as a scale distance that increases with storm magnitude and can be explicitly incorporated into the expression to indicate runoff response to simple events as a function of storm size, storm duration, slope length and gradient. The model has also been applied to a 10‐year rainfall record, using both hourly and daily time steps, and the implications explored for coarser scale models. Initial trails incorporating erosion continuously update topography and suggest that successive storms produce an initial increase in erosion as rilling develops, while runoff totals are only slightly modified. Other factors not yet considered include the dynamics of soil crusting and vegetation growth. Copyright © 2014 John Wiley & Sons, Ltd.     

Water percolation process studies in a Mediterranean karst area

To investigate processes of water percolation, the drip response of stalactites in a karstic cave below a 143 m² sprinkling plot was measured. The experiment was conducted in Mount Carmel, Israel, at the end of the dry season and intended to simulate a series of two high‐intensity storms on dry and wet soils. In addition to hydrometric measurements (soil moisture, surface runoff, stalactite dripping rates), two types of tracers (electrical conductivity and bromide) were used to study recharge processes, water origin and mixing inside a 28‐m vadose zone. Results suggested that slow, continuous percolation through the rock matrix is of minor importance and that percolating water follows a complicated pattern including vertical and horizontal flow directions. While bromide tracing allowed identification of quick direct flow paths at all drips with maximum flow velocities of 4·3 m/h, mixing analysis suggested that major water fractions were mobilized by piston flow, pushing out water stored in the unsaturated zone above the cave. Under dry preconditions, 80 mm of artificial rainfall applied in less than 7 h was not enough to initiate significant downward water percolation. Most water was required to fill uppermost soil and rock storages. Under wet preconditions during the second day sprinkling, higher water contents in soils and karst cavities facilitated piston flow effects and a more intense response of the cave drips. Results indicate that in Mediterranean karst regions, filling of the unsaturated zone, including soil and rock storages, is an important precondition for the onset of significant water percolation and recharge. This results in a higher seasonal threshold for water percolation than for the generation of surface runoff. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

A model to couple overland flow and infiltration into macroporous vadose zone

Most vegetated land surfaces contain macropores that may have a significant effect on the rate of infiltration of water under ponded conditions on the ground surface. Owing to the small-scale variations of the land topography (microtopography), only portions of the land area may get ponded during the process of overland flow. As the macropores transmit water at much higher rates than the primary soil matrix, higher macropore activation in ponded areas produces larger effective infiltration rates into the soil. Therefore, overland flow and infiltration into the macroporous vadose zone are interrelated. Representing the microtopographic variation of the land surface by a simple sine wave function, a method was developed to relate the ponding area to the average ponding depth which was determined by overland flow. A numerical model coupling overland flow and infiltration into the macroporous vadose zone was developed. Overland flow was simulated using the St. Venant equations with the inertia terms neglected. A single macropore model was used to simulate the infiltration into the macroporous vadose zone. The interaction between overland flow and the infiltration into the macroporous vadose zone was analyzed for a hypothetical watershed. The sensitivity analysis revealed that the interaction of macropore flow and overland flow is significant. For the conditions tested, the macropore flow and the overland flow were found to be more sensitive to the macroporosity and less sensitive to the microtopographic surface variation.