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

    氣候變化下紅茴砂在中國潛在適生區(qū)的最大熵生態(tài)位模型預(yù)測(cè)

    2021-12-23 16:55:29唐夢(mèng)詩袁淑娜余文剛施琦侯天澤吳君楠劉子毓
    熱帶作物學(xué)報(bào) 2021年11期
    關(guān)鍵詞:氣候變化預(yù)測(cè)

    唐夢(mèng)詩 袁淑娜 余文剛 施琦 侯天澤 吳君楠 劉子毓

    摘 ?要:紅茴砂是姜科茴香砂仁屬多年生草本植物,我國僅在海南地區(qū)有自然分布,是具有很高藥用價(jià)值的瀕危植物。為了解紅茴砂在未來氣候變化下的潛在適生區(qū)及其主要環(huán)境影響因子,以期對(duì)紅茴砂的保育提供科學(xué)指導(dǎo),將45條紅茴砂分布位點(diǎn)數(shù)據(jù)與20個(gè)環(huán)境因子相結(jié)合,運(yùn)用最大熵生態(tài)位模型MaxEnt和ArcGIS軟件模擬了當(dāng)前氣候和未來2050年RCP2.6和RCP8.5兩種不同氣候條件下紅茴砂在中國的潛在分布區(qū),并分析影響紅茴砂分布的主導(dǎo)環(huán)境因子。結(jié)果表明:當(dāng)前氣候條件下,紅茴砂適生區(qū)總范圍在18°~32° N、27°~122° E,面積約為1.24×106 km2,主要集中分布在海南、貴州、福建、廣東、廣西、云南等地;影響紅茴砂地理分布的主導(dǎo)環(huán)境因子為最濕月降水量、最暖季度降水量、年均溫度變化范圍、最冷月份最低溫、海拔和最熱季度均溫,累計(jì)貢獻(xiàn)率達(dá)87%;在未來2050年RCP2.6和RCP8.5兩種不同氣候情景下,紅茴砂適生區(qū)域喪失面積均達(dá)到95%以上,潛在分布區(qū)縮小到云南、四川、西藏、貴州和臺(tái)灣這5個(gè)省區(qū),新增區(qū)域主要在西藏東南部和四川中南部,同時(shí)潛在總適生區(qū)和高適生區(qū)的質(zhì)心有向西北方遷移的趨勢(shì)。

    關(guān)鍵詞:紅茴砂;最大熵生態(tài)位模型;中國潛在適生區(qū);氣候變化;預(yù)測(cè)

    中圖分類號(hào):Q948 ? ? ?文獻(xiàn)標(biāo)識(shí)碼:A

    Abstract: Zingiberaceae Etlingera littoralis (J. Koenig) Giseke is a endangered perennial herb with high medicinal val-ue,the natural distribution in China is Hainan. Predicting the potential suitable areas and key environmental factors are critical to the conservation of E. littoralis. Based on the data of 45 current distribution points of E. littoralis in global and 20 environmental factors the key environmental factors and the potential suitable areas under current and future two scenarios (2050 s RCP 2.6, 2050 s RCP 8.5) in China were dtudied using the MaxEnt and ArcGIS, and the potential distribution changes of E. littoralis in China was analyzed, using the SDM toolbox of ArcGIS software, Under the present climatic conditions, E.littoralis is mainly distributed from North latitude 18° to 32°, and East longitude from 27° to 122°, with a total suitable area about 1.24×106 km2. The current potential distribution areas are Hainan, Guizhou, Fujian, Guangdong, Guangxi, Yunnan, Sichuan, Chongqing, Tibet, Hunan, Hubei, Taiwan. Dominant envi-ronmental factors affecting the distribution of E. littorali. were precipitation of the wettest month, precipitation of the warmest quarter, temperature annual range, mean temperature of the coldest month, altitude and mean temperature of the warmest quarter, with a cumulative contribution rate of 87%. Under future climate scenarios of RCP2.6 and RCP8.5 by the 2050 s, the loss of potential suitable areas of E. littoralis would be up to 95%, and the potential distribution areas would be limited to Yunnan, Sichuan, Tibet, Guizhou and Taiwan. The increased areas are in southeastern Tibet and south-central Sichuan. The centers of total suitable areas and most suitable areas would be moved towards northwest.

    Keywords: Etlingera littoralis (J. Koenig) Giseke; MaxEnt; potential suitable area in China; climate change; predic-tion

    DOI: 10.3969/j.issn.1000-2561.2021.11.040

    物種分布區(qū)是物種重要的生態(tài)學(xué)和進(jìn)化學(xué)特征,可作為物種瀕危狀況的一個(gè)指示特征[1],而氣候是影響物種分布的最主要因素之一,氣候變化可以直接或間接改變植物的物種分布范圍,分布格局,物種豐富度和多樣性[2-3],近百年來,在地球表面溫度持續(xù)呈現(xiàn)上升的趨勢(shì)下,我國珍稀瀕危物種表現(xiàn)出脆弱性,氣候變暖使大部分物種空間分布改變,加劇了分布范圍縮減和破碎化而導(dǎo)致的物種滅絕風(fēng)險(xiǎn)[4],提前了解物種地理分布在未來氣候環(huán)境中的變化,有利于指導(dǎo)瀕危物種保育策略的制定。生態(tài)位模型是以生態(tài)位理論為基礎(chǔ)的新興研究領(lǐng)域,通過模型預(yù)測(cè)物種的實(shí)際分布和潛在分布,包括生物氣候分析系統(tǒng)(BIOCLIM)、生態(tài)位因子分析模型(ENFA)、基于規(guī)則集的遺傳算法模型(GARP)、最大熵模型(Maxent)等多個(gè)模型算法[5],其中最大熵生態(tài)位MaxEnt模型近年來被廣泛應(yīng)用于物種保育學(xué)[6-7]、生物入侵學(xué)[8-9]、全球生物變化和生物進(jìn)化學(xué)[10-11]等領(lǐng)域研究,該模型基于最大熵原理利用物種已知的地理分布數(shù)據(jù),進(jìn)行物種生態(tài)位建模和分布分析,在當(dāng)前諸多預(yù)測(cè)模型中,效果突出,即使在物種的分布數(shù)據(jù)較少情況下,仍然能得到滿意的結(jié)果[12-13]。

    紅茴砂[Etlingera littoralis (J. Koenig) Giseke],為姜科茴香砂仁屬多年生草本植物,主要分布在馬來西亞、泰國、越南、印度尼西亞等熱帶地區(qū)和部分亞熱帶地區(qū),生于海拔90~300 m的河道旁、林下、林緣等地,我國僅海南地區(qū)有自然分布。紅茴砂葉和根莖部位含有大量成分豐富的揮發(fā)油,主要為的單萜(monoterpenoids)類和苯基丙烷類化合物,而且紅茴砂提取的揮發(fā)油含有在姜科植物中2種未見報(bào)的二萜醇化合物——8(14), 15-sandaracopimaradien-3b-ol和7,15-isopimaradien- 3b-ol[14],其中(E)-甲基異丁香酚[(E)-methyl isoeugenol)含量最高達(dá)58.1%,研究表明這種天然食品香料具有抗焦慮和抗抑郁作用[15],同時(shí)也對(duì)空腸彎曲桿菌、大腸桿菌彎曲桿菌和弧菌彎曲桿菌等菌株具有明顯抑菌效果[16]。紅茴砂株高2~3 m,花期為3—6月,花色艷紅奪目,枝葉繁茂,在廣東、福建、海南等地區(qū)終年常綠可作為優(yōu)秀的園林背景植物[17]。目前由于紅茴砂生境條件的變遷,其自然種群的延續(xù)受到威脅,現(xiàn)已被列入《中國高等植物受威脅物種名錄》,受威脅等級(jí)為EN[18],但紅茴砂當(dāng)前的研究多集中在其揮發(fā)精油的成分與活性分析[14, 19]、內(nèi)生真菌[20]和紅茴砂種子萌發(fā)、傳粉機(jī)制等方面[21-22],而未有影響其地理分布的主要限制性氣候因子的研究,同時(shí)未來不同氣候變化情景下的紅茴砂潛在分布格局的變遷也尚不明確。

    在全球氣候持續(xù)變暖背景下,本研究運(yùn)用最大熵生態(tài)位模型MaxEnt和ArcGIS軟件,分析紅茴砂目前地理分布區(qū)的氣候特點(diǎn),預(yù)測(cè)紅茴砂當(dāng)前的分布范圍及未來潛在適生區(qū)空間分布格局的變遷,確定影響紅茴砂植物分布的主要環(huán)境因子,旨在為紅茴砂的保育及利用提供支撐。

    1 ?材料與方法

    1.1 ?物種分布點(diǎn)信息來源及篩選

    通過實(shí)地考察,查閱公開發(fā)表的學(xué)術(shù)論文,檢索中國數(shù)字植物標(biāo)本館、NSII-中國國家標(biāo)本資源平臺(tái)和世界生物多樣性信息機(jī)構(gòu)網(wǎng)(GBIF,http://www.gbif.org/)獲得紅茴砂分布數(shù)據(jù),部分只有采集地而無具體經(jīng)緯度信息的分布點(diǎn)結(jié)合GoogleEarth查找相應(yīng)的經(jīng)緯度信息。

    采用緩沖區(qū)分析法對(duì)紅茴砂分布點(diǎn)進(jìn)行校對(duì)、篩選,排除過擬合模擬對(duì)研究結(jié)果的影響。方法:因本研究環(huán)境變量空間分辨率為2.5 minutes(約4.5 km),所以設(shè)緩沖區(qū)為3 km,當(dāng)2個(gè)分布點(diǎn)處于同一個(gè)緩沖區(qū)類,則僅保留其中的1個(gè)點(diǎn)[23]。最終獲得45個(gè)紅茴砂分布點(diǎn)進(jìn)行分析。最后根據(jù)Maxent 軟件要求將分布點(diǎn)數(shù)據(jù)生成CSV格式文件。

    1.2 ?環(huán)境變量和地圖數(shù)據(jù)來源

    本研究采用的地形變量和氣候變量作為本研究的環(huán)境變量,均從世界氣候數(shù)據(jù)庫World-Clim (http://www.worldclim.org/)中獲取,包括3份全球的19個(gè)生物氣候變量圖層 (1970—2000年、2050年RCP2.6和2050年RCP8.5)以及全球海拔變量圖層,選取的空間分辨率均為2.5 minutes。利用2050年RCP2.6和RCP8.5兩個(gè)氣候變化情景,模擬2050年CO2濃度上升對(duì)未來氣候影響的最低和最高情況[24]。本研究繪圖使用底圖數(shù)據(jù)均來自國家基礎(chǔ)地理信息系統(tǒng)網(wǎng)站(http://nfgis. nsdi.gov.cn)所提供的中國行政區(qū)劃圖和標(biāo)準(zhǔn)地圖。

    1.3 ?分析軟件信息

    本研究所用軟件為 Maxent 3.4.1 版本的最大熵模型預(yù)測(cè)軟件(http://Biodiversityinformatics. amnh.org/open_source/maxent/)和ArcGis 10.4軟件(http://desktop.arcgis.com /zh-cn/)。

    1.4 ?紅茴砂適生區(qū)最大熵模型預(yù)測(cè)

    將CSV格式的紅茴砂植物分布數(shù)據(jù)分別與ASCII格式的當(dāng)前時(shí)段1970—2000年及未來2050年RCP2.6和RCP8.5情景下的20個(gè)環(huán)境因子數(shù)據(jù)在MaxEnt軟件中加載,設(shè)置25%的分布點(diǎn)作為測(cè)試集(testing data),75%的分布點(diǎn)作為訓(xùn)練集(training data),利用刀切法(Jackknife)分析每個(gè)氣候因子在預(yù)測(cè)中對(duì)模型的貢獻(xiàn)率。并用ArcGis軟件對(duì)MaxEnt 軟件輸出結(jié)果進(jìn)行可視化表達(dá)和重分類。采用自然間斷點(diǎn)分級(jí)法將其劃分為4個(gè)適生等級(jí)[25],分別是非適生區(qū)(0~0.08)、低適生區(qū)(0.08~0.26)、中適生區(qū)(0.26~0.49)和高適生區(qū)(0.49~0.98)。

    1.5 ?紅茴砂不同時(shí)期適生區(qū)變化分析

    在ArcGis軟件中利用SDMtoolbox工具對(duì)MaxEnt軟件輸出的3份紅茴砂適生區(qū)預(yù)測(cè)圖進(jìn)行計(jì)算當(dāng)前和未來紅茴砂適生區(qū)收縮,擴(kuò)張等變化的面積[26]。

    2 ?結(jié)果與分析

    2.1 ?MaxEnt模型預(yù)測(cè)精度評(píng)估

    本研究基于紅茴砂實(shí)際分布數(shù)據(jù)和20個(gè)綜合性氣候指標(biāo),運(yùn)用最大熵模型進(jìn)行建模,采用ROC曲線分析法來評(píng)價(jià)分布預(yù)測(cè)模型。當(dāng)AUC值大于0.9時(shí)預(yù)測(cè)精度和診斷價(jià)值較高[27]。因此,計(jì)算當(dāng)前紅茴砂潛在分布區(qū)域時(shí)重復(fù)運(yùn)行了10次建模,最終得到平均訓(xùn)練集的AUC值為0.982(圖1)。而在未來2050年RCP2.6和RCP8.5情景下測(cè)試圖的訓(xùn)練集和驗(yàn)證集的AUC值均大于0.98(圖1)說明模型預(yù)測(cè)效果比較準(zhǔn)確,可用于紅茴砂在中國地區(qū)分布情況與氣候關(guān)系的研究。

    2.2 ?影響紅茴砂分布的主要?dú)夂蛞蜃?/p>

    在Max Ent模型中評(píng)估20個(gè)環(huán)境因子在10次重復(fù)預(yù)測(cè)中的重要性,對(duì)模型預(yù)測(cè)貢獻(xiàn)率較高的變量依次為Bio16、Bio18、Bio2、Bio4、Bio6、Altitude、Bio11、Bio10 等個(gè)8因子,累計(jì)貢獻(xiàn)率達(dá)87%。而刀切檢驗(yàn)獲得的模型正則化訓(xùn)練增益最高的8個(gè)環(huán)境變量如圖2依次為Bio7、Bio16、Bio13、Bio12、Bio6、Bio18、Bio4。同時(shí)對(duì)這20個(gè)氣候變量對(duì)其進(jìn)行Pearson相關(guān)性分析,貢獻(xiàn)率和增益較高的變量中存在Bio13和Bio16、Bio7和Bio4、Bio6和Bio11這3組氣候變量相關(guān)性>±0.9,所以每組中選取一個(gè)對(duì)紅茴砂生物意義更高的作為其關(guān)鍵影響因子,綜合分析,篩選出Bio13、Bio18、Bio7、Bio6、altitude、Bio10等6個(gè)變量,為影響紅茴砂地理分布的主要?dú)夂蜃兞?,其貢獻(xiàn)率和主要閾值范圍如表1所示。

    2.3 ?當(dāng)前情景下紅茴砂在中國的潛在分布

    在Maxent模型中,紅茴砂當(dāng)前在我國的潛在分布區(qū)域在海南、云南、四川、臺(tái)灣、廣西、福建、廣東、貴州、重慶的大部分地區(qū),湖南和西藏也有小部分地區(qū)適宜紅茴砂的生長(zhǎng)(表2)。其中,高適生區(qū)(存在概率0.49~0.98)面積為157 068.28 km2,占國土面積1.63%,集中在海南、云南東南部、廣西西部,四川南部小部分地區(qū);中適生區(qū)(存在概率0.26~0.49)面積為313 544.22 km2,占國土面積3.26%,包括云南、海南、廣西、貴州、四川、廣東、西藏和臺(tái)灣等地;低適生區(qū)(存在概率0.08~ 0.26)面積為769 227.15 km2,占國土面積7.99%,分布在廣東、廣西、福建、云南、四川、重慶、西藏、湖南、貴州和臺(tái)灣等地。當(dāng)前情景下,紅茴砂在我國的潛在適生區(qū)總面積約為1.2×106 km2。

    2.4 ?在未來氣候條件下紅茴砂的地理分布格局變遷

    在未來2050年RCP2.6和RCP8.5兩種氣候模型下,紅茴砂適宜分布界線均從華南地區(qū)向西南地區(qū)遷移,但適生區(qū)總體面積喪失極大,如表3所示2種未來氣候模式下紅茴砂適生區(qū)域喪失面積均達(dá)到95%以上,華南地區(qū)(海南、廣東、廣西)、華中地區(qū)(湖南)、西南地區(qū)(重慶、四川東部、貴州中部和東部、云南南部)和華東地區(qū)(福建)基本喪失全部分布,新增區(qū)域主要在西藏東南部和四川中南部。潛在分布區(qū)域減少為云南、四川、西藏、貴州和臺(tái)灣5個(gè)地區(qū),主要集中在滇中北部和川南,進(jìn)一步收縮在橫斷山脈、云貴高原等高海拔地區(qū),呈現(xiàn)向西藏延伸分布的現(xiàn)象,臺(tái)灣地區(qū)適生區(qū)也向中央山脈高海拔地區(qū)聚攏。當(dāng)前紅茴砂潛在總適生區(qū)質(zhì)心分布在貴州南部獨(dú)山縣內(nèi),高適生區(qū)質(zhì)心在紅河哈尼族彝族自治州,PCR2.6和PCR8.5兩種未來情景模式紅茴砂潛在總適生區(qū)質(zhì)心分別在云南省楚雄彝族自治州和云南迪慶藏族自治州,但高適生區(qū)質(zhì)心均在云南貢山獨(dú)龍族怒族自治區(qū)。不論是紅茴砂潛在總適生區(qū)質(zhì)心還是高適區(qū)適生區(qū)質(zhì)心都向西北方遷移。

    3 ?討論

    3.1 ?模型分析

    本研究中紅茴砂MaxEnt生態(tài)位模型是以最大熵理論為基礎(chǔ),將紅茴砂分布數(shù)據(jù)和全球變化下氣候環(huán)境變量相結(jié)合模擬紅茴砂分布,同時(shí)轉(zhuǎn)移、投射至其他區(qū)域,來預(yù)測(cè)紅茴砂在氣候變化下的潛在分布,為其保育提供支持。已有研究表明MaxEnt生態(tài)位模型雖然具有運(yùn)算穩(wěn)定、運(yùn)算時(shí)間短、結(jié)果精度高等優(yōu)點(diǎn),但MaxEnt模型對(duì)采樣偏差敏感,在模擬物種潛在分布時(shí)易產(chǎn)生過度擬合的問題,使模型轉(zhuǎn)移能力降低[28],為提高紅茴砂MaxEnt生態(tài)位模型轉(zhuǎn)移能力,本研究從分布點(diǎn)篩選和環(huán)境因子選擇這兩個(gè)方面進(jìn)行了控制。在分布點(diǎn)數(shù)據(jù)處理上采用緩沖區(qū)分析法進(jìn)行篩選、校對(duì)。環(huán)境因子選擇時(shí)考慮到紅茴砂為一種林下植物,其光照環(huán)境與植物群落息息相關(guān),所以模型建立未選擇太陽輻射因子[29],最后選取了影響植物生長(zhǎng)和分布的溫度、降雨和海拔等主要非生物環(huán)境因子運(yùn)用Pearson相關(guān)性分析,在相關(guān)性較高的環(huán)境因子中選取對(duì)紅茴砂生物意義更高的作為其關(guān)鍵影響因子進(jìn)行建模,減少因過度擬合產(chǎn)生的誤差。結(jié)果顯示本研究紅茴砂三個(gè)MaxEnt模型檢驗(yàn)AUC值均達(dá)0.98以上,表明模型模擬效果好,可信度高。

    然而,本研究也存在一定的局限性。除已選擇的非生物環(huán)境因子變量外,還有其他一些環(huán)境因素也影響紅茴砂分布,比如土壤因子、地形特征和紅茴砂傳播能力、種間競(jìng)爭(zhēng)力、生物相互作用等生物環(huán)境因子[30],基于以上問題,本模型預(yù)測(cè)的基礎(chǔ)生態(tài)位比紅茴砂實(shí)際占據(jù)的生態(tài)位要寬。所以應(yīng)盡可能收集全面且準(zhǔn)確的分布數(shù)據(jù)、環(huán)境數(shù)據(jù)、以及其他因素?cái)?shù)據(jù),加強(qiáng)紅茴砂生理生態(tài)應(yīng)答與環(huán)境要素的耦合關(guān)系研究,有助于了解未來全球氣候變暖背景下紅茴砂可能發(fā)生的分布變化,對(duì)其物種保護(hù)具有現(xiàn)實(shí)意義。

    3.2 ?氣候因子對(duì)紅茴砂分布的影響

    根據(jù)環(huán)境因子貢獻(xiàn)率的結(jié)果,Bio13(最濕月降水量),Bio18(最暖季度降水量)2個(gè)氣候環(huán)境因子對(duì)紅茴砂分布影響最大,其貢獻(xiàn)率累計(jì)高達(dá)57.3%,最濕月降水量454 mm,最暖季度降水量600 mm時(shí)為其最佳生長(zhǎng)條件,表明紅茴砂是對(duì)降水最為敏感的植物,適合在濕潤(rùn)環(huán)境生存,需要足夠水分滿足其生長(zhǎng)需求。同時(shí)溫度對(duì)紅茴砂也是影響紅茴砂分布的主要限制因素,最冷季度平均溫度在3 ℃以上才能使其正常生長(zhǎng),最熱季度均溫適宜范圍為19.6~31.9 ℃,年均溫差不超過15 ℃,這也與紅茴砂當(dāng)前分布較廣的熱帶和亞熱帶潮濕常綠林生境高度吻合[31]。因此寒潮少,熱量資源豐富,降雨充沛,常年濕潤(rùn)的區(qū)域利于紅茴砂分布,可為紅茴砂資源的保護(hù)和合理的開發(fā)利用提供參考。

    3.3 ?氣候變化背景下紅茴砂的分布區(qū)

    2050年RCP2.6和RCP8.5兩種未來氣候情景變化下紅茴砂在華中和華南等區(qū)域潛在分布區(qū)全部喪失,可能與全球氣候持續(xù)變暖趨勢(shì)下這兩個(gè)地區(qū)降水量的減少有關(guān)。未來21世紀(jì)華南地區(qū)溫度整體處于上升趨勢(shì),但年降水持續(xù)減少[32],雖然夏季強(qiáng)降水頻率增大但冬旱現(xiàn)象加劇[33],而紅茴砂作為對(duì)降水最為敏感的植物,廣東、廣西和海南等華南地區(qū)將不再適宜其生存,所以紅茴砂適生區(qū)在未來氣候下逐漸向西北方轉(zhuǎn)移,以減緩極端高溫和冬旱所帶來的危害。

    2種未來氣候環(huán)境下紅茴砂潛在分布收縮集中在川南、滇中北地區(qū)碎片化分布,并呈現(xiàn)沿西藏延伸分布現(xiàn)象,此地區(qū)處于過渡帶,地貌復(fù)雜,高程落差大[34],形成了氣候獨(dú)特,穩(wěn)定的片段化生境,可抵抗全球變暖帶來的極端氣候等不利影響,同時(shí)怒江、瀾滄江、金沙江、雅礱江、大渡河等主要河流在該區(qū)域也呈南北走向,一直延伸到青藏高原內(nèi)部,為紅茴砂向西藏和川北地區(qū)延伸分布提供條件[35] 。中國西南地區(qū)未來21世紀(jì)氣候大部分地區(qū)呈暖濕型發(fā)展的趨勢(shì)[36],區(qū)域溫度逐漸升高、降水量和相對(duì)濕度也逐步增加,許多學(xué)者也對(duì)此地區(qū)潛在植被分布進(jìn)行了分析,21世紀(jì)未來氣候下與紅茴砂潛在分布區(qū)相吻的植被類型主要為亞熱濕潤(rùn)常綠闊葉林、暖熱潮濕落葉-常綠闊葉林、亞熱潮濕闊葉林類、溫帶常綠針葉灌木和部分熱帶森林草地與[37-38],與當(dāng)前紅茴砂分布點(diǎn)植被類型相似,由此可見在全球氣候持續(xù)變暖趨勢(shì)下,西南地區(qū)滇中北和川南等地區(qū)能逐步滿足紅茴砂分布對(duì)地形、植被、熱量資源豐富、水分需求高的要求,成為紅茴砂新的避難所。此現(xiàn)象與前人關(guān)于氣候變化影響植物分布變化的研究相符,認(rèn)為當(dāng)前物種分布范圍及范圍大小會(huì)隨著全球氣溫的升高逐漸向高緯度和高海拔遷移[39-40]。PCR8.5氣候模型情景相較于PCR2.6氣候情景模型紅茴砂潛在適宜分布區(qū)域遷移更遠(yuǎn),向西南地區(qū)新增適生面積更大,即未來全球升溫閾值越高,紅茴砂適合生分布范圍向西南地區(qū)遷移現(xiàn)象更明顯。

    3.4 ?紅茴砂的保育策略

    本研究結(jié)果對(duì)紅茴砂的保護(hù)和生物地理學(xué)研究都具有重要的意義。由前所述紅茴砂應(yīng)對(duì)氣候變暖的一條逃生路線是向高海拔山地遷移,但是目前紅茴砂在中國的自然分布區(qū)僅為海南島,分布范圍十分狹窄,如果紅茴砂持續(xù)向高緯度和高海拔地區(qū)遷移,在遷移途中不可避免遇到海洋等自然地理障礙阻隔導(dǎo)致局部滅絕[41]。所以需要人為干預(yù)降低其在中國滅絕的風(fēng)險(xiǎn)。可根據(jù)本文對(duì)紅茴砂適生區(qū)預(yù)測(cè)情況結(jié)合其瀕危機(jī)制研究,綜合考慮多方面的因素,制定合理保育策略。對(duì)于紅茴砂當(dāng)前分布區(qū)可采取就地保護(hù)措施,例如可在海南島中南部山地地區(qū)建立紅茴砂保育關(guān)鍵區(qū),通過定點(diǎn)觀察、動(dòng)態(tài)監(jiān)測(cè)結(jié)合人工復(fù)壯,最大確保其遺傳多樣性。而在新增適合生境方面,可以考慮遷地保育方法,如紅茴砂潛在分布中心云南地區(qū),可進(jìn)行人工引種栽培,制定合理的生態(tài)調(diào)控措施,避免紅茴砂遺傳資源喪失風(fēng)險(xiǎn),同時(shí)減少對(duì)當(dāng)?shù)厣鷳B(tài)系統(tǒng)平衡影響。

    參考文獻(xiàn)

    [1] 李利平, 何思源, 蔣樣明, 等. 物種分布區(qū)特征及其對(duì)生物多樣性保育的意義[J]. 中國科學(xué): 生命科學(xué), 2019, 49(8): 929-937.

    [2] 朱耿平, 劉國卿, 卜文俊, 等. 生態(tài)位模型的基本原理及其在生物多樣性保護(hù)中的應(yīng)用[J]. 生物多樣性, 2013, 21(1): 90-98.

    [3] Fitzpatrick M C, Gove A D, Sanders N J, et al. Climate change, plant migration, and range collapse in a global bio-diversity hotspot: the Banksia (Proteaceae) of Western Aus-tralia[J]. Global Change Biology, 2008, 14(6): 1337-1352.

    [4] 吳建國, 呂佳佳, 周巧富. 我國珍稀瀕危物種適應(yīng)氣候變化的對(duì)策探討[J]. 中國人口·資源與環(huán)境, 2011, 21(S1): 566-570.

    [5] 喬慧捷, 胡軍華, 黃繼紅. 生態(tài)位模型的理論基礎(chǔ)、發(fā)展方向與挑戰(zhàn)[J]. 中國科學(xué): 生命科學(xué), 2013, 43(11): 915-927.

    [6] 楊 ?楠, 馬東源, 鐘 ?雪, 等. 基于MaxEnt模型的四川王朗國家級(jí)自然保護(hù)區(qū)藍(lán)馬雞棲息地適宜性評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào), 2020, 40(19): 7064-7072.

    [7] 袁智文, 徐愛春, 俞平新, 等. 浙江清涼峰國家級(jí)自然保護(hù)區(qū)華南梅花鹿棲息地適宜性評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào), 2020, 40(18): 6672-6677.

    [8] 侯沁文, 白海艷, 李云玲, 等. 馬鈴薯甲蟲在中國的適生區(qū)[J]. 生態(tài)學(xué)雜志, 2020, 39(10): 3311-3319.

    [9] 張靜涵, 張軒波, 沈美亞, 等. 互花米草在中國的潛在生存區(qū)域預(yù)測(cè)[J]. 濕地科學(xué)與管理, 2020, 16(3): 42-46.

    [10] 趙金鵬, 王閆利, 陸興利, 等. 軟棗獼猴桃在中國的適生區(qū)分析及對(duì)未來氣候變化的響應(yīng)[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2020, 28(10): 1523-1532.

    [11] 檀逸虹, 張喜娟, 原樹生, 等. 全球氣候變化背景下秦艽生態(tài)適宜性預(yù)測(cè)[J]. 生態(tài)學(xué)雜志, 2020, 39(11): 3766-3773.

    [12] Kruijer H J D, Raes N, Stech M. Modelling the distribution of the moss species Hypopterygium tamarisci (Hypoptery-giaceae, Bryophyta) in Central and South America[J]. Nova Hedwigia, 2010, 91(3-4): 399-420.

    [13] 李 ?響, 張成福, 賀 ?帥, 等. MaxEnt模型綜合應(yīng)用研究進(jìn)展分析[J]. 綠色科技, 2020(14): 14-17.

    [14] Wong K C, Sivasothy Y, Boey P L, et al. Essential oils of Etlingera elatior (Jack) R. M. Smith and Etlingera littoralis (Koenig) giseke[J]. Journal of Essential Oil Research, 2010, 22(5): 461-466.

    [15] Fajemiroye J O, Galdino P M, De Paula J A M, et al. Anxi-olytic and antidepressant like effects of natural food flavour (E)-methyl isoeugenol[J]. Food & Function, 2014, 5(8): 1819-1828.

    [16] Rossi P G, Bao L, Luciani A, et al. (E)-Methy?lisoeugenol and elemicin: antibacterial components of Daucus carota L. essential oil against Campylobacter jejuni[J]. Journal Agri-cultural Food Chemistry, 2007, 55(18): 7332-7336.

    [17] 鐘玉玲. 姜科植物景觀園規(guī)劃設(shè)計(jì)[D]. 廣州: 仲愷農(nóng)業(yè)工程學(xué)院, 2017: 21.

    [18] 覃海寧, 楊 ?永, 董仕勇, 等. 中國高等植物受威脅物種名錄[J]. 生物多樣性, 2017, 25(7): 696-744.

    [19] Norhayat S N, Mohd. Sirat H. Chemical constituents of Etlingera littoralis (konig.) giseke (red flower variety)[J]. Zamm Journal of Applied Mathematics & Mechanics, 2010, 64(3): 200.

    [20] Lutfia A, Munir E, Yurnaliza Y, et al. Screening of antago-nistic fungi from Etlingera littoralis (J. Knig) Giseke rhi-zome in Sibolangit Forest, North Sumatra[J]. Journal of Physics: Conference Series, 2019, 1351: 012025.

    [21] Kittipanangkul N, Ngamriabsakul C. Pollen and Pollinator Limitation of Seed Initiation in Etlingera littoralis (J. K?nig) Giseke (Zingiberaceae) in Klong Klai Basin, Khao Nan Na-tional Park, Thailand[J]. Walailak Journal of Science & Technology, 2006, 3(2): 207-217.

    [22] 歐云鋒, 檀廣兵, 王英強(qiáng). 姜科熱帶分布與高山分布類群傳粉生物學(xué)的比較研究——以茴香砂仁屬和象牙參屬的代表種為例[C]//廣東省植物學(xué)會(huì). 廣東省植物學(xué)會(huì)第十九期學(xué)術(shù)研討會(huì)論文集. 廣州: 廣東省植物學(xué)會(huì), 2010: 27-28.

    [23] 王茹琳, 李 ?慶, 封傳紅, 等. 基于MaxEnt的西藏飛蝗在中國的適生區(qū)預(yù)測(cè)[J]. 生態(tài)學(xué)報(bào), 2017, 37(24): 8556-8566.

    [24] Riahi K, Rao S, Krey V, et al. RCP 8.5—A scenario of comparatively high greenhouse gas emissions[J]. Climatic Change, 2011, 109(1-2): 33.

    [25] 王小軍, 劉光旭, 肖 ?彤. 氣候變化情景下油茶生長(zhǎng)的適宜性特征[J]. 熱帶地理, 2020, 40(5): 868-880.

    [26] 陳衍如, 謝慧敏, 羅火林, 等. 氣候變化對(duì)寒蘭分布的影響及其分布格局模擬[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2019, 30(10): 3419-3425.

    [27] Swets J. Measuring the accuracy of diagnostic systems[J]. Science ,1988, 240(4857): 1285-1293.

    [28] 朱耿平, 喬慧捷. Maxent模型復(fù)雜度對(duì)物種潛在分布區(qū)預(yù)測(cè)的影響[J]. 生物多樣性, 2016, 24(10): 1189-1196.

    [29] 胡 ?秀, 郭 ?微, 吳福川, 等. MaxEnt生態(tài)學(xué)模型在野生植物近自然林引種區(qū)劃中的應(yīng)用——以紅姜花為例[J]. 廣西植物, 2015, 35(3): 325-330.

    [30] 張殷波, 高晨虹, 秦 ?浩. 山西翅果油樹的適生區(qū)預(yù)測(cè)及其對(duì)氣候變化的響應(yīng)[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2018, 29(4): 1156-1162.

    [31] Chongkraijak W, Ngamnriabsakul C, Poulsen A D. Mor-phological diversity and distribution of Etlingera littoralis (K?nig) Giseke (Zingiberaceae) in Southern Thailand[J]. Walailak Journal of Science & Technology, 2013, 10(6): 643-656.

    [32] 黃曉瑩, 溫之平, 杜堯東, 等. 華南地區(qū)未來地面溫度和降水變化的情景分析[J]. 熱帶氣象學(xué)報(bào), 2008(3): 254-258.

    [33] 陳曉晨, 徐 ?影, 姚 ?遙. 不同升溫閾值下中國地區(qū)極端氣候事件變化預(yù)估[J]. 大氣科學(xué), 2015, 39(6): 1123-1135.

    [34] 尤麗鈺, 霍義強(qiáng), 陳隆勛. 云南雨季開始和大氣環(huán)流季節(jié)變化的關(guān)系[J]. 云南大學(xué)學(xué)報(bào)(自然科學(xué)版), 1981(1): 96-112.

    [35] 李寧寧, 張愛平, 張 ?林, 等. 氣候變化下青藏高原兩種云杉植物的潛在適生區(qū)預(yù)測(cè)[J]. 植物研究, 2019, 39(3): 395-406.

    [36] 劉曉冉, 程炳巖, 楊 ?茜, 等. 我國西南地區(qū)21世紀(jì)氣候變化的情景預(yù)估分析[J]. 西南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 34(9): 82-89.

    [37] 任正超, 朱華忠, 史 ?華, 等. 最后間冰期至未來2070s中國潛在自然植被時(shí)空分布格局及其對(duì)氣候變化的響應(yīng)[J]. 自然資源學(xué)報(bào), 2020, 35(6): 1484-1498.

    [38] 陸雙飛, 殷曉潔, 韋晴雯, 等. 氣候變化下西南地區(qū)植物功能型地理分布響應(yīng)[J]. 生態(tài)學(xué)報(bào), 2020, 40(1): 310-324.

    [39] Mccarty J P. Ecological consequences of recent climate change[J]. Conservation Biology, 2001, 15(2): 320-33.

    [40] 賴云云, 馮建孟, 袁媛媛. 氣候變化對(duì)尼泊爾地區(qū)熱帶植物海拔分布格局的影響[J]. 信陽師范學(xué)院學(xué)報(bào)(自然科學(xué)版), 2018, 31(2): 233-239.

    [41] 黎 ?磊, 陳家寬. 氣候變化對(duì)野生植物的影響及保護(hù)對(duì)策[J]. 生物多樣性, 2014, 22(5): 549-563.

    責(zé)任編輯:崔麗虹

    猜你喜歡
    氣候變化預(yù)測(cè)
    北極熊如何適應(yīng)氣候變化
    英語世界(2023年10期)2023-11-17 09:19:12
    無可預(yù)測(cè)
    黃河之聲(2022年10期)2022-09-27 13:59:46
    選修2-2期中考試預(yù)測(cè)卷(A卷)
    選修2-2期中考試預(yù)測(cè)卷(B卷)
    選修2—2期中考試預(yù)測(cè)卷(A卷)
    《應(yīng)對(duì)氣候變化報(bào)告(2022)》發(fā)布
    保護(hù)生物多樣性 積極應(yīng)對(duì)氣候變化
    不可預(yù)測(cè)
    氣候變化與環(huán)保法官
    氣候變化:法官的作用
    国产精品秋霞免费鲁丝片| 亚洲成人国产一区在线观看 | 麻豆乱淫一区二区| 午夜久久久在线观看| 最新在线观看一区二区三区 | 国产成人av激情在线播放| 成人18禁高潮啪啪吃奶动态图| 多毛熟女@视频| 日韩 亚洲 欧美在线| 一级毛片女人18水好多 | 国产无遮挡羞羞视频在线观看| 大片电影免费在线观看免费| 精品国产乱码久久久久久小说| 夫妻性生交免费视频一级片| 国产免费又黄又爽又色| 久久久精品区二区三区| 午夜两性在线视频| 午夜两性在线视频| 99热全是精品| 一级黄片播放器| 国产精品国产三级国产专区5o| 下体分泌物呈黄色| 50天的宝宝边吃奶边哭怎么回事| 香蕉国产在线看| 大话2 男鬼变身卡| 国产亚洲欧美在线一区二区| a级片在线免费高清观看视频| 精品人妻一区二区三区麻豆| 亚洲精品第二区| 久久性视频一级片| 亚洲av日韩精品久久久久久密 | 久久av网站| 国产又爽黄色视频| 久久久久久人人人人人| 一级a爱视频在线免费观看| 99国产精品99久久久久| 国产片内射在线| av片东京热男人的天堂| 两性夫妻黄色片| 国产成人精品久久二区二区91| 国产在线免费精品| 久久鲁丝午夜福利片| 啦啦啦 在线观看视频| 黄色视频不卡| www.av在线官网国产| 亚洲精品久久成人aⅴ小说| 久热这里只有精品99| 巨乳人妻的诱惑在线观看| 丝袜喷水一区| 黄网站色视频无遮挡免费观看| 中文字幕精品免费在线观看视频| 午夜免费鲁丝| 欧美国产精品一级二级三级| 青春草亚洲视频在线观看| 在线观看国产h片| www.精华液| 美女扒开内裤让男人捅视频| 亚洲精品国产一区二区精华液| bbb黄色大片| 777米奇影视久久| 国产亚洲午夜精品一区二区久久| 日本猛色少妇xxxxx猛交久久| 亚洲成国产人片在线观看| 久久久久国产一级毛片高清牌| 首页视频小说图片口味搜索 | 国产成人a∨麻豆精品| 亚洲情色 制服丝袜| 久久热在线av| av一本久久久久| 黄色a级毛片大全视频| av天堂久久9| 多毛熟女@视频| 免费久久久久久久精品成人欧美视频| 国产男人的电影天堂91| 精品熟女少妇八av免费久了| 欧美日韩成人在线一区二区| 免费看十八禁软件| 亚洲国产精品一区二区三区在线| 丝袜美腿诱惑在线| 成年人免费黄色播放视频| 欧美黑人欧美精品刺激| 欧美黑人精品巨大| 日韩电影二区| 男女无遮挡免费网站观看| 亚洲第一青青草原| 亚洲国产精品国产精品| 亚洲成人免费电影在线观看 | 19禁男女啪啪无遮挡网站| 在线观看免费视频网站a站| 国产一区有黄有色的免费视频| 精品亚洲乱码少妇综合久久| 热re99久久精品国产66热6| 黄色视频在线播放观看不卡| 中文字幕精品免费在线观看视频| 午夜日韩欧美国产| 啦啦啦在线免费观看视频4| 50天的宝宝边吃奶边哭怎么回事| 青春草亚洲视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 精品久久蜜臀av无| av又黄又爽大尺度在线免费看| 免费久久久久久久精品成人欧美视频| 久久久久视频综合| 欧美中文综合在线视频| 亚洲 欧美一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩精品网址| 亚洲国产欧美日韩在线播放| 一边摸一边做爽爽视频免费| 在现免费观看毛片| 精品国产国语对白av| 我要看黄色一级片免费的| 日本欧美国产在线视频| 久热爱精品视频在线9| 亚洲国产成人一精品久久久| 亚洲成人免费av在线播放| 在线观看免费视频网站a站| 免费高清在线观看日韩| 美女主播在线视频| 男人操女人黄网站| 日本vs欧美在线观看视频| 大片电影免费在线观看免费| 国产亚洲欧美精品永久| 永久免费av网站大全| 久久久久久久久久久久大奶| 成年人黄色毛片网站| 精品一区二区三区四区五区乱码 | 国产99久久九九免费精品| 天天躁日日躁夜夜躁夜夜| 精品福利观看| 亚洲伊人色综图| 美女午夜性视频免费| 99热全是精品| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产最新在线播放| 亚洲精品国产一区二区精华液| 搡老乐熟女国产| 亚洲av男天堂| 亚洲五月婷婷丁香| 99香蕉大伊视频| 亚洲欧美一区二区三区黑人| 免费人妻精品一区二区三区视频| 色94色欧美一区二区| 一二三四社区在线视频社区8| 久久亚洲精品不卡| 成人国产一区最新在线观看 | 人人妻人人澡人人爽人人夜夜| 国产精品久久久久久人妻精品电影 | 捣出白浆h1v1| 亚洲精品久久成人aⅴ小说| 久久天堂一区二区三区四区| 国产深夜福利视频在线观看| 丁香六月天网| 老汉色∧v一级毛片| 国产在视频线精品| 亚洲国产精品成人久久小说| 国产成人欧美| 少妇被粗大的猛进出69影院| 老司机影院毛片| 久久精品亚洲av国产电影网| 可以免费在线观看a视频的电影网站| 国产精品 欧美亚洲| 91精品伊人久久大香线蕉| 曰老女人黄片| 国产深夜福利视频在线观看| 国产午夜精品一二区理论片| 老汉色av国产亚洲站长工具| 51午夜福利影视在线观看| 99香蕉大伊视频| 国产成人欧美在线观看 | 午夜免费鲁丝| 秋霞在线观看毛片| 考比视频在线观看| 狠狠婷婷综合久久久久久88av| 国产精品偷伦视频观看了| www.精华液| 精品少妇内射三级| 欧美 日韩 精品 国产| 免费在线观看完整版高清| 日韩av在线免费看完整版不卡| 男女边吃奶边做爰视频| netflix在线观看网站| 国产高清视频在线播放一区 | 国产成人一区二区在线| 欧美精品人与动牲交sv欧美| 精品久久蜜臀av无| 天天躁日日躁夜夜躁夜夜| 久久久久久人人人人人| 亚洲精品美女久久久久99蜜臀 | 午夜福利视频精品| 我要看黄色一级片免费的| 亚洲欧美一区二区三区国产| 免费看av在线观看网站| 亚洲国产成人一精品久久久| av天堂在线播放| 国产精品麻豆人妻色哟哟久久| tube8黄色片| 久久狼人影院| 亚洲国产精品999| 90打野战视频偷拍视频| 99国产综合亚洲精品| 亚洲av欧美aⅴ国产| 久久九九热精品免费| 国产欧美日韩综合在线一区二区| 婷婷丁香在线五月| 91老司机精品| 人妻人人澡人人爽人人| 欧美日韩亚洲高清精品| 免费少妇av软件| 国产高清不卡午夜福利| 亚洲五月婷婷丁香| 男人操女人黄网站| 国产熟女午夜一区二区三区| 日本vs欧美在线观看视频| 亚洲精品国产一区二区精华液| 国产成人av激情在线播放| 王馨瑶露胸无遮挡在线观看| 美女国产高潮福利片在线看| 999精品在线视频| 日本五十路高清| 欧美日本中文国产一区发布| 欧美xxⅹ黑人| 中文欧美无线码| 少妇的丰满在线观看| 国产精品成人在线| 久久九九热精品免费| 国产成人免费观看mmmm| 一级毛片电影观看| 国产成人精品久久二区二区免费| 欧美精品啪啪一区二区三区 | 中国国产av一级| 欧美乱码精品一区二区三区| 久热爱精品视频在线9| 好男人电影高清在线观看| 国产在线观看jvid| 久久av网站| 亚洲九九香蕉| 亚洲精品日韩在线中文字幕| 欧美人与性动交α欧美精品济南到| 日韩伦理黄色片| 18禁裸乳无遮挡动漫免费视频| 搡老乐熟女国产| 无限看片的www在线观看| 蜜桃国产av成人99| 男女无遮挡免费网站观看| 大香蕉久久网| 婷婷色麻豆天堂久久| 精品久久久久久电影网| 亚洲情色 制服丝袜| 中文字幕av电影在线播放| 在现免费观看毛片| 悠悠久久av| www.av在线官网国产| 亚洲熟女精品中文字幕| 日本av手机在线免费观看| 一级黄色大片毛片| av线在线观看网站| 一本一本久久a久久精品综合妖精| 美女高潮到喷水免费观看| 最新的欧美精品一区二区| 后天国语完整版免费观看| 久久久精品国产亚洲av高清涩受| 国产亚洲午夜精品一区二区久久| 久久99热这里只频精品6学生| 亚洲精品国产av成人精品| 三上悠亚av全集在线观看| 国产又色又爽无遮挡免| 亚洲国产精品一区二区三区在线| 爱豆传媒免费全集在线观看| 91精品伊人久久大香线蕉| 我要看黄色一级片免费的| 每晚都被弄得嗷嗷叫到高潮| 午夜免费观看性视频| 人人妻人人澡人人爽人人夜夜| 麻豆av在线久日| 精品国产一区二区久久| 国产人伦9x9x在线观看| 黄色毛片三级朝国网站| 老司机深夜福利视频在线观看 | 久久人人爽人人片av| 在线观看一区二区三区激情| 久久国产精品大桥未久av| 午夜av观看不卡| av天堂久久9| 99热全是精品| 精品国产超薄肉色丝袜足j| 下体分泌物呈黄色| 黄色怎么调成土黄色| 丰满人妻熟妇乱又伦精品不卡| 男女国产视频网站| 国精品久久久久久国模美| 欧美精品啪啪一区二区三区 | 亚洲av电影在线进入| av福利片在线| 人妻人人澡人人爽人人| 欧美日韩亚洲高清精品| 成人午夜精彩视频在线观看| 美女高潮到喷水免费观看| 成人手机av| 校园人妻丝袜中文字幕| 一边摸一边做爽爽视频免费| 丝袜喷水一区| 男女之事视频高清在线观看 | 亚洲欧美日韩高清在线视频 | 国产男女内射视频| 亚洲专区中文字幕在线| 母亲3免费完整高清在线观看| 久久影院123| 亚洲av综合色区一区| 亚洲,欧美,日韩| a级毛片在线看网站| 久久久精品94久久精品| 亚洲熟女毛片儿| 捣出白浆h1v1| 观看av在线不卡| 久久国产亚洲av麻豆专区| 国产爽快片一区二区三区| 久久人妻福利社区极品人妻图片 | 国产欧美日韩一区二区三区在线| 色网站视频免费| 亚洲av男天堂| 看免费av毛片| 国产一级毛片在线| 欧美激情极品国产一区二区三区| 日本一区二区免费在线视频| 国产成人精品久久二区二区免费| 人人妻,人人澡人人爽秒播 | 中文欧美无线码| 黄片播放在线免费| 老司机靠b影院| 欧美激情极品国产一区二区三区| 视频区图区小说| 好男人视频免费观看在线| av线在线观看网站| 色94色欧美一区二区| 精品久久久久久电影网| 人人澡人人妻人| 国产亚洲一区二区精品| 两个人免费观看高清视频| 男的添女的下面高潮视频| 一本—道久久a久久精品蜜桃钙片| 国语对白做爰xxxⅹ性视频网站| 欧美日韩亚洲高清精品| 成年人免费黄色播放视频| 亚洲色图 男人天堂 中文字幕| 黄片播放在线免费| 丰满饥渴人妻一区二区三| 国产精品国产av在线观看| 亚洲一区二区三区欧美精品| 多毛熟女@视频| 国产黄色视频一区二区在线观看| 久久久精品免费免费高清| 国产黄色免费在线视频| 久久午夜综合久久蜜桃| 久久久精品区二区三区| 国产伦理片在线播放av一区| 国产黄频视频在线观看| 国产一区二区 视频在线| 欧美成人午夜精品| 国产精品久久久人人做人人爽| 精品免费久久久久久久清纯 | 婷婷色综合大香蕉| 久久精品久久精品一区二区三区| 久久久久久久大尺度免费视频| 一二三四社区在线视频社区8| 香蕉国产在线看| 高潮久久久久久久久久久不卡| 精品国产超薄肉色丝袜足j| 欧美激情极品国产一区二区三区| 国产精品秋霞免费鲁丝片| 好男人视频免费观看在线| 极品少妇高潮喷水抽搐| 日本a在线网址| 婷婷丁香在线五月| 美女国产高潮福利片在线看| 欧美在线黄色| 亚洲av欧美aⅴ国产| 国产一区亚洲一区在线观看| 高清欧美精品videossex| 亚洲精品久久成人aⅴ小说| 大型av网站在线播放| 精品一区在线观看国产| 亚洲五月色婷婷综合| 99re6热这里在线精品视频| 国产熟女午夜一区二区三区| 啦啦啦在线免费观看视频4| 亚洲 欧美一区二区三区| 人人澡人人妻人| 国产欧美日韩一区二区三区在线| 99国产综合亚洲精品| 高清av免费在线| 国产视频一区二区在线看| 亚洲成av片中文字幕在线观看| 18在线观看网站| 国产日韩一区二区三区精品不卡| 天天躁夜夜躁狠狠躁躁| 亚洲国产精品一区二区三区在线| svipshipincom国产片| 婷婷成人精品国产| 精品免费久久久久久久清纯 | 久久影院123| 超碰成人久久| 国产淫语在线视频| 欧美精品人与动牲交sv欧美| 久久久久久久久免费视频了| 亚洲图色成人| 母亲3免费完整高清在线观看| 人人澡人人妻人| 婷婷色av中文字幕| 日韩熟女老妇一区二区性免费视频| e午夜精品久久久久久久| 一边摸一边抽搐一进一出视频| h视频一区二区三区| 亚洲精品自拍成人| 人人妻人人添人人爽欧美一区卜| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲综合一区二区三区_| 亚洲成国产人片在线观看| 国产日韩欧美视频二区| 啦啦啦 在线观看视频| 精品久久久久久电影网| 久久精品亚洲熟妇少妇任你| 熟女av电影| 菩萨蛮人人尽说江南好唐韦庄| 中文乱码字字幕精品一区二区三区| 亚洲第一青青草原| 黄色片一级片一级黄色片| 热99国产精品久久久久久7| 久久毛片免费看一区二区三区| 青草久久国产| 亚洲国产最新在线播放| 两个人免费观看高清视频| 只有这里有精品99| 美女扒开内裤让男人捅视频| 久久午夜综合久久蜜桃| 欧美另类一区| 91字幕亚洲| 最近最新中文字幕大全免费视频 | 欧美在线黄色| 日韩av不卡免费在线播放| 啦啦啦在线免费观看视频4| 欧美精品啪啪一区二区三区 | 国产片特级美女逼逼视频| 精品人妻在线不人妻| 1024视频免费在线观看| 午夜日韩欧美国产| 80岁老熟妇乱子伦牲交| 午夜老司机福利片| 午夜两性在线视频| 免费观看人在逋| 亚洲精品av麻豆狂野| 精品人妻在线不人妻| 精品久久久久久电影网| 亚洲熟女毛片儿| 欧美激情高清一区二区三区| 亚洲欧美日韩高清在线视频 | 伦理电影免费视频| 精品亚洲成国产av| 一区二区三区精品91| 欧美国产精品va在线观看不卡| 岛国毛片在线播放| 亚洲欧洲精品一区二区精品久久久| videos熟女内射| 日本一区二区免费在线视频| 国产精品久久久久久精品电影小说| 一级毛片黄色毛片免费观看视频| 国产免费视频播放在线视频| 免费看不卡的av| 性少妇av在线| 热re99久久精品国产66热6| 国产成人精品无人区| 在线观看免费高清a一片| 80岁老熟妇乱子伦牲交| 七月丁香在线播放| 晚上一个人看的免费电影| 久久久亚洲精品成人影院| 国产在视频线精品| 高清视频免费观看一区二区| 婷婷丁香在线五月| av在线app专区| 色婷婷av一区二区三区视频| 97在线人人人人妻| 在线av久久热| 久久久国产欧美日韩av| 亚洲av男天堂| 18禁观看日本| 少妇人妻 视频| www.自偷自拍.com| 一区在线观看完整版| 热re99久久精品国产66热6| 国产色视频综合| 建设人人有责人人尽责人人享有的| 精品福利永久在线观看| 91九色精品人成在线观看| 精品久久久久久久毛片微露脸 | 亚洲欧美成人综合另类久久久| 啦啦啦啦在线视频资源| 又大又爽又粗| 亚洲专区国产一区二区| 美女脱内裤让男人舔精品视频| 黄色毛片三级朝国网站| www.999成人在线观看| 亚洲欧洲日产国产| 精品亚洲成a人片在线观看| 国产极品粉嫩免费观看在线| 妹子高潮喷水视频| 麻豆av在线久日| 国产精品99久久99久久久不卡| 久久女婷五月综合色啪小说| 国产成人av激情在线播放| 午夜福利乱码中文字幕| 日韩大片免费观看网站| 在线亚洲精品国产二区图片欧美| 只有这里有精品99| 麻豆乱淫一区二区| 这个男人来自地球电影免费观看| 欧美日韩精品网址| 免费女性裸体啪啪无遮挡网站| 亚洲精品国产av成人精品| 国产精品av久久久久免费| 王馨瑶露胸无遮挡在线观看| 99re6热这里在线精品视频| 国产高清videossex| 午夜福利乱码中文字幕| 国产成人欧美在线观看 | 精品人妻在线不人妻| 国产一区亚洲一区在线观看| 啦啦啦啦在线视频资源| 欧美成狂野欧美在线观看| 在线观看www视频免费| 男人操女人黄网站| 女人被躁到高潮嗷嗷叫费观| 99国产精品免费福利视频| 亚洲综合色网址| 亚洲国产看品久久| 自拍欧美九色日韩亚洲蝌蚪91| 国产野战对白在线观看| 色婷婷av一区二区三区视频| 男女午夜视频在线观看| 亚洲国产av新网站| 亚洲图色成人| 少妇人妻 视频| 亚洲伊人色综图| 国产精品秋霞免费鲁丝片| 黑人巨大精品欧美一区二区蜜桃| 亚洲人成电影免费在线| 18禁观看日本| 日本一区二区免费在线视频| 韩国高清视频一区二区三区| 两性夫妻黄色片| 国产一级毛片在线| 十八禁人妻一区二区| 亚洲精品一卡2卡三卡4卡5卡 | 国产av一区二区精品久久| √禁漫天堂资源中文www| 国产免费又黄又爽又色| 一边亲一边摸免费视频| 亚洲国产看品久久| 日韩av免费高清视频| 午夜福利乱码中文字幕| 成人影院久久| 精品少妇黑人巨大在线播放| 国产精品一区二区免费欧美 | 欧美日本中文国产一区发布| xxxhd国产人妻xxx| 国产成人精品在线电影| 王馨瑶露胸无遮挡在线观看| 久久人妻福利社区极品人妻图片 | 日韩熟女老妇一区二区性免费视频| av在线播放精品| 久久精品国产亚洲av涩爱| 亚洲 欧美一区二区三区| 性高湖久久久久久久久免费观看| 精品国产国语对白av| 国产熟女欧美一区二区| 欧美 亚洲 国产 日韩一| 日韩人妻精品一区2区三区| 建设人人有责人人尽责人人享有的| 观看av在线不卡| 日韩制服丝袜自拍偷拍| 色婷婷av一区二区三区视频| 女人高潮潮喷娇喘18禁视频| 婷婷丁香在线五月| 亚洲视频免费观看视频| 50天的宝宝边吃奶边哭怎么回事| 日韩,欧美,国产一区二区三区| 亚洲中文日韩欧美视频| 夫妻性生交免费视频一级片| 王馨瑶露胸无遮挡在线观看| 一级毛片黄色毛片免费观看视频| 9热在线视频观看99| 中文欧美无线码| 精品亚洲乱码少妇综合久久| 亚洲综合色网址| 黄片小视频在线播放| 免费一级毛片在线播放高清视频 | 国产极品粉嫩免费观看在线| tube8黄色片| 免费av中文字幕在线| 国产伦理片在线播放av一区| 国产av国产精品国产| 叶爱在线成人免费视频播放| 国产片特级美女逼逼视频| 嫁个100分男人电影在线观看 | 午夜av观看不卡| 久久久久久久国产电影| 成人国语在线视频| 久久精品aⅴ一区二区三区四区| 9色porny在线观看| 90打野战视频偷拍视频| 欧美精品啪啪一区二区三区 | 午夜免费男女啪啪视频观看| 久久热在线av| 大香蕉久久网| 一区二区三区乱码不卡18|