梁大林,唐海萍

北京师范大学地理科学学部自然资源学院, 北京 100875

青藏高原草地资源占全国总草地资源的1/3,是我国面积最大的天然草地[1—2],研究青藏高原不同草地植被变化对草地资源保护与调查具有重要意义。青藏高原也是我国重要的生态功能区,是全球气候变化最敏感的地区之一[3—4],其草地植被对气候变化响应研究一直是众多学者关注的焦点问题之一[5—9]。

上世纪以来青藏高原地区年平均气温总体呈现出升高的趋势,年降水量在近40年呈现增加趋势[10],这种以变暖变湿为主要特征的气候变化对青藏高原植被覆盖具有重要影响,高寒草地作为青藏高原主要植被类型之一,受气候变化影响较为显著[11—12]。植被NPP是衡量植物在自然环境条件下的生产能力,表征全球变化对陆地生态系统影响及相应关键参数[13—14],不少学者基于植被NPP研究植被变化[15—17],从时空尺度探讨植被NPP变化特征及其与气候变化之间的关系[18]。张镱锂等[19]利用CASA模型估算青藏高原高寒草地1982—2009年间植被净初级生产力,分析了高寒草地植被NPP时空格局与变化特征,对青藏高原草地植被NPP变化有整体认识；Zheng等[20]研究青藏高原2001—2015年植被NPP时空格局及其对气温、降水和太阳辐射的响应,从总体上对青藏高原植被NPP与气候因素进行分析；韩炳宏等[21]通过植被NDVI研究了青藏高原2000—2018年植被覆盖变化及其与气候因素的关系,得出青藏高原植被整体呈稳定恢复状态,与降水和温度呈正相关。前人在青藏高原植被NPP变化方面做了很多研究,取得众多研究成果,但对青藏高原不同地区具体植被类型NPP变化与气候变化响应差异研究相对欠缺。周秉荣等[22]研究发现三江源区近54年植被NPP整体呈显著增加趋势,气温是影响该地区植被NPP增加的主要气象因素,该研究只是针对三江源植被,气温是否也是影响青藏高原其他地区植被NPP增加的主要气象因素需进一步研究； 陈舒婷等[23]通过相关分析和重心模型等方法对2000—2015年整个青藏高原植被NPP时空变化格局及其驱动机理进行研究,从整体上区分气候变化和人类活动对植被变化的影响。而戴黎聪等[24]只是分析了青藏高原高寒小嵩草草甸和高寒金露梅灌丛两种植被NPP变化和气候相关性,对全球变暖背景下这两种植被NPP变化差异有进一步认识,而对于高寒草甸和高寒灌丛是否也成立,需要进一步研究。

现有研究主要集中在青藏高原整体植被NPP时空分布格局及其对气候响应,为青藏高原区域尺度上长时间序列植被NPP时空分异格局与变化规律研究奠定了基础[25],但对青藏高原内部不同区域、不同草地类型之间植被NPP与气候变化相关性差异研究较少,缺乏对青藏高原不同草地植被生长与气候制约因子比较研究。而高寒草甸和高寒草原作为青藏高原两种典型植被类型,是青藏高原高寒草地生态系统的主体[26—27],在不同地区、不同气候环境下,二者对气候变化的响应也可能存在显著差异,目前关于两者对比研究相对欠缺。

基于此,本文以三江源高寒草甸和阿里高寒草原为研究对象,利用研究区2000—2017年草地植被月最大值NPP变化表征植被变化(下文所提NPP均代表月最大值NPP),根据气象站点数据插值得到研究区气温、降水时空变化,研究两种高寒草地植被变化与气象因子相关性,得出气象因子对两种高寒草地植被变化影响差异,为青藏高原不同草地植被对气候变化响应研究提供参考。

1 研究区与研究方法

1.1 研究区概况

三江源地区位于我国青海省南部(图1),是长江、黄河和澜沧江的发源地,平均海拔3500—4800 m。该区域总面积达38.24×104km2,年平均气温低于-1℃,属于典型的高原大陆性气候,冷热两季交替、干湿季节分明,且植被类型以高寒草甸为主[28]。

阿里地区位于青藏高原北部西藏自治区西部(图1),平均海拔4500 m以上,素有“世界屋脊”的屋脊之称。该地区总面积约为33.7×104km2,年平均气温在0.4—3.6℃之间,属亚寒带干旱高原气候区,植被类型以高寒草原为主[29—30]。

图1 研究区采样点、气象站点和高程分布Fig.1 Sampling points, meteorological stations and elevation distribution in the study area

1.2 数据来源与预处理

本研究所采用数据和预处理过程如下

(1)研究区草地类型数据源自1∶100万植被图(图1)。(2)2000年3月—2018年2月降水、气温月值数据集均来自国家气象科学数据中心(http://www.nmic.cn/),包括青藏高原及周边地区102个气象站点(图1)；利用专业气象插值软件ANUSPLIN对各年月值气象数据插值,得到青藏高原2000—2017年250 m分辨率气温与降水栅格数据。ANUSPLIN插值是一种广泛应用于青藏高原气候因子模拟的方法,插值时加入高程协变量能较好地拟合青藏高原气象要素[7]。(3)三江源与阿里地区各县牛羊肉产量数据来源于中国知网中青海省统计年鉴(https://www.cnki.net/)。(4)NPP数据来源于国家青藏高原科学数据中心(https://data.tpdc.ac.cn/zh-hans/)青藏高原生态资产评估遥感反演净初级生产力数据集[31—32],空间分辨率250 m×250 m。本研究最后将所有数据重采样为250 m×250 m栅格数据。

1.3 NPP数据评估

本研究分析所用NPP数据均来自青藏高原生态资产评估遥感反演净初级生产力数据集,时间跨度2000—2017年,该数据集基于CASA模型计算得到植被月最大值NPP[31],代表植被一年中生长最好月份状况。为评估实验数据在三江源与阿里地区的可靠性,本研究基于Liang等人[33]发表文献中实测生物量数据评估三江源区数据准确性,利用国家青藏高原科学数据中心公开的藏北样带实测生物量数据评估阿里区数据准确性。通过对研究区NPP数据与实测生物量数据进行相关与线性回归分析,结果表明,三江源区高寒草甸植被NPP与实测生物量数据相关系数达0.863(P<0.001),线性回归拟合优度R2为0.749；阿里地区高寒草原NPP与实测生物量数据相关系数达0.833(P<0.05),线性回归拟合优度R2为0.695(图2)。由此可见,本研究所用产品中两种高寒草地植被NPP计算结果可靠,可进行下步相关研究。

图2 研究区NPP数据与实测生物量数据线性拟合Fig.2 Linear fit between NPP data and measured biomass data in the study area

1.4 研究方法

1.4.1Sen+M-K趋势分析

对NPP及气象数据趋势分析一般采用回归分析或Mann-Kendall(M-K)趋势分析,回归分析要求时间序列符合正态分布,且结果容易受噪声干扰[34]。相比之下Theil-Sen Median(Sen)趋势度分析是计算序列中的中位数,可以很好地减少噪声干扰,是一种稳健的非参数统计的趋势计算方法,但其本身不能实现序列趋势显著性判断。这两种方法结合使用在时间序列趋势分析中具有很大优势,不仅可以有效增强抗噪性,且更准确地检验时间序列趋势及其显著性。本文利用Sen趋势度对时间序列趋势进行判断,利用M-K趋势检验法对趋势显著性进行检验。

Sen趋势计算公式：

(1)

式中,xj和xi为时间序列数据,Median()表示序列中位数,当Sen>0表示时间序列呈现上升趋势；当Sen<0表示时间序列呈下降趋势[35]。

M-K趋势显著性检验方法如下：

(2)

(3)

式中,xj和xi为时间序列数据,n为数据个数；当n≥8时,(3)式检验统计量S近似为正态分布。

根据统计量S标准化检验统计量Z：

(4)

原假设H0表示时间序列(x1、x2、x3、…,xn)是独立同分布样本数据,不存在时间序列趋势,备择假设H1表示序列存在上升或下降趋势。在给定显著性水平α下,若|Z|≥Z1-α/2,则拒绝原假设H0,在置信水平α上采用备择假设,即时间序列数据存在显著上升或下降趋势。在本文检验中Sen>0表示呈上升趋势,Sen<0表示下降趋势；当Z的绝对值大于1.65、1.96和2.58时,表示趋势通过的信度分别为90%、95%和99%的显著性检验。

M-K方法还可以用于突变检验,假定有n个样本的时间序列x1,x2,…,xn,构造一秩序列：

(5)

式中,xi>xj时,ri=+1；xi

(6)

1.4.2相关与偏相关性计算

本文对阿里和三江源地区草地NPP与气象因子分别采用逐像元相关与偏相关分析。由于年降水量和年平均温度都对草地NPP有影响,控制其中一种气象要素不变,研究另一种气象要素对草地NPP的影响[36],计算原理见参考文献 [36]。

Pearson相关分析计算公式如下式(7)：

(7)

偏相关分析计算公式如下式(8)

(8)

式中,n表示时间序列长度,本研究时间序列为18,rxy表示x与y之间相关系数,rxy·z表示控制变量为z时,x与y之间偏相关系数[37]。

本研究数据处理使用了ANUSPLIN和Matlab软件, 结果分析与统计地图制作使用了Origin和ArcMap软件。

2 结果与分析

2.1 三江源与阿里地区气温与降水变化

通过对气象站点进行ANUSPLINE气象插值得到三江源地区与阿里地区2000—2017年气温与降水栅格数据,并对数据进行趋势分析。计算结果表明三江源地区与阿里地区年平均气温总体均呈现上升趋势,但二者年际升温幅度差异较大。阿里地区从2000年的2.59℃上升至2017年的4.12℃,平均年际增长为0.085℃/a(图3、图4),而三江源地区升温相对不稳定,三江源地区年际平均气温变化大,在2000—2013年间,气温上下波动较大,从2014年开始呈现明显上升趋势, 2000—2017年平均增长速率为0.084℃/a。阿里地区年平均气温M-K趋势分析表明,2000—2017年气温均为上升趋势,其中2007年、2012年和2017年上升趋势显著(图4)。即2000—2017年阿里地区年平均气温相对稳步上升,而三江源地区2000—2014年气温波动较大,2014—2017年气温上升显著。

图3 研究区年平均气温与累计年降水变化Fig.3 Annual average temperature and cumulative annual precipitation changes in the study area

图4 研究区年平均气温与累计年降水变化趋势M-K突变性检验Fig.4 M-K mutation test of the variation trend of annual average temperature and cumulative annual precipitation in the study areaUF:正序统计量Positive sequence statistics；UB:逆序统计量Inverse sequence statistics

在降水方面,三江源地区与阿里地区年降水量均呈现不同程度的上升趋势。三江源地区年降水量从2000年363.84 mm上升至2017年的433.61 mm；而阿里地区年降水量从2000年的163.93 mm上升至2017年的204.12 mm,两地降水均呈现不同程度波动增长。降水M-K趋势分析表明,2008—2014年三江源地区降水呈显著上升趋势(α<0.05),阿里地区降水上升趋势不显著(α>0.05)。

2.2 研究区两种高寒草地植被变化时间对比分析

从2000年到2017年,三江源和阿里地区高寒草地植被平均NPP均呈现轻微上升趋势,其中三江源地区高寒草甸植被NPP平均年增长率为0.49gC m-2a-1,阿里地区高寒草原植被NPP平均年增长率为0.0454gC m-2a-1(图5)。分时间段看,两种高寒草地植被NPP变化具有一定差异,三江源地区 2000—2005年高寒草甸植被NPP出现轻微下降趋势；而阿里地区2000—2005年高寒草原植被NPP出现轻微上升趋势。从2007年到2017年,三江源地区高寒草甸植被NPP出现明显波动上升趋势,而阿里地区高寒草原植被NPP在2010—2017年呈现增长趋势。

图5 两种高寒草地植被NPP年际变化Fig.5 Interannual changes in NPP of vegetation in two types of alpine grasslands

2.3 研究区两种高寒草地NPP空间分布及变化趋势

空间上两种高寒草地植被NPP均呈现由东南向西北逐渐递减格局(图6)；阿里地区高寒草原NPP整体处于相对较低水平,月最大值NPP不超过30gC m-2a-1,NPP相对高值区集中在普兰县、措勤县和改则县南部,而日土县及改则县西北部地区,高寒草原NPP均处于最低值。三江源地区高寒草甸NPP整体比阿里地区高寒草原高,空间上高值区分布也相对更集中。如三江源地区NPP高值集中分布东北部库泽县和河南蒙古自治县以及玉树县和囊谦县连片区92.7—139.1gC m-2a-1,在玛多县和称多县地区NPP值处于中间水平,而格尔木市及治多县西北部NPP最低0—26.2gC m-2a-1；相比阿里地区,三江源地区高寒草原NPP低值集中分布在北部及西部边缘地区。

图6 研究区NPP空间分布Fig.6 Patial distribution of NPP in the study area

研究区两种高寒草地2000—2017年NPP变化趋势差异明显(图7和表1),高寒草甸NPP增长趋势较高寒草原明显。高寒草甸NPP呈上升趋势部分占总面积57.7%,下降趋势占总面积42.3%,其中呈显著上升趋势占比13%,显著下降趋势占比6.49%,高寒草甸NPP总体呈上升趋势。而高寒草原NPP上升趋势略比下降趋势占比高,上升趋势占比50.96%,下降趋势占比49.04%,总体NPP上升部分相对较少。空间上,位于三江源地区的高寒草甸NPP增长区主要集中在西南部的格尔木市、杂多县和治多县,这些地区均为气温和降水相对较低区域。阿里地区高寒草原NPP增长区主要集中在改则县和日土县北部及札达县,其中改则县北部上升趋势最显著。通过对比发现,三江源与阿里地区高寒草地NPP增长显著区均位于降水和温度相对低值区域。

图7 研究区NPP年际变化Sen+M-K趋势检验结果Fig.7 M-K trend test results of interannual changes in NPP in the study areaSen>0,|Z|>1.65表示NPP序列呈显著上升趋势；Sen>0,|Z|<1.65表示NPP序列呈上升趋势但不显著；Sen<0,|Z|>1.65表示NPP序列呈显著下降趋势,Sen<0,|Z|<1.65, 表示NPP序列呈下降趋势但不显著；以上通过显著性检验均为90%置信区间

表1 两种高寒草地NPP变化趋势M-K检验统计

2.4 放牧与高寒草地植被变化

放牧等人类活动可能对两种高寒草地植被NPP变化产生一定影响。基于三江源和阿里地区各县统计年鉴(2000—2016)中牛羊肉产量代表各县放牧强度,利用各县域草地NPP均值代表植被生长状况,统计结果表明两种高寒草地牛羊肉产量年际变化与草地NPP变化不存在明显负相关性。在三江源区,牛羊肉产量呈明显逐年上升趋势,而高寒草甸植被NPP并没有出现逐年下降趋势,而是出现一定波动变化,二者线性回归拟合优度R2为0.0352,相关系数-0.187(P>0.05)。在阿里地区,2000—2014年牛羊肉产量呈现波动上升趋势,2014—2016年牛羊肉产量出现下降趋势,而该地区高寒草原植被NPP变化并没有出现相反的变化趋势,部分年份(2014—2016年)变化趋势与牛羊肉产量变化趋同,二者线性回归拟合优度R2为0.0256,相关系数-0.16(P>0.05)(图8)。从相关系数符号可以看出,放牧强度对两种高寒草地植被NPP均具有负作用,但影响均不明显(r<-0.2,P>0.05),不是影响两种高寒草地植被变化的主要因素。

图8 两种高寒草地牛羊肉产量与植被NPP变化Fig.8 Changes of beef and mutton yield and vegetation NPP in two alpine grasslands

2.5 两种高寒草地NPP与气温、降水相关性分析

分别对两种高寒草地NPP与气温和降水变化进行相关和偏相关分析,相关与偏相关系数大小反映温度与降水对草地植被NPP的影响,该方法被广泛用于研究气象因素对植被生长的影响[36—39], 2000—2017年两种高寒草地NPP对气温和降水的响应差异显著,并且偏相关性分析效果明显好于相关性分析(图9)。

相关性分析结果表明,高寒草甸NPP与气温呈显著正相关区集中在三江源地区西部格尔木市和杂多县,整体呈正相关区占比48.14%,而在东南部班玛县、久治县和甘德县则呈现负相关。高寒草原NPP与气温呈显著正相关区集中在阿里地区东北部即改则县北部(图9)。气温对两种草地NPP的影响存在空间差异,在高寒草甸区气温越低其与NPP正相关系数越大,而在高寒草原区气温相对较高地区其与NPP正相关系数越大。NPP与降水相关性分析表明,高寒草原区NPP与降水呈正相关性面积占比49.01%,而高寒草甸区NPP与降水呈正相关面积占比20.37%。

图9 研究区2000—2017年两种高寒草地NPP与气温、降水相关性分析Fig.9 Correlation analysis of NPP, temperature and precipitation of two typical grasslands in the study area from 2000 to 2017

偏相关性分析结果表明,在控制降水不变,高寒草甸NPP与气温呈显著正相关,面积占比达64.9%,比相关性分析48.14%提高16.76%(表3和表2)。偏相关性分析对控制降水研究温度对高寒草地NPP的影响更敏感,尤其在三江源地区东南部甘德县、达日县、久治县及其周边县,气温与高寒草甸NPP呈显著正相关。相比之下,阿里地区改则县北部、札达县和革吉县中部部分地区气温与草地NPP呈显著正相关,占比达49.91%(表2)。同样,控制气温不变对降水与草地植被NPP进行偏相关性分析,结果表明在三江源大部分区域,降水与高寒草甸植被NPP不呈显著正相关,相反在东部及中部大部分县均呈负相关性,占比达77.92%(表3)；而在西南部格尔木市及杂多县西北部降水相对较少地区呈正相关,总体呈正相关性面积占比22.89%(表3)。阿里地区降水与草地植被NPP偏相关分析表明札达县、革吉县中部呈显著正相关性,改则县西北部降水与草地NPP呈正相关性,总体呈正相关性面积占比47.76%(表3)。

表2 两种高寒草地NPP与气温相关性统计

表3 两种高寒草地NPP与降水相关性统计

3 讨论

3.1 两种高寒草地气温、降水变化差异分析

三江源与阿里地区所处纬度相当,且均处于暖湿化趋势,二者年平均气温分别以0.085℃/a和0.084℃/a的速率上升,相比张镱锂等人[8]计算青藏高原1980—2010年平均每年上升0.06℃略高；而年平均降水量分别以平均每年3.87 mm和2.23 mm速率增加,相比陆晴等[40]计算整个青藏高原1982—2013年平均每年降水增加0.96 mm高,这表明2000—2017年阿里与三江源地区暖湿化速率高于青藏高原平均水平。两地升温速率相当,但三江源高寒草甸区增湿速率明显较阿里高寒草地区快,可能原因是2000—2017年三江源高寒草甸区同时受中纬度西风带和南亚季风共同影响增强,导致降水增加相对较多[41]。

3.2 两种高寒草地放牧等人类活动与植被变化

研究中采用两种高寒草地牛羊肉产量表示放牧情况,一定程度反映当地放牧活动强弱。根据结果分析,三江源区放牧强度比阿里地区强,且三江源区东部人类活动相对密集,放牧等人类活动比阿里地区强度更大,但两地高寒草地植被NPP变化主要受气候变化影响,放牧活动影响较小。这与一些学者研究结果相近,三江源区位于气候变化主导作用区[23,42]。在三江源区同时有一系列生态保护工程,对高寒草甸植被NPP具有一定影响,主要集中在三江源国家公园区域内,彭凯峰等人[17]对该区域内人类活动和气候变化对植被变化影响进行对比研究,发现降水和温度预测的植被覆盖度与观测到的植被覆盖度变化趋势基本吻合,表明三江源国家公园生态保护区植被变化主要受气候变化影响,人类活动影响较小。两种高寒草地植被NPP变化受放牧等人类活动影响较小,受气候变化影响较大。

3.3 两种高寒草地植被NPP变化与气候变化相关性差异分析

通过相关与偏相关分析得出气温是促进高寒草甸植被生长主要因素,而降水对大部分高寒草甸植被NPP促进作用不明显,甚至在降水相对丰沛地区出现显著负相关,该研究结果与周秉荣等[22]、王军邦等[43]、朱文会等[44]和Xiong等[45]研究得到气温是影响三江源区植被NPP增加的关键气象因子结论相似。Bai等[42]研究直接得出气候变暖有助于三江源高寒植被增长,与本研究结果一致。本研究结果中降水相对丰沛区降水与高寒草甸植被NPP 呈显著负相关,可能原因是高寒草甸在满足一定降水量后降水的增加意味着日照时数的减小,因此太阳有效辐射的减少,造成气温的下降,而气温是高寒草甸植被生长的主要促进因素,降水增多一定程度上会对高寒草甸植被生长产生负作用。以上温度和降水对高寒草地植被的影响在阿里高寒草原区却没有得到相似的结果,温度并不是高寒草原植被NPP增长关键因子,相反温度和降水对高寒草原植被NPP增长的影响相当(呈正相关占比均为49%左右)。一方面可能是阿里地区温度相对更高,不像三江源大部分地区温度成为植被生长关键限制因素；另一方面阿里地区降水相对匮乏,植被生长受到降水不足的胁迫,不能与相对充足的温度协调促进植被生理生态过程。因此在三江源地区高寒草甸恢复与保护应该更加注重对耐低温植物培育与保护,在阿里高寒草原区应该对耐低温耐旱植物培育与保护。

本文研究了三江源高寒草甸与阿里高寒草原植被NPP变化及其气象驱动因子差异,为青藏高原气候变化与不同草地响应研究提供参考资料。同时本研究也存在一些不足,以气温和降水表征气候影响因素,认为太阳辐射差异较小,未考虑太阳辐射差异对植被生长影响；在分析放牧强度对高寒草地植被的影响中未能将县域尺度产肉量分解到栅格尺度与植被NPP对比研究,存在一定不确定性；对于三江源与阿里地区不同生态保护工程实施对两种高寒草地NPP产生的影响缺乏研究。接下来的研究应该更深入考虑气候变化、生态保护工程共同对两种高寒草地NPP变化的影响。

4 结论

本研究主要结论如下：

(1)三江源高寒草甸与阿里高寒草原区均处于暖湿化趋势,且三江源地区增湿速率较阿里地区快。(2)三江源高寒草甸和阿里高寒草原植被NPP均呈现东南向西北递减空间格局；三江源及阿里地区高寒草地NPP增长显著区均位于降水和温度相对低值区域,且高寒草甸植被NPP增长趋势总体较高寒草原区明显。(3)两种高寒草地植被变化受放牧影响不大,受气候变化影响较大。(4)两种高寒草地植被NPP增长的气象驱动因子差异明显,气温升高是促进高寒草甸植被NPP增长关键因子,降水增多对三江源东北部高寒草甸植被增长具有一定抑制作用；而对于高寒草原植被,气温和降水增加共同促进高寒草原植被NPP增长。

致谢：NPP数据来源于“国家青藏高原科学数据中心”(http://data.tpdc.ac.cn)