基于遥感反演的1982-2015年中国北方温带和青藏高原高寒草地地上生物量空间数据集

焦翠翠，于贵瑞，陈智，何念鹏

1. 四川轻化工大学经济学院，四川宜宾 644000

2. 中国科学院地理科学与资源研究所生态系统网络观测与模拟重点实验室，北京100101

3. 中国科学院大学资源与环境学院，北京 100190

数据库（集）基本信息简介

数据库（集）组成数据集包括3 8个文件，其中3 4个文件为中国北方温带和青藏高原高寒草地1 9 8 2-2 0 1 5年逐年的A G B空间数据。命名规则为N G T P_Y Y Y Y_A G B.t i f，其中N G T P代表中国北方温带和青藏高原高寒草地，Y Y Y Y代表年份，A G B代表地上生物量；另外4个文件为N G T P_1 9 8 0 s A G B.t i f，N G T P_1 9 9 0 s A G B.t i f，N G T P_2 0 0 0 s A G B.t i f和N G T P_2 0 1 0 s A G B.t i f， 分别代表1 9 8 0 s（1 9 8 2-1 9 8 9年），1 9 9 0 s（1 9 9 0-1 9 9 9年），2 0 0 0 s（2 0 0 0-2 0 0 9年）和2 0 1 0 s（2 0 1 0-2 0 1 5年）A G B的平均值，所有数据单位均为g·m-2。

引 言

草地生物量是指某一时刻单位面积内草地中实际存有的有机物质总量。其中，地上生物量（Aboveground biomass，AGB）不仅能够反映植被的生长状况、表征草地载畜能力[1-2]，还与生态系统质量密切相关[3]。另外，AGB是草地生态系统碳库的重要组成部分[4]，虽然它在草地碳储量中所占比重不大，但它却与草地生态系统的物质循环和能量流动过程紧密相连[5]。中国草地的总面积约为3.9×106km2，居世界第二位。北方温带草地和青藏高原高寒草地面积约占中国草地总面积的79%，是全国重要的畜牧业基地，又是北方和京津地区的重要生态屏障[6]。因此，评估中国这两个区域草地AGB并揭示其动态变化，对我国草地资源的保护、畜牧业的发展、水土保持和生态系统的可持续发展具有重要意义[1-2]。

在草地AGB的众多测定方法中，样地收获法被认为是简单、最准确的一种方法，但该方法耗时、费力，通常只适用于样点及样带尺度的研究[1,7-8]。近年的大量研究表明，地面生物量调查数据与卫星遥感技术相结合，是评估区域尺度草地AGB的有效途径[8-9]。随着过去40年间中国社会与经济的快速发展，北方温带草地和青藏高原高寒草地所承受的扰动和压力也在不断变化。AGB的长期变化不仅能够很好地反映这种变化，而且能够为后续草地的可持续管理提供可借鉴经验。然而，过去对中国这两个区域草地AGB的评估多关注于某一时期的平均状态[10]，较少关注AGB动态变化的研究，探讨的时间尺度也多集中在10-20年间[11]。关于1982-2015年期间中国北方温带和青藏高原高寒草地的AGB数据，到目前为止尚未公开共享的空间数据集，一定程度上阻碍了相关研究的进展。

有鉴于此，本研究以大量AGB实测数据为基础，借助美国国家航空航天局（National Aeronautics and Space Administration，NASA）全球监测与模型研究组（Global Inventory Monitoring And Modeling Studies，GIMMS）发布的长时间序列归一化植被指数（Normalized Difference Vegetation Index，NDVI）数据，科学地构建估算AGB的遥感模型，对1982-2015年中国北方温带和青藏高原高寒草地的AGB进行了评估，在此基础上分别对1980s（1982-1989年），1990s（1990-1999年），2000s（2000-2009年）和2010s（2010-2015年）的平均AGB进行了计算，生产了一套1982-2015年中国北方温带和青藏高原高寒草地AGB的空间数据集。期望通过公开共享模式，为中国草地生态系统（草地资源、畜载力、承载力和碳储量等）动态变化的后续研究提供科学数据。

1 数据采集和处理方法

1.1 基础数据

1.1.1 研究区域概况

中国北方温带和青藏高原高寒草地空间分布范围广，南北跨越23个纬度，东西跨越50个经度。综合过去对北方草地的相关探讨，本研究主要考了内蒙古自治区、青海省、西藏自治区、新疆维吾尔族自治区、甘肃省和宁夏回族自治区6个北方省级行政区范围内的草地，没有包括辽东湾北端向西和陕西一带等范围[11-12]。北方温带草地主要分布在大、小兴安岭向西、西南直至新疆西部国境线的区域。青藏高寒草地主要包括西藏、青海大部分区域、甘肃和新疆的部分地区[12]。根据1:100万中国植被图[13]，我国这两个区域的草地主要包括6种草地类型：高寒草原、高寒草甸、山地草甸、草甸草原、典型草原和荒漠草原（图1）。

本研究将内蒙古自治区、青海省、西藏自治区、新疆维吾尔族自治区、甘肃省和宁夏回族自治区6个省级行政区的分布范围与中国1:100万植被类型数据[13]中草地类型的空间分布范围相重叠的部分作为中国北方温带和青藏高原高寒草地的总体分布范围。将青藏高原的空间分布范围（数据来源详见1.1.4部分介绍）与中国1∶100万植被类型数据中草地类型的空间分布范围相重叠的部分作为青藏高原高寒草地的分布范围；将中国北方温带和青藏高原高寒草地的总体分布范围中除去青藏高原高寒草地之外的部分作为北方温带草地的分布范围。

图1 中国北方温带和青藏高原高寒草地的植被类型和生物量实测数据样点的空间分布

1.1.2 生物量实测数据

本研究中使用的地上生物量实测数据（AGBobs）是采用经典的样地收获法测定得到，主要有以下 2个数据来源：① 通过在中国知网（http://epub.cnki.net）和 Web of Science（www.Webofknowledge.com）以“生物量”“碳储量”“生产力”及“产量”为关键词收集到已公开发表的学术论文230篇，得到1982-2015年期间787个调查点，953个调查点年的AGBobs数据；②相关研究者提供了608个调查点，630个调查点年的AGBobs数据。通过以上2个途径共得到1982-2015年期间中国北方温带和青藏高原高寒草地区域内1395个调查点，1583个调查点年的AGBobs，数据的详细介绍见Jiao[14-15]、Ma[10,16]和Xu[17]等的相关研究。

采用以下3个标准对初步得到的AGBobs数据进行了筛选：① 剔除地理位置信息记录不完整的样点数据；② 剔除缺失采样时间的样点数据；③ 以草地类型为基础，将平均值±3倍标准差作为标准，剔除数值过低或者过高的样点数据。筛选之后最终得到了一个包括 1104个调查点，1259个调查点年的AGBobs数据集（图1）。

1.1.3 归一化植被指数（NDVI）

本研究中使用了1982-2015年的NDVI数据，该数据来源于GIMMS生产的NDVI第三代数据产品NDVI3g.v1[18]。GIMMS的NDVI3g.v1[19]是以搭载在NOAA（National Oceanic and Atmospheric Administration）卫星上的AVHHR（Advanced Very High Resolution Radiometer）传感器采集的卫星数据为基础，通过运用最大值（Maximum Value Composite，MVC）合成方法生产得到的空间分辨率为0.083°×0.083°（～8 km×8 km）[18,20-21]，时间分辨率为 15 d 的 NDVI产品[18]。

1.1.4 空间矢量数据

植被类型分布图采用的是中国科学院中国植被图编辑委员会绘制的1:100万植被类型图[13]，从中国科学院资源环境科学数据中心下载得到[22]。青藏高原空间范围数据也是从中国科学院资源环境科学数据中心下载得到[22]。

1.2 生物量空间数据集生产方法

1.2.1 数据预处理

在构建AGB的遥感估算模型之前，首先通过以下3个步骤对数据进行预处理：

① 以时间分辨率为15 d的NDVI数据为基础，运用MVC方法[23]（公式1）计算得到时间分辨率为1个月的NDVI值。

② 在月NDVI值的基础上计算参与构建AGB遥感估算模型的年NDVI值。通过构建NDVI数据与AGB实测数据之间的经验关系来模拟区域尺度的草地AGB时，通常选用年最大NDVI值或生长季的平均NDVI值。草地生态系统的生长季通常开始于每年的4月、结束于10月，但是表征AGB变化的最优NDVI合成时相会因研究区域内部的异质性，如气候、植被条件等的差异有所不同[24-25]，可能是4-10月、5-10月、4-9月、6-8月、7-8月等不同时期的NDVI值。

在本研究中为了系统地选择表征AGB变化的最优NDVI合成时相，年NDVI的计算时期可能开始于4月、5月、6月或者7月，结束时期可能是8月、9月或者10月，然后以1个月为时间步长进行组合，共计得到12种不同时相组合的生长季平均NDVI值：4-10月、5-10月、6-10月、7-10月；4-9月、5-9月、6-9月；4-8月、5-8月、6-8月；7-8月、7-9月。再加上年NDVI最大值，这样共计得到每年13种不同时相组合的年NDVI值（图2，附表1）。

③ 根据AGBobs的采样年份和地理位置信息分别提取每个样点对应年份由不同时相组合得到的13种年NDVI值，得到AGBobs与对应年NDVI数据集。采用随机抽样的方法选取约占总数75%的样点数据为测试样本，用于模型的构建，约25%的样点数据为验证样本，用于模型的验证。在此需要说明，在年最大NDVI值小于0.1的区域，因植被覆盖率低，NDVI数据受下垫面影响比较大，在本研究中没有考虑。

1.2.2 AGB遥感估算模型的构建

基于预处理后的AGBobs和NDVI数据，分别拟合AGBobs与对应13种年NDVI值之间的回归关系（线性、指数、乘幂及对数），共计得到52种回归模型（附表1）。本研究选取决定系数（R2，公式2）和均方根误差（RMSE，公式3）作为模型精度的评价指标。

公式（2）和（3）式中，R2为模型估算值AGBmod和地面实测值AGBobs之间的决定系数，能够指示AGBmod和AGBobs两组数据之间分布的相似性，表征模型能够解释实测AGBobs变异的程度；RMSE为模型估算值AGBmod和地面实测值AGBobs之间的均方根误差，表征AGBmod与AGBobs之间的偏差程度；n为样点个数。

在构建的52个回归模型中（附表1），模型精度评价结果同时符合RMSE最小和R2最大这两个条件的模型，被认为是AGB的最优估算模型[26-27]。如果模型的精度评价结果表明同时满足RMSE最小和最大这两个条件的模型不存在，那么需要优化现有的模型[26]。在本研究中，将分别满足RMSE最小的模型和最大的模型在像元尺度的估算值取平均值作为AGB的最优估算结果。同时，将符合RMSE最小和最大的模型在像元尺度取平均的结果作为AGB的最优估算模型。

将上述的建模思路分别应用于北方温带草地和青藏高原高寒草地，然后用预留的验证数据分别对两个草地区域各自构建的 52个回归模型进行精度评价。评价结果表明（附表1），对于北方温带草地和青藏高原高寒草地来说，都不存在任何一个模型同时满足最大和RMSE最小这两个条件。根据上述模型构建方法的介绍，在北方温带草地和青藏高原高寒草地的AGB最优估算模型都是由各自最大和RMSE最小的模型（附表 1中加粗显示）在像元尺度取均值的平均模型。AGB遥感估算模型的构建思路参考Jiao等[14]的研究，主要的技术流程如图2所示。

图2 AGB空间数据集生产的技术路线图

根据上述技术流程，构建中国北方温带草地AGB的遥感估算模型AGB-RSMNG（公式4）和青藏高原高寒草地AGB的遥感估算模型AGB-RSMTP（公式5）如下：

2 数据样本描述

经过一系列处理，获取了1982-2015年期间中国北方温带和青藏高原高寒草地逐年的AGB空间数据，并且以此为基础分别计算了1980s（1982-1989年）、1990s（1990-1999年）、2000s（2000-2009年）和2010s（2010-2015年）不同年代的AGB，共同构成了1982-2015年间中国北方温带和青藏高原高寒草地AGB的空间数据集，空间分辨率为8 km，单位为g·m-2，格式为ARCGIS TIFF。以1980s、1990s、2000s和2010s的数据为例，展示AGB的空间数据情况（图3）。

3 数据质量控制和评估

本数据集的生产从基础数据的获取与预处理、遥感估算模型的构建与验证都有严格的质量控制和评估标准，保证了数据的可靠性。AGB实测数据是通过样地收获法得到的，数据获取途径包括文献调研和野外调查两个方面。对于文献调研的数据，检索关键词的确定、数据获取与处理的方法均经过专家的论证和认可。在数据收集完毕之后，由不同人员对数据进行交叉校对，对数据进行单位统一、异常值剔除等后续处理。对于野外调查数据，科研人员在样方设置和样品处理都严格遵循了相应的调查技术规范。

图3 1980s、1990s、2000s和2010s中国北方温带和青藏高原高寒草地地上生物量的空间数据

AGB遥感估算模型所需的NDVI数据来源于GIMMS发布的NDVI3g.v1，该数据集在生产的过程中已经过一系列的预处理：如辐射校正、去云、传感器退化性订正及几何纠正等处理。GIMMS NDVI3g.v1数据集中还提供了NDVI数据的质量控制文件，目前已经在全球及区域尺度生态系统生物量及生产力的空间分异和动态变化研究中被公认和广泛应用[28-30]。

在运用NDVI数据估算AGB的研究中，多数是通过选用预先设置好的某一个固定时期（如每年的5-8月、5-9月、4-9月、4-10月等）的NDVI值与实测AGB数据来构建AGB的遥感估算模型[31-33]。本数据集在生产过程中充分考虑了区域内植被及气候条件差异对表征AGB变化的最优NDVI时相的影响作用，将每年4-10月的月NDVI值按照不同的时间窗口进行组合，得到每年13种不同时相组合的年NDVI值。在此基础上，较系统地筛选出表征北方温带草地和青藏高原高寒草地AGB变化的最优NDVI时相，从一定程度上优化了这类传统模型的计算方法。

为了更加直观地表达本数据集的精度，我们用之前预留的约25%样点的AGB实测数据作为验证数据，与AGB-RSMNG和AGB-RSMTP的模拟结果进行了对比（图4）。对比结果表明，AGB-RSMNG和AGB-RSMTP的模拟结果能够很好地表征中国北方温带和青藏高原高寒草地AGB的变异。对于北方温带草地来说（图4a），AGB-RSMNG的模拟结果与AGB实测数据之间的R2和RMSE分别为0.63和55.38 g·m-2。对于青藏高原高寒草地来说（图4b），R2和RMSE则为0.64和56.43 g·m-2。

图4 由AGB-RSM 的估算结果与AGB实测数据对比验证的结果

4 数据使用方法和建议

基于遥感反演的1982-2015年中国北方温带和青藏高原高寒草地地上生物量的空间数据集，可以通过CERN综合中心数据资源服务网站（http://www.cnern.org.cn）下载得到。用户登录系统后，在首页点击“数据论文数据”图标或在“数据资源”栏目选择“数据论文数据”中的“碳氮水通量观测专题”，进入相应页面后可以下载完整数据。本数据集的生产是在ArcGIS软件平台下完成的，数据格式为ArcGIS TIFF格式，用户需要借助ArcGIS软件才能打开。用户若需要其他格式，可以在ArcGIS软件中进行格式转换。用户使用时需要注意数据的单位，本数据集的单位统一为g·m-2。

在此需要说明的是，本数据集可能存在以下几个方面的不确定性：（1）因本数据集的时间尺度比较长（1982-2015年），为了充分运用已有的调查数据，在构建AGB-RSM时采用的AGB实测数据来源于多个调查者、多次野外采样获得的数据集。尽管这些野外数据在获取时都严格遵守了野外观测与技术规范，但是多次的野外调查规范没有统一，可能会对AGB实测数据的可比性产生一定的影响，也是本文公开数据不确定性的重要来源之一。（2）因考虑到AGB实测数据的可得性，本研究在构建AGB-RSM时，只考虑了因气候不同而形成的北方温带草地和青藏高原高寒草地两大区域之间的差异。对于研究区域的6种草地类型，高寒草原和高寒草甸主要分布在青藏高原高寒草地区域，运用的是同一个AGB估算模型（AGB-RSMTP）；其他4种草地类型主要分布在北方温带草地区域（AGB-RSMNG），运用的是同一个AGB估算模型，尚未针对各种草地类型的差异角度来考虑对模型的影响，这也可能会对AGB-RSM模型的构建造成一定的不确定性。