胡 珂，何 鵬，高文波

（1四川省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)信息與農(nóng)村經(jīng)濟(jì)研究所，成都 610066； 2成都理工大學(xué)旅游與城鄉(xiāng)規(guī)劃學(xué)院，成都 610059；3四川省儀隴中學(xué)校，四川儀隴 637600；4四川省農(nóng)業(yè)科學(xué)院大數(shù)據(jù)中心，成都 610066）

0 引言

長江流域是中國最重要的油菜產(chǎn)區(qū)，冬油菜常年種植面積和產(chǎn)量均占全國的80%以上[1]。然而受季風(fēng)氣候影響，中國長江流域油菜產(chǎn)區(qū)秋、冬生長季節(jié)階段性持續(xù)降水偏多，連陰雨出現(xiàn)概率高、持續(xù)時間長[2]，極易發(fā)生漬害，嚴(yán)重影響油菜的生長發(fā)育和產(chǎn)量品質(zhì)。漬害是指土壤含水量達(dá)到飽和或者過飽和，導(dǎo)致植株根系的生長環(huán)境缺少氧氣，影響植株正常的生理代謝，最終引起農(nóng)作物產(chǎn)量明顯下降甚至絕收的一種農(nóng)業(yè)氣象災(zāi)害[3]。油菜耐漬能力較差，萌發(fā)期和幼苗期漬水后，會降低出苗率和成苗質(zhì)量[4]，削弱油菜群體效應(yīng)[5]；蕾薹期遭受漬害后會抑制油菜一次分枝的形成，開花期、角果期的漬害將影響有效角果數(shù)、千粒質(zhì)量[6]，最終都將影響油菜正常的生長發(fā)育，導(dǎo)致減產(chǎn)。近年來，油菜漬害已逐漸引起國內(nèi)外學(xué)者的關(guān)注，并從油菜漬害的生理響應(yīng)[7-8]、漬害對油菜產(chǎn)量的影響[9-10]以及耐濕性品種篩選[11-12]等方面做了大量研究。作為一種區(qū)域性頻發(fā)的農(nóng)業(yè)氣象災(zāi)害，不少學(xué)者也開展了油菜濕漬害氣象指標(biāo)以及氣候變化對油菜漬害影響等方面的研究。霍治國等[13-14]選取了降水量、降水日數(shù)以及日照時數(shù)等氣象因子，構(gòu)建了油菜澇漬指數(shù)計算模型，形成了適用于中國南方地區(qū)油菜澇漬澇漬等級行業(yè)標(biāo)準(zhǔn)。在此基礎(chǔ)上，劉瑞娜等[15]利用逐旬氣象要素資料，參照行業(yè)標(biāo)準(zhǔn)中生育期的劃分，對安徽省油菜不同生育期澇漬災(zāi)害的特征和時空變化特征進(jìn)行了系統(tǒng)研究；鞠英芹等[16]構(gòu)建了油菜全生育期的綜合澇漬指數(shù)，利用了經(jīng)驗正交函數(shù)分解和小波分析，研究了江淮地區(qū)油菜漬害的時空分布特征；郭翔等[17]分析了四川省盆地區(qū)1981—2020 年4 個年代際油菜播種期澇漬指數(shù)、輕/中度澇漬害發(fā)生頻率、澇漬害危險性的時空變化特征。

受華西秋雨影響，位于長江上游的西南地區(qū)冬油菜種植地區(qū)苗期的漬害出現(xiàn)概率大，且造成的后果尤為嚴(yán)重[18]。在氣候變化背景下開展西南地區(qū)冬油菜關(guān)鍵生育期氣象災(zāi)害分布特征的研究，對于優(yōu)化區(qū)域冬油菜生產(chǎn)布局、提升防災(zāi)減災(zāi)、保障油料供給安全有重要意義。然而，要準(zhǔn)確分析跨地區(qū)的作物關(guān)鍵生育期災(zāi)害危險性的分布特征，必須考慮不同地理環(huán)境（經(jīng)緯度和海拔）條件下作物關(guān)鍵生育期起止時間的差異。迄今為止，結(jié)合區(qū)域油菜苗期分布特征，系統(tǒng)開展西南地區(qū)油菜苗期漬害危險性分布特征的研究鮮有報道。本研究以西南地區(qū)四川、重慶、貴州、云南四省（市）為研究區(qū)域，利用區(qū)域內(nèi)農(nóng)業(yè)氣象站點數(shù)據(jù)系統(tǒng)分析了研究區(qū)域內(nèi)甘藍(lán)型冬油菜出苗期的時空分布特征，結(jié)合行業(yè)標(biāo)準(zhǔn)QX/T107—2009《冬小麥、油菜澇漬等級》中油菜漬害指標(biāo)的定義與計算方法，開展了西南地區(qū)冬油菜苗期漬害危險性的時空分布特征研究，以期為西南地區(qū)油菜種植布局優(yōu)化和防災(zāi)減災(zāi)工作提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)域

本研究主要的研究對象為甘藍(lán)型冬油菜，為了保證數(shù)據(jù)的一致性，通過統(tǒng)計年鑒及文獻(xiàn)資料查閱，結(jié)合前人對中國油菜種植區(qū)劃的研究[19-20]，排除部分以春油菜為主的種植區(qū)域，最終篩選出有冬油菜種植的397個縣級行政單位所在范圍作為研究區(qū)域，即研究區(qū)域包括云南省、貴州省、重慶市和四川省除西北高原以外的地區(qū)，具體見圖1。

圖1 研究區(qū)域及站點分布圖

1.2 資料來源

氣象數(shù)據(jù)來自中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)，選取符合條件的87 個國家氣象站點1971—2020年9月—翌年3月的逐日數(shù)據(jù)，包括降雨量和日照時數(shù)等指標(biāo)。氣象站點篩選條件和處理方式如下：（1）剔除缺失率大于1%的站點，對缺測序列小于1%的站點利用多年均值替代法進(jìn)行插補(bǔ)；（2）考慮站點遷移造成的非均一性影響，排除遷移信息不明確站點，運(yùn)用多元回歸模型插值方法對遷移站點數(shù)據(jù)進(jìn)行校正[21]。冬油菜生育期資料來源于以及中國氣象數(shù)據(jù)網(wǎng)中國農(nóng)作物生長發(fā)育狀況資料數(shù)據(jù)集和前人研究論文[22-27]，整理選取西南地區(qū)四川、云南、貴州、重慶四個省市內(nèi)50個農(nóng)業(yè)氣象觀測點的油菜生育觀測資料，其中包括油菜發(fā)育期名稱、發(fā)育日期、發(fā)育程度、發(fā)育期距平等數(shù)據(jù)。地理空間數(shù)據(jù)來源于地理空間數(shù)據(jù)云(http://www.gscloud.cn/search)，獲取了研究區(qū)域90 m×90 m 數(shù)字高程模型數(shù)字高程模型。研究區(qū)域、氣象站點和農(nóng)業(yè)氣象站點分布和數(shù)字高程模型見圖1。

1.3 研究方法

1.3.1 漬害指數(shù)計算 根據(jù)國家規(guī)范的氣象行業(yè)標(biāo)準(zhǔn)QX/T107—2009《冬小麥、油菜澇漬等級》建立油菜苗期的澇漬指標(biāo)[14]，油菜的澇漬指標(biāo)公式為：

其中，Qw為澇漬指數(shù)；R為旬降水量(mm)；Rmax為近3個年代的旬最大降水量(mm)；DR為旬降水日數(shù)（天）；D為旬天數(shù)（天）；S為旬日照時數(shù)(h)；Smax為旬可能日照時數(shù)(h)；b1為降水量；b2為降水日數(shù)；b3為日照時數(shù)對澇漬影響所占權(quán)重，本研究結(jié)合以往研究[28]和專家意見，對3個系數(shù)分別取值為b1=0.75，b2=0.75，b3=0.6。

1.3.2 漬害等級指標(biāo) 按照行業(yè)標(biāo)準(zhǔn)QX/T107—2009《冬小麥、油菜澇漬等級》逐旬計算87個站點歷年苗期澇漬指數(shù)Qw，計算2旬澇漬指數(shù)均值Qw’，將油菜苗期澇漬害分輕、中、重3個等級（表1）[28]。

表1 油菜苗期澇漬等級指標(biāo)

1.3.3 油菜苗期漬害危險性評估模型的建立 參考前人相關(guān)研究[17,21]，致災(zāi)因子危險性同時取決于致災(zāi)因子強(qiáng)度和致災(zāi)因子發(fā)生的概率，如式(2)所示。

式中，R為致災(zāi)因子危險性，F(xiàn)表示致災(zāi)因子發(fā)生的強(qiáng)度，P表示致災(zāi)因子發(fā)生的概率，以某站點多年發(fā)生災(zāi)害的頻率表示。

研究區(qū)域苗期漬害的強(qiáng)度是受一次災(zāi)害發(fā)生的致災(zāi)因子范圍和程度共同作用的結(jié)果，如式(3)所示。

式中，F(xiàn)為冬油菜苗期漬害強(qiáng)度。A表示漬害影響范圍，不同的漬害等級影響范圍以研究時間段內(nèi)發(fā)生過輕、中、重澇漬災(zāi)害的站點數(shù)量占總站數(shù)的百分比表示。I為漬害影響程度，結(jié)合以往研究[29]和專家意見，以數(shù)值1、2、3表征輕、中、重3種等級漬害影響程度。

漬害發(fā)生概率可用站點發(fā)生漬害年份數(shù)占總樣本年份數(shù)量之比來表示，如式(4)所示。

式中，P為冬油菜苗期漬害發(fā)生概率。n表示冬油菜苗期發(fā)生過漬害（即澇漬指數(shù)Qw≥0.8）的年數(shù)，N表示總年數(shù)。

考慮到漬害指數(shù)Qw計算中旬降水量和旬降水日數(shù)的年際變化較大，本研究引入降水量變異系數(shù)CV來表征漬害致災(zāi)因子的變異程度。變異系數(shù)CV(Coefficient of Variation)是標(biāo)準(zhǔn)差與平均值之比，表示了樣本降水量的變化幅度，其值表征了實際數(shù)據(jù)偏離平均值的程度。計算公式如式(5)所示。

當(dāng)前國內(nèi)質(zhì)量控制工具的利用尚不成熟,仍處于學(xué)習(xí)階段。而在美國的筆譯項目中,術(shù)語管理貫穿整個項目,質(zhì)量控制更是一個動態(tài)過程,并緊密結(jié)合各類計算機(jī)軟件技術(shù)。

綜上，本研究構(gòu)建了冬油菜苗期澇漬災(zāi)害致災(zāi)因子危險性模型，如式(6)所示。

式中，Ri表示各個子區(qū)域冬油菜苗期漬害的危險性指數(shù)，F(xiàn)a表示輕度漬害的強(qiáng)度，F(xiàn)b表示中度漬害的強(qiáng)度，F(xiàn)c表示重度漬害的強(qiáng)度，Pa表示輕度漬害的頻率，Pb表示中度漬害的頻率，Pc表示重度漬害的頻率。

1.3.4 空間特征分析方法 對于不同的地理因子，為了保證高精度，采取的插值方法各有不同。本研究根據(jù)研究對象的不同特征，分別采用了多元回歸殘差反距離訂正法[30]和反距離權(quán)重插值法進(jìn)行空間分析與制圖表達(dá)。

2 結(jié)果與分析

2.1 冬油菜苗期空間分布特征

研究區(qū)域跨度范圍大，不同區(qū)域冬油菜苗期有差異。為準(zhǔn)確分析研究區(qū)域苗期漬害危險性特點，需對研究區(qū)域苗期分布特征進(jìn)行研究。本研究以研究區(qū)域內(nèi)農(nóng)業(yè)氣象站點穩(wěn)定出苗期的日序數(shù)(Day Of Year,DOY)為因變量y，以農(nóng)業(yè)氣象站點的緯度、經(jīng)度和海拔高程為自變量x1、x2、x3，分別考慮高程(x3)、緯度+高程(x1，x3)、緯度+經(jīng)度(x1，x2)、經(jīng)度+高程(x2，x3)、緯度+經(jīng)度+高程(x1，x2，x3)等5 種方式構(gòu)建出苗期日序數(shù)(y)的多元回歸方程，比較油菜出苗期日序數(shù)與不同方式組合的相關(guān)程度，結(jié)果見表2。

表2 多元回歸5種模型對比

由表2可知，高程x3、經(jīng)度x2+高程x3、緯度x1+經(jīng)度x2+高程x3與出苗期日序數(shù)構(gòu)建的多元回歸方程的P值均小于0.01，表明出苗期日序數(shù)與以上三種回歸方程模式呈極顯著相關(guān)的特征。其中，以（經(jīng)度x2+高程x3）為自變量與油菜出苗期日序數(shù)y構(gòu)建的多元回歸方程的決定系數(shù)R2、F值、P值最優(yōu)。因此，本研究以經(jīng)度為x2、海拔為x3，苗期開始的日期的日序數(shù)(DOY)為y值，常數(shù)項b，x2系數(shù)a2，x3系數(shù)a3構(gòu)建出苗期與經(jīng)度、海拔的多元回歸方程，如式(7)所示。

式(7)即為多元回歸殘差反距離訂正法中所應(yīng)用的多元回歸方程。

在ArcGIS 軟件中利用多元回歸殘差反距離訂正法分析了研究區(qū)域苗期的空間分布特征，得出的研究區(qū)域冬油菜分布特征圖結(jié)果見圖2a。

圖2 研究區(qū)域甘藍(lán)型冬油菜出苗期分區(qū)圖

由圖2a可知，研究區(qū)域東北部的大部分地區(qū)苗期開始的日序數(shù)在254～279 d 區(qū)間，西南部的地區(qū)在254～290 d區(qū)間，中部地區(qū)以291～317 d區(qū)間為主，西部地區(qū)苗期開始最晚，日序數(shù)在303～383 d。根據(jù)苗期開始的空間特征，結(jié)合研究區(qū)域甘藍(lán)型冬油菜苗期平均為93 d 這一特征將研究區(qū)域劃分為4 個子區(qū)域（見圖2b），分別是：Ⅰ區(qū)，是研究區(qū)域冬油菜最早進(jìn)入苗期的地區(qū)，大部分油菜種植區(qū)域9月中下旬進(jìn)入苗期，主要包括四川盆地和貴州東部，該區(qū)域雖然日照較少，但平均海拔較低，四面環(huán)山，寒潮不易入侵，熱量條件較好，降水較為適中。Ⅱ區(qū)，即大部分油菜種植區(qū)域10月上中旬進(jìn)入苗期，該區(qū)域主要是云南的南部，氣候干濕兩季分明，南部區(qū)域的降水、熱量條件都較為豐富。Ⅲ區(qū)，即大部分油菜種植區(qū)域10 月下旬進(jìn)入苗期，主要包括云南省的東北部、中部和貴州西部地區(qū)，該區(qū)域海拔較高，熱量和水分相較Ⅱ區(qū)要少。Ⅳ區(qū)，是研究區(qū)域內(nèi)最晚進(jìn)入苗期的油菜種植區(qū)域，即油菜種植區(qū)域11月上中旬進(jìn)入苗期，主要包括四川省西南山地和云南省北部高山區(qū)，區(qū)域內(nèi)平均海拔較高，熱量條件和降水條件較差。

2.2 冬油菜苗期漬害強(qiáng)度結(jié)果

圖3 1971—2020年油菜苗期漬害發(fā)生臺站數(shù)年際統(tǒng)計圖

利用式(3)計算研究區(qū)域不同等級漬害強(qiáng)度。根據(jù)油菜苗期漬害等級指標(biāo)（表1）計算結(jié)果，梳理出研究區(qū)域氣象站點中發(fā)生過輕度漬害、中度漬害、重度漬害的站點數(shù)量，將發(fā)生過漬害的站點數(shù)量與總站點數(shù)量的比值作為漬害在研究區(qū)域的影響范圍，以數(shù)值1、2、3表征輕、中、重3種等級的漬害的影響程度，綜合考慮范圍和程度計算油菜苗期漬害強(qiáng)度，結(jié)果如表3所示。

表3 漬害強(qiáng)度表

由表2可知，1971—2020年間，研究區(qū)域內(nèi)冬油菜苗期有69%的站點發(fā)生輕度漬害，有64%的站點發(fā)生過中度漬害，55%的站點發(fā)生過重度災(zāi)害。1971—2020年間，研究區(qū)域內(nèi)冬油菜苗期輕度、中度、重度漬害強(qiáng)度指數(shù)分別為0.69、1.28和1.65。

2.3 冬油菜苗期不同漬害等級發(fā)生頻率

根據(jù)式(3)對研究區(qū)域發(fā)生漬害的等級頻率進(jìn)行計算，并將計算結(jié)果進(jìn)行空間可視化，得到不同等級的漬害發(fā)生的頻率如圖4。

圖4 不同等級漬害發(fā)生概率分布圖

由圖4 可知，研究區(qū)域內(nèi)輕度、中度、重度漬害發(fā)生的頻率區(qū)間分別為[0，0.33]、[0，0.20]、[0，0.16]，表明區(qū)域內(nèi)發(fā)生輕度漬害的概率較高，而中度和重度漬害發(fā)生的概率相對較低。從空間分布來看，研究區(qū)域冬油菜苗期漬害主要發(fā)生在東部地區(qū)，四川西南地區(qū)和云南省不同等級的苗期漬害發(fā)生的頻率均極低。

從不同等級漬害發(fā)生頻率和分布特征來看，冬油菜苗期輕度漬害發(fā)生的高頻區(qū)集中在四川省東部、重慶市和貴州省，其中四川東南部地區(qū)、重慶西部和貴州北部的輕度漬害發(fā)生概率較高，多年最大頻率值在[0.11，0.33]之間。中度漬害和重度漬害發(fā)生的高頻區(qū)間面積相對輕度漬害面積均有大幅減小。其中，中度漬害高頻區(qū)集中在四川東部的瀘州市、重慶市南部的江津區(qū)、綦江區(qū)和貴州省的遵義、畢節(jié)和貴陽等地區(qū)；重度漬害發(fā)生的高頻區(qū)在四川東南部瀘州敘永，貴州省的畢節(jié)、織金安順等地。

2.4 冬油菜苗期漬害危險性分布特征

利用式(5)計算出研究區(qū)域各氣象站點1971—2020年甘藍(lán)型冬油菜苗期漬害的危險性指數(shù)值，并利用ArcGIS 軟件進(jìn)行空間可視化，根據(jù)取值范圍利用GIS中的自然斷點法對危險性指數(shù)劃分了澇漬危險性等級，結(jié)果如圖5、表4所示。

表4 油菜苗期漬害危險性分級及面積統(tǒng)計

由圖5可知，研究區(qū)域冬油菜苗期漬害呈現(xiàn)“東高西低”的分布特征。結(jié)合表4和圖5可知，研究區(qū)域甘藍(lán)型冬油菜苗期漬害危險性較低的區(qū)域分布在云南地區(qū)、四川西南地區(qū)和北部地區(qū)，占整個區(qū)域面積的近60%；甘藍(lán)型冬油菜苗期漬害危險性中等及以上區(qū)域集中在四川東部地區(qū)、重慶地區(qū)和貴州地區(qū)，占比40.4%，并呈現(xiàn)出明顯的圈層特點，即由苗期漬害高風(fēng)險地區(qū)向四周風(fēng)險逐漸降低，冬油菜苗期漬害危險性高的區(qū)域集中在四川省、重慶市、和貴州省結(jié)合部，占研究區(qū)域面積約8%，包括四川省的宜賓、瀘州等地，重慶市的江津區(qū)、綦江區(qū)、永川區(qū)、南川區(qū)等地和貴州省的遵義、畢節(jié)、貴陽等地。

3 結(jié)論

本研究根據(jù)西南地區(qū)農(nóng)業(yè)氣象站點多年生育期數(shù)據(jù)，利用GIS 空間分析技術(shù)開展了甘藍(lán)型冬油菜苗期空間分布特征研究，結(jié)果表明西南地區(qū)甘藍(lán)型冬油菜苗期分布存在明顯的區(qū)域特征，因此要準(zhǔn)確分析冬油菜苗期漬害，有必要結(jié)合不同區(qū)域的冬油菜苗期起止時間來計算區(qū)域內(nèi)站點漬害指數(shù)。研究結(jié)果還發(fā)現(xiàn)，西南地區(qū)甘藍(lán)型冬油菜苗期漬害呈現(xiàn)“東高西低”的分布特征，漬害危險性較高的區(qū)域主要集中在四川省中東部、重慶市和貴州省，尤其是四川省宜賓市、瀘州市，重慶市江津區(qū)、綦江區(qū)、永川區(qū)、南川區(qū)，貴州省遵義市、畢節(jié)市、貴陽市等地冬油菜苗期漬害危險性最高。研究結(jié)果與區(qū)域連陰雨分布特征情況基本吻合[31-34]，表明秋冬季節(jié)連陰雨過程是西南地區(qū)甘藍(lán)型冬油菜苗期漬害形成的最重要因素，這與相似研究的結(jié)果[17,35]以及實際種植經(jīng)驗相一致。因此建議在上述區(qū)域開展諸如冬油菜品種改良和布局優(yōu)化等措施來減小冬油菜苗期漬害帶來的可能影響。油菜苗期前的降水量以及土壤墑情等因素也會一定程度上也會影響苗期漬害形成和漬害嚴(yán)重程度。在今后的研究中，應(yīng)將結(jié)合苗期前的降水、土壤墑情特征和苗期降水量、日照時數(shù)等數(shù)據(jù)，對油菜苗期漬害危險性及其影響做進(jìn)一步研究。