摘要
准确的高山植被物候遥感产品是维持生态平衡、正确反映地表能量平衡以及冰冻圈稳定的基础性数据。然而,目前常用于提取物候的遥感数据空间分辨率多样,基于不同空间尺度提取的物候可能存在尺度效应,但该问题目前并未受到广泛关注。本研究利用10m分辨率Sentinel-2,500m以及5600m分辨率 MODIS(Moderate Resolution Imaging Spectroradiometer)遥感数据,使用NDGI(Normalized Difference Greenness Index)植被指数,基于GEE(Google earth engine)云平台提取了2022年青藏高原地区相应空间分辨率的物候产品,探讨了不同空间分辨率物候差异,定量分析了物候产品的空间尺度效应,并阐明了其随地形复杂度的变化。研究发现,高分辨率遥感产品(10m)更能清楚地反映精细物候的变化。另外,使用不同空间分辨率遥感产品探测草地物候时,尺度效应普遍存在,其中,较粗空间分辨率春季物候存在滞后现象,且秋季物候的尺度效应小于春季物候。不同空间分辨率差异形成的尺度效应对地形变化的敏感度不同,总体上,物候的尺度效应随着地形复杂度的增加而增加。
背景
青藏高原具有集中且丰富的生物气候垂直带谱,在生物多样性,调节区域气候和涵养水源等方面扮演着重要的角色。植被的存在改变了局地气候,为青藏高原区域提供了丰富的生境和栖息地,是维持该生态系统中生物多样性不可或缺的重要存在,监测青藏高原植被物候可以及时发现和评估潜在的自然灾害风险,反馈冰冻圈的稳定性,保护青藏高原生态系统的多样性。然而,目前用于监测植被动态遥感产品的空间分辨率范围从米级到公里级不等,通常基于不同尺度遥感产品探测到的植被动态可能存在差异,以前的研究已充分地证明了这一点。但是目前为止,使用不同空间尺度数据进行物候提取时是否存在空间尺度效应尚不明确,这可能会造成使用单一空间分辨率遥感产品进行物候提取时得出的结论不够准确。因此,准确认识物候提取中的空间尺度效应至关重要。鉴于归一化差分绿度指数(Normalized Difference Greenness Index,NDGI)在提取植被物候时精度较高,本研究分别利用10m Sentinel-2、500m MOD09A.061以及5600m MOD09CMG.061遥感数据计算NDGI指数,并根据NDGI的时间序列提取青藏高原的草地物候,分析了不同空间尺度的物候差异,并在该结果的基础上计算了物候的空间尺度效应。另外,为了解青藏高原区域植被物候空间尺度效应的分布特征,探讨了空间尺度效应随地形的变化情况。
主要结果
由图1, 虽然不同空间尺度NDGI的空间分布特征相同,但春秋季物候SOS值本身的大小存在较大差异,尤其是公里级遥感产品(5600m)MOD09CMG.061提取的春季物候相较10m和500m春季物候较晚。这可能是因为粗空间分辨率像元涵盖的空间范围较大,而物候可能发生在较小的区域内,这一精细变化在聚合到粗空间分辨率时会被平均,因此,粗空间分辨率不能及时探测到物候的变化。随着时间的推迟,当变绿更加明显时,粗空间分辨率像元才会将此时植被绿度的变化认别为生长季的开始。对比不同空间分辨率EOS值的分布情况,发现公里级和百米级遥感产品秋季物候值的分布非常相似,即公里级空间分辨率遥感产品(5600m)探测到EOS的滞后现象不明显,这可能归因于以下两方面:(1)在植被绿度下降阶段,草地叶子的凋落和颜色变化较为明显,在粗空间分辨率上也能捕捉到;(2)秋季青藏高原较为干旱,云量较少,遥感影像的时间序列较为完整。
图1. 不同空间分辨率春季(SOS)和秋季(EOS)物候的空间分布,第一行(a-c)和第二行(d-f)分别指SOS和EOS的空间分布,第一列(a,d),第二列(b, e)和第三列(c, f)分别指的是10m(Sentinel-2),500m(MOD09A1.061)和5600m NDGI(MOD09CMG.061)物候的空间分布,Area指的是对应空间分辨率NDGI遥感产品提取的物候像元所占的面积,单位‘d’为‘day’的缩写,表示‘天’
由图2,使用不同空间分辨率遥感产品进行物候提取时尺度效应普遍存在。其中,500m和10m的尺度误差在青藏高原大部分区域分布较为均匀,但在青藏高原西南部区域,尺度误差较大。由图2b和图2c,5600m与10m以及5600m与500m的尺度误差在空间上分布具有相似性,这主要是因为10m和500m遥感产品提取的SOS差异较小的缘故。
图2. 不同空间分辨率遥感产品尺度误差的空间分布,第一行(a-c)和第二行(d-f)分别指SOS和EOS尺度误差的空间分布,第一列(a, d),第二列(b, e)和第三列(c, f)分别指500m和 10m,5600m和10m以及5600m和500m尺度误差的空间分布,单位‘d’为‘day’的缩写,表示‘天’。
对比不同空间分辨率尺度误差值的分布情况(图3),发现500m和10m SOS尺度误差最小,约75%的区域尺度误差小于等于12d,当空间分辨率差异为5600m和10m以及5600m和500m时,约50%的尺度误差位于[0d, 12d]的区间范围内。这一现象可能表明了空间尺度越接近,春季物候的空间尺度就越小。基于不同空间分辨率EOS的尺度误差在青藏高原中南部区域最大,其他大部分区域尺度误差较小。三种空间分辨率差异引起的尺度误差均有超过75%的区域小于12d,因此,与SOS相比,EOS物候的空间尺度效应较小。
图3. 尺度误差在不同数值区间上的概率统计,(a)图和(b)图分别指的是SOS和EOS尺度误差的概率统计情况,第一行到第三行分别指空间分辨率差异为500m和10m,5600m和10m以及5600m和500m尺度误差的分布。
由图4,由不同空间尺度差异引起的尺度误差随CV增大显著增加,其中5600m与500m空间分辨率SOS的尺度误差对CV敏感性较小,CV每增加0.1,SOS的尺度误差增加0.34d。5600m与10m SOS的尺度误差对CV变化最敏感,CV每增加0.1,SOS的尺度误差便增加0.74d。同SOS尺度误差随CV的变化一致,5600m与500m空间尺度差异引起的EOS尺度误差对CV的敏感性最弱,CV每增加0.1,尺度误差则增加0.49d。500m与10m EOS的尺度误差对地形最敏感,CV每增加0.1,物候偏差增加0.78d。尺度效应随地形复杂度增大而增大的原因是复杂的地形形成了多种微气候环境,使得植被生长环境各异,导致不同区域植被物候差异较大,这种差异能够在精细尺度的遥感产品中表现出来,但是当遥感产品分辨率较低时,这种变化就会在聚合的过程中被平均。当地形较为平坦,即地形复杂度较小时,植被异质性较小,基于不同空间尺度得到的物候较为相近,因此,这种情况下物候尺度效应较小。
图4. 由空间分辨率不同引起的物候偏差随地形复杂度的变化。第一到第三列(a-c)分别是500m与10m, 5600m与10m以及5600m与500m的尺度误差随地形复杂度的变化。蓝色和橘色箱线图分别指的是SOS和EOS的在不同地形复杂度区间内统计情况,箱线图上的黑色圆圈对应每个柱子的均值,蓝色和橘色虚线是对SOS和EOS均值线性拟合的结果。slope代表的是拟合线的趋势,p指的是显著性水平。
文章信息:马杜娟,谢江流,陈瑞,等.青藏高原草地物候的空间尺度效应[J].测绘学报,2024,53(5):860-868. DOI:10.11947/j.AGCS. 2024.20230256.
供稿:马杜娟