Regional assimilation of in situ observed soil moisture into the VIC model considering spatial variability

ABSTRACT

In order to improve the soil moisture (SM) modelling capacity, a regional SM assimilation scheme based on an empirical approach considering spatial variability was constructed to assimilate in situ observed SM data into a hydrological model. The daily variable infiltration capacity (VIC) model was built to simulate SM in the Upper Huai River Basin, China, with a resolution of 5 km × 5 km. Through synthetic assimilation experiments and validations, the assimilated SM was evaluated, and the assimilation feedback on evapotranspiration (ET) and streamflow are analysed and discussed. The results show that the assimilation scheme improved the SM modelling capacity, both spatially and temporally. Moreover, the simulated ET was continually affected by changes in SM simulation, and the streamflow predictions were improved after applying the SM assimilation scheme. This study demonstrates the potential value of in situ observations in SM assimilation, and provides valuable ways for improving hydrological simulations.     

Regional scale spatial and temporal variability of soil moisture in a prairie region

Understanding the dynamics of spatial and temporal variability of soil moisture at the regional scale and daily interval, respectively, has important implications for remote sensing calibration and validation missions as well as environmental modelling applications. The spatial and temporal variability of soil moisture was investigated in an agriculturally dominated region using an in‐situ soil moisture network located in central Saskatchewan, Canada. The study site evaluated three depths (5, 20, 50 cm) through 139 days producing a high spatial and temporal resolution data set, which were analysed using statistical and geostatistical means. Processes affecting standard deviation at the 5‐cm depth were different from the 20‐cm and 50‐cm depths. Deeper soil measurements were well correlated through the field season. Further analysis demonstrated that lag time to maximum correlation between soil depths increased through the field season. Temporal autocorrelation was approximately twice as long at depth compared to surface soil moisture as measured by the e‐folding frequency. Spatial correlation was highest under wet conditions caused by uniform rainfall events with low coefficient of variation. Overall soil moisture spatial and temporal variability was explained well by rainfall events and antecedent soil moisture conditions throughout the Kenaston soil moisture network. It is expected that the results of this study will support future remote sensing calibration and validation missions, data assimilation, as well as hydrologic model parameterization for use in agricultural regions. Copyright © 2016 John Wiley & Sons, Ltd.     

The role of landscape morphology on soil moisture variability in semi-arid ecosystems

Previous studies on semi-arid ecosystems have shown high values of soil moisture variability (SMV) primarily induced by the combined effects of non-uniform precipitation, incoming solar radiation, and soil and vegetation properties. However, the relative impact of these various factors on SMV has been difficult to evaluate due to limited availability of field data. In addition, only a limited number of studies have analysed the role of landscape morphology on SMV. Here we use numerical simulations of a simple hydrological model, the Bucket Grassland Model, to systematically analyse the effect of each contributing factor on SMV on two different landscape morphologies. The two different landform morphologies represent landscapes dominated respectively by either diffusive erosion or fluvial erosion processes. We conducted various simulations driven by a stochastically generated 100-year climate time series, which is long enough to capture climatic fluctuations, in order to understand the effect of various soil moisture controlling factors on the spatiotemporal SMV. Our modelling results show that the fluvial dominated landscapes promote higher spatial SMV than the diffusive dominated ones. Further, the role of landform morphology on SMV is more pronounced in regions where the spatial variability of incoming solar radiation and precipitation is high.     

Factors affecting soil moisture spatial variability for a humid forest hillslope

Soil moisture is an important variable in explaining hydrological processes at hillslope scale. The distribution of soil moisture along a hillslope is related to the spatial distribution of the soil properties, the topography, the soil depth, and the vegetation. In order to investigate the factors affecting soil moisture, various environmental data were collected from a humid forest hillslope in this study. Several factors (the wetness index; the contributing area; the local slope; the soil depth; the composition of sand, silt, and clay; the scaling parameter; the hydraulic conductivity; the tree diameter at breast height; and the total weighted basal area) were evaluated for their effect on soil moisture and its distribution over the hillslope at depths of 10, 30, and 60 cm. Both linear correlation analysis and empirical orthogonal function analysis indicated that the soil texture was a dominant factor in soil moisture distribution. The impact of soil hydraulic conductivity was important for all soil moisture ranges at a depth of 30 cm, but those at 10 and 60 cm were limited to very wet and dry conditions, respectively. The relationships of the various factors with the spatial variability of soil moisture indicated the existence of a threshold soil moisture that is related to the composition of the soil and the factors related to the distribution of water in the study area.     

Role of time variant and invariant terms on soil moisture spatio-temporal variability in a subhumid coastal wetland

The knowledge of soil moisture spatio-temporal variability is highly relevant for water resources management. This paper reports an analysis of the spatial–temporal variability of soil moisture data for a small to medium-scale soil-sensors network in a coastal wetland of southwestern Spain. Measurements were taken from five sites located in the Doñana National Park over the time-period of one hydrological year from September 2017 to September 2018. The total area of the soil-sensors network shows an extension about 25 × 3 km. Soil moisture data was separated into time invariant (the temporal mean of the whole period at each site) and time-variant terms (the deviations of soil moisture from the mean, or anomalies). The time-invariant component was generally the main contributor to the total spatial variance of soil moisture and it was mostly controlled by the groundwater levels in the area. Nevertheless, the time variant terms have a huge effect on soil moisture variability in very dry states. Characteristic convex time-dependent patterns for this field site were found between spatially averaged soil moisture and its variability. This information could be used for the up and downscaling of soil moisture from satellite data. Those patterns of relation between spatial mean and variability of soil moisture were only affected by heavy rainfalls giving rise to hysteretic behaviour. This study shows that even though groundwater level is a time-variant variable, it significantly affects soil moisture's time-variant but also time-invariant terms due to the different average groundwater level depths at the different sites.     

Spatial patterns and temporal stability of soil moisture across a range of scales in a semi‐arid environment

The antecedent soil moisture status of a catchment is an important factor in hydrological modelling. Traditional Hortonian infiltration models assume that the initial moisture content is constant across the whole catchment, despite the fact that even in small catchments antecedent soil moisture exhibits tremendous spatial heterogeneity. Spatial patterns of soil water distribution across three transects (two in a burnt area and one in an unburnt area) in a semi‐arid area were studied. At the transect scale, when the factors affecting soil moisture were limited to topographical position or local topography, spatial patterns showed time stability, but when other factors, such as vegetation, were taken into account, the spatial patterns became time unstable. At the point scale, and in the same areas, topographical position was the main factor controlling time stability. Scale dependence of time stability was studied and local topography and vegetation presence were observed to play an important role for the correlation between consecutive measures depending on the scale. Copyright © 2000 John Wiley & Sons, Ltd.     

Spatial variability in soil moisture as predicted from airborne thematic mapper (ATM) data

In the first part of this project, the extent to which moisture content of alluvial soils could be predicted from imagery derived from an airborne thematic mapper (ATM) was investigated. From sampling done on the same day as the flight, it was found that digital numbers derived from the thermal channel (waveband 11) were strongly correlated with gravimetric moisture content. From sampling three fields of contrasting land cover, the relationship between waveband 11 values and moisture content was found to be independent of land cover type. Spatial variation in waveband 11 values and thus moisture content were related to palaeochannel patterns on the alluvial land. This was investigated by deriving variograms for long transects from each of the three investigated fields. The range and sills of the variograms are shown to express the nature and pattern of palaeochannels. By the application of such geostatistical techniques, high resolution imagery can thus be used to quantify palaeochannel characteristics on alluvial land.     

Spatial soil moisture scaling structure during Soil Moisture Experiment 2005

Soil moisture state and variability control many hydrological and ecological processes as well as exchanges of energy and water between the land surface and the atmosphere. However, its state and variability are poorly understood at spatial scales larger than the fields (i.e. 1 km²) as well as the ability to extrapolate field scale to larger spatial scales. This study investigates soil moisture profiles, their spatial organization, and physical drivers of variability within the Walnut Creek watershed, Iowa, during Soil Moisture Experiment 2005 and relates the watershed scale findings to previous field‐scale results. For all depths, the watershed soil moisture variability was negatively correlated with the watershed mean soil moisture and followed an exponential relationship that was nearly identical to that for field scales. This relationship differed during drying and wetting. While the overall time stability characteristics were improved with observation depth, the relatively wet and dry locations were consistent for all depths. The most time stable locations, capturing the mean soil moisture of the watershed within ± 0·9% volumetric soil moisture, were typically found on hill slopes regardless of vegetation type. These mild slope locations consistently preserve the time stability patterns from field to watershed scales. Soil properties also appear to impact stability but the findings are sensitive to local variations that may not be well defined by existing soil maps. Copyright © 2010 John Wiley & Sons, Ltd.     

新疆艾比湖湿地土壤水盐空间变异性分析

为揭示艾比湖湿地土壤退化程度空间分布特征,在离艾比湖湖滨5~15 km,绕湖一周160 km范围内,以湖心质点为中心,将艾比湖划分为东北、东南、西南、西北4个区域,采用传统统计学和地统计学相结合的方法对表层(0~20 cm)土壤盐分、含水量与p H的空间分异特征进行研究.结果表明:绕湖一周不同区域的土壤盐分均属中等变异强度;土壤含水量在西北部属强变异性,而东北、东南和西南部均属中等变异强度;土壤p H在不同区域内均属弱变异强度.绕湖一周除西北部土壤盐分的半方差理论模型较符合球状模型外,其它区域土壤盐分、含水量和p H均符合高斯模型;受结构性因素影响,不同区域土壤盐分、含水量和p H均具有较强的空间相关性;西南部土壤盐分、含水量和p H的Moran's I系数比其它区域的波动大,表明空间相关性较强.艾比湖湿地常年大风、干旱、缺水及沙化盐化的自然因素与引水围堰、种植耐盐碱植物的人为活动造成了采样区表层土壤盐分、含水量和p H的空间分布多呈现不规则条带状格局.艾比湖湿地土壤以盐土为主,重度盐化土次之,土壤盐渍化日益严重.     

Temporal and spatial features of the soil moisture in boreal spring in eastern China

Soil moisture data of 45 years from European Center for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40) and the in situ observational data are used to study the temporal and spatial characteristics of the soil moisture in boreal spring in the area to the east of 100°E in China. Results show that ERA-40 soil moisture well reproduces the temporal and spatial features of observations. ERA-40 data capture the spatial pattern that the soils in Northeast China and Southwest China are wetter than those...     

Effects of spatial aggregation of soil spatial information on watershed hydrological modelling

Many researchers have examined the impact of detailed soil spatial information on hydrological modelling due to the fact that such information serves as important input to hydrological modelling, yet is difficult and expensive to obtain. Most research has focused on the effects at single scales; however, the effects in the context of spatial aggregation across different scales are largely missing. This paper examines such effects by comparing the simulated runoffs across scales from watershed models based on two different levels of soil spatial information: the 10‐m‐resolution soil data derived from the Soil‐Land Inference Model (SoLIM) and the 1:24000 scale Soil Survey Geographic (SSURGO) database in the United States. The study was conducted at three different spatial scales: two at different watershed size levels (referred to as full watershed and sub‐basin, respectively) and one at the model minimum simulation unit level. A fully distributed hydrologic model (WetSpa) and a semi‐distributed model (SWAT) were used to assess the effects. The results show that at the minimum simulation unit level the differences in simulated runoff are large, but the differences gradually decrease as the spatial scale of the simulation units increases. For sub‐basins larger than 10 km² in the study area, stream flows simulated by spatially detailed SoLIM soil data do not significantly vary from those by SSURGO. The effects of spatial scale are shown to correlate with aggregation effect of the watershed routing process. The unique findings of this paper provide an important and unified perspective on the different views reported in the literature concerning how spatial detail of soil data affects watershed modelling. Different views result from different scales at which those studies were conducted. In addition, the findings offer a potentially useful basis for selecting details of soil spatial information appropriate for watershed modelling at a given scale. Copyright © 2011 John Wiley & Sons, Ltd.     

A mathematical model of soil moisture spatial distribution on the hill slopes of the Loess Plateau

Based on important factors that affect soil moisture spatial distribution, such as the slope gradients, land use, vegetation cover, and surface water diffusion characteristics together with field measurements of soil moisture data obtained from the surface soil under different land use structures, a soil moisture spatial distribution model was established. The diffusion degree coefficient of surface water for different vegetations was estimated from soil moisture values obtained from field measurements. The model can be solved using the finite unit method. The soil moisture spatial distribution on the hill slopes in the Loess Plateau were simulated by the model. A comparison of the simulated values with measurement data shows that the model is a good fit.     

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

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

Soil moisture temporal stability and spatio-temporal variability about a typical subalpine ecosystem in northwestern China

Knowledge of the spatial–temporal variability of soil water content is critical for water management and restoration of vegetation in semi-arid areas. Using the temporal stability method, we investigated soil water relations and spatial–temporal variability of volumetric soil water content (VSWC) in the grassland–shrubland–forest transect at a typical semi-arid subalpine ecosystem in the Qilian Mountains, northwestern China. The VSWC was measured on 48 occasions to a depth of 70 cm at 50 locations along a 240-m transect during the 2016–2017 growing seasons. Results revealed that temporal variability in VSWC in the same soil layer in the three vegetation types and averaged across vegetation types tended to exhibit similar patterns of a decrease with increasing soil depth. Temporal stability in each vegetation type was stronger with an increase in soil depth. However, the results of temporal stability determined with standard deviation of relative difference (SDRD) disagreed with those based on the Spearman's rank correlation coefficient; the forest site had the highest Spearman rank correlation coefficient while the shrubland—the smallest SDRD in the 0–20 cm soil layer. Correlation analyses of VSWCs between two vegetation types indicated that soil water was related among all three vegetation types at the 0–20, and 0–70 cm soil layer, but in the 20–40 and 40–70 cm soil layers, significant correlation (p < .01) occurred only between adjacent vegetation types. In the upper soil layer (0–20 cm), soil water relations were mainly affected by surface runoff. In the lower soil layer (20–40 and 40–70 cm), soil water relations among the three vegetation types were highly complex, and probably resulting from a combination of root distribution and activity, interflow, and the impact of deep soil freeze–thaw dynamics. These results suggest that the factors affecting soil water are complex, and further research should address the relative importance of and interactions among different determining factors.     

Time stability and variability of Electronically Scanned Thinned Array Radiometer soil moisture during Southern Great Plains hydrology experiments

Variability and time‐stability analysis for field‐scale (800 m) Electronically Scanned Thinned Array Radiometer soil moisture within a satellite scale footprint (∼ 50 km) were quantified using observations from the Southern Great Plains Hydrology Experiment 1997 and 1999 (SGP97 and SGP99). The pixels' time‐stability properties were examined with respect to soil, vegetation and topographic parameters in order to determine which physical parameters can be used to identify good candidate observation locations for validating soil moisture from satellite observations and global‐scale model output. The results show that the time‐stability concept remains valid at the satellite scale. The root mean square error values were 1·47, 1·51, 1·93 and 2·32% for the 1st, 2nd, 50th and 100th most stable fields, respectively. The most stable locations had sand and clay percentages consistent with sandy loam soils and moderate to high normalized difference vegetation index values. Neither land cover nor topography properties could be used to identify potentially stable fields in the study region. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Spatial scale effect on the upper soil effective parameters of a distributed hydrological model

Nonlinear dynamics and spatial variability in hydrological systems make the formulation of scaling theories difficult. Therefore, the development of knowledge related to scale effects, scaling techniques, parameterization and linkages of parameters across scales is highly relevant. The main purpose of this work is to analyse the spatial effect of the static storage capacity parameter H_u and the saturated hydraulic conductivity parameter k_s from microscale (sub‐grid level) to mesoscale (grid level) and its implication to the definition of an optimum cell size. These two parameters describe the upper soil water characteristics in the infiltration process conceptualization of the TETIS hydrological model. At microscale, the spatial heterogeneity of H_u and k_s was obtained generating random parameter fields through probability distribution functions and a spatial dependence model with pre‐established correlation lengths. The effective parameters at mesoscale were calculated by solving the inverse problem for each parameter field. Results indicate that the adopted inverse formulation allows transferring the nonlinearity of the system from microscale to the mesoscale via non‐stationary effective parameters. Their values at each cell and time step are in the range of zero to the mean value of the parameter at microscale. The stochastic simulations showed that the variance of the estimated effective parameters decreases when the ratio between mesoscale cell size and correlation length at microscale increases. For a ratio greater than 1, we found cell sizes having the characteristics of a representative elementary area (REA); in such case, the microscale variability pattern did not affect the system response at mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.     

Configuration of hydrological processes in a steep hillslope through causality analysis of soil moisture measurements

The vertical and lateral profiles of temporal variations in soil moisture are important for understanding the hydrological process along hillside transects. In this study, relationships among measured soil moistures were explored to configure the hydrological contributions of different flowpaths. All the measured soil moistures included a common stochastic structure because rainfall, the hydrometeological driver, is the main factor that determines the soil moisture response feature, and the infiltration process through the topsoil at a shallow depth is also common in all measured soil moisture histories. Therefore, the relationships between the measured series are also affected by both rainfall and topsoil infiltration. The common stochastic structure of the soil moisture series was removed via a prewhitening procedure. A systematic analysis procedure is presented to delineate the exclusive causal relationships among multiple soil moisture measurements. A monitoring system based on multiplexed time domain reflectometry was used to obtain soil moisture time series along two transects on a steep hillslope during the rainy season. The application of the proposed method for monitoring points in two adjacent locations provided 8, 12, 14, and 13, 16, 22 causal relationships for vertical, lateral in parallel, and diagonal directions, respectively, along the two transects. The point‐based contributions of the internal flowpath can be evaluated as the correlation is normalized in the context of inflow and outflow. The hydrological processes in the soil layer, vertical flow, lateral flow, downslope recharge, and return flow were quantified, and the relative importance of each hydrological component was determined to improve our understanding of the hydrological processes along the two transects of the study area. Copyright © 2011 John Wiley & Sons, Ltd.