魏忠光 崔修來 孫瑤 王曉東

(1 安徽海峰分析測試科技有限公司，合肥 230031；2 遼寧省營口市氣象局，營口 115001；3 安徽省農(nóng)業(yè)氣象中心，合肥 230031)

引言

澇漬災(zāi)害主要是指由于長時間維持陰雨天氣，降水量過多，地下水位過高，地面排水不暢，造成植物根層土壤水分長時間處于飽和或過飽和、土壤中空氣含量不足，危害植物正常生長的一種土壤水分過多的農(nóng)業(yè)氣象災(zāi)害。其中田間作物澇漬主要表現(xiàn)出土壤相對濕度(土壤含水量占田間持水量的百分比)長時間維持在85%以上、植株間的空氣相對濕度大，導(dǎo)致作物生長不良和病蟲草害嚴(yán)重等[1-5]。確定澇漬災(zāi)害指標(biāo)的方法前人已作過較多研究，也建立了許多定量指標(biāo)和模型，針對夏玉米澇漬這種特殊的災(zāi)害，國內(nèi)在指標(biāo)方法的確定以及定量指標(biāo)和模型等方面開展了相關(guān)的研究，很好的揭示了不同區(qū)域玉米澇漬災(zāi)害時空發(fā)生規(guī)律，對科學(xué)開展區(qū)域種植布局具有重要參考意義[6-12]。另外采用風(fēng)險評估指標(biāo)(包括受災(zāi)減產(chǎn)率出現(xiàn)的概率、風(fēng)險指數(shù)、變異系數(shù)等)來研究澇漬害產(chǎn)量災(zāi)損程度，并利用風(fēng)險分析技術(shù)探討作物產(chǎn)量災(zāi)損規(guī)律，對制定防災(zāi)減災(zāi)決策，降低災(zāi)損風(fēng)險程度，實現(xiàn)穩(wěn)產(chǎn)增收具有重要意義[13-17]。

由于安徽省地處暖溫帶和北亞熱帶兩大氣候帶的過渡區(qū)，過渡型氣候特征明顯，南北冷暖氣團(tuán)交匯頻繁，天氣多變，降水年際變率大，旱澇自然災(zāi)害頻發(fā)，降水集中時段與夏玉米的生長季基本吻合，期間作物水分基本處于供大于求的狀況，澇漬災(zāi)害發(fā)生概率較大，嚴(yán)重影響玉米的籽粒灌漿、結(jié)實粒數(shù)和千粒重等產(chǎn)量構(gòu)成因素，因此澇漬是安徽省夏玉米穩(wěn)產(chǎn)、高產(chǎn)的主要限制因子[18-21]。本文綜合考慮夏玉米澇漬孕災(zāi)環(huán)境、澇漬致災(zāi)因子(澇漬害發(fā)生時期的降水量、雨日、日照時數(shù))等因素，利用夏玉米澇漬指數(shù)模型對澇漬的時空分布特征進(jìn)行了詳細(xì)的分析，同時參考國家氣象行業(yè)標(biāo)準(zhǔn)《冬小麥油菜澇漬等級》[22]，確定了適合安徽省的夏玉米澇漬害分級指標(biāo);另外為揭示夏玉米澇漬可能造成的災(zāi)損風(fēng)險, 采用災(zāi)損風(fēng)險指數(shù)模型計算得到了安徽省夏玉米產(chǎn)量澇漬災(zāi)損風(fēng)險分布[23-24]，研究結(jié)果對開展區(qū)域農(nóng)業(yè)氣象防災(zāi)減災(zāi)有很好的借鑒意義。

1 研究資料與方法

1.1 研究資料

本文所用氣象資料為安徽省1971—2018年78個氣象觀測站的逐日平均氣溫、最高氣溫、最低氣溫、降水量、日照時數(shù)、水汽壓、平均風(fēng)速、相對濕度等;夏玉米發(fā)育期及土壤水分資料均來源于安徽省氣象信息中心；夏玉米產(chǎn)量資料來源于安徽省統(tǒng)計部門，包括78個站點1981—2017年的夏玉米分縣產(chǎn)量數(shù)據(jù)。研究區(qū)內(nèi)根據(jù)氣候、地形和農(nóng)業(yè)種植等可將其分為淮北地區(qū)、沿淮地區(qū)、江淮地區(qū)、大別山區(qū)、沿江地區(qū)和江南地區(qū)共6個不同區(qū)域，每個區(qū)域分別包含16個、12個、14個、4個、19個和13個氣象站點(圖1)[25]。

圖1 安徽省主要種植分區(qū)及氣象站點分布

1.2 研究方法

1.2.1 夏玉米水分敏感系數(shù)的計算

作物各生育階段對水分脅迫的影響程度不同，一般用水分敏感系數(shù)來表示，主要由作物實際產(chǎn)量、最大產(chǎn)量、實際蒸散量、參考作物蒸散量、作物系數(shù)和土壤水分系數(shù)等參數(shù)決定，各階段水分敏感系數(shù)的計算方法詳見文獻(xiàn)[26]。安徽省夏玉米一般6 月上中旬播種，9月中下旬收獲，播種較早的江南地區(qū)和播種較晚的淮北地區(qū)生育進(jìn)程普遍相差1旬左右，充分考慮夏玉米南北生育進(jìn)程的差異，本文首先計算獲得了安徽省夏玉米各站點的逐日水分敏感系數(shù)，基于逐日水分敏感系數(shù)通過統(tǒng)計分析得到了分區(qū)逐旬作物水分敏感系數(shù)，最后將分旬水分敏感系數(shù)作為各區(qū)計算全生育期夏玉米澇漬指數(shù)的分旬權(quán)重，具體的分區(qū)分旬水分敏感系數(shù)值見表1。

表1 安徽省夏玉米分區(qū)分旬水分敏感系數(shù)

1.2.2 夏玉米澇漬指數(shù)Qw

在氣象行業(yè)標(biāo)準(zhǔn)《冬小麥、油菜澇漬等級》[22]中澇漬指數(shù)計算方法的基礎(chǔ)上,結(jié)合安徽省的實際情況,通過優(yōu)化和改進(jìn)建立夏玉米澇漬指數(shù)的計算公式。

(1)

其中,Qwi為第i時間段的澇漬指數(shù)，Ri、Rdi、Si分別為第i時段的降水量、降水日數(shù)和日照時數(shù)，Rmax、Rdmax、Smax分別為序列內(nèi)的降水量、降水日數(shù)和日照時數(shù)的最大值。

本文中計算的時間尺度以旬作為基本單位，分別計算全省各站點1971—2018年6月上旬至9月下旬逐旬的澇漬指數(shù)，例如計算合肥站2000年的8月上旬的澇漬指數(shù)，其中Ri為8月上旬的降水量，Rmax為1971—2018年共48年每年8月上旬降水量的最大值，降水日數(shù)和日照時數(shù)情況類似。參數(shù)a1、a2、a3為表征降水量、雨日和日照時數(shù)對形成澇漬災(zāi)害貢獻(xiàn)的經(jīng)驗系數(shù)，參考澇漬等級的氣象行業(yè)標(biāo)準(zhǔn)并結(jié)合本地區(qū)澇漬發(fā)生的實際情況，a1、a2、a3的值分別取1、1、0.75。

1.2.3 夏玉米趨勢產(chǎn)量的分解

影響作物最終產(chǎn)量形成的因素很多，一般把作物最終產(chǎn)量分解為趨勢產(chǎn)量和氣象產(chǎn)量，其中農(nóng)業(yè)技術(shù)措施類包括施肥、經(jīng)營管理、病蟲害控制、品種改良以及其它增產(chǎn)措施等，它反映了一定歷史時期的社會經(jīng)濟(jì)技術(shù)發(fā)展水平；相應(yīng)的產(chǎn)量分量稱為時間技術(shù)，簡稱趨勢產(chǎn)量。而由于年際間氣象條件的差異造成作物產(chǎn)量的波動，相應(yīng)的產(chǎn)量分量稱為氣象產(chǎn)量。本文采用直線滑動平均法模擬趨勢產(chǎn)量，該方法是將產(chǎn)量的時間序列在某個階段內(nèi)的變化看作線性函數(shù)，呈一直線；隨著階段的連續(xù)滑動，直線不斷改變位置、后延滑動而反映產(chǎn)量歷史演變趨勢的連續(xù)變化；依次求取各個階段內(nèi)的回歸模型，而各時間點上各直線滑動回歸模擬值的均值, 即為其趨勢產(chǎn)量值[27]。

y=yt+yw

(2)

式中，y為作物的實際產(chǎn)量，單位為kg/hm2；yt為作物的趨勢產(chǎn)量，單位為kg/hm2；yw為作物產(chǎn)量的氣象分量，又稱為氣象產(chǎn)量，單位為kg/hm2。

1.2.4 產(chǎn)量澇漬災(zāi)損風(fēng)險指數(shù)

災(zāi)損減產(chǎn)風(fēng)險指數(shù)[13]是指不同減產(chǎn)率與相應(yīng)出現(xiàn)概率乘積的總和，是反映減產(chǎn)風(fēng)險程度的一個指標(biāo)，指數(shù)越大，風(fēng)險就越大。玉米產(chǎn)量災(zāi)損風(fēng)險指數(shù)(M)由3方面因素綜合構(gòu)成，其中產(chǎn)量災(zāi)損風(fēng)險強(qiáng)度指數(shù)(Q)是災(zāi)損風(fēng)險在強(qiáng)度上的表示；減產(chǎn)變異系數(shù)(Cv)可反映實際災(zāi)損的年際波動情況，是災(zāi)損風(fēng)險振幅的體現(xiàn)；澇漬年發(fā)生概率(F)是澇漬在時域尺度上的度量；3者分別代表實際災(zāi)損風(fēng)險發(fā)生的強(qiáng)度、幅度和頻度，均與實際災(zāi)損呈正比。

玉米減產(chǎn)變異系數(shù)(Cv)指減產(chǎn)年(氣象產(chǎn)量為負(fù))玉米產(chǎn)量的波動情況。

(3)

采用概率分布函數(shù)計算玉米不同等級的減產(chǎn)強(qiáng)度及其發(fā)生概率，定義不同等級的減產(chǎn)強(qiáng)度與其發(fā)生概率之積的總和為災(zāi)損的風(fēng)險強(qiáng)度指數(shù)(Q)。

(4)

式中,Ei為第i個澇漬等級的減產(chǎn)強(qiáng)度，共3個等級的澇漬減產(chǎn)強(qiáng)度，減產(chǎn)5%～10%為輕度，減產(chǎn)10%～20%為中度，減產(chǎn)20%以上為重度；Pi為第i個減產(chǎn)強(qiáng)度出現(xiàn)的概率。

基于產(chǎn)量減產(chǎn)序列獲取因澇漬導(dǎo)致的減產(chǎn)年，并確定澇漬減產(chǎn)年發(fā)生概率(F):

(5)

式中,d為因澇漬導(dǎo)致的減產(chǎn)年數(shù)；D為統(tǒng)計的總年數(shù)，1971—2018年共48年。

則玉米災(zāi)損風(fēng)險指數(shù)M的表達(dá)式為

(6)

2 結(jié)果與分析

2.1 夏玉米澇漬指數(shù)的時間變化特征

安徽省1971—2018年夏玉米全生育期澇漬指數(shù)多年平均值為0.239，全省平均值隨年代變化呈顯著上升趨勢，其上升趨勢為0.027/10a，通過α=0.1的顯著性檢驗(圖2)。

圖2 1971—2018年安徽省夏玉米全生育期澇漬指數(shù)變化

全省分旬澇漬指數(shù)的時間變化特征不盡相同，其中6月上旬和8月上旬澇漬指數(shù)均呈明顯的上升趨勢，通過α=0.1的顯著性檢驗，其余各旬變化不顯著。其中澇漬指數(shù)最大的是6月下旬，其次是7月上旬和7月中旬，期間正值我省的梅雨期，降水較為集中，易發(fā)生澇漬害；澇漬指數(shù)最小的是9月中旬和9月下旬，期間受副熱帶高壓控制，正值我省晴熱少雨的時期，出現(xiàn)澇漬害的可能性也較低(表2)。

從夏玉米全生育期澇漬指數(shù)的年代際變化特征來看，20世紀(jì)70年代澇漬指數(shù)最小，20世紀(jì)80年代澇漬指明顯增大，20世紀(jì)90年代澇漬指數(shù)有所減小，進(jìn)入21世紀(jì)澇漬指數(shù)也是明顯增大，其中20世紀(jì)80年代和21世紀(jì)分別為玉米澇漬指數(shù)的2個高值時期，這與夏玉米典型澇漬災(zāi)害年多發(fā)生在這2個時期的實際情況相吻合。隨機(jī)選取了幾個代表旬分析了其澇漬指數(shù)年代的特征，存在的共同的特征為20世紀(jì)70年代普遍為澇漬指數(shù)的低值區(qū)(表3)。

表2 1971—2018年安徽省夏玉米分旬澇漬指數(shù)及其變化特征

表3 1971—2018年安徽省夏玉米旬及全生育期各年代澇漬指數(shù)

2.2 夏玉米澇漬指數(shù)的空間變化特征

圖3為安徽省1971—2018年各站點夏玉米全生育期澇漬指數(shù)的空間分布特征，全省各站點澇漬指數(shù)的變化幅度在0.08～0.38之間，其中大別山區(qū)、江南中東部和淮北局部是夏玉米澇漬指數(shù)的高值區(qū)(Qw>0.25)，淮北北部、江淮北部、沿江西部以及皖南南部是夏玉米澇漬指數(shù)的低值區(qū)(Qw<0.2)。

圖3 1971—2018年安徽省夏玉米全生育期 澇漬指數(shù)空間分布

安徽省1971—2018 年各站點夏玉米分月澇漬指數(shù)空間分布特征顯示：6月正值安徽省夏玉米的苗期生長階段，全省自南向北陸續(xù)進(jìn)入汛期，降水南多北少，澇漬指數(shù)呈明顯的緯向分布特征，沿淮淮北地區(qū)澇漬指數(shù)普遍在0.05～0.15，其它大部分地區(qū)0.2～0.5，其中江南部分地區(qū)在0.5以上(圖4a)；7月正值安徽省夏玉米的拔節(jié)抽雄期，也是夏玉米需水關(guān)鍵期，期間雨帶北移，北部的澇漬指數(shù)普遍高于南部，沿淮淮北大部澇漬指數(shù)普遍在0.3～0.4，其他大部分地區(qū)普遍在0.1～0.3(圖4b)；8月安徽省夏玉米普遍處于抽雄至開花吐絲期，澇漬指數(shù)普遍在0.1～0.4，南北差異不是很明顯，2個高值區(qū)分別位于大別山區(qū)和皖南中東部地區(qū)(圖4c)；9月安徽省夏玉米普遍處于灌漿至乳熟期，澇漬指數(shù)普遍在0.05～0.25，南部的澇漬指數(shù)普遍高于北部(圖4d)。

2.3 夏玉米澇漬災(zāi)害等級指標(biāo)確定

參考國家行業(yè)標(biāo)準(zhǔn)《冬小麥、油菜澇漬等級》[22]，結(jié)合因澇漬導(dǎo)致的夏玉米減產(chǎn)年，以各旬夏玉米澇漬指數(shù)(以兩旬的平均值作為判別的依據(jù))確定不同減產(chǎn)程度下的夏玉米澇漬等級指標(biāo)；其中分別確定了各生育階段的澇漬等級指標(biāo)，將夏玉米澇漬分為輕度、中度、重度3級(表4)。

為了確定夏玉米澇漬等級指標(biāo)在安徽省的適用性，對其進(jìn)行了指標(biāo)驗證，同時為了排除其它災(zāi)害導(dǎo)致夏玉米減產(chǎn)的樣本，通過統(tǒng)計年鑒并結(jié)合當(dāng)年災(zāi)害實際發(fā)生的情況，對災(zāi)害年樣本進(jìn)行篩選，挑選出安徽省1981-2010年主要由澇漬引起減產(chǎn)的年份作為夏玉米澇漬災(zāi)損的樣本，全省各縣共計467個樣本，對災(zāi)害評估指標(biāo)(表5)進(jìn)行檢驗。檢驗結(jié)果顯示，計算得到的澇漬等級對應(yīng)的減產(chǎn)率與實際由于澇漬災(zāi)害導(dǎo)致的減產(chǎn)率一致的比例為82.2%，其中中度澇漬樣本達(dá)到一致的比例最高，為85.5%，澇漬等級與實際減產(chǎn)率相差1個等級的比例為15.2%，兩者之和為97.4%，可以看出夏玉米澇漬災(zāi)害指標(biāo)的評價結(jié)果與實際澇漬減產(chǎn)的結(jié)果較為一致。

表4 安徽省夏玉米澇漬等級指標(biāo)

表5 1981—2010年安徽省夏玉米澇漬災(zāi)害評估指標(biāo)檢驗結(jié)果

2.4 安徽省夏玉米產(chǎn)量澇漬災(zāi)損風(fēng)險分布特征

圖5 安徽省夏玉米產(chǎn)量澇漬災(zāi)損風(fēng)險分布

基于產(chǎn)量災(zāi)損風(fēng)險理論，分別計算安徽省各個站點的夏玉米產(chǎn)量澇漬災(zāi)損風(fēng)險，采用ArcGIS中的克里格插值方法對各站的澇漬災(zāi)損風(fēng)險指數(shù)進(jìn)行插值，得到安徽省夏玉米產(chǎn)量澇漬災(zāi)損風(fēng)險分布(圖5)，可知安徽省夏玉米產(chǎn)量澇漬災(zāi)害損失的低風(fēng)險區(qū)主要分布在淮北大部以及江淮的東北部地區(qū)，其值普遍在0.05～0.16;中風(fēng)險區(qū)主要分布在沿淮中西部、江淮大部和沿江江南中西部地區(qū)，其值普遍在0.16～0.24;高風(fēng)險區(qū)主要分布大別山區(qū)、江淮東部和江南東部地區(qū)，其值普遍在0.24～0.32，夏玉米產(chǎn)量澇漬災(zāi)損風(fēng)險指數(shù)與夏玉米全生育期澇漬指數(shù)空間分布特征存在一定的差異，其中大別山區(qū)較為吻合，為澇漬指數(shù)的高值區(qū)和產(chǎn)量澇漬災(zāi)害損失風(fēng)險的高值區(qū)，而江南中東部作為夏玉米全生育期澇漬指數(shù)的高值區(qū)，實際上僅江南東部為夏玉米產(chǎn)量澇漬災(zāi)害損失的高風(fēng)險區(qū)，另外江淮東部為山地丘陵區(qū)，受地形等因素的影響，其也成為夏玉米產(chǎn)量澇漬災(zāi)害損失的高風(fēng)險區(qū)，因此澇漬災(zāi)害風(fēng)險的高低與當(dāng)?shù)睾０屋^高，土壤類型和當(dāng)?shù)氐膶嶋H作物抗災(zāi)能力存在較大的關(guān)系。

3 結(jié)論與討論

(1) 安徽省夏玉米澇漬的時空特征顯示：夏玉米全生育期澇漬指數(shù)以0.027/10a的速率呈顯著上升趨勢，其中6月上旬和8月上旬上升趨勢最為明顯，澇漬指數(shù)的年代際變化總體呈現(xiàn)上升的趨勢，僅90年代澇漬指數(shù)有所減小；澇漬指數(shù)的空間分布其值變幅在0.08～0.38之間，其中大別山區(qū)、江南中東部和淮北局部是夏玉米澇漬指數(shù)的高值區(qū)，分月的澇漬指數(shù)空間分布隨著夏季雨帶的南北擺動呈現(xiàn)出不同的特征。

(2) 根據(jù)國家氣象行業(yè)標(biāo)準(zhǔn)《冬小麥、油菜澇漬等級》，結(jié)合安徽夏玉米全生育期澇漬發(fā)生特征，確定了安徽省夏玉米澇漬災(zāi)害等級指標(biāo)；指標(biāo)驗證結(jié)果與實際因澇漬災(zāi)害而減產(chǎn)的結(jié)果較為一致，其中澇漬等級與實際減產(chǎn)率一致的比例為82.2%，一致和相差1個等級的比例大于95%。

(3) 根據(jù)夏玉米產(chǎn)量澇漬災(zāi)損綜合風(fēng)險指數(shù)將安徽省夏玉米劃分為高、中、低3 個風(fēng)險區(qū)，低風(fēng)險區(qū)主要分布在淮北大部以及江淮的東北部地區(qū)；中風(fēng)險區(qū)主要分布在沿淮中西部、江淮大部和沿江江南中西部地區(qū)，高風(fēng)險區(qū)主要分布大別山區(qū)、江淮東部和江南東部地區(qū)。

由于作物災(zāi)害風(fēng)險分析尚處在研究階段，在風(fēng)險評估指標(biāo)的選擇上，僅選擇了災(zāi)損風(fēng)險發(fā)生的強(qiáng)度、幅度和頻度等常用的指標(biāo)，有待作更深一步的研究確定更精細(xì)和更合理的風(fēng)險評估指標(biāo)；另外在確定漬指數(shù)等級分類指標(biāo)時，存在很大的經(jīng)驗性，還需更進(jìn)一步對指標(biāo)進(jìn)行驗證，精確分類指標(biāo)。