不同格网尺度下生态系统服务价值空间分异及其影响因素差异——以大南昌都市圈为例.pdf

第 卷第 期 年 月生态学报 .:/.基金项目:国家自然科学基金资助项目()收稿日期:网络出版日期:通讯作者 .:.:./.危小建辛思怡张颖艺龙英豪张茜.不同格网尺度下生态系统服务价值空间分异及其影响因素差异 以大南昌都市圈为例.生态学报():.:.():.不同格网尺度下生态系统服务价值空间分异及其影响因素差异 以大南昌都市圈为例危小建辛思怡张颖艺龙英豪张 茜 东华理工大学测绘工程学院 南昌 自然资源部环鄱阳湖区域矿山环境监测与治理重点实验室 南昌摘要:在江西省持续推进生态文明建设的背景下生态系统服务价值()空间分异及其影响因素的探析对生态环境保护与改善、促进区域性可持续发展等方面具有重要意义 以大南昌都市圈为例利用当量因子法、空间自相关分析等方法分析 、及 各个格网尺度下地均生态系统服务价值的空间分布特征并利用地理探测器和空间回归模型研究不同尺度下 空间异质性的影响因素及其尺度差异性 研究结果表明:()分布总体呈现西北部高东南部低的特点且各种格网尺度下均存在显著空间正相关性和空间集聚效应但随着格网尺度增大其集聚效应减弱()空间异质性受自然、社会的协同作用其中人为影响指数的贡献最大且任意双因子都比单一因子对 空间异质性的解释力高但随着格网尺度增大各因子及因子间的耦合协调作用对 的解释力都呈下降趋势()随着格网尺度的增大空间回归模型的拟合度下降且不同格网尺度下影响 空间异质性的影响因素的作用强度不同作用方向也有发生变化关键词:多尺度生态系统服务价值影响因素空间回归地理探测器大南昌都市圈 :().:():/.().().:生态系统服务是指人类能够从生态系统获取的所有价值包括可以保障其赖以生存和发展的环境资源及其他效用 随着经济社会发展、城市化水平提升、民众生活水平提升等带动生态系统服务消费增长生态系统承受着巨大的压力 为提高生态系统服务消费供给水平与经济社会发展的协调水平急需对现有的生态系统进行合理调整 在此背景下大量学者开始关注生态系统服务价值的研究并进一步分析和探讨其演变趋势及影响因素美国学者 等于 年全面系统地介绍了生态系统服务价值()的定义后受到了国际学界的普遍重视引发了研究 的兴趣浪潮 目前 的核算方法一般包括二大类:一是由 等建议的当量因子法该方式核算具有规范性结果具有横向水平一致性但不能完全反映某地区内的具体生态特点较适用于评价地区及世界范围的 二是功能价值法该方法能够专门核算某种土地利用类型的生态价值并应用于国内各大科学研究领域 谢高地等于 年研究中国生态系统及经济社会发展现状并在 等观点的基础上对当量因子法加以修正和完善为我国学者进一步深入研究 提供了借鉴 现如今国内学者主要以行政区划或格网为基本单元以当量因子法、空间自相关模型、地理探测器模型为研究方法测算研究区域的、分析 的时空演变及探讨其空间异质性的影响因素 例如黄木易等以 为基本单元、张鹏岩等以 的格网为基本单元核算了研究区域的 并利用地理探测器分析影响 空间异质性的影响因素赵育恒等以城市为单元估算 年武汉城市圈的 并采用空间回归模型探究其影响因素 在研究 的影响因素过程中空间尺度非常重要在不同尺度中生态系统格局呈现不同特征 尺度较大时结果不够精确容易忽略单元内各因素的影响尺度较小时所得结果难以推广至大范围生态调节应用之中 因此选择合适的尺度有利于对 进行研究 黄木易等以多尺度格网、乡镇行政区划和子流域为基本单元分析不同尺度下 的空间分布及其异质性并利用地理探测器探究其影响因素 相关研究已取得较为显著的成果但仍存在一定不足有关 影响因素的探究中大多以单一尺度为主没有关注到在不同格网尺度下同一地区的 影响因素可能存在差异 此外同时利用地理探测器和空间回归模型进行影响因素对 的正负影响及两两因子间的交互作用力的研究的相关文章并不多见本文选取大南昌都市圈为研究区域利用 年土地覆盖数据采用 的 工具得到生 态 学 报 卷:/.、及 等多种空间尺度格网结合 和 软件采用当量因子法、空间自相关分析等方法分析地均 的空间分布特征进而运用空间回归分析和地理探测器两种模型研究 空间异质性的主导因子和各因子交互作用强度及其不同尺度下影响因子的差异性 本文旨在揭示大南昌都市圈内各影响因素对生态系统服务价值空间异质性的作用为改善生态环境质量、促进绿色发展提供科学依据 研究区概况及数据来源.研究区概况大南昌都市圈(图)地处中国江西省北部区域涵盖了南昌市、九江市、抚州市(临川区、东乡区)、宜春市(丰城市、樟树市、高安市、靖安县、奉新县)、上饶市(鄱阳县、余干县、万年县)以及国家级新区赣江新区等 个县市国土面积.万 由于位于我国城镇化的战略布局长江横轴线与京九地区发展轴线的交汇处大南昌都市圈在我国地区发展布局中有着沟通南北、承东启西的关键战略地位 同时该区域北临长江西依幕阜山和九岭山东含鄱阳湖和庐山包括了山区、丘陵、平原和河湖等各类山势地貌西边为山地、中部为山丘东边则为平原和湖泊 其中鄱阳湖作为中国面积最大的淡水湖它承载了调控径流、净化水体、调控气候等多重生态功能有着大量的鱼类、鸟类等生物资源在维护地球生物多样性上有着难以取代的作用 该区域内绿色生态显著是推动中国国家生态文明实验区(江西)建设的重点地区图 大南昌都市圈土地利用分布图.数据来源与处理本文使用了谷歌地球引擎上 年的 影像结合 提取的稳定样本制作的土地利用数据空间分辨率为 数据来源于:/./根据中科院“遥感监测数据分类体系”综合考虑大南昌都市圈的土地资源特征将研究区域的 数据分为耕地、林地、草地、水域、建 期 危小建 等:不同格网尺度下生态系统服务价值空间分异及其影响因素差异:/.设用地和未利用地 类 粮食产量数据来源于 年的江西统计年鉴和中国统计年鉴粮食价格根据 年中国农产品价格调查年鉴获得本文选用的影响 空间异质性的指标包括:高程、归一化植被指数()、气温、降水、地均、人口密度数据来源如表 所示坡度利用高程数据在 中计算所得城镇用地从 数据提取人为影响指数()利用公式计算获得 运用 的自然间断点分级法对各因子进行离散化处理表 数据来源与说明 数据名称 数据类型 时间 分辨率数据来源 网址高程格网 美国对地观测全球高程数据:/./归一化植被指数 格网 中国科学院地理科学与资源研究所资源环境数据云平台:/./气温格网 国家地球系统科学数据中心 年中国 分辨率逐月平均气温数据集:/./降水格网 国家地球系统科学数据中心 年中国 分辨率逐月降水量数据集:/./地均 格网 中国科学院地理科学与资源研究所资源环境数据云平台:/./人口密度 格网 国家地球系统科学数据中心中国 格网人口数据集:/./在查阅相关文献后发现采用格网尺度做此类研究选用的最小尺度为 的格网 此外人口密度和 等影响因素的分辨率为 考虑到数据处理的易操作性以及结果差异性决定研究的尺度选取 、及 格网 研究方法.生态系统服务价值测算以“中国生态系统单位面积生态服务价值当量”为基础计算大南昌都市圈生态系统服务价值按单位面积耕地生态系统服务价值相当于平均粮食作物单产市场经济价值的七分之一计算通过查阅江西统计年鉴、中国统计年鉴、和国家统计局关于 年粮食产量的公告等资料经测算后确定大南昌都市圈粮食作物单位面积产量为./该地区主要农产品为水稻、大豆和油菜籽选取 年水稻、大豆、油菜籽单价的平均值作为实际计算的粮食单价通过查阅中国农产品价格调查年鉴计算得出大南昌都市圈粮食单价为.元/参照前人研究将建设用地的价值系数假定为 进而确定大南昌都市圈的单位面积生态系统服务价值系数(表)根据公式()测算研究区内生态系统服务价值总量 ()()式中 为生态系统服务价值元为第 种生态系统类型面积为第 种生态系统服务价值系数元/生 态 学 报 卷:/.表 大南昌都市圈生态系统服务价值系数表 生态系统类型 生态系统服务价值 耕地林地草地水域建设用地 未利用地 食物生产 .原材料 .气体调节 .气候调节 .水文调节 .废物处理 .保持土壤 .生物多样性.美学景观 .合计.利用敏感性指数解析生态系统服务价值相对生态系统价值系数的依赖性程度(式)()/()/()式中 为敏感性指数、分别表示初始值和当量系数调整后的值.生态系统服务价值的人为影响指数人为影响指数()反映人类的社会活动对生态系统及自然景观结构变迁的影响程度利用 指数能够反映大南昌都市圈的 与人类行为影响间的相互关系其公式如下:(/)()式中 为人为影响指数为第 种土地利用类型的总面积为第 种土地利用类型所反映的人为影响强度参数 为研究区总面积 根据前人研究成果并综合大南昌都市圈土地利用现状特征运用 法确定各土地利用类型的 赋值:耕地(.)、林地(.)、草地(.)、水域(.)、建设用地(.)、未利用地(.).探索性空间数据分析.全局空间自相关分析选用全局 统计量对不同尺度下的地均 进行全局空间自相关检验公式如下:()()()()式中 为网格单元总数()为网格单元()的测度值()为第 个网格单元上测度值与平均值的偏差 为标准化的空间权重矩阵 为方差下同.局部空间自相关分析通过局部空间自相关可分析大南昌都市圈内任一区域与其相邻区域对该属性值的集聚趋势空间异质性常见方法有 统计量、空间关联局部指标()、散点图等 本文利用 集聚图对不同尺度下的地均 进行局部空间自相关检测公式如下:()()()期 危小建 等:不同格网尺度下生态系统服务价值空间分异及其影响因素差异:/.()().地理探测器地理探测器()是通过其表现处理的空间分异性来探测各因子对因变量的驱动力包括四个探测器本研究通过因子探测器和交互探测器判断影响大南昌都市圈 空间异质性的主导因子及其交互作用以及在不同尺度下各因子影响程度差异公式如下:()式中 为某因子对 空间异质性的影响力取值为 为变量 或因子 的分层 和 分别为全区和 区单元数s和s分别为全区和 区 值的方差.空间回归模型地理探测器仅仅研究因子对 的影响程度而无法确定其正负相互影响为了更全面的研究各种因素对 的驱动影响揭示各因子和 之间的正负相关性构建了空间回归模型其基本形式如下:()()式中()分别代表被解释变量(一个 维向量)、解释变量()、误差项()与随机扰动项()的空间权重矩阵 则表示空间回归系数通过逐步线性回归来剔除共线性的变量然后构建 邻接的空间权重矩阵利用空间计量回归的方法进一步研究 与各因素之间的内在驱动关系 结果与分析.生态系统服务价值的空间变化.生态系统服务价值的空间布局由图 可以看出在不同格网尺度下地均 的空间分布较为一致总体呈现为西部和北部较高、东部和南部较低的特点 高价值区主要集聚于鄱阳湖水域及其周围的湿地较高价值区主要位于西部的幕阜山和九岭山等林地中等价值区散乱分布于林地内部低价值区和较低价值区则主要分布于各县市城区、乡村居民点用地和耕地通过敏感性分析发现研究区内各地类 的 值均低于 表明大南昌都市圈 相对于价值系数而言是比较稳定的因此本文所用的价值系数适用于大南昌都市圈.生态系统服务价值的空间自相关分析由表 可知不同格网尺度下地均 分布均呈现显著的空间集聚效应和空间正相关性即地均 低值区与高值区都趋于邻近 随着格网尺度的增大 值也呈现了逐步降低的特点表明研究区内地均 高值和低值的集聚效应随着格网尺度的增大在逐步降低 格网尺度较小时大多数格网中地类结构相似其地均 相近集聚效应明显随着格网尺度的增大各格网中的地类结构将发生不同的变化其地均 也趋于不同集聚效应减弱 因此生态服务价值空间差异对评价单元尺度的变化具有一定的敏感性由图 可以发现在不同格网尺度下地均 的空间集聚分布都比较一致高高集聚区都主要分布于鄱阳湖水域以及武宁县、修水县、奉新县、靖安县之间的山区因为此区域内包含大量的生态资源是大南昌都市圈 的重点来源区域高低和低高集聚区分布相对较少且散乱且大多处在高高和低低集聚地附近低低集聚区则大多散布于研究区的中部及南部该区域多为建设用地及耕地地均 均较低生 态 学 报 卷:/.图 不同格网尺度下地均 的空间分布.表 不同格网尺度下地均 全局 统计量 格网尺度 ()()格网 .格网 .格网 .格网 .()为理论期望()为标准差 为概率即所观测到的空间模式是由某一随机过程创建而成的概率.生态系统服务价值空间分异的地理探测分析.因子探测从表 可以看出大南昌都市圈 空间异质性受到自然因素与社会经济因素的共同作用且不同格网尺度下各因子对 的影响程度差异较大 通过 统计量发现 的贡献力最强其中在 格网的贡献力达到了.表明人类活动对生态环境造成了极大的干扰在 、格网下 的 统计量逐步减小这也表明格网尺度越小所反映的空间分异特性就越强 格网的坡度贡献力为 说明地形起伏情况对 的空间分异也有着极大的影响且坡度的 统计值也是随着格网的增大而逐步减小的贡献力在 格网中最高为.表明其在 格网中空间结构异质性最大该因子对研究区 空间异质性的影响很强高程和气温在各格网尺度的贡献力为.对研究区 空间异质性的影 期 危小建 等:不同格网尺度下生态系统服务价值空间分异及其影响因素差异:/.图 不同格网尺度下地均 的 集聚图.响较强这两种因子在 和 格网中的 统计量都较大 和 中则相对较小人口密度在 格网中贡献度最大为.格网次之地均、降水和距城市距离的各尺度贡献力均在.及以下对研究区 空间异质性的影响不大表 不同格网尺度下 空间分异的因子探测 影响因素 格网 格网 格网 格网 统计量 统计量 统计量 统计量 高程.坡度 .归一化植被指数 .气温.降水.距城市距离 .地均 .人口密度 .人为影响指数 .生 态 学 报 卷:/.因子交互探测由图 可知因子交互作用对 的影响大于单一因子的影响 在 格网中 和坡度之间的交互作用对研究区域 的空间异质性影响最大 值为.在 、和 格网 和降水之间的交互作用对研究区域 的空间异质性影响最大 值依次为.、.和.双因子对 空间异质性的解释力超过 的在 格网中有:(.)、高程(.)、气温 期 危小建 等:不同格网尺度下生态系统服务价值空间分异及其影响因素差异:/.图 不同格网尺度下 空间分异的因子交互探测.(.)、降水(.)、地均(.)、距城市距离(.)、人口密度(.)在 和 格网中 和其他因素之间的相互作用对 的空间异质性也有超过 的解释力 其余双因子在各尺度上的交互作用对 的空间异质性的解释力不足 但其交互作用强于单一因子.生态系统服务价值的空间回归分析利用 软件对不同格网尺度 与各影响因素分别进行逐步线性回归选取与研究区 拟合最好的模型构建线性回归方程 根据表 可知 格网回归结果最好格网尺度越大模型拟合结果越差且不同格网尺度下影响 空间异质性的影响因素的作用强度不同作用方向也可能发生变化 具体来看在各个格网尺度下对 解释力最强的均为 呈现显著负作用在 格网中高程、气温、降水、距城市距离、地均、人口密度对研究区 呈正向作用坡度、与研究区 呈负向作用 格网中的人表 不同格网尺度下 逐步线性回归结果 变量 经典最小二乘 格网 格网 格网 格网 常数项 .高程.坡度 .归一化植被指数 .气温.降水.距城市距离 .地均 .人口密度 .人为影响指数 .、依次表示通过、的显著性水平检验生 态 学 报 卷:/.口密度因素、格网中的降水因素、和 以及 格网中的距城市距离因素与 呈负向作用作用方向与其在 格网中不同其他因素对 的作用方向和 格网一致从空间回归结果来看(表)模型的拟合效果较线性回归有所改进总体上仍表现为格网尺度越大模型的拟合效果就越差且在各格网尺度中 对 的解释仍最强呈负向作用 与线性回归结果比较考虑了空间效应在各个格网尺度下各影响因素的作用强度甚至作用方向都存在不同 具体来看与线性回归相比 格网中的高程因素和 格网中的距城市距离因素对 的作用明显发生了变化分别呈负向和正向且均未通过显著性检验 格网中的距城市距离因素和 格网中的降水因素对 的作用均呈正向 格网中的坡度因素和人口密度因素、格网中的人口密度因素、格网中的气温因素和地均 因素均未通过显著性检验其他因素与线性回归结果基本一致表 不同格网尺度下 空间误差回归 变量 空间误差模型 格网 格网 格网 格网 常数项 .高程.坡度 .归一化植被指数 .气温.降水.距城市距离 .地均 .人口密度 .人为影响指数 .系数.、依次表示通过、的显著性水平检验 结论与讨论.讨论对比分析各尺度地理探测器和空间回归结果可以发现空间回归模型与地理探测器揭示出 与各影响因素之间的不同关系 比如 格网中高程在地理探测器因子探测中对 的影响贡献力排序为第四但逐步线性回归分析却建议将其剔除距城市距离在地理探测器因子探测中对 影响力排序为第八而逐步线性回归分析却建议将其输入 经深入分析发现在 尺度上坡度与高程呈现高度相关性(.)而地理探测器中的因子探测并没有考虑多重共线性问题 此外相比于地理探测器空间回归模型能揭示影响因素对 的不同作用方向 因此将地理探测器与空间回归模型结合用于生态系统服务价值空间分异的探测结果较好的展示了其空间分异的主导因子、因子交互特点以及正负作用情况揭示了生态系统服务价值空间分异的影响机制 受多种因素共同作用的影响 其中人为影响指数是 空间分布的主导因子其原因是人类活动如城市扩张和开垦农田等的增多导致土地利用结构发生变化从而影响了 的空间分布 高程、坡度、温度、降水等自然因素对 的空间分布也有着较大的影响其中水热条件直接影响着植被的生长由气候、地质活动、生物等共同作用形成的高程、坡度等地形条件的变化对人类的生产生活和土地利用有着重要影响而 能够反应植被的分布状态及生长潜力是人类改造利用林地的重要依据之一 人口密度、等社会经济因素对 的空间分布影响较小因为研究区内林地、水域所占面积较大境内生态系统稳定性 期 危小建 等:不同格网尺度下生态系统服务价值空间分异及其影响因素差异:/.强受社会经济影响小 各影响因素之间的相互耦合作用加大了对研究区域 空间异质性的影响 与影响因素的关系呈现出明显的尺度依赖特征 随着格网尺度增大格网内各地类的比例变小且分布越分散格网内的同质性变低 的空间相关性也变低 空间自相关低或不显著的区域主要集中在不同地类的边界处格网尺度越大空间自相关低或不显著的区域越大 不同的影响因素也受到了格网尺度变化的影响随着格网尺度增加数据聚合程度越高呈现出更强的空间异质性 随着格网尺度的增大 与影响因素同时受到尺度的影响空间异质性增强使得 与影响因素的拟合程度变低 且由于聚合过程的复杂性不同格网尺度下影响因素是否进入回归方程以及作用方向都可能不同 空间异质性的影响因素差异较大 格网尺度越小影响因素间交互作用越强且双因子交互作用排序也发生改变 因此改善生态系统的过程中需要选取合适的分析尺度且不仅需要考虑不同格网尺度下影响因素的不同还需考虑各影响因素间的交互作用改善生态环境应精准调控生态环境空间积极落实生态保护政策持续推进生态环境治理 在规划改善生态系统服务空间结构时可以考虑其空间布局以及影响程度较大的利益群体 本文研究结果表明人类活动对生态系统服务具有明显抑制作用表明城市扩张、开垦农田等人类活动造成较高价值生态系统服务向较低价值生态系统服务转移导致生态系统向负向发展 因此建议应协调城乡发展合理用地实现高质量的经济发展建设避免城市无序扩张 此外管理政策和措施需要考虑尺度变化的影响选择合适的尺度影响因素 尺度太小其结果难以推广至大范围生态调节应用之中 基于生态环境保护的易操作性及区域性可持续发展的可行性结合研究区的实际概况考虑到数据量的大小本文认为 格网的研究结果更具参考价值 此外需要指出的是格网尺度的分析有助于发现适合的分析尺度但是不如行政区划单元样便于生态环境保护管控和治理 且本文所考虑的影响因素可能不够全面在今后的研究中还需进一步加强与完善.结论以大南昌都市圈为研究对象结合地理探测器和空间回归模型从多尺度角度全面分析大南昌都市圈 的影响因素并探究 空间异质性的主导因子和各因子交互作用以及不同尺度下影响因素作用的差异性结果表明:()在不同格网尺度下地均 分布总体上较为一致均为西部和北部的地均 较高东部和南部的地均 较低 不同格网尺度下地均 呈现一定的空间异质性特点且各格网尺度下均存在显著空间正相关性和空间集聚效应 随着格网尺度的增大地均 的高值和低值集聚效应均在逐渐弱化即生态价值空间集聚效应逐渐减弱()地理探测器结果表明大南昌都市圈 空间异质性受到自然环境和经济社会发展的协同作用不同因子对 的影响程度差异显著在不同格网尺度下 贡献度均为最大是影响大南昌都市圈空间异质性的主导因子且任意两两因子都较单一因子对 空间异质性的解释力高 随着格网尺度的增大因子探测与交互探测统计量明显下降即因子及因子间对 空间异质性的解释力逐渐下降()空间回归的结果显示随着格网尺度增大减小空间回归模型的拟合度降低 格网尺度不同各因子对 空间异质性的影响程度和正负作用方向存在差异 在不同格网尺度下 均对 解释力最强且呈明显负作用说明人类活动对 有较强抑制性参考文献():.:.():.:.:():.():.王思博 龚晶 陈俊红 李冬冬 李婷伟.生态系统服务消费与经济增长、环保投入相关性研究 基于居民旅游对生态景观的需求视角.价格理论与实践():.生 态 学 报 卷:/.():.严恩萍 林辉 王广兴 夏朝宗.年三峡库区生态系统服务价值演变及驱动力.生态学报 ():.():.高登.大都市区土地利用变化及其生态系统服务价值响应 以济南市为例.中国资源综合利用 ():.:.:.:.():.赵育恒 曾晨.武汉城市圈生态服务价值时空演变分析及影响因素.生态学报 ():.彭保发 陈端吕.常德市土地覆被的生态服务价值空间变异分析.经济地理 ():.:.():.():.薛明皋 邢路 王晓艳.中国土地生态系统服务当量因子空间修正及价值评估.中国土地科学 ():.孙新章 周海林 谢高地.中国农田生态系统的服务功能及其经济价值.中国人口资源与环境 ():.王景升 李文华 任青山 刘某承.西藏森林生态系统服务价值.自然资源学报 ():.:.():.韩增林 赵文祯 闫晓露 钟敬秋 孟琦琦.基于生态系统服务价值损益的生态安全格局演变分析 以辽宁沿海瓦房店市为例.生态学报 ():.江波 欧阳志云 苗鸿 郑华 白杨 庄长伟 方瑜.海河流域湿地生态系统服务功能价值评价.生态学报 ():.谢高地 张彩霞 张雷明 陈文辉 李士美.基于单位面积价值当量因子的生态系统服务价值化方法改进.自然资源学报 ():.谢高地 鲁春霞 冷允法 郑度 李双成.青藏高原生态资产的价值评估.自然资源学报 ():.陈睿 杨灿 杨艳 董祥泽.洞庭湖生态经济区生态系统服务价值的时空演绎及其驱动因素.应用生态学报 ():.危小建 陈竹安 张蕾 江平 吴芳.引入城市扩张干扰效应的生态服务价值化方法改进.长江流域资源与环境 ():.叶妍君 齐清文 姜莉莉 张岸.基于地理探测器的黑龙江垦区农场粮食产量影响因素分析.地理研究 ():.王劲峰 徐成东.地理探测器:原理与展望.地理学报 ():.:.黄木易 方斌 岳文泽 冯少茹.近 来巢湖流域生态服务价值空间分异机制的地理探测.地理研究 ():.张鹏岩 耿文亮 杨丹 李颜颜 张宇 秦明周.黄河下游地区土地利用和生态系统服务价值的时空演变.农业工程学报 ():.黄木易 岳文泽 方斌 冯少茹.年大别山区生态服务价值尺度响应特征及地理探测机制.地理学报 ():.():.乔斌 祝存兄 曹晓云 肖建设 扎西拉旦 颜玉倩 陈国茜 史飞飞.格网尺度下青海玛多县土地利用及生态系统服务价值空间自相关分析.应用生态学报 ():.侯孟阳 姚顺波 邓元杰 丁振民 鲁亚楠 郑雪 李雅男.格网尺度下延安市生态服务价值时空演变格局与分异特征 基于退耕还林工程的实施背景.自然资源学报 ():.谢高地 甄霖 鲁春霞 肖玉 陈操.一个基于专家知识的生态系统服务价值化方法.自然资源学报 ():.姚小薇 曾杰 李旺君.武汉城市圈城镇化与土地生态系统服务价值空间相关特征.农业工程学报 ():.刘桂林 张落成 张倩.长三角地区土地利用时空变化对生态系统服务价值的影响.生态学报 ():.刘祗坤 吴全 苏根成.土地利用类型变化与生态系统服务价值分析 以赤峰市农牧交错带为例.中国农业资源与区划 ():.陈浮 葛小平 陈刚 彭补拙.城市边缘区景观变化与人为影响的空间分异研究.地理科学 ():.侯孟阳 姚顺波 邓元杰 丁振民 鲁亚楠 郑雪 李雅男.格网尺度下延安市生态服务价值时空演变格局与分异特征 基于退耕还林工程的实施背景.自然资源学报 ():.雷金睿 陈宗铸 吴庭天 李苑菱 杨琦 陈小花.海南岛东北部土地利用与生态系统服务价值空间自相关格局分析.生态学报():.马晓熠 裴韬.基于探索性空间数据分析方法的北京市区域经济差异.地理科学进展 ():.彭珏 何金廖.电商粉丝经济的地理格局及其影响因子探析 以抖音直播带货主播为例.地理科学进展 ():.期 危小建 等:不同格网尺度下生态系统服务价值空间分异及其影响因素差异