以研究区1990年、2000年、2010年、2015年的的四期Landsat TM/ETM遥感影像为数据源,运用ERDAS 9.3软件对四期遥感影像进行波段组合、几何校正、投影转换、图像镶嵌、区域裁剪等预处理的基础上,按照2007年8月颁布的《土地利用现状分类》国家标准,结合研究区土地利用特征及遥感影像特点,将研究区土地利用类型分为6个一级类型:耕地、林地、草地、水域、建设用地(城镇用地、工矿用地)、未利用地。利用监督分类和人工目视判读相结合的方法进行遥感影像的解译,生成四期土地利用数据(图1)。为保证解译质量和数据一致性,利用汉中市1:5万地形图以及Google Earth高分辨率影像对土地利用解译结果进行了选点验证,评价结果显示各土地利用类型的解译精度均超过85%。最后应用GIS技术对4期数据进行空间叠置,最终得出3个时段土地利用类型变化数据。

汉中市土地利用以林地和耕地为主,分别占总面积的73%和22%以上(表2)。1990—2015年汉中市土地利用变化总体趋势表现为草地、耕地、水域和未利用地减少,而林地和建设用地面积增加,其中耕地减少面积最多,达54 769.24 hm²,由1990年的24.19%下降到2015年的22.18%,草地面积减少了3 069.95 hm²,减幅为23.82%,水域面积减少了5 253.75 hm²,减幅为23.32%,未利用地减少97.61 hm²,减幅为11.29%; 面积增加最多的为林地,25 a间增加了60 308.86 hm²,增幅3.02%,建设用地面积增加3 757.30 hm²,增幅16.26%。1990—2000年,林地、建设用地和未利用地面积增加,草地、耕地、水域面积减少,建设用地增幅最大达8.38%,水域减幅最大达22.33%。2000—2010年耕地和未利用地分别减少了4.44%和6.40%,其他地类面积均增加,增幅最大的为建设用地7.65%。2010—2015年,除林地面积增加外,其他地类面积均减少,减幅最大的为草地和未利用地,分别减少了9.27%和7.29%。这与城镇化进程中,耕地被大量建设成住宅和公共设施有关,另外,退耕还林还草政策是导致汉中市耕地减少,林地增加的主要原因。

1990—2015年汉中市各县区土地利用类型面积的变化见表3,每个县区的林地面积都在增加,耕地面积均在减少,林地面积增加和耕地面积减少最多的是位于山区的镇巴县、宁强县、略阳县和勉县,位于盆地中心的汉台区增减幅度不大,这与退耕还林主要在坡度较大的山区实施有关。所有县区的草地面积均在减少,其中减少较多的为勉县和西乡县,略阳县、佛坪县和洋县面积变化不大,其他县区减少均在200 hm²左右,从时间序列来看,位于山区的镇巴县、宁强县、西乡县、勉县和留坝县草地面积减少主要发生在1990—2000年,而位于盆地中心的南镇县、汉台区和城固县则主要发生在2010—2015年,总体上所有县区在2000—2010年草地面积变化不大。城固县、勉县、洋县、西乡县、南镇县、宁强县水域面积减少较多,其他县区减少不明显。每个县区的建设用地面积都在增加,增加最快的是汉台区,其次是勉县、洋县、城固县、西乡县和南镇县,其他县略有增加。城固县和佛坪县未利用地面积有所增加,其他县区均在减少,面积变化主要集中在镇巴县、城固县、西乡县和勉县。

在3个研究时段内,2010—2015年汉中市土地利用类型变化最为强烈(0.426),而1990—2000年和2000—2010年分别仅为0.144,0.103(表4),说明2010年以来大多数土地类型受人类活动影响明显加强。1990—2015年各土地利用类型变化的速度由大到小的的顺序为:草地>水域>建设用地>未利用地>耕地>林地,增加幅度最大的为建设用地,减少幅度最大的为草地和水域。

表2 1990-2015年汉中市土地利用类型面积变化

运用CS计算模型,按50%的幅度上下调整各土地利用类型的价值指数,计算了汉中市不同年份的生态系统服务价值敏感性指数(表5)。结果显示,各地类的敏感性指数均小于1; ESV相对于VC是缺乏弹性的,说明所采用的生态服务价值系数适合研究区的情况,计算的生态服务价值是可信的。1990年以来,林地的敏感性指数呈增长的趋势,草地和耕地的的敏感性指数呈减小的趋势,说明林地生态服务价值变化对汉中市生态系统服务总价值具有放大作用,而草地和耕地则产生了缩小作用。未利用地的敏感性指数接近于0,表明未利用地生态系统服务价值变化对汉中市生态系统服务总价值的变化影响不大。水域敏感性指数呈现减小—增大—减小的变化,表明水域的生态系统服务价值变化对汉中市生态系统服务总价值的变化产生先缩小,后放大,再缩小的变化。敏感性指数最高的为林地,在0.884 2～0.898 7,说明林地对服务价值的贡献率最大,其次为耕地和水域,分别在0.083 4～0.092 2; 0.016 4～0.021 6,草地和未利用地均在0.002以下。

汉中市1990年、2000年、2010年和2015年的生态服务价值分别为368.15亿元、369.34亿元、372.36亿元和373.13亿元,总体上呈现逐年增加的趋势,25 a间增加了4.98亿元(表6)。林地的生态价值增加了9.82亿元,耕地和水域生态价值减少较多,分别减少了2.82亿元和1.86亿元; 草地和未利用土地生态价值略有减少,林地生态服务价值的增加足以弥补其他地类生态价值的减少。在各种土地利用类型中林地产生的生态价值最大,占总价值的88%以上,其次是耕地,占总价值的8.34%～9.22%,最少的为未利用地仅占总价值的0.01%左右。水域面积约占总土地面积的0.7%左右,但其产生的生态价值却占总生态价值的1.6%以上。可见,水域的生态效应与景观效应是十分明显的。

表3 1990-2015年汉中市各县区不同土地类型面积变化hm2

表4 1990-2015年汉中市土地利用类型的动态度

表5 汉中市各土地利用类型生态服务价值的敏感性指数

1990—2015年汉中市各县区土地的生态服务价值见表7,生态服务价值较大的为西乡县、洋县、宁强县、略阳县和南郑县,均在40亿元以上,其次是勉县、留坝县和城固县,均在30亿元以上,汉台区最小,仅有5.6亿元。镇巴县生态服务价值呈先增加后减少的变化,南镇县、汉台区、城固县和洋县则相反,其他县呈持续增加的趋势。总体上生态服务价值增加最快的为西部的宁强县、勉县和略阳县,增幅均在2.8%以上,其次为北部的留坝县和佛坪县,增幅在1%左右,东、南部普遍增幅较小,均在0.5%以下,处于盆地中部的汉台区和城固县则有减小趋势,特别是汉台区减少幅度较大,达2.03%,主要是耕地和水域大面积减少,而建设用地大幅增加所致。

表6 1990-2015年汉中市土地利用类型生态服务价值及其变化

表7 1990-2015年汉中市各县区生态服务价值及其变化106元

汉中市生态系统各项服务功能产生的价值由大到小依次为:土壤形成与保护、气体调节、水源涵养、生物多样性保护、气候调节、原材料生产、废物处理、娱乐文化、食物生产(表8),其中价值最大的4项服务功能始终占总价值的61%以上,这恰好反映了汉中市生态系统类型的分布格局和服务功能的构成情况,汉中市以林地生态系统为主,而林地的土壤形成与保护、气体调节和生物多样性保护价值远大于其他生态系统。近25 a来,废物处理和食物生产2项服务功能价值减少,其他7项服务功能价值均有不同程度的增加,其中增加较多的为土壤形成与保护、气体调节、原材料生产和生物多样性保护等4项功能,主要原因是产生较大土壤形成与保护、气体调节、原材料生产和生物多样性保护服务功能的林地面积不断扩大,而产生最大食物生产、废物处理服务价值的耕地和水域面积持续减少。汉中市各生态服务功能价值构成及其变化趋势与陕西省其他地区不同^[20],这与研究区所处的地理位置和生态条件等有关。

表8 1990-2015年汉中市生态系统各类服务价值构成及其变化106元

1990—2000年、2000—2010年、2010—2015年3个时段,林地对生态服务价值变化的贡献率最大(72.02%),其次为耕地(24.59%)和水域(25.19%)(表9)。1990—2015年林地、耕地和的生态贡献率分别为66.95%,12.67%和19.23%,三者之和超过98%。3个时期,林地、耕地、水域的生态贡献率变化较大,最大变化分别为13.64%,9.9%和22.02%,说明林地、耕地和水域产生的生态服务价值变化量对区域总服务价值变化影响较大,因而是主要贡献因子和敏感因子。

表9 1990-2015年汉中市不同土地利用类型生态服务价值变化的贡献率%