Information and complexity measures applied to observed and simulated soil moisture time series

Abstract

Time series of soil moisture-related parameters provide important insights into the functioning of soil water systems. Analysis of patterns within such time series has been used in several studies. The objective of this work was to compare patterns in observed and simulated soil moisture contents to understand whether modelling leads to a substantial loss of information or complexity. The time series were observed at four plots in sandy soils within the USDA-ARS OPE3 experimental watershed, for a year; precipitation and evapotranspiration (ET) were measured and estimated, respectively, and used for soil water flow simulation with the HYDRUS-1D software. The information content measures are the metric entropy and the mean information gain, and complexity measures are the fluctuation complexity and the effective measure complexity. These measures were computed based on the binary encoding of soil moisture time series, and used probabilities of patterns, i.e. probabilities of joint or sequential appearance of symbol sequences. The information content of daily soil moisture time series was much smaller than that of rainfall data, and had higher complexity, indicating that soil worked essentially as an information filter. Information content and complexity decreased and increased with depth, respectively, demonstrating the increase in the information filtering action of soil. The information measures of simulated soil moisture content were close to those of the measurements, indicating the successful simulation of patterns in the data. The spatial variability of the information measures for simulated soil moisture content at all depths was less pronounced than the one of measured time series. Compared with precipitation and estimated ET, soil moisture time series had more structure and less randomness in this work. The information measures can provide useful complementary knowledge about model performance and patterns in observation and modelling results.

Citation Pan, F., Pachepsky, Y. A., Guber, A. K., & Hill, R. L. (2011) Information and complexity measures applied to observed and simulated soil moisture time series. Hydrol. Sci. J. 56(6), 1027–1039.     

Estimation of evapotranspiration and its components from an apple orchard in northwest China using sap flow and water balance methods

Field experiments were conducted to investigate the effects of leaf area index and soil moisture content on evapotranspiration and its components within an apple orchard in northwest China for 2 years. Evapotranspiration in the non‐rainfall period was estimated using two approaches: the soil water balance method based on tube‐type time‐domain reflection measurements, and sap flow plus micro‐lysimeter methods. The two methods were in good agreement, with differences usually less than 10%. The components of evapotranspiration varied with canopy development. During spring and autumn, soil evaporation was dominating as result of low leaf area index. In summer, plant transpiration became significant, with an average transpiration to evapotranspiration ratio of 0·87. The crop coefficient K_c showed a strong linear dependence on leaf area index. The water stress coefficient K_s was around 1·0 when soil moisture was above 23% and started to decrease linearly after that. This study demonstrates that prediction of evapotranspiration in apple orchards can be made using the Food and Agriculture Organization's crop coefficient method from commonly available meteorological data in the area. Copyright © 2006 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.     

Characteristics of soil moisture in permafrost observed in East Siberian taiga with stable isotopes of water

Soil moisture and its isotopic composition were observed at Spasskaya Pad experimental forest near Yakutsk, Russia, during summer in 1998, 1999, and 2000. The amount of soil water (plus ice) was estimated from volumetric soil water content obtained with time domain reflectometry. Soil moisture and its δ¹⁸O showed large interannual variation depending on the amount of summer rainfall. The soil water δ¹⁸O decreased with soil moisture during a dry summer (1998), indicating that ice meltwater from a deeper soil layer was transported upward. On the other hand, during a wet summer (1999), the δ¹⁸O of soil water increased due to percolation of summer rain with high δ¹⁸O values. Infiltration after spring snowmelt can be traced down to 15 cm by the increase in the amount of soil water and decrease in the δ¹⁸O because of the low δ¹⁸O of deposited snow. About half of the snow water equivalent (about 50 mm) recharged the surface soil. The pulse of the snow meltwater was, however, less important than the amount of summer rainfall for intra‐annual variation of soil moisture. Excess water at the time just before soil freezing, which is controlled by the amount of summer rainfall, was stored as ice during winter. This water storage stabilizes the rate of evapotranspiration. Soil water stored in the upper part of the active layer (surface to about 120 cm) can be a water source for transpiration in the following summer. On the other hand, once water was stored in the lower part of the active layer (deeper than about 120 cm), it would not be used by plants in the following summer, because the lower part of the active layer thaws in late summer after the plant growing season is over. Copyright © 2002 John Wiley & Sons, Ltd.     

Effects of surface characteristics on infiltration patterns in an arid shrub desert

Precipitation is often the sole source of water replenishment in arid and semi‐arid areas and, thus, plays a pertinent role in sustaining desert ecosystems. Revegetation over 40 years using mainly Artemisia ordosica and Caragana korshinskii at Shapotou Desert Experimental Research Station near Lanzhou, China, has established a dwarf‐shrub and microbiotic soil crust cover on the stabilized sand dunes. The redistribution of infiltrated moisture through percolation, root extraction, and evapotranspiration pathways was investigated. Three sets of time‐domain reflectometry (TDR) probes were inserted horizontally at 5, 10, 15, 20, 30 and 40 cm depths below the ground surface in a soil pit. The three sets of TDR probes were installed in dwarf‐shrub sites of A. ordosica and C. korshinskii community with and without a microbiotic soil crust cover, and an additional set was placed in a bare sand dune area that had neither vegetation nor a microbiotic soil crust present. Volumetric soil moisture content was recorded at hourly intervals and used in the assessment of infiltration for the different surface covers. Infiltration varied greatly, from 7·5 cm to more than 45 cm, depending upon rainfall quantity and soil surface conditions. In the shrub community area without microbiotic soil crust cover, infiltration increased due to preferential flow associated with root tunnels. The microbiotic soil crust cover had a significant negative influence on the infiltration for small rainfall events (～10 mm), restricting the infiltration depth to less than 20 cm and increasing soil moisture content just beneath the soil profile of 10 cm, whereas it was not as strong or clear for larger rainfall events (～60 mm). For small rainfall events, the wetting front depth for the three kinds of surface cover was as follows: shrub community without microbiotic soil crust > bare area > shrub community with microbiotic soil crust. In contrast, for large rainfall events, infiltration was similar in shrub communities with and without microbiotic soil crust cover, but significantly higher than measured in the bare area. Soil water extraction by roots associated with evapotranspiration restricted the wetting front penetration after 1 to 3 h of rainfall. Copyright © 2006 John Wiley & Sons, Ltd.     

Understanding hydrological trends by combining the Budyko hypothesis and a stochastic soil moisture model

Abstract

The actual evapotranspiration and runoff trends of five major basins in China from 1956 to 2000 are investigated by combining the Budyko hypothesis and a stochastic soil moisture model. Based on the equations of Choudhury and Porporato, the actual evapotranspiration trends and the runoff trends are attributed to changes in precipitation, potential evapotranspiration, rainfall depth and water storage capacity which depends on the soil water holding capacity and the root depth. It was found that the rainfall depth increased significantly in China during the past 50 years, especially in southern basins. Contributions from changes in the water storage capacity were significant in basins where land surface characteristics have changed substantially due to human activities. It was also observed that the actual evapotranspiration trends are more sensitive to precipitation trends in water-limited basins, but more sensitive to potential evapotranspiration trends in energy-limited basins.

Editor D. Koutsoyiannis; Associate editor A. Porporato     

Temporal dynamics of soil water balance components in a karst range in southeastern Spain: estimation of potential recharge

Abstract

This paper analyses the temporal dynamics of soil water balance components in a representative recharge area of the Sierra de Gádor (Almeria, southeastern Spain) in two hydrological years. Two approaches are used to estimate daily potential recharge (PR): Approach 1 based on deriving PR from the water balance as the difference between measurements of rainfall (P) and actual evapotranspiration (E) obtained by eddy covariance; and Approach 2 with PR obtained from the dynamic pattern of the soil moisture (θ) recorded at two depths in the site's thin soil (average 0.35 m thickess). For the hydrological year 2003/04, which was slightly drier than the 30-year average, E accounted for 64% of rainfall and occurred mainly in late spring and early summer. The PR estimated by Approach 1 was 181 ± 18 mm year^-1 (36% of rainfall), suggesting an effective groundwater recharge in the study area. In the unusually dry hydrological year 2004/05, E was about 215 mm year^-1, close to the annual rainfall input, and allowing very little (8 ± 12 mm year^-1) PR according to Approach 1. Estimation of PR based on Approach 2 resulted in PR rates lower than those found by Approach 1, because Approach 2 does not take into account the recharge that occurs through preferential flow pathways (cracks, joints and fissures) which were not monitored with the θ probes. Moreover, using Approach 2, the PR estimates differed widely depending on the time scale considered: with daily mean θ data, PR estimation was lower, especially in late spring, while θ data at 30 min resolution yielded a more reliable prediction of the fraction of total PR resulting from the downward movement of soil water by gravity.

Citation Cantón, Y., Villagarcía, L., Moro, M. J., Serrano-Ortíz, P., Were, A., Alcalá, F. J., Kowalski, A. S., Solé-Benet, A., Lázaro, R. & Domingo, F. (2010) Temporal dynamics of soil water balance components in a karst range in southeastern Spain: estimation of potential recharge. Hydrol. Sci. J. 55(5), 737–753.     

Modelling the spatial variation in soil moisture at the landscape scale: an application to five areas of ecological interest in the UK

A simple conceptual hydrological model that explicitly includes the lateral movement of soil water and operates efficiently at the landscape scale is outlined. It is applied to five areas of ecological interest in the UK to provide distributed mean monthly soil moisture on a 50 m grid. As the model's driving variables—daily rainfall and potential evapotranspiration—are assumed constant over each of the tracts of land, the variability in soil moisture is due to different soil types and to topographic effects. Box plots of the mean monthly simulated soil moisture clearly show the spread of values occasioned by modelling the lateral water movement down the hillslope. The general magnitude of the results are compared with published data wherever possible and there is some discussion of the form of the curve used in the model to describe the attenuation of evapotranspiration with decreasing soil moisture. Copyright © 2000 John Wiley & Sons, Ltd.     

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

Abstract

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

Editor D. Koutsoyiannis; Associate editor D. Hughes

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

Estimating water balance components in irrigated agriculture using a combined approach of soil moisture and energy balance monitoring,and numerical modelling

Information on water balance components such as evapotranspiration and groundwater recharge are crucial for water management. Due to differences in physical conditions, but also due to limited budgets, there is not one universal best practice, but a wide range of different methods with specific advantages and disadvantages. In this study, we propose an approach to quantify actual evapotranspiration, groundwater recharge and water inflow, i.e. precipitation and irrigation, that considers the specific conditions of irrigated agriculture in warm, arid environments. This approach does not require direct measurements of precipitation or irrigation quantities and is therefore suitable for sites with an uncertain data basis. For this purpose, we combine soil moisture and energy balance monitoring, remote sensing data analysis and numerical modelling using Hydrus. Energy balance data and routine weather data serve to estimate ET₀. Surface reflectance data from satellite images (Sentinel-2) are used to derive leaf area indices, which help to partition ET₀ into energy limited evaporation and transpiration. Subsequently, first approximations of water inflow are derived based on observed soil moisture changes. These inflow estimates are used in a series of forward simulations that produce initial estimates of drainage and ET_act, which in turn help improve the estimate of water inflow. Finally, the improved inflow estimates are incorporated into the model and then a parameter optimization is performed using the observed soil moisture as the reference figure. Forward simulations with calibrated soil parameters result in final estimates for ET_act and groundwater recharge. The presented method is applied to an agricultural test site with a crop rotation of cotton and wheat in Punjab, Pakistan. The final model results, with an RMSE of 2.2% in volumetric water content, suggest a cumulative ET_act and groundwater recharge of 769 and 297 mm over a period of 281 days, respectively. The total estimated water inflow accounts for 946 mm, of which 77% originates from irrigation.     

On the estimation of antecedent wetness conditions in rainfall–runoff modelling

The analysis of the physical processes involved in a conceptual model of soil water content balance is addressed with the objective of its application as a component of rainfall–runoff modelling. The model uses routinely measured meteorological variables (rainfall and air temperature) and incorporates a limited number of significant parameters. Its performance in estimating the soil moisture temporal pattern was tested through local measurements of volumetric water content carried out continuously on an experimental plot located in central Italy. The analysis was carried out for different periods in order to test both the representation of infiltration at the short time‐scale and drainage and evapotranspiration processes at the long time‐scale. A robust conceptual model was identified that incorporated the Green–Ampt approach for infiltration and a gravity‐driven approximation for drainage. A sensitivity analysis was performed for the selected model to assess the model robustness and to identify the more significant parameters involved in the principal processes that control the soil moisture temporal pattern. The usefulness of the selected model was tested for the estimation of the initial wetness conditions for rainfall–runoff modelling at the catchment scale. Specifically, the runoff characteristics (runoff depth and peak discharge) were found to be dependent on the pre‐event surface soil moisture. Both observed values and those estimated by the model gave good results. On the contrary, with the antecedent wetness conditions furnished by two versions of the antecedent precipitation index (API), large errors were obtained. Copyright © 2007 John Wiley & Sons, Ltd.     

Water balance of a tropical woodland ecosystem, Northern Australia: a combination of micro-meteorological, soil physical and groundwater chemical approaches

A combination of micro-meteorological, soil physical and groundwater chemical methods enabled the water balance of a tropical eucalypt savanna ecosystem in Northern Australia to be estimated. Heat pulse and eddy correlation were used to determine overstory and total evapotranspiration, respectively. Measurements of soil water content, matric suction and water table variations were used to determine changes in soil moisture storage throughout the year. Groundwater dating with chlorofluorocarbons was used to estimate net groundwater recharge rates, and stream gauging was used to determine surface runoff. The wet season rainfall of 1585 mm is distributed as: evapotranspiration 810 mm, surface runoff (and shallow subsurface flow) into the river 410 mm, groundwater recharge 200 mm and increase in soil store 165 mm. Of the groundwater recharge, 160 mm enters the stream as baseflow in the wet season, 20 mm enters as baseflow in the dry season, and the balance (20 mm) is distributed to and used by minor vegetation types within the catchment or discharges to the sea. In the dry season, an evapotranspiration of 300 mm comprises 135 mm rainfall and 165 mm from the soil store. Because of the inherent errors of the different techniques, the water balance surplus (estimated at 20 mm) cannot be clearly distinguished from zero. It may also be as much as 140 mm. To our knowledge, this is the first time that such diverse methods have been combined to estimate all components of a catchment's water balance.     

Sensitivity of soil moisture field evolution to rainfall forcing

In this paper the temporal behaviour of soil moisture is modelled and statistically characterized by use of the zero‐dimensional model for soil moisture dynamics and the rectangular pulses Poisson process model for rainfall forcing. The mean, covariance and spectral density function of soil moisture (both instantaneous and locally averaged cases) are analytically derived to evaluate its sensitivity to the model parameters. Finally, the probability density function of soil moisture is derived to evaluate the effect of rainfall forcing. All the model parameters used have been tuned to the Monsoon '90 data. Results can be summarized as follows. (1) Only the soil moisture model parameters (η and nZ_r) are found to affect the autocorrelation function in a distinguishable manner. On the other hand, both the rainfall model parameter (θ) and the effective soil depth (nZ_r) are found to be of impact to the soil moisture spectrum. However, as the smoothing (or damping) effect of soil is so dominant, about ±20% variation of one parameter seems not to affect significantly the second‐order statistics of soil moisture. (2) More difference can be found by applying a longer averaging time, which is found to obviously decrease the variance but increase the correlation even though no overlapping between neighbouring soil moisture data was allowed. (3) Among rainfall model parameters, the arrival rate (λ) was found to be most important for the soil moisture evolution. When increasing the arrival rate of rainfall, the histogram of soil moisture shifts its peak to a certain value as well as becomes more concentrated around the peak. However, by decreasing the arrival rate of rainfall, a much smaller (almost to zero) mean value of soil moisture was estimated, even though the total volume of rainfall remained constant. This indicates that desertification may take place without decreasing the total volume of rainfall. Copyright © 2005 John Wiley & Sons, Ltd.     

The role of rainfall temporal and spatial averaging in seasonal simulations of the terrestrial water balance

The partitioning of rainfall into surface runoff and infiltration influences many other aspects of the hydrologic cycle including evapotranspiration, deep drainage and soil moisture. This partitioning is an instantaneous non-linear process that is strongly dependent on rainfall rate, soil moisture and soil hydraulic properties. Though all rainfall datasets involve some degree of spatial or temporal averaging, it is not understood how this averaging affects simulated partitioning and the land surface water balance across a wide range of soil and climate types. We used a one-dimensional physics-based model of the near-surface unsaturated zone to compare the effects of different rainfall discretization (5-min point-scale; hourly point-scale; hourly 0.125° gridded) on the simulated partitioning of rainfall for many locations across the United States. Coarser temporal resolution rainfall data underpredicted seasonal surface runoff for all soil types except those with very high infiltration capacities (i.e., sand, loamy sand). Soils with intermediate infiltration capacities (i.e., loam, sandy loam) were the most affected, with less than half of the expected surface runoff produced in most soil types when the gridded rainfall dataset was used as input. The impact of averaging on the water balance was less extreme but non-negligible, with the hourly point-scale predictions exhibiting median evapotranspiration, drainage and soil moisture values within 10% of those predicted using the higher resolution 5-min rainfall. Water balance impacts were greater using the gridded hourly dataset, with average underpredictions of ET up to 27% in fine-grained soils. The results suggest that “hyperresolution” modelling at continental to global scales may produce inaccurate predictions if there is not parallel effort to produce higher resolution precipitation inputs or sub-grid precipitation parameterizations.     

Hydrological sciences—society's needs

Abstract

The hydrological regime of a mountainous catchment, in this instance the Mesochora catchment in Central Greece, was simulated for altered climates resulting when using the Goddard Institute for Space Studies (GISS) model for carbon dioxide doubling. The catchment snow water equivalent was predicted on the basis of the snow accumulation and ablation model of the US National Weather Service River Forecast System (NWSRFS), while the catchment runoff, as well as actual evapotranspiration and soil moisture storages, were simulated through application of the soil moisture accounting model of NWSRFS. Two scenarios of monthly climate change were drawn from the GISS model, one associated with temperature and precipitation changes, while the other referred to temperature changes alone. A third hypothetical scenario with temperature and precipitation changes similar to those corresponding to the mean monthly GISS scenarios was used to test the sensitivity of the monthly climate change of the hypothetical case on catchment hydrology. All three scenarios projected decreases in average snow accumulations and in spring and summer runoff and soil moisture, as well as increases in winter runoff and soil moisture storage and spring evapotranspiration.     

Comparative application of two mathematical models to predict sedimentation in Yermasoyia Reservoir,Cyprus

The Yermasoyia Reservoir is located northeast of the town of Limassol, Cyprus. The storage capacity of the reservoir is 13·6 × 10⁶ m³. The basin area of the Yermasoyia River, which feeds the reservoir, totals 122·5 km². This study aims to estimate the mean annual deposition amount in the reservoir, which originates from the corresponding basin. For the estimate of the mean annual sediment inflow into the reservoir, two mathematical models are used alternatively. Each model consists of three submodels: a rainfall‐runoff submodel, a soil erosion submodel and a sediment transport submodel for streams. In the first model, the potential evapotranspiration is estimated for the rainfall‐runoff submodel, and the soil erosion submodel of Schmidt and the sediment transport submodel of Yang are used. In the second model, the actual evapotranspiration is estimated for the rainfall‐runoff submodel, and the soil erosion submodel of Poesen and the sediment transport submodel of Van Rijn are used. The deposition amount in the reservoir is estimated by means of the diagram of Brune, which delivers the trap efficiency of the reservoir. Daily rainfall data from three rainfall stations, and daily values of air temperature, relative air humidity and sunlight hours from a meteorological station for four years (1986–89) were available. The computed annual runoff volumes and mean annual soil erosion rate are compared with the respective measurement data. Copyright © 2006 John Wiley & Sons, Ltd.     

Sensitivity of hydrological variables to climate change in the Haihe River basin,China

Using the defined sensitivity index, the sensitivity of streamflow, evapotranspiration and soil moisture to climate change was investigated in four catchments in the Haihe River basin. Climate change contained three parts: annual precipitation and temperature change and the change of the percentage of precipitation in the flood season (P_f). With satisfying monthly streamflow simulation using the variable infiltration capacity model, the sensitivity was estimated by the change of simulated hydrological variables with hypothetical climatic scenarios and observed climatic data. The results indicated that (i) the sensitivity of streamflow would increase as precipitation or P_f increased but would decrease as temperature increased; (ii) the sensitivity of evapotranspiration and soil moisture would decrease as precipitation or temperature increased, but it to P_f varied in different catchments; and (iii) hydrological variables were more sensitive to precipitation, followed by P_f, and then temperature. The nonlinear response of streamflow, evapotranspiration and soil moisture to climate change could provide a reference for water resources planning and management under future climate change scenarios in the Haihe River basin. Copyright © 2011 John Wiley & Sons, Ltd.     

A method for estimating soil moisture storage in regions under water stress and storage depletion: a case study of Hai River Basin,North China

Soil moisture is a consideration for soil conservation, agricultural production and climate modelling. This article presents a simple method for estimating soil moisture storage under water stress and storage depletion conditions. The method is driven by the common agro‐hydrologic variables of precipitation (PPT), irrigation (IRR) and evapotranspiration (ET). The proposed method is successfully tested for the 152 000 km² floodplain region of Hai River Basin using 48 consecutive months (2003–2006) of data. Soil moisture data from global land data assimilation system/Noah land surface model are validated with ground‐truth data from 102 soil moisture monitoring sites. The validated soil moisture is used in combination with in situ groundwater data to quantify total water storage change (TWSC) in the region. The estimated storage change is in turn compared with gravity recovery and climate experiment‐derived TWSC for the study area. The soil moisture and TWSC terms show favourable agreements, with discrepancies of < 10% on the average. While there is no consistent seasonal trend in soil moisture, TWSC shows a strong seasonality. It is low in spring and high in summer. This trend corresponds with the IRR–PPT season in the study area. Change in groundwater and total water storage indicates storage depletion in the basin. Storage depletion in the region is driven mainly by groundwater IRR and ET loss. Despite the low PPT and high ET, there is narrowing seasonal trend in soil moisture. This is achieved at the expense of groundwater storage. IRR pumping has induced extensive groundwater depletion in the basin. It is therefore vital to develop cultivation strategies that aim at limiting IRR pumping and ET loss. Water management practices that not only reduce waste but also ensure high productivity and ecological sustainability could also mitigate storage depletion in the region. These measures could reduce further not only the seasonal trend in soil moisture but also that in groundwater storage. Copyright © 2011 John Wiley & Sons, Ltd.