陳小敏 李偉光 陳匯林 劉少軍 張京紅 吳青蘭

摘要：建立海南儋州地區(qū)橡膠割膠期氣候適宜性評價指標和模型。對比絕對值法、歸一化法、相關(guān)系數(shù)法等3種方法確定權(quán)重系數(shù)建立氣候適宜指數(shù)的優(yōu)劣;建立產(chǎn)量豐歉指數(shù)和氣候適宜度指數(shù)的關(guān)系模型，確定氣候適宜性評價指標。研究結(jié)果表明，相關(guān)系數(shù)法建立氣候適宜度指數(shù)效果最優(yōu);產(chǎn)量豐歉指數(shù)與氣候適宜度指數(shù)整體呈顯著或極顯著相關(guān)關(guān)系;各個時段指標回代檢驗準確率在70%以上。研究結(jié)果可用于評價氣候條件對橡膠割膠和產(chǎn)量形成的適宜程度。

關(guān)鍵詞：天然橡膠;割膠;氣候適宜性;評價指標;儋州

中圖分類號： S162.5

文獻標志碼： A

文章編號：1002-1302（2019）15-0278-04

天然橡膠是國際戰(zhàn)略物質(zhì)，直接關(guān)系到國家經(jīng)濟發(fā)展、政治穩(wěn)定和國防安全。海南島天然橡膠種植規(guī)模占全國的1/2左右，至2016年，海南島橡膠種植面積達54.1萬hm2，干膠年產(chǎn)量35萬t[1]。而儋州市位于海南島西北部，是橡膠海南引種的第1站，也是全省種植面積和產(chǎn)量排名第1的市（縣），種植面積和產(chǎn)量占全島的1/5左右[1]，已成為該地區(qū)農(nóng)戶收入的主要來源。橡膠產(chǎn)量和質(zhì)量容易受氣候生態(tài)環(huán)境的影響和制約，如2005年“達維”臺風等造成儋州市橡膠林不同程度的折枝或倒伏[2-5];2008年年初海南島持續(xù)24 d的低溫陰雨過程，導致西北部內(nèi)陸橡膠樹出現(xiàn)嚴重寒害災害[6-9];2010年海南出現(xiàn)嚴重干旱，導致大量橡膠樹葉片枯黃脫落，嚴重影響割膠[10]。可見橡膠樹生長發(fā)育、割膠作業(yè)和產(chǎn)量與氣候條件適宜性密切相關(guān)[11-14]。因此建立橡膠割膠氣候適宜性評價指標，對科學指導和調(diào)整橡膠種植產(chǎn)區(qū)，合理利用氣候資源，最大限度地避免、抵御不利氣候條件的影響，提高橡膠產(chǎn)量具有重要意義。

近年來，隨著現(xiàn)代農(nóng)業(yè)生產(chǎn)水平的發(fā)展，農(nóng)業(yè)生產(chǎn)對農(nóng)業(yè)氣象業(yè)務提出了定量化和精細化的要求，如利用作物對日照、溫度、水分的適宜情況來研究作物適宜種植區(qū)，評價氣象要素對作物生長發(fā)育和產(chǎn)量形成的優(yōu)劣或適宜程度，獲得了很多研究成果[15-17]。也有專家利用氣候適宜度模型，結(jié)合作物氣象產(chǎn)量進行統(tǒng)計分析，確定評價氣候適宜度診斷指標的閾值，構(gòu)建適宜度診斷指標[18-20]。關(guān)于天然橡膠的氣候適宜度研究鮮有報道，劉少軍等利用溫度、降水量和降水日數(shù)、日照、風速等氣象要素建立模型，定量評價海南橡膠在葉片生長期的氣候適宜性[21]，而對割膠期的氣候適宜性分析未見報道。因此，本研究通過建模的方法，構(gòu)建適宜性評價指標，對儋州地區(qū)橡膠割膠氣候適宜性進行定量分析，以期為橡膠產(chǎn)量豐歉提供綜合定量分析方法，為政府、橡膠期貨和農(nóng)戶的橡膠生產(chǎn)、產(chǎn)量估算提供參考。

1 材料與方法

1.1 資料來源

儋州市位于海南省西北部，地處東亞大陸季風氣候的南緣，屬熱帶季風氣候，夏無酷暑，冬無嚴寒，陽光充足，雨量充沛。常年平均氣溫為23.5 ℃，最熱月（7月）的平均溫度為 28.0 ℃，最冷月（1月）的平均氣溫為17.8 ℃;年平均降水量為1 815 mm，年均日照時數(shù)為1 967 h。

本研究采用海南省氣象信息中心提供的儋州市氣象站（109.5833°E、19.5167°N）氣象資料，選取1961—2017年逐日氣象要素，包括平均氣溫、最低氣溫、降水量、降水日數(shù)、日照時數(shù)和風速等進行分析。1989—2017年橡膠單產(chǎn)數(shù)據(jù)來源于海南省統(tǒng)計局。

1.2 數(shù)據(jù)處理方法

橡膠單產(chǎn)豐歉指數(shù)的計算方法[22]，

式中：Kyi是第i年產(chǎn)量豐歉指數(shù);yi是第i年的實際單產(chǎn)值;yi是近5年單產(chǎn)的滑動平均值。

2 橡膠割膠氣候適宜度模型的構(gòu)建

2.1 單要素適宜度模型

參考侯英雨等的方法[16-21]，分別建立橡膠割膠溫度、日照時數(shù)、降水量、降水日數(shù)、風速的適宜度模型。

溫度適宜度模型為

式中：T表示日割膠溫度（由于割膠作業(yè)通常安排在02：00—09：00，因此T值處理為日最低氣溫和日平均氣溫值的平均值）;T1、T2、T0分別為研究時間段內(nèi)橡膠生長的最低溫度、最高溫度、最適宜溫度;ST表示溫度為T時的溫度適宜度;B表示最高溫度和最適宜溫度的差值與最適宜溫度和最低溫度差值之比。

2.2 多要素氣候適宜度模型

氣候適宜度模型綜合考慮了溫度、日照時數(shù)、降水量、降水日數(shù)和風速等多個要素對割膠的影響，在文獻[16-17]的基礎(chǔ)上，采用幾何平均和綜合乘積的方法，建立橡膠割膠期綜合氣候適宜度模型：

2.3 割膠適宜度模型中指標的確定

參照華南熱帶作物研究院等的相關(guān)文獻[11-14]，根據(jù)橡膠樹生長對溫度、降水、光照、風速條件的要求，確立橡膠割膠適宜度模型中各個割膠指標（表1）。

2.4 割膠期氣候適宜度模型的構(gòu)建

天然橡膠的產(chǎn)量高低與其生長過程中的光、溫、水、風條件配置有緊密關(guān)系，而且在不同生長時段對光、溫、水、風的需求和敏感度存在很大差異。橡膠割膠期的氣候適宜度指數(shù)，主要由割膠期間各月氣候適宜度的加權(quán)集成構(gòu)成，加權(quán)過程中權(quán)重系數(shù)的確定對指數(shù)的構(gòu)建及最終產(chǎn)量預報的準確性具有重要意義。其中，權(quán)重系數(shù)采用絕對值法、歸一化法和相關(guān)系數(shù)法等3種方法進行計算[23-24]。

參照氣象行業(yè)標準《主要糧食作物產(chǎn)量年景等級》[22]可將年份分為豐年、偏豐年、持平略增、持平略減、偏歉年、歉年等6種類型，對應的產(chǎn)量豐歉指數(shù)分別為Kyi>11%、4%

[10]李海亮，戴聲佩，陳幫乾，等. 基于HJ-1A/1B數(shù)據(jù)的天然橡膠干旱監(jiān)測[J]. 農(nóng)業(yè)工程學報，2016，32（23）：176-182.

[11]郭玉清，張 汝. 氣象條件與橡膠樹產(chǎn)膠量的關(guān)系[J]. 云南熱作科技，1980（1）：8-11.

[12]楊 銓. 幾種氣象因子與產(chǎn)膠量的關(guān)系[J]. 中國農(nóng)業(yè)氣象，1987（1）：42-44.

[13]李國堯，王權(quán)寶，李玉英，等. 橡膠樹產(chǎn)膠量影響因素[J]. 生態(tài)學雜志，2014，33（2）：510-517.

[14]華南熱帶作物學院. 橡膠栽培學[M]. 2版.北京：農(nóng)業(yè)出版社，1989.

[15]譚方穎，宋迎波，毛留喜，等. 東北地區(qū)玉米氣候適宜評價指標的確定與驗證[J]. 干旱地區(qū)農(nóng)業(yè)研究，2016，34（5）：234-239.

[16]侯英雨，張艷紅，王良宇，等. 東北地區(qū)春玉米氣候適宜度模型[J]. 應用生態(tài)學報，2013，24（11）：3207-3212.

[17]金志鳳，葉建剛，楊再強，等. 浙江省茶葉生長的氣候適宜性[J]. 應用生態(tài)學報，2014，25（4）：967-973.

[18]宋迎波，王建林，李昊宇，等. 冬小麥氣候適宜診斷指標確定方法探討[J]. 氣象，2013，39（6）：768-773.

[19]張建軍，馬曉群，許 瑩. 安徽省一季稻生長氣候適宜性評價指標的建立與試用[J]. 氣象，2013，39（1）：88-93.

[20]劉琰琰，李海燕，陳 超，等. 攀西地區(qū)烤煙氣候適宜性評價指標建立及應用[J]. 四川農(nóng)業(yè)大學學報，2015，33（3）：299-305.

[21]劉少軍，房世波. 海南島天然橡膠氣候適宜性及變化趨勢分析——以第一蓬葉生長期為例[J]. 農(nóng)業(yè)現(xiàn)代化研究，2015，36（6）：1062-1066.

[22]中國氣象局. 主要糧食作物產(chǎn)量年景等級： QX/T 335—2016[S]. 北京：中國氣象出版社，2017.

[23]黃 璜. 中國紅黃壤地區(qū)作物生產(chǎn)的氣候生態(tài)適應性研究[J]. 自然資源學報，1996，11（4）：340-346.

[24]邱美娟，宋迎波，王建林，等. 農(nóng)業(yè)氣象產(chǎn)量預報方法中權(quán)重系數(shù)的確定[J]. 氣象與環(huán)境學報，2016，32（2）：106-111.于 芳，李海明. 基于地統(tǒng)計學和GIS的湖北省土壤有機質(zhì)空間變異性研究[J]. 江蘇農(nóng)業(yè)科學，2019，47（15）：282-286.