科尔沁沙地土壤水分特征曲线传递函数的构建与评估

为了快速准确的获取某一区域的水力性质,该文以科尔沁沙地典型沙丘-草甸相间地区为研究区,在对该区49个不同地貌类型采样点土壤水分特征曲线与土壤基本物理化学特性参数测试分析的基础上,采用函数参数非线性规划法构建了土壤干容重、粒径分布、有机质、pH值、电导率值等基本参数与水分特征曲线之间的传递函数,并进行了精度评估与分析。结果表明:1)研究区土壤水分特征曲线的坡度陡峭,不同地貌类型与土地利用方式下,土壤水分特征曲线有较大差异,相同负压下,土壤持水量按照流动沙丘-半固定沙丘-固定沙丘-农田-草甸的顺序递增;土壤的供水能力按照流动沙丘-半固定沙丘-固定沙丘-草甸-农田的顺序递减;2)利用土壤基本物化特性参数通过函数参数非线性规划法建立了研究区土壤水分特征曲线的传递函数,干容重和砂粒含量是预测土壤水分特征曲线模型参数的主要变量,增加土壤的理化指标可以提高预测精度,然而有机质含量、pH值、电导率值对本区土壤水分常数的影响并不大,对水分特征曲线模型3个参数的影响略微增加;3)通过对传递函数的检验与精度评估分析,各参数的平均误差均在0附近;饱和含水率、残余含水率的均方根误差分别为0.017、0.023;土壤水分常数的相关系数在0.95附近,饱和含水率、残余含水率的误差比的几何标准偏分别为1.04、1.27。表明所建土壤水分特征曲线传递函数的精度较高,可用于该区土壤水分特性研究。该研究可为该区水分、溶质运移、水-热-盐耦合运移模拟提供技术支持和理论保证。     

轻质土壤水分特征曲线估计的简便方法

以黄淮海平原封丘地区的潮土和风沙土为研究对象,根据大量的土壤基本物理性质和土壤持水数据,利用多元逐步回归分析方法,建立了轻质土壤在不同基质势下土壤含水量(θ)的传递函数模型,并进行了模型验证。结果表明,利用轻质土壤的基本物理性质估计其水分特征曲线是一种简便可行的方法,并且,在回归方程中,增加-30kPa含水量项可提高-30kPa以上土壤含水量的估计精度;增加-1500kPa含水量项可以明显提高-100kPa至-1500kPa间土壤含水量的估计精度。     

河北省土壤容重的传递函数研究

由于非破坏性取样测定土壤容重是一项困难、昂贵、费时和不符合实际的工作。为克服这些困难,国外学者已经研究出多种与土壤有机质或有机碳含量和土壤质地有关的传递函数去估计其研究范围内的土壤容重。本文利用河北省第二次土壤普查资料中的土壤有机质和质地数据,对已有的13种土壤传递函数模拟的结果进行对比,确定最佳传递函数形式后,用SPSS统计分析软件进行回归分析,求得适宜河北省土壤容重的传递函数表达式。     

中国主要土壤类型的土壤容重传递函数研究

由于我国大多数土壤数据库缺失或部分缺失土壤容重数据,利用土壤其他属性来预测土壤容重具有重要意义。我国土壤类型多样,已有的土壤容重传递函数在不同土壤类型中的适用性也值得进一步探讨。本文基于我国现有的土壤数据库,对已有的两种土壤容重传递函数在不同土壤类型(土纲)中的预测精度与适用性进行评估,最后通过SPSS进行回归分析建立我国主要土壤类型最适宜的容重传递函数。研究结果表明,已有的两种土壤容重传递函数应用于部分土壤类型时预测精度不高,但基于土壤系统分类的数据分组后建立的容重传递函数能够明显提高预测精度。新建的容重传递函数对有机土、铁铝土、潜育土、均腐土、富铁土、淋溶土、雏形土、新成土和变性土土壤容重的预测精度较高,但人为土、盐成土和干旱土土壤容重传递函数的预测精度整体偏低,在应用时需要慎重。在利用土壤容重传递函数时一定要注意研究区及适用范围。     

灰色关联及非线性规划法构建传递函数估算黑土水力参数

土壤水分特征曲线和饱和导水率是重要的水力参数,为了简便准确获取这些参数,以松嫩平原黑土区南部为研究区域,采集136个采样点土样用于测定不同土层土壤水分特征曲线、饱和导水率以及土壤理化性质,并运用灰色关联分析确定影响土壤水力参数的主要土壤理化性质,采用非线性规划构建土壤分形维数、有机质、干容重、土壤颗粒组成与土壤水分特征曲线、饱和导水率之间的土壤传递函数,并通过与现有土壤传递函数对比分析进行精度验证。结果表明:1)土壤分形维数是估算土壤水分特征曲线模型参数和饱和导水率的主要参数之一,同时,干容重和有机质含量也在不同土层土壤传递函数中起到重要的作用;2)通过验证分析,不同土层各参数平均绝对误差接近于0,均方根误差值也都较小,其中在不同土层土壤传递函数估算的土壤含水率均方根误差分别为0.022、0.017cm~3/cm~3;3)对比分析其他已存的土壤水分特征曲线和饱和导水率的土壤传递函数,该文构建的土壤传递函数均方根误差值均较小,决定系数值都在0.66以上,表明估算精度较高,均好于其他方法估算精度,具有良好的区域适应性。综上,所构建的土壤水分特征曲线和饱和导水率土壤传递函数可以用于松嫩平原黑土区土壤水力参数估算。     

科尔沁沙丘-草甸相间地区表土饱和导水率的土壤传递函数研究

以科尔沁沙地典型坨-甸相间地区为研究区,野外布设240个采样点,对流动沙丘、半固定沙丘、固定沙丘、沙丘区杨树林、沙丘区耕地、低覆盖度草甸、高覆盖度草甸、草甸区耕地、撂荒地9种地貌类型下的表层土壤进行了采样,测定了其含水率、干容重、有机质、饱和导水率等理化特性,分析了不同地貌类型下表层土壤理化参数差异。选取Campbell、Cosby、Wosten等、Saxton等4种土壤饱和导水率传递函数,对该地区表土饱和导水率进行了预测。结果显示这几种土壤传递函数预测值与实测值偏差较大,相关系数均小于0.3,精度难以满足本地区应用。在此基础上,选取土壤容重、有机质含量、饱和含水率、平均粒径、粒径标准偏差5种土壤特性参数作为输入变量,采用主成分分析与非线性回归分析相结合的方法,重新建立了预测本地区表土饱和导水率的土壤传递函数,结果显示预测值与实测值相关系数为0.661,该传递函数可用于科尔沁沙地表层土壤饱和导水率的预测。     

黑土区水蚀坡耕地土壤穿透阻力时空变异特征及传递函数

[目的] 为探明侵蚀坡耕地黑土硬度的时空特征及关键驱动因素。[方法] 基于典型水蚀坡耕地110个样点土壤穿透阻力动态监测及其相关因素的测定,在阐明土壤穿透阻力时空特征的基础上,利用多元线性回归(MLR)和随机森林(RFR)模型分析土壤穿透阻力的影响因素,并进一步构建其传递函数。[结果] 土壤穿透阻力时空变异性受土层深度、水分条件及农事活动等因素的综合影响,耕作层土壤穿透阻力空间上的异质性显著低于亚表层,变异系数分别为17.4%和26.3%。随土壤由湿变干和距离翻耕的时间增加,空间异质性呈增强趋势。土壤侵蚀影响土壤穿透阻力,强烈侵蚀区土壤穿透阻力高于沉积区,特别是在湿润条件下,差异更为显著(p<0.05)。与自然林地相比,所研究坡耕地73%样点的土壤穿透阻力增加,尤其在强侵蚀区。土壤含水量、容重和有机碳(SOC)是土壤穿透阻力的主要影响因素,Pearson系数分别为-0.69,0.58,-0.54,三者共同解释土壤穿透阻力88%的变异。RFR模型在土壤穿透阻力预测中表现优于MLR模型,在生长季均值的预测中,R²可达0.91,RMSE仅为91.2 kPa。[结论] 研究结果加深对侵蚀背景下黑土硬度空间分异特征的理解,为黑土压实管控提供理论依据。     

基于传递函数的土壤数据库缺失数据的填补研究

数据缺失在土壤调查研究中是一个非常普遍的现象,处理不当一定程度上会影响研究结果的可靠性。土壤转换函数(pedotransfer functions,PTFs)是简单、快速、大批量填补土壤数据库缺失信息的有效手段。但目前分析和厘定我国土壤数据库缺失数据特征的研究较少,针对土壤数据库缺失数据的填补方法也亟待规范。本文对我国第二次土壤普查数据库进行分析,探讨该数据库的数据缺失特征,并对数据缺失严重的土壤属性进行预测,以期为今后的土壤数据库缺失数据填补工作提供参考。总体来看,质地(砂粒、粉粒和黏粒含量)、pH、有机质、全氮、全磷、全钾是土壤普查中最基础的调查项目,这些土壤属性信息的完整性最好。有效磷、速效钾和阳离子交换量数据有一定的缺失。碱解氮、容重、砾石含量、各种类型氧化铁数据缺失严重。在填补缺失数据时,建议首先考虑模型的稳定性,尽量使用那些相对稳定且数据完整性好的土壤属性来预测缺失数据。我国第二次土壤普查数据库基本都缺少空间属性信息,在填补缺失数据时最好采用简单而相对稳定的回归模型。利用回归分析得到的土壤传递函数可以较好地实现容重、碱解氮和部分阳离子交换量缺失数据的填补工作。尽管如此,由于部分土壤属性信息有一定的时效性,应用传递函数时要注意数据源的历史背景。     

科尔沁沙丘-草甸相间地区表土饱和导水率的土壤传递函数研究

本研究以科尔沁沙地典型坨-甸相间地区为研究区,野外布设240个采样点,对流动沙丘、半固定沙丘、固定沙丘、沙丘区杨树林、沙丘区耕地、低覆盖度草甸、高覆盖度草甸、草甸区耕地、撂荒地九种地貌类型下的表层土壤进行了采样,测定了这些样点的含水率、干容重、有机质、饱和导水率等理化特性参数,分析了不同地貌类型下表层土壤理化参数差异。选取Campbell、Cosby、Wosten1997/1999、Saxton四种土壤传递函数,对该地区表土饱和导水率进行了预测,结果显示这几种土壤传递函数预测值与实测值偏差较大,相关系数均小于0.3,由此可见这几种传统的土壤传递函数在本地区应用具有一定的局限性。在此基础上,选取土壤容重、有机质含量、饱和含水率、平均粒径、粒径标准偏差五种土壤特性参数作为输入变量,采用逐步多元统计回归、主成分分析及非线性回归分析相结合的方法,重新建立了预测本地区表土饱和导水率的土壤传递函数,结果显示预测值与实测值相关系数为0.648,该传递函数可应用于科尔沁沙地表层土壤饱和导水率的预测。     

Methods for predicting the optimum and the range of soil water contents for tillage based on the water retention curve

Information is needed on the range of soil water contents for tillage. The objective of the work was to develop methods for the prediction of the soil water contents at which tillage may be done satisfactorily. Three water contents are considered: the lower (dry) limit, the optimum water content, and the upper (wet) limit. This paper makes a synthesis of published results from tillage and soil physics experiments and also includes some new experimental results. The effects of tillage are considered in relation to some “fixed points” including the lower plastic limit, field capacity and a new fixed point “the inflection point”. These considerations lead to methods for prediction of the lower (dry) tillage limit, the optimum water content, and the upper (wet) tillage limit in terms of the parameters of the van Genuchten equation for soil water retention. Predictions can be made in terms of soil composition through the use of pedotransfer functions for the parameters of the van Genuchten equation. The new methods will enable the effects of soil degradation and climate change on tillage work days to be estimated. The results are potentially mappable using geographic information systems.     

非饱和土壤水分运动与热力学函数关系初探

在不同的温度条件下,研究了土壤水势对水分运动的影响,结果表明,在同一温度条件下,提高土壤水势可增加土壤非饱和导水率,呈现黄绵土＞lou土,在相同的土壤含水量条件下,增加温度可提高土壤非饱和导水率,其导水率温度效应值（dk/kt)lou土＞黄绵土。土壤含水量一定时,随着相对偏摩尔自由能亦[△（△-/G）]和相对偏摩尔焓变[△（△-/H）]增大,土壤非饱和导水率也增大,并且呈出黄绵土＞lou土,拟合得出的相对偏摩尔自由能变[△（△-/G）]和相对偏摩尔焓变[△（△-/H）]与土壤非饱和导水率方程,具有较好的适应性。     

Applicability of site-specific pedotransfer functions and rosetta model for the estimation of dynamic soil hydraulic properties under different vegetation covers

Background, Aims, and Scope  During the last decades, different methods have been developed to determine soil hydraulic properties in the field and laboratory. These methodologies are frequently time-consuming and/or expensive. An indirect method, named Pedotransfer Functions (PTFs), was developed to predict soil hydraulic properties using other easily measurable soil (physical and chemical) parameters. This work evaluates the use of the PTFs included in the Rosetta model (Schaap et al. 2001) and compares them with PTFs obtained specifically for soils under two different vegetation covers. Methods  Rosetta software includes two basic types of pedotransfer functions (Class PTF and Continuous PTF), allowing the estimation of van Genuchten water retention parameters using limited (textural classes only) or more extensive (texture, bulk density and one or two water retention measurements) input data. We obtained water retention curves from undisturbed samples using the ‘sand box’ method for potentials between saturation and 20 kPa, and the pressure membrane method for potentials between 100 and 1500 kPa. Physical properties of sampled soils were used as input variables for the Rosetta model and to determine site-specific PTFs. Results  The Rosetta model accurately predicts water content at field capacity, but clearly underestimates it at saturation. Poor agreement between observed and estimated values in terms of root mean square error were obtained for the Rosetta model in comparison with specific PTFs. Discrepancies between both methods are comparable to results obtained by other authors. Conclusions  Site-specific PTFs predicted the van Genuchten parameters better than Rosetta model. Pedotransfers functions have been a useful tool to solve the water retention capacity for soils located in the southern Pyrenees, where the fine particle size and organic matter content are higher. The Rosetta model showed good predictions for the curve parameters, even though the uncertainty of the data predicted was higher than for the site-specific PTFs. Recommendations and Perspectives  The Rosetta model accurately predicts the retention curve parameters when the use is related with wide soil types; nevertheless, if we want to obtain good predictors using a homogenous soil database, specific PTFs are required. ESS-Submission Editor: Prof. Zhihong Xu, PhD (zhihong.xu@griffith.edu.au)     

土壤持水曲线van Genuchten模型求参的Matlab实现

土壤持水曲线是研究土壤水力学性质必不可少的 ,在已经建立的众多数学模型中 ,vanGenuchten模型是目前运用最广泛的模型 ,而运用该模型的关键是其 4个参数的求解。为此 ,本文对同一组东北褐土的土壤水吸力和对应的土壤含水量数据较详细地介绍了Matlab软件的非线性拟合和非线性回归函数的运用 ,得出了该土壤vanGenuchten模型的 4个参数值 ,分别建立了该土壤的vanGenuchten模型 ,并利用Mat lab强大的绘图功能对它们进行了直观比较。最后运用方差分析和残差分析对该模型的计算值与实测数据进行了分析 ,结果表明 :非线性拟合和非线性回归函数求参结果的显著水平均达到p<0 0 0 0 1,残差平方和均小于 0 0 0 0 5 ,其中非线性回归函数的求参结果较非线性拟合好。因此 ,运用Matlab软件的非线性拟合和非线性回归函数对土壤持水曲线的vanGenuchten模型进行求参是切实可行的 ,从而为土壤学工作者寻求出了一条运用数值计算方法的新途径。     

Change in the hydraulic properties of a Brazilian clay Ferralsol on clearing for pasture

Ferralsols under native vegetation have a weak to moderate macrostructure and a well-developed microstructure corresponding to subrounded microaggregates that are usually 80 to 300 μm in size. The aim of this study was to analyze how the hydraulic properties of a clay Ferralsol were affected by a change of structure when the native vegetation is cleared for pasture. We studied the macrostructure in the field and microstructure in scanning electron microscopy. The water retention properties were determined by using pressure cell equipment. We determined the saturated hydraulic conductivity, K_s, by applying a constant hydraulic head to saturated core samples, and the unsaturated hydraulic conductivity, K(Ψ), by applying the evaporation method to undisturbed core samples. Results showed a significant decrease in the water retained at −1 and −10 hPa from 0- to 40-cm-depth when the native vegetation is cleared for pasture. That decrease in the water retained was related to a smaller development of microaggregation and greater proportion of microaggregates in close packing. For smaller water potential, there was no difference of water retained at every depth between native vegetation and pasture. Pedotransfer functions established earlier for Brazilian Ferralsols and using clay content as single predictor gave pretty good results but the precision of the estimation decreased when the water potential increased. This decrease in the precision was related to the lack of predictor taking structure into account. K_s and K(Ψ) showed an upward trend with depth under native vegetation and pasture. Except at 0–7-cm depth between the Brachiaria clumps in the pasture where smaller K_s and K(Ψ) than at the other depth was recorded whatever land use, we did not record any significant difference of K_s and K(Ψ) at every depth between native vegetation and pasture. The upward trend shown by the hydraulic conductivity with depth was related to the increase in the development of microaggregation with depth.     

Estimating the water retention curve from soil properties: comparison of linear, nonlinear and concomitant variable methods

The unsaturated soil hydraulic functions involving the soil–water retention curve (SWRC) and the hydraulic conductivity provide useful integrated indices of soil quality. Existing and newly devised methods were used to formulate pedotransfer functions (PTFs) that predict the SWRC from readily available soil data. The PTFs were calibrated using a large soils database from Hungary. The database contains measured soil–water retention data, the dry bulk density, sand, silt and clay percentages, and the organic matter content of 305 soil layers from some 80 soil profiles. A three-parameter van Genuchten type function was fitted to the measured retention data to obtain SWRC parameters for each soil sample in the database. Using a quasi-random procedure, the database was divided into “evaluation” (EVAL) and “test” (TEST) parts containing 225 and 80 soil samples, respectively. Linear PTFs for the SWRC parameters were calculated for the EVAL database. The PTFs used for this purpose particle-size percentages, dry bulk density, organic matter content, and the sand/silt ratio, as well as simple transforms (such as logarithms and products) of these independent variables. Of the various independent variables, the eight most significant were used to calculate the different PTFs. A nonlinear (NL) predictive method was obtained by substituting the linear PTFs directly into the SWRC equation, and subsequently adjusting the PTF parameters to all retention data of the EVAL database. The estimation error (SSQ) and efficiency (EE) were used to compare the effectiveness of the linear and nonlinearly adjusted PTFs. We found that EE of the EVAL and the TEST databases increased by 4 and 7%, respectively, using the second nonlinear optimization approach. To further increase EE, one measured retention data point was used as an additional (concomitant) variable in the PTFs. Using the 20 kPa water retention data point in the linear PTFs improved the EE by about 25% for the TEST data set. Nonlinear adjustment of the concomitant variable PTF using the 20 kPa retention data point as concomitant variable produced the best PTF. This PTF produced EE values of 93 and 88% for the EVAL and TEST soil data sets, respectively.     

耕作措施对土壤物理性状的影响

耕作措施对0~20cm土壤容重影响较大,深层影响较小。稳定入渗率呈现出强烈的时间变异性,从整个生育期看,翻耕都表现出较高的稳定入渗率,随时间的变化,差异减小;不同负压下,翻耕稳定入渗率最高,免耕次之,铁茬最低,不同时期变化趋势一致;翻耕大孔隙较其它两耕作措施多,翻耕后,土壤大孔隙增多,比免耕和铁茬处理高65%左右,达到5%水平差异显著,免耕与铁茬差异不显著。随着时间的推移,不同耕作措施大孔隙比例有所下降,中小孔隙比例增加,大孔隙数目翻耕仍最高,从显著性看,与免耕和铁茬在5%水平显著。