关键词:
滨海盐沼湿地
涡度协方差技术
总初级生产力
光学卫星遥感
机器学习方法
摘要:
总初级生产力(Gross primary productivity,GPP)是衡量生态系统功能的重要参数之一,准确估算滨海盐沼湿地GPP对阐明滨海湿地碳收支情况及碳汇能力至关重要。然而,由于滨海湿地具有面积小、破碎化等特征并且受潮汐淹水的影响,区域尺度滨海湿地GPP估算依然具有较大的不确定性,缺乏适用的高分辨率滨海盐沼湿地GPP估算方法。本研究以长江口典型盐沼湿地九段沙湿地与崇明西沙湿地为研究区域,基于线性回归模型和机器学习随机森林回归模型,通过整合涡度协方差(Eddy covariance,EC)、陆地卫星(Landsat)和中分辨率成像光谱仪(Moderate-resolution imaging spectroradiometer,MODIS)影像等多源数据,开展长江口盐沼湿地的高空间分辨率GPP估算模型的构建与评价、影响因素探究以及空间格局的分析。本研究基于不同分辨率的卫星遥感数据源,使用卫星衍生的光合有效辐射(Photosynthetically active radiation,PAR)、增强植被指数(Enhanced vegetation index,EVI)、归一化差异植被指数(Normalized difference vegetation index,NDVI)和归一化差异水指数(Normalized difference water index,NDWI)单独或不同的指标组合方式为驱动因素来实现两种回归模型,并研究它们在估算GPP方面的表现。主要研究结果如下:(1)本研究通过采用完全由遥感数据驱动的VI×PAR模型结构,实现了长江口典型盐沼湿地GPP估算。结果表明,采用线性回归模型进行EVI与PAR的整合通常具有最大的GPP估算潜力,其表现达到了R2>0.75,RMSE<6.80μmol m-2 s-1。在几乎不受潮汐和植物凋落物影响的滨海盐沼湿地中,NDVI的表现超过了EVI(R2>0.78)。与线性回归方法相比,随机森林方法在使用EVI、NDVI、NDWI和PAR组合作为输入变量时,在估算近瞬时尺度的GPP变化时具有最佳性能。此外,使用30米分辨率的Landsat数据比使用500米分辨率的MODIS数据或10米分辨率的Sentinel-2数据更适合于产生准确的GPP估算结果。(2)探讨环境因素特别是潮汐淹水对使用的VI×PAR估算形式估算GPP的影响,结果表明,考虑潮汐淹没的影响后,在剔除NDWI>0.2的数据的情况下,VI×PAR与实测GPP之间的简单线性关系强度普遍有一定程度的提高,EVI×PAR的R2可以达到0.81。并且在随机森林回归模型中,使用EVI、NDVI、NDWI和PAR的组合表现最好,具有最高的测试数据的R2(0.71)和最低的RMSE(6.64μmol m-2 s-1)值。这表明使用机器学习算法,以及增加反映潮汐影响的遥感指数,对于滨海盐沼的GPP估算是有利的。(3)本研究对以长江口典型盐沼湿地为研究区域进行了植被分类和区域尺度的GPP估算。估算结果表明,2005年至2020年期间,长江口盐沼湿地的植被面积整体呈现上升趋势,整体上长江口盐沼植被面积从2005年的121.28km2到2020年的101.07 km2,净损失20.21 km2。不同种类的长江口盐沼植被的面积变化量差异很大,其中C4植物互花米草的GPP值明显高于其他类型植被,例如基于线性回归模型的GPP估算结果,在2020年8月16日长江口互花米草湿地的空间平均GPP为20.96μmol m-2 s-1,比芦苇的GPP高出约49%。综上所述,本研究结果表明,在进行滨海湿地GPP近瞬时尺度估算及其空间格局分析时,运用高空间分辨率的数据,准确使用遥感指数,并结合机器学习算法,将有利于获得更精确、更全面的估算结果。