Evaluation of porous medium hydraulic properties using experimental methods and RETC code

Two experimental procedures were used to determine both hydraulic properties, soil water retention θ(h) curve and unsaturated hydraulic conductivity K(θ), of a sand sample. Knowledge of hydraulic properties is essential, since they generally control soil water dynamics. A steady-state laboratory method was used for the simultaneous determination of θ(h) and K(θ). A one-step outflow method was used for the determination of diffusivity D(θ) and subsequently K(θ) from soil water retention data which were measured independently on the same sample and using the same apparatus. The comparison of K(θ) measured values from the above-mentioned methods showed very good agreement of the results. Also, the comparison between the experimental K(θ) and θ(h) functions and the predictions obtained using retention curve (RETC) code by simultaneous fit of experimental soil water retention and hydraulic conductivity data from outflow data, assuming the Mualem-van Genuchten model, showed very good agreement. It is noted that the main disadvantage of the one-step outflow method is the weakness to predict K(θ) values near saturation. This disadvantage could be overcome using RETC code with the above procedures, since the K(θ) values between the predictive approach and the steady-state method were similar.     

Determination of chemical transport properties for different textures of undisturbed soils

Agrichemicals usually contaminate groundwater via preferential flow, therefore determination of the preferential flow characteristics of soil is needed. One model that predicts solute transport due to preferential flow is the mobile–immobile (MIM) solute-transport model, which partitions total water content (θ; m³ m^?3) into mobile (θ_m) and immobile fractions (θ_im). In undisturbed soils, a method is proposed for determining the MIM model parameters, i.e. immobile water fraction (θ_im), mass transfer coefficient (α) and hydrodynamic dispersion coefficient (D _h). Breakthrough curves were obtained for five different soil textures in three replicates, by miscible displacement of Cl^? in undisturbed soil columns. Cl^? breakthrough curves were evaluated in terms of the MIM model. Analysis suggests that the values of D _h and α increased with lighter soil textures and θ_im increased with heavier soil textures. The values of θ_im ranged from 5.31 to 14.28% in different soil textures. Furthermore, values of θ_im were found to be related to soil clay content. Values of α ranged from 0.0257 to 0.32 h^?1 and values of D _h ranged from 0.36 to 11.2 cm² h^?1 in different soil textures. A significant linear correlation was obtained between α, θ_im, D _h and soil saturated hydraulic conductivity (K _s) and pore water velocity (v). A multivariate pedotransfer function was developed to estimate α, θ_im and D _h based on the geometric mean (d _g) and the standard deviation (σ_g) of the diameter of soil particles and soil organic matter content. The pedotransfer functions for D _h, θ_im and α were validated by independent data sets from other investigators.     

陕北水蚀风蚀交错区两种生物结皮对土壤饱和导水率的影响

该文以陕北水蚀风蚀交错区普遍发育的地表和地上两种生物结皮为研究对象，分别以3种非生物结皮（无结皮、物理结皮、去除生物结皮）为对照，使用盘式入渗仪测定其饱和导水率。结果表明：与无结皮土壤相比，两种类型生物结皮均可极显著降低土壤饱和导水率；与去除生物结皮土壤相比，两种类型生物结皮对土壤饱和导水率的降低均不显著；与有物理结皮发育的土壤相比，地表生物结皮对土壤饱和导水率的降低不显著，而地上生物结皮对土壤饱和导水率的降低显著。一方面，两种生物结皮对土壤饱和导水率均有明显降低作用，预示生物结皮在降雨活动中可能会增加径流、降低入渗，阻碍研究区水分亏缺条件下的植被恢复和生态与环境建设。另一方面，与不同的对照相比，生物结皮对土壤饱和导水率的影响截然不同，该结论可在一定程度上解释当前有关生物结皮影响土壤水分入渗方面所存在的分歧。     

Spatio-temporal variation of anisotropy of saturated hydraulic conductivity in a tilled sandy loam soil

Knowledge on anisotropy of saturated hydraulic conductivity can improve the understanding of transport phenomena in soil. We hypothesized that saturated hydraulic conductivity (K_s) in the upper part of the root zone of an agricultural sandy loam soil was anisotropic at different soil depths and times after tillage. K_s was measured on undisturbed 100 cm³ core samples taken in the horizontal and vertical directions in up to four soil layers (Surf: surface layer (0–5 cm); Top: topsoil (10–15 cm); Trans: transition layer between topsoil and subsoil; Sub: subsoil (40–60 cm)) 1, 8 and 32 months, respectively, after mouldboard ploughing and drilling. The ratio between estimated geometric mean values for K_s in the vertical and the horizontal directions (K_ms,v/K_ms,h) was used to test the hypotheses. A total of 669 soil samples were analysed.K_ms,v/K_ms,h varied with time after tillage and between soil layers. One month after ploughing, K_ms,v/K_ms,h was <0.23 (P = 0.975) in the Trans layer with an average value of 0.084, i.e. K_ms,h was 12 times larger than K_ms,v. Anisotropy could not be documented in this layer 8 or 32 months after ploughing, i.e. K_ms,v/K_ms,h was not significantly different from 1.0. For the Surf and Top layers 32 months after ploughing, K_ms,v/K_ms,h was in the intervals 1.4–50 and 3.1–77, respectively, (P = 0.95) with average values of 8.4 and 15, respectively. Thus, K_ms,v was 8.4 respectively 15 times larger than K_ms,h in the two layers. Anisotropy was not found in these layers 1 or 8 months after tillage. Strong anisotropy was found in the Sub layer with K_ms,v/K_ms,h averaging to 14 and 32, respectively, 8 and 32 months after tillage. K_ms,v and K_ms,h generally decreased with time in the Surf, Top and Trans layers, except in the vertical direction in the Top layer between 8 and 32 months after ploughing, and in the Trans layer between 1 and 8 months after ploughing. Overall, the geometric means of K_s varied between 10^−4.0 and 10^−7.1 m s⁻¹.The results may reflect systematic effects of soil settlement and drying/wetting phenomena coupled with biological activity and the existence of stable, vertically oriented biopores in the subsoil. It appears to be necessary to consider anisotropy of K_s and its variation in the analysis and modelling of water flow and chemical transport in agricultural soils, particularly to explain heterogeneous flow phenomena at the plot and field scales.     

Mualem模型中的饱和导水率修正研究

[目的]提高Mualem模型计算非饱和导水率的准确性。[方法]采用理论推导结合数据统计的方法研究该模型中的饱和导水率ks修正问题。基于Brooks—Corey土壤水分特征曲线模型,建立修正导水率ko与土壤水分特征曲线之间的理论关系式;通过回归分析得到ko与土壤水分特征曲线之间的理论关系式中相关参数。[结果]利用原状黄土的非饱和渗透试验数据,对考虑修正导水率ko的Mualem模型的准确性进行了验证,得到了比较满意的结果。[结论]研究成果可用于依据土壤水分特征曲线直接确定非饱和土导水率,对非饱和导水率预测研究具有一定参考价值。     

黄土丘陵区典型植被土壤物理性质差异及其对导水特性影响

选取黄土丘陵区12种典型植被样地,通过测定各样地不同土层植物残体生物量、土壤容重、毛管孔隙度、非毛管孔隙度及饱和导水率,研究各指标随土层深度和植被类型的变化规律及其对土壤饱和导水率的影响。结果表明:(1)除容重随土层深度增加外,植物残体、毛管孔隙度、非毛管孔隙度和饱和导水率均随土层深度减少,其中植物残体大多集中于表层土壤(0—10 cm),占总残体生物量的51.4%～85.7%。(2)不同植被类型其植物残体及土壤物理性质存在显著差异,乔木林地植物残体、农耕地土壤容重、灌木林地非毛管孔隙度及饱和导水率均最大,而毛管孔隙度与不同土地利用类型间无显著差异。(3)饱和导水率随植物残体生物量密度(0—10 cm)和土壤容重呈幂函数减小,随毛管孔隙度和非毛管孔隙度呈幂函数增大;土壤容重(BD)和非毛管孔隙度(NCP)是影响土壤饱和导水率(K_s)的主要因素,且土壤饱和导水率可表示为两者的综合非线性方程(K_s=0.6BD~(-4.717)NCP~(0.203),P0.01,R~2=0.63,NSE=0.50)。此外,沙棘灌木林地平均饱和导水率最大,有利于降雨过程中土壤水分入渗,具有较强的水土保持功能。本研究结果可为黄土高原植被恢复生态水文效益评价提供理论依据。     

Pedotransfer Functions for the Estimation of Saturated Hydraulic Conductivity for Some Indian Sandy Soils

Determination of the saturated hydraulic conductivity (k_s) is needed in many studies and applications related to irrigation, drainage, water movement and solute transport in the soil. Although many advances are made for direct measurements of k_s, they are usually time consuming and costly. Some attempts have been made to indirectly predict the saturated hydraulic conductivity from the more easily or readily available basic soil properties. The objective of this study was to develop and validate Pedotransfer Functions (PTFs) for estimation of saturated hydraulic conductivity using multiple non-linear regression technique. One hundred and one soil samples were collected from agricultural and forest soils at different depths, at different locations in the Pavanje River basin that lies in the southern coastal region of Karnataka, India. Saturated hydraulic conductivity was measured, by variable falling head method through Permeameter in the laboratory. Prediction accuracies were evaluated using coefficient of determination, root mean square error, mean error, geometric mean error ratio and geometric standard deviation of the error ratio between measured and predicted values. The results show that, the PTFs for the estimation of saturated hydraulic conductivity could be used appropriately for the soils with loamy sand and sandy loam textures falling in this area of the coastal region of southern India.     

Estimation of saturated hydraulic conductivity from double‐ring infiltrometer measurements

This research aims to determine soil vertical saturated hydraulic conductivity (K_s) in situ from the measured steady infiltration rate (I), initial soil properties and double‐ring infiltrometer (DRI) test data. Characterizing the effects of these variables on the measured steady infiltration rate will enable more accurate prediction of K_s. We measured the effects of the ring diameter, head of ponding, ring depth, initial effective saturation and soil macroscopic capillary length on measured steady infiltration rates. We did this by simulating 864 DRI tests with the finite element program HYDRUS‐2D and by conducting 39 full‐scale in situ DRI tests, 30 Mini‐Disk infiltrometer experiments and four Guelph Permeameter tests. The M5′ model trees and genetic programming (GP) methods were applied to the data to establish formulae to predict the K_s of sandy to sandy‐clay soils. The nine field DRI tests were used to verify the computer models. We determined the accuracy of the methods with 30% of the simulated DRI data to compare I/K_S values of the finite element models with estimates from the suggested formulae. We also used the suggested formulae to predict the K_s values of 30 field DRI experiments and compared them with values measured by Guelph Permeameter tests. Compared with the GP method, the M5′ model was better at predicting K_S, with a correlation coefficient of 0.862 and root mean square error (RMSE) of 0.282 cm s^?1. In addition, the latter method estimated K_svalues of the field experiments more accurately, with an RMSE of 0.00346 cm s^?1.     

Variability in Soil Hydraulic Conductivity and Soil Hydrological Response Under Different Land Covers in the Mountainous Area of the Heihe River Watershed,Northwest China

Understanding the variability in soil hydraulic conductivity in the mountainous headwaters is critical to the modeling of mountainous runoff and the water resources management of river basins, especially in the arid and semiarid areas. In this study, a total of 32 soil profiles with five layers within 0–70 cm were sampled under different land cover types: forest, meadow, high coverage grassland (HCG), medium coverage grassland (MCG), and barren land in the upper stream of the Heihe river watershed, Northwest China. Saturated hydraulic conductivity (K _S) was measured for each sample. The vertical variation of K _S and soil hydrological response under different land covers were analyzed. Results show that K _S value in layer 5 was significantly lower than the values of above four layers. K _S decreased in the order of forest, meadow, HCG, MCG, and barren land, corresponding to the degree of vegetation degradation. The K _S decreased with depth under forest, HCG, and barren land, but increased first and then decreased under meadow and MCG. The dominant stormflow paths for different land covers were different: forest was dominated by deep percolation, HCG was dominated by subsurface flow (SSF), meadow was prevailed by Hortonian overland flow and had no SSF, while MCG and barren land were also dominated by Hortonian overland flow, but still formed SSF. This result provides important information for improving the accuracy of mountainous hydrological modeling and, in turn, leading to sustainable management of water resources in the study watershed. Copyright © 2016 John Wiley & Sons, Ltd.     

Prüfung von 3 Verfahren zur Vorhersage der hydraulischen Leitfähigkeit ungesättigter Böden aus Wasserretentionsdaten oder aus der Bodenart

Testing of three methods to predict unsaturated soil hydraulic conductivity from water retention data or from texture class Using 60 soils taken from UNSODA (Leij et al., 1996) the method proposed by Renger et al. (1999) and the prediction according to Mualem (1976)/van Genuchten (1980) were tested. The parameters of the Mualem/vanGenuchten model were estimated either from water retention data or from a table of reference values. Using the reference values requires only the knowledge of texture class (German classification system). An advantage of the method proposed by Renger et al. (1999) is its capability to predict saturated conductivity too. The model of Mualem/vanGenuchten using reference values of parameters yields the best results. The standard deviation between observed and predicted values of unsaturated hydraulic conductivity was 0.93 lg (K) for the Mualem/vanGenuchten model and 1.3 lg (K) for the Renger et al. (1999) predictions.     

Prediction of saturated hydraulic conductivity of compacted soils using empirical scaling factors

A new empirical-based scaling method is introduced to predict saturated hydraulic conductivity (K _s ) of compacted soils. This method is an improvement of the former non-similar media concept (NSMC) model that is generalized for tilled and untilled conditions. In this method, geometric mean particle size diameter (d_g ), geometric standard deviation (σ _g ) and saturated soil water content (total porosity) are successfully incorporated in the empirical-based scaling factor of K _s . Results showed that the scaled model overestimated K _s by ～18%, whereas the NSMC model underestimated K _s by ～21%. However, the scaled model based on the similar media concept (SMC) failed to predict K _s . Because of the complexity and high uncertainty in determining the shape factor parameter in the NSMC model, it is suggested that the new scaled model might be used reliably in practical cases to predict K _s in the various layers of compacted soils irrespective of the tillage condition. Further assessment of the new scaling model in other areas, in which new collected data are available, is recommended.     

Measurement of solid phase permittivity for volcanic soils by time domain reflectometry

Dielectric properties of soils are widely used to estimate their water content. Andisols are unique soils in terms of aggregate structure and show dielectric properties different from other mineral soils. To understand the relationship between the dielectric properties and soil physical properties, multi‐phase dielectric mixing models are often used. However, solid phase permittivity (?_s) for Andisols, which is required for calculating the model output, has not been directly measured yet. Therefore, the objectives for this study were to measure ?_s for Andisols in Japan. In addition, the measured ?_s values were compared with those calculated from the traditional method, which applies two‐phase mixing models to estimate ?_s values based on effective permittivities of repacked soil samples. The applicability of the traditional method to estimate such values for Andisols is also discussed. The effects of bound water and aggregate structure on measured ?_s values were evaluated prior to their measurements. We found that the aggregate structural effects were negligible. However, the amount of bound water caused overestimates of ?_s. Four Andisols from the A and B horizons of soils in Memuro Hokkaido (northern part of Japan), from the A horizon of a soil in Tsukuba Ibaraki (central part of Japan), and from the A horizon of a soil in Koshi Kumamoto (southern part of Japan) were used in this study. The ?_s values obtained fall between 5.6 and 6.1, and deviated from the estimated ?_s values derived from the traditional method. Therefore, the traditional method is probably unsuitable to estimate ?_s values for Andisols.     

Ploughing frequency and compost application effects on soil infiltrability in a cotton–maize (Gossypium hirsutum–Zea mays L.) rotation system on a Ferric Luvisol and a Ferric Lixisol in Burkina Faso

One of the key issues to increase soil productivity in the Sahel is to ensure water infiltration and storage in the soil. We hypothesised that reducing tillage from annual to biennial ploughing and the use of organic matter, like compost, would better sustain soil hydraulic properties. The study had the objective to propose sustainable soil fertility management techniques in the cotton–maize cropping systems. The effects of reduced tillage (RT) and annual ploughing (AP) combined with compost application (Co) on soil infiltration parameters were assessed on two soil types. Topsoil mean saturated hydraulic conductivities (K_s) were between 9 and 48 mm h⁻¹ in the Luvisol, while in the Lixisol they were between 18 and 275 mm h⁻¹. In the two soil types compost additions with reduced tillage or with annual ploughing had the largest effect on K_s. Soil hydraulic behaviour was in reasonable agreement with soil pore size distribution (mean values varied from 19.5 to 237 μm) modified by tillage frequency and organo-mineral fertilization. Already the first 3 years of this study showed that use of organic matter, improved soil infiltration characteristics when annual ploughing was used. Also biennial ploughing showed promising results and may be a useful strategy for smallholders to manage these soils.     

Effect of zeolite on saturated hydraulic conductivity and crack behavior of silty clay paddled soil

The addition of zeolites to soil modifies soil physical and chemical properties. The objectives of this research were to study the effect of zeolite on saturated hydraulic conductivity, K _s, and crack behavior in a silty clay paddled soil. Soil samples were mixed with 0, 4, 8 and 12 g kg^?1 of zeolite for K _s determination, and 0, 2, 4, 8 and 12 g kg^?1 for soil crack measurements. Saturated hydraulic conductivity was measured using the constant head method. The results indicated that K _s was significantly increased at a zeolite application rate of 8 g kg^?1. Furthermore, an increase in zeolite content up to 8 g kg^?1 decreased soil crack area after paddling and first rewetting. Higher amounts of zeolite (e.g. 12 g kg^?1) increased crack density after the second rewetting. However, a 50% reduction in crack depth occurred with zeolite application rates of 8 and 12 g kg^?1 in comparison with controls. Thus, a zeolite application rate of 8 g kg^?1 is recommended for soil crack reduction in intermittent-flood irrigation. Furthermore, a relationship was obtained between crack area density (Ln), gravimetric soil water content and zeolite application rate. After the second irrigation, a relationship was also obtained between crack depth, gravimetric soil water content and zeolite application rate. Crack depth showed a positive and highly significant linear correlation with crack width.     

Spatial variability of water retention parameters and saturated hydraulic conductivity in a calcareous Inceptisols (Khuzestan province of Iran) under sugarcane cropping

The objective of this study was to quantify inherent spatial variability and spatial cross-correlation of the van Genuchten retention parameters and saturated hydraulic conductivity (K_s) of surface and subsurface layers in a calcareous Inceptisols (Khuzestan province, Iran) under sugarcane cropping. Measurements were performed on 100-cm³ undisturbed soil cores collected at 94 locations along a 30-m-long transect with horizontal sampling distance intervals of 0.3 and 1 m at soil depths of 0–40 and 40–80 cm, respectively. Spatial variability was investigated using conventional statistics and geostatistical techniques. Coefficient of variation (CV) varied from 8.2% (for shape parameter, n at 40–80 cm depth) to 256.7% (for K_s at 0–40 cm depth). The n parameter and saturated water content, θ_s, showed a small-scale spatial heterogeneity with a maximum CV of 11.3% for the first depth and 9.2% for the second depth. Most of the hydraulic parameters at both depths showed a spatial structure and convex experimental semivariograms with dominant spherical models with the influence range of 3.2–41 m. In most cases, the extent of spatial correlation scales of cross-semivariograms for pairs of cross-correlated hydraulic variables was found to be different with reference to those relating to the direct semivariograms of correlated variables.     

渝东北紫色土饱和导水率传递函数研究

为研究渝东北紫色土理化性质在垂直空间上的分布情况以及对饱和导水率的影响,进而建立饱和导水率与各理化性质间的关系函数,推求饱和导水率的传递函数,选择渝东北开州区、云阳县等7个区县内45个紫色土典型田块为研究区域,运用Excel 2013和Matlab 2015b软件统计分析后,利用多元非线性回归法推求并验证了渝东北紫色土饱和导水率传递函数模型和模型参数。研究表明:①研究区土壤饱和导水率变化范围在0.16～195.68 cm/d,变化范围广,空间变异系数大,变异性较强;同一采样点深度越大,饱和导水率越小;②土壤饱和导水率与有机质含量有显著的指数函数关系,与饱和含水量有较强的二次函数关系,与土壤容重和土壤颗粒的相关性不大;③本次试验建立的土壤饱和导水率传递函数模型及模型系数检验合格,预测值与实际测算值误差较小,精度良好,可用于渝东北紫色土饱和导水率的预测工作。     

青海省东部不同土地利用方式下土壤饱和导水率分布及其影响因素

研究高寒地区土壤饱和导水率分布特征及其影响因素可为评估脆弱生态系统水源涵养能力和构建区域水文模型提供参数.通过测定青海省东部南北样线24个样点(0—30 cm)土壤基本理化性质和饱和导水率(Ks),分析了不同土地利用方式下Ks分布特征及其影响因素.结果表明:Ks均值表现为林地(1.89 cm/h)>草地(1.62 cm...     

Estimation of the diffusion coefficients of adsorbed and non-adsorbed solutes in soil

A method is proposed which follows Darrah's experimental procedure and takes advantage of a mathematical solution provided by Carslaw & Jaeger to estimate the diffusion coefficients of adsorbed and non-adsorbed solutes in soil. The method requires only the values of the concentration of the solute at the input face of a uniform column of soil, C^s, and of the total amount, Q_t, that has entered the soil after a specified time during which the surface of the block is in contact with a thin porous pad containing a known initial amount of solute, Q₀, at concentration C₀, expressed in the same units as C^s. In the C^s/C₀ vs. Q_t/Q₀ space there is a unique relationship between the effective diffusion coefficient, D_e, of the solute in the soil and the contact conductance for this solute, h, between the pad and the soil surface. The proposed procedure is firstly to determine D_e, and h for a non-adsorbed solute in the experimental soil using the experimental values of C^s/C₀ and Q/_Q for that solute. This value of D_e, gives the diffusion impedance factor for the solute in the soil, f, which is assumed also to apply to adsorbed solutes. A first estimate of the effective diffusion coefficient of an adsorbed solute, ¹D_ea, is then made using f and the diffusion coefficient of the free solute in water, D_L, obtained from the literature (i.e. ¹D_ea= D_Lf). Only if the solute is weakly adsorbed will the values of C^s/C₀, and Q_t/Q₀ lie in C^s/C₀, vs. Q_t/Q₀, space as defined by ¹D_ea and the contact conductance, h. Instead a second space relating C^s/C₀ and Q_t/Q₀, is now constructed from nominated values of h and D_e, where D_e, is defined in terms of ¹D_ea, the adsorption coefficient, F , and the volumetric moisture content of the soil, θ. The position of the experimental values of C^s/C₀, and Q_t/Q₀ within this new space defines h and the actual D_e, and F of the solute as it diffuses and is adsorbed in the soil. The advantages and limitations of the method are discussed. In particular, the method assumes that the adsorption process is linear and reversible.     

Exploring the variation in soil saturated hydraulic conductivity under a tropical rainforest using the wavelet transform

Saturated hydraulic conductivity (K_s) of the soil is a key variable in the water cycle. For the humid tropics, information about spatial scales of K_s and their relation to soil types deduced from soil map units is of interest, as soil maps are often the only available data source for modelling. We examined the influence of soil map units on the mean and variation in K_s along a transect in a tropical rainforest using undisturbed soil cores at 0–6 and 6–12 cm depth. The K_s means were estimated with a linear mixed model fitted by residual maximum likelihood (REML), and the spatial variation in K_s was investigated with the maximum overlap discrete wavelet packet transform (MODWPT). The mean values of K_s did not differ between soil map units. The best wavelet packet basis for K_s at 0–6 cm showed stationarity at high frequencies, suggesting uniform small‐scale influences such as bioturbation. There were substantial contributions to wavelet packet variance over the range of spatial frequencies and a pronounced low frequency peak corresponding approximately to the scale of soil map units. However, in the relevant frequency intervals no significant changes in wavelet packet variance were detected. We conclude that near‐surface K_s is not dominated by static, soil‐inherent properties for the examined range of soils. Several indicators from the wavelet packet analysis hint at the more dominant dynamic influence of biotic processes, which should be kept in mind when modelling soil hydraulic properties on the basis of soil maps.     

重庆市四面山不同土地利用类型饱和导水率

[目的]探讨不同土地利用类型和土壤理化性质对饱对导水率的影响。[方法]采用定水头法测定四面山不同土地利用类型的饱和导水率,并运用回归分析,相关分析和主成分分析法分析其与土壤物理因子和有机质的关系,以及影响饱和导水率的主导因子。[结果]各土地利用类型的平均饱和导水率均高于荒地,其顺序为:林地农地草地,林地中天然林饱和导水率大于人工林;饱和导水率随土层深度的增加呈负指数递减规律;饱和导水率与容重呈幂函数关系,与孔隙度呈正相关,与黏粒含量呈负相关;有机质含量的提高对饱和导水率有积极的促进作用。[结论]影响饱和导水率的主导因子是容重、有机质、非毛管孔隙度和毛管孔隙度,其次,土壤机械组成对其也有一定影响。