基于SW物候模型“牛媽媽棗”收獲期預(yù)測

劉振宇 曹曉霞 史有瑜 陳秀影

[摘 要]基于唐山豐潤區(qū)國家一般氣象站2013年1月～2017年8月逐日平均氣溫資料和2017年8月逐日平均氣溫滾動預(yù)報資料，應(yīng)用SW物候模型，對2017年豐潤山緣景區(qū)“牛媽媽棗”收獲（采摘）期進行預(yù)測，取得了較好的效果。研究結(jié)果表明：應(yīng)用SW物候模型預(yù)測“牛媽媽棗”收獲（采摘）期，預(yù)測結(jié)果與實況接近，該模型適用性較高，可進行進一步的業(yè)務(wù)試用。

[關(guān)鍵詞]物候模型;牛媽媽棗;收獲期;預(yù)測

[中圖分類號]S685.99 [文獻標識碼]A

“牛媽媽棗”是20世紀50年代在河北省唐山市豐潤區(qū)王官營鎮(zhèn)下水路村發(fā)現(xiàn)的一種優(yōu)良鮮食棗單株，因其所結(jié)果實呈長圓錐形，胴部具有隘痕，似牛乳頭而被命名為“牛媽媽棗”。“牛媽媽棗”果實皮薄肉細，質(zhì)地細脆，汁液多，酸甜適口，2003年9月，經(jīng)河北科技師范學(xué)院、河北省農(nóng)林科學(xué)院、唐山市林業(yè)局專家實地考察和品嘗鑒定，河北農(nóng)業(yè)大學(xué)中國棗研究中心對鮮果分析測定，認為品質(zhì)上等;2005年唐山市豐潤區(qū)政府將其定義為豐潤區(qū)的特色果品，并注冊了“峪泉”商標;2016年，豐潤“牛媽媽棗”通過了國家有機食品認證;2017年，豐潤山緣景區(qū)成功申報國家地理標志認證，確定山緣景區(qū)為“牛媽媽棗”的原產(chǎn)地。

豐潤區(qū)隸屬于河北省唐山市，地處河北省東部，北依燕山，中南部偎冀東平原。豐潤山緣景區(qū)為地勢起伏不大的山地，林果生長旺盛，現(xiàn)已建成集生態(tài)特種養(yǎng)殖、旅游觀光、休閑度假、果蔬采摘、鄉(xiāng)村體驗、苗木栽培、科技農(nóng)業(yè)研發(fā)推廣為一體的現(xiàn)代休閑農(nóng)業(yè)示范園區(qū)。景區(qū)自2016年開始承辦豐潤區(qū)大棗文化節(jié)，為保證大棗文化節(jié)成功舉辦，準確預(yù)報“牛媽媽棗”適宜收獲期，并以此為參考提前確定文化節(jié)舉辦日期至關(guān)重要，對提高文化節(jié)期間企業(yè)的經(jīng)濟效益和宣傳效果有很大幫助。本文基于SW（Spring Warming）物候模型利用有效積溫指標對豐潤山緣景區(qū)“牛媽媽棗”2013-2016年收獲期（本文指始采摘期，下同）的觀測資料進行模擬檢驗，并對2017年收獲期進行試預(yù)報，以期推動和提高“牛媽媽棗”采摘期預(yù)報服務(wù)水平，從而更好地為園區(qū)提供服務(wù)。

1 研究資料與方法

1.1 研究資料

以唐山豐潤山緣景區(qū)“牛媽媽棗”為研究對象，物候資料為該植物2013～2017年的成熟可采期觀測資料;氣象資料取自河北省唐山市豐潤區(qū)氣象局國家一般氣象站2013年1月1日～2017年8月31日逐日平均氣溫和最高氣溫觀測資料，觀測資料經(jīng)過了嚴格的質(zhì)量控制。

1.2 研究方法

本文采用SW物候模型對“牛媽媽棗”收獲期進行預(yù)測，積溫使用有效積溫。有效積溫預(yù)測模型如下：

式中，y是預(yù)測的“牛媽媽棗”收獲期（日序），xt是第t天的日平均氣溫;R（xt）是高于某一界限溫度值的有效溫度，Tb是界限溫度，本文依據(jù)從唐山市豐潤區(qū)林業(yè)局獲得的“牛媽媽棗”生長發(fā)育指標，將Tb取13.0℃;F為有效積溫閾值，本文統(tǒng)計值為1340℃;t0是有效積溫開始累積的時間，以日序來表示，起始日期取每年日平均氣溫穩(wěn)定通過13.0℃的初日。由于高溫天氣對“牛媽媽棗”的生長發(fā)育影響很大，本文在計算累積有效溫度時，將連續(xù)2天以上高于35. 0℃日期的有效溫度扣除。

當收獲期模擬日期與實際日期差值≤1d評定為正確，2～3為基本正確，≥4d為不正確。

2 結(jié)果與分析

對唐山市豐潤區(qū)國家一般氣象站2013～2016年逐日平均氣溫觀測資料進行5日滑動平均，得到各年日平均氣溫穩(wěn)定通過13.0℃的初日，如表1所示。

2.1 2013～2016年收獲期模擬檢驗

以日平均氣溫穩(wěn)定通過13.0℃的初日為起始日期，用SW物候模型計算有效積溫，得到2013～2016年“牛媽媽棗”達到有效積溫閾值的日期，其模擬日期與觀測日期如表2所示。

可見，預(yù)測模型對收獲期擬合效果較好，4年均達到正確和基本正確的評價標準，擬合值與實際值的相關(guān)系數(shù)為0.9825，通過0.05的顯著性檢驗。

2.2 2017年收獲期預(yù)報及檢驗

目前業(yè)務(wù)中可獲得時效達10天的ECMWF細網(wǎng)格預(yù)報產(chǎn)品，實際業(yè)務(wù)中通常做提前一周的預(yù)報，本文以10天預(yù)報時效進行檢驗。進行收獲期預(yù)報的具體做法是從8月1日起每天進行溫度預(yù)報，預(yù)報日的日平均氣溫用預(yù)報的日最高氣溫和日最低氣溫的平均值代替;如預(yù)報連續(xù)兩日以上最高氣溫大于等于35℃時，則將該連續(xù)日平均溫度置為0。首先計算出當年穩(wěn)定通過13℃的初日，再計算初日至預(yù)報起點日后的有效積溫，當累積有效積溫達到閾值時，即確定該日為預(yù)測的收獲日（始采摘日）。

用唐山市豐潤區(qū)國家一般氣象站2017年的逐日平均氣溫實況值進行統(tǒng)計，2017年穩(wěn)定通過13℃的初日為4月2日，2017年4月2日至8月1日期間有效積溫為1091.6℃，有效積溫達到或超過閾值（1340℃）取決于預(yù)報日之后逐日平均氣溫的預(yù)報結(jié)果。將2017年8月1日、8月8日、8月9日、8月12日對未來10天的預(yù)報結(jié)果分別列于表3和表4。從表3可以看出，根據(jù)8月1日的預(yù)報結(jié)果，8月11日有效積溫值為1232.1℃，距離閾值還尚需些時日;8月8日滾動預(yù)報的8月18日有效積溫值為1328.9℃，已接近閾值。繼續(xù)滾動訂正預(yù)報結(jié)果，從表4可以看出，根據(jù)8月9日預(yù)報結(jié)果，8月19日有效積溫值為1340.9℃，剛好達到閾值;繼續(xù)滾動預(yù)報加以印證，8月10日預(yù)報（表略）8月19日有效積溫為1339.9℃，8月20日有效積溫為1351.4℃;8月11日預(yù)報（表略）8月19日有效積溫為1342.4℃;8月12日預(yù)報（表4）8月20日有效積溫為1344.5℃。因此，在8月9日預(yù)測的大約8月19日進入收獲（始采摘）期，根據(jù)后續(xù)幾天滾動預(yù)報結(jié)果，結(jié)論也基本一致。

2017年實況是，物候觀測8月19日進入采摘期，氣象實況資料8月20日累積有效積溫達到1344.8℃，與預(yù)測結(jié)果基本一致。

3 討論

本文采用SW物候模型，在前期氣溫實際觀測資料以及后期逐日氣溫滾動預(yù)報的基礎(chǔ)上，進行了2017年“牛媽媽棗”收獲期的預(yù)測，為“牛媽媽棗”適時采摘和園區(qū)“采摘文化節(jié)”的籌備工作提供了科學(xué)的參考依據(jù)，收到了較好的效果。采用一周（7日）滾動預(yù)報結(jié)果所預(yù)測的始采摘日也在19-20日，具體為8月11日可預(yù)測19日累積有效積溫可達閾值，8月12-15日均可預(yù)測8月20日累積有效積溫達到閾值。實際運用中可根據(jù)服務(wù)對象的需求進行提前10天或提前一周的預(yù)測服務(wù)，通常情況下提前一周的效果會更好。

鑒于目前用于檢驗的樣本較少，本文的研究結(jié)果還不夠全面，之后將隨著樣本的增多進行更深入的檢驗和業(yè)務(wù)應(yīng)用。

[參考文獻]

[1] 翟君芳.牛媽媽棗生物學(xué)特性觀察[J].河北果樹，2011（05）.

[2] 王連芹.優(yōu)良鮮食棗品系——“牛媽媽棗”[J].現(xiàn)代農(nóng)村科技，2009（15）.

[3] 戴君虎，陶澤興，王煥炯，等.物候模型在氣候定量重建中的應(yīng)用[J].第四紀研究，2016（03）.

[4] 董海濤，譚麗靜，劉洪林，左曉強，于文革，孟鑫.基于積溫模型丹東地區(qū)桃樹盛花期預(yù)測研究[J].氣象與環(huán)境學(xué)報，2018（01）.

[5] 張愛英，張建華，高迎新，等.SW物候模型在北京櫻花始花期預(yù)測中的應(yīng)用[J].氣象科技，2015（02）.