• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2022-10-25 06:20:44陳彥希婁運(yùn)生任麗軒湯麗玲楊建洲
    中國(guó)農(nóng)業(yè)氣象 2022年10期
    關(guān)鍵詞:日照時(shí)數(shù)海南島區(qū)劃

    陳彥希,婁運(yùn)生,任麗軒,蘇 磊,湯麗玲,楊建洲

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

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

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

    基于海南島主島19個(gè)氣象站點(diǎn)連續(xù)40a(1980-2019年)逐日氣象數(shù)據(jù)及地形等要素,建立空間分析模型,綜合確定海南島綠橙種植氣候適宜區(qū)指標(biāo)。在氣候適宜度模型基礎(chǔ)上構(gòu)建海南島綠橙種植溫度、日照、降水以及綜合適宜度模型,采用地理信息系統(tǒng)(GIS)技術(shù)綜合分析,利用自然點(diǎn)斷法對(duì)綠橙種植氣候適宜性進(jìn)行精細(xì)化區(qū)劃,并利用RCP4.5未來(lái)氣候情景模式模擬數(shù)據(jù)探討未來(lái)30a(2020-2049年)氣候變化情景下海南島綠橙種植適宜區(qū)變化趨勢(shì)。結(jié)果表明:現(xiàn)階段海南島綠橙種植最適宜區(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。未來(lái)氣候變化情景下海南島氣溫、降水分布呈較大變化,溫度適宜區(qū)由四周向中部逐漸縮小,降水適宜區(qū)由東部逐漸遷移至中部地區(qū)。未來(lái)綠橙種植氣候最適宜區(qū)主要分布于瓊中、屯昌和保亭縣大部分地區(qū)、萬(wàn)寧市西部及白沙縣中東部等地。

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

    綠橙是海南島知名特色水果之一,前身為紅江橙,20世紀(jì)80年代由廣東廉江引入海南后,表現(xiàn)出比原產(chǎn)地更強(qiáng)的適應(yīng)性,具有結(jié)實(shí)早、豐產(chǎn)快等特點(diǎn)[1?2]。海南島綠橙已獲農(nóng)業(yè)農(nóng)村部無(wú)公害農(nóng)產(chǎn)品及綠色食品認(rèn)證,是海南優(yōu)質(zhì)名牌農(nóng)產(chǎn)品[3],種植面積約占該省柑橘橙柚類種植面積的36.4%。

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

    1 資料與方法

    1.1 研究區(qū)概況

    海南島位于中國(guó)最南端,下屬18個(gè)縣市,陸地面積約3.39×104km2[17]。該地氣候?yàn)闊釒Ъ撅L(fēng)氣候,年均氣溫23.2~27.1℃,總降水量1009~2367.7mm,日照時(shí)數(shù)1776~2783h,島內(nèi)氣象觀測(cè)站點(diǎn)分布見(jiàn)圖1。

    圖1 海南島主島19個(gè)氣象觀測(cè)站點(diǎn)分布

    1.2 數(shù)據(jù)來(lái)源

    (1)1980-2019年海南島主島19個(gè)氣象站點(diǎn)逐日氣象資料,包括日平均氣溫、日總降水量、日照等氣象要素,來(lái)自國(guó)家氣象數(shù)據(jù)網(wǎng)(http://data. cma.cn/)。

    (2)海南島主島DEM數(shù)字高程模型數(shù)據(jù)。利用DEM數(shù)據(jù)獲取坡度、坡向、海拔高度等數(shù)據(jù)(30m×30m),來(lái)自國(guó)家基礎(chǔ)地理信息中心(http:// www. ngcc.cn/ngcc/)。

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

    (4)海南島未來(lái)氣候變化情景數(shù)據(jù),即RCP4.5排放情景下CMIP5(Coupled Model Intercomparison Project Phase 5)模式模擬數(shù)據(jù),包括2020-2049年溫度數(shù)據(jù)和降水?dāng)?shù)據(jù)來(lái)自http://data.ceda.ac.uk/。

    (5)海南島柑橘橙柚、綠橙農(nóng)產(chǎn)品種植與生產(chǎn)情況相關(guān)數(shù)據(jù)來(lái)自海南省統(tǒng)計(jì)年鑒(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,相對(duì)濕度85.3%,干燥度0.75,降水充足,但季節(jié)降水分布不均,多年平均降水量冬季90mm,春季321mm,夏季721mm,秋季676.8mm;年均總?cè)照諘r(shí)數(shù)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 氣象數(shù)據(jù)空間分析模型建立

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

    式中,Y為氣候區(qū)劃指標(biāo)因子,?、λ、h、β和θ分別為緯度(°)、經(jīng)度(°)、海拔高度(m)、坡向(°)、坡度(°)等地理因子,e為余差項(xiàng),稱為地理殘差。

    1.3.3 氣候適宜度模型建立

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

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

    1.3.3.2 溫度適宜度模型

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

    式中,ST為作物溫度適宜度(℃),t為實(shí)際溫度(℃),t0、tl、th分別為作物生長(zhǎng)最適、下限和上限溫度(℃),具體取值見(jiàn)表2。

    1.3.3.3 日照適宜度模型

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

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

    1.3.3.4 降水適宜度模型

    水分影響綠橙果實(shí)大小、質(zhì)地和風(fēng)味等品質(zhì)指標(biāo)[28],綠橙種植最佳年降水量為2000~2200mm[20]。作物降水適宜度模型為[27?29]

    式中,SP為作物降水適宜度,Sr為作物降水量適宜度,Sd為作物降水日數(shù)適宜度。

    表2 溫度適宜度指標(biāo)

    式中,R為實(shí)際降水量(mm),R0為作物生長(zhǎng)最適降水量(mm),取值2100mm。

    式中,d為實(shí)際年降水日數(shù)(d),dl、dh分別為作物生長(zhǎng)適宜降水日數(shù)的下限和上限(d),取值分別為140d和200d。

    1.3.3.5 綜合適宜度模型

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

    等權(quán)重求和法計(jì)算式為

    幾何平均法計(jì)算式為

    不等權(quán)求和法計(jì)算式為

    式中,S為作物綜合適宜度,WT、WS和WP分別為作物溫度、日照時(shí)數(shù)以及降水的適宜度權(quán)重,均采用層次分析法(AHP)所確定[23]。

    2 結(jié)果與分析

    2.1 未來(lái)海南島綠橙種植關(guān)鍵氣候因子時(shí)空變化

    2.1.1 氣候因子空間分析模型

    利用1980-2019年海南島19個(gè)氣象站點(diǎn)年平均氣溫、年總?cè)照諘r(shí)數(shù)、年總降水量3個(gè)氣象要素,分別與經(jīng)度、緯度、海拔高度、坡向和坡度5個(gè)地理信息因子進(jìn)行回歸分析,構(gòu)建空間分析模型[32?34],結(jié)果見(jiàn)表3。

    由表3可見(jiàn),各模型的相關(guān)系數(shù)均通過(guò)了0.01水平的顯著性檢驗(yàn)。其中,年平均氣溫與緯度、經(jīng)度、海拔高度、坡度和坡向均呈負(fù)相關(guān)關(guān)系;年總?cè)照諘r(shí)數(shù)與緯度、經(jīng)度、坡度和坡向呈負(fù)相關(guān)關(guān)系,與海拔高度呈正相關(guān)關(guān)系;年總降水量與緯度和海拔高度呈負(fù)相關(guān)關(guān)系,與經(jīng)度、坡度和坡向呈負(fù)相關(guān)關(guān)系;年≥10℃積溫與緯度、經(jīng)度、海拔高度、坡度和坡向均呈負(fù)相關(guān)關(guān)系。

    2.1.2 關(guān)鍵氣象因子時(shí)空變化

    2.1.2.1 溫度

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

    2.1.2.2 降水量

    海南島降水強(qiáng)度呈現(xiàn)東高西低、四周向中間漸增趨勢(shì)。從1980-2019年40a平均狀況看(圖3a),中部部分地區(qū)年均總降水量>2400mm,西部沿海地區(qū)年均總降水量<1800mm。從RCP4.5情景下2020-2049年30a平均狀況看(圖3b),年均總降水量>2400mm區(qū)域與過(guò)去40a相比,其范圍由中部部分地區(qū)擴(kuò)至中部大部分地區(qū),面積增加0.43×103km2,約為原面積的14%;東部大部分區(qū)域年降水量均>2200mm??傮w上,未來(lái)30a海南島平均年降水量較過(guò)去40a均有不同程度增加,且年降水量<1800mm區(qū)域范圍縮小。由于綠橙生長(zhǎng)適宜年平均總降水量在2000~2200mm,因此,綠橙種植降水適宜區(qū)范圍由東部逐漸遷移至中部地區(qū)。

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

    注:**表示相關(guān)系數(shù)通過(guò)0.01 水平顯著性檢驗(yàn)。Φ,λ,h,β,θ分別為緯度坐標(biāo)值、經(jīng)度坐標(biāo)值、海拔高度(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 海南島過(guò)去40a(a)和未來(lái)RCP4.5情景下(b)年平均氣溫空間分布

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

    2.2 未來(lái)海南島綠橙種植氣候適宜性變化

    利用氣候適宜度模型按照綠橙生長(zhǎng)適宜度評(píng)價(jià)指標(biāo)進(jìn)行計(jì)算,在ArcGIS空間分析模塊中生成單個(gè)氣候適宜度柵格圖,即形成氣象要素單因子氣候適宜度空間分布圖(圖4?圖6)。

    (1)溫度適宜度:基于海南島近40a氣象數(shù)據(jù)計(jì)算當(dāng)前溫度適宜度,繪制空間分布圖。由圖4a可見(jiàn),40a內(nèi)全島溫度適宜度總體呈四周高中間低的分布趨勢(shì),溫度適宜度范圍在?1~1,海南島綠橙生長(zhǎng)溫度最適宜區(qū),分布于中部,以瓊中、屯昌及周邊部分地區(qū)為主,該區(qū)年平均氣溫在22~24℃;溫度次適宜區(qū),一部分位于北部及東北部,為???、澄邁、儋州、定安等地,該區(qū)年平均氣溫在24~25℃,另一部分位于白沙、昌江等局部區(qū)域,該地年平均氣溫在21~22℃;溫度不適宜區(qū),一部分位于東部和北部,以東方、三亞、樂(lè)東等地,該區(qū)年平均氣溫>25℃,另一部分位于昌江、五指山、保亭等地,該地年平均氣溫<21℃,不適宜綠橙生長(zhǎng)。

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

    基于海南島未來(lái)30a氣象數(shù)據(jù),計(jì)算得到未來(lái)溫度適宜度,繪制空間分布圖。由圖4b可見(jiàn),全島未來(lái)溫度適宜度總體呈中間高四周低的分布趨勢(shì),未來(lái)溫度適宜度范圍在?1~1,綠橙生長(zhǎng)溫度最適宜區(qū)范圍相較于歷史時(shí)期(1980?2019年)有縮小趨勢(shì),不適宜區(qū)范圍也呈擴(kuò)大趨勢(shì),主要分布于海南島西部及北部沿海地區(qū)。

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

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

    (3)降水適宜度:由海南島過(guò)去40a氣象數(shù)據(jù)計(jì)算當(dāng)前降水適宜度,繪制空間分布圖。由圖6a可見(jiàn),40a內(nèi)全島降水適宜度呈現(xiàn)東高西低趨勢(shì),降水適宜度范圍在0.45~1.00,綠橙生長(zhǎng)發(fā)育降水最適宜區(qū),位于東部及中部局部地區(qū),以???、文昌、定安等縣市為主,年降水量2000~2200mm;次適宜區(qū),一部分位于北部及中部局部地區(qū),包括澄邁、瓊中等地,年降水量1800~2000mm,另一部分位于中部零散區(qū)域,年降水量2200~2400mm;不適宜區(qū),一部分位于西部及西北部,以儋州、昌江、東方、樂(lè)東等市縣及周邊為主,年降水量<1800mm,另一部分位于中部局部地區(qū),年降水量>2400mm。由海南島未來(lái)30a氣象數(shù)據(jù)計(jì)算未來(lái)降水適宜度,繪制空間分布圖,由圖6b可見(jiàn),未來(lái)降水適宜度總體也呈東高西低的分布趨勢(shì),未來(lái)降水適宜度范圍在0.66~1.00,綠橙生長(zhǎng)降水最適宜區(qū)范圍,相較于圖6a有縮小趨勢(shì),由東部逐漸遷移至中部地區(qū),不適宜區(qū)范圍呈擴(kuò)大趨勢(shì),主要分布于海南島西部及北部沿海地區(qū)。

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

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

    (1)最適宜區(qū)

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

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

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

    (2)適宜區(qū)

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

    (3)次適宜區(qū)

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

    (4)不適宜區(qū)

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

    3 結(jié)論與討論

    3.1 結(jié)論

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

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

    3.2 討論

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

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

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

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

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

    [2] 張業(yè)忠.瓊中縣種植綠橙的氣候條件分析[J].熱帶農(nóng)業(yè)科學(xué),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] 吳海虹,廖孝文,項(xiàng)秋蘭,等.海南瓊中綠橙產(chǎn)業(yè)發(fā)展現(xiàn)狀及對(duì)策[J].中國(guó)園藝文摘,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] 王玉潔,周波濤,任玉玉,等.全球氣候變化對(duì)我國(guó)氣候安全影響的思考[J].應(yīng)用氣象學(xué)報(bào),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] 沈兆敏,張伯雍,何天富,等.我國(guó)柑桔的生態(tài)適宜性區(qū)劃研究[J].中國(guó)農(nóng)業(yè)科學(xué),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] 沈兆敏.中國(guó)柑桔區(qū)劃與柑桔良種[M].北京:中國(guó)農(nóng)業(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].西南農(nóng)業(yè)學(xué)報(bào),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的大英縣甜橙種植農(nóng)業(yè)氣候區(qū)劃研究[J].中國(guó)農(nóng)學(xué)通報(bào),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的重慶市柑橘農(nóng)業(yè)氣候區(qū)劃[D].重慶:重慶師范大學(xué),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] 王迅,熊博,李啟權(quán),等.四川省晚熟柑橘生態(tài)氣候適宜性區(qū)劃研究[J].中國(guó)農(nóng)學(xué)通報(bào),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] 楊愛(ài)萍,杜筱玲,王保生,等.江西省多氣象要素的柑橘凍害指標(biāo)[J].應(yīng)用氣象學(xué)報(bào),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] 謝遠(yuǎn)玉,王培娟,朱凌金,等.基于氣象因子的贛南臍橙氣候品質(zhì)指標(biāo)評(píng)價(jià)模型[J].生態(tài)學(xué)雜志,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].熱帶農(nóng)業(yè)科學(xué),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] 劉少軍,周廣勝,房世波.中國(guó)橡膠樹種植氣候適宜性區(qū)劃[J].中國(guó)農(nóng)業(yè)科學(xué),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的茶葉生產(chǎn)適宜性區(qū)劃研究:以旺蒼縣為例[J].西南大學(xué)學(xué)報(bào)(自然科學(xué)版),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] 張亞杰,陳升孛,張京紅,等.屯昌縣黃金蜜柚精細(xì)化農(nóng)業(yè)氣候區(qū)劃[J].中國(guó)農(nóng)業(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] 李丹陽(yáng),黃霖霄,楊毅.柑桔黃龍病菌濃度對(duì)瓊中綠橙果實(shí)品質(zhì)的影響[J].中國(guó)南方果樹,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].中國(guó)農(nóng)學(xué)通報(bào),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] 劉聰,李凱偉,張繼權(quán),等.基于氣候適宜度的南方柑橘種植精細(xì)化氣候區(qū)劃[J].應(yīng)用氣象學(xué)報(bào),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] 李麗純,周廣勝,陳惠,等.臺(tái)灣青棗在福建主產(chǎn)區(qū)的氣候適宜度[J].應(yīng)用生態(tài)學(xué)報(bào),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] 林正雨,陳強(qiáng),鄧良基,等.中國(guó)柑橘生產(chǎn)空間變遷及其驅(qū)動(dòng)因素[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技術(shù)的重慶市巴南區(qū)柑橘種植氣候區(qū)劃[J].重慶師范大學(xué)學(xué)報(bào)(自然科學(xué)版),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)質(zhì)豐產(chǎn)栽培300問(wèn)[M].北京:金盾出版社,2008.

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

    [27] 陳小敏,李偉光,陳匯林,等.海南島橡膠割膠氣候適宜度評(píng)價(jià)指標(biāo)的建立及應(yīng)用:以儋州市為例[J].江蘇農(nóng)業(yè)科學(xué),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)學(xué)[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] 邱美娟,劉布春,劉園,等.中國(guó)北方蘋果種植需水特征及降水適宜性[J].應(yīng)用氣象學(xué)報(bào),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].應(yīng)用生態(tài)學(xué)報(bào),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] 肖秀珠,劉佳明,孫朝鋒,等.福建省藍(lán)莓種植氣候適宜性區(qū)劃[J].中國(guó)農(nóng)業(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] 王苗,劉敏,任永建.基于高分辨率模擬數(shù)據(jù)RCP4.5情景下的華中區(qū)域氣候變化預(yù)估[J].氣象與環(huán)境學(xué)報(bào),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].中國(guó)農(nóng)業(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].中國(guó)農(nóng)業(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

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

    2021?12?15

    中國(guó)地質(zhì)調(diào)查局地質(zhì)調(diào)查項(xiàng)目(DD20190305);國(guó)家自然科學(xué)基金項(xiàng)目(41875177)

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

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

    猜你喜歡
    日照時(shí)數(shù)海南島區(qū)劃
    福州市近70年日照變化趨勢(shì)分析
    西昌近60年日照時(shí)數(shù)的變化特征分析
    1961~2020年曲麻萊縣日照時(shí)數(shù)變化特征
    林芝市雷電風(fēng)險(xiǎn)區(qū)劃
    基于“三線一單”的環(huán)境功能區(qū)劃實(shí)證研究
    川渝傳統(tǒng)民居區(qū)劃的環(huán)境蘊(yùn)意
    1980年~2017年大冶市日照時(shí)數(shù)變化特征分析
    是海南省還是海南島?
    基于GIS的陜西省烤煙氣候生態(tài)適宜性區(qū)劃
    關(guān)于海南島戰(zhàn)役作戰(zhàn)方針的考察
    軍事歷史(1994年6期)1994-08-15 08:56:38
    麻豆一二三区av精品| bbb黄色大片| 亚洲精品乱码久久久v下载方式 | 国产精品av视频在线免费观看| 日韩欧美精品v在线| 国产精品爽爽va在线观看网站| 天堂√8在线中文| 日日干狠狠操夜夜爽| 男人的好看免费观看在线视频| 成年免费大片在线观看| 成人鲁丝片一二三区免费| 色尼玛亚洲综合影院| 午夜日韩欧美国产| 日韩国内少妇激情av| 欧美三级亚洲精品| 欧美一区二区精品小视频在线| 在线观看av片永久免费下载| 老司机在亚洲福利影院| 无限看片的www在线观看| 亚洲无线观看免费| 亚洲精品一卡2卡三卡4卡5卡| 欧美在线黄色| 久久亚洲真实| 天堂影院成人在线观看| 午夜福利在线观看免费完整高清在 | 人妻丰满熟妇av一区二区三区| 又紧又爽又黄一区二区| 免费电影在线观看免费观看| 亚洲最大成人手机在线| 脱女人内裤的视频| 中文字幕av成人在线电影| 国产高清视频在线观看网站| 久久久色成人| 午夜免费男女啪啪视频观看 | 女人被狂操c到高潮| 日韩欧美在线乱码| 亚洲成人精品中文字幕电影| 99国产综合亚洲精品| 少妇人妻精品综合一区二区 | 亚洲五月天丁香| 一个人看的www免费观看视频| 欧美最黄视频在线播放免费| xxx96com| 淫妇啪啪啪对白视频| 国产成人av激情在线播放| 久久精品亚洲精品国产色婷小说| 少妇丰满av| 亚洲av第一区精品v没综合| 国产亚洲欧美98| 久久久久久久久中文| 国产乱人视频| 国产精品女同一区二区软件 | 岛国在线免费视频观看| 亚洲精品456在线播放app | 窝窝影院91人妻| 午夜福利在线观看免费完整高清在 | 十八禁网站免费在线| 国产精华一区二区三区| 美女被艹到高潮喷水动态| 悠悠久久av| 亚洲18禁久久av| 久久性视频一级片| 一级毛片女人18水好多| 精品人妻1区二区| 舔av片在线| 免费在线观看亚洲国产| 天天一区二区日本电影三级| 欧美黄色片欧美黄色片| 好男人在线观看高清免费视频| 欧美成人免费av一区二区三区| 黑人欧美特级aaaaaa片| 12—13女人毛片做爰片一| 日本免费一区二区三区高清不卡| 最近视频中文字幕2019在线8| 国产成人aa在线观看| 丁香六月欧美| 日本五十路高清| 乱人视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 手机成人av网站| 午夜a级毛片| 亚洲黑人精品在线| 国产色爽女视频免费观看| 久久精品91蜜桃| 精品乱码久久久久久99久播| 成人av在线播放网站| 日日夜夜操网爽| 国产精品久久久久久亚洲av鲁大| 欧美激情在线99| 精品99又大又爽又粗少妇毛片 | 久久精品国产综合久久久| 亚洲精品在线观看二区| 99久久无色码亚洲精品果冻| 精品久久久久久成人av| 国产99白浆流出| 99在线人妻在线中文字幕| 一进一出抽搐动态| 天天一区二区日本电影三级| 亚洲av免费在线观看| 日韩亚洲欧美综合| 亚洲无线观看免费| 黄色女人牲交| 久99久视频精品免费| 在线播放无遮挡| 久久久久久九九精品二区国产| 国产精品久久久久久精品电影| 精品久久久久久成人av| 国产午夜精品论理片| 欧美最新免费一区二区三区 | 一级毛片女人18水好多| 夜夜看夜夜爽夜夜摸| 一级黄色大片毛片| 亚洲精品亚洲一区二区| 婷婷六月久久综合丁香| 日日干狠狠操夜夜爽| 97超级碰碰碰精品色视频在线观看| 99精品在免费线老司机午夜| 国产亚洲精品久久久久久毛片| 午夜影院日韩av| 国产精品自产拍在线观看55亚洲| 国产精品嫩草影院av在线观看 | 亚洲专区中文字幕在线| 99久久精品热视频| 国产高清激情床上av| 午夜福利免费观看在线| 亚洲精品粉嫩美女一区| 一区二区三区免费毛片| 老司机福利观看| 波多野结衣巨乳人妻| 淫妇啪啪啪对白视频| 又黄又粗又硬又大视频| 亚洲一区二区三区不卡视频| 小说图片视频综合网站| 又粗又爽又猛毛片免费看| 搡老岳熟女国产| 乱人视频在线观看| 美女高潮的动态| 99热这里只有是精品50| 亚洲一区二区三区不卡视频| 又紧又爽又黄一区二区| 亚洲欧美日韩无卡精品| 床上黄色一级片| 免费搜索国产男女视频| 欧美丝袜亚洲另类 | 女生性感内裤真人,穿戴方法视频| 久久亚洲真实| 亚洲aⅴ乱码一区二区在线播放| 蜜桃亚洲精品一区二区三区| 日日摸夜夜添夜夜添小说| 精品福利观看| 女人高潮潮喷娇喘18禁视频| 97人妻精品一区二区三区麻豆| 噜噜噜噜噜久久久久久91| www日本黄色视频网| 97超视频在线观看视频| 俄罗斯特黄特色一大片| 久久久久国内视频| x7x7x7水蜜桃| 人妻丰满熟妇av一区二区三区| 精品福利观看| 久久久久亚洲av毛片大全| 久久婷婷人人爽人人干人人爱| 亚洲激情在线av| 香蕉丝袜av| 欧美色视频一区免费| 欧美日韩亚洲国产一区二区在线观看| 亚洲中文日韩欧美视频| 我的老师免费观看完整版| 男人舔女人下体高潮全视频| 欧美bdsm另类| 久久伊人香网站| 看免费av毛片| 熟女电影av网| 国产伦精品一区二区三区视频9 | 国产亚洲精品久久久com| 久9热在线精品视频| 老司机福利观看| 中文字幕av在线有码专区| 国产精品98久久久久久宅男小说| 久久久久久九九精品二区国产| 久久久久久久久久黄片| 日韩欧美精品免费久久 | 久久精品国产99精品国产亚洲性色| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久久电影 | av女优亚洲男人天堂| 午夜福利成人在线免费观看| 香蕉久久夜色| 国产三级在线视频| 色吧在线观看| 狠狠狠狠99中文字幕| 色综合欧美亚洲国产小说| 好看av亚洲va欧美ⅴa在| 日本黄大片高清| 岛国视频午夜一区免费看| a在线观看视频网站| 亚洲一区高清亚洲精品| xxxwww97欧美| 美女被艹到高潮喷水动态| 国产一区二区三区视频了| 在线观看一区二区三区| 久久精品91无色码中文字幕| 亚洲av成人精品一区久久| 小说图片视频综合网站| 最新中文字幕久久久久| 男女做爰动态图高潮gif福利片| www.999成人在线观看| 日韩av在线大香蕉| 女人被狂操c到高潮| 一本综合久久免费| 久久人人精品亚洲av| 床上黄色一级片| 男人和女人高潮做爰伦理| 一级毛片女人18水好多| 在线观看av片永久免费下载| 九九久久精品国产亚洲av麻豆| 精品免费久久久久久久清纯| 少妇人妻精品综合一区二区 | 国产精品久久久人人做人人爽| 亚洲五月天丁香| 久久人妻av系列| 在线观看免费视频日本深夜| 久久亚洲精品不卡| 国产一区二区激情短视频| www.www免费av| 欧美大码av| 日韩免费av在线播放| 俄罗斯特黄特色一大片| 欧美日韩一级在线毛片| 国产 一区 欧美 日韩| 亚洲天堂国产精品一区在线| 成年女人毛片免费观看观看9| 精品国内亚洲2022精品成人| 国产美女午夜福利| 亚洲在线观看片| 午夜两性在线视频| 亚洲欧美日韩无卡精品| 久久久精品欧美日韩精品| 国产成人欧美在线观看| 国产老妇女一区| 国产高清激情床上av| 日本 欧美在线| 精品久久久久久久久久免费视频| 日韩国内少妇激情av| 久久人人精品亚洲av| 欧美性猛交╳xxx乱大交人| 欧美成狂野欧美在线观看| 真实男女啪啪啪动态图| 国产亚洲精品久久久久久毛片| 五月玫瑰六月丁香| 女警被强在线播放| 国产激情偷乱视频一区二区| 日韩欧美国产一区二区入口| 18禁美女被吸乳视频| 99久国产av精品| 成人高潮视频无遮挡免费网站| 欧美在线黄色| 欧美乱妇无乱码| 国产高清三级在线| 免费在线观看影片大全网站| 18禁在线播放成人免费| 欧美中文日本在线观看视频| 亚洲欧美日韩东京热| 免费看a级黄色片| 美女高潮的动态| 亚洲欧美日韩东京热| 最好的美女福利视频网| 午夜福利18| 九九热线精品视视频播放| 亚洲成av人片免费观看| 国产精品永久免费网站| 国内精品一区二区在线观看| 亚洲自拍偷在线| 免费大片18禁| 成年免费大片在线观看| 特级一级黄色大片| 男女下面进入的视频免费午夜| 国产精品久久久久久人妻精品电影| 欧洲精品卡2卡3卡4卡5卡区| 久久国产精品影院| 亚洲,欧美精品.| 免费观看的影片在线观看| 99热精品在线国产| 精品一区二区三区人妻视频| 长腿黑丝高跟| 熟女人妻精品中文字幕| a级一级毛片免费在线观看| 国产男靠女视频免费网站| 超碰av人人做人人爽久久 | 尤物成人国产欧美一区二区三区| 欧美三级亚洲精品| 久久久久久大精品| 99久久精品国产亚洲精品| 欧美日韩精品网址| 少妇丰满av| 一边摸一边抽搐一进一小说| 免费一级毛片在线播放高清视频| 午夜福利18| 嫩草影院入口| 日韩大尺度精品在线看网址| 久久国产乱子伦精品免费另类| 色播亚洲综合网| 久久伊人香网站| 日韩国内少妇激情av| 久久亚洲精品不卡| 两个人视频免费观看高清| 成人18禁在线播放| 变态另类成人亚洲欧美熟女| 亚洲精品美女久久久久99蜜臀| 在线免费观看不下载黄p国产 | 成人无遮挡网站| 一级毛片女人18水好多| 村上凉子中文字幕在线| 国产精品日韩av在线免费观看| 国产成人av教育| 一级毛片女人18水好多| 免费看a级黄色片| 国产黄色小视频在线观看| bbb黄色大片| 欧美绝顶高潮抽搐喷水| 亚洲美女黄片视频| 国产老妇女一区| 美女黄网站色视频| 亚洲人与动物交配视频| 亚洲无线观看免费| 亚洲精品亚洲一区二区| 村上凉子中文字幕在线| 欧美在线一区亚洲| 窝窝影院91人妻| 一a级毛片在线观看| 成人av在线播放网站| 亚洲av二区三区四区| 女生性感内裤真人,穿戴方法视频| 精品日产1卡2卡| 国产精品综合久久久久久久免费| 听说在线观看完整版免费高清| 国产精品永久免费网站| 身体一侧抽搐| 99精品欧美一区二区三区四区| 色视频www国产| 欧美三级亚洲精品| 亚洲欧美一区二区三区黑人| 高清在线国产一区| 又粗又爽又猛毛片免费看| 久久久久九九精品影院| 在线播放无遮挡| 久久婷婷人人爽人人干人人爱| 黄色片一级片一级黄色片| 无限看片的www在线观看| avwww免费| 狠狠狠狠99中文字幕| 国产成人啪精品午夜网站| 国产成+人综合+亚洲专区| 一级a爱片免费观看的视频| 日韩av在线大香蕉| 中文在线观看免费www的网站| 国产99白浆流出| 偷拍熟女少妇极品色| 亚洲第一电影网av| 亚洲专区中文字幕在线| av在线天堂中文字幕| 尤物成人国产欧美一区二区三区| 日韩精品青青久久久久久| 久久6这里有精品| 久久精品国产99精品国产亚洲性色| 亚洲国产色片| 成人精品一区二区免费| 欧美日本亚洲视频在线播放| 国产精品久久久久久久电影 | 欧美日韩乱码在线| 老司机福利观看| 一区二区三区激情视频| 欧美性猛交╳xxx乱大交人| 狂野欧美激情性xxxx| 日韩精品中文字幕看吧| 午夜免费观看网址| 免费在线观看亚洲国产| 成人永久免费在线观看视频| 人人妻人人澡欧美一区二区| 99国产精品一区二区蜜桃av| 欧美极品一区二区三区四区| 国产精品亚洲一级av第二区| 综合色av麻豆| 精品99又大又爽又粗少妇毛片 | 亚洲熟妇熟女久久| 一级作爱视频免费观看| 18+在线观看网站| 五月伊人婷婷丁香| 日韩欧美精品免费久久 | 校园春色视频在线观看| 老司机深夜福利视频在线观看| 欧美在线黄色| 十八禁人妻一区二区| 亚洲人成网站高清观看| 成年女人看的毛片在线观看| 男女床上黄色一级片免费看| 天天添夜夜摸| 男女视频在线观看网站免费| 亚洲人成网站在线播| 亚洲av二区三区四区| 噜噜噜噜噜久久久久久91| 久久久久久久亚洲中文字幕 | 亚洲精品亚洲一区二区| 国产精品日韩av在线免费观看| 亚洲人成网站在线播| 午夜福利高清视频| 免费观看的影片在线观看| 757午夜福利合集在线观看| 久久久久久大精品| 国产精品一区二区免费欧美| 国产精品香港三级国产av潘金莲| 欧美色视频一区免费| netflix在线观看网站| АⅤ资源中文在线天堂| 国产视频内射| 老汉色∧v一级毛片| 亚洲av免费在线观看| 别揉我奶头~嗯~啊~动态视频| 人妻丰满熟妇av一区二区三区| 人人妻人人澡欧美一区二区| 午夜亚洲福利在线播放| 久久久久性生活片| 婷婷精品国产亚洲av| 国产精品,欧美在线| 制服丝袜大香蕉在线| 日本黄色片子视频| 亚洲精品一卡2卡三卡4卡5卡| 九九久久精品国产亚洲av麻豆| 国内精品久久久久久久电影| 婷婷亚洲欧美| 九色成人免费人妻av| 狂野欧美激情性xxxx| 757午夜福利合集在线观看| 国产蜜桃级精品一区二区三区| 国产美女午夜福利| 国产精品一区二区三区四区免费观看 | 亚洲不卡免费看| 午夜福利18| 国产亚洲av嫩草精品影院| 亚洲va日本ⅴa欧美va伊人久久| 久久精品夜夜夜夜夜久久蜜豆| 欧美黑人巨大hd| 91九色精品人成在线观看| 美女大奶头视频| 一进一出抽搐动态| 一级黄色大片毛片| 十八禁人妻一区二区| 老汉色∧v一级毛片| 99热6这里只有精品| 亚洲乱码一区二区免费版| 亚洲avbb在线观看| 美女高潮喷水抽搐中文字幕| 日韩欧美在线乱码| 免费人成视频x8x8入口观看| 日本黄大片高清| 国内精品久久久久久久电影| 精华霜和精华液先用哪个| 久久久精品欧美日韩精品| 亚洲一区二区三区不卡视频| 国产精品久久久久久人妻精品电影| 中文亚洲av片在线观看爽| www.www免费av| 91久久精品电影网| 久久草成人影院| 亚洲激情在线av| 99久久精品一区二区三区| 最近视频中文字幕2019在线8| 国内精品一区二区在线观看| 毛片女人毛片| 免费看a级黄色片| 内地一区二区视频在线| 天堂√8在线中文| 国产亚洲精品一区二区www| 国产精品 国内视频| 国产一级毛片七仙女欲春2| 狂野欧美白嫩少妇大欣赏| 亚洲美女黄片视频| 亚洲天堂国产精品一区在线| 亚洲av日韩精品久久久久久密| 日本熟妇午夜| 亚洲成av人片免费观看| 国产野战对白在线观看| 国产精品 欧美亚洲| 男女床上黄色一级片免费看| 亚洲在线自拍视频| 久久香蕉国产精品| 一区二区三区国产精品乱码| av视频在线观看入口| 久久香蕉精品热| 无人区码免费观看不卡| 欧美bdsm另类| 亚洲在线自拍视频| 在线观看美女被高潮喷水网站 | 麻豆久久精品国产亚洲av| 夜夜夜夜夜久久久久| 91av网一区二区| 亚洲av不卡在线观看| 村上凉子中文字幕在线| 欧美中文日本在线观看视频| 高潮久久久久久久久久久不卡| 日韩人妻高清精品专区| 日韩亚洲欧美综合| 成年女人看的毛片在线观看| 一边摸一边抽搐一进一小说| 亚洲av二区三区四区| 夜夜爽天天搞| 欧美日韩精品网址| 欧美zozozo另类| 老司机深夜福利视频在线观看| 免费观看人在逋| 欧美乱色亚洲激情| 香蕉丝袜av| 国产精品乱码一区二三区的特点| 免费人成视频x8x8入口观看| 91久久精品电影网| 亚洲人成网站在线播| 99在线人妻在线中文字幕| 天天一区二区日本电影三级| 欧美色视频一区免费| 亚洲成人中文字幕在线播放| 日韩大尺度精品在线看网址| 岛国在线免费视频观看| 日韩欧美 国产精品| 婷婷丁香在线五月| 亚洲熟妇中文字幕五十中出| 成人高潮视频无遮挡免费网站| 日韩亚洲欧美综合| 国产一区在线观看成人免费| 免费人成在线观看视频色| av视频在线观看入口| 12—13女人毛片做爰片一| 老司机福利观看| АⅤ资源中文在线天堂| 少妇高潮的动态图| 亚洲国产色片| 欧美日本亚洲视频在线播放| 亚洲人成网站高清观看| 51午夜福利影视在线观看| 国产99白浆流出| 久久草成人影院| 国产av不卡久久| 亚洲最大成人中文| 精品欧美国产一区二区三| 五月伊人婷婷丁香| 国产精品野战在线观看| 日本 欧美在线| 久久久久亚洲av毛片大全| 中文字幕人妻熟人妻熟丝袜美 | 亚洲在线自拍视频| 91在线精品国自产拍蜜月 | 蜜桃亚洲精品一区二区三区| 久久精品影院6| 老熟妇乱子伦视频在线观看| 热99re8久久精品国产| 男女做爰动态图高潮gif福利片| 看黄色毛片网站| 免费看十八禁软件| 91久久精品电影网| 色尼玛亚洲综合影院| 嫩草影院精品99| svipshipincom国产片| 国产精品日韩av在线免费观看| 亚洲av二区三区四区| 亚洲欧美激情综合另类| 日韩欧美三级三区| www日本在线高清视频| 熟女人妻精品中文字幕| www.色视频.com| 香蕉久久夜色| 麻豆国产av国片精品| 久久精品夜夜夜夜夜久久蜜豆| 看黄色毛片网站| 嫁个100分男人电影在线观看| 亚洲乱码一区二区免费版| 天堂网av新在线| 免费人成视频x8x8入口观看| 欧美又色又爽又黄视频| 观看免费一级毛片| 两个人视频免费观看高清| 国产乱人视频| 国产av一区在线观看免费| 国内毛片毛片毛片毛片毛片| 悠悠久久av| 日本三级黄在线观看| 久久久精品欧美日韩精品| 国产99白浆流出| 亚洲男人的天堂狠狠| 九九热线精品视视频播放| 一区二区三区激情视频| 人妻久久中文字幕网| 综合色av麻豆| 国产午夜福利久久久久久| 日本撒尿小便嘘嘘汇集6| 精品乱码久久久久久99久播| 岛国在线观看网站| 99热6这里只有精品| 深夜精品福利| 欧美激情久久久久久爽电影| 精品无人区乱码1区二区| 国产成年人精品一区二区| 白带黄色成豆腐渣| 日本撒尿小便嘘嘘汇集6| 天堂动漫精品| 欧美精品啪啪一区二区三区| 可以在线观看的亚洲视频| 国产精品免费一区二区三区在线| 99国产精品一区二区三区| 51国产日韩欧美| 国产精品,欧美在线| 又黄又爽又免费观看的视频| 国产中年淑女户外野战色| 天堂动漫精品| av天堂中文字幕网| 成人鲁丝片一二三区免费| 国产亚洲欧美在线一区二区|