基于样条插值方法的干旱半干旱区植被根系层土壤水分特性的研究(新2)精选.doc

基于样条插值方法的干旱半干旱区 植被根系层土壤水分特性的研究 徐艳野金花 （黑龙江八一农垦大学文理学院大庆）内蒙古自治区鄂尔多斯 关键词：；；； 土壤水分是影响干旱半干旱地区植物生长和发育的重要因素，在这些地区开展植被根系层土壤水分特性研究对植被种植和生态恢复的重要性是不言而喻的。土壤水分特性描述了土壤含水量与土壤基质势之间的函数关系，是模拟水分和化学物质在非饱和土壤中运移的关键参数。目前测定土壤水分特征曲线的方法主要分为三种：第一种是直接测定法，如张力计法、离心机法、压力膜法、砂芯漏斗法、平衡水汽压法等。张力计法简单且便于操作，在室内和田间都得到了广泛应用；第二种是参数估计法，由于直接法在实际应用中存在诸多问题，特别是在区域尺度上进行实际问题研究时，这类方法多数是不可行的，甚至是不可能的，因此许多土壤物理学家尝试着用数学表达（经验公式）来描述水分特征曲线，通过估计表达式中的参数来确定土壤水分特征曲线，这种方法称为参数估计法。目前比较常用的经验公式有Brooks-Corey(1964)模型、Gardner(1970)模型、van Genuchten(1980)模Gardner-Russo(1988)模型等，近几年Ferdinand J Kastanek 等人又提供了一种新的用经典的三次样条插值函数来描绘土壤水分特性的方法,，这种方法较目前土壤水文学中应用的其它方法有许多优点；第三种是间接推求法。间接法主要是利用易于测定的土壤物理特性来推求土壤水分特征曲线，Arya．Paris(1999) 利用土壤颗粒组成来估算土壤水分特征曲线，近年夏达忠等也发表了《基于土壤物理特性的土壤水分特征曲线推求法》。这是近年来发展的利用土壤物理特性确定土壤水分特征曲线的新方法。本文主要采用张力计法对所选干旱半干旱试验区进行土壤水分含量的测定，利用所得数据我们建立了基于三次样条函数插值的土壤水分特征曲线模型，并利用所得结果分析出该地区土壤水分的变化规律，从而为干旱半干旱地区植被种植、退耕还林等方面提出了更有效的依据。 1 试验方案与研究方法 1.1 取样 我们选取具有代表性的内蒙古自治区科尔沁沙地为主要研究区，该地区属于温带半干旱气候类型选择科尔沁沙地典型人工植被作为研究对象分别以杨树林地 (株行距 1 . 5 m ×2 . 0 m年龄 12年 )、 樟子松林地 (株行距 1 . 5 m ×2 . 0 m年龄 25年 )15年生小叶锦鸡儿固定沙丘 (株行距 1 . 0 m ×1 . 0 m)天然草场为试验地。(体积为100 cm3)取土，每个试验区分3层，每层取3个土样。取样时同时记录毛根密集分布位置、个别粗根直径及分布位置、植物根系总体分布等特征信息。 1.2 试验数据及研究方法 以上四个试验区的土壤水分特征曲线的测定主要采用张力计法。用张力计(负压计)测定土壤负压h，土壤水分计算公式为：。 为了保证装土的初始含水率均匀、密度一致，土样经过粉碎和过筛处理，然后按一定的密度装填土样，装好的土样经过充分饱和后开始脱湿试验。受实验条件的限制，只测定了低吸力条件下的土壤负压。用张力计实测的低吸力条件下各个土样的负压h和用烘干法实测的各个土样的含水率数据见表1。 表1 试验区土壤水分特征曲线实测数据 杨树林地 樟子松地 小叶锦鸡儿固定沙丘h h h h 0.00 0.540 0.00 0.398 0.00 0.392 0.00 0.568 26.3 0.476 22.12 0.373 34 0.378 25.2 0.523 52.04 0.435 47.63 0.351 55.7 0.362 45.24 0.509 63.59 0.406 99.38 0.327 82.86 0.35 61.2 0.501 81.6 0.398 129.20 0.318 102.00 0.329 76.16 0.486 95.4 0.382 157.76 0.305 121.76 0.313 94.2 0.475 108.6 0.356 190.40 0.296 146.24 0.304 112.16 0.466 126.48 0.341 244.81 0.288 171.52 0.297 136.02 0.452 159.24 0.312 299.20 0.282 192.11 0.289 148.60 0.429 2 土壤水分特征曲线模型建立与求解 2.1 基于三次样条插值法模型的建立 根据表1中的数据，我们利用三次样条插值函数建立了试验区土壤水分特征曲线模型，其原理如下： 假设在闭区间 有一组结点：和一组对应的函数值：，若函数在子区间上是一个不超过三次的多项式，且在每一点上满足：，且一阶和二阶可导，如下式：