陳彥希，婁運生，任麗軒，蘇 磊，湯麗玲，楊建洲

未來氣候變化情景下海南島綠橙種植氣候適宜性分析*

陳彥希1,2,3，婁運生1,3**，任麗軒4,5,6**，蘇 磊7**，湯麗玲7，楊建洲7

(1. 南京信息工程大學氣象災害預報預警與評估協(xié)同創(chuàng)新中心，南京 210044；2. 海南省氣象探測中心，?？?570203；3. 南京信息工程大學江蘇省農業(yè)氣象重點實驗室，南京 210044；4.南京農業(yè)大學資源與環(huán)境科學學院，南京 210095；5. 農業(yè)農村部長江中下游植物營養(yǎng)與肥料重點實驗室，南京 210095；6. 江蘇省有機固體廢棄物資源化協(xié)同創(chuàng)新中心，南京 210095；7. 中國地質科學院地球物理地球化學勘查研究所，廊坊 065000)

基于海南島主島19個氣象站點連續(xù)40a（1980-2019年）逐日氣象數據及地形等要素，建立空間分析模型，綜合確定海南島綠橙種植氣候適宜區(qū)指標。在氣候適宜度模型基礎上構建海南島綠橙種植溫度、日照、降水以及綜合適宜度模型，采用地理信息系統(tǒng)（GIS）技術綜合分析，利用自然點斷法對綠橙種植氣候適宜性進行精細化區(qū)劃，并利用RCP4.5未來氣候情景模式模擬數據探討未來30a（2020-2049年）氣候變化情景下海南島綠橙種植適宜區(qū)變化趨勢。結果表明：現階段海南島綠橙種植最適宜區(qū)主要分布于中部地區(qū)，面積0.87×104km2，氣候適宜度為0.9～1.0；適宜區(qū)主要分布于東部局部地區(qū)及中西部大部地區(qū)，面積1.83×104km2，氣候適宜度為0.7～0.9；次適宜區(qū)主要分布于西部沿海及中西部海拔較高地區(qū)，面積0.51×104km2，氣候適宜度為0.4～0.7；不適宜區(qū)主要分布于中部山區(qū)，面積0.17×104km2，氣候適宜度為0～0.4。未來氣候變化情景下海南島氣溫、降水分布呈較大變化，溫度適宜區(qū)由四周向中部逐漸縮小，降水適宜區(qū)由東部逐漸遷移至中部地區(qū)。未來綠橙種植氣候最適宜區(qū)主要分布于瓊中、屯昌和保亭縣大部分地區(qū)、萬寧市西部及白沙縣中東部等地。

綠橙；適宜性；氣候區(qū)劃；GIS；未來氣候變化情景

綠橙是海南島知名特色水果之一，前身為紅江橙，20世紀80年代由廣東廉江引入海南后，表現出比原產地更強的適應性，具有結實早、豐產快等特點[1?2]。海南島綠橙已獲農業(yè)農村部無公害農產品及綠色食品認證，是海南優(yōu)質名牌農產品[3]，種植面積約占該省柑橘橙柚類種植面積的36.4％。

IPCC第五次報告指出，全球氣溫近100a上升了0.85～0.89℃[4]。氣候是影響作物生長、產量及品質的重要因素，因此，氣候變化背景下作物氣候適宜性及其變化趨勢日益受到關注[5]。目前，已有學者采用模糊數學法、層次分析法、加權指數法、聚類分析法、空間分析法等方法對作物氣候適宜性區(qū)劃開展了研究[6?7]。沈兆敏等[8?9]利用氣溫、≥10℃積溫、極端最低氣溫等氣溫指標繪制了全國柑橘氣候區(qū)劃圖；池再香等[10]運用模糊聚類法研究了貴州東南部椪柑氣候區(qū)劃并將椪柑分為四個種植適宜區(qū)；林紅等[11]研究了四川省大英縣甜橙種植氣候區(qū)劃；也有學者在柑橘生產氣候適宜性研究中考慮了氣象災害因子[12]，如王迅等[13]報道了四川省晚熟柑橘成熟所需光熱條件，在原區(qū)劃指標基礎上采用特爾斐法疊加冬春氣象災害因子，解決了晚熟柑橘冬春季產量低的問題。前人對柑橘橙柚類果樹的生產種植、氣候適宜性區(qū)劃等方面已開展相關研究[14?15]，但迄今有關海南島綠橙種植氣候適宜性研究尚不多見[16]，現有相關研究存在所用氣象數據時間序列較短較早，所用站點數據分布相對于綠橙種植空間差異性要求的精度偏低的問題，難以滿足綠橙種植精細化區(qū)劃要求。此外，隨著氣候變化及綠橙栽培技術進步，傳統(tǒng)的氣候適宜性區(qū)劃已不能滿足綠橙優(yōu)質穩(wěn)產的要求。海南島地處中國最南端，具有獨特的氣候和地理條件，為合理布局綠橙種植，有必要采用適當的適宜度模型的評價方法開展綠橙氣候適宜性精細化區(qū)劃研究。針對以上問題，本研究基于海南島19個主要氣象臺站近40a（1980-2019年）逐日氣象觀測資料，結合多點地形等數據要素，運用數理統(tǒng)計方法構建空間分析模型，根據海南島綠橙種植氣候條件確定區(qū)劃指標，在氣候適宜度函數基礎上構建海南綠橙溫度、日照、降水以及綜合適宜度模型。在分析現階段綠橙種植氣候適宜性基礎上，利用未來氣候變化情景（RCP4.5）模式模擬氣象數據，探討未來30a（2020-2049年）海南綠橙種植適宜區(qū)對氣候變化的響應?；贕IS技術利用自然點斷法，開展海南島綠橙種植氣候適宜性精細區(qū)劃，評估復雜地形下綠橙種植氣候適宜性，旨在為海南島綠橙種植優(yōu)化空間布局及農業(yè)氣候資源開發(fā)利用提供生產參考依據。

1 資料與方法

1.1 研究區(qū)概況

海南島位于中國最南端，下屬18個縣市，陸地面積約3.39×104km2[17]。該地氣候為熱帶季風氣候，年均氣溫23.2～27.1℃，總降水量1009～2367.7mm，日照時數1776～2783h，島內氣象觀測站點分布見圖1。

圖1 海南島主島19個氣象觀測站點分布

1.2 數據來源

（1）1980-2019年海南島主島19個氣象站點逐日氣象資料，包括日平均氣溫、日總降水量、日照等氣象要素，來自國家氣象數據網（http://data. cma.cn/）。

（2）海南島主島DEM數字高程模型數據。利用DEM數據獲取坡度、坡向、海拔高度等數據（30m×30m），來自國家基礎地理信息中心（http:// www. ngcc.cn/ngcc/）。

（3）海南島1:25萬行政區(qū)劃矢量圖層，包括海南島邊界、18市縣邊界等，來自資源環(huán)境科學與數據中心（https://www.resdc.cn/）。

（4）海南島未來氣候變化情景數據，即RCP4.5排放情景下CMIP5（Coupled Model Intercomparison Project Phase 5）模式模擬數據，包括2020-2049年溫度數據和降水數據來自http://data.ceda.ac.uk/。

（5）海南島柑橘橙柚、綠橙農產品種植與生產情況相關數據來自海南省統(tǒng)計年鑒（2020）。

1.3 研究方法

1.3.1 氣候資源分布

海南島綠橙種植區(qū)主要分布于中部及東北部地區(qū)（表1），種植區(qū)年平均氣溫24.2℃，年平均極端最低溫度7.8℃，≥10℃平均年積溫8851.1℃·d，1月平均氣溫18.1℃，7月平均氣溫28.3℃；年均降水量2023.5mm，相對濕度85.3％，干燥度0.75，降水充足，但季節(jié)降水分布不均，多年平均降水量冬季90mm，春季321mm，夏季721mm，秋季676.8mm；年均總日照時數1964.8h，日照充足。

表1 綠橙種植主要市縣1980?2019年氣象要素特征

Note: T is mean temperature, D is annual mean cumulative temperature, R is annual mean precipitation, ST is sunshine hours.

1.3.2 氣象數據空間分析模型建立

基于海南島主島19個氣象站點近40a（1980-2019年）氣溫、降水及日照等逐日氣象數據，運用數理統(tǒng)計方法構建空間分析模型。因氣候與經緯度、海拔高度、坡向和坡度等地理因子相關，故區(qū)劃指標與地理因子多元回歸關系模型可表示為[18?19]

式中，Y為氣候區(qū)劃指標因子，?、λ、h、β和θ分別為緯度（°）、經度（°）、海拔高度（m）、坡向（°）、坡度（°）等地理因子，e為余差項，稱為地理殘差。

1.3.3 氣候適宜度模型建立

1.3.3.1 關鍵氣候指標篩選及氣候適宜度區(qū)劃指標閾值確定

參考前人研究，分析海南島綠橙生長發(fā)育所需氣候條件，篩選出影響綠橙生長發(fā)育的關鍵氣候指標[2,16,20]，包括年平均氣溫、年日照時數和年降水量。為定量分析研究區(qū)氣候條件對綠橙種植適宜度的影響，結合相關文獻和研究，分別建立綠橙溫度、日照時數、降水的適宜度模型。結合相關文獻[2,14,20-28]，基于海南島19個氣象站點氣象數據及地理數據，以及綠橙種植氣候條件要求，確定綠橙適宜度模型各指標。

1.3.3.2 溫度適宜度模型

氣溫是影響綠橙生長發(fā)育和產量的決定因素[16]。張業(yè)忠[2]通過農業(yè)氣候分析提出，最適宜綠橙生長的年平均氣溫為22～24℃，溫度過高或過低均會使果實的品質受到影響[12]?？紤]綠橙不同生育期溫度差異，選用1月平均氣溫（最冷月）及7月平均氣溫（最熱月）數據，作為溫度適宜度區(qū)劃指標。作物溫度適宜度模型多采用與作物三基點溫度相匹配的非線性beta函數[21]，溫度適宜度模型為[22?23]

式中，ST為作物溫度適宜度（℃），t為實際溫度（℃），t0、tl、th分別為作物生長最適、下限和上限溫度（℃），具體取值見表2。

1.3.3.3 日照適宜度模型

日照影響果實質量、色澤和品質，日照充足可增加果實質量，使果實飽滿、著色良好[24]。日照時數過低或過高均導致果實發(fā)育不良，特別是炙熱的夏季，光照過強會導致綠橙日灼病、果實呈海綿狀等疾病[25]。當日照時數達到可照時數的55％時，作物對光照反應達到最適宜狀態(tài)[26]。日照適宜度模型一般采用自然數e為底數的指數函數分段函數。日照適宜度函數模型為[21?22]

式中，SS為作物日照適宜度，S為實際日照時數（h），S0為特定地區(qū)特定時間可照時數（h），取值在4.9～5.5，b為常數，取值為5.1[27]。

1.3.3.4 降水適宜度模型

水分影響綠橙果實大小、質地和風味等品質指標[28]，綠橙種植最佳年降水量為2000～2200mm[20]。作物降水適宜度模型為[27?29]

式中，SP為作物降水適宜度，Sr為作物降水量適宜度，Sd為作物降水日數適宜度。

表2 溫度適宜度指標

式中，R為實際降水量（mm），R0為作物生長最適降水量（mm），取值2100mm。

式中，d為實際年降水日數（d），dl、dh分別為作物生長適宜降水日數的下限和上限（d），取值分別為140d和200d。

1.3.3.5 綜合適宜度模型

根據單因子氣候適宜度區(qū)劃結果，利用等權重求和法、幾何平均法和不等權求和法，分別建立綜合氣候適宜度模型[30?31]，取其均值。

等權重求和法計算式為

幾何平均法計算式為

不等權求和法計算式為

式中，S為作物綜合適宜度，WT、WS和WP分別為作物溫度、日照時數以及降水的適宜度權重，均采用層次分析法（AHP）所確定[23]。

2 結果與分析

2.1 未來海南島綠橙種植關鍵氣候因子時空變化

2.1.1 氣候因子空間分析模型

利用1980-2019年海南島19個氣象站點年平均氣溫、年總日照時數、年總降水量3個氣象要素，分別與經度、緯度、海拔高度、坡向和坡度5個地理信息因子進行回歸分析，構建空間分析模型[32?34]，結果見表3。

由表3可見，各模型的相關系數均通過了0.01水平的顯著性檢驗。其中，年平均氣溫與緯度、經度、海拔高度、坡度和坡向均呈負相關關系；年總日照時數與緯度、經度、坡度和坡向呈負相關關系，與海拔高度呈正相關關系；年總降水量與緯度和海拔高度呈負相關關系，與經度、坡度和坡向呈負相關關系；年≥10℃積溫與緯度、經度、海拔高度、坡度和坡向均呈負相關關系。

2.1.2 關鍵氣象因子時空變化

2.1.2.1 溫度

由圖2a可見，海南島年平均氣溫總體呈四周高、中間低的分布格局，其中西南地區(qū)溫度最高。從40a平均狀況看，中部地區(qū)年均溫度低于21.0℃，西部及北部地區(qū)年均溫度高于25.0℃。從RCP4.5情景下2020?2049年30a平均狀況看（圖2b），中部地區(qū)年均低于21.0℃的區(qū)域較過去40a范圍更廣，面積增加0.33×103km2，約為原面積的23％，且四周沿海大部分地區(qū)年均溫度均高于25℃，較圖2a范圍有大幅增加，面積增加1.14×103km2，約為原面積的35％。總體上，未來30a海南島平均氣溫較過去40a均有不同程度增加，但海拔較高地區(qū)則相反，低溫地區(qū)也有擴展趨勢。因綠橙生長適宜年平均氣溫為22～24℃，因此，綠橙種植溫度適宜區(qū)范圍呈縮小趨勢。

2.1.2.2 降水量

海南島降水強度呈現東高西低、四周向中間漸增趨勢。從1980-2019年40a平均狀況看（圖3a），中部部分地區(qū)年均總降水量＞2400mm，西部沿海地區(qū)年均總降水量＜1800mm。從RCP4.5情景下2020-2049年30a平均狀況看（圖3b），年均總降水量＞2400mm區(qū)域與過去40a相比，其范圍由中部部分地區(qū)擴至中部大部分地區(qū)，面積增加0.43×103km2，約為原面積的14％；東部大部分區(qū)域年降水量均＞2200mm。總體上，未來30a海南島平均年降水量較過去40a均有不同程度增加，且年降水量＜1800mm區(qū)域范圍縮小。由于綠橙生長適宜年平均總降水量在2000～2200mm，因此，綠橙種植降水適宜區(qū)范圍由東部逐漸遷移至中部地區(qū)。

表3 海南島氣象要素空間分析模型

注：**表示相關系數通過0.01 水平顯著性檢驗。Φ,λ，h，β，θ分別為緯度坐標值、經度坐標值、海拔高度(m)、坡度(°)、坡向（°）。

Note:**is P＜0.01.φ is latitude coordinate value of kilometer network，λis longitude coordinate value of kilometer network，h is Altitude (m)，β is slope (°)，θ is slope direction (°).

圖2 海南島過去40a(a)和未來RCP4.5情景下(b)年平均氣溫空間分布

圖3 海南島過去40a(a)和未來RCP4.5情景下(b)年平均降水量空間分布

2.2 未來海南島綠橙種植氣候適宜性變化

利用氣候適宜度模型按照綠橙生長適宜度評價指標進行計算，在ArcGIS空間分析模塊中生成單個氣候適宜度柵格圖，即形成氣象要素單因子氣候適宜度空間分布圖（圖4?圖6）。

（1）溫度適宜度：基于海南島近40a氣象數據計算當前溫度適宜度，繪制空間分布圖。由圖4a可見，40a內全島溫度適宜度總體呈四周高中間低的分布趨勢，溫度適宜度范圍在?1～1，海南島綠橙生長溫度最適宜區(qū)，分布于中部，以瓊中、屯昌及周邊部分地區(qū)為主，該區(qū)年平均氣溫在22～24℃；溫度次適宜區(qū)，一部分位于北部及東北部，為?？?、澄邁、儋州、定安等地，該區(qū)年平均氣溫在24～25℃，另一部分位于白沙、昌江等局部區(qū)域，該地年平均氣溫在21～22℃；溫度不適宜區(qū)，一部分位于東部和北部，以東方、三亞、樂東等地，該區(qū)年平均氣溫＞25℃，另一部分位于昌江、五指山、保亭等地，該地年平均氣溫＜21℃，不適宜綠橙生長。

圖4 海南島過去40a(a)和未來RCP4.5情景下(b)綠橙種植溫度適宜度空間分布

基于海南島未來30a氣象數據，計算得到未來溫度適宜度，繪制空間分布圖。由圖4b可見，全島未來溫度適宜度總體呈中間高四周低的分布趨勢，未來溫度適宜度范圍在?1～1，綠橙生長溫度最適宜區(qū)范圍相較于歷史時期（1980?2019年）有縮小趨勢，不適宜區(qū)范圍也呈擴大趨勢，主要分布于海南島西部及北部沿海地區(qū)。

（2）日照適宜度：基于海南島過去40a氣象數據，計算當前日照適宜度，繪制空間分布圖。由圖5可見，40a內全島日照適宜度總體呈現由西向東遞減趨勢，日照適宜度范圍在0.63～1.00，綠橙生長日照時數最適宜區(qū)，分布于中部及中東部，以澄邁、定安、屯昌、瓊中等市縣及周邊局部為主，該區(qū)年日照時數1800～2000h；日照時數次適宜區(qū)，一部分位于東北部，另一部分由中部向南、北發(fā)散，包括文昌、臨高、儋州等市縣，年日照時數1600～1800h或2000～2200h；日照時數不適宜區(qū)，小部分位于中部河流沿岸，大部分位于西部，包括東方、樂東等市縣及周邊區(qū)域，該區(qū)年均日照時數＞2200h，不適宜綠橙生長。因未來氣候情景數據缺少日照數據，故僅對當前綠橙種植氣候適宜度進行分析。

圖5 海南島過去40a(1980?2019年)綠橙種植日照適宜度空間分布

（3）降水適宜度：由海南島過去40a氣象數據計算當前降水適宜度，繪制空間分布圖。由圖6a可見，40a內全島降水適宜度呈現東高西低趨勢，降水適宜度范圍在0.45～1.00，綠橙生長發(fā)育降水最適宜區(qū)，位于東部及中部局部地區(qū)，以?？?、文昌、定安等縣市為主，年降水量2000～2200mm；次適宜區(qū)，一部分位于北部及中部局部地區(qū)，包括澄邁、瓊中等地，年降水量1800～2000mm，另一部分位于中部零散區(qū)域，年降水量2200～2400mm；不適宜區(qū)，一部分位于西部及西北部，以儋州、昌江、東方、樂東等市縣及周邊為主，年降水量＜1800mm，另一部分位于中部局部地區(qū)，年降水量＞2400mm。由海南島未來30a氣象數據計算未來降水適宜度，繪制空間分布圖，由圖6b可見，未來降水適宜度總體也呈東高西低的分布趨勢，未來降水適宜度范圍在0.66～1.00，綠橙生長降水最適宜區(qū)范圍，相較于圖6a有縮小趨勢，由東部逐漸遷移至中部地區(qū)，不適宜區(qū)范圍呈擴大趨勢，主要分布于海南島西部及北部沿海地區(qū)。

2.3 海南島綠橙種植氣候適宜性分區(qū)評估和未來可能變化分析

考慮氣候條件、地形因素及土地利用類型等多因子進行區(qū)劃分級[22]，即得海南島當前和未來RCP4.5情景下綠橙種植最適宜區(qū)、適宜區(qū)、次適宜區(qū)及不適宜區(qū)（圖7），因未來氣候情景數據缺少日照時數，故未來氣候變化情景下綠橙氣候適宜性區(qū)劃僅用溫度、降水適宜度進行綜合分析。

（1）最適宜區(qū)

由圖7a可知，過去40a最適宜區(qū)主要分布于海南島中部區(qū)域，包括?？谑心喜?、?？谑协偵絽^(qū)、澄邁縣南部、定安縣南部區(qū)域、文昌東北部小部分區(qū)域、瓊海市西北部分區(qū)域、屯昌縣大部分地區(qū)、儋州市南部小部分區(qū)域、萬寧市西部、瓊中縣大部分地區(qū)、白沙縣中部及南部地區(qū)、五指山市小部分區(qū)域等，面積約0.87×104km2。該區(qū)年均氣溫22.8～24.7℃，≥10℃年均積溫＞8100℃·d，年均日照時數1782～2064h，年均降水量1784～2367mm。該區(qū)氣溫適宜、地勢較低、坡度較小、水分充足，可滿足綠橙生長發(fā)育所需大部分條件。但是，因?？?、瓊海等地城市化較快，故綠橙實際種植面積較少。由圖7b可知，未來（2020-2049年）綠橙種植氣候最適宜區(qū)與現階段相比，區(qū)劃面積無較大變化，但分布區(qū)域總體呈下移趨勢，主要分布于瓊中縣、屯昌縣和保亭縣大部分地區(qū)、萬寧市西部、白沙縣中部及東部等地，總面積約0.93×104km2。

圖6 海南島過去40a(a)和未來RCP4.5情景下(b)綠橙種植降水適宜度空間分布

圖7 海南島過去40a(a)和未來RCP4.5情景下(b)綠橙種植氣候適宜性區(qū)劃圖

（2）適宜區(qū)

由圖7a可知，過去40a適宜區(qū)主要分布于海南島東部區(qū)域及北部大部分地區(qū)，包括?？谑兄胁?、文昌市中部及東部沿海大部分區(qū)域、定安縣中部地區(qū)、瓊海市中部及東部沿海地區(qū)，澄邁縣西部和中部、臨高縣、儋州市大部分地區(qū)、白沙縣北部區(qū)域、東方市東部小部分地區(qū)、樂東縣東部小部分區(qū)域、三亞市北部小部分地區(qū)、陵水縣北部地區(qū)等，面積約1.83×104km2。該區(qū)年均氣溫21.3～25.6℃，≥10℃年均積溫＞7600℃·d，年均日照時數1624～2271h，年均降水量1309～2465mm。但三亞市、陵水縣、保亭縣等地位于赤道北緣，氣候炎熱，且三亞城市化發(fā)展較快，故不適宜綠橙種植。由圖7b可知，未來（2020-2049年）綠橙種植氣候適宜區(qū)與現階段相比，區(qū)劃面積明顯減少，分布區(qū)域減少，主要分布于?？诤臀牟写蟛糠值貐^(qū)、定安縣、瓊海市中部及東部、萬寧市東部地區(qū)、三亞市北部及陵水縣等地，總面積約1.52×104km2。

（3）次適宜區(qū)

由圖7a可知，過去40a次適宜區(qū)主要分布于海南島西部沿海及中西部海拔較高區(qū)域，包括昌江縣西部和北部小部分區(qū)域、東方市西部沿海大部分區(qū)域及中部、樂東縣沿海西部和南部區(qū)域、五指山市邊緣及中部、保亭縣北部及西部零星區(qū)域、萬寧市東部沿海小部分區(qū)域、瓊中縣邊緣小部分地區(qū)等，面積約0.51×104km2。該區(qū)年均氣溫19.8～26.7℃，≥10℃年均積溫分布不穩(wěn)，年均日照時數1394～2542h，降水量較少。該區(qū)降水較少，平均氣溫較低，部分地區(qū)氣溫偏高，部分地區(qū)日照不足。但綜合考慮，如管理得當，也可種植綠橙。由圖7b可知，未來（2020-2049年）綠橙種植氣候次適宜區(qū)與現階段相比，區(qū)劃面積明顯增加，其分布區(qū)域從海南島西部、南部沿海地區(qū)遷移至西部、北部沿海地區(qū)，主要分布于臨高縣和儋州市大部分地區(qū)、昌江縣北部、東方市西部、樂東縣南部等地，總面積約0.69×104km2。

（4）不適宜區(qū)

由圖7a可知，過去40a不適宜區(qū)主要分布于海南島中部海拔較高地區(qū)，包括昌江縣南部部分地區(qū)、瓊中縣東部及西部較小區(qū)域、陵水縣北部小部分區(qū)域、樂東縣西部沿海部分地區(qū)、五指山市中部及北部小部分地區(qū)、白沙縣中部及西部零星區(qū)域等，面積約0.17×104km2。該區(qū)海拔高度均400m以上，地形復雜，年均氣溫＜21℃，降水量較高，不適宜綠橙生長發(fā)育。由圖7b可知，未來（2020-2049年）綠橙種植氣候不適宜區(qū)與現階段相比，區(qū)劃面積和區(qū)域無較大變化，主要分布于東方市西部沿海地區(qū)、昌江縣南部地區(qū)、白沙縣南部小部分地區(qū)、五指山市北部小部分地區(qū)及瓊中縣南部零星地區(qū)等地，總面積約0.24×104km2。

3 結論與討論

3.1 結論

海南島綠橙種植最適宜區(qū)，主要位于島內中部地區(qū)，氣候適宜度為0.9～0.1；適宜區(qū)主要位于東部及北部大部地區(qū)，氣候適宜度為0.9～0.7；次適宜區(qū)主要分布于西部沿海及中西部海拔較高區(qū)域，氣候適宜度為0.4～0.7；不適宜區(qū)主要分布于中部海拔較高區(qū)域，氣候適宜度為0～0.4。

未來氣候變化情景下（RCP 4.5）2020-2049年氣溫、降水分布有較大變化，溫度適宜區(qū)由四周向中部逐漸縮小，降水適宜區(qū)由東部逐漸遷移至中部地區(qū)。未來綠橙種植氣候最適宜區(qū)域，總體呈下移趨勢；適宜區(qū)域減少；次適宜區(qū)域增加，從海南島西部、南部沿海地區(qū)遷移至海南島西部、北部沿海地區(qū)；不適宜區(qū)域無較大變化。

3.2 討論

選取區(qū)劃指標是進行果樹氣候適宜性區(qū)劃的關鍵[35]，但目前有關綠橙種植區(qū)劃指標尚未有統(tǒng)一標準。沈兆敏等[9]利用氣溫、≥10℃積溫、極端最低氣溫等氣溫指標研究了全國柑橘種植氣候區(qū)劃；趙彤[12]利用年平均氣溫、年總輻射量、空氣相對濕度、年日照時數開展了重慶不同柑橘品種種植分區(qū)研究；陳小敏等[16]利用年平均溫度、年日照時數、年降水量、海拔高度、坡度等指標劃分了海南瓊中綠橙適宜種植區(qū)；劉聰等[22]利用柑橘發(fā)育期溫度、降水、日照適宜度模型對南方柑橘進行前、后30a適宜性氣候區(qū)劃分析研究。本研究在前人研究的基礎上[21?23]，基于氣候適宜度函數構建海南綠橙溫度、日照、降水適宜度模型，對海南島綠橙進行氣候適宜性精細化區(qū)劃分析。為了區(qū)劃結果的精確性，該研究還結合了等權重求和法、幾何平均法和不等權求和法共同建立綜合氣候適宜度模型，在GIS技術的支持下，采用自然點斷法將海南島綠橙種植氣候適宜性分為最適宜區(qū)、適宜區(qū)、次適宜區(qū)和不適宜區(qū)四個區(qū)域。

目前海南島綠橙種植區(qū)，主要集中于瓊中縣、白沙縣、澄邁縣、儋州縣、臨高縣、?？诘仁锌h，均處于現階段綠橙種植適宜區(qū)至最適宜區(qū)范圍內，表明本研究區(qū)劃結果較可靠。本研究也發(fā)現一些綠橙潛在最適宜或適宜種植區(qū)，包括萬寧市北部小部分地區(qū)、定安縣大部地區(qū)、瓊海市中部及沿海東部地區(qū)等地，這些地區(qū)可結合當地經濟社會及果樹生產適當發(fā)展綠橙種植業(yè)。農業(yè)氣候適宜性區(qū)劃只能從宏觀上對綠橙種植業(yè)提出指導建議，因海南島地形復雜，還需綜合考慮當地實際生產情況，開展小范圍引種試驗，再根據試驗情況進行綠橙規(guī)?；N植。

通過對未來（2020-2049年）海南綠橙種植氣候適宜性變化分析，明確了綠橙種植氣候適宜區(qū)的變化趨勢，對綠橙種植基地選址及綠橙產業(yè)發(fā)展，均具有參考價值。本研究主要考慮了氣候和地理因子，但未涉及氣象災害、品種、土壤等因素。因此，未來研究中可結合氣象災害等因素，進一步開展氣候變化背景下海南島綠橙種植氣候適宜性研究，為農業(yè)可持續(xù)生產提供參考依據。

[1] 海南省統(tǒng)計局.2020年海南省統(tǒng)計年鑒[M].北京:中國統(tǒng)計出版社,2020.

Hainan Provincial Bureau of Statistics.Hainan statistical yearbook 2020[M].Beijing:China Statistics Press,2020.(in Chinese)

[2] 張業(yè)忠.瓊中縣種植綠橙的氣候條件分析[J].熱帶農業(yè)科學,2017,37(12):44-47.

Zhang Y Z.Analysis of climatic conditions for green orange plantedin Qiongzhong,Hainan[J].Chinese Journal of Tropical Agriculture,2017,37(12):44-47.(in Chinese)

[3] 吳海虹,廖孝文,項秋蘭,等.海南瓊中綠橙產業(yè)發(fā)展現狀及對策[J].中國園藝文摘,2015,31(3):55-58.

Wu H H,Liao X W,Xiang Q L,et al.Development status and countermeasures of green orange industry in Qiongzhong, Hainan[J].Chinese Horticultural Abstracts,2015,31(3):55- 58.(in Chinese)

[4] IPCC.Climate change 2013:the physical science basis[M]. Cambridge:Cambridge University Press,2013.

[5] 王玉潔,周波濤,任玉玉,等.全球氣候變化對我國氣候安全影響的思考[J].應用氣象學報,2016,27(6):750-758.

Wang Y J,Zhou B T,Ren Y Y,et al.Impacts of global climate change on china's climate security[J].Journal of Applied Meteorological Science,2016,27(6):750-758.(in Chinese)

[6] Aparecido L E d O,Moraes J R d S C d,Rolim G d S,et al.Neural networks in spatialization of meteorological elements and their application in the climatic agricultural zoning of bamboo[J].Int J Biometeorol,2018,62:1955- 1962.

[7] de Oliveira Aparecido L E,da Silva Cabral de Moraes J R,de Meneses K C,et al.Agricultural zoning as tool for expansion of cassava in climate change scenarios[J].Theor Appl Climatol,2020,142:1085-1095.

[8] 沈兆敏,張伯雍,何天富,等.我國柑桔的生態(tài)適宜性區(qū)劃研究[J].中國農業(yè)科學,1984(2):1-7.

Shen Z M,Zhang B Y,He T F,et al.Study on ecological suitability zoning of citrus in China[J].Scientia Agricultura Sinica,1984(2): 1-7.(in Chinese)

[9] 沈兆敏.中國柑桔區(qū)劃與柑桔良種[M].北京:中國農業(yè)科技出版社,1988.

Shen Z M.Citrus regionalization and improved citrus varieties in China[M].Beijing:China Agricultural Science and Technology Press,1988.(in Chinese)

[10] 池再香,白慧,彭興德,等.黔東南山區(qū)椪柑生態(tài)氣候適宜性分析及區(qū)劃[J].西南農業(yè)學報,2008,4(4):1063-1068.

Chi Z X,Bai H,Peng X D,et al.Analysis of ecological climate adaptability of citrus and its division in mountainous areas of southeast Guizhou[J].Southwest China Journal of Agriculture Sciences,2008,4(4):1063- 1068.(in Chinese)

[11] 林紅,楊麗霞,張玉芳,等.基于GIS的大英縣甜橙種植農業(yè)氣候區(qū)劃研究[J].中國農學通報,2016,32(32):164-168.

Lin H,Yang L X,Zhang Y F,et al.Agricultural climate division for orange-growing based on GIS in Daying[J]. Chinese Agricultural Science Bulletin,2016,32(32):164- 168.(in Chinese)

[12] 趙彤.基于GIS的重慶市柑橘農業(yè)氣候區(qū)劃[D].重慶:重慶師范大學,2018.

Zhao T.The study on agricultural climate regionalization of citrus in Chongqing Province based on GIS[D].Chongqing: Chongqing Normal University,2018.(in Chinese)

[13] 王迅,熊博,李啟權,等.四川省晚熟柑橘生態(tài)氣候適宜性區(qū)劃研究[J].中國農學通報,2021,37(13):94-101.

Wang X,Xiong B,Li Q Q,et al.Analyzing and zoning of eco-climate suitability of late-maturing citrus in Sichuan[J]. Chinese Agricultural Science Bulletin,2021,37(13):94- 101.(in Chinese)

[14] 楊愛萍,杜筱玲,王保生,等.江西省多氣象要素的柑橘凍害指標[J].應用氣象學報,2013,24(2):248-256.

Yang A P,Du X L,Wang B S,et al.Citrus freezing injury index of multi meteorological elements in Jiangxi province[J]. Journal of Applied Meteorological Science,2013,24(2): 248-256.(in Chinese)

[15] 謝遠玉,王培娟,朱凌金,等.基于氣象因子的贛南臍橙氣候品質指標評價模型[J].生態(tài)學雜志,2019,38(7):2265- 2274.

Xie Y Y,Wang P J,Zhu L J,et al.Climate quality evaluation model for navel orange in Ganzhou[J].Chinese Journal of Ecology,2019,38(7):2265-2274.(in Chinese)

[16] 陳小敏,陳匯林,張業(yè)忠,等.基于GIS的瓊中縣綠橙種植氣候區(qū)劃研究[J].熱帶農業(yè)科學,2012,32(12):100-102.

Chen X M,Chen H L,Zhang Y Z,et al.Agro-climatic division for growing Qiongzhong green orange cultivar by GIS technology[J].Chinese Journal of Tropical Agriculture, 2012,32 (12):100-102.(in Chinese)

[17] 劉少軍,周廣勝,房世波.中國橡膠樹種植氣候適宜性區(qū)劃[J].中國農業(yè)科學,2015,48(12):2335-2345.

Liu S J,Zhou G S,Fang S B.Climatic suitability regionalization of rubber plantation in China[J].Scientia Agricultura Sinica,2015,48(12):2335-2345.(in Chinese)

[18] 鄒雨伽,楊柳,張秀瓊,等.基于GIS的茶葉生產適宜性區(qū)劃研究:以旺蒼縣為例[J].西南大學學報(自然科學版),2021, 43(3):53-59.

Zou Y J,Yang L,Zhang X Q,et al.Study on Suitability zoning of tea production based on GIS:take Wangcang County as an example[J].Journal of Southwest University (Natural Science Edition),2021,43(3):53-59.(in Chinese)

[19] 張亞杰,陳升孛,張京紅,等.屯昌縣黃金蜜柚精細化農業(yè)氣候區(qū)劃[J].中國農業(yè)資源與區(qū)劃,2021,42(3):167-172.

Zhang Y J,Chen S B,Zhang J H,et al.Fine agroclimatic regionalization of golden pomelo in Tunchang County[J].Chinese Journal of Agricultural Resources and Regional Planning,2021,42(3):167-172.(in Chinese)

[20] 李丹陽,黃霖霄,楊毅.柑桔黃龍病菌濃度對瓊中綠橙果實品質的影響[J].中國南方果樹,2021,50(2):26-29.

Li D Y,Huang L X,Yang Y.Effect of citrus yellow dragon disease concentration on fruit quality of green oranges in Qiongzhong[J].Southern China Fruit Trees,2021,50(2): 26-29.(in Chinese)

[21] 劉武,莫家堯,李政,等.廣西柑橘氣候適宜度模型[J].中國農學通報,2021,37(25):109-114.

Liu W,Mo J Y,Li Z,et al.Climate suitability model of citrus in Guangxi[J].Chinese Agricultural Science Bulletin,2021, 37(25):109-114.(in Chinese)

[22] 劉聰,李凱偉,張繼權,等.基于氣候適宜度的南方柑橘種植精細化氣候區(qū)劃[J].應用氣象學報,2021,32(4):421-431.

Liu C,Li K W,Zhang J Q,et al.Refined climatic zoning for citrus cultivation in Southern China based on climate suitability[J].Journal of Applied Meteorological Science, 2021,32(4):421-431.(in Chinese)

[23] 李麗純,周廣勝,陳惠,等.臺灣青棗在福建主產區(qū)的氣候適宜度[J].應用生態(tài)學報,2020,31(4):1250-1258.

Li L C,Zhou G S,Chen H,et al.Climate suitability ofin the main planting areas of Fujian province,China[J].Chinese Journal of Applied Ecology, 2020,31(4):1250-1258.

[24] 林正雨,陳強,鄧良基,等.中國柑橘生產空間變遷及其驅動因素[J].熱帶地理,2021,41(2):374-387.

Lin Z Y,Chen Q,Deng L J,et al.Spatial pattern changes and driving factors of citrus production in China[J].Tropical Geography,2021,41(2):374-387.(in Chinese)

[25] 胡蕓蕓,楊世琦,陳志軍,等.基于GIS技術的重慶市巴南區(qū)柑橘種植氣候區(qū)劃[J].重慶師范大學學報(自然科學版),2010,27(6):79-82.

Hu Y Y,Yang S Q,Chen Z J,et al.Climatic regionalization of citrus planting in Banan district of Chongqing based on GIS technology[J].Journal of Chongqing Normal University (Natural Science),2010,27(6):79-82.(in Chinese)

[26] 胡正月.柑橘優(yōu)質豐產栽培300問[M].北京:金盾出版社,2008.

Hu Z Y.Citrus quality and high yield cultivation 300 questions[M].Beijing:Jindun Publishing House,2008.(in Chinese)

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

Chen X M,Li W G,Chen H L,et al.Establishment and application of climate suitability evaluation index for rubber tapping in Hainan island:take Danzhou city as an example[J].Jiangsu Agricultural Sciences,2019,47(15):278- 281.(in Chinese)

[28] 蔣高明,常杰,高玉葆,等.植物生理生態(tài)學[M].北京:高等教育出版社,2004:150-154.

Jiang G M,Chang J,Gao Y B,et al.Plant physiological ecology[M].Beijing:Higher Education Press,2004:150-154. (in Chinese)

[29] 邱美娟,劉布春,劉園,等.中國北方蘋果種植需水特征及降水適宜性[J].應用氣象學報,2021,32(2):175-187.

Qiu M J,Liu B C,Liu Y,et al.Water requirement and precipitation suitability of apple planting in Northern China[J].Journal of Applied Meteorological Science,2021, 32(2):175-187.(in Chinese)

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

Hou Y Y, Zhang Y H,Wang L Y,et al.Climatic suitability model for spring maize in Northeast China[J].Chinese Journal of Applied Ecology,2013,24(11):3207-3212.(in Chinese)

[31] Pan Y X. GIS-based niche suitability model of rice breeding and cultivation and evaluation of multi-suitability of cultivated land[J]. Arabian Journal of Geosciences, 2021, 14(11): 932.

[32] 肖秀珠,劉佳明,孫朝鋒,等.福建省藍莓種植氣候適宜性區(qū)劃[J].中國農業(yè)氣象,2020,41(8):520-528.

Xiao X Z,Liu J M,Sun C F,et al.Climate suitability zoning of blueberry in Fujian province[J].Chinese Journal of Agrometeorology,2020,41(8):520-528.(in Chinese)

[33] 王苗,劉敏,任永建.基于高分辨率模擬數據RCP4.5情景下的華中區(qū)域氣候變化預估[J].氣象與環(huán)境學報,2021, 37(3):65-72.

Wang M,Liu M,Ren Y J.Projection of climate change over Centre China by high-resolution simulation data under RCP4.5 scenario[J].Journal of Meteorology and Environment, 2021,37(3):65-72.(in Chinese)

[34] 左琛,宋鋒惠,史彥江,等.新疆平歐雜種榛種植氣候區(qū)劃[J].中國農業(yè)氣象,2020,41(3):146-155.

Zuo C,Song F H,Shi Y J,et al.Climate regionalization of hybrid hazel planting in Xinjiang[J].Chinese Journal of Agrometeorology,2020,41(3):146-155.(in Chinese)

[35] 薛麗芳,申雙和,王春林.基于GIS的廣東香蕉種植氣候適宜性區(qū)劃[J].中國農業(yè)氣象,2010,31(4):575-581.

Xue L F,Shen S H,Wang C L.Climatic suitability and regionalization of banana planting in Guangdong based on GIS[J].Chinese Journal of Agrometeorology,2010,31(4): 575-581.(in Chinese)

Climate Suitability Analysis of Green Orange Cultivation in Hainan Island under Future Climate Change Scenarios

CHEN Yan-xi1, 2, 3, LOU Yun-sheng1, 3, REN Li-xuan4, 5, 6, SU Lei7, TANG Li-ling7, YANG Jian-zhou7

(1. Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science and Technology, Nanjing 210044, China; 2. Hainan Meteorological Detection Center, Haikou 570203; 3. Jiangsu Key Laboratory of Agricultural Meteorology/School of Applied Meteorology, Nanjing University of Information Science and Technology, Nanjing 210044; 4. School of Resources and Environmental Science, Nanjing Agricultural University, Nanjing 210095; 5. Key Laboratory of Plant Nutrition and Fertilization in the Low-Middle Reaches of the Yangtze River, Ministry of Agriculture and Rural Affairs, Nanjing 210095; 6. Jiangsu Collaborative Innovation Center for Solid Organic Waste Resource Utilization, Nanjing 210095; 7. Institute of Geophysical and Geochemical Exploration, Chinese Academy of Geological Sciences, Langfang 065000)

Green orange is a characteristic fruit and plays an important role in local agricultural production in Hainan island. However, few researches have been available regarding the climatic suitability of green orange planting in the island. Based on the last 40a (1980-2019) daily meteorological data, topographical and other factors from 19 meteorological stations in the main island of Hainan island, a spatial analysis model is established to comprehensively determine the climate suitability index for green orange planting in the island. With climate suitability model, this study constructed the suitability models of temperature, sunshine and precipitation, as well as comprehensive suitability model for green orange planting. The methods of geographic information system (GIS) and natural point break were used in finely zoning green orange planting climate suitability. With the model simulation data of RCP4.5 climate scenario, this study assessed the trend of green orange planting suitability zone in the next 30y (2020-2049) under future climate change scenario. The results show that, the most suitable area for green orange planting is mainly distributed in the central region, with the subtotal area being 0.87×104km2, and the climate suitability index ranging from 0.9-1.0; the suitable area is mainly distributed in the eastern local region and most of the central and western regions, with the subtotal area being 1.83×104km2, and the climate suitability index being 0.7-0.9; the subsuitable area is mainly located in the western coastal, central and western high altitude regions, with subtotal area being 0.51×104km2, and the climate suitability index varying from 0.4-0.7; the unsuitable area is mainly distributed in the central mountainous region, with the subtotal area being 0.17×104km2, and the climate suitability index ranging from 0-0.4. Under the future climate scenario, the area with suitable temperature and precipitation will have large change, namely, the area with suitable temperature will gradually shrink from the surrounding to the central regions, the area with suitable precipitation will gradually move from the eastern to the central regions. The most climatic suitable area for future green orange planting will mainly distribute in the most parts of Qiongzhong, Tunchang and Baoting counties, western parts of Wanning, central and eastern parts of Baisha county in the island.

Green orange; Suitability; Climate zones; GIS; Future climate scenarios

10.3969/j.issn.1000-6362.2022.10.002

陳彥希,婁運生,任麗軒,等.未來氣候變化情景下海南島綠橙種植氣候適宜性分析[J].中國農業(yè)氣象,2022,43(10):786-797

2021?12?15

中國地質調查局地質調查項目（DD20190305）；國家自然科學基金項目（41875177）

婁運生，教授，主要從事氣候變化與氣候資源利用、農業(yè)氣象研究，E-mail: yunshlou@163.com；任麗軒，副教授，研究方向農業(yè)資源高效利用，E-mail: lxren@njau.edu.cn；蘇磊，工程師，研究方向生態(tài)地質調查，E-mail: sulei@igge.cn

陳彥希，E-mail:yumcyC@outlook.com