清鎮市櫻花始花期預測分析

楊濤 陳宇 張皓 蔡志穎

摘 要 利用貴州省清鎮市2013—2020年逐日氣象數據及櫻花花期觀測資料，對清鎮櫻花始花期特征及其氣象影響因子進行了分析。結果表明：清鎮櫻花始花期平均日序為68.6 d，對應日期為3月8—9日，平均花期可達17.6 d，始花期變化幅度較大，最大可達9 d;櫻花花期與各階段氣溫、最高溫度、日照時數緊密相關，尤其是與2月下旬最高氣溫存在極顯著負相關;采用逐步回歸分析方法，得到了清鎮櫻花始花期預測模型并進行了歷年回測，發現預測始花期與實際始花期平均絕對誤差為1.0 d，最大誤差為3 d，模型檢驗效果較好。采用積溫分析法得到了清鎮櫻花始花期預測指標為1月15日起，日最高溫度累積值≥611.8 ℃時，即為清鎮櫻花花開日;用該指標同樣進行了歷年回測，發現始花期平均預測誤差天數為1.5 d，差距為2 d之內的預測準確率達88%。所建立的預測模型能夠較準確地預測出清鎮櫻花始花期。

關鍵詞 櫻花;花期;氣象因子;預測模型;貴州省清鎮市

中圖分類號：S162.5+5 文獻標志碼：B DOI：10.19415/j.cnki.1673-890x.2021.31.043

觀賞性植物的始花期預測已成為近年來倍受關注的焦點問題之一，植物物候期的變化與氣象條件密切相關。研究表明，諸多環境要素（溫度、降水、光照、風速、養分等）都會影響植物物候期，其中溫度作用最大[1]。饒紅欣等利用谷城櫻花開花前1～2月的最高溫度建立了始花期線性回歸方程，結果表明當年1—2月的平均最高氣溫升高1℃，櫻花始花期可提前2.4 d[2]。肖丹等對眉山晚櫻盛花期和氣溫進行了相關性分析發現，櫻花盛花期受上年11月、12月氣溫影響較小，與當年1—3月平均氣溫的相關性最高[3]。櫻花花期與溫度呈顯著負相關，開花前一段時間內溫度越高，花期時序越短。陳正洪等用武漢大學櫻花園連續62年的花期資料研究了櫻花始花期、落花期、持續天數的氣候變化特征，并建立了始花期預報模型，發現始花期與冬季及2月平均氣溫密切相關，其中基于2月平均氣溫的非線性模型對異常花早開有較好的模擬效果[4]。魏華兵等利用積溫方法研究了櫻花花期，發現下限溫度和不同起始日期都會影響花期積溫指標[5]。

本文研究清鎮氣象條件對櫻花始花期的影響，以期為當地有關部門提供服務決策依據，指導人們合理安排時間觀賞櫻花，助力鄉村旅游發展。

1 ?資料與方法

1.1 ?研究區概況

清鎮市位于貴州省中部，東經106°07′～106°33′，北緯26°21′～26°59′，總面積1 386.6 km2，屬北亞熱帶季風濕潤氣候，平均海拔1 200 m，年平均氣溫14.3 ℃，其中夏季平均氣溫22.7 ℃，年均降雨量1 168.9 mm，年均日照時數1 228.8 h。氣候溫和，舒適宜人，季風明顯，總儲水量近20億立方米。清鎮櫻花園位于清鎮市紅楓湖上游，緊連貴安新區，共同占地1 300余公頃的農場，種植了近70萬株名貴櫻花，其面積之廣，密度之大，享有“貴州最佳櫻花觀賞區”的美譽。每年3月櫻花盛開時期都會吸引上百萬人參觀游覽，為當地經濟發展作出了巨大貢獻。

1.2 ?數據來源

氣象數據來源于貴州省信息中心，選用清鎮市2013—2020年逐日平均氣溫、最高氣溫、降水、日照時數、平均風速等氣象數據，并采取調查及專業觀測方法，參照《農業氣象觀測規范》獲取了清鎮櫻花2013—2020年花期的相關資料。

1.3 ?研究方法

1.3.1 ?數據處理

首先將始花期、落花期轉換為日序數（每年1月1日記為1，1月2日記為2，……以此類推），其次計算始花期至落花期的持續天數作為花期天數，從而得到了始花期、落花期、花期天數等3個物候期的完整時序資料。

1.3.2 ?積溫計算方法

櫻花花芽的發育與氣溫有關，櫻花花期的相關研究表明櫻花生長過程需要一定的總熱量，總熱量常用總積溫或有效積溫來表示[6]。對櫻花花開前期總熱量的表述，除了用日均有效積溫表示外，日本預測櫻花多利用日最高氣溫累加值的方法來表示并分析櫻花花期影響指標。本研究采用一定溫度條件下的日平均氣溫或日最高氣溫表示，如下：

[Q=i=1mTT=0ttc] ? （1）

式（1）中，Q為花開前期所需一定條件的階段積溫，m為一定溫度條件下的發育歷期，T為發育歷期間的階段溫度值，c為影響發育的溫度閾值。

1.3.3 ?統計分析

應用SPSS26.0、Excel軟件對數據進行統計分析，用相關分析法分析清鎮櫻花始花期與氣溫、日照、風速等氣象因子的相關性，利用逐步回歸分析法找出對櫻花始花期有顯著影響的關鍵氣象因子并建立預報模型，其中2013—2020年的資料用于模型建立與回測，2021年的資料則作為獨立樣本進行了預測檢驗。

2 ?結果與分析

2.1 ?清鎮櫻花花期特征分析

根據2013—2020年清鎮櫻花的物候資料，發現清鎮櫻花始花期平均日序為68.6 d，對應日期為3月8—9日，平均盛花期為3月14日，平均落花期為3月26日，平均花期可達17.6 d。清鎮櫻花始花期最早為2013年3月5日，最晚為2017年3月14日，變幅長達9 d。借助均方差的大小判定始花期的穩定性程度[7]，發現櫻花始花期日序數均方差為3.2，超過63%的年份均偏離平均始花期2 d以上，這表明清鎮櫻花始花期近8年來并不穩定。

2.2 ?櫻花始花期與氣象因子的關系

2.2.1 ?影響櫻花始花期的氣象因子

將清鎮歷年1、2 月的日平均氣象資料換算成逐旬并逐一進行相關性分析。表1為清鎮櫻花始花期序列與各氣象因子的相關系數，可以看出氣溫是影響櫻花花期的主要氣象因子，其中1月中旬、2月下旬平均氣溫，1月下旬日照時數，2月最高氣溫與櫻花始花期呈顯著負相關關系，相關系數分別為-0.72、-0.83、-0.73、-0.76，尤其是2月下旬最高氣溫與櫻花花期的相關性達極顯著水平（a=0.01），這表明清鎮櫻花開花前，2月下旬的氣溫對花期影響最大，即2月下旬的最高氣溫越高，櫻花開花越早。從表中可以看出降水量、平均風速和最大風速均未通過顯著性檢驗，這表明此三項氣象因子對清鎮櫻花的始花期影響較小。

2.2.2 ?清鎮櫻花花期預報模型建立

通過以上分析，從各氣象因子中挑選出與始花期日序數最相關的2個氣象因子，利用SPSS 26.0軟件對始花期日序數與氣象因子進行逐步回歸分析，從而得到了清鎮始花期預報模型：

Y= 79.564-0.114×X1-0.763×X2 （2）

公式（2）中，Y表示始花期日序數，X1表示2月下旬平均氣溫，X2表示2月下旬最高氣溫。方程通過了0.05顯著性檢驗，即回歸方程有效。

2.2.3 ?模型檢驗與預測檢驗

用2013—2020年氣象資料對所建立的花期預報方程進行回測，以驗證模型的準確性。結果表明模型擬合值與實際觀測值差異較小，絕對誤差平均為1.0 d，最大誤差為3 d（見圖1）。利用清鎮櫻花始花期預測模型，將2021年的各預測因子代入得到Y值為65，換算成日期，得出2021年清鎮櫻花始花期為3月6日，與實際觀測日期3月4日相差2 d，在允許誤差范圍內。

2.3 ?櫻花始花期與前期積溫

2.3.1 ?清鎮櫻花花期預報模型建立

大多數植物開花前的積溫會在一定程度上影響花期早晚，而植物不同生長階段所需要的積溫是相對固定的，積溫模型預測法就是利用櫻花始花期所需要的積溫，檢驗是否達到閾值。不同地方、不同植物始花期所需積溫閾值各不相同。利用2013—2020年間清鎮逐日氣象數據，選取下限溫度0、3、5 ℃和日最高氣溫，起算時間為1月1日、1月15日、2月1日至當年始花期日，分別計算從起始日期到始花期的活動積溫值。表2為清鎮櫻花始花期序列的各積溫平均值。方差大小可表示花期的穩定性程度[8]，差值越小，表明花期與平均花期差距越小，結果越穩定。研究發現1月15日最高溫度累積平均值方差最小為542.3。因此采用起算日為1月15日、日最高溫度累積溫度來確定櫻花始花期更為合適。具體指標為1月15日起，當日最高溫度累積值≥611.8 ℃時，這天即為清鎮櫻花開花日。

2.3.2 模型檢驗與預測檢驗

用2013—2020年清鎮逐日氣象數據，將每年1月15日至始花期序列對各積溫模型進行回測，算出歷年櫻花期的擬合日期，將擬合日期減去實際日期的絕對值，得出模型的擬合差異。結果表明（見表3），1月15日至櫻花始花期的0、3、5 ℃有效積溫及最高溫度累積預報的始花期平均絕對誤差分別為2.8 d、3.8 d、2.1 d、1.5 d，其中最高溫度累積模型預報誤差天數最小。從圖2可以看出，其擬合結果與實際結果誤差最大值為3 d，差距為2 d之內的預測準確率達88%。這表明最高溫度積溫預報模型對清鎮櫻花始花期預測有良好的擬合效果。利用清鎮櫻花始花期預測模型，預測2021年櫻花始花期為3月4日，與實際觀測日期3月4日符合，預測效果良好。

2.4 ?特早、特遲開花年份的氣候分析

通過對清鎮櫻花開花特早年（2013年）、特遲年份（2017年）與近30年（1991—2020年）氣象條件分析，發現2013年1—3月的平均氣溫和最高氣溫明顯偏高，尤其是2月的平均氣溫、最高氣溫、日照時數較歷史同期偏高32%、22%、19.7%，降水偏少47.6%，有利于清鎮櫻花始花期提前。分析2017年的氣象資料，發現2 月下旬的平均氣溫、最高氣溫、日照時數較歷史同期偏低35.8%、26.2%、50.3%，降水偏多，濕度偏大，從而導致櫻花花期推遲。

3 ?小結

本文對清鎮櫻花始花期與氣象因子的相關性進行了分析，建立了櫻花始花期預測的線性回歸方程，同時利用積溫方法建立櫻花始花期預測指標，并進行了回測檢驗及2021年始花期預測。

1）清鎮櫻花始花期平均日序為68.6 d，對應日期為3月8—9日，平均盛花期為3月14日，平均落花期為3月26日。清鎮櫻花始花期最早為3月5日，最晚為3月14日，變幅長達9 d。

2）清鎮櫻花始花期與1月中旬、2月下旬平均氣溫，1月下旬日照時數和2月下旬、2月最高氣溫呈顯著負相關關系;利用相關氣象因子建立了回歸模型，模型擬合值與實際觀測值差異較小，平均絕對誤差為1.0 d，最大誤差為3 d。

3）清鎮櫻花始花期積溫模型具體指標為自1月15日起，當日最高溫度累積值≥611.8 ℃時，這天即為清鎮櫻花始花期，利用該指標驗證2013—2020年清鎮櫻花始花期預報平均誤差天數為1.5 d，誤差較小。因此，該指標可在未來清鎮櫻花始花期預報服務中使用。

4）通過對清鎮櫻花開花特早年和特遲年的氣候分析，發現特早年1—3月平均氣溫和最高氣溫明顯偏高，特別是2月氣溫偏高、日照偏多，有利于櫻花始花期提前;特遲年2月下旬氣溫偏低、日照偏少、降雨偏多，導致花期推遲。

櫻花始花期受到多種因素的影響，除氣象因素外，櫻花生長環境也有很大的影響。如果預測花期前一周的氣溫發生劇烈波動，花期不穩定性將大大增加。目前已有的清鎮櫻花物候觀測資料不多，需進一步觀測、收集更多櫻花花期樣本資料，繼續優化預測模型，提高花期預報準確率。

參考文獻：

[1] ? 劉亞辰，王煥炯，戴君虎，等.物候學方法在歷史氣候變化重建中的應用[J].地理研究，2014，33（4）：603-613.

[2] ? 韓路，徐星.谷城櫻花谷櫻花花期預報模型建設探討[J].現代農業科技，2019（12）：173，176.

[3] ? 肖丹，徐衛紅，黃保云，等.眉山市東坡區櫻花盛花日預測[J].南方農業，2017，11（22）：20-23.

[4] ? 陳正洪，肖玫，陳璇.櫻花花期變化特征及其與冬季氣溫變化的關系[J].生態學報，2008（11）：5209-5217.

[5] ? 魏華兵，李磊，羅昱.咸寧市葛仙山景區原生態櫻花群落花期的分析研究[J].湖北農業科學，2020，59（8）：103-106，120.

[6] ? 成青.竺可楨的物候學研究與影響[D].杭州：浙江工業大學，2019.

[7] ? 郭連云，趙年武.貴德縣梨樹始花期與氣象因子的相關分析及預報模型[J].中國農學通報，2016，32（7）：147-151.

[8] ? 饒莉娟，陳清峰，李慶寶.青島大珠山杜鵑花初花期預報方法研究[J].江西農業學報，2019，31（11）：90-97.

（責任編輯：易 ?婧）

收稿日期：2020-09-03;

基金項目：貴州省氣象局項目“清鎮櫻花始花期預報”（黔氣科登〔2020〕08-12）。

作者簡介：楊濤（1996—），男，貴州銅仁人，本科，助理工程師，主要從事綜合氣象服務。E-mail： gzqx9606@163.com。

為通信作者，E-mail： hao.0210@163.com。