乔雪丽¹ , 郑江华^1,2 , 穆晨³     Zheng_jianghua@126.com.

QIAO Xueli¹ , ZHENG Jianghua^1,2 , MU Chen³     1-2]。草地净初级生产力(NPP)是草地碳固定能力的重要表征, 草地NPP研究通过对草地生态系统生产力的模拟, 定量分析其生产力的时空变化, 正确评价在自然条件下草地的生产能力, 从而对系统分析草地在全球气候变化中的生态价值和贡献、研究陆地碳循环机制等方面的研究提供了依据^[3]。草地NPP的研究对于合理利用草地资源, 最大限度地提高草地产量以及对草地农业生态系统中的碳循环研究都具有重要的指导意义^[4-6]。

自20世纪, 随着遥感技术的迅速发展和基于遥感观测生理生态理论研究的进展, 利用遥感影像反演陆地植被净初级生产力(NPP)成为模拟NPP的主要发展方向, 这使得多尺度生态系统草地NPP的空间格局和过程的定量、动态观测成为可能^[7]。基于遥感数据的CASA(Carnegie-Ames-Stanford Approach)模型建立在植物光合作用过程和Monteith^[8]提出的光能利用率概念上, 因其能够在全球和区域尺度上很好地评估NPP的动态变化和时空变异性而得到广泛的应用^[9], 例如周夏飞等利用CASA模型分析了2001—2015年青藏高原草地碳源/汇时空及其与气候因子关系^[10], 黄露等利用CASA模型研究土地利用/覆被变化对鄂尔多斯市草地生态系统净初级生产力的影响^[11], 韩王亚等则利用CASA模型分析了2000—2015年拉萨河流域NPP时空变化及驱动因子^[12]。

随着遥感技术的迅速发展, 中高空间分辨率遥感数据类型逐渐增多, 为区域NPP估算提供了多源遥感数据。由于受时间分辨率及云雨天气的影响, 难以获取完整时间序列的中高空间分辨率遥感影像。为解决中高空间分辨率卫星数据缺失的问题, 基于多源遥感数据时空融合技术得到发展, 并已在植被、土地覆盖与土地利用中得到较多的应用^[13-17]。近年来, 大多数学者从不同的插值方法探究其对NPP精度的影响^[18-19], 但研究不同数据源对草地NPP精度影响的学者较少。

基于以上背景, 本文将处于亚洲大陆腹地的青海省海北藏族自治州作为研究区域, 探讨多源遥感数据在反演草地NPP时的精度, 并进行精度检验, 为改善与提高该区域生态系统生产力的估算与模拟提供有效的数据参考与支撑。

1 研究区概况

本研究区域海北藏族自治州位于青海省东北部, 北面与甘肃省毗邻。全州土地面积3.4万km², 横跨98°05′—102°41′E, 36°44′—39°05′N如图 1。境内最高海拔5287 m, 最低海拔2180 m, 海拔超过3000 m的高原面积约占全州土地总面积的85%为草原地带。海北藏族自治州位于亚洲大陆腹地, 属高原大陆性气候, 东南季风及西南季风微弱。受高海拔条件的制约, 气温极低, 按气候四季的标准划分, 这里全年皆冬, 无明显四季之分, 仅有冷暖二季之别, 干湿季分明^[20]。全年日照时数在2440—3140 h。年降水量300—500 mm, 且集中在7、8月份, 土壤以高山草甸土和山地草甸土为主, 兼有黑钙土、栗钙土、灰褐土等, 有机质含量丰富, 有利于农作物和牧草的生长。

2 数据来源及预处理 2.1 遥感数据 2.1.1 Landsat 8 OLI数据

研究中使用了Landsat 8 OLI(Operational land imager, 陆地成像仪)影像(http://www.gscloud.cn/), Landsat 8卫星的成功发射(2013年2月在美国加州)使得Landsat对地观测得以继续^[21]。在性能上相比之前的Landsat卫星拥有明显提升, 除了保持前阶段Landsat系列的基本特征外, 还在波段数量、光谱范围及辐射分辨率上进行了改进, 这些优点将会对全球生态变化检测产生积极作用^[22]。本地区选取的5景影像的成像时间都在2017年3月份, 且影像质量较好。

2.1.2 MODIS数据

本研究首先将下载得到的MODIS16天合成产品NDVI影像(https://modis.gsfc.nasa.gov), 一共2景影像。利用专业处理软件MRT TOOLS对数据进行格式转换、投影转换将HDF格式转换成Tiff格式, 之后利用ENVI 5.3完成影像裁剪和拼接, 以备后续工作使用。

2.1.3 天宫二号宽波段成像仪数据

研究中使用了天宫二号对地观测数据的宽波段成像仪数据(http://www.msadc.cn/sy/), 天宫二号空间实验室于2016年9月15日成功发射, 设计在轨寿命2年, 是中国第一个真正意义上的空间实验室, 将用于进一步验证空间交会对接技术及进行一系列空间试验。本研究使用数据为宽波段成像仪的可见光近红外数据的二级产品, 已经过辐射定标和几何校正。研究中所使用的遥感数据介绍如表 1。

2.2 气象数据

模型中需要月均温(0.1℃)、月降水量(0.1 mm)和月均辐射(0.01 MJ/m²)参数, 气象数据从国家气象信息中心获取(http://data.cma.cn)。本研究月均温、月降水量气象站点是基于海北藏族自治州300 km的缓冲区区域选取的33个气象站点。月均辐射数据站点的选取是基于海北藏族自治州500 km的缓冲区区域选取的13个气象站点。气象要素栅格图像的生成过程如下：①根据所获取的站点的经纬度信息, 将其转换成空间矢量数据, 定义空间投影与影像信息保持一致的WGS84坐标系；②进行Kriging插值, 利用掩膜工具裁剪海北藏族自治州的栅格图像。以上处理基于ArcGIS 10.2软件实现。

2.3 土地利用类型数据

本研究采用的海北藏族自治州植被覆盖分类方法已发表在《地理学报》英文版上^[3], 数据来源于中国科学院资源环境科学数据中心(http://www.resdc.cn)的2015年中国土地利用现状遥感监测数据库数据集。将所获取的土地利用类型图像在ArcGIS中利用栅格计算器提取高覆盖度草地、中覆盖度草地和低覆盖度草地, 并利用全国基础数据库(2017融合版)裁剪青海省海北藏族自治州地区的草地利用类型图(图 2)。

3 方法 3.1 Gram-Schmidt图像融合

基于格拉姆-施密特(Gram-Schmidt)变换的融合方法是一种多光谱数据融合方法^[23], 其融合效果较好, 可使空间信息明显增强, 光谱保真度较高^[24]。图像融合方法如图 3。

本研究为验证天宫二号宽波段成像仪数据在估算草原NPP方面的精度, 将高分辨影像Landsat 8 OLI影像与经过几何校正的天宫二号宽波段成像仪影像进行Gram-Schmidt图像融合。

3.2 NPP估算模型

已被全球1900多个实测站点校准的CASA模型, 是由遥感、气象、植被以及土壤类型数据共同驱动的光能利用率模型。采用朱文泉等改进的CASA模型计算海北藏族自治州的草地NPP^[25], CASA遥感估算模型为：

(1)

式中, P_PAR为光合有效辐射, 单位为MJ/m²；F_FPAR为植被层对入射光合有效辐射(PAR)的吸收分量, CASA模型使用归一化差异植被指数(Normalized Difference Vegetation Index, NDVI)代替F_FPAR进行估算, 无单位；Ts为温度对光能利用率的影响, 无单位；Ws为水分对光能利用率的影响, 无单位；E_max为最大光能利用率, 单位为gC/MJ；SR代表比值植被指数, PAR是根据SR和NDVI估算得到的。根据查阅文献可知3月该地区的最大光能利用率为0.115、0.085和0.07^[3]。CASA模型中静态参数如表 2。

3.3 检验方法

为检验不同数据源数据对草原NPP的预测精度, 在本研究中将高空间分辨率影像Landsat 8 OLI反演产生的每种NPP类型百分比作为参照物。将天宫二号宽波段成像仪数据、影像融合数据和MODIS数据分别与Landsat 8 OLI数据进行计算均方根误差(Root mean squared error, RMSE)、平均绝对误差(Mean absolute error, MAE)、平均绝对百分比误差(Mean absolute percentage error, MAPE)和相关系数(Correlation coefficient)R_xy, 将其作为检验不同数据源精度的标准^[19]。MAE由于将不同数据源草地NPP值与Landsat草地NPP值的差值绝对值化, 不会出现正负相抵消的情况, 因此能更好地反映偏离的实际情况；RMSE表示不同类型数据源草地NPP值偏离程度, 可反映离散程度；MAPE描述了偏离Landsat草地NPP值的相对程度^[26]；采用逐像元的皮尔逊相关系数分析对海北藏族自治州的不同数据源数据的NPP值与Landsat 8 OLI的草地NPP之间的相关性进行定量研究, 并进行F检验, 当P＜0.05时, 相关性显著, R_xy揭示不同遥感源的草地NPP值与Landsat 8 OLI草地NPP值之间的相关性。

(2)

(3)

(4)

(5)

式中, n为分类数；Z⁰(x_i)为第i种类型NPP值的百分比；Z(x_i)为Landsat 8 OLI影像反演结果的第i种类型NPP值百分比；x和y分别表示Landsat 8 OLI与天宫二号、融合影像和MODIS影像反演的草地NPP样本值的平均值。

技术路线如图 4。

4 结果与分析 4.1 不同数据源数据反演NPP值

利用不同数据源的影像反演的青海省海北藏族自治州的3月份草地NPP值存在显著差异。从表 3可以看出：Landsat 8 OLI的NPP值位于150—200 g C m^-2 a^-1所占比例最高为30.14%, 所对应比例最接近的数据类型为MODIS数据为18.55%；Landsat 8 OLI反演NPP值大于200 g C m^-2 a^-1的比例次之为25.71%, 所对应比例最接近的数据类型为天宫二号宽波段成像仪影像数据类型；Landsat 8 OLI反演NPP值100—150 g C m^-2 a^-1的比例为23.49%, 所对应比例最接近的数据类型为天宫二号宽波段成像仪影像数据；Landsat 8 OLI反演NPP值50—100 g C m^-2 a^-1的比例为15.90%, 所对应比例最接近的数据类型为MODIS数据；Landsat 8 OLI反演NPP值小于等于50 g C m^-2 a^-1的比例最小为4.76%, 所对应比例最接近的数据类型为天宫二号宽波段成像仪影像数据类型。从图 5中可以初步得到, 天宫二号宽波段成像仪影像数据反演的NPP值最为接近Landsat 8 OLI数据反演的NPP值, 反演NPP精度排序为天宫二号宽波段成像仪影像>融合影像>MODIS。

利用CASA模型反演的NPP值如图 5所示。

4.2 不同数据源数据反演NPP值检验

统计不同数据源数据的NPP值分布比例(图 6), 将得到的NPP像元百分比值代入公式(2)、(3)、(4)和(5)中, 计算其MAE、RMSE、MAPE和R_xy, 结果如表 4：

(1) 通过观察MAE值, 发现MAE(天宫二号) < MAE(融合影像) < MAE(MODIS), 此处的MAE值代表 3种数据源数据NPP值与Lnadsat 8 OLI草地NPP值偏离的实际情况。可知MODIS与Landsat 8 OLI草地NPP值的MAE最大, 反映了MODIS数据反演NPP值精度较差；天宫二号宽波段成像仪影像数据的草地NPP值与Landsat 8 OLI数据的NPP值最为接近, 精度要高于另外两种数据。

(2) 比较表 4中的RMSE值, 发现RMSE(天宫二号) < RMSE(融合影像) < RMSE(MODIS), 此处的RMSE值代表 3种数据源数据NPP值与Lnadsat 8 OLI草地NPP值偏离的实际情况。可知天宫二号宽波段成像仪影像数据相较于融合图像和MODIS影像与Landsat 8 OLI的草地NPP值偏离较小, 精度最高, MODIS数据精度效果最差；

(3) 观察MAPE值, 发现MAPE(天宫二号) < MAPE(融合影像) < MAPE(MODIS), 此处的MAPE值代表 3种数据源数据NPP值与Lnadsat 8 OLI草地NPP值偏离的相对程度。天宫二号宽波段成像仪影像数据反演草地NPP值精度最高。

(4) R_xy结果值为R_xy(天宫二号)>R_xy(融合影像)>R_xy(MODIS), 此处的R_xy值代表 3种数据源数据NPP值与Lnadsat 8 OLI草地NPP值的相关性。同样天宫二号宽波段成像仪影像数据反演草地NPP值与Lnadsat 8 OLI草地NPP值的相关度最高, 且相关性显著。

5 结论与讨论

本文基于CASA模型估算了青海省海北藏族自治州的草地植被净初级生产力, 同时从不同数据源数据出发, 分析了其对应的反演结果精度, 初步得到以下结论：

(1) 青海省海北藏族自治州2017年3月的草地NPP值受不同数据源的影响而呈现不同的空间分布。Landsat 8 OLI数据的NPP值位于150 g—200 g C m^-2 a^-1所占比例最高；天宫二号宽波段成像仪影像数据和融合后影像的NPP值位于50 g—100 g C m^-2 a^-1所占比例最高；MODIS数据的NPP值位于小于50 g C m^-2 a^-1比例最高；

(2) 天宫二号宽波段成像仪影像数据较融合后影像和MODIS数据更能反映海北藏族自治州的草地净初级生产力的时空分布, MODIS精度要低于天宫二号和融合影像反演的NPP值的精度；

(3) 天宫二号作为2016年9月15日成功发射的空间实验室, 适用于反演青海省海北藏族自治州的草地NPP。

CASA模型是光能利用模型, 其估算NPP值是基于植被光能利用效率和光合有效辐射, 估算结果能够表征未干扰天然草地的发育状况和演变过程, 但该模型参数中并未考虑人类活动, 如放牧强度和人工恢复等对于草地NPP的影响, 因此, 在CASA模型的改进和NPP的定量化估算中, 如何增加人为因素影响将是今后研究的方向之一^[27]。

如今, 研究者主要利用MODIS影像产品数据、Landsat影像数据或GIMMS和SPOT VGT两种数据集的NDVI数据来反演植被的净初级生产力, 由于所研究对象空间尺度较大, 那么如何获取既具备高空间分辨率和高时间分辨率的影像仍是一个方法性的难点, 这也是今后的研究工作之一。

近年来, 越来越多的学者利用遥感数据参与到气候变化的研究中, 本团队也进行了相关研究^[28], 未来的工作中应继续深入探讨不同数据源与气候因子在NPP值之间的相关性和不同数据源对草地NPP的估算精度。基于不同数据源利用CASA模型对草地NPP估算的影响因素较多, 应扩大时间尺度、变换插值方法、探究草地类型对NPP的影响。研究中发现天宫二号宽波段成像仪影像数据反演的草地NPP精度高于利用Gram-Schmidt算法融合的影像数据, 误差原因的可能性有多种, 例如：融合算法、模型参数或评价方法等, 应加大对此方向的研究, 从而对草业检测提供参考。

致谢: 感谢载人航天工程提供天宫二号宽波段成像仪数据产品。