Responses of soil water dynamic processes and groundwater recharge to irrigation intensity and antecedent moisture in the vadose zone

In this study, a field experiment was conducted to investigate the soil water dynamics and water percolation through the deep vadose zone. A calibrated HYDRUS‐1D model was used to simulate the process of soil water movement and the water budget. Based on the measured volumetric soil water contents, the model was well calibrated and validated. Then, we conducted scenario analyses to determine the combined effects of irrigation amount (IA), antecedent soil moisture (AM), crop evapotranspiration, and deep percolation (DP) in an irrigation event. Four IAs (5, 10, 15, and 20 cm) and three AM conditions (AM‐1, AM‐2, and AM‐3) were controlled in the scenario analyses. The results indicate that according to the S_e's (effective saturation) values status and the observed or simulated depth, there could be different conclusions on the influence of DP. Under different IAs in dry (AM‐1) and medium (AM‐2) AM status, DP changed slightly; it was 0.39 and 2.47 cm in AM‐1 and 0.40 and 2.48 cm in AM‐2 for the summer maize and winter wheat crop, respectively; the AM had a crucial contribution to DP. While under the condition of wet AM (AM‐3) or small observation depth, the water inputs could have a significant effect on DP. According to increasing irrigation intensity, the higher values of S_e (>0.6) in the whole profile were only displayed between 70 and 300 cm at AM‐1, 70–500 cm at AM‐2, and 70‐below 600 cm at AM‐3, which were gradually extended and moved down with increasing AM. Hence, the IA significantly affected the water percolation at a depth of 200 cm, whereas there was a weak influence at 600 cm except in AM‐3. Furthermore, in the higher values of the S_e (>0.65) domain, the correlation between IA and DP was an exponential function and significantly under P < 0.05. In addition, DP began to occur when the soil water content was equal to or greater than 0.75 times that of the field water capacity or the S_e > 0.65. When the coarse silt layer became embedded in the silt clay soil profile, it lagged the process of water transport but did not affect permeability in the end.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 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.     

DIELECTRIC PROPERTIES OF WET SOILS IN THE FREQUENCY RANGE 1–3000 MHz1

The effective relative dielectric constant ?_{e, r} and the effective conductivity σ_e have each been determined as a function of frequency in the range 1–3000 MHz at volumetric water contents of up to approximately 0.74 for clays, 0.83 for a peat and 0.56 for a silt. At frequencies above about 25 MHz (depending on soil type), ?_{e, r}increases with water content for all samples. However, at lower frequencies, ?_{e, r}only increases with water content as long as the wet density also increases, which is the case for water contents up to a critical value lying between 0.35 and 0.48. At higher water contents, ?_{e, r}and the wet density decrease with increasing water content. Consequently, curves of ?_{e, r}versus frequency for two wet samples with different water contents, at least one of them higher than the critical value, are seen to cross at about 25 MHz. Below the critical value the curve of the sample with the lower water content is below the other curve at all freqencies applied. At a given frequency, σ_e has a maximum as a function of water content. This is tentatively explained by assuming that σ_e is the sum of pore water conductivity (increasing with water content until all salts in the soil are dissolved into the water and then decreasing) and surface water conductivity (increasing with wet density and therefore increasing with water content up to the critical value and then decreasing). At frequencies higher than 1000 MHz, ?_{e, r}depends only weakly on salinity (which is represented by the measured conductivity). It shows an increasing dependence if the frequency is decreased towards 1 MHz. The highest values of ?_{e, r}and σ_e, measured in this work, occur for a sample of wet, nearly saturated silt originating from a location below sea-level near to the Dead Sea, Israel: ?_{e, r}decreases continuously from a value of about 10⁴ at 3 MHz to about 10² at 200 MHz, while σ_e rises from about 4 S/m to 5 S/m at these respective frequencies. The dependence of the wavelength on the loss-tangent is strong and the wavelength is considerably smaller than it would be in a dielectric. This is the only sample for which σ_e increases with water content, even if the latter is above its critical value. Therefore it is assumed that the pore water conductivity is greater than the surface water conductivity if the volumetric water content is lower than 0.564, the maximum value applied. The measurements give evidence for the presence of a relaxation at about 3 MHz for all samples examined.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Variation in the vertical zonality of erodibility and critical shear stress of rill erosion in China's Hengduan Mountains

Soil erodibility is an essential parameter used in soil erosion prediction. This study selected the Liangshan town watershed to quantify variation in the vertical zonality of rill erodibility (k_r) in China's ecologically fragile Hengduan Mountains. Soil types comprised of yellow–brown (soil A), purple (soil B), and dry-red (soil C) in a descending order of occurrence from the summit to the valley, which roughly corresponds to the vertical climate zone (i.e. cool-high mountain, warm-low mountain, and dry-hot valley sections) of the study area. With elevation, vertical soil zonality varied in both soil organic matter (SOM) content and soil particle-size fractions. A series of rill erosion-based scour experiments were conducted, using water discharge rates of 100, 200, 300, 400, 500, and 600 mL min^-1. Additionally, detachment rates (Dr) were measured under three hydrological conditions (the drainage, saturation, and seepage treatments). Results show that both Dr and flow shear stress (Ʈ) values increased as discharge increased. As elevation increased, the k_r values decreased, while the vertical zonality of critical shear stress (Ʈ_c) values showed no obvious variation. The highest k_r values were observed during the seepage treatment, followed by the saturation treatment then drainage treatment, indicating that variation in vertical hydraulic gradients could significantly alter k_r values. This study also found that land-use types could also alter k_r and Ʈ_c values. Further research, however, is necessary to better quantify the effects of subsurface hydrological conditions and land-use types on k_r under different soil zonalities in China's Hengduan Mountains. © 2018 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.     

A soil‐water‐balance approach to quantify groundwater recharge from irrigated cropland in the North China Plain

Rapidly depleting unconfined aquifers are the primary source of water for irrigation on the North China Plain. Yet, despite its critical importance, groundwater recharge to the Plain remains an enigma. We introduce a one‐dimensional soil‐water‐balance model to estimate precipitation‐ and irrigation‐generated areal recharge from commonly available crop and soil characteristics and climate data. To limit input data needs and to simplify calculations, the model assumes that water flows vertically downward under a unit gradient; infiltration and evapotranspiration are separate, sequential processes; evapotranspiration is allocated to evaporation and transpiration as a function of leaf‐area index and is limited by soil‐moisture content; and evaporation and transpiration are distributed through the soil profile as exponential functions of soil and root depth, respectively. For calibration, model‐calculated water contents of 11 soil‐depth intervals from 0 to 200 cm were compared with measured water contents of loam soil at four sites in Luancheng County, Hebei Province, over 3 years (1998–2001). Each 50‐m² site was identically cropped with winter wheat and summer maize, but received a different irrigation treatment. Average root mean‐squared error between measured and model‐calculated water content of the top 180 cm was 4·2 cm, or 9·3% of average total water content. In addition, model‐calculated evapotranspiration compared well with that measured by a large‐scale lysimeter. To test the model, 12 additional sites were simulated successfully. Model results demonstrate that drainage from the soil profile is not a constant fraction of precipitation and irrigation inputs, but rather the fraction increases as the inputs increase. Because this drainage recharges the underlying aquifer, improving irrigation efficiency by reducing seepage will not reverse water‐table declines. Copyright © 2003 John Wiley & Sons, Ltd.     

TDR pressure cell for monitoring water content retention and bulk electrical conductivity curves in undisturbed soil samples

The water retention curve (θ(ψ)), which defines the relationship between soil volumetric water content (θ) and matric potential (ψ), is of paramount importance in characterizing the hydraulic behaviour of soils. However, few methods are so far available for estimating θ(ψ) in undisturbed soil samples. We present a new design of TDR‐pressure cell (TDR‐Cell) for estimating θ(ψ) in undisturbed soil samples. The TDR‐Cell consists of a 50‐mm‐long and 50‐mm internal diameter stainless steel cylinder (which constitutes the outer frame of a coaxial line) attached to a porous ceramic disc and closed at the ends with two aluminium lids. A 49‐mm‐long and 3‐mm‐diameter stainless steel rod, which runs longitudinally through the centre of the cylinder, constitutes the inner rod of a coaxial TDR probe. The TDR‐Cell was used to determine the θ(ψ) curves of a packed sand and seven undisturbed soil samples from three profiles of agricultural soils. These θ(ψ) curves were subsequently compared to those obtained from the corresponding 2‐mm sieved soils using the pressure plate method. Measurements of bulk electrical conductivity, σ_a, as a function of the water content, σ_a(θ), of the undisturbed soil samples were also performed. An excellent correlation (R² = 0·988) was found between the θ values measured by TDR on the different undisturbed soils and the corresponding θ obtained from the soil gravimetric water content. A typical bimodal θ(ψ) function was found for most of the undisturbed soil samples. Comparison between the θ(ψ) curves measured with the TDR‐Cell and those obtained from the 2‐mm sieved soils showed that the pressure plate method overestimates θ at low ψ values. The σ_a(θ) relationship was well described by a simple power expression (R² > 0·95), in which the power factor, defined as tortuosity, ranged between 1·18 and 3·75. Copyright © 2011 John Wiley & Sons, Ltd.     

Alternative analysis of transient infiltration experiment to estimate soil water repellency

The repellency index (RI) defined as the adjusted ratio between soil‐ethanol, S_e, and soil‐water, S_w, sorptivities estimated from minidisk infiltrometer experiments has been used instead of the widely used water drop penetration time and molarity of ethanol drop tests to assess soil water repellency. However, sorptivity calculated by the usual early‐time infiltration equation may be overestimated as the effects of gravity and lateral capillary are neglected. With the aim to establish the best applicative procedure to assess RI, different approaches to estimate S_e and S_w were compared that make use of both the early‐time infiltration equation (namely, the 1 min, S1, and the short‐time linearization approaches), and the two‐term axisymmetric infiltration equation, valid for early to intermediate times (namely, the cumulative linearization and differentiated linearization approaches). The dataset included 85 minidisk infiltrometer tests conducted in three sites in Italy and Spain under different vegetation habitats (forest of Pinus pinaster and Pinus halepensis, burned pine forest, and annual grasses), soil horizons (organic and mineral), postfire treatments, and initial soil water contents. The S1 approach was inapplicable in 42% of experiments as water infiltration did not start in the first minute. The short‐time linearization approach yielded a systematic overestimation of S_e and S_w that resulted in an overestimation of RI by a factor of 1.57 and 1.23 as compared with the cumulative linearization and differentiated linearization approaches. A new repellency index, RI_s, was proposed as the ratio between the slopes of the linearized data for the wettable and hydrophobic stages obtained by a single water infiltration test. For the experimental conditions considered, RI_s was significantly correlated with RI and WDPT. Compared with RI, RI_s includes information on both soil sorptivity and hydraulic conductivity and, therefore, it can be considered more physically linked to the hydrological processes affected by soil water repellency.     

A bimodal four‐parameter lognormal linear model of soil water repellency persistence

Soil water repellency may be characterized in terms of the delayed infiltration time of a water droplet resting on the soil surface, which is, water drop penetration time (WDPT), or repellency persistence. Such repellency persistence varies nonlinearly with soil water content (θ_g), although no models have been proposed to reproduce the variation of WDPT with θ_g in soils. Dynamic factor analysis (DFA) is used to identify two common patterns of unexplained variability in a scattered dataset of WDPT versus θ_g measurements. A four‐parameter lognormal distribution was fitted to both common patterns obtained by DFA, and these were combined additively in a weighted multiple linear bimodal model. We show how such an empirical model is capable of reproducing a large variety of WDPT versus θ_g curve shapes (N = 80) both within a wide range of measured WDPTs (0–17 000 s) and for samples with organic matter content ranging from 21·7 to 80·6 g (100 g)^?1. Copyright © 2009 John Wiley & Sons, Ltd.     

Recharge,upflux and water table response for shallow water table conditions in southwest Florida

A disproportionate increase or decrease in water table in response to minor water input or drainage is observed in shallow water table conditions inside drainage lysimeters. This increase happens because the capillary fringe of the shallow water table reaches up to or near the surface (Wieringermeer effect). The correlations between water table level changes and rainfall, seepage irrigation, drip irrigation, and drainage were analysed. Correlations with rainfall, seepage irrigation, and drainage were high (R² ranged from 0·46 to 0·97). Drip irrigation had low correlations due to the low rates of application (R² ranged from 0·26 to 0·44). Conventional methods of calculating recharge, such as multiplying the specific yield with the water table fluctuations, cannot be used for Wieringermeer effect situations. A method using water balance data and soil moisture at different depths in the lysimeters was developed to estimate recharge and upflux. The recharge results were used to develop the apparent specific yield S_ya, which could be used to calculate consequent recharge events from water table fluctuation data. Combining the water table fluctuation relationships developed with the S_ya value will allow the prediction of recharge from rainfall and irrigation events without the need for soil moisture equipment. Copyright © 2005 John Wiley & Sons, Ltd.     

Modification of the Penman method for computing bare soil evaporation

The Penman equation, which calculates potential evaporation, was modified by Staple (1974, Soil Science Society of America Proceedings 38 : 837) to include in it the relative vapour pressure h_s of an unsaturated soil to predict actual evaporation from a soil surface. This improved the prediction when the difference between the temperature of the soil surface and ambient air is relatively small. The objectives of this study were (i) to revise it further using the actual temperature of the soil surface and air to provide the upper boundary condition in computing evaporative flux from the soil surface and (ii) to determine the range of water content for which the modified form of the Penman equation is applicable. The method adopted was tested by a series of outdoor experiments with a clay soil. The method of Staple (1974) overestimated the rate of evaporation above the water content 0·14 m³ m^?3 (up to 30% deviation), whereas the new method agreed well with the measured rates (maximum 7% deviation). Below 0·14 m³ m^?3 water content, both methods underestimated, but the Staple (1974) method deviated more from the measured values: the deviations were above 70% and around 30% for the Staple (1974) and the new methods respectively. Although the new method provided accurate solutions for a wider range of water content from saturation to the lower limit of the liquid phase of a particular soil, the modification did not respond to the vapour phase of the soil moisture. Therefore, in the dry range (i.e. in the vapour phase in which the flow was entirely as vapour), either resistance models or a Fickian equation should be used. Although the effect of salinity on the measured rates was significant, the model erroneously calculated the same rates for both saline and non‐saline conditions. The effect of soil texture can easily be accounted by defining appropriate matric potential water content ψ_m(θ) and soil relative humidity water content h_s(θ) relationships. Copyright © 2007 John Wiley & Sons, Ltd.     

Coupling two radar backscattering models to assess soil roughness and surface water content at farm scale

Abstract

Remote sensing techniques are useful for agro-hydrological monitoring at the farm scale because the availability of spatially and temporally distributed data improves agricultural models for irrigation and crop yield optimization under water scarcity conditions. This research focuses on the surface water content retrieval using active microwave data. Two semi-empirical models were chosen as these showed the best performances in simulating cross and co-polarized backscatter. Thus, these models were coupled to obtain reliable assessments of both soil water content and soil roughness. The use of the coupled model enables one to avoid using roughness measured in situ. Remote sensing images and in situ data were collected between April and July 2006 within the European Space Agency-funded project AgriSAR 2006. The images data set includes L-band in HH, VV and VH polarizations acquired from the airborne E-SAR sensor, operated by the German Aerospace Centre. Results were validated using in situ soil water content and roughness measurements. The results show that reliable assessment of both soil roughness (r ² up to ?0.8) and soil water content (r ² ? 0.9) can be retrieved in fields characterized by low fractional coverage.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Capodici, F., Maltese, A., Ciraolo, G., La Loggia, G., and D’Urso, G., 2013. Coupling two radar backscattering models to assess soil roughness and surface water content at the farm scale. Hydrological Sciences Journal, 58 (8), 1677–1689.     

Development of pedotransfer functions for soil hydraulic properties in the critical zone on the Loess Plateau,China

Soil hydraulic properties (SHPs) including the soil water retention curve and saturated soil hydraulic conductivity (Ks) are crucial input data for simulations of soil water and solute transport in the Earth's critical zone. However, obtaining direct measurements of SHPs at a wide range of scales is time consuming and expensive. Pedotransfer functions (PTFs) are employed as an alternative method for indirectly estimating these parameters based on readily measured soil properties. However, PTFs for SHPs for the deep soil layer in the Earth's critical zone are lacking. In this study, we developed new PTFs in the deep soil profile for Ks and soil water retention curve on the Loess Plateau, China, which were fitted with the van Genuchten equation. In total, 206 data sets comprising the hydraulic and basic soil properties were obtained from three typical sites. Samples were collected from the top of the soil profile to the bedrock by soil core drilling. PTFs were developed between the SHPs and basic soil properties using stepwise multiple linear regression. The PTFs obtained the best predictions for Ks (R_adj² = 0.561) and the worst for van Genuchten α (R_adj² = 0.474). The bulk density and sand content were important input variables for predicting Ks, α, and θs, and bulk density, clay content, and soil organic carbon were important for n. The PTFs developed in this study performed better than existing PTFs. This study contains the first set of PTFs of SHPs to be developed for the deep profile on the Loess Plateau, and they may be applicable to other regions.     

Seasonal climate patterns and their influence on calibration of the Hargreaves-Samani equation

Abstract

Annual patterns in climate parameters were studied to evaluate how these influence the quality of reference evapotranspiration (ETo) estimates obtained from the Hargreaves-Samani (HS) equation, since the method only uses the measured temperature directly. The work evaluates how these patterns can be used to improve the HS ETo estimates. Ten-year moving averages from a set of California Irrigation Management Information System (CIMIS) stations were used to evaluate the relationships between solar radiation (R_s), temperature (T) and ETo. The results indicate that T treads behind solar radiation and its value peaks some 25 days later. Thus, the main irrigation season in the Mediterranean climate (1 May–30 September) can be divided into three phases: increasing R_s and T; decreasing R_s with increasing T; and decreasing R_s and T. Non-univocal annual cycles were observed between R_s and T, ETo and R_s, and ETo and T. These annual patterns result in important seasonal changes in the ratio between the HS and Penman-Monteith (FAO PM) ETo estimates. The changes are particularly important during the irrigation season, where the FAO PM initially calculates greater ETo values than the HS methodology, and from the end of May to early September, where the HS equation overestimates the ETo values (by 17 mm, or 3%). These patterns obtained from 2000–2009 data were used to calibrate and improve HS ETo estimates at new sites for the 2010–2011 period. Calibration based on the proposed seasonal region-wide FAO PM/HS ETo ratios improved both the bias (decreased from 0.40 to 0.36 mm d^-1) and r² (increased from 0.67 to 0.87) of the ETo estimates for the irrigation season. The proposed methodology can be easily applied to other regions, even when the existing weather stations are sparse.

Editor Z.W. Kundzewicz     

Water flux measurement and prediction in young cashew trees using sap flow data

Measurements of sap flow, meteorological parameters, soil water content and tension were made for 4 months in a young cashew (Anacardium occidentale L.) plantation during the 2002 rainy season in Ejura, Ghana. This experiment was part of a sustainable water management project in West Africa. The Granier system was used to measure half‐hourly whole‐tree sap flow. Weather variables were observed with an automatic weather station, whereas soil moisture and tension were measured with a Delta‐T profile probe and tensiometers respectively. Clearness index (CI), a measure of the sky condition, was significantly correlated with tree transpiration (r² = 0·73) and potential evaporation (r² = 0·86). Both diurnal and daily stomata conductance were poorly correlated with the climatic variables. Estimated daily canopy conductance g_c ranged from 4·0 to 21·2 mm s⁻¹, with a mean value of 8·0 ± 3·3 mm s⁻¹. Water flux variation was related to a range of environmental variables: soil water content, air temperature, solar radiation, relative humidity and vapour pressure deficit. Linear and non‐linear regression models, as well as a modified Priestley–Taylor formula, were fitted with transpiration, and the well‐correlated variables, using half‐hourly measurements. Measured and predicted transpiration using these regression models were in good agreement, with r² ranging from 0·71 to 0·84. The computed measure of accuracy δ indicated that a non‐linear model is better than its corresponding linear one. Furthermore, solar radiation, CI, clouds and rain were found to influence tree water flux. Copyright © 2005 John Wiley & Sons, Ltd.     

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

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

Impacts of different vegetation restoration options on gully head soil resistance and soil erosion in loess tablelands

Vegetation restoration is identified as an effective approach to control soil erosion and affects soil detachment and resistance to concentrated flow on the Loess Plateau. However, the effects of vegetation restoration at gully heads in loess-tableland remains unclear. This study was performed to investigate the effects of nine vegetation restoration types at gully heads on soil detachment rate (Dr) and soil resistance to concentrated flow (i.e. soil erodibility, K_r and critical shear stress, τ_c). Undisturbed soil samples were collected from nine vegetation-restored lands and one slope cropland (as the control) and were subjected to a hydraulic flume to obtain Dr values of gully heads under six inflow discharges (0.5–3.5 L s^-1). The results showed that the Dr values of nine revegetated gully heads were 77.11% to 95.81% less than that of slope cropland, and the grassland dominated by Cleistogenes caespitosa and the shrubland dominated by Hippophae rhamnoides had a relatively greater decrease in Dr than those of other seven restoration types. The Dr value of nine revegetated gully heads could be better simulated by stream power than by flow velocity and shear stress and was also significantly affected by soil disintegration rate (positively), soil bulk density, saturated hydraulic conductivity, organic matter content, and water-stable aggregate stability (negatively). Additionally, roots with diameters of 0 to 0.5 mm showed a greater effect on Dr than those with larger diameters. Compared to cropland, the nine restored types reduced K_r by 76.26% to 94.26% and improved τ_c by 1.51 to 4.68 times. The decrease in K_r and the increase in τ_c were significantly affected by organic matter content, water-stable aggregate, mean weight diameter of aggregate and root mass density. The combination of grass species (Cleistogenes caespitosa) and shrub (Hippophae rhamnoides) could be considered the best vegetation restoration types for improving soil resistance of gully heads to concentrated flow. © 2019 John Wiley & Sons, Ltd.