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

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

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

關鍵詞：氣候風險度;區(qū)劃;小麥;達州

中圖分類號：S16文獻標識碼：ADOI：10.19754/j.nyyjs.20200815044

收稿日期：2020-07-02

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

引言

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

1材料與方法

1.1研究區(qū)概況

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

1.2數據來源

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

1.3研究方法

1.3.1氣候風險度指標的建立

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

式中，Rt、Rr、Rs、R分別表示同一時期內溫度風險度、降水風險度、日照風險度和氣候風險度;At、Ar、As、A分別表示同一時期內溫度適宜度、降水適宜度、日照適宜度和氣候適宜度的標準差;St、Sr、Ss、S分別表示某一時期內的平均溫度適宜度、降水適宜度、日照適宜度和氣候適宜度。

1.3.2氣候適宜度模型的建立

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

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

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

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

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

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

1.3.3傾向率

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

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

1.3.4風險度與產量的相關性檢驗

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

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

1.3.5區(qū)劃方法

氣候要素與當地的經度、緯度、海拔等地理因子的關系模型可表示為：

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

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

1.3.6統(tǒng)計分析

采用Excel和SPSS進行數據與圖形處理與分析。

2結果與分析

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