-
生态系统服务价值是生态系统服务功能的货币化形式。区域生态系统服务价值越高在一定程度上表明区域生态服务功能越好,人们的生态意识越强,越有利于区域生态安全。20世纪90年代以来,随着人类对生态系统功能认识的加强,生态系统服务价值研究受到广泛关注,最具代表性的是Costanza等[1]提出的“生态系统服务价值与自然资本”理论,然而其研究中的“全球静态部分平衡模型”是基于全球生态系统服务价值均值化的处理,在地域性生态系统服务价值评估中会产生较大误差。谢高地等[2]在Costanza等[1]研究基础上结合中国的自然特征制定了青藏高原单位面积生态系统服务价值当量因子, 成为国内生态系统服务价值研究的理论基础和主要方法参考。近30年来,国内外学者从景观尺度、单一生态系统、物种和生物多样性等方面开展了生态系统服务价值多角度研究[3-8], 不仅夯实了生态系统功能及价值评估的理论基础,同时也较有力的推动了该领域研究方法的创新[9-15]。森林是陆地生态系统的重要组成部分,在维护区域生态平衡中具有关键作用。国外研究侧重于林地的单一生态服务功能的测评及森林生态系统服务与局地生态环境的关系分析等方面[16-20],国内研究则侧重于从土地利用角度评估区域生态系统服务价值,王丽等[21]从涵养水源价值、土壤保持价值、固定CO2价值、净化空气价值、社会与文化价值5个方面评估了青城生态系统服务价值,郑江坤等[22]利用马尔科夫链模型评估了潮白河流域森林生态服务价值的动态变化,马国军等[23]和汪有奎等[24]对干旱和半干旱区的森林生态系统服务价值进行了评价。
干热河谷是中国西南地区一种局地特殊的地理景观和气候类型,气候炎热、干燥,水分蒸发比极为失衡,土壤侵蚀和水土流失严重,植被类型为“河谷型萨瓦纳植被”,特殊的自然生态环境使干热河谷成为西南地区典型的生态脆弱区。林地是干热河谷山地生态系统的重要组成部分,具有水源涵养和气候调节等多种功能,是维护干热河谷区生态系统功能的主要地类,林地的数量及时空分布对干热河谷区域生态安全的影响显著。现有的干热河谷林地研究侧重于林木育苗[25-26]、种群结构[27]和植被恢复[28-29]等方面,这些研究虽从单一角度部分探讨了林地的生态功能,但是无法有效系统量化林地生态系统服务价值。尽管有学者曾研究了林地生态系统服务静态价值[30],但没能深入评估林地生态系统动态价值的变化,不能揭示社会、经济和环境因素对干热河谷林地生态系统服务价值的影响。鉴于此,本文以具有干热河谷典型环境特征的元谋县为案例,研究干热河谷林地的时空格局及生态系统服务静态价值和动态价值的变化特征,辨明影响林地生态系统服务价值变化的主要因素,以期为干热河谷林地可持续发展及区域生态环境保护提供科学依据。
HTML
-
本研究的数据源为:(1)遥感数据,选取2008年Landsat7 ETM及2016年Landsat8 OLI卫星遥感影像,过境日期分别为2008年1月28日和2016年1月4日,分辨率为30 m。两景遥感影像时段均处于研究区旱季。数据由中国科学院地理空间数据云平台提供。(2)地理信息系统数据,选取研究区1∶100 000的国家地理基础信息系统数据,包括水系、道路和居民点等基本信息。数据由国家基础地理信息中心提供。(3)DEM数据,选取ASTER GDEM V2产品,分辨率为30 m。研究区共涉及到4个图幅(ASTGTM2 _N26E 101、ASTGTM2_N26E 102、ASTGTM2_N25E 101、ASTGTM2_N25E102)。数据由中国科学院地理空间数据云平台提供。
在ArcGIS中利用Mosaic工具对4个图幅的DEM进行拼接,拼接图幅包含了研究区区域范围。利用水系和高程点为控制点完成拼接的DEM图幅与1∶100 000基础地理信息系统数据的空间配准,以研究区范围的矢量数据裁剪拼接的DEM图幅获得研究区范围DEM,利用ArcGIS栅格计算功能提取研究区4个垂直自然带的空间范围,利用ENVI软件完成两景遥感影像的去云、去条带、标准假彩色合成及几何精校正等图像处理工作。完成2个时段遥感影像与1∶100 000地形图的空间配准,并以研究区矢量边界为掩膜对裁切遥感影像获得了研究区范围内的遥感影像。
-
根据《国家森林资源连续清查技术规定》地类划分标准和研究区遥感影像特征,结合前人研究成果[30-33]和实际研究需要,将研究区分为林地和非林地2个大类别,非林地包括草地、耕地、建设用地、未利用地和水域5种一级土地利用类型。
-
综合运用非监督分类、监督分类及目视解译等方法进行遥感图像解译。以土地利用类型历史图件与GPS野外采点调研验证了土地利用类型的解译精度。2期遥感图像的解译精度均在88%以上,符合解译精度要求(表 1)。在ENVI软件中以像元二分法提取研究区2008年和2016年的植被覆盖度。利用ArcGIS软件提取主要土地利用类型的基本统计数据,计算林地土地利用指数和土地状态指数[34],利用Dissolve和Intersect工具对2008年和2016年2个时段土地利用类型进行叠置分析,得到研究区土地利用转移矩阵。
年份Year 林地
Forest land草地
Grassland耕地
Arable land建设用地
Construction land未利用地
Unused land水域
Water land综合精度
Synthetic accuracy2008 86.1 86.4 90.6 86.3 88.4 93.2 88.5 2016 88.9 86.9 91.4 85.7 87.2 94.5 89.1 均值Mean 88.8 Table 1. Interpretation accuracy of land use types
% -
谢高地等[2]的研究以一级土地利用类型为研究基础,将不同的林地类型归为一种,不能有效揭示植被种类之间的差异性导致的生态系统服务价值的差别。李晓赛等[35]在方精云等[36]对植被生物量差异性研究的基础上,基于材积源生物量法提出了有林地、灌木林和其它林地的生态系统服务当量因子系数的调整方案,取3种林地的功能性调整系数平均值为林地生态系统服务静态价值当量因子调整系数,即1.032(表 2)。研究区的林地主要类型为有林地、灌木林和其它林地。为了更好地反映干热河谷不同林地的生态系统服务价值的差异性,本研究采用李晓赛等[35]的研究方法设定研究区林地生态系统服务当量因子的调整系数为1.032。
-
单位面积食物生产功能价值当量因子是指1 hm2农田每年自然粮食的平均产量的产值, 林地生态系统服务价值当量因子则以林地生态系统产生该生态服务相对于农田食物生产服务贡献的大小:
式中:Ea为单位面积农田生态系统提供食物生产服务功能的经济价值(元·hm-2);Ti是研究区每年粮食作物的经济产值(元·a-1);Mi是研究区粮食作物的总面积。1/7是指在没有人工干预的自然生态系统提供的经济价值与现有单位面积农田提供的食物生产服务经济价值的比值。以2008年和2016年研究区农田粮食作物(水稻、玉米和杂粮)的价格为生态系统服务价值系数计算的基本依据,依据统计资料计算得到2个年份的单价分别为979.36、1 142.55元·hm-2,取二者的平均值1 060.96元·hm-2作为单位面积农田食物生产功能价值。结合调整后的林地生态系统服务当量因子,通过公式(1)计算研究区林地及其它地类单位面积生态系统服务价值。
-
经济学研究表明, 对物品的支付水平与区域内某时段人们的支付能力和基于物品的需求程度的支付意愿有关[35],同时由于研究区所提供生态资源并非一成不变,而是受生态环境变化的影响,因此,本研究将支付能力、支付意愿[37]和环境能力等指标纳入生态系统服务动态价值评估体系,对动态价值进行系数调整,获得研究区生态系统服务动态价值系数。
-
支付能力是个人经济能力的直接表现。个人支付能力越强,越有可能承受对生态系统服务功能的支付。本研究将研究区人均国内生产总值与当年中国人均国内生产总值的比值作为支付能力指数,其计算公式如下:
式中:GDPi表示第i年元谋县人均国内生产总值;GDPimean表示第i年中国人均国内生产总值。利用统计数据通过公式(2)计算研究区2008年和2016年生态系统服务价值支付能力指数。
-
人们对生态系统价值的认识和需求通常与社会发展阶段相关[35]。社会发展阶段越高,人们的生态环境意识及对生态服务功能的需求越多(直至饱和状态),支付意愿也越强。因此,本研究对支付意愿指数(Ni)进行调整以反映社会发展阶段对生态系统服务价值的动态影响。支付意愿通常用Logistics生长曲线模型来刻画,其表达式为:
式中:Ni为代表生长特性的参数,在此表示第i年社会发展阶段系数;t在此表示社会经济发展阶段;a、b为常数取值为1,e为自然对数。当t值很小时,即社会发展水平很低,Ni值趋于0;当t值很大时,社会发展水平很高,Ni值趋于饱和值[35]。社会发展阶段的计算通常与反映生活支出结构的恩格尔系数(En)的倒数对应起来,表达式如下:
利用统计数据由公式(4)计算得到研究区2008年和2016年的t值分别为-0.46、-0.32,以公式(3) 计算研究区2008年和2016年生态系统服务价值支付意愿指数。
-
植被的生物量不仅与植被的类型有关,同时受当年水分和热量等环境条件的影响。植被覆盖度是反映植被生长状况的主要指标,能较客观的反映植被的生产力和生物量的变化。本研究以植被覆盖度作为基本参数构建环境能力指数(Si)反映研究区生态系统服务价值因环境条件而产生的动态变化,其公式为:
式中:Sd表示该年度研究区植被覆盖度,Smean表示研究区多年植被覆盖度均值。根据遥感图像运算结果,以公式(5)计算研究区2008年和2016年的生态系统服务价值环境能力指数。
基于上述研究成果,同时参考相关[37-40]研究,提出生态系统服务动态价值调整系数(Mij),在此基础上构建研究区生态系统服务动态价值评估模型(Qij),其具体模型如下:
式中:i为第i个年份,j为第j种土地利用类型;Pij表示居民在第i年第j种土地利用类型生态系统服务价值的支付能力指数;Nij表示居民在第i年第j种土地利用类型生态系统服务价值的支付意愿指数;Sij表示第i年第j种土地利用类型生态系统服务价值的环境能力指数;Qij为研究区第i年第j种土地利用类型生态系统服务动态价值;Cij为第i年第j种土地利用类型的面积;Eij是第i年第j种土地利用类型单位面积生态系统服务静态价值(元·hm-2);通过公式(6)计算2008年和2016年研究区生态系统服务价值动态调整系数值,林地和其它地类被赋予相同的生态系统服务动态价值调整系数值。利用公式(7)分别计算研究区2008和2016年及林地和其它地类生态系统服务动态价值及区域生态系统服务总动态价值。
2.1. 数据源及数据预处理
2.2. 土地利用分类体系
2.3. 土地利用数据的提取
2.4. 林地生态系统服务静态价值系数调整
2.4.1. 林地生态系统服务当量因子系数调整
2.4.2. 单位面积食物生产功能价值当量因子
2.5. 林地生态系统动态服务价值系数调整
2.5.1. 支付能力指数(Pi)
2.5.2. 支付意愿指数(Ni)
2.5.3. 环境能力指数(Si)
-
2008—2016年期间,研究区林地空间分布的基本格局未发生明显变化,即林地主要分布在南部和东部中高山,在西部、西南部和北部、东部边缘中高山地带有少许分布,在河谷坝区也有零星分布(图 2)。林地是东部和南部中高山的基质景观,与部分耕地和草地混交在一起,在其它区域有少许草地和林地混交。研究区海拔1 600 m以上属于中高山地带,从亚热带山区逐渐过渡到暖温带山区。海拔的升高明显改善了水热条件,为有林地、灌木林及其它林地的发育提供了条件。林木类型从云南松(Pinus yannanensis Franch)等滇中高原树种逐渐过渡到亚热带半湿润常绿阔叶林及针叶林。海拔1 6002 000 m的林木类型为云南松和旱冬瓜(Alnus nepalensis D.Don)等乔灌木,海拔2 500 m以上地区的林木类型以云南松为主,其次为滇栲(Castanopsis diversifolia Franch)和马缨花(Rhododendron delavayi Franch)等针阔叶混交中幼林[21]。从研究区林地的空间分布变化看,2016年金沙江沿岸的东部中高山及南部中高山部分耕地转化为林地,表明退耕还林等生态保护工程对林地恢复有积极的促进作用;然而,与2008年相比,2016年中这2个区域的林地斑块增多,破碎化程度更明显,主要是因为人为干扰的加强(特别是人工林的建设)使林地破碎化程度加深。
-
通过本研究的数据源获得的研究区范围面积为2 032.84 km2。2008—2016年,研究区林地的面积从2008年的562.08 km2增长至2016年的626.32 km2,增长了64.24 km2,变化幅度达到11.42%,年变化率达1.42%。林地占研究区总面积比例由2008年的27.65%增至30.81%,增长了3.16%,增幅较为明显。研究区林地的土地利用综合指数由2008年的55.30增至60.20,增长了4.90,表明林地土地利用指数有所增长且增幅较大。
从表 3中可以看出:林地与非林地之间相互转化较为频繁。林地向非林地的转移面积为90.53 km2,由非林地转入为林地的面积为154.77 km2。林地主要与非林地中的草地、耕地和未利用地和建设用地相互转化。林地转为草地的面积最多,达到73.61 km2;其次是林地向耕地和建设用地的转移,分别达到14.02、2.9 km2。非林地向林地的转移中,草地的转移面积最大,为123.11 km2;其次是耕地和未利用地,分别向林地转移27.22、4.44 km2。林地虽与草地相互转化量最大,但林地以转入为主,有49.5 km2的净增长,主要原因在于研究区林地以灌木林为主,有林地少。灌木林和灌草丛的生境相似,易受水分条件的波动及人为干扰的影响而相互转换,发展人工林使更多的草地转为林地;林地与耕地的相互转化中,林地以转出为主,净减少13.2 km2,主要受耕地空间分布变化的影响。退耕还林和天然林工程及水土保持工作的开展使山区坡度≥25°的耕地转为林地,而坝周低山区和河谷坝区地势较平,坡度小,土地开发难度相对较小,农业经济发展对耕地的强烈需求使更多的林地开发成耕地,因而林地在与耕地的相互转化中表现为净减少;林地与未利用地的相互转化中,林地以转入为主,净增长4.44 km2,主要是由于经济利益驱使的土地开发和生态恢复工程的实施使未利用地向林地转化,因而林地获得一定净增长;林地与建设用地的相互转化中,林地以转出为主,净减少2.9 km2,这主要是由于交通运输的发展使部分林地转化为道路,居民用地增长使部分林地转为城乡居民用地,因而林地表现为一定的净减少。林地的土地利用状态指数为0.26,指数为正且数值较小,表明林地以转入为主,且转入转出较为频繁,这与林地转移状况的研究结果一致。
土地利用类型
Land use type林地
Forest land草地
Grassland耕地
Arable land建设用地
Construction land未利用地
Unused land水域
Water land2008年面积
Area in 2008林地Forest land 471.55 73.61 14.02 2.90 0.00 0.00 562.08 草地Grassland 123.11 812.12 70.39 7.34 3.66 0.00 1 016.62 耕地Arable land 27.22 25.66 192.86 10.17 2.87 0.00 258.78 建设用地Construction land 0.00 0.00 0.00 82.34 2.23 0.00 84.57 未利用地Unused land 4.44 6.26 18.91 1.53 42.27 4.04 77.45 水域Water land 0.00 0.38 0.00 0.00 5.08 27.88 33.34 2016年面积Area in 2016 626.32 918.03 296.18 104.28 56.11 31.92 2 032.84 Table 3. Land use transfer matrix
km2 -
通过对林地生态系统服务当量因子系数的调节,结合谢高地等[2]的生态系统服务当量因子基础值确定研究区林地生态系统服务当量因子及林地单位面积生态系统服务价值(表 4)。
生态系统服务功能
Ecosystem service function气体调节
Gasadjust气候调节
Climateadjust水源涵养
Water sourceconservation土壤形成与保护
Soil formationWith protection废物处理
wastedeal生物多样性保护
Biological diversityprotection食物生产
Foodproduce原材料
Raw materials娱乐文化
Entertainment合计
Total当量因子
Equivalent factor3.61 2.79 3.30 4.02 1.35 3.36 0.10 2.68 1.32 价值系数
Value factor/×104(元·hm-2)0.383 0.296 0.350 0.427 0.143 0.356 0.011 0.284 0.140 2.390 Table 4. Equivalent factor of forest ecosystem service and forest ecosystem service value per unit area
综合研究区土地利用变化(表 3)及林地单位面积生态系统服务价值研究(表 4),分析2008年和2016年研究区生态系统服务静态价值的变化状况(表 5)。2008年和2016年研究区生态系统服务总静态价值分别为25.98亿元和26.54亿元,2016年较2008年增长了0.56亿元,增长了2.16%,平均增幅为0.27%,总价值增长幅度不大。2016年林地生态系统服务静态价值为14.97亿元,较2008年的13.44亿元增长了1.53亿元。2016年林地在研究区生态系统服务静态服务总价值比例为56.41%,较2008年的51.73%提高了4.68%。2个年份林地在研究区生态系统服务静态价值的比重都在50%以上,远高于草地及其它地类(表 5)。林地在区域生态系统服务静态价值的比重大且增幅明显,表明林地是干热河谷生态安全服务价值的关键类型,且其作用力不断增强,主要原因在于林地自身的单位面积生态系统服务静态价值高,林地面积的增长使其生态系统服务静态价值增长更为明显,在区域生态系统服务价值的比重也提升较快。自20世纪50年代以来,人类大幅砍伐森林用于薪材及其它生活和生产活动,导致20世纪80年代研究区的森林覆盖率不足10%,森林的生态功能遭到毁灭性的破坏,极大的危害了长江流域的水土安全,引起了国家和各级政府的高度关注。作为长江上游水土保持重点区域之一,元谋强化了天然林工程、退耕还林还草工程及其它植被恢复措施的实施,林地比例的不断增加促进了元谋干热河谷生态环境的恢复。
土地利用类型
Land use type2008年 2016年 变化值
Change value/×109元变化率
Change rate/%年变化率
Annual change rate/%价值
Value/×109元比例
Percentage/%价值
Value/×109元比例
Percentage/%林地Forest land 13.44 51.73 14.97 56.41 1.53 11.38 1.42 草地Grassland 9.83 37.84 8.87 33.41 -0.96 -9.77 -1.22 耕地Arable land 1.90 7.31 2.17 8.18 0.27 14.21 1.78 建设用地Construction land -0.86 -3.31 -1.06 -3.99 -0.20 -23.26 2.91 未利用地Unused land 1.63 6.28 1.56 5.88 -0.07 -4.29 -0.54 水域Water land 0.04 0.15 0.03 0.11 -0.01 -25.00 -3.13 合计Total 25.98 100 26.54 100 0.56 2.16 0.27 Table 5. Static value of ecosystem service of all land use types in different years
林地9个单项生态系统服务功能按价值量的大小排序依次为: 土壤形成与保护>气体调节>生物多样性保护>水源涵养>气候调节>原材料>废物处理>娱乐文化>食物生产(表 6)。林地在土壤形成与保护、气体调节和生物多样性保护的作用最大,食物生产价值最弱。由于研究区林地以灌木林为主,灌木林是研究区的景观基质,较好的适应了当地的干热气候。林木根系在土壤中的延伸为土壤的形成和保护提供了丰富的养分和水分来源,是固土和水土保持的关键因素,因而林地的土壤形成与保护单向价值高。林木的光合作用在固定CO2和释放O2方面具有重要作用,故气体调节功能价值高。林地植被覆盖度高,生态环境质量好,为研究区动植物生存提供了良好的生境,因而生物多样性保护价值高。灌木林与有林地相比生物量偏低,食物产出少,以灌木林为主的林地结构导致林地食物生产价值弱。各单项生态系统服务价值虽有一定差异,但林地总体具有较高的生态系统服务价值。
年份
year生态系统服务功能
Ecosystem service function合计
Total气体调节
Gasadjust气候调节
Climateadjust水源涵养
Water sourceconservation土壤形成与保护
Soil formationWith protection废物处理
wastedeal生物多样性保护
Biological diversityprotection食物生产
Foodproduce原材料
Raw materials娱乐文化
Entertainment2008年 2.15 1.66 1.97 2.40 0.81 2.00 0.06 1.60 0.79 13.44 2016年 2.40 1.85 2.19 2.67 0.90 2.23 0.07 1.78 0.88 14.97 变化值Change value 0.25 0.19 0.22 0.27 0.09 0.23 0.01 0.18 0.09 1.53 Table 6. Individual static value of forest ecosystem service in different years
×109元 -
2008年和2016年生态系统服务动态价值调整系数值分别为0.354和0.603(表 7)。2008年和2016年研究区林地生态系统服务动态价值分别为4.758亿元和9.027亿元,仅为同期林地生态系统服务静态价值的35.40%和60.30%,这主要是因为研究区经济发展水平较落后,人均GDP与全国平均水平的差异较大,恩格尔系数高,因此人们对林地生态系统服务价值的支付能力和支付意愿较低,导致林地生态系统服务动态价值总体偏低。2008年和2016年研究区生态系统服务总动态价值分别为9.197亿元和16.004亿元,林地生态系统服务动态价值占研究区生态系统服务总动态价值的比例分别为51.73%和56.41%,表明林地是区域生态系统服务价值的关键类型。
年份
Year支付能力
Payment ability支付意愿
Payment willingness环境能力
Environment capacity调整系数
Adjustment coefficient2008 0.456 0.774 1.004 0.354 2016 0.712 0.840 1.008 0.603 Table 7. Adjustment coefficient of dynamic value of forest ecosystem service
从林地生态系统服务动态价值变化状况看,8年间其值增长了4.269亿元, 增幅达89.72%,远超过同期林地生态系统静态服务价值增长率(11.38%)。这是由于较2008年相比,2016年研究区人均GDP增长非常迅速, 与全国平均GDP水平的差异缩小了25.6%,支付能力提高明显。社会和经济的迅速发展较大的提升了社会发展阶段,恩格尔系数降低,逐渐接近0.40的拐点,支付意愿得到提升,因而生态系统服务动态价值显著提升。随着支付能力和支付意愿的迅速提升,研究区林地生态系统服务动态价值将会显著增加。林地生态系统服务动态价值增长速度(89.72%)高于区域生态系统服务动态价值的增长均值(74.01%),表明高生态价值的土地利用类型对社会和环境调节能力增长更为敏感。