基于协同克立格的地下水盐分插值方法研究.docxVIP

基于协同克立格的地下水盐分插值方法研究 水盐和水位对地下水的使用和土壤盐分的测定有一定影响（ka混淆等人，2005年），土壤养分和盐渍对农业的可持续生产有重大威胁（centin等人，2003年）。土壤盐分的增加不仅降低了作物产量，而且限制了作物种类的选择（siangh等人，2009），尤其是干旱半干旱地下水和沿海低地区的浅埋区（il’ichev等人，2008；guswa，2002）。对地下水和土壤盐分含量的准确估计和分布研究在区域和全球水土环境及相关研究中是一项重要的课题(Liu et al.2005)。对于水土而言其可溶盐的分布在一定程度上受到地形的影响(Salama et al.,1999),地形特征对地下水及土壤形成与利用具有重要作用(Lin et al.2007;Wu et al.,2008)。因此在地下水和土壤盐分含量的插值估算中区域地形因素的影响不容忽视。 协同克立格(Co Kriging)插值方法充分利用了待求变量和易测变量的相关性(Odeh et al.,1995),借助辅助变量提高主变量的预测精度,在地下水和土壤盐分空间插值中得到广泛的应用,是公认的有效插值方法(Triantafilis et al.,2001;Darwish et al.2007)。Kim等应用水化学、Kriging与Co Kriging方法对首尔以南的五松地区地下水硝酸盐污染进行了评价,阐明了硝酸盐的空间分布及与其它物理化学参数的关系,并解释了Co Kriging在复杂地下水化学数据插值中的优越性(Kim et al.,2009)。Eldeiry等综合应用地面验证和遥感数据对美国阿肯色流域土壤盐分含量的Kriging和Co Kriging插值方法进行了比较,认为Co Kriging方法对玉米地的土壤盐分估计精度最高,其次是小麦和苜蓿(Eldeiry et al.2010)。José等应用Co Kriging方法对西班牙东南部塞奎拉河域土地,以土壤水电导率为辅助变量进行了50 cm深度的盐分预测研究,较好地反应了作物根系层的盐渍化状况(Joséet al.,2010)。在国内相关的研究以西北干旱内陆和黄河三角洲为主,结果表明与普通克立格相比Co Kriging既能够减少样点数量又能提高估值精度(赵成义等,2003;姚荣江等2006)。 针对整个环渤海低平原水土盐分问题,进行大尺度区域布点,结合地形因素对水土盐分的空间分布格局进行研究的尚不多见。为此,本文以环渤海低平原为研究对象,以研究区高程为辅助变量应用协同克立格插值方法对水土盐分及水位埋深的空间分布进行研究,为提高水土盐分含量的插值估算精度,为地下水和土地资源管理利用和农业生产的合理布局提供科学依据。 1 渤海低平生态农业特点 环渤海低平原位于华北东部,包括河北、山东内陆低平原和滨海低平原两大类型,含天津、沧州等80多个县市。由黄河、海河、滦河冲积而成,地势低平,大部分海拔50 m以下,滨海区10 m左右。总面积9.6×104 km2;总耕地面积6000多万亩,是我国重要的粮棉、果蔬产区。 环渤海低平原属于欧亚大陆东岸暖温带半干旱季风气候区,冬春寒冷干燥,夏季炎热多雨。多年平均气温12.2℃,平均年降水量500~600 mm,全区年内降水分配不均,主要集中在6—9月,占全年降水量的60%~80%,年均蒸发量900~1400 mm,干燥度达1.5左右。 1.1 土样、土壤、内盐量测定方法 根据该区的土地利用类型和地下水状况,全区共测量水井130个,布设水样点128个、土样点127个(图1),取样点间距约20 km。取样范围为(N36°03′~N39°35′;E114°36′~E119°28′),取样面积8.97×104km2。取样时间为2010年4—5月,每个点采用GPS定位,土样深度为0~20 cm。 水位埋深用皮尺与测绳测定,水样矿化度采用重量法测定。土样全盐量测定过程如下,称取过2 mm筛的风干土试样50~100 g,按土水比1:5配制浸出液。采用蒸干法,添加15%双氧水溶液与2%碳酸钠溶液进行测定,详细操作方法参见《土壤农业化学分析方法》(鲁如坤等,1999)。 以1:5万地形图(等高距20 m)与SRTM高程数据为基础,用实测高程点进行高程验证(图2),最终确定高程点1120个,平均间距约6.7 km,覆盖研究区主要地貌单元。利用均方根误差描述高程精度,Zi为高程数据真值,zi为计算值,n为误差个数,计算得RMSE=6.52,说明所得到的高程空间插值结果满足数据分析要求。 1.2 小区域化变量和多随机变量协同共进模型 在地统计学中,半方差函数的一些重要参数如块金值、基台值和变程等可以用来表示区域化变量在一定尺度上的空间变异和相关程度,它是研究土壤特性空间变异的关键,也是克立格插