诊断陡峭山地森林中GEDI冠层高度误差：地面高程表征与地理定位不确定性

《Forests》：Diagnosing GEDI Canopy Height Errors in Steep Mountainous Forests: Ground Elevation Representation and Geolocation Uncertainty

　　

星载激光雷达（LiDAR）提供了森林垂直结构的重要观测，然而在陡峭山地森林中，冠层高度反演仍易受地形相关误差影响，此类区域的地形异质性使得误差难以归因于足迹地理定位与地面高程表征。为区分这两种效应，研究人员以中国九寨沟地区1 m机载激光扫描（ALS）衍生的数字地形模型和冠层高度模型（DTM与CHM）作为局部参考。通过2 × 2诊断设计，在四个坡度等级中独立变化足迹地理定位与地面参考。对于2102个严格筛选的足迹，地理定位细化使均方根误差（RMSE）从13.36 m略增至13.47 m，而基于中位数的地面参考替换使RMSE降至9.98 m，降幅达25.3%。GEDI隐含地面在坡度低于35°时最接近足迹中位数，但在最陡地形中转向P₂₀–P₃₀。垂直参考敏感性分析显示，地面诊断的最优参考在P₂₀与P₃₀之间变化，而足迹中位数始终使校正后RH₉₈与ALS CHM P₉₈的RMSE最小。ATL03对比同样表明，在其他光子处理条件完全相同的情况下，以ALS DTM替换哥白尼DEM（Copernicus DEM）作为地形支持输入可显著改善冠层高度一致性。这些结果确定地面高程表征（而非所测试的水平地理定位细化）是陡峭地形中主导GEDI–ALS冠层高度分歧的地形相关因素，并表明地面诊断与RH₉₈校正需要不同的地形参考选择。

**陡峭山地森林中GEDI冠层高度误差的诊断：地面高程表征与地理定位不确定性**

**1. 研究背景与问题**

星载激光雷达（LiDAR）为量化森林垂直结构提供了前所未有的观测能力，其中全球生态系统动力学调查（GEDI）与冰、云和陆地高程卫星二号（ICESat-2）是当前两种主要传感器。然而，在陡峭山地森林中，地形起伏大、冠层异质性强，导致GEDI全波形信号与ICESat-2光子计数信号在解释地面和冠层返回时面临显著挑战。已有研究多将GEDI冠层高度误差归因于足迹地理定位不准或地面参考不一致，但很少有研究在相同严格筛选样本内定量区分二者的相对贡献。由于单个GEDI足迹（直径约25 m）内高程范围可能很大，波形推导出的地面高程未必对应足迹内地形分布的简单中心趋势，因此尤其需要明确：是水平定位误差还是地面高程表征误差主导了冠层高度反演偏差。为此，研究人员以中国四川九寨沟陡峭山地森林为研究区，利用高精度机载激光扫描（ALS）参考数据，设计了严格的诊断框架，以独立评估地理定位细化与地面参考替换对GEDI冠层高度精度的影响。该研究发表于《Forests》。

**2. 研究内容与主要技术方法**

研究人员利用2020年采集的1 m分辨率ALS数字地形模型（DTM）与冠层高度模型（CHM）作为局部参考，匹配经过严格质量筛选的2021年GEDI L2A足迹（共2102个）及少量ICESat-2 ATL03数据。主要方法包括：构建“Raw、Geo-only、Ground-only、Geo+Ground”四种诊断状态，分别表示原始状态、仅地理定位细化、仅地面参考替换、两者联合；采用基于地面高程匹配的两阶段地理定位细化流程；使用ALS DTM的六种足迹尺度统计量（最小值、P₀₅、P₁₀、P₂₀、P₃₀、中位数）作为候选地面参考；对ICESat-2光子数据实施双地形支持（哥白尼DEM与ALS DTM）的配对处理，并对全部结果进行坡度分层与稳健性检验。

**3. 研究结果**

**3.1 GEDI与ALS在不同坡度等级下的原始表现**
原始GEDI RH₉₈相对ALS CHM P₉₈的整体偏差和RMSE分别为10.03 m与13.36 m，且误差随坡度增大而显著增加，>35°坡度时RMSE高达16.24 m，表明陡峭地形中存在固有的坡度相关残余误差。

**3.2 地理定位细化对GEDI冠层高度一致性的影响**
利用地面匹配细化足迹位置后，整体RMSE从13.36 m略增至13.47 m，偏差与R2均未改善，各坡度等级亦无一致减差效果。即使改用P₂₀或P₃₀作为搜索目标，结论仍然不变，说明单纯的水平位置调整无法消除地形引起的冠层高度误差。

**3.3 地理定位细化与地面参考替换的相对效应**
地面参考替换显著改善了GEDI–ALS一致性：Ground-only将RMSE降至9.98 m，降幅25.3%，而Geo+Ground的RMSE为12.37 m，劣于Ground-only。坡度分层显示，地面参考替换在>35°时RMSE降低4.97 m，而地理定位细化在陡坡上反而使误差增大。平衡坡度样本与连续坡度分析均支持这一排序，且不受冠层结构混杂影响。

**3.4 DTM参考对GEDI地面诊断与RH₉₈校正的影响**
在名义垂直归一下，P₃₀在全体样本中最接近GEDI隐含地面，而在最陡坡度类中P₂₀–P₃₀更优。然而，校正RH₉₈时，足迹中位数始终产生最低RMSE（9.98 m），且这一优势在垂直偏移敏感性分析中保持稳定。这表明地面诊断与冠层高度校正的最优地形参考并不一致；低分位数参考更适用于解释波形-地形交互，而中位数更适合作为经验性校正参考。

**3.5 不同地形支持下的ICESat-2 ATL03比较**
在相同光子处理流程下，以ALS DTM替代哥白尼DEM作为地形支持后，ATL03 RH₉₈与ALS CHM P₉₈的RMSE从15.37 m降至9.00 m，R2从0.010升至0.697，陡坡上改善尤为明显，证实地形支持质量对光子计数激光雷达同样重要。

**3.6 稳健性检验与退化GEDI足迹的影响**
放宽degrade_flag筛选后（n=5328），地理定位细化仍无改善，但地面参考替换失效，整体RMSE增至15.29 m。这归因于degrade_flag=70的足迹具有不稳定的波形地面高程（标准差20.54 m），将噪声传播至校正后的RH₉₈。因此，地面参考替换的有效性依赖GEDI自身波形质量，严格的质量筛选不可或缺。

**4. 讨论总结与结论**

讨论指出，在严格筛选样本中，地面参考替换对GEDI–ALS冠层高度分歧的影响远大于所测试的地理定位细化，但这一结论受限于诊断性位移而非真实定位校正。GEDI隐含地面随坡度增大向低分位数移动，而中位数作为RH₉₈校正参考更为稳健。局限性包括：单一研究区、ALS参考并非真值、垂直基准采用EGM2008近似转换、时间间隔可能引入植被生长影响等。未来需在多种森林类型与地形条件下扩展，并引入实测足迹位置与野外冠层高度数据。研究结论如下：在九寨沟陡峭山地森林中，以ALS DTM和CHM为局部参考，研究人员评估了与地形相关的GEDI RH₉₈误差来源，并辅以ICESat-2 ATL03对比。在严格筛选样本中，所测试的基于地面匹配的地理定位细化未能改善冠层高度一致性，并略微增加了整体RMSE；而基于中位数的地面参考替换将RMSE从13.36 m大幅降至9.98 m。四种诊断状态的相对排序经聚类配对自举重采样支持，且在平衡坡度样本和控制冠层结构后保持不变。GEDI地面参考行为呈强烈坡度依赖性：在低至中坡度地形，GEDI地面高程与足迹中位数更一致；在最陡地形中，则偏向低分位数ALS DTM统计量（尤其是P₂₀–P₃₀）。然而，最匹配GEDI地面高程的DTM统计量并非使校正后冠层高度误差最小的统计量。对于RH₉₈校正，足迹中位数在该数据集中提供了最有效的经验参考，表明地面诊断与冠层高度校正需采用相关但不同的地形参考选择。ICESat-2结果进一步支持地形支持的重要性，但该实验并非平行的四状态分解，其结果显示地形支持输入的改变会显著影响光子计数冠层高度反演，尤其在陡峭地形中。放宽筛选的分析表明，地面参考替换无法可靠推广至波形地面高程不稳定的退化足迹，凸显了严格波形质量筛选的必要性。总之，陡峭地形的星载激光雷达冠层高度分析应明确考虑地形坡度、波形或光子质量、地形支持输入以及足迹尺度参考构建。低分位数DTM统计量有助于解释波形-地形交互，而足迹中位数在类似陡峭茂密森林中可作为RH₉₈校正的保守经验参考。