農(nóng)業(yè)氣候年景評估指數(shù)構(gòu)建及在江漢平原的應(yīng)用

秦鵬程 夏智宏 陳偉亮

（1 武漢區(qū)域氣候中心，武漢 430074；2 荊州市氣象局，荊州 434020；3 江漢平原生態(tài)氣象遙感監(jiān)測技術(shù)協(xié)同創(chuàng)新中心, 荊州 434025）

0 引言

氣候年景評估是對一年（或一個生產(chǎn)周期）氣候要素異常程度及其對國民經(jīng)濟影響利弊的綜合評價，是開展氣候以及氣候變化對敏感行業(yè)影響評估的基礎(chǔ)[1]。農(nóng)業(yè)是一個高度依賴于天氣氣候條件的敏感性行業(yè)，農(nóng)業(yè)氣候年景的評估是氣候評估業(yè)務(wù)的重要內(nèi)容之一，而作物播種前開展的農(nóng)業(yè)氣候年景預(yù)測和展望可為農(nóng)業(yè)生產(chǎn)布局安排、資源配置和災(zāi)害防御提供重要參考。

20世紀90年代以來，為了探索定量化的作物年景評價方法，相關(guān)學者相繼以作物產(chǎn)量豐歉、綜合氣象要素、氣象災(zāi)害災(zāi)情等為指標建立了評價模型和指標。如王建林等[2]綜合逐旬氣象要素建立了棉花年景評價方法，婁秀榮等[3]以各生育期氣象要素距平為指標建立了晚稻氣候年景評價方法，孫家民等[4]利用降水距平絕對值累計方法，從旱澇災(zāi)害等級的角度開展了中國農(nóng)業(yè)氣候年景評估研究，吳義城等[5]、劉敏等[6]先后對湖北省旱澇與中稻、早稻、棉花年景進行了系統(tǒng)分析，吳義城等[7]以作物生育期內(nèi)的氣溫、積溫、日照及降水為指標，從氣候異常度角度建立了農(nóng)作物氣候年型劃分方法，石磊等[1]、王勝等[8]分別以氣象災(zāi)害和氣候適宜度指標構(gòu)建了安徽省氣候年景和作物年景評估指標。上述研究為開展農(nóng)業(yè)氣候年景評價提供了重要依據(jù)，但各類方法局限于單方面因子的年際變異描述，以產(chǎn)量豐歉為依據(jù)評價氣候年景，忽略了農(nóng)業(yè)生產(chǎn)技術(shù)水平提高及防災(zāi)救災(zāi)措施對不利氣象條件的補償作用，可能低估氣候的不利影響；以氣象災(zāi)害災(zāi)情等級為基礎(chǔ)的農(nóng)業(yè)氣候年景評價，考慮的災(zāi)害類型尚不全面，且未反映光、溫、水等氣候資源的優(yōu)劣。

綜合前人研究成果，農(nóng)業(yè)氣候年景的內(nèi)涵應(yīng)至少包括農(nóng)業(yè)氣候資源優(yōu)劣和氣象災(zāi)害發(fā)生程度兩個方面，為此，本文以光、溫、水氣候適宜度為指標，建立作物生長季農(nóng)業(yè)氣候資源等級評價指標，以旱、澇、冷、熱等氣象災(zāi)害指數(shù)為指標，建立農(nóng)業(yè)氣象災(zāi)害等級評價指標，并通過有效綜合建立農(nóng)業(yè)氣候年景綜合指數(shù)。最后，以江漢平原為例，對農(nóng)業(yè)氣候年景評估指數(shù)在氣候變化評估以及業(yè)務(wù)服務(wù)中的應(yīng)用進行了案例分析。

1 農(nóng)業(yè)氣候年景指數(shù)構(gòu)建

1.1 農(nóng)業(yè)氣候資源評價

1.1.1 氣候適宜度模型

根據(jù)農(nóng)作物生長發(fā)育對氣候資源的需求和生理特點，可采用模糊數(shù)學理論中的隸屬函數(shù)對農(nóng)作物氣候適宜度及氣候資源有效性進行量化，參考已有研究成果[8-12]，分別建立氣溫、降水和日照的氣候適宜度模型如下

式中，S(T)為作物生長發(fā)育的溫度適宜度，T為日平均氣溫，Ttop、Tbase、Topt分別為作物生長發(fā)育所需的上限溫度、下限溫度和最適溫度。B為形狀參數(shù)，計算公式為

式中，S(R)為作物生長發(fā)育的降水適宜度，R0為作物需水量的臨界值，Re為當前日有效降水量，是對前期降水量的加權(quán)累積值，計算公式為[13]

式中，n為距離當前的日數(shù)，Rn為當前日期前第n天的日降水量，wn為權(quán)重系數(shù)，可近似表示為為降水衰減參數(shù)，取值0.95。

式中，S(S)為作物生長發(fā)育的日照適宜度，S0為作物生長發(fā)育達到適宜狀態(tài)所需的臨界日照時數(shù)，b為形狀參數(shù)。

S(T)、S(R)、S(S)的形狀分布如圖1所示，其值域界于0～1，因此不同氣候要素的適宜度具有可比性。受作物三基點溫度生理特性的影響，當平均氣溫低于或高于最適溫度時，熱量條件的有效性下降，因此溫度適宜度曲線呈拋物線形式，降水和日照適宜度則呈現(xiàn)逐漸增加并趨于穩(wěn)定的形式，當降水或日照低于最適需求量時其適宜度為小于1的數(shù)值，并逐漸升高，達到最適需求量以上時其適宜度為1。

圖1 農(nóng)作物溫度、降水、日照氣候適宜度曲線Fig.1 The climatic suitability curves of crop for temperature, precipitation, and sunshine duration

1.1.2 氣候資源評價

為了綜合評價光、溫、水氣候資源及其匹配程度，采用幾何平均法對各氣候要素單因子適宜度求取平均值構(gòu)建綜合氣候適宜度模型[11]，計算公式為

為了綜合評價某地一年氣候資源總體狀況，根據(jù)該地區(qū)農(nóng)業(yè)種植結(jié)構(gòu)，選擇主要農(nóng)作物及其生長發(fā)育時段，分作物分別計算各作物每日溫度、降水、日照適宜度及其綜合氣候適宜度，然后以作物種植面積比例為權(quán)重，計算該地每日的氣候適宜度，最后將一年內(nèi)各日氣候適宜度平均值作為該年度的氣候適宜度。

通過對歷年氣候適宜度進行排序，采用百分位法建立氣候資源等級劃分依據(jù)，分別按照25%和75%百分位數(shù)對應(yīng)的適宜度值，將氣候資源年景劃分為偏差、一般、偏好年景。

1.2 農(nóng)業(yè)氣象災(zāi)害年景評價

綜合考慮作物生長季內(nèi)影響生長發(fā)育、產(chǎn)量形成以及農(nóng)事活動開展的旱、澇、冷、熱等氣象災(zāi)害類型，利用相關(guān)農(nóng)業(yè)氣象災(zāi)害指標或模型建立主要農(nóng)業(yè)氣象災(zāi)害發(fā)生程度的量化表達指數(shù)（表1），統(tǒng)計歷年各類氣象災(zāi)害過程強度（年內(nèi)多次過程取各過程強度之和），通過經(jīng)驗累積頻率轉(zhuǎn)換或離差標準化將其歸一化為0～1的無量綱指標，以歷年實際災(zāi)情統(tǒng)計資料為基礎(chǔ)，根據(jù)農(nóng)作物受災(zāi)面積確定不同災(zāi)害的權(quán)重系數(shù)，最后將歸一化后的各災(zāi)害指數(shù)進行加權(quán)平均，構(gòu)建農(nóng)業(yè)氣象災(zāi)害年景指數(shù)。采用百分位法建立農(nóng)業(yè)氣象災(zāi)害年景等級劃分依據(jù)，分別按照25%、75%百分位數(shù)將氣象災(zāi)害年景劃分為偏輕、一般和偏重年景。

表1 主要農(nóng)業(yè)氣象災(zāi)害及氣象指標Table 1 Main agrometeorological disasters and meteorological indices

1.3 農(nóng)業(yè)氣候年景綜合評價

對農(nóng)業(yè)氣候資源年景和氣象災(zāi)害年景等級分別進行打分：一般年景為0，偏好或偏輕年景為1，偏差或偏重年景為－1，將兩類年景得分累積，若累積得分等于2，則認定該年度農(nóng)業(yè)氣候年景等級為較好年景，累積得分等于1為一般偏好年景，累積得分等于0為一般年景，累積得分等于－1為一般偏差年景，否則為偏差年景。

2 江漢平原農(nóng)業(yè)氣候年景評估應(yīng)用

江漢平原地處長江中游，由長江及其支流漢江沖積而成，土壤肥沃，水源充足，氣候溫和，農(nóng)業(yè)生產(chǎn)自然條件優(yōu)越，是我國重要的糧、棉、油、魚生產(chǎn)基地，素有“魚米之鄉(xiāng)”之稱。同時，因受季風氣候以及地形、地理位置制約，雨季降水集中、洪澇災(zāi)害多發(fā)，盛夏在副熱帶高壓控制下，常出現(xiàn)持續(xù)晴熱高溫天氣，造成農(nóng)作物高溫熱害，春、秋季易受冷空氣南下影響，常發(fā)生陰雨濕漬災(zāi)害，農(nóng)業(yè)年景因氣候資源和氣象災(zāi)害的程度不同而呈現(xiàn)較大的年際變異，是制約糧食生產(chǎn)穩(wěn)定性的主要限制因子之一。

根據(jù)江漢平原分縣作物種植面積及產(chǎn)量統(tǒng)計資料，選擇油菜、小麥、早稻、中稻、晚稻、玉米、棉花等主要農(nóng)作物為農(nóng)業(yè)氣候年景評估對象，利用區(qū)域內(nèi)農(nóng)業(yè)氣象觀測站作物發(fā)育期觀測資料確定各作物的平均發(fā)育期，結(jié)合《大宗作物氣象服務(wù)手冊》[19]、《湖北省公共氣象服務(wù)手冊》[17]以及前人相關(guān)研究成果[8-11]，確定各作物氣候適宜度模型的參數(shù)以及相關(guān)災(zāi)害指標，利用區(qū)域內(nèi)16個氣象臺站1961—2017年逐日氣象觀測資料，分別計算歷年氣候資源和氣象災(zāi)害年景指數(shù)，在此基礎(chǔ)上通過區(qū)域平均和等級劃分建立江漢平原農(nóng)業(yè)氣候年景序列。

2.1 農(nóng)業(yè)氣候資源評價

圖2為江漢平原1961—2017年逐年溫度、降水、日照適宜度序列，從圖中可以看出，適宜度指數(shù)的年際間具有明顯波動，能夠反映氣候資源優(yōu)劣的年際差異，其中溫度適宜度數(shù)值介于0.58～0.7，降水適宜度介于0.6～0.9，日照適宜度介于0.45～0.65，三者中以降水適宜度的變化幅度最大，這是由于降水的年際變異大，且降水資源不足引起水分虧缺通常會持續(xù)一定的時間，日照和溫度適宜度年際變化相對較小，表明本地區(qū)熱量和光照資源相對穩(wěn)定。從圖中還可以看出，溫度和日照適宜度在近60年中呈現(xiàn)明顯的趨勢性變化，這與近百年全球氣候變化的背景有關(guān)，研究表明近60年來江漢平原氣溫呈顯著上升趨勢，與此同時熱量資源有所改善，但地表總輻射呈顯著下降趨勢，導(dǎo)致光照資源有減少趨勢。綜合光、溫、水及其匹配狀況的氣候適宜度指數(shù)序列（如圖3a所示），并基于百分位數(shù)劃分農(nóng)業(yè)氣候年景等級（圖3b），統(tǒng)計表明1961—2017年農(nóng)業(yè)氣候資源共計有偏好年景16年，偏差年景15年，一般年景26年，不同年代間20世紀70年代和21世紀以來，以偏差或一般年景為主，20世紀90年代以偏好年景為主，20世紀80年代氣候資源變化幅度最大。可見，在氣候變暖背景下，盡管江漢平原熱量資源有所增加，但由于光照資源的下降以及降水資源的年際變異，導(dǎo)致農(nóng)業(yè)氣候資源總量及其匹配協(xié)調(diào)的有效性下降，表明氣候變化對江漢平原農(nóng)業(yè)氣候資源的影響總體以不利為主。

圖2 1961—2017年江漢平原農(nóng)作物溫度、降水、日照適宜度變化Fig.2 Temporal changes in climatic suitability of temperature, precipitation, and sunshine for the Jianghan Plain from 1961 to 2017

圖3 1961—2017年江漢平原（a）農(nóng)業(yè)氣候資源綜合指數(shù)時間演變及（b）年景等級劃分Fig.3 Temporal change of (a) comprehensive index of agroclimate resources and (b) classification of patterns by year for the Jianghan Plain from 1961 to 2017

2.2 農(nóng)業(yè)氣象災(zāi)害年景

圖4是基于表1農(nóng)業(yè)氣象災(zāi)害指標和模型反演的1961—2017年江漢平原主要農(nóng)業(yè)氣象災(zāi)害時間序列，從圖中可以看出，江漢平原地區(qū)發(fā)生頻率較高的農(nóng)業(yè)氣象災(zāi)害主要有暴雨洪澇、干旱、濕漬害、連陰雨、高溫熱害、倒春寒等，幾乎每年均有發(fā)生，五月寒、寒露風、盛夏低溫及低溫凍害發(fā)生頻率相對較低。近幾十年來在氣候變暖背景下，低溫冷害、凍害的發(fā)生頻率明顯減少，高溫熱害的頻率和強度有所增加，春、秋季連陰雨和濕漬害的強度在21世紀以來有所減弱，暴雨洪澇和干旱的發(fā)生頻率和強度沒有明顯的趨勢性變化。圖5為農(nóng)業(yè)氣象災(zāi)害綜合指數(shù)的時間序列及其年景等級，從圖中可以看出，20世紀60年代至70年代前期，氣象災(zāi)害年景以一般和偏差年景為主，期間寒露風、低溫凍害和連陰雨災(zāi)害較重，20世紀70年代后期至90年代前期，氣象災(zāi)害年景的年際變幅較大，期間出現(xiàn)了1961年以來最嚴重的暴雨洪澇（1983年）和濕漬害（1991年），21世紀以來農(nóng)業(yè)氣象災(zāi)害年景以一般或偏好年景為主，這與低溫冷害、凍害及濕漬害、連陰雨的強度減弱有關(guān)，但部分年份（如2016年）由于個別氣象災(zāi)害異常偏重，氣象災(zāi)害年景總體偏差，且強度居歷史前列。

2.3 農(nóng)業(yè)氣候年景綜合評價

圖4 1961—2017年江漢平原主要農(nóng)業(yè)氣象災(zāi)害時間演變Fig.4 Temporal changes of main agrometeorological disasters

圖5 1961—2017年（a）江漢平原農(nóng)業(yè)氣象災(zāi)害綜合指數(shù)時間演變及（b）年景等級劃分Fig.5 Temporal change of (a) comprehensive index of agrometeorological disasters and (b) classification of patterns by year for the Jianghan Plain from 1961 to 2017

圖6 1961—2017年江漢平原農(nóng)業(yè)氣候年景等級Fig.6 Agroclimatic patterns by year for the Jianghan Plain from 1961 to 2017

綜合農(nóng)業(yè)氣候資源和農(nóng)業(yè)氣象災(zāi)害年景建立農(nóng)業(yè)氣候年景等級（如圖6所示），江漢平原1961—2017年較好年景有3年（1965，1975，1993年），上述年份農(nóng)業(yè)氣候資源總量和匹配程度較好，農(nóng)業(yè)氣象災(zāi)害發(fā)生程度總體偏輕，較差年份有3年（1988，1991，2003年），上述年份農(nóng)業(yè)氣候資源總量和匹配程度較差，且農(nóng)業(yè)氣象災(zāi)害發(fā)生程度偏重，一般偏好和一般偏差年景各有15年，一般年景有21年。以2017年為例，綜合圖2～6可以看出，該年度溫度適宜度處于歷史較高水平，熱量條件較好，降水適宜度為一般偏好等級，基本滿足農(nóng)業(yè)用水需求，但日照適宜度接近歷史最低水平，日照條件異常偏差，農(nóng)業(yè)氣候資源總體屬一般年景，年內(nèi)高溫熱害和連陰雨災(zāi)害較為突出，部分時段還出現(xiàn)了暴雨洪澇、濕漬害、干旱、倒春寒、五月寒、寒露風和盛夏低溫災(zāi)害，但發(fā)生程度相對較輕，農(nóng)業(yè)氣象災(zāi)害年景總體屬一般年景，綜合農(nóng)業(yè)氣候資源狀況和氣象災(zāi)害發(fā)生程度，2017年江漢平原農(nóng)業(yè)氣候年景為一般年景。

3 結(jié)論與討論

本文從農(nóng)業(yè)氣候資源優(yōu)劣和農(nóng)業(yè)氣象災(zāi)害發(fā)生程度兩個方面，基于氣候適宜度和氣象災(zāi)害指標分別構(gòu)建指標體系，并采用百分位數(shù)進行等級劃分，建立了農(nóng)業(yè)氣候年景評估指標體系和綜合指數(shù)。該指數(shù)綜合反映了農(nóng)業(yè)生產(chǎn)的氣候資源總量和匹配程度，以及不利氣象災(zāi)害的影響程度，物理意義清晰、易于擴展和應(yīng)用，可作為農(nóng)業(yè)氣象及生態(tài)氣象業(yè)務(wù)服務(wù)指標。

以江漢平原為例，利用農(nóng)業(yè)氣候年景評估指標體系對1961—2017年江漢平原農(nóng)業(yè)氣候資源、農(nóng)業(yè)氣象災(zāi)害以及農(nóng)業(yè)氣候綜合年景的時間演變特征分析表明，江漢平原農(nóng)業(yè)氣候資源、氣象災(zāi)害發(fā)生程度具有明顯的年際和年代際變化，在氣候變暖背景下，江漢平原熱量資源有所增加，降水資源變化不明顯，但光照資源明顯減少，農(nóng)業(yè)氣候資源總量及其匹配協(xié)調(diào)的有效性下降；隨著氣候變暖，江漢平原倒春寒、五月寒、寒露風、盛夏低溫及低溫凍害等農(nóng)業(yè)氣象災(zāi)害發(fā)生頻率明顯減少，高溫熱害的頻率和強度有所增加；由于氣候變化對農(nóng)業(yè)氣候資源和氣象災(zāi)害影響的復(fù)雜性，近60年來江漢平原農(nóng)業(yè)氣候綜合年景雖有波動，但未呈現(xiàn)明顯的趨勢性變化。以2017年為例，對開展任意年份農(nóng)業(yè)氣候年景評估示例，該年度熱量條件較好，降水資源一般，光照資源偏差，農(nóng)業(yè)氣候資源屬一般年景，年內(nèi)出現(xiàn)嚴重的高溫熱害和連陰雨災(zāi)害，但暴雨洪澇、濕漬害、干旱等其他災(zāi)害發(fā)生程度較輕，農(nóng)業(yè)氣象災(zāi)害總體為一般年景，綜合農(nóng)業(yè)氣候資源狀況和氣象災(zāi)害發(fā)生程度，2017年江漢平原農(nóng)業(yè)氣候年景為一般年景。

如何界定農(nóng)業(yè)氣候年景的優(yōu)劣尚沒有統(tǒng)一的標準[4,20-21]。以往研究通常以產(chǎn)量豐歉作為年景劃分的依據(jù)[21-22]，但糧食產(chǎn)量是氣候條件與農(nóng)業(yè)生產(chǎn)管理多方因素綜合的結(jié)果，因此農(nóng)業(yè)氣候年景不等同于農(nóng)業(yè)年景，此外，不同農(nóng)作物產(chǎn)量水平和經(jīng)濟價值存在差異，當用于評價某地某年氣候年景的總體狀況時，如何確定不同作物的權(quán)重是年景劃分的難點，也是影響年景評價結(jié)果科學性和公眾理解、接受程度的關(guān)鍵。本文從農(nóng)業(yè)氣候年景的內(nèi)涵出發(fā)，基于氣象指標建立了年景評估和劃分方法，是對農(nóng)業(yè)氣候年景評估技術(shù)體系的有效補充，由于缺乏必要的驗證，其應(yīng)用效果尚需要在業(yè)務(wù)應(yīng)用和實踐中逐步檢驗。本文提出的指標體系不僅適用于年度氣候影響評價，同時可應(yīng)用于單作物的氣候年景評價和產(chǎn)量預(yù)報，在不同地區(qū)應(yīng)用時根據(jù)需要可進一步對指標體系進行細化和完善。