1982—2017年黑河流域植被指数动态及其对气候变化与生态建设工程的响应
1.
2.
Vegetation dynamics and their response to climate change and ecological protection projects in the Heihe River Basin from 1982 to 2017
1.
2.
Received: 2022-10-05 Revised: 2022-11-08
作者简介 About authors
韩子言(1997—),女,天津人,博士研究生,研究方向为干旱区社会生态系统的耦合机制E-mail:
关键词:
Keywords:
本文引用格式
韩子言, 蒙吉军, 邹易, 朱利凯.
Han Ziyan, Meng Jijun, Zou Yi, Zhu Likai.
0 引言
在气候变化和人类活动的共同作用下,全球陆地生态系统表现出不同程度的退化[1-2],并且受人类活动的影响较大[3-4]。特别是干旱和半干旱地区,因气候变暖、干旱加剧和人口增长的综合效应,未来生态退化、沙漠化风险将进一步加剧[5],阻碍社会和经济发展,并威胁到生态系统的稳定性。植被作为陆地生态系统的关键组成部分,通过调节物质和能量循环维持生态系统功能,提供各种生态系统服务[6-7],其退化对以农牧业为主的干旱和半干旱地区的持续发展将产生巨大影响[8]。在过去几十年,随着对植被退化、生态恢复的重视,许多国家开始实施生态建设项目,旨在不降低基本要求的前提下,实现社会经济发展和生态环境保护的双赢,改善人类福祉[9-10]。
归一化差异植被指数(Normalized Difference Vegetation Index,NDVI)具有算法简单、覆盖范围广等优势,已被广泛应用于环境和气候变化的研究[11-13]。人类活动对植被NDVI产生重要影响。一方面,耕地开垦、工业布局、城镇扩建等经济活动将导致NDVI下降[14];另一方面,生态建设工程如退耕还林、开荒造林等直接改变土地覆被类型,或者通过封山育林、森林抚育等对提高NDVI具有重要作用[15]。近年来,一些学者采用相关分析[16]、回归分析[17]研究气候变化和人类活动影响下植被动态变化。而植被对外界变化的响应是非线性和非稳态的,故此类研究存在一定误差。基于实验机理的过程模型研究植被变化的驱动机制虽然较为精确,但需大量基础数据,并且模型验证也受到限制[18]。考虑时间尺度效应和空间异质性,残差分析可通过模拟无人类活动干扰下的植被变化与实际植被的差异来间接衡量人类活动的影响,进而有效衡量人类活动对植被的动态影响[19]。另外,还有学者采用人工神经网络获取生态建设基准期和改变期的插值量化分析生态建设工程对植被变化的影响[20-21]。一般来说,生态建设工程的实施会促使土地覆被向林地或草地转化[22],因此,本研究从土地覆被角度出发,通过揭示植被动态变化分析生态建设工程的实施效益。
黑河流域位于中国西北内陆干旱与半干旱地带,是荒漠集中、植被稀少、土壤侵蚀、沙漠化强烈的典型生态脆弱区[23]。多年来,在气候变化的背景下,随着中游绿洲农业和城镇化建设对水土资源需求的激增,黑河流域生态系统的敏感性和脆弱性进一步加剧,尤其是下游出现绿洲萎缩、胡杨林枯死等生态退化现象。为改善流域生态状况,2000年后中国先后实施了包括退耕还林还草、生态补偿、重点防护林建设计划、野生动物保护和自然保护区发展计划以及“黑河甘蒙分水方案”等多项生态建设工程。截至目前,这些生态建设工程产生了积极的效应,并在一定程度上缓解了生态退化的趋势[24]。与此同时,随着黑河流域社会经济的快速发展,耕作放牧、城镇扩张等人类活动又会导致植被出现衰退趋势。因此,剥离气候影响,量化不同土地覆被类型下植被对人类活动的响应,对评估黑河流域生态建设工程的有效性具有重要意义。
目前,黑河流域植被面临气候变化、社会经济高速发展和生态建设工程的多重影响,区域生态恢复和生态系统的管理具有很大不确定性。为此,本文基于遥感数据,采用趋势分析法、残差分析法以及地统计方法探究生态建设工程实施前后黑河流域植被动态变化,揭示植被对气候变化和人类活动的响应,进而在一定程度上衡量生态建设工程的效果。研究结果有望增强对气候变化背景下人类活动对植被变化过程的理解,为干旱区生态恢复制定有针对性的生态系统管理措施提供科学依据。
1 研究区概况与研究方法
1.1 研究区概况
黑河流域位于中国西北内陆地区(37°50′—42°40′N、97°—102°E),涉及青海、甘肃和内蒙古三省区,跨越半干旱、干旱与极度干旱气候区,多年平均气温由北向南递减、多年平均降水量由北向南递增,是西北内陆重要的生态屏障区。研究区总面积为13.4万km2,地势南高北低(图1),上游祁连山区大部分为中高山地,植被覆盖度高,人烟稀少,是黑河的水源涵养与径流补给区;中游为河西走廊平原区,地势较为平缓,绿洲、荒漠、戈壁、沙漠断续相间分布,光热资源丰富,是黑河流域人口分布最为密集的区域,由于农业灌溉和城镇生产生活用水,也是黑河径流消耗的主要区域;下游额济纳盆地包括甘肃省金塔县和内蒙古额济纳旗,以戈壁沙漠为主(75%以上),生态环境极为脆弱。
图1
1.2 数据来源和处理
为延长对区域植被变化的研究时间,本研究采用的植被NDVI数据来源于GIMMS NDVI和SPOT/VEGETATION NDVI数据集。其中,GIMMS NDVI数据来自美国航空航天局(
研究使用的气候数据来自国家青藏高原科学数据中心(
1.3 研究方法
1.3.1 趋势分析
式中:slope为NDVI变化趋势,slope>0 表明NDVI呈增加的趋势,slope<0 表明NDVI呈减少的趋势;n为监测年数;j代表时间序列。本文根据上述显著性检验结果,将NDVI的趋势变化分为6类:极显著下降(P<0.01,slope <0)、显著下降(P<0.05,slope <0)、不显著下降(P>0.05,slope <0)、不显著增加(P>0.05,slope >0)、显著增加(P<0.05,slope >0)、极显著增加(P<0.01,slope >0)。
1.3.2 偏相关分析
偏相关分析可以反映植被NDVI与环境要素之间的相互关系[29]。当某一环境要素的变化引起其他要素变化的时候,偏相关分析可以将其他要素的影响剔除,计算公式如下:
式中:
1.3.3 残差分析
式中:
根据残差分析结果,分别得到气候变化和人类活动对NDVI的相对贡献率:
式中:
2 结果
2.1 NDVI 和气候因子的变化趋势
2.1.1 区域尺度 NDVI 和气候因子变化趋势
如图2所示,在流域尺度上,黑河流域年平均气温和年降水量呈现显著的增加趋势,年变化速率分别为0.0518 ℃和1.8915 mm,表明近40年黑河流域气候朝着温暖湿润的方向发展。NDVI呈波动上升趋势,NDVI为0~0.9,多年平均值为0.19。其中,1982—1999年NDVI上升速率为0.0009;2000年上升速度加快,达0.0018,是之前的2倍。
图2
图2
1982—2017年黑河流域植被指数(NDVI)和气温、降水量的年际变化趋势
Fig.2
Trends of NDVI (A), temperature (B) and precipitation (C) in Heihe River Basin during 1982-2017
2.1.2 栅格尺度的 NDVI 变化趋势
1982—2017年,黑河流域NDVI变化趋势表现出显著的空间异质性(图3C和图3F),呈现南高北低的态势。从表1看出,49.35%的区域NDVI得到改善,其中极显著增加占27.36%,显著增加占3.11%;50.65%的区域NDVI表现出降低趋势,其中极显著降低占22.13%,显著降低占4.88%, 主要分布在下游荒漠区。2000年前,除下游荒漠地区外,植被普遍呈现缓慢变绿趋势(图3A和图3D),71.86%的区域植被NDVI呈现增加趋势,其中极显著增加的占6.93%,显著增加的占36.11%;28.14%的区域植被NDVI呈现降低趋势,只有7.81%的面积植被NDVI呈现出显著降低趋势。2000年以后,研究区NDVI的变绿趋势加速(图3B和图3E)。除未利用地和建设用地,其他土地覆被类型均呈现NDVI增加趋势的区域大于减少趋势的区域(表2),95.73%的耕地NDVI得到较大改善,其次是林地和草地,占比分别为90.30%和64.13%。
图3
图3
不同阶段黑河流域植被NDVI变化趋势(A~C)及显著性检验(D~F)的空间分布
Fig.3
Spatial patterns of NDVI trend (A,B,C) and statistical significance (D,E,F) in Heihe River Basin during 1982-2017
表1 不同阶段 NDVI 变化趋势统计(%)
Table 1
阶段 | NDVI增加 | NDVI降低 | ||||
---|---|---|---|---|---|---|
极显著 | 显著 | 不显著 | 极显著 | 显著 | 不显著 | |
1982—1999年 | 6.93 | 36.11 | 28.83 | 1.80 | 6.01 | 20.33 |
2000—2017年 | 17.15 | 5.58 | 13.50 | 21.13 | 17.19 | 25.43 |
1982—2017年 | 27.36 | 3.11 | 18.88 | 22.13 | 4.88 | 23.64 |
表2 1982—2017年不同土地覆被类型 NDVI 变化趋势面积占比统计(%)
Table 2
土地覆被类型 | 耕地 | 林地 | 草地 | 建设用地 | 未利用地 |
---|---|---|---|---|---|
NDVI增加区域面积占比 | 95.73 | 90.30 | 64.13 | 33.33 | 34.18 |
NDVI降低区域面积占比 | 4.27 | 9.70 | 35.87 | 66.67 | 65.82 |
2.2 NDVI 对气候因子的响应
图4
图4
黑河流域NDVI与气候因子之间的偏相关性(A、B)及显著性检验(C、D)的空间分布
Fig.4
Spatial patterns of partial correlation coefficients (A,B) and statistical significance (C,D) between NDVI and climatic factors in Heihe River Basin
1982—2017年,气候变化对植被变异改善的贡献率为75.70%,对植被退化归因占82.97%(图5)。
图5
图5
气候变化对黑河流域NDVI的相对贡献率的空间分布
Fig.5
Spatial pattern of the relative contribution of climate change to NDVI in the Heihe River Basin
2.3 NDVI 对生态建设工程的响应
图6
图6
黑河流域1982—2017年NDVI残差趋势(A)及显著性检验(B)的空间分布
Fig.6
Spatial distribution of NDVI residual trends (A) and statistical significance (B) in the Heihe River Basin
人类活动对NDVI的相对贡献与气候变化的相对贡献率具有明显的空间差异性(图7)。1982—2017年,共有24.30%植被改善归因于生态建设工程的实施;而基础设施建设、过度放牧、耕地开垦等人类活动对植被的负向影响占17.03%。此外,78.12%的区域受气候变化的影响大于人类活动的影响。
图7
图7
人类活动对黑河流域NDVI的相对贡献率的空间分布
Fig.7
Spatial pattern of relative contribution of human activities to NDVI in the Heihe River Basin
图8
图8
不同时期气候变化和生态建设工程对植被动态影响的空间分布(A~C)及不同土地覆被类型统计(D~F)
Fig.8
Spatial patterns of the distinguished impacts of climate change and ecological project on vegetation in different periods (A,B,C) and the summarized results by land cover (D,E,F)
为了进一步探讨气候变化和生态建设工程对NDVI的潜在影响,基于对应的土地覆被类型计算气候变化和生态建设工程的相对贡献率。在不同的土地覆被类型中,生态建设工程和气候变化对NDVI变化的相对贡献率不同(图8D~F)。1982—1999年,植被NDVI的改善主要受气候影响,尤其是耕地、林地和草地中气候变化产生的正向影响分别占88.13%、79.37%和56.75%,说明植被NDVI的改善主要归因于气候变化;2000—2017年,生态建设工程对植被的改善效果更为显著,在耕地和林地的表现尤其突出,生态建设工程带来的改善分别占75.55%和72.42%,此外,生态建设工程对建设用地植被改善也较为明显,有56.67%的建设用地受到生态建设工程的正向影响,同时有43.81%的草地受到生态建设工程的正向影响,未利用地中生态建设工程的影响最小,仅有11.23%的未利用地改善归因于生态建设工程。
3 讨论
通过遥感持续监测植被动态对生态恢复具有重要意义[32-33]。本研究发现,黑河流域近40年来NDVI在气候变化和人类活动的共同影响下总体呈波动上升趋势。这与Han等[34]发现黑河流域生态环境质量近年来明显改善的结论基本一致。黑河流域自然地理条件差异显著,近年来气候呈变暖和变湿的趋势,而气温和降水是影响植被分布的主要因素[35],呈现明显区域差异。根据偏相关结果可以发现,气候变化主导的NDVI正向影响以黑河中游地区最为明显,主要分布在中游绿洲的耕地与草地,其次是上游地区,高值区主要分布在上游草地,这与Spano等[36]发现干旱与半干旱区植被生长开始期和持续期主要受降水控制的结论一致。在干旱与半干旱区,虽然气候变暖可能会加剧土壤水分缺失,但气温升高也可提高植被光合效率,改善植被生长状况[37]。
气候变化和人类活动相互作用,共同影响干旱区植被生态系统的生长。一些研究发现,在气候因子的主导下,20世纪80年代以来植被NDVI呈现上升趋势,气温的升高和降水的增加延长了植被生长季节的长度并促进光合作用,促进了植被生长[38],与本研究发现1984—1999年NDVI的增加主要是由于气候变化主导产生的正向影响一致。中游和下游部分地区出现生态退化现象,则归因于人类活动,包括人口快速增长、农牧业发展和城市化导致水资源的过度消耗,使得植被恢复的生态需水受到影响,进而出现绿洲萎缩、胡杨林枯死的现象[39]。NDVI的残差分析表明人类活动对流域植被的影响程度呈显著增加态势。研究发现,2000年以来由于生态建设工程的实施,黑河流域上游以及中下游绿洲地带NDVI增长是显而易见的,与You等[23]发现黑河流域生态系统的退化在生态建设工程的实施下得以缓解的研究结果一致。2001年以来,甘肃祁连山国家级自然保护区先后实施天然林保护与生态公益林建设生态工程,上游林地、灌丛和草地等生态空间的大幅提升,对大部分地区植被状况产生正向影响。经过国家两期“三北防护林建设工程”的实施,中游的人工绿洲内部90%以上的农田均得到保护,通过建设完整的人工绿洲防护林体系庇护中游绿洲与农田的生态环境,生态建设工程对中游林地与耕地植被状况产生正向影响。“黑河甘蒙分水方案”则通过约束中游水资源利用来保障下游生态需水,使得下游绿洲NDVI呈现显著提高。Hao等[40]也发现下游绿洲区地下水位的显著提高促进了当地植被的生长。
本研究发现不同土地覆被类型对人类活动的响应不相同,表现为耕地、林地和建设用地较为明显。1980—2017年,不同土地覆盖类型的面积并未发生显著变化,说明植被恢复主要源于植被改良,而不是土地覆被类型的转换。特别是2000年后,林地变化和草地变化主要来自生态建设工程产生的正向影响,分别占79.37%和56.75%,说明政策驱动的大规模生态建设工程对黑河流域植被恢复具有积极影响。生态建设工程主要通过直接退耕还林或荒地造林,从耕地和荒地转入植被,或者通过间接天然林保护计划、封山育林增加原有植被这两种方式作用于生态系统。通过已有研究总结2000年以来中国实施的重大生态工程,其中发生在黑河流域的包括三北防护林体系建设工程、退耕还林工程、天然林资源保护工程(上游)、退牧还草工程(上游、下游)[41],分别对流域上游、中游和下游展开分析。可以发现,2000年以后上游地区生态建设工程对NDVI的正向影响明显提升,占20.11%,在草地的表现尤为突出(占14.53%),其次为林地(占3.03%),反映出上游退耕还林、天然林资源保护与退牧还草等工程对林草植被改善起到的积极作用。中游地区生态建设工程对NDVI的正向影响占17.60%,以耕地和草地最为明显,分别占7.67%与2.71%,此外,建设用地的改善也相对上游与下游较高,占0.07%,这是由于中游是耕地和建设用地集中分布区,耕地农作物变化周期较短,响应速度较快,高度依赖于灌溉活动[42];而在人口增长和城市化进程下,建设用地一直是对人类活动影响最强烈的土地利用类型之一[43]。下游地区的生态建设工程对NDVI的正向主要体现在草地(占5.72%),其次为林地(占0.38%),反映出退牧还草、天然林资源保护工程等对下游植被恢复的积极影响。本研究基于遥感监测分析区域土地覆被和植被覆盖的空间分布和变化,衡量生态建设工程的有效性,确保生态系统的可持续性。残差分析对比不同土地覆被类型对应NDVI对气候变化和生态建设工程的响应,未考虑气候变化和生态建设工程所带来的滞后效应,可能会产生一定误差。
4 结论
在生态环境极为脆弱的黑河流域,由于绿洲农业和城镇化发展对水土资源需求的激增,对植被产生了重要的影响。为改善生态系统功能,缓解植被退化,中国自2000年起在黑河流域实施一系列生态建设工程。研究发现,在气候变暖变湿和生态建设工程的双重影响下,1982—2017年黑河流域植被NDVI总体呈现波动上升趋势。尤其是2000年后,NDVI增长速率上升,且上升态势呈现南高北低的分布格局。气温和降水显著影响NDVI的变化,偏相关分析结果表明降水对NDVI的空间异质性影响更为显著。此外,残差分析结果表明黑河上游祁连山区、中游城镇区以及下游绿洲区存在显著的正残余趋势,受人类活动影响较大。特别是2000年后生态建设工程对植被变化的影响远超过气候变化,30.18%的研究区受到生态建设工程的正向影响,其中在耕地和林地的表现尤为突出。由此可见,生态建设工程对黑河流域植被变化起到积极影响。本研究旨在提高对黑河流域植被动态变化驱动力的理解,为干旱区生态系统管理和可持续发展提供理论指导。
参考文献
Millennium Ecosystem Assesment.Ecosystems and human well-being:synthesis
[J].
Summary for Policymakers of the IPBES Regional Assessment Report on Biodiversity and Ecosystem Services for Europe and Central Asia
[R].
Human-induced greening of the northern extratropical land surface
[J].
Global land change from 1982 to 2016
[J].
Mapping and modelling trade-offs and synergies between grazing intensity and ecosystem services in rangelands using global-scale datasets and models
[J].
Accelerating net terrestrial carbon uptake during the warming hiatus due to reduced respiration
[J].
Satellites reveal contrasting responses of regional climate to the widespread greening of Earth
[J].
Achieving a fit between social and ecological systems in drylands for sustainability
[J].
Variation in actual evapotranspiration following changes in climate and vegetation cover during an ecological restoration period (2000-2015) in the Loess Plateau,China
[J].
Increased plant growth in the northern high latitudes from 1981 to 1991
[J].
Changes in satellite‐derived vegetation growth trend in temperate and boreal Eurasia from 1982 to 2006
[J].
Global urban expansion offsets climate-driven increases in terrestrial net primary productivity
[J].
Drought reduces the effectiveness of ecological projects:perspectives from the inter-annual variability of vegetation index
[J].
Determining the contributions of climate change and human activities to vegetation dynamics in agro-pastural transitional zone of northern China from 2000 to 2015
[J].
Impacts of climate change and human activities on grassland vegetation variation in the Chinese Loess Plateau
[J].
Response of wind erosion dynamics to climate change and human activity in Inner Mongolia,China during 1990 to 2015
[J].
NDVI-based vegetation dynamics and its response to climate changes at Amur-Heilongjiang River Basin from 1982 to 2015
[J].
Impacts of ecological restoration projects on the ecosystem carbon storage of inland river basin in arid area,China
[J].
Isolating the impacts of land-use/cover change and climate change on the GPP in the Heihe River Basin of China
[J].
Dynamic analysis of ecological environment combined with land cover and NDVI changes and implications for sustainable urban-rural development:the case of Mu Us Sandy Land,China
[J].
Comparison and conversion of AVHRR GIMMS and SPOT VEGETATION NDVI data in China
[J].
Spatio-temporal analysis of vegetation variation in the Yellow River Basin
[J].
Discrimination between climate and human-induced dryland degradation
[J].
What drives the vegetation restoration in Yangtze River Basin,China:climate change or anthropogenic factors?
[J].
North American vegetation dynamics observed with multi-resolution satellite data
[J].
Evaluation of earth observation based long term vegetation trends-Intercomparing NDVI time series trend analysis consistency of Sahel from AVHRR GIMMS,Terra MODIS and SPOT VGT data
[J].
Quantifying trade-offs of land multifunctionality evaluated by set pair analysis in ecologically vulnerable areas of northwestern China
[J].
Phenological stages of natural species and their use as climate indicators
[J].
Northern hemisphere photosynthetic trends 1982-99
[J].
Changes in riparian woods over space and time:influence of environment and land use
[J].
Impacts of ecological water conveyance on groundwater dynamics and vegetation recovery in the lower reaches of the Tarim River in northwest China
[J].
Climatic change and human activities link to vegetation dynamics in the Aral Sea Basin using NDVI
[J].
/
〈 |
|
〉 |
