太行山区植被NPP时空变化特征及其驱动力分析*

李晓荣, 高 会, 韩立朴, 刘金铜**



太行山区植被NPP时空变化特征及其驱动力分析*

李晓荣1,2, 高 会1,2, 韩立朴1, 刘金铜1**

(1. 中国科学院遗传与发育生物学研究所农业资源研究中心 石家庄 050022; 2. 中国科学院大学 北京 100049)

本文基于2000—2014年MODIS NPP数据, 结合同期土地利用变化、气温、降水和DEM数据, 运用趋势分析法、相关系数法及分区统计法等方法, 研究了太行山区2000—2014年植被NPP时空变化特征, 分析了气温、降水等气候因素和人为因素对植被NPP变化的影响, 为太行山区植被资源管理及生态环境调控提供参考。研究结果表明: (1)太行山区植被NPP多年平均值为284.0 g(C)·m-2·a-1, 耕地、林地和草地的NPP均值分别为302.5 g(C)·m-2·a-1、258.1 g(C)·m-2·a-1、286.5 g(C)·m-2·a-1。(2)2000—2014年太行山区植被NPP整体呈上升趋势, 但大部分植被NPP变化未达到显著水平; 16.17%的植被NPP显著或极显著升高, 主要分布在太行山区西侧; 0.88%的植被NPP显著或极显著降低, 零散分布在研究区内。(3)不同植被类型NPP变化速率为草地>耕地>林地。(4)基于区域平均计算, 太行山区植被NPP与降水显著正相关(<0.05), 与气温负相关(>0.05)。基于像元计算, 植被NPP与降水显著或极显著正相关区面积比例为23.82%, 主要分布在太行山区北段, 几乎没有显著负相关区; 植被NPP与气温显著或极显著负相关区面积比例为8.42%, 主要分布在太行山区西侧, 显著或极显著正相关区面积比例为0.81%, 主要分布在太行山区最北端。(5)研究期内气候因子对植被NPP的升高整体上表现为促进作用, 而人为因素主要表现为抑制作用。太行山区生态环境保护仍应以减少人为干扰为主。

太行山区; 净初级生产力(NPP); 土地利用类型; 植被变化; MODIS; 驱动力

植被净初级生产力(net primary productivity, NPP)是指植被在单位时间、单位面积上由光合作用产生的有机物质总量扣除自养呼吸后的部分[1], 是反映植被生长变化的重要指标。植被NPP作为生态系统物质循环和能量流动的关键环节, 是生物圈的功能基础[2]。当前, 气候和环境问题日益突出, 及时准确了解区域或全球植被NPP, 并分析其时空变化特征, 理解其与各种影响因素之间的关系对于全球变化及可持续发展研究有重要意义, 同时也能为区域合理开发利用自然资源提供科学依据。

从20世纪80年代中后期开始, 学者们针对NPP经验统计模型的不足, 在考虑植物生理生态特征的基础上, 提出了植被NPP过程模型; 遥感(RS)、地理信息系统(GIS)和全球定位系统(GPS)技术的发展则为NPP过程模型的广泛应用提供了重要的技术平台[3]。20世纪末期以来, NOAA/AVHRR、EOS/ MODIS等遥感数据的引入, 有力推动了我国植被NPP遥感模型的研究和相关数据产品的应用。孙睿等[4]、朴世龙等[5]、陈利军等[6]、朱文泉等[7]先后对我国植被NPP时空变化特征进行了研究; 并且随着相关遥感数据的积累, 长时间序列遥感产品在区域植被变化研究中发挥的作用越来越大。

太行山区是京津冀地区的重要生态屏障, 又是生态敏感区, 生态环境极易受到各种因素的干扰; 历朝历代的开发对太行山天然植被造成了严重破坏, 使该区域水土流失严重, 旱涝灾害频繁, 影响了当地社会的可持续发展[8]。解放后, 对太行山区多半个世纪的生态恢复绿化了荒山秃岭, 产生了显著的综合效益[9]。当前气候变化和市场经济背景下, 中国陆地植被波动增大, 影响了区域生态系统服务功能的正常发挥。目前, 有很多植被变化方面的研究涉及太行山地区, 但相关研究中以NPP为中心的还很少, 且这些研究都只涉及了太行山区的部分区域。如刘征等[10]对河北省山区植被NPP空间格局进行研究, 发现NPP在河北省太行山区呈带状分布; 张莎[11]对河北省植被NPP时空变化特征的研究中发现太行山迎风坡在生长季降水较多, 形成NPP高值带; 王新闯等[12]对河南省NPP的研究中发现太行山河南部分为相对NPP低值区。目前, 将整个太行山区的植被作为研究对象的还鲜见报道。

太行山地处我国第二阶梯的东沿, 是平原到山地高原的过渡带, 同时, 它又是东部经济发达地区到中西部经济欠发达地区的过渡带。在太行山整体尺度上研究植被NPP变化, 可以增进人们对全球变化背景下生态敏感地区植被反馈特征的认识, 这对于典型山区自然资源的合理开发利用、社会的可持续发展有重要的现实意义。本研究借助3S技术对太行山区植被变化进行研究, 了解植被变化趋势, 并对导致植被变化的驱动力进行分析, 可以帮助管理者改进现有的植被资源管理策略, 更有效地调控太行山区的生态环境。

1 研究区概况

关于太行山区范围的说法有很多, 本文参考范晓[13]的说法, 将太行山区作为和西部吕梁山区、东部华北平原对等的地理单元; 同时参考古代“太行八陉”的说法[14], 确定太行山南北界。这样, 太行山自然区域北起永定河, 南止黄河, 西自汾河谷地, 东至华北平原, 跨越北京、山西、河南、河北4省(市), 115个县级行政区。本文研究范围(110°14′E~116°35′E, 34°34′N~40°47′N)特指太行山过境的所有县级行政区, 总面积约13.7×104km2, 平均海拔861 m, 地势西北高东南低。区内气候属温带大陆性季风气候, 四季分明, 年均温11.1 ℃, 年均降水511.7 mm, 全年一半以上的降水集中在夏季。研究区内植被状况随海拔、经度、纬度的变化而呈现出显著差异, 沟谷、盆地内平坦地区多为农田或建设用地, 崎岖陡坡地段多为林地或草地(图1)。

2 材料与方法

2.1 数据来源及处理

DEM数据来源于中国科学院计算机网络信息中心地理空间数据云平台(http://www.gscloud.cn), 为日本METI和美国NASA联合研制的ASTER GDEM V2数据, 基于“先进星载热发射和反辐射计(ASTER)”数据计算生成, 空间分辨率30 m。

植被NPP数据为美国蒙大拿大学发布的2000— 2014年的MOD17A3数据(https://ipdaac.usgs. gov), 空间分辨率为1 km, 时间分辨率为1 a。该产品基于MODIS/TERRA卫星遥感数据, 通过参考BIOME- BGC模型计算得到植被NPP数据。BIOME-BGC模型通过模拟植被生理过程来估算NPP, 估算精度相对较高[15], MODIS NPP数据的相关检验已有很多[16-17]。用MRT (MODIS reprojection tools)软件对图像进行拼接、转投影、重采样等操作。用ArcGIS 10.2软件进行有效值提取。

土地利用数据基于2010年TM影像、30 m DEM、Google高分影像, 进行目视解译生成, 数据来源于中国科学院山地灾害与环境研究所。

气象数据为2000—2014年太行山区及周边北京、山西、河南、河北4省(市)58个气象站点逐年的降水、气温数据, 数据来源于中国气象科学数据共享服务网(http://cdc.cma.gov.cn)。气象数据采用反距离加权插值法(inverse distance weighted, IDW)进行空间插值。

所有图像数据均采用WGS84/Albers Equal Area投影, 由ArcGIS 10.2软件进行进一步的处理分析。

2.2 分析方法

2.2.1 植被NPP变化趋势分析

用一元线性回归分析法, 基于像元分析太行山区2000—2014年植被NPP变化, 线性回归的趋势斜率即为NPP变化速率[16,18], 公式为:

式中:为NPP变化速率,=15为年数,x为第年的NPP值,为年份。基于ArcGIS的栅格计算功能, 可得到植被NPP在2000—2014年间的变化趋势图。当>0, 说明植被NPP呈增加趋势, 反之为减少趋势。

用检验法检验NPP变化趋势的显著性, 该检验结果仅代表趋势变化可置信程度的高低, 与变化速率大小无关[16,19]。检验公式为:

2.2.2 植被NPP变化与气候因素的相关分析

采用基于像元计算的偏相关系数分析NPP与降水或气温的相关性[16]。先计算简单相关系数, 进而计算偏相关系数。简单相关系数计算公式如下:

式中:R为、两变量的相关系数;x和y分别表示第年NPP和气候因子值;和分别表示多年平均NPP和多年平均气候因子值;=15, 为年数。

偏相关系数计算公式如下:

用检验法对偏相关系数进行显著性检验[20], 公式如下:

式中:为偏相关系数,=15为年数。查询分布表可判断偏相关系数的显著性, 结合偏相关系数正负情况, 可将其分为5个等级: 极显著负相关(<0,<0.01)、显著负相关(<0, 0.01<<0.05)、相关性不显著(>0.05)、显著正相关(>0, 0.01<<0.05)和极显著正相关(>0,<0.01)。

2.2.3 区域统计分析

将研究区分别按海拔梯度和坡度梯度分区, 结合土地利用数据、多年平均NPP、NPP变化分析结果、NPP与降水或气温的相关分析结果, 采用ArcGIS分区统计的方法, 研究不同区域植被NPP变化特征, 分析人为因素和气候因素对植被NPP变化的影响。

3 结果与分析

3.1 植被NPP的空间格局

太行山大体呈东北-西南走向, 跨越6个纬度, 山体东西两侧受季风影响差异较大, 且太行山区局部地形复杂多变, 因此区域内生态环境条件差异较大, 植被NPP空间分布异质性较高。太行山区植被NPP多年平均值为284.0 g(C)·m-2·a-1。其中, 耕地、林地和草地的NPP均值分别为302.5 g(C)·m-2·a-1、258.1 g(C)·m-2·a-1、286.5 g(C)·m-2·a-1。从整个区域来看(图2), 太行山区植被NPP空间分布呈现出北段低值、南段高值、中段大部分区域高低值交叉分布的特征。这一方面与研究区整体尺度上水热条件自北向南变好有关, 另一方面也与不同地理区域的水分条件差异有关。受季风气候和地形影响, 太行山区一般迎风坡和中高山区降水多、沟谷盆地降水少; 且地形对太阳辐射的分配使得阳坡比阴坡的蒸散量高, 植被干旱胁迫更严重, 因此, 当地人工造林多植于阴坡。耕地条件较高, 植被长势一般可以比同生长型自然植被更好。部分中高山NPP红色低值区分布突出, 对比土地利用图可发现, NPP低值区主要为林地。对比其他研究中林地NPP值为500 g(C)·m-2·a-1[16], 甚至更高的生产水平[21], 太行山区林地的生产力较差。结合太行山区林地NPP整体较低的特点, 太行山区林地植被生长受到很大的环境胁迫。主要由于太行山区林地分布的平均海拔为1 127 m, 平均坡度达21.4°, 年均降水507.7 mm, 自然条件比农田和草地的差; 另一方面, 山地生态系统比平原区更脆弱, 历史开发引起的水土流失使山体土层瘠薄, 很多地方土层厚度小于15 cm[22], 不少地方岩石大面积裸露, 只适合生长杂草, 这严重限制了植被的发育。

3.2 植被NPP年际变化分析

2000—2014年, 太行山区植被NPP整体呈波动上升的趋势, 不同植被类型NPP年际波动与研究区NPP整体波动基本一致(图3)。这说明影响太行山区植被NPP整体变化的因素是一种能作用于大尺度区域的全局性因素。2001年, 中国北方地区遭遇干旱, 研究区NPP值较低; 2003年, 研究区降水丰沛, 形成NPP高值, 并改善了来年植被长势。2000年左右, “退耕还林工程”启动, 封山禁牧等生态保护措施, 以及农村人口向城市转移等进程, 对植被改善产生了积极影响。

在植被NPP变化速率空间分布图中可看出 (图4a), NPP变化速率大部分介于-5~5 g(C)·m-2·a-1, 正值区面积远大于负值区; NPP增长较快的区域主要位于太行山西侧, NPP降低区除部分在太行山区南段林地成片分布外, 大多零散分布在各地。不同植被类型NPP变化速率有较大差别, 草地、耕地和林地的NPP变化速率平均值分别为2.77 g(C)·m-2·a-1、2.19 g(C)·m-2·a-1、0.83 g(C)·m-2·a-1。从NPP变化显著性检验结果(图4b)中可知, 研究区内大部分植被NPP变化不显著, NPP显著或极显著升高区面积比例为16.17%, 主要分布在太行山区西侧; NPP显著或极显著降低区面积比例只有0.88%, 零散分布在整个研究区域内。不同植被类型NPP显著或极显著升高的比例: 草地>耕地>林地, 不同植被NPP显著或极显著降低的比例: 林地>耕地>草地(表1)。这除了与林地生境较差有关外, 也与不同植被类型的生态学特征有关。林地生态系统结构复杂, 稳定性相对较高, 在环境变化中植被波动往往较小; 而草地和耕地以草本植物为主, 对环境变化敏感, 尤其草地植物适应性强, 在环境改善后更容易有高的NPP变化速率值。林地植被NPP显著降低则可能与人类活动有较大关联。

从植被NPP变化显著性检验结果的空间分布来看, NPP显著或极显著升高区和NPP变化速率高值区分布一致, 说明该区植被在研究时段内有明显改善; NPP显著或极显著降低区面积占全区面积的比例极小, 呈零星分布, 说明NPP变化速率负值区植被大部分并未明显恶化。NPP显著或极显著升高区主要为草地和耕地; NPP显著或极显著降低区主要为建设用地、耕地和林地。其中, 建设用地植被NPP降低是因为城市化过程中占用了原有植被用地; 耕地NPP降低部分多分布在建设用地周边, 说明人为因素是耕地区植被恶化的主要影响因素; 林地NPP降低部分多位于太行山区南段, 大多与NPP低值区、NPP变化速率低值区重合, 说明这部分林地中存在长期作用的负面影响因素或机制。气候因素作用范围大, 但林地植被显著恶化区孤立分布, 说明这些区域的植被恶化不是气候因素导致的, 结合近些年太行山区南段旅游业快速发展的事实, 推测这可能是人为破坏的结果。

3.3 植被NPP与气候因素的相关分析

基于空间平均计算, 2000—2014年, 研究区降水与NPP显著正相关, 相关系数为0.530; 气温与NPP负相关, 相关系数为-0.280。由偏相关系数显著性检验结果(图5)可知, 降水与植被NPP主要表现为正相关关系, 降水与NPP显著或极显著相关区的面积比例为23.82%, 几乎全部表现为正相关, 主要分布在太行山北段; 气温与植被NPP的关系主要为负相关关系, 其显著或极显著负相关区面积比例为8.42%, 主要分布在太行山区西侧; 气温与NPP显著或极显著正相关区面积比例只有0.81%, 这部分植被主要是分布在受热量限制的太行山区最北端的耕地和林地。

图b中, A为变化不显著, B为显著降低, C为显著升高, D为极显著降低, E为极显著升高。In figure b, A means not significant change, B means significant reduction, C means significant increase, D means extremely significant reduction, E means extremely significant increase.

A: 极显著负相关; B: 显著负相关; C: 相关性不显著; D: 显著正相关; E: 极显著正相关。A: extremely significant negative correlation; B: significant negative correlation; C: not significant correlation; D: significant positive correlation; E: extremely significant positive correlation.

表1 2000—2014年太行山区不同植被类型NPP变化显著性检验结果面积百分比

不同植被类型NPP与年降水量、年均温的相关性检验结果与全区平均状况基本一致(图6)。对比NPP变化显著性检验结果可发现, NPP显著变化区与NPP-降水显著相关区、NPP-气温显著相关区有较多重叠部分, 说明一部分植被NPP的显著变化与气候因素有较强关联。

以上分析中, 降水对植被NPP的增加主要表现为促进作用, 气温对植被NPP的增加主要表现为抑制作用, 这两种作用形式正是干旱胁迫的表现, 说明干旱是限制太行山区植被NPP增长的重要原因。一些学者也认为太行山区植被面对的干旱胁迫较为严重[8,22-23]。但降水只在全区平均水平上显著影响NPP变化, 在像元尺度的分析中, 其显著影响区的面积比例很小, 温度对植被NPP的影响力相对更小。这说明太行山区植被NPP的变化原因是复杂的, 降水和气温只是多种影响因素中的两种。

A: 极显著负相关; B: 显著负相关; C: 相关性不显著; D: 显著正相关; E: 极显著正相关。A: extremely significant negative correlation; B: significant negative correlation; C: not significant; D: significant positive correlation; E: extremely significant positive correlation.

3.4 地理梯度上植被特征的变化

基于海拔梯度和坡度梯度进行分区统计, 太行山区不同植被类型的面积百分比、年均NPP、NPP变化速率、NPP-降水偏相关系数和NPP-气温偏相关系数均产生变化。面积百分比表现为: 随着海拔上升、坡度增大, 耕地面积比例下降, 草地和林地面积比例上升, 且林地在中、高海拔和大坡度地区所占比例较大(图7a, b)。说明林地和草地的地理条件一般比耕地差。NPP的变化表现为: 随着海拔上升、坡度增大, 林地和耕地NPP趋于下降, 且林地下降幅度更大; 草地NPP则变化不大(图7c, d), 表现出了对不同地理环境的适应性。NPP变化速率方面, 随海拔和坡度增加, 耕地和草地均呈现出了先升高后降低的趋势, 而林地波动下降(图7e, f)。这表明林地植被长势受地形影响更大。降水与NPP的相关性表现为: 随林地在海拔梯度上波动下降(图7g, h), 表明随着海拔上升降水对林地植被生长的影响减小; 草地在坡度梯度上波动下降, 表明随着坡度增大, 降水对草地植被生长的影响减小; 但其中又有很多不规则的变化情况, 表明地理梯度上降水对植被生长的影响是复杂的。气温与NPP的相关性表现为, 不同植被类型NPP与气温的相关系数在地理梯度上均表现出了先降低后升高的趋势(图7i, j)。说明在不同地理梯度上, 气温所起的作用不同, 其中可能发生了主导因子的转变。除林地外, 气候因素与NPP的相关系数曲线大多在中等海拔、中等坡度部位表现凸出, 这与NPP变化速率曲线表现相似。一般而言, 低海拔、小坡度地区, 人类生产活动强烈, 人为因素是植被变化的主导因素, 气候因子作用相对较弱; 高海拔、大坡度地区, 植被NPP变化可能更多受到土壤等生态系统内部因素的影响。这反映出气候因素作用相对强烈的地方, 植被生长得更好。而人为因素作用强烈的低海拔、小坡度地区, 人为因素降低了耕地和草地的改善幅度, 尽管人类重视林地植被建设, 但这些地区林地比例较小, 人为因素主要表现为抑制植被改善。

4 结论和讨论

4.1 结论

本文基于2000—2014年MODIS NPP数据, 结合土地利用、气温、降水和DEM数据, 运用趋势分析、相关系数及分区统计等方法, 研究了太行山区2000—2014年植被NPP时空变化特征, 并对影响植被NPP变化的因素进行了分析, 结论如下:

(1)2000—2014年, 太行山区植被NPP多年平均值为284.0 g(C)·m-2·a-1。其中, 耕地、林地和草地的NPP均值分别为302.5 g(C)·m-2·a-1、258.1 g(C)·m-2·a-1、286.5 g(C)·m-2·a-1。

(2)2000—2014年, 太行山区植被NPP整体有上升趋势, 但大部分植被NPP变化达不到显著水平, 16.17%的植被NPP显著或极显著升高, 0.88%的植被NPP显著或极显著降低; 植被NPP显著或极显著升高区主要分布在太行山区西侧, 植被NPP显著或极显著降低区零散分布在研究区内。

(3)不同植被类型NPP变化速率为草地>耕地>林地。草地适应性相对较强, 更多利用了环境变化的机会; 林地大部分分布区生境较差, 制约了NPP的升高。

(4)基于区域平均计算, 太行山区植被NPP与降水显著正相关(<0.05), 与气温负相关(>0.05)。基于像元计算, 植被NPP与降水显著或极显著正相关区面积比例为23.82%, 主要分布在太行山区北段, 几乎没有显著负相关区; 植被NPP与气温显著或极显著负相关区面积比例为8.42%, 主要分布在太行山区西侧, 显著或极显著正相关区面积比例为0.81%, 主要分布在太行山区最北端。

(5)研究期内气候因子对植被NPP的升高整体上表现为促进作用, 而人为因素主要表现为抑制作用。

4.2 讨论

本研究中, 地理位置、地形、开发历史和人为管理等通过水、热、土壤等因子影响太行山区植被NPP的空间分布, 但水分条件是NPP增长的主要限制因子。太行山区受季风气候影响, 雨热同期, 但潜在蒸发量多在1 600 mm左右(根据MODIS蒸散数据计算多年平均得到), 远大于实际降水, 易形成干旱胁迫; 长期水土流失使土层瘠薄, 水源涵养功能降低, 植被水分状况进一步恶化。耕地植被受人类照料, 生长较好; 草地植被NPP低于耕地, 但因其适应性较强, NPP高于林地。本研究中太行山区林地NPP较低, 这与林地生境较差有关, 但也与不适宜的造林活动有一定关系。自20世纪80年代开始, 太行山区相继实施了“太行山绿化工程”、“退耕还林工程”等多项生态工程。但各地大规模造林过程中, 除在宜林地造林外, 因为指标问题, 难免会在不适宜造林的地方造林。这一方面增加了造林成本, 另一方面也违背了“因地制宜”的原则, 影响了林木生长, 使林分成为低质低效林[24-25]。张莎[11]对河北省2001—2010年NPP的研究中发现, 河北省NPP绝大部分集中在200~400 g(C)·m-2·a-1, 且耕地NPP大于林地和草地, 与本研究一致。王新闯等[12]对2000—2010年河南省NPP的研究中, 河南省年NPP变化范围为308~430 g(C)·m-2·a-1, 耕地NPP在300~500 g(C)·m-2·a-1, 豫西和豫东南林区NPP多在500 g(C)·m-2·a-1以上, 其植被长势整体上优于太行山区植被。

本研究中, 太行山区植被大体呈好转趋势, 这与国内很多学者的研究结果一致。杜加强等[26]对新疆植被NDVI的研究中, 发现植被在1982—2012年间呈极显著增加趋势; 李辉霞等[27]基于NDVI对三江源地区的研究, 发现植被在2000—2014年间呈好转趋势; 王新闯等[12]对河南的研究也发现植被NPP在2000—2010年间以增加为主。但在其他区域的研究中, 也有一些区域植被未改善, 甚至恶化。张莎[11]对河北植被NPP的研究中发现, 2000—2014年间河北大部分地区植被呈减少趋势。有学者指出, 植被变化有明显的区域差异性, 与气候因子关系的空间异质性也较强[20]。因此, 即便在大的环境变化背景下, 不同地区植被的响应也可能有很大差异。

对于太行山区植被的变化, 不同学者的研究结果也有所不同。孙雷刚等[28]对太行山中部滹沱河流域植被的研究中发现, 1984—2014年植被覆盖度总体呈增加趋势, 但海拔相对较低的区域表现为植被退化。本文中该区域虽然也存在植被退化, 但显著退化区比例较小。考虑到孙雷刚等[28]的研究中包括了石家庄平原区(本文研究范围未包含), 且研究时段更长, 而该区域又是城市化较快的地区, 这种差异是很可能出现的。贾宝全[29]基于1987年和2009年两期TM卫星影像对北京植被覆盖变化进行研究, 发现植被质量以改善为主, 且认为植被变化的驱动力主要是人为生态工程因素。本文也认为封育禁牧等生态措施对植被改善有利, 在低海拔区植被变化受人为因素主导, 但本文中低海拔区人类行为多表现为负面影响, 而非正面影响; 本研究缺乏林业相关数据分析及对2000年之前植被状况的认识, 贾宝全[29]的研究中北京造林绿化在2000年后已趋于饱和, 且MODIS NPP数据在建设用地缺失, 这些都会导致结果的差异。Chen等[30]基于AVHRR/NDVI数据对海河流域植被的研究中发现, 1981—2000年, 北京NDVI显著减少; 这与贾宝全[29]的研究结果有一定冲突, 但两者时空尺度和研究方法有较大差别,这可能影响了最终的结果。从学者们的研究中可看到, 方法差异、时空差异都可能导致研究结果的不同。

目前, 本研究还基本停留在对MODIS NPP数据的应用分析上, 缺乏对实地状况的认识, 与社会经济数据的结合还很少。而对人类行为只是定性分析, 其实际应用意义将大打折扣, 这是后续研究需要改进的。MODIS数据本身存在一定的精度问题, 本土化过程中需要多方考虑。当前, 各种植被生产力模型精度和准度均有较大差别, 模型应用结果可比性差; 有研究认为基于多种算法的集合预估方法能有效提高模拟精度, 产生显著优于单个算法的模拟结果[31], 这可能是今后的一个发展方向。

[1] Lieth H, Whittaker R H. Primary Productivity of the Biosphere[M]. Berlin Heidelberg: Springer, 1975

[2] 周广胜, 张新时. 自然植被净第一性生产力模型初探[J]. 植物生态学报, 1995, 19(3): 193–200 Zhou G S, Zhang X S. A natural vegetation NPP model[J]. Acta Phytoecologica Sinica, 1995, 19(3): 193–200

[3] 孙力炜. 祁连山区植被净第一性生产力的时空分布特征及气候变化和人类活动的影响[D]. 兰州: 西北师范大学, 2013 Sun L W. Temporal and spatial variation of vegetable net primary productivity and influence of climate change and human activity in Qilian Mountains[D]. Lanzhou: Northwest Normal University, 2013

[4] 孙睿, 朱启疆. 气候变化对中国陆地植被净第一性生产力影响的初步研究[J]. 遥感学报, 2001, 5(1): 58–61 Sun R, Zhu Q J. Effect of climate change of terrestrial net primary productivity in China[J]. Journal of Remote Sensing, 2001, 5(1): 58–61

[5] 朴世龙, 方精云, 郭庆华. 利用CASA模型估算我国植被净第一性生产力[J]. 植物生态学报, 2001, 25(5): 603–608 Piao S L, Fang J Y, Guo Q H. Application of CASA model to the estimation of Chinese terrestrial net primary productivity[J]. Acta Phytoecologica Sinica, 2001, 25(5): 603–608

[6] 陈利军, 刘高焕, 励惠国. 中国植被净第一性生产力遥感动态监测[J]. 遥感学报, 2002, 6(2): 129–135 Chen L J, Liu G H, Li H G. Estimating net primary productivity of terrestrial vegetation in China using remote sensing[J]. Journal of Remote Sensing, 2002, 6(2): 129–135

[7] 朱文泉, 潘耀忠, 张锦水. 中国陆地植被净初级生产力遥感估算[J]. 植物生态学报, 2007, 31(3): 413–424 Zhu W Q, Pan Y Z, Zhang J S. Estimation of net primary productivity of Chinese terrestrial vegetation based on remote sensing[J]. Journal of Plant Ecology, 2007, 31(3): 413–424

[8] 刘欣, 葛京凤, 冯现辉. 河北太行山区土地资源生态安全研究[J]. 干旱区资源与环境, 2007, 21(5): 68–74 Liu X, Ge J F, Feng X H. Study on ecological security of land resources in Taihang Mountain Hebei[J]. Journal of Arid Land Resources and Environment, 2007, 21(5): 68–74

[9] 曾宪芷, 杨跃军, 张国红, 等. 对太行山绿化工程建设的思考[J]. 林业经济, 2010, (7): 52–54 Zeng X Z, Yang Y J, Zhang G H, et al. On the greening project of Taihang mountains[J]. Forestry Economics, 2010, (7): 52–54

[10] 刘征, 赵旭阳, 米林迪. 基于3S技术的河北省山区植被净初级生产力估算及空间格局研究[J]. 水土保持研究, 2014, 21(4): 143–147 Liu Z, Zhao X Y, Mi L D. Study on estimation and spatial pattern of vegetation net primary productivity in mountainous area of Hebei Province based on 3S[J]. Research of Soil and Water Conservation, 2014, 21(4): 143–147

[11] 张莎. 基于长时间序列多源遥感数据的河北省植被NPP时空分布研究[D]. 石家庄: 河北师范大学, 2015 Zhang S. Research on spatio-temporal distribution of vegetation net primary productivity in Hebei Province based on long time-series multi-source remote sensing data[D]. Shijiazhuang: Hebei Normal University, 2015

[12] 王新闯, 王世东, 张合兵. 基于MOD17A3的河南省NPP时空格局[J]. 生态学杂志, 2013, 32(10): 2797–2805 Wang X C, Wang S D, Zhang H B. Spatiotemporal pattern of vegetation net primary productivity in Henan Province of China based on MOD17A3[J]. Chinese Journal of Ecology, 2013, 32(10): 2797–2805

[13] 范晓. 太行山——高原向平原的转折很壮丽[J]. 中国国家地理, 2011, (5): 15–25 Fan X. Taihang Mountain — The turning from the plateau to the plain is very magnificent[J]. Chinese National Geography, 2011, (5): 15–25

[14] 王尚义. 刍议太行八陉及其历史变迁[J]. 地理研究, 1997, 16(1): 68–76 Wang S Y. A brief discussion on the eight Taihang passages and their vicissitudes in the history[J]. Geographical Research, 1997, 16(1): 68–76

[15] 李登科, 范建忠, 王娟. 基于MOD17A3的陕西省植被NPP变化特征[J]. 生态学杂志, 2011, 30(12): 2776–2782 Li D K, Fan J Z, Wang J. Variation characteristics of vegetation net primary productivity in Shaanxi Province based on MOD17A3[J]. Chinese Journal of Ecology, 2011, 30(12): 2776–2782

[16] 穆少杰, 李建龙, 周伟, 等. 2001—2010年内蒙古植被净初级生产力的时空格局及其与气候的关系[J]. 生态学报, 2013, 33(12): 3752–3764 Mu S J, Li J L, Zhou W, et al. Spatial-temporal distribution of net primary productivity and its relationship with climate factors in Inner Mongolia from 2001 to 2010[J]. Acta Ecologica Sinica, 2013, 33(12): 3752–3764

[17] 朱锋, 刘志明, 王宗明, 等. 东北地区农田净初级生产力时空特征及其影响因素分析[J]. 资源科学, 2010, 32(11): 2079–2084 Zhu F, Liu Z M, Wang Z M, et al. Temporal-spatial characteristics and factors influencing crop NPP across northeastern China[J]. Resources Science, 2010, 32(11): 2079–2084

[18] 吴珊珊, 姚治君, 姜丽光, 等. 基于MODIS的长江源植被NPP时空变化特征及其水文效应[J]. 自然资源学报, 2016, 31(1): 39–51 Wu S S, Yao Z J, Jiang L G, et al. The spatial-temporal variations and hydrological effects of vegetation NPP based on MODIS in the source region of the Yangtze River[J]. Journal of Natural Resources, 2016, 31(1): 39–51

[19] 阿斯姆古丽·阿纳耶提, 师庆东, 刘曼, 等. 长时间序列遥感数据回归分析支持下的植被变化分析[J]. 新疆大学学报: 自然科学版, 2014, 31(3): 341–344 Anayeti A, Shi Q D, Liu M, et al. The vegetation cover changes of regress analysis on using time-serial images of remote sensing[J]. Journal of Xinjiang University: Natural Science Edition, 2014, 31(3): 341–344

[20] 阿多, 赵文吉, 宫兆宁, 等. 1981—2013华北平原气候时空变化及其对植被覆盖度的影响[J]. 生态学报, 2017, 37(2): 576–592 A D, Zhao W J, Gong Z N, et al. Temporal analysis of climate change and its relationship with vegetation cover on the north china plain from 1981 to 2013[J]. Acta Ecologica Sinica, 2017, 37(2): 576–592

[21] 陈福军, 沈彦俊, 李倩, 等. 中国陆地生态系统近30年NPP时空变化研究[J]. 地理科学, 2011, 31(11): 1409–1414 Chen F J, Shen Y J, Li Q, et al. Spatio-temporal variation analysis of ecological systems NPP in China in past 30 years[J]. Scientia Geographica Sinica, 2011, 31(11): 1409–1414

[22] 李保国. 生态经济沟建设理论与技术——以太行山为例[M]. 北京: 科学出版社, 2015 Li B G. Theory and Technology of Ecological Economic Ditch Construction — A case of Taihang Mountain[M]. Beijing: Science Press, 2015

[23] 邓伟, 戴尔阜, 贾仰文, 等. 山地水土要素时空耦合特征、效应及其调控[J]. 山地学报, 2015, 33(5): 513–520 Deng W, Dai E F, Jia Y W, et al. Spatiotemporal coupling characteristics, effects and their regulation of water and soil elements in mountainous area[J]. Mountain Research, 2015, 33(5): 513–520

[24] 杨海青. 南太行山低山丘陵区主要灌木林生态服务功能及可持续经营研究[D]. 郑州: 河南农业大学, 2013: 5Yang H Q. The ecological services and sustainable management of the main shrubs in the hill areas of south Taihang Mountains[D]. Zhengzhou: Henan Agricultural University, 2013: 5

[25] 罗晓华, 何成元, 刘兴良, 等. 国内低效林研究综述[J]. 四川林业科技, 2004, 25(2): 31–36 Luo X H, He C Y, Liu X L, et al. A study of the domestic low-benefit stand[J]. Journal of Sichuan Forestry Science and Technology, 2004, 25(2): 31–36

[26] 杜加强, 贾尔恒·阿哈提, 赵晨曦, 等. 1982—2012年新疆植被NDVI的动态变化及其对气候变化和人类活动的响应[J].应用生态学报, 2015, 26(12): 3567–3578 Du J Q, Ahati J, Zhao C X, et al. Dynamic changes in vegetation NDVI from 1982 to 2012 and its responses to climate change and human activities in Xinjiang, China[J]. Chinese Journal of Applied Ecology, 2015, 26(12): 3567–3578

[27] 李辉霞, 刘国华, 傅伯杰. 基于NDVI的三江源地区植被生长对气候变化和人类活动的响应研究[J]. 生态学报, 2011, 31(19): 5495–5504 Li H X, Liu G H, Fu B J. Response of vegetation to climate change and human activity based on NDVI in the Three-River Headwaters region[J]. Acta Ecologica Sinica, 2011, 31(19): 5495–5504

[28] 孙雷刚, 郑振华. 基于RS的近30年滹沱河流域植被覆盖度动态变化研究[J]. 地理与地理信息科学, 2014, 30(6): 36–40 Sun L G, Zheng Z H. RS-based study on dynamic change of vegetation coverage in Hutuo River Watershed in the past 30 years[J]. Geography and Geo-Information Science, 2014, 30(6): 36–40

[29] 贾宝全. 基于TM卫星影像数据的北京市植被变化及其原因分析[J]. 生态学报, 2013, 33(5): 1654–1666 Jia B Q. Driving factor analysis on the vegetation changes derived from the Landsat TM images in Beijing[J]. Acta Ecologica Sinica, 2013, 33(5): 1654–1666

[30] Chen F J, Shen Y J, Hu Q L, et al. Responses of NDVI to climate change in the Hai Basin[J]. Journal of Remote Sensing, 2011, 15(2): 401–414

[31] 袁文平, 蔡文文, 刘丹, 等. 陆地生态系统植被生产力遥感模型研究进展[J]. 地球科学进展, 2014, 29(5): 541–550 Yuan W P, Cai W W, Liu D, et al. Satellite-based vegetation production models of terrestrial ecosystem: An overview[J]. Advances in Earth Science, 2014, 29(5): 541–550

Spatio-temporal variations in vegetation NPP and the driving factors in Taihang Mountain Area*

LI Xiaorong1,2, GAO Hui1,2, HAN Lipu1, LIU Jintong1**

(1. Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang 050022, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China)

Net primary productivity (NPP) is an important indicator for the condition of vegetation in a given region. Research on NPP is not only important for the management of vegetation resources, but also a key element of global change. Technological development (such as remote sensing, geographic information system and global positioning system) had created the conditions for the establishment of complex process-based NPP models. On this basis, global NPP products today that continuously release data for long periods of time (e.g., MODIS NPP data appeared) greatly enhance research on NPP of regional vegetation. The objective of this study was to analyze the spatio-temporal variations in NPP of vegetation in Taihang Mountain Area for the period 2000–2014 using MODIS NPP data. At the same time, the study investigated the effects of climatic factors (e.g., temperature and precipitation) and of human factors (e.g., farming) on the change in NPP of vegetation in the region. Also trend analysis, correlation coefficient and zonal statistics were used in the study to analyze the various dataset (e.g., LUCC, temperature, precipitation and DEM). Results showed that the average NPP values of the study area was 284.0 g(C)·m-2·a-1, while those of farmland, forest and grassland were302.5 g(C)·m-2·a-1, 258.1 g(C)·m-2·a-1and 286.5 g(C)·m-2·a-1, respectively. Geographical location, topography, development history and human management influenced the distribution of the NPP of vegetation in the study area. Poor geographical environment was the main reason for the small NPP of forest vegetation in the region. The NPP of vegetation generally showed an upward trend, but was not significant for most of study area. About 16.17% area had significantly or extremely significantly increased vegetation NPP, which was mainly in the west of study area. In another 0.88% area, NPP significantly or extremely significantly dropped. The order of the NPP change rate of different vegetation types was — grassland > farmland > forest. Grassland vegetation was more effective for environmental rehabilitation because it had a better adaptability to the local conditions. Based on calculated regional averages, the vegetation NPP was significantly positively correlated with precipitation (< 0.05), but negatively correlated with temperature (> 0.05). About 23.82% of the study area had a significantly or extremely significantly positive correlation between NPP of vegetation and precipitation, which mainly distributed in the northern section of Taihang Mountain. No significantly negative correlation was observed. Furthermore, about 8.42% of the study area had a significantly or extremely significantly negative correlation between vegetation NPP and temperature, and this was mainly in the west side of Taihang Mountain. In another 0.81% area, the vegetation NPP was significantly or extremely significantly positively correlated with temperature, which was mainly distributed in the extreme north of the study area. Also the rate of NPP change and the correlation coefficient between NPP and climatic factors was positively correlated with altitude and slope gradients — both relatively smaller at low altitude and small slope, where human activity intensity was relatively higher. The area with significantly or extremely significantly reduced vegetation NPP was across the study area, especially around construction lands, which was not as a result of climatic factors. It was therefore suggested that while climatic factors generally enhanced vegetation NPP, human factors mainly inhibited vegetation NPP in the study area during the study period.

Taihang Mountain; NPP; Land use type; Vegetation change; MODIS; Driving factor

10.13930/j.cnki.cjea.160780

Q149

A

1671-3990(2017)04-0498-11

2016-09-02

2016-10-11

Sep. 2, 2016; accepted Oct. 11, 2016

* 国家重点基础研究发展计划(973计划)课题(2015CB452705)资助

* Supported by the National Program on Key Basic Research Project of China (2015CB452705)

**Corresponding author, E-mail: jtliu@sjziam.ac.cn

**通讯作者: 刘金铜, 主要研究方向为生态系统可持续管理与生态工程。E-mail: jtliu@sjziam.ac.cn

李晓荣, 主要研究方向为山地生态系统。E-mail: lxrlvge@163.com

李晓荣, 高会, 韩立朴, 刘金铜. 太行山区植被NPP时空变化特征及其驱动力分析[J]. 中国生态农业学报, 2017, 25(4): 498-508

Li X R, Gao H, Han L P, Liu J T. Spatio-temporal variations in vegetation NPP and the driving factors in Taihang Mountain Area[J]. Chinese Journal of Eco-Agriculture, 2017, 25(4): 498-508