清鎮(zhèn)市櫻花始花期預(yù)測分析

楊濤 陳宇 張皓 蔡志穎

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

關(guān)鍵詞 櫻花;花期;氣象因子;預(yù)測模型;貴州省清鎮(zhèn)市

中圖分類號：S162.5+5 文獻(xiàn)標(biāo)志碼：B DOI：10.19415/j.cnki.1673-890x.2021.31.043

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

本文研究清鎮(zhèn)氣象條件對櫻花始花期的影響，以期為當(dāng)?shù)赜嘘P(guān)部門提供服務(wù)決策依據(jù)，指導(dǎo)人們合理安排時間觀賞櫻花，助力鄉(xiāng)村旅游發(fā)展。

1 ?資料與方法

1.1 ?研究區(qū)概況

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

1.2 ?數(shù)據(jù)來源

氣象數(shù)據(jù)來源于貴州省信息中心，選用清鎮(zhèn)市2013—2020年逐日平均氣溫、最高氣溫、降水、日照時數(shù)、平均風(fēng)速等氣象數(shù)據(jù)，并采取調(diào)查及專業(yè)觀測方法，參照《農(nóng)業(yè)氣象觀測規(guī)范》獲取了清鎮(zhèn)櫻花2013—2020年花期的相關(guān)資料。

1.3 ?研究方法

1.3.1 ?數(shù)據(jù)處理

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

1.3.2 ?積溫計(jì)算方法

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

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

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

1.3.3 ?統(tǒng)計(jì)分析

應(yīng)用SPSS26.0、Excel軟件對數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，用相關(guān)分析法分析清鎮(zhèn)櫻花始花期與氣溫、日照、風(fēng)速等氣象因子的相關(guān)性，利用逐步回歸分析法找出對櫻花始花期有顯著影響的關(guān)鍵氣象因子并建立預(yù)報模型，其中2013—2020年的資料用于模型建立與回測，2021年的資料則作為獨(dú)立樣本進(jìn)行了預(yù)測檢驗(yàn)。

2 ?結(jié)果與分析

2.1 ?清鎮(zhèn)櫻花花期特征分析

根據(jù)2013—2020年清鎮(zhèn)櫻花的物候資料，發(fā)現(xiàn)清鎮(zhèn)櫻花始花期平均日序?yàn)?8.6 d，對應(yīng)日期為3月8—9日，平均盛花期為3月14日，平均落花期為3月26日，平均花期可達(dá)17.6 d。清鎮(zhèn)櫻花始花期最早為2013年3月5日，最晚為2017年3月14日，變幅長達(dá)9 d。借助均方差的大小判定始花期的穩(wěn)定性程度[7]，發(fā)現(xiàn)櫻花始花期日序數(shù)均方差為3.2，超過63%的年份均偏離平均始花期2 d以上，這表明清鎮(zhèn)櫻花始花期近8年來并不穩(wěn)定。

2.2 ?櫻花始花期與氣象因子的關(guān)系

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

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

2.2.2 ?清鎮(zhèn)櫻花花期預(yù)報模型建立

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

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

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

2.2.3 ?模型檢驗(yàn)與預(yù)測檢驗(yàn)

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

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

2.3.1 ?清鎮(zhèn)櫻花花期預(yù)報模型建立

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

2.3.2 模型檢驗(yàn)與預(yù)測檢驗(yàn)

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

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

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

3 ?小結(jié)

本文對清鎮(zhèn)櫻花始花期與氣象因子的相關(guān)性進(jìn)行了分析，建立了櫻花始花期預(yù)測的線性回歸方程，同時利用積溫方法建立櫻花始花期預(yù)測指標(biāo)，并進(jìn)行了回測檢驗(yàn)及2021年始花期預(yù)測。

1）清鎮(zhèn)櫻花始花期平均日序?yàn)?8.6 d，對應(yīng)日期為3月8—9日，平均盛花期為3月14日，平均落花期為3月26日。清鎮(zhèn)櫻花始花期最早為3月5日，最晚為3月14日，變幅長達(dá)9 d。

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

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

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

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

參考文獻(xiàn)：

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

[2] ? 韓路，徐星.谷城櫻花谷櫻花花期預(yù)報模型建設(shè)探討[J].現(xiàn)代農(nóng)業(yè)科技，2019（12）：173，176.

[3] ? 肖丹，徐衛(wèi)紅，黃保云，等.眉山市東坡區(qū)櫻花盛花日預(yù)測[J].南方農(nóng)業(yè)，2017，11（22）：20-23.

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

[5] ? 魏華兵，李磊，羅昱.咸寧市葛仙山景區(qū)原生態(tài)櫻花群落花期的分析研究[J].湖北農(nóng)業(yè)科學(xué)，2020，59（8）：103-106，120.

[6] ? 成青.竺可楨的物候?qū)W研究與影響[D].杭州：浙江工業(yè)大學(xué)，2019.

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

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

（責(zé)任編輯：易 ?婧）

收稿日期：2020-09-03;

基金項(xiàng)目：貴州省氣象局項(xiàng)目“清鎮(zhèn)櫻花始花期預(yù)報”（黔氣科登〔2020〕08-12）。

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

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