婁仲山

摘要?利用1961—2018年青海省共和縣氣象局觀測的逐日平均氣溫資料，采用線性趨勢、滑動平均、Mann-Kendall法等數(shù)理統(tǒng)計方法分析了春小麥生長季平均氣溫、氣溫穩(wěn)定通過0、5、10 ℃積溫的變化趨勢及突變特征。結果表明，近58年共和地區(qū)春小麥生長季平均氣溫以0.37 ℃/10 a的速率呈極顯著升溫趨勢。春小麥生長季氣溫穩(wěn)定通過0、5、10 ℃積溫也呈極顯著增溫趨勢，近58年氣溫穩(wěn)定通過0、5、10 ℃積溫分別增加418.0、449.8、568.2 ℃·d，且2000年以后熱量資源增加更顯著。春小麥生長季平均氣溫的突變時間從1991年開始，突變后平均氣溫較突變前升高了1.2 ℃。春小麥生長季氣溫穩(wěn)定通過0、5 ℃積溫的突變時間分別為1993和1994年，氣溫穩(wěn)定通過10 ℃積溫的突變時間為1988年。

關鍵詞?春小麥;生長季;熱量資源;變化特征;共和地區(qū)

中圖分類號?S162文獻標識碼?A

文章編號?0517-6611（2020）02-0228-04

doi：10.3969/j.issn.0517-6611.2020.02.066

開放科學（資源服務）標識碼（OSID）：

Variation Characteristics of Heat Resources in Spring Wheat Growing Season in Gonghe Region

LOU Zhong-shan?（Gonghe County Meteorological Bureau，Gonghe，Qinghai 813099）

Abstract?Using the daily average temperature data observed by the Meteorological Bureau of Gonghe County， Qinghai Province from 1961 to 2018， the linear trend， the moving average， the Mann-Kendall method and other mathematical statistics methods were used to analysis the variation trend and mutation characteristics of the average temperature and the accumulated temperature at 0， 5 and 10 ℃ in spring wheat growing season.The results showed that the average temperature of spring wheat growing season in Gonghe region from 1961 to 2018 showed extremely significant warming trend at a rate of 0.37 ℃/10 a.The temperature stability through 0， 5 and 10 ℃ accumulated temperature also showed extremely significant temperature increase trend in spring wheat growing season. In the past 58 years， the temperature stability through 0， 5 and 10 ℃ accumulated temperature had increased by 418.0， 449.8， 568.2 ℃·d， respectively， and after 2000， the increase in thermal resources was more significant.The sudden change of the average temperature in spring wheat growing season began in 1991， the average temperature after the mutation increased by 1.2 ℃ before the mutation.The mutation time of the temperature stability through 0 and 5 ℃ accumulated temperature was 1993 and 1994， respectively， and the mutation time of the temperature stability through 10 ℃ accumulated temperature was 1988.

Key words?Spring wheat;Growing season;Heat resources;Variation characteristics;Gonghe region

氣候變暖已經(jīng)是一個不爭的事實，IPCC第五次評估報告指出，過去130年全球平均溫度增加了0.85 ℃[1]。作物的生長發(fā)育不僅需要在一定的溫度條件下進行，而且只有當熱量累積到一定程度，才能完成其全生育期過程并獲得產(chǎn)量[2-3]。光照資源和熱量資源等農(nóng)業(yè)氣候資源直接參與農(nóng)業(yè)生產(chǎn)過程，其數(shù)量、組成及其分配狀況在一定程度上決定了一個地區(qū)的農(nóng)業(yè)格局和農(nóng)業(yè)生產(chǎn)潛力[4]。李英年等[5]研究了氣候變暖對青海農(nóng)業(yè)生產(chǎn)格局的影響。伏洋等[6]研究了氣候變化對青海省種植業(yè)的影響及適應對策。劉彩虹等[7]利用1961—2009年青海海東農(nóng)業(yè)區(qū)有關氣象、農(nóng)業(yè)觀測資料，對農(nóng)業(yè)區(qū)氣候變化特征及氣候變化對農(nóng)業(yè)的影響進行了分析，預估了未來氣候變化對海東農(nóng)業(yè)產(chǎn)生的可能影響。劉義花等[8]利用湟水河地區(qū)9個氣象站1961—2012年逐日氣溫、降水、積溫等氣象數(shù)據(jù)以及干旱、暴雨、干熱風、霜凍等災害資料，應用數(shù)理統(tǒng)計方法，研究了近52年湟水河流域氣候資源的變化趨勢。此外，不同學者對全國不同區(qū)域農(nóng)業(yè)氣候資源的時空變化特征進行了研究[9-13]。

共和縣隸屬青海省海南藏族自治州，地處青藏高原東北角，平均海拔3 200 m，全縣有可利用草場125.08萬hm2，耕地資源3.13萬hm2，農(nóng)作物總播種面積2.96萬hm2，其中糧食作物播種面積1.70萬hm2，油料作物播種面積9 126.7 hm2。

在全球氣候變暖的背景下，共和地區(qū)19世紀80年代中期以后氣溫升高趨勢顯著，尤其是21世紀以來，變暖趨勢更加明顯[14]。根據(jù)研究，共和地區(qū)活動積溫、日照時數(shù)等農(nóng)業(yè)氣候資源均發(fā)生了一定的時空變化，將會對共和地區(qū)的種植制度和農(nóng)業(yè)布局產(chǎn)生一定的影響[15-16]。因此分析共和地區(qū)氣候變化狀況下熱量資源變化特征，對充分利用熱量資源、搭配適宜的作物和品種，爭取農(nóng)業(yè)的高產(chǎn)和穩(wěn)產(chǎn)，提高熱量資源利用率均具有重要意義。

1?資料與方法

1.1?資料來源

氣象資料來源于共和縣氣象局的觀測資料，包括1961—2018年逐月溫度、日平均氣溫穩(wěn)定通過0、5、10 ℃積溫等氣象資料。多年平均值采用1981—2010年的平均值，春小麥生長季為4—9月。

1.2?研究方法

1.2.1?線性趨勢。

用xi表示樣本量為n的某一氣候變量，用ti表示xi所對應的時間，建立xi與ti之間的一元線性回歸方程[17]：

xi=a+bti（i=1，2，3，…，n）（1）

式中，a為回歸常數(shù);b為回歸系數(shù)，b和a用最小二乘法進行估計。以b×10為氣候傾向率，b>0表示直線遞增，b<0表示直線遞減，并對結果分別進行置信度為95%和99%的顯著性檢驗。

1.2.2?日滑動平均法。采用5日滑動平均法確定春小麥生長季內(nèi)≥0 ℃、≥5 ℃、≥10 ℃界線溫度的起止日期，計算穩(wěn)定通過0、5 10 ℃的積溫[18-20]。

1.2.3?Mannan-Kendall檢驗法。

Mannan-Kendall檢驗法（簡稱M-K法）是世界氣象組織推薦的用于提取序列變化趨勢的有效工具，其優(yōu)點在于不受個別異常值的干擾，能夠客觀反映時間序列趨勢[21-22]。目前已廣泛用于氣候參數(shù)和水文序列的分析中[23-24]。M-K法可以根據(jù)輸出的2個序列（UF和UB）明確突變的時段區(qū)域。若UF>0且超過臨界直線，表明時間序列呈顯著上升趨勢（P<0.05），反之，若UF<0且超過臨界直線，表明時間序列呈顯著下降趨勢（P<0.05）。如果UF和UB 2條曲線相交點在臨界直線之間，則突變開始的時間即為交點所對應的時刻，具體計算方法及參數(shù)說明見參考文獻[21-24]。

1.3?數(shù)據(jù)處理?應用DPS 7.05、Excel 2003軟件對數(shù)據(jù)進行統(tǒng)計處理、分析及繪圖。

2?結果與分析

2.1?春小麥生長季光熱資源年際變化趨勢分析

2.1.1?氣溫變化。

氣溫是影響春小麥生長的重要指標之一，春小麥生長季適宜的溫度與其產(chǎn)量密切相關。由圖1可見，1961—2018年共和地區(qū)春小麥生長季平均氣溫在波動中呈升溫趨勢，氣候傾向率為0.37 ℃/10 a，相關系數(shù)為0.807，通過0.001信度顯著水平檢驗，即近58年增溫2.2 ℃，增溫趨勢

極顯著，并呈現(xiàn)出階段性變化。春小麥生長季平均氣溫年際變化振幅較大，最高值出現(xiàn)在1998和2013年，為13.8 ℃;最低值出現(xiàn)在1976年，為10.6 ℃，生長季平均氣溫的年極差達3.2 ℃。在1961—1997年的37年間只有6年高于或等于多年值，其他31年平均氣溫均低于多年值。在1998—2018年的21年中春小麥生長季平均氣溫均高于多年值。由表1可見，共和地區(qū)春小麥生長季多年平均氣溫為12.2 ℃，20世紀60、70、80年代均低于平均值，屬于氣溫偏低期;從90年代至今氣溫高于平均值，呈現(xiàn)出氣溫快速升高趨勢。

2.1.2?積溫變化趨勢。

熱量條件與作物的生長、發(fā)育和產(chǎn)量形成密切相關，其中積溫是作物生長發(fā)育的重要指標之一。圖2顯示，1961—2018年共和地區(qū)春小麥生長季穩(wěn)定通過0、5、10 ℃的積溫均呈增加趨勢，氣候傾向率分別為72.0、77.5、98.0 ℃·d/10 a，春小麥生長季氣溫穩(wěn)定通過0、5、10 ℃的積溫與年份之間的線性相關系數(shù)分別為0.776、0.775和0.687，均通過0.001信度顯著水平檢驗，即近58年春小麥生長季氣溫穩(wěn)定通過0、5、10 ℃的積溫分別增溫418.0、449.8、568.2 ℃·d，增溫趨勢極顯著。

由表1可見，1961—2018年共和地區(qū)春小麥生長季氣溫穩(wěn)定通過0、5、10 ℃積溫平均值分別為2 252.9、2 165.3、1 706.7 ℃·d，20世紀60、70、80年代均低于平均值，屬于積溫偏低期;從20世紀90年代至今積溫高于平均值，呈現(xiàn)出積溫快速增多趨勢。春小麥生長季氣溫穩(wěn)定通過0、5、10 ℃積溫的年際變化振幅很大，最高值分別為2 543.3 ℃·d（2016年）、2 497.4 ℃·d（1998年）、2 099.8 ℃·d（2008年），最低值分別為1 947.0 ℃·d（1976年）、1 763.2 ℃·d（1982年）、1 139.0 ℃·d（1967年），春小麥生長季氣溫穩(wěn)定通過0、5、10 ℃積溫的年極差分別達596.3、734.2、960.8 ℃·d。積溫的增加將對春小麥生長季的延長生長起到顯著的促進作用。

2.2?春小麥生長季熱量資源不同時段變化趨勢分析

將1961—2018年劃分為3個時段（1961—1980、1981—2000、2001—2018年），計算每個時段內(nèi)共和地區(qū)春小麥生長季熱量資源平均值。從表2可以看出，共和地區(qū)春小麥生長季平均氣溫在不同時段呈現(xiàn)出升高趨勢，P2時段（1981—2000年）比P1時段（1961—1980年）升高了0.5 ℃，P3時段（2001—2018年）較P2時段升高了0.9 ℃。春小麥生長季氣溫穩(wěn)定通過0、5、10 ℃積溫P2時段比P1時段分別增加了90.1、86.7、178.2 ℃·d，P3時段較P2時段分別增加了197.3、197.5、195.0 ℃·d?？梢?，春小麥生長季平均氣溫和積溫在P2～P3時段的增加值高于P1～P2時段的增加值。

2.3?春小麥生長季光熱氣象要素的突變分析

2.3.1?平均氣溫的突變趨勢。從圖3可以看出，近58年來共和地區(qū)春小麥生長季平均氣溫在1969—1981、1988—2018年呈急劇升高趨勢，1961—1968年呈現(xiàn)波動中升降，1982—1987年呈現(xiàn)明顯下降趨勢， UF線于1979年超過+1.96信度水平線，在0.05信度水平下，UF線和UB線相交于1991年，氣溫升高趨勢顯著，可以確定1961—2018年共和地區(qū)春小麥生長季平均氣溫的突變時間為1991年，突變后平均氣溫較突變前升高了1.2 ℃。

2.3.2?氣溫穩(wěn)定通過0、5、10 ℃積溫的突變趨勢。從圖4可以看出，1961—2018年共和地區(qū)春小麥生長季氣溫穩(wěn)定通過0、5 ℃的積溫變化趨勢基本一致，在1969—1981、1988—2018年呈急劇升高趨勢，1961—1968年呈現(xiàn)波動中下降趨勢，1982—1987年呈現(xiàn)明顯下降趨勢，UF線和UB線分別相交于1993和1994年，0、5 ℃的積溫增加趨勢顯著，可以確定1961—2018年共和地區(qū)春小麥生長季氣溫穩(wěn)定通過0、5 ℃積溫的突變時間分別為1993和1994年。

近58年共和地區(qū)春小麥生長季氣溫穩(wěn)定通過10 ℃積溫變化趨勢與0、5 ℃積溫不盡相同，在1961—1965年呈急劇下降趨勢，1966—1979年波動中持續(xù)升高，1981—1987年呈持續(xù)下降趨勢，從1988年開始積溫顯著升高。UF和UB曲線相交于1988年，1992年UF曲線突破信度0.05水平線，可以確定1961—2018年共和地區(qū)春小麥生長季氣溫穩(wěn)定通過10 ℃積溫的突變時間為1988年。

3?結論與討論

（1）1961—2018年共和地區(qū)春小麥生長季平均氣溫以0.37 ℃/10 a的速率呈極顯著升溫趨勢，近58年增溫2.2 ℃。春小麥生長季氣溫穩(wěn)定通過0、5、10 ℃積溫也呈極顯著增溫趨勢，近58年氣溫穩(wěn)定通過0、5、10 ℃積溫分別增加418.0、449.8、568.2 ℃·d，且2000年以后熱量資源增加更顯著。春小麥生長季平均氣溫和積溫在P2（1981—2000年）～P3（2001—2018年）時段的增加值高于P1（1961—2000年）～P2時段的增加值。劉義花等[8]研究指出，氣候變暖使湟水河流域0、10 ℃的積溫總體呈增加趨勢，熱量資源更加豐富，與該研究結果一致，張曉云等[25]也有類似的研究結果。劉彩虹等[7]研究表明，溫度上升加快了春小麥的生理發(fā)育速度，營養(yǎng)生長階段提前，全生育期延長，發(fā)育期縮短。積溫的增加，使得春播作物播種期提早，低溫冷害的不良影響有所緩解。

（2）春小麥生長季平均氣溫的突變時間為1991年，突變后平均氣溫較突變前升高了1.2 ℃。春小麥生長季氣溫穩(wěn)定通過0、5和10 ℃積溫的突變時間分別為1993、1994、1988年。

參考文獻

[1] IPCC.Climate change 2013：The physical science base[M].Cambridge，UK：Cambridge University Press，2013：1-300.

[2] 李樹巖，方文松，馬志紅.河南省夏玉米生長季農(nóng)業(yè)氣候資源變化分析[J].河南農(nóng)業(yè)科學，2012，41（7）：21-26.

[3] 王位泰，張?zhí)旆?，蒲金涌，?黃土高原中部冬小麥生長對氣候變暖和春季晚霜凍變化的響應[J].中國農(nóng)業(yè)氣象，2011，32（1）：6-11.

[4] 葉彩華，欒慶祖，胡寶昆，等.北京農(nóng)業(yè)氣候資源變化特征及其對不同種植模式玉米各生育期的影響[J].自然資源學報，2010，25（8）：1350-1364.

[5] 李英年，王啟基.氣候變暖對青海農(nóng)業(yè)生產(chǎn)格局的影響[J].西北農(nóng)業(yè)學報，1999，8（2）：102-107.

[6] 伏洋，張國勝，顏亮東.氣候變化對青海省種植業(yè)的影響及適應對策[J].中國農(nóng)業(yè)氣象，2004，35（3）：11-14.

[7] 劉彩虹，祁貴明，戴升.氣候變化對青海海東農(nóng)業(yè)區(qū)的影響評估與適應對策[J].安徽農(nóng)業(yè)科學，2011，39（11）：6608-6610，6613.

[8] 劉義花，周強，魯延榮，等.湟水河流域近50年來農(nóng)業(yè)氣候資源變化[J].中國農(nóng)學通報，2016，32（12）：163-170.

[9] 孫楊，張雪芹，鄭度.氣候變暖對西北干旱區(qū)農(nóng)業(yè)氣候資源的影響[J].自然資源學報，2010，25（7）：1153-1162.

[10] 龔強，汪宏宇，張運福，等.氣候變化背景下遼寧省氣候資源變化特征分析[J].資源科學，2010，32（4）：671-678.

[11] 董旭光，李勝利，石振彬，等.近50年山東省農(nóng)業(yè)氣候資源變化特征[J].應用生態(tài)學報，2015，26（1）：269-277.

[12] 陳超，龐艷梅，潘學標，等.氣候變化背景下四川省氣候資源變化趨勢分析[J].資源科學，2011，33（7）：1310-1316.

[13] 郭良才，岳虎，王強，等.河西走廊干旱區(qū)農(nóng)業(yè)氣候資源變化特征[J].干旱地區(qū)農(nóng)業(yè)研究，2008，26（3）：14-22.

[14] 劉運華，毛萬珍，李樂，等.近53年共和地區(qū)作物生長季氣候變化特征分析[J].農(nóng)業(yè)災害研究，2016，6（5）：31-33.

[15] 趙恒和，郭連云，趙年武，等.共和盆地界限溫度初終日和積溫對氣候變化的響應[J].水土保持研究，2012，19（4）：207-211.

[16] 屈柳燕，王靖，劉澤華，等.青海省共和縣小麥和油菜物候?qū)夂蜃兓捻憫猍J].青海草業(yè)，2017，26（3）：31-36.

[17] 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預測技術[M].北京：氣象出版社，2007：37-59.

[18]馮定原.農(nóng)業(yè)氣象預報和情報方法[M].北京：氣象出版社，1988：101-120.

[19] 高紹風，陳萬隆，朱超群，等.應用氣候?qū)W[M].北京：氣象出版社，2004：40-46.

[20] 馮秀藻，陶炳炎.農(nóng)業(yè)氣象學原理[M].北京：氣象出版社，1991：72-127.

[21] 符淙斌，王強.氣候突變的定義和檢測方法[J].大氣科學，1992，16（4）：482-493.

[22] 吳洪寶，吳蕾.氣候變率診斷和預測方法[M].北京：氣象出版社，2005：208-244.

[23] 張涵丹，衛(wèi)偉，薛萐.基于R/S分析和Mann-Kendall檢驗的定西市氣溫降水變化特征[J].水土保持研究，2015，22（6）：183-189.

[24] 康淑媛，張勃，柳景峰，等.基于Mann-Kendall法的張掖市降水量時空分布規(guī)律分析[J].資源科學，2009，31（3）：501-508.

[25] 張曉云，孔祥萍，朱元福.1961—2017年貴南地區(qū)氣溫穩(wěn)定通過10 ℃界限溫度變化特征[J].中國農(nóng)學通報，2018，34（28）：101-105.