草地是陆地生态系统的重要组成部分, 在调节气候、保持水土、改良土壤和维持生物多样性等方面发挥着重要作用。植被覆盖度及其变化深刻影响着生态系统的植被生产能力、土壤理化性质、产汇流过程和生态系统功能, 从而影响着土壤流失等地表物质的传输和平衡, 是土壤侵蚀评价中的重要影响因素^[1]。内蒙古地区位于中国北方, 气候主要为干旱、半干旱和半湿润区, 是我国北方草地的主体部分和重要的农牧业生产基地, 生态环境脆弱, 土地退化严重^[2]。研究内蒙古不同类型草地植被在年、季、月等尺度上对降水的响应和敏感程度, 有助于揭示草地覆盖变化影响因素, 提高区域土壤侵蚀动态评价效率。

很多研究结果认为, 近年来内蒙古植被覆盖度总体呈上升趋势^[3], 但不同地区植被覆盖度变化趋势表现不同。如彭飞等^[4]、李林叶等^[5]发现2000年后的15年中呼伦贝尔草原植被覆盖度整体下降；焦全军等^[6]、张圣微等^[7]认为, 锡林郭勒草原植被覆盖2001—2014年总体呈增加趋势, 不同年份各类型草原植被覆盖度表现为由东向西递减的规律；而希拉穆仁草原近年来则处于中度退化状态^[8]。在植被覆盖度对降水的响应方面, 大量研究结果认为, 降水是影响植被变化最主要的自然因素, 是制约中国北方温带草原植被生长的根本原因。孙艳玲等^[9]、张清雨等^[10]认为, 1982—2000年内蒙古地区植被NDVI与降水有很好的相关性, 植被受降水的影响较大。许旭等^[11]、穆少杰等^[3]认为, 2001—2010年内蒙古草原生态区植被覆盖度在年际水平上主要受降雨影响, 月植被覆盖度对降雨量的响应存在时滞效应。上述研究均从植被生态角度, 探讨了植被覆盖度的变化及降水对植被生长的一般影响, 而对不同类型草原植被覆盖度的多年变化状态, 对降水量变化的敏感程度和定量关系, 植被覆盖度对不同时段降水量的响应过程和密切程度, 最佳模拟式等等问题还有待研究。

植被覆盖度(Fractional Vegetation Cover, FVC)是植被的叶、茎、枝的垂直投影面积占地面总面积的比例, 是衡量地表植被状况的重要指标^[12]。地表植被按照功能属性可以分为绿色光合植被(Photosynthetic Vegetation, PV)和非光合植被(Non-photosynthetic Vegetation, NPV)。在以往像元二分模型基础上, Guerschman等^[13]在2009年提出了像元三分模型, 将光合植被(PV)、非光合植被(NPV)和裸土(BS)作为混合像元的3种组分进行线性分离, 经过澳大利亚稀树草原区域的验证、评估和改进, 推出了基于MODIS数据(MCD43A4, MOD09A1)的全球产品^[13]。在我国北方草原地区的应用^[14-16]表明, 像元三分模型中的PV覆盖度精度可以达到91.2%, 而NPV、BS的反演还存在很大不确定性。

本文将使用该数据集的光合植被(PV)数据部分, 深入分析内蒙古不同类型草原年、月植被覆盖度变化, 及其对不同时间尺度降水量的响应过程和特征, 为区域土壤侵蚀动态评价和生态环境保护提供理论依据。

1 研究区和样地选择

内蒙古自治区位于我国的北部边疆, 由东北向西南斜伸, 总面积118.3万km²。其地貌以蒙古高原为主体, 主要包括由呼伦贝尔、锡林郭勒、巴彦淖尔—阿拉善及鄂尔多斯等高平原, 平均海拔1000 m左右。以温带大陆性季风气候为主, 平均降水量50—450 mm, 平均气温-3.7—11.2 ℃。地带性土壤的分布从东部的黑土, 向西依次为暗棕壤、黑钙土、栗钙土、棕壤土、黑垆土、灰钙土、风沙土和灰棕漠土。按照降雨量和温度的梯度变化, 植被类型沿东北西南一线依次为森林、草原和荒漠^[17]。随气候条件的地带性差异形成了草原植被的3个亚类型, 即草甸草原、典型草原和荒漠草原。

本文结合中国植被图和Google Earth以内蒙古3种草原类型, 即呼伦贝尔市鄂温克族自治旗的草甸草原区、锡林郭勒盟西乌珠穆沁旗的典型草原区、包头市达尔罕茂明安联合旗的荒漠草原区为研究对象, 各草原类型区分别选取3个代表性研究样地, 每个样地均包括9个像元, 即每个草原区的样地面积均为6.75 km², 取像元均值作为该类型草原的月、年植被覆盖度值。草原类型的分区由《中华人民共和国植被图 1:100万》^[18]提供。研究区和样地位置如图 1所示, 基本情况特征和降水组合方式如表 1。

2 数据来源和方法 2.1 数据来源

采用MODIS产品(MCD43A4, MOD09A1, 2002—2016年), 使用线性分离方法将每个像元分为PV、NPV和BS 3种组分^[13], 得到空间分辨率为500 m的月植被覆盖产品数据集^[13]。分析该数据集中的PV覆盖度(f_PV)部分并对数据进行了质量检查。采用6—9月生长季f_PV值的平均值来表征内蒙古不同类型草原年植被覆盖度。

同期降水数据来源于中国气象数据共享服务网(http://data.cam.gov.cn)。分别选取距离研究样点最近的3个气象站, 即海拉尔、西乌珠穆沁旗、希拉穆仁的日降水数据, 年降水量是指1—12月的月降水量之和。

2.2 分析方法

(1)变化趋势

为研究植被覆盖度和降水量的变化趋势, 对数据和年份之间进行线性回归, 计算其拟合斜率b和相关系数r^[19]。如果斜率大于零, 说明要素变化趋于上升趋势, 反之则说明要素变化趋于下降。如果r值通过0.05的显著性水平(P < 0.05), 则认为要素呈显著增加或显著下降趋势。

(1)

(2)

式中, b为斜率；x_i 为第i年的降水量(植被覆盖度)；t_i 为第i年份, 在2002—2016间取值；n为研究时间段, 取15。

(2)相关性分析

植被覆盖度与降水量之间相关系数的计算公式如下：

(3)

式中, R_xy 为两个变量的相关系数, x_i 为第i年/月的植被覆盖度, y_i 为第i年/月的降水量, x 为多年/月植被覆盖度的平均值, y为多年/月降水的平均值, n为样本数。并建立内蒙古草原f_PV与降水的回归模型, 利用F检验对回归模型精度进行验证。

3 结果分析 3.1 光合植被覆盖度f_PV和降水量的年际变化特征

内蒙古草原生长季2002—2016年多年平均植被覆盖度, 以草甸草原样地最高, 为46.5%, 其次为典型草原, 为36.3%, 最低为荒漠草原22.4%。对3个草原区f_PV时间序列进行分析发现, 草甸草原f_PV随年际变化呈不显著增长趋势, 年变化斜率为0.29%/a, 而典型草原和荒漠草原f_PV呈不显著下降趋势, 年变化斜率分别为-0.04%/a和-0.21%/a, 如图 2所示。

2002—2016年3个草原区多年平均降水量分别为341.0、299.7、295.4 mm, 降水量随年际变化都呈现不显著波动上升趋势。其中海拉尔的年降水量增长速率最大, 年变化斜率为8.5 mm/a, 希拉穆仁的降水量上升趋势最弱, 年变化斜率为2.0 mm/a, 西乌珠穆沁旗的年降水量变化居中, 为6.1 mm/a。

3.2 光合植被覆盖度f_PV和降水量的年内变化特征

图 3是各草原区样地月均植被覆盖度变化曲线(源数据缺少11、12月和1月数据)。受季风气候影响, 各草原区月均植被覆盖度均表现出2—3月份最低, 3—7月持续上升, 7月均达到植被覆盖度峰值, 分别为61.4%、50.2%、25.0%。草甸草原和典型草原植被覆盖度在7—10月呈现快速下降过程, 而荒漠草原7—9月植被覆盖度缓慢降低, 10月份植被均处于荒芜状态, f_PV值小于10%。

从3类草原区的多年月平均降水曲线(图 4)明显看出, 降水量最大的月份均为7月份, 草甸草原区的海拉尔站7月降水量(96.4 mm)明显大于典型草原区的西乌珠穆沁旗站降水量(70.7 mm)和荒漠草原区的希拉穆仁站降水量(66.8 mm)。3类草原的夏季(6、7、8月)降水占全年的比例分别为63.5%、61.7%、60.0%左右。3类草原区月均降水量的分布格局与月植被覆盖度的分布十分相似。

3.3 月植被覆盖度f_PV对降水的响应

对3类草原区月植被覆盖度值与月降水量进行相关性分析, 计算生长季逐月f_PV与当月和前几月降水的相关系数, 探讨月降水量对月植被生长影响的滞后效应, 前期累积降水量对月植被生长的影响程度, 如表 2、表 3所示。

从月植被覆盖度与逐月降水的相关系数(表 2)可知, 内蒙古3种类型草原在6—9月的当月植被覆盖度与前一个月、当月降水量关系最为密切；5月植被覆盖度受降水的影响, 仅在荒漠草原与5月降水达到显著相关, 在草甸草原和典型草原不显著。表 2还可以看出, 3类草原月植被覆盖度对月降水量的响应均存在滞后效应, 尤其1—2个月的滞后效应最明显, 且表现出越干旱的草原类型滞后效应越明显的特征。

从月植被覆盖度f_PV与前期累积降水量的相关系数(表 3)看出, 内蒙古草原月累积降水对植被的作用在各月间以及不同类型草原区差异较大。5月植被覆盖度受累积降水的影响在各类型草原均不显著, 整体以6至10月的植被覆盖度与月累积降水相关性较好, 表现为越干旱的地区越受更长时间的前期降水影响的特征。总体上, 草甸草原植被覆盖度与当月+前1个月、当月+前2个月降水量最密切, 典型草原植被覆盖度与当月+前2个月、当月+前3个月降水量最密切, 荒漠草原植被覆盖度与当月+前3个月、当月+前4个月降水量最密切。

从草原类型来看, 草甸草原和典型草原植被覆盖度与累积降水相关性高的月份主要是6、7、8月；荒漠草原在6至10月植被覆盖度与累积降水相关性整体都较好, 同时, 荒漠草原植被覆盖度与累积降水的相关性较草甸草原和典型草原显著, 因为荒漠草原区降水量低, 植被对降水更加敏感。

总的来看, 5月的植被覆盖度和降水量, 除了典型草原与前第2个月的降水量、荒漠草原植被覆盖度与当月降水量为显著相关外, 其余均无显著相关, 受累积降水量的影响在各类型草原也均不显著。说明研究区5月份植被覆盖度开始增加, 与当月降水量、前期各月降水量、累积降水量关系不明显, 但有草原类型越干旱, 受前期降水量、累积降水量影响越大的趋势；内蒙古草原6—9月生长季植被对降水的响应较为强烈。

3.4 年植被覆盖度f_PV对降水的响应

用生长季6—9月植被覆盖度均值(f_PV)代表年植被覆盖度, 探讨年植被覆盖度与生长季各月、不同月累积降水量间的响应密切程度, 以及回归关系(表 4)；年植被覆盖度与当年、年不同时段降水量间的响应密切程度, 及其回归关系, 如当年降水量、生长季(6—9月)降水量、前期(11—5月)降水量, 结果如表 5。

从年植被覆盖度(6—9月植被覆盖度均值, f_PV)与生长季各月(6—9月)降水量间的响应程度来看(表 4), 单月降水因子对草甸草原和典型草原植被覆盖度主要是与7月降水量相关性最好, 7、8月降水对荒漠草原植被影响最大, 其相关系数为显著程度。将生长季降水量进行逐月累积后发现, 3类草原年植被覆盖度与生长季各月组合降水的相关系数都达到显著甚至极显著, 总体上与6—8月总降水量最相关。从年植被覆盖度与生长季各月降水量间的多元关系来看, 随着月份的增加, 3类草原年植被覆盖度与各月的关系也越密切, 与6、7、8、9月的多元关系最优, 其方程决定系数分别为0.746、0.407、0.818。说明3类草原年植被覆盖度与整个生长季降水都相关, 尤其是6—8月的总降水量的影响最显著。

将年植被覆盖度与当年降水量、生长季(6—9月)降水量、前期降水量(前一年11月至当年5月)、前一年降水量等要素进行回归(表 5), 发现年植被覆盖度与当年降水量、生长季降水量的关系都达到极显著相关。内蒙古草原为长冬短夏, 冬季开始时间从当年10月初到次年4月底, 持续7个月时间, 夏季为7月底到8月初的半月时间。表 5中进行了植被覆盖度与前期降水(与冬季时间几乎一致)的相关分析, 发现不同草原类型年植被覆盖度与冬季降水量均为正相关关系, 但在草甸草原和典型草原均不显著, 说明在草甸草原区和典型草原区冬季降水量可以促进植被覆盖度的增加, 但作用不显著。而对较干旱的荒漠草原植被覆盖度起到显著促进作用。将各时段降水因素与年植被覆盖度进行回归分析结果显示, 草甸草原和典型草原植被覆盖度均与当年、前期降水的回归效果最好, 而荒漠草原植被与生长季、前期降水拟合结果最优。

总体上, 内蒙古3种类型草原年植被覆盖度对降水量的响应均表现出在年、季、月尺度上分别受当年降水量、生长季降水量以及6、7、8月份降水量显著影响的特征。

4 讨论

从样地光合植被覆盖度以及降水年际变化趋势可以看出(图 2), 不同类型草原区植被覆盖度与降水量均具有很好的相关关系, f_PV值变化曲线都与降水量变化曲线波动情况基本一致, 也发现某些年份植被覆盖度并不能与降水量很好地吻合, 可能是草原植被覆盖度对降水存在一定滞后效应或相对于降水的变化而言植被结构更趋于稳定。

对3个草原多年f_PV月平均值与多年月平均降水量进行相关性分析发现：草原区月植被覆盖度对降水量的响应存在时滞效应, 这与穆少杰、张戈丽等人^{[3, 20]}研究结果基本一致。不管是在年尺度还是月尺度上, 相比草甸草原区和典型草原区, 荒漠草原植被覆盖度对降水更加敏感, 这与许旭研究结果一致^[11]。

本文对植被覆盖度与降水的相关性做了分析, 比较细致地探讨了植被覆盖度受不同时间降水影响的规律。以往文献^[21-22]很多都认为年际尺度上, 植被覆盖度主要受降水影响, 而月季尺度上, 植被受温度和降水综合因素影响大于单因子的影响。作为重要的气象因素, 温度对于植被的影响不容忽视, 但经过分析, 内蒙古3种类型草原的5—10月逐月植被覆盖度, 与前期月均温, 年植被覆盖度与生长季、非生长季的均温在多数情况下关系不大。也就是说, 对于内蒙古草原植被, 降水量对植被覆盖度起着决定性的作用^[23-24]。这与孟梦^[25]研究结果较吻合, 她认为内蒙古NDVI与气温和降水之间的关系表现出空间差异, 17.6%的面积与降水显著相关, 仅有0.4%的区域与气温显著相关, 降水对NDVI的影响比气温更为显著。因此, 本文以降水量作为主要气象因素, 分析了其对植被覆盖度的影响。

文中重点研究了光合植被覆盖度对降水的定量响应关系, 可为土壤侵蚀动态监测中消除降水波动的影响提供理论依据和具体算法。在地表覆盖对土壤侵蚀的影响中, 非光合植被覆盖度也起着重要作用, 需要进一步研究。总植被覆盖度(光合植被覆盖度+非光合植被覆盖度)需要根据具体研究区不同地物特征的光谱反射机理出发改进像元三分模型, 增强非光合植被覆盖度反映地物信息的准确程度, 进一步定量化光合植被覆盖度、非光合植被覆盖度和裸土以及三者的交互作用, 为区域土壤侵蚀动态监测与评估提供依据。

5 结论

本文通过对内蒙古3种草原类型2002—2016年间的光合植被覆盖度和同期降水进行分析, 得出如下结论：

(1) 2002—2016年间草甸草原、典型草原和荒漠草原多年平均光合植被覆盖度f_PV分别为46.5%、36.3%和22.4%；草甸草原f_PV随时间变化呈不显著增长趋势(斜率为0.29%/a), 典型草原和荒漠草原f_PV呈不显著下降趋势(斜率分别为-0.04%/a和-0.21%/a)。

(2) 2002—2016年3个草原区多年平均降水量分别为341.0、299.7、295.4 mm, 降水量随年际变化都呈现不显著波动上升趋势, 年内降水差异较大。

(3) 内蒙古草原的月植被覆盖度对月降水量存在明显的1—2个月滞后效应和显著的累积效应, 且表现出越干旱的草原类型滞后效应越明显的特征；相比草甸草原和典型草原, 荒漠草原植被对降水量变化更加敏感。

(4) 3类草原年植被覆盖度对降水量的响应, 均表现出在年、季、月尺度上分别受当年降水量、生长季降水量以及6、7、8月份降水量显著影响的特征；3类草原年植被覆盖度与生长季降水线性拟合结果都较好。