達州小麥氣候風(fēng)險評估與區(qū)劃

何希望 張利平 鄭文全 覃紹彬 伍文玲

摘要：在小麥氣候適宜度基礎(chǔ)上，構(gòu)建風(fēng)險度模型，對達州小麥的氣候風(fēng)險度進行評估。結(jié)果表明，近幾十年來，達州氣溫和降水以及小麥氣候風(fēng)險度變化不大，趨勢較為平穩(wěn);日照風(fēng)險度偏高，變幅較大且有增加趨勢;小麥全生育期的日照風(fēng)險度最高，但對產(chǎn)量影響不顯著，降水風(fēng)險度次之，但對產(chǎn)量影響明顯，氣溫風(fēng)險度維持在較低水平且有下降趨勢。因此，在生產(chǎn)中可通過合理安排種植制度和種植區(qū)域，選擇耐濕且對日照不敏感的品種、采取合適的農(nóng)技措施，調(diào)控水分、改善田間通光條件，從而規(guī)避小麥種植氣候風(fēng)險。

關(guān)鍵詞：氣候風(fēng)險度;區(qū)劃;小麥;達州

中圖分類號：S16文獻標(biāo)識碼：ADOI：10.19754/j.nyyjs.20200815044

收稿日期：2020-07-02

作者簡介：何希望（1987-），男，本科，助理工程師。研究方向：農(nóng)業(yè)氣象和公共氣象服務(wù)。

引言

針對氣候因子對小麥生長發(fā)育的影響國內(nèi)外已經(jīng)做過大量研究工作，這些研究[3-5]主要側(cè)重于作物產(chǎn)量的統(tǒng)計分析、作物生理生態(tài)模擬及作物響應(yīng)全球氣候變化的試驗研究等方面。根據(jù)作物生理氣候需求曲線和研究區(qū)域的氣候特點，確定小麥各生育期參數(shù)值，應(yīng)用生態(tài)適宜性理論，構(gòu)建風(fēng)險度模型，研究歷年平均氣候條件下作物生長的風(fēng)險度和氣候變化對小麥農(nóng)業(yè)生態(tài)的影響，并對研究區(qū)域進行風(fēng)險度區(qū)劃，以期為本地區(qū)農(nóng)業(yè)的種植制度調(diào)控提供重要的參考。

1材料與方法

1.1研究區(qū)概況

研究區(qū)位于四川省東北部，屬中北亞熱帶濕潤季風(fēng)氣候，主要氣候特征是四季分明、雨熱同季，山地氣候特征明顯。2015年度達州小麥播種面積6.53萬hm2，總產(chǎn)19.90萬t，平均單產(chǎn)0.203t·hm-2，各地播種面占比見圖1。

1.2數(shù)據(jù)來源

小麥生育期資料來自宣漢國家農(nóng)業(yè)氣象觀測站1990—2016年的觀測數(shù)據(jù)（訂正）;各站點1961—2016年55a逐日溫度、降水量和日照時數(shù)等資料來源于四川省氣象局?jǐn)?shù)據(jù)中心，小麥播種面積和產(chǎn)量資料來源于達州市統(tǒng)計局。1∶5萬GIS資料來源于四川省氣候中心。

1.3研究方法

1.3.1氣候風(fēng)險度指標(biāo)的建立

本研究基于氣候適宜度模型構(gòu)建氣候風(fēng)險度指標(biāo)。綜合考慮溫度、降水量、日照時數(shù)、氣候條件適宜度水平的高低及達州小麥生長發(fā)育的實際情況，建立小麥氣候風(fēng)險度指標(biāo)，其計算分別見公式（1）～（4）：

式中，Rt、Rr、Rs、R分別表示同一時期內(nèi)溫度風(fēng)險度、降水風(fēng)險度、日照風(fēng)險度和氣候風(fēng)險度;At、Ar、As、A分別表示同一時期內(nèi)溫度適宜度、降水適宜度、日照適宜度和氣候適宜度的標(biāo)準(zhǔn)差;St、Sr、Ss、S分別表示某一時期內(nèi)的平均溫度適宜度、降水適宜度、日照適宜度和氣候適宜度。

1.3.2氣候適宜度模型的建立

根據(jù)黃璜等[6，7]的研究，結(jié)合達州實際情況，以旬為時間單位，建立溫度、降水和日照時數(shù)的適宜度模型，分別見公式（5）～（7）：

式中，St、Sr、Ss分別為小麥生育期間溫度、降水、日照時數(shù)的適宜度;T、R、Sun為溫度、降水量和日照時數(shù)的觀測值;Tmin、Tmax、T0分別為小麥在該時段內(nèi)的下限平均溫度、上限平均溫度和最適宜平均溫度。

R0為小麥旬需水量。一般而言，達州小麥濕害偏重，特別是播種-出苗期和成熟期，存在“兩頭濕、中間干”的氣候特點十分明顯。因此，認為此2個時期的旬降水少于適宜旬平均降水量時，降水處于適宜狀態(tài)，取值1。

Sun0為可照時數(shù)（單位為日）的50%，b為經(jīng)驗常數(shù)（b=0.538Sun+0.013）。受山地自然遮擋和冬半年云霧天氣較多等因素影響，達州日照條件較差，年日照百分率只有30%左右，且冬半年不足20%。按原日照模型參數(shù)進行計算，指數(shù)結(jié)果出現(xiàn)極化，因此，將50%的日照百分率作為臨界點，認為當(dāng)日照百分率達到50%以上時，達州小麥對光照條件的反應(yīng)達到適宜狀態(tài)，取值1。

各生育期權(quán)重采用專家打分法與層次分析法綜合確定，各旬及全生育期氣候適宜度的采用幾何平均法進行計算，見公式（8）：

式中，S為某發(fā)育期光溫水綜合適宜度。

1.3.3傾向率

傾向率表示氣候因子的變化趨勢，以下式計算：

式中，Xi是某一氣候因子;ti是時間;n是樣本數(shù)量;a是常數(shù);b是一元回歸系數(shù)。

1.3.4風(fēng)險度與產(chǎn)量的相關(guān)性檢驗

將小麥產(chǎn)量分解為趨勢產(chǎn)量和氣象產(chǎn)量，誤差一般忽略不計，然后計算相對氣象產(chǎn)量，用最小二乘法計算相對氣象產(chǎn)量和風(fēng)險度的相關(guān)系數(shù)，進行兩者間的相關(guān)性檢驗。

式中，y、yt、yq、ys、ε分別表示小麥的實際產(chǎn)量、趨勢產(chǎn)量、氣象產(chǎn)量、相對產(chǎn)量和誤差，趨勢產(chǎn)量用5a滑動平均法計算。

1.3.5區(qū)劃方法

氣候要素與當(dāng)?shù)氐慕?jīng)度、緯度、海拔等地理因子的關(guān)系模型可表示為：

式中，Y為生態(tài)氣候要素，φ、λ和h分別代表緯度、經(jīng)度和海拔高度等地理因子，函數(shù)f（φ，λ，h）為區(qū)劃指標(biāo)的方程表達式，ε為殘差，即：

利用達州、巴中和廣元3市15個氣象站的氣候觀測資料（因變量）與對應(yīng)的經(jīng)度、緯度和海拔高度等地理信息數(shù)據(jù)（自變量），計算各站的氣候風(fēng)險度并建立空間線性模型（參數(shù)見表1），運用Arcgis10進行空間分析，殘差ε采用反距離插值法進行空間內(nèi)插。

1.3.6統(tǒng)計分析

采用Excel和SPSS進行數(shù)據(jù)與圖形處理與分析。

2結(jié)果與分析

本研究用5a滑動平均法進行趨勢分析。由圖2可知，近50a來達州各地小麥氣候風(fēng)險度變化較大，受各地小氣候變化差異影響，北部明顯高于南部，但變化趨勢基本一致。本文以種植面積最大的渠縣為例，對達州小麥氣候風(fēng)險度進行描述。