張祖蓮，王命全，張山清，姚艷麗*，木沙江·艾代吐力，吐瑪熱色·吾拉木

（1.新疆大學(xué)，新疆 烏魯木齊 830017；2.新疆興農(nóng)網(wǎng)信息中心/新疆農(nóng)業(yè)氣象臺，新疆 烏魯木齊 830002；3.新疆教育管理信息中心，新疆 烏魯木齊 830049）

霜凍是指作物生長季土壤表面和植物表面溫度下降到0 ℃及以下，足以使植物遭受傷害甚至死亡的氣象災(zāi)害。春季終霜凍出現(xiàn)愈遲，對農(nóng)業(yè)的危害就愈嚴重[1]。霜凍是北疆主要的農(nóng)業(yè)氣象災(zāi)害，其發(fā)生頻次多、影響范圍大，對農(nóng)業(yè)危害嚴重[2-3]，常造成嚴重的經(jīng)濟損失。國內(nèi)眾多學(xué)者對霜凍給農(nóng)業(yè)帶來的災(zāi)害進行研究[4-6]。北疆地區(qū)終霜凍對棉花苗期的危害要大于南疆，1989—2018 年新疆生產(chǎn)建設(shè)兵團八師墾區(qū)棉花播種期70%集中在4 月中旬，終霜風(fēng)險達36.3%～41.4%，即使將播種期推至4 月26 日后，仍存在3.4%～8.6%的終霜凍風(fēng)險[7-8]。終霜凍已成為制約北疆棉花穩(wěn)產(chǎn)高產(chǎn)的主要限制性因素之一[9-12]。以上研究的基礎(chǔ)氣象數(shù)據(jù)均使用的是國家氣象站數(shù)據(jù)。

新疆地域遼闊、地形地貌復(fù)雜，氣候類型多樣，已有的國家氣象站均相隔幾百公里及以上，且大部分布在城市區(qū)，不能代表農(nóng)區(qū)氣候要素，尤其不能滿足分區(qū)域、分作物的農(nóng)業(yè)氣象精細服務(wù)需求。例如，2021 年4 月22—26 日，受西西伯利亞南下強冷空氣影響，全疆大部出現(xiàn)了一場雨雪、大風(fēng)及強降溫天氣過程，多地達到寒潮標準。4 月25—26 日凌晨，北疆大部最低氣溫降至0 ℃以下，對早播出苗的棉花、玉米及部分林果開花授粉造成較大影響。由于北疆棉區(qū)持續(xù)低溫，部分棉區(qū)棉花遭受低溫冷害危害。

為準確了解此次霜凍、冷害影響范圍及危害程度，本文選用中國氣象局陸面數(shù)據(jù)同化系統(tǒng)實時產(chǎn)品數(shù)據(jù)集（Chinese Land Data Assimilation System，CLDAS），該數(shù)據(jù)集為逐小時數(shù)據(jù)且空間分辨率為5 km×5 km，產(chǎn)品優(yōu)勢在于結(jié)合地面、衛(wèi)星等多源資料融合[13]。國內(nèi)很多學(xué)者對CLDAS 實況數(shù)據(jù)涉及新疆區(qū)域的氣溫進行過評估檢驗[14]，劉瑩等[15]利用全國地面自動氣象站對比分析2017 年10 月、2018 年1、4、7 月CLDAS-V2.0 氣溫，結(jié)果顯示：CLDAS 氣溫數(shù)據(jù)在中國區(qū)域整體上具有較高的可信度，其中，西北地區(qū)非獨立性檢驗為0.972，獨立性檢驗為0.950，平均偏差為-0.271 ℃，均方根誤差為2.406 ℃，平均絕對誤差為1.588 ℃，存在極端氣候事件發(fā)生強度的刻畫與站點觀測存在一定差異。本文檢驗CLDAS 數(shù)據(jù)在北疆4 月極端災(zāi)害天氣中是否滿足霜凍災(zāi)害監(jiān)測，利用北疆地面自動氣象站和CLDAS溫度數(shù)據(jù)對災(zāi)害天氣中溫度均方根誤差、氣溫預(yù)報準確率、最低氣溫、平均氣溫、降溫幅度及棉花受災(zāi)各溫度指標等分別對比檢驗分析。并通過對生產(chǎn)建設(shè)兵團六師、七師和八師的棉花、林果等進行受災(zāi)情況實地檢驗，主要驗證在沒有地面自動觀測站的棉區(qū)且CLDAS 數(shù)據(jù)顯示該區(qū)域達到氣象災(zāi)害指標實際情況是否受災(zāi)，已受災(zāi)棉區(qū)中CLDAS 溫度是否達到溫度氣象災(zāi)害指標。結(jié)合實地調(diào)研檢驗的結(jié)果，驗證CLDAS 多源數(shù)據(jù)產(chǎn)品與國家氣象站實際觀測數(shù)據(jù)的吻合度，棉花氣象災(zāi)害指標是否完善，為精細農(nóng)業(yè)氣象服務(wù)、挖掘高分辨率多源數(shù)據(jù)的應(yīng)用提供依據(jù)。

1 資料和方法

1.1 資料

研究區(qū)域位于新疆北部（圖1），42.30°～49.15°N，79.90°～96.05°E，區(qū)域內(nèi)共有1 120 個氣象站點，其中國家觀測氣象站53 個、自動氣象站1 067 個。研究時間為2021 年4 月20—27 日。

圖1 北疆地理位置與地形高程

研究資料選用中國氣象局陸面數(shù)據(jù)同化系統(tǒng)（CLDAS-V2.0）小時實況數(shù)據(jù)產(chǎn)品（0.05°×0.05°，GRIB 格式）。研究區(qū)域按5 km 格點共計20 447 個格點，其中發(fā)生霜凍的格點為18 534 個。

1.2 方法

1.2.1 插值方法

對于CLDAS 的5 km 格點的氣溫，本文采用雙線性插值獲得相應(yīng)位置的氣象要素值。雙線性插值，其核心思想是獲取最鄰近4 個格點的位置，在兩個方向分別進行一次線性插值，按照其距目標點的距離賦予相應(yīng)的權(quán)重進行線性插值。要獲得未知函數(shù)f 在點P =（x，y）的值，最鄰近4 個格點分別為Q11（x1，y1）、Q12（x1，y2）、Q21（x2，y1）以及Q22（x2，y2）。其公式如下：

式中：（x，y）表示經(jīng)緯度位置。若x=x1，y=y1則f（x，y）=f（Q11）；若x=x1，y=y2則f（x，y）=f（Q12）；若x=x2，y=y1則f（x，y）=f（Q21）；若x=x2，y=y2則f（x，y）=f（Q22）。

1.2.2 計算方法

（1）最大降溫幅度（Tmd）：降溫幅度越大，對出苗期棉花致災(zāi)越嚴重。此次天氣過程中逐站點或格點的最高溫度與最低溫度之差作為該站點或格點的降溫幅度，且最大溫度時刻早于最小溫度時刻。

（2）溫度≤0 ℃及累計小時數(shù)：終霜凍是決定棉花播種期的主要因子，早播出苗后往往遭霜凍襲擊，使棉苗嚴重受凍甚至凍死。低溫持續(xù)時間越長凍害越嚴重。統(tǒng)計每個站點或格點溫度≤0 ℃所有時刻的平均溫度及累計小時總數(shù)。

（3）溫度≤12 ℃及累計小時數(shù)：棉花在播種后至出苗前遇低溫陰雨天氣發(fā)生爛種爛芽的現(xiàn)象稱為苗前冷害。根據(jù)研究，隨著日平均氣溫≤12 ℃低溫日數(shù)和強度的增加，棉花苗期冷害加重。播種至出苗期間，從棉花開始發(fā)芽到胚軸顯著伸長，是棉花苗期冷害敏感期。統(tǒng)計北疆此次天氣中是否達到棉花苗期冷害氣象指標，并逐時刻統(tǒng)計每一個站點或格點溫度≤12 ℃的所有時刻的平均溫度及累計小時數(shù)。

在霜凍災(zāi)害分析中分別計算站點和格點的霜凍、棉花播種及出苗期間的各溫度指標對應(yīng)指標值，并在區(qū)域中對比檢驗分析。

1.2.3 檢驗方法

采用均方根誤差RMSE、氣溫準確率TT（溫度絕對誤差≤1 ℃，記為TT1℃；溫度絕對誤差≤2 ℃記為TT2℃）。其中，均方根誤差：

式中：ft為從多源數(shù)據(jù)插值到氣象站點獲取的氣象要素值，at為氣象站點值，N 為滑動時刻數(shù)或總站數(shù)。

氣溫準確率TT公式：

式中：NA為溫度絕對誤差≤1 ℃或2 ℃的格點（次）數(shù)，NB為溫度絕對誤差＞1 ℃或2 ℃的格點（次）數(shù)。

對CLDAS 溫度數(shù)據(jù)進行RMSE、TT1℃、TT2℃檢驗時，首先將CLDAS 的5 km×5 km 的雙線性插值到所有地面氣象站點位置，與氣象站點的值分別進行對比檢驗。

2 結(jié)果與分析

2.1 CLDAS 數(shù)據(jù)分時效對比檢驗分析

由圖2a 可知，4 月20 日00：00—27 日23：00，CLDAS 溫度數(shù)據(jù)檢驗的RMSE為0.77～3.92 ℃，其平均值為1.19 ℃；TT1℃為23.69%～90.9%，平均值為80.48%；TT2℃為43.82%～98.08%，平均值為93.95%，其中RMSE＞2 ℃有11 個時刻，分布在20 日01：00、22日22：00、23 日23：00、24 日00：00—03：00、25 日13：00—14：00、26 日19：00、27 日08：00。剩下180 個時刻，RMSE、TT1℃、TT2℃平均值分別為1.06 ℃、83.12%、96.09%。

圖2 2021 年4 月20—27 日逐小時的CLDAS 數(shù)據(jù)均方根誤差、溫度準確率（a）和觀測站點平均溫度（b）的逐日逐時刻變化

由圖2b 可知，-3.2 ℃≤平均氣溫Tavg≤-2 ℃分布于24 日05：00—08：00，25 日06：00—07：00。平均氣溫＜0 ℃在22 日06：00—09：00，23 日04：00—10：00，24 日00：00—10：00，25 日01：00—08：00。

當(dāng)氣象站點降溫或升溫速度（斜率絕對值）越大，格點的均方根誤差越大，相應(yīng)溫度準確率也最低。

2.2 CLDAS 數(shù)據(jù)區(qū)域?qū)Ρ葯z驗分析

由圖3 可知，RMSE、TT1℃、TT2℃平均值分別為1.04 ℃、80.98%、94.01%，其中RMSE＞2 ℃平均海拔高度為1 575.83 m，共有62 個站點，阿勒泰地區(qū)7 個，博爾塔拉蒙古自治州4 個，昌吉回族自治州20 個，哈密市巴里坤3 個，75 團1 個，塔城地區(qū)6 個，烏魯木齊市4 個，伊犁州16 個，五家渠市1 個。RMSE≤2 ℃共有981 個站點，平均海拔高度為966.81 m，RMSE、TT1℃、TT2℃平均值分別為0.86 ℃、84.21%、96.33%。

圖3 檢驗CLDAS 數(shù)據(jù)的均方根誤差（a）和1 ℃、2 ℃氣溫準確率（b、c）區(qū)域分布

在時效檢驗中94.24%的時刻RMSE、TT1℃、TT2℃平均值分別為1.06 ℃、83.12%、96.09%。區(qū)域檢驗中94.15%的站點區(qū)域RMSE、TT1℃、TT2℃平均值分別為0.86 ℃、84.21%、96.33%。經(jīng)檢驗，CLDAS 格點實況數(shù)據(jù)與氣象站數(shù)據(jù)吻合度較高。

2.3 霜凍災(zāi)害天氣分析

2.3.1 平均氣溫（Average Temperature，Tavg）和最低氣溫（Minimum Temperature，Tmin）

由圖4a、4b 可知，此次霜凍災(zāi)害天氣過程中格點平均溫度分布和站點平均溫度等級分布基本一致。CLDAS 溫度數(shù)據(jù)中平均氣溫≤0 ℃的區(qū)域主要分布在海拔1 200 m 以上的阿勒泰北部及和布克賽爾區(qū)域，1 800 m 及以上的巴里坤山周邊，海拔2 000 m 以上的巴爾魯克山區(qū)域、博格達峰周邊、依連哈比爾尕山、博羅科努山、哈爾克他烏山。

圖4 站點和格點平均氣溫（a～b）、最低氣溫（c～d）、最大降溫幅度（e～f）分布

由圖4c、4d 可知，格點最低溫度分布和站點最低溫度等級分布基本一致。CLDAS 溫度數(shù)據(jù)的最低氣溫≤-12 ℃區(qū)域主要集中在海拔2 000 m 以上的中高山帶，共有4 656 個格點；-12 ℃≤最低氣溫≤-8 ℃集中在低山帶，共有1 925 個格點；最低氣溫＞0 ℃共有1 913 個格點，主要分布在哈密北部三塘湖、淖毛湖一帶以及巴里坤盆地，在準噶爾盆地腹地等局部區(qū)域也有零星分布。

2.3.2 降溫幅度分析

由圖4e、4f 可知，此次霜凍災(zāi)害天氣過程格點最大降溫幅度分布和站點的分布等級略有差別。CLDAS 溫度數(shù)據(jù)降溫幅度≤25 ℃格點區(qū)域少于站點分布區(qū)域，阿勒泰區(qū)域大部、伊犁地區(qū)鞏留和新源以南、巴里坤盆地降溫幅度為12～20 ℃；降溫幅度為0～10 ℃的情況較少，僅在阿勒泰地區(qū)以及天山山區(qū)有零星出現(xiàn)。以上有差別區(qū)域的格點降溫區(qū)域等級均大于站點降溫1 個等級。

最大降溫幅度＞30 ℃共有53 個氣象站。其中，伊犁哈薩克自治州27 個，博爾塔拉蒙古自治州2個，昌吉回族自治州6 個，哈密北部1 個，克拉瑪依市2 個，十三師1 個、四師6 個，塔城地區(qū)3 個，烏魯木齊3 個，五家渠市1 個，石河子市1 個。最大降溫幅度≤20 ℃共有136 個站，其中阿勒泰地區(qū)有91個，其余為零星分布。

由圖4e、4f 可知，降溫幅度最少的區(qū)域位于阿勒泰地區(qū)，降溫幅度≤25 ℃；整個天山山脈降溫幅度＞20 ℃，天山高海拔區(qū)域降溫在20～25 ℃，天山北半坡和伊犁地區(qū)降溫幅度均＞25 ℃且局部區(qū)域＞30 ℃。

2.3.3 溫度≤0 ℃及累計小時數(shù)分析

由圖5a、5b 可知，此次天氣過程中站點和格點T≤0 ℃所有時刻平均溫度分布基本一致，有差別的區(qū)域為：格點在沙灣、呼圖壁、精河、新源中部區(qū)域沒有顯示出氣溫≤0 ℃的情況，但這些區(qū)域的大部分氣象站點溫度在-1～0 ℃，且持續(xù)1～5 h。由氣象站數(shù)據(jù)（圖5c、5d）可知，1 h≤累計小時數(shù)≤5 h 主要集中在烏蘇、沙灣、呼圖壁、鞏留、炮臺、莫所灣站附近的站點。5 h＜累計小時數(shù)≤12 h 主要集中在阜康、奇臺、烏蘇、伊寧等區(qū)域。累計小時數(shù)＞72 h 主要集中在沿天山山脈海拔2 000 m 及以上、托里南部和布克賽爾海拔1 800 m 以上、哈巴河北部喀拉斯湖海拔1 500 m 及以上的區(qū)域，這些區(qū)域農(nóng)牧業(yè)較少。逐格點統(tǒng)計中，≤0 ℃的格點數(shù)共18 534 個，其中，霜凍小時數(shù)≤5 h 的格點共有2 119 個，小時數(shù)≤12 h的格點共有4 283 個，小時數(shù)≤24 h 的格點共有6 086 個，小時數(shù)≤48 h 的格點共有8 813 個，小時數(shù)≤72 h 的格點共有11 573 個。主要棉區(qū)在此次降溫至溫度≤0 ℃的持續(xù)時間均≤48 h。

圖5 站點（Stations）和CLDAS 格點（Grids）≤0 ℃的平均溫度（a、b）、≤0 ℃的累計小時數(shù)（c、d）、≤12 ℃的平均溫度（e、f）、≤12 ℃的累計小時數(shù)（g、h）分布

2.3.4 溫度≤12 ℃及累計小時數(shù)分析

由圖5e、5f 可知，此次霜凍災(zāi)害天氣過程中，所有氣象站點均有溫度≤12 ℃的情況。其中T≤12 ℃的所有時刻且平均溫度＜0 ℃的站點共有141 個。由圖5g、5h 可知，累計小時數(shù)為72～120 h 的T≤12 ℃的站數(shù)共296 個（昌吉回族自治州77 個，伊犁哈薩克自治州75 個，塔城地區(qū)23 個，哈密市北站21 個，石河子19 個，烏魯木齊市18 個，博爾塔拉蒙古自治州12 個，克拉瑪依市13 個，阿勒泰地區(qū)9 個，四師9 個，八師1 個，六師1 個，七師6 個，十三師1 個，五師4 個，五家渠市5 個）。由圖5g 可知，T≤12 ℃累計小時數(shù)站數(shù)＞120 h 的，阿勒泰地區(qū)、塔城地區(qū)、伊犁哈薩克自治州、烏魯木齊市明顯較多。由圖5 h可知，CLDAS 數(shù)據(jù)有差別的區(qū)域是站點在累計小時數(shù)為72～120 h 的中天山北半坡格點分布上沒有顯示完全。

此次天氣過程CLDAS 數(shù)據(jù)中所有格點均有T≤12 ℃的站點，其中≤12 ℃格點累計小時數(shù)＜72 h的站點共有671 個，其相應(yīng)的平均溫度為9.76 ℃。T≤12 ℃且累計小時數(shù)≥72 h 的格點占總格點數(shù)的96.72%。72～120 h 的格點數(shù)共有1 052 個，平均溫度為7.13 ℃。T≤12 ℃且累計小時數(shù)在120～150 h 的站點共有6 215 個，平均溫度為5.40 ℃。

對比站點和格點區(qū)域分布圖發(fā)現(xiàn)，T≤12 ℃的站點平均溫度分布和格點平均溫度分布等級基本一致。在累計小時數(shù)上，格點中有96.72%的格點溫度≤12 ℃，持續(xù)時間≥72 h。

2.4 應(yīng)用

據(jù)調(diào)查，受2021 年4 月22—26 日低溫雨雪天氣影響，北疆處于萌芽、出苗期的棉田均遭受不同程度的霜凍害、低溫冷害影響。但烏蘇市、瑪納斯縣、呼圖壁縣大部棉區(qū)由于棉花播種后未滴水或滴水晚（19 日后滴水），棉花尚未出苗，損失較輕。在沙灣、呼圖壁局部區(qū)域出現(xiàn)弱霜凍現(xiàn)象，CLDAS 溫度數(shù)據(jù)沒有≤0 ℃，此區(qū)域CLDAS 氣溫為0～1 ℃。

4 月25—26 日北疆棉區(qū)日最低氣溫降至-6.2～0.4 ℃，0 ℃以下低溫持續(xù)時間為1～9 h，其中七師124 團下雙河片區(qū)日最低氣溫為-6.2 ℃，0 ℃以下低溫持續(xù)9 h。對于剛出苗的棉花，因其生長勢弱、抗凍能力差，當(dāng)氣溫＜2 ℃（地表溫度≤0 ℃）時棉花子葉就會受凍。氣溫＜0 ℃，會導(dǎo)致棉花子葉、生長點嚴重受凍，致使棉苗死亡，并且隨低溫強度的增強和持續(xù)時間的增長死苗率增加。多源數(shù)據(jù)結(jié)果與調(diào)查實況吻合，多源融合資料（圖5e～5h）顯示，北疆棉區(qū)22—26 日，連續(xù)5～6 d 日平均氣溫＜12 ℃，其中博州棉區(qū)持續(xù)6 d，北疆其余棉區(qū)持續(xù)5 d，持續(xù)低溫雨雪天氣使棉花遭受不同程度的苗前或苗期低溫冷害。

實地調(diào)查發(fā)現(xiàn)：（1）棉花苗期遭受霜凍和低溫冷害的程度、范圍，除了受人工管理如滴水時間、封土等影響外，和當(dāng)?shù)貧鉁赜兄苯拥年P(guān)系。（2）棉花播種苗期冷害指標，隨著“干播濕出”和種子“精選包衣”技術(shù)的廣泛應(yīng)用，現(xiàn)在棉花種子萌發(fā)出苗期抗低溫的能力明顯增強，棉花播種苗期冷害的氣象指標，需要修訂完善。（3）由于大多數(shù)棉田沒有自動氣象站，不能及時準確獲取棉田氣溫變化，農(nóng)戶特別需要格點數(shù)據(jù)進行預(yù)報和監(jiān)測。多源精細化的氣溫監(jiān)測數(shù)據(jù)可以第一時間關(guān)注到受風(fēng)險地塊，在降溫持續(xù)時間還不長的情況下，及時采取補救措施，對挽救棉苗受災(zāi)至關(guān)重要。

3 結(jié)論

選擇CLDAS 實況格點數(shù)據(jù)對2021 年4 月春季北疆霜凍災(zāi)害性天氣進行分析，利用氣象站點數(shù)據(jù)對CLDAS 進行多要素檢驗。對北疆棉區(qū)進行實地災(zāi)害調(diào)查來檢驗多源融合數(shù)據(jù)的適用性，得出如下結(jié)論：

（1）北疆地形復(fù)雜，通過對CLDAS 的RMSE、TT1、TT2檢驗分析，CLDAS 格點實況數(shù)據(jù)與地面氣象站數(shù)據(jù)吻合度較高。

（2）格點數(shù)據(jù)中平均溫度、最低溫分布和站點平均溫度、最低溫度等級分布基本一致。在天山北波經(jīng)濟帶降溫幅度為20～25 ℃的局部區(qū)域CLDAS 數(shù)據(jù)為25～30 ℃；整個研究區(qū)域中共有14 個站點的降溫＜10 ℃在CLDAS 數(shù)據(jù)中沒有顯示出來。

（3）溫度≤0 ℃及累計小時數(shù)分析中，沙灣、呼圖壁、精河、新源中局部區(qū)域降溫至-1～0 ℃，區(qū)域內(nèi)CLDAS 溫度數(shù)據(jù)為0～1 ℃。溫度≤12 ℃及累計小時數(shù)分析中，棉花主產(chǎn)區(qū)溫度≤12 ℃時累計小時數(shù)均＞72 h，達到低溫冷害災(zāi)害指標。在累計持續(xù)小時數(shù)≤120 h 的天山北半坡，個別站點CLDAS 溫度數(shù)據(jù)表征持續(xù)小時數(shù)＞120 h。

通過此次霜凍災(zāi)害天氣過程分析和實地調(diào)查檢驗，在沒有地面觀測站點區(qū)域，可參考CLDAS 溫度數(shù)據(jù)代表當(dāng)?shù)販囟?。在關(guān)鍵農(nóng)事季節(jié)時，農(nóng)戶需做好農(nóng)作物人工管理并及時關(guān)注天氣預(yù)報。在極端降溫天氣情況下，尤其是作物播種出苗期間，農(nóng)戶需關(guān)注溫度是否≤0 ℃、≤12 ℃的持續(xù)時間及降溫幅度，CLDAS 溫度降至2～3 ℃時，已達作物開始受凍上限溫度，建議立即采取預(yù)防措施。