中國農(nóng)業(yè)氣象 Chinese Journal of Agrometeorology 2020 年 doi 10 3969 j issn 1000 6362 2020 12 005 張利平 周澤英 彭雲(yún) 等 四川渠江中上游羊肚菌低溫冷害風(fēng)險分析 J 中國農(nóng)業(yè)氣象 2020 41 12 794 806 四川渠江中上游羊肚菌低溫冷害風(fēng)險分析 張利平 1 2 周澤英 2 彭 雲(yún) 2 王慶莉 3 賴自力 4 王明田 1 5 6 1 中國氣象局成都高原氣象研究所 成都 610072 2 達州市氣象局 達州 635000 3 甘孜藏族自治州氣象局 康定 626000 4 瀘州市氣象局 瀘州 646000 5 南方丘區(qū)節(jié)水農(nóng)業(yè)研究四川省重點實驗室 成都 610066 6 四川省氣象臺 成都 610072 摘要 利用達州 巴中及周邊 21 個站點 1991 2020 年氣象資料 選取最低氣溫和積寒作為區(qū)劃指標(biāo) 分析渠 江中上游海拔 600m 以下地區(qū) 2 3 月羊肚菌子實體期間冷害指標(biāo)的時空分布 建立冷害指標(biāo)的地理空間分布 模型 基于氣象災(zāi)害風(fēng)險評價原理 運用 ArcGIS10 1 進行冷害風(fēng)險等級區(qū)劃 結(jié)果表明 1 羊肚菌冷害風(fēng) 險指標(biāo)為日最低氣溫 2 冷害風(fēng)險預(yù)警基礎(chǔ)指標(biāo)為日平均氣溫 6 0 2 羊肚菌冷害頻率為 10 年 7 遇 以上 主要發(fā)生在 2 月上旬 3 月上旬 發(fā)生頻次和強度逐旬減小 其中 2 月上中旬約占 70 冷害過程持續(xù) 日數(shù)以 1 3d 為主 占總過程數(shù)的 82 89 最長過程日數(shù) 6 8d 最大積寒 3 64 d 3 羊肚菌冷害風(fēng) 險等級自南向北 自低海拔向高海拔增大 以中低風(fēng)險區(qū)為主 中低風(fēng)險區(qū)主要分布在渠江中上游中南部海 拔 520m 以下的淺丘河谷平壩區(qū) 面積約占 2 3 次高和高風(fēng)險區(qū)分布零散 主要分布在中北部海拔 520 600m 的低山河谷區(qū) 面積約占 1 3 4 將羊肚菌原基分化期調(diào)控至 2 月下旬以后 冷害頻次降至約 3 年 2 遇 冷害強度明顯降低 其危害明顯減小 5 利用黑色遮陽網(wǎng)搭建的栽培設(shè)施會降低環(huán)境最低溫度約 0 3 加重冷害過程對羊肚菌的損害 渠江中上游低海拔地區(qū)羊肚菌順季栽培時應(yīng)控制在海拔 520m 以下地區(qū) 并 根據(jù)天氣狀況 將原基分化期調(diào)控至 2 月下旬以后 注意防范冷害過程中局地小氣候和黑色鏤空大棚的不利 影響 以避免或減輕冷害風(fēng)險及危害 關(guān)鍵詞 羊肚菌 低溫冷害 氣候風(fēng)險 渠江中上游 Risk Analysis of Chilling Injury of Morchella in Middle and Upper Reaches of Qu River in Sichuan Province ZHANG Li ping 1 2 ZHOU Ze ying 2 PENG Yun 2 WANG Qing li 3 LAI Zhi li 4 WANG Ming tian 1 5 6 1 Institute of Plateau Meteorology China Meteorological Administration Chengdu 6l0072 China 2 Dazhou Meteorological Bureau Dazhou 635000 3 Meteorological Bureau of Ganzi Tibetan Autonomous Prefecture Kangding 626000 4 Luzhou Meteorological Bureau Luzhou 646000 5 Water Saving Agriculture in Southern Hill Area Key Laboratory of Sichuan Province Chengdu 6l0066 6 Sichuan Metorological Observatory Chengdu 610072 Abstract Chilling injury is one of the main climatic problems affecting the cultivation and production of Morchella However there are few reports about the mechanism of chilling injury meteorological index and climate risk zoning of Morchella and there is lack of necessary technical support to carry out the weather services for Morchella chilling injury It provides scientific basis for the production planning of Morchella in the study area the dynamic assessment of the risk of chilling injury and the meteorological service of disaster prevention and reduction and facilitates the development of agricultural industry in mountainous areas Taking the upper and middle reaches of the Qu River below 600m above sea level as the study area DEM data of this area and climatic data of 21 stations in and around 收稿日期 2020 06 22 通訊作者王明田 E mail Wangmt0514 基金項目 高原與盆地暴雨旱澇災(zāi)害四川省重點實驗室科技發(fā)展基金項目 省重實驗室 2018 重點 05 05 第一作者聯(lián)系方式 張利平 E mail 854874815 第 12 期 張利平等 四川渠江中上游羊肚菌低溫冷害風(fēng)險分析 795 the basin from 1991 to 2020 and experimental data of Morchella production from 2018 to 2019 were used based on the chilling injury types and disaster causing factors of Morchella extreme minimum air temperature and accumulated cold were selected as zoning indices to analyze the relationship between air temperature and ground temperature from February and March in the study area to determine the threshold of extreme minimum air temperature to statistically analyze the temporal and spatial distributions of extreme minimum air temperature cold accumulation chilling injury frequency and its lasting days The spatial distribution data of index factors were standardized by establishing the spatial distribution model of chilling damage index Analytic hierarchy process and natural breakpoint method were used to determine the influence weight and risk level thresholds of chilling injury indexes respectively Based on the principle of comprehensive meteorological disaster risk assessment ArcGIS10 1 was used to make a classification in Morchella chilling risk The results showed that 1 the average difference between air and ground temperature from February and March was about 1 0 0 7 1 4 in the study area smaller in the south and larger in the north The risk index of chilling injury of Morchella was the daily minimum temperature that was less than or equal to 2 0 and the basic index of early warning of chilling injury risk was daily average temperature that was less than 6 0 2 The chilling injury frequency of Morchella was more than 7 times in 10 years mainly occurred from early February to early March and the occurrence frequency and intensity decreased every 10 days among which 70 occurred in early and middle February The duration of chilling injury process was 1 3 days accounting for 82 89 of the total number of processes the longest duration was 6 8 days and the maximum accumulated cold was 3 64 d 3 The risk levels of chilling injury of Morchella increased from south to north and from low altitude to high altitude mainly in the middle and low risk areas The middle and low risk areas are mainly distributed in the shallow hills and valleys below 520m above sea level in the middle and upper reaches of the Qu River accounting for about 2 3 of the total area The secondary high and high risk areas were scattered mainly distributed in the low mountain valley area with an altitude of 520 600m in the central and northern part accounting for about 1 3 of the area 4 After controlling the differentiation stage of Morchella to late February the chilling injury frequency was reduced to about 2 times in 3 years the chilling injury intensity was significantly reduced and the harm was significantly reduced 5 The planting site selection and the establishment of greenhouses had different effects on the chilling injury of Morchella the influence of local topography on microclimate can be judged comprehensively according to the altitude and whether it is conducive to the southward cold air flow and aggregation the cultivation facilities built with black sunshade net can reduce the environmental minimum temperature by about 0 3 and aggravate the chilling damage to Morchella In the study area the suitable season cultivation of Morchella should be controlled below 520m above sea level According to the weather conditions the primary basal differentiation period should be adjusted after late February Attention should be paid to the adverse effects of local microclimate and black hollowed out greenhouse in the process of chilling injury so as to avoid or reduce the risk and harm of chilling injury Key words Morchella Chilling injury Climate risk The middle and upper reaches of the Qu River 羊肚菌 Morchella 是一種珍稀的野生食用菌 1 有豐富的營養(yǎng)價值和很高的保健功能 中國 20 世紀(jì) 80 年代開始羊肚菌人工栽培馴化研究 2 3 2000 年 四川省林業(yè)科學(xué)院首次采用外營養(yǎng)添加技術(shù) 突破 羊肚菌的室外栽培技術(shù)難題 4 7 之后 羊肚菌優(yōu)良 菌株的引進與篩選 人工栽培技術(shù)的改進和商業(yè)化 應(yīng)用等得到快速發(fā)展 相關(guān)仿生態(tài)栽培技術(shù)已在國 內(nèi)多省 市 自治區(qū) 得到大面積推廣 在實際栽 培生產(chǎn)中 不同品種 不同地區(qū) 不同年份的羊肚 菌產(chǎn)量表現(xiàn)極不穩(wěn)定 主要原因有菌種不穩(wěn)定 抗 逆性差以及對環(huán)境條件的依賴性強等 其中最大的 制約因素是天氣與氣候條件的不穩(wěn)定性 4 羊肚菌主 要氣象災(zāi)害有冷害 熱害 干旱和低溫連陰雨等 其中影響最大的是冬末春初羊肚菌原基分化 子囊 中 國 農(nóng) 業(yè) 氣 象 第 41 卷 796 果生長期的低溫冷害 如四川省甘孜州康定市 2016 年春季發(fā)生的強寒潮低溫冷凍災(zāi)害 導(dǎo)致當(dāng)年羊肚 菌絕收 羊肚菌冷害是指當(dāng)環(huán)境溫度低于某一臨界值 時 子實體生長發(fā)育不良甚至受凍死亡的減產(chǎn)現(xiàn)象 國外對冷害機理研究較多 8 11 國內(nèi)冷害研究主要集 中在水稻 東北玉米 新疆棉花以及熱帶經(jīng)濟林木 王遠皓等 12 將冷害指標(biāo)歸納為生長季溫度距平 13 生長季積溫 14 生長發(fā)育關(guān)鍵期冷積溫 15 作物發(fā) 育期距平 16 熱量指標(biāo) 17 和綜合指標(biāo) 18 等六類 霍 治國等 19 21 分別研究確定了香蕉和荔枝 龍眼 蓮 霧苗的低溫寒害指標(biāo)等級 王春乙等 22 24 基于自然 災(zāi)害風(fēng)險理論對作物冷害風(fēng)險評估方法進行了深入 研究 梁立江等 25 29 利用 GIS 技術(shù)分別對東北水稻 湖北馬鈴薯 華南香蕉和荔枝 福建芒果等開展了 冷害風(fēng)險區(qū)劃 羊肚菌生產(chǎn)是一種新興的特色與設(shè) 施農(nóng)業(yè) 各地栽培歷史短 栽培總面積較小 目前 關(guān)于羊肚菌冷害的致害機理 氣象指標(biāo)和氣候風(fēng)險 與種植氣候區(qū)劃的研究少見報道 開展冷害氣象服 務(wù)缺乏必要的技術(shù)支撐 近年來 因產(chǎn)品需求旺盛 經(jīng)濟效益顯著 渠 江中上游羊肚菌栽培面積逐年增加 是山區(qū)農(nóng)戶脫 貧和農(nóng)村產(chǎn)業(yè)發(fā)展的主要項目之一 但冷害是影響 其安全生產(chǎn)和快速發(fā)展的主要氣候問題 本研究基 于現(xiàn)有糧經(jīng)作物冷害分析和風(fēng)險區(qū)劃研究方法與 成果 在不考慮社會經(jīng)濟與管理水平等影響的前提 下 以四川渠江中上游羊肚菌主要栽培區(qū)出菇期冷 害為研究對象 結(jié)合田間試驗和生產(chǎn)調(diào)查 研究羊 肚菌冷害致災(zāi)因子及風(fēng)險指標(biāo) 運用 GIS 技術(shù)開展 羊肚菌冷害氣候風(fēng)險精細區(qū)劃 以期為本地區(qū)羊肚 菌生產(chǎn)規(guī)劃 冷害動態(tài)評估和防災(zāi)減災(zāi)對策提供科 學(xué)依據(jù) 促進山區(qū)農(nóng)戶脫貧致富和農(nóng)村產(chǎn)業(yè)發(fā)展 1 材料與方法 1 1 研究區(qū)域與時段 渠江中上游地處四川東北部 大巴山 米倉山南 麓 106 35 108 54 E 30 32 32 74 N 海拔 222 0 2512 9m 主要包括達州和巴中兩市 地勢北高 大 巴山區(qū) 南低 盆地丘區(qū) 以山地和丘陵為主 屬 亞熱帶濕潤季風(fēng)氣候類型 由于地形復(fù)雜 山地氣 候特征明顯 海拔 800m 以下的低山 丘陵 河谷地 區(qū) 氣候溫和 四季分明 無霜期長 海拔 800 1000m 的低 中山區(qū) 氣候溫涼 陰濕 回春遲 夏日熱 秋涼早 冬寒長 海拔 1000m 以上的中山區(qū) 光熱 資源不足 寒冷期較長 春寒和秋霜突出 羊肚菌一般從低海拔的平原地區(qū)到 3200m 的高 海拔地區(qū)都有分布 3 但各地生育期不同 據(jù)調(diào)查和 試驗 表 1 目前渠江中上游區(qū)域內(nèi)羊肚菌主要種 植于海拔 300 600m 的河谷 低山淺丘及平壩區(qū) 600m 以上山區(qū)種植極少 且子實體生長季明顯延后 該區(qū)域內(nèi)羊肚菌栽培季一般為 10 中下旬 翌年 3 月 底或 4 月初 一般在 10 月中下旬播種 1 月下旬 2 月上旬為原基分化期 3 月上旬隨著子實體生長成 熟 開始采收 4 月上中旬結(jié)束 采收期 30d 左右 可見 2 3 月為羊肚菌生長的關(guān)鍵期 據(jù)此 確定研 究區(qū)域為渠江中上游海拔 600m 以下地區(qū) 以下稱 研究區(qū) 研究統(tǒng)計時段為 2 3 月 研究區(qū)地形與 資料站點分布見圖 1 表 1 研究區(qū)內(nèi)羊肚菌生育期調(diào)查結(jié)果 Table 1 Investigation result of the growth period of Morchella in research station 地點 Site 海拔 Altitude m 播種期 Sowing stage 原基分化 Primordium differentiation 收獲始期 Begin of harvest 收獲盛期 Bloom of harvest 收獲末期 End of harvest 渠縣 Quxian 295 1 10 月下旬末 Late in L Oct 1 月下旬初 Early in L Jan 3 月上旬初 Early in E Mar 3 月中旬初 Early in M Mar 3 月末 End in Mar 宣漢 Xuanhan 389 4 10 月下旬初 Early in L Oct 2 月上旬初 Early in E Feb 3 月上旬中 Middle in E Mar 3 月中 M Mar 4 月初 E Apr 南江 Nanjiang 579 3 10 月中旬中 Middle in M Oct 2 月上旬末 Late in E Feb 3 月上旬末 Late in E Mar 3 月中旬末 Late in M Mar 4 月上旬末 Late in E Apr 萬源 Wanyuan 674 0 10 月中旬初 Early in M Oct 2 月下旬初 Early in L Feb 3 月中 M Mar 3 月下旬末 Late in L Mar 4 月中 M Apr Note E is the first ten day of a month M is the middle ten day of a month L is the last ten day of a month The same as below 第 12 期 張利平等 四川渠江中上游羊肚菌低溫冷害風(fēng)險分析 797 圖1 研究區(qū)域 海拔600m以下 地形與站點分布 Fig 1 Topography and site distribution in the study area altitude 600m 1 2 資料及其來源 DEM 高程數(shù)據(jù) 分辨率 30m 的 DEM 高程數(shù)據(jù) 來源于地理空間云 1 100 萬行政區(qū)劃邊界 shp 2017 版 數(shù)據(jù)來源于國家基礎(chǔ) 地理信息中心 用于 空間分析 氣候數(shù)據(jù) 渠江中上游及其周邊 21 個縣級氣象 站 1991 2020 年逐日地面氣象觀測資料來源于四川 省氣象探測數(shù)據(jù)中心 重慶市氣象科學(xué)研究所和陜 西省漢中市氣象局 用于對研究區(qū)域內(nèi)典型站點冷 害指標(biāo)的統(tǒng)計分析 試驗數(shù)據(jù) 2018 年 10 月 2019 年 4 月 在宣漢 縣雙河鎮(zhèn)方斗村 31 463 N 107 772 E 494 0m 和 萬源縣井溪鎮(zhèn)鹽井壩村 31 772 N 108 229 E 500 0m 同時開展羊肚菌生產(chǎn)試驗 方斗村為丘陵 臺地 鹽井壩村為河谷 土壤類型均為水稻土 試 驗在遮光率為 75 黑色遮陽網(wǎng)搭建的 2m 高簡易大 棚中進行 試驗品種為六妹羊肚菌 M sextelata 栽 培種 每個小區(qū)面積為 10m 2 設(shè) 1 個重復(fù) 5 個播期 播期間隔 5d 首播期為 10 月 29 日 試驗區(qū)外生產(chǎn) 對照地塊播期為 11 月 8 日 棚內(nèi)安裝便攜式自動氣 象觀測儀 高度 1 5m 和土壤多參數(shù)速測儀 外接 電源 主要觀測氣溫 空氣濕度 日照時數(shù)以及 0 10cm 土壤溫濕度 棚外 800m 范圍內(nèi)建有 1 個區(qū) 域自動氣象站 自動定時收集氣象資料 結(jié)合羊肚 菌生育期觀測 分析氣象條件和栽培設(shè)施內(nèi)外溫差 及其變化對羊肚菌生長發(fā)育的影響 1 3 冷害指標(biāo)選擇 羊肚菌子實體期發(fā)生冷害 原基和子囊果會因 凍而受損或死亡 受害程度與地面最低氣溫直接相 關(guān) 還與降溫幅度和持續(xù)時間等密切相關(guān) 根據(jù)降 溫天氣條件的不同 冷害可分成輻射型冷害 平流 型冷害和混合型冷害 30 其中 輻射型冷害主要由 極端最低氣溫造成 平流型冷害主要由持續(xù)低溫造 成 混合型冷害主要是由平流型冷害轉(zhuǎn)為輻射型冷 害所致 冷害致災(zāi)因子可以概括為極端最低氣溫 最大降溫幅度 持續(xù)日數(shù)和有害積寒 降溫幅度和 極端最低氣溫可較好地表示強冷空氣入侵引起的上 述 3 種冷害類型的劇烈程度 低溫持續(xù)日數(shù)和過程 有害積寒可較好地表示中弱冷空氣多次補充造成的 平流型冷害的累積作用 據(jù)研究 27 過程有害積寒與過程持續(xù)日數(shù) 最 低溫度和降溫幅度呈顯著相關(guān) 能綜合反應(yīng)不同類 型冷害過程的時長與強度 因此 選擇極端最低氣 溫和過程有害積寒作為羊肚菌的冷害指標(biāo) 1 4 統(tǒng)計口徑與計算方法 1 羊肚菌低溫冷害過程 羊肚菌原基分化 子 實體期間 當(dāng)日最低氣溫 T min 冷害臨界溫度 T c 且持續(xù) d 或以上時 定義為 1 次羊肚菌低溫 冷害過程 2 冷害頻率 P 統(tǒng)計年 N 內(nèi)冷害發(fā)生次 數(shù) n 的百分比 P n N 100 1 3 冷害過程日平均氣溫 T k 冷害過程中日 平均氣溫的平均值 T k T m 2 式中 T 為最低氣溫 界限溫度之日的平均氣 溫 m 為最低氣溫 界限溫度的日數(shù) d 4 過程有害積寒 ACHT accumulated cold harmful temperature 簡稱積寒 指冷害過程中 低 于臨界溫度的逐時溫度與臨界溫度之差的絕對值累 積量 計算方法 根據(jù)確定的羊肚菌發(fā)生低溫冷害臨 界溫度指標(biāo) 統(tǒng)計低溫冷害過程數(shù)據(jù) 由于沒有逐 時溫度資料 參考文獻 19 中公式 計算每次低溫冷 中 國 農(nóng) 業(yè) 氣 象 第 41 卷 798 害過程的積寒值 計算式為 2 min c 2 x cmin min c m1 mmin 0 T T A TT 1 T T 4TT 采用層次分析 法確定 1 7 空間分布數(shù)理模型建立 運用多元回歸法建立各指標(biāo)或要素的空間分布 數(shù)理模型 即 ii i i i Aa b ch d 8 式中 i A 為第 i 項氣象災(zāi)害指標(biāo)的回歸頻次 h 分別代表經(jīng)度 緯度和海拔 i a i b i c 為 回歸系數(shù) d i 為常數(shù) 為殘差 利用 ArcGIS10 1 繪圖 殘差采用反距離權(quán)重法 插值訂正 1 8 風(fēng)險等級劃分與區(qū)劃圖繪制 利用 ArcGIS10 1 的空間分析模塊制作區(qū)劃圖 采用層次分析法確定冷害指標(biāo)的權(quán)重系數(shù) 分別對 冷害指標(biāo)圖層進行標(biāo)準(zhǔn)化處理 按權(quán)重對指標(biāo)圖層 進行疊加計算 繪制風(fēng)險柵格圖層 按照自然斷點 法確定低 次低 中 次高 高 5 級風(fēng)險分級閾值 對風(fēng)險柵格圖層分級 最后 疊加行政邊界和主要 河流等地圖信息 繪制冷害風(fēng)險區(qū)劃圖 2 結(jié)果與分析 2 1 羊肚菌低溫冷害指標(biāo)確定 2 1 1 地溫臨界值 羊肚菌是低溫菌 不同品種的適宜生長發(fā)育溫 度范圍略有差異 2 4 31 32 一般而言 菌絲體發(fā)育溫 度在 5 0 30 0 均可 最適溫度為 15 0 22 0 能耐 0 0 以下低溫 當(dāng)溫度穩(wěn)定在 8 0 以上 3d 后 原基開始分化 子囊果含水量在 90 以上 環(huán)境溫 度在 6 0 25 0 區(qū)間均能生長 低于 15 0 品質(zhì)較 好 原基分化 子囊果初期是低溫敏感期 研究區(qū)多 處于冬末春初 該時段內(nèi)環(huán)境溫度為 0 0 6 0 時 原基與子囊果生長緩慢甚至停止分化 生長 當(dāng)環(huán) 境溫度 0 0 時 原基或子囊果會因受凍結(jié)冰而死 亡 嚴(yán)重時導(dǎo)致羊肚菌大量減產(chǎn)甚至絕收 羊肚菌 主要生長在地表 0 15cm 范圍內(nèi) 可用地表溫度表示 其生長的環(huán)境溫度 因此 確定 2 3 月羊肚菌低溫 冷害的地面溫度臨界值為 0 0 2 1 2 氣溫臨界值 地溫受環(huán)境影響因素較多 加之一般種植地區(qū) 并沒有觀測或預(yù)報 所以實際應(yīng)用中地溫不如氣溫 H 1 5m 方便 因此 需將 0 0 的地面溫度臨界 值轉(zhuǎn)換為對應(yīng)的氣溫值 1 以南江 宣漢和大竹 3 個國家氣象觀測站 分別代表研究區(qū)的北 中 南 3 個區(qū)位 選取 2010 2019 年 2 3 月日最低氣溫 T min 3 0 觀測 資料與同日地面最低溫度 T 0min 觀測資料為一個 數(shù)據(jù)組 T m 3 0 的頻率在 98 以上 利用上述資 料進行線性回歸分析 結(jié)果見表 2 第 12 期 張利平等 四川渠江中上游羊肚菌低溫冷害風(fēng)險分析 799 由表 2 可見 研究站點最低氣溫與地面最低溫 度呈極顯著正相關(guān) 但回歸系數(shù)較小 0 1 常數(shù) 項 平均溫差 在 0 90 1 35 從南至北增大 最 低氣溫與地面最低溫度之差相對恒定 2 統(tǒng)計同日最低氣溫 T min 與地面最低溫度 T 0min 之差 T m 和當(dāng) T min 臨界值 臨界值取 0 0 1 0 2 0 和 3 0 時 T 0min 0 0 的頻率 臨界值 結(jié)果見表 3 由表中可見 日最低氣溫與地面最低溫 度關(guān)系較復(fù)雜 二者平均差值在 0 7 1 4 溫差范 圍在 1 9 3 5 其中最低氣溫高于地面最低溫度 的頻率為 82 2 94 1 當(dāng)最低氣溫 0 0 時 一 般都有冷害發(fā)生 最低氣溫 1 0 時 冷害頻率 50 0 91 0 最低氣溫 2 0 時 冷害頻率 41 2 77 3 最低氣溫 3 0 時 冷害頻率 31 1 59 3 3 統(tǒng)計當(dāng) T 0min T min 時 T min 小于 1 0 2 0 和 3 0 等臨界值的頻率 臨界值 和溫差 T m 結(jié)果見 表 4 由表中可見 當(dāng)冷害發(fā)生時 最低氣溫高于地 面最低溫度的頻率在 90 以上 平均溫差為 1 2 1 6 最低氣溫 2 0 的頻率均在 95 0 以上 從地理區(qū)位上看 隨著緯度北移和海拔升高 研究區(qū)內(nèi)冷害日數(shù)明顯增加 且最低氣溫與地面最 低溫度之差增大 表 2 代表站點 2 3 月日最低氣溫 T min 與地面最低溫度 T 0min 間回歸分析結(jié)果 2010 2019 年 Tble 2 Regression analysis between the minimum air temperature T min and minimum ground temperature T 0min at three representative stations from Feb to Mar 2010 2019 站點 Station 數(shù)據(jù)組 Data array 0 a 1 a R F 南江 北部 Nanjiang northern 101 1 35 0 09022 0 847 257 2 2 2 2 0 宣漢 中部 Xuanhan middle 62 1 11 0 06608 0 730 68 49 2 0 1 3 大竹 南部 Dazhu southern 45 0 90 0 06127 0 678 36 63 1 6 1 8 注 min 0 1 0min TaaT 表示方程通過 0 01 水平顯著性檢驗 下同 Note min 0 1 0min TaaT indicated that the equation passed the significance test of P 0 01 The same as below 表 3 不同臨界氣溫水平下地溫 T 0min 0 的頻率及不同條件下的氣 地溫差統(tǒng)計 2010 2019 年 Table 3 Statistics of the ground air temperature difference and frequency of T 0min 0 under different condition 2010 2019 當(dāng) T min 不同值時 T 0min 0 的頻率 Frequency of T 0min 0 when T min different value T min T 0min T min T 0min 站點 Station 海拔 Altitude m 數(shù)據(jù)組 Data array P 0 0 P 1 0 P 2 0 P 3 0 溫差 標(biāo)準(zhǔn)差 T m SD P 1 T m1 1 P 2 T m2 2 南江 Nanjiang 579 3 101 100 0 91 0 77 3 59 3 1 4 0 84 94 1 1 4 3 9 0 6 宣漢 Xuanhan 389 4 62 100 0 77 8 52 3 36 5 0 9 0 85 83 9 1 2 16 1 0 4 大竹 Dazhu 418 0 45 80 0 50 0 41 2 31 1 0 7 0 84 82 2 0 9 17 8 0 4 注 T m T min T 0min T m 為同日最低氣溫 T min 與地面最低溫度 T 0min 之差的平均值即氣 地溫差 T m1 和 T m2 為不同條件 下的溫差 P P 1 和 P 2 為相應(yīng)的頻率 下同 Note T m T min T 0min T m is the average difference between T min and T 0min i e temperature difference and T m1 and T m2 are mean temperature difference under different condition P P 1 and P 2 are their frequency SD is standard deviation The same as below 表 4 冷害發(fā)生時的氣 地溫差 和不同臨界氣溫下 T 0min 0 的頻率 2010 2019 年 Table 4 Statistics of the ground air temperature difference and frequency under different critical air temperature when chilling injury occurs in 2010 2019 T min T 0min T min T 0min 站點 Station 數(shù)據(jù)組 Data array P 1 0 P 2 0 P 3 0 溫差 標(biāo)準(zhǔn)差 T m SD P 1 T m1 P 2 T m2 南江 Nanjiang 43 77 7 95 1 100 0 1 5 0 83 94 8 1 6 5 2 0 3 宣漢 Xuanhan 21 61 9 95 2 100 0 1 4 0 82 90 5 1 5 9 5 0 3 大竹 Dazhu 14 64 3 100 0 100 0 1 0 0 81 92 9 1 2 7 1 0 7 中 國 農(nóng) 業(yè) 氣 象 第 41 卷 800 2 1 3 低溫冷害指標(biāo)確定 統(tǒng)計方斗和鹽井壩站羊肚菌試驗大棚內(nèi)外 2019 年 2 月的最低氣溫和地溫情況 結(jié)果見表 5 由表中 可見 兩個試驗大棚內(nèi)外 2 月平均最低氣溫與平均 最低地面溫度之差 0 6 0 9 極端最低氣溫與極端 最低地面溫度之差 1 4 1 7 棚內(nèi)外溫差變化方向 與此一致但有差異 2 月 19 20 日 渠江中上游發(fā)生了一次輻射型 冷害天氣過程 18 日下午天氣轉(zhuǎn)晴以后 19 日試驗 大棚內(nèi)外平均氣溫 5 0 左右 比前一天降低約 1 5 而極端最低氣溫由前一天的 2 0 左右下降 到 0 0 以下 地面極端最低溫度降至 2 5 1 7 差值 1 4 1 7 20 日 方斗站冷害結(jié)束 鹽井壩 站平均氣溫回升約 1 0 極端最低氣溫僅回升 0 2 0 4 棚內(nèi)地面極端最低溫度比棚外偏低 0 4 0 3 宣漢站 2019 年 2 月 19 日冷害日的氣溫和地溫 實時變化情況見圖 2 由圖可見 受云天狀況變化 影響 當(dāng)日 2 00 后地溫和氣溫迅速下降 4 00 5 00 先后降至最低 然后逐步回升 8 00 后受太 陽輻射影響 地溫迅速上升 日極端最低氣溫 0 8 與極端最低地面溫度 1 0 之差為 1 8 綜上 研究區(qū) 2 3 月最低氣溫與地面最低溫度 的平均差值在 0 7 1 4 當(dāng)最低氣溫 2 0 時 雖不一定出現(xiàn)冷害 但當(dāng)有冷害出現(xiàn)時 最低氣溫 一般在 2 0 以下 綜合考慮羊肚菌風(fēng)險管理 栽培 環(huán)境復(fù)雜性等 確定 2 3 月羊肚菌冷害風(fēng)險指標(biāo)為 日最低氣溫 2 0 2 2 羊肚菌冷害發(fā)生情況統(tǒng)計 按照日最低氣溫 2 0 標(biāo)準(zhǔn) 統(tǒng)計羊肚菌 2 3 月冷害發(fā)生情況 結(jié)果見圖 3 表 6 表 8 由圖 3 和 表 6 可見 不同年份各地冷害差異明顯 除個別年 表 5 2019 年 2 月兩個試驗站簡易大棚內(nèi)外氣溫與地面溫度 Table 5 Air temperature and ground temperature 0cm inside an