雅鲁藏布大峡谷位于藏东南地区, 北起林芝市米林县, 经排龙乡的雅鲁藏布江大拐弯, 南到林芝市墨脱县。 选取以雅鲁藏布江大拐弯(94° 52'E, 29° 30'N)为中心的3 200 km²范围为研究区(图1), 该区域海拔673~7 198 m, 原始森林覆盖率达到56%以上, 其中, 高山松、 林芝云杉和川滇高山栎是分布最广泛的树种, 分别占森林总面积的54.11%, 35.14%和10.74%。 由于大峡谷地区几乎未受到人为活动的干扰, 且乔木树龄大, 林分郁闭度高, 枯落物丰富, 再加以充足的水汽条件, 使得该地森林具有极高的水源涵养价值。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 研究区地表概况及采样点位置分布Fig.1 Surface profile and sampling points distribution in the study area

高山松和林芝云杉属于常绿针叶树种, 主要分布在雪线灌丛以下较高海拔区域, 川滇高山栎为常绿硬叶阔叶树种, 多分布于高山松林的下坡位。 为保证枯落物的蓄积量保持在较高水平, 于2020年8月下旬正午12:00—14:00, 在雅鲁藏布大峡谷对高山松、 林芝云杉、 川滇高山栎纯林进行野外采样。 样本采集过程如下, (1)建模组: 按照地势条件和土壤类型均相当的原则在高山松林和林芝云杉林选择60棵成年乔木, 在川滇高山栎林选择30棵成年乔木, 剪取每棵乔木不同高度处的健康叶片5张(针叶为5丛)。 在已被采集叶片的乔木下设置小样方采集枯落物样本, 将未分解层和半分解层分开置于取样箱, 准确标记样本编号以确保叶片与枯落物一一对应。 (2)验证组: 在高山松林和林芝云杉林设置15个20 m× 20 m的大样方, 川滇高山栎林设置10个20 m× 20 m的大样方, 每个样方内随机剪取10张(丛)健康成熟的叶片, 在样方四角及中间位置设置20 cm× 20 cm的小样方采集5个枯落物样本, 最后以叶片和枯落物的平均值代表该大样方的实际水平。

1.3.1 叶片高光谱反射率测定

使用ASD Fieldspec 3多功能光谱仪测定叶片光谱, 波长范围为350~2 500 nm, 光谱分辨率3 nm, 光谱采样间隔为1 nm。 采样时, 模拟正午太阳入射光设置线性光源, 视场角为30° , 探头垂直于叶片, 距离约5 cm; 为尽量降低环境和仪器系统误差的干扰, 每15 min对光谱仪进行一次白板校正。

1.3.2 曲线平滑与降噪

光谱信息在采集过程中受大气水分和环境噪声的影响, 使光谱曲线出现波动, 采用移动平均法对光谱曲线进行降噪和平滑处理, 见式(1)

R'i=12k+1(Ri-k+Ri-k+1+…+Ri+…+Ri+k-1+Ri+k)(1)

式(1)中, R'_i为经过平滑处理后的光谱信息值; R_i为原始采集的光谱信息值; i=1, 2, 3, …, n; k=5。 经过上式的运算后, 光谱曲线的每一个点都由该点前后共10个点的平均值代替, 从而达到降低偶然误差和平滑曲线的目的。

哨兵2号(Sentinel-2)卫星携带的MSI多光谱成像仪包含可见光至短波红外共13个光谱波段, 波段空间分辨率分别为10, 20和60 m, 过境周期为5 d。 相对于其他成像卫星产品而言, 哨兵2号数据在红边范围有三个窄波段, 对于植被信息的提取非常有利。 从欧空局哥白尼数据中心(https://scihub.copernicus.eu/dhus/#/home)下载大峡谷地区2020年8月下旬的遥感影像, 经拼接裁剪、 辐射定标、 大气校正、 几何校正和重采样得到可用于计算研究区植被指数的图像。 此外, 为消除不同树种枯落物持水过程的差异性对反演精度的影响, 分别建立高山松、 林芝云杉和川滇高山栎的枯落物有效拦蓄量反演模型, 最后将三类树种的反演结果进行拼接得到大峡谷地区森林枯落物水源涵养能力的空间分布。

考虑到叶片的组织结构、 自身含水量和易分解程度是影响枯落物持水性能的主要因素^{[8, 9, 10, 11]}, 筛选出能反映上述特征的植被指数进行Person相关系数分析: (1)归一化植被指数(NDVI)、 增强植被指数(EVI)、 红边指数(VOG): 表征植被长势、 健康状态、 环境胁迫程度等; (2)叶面叶绿素指数(LCI): 反映叶片叶绿素含量; (3)全球植被水分指数(GVMI): 监测叶片含水量; (4)氮素比值植被指数(NRI): 反映叶片氮素密度; (5)植物蜡质参数(PWP): 与叶片表皮蜡质含量相关性高; (6)植被衰减指数(PSRI): 最大限度地提高类胡萝卜素与叶绿素比率的灵敏度, 预示植被衰老的程度。 上述植被指数的计算方法如表1, 建模时采用ASD光谱仪采集的高光谱数据, 并取各波段中心波长处的反射率计算植被指数。

表1

Table 1

表1(Table 1)

表1 植被指数计算方法 Table 1 Characteristic spectral parameters and calculation method 植被指数 计算方法 NDVI ρNIR-ρRρNIR+ρR EVI 2.5× ρNIR-ρRρNIR+(6×ρR-7.5×ρB)+1 VOG ρRE2ρRE1 LCI ρNIR-ρRE1ρNIR-ρR GVMI (ρNIR+0.1)-(ρMIR+0.02)(ρNIR+0.1)+(ρMIR+0.02) NRI ρGρR PWP ρGρG-ρR PSRI ρR-ρBρRE2 注: ρ NIR, ρ R, ρ B, ρ RE1, ρ RE2, ρ MIR和ρ G分别代表近红外波段、 红光波段、 蓝光波段、 红边第一个窄波段、 红边第二个窄波段、 中红外波段和绿光波段

表1 植被指数计算方法 Table 1 Characteristic spectral parameters and calculation method

取回枯落物样本立即称其自然鲜重, 然后放入85 ℃恒温箱中烘干并称其干重, 计算蓄积量和自然含水率; 将烘干后的枯落物装入垫有尼龙布的土壤筛后浸水, 完全浸泡24 h后取出悬挂, 待其不再滴水时迅速称重, 计算最大持水率(量)和有效拦蓄率(量), 计算公式如式(2)—式(6)。

R0=(M0-M)/M×100%(2)

Rm=(M24-M)/M×100%(3)

Hm=Rm×M(4)

Ru=0.85×Rm-R0(5)

Su=Ru×M(6)

式中, M₀为枯落物单位面积的自然鲜重(t· ha^-1), M为烘干后单位面积的干重(t· ha^-1), 即蓄积量; R₀为自然含水率, R_m为最大持水率, H_m为最大持水量(t· ha^-1); R_u为有效拦蓄率, S_u为有效拦蓄量(t· ha^-1)。

2.1.1 光谱反射率特征

在View Spectral Pro软件中显示ASD Fieldspec 3多功能光谱仪测得的样本数据, 对曲线进行降噪和平滑处理, 而后分别绘制三类树种反射率的平均值曲线, 如图2所示。 大峡谷地区高山松、 林芝云杉和川滇高山栎的光谱反射率总体呈现一致的变化规律, 但各树种的反射波峰(谷)、 高反射平台等特征波段差异显著。 400~700 nm波段(可见光)是植物叶片的强吸收波段, 三组树种的反射率都很低; 但川滇高山栎在< 400 nm的紫外波段表现出较高的反射率, 分析认为川滇高山栎在高海拔环境下生长出较厚的角质膜以抵御强辐射和水分的损失, 角质膜的存在使得叶片对紫外线极为敏感; 林芝云杉和高山松在560 nm(绿光)处的反射峰值非常接近, 但在680 nm(红光)处林芝云杉的吸收性更强。 700~780 nm是植被的反射率红边, 叶片光谱从强吸收红光波段过渡到高反射近红外波段。 780~1 350 nm(近红外)的光谱特征由叶片细胞结构决定, 大量透明的色素和纤维素使得反射率变高, 形成“ 近红外反射肩” ; 川滇高山栎的反射肩最高, 分析认为阔叶树种的叶片结构比针叶复杂, 因此光线在叶片内部散射的次数越多反射率就越高。 此外, 在此波段出现的两个微弱的吸收谷(970, 1 200 nm)是叶片水分所致。 1 350~2 500 nm(中红外)是植物叶片的水分强吸收带, 三组样本均在水分吸收谷(1 450, 1 940 nm)表现出明显的吸收性。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 叶片光谱反射率Fig.2 Spectral reflectance of leaves

2.1.2 光谱一阶导数特征

通过求导可见光范围的光谱曲线可以确定曲线的蓝边、 黄边和红边等特征点, 如图3(a)所示。 通常, 蓝边是光谱反射率在490~530 nm区间内一阶导数峰值对应的波长, 黄边对应的是550~582 nm波段内反射率一阶导数最小值的位置, 红边是680~750 nm范围内反射率一阶导数最大值所在的波长。 植被的蓝边、 黄边和红边位置可以反映植被的生长态势和营养情况; 对红边而言, 植被长势越旺盛、 健康状况越好, 红边位置越向长波方向移动; 反之, 则越短波方向移动。 从表2可以看出, 高山松、 林芝云杉和川滇高山栎的三边位置十分接近, 说明三类树种在光谱反射率的特征点具有共性, 且长势和健康状况处于同等水平。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 光谱一阶导数光谱 (a): 可见光波段; (b): 红边波段Fig.3 First derivative of spectrum (a): The visible light band; (b): The red-edged band

表2

Table 2

表2(Table 2)

表2 蓝边、 黄边、 红边参数 Table 2 Parameter of blue edge, yellow edge and red edge 树种 蓝边位 置/nm 黄边位 置/nm 红边位 置/nm 红边 斜率 红边峰 值面积 高山松 524 567 721 0.009 5 0.404 1 林芝云杉 522 572 717 0.006 3 0.307 6 川滇高山栎 524 568 721 0.011 9 0.523 0

表2 蓝边、 黄边、 红边参数 Table 2 Parameter of blue edge, yellow edge and red edge

如图3(b)所示, 由于植被在红边区域的光谱反射率变化最剧烈, 且该波段蕴含的光谱信息最多, 因此红边的应用还包含了以下几个参数: (1)红边斜率: 红边区域内曲线斜率的最大值, 与叶片叶绿素含量和叶面积指数呈显著正相关; (2)红边峰值面积: 一阶导数曲线在红边范围内与坐标轴围成的面积, 主要受限于叶片叶绿素含量。 结合表2的统计数值, 尽管高山松、 林芝云杉和川滇高山栎的红边位置区别不大, 但红边斜率和红边峰值面积差异明显; 川滇高山栎的红边斜率为0.011 9, 在三类树种中值最大, 高山松次之为0.009 5, 林芝云杉的红边斜率最小, 为0.006 3。 红边峰值面积与红边斜率表现出相同的趋势, 这一趋势也与反射肩平台的反射率规律一致。

2.2.1 枯落物蓄积量和自然含水率

将采集的样本按树种分类, 计算每类树种枯落物的未分解层和半分解层占比、 蓄积总量以及自然含水率的平均值。 从表3看出, 三类树种的半分解层占比由高到低次序为: 林芝云杉> 川滇高山栎> 高山松, 半分解层占比越高说明枯落物分解程度越高, 反之, 则越低。 川滇高山栎作为阔叶树种却没有林芝云杉分解效率高, 可能是因为高山栎叶片表面的蜡质膜阻碍了微生物的分解作用; 高山松因自身含有较多油脂所以最不易分解。 林芝云杉枯落物的蓄积总量(26.46 t· ha^-1)和自然含水率(69.06%)在三类树种中最为突出, 这是因为林芝云杉叶片易于分解且自身含水量较高, 所以凋落物在林下快速腐化并堆积; 另外, 经过分解的枯落物对水分的吸收和贮存能力大于未分解层, 因此林芝云杉的自然含水率也最高。 高山松的蓄积量最低, 为17.53 t· ha^-1, 分析认为未来得及分解的凋落物更易受到外力的冲刷而转移到别处, 从而存有量偏少; 此外, 高山松林枯落物的自然含水率(30.84%)明显低于林芝云杉(69.06%)和川滇高山栎(62.89%), 一方面是因为高山松作为针叶树种, 其林冠层截留量大于阔叶树种, 致使林内枯落物相对干燥, 另一方面是高山松叶片的自身含水率低于同为针叶树种的林芝云杉, 且相对低的分解程度不利于其保留水分。

表3

Table 3

表3(Table 3)

表3 枯落物蓄积量和自然含水率 Table 3 Litter Stock volume and natural water content 树种 蓄积量 自然 含水 率/% 未分解层 占比/% 半分解层 占比/% 总量 /(t· ha-1) 高山松 53.68 46.32 17.53 30.84 林芝云杉 29.18 70.82 26.46 69.06 川滇高山栎 41.17 58.83 25.14 62.89

表3 枯落物蓄积量和自然含水率 Table 3 Litter Stock volume and natural water content

2.2.2 枯落物持水量和拦蓄量

根据枯落物浸水后的质量计算三类枯落物的最大持水量(率)和有效拦蓄量(率), 如表4。 可见, 林芝云杉的最大持水率和最大持水量均为最高, 高山松和川滇高山栎最大持水率相当, 但由于川滇高山栎的蓄积量更丰富, 所以最大持水量高于高山松。 但是, 最大持水率和最大持水量只能说明枯落物持水能力的强弱, 有效拦蓄率和有效拦蓄量才能表征枯落物对降水的实际拦蓄效果。 对本研究的三类树种而言, 有效拦蓄率和有效拦蓄量的排序分别为: 高山松> 林芝云杉> 川滇高山栎, 林芝云杉> 川滇高山栎> 高山松。 高山松的枯落物在自然状态下最干燥, 可以用来吸收外界水量的空间最大, 造成其有效拦蓄率排名最靠前, 但由于较少的蓄积量, 导致有效拦蓄量排名滞后。 分解程度最高的林芝云杉枯落物拥有和高山松相当的有效拦蓄率, 自身蓄积量也大, 因此对降水的有效拦蓄量最大。 森林枯落物的水源涵养能力是对枯落物层整体功能的描述, 包含对枯落物蓄积量和吸水性能的综合定量评价, 因此, 在描述持水拦蓄性能的各项参数中, 有效拦蓄量最能直观地反映出枯落物的水源涵养能力。

表4

Table 4

表4(Table 4)

表4 枯落物持水拦蓄性能 Table 4 Water holding capacity of litter 树种 最大持水 率/% 最大持水量 /(t· ha-1) 有效拦蓄 率/% 有效拦蓄量 /(t· ha-1) 高山松 253.19 44.38 184.37 32.32 林芝云杉 296.25 78.39 182.75 48.36 川滇高山栎 257.61 64.76 156.08 39.24

表4 枯落物持水拦蓄性能 Table 4 Water holding capacity of litter

2.3.1 植被指数筛选

分别计算每类树种拟选定的8种植被指数(NDVI, EVI, VOG, LCI, GVMI, NRI, PWP和PSRI)与枯落物有效拦蓄量的Person相关系数, 结果如表5。

表5

Table 5

表5(Table 5)

表5 Person相关系数分析 Table 5 Person correlation coefficient analysis 植被指数 Person相关系数 高山松 林芝云杉 川滇高山栎 NDVI 0.481* 0.277 0.204 EVI 0.288 0.331* 0.231 VOG 0.345 0.312* 0.358* GVMI 0.507* 0.546* * 0.407* LCI 0.396* 0.235 0.375* NRI 0.194 0.363* 0.289 PWP -0.623* * -0.437* -0.565* * PSRI -0.519* * -0.544* -0.442* * 注: * 显著相关(p< 0.05), * * 极显著相关(p< 0.01)

表5 Person相关系数分析 Table 5 Person correlation coefficient analysis

筛选出与有效拦蓄量显著相关和极显著相关的植被指数作为回归模型的自变量, 按照系数由大到小的顺序排列, 各树种参与回归的自变量如下: (1)高山松: PWP, PSRI, GVMI, NDVI, LCI; (2)林芝云杉: GVMI, PWP, PSRI, NRI, EVI, VOG; (3)川滇高山栎: PWP, PSRI, GVMI, LCI, VOG。 从保留的植被指数中可以看出, 三类树种林下枯落物的有效拦蓄量都与叶片蜡质参数(PWP)、 植被长势(NDVI、 EVI、 VOG)、 衰减程度(PSRI)、 叶片含水量(GVMI)相关; 此外, 林芝云杉的拦蓄能力还与叶片氮素含量(NRI)相关, 高山松和川滇高山栎的拦蓄能力还与叶绿素含量(LCI)相关。

2.3.2 水源涵养反演模型

将有效拦蓄量作为因变量, 筛选后的植被指数作为自变量代入SPSS22.0软件中逐步回归, 得到枯落物水源涵养能力的反演模型(表6)。 高山松、 林芝云杉和川滇高山栎的模型拟合优度(R²)分别为0.691, 0.779和0.743, 表明模型拟合效果好, 可以较好地反映枯落物水源涵养能力与植被指数的关系。 此外, 植被叶片蜡质参数和衰减指数与有效拦蓄量呈负向相关, 说明叶片蜡质含量越高则越不易分解, 因此蓄水能力越弱, 同理, 叶片衰老程度越高其细胞机能越低, 吸水能力越弱。 其余植被指数与有效拦蓄量呈正向相关, 说明植被长势越好其枯落物的水源涵养能力越强。

表6

Table 6

表6(Table 6)

表6 各树种枯落物层水源涵养能力反演模型 Table 6 Inversion model of water conservation in litter layer of different tree species 植被类型 反演模型 R2 高山松 WL=-13.78PWP* * -10.69PSRI* +8.24GVMI* +5.51NDVI* +2.82LCI* +69.03 0.691 林芝云杉 WL=21.43GVMI* * -15.77PWP* * -10.81PSRI* +9.38NRI* +7.88EVI* +2.51VOG* +20.14 0.779 川滇高山栎 WL=-15.16PWP* * -12.73PSRI* +12.62GVMI* +9.64LCI* +6.79VOG* +72.11 0.743 注: WL为枯落物层的水源涵养量(t· ha-1), * 显著相关(p< 0.05), * * 极显著相关(p< 0.01)

表6 各树种枯落物层水源涵养能力反演模型 Table 6 Inversion model of water conservation in litter layer of different tree species

以Sentinel-2影像为底图, 基于上述模型分别反演高山松、 林芝云杉和川滇高山栎的枯落物水源涵养能力, 经空间拼接处理后得到图4。 水源涵养的高值主要出现在林芝云杉林区内部, 少数出现在川滇高山栎林; 从各林种内部来看, 海拔越低, 水源涵养能力越强, 在近河谷的两侧, 水源涵养值达到最高。 总体而言, 大峡谷研究区内森林枯落物层的水源涵养情况沿植被类型和地形大致呈带状分布, 不同植被枯落物的水源涵养能力差异明显。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 大峡谷森林枯落物水源涵养能力分布情况Fig.4 Water conservation and distribution of forest litter in Grand Canyon

将验证组的叶片光谱反射信息代入各自树种的反演模型进行计算, 得到枯落物的预测水源涵养量, 并与实际测得的有效拦蓄量进行比对, 模型的精度评价结果分别如图5(a, b, c)所示。 高山松、 林芝云杉和川滇高山栎的检验点R²分别为0.943, 0.815和0.812, RMSE分别为1.597, 2.270和1.953, 表明模型的预测效果良好, 可以用于大峡谷地区森林枯落物水源涵养能力的预测分布研究。

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 枯落物水源涵养能力反演模型精度评价 (a): 高山松; (b): 林芝云杉; (c): 川滇高山栎Fig.5 Precision evaluation of water conservation capacity for different species litter layer in model (a): Pinus densata; (b): Picea Linzhi; (c): Quercusaqui folioides