壽寧縣高山茶春季霜凍特征及預報模型構(gòu)建

王政 王明輝 謝華清

摘要 利用壽寧縣1980—2019年3—4月逐日最低氣溫、14：00本站氣壓及其24 h變壓、14：00氣溫及其24 h變溫、14：00相對濕度、風向和風速等常規(guī)氣象資料，運用線性傾向估計法、t檢驗法和Pearson相關(guān)系數(shù)法等氣候統(tǒng)計診斷方法，分析春季霜凍害氣候特征;使用SPSS 19.0軟件對日最低氣溫主要相關(guān)因子進行多元回歸分析，建立日最低氣溫預報數(shù)學模型。結(jié)果表明，壽寧縣春季霜凍呈始日提前、終日延后的變化趨勢，預示未來春季霜凍期將延長，茶葉春季霜凍害風險增大。經(jīng)檢驗模型定性預報準確率達91.0%，預報效果較好，可應(yīng)用于氣象預報業(yè)務(wù)。

關(guān)鍵詞 茶葉;春季霜凍;氣象預報;數(shù)學模型

中圖分類號：S425 文獻標識碼：B 文章編號：2095–3305（2021）12–0038–03

茶葉是原產(chǎn)于中國的特色經(jīng)濟作物，因其生態(tài)綠色、優(yōu)質(zhì)高效的優(yōu)勢而在中國農(nóng)業(yè)生產(chǎn)中占有重要地位[1]。隨著全球氣候變暖加劇，壽寧縣各類氣象災害頻發(fā)、重發(fā)，對農(nóng)業(yè)生產(chǎn)造成嚴重影響。在眾多氣象災害中，春季霜凍害對茶葉生產(chǎn)影響尤為嚴重，春季正是茶樹新葉萌芽生長期，茶葉芽和嫩梢抗低溫能力較弱，此時遭遇霜凍害，輕則造成茶樹新芽葉尖受凍損傷，茶葉產(chǎn)量和品質(zhì)下降;重則造成茶樹枝枯葉落，甚至死亡。因此提高春季霜凍氣象的預報準確率，降低災害造成的損失，對茶葉生產(chǎn)具有十分重要的意義。

20世紀70年代以來，國內(nèi)許多學者對茶葉霜凍害做過大量研究。吳楊等[2]以日最低氣溫≤4℃為霜凍指標，根據(jù)茶葉受災癥狀進行分級，確定了不同等級的氣象指標。唐力生等[3]通過分析氣溫直減率、輻射與溫度的關(guān)系，建立了寒害過程日平均氣溫分布式模型。這些研究成果和方法對茶葉春季霜凍災害的研究提供了有益的借鑒。通過研究福建省壽寧縣春季霜凍害的氣候規(guī)律特征，構(gòu)建春季霜凍預報數(shù)學模型，建立茶葉生產(chǎn)與霜凍災害的氣象預警機制，為提高茶葉種植效益和防災水平提供技術(shù)支撐，為鄉(xiāng)村振興提供精準專業(yè)氣象服務(wù)。

1 資料與方法

1.1 研究區(qū)概況

福建省壽寧縣位于27°16′N～ 27°41′N、119°12′E～119°44′E，地處鷲峰山系洞宮山脈東麓，距離太平洋西海岸約83 km，地勢從西北向東南傾斜，境內(nèi)層巒疊嶂、地勢起伏，海拔落差大，海洋性季風氣候明顯。該縣氣候溫暖，四季分明，雨量充沛，光照充足，無霜期長，屬中亞熱帶海洋性季風氣候，常年平均氣溫15.1℃，年日照時數(shù)1 663.3 h，年降水量1 976.9 mm，年無霜期238 d。

1. 2 資料來源與方法

選取壽寧縣氣象局氣象資料數(shù)據(jù)庫中1980—2019年3—4月的逐日最低氣溫和在14：00的氣溫、本站氣壓、相對濕度、風向、風速等氣象觀測數(shù)據(jù)，所有資料序列完整、連續(xù)、均一。運用線性傾向估計法、t檢驗法和Pearson相關(guān)系數(shù)法等氣候統(tǒng)計診斷方法，分析春季霜凍害氣候特征，使用SPSS 19.0軟件對日最低氣溫的主要相關(guān)因子進行多元回歸分析，并建立日最低氣溫預報數(shù)學模型[4]。

2 結(jié)果與分析

2. 1 春季霜凍的氣候特征

2.1.1 春季霜凍日數(shù)的年際和年代際變化 1980—2019年壽寧縣春季霜凍日數(shù)線性擬合統(tǒng)計結(jié)果顯示，近40年壽寧春季霜凍日數(shù)呈不顯著減少趨勢，減少速率為0.423 d/10 a。年平均春季霜凍日數(shù)10.63 d，其最多年份達23 d（2011年），最少年份僅4 d，為1992年、1997年、2002年，說明春季霜凍日數(shù)年際分布差異大。從春季霜凍日數(shù)年代平均值曲線看，20世紀80年代春季霜凍日數(shù)最多，為12.1 d，其次是21世紀10年代（10.6 d）、20世紀90年代（10.0 d）和21世紀00年代（9.8 d），春季霜凍日數(shù)年代際振幅為2.3 d，年代際分布較為均衡，年代際差異不明顯（圖1）。

2.1.2 春季霜凍日的起止時間變化 每年春季首次出現(xiàn)的霜凍日稱為春季霜凍始日，最后一次出現(xiàn)的霜凍日稱為春季霜凍終日。近40年壽寧縣春季霜凍始日最早在3月1日，最晚在4月23日，5月無霜凍日。將3月1日的日序記為1，3月2日的日序記為2……以此類推，4月30日的日序記為61，以此法建立1980—2019年春季霜凍始日和終日的日期序列。

1980—2019年壽寧縣春季霜凍始日和終日日序的線性擬合統(tǒng)計結(jié)果顯示，近40年來春季霜凍始日平均為3月6日，最早在3月1日，最晚是3月15日;春季霜凍終日平均出現(xiàn)在4月1日，最早在3月5日，最晚在4月23日（圖2）。春季霜凍始日呈不顯著提前趨勢，提前速率為0.341 d/10 a ;春季霜凍終日呈不顯著延后趨勢，延后速率為0.591 d/10 a ，且終日延后速率大于始日提前速率，說明未來春季霜凍期將延長，茶葉春季霜凍害風險增大。

2.2 春季霜凍預報

2.2.1 霜凍預報因子的選取 霜凍及其強度的預報，即對日最低氣溫的預報。選擇日最低氣溫值作為霜凍預報量，根據(jù)霜凍形成機理和當?shù)靥鞖忸A報工作經(jīng)驗，選取8個與日最低氣溫相關(guān)且物理意義明顯的氣象要素作為預報因子，分別是：霜凍日前一天的日最低氣溫、14：00本站氣壓及其24 h變壓、14：00氣溫及其24 h變溫、14：00相對濕度、風向和風速。

統(tǒng)計分析1980—2009年3—4月的逐日日最低氣溫與8個預報因子的相關(guān)性，結(jié)果顯示，8個預報因子中，僅風向和風速2個預報因子相關(guān)系數(shù)未通過顯著性檢驗，其他6個預報因子相關(guān)系數(shù)均通過0.01或0.05的顯著性檢驗（表1）。表明這6個預報因子與預報量相關(guān)性顯著，可用于預測日最低氣溫值。6個預報因子中，本站氣壓和24 h變壓2個預報因子與預報量呈負相關(guān)，其他4個預報因子與預報量均呈正相關(guān)。說明前一日14：00本站氣壓值及其24 h變壓值越大，則次日的日最低氣溫越低;反之，則次日的日最低氣溫越高。而前一日的日最低氣溫值、14：00氣溫值及其24 h變溫值、14：00相對濕度值越大，則次日的日最低氣溫越高;反之，則次日的日最低氣溫越低。

2.2.2 霜凍預報方程的建立 根據(jù)以上分析，對選取的6個預報因子，運用多元線性回歸方法建立日最低氣溫預報的數(shù)學模型，其回歸方程為：

Y=0.192X1-0.307X2+0.555X3- 0.097X4+0.328X5+0.07X6-184.952? ? ? ? （1）

式（1）中，Y為預報日的日最低氣溫，X1、X2、X3、X4、X5、X6分別為預報日前一日14：00的本站氣壓、24 h變壓、氣溫、24 h變溫、日最低氣溫、相對濕度。方程的方差檢驗F= 398.924，通過了α=0.01（F0.01=29.343）顯著性檢驗。

2.2.3 預報方程檢驗 利用公式（1）對2010—2019年3—4月逐日最低氣溫值進行回報檢驗，結(jié)果顯示，公式（1）預報是否出現(xiàn)霜凍的準確率達91.0 %，定性預報準確率滿意（表2）。

3 結(jié)論

（1）1980—2019年壽寧縣年平均春季霜凍日數(shù)10.63 d。近40年來壽寧縣春季霜凍日數(shù)呈略減少趨勢，減少速率為0.423 d/10 a。春季霜凍始日呈不顯著提前趨勢，提前速率為0.341 d/10 a，春季霜凍終日趨于略延后，延后速率為0.591 d/10 a?？梢娢磥韷蹖幙h春季霜凍期將延長，茶葉春季霜凍害風險增大。

（2）建立春季日最低氣溫與前一日14：00本站氣壓、24 h變壓、14：00氣溫、24 h變溫、日最低氣溫、14：00相對濕度等6個氣象要素的多元回歸方程，經(jīng)回報檢驗，方程定性預報準確率達91.0%，預報效果較好，可以作為壽寧縣春季霜凍預報模型在氣象業(yè)務(wù)中運用。

參考文獻

[1] 金志鳳，胡波，嚴甲真，等.浙江省茶葉農(nóng)業(yè)氣象災害風險評價[J].生態(tài)學雜志，2014，33（3）：771-777.

[2] 吳楊，金志鳳，葉建剛，等.浙江茶樹春霜凍發(fā)生規(guī)律及其與太平洋海溫的遙相關(guān)分析[J].中國農(nóng)業(yè)氣象，2014，35 （4）：434-439.

[3] 唐力生，杜堯東，陳新光，等.廣東寒害過程溫度動態(tài)監(jiān)測模型[J].生態(tài)學雜志，2009，28 （2）：366-370.

[4] 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預測技術(shù)[M].北京：氣象出版社，1999.

責任編輯：黃艷飛

Frost Characteristics and Prediction Model Construction of Shouning Gaoshan Tea in Spring

WANG Zheng et al（Shouning County Meteorological Bureau， Fujian Province， Shouning， Fujian 355500）

Abstract Conventional meteorological data including daily minimum temperature， 14：00 air pressure and its 24 h variable pressure， 14：00 air temperature and its 24 h variable temperature， 14：00 relative humidity， wind direction and wind speed were used in Shouning county from 1980 to 2019; and SPSS software was used to make multiple regression analysis on the main correlation factors of daily minimum temperature， and the mathematical model of daily minimum temperature forecast was established. The results showed that the beginning date of spring frost in Shouning county was earlier than the beginning date of spring frost， but the end date of spring frost was later， which indicated that the frost period in spring would be prolonged in the future， and the risk of tea frost injury in spring would be increased. The qualitative prediction accuracy of the tested model is 91.0%， and the prediction effect was good， which could be used in meteorological forecasting.

Key words Tea; Spring frost; Weather forecast; Mathematical model