謝佰承 郭凌曜 杜東升 譚 儼 王國(guó)棟
(1.西北農(nóng)林科技大學(xué)理學(xué)院 楊凌 712100; 2.氣象防災(zāi)減災(zāi)湖南省重點(diǎn)實(shí)驗(yàn)室 長(zhǎng)沙 410118; 3.湖南省氣候中心 長(zhǎng)沙 410118)
油茶(Camelliaoleifera)屬山茶科(Theaceae)山茶屬(Camellia)常綠小喬木,為世界四大木本油料樹種之一,具有生長(zhǎng)快、結(jié)果早、適應(yīng)性強(qiáng)、用途廣等特點(diǎn)(莊端林, 2008)。目前,我國(guó)油茶種植面積超過430萬hm2,主要分布在長(zhǎng)江流域及其以南的14個(gè)省份,油茶籽年產(chǎn)量約210萬t,年產(chǎn)茶油約52萬t。其中,湖南油茶種植面積已超過140萬hm2,茶油產(chǎn)量和產(chǎn)值均居全國(guó)首位,是當(dāng)?shù)亟?jīng)濟(jì)發(fā)展的優(yōu)勢(shì)產(chǎn)業(yè)。
長(zhǎng)期以來,我國(guó)油茶存在低產(chǎn)低效問題,除品種、管理水平和社會(huì)因素外,氣象災(zāi)害也是一個(gè)重要制約因素。在氣象因子對(duì)油茶產(chǎn)量的影響方面已有大量研究,如康志雄等(1993)在連續(xù)11年調(diào)查516株固定樣株產(chǎn)量的基礎(chǔ)上,運(yùn)用灰色關(guān)聯(lián)度方法分析17個(gè)氣象因子的作用,結(jié)果發(fā)現(xiàn)上年度11月的日照量、7—8月的蒸發(fā)量、1月的平均氣溫和年日照量等對(duì)油茶產(chǎn)量影響很大; 左繼林等(2010)研究得出,油茶優(yōu)良無性系的干出籽率、鮮果含油率和種仁含油率等經(jīng)濟(jì)性狀指標(biāo)受年內(nèi)不同時(shí)段以及不同年份的氣象因子影響較大,其中最低氣溫、日降雨量、風(fēng)速和日均溫等因子對(duì)油茶鮮出籽率影響很大; 曾燕如等(2010)研究表明,油茶花期的降雨和霜凍等氣象條件會(huì)影響坐果率,進(jìn)而影響產(chǎn)量; 夏尚光等(2011)黃拯等(2017)研究認(rèn)為,特定環(huán)境因子會(huì)影響凈光合速率和蒸騰速率,從而影響油茶生長(zhǎng)發(fā)育; 彭嘉棟等(2016)統(tǒng)計(jì)分析不同物候期氣象因子對(duì)油茶產(chǎn)量的影響,并建立了產(chǎn)量預(yù)測(cè)模型; 丁少凈等(2017)探討干旱脅迫對(duì)油茶花苞生長(zhǎng)及油茶籽和茶油產(chǎn)量的影響,結(jié)果發(fā)現(xiàn)油脂轉(zhuǎn)化期的干旱脅迫影響較大,可使茶油產(chǎn)量下降45.73%,但采用覆蓋和松土措施后僅下降24.14%和27.40%,表明相應(yīng)的管理措施可一定程度緩解干旱危害。
目前,研究主要集中在不同油茶生長(zhǎng)時(shí)期的氣溫、日照、降水量、降水日數(shù)等氣象指標(biāo)與油茶產(chǎn)量的關(guān)系方面,(王道藩, 1983; 郭文楊等,1987; 韋金霖, 2013; 周定杰等, 2013; 袁昌選等, 2014; 鐘飛霞等, 2015; 左繼林, 2014; 宋英強(qiáng)等, 2015; 許彥明,2016)高溫干旱是影響油茶產(chǎn)量和品質(zhì)的重要?dú)庀鬄?zāi)害,如2003年江西省的歷史罕見夏秋高溫干旱使油茶總產(chǎn)量較前3年均值減少12%; 2013年湖南省的夏秋季節(jié)性干旱導(dǎo)致全省油茶產(chǎn)量減少3成以上。本研究探討氣象因子對(duì)油茶產(chǎn)量的影響,確定油茶關(guān)鍵生長(zhǎng)時(shí)期的高溫干旱指標(biāo)及其影響關(guān)系,以期為油茶高溫干旱災(zāi)害提供參考。
油茶產(chǎn)量源于湖南省林業(yè)科學(xué)院提供的2007—2015年省內(nèi)18個(gè)縣(市、區(qū))油茶樣地?cái)?shù)據(jù),基于油茶樣地鮮果采摘量計(jì)算的鮮果單產(chǎn)數(shù)據(jù)。為剔除油茶樹齡、品種和土壤等因素差異的影響,數(shù)據(jù)分析時(shí)首先對(duì)油茶產(chǎn)量數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理: 1) 因良種林和非良種林(指其他未審定的品種)油茶產(chǎn)量差異較大,故剔除非良種林油茶產(chǎn)量數(shù)據(jù),統(tǒng)一選擇良種林油茶產(chǎn)量數(shù)據(jù); 2) 油茶種植后進(jìn)入豐產(chǎn)期(指油茶林進(jìn)入豐產(chǎn)、穩(wěn)產(chǎn)的階段,一般指油茶良種造林后第8年)需7~8年時(shí)間,進(jìn)入豐產(chǎn)期前產(chǎn)量隨樹齡增大而增加,故剔除未進(jìn)入豐產(chǎn)期的油茶產(chǎn)量數(shù)據(jù); 3) 不同管理措施使油茶產(chǎn)量差異較大,高產(chǎn)品種年均產(chǎn)量達(dá)13 395 kg·hm-2,而低產(chǎn)良種僅1 455 kg·hm-2,故對(duì)每塊良種林樣地的油茶產(chǎn)量進(jìn)行標(biāo)準(zhǔn)化處理:
(1)
式中:Xi為單塊樣地某年的產(chǎn)量;X為單塊樣地進(jìn)入豐產(chǎn)期后的多年平均產(chǎn)量;σ為進(jìn)入豐產(chǎn)期后產(chǎn)量的標(biāo)準(zhǔn)差;i為年序; XSi為某年的標(biāo)準(zhǔn)化產(chǎn)量。
剔除低產(chǎn)林、未進(jìn)入豐產(chǎn)期的樣地?cái)?shù)據(jù)后,共有48組標(biāo)準(zhǔn)化油茶產(chǎn)量數(shù)據(jù),排除了油茶樹齡、品種和土壤條件的干擾,可認(rèn)為油茶產(chǎn)量?jī)H受氣象因素影響。油茶標(biāo)準(zhǔn)化產(chǎn)量數(shù)據(jù)見表1。
表1 2007—2015年油茶標(biāo)準(zhǔn)化產(chǎn)量Tab.1 Standardized production of C. oleifera from 2007 to 2015
氣象數(shù)據(jù)來自湖南省氣象信息中心檔案館,采用距油茶樣地最近的氣象站的逐日數(shù)據(jù)。不同物候期主要包括果實(shí)膨大高峰期(6月1日—8月5日)和油脂轉(zhuǎn)化期(8月11日—10月10日),物候期的氣象因子數(shù)據(jù)見表2; 不同物候期≥30 ℃熱積溫和≥30 ℃高溫日數(shù)數(shù)據(jù)見表3。
表2 油茶不同物候期氣象因子統(tǒng)計(jì)Tab.2 Statistics of meteorological factors of C. oleifera in different phenological periods
表3 油茶不同物候期≥30 ℃熱積溫和≥30 ℃高溫日數(shù)Tab.3 The heat accumulated temperature ≥30 ℃ and days of high temperature ≥30 ℃ of C. oleifera in different phenological periods
油茶雖然耐旱,但果實(shí)膨大高峰期(6月1日—8月5日)和油脂轉(zhuǎn)化期(8月11日—10月10日)的持續(xù)高溫對(duì)油茶產(chǎn)量和品質(zhì)具有影響。本研究以日均氣溫≥30 ℃作為高溫指標(biāo),當(dāng)日均氣溫≤30 ℃時(shí),則日熱積溫為0:
∑Td=∑[Km(Tm-T0)]。
(2)
式中:Td為日熱積溫(℃·d-1);Tm為日均最高氣溫(℃);T0為影響油茶產(chǎn)量和品質(zhì)的臨界溫度(30 ℃);Km=1,為訂正系數(shù)。
統(tǒng)計(jì)各油茶樣地不同生長(zhǎng)階段的≥30 ℃熱積溫和≥30 ℃高溫日數(shù),用于分析果實(shí)膨大高峰期和油脂轉(zhuǎn)化期的氣象因素對(duì)油茶產(chǎn)量的影響。采用相關(guān)分析研究氣象因子與油茶產(chǎn)量的關(guān)系,應(yīng)用逐步回歸建立油茶氣象產(chǎn)量預(yù)測(cè)模型,利用SPSS軟件分析。
基于2007—2013年的產(chǎn)量和氣象數(shù)據(jù)建立產(chǎn)量模型,基于2014—2015年的數(shù)據(jù)驗(yàn)證模型精度:
(3)
2.1.1 果實(shí)膨大高峰期氣象因子影響 由圖1可知,油茶標(biāo)準(zhǔn)化產(chǎn)量與平均氣溫呈顯著負(fù)相關(guān)(P<0.05,-0.433 9),與降水日數(shù)呈顯著正相關(guān)(P<0.05,0.591 5),與降水量呈正相關(guān),與日照時(shí)數(shù)呈負(fù)相關(guān)。這表明,果實(shí)膨大高峰期降水量和降水日數(shù)較多利于油茶產(chǎn)量提高。油茶生長(zhǎng)除需要一定降水量外,果實(shí)膨大高峰期降水的合理均衡分布也是油茶正常生理活動(dòng)和油脂合成轉(zhuǎn)化的重要條件。
圖1 果實(shí)膨大高峰期氣象因子與油茶標(biāo)準(zhǔn)化產(chǎn)量的關(guān)系Fig. 1 The relationship between meteorological factors and the standardized yield in the peak period of fruit expansion
2.1.2 油脂轉(zhuǎn)化期氣象因子的影響 由圖2可知,油茶標(biāo)準(zhǔn)化產(chǎn)量與平均氣溫、降水量、降水日數(shù)呈負(fù)相關(guān),與日照時(shí)數(shù)呈正相關(guān)。雖然油脂轉(zhuǎn)化期高溫干旱不利于果實(shí)增長(zhǎng),但充足的日照能保障油茶光合作用,利于果實(shí)增長(zhǎng)和油脂積累,因此連續(xù)降水會(huì)限制油茶產(chǎn)量增加。
圖2 油脂轉(zhuǎn)化期氣象因子與油茶標(biāo)準(zhǔn)化產(chǎn)量的關(guān)系Fig. 2 The relationship between meteorological factors and the standardized yield in oil transformation period
2.2.1 果實(shí)膨大高峰期≥30 ℃熱積溫和≥30 ℃高溫日數(shù)的影響 由圖3A可知,油茶標(biāo)準(zhǔn)化產(chǎn)量與≥30 ℃熱積溫呈顯著負(fù)相關(guān)(P<0.05),其統(tǒng)計(jì)關(guān)系為y=-0.040 5x+1.486 6(r=0.759 8,P<0.05),當(dāng)≥30 ℃熱積溫達(dá)到36.7 ℃時(shí)油茶標(biāo)準(zhǔn)產(chǎn)量為0,表明若≥30 ℃熱積溫繼續(xù)增加將會(huì)抑制油茶果實(shí)生長(zhǎng),使標(biāo)準(zhǔn)化產(chǎn)量為負(fù),即實(shí)際產(chǎn)量在低于多年平均產(chǎn)量的范圍內(nèi)繼續(xù)降低。油茶雖為耐旱樹種,但果實(shí)膨大高峰期持續(xù)高溫和缺少降水會(huì)影響產(chǎn)量,難獲豐產(chǎn)。
圖3 果實(shí)膨大高峰期日熱積溫及高溫日數(shù)與標(biāo)準(zhǔn)產(chǎn)量的關(guān)系Fig. 3 The relationship between the standardized yield and the daily accumulated temperature and days of high temperature in the peak period of fruit expansion period
由圖3B可知,油茶標(biāo)準(zhǔn)化產(chǎn)量與的≥30 ℃高溫日數(shù)呈顯著負(fù)相關(guān)(P<0.05),其統(tǒng)計(jì)關(guān)系為y=-0.071 7x+1.885 9(r=0.650 9,P<0.05)。當(dāng)≥30 ℃高溫日數(shù)達(dá)到26天時(shí)油茶標(biāo)準(zhǔn)化產(chǎn)量為0,表明若≥30 ℃高溫日數(shù)繼續(xù)增加將會(huì)抑制油茶果實(shí)生長(zhǎng),使實(shí)際產(chǎn)量在低于多年平均產(chǎn)量的范圍內(nèi)繼續(xù)降低。因此,當(dāng)≥30 ℃高溫日數(shù)超過26天后,需采取防災(zāi)減災(zāi)措施,如在油茶根部灌溉補(bǔ)水,以減輕水分虧缺并降低災(zāi)害減產(chǎn)。
從油茶標(biāo)準(zhǔn)化產(chǎn)量與≥30 ℃熱積溫和≥30 ℃高溫日數(shù)的相關(guān)系數(shù)可知,油茶標(biāo)準(zhǔn)化產(chǎn)量對(duì)果實(shí)膨大高峰期≥30 ℃熱積溫的指標(biāo)更敏感。
2.2.2 油脂轉(zhuǎn)化期≥30 ℃熱積溫和≥30 ℃高溫日數(shù)的影響 由圖4A可知,油茶標(biāo)準(zhǔn)化產(chǎn)量與≥30 ℃熱積溫呈顯著負(fù)相關(guān)(P<0.05),其統(tǒng)計(jì)關(guān)系為y=-0.068 8x+0.936 9(r=0.725 1,P<0.05),當(dāng)≥30 ℃熱積溫達(dá)到13.6 ℃時(shí)油茶≥30 ℃標(biāo)準(zhǔn)化產(chǎn)量為0,表明若≥30 ℃熱積溫繼續(xù)增加將會(huì)使油茶標(biāo)準(zhǔn)產(chǎn)量化為負(fù),即實(shí)際產(chǎn)量在低于多年平均產(chǎn)量的范圍內(nèi)繼續(xù)降低。
由圖4B可知,油茶標(biāo)準(zhǔn)化產(chǎn)量與≥30 ℃高溫日數(shù)呈顯著負(fù)相關(guān)(P<0.05),其統(tǒng)計(jì)關(guān)系為y=-0.1270x+1.226 7(r=0.613 0,P<0.05),當(dāng)≥30 ℃高溫日數(shù)達(dá)到10天時(shí)油茶標(biāo)準(zhǔn)化產(chǎn)量為0,表明若≥30 ℃高溫日數(shù)繼續(xù)增加將會(huì)使油茶標(biāo)準(zhǔn)化產(chǎn)量為負(fù),即持續(xù)高溫不利油茶鮮果產(chǎn)量形成,實(shí)際產(chǎn)量會(huì)在低于多年平均產(chǎn)量的范圍內(nèi)繼續(xù)降低。這主要是由于高溫干旱不利于油茶光合作用,造成油茶果的果徑、果高、果形指數(shù)和體積等參數(shù)降低(王瑞輝等, 2014)。
圖4 油脂轉(zhuǎn)化期熱積溫和高溫日數(shù)與標(biāo)準(zhǔn)化產(chǎn)量的關(guān)系Fig. 4 The relationship between standardized yield and the daily accumulated temperature and days of high temperature in oil transformation period
與油茶果實(shí)膨大高峰期(6月1日—8月5日)相比,油脂轉(zhuǎn)化期(8月11日—10月10日)影響產(chǎn)量的≥30 ℃熱積溫和≥30 ℃高溫日數(shù)的閾值均較小,主要是因?yàn)榇藭r(shí)已逐漸進(jìn)入秋季。與果實(shí)膨大高峰期類似,油茶標(biāo)準(zhǔn)化產(chǎn)量對(duì)≥30 ℃熱積溫的敏感性大于對(duì)高溫日數(shù)的敏感性。
根據(jù)上述結(jié)果,選用果實(shí)膨大高峰期和油脂轉(zhuǎn)化期及果實(shí)膨大高峰期-油脂轉(zhuǎn)化期≥30 ℃熱積溫、≥30 ℃高溫日數(shù)建立油茶標(biāo)準(zhǔn)化產(chǎn)量預(yù)測(cè)模型(表4),分析可知,在果實(shí)膨大高峰期-油脂轉(zhuǎn)化期建立的回歸模型擬合效果較好,Y=-0.024x1-0.012X2+1.611,R2=0.629,達(dá)到最優(yōu),根據(jù)方程可知,≥30 ℃熱積溫和≥30 ℃高溫日數(shù)與油茶標(biāo)準(zhǔn)化產(chǎn)量均呈負(fù)相關(guān)關(guān)系。
表4 不同物候期油茶標(biāo)準(zhǔn)化產(chǎn)量預(yù)測(cè)模型①Tab.4 Standardized yield prediction model of C. oleifera in different phenological periods
由圖5可知,基于2014—2015年湖南省24個(gè)油茶種植區(qū)標(biāo)準(zhǔn)化產(chǎn)量進(jìn)行檢驗(yàn)時(shí),擬合產(chǎn)量與實(shí)際產(chǎn)量相關(guān)性較好,其r=0.793 0, RMSE=0.775 3,P<0.05。因此,選用果實(shí)膨大高峰期-油脂轉(zhuǎn)化期建立的預(yù)測(cè)模型,結(jié)合長(zhǎng)期氣象趨勢(shì),可進(jìn)行當(dāng)年或來年油茶產(chǎn)量預(yù)測(cè)。
圖5 湖南省不同站點(diǎn)油茶實(shí)際產(chǎn)量與擬合產(chǎn)量Fig. 5 The actual yield and fitting yield of C. oleifera at different sites in Hunan Province
高溫干旱是油茶林的主要災(zāi)害,油茶林生長(zhǎng)初期,樹齡越小,抗旱能力越弱,高溫干旱災(zāi)害越嚴(yán)重。1年生新造油茶林平均死亡率達(dá) 52.5%,正常株僅13.4%,嚴(yán)重時(shí)林木全部死亡(陳永忠等, 2014)。本研究首次采用≥30 ℃熱積溫和≥30 ℃高溫日數(shù)作為關(guān)鍵氣象因子指標(biāo)對(duì)油茶關(guān)鍵生長(zhǎng)時(shí)期(果實(shí)膨大高峰期、油脂轉(zhuǎn)化期)的高溫干旱災(zāi)害進(jìn)行提前預(yù)測(cè),通過統(tǒng)計(jì)分析,明確給出油茶果實(shí)膨大高峰期≥30 ℃熱積溫達(dá)到36.7 ℃或≥30 ℃高溫日數(shù)達(dá)到26天,以及油脂轉(zhuǎn)化期≥30 ℃熱積溫達(dá)到13.6 ℃或≥30 ℃高溫日數(shù)達(dá)到10天,為高溫干旱出現(xiàn)的閾值,可以指導(dǎo)油茶生產(chǎn)中采取相應(yīng)抗旱措施,為穩(wěn)產(chǎn)高產(chǎn)提供保障。
以往研究多采用油茶樹外觀形態(tài)和生理指標(biāo)指示油茶干旱脅迫,如葉卷曲變褐、頂梢枯萎、老葉片脫落、枝條枯萎等(曹永慶等, 2014); 成熟葉片變薄變小,氣孔面積、氣孔周長(zhǎng)、氣孔開度減小,主脈厚度增加等(曹林青等, 2018); 果徑、果高、果形指數(shù)和體積等參數(shù)變小等(王瑞輝等, 2014)。僅采用外觀形態(tài)和生理指標(biāo)指示干旱脅迫時(shí),因油茶生長(zhǎng)周期已基本完成,很難通過其他管理或補(bǔ)救措施避免油茶產(chǎn)量降低,因此,采用關(guān)鍵生長(zhǎng)時(shí)期氣象指標(biāo)進(jìn)行提前診斷,可提前做好災(zāi)害預(yù)防,以確保油茶產(chǎn)量形成。
高溫干旱一直是影響南方油茶產(chǎn)量的關(guān)鍵因子,隨著氣候變暖,高溫干旱出現(xiàn)頻次逐年增加,7—9月是果實(shí)增長(zhǎng)和油脂形成的關(guān)鍵時(shí)期,持續(xù)高溫、干旱缺水會(huì)限制油茶果實(shí)生長(zhǎng)和油脂轉(zhuǎn)化,不利于果實(shí)形成或形成飽滿果實(shí),不利于優(yōu)質(zhì)豐產(chǎn),含油量低、品質(zhì)差(王瑞輝等, 2014)。在果實(shí)膨大高峰期和油脂轉(zhuǎn)化期,隨著≥30 ℃高溫日數(shù)增加,油茶產(chǎn)量減少,在前人研究中也得到印證(蔣元華等, 2015)。
此外,以往國(guó)內(nèi)學(xué)者大多選用油茶產(chǎn)量與整個(gè)物候期氣象因子(氣溫、降水、日照時(shí)數(shù))進(jìn)行逐步相關(guān)分析,篩選相關(guān)因子建立油茶產(chǎn)量回歸模型(彭嘉棟等, 2016),結(jié)果發(fā)現(xiàn)大多數(shù)常規(guī)氣象因子與油茶產(chǎn)量的相關(guān)性不強(qiáng)。本研究主要考慮油茶果實(shí)膨大高峰期和油脂轉(zhuǎn)化期≥30 ℃熱積溫和≥30 ℃高溫日數(shù)作為油茶生長(zhǎng)的脅迫氣象因子,建立油茶產(chǎn)量預(yù)測(cè)模型,適用性較強(qiáng),便于運(yùn)用,可為油茶產(chǎn)量預(yù)測(cè)提供借鑒參考。隨著不斷增多的極端氣候事件及其對(duì)油茶產(chǎn)量的影響,需要更加精細(xì)地預(yù)測(cè)油茶產(chǎn)量,這將成為油茶高產(chǎn)穩(wěn)產(chǎn)的氣象保障服務(wù)的重要內(nèi)容,可有效提升油茶防災(zāi)減災(zāi)的科學(xué)性。
當(dāng)前,人工管理和撫育對(duì)油茶產(chǎn)量有較大影響,因此油茶產(chǎn)量預(yù)測(cè)模型的準(zhǔn)確度十分關(guān)鍵,多年生油茶產(chǎn)量與氣候相關(guān)更為緊密。積溫是某一時(shí)段內(nèi)逐日平均氣溫之和,是反映作物發(fā)育對(duì)熱量的要求和地區(qū)熱量資源的評(píng)價(jià)指標(biāo),是影響植物生長(zhǎng)的重要因素之一,本研究?jī)H采用≥30 ℃熱積溫作為指標(biāo),考慮果實(shí)膨大高峰期和油脂轉(zhuǎn)化期熱害累積,如采用≥30 ℃小時(shí)當(dāng)量熱積溫建立模型,或許更能反映出油茶生長(zhǎng)受高溫?zé)岷Φ膹?qiáng)度,建立的預(yù)測(cè)模型更加趨于穩(wěn)定和準(zhǔn)確。
果實(shí)膨大高峰期和油脂轉(zhuǎn)化期是影響油茶產(chǎn)量和品質(zhì)的關(guān)鍵時(shí)期,油茶產(chǎn)量與≥30 ℃熱積溫和≥30 ℃高溫日數(shù)均呈顯著負(fù)相關(guān),表明這2個(gè)關(guān)鍵物候期的高溫干旱會(huì)抑制果實(shí)生長(zhǎng),造成油茶產(chǎn)量減產(chǎn)。導(dǎo)致油茶產(chǎn)量低于多年均值的≥30 ℃熱積溫閾值在果實(shí)膨大高峰期為36.7 ℃,在油脂轉(zhuǎn)化期為13.6 ℃,相應(yīng)的≥30 ℃高溫日數(shù)閾值為26天和10天,油茶減產(chǎn)對(duì)果實(shí)膨大高峰期≥30 ℃熱積溫最為敏感。利用果實(shí)膨大高峰期-油脂轉(zhuǎn)化期的≥30 ℃熱積溫和高溫日數(shù)氣象指標(biāo)可較好預(yù)測(cè)油茶產(chǎn)量,并及時(shí)發(fā)出高溫干旱災(zāi)害預(yù)警,能夠指導(dǎo)油茶種植區(qū)及時(shí)采取預(yù)防高溫干旱災(zāi)害的管理措施,如人工降雨或灌溉補(bǔ)水,以確保油茶產(chǎn)量和經(jīng)濟(jì)效益。