0 引言
1 研究区与研究方法
1.1 研究区概况
1.2 估算流程
2 数据获取和处理
2.1 遥感数据获取和处理
2.1.1 无人机多光谱数据获取
2.1.2 波段选择
2.1.3 乔木、灌丛优势种的单木识别和分割
2.2 样方设置
2.3 地面数据获取和处理
表1 优势树种冠幅与胸径(基径)参数Tab.1 The parameters for crown and diameter at breast height(branch) diameter of dominant plant species |
物种 | 样本量 | 胸径或基径/cm | 冠幅/m2 | |||
---|---|---|---|---|---|---|
范围 | 标准差 | 范围 | 标准差 | |||
台湾相思 | 15 | 4.80~18.90 | 4.06 | 3.00~48.00 | 14.89 | |
椿叶花椒 | 15 | 6.30~23.00 | 5.75 | 4.00~70.68 | 22.73 | |
黑松 | 15 | 6.80~11.70 | 5.15 | 4.94~16.38 | 4.18 | |
天仙果 | 15 | 6.80~10.70 | 0.93 | 8.78~19.32 | 2.74 | |
野梧桐 | 15 | 1.30~8.77 | 1.85 | 47.75~11.40 | 1.94 | |
滨柃 | 15 | 6.90~9.10 | 0.65 | 4.32~7.83 | 1.23 |
表2 灌丛优势种基径与生物量Tab.2 The branch diameter and biomass of dominant shrubs |
类型 | 物种 | 样本量 | 基径/cm | 生物量/kg | |||
---|---|---|---|---|---|---|---|
范围 | 标准差 | 范围 | 标准差 | ||||
灌丛 | 野梧桐 | 7 | 2.20~6.80 | 1.38 | 0.08~1.33 | 0.44 | |
天仙果 | 8 | 2.23~5.06 | 0.89 | 0.28~1.51 | 0.65 | ||
滨柃 | 12 | 1.08~5.38 | 1.51 | 0.05~2.36 | 0.81 |
表3 草丛样方的生物量与归一化植被指数Tab.3 The biomass and normalized difference vegetation index(NDVI) of herb plots |
样地编号 | 生物量/kg | NDVI | 样地编号 | 生物量/kg | NDVI | 样地编号 | 生物量/kg | NDVI |
---|---|---|---|---|---|---|---|---|
H13 | 1.04 | 0.88 | H23 | 0.32 | 0.34 | H33 | 0.79 | 0.73 |
H14 | 0.15 | 0.23 | H24 | 0.83 | 0.71 | H34 | 0.69 | 0.66 |
H15 | 0.22 | 0.28 | H25 | 0.94 | 0.79 | H35 | 0.24 | 0.34 |
H16 | 0.20 | 0.26 | H26 | 0.53 | 0.49 | H36 | 0.37 | 0.34 |
H17 | 0.17 | 0.24 | H27 | 0.31 | 0.37 | H37 | 0.13 | 0.26 |
H18 | 0.43 | 0.44 | H28 | 0.37 | 0.37 | H38 | 0.21 | 0.22 |
H19 | 0.76 | 0.67 | H29 | 0.15 | 0.25 | H39 | 0.30 | 0.40 |
H20 | 0.25 | 0.29 | H30 | 0.57 | 0.51 | H40 | 0.71 | 0.70 |
H21 | 0.40 | 0.42 | H31 | 0.59 | 0.52 | H41 | 0.17 | 0.31 |
H22 | 0.08 | 0.16 | H32 | 0.30 | 0.31 | H42 | 0.68 | 0.53 |
2.4 生物量反演
表4 回归方程方差分析Tab.4 Variance analysis of regression equations |
回归方程 | a | b | R2 | F值 | SEE | p值 | ||
---|---|---|---|---|---|---|---|---|
生物量 反演 | 天仙果 | Y=a+bX | -1.374 1 | 0.672 3 | 0.847 0 | 33.204 0 | 0.276 0 | 0.001 0 |
Y=aXb | 0.034 2 | 2.565 7 | 0.864 5 | 38.287 0 | 0.284 0 | 0.001 0 | ||
野梧桐 | Y=a+bX | -0.795 5 | 0.467 3 | 0.756 0 | 15.521 0 | 0.438 0 | 0.011 0 | |
Y=aXb | 0.030 3 | 2.391 5 | 0.938 0 | 91.717 0 | 0.340 0 | 0.000 2 | ||
滨柃 | Y=a+bX | -0.462 8 | 0.275 3 | 0.729 0 | 30.594 0 | 0.228 0 | 0.000 3 | |
Y=aXb | 0.027 0 | 2.154 4 | 0.738 0 | 31.913 0 | 0.430 0 | 0.000 2 | ||
中华结缕草 | Y=a+bX | -158.480 0 | 1 352.445 0 | 0.967 0 | 809.549 0 | 49.863 0 | 0.000 0 | |
Y=aXb | 1 384.025 0 | 1.483 0 | 0.934 0 | 394.064 0 | 0.177 0 | 0.000 0 | ||
胸径(基径)与冠幅 | Y=a+bX | 5.837 2 | 0.252 0 | 0.845 0 | 484.655 0 | 1.490 0 | 0.000 0 | |
Y=aXb | 3.578 0 | 0.382 0 | 0.526 0 | 97.508 0 | 0.258 0 | 0.000 0 |
表5 大竹峙岛植被优势种生物量方程Tab.5 Biomass equations for dominant species on Dazhuzhi Island |
植被类型 | 物种名称 | 生物量方程 | 方程来源 |
---|---|---|---|
乔木 | 台湾相思/ 椿叶花椒 | W干=0.089 5×D2.452 1; W枝=0.020 5×D2.505 9; W叶=0.021 5×D2.039 3; W根=0.006 7×D2.877 4; W=W干+W枝+W叶+W根 | 文献[14] |
灌丛 | 黑松 | W干=0.067×D2.409; W枝=0.022×D2.27; W叶=0.025×D2.124; W根=0.038×D2.165; W=W干+W枝+W叶+W根 | 文献[14] |
灌丛 | 天仙果 | W=0.034 2×D2.565 7 | 实测构建 |
滨柃 | W=0.027 0×D2.154 4 | 实测构建 | |
野梧桐 | W=0.030 3×D2.391 5 | 实测构建 | |
草丛 | 中华结缕草 | W=1 352.445 0×NDVI-158.480 0 | 实测构建 |
注:W表示单个植株生物量,D表示植株胸径(基径),D=0.252 0×P+5.837 2,P表示植株冠幅。 |
2.5 植被碳储量计算
3 结果与分析
3.1 大竹峙岛植被碳储量
表6 大竹峙岛乔木、灌丛和草丛碳储量估算结果Tab.6 Carbon storage estimation results of arbors, shrubs and herbs on Dazhuzhi Island |
物种 名称 | 株数 | 面积/hm2 | 总生物量/t | 总碳储量/t | 平均碳密度/ (t·hm-2) |
---|---|---|---|---|---|
乔木 | 28 856 | 11.69 | 667.47 | 300.36 | 25.69 |
灌丛 | 22 597 | 7.98 | 105.75 | 47.59 | 5.96 |
草丛 | 9.07 | 47.98 | 21.59 | 2.38 |
3.2 估算结果分析
3.2.1 无人机遥感估算结果精度验证
图8 样方内无人机遥感估算和样方调查结果的相关性(图8a表示9个乔木和灌丛样方内无人机遥感估算株数和实地调查结果的相关性;图8b~图8d分别表示3个乔木样方、6个灌丛样方和12个草丛样方内无人机遥感估算生物量和实地测定结果的相关性。) Fig.8 Correlation between UAV multispectral remote sensing and quadrat measurement in sample plots (Fig.8a shows the correlation of the plant numbers between UAV multispectral remote sensing and quadrat measurement in shrub and arbor sample plots; Fig.8b to Fig.8d show the correlation of biomass between UAV multispectral remote sensing and quadrat measurement in 3 arbor sample plots, 6 shrub sample plots and 12 herb sample plots, respectively.) |
表7 乔木、灌丛植株数、生物量的遥感反演和样方调查比较Tab.7 Comparison of arbor and shrub plant number and biomass by remote sensing inversion and quadrat measurement |
植被 类型 | 样方 编号 | 株数 | 生物量 | |||||
---|---|---|---|---|---|---|---|---|
实测 /株 | 单木分 割/株 | 样方调 查/kg | 遥感估 算/kg | 差值/kg | 相对偏 差/% | |||
乔木 | T1 | 86 | 80 | 6 397.01 | 5 231.35 | -1 165.66 | 18.22 | |
T2 | 91 | 95 | 2 171.10 | 2 420.34 | 249.24 | 11.48 | ||
T3 | 46 | 39 | 1 006.23 | 879.42 | -126.81 | 12.60 | ||
灌丛 | S1 | 09 | 11 | 35.28 | 27.96 | -7.32 | 20.75 | |
S2 | 7 | 9 | 219.99 | 192.22 | -27.77 | 12.62 | ||
S3 | 7 | 5 | 34.95 | 23.48 | -11.47 | 32.82 | ||
S4 | 15 | 5 | 36.52 | 25.09 | -11.43 | 31.30 | ||
S5 | 6 | 5 | 14.28 | 10.03 | -4.25 | 29.76 | ||
S6 | 6 | 5 | 86.41 | 63.18 | -23.23 | 26.88 |
表8 草丛生物量的遥感反演和样方调查比较Tab.8 Comparison between remote sensing inversion and quadrat measurement of herb biomass |
样地 编号 | NDVI | 生物量 | |||
---|---|---|---|---|---|
样方调查/kg | 遥感估算/kg | 差值/kg | 相对偏差/% | ||
H1 | 0.27 | 0.21 | 0.21 | 0.00 | 0.00 |
H2 | 0.25 | 0.18 | 0.17 | -0.01 | 5.56 |
H3 | 0.25 | 0.18 | 0.18 | 0.00 | 0.00 |
H4 | 0.22 | 0.12 | 0.14 | 0.02 | 16.67 |
H5 | 0.40 | 0.37 | 0.39 | 0.02 | 5.41 |
H6 | 0.31 | 0.24 | 0.26 | 0.02 | 8.33 |
H7 | 0.45 | 0.44 | 0.45 | 0.01 | 2.27 |
H8 | 0.38 | 0.33 | 0.35 | 0.02 | 6.06 |
H9 | 0.69 | 0.74 | 0.77 | 0.03 | 4.05 |
H10 | 0.24 | 0.22 | 0.17 | -0.05 | 22.73 |
H11 | 0.65 | 0.79 | 0.72 | -0.07 | 8.86 |
H12 | 0.68 | 0.84 | 0.76 | -0.08 | 9.52 |