基于样条插值方法的干旱半干旱区植被根系层土壤水分特性的研究新2Word文档格式.docx
《基于样条插值方法的干旱半干旱区植被根系层土壤水分特性的研究新2Word文档格式.docx》由会员分享，可在线阅读，更多相关《基于样条插值方法的干旱半干旱区植被根系层土壤水分特性的研究新2Word文档格式.docx（13页珍藏版）》请在冰豆网上搜索。
以上四个试验区的土壤水分特征曲线的测定主要采用张力计法。
用张力计（负压计）
测定土壤负压h，用烘干法测定相应的土壤含水率
，土壤水分计算公式为：
。
为了保证装土的初始含水率均匀、密度一致，土样经过粉碎和过筛处理，然后按一定的密度装填土样，装好的土样经过充分饱和后开始脱湿试验。
受实验条件的限制，只测定了低吸力条件下的土壤负压。
用张力计实测的低吸力条件下各个土样的负压h和用烘干法实测的各个土样的含水率
数据见表1。
表1试验区土壤水分特征曲线实测数据
杨树林地
樟子松地
小叶锦鸡儿固定沙丘
天然草场
h
0.00
0.540
0.398
0.392
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
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中的数据，我们利用三次样条插值函数建立了试验区土壤水分特征曲线模型，其原理如下：
假设在闭区间
有一组结点
：
和一组对应的函数值
，若函数
在子区间
上是一个不超过三次的多项式，且在每一点上满足：
，且
一阶和二阶可导，如下式：
（1）
（2）
（3）
那么就称
为三次样条插值函数。
用三次样条插值意味着在左右两子区间上的三次多项式有不同的系数，但是在连接点处的数值、一阶和二阶导数分别相等，这样就使得相邻两区间可以平滑过渡。
（4）
此外，
（5）
（6）
并且对于自然三次样条，
这样解三次多项式所需的参数都可以通过解一系列的线性方程来求得。
2.2模型求解
利用Matlab优化工具箱对上述模型进行求解，得到的四个试验区土壤水分特征曲线如图1、图2、图3、图4所示。
将计算所得的土壤水分特征曲线与实测值进行比较，结果表明基于三次样条插值法的土壤水分特征曲线模型能够较好的反映不同植被类型根系层土壤水分含量与基质势之间的函数关系。
2.3结果分析
（1）研究表明，不同植被类型土壤含水量有着显著差异。
杨树林地、天然草地土壤含水量能够保持在较高水平。
这两种植被样地由于地下水位较高,土壤类型主要为草甸土,植被根系周围土壤能够得到地下水的补充,故土壤含水量能够维持较高水平。
樟子松林地、小叶锦鸡儿固定沙丘含水量较低,这两种植被样地均处在沙丘的中上部,土壤保水能力差,地下水位低,上层土壤得不到地下水的补充,完全依靠降雨的补充,土壤水分含量波动大。
（2）由上图可以看出，不同层次土壤水分特征曲线差异较大，而且种植植被不同，各层含水率的变化趋势也不一致。
在0—60
吸力范围内，樟子松地3个土壤剖面层次的土壤水分特征曲线变化较平缓，这表明该试验区土壤的重力水含量少。
而杨树林地土壤水分特征曲线变化很陡，表明土壤中大部分水分保持在这一范围内，重力水移动性大，但当降雨时，部分水量会迅速下渗到土壤深层，因此保持这部分水量使之变为土壤深层储水，对于改善土壤水状况以及提高水分利用效率有重要作用。
在60—200
吸力范围内，四种植被土壤水分特征曲线变化都比较平缓，表明土壤持水能力增强。
（3）通过对四种植被根系层土壤水分含量的分析发现，根系集中生长且发达的地方,植物水势高于土水势,这样土壤水分向着植物根系附近运动,使植物获得生存的条件，所以在该地区土壤水分条件明显优于其它地方。
3.1问题一的建模及求解
为分析干旱半干旱地区的土壤水分特征，我们通过实地考察，选取具有代表性的内蒙古自治区科尔沁沙地为主要研究区，该地区属于温带半干旱气候类型，选择科尔沁沙地典型人工植被作为研究对象，分别以杨树林地（株行距1.5m×
1.0m）、流动沙丘、丘间地、天然草场为试验地。
土壤水分测定采用：
在每个类型区随机设定4m×
4m样方3个，每个样方设置1个取样点，测定深度200cm（部分样地由于地下水位较高，取样深度随机变化），每10cm深度取样1层，每层3个重复。
针对干旱、半干旱区植被根系层土壤水分特性这一问题，由于土壤的水分是一个间接获得的信息量，所以通过建立模型Ⅰ：
烘干法模型，来测量土壤水分；
进而我们通过建立的模型Ⅱ：
基于三次样条插值法的模型，来从不同角度描述土壤的水分特征。
3.1.1模型Ⅰ的建立
根据实地考察，为降低实验的偶然性以及数据间的不确定性，我们在同一时间、使用同样地工具获得不同类型区域的样本土壤，为得到不同类型区间的土壤水分，我们建立了烘干法模型来测量土壤的水分：
（1）用已知重量的铝盒在天平上称取欲测土样15—20克。
（2）将盛土样的铝盒放入烘箱内，打开盖，在105～110℃温度条件下连续烘6小时，取出后，放入干燥器内冷却。
（3）将铝盒盖盖上，从干燥器中取出，称量。
（4）称后再将盖打开，放入105～110℃温度的烘箱中烘2小时，取出称重，如此连续烘至恒重（两次差数小于0.05克）。
（5）计算：
3.1.2模型Ⅰ的求解
通过烘干法的原理，我们对选取的对象（小叶锦鸡儿和榆树）进行水分测量，同时为了更具说服性，我们选取不同深度的根茎处的植物进行测量，得到如下水分数据（见表1和表2）：
小叶锦鸡儿
含水量（%）
根茎深度（cm）
1
0.03128–0.03375
33.57–26.12
2
0.03375–0.0428
26.12–21.19
3
0.0428–0.09628
21.19–14.12
4
0.09628–0.2304
14.12–11.09
5
0.2304–0.3275
11.09–10.33
6
0.3275–0.3316
10.33–0.3526
7
0.3316–0.3316
0.3526–0
表1小叶锦鸡儿含水量
榆树
含水量（%）
2.9893-2.5643
0-20
2.5643-3.59
30-40
1.73-3.59
40-50
1.73-2.6126
50-70
2.4367-3.0218
70-100
3.0218-2.12459
100-130
130-160
表2榆树含水量
由于所选取的植物对象出自同一地点，为了更形象且具体的体现土壤水分特征，我们又针对不同土壤中植物所含水分的多少，建立出vanGenuchten模型，从而计算出不同土壤所对应的含水量，根据公式：
式中：
为体积含水量
；
为滞留含水量
为饱和含水量
h为土壤吸力（cm）；
为土壤进气吸力（cm）；
，m，n，a为拟合参数，m=1-1/n。
其计算结果见下表3和折线图1所示：
土样
a
n
m
壤质沙
0.0657
0.4862
0.0124
1.4627
0.3163
0.984
西安土
0.0627
0.5047
0.0130
1.5003
0.3335
0.964
沙壤土
0.0359
0.4900
0.0181
1.6473
0.3929
0.976
粉壤土
0.0335
0.4988
0.0193
1.6636
0.3989
0.993
粉粘壤土
0.0317
0.3809
0.0426
1.8160
0.4493
0.992
表3不同土壤水分含量
图1不同土壤水分含量折线图
通过以上的数据及图像显示，可以看出：
（1）对小叶锦鸡儿而言，根茎浅处所含的水分较多，越往深处水分就相对较少；
而对榆树来说，水分含量分布较均匀，但以深处较多。
由此可知：
植物根茎的含水量与植物根茎深度有着很好的对应关系，出现土壤水分空间分布变化的原因，主要由于植物吸收水分，使土壤水分向植物根系附近运动，并且根系附近土壤疏松，毛管作用强，保持水性大，这样在植物根系发育范围的土壤含水率较大。
（2）根据不同土壤植物水分含量的数据，我们可以发现：
不同土壤的含水率是不同的，对其研究可以更好的针对不同土壤种植不同的植物，以达到人工植树造林的目的，例如：
人工种植小叶锦鸡儿，应充分考虑小叶锦鸡儿的根系状况，在灌丛空地可以种植浅根性草类，一方面它与小叶锦鸡儿竞争水分，另一方面又可以增加地表覆盖率以减少水分蒸发，这样使整个生态环境中的植被能够充分有效的利用土壤水分，促使生态环境进行正向演替。
3.1.3模型Ⅱ的建立与求解
根据模型Ⅰ所得的水分值，为详细的分析出干旱半干旱地区的土壤水分特征，我们建立了基于三次样条法的模型来进行求解，其原理如下：
根据三次样条插值法，我们可以更加精确方便的得出干旱半干旱地区土壤水分动态变化的环境因子的影响，呈现出非常复杂的动态变化，其动态变化呈现如下：
一、土壤水分的时间变化
降水、气温、光照强度和蒸发量等气候因子的季节性变化，会引起土壤水分发生相应变化。
植物的生长具有季节性，随着季节的变化植物对土壤水分的利用和植被覆盖地表的情况随之变化，对土壤水分含量也有一定影响。
据降水情况，可粗略地将一年分为雨季和旱季雨季前土壤含水量较低，随着降雨增加进入雨季，土壤含水量逐步增加并在一段时期内相对稳定在一个较高的水平，接着降雨量减少而进入旱季，土壤水分也随之减少并在一段时期内相对稳定在一个较低的水平。
总之，土壤水分季节性变化受该地区降雨的强烈影响，从总趋势上看，土壤水分季节性变化与当地气候的季节性变化，尤其是降雨的季节性变化基本是一致的。
2、降水
在干旱半干旱区降水对土壤水分动态有着强烈影响。
土壤的干湿季与气候的干湿季基本一致，土壤水分的变化趋势一般情况下与降水的变化趋势也基本一致，由于降水发生时间与水分入渗到相应土层有一定时间差，土壤含水量的峰值出现时间要滞后于降雨量峰值出现时间。
一次降雨的降雨量的大小对土壤水分的影响效果是不一样的。
在内蒙古自治区鄂尔多斯灌木林地，降雨量小于10mm会从地表迅速蒸发而不能补充土壤水分，故对于土壤水分补充来讲是无效降水。
土壤水分的丰水年和欠水年的形成与该年度降水量大于10mm的降水次数直接相关。
就降雨量相同的降水过程而言，不同的降水强度对土壤水分补充效果也不同。
当降雨量相同时，土壤含水量受中等强度降雨的影响最大，原因在于强度较大的降雨雨滴击溅易引起土壤空隙阻塞，造成下渗水量减少，而小强度历时长的降雨，对气温影响小，降雨期间蒸发量大，降雨对土壤水分的补充作用减弱。
降雨对土壤水分的影响是一个复杂的过程，既与降雨量有关，还与降雨强度和降雨历时长短等诸多因素有关。
3、气温和太阳辐射
蒸散损失是土壤水分损失的主要形式，它主要受气温、太阳辐射和风速的影响。
干旱和半干旱地区土壤平均含水量低，上层土壤温度在一天中的波动很大，土壤水分蒸散随土壤温度的波动变化很大。
干旱和半干旱地区光照充足，太阳辐射强烈，光照强度高和日照时数长，这样导致年蒸发量往往是年降雨量的数倍至数十倍，这些地区山地的坡向就成为影响太阳辐射强度和日照长度的关键性因子，也就是影响土壤水分的蒸散的关键性因子。
一般北坡（阴坡）的太阳辐射强度远低于南坡（阳坡），土壤含水量则明显高于南坡，这也是干旱和半干旱地区山地南坡造林难度大的主要原因。
通过三次样条差值法研究的干旱半干旱地区的土壤水分变化，除了与以上因素有关外，还与土壤本身的特性、坡向、坡度以及土地利用方式等都有一定的关系。
通过三次样条差值的描述，同时加入虚拟点的使用，使得曲线更加光滑，切合实际，将干旱半干旱地区土壤水分变化特征描述的恰当好处。
3.2问题二的求解
在问题二的基础上，根据最后的结论：
所选取的植被主要依靠降水补充，土壤水分波动较大，在半干旱风沙区的环境条件下，影响土壤水分变化的因素很多，通过实验和数据显示可知，土壤水分状况主要受到大气降水和土壤地下水供给的影响。
在地下水埋较浅的区内，虽然有植被存在，消耗了较多的水分，但土壤整体水分状况依然较好，能够满足植物生长所需的水分，植被可持续稳定生长存在。
在地下水埋较深的地区，特别是沙丘区，土壤水分的供给主要依靠降水的补充，由于植被生长的影响，植被区土壤水分低于流动沙丘区，土壤水分是限制植物生长和生存的主要因素之一，此类型植被区土壤水分受降水影响显著。
我们可以得出以下对干旱半干旱地区防风固沙、退耕还林的方法：
（1）在一定降水条件下，可通过植物种类和栽植密度调整方式，提出沙丘区最适宜的植被种类和密度管理方式。
在半干旱风沙区条件下，固定流动沙丘的植被应该采用耗水少的低矮灌木种类，同时应严格控制栽植密度，或通过间伐方式控制密度，才能保证土壤水分的供需平衡。
对半干旱风沙区受地下水影响的区域，如平坦的丘间地或甸子地等类型区域，在人工植被建立过程中，由于土壤有地下水的水分供给，可以根据生产需要确定合理的物种和栽植密度。
（2）应用草方格工程措施结合播种小叶锦鸡儿植物，建立沙地人工固沙植物群落是固定流动沙丘的一种重要途径。
它不但可以通过植物在不断生长过程中对沙丘表面加以固定，还可以通过植物地下庞大的根系和地上部的枯枝落叶对沙地土壤养分及化学性质具有一定的改良与提高作用。
（3）封育措施是恢复沙地植被的有效方法之一。
与未封育的对照相比，经过2~3年围封后，3种柳树灌丛的平均高度与平均地径有明显增加、生物量成倍增长；
灌丛内的草本植物盖度、平均高及产草量均有明显提高，总生物量（灌木层和草本层）也明显提高。
4模型评价
4.1模型的优点
（1）应用虚拟数据点得到的曲线与实测点之间拟合性很好；
（2）三次样条插值法对土壤水分特征曲线没有限制要求；
（3）烘干法直接测量，成本低，简单易行，准确性高。
4.2模型的缺点
（1）三次样条插值法需要很多的参数，不能用袖珍计算器来处理；
（2）烘干法费时费力。
参考文献
[1]BaveyeP,SpositoG.Macroscopicbalanceequationsinsoilsandaquifers:
Thecaseofspaceandtime-dependentinstrumentalresponse_WaterResource,1985（21）.
[2]苏敏,卢宗凡,李够霞.陕北丘陵沟壑区主要农作物水分利用与平衡[J].水土保持研究,1996,3
（2）:
36～451.
[3]阿拉木萨,裴铁璠,蒋德明.科尔沁沙地人工固沙林土壤水分与植被建设适宜度探讨
[J].水科学进展,2005,16（3）:
426-431.
[4]史小红,刘廷玺,李畅游.科尔沁沙地低缓沙丘主要植物根系层土壤水分特征分析[J].内蒙古农业大学学报,2003,24（3）:
23-29.
[5]高琼,董学军,梁宁.基于土壤水分平衡的沙地草地最优植被覆盖率的研究[J].生态
学报,1996,16
（1）:
33-39.
[6]牛海,李和平,赵萌莉.毛乌素沙地不同水分梯度根系垂直分布与土壤水分关系的研
究[J].干早区资源与环境,2008,22
（2）:
157-163.
升级会员 