A land surface soil moisture data assimilation framework in consideration of the model subgrid-scale heterogeneity and soil water thawing and freezing

The Ensemble Kalman Filter (EnKF) is well known and widely used in land data assimilation for its high precision and simple operation. The land surface models used as the forecast operator in a land data assimilation system are usually designed to consider the model subgrid-heterogeneity and soil water thawing and freezing. To neglect their effects could lead to some errors in soil moisture assimilation. The dual EnKF method is employed in soil moisture data assimilation to build a soil moisture data as- similation framework based on the NCAR Community Land Model version 2.0 (CLM 2.0) in considera- tion of the effects of the model subgrid-heterogeneity and soil water thawing and freezing: Liquid volumetric soil moisture content in a given fraction is assimilated through the state filter process, while solid volumetric soil moisture content in the same fraction and solid/liquid volumetric soil moisture in the other fractions are optimized by the parameter filter. Preliminary experiments show that this dual EnKF-based assimilation framework can assimilate soil moisture more effectively and precisely than the usual EnKF-based assimilation framework without considering the model subgrid-scale heteroge- neity and soil water thawing and freezing. With the improvement of soil moisture simulation, the soil temperature-simulated precision can be also improved to some extent.     

Development of the universal and simplified soil model coupling heat and water transport

It is very important to develop a universal soil model with higher simplicity and more accuracy, which can be widely applied to very general cases such as wet or dry soil, frozen or unfrozen soil and homogeneous or heterogeneous soil. Firstly in this study, based on analysis of both magnitude order and the numerical simulation results, the universal and simplified soil model (USSM) coupling heat and mass transport processes is developed. Secondly, in order to avoid the greater uncertainty caused by the phase change term in numerical iteration process for the model solution obtaining, new version of the universal simplified soil model (NUSSM) is further derived through variables transformation, and accordingly a more efficient numerical scheme for the new version is designed well. The simulation results from the NUSSM agree with the results from more complicated and accurate soil model very well, also reasonably reproduce the observed data under widely real conditions. The new version model, because of its simplicity, will match for the development of land surface model.     

A land surface model incorporated with soil freeze/thaw and its application in GAME/Tibet

Land surface process is of great importance in global climate change, moisture and heat exchange in the interface of the earth and atmosphere, human impacts on the environment and eco- system, etc. Soil freeze/thaw plays an important role in cold land surface processes. In this work the diurnal freeze/thaw effects on energy partition in the context of GAME/Tibet are studied. A sophisti- cated land surface model is developed, the particular aspect of which is its physical consideration of soil freeze/thaw and vapor flux. The simultaneous water and heat transfer soil sub-model not only reflects the water flow from unfrozen zone to frozen fringe in freezing/thawing soil, but also demon- strates the change of moisture and temperature field induced by vapor flux from high temperature zone to low temperature zone, which makes the model applicable for various circumstances. The modified Picard numerical method is employed to help with the water balance and convergence of the numerical scheme. Finally, the model is applied to analyze the diurnal energy and water cycle char- acteristics over the Tibetan Plateau using the Game/Tibet datasets observed in May and July of 1998. Heat and energy transfer simulation shows that: (i) There exists a negative feedback mechanism between soil freeze/thaw and soil temperature/ground heat flux; (ii) during freezing period all three heat fluxes do not vary apparently, in spite of the fact that the negative soil temperature is higher than that not considering soil freeze; (iii) during thawing period, ground heat flux increases, and sensible heat flux decreases, but latent heat flux does not change much; and (iv) during freezing period, soil temperature decreases, though ground heat flux increases.     

Canopy cover effects on local soil water dynamics in a tropical agroforestry system: Evaporation drives soil water isotopic enrichment

Despite the widely held assumption that trees negatively affect the local water budget in densely planted tree plantations, we still lack a clear understanding of the underlying processes by which canopy cover influences local soil water dynamics in more open, humid tropical ecosystems. In this study, we propose a new conceptual model that uses a combination of stable isotope and soil moisture measurements throughout the soil profile to assess potential mechanisms by which evaporation (of surface soil water and of canopy‐intercepted rainfall) affects the relationship between surface soil water isotopic enrichment (lc‐excess) and soil water content. Our conceptual model was derived from soil water data collected under deciduous and evergreen plants in a shade grown coffee agroforestry system in Costa Rica. Reduced soil moisture under shade trees during the “drier” season, coinciding when these trees were defoliated, was largely the result of increase soil water evaporation as indicated by the positive relationship between soil water content and lc‐excess of surface soil water. In contrast, the evergreen coffee shrubs had a higher leaf area index during the “drier” season, leading to enhanced rainfall interception and a negative relationship between lc‐excess and soil water content. During the wet season, there was no clear relationship between soil water content and between lc‐excess of surface soil water. Greater surface soil water under coffee during the dry season may, in part, explain greater preferential flow under coffee compared with under trees in conditions of low rainfall intensities. However, with increasing rainfall intensities during the wet season, there was no obvious difference in preferential flow between the two canopy covers. Results from this study indicate that our new conceptual model can be used to help disentangling the relative influence of canopy cover on local soil water isotopic composition and dynamics, yet also stresses the need for additional measurements to better resolve the underlying processes by which canopy structure influences local water dynamics.     

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.     

Evapotranspiration and runoff from large land areas: Land surface hydrology for atmospheric general circulation models

A land surface hydrology parameterization for use in atmospheric GCMs is presented. The parameterization incorporates subgrid scale variability in topography, soils, soil moisture and precipitation. The framework of the model is the statistical distribution of a topography-soils index, which controls the local water balance fluxes, and is therefore taken to represent the large land area. Spatially variable water balance fluxes are integrated with respect to the topography-soils index to yield our large scale parameterizations: water balance calculations are performed for a number of intervals of the topography-soils distribution, and interval responses are weighted by the probability of occurrence of the interval. Grid square averaged land surface fluxes result. The model functions independently as a macroscale water balance model. Runoff ratio and evapotranspiration efficiency parameterizations are derived and are shown to depend on the spatial variability of the above mentioned properties and processes, as well at the dynamics of land surface-atmosphere interactions.     

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.     

Experimental study and mathematical modelling of soluble chemical transfer from unsaturated/saturated soil to surface runoff

A one‐dimensional, two‐layer solute transport model is developed to simulate chemical transport process in an initially unsaturated soil with ponding water on the soil surface before surface runoff starts. The developed mathematical model is tested against a laboratory experiment. The infiltration and diffusion processes are mathematically lumped together and described by incomplete mixing parameters. Based on mass conservation and water balance equations, the model is developed to describe solute transport in a two‐zone layer, a ponding runoff zone and a soil mixing zone. The two‐zone layer is treated as one system to avoid describing the complicated chemical transport processes near the soil surface in the mixing zone. The proposed model was analytically solved, and the solutions agreed well with the experimental data. The developed experimental method and mathematical model were used to study the effect of the soil initial moisture saturation on chemical concentration in surface runoff. The study results indicated that, when the soil was initially saturated, chemical concentration in surface runoff was significantly (two orders of magnitude) higher than that with initially unsaturated soil, while the initial chemical concentrations at the two cases were of the same magnitude. The soil mixing depth for the initially unsaturated soil was much larger than that for the initially saturated soil, and the incomplete runoff mixing parameter was larger for the initially unsaturated soil. The higher the infiltration rate of the soil, the greater the infiltration‐related incomplete mixing parameter. According to the quantitative analysis, the soil mixing depth was found to be sensitive for both initially unsaturated and saturated soils, and the incomplete runoff mixing parameter was sensitive for initially saturated soil but not for the initially unsaturated soil; the incomplete infiltration mixing parameter behaved just the opposite. Some suggestions are made for reducing chemical loss from runoff. Copyright © 2010 John Wiley & Sons, Ltd.     

Hillslope soil water flowpaths and the dynamics of roadside soil cation pools influenced by road deicers

Over the past 60 years, road deicers (i.e. road salt) have been applied to roadways in high latitudes to improve road conditions in winter weather. However, the dissolution of road deicers in highway runoff creates waters with high concentrations of sodium, which can mobilize soil metals via soil cation‐exchange reactions. While several studies have detailed the interactions of road salt‐rich solutions and surface and ground waters, less attention has been given to how local hydrologic flowpaths can impact the delivery of these solutions to near‐road soils. Between 2013 and 2014, soil water samples were collected from a roadside transect of lysimeter nests in Pittsburgh, Pennsylvania (USA). Soil water samples were analysed for metal concentrations and resulting data used to examine cation dynamics. While patterns in soil water calcium and magnesium concentrations follow patterns in soil water sodium concentrations, additional processes influence patterns in soil water potassium concentrations. Specifically, we observe the highest calcium and magnesium concentrations in the deepest lysimeters, suggesting divalent cations are mobilized to, and potentially accumulate in, deeper soil horizons. In contrast, soil water potassium concentrations do not follow this pattern. Additionally, in all examined elements (Ca, Mg, K, Na, and Cl), the timing of concentration peaks appears be influenced by a combination of both distance from the roadside and sampling depth. These relationships not only suggest that multiple soil water flowpaths interact with our study transect but also confirm that road salt plumes persist and migrate following the road salting season. Characterizing the interactions of sodium‐rich solutions and roadside soil cation pools clarifies our understanding of metal dynamics in the roadside environment. A deeper understanding of these processes is necessary to effectively restore and manage watersheds as high total dissolved solid solutions (e.g. road deicing melt, unconventional natural gas brines, and marginal irrigation water) continue to influence hydrological systems. Copyright © 2016 John Wiley & Sons, Ltd.     

Modeling soil erosion using a spatially distributed model in a karst catchment of northwest Guangxi,China

Reliable assessment of the spatial distribution of soil erosion is important for making land management decisions, but it has not been thoroughly evaluated in karst geo‐environments. The objective of this study was to modify a physically based, spatially distributed erosion model, the revised Morgan, Morgan and Finney (RMMF) model, to estimate the superficial (as opposed to subsurface creep) soil erosion rates and their spatial patterns in a 1022 ha karst catchment in northwest Guangxi, China. Model parameters were calculated using local data in a raster geographic information system (GIS) framework. The cumulative runoff on each grid cell, as an input to the RMMF model for erosion computations, was computed using a combined flow algorithm that allowed for flow into multiple cells with a transfer grid considering infiltration and runoff seepage to the subsurface. The predicted spatial distributions of soil erosion rates were analyzed relative to land uses and slope zones. Results showed that the simulated effective runoff and annual soil erosion rates of hillslopes agreed well with the field observations and previous quantified redistribution rates with caesium‐137 (¹³⁷Cs). The estimated average effective runoff and annual erosion rate on hillslopes of the study catchment were 18 mm and 0.27 Mg ha^?1 yr^?1 during 2006–2007. Human disturbances played an important role in accelerating soil erosion rates with the average values ranged from 0.1 to 3.02 Mg ha^?1 yr^?1 for different land uses. The study indicated that the modified model was effective to predict superficial soil erosion rates in karst regions and the spatial distribution results could provide useful information for developing local soil and water conservation plans. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Method development for estimating soil organic carbon content in an alpine region using soil moisture data

The high soil organic carbon(SOC) content in alpine meadow can significantly change soil hydrothermal properties and further affect the soil temperature and moisture as well as the surface water and energy budget. Therefore, this study first introduces a parameterization scheme to describe the effect of SOC content on soil hydraulic and thermal parameters in a land surface model(LSM), and then the SOC content is estimated by minimizing the difference between observed and simulated surface-layer soil moisture. The accuracy of the estimated SOC content was evaluated using in situ observation data at a soil moisture and temperature-measuring network in Naqu, central Tibetan Plateau. Sensitivity experiments show that the optimum time window for stabilizing the estimation results cannot be shorter than three years. In the experimental area, the estimated SOC content can generally reflect the spatial distribution of the measurements, with a root mean square error of 0.099 m^3 m^-3, a mean bias of 0.043 m^3 m^-3, and a correlation coefficient of 0.695. The estimated SOC content is not sensitive to the temporal frequency of the soil moisture data input. Even if the temporal frequency is as low as that of current soil moisture products derived from passive microwave satellites, the estimation result is still stable. Therefore, by combining a high-quality satellite soil moisture product and a parameter optimization method, it is possible to obtain grid-scale effective parameter values, such as SOC content,for an LSM and improve the simulation ability of the LSM.     

Multi‐scale assimilation of root zone soil water predictions

When hydrology model parameters are determined, a traditional data assimilation method (such as Kalman filter) and a hydrology model can estimate the root zone soil water with uncertain state variables (such as initial soil water content). The simulated result can be quite good. However, when a key soil hydraulic property, such as the saturated hydraulic conductivity, is overestimated or underestimated, the traditional soil water assimilation process will produce a persistent bias in its predictions. In this paper, we present and demonstrate a new multi‐scale assimilation method by combining the direct insertion assimilation method, particle swarm optimisation (PSO) algorithm and Richards equation. We study the possibility of estimating root zone soil water with a multi‐scale assimilation method by using observed in situ data from the Wudaogou experiment station, Huaihe River Basin, China. The results indicate there is a persistent bias between simulated and observed values when the direct insertion assimilation surface soil water content is used to estimate root zone soil water contents. Using a multi‐scale assimilation method (PSO algorithm and direct insertion assimilation) and an assumed bottom boundary condition, the results show some obvious improvement, but the root mean square error is still relatively large. When the bottom boundary condition is similar to the actual situation, the multi‐scale assimilation method can well represent the root zone soil water content. The results indicate that the method is useful in estimating root zone soil water when available soil water data are limited to the surface layer and the initial soil water content even when the soil hydraulic conductivities are uncertain. Copyright © 2011 John Wiley & Sons, Ltd.     

Insight into sediment transport processes on saline rangeland hillslopes using three‐dimensional soil microtopography changes

In arid and semi‐arid rangeland environments, an accurate understanding of runoff generation and sediment transport processes is key to developing effective management actions and addressing ecosystem response to changes. Yet, many primary processes (namely sheet and splash and concentrated flow erosion, as well as deposition) are still poorly understood due to a historic lack of measurement techniques capable of parsing total soil loss into these primary processes. Current knowledge gaps can be addressed by combining traditional erosion and runoff measurement techniques with image‐based three‐dimensional (3D) soil surface reconstructions. In this study, data (hydrology, erosion and high‐resolution surface microtopography changes) from rainfall simulation experiments on 24 plots in saline rangelands communities of the Upper Colorado River Basin were used to improve understanding on various sediment transport processes. A series of surface change metrics were developed to quantify and characterize various erosion and transport processes (e.g. plot‐wide versus concentrated flow detachment and deposition) and were related to hydrology and biotic and abiotic land surface characteristics. In general, erosivity controlled detachment and transport processes while factors modulating surface roughness such as vegetation controlled deposition. The extent of the channel network was a positive function of slope, discharge and vegetation. Vegetation may deflect runoff in many flow paths but promoted deposition. From a management perspective, this study suggests that effective runoff soil and salt load reduction strategies should aim to promote deposition of transported sediments rather than reducing detachment which might not be feasible in these resource‐limited environments. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamics of soil erosion rates and controlling factors in the Northern Ethiopian Highlands – towards a sediment budget

This paper analyses the factors that control rates and extent of soil erosion processes in the 199 ha May Zegzeg catchment near Hagere Selam in the Tigray Highlands (Northern Ethiopia). This catchment, characterized by high elevations (2100–2650 m a.s.l.) and a subhorizontal structural relief, is typical for the Northern Ethiopian Highlands. Soil loss rates due to various erosion processes, as well as sediment yield rates and rates of sediment deposition within the catchment (essentially induced by recent soil conservation activities), were measured using a range of geomorphological methods. The area‐weighted average rate of soil erosion by water in the catchment, measured over four years (1998–2001), is 14·8 t ha^?1 y^?1, which accounts for 98% of the change in potential energy of the landscape. Considering these soil loss rates by water, 28% is due to gully erosion. Other geomorphic processes, such as tillage erosion and rock fragment displacement by gravity and livestock trampling, are also important, either within certain land units, or for their impact on agricultural productivity. Estimated mean sediment deposition rate within the catchment equals 9·2 t ha^?1 y^?1. Calculated sediment yield (5·6 t ha^?1 y^?1) is similar to sediment yield measured in nearby catchments. Seventy‐four percent of total soil loss by sheet and rill erosion is trapped in exclosures and behind stone bunds. The anthropogenic factor is dominant in controlling present‐day erosion processes in the Northern Ethiopian Highlands. Human activities have led to an overall increase in erosion process intensities, but, through targeted interventions, rural society is now well on the way to control and reverse the degradation processes, as can be demonstrated through the sediment budget. Copyright © 2007 John Wiley & Sons, Ltd.     

Ability of soil bacterial composition as an indicator of levels of soil erosion in a badland

Calanchi(plural of calanco) are typical Italian badlands created by a combination of morphogenetic processes(rill and interrill erosion, gullying, piping, and mass movements) mainly originated by the effect of water. Calanchi are characterized by the sparse and patchy distribution of vegetation, and, in interplant areas, the soil surface is colonized by an association of organisms known as biological soil crust(BSC). A morphometric analysis of 45 basins in the studied calanchi area, based on a h...     

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.     

An improved process-based representation of stream solute transport in the soil and water assessment tools

Hydrological models have long been used to study the interactions between land, surface and groundwater systems, and to predict and manage water quantity and quality. The soil and water assessment tool (SWAT), a widely used hydrological model, can simulate various ecohydrological processes on land and subsequently route the water quality constituents through surface and subsurface waters. So far, in-stream solute transport algorithms of the SWAT model have only been minimally revised, even though it has been acknowledged that an improvement of in-stream process representation can contribute to better model performance with respect to water quality. In this study, we aim to incorporate a new and improved solute transport model into the SWAT model framework. The new process-based model was developed using in-stream process equations from two well established models—the One-dimensional Transport with Inflow and Storage model and the Enhanced Stream Water Quality Model. The modified SWAT model (Mir-SWAT) was tested for water quality predictions in a study watershed in Germany. Compared to the standard SWAT model, Mir-SWAT improved dissolved oxygen (DO) predictions by removing extreme low values of DO (<6 mg/L) simulated by SWAT. Phosphate concentration peaks were reduced during high flows and a better match of daily predicted and measured values was attained using the Mir-SWAT model (R² = 0.17, NSE = −0.65, RSR = 1.29 with SWAT; R² = 0.28, NSE = −0.04, RSR = 1.02 with Mir-SWAT). In addition, Mir-SWAT performed better than the SWAT model in terms of Chlorophyll-a content particularly during winter months, improving the NSE and RSR for monthly average Chl-a by 74 and 42%, respectively. With the new model improvements, we aim to increase confidence in the stream solute transport component of the model, improve the understanding of nutrient dynamics in the stream, and to extend the applicability of SWAT for reach-scale analysis and management.     

Simulated soil erosion from a semiarid typical steppe watershed using an integrated aeolian and fluvial prediction model

Total soil erosion is a result of both aeolian and fluvial processes, which is particularly true in semiarid regions. However, although physically interrelated, these two processes have conventionally been studied and modelled independently. Recently, a few researchers highlighted the importance and need of considering both processes in concert as well as their interactions, but they did not give specific modelling approaches or algorithms. The objectives of this study were to (1) formulate an integrated aeolian and fluvial prediction (IAFP) model, (2) parameterize the IAFP model for a semiarid steppe watershed located in northeastern China by using literature and historical data and (3) use the model to predict soil erosion in the watershed and assess the sensitivity of predicted erosion to environmental factors such as soil moisture and vegetation coverage. The results indicated that the IAFP model can capture the dynamic interactions between aeolian and fluvial erosion processes. For the study watershed, the model predicted a higher occurrence frequency of fluvial erosion than that of aeolian erosion and showed that these two processes almost equivalently contributed to the average total erosion of 0.07 mm year^?1 across the simulation period. The ‘existing’ vegetation cover can provide an overall good protection of the soils, although the vegetation cover was predicted to play a larger role in a drier than a wetter year as well as in controlling aeolian than fluvial erosion. In addition, soil erosion was predicted to be more sensitive to soil moisture than land coverage. A soil moisture level of 0.23–0.25 was determined to be the probable switch point from aeolian‐to fluvial‐dominant process or vice versa. Copyright © 2012 John Wiley & Sons, Ltd.