土壤水分空间变异性研究进展.docxVIP

土壤水分空间变异性研究进展 水是土壤中蒸发、蒸发、渗透和径流过程的重要参与者。其变化和区域差异对区域水文、水土保护、生态和环境有显著影响。在大范围土壤水的研究中,通常需要布设监测站点或进行高密度田间采样。由于监测站点的维护和管理需要耗费大量的人力和物力,且土壤特性田间变异程度大,传统的田间采样方法比较费时、费力,常无法采集足够数量的样品。所以,如何利用较少的资料来获取比较多的信息,或者采用较易获取的资料来估算难以获取的水分信息成为研究的热门问题。 研究表明,地统计学是研究土壤属性或其他农田特征变异定量化的有效方法,已有的研究涉及到土壤水分特征及状态参数、物理性质、化学性质、重金属及其他元素。自20世纪80年代以来,Burgess及Webster等将区域化变量理论和普通克立格和协同克立格估值方法引入土壤性质空间变异性的研究中,推动了研究向前发展。随后一些研究者结合地统计学探讨了土壤水分的空间格局与土壤水文过程的联系以及土壤水分的时空变异与环境因素的关系。国内潘成忠等探讨了黄土半干旱丘陵区陡坡坡地土壤水分的空间变异特征,刘贤赵等研究了渭北旱塬不同土地利用条件下土壤水分的时空变异特征,赵学勇等揭示了科尔沁不同类型沙地土壤水分在降水后的空间变异性,张伟等分析了喀斯特地区典型峰丛洼地旱季表层土壤水分的空间变异性。至今,国内对海涂区土壤水分空间变异性的研究较少。本文以江苏北部滨海滩涂围垦区为例,研究了土壤含水量的空间变异性,探讨了田间土壤水分的普通克立格和协同克立格估值方法,旨在为该地区土壤水分分区管理提供科学依据,并为滨海滩涂地区土壤水盐优化调控提供一定的理论参考。 1 试验材料和方法 1.1 采样点布设和样品采集 研究区域为江苏省大丰市金海农场,地理位置介于32°59′30″~33°0′31″N、120°49′40″~120°51′4″E,总面积约180 hm2,东距黄海约4 km,西临大丰麋鹿国家级自然保护区。该区地处北亚热带季风气候区,具有明显的过渡性、海洋性和季风性,四季分明,年均降水量1058.4 mm,主要集中在6—8月份的雨季。由于当地缺乏蓄水工程,雨季水量丰富但不能调蓄,非雨季淡水资源紧缺,农田灌溉一般采用当地微咸水。在研究区域东部的种植制度主要为“水稻-油菜”的一年两熟制,西部则是“棉花-大麦”的一年两熟制。从图1显示的研究区域来看,受种植模式的影响,研究区东部和西部的排渠密度存在一定差异,表现为东部水稻种植区域的排渠密度(平均间距约50 m)明显大于棉花种植区域(平均间距约100 m)。 在研究区范围内进行网格布点采样,采样点平均间距大约200 m。布设的步骤是:首先从该研究区最南端的边界作一东西向的平行线,再以此线为起点线由西向东作平行线(至最北端边界),各平行线之间的间隔大约200 m;再由南向北作间隔200 m垂直于东西线的垂直线,最后以平行线与垂直线相交的“点”为采样点。规则网格数据的特点是样点间距离相等或成倍数,这种结构可以在研究区域上进行全部取样,也可以在网格上随机取样以减少工作量。在研究区内共挖取60个采样点剖面,其中水稻田25个,棉花地29个,盐蒿杂草地6个,实际田间采样点分布如图1所示。每个采样点剖面均对0~20、20~40、40~60和60~80 cm进行分层采样。各采样点的相对坐标采用差分GPS定位技术确定,样品采集于2007年10月下旬,正处于棉铃吐絮期的晚期(拔秆期)。土壤含水量的测定采用烘干法,将装有土壤样品的铝盒于105℃下烘至恒重,称重并计算出各相应采样土层的土壤质量含水量。 1.2 预测精度的表征 交叉验证用来评价普通克立格法和协同克立格法的预测精度。其基本思路是依次假设每一个实测数据点未测定,由所选定的变异模型,根据其他测定点数据用特定的克里格法估算这个点的值。均方根误差(RMSE)、平均标准误差(ASE)和预测值与实测值的相关系数r用来表征预测精度。RMSE和ASE越小且越接近、r越大,则预测的精度越高。本研究采用协同克立格预测的RMSE、ASE与普通克立格预测的RMSE、ASE减少的百分数(分别为RMSEa、ASEa)来表示预测精度的提高程度,并利用ra来表示协同克立格法对普通克立格法的相关系数的提高程度。 式中RMSEOK、ASEOK分别表示普通克立格法预测的均方根误差、平均标准误差;RMSECK、ASECK分别表示协同克立格法预测的均方根误差、平均标准误差;rCK和rOK分别表示协同克立格法和普通克立格法的预测值和实测值之间的相关系数。 2 结果与分析 2.1 含水量的同步控制 表1列出了各深度土壤含水量的统计特征值。从平均值来看,含水量变化范围在277.9~324.3 g/kg之间,虽然采样期正值旱季,但土壤含水量普遍较高,滨海滩涂区低地势及浅地下水埋深条件是造成这种现