Estimation of soil hydraulic properties and their uncertainty through the Beerkan infiltration experiment

The Beerkan method based on in situ single‐ring water infiltration experiments along with the relevant specific Beerkan estimation of soil transfer parameters (BEST) algorithm is attractive for simple soil hydraulic characterization. However, the BEST algorithm may lead to erroneous or null values for the saturated hydraulic conductivity and sorptivity especially when there are only few infiltration data points under the transient flow state, either for sandy soil or soils in wet conditions. This study developed an alternative algorithm for analysis of the Beerkan infiltration experiment referred to as BEST‐generalized likelihood uncertainty estimation (GLUE). The proposed method estimates the scale parameters of van Genuchten water retention and Brooks–Corey hydraulic conductivity functions through the GLUE methodology. The GLUE method is a Bayesian Monte Carlo parameter estimation technique that makes use of a likelihood function to measure the goodness‐of‐fit between modelled and observed data. The results showed that using a combination of three different likelihood measurements based on observed transient flow, steady‐state flow and experimental steady‐state infiltration rate made the BEST‐GLUE procedure capable of performing an efficient inverse analysis of Beerkan infiltration experiments. Therefore, it is more applicable for a wider range of soils with contrasting texture, structure, and initial and saturated water content. Copyright © 2015 John Wiley & Sons, Ltd.     

Determining hydraulic properties of a loam soil by alternative infiltrometer techniques

Testing infiltrometer techniques to determine soil hydraulic properties is necessary for specific soils. For a loam soil, the water retention and hydraulic conductivity predicted by the BEST (Beerkan Estimation of Soil Transfer parameters) procedure of soil hydraulic characterization was compared with data collected by more standard laboratory and field techniques. Six infiltrometer techniques were also compared in terms of saturated soil hydraulic conductivity, K_s. BEST yielded water retention values statistically similar to those obtained in the laboratory and K_s values practically coinciding with those determined in the field with the pressure infiltrometer (PI). The unsaturated soil hydraulic conductivity measured with the tension infiltrometer (TI) was reproduced satisfactorily by BEST only close to saturation. BEST, the PI, one‐potential experiments with both the TI and the mini disk infiltrometer (MDI), the simplified falling head (SFH) technique and the bottomless bucket (BB) method yielded statistically similar estimates of K_s, differing at the most by a factor of three. Smaller values were obtained with longer and more soil‐disturbing infiltration runs. Any of the tested infiltration techniques appears usable to obtain the order of magnitude of K_s at the field site, but the BEST, BB and PI data appear more appropriate to characterize the soil at some stage during a rainfall event. Additional investigations on both similar and different soils would allow development of more general procedures to apply infiltrometer techniques for soil hydraulic characterization. Copyright © 2015 John Wiley & Sons, Ltd.     

Synergistic effects of water repellency and macropore flow on the hydraulic conductivity of a burned forest soil,south‐east Australia

Research shows that water repellency is a key hydraulic property that results in reduced infiltration rates in burned soils. However, more work is required in order to link the hydrological behaviour of water repellent soils to observed runoff responses at the plot and hillslope scale. This study used 5 M ethanol and water in disc infiltrometers to quantify the role of macropore flow and water repellency on spatial and temporal infiltration patterns in a burned soil at plot (<10 m²) scale in a wet eucalypt forest in south‐east Australia. In the first summer and winter after wildfire, an average of 70% and 60%, respectively, of the plot area was water repellent and did not contribute to infiltration. Macropores (r > 0·5 mm), comprising just 5·5% of the soil volume, contributed to 70% and 95%, respectively, of the field‐saturated and ponded hydraulic conductivity (K_p). Because flow occurred almost entirely via macropores in non‐repellent areas, this meant that less than 2·5% of the soil surface effectively contributed to infiltration. The hydraulic conductivity increased by a factor of up to 2·5 as the hydraulic head increased from 0 to 5 mm. Due to the synergistic effect of macropore flow and water repellency, the coefficient of variation (CV) in K_p was three times higher in the water‐repellent soil (CV = 175%) than under the simulated non‐repellent conditions (CV = 66%). The high spatial variability in K_p would act to reduce the effective infiltration rate during runoff generation at plot scale. Ponding, which tend to increase with increasing scale, activates flow through macropores and would raise the effective infiltration rates at larger scales. Field experiments designed to provide representative measurements of infiltration after fire in these systems must therefore consider both the inherent variability in hydraulic conductivity and the variability in infiltration caused by interactions between surface runoff and hydraulic conductivity. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydraulic characterization and hydrological behaviour of a pilot permeable pavement in an urban centre,Brazil

Considering all the alterations on hydrology and water quality that urbanization process brings, permeable pavement (PP) is an alternative to traditional impermeable asphalt and concrete pavement. The goal of the PP and other low impact development devices is to increase infiltration and reduce peak runoff flows. These structures are barely used in Brazil aiming stormwater management, one of the big hydrological issues in cities throughout the country, with increasing urbanization rates. The main objective of this paper is the hydraulic characterization of a PP and the assessment of its hydrological efficiency from the point of view of the infiltration process. The study focuses on a pilot area in a parking lot in an urban area (Recife, Brazil). Soil elements filling the voids between concrete elements were sampled (particle size density, water contents) and tested with water infiltration experiments at several points of the 3 m × 1.5 m surface pilot area. Beerkan Estimation of Soil Transfer parameters algorithm was applied to the infiltration experiment data to obtain the hydraulic characteristics of the soil composing the PP surface layer, the concrete grid pavers (with internal voids filled with natural soil) permeability being neglected. Results show that the soil hydraulic characteristics vary spatially within the pilot area and that the soil samples have different hydraulic behaviours. The hydraulic characteristics derived from Beerkan Estimation of Soil Transfer parameters analysis were implemented into Hydrus code to simulate runoff, infiltration and water balance over a year. The numerical simulation showed the good potential of the PP for rainfall–runoff management, which demonstrates that PP can be used to retrofit existing parking infrastructure and to promote hydrological behaviour close to natural soils. Copyright © 2016 John Wiley & Sons, Ltd.     

Solute transport scales in an unsaturated stony soil

Solute transport parameters are known to be scale-dependent due mainly to the increasing scale of heterogeneities with transport distance and with the lateral extent of the transport field examined. Based on a transect solute transport experiment, in this paper we studied this scale dependence by distinguishing three different scales with different homogeneity degrees of the porous medium: the observation scale, transport scale and transect scale. The main objective was to extend the approach proposed by van Wesenbeeck and Kachanoski to evaluating the role of textural heterogeneities on the transition from the observation scale to the transport scale. The approach is based on the scale dependence of transport moments estimated from solute concentrations distributions. In our study, these moments were calculated starting from time normalized resident concentrations measured by time domain reflectometry (TDR) probes at three depths in 37 soil sites 1 m apart along a transect during a steady state transport experiment. The Generalized Transfer Function (GTF) was used to describe the evolution of apparent solute spreading along the soil profile at each observation site by analyzing the propagation of the moments of the concentration distributions. Spectral analysis was used to quantify the relationship between the solid phase heterogeneities (namely, texture and stones) and the scale dependence of the solute transport parameters. Coupling the two approaches allowed us to identify two different transport scales (around 4-5 m and 20 m, respectively) mainly induced by the spatial pattern of soil textural properties. The analysis showed that the larger transport scale is mainly determined by the skeleton pattern of variability. Our analysis showed that the organization in hierarchical levels of soil variability may have major effects on the differences between solute transport behavior at transport scale and transect scale, as the transect scale parameters will include information from different scales of heterogeneities.     

Determining short-term changes in the hydraulic properties of a sandy-loam soil by a three-run infiltration experiment

ABSTRACT

Soil structure-dependent parameters can vary rapidly as a consequence of perturbing events such as intense rainfall. Investigating their short-term changes is therefore essential to understand the general behaviour of a porous medium. The aim of this study is to gain insight into the effects of wetting, perturbation and recovery processes through different sequences of Beerkan infiltration experiments performed on a sandy-loam soil. Two different three-run infiltration experiments (LHL and LLL) were carried out by pouring water at low (L, non-perturbing) and high (H, perturbing) heights above the soil surface and at short time intervals (hours, days). The results demonstrate that the proposed method allows one to capture short-term variations in soil structure-dependent parameters. The developed methodology is expected to simplify the parameterization of hydrological models with temporally variable soil hydraulic properties.     

Variability in surface moisture content along a hillslope transect: Rattlesnake Hill, Texas

Surface soil moisture content exhibits a high degree of spatial and temporal variability. The purpose of this study was (a) to characterize variations in moisture content in the 0–5 cm surface soil layer along a hillslope transect by means of intensive sampling in both space and time; and (b) to make inferences regarding the environmental factors that influence this variability. Over a period of seven months, soil moisture content was measured (gravimetric method) on a near-daily basis at 10 m intervals along a 200 m downslope transect at the Rattlesnake Hill field site in Austin, Texas. Results indicate that significant variability in soil moisture content exists along the length of the transect; that variability decreases with decreasing transect-mean moisture content as the hillslope dries down following rain events; and that the dominant influences on moisture content variability are dependent upon the moisture conditions on the hillslope. While topographic and soil attributes operate jointly to redistribute soil water following storm events, under wet conditions, variability in surface moisture content is most strongly influenced by porosity and hydraulic conductivity, and under dry conditions, correlations are strongest to relative elevation, aspect and clay content. Consequently, the dominant influence on soil moisture variability gradually changes from soil heterogeneity to joint control by topographic and soil properties as the transect dries following significant rain events.     

Functional evaluation of a simplified scaling method for assessing the spatial variability of soil hydraulic properties at the hillslope scale

Abstract

Mapping soil hydraulic parameters with traditional scaling methods that use laboratory-determined hydraulic characteristics (the LAB method) is not always feasible as it involves expensive, time-consuming and sophisticated measurements on soil samples collected in several locations of the study area. An alternative scaling method (the AP method) has been recently proposed to indirectly retrieve the soil hydraulic properties following the Arya-Paris physico-empirical pedotransfer function, which makes use of particle-size distribution and bulk density values. In this synthetic study we verify the performance of the AP method from a functional perspective, by evaluating the differences in the simulated soil water budget through a Monte Carlo approach. Notwithstanding that the AP method can provide soil hydraulic property patterns with faster experimental procedures and minor costs, we observe significant bias in the predicted spatially-averaged soil water budget due to a poor parametric calibration of the AP method and an imprecise identification of the spatial correlation structure of the AP-estimated scaling factors.

Citation Nasta, P., Romano, N., and Chirico, G.B., 2013. Functional evaluation of a simplified scaling method for assessing the spatial variability of soil hydraulic properties at the hillslope scale. Hydrological Sciences Journal, 58 (5), 1059–1071.     

Identifying effects of local and nonlocal factors of soil water storage using cyclical correlation analysis

There are various factors governing the spatial and temporal variability of soil water storage including soil properties, topography and vegetation. Some factors act locally, whereas others act nonlocally, which means that a factor measured at one location has effect on soil water storage at another location. The objective of this study was to examine the effects of local and nonlocal controls of soil water storage in a hummocky landscape using cyclical correlation analysis. Soil water storage, soil properties and terrain indices were measured along a 128‐point transect of 576 m long from the semiarid, hummocky, prairie pothole region of North America. There are large coefficients of determination (r²) between soil water storage and sand content (r² = 0.32–0.53), organic carbon content (r² = 0.22–0.56), depth to carbonate layer (r² = 0.13–0.63), wetness index (r² = 0.25–0.45) and other variables at the measurement scale at different times, indicating strong local effects from these variables. The correlation coefficients were also calculated by physically shifting the spatial series of soil water storage with respect to that of controlling factors. The shifting improves the correlation between the spatial series, and the length of shifting indicated the difference in the response of soil water to its controlling factors. For example, the value of r² increased more than eightfold (r² = 0.47–0.64) after shifting the spatial series of soil water storage by 54 m, almost equal to the average length of existing slope, compared with the very weak correlation (r² = 0.02–0.08) at the measurement scale. This indicated the nonlocal effect from the relative elevation. The identification of nonlocal effects from factors improves the prediction of soil water storage. Copyright © 2011 John Wiley & Sons, Ltd.     

LISEM: A SINGLE-EVENT,PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. II: SENSITIVITY ANALYSIS,VALIDATION AND APPLICATION

A new hydrological and soil erosion model has been developed and tested: LISEM, the Limburg soil erosion model. The model uses physically based equations to describe interception, infiltration and soil water transport, storage in surface depressions, splash and flow detachment, transport capacity and overland and channel flow. From the validation results it is clear that, although the model has several advantages over other models, the results of LISEM 1.0 are far from perfect. Based on the sensitivity analysis and field observations, the main reasons for these differences seems to be the spatial and temporal variability of the soil hydraulic conductivity and the initial pressure head at the basin scale. Another reason for the differences between measured and simulated results is our lack or understanding of the theory of hydrological and soil erosion processes.     

Impact of land use on the hydraulic properties of the topsoil in a small French catchment

The hydraulic properties of the topsoil control the partition of rainfall into infiltration and runoff at the soil surface. They must be characterized for distributed hydrological modelling. This study presents the results of a field campaign documenting topsoil hydraulic properties in a small French suburban catchment (7 km²) located near Lyon, France. Two types of infiltration tests were performed: single ring infiltration tests under positive head and tension‐disk infiltration using a mini‐disk. Both categories were processed using the BEST—Beerkan Estimation of Soil Transfer parameters—method to derive parameters describing the retention and hydraulic conductivity curves. Dry bulk density and particle size data were also sampled. Almost all the topsoils were found to belong to the sandy loam soil class. No significant differences in hydraulic properties were found in terms of pedologic units, but the results showed a high impact of land use on these properties. The lowest dry bulk density values were obtained in forested soils with the highest organic matter content. Permanent pasture soils showed intermediate values, whereas the highest values were encountered in cultivated lands. For saturated hydraulic conductivity, the highest values were found in broad‐leaved forests and small woods. The complementary use of tension‐disk and positive head infiltration tests highlighted a sharp increase of hydraulic conductivity between near saturation and saturated conditions, attributed to macroporosity effect. The ratio of median saturated hydraulic conductivity to median hydraulic conductivity at a pressure of − 20 mm of water was about 50. The study suggests that soil texture, such as used in most pedo‐transfer functions, might not be sufficient to properly map the variability of soil hydraulic properties. Land use information should be considered in the parameterizations of topsoil within hydrological models to better represent in situ conditions, as illustrated in the paper. Copyright © 2010 John Wiley & Sons, Ltd.     

Spatial heterogeneity and temporal dynamics of soil water tension in a mature Norway spruce stand

In ecosystem research great effort is made in measuring soil water tension, because this is a critical calibration variable for modelling soil water fluxes. In this paper the spatial heterogeneity and temporal dynamics of soil tensions and their consequences for the determination of water fluxes are investigated. Studies were carried out at a Norway spruce stand in the Fichtelgebirge (NE Bavaria). Standard tensiometers were installed at three soil depths (20 each) on the whole experimental plot, as well as 45 microtensiometers as a dense grid in a small soil pit. Microtensiometry at the centimetre scale showed that, depending on rain intensity and initial soil water tension, even a soil without discernible macrostructure may show preferential water infiltration. At the stand scale the variability of soil hydraulic properties and tree root distribution causes substantial heterogeneity of soil water tension, as observed by standard tensiometers. A functional relationship between increasing spatial heterogeneity of tensiometer readings and increasing soil water tension was found, which was particularly pronounced after longer dry periods. Also at low soil water tension, where spatial heterogeneity was low, the calculation of water fluxes from tensiometer values was critical, owing to the fact that small differences in measuring soil water tension resulted in big differences in calculated water fluxes. At high soil water tension in summer the spatial heterogeneity of tensiometer readings was extremely high. At our experimental site, since 30% of the total rain in summer falls in events having a precipitation rate greater than 5 mm h⁻¹, preferential water and solute flow was an important phenomenon. We conclude that the validation of calculated water fluxes using measured soil water tension at the stand scale is not an appropriate tool, because of measurement difficulties, considerable spatial heterogeneity, especially in dry periods, and the great variability of soil hydraulic properties. © 1998 John Wiley & Sons, Ltd.     

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.     

Spatial scale effect on the upper soil effective parameters of a distributed hydrological model

Nonlinear dynamics and spatial variability in hydrological systems make the formulation of scaling theories difficult. Therefore, the development of knowledge related to scale effects, scaling techniques, parameterization and linkages of parameters across scales is highly relevant. The main purpose of this work is to analyse the spatial effect of the static storage capacity parameter H_u and the saturated hydraulic conductivity parameter k_s from microscale (sub‐grid level) to mesoscale (grid level) and its implication to the definition of an optimum cell size. These two parameters describe the upper soil water characteristics in the infiltration process conceptualization of the TETIS hydrological model. At microscale, the spatial heterogeneity of H_u and k_s was obtained generating random parameter fields through probability distribution functions and a spatial dependence model with pre‐established correlation lengths. The effective parameters at mesoscale were calculated by solving the inverse problem for each parameter field. Results indicate that the adopted inverse formulation allows transferring the nonlinearity of the system from microscale to the mesoscale via non‐stationary effective parameters. Their values at each cell and time step are in the range of zero to the mean value of the parameter at microscale. The stochastic simulations showed that the variance of the estimated effective parameters decreases when the ratio between mesoscale cell size and correlation length at microscale increases. For a ratio greater than 1, we found cell sizes having the characteristics of a representative elementary area (REA); in such case, the microscale variability pattern did not affect the system response at mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatial patterns and temporal stability of soil moisture across a range of scales in a semi‐arid environment

The antecedent soil moisture status of a catchment is an important factor in hydrological modelling. Traditional Hortonian infiltration models assume that the initial moisture content is constant across the whole catchment, despite the fact that even in small catchments antecedent soil moisture exhibits tremendous spatial heterogeneity. Spatial patterns of soil water distribution across three transects (two in a burnt area and one in an unburnt area) in a semi‐arid area were studied. At the transect scale, when the factors affecting soil moisture were limited to topographical position or local topography, spatial patterns showed time stability, but when other factors, such as vegetation, were taken into account, the spatial patterns became time unstable. At the point scale, and in the same areas, topographical position was the main factor controlling time stability. Scale dependence of time stability was studied and local topography and vegetation presence were observed to play an important role for the correlation between consecutive measures depending on the scale. Copyright © 2000 John Wiley & Sons, Ltd.     

Comparing Beerkan infiltration tests with rainfall simulation experiments for hydraulic characterization of a sandy‐loam soil

Saturated soil hydraulic conductivity, K _s , data collected by ponding infiltrometer methods and usual experimental procedures could be unusable for interpreting field hydrological processes and particularly rainfall infiltration. The K _s values determined by an infiltrometer experiment carried out by applying water at a relatively large distance from the soil surface could however be more appropriate to explain surface runoff generation phenomena during intense rainfall events. In this study, a link between rainfall simulation and ponding infiltrometer experiments was established for a sandy‐loam soil. The height of water pouring for the infiltrometer run was chosen, establishing a similarity between the gravitational potential energy of the applied water, E _p , and the rainfall kinetic energy, E _k . To test the soundness of this procedure, the soil was sampled with the Beerkan estimation of soil transfer parameters procedure of soil hydraulic characterization and two heights of water pouring (0.03 m, i.e., usual procedure, and 0.34 m, yielding E _p  = E _k ). Then, a comparison between experimental steady‐state infiltration rates, i _sR , measured with rainfall simulation experiments determining runoff production and K _s values for the two water pouring heights was carried out in order to discriminate between theoretically possible (i _sR  ≥ K _s ) and impossible (i _sR  < K _s ) situations. Physically possible K _s values were only obtained by applying water at a relatively large distance from the soil surface, because i _sR was equal to 20.0 mm h^?1 and K _s values were 146.2–163.9 and 15.2–18.7 mm h^?1 for a height of water pouring of 0.03 and 0.34 m, respectively. This result suggested the consistency between Beerkan runs with a high height of water pouring and rainfall simulator experiments. Soil compaction and mechanical aggregate breakdown were the most plausible physical mechanisms determining reduction of K _s with height. This study demonstrated that the height from which water is poured onto the soil surface is a key parameter in infiltrometer experiments and can be adapted to mimic the effect of high intensity rain on soil hydraulic properties.     

The scaling characteristics of soil parameters: From plot scale heterogeneity to subgrid parameterization

The variation in soil texture, surface moisture or vertical soil moisture gradient in larger scale atmospheric models may lead to significant variations in simulated surface fluxes of water and heat. The parameterization of soil moisture fluxes at spatial scales compatible with the grid size of distributed hydrological models and mesoscale atmospheric models ( 100 km²) faces principal problems which relate to the underlying microscopic or field scale heterogeneity in soil characteristics.

The most widely used parameterization in soil hydrology, the Darcy-Richards (DR) equation, is gaining increasing importance in mesoscale and climate modelling. This is mainly due to the need to introduce plant-interactive soil water depletion and stomatal conductance parameterizations and to improve the calculation of deep percolation and runoff. Covering a grid of several hundreds of square kilometres, the DR parameterization in soil-vegetation-atmosphere-transfer schemes (SVATs) is assumed to be scale-invariant. The parameters describing the non-linear, area-average soil hydraulic functions in this scale-invariant DR-equation should be treated as calibration-parameters, which do not necessarily have a physical meaning. The saturated hydraulic conductivity is one of the soil parameters to which the models show very high sensitivity. It is shown that saturated hydraulic conductivity can be scaled in both vertical and horizontal directions for large flow domains.

In this paper, a distinction is made between effective and aggregated soil parameters. Effective parameters are defined as area-average values or distributions over a domain with a single, distinct textural soil type. They can be obtained by scaling or inverse modelling. Aggregated soil parameters represent grid-domains with several textural soil types. In soil science dimensional methods have been developed to scale up soil hydraulic characteristics. With some specific assumptions, these techniques can be extrapolated from classical field-scale problems in soil heterogeneity to larger domains, compatible with the grid-size of large scale models. Particularly promising is the estimation of effective soil hydraulic parameters from area averaging measurements through inverse modelling of the unsaturated flow.

Techniques to scale and aggregate the soil characteristics presented in this paper qualify for direct or indirect use in large scale meteorological models. One of the interesting results is the effective behaviour of the reference curve, which can be obtained from similar media scaling. If the conclusions of this paper survive further studies, a relatively simple method will become available to parameterize soil variability at large scales. The inverse technique is found to provide effective soil parameters which perform well in predicting both the area-average evaporation and the area-average soil moisture fluxes, such as subsurface runoff. This is not the case for aggregated soil parameters. Obtained from regression relationships between soil textural composition and hydraulic characteristics, these aggregated parameters predict evaporation fluxes well, but fail to predict water balance terms such as percolation and runoff. This is a serious drawback which could eventually hamper the improvement of the representation of the hydrological cycle in mesoscale atmospheric models and in GCMs.     

Interactions and connectivity between runoff generation processes of different spatial scales

Monitoring runoff generation processes in the field is a prerequisite for developing conceptual hydrological models and theories. At the same time, our perception of hydrological processes strongly depends on the spatial and temporal scale of observation. Therefore, the aim of this study is to investigate interactions between runoff generation processes of different spatial scales (plot scale, hillslope scale, and headwater scale). Different runoff generation processes of three hillslopes with similar topography, geology and soil properties, but differences in vegetation cover (grassland, coniferous forest, and mixed forest) within a small v‐shaped headwater were measured: water table dynamics in wells with high spatial and temporal resolution, subsurface flow (SSF) of three 10 m wide trenches at the bottom of the hillslopes subdivided into two trench sections each, overland flow at the plot scale, and catchment runoff. Bachmair et al. ( 2012 ) found a high spatial variability of water table dynamics at the plot scale. In this study, we investigate the representativity of SSF observations at the plot scale versus the hillslope scale and vice versa, and the linkage between hillslope dynamics (SSF and overland flow) and streamflow. Distinct differences in total SSF within each 10 m wide trench confirm the high spatial variability of the water table dynamics. The representativity of plot scale observations for hillslope scale SSF strongly depends on whether or not wells capture spatially variable flowpaths. At the grassland hillslope, subsurface flowpaths are not captured by our relatively densely spaced wells (3 m), despite a similar trench flow response to the coniferous forest hillslope. Regarding the linkage between hillslope dynamics and catchment runoff, we found an intermediate to high correlation between streamflow and hillslope hydrological dynamics (trench flow and overland flow), which highlights the importance of hillslope processes in this small watershed. Although the total contribution of SSF to total event catchment runoff is rather small, the contribution during peak flow is moderate to substantial. Additionally, there is process synchronicity between spatially discontiguous measurement points across scales, potentially indicating subsurface flowpath connectivity. Our findings stress the need for (i) a combination of observations at different spatial scales, and (ii) a consideration of the high spatial variability of SSF at the plot and hillslope scale when designing monitoring networks and assessing hydrological connectivity. Copyright © 2013 John Wiley & Sons, Ltd.     

Multi‐scale assimilation of root zone soil water predictions

When hydrology model parameters are determined, a traditional data assimilation method (such as Kalman filter) and a hydrology model can estimate the root zone soil water with uncertain state variables (such as initial soil water content). The simulated result can be quite good. However, when a key soil hydraulic property, such as the saturated hydraulic conductivity, is overestimated or underestimated, the traditional soil water assimilation process will produce a persistent bias in its predictions. In this paper, we present and demonstrate a new multi‐scale assimilation method by combining the direct insertion assimilation method, particle swarm optimisation (PSO) algorithm and Richards equation. We study the possibility of estimating root zone soil water with a multi‐scale assimilation method by using observed in situ data from the Wudaogou experiment station, Huaihe River Basin, China. The results indicate there is a persistent bias between simulated and observed values when the direct insertion assimilation surface soil water content is used to estimate root zone soil water contents. Using a multi‐scale assimilation method (PSO algorithm and direct insertion assimilation) and an assumed bottom boundary condition, the results show some obvious improvement, but the root mean square error is still relatively large. When the bottom boundary condition is similar to the actual situation, the multi‐scale assimilation method can well represent the root zone soil water content. The results indicate that the method is useful in estimating root zone soil water when available soil water data are limited to the surface layer and the initial soil water content even when the soil hydraulic conductivities are uncertain. Copyright © 2011 John Wiley & Sons, Ltd.     

Monitoring spatial variability of forest interception

Field instrumentation was designed and installed to quantify the influence of forest interception on the spatial and temporal distribution of water flux onto and into the forest soil at the plot scale. An application is presented which demonstrates that the instrumentation has the required resolution to monitor the spatial variability and dynamics of the flux processes. The observations show that spatial variability of interception may play an important role, not only in small scale soil moisture heterogeneity, but also in the hydrological response of a forested catchment at the hillslope scale. They also highlight the need of gathering more field information on the effects of vegetation on the spatial variability of soil surface water input.