格陵兰西北部Qaanaaq冰盖冰川藻类物候动态及其对冰面暗化过程的高分辨率遥感研究

《Journal of Glaciology》：Glacier algae phenology on the Qaanaaq Ice Cap (Northwest Greenland)

本研究针对格陵兰冰盖暗化现象中冰川藻类的关键作用，通过多源遥感与实地观测相结合，建立了基于Sentinel-2波段比（B5/B4）的藻类丰度反演模型，揭示了2016-2023年间Qaanaaq冰盖藻类时空异质性分布规律，发现持续正温与降雪间歇是藻华爆发的关键驱动因素，为评估生物效应对极地冰盖物质平衡的影响提供了重要方法论支撑。

　　

随着全球气候变暖加剧，格陵兰冰盖（GrIS）正经历着显著的表面暗化现象，这种现象直接导致冰面反照率降低，从而加速冰体消融。过去几十年间，格陵兰冰盖特别是其西南部所谓的“暗带”区域，表面反照率持续下降，这种现象被归因于多种因素的综合作用：降雪频率减少、太阳辐射增强导致雪粒老化、黑碳和矿物粉尘等吸光性颗粒物的积累，以及生物活动的增强。其中，冰川藻类因其含有深色色素（如酚类红紫棓精衍生物），能够强烈吸收太阳辐射，成为冰面暗化的重要生物驱动因子。藻类生长依赖于冰雪融化提供的水分和养分，进而形成融化-反照率正反馈机制：温度升高和融化促进藻华爆发，藻华又进一步降低反照率，加剧融化。在这一背景下，冰面风化壳的形成为藻类提供了理想栖息地，而特定蓝藻则促进有机质形成，导致冰面上出现深色聚集体（即冰尘）。尽管前人研究已认识到冰川藻类在冰盖暗化中的重要作用，但对于高纬度地区如格陵兰西北部Qaanaaq冰盖的藻类时空动态及其与气象、地形因子的关系，仍缺乏系统性的高分辨率研究。

为深入探究冰川藻类在Qaanaaq冰盖的物候特征及其环境驱动机制，研究团队整合了多源遥感数据与实地观测，通过建立藻类丰度遥感反演模型，分析了2016-2023年间藻类的时空变化规律，并探讨了其与气象、地形因子的关联。本研究主要依托以下关键技术方法：首先，基于2014年和2023年野外采集的53个冰面样本，通过显微镜直接计数法测定藻类细胞浓度（cells mL^-1），并采用燃烧法量化粉尘含量；同时利用野外光谱仪（MS-720和RoX）测量地表反射率，波段范围覆盖400-865 nm。其次，通过相关性矩阵分析筛选出对藻类丰度最敏感的反射率比值（695/681 nm或695/687 nm），进而建立基于Sentinel-2波段5（705±15 nm）与波段4（665±30 nm）比值的藻类丰度反演模型。再次，利用Sentinel-2影像（10 m分辨率）提取冰盖反照率及藻类丰度，并结合ERA5-Land再分析数据（9 km分辨率）和ArcticDEM地形数据（32 m分辨率），通过物理降尺度方法获取气象要素（气温、降雪）的高分辨率分布。最后，采用物候学分析方法，量化了裸冰期始末（SOS/EOS）、藻华持续时间（LObS）等关键参数，并在Qaanaaq冰川不同海拔设置关注区（AOI），分析藻类丰度与气象、地形因子的关系。

Glacier algae and dust abundance from field observations

通过对比2014年和2023年野外样本发现，藻类丰度存在显著差异：2023年平均藻类浓度（8.1±5.5×10⁵ cells mL^-1）比2014年（2.2±2.4×10⁴ cells mL^-1）高出一个数量级，且两年均以Ancylonema nordenskioldii为优势物种（占比＞82%）。藻类分布沿海拔呈现规律性变化，中海拔区域（400-700 m）浓度最高，而冰川末端因坡度较陡（＞10°）、表层径流冲刷作用强，藻类丰度较低。粉尘含量方面，无机粉尘占总颗粒物的86-96%，有机质占比仅4-14%，表明藻类生长与粉尘积累存在耦合关系，但模型分析显示反射率比值对粉尘不敏感。

Field spectroscopy data and reflectance ratio

野外光谱测量显示，2023年冰面反照率（最低0.13）显著低于2014年，且光谱曲线在680 nm附近出现明显的叶绿素a吸收谷。通过216,225个反射率比值与藻类丰度的相关性分析，确定695/681 nm或695/687 nm比值与藻类浓度相关性最高（R²=0.90）。Sentinel-2波段5/波段4比值与藻类丰度的线性回归模型为：AA = (8.04×10⁶ × Ratio) - 7.56×10⁶（R²=0.84），该模型对粉尘无显著响应（R²=0.03, p>0.05），验证了其特异性。

Interannual variability over the Qaanaaq Ice Cap by means of glacier phenology metrics

2016-2023年物候分析表明，藻类丰度存在显著年际差异：2019、2020和2023年为藻华高发年，夏季平均藻类丰度（mean(AA)）＞1×10⁶ cells mL^-1，相当于等效碳储量最高达2114 kg C（2020年）；而2017、2018和2021年因持续降雪和低温，藻华规模较小。空间分布上，东部分水岭区域藻类丰度较西部分水岭高0.5-1.7×10⁶ cells mL^-1，且中海拔区域（400-500 m）是藻华热点区。反照率最低值（min(α)＜0.20）出现在2019年，与藻类高丰度区域吻合。

Role of meteorological and glaciological parameters on the algae abundance variations at Qaanaaq Glacier

气象因子分析显示，连续正温天数与降雪间歇是藻华爆发的关键驱动因素。2019和2023年分别出现57-64天连续正温期，藻类丰度峰值超过10×10⁵ cells mL^-1；而2017、2018年正温期较短（约30天），且降雪事件频繁，藻华规模受限。线性回归表明，在400-600 m海拔区间，藻类丰度与反照率（R²=0.70）和气温（R²=0.40）显著相关，但与降雪量的关系不显著。地形分析进一步证实，陡坡区域（如冰川末端）因表层径流冲刷作用，藻类丰度较低。

Limitations on glacier algae abundance estimation

本研究提出的反演模型仅适用于以Ancylonema nordenskioldii为优势种的冰面环境，且对藻类浓度较低（反射率比值＜0.94）区域的敏感性不足。Sentinel-2数据的时空分辨率受高纬度云覆盖影响，可能引入不确定性。未来可结合更高光谱分辨率卫星（如PRISMA、EnMAP）或每日重访卫星（如Sentinel-3、GCOM-C）提升监测能力。

本研究系统揭示了Qaanaaq冰盖冰川藻类的物候特征及其环境驱动机制，建立了基于高分辨率遥感数据的藻类丰度定量反演方法，填补了高纬度地区生物暗化研究的空白。研究发现，持续正温与降雪间歇是藻华爆发的关键条件，而地形坡度通过影响表层水文过程调节藻类空间分布。这些结论不仅深化了对冰川生物-物理耦合过程的理解，也为预测未来气候变暖背景下极地冰盖物质平衡变化提供了科学依据。论文所提出的遥感物候分析方法可扩展至整个格陵兰冰盖，为评估生物效应对全球海平面上升的贡献提供方法论支持。