［提要］　土地利用及生态系统服务价值（ESV）是衡量经济发展和生态保护的重要指标，该研究有利于促进社会经济和生态环境协调发展。以广州市为研究对象，以土地利用变更数据为基础，对研究区生态系统服务价值进行评估，并且分析其时序变化。结果表明：2000～2020年广州市生态系统服务价值总体呈现下降趋势，林地提供最大比例的生态系统服务价值；广州市esv变化幅度逐渐收窄。当前广州市生态状况较为稳定，但未来仍然要注意保持生态与经济的协调发展，合理规划用地，加强生态保护。

关键词：广州；土地利用；生态系统服务价值

中图分类号：F127；X171.1 文献标识码：A

收录日期：2022年2月28日

生态系统作为人类在生存发展进程中无可替代的自然资产，其服务是通过生态系统的结构、过程和功能直接或间接得到的生态支持产品和服务，关系着生态系统和人类福祉。其中，生态系统服务价值（ESV）是量化和评价生态系统服务能力的重要方法，土地利用作为生态系统服务价值的重要体现，近年来城镇化高速发展使土地利用剧烈变化，土地盲目开发、水土流失、空气污染等环境问题频发，导致城市生态环境遭受破坏，损害了城市生态系统服务价值。通过探究土地利用变化与生态系统服务价值之间的关联，能够更好地从生态系统的整体角度保护生态系统结构、功能与资源的可持续供给能力，对促进社会经济和生态环境协调发展具有重要意义。

近年来，国内外学者对土地利用转移和生态系统服务价值进行了广泛的研究，逐渐聚焦于两者的耦合研究。由于不同区域生态系统组成的多样性，使土地利用变化对生态系统服务价值影响的强度、方式及作用效果等各不相同，研究视角也从区域向地方转变，较小尺度和典型区域的分析让研究问题更加聚焦。另外，在生态系统服务评价理论和评价计算方法上取得了重要进展。余建华等在Costanza研究基础上，构造了土地利用变化生态效应指数等指标定量测度研究区土地利用变化的生态效应。谢高地等结合中国实际情况对生态系统服务价值的评估方法进行了多次修正，得出了适合中国国情的生态系统服务价值的基础当量表，得到了广泛应用。

基于此，本文以广州市为研究区域，以2000～2020年为研究时段，基于区域土地利用变化评估其生态系统服务价值变化，分析土地利用与生态系统服务价值的协同发展状况及时序变化差异，为广州市生态环境保护、生态功能区划以及经济社会可持续发展提供建议，同时也为同类城镇化快速发展的城市生态系统服务价值评估提供借鉴。

一、研究区概况和数据准备

（一）研究区概况。广州市位于珠江三角洲北缘（112°57′E～114°3′E，22°26′N～23°56′N），地势东北高、西南低，属于海洋性亚热带季风气候，土地类型多样，适宜性广。境内河流水系发达，河流众多，水域面积广阔。全市土地总面积为7，434.4平方千米，现辖越秀、海珠、荔湾、天河、白云、黄埔、花都、番禺、南沙、从化、增城共11个区。自21世纪以来，广州市经济实现飞速发展，地区生产总值（GDP）从2000年的2，375.9亿元上升至2021年的28，231.9亿元，加上城镇化进程的不断推进，2021年常住人口高达了1，867.66万人。随着经济社会的发展，广州市建设用地面积从2000年的825.15平方千米增至2020年的1，633.9平方千米，年均增长率为6.2%。建设用地的不断扩张，使得生态用地结构逐渐缩小，生态环境问题成为当前广州市经济进一步发展的瓶颈。

（二）数据来源与处理。选用2000版、2010版和2020版GlobeLand30全球地表覆盖数据，该数据是由国家地理信息中心提供的精度为30米全球地表土地利用数据，采用的是WGS84坐标系统。分类体系参照《土地利用现状分类》（GB/T21010-2017），并根据广州市实际土地利用使用情况，分为耕地、林地、草地、水域、建设用地和未利用地共6种地类，得到广州市2000年、2010年、2020年的土地利用分类图，分类后经Kappa系数和Google Earth检验，解译结果较为理想。

二、研究方法

（一）土地利用转移矩阵。土地利用转移矩阵用于反映在研究时段初至末广州各用地类型面积间相互转化的动态过程。公式如下：

Sij=■

式中，S为面积；n为广州地区土地利用类型数；ij（i，j=1，2，…，n）分别是发生土地转移前后所代表的土地类型；Sij为i与j地类间转换的面积量。

（二）生态系统服务价值估算。本研究参照Costanza等和谢高地等的研究成果，估算研究区生态系统服务价值，划分一级类型4个，分别是供给服务、调节服务、支持服务和文化服务；分支的二级类型10个，分别是食物生产、原材料生产、水资源供给、气体调节、气候调节、水文调节、废物处理、保持土壤、维持生物多样性、提供美学景观等10项生态系统服务功能价值，主要来源于2000～2020年《广州市统计年鉴》《广东省统计年鉴》《广州市国民经济和社会发展统计公报》。考虑到地域和时间的差异，根据研究区实际情况对谢高地等的“中国生态系统单位面积生态系统服务价值当量表（2007）”进行修正，将建设用地生态系统服务价值定义为0，最终确定符合广州城市发展实情与用地特征的单位面积生态系统服务价值表。（表1）

三、结果分析

（一）土地利用变化分析。2000～2020年各类用地类型相互转换，转出量最高的是耕地，其次是林地，两者转出的主要方向均为建设用地，其中约有82.31%的耕地面积转换为建设用地，林地的转出量仅次于耕地，为138.23km2。建设用地是研究时段内转入量最大的地类，除了耕地和林地外，水域土地也是建设用地面积增加的主要贡献者之一，其转入建设用地的面积占该地类面积的40.58%，草地对建设用地面积增加的贡献程度仅次于水域用地，草地向建设用地转出面积约占草地面积的50.23%。除建设用地外，其余用地类型的转出量明显大于转入量，呈现“入不敷出”的态势。（表2）

（二）生态系统服务价值变化分析。根据广州市单位面积生态系统服务价值表，估算出广州市2000～2020年各类用地的生态系统服务价值，并对生态系统服务价值时空变化特征进行分析。（表3）

广州市2000～2020年生态系统服务价值（下文简称ESV）总体呈下降趋势。2000～2010年广州市ESV的变化幅度较大，总值有明显减少。2010～2020年广州市ESV的降幅收窄，但与2000年相比，ESV总值下降约5%。与2000年相比，2010年的ESV下降主要是因为各地类均有较大面积转入到建设用地，各项用地的ESV下降明显；2020年变化原因与2010年相似，除了林地和草地保持稳定，其他地类的ESV均下降。由表3可知，耕地和林地ESV的变化是影响整个研究时段ESV变化的主要因素。在所有地类中，林地是广州市ESV构成的主体，其ESV占全市总量的43%；其次是水域，约占全市总量的33%；第三位的是耕地，约占全市总量的14%，其余地类的ESV占比较小。2000～2020年间，耕地、未利用地的ESV降幅最大，而林地、草地和水域的ESV整体保持平稳，这一变化与土地利用变化状况基本吻合。

从各类生态系统服务价值类型来看，广州市的ESV以水文调节和气候调节为主，最少的服务是原材料生产和提供美学景观。20年间，广州市ESV持续下降，各项服务类型均有变化，其中单项ESV变化最大的是水文调节，减少了1.8亿元；其次是水资源供给，减少了1.4亿元，整体而言，水域和林地是广州市ESV的主要提供者，广州市水系密布，流域面积稳定，还有丘陵地区较高的林地覆盖率，都较好地维持了区域ESV的稳定。（表4）

四、结论与讨论

（一）结论。研究利用GIS相关知识技术，借鉴国内外相关研究，通过土地利用转移矩阵和生态系统服务价值评价方法，对2000～2020年广州市各类型用地进行分析，结果表明：（1）2000～2020年，广州市建设用地面积显著增加，其余用地类型面积减少趋势稳定，其中耕地和林地面积负向变化最为显著。各用地类型间均存在类型转换，其中转出量最高的是草地，转入量最高的则是建设用地，除建设用地外，其余用地类型均呈现“入不敷出”态势。（2）2000～2020年，广州市ESV总价值呈下降趋势，下降速度由迅速到缓慢的变化，调节服务是广州市最主要的生态系统服务功能，其次则是支持服务功能。2000～2010年，水域、林地和耕地的生态系统服务价值变化度均呈负向变化，其中未利用土地的变化度最大，耕地次之。2010～2020年，林地和未利用土地面积减少的趋势得到缓和。

（二）讨论。研究中方法的确定及结果的分析均具有合理性与正确性，基于生态系统服务价值基础上确定广州市各类用地的生态系统服务价值，分析其演变结果与前期学者的研究结果相一致，研究结果为城镇化快速发展及类似的沿海城市的社会经济发展提供参考和借鉴。

研究并未对区域土地利用转型驱动力进行深入研究，后续研究中需进一步基于区域性，以不同空间、时间尺度来探讨土地利用转型的驱动机制，更好地揭示并解决相关问题，并进一步完善生态系统服务价值评估体系，进行更精准的评估。

（作者单位：广东工业大学管理学院）

主要参考文献：

［1］谭术魁，韩思雨，周敏.土地城市化背景下武汉市资源环境承载力仿真研究［J］.长江流域资源与环境，2017.26（11）.

［2］谢高地，张彩霞，张昌顺，肖玉，鲁春霞.中国生态系统服务的价值［J］.资源科学，2015.37（09）.

［3］胡和兵，刘红玉，郝敬锋，安静.城市化流域生态系统服务价值时空分异特征及其对土地利用程度的响应［J］.生态学报，2013.33（08）.

［4］吴建寨，李波，张新时，夏艳玲，崇洁.天山北坡土地利用/覆被及生态系统服务功能变化［J］.干旱区地理，2007（05）.

［5］刘永强，廖柳文，龙花楼，秦建新.土地利用转型的生态系统服务价值效应分析——以湖南省为例［J］.地理研究，2015.34（04）.

［6］刘兴诏，林丽丽，董建文，陈颐，王齐，张增可，黄柳菁.国家森林公园生态系统服务评估——以福州国家森林公园为例［J］.石河子大学学报（自然科学版），2019.37（05）.

［7］肖红克，李洪忠，王莉，陈劲松，韩宇.粤港澳大湾区土地利用及生态系统服务价值的变化——以广佛肇为例［J］.水土保持研究，2020.27（01）.

［8］李斐，申丽娟，王三，谢德体.土地规划中建设用地变化对耕地保护的潜在影响——以重庆市北碚区为例［J］.西南大学学报（自然科学版），2017.39（09）.

［9］李锋，叶亚平，宋博文，王如松.城市生态用地的空间结构及其生态系统服务动态演变——以常州市为例［J］.生态学报，2011.31（19）.

［10］孙贤斌，孙良萍，王升堂，李光耀，张欣.基于GIS的跨流域生态补偿模型构建及应用——以安徽省大别山区为例［J］.中国生态农业学报（中英文），2020.28（03）.

［11］喻建华，高中贵，张露，彭补拙.昆山市生态系统服务价值变化研究［J］.长江流域资源与环境，2005（02）.

［12］吕乐婷，张杰，彭秋志，任斐鹏，江源.东江流域景观格局演变分析及变化预测［J］.生态学报，2019.39（18）.

［13］翟浩然，唐新明，王光辉，李佳雨，邢丽玮.中国超大城市建设用地遥感监测与时空格局［J］.遥感信息，2019.34（04）.

［14］Costanza R，d'Arge R，de Groot，Farber S，Grasso M，Hannon B，Limburg K，Nacem S.The value of the world's ecosystem services and natural capital［J］.Nature，1997（387）.

［15］李薇，谈明洪.太行山区不同坡度NDVI变化趋势差异分析［J］.中国生态农业学报，2017.25（04）.

［16］谢高地，甄霖，鲁春霞，肖玉，陈操.一个基于专家知识的生态系统服务价值化方法［J］.自然资源学报，2008（05）.