-
基于RS与GIS技术的应用为分析沙化土地时空变化趋势及其分布格局提供了技术手段,也为利用转移矩阵方法分析沙化土地内部不同沙化类型的转移情况提供了技术保障。以理塘县沙化治理重点工程区为研究区,通过对影像数据进行解译,构建了沙化土地类型转移矩阵,分析其在5年内沙化土地结构的变化过程及其动态特征,并将主要治理模式与动态变化趋势进行空间叠加分析,宏观评估了各模式的治理成效,为理解沙化土地时空动态变化提供了新思路,对理塘县乃至川西北沙化土地治理工程宏观成效评估具有现实意义。研究区在2007—2012年间针对不同沙地类型主要实施了人工种草(低用量)、人工种草(高用量)、植灌种草、客土+人工种草、客土+植灌种草、沙障+植灌种草等6大治理模式(见表1),在部分冲蚀沟设置了挡沙墙、治理区进行了围栏封禁。2013—2017年,研究区主要开展了以前期沙化治理成果巩固为目标的防治工作[1]。
主要模式类型 模式内容 实施面积/hm2 人工种草(低用量) 混合牧草(4披碱草:4老芒麦:2燕麦)30~45 kg/hm2+有机肥(腐熟的牛羊粪)1.5—3 t/ hm2+围栏管护+补播 440.1774 人工种草(高用量) 混合牧草(4披碱草:4老芒麦:2燕麦)60~90 kg/hm2+有机肥(腐熟的牛羊粪)1.5—3 t/ hm2+围栏管护+补播 493.6198 植灌种草 高山柳(沙棘)2 500—5 000株/ hm2+混合牧草(4披碱草:4老芒麦:2燕麦)30—60 kg/hm2+有机肥(腐熟的牛羊粪)1.5~3 t/ hm2+围栏管护+补植补播 178.1265 客土+人工种草 客土500—1 000 m3/ hm2+混合牧草(4披碱草:4老芒麦:2燕麦)120—150 kg/hm2+有机肥(腐熟的牛羊粪)1.5~3 t/ hm2+围栏管护+补播 245.6159 客土+植灌种草 客土500~1 000 m3/ hm2+高山柳(沙棘)2 500—5 000株/ hm2+混合牧草(4披碱草:4老芒麦:2燕麦)90—120 kg/hm2+有机肥(腐熟的牛羊粪)1.5~3 t/ hm2+围栏管护+补植补播 116.5253 沙障+植灌种草 高山柳沙障(2 m*4 m)+高山柳5 000株/ hm2+混合牧草(4披碱草:4老芒麦:2燕麦)90~120 kg/hm2+有机肥(腐熟的牛羊粪)12 t/ hm2+围栏管护+补植补播 36.1576 Table 1. Main desertification land control measures in the study area
HTML
-
研究所用Landsat 8遥感影像来自“地理空间数据云”网站,数据源包括2013年Landsat 8遥感影像和2018年Landsat 8遥感影像各1景(见表2)。遥感影像选择原则:(1)以影像质量完好,确保研究区范围云层覆盖量尽可能少为原则(2018年的影像云量11.05%,但经过判断不影响研究区域)(2)选取研究区植被生长季节期内遥感影像,有利于沙漠化信息的识别与提取。基础数据还包括研究区1:5万地形图,以及影像解译辅助分类相关数据,包括:(1)2013年和2018年研究区的谷歌地图;(2)2013年涉及研究区的76个野外典型沙地类型调查样地资料;(3)2018年涉及研究区的72个野外典型沙地类型调查样地资料。同时,也收集了2007—2017年间研究区相关沙化治理模式布局图等。
传感器 影像时间 条带号 云量 Landsat-8 OLI_TIRS 2013.5.27 132/39 3.74% 2018.5.25 132/39 11.05% Table 2. Main information of remote sensing images
-
遥感影像数据的预处理及解译工作全部在Arcgis10.2软件中进行,包括投影转换、影像裁剪、影像几何校正(通过1:5万地形图进行)等,采用监督分类中的最大似然法进行分类,波段组合选择为654。再将波段组合后的数据(分辨率30 m)与第8波段(分辨率15 m的全色波段)进行影像融合,得到分辨率约15 m的假彩色影像初步数据。参照国内沙漠化遥感监测分类体系标准相关成果[2-5],根据地表流沙面积百分比、植被覆盖度以及地表景观特征将沙化类型确定分为轻度沙地、中度沙地、重度沙地、极重度沙地4大类(见表3),结合2013年和2018年的谷歌高清影像和随机抽取的50%野外典型沙地类型调查样地资料创建了4类沙化土地类型的训练样本,采用最大似然法进行解译,并对解译后的初步成果进行影像的平滑相关处理,并消除3个像元及以下的“孤岛”,得到解译成果栅格数据,并将栅格数据转化为矢量数据。再用剩余50%野外典型沙地类型调查样地资料对解译成果进行对比分析,解译准确度分别89.5%和88.9%,均在80%以上,满足研究要求。
沙化类型 植被盖度 流沙面积 景观特征 轻度 >50% <5% 土壤表层主要为土质,有斑点状流沙出露(<5%) 中度 30%~50% 5~25 地表分布有块状、片状流沙,但风沙流活动不明显,地表稳定或基本稳定。植被分布密度较低。 重度 10%~30% 25~50 原始地表形态基本破坏,沙地呈半流动状态,风沙流活动明显或流沙纹理明显可见,植被显小块状零星分布。 极重度 <10% >50% 流沙大面积分布,几乎无植被或仅有少量植被分布 Table 3. Characteristics of desertification land types in the study area
-
在Arcgis10.2对中两期沙化土地类型图进行叠加分析对获得的交集图层属性表赋值并计算其面积,再导出数据,并在EXCEL中通过数据透视表分析提取沙化土地类型转移矩阵。土地类型转移矩阵动态变化指标通常用动态度(D)反映[6-7]。单一沙化土地类型动态度表达的是研究区域一定时段范围内某土地利用类型的数量变化情况,其表达式为:
其中,LC代表在研究时间内某一种沙化土地类型动态度;Ua、Ub分别表示研究刚开始和研究结束时某一种沙化土地类型的总面积;T表示研究期的总时长,把T的单位设定为年时,LC就可以表示某一种沙化土地类型的年变化率。本研究通过年变化率来定量描述区域沙化土地类型变化的速度,比较沙化土地类型变化的差异。
1.1. 数据源
1.2. 遥感影像解译
1.3. 沙化土地类型转化分析
-
2013年研究区沙化土地总面积1 510.2225 hm2,其中轻度、中度、重度和极重度沙化土地面积分别为753.2775 hm2,229.9725 hm2,346.9725 hm2和180 hm2,分别占研究区沙化土地总面积的49.88%、15.23%、22.97%和11.92%,轻度沙化土地面积分布最广,其次是重度、中度、极重度(见图1)。极重度、重度沙地主要分布在丘陵中下部,呈相对集中连片、连块分布,主要原因是上世纪70年代在该区域大量成片取草皮,导致研究区草皮层下砂质土壤出露成沙;重度沙地和中度沙地多分布在极重度沙地附近,主要是过度放牧加之夏季降水集中,冲蚀加剧;轻度沙地分布广泛,主要分布在高原丘陵的中上部,人为干扰相对较小。
2018年研究区沙化土地总面积1510.2225 hm2,其中轻度、中度、重度和极重度沙化土地面积分别为952.965 hm2,116.3925 hm2,363.06 hm2和77.805 hm2,分别占研究区沙化土地总面积的63.10%、7.71%、24.04%和5.15%,轻度沙化土地面积分布最广,其次是重度、中度、极重度(见图1)。与2013年相比,极重度沙地和中度沙地分布总体上呈现“收缩”状态,而重度沙地和轻度沙地分布呈现“扩展”状态,主要原因是相关治理和管护模式使得部分极重度沙地逆转为重度沙地,部分中度沙地逆转为轻度沙地。
-
总体来看,2013—2018年的5年研究区沙化土地总体呈逆转的趋势,表现在:极重度面积减少102.195 hm2,年均减少11.36%;中度面积减少113.58 hm2,年均减少9.88%;轻度面积增加199.6875 hm2,年均增加5.30%;重度面积变化不大,年均增加0.93%。从各类沙化土地面积看,沙化程度相对严重的中度、重度、极重度沙地总面积由2013年的756.945 hm2下降到2018年的557.2575 hm2,下降幅度达26.38%,逆转趋势明显。
由表4、表5可知:(1)极重度面积由2013年的180 hm2减少到77.805 hm2,共减少102.195 hm2,年均减少达11.36%,主要逆转为重度和中度沙地。具体动态变化包括:①减少面积及来源结构:共减少面积110.5857 hm2,主要逆转为重度和中度;②增加面积及来源结构:共增加面积8.3925 hm2,主要来源于重度和中度沙地恶化。(2)重度面积由2013年的346.9725 hm2增加到363.06 hm2,共增加16.0875 hm2,年均增加0.93%,主要是一部分极重度逆转为了重度沙地,还有一部分是由中度、轻度沙地恶化所致。具体动态变化包括:①减少面积及来源结构:共减少面积149.265 hm2,主要来源于中度、轻度的沙地逆转;②增加面积及来源结构:共增加面积165.3525 hm2,主要来源于极重度逆转和中度沙地,少部分来源于中度、轻度恶化。(3)中度面积由2013年的229.9725 hm2减少到116.3925 hm2,共减少113.58 hm2,年均减少达9.88%,主要逆转为轻度沙地。具体动态变化包括:①减少面积及来源结构:共减少面积220.7025 hm2,主要来源于轻度沙地逆转以及恶化为重度沙地;②增加面积及来源结构:共增加面积107.1225 hm2,主要来源于重度和极重度沙地逆转。(4)轻度面积由2013年的753.2775 hm2增加到952.965 hm2,共增加199.6875 hm2,年均增加5.30%,主要是中重、重度沙地逆转为了轻度沙地。具体动态变化包括:①减少面积及来源结构:共减少面积20.205 hm2,主要是恶化为中度、重度、极重度沙地;②增加面积及来源结构:共增加面积219.8925 hm2,主要来源于中度、重度逆转为轻度沙地。
沙化土地类型 2013年/hm2 2018年/hm2 5年间沙地动态变化/hm2 沙地类型年变化率 轻度 753.2775 952.965 199.6875 5.30% 中度 229.9725 116.3925 –113.58 –9.88% 重度 346.9725 363.06 16.0875 0.93% 极重度 180 77.805 –102.195 –11.36% Table 4. Overall change of desertification land types in the study area
沙化土地类型 2018年(面积/hm2) 2013年面积比例 轻度 中度 重度 极重度 2013年合计 2013年(面积/hm2) 轻度 733.0725 5.895 13.8375 0.4725 753.2775 49.88% 中度 163.935 9.27 53.8875 2.88 229.9725 15.23% 重度 54.36 89.865 197.7075 5.04 346.9725 22.97% 极重度 1.5975 11.3625 97.6275 69.4125 180 11.92% 2018年合计 952.965 116.3925 363.06 77.805 1 510.2225 – 2018年面积比例 63.10% 7.71% 24.04% 5.15% – – Table 5. Transition matrix of desertification land types in the study area
-
将研究区内6大治理模式GIS图层与2013年的遥感解译的沙化土地类型矢量图层进行空间叠加分析,提取了各类模式实施面积及与沙地面积动态变化趋势的关系表(见表6)。
模式 沙化类型 面积/比例 总计 沙地变化趋势 逆转 无明显变化 恶化 人工种草(低用量) 轻度 实施面积/hm2 405.4572 390.6984 14.7588 面积比例/% 96.36 3.64 中度 实施面积/hm2 24.0668 2.6191 0.5768 20.8708 面积比例/% 10.88 2.40 86.72 重度 实施面积/hm2 9.7831 0.4989 9.1603 0.1239 面积比例/% 5.10 93.63 1.27 极重度 实施面积/hm2 0.8704 0.0585 0.8119 面积比例/% 6.72 93.28 人工种草(高用量) 轻度 实施面积/hm2 153.0564 147.7424 5.3140 面积比例/% 96.53 3.47 中度 实施面积/hm2 174.2728 137.4076 6.0129 30.8523 面积比例/% 78.85 3.45 17.70 重度 实施面积/hm2 157.2788 53.7885 101.0452 2.4451 面积比例/% 34.20 64.25 1.55 极重度 实施面积/hm2 9.0118 0.3343 8.6775 面积比例/% 3.71 96.29 植灌种草 轻度 实施面积/hm2 168.8379 168.7058 0.1322 面积比例/% 99.92 0.08 中度 实施面积/hm2 5.3870 3.9541 0.3843 1.0487 面积比例/% 73.40 7.13 19.47 重度 实施面积/hm2 3.3683 0.8571 2.3704 0.1407 面积比例/% 25.45 70.37 4.18 极重度 实施面积/hm2 0.5333 0.1274 0.4059 面积比例/% 23.88 76.12 客土+人工种草 轻度 实施面积/hm2 22.6106 22.6106 面积比例/% 100.00 中度 实施面积/hm2 11.5535 7.4577 1.4202 2.6756 面积比例/% 64.55 12.29 23.16 重度 实施面积/hm2 123.8280 67.1754 55.0106 1.6420 面积比例/% 54.25 44.43 1.33 极重度 实施面积/hm2 87.6238 84.7274 2.8964 面积比例/% 96.69 3.31 客土+植灌种草 轻度 实施面积/hm2 3.3154 3.3154 面积比例/% 100.00 中度 实施面积/hm2 14.6924 12.4964 0.8758 1.3202 面积比例/% 85.05 5.96 8.99 重度 实施面积/hm2 52.3028 21.8151 29.7994 0.6883 面积比例/% 41.71 56.97 1.32 极重度 实施面积/hm2 46.2147 21.9611 24.2536 面积比例/% 47.52 52.48 沙障+植灌种草 重度 实施面积/hm2 0.4116 0.0900 0.3216 面积比例/% 21.87 78.13 极重度 实施面积/hm2 35.7460 3.3788 32.3672 面积比例/% 9.45 90.55 Table 6. Relationship between main desertification control measures and dynamic change trend of sandy land types in the study area
由表6可知:(1)人工种草(低用量)模式可以明显遏制轻度沙地退化,其恶化面积比例仅3.64%;但对中度、重度、极重度治理效果较差,逆转面积比例分别仅为10.88%、5.10%、6.72%,且中度沙地的恶化比例达86.72%,重度、极重度沙地的现状保持面积比例分别高达93.63%、93.28%,说明该模式无法遏制中度沙地的进一步退化,对重度和极重度沙地治理基本没有效果。(2)人工种草(高用量)模式对轻度、中度沙地治理效果较好,可以明显遏制轻度沙地退化,其恶化面积比例仅3.47%;对中度沙地逆转比例为78.85%;对重度沙地有一定效果,逆转比例为34.20%。(3)植灌种草模式对轻度、中度沙地治理效果较好,对重度、极重度沙地也有一定效果,其中:轻度沙地恶化面积比例仅0.08%,对中度沙地逆转比例达73.40%;对重度、极重度沙地有一定效果,逆转比例分别为25.45%、23.88%。(4)客土+人工种草模式,主要针对中度、重度、极重度沙地,治理效果较明显,对中度沙地的逆转率为64.55%;对重度沙地的逆转率为54.25%;对极重度沙地逆转率为96.69%,极重度沙地的逆转率反而高于中度和重度的主要原因,是该模式所在的极重度沙地治理区周边大多配套有挡沙墙工程措施,强化了治沙成效。(5)客土+植灌种草模式主要针对中度、重度、极重度沙地,治理效果较明显,对中度沙地的逆转率为85.05%;对重度沙地的逆转率为41.71%;对极重度沙地逆转率为47.52%,极重度沙地的逆转率反而高于重度的主要原因,是该模式所在的极重度沙地治理区部分配套有挡沙墙工程措施固沙,强化了治沙成效。(6)沙障+植灌种草模式主要针对重度、极重度沙地,而且是位于沙化冲蚀沟的沟口,是沙化最严重的地区,该模式能够有效固沙,且能有效遏制重度沙地进一步退化,也能起到一定逆转效果(逆转率为21.87%)。