Improving the estimation of hydrothermal state variables in the active layer of frozen ground by assimilating <Emphasis Type="Italic">in situ</Emphasis> observations and SSM/I data

The active layer of frozen ground data assimilation system adopts the SHAW (Simulteneous Heat and Water) model as the model operator. It employs an ensemble kalman filter to fuse state variables predicted by the SHAW model with in situ observation and the SSM/I 19 GHz brightness temperature for the purpose of optimizing model hydrothermal state variables. When there is little water movement in the frozen soil during the winter season, the unfrozen water content depends primarily on soil temperature. Thus, soil temperature is the crucial state variable to be improved. In contrast, soil moisture is heavily influenced by precipitation during the summer season. The simulation accuracy of soil moisture has a strong and direct impact on the soil temperature. In this case, the crucial state variable to be improved is soil moisture. One-dimensional assimilation experiments that have been carried out at AMDO station show that land data assimilation method can improve the estimation of hydrothermal state variables in the soil by fusing model information and observation information. The reasonable model error covariance matrix plays a key role in transferring the optimized surface state information to the deep soil, and it provides improved estimations of whole soil state profiles. After assimilating the 4-cm soil temperature by in situ observation, the soil temperature RMSE (Root Mean Square Error) of each soil layer decreased by 0.96°C on average relative to the SHAW simulation. After assimilating the 4-cm soil moisture in situ observation, the soil moisture RMSE of each soil layer decreased by 0.020 m³·m⁻³. When assimilating the SSM/I 19 GHz brightness temperature, the soil temperature RMSE of each soil layer during the winter decreased by 0.76°C, while the soil moisture RMSE of each soil layer during the summer decreased by 0.018 m³·m⁻³.     

Validation of ERS scatterometer‐derived soil moisture data in the central part of the Duero Basin,Spain

The objective of this study was to validate the soil moisture data derived from coarse‐resolution active microwave data (50 km) from the ERS scatterometer. The retrieval technique is based on a change detection method coupled with a data‐based modelling approach to account for seasonal vegetation dynamics. The technique is able to derive information about the soil moisture content corresponding to the degree of saturation of the topmost soil layer (∼5 cm). To estimate profile soil moisture contents down to 100 cm depth from the scatterometer data, a simple two‐layer water balance model is used, which generates a red noise‐like soil moisture spectrum. The retrieval technique had been successfully applied in the Ukraine in a previous study. In this paper, the performance of the model in a semi‐arid Mediterranean environment characterized by low annual precipitation (400 mm), hot dry summers and sandy soils is investigated. To this end, field measurements from the REMEDHUS soil moisture station network in the semi‐arid parts of the Duero Basin (Spain) were used. The results reveal a significant coefficient of determination (R² = 0·75) for the averaged 0–100 cm soil moisture profile and a root mean square error (RMSE) of 2·2 vol%. The spatial arrangement of the REMEDHUS soil moisture stations also allowed us to study the influence of the small‐scale variability of soil moisture within the ERS scatterometer footprint. The results show that the small‐scale variability in the study area is modest and can be explained in terms of texture fraction distribution in the soil profiles. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulated water budget of a small forested watershed in the continental/maritime hydroclimatic region of the United States

Annual streamflows have decreased across mountain watersheds in the Pacific Northwest of the United States over the last ~70 years; however, in some watersheds, observed annual flows have increased. Physically based models are useful tools to reveal the combined effects of climate and vegetation on long‐term water balances by explicitly simulating the internal watershed hydrological fluxes that affect discharge. We used the physically based Simultaneous Heat and Water (SHAW) model to simulate the inter‐annual hydrological dynamics of a 4 km² watershed in northern Idaho. The model simulates seasonal and annual water balance components including evaporation, transpiration, storage changes, deep drainage, and trends in streamflow. Independent measurements were used to parameterize the model, including forest transpiration, stomatal feedback to vapour pressure, forest properties (height, leaf area index, and biomass), soil properties, soil moisture, snow depth, and snow water equivalent. No calibrations were applied to fit the simulated streamflow to observations. The model reasonably simulated the annual runoff variations during the evaluation period from water year 2004 to 2009, which verified the ability of SHAW to simulate the water budget in this small watershed. The simulations indicated that inter‐annual variations in streamflow were driven by variations in precipitation and soil water storage. One key parameterization issue was leaf area index, which strongly influenced interception across the catchment. This approach appears promising to help elucidate the mechanisms responsible for hydrological trends and variations resulting from climate and vegetation changes on small watersheds in the region. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling soil moisture in a high‐latitude landscape using LiDAR and soil data

Soil moisture has a pronounced effect on earth surface processes. Global soil moisture is strongly driven by climate, whereas at finer scales, the role of non‐climatic drivers becomes more important. We provide insights into the significance of soil and land surface properties in landscape‐scale soil moisture variation by utilizing high‐resolution light detection and ranging (LiDAR) data and extensive field investigations. The data consist of 1200 study plots located in a high‐latitude landscape of mountain tundra in north‐western Finland. We measured the plots three times during growing season 2016 with a hand‐held time‐domain reflectometry sensor. To model soil moisture and its temporal variation, we used four statistical modelling methods: generalized linear models, generalized additive models, boosted regression trees, and random forests. The model fit of the soil moisture models were R² = 0.60 and root mean square error (RMSE) 8.04 VWC% on average, while the temporal variation models showed a lower fit of R² = 0.25 and RMSE 13.11 CV%. The predictive performances for the former were R² = 0.47 and RMSE 9.34 VWC%, and for the latter R² = 0.01 and RMSE 15.29 CV%. Results were similar across the modelling methods, demonstrating a consistent pattern. Soil moisture and its temporal variation showed strong heterogeneity over short distances; therefore, soil moisture modelling benefits from high‐resolution predictors, such as LiDAR based variables. In the soil moisture models, the strongest predictor was SAGA (System for Automated Geoscientific Analyses) wetness index (SWI), based on a 1 m² digital terrain model derived from LiDAR data, which outperformed soil predictors. Thus, our study supports the use of LiDAR based SWI in explaining fine‐scale soil moisture variation. In the temporal variation models, the strongest predictor was the field‐quantified organic layer depth variable. Our results show that spatial soil moisture predictions can be based on soil and land surface properties, yet the temporal models require further investigation. Copyright © 2017 John Wiley & Sons, Ltd.     

Assessing the impact of interannual variability of precipitation and potential evaporation on evapotranspiration

The impact of interannual variability of precipitation and potential evaporation on the long-term mean annual evapotranspiration as well as on the interannual variability of evapotranspiration is studied using a stochastic soil moisture model within the Budyko framework. Results indicate that given the same long-term mean annual precipitation and potential evaporation, including interannual variability of precipitation and potential evaporation reduces the long-term mean annual evapotranspiration. This reduction effect is mostly prominent when the dryness index (i.e., the ratio of potential evaporation to precipitation) is within the range from 0.5 to 2. The maximum reductions in the evaporation ratio (i.e., the ratio of evapotranspiration to precipitation) can reach 8–10% for a range of coefficient of variation (CV) values for precipitation and potential evaporation. The relations between the maximum reductions and the CV values of precipitation and potential evaporation follow power laws. Hence the larger the interannual variability of precipitation and potential evaporation becomes, the larger the reductions in the evaporation ratio will be. The inclusion of interannual variability of precipitation and potential evaporation also increases the interannual variability of evapotranspiration. It is found that the interannual variability of daily rainfall depth and that of the frequency of daily rainfall events have quantitatively different impacts on the interannual variability of evapotranspiration; and they also interact differently with the interannual variability of potential evaporation. The results presented in this study demonstrate the importance of understanding the role of interannual variability of precipitation and potential evaporation in land surface hydrology under a warming climate.     

Estimation of soil hydraulic properties and their uncertainty: comparison between laboratory and field experiment

Often the soil hydraulic parameters are obtained by the inversion of measured data (e.g. soil moisture, pressure head, and cumulative infiltration, etc.). However, the inverse problem in unsaturated zone is ill‐posed due to various reasons, and hence the parameters become non‐unique. The presence of multiple soil layers brings the additional complexities in the inverse modelling. The generalized likelihood uncertainty estimate (GLUE) is a useful approach to estimate the parameters and their uncertainty when dealing with soil moisture dynamics which is a highly non‐linear problem. Because the estimated parameters depend on the modelling scale, inverse modelling carried out on laboratory data and field data may provide independent estimates. The objective of this paper is to compare the parameters and their uncertainty estimated through experiments in the laboratory and in the field and to assess which of the soil hydraulic parameters are independent of the experiment. The first two layers in the field site are characterized by Loamy sand and Loamy. The mean soil moisture and pressure head at three depths are measured with an interval of half hour for a period of 1 week using the evaporation method for the laboratory experiment, whereas soil moisture at three different depths (60, 110, and 200 cm) is measured with an interval of 1 h for 2 years for the field experiment. A one‐dimensional soil moisture model on the basis of the finite difference method was used. The calibration and validation are approximately for 1 year each. The model performance was found to be good with root mean square error (RMSE) varying from 2 to 4 cm³ cm^?3. It is found from the two experiments that mean and uncertainty in the saturated soil moisture (θ_s) and shape parameter (n) of van Genuchten equations are similar for both the soil types. Copyright © 2010 John Wiley & Sons, Ltd.     

Particle Filter-based assimilation algorithms for improved estimation of root-zone soil moisture under dynamic vegetation conditions

In this study, we implement Particle Filter (PF)-based assimilation algorithms to improve root-zone soil moisture (RZSM) estimates from a coupled SVAT-vegetation model during a growing season of sweet corn in North Central Florida. The results from four different PF algorithms were compared with those from the Ensemble Kalman Filter (EnKF) when near-surface soil moisture was assimilated every 3 days using both synthetic and field observations. In the synthetic case, the PF algorithm with the best performance used residual resampling of the states and obtained resampled parameters from a uniform distribution and provided reductions of 76% in root mean square error (RMSE) over the openloop estimates. The EnKF provided the RZSM and parameter estimates that were closer to the truth than the PF with an 84% reduction in RMSE. When field observations were assimilated, the PF algorithm that maintained maximum parameter diversity offered the largest reduction of 16% in root mean square difference (RMSD) over the openloop estimates. Minimal differences were observed in the overall performance of the EnKF and PF using field observations since errors in model physics affected both the filters in a similar manner, with maximum reductions in RMSD compared to the openloop during the mid and reproductive stages.     

Effect of simultaneous state–parameter estimation and forcing uncertainties on root-zone soil moisture for dynamic vegetation using EnKF

In this study, an EnKF-based assimilation algorithm was implemented to estimate root-zone soil moisture (RZSM) using the coupled LSP–DSSAT model during a growing season of corn. Experiments using both synthetic and field observations were conducted to understand effects of simultaneous state–parameter estimation, spatial and temporal update frequency, and forcing uncertainties on RZSM estimates. Estimating the state–parameters simultaneously with every 3-day assimilation of volumetric soil moisture (VSM) observations at 5 depths lowered the average standard deviation (ASD) and the root mean square error (RMSE) for RZSM by approximately 1.77% VSM (78%) and 2.18% VSM (93%), respectively, compared to the open-loop ASD where as estimating only states lowered the ASD by approximately 1.26% VSM (56%) and the RMSE by 1.66% VSM (71%). The synthetic case obtained RZSM estimates closer to the observations than the MicroWEX-2 case, particularly after precipitation/irrigation events. The differences in EnKF performance between MicroWEX-2 and synthetic observations may indicate other sources of errors in addition to those in parameters and forcings, such as errors in model biophysics.     

Modelling spatio-temporal soil moisture dynamics in mountain tundra

Soil moisture has a fundamental influence on the processes and functions of tundra ecosystems. Yet, the local dynamics of soil moisture are often ignored, due to the lack of fine resolution, spatially extensive data. In this study, we modelled soil moisture with two mechanistic models, SpaFHy (a catchment-scale hydrological model) and JSBACH (a global land surface model), and examined the results in comparison with extensive growing-season field measurements over a mountain tundra area in northwestern Finland. Our results show that soil moisture varies considerably in the study area and this variation creates a mosaic of moisture conditions, ranging from dry ridges (growing season average 12 VWC%, Volumetric Water Content) to water-logged mires (65 VWC%). The models, particularly SpaFHy, simulated temporal soil moisture dynamics reasonably well in parts of the landscape, but both underestimated the range of variation spatially and temporally. Soil properties and topography were important drivers of spatial variation in soil moisture dynamics. By testing the applicability of two mechanistic models to predict fine-scale spatial and temporal variability in soil moisture, this study paves the way towards understanding the functioning of tundra ecosystems under climate change.     

Evapotranspiration models of different complexity for multiple land cover types

A comparison between half‐hourly and daily measured and computed evapotranspiration (ET) using three models of different complexity, namely, the Priestley–Taylor (P‐T), the reference Penman–Monteith (P‐M) and the Common Land Model (CLM), was conducted using three AmeriFlux sites under different land cover and climate conditions (i.e. arid grassland, temperate forest and subhumid cropland). Using the reference P‐M model with a semiempirical soil moisture function to adjust for water‐limiting conditions yielded ET estimates in reasonable agreement with the observations [root mean square error (RMSE) of 64–87 W m^?2 for half‐hourly and RMSE of 0.5–1.9 mm day^?1 for daily] and similar to the complex Common Land Model (RMSE of 60–94 W m^?2 for half‐hourly and RMSE of 0.4–2.1 mm day^?1 for daily) at the grassland and cropland sites. However, the semiempirical soil moisture function was not applicable particularly for the P‐T model at the forest site, suggesting that adjustments to key model variables may be required when applied to diverse land covers. On the other hand, under certain land cover/environmental conditions, the use of microwave‐derived soil moisture information was found to be a reliable metric of regional moisture conditions to adjust simple ET models for water‐limited cases. Further studies are needed to evaluate the utility of the simplified methods for different landscapes. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical investigation of the impact of uncertainties in satellite rainfall estimation and land surface model parameters on simulation of soil moisture

This study aims at evaluating the uncertainty in the prediction of soil moisture (1D, vertical column) from an offline land surface model (LSM) forced by hydro-meteorological and radiation data. We focus on two types of uncertainty: an input error due to satellite rainfall retrieval uncertainty, and, LSM soil-parametric error. The study is facilitated by in situ and remotely sensed data-driven (precipitation, radiation, soil moisture) simulation experiments comprising a LSM and stochastic models for error characterization. The parametric uncertainty is represented by the generalized likelihood uncertainty estimation (GLUE) technique, which models the parameter non-uniqueness against direct observations. Half-hourly infra-red (IR) sensor retrievals were used as satellite rainfall estimates. The IR rain retrieval uncertainty is characterized on the basis of a satellite rainfall error model (SREM). The combined uncertainty (i.e., SREM + GLUE) is compared with the partial assessment of uncertainty. It is found that precipitation (IR) error alone may explain moderate to low proportion of the soil moisture simulation uncertainty, depending on the level of model accuracy—50–60% for high model accuracy, and 20–30% for low model accuracy. Comparisons on the basis of two different sites also yielded an increase (50–100%) in soil moisture prediction uncertainty for the more vegetated site. This study exemplified the need for detailed investigations of the rainfall retrieval-modeling parameter error interaction within a comprehensive space-time stochastic framework for achieving optimal integration of satellite rain retrievals in land data assimilation systems.     

Soil moisture dynamics across landscape types in an arid inland river basin of Northwest China

A combination of field measurements, continuous monitoring and numerical modelling was used to evaluate soil moisture regimes at four sites across a landscape gradient of the Heihe River Basin. Recorded data of precipitation, irrigation and floods were used to build the model, and an optimization technique was employed to calibrate the parameters. Based on the optimized parameters and estimates of future scenarios, the modelling structure was employed to predict the changes in the growing season soil moisture regimes due to climate change and intensive management. The results suggest that the upper‐reach Yeniugou and Xishui sites will become wetter because of alterations in the precipitation regime, and this trend could be strengthened by the expected amplified interannual variability. Precipitation features at middle‐reach Linze and lower‐reach Ejina, however, are not expected to change in the future. We assumed that a more water‐saving irrigation system will replace the current traditional one, and hence, the soil moisture probability density function at the Linze site would tend to be narrowed to ranges around either the wilting point or the point of incipient water stress, depending on how the intervention point and target level are settled. Ejina is expected to experience the most extreme ecological conversion effects in the future, and soil moisture would be more frequently recharged by water delivery. However, the soil moisture regime would not change much because of the poor water‐holding capacity and intensive evaporation. The revealed patterns and predicted shifts in soil moisture dynamics could provide a useful reference for identifying robust long‐term water resource management strategies for the Heihe River Basin. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimating growing-season root zone soil moisture from vegetation index-based evapotranspiration fraction and soil properties in the Northwest Mountain region,USA

A soil moisture retrieval method is proposed, in the absence of ground-based auxiliary measurements, by deriving the soil moisture content relationship from the satellite vegetation index-based evapotranspiration fraction and soil moisture physical properties of a soil type. A temperature–vegetation dryness index threshold value is also proposed to identify water bodies and underlying saturated areas. Verification of the retrieved growing season soil moisture was performed by comparative analysis of soil moisture obtained by observed conventional in situ point measurements at the 239-km² Reynolds Creek Experimental Watershed, Idaho, USA (2006–2009), and at the US Climate Reference Network (USCRN) soil moisture measurement sites in Sundance, Wyoming (2012–2015), and Lewistown, Montana (2014–2015). The proposed method best represented the effective root zone soil moisture condition, at a depth between 50 and 100 cm, with an overall average R² value of 0.72 and average root mean square error (RMSE) of 0.042.     

Evaluation of the PESERA model in two contrasting environments

The performance of the Pan‐European Soil Erosion Risk Assessment (PESERA) model was evaluated by comparison with existing soil erosion data collected in plots under different land uses and climate conditions in Europe. In order to identify the most important sources of error, the PESERA model was evaluated by comparing model output with measured values as well as by assessing the effect of the various model components on prediction accuracy through a multistep approach. First, the performance of the hydrological and erosion components of PESERA was evaluated separately by comparing both runoff and soil loss predictions with measured values. In order to assess the performance of the vegetation growth component of PESERA, the predictions of the model based on observed values of vegetation ground cover were also compared with predictions based on the simulated vegetation cover values. Finally, in order to evaluate the sediment transport model, predicted monthly erosion rates were also calculated using observed values of runoff and vegetation cover instead of simulated values. Moreover, in order to investigate the capability of PESERA to reproduce seasonal trends, the observed and simulated monthly runoff and erosion values were aggregated at different temporal scale and we investigated at what extend the model prediction error could be reduced by output aggregation. PESERA showed promise to predict annual average spatial variability quite well. In its present form, short‐term temporal variations are not well captured probably due to various reasons. The multistep approach showed that this is not only due to unrealistic simulation of cover and runoff, being erosion prediction also an important source of error. Although variability between the investigated land uses and climate conditions is well captured, absolute rates are strongly underestimated. A calibration procedure, focused on a soil erodibility factor, is proposed to reduce the significant underestimation of soil erosion rates. Copyright © 2009 John Wiley & Sons, Ltd.     

A soil moisture assimilation scheme based on the ensemble Kalman filter using microwave brightness temperature

This study presents a soil moisture assimilation scheme, which could assimilate microwave brightness temperature directly, based on the ensemble Kalman filter and the shuffled complex evolution method (SCE-UA). It uses the soil water model of the land surface model CLM3.0 as the forecast operator, and a radiative transfer model (RTM) as the observation operator in the assimilation system. The assimilation scheme is implemented in two phases: the parameter calibration phase and the pure soil moisture assimilation phase. The vegetation optical thickness and surface roughness parameters in the RTM are calibrated by SCE-UA method and the optimal parameters are used as the final model parameters of the observation operator in the assimilation phase. The ideal experiments with synthetic data indicate that this scheme could significantly improve the simulation of soil moisture at the surface layer. Furthermore, the estimation of soil moisture in the deeper layers could also be improved to a certain extent. The real assimilation experiments with AMSR-E brightness temperature at 10.65 GHz (vertical polarization) show that the root mean square error (RMSE) of soil moisture in the top layer (0–10 cm) by assimilation is 0.03355 m³ · m⁻³, which is reduced by 33.6% compared with that by simulation (0.05052 m³ · m⁻³). The mean RMSE by assimilation for the deeper layers (10–50 cm) is also reduced by 20.9%. All these experiments demonstrate the reasonability of the assimilation scheme developed in this study.     

Determination of bare surface soil moisture from combined temporal evolution of land surface temperature and net surface shortwave radiation

Land surface soil moisture (SSM) is an important variable for hydrological, ecological, and meteorological applications. A multi‐linear model has recently been proposed to determine the SSM content from the combined diurnal evolution of both land surface temperature (LST) and net surface shortwave radiation (NSSR) with the parameters TN (the LST mid‐morning rising rate divided by the NSSR rising rate during the same period) and t_d (the time of daily maximum temperature). However, in addition to the problem that all the coefficients of the multi‐linear model depend on the atmospheric conditions, the model also suffers from the problems of the nonlinearity of TN as a function of the SSM content and the uncertainty of determining the t_d from the diurnal evolution of the LST. To address these problems, a modified multi‐linear model was developed using the logarithm of TN and normalizing t_d by the mid‐morning temperature difference instead of using the TN and t_d. Except for the constant term, the coefficients of all other variables in the modified multi‐linear model proved to be independent of the atmospheric conditions. Using the relevant simulation data, results from the modified multi‐linear model show that the SSM content can be determined with a root mean square error (RMSE) of 0.030m³/m³, provided that the constant term is known or estimated day to day. The validation of the model was conducted using the field measurements at the Langfang site in 2008 in China. A higher correlation is achieved (coefficient of determination: R² = 0.624, RMSE = 0.107m³/m³) between the measured SSM content and the SSM content estimated using the modified multi‐linear model with the coefficients determined from the simulation data. Another experiment is also conducted to estimate the SSM content using the modified model with the constant term calibrated each day by one‐spot measurements at the site. The estimation result has a relatively larger error (RMSE = 0.125m³/m³). Additionally, the uncertainty of the determination of the coefficients is analysed using the field measurements, and the results indicate that the SSM content obtained using the modified model accurately characterizes the surface soil moisture condition. Copyright © 2012 John Wiley & Sons, Ltd.     

Improved SAL method and its application to verifying regional soil moisture forecasting

The regional verification of soil moisture is a vital step in evaluating and improving numerical model performance and utilizing forecast results. Currently, even with improved spatial and temporal resolutions of numerical model, verification methods for soil moisture data still rely on the traditional intensity verification parameters, such as mean error (ME) and root-mean-squared error (RMSE). Those methods provide only incomplete and sometimes inaccurate messages and thus hinder a proper evaluation of a forecast model. The SAL method is an object-based regional verification method with respect to precipitation forecasts. Based on the SAL method, a novel object-based method (SAL-DN) is proposed here, which can be used to test regional soil moisture. Both the ideal experiment and real experiment show that the SAL-DN method can reveal the differences between the observed and forecast soil moisture in three aspects: structure, amplitude, and location, and the results can reflect the actual situation. Furthermore, compared with the SAL method, the SAL-DN method is also capable of verifying physical quantities with high-value and low-value centers like temperature. Therefore, the SAL-DN method enhances verification accuracy and can be applied widely.     

Improving the representation of soil moisture by using a semi‐analytical infiltration model

Soil moisture is widely recognized as a fundamental variable governing the mass and energy fluxes between the land surface and the atmosphere. In this study, the soil moisture modelling at sub‐daily timescale is addressed by using an accurate representation of the infiltration component. For that, the semi‐analytical infiltration model proposed by Corradini et al. (1997) has been incorporated into a soil water balance model to simulate the evolution in time of surface and profile soil moisture. The performances of this new soil moisture model [soil water balance module‐semi‐analytical (SWBM‐SA)] are compared with those of a precedent version [SWBM‐Green–Ampt (GA)] where the GA approach was employed. Their capability to reproduce in situ soil moisture observations at three sites in Italy, Spain and France is analysed. Hourly observations of quality‐checked rainfall, temperature and soil moisture data for a 2‐year period are used for testing the modelling approaches. Specifically, different configurations for the calibration and validation of the models are adopted by varying a single parameter, that is, the saturated hydraulic conductivity. Results indicate that both SWBMs are able to reproduce satisfactorily the hourly soil moisture temporal pattern for the three sites with root mean square errors lower than 0.024 m³/m³ both in the calibration and validation periods. For all sites, the SWBM‐SA model outperforms the SWBM‐GA with an average reduction of the root mean square error of ~20%. Specifically, the higher improvement is observed for the French site for which in situ observations are measured at 30 cm depth, and this is attributed to the capability of the SA infiltration model to simulate the time evolution of the whole soil moisture profile. The reasonable models performance coupled with the need to calibrate only a single parameter makes them useful tools for soil moisture simulation in different regions worldwide, also in scarcely gauged areas. Copyright © 2013 John Wiley & Sons, Ltd.