受气候变化影响, 祁连山东部温度和降雨均逐年上升[1], 为适应其引起的立地条件变化, 在立地条件改善区进行植树造林十分必要[2]。目前已有研究揭示在耕地和裸地进行造林的效果差异[3—6], 但鲜有在灌木林地上进行植树造林的研究, 本次研究区位于祁连山甘沟小流域, 原有土地利用类型为灌草地, 成为潜在的灌草转为森林的区域, 为探究该地区植被转换生态环境效益, 本研究对造林地和原有灌草地两种植被类型进行比较。

传统造林效果评价的研究方法集中在相同区域造林地与未经造地土壤理化性质、物种多样性、生物量等指标的比较[3, 5—10], 此方法需要耗费大量人力和物力且无法形成长时间大尺度观测, 而遥感技术的发展弥补了传统调查的不足, 以其高时空分辨率、全球尺度监测的优势被广泛应用于植被调查[1, 11—12], 植被水分监测[13, 14]。特别是具有高时空分辨率以及其易于监测植被健康的红边波段[15, 16]的哨兵2号(Sentinel-2)数据的出现, 促进了遥感在高精度植被健康状况监测的应用[11, 13, 16—18]。目前以有大量用于监测植被状况的遥感指数如归一化植被指数(NDVI), 增强植被指数(EVI), 水分胁迫指数(MSI), 植被水分指数(NDMI), 另外由哨兵数据红外波段衍生的准确度较高的陆地叶绿素指数(MTCI)[17, 19], 和叶绿素红外指数(CI)[18]也被验证和应用。但传统手段对遥感数据下载、处理等环节较为繁琐, 限制其广泛应用, 近年来随着Google earth engine(GEE)遥感云处理平台的发展, 其以海量数据库、云端处理数据无需下载及可重复的优势[20], 加速了遥感在植被监测中的应用。

本文基于GEE平台的高精度sentinel-2数据计算获取的植被指数并结合传统实地调查方法对门源地区灌草地造林与未造林地两种植被类型在短时间内生态环境状况进行比较, 揭示短时间造林地与原有灌草地之间的差异, 为森林植树管理及监测提供借鉴, 并探讨遥感在森林监测中的应用不足和前景。

本次研究位于青海省东部门源县甘沟小流域(37.5°—38°N, 101°—102°E), 海拔2 800—3 100 m, 属高原大陆性气候, 光照时间长, 日照强烈, 多年平均气温2.8—7.9℃、降水量360—540 mm (集中在6—9月)、蒸发量1100—1800mm。土壤垂直地带分布以栗钙土为主。甘沟地区原先有稀疏乔木, 但被当地居民破坏形成灌草地用于放牧, 为改善当地生态环境, 于2019年开始对灌木地区进行幼苗栽植造林, 造林树种有青海云杉(Picea crassifolia), 白桦(Betula platyphylla Suk), 小叶杨(Populus simonii)；原有灌木种类为金露梅(Potentilla fruticosa), 银露梅(Potentilla glabra), 沙棘(Hippophae rhamnoides), 小檗(Berberis thunbergii), 金银忍冬(Lonicera maackii), 鲜卑花(Sibiraea laevigata)。

本次调查于2021年7月份进行, 共调查13个样方, 6个灌草地和7个乔木林地(见表 1), 样地大小为24m×24m。所有样地均设置在阴坡, 并用GPS定位每一个样地的坐标、海拔、坡向且以5m等高线进行校准, 地理坐标系为WGS1984(EPSG：4326)。每个样地随机取两个5m×5m的灌木样方对灌木的体积进行测算, 五个1m×1m样方进行草本植物多样性调查并测定土壤水分, 一个区域对0—20, 20—40, 40—60cm三层取土分别测量其有机质含量。

土壤水分用TRIME-PICO 64/32 TDR便携式土壤水分测量仪测量, 采样土壤采用GB7858-87标准进行土壤有机质测量。

基于Google Earth Engine(GEE)平台获取的Sentinel-2数据ID: (COPERNICUS/S2_SR), 具有5天的重返周期, 且大部分波段精度为10m(见表 2), 为小尺度范围内监测提供了更多可能。数据经CFMask算法[21]生成的QA波段进行掩膜去云。该地区7—8月为植物生长季降水和温度均最高, 4—5月为植物发芽季, 因此采用NDMI和MSI在每年4—5月最小和最大值, 在7—8月最大和最小值的三年六个值对植物水分状况进行监测, 其余指数均采用最大值合成法(Maximum Value Composite, MVC)消除误差获取2019、2020、2021三年的年最大值。

首先计算草本Margalef丰富度指数、Shannon多样性指数与Pielou均匀度指数用来指示α多样性的分布, 计算公式如下：

Margalef丰富度指数

Shannon多样性指数

Pielou均匀度指数

R为Margalef丰富度指数, H′为Shannon多样性指数, E为Pielou均匀度指数, N为样地内个体数量, S为样地内的总物种数, pi为第i物种所占比例。

在GEE平台上基于地理坐标系为WGS1984(EPSG：4326)获取的实地勘测的坐标进行掩膜取值, 获得15个样地的三个类别6个植被指数, 分别为用以表征植被覆盖情况的结构指数：植被归一化指数(NDVI)[22], 增强植被指数(EVI)[23]；用以表征植被水分状况的水分指数：植被水分指数(NDMI)[24], 水分胁迫指数(MSI)[25]；用以表征植被健康状况的红外波段衍生的理化指数：叶绿素红外指数(CI) 用于观测冠层内的叶绿素含量[18, 26]和用以观测绿素含量的陆地叶绿素指数(MTCI)[19], 计算公式见表 3。

分别比较乔木林地和灌草地在土壤含水量, 土壤有机质含量, 草本植物多样性, NDVI, EVI, NDMI, MSI, CI, MTCI上的差异性(P < 0.05), 利用shapiro (P>0.05)检验其正态性, 通过检验使用单因素方差分析(One-way ANOVA), 未通过使用Wilcoxon分析, 数据以箱线图和柱状图并标有P值和使用方法呈现。以上数据的处理及画图和分析均在R中stats包[27] tidyverse包[28]中进行。

比较灌草地和乔木林近三年的植被指数算术平均值, 仅有水分指数在两种植被类型比较中出现显著性差异, 其中4—5月份MSI灌木林大于(P < 0.01)乔木林, 分别为1.69, 1.55；NDMI指数中乔木林大于灌草地(P < 0.01)分别为-0.213, 灌木为-0.253。7—8月份水分指数灌草地和乔木林基本相同, NDMI分别为0.322, 0.325, MSI分别为0.514, 0.510。结构指数和叶绿素指数中灌草与乔木林之间无显著差异(P>0.05), 乔木林和灌木林的NDVI分别为0.836, 0.840, EVI分别为0.693, 0.714, 两植被类型的结构指数均属于植被覆盖高的范围, 但乔木还有较大的生长空间, 而灌草地基本饱和。叶绿素指数中, MTCI乔木林为-0.702, 灌草地为-0.688, CI乔木为-0.879, 灌木为-0.844。其分布见图 1。

比较灌草地和乔木林的实地调查指标的算术平均值, 含水率灌草地显著(P < 0.01)高于林地, 林地和灌草地分别为33.35%, 31.60%。土壤有机质灌草地略微高于乔木林(P>0.05)灌草地和乔木林分别为114.48g/kg, 为93.51g/kg。两种植被类型土壤有机质在0—20cm, 20—40cm, 40—60cm分别为灌草地(137.43 g/kg, 101.15 g/kg, 81.69 g/kg)和乔木林(103.28 g/kg, 101.15 g/kg, 81.69 g/kg), 灌草地均高于乔木林且均随着深度的增加而减小。Margalef丰富度指数, Shannon多样性指数和Pielou均匀度指数草地和灌草地均无显著差异, 灌草地和乔木林的Margalef丰富度指数分别为1.76, 1.80；shannon多样性指数分别为1.67, 1.77；Pielou均匀度指数分别为0.69, 0.73。其分布见图 2。

所有指标相关性分析表明(图 3), 植被覆盖度与MSI成负相关与NDMI成正相关(P < 0.05)。另外CI指数与土壤有机质和Shannon多样性指数均显著正相关为0.61(P < 0.05)。Shannon多样性指数与MSI负相关为-0.60(P < 0.05), 其余植被指数与实地指标不相关(P>0.05)。通过筛查各指数发现, EVI和NDVI与植被覆盖度其相关系数较低是由于覆盖度较低的灌草地具有较高的植被指数值(图 4)。

水分是限制树木生长的重要因素, 目前大量造林失败由于土壤水分的缺失[1, 3], 本文通过遥感指数观测的植被水分发现, 4—5月份, 灌木的含水量显著低于乔木, 而7—8月份两种植被类型水分基本相同, 这可能是由于, 4—5月份降雨不充沛条件下, 灌草比乔木耗水少[29]。7月份的实地调查发现灌木的土壤水分显著高于乔木林, 但其植被水分基本相同, 造成这种现象可能由于两个原因：一是灌丛相对于乔木林倾向于使用深层地下水[30], 而本次只对土壤浅层水分进行了测定, 二是可能由于乔木相对于灌木需要更多水分生长或者是因为乔木林的盖度普遍高于灌木林, 耗水较多, 但王凯博等[4]发现灌草地的土壤水分较乔木林更低, 与本研究结果相悖, 可能由于其研究的森林郁闭度较高蓄存更多的水分, 而本次研究树木尚未郁闭, 太阳辐射强, 导致土壤水分较低。另外也有研究表明过度造林会造成土壤水分的干化并且随着时间推移效果会加剧[3, 31], 本次研究仅发现乔木林浅层土壤水较低, 且并未发现干化现象, 表明现有短时间造林下土壤水分并未成为造林区树木生长的限制条件, 且植被也并未遭受水分胁迫。

土壤是植物生长的基础, 也是植物对其生长环境的反馈, 良好的土壤环境能实现造林的健康生长。已有研究表明植树造林可以提高土壤有机质的含量[5—6], 但其研究多集中在原有农田或草地等土壤相对贫瘠地区。在天然的森林和灌木中两种植被类型土壤有机质含量并不明确且受诸多环境因素影响[8, 32], Hong等[7]发现植树造林对土壤有机质的效应由原有区域土壤有机质含量决定, 一般表现为土壤贫瘠地区增加而土壤肥沃地区减少, 本次研究区原有灌草地土壤有机质算术平均值为114.48g/kg属于有机质丰富地区, 造林后出现有机质减少与其研究结果一致, 但造林对土壤有机质的影响无法抵消其将来对整个生态环境的积极作用。

草本多样性受诸多因素影响, Bremer等[33]通过全球荟萃分析发现在非退化土地进行造林会降低其草本多样性, 但其影响还受树木年龄、气候条件、造林措施等诸多因素影响。本研究树木尚处于幼林阶段未形成郁闭对草本获取资源的影响较小[34], 从而草本多样性未发生显著降低, 这也从另一方面证明该地区目前资源条件可以支持乔木生长而不影响其他植被生长。遥感监测中, 反应植被健康的遥感指数并未发生显著性差异, 结构性指数中两者基本相同, 而叶绿素灌草地略高于乔木林(P>0.05), 可能由于乔木林尚未郁闭木质部分占比过大影响监测效果。因此从以上结果可以得出现阶段短时间造林对原有植被和环境并未产生负效应。

遥感技术在森林状况监测的发展克服了传统森林调查的不足, 能够实现及时大尺度的森林监测, 如森林退化[12]和水分的胁迫[13], 但由于遥感数据时空精度的问题限制其在大尺度范围的应用, 但随着GEE的发展减化了遥感应用的繁琐步骤, 以及sentinel-2高精度开源数据的出现提高了小尺度遥感监测的精确性, 为小尺度监测提供了便利。本次以单样地进行皮尔逊相关分析发现CI与土壤有机质和Shannon多样性指数均成显著正相关, 可能由于土壤有机质和生物多样性能够为植被提供较好的生存环境因而植被具有较强的光合作用能力, 但由于本次研究样本量和区域性限制其关系还需进一步验证。另外EVI和NDVI高估了灌草地植被盖度, 可能由于两方面原因, 一方面乔木年龄过小, 叶面积指数较小, 造成了其主要以树干为主的生物量难以被监测, 另一方面本次研究对草本的生物量没有进行准确估测, 灌草地内大量的草本造成EVI和NDVI偏大。NDVI和EVI同为植被结构指数中两种植被类型比较发现, 中位数表现有差异, 这可能由于NDVI在植被茂密地区易发生过饱和现象, 而EVI能在饱和区域发现差异[23], 灌草地的郁闭度较低, 无法充分利用EVI的特点从而造成此差异, 因此在不同的植被类型和生长阶段应慎重考虑不同指数的应用特性再进行选择。在水分指数中均未发现土壤水分与遥感监测植被水分的关系, 可能由于该地区目前水分条件较好, 土壤水分的变化对植被水分影响较小[25]。综上可得出, 遥感指数对植被的监测受植被类型和土壤环境上及其他因素仍有较多的不确定性[17, 35], 在小尺度上的应用分析仍需进一步研究。

本论文结合遥感指数和实地调查指标对灌草地和灌草地栽植的乔木林进行比较, 为潜在植被类型转换地区植树造林提供了借鉴。研究区除水分相关的遥感指数和调查指标存在显著性差异外, 其余指标基本相同, 通过植被叶绿素和结构指数发现乔木林生长并没有受到影响, 草本多样性未发生降低也证明该地区资源条件能够支持乔木生长, 表明该地区短时间造林通过牺牲一部分土壤水分提高了植被盖度, 并维持了生态系统的稳定, 但其对生态环境改善并不明显。

本次研究验证了遥感技术在小尺度范围内监测植被状况的可行性, 对于遥感监测的不确定性进行了分析, 为后续遥感植被监测提供经验。以后遥感植被监测应充分发挥哨兵系列等高频率、高精度数据以及GEE等云处理平台, 结合不同指标降低干扰以便能够更加准确的对植被进行监测。