上海市土地利用变化对生态系统服务价值的影响研究.pdf

上海国土资源Shanghai Land&Resources 2023Vol.44.3 79doi:10.3969/j.issn.2095-1329.2023.03.013上海市土地利用变化对生态系统服务价值的影响研究王丹妮（上海市地质调查研究院，上海 200072）摘 要：面对资源约束趋紧、环境污染严重、生态系统退化的严峻形势，建设生态文明是当前改善生态环境的迫切需要。研究土地利用变化对地区生态系统服务价值的影响，能够为土地资源的可持续利用与生态环境保护建设提供科学依据。本文以上海市（陆域）为研究对象，基于 20102020 年土地利用变更数据，采用当量因子法和 GIS 空间分析技术，分析上海市 2010、2020 年两期的土地利用变化特征与生态系统服务价值时空变化的关系。结果表明：（1）20102020 年，上海市土地利用格局变化明显，除草地和林地的面积增长外，其他各类用地面积均有不同程度的减少，其中水域面积减少最多，林地面积增加最多。（2）20102020 年间上海市生态系统服务价值减少 52.65 亿元，主要是因为水域面积减少所致。（3）在 500 m500 m 格网尺度下，上海市生态系统服务价值在空间上总体表现为高高、低低集聚，局部呈现高低、低高集聚的特点，人类活动是其主要驱动因素。关键词：土地利用；生态系统；服务价值；空间分析；驱动因素中图分类号：F293.22；X171.1 文献标志码：A 文章编号：2095-1329(2023)03-0079-06生态系统服务是指生态系统与生态过程中所形成及所维持的人类赖以生存的自然环境条件与效用。量化生态系统服务价值（ecosystem service value,ESV）有利于了解生态系统的供给情况与生态效益的大小。近年来，随着人类活动的不断加剧，土地不合理利用导致的生态环境、资源利用等严重问题不断涌现，也对生态系统服务价值产生了一系列影响。随着 2001 年联合国组织千年生态系统评估项目的启动，国内外掀起了一股生态系统服务价值相关的研究热潮1-3。目前，国内外对于土地利用变化与生态系统服务价值的动态研究多从土地利用方式、类型、格局变化的响应角度出发4，各研究者之间仅研究尺度略有不同，从国家、流域到省市、县域5-7。在城市方面，有学者基于乐山市 2010、2017 年遥感影像，探析乐山市土地利用变化对 ESV 空间分布差异的影响，结果表明乐山市 ESV 的变化呈现极化效应，需进一步加强区域生态管控8；一些学者发现在 1 km1 km 格网尺度下，南昌市 ESV 在空间上表现为高高、低低集聚，自然资源禀赋与人类活动是其主要驱动因素9。对于上海市，有学者基于 Costanza 的生态系统服务单价和我国国情，评估了上海中心城区 19472006 年土地利用变化对区域生态系统服务功能的影响，结果表明 60 年间生态系统服务功能呈现 4 个发展阶段，除水域和耕地的生态系统服务功能显著下降外，其余地类均有上升10；另一些学者同样基于 Costanza 模型，结合谢高地等的研究成果，计算了 20032007 年间上海市 ESV 的变化情况，结果显示ESV 对价值系数的变化是缺乏弹性的11。这两项针对上海市的研究均缺乏土地利用变化与ESV的空间对应研究，故本文以上海市（陆域）为研究对象，基于 20102020年土地利用变更数据，采用当量因子法和 GIS 空间分析技术，分析上海市 2010、2020 年两期的土地利用变化特征与生态系统服务价值时空变化的关系，旨在为土地资源的可持续利用与生态环境保护建设提供科学依据。1 数据来源与研究方法1.1 研究区概况上海市位于中国华东地区，平均海拔高度 2.19 m。上海市包括黄浦区、徐汇区、长宁区、静安区、普陀区、虹口区、杨浦区、闵行区、宝山区、嘉定区、浦东新区、收稿日期：2023-06-30修回日期：2023-08-08作者简介：王丹妮，1988 年生，女，硕士，工程师，主要从事土地质量监测与评价等研究。电子邮箱：基金项目：上海市规划和自然资源局科研项目“上海市土地质量调查与监测”“上海市土地质量环境质量监测网优化研究”上海国土资源 Shanghai Land&Resources80 2023Vol.44.3 金山区、松江区、青浦区、奉贤区、崇明区的 16 个区，常住人口约为 2475.89 万人。1.2 数据来源研究数据来源于上海市 2010 年、2020 年土地利用变更数据，依据中国科学院土地利用/覆被分类标准，将各细分土地利用类型归并为耕地、林地、草地、水域、建设用地和未利用地 6 种土地利用类型（图 1）。社会经济数据主要来自上海统计年鉴（20112021 年）。1.3 研究方法（1）生态系统服务价值当量的确定采用我国学者提出的中国生态系统价值系数12，结合上海市具体情况，依据上海市 20102020 年粮食作物产值、粮食单产数据，综合相关文献和区域自然情况，采用一个标准生态服务当量子的经济价值等于研究区单位面积农田粮食产值的 1/7 的中国生态系统服务系数修正方法13，得出上海市一个生态系统服务价值当量为 2491.64元/hm2。（2）基于格网的 ESV 空间表达结合上海市地形特征、面积大小等因素，采用 500 m500 m 格网作为评价单元，利用 ArcGIS 10.6 软件确立覆盖上海市陆域的格网。格网确定后，以每个格网内土地利用类型数据为基础，计算出各个格网中的 ESV：nmnmEALU=（1）1ksnnEE=（2）式中：nE 是每一个格网单元内第 n 类土地利用类型的生态系统服务价值；sE 是每一个格网单元内生态系统服务价值总量；mnA为第 m 个网格中第 n 种土地利用类型的面积（2hm）；mLU 为第 m 种土地利用类型的单位面积生态系统服务价值；k 为土地利用类型个数。（3）空间自相关分析本文中的全域空间自相关是利用 Morans I 指数来判断上海市 ESV 在空间上是否有聚集特性存在，局域空间自相关是利用 Anselin Local Morans I 指数判断上海市ESV 在空间上的聚集空间单元及其空间自相关性的显著性。2 土地利用结构变化2.1 土地利用数量变化表 1 反映出 20102020 年上海市主要土地利用类型一直是建设用地、耕地，2010 年、2020 年这两种地类总面积分别占上海市总面积的 70.81%、66.82%。从变化量可知，林地和草地面积分别增加 39270.09 hm2、17263.71 hm2，林地增幅为 57.29%，草地增长了 10.52 倍；水域和耕地面积大幅减少，面积分别缩减 28287.43 hm2、19245.35 hm2，降幅分别为 21.61%、10.14%，建设用地面积缩减 8474.49 hm2，降幅达 2.81%。2.2 土地利用转移分析由表 2 可知，20102020 年间上海市耕地的转出面积最多，转出 67010.97 hm2，主要转变为林地、建设用地，面积分别为 35686.95 hm2、21318.34 hm2；其次是建设用地，转出 49405.24 hm2，主要转变为林地、耕地，面积分别为 18502.13 hm2、14375.76 hm2；水域转出面积达 45698.66 hm2，主要转变为耕地、建设用地、林地，面积 分 别 为 20202.66 hm2、10709.78 hm2、10407.45 hm2。林地的转入面积最多，为 65164.72 hm2，主要由耕地、建表 1 20102020 年上海市土地利用面积和变化率Table 1 The land use area and rate of change in Shanghai from 2010 to 2020 时间项目耕地林地草地水域建设用地未利用地2010 年面积（hm2）189739.2168544.001640.19130897.38301268.741302.26占比（%）27.369.890.2418.8843.450.192020 年面积（hm2）170493.86107814.0918903.91102609.96292794.24775.72占比（%）24.5915.552.7314.8042.230.1120102020 年面积（hm2）-19245.3539270.0917263.71-28287.43-8474.49-526.55变化率（%）-10.1457.291052.54-21.61-2.81-40.43图 1 不同时期上海市土地利用现状图Fig.1 Land use status map of Shanghai in different periods上海国土资源Shanghai Land&Resources 2023Vol.44.3 81设用地、水域转入，面积分别为 35686.95 hm2、18502.13 hm2、10407.45 hm2；其次是耕地，转入面积为 47765.62 hm2，主要由水域、建设用地、林地转入，面积分别为20202.66 hm2、14375.76 hm2、12512.20 hm2。其中，耕地林地是永久基本农田划定过程中耕地从严判定的体现，建设用地林地、耕地是上海市近 5 年大力推进集建区外低效建设用地减量化，将低效建设用地转为农用或生态用地功能的行动结果，水域耕地、林地多为沿海滩涂围垦造田、郊区以坑塘水面为代表的水域土地整治的结果，而耕地、水域建设用地是上海市城市化进程中集建区建设开发的结果。3 生态系统服务价值变化3.1 生态系统服务价值数量变化通过确立的上海市生态系统服务价值当量（2491.64元/hm2），由此计算出上海市 6 种土地利用类型的单位面积生态系统服务价值（见表 3）。在确立各种土地利用类型的单位面积生态服务价值的基础上，采用生态系统服务价值计算公式，计算上海市土地生态系统服务价值：1niiiEVCA=（3）式中：E 代表研究区生态系统服务价值总量；iVC 代表第i 类土地利用类型的单位面积生态系统服务价值(当量/hm2)；iA 代表第 i 类土地利用类型的面积(hm2)；n 为研究区所有土地利用类型个数。表 4 反映的是 20102020 年上海市土地利用的 ESV变化情况。由表可知，10 年间 ESV 共减少了 52.65 亿表 2 20102020 年上海市土地利用和生态系统服务价值转移矩阵Table 2 Conversion matrix of the change of land use and ESV from 2010 to 2020类别2010 年耕地林地草地水域建设用地未利用地转入2020 年耕地面积（hm2）12512.20315.5420202.6614375.76359.4647765.62ESV（万元）-59421.29-592.80-538624.0413933.69330.50林地面积（hm2）35686.95355.2010407.4518502.13212.9965164.72ESV（万元）169479.681019.56-228047.92105800.981207.36草地面积（hm2）1243.143110.994372.759718.67141.4218586.97ESV（万元）2335.48-8929.71-108366.9727678.22395.71水域面积（hm2）8745.592382.64142.456060.0580.5117411.24ESV（万元）233166.6652208.313530.36167441.052220.20建设用地面积（hm2）21318.347887.28510.0610709.78505.2840930.74ESV（万元）-20662.77-45101.90-1452.63-295914.49-25.18未利用地面积（hm2）16.951.520.006.02748.63773.12ESV（万元）-15.59-8.62-165.9837.31转出面积（hm2）67010.9725894.631323.2545698.6649405.241299.66表 4 20102020 年上海市土地利用的 ESV 变化情况Table 4 Changes of ESV about land use in Shanghai from 2010 to 2020类型ESV（亿元）ESV变化量（亿元）ESV变化率（%）2010 年2020 年2010-2020 年2010-2020 年耕地18.39016.525-1.865-10.14林地39.19661.65122.45557.29草地0.4675.3844.9171052.89水域361.673283.514-78.159-21.61建设用地0.000.000未利用地0.0060.004-0.002-33.33总计419.733367.079-52.654-12.54表 3 20102020 年上海市 6 种土地利用类型的单位面积生态系统服务价值Table 3 The ESV per unit of six types of land use in Shanghai from 2010 to 2020服务类型评价指标耕地（元/hm2）林地（元/hm2）草地（元/hm2）水域（元/hm2）建设用地（元/hm2）未利用地（元/hm2）供给服务食物生产3388.50722.58548.161848.800.000.00原料生产224.321644.48822.24707.630.000.00水资源供给-6552.35847.16448.5017815.230.000.00调节服务气体调节2765.615406.862840.472481.670.0049.83气候调节1420.1816195.667524.756358.670.000.00净化环境423.564808.872491.6412856.860.00249.16水文调节6776.6511810.375506.52215870.710.0074.75支持服务土壤保持25.176602.853463.383004.920.0049.83维持养分循环473.39498.33274.08229.230.000.00生物多样性523.226004.853164.389004.790.0049.83文化服务美学景观224.242641.141395.326124.450.0024.92上海国土资源 Shanghai Land&Resources82 2023Vol.44.3 元，降幅达 12.54%。从土地利用类型上看，水域是上海市 ESV 的主体，在 2010 年 ESV 达 361.673 亿元，占总价值量的 86.17%；其次是林地和耕地，2010 年 ESV 分别为 39.196 亿元、18.390 亿元，占比到达 9.34%、4.38%。到了 2020 年，六种土地利用类型中，林地和草地的 ESV分别增长了 22.455 亿元、4.917 亿元；水域和耕地的 ESV分别减少了 78.159 亿元、1.865 亿元。比较不同生态系统服务类型的 ESV 可知，水文调节的 ESV 占比最高，2010、2020 年分别达到 72.34%、67.24%。主要原因一是提供水文调节的水域基础当量远远大于其他地类，二是在上海市土地利用类型中水域面积较大（见表 5）。比较生态系统服务价值变化量（见表 2）可知，在土地利用转移过程中，耕地水域、耕地林地、建设用地水域、建设用地林地使得 ESV 增加，分别增加 23.317亿元、16.948 亿元、16.744 亿元、10.580 亿元，共占土地利用类型转移过程中 ESV 增加量的 86.57%；水域耕地、水域建设用地、水域林地、水域草地会导致ESV 减少，分别减少 53.862 亿元、29.591 亿元、22.805亿元、10.837 亿元，减少量占土地利用类型转移过程中ESV 减少量的 89.57%。可知，20102020 年间，上海市水域转向其他用地类型导致了 ESV 的减少，而耕地、建设用地转向水域、林地导致了 ESV 的增加，水域的转入转出是引起 ESV 较大变化的主要原因。3.2 生态系统服务价值空间分布变化以 500 m500 m 的格网作为研究单元，计算每个格网的 ESV，并按照自然断点法将上海市 ESV 分为五个等级，一级最低，五级最高。图 2 反映的是 20102020 年上海市 ESV 分布情况，由图可知：（1）2010 年上海市ESV 高值区主要分布于崇明岛的长江北港北支、北湖、东滩、东风西沙、青草沙，宝山区的陈行宝钢水库、青浦区的淀山湖金泽水库一带、浦东新区的南汇沿岸新围垦滩地、奉贤区的杭州湾北岸沿海滩涂以及黄埔江等水域占比高的沿河沿岸地区；2020 年 ESV 高值区在南汇沿岸、崇明地区面积大量缩小，仅在东风西沙地区有所扩大。（2）2010 年较高值区除沿高值区外围分布外，主要集中分布于崇明岛、青浦以及奉贤等郊区中河流水面、坑塘水面密集分布的地区，2020 年较高值区在青浦、奉贤区范围缩小，在崇明岛内整体向岛内、东部转移。（3）中值区与较低值区交错分布，主要分布于郊区的城郊结合带地区，2020 年分布情况基本一致。（4）上海市ESV 低值区主要位于中心城区、郊区集建区等建设用地占比高的城镇化地区，2020 年上海市 ESV 低值区范围进一步扩大，沿已有集建区外围扩张，低值区的增长主要集中在金山南部、浦东新区临港新城地区。3.3 生态系统服务价值空间自相关变化分析表 6 反映的是上海市 ESV 各年份 Morans I 统计值计算结果，由表可知：2010 年、2020 年上海市 ESV 的Morans I 值均大于 0，P 值均小于 0.001，表明上海市表 5 20102020 年上海市生态系统服务价值变化情况 Table 5 Changes of ESV in Shanghai from 2010 to 2020 一级分类二级分类ESV（亿元）ESV变化量（亿元）ESV变化率（%）2010 年2020 年20102020 年20102020年供给服务食物生产9.3538.558-0.795-8.50原料生产2.4923.0370.54521.87水资源供给11.4768.107-3.369-29.36调节服务气体调节12.24913.6281.37911.26气候调节22.24227.8295.58725.12净化环境20.97319.572-1.401-6.68水文调节303.614246.834-56.780-18.70支持服务土壤保持8.56510.9012.33627.27维持养分循环1.5451.6310.0865.57生物多样性16.94917.2040.2551.50文化服务美学景观10.2759.778-0.497-4.84图 2 20102020 年上海市生态系统服务价值分布图Fig.2 Distribution of ESV in Shanghai from 2010 to 2020表 6 上海市生态系统服务价值各年份 Morans I 统计值Table 6 Morans I statistics of ESV in Shanghai from 2010 to 2020年份Morans IZ 值P 值2010 年0.711237.9960.0002020 年0.689230.5600.000上海国土资源Shanghai Land&Resources 2023Vol.44.3 83ESV 整体上呈显著的空间正相关关系；2020 年 Morans I值略小于 2010 年，表明上海市近 10 年来，生态系统服务价值集聚程度降低，相关性逐渐减弱，在空间上表现为整体较为分散、区域集聚程度降低的特点。图 3 是上海市生态系统服务价值 LISA 分布图。由图可知：2010 年上海市 ESV 高高值聚集区除集中于 ESV高值区（五级区）外，还分布于水网密集的横沙岛红星港新民港地区、奉贤区东北部四团港南大港三团港地区以及东南部焦墩港南门港南横港（河）一带。低低值聚集区基本与 ESV 低值成片区分布一致，主要集中于中心城区、崇明岛城桥镇、长兴岛船舶制造基地（东、西块）、宝山区东部及蕴藻浜两岸地区、嘉定区南部、闵行区黄浦江西岸、松江区东北中部地区、奉贤区上海化学工业区、金山区南部杭州湾沿岸、浦东新区西部黄浦江沿岸；2020 年上海市 ESV 空间分布上高高值聚集区范围大幅缩小，低低值聚集区仅在浦东新区、闵行区及嘉定区范围有所缩小。上海市 20102020 的十年间 ESV 聚集情况发生较大转变，主要有以下四种变化方式：（1）非显著区转变成高高值聚集区的 ESV 增值区面积仅占上海市 ESV 增加区总面积的 2.06%，增加量占到了总增加量的 10.44%。此类区域主要分布于崇明岛北湖与东滩东侧向岛地区、东风西沙沿江地区、黄浦江沿岸、奉贤金山北杭州湾沿岸地区以及青浦环淀山湖水网密集区。崇明岛内主要以耕地转化为水域和林地为主；黄浦江沿岸土地多以建设用地转为水域、林地和草地为主；奉贤金山北杭州湾沿岸地区多以建设用地转为水域、林地和草地为主；青浦环淀山湖水网密集区多以耕地、建设用地转化为水域和林地为主（图 4a）。（2）低低值聚集区转变成非显著区中的 ESV 增值区面积占 ESV 增加区总面积的 8.31%，增加量占总增加量的 12.20%。从图 4b 可知，上述区域集中成片区基本与近 10 年上海市重要生态空间建成区重合。浦东新区东北角沿上海绕城高速沿线上海浦东国际机场西侧、浦东新区外环高速三林康桥沿线出现的低低值聚集区转变成非显著区主要是建设用地转化成林地、草地，建成了环城绿带、近郊绿环区；宝山江杨北路泰和路铁峰路兰岗路合围区以建设用地转为转成林地、水域为主，形成吴淞楔形绿地；川杨河、申嘉湖、赵家沟地区形成的聚集区主要是建设用地转化成林地、水域，形成郊区骨干河道两侧建设生态林带；宝山罗店美兰湖、嘉定远香湖公园及浦东新区森兰公园附近均以建设用地转成水域、草地与林地为主，成为主要城市公园。（3）非显著区转变成低低值聚集区中 ESV 值降值区占此类区域总面积的 71.25%，成片分布于工业园的建设区、大型居住社区新建区和市政大型工程建设区(图4c)，其中工业园区有依托临港发展的临港重装备产业基地临港物流园区（奉贤园区）海港综合开发区、上海化学工业区奉贤分区（基地）杭州湾沿岸、松江工业区东区、嘉定工业园区（北区）；大型居住社区新建区有浦东新区航头镇航三路航海路航浦路沪南公路合围区、闵行区吴泾镇申嘉湖高速樱桃河淡水河图 3 上海市生态系统服务价值 LISA 分布图Fig.3 LISA distribution of ESV in Shanghai（a:非显著区转变成高高值聚集区；b:低低值聚集区转变成非显著区；c:非显著区转变成低低值聚集区；d:高高值聚集区转变成非显著区）图 4 20102020 年生态系统服务价值 LISA 转变分布图Fig.4 Distribution of LISA transformation of ESV from 2010 to 2020上海国土资源 Shanghai Land&Resources84 2023Vol.44.3 剑川路合围区、滴水湖临港大道申港大道地区；市政大型工程建设区主要在浦东机场东侧临海区。工业园的建设区以林地、耕地转化为建设用地为主，大型居住社区新建区多以耕地、水域转化为建设用地为主，浦东机场东侧临海区多以耕地、水域转化为草地、林地为主。（4）高高值聚集区转变成非显著区中的 ESV 降值区域面积仅占上海市 ESV 减少区总面积的 6.73%，而ESV 值减少量占总减少量的 30.99%。主要分布于崇明岛的长江北港北支与东滩地区、浦东新区的南汇沿岸新围垦滩地、青浦区南部、松江区西南部以及奉贤区南部地区。其中，崇明岛的长江北港北支与东滩地区、浦东新区的南汇沿岸地区均以水域转成耕地为主，青浦、松江及奉贤区内多以水域转成耕地、林地为主（图 4d）。4 结论本文基于 20102020 年土地利用变更数据，采用当量因子法和 GIS 空间分析技术，分析了上海市（陆域）2010、2020 年两期的土地利用变化及其 ESV 的转换矩阵、ESV 时空变化的关系。结果表明：（1）20102020年间，上海市土地利用格局变化明显，林地、草地面积增加，而水域、耕地和建设用地面积减少，反映了上海市近 10 年城市化进程背景下，构建“双环、九廊、十区”市域生态空间结构，以保障城市生态安全，提升城市生态环境品质的建设结果。期间上海市 ESV 减少 52.65 亿元，主要是因为水域面积减少所致。（2）在 500 m500 m 格网尺度下，上海市生态系统服务价值在空间上总体表现为高高、低低集聚，局部呈现高低、低高集聚，上海市 ESV 高值区主要集中于崇明岛、淀山湖、杭州湾、南汇边滩等大片生态红线保护区域；ESV 低值区主要集中于市中心城区。（3）崇明岛长江北港北支与东滩地区、浦东新区的南汇沿岸地区的围垦造田和近郊青浦、松江及奉贤区的土地整治以及工业园区、城镇化大型居住社区、市政重大工程的重大城市建设是近 10 年间上海市 ESV 下降的重要因素；近10年间崇明岛（退耕还湿还林工程）、杭州湾（恢复生态用地工程）、淀山湖（湖泊生态治理工程）为代表的重要生态保护区上海市域环城绿带、近郊绿带连接两者的楔形绿地骨干河道生态带（滨水生态空间建设工程）区级公园不同层级生态空间建设的重要成效区，使上海市 ESV 有一定的提升，表明上海市现有生态整治工程和手段效果明显。人类活动是上海市 ESV 变化的主要驱动因素。参考文献(References)1 DAILY G C.The value of nature and the nature of valueJ.Science,2000,289(5478):395-396.2 EGOH B,ROUGET M,REYERS B,et al.Integrating ecosystem services into conservation assessment:A reviewJ.Ecological Economics,2007,63:714-721.3 LAUTENBACH S,KUGEL C,LAUSCH A,et al.Analysis of historic changes in regional ecosystem service provisioning using land use dataJ.Ecological Indicators,2011,11:676-687.4 李雅,刘玉卿.滩涂湿地生态系统服务价值评估研究综述 J.上海国土资源,2017,38(4):86-92.LI Y,LIU Y Q.Review of ecosystem service evaluation of a tidal flat wetlandJ.Shanghai Land&Resources,2017,38(4):80-92.5 黄耀文,张晓平.国家重点生态功能区生态系统服务价值时空变化研究以江西省大余县为例 J.上海国土资源,2022,43(3):99-104,123.HUANG Y W,ZHANG X P.Study on the spatio-temporal changes of ecosystem service value in national key ecological function areas:Take Dayu County of Jiangxi Province as an exampleJ.Shanghai Land&Resources,2022,43(3):99-104,123.6 王双,刘永强.浙江省土地景观格局演变及其生态服务价值效应分析 J.上海国土资源,2021,42(3):63-70.WANG S,LIU Y Q.Analysis of the evolution of land use landscape pattern and value effect of ecological services in Zhejiang provinceJ.Shanghai Land&Resources,2021,42(3):63-70.7 王丹妮.基于上海市土地利用方式变化下的土壤有机碳储量核算研究 J.上海国土资源,2021,42(3):1-6.WANG D N.Soil organic carbon storage based on land use change in ShanghaiJ.Shanghai Land&Resources,2021,42(3):1-6.8 舒波,李雨哲,王玲,等.乐山市土地利用与生态系统服务价值时空动态变化分析 J.生态科学,2022,41(1):159-168.SHU B,LI Y Z,WANG L,et al.Spatial and temporal dynamic changes of land use and ecosystem service value in Leshan CityJ.Ecological Science,2022,41(1):159-168.9 杨光宗,吕凯,李峰.基于格网尺度的南昌市土地利用变化及生态系统服务价值时空相关性分析 J.中国土地科学,2022,36(8):121-130.YANG G Z,LYU K,LI F.Spatial and temporal correlation analysis of land use change and ecosystem service value in Nanchang City based on grid scaleJ.China Land Science,2022,36(8):121-130.10 程江,杨凯,赵军,等.基于生态服务价值的上海土地利用变化影响评价 J.中国环境科学,2009,29(1):95-100.CHENG J,YANG K,ZHAO J,et al.Impact assessment of land use change in center district of Shanghai based on ecosystem services valueJ.China Environmental Science,2009,29(1):95-100.11 鹿亚楠,石登荣,冯艳之,等.上海市土地利用与生态服务价值变化研究 J.北方环境,2011,23(7):198-199.LU Y N,SHI D R,FENG Y Z,et al.Land-use on and the ecosystem service value changes in Shanghai cityJ.Northern Environment,2011,23(7):198-199.12 谢高地,甄霖,鲁春霞,等.一个基于专家知识的生态系统服务价值化方法 J.自然资源学报,2008,23(5):911-919.XIE G D,ZHEN L,LU C X,et al.Expert knowledge based valuation method of ecosystem services in ChinaJ.Journal of Natural Resources,2008,23(5):911-919.13 闫安安,刘雁.大同市土地利用变化对生态系统服务价值的影响 J.农村经济与科技,2022,33(3):49-52.YAN A A,LIU Y.Effects of land use change on ecosystem service value in Datong cityJ.Rural Economy and Science-Technology,2022,33(3):49-52.（下转第 111 页）上海国土资源Shanghai Land&Resources 2023Vol.44.3 11125 肖 云 飞,陈 文 业,王 斌 杰,等.祁 连 山 国 家 级 自 然 保 护区土地利用时空变化及与气候因子关系研究 J.草地学报,2021,29(9):2049-2057.XIAO Y F,CHEN W Y,WANG B J,et al.Study on temporal and spatial change of land use and its relationship with climate factors in oilian mountain national nature reserveJ.Acta Agrestia Sinica,2021,29(9):2049-2057.26 冯姗姗,罗萍嘉.面向空间规划体系变革的资源型城市生态转型规划策略研究以徐州市为例 J.上海城市规划,2020,30(1):55-60.FENG S S,LUO P J.Planning strategies of ecological transformation of resource-based cities oriented to spatial planning system reform:a case study of XuzhouJ.Shanghai Urban Planning Review,2020,30(1):55-60.27 廖谌婳,吴克宁.矿区农田景观生态适宜性评价以徐州市潘安采煤塌陷区为例.资源与产业,2012,14(03):36-42.LIAO S H,WU K N.A case study on Xuzhous Panan coal mine subsidence:assessment of farmland landscape ecological suitability in mining area J.Resources Industries,2012,14(03):36-42.28 董光,何兰,王艺积,等.19902017 年河北省易县土地利用时空格局动态变化 J.中国农业资源与区划,2020,41(1):242-249.DONG G,HE L,WANG Y J,et al.Study on spatio-temporal pattren of land use change in yi county,Hebei Province from 1990 to 2017J.Chinese Journal of Agricultural Resources and Regional Planning,2020,41(1):242-249.29 毕雪昊,凌宇,杨小艳,等.徐州市农村居民点用地时空变化及其影响因素 J.江苏农业科学,2020,48(10):263-268.BI X H,LING Y,YANG X Y,et al.Spatial-temporal changes of rural residential land use in Xuzhou City and its influencing factorsJ.Jiangsu Agricultural Sciences,2020,48(10):263-268.30 崔秀萍,郭永新.近 40a 西部城市群土地利用结构时空动态演化分析 J.干旱区资源与环境,2022,36(2):16-24.CUI X P,GUO Y X.Analysis on the spatio-temporal dynamic evolution of land use structure of western urban agglomerations in the past 40 yearsJ.Journal of Arid Land Resources and Environment,2022,