GIS与地统计学的土壤水分空间插值方法.docVIP

GIS与地统计学的土壤水分空间插值方法 土壤是由多种因素相互作用而形成的，其属性在空间分布上有一定的差异[1] . 实现精确农业养分管理和解决全球变化等环境问题的关键在于准确掌握土壤理化性状的空间变异规律[2] . 土壤干湿状况对监测土壤退化、土地干旱有着巨大作用，是监测土地状况的重要指标[3] . 土壤水分的时空分布研究对农业干旱监测和生态环境监测具有重大意义. 目前，揭示土壤属性空间分布特征的主要手段仍是通过野外采样和室内测定. 但其问题在于，无论采样密度多大，均无法得到所有区域的土壤属性值. 而空间插值技术可以将离散的数据点转化为连续的数据曲面从而实现研究区土壤属性的全覆盖. 作为土壤属性变化的时空定量监测方法，土壤属性空 间插值方法及其精度是数字土壤研究的重点. 姚雪玲等[4] 以黄土丘陵沟壑区洋圈沟小流域为例，介绍了基于ArcGIS 与地统计学的回归克里格模型插值方法，详细介绍了基于虚拟变量技术将分类变量引 入模型的过程，同时将回归克里格模型、回归线性模型和普通克里格模型的插值精度进行对比分析. 史文娇等[5] 介绍了土壤属性空间插值的常用方法与精度验证的方法和指标，并且总结了能够提高土壤属性插值精度的6 种途径. 陈思萱等[6] 针对广东省某地区土壤中采样样品点的重金属污染物砷(As)，选取统计学中具有代表性的3 种空间插值方法，对比分析不同空间插值方法下土壤重金属浓度估算精度和污染格局空间区划差异. 龙军等[1] 根据不同地貌类型特征在福建省各地级市选取了9 个典型县，利用2008 年采集的样点数据评价不同空间插值方法对土壤有机质含量插值精度的影响. 关键词：地统计学; GIS; 土壤水分; 空间插值 本文以山地平原过渡地区四川省绵竹市为研 究区，选择了几种常用的空间插值方法，包括反距离权重法(IDW)、径向基函数空间插值方法、普通克里格插值方法以及结合线性回归模型和普通克 里格的回归克里格插值方法，探究不同插值方法对 山地平原过渡地带土壤水分的影响. 1 研究区概况 绵竹市位于四川盆地西北部，位于山地平原过渡地区，地，1035410420E之间，面积1 245. 3 km 2 ，辖19 镇2 乡，总人口52万. 本文研究区为绵竹市部分区域(图1 ). 该区域属四川盆地中亚热带湿润气候区，气候温和，年均温15. 7℃ ，年降水量为1 053. 2 mm，雨量充沛，四季分明，大陆季风性气候显著. 地貌形态区域分异十分明显，西北部为山地，东南部为平原，地势西北高，东南低，由西北向东南倾斜. 森林占幅员面积49. 8% . 土地利用类型主要以林地、耕地、草地、建设用地、水体和未利用地6 类为主. 境内平原区土壤包括13 个土属，8 个亚类，土地肥沃，物产丰富. 2 数据准备及预处理 图1 研究区位置图 图上创建2 km 2 km 的网格，以各个格网中心作为预采样点，然后叠合土地利用数据，对预采样点 研究涉及数据主要包括：绵竹市行政区划矢量 图，来源于四川省测绘地理信息局;绵竹市2016 年土地利用类型数据为Landsat 8 OLI 解译数据，解译为耕地、林地、草地、水域、建设用地和未利用地共6 类[7] ，经实地样点验证，一级土地利用解译精度满足后续分析需要(总体精度达94. 3% );绵竹市数字地面模型(DEM)来自USGS(https：/ / www. usgs. gov / )，空间分辨率为30 m 30 m. 绵竹市土壤含水量样点数据通过野外采集获 取，考虑到样点空间分布对插值结果的影响，在样点选取时按以下思路进行. 首先在研究区行政区划进行适当调节以规避诸如水体、城市不透水面等区 域. 此外，为适应研究区西北部复杂地形条件下土壤水分空间分布的潜在异质性，采样时对样点进行 了适当加密，最终确定在空间上尽可能均匀分布的样点共184 个. 结合手持GPS，土壤采样工作在 2016 年4 ～ 5 月间进行，采样时选取有代表性的新鲜土壤，采样深度为0 ～ 20 cm. 取样时，先铲出一个耕层断面，再平行于断面下铲取土，混合土样太多时用四分法将多余土壤弃去[8] ，采样的同时记录样点的土地利用类型. 采样当天即采用烘干法进行土样含水量测定，具体过程：将鲜土捏碎放入已知准确质量的铝盒中，对其进行称量，精确至0. 001 g;然后，将其置于已预热至(105 2) ℃ 的烘箱中烘烤12 h，取出冷却至室温后立即称重，每个样点做3个平行测定. 若第一次测定结果中3 个平行样本之间的差距过大(超过0. 05% )，则进行二次烘干再次称重[9] . 利用ArcGIS 对土壤含水量样点数据进行数据子集提取，为满足不同插值