Theory, Construction and Operation of Simple Tensiometers

The tensiometer, introduced by Richards and Gardner (1936) has been a tool often ignored in unsaturated zone investigations. Hazardous waste disposal studies that require values of unsaturated zone matric potentials will expand tensiometer use. Familiarity with operation and installation principles is necessary to ensure quality data collection.
The essential elements of a tensiometer are a porous cup connected with tubing to a vacuum gauge; all are filled with water. When the porous cup is buried in unsaturated soil, the soil draws water through the cup wall, until the water tension inside the cup equals the surrounding soil-water tension. The tension is transmitted through the water filled tube to an accessible vacuum gauge (manometer or transducer) and is recorded manually or automatically.
Although specialized tensiometers exist, the tensiometer presented here in detail is suited to diverse on-site applications. Constructed from readily available, inexpensive parts, it can measure as much as 0.85 bar of tension. Design features include a flushing system for removal of entrapped air or mercury, and an easily maintained modular network of nylon manometers and water-supply tubes.
Good hydraulic connection between porous cup and soil and meticulous backfill procedures are necessary for accurate measurement of representative matric potentials. Proper and consistent use of tensiometers should yield reliable values of matric potentials, and lead to a better understanding of unsaturated zone processes.     

The role of vegetation patterns in structuring runoff and sediment fluxes in drylands

The dynamics of vegetation‐driven spatial heterogeneity (VDSH) and its function in structuring runoff and sediment fluxes have received increased attention from both geomorphological and ecological perspectives, particularly in arid regions with sparse vegetation cover. This paper reviews the recent findings in this area obtained from field evidence and numerical simulation experiments, and outlines their implications for soil erosion assessment. VDSH is often observed at two scales, individual plant clumps and stands of clumps. At the patch scale, the local outcomes of vegetated patches on soil erodibility and hydraulic soil properties are well established. They involve greater water storage capacity as well as increased organic carbon and nutrient inputs. These effects operate together with an enhanced capacity for the interception of water and windborne resources, and an increased biological activity that accelerates breakdown of plant litter and nutrient turnover rates. This suite of relationships, which often involve positive feedback mechanisms, creates vegetated patches that are increasingly different from nearby bare ground areas. By this way a mosaic builds up with bare ground and vegetated patches coupled together, respectively, as sources and sinks of water, sediments and nutrients. At the stand scale within‐storm temporal variability of rainfall intensity controls reinfiltration of overland flow and its decay with slope length. At moderate rainfall intensity, this factor interacts with the spatial structure of VDSH and the mechanism of overland flow generation. Reinfiltration is greater in small‐grained VDSH and topsoil saturation excess overland flow. Available information shows that VDSH structures of sources and sinks of water and sediments evolve dynamically with hillslope fluxes and tune their spatial configurations to them. Rainfall simulation experiments in large plots show that coarsening VDSH leads to significantly greater erosion rates even under heavy rainfall intensity because of the flow concentration and its velocity increase. Copyright © 2005 John Wiley & Sons, Ltd.     

Beerkan multi-runs for characterizing water infiltration and spatial variability of soil hydraulic properties across scales

A method is presented for characterizing the spatial variability of water infiltration and soil hydraulic properties at the transect and field scales. The method involves monitoring a set of 10 Beerkan runs distributed over a 1-m length of soil, and running BEST (Beerkan estimation of soil transfer parameters) methods to derive hydraulic parameters. The Beerkan multi-runs (BMR) method provides a significant amount of data at the transect scale, allowing the determination of correlations between water infiltration variables and hydraulic parameters, and the detection of specific runs affected by preferential flow or water repellence. The realization of several BMRs at several transects on the same site allows comparison of the variation between locations (spatial variability at the field scale) and at the transect scale (spatial variability at the metre scale), using analysis of variance. From the results, we determined the spatial variability of water infiltration and hydraulic parameters as well as its characteristic scale (transect versus field).     

Is the soil moisture precipitation feedback enhanced by heterogeneity and dry soils? A comparative study

The interaction between the land surface and the atmosphere is a crucial driver of atmospheric processes. Soil moisture and precipitation are key components in this feedback. Both variables are intertwined in a cycle, that is, the soil moisture – precipitation feedback for which involved processes and interactions are still discussed. In this study the soil moisture – precipitation feedback is compared for the sempiternal humid Ammer catchment in Southern Germany and for the semiarid to subhumid Sissili catchment in West Africa during the warm season, using precipitation datasets from the Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), from the German Weather Service (REGNIE) and simulation datasets from the Weather Research and Forecasting (WRF) model and the hydrologically enhanced WRF-Hydro model. WRF and WRF-Hydro differ by their representation of terrestrial water flow. With this setup we want to investigate the strength, sign and variables involved in the soil moisture – precipitation feedback for these two regions. The normalized model spread between the two simulation results shows linkages between precipitation variability and diagnostic variables surface fluxes, moisture flux convergence above the surface and convective available potential energy in both study regions. The soil moisture – precipitation feedback is evaluated with a classification of soil moisture spatial heterogeneity based on the strength of the soil moisture gradients. This allows us to assess the impact of soil moisture anomalies on surface fluxes, moisture flux convergence, convective available potential energy and precipitation. In both regions the amount of precipitation generally increases with soil moisture spatial heterogeneity. For the Ammer region the soil moisture – precipitation feedback has a weak negative sign with more rain near drier patches while it has a positive signal for the Sissili region with more rain over wetter patches. At least for the observed moderate soil moisture values and the spatial scale of the Ammer region, the spatial variability of soil moisture is more important for surface-atmosphere interactions than the actual soil moisture content. Overall, we found that soil moisture heterogeneity can greatly affect the soil moisture – precipitation feedback.     

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.     

Unsaturated flows in soils with self-similar hydraulic conductivity distribution

In this article, we are concerned with the statistics of steady unsaturated flow in soils with a fractal hydraulic conductivity distribution. It is assumed that the spatial distribution of log hydraulic conductivity can be described as an isotropic stochastic fractal process. The impact of the fractal dimension of this process, the soil pore-size distribution parameter, and the characteristic length scale on the variances of tension head and the effective conductivity is investigated. Results are obtained for one-dimensional and three-dimensional flows. Our results indicate that the tension head variance is scale-dependent for fractal distribution of hydraulic conductivity. Both tension head variance and effective hydraulic conductivity depend strongly on the fractal dimension. The soil pore-size distribution parameter is important in reducing the variability of the unsaturated hydraulic conductivity and of the fluxes.     

The effect of beech stemflow on spatial patterns of soil solution chemistry and seepage fluxes in a mixed beech/oak stand

Stemflow of beech (Fagus sylvatica L.) represents a significant input of water and elements to the soil and might influence the spatial patterns and the rate of seepage fluxes at the stand scale. We investigated the soil solution chemistry at different depths and distances from the stem and the element fluxes with stemflow, throughfall and seepage in proximal and distal stem areas of a 130‐year‐old beech/oak forest in Steigerwald (northern Bavaria, Germany). The proximal stem area (in total 286 m² ha⁻¹) was defined as a 1 m², 60 cm deep cylinder around the beech stem. Seepage fluxes were calculated by a soil hydrological model for 1996 using measured soil matrix potentials and tree xylem flow data for calibration. Stemflow represented 6·6% of the annual soil water input. With the exception of H⁺ fluxes, less than 10% of the total element fluxes with throughfall and stemflow reached the soil via stemflow. The volume‐weighted concentrations of H⁺, K⁺ and SO₄²⁻ in stemflow were higher than those in throughfall, while other elements had similar concentrations. Soil solution K⁺ concentrations decreased with stem distance, but the Na⁺, Mg²⁺, Cl⁻ and SO₄²⁻ concentrations increased. Gradients for other elements were not statistically significant. Stemflow had a strong influence on the spatial patterns of element fluxes with seepage. The water fluxes through the soil of the proximal stem areas at a depth of 60 cm contributed 13·5% to the total seepage at the stand scale. Proximal to the stems about 20% of total seepage for K⁺, Mn²⁺, Alⁿ⁺, dissolved organic N and dissolved organic C were concentrated, but only 8–10% for Na⁺, Mg²⁺ and Ca²⁺. The loss of acid‐neutralizing capacity calculated from the flux balance was about four times higher proximal to the stems compared with distal areas, indicating high rates of soil acidification proximal to the stems. Our results confirm the concept of a microsite around beech stems, characterized by high element and water fluxes in comparison with distal stem areas. Calculations of seepage fluxes and element budgets in beech stands have to consider the spatial heterogeneity of fluxes induced by stemflow. Copyright © 2000 John Wiley & Sons, Ltd.     

Physical controls on septic leachate movement in the vadose zone at the hillslope scale,Putnam County,New York,USA

The fate and transport of contaminants in the vicinity of septic fields remains poorly understood in many hydrogeomorphological environments. We report hydrometric data from an intensive hillslope‐scale experiment conducted between 29 August and 11 November 1998 at a residential leach field in New York State. The objective of our study was to characterize water flux within the vadose zone, understand the physical controls on the flux, and predict how this ultimately will affect subsurface water quality. Soil‐water flux was calculated using matric potential measurements from a network of 25 tensiometer nests, each nest consisting of three tensiometers installed to depths of 10, 50 and 130 cm. Unsaturated hydraulic conductivity curves were derived at each depth from field‐determined time‐domain reflectometry–tensiometry moisture‐release curves and borehole permeametry measurements. Flownets indicated that a strong upward flux of soil water occurred between rainstorms. Following the onset of (typically convective) rainfall, low near‐surface matric potentials were rapidly converted to near‐saturated and saturated conditions, promoting steep vertical gradients through the near‐surface horizons of the hillslope. Lateral hydraulic gradients were typically 10 times smaller than the vertical gradients. Resultant flow vectors showed that the flux was predominantly vertical through the vadose zone, and that the flux response to precipitation was short‐lived. The flux response was controlled primarily by the shape of the unsaturated hydraulic conductivity curves, which indicated a rapid loss of conductivity below saturation. Thus, soil water had a very high residence time in the vadose zone. The absence of rapid wetting at 130 cm and the delayed and small phreatic zone response to rainfall indicated that water movement through macropores did not occur on this hillslope. These results are consistent with a Cl tracing experiment, which demonstrated that the tracer was retained in the vadose zone for several months after injection to the system. Copyright © 2002 John Wiley & Sons, Ltd.     

Remote sensing and ground‐based measurements of evapotranspiration in an extreme cold Patagonian desert

Accurate estimates of seasonal evapotranspiration (ET) at different temporal and spatial scales are essential for understanding the biological and environmental determinants of ecosystem water balance in arid regions and the patterns of water utilization by the vegetation. For this purpose, remote sensing ET estimates of a Patagonian desert in Southern Argentina were verified with field measurements of soil evaporation and plant transpiration using an open top chamber. Root distribution and seasonal variation in soil volumetric water content were also analysed. There was a high correlation between remote sensing and field measurements of ecosystem water fluxes. A substantial amount of the annual ET occurred in spring and early summer (73.4 mm) using winter rain stored in the soil profile and resulting in water content depletion of the upper soil layers. A smaller amount of annual ET was derived from few rainfall events occurring during the mid or late summer (41.4 mm). According to remote sensing, the 92.9% of the mean annual precipitation returns to the atmosphere by transpiration or evaporation from the bare soil and by canopy interception. Only 7.1% infiltrates to soil layers deeper than 200 cm contributing to the water table recharge. Fourier time series analysis, cross‐correlation methods and multiple linear regression models were used to analyse 11 years of remote sensing data to assess determinants of water fluxes. A linear model predicts well the variables that drive complex ecosystem processes such as ET. Leaf area index and air temperature were not linearly correlated to ET because of the multiple interaction among variables resulting in time lags with ET variations and thus these two variables were not included in the linear model. Soil water content, the fraction of photosynthetic active radiation and precipitation explained 86% of the ET monthly variations. The high volumetric water content and the small seasonal variations at 200‐cm depth were probably the result of little water uptake from deeper soil horizons by roots with low hydraulic conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Monitoring over‐winter soil water dynamics in a freezing and snow‐covered environment using a thermally insulated tensiometer

Water potential below a frozen soil layer was continuously monitored over an entire winter period (using thermally insulated tensiometers sheltered in a heated chamber) along with other soil, snow and atmospheric variables. In early winter, the freezing front advanced under a thin snow cover, inducing upward soil water flow in the underlying unfrozen soil. The freezing front started to retreat when the snow cover became thick enough to insulate the soil, resulting in the reversal of the flow direction in the unfrozen zone. These data provide a clear illustration of soil water dynamics, which have rarely been monitored with a tensiometer. Copyright © 2005 John Wiley & Sons, Ltd.     

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

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

Modelling the impacts of spatial heterogeneity in soil hydraulic properties on hydrological process in the upper reach of the Heihe River in the Qilian Mountains,Northwest China

Spatial heterogeneity of soil has great impacts on dynamic processes of the hydrological systems. However, it is challenging and expensive to obtain spatial distribution of soil hydraulic properties, which often requires extensive soil sampling and observations and intensive laboratory analyses, especially in high elevation, hard to access mountainous areas. This study evaluates the impacts of soil heterogeneity on hydrological process in a high elevation, topographically complex watershed in Northwest China. Two approaches were used to derive the spatial heterogeneity of soil properties in the study watershed: (1) the spatial clustering method, Full‐Order‐CLK was used to determine five soil heterogeneous clusters (configurations 97, 80, 60, 40 and 20) through large number of soil sampling and in situ observations, and (2) the average values of soil hydraulic properties for each soil type were derived from the coarse provincial soil data sets (Gansu Soil Handbook at 1 : 1 000 000 scale). Subsequently, Soil and Water Assessment Tool model was used to quantify the impact of the spatial heterogeneity of soil hydraulic properties on hydrological process in the study watershed. Results show the simulations by Soil and Water Assessment Tool with the spatially clustered soil hydraulic information from the field sampling data had much better representation of the soil heterogeneity and had more accurate performance than the model using the average soil property values for each soil type derived from the coarse soil data sets. Thus, incorporating detailed field sampling, soil heterogeneity data greatly improve performance in hydrological modelling. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatial and temporal analysis of groundwater recharge with application to sampling design

A numerical experiment of flow in variably saturated porous media was performed in order to evaluate the spatial and temporal distribution of the groundwater recharge at the phreatic surface for a shallow aquifer as a function of the input rainfall process and soil heterogeneity. The study focused on the groundwater recharge which resulted from the percolation of the excess rainfall for a 90-days period of an actual precipitation record. Groundwater recharge was defined as the water flux across the moving phreatic surface. The observed spatial non-uniformity of the groundwater recharge was caused by soil heterogeneity and is particularly pronounced during the stage of recharge peak (substantial percolation stage). During that stage the recharge is associated with preferential flow paths defined as soil zones of locally higher hydraulic conductivity. For the periods of low percolation intensity the groundwater recharge was exhibiting more uniform spatial characteristics. The temporal distribution of the recharge was found to be a function of the frequency and intensity of the rainfall events. Application of sampling design demonstrates the joint influence of the spatial and temporal recharge variability on the cost-effective monitoring of groundwater potentiometric surfaces.     

How much does dry‐season fog matter? Quantifying fog contributions to water balance in a coastal California watershed

The seasonally‐dry climate of Northern California imposes significant water stress on ecosystems and water resources during the dry summer months. Frequently during summer, the only water inputs occur as non‐rainfall water, in the form of fog and dew. However, due to spatially heterogeneous fog interaction within a watershed, estimating fog water fluxes to understand watershed‐scale hydrologic effects remains challenging. In this study, we characterized the role of coastal fog, a dominant feature of Northern Californian coastal ecosystems, in a San Francisco Peninsula watershed. To monitor fog occurrence, intensity, and spatial extent, we focused on the mechanisms through which fog can affect the water balance: throughfall following canopy interception of fog, soil moisture, streamflow, and meteorological variables. A stratified sampling design was used to capture the watershed's spatial heterogeneities in relation to fog events. We developed a novel spatial averaging scheme to upscale local observations of throughfall inputs and evapotranspiration suppression and make watershed‐scale estimates of fog water fluxes. Inputs from fog water throughfall (10–30 mm/year) and fog suppression of evapotranspiration (125 mm/year) reduced dry‐season water deficits by 25% at watershed scales. Evapotranspiration suppression was much more important for this reduction in water deficit than were direct inputs of fog water. The new upscaling scheme was analyzed to explore the sensitivity of its results to the methodology (data type and interpolation method) employed. This evaluation suggests that our combination of sensors and remote sensing allows an improved incorporation of spatially‐averaged fog fluxes into the water balance than traditional interpolation approaches.     

Impact of soil moisture heterogeneity length scale and gradients on daytime coupled land‐cloudy boundary layer interactions

Using a coupled large‐eddy simulation–land surface model framework, the impact of two‐dimensional soil moisture heterogeneity on the cloudy boundary layer under varied free‐atmosphere stabilities is investigated. Specifically, the impacts of soil moisture heterogeneity length scale and heterogeneity in terms of soil moisture gradients on micrometeorological states, surface fluxes, boundary layer characteristics, and cloud development are examined. The results show that mesoscale circulations due to surface heterogeneity in soil moisture play an important role in transferring water vapour within the boundary layer and in regulating cloud distribution at the entrainment zone, which, in turn, provides feedbacks on boundary layer/surface energy budgets. The initial domain‐averaged soil moisture is identical for all homogenous and heterogeneous cases; however, the soil moisture heterogeneity in gradient and length scale between dry and wet regions has a significant impact on the estimates of near‐surface micrometeorological properties and surface fluxes, which further affect the boundary layer states and characteristics. Both liquid water potential temperature and liquid water mixing ratio increase with an increasing soil moisture gradient, whereas the amount of specific humidity decreases. Heterogeneity length scale and free atmosphere stability also amplify these impacts on the boundary layer structure and cloud formation. In a low atmospheric stability condition that potentially allows for a deeper boundary layer and a higher entrainment rate, cloud base height and cloud thickness significantly increase as the soil moisture gradient and length scale increase. Analysis to differentiate the influences of surface heterogeneity type (i.e. length scale vs gradient) shows that in general soil moisture gradient provides a larger impact than heterogeneity length scale, although the heterogeneity length scale is large enough to initiate circulation features responsible for differences in the coupled system between homogeneous and heterogeneous soil moisture cases. Copyright © 2012 John Wiley & Sons, Ltd.     

EVALUATION OF VERTICAL AND LATERAL FLOW THROUGH AGRICULTURAL LOESSIAL HILLSLOPES USING A TWO-DIMENSIONAL COMPUTER SIMULATION MODEL

On four hillslopes in the loess region of The Netherlands pressure heads were monitored using stand alone measuring devices. During rain events pressure heads were measured regularly with time intervals of five minutes. Rainfall was recorded with a tipping bucket rainfall gauge. Two-dimensional simulations were executed to simulate water flow through these hillslopes during erosive rain events in summer and winter periods. From these simulations it appeared that vertical flow is dominant during rain events, and that lateral water transport is of minor importance despite distinct layer stratifications. During the selected rain showers, average lateral water movement varied between 1.6 and 4.7% of the total water displacement through the hillslopes. Therefore, it was decided to incorporate a one-dimensional water flow module into the event-based distributed soil erosion and hydrological model LISEM.     

Multi-scale heterogeneity of soil moisture following snow thawing in Haloxylon ammodendron shrubland

Spatial variation of soil moisture after snow thawing in South Gurbantunggut was quantitatively studied using ANOVA and geostatistics at various scales. The results show that the soil moisture heterogeneity varies along with spatial scales. At the shrub individual scale, there is a gradient in soil moisture from shrub-canopied area to canopy margin and to the interspaces between shrubs. At the community scale, soil moisture is highly autocorrelated and the semivariogram is fitted as spherical model, with an 89.6% structural variance and a range of 4.02 m. In addition, Kringing map indicates that the soil moisture distribution pattern after snow thawing is highly consistent with the shrub patch pattern. At the typical inter-dune transect scale, soil moisture presents a pattern of high value at inter-dune depression and low value at dune, and this variation is fitted as Gaussian model with a structural variance of 95.8% and a range of 66.16 m. The range is comparable with the scale of topography zoning, suggesting that the topography pattern controls the pattern of snowmelt at this scale. The evidence indicates that the heterogeneity of soil moisture at various scales is controlled by various land surface processes after snow thawing. For Gurbantunggut Desert, the spatial heterogeneity of snowmelt at various scales is ecologically valuable, because it promotes the utilization efficiency of the snowmelt for the desert vegetation.     

Groundwater recharge and evapotranspiration for two natural ecosystems covered with oak and heather

The evapotranspiration and groundwater recharge from two natural areas with high (oak) and low (heath) vegetation were estimated by calibrating a semi-physical numerical soil water and heat model to fit 8 and 7 years of TDR-measurements of water content, respectively. The measurements were made between the surface and 7 m depth. For the oak stand, the estimated annual recharge for the years 1992–1999 is 390 mm, the evaporation from soil and interception is 205 mm, and the transpiration is 285 mm. For the heath area estimation was carried out for the years 1993–1999. However, the heath was struck by a heavy beetle attack in 1994, which strongly affected the vegetation and thus the water balance for the following 3 years. For years not affected, the estimated recharge is 733 mm (about 50% larger than for the oak stand for the same years), the evaporation is 316 mm, and the transpiration is 128 mm. The estimated recharge values compare fairly well to estimates obtained from bromide tracer experiments. However, the recharge estimates obtained from the tracer experiments are very uncertain. The uncertainty is mainly due to spatial heterogeneity making the three replicate samples taken here for each time and depth insufficient.

The analyses of TDR-measurements and tracer data showed that water front movement depends on the antecedent soil water content. Some layers are bypassed, especially at low water contents, and at high soil water contents preferential flow was observed at the heath site.     

Multi-scale heterogeneity of soil moisture following snow thawing in <Emphasis Type="Italic">Haloxylon ammodendron</Emphasis> shrubland

Spatial variation of soil moisture after snow thawing in South Gurbantunggut was quantitatively studied using ANOVA and geostatistics at various scales. The results show that the soil moisture heterogeneity varies along with spatial scales. At the shrub individual scale, there is a gradient in soil moisture from shrub-canopied area to canopy margin and to the interspaces between shrubs. At the community scale, soil moisture is highly autocorrelated and the semivariogram is fitted as spherical model, with an 89.6% structural variance and a range of 4.02 m. In addition, Kringing map indicates that the soil moisture distribution pattern after snow thawing is highly consistent with the shrub patch pattern. At the typical inter-dune transect scale, soil moisture presents a pattern of high value at inter-dune depression and low value at dune, and this variation is fitted as Gaussian model with a structural variance of 95.8% and a range of 66.16 m. The range is comparable with the scale of topography zoning, suggesting that the topography pattern controls the pattern of snowmelt at this scale. The evidence indicates that the heterogeneity of soil moisture at various scales is controlled by various land surface processes after snow thawing. For Gurbantunggut Desert, the spatial heterogeneity of snowmelt at various scales is ecologically valuable, because it promotes the utilization efficiency of the snowmelt for the desert vegetation.     

Spatio‐temporal heterogeneity in soil water stable isotopic composition and its ecohydrologic implications in semiarid ecosystems

Spatio‐temporal heterogeneity in soil water content is recognized as a common phenomenon, but heterogeneity in the hydrogen and oxygen isotope composition of soil water, which can reveal processes of water cycling within soils, has not been well studied. New advances are being driven by measurement approaches allowing sampling with high density in both space and time. Using in situ soil water vapour probe techniques, combined with conventional soil and plant water vacuum distillation extraction, we monitored the hydrogen and oxygen stable isotopic composition of soil and plant waters at paired sites dominated by grasses and Gambel's oak (Quercus gambelii) within a semiarid montane ecosystem over the course of a growing season. We found that sites spaced only 20 m apart had profoundly different soil water isotopic and volumetric conditions. We document patterns of depth‐ and time‐explicit variation in soil water isotopic conditions at these sites and consider mechanisms for the observed heterogeneity. We found that soil water content and isotopic variability were damped under Q. gambelii, perhaps due in part to hydraulic redistribution of deep soil water or groundwater by Q. gambelii in these soils relative to the grass‐dominated site. We also found some support for H isotope discrimination effects during water uptake by Q. gambelii. In this ecosystem, the soil water content was higher than that at the neighbouring Grass site, and thus, 25% more water was available for transpiration by Q. gambelii compared with the Grass site. This work highlights the role of plants in governing soil water variation and demonstrates that they can also strongly influence the isotope ratios of soil water. The resulting fine‐scale heterogeneity has implications for the use of isotope tracers to study soil hydrology and evaporation and transpiration fluxes to improve understanding of water cycling through the soil–plant–atmosphere continuum.