劉宏舉 姚樹然 董航宇

摘要：選取河北省24個農(nóng)業(yè)氣象觀測站1960—2016年逐日氣象資料，以日最高氣溫、14：00相對濕度和14：00風(fēng)速等3個氣象因素作為分析依據(jù)，得到全省冬麥區(qū)近57年來干熱風(fēng)發(fā)生日數(shù)，通過線性趨勢、Mann-Kendall檢驗(yàn)法和滑動t檢驗(yàn)法等方法，研究河北省冬小麥干熱風(fēng)時(shí)空分布規(guī)律。結(jié)果表明，近57年來河北省干熱風(fēng)發(fā)生整體呈下降趨勢，其中重度干熱風(fēng)下降趨勢顯著;輕度干熱風(fēng)和干熱風(fēng)總次數(shù)在近57年來并沒有發(fā)生突變，而重度干熱風(fēng)的發(fā)生在20世紀(jì)80年代中期出現(xiàn)了突變性下降;干熱風(fēng)的空間分布呈東北部低、南部高的趨勢，河北省邯鄲市、邢臺市的太行山山前地帶以及黑龍港流域的中西部為干熱風(fēng)的重災(zāi)區(qū)。

關(guān)鍵詞：干熱風(fēng);河北省;冬小麥;M-K檢驗(yàn);時(shí)空分布;建議

中圖分類號： S429 文獻(xiàn)標(biāo)志碼： A 文章編號：1002-1302（2019）01-0242-04

干熱風(fēng)是冬小麥生育后期出現(xiàn)的一種高溫低濕并伴有一定風(fēng)力的綜合性農(nóng)業(yè)氣象災(zāi)害。干熱風(fēng)危害的實(shí)質(zhì)是高溫、低濕引起農(nóng)作物生理干旱，風(fēng)加重了危害的程度。干熱風(fēng)是河北省冬小麥生產(chǎn)的主要?dú)庀鬄?zāi)害之一，主要出現(xiàn)在冬小麥的灌漿期，抑制小麥旗葉的光合蒸騰[1-2]，導(dǎo)致冬小麥灌漿速率不同程度下降[3]，千粒質(zhì)量下降[4]，一般年份可造成小麥減產(chǎn)1%～2%，偏重年份可減產(chǎn)3%以上[5]。干熱風(fēng)主要有高溫低濕型、雨后熱枯型、旱風(fēng)型等3種類型[6]，其中，高溫低濕型是河北麥區(qū)干熱風(fēng)的主要類型，在小麥開花灌漿過程均可發(fā)生，其特點(diǎn)是高溫低濕，干熱風(fēng)發(fā)生時(shí)溫度升高，空氣濕度劇降，風(fēng)力在3～4 m/s以上;這類干熱風(fēng)發(fā)生的區(qū)域廣，使小麥芒尖干枯炸芒，穎殼呈灰白色或青灰色，葉片卷曲凋萎，從而造成小麥大面積干枯逼熟死亡，對小麥產(chǎn)量威脅很大。

近年來，關(guān)于干熱風(fēng)的研究較多，也取得了不少成果。喇永昌等利用1981—2014年寧夏灌區(qū)16個氣象站的氣象資料，分析了當(dāng)?shù)卮盒←湼邷氐蜐裥透蔁犸L(fēng)日數(shù)和天氣過程次數(shù)的時(shí)空分布特征[7];屈振江等分析了陜西省冬小麥主產(chǎn)區(qū)干熱風(fēng)的年際變化規(guī)律和空間分布特點(diǎn)[8];成林等對河南省干熱風(fēng)發(fā)生時(shí)氣象要素的變化、發(fā)生范圍、天數(shù)及頻率的變化進(jìn)行了詳細(xì)分析[9-10];楊霏云等基于冬小麥干熱風(fēng)發(fā)生指標(biāo)，構(gòu)建了冬小麥干熱風(fēng)的發(fā)生強(qiáng)度指標(biāo)，反演了華北平原冬小麥主產(chǎn)區(qū)1961年以來干熱風(fēng)的發(fā)生規(guī)律[11];鄔定榮等利用華北平原48個站點(diǎn)1961—2008年的逐日氣象資料，分析了近50年來各站點(diǎn)的干熱風(fēng)日數(shù)及區(qū)域分異規(guī)律[12];趙俊芳等研究了過去黃淮海地區(qū)冬小麥干熱風(fēng)發(fā)生的時(shí)空演變規(guī)律，并基于干熱風(fēng)危害指數(shù)，對黃淮海地區(qū)冬小麥干熱風(fēng)災(zāi)損進(jìn)行了評估[13-15]。但在現(xiàn)有研究中，針對河北省內(nèi)干熱風(fēng)分布的研究相對較少，且歷史較久遠(yuǎn)，所采用的干熱風(fēng)指標(biāo)與后來統(tǒng)一的中華人民共和國氣象行業(yè)標(biāo)準(zhǔn)（QX/T 82—2007《小麥干熱風(fēng)災(zāi)害等級》）也有一定差異[16]，不能完全反映出河北省干熱風(fēng)發(fā)生的實(shí)際情況。

河北省作為全國的小麥生產(chǎn)大省之一，小麥生產(chǎn)的豐歉，不僅關(guān)系到農(nóng)村社會經(jīng)濟(jì)發(fā)展和農(nóng)民的切身利益，還關(guān)系到全國糧食市場的供求平衡和糧食價(jià)格的穩(wěn)定;而河北省邯鄲市、邢臺市的太行山山前地帶以及黑龍港流域的中西部為干熱風(fēng)的重災(zāi)區(qū)。本研究以河北省24個農(nóng)業(yè)氣象觀測站1960—2016年逐日日最高氣溫、14：00相對濕度和14：00風(fēng)速等3個氣象因素作為分析依據(jù)，通過線性趨勢、Mann-Kendall檢驗(yàn)法和滑動t檢驗(yàn)法等方法，研究河北省冬小麥干熱風(fēng)時(shí)空分布規(guī)律，以提高干熱風(fēng)預(yù)警的準(zhǔn)確性和可靠性，降低干熱風(fēng)災(zāi)害造成的經(jīng)濟(jì)損失，也為后期進(jìn)一步建立冬小麥干熱風(fēng)監(jiān)測預(yù)警模型提供一定的理論和應(yīng)用基礎(chǔ)。

1 材料與方法

1.1 材料來源

河北省大部地區(qū)冬小麥在5月中下旬開始進(jìn)入灌漿期，因而選取每年5月中旬至6月上旬（5月10日—6月10日）的逐日氣象資料，系統(tǒng)分析干熱風(fēng)對冬小麥灌漿至成熟期的影響[4]。選取定州、藁城、南宮、黃驊、涿州等24個農(nóng)業(yè)氣象觀測站1960—2016年逐日氣象資料，以日最高氣溫、14：00相對濕度、14：00風(fēng)速等3個氣象因素作為分析依據(jù)，具體氣象站點(diǎn)空間分布見圖1。干熱風(fēng)災(zāi)害氣象指標(biāo)選用中華人民共和國氣象行業(yè)標(biāo)準(zhǔn)（QX/T 82—2007《小麥干熱風(fēng)災(zāi)害等級》）規(guī)定的相應(yīng)等級，具體見表1。

采用非參數(shù)Mann-Kendall突變檢驗(yàn)法和滑動t檢驗(yàn)法檢測干熱風(fēng)發(fā)生的突變性。Mann-Kendall檢驗(yàn)中，UF、UB分別表示正向、反向統(tǒng)計(jì)量，并給出顯著性水平α=0.05的臨界線（±1.96）。如果UF>0，則表示序列呈上升趨勢，反之則表示序列呈下降趨勢;當(dāng)超過±1.96時(shí)，表明上升或下降趨勢達(dá)到顯著水平。如果UF、UB等2條曲線出現(xiàn)交點(diǎn)，且交點(diǎn)在臨界線之間，那么交點(diǎn)對映的時(shí)刻便是突變開始的時(shí)間，當(dāng)出現(xiàn)多個突變點(diǎn)時(shí)，使滑動t檢驗(yàn)對突變點(diǎn)進(jìn)一步確定[17]。

在ArcGIS軟件下通過反距離權(quán)重插值法（inverse distance weighted，簡稱IDW）定量分析河北省冬麥區(qū)輕度、重度干熱風(fēng)災(zāi)害的時(shí)空變化特征。

2 結(jié)果與分析

2.1 時(shí)間變化特征

2.1.1 年際變化特征 由圖2-a可知，河北省輕度干熱風(fēng)發(fā)生最大值出現(xiàn)在1967年，發(fā)生次數(shù)為120站次;最小值出現(xiàn)在1963年，發(fā)生12站次;57年來輕度干熱風(fēng)發(fā)生整體呈下降趨勢，但年際間波動幅度較大，無顯著差異（P>0.05）;在各時(shí)段發(fā)生趨勢也不盡相同，其中1960—2000年輕度干熱風(fēng)發(fā)生站次整體呈波動中下降趨勢，2001年以后又呈現(xiàn)上升趨勢。重度干熱風(fēng)（圖2-b）在1960—1985年是高發(fā)期，1985年后，除2001年出現(xiàn)1次極值和2014年1次相對高值外，其他年份發(fā)生均較少;其整體發(fā)生趨勢與輕度干熱風(fēng)均呈下降趨勢，但下降趨勢更為明顯，有顯著差異（P<0.05）。干熱風(fēng)總站次（圖2-c）發(fā)生趨勢與輕度干熱風(fēng)的變化規(guī)律更為接近，同樣呈下降趨勢，無顯著差異（P>0.05）。近57年來，僅1968、1972、2014年等3年重度干熱風(fēng)占比超過50%，其他年份均在50%以下（圖2-d），最低的2006年僅占干熱風(fēng)總站次的1.3%;從整體上看，重度干熱風(fēng)所占比重呈下降趨勢，但同樣無顯著差異（P>0.05）。

2.1.2 干熱風(fēng)的突變分析 對所選站點(diǎn)近57年來干熱風(fēng)發(fā)生站次變化進(jìn)行M-K檢驗(yàn)及滑動t檢驗(yàn)，結(jié)果見圖3。由圖3-a可見，在1960—2000年，河北省輕度干熱風(fēng)發(fā)生呈波動下降趨勢;其中，1960—1970年波動幅度較大，1971—1979年整體較平緩，1982年后開始迅速下降，2000年達(dá)到最低點(diǎn)，但沒有超過-1.96（P<0.05）的臨界值。2001年UF逐漸上升，說明2001年開始輕度干熱風(fēng)加重發(fā)生，一直到2015年，但UF依然沒有超過1.96（P<0.05）的臨界值。因此，輕度干熱風(fēng)的發(fā)生在近57年來并沒有發(fā)生突變。圖3-b為滑動t檢驗(yàn)的曲線圖驗(yàn)證了前面的觀點(diǎn)，輕度干熱風(fēng)的下降或者上升均沒有達(dá)到顯著性水平。重度干熱風(fēng)UF變化與輕度干熱風(fēng)稍有差異（圖3-c），1972年開始，UF有明顯的下降趨勢，1985年后更為明顯，并在1997年超過-1.96（P<0.05）的臨界值，且UF、UB曲線相交于1985年，說明重度干熱風(fēng)的發(fā)生在1985年出現(xiàn)了突變性下降。圖3-d也驗(yàn)證了前面的觀點(diǎn)，滑動t值在1986年前后達(dá)到最大值，且超過了臨界值。2001年開始，干熱風(fēng)發(fā)生出現(xiàn)加重趨勢，這與前人的研究結(jié)果[4]一致，但這種趨勢并不明顯，特別是在2016年，全省干熱風(fēng)災(zāi)害發(fā)生較輕，重度干熱風(fēng)僅出現(xiàn)1站次，發(fā)生總站次也屬于歷史較少年份。1996年開始這種減少趨勢，并通過95%的信度檢驗(yàn)。UF與UB曲線相交于1985年，可以認(rèn)為該年為干熱風(fēng)出現(xiàn)站次發(fā)生突變的年份。干熱風(fēng)發(fā)生總數(shù)M-K檢驗(yàn)結(jié)果如圖3-e所示，UF、UB曲線相交于1976年，但在滑動t檢驗(yàn)中，該時(shí)段滑動t值并沒有出現(xiàn)峰值（圖3-f），由此認(rèn)為該點(diǎn)并不是真正的突變點(diǎn)，干熱風(fēng)發(fā)生總數(shù)在近57年來沒有發(fā)生突變。

2.1.3 發(fā)生時(shí)段分布 干熱風(fēng)發(fā)生時(shí)段對冬小麥灌漿的影響也各不相同[18]，張志紅等發(fā)現(xiàn)，干熱風(fēng)對冬小麥生理機(jī)能的危害程度表現(xiàn)為灌漿中期>灌漿前期[2]。5月中旬、5月下旬、6月上旬分別對應(yīng)河北省冬小麥的灌漿前期、中期和后期。因此，從旬尺度分析研究區(qū)域內(nèi)干熱風(fēng)發(fā)生時(shí)段分布，對綜合評估干熱風(fēng)發(fā)生對冬小麥灌漿的影響有重要意義。由圖4可見，隨著時(shí)間的推遲，輕、重干熱風(fēng)的發(fā)生均呈上升趨勢，重度干熱風(fēng)上升趨勢更為明顯。5月下旬輕、重干熱風(fēng)發(fā)生次數(shù)均最少，分別占總站次的18.9%、11.8%;5月中旬輕、重干熱風(fēng)發(fā)生次數(shù)分別占總站次的40.2%、32.3%;6月上旬輕、重干熱風(fēng)發(fā)生站次均為時(shí)段內(nèi)最高，分別占總站次的40.9%、55.9%。超過50%的重度干熱風(fēng)發(fā)生在冬小麥的灌漿后期，灌漿中期次之，灌漿前期不管是輕度干熱風(fēng)還是重度干熱風(fēng)，發(fā)生的概率均最小。

2.2 空間分布

根據(jù)表1中干熱風(fēng)氣象行業(yè)標(biāo)準(zhǔn)（QX/T 82—2007《小麥干熱風(fēng)災(zāi)害等級》）計(jì)算出各站點(diǎn)年均干熱風(fēng)日數(shù)，結(jié)果見表2。輕、重干熱風(fēng)的分布具有比較高的一致性，最大值均出現(xiàn)在南宮市，年均輕、重干熱風(fēng)發(fā)生次數(shù)分別為3.5、2.2次，同區(qū)域的肥鄉(xiāng)、內(nèi)丘、涉縣等地也是干熱風(fēng)發(fā)生的重災(zāi)區(qū)。東北部為干熱風(fēng)發(fā)生較少區(qū)域，其中唐海輕、重干熱風(fēng)年均發(fā)生次數(shù)僅為0.5、0.1次，同區(qū)域的昌黎、青龍等干熱風(fēng)發(fā)生也均較少。

根據(jù)各站點(diǎn)的計(jì)算結(jié)果，按其站點(diǎn)所在的經(jīng)緯度，在ArcGIS中以反距離權(quán)重法進(jìn)行插值，得到所選區(qū)域的年均干熱風(fēng)日數(shù)[圖5（承德、張家口為非種植區(qū)）]。河北省冬麥區(qū)年均輕度（圖5-a）、重度（圖5-b）干熱風(fēng)日數(shù)分布一致，大體呈東北部低、南部高的趨勢。河北省邯鄲、邢臺的太行山山前地帶以及黑龍港流域的中西部為干熱風(fēng)的重災(zāi)區(qū)，東北部的唐山、秦皇島地區(qū)為低值區(qū);南宮、肥鄉(xiāng)、黃驊一帶是干熱風(fēng)的高發(fā)區(qū)域。本研究所選24個農(nóng)業(yè)氣象觀測站點(diǎn)57年來總計(jì)發(fā)生輕度干熱風(fēng)2 992站次，平均每年每站2.18次;總計(jì)發(fā)生重度干熱風(fēng)1 399站次，平均每年每站1.02次。

3 結(jié)論與討論

近57年來，河北省干熱風(fēng)發(fā)生整體呈下降趨勢，重度干熱風(fēng)的下降趨勢更為明顯，通過了顯著性檢驗(yàn)（P<0.05）。輕度干熱風(fēng)和干熱風(fēng)總次數(shù)的下降并沒有通過顯著性檢驗(yàn)（P>0.05）。從年際尺度上看，河北省1960—1990前期是干熱風(fēng)的高發(fā)時(shí)期，80年代中后期到90年代末，干熱風(fēng)發(fā)生較少，2001年出現(xiàn)歷史極大值，發(fā)生干熱風(fēng)共231站次，這次干熱風(fēng)災(zāi)害發(fā)生突然性強(qiáng)，波及范圍廣。此后，干熱風(fēng)的發(fā)生進(jìn)入一個短暫的低谷期，之后又呈緩慢上升的趨勢。

根據(jù)M-K檢驗(yàn)及滑動t檢驗(yàn)結(jié)果，輕度干熱風(fēng)和干熱風(fēng)總數(shù)在近57年來并沒有發(fā)生突變，而重度干熱風(fēng)的發(fā)生在20世紀(jì)80年代中期出現(xiàn)了突變性下降。

從區(qū)域上看，河北省麥區(qū)各站點(diǎn)的年均輕度、重度干熱風(fēng)

日數(shù)分布一致，大體呈東北部低、南部高的趨勢。河北省邯鄲、邢臺的太行山山前地帶以及黑龍港流域的中西部為干熱風(fēng)的重災(zāi)區(qū)，東北部的唐山、秦皇島地區(qū)為低值區(qū);南宮、肥鄉(xiāng)、黃驊一帶是干熱風(fēng)的高發(fā)區(qū)域。

氣象指標(biāo)雖然可以作為判別干熱風(fēng)發(fā)生的依據(jù)，但其危害程度卻不僅僅是由氣象指標(biāo)決定的。還與土壤墑情、小麥發(fā)育期、品種、干熱風(fēng)持續(xù)時(shí)間、農(nóng)田小氣候等有密切聯(lián)系[19-21]。在小麥籽粒中積累的干物質(zhì)中，有70%～80%來自于開花期至灌漿成熟期間的生育后期[22]，通過對河北省冬小麥干熱風(fēng)時(shí)空分布規(guī)律的研究，可以提高干熱風(fēng)預(yù)警的準(zhǔn)確性和可靠性，對于冬小麥干熱風(fēng)災(zāi)害的防災(zāi)減災(zāi)有著重要的現(xiàn)實(shí)意義。

參考文獻(xiàn)：

[1]趙風(fēng)華，居 輝，歐陽竹. 干熱風(fēng)對灌漿期冬小麥旗葉光合蒸騰的影響[J]. 華北農(nóng)學(xué)報(bào)，2013，28（5）：144-148.

[2]張志紅，成 林，李書嶺，等. 干熱風(fēng)天氣對冬小麥的生理影響[J]. 生態(tài)學(xué)雜志，2015，34（3）：712-717.

[3]成 林，張志紅，方文松. 干熱風(fēng)對冬小麥灌漿速率和千粒重的影響[J]. 麥類作物學(xué)報(bào)，2014（2）：248-254.