海南島天然橡膠產(chǎn)量和氣候適宜度相關(guān)性研究

佟金鶴, 張衛(wèi)紅, 劉少軍, 甘業(yè)星

海南島天然橡膠產(chǎn)量和氣候適宜度相關(guān)性研究

佟金鶴1, *, 張衛(wèi)紅2, 劉少軍1, 甘業(yè)星1

1. 海南省氣象科學研究所/海南省南海氣象防災(zāi)減災(zāi)重點實驗室, ?？?570203 2. 河北省氣象科學研究所, 石家莊 050022

使用氣候適宜度指數(shù)方法對海南島橡膠樹氣候適宜性進行評價, 并選擇去趨勢互相關(guān)分析(DCCA)方法對橡膠樹單產(chǎn)和氣候適宜度進行相關(guān)性分析。結(jié)果表明: 海南島橡膠單產(chǎn)在東部和西部較低, 最低值出現(xiàn)在文昌, 為657.30 kg·hm–2; 最高值為三亞的1376.54 kg·hm–2。海南島中部單產(chǎn)的變化幅度較小, 其余地區(qū)單產(chǎn)均呈明顯增加趨勢, 線性傾向率普遍超過10 kg·hm–2·a–1。海南島中部和東部綜合適宜度水平較高, 適宜度指數(shù)超過0.70, 西部適宜度水平較低, 為0.39; 溫度和光照適宜度指數(shù)空間差異小, 整體在0.80和0.70以上; 降水適宜度指數(shù)在中部和東部較高, 西部南部偏低, 數(shù)值在0.42—0.77之間; 風速適宜度在中部和南部超過0.75, 在西部和北部適宜度指數(shù)偏低, 最低達到0.18。DCCA分析結(jié)果顯示, 氣候適宜度與橡膠單產(chǎn)存在明顯的長程相關(guān)關(guān)系, 降水和產(chǎn)量的相關(guān)性最高, 風速相關(guān)性最低; 空間上, 東部相關(guān)性較低, 西部、北部和和南部相關(guān)性稍高。

DCCA; 天然橡膠產(chǎn)量; 氣候適宜度指數(shù)

0 前言

天然橡膠是在交通、醫(yī)藥、日常生活等多領(lǐng)域起重要作用的戰(zhàn)略性資源, 隨著經(jīng)濟的增長, 我國對天然橡膠的需求量節(jié)節(jié)攀升[1]。海南作為中國天然橡膠的主產(chǎn)區(qū), 其橡膠生長氣候適宜性偏低, 過高或過低的氣溫, 較大的風速等不利氣象條件均在不同程度上影響了海南橡膠的高產(chǎn)穩(wěn)產(chǎn)[2, 3]。橡膠樹在高溫、高濕、微風的環(huán)境下才能保持旺盛的生長, 為定量衡量溫度、光照、降水、風速等氣象因子對橡膠樹生長發(fā)育、產(chǎn)量形成、品質(zhì)指標等的適宜程度, 佟金鶴等使用氣候適宜度方法評價了中國橡膠樹種植區(qū)的氣候條件[4], 亦有研究探討了橡膠樹第一蓬葉期和割膠期的氣候適宜度分布及變化特征[5, 6]。陳小敏等進一步使用相關(guān)系數(shù)法分析了割膠期各月氣候適宜度的加權(quán)系數(shù)構(gòu)成, 建立了橡膠樹割膠期的氣候適宜度評價指標[7]。橡膠產(chǎn)量除受到氣象因子的影響外, 還與橡膠樹品種、田間管理水平、橡膠市場行情等多種因素關(guān)系密切[8], 傳統(tǒng)的統(tǒng)計方法難以消除非平穩(wěn)、復(fù)雜性等特征導致的偽相關(guān)現(xiàn)象。DCCA(Detrended Cross-Correlation Analysis)是一種基于去趨勢協(xié)方差的互相關(guān)性分析方法[9], 通過系統(tǒng)性地濾去各階趨勢成分、消除原始序列中數(shù)據(jù)非平穩(wěn)性的影響, 可以有效避免由于數(shù)據(jù)非平穩(wěn)性所導致的序列之間的偽相關(guān), 從而檢測出含有噪聲且疊加有多項式趨勢信號的長程相關(guān)性[10, 11]。目前已逐步應(yīng)用于大氣環(huán)境[12, 13]、金融市場[14]、水文氣象[15]、旅游氣象[16]等多個領(lǐng)域。

目前尚未見DCCA方法在橡膠氣候適宜度與產(chǎn)量相關(guān)關(guān)系分析中的應(yīng)用。應(yīng)用去趨勢互相關(guān)分析(DCCA)模型探討氣候適宜度和橡膠產(chǎn)量的相關(guān)性, 對進一步探索氣候條件對橡膠樹的影響機制, 提升橡膠樹防災(zāi)減災(zāi)能力, 保障我國天然橡膠生產(chǎn)具有重要的參考意義。

1 數(shù)據(jù)與方法

1.1 數(shù)據(jù)

選取國家氣象信息中心《中國國家級地面氣象站基本氣象要素日值數(shù)據(jù)集(V3.0)》中海南島1989—2010年10個測站的氣溫、降水、日照時數(shù)、風速等氣象因子, 通過計算獲得22年的氣象因子適宜度及綜合適宜度序列, 并將其與所在市縣的橡膠單產(chǎn)序列進行DCCA分析。其中橡膠產(chǎn)量數(shù)據(jù)來自海南省統(tǒng)計廳。使用Matlab 函數(shù)adftest()對各市縣的橡膠單產(chǎn)(單位面積產(chǎn)量, 下同)及氣候適宜度序列進行ADF檢驗, 結(jié)果顯示研究時段內(nèi)各市縣橡膠產(chǎn)量序列及適宜度序列均存在單位根, 為非平穩(wěn)序列, 因此分析海南橡膠產(chǎn)量和氣候適宜度相關(guān)性適用DCCA方法。

1.2 方法

選取對橡膠生長具有重要影響的溫度、降水、光照、風速四個主要影響因子計算逐月適宜度函數(shù)[4–6]。隨后將逐月四個氣象因子和綜合氣候適宜度序列分別與橡膠單產(chǎn)序列進行DCCA分析。具體方法如下:

月溫度適宜度S為逐日氣溫適宜度S的月平均值,S計算方法為:

其中:為平均氣溫(℃),T=15 ℃、T=39 ℃、T= 26 ℃ 分別為橡膠生長的最低溫度、最高溫度和最適溫度[14]。若t< T, 則S=0。

月降水適宜度S:

式中,S表示月降水量適宜度, 而 S為月降水日數(shù)適宜度。R=150 mm為月適宜降水量,為月降水量(mm)。d=10 d、d=16 d為橡膠月降水適宜日數(shù)的上限和下限,為月降水日數(shù)(d)。

月光照適宜度S為逐日光照適宜度S的月平均值,S計算方法為:

式中為實際日照時數(shù)(h),=5.1為常數(shù),S為日照百分率為55%的日照時數(shù)(h), 由測站經(jīng)緯度及日序確定。

月風速適宜度S為逐日風速適宜度S的月平均值,S計算方法為:

為日平均風速(m/s),W=1.9 m/s、W=2.9 m/s分別為橡膠適宜風速的下限和上限。

月綜合適宜度S定義為:

式中,S、S、S、S分別為月溫度適宜度、月降水適宜度、月光照適宜度和月風速適宜度。

適宜度年值為各自月值的平均。

DCCA方法具體如下[10, 17]:

隨后改變時間尺度, 重復(fù)進行以上步驟得到不同時間尺度下的, 若和呈現(xiàn)線性則其線性傾向率為DCCA標度指數(shù)。=0.5時, 表明分析的兩組時間序列中不存在互相關(guān)性; 當>0.5時, 表示兩組時間序列之間存在長程正相關(guān)性,<0.5則表示兩組序列具有反持續(xù)性的長程互相關(guān)性, |-0.5|值越大, 相關(guān)性越強[11]。

2 結(jié)果與分析

2.1 海南島橡膠產(chǎn)量分析

1989—2010年間, 海南島18市縣橡膠平均單產(chǎn)和單產(chǎn)變化趨勢如下圖:

如圖1所示海南島橡膠產(chǎn)量東部和西部較低, 文昌、東方、樂東等地橡膠單產(chǎn)低于800 kg·hm–2, 其中文昌的橡膠單產(chǎn)最低, 為657.30 kg·hm–2; 澄邁、昌江-白沙-瓊中以及三亞橡膠單產(chǎn)較高, 普遍在1000 kg·hm–2以上, 單產(chǎn)最高的地區(qū)為三亞, 達到1376.54 kg·hm–2, 其次為澄邁1278.65 kg·hm–2。橡膠單產(chǎn)以增加為主, 僅樂東和保亭橡膠單產(chǎn)略有下降; 單產(chǎn)增加最明顯的地區(qū)為瓊海, 增加趨勢為46.85 kg·hm–2·a–1, 其次為臨高、文昌、三亞、定安等地, 增產(chǎn)幅度均超過30 kg·hm–2·a–1。除單產(chǎn)下降的樂東、保亭及增產(chǎn)幅度最小的萬寧、東方, 研究時段產(chǎn)量序列均通過了α=0.05的顯著性檢驗。

圖1 海南橡膠單產(chǎn)(kg·hm–2)及單產(chǎn)變化趨勢(kg·hm–2·a–1) (填色圖表示橡膠單產(chǎn), △表示橡膠單產(chǎn)變化趨勢)

Figure 1 The rubber tree yields (kg·hm–2) and linear trend (kg·hm–2·a–1) in Hainan (color map indicated rubber tree yields , △ indicated their linear trends )

2.2 海南島氣候適宜度分布及變化狀況

1989—2010年間, 海南島氣候適宜度分布和變化情況如圖2:

海南島溫度適宜度最高, 全島均在0.80以上, 基本呈現(xiàn)南高北低中部山區(qū)低于同緯度周邊地區(qū)的分布形式。降水適宜度空間差異較大, 其值在0.42—0.77之間, 適宜度最高的地區(qū)是瓊中-屯昌-瓊海一帶, 適宜度最低是西南部的三亞-東方一帶。光照適宜度整體較高, 南部高于北部, 西部高于東部, 三亞和東方是光照適宜度最高的地區(qū), 分別為0.83和0.84。風速適宜度中部較高, 北部、東部、西部適宜度均較低, 其中西部的東方風速適宜度最低, 僅為0.18。綜合適宜度以中部最高, 在0.70以上, 其次為東部, 適宜度最低的地區(qū)位于海南島的西部, 其值為0.39。由于溫度和光照適宜度空間差異較小, 海南島的綜合適宜度差異主要是由降水和風速貢獻。

圖2 海南島橡膠溫度(1)、降水(2)、光照(3)、風速(4)及綜合(5)適宜度指數(shù)分布及變化(填色圖表示氣候適宜度, △表示橡膠適宜度變化趨勢)

Figure 2 The distribution and changement of temperature(1), precipitation(2), sunshine duration(3), wind speed(4)and comprehensive(5) suitability indexes of rubber tree in Hainan Island (color map indicated rubber tree yields , △ indicated their linear trends )

在變化方面, 溫度適宜度在全島基本呈現(xiàn)一致的升高趨勢, 僅三亞出現(xiàn)下降。降水適宜度在全島呈現(xiàn)一致的下降趨勢, 下降最明顯的地區(qū)是東方, 線性傾向率為0.003 a–1。光照適宜度以下降為主, 海南島東部沿海地區(qū)光照適宜度下降最為明顯, 僅三亞和屯昌地區(qū)出現(xiàn)了升高。風速適宜度在四個因子中變化幅度最大, 在?？谌齺喅霈F(xiàn)明顯的下降, 降幅分別為–0.12和–0.08 a–1, 在五指山地區(qū)明顯上升, 幅度達到0.014 a–1。綜合適宜度在?？诤腿齺喌貐^(qū)出現(xiàn)了較為明顯的下降, 其他地區(qū)變化幅度較小。整體上, 海南島氣候適宜度變化幅度小, 大多未通過顯著性檢驗。

2.3 橡膠產(chǎn)量與氣候適宜度的去趨勢互相關(guān)分析

DCCA結(jié)果顯示, 橡膠單產(chǎn)與四個氣象因子及綜合適宜度的標度系數(shù)均在0.84以上, 超過0.5的閾值, 表明氣候適宜度與橡膠單產(chǎn)存在明顯的長程相關(guān)關(guān)系。圖3為各測站DCCA標度系數(shù)的平均值。

如圖3-1所示, 橡膠單產(chǎn)和降水適宜度的標度系數(shù)最高, 其次為溫度和光照, 風速適宜度最低,顯示橡膠單產(chǎn)對降水適宜度的變化最為敏感。相較之下, 風速適宜度對橡膠單產(chǎn)的影響稍小?？臻g分布上, DCCA標度系數(shù)在海南島東部尤其是文昌地區(qū)稍低, 在西部、北部和南部較高(圖3-2)。

3 討論

橡膠樹適宜生長在高溫、高濕、微風、光照充足的環(huán)境, 海南島光熱條件優(yōu)越、降水充分[18], 是我國最適宜種植橡膠樹的區(qū)域之一。但盛夏季節(jié)過高的氣溫易造成膠水過快凝結(jié); 過于頻繁的降水會增加橡膠病害風險[19], 可能導致大雨沖膠、產(chǎn)量降低; 微風可以促進橡膠樹蒸騰, 進而促進代謝活動, 同時有助于保持樹皮干燥, 減輕割皮病害, 但過高的風速會吹裂嫩葉, 橡膠林中年平均風速大于3 m·s-1時, 橡膠樹不能正常生長[8]。使用氣候適宜度方法評價氣象條件, 綜合考慮了作物多個關(guān)鍵氣象因子的最高、最低、最適氣象條件, 能較充分地反映橡膠樹生產(chǎn)實際[20, 21]。適宜度的研究結(jié)果顯示海南島溫度和光照適宜度較高, 降水和風速空間差異較大, 綜合適宜度以中部為最高, 西部最低。其中西部代表測站東方站位置臨海, 局地風速偏高是導致該地適宜度明顯偏低于周圍地區(qū)的主要原因。從橡膠產(chǎn)量上看, 文昌—陵水一帶橡膠單產(chǎn)較低, 可能是由于該地相對頻繁的臺風會導致橡膠斷倒, 及其他次生病害, 可對成林產(chǎn)膠量造成長期影響[22], 對斷倒比例高的膠林, 氣候適宜度的提升不足以彌補產(chǎn)量損失。三亞橡膠產(chǎn)量為全島最高, 主要是由于三亞地區(qū)光熱條件最好, 受臺風影響較少[18], 降水雖較少但可滿足橡膠樹生長需要[23]。通過與橡膠生產(chǎn)實際情況的比較分析, 認為對于海南島, 上述氣候適宜度指數(shù)方法存在對大風估計不足, 降水權(quán)重偏高的問題。去趨勢互相關(guān)分析的結(jié)果顯示, 橡膠氣候適宜度與橡膠單產(chǎn)呈明顯的長程正相關(guān)關(guān)系, 這與之前的研究結(jié)果一致[5–6]。氣候適宜度雖不能完全解釋產(chǎn)量的空間差異, 但氣候適宜度指數(shù)的提升對產(chǎn)量的提升有明顯的正向促進作用。

圖3 海南島橡膠單產(chǎn)與氣候適宜度DCCA相關(guān)關(guān)系統(tǒng)計分析(1)及空間分布(2)

Figure 3 Statistic Analysis(1) and spatial distribution(2) of the relationship between rubber tree yield and climate suitability in Hainan Island by DCCA

4 結(jié)論

1. 海南島橡膠單產(chǎn)在東部和西部較低, 最低值出現(xiàn)在文昌, 為657.30 kg·hm–2; 單產(chǎn)高值在三亞和澄邁, 分別達到1376.54和1278.65 kg·hm–2。海南島中部單產(chǎn)的變化幅度較小, 其余地區(qū)單產(chǎn)均呈明顯增加趨勢, 線性傾向率普遍超過10 kg·hm–2·a–1。

2. 海南島溫度和光照適宜度水平較高, 分別在0.80和0.70以上, 空間差異較小, 基本呈現(xiàn)南高北低的分布形式; 降水適宜度空間差異較大, 海南島中部和東部適宜度較高在0.70以上, 西部和南部適宜度較低, 普遍在0.60以下; 風速的區(qū)域性更加明顯, 海南島中部和南部適宜度高于0.75, 而西部和北部適宜度偏低, 最低值位于東方, 適宜度指數(shù)僅為0.18; 綜合適宜度以中部最高, 在0.70以上, 其次為東部, 適宜度最低的地區(qū)位于海南島的西部, 其值為0.39。

3. DCCA分析結(jié)果顯示, 氣候適宜度與橡膠單產(chǎn)存在明顯的長程相關(guān)關(guān)系。其中, 降水和產(chǎn)量的相關(guān)性最高, 風速相關(guān)性最低, 溫度和光照相關(guān)性相近?？臻g上, 海南島東部相關(guān)性較低, 西部、北部和和南部相關(guān)性相對稍高。

[1] 鞠巖峰, 張劍, 吳潤, 等. 天然橡膠種植現(xiàn)狀及市場需求預(yù)測分析[J]. 林業(yè)資源管理, 2014, (1): 152–157.

[2] 劉世紅, 田耀華. 橡膠樹抗寒性研究現(xiàn)狀與展望[J]. 廣東農(nóng)業(yè)科學, 2009, 236(11): 26–28.

[3] 王兵, 鄭璟, 杜堯東, 劉錦鑾, 王志春. 廣東橡膠風害等級標準及風險區(qū)劃研究[J]. 自然災(zāi)害學報, 2019, 28(5): 189–197.

[4] 佟金鶴, 劉少軍, 陳小敏, 等. 中國橡膠樹氣候適宜度分布特征研究[J]. 生態(tài)科學, 2021, 40 (1): 162—168.

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

[6] 陳小敏, 劉少軍, 張京紅, 等. 海南島橡膠割膠期的氣候適宜度變化特征分析[J]. 氣象與環(huán)境科學, 2019, 42(2): 35–41.

[7] 陳小敏, 李偉光, 陳匯林, 等. 海南島橡膠割膠氣候適宜度評價指標的建立及應(yīng)用——以儋州市為例[J]. 江蘇農(nóng)業(yè)科學, 2019b, 47(15): 278–281.

[8] 鄧須軍, 李玉鳳. 海南天然橡膠產(chǎn)業(yè)發(fā)展研究[M]. 北京: 中國農(nóng)業(yè)出版社, 2009.

[9] 史凱, 劉春瓊, 吳生虎. 基于DCCA方法的成都市市區(qū)與周邊城鎮(zhèn)大氣污染長程相關(guān)性分析[J]. 長江流域資源與環(huán)境, 2014, 23(11): 1633–1640.

[10] 李思川. 典型城市近地面O3濃度演化的自組織動力機制[D].吉首: 吉首大學, 2015.

[11] 劉春瓊, 劉萍, 吳生虎, 等. 基于DCCA方法分析氣候變化對四川省糧食產(chǎn)量的影響[J]. 中國農(nóng)業(yè)氣象, 2016, 37(1): 43–50.

[12] 吳波, 劉春瓊, 張嬌, 等. COVID–19期間區(qū)域大氣高污染發(fā)生的非線性動力機制[J]. 中國環(huán)境科學, 2021, 41(5): 2028–2039.

[13] 謝煥麗, 何紅弟. 香港港口PM2.5和PM10的多重分形特征[J]. 大氣與環(huán)境光學學報, 2019, 14 (3): 179–190.

[14] 曹廣喜, 謝文浩. 加密貨幣量價關(guān)系研究——基于去趨勢交叉相關(guān)分析和分位數(shù)回歸的方法[J]. 經(jīng)濟與管理研究, 2021, 42 (3): 45–63.

[15] 韋曉偉, 張洪波, 辛琛, 楊建濤, 李慈. 變化環(huán)境下流域氣象水文要素的相關(guān)性演化[J]. 南水北調(diào)與水利科技(中英文), 2020, 18 (6): 17–26.

[16] 歐陽文言. 耦合EEMD和DCCA統(tǒng)計方法分析典型旅游城市空氣質(zhì)量對旅游人數(shù)的響應(yīng)關(guān)系[D]. 吉首: 吉首大學, 2019.

[17] PODOBNIK B, STANLEY H E. Detrended Cross-Correlation Analysis: A New Method for Analyzing Two Non-stationary Time Series[J]. Physical Review Letters, 2008, 100(8): 84–102.

[18] 王春乙. 海南氣候[M]. 北京: 氣象出版社, 2014.

[19] 李晗, 冉茂, 陳海濤, 等. 植物棒孢霉葉斑病的發(fā)生及防治研究進展[J]. 植物醫(yī)生, 2020, 195(1): 15–20.

[20] 魏瑞江, 王鑫. 氣候適宜度國內(nèi)外研究進展及展望[J]. 地球科學進展, 2019, 296(6): 584–595.

[21] 邱美娟, 宋迎波, 王建林, 等. 耦合土壤墑情的氣候適宜度指數(shù)在山東省冬小麥產(chǎn)量動態(tài)預(yù)報中的應(yīng)用[J]. 中國農(nóng)業(yè)氣象, 2015, 36(2): 187–194.

[22] 胡真臻, 李增平, 單金雪, 等. 橡膠樹靈芝莖腐病病原菌鑒定及其生物學特性測定[J]. 熱帶作物學報, 2021, 42(2): 488–494.

[23] 凌禎. 西雙版納橡膠林蒸散量時空變異特征及其預(yù)報模型研究[D].昆明: 云南師范大學, 2021.

Correlation between natural rubber yield and climate suitability in Hainan Island

TONG Jinhe1, *, ZHANG Weihong2, LIU Shaojun1, GAN Yexing

1. Institute of Meteorological Science of Hainan Province / Key Laboratory of South China Sea Meteorological Disaster Prevention and Mitigation of Hainan Province, Haikou 570203, China 2. Institute of Meteorological Science in Hebei, Shijiazhuang 050022, China

The climate suitability index method was used to evaluate the climate suitability of rubber trees in Hainan Island, and the detrending cross-correlation analysis (DCCA) method was selected to analyze the correlation between rubber tree yield and climate suitability.The results showed that rubber tree yield in Hainan Island was lower in the east and west and the lowest value was 657.30 kg·hm–2which appeared at Wenchang, and the highest value was 1376.54 kg·hm–2in Sanya. Per unit area yield in the central part of Hainan Island had a small range of change, and in other areas showed a significant increase, with linear tendency rates generally exceeding 10 kg·hm–2·a–1. The central and eastern parts of Hainan Island had a relatively high level of comprehensive suitability, with suitability indexes exceeding 0.70, and a lower suitability level of 0.39 in the west. The temperature and sunshine duration suitability indexes were generally above 0.80 and 0.70, which had small differences in spatial. Precipitation suitability was higher in the central and eastern parts, and was lower in the western and southern parts, overall range from 0.42 to 0.77. The wind speed suitability index was more than 0.75 in the central and southern parts, and was lower in the western and northern parts, with the lowest being 0.18. The results of DCCA analysis showed that there was an obvious long-range correlation between climate suitability and rubber tree yield. The correlation level between precipitation and yield was the highest, and the wind speed was the lowest. Spatially, the correlation in east was weak, and the correlation in the west, north and south was slightly stronger.

DCCA; natural rubber yield; climate suitability index

10.14108/j.cnki.1008-8873.2024.01.018

P464

A

1008-8873(2024)01-154-06

2021-08-29;

2022-01-07基金項目:國家自然科學基金(41765007); 海南省基礎(chǔ)與應(yīng)用基礎(chǔ)研究計劃(自然科學領(lǐng)域)高層次人才項目(2019RC359)

佟金鶴(1991—) , 女, 碩士, 工程師, 主要從事橡膠氣象服務(wù), E-mail: tongjinhesuns@163.com

通信作者:佟金鶴

佟金鶴, 張衛(wèi)紅, 劉少軍, 等. 海南島天然橡膠產(chǎn)量和氣候適宜度相關(guān)性研究[J]. 生態(tài)科學, 2024, 43(1): 154–159.

TONG Jinhe, ZHANG Weihong, LIU Shaojun, et al. Correlation between natural rubber yield and climate suitability in Hainan Island[J]. Ecological Science, 2024, 43(1): 154–159.