吴爱林, 陈 燕, 燕彩霞, 罗 婵, 宋新山

(东华大学 环境学院, 上海 201620)

长江三角洲生态系统服务价值分析及趋势预测

吴爱林, 陈 燕, 燕彩霞, 罗 婵, 宋新山

(东华大学环境学院,上海201620)

[目的] 对长江三角洲生态系统服务价值进行分析及趋势预测，以期为该区未来城市发展以及制定生态环境建设补偿政策提供支持。 [结果] 以长江三角洲16个地市2000—2014年的土地利用数据和相关统计资料为数据源，采用Costanza生态系统服务价值计算方法，参照修正后的生态系统单位面积生态服务价值系数，结合敏感性分析，探讨生态系统服务价值的时空演化规律，并利用马尔科夫模型对其未来3 a的发展趋势进行了预测分析。 [结果] 区域生态系统总服务价值从2000年的1 303.89亿元减少到2014年的1 294.25亿元，变化率为-0.74%，水源涵养、废物处理和娱乐文化生态系统服务功能价值上升；在空间上，常州、南京、宁波、泰州、台州、扬州的生态系统服务价值不断增加，而其他城市的则处于下降趋势；对长江三角洲2015和2017年的土地利用构成和生态系统服务价值进行预测，呈现耕地、林地持续减少，建设用地持续增加的趋势，生态系统总服务价值将减少276.65亿元。 [结论] 在未来土地调整中应更多地考虑到生态服务价值的重要意义，格外注意对耕地、水域的保护，控制建设用地占用耕地和林地，继续加大植树造林力度，提高土地利用集约度，维持生态系统的功能。

长江三角洲; 土地利用; 生态系统服务价值; 敏感性分析; 马尔科夫模型

文献参数： 吴爱林, 陈燕, 燕彩霞, 等.长江三角洲生态系统服务价值分析及趋势预测[J].水土保持通报，2017,37(4)：254-259.DOI:10.13961/j.cnki.stbctb.2017.04.043； Wu Ailin, Chen Yan, Yancai Xia, et al. Ecosystem service value analysis and trend prediction in Yangtze River delta[J]. Bulletin of Soil and Water Conservation, 2017,37(4)：254-259.DOI:10.13961/j.cnki.stbctb.2017.04.043

生态系统服务功能，是生态系统、生态过程形成的自然环境条件及其效用[1]，是维持人类生存的根本。近几年，生态服务价值已经成为可持续发展研究的热点，普遍研究者已将生态服务价值作为评价生态环境变化的主要指标之一[2-3]，生态系统服务价值(ecosystem service value， ESV)的定量评估已成为生态经济学和环境经济学的研究重点[4-5]。对生态系统服务价值及其变化进行评估，不仅反映了人类对自然资源周期性利用情况，更重要是突出在经济发展过程中人类对生态系统平衡的重要性的认知[6]。生态系统服务价值是指通过生态系统的结构、过程和功能直接或间接提供的生命支持产品(如原材料和食物等)和服务(如栖息地提供等)，包括人类生活所必需的生态产品和保证人类生活质量的生态功能两部分[7]。国内外学者已经分别从全球、国家、区域等不同尺度以及单个生态系统服务价值的评估进行了研究，其中Costanza等发表的“全球生态系统服价值和自然资本”对之后的相关研究产生了深远的影响，使生态系统服务价值评估的原理与方法从科学意义上得以明确。本文基于ENVI软件平台，采用转移矩阵、Costanza等[1]提出的生态系统服务价值计算公式和马尔科夫预测方法，结合敏感度分析，探讨了2000—2014年15 a间不同空间尺度上的生态系统服务价值变化，并对2015—2025年的生态系统服务价值进行预测，以期为长江三角洲未来城市发展以及制定生态环境建设补偿政策提供支持。

1 研究区概况

长江三角洲地区包括杭州、苏州、台州、绍兴、宁波、上海、南通、扬州、湖州、南京、无锡、泰州、常州、镇江、嘉兴和舟山两省一市共16个城市。是中国综合实力最强的经济中心。截至2014年末，长江三角洲国内生产总值是2000年的6.62倍，达到了128 802.8亿元，占全国的20.25%。经济的快速发展，城镇化的快速推进，城市用地类型及空间结构也随之发生变化。研究期间，耕地、林地、草地面积减少，其余类型面积呈增加趋势。其中，耕地的变化幅度最大，面积减少了8 129.45 km2，其次为城镇用地，城镇用地净增加了5 162.18 km2，从而严重影响了区域内的生态系统服务价值。

2 研究方法与数据来源

2.1 研究方法

2.1.1 生态系统服务价值评价方法 本研究以谢高地等[8]提出的“中国生态系统单位面积生态服务价值当量”为基础，结合长三角的实际情况进行修订，根据长三角2000—2014年统计年鉴的相关数据，计算得到长三角2000—2014年的平均粮食产量为2 087.42 kg/hm2，2014年间的平均粮食价格为2.604 6元/kg，参考谢高地[8]关于单位生态服务价值当量因子的经济价值量等于当年全国平均粮食产量市场价值的1/7的方法，计算出长三角一个当量因子的经济价值约为776.82元/hm2。由此计算出长三角不同生态系统单位面积生态服务价值(表1)。结合相关土地利用类型面积，根据表1对长三角生态系统服务价值进行评估。其中耕地和农田对应，林地和森林对应，未利用土地和荒漠对应，城镇用地、农村居民点和其他建设用地赋值为零。

(1)

(2)

式中：ESV， ESVf——生态系统服务价值和第f项服务价值；Ai——土地利用类型i的面积； VCfi——单位面积土地利用类型i的第f项服务功能价值系数； VCi——生态系统价值系数。

2.1.2 敏感性分析方法 生态系统价值敏感性指数(coefficient of sensitivity， CS)是为了验证生态系统类型对于土地利用类型的代表性以及生态价值系数的准确性而引入的[9-10]。该指数主要是用来确定ESV对VC的依赖程度，当CS<1，表明ESV对于VC是缺乏弹性的，当CS>1，表明ESV对于VC是富有弹性的。CS越大，表明VC的准确性越关键。本文通过对各土地利用类型的VC进行上下调整50%来说明ESV的变化情况。其计算公式为：

(3)

式中：b,a——初始价值和生态价值系数调整后的价值。

表1 长江三角洲地区生态系统单位面积生态服务价值

2.1.3 马尔科夫预测方法 马尔科夫预测过程是一种具有“无后效应”的特殊随机过程[11]，它表示随机过程X(t)时刻(n+(1)的状态概率分布只与时刻n的状态概率分布有关，而与时刻n之前的状态概率分布无关[12]。可用转移概率矩阵P和初始分布S(0)并利用如下公式对土地利用状态进行预测[13-16]：

(4)

式中：Pij——i类土地利用类型转化为j类土地利用类型的概率，Pij满足2个条件：

平稳马尔柯夫链在时刻T+1的状态采用下式计算

S(T+(1)=S(T)P

(5)

2.2 数据来源

以长江三角洲为研究区域，以2000和2014年Landsat TM卫星影像为主要数据源，数据预处理包括数据读取、辐射定标、大气校正(FLAASH模型)、波段合成、影像镶嵌、研究区裁剪等。将长江三角洲的土地利用划分为6大类17小类，分别是耕地(包括水田和旱地)；林地(包括有林地、灌木林、疏林地和其他林地)；草地(包括高覆盖度草地、中覆盖度草地和地覆盖度草地)、水域(包括河渠、湖泊、水库坑塘、滩涂和滩地)；城乡、工矿、居民用地(包括城镇用地、农村居民点和其他建设用地)；未利用土地。解译精度约为89.4%，并以当年的植被图、地形图、地貌图、等专题图为依据进行数据校验。其他数据包括长江三角行政边界、社会经济数据等。图像处理及数据分析软件分采用ArcGIS 10.2，ENVI5.1，Matlab等。

3 结果与分析

3.1 生态系统服务价值计算

根据公式(1)，可计算得到2000和2014年长三角生态系统服务总价值量(表2)。结果表明：长三角15 a生态系统服务总价值减少了9.65亿元，下降了约0.74%，2000年的生态系统服务总价值为1 303.89亿元，2014年为1 294.25亿元。其中，林地对生态系统服务价值总量的贡献最大，分别为41.15%和41.01%，而未利用地的贡献最小。从各类用地生态系统服务价值的总体变化情况来看，水域和未利用土地的生态系统服务价值呈增长趋势，水域的增加量明显，占增加总量的93.95%；耕地、林地和草地生态系统服务价值呈降低趋势，以耕地的降低量最大，占三者降低总量的86.56%，其次是林地，草地最低。研究期内，水域生态系统服务价值占长江三角生态系统服务总价值的比例增加了3.4%，其生态系统服务价值所占比例均达30%以上，但其面积比例却远小于此，2000，2014年长江三角洲水域面积比例分别为11.29%和12.18%。表明生态价值系数较高的水域面积增加已经弥补不了由于耕地和林地减少引起的生态系统服务总价值的下降，使得长江三角洲生态系统服务总价值近年来有所下降。长江三角洲在今后土地利用中应注重对林地和耕地生态系统的保护，继续加强生态环境建设。

表2 长江三角洲地区生态系统服务价值(ESV)变化

3.2 单项生态系统服务总价值计算

根据公式(2)，可计算得到2000和2014年长三角单项生态系统服务价值量(表3)。可以看出研究期内生态系统各项服务价值的结构顺序为:水源涵养>废物处理>土壤形成与保护>生物多样性保护>气候调节>气体调节>娱乐文化>原材料>食物生产。前4项生态系统服务功能都占据各时期内所有功能的10%以上，并且总和也超过了生态系统服务功能的60%，说明水源涵养、废物处理、土壤形成与保护、生物多样性保护是长三角地最主要的生态系统服务功能。水源涵养、废物处理和娱乐文化呈不断增加的趋势，尤其是水源涵养功能增加程度最显著，达到13.3亿元，占增加量的60.56%。这与长三角水系发达、水量充沛密切相关。土壤形成与保护价值损失最大(10.5亿元)，气候调节次之(6.05亿元)，其损失均主要由林地和高覆盖地草地退化所致。

表3 长江三角洲地区生态系统单项服务功能价值(ESV)变化

3.3 生态系统服务价值空间分异变化

以长三角16地级市为单位，对生态系统服务价值变化量和变化率进行计算和比较，结果详见表4。可以看出，在2000—2014年14 a间，杭州、湖州等10个市的生态服务价值处于下降趋势，其中杭州市减少最多，达到5.57亿元，其次是上海、绍兴、无锡等市。其他地市生态系统服务价值均有不同程度的增加。其中，宁波市的增加量最大，为9.92亿元，占全省生态系统服务价值总增加量的53.94%，年均增加0.71亿元，其次是台州市，增加了4.97亿元，占全省生态系统服务价值总增加量的27.03%，年均增加0.36亿元，然后依次为泰州市、南京市、扬州市，常州市的增加量最少，为0.02亿元从变化幅度来看，生态系统服务价值变化率下降最快的是嘉兴市，达到了-9.07%，其次是上海市、无锡市，分别达到了-5.46%，-5.17%；其余地级市的变化幅度均低于±5%。

3.4 敏感性分析

根据式(3)，通过对各土地利用类型的VC进行上下调整50%的生态系统价值系数来计算敏感性指数，结果详见表5。由表5可以看出：ESV对VC的敏感性指数都小于1，由高到低依次为林地、水域、耕地、草地、未利用土地，其中，草地敏感性指数约为0.01，说明当草地的VC增加1%时，ESV将增加0.01%。这表明，研究区域ESV对VC缺乏弹性，研究结果是可信的，可为长三角土地利用规划提供借鉴和参考。

3.5 生态系统服务价值的马尔科夫预测

运用马尔科夫过程进行预测的关键在于转移概率的确定。本研究先用ENVI 5.1软件中的Class Statistics工具统计长三角2014年土地利用类型面积，形成初始状态矩阵，再用ENVI5.1软件对长三角2000和2014年的遥感影像分类结果图进行转移矩阵计算，获得了各种土地利用类型之间的转化关系(表6)。由表6可以看出，2014年长江三角洲城镇用地面积共增加了5 345.60 km2。有4 166.05 km2来自耕地，92.52 km2来自林地，11.26 km2来自草地，134.35 km2来自水域，826.60 km2来自农村居民地，还有112.76，2.06 km2分别来自其他建设用地和未利用土地。同样了解到林地、草地、水域等土地利用类型的转变及其构成。其中2014年耕地面积减少最多，主要转化为城镇用地和农村居民用地。利用已求出的2000—2014年内的土地利用各类型面积转移矩阵，求出该时间段内各个土地利用类型的年平均转移概率，构成1个年平均转移概率矩阵，该矩阵即为初始状态转移矩阵。利用P(n)=Pn(n=1,2,3)求得最终转移概率矩阵。最后用初始状态矩阵与转移概率矩阵相乘(运算过程采用Matlab软件进行)，即可得最终状态矩阵，即2015，2020，2025年预测的长三角各土地类型面积，并求得未来10 a生态系统总服务价值。

表4 长江三角洲地区各市生态系统服务价值空间分异

表5 调整生态价值系数后生态系统服务总价值(ESV)的变化情况(VC)及敏感性指数(CS)

表6 长江三角洲地区2000-2014年土地利用各类型转化情况 km2

根据表7对长三角2015-2017年各类土地利用面积变化进行预测，可以看出城镇用地、农村居民点、未利用土地面积出现不同程度增加，变化率分别达到42.87%，15.20%，21.02%。这表明随着经济高速发展、城市化进程加快，长江三角洲建设用地占用土地的强度还将不断加大。

由预测结果可以看出，长江三角洲2015—2017年总的生态系统服务价值呈下降的趋势，从1 223.5亿元减少至946.85亿元，共减少276.65亿元。长江三角洲不同土地利用类型的生态系统服务价值变化规律略有不同:2015—2017年内，仅有未利用土地的生态系统服务价值不变，其他的均出现下降。其中，耕地的生态系统服务价值量减少最多，达到46.20亿元，占减少总量的71.74%；其次是林地，减少量达到11.95亿元，占减少总量的18.56%；水域、草地的生态系统服务价值分别减少了5.03,1.22亿元，所占土地生态服务价值减少总量的比例分别为7.81%和1.89%。

表7 研究区2015,2020，2025年各类土地利用面积及生态系统服务价值预测

4 结 论

2000—2014年长三角土地生态系统服务价值总体呈不断下降的时序变化趋势，年下降率为0.74%。这与长三角耕地、林地、草地面积的大幅度减少，引起土壤形成与保护、食物生产、气候调节等价值快速下降密切相关。即使高生态价值的水域面积增加也弥补不了由于耕地、林地、草地减少引起的生态系统服务总价值的下降。在空间上，杭州市、苏州市、台州市的生态系统服务价值在各地级市排序中排前3位，镇江市、嘉兴市、舟山市的生态系统服务价值最低，尤其是嘉兴市在研究期间生态系统服务价值的下降率最高，迫切需要保护耕地、水域，继续加大植树造林力度，维持生态系统的功能。

根据马尔柯夫过程对2015—2017年长三角的土地利用结构和生态系统服务价值进行了预测。未来土地利用生态系统服务价值变化趋势与研究区主要构成者耕地、林地、水域的变化趋势保持一致，研究区域生态系统总服务价值将减少276.65亿元。因此，在未来土地调整中应更多地考虑到生态服务价值的重要意义，格外注意对森林生态系统的保护，控制建设用地占用耕地和林地，提高土地利用集约度。

[1] Daily G C, Myers J P, Reichert J. Nature’s service:societal dependence on natural ecosystems[M].Washington: Island Press, 1997.

[2] Seppelt R, Dormann C F, Eppink F V, et al. A quantitative review of ecosystem service studies: Approaches,shortcomings and the road ahead [J].Journal of Applied Ecology, 2011,48(3):630-636.

[3] Fisher B, Turner R K, Morling P. Defining and classifying ecosystem services for decision making [J].Ecological Economics, 2009,68(3):643-653.

[4] Chen Zhongxin, Zhang Xinshi. Value of ecosystem services in China[J]. Chinese Science Bulletin, 2000,45(10):870-876.

[5] Pearce D. Auditing the Earth: The value of the world’s ecosystem services and natural capital[J]. Environment, 1998,40(2):23-28.

[6] 黄湘,陈亚宁,马建新.西北干旱区典型流域生态系统服务价值变化[J].自然资源学报,2011,26(8):1364-1376.

[7] Costanza R, Arge R, Groot R, et al. The value of the world’s ecosystem services and natural capital[J]. Nature, 1997,387(7):253-260.

[8] 谢高地,鲁春霞,冷允法,等.青藏高原生态资产的价值评估[J].自然资源学报,2003,18(2):189-196.

[9] Li Renqiang, Dong Ming, Cui Jianyong, et al. Quantification of the impact of land-use changes on ecosystem services: A case study in Pingbian County, China[J]. Environmental Monitoring & Assessment, 2007, 128(4):503-510.

[10] Kreuter U P, Harris H G, Matlock M D, et al. Change in ecosystem service values in the San Antonio area, Texas[J]. Ecological Economics, 2001,39(3)：333-346.

[11] Hulst R. On the dynamics of vefetation: Markovchains as models of succession[J]. Vegetation, 1979,40(1):3-14.

[12] 胡汝晓,赵松义,谭周进,等.烟草连作对稻田土壤微生物及酶的影响[J].核农学报,2007,21(5):494-497.

[13] 王晓峰,任志远,黄青.农牧交错区县城土地利用变化及驱动力分析:以陕北神木县为例[J].干旱区地理,2003,26(4):402-407.

[14] 郜鲁豪,徐旌.基于Markov模型的安宁市土地利用预测[J].云南地理环境研究,2010,22(2):87-89.

[15] 蓝永超,丁永建,康尔酒,等.黑河流域水资源动态变化及其趋势的灰色Markov链预测[J].中国沙漠,2003,24(3):435-440.

[16] 李黔湘,王华斌.基于马尔科夫模型的涨渡湖流域土地利用变化预测[J].资源科学,2008,30(10):1542-1543.

Ecosystem Service Value Analysis and Trend Prediction in Yangtze River Delta

WU Ailin, CHEN Yan, Yancai Xia, LUO Chan, SONG Xinshan

(College of Environment, Donghua University, Shanghai 201620, China)

[Objective] This paper examined the spatio-temporal variations in the value of ecosystem services and then predicted it to provide support for formulating compensation policy in ecological environment construction and future urban development. [Methods] Based on land use data and relevant statistical data in the Yangtze River delta from 2000 to 2014, with the method proposed by Costanza wherein the values of ecosystem services are calculated with regard to the value coefficients of ecosystem services per unit area in the Chinese terrestrial ecosystem and the sensitivity analysis, this paper examined the spatio-temporal variations in the value of ecosystem services and then predicted it by Markov model. [Results] The total value of ecosystem services decreased from 130 389 million yuan in 2000 to 129 425 million yuan, reduced by 0.74%. In addition, the ecosystem service value of water supply, recreation and culture function and waste treatment increased while other ecosystem functions declined. From the perspective of spatial analysis, the ecosystem service values of Changzhou, Nanjing, Ninbo, Taizhou, Taizhou, Yangzhou City increased, while values of other cities decreased. The land use structures of the Yangtze River delta in 2015 and 2017 were predicted using the Markov model. The trends of woodland and cultivated land will maintain a continual increase while the construction land will be in a continual decrease in the coming three years. The total ecosystem services value will have a decrement of 276.65 million yuan. [Conclusion] Therefore, the importance of ecological service value should be given more consideration in future land adjustment, especially for the protection of cultivated land and water and improvement of land utilization intensity. it is very important to control the occupancies of cultivated and forest land by construction, to continually increase the effort of afforestation and maintain the function of ecosystem.

YangtzeRiverDelta;landuse;ecosystemservicevalue;sensitivityanalysis;Markovmodel

A

： 1000-288X(2017)04-0254-06

： F301.24， Q149

2016-10-15

：2016-12-22

国家自然科学基金项目“长江三角洲区域生态安全时空演变及调控图谱研究”(41471089); 中央高校基本科研业务费项目

吴爱林(1990—),女(汉族),江西省玉山县人,硕士研究生,研究方向为长江三角洲的土地生态敏感性评价。E-mail:270134083@qq.com。

陈燕(1983—),女(汉族),上海市人，博士，副教授，主要从事生态安全、纳米技术(材料)在水处理中的应用以及环境遥感与环境信息系统。E-mail：chenyanhjxy@dhu.edu.cn。