摘要：本数据集提供了2000年至2024年全球逐年净生态系统生产力（NEP）数据，空间分辨率为1000米。数据集基于多源遥感数据、气象再分析数据和生态系统过程模型生成，反映了在标准气候与环境驱动场条件下陆地生态系统的碳源汇时空动态。数据经过了一致性处理和误差分析等质量控制，确保了其在长时间序列上的可比性与可靠性。该数据集对于评估全球及区域碳收支、研究气候变化对生态系统的影响以及支持相关环境政策制定具有重要的利用价值。

关键词：净生态系统生产力；NEP；CASE模型；30米分辨率

引  言

陆地生态系统作为全球碳循环的关键组成部分，其碳源汇格局的动态变化深刻影响着大气二氧化碳浓度的增长速率，进而对气候变化产生重要反馈。精确量化陆地生态系统碳收支的时空变异，是当前全球变化生态学和地球系统科学研究的核心问题之一，对于理解和预测气候变化、制定科学合理的碳管理政策具有重要意义。

传统的站点观测虽然能够提供精确的生态系统通量数据，但在空间代表性上存在局限，难以直接扩展到区域和全球尺度。随着遥感技术的发展，利用卫星遥感数据驱动生态系统模型或基于光能利用率原理的反演模型，已成为估算大尺度、长时间序列生态系统生产力（如总初级生产力GPP和净初级生产力NPP）的主流方法。净生态系统生产力（NEP）作为GPP扣除自养呼吸和异养呼吸后的剩余部分，直接表征了生态系统与大气之间的净碳交换量，是判断区域为碳源或碳汇的关键指标。

为满足对长时间序列、高空间分辨率全球NEP数据的需求，我们在整合多源遥感数据（如NDVI、EVI等）、气象再分析数据以及生态系统过程模型的基础上，开发了本数据集。它延续并拓展了先前在区域碳循环模拟和基于遥感的生产力估算等方面的研究工作。本数据集提供了2000-2024年间逐年、30米分辨率的全球NEP数据，能够清晰地揭示近20年来全球碳源汇的时空分异特征与演变趋势。

本数据集可为全球变化研究、碳循环模型验证、生态脆弱性评估以及国家尺度碳汇核算提供关键的数据支撑，具有广泛的重用价值。

1  数据采集和处理方法

1.1数据采集方法

本数据集计算使用了归一化植被指数（NDVI），土地利用，太阳辐射,降水，温度。分为预处理和计算两个部分，预处理部分主要通过裁剪，投影的方式统一数据的范围，行列数，投影，分辨率。计算部分通过公式分步计算。

1.1.1归一化植被指数数据（NDVI）数据

2000-2024年数据来源主要依靠GEE（Google Earth Engine）平台，2000-2011年份使用LANDSAT/LT05/C02/T1_TOA数据集。2012-2013年使用LANDSAT/LE07/C02/T1_TOA数据集。由于LANDSAT7影像的扫描行校正器（Scan Line Corrector, SLC）的永久性故障。因此，自2003年5月31日之后获取的LANDSAT 7影像被称为 “SLC-off”数据，在SLC-off影像中，每一景图像都会出现平行的、交错的数据缺失条带。图像中间部分几乎完整，但越往边缘，条带缺失越严重，导致边缘部分约有22%的数据完全丢失。所以2012年和2013年的数据采用了核函数的方法进行了修补。2014-2018年使用LANDSAT/LC08/C02/T1_TOA数据集。2019-2024年使用LANDSAT/LC08/C02/T1_TOA和LANDSAT/LC09/C02/T1_TOA数据集。

1.1.2土地利用数据

本数据来源于中国科学院资源环境科学数据中心 (http://www.resdc.cn)的中国土地利用/土地覆盖数据数据。

中国土地利用/土地覆盖数据数据的分类系统采用三级分类系统（表 1）：一级分为 6 类，主要根据土地资源及其利用属性，分为耕地、林地、草地、水域、建 设用地和未利用土地；二级主要根据土地资源的自然属性，分为 25 个类型；三 级类型 8 个，主要根据耕地的地貌部位，具体分类如下：耕地分水田和旱地（二级类型），水田根据其所处的地貌位置又分为四个三级类型：山地水田（111）、 丘陵水田（112）、平原水田（113）和大于 25°坡地水田（114），旱地根据其所处的地貌位置又分为四个三级类型：山地旱地（121）、丘陵旱地（122）、平原旱地 （123）、大于 25°坡地旱地（124）。

1.1.3 太阳辐射数据

1.1.4 降水数据

1.1.5 温度数据

1.2 数据处理

1.2.1 土地利用数据处理

土地利用数据为每五年一幅的TIFF影像，分辨率为30米，首先使用ARCGIS 10.6软件将其投影，处理1985-2020年间共8幅影像。由于不同土地利用类型的最大光能利用效率（ε_max）不同，目前的利用类型不能与其完全呼应，所以需要将土地利用影像进行重分类。

1.2.2 归一化植被指（NDVI）数据处理

对下载好2000-2024年的逐月NDVI最大值合成数据使用ARCGIS 10.6软件进行投影，设置分辨率30米，并使用裁剪工具，以土地利用数据为基础进行行列数匹配，或者行列数、投影、分辨率均相同的影像。

1.2.3 太阳辐射数据处理

通过 ArcGis软件插值生成的 1985年以来各月度 1000 米分辨率地表太阳辐射栅格数据，单位为MJ。再进行与NDVI同样的投影裁剪操作，获得与土地利用行列数，投影，分辨率相同的影像。

1.2.4 降水数据处理

基于 Anusplin 软件插值生成的 1985年以来各月度 1000 米分辨率降水栅格数据，插值过程中引入高程作为协变量,单位为mm。使用ARCGIS 10.6软件进行与NDVI同样的投影裁剪操作，获得与土地利用行列数，投影，分辨率相同的影像。

1.2.5 温度数据处理

基于 Anusplin 软件插值生成的 1985年以来各月度 1000 米分辨率气温栅格数据，插值过程中引入高程作为协变量，单位为℃。使用ARCGIS 10.6软件进行与NDVI同样的投影裁剪操作，获得与土地利用行列数，投影，分辨率相同的影像。

1.2.6 NPP计算

CASA（Cainegie-Ames-Stanford Approach）模型计算公式如下：

NPP(x,t)=APAR(x,t)×ε(x,t)

其中，NPP(x,t)表示像元x在t月的净初级生产力（gC/m²）；APAR(x,t)表示像元x在t月吸收的光合有效辐射（单位：gC·m-2·month -1）；ε(x,t)表示单个像元x在t月的实际光能利用率(单位：gC·MJ -1 )。

APAR(x,t)= SOL(x,t) ×0.5 ×FPAR(x,t)

式中:SOL(x,t)是时间为t月象元x处的太阳总辐射量(单位：MJ۰m-2۰month -1 ), SOL(x,t)应用经验公式进行计算。0.5指的是植被所能利用的太阳有效辐射(波长是0.4~0.7μm)占太阳总辐射的比例；FPAR(x, t)为植被吸收光合有效辐射的吸收比例；FPAR与归一化植被指数（NDVI）、比值植被指数（SR）皆存在一定的线性关系。

植被吸收的光和有效辐射取决于太阳总辐射和植被本身特性,太阳总辐射可由大气上界太阳辐射量和日照百分率计算。

式中：SOL为陆表短波辐射，常被称为陆表太阳辐射（MJ·m^-2·d^-1）；n为实际日照时数（单位：h），由气象资料提供；N为最大日照时数（单位：h）；n/N为日照百分率，S₀为大气外界辐射量（单位：MJ·m^-2·d^-1）；a_n，b_n表示晴天最大日照时数时到达地面的大气外界辐射分量，参考已有研究选取a_n=0.207；b_n=0.725。

式中：dr为大气外界相对日地距离；Q₀为太阳常数（取0.0820MJ·m^-2·min^-1）；d为赤纬（单位：rad），W_s为太阳时角（单位：rad）；j为纬度（单位:rad）。d，dr，W_s由下式计算：

式中：J为该年中所处的天数。

植被层对光合有效辐射的吸收比例（FPAR)主要受到地表植被覆盖类型和覆盖程度影响。其计算公式如下:

FPAR（x,t）=αFPAR_NDVI+（1-α）FPAR_SR

其中α = 0.5，