秦鵬程 夏智宏 陳偉亮

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

0 引言

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

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

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

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

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

1.1.1 氣候適宜度模型

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

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

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

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

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

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

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

1.1.2 氣候資源評(píng)價(jià)

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

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

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

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

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

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

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

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

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

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

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

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

圖2為江漢平原1961—2017年逐年溫度、降水、日照適宜度序列，從圖中可以看出，適宜度指數(shù)的年際間具有明顯波動(dòng)，能夠反映氣候資源優(yōu)劣的年際差異，其中溫度適宜度數(shù)值介于0.58～0.7，降水適宜度介于0.6～0.9，日照適宜度介于0.45～0.65，三者中以降水適宜度的變化幅度最大，這是由于降水的年際變異大，且降水資源不足引起水分虧缺通常會(huì)持續(xù)一定的時(shí)間，日照和溫度適宜度年際變化相對(duì)較小，表明本地區(qū)熱量和光照資源相對(duì)穩(wěn)定。從圖中還可以看出，溫度和日照適宜度在近60年中呈現(xiàn)明顯的趨勢(shì)性變化，這與近百年全球氣候變化的背景有關(guān)，研究表明近60年來(lái)江漢平原氣溫呈顯著上升趨勢(shì)，與此同時(shí)熱量資源有所改善，但地表總輻射呈顯著下降趨勢(shì)，導(dǎo)致光照資源有減少趨勢(shì)。綜合光、溫、水及其匹配狀況的氣候適宜度指數(shù)序列（如圖3a所示），并基于百分位數(shù)劃分農(nóng)業(yè)氣候年景等級(jí)（圖3b），統(tǒng)計(jì)表明1961—2017年農(nóng)業(yè)氣候資源共計(jì)有偏好年景16年，偏差年景15年，一般年景26年，不同年代間20世紀(jì)70年代和21世紀(jì)以來(lái)，以偏差或一般年景為主，20世紀(jì)90年代以偏好年景為主，20世紀(jì)80年代氣候資源變化幅度最大?？梢?jiàn)，在氣候變暖背景下，盡管江漢平原熱量資源有所增加，但由于光照資源的下降以及降水資源的年際變異，導(dǎo)致農(nóng)業(yè)氣候資源總量及其匹配協(xié)調(diào)的有效性下降，表明氣候變化對(duì)江漢平原農(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ù)時(shí)間演變及（b）年景等級(jí)劃分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)害指標(biāo)和模型反演的1961—2017年江漢平原主要農(nóng)業(yè)氣象災(zāi)害時(shí)間序列，從圖中可以看出，江漢平原地區(qū)發(fā)生頻率較高的農(nóng)業(yè)氣象災(zāi)害主要有暴雨洪澇、干旱、濕漬害、連陰雨、高溫?zé)岷Α⒌勾汉?，幾乎每年均有發(fā)生，五月寒、寒露風(fēng)、盛夏低溫及低溫凍害發(fā)生頻率相對(duì)較低。近幾十年來(lái)在氣候變暖背景下，低溫冷害、凍害的發(fā)生頻率明顯減少，高溫?zé)岷Φ念l率和強(qiáng)度有所增加，春、秋季連陰雨和濕漬害的強(qiáng)度在21世紀(jì)以來(lái)有所減弱，暴雨洪澇和干旱的發(fā)生頻率和強(qiáng)度沒(méi)有明顯的趨勢(shì)性變化。圖5為農(nóng)業(yè)氣象災(zāi)害綜合指數(shù)的時(shí)間序列及其年景等級(jí)，從圖中可以看出，20世紀(jì)60年代至70年代前期，氣象災(zāi)害年景以一般和偏差年景為主，期間寒露風(fēng)、低溫凍害和連陰雨災(zāi)害較重，20世紀(jì)70年代后期至90年代前期，氣象災(zāi)害年景的年際變幅較大，期間出現(xiàn)了1961年以來(lái)最嚴(yán)重的暴雨洪澇（1983年）和濕漬害（1991年），21世紀(jì)以來(lái)農(nóng)業(yè)氣象災(zāi)害年景以一般或偏好年景為主，這與低溫冷害、凍害及濕漬害、連陰雨的強(qiáng)度減弱有關(guān)，但部分年份（如2016年）由于個(gè)別氣象災(zāi)害異常偏重，氣象災(zāi)害年景總體偏差，且強(qiáng)度居歷史前列。

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

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

圖5 1961—2017年（a）江漢平原農(nóng)業(yè)氣象災(zāi)害綜合指數(shù)時(shí)間演變及（b）年景等級(jí)劃分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è)氣候年景等級(jí)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è)氣候年景等級(jí)（如圖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可以看出，該年度溫度適宜度處于歷史較高水平，熱量條件較好，降水適宜度為一般偏好等級(jí)，基本滿足農(nóng)業(yè)用水需求，但日照適宜度接近歷史最低水平，日照條件異常偏差，農(nóng)業(yè)氣候資源總體屬一般年景，年內(nèi)高溫?zé)岷瓦B陰雨災(zāi)害較為突出，部分時(shí)段還出現(xiàn)了暴雨洪澇、濕漬害、干旱、倒春寒、五月寒、寒露風(fēng)和盛夏低溫災(zāi)害，但發(fā)生程度相對(duì)較輕，農(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ā)生程度兩個(gè)方面，基于氣候適宜度和氣象災(zāi)害指標(biāo)分別構(gòu)建指標(biāo)體系，并采用百分位數(shù)進(jìn)行等級(jí)劃分，建立了農(nóng)業(yè)氣候年景評(píng)估指標(biāo)體系和綜合指數(shù)。該指數(shù)綜合反映了農(nóng)業(yè)生產(chǎn)的氣候資源總量和匹配程度，以及不利氣象災(zāi)害的影響程度，物理意義清晰、易于擴(kuò)展和應(yīng)用，可作為農(nóng)業(yè)氣象及生態(tài)氣象業(yè)務(wù)服務(wù)指標(biāo)。

以江漢平原為例，利用農(nóng)業(yè)氣候年景評(píng)估指標(biāo)體系對(duì)1961—2017年江漢平原農(nóng)業(yè)氣候資源、農(nóng)業(yè)氣象災(zāi)害以及農(nóng)業(yè)氣候綜合年景的時(shí)間演變特征分析表明，江漢平原農(nóng)業(yè)氣候資源、氣象災(zāi)害發(fā)生程度具有明顯的年際和年代際變化，在氣候變暖背景下，江漢平原熱量資源有所增加，降水資源變化不明顯，但光照資源明顯減少，農(nóng)業(yè)氣候資源總量及其匹配協(xié)調(diào)的有效性下降；隨著氣候變暖，江漢平原倒春寒、五月寒、寒露風(fēng)、盛夏低溫及低溫凍害等農(nóng)業(yè)氣象災(zāi)害發(fā)生頻率明顯減少，高溫?zé)岷Φ念l率和強(qiáng)度有所增加；由于氣候變化對(duì)農(nóng)業(yè)氣候資源和氣象災(zāi)害影響的復(fù)雜性，近60年來(lái)江漢平原農(nóng)業(yè)氣候綜合年景雖有波動(dòng)，但未呈現(xiàn)明顯的趨勢(shì)性變化。以2017年為例，對(duì)開(kāi)展任意年份農(nóng)業(yè)氣候年景評(píng)估示例，該年度熱量條件較好，降水資源一般，光照資源偏差，農(nóng)業(yè)氣候資源屬一般年景，年內(nèi)出現(xiàn)嚴(yán)重的高溫?zé)岷瓦B陰雨災(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)劣尚沒(méi)有統(tǒng)一的標(biāo)準(zhǔn)[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)濟(jì)價(jià)值存在差異，當(dāng)用于評(píng)價(jià)某地某年氣候年景的總體狀況時(shí)，如何確定不同作物的權(quán)重是年景劃分的難點(diǎn)，也是影響年景評(píng)價(jià)結(jié)果科學(xué)性和公眾理解、接受程度的關(guān)鍵。本文從農(nóng)業(yè)氣候年景的內(nèi)涵出發(fā)，基于氣象指標(biāo)建立了年景評(píng)估和劃分方法，是對(duì)農(nóng)業(yè)氣候年景評(píng)估技術(shù)體系的有效補(bǔ)充，由于缺乏必要的驗(yàn)證，其應(yīng)用效果尚需要在業(yè)務(wù)應(yīng)用和實(shí)踐中逐步檢驗(yàn)。本文提出的指標(biāo)體系不僅適用于年度氣候影響評(píng)價(jià)，同時(shí)可應(yīng)用于單作物的氣候年景評(píng)價(jià)和產(chǎn)量預(yù)報(bào)，在不同地區(qū)應(yīng)用時(shí)根據(jù)需要可進(jìn)一步對(duì)指標(biāo)體系進(jìn)行細(xì)化和完善。