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

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

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

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

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

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

0 前言

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

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

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

1.1 數(shù)據(jù)

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

1.2 方法

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

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

其中:為平均氣溫(℃),T=15 ℃、T=39 ℃、T= 26 ℃ 分別為橡膠生長(zhǎng)的最低溫度、最高溫度和最適溫度[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計(jì)算方法為:

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

月風(fēng)速適宜度S為逐日風(fēng)速適宜度S的月平均值,S計(jì)算方法為:

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

月綜合適宜度S定義為:

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

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

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

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

2 結(jié)果與分析

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

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

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

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

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。風(fēng)速適宜度中部較高, 北部、東部、西部適宜度均較低, 其中西部的東方風(fēng)速適宜度最低, 僅為0.18。綜合適宜度以中部最高, 在0.70以上, 其次為東部, 適宜度最低的地區(qū)位于海南島的西部, 其值為0.39。由于溫度和光照適宜度空間差異較小, 海南島的綜合適宜度差異主要是由降水和風(fēng)速貢獻(xiàn)。

圖2 海南島橡膠溫度(1)、降水(2)、光照(3)、風(fēng)速(4)及綜合(5)適宜度指數(shù)分布及變化(填色圖表示氣候適宜度, △表示橡膠適宜度變化趨勢(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)一致的升高趨勢(shì), 僅三亞出現(xiàn)下降。降水適宜度在全島呈現(xiàn)一致的下降趨勢(shì), 下降最明顯的地區(qū)是東方, 線性?xún)A向率為0.003 a–1。光照適宜度以下降為主, 海南島東部沿海地區(qū)光照適宜度下降最為明顯, 僅三亞和屯昌地區(qū)出現(xiàn)了升高。風(fēng)速適宜度在四個(gè)因子中變化幅度最大, 在?？谌齺喅霈F(xiàn)明顯的下降, 降幅分別為–0.12和–0.08 a–1, 在五指山地區(qū)明顯上升, 幅度達(dá)到0.014 a–1。綜合適宜度在?？诤腿齺喌貐^(qū)出現(xiàn)了較為明顯的下降, 其他地區(qū)變化幅度較小。整體上, 海南島氣候適宜度變化幅度小, 大多未通過(guò)顯著性檢驗(yàn)。

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

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

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

3 討論

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

圖3 海南島橡膠單產(chǎn)與氣候適宜度DCCA相關(guān)關(guān)系統(tǒng)計(jì)分析(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)高值在三亞和澄邁, 分別達(dá)到1376.54和1278.65 kg·hm–2。海南島中部單產(chǎn)的變化幅度較小, 其余地區(qū)單產(chǎn)均呈明顯增加趨勢(shì), 線性?xún)A向率普遍超過(guò)10 kg·hm–2·a–1。

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

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

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

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

[3] 王兵, 鄭璟, 杜堯東, 劉錦鑾, 王志春. 廣東橡膠風(fēng)害等級(jí)標(biāo)準(zhǔn)及風(fēng)險(xiǎn)區(qū)劃研究[J]. 自然災(zāi)害學(xué)報(bào), 2019, 28(5): 189–197.

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

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

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

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

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

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

[10] 李思川. 典型城市近地面O3濃度演化的自組織動(dòng)力機(jī)制[D].吉首: 吉首大學(xué), 2015.

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

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

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

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

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

[16] 歐陽(yáng)文言. 耦合EEMD和DCCA統(tǒng)計(jì)方法分析典型旅游城市空氣質(zhì)量對(duì)旅游人數(shù)的響應(yīng)關(guān)系[D]. 吉首: 吉首大學(xué), 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ìn)展[J]. 植物醫(yī)生, 2020, 195(1): 15–20.

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

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

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

[23] 凌禎. 西雙版納橡膠林蒸散量時(shí)空變異特征及其預(yù)報(bào)模型研究[D].昆明: 云南師范大學(xué), 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基金項(xiàng)目:國(guó)家自然科學(xué)基金(41765007); 海南省基礎(chǔ)與應(yīng)用基礎(chǔ)研究計(jì)劃(自然科學(xué)領(lǐng)域)高層次人才項(xiàng)目(2019RC359)

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

通信作者:佟金鶴

佟金鶴, 張衛(wèi)紅, 劉少軍, 等. 海南島天然橡膠產(chǎn)量和氣候適宜度相關(guān)性研究[J]. 生態(tài)科學(xué), 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.