留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于机载P波段全极化SAR数据的森林地上生物量估测

姬永杰 杨丛瑞 张王菲 曾鹏 张甫香 屈亚妮

朱文见, 张慧, 王懿祥. 采伐对森林土壤呼吸影响的研究进展[J]. 浙江农林大学学报, 2021, 38(5): 1000-1011. DOI: 10.11833/j.issn.2095-0756.20210365
引用本文: 姬永杰, 杨丛瑞, 张王菲, 等. 基于机载P波段全极化SAR数据的森林地上生物量估测[J]. 浙江农林大学学报, 2022, 39(5): 971-980. DOI: 10.11833/j.issn.2095-0756.20220111
ZHU Wenjian, ZHANG Hui, WANG Yixiang. Research progress on effects of cutting on forest soil respiration[J]. Journal of Zhejiang A&F University, 2021, 38(5): 1000-1011. DOI: 10.11833/j.issn.2095-0756.20210365
Citation: JI Yongjie, YANG Congrui, ZHANG Wangfei, et al. Forest above ground biomass estimation using airborne P band polarimetric SAR data[J]. Journal of Zhejiang A&F University, 2022, 39(5): 971-980. DOI: 10.11833/j.issn.2095-0756.20220111

基于机载P波段全极化SAR数据的森林地上生物量估测

DOI: 10.11833/j.issn.2095-0756.20220111
基金项目: 国家自然科学基金资助项目(32160365,31860240,42161059);云南省万人计划青年拔尖人才项目(80201444)
详细信息
    作者简介: 姬永杰(ORCID: 0000-0001-8012-4115),讲师,博士,从事SAR、LiDAR数据资源环境反演研究。E-mail: jiyongjie@live.cn
    通信作者: 张王菲(ORCID: 0000-0002-2147-5246),教授,博士,从事农林业遥感研究。E-mail: mewhff@163.com
  • 中图分类号: S757

Forest above ground biomass estimation using airborne P band polarimetric SAR data

  • 摘要:   目的  森林生物量的空间精准量化对了解陆地碳储量、碳收支、碳平衡,以及揭示森林碳储量与全球气候变化的影响过程具有重要意义。P波段波长较长,在森林中具有更高的穿透能力,研究机载P波段SAR数据提高森林地上生物量(AGB)估测精度的可行性。  方法  以机载P波段全极化合成孔径雷达(SAR)数据和高精度激光雷达(LiDAR)数据估测的森林AGB抽样点为基础,提取20个极化SAR特征,并分别与森林AGB变化作敏感性响应情况分析。采用多元线性回归模型(MLR)、K近邻方法(KNN)、支持向量回归(SVR)和随机森林(RF)4种估测方法,探究机载P波段SAR数据的森林AGB估测精度。  结果  在较低森林AGB(均值约45 t·hm−2)的森林覆盖区中,P波段的同极化后向散射系数、Freeman-Durden和Yamaguchi分解中的表面和二次散射分量对森林AGB变化敏感;此外H-A-ALPHA极化分解的散射角(alpha)、拓展极化参数极化辨别率参数(PDR)也对森林AGB变化敏感。4种方法估测的森林AGB相对误差均约30%,其中MLR估测结果精度最低,估测精度为63.55%,均方根误差(RMSE)为19.16 t·hm−2;RF估测结果精度最高,估测精度为72.97%,RMSE为15.98 t·hm−2;KNN和SVR估计结果差别不明显,RMSE分别为17.04和17.09 t·hm−2  结论  P波段SAR数据对估测森林AGB具有一定潜力,3种非参数方法的估测结果明显优于MLR参数方法。此外,P波段的森林AGB估测精度受到待估森林AGB水平高低的影响明显,在森林AGB水平较高的分组中估测精度较高。在森林AGB均值为45 t·hm−2,最大值为120 t·hm−2的森林覆盖区,以50 t·hm−2将森林AGB样点分为2组时,高森林AGB组的估测精度高出低AGB组约6%。图5表3参34
  • 全球土壤碳储量约为1 500 Pg,超过全球陆地植被碳储量和大气碳储量之和[1]。土壤呼吸年均释放80~95 Pg二氧化碳-碳(CO2-C)到大气中[2-3],是化石燃料燃烧产生二氧化碳排放量的11倍以上[4-5],是陆地生态系统第二大碳通量。全球范围内,森林在减缓气候变化方面发挥着重要的作用[6]。作为地下生态过程的土壤呼吸显著影响着陆地生态系统的碳循环[7],其通量过程已成为全球变化生态学研究的核心和焦点之一。一方面,大气中CO2等温室气体的增加是导致全球气候变化的主要原因。另外一方面,全球气候变化也会加速土壤呼吸速率,进一步增加CO2年排放量。CO2排放与全球气候变化的正反馈作用将放大全球气候变化对陆地生态系统的影响,因此备受学术界和各国政府关注[5, 8]。森林是陆地生态系统中最大的碳库,其面积约占陆地面积的1/3,对全球碳收支有着重要影响。森林土壤碳储量约占森林生态系统碳储量的2/3,约占全球土壤碳储量的39%[9-10]。森林土壤不仅是植被生长的基础,也是CO2的源、汇地之一,通过土壤呼吸排放到大气中的CO2是大气的重要碳源[11]。在碳中和背景下,被关注的重点是非自然变动引起的森林土壤呼吸的增加或减少,这种变化量才是森林生态系统的有效碳源或碳汇。森林生态系统中的人为干扰(如森林经营活动)能在很大程度上影响土壤CO2排放[12]。其中,森林采伐作为最重要的经营措施及干扰程度最大的人为干扰活动之一,通过改变植被组成、林内光照、凋落物质量、数量及土壤温湿度等进而影响着土壤CO2排放。学者们针对不同气候带的森林开展了多种采伐方式对土壤呼吸影响的研究,但是结论并不一致,存在很大的不确定性。科学认识采伐干扰下森林土壤CO2排放的特征,探讨减少土壤呼吸的森林经营措施对于增强森林的固碳减排功能具有重要的科学意义和实践价值。为此,本研究综述了不同采伐方式对森林土壤呼吸的影响及其机制,主要包括不同采伐方式处理下、不同森林类型对于森林土壤呼吸总量、土壤呼吸组分及其温度敏感性(Q10)的影响,并总结了采伐对土壤呼吸影响的调控因子,在此基础上,提出了该领域的研究前景,以期为中国选择合理的采伐方式,降低森林土壤CO2排放,2060年实现碳中和提供参考。

    土壤释放CO2的过程称为土壤呼吸,包括3个生物学过程和1个非生物学过程[13]。3个生物学过程分别是自养呼吸、土壤微生物异养呼吸和土壤动物异养呼吸。植物根系与根际呼吸产生的CO2排放,称为自养呼吸;微生物分解土壤有机质产生的CO2排放,称为土壤微生物异养呼吸;土壤动物呼吸产生的CO2排放,称为土壤动物异养呼吸[13]。非生物学过程是指土壤含碳矿物质化学氧化产生的CO2排放[13],其产生的CO2量远少于生物学过程而通常被忽略不计。

    土壤呼吸组分因其产生途径、产生部位和所利用碳源的不同有着不同的术语表达,且经常存在土壤呼吸组分术语混用的问题[14]。在分析森林采伐对土壤呼吸的影响时,可以以采伐影响土壤呼吸的产生途径、产生部位和碳源等某一方面为主进行分析。从土壤CO2排放的产生途径来分析,可以分为自养呼吸(autotrophic respiration)和异养呼吸(heterotrophic respiration)[8, 15]。从土壤CO2排放的产生部位来分析,可分为根际区、无根系影响的土壤和凋落物层3个部位[16]。从土壤CO2排放所利用的碳源来分析,可以分为土壤有机质源CO2和植物源CO2(包括凋落物源、死根源、活根源)[17-19]

    森林采伐是一种非常普遍的经营作业方式,一般分为针对成熟林或过熟林的皆伐、择伐和渐伐等主伐、针对中幼龄林的间伐以及针对防护林的更新采伐。皆伐是将伐区上的林木一次性全部伐除或几乎伐除(保留部分母树)的主伐方式。择伐、渐伐、间伐、更新采伐都是仅将伐区上的林木移除一部分,为方便叙述,本研究统一称它们为部分采伐。森林采伐要砍伐林分中的所有或部分林木,势必会降低冠层覆盖,去除林分或改变林分结构,影响各种环境因子,进而影响土壤呼吸。

    目前,关于皆伐影响土壤总呼吸的研究有很多,结果并不一致(表1),可以分为增加、不变、减少3种结论。通常认为皆伐短期内会增加土壤总呼吸[20]:锐齿栎Quercus aliena皆伐4个月后土壤总呼吸增加5%[21];挪威云杉Picea abies林皆伐后第2年土壤总呼吸增加29%,第3年增加52%[22];云杉Picea asperata林皆伐后2 a土壤总呼吸增加50%[23]。其主要原因有:①土壤温度升高提升了异养呼吸速率。林地皆伐后土壤受阳光直射,其温度会发生剧烈的变化[24],从而提升了土壤有机质的分解速率和土壤微生物异养呼吸[21],大量研究表明土壤温度提升可以解释85%~98%的土壤呼吸变化[25-29]。②土壤有机质增加。皆伐林地内残留的死根、凋落物和伐木残留物的丰富和矿化导致土壤呼吸在皆伐后几年内增加[23]。③土壤理化性质变化。皆伐会通过影响土壤理化性质,进而影响土壤呼吸。皆伐影响土壤氮含量,土壤氮能加速植物生长,影响土壤根呼吸,同时土壤氮也是土壤微生物的重要影响因子;皆伐还会影响土壤pH,土壤pH通过调控土壤中化学反应的进程和土壤酶活性来间接影响土壤呼吸[30]。还有研究表明皆伐会影响土壤全碳、全氮、碳氮比、速效氮磷钾和土壤容重等,而这些都是土壤呼吸的影响因子[31-35]

    表 1  土壤呼吸及其组分对皆伐的响应
    Table 1  Response of soil respiration and its components to clear cutting
    地点气候带皆伐更新
    方式
    剩余物
    处理方式
    伐后时
    间/a
    观察时间森林类型总呼吸/
    %
    自养呼吸/
    %
    异养呼吸/
    %
    Q10/%参考文献
    中国福建省 亚热带 1 5 cm以上收集,
    以下归堆清理
    5~6 整年 杉阔混交林 −37 −48 −34 −17 [37]
    中国黑龙江省 温带  1 1 生长季 白桦沼泽 −6 [58]
    中国吉林省 温带  1 主干移除
    枝叶未清
    12~13 生长季 阔叶红松林 −25 −35 [72]
    美国加利福尼亚州 温带  1 1~2 整年 云杉林 −29 [38]
    中国甘肃省 暖温带 1 1 4个月后整年 锐齿栎 5 [21]
    俄罗斯莫斯科州 温带  1 凋落物保留
    剩余物保留
    1~2 生长季 云杉林 50
    50
    [23]
    芬兰 温带  1 保留 1 整年 挪威云杉 −16 16
    17
    25
    [22]
    2 整年 29
    3 整年 52
    美国密苏里州 热带  主干移除 2~4 整年 栎-山核桃林 −18 [73]
    全移除 −17
    芬兰 温带  1 全部移除 1 整年 苏格兰松 23 [64]
    2 整年 −16
    3 整年 −20
    加拿大魁北克省 寒带  1 6~7 整年 黑云杉 16 [74]
    2 9
    加拿大新斯科舍省 温带  3 3~4 整年 混合杉木林 −1 [43]
    中国浙江省 亚热带 1 25~26 整年 杉木林 17 −15 [75]
    瑞典乌普萨拉省 温带  4 树干树桩收获
    树冠枝条保留
    21~22 整年 苏格兰松
    挪威云杉
    −10 [65]
    五大湖流域 温带  1 生长季 糖枫 −7 [57]
    日本 温带  4 保留竹类 1~3 整年 寒温带针阔
    混交林
    17 [76]
    日本 温带  4 保留竹类 1~10 整年 寒温带针阔
    混交林
    61 [77]
    马来西亚 热带  1 树干收获,
    其余保留
    1~2 5个月 重红婆罗双林
    龙脑香林
    不变
    不变
    [44]
    韩国 温带  4 1 整年 红松林 41 [78]
    中国浙江省 亚热带 1 移除 1 整年 杉木林 −15 −20 [79]
    5 火烧 −27 −27
    中国浙江省 亚热带 1 保留 2 整年 杉木林 13 −10 [79]
    1 保留且翻土 32 −11
    中国浙江省 亚热带 1 保留 3 整年 杉木林 16 −10 [79]
    1 保留且翻土 30 −12
    英国英格兰 温带  1 1 整年 云杉 −22 [80]
    2 −42
    3 −30
    4 −10
    马来西亚 热带  1 1~9 隔4周测2周 阔叶混交林 13 [81]
    日本 亚热带 1 清除 2 每年5−10月 天然混交林 16 14 [82]
    3 11 33
    4 20 48
    5 5 57
    6 5 67
    7 20 29
    8 4 38
      说明:皆伐更新方式中1表示皆伐后自然恢复,2表示皆伐后翻土,3表示皆伐后喷洒除草剂,4表示皆伐后人工种植,5表示皆伐后火烧。     栎Quercus spp.,山核桃Carya spp.,黑云杉Picea mariana,重红婆罗双Shorea spp.,龙脑香Dipterocarpus spp.,红松      Pinus koraiensis。空白表示无此项观测记录
    下载: 导出CSV 
    | 显示表格

    也有少数研究认为,皆伐造成的根呼吸降低大于采伐造成的异养呼吸增加,因此皆伐造成土壤总呼吸的降低[36]。杉阔混交林皆伐第5年土壤呼吸减少48%[37]。云杉林皆伐1 a后土壤呼吸减少29%[38]。皆伐减少土壤呼吸的原因主要有:①皆伐后土壤自养呼吸显著下降。根呼吸占土壤呼吸的50%[39],皆伐迹地植被活根的减少会导致土壤自养呼吸速率下降[40],当自养呼吸下降幅度大于异养呼吸的增加幅度时土壤总呼吸速率表现为降低[21]。②皆伐后采伐剩余物的清除方式。皆伐后火烧或清除采伐剩余物、清理凋落物等都会减少土壤有机质输入,从而减少碳输入[41],微生物的异养呼吸会在一段时间后消耗掉大量的土壤碳[42],减少皆伐迹地土壤碳含量,进而降低土壤呼吸。③皆伐迹地植被恢复的时间不同。从皆伐后立即开展研究到皆伐后若干年开展研究,观察到的皆伐迹地恢复阶段不统一,导致相同气候和人为干扰措施可能因为不同植被恢复阶段而得到不同的研究结论。

    还有研究发现皆伐对土壤呼吸无显著影响。例如:杉木Cunninghamia lanceolata林皆伐后第25年土壤呼吸未发生明显变化[43]。杨玉盛等[25]发现杉木林皆伐后土壤呼吸的变化不显著。皆伐后土壤呼吸变化不大的原因可能有:①土壤异养呼吸的增加弥补了根呼吸的减少导致了土壤总呼吸基本不变。皆伐后根系呼吸的下降和物质输入的消失可降低土壤自养呼吸,而采伐剩余物的分解增加及新近死亡的根系分解可能促进土壤异养呼吸,两方面综合作用可能导致土壤总呼吸的不变[25]。也有研究表明,皆伐后土壤微生物呼吸的增加与根呼吸的减少相抵消,从而使得土壤总呼吸未发生明显变化[44-45]。②研究区微地形的影响和地下潜在因素众多,尤其是皆伐后林区排水能力的变化影响地下水位,进一步影响微生物活性,本应增加的土壤微生物呼吸未发生明显变化,导致土壤呼吸未发生明显变化[43]

    综上可见,皆伐对土壤呼吸影响的效果因皆伐措施的不同、森林类型的不同和伐后恢复时间的不同呈现显著的时空和地域异质性[46-47]

    部分采伐对森林土壤呼吸影响的研究相对于皆伐较少[48]。部分采伐收获了部分林木,对林分及其土壤的干扰程度相比皆伐较低。从目前的研究情况(表2)来看,部分采伐对土壤呼吸影响的研究结果也不一致,有增加[49-50],减少[36, 51]和基本不变[52-54]共3类。有关部分采伐对土壤呼吸影响的研究常聚焦于不同采伐强度的影响上。如马尾松Pinus massoniana林间伐15%和间伐70%后1 a内土壤呼吸分别为保持不变和增加17%[55];杉阔混交林间伐35%、49%和68%第5年土壤呼吸分别增加15%、增加16%和减少10%[37]。毛竹Phyllostachys edulis林择伐24%第3~8个月土壤呼吸减少16%[56]。糖枫Acer saccharum林间伐35%第5~10个月土壤呼吸减少19%[57]。白桦Betula platyphylla沼泽林渐伐45%第8~13个月土壤呼吸减少15%[58]

    表 2  土壤呼吸及其组分对部分采伐的响应
    Table 2  Response of soil respiration and its components to partial cutting
    地点气候带部分采伐
    强度/%
    剩余物
    处理方式
    伐后时
    间/a
    观察时间森林类型总呼吸/
    %
    自养呼吸/
    %
    异养呼吸/
    %
    Q10/%参考文献
    中国湖北省 亚热带 除灌 清理 1 整年 马尾松林 −17 −17 −18 [56]
    15 移除树干 −14 11
    70 移除树干 17 11 22
    中国山西省 温带  20 清除 1 生长季 油松人工林 −4 18 −6 6 [69]
    30 23 64 19 −30
    40 52 290 30 −13
    中国湖北省 亚热带 24 1 生长季 毛竹林 −16 28 −29 9 [56]
    中国福建省 亚热带 35 5 cm以上收集,
    以下归堆
    5~6 整年 杉阔混交林 15 14 15 52 [37]
    49 16 13 17 34
    68 −10 −5 −12 −1
    中国黑龙江省 温带  45 1 生长季 白桦沼泽 −15 [58]
    中国陕西省 温带  15 清除采伐剩余物 3~4 生长季 华北落叶松 −5 47 [84]
    35 16 3
    50 −3 15
    中国陕西 温带  12 1~4 生长季 华北落叶松 [28]
    32 17
    47
    斯洛文尼亚 温带  50
    100
    1~3 生长季 山毛榉林 47
    69
    [85]
    中国黑龙江省 温带  20 堆腐 1~4 生长季 针阔混交林 23 [59]
    39 22
    52 24
    62 27
    71 22
    挪威 寒带  41 32~33 夏季 挪威云杉 13 [86]
    55 17
    加拿大安大略省 温带  50 2 生长季 耐寒阔叶林 54 [57]
    爱尔兰 温带  42 1~2 整年 云杉 13 [87]
    日本 亚热带 50 2~4 整年 日本雪松林 46 [88]
    加拿大安大略省 温带  28 1 生长季 杉阔混交林 17 −25 [50]
    2 18 −6
    3 16 19
    中国湖北省 亚热带 23 手工除草为对照
    除草剂除草为处理
    1 整年 毛竹林 −7 20 −13 3 [83]
    斯洛文尼亚 温带  50 1~3 生长季 云杉林/冷杉林 26 [85]
    100 48
    中国湖北省 亚热带 15 清除 1~3 全年 马尾松林 29 14 39 [89]
    70 42 19 59
      说明:日本雪松Cryptomeria japonica。空白表示无此项观测记录;−表示减少
    下载: 导出CSV 
    | 显示表格

    部分采伐增加土壤呼吸的原因有:①部分采伐减小了森林郁闭度,林下光照强度增加导致土壤温度增加,促进土壤有机质分解,从而增加土壤异养呼吸,同时也促进植物根系的生长,增加土壤自养呼吸[49];②部分采伐后采伐剩余物例如木屑和树枝树叶等进入土壤,为土壤微生物活动提供底物,增加土壤异养呼吸[59]。部分采伐降低土壤呼吸可以归因为:①部分采伐时整株植物被移除,凋落物减少,碳底物供应下降导致土壤呼吸减弱[56]。②部分采伐后乔木层减少,树木蒸腾作用减弱,地下水位上升,土壤孔隙减少,导致土壤呼吸减小[58]。部分采伐对土壤呼吸无显著影响可能是因为:①部分采伐提高了土壤异养呼吸,但又同时降低了根呼吸,综合作用下部分采伐对土壤呼吸无影响[55]。②部分采伐后林地凋落物储量、有机碳储量、土壤总孔隙度及细根生物量仍能维持较高的水平,与对照相比土壤呼吸未发生显著变化[37]

    总体上,部分采伐对土壤湿度、细根生物量和土壤碳储量(包括土壤总碳含量、土壤有机碳和微生物量碳)无显著影响。但是部分采伐会导致凋落物等显著减少,土壤温度升高,土壤总呼吸上升。轻度和中度部分采伐显著增加土壤呼吸,尤其是在植被恢复的早期阶段(≤2 a)[60]

    虽然近些年来对土壤呼吸组分的研究大幅度增加(表1表2),但是与采伐对森林土壤呼吸影响的研究相比,采伐对土壤呼吸组分影响的研究要少得多。土壤自养呼吸和土壤异养呼吸受到土壤温度、土壤湿度和细根生物量等一系列因素的影响[44]

    皆伐导致细根大量死亡,土壤自养呼吸显著下降[37]。皆伐后森林乔木层消失,太阳直射地表导致土壤温度升高,地表水分加速蒸发[61]。地表温度的上升促进了枯枝落叶层和表层土壤有机质的分解[29];皆伐带来的新鲜采伐剩余物为土壤微生物提供了大量的碳源[62],以上2点原因导致了皆伐后土壤异养呼吸增加[63]。但此部分碳源分解较快,长时间土壤异养呼吸下降会导致土壤异养呼吸短时间内增加长时间内减少,其他研究也佐证了这一结论。例如苏格兰松Pinus sylvestris皆伐第1年土壤异养呼吸增加23%,第2年减少16%,第3年减少20%[64]。这是因为皆伐时产生的碎木屑进入土壤,增加了土壤微生物呼吸的底物,导致了土壤异养呼吸的增加,但是这部分底物很少,在第2年和第3年时底物分解殆尽,土壤异养呼吸下降。杉阔混交林皆伐第5年土壤自养呼吸减少48%,土壤异养呼吸减少34%[37]。这是因为皆伐收获了林木,植物根大量死亡,土壤自养呼吸显著下降,同时皆伐后林地凋落物、土壤总孔隙度和土壤有机质都出现了明显的下降,土壤异养呼吸显著下降。苏格兰松和挪威云杉在皆伐第22年土壤异养呼吸减少10%[65]。而这可能是因为此研究采用挖掘机收获伐桩,比起用带有刀片的推土机,对土壤的扰动更小,不同收获方式导致土壤呼吸的变化不同。

    总体来看,与对照组相比,皆伐破坏了森林地上植被,导致根系死亡,土壤自养呼吸下降;皆伐后保留采伐剩余物短时间内土壤异养呼吸增加,长时间后则土壤异养呼吸会下降。这是因为保留采伐剩余物为土壤微生物呼吸和土壤动物呼吸提供了碳源,但是这种碳源易分解,短时间内会释放大量CO2,长时间后则易分解有机质减少,土壤异养呼吸下降。同时皆伐砍伐灌木、清除草本和根系分解可能补偿根系和根际呼吸的减少[66]

    部分采伐主要通过以下两方面影响土壤呼吸组分:①不同的采伐剩余物处理方法对土壤微生物底物的供应不同,影响土壤微生物呼吸,从而影响土壤异养呼吸。②部分采伐强度不同,对植物根的破坏程度不同,对土壤自养呼吸的影响也不同。例如,马尾松林间伐15%和70%在1 a内(仅移除树干)土壤自养呼吸分别减少14%和增加11%,土壤异养呼吸分别增加11%和22%。这是因为15%间伐清除了林下灌木和部分林下树种,这些植被细根比例大且分布较浅,清除后可能会显著降低表层土壤根系生物量,导致土壤自养呼吸减少[55];70%间伐导致地上植被减少,但是充足的养分会促进剩余植被的生长,导致根系生物量增加,进而增加根呼吸,原本应减少的根呼吸无显著变化[55];2种强度的采伐后林地残留的伐根死亡为土壤异养呼吸增加了底物,同时活立木的减少改变了林木微环境,为土壤微生物活动创造了适宜的条件,导致土壤异养呼吸增加[67-68]。油松Pinus tabulaeformis人工林择伐20%、30%和40%第2~7个月(采伐剩余物清除)土壤自养呼吸分别增加18%、64%和290%,土壤异养呼吸分别减少6%、增加19%和增加30%[69]。此研究中随着林分密度的递减,林地总活根量密度增大,而总活根量在一定程度上决定根呼吸,故随采伐强度增加,土壤自养呼吸越强。随着采伐强度的增加,进入土壤的枯枝落叶增加,而枯枝落叶层的覆盖对土壤CO2的排放有一定的阻碍[70],故对照组异养呼吸低于处理组。毛竹林间伐24%第3~8个月土壤异养呼吸增加28%,土壤自养呼吸减少29%[56]。这是因为采伐后林地表面温度升高,地上碳供应减少,根基分泌物减少,导致土壤有机碳分解增加,土壤矿质呼吸增加,而根呼吸的下降可能是因为底物供应的下降[71]。杉阔混交林择伐35%、49%和68%第5年(采伐剩余物长度5 cm以上收集以下归堆清理)土壤自养呼吸分别增加14%、增加13%和减少5%,土壤异养呼吸分别增加15%、增加17%和减少12%[37],而这些差异在统计学上并不显著。这是因为择伐后林地凋落物储量、土壤总孔隙度、有机碳储量、有机质和细根生物量仍维持在较高的水平,土壤呼吸组分未发生显著变化。

    可以看出,部分采伐对土壤呼吸组分的影响会随着采伐剩余物处理方式的不同而发生显著的变化,保留采伐剩余物短时间内通常会增加土壤异养呼吸;同时林分根系的生长也会随着伐后恢复的程度而得到增强,伐后恢复时间越久,部分采伐对土壤呼吸组分的影响越小。

    土壤温度是影响土壤呼吸的重要环境因子,土壤呼吸的温度敏感性用Q10来表示,是指土壤呼吸随温度每升高10 ℃所增加的倍数。Q10值不仅随地理位置、森林生态系统的不同而不同,也会受到人为干扰活动如采伐的影响。

    皆伐对土壤呼吸温度敏感性的影响主要取决于皆伐迹地植被恢复的时间。例如欧洲云杉皆伐1~3 a Q10连年上升,第1年增加16%,第2年增加17%,第3年增加25%[22],阔叶红松林皆伐13 a后生长季Q10减少35%[72],但杉木林皆伐1~3 a无论是移除还是保留采伐剩余物Q10皆下降[79]。而杉阔混交林皆伐5 a后Q10减少17%[37]。由于皆伐后采伐剩余物管理方式的不同,进入土壤的易分解有机质有多有少,短期内Q10也表现出不同的变化规律,但长期后因为皆伐迹地植被的恢复,土壤温度敏感性基本呈现下降的趋势。

    部分采伐对土壤温度敏感性的影响主要取决于部分采伐的强度,但是不同研究的结果并不统一。低强度部分采伐下,短时间内Q10通常增加,毛竹林23%间伐1 a后Q10增加3%[83],油松人工林20%间伐1 a后生长季Q10增加6%[69],毛竹林24%间伐1 a后生长季Q10增加9%[56],杉阔混交林35%和49%间伐5~6 a内Q10分别增加52%和34%[37],华北落叶松15%间伐3~4 a内生长季Q10增加47%[84]。但是也有结果相反的研究,例如杉阔混交林28%间伐1、2 a后Q10分别减少25%和6%,这和此研究中夏季降雨量减小有关。高强度采伐后Q10的变化并不统一,例如油松人工林40%间伐1 a后生长季Q10减少13%[69],杉阔混交林间伐68% 5~6 a内Q10减少1%[37],华北落叶松50%间伐3~4 a Q10增加15%[74]。这可能是因为高强度部分采伐后林窗面积增大,其他植物荫蔽林窗的能力受到当地气候等因素的影响。从以上研究中可以看出,一部分研究结果呈现轻度、中度部分采伐短时间内Q10增加的趋势,随着植被的恢复,Q10也逐渐接近对照林。但是也有部分研究受到其它因素例如降雨量变化的影响,结果与上述研究相反。

    总体上皆伐会破坏森林植被,造成植物根系大量死亡,土壤自养呼吸降低,同时皆伐将更多的枯枝落叶带入土壤,加上死亡的植物根系,土壤异养呼吸增加。两者共同作用决定了土壤总呼吸的变化,如果皆伐后对皆伐迹地进行清理,土壤总呼吸往往会下降,如果皆伐迹地内采伐剩余物较多,土壤总呼吸可能会先上升后下降。与皆伐相比,部分采伐对森林的干扰程度不同,一定强度的部分采伐可能会增加土壤总呼吸,随着部分采伐强度的增大,土壤呼吸的变化接近皆伐迹地内土壤呼吸的变化。

    森林土壤呼吸是陆地生态系统碳循环的重要组成部分,在全球气候变化中起着重要的作用。皆伐或部分采伐作为重要的人为干扰经营措施,对森林林冠、覆盖率、枝叶雨水截流、土壤温度、土壤湿度等土壤理化性质和土壤呼吸有着显著的影响。森林不同强度部分采伐对伐后植被不同恢复阶段土壤呼吸和土壤碳储量的影响尚不清晰,建议加强土壤呼吸组分对部分采伐强度响应的长期研究。除此之外,森林采伐和林下除灌、除草、定期打枝等其他经营措施的交互作用以及全球大气CO2浓度上升等全球变化因子对区域森林变化也应纳入考量中。

  • 图  1  P波段SAR数据(A)及其覆盖区样点分布(B)、林分概况(C)

    Figure  1  P band SAR data (A), the distribution of samples (B) and stand examples (C)

    图  2  LiDAR衍生数据

    Figure  2  Lidar derived data

    图  3  P波段森林散射机制分析

    Figure  3  Forest scattering mechanism analysis at P band

    图  4  4种模型森林AGB估测与LiDAR抽样点散点图

    Figure  4  Scatter plots of estimated and LiDAR forest AGB of four models

    图  5  森林AGB估测结果与LiDAR抽样点差值直

    Figure  5  Histograms of difference between forest AGB estimation results and LiDAR sampling points

    表  1  P波段SAR极化特征与森林AGB相关性分析

    Table  1.   Correlation coefficients of polarization feature and forest AGB

    特征RP特征RP特征RP
    HH 0.401** 0 F_VOL −0.129 0.173 alpha −0.734** 0
    HV −0.249** 0.008 Y_DBL 0.719** 0 beta 0.037 0.696
    VV −0.322** 0 Y_ODD 0.323** 0 delta 0.221* 0.018
    FD1 −0.453** 0 Y_VOL −0.213* 0.024 gamma 0.314** 0.001
    FD2 −0.448** 0 Y_HLX 0.187* 0.047 PDR 0.629** 0
    F_DBL 0.650** 0 entropy −0.136 0.151 RVI −0.374** 0
    F_ODD 0.303** 0.001 anisotropy 0.436** 0
      说明:R代表Pearson相关系数;P代表显著性水平,**代表0.01水平上显著;*代表0.05水平上显著
    下载: 导出CSV

    表  2  基于4种模型的P波段SAR 森林AGB估测结果

    Table  2.   P band SAR forest AGB inversion using four models

    波段反演模型模型参数R2RMSEAcc/%MAE
    P MLR 显著性为0.0008 0.43 19.16 63.55 16.99
    KNN K值为11,欧氏距离 0.54 17.04 71.18 13.44
    SVR 惩罚系数7.0 0.54 17.09 71.15 13.45
    RF 决策树数为100,树深度为11 0.60 15.98 72.97 12.60
    下载: 导出CSV

    表  3  基于4种模型的P波段SAR 森林AGB反演情况     

    Table  3.   P band SAR forest AGB inversion based on four models

    序号AGB分段/
    (t·hm−2)
    R2RMSE/
    (t·hm−2)
    Acc/%
    10~300.554.9276.12
    30~600.137.9581.85
    >600.2514.8981.57
    20~500.478.4672.15
    >500.3015.8078.19
    30~400.486.5274.36
    >400.3616.2174.82
    下载: 导出CSV
  • [1] CARTUS O, SANTORO M, WEGMÜLLER U, et al. Benchmarking the retrieval of biomass in boreal forests using P band SAR backscatter with multi-temporal C- and L-band observations[J/OL]. Remote Sensing, 2019, 11(14): 1695 [2022-01-01]. doi: 10.3390/rs11141695.
    [2] SANTI E, PALOSCIA S, PETTINATO S, et al. Machine-learning applications for the retrieval of forest biomass from airborne P band SAR data[J/OL]. Remote Sensing, 2020, 12(5): 804 [2022-01-05]. doi: 10.3390/rs12050804.
    [3] LIAO Zhanmang, HE Binbin, QUAN Xingwen. Potential of texture from SAR tomographic images for forest aboveground biomass estimation[J/OL]. Int J Appl Earth Obs Geoinf, 2020, 88: 102049 [2022-01-05]. doi: 10.1016/j.jag.2020.102049.
    [4] QUAN Xingwen, HE Binbin, YEBRA M, et al. A radiative transfer model-based method for the estimation of grassland aboveground biomass [J]. Int J Appl Earth Obs Geoinf, 2017, 54: 159 − 168.
    [5] JI Yongjie, XU Kunpeng, ZENG Peng, et al. GA-SVR algorithm for improving forest above ground biomass estimation using SAR data [J]. IEEE J Selected Top Appl Earth Obs Remote Sensing, 2021, 14: 6585 − 6595.
    [6] LIU Y Y, DIJK A, RAMD J, et al. Recent reversal in loss of global terrestrial biomass [J]. Nat Clim Change, 2015, 5(5): 470 − 474.
    [7] BOUDREAU J, NELSON R F, MARGOLIS H A, et al. Regional aboveground forest biomass using airborne and spaceborne LiDAR in Québec [J]. Remote Sensing Environ, 2008, 112(10): 3876 − 3890.
    [8] JI Yongjie, HUANG Jimao, JU Yilin, et al. Forest structure dependency analysis of L-band SAR backscatter[J/OL]. PeerJ, 2020, 8: e10055 [2022-01-05]. doi: 10.7717/PEERJ.10055.
    [9] 李增元, 赵磊, 李堃, 等. 合成孔径雷达森林资源监测技术研究综述[J]. 南京信息工程大学学报(自然科学版), 2020, 12(2): 150 − 158.

    LI Zengyuan, ZHAO Lei, LI Kun, et al. A survey of developments on forest resources monitoring technology of synthetic aperture radar [J]. J Nanjing Univ Inf Sci Technol Nat Sci Ed, 2020, 12(2): 150 − 158.
    [10] 冯琦, 陈尔学, 李增元, 等. 基于机载P波段全极化SAR数据的复杂地形森林地上生物量估测方法[J]. 林业科学, 2016, 52(3): 10 − 22.

    FENG Qi, CHEN Erxue, LI Zengyuan, et al. Forest above-ground biomass estimation method for rugged terrain based on airborne P-band PolSAR data [J]. Sci Silv Sin, 2016, 52(3): 10 − 22.
    [11] DOBSON M C, ULABY F T, LETOAN T, et al. Dependence of radar backscatter on coniferous forest biomass [J]. IEEE Trans Geosci Remote Sensing, 1992, 30(2): 412 − 415.
    [12] TOAN T L, BEAUDOIN A, RIOM J, et al. Relating forest biomass to SAR data [J]. IEEE Trans Geosci Remote Sensing, 1992, 30(2): 403 − 411.
    [13] TOAN T L, QUEGAN S, WOODWARD I, et al. Relating radar remote sensing of biomass to modelling of forest carbon budgets [J]. Clim Change, 2004, 67(2): 379 − 402.
    [14] LIAO Zhanmang, HE Binbin, QUAN Xingwen, et al. Biomass estimation in dense tropical forest using multiple information from single-baseline P band PolInSAR data [J]. Remote Sensing Environ, 2019, 221: 489 − 507.
    [15] 黄国满. 机载多波段多极化干涉SAR测图系统——CASMSAR[J]. 测绘科学, 2014, 39(8): 111 − 115.

    HUANG Guoman. An airborne interferometric SAR mapping system with multi-band and multi-polarization: CASMSAR [J]. Sci Surv Mapp, 2014, 39(8): 111 − 115.
    [16] ZHAO Lei, CHEN Erxue, LI Zengyuan, et al. Three-step semi-empirical radiometric terrain correction approach for PolSAR data applied to forested areas[J/OL]. Remote Sensing, 2017, 9(3): 269 [2022-01-05]. doi: 10.3390/rs9030269.
    [17] 穆喜云. 森林地上生物量遥感估测方法研究[D]. 呼和浩特: 内蒙古农业大学, 2015.

    MU Xiyun. A Study on the Estimating Method of Forest Above Ground Biomass Based on Remote Sensing Data[D]. Huhhot: Inner Mongolia Agricultural University, 2015.
    [18] 陈传国. 东北主要林木生物量手册[M]. 北京: 中国林业出版社, 1989.

    CHEN Chuanguo. Woody Biomass Manual of Typical Species in the Northeast of China[M]. Beijing: China Forestry Publishing House, 1989.
    [19] CLOUDE S R, POTTIER E. A review of target decomposition theorems in radar polarimetry [J]. IEEE Trans Geosci Remote Sensing, 1996, 34(2): 498 − 518.
    [20] 张王菲, 姬永杰. 极化与干涉SAR植被参数反演[M]. 北京: 中国林业出版社, 2019.

    ZHANG Wangfei, JI Yongjie. Retrieval of Vegetation Parameters from Polarimetric and Interferometric SAR[M]. Beijing: China Forestry Publishing House, 2019.
    [21] 姬永杰. 多频极化SAR技术森林地上生物量反演研究[D]. 昆明: 西南林业大学, 2021.

    JI Yongjie. Retrieval of Forest above Ground Biomass Using Multi-frequency Polarimetric SAR[D]. Kunming: Southwest Forestry University, 2021.
    [22] 郭颖. 森林地上生物量的非参数化遥感估测方法优化[D]. 北京: 中国林业科学研究院, 2011.

    GUO Ying. Optimum Non-Parametric Method for Forest Above Ground Biomass Estimation Based on Remote Sensing Data[D]. Beijing: Chinese Academy of Forestry, 2011.
    [23] TIAN Xin, SU Zhongbo, CHEN Erxue, et al. Estimation of forest above-ground biomass using multi-parameter remote sensing data over a cold and arid area [J]. Int J Appl Earth Obs Geoinf, 2011, 17: 102 − 110.
    [24] LU Dengsheng, CHEN Qi, WANG Guangxing, et al. A survey of remote sensing-based aboveground biomass estimation methods in forest ecosystems [J]. Int J Digital Earth, 2016, 9(1): 63 − 105.
    [25] 韩宗涛, 江洪, 王威, 等. 基于多源遥感的森林地上生物量KNN-FIFS估测[J]. 林业科学, 2018, 54(9): 70 − 79.

    HAN Zongtao, JIANG Hong, WANG Wei, et al. Forest above-ground biomass estimation using KNN-FIFS method based on multi-source remote sensing data [J]. Sci Silv Sin, 2018, 54(9): 70 − 79.
    [26] 方匡南, 吴见彬, 朱建平, 等. 随机森林方法研究综述[J]. 统计与信息论坛, 2011, 26(3): 32 − 38.

    FANG Kuangnan, WU Jianbin, ZHU Jianping, et al. A review of technologies on random forests [J]. Stat Inf Forum, 2011, 26(3): 32 − 38.
    [27] 魏晶昱, 范文义, 于颖, 等. GF-3全极化SAR数据极化分解估算人工林冠层生物量[J]. 林业科学, 2020, 56(9): 174 − 183.

    WEI Jingyu, FAN Wenyi, YU Ying, et al. Polarimetric decomposition parameters for artificial forest canopy biomass estimation using GF-3 fully polarimetric SAR data [J]. Sci Silv Sin, 2020, 56(9): 174 − 183.
    [28] SAATCHI S, MARLIER M, CHAZDON R L, et al. Impact of spatial variability of tropical forest structure on radar estimation of aboveground biomass [J]. Remote Sensing Environ, 2011, 115(11): 2836 − 2849.
    [29] TOAN T L, QUEGAN S, DAVIDSON M W J, et al. The BIOMASS mission: mapping global forest biomass to better understand the terrestrial carbon cycle [J]. Remote Sensing Environ, 2011, 115(11): 2850 − 2860.
    [30] NEUMANN M, SAATCHI S S, ULANDER L M H, et al. Assessing performance of L and P band polarimetric interferometric SAR data in estimating boreal forest above-ground biomass [J]. IEEE Trans Geosci Remote Sensing, 2012, 50(3): 714 − 726.
    [31] KASISCHKE E S, CHRISTENSEN N L, BOURGEAU-CHAVEZ L L. Correlating radar backscatter with components of biomass in loblolly pine forests [J]. IEEE Trans Geosci Remote Sensing, 1995, 33(3): 643 − 659.
    [32] RODRÍGUEZ-VEIGA P, QUEGAN S, CARREIRAS J, et al. Forest biomass retrieval approaches from earth observation in different biomes [J]. Int J Appl Earth Obs Geoinf, 2019, 77: 53 − 68.
    [33] GOLSHANI P, MAGHSOUDI Y, SOHRABI H. Relating ALOS-2 PALSAR-2 parameters to biomass and structure of temperate broadleaf Hyrcanian forests [J]. J Ind Soc Remote Sensing, 2019, 47(5): 749 − 761.
    [34] ENGLHART S, KEUCK V, SIEGERT F. Modeling aboveground biomass in tropical forests using multi-frequency SAR data: a comparison of methods [J]. IEEE J Selected Top Appl Earth Obs Remote Sensing, 2012, 5(1): 298 − 306.
  • [1] 卢佶, 张国威, 吴昊.  基于多时相光学和雷达遥感的太平湖生态保护区森林地上生物量反演 . 浙江农林大学学报, 2023, 40(5): 1082-1092. doi: 10.11833/j.issn.2095-0756.20220682
    [2] 黄屹杰, 张加龙, 胡耀鹏, 程滔.  高山松地上生物量遥感估算的不确定性分析 . 浙江农林大学学报, 2022, 39(3): 531-539. doi: 10.11833/j.issn.2095-0756.20210473
    [3] 杨绍钦, 王翔, 许澄, 商天其.  基于MODIS时间序列数据的竹林地上生物量估算 . 浙江农林大学学报, 2022, 39(4): 734-741. doi: 10.11833/j.issn.2095-0756.20210431
    [4] 卢腾飞, 周律, 胡中岳, 欧光龙, 胥辉.  基于混合效应模型的迪庆云冷杉林地上生物量遥感估测 . 浙江农林大学学报, 2021, 38(3): 510-522. doi: 10.11833/j.issn.2095-0756.20200327
    [5] 庾露, 黄艳霞, 刘警鉴, 段炼.  降雨影响下的水稻长势极化特征提取 . 浙江农林大学学报, 2020, 37(5): 992-998. doi: 10.11833/j.issn.2095-0756.20190605
    [6] 谢福明, 字李, 舒清态.  基于优化k-NN模型的高山松地上生物量遥感估测 . 浙江农林大学学报, 2019, 36(3): 515-523. doi: 10.11833/j.issn.2095-0756.2019.03.012
    [7] 王艳芳, 刘领, 悦飞雪, 李冬, 上官周平.  退耕还林工程对河南省森林地上碳储量的影响 . 浙江农林大学学报, 2019, 36(3): 507-514. doi: 10.11833/j.issn.2095-0756.2019.03.011
    [8] 吕常笑, 邓华锋, 王少杰, 陈振雄, 王雪军.  马尾松不同区域相容性立木材积和地上生物量模型 . 浙江农林大学学报, 2016, 33(5): 790-797. doi: 10.11833/j.issn.2095-0756.2016.05.010
    [9] 冉啟香, 邓华锋, 黄国胜, 王雪军, 陈振雄.  云南松地上生物量模型研究 . 浙江农林大学学报, 2016, 33(4): 605-611. doi: 10.11833/j.issn.2095-0756.2016.04.008
    [10] 王金亮, 程鹏飞, 徐申, 王小花, 程峰.  基于遥感信息模型的香格里拉森林生物量估算 . 浙江农林大学学报, 2013, 30(3): 325-329. doi: 10.11833/j.issn.2095-0756.2013.03.003
    [11] 商珍珍, 周国模, 杜华强.  毛竹林地上生物量与胸径的分形关系 . 浙江农林大学学报, 2013, 30(3): 319-324. doi: 10.11833/j.issn.2095-0756.2013.03.002
    [12] 季碧勇, 陶吉兴, 张国江, 杜群, 姚鸿文, 徐军.  高精度保证下的浙江省森林植被生物量评估 . 浙江农林大学学报, 2012, 29(3): 328-334. doi: 10.11833/j.issn.2095-0756.2012.03.002
    [13] 王晓宁, 徐天蜀, 李毅.  利用ALOS PALSAR双极化数据估测山区森林蓄积量模型 . 浙江农林大学学报, 2012, 29(5): 667-670. doi: 10.11833/j.issn.2095-0756.2012.05.005
    [14] 苏文会, 范少辉, 刘亚迪, 彭颖, 封焕英.  车筒竹地上生物量分配格局及秆形特征 . 浙江农林大学学报, 2011, 28(5): 735-740. doi: 10.11833/j.issn.2095-0756.2011.05.008
    [15] 杨同辉, 达良俊, 李修鹏.  浙江天童国家森林公园常绿阔叶林生物量研究(Ⅱ)群落生物量及其分配规律 . 浙江农林大学学报, 2007, 24(4): 389-395.
    [16] 杨同辉, 达良俊, 宋永昌, 杨永川, 王良衍.  浙江天童国家森林公园常绿阔叶林生物量研究(Ⅰ)群落结构及主要组成树种生物量特征 . 浙江农林大学学报, 2005, 22(4): 363-369.
    [17] 林新春, 方伟, 俞建新, 余学军, 胡超宗, 周林.  苦竹各器官生物量模型 . 浙江农林大学学报, 2004, 21(2): 168-171.
    [18] 姜培坤, 徐秋芳, 俞益武.  土壤微生物量碳作为林地土壤肥力指标 . 浙江农林大学学报, 2002, 19(1): 17-19.
    [19] 贺东北, 骆期邦, 曾伟生.  立木生物量线性联立模型研究 . 浙江农林大学学报, 1998, 15(3): 298-303.
    [20] 斯金平, 姚荣明, 陈德标, 吴长辉.  厚朴人工林生物量的研究 . 浙江农林大学学报, 1993, 10(2): 162-168.
  • 期刊类型引用(1)

    1. 吴宁, 肖瑞, 许艳萍, 杜官本, 李晓平, 孙飞. 生长期和植株性别对工业大麻秆“三大素”的影响. 浙江农林大学学报. 2015(05): 776-782 . 本站查看

    其他类型引用(0)

  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20220111

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2022/5/971

图(5) / 表(3)
计量
  • 文章访问数:  927
  • HTML全文浏览量:  131
  • PDF下载量:  109
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-01-10
  • 修回日期:  2022-04-25
  • 录用日期:  2022-07-08
  • 网络出版日期:  2022-09-22
  • 刊出日期:  2022-10-20

基于机载P波段全极化SAR数据的森林地上生物量估测

doi: 10.11833/j.issn.2095-0756.20220111
    基金项目:  国家自然科学基金资助项目(32160365,31860240,42161059);云南省万人计划青年拔尖人才项目(80201444)
    作者简介:

    姬永杰(ORCID: 0000-0001-8012-4115),讲师,博士,从事SAR、LiDAR数据资源环境反演研究。E-mail: jiyongjie@live.cn

    通信作者: 张王菲(ORCID: 0000-0002-2147-5246),教授,博士,从事农林业遥感研究。E-mail: mewhff@163.com
  • 中图分类号: S757

摘要:   目的  森林生物量的空间精准量化对了解陆地碳储量、碳收支、碳平衡,以及揭示森林碳储量与全球气候变化的影响过程具有重要意义。P波段波长较长,在森林中具有更高的穿透能力,研究机载P波段SAR数据提高森林地上生物量(AGB)估测精度的可行性。  方法  以机载P波段全极化合成孔径雷达(SAR)数据和高精度激光雷达(LiDAR)数据估测的森林AGB抽样点为基础,提取20个极化SAR特征,并分别与森林AGB变化作敏感性响应情况分析。采用多元线性回归模型(MLR)、K近邻方法(KNN)、支持向量回归(SVR)和随机森林(RF)4种估测方法,探究机载P波段SAR数据的森林AGB估测精度。  结果  在较低森林AGB(均值约45 t·hm−2)的森林覆盖区中,P波段的同极化后向散射系数、Freeman-Durden和Yamaguchi分解中的表面和二次散射分量对森林AGB变化敏感;此外H-A-ALPHA极化分解的散射角(alpha)、拓展极化参数极化辨别率参数(PDR)也对森林AGB变化敏感。4种方法估测的森林AGB相对误差均约30%,其中MLR估测结果精度最低,估测精度为63.55%,均方根误差(RMSE)为19.16 t·hm−2;RF估测结果精度最高,估测精度为72.97%,RMSE为15.98 t·hm−2;KNN和SVR估计结果差别不明显,RMSE分别为17.04和17.09 t·hm−2  结论  P波段SAR数据对估测森林AGB具有一定潜力,3种非参数方法的估测结果明显优于MLR参数方法。此外,P波段的森林AGB估测精度受到待估森林AGB水平高低的影响明显,在森林AGB水平较高的分组中估测精度较高。在森林AGB均值为45 t·hm−2,最大值为120 t·hm−2的森林覆盖区,以50 t·hm−2将森林AGB样点分为2组时,高森林AGB组的估测精度高出低AGB组约6%。图5表3参34

English Abstract

朱文见, 张慧, 王懿祥. 采伐对森林土壤呼吸影响的研究进展[J]. 浙江农林大学学报, 2021, 38(5): 1000-1011. DOI: 10.11833/j.issn.2095-0756.20210365
引用本文: 姬永杰, 杨丛瑞, 张王菲, 等. 基于机载P波段全极化SAR数据的森林地上生物量估测[J]. 浙江农林大学学报, 2022, 39(5): 971-980. DOI: 10.11833/j.issn.2095-0756.20220111
ZHU Wenjian, ZHANG Hui, WANG Yixiang. Research progress on effects of cutting on forest soil respiration[J]. Journal of Zhejiang A&F University, 2021, 38(5): 1000-1011. DOI: 10.11833/j.issn.2095-0756.20210365
Citation: JI Yongjie, YANG Congrui, ZHANG Wangfei, et al. Forest above ground biomass estimation using airborne P band polarimetric SAR data[J]. Journal of Zhejiang A&F University, 2022, 39(5): 971-980. DOI: 10.11833/j.issn.2095-0756.20220111
  • 森林在陆地生态系统中发挥着重要的碳汇作用。森林生物量的空间精准量化对了解陆地碳储量、碳收支、碳平衡,以及它们与全球气候变化之间的关系具有重要意义[1-3]。森林地上生物量(AGB)为林木干、枝、叶的有机体干质量,不含林木根部、灌木、草本及林下枯枝落叶,是森林生物量的主体部分。近年来,采用遥感技术反演森林AGB展露出巨大的潜力,光学、激光雷达(LiDAR)、被动微波雷达和合成孔径雷达(SAR)等多种遥感数据用于森林AGB的估测,使得森林AGB估测的区域尺度性、经济性及精度等不断提升[3-7]。其中,SAR不受天时、天气的影响,且微波波长较长,在森林中具有较强的穿透能力,因此在获取森林密度、树高、AGB等森林垂直结构因子方面极具潜力[2, 8-10]。SAR估测森林AGB的能力有赖于其工作频率的高低,大量研究证实了低频波段后向散射对森林AGB变化的敏感性更强。在常用于森林监测的微波波段中,频率较低的P波段SAR后向散射系数对森林AGB变化最为敏感[11-13]。目前P波段的研究多基于机载数据展开,且研究区多位于国外,而国内相关研究则开展较少。CARTUS等[1]基于AfriSAR、BioSAR和TropiSAR机载飞行试验,使用P、C、L波段联合进行了森林AGB反演,结果表明P波段与森林结构因子的树干、树枝的相关性较高,P波段的加入有力地提高了反演的准确度和精度。LIAO等[14]使用TropiSAR机载P波段数据,采用层析方法(TomoSAR)将相干幅度、干涉相位和后向散射特征建模对法属圭亚那热带雨林AGB进行了反演,结果表明树高特征的引入可有效提高森林AGB反演精度,决定系数(R2)最高可达0.7。冯琦等[10]使用国产机载P波段数据结合坡度因子,在考虑当地入射角和坡度的情况下建立对数统计模型,对内蒙古根河市生态站的寒温带针叶林进行森林AGB的反演,最高反演精度R2为0.634、均方根误差(RMSE)为12.07 t·hm−2

    现有采用P波段SAR数据进行的森林AGB估测,多集中在采用P波段的后向散射信息、相位和相干性信息,而对于极化信息的利用则较少。由于SAR信息的差异,使得其用于森林AGB估测的方法也差异明显,如采用后向散射信息估测森林地上生物量,多采用线性回归参数模型;而采用相位和相干性信息利用层析技术进行森林AGB估测则多基于植被微波散射模型或电磁波信号模型进行反演。线性回归模型简单灵活,但通常无法表征P波段特征与森林AGB变化之间的复杂关系;植被散射模型或者电磁波信号模型能够体现森林与P波段电磁波之间的部分物理作用机制,但模型较复杂,应用推广困难。

    近年来,随着极化SAR数据的丰富,可提取的极化SAR特征涌现,非参数模型被广泛应用于SAR极化特征农作物生长参数的定量反演,并表现出较强的反演能力和较好的应用推广性。鉴于P在森林监测中的潜力,其极化特征在森林AGB估测中并未深入探索,参数模型过于简单、植被微波散射模型理解过于困难,本研究以中国北方典型寒温带森林作为研究对象,使用机载P波段SAR数据,提取多种极化特征,在分析其P波段极化散射特征的基础上,探索采用参数和非参数模型进行森林AGB估测的可行性,旨在明确P波段森林的极化散射特征,探索采用P波段极化SAR数据进行森林AGB估测具体应用的有效方法。

    • 研究区为位于内蒙古呼伦贝尔盟的根河市大兴安岭森林生态系统国家野外科学观测研究站(大兴安岭生态站),50°49′~50°51′N,121°30′~121°31′E。该研究站是中国纬度最高的森林生态研究站,面积为102 km2。研究区气候类型为典型的寒温带大陆性季风气候。研究区地势相对平缓,区域内80%的坡度小于15°,海拔高度分布在800~1 200 m,森林覆盖率大于75%,主要树种为兴安落叶松Larix gmelinii、白桦Betula platyphylla、樟子松Pinus sylvestris var. mongolica等。

    • 采用的P波段SAR数据由机载CAMSAR系统获得[15]。通过协议方式获取根河实验区机载P波段全极化(HH/HV/VV/VH) SLC数据。该数据是2013年9月以“奖状Ⅱ”飞机平台,以CASMSAR系统右视观测获取的全极化SAR数据。获取过程中飞行方向为自西向东,飞行高度为5 807 m,获取时间为2013年9月13—16日,极化方式为HH/HV/VH/VV,产品模式为SLC,幅宽为6 km × 7 km,中心入射角为55.058°,距离向分辨率为0.666 m,方位向分辨率为0.625 m。机载P波段SAR数据的预处理关键步骤包括多视处理、正射校正、入射角校正和极化方位角校正。多视视数在距离向和方位向均为3,数据的正射校正、入射校正参考文献[10];极化方位角校正过程和算法见文献[16]。P波段SAR数据预处理结果、森林AGB抽样点及林分概况见图1

      图  1  P波段SAR数据(A)及其覆盖区样点分布(B)、林分概况(C)

      Figure 1.  P band SAR data (A), the distribution of samples (B) and stand examples (C)

    • LiDAR获取的数字表面模型(DSM)、数字高程模型(DEM)用于P波段SAR数据的地理编码,由冠层高度模型(CHM)获取的LiDAR森林AGB数据用于反演模型的训练及验证。本研究获取的机载LiDAR数据是将Leica机载雷达系统荷载于“运-5”飞机平台上,于2012年8—9月在根河实验区开展飞行任务。该原始数据密度为5.6个·m−2的点云数据,激光中心波谱值为1550 nm,在初始点云数据的基础上,提取了研究区域高精度的DEM (图2A)、CHM (图2B)和森林AGB(图2C)等衍生产品。DSM、DEM、CHM数据的详细生成方法参考文献[17],高精度LiDAR森林AGB的详细提取过程与方法参考文献[10,18]。

      图  2  LiDAR衍生数据

      Figure 2.  Lidar derived data

      为保证建模样本和验证样本能够代表整个研究区的森林AGB水平,以高精度LiDAR森林AGB图为基础,按照750 m的空间采样间隔,在ArcGIS中采用交互人工干预(去除道路及裸露地物)的方法选取113个样点(图1B)作为森林AGB反演模型的训练与验证。机载P波段SAR数据覆盖区森林林相不均匀,平均AGB较低,约46.7 t·hm−2,且大于100 t·hm−2的采样点仅有5个。113个样点AGB以10 t·hm−2间隔得分布情况(图1B):0~10 t·hm−2,3个;>10~20 t·hm−2,14个;>20~30 t·hm−2,17个;>30~40 t·hm−2,17个;>40~50 t·hm−2,18个;>50~60 t·hm−2,14个;>60~70 t·hm−2,10个;>70~80 t·hm−2,8个;>80~90 t·hm−2,7个;>90 t·hm−2,5个。

    • 极化目标分解方法是从全极化数据中提取地物极化信息的有效方法,目前多种极化分解参数在森林类型识别中表现出巨大的潜力。本研究采用目前常用的极化后向散射系数、常用的3种极化分解方法提取极化SAR特征,并基于此分析森林的极化散射特征。提取P波段HH、HV和VV等3个极化的后向散射系数,基于3个后向散射系数的雷达植被指数(RVI)、极化辨别率参数(PDR);基于Freeman-Durden三分量分解的体散射分量(FVOL)、单次散射分量(FODD)、二次散射分量(FDBL)、体地散射比分量(FD1/FD2,1表示地散射分量是ODD和DBL的和,2表示地散射分量仅为ODD);基于Yamaguchi的体散射分量(YVOL)、单次散射分量(YODD)、二次散射分量(YDBL)、螺旋体散射分量(YHLX);基于H-A-ALPHA极化分解的极化散射熵(entropy)、反熵(anisotropy)、散射角(alpha)、目标方位向角(beta)、相位差角1(gamma)、相位差角2(delta)。3种极化分解方法参考文献[19-20]。由于本次飞行试验时并未布设角反射器用于定标,因此本研究使用的后向散射系数值仅有相对含义。

      森林在P波段极化特征响应分析的目的是为了确定P波段对森林AGB动态变化敏感的极化特征参数,从而确定有效的极化特征参数进行森林AGB的估测。将研究区的森林AGB划分为A (表示生物量在0~30 t·hm−2变化时对应后向散射系数的箱线变化,均值约20 t·hm−2),B (31~50 t·hm−2,40 t·hm−2),C (51~70 t·hm−2,60 t·hm−2),D (71~90 t·hm−2,80 t·hm−2),E (>91 t·hm−2,100 t·hm−2)等5个变化等级,分别制作各等级相应P波段SAR提取参数值的箱线图,分析研究区各极化特征参数对森林AGB动态变化的响应,进而分析其极化散射特征。为了定量的分析各极化特征与森林AGB变化的关系,计算了它们之间的皮尔逊相关系数及显著性水平。

    • 多元线性逐步回归模型(MLR)是森林AGB估测中最常用、最经典的方法之一。与常规的线性回归模型相比,MLR可同时完成模型输入参数的优选,进而提高森林AGB估测的效率和精度。MLR是将自变量逐个引入,每次判断自变量对因变量影响的显著性,并对模型中的自变量进行检验,逐个从模型中剔除不显著的变量,从而得到最优模型。即在保证显著性值在0.05以下的情况下,筛选相关性高的特征变量,进而得到因变量的最优估计。MLR的实现算法详见文献[21]。

    • KNN[22-25]、SVR[22, 24]和RF[24, 26]是森林AGB估测中常用的非参数模型,与参数模型相比,无固定的模型结构,通常通过数据驱动的方法来确定模型结构并用于森林AGB的估测,因此也称为机器学习方法。在森林AGB估测中,这3种方法各有优势,因此选取这3种方法来探究非参数方法在P波段森林AGB估测中的潜力。

    • 反演结果精度的定量评价通过反演结果与真值之间的决定系数(R2)、均方根误差(RMSE)、平均绝对误差(MAE)、估测精度(Acc)来表征。这4个参数的计算公式参考文献[27]。

    • 本研究获取的P波段中心波长为64 cm,该波段在森林中的穿透性较强。由图3可知:P波段20个SAR特征均表现出对研究区森林AGB变化的敏感性。P波段3种极化方式后向散射系数对森林AGB在0~110 t·hm−2内动态变化的响应与已有研究[28]相似,即后向散射系数的值均随着森林AGB的升高而缓慢增长。在各森林AGB水平中,HH、HV和VV散射能量强度相差并不明显,并且从箱线图(图3)的中值可以看出:HH的单调增长趋势较其他2个极化明显,且异常值较少。这可能是由于P波段波长较长,比较粗树干等圆柱体形状的散射体为森林中的散射体,使得去极化特征明显降低,因此同极化的HH和VV后向散射能量对AGB变化的敏感性明显高于HV[29]。此外,由后向散射系数计算的RVI也随着森林AGB增加而单调增加且未出现饱和现象,并且通过箱线图的宽度变化可以看出:各水平RVI箱宽变化不明显,说明该值对森林结构的变化不敏感,而对森林AGB的变化比较敏感。相比RVI,PDR非单调增长,在不同森林AGB水平,占主导散射的散射体会发生明显变化[30-31],该现象也可由Freeman-Durden中各个参数随森林AGB变化的响应加以说明。P波段Freeman-Durden分解相关特征中,除FODD和FDBL外,均随着森林AGB的增加而增加,并且2种体地散射比的值均远远大于1。Yamaguchi分解特征对森林AGB的响应趋势与Freeman-Durden分解相关特征的趋势基本一致。P波段H-A-ALPHA分解提取的参数中,散射熵(entropy)随着森林AGB增加而单调增加,而alpha则敏感性不明显。但beta、gamma和delta对森林AGB变化的敏感性则高于散射角alpha。P波段森林的极化散射特征随森林AGB的变化与已有研究中短波长的研究差异明显,说明现有的极化分解建模方法与该森林AGB水平P波段的散射特征有明显差异。

      图  3  P波段森林散射机制分析

      Figure 3.  Forest scattering mechanism analysis at P band

      表1可知:极化分解中表征体散射分量的极化特征与森林AGB变化的显著性较低或不存在显著性,而二次散射特征分量则与森林AGB变化具有强相关性且显著性水平最高。这说明在本研究区,树干与地表形成的二次散射对森林AGB的变化最为敏感,且目前极化分解中的体散模型并不适合P波段研究区森林AGB水平下森林的散射机制。表1的定量分析结果与图4的定性分析结果对比可知:由于各极化参数的值域范围差异明显,图3对森林AGB敏感性较明显的参数在定量分析时相关性并不一定最高,因此图4的分析仅可作为初步参考,对于森林AGB敏感的极化特征参数的分析,仍需要采用皮尔逊系数进行相关性定量分析。

      表 1  P波段SAR极化特征与森林AGB相关性分析

      Table 1.  Correlation coefficients of polarization feature and forest AGB

      特征RP特征RP特征RP
      HH 0.401** 0 F_VOL −0.129 0.173 alpha −0.734** 0
      HV −0.249** 0.008 Y_DBL 0.719** 0 beta 0.037 0.696
      VV −0.322** 0 Y_ODD 0.323** 0 delta 0.221* 0.018
      FD1 −0.453** 0 Y_VOL −0.213* 0.024 gamma 0.314** 0.001
      FD2 −0.448** 0 Y_HLX 0.187* 0.047 PDR 0.629** 0
      F_DBL 0.650** 0 entropy −0.136 0.151 RVI −0.374** 0
      F_ODD 0.303** 0.001 anisotropy 0.436** 0
        说明:R代表Pearson相关系数;P代表显著性水平,**代表0.01水平上显著;*代表0.05水平上显著

      图  4  4种模型森林AGB估测与LiDAR抽样点散点图

      Figure 4.  Scatter plots of estimated and LiDAR forest AGB of four models

    • 表2可知:4种模型中,MLR模型估测结果精度最低,RF精度最高;而KNN和SVR模型的估测精度则相近,R2相同,RMSE、MAE和Acc也基本相同,比RF估测结果的最高值仅低约2%。

      表 2  基于4种模型的P波段SAR 森林AGB估测结果

      Table 2.  P band SAR forest AGB inversion using four models

      波段反演模型模型参数R2RMSEAcc/%MAE
      P MLR 显著性为0.0008 0.43 19.16 63.55 16.99
      KNN K值为11,欧氏距离 0.54 17.04 71.18 13.44
      SVR 惩罚系数7.0 0.54 17.09 71.15 13.45
      RF 决策树数为100,树深度为11 0.60 15.98 72.97 12.60

      图4可看出:在整个森林AGB分布中,MLR估测结果分布较其他3种非参数方法分散,在森林AGB大于80 t·hm−2时出现了明显的饱和低估现象。其他3种方法尽管也有低估的现象,但是饱和现象并不明显。此外,本研究中4种方法估测结果的相对误差约30%,而以往区域性森林AGB的估测误差为37%~67%[32],采用P波段HV后向散射系数的估测结果中,同质性森林地区的相对误差约13%,而异质性地区相对误差则约60%[28]。在采用L-波段极化分解参数进行森林AGB反演研究中,在森林AGB水平低于120 t·hm−2时,后向散射的特征优于3种极化分解的特征,这与本研究中P波段的研究结果基本一致[33]。4种方法的估测结果中均出现了低值高估和高值低估的现象,这与以往采用不同遥感数据源进行森林AGB估测的结果一致[32]图5中,RF的残差分布最接近高斯分布,尽管峰值出现在10 t·hm−2左右,由于其值分布较窄,因此具有较高的估测精度。KNN和SVR的残差分布图均较为连续,并且分布区间明显高于MLR,但由于正值占比较大,因此总体上呈现高估现象。MLR估测结果的残差分布图在10和20 t·hm−2出现了2个明显的峰值,且所有残差值的分布范围较宽,解释了其估测精度低于其他3种方法的原因。在ENGHART等[34]基于C和L-波段参数和非参数方法森林AGB估测的研究中,也发现MLR的估测结果要略低于非参数模型。然而参数和非参数模型的适用性还受到训练样本大小的影响,因此在训练样本较大时可选择非参数模型进行估测,当样本较小时则可优先选择MLR方法进行森林AGB估测[32]

      图  5  森林AGB估测结果与LiDAR抽样点差值直

      Figure 5.  Histograms of difference between forest AGB estimation results and LiDAR sampling points

      为了进一步分析P波段极化特征对森林AGB估测的潜力,本研究将研究区森林AGB划分为不同的等级,然后采用4种估测模型中估测结果最优的RF模型进行估测,估测结果见表3。在本研究中将森林AGB划分为3组,其中第1组分别以30、60 t·hm−2为边界划分为3个子组;而第2组和第3组分别以40和50 t·hm−2为界划分为2个子组。由表3可知:分组界限不同,估测精度有明显的差异,分组划分越详细,估测精度(Acc)越高;此外,3种分组情况的估测结果均表明:在森林AGB平均值约45 t·hm−2,最高值不超过120 t·hm−2时,P波段在森林AGB水平较高的分组估测精度较高,如在第1组中,森林AGB大于30 t·hm−2分组的估测精度比小于30 t·hm−2分组的估测精约高5%;而在以50 t·hm−2为分组界限的2组中,森林高AGB组的估测精度比低AGB组的估测精度约高出6%。表3的结果分析表明:待估森林的AGB水平对P波段极化特征进行森林AGB估测的估测精度有明显影响,且P波段极化信息更适合森林AGB较高区域森林AGB的估测。

      表 3  基于4种模型的P波段SAR 森林AGB反演情况     

      Table 3.  P band SAR forest AGB inversion based on four models

      序号AGB分段/
      (t·hm−2)
      R2RMSE/
      (t·hm−2)
      Acc/%
      10~300.554.9276.12
      30~600.137.9581.85
      >600.2514.8981.57
      20~500.478.4672.15
      >500.3015.8078.19
      30~400.486.5274.36
      >400.3616.2174.82
    • 针对P波段极化SAR数据在森林AGB估测中的潜力,提取了20个P波段极化SAR参数,探索了森林AGB估测中常用的MLR、KNN、SVR和RF方法在使用P波段极化SAR数据进行森林AGB估测的潜力。结果表明:P波段极化SAR信息在森林AGB估测中具潜力,但估测精度受到待估区域森林AGB水平高低的影响;4种估测方法中,非参数方法的估测结果明显优于MLR估测。P波段森林AGB估测结果中,同样存在AGB低值高估和高值低估的现象,其原因仍需要进一步探索。此外由于本研究区的森林AGB均值为45 t·hm−2,最高值低于120 t·hm−2,所以P波段极化信息在AGB高于120 t·hm−2的森林覆盖区中对森林AGB的估测能力仍有待进一步研究。

参考文献 (34)

目录

/

返回文章
返回