An analysis of infiltration with moisture content distribution in a two‐dimensional discretized water content domain

On the basis of unsaturated Darcy's law, the Talbot–Ogden method provides a fast unconditional mass conservative algorithm to simulate groundwater infiltration in various unsaturated soil textures. Unlike advanced reservoir modelling methods that compute unsaturated flow in space, it only discretizes the moisture content domain into a suitable number of bins so that the vertical water movement is estimated piecewise in each bin. The dimensionality of the moisture content domain is extended from one dimensional to two dimensional in this study, which allows us to distinguish pore shapes within the same moisture content range. The vertical movement of water in the extended model imitates the infiltration phase in the Talbot–Ogden method. However, the difference in this extension is the directional redistribution, which represents the horizontal inter‐bin flow and causes the water content distribution to have an effect on infiltration. Using this extension, we mathematically analyse the general relationship between infiltration and the moisture content distribution associated with wetting front depths in different bins. We show that a more negatively skewed moisture content distribution can produce a longer ponding time, whereas a higher overall flux cannot be guaranteed in this situation. It is proven on the basis of the water content probability distribution independent of soil textures. To illustrate this analysis, we also present numerical examples for both fine and coarse soil textures. Copyright © 2014 John Wiley & Sons, Ltd.     

Water condition control of in situ soil water repellency: an observational study from a hillslope in a Japanese humid‐temperate forest

To evaluate the relationship between the moisture conditions and the water repellency of soil surfaces in situ, we periodically conducted field surveys for more than a year in a humid‐temperate forest in Japan. Measurements were made in four plots with varying soil physicochemical properties and under different topographic conditions across a hillslope. Each plot contained permanent quadrats with measurement points in a grid pattern. At each point, we measured the volumetric water content at 0‐ to 5‐cm depths and the water repellency at soil surfaces approximately twice a month. The repeated measurements enabled us to estimate the critical water content (CWC) below which soils repelled water at each point. We defined the representative CWC (RCWC) of a plot as the median of all CWCs in a plot and estimated the representative critical water potential (RCWP) on the basis of the RCWC using the water retention curve. The RCWC values differed among plots, but the corresponding RCWP values were similar (pF = 3.5–3.9). The relationship of the areal fraction showing water repellency against soil water potentials was similar across plots, but the relationship differed among plots against the soil moisture content. These results suggest that soil water potential is more indicative of the spatial occurrence of water repellency than moisture content on a hillslope where soil physicochemical properties vary. Plots located on ridge crests frequently exhibited lower water potentials and showed a higher areal fraction of water repellency, implying a greater chance of generating surface runoff by rainfall events. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimation of tree water stress from stem and soil water monitoring with time‐domain reflectometry in two small forested basins in Spain

Soil‐tree water relationships were studied using time domain reflectometry (TDR) in two small forested basins in Spain. The stem water content of two Mediterranean Quercus species (Quercus pyrenaica and Quercus rotundifolia) was measured using previously constructed species‐specific equations. To monitor soil moisture, a TDR station network was used in both cases. Sixteen Q. pyrenaica and six Q. rotundifolia individuals were selected to install two TDR probes in their trunks (at 20 and 120 cm above the ground) to monitor stem water content. Stem and soil water contents were measured fortnightly. The stem water content of both species showed a similar temporal trend for the period studied. A spring maximum (0·654 cm³ cm^?3 for Q. pyrenaica and 0·568 cm³ cm^?3 for Q. rotundifolia) was found to be associated with high transpiration and no soil moisture deficit, and a late‐summer minimum (0·520 cm³ cm^?3 for Q. pyrenaica and 0·426 cm³ cm^?3 for Q. rotundifolia) was associated with the end of the dry season. This drop in stem water content occurs when the available water in the soil decreases. This seasonal difference presumably reflects water withdrawn from stem storage to support the transpirational demands of the tree. Since plant water stress results in reduced stem water content and since this drop can be measured by TDR, it may be concluded that this technology offers a suitable tool for detecting plant water stress. Copyright © 2007 John Wiley & Sons, Ltd.     

First approximations of soil moisture retention curves using the filter‐paper method

Relationships between gravimetric soil moisture content (w) and matric potential (ϕ), and between volumetric soil moisture content (θ_v) and pressure head (h) were approximated for the unsaturated zone on Long Island, New York. Soil samples were collected from two sites using a hand auger. The soil moisture content was determined using the filter‐paper (w_f) and gravimetric (w) methods, respectively. The w_f was then used in an empirical equation to estimate ϕ_m. Each set of ϕ_m and w was combined with a straight‐line empirical model to obtain a w(ϕ_m) relationship. Soil ϕ_m was converted to h, and w to the volumetric moisture content θ_v, in order to produce a θ_v(h) curve. Graphical and statistical comparison showed that the resulting θ_v(h) curves fell within one order of magnitude of similar curves generated by a more sophisticated non‐linear model developed previously. The simplicity and low cost of the filter‐paper approach described in this study recommends it for preliminary studies of hydraulic properties in the unsaturated zone. Copyright © 2000 John Wiley & Sons, Ltd.     

Hillslope‐scale probabilistic characterization of soil moisture dynamics and average water balance

Probabilistic water balance modelling provides a useful framework for investigating the interactions between soil, vegetation, and the atmosphere. It has been used to estimate temporal soil moisture dynamics and ecohydrological responses at a point. This study combines a nonlinear rainfall–runoff theory with probabilistic water balance model to represent varied source area runoff as a function of rainfall depth and a runoff coefficient at hillslope scale. Analytical solutions of the soil‐moisture probability density function and average water balance model are then developed. Based on a sensitivity analysis of soil moisture dynamics, we show that when varied source area runoff is incorporated, mean soil moisture is always lower and total runoff higher, compared with the original probabilistic water balance model. The increased runoff from the inclusion of varied source area runoff is mainly because of a reduction in leakage when the index of dryness is less than one and evapotranspiration when the index of dryness is greater than one. Inclusion of varied source area runoff in the model means that the actual evapotranspiration is limited by less available water (i.e. water limit), which is stricter than Budyko’s and Milly’s water limit. Application of the model to a catchment located in Western Australia showed that the method can predict annual value of actual evapotranspiration and streamflow accurately. Copyright © 2012 John Wiley & Sons, Ltd.     

TDR measurement of stem and soil water content in two Mediterranean oak species

Abstract

Since 1990s, time domain reflectometry (TDR) has been applied to estimate the stem water content of living trees. Here, new calibration equations relating the apparent dielectric constant (K_a ) to the volumetric water content (θ) were developed for two Mediterranean oak species. Our calibration equations differ from those previously calculated for other species, suggesting that stem water contents could be monitored more accurately using species-specific curves. The stem water content in the trees of these species and the surrounding soil were monitored with TDR to examine the feasibility of this technology for recording changes in trunk water storage. The average stem water contents of the oaks reflect the soil water contents, and the temporal differences observed (17%) point to the importance of trunk water for coping with soil water deficit. Although it would be very useful to obtain a single function to estimate the stem water content of trees, it remains necessary to obtain the results in more species.     

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.     

Soil moisture states,lateral flow,and streamflow generation in a semi‐arid,snowmelt‐driven catchment

Hydraulic connectivity on hillslopes and the existence of preferred soil moisture states in a catchment have important controls on runoff generation. In this study we investigate the relationships between soil moisture patterns, lateral hillslope flow, and streamflow generation in a semi‐arid, snowmelt‐driven catchment. We identify five soil moisture conditions that occur during a year and present a conceptual model based on field studies and computer simulations of how streamflow is generated with respect to the soil moisture conditions. The five soil moisture conditions are (1) a summer dry period, (2) a transitional fall wetting period, (3) a winter wet, low‐flux period, (4) a spring wet, high‐flux period, and (5) a transitional late‐spring drying period. Transitions between the periods are driven by changes in the water balance between rain, snow, snowmelt and evapotranspiration. Low rates of water input to the soil during the winter allow dry soil regions to persist at the soil–bedrock interface, which act as barriers to lateral flow. Once the dry‐soil flow barriers are wetted, whole‐slope hydraulic connectivity is established, lateral flow can occur, and upland soils are in direct connection with the near‐stream soil moisture. This whole‐slope connectivity can alter near‐stream hydraulics and modify the delivery of water, pressure, and solutes to the stream. Copyright © 2005 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.     

Influence of antecedent soil moisture on hydraulic conductivity in a series of texture‐contrast soils

Antecedent soil moisture significantly influenced the hydraulic conductivity of the A1, A2e and B21 horizons in a series of strong texture‐contrast soils. Tension infiltration at six supply potentials demonstrated that in the A1 horizon, hydraulic conductivity was significantly lower in the ‘wet’ treatment than in the ‘dry’ treatment. However in the A2e horizon, micropore and mesopore hydraulic conductivity was lower in the ‘dry’ treatment than the ‘wet’ treatment, which was attributed to the precipitation of soluble amorphous silica. In the B21 horizon, desiccation of vertic clays resulted in the formation of shrinkage cracks which significantly increased near‐saturated hydraulic conductivity and prevented the development of subsurface lateral flow in the ‘dry’ treatment. In the ‘wet’ treatment, the difference between the hydraulic conductivity of the A1 and B21 horizons was reduced; however, lateral flow still occurred in the A1 horizon due to difficulty displacing existing soil water further down the soil profile. Results demonstrate the need to account for temporal variation in soil porosity and hydraulic conductivity in soil‐water model conceptualisation and parameterisation. Copyright © 2011 John Wiley & Sons, Ltd.     

Improving soil moisture accounting and streamflow prediction in SWAT by incorporating a modified time‐dependent Curve Number method

The objective of this study is to incorporate a time‐dependent Soil Moisture Accounting (SMA) based Curve Number method (SMA_CN) in Soil and Water Assessment Tool (SWAT) and compare its performance with the existing CN method in SWAT by simulating the hydrology of two agricultural watersheds in Indiana, USA. Results show that fusion of the SMA_CN method causes decrease in runoff volume and increase in profile soil moisture content, associated with larger groundwater contribution to the streamflow. In addition, the higher amount of moisture in the soil profile slightly elevates the actual evapotranspiration. The SMA‐based SWAT configuration consistently produces improved goodness‐of‐fit scores and less uncertain outputs with respect to streamflow during both calibration and validation. The SMA_CN method exhibits a better match with the observed data for all flow regimes, thereby addressing issues related to peak and low flow predictions by SWAT in many past studies. Comparison of the calibrated model outputs with field‐scale soil moisture observations reveals that the SMA overhauling enables SWAT to represent soil moisture condition more accurately, with better response to the incident rainfall dynamics. While the results from the modification of the CN method in SWAT are promising, more studies including watersheds with various physical and climatic settings are needed to validate the proposed approach. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of near surface soil moisture in a sloping terrain of a Himalayan watershed using ENVISAT ASAR multi‐incidence angle alternate polarisation data

Soil moisture is one of the important input variables in hydrological and water erosion models. The extraction of information on near surface soil moisture from synthetic aperture radar (SAR) is well established mostly for flat terrain and using low incidence angle single polarisation data. The ENVISAT advanced SAR (ASAR) data available in multiple incidence angles and alternate polarisation modes were investigated in this study for soil moisture estimation in sloping terrain. The test site was Sitla Rao watershed in the Lesser Himalayas of northern India. Empirical models were developed to estimate near surface soil moisture in bare agricultural fields using alternate polarisation ASAR data. Both soil moisture and surface roughness field measurements were performed during the satellite passes. Backscatter from medium incidence angle (IS‐4) and vertical‐vertical (VV) polarisation signal is correlated better with volumetric soil moisture content compared to other incidence angles. The model parameters were further improved, and soil moisture estimation was refined by combining medium incidence angle (IS4) vertical‐horizontal polarisation response as another variable along with VV polarisation response. The effect of slope on the radar backscatter was minimized by incorporating local incidence angles derived from an ASTER DEM. Copyright © 2012 John Wiley & Sons, Ltd.     

Complexity‐reduction modelling for assessing the macro‐scale patterns of historical soil moisture in the Euro‐Mediterranean region

Complexity‐reduction modelling can be useful for increasing the understanding of how the climate affects basin soil moisture response upon historical times not covered by detailed hydrological data. For this purpose, here is presented and assessed an empirical regression‐based model, the European Soil Moisture Empirical Downscaling (ESMED), in which different climatic variables, easily available on the web, are addressed for simplifying the inherent complexity in the long‐time studies. To accommodate this simplification, the Palmer Drought Severity Index, the precipitation, the elevation and the geographical location were used as input data in the ESMED model for predicting annual soil moisture budget. The test area was a large region including central Europe and Mediterranean countries, and the spatial resolution was initially set at 50 km. ESMED model calibration was made according to the soil moisture values retrieved from the Terrestrial Water Budget Data archive by selecting randomly 285 grid points (out of 2606). Once parameterized, ESMED model was performed at validation stage both spatially and temporally. The spatial validation was made for the grid points not selected in the calibration stage while the comparison with the soil moisture outputs of the Global Land Data Assimilation System–NOAH10 simulations upon the period 1950–2010 was carried out for the temporal validation. Moreover, ESMED results were found to be in good agreement with a root‐zone soil moisture product obtained from active and passive microwave sensors from various satellite missions. ESMED model was thus found to be reliable for both the temporal and spatial validations and, hence, it might represent a useful tool to characterize the long‐term dynamics of soil moisture–weather interaction. Copyright © 2013 John Wiley & Sons, Ltd.     

Synergistic effect of vegetation and air temperature changes on soil water content in alpine frost meadow soil in the permafrost region of Qinghai‐Tibet

Seasonal changes over 2 years (2004–2006) in soil moisture content (θ_v) of frozen alpine frost meadow soils of the Qinghai‐Tibet plateau permafrost region under three different levels of vegetation cover were investigated. Vegetation cover and air temperature changes had significant effects (synergistic effect) on θ_v and its distribution in the soil profile. During periods of soil freezing or thawing, the less the vegetation cover, the quicker the temperature drop or rise of soil water, and the shorter the duration of the soil water freeze–thaw response in the active soil layer. Under 30% and 65% vegetation cover the amplitude of variation in θ_v during the freezing period was 20–26% greater than that under 93% cover, while during the thawing period, it was 1·5‐ to 40·5‐fold greater. The freezing temperature of the surface soil layer, ^fT_s, was 1·6 °C lower under 30% vegetation cover than under 93% vegetation cover. Changes in vegetation cover of the alpine frost meadow affected θ_v and its distribution, as well as the relationship between θ_v and soil temperature (T_s). As vegetation cover decreased, soil water circulation in the active layer increased, and the response to temperature of the water distribution across the soil profile was heightened. The quantity of transitional soil phase water at different depths significantly increased as vegetation cover decreased. The influence of vegetation cover and soil temperature distribution led to a relatively dry soil layer in the middle of the profile (0·70–0·80 m) under high vegetation cover. Alpine meadow θ_v and its pattern of distribution in the permafrost region were the result of the synergistic effect of air temperature and vegetation cover. Copyright © 2007 John Wiley & Sons, Ltd.     

Soil sealing and soil water content under no‐tillage and conventional tillage in irrigated corn: Effects on grain yield

The main objective of this research was to analyse the effect of soil management on soil sealing and on soil water content under contrasting tillage practices and its influence on corn yield. The experimental research was carried out in a field cultivated with irrigated corn differentiated into three zones representing a gradient of soil texture (Z1, Z2, and Z3, i.e., increasingly coarser). Two plots under different soil management practices (conventional intensive tillage, CT, and no‐tillage, NT) were selected in each zone. The susceptibility to sealing of each soil and the steady infiltration rates were evaluated in the laboratory subjecting the soils to rainfall simulation applied at an intensity of 25 mm h^?1. In addition, soil porosity under each treatment was quantified. Soil water content (0–90 cm depth) was determined gravimetrically at the beginning and the end of the growing cycle and at the surface (0–5 cm) during three growing seasons and continuously at two depths (5–15 and 50–60 cm) during the last growing cycle. Soil water content was simulated using the SIMPEL model, which was calibrated for the experimental conditions. Corn yield and above‐ground biomass were also analysed. Significant differences in soil sealing among zones, with decreasing soil sealing for coarser textures, and treatments were observed with infiltration rates that were near twice in NT than in CT, being the effect of soil cover significant in the reduction of soil detachment and soil losses. NT showed higher soil water content than CT, especially in the surface layers. Above‐ground biomass production was smaller in CT than in NT, and in the areas with higher sealing susceptibility was 30% to 45% smaller than in other zones, reaching the smallest values in Z1. A similar reduction in corn yield was observed between treatments being smaller in CT than in NT. No‐tillage has been confirmed as an effective technique that benefits soil physical properties as well as crop yields in relation to CT, being its impact greater in soils susceptible to sealing.     

Mesoscale spatial variability of soil‐water content in an alpine meadow on the northern Tibetan Plateau

Soil water is an important limiting factor for restoring alpine meadows on the northern Tibetan Plateau. Field studies of soil‐water content (SWC), however, are rare due to the harsh environment, especially in a mesoscale alpine‐meadow ecosystem. The objective of this study was to assess the spatial variability of SWC and the temporal variation of the spatial variability in a typical alpine meadow using a geostatistical approach. SWC was measured using a neutron probe to a depth of 50 cm at 113 locations on 22 sampling occasions in a 33.5‐hm² alpine meadow during the 2015 and 2016 growing seasons. Mean SWC in the study plot for the two growing seasons was 18.7, 14.0, 13.9, 14.3, and 14.8% for depths of 10, 20, 30, 40, and 50 cm, respectively, and SWC was significantly larger at 10 cm than at other depths. SWC was negatively correlated with its spatial variability, and the spatial variability was higher when SWC was lower. Thirty‐three sampling locations in this study plot met the requirement of accuracy of the central limit theorem. A Gaussian model was the best fit for SWC semivariance at depths of 10, 20, and 30 cm, and the spatial structural ratio was between 0.997 and 1, indicating a strong spatial dependence of SWC. The sill and range fluctuated temporally, and the nugget and spatial structural ratio did not generally vary with time. The sill was significantly positively correlated with SWC and was initially stable and then tend to increase with SWC. The nugget, range, and spatial structure ratio, however, were not correlated with SWC. These results contribute to our understanding of SWC spatial distribution and variation in alpine meadows and provide basic empirical SWC data for mesoscale model simulations, optimizing sampling strategies and managing meadows on the Tibetan Plateau.     

Characterizing solute transport in undisturbed soil cores using electrical and X‐ray tomographic methods

Solute transport in undisturbed soil is a complex process and detailed information on the transport characteristics is needed to provide fundamental understanding of the processes involved. X‐ray computer tomography (CT) and electrical resistivity tomography (ERT) have been used to gain information on the transport characteristics. Both methods are non‐intrusive and do not disturb the soil, in contrast to other methods. CT provides high resolution information on bulk density and macropores, while ERT provides a three‐dimensional image of the internal resistivity structure. By adding a suitable solute under steady‐state flow, the internal resistivity changes can be interpreted as a change in resident concentrations. In our experiment two cores from different field sites were investigated. The ERT measurements revealed two transport modes (one fast and one slow) in one of the cores and only one mode in the other. This was consistent with the results of transfer function modelling on the independently measured breakthrough curves (BTCs). The fast transport mode is perhaps a result of many connected macropores, detected by CT, but this could not be verified with the ERT measurements because of the coarser resolution. However, with ERT in both cases we were able to explain the observed BTC qualitatively. Copyright © 1999 John Wiley & Sons, Ltd.