徐韻佳,葛全勝,戴君虎,王煥炯,*

1 中國科學(xué)院地理科學(xué)與資源研究所, 陸地表層格局與模擬重點實驗室, 北京 100101 2 中國科學(xué)院大學(xué), 北京 100049

近幾十年來,北半球溫帶地區(qū)大部分植物春季物候期明顯提前,而秋季物候期的變化空間差異較大[1-3]。與其他氣候因素(日照、降水量等)相比,溫度升高是造成木本植物春季物候期提前的最主要因素[4-5]。因此,春季物候期對溫度變化的響應(yīng)程度受到廣泛關(guān)注。溫度每升高1℃,植物物候期變化的天數(shù)定義為溫度敏感度[6-7]。研究表明,春季物候溫度敏感度高的物種,在群落中的覆蓋度和生物量占據(jù)優(yōu)勢[8]。因此,研究物候的溫度敏感度變化有助于評估不同植物對氣候變化的響應(yīng)與適應(yīng)能力,從而掌握未來生態(tài)系統(tǒng)中的種群動態(tài)[8-9]。

很多學(xué)者研究了植物春季物候溫度敏感度在不同時期的變化特征。例如,在瑞士地區(qū)的研究表明,植物春季物候的溫度敏感度在過去300年總體呈升高態(tài)勢,而在1890—1950和1970—2007年兩個典型暖期出現(xiàn)了降低趨勢[10]。在中國,Dai等[11]分析了牡丹江地區(qū)40個木本植物1978—2014年春季展葉始期溫度敏感度的變化特征,發(fā)現(xiàn)75%的植物溫度敏感度呈上升趨勢,25%呈降低趨勢。在歐洲地區(qū),Fu等[6]分析了1980—1994年和1999—2013年兩個時段7種廣泛分布的木本植物展葉始期對溫度的敏感度變化,發(fā)現(xiàn)后時段的平均溫度敏感度相較于前時段下降了40%。但當比較30年變化時,兩時段(1951—1980年和1984—2013年)溫度敏感度不存在顯著差異[12]。

造成上述溫度敏感度時間變化的原因仍有爭議。一些學(xué)者認為植物春季物候期溫度敏感度變化與植物在冬季受到的冷激量有關(guān)[6,13]。當冬季溫度升高時,植物所經(jīng)歷的冷激量降低,從而使萌芽所需的積溫量增加,降低春季物候期的溫度敏感度[14-16]。也有學(xué)者認為光周期與植物春季物候期溫度敏感度有關(guān)。一些樹種對日照長度較為敏感,當展葉始期提前到一定程度,達到了植物開始發(fā)育晝長的臨界值,展葉始期將不再提前[17]。另外,有研究表明春季氣溫變率與植物展葉始期的敏感度有關(guān)[18-20]。春季氣溫變率較大的地區(qū)更易發(fā)生霜凍害事件[21],植物通過降低溫度敏感度這一適應(yīng)策略來避免遭受霜凍災(zāi)害[22-23]。

國內(nèi)對溫度敏感度的研究主要集中在個別站點的種間差異和時間變化特征[11,24-25]或少數(shù)物種的空間差異方面[26-27],而對于中國近50年溫度敏感度的變化情況尚無系統(tǒng)報導(dǎo)。為增強對未來物候變化的預(yù)測能力,亟需對過去50年來中國地區(qū)植物春季物候期的溫度敏感度變化特征及影響因素進行深入研究。本文選擇1963—2014年中國中東部地區(qū)10個站點163種植物的展葉始期資料,利用滑動分析法計算了每15a各植物的展葉始期溫度敏感度。在此基礎(chǔ)上,分析了溫度敏感度的變化趨勢及空間格局,并探討了導(dǎo)致展葉始期溫度敏感度變化的可能原因。

1 數(shù)據(jù)與方法

1.1 物候與氣象數(shù)據(jù)

本文所用物候觀測數(shù)據(jù)來自“中國物候觀測網(wǎng)”(China Phenological Observation Network, http://www.cpon.ac.cn/)。根據(jù)觀測序列較為連續(xù)、物種較豐富的原則,研究選取了中東部地區(qū)哈爾濱、北京、長沙等10個站點(表1)。物候數(shù)據(jù)時間跨度為1963—2014年,各站點觀測年數(shù)在21—43a之間。為研究溫度敏感度的時間變化,篩選出各站點觀測年數(shù)超過20年的木本植物,共計163種、313條序列。這些植物隸屬48個科,103個屬,包括66種灌木,97種喬木。由于觀測中斷,1969—1972和1997—2002年兩時段的物候數(shù)據(jù)缺失,所選序列主要分布于1963—1968、1973—1996和2003—2014年(圖1)。各觀測站點對應(yīng)的氣象資料從中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.gov.cn/)獲取,包括1963—2014年的逐日平均氣溫數(shù)據(jù)。對由于遷站等原因造成氣溫數(shù)據(jù)不均一的站點(如貴陽、西安等),使用中國氣象數(shù)據(jù)網(wǎng)提供的“中國近50年均一化歷史氣溫數(shù)據(jù)集”。

表1 研究站點地理位置及物候資料概況

圖1 1963—2014年中國10個站點展葉始期逐年觀測記錄數(shù)量Fig.1 The number of the observation records for leaf unfolding date (LUD) at 10 sites of China from 1963 to 2014

1.2 研究方法

1.2.1溫度敏感度計算方法

為估計溫度敏感度,首先需要確定溫度對植物展葉始期影響最為顯著的時段[28-29],即在多年平均展葉始期前一段時間內(nèi),以15d為單位滑動計算每一時段(即展葉始期前15d,展葉始期前30d,直至前120d)的平均溫度與展葉始期的相關(guān)系數(shù),相關(guān)系數(shù)絕對值最大的時段作為響應(yīng)時段。展葉始期與響應(yīng)時段內(nèi)平均氣溫的回歸方程斜率,可用來表示展葉始期對溫度變化響應(yīng)的敏感度。

為評估展葉始期溫度敏感度的變化趨勢,本研究首先選擇連續(xù)的15a時間窗口(例如1963—1977,1974—1988,…,2000—2014),然后對時間窗口內(nèi)展葉始期序列與響應(yīng)時段平均氣溫進行線性回歸分析,回歸系數(shù)用來表征每15a植物的展葉始期溫度敏感度[6]。當某一時間窗口內(nèi)觀測年數(shù)不足8a時,不予計算。其次,對溫度敏感度與時間窗口的終止年份進行線性回歸分析,回歸系數(shù)即為該物種展葉始期溫度敏感度的變化速率。由于溫度敏感度在大多數(shù)情況下為負值,因此溫度敏感度在數(shù)值上越小意味著絕對值越大,敏感度越高。

1.2.2溫度敏感度變化的歸因分析

根據(jù)引言部分的文獻綜述,植物展葉始期溫度敏感度可能與3個因子有關(guān)：植物在冬季受到的冷激量,光周期和春季氣溫變率。本研究分別用以下變量來表征這3個因子：

(1)冷激天數(shù)(CD),用前一年11月1日至平均展葉始期間的平均溫度介于0到5℃之間的日數(shù)總和來表示[30]。根據(jù)本研究的假設(shè),植物在冬季接受的冷激天數(shù)降低會導(dǎo)致植物展葉始期溫度敏感度降低。

(2)平均展葉始期(LUD)。在固定站點,晝夜長短只與日期有關(guān)。當物候期提前到一定日期,晝長縮短,展葉可能受到限制而不再提前。因此,可以用植物平均展葉始期這一變量來表征植物物候發(fā)生時的光周期情況。根據(jù)本研究的假設(shè),植物展葉始期提前會導(dǎo)致溫度敏感度降低。

(3)春季氣溫標準差(SD),通過計算響應(yīng)時段內(nèi)氣溫標準差的多年平均值來計算[18,20]。根據(jù)本研究的假設(shè),春季氣溫標準差增大,植物溫度敏感度降低。

對每條展葉始期序列采用隨機抽樣法進行歸因分析。具體做法是：對每條序列隨機選擇15a形成一組春季物候期和溫度序列,計算該隨機序列的溫度敏感度以及3個變量(CD,LUD和SD)。將隨機抽樣實驗重復(fù)1000次,對隨機形成的各時段溫度敏感度和各變量進行Pearson偏相關(guān)分析。通過分析偏相關(guān)系數(shù)的顯著度和正負號,來判斷各因子對溫度敏感度的影響是否具有顯著性和一致性。

2 結(jié)果分析

2.1 展葉始期溫度敏感度變化的總體趨勢

圖2 中國地區(qū)313條展葉始期序列溫度敏感度變化趨勢的頻率分布(1963—2014年)Fig.2 Frequency distribution of the trends in temperature sensitivity (ST) for 313 time series of leaf unfolding date (LUD) in China (1963—2014) 虛線表示變化趨勢的中值(-0.15d ℃-1 (10a)-1)

從植物展葉始期溫度敏感度變化的總體趨勢來看,溫度敏感度升高的序列明顯多于降低的序列。313條展葉始期溫度敏感度變化序列中,60.1%的序列呈升高趨勢,其中顯著的占40.0%(P<0.05);39.9%的溫度敏感度序列降低,其中顯著降低的占28.4%(P<0.05)(圖2)。就變化幅度而言,展葉始期溫度敏感度的變化趨勢在-7.03d ℃-1(10a)-1(武漢的馬尾松Pinusmassoniana)到4.96d ℃-1(10a)-1(長沙的楓香樹Liquidambarformosana)之間,中值為-0.15d ℃-1(10a)-1。但81.2%的序列溫度敏感度變化趨勢絕對值小于1d ℃-1(10a)-1(圖2)。

通過對所有展葉始期序列在各個時間窗口內(nèi)的統(tǒng)計分析發(fā)現(xiàn),1977—2014年展葉始期溫度敏感度的平均值和中值呈波動變化。2004年以前主要呈升高趨勢,溫度敏感度中值從-2.43d/℃增加至-4.08d/℃。2004年后主要呈降低趨勢,中值重新降低到約-3.46d/℃(圖3)。從線性趨勢來看,所有展葉始期序列溫度敏感度的中值顯著降低,但變化速率較小,僅-0.3d ℃-1(10a)-1(R2=0.64,P<0.05)。

圖3 1977—2014年中國地區(qū)313條展葉始期序列15年滑動溫度敏感度變化Fig.3 Temperature sensitivity of leaf unfolding date for each 15-year moving window of 313 time series in China (1977—2014) 箱底部和頂部代表四分位數(shù);箱內(nèi)橫線代表中值;箱內(nèi)圓點代表平均值;上下兩個須的終端分別代表10和90分位數(shù);紅線為中值的線性擬合線：y=-0.03x+53.8, R2=0.64, P<0.05

2.2 展葉始期溫度敏感度變化的站點差異

圖4 中國10個站點植物展葉始期溫度敏感度變化的空間分布Fig.4 The spatial patterns of trends in temperature sensitivity of leaf unfolding date at 10 sites in China

從植物展葉始期溫度敏感度變化的空間格局來看,分布在溫帶地區(qū)的6個站點溫度敏感度升高的物種數(shù)量均超過60%(圖4)。其中,北京植物展葉始期溫度敏感度升高最為普遍,升高的物種比例達到89.3%。其中顯著升高的物種占75.0%(P<0.05),平均變化速率為-0.43d ℃-1(10a)-1。而分布在亞熱帶溫度帶的4個站點,除合肥植物展葉始期溫度敏感度升高(-0.51d ℃-1(10a)-1)外,其他站點(武漢、長沙和貴陽)均有超過65%的物種展葉始期溫度敏感度降低,且顯著降低的物種比例在45.7%到68.4%之間。

2.3 展葉始期溫度敏感度變化的影響因子

對溫帶地區(qū)站點而言,有66%的展葉始期序列溫度敏感度與冷激天數(shù)呈負相關(guān)關(guān)系,其中偏相關(guān)系數(shù)顯著的序列占58.3%(P<0.05);只有10%的展葉始期序列溫度敏感度與冷激天數(shù)呈顯著正相關(guān)關(guān)系(P<0.05)。這表明,對溫帶地區(qū)大部分植物而言,當冷激天數(shù)減少時,展葉始期溫度敏感度降低,符合預(yù)期假設(shè)。對亞熱帶地區(qū)的站點,則有超過一半(58.9%)的序列溫度敏感度與冷激天數(shù)呈正相關(guān)關(guān)系,顯著的序列達46.8%(P<0.05);只有38.7%序列溫度敏感度與冷激天數(shù)呈顯著負相關(guān)關(guān)系(P<0.05)。這表明,對亞熱帶地區(qū)大部分植物而言,隨著冷激天數(shù)減少,展葉始期溫度敏感度反而升高,不符合預(yù)期假設(shè)(圖5)。

圖5 中國地區(qū)全部站點以及分區(qū)的植物展葉始期溫度敏感度與3個因子的偏相關(guān)系數(shù)頻率分布Fig.5 Frequency distributions of the partial correlation coefficients between temperature sensitivity of leaf unfolding date and three factors

在溫帶和亞熱帶地區(qū),分別有63.7%和70.2%的序列展葉始期溫度敏感度與平均展葉始期呈正相關(guān)關(guān)系,其中顯著的序列分別有57.8%和66.1%(P<0.05)。兩個地區(qū)分別只有25.1%和24.2%的序列溫度敏感度與平均展葉始期呈顯著負相關(guān)(P<0.05)。因此,無論是在溫帶地區(qū)還是亞熱帶地區(qū),對大部分植物而言,展葉始期越早,溫度敏感度越高,不符合本研究的預(yù)期假設(shè)(圖5)。

與氣溫標準差的偏相關(guān)分析表明,在溫帶和亞熱帶地區(qū),分別有65.8%和63.7%的序列展葉始期溫度敏感度與春季氣溫標準差呈正相關(guān)關(guān)系,其中偏相關(guān)系數(shù)顯著的序列數(shù)分別有52.4%和57.2%(P<0.05)。兩個地區(qū)分別只有24.6%和25%的序列展葉始期溫度敏感度與春季氣溫標準差呈顯著負相關(guān)關(guān)系(P<0.05)。該結(jié)果表明,當春季氣溫變率增大時,大部分植物展葉始期溫度敏感度降低,符合預(yù)期假設(shè),并且在兩個區(qū)域的結(jié)果一致(圖5)。

3 討論

本研究系統(tǒng)探討了1963—2014年中國中東部地區(qū)10個站點313條展葉始期序列的溫度敏感度變化趨勢及可能原因。其中,60.1%的序列溫度敏感度呈升高趨勢,39.9%的溫度敏感度序列呈降低趨勢。該結(jié)果與歐洲地區(qū)的研究結(jié)果不同。歐洲7種廣布種展葉始期平均溫度敏感度1999—2013年相較于1980—1994年降低了40%[6]。這表明植物春季物候期的溫度敏感度的變化存在顯著的區(qū)域差異。

研究結(jié)果表明,展葉始期溫度敏感度與冷激天數(shù)的相關(guān)關(guān)系存在區(qū)域差異。在溫帶地區(qū),66%的植物展葉始期序列的溫度敏感度與冷激天數(shù)呈負相關(guān)關(guān)系,表明在溫帶地區(qū),冷激量的變化能夠解釋大部分植物溫度敏感度變化。然而,在亞熱帶地區(qū),大部分植物可能不需要或只需要很小的冷激量即可解除芽在秋冬季的休眠[31],冬季低溫不足不會使展葉始期溫度敏感度降低,表現(xiàn)為58.9%的序列展葉始期溫度敏感度與冷激天數(shù)呈正相關(guān)關(guān)系。對于冷激條件變化的響應(yīng)差異可能是導(dǎo)致溫帶和亞熱帶植物展葉始期溫度敏感度變化趨勢不同的原因。當前,冷激條件對植物葉芽發(fā)育的影響仍缺乏深刻認識,需要進一步研究。

根據(jù)光周期限制的假設(shè),在氣溫較高的年份,展葉始期對光周期敏感的植物由于晝長不足的限制而停止提前。因此,展葉始期提前,溫度敏感度應(yīng)降低。但僅有36.4%的溫帶序列和29.8%的亞熱帶序列支持這一結(jié)論,對于其余63.7%和70.2%的序列而言,隨著展葉始期提前,溫度敏感度反而升高。因此,光周期限制不能解釋大部分植物春季展葉始期溫度敏感度的變化。可能的原因是所選植物的春季物候期對光周期變化并不敏感。另一個可能原因是在當前氣候狀態(tài)下,大部分植物春季展葉始期尚未受到光周期限制。在未來氣候持續(xù)升溫的情景下,對光周期敏感的植物可能因到達晝長的閾值而不再提前[32]。

本研究結(jié)果還表明無論在溫帶還是亞熱帶地區(qū),超過一半的展葉始期序列溫度敏感度與春季氣溫標準差呈顯著正相關(guān)關(guān)系,說明春季氣溫變率對展葉始期溫度敏感度有顯著影響,且區(qū)域差異性較小?？赡艿臋C制是春季溫度波動較大的區(qū)域更易發(fā)生霜凍事件,植物通過降低對氣候變化響應(yīng)的敏感度以減小遭受霜凍災(zāi)害的風(fēng)險[18]。該結(jié)果也與溫度敏感度空間分布的影響因素一致。Wang等[18]研究發(fā)現(xiàn)歐洲不同站點間植物春季物候期溫度敏感度與春季氣溫標準差存在顯著正相關(guān)關(guān)系,在春季氣溫變率較大的地區(qū)植物春季物候敏感度較低。Zhang等[20]通過分析中國59個站點植物花期溫度敏感度分布規(guī)律也得到類似結(jié)論。

4 結(jié)論

本研究定量分析了近50年中國中東部地區(qū)10個站點163種植物展葉始期溫度敏感度的變化趨勢及空間格局,并探討了導(dǎo)致展葉始期溫度敏感度變化的可能原因,主要結(jié)論為：

(1)在全部313條展葉始期時間序列中,60.1%的序列溫度敏感度呈升高趨勢,其中顯著升高的占40.0%(P<0.05)。39.9%的序列溫度敏感度降低,其中顯著降低的占28.4%。所有植物溫度敏感度的中值在1963—2014年顯著升高。

(2)在空間分布上,溫帶地區(qū)的6個站點展葉始期溫度敏感度平均呈升高趨勢。其中,北京植物展葉始期溫度敏感度升高最為普遍,顯著升高的物種比例達到75.0%(P<0.05)。而亞熱帶站點(除合肥外)的展葉始期溫度敏感度主要呈降低趨勢。其中,長沙植物展葉始期溫度敏感度顯著降低的物種比例最高,達68.4%。

(3)冬季冷激量和春季氣溫變率是影響植物展葉始期溫度敏感度隨時間變化的主要因素。冬季冷激量降低將導(dǎo)致植物展葉始期溫度敏感度降低,而春季氣溫變率降低將導(dǎo)致植物展葉始期溫度敏感度升高。