中国陆地生态系统净初级生产力时空变化特征及影响因素_涂海洋.pdf

第 卷第 期 年 月生态学报 ，：基金项目：第三次新疆综合科学考察项目（）；中国科学院 类战略性先导科技专项（）；青海省科技成果转化专项（）收稿日期：；网络出版日期：通讯作者 ：涂海洋，古丽加帕尔，于涛，李旭，陈柏建中国陆地生态系统净初级生产力时空变化特征及影响因素生态学报，（）：，（）：中国陆地生态系统净初级生产力时空变化特征及影响因素涂海洋，古丽加帕尔，于 涛，李 旭，陈柏建，中国科学院新疆生态与地理研究所，荒漠与绿洲生态国家重点实验室，乌鲁木齐 中国科学院中亚生态与环境研究中心，乌鲁木齐 中国科学院大学，北京 摘要：为揭示气候变化背景下我国各陆地生态系统净初级生产力（）的时空分布特征与驱动机制，引入重心模型分析 年我国 的空间分布格局变化，并利用相关分析方法结合 模型定量区分气候变化与人类活动影响 的相对作用。结果表明：（）年全国 均值为 ，整体呈现出南方高北方低，东南向西北逐渐递减的特点。（）近 年全国与各陆地生态系统 均呈现增长趋势，全国 增长速率为 ，总净增加约 。空间上全国与森林、草地、荒漠生态系统的 重心向东北方向移动，农田与城市生态系统的 重心向西北方向移动，表明 在该方向上的增速和增量最大。（）全国 在华北、西北地区与四川盆地主要受降水的影响，在青藏高原与云贵高原的东部主要受气温的影响，各陆地生态系统之间城市生态系统 对降水响应的敏感度相对最高，荒漠生态系统 对温度响应的敏感度相对最高。（）气候变化和人类活动对全国 变化的贡献分别约占 和，各陆地生态系统的 也表现为受气候变化的影响大于人类活动的影响。此外气候变化对荒漠生态系统 的影响最大，对森林生态系统的影响最小，人类活动对各陆地生态系统 的影响大小则相反。研究结果可为我国生态环境保护、资源的合理配置以及制定环境、社会经济发展等战略决策提供科学依据。关键词：陆地生态系统；净初级生产力；重心模型；模型；驱动机制 ，：（），：（），：，（），（），（），：；植被作为陆地生态系统的主体，对气候变化具有很重要的反馈作用。净初级生产力（，）是指单位时间、单位面积上绿色植物通过光合作用产生的有机质总量减去自身呼吸消耗后所储存下来的部分，即绿色植物所积累的有机干物质总量，它受降水、温度、沉降、浓度等诸多因素共同作用影响。能够用于判定生态系统的健康状况和可持续发展水平，其时空变化能够直接体现区域陆地碳循环的变化，是揭示区域碳源 汇过程的重要指标。近年来，随着遥感观测技术的快速发展以及对生态系统过程机理认识的不断加深，已成为表征生态系统受全球变化影响的重要指标和核心内容。部分学者借助 揭示区域生态环境质量的变化情况，还有一些学者利用 对生态系统服务价值进行评价。但目前 的研究多集中于区域、省份和生态脆弱区，侧重于分析 的时空分布特征和趋势，并且大多没有区分或只是针对单一的生态系统进行分析，对于全国尺度和多个生态系统的定量研究以及 在空间与时间上的变化轨迹研究较少。同时全国尺度上也仅研究 的空间分布格局或探讨 对气候变化与人为活动的响应机制，没有将两者有机结合起来。此外，气候变化与人类活动被认为是影响生态环境质量的主要因素，传统的方法主要从统计学的角度分析生态环境质量变化的原因，没有基于生态过程考虑其变化的机理。是反映生态系统功能的重要参数，目前研究通过比较潜在（）与实际（）之间的差异分析气候变化与人类活动影响 的相对作用，进而区分气候变化与人类活动对生态环境质量的影响。该方法不仅考虑到植被的生理生态过程，并能够在像元尺度体现出气候变化与人类活动对 的影响，利于为下一步的生态系统管理提供指导。本研究通过重心模型，结合趋势分析方法分析全国及各陆地生态系统（农田、森林、草地、城市、荒漠）大生 态 学 报 卷：尺度、长时间序列的 时空变化特征和变化轨迹，利用相关性分析与 模型定量区分 变化的驱动机制。结果有利于全面掌握我国各陆地生态系统的健康状况，为生态治理工程效益的验证提供参考。同时分析在全球气候背景下气候变化与人类活动对我国陆地生态系统 的影响，进而为我国生态环境保护、资源的合理配置以及制定环境、社会经济发展等战略决策提供科学依据。数据源 产品来源于 数据集（：），该数据集已通过第三阶段的验证，空间分辨率为。由于该数据集在研究区域内有填充值，不具有 数值大小的意义，因此本研究将其设置为空值，不作考虑。此外，已有研究论证 产品在我国东南部地区、黄土高原和天山中段等区域的应用是可靠的，本研究也利用前人的实测 数据对该产品进行验证，发现两者之间、，综上说明 产品用于本文研究可行。因此本文基于该 产品进行了后续的分析。中国土地利用现状遥感监测数据来源于“国家青藏高原科学数据中心”（：），空间分辨率为 ，是我国目前运用精度较高的土地利用遥感监测产品，产品缺少南海部分岛屿数据。本研究涉及、年共四期监测数据，并根据一级类的范围划分农田、森林、草地、城市、荒漠五类生态系统。为排除土地利用类型变化对研究结果的影响，后续分析仅考虑 年间土地利用类型没有发生改变的区域。通过与 数据集叠加，去除 空值所在区域的土地利用现状数据，得到剩余所有土地利用类型的面积约为 万，其中土地利用类型没有发生改变的区域面积约为 万，占比约。中国省级行政边界矢量数据来源于中国科学院资源环境科学与数据中心（：）。年降水量与平均气温数据来源于国家科技基础条件平台国家地球系统科学数据中心（：）。时间分辨率为月，空间分辨率 ，数据经 个独立气象观测点数据进行验证，结果可信。为保证数据的一致性，便于后续分析，以上数据均统一为以 为基准面的 等面积圆锥投影，空间分辨率 。研究方法 重心模型“重心”在地理学中被定义为使得区域诸多要素在空间上平衡的点，重心坐标能够清晰的反映出地理现象在空间与时间上的变化轨迹、程度以及趋势。通过研究 重心的移动轨迹，一方面能够反映 的空间分布在时间上发展的不均衡性和偏向性，另一方面，通过进一步探讨发展不均衡的原因，能够辅助区域生态环境保护相关的决策制定。由此引入区域重心的概念。假设区域由 个平面空间单元所组成，为第 个空间单元的属性值，（，）为该单元的重心坐标，那么这一区域 个空间单元的重心坐标的均值就构成了区域的重心坐标（，），公式如下：（）（）期 涂海洋 等：中国陆地生态系统净初级生产力时空变化特征及影响因素：趋势分析使用最小二乘法逐像元计算从起始年到终止年，各年同一时段的 在该时段内随年份变化的一元线性方程的斜率。一般认为斜率为正表示多年来 增加，斜率为负表示 减少，斜率等于 表示 不变。相关性分析由于各陆地生态系统净初级生产力（）往往受到以降水、气温等为主的气候因子的影响。因此本研究将降水、气温作为 的气候影响因素，并运用偏相关与复相关分析方法，逐像元探讨全国与各陆地生态系统 与降水、气温的相关性和显著水平，从而揭示 对气候因子的响应情况。以、为临界值划分 与降水、气温之间的偏相关、复相关系数区间，表示高度线性相关，表示中度线性相关，表示低度线性相关，表示两者不呈线性相关。为区分不同区域、陆地生态系统的 受年降水量、气温因子的影响，在前人研究基础上，定义降水、气温因子影响 的相对作用的四种类型（表）：表 降水、气温因子影响 的相对作用类型 驱动因子 强共同作用 弱共同作用 降水作用 气温作用 ：净初级生产力 ；：与降水、气温之间复相关系数的 显著性检验结果 ；：与降水之间偏相关系数的 显著性检验结果 ；：与气温之间偏相关系数的 显著性检验结果 模型 模型考虑的气候因素较为简单且能够较好反映影响植物生长发育的关键因子，其利用潜在蒸散模型对 模型进行改进，与气候变化作用的潜在（）拟合效果更好，因此被广泛用于计算潜在。目前国内利用 模型估算潜在 的相关工作已开展较多，且均取得了较好的结果，表明其在中国区域适用性较好。其计算公式如下：（）（）其中：（）（）式中，是年均实际蒸散量（）；是年均降水量（）；是年均最大蒸散量（）；是年均温度（）。人为（）是指受人类活动影响产生的，一般用潜在 与实际 的差值来表示。（）年 的变化量由 表示，其计算公式如下：（）式中，是 时间序列的斜率，包括（实际）、（潜在）、（人为），为除基准年外的总年份数。本研究利用、，区分气候变化与人类活动对全国与各陆地生态系统 的影响。其中：表生 态 学 报 卷：明 逐年增加，表明 逐年减少；表示气候变化促进植被生长，表示气候变化导致植被退化；表示人类活动导致植被退化，表示人类活动促进植被生长。根据各 的斜率组合，定义 变化的 种情景，（表）。表 气候变化与人类活动对 增加 减少的相对作用情景 变化情况 情景影响因素的贡献率 气候变化 人类活动 增加情景 （）情景 情景 减少情景 （）情景 情景 ：实际 的斜率 ；：潜在 的斜率 ；：人为 的斜率 ；：年潜在 的变化量 ；：年人为 的变化量 结果与分析 空间变化特征图 年全国平均 空间分布 年全国平均 空间分布特征 年全国 分布具有明显的区域分异规律，整体呈现出南方高北方低，东南向西北逐渐递减的特点（图）。近 年全国平均 为 。小于 的面积最大，占比为，主要分布在我国的西北和青藏高原地区，该地区很大一部分面积常年被沙漠或积雪覆盖，土地贫瘠、植被稀疏，大多低于 。之间的面积占比为，主要分布在我国的东北、华北、华中以及华东地区北部的部分省份。的区域占比，主要分布在华东地区的南部、西南地区的东北部、华南地区以及海南岛的北部。大于 的区域占比为，主要分布在台湾的东部、海南岛的南部、云南省的西南部以及青藏高原的东南部。年 时空分布迁移特征本研究按照时间序列连接 年各年 重心，构建重心演变曲线（图）。整体上看，年全国 重心由湖北省十堰市西南部迁移至陕西省安康市东北部，呈现出向东北方向移动的趋势，共移动 。农田生态系统 重心向西北方向移动 ，森林生态系统 重心向东北方向移动，草地生态系统 重心向东北方向移动 ，城市生态系统 重心向西北方向移动，荒漠生态系统 重心向东北方向移动 。可以看出 年农田与城市生态系统 期 涂海洋 等：中国陆地生态系统净初级生产力时空变化特征及影响因素：在西北方向的增速和增量最大，全国与森林、草地、荒漠生态系统 在东北方向的增速和增量最大。图 全国与各陆地生态系统 空间分布格局变化 图 年全国 变化趋势空间分布 时间变化特征 年全国 增加的区域面积占比达到了，减少的区域仅占，增加的区域在全国分布较广，而减少的区域主要分布在我国南方地区（图）。此外，年全国 表现为净的增加，其净增加量约为 ，平均每年净增加 ，期间 总的增加量为 ，总的减少量为 。各陆地生态系统 增加与减少的区域占比分别为：农田（，）、森林（，）、草地（，）、城市（，）、荒漠（，）。且各陆地生态系统 近 年来均表现为净增加，其中农田生态系统 净增加 （），森林生态系统 净增加 （），草地生态系统 净增加 （），城市生态系统 净增加 （），荒漠生态系统 净增加 （）。由此说明在不改变土地利用类型的情况下全国植被状况具有明显的改善趋势，全国 的增长主要是由农田、森林、草地生态系统贡献的（共占全国 增加量的）。近 年来全国与各陆地生态系统的 增长速率从大到小依次为农田（）、全国（）、城市（）、森林（）、草地（）、荒漠（）。平均 的变化也表明（图），年全国与各陆地生态系统的 整体上生 态 学 报 卷：均呈现上升的趋势，平均 的最低值均出现在 年，除农田与城市生态系统的最高值出现在 年之外，其余最高值均出现在 年。单位面积 均值范围从大到小依次为：森林（）、农田（）、全国（）、城市（）、草地（）、荒漠（），该结果与前人研究结果一致。图 年全国与各陆地生态系统平均 年际变化趋势 与气候因子相关性分析 全国及各陆地生态系统 与降水、气温因子相关性图 结果显示，年全国 与降水呈正相关的区域占比，主要分布于我国东北地区的中部、华北、华中和西北地区以及四川盆地、华南地区的西部和华东地区的北部；呈负相关的区域占比，主要分布于我国西南地区以及华东地区南部、东北地区的北部、南部等地，这些地区由于降水、土壤水分较为充沛或是冷湿地区主要受温度影响等原因，受降水影响较小，而降水增加导致的太阳辐射减少，温度降低等现象，会限制植被的生长。与气温呈正相关的区域占比，主要分布于我国的西北、西南、华南地区以及华东地区的北部、华北地区的东部和东北地区的南部；呈负相关的区域占比，主要分布于我国的华中地区以及东北地区的北部、华东地区的南部。年全国 与降水、气温的复相关系数处于 之间，相关性较高的区域主要分布于华北、西南地区。年各陆地生态系统 与降水 气温呈正相关的面积均大于负相关所占的面积（表），且以低、中度相关为主。与降水偏相关的 显著性检验结果相比，温度呈显著关系的面积占比有所下降，两者分布范围与中 高度正相关系数范围基本一致。各陆地生态系统 与降水、气温的复相关系数均处于 之间，通过、显著性水平检验的范围与相关系数高的区域基本一致。降水、气温因子作用的空间分布图 结果显示，在全国范围内，受降水、气温影响显著的区域可分为两部分。第一个部分包括华北地区、四川盆地以及新疆伊犁地区、塔城地区、天山北坡经济带及塔里木河流域绿洲地带。在气温不变的情况下，随降水的增加而增加，当降水不变时，随气温升高而增加的区域仅沿着内蒙古东南方向与其他省份的交界处以及华北平原的部分区域分布，其它区域气温升高而 基本不变，可知该区域主要受降水条件影响；在青藏高原与云贵高原的东部等地，当气温不变时，降水增加导致 减少，当降水不变时，与气温变化的趋势保持一致，可知该区域主要受气温条件的影响。结合图 与图 可以看出，农田生态系统主要受降水的影响（），分布在东北平原、华北平原的北部和四川盆地以及塔里木盆地边缘；森林生态系统受降水、气温因素影响（），受降水影响的区域主要 期 涂海洋 等：中国陆地生态系统净初级生产力时空变化特征及影响因素：图 与降水、气温因子之间相关性 分布在大兴安岭、太行山周边，受气温影响的区域主要分布在横断山脉与云贵高原的东部；草地生态系统受降水、气温因素影响（），受降水影响的地区主要分布在华北地区的东北部以及新疆的伊犁地区、塔城地区以及天山北坡经济带，受气温影响的区域主要分布在青藏高原；城市生态系统主要受降水的影响（），影响范围主要包括东北平原、华北平原以及内蒙古的中部；荒漠生态系统受降水、气温因素的影响（），受降水影响的范围主要包括内蒙古的中部与东部、天山北麓以及新疆两大盆地的边缘，受气温影响的区域主要分布在三江源地区。受降水、气温影响显著的区域中，除荒漠生态系统外，全国及其余陆地生态系统 受降水影响的面积最大。各陆地生态系统之间城市生态系统 对降水响应的敏感度相对最高，其受降水主导的面积占比达，其次是农田（）、草地（）、森林（）与荒漠（）。荒漠生态系统 对温度响应的敏感度相对最高，其受温度影响的面积占比达，其次是草地（）、森林（）、农生 态 学 报 卷：田（）和城市（）。表 与降水、气温的偏相关、复相关结果面积占比 区域降水偏相关 气温偏相关 复相关 正相关负相关 显著（）显著（）正相关负相关 显著（）显著（）显著（）显著（）全国 农田 森林 草地 城市 荒漠 图 降水、气温因素影响 的相对作用空间分布 期 涂海洋 等：中国陆地生态系统净初级生产力时空变化特征及影响因素：图 降水、气温因子影响 相对作用的面积占比 图 潜在 与人为 斜率的空间分布 气候变化和人类活动对 的影响 年全国潜在 时间序列的斜率处于 之间，人为 时间序列的斜率处于 之间（图）。结合气候变化与人类活动对实际 的贡献率的空间分布结果（图），气候变化引起 增加的区域占比为，增加 ，主要分布在我国黑龙江河流域、东北平原、华北平原、内蒙古高原的东部、长江中下游平原、四川盆地、三江源地区以及祁连山和天山山脉；人类活动引起 增加的区域占比为，增加 ，主要分布在青藏高原、横断山脉以及新疆北部和秦岭华北平原中部一带；气候变化导致 减少的区域占比，减少 ，主要分布在青藏高原东南部；人类活动导致 减少的区域占比为，减少 ，主要分布在我国华东地区南部、华南地区的东部以及云南省的南部。通过计算气候变化和人类活动对 变化的贡献生 态 学 报 卷：图 气候变化与人类活动对 的贡献率空间分布 率得出，气候变化和人类活动对全国 变化的贡献分别约占 和，其中对 增加的贡献分别约为 和，对 减少的贡献分别约为 和，各自促使 净增加约、。在不改变土地利用类型的情况下，年各陆地生态系统气候变化影响 的面积占比均大于人类活动影响的占比（图）。此外，近 年各陆地生态系统 产生的变化中气候变化与人类活动所占的贡献分别为：农田（，）；森林（，）；草地（，）；城市（，）；荒漠（，）。以上结果表明，各陆地生态系统 受气候变化的影响要大于人类活动的影响。结合 的时间变化分析，气候变化与人类活动的影响都引起 的净增长，增长量分别为：农田生态系统增加 （）、（），森林生态系统增加 （）、（），草地生态系统增加 （）、（），城市生态系统增加 （）、（），荒漠生态系统增加 （）、（）。图 气候变化与人类活动影响 的相对作用面积占比 期 涂海洋 等：中国陆地生态系统净初级生产力时空变化特征及影响因素：依据各陆地生态系统中 受气候变化与人类活动的贡献占比，在生态系统之间进行对比。近 年各陆地生态系统中 受气候变化影响最大的是荒漠生态系统，气候变化的贡献达到了，其次是农田（）、草地（）、城市生态系统（），最小的是森林生态系统（），受人类活动影响最大的是森林生态系统，其次是城市、草地、农田生态系统，最小的是荒漠生态系统。讨论 年全国的植被状况整体上具有改善趋势，研究期内 的最大值主要出现在 年与 年，该结果与前人研究结果一致。有一部分原因可以理解为是由气候变化的自然原因所引起的。世界气象组织发布的 年大数据显示，年与 年并列为 年以来的第六最暖年，全球冰川呈消融退缩的趋势，短期内增加了冰川融水径流补给，并在一定程度上对依赖河川径流补给为主的生态系统提供了更多的水源保证。此外，中国气候公报显示，年与 年我国主要粮食作物产区光温水总体匹配较好，气候灾害偏轻，气候条件对植被生长也较为有利。年全国与森林、草地和荒漠生态系统 在东北方向的增速与增量最大。已有研究表明，国家实施的一系列生态恢复工程促进了森林（黄土高原与东北地区）与草地（内蒙古）的恢复生长，同时受干旱区暖湿化影响，北方荒漠植被 呈增加态势，而东南地区由于城市化初期的快速发展，虽近十年 逐步增加，但整体增速相对不及北方。城市与农田生态系统 在西北方向的增速与增量最大。近年来由于耕地占补平衡等原因，西北地区耕地扩张使耕地的聚集度增加，提高了农田灌溉效率与机械化水平，促使耕地质量向好的方向发展。同时在西部地区的自然条件下，我国西部的城市化会在一定程度上持续改善植被覆盖情况。而东部城市化快速发展会减少部分绿色植被覆盖，城市排放的污染物也会对周边耕地产生一些不利影响。全国 变化与气候因子呈显著相关的区域中，华北、西北地区与四川盆地主要受降水影响。其中华北地区与四川盆地主要分布着农田与草地，草地对降水的敏感性较高，作物生长主要取决于区域水热情况和灌溉条件，西北地区由于区域降水较少，使得降水成为该地区植被生长的主控因子。青藏高原与云贵高原的东部主要受气温的影响。因为该区域海拔较高，降水充沛，但低温天气一般意味着日照时长与太阳辐射减少，从而抑制植被的生长，所以温度升高，有利于植被的光合作用，进而促进植被生长。各陆地生态系统之间，城市生态系统 对降水响应的敏感度相对最高，是因为全球气候变化和城市热岛效应加剧了城市热胁迫，温度上升使植被的蒸腾作用增强、土壤水分减少，从而阻碍植被生长，而充足的降雨在满足植被生长需要的同时也可以缓解因高温导致的强蒸散发对植被生长的影响。荒漠生态系统 受降水、气温影响的面积分别占比、，这一结果可能与植被类型有关。荒漠植被中深根灌木能通过其发达的根系从地下水中获取水分补充，相比之下更易受气温的影响；而非深根灌丛的水分补充主要依赖地表水，因此其受水分条件制约。同时，过去 年干旱半干旱区降水以 的速度减少，降水变化幅度很小，因此荒漠生态系统 对降水响应的敏感度相对较小。在不改变土地利用类型的情况下，气候变化对全国与各陆地生态系统 变化的贡献均要大于人类活动，且两者都引起 的净增长。其中气候变化主要表现为区域水热条件的改善延长了植被的生长周期，加速土壤有机质分解与营养元素释放，从而促进植被生长与 的积累。人类活动主要表现为农业机械水平、技术、投入的提升以及高投入生态修复工程的实施等。尤其生态修复工程的实施能够持续增强植被的水源涵养功能，减少土壤侵蚀，改善植被生长环境并促进植被恢复。各陆地生态系统之间 受气候变化影响最大的是荒漠生态系统，最小的是森林生态系统，人类活动对各陆地生态系统的影响则相反。主要原因是荒漠生态系统受人为活动干扰较弱而其他生态系统与人类的生产生活活动息息相关，受人类活动的影响深刻。如京津源风沙治理工程、三北防护林工程、退耕还林还草、退牧还草、围封转移等，同时各生态系统也受到城市化扩张、森林乱砍滥伐、过度放牧等人为活动的不利影响。虽然各陆地生态系统受气候变化与生 态 学 报 卷：人类活动的影响整体上均表现为利于 的增加，但仍有一些区域如青藏高原的东南部与我国东南沿海的部分省份受到人类活动的负面影响，因此在这些区域未来的环境保护、资源配置规划过程中，应充分考虑人类活动如过度放牧、城市扩张占用耕地等产生的不利影响，并制定相关政策指导，以期促使区域的环境状况得到改善。不足与展望本研究利用 产品分析全国与五类陆地生态系统长时间序列的 时空变化，同时定量分析 变化的影响因素，结果具有一定的参考价值，但研究仍存在一些不足。就产品本身而言。由于 产品在我国西北干旱区与青藏高原的部分地区存在数据缺失的情况，因此可能会增加研究结果的不确定性；其次，由于生产该产品的 数据的时间覆盖不完整，也在一定程度上造成产品的估算误差，并对研究结果产生一定影响。就方法而言。虽然 模型对潜在 的拟合效果较好，但未考虑到不同生态系统在相同气象条件下的生产力可能不同；此外，在利用潜在 与实际 的差值区分气候变化与人类活动对 的影响时，对水分利用率等其他相关要素的影响考虑有所欠缺。综上所述，利用长时序产品数据来模拟，并将其作为模型的输入数据估算 可以更全面地反映全国的植被长势情况。同时融合多学科知识，加深了解 的形成机理与影响因素，深化模型的生理生态过程机理研究，这将是下一步研究的重点。结论（）年全国 分布具有明显的区域分异规律，整体呈现出南方高北方低，东南向西北逐渐递减的特点，全国 均值为 ，各陆地生态系统的 均值从大到小依次为：森林、农田、城市、草地和荒漠。（）全国和森林、草地、荒漠生态系统的 重心向东北方向移动，农田与城市生态系统的 重心向西北方向移动，表明 在该方向上的增速和增量最大。（）年全国植被状况整体上具有改善趋势。全国 的平均增长速率为 ，各陆地生态系统的 增长速率从大到小依次为农田、城市、森林、草地和荒漠。近 年来全国 总的净增加量约为 ，全国 的增长主要是由农田、森林、草地生态系统贡献的（共占全国 净增加量的）。（）全国的 在华北、西北地区与四川盆地主要受降水的影响，青藏高原与云贵高原的东部主要受气温的影响。各陆地生态系统之间，城市生态系统 对降水响应的敏感度相对最高，其次是农田、草地、森林与荒漠。荒漠生态系统 对温度响应的敏感度相对最高，其次是草地、森林、农田和城市。（）在不改变土地利用类型的情况下，气候变化和人类活动对全国 变化的贡献分别约占 和，各陆地生态系统 也表现为受气候变化的影响大于人类活动的影响，且两者都引起 的净增长。各陆地生态系统之间气候变化对荒漠生态系统 的影响最大，对森林生态系统的影响最小，人类活动对各陆地生态系统 的影响大小则相反。参考文献（）：孙金珂，牛海鹏，袁鸣 中国陆地植被生态系统 空间格局变迁分析 农业机械学报，（）：董晓宇，姚华荣，戴君虎，朱梦瑶 年内蒙古荒漠草原植被物候变化及对净初级生产力的影响地理科学进展，（）：，（）：期 涂海洋 等：中国陆地生态系统净初级生产力时空变化特征及影响因素：刘洁，孟宝平，葛静，高金龙，殷建鹏，侯蒙京，冯琦胜，梁天刚基于 模型和 数据的甘南草地 时空动态变化研究草业学报，（）：，（）：马良，朱再春，曾辉 评估过程模型应用研究进展中国沙漠，（）：肖文发，徐德应 森林能量利用与产量形成的生理生态基础 北京：中国林业出版社，汲玉河，周广胜，王树东，王丽霞，周梦子 年秦岭地区植被生态质量演变特征及驱动力分析植物生态学报，（）：徐洁，谢高地，肖玉，李娜，江源，陈文辉国家重点生态功能区生态环境质量变化动态分析生态学报，（）：牛丽楠，邵全琴，宁佳，黄海波西部地区生态状况变化及生态系统服务权衡与协同地理学报，（）：吴英迪，蒙吉军中国自然资源生态服务重要性评价与空间格局分析自然资源学报，（）：安妮，宁小莉，海全胜，丁蕾基于 数据的近 年浑善达克沙地植被净初级生产力时空分布研究干旱区资源与环境，（）：尹小君，祝宏辉，高军，郭丽洁，王娇娇气候变化和人类活动对天山北坡净初级生产力变化的影响农业工程学报，（）：刘旻霞，焦骄，潘竟虎，宋佳颖，车应弟，李俐蓉青海省植被净初级生产力（）时空格局变化及其驱动因素生态学报，（）：，：，（）：，：？，：，：，（）：，（）：，：，：，：，（）：崔林丽，杜华强，史军，陈昭，郭巍中国东南部植被 的时空格局变化及其与气候的关系研究地理科学，（）：谢宝妮，秦占飞，王洋，常庆瑞 黄土高原植被净初级生产力时空变化及其影响因素 农业工程学报，（）：乔郭亮，金晓斌，顾铮鸣，杨绪红，徐伟义，尹延兴，周寅康 年天山中段高海拔草地暖季承载力农业工程学报，（）：，（）：安尼瓦尔买买提，杨元合，郭兆迪，方精云新疆天山中段巴音布鲁克高山草地碳含量及其垂直分布植物生态学报，（）：，（）：杨潇，郭兵，韩保民，陈舒婷，杨飞，范业稳，何田莉，刘悦青藏高原 时空演变格局及其驱动机制分析长江流域资源与环境，（）：陈舒婷，郭兵，杨飞，韩保民，范业稳，杨潇，何田莉，刘悦，杨雯娜 年青藏高原植被 时空变化格局及其对气候变化的响应 自然资源学报，（）：生 态 学 报 卷：刘洋洋，章钊颖，同琳静，王倩，周伟，王振乾，李建龙中国草地净初级生产力时空格局及其影响因素生态学杂志，（）：李登科，王钊基于 的中国陆地植被 变化特征分析生态环境学报，（）：刘刚，孙睿，肖志强，崔天翔 年中国植被净初级生产力时空变化及其与气象因素的关系生态学报，（）：苏胜涛，曾源，赵旦，郑朝菊，吴兴华基于 数据的中国陆地植被 估算模型优化与分析 基于中国生态系统研究网络数据生态学报，（）：，：，（）：赵东升，郭彩赟，郭义强，刘磊，高璇，张家诚，朱瑜，张雪梅气候变化对“山水林田湖草”重大生态工程的影响生态学报，（）：王玉涛 年以来人类活动对中国植被净初级生产力的影响 兰州：西北师范大学，欧阳玲，马会瑶，王宗明，王昭伟，于显双气候变化与人类活动对内蒙古东部草地净初级生产力的影响生态学报，（）：刘雪佳，董璐，赵杰，杜自强，张红我国荒漠植被生产力动态及其与水热因子的关系干旱区研究，（）：，：黄海潮，温良友，孔祥斌，陈文广，孙晓兵中国耕地空间格局演化对耕地适宜性的影响及政策启示中国土地科学，（）：，：赵安周，刘宪锋，裴韬，王金杰，张安兵，宋辞京津冀地区不透水表面覆盖率的时空演变及其对植被的影响地理研究，（）：陈印军，肖碧林，方琳娜，马宏岭，杨瑞珍，易小燕，李倩倩中国耕地质量状况分析中国农业科学，（）：，：孙雷刚，王绍强，徐全洪，王苗苗京津冀地区生态系统生产力估算及时空格局分析地理与地理信息科学，（）：卢阳，杨亦辰，吴昕燃，王磊，万梓文，郑哲俊，徐家平，曹畅南京城市街区湿球温度时空变化特征及其影响因子应用生态学报，（）：温旭丁，罗赵慧，符良刚气候与土地利用变化对粤港澳大湾区 的影响广西科学，（）：方霞中国干旱半干旱生态系统生产力和碳动态对气候变化响应的数值模拟研究 乌鲁木齐：新疆大学，李净，刘红兵，李彩云，李龙基于 的近 年中国北部植被生长季始期变化研究地理科学，（）：卢龙辉，陈福军，许月卿，黄安，黄玲京津冀“生态系统服务转型”及其空间格局自然资源学报，（）：刘宪锋，朱秀芳，潘耀忠，李宜展，赵安周 年中国植被覆盖时空变化特征生态学报，（）：王耀，张昌顺，刘春兰，甄霖三北防护林体系建设工程区森林水源涵养格局变化研究生态学报，（）：张彪，王爽，李庆旭，谢高地京津风沙源治理工程区水源涵养功能时空变化分析生态学报，（）：刘文超，刘纪远，匡文慧 陕北地区退耕还林还草工程土壤保护效应的时空特征 地理学报，（）：洪长桥，金晓斌，陈昌春，王慎敏，杨绪红，项晓敏集成遥感数据的陆地净初级生产力估算模型研究综述地理科学进展，（）：期 涂海洋 等：中国陆地生态系统净初级生产力时空变化特征及影响因素