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

2.
Our understanding of hydraulic properties of peat soils is limited compared with that of mineral substrates. In this study, we aimed to deduce possible alterations of hydraulic properties of peat soils following degradation resulting from peat drainage and aeration. A data set of peat hydraulic properties (188 soil water retention curves [SWRCs], 71 unsaturated hydraulic conductivity curves [UHCs], and 256 saturated hydraulic conductivity [Ks] values) was assembled from the literature; the obtained data originated from peat samples with an organic matter (OM) content ranging from 23 to 97 wt% (weight percent; and according variation in bulk density) representing various degrees of peat degradation. The Mualem‐van Genuchten model was employed to describe the SWRCs and UHCs. The results show that the hydraulic parameters of peat soils vary over a wide range confirming the pronounced diversity of peat. Peat decomposition significantly modifies all hydraulic parameters. A bulk density of approximately 0.2 g cm?3 was identified as a critical threshold point; above and below this value, macroporosity and hydraulic parameters follow different functions with bulk density. Pedotransfer functions based on physical peat properties (e.g., bulk density and soil depth) separately computed for bog and fen peat have significantly lower mean square errors than functions obtained from the complete data set, which indicates that not only the status of peat decomposition but also the peat‐forming plants have a large effect on hydraulic properties. The SWRCs of samples with a bulk density of less than 0.2 g cm?3 could be grouped into two to five classes for each peat type (botanical composition). The remaining SWRCs originating from samples with a bulk density of larger than 0.2 g cm?3 could be classified into one group. The Mualem‐van Genuchten parameter values of α can be used to estimate Ks if no Ks data are available. In conclusion, the derived pedotransfer functions provide a solid instrument to derive hydraulic parameter values from easily measurable quantities; however, additional research is required to reduce uncertainty.  相似文献   

3.
  相似文献   

4.
The Beerkan method based on in situ single‐ring water infiltration experiments along with the relevant specific Beerkan estimation of soil transfer parameters (BEST) algorithm is attractive for simple soil hydraulic characterization. However, the BEST algorithm may lead to erroneous or null values for the saturated hydraulic conductivity and sorptivity especially when there are only few infiltration data points under the transient flow state, either for sandy soil or soils in wet conditions. This study developed an alternative algorithm for analysis of the Beerkan infiltration experiment referred to as BEST‐generalized likelihood uncertainty estimation (GLUE). The proposed method estimates the scale parameters of van Genuchten water retention and Brooks–Corey hydraulic conductivity functions through the GLUE methodology. The GLUE method is a Bayesian Monte Carlo parameter estimation technique that makes use of a likelihood function to measure the goodness‐of‐fit between modelled and observed data. The results showed that using a combination of three different likelihood measurements based on observed transient flow, steady‐state flow and experimental steady‐state infiltration rate made the BEST‐GLUE procedure capable of performing an efficient inverse analysis of Beerkan infiltration experiments. Therefore, it is more applicable for a wider range of soils with contrasting texture, structure, and initial and saturated water content. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

5.
  相似文献   

6.
The calculation of the relative hydraulic conductivity function based on water retention data is an attractive and widely used approach, since direct measurements of unsaturated conductivities are difficult. We show theoretically under which conditions an air-entry value for water retention data is definitely required when using the statistical approach of Mualem. Moreover we rigorously specify the conditions for which the classical van Genuchten–Mualem model leads to wrong predictions of relative hydraulic conductivity and, hence, an alternative formulation including an air-entry value should be used. Significant consequences are demonstrated for the inverse parameter estimation based on multistep outflow experiments. Furthermore it is shown that the use of a physically correct formulation of the water retention curve including an air-entry value and the derived hydraulic conductivity function influences not only the stability of numerical simulations but also their final results. This is especially grave as simulations with van Genuchten–Mualem parameters are frequently used to compare experiments and simulations and to draw conclusions on the correctness of Richards’ equation.  相似文献   

7.
The analytic element method is well suited for the Gardner hydraulic conductivity function, but is limited in describing real soils. Therefore, parameter equivalence between the van Genuchten and Gardner hydraulic conductivity functions is explored for the case of steady vertical flow through a homogeneous medium with a single inclusion, i.e., a binary soil. The inclusion has different hydraulic parameters than the background medium. Equivalence is established using three methods: (1) effective capillary drive; (2) capillary length; (3) and a least-squares optimization method that aims to fit a Gardner function to a corresponding van Genuchten function by minimizing the difference in log conductivity over a specified pressure range. Comparisons between hydraulic models are made based on scatterplots of pressure head and the vertical Darcian flux obtained using a finite-element numerical solution with both constitutive relations. For applicability of an equivalent Gardner function over a broad range of pressure heads, the crossover pressure must be maintained between the two parametric functions. The crossover pressure is defined as the pressure in which the hydraulic conductivity of the inclusion is equal to the background. It can be shown that a hybrid methodology of preserving the crossover pressure exactly and using the effective capillary drive will result in hydraulic parameters that are easily obtained and provide good agreement between the conductivity functions of the GR model to the VG model.  相似文献   

8.
Consistent parameter constraints for soil hydraulic functions   总被引：2，自引：0，他引：2  
Parameters of functions to describe soil hydraulic properties are derived from measurements by means of parameter estimation. Of crucial importance here is the choice of correct constraints in the parameter space. Often, the parameters are mere shape parameters without physical meaning, giving flexibility to the model. A fundamental requirement is that the hydraulic functions are monotonic: the retention function and the conductivity function can only decrease as the capillary suction increases. A stricter physical requirement for the conductivity function is that its decrease with respect to saturation is at least linear. This linear decrease would occur if all pores of a capillary bundle had an equal radius. In the first part of this contribution, we derive constraints for the so-called tortuosity parameter of the Mualem conductivity model, which allow highest possible flexibility on one hand and guarantee physical consistency on the other hand. In combination with the retention functions of Brooks and Corey, van Genuchten, or Durner, such a constraint can be expressed as a function of the pore-size distribution parameters. In the second part, we show that a common modification of retention models, which is applied to reach zero water content at finite suction, can lead to the physically unrealistic case of increasing water content with increasing suction. We propose a solution for this problem by slightly modifying these models and introducing a correct parameter constraint.  相似文献   

9.
Evaporation from mosses and lichens can form a major component of the water balance, especially in ecosystems where mosses and lichens often grow abundantly, such as tundra, deserts and bogs. To facilitate moss representation in hydrological models, we parameterized the unsaturated hydraulic properties of mosses and lichens such that the capillary water flow through moss and lichen material during evaporation could be assessed. We derived the Mualem‐van Genuchten parameters of the drying retention and the hydraulic conductivity functions of four xerophilous moss species and one lichen species. The shape parameters of the retention functions (2.17 < n < 2.35 and 0.08 < α < 0.13 cm?1) ranged between values that are typical for sandy loam and loamy sand. The shapes of the hydraulic conductivity functions of moss and lichen species diverged from those of mineral soils, because of strong negative pore‐connectivity parameters (?2.840 < l < ?2.175) and low hydraulic conductivities at slightly negative pressure heads (0.016 < K0 < 0.280 cm/d). These K0 values are surprisingly low, considering that mosses are very porous. However, during evaporation, large pores and voids were air filled and did not participate in capillary water flow. Small K0 values cause mosses and lichens to be conservative with water during wet conditions, thus tempering evaporation compared to mineral soils. On the other hand, under dry conditions, mosses and lichens are able to maintain a moisture supply from the soil, leading to a higher evaporation rate than mineral soils. Hence, the modulating effect of mosses on evaporation possibly differs between wet and dry climates. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

10.
The capillary-sorption potential and the unsaturated hydraulic conductivity of soils as functions of soil water content are derived for forest, meadow, and agricultural ecosystems of the Ivankovo Reservoir watershed. The parameters of van Genuchten-Mualem equations were determined for the same soils based on data on soil particle size distribution and density. Computer code HYDRUS-1D and the obtained data were used to calculate the values of field capacity of soils in forest, meadow, and agricultural ecosystems in the drainage basin of the Ivankovo Reservoir.  相似文献   

11.
The selective radius shift model was used to relate changes in mineral volume due to precipitation/dissolution reactions to changes in hydraulic properties affecting flow in porous media. The model accounts for (i) precipitation/dissolution taking place only in the water-filled part of the pore space and further that (ii) the amount of mineral precipitation/dissolution within a pore depends on the local pore volume. The pore bundle concept was used to connect pore-scale changes to macroscopic soil hydraulic properties. Precipitation/dissolution induces changes in the pore radii of water-filled pores and, consequently, in the effective porosity. In a time step of the numerical model, mineral reactions lead to a discontinuous pore-size distribution because only the water-filled pores are affected. The pore-size distribution is converted back to a soil moisture characteristic function to which a new water retention curve is fitted under physically plausible constraints. The model equations were derived for the commonly used van Genuchten/Mualem hydraulic properties. Together with a mixed-form solution of Richards’ equation for aqueous phase flow, the model was implemented into the geochemical modelling framework PHREEQC, thereby making available PHREEQC’s comprehensive geochemical reactions. Example applications include kinetic halite dissolution and calcite precipitation as a consequence of cation exchange. These applications showed marked changes in the soil’s hydraulic properties due to mineral precipitation/dissolution and the dependency of these changes on water contents. The simulations also revealed the strong influence of the degree of saturation on the development of the saturated hydraulic conductivity through its quadratic dependency on the van Genuchten parameter α. Furthermore, it was shown that the unsaturated hydraulic conductivity at fixed reduced water content can even increase during precipitation due to changes in the pore-size distribution.  相似文献   

12.
Abstract. A simple closed-form expression relating saturated hydraulic conductivity to the van Genuchten capillary retention model parameters is derived. Application of this equation to an experimental data set shows reasonable agreement between measured and predicted saturated conductivity values. The proposed equation provides a consistent theoretical basis for estimating both saturated and unsaturated hydraulic conductivity from statistical pore structure models.  相似文献   

13.
The Mualem and the Burdine hydraulic conductivity prediction models are considered in combination with the van Genuchten analytical retention curve, as well as the Brooks and Corey prediction model. An equivalence is presented between the retention curves of these models. A comparative study follows between hydraulic conductivities that are based on equivalent retention curves. A unified presentation of prediction models provides a framework for the whole analysis. The treatment of the equivalence problem consists in a minimization procedure characterized by uncoupling of the parameters and analytical evaluation of the objective function. Exact analytical equivalence relations are given for significant parts of the parameter ranges, and, for the remaining parts, analytical approximations are proposed. The comparisons between hydraulic conductivities are carried out via an inequality analysis. It is shown that the hydraulic conductivity of the Burdine model is less than that of the other models for extended ranges of equivalent parameters.  相似文献   

14.
Soil hydraulic parameter values for large‐scale modelling cannot be obtained by direct methods. Pedotransfer functions (PTFs) that relate soil hydraulic properties (SHPs) to generally available soil texture data may provide an alternative. A considerable number of PTF models has been developed, the application of three recent PTFs is evaluated. As a first step sets of SHPs derived from the PTFs are compared with measured sets of SHPs for three sites. No good agreement was found statistically between measured and PTF results or between PTF results. As a second step and from a practical point of view results for three hydrologically functional variables were compared and evaluated. The three selected functional variables are saturated hydraulic conductivity, k0, in relation to infiltration excess runoff, available soil water amounts for evapotranspiration and water table depth for a specified upward flux or capillary rise. Derived k0 distributions from PTFs show substantially less variance than from the measured data at all three sites. This can have a considerable impact on infiltration excess runoff, depending on the actual rainfall regime. Simulated available soil water amounts for evapotranspiration for some combinations of PTFs and sites are close to those obtained for measured SHPs, however, no consistency in results can be detected. Water table depths for specified upward flux densities using PTF derived SHPs are generally deeper than those based on measured SHPs and means a potentially higher water availability. Overall, differences in capillary rise among the selected PTFs and between measured and PTF based results are again inconsistent and show no clear relationship with soil texture. Finally, as a third step, effective SHPs were calculated by using spatially averaged texture as PTF input representing areal average behaviour. For these effective SHPs the calculated effective values for the three selected functional variables appear to be close to the areally averaged values obtained with step 2. The selected functional variables thus appear to depend linearly on the PTFs over the range for which the data are representative. This suggests that for our specific PTFs areal mean or effective values for the three functional variables can be obtained fairly accurately from a single measurement of a bulk collection of soil samples as input. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

15.
Pedotransfer functions (PTFs) have become a topic drawing increasing interest within the field of soil and environmental research because they can provide important soil physical data at relatively low cost. Few studies, however, explore which contributions PTFs can make to land‐use planning, in terms of examining the expected outcome of certain changes in soil and water management practices. This paper describes three scenario studies that show some aspects of how PTFs may help improve decision making about land management practices. We use an exploratory research approach using simulation modelling to explore the potential effect of alternative solutions in land management. We: (i) evaluate benefits and risks when irrigating a field, and the impact of soil heterogeneity; (ii) examine which changes can be expected (in terms of soil water balance and supply) if organic matter content is changed as a result of an alternative management system; (iii) evaluate the risk of leaching to deeper horizons in some soils of Hungary. Using this research approach, quantitative answers are provided to ‘what if?’ type questions, allowing the distinction of trends and potential problems, which may contribute to the development of sustainable management systems. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

16.
Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten–Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.  相似文献   

17.
Measurement uncertainty is a key hindrance to the quantification of water fluxes at all scales of investigation. Predictions of soil‐water flux rely on accurate or representative measurements of hydraulic gradients and field‐state hydraulic conductivity. We quantified the potential magnitude of errors associated with the parameters and variables used directly and indirectly within the Darcy – Buckingham soil‐water‐flux equation. These potential errors were applied to a field hydrometric data set collected from a forested hillslope in central Singapore, and their effect on flow pathway predictions was assessed. Potential errors in the hydraulic gradient calculations were small, approximately one order of magnitude less than the absolute magnitude of the hydraulic gradients. However, errors associated with field‐state hydraulic conductivity derivation were very large. Borehole (Guelph permeameter) and core‐based (Talsma ring permeameter) techniques were used to measure field‐saturated hydraulic conductivity. Measurements using these two approaches differed by up to 3\9 orders of magnitude, with the difference becoming increasingly marked within the B horizon. The sensitivity of the shape of the predicted unsaturated hydraulic conductivity curve to ±5% moisture content error on the moisture release curve was also assessed. Applied moisture release curve error resulted in hydraulic conductivity predictions of less than ±0\2 orders of magnitude deviation from the apparent conductivity. The flow pathways derived from the borehole saturated hydraulic conductivity approach suggested a dominant near‐surface flow pathway, whereas pathways calculated from the core‐based measurements indicated vertical percolation to depth. Direct tracer evidence supported the latter flow pathway, although tracer velocities were approximately two orders of magnitude smaller than the Darcy predictions. We conclude that saturated hydraulic conductivity is the critical hillslope hydrological parameter, and there is an urgent need to address the issues regarding its measurement further. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

18.
Antecedent soil moisture significantly influenced the hydraulic conductivity of the A1, A2e and B21 horizons in a series of strong texture‐contrast soils. Tension infiltration at six supply potentials demonstrated that in the A1 horizon, hydraulic conductivity was significantly lower in the ‘wet’ treatment than in the ‘dry’ treatment. However in the A2e horizon, micropore and mesopore hydraulic conductivity was lower in the ‘dry’ treatment than the ‘wet’ treatment, which was attributed to the precipitation of soluble amorphous silica. In the B21 horizon, desiccation of vertic clays resulted in the formation of shrinkage cracks which significantly increased near‐saturated hydraulic conductivity and prevented the development of subsurface lateral flow in the ‘dry’ treatment. In the ‘wet’ treatment, the difference between the hydraulic conductivity of the A1 and B21 horizons was reduced; however, lateral flow still occurred in the A1 horizon due to difficulty displacing existing soil water further down the soil profile. Results demonstrate the need to account for temporal variation in soil porosity and hydraulic conductivity in soil‐water model conceptualisation and parameterisation. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

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

20.
One‐dimensional flow simulations were conducted at four locations of the shallow alluvial aquifer of the upper Rhine River (at the Erstein polder) to quantify the time‐dependent moisture distribution, the water flux and the water volume infiltrated in the unsaturated zone as a function of soil heterogeneities during a five‐day‐long flooding event. Three methods of estimating the hydraulic parameters of soil in the vadose zone were tested. They are based on the following: (1) experimental data, (2) soil particle‐size distribution and (3) pedology information on soils. Water fluxes calculated from modelling approaches 2 and 3 were compared with those of the experiment‐based values and the effect of these differences on the arrival time and velocity of water at the water table were analysed. Major differences in water fluxes were found among the methods of estimating the hydrodynamic parameters. At the Terrace location, the groundwater recharge predicted using soil data from methods 1 and 2 are approximately 4500 and 2400 mm, respectively. Flow simulations using soil data and the experiment‐based method show the highest velocities of infiltrating water at the soil surface and largest volume of groundwater infiltration but result in the lowest centres of the moisture content mass. The results obtained using soil data based on the pedological method are similar to those calculated using soil parameters based on the particle‐size distribution of extracted soil samples. Water pressure profiles calculated on Terrace and Channel location, 3 and 7 days after the inundation event agreed reasonably well with those observed when using hydrodynamic parameters from the experiment‐based method. However, the flow model using the pedology‐based parameters largely underestimates the time needed to achieve hydrostatic conditions of the soil water profile once water flooding at the soil surface stops. This can be mainly attributed to the low values of estimated van Genuchten parameter α. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号