基于无人机平台的多角度植被冠层遥感观测前沿研究

《Forests》：Advanced Multi-Angle Remote Sensing Observation of Vegetation Canopy Leveraging UAV Platform

　　

山地集水区整合了环境梯度、干扰遗留效应及多尺度过程，这些过程塑造了生态系统功能在时空上的变异。为更好理解此类空间过程的作用，研究人员分析了巴塔哥尼亚（Patagonia）源头集水区中生态系统功能关系是否按单一多功能梯度组织，或是否对生态系统状况与结构具有尺度依赖性。研究人员基于三个空间分辨率——集水区（n=12）、森林簇（forest cluster，n=63）和森林样地（forest plot，n=175）——的野外实测生态系统功能，拟合了三个专家约束贝叶斯网络（Bayesian network，BN）。预测性能在尺度与功能间差异显著：集水区尺度下薪材体积预测最佳（r=0.79），而死木质碳储量、活树碳储量、维管植物丰富度、用材体积和土壤侵蚀呈负向预测相关（r分别为?0.14、?0.15、?0.33、?0.44和?0.80）；森林簇尺度下薪材体积、活树碳储量、土壤侵蚀和死木质碳储量预测较强（r=0.96、0.81、0.77、0.69）；样地尺度下薪材体积、树碳储量、死木质碳储量和林下植物多样性局部信号最强（r=0.91、0.67、0.67、0.60）。贝叶斯网络揭示出反复出现的木材–碳正向结构，联结树碳储量、薪材体积、死木质碳储量、树碳固存及部分土壤响应。相比之下，林下植物多样性更多与海拔、郁闭度、坡度和权属等局部驱动因子关联，侵蚀调控随尺度和森林发育阶段改变方向。结果表明，监测应跨嵌套空间尺度同步评估碳与木材生产相关功能及生物多样性与土壤相关功能，因为权衡与协同作为由共享环境与经营驱动因子塑造的情境与尺度依赖关系而涌现。

研究背景与立项依据：山地源头集水区将气候、地形、植被、土壤与生物过程耦合于邻近空间，森林生态系统功能（如碳储存、生物量生产、土壤发育、侵蚀调控、死木质累积与植物多样性）并非由单因子控制，而是同时受气候、地形、森林结构、土壤条件与干扰影响。已有研究指出，单一功能高低不足以推断系统整体运行状态，功能间可呈协同或权衡，且简单两两相关无法区分“直接功能联结”与“共同驱动介导关联”。现有森林多功能研究多依赖制图或模型指标，缺乏基于嵌套野外实测、跨集水区–森林簇–样地三尺度并同时纳入环境/干扰驱动与条件依赖结构的分析。为此，研究人员以智利巴塔哥尼亚艾森（Aysén）地区12个源头集水区为对象，采用专家约束贝叶斯网络刻画生态系统功能关系的尺度依赖性。

主要技术方法概括：研究人员沿用Moreno-Meynard公布的系统性嵌套野外清查数据，在干湿两种气候背景与参照/受扰（火、择伐、放牧）两种遗留下选取12个集水区（面积0.5–2.8 km²），每集水区内布设森林簇（含至多3个500 m²样地及50 m线性样线），最终集水区尺度n=12、森林簇尺度n=63、森林样地尺度n=175。八项共有功能节点为树碳储量、树碳固存、死木质碳储量、用材体积、薪材体积、土壤形成、土壤侵蚀（露根代理反转作侵蚀调控）、维管植物丰富度；协变量含生态系统背景、海拔、坡度、郁闭度、南向比例、森林发育阶段、干扰类型与强度、权属等。结构学习采用bnlearn包爬山算法（hill-climbing），500次bootstrap估弧强，按尺度设生态黑名单（如禁功能节点指向预测因子、禁维管丰富度→死木质碳）；连续节点作log(x+1)变换与标准化，预测性能以10折交叉验证的预测相关（predictive correlation, r）评价，条件协同/权衡通过固定条件因子下树碳储量梯度（0–1归一）的响应剖面bootstrap重采样获得。

研究结果

3.1 模型1集水区尺度生态系统功能网络：最终保留17条弧，预测异质性强；薪材体积r=0.79最高，土壤形成r=0.45、树碳固存r=0.36次之，死木质碳r=?0.14、树碳r=?0.15、维管丰富度r=?0.33、用材体积r=?0.44、土壤侵蚀r=?0.80为负。干集水区树碳高于湿集水区（+51.94 t C ha^?1，std. β=+0.52）但土壤形成更低（?377.79 m³ ha^?1）；受扰集水区土壤形成低于参照（?661.87 m³ ha^?1）。树碳与薪材体积、死木质碳呈正局部效应，与维管丰富度负向；树碳固存主要受薪材体积正向驱动，与树碳/土壤形成负向；土壤形成受死木质碳、用材体积、土壤侵蚀正向联结。结论：集水区尺度网络可识别广域功能结构，但多数功能预测不可靠，宜作受限概率综合而非普适结构。

3.2 模型2森林簇尺度生态系统功能网络：保留11条弧，预测显著提升，薪材体积r=0.96、树碳r=0.81、土壤侵蚀r=0.77、死木质碳r=0.69；维管丰富度仅r=0.30且由海拔等环境链入。干簇土壤形成低于湿簇（?843.99 m³ ha^?1）；海拔每升100 m死木质碳降3.04 t C ha^?1、维管种数降1.15。森林发育阶段分类显示老熟林树碳（170.74 t C ha^?1）高于成年林（204.60）系原文数值分布所致（次生林127.01、幼树簇3.17）。木材–碳块清晰：树碳↑→薪材体积↑（β=+0.90）、树碳固存↑（β=+0.29）、土壤侵蚀↑（β=+1.51）；薪材体积↑→死木质碳↑（β=+0.43）。结论：森林簇尺度在生态分辨率与预测稳定性间最佳，是监测木–碳–土模块首选。

3.3 模型3样地尺度生态系统功能网络：保留19条弧，薪材体积r=0.91、死木质碳r=0.67、树碳r=0.67、维管丰富度r=0.60；土壤侵蚀r=0.45、用材r=0.41、土壤形成r=0.29、固存r=0.27。维管丰富度随郁闭度（每+1%减0.064种）、海拔（每+100 m减0.77种）、坡度（每+1%减0.041种）下降，小地主权属最高（11.40种）、大地产权最低（2.83种）。海拔↑→树碳↑7.10 t C ha^?1但死木质碳↓5.49；郁闭度↑→土壤形成↑6.27 m³ ha^?1。样地木材块正向：树碳↑→固存↑、用材↑、薪材↑；薪材↑→死木质碳↑→土壤形成↑；但薪材与用材局部负向（β=?0.16），固存与用材负向。结论：样地模型为敏感性分析，确认木–碳块跨尺度稳定，同时显化维管丰富度局部驱动。

3.4 跨尺度综合：薪材体积与树碳为多尺度高连接核心；死木质碳–土壤形成在簇/样地稳为正向；维管丰富度不入核心碳生产块而走环境链。系统未坍缩为单一多功能轴，而呈部分联结子系统。

3.5 沿树碳梯度的条件权衡与协同：集水区尺度薪材、死木质碳、固存随树碳↑而↑，维管丰富度↓，侵蚀调控多↑（干参照略↓）；森林簇尺度薪材一致↑，侵蚀调控在次生/成年林↓、老熟林↑；样地尺度薪材、死木质碳、固存、用材、土壤形成↑，维管平坦或老熟林略↓，侵蚀调控次生林强↓、老熟林↑。证明协同/权衡非固定配对，而依尺度、发育阶段与条件因子变向。

讨论部分总结：研究人员指出，功能间关联须区分直接链路与共享驱动（气候、发育、海拔、郁闭、干扰）介导，两两相关会误判；最稳协同为树碳–薪材体积，死木质碳随尺度并入木–碳–土块，固存与储量不等价（速率为峰值可移），维管丰富度与碳块呈条件背离而非普适权衡，侵蚀调控随发育阶段变向反映次生林恢复未及老熟林根–枯落–死木防护结构。尺度上集水区过粗（n=12）仅作综合，森林簇最宜监测 recurrent 木–碳–土模块，样地补局部驱动。管理含义为碳/木材指标不能替代生物多样性与土壤功能直接监测，目标决定尺度与指标选择。局限含集水区样本小、观察性非因果、嵌套非独立、未测微气候/管理强度等潜变量。

结论部分翻译：最终模型表明，巴塔哥尼亚源头集水区生态系统功能关系具条件性、驱动介导性与尺度依赖性。碳与木材相关功能构成最反复出现的正向模块，而维管植物丰富度与侵蚀调控呈更强情境依赖。结果支持对权衡与协同的审慎解读：不应作为固定两两相关，而应视为依赖驱动因子、条件因素与观测空间尺度的关系。对管理而言，碳与木材生产指标可描述森林功能的重要部分，但若目标是维持多功能山地集水区，须与生物多样性及土壤功能同步评估。监测方案须使数据类型与尺度匹配管理情境与监测目标。

（全文登于《Forests》）