李晓婷, 郭 伟,2,*, 倪向南, 卫晓依

1 西安交通大学地球环境科学系, 西安 710049 2 西安交通大学全球环境变化研究院, 西安 710049

植物物候及其对气候变化的响应一直是全球变化研究的热点。植物物候是指植物为适应其生长环境而呈现的规律性变化[1],被视为全球变化的“诊断指纹”[2]。植物对气候变化的响应首先表现为物候期的提前或推迟,带来物种水平上植物生长期的改变,进而驱动物种分布和群落组合[3],影响陆地生态系统的结构和功能。生态系统的物质能量变化过程,如碳循环(植物生长及其生产力)、水循环(蒸发蒸腾和径流)和养分循环(分解和矿化)等也直接或间接受到植物物候变化的影响[4- 6]。因此,研究植物物候与气候变化之间的关系对认识全球变化区域响应以及指导地方农牧业可持续发展都有着重要的科学价值和实际意义。

尽管植物物候对气候要素变化的响应呈现出地域性和个体差异性,温度变化仍是影响全球植物物候的主导因素,尤其在高原山地气候区[7]。物候变化反映植物生理过程对温度变化的响应机制,在不同温度区间,种子发育速率对温度的响应不同,一般认为,低温区间发育速率呈现指数型增长,适中温度区间为线性增长,而接近高温区间的极限温度阈值时,植物发育速率迅速下降[8]。因此,气候变暖背景下,全球范围内大多数植物返青提前,如北美[9-12]、欧洲[13-14]、中国[15-19]等,植被枯黄期对气候变化响应无明显一致性趋势[20- 22]。然而,植物滞育或休眠的结束不仅与早春的热量积累有关,而且与冬季的低温有关,研究表明冬季低温有利于诱导植物自然休眠和低温需冷量的积累[23- 26]。虽然温暖的冬季可以减少许多植物的越冬死亡率,但它通常也减少了打破休眠的低温需冷量的积累,从而延迟了进入热量积累的时间使得返青推迟。近年来,研究表明植物物候响应温度升高的变化趋势逐渐减缓,甚至出现相反的变化。例如内蒙古克氏针茅草原植被返青普遍推迟[27],青藏高原高寒植被出现返青推迟、枯黄提前的现象[28-29],北京果树开花也受冬季温度升高而推迟[30],这表明冬季温度升高同样对春季物候产生了显著的影响。此外,随着对植物生长发育过程的认识,植物物候模型由线性模型逐渐发展为更精确的非线性物候过程模型,不同研究基于模型应用定量分析了冬季低温对植物返青的影响[24,31]。如Shi等对比不同物候模型,发现考虑低温需冷积累显著提高了对樱花开花期的模拟精度,冬季最低温度升高延迟樱花树开花,但同时随后极度低温日数的减少又会抵消升温带来的推迟效应[24]。因此,有必要区别讨论冬季低温与春季高温对植物物候的影响。目前,已有研究发现植物物候受临近月份的温度变化的影响最为强烈[32],其中昼夜温度对植物物候的作用也存在差异[33- 35],但基于站点监测研究高寒植物物候对冬春季最低最高温度变化的响应很少,这同时也限制了高寒地区植物物候模型的发展。

高寒草地在陆地生态系统中占据重要地位,草地植物生长发育受到温度的强烈限制。青海湖地区孕育于特殊的高寒干旱环境,属于全球变化的敏感区和生态系统典型脆弱地区。近50年来,青海湖地区温度显著升高,并且日最高最低温度变化并非同步,其中日最低温度升高程度更加剧烈[36]。在此背景下,对比研究高原植物物候如何响应最高最低温度的变化,可以进一步提高对高寒植物物候与温度变化相互关系的认识。目前,青藏高原地区植物物候研究主要集中于实测或遥感反演的植物物候变化与环境因子相关分析方面[37-39],其中环境因子大多考虑年、月均温,较少研究从日分辨率上考虑最高最低温度变化对高寒植物物候的作用,也较少在植物生理层面研究其对气候变化的响应机理,导致物候在未来预测上存在很大的不确定性。因此,本文以青藏高原东北部青海湖地区高寒草甸优势种矮嵩草为研究对象,结合青海湖海北牧业气象站物候实地观测数据和同步气象监测数据,运用偏最小二乘(PLS)回归等数理统计方法定量分析1997—2010年植物物候期(返青期和枯黄期)变化规律及其与温度的关系,以期：(1)查明温度影响高寒矮嵩草物候变化的关键时期;(2)分析关键时期内矮嵩草物候对温度变化的响应规律。

1 材料与方法

1.1 研究区概况

研究区域选择为青藏高原东部青海海北牧业气象站(海北站,(36.98°N, 100.98°E))。海北站地处高寒地区,海拔约3200 m,属于典型的高原大陆性气候,气象监测数据显示1995年至2013年,年平均气温1.4℃,冬季寒冷而漫长,霜冻情况严重,年平均降水量418.9 mm,降水集中在植物生长季(5—9月)。研究区域年蒸发量为1439.7 mm,约为年降水量的3倍。研究区域优势种为典型的矮嵩草草甸,主要由多年生草本植物组成,优势种为矮嵩草 (Kobresiahumilis),次优势种有羊茅(Festucaovina)、异针茅(Stipaaliena)和垂穗披碱草(Elymusnutans)等[40]。

1.2 资料来源及方法

本研究气象数据来自青海湖海北牧业气象站,包含1997—2010年日均温(Tmean)、日最高气温(Tmax)、日最低气温(Tmin)等气象数据。高山矮嵩草物候人工观测包括返青、花序形成、开花、果实成熟和枯黄等植物生长关键时期。物候观测按照《中国物候观测方法》进行,并按照植物物候现象出现比例确定某一物候现象的始、盛、末期。对得到的1997—2010年矮嵩草物候观测结果进行儒略历(Julian days)日序换算[41],统计矮嵩草各物候发生日期。

本研究应用偏最小二乘法(Partial Least Squares, PLS)分析植物物候与温度的相关关系。PLS回归分析在处理样本数远小于变量数而自变量存在严重多重相关性的问题时具有明显优势[42],近年来在生态和气象统计等方面已发挥重要作用[43-44]。PLS分析主要的两项输出即为变量投影重要性指标VIP(Variable Importance in Projection)和标准化模型系数。变量投影重要性指标VIP反映自变量对因变量的解释能力,其阈值通常设置为0.8,表示自变量对于因变量具有显著的解释意义[45]。而标准化模型系数指示每个变量在PLS模型中的影响的强度和方向。本文采用PLS回归方法实现高原植物矮嵩草与1997—2010年气温的回归分析,具体通过R语言(R i386 3.3.2)[46]加载pls[47]和chillR[48]扩展包编程实现。进一步相关分析及显著性检验通过SPSS 24.0计算实现。

2 结果与分析

2.1 矮嵩草物候变化特征

矮嵩草是高寒草甸中的早花植物,返青时间较早,生长季为4月至10月。1997—2010年,矮嵩草平均返青期、枯黄期分别为4月18日和10月2日(表1),平均生长季长度为166天。就14年总体来看,矮嵩草返青期2010年较1997年推迟8天,而矮嵩草枯黄期提前12天。由于矮嵩草返青期推迟,枯黄期提前,生长季长度出现缩短趋势。1997—2010年内,返青期主要集中于4月中下旬,仅在2007年于4月上旬最早发生;枯黄期主要集中于10月上旬,仅在2010年于9月中旬提前结束。总体来看,矮嵩草返青和枯黄出现的最早和最晚日期相差约1个月,返青期最早和最晚发生时枯黄期同向提前和推迟,而枯黄期最早和最晚发生时对应的返青期反向推迟和提前,造成生长季长度在14 a中差异最大值达23 d。

表1 矮嵩草物候期特征值

对矮嵩草物候期年际变化特征进一步分析发现(图1),返青期在2002年前后均呈现推迟趋势,且2002年之后推迟趋势减缓,在1997—2010年整体上返青期呈微弱推迟趋势(0.37 d/10 a);1997—2010年矮嵩草枯黄期整体呈现提前趋势(8.1 d/10 a),枯黄期也在2002年出现转折,1997—2002年矮嵩草枯黄期显著推迟(P<0.05),而2002年之后枯黄期整体呈现出提前的趋势。植物返青和枯黄的年际变化带来生长季长度的改变,1997—2002年,返青期和枯黄期均表现为推迟,其中枯黄期显著推迟,使矮嵩草生长季长度略微延长;而2002—2010年返青期推迟,枯黄期提前,造成生长季逐渐缩短;整体上,1997—2010年矮嵩草返青期推迟,枯黄期提前,生长季长度呈缩短趋势(8.4 d/10 a)。

图1 1997—2010年植物物候及生长季长度变化趋势Fig.1 Temporal variation of vegetation phenology and the length of growing season from 1997 to 2010

2.2 研究区温度变化特征

研究区域1997—2010年平均气温总体上升,倾向率为0.5℃/10a,与整个全球和青藏高原气候变暖的大背景相一致。最低气温呈现显著升高趋势,倾向率为0.7℃/10a,而最高气温在波动中呈现略微升高趋势(图2)。进一步分析发现,研究区气温同样存在阶段性变化特征,1997—2002年,平均气温、最低气温和最高气温均无明显变化趋势,而2002年之后,平均气温、最低气温和最高气温均呈现出明显上升趋势,其中平均气温显著上升(1.2℃/10a,P<0.05),最高气温和最低气温也逐步升高。由此可见,研究区域最高温度和最低温度呈现出非对称型变化,最低气温升高是该地区增温的关键。

图2 1997—2010年气温变化特征Fig.2 Temporal variations of temperature from 1997 to 2010

以2002年作为矮嵩草物候期和区域气温变化趋势节点分段,进一步分析年际温度与矮嵩草物候期的相关关系(图3),结果表明仅2002—2010年年均最高温度与矮嵩草枯黄期显著负相关,因此年际温度变化对矮嵩草返青期和枯黄期影响并不显著,本研究利用PLS回归在日分辨率尺度上研究影响矮嵩草物候期的关键时期,以及植物物候响应期间温度变化的规律。

图3 1997—2002年矮嵩草物候期与同期温度的变化趋势Fig.3 Changing trends of vegetation phenology and temperature from 1997 to 2010

2.3 矮嵩草返青期对温度变化的响应

运用偏最小二乘(PLS)回归法对植物返青期和上年度5月和下年度4月之间的365日气温回归(图4),日最高温度和日最低温度对应的返青期均方根误差(RMSE)分别为1.07 d和1.09 d,表明模型对数据是适用的。基于变量投影重要性指标(VIP)和标准化模型系数,可以确定温度变化与矮嵩草返青日期具有显著相关性的关键时期分别为1月和3—4月。

1月最高温度、最低温度与返青期回归模型系数大部分为正值且重要性指标超过0.8(变量投影重要性指标阈值),表明温度越高矮嵩草返青期越晚的关系。

3—4月最低温度与返青期回归模型系数绝大部分为负值且重要性指标超过0.8(变量投影重要性指标阈值),表明最低温度越高矮嵩草返青期越早的关系,最低温度升高表现出对返青的促进作用。然而,3—4月最高温度与返青期无显著相关关系,期间最高温度变化对矮嵩草返青的影响无连续和统一的作用。

图4 1997—2011年矮嵩草返青期与15日滑动平均气温PLS回归Fig.4 The results of Partial Least Squares (PLS) regression correlating green-up dates for Kobresia humilis in Qinghai lake basin with 15-day average daily air temperature from May to April上图黑色代表变量投影重要性指标(VIP)超过0.8; 中间图黑色代表相应的标准化模型系数; 下图黑色曲线代表温度均值,灰色和黑色区域代表每日温度标准偏差

矮嵩草返青期响应关键时期气温变化的敏感系数见表2。1月的平均最低温度与返青期呈显著正相关,1月平均最低温度每升高1℃,返青期显著推迟1.7 d;而3月至4月平均最低温度与返青期呈极显著负相关,最低温度每升高1℃,返青期显著提前3.9 d。而关键时期内最高温度对矮嵩草春季返青的影响均未通过显著性检验。由此可见,最低温度变化是影响高寒草地优势种矮嵩草返青变化的关键因子。

综上,1月和3—4月是影响高寒草地植物矮嵩草返青变化的关键时期,期间最低温度变化是影响返青变化的关键因子。回归分析发现1月和3月至4月平均最低温度变化对返青期存在相反作用,1月最低温度升高不利于矮嵩草返青,而3—4月的升温反而促进返青的发生。

2.4 矮嵩草枯黄期对温度变化的响应

运用偏最小二乘回归法对植物枯黄期和上年度11月和下年度10月之间的365日气温回归(图4),日最高温度和日最低温度对应的枯黄期均方根误差(RMSE)分别为1.55 d和2.83 d,表明模型对数据是适用的。本研究仅考虑相应物候期发生前临近月份的影响,且变量投影重要性指标(VIP)需要在一定连续影响天数(>15 d)中超过0.8,因此确定7月上中旬和8月是影响矮嵩草枯黄期的关键时期。

7月上中旬仅最低温度与枯黄期回归模型系数大部分为负值且重要性指标超过0.8,最高温度与枯黄期回归无明显相关性。8月最高温度、最低温度与枯黄期回归模型系数大部分为负值且重要性指标超过0.8(变量投影重要性指标阈值),表明温度越高矮嵩草枯黄期越早的关系。

图5 1997—2011年矮嵩草枯黄期与15日滑动平均气温PLS回归Fig.5 The results of Partial Least Squares (PLS) regression correlating coloring dates for Kobresia humilis in Qinghai lake basin with 15-day average daily air temperature from November to October上图黑色代表变量投影重要性指标(VIP)超过0.8;中间图黑色代表相应的标准化模型系数;下图黑色曲线代表温度均值,灰色和黑色区域代表每日温度标准偏差

进一步分析关键时期内矮嵩草枯黄期对温度变化的响应发现,7月上中旬的平均最低温度与枯黄期呈显著负相关,最低温度每升高1℃,枯黄期显著提前2.9 d;而期间平均最高温度对枯黄期的影响未通过显著性检验。8月平均最低温度和最高温度每升高1℃,枯黄期显著提前3.8 d(P=0.016)、5.6天(P=0.03),最低温度升高对枯黄期提前的影响更显著(表2)。由此可见,关键时期内最低温度是影响高寒草地矮嵩草枯黄变化的关键因子。

7、8月是夏末秋初时期,增温可能导致的土壤水分蒸发和植物蒸腾作用提高,导致土壤水分胁迫增强,因此植物可能无法有效利用长时段的温热条件,可利用水分减少,7、8月温度升高使植物加快进入枯黄。结合温度变化来看,本文研究区域7月上中旬和8月温度14年来呈现出上升趋势,升高速率分别为1.3℃/10a和1.1℃/10a,气温越高,植株散失水分越快,则矮嵩草进入枯黄期就越早。

表2 矮嵩草物候期对关键时期最高温和最低温变化的敏感系数

显著性：***、**分别代表显著性水平通过0.01和0.05。无星号标记表示未通过显著性检验(P>0.10)

3 讨论

本研究以典型高寒草地优势种矮嵩草为研究对象,基于实地观测物候数据,分析了1997—2010年矮嵩草物候期对日最低气温、最高气温的响应特征,研究发现近14年来青海湖地区矮嵩草返青期推迟,枯黄期提前,生长季长度缩短,这与全球气候变化下北半球高纬地区植被物候变化趋势相反[49- 50]。然而在青藏高原这一高海拔地区,已有研究发现在1998年左右,青藏高原植被返青期提前趋势减弱甚至推迟[51- 52],不同学者结合遥感数据研究植被物候变化特征也发现在青藏高原地区植被枯黄提前[22,28]。这一发现表明,温度升高并不总是有利于植物生长季长度延长,在不同时期温度对植物关键物候期的影响机制存在差异。

本研究发现3—4月份温度升高有利于矮嵩草返青,而1月份温度升高反而使矮嵩草返青推迟,冬春两季温度对植物休眠进程和种子发育阶段的相反效应在高寒草甸植物中也得到证实。具体来说,1月正处于青海湖地区隆冬时期,最低气温极低,高寒植物为避免霜冻损害进入休眠期。因此1月正处于矮嵩草休眠期低温积累时期,此时温度升高不利于矮嵩草返青需冷量的积累,使得低温积累阶段时间延长,从而推迟返青期的到来。而接下来3月至4月土壤表面解冻,植物进入热量累积阶段准备复苏。此时温度升高有利于矮嵩草积累热量打破休眠并恢复生长,温度表现出对植物返青的促进作用。因此,高寒草甸植物返青同样不仅与早春的热量积累有关,而且与冬季的低温有关。以往研究多利用遥感数据反演青藏高原地区植被物候变化趋势,发现青藏高原植被返青在1998年左右发生转折,由显著提前趋势转变为微弱推迟趋势[22,51- 52],但由于应用不同方法对不同类型遥感数据提取物候期的不确定性,目前,关于青藏高原植物初期物候指标转折之后趋势仍具有较大争议。本研究基于14年实地监测植物物候数据,发现1997年之后高寒植物返青期在隆冬温度和春季温度综合影响下仅呈现自然波动,无明显提前或推迟趋势。因此,未来植物返青期的变化趋势很大程度取决于隆冬和春季温度变化对植物生长的综合影响,而这一影响随着两季温度变化的相关性可以表现为抵消作用(冬春两季温度变化正相关)或同向作用(冬春两季温度变化负相关)。

对于生长季末期物候指标,本研究发现7月上中旬、8月是影响矮嵩草枯黄期到来的关键时期。7、8月是植物生长旺盛期,在这一时期结束后植物开始枯黄,因此期间温度变化是影响植物枯黄期变化的最相关月份。与中国温带植物秋季物候对夏秋季升温的响应相反[21,53],本研究发现青藏高原高寒草甸优势种矮嵩草随夏季(7月上中旬、8月)温度升高,枯黄期提前。因此,夏秋季升温对植物秋季物候的影响表现出地理差异性,其推迟效应在青藏高原高寒地区更为明显。7、8月是夏末秋初时期,增温可能导致土壤水分胁迫增强,而有研究表明青藏高原年平均降水量在减少,尤其是夏季降水量显著减少[54],植物可利用水分减少可能加速其进入枯黄期。而在青藏高原植被秋季物候对水分的响应较为复杂,有研究表明8月限制植物生长的主要因素为降水,降水增多促进生长季末期物候指标推迟[22],也有研究通过人工增雪实验发现矮生嵩草对水分添加的响应不明显[55]。由于本实验主要研究温度对高原植物物候的影响,降水和土壤水分对植物物候尤其是枯黄期的影响会在后续实验和观测研究中继续深入。

本研究发现最低温度是影响高寒草甸植物物候的关键因子,最低温度的变化显著影响矮嵩草的返青与枯黄。海北站位于青藏高原东北部,是典型的高寒干旱半干旱地区,此地植物生长发育处于生物极限水平的水热条件,因此高寒植物对温度变化尤其是低温的扰动可能更为敏感。极低的夜间最低温度和霜冻灾害威胁着高寒植物的生长,也对植物返青前的休眠解除、发芽和展叶等生理过程起到抑制作用[56]。在鲜有降水的冬季,高寒植物返青所需水分主要依靠土壤融雪解冻,而较高的最低温度有利于土壤解冻,从而为植物返青创造有利的水热条件[57-58]。而适应常年低温的高寒植物对夏秋季最低温度的变化同样更敏感。此外,对研究区域温度年际变化趋势分析发现,青海湖流域14年来温度年际、月份变化基本呈现上升趋势,但仅全年平均最低温度、1月、7月、9月最低温度呈现显著升高趋势,而最高温度均无显著变化,这与Liu等[36]对青藏高原近50年来最高温度、最低温度变化特征认识一致,最低温度的变化较最高温度明显,对高寒植物物候潜在影响更大。

4 结论

温度对高寒草甸植物物候的影响十分复杂,温度升高并非总是有利于植物返青、延缓枯黄和生长季延长。本研究发现,1月和3—4月的温度变化对高寒草甸优势种矮嵩草返青的作用方向不同,1月温度升高影响高寒植物休眠进程进而延迟返青,而3—4月温度升高有利于热量积累使返青提前。这说明植物春季物候对冬春两季温度升高的响应机制具有普适性,高寒地区植物返青也遵循低温积累和高温积累两段时期的温度调控。此外,影响矮嵩草枯黄的关键时期为7月上中旬和8月,期间温度升高使矮嵩草枯黄期提前。最低温度是影响高寒植物矮嵩草物候期的关键因子,其变化对高寒草甸植物物候潜在影响更为显著,这将提高对高寒植物物候影响因子的认识,优化物候过程模型的预测精度。