2017年北美日食期间光谱植被指数与太阳诱导荧光对光照强度与几何变化响应的研究综述

本研究的主要作者为 Cheryl A. Rogers（第一作者）、Jing M. Chen、Ting Zheng、Holly Croft、Alemu Gonsamo、Xiangzhong Luo 和 Ralf M. Staebler，其所属机构分别来自加拿大多伦多大学、美国威斯康星大学麦迪逊分校、英国谢菲尔德大学、加拿大麦克马斯特大学、美国劳伦斯伯克利国家实验室以及加拿大环境与气候变化部。该研究成果于2020年发表在《Journal of Geophysical Research: Biogeosciences》期刊上。

一、 研究背景与目的

本研究隶属于遥感科学与全球变化生态学的交叉领域。准确监测陆地生态系统光合作用是理解全球碳循环和气候变化减缓潜力的关键。遥感技术，特别是各类光谱植被指数（Spectral Vegetation Indices），已成为获取连续空间尺度植被信息（如健康状况、生物物理性状）并进而模拟光合作用的核心手段。例如，归一化差异植被指数（NDVI）、光化学反射指数（PRI）和太阳诱导荧光（SIF）等指数，分别被用作植被结构、光合生理状态和光合作用本身的替代指标。然而，这些光谱信号在实际观测中受到光照强度（由云、气溶胶、日变化等引起）和光照/观测几何角度变化的共同影响，这两种因素通常在自然条件下相互耦合，难以分离，尤其是在整个森林冠层尺度上。这种耦合效应严重干扰了基于光谱指数准确反演植被生理、生化及结构参数的能力。因此，厘清光照强度与几何角度对光谱指数的独立影响，对于提高遥感数据的解释精度至关重要。

基于此背景，本研究团队创新性地利用了2017年8月21日发生的北美日食这一独特的自然事件，作为在冠层尺度上进行的一次“可控”遮光实验。日食过程在短时间内造成光照强度大幅减弱，但太阳在天空中的位置（即光照几何角度）变化极小。这为研究者提供了一个前所未有的机会，能够将光照强度的影响从光照/观测几何角度的复杂影响中剥离出来。研究设置了三天观测期：日食前一日（晴朗）、日食当日（晴朗，含日食过程）和日食后一日（多云），旨在对比不同光照条件下（晴朗、日食遮光、多云扩散光）光谱指数的响应。

本研究的具体目标包括：1) 确定日食期间及整个研究期内环境条件和植被的响应；2) 评估生理相关指数（PRI和SIF）在日食遮光以及光照随几何角度同步变化时，对光照减弱的响应；3) 评估用于估计冠层结构和叶片生物化学的指数在日食、日出日落及多云造成的光照变化下的稳定性；4) 探索冠层空间异质性和天空状况对光谱指数变异性的影响。

二、 详细研究方法与流程

1. 研究地点与数据收集： 研究地点位于加拿大安大略省的博登森林研究站，该站点设有一座装备齐全的涡度相关通量塔，提供长期的二氧化碳、能量和水汽通量数据。森林为混交林，主要树种包括红枫、北美乔松、大齿杨等。数据收集于2017年8月20日至22日，涵盖了一个晴朗日、一个包含日食的晴朗日和一个多云日。

2. 环境与植被生理数据获取： 从通量塔获取了核心气象与生理数据，包括：光合有效辐射（PAR）、空气温湿度、风速风向、二氧化碳净通量（NEP）。此外，还使用热比率法测量了临近塔的红枫树的树干液流速率，作为植被蒸腾和水分状况的指示。

3. 冠层光谱观测与光谱指数计算： 研究的核心数据来自安装在塔上38米高处的双通道光谱仪（UniSpec-DC）。该设备的一个探测器向上测量天空辐照度，另一个探测器向下，通过一个云台系统进行多角度观测，测量植被冠层的反射辐射。观测周期为每15分钟完成一轮四个不同天顶角（43°, 48°, 53°, 58°）和一系列方位角的测量循环，整个研究期间共获得超过33,000条光谱反射率测量数据。

数据处理流程包括：基于温度的暗电流校正、上下探测器的交叉定标（确保两个通道的响应一致性）、利用高光谱数据计算各个波段的反射率。研究基于校正后的光谱数据计算了六种植被指数，并将其分为三类： * 生理指数：光化学反射指数（PRI，旨在追踪光胁迫下的叶黄素循环变化）和太阳诱导荧光（SIF，利用氧气A吸收带760 nm附近的重填效应计算，作为光合作用的直接指示）。 * 结构指数：归一化差异植被指数（NDVI）、增强型植被指数（EVI，对高生物量区更敏感）和植被近红外反射率（NIRv，被证明与SIF和GPP密切相关）。 * 叶绿素指数：选择Maccioni等人（2001）开发的指数（MACC01），该指数在林冠尺度叶绿素含量反演中表现最佳。

对于SIF的计算，研究采用了基于UniSpec-DC光谱仪在氧气A吸收带（约760 nm）及其临近波段（756.7 nm）的简化反演算法，以获得SIF的相对变化值。虽然该仪器光谱分辨率有限，但此前研究已证明其可用于追踪SIF的日变化趋势。

4. 数据分析与建模流程： 分析遵循以下关键步骤： * 时间序列分析：将每15分钟完成的所有角度观测结果进行平均，以减弱冠层空间异质性的影响，然后分析各光谱指数随时间及PAR的变化趋势。 * 回归分析：针对晴朗、日食期间、多云三种天空状况，分别建立各光谱指数与PAR之间的线性回归模型，计算决定系数（R²）和斜率，以量化光照强度对不同天空条件下光谱指数的影响程度和模式。 * 光谱差异分析：计算日食开始和日食峰值时刻平均冠层反射光谱的差值，分析在整个可见光-近红外波段范围内，遮光引起的反射率变化的谱形特征，以探究驱动光谱指数变化的物理机制。 * 几何角度效应分离：为分离观测几何效应，计算了每次观测时太阳与传感器视线方向之间的相位角。首先，分析平均后的光谱指数与相位角的关系。其次，进行更精细的分析：将未平均的原始单次测量数据（3秒一次）按光束分数（光束辐射占总PAR的比例，反映天空漫射程度）和传感器观测方位角进行分组。然后，在每个分组内建立光谱指数与PAR的回归模型，并绘制所有分组回归的R²和斜率图，以可视化空间异质性（不同观测方向对应不同冠层区域）和天空状况如何共同影响光谱指数对光照的响应。最后，在模型中额外加入相位角作为解释变量，检验其是否显著提高了对光谱指数变异的解释能力。

三、 主要研究结果

1. 环境与生理响应： 日食期间，PAR急剧下降，导致空气温度降低、相对湿度升高。二氧化碳净吸收通量（NEP）和树干液流速率均出现明显下降，证实了遮光对森林光合作用和蒸腾作用的生理抑制效应。

2. 光谱指数的时间动态与对PAR的响应： * 生理指数（PRI和SIF）：如预期所示，它们对光照变化高度敏感。日食遮光期间，SIF随PAR下降而显著降低，SIF产额（SIF/PAR）则略有上升；PRI则显著增加（表明光胁迫缓解）。回归分析表明，在日食期间，PAR可以解释PRI变异的99%和SIF变异的97%，相关性极强。然而，在晴朗日，PRI与PAR的关系呈现明显的滞后现象：早晨的PRI值远高于傍晚相同光照水平下的值，导致整体相关性较弱。这种滞后在日食日有所减弱，可能与日食缓解了午后典型的高温和水分胁迫有关。 * 结构指数（NDVI, EVI, NIRv）：在日食期间，这些指数对PAR变化响应微弱，保持相对稳定，符合其反映相对稳定冠层属性的预期。但在晴朗日，NDVI与PAR呈显著正相关，这主要是由日出日落时极低的PAR和NDVI值驱动，分析表明这更可能与此时观测几何接近“热点”方向（观测方向与太阳光照方向一致）导致阴影比例最小、NDVI降低的几何效应有关。 * 叶绿素指数（MACC01）：在日食期间，MACC01与PAR呈显著负相关（R²=0.91），但其绝对值变化范围很小。在晴朗和多云条件下，其与PAR的关系模式不同且相关性较弱，表明其受光照角度影响较大。

3. 光谱变化机制分析： 对日食始末反射光谱差异的分析揭示了关键发现。反射率的变化并非均匀分布在所有波长。在整个可见光波段（450-700 nm），反射率普遍增加，但在绿光波段（约550 nm，叶绿素吸收谷）增加幅度小于红光和蓝光波段。在近红外波段（700-800 nm），反射率普遍下降，但在760 nm（SIF探测波段）附近的下降幅度小于周边波段。 * 最重要的发现是：并未观察到预期中由叶黄素循环引起的、以531 nm为中心的狭窄反射峰。相反，观测到一个以550 nm为中心的较宽光谱特征。这意味着，日食期间PRI的升高，主要是由于其在570 nm处的参考波段反射率相对于531 nm信号波段下降得更明显所致，而非531 nm处反射率的生理性上升。这一结果表明，冠层尺度的PRI信号变化可能更多地由冠层内多次散射的物理过程变化驱动，而非纯粹的叶片水平叶黄素色素转换的生理响应。 * 光谱变化模式也解释了其他指数的行为：NDVI、EVI、NIRv所使用的波段（如680 nm和800 nm）受到日食遮光的相似影响，故指数稳定；而MACC01使用的波段（680, 710, 780 nm）受影响程度不同，尤其在680-710 nm间光谱斜率变化大，导致其对光照变化表现出（尽管幅度小）可探测的敏感性。

4. 几何角度与空间异质性的影响： * 相位角分析证实，在晴朗日，NDVI和SIF与相位角存在明显关系，但在日食期间（相位角变化°），这种关系被最小化。 * 分组回归分析（按光束分数和观测方位角分组）显示，对于SIF和PRI，模型拟合优度（R²）在日食条件下最高且最稳定。在晴朗和多云条件下，拟合优度和回归斜率随观测方位角（即冠层不同区域）变化较大，表明冠层空间异质性（树种组成、空隙等）显著影响了光谱指数对光照的响应关系。添加相位角作为解释变量通常能提高模型的R²，但其重要性在不同观测方向和天空条件下表现不一。

四、 研究结论与价值

本研究通过巧妙利用日食这一自然实验，首次在真实森林冠层尺度上，将光照强度的影响从光照/观测几何角度的耦合效应中有效分离。主要结论如下：

1. 生理相关光谱指数（PRI和SIF）对光照强度变化高度敏感，但这种敏感性在自然条件下常被几何角度效应和冠层异质性所混淆。日食提供了最“纯净”的强度-响应关系。
2. 冠层尺度PRI信号对遮光的响应机制可能与预期不同。研究发现，日食引起的PRI升高主要由570 nm参考波段反射率的相对变化驱动，而非531 nm信号波段的生理性上升。这揭示了冠层多次散射的物理过程在塑造PRI观测信号中的重要作用，这对准确解释PRI作为光能利用效率（LUE）指标的可靠性提出了重要警示。
3. 结构指数（NDVI, EVI, NIRv）在光照强度独立变化时（日食）相对稳定，验证了它们反映较稳定冠层属性的能力。但其在自然日变化中表现出的变化，更多地归因于随太阳位置变化的观测几何效应（如阴影比例变化），而非强度本身。
4. 漫射光照条件（多云）引入了额外的复杂性。在这种条件下，光在冠层中穿透更深，使得光谱指数与PAR的关系变得更弱、更分散，表明漫射光改变了冠层内光分布和叶片贡献比例，其效应无法仅用PAR强度来解释。
5. 冠层空间异质性是导致光谱指数响应关系随观测角度变化的重要因素，在利用点尺度测量推演区域信息时需谨慎考虑。

本研究的科学价值在于，它通过一个独特且不可多得的自然实验，深化了对遥感植被指数响应机制的理解，特别是揭示了冠层物理过程（多次散射）对生理指数信号的潜在干扰。这对于改进基于PRI、SIF等指数的大尺度光合作用遥感估算模型、开发角度校正算法、以及理解不同天空条件下遥感数据的可比性具有重要的理论意义。同时，研究也展示了利用未来日食事件（如2024年北美日全食）进行进一步验证和探索的潜力。

五、 研究亮点

1. 新颖的实验设计：创造性地利用日食作为“自然遮光实验”，在真实森林冠层尺度上实现了光照强度与几何角度的“准可控”分离，这是传统人工遮光实验或模型模拟难以实现的。
2. 重要的机制发现：通过高光谱差异分析，首次在冠层尺度上明确指出，日食导致的PRI变化主要源于参考波段（570 nm）的物理响应，而非信号波段（531 nm）的生理响应，挑战了对PRI生理信号纯度的传统认知，强调了冠层辐射传输过程（多次散射）的关键作用。
3. 系统性的对比分析：在同一地点、短时间内，对比了晴朗、日食遮光、多云三种典型光照条件下的多种光谱指数响应，全面评估了不同类别指数对光照变化和几何效应的敏感性。
4. 精细化的数据解析：采用多角度观测并结合光束分数和空间分组回归分析，深入剖析了冠层异质性和天空状况对指数-光照关系的影响，提供了更细致的见解。