Evaluation of seasonal patterns of hydraulic redistribution in a humid subtropical area,East China

Antecedent soil moisture or soil moisture status has a great impact on hydrological processes. Hydraulic redistribution (HR), a widely observed phenomenon, is defined as water distributed (typically at night) from moist soil to drier soil via plant roots, which plays an important role in soil moisture replenishment. Knowledge on seasonal patterns of HR and on the relationship between HR and soil water use is not fully understood. We investigated temporal variations in HR and total daily water use (Δθ) at stands of camphor and peach by monitoring soil moisture content in a humid region in eastern China. HR at the three locations reached its maximum values in summer (0.68 mm d⁻¹ to 1.15 mm d⁻¹) at depths of 15 cm and 35 cm. Redistributed water replenished 41% of water depleted in the soil at a 5–45-cm depth. Interestingly, normalized HR (i.e., HR/Δθ) showed a constant pattern during the growing season implying it is independent of seasonal climate alterations. This also indicated that HR had a stable effect on the replenishment of daily water use. Positive linear relationships between HR and Δθ were found at three measuring locations (camphor: R² = .35, p < .01; peach1: R² = .57, p < .01; peach2: R² = .63, p < .01), suggesting a relatively stable inherent linkage between HR and Δθ. This study suggested that hydrological processes involving soil moisture status or antecedent soil moisture, needs to take the HR effect into account across timescales from intraday to seasonal.     

Effective hydraulic conductivity of nonstationary aquifers

Assuming that the ln hydraulic conductivity in an aquifer is mathematically approximated by a spatial deterministic surface, or trend, plus a stationary random noise, we treat the problem of finding what the effective hydraulic conductivity of that aquifer is. This problem is tackled by spectral methods applied to a type of diffusion equation of groundwater flow, together with suitable coordinate transformations. Analytical (exact) solutions in terms of elementary functions are presented for one- and three-dimensional finite and infinite domains. Stability criteria are obtained for the solutions, in terms of a critical parameter, that turns out to involve the product of correlation scale and trend gradient. For the case of finite and symmetrical domains, additional provisions to insure the stability of numerical calculations of effective hydraulic conductivity are provided. Effective hydraulic conductivity is an important property, with potential applications in the calibrations of groundwater and transport numerical models.     

Effective hydraulic conductivity of nonstationary aquifers

Assuming that the ln hydraulic conductivity in an aquifer is mathematically approximated by a spatial deterministic surface, or trend, plus a stationary random noise, we treat the problem of finding what the effective hydraulic conductivity of that aquifer is. This problem is tackled by spectral methods applied to a type of diffusion equation of groundwater flow, together with suitable coordinate transformations. Analytical (exact) solutions in terms of elementary functions are presented for one- and three-dimensional finite and infinite domains. Stability criteria are obtained for the solutions, in terms of a critical parameter, that turns out to involve the product of correlation scale and trend gradient. For the case of finite and symmetrical domains, additional provisions to insure the stability of numerical calculations of effective hydraulic conductivity are provided. Effective hydraulic conductivity is an important property, with potential applications in the calibrations of groundwater and transport numerical models.     

Identification of optimal soil hydraulic functions and parameters for predicting soil moisture

Abstract

The accuracy of six combined methods formed by three commonly-used soil hydraulic functions and two methods to determine soil hydraulic parameters based on a soil hydraulic parameter look-up table and soil pedotransfer functions was examined for simulating soil moisture. A novel data analysis and modelling approach was used that eliminated the effects of evapotranspiration so that specific sources of error among the six combined methods could be identified and quantified. By comparing simulated and observed soil moisture at six sites of the USDA Soil Climate Analysis Network, we identified the optimal soil hydraulic functions and parameters for predicting soil moisture. Through sensitivity tests, we also showed that adjusting only the soil saturated hydraulic conductivity, K_s , is insufficient for representing important effects of macropores on soil hydraulic conductivity. Our analysis illustrates that, in general, soil hydraulic conductivity is less sensitive to K_s than to the soil pore-size distribution parameter.

Editor D. Koutsoyiannis; Associate editor D. Hughes

Citation Pan, F., McKane, R.B. and Stieglitz, M., 2012. Identification of optimal soil hydraulic functions and parameters for predicting soil moisture. Hydrological Sciences Journal, 57 (4), 723–737.     

Year‐round estimation of soil moisture content using temporally variable soil hydraulic parameters

Soil moisture plays a key role in the hydrological cycle as it controls the flux of water between soil, vegetation, and atmosphere. This study is focused on a year‐round estimation of soil moisture in a forested mountain area using the bucket model approach. For this purpose, three different soil moisture models are utilised. The procedure is based on splitting the whole year into two complement periods (dormant and vegetation). Model parameters are allowed to vary between the two periods and also from year to year in the calibration procedure. Consequently, two sets of average model parameters corresponding to dormant and vegetation seasons are proposed. The process of splitting is strongly supported by the experimental data, and it enables us to variate saturated hydraulic conductivity and pore‐size characterisation. The use of the two different parameter sets significantly enhances the simulation of two (Teuling and Troch model and soil water balance model‐green–ampt [SWBM‐GA]) out of three models in the 6‐year period from 2009 to 2014. For these two models, the overall Nash‐Sutcliffe coefficient increased from 0.64 to 0.79 and from 0.55 to 0.80. The third model (the Laio approach) proved to be insensitive to parameter changes due to its insufficient drainage prediction. The variability of the warm and cold parameter sets between particular years is more pronounced in the warm periods. The cold periods exhibited approximately similar character during all 6 years.     

Measurement of soil moisture change in spatially heterogeneous weathered soils using a capacitance probe

A capacitance probe has been employed to measure the spatial and temporal change of soil moisture in a weathered soil profile with heterogeneous texture characteristics. Measurements were made to a depth of 2 m at intervals of 2 cm. The field work was carried out during a wetting-up period in a subtropical monsoonal environment in the Middle Hills of Nepal. Calibrations of the frequency readings against gravimetric and volumetric soil moisture are compared. It was found that the frequency reading correlates much better with the volumetric soil moisture. Although the soil is spatially heterogeneous in its particle size, porosity and bulk density, it was found that a single regression equation can be used to represent the relationship between the frequency reading of the capacitance probe and volumetric soil moisture. The results demonstrate that the capacitance probe enables the investigation of soil moisture variation temporally and spatially at high resolution and can measure a wide moisture range. © 1998 John Wiley & Sons, Ltd.     

Climate,vegetation, and soil controls on hydraulic redistribution in shallow tree roots

Hydraulic redistribution defined as the translocation of soil moisture by plant root systems in response to water potential gradients is a phenomenon widely documented in different climate, vegetation, and soil conditions. Past research has largely focused on hydraulic redistribution in deep tree roots with access to groundwater and/or winter rainfall, while the case of relatively shallow (i.e., ≈1–2 m deep) tree roots has remained poorly investigated. In fact, it is not clear how hydraulic redistribution in shallow root zones is affected by climate, vegetation, and soil properties. In this study, we developed a model to investigate the climate, vegetation, and soil controls on the net direction and magnitude of hydraulic redistribution in shallow tree root systems at the growing season to yearly timescale. We used the model to evaluate the effect of hydraulic redistribution on the water stress of trees and grasses. We found that hydraulic lift increases with decreasing rainfall frequency, depth of the rooting zone, root density in the deep soil and tree leaf area index; at the same time for a given rainfall frequency, hydraulic lift increases with increasing average rainstorm depth and soil hydraulic conductivity. We propose that water drainage into deeper soil layers can lead to the emergence of vertical water potential gradients sufficient to explain the occurrence of hydraulic lift in shallow tree roots without invoking the presence of a shallow water table or winter precipitation. We also found that hydraulic descent reduces the water stress of trees and hydraulic lift reduces the water stress of grass with important implications on tree–grass interactions.     

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.     

The effect of conservation practices in sloped croplands on soil hydraulic properties and root‐zone moisture dynamics

Rain‐induced erosion and short‐term drought are the two factors that limit the productivity of croplands in the red soil region of subtropical China. The objective of this study was to estimate the effects of conservation practices on hydraulic properties and root‐zone water dynamics of the soil. A 3‐year experiment was performed on a slope at Xianning. Four treatments were evaluated for their ability to reduce soil erosion and improve soil water conditions. Compared with no practices (CK) and living grass strips (GS), the application of polyacrylamide (PAM) significantly reduced soil crust formation during intense rainfall, whereas rice straw mulching (SM) completely abolished soil crust formation. The SM and PAM treatments improved soil water‐stable aggregates, with a redistribution of micro‐aggregates into macro‐aggregates. PAM and SM significantly increased the soil water‐holding capacity. These practices mitigated the degradation of the soil saturated hydraulic conductivity (Ks) during intense rainfalls. These methods increased soil water storage but with limited effects during heavy rainfalls in the wet period. In contrast, during the dry period, SM had the highest soil water storage, followed by PAM and CK. Grass strips had the lowest soil water storage because of the water uptake during the vigorous grass growth. A slight decline in the soil moisture resulted in a significant decrease in the unsaturated hydraulic conductivity (Ku) of the topsoil. Therefore, the hydraulic conductivity in the field is governed by soil moisture, and the remaining soil moisture is more important than improving soil properties to resist short‐term droughts. As a result, SM is the most effective management practice when compared with PAM and GS, although they all protect the soil hydraulic properties during wet periods. These results suggest that mulching is the best strategy for water management in erosion‐threatened and drought‐threatened red soils. Copyright © 2014 John Wiley & Sons, Ltd.     

Soil water content vertical profiles under natural conditions: matching of experiments and simulations by a conceptual model

The prediction of soil moisture content, θ, as a function of depth, z, and time, t, is of fundamental importance for applications in many hydrological processes. The main objective of this paper is to provide an approach to solve this problem at a local scale in soils with vegetation. The matching of soil moisture vertical profiles observed under natural conditions in grassy plots and their simulations by a conceptual model is presented. Experimental measurements were performed in a plot located in Central Italy, complete with hydrometeorological sensors specifically set up and equipped with a time domain reflectometry system providing the water content, θ_e(z, t). A conceptual model framework earlier proposed for two‐layered soil vertical profiles was modified and adopted for simulations. The changes concern the incorporation of evapotranspiration, the reduction of the original model for applications also to homogeneous soil vertical profiles, and a correction for the differences existing between assumed and observed initial moisture contents. In the model calibration, it was found that the effects of vegetation could be represented adequately by a fictitious soil vertical profile with a more permeable upper layer of saturated hydraulic conductivity, K_s, independent of time. Then, for the validation events, the model simulations in the stages of both infiltration and redistribution/evapotranspiration reproduced appropriately θ_e(z, t) with typical values of root mean square error in the range 0.0017–0.0657. Similar results were obtained by applying the modified two‐layered model for simulations of experimental data observed in three other plots located in Northern Italy and Germany. For all four vegetated sites, the two‐layer profile better matched the experimental data than the assumption of a homogeneous profile. Thus, the conceptual approach based on a two‐layered scheme for representing θ(z, t) in soils with vegetation appears to be appropriate for many hydrological applications. Copyright © 2013 John Wiley & Sons, Ltd.     

Impact of evaporation on field capacity during water drainage redistribution in a soil

To this day, field capacity (FC) is rarely defined in the context of soil properties, and the use of non-physical simplistic models is the common way to normalize water content at FC. In this study, the problem of water drainage redistribution in a soil column with and without the presence of evaporation (EV) was extensively studied. Analytical solutions for the Richards equation were established for the case of water drainage redistribution through a deeply wetted soil water column with and without EV at FC conditions. Water retention and depth evolution curves were plotted first, using different EV values of (2 $\frac{\mathbf{mm}}{\mathbf{day}}$, 5 $\frac{\mathbf{mm}}{\mathbf{day}}$ and 8 $\frac{\mathbf{mm}}{\mathbf{day}}$) and second, for different drainage redistribution durations of (1 day, 4 days and 6 days) where EV was set to zero for the case with no EV or to a fixed value of 5 $\frac{\mathbf{mm}}{\mathbf{day}}$ for the case with EV. The results suggest that EV plays a significant role in soil water drainage suggesting that, in the presence of EV, the FC drying front reaches much higher depths in the soil water profile than if EV is turned off. It was also concluded that FC reaches deeper depths faster the stronger EV is acting at the surface of a soil water column. Additionally, the results suggest that the texture of the soil receiving drainage controls the amount of water available for EV and as a result, EV was found to play a stronger role the smaller the hydraulic conductivity of the soil is.     

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.     

Relationships between hydraulic parameters in a small stream under varying flow and seasonal conditions

Twenty conservative tracer injections were carried out in the same reach of a small woodland stream in order to determine how variation in discharge and leaf accumulation affect stream hydraulic parameters. The injections were made at various discharge rates ranging from 2·6 to 40 l/s. Five of the injections were made during late autumn, when there were large accumulations of leaves in the stream. Estimates of hydraulic parameters were made by fitting a transient storage solute transport model to conservative tracer concentration profiles. Velocity increased almost linearly with increasing discharge, indicating a decline in the Darcy friction factor. Dispersion also increased with increasing discharge, especially for the lower flow injections. The relative size of the storage zone was small (∽0·1). There was no definable relationship between discharge and the relative storage zone size, but the rates of exchange between the storage zone and the main channel increased markedly with increasing discharge. The presence of large accumulations of leaves had a clear effect on the hydraulic characteristics of the stream, producing much higher friction factors, larger storage zone sizes and lower velocity than would have been predicted by discharge alone. Copyright © 1999 John Wiley & Sons, Ltd.     

Hydraulic redistribution and hydrological controls on aspen transpiration and establishment in peatlands following wildfire

In the sub‐humid Western Boreal Plains of Alberta, where evapotranspiration often exceeds precipitation, trembling aspen (Populus tremuloides Michx.) uplands often depend on adjacent peatlands for water supply through hydraulic redistribution. Wildfire is common in the Boreal Plains, so the resilience of the transfer of water from peatlands to uplands through roots immediately following wildfire may have implications for aspen succession. The objective of this research was to characterize post‐fire peatland‐upland hydraulic connectivity and assess controls on aspen transpiration (as a measure of stress and productivity) among landscape topographic positions. In May 2011, a wildfire affected 90,000 ha of north central Alberta, including the Utikuma Region Study Area (URSA). Portions of an URSA glacio‐fluval outwash lake catchment were burned, which included forests and a small peatland. Within 1 year after the fire, aspen were found to be growing in both the interior and margins of this peatland. Across recovering land units, transpiration varied along a topographic gradient of upland midslope (0.42 mm hr^?1) > upland hilltop (0.29 mm hr^?1) > margin (0.23 mm hr^?1) > peatland (0.10 mm hr^?1); similar trends were observed with leaf area and stem heights. Although volumetric water content was below field capacity, P. tremuloides were sustained through roots present, likely before fire, in peatland margins through hydraulic redistribution. Evidence for this was observed through the analysis of oxygen (δ¹⁸O) and hydrogen (δ²H) isotopes where upland xylem and peat core signatures were ?10.0‰ and ?117.8‰ and ?9.2‰ and ?114.0‰, respectively. This research highlights the potential importance of hydraulic redistribution to forest sustainability and recovery, in which the continued delivery of water may result in the encroachment of aspen into peatlands. As such, we suggest that through altering ecosystem services, peatland margins following fire may be at risk to aspen colonization during succession.     

On the dependence of the saturated hydraulic conductivity upon the effective porosity through a power law model at different scales

We study the scale dependence of the saturated hydraulic conductivity K_s through the effective porosity n_e by means of a newly developed power‐law model (PLM) which allows to use simultaneously measurements at different scales. The model is expressed as product between a single PLM (capturing the impact of the dominating scale) and a characteristic function κ^? accounting for the correction because of the other scale(s). The simple (closed form) expression of the κ^?‐function enables one to easily identify the scales which are relevant for K_s. The proposed model is then applied to a set of real data taken at the experimental site of Montalto Uffugo (Italy), and we show that in this case two (i.e. laboratory and field) scales appear to be the main ones. The implications toward an important application (solute transport) in Hydrology are finally discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

Temporal sampling and role of flux measurements for subsurface heterogeneous characterization in groundwater basins using hydraulic tomography

Accurate characterization of heterogeneity in groundwater basins is crucial to the sustainable management of groundwater resources. This study explores the temporal sampling issues and the role of flux measurements in the characterization of heterogeneity in groundwater basins using numerical experiments. The experiments involve a digital basin imitating the groundwater basin of the North China Plain (NCP), where the groundwater exploitation reduction program is ongoing. Using the experiments, we champion that the reduction program could collect groundwater level information induced by operational variations of existing pumping wells at different locations in the basin. Such a dataset could serve as a basin-scale hydraulic tomography (HT) to characterize the basin-scale heterogeneity cost-effectively. Both steady-state and transient-state inversion experiments demonstrate the advantage of HT surveys in characterizing basin-scale heterogeneity over conventional pumping tests at fixed well locations. Additionally, head data at the early, intermediate, and late time from well hydrographs should be selected for the HT analysis to maximize HT's power and save computational costs. When accurate geological zones are incorporated in prior information, flux measurements significantly improve parameter estimates based on conventional pumping tests. However, their effects are less noticeable for long-term HT surveys in such basin-scale aquifers without fissures or fractures. This basin-scale tomographic survey example serves a guide for field data collection and optimization of the analysis of future basin-scale HT.     

Establishing the role of pedogenesis in changing soil hydraulic properties

Chronosequences provide suitable sources for the study of changes in soil hydraulic behaviour as a result of long-term pedogenesis. For a podsol chronosequence in the Scottish Highlands, data are presented to indicate the changes that have occurred over 13000 years in the saturated hydraulic conductivity (K_sat) in each horizon. As the soil profile has evolved into a differentiated sequence of three horizons, the resulting hydrological changes can be both measured and quantified by relating K_sat to textural properties and bulk density. The results are significant for interpretation of changing runoff processes and slope stability.     

Obtaining basic simulation data for a heterogeneous field with stratified marine soils

Thinly stratified sedimentary deposits in a heterogeneous field were investigated to obtain basic physical data for the simulation of water flow. A procedure is described which translates a thinly stratified soil profile into a number of functional layers using functional hydrological properties. A functional layer is defined as a combination of one or more soil horizons and should (i) be recognizable during a soil survey using an auger and (ii) show significantly different functional hydrological properties when compared with another functional layer. This procedure gave three easily recognizable functional layers. Sets of hydrological characteristics of these three functional layers were obtained by physical measurements of the soil and by estimation, using textural data for classification into a standard Dutch series. The performance of several combinations of these sets was tested by comparing simulated and measured soil matric potentials for seven plots during one year. The best simulation results were obtained if measured soil hydraulic characteristics were used for relatively homogeneous functional layers and if the soil hydraulic characteristics were estimated at each location for the most heterogeneous layer.     

On estimating the hydraulic properties of soil,Part 2. A new empirical equation for estimating hydraulic conductivity for sands

A new empirical equation to estimate hydraulic conductivity is proposed, based on a large set of measured data for hydraulic properties of soil. The equation is simpler and more accurate than the series-parallel model. Under conditions of insufficient data, the new equation provides a good estimation of hydraulic conductivity for sands. For the same class of soils, another empirical equation is proposed to estimate the power N in the Averjanov-Irmay function.