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

    基于MaxEnt生態(tài)軟件劃分澳洲堅果的潛在地理適生區(qū)

    2017-06-19 18:29:34田榮榮那立妍許新橋
    林業(yè)科學(xué)研究 2017年3期
    關(guān)鍵詞:適生區(qū)氣候因子堅果

    王 偉,田榮榮,那立妍,,楊 穎,許新橋*

    (1.林木遺傳育種國家重點實驗室,國家林業(yè)局林木培育重點實驗室,中國林業(yè)科學(xué)研究院林業(yè)研究所,北京 100091;2.北京城市學(xué)院生物醫(yī)藥學(xué)部,北京 100094; 3.山東省食品藥品檢驗研究院,山東 濟南 250101)

    基于MaxEnt生態(tài)軟件劃分澳洲堅果的潛在地理適生區(qū)

    王 偉1,田榮榮2,那立妍1,2,楊 穎3,許新橋1*

    (1.林木遺傳育種國家重點實驗室,國家林業(yè)局林木培育重點實驗室,中國林業(yè)科學(xué)研究院林業(yè)研究所,北京 100091;2.北京城市學(xué)院生物醫(yī)藥學(xué)部,北京 100094; 3.山東省食品藥品檢驗研究院,山東 濟南 250101)

    [目的]對澳洲堅果在中國的潛在分布區(qū)進行預(yù)測,并對其適生區(qū)進行分析和劃分。[方法] 通過收集澳洲堅果的地理信息數(shù)據(jù),利用最大熵模型(MaxEnt)與地理信息系統(tǒng)(ArcGIS),綜合相關(guān)19項氣候因子,預(yù)測劃定澳洲堅果在世界以及我國的潛在地理分布區(qū)。[結(jié)果]該物種生長區(qū)域狹窄,對環(huán)境要求苛刻。世界范圍內(nèi),澳洲堅果的較適宜生長區(qū)在澳洲東部、南美洲東南部和馬達(dá)加斯加島東部以及亞洲地區(qū)23°26′~30°N,73°~122°E范圍內(nèi)。在我國,澳洲堅果適宜分布區(qū)主要集中在西藏、臺灣、廣西、廣東和云南等地,其高適宜區(qū)面積依次為西藏(15 359 km2),臺灣(14 054 km2),廣西(7 372 km2),廣東(6 147 km2)和云南(3 776 km2)。刀切法(Jackknife)分析顯示, 澳洲堅果分布主要受到極端最高溫、年均氣溫變化范圍、最干月降雨量、溫度季節(jié)性變化和等溫性等氣象因子的影響。[結(jié)論] 本研究用MaxEnt模擬澳洲堅果的潛在地理分布有一定的準(zhǔn)確性,劃分出了澳洲堅果基本的地理分布格局和潛在分布區(qū)域,并闡明了主導(dǎo)其地理分布的生物氣候因子,為澳洲堅果尤其是在我國的引種和推廣應(yīng)用提供了參考。

    澳洲堅果;潛在適生區(qū);MaxEnt;ArcGIS

    澳洲堅果(Macadamiaintegrifolia)又稱夏威夷果、澳洲胡桃、昆士蘭栗、昆士蘭果等,澳洲堅果作為常綠喬木,雙子葉植物,隸屬山龍眼科、澳洲堅果屬。在眾多的干果之中,澳洲堅果具有較高的經(jīng)濟價值[1],有“干果之王”的美譽。澳洲堅果其原產(chǎn)于澳洲,適宜生長在溫和、濕潤、風(fēng)力小的地區(qū)[2]。目前主要分布于澳大利亞東部、新喀里多尼亞、印度尼西亞蘇拉威西島。

    在我國,澳洲堅果作為外來樹種,引種栽培不僅可以豐富我國的堅果市場,而且可以增加當(dāng)?shù)氐慕?jīng)濟效益。我國在1910年即開始引種實生苗,而直至1979年才開始引入澳洲堅果無性系,當(dāng)前在云南、廣東、廣西、海南以及福建等地發(fā)展迅速,尤其在云南栽培面積高達(dá)6.7 hm-2[2]。在我國,盡管開展了一些初步的引種試驗,引種結(jié)果也為我國澳洲堅果新興產(chǎn)業(yè)的發(fā)展范圍和空間提供了一定的參考,但總體上缺乏對其適生區(qū)的劃分,引種和推廣地區(qū)的選擇較盲目,導(dǎo)致部分人力、物力和財力的浪費。明確劃定其潛在的適生區(qū)域,對謀劃我國澳洲堅果產(chǎn)業(yè)的發(fā)展極其重要。本文利用MaxEnt生態(tài)學(xué)模型根據(jù)原產(chǎn)地的氣象條件劃分我國的潛在適生區(qū),旨在為澳洲堅果在我國的引種和推廣提供基本的數(shù)據(jù)支持。

    MaxEnt生態(tài)學(xué)模型最大熵理論[3]是眾多生物地理分布預(yù)測軟件的一種,其是基于有限的已知信息進行無偏推斷預(yù)測未知分布的方法。與同類預(yù)測模型如GARP[4]、ENFA和BIOCLIM等相比,利用AUC分析顯示MaxEnt預(yù)測結(jié)果更優(yōu)[5-7],AUC值越接近1說明預(yù)測準(zhǔn)確度越高,通過比較發(fā)現(xiàn)MaxEnt的AUC值最高。尤其對物種潛在分布區(qū)劃分的預(yù)測,MaxEnt往往能得到較滿意的結(jié)果[8-9]。目前,用該模型已成功預(yù)測了紫花苜蓿[10]、南丹參[11]、金錢松[12]、普洱姜花[13]、加拿大一枝黃花[14]、檀香[15]等在我國的潛在生境分布。

    1 收集數(shù)據(jù)

    1.1 澳洲堅果數(shù)據(jù)的收集

    本文通過全球生物多樣性信息平臺(GBIF)(http://www.gbif.org/)共獲得澳洲堅果地理信息數(shù)據(jù)693個(中國不是原始分布區(qū),故沒有中國的地理信息),每個數(shù)據(jù)都含有相應(yīng)的經(jīng)緯度,部分含有海拔信息。這些地理信息數(shù)據(jù)主要分布在澳大利亞東部,美國有零星分布。將這些地理信息數(shù)據(jù)轉(zhuǎn)換為統(tǒng)一的的經(jīng)緯度坐標(biāo),所有經(jīng)緯度數(shù)據(jù)默認(rèn)東經(jīng)與北緯為正,西經(jīng)與南緯為負(fù)。根據(jù)MaxEnt軟件的需要,將澳洲堅果已知分布點通過EXCEL轉(zhuǎn)化為軟件支持的.csv格式文件(原格式.xls),表格中包含物種名稱、經(jīng)度、緯度等3列內(nèi)容。

    1.2 氣候因子數(shù)據(jù)收集以及相關(guān)地理信息數(shù)據(jù)

    本文使用的氣候因子數(shù)據(jù)來源于世界氣候數(shù)據(jù)庫(http://www.worldclim.org/),選擇下載,精度2.5弧分(2.5 km),得到研究所需的19個生物氣候因子。分別為年均氣溫(bio1)、平均周溫度變化范圍(bio2)、年均氣溫變化范圍(bio3)、極端最高氣溫(bio4)、極端最低氣溫(bio5)、等溫性(bio6)、溫度季節(jié)性變化(bio7)、最濕季平均氣溫(bio8)、最干季平均氣溫(bio9)、最熱季平均氣溫(bio10)、最冷季平均氣溫(bio11)、年降水量(bio12)、最濕月降水量(bio13)、最干月降水量(bio14)、降水量的季節(jié)性變化(bio15)、最干季降水量(bio16)、最濕季降水量(bio17)、最熱季降水量(bio18)、最冷季降水量(bio19)。在DIVA-GIS官方網(wǎng)站(http://www.diva-gis.org/Data)下載全球矢量地圖,在國家地理信息系統(tǒng)網(wǎng)站(http://nfgis.nsdi.gov.cn/)下載1∶400萬中國地圖及中國行政區(qū)劃圖用作研究。

    2 建立模型

    通過ArcGIS(版本10.2)軟件[4]和MaxEnt(版本3.3.3k)[3]最大熵模型對澳洲堅果進行預(yù)測。前期數(shù)據(jù)均由ArcGIS處理,MaxEnt僅用來預(yù)測數(shù)據(jù)。具體見圖1。

    圖1 ArcGIS軟件和MaxEnt軟件預(yù)測流程Fig.1 Based on ArcGIS, using MaxEnt to predict the geographical distribution of Technology Roadmap

    2.1 采用ArcGIS軟件的預(yù)測

    將.csv格式文件與全球矢量地圖導(dǎo)入軟件,通過ArcGIS中Raster to ASCII功能可將柵格研究所需的生物氣候因子轉(zhuǎn)化為MaxEnt所需的.asc格式文件。將.csv格式數(shù)據(jù)導(dǎo)入ArcGIS并導(dǎo)入世界地圖,同樣通過Raster to ASCII功能將其轉(zhuǎn)化為MaxEnt最大熵模型需要的.asc格式文件。

    在數(shù)據(jù)預(yù)測完成后通過ArcGIS中ASCII to Raster功能導(dǎo)入適生區(qū)預(yù)測地圖(.asc),利用Raster Reclassify功能可計算出澳洲堅果在各個國家的高適宜區(qū)面積、較適宜生長區(qū)面積、邊緣適生區(qū)面積和低適生區(qū)面積。

    2.2 MaxEnt最大熵模型的預(yù)測

    MaxEnt最大熵模型,原理即對一個隨機事件的概率分布進行預(yù)測時,預(yù)測應(yīng)當(dāng)滿足全部已知的條件,而對未知的情況不要做任何主觀假設(shè)。在這種情況下,概率分布最均勻,預(yù)測的風(fēng)險最小。因為這時概率分布的信息熵最大,所以人們稱這種模型叫“最大熵模型”。

    MaxEnt軟件分析,將.asc格式的澳洲堅果已知分布點添加進Samples中,在Environmental layers里添加19個生物氣候因子(1950—2000年間的氣候數(shù)據(jù))。設(shè)置75%的分布數(shù)據(jù)作為訓(xùn)練數(shù)據(jù),25%的分布數(shù)據(jù)作為檢驗數(shù)據(jù),將ArcGIS中導(dǎo)出的.csv文件添加到設(shè)置欄下的Bias file中,其他參數(shù)為軟件默認(rèn)值。

    3 結(jié)果與分析

    3.1 澳洲堅果在世界的潛在分布與分析

    利用MaxEnt軟件分析得到澳洲堅果在世界的適生分布區(qū)(表1)。表1顯示:澳洲堅果生長區(qū)域狹窄,對環(huán)境要求苛刻,主要分布在熱帶和亞熱帶的部分地區(qū),主要集中在澳洲堅果的原產(chǎn)地澳洲東部以及南美洲東南部和馬達(dá)加斯加島東部,亞洲23°26′~30°N,73°~122°E是其大面積的適生區(qū),因此,澳洲堅果在這個區(qū)間應(yīng)該有較大的發(fā)展空間。澳洲堅果在亞洲高適宜區(qū)屬于熱帶、亞熱帶地區(qū),這些地區(qū)降水量豐富,空氣濕度大,有利于澳洲堅果的生長。尤其在中國的南部地區(qū)、印度的北部地區(qū)、尼泊爾自西向東的中部地區(qū)、老撾的中部及北部地區(qū)以及緬甸的北部及東北部地區(qū)。澳洲堅果在印度、緬甸、中國、越南、尼泊爾等地均適宜引種種植。

    表1 世界23個國家不同等級適生區(qū)面積

    澳洲堅果在非洲地區(qū)適生范圍較小,主要在馬達(dá)加斯加島的東部、津巴布韋的東部、南非的南部和東部,少量分布在東北部地區(qū)。在莫桑比克,僅在中東部適宜栽培澳洲堅果且面積較小。這些地區(qū)均處在熱帶地區(qū),地形多為平坦的高原或平原,且降水量充沛,年降水量最少可達(dá)1 000 mm。

    在南、北美洲地區(qū),僅厄瓜多爾和墨西哥有最適宜澳洲堅果種植的高適宜區(qū),并且墨西哥的高適宜區(qū)少,大部分屬于低適生區(qū)。澳洲堅果適生區(qū)大多分布在墨西哥的東部沿海地區(qū)、厄瓜多爾的東北部地區(qū)、巴西的南部地區(qū),少量分布于美國的佛羅里達(dá)州、伯利茲的南部地區(qū)、危地馬拉的中部及東部地區(qū)等地;而在大洋洲,除澳洲堅果的原產(chǎn)地澳大利亞外,其適生區(qū)只有新喀里多尼亞。新喀里多尼亞屬熱帶氣候,年平均降水量可達(dá)1 500 mm,整個國家大部分地區(qū)適宜澳洲堅果的栽培種植。

    3.2 澳洲堅果在我國潛在分布區(qū)與分析

    以澳洲堅果的世界適生區(qū)分布圖(.asc格式)為底圖,利用ArcGIS的掩膜提取功能提取中國分布圖以及中國境內(nèi)各省(市、自治州)的分布圖。通過ArcGIS中Reclassify功能對各國以及中國境內(nèi)各省(市、自治州)的適生區(qū)統(tǒng)計適生面積(表2)。通過ArcGIS計算,澳洲堅果在中國的高適宜區(qū)面積為47 216 km2,較適宜生長區(qū)面積為90 391 km2(表1),高適宜區(qū)和較適宜生長區(qū)的累計面積是全球最大的,說明澳洲堅果在我國的發(fā)展空間很大。

    表2顯示:澳洲堅果適宜區(qū)主要集中在我國南部地區(qū),尤其是南部沿海地區(qū),北部地區(qū)沒有澳洲堅果的適生區(qū)。在我國西藏、臺灣、廣西、廣東和云南,澳洲堅果高適宜區(qū)的面積最大。西藏地區(qū)的高適宜區(qū)面積最大,高達(dá)15 359 km2;就高適宜區(qū)占各省份或地區(qū)面積的比例看,中國臺灣最適宜,高適宜生長區(qū)面積為14 054 km2,占中國臺灣總面積的38%,高適宜區(qū)基本上遍布整個地區(qū),而總適生區(qū)占其面積的51%。較適宜生長區(qū)面積較大的主要是廣西(24 355 km2)、廣東(22 398 km2)和云南(20 511 km2)。

    表2 我國8個行政區(qū)不同等級適生區(qū)面積

    由圖2所示:在我國,澳洲堅果較適宜生長區(qū)主要分布在廣東東部、廣西西北和中南部、海南東北部、福建中部和南部、云南中南部以及貴州東南和西南部地區(qū),而在廣東與廣西的澳洲堅果總適生區(qū)面積則分別是其省總面積的69%和71%。

    表1表明:在我國南部地區(qū),澳洲堅果的高適宜區(qū)總面積為47 216 km2,較適宜生長區(qū)總面積為90 391 km2,邊緣適生區(qū)總面積為134 757 km2,低適生區(qū)總面積為444 335 km2,不同適生區(qū)總面積為716 699 km2。該表只列出了預(yù)測分布總面積,分布區(qū)面積內(nèi)的實際可利用土地,會遠(yuǎn)低于預(yù)測面積。

    3.3 預(yù)測結(jié)果評價

    經(jīng)受試者工作曲線檢驗,MaxEnt模型的ROC曲線下面積訓(xùn)練集和測試集的AUC值依次為0.994和0.993,接近于l,表明MaxEnt模型對澳洲堅果潛在分布區(qū)的預(yù)測效果非常好。

    經(jīng)MaxEnt最大熵模型ROC曲線檢驗,澳洲堅果ROC曲線的AUC值為0.993(圖3)。AUC的范圍值為0~1。一般認(rèn)為,當(dāng)AUC≤0.7時,表明預(yù)測價值較低;當(dāng)0.8≤AUC<0.9時,表明預(yù)測價值較好;當(dāng)0.9≤AUC<1時,則表明預(yù)測價值較高。因此,在生態(tài)位模型對物種進行適生區(qū)預(yù)測時,AUC值越接近1代表預(yù)測結(jié)果價值越高,生態(tài)位模型預(yù)測越準(zhǔn)確,其預(yù)測效果越好。

    3.4 影響澳洲堅果分布的不同氣候因子的評價

    用刀切法對環(huán)境氣候因子對分布所作貢獻(xiàn)進行檢測,由圖4可知:影響澳洲堅果分布情況的環(huán)境氣候因子由高到低依次為bio4(極端最高溫)、bio3(年均溫度變化范圍)、bio14(最干月降水量)、bio7(溫度季節(jié)性變化)、bio6(等溫性)、bio18(最熱季降水量)、bio11(最冷季平均溫度)、bio17(最濕季降水量)、bio9(最干季平均氣溫)、bio1(年均氣溫)、bio8(最濕季平均氣溫)、bio19(最冷季降水量)、bio13(最濕月降水量)、bio15(降水量的季節(jié)性變化)、bio16(最干季降水量)、bio5(極端最低氣溫)、bio10(最熱季平均氣溫)、bio12(年降水量)、bio2(平均周溫度變化范圍)。結(jié)果表明:極端最高氣溫、最干月降水量、年均氣溫變化范圍、溫度季節(jié)性變化和等溫性對澳洲堅果分布的影響較大。

    圖2 澳洲堅果在我國主要適生省份或地區(qū)分布圖Fig.2 Distribution of macadamia nuts in suitable provinces of China

    圖3 澳洲堅果ROC曲線Fig.3 Macadamia ROC curve

    圖4 刀切法對環(huán)境氣候因子檢測Fig.4 Jackknife test for environmental variable significance performed by Maxent

    4 討論

    利用MaxEnt模型與GIS的定量分析,對澳洲堅果在世界范圍內(nèi)以及在我國的潛在適生區(qū)進行預(yù)測,直觀且定量的獲得了澳洲堅果的潛在分布區(qū),預(yù)測結(jié)果經(jīng)ROC曲線分析評價,得到模型的AUC值為0.993,預(yù)測效果非常好。預(yù)測結(jié)果顯示:亞洲地區(qū)適宜引種種植澳洲堅果的國家為印度、緬甸、中國、越南、尼泊爾等地。印度[16]和中國[17]等地已成功引種,說明分析預(yù)測的準(zhǔn)確性。預(yù)測結(jié)果還表明,澳洲堅果在我國的高適宜區(qū)和較適宜生長區(qū)累計面積是全球最大的,說明澳洲堅果在我國的發(fā)展空間潛力最大,尤其是在南部沿海地區(qū)。預(yù)測范圍內(nèi)的廣西桂中[18]、云南普洱[17]等地的成功引種澳洲堅果,再次證實了該預(yù)測的準(zhǔn)確性。

    目前,澳洲堅果在美國的夏威夷、加州,澳大利亞,南非,巴西,哥斯達(dá)黎加,以色列等國家和地區(qū)有澳洲堅果的商業(yè)生產(chǎn)地。這些生產(chǎn)區(qū)也基本在本研究劃定范圍內(nèi),成功預(yù)測的美國佛羅里達(dá)州,將澳洲堅果作為庭院植物[19];但預(yù)測并沒有把美國夏威夷和加州等實際商業(yè)生產(chǎn)地預(yù)測在內(nèi),說明預(yù)測與現(xiàn)實存在一定的誤差,抑或通過園藝措施可以彌補部分氣候的影響,擴大了其適生的范圍,導(dǎo)致差異的產(chǎn)生。這些結(jié)果和現(xiàn)象說明預(yù)測的范圍窄于其可能的適生范圍,以及通過園藝措施改善一些氣象條件,可以擴大其適生范圍。

    5 結(jié)論

    澳洲堅果生長區(qū)域相對狹窄,對環(huán)境要求苛刻。在世界范圍內(nèi),在澳洲東部、南美洲東南部和馬達(dá)加斯加島東部較適宜澳洲堅果生長,較大面積的高適宜區(qū)位于亞洲23°26′~30°N, 73°~122°E范圍內(nèi)。

    在我國西藏、臺灣、廣西、廣東和云南澳洲堅果的高適宜區(qū)和較適宜生長區(qū)面積較大,高適宜區(qū)面積為47 216 km2,較適宜生長區(qū)面積為90 391 km2,高適宜區(qū)和較適宜生長區(qū)累計面積是全球最大的地區(qū)。澳洲堅果在我國的發(fā)展空間很大。

    通過MaxEnt最大熵模型對19個環(huán)境氣候因子進行檢測發(fā)現(xiàn),平均周溫度變化對澳洲堅果分布影響最小,極端最高氣溫、降水量、平均氣溫變化范圍、溫度季節(jié)性變化和等溫性這5個氣候因子對澳洲堅果分布影響最大。

    MaxEnt最大熵模型對澳洲堅果預(yù)測結(jié)果進行ROC曲線檢驗,檢驗結(jié)果其AUC>0.9。說明模型對澳洲堅果的預(yù)測結(jié)果較高,若在國內(nèi)引種和推廣澳洲堅果,應(yīng)在高適宜區(qū)和較適宜生長區(qū)進行引種試驗或者推廣,低適生區(qū)和不適生區(qū)需謹(jǐn)慎引種。

    綜上所述,研究采集了693個澳洲堅果的分布數(shù)據(jù),基本考慮到絕大多數(shù)澳洲堅果的實際生態(tài)位,明確劃定了澳洲堅果的具體適生區(qū),為澳洲堅果引種,尤其為我國澳洲堅果引種和推廣提供了參考。由于實際情況的復(fù)雜性,研究中僅考慮了氣候和地區(qū)的重要影響因素,而沒有考慮澳洲堅果對光照以及土壤等條件的需求。在引種澳洲堅果時,既要參考本文的區(qū)劃結(jié)果也應(yīng)考察土壤條件等,以增加引種的成功率。另外,澳洲堅果根系不粗壯,抓地不牢,不耐臺風(fēng),在廣西等有較強臺風(fēng)的地區(qū)也限制其發(fā)展。綜合考慮多方面的因素,結(jié)合本論文的分布劃分,可以提高引種和推廣的成功率[2]。模型預(yù)測劃分的潛在地理分布范圍之外,盡管引種成功的概率低,但通過采取相應(yīng)的措施不排除實現(xiàn)成功引種栽培的可能。

    [1] 劉建福,黃 莉. 澳洲堅果的營養(yǎng)價值及其開發(fā)利用[J]. 中國食物與營養(yǎng), 2005, 16(2): 25-26.

    [2] 葉振喜,陳作泉,李仍然,等. 澳洲堅果引種試種研究初報[J]. 熱帶作物學(xué)報, 1995, 16(2): 70-77.

    [3] Phillips S J, Anderson R P, Schapire R E. Maximum entropy modeling of species geographic distributions[J]. Ecological modelling, 2006, 190(3): 231-259.

    [4] Townsend Peterson A, Papes M, Eaton M. Transferability and model evaluation in ecological niche modeling: a comparison of GARP and Maxent[J]. Ecography, 2007, 30(4): 550-560.

    [5] 張海娟,陳 勇,黃烈健,等. 基于生態(tài)位模型的薇甘菊在中國適生區(qū)的預(yù)測[J]. 農(nóng)業(yè)工程學(xué)報, 2011, 27(S1): 413-418.

    [6] 馬松梅,張明理,張宏祥,等. 利用最大熵模型和規(guī)則集遺傳算法模型預(yù)測孑遺植物裸果木的潛在地理分布及格局[J]. 植物生態(tài)學(xué)報, 2010, 34(11): 1327-1335.

    [7] 曹向鋒,錢國良,胡白石,等. 采用生態(tài)位模型預(yù)測黃頂菊在中國的潛在適生區(qū)[J]. 應(yīng)用生態(tài)學(xué)報, 2010, 21(12): 3063-3069.

    [8] 馬松梅,聶迎彬,段 霞,等. 蒙古扁桃植物的潛在地理分布及居群保護優(yōu)先性[J]. 生態(tài)學(xué)報, 2015, 35(9): 2960-2966.

    [9] 龔 維,夏 青,陳紅鋒,等. 珍稀瀕危植物伯樂樹的潛在適生區(qū)預(yù)測[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報, 2015, 36(4): 98-104.

    [10] 鄧 飛,李曉兵,王 宏,等. 基于MaxEnt模型評價紫花苜蓿在錫林郭勒盟的分布適宜性及主導(dǎo)因子[J]. 草業(yè)科學(xué), 2014, 31(10): 1840-1847.

    [11] 何淑婷,白碧玉,但佳惠,等. 基于MaxEnt的南丹參在中國的潛在分布區(qū)預(yù)測及適生性分析[J]. 安徽農(nóng)業(yè)科學(xué), 2014, 42(8): 2311-2314.

    [12] 王雷宏,楊俊仙,徐小牛. 基于MaxEnt分析金錢松適生的生物氣候特征[J]. 林業(yè)科學(xué), 2015, 51(1): 127-131.

    [13] 胡 秀,楊劍文,高麗霞,等. 基于 MaxEnt 生態(tài)學(xué)模型對普洱姜花潛在的園林引種區(qū)進行預(yù)測[J]. 廣東園林, 2014, 36(3): 65-68.

    [14] 雷軍成,徐海根. 基于MaxEnt的加拿大一枝黃花在中國的潛在分布區(qū)預(yù)測[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2010, 26(2): 137-141.

    [15] 胡 秀,吳福川,郭 微,等. 基于MaxEnt生態(tài)學(xué)模型的檀香在中國的潛在種植區(qū)預(yù)測[J]. 林業(yè)科學(xué), 2014, 50(5): 27-33.

    [16] Singh R V, Chand D, Tyagi V,etal. Priorities for Introduction of fruit crops in India[J]. Indian Journal of Plant Genetic Resources, 2005, 18(1): 67-68.

    [17] 周 元. 試論澳洲堅果在思茅引種試種栽培前景[J]. 資源科學(xué), 1999, 21(3): 83-86.

    [18] 陳顯國,黃 華,林玉虹,等. 澳洲堅果在廣西桂中地區(qū)的適應(yīng)性研究[J]. 熱帶農(nóng)業(yè)科學(xué), 2013, 33(7): 33-37.

    [19] Boning C. Florida's Best Fruiting Plants: Native and Exotic Trees, Shrubs, and Vines[Z]. Florida: Pineapple Press Inc. ISBN, 2006.

    (責(zé)任編輯:詹春梅)

    Predicting Potential Geographic Suitable Regions ofMacadamiaintegrifoliaBased on MaxEnt

    WANG Wei1, TIAN Rong-rong2, NA Li-yan1,2, YANG Ying3, XU Xin-qiao1

    (1.State Key Laboratory of Tree Genetics and Breeding, Key Laboratory of Tree Breeding and Cultivation of State Forestry Administration, Research Institute of Forestry, Chinese Academy of Forestry, Beijing 100091, China; 2.Biomedicine Department, Beijing City University, Beijing 100094, Chinia; 3.Shandong Institute for Food and Drug Control, Ji’nan 250101, Shandong, China)

    [Objective]To predict the potential distribution area ofMacadamiaintegrifoliain China and analyze its suitability. [Method]The distribution ofM.integrifoliaand 19 related environmental factors were collected. The potential suitable distribution range ofM.integrifoliain the world was analyzed by the softwares MaxEnt and ArcGIS. [Result]The results showed that the growth areas is very narrow,and this species is demanding on the growth environment.In the world, macadamia nuts is suitable to growth in eastern Australia, southeast of South America and the eastern part of Madagascar Island, but the suitable area are mainly located in the range of latitude 23°26′~30° N and longitude 73°~122°E in Asia. In China, the suitable distribution areas forM.integrifoliaare mainly concentrated in Xizang, Taiwan, Guangdong, Guangxi and Yunnan. The highly suitable area in the province mentioned above are: Xizang (15 359 km2), Taiwan (14 054 km2), Guangdong (7 372 km2), Guangxi (6 147 km2) and Yunnan (3 776 km2). Jackknife analysis showed that five environmental factors (maximum temperature of the warmest month, mean annual temperature range, precipitation of the driest month, temperature seasonality, Isothermality) have obvious influence onM.integrifoliadistributions. [Conclusion]In general, MaxEnt accurately simulated the geographical distribution ofM.integrifolia. It showed the basic pattern of geographic distribution and the potential distribution areas, and clarified the dominant bioclimatic factors to geographic distribution ofM.integrifolia. The results provide references forM.integrifoliaintroduction and extension in China.

    Macadamiaintegrifolia; potential geographical distribution; MaxEnt; ArcGIS

    10.13275/j.cnki.lykxyj.2017.03.012

    2015-12-14

    國家自然科學(xué)基金(41501059);中央級公益性科研院所基本科研業(yè)務(wù)費專項資金(CAFYBB2016QB004)

    王偉(1980—),男,山東臨沂人,助理研究員.主要研究方向:經(jīng)濟林栽培與加工.Email: greatkingww@126.com

    * 通訊作者:許新橋(1972—),男,山東菏澤人,教授級高級工程師.主要研究方向: 經(jīng)濟林.電話:010-62889603,E-mail: greatkingww@163.com

    S727.3

    A

    1001-1498(2017)03-0444-06

    猜你喜歡
    適生區(qū)氣候因子堅果
    堅果不見了
    青海省果洛州高原鼠兔宜生區(qū)分布及其對氣候變化的響應(yīng)
    氣候變化下中國蒟蒻薯科箭根薯的地理分布格局預(yù)測
    未來氣候條件下當(dāng)歸適生區(qū)預(yù)測及時空變化分析
    搖一搖堅果罐
    氣候變化下瀕危植物半日花在中國的潛在分布
    巴拉圭瓜多竹適生區(qū)分布研究
    愛吃堅果的小松鼠
    啟蒙(3-7歲)(2018年10期)2018-10-13 02:19:26
    四川水稻產(chǎn)量及其構(gòu)成要素對不同生育期氣候因子的響應(yīng)分析
    天山北坡NDVI對氣候因子響應(yīng)的敏感性分析
    国产精品香港三级国产av潘金莲 | 亚洲精品美女久久av网站| 亚洲,欧美,日韩| 国产在线视频一区二区| 午夜福利在线免费观看网站| 亚洲人成电影观看| 精品一区二区三区四区五区乱码 | 亚洲欧洲日产国产| 久久人人爽av亚洲精品天堂| 又黄又粗又硬又大视频| 麻豆精品久久久久久蜜桃| 免费黄网站久久成人精品| 国产又爽黄色视频| 成年人免费黄色播放视频| 男女之事视频高清在线观看 | 纵有疾风起免费观看全集完整版| 99热国产这里只有精品6| 亚洲av福利一区| 国产亚洲欧美精品永久| 国产亚洲欧美精品永久| 久久性视频一级片| 老熟女久久久| 午夜影院在线不卡| 午夜福利免费观看在线| 日韩 亚洲 欧美在线| 亚洲av国产av综合av卡| 午夜福利视频在线观看免费| 日韩一卡2卡3卡4卡2021年| 天天躁夜夜躁狠狠躁躁| 1024视频免费在线观看| 黑人欧美特级aaaaaa片| 国产精品久久久av美女十八| 日本av免费视频播放| 日日啪夜夜爽| 在线观看国产h片| 纯流量卡能插随身wifi吗| 精品国产乱码久久久久久男人| 一边亲一边摸免费视频| 亚洲精品一二三| 青春草亚洲视频在线观看| 久久久欧美国产精品| 久久性视频一级片| 啦啦啦啦在线视频资源| 日韩制服骚丝袜av| 亚洲精品国产av蜜桃| 亚洲美女黄色视频免费看| 日韩欧美一区视频在线观看| 美女扒开内裤让男人捅视频| 少妇人妻久久综合中文| 亚洲av欧美aⅴ国产| 国产午夜精品一二区理论片| 精品亚洲成国产av| 妹子高潮喷水视频| 国产在视频线精品| 日本一区二区免费在线视频| 国产福利在线免费观看视频| 在线亚洲精品国产二区图片欧美| 久久综合国产亚洲精品| 精品国产乱码久久久久久小说| 欧美日本中文国产一区发布| 在线天堂最新版资源| 国产免费福利视频在线观看| 香蕉丝袜av| 捣出白浆h1v1| 欧美av亚洲av综合av国产av | 国产精品 国内视频| 国产av一区二区精品久久| 久久人人97超碰香蕉20202| a 毛片基地| 另类精品久久| 亚洲欧美日韩另类电影网站| 亚洲精品中文字幕在线视频| 久久久国产欧美日韩av| www.自偷自拍.com| 国产日韩欧美视频二区| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费现黄频在线看| 国产深夜福利视频在线观看| 丰满乱子伦码专区| 国产精品一国产av| 成年人免费黄色播放视频| 亚洲第一av免费看| 天天添夜夜摸| 国产国语露脸激情在线看| 色吧在线观看| 最近手机中文字幕大全| 99热全是精品| 丝袜脚勾引网站| 中文字幕av电影在线播放| 91aial.com中文字幕在线观看| 久久毛片免费看一区二区三区| 少妇人妻精品综合一区二区| 看免费av毛片| 无限看片的www在线观看| 黄色视频在线播放观看不卡| 国产乱来视频区| 人妻 亚洲 视频| 久久久久精品人妻al黑| 免费观看人在逋| 亚洲国产av影院在线观看| 国产亚洲欧美精品永久| av在线播放精品| 国语对白做爰xxxⅹ性视频网站| 日本欧美国产在线视频| 建设人人有责人人尽责人人享有的| 久久婷婷青草| 午夜福利在线免费观看网站| 国产熟女午夜一区二区三区| 国产免费视频播放在线视频| 国产亚洲一区二区精品| 一个人免费看片子| 国产成人一区二区在线| av免费观看日本| 亚洲一级一片aⅴ在线观看| 一级毛片我不卡| 最近的中文字幕免费完整| 午夜影院在线不卡| 亚洲中文av在线| 大香蕉久久网| 亚洲色图 男人天堂 中文字幕| 99热网站在线观看| 啦啦啦在线观看免费高清www| 午夜福利影视在线免费观看| 久久久久久久大尺度免费视频| 无遮挡黄片免费观看| 深夜精品福利| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产色婷婷电影| 精品卡一卡二卡四卡免费| 一级黄片播放器| 黑人巨大精品欧美一区二区蜜桃| 欧美激情极品国产一区二区三区| 18在线观看网站| 免费不卡黄色视频| 最近手机中文字幕大全| av卡一久久| xxxhd国产人妻xxx| 国产激情久久老熟女| 久久久亚洲精品成人影院| 国产一区二区三区综合在线观看| 国产成人系列免费观看| 亚洲精品中文字幕在线视频| 狂野欧美激情性bbbbbb| 在现免费观看毛片| 日本爱情动作片www.在线观看| 观看av在线不卡| 国产精品久久久久久人妻精品电影 | xxxhd国产人妻xxx| 精品一区二区三区av网在线观看 | 亚洲av电影在线观看一区二区三区| 国产精品成人在线| 亚洲av欧美aⅴ国产| 成人毛片60女人毛片免费| 最新的欧美精品一区二区| 满18在线观看网站| 在线观看免费视频网站a站| 欧美精品亚洲一区二区| 国产精品一区二区在线不卡| 欧美日韩亚洲国产一区二区在线观看 | 人人澡人人妻人| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美精品综合一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 国产在线视频一区二区| 成年美女黄网站色视频大全免费| 国产免费福利视频在线观看| 精品久久蜜臀av无| 成人影院久久| 如日韩欧美国产精品一区二区三区| 日韩中文字幕欧美一区二区 | 亚洲国产欧美日韩在线播放| 人人妻人人爽人人添夜夜欢视频| bbb黄色大片| 国产人伦9x9x在线观看| 婷婷成人精品国产| 久热这里只有精品99| 国产男女超爽视频在线观看| 日本色播在线视频| 欧美国产精品一级二级三级| www日本在线高清视频| 妹子高潮喷水视频| 色网站视频免费| 午夜av观看不卡| 国产精品欧美亚洲77777| 狂野欧美激情性bbbbbb| 成年动漫av网址| 日韩一区二区视频免费看| 日韩av不卡免费在线播放| 18禁裸乳无遮挡动漫免费视频| 99久国产av精品国产电影| 国产欧美日韩一区二区三区在线| 亚洲久久久国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦中文免费视频观看日本| 999久久久国产精品视频| 亚洲国产精品999| 五月开心婷婷网| 最新的欧美精品一区二区| 欧美日韩亚洲国产一区二区在线观看 | 免费黄色在线免费观看| 国产欧美亚洲国产| 91国产中文字幕| 色吧在线观看| 日本欧美国产在线视频| 999精品在线视频| 免费看av在线观看网站| 国产成人欧美| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区四区五区乱码 | 亚洲国产av新网站| 女性被躁到高潮视频| 纯流量卡能插随身wifi吗| 亚洲第一区二区三区不卡| 免费不卡黄色视频| 99久久综合免费| 老司机靠b影院| 日韩一区二区视频免费看| 亚洲人成77777在线视频| 日韩伦理黄色片| 天堂8中文在线网| 精品少妇一区二区三区视频日本电影 | 精品人妻熟女毛片av久久网站| 1024视频免费在线观看| 日韩一卡2卡3卡4卡2021年| 丁香六月欧美| 婷婷成人精品国产| 国产伦理片在线播放av一区| 国产97色在线日韩免费| 99久国产av精品国产电影| 日本91视频免费播放| 久久久国产精品麻豆| av片东京热男人的天堂| av电影中文网址| av线在线观看网站| 久久女婷五月综合色啪小说| 免费在线观看黄色视频的| 我要看黄色一级片免费的| 亚洲伊人色综图| 午夜福利免费观看在线| 亚洲美女视频黄频| 亚洲熟女毛片儿| 国产午夜精品一二区理论片| 十分钟在线观看高清视频www| 啦啦啦视频在线资源免费观看| 亚洲国产毛片av蜜桃av| 中文天堂在线官网| 亚洲情色 制服丝袜| videos熟女内射| 久久国产精品大桥未久av| 99热国产这里只有精品6| xxxhd国产人妻xxx| 精品一区二区三卡| 日韩人妻精品一区2区三区| 日韩av免费高清视频| 亚洲av电影在线观看一区二区三区| 老司机在亚洲福利影院| 欧美日韩亚洲高清精品| tube8黄色片| a级片在线免费高清观看视频| www.av在线官网国产| 精品国产露脸久久av麻豆| 中文字幕人妻丝袜制服| 久久免费观看电影| 精品国产一区二区久久| 亚洲国产精品一区三区| 超碰97精品在线观看| 亚洲伊人色综图| 欧美97在线视频| 中文字幕人妻丝袜制服| 成人黄色视频免费在线看| 99香蕉大伊视频| 日日撸夜夜添| 亚洲国产精品一区三区| 老司机影院成人| 国产精品一二三区在线看| 亚洲精品国产区一区二| 捣出白浆h1v1| 久久精品国产亚洲av高清一级| 黄片无遮挡物在线观看| 熟女少妇亚洲综合色aaa.| 免费av中文字幕在线| 日日撸夜夜添| 亚洲熟女精品中文字幕| 日韩人妻精品一区2区三区| 丝袜人妻中文字幕| 免费观看a级毛片全部| 久久精品亚洲av国产电影网| 亚洲精品视频女| 中文欧美无线码| 亚洲一级一片aⅴ在线观看| 我要看黄色一级片免费的| 国产精品久久久av美女十八| 欧美日韩成人在线一区二区| 久久精品国产综合久久久| 十八禁高潮呻吟视频| 精品少妇黑人巨大在线播放| 黑人巨大精品欧美一区二区蜜桃| 成人亚洲精品一区在线观看| 久久精品久久久久久噜噜老黄| 久久国产亚洲av麻豆专区| 女人被躁到高潮嗷嗷叫费观| 叶爱在线成人免费视频播放| 国产成人啪精品午夜网站| 国产精品99久久99久久久不卡 | 91aial.com中文字幕在线观看| av在线播放精品| 一个人免费看片子| 汤姆久久久久久久影院中文字幕| 亚洲国产精品一区三区| 婷婷成人精品国产| 亚洲av中文av极速乱| 久久午夜综合久久蜜桃| 欧美成人精品欧美一级黄| 男女下面插进去视频免费观看| 男女高潮啪啪啪动态图| 国产精品 欧美亚洲| 久久 成人 亚洲| 久久久久精品性色| 日韩一区二区三区影片| 男女午夜视频在线观看| 精品视频人人做人人爽| 考比视频在线观看| av有码第一页| 欧美xxⅹ黑人| 超色免费av| 色网站视频免费| 久久女婷五月综合色啪小说| 999久久久国产精品视频| 免费少妇av软件| 免费黄频网站在线观看国产| 成年动漫av网址| 尾随美女入室| 国产av码专区亚洲av| 大陆偷拍与自拍| 国产精品一国产av| 又粗又硬又长又爽又黄的视频| 久久青草综合色| 黄色视频在线播放观看不卡| 欧美xxⅹ黑人| 亚洲一卡2卡3卡4卡5卡精品中文| 一级片'在线观看视频| 日韩中文字幕视频在线看片| 在线观看免费高清a一片| 国产日韩欧美视频二区| 亚洲第一av免费看| 午夜福利一区二区在线看| 国产一区二区激情短视频 | 校园人妻丝袜中文字幕| 9热在线视频观看99| 飞空精品影院首页| 宅男免费午夜| 99国产综合亚洲精品| 1024视频免费在线观看| 国产成人啪精品午夜网站| 美女福利国产在线| 在线免费观看不下载黄p国产| 亚洲少妇的诱惑av| 欧美av亚洲av综合av国产av | 夫妻性生交免费视频一级片| 又大又爽又粗| www.熟女人妻精品国产| 国产精品久久久久久精品古装| 久久精品熟女亚洲av麻豆精品| 捣出白浆h1v1| 一区二区av电影网| 亚洲国产欧美网| 亚洲图色成人| 五月开心婷婷网| 99热国产这里只有精品6| 成人漫画全彩无遮挡| 欧美日韩亚洲国产一区二区在线观看 | 免费少妇av软件| 黄色视频在线播放观看不卡| 狂野欧美激情性bbbbbb| 久久精品aⅴ一区二区三区四区| 免费女性裸体啪啪无遮挡网站| av又黄又爽大尺度在线免费看| 国产成人a∨麻豆精品| 国产亚洲av片在线观看秒播厂| 男女下面插进去视频免费观看| 日韩成人av中文字幕在线观看| 水蜜桃什么品种好| 老汉色av国产亚洲站长工具| 性少妇av在线| 欧美日韩av久久| 久久久久精品人妻al黑| 日韩大片免费观看网站| 韩国高清视频一区二区三区| 午夜福利乱码中文字幕| 男女床上黄色一级片免费看| 欧美国产精品va在线观看不卡| 青春草国产在线视频| 亚洲精品国产av成人精品| 18禁动态无遮挡网站| 久久av网站| 成人毛片60女人毛片免费| 精品国产一区二区久久| xxx大片免费视频| 日日摸夜夜添夜夜爱| 亚洲三区欧美一区| 日韩欧美一区视频在线观看| videosex国产| 丝袜人妻中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 天堂俺去俺来也www色官网| 亚洲精华国产精华液的使用体验| 日本av免费视频播放| 国产免费视频播放在线视频| 午夜福利一区二区在线看| 久久精品久久久久久噜噜老黄| 亚洲精品日本国产第一区| 亚洲精品美女久久久久99蜜臀 | 国产免费又黄又爽又色| 久久性视频一级片| 婷婷色综合www| 亚洲欧美成人精品一区二区| 日本wwww免费看| 亚洲欧美激情在线| 国产成人精品在线电影| 青春草亚洲视频在线观看| 新久久久久国产一级毛片| 亚洲激情五月婷婷啪啪| 91老司机精品| 亚洲伊人色综图| 久久免费观看电影| 亚洲国产欧美网| 久久久欧美国产精品| 久久精品久久久久久噜噜老黄| 最近中文字幕高清免费大全6| 国产乱人偷精品视频| 亚洲精品久久久久久婷婷小说| 午夜91福利影院| 九九爱精品视频在线观看| 天天影视国产精品| 欧美日韩精品网址| 欧美日韩亚洲高清精品| 涩涩av久久男人的天堂| 99国产精品免费福利视频| 午夜日本视频在线| 又粗又硬又长又爽又黄的视频| 女人爽到高潮嗷嗷叫在线视频| 如何舔出高潮| 亚洲第一青青草原| 免费观看a级毛片全部| 三上悠亚av全集在线观看| 别揉我奶头~嗯~啊~动态视频 | av一本久久久久| 热99久久久久精品小说推荐| 黄片播放在线免费| www.av在线官网国产| 99久久99久久久精品蜜桃| 免费观看人在逋| a级毛片黄视频| 亚洲一区二区三区欧美精品| 国产一级毛片在线| 日本vs欧美在线观看视频| 欧美精品高潮呻吟av久久| 亚洲国产欧美在线一区| 午夜日本视频在线| 国产片特级美女逼逼视频| 欧美日韩国产mv在线观看视频| 中文字幕制服av| 亚洲精品自拍成人| 精品国产国语对白av| 人人妻人人爽人人添夜夜欢视频| 国产精品欧美亚洲77777| 国产毛片在线视频| 日本猛色少妇xxxxx猛交久久| 日韩中文字幕欧美一区二区 | 一区福利在线观看| 母亲3免费完整高清在线观看| 成人国产av品久久久| 久久久久久久久久久久大奶| 99香蕉大伊视频| 两个人看的免费小视频| 国产精品一二三区在线看| 精品人妻熟女毛片av久久网站| 欧美日韩视频精品一区| 黄色怎么调成土黄色| 老司机靠b影院| www.av在线官网国产| av福利片在线| 99久久99久久久精品蜜桃| 国产精品久久久av美女十八| 国产福利在线免费观看视频| 日本猛色少妇xxxxx猛交久久| 一区二区三区激情视频| 久久精品人人爽人人爽视色| 国产精品成人在线| 最近中文字幕高清免费大全6| 十八禁高潮呻吟视频| 国产探花极品一区二区| 国产欧美日韩综合在线一区二区| 综合色丁香网| 亚洲av在线观看美女高潮| 你懂的网址亚洲精品在线观看| 欧美日本中文国产一区发布| 成人亚洲欧美一区二区av| 精品一区二区免费观看| 美女福利国产在线| 免费少妇av软件| 国产女主播在线喷水免费视频网站| av网站在线播放免费| 亚洲精品第二区| 午夜福利乱码中文字幕| 免费人妻精品一区二区三区视频| 亚洲精品乱久久久久久| 久久人人爽人人片av| 国产精品av久久久久免费| 中文字幕av电影在线播放| 国产不卡av网站在线观看| 一级片'在线观看视频| videos熟女内射| 免费在线观看黄色视频的| 欧美 亚洲 国产 日韩一| 99国产精品免费福利视频| 激情五月婷婷亚洲| 国产免费又黄又爽又色| 成人国产av品久久久| 国产一区亚洲一区在线观看| 超碰成人久久| 成人亚洲欧美一区二区av| 国产精品99久久99久久久不卡 | 国产精品女同一区二区软件| 久久狼人影院| 久久影院123| 国产精品久久久人人做人人爽| 2018国产大陆天天弄谢| 日本欧美视频一区| 日韩一卡2卡3卡4卡2021年| 免费少妇av软件| 日韩一卡2卡3卡4卡2021年| 欧美精品亚洲一区二区| 精品福利永久在线观看| 一级片'在线观看视频| 丰满乱子伦码专区| 亚洲专区中文字幕在线 | 别揉我奶头~嗯~啊~动态视频 | 国产一区二区 视频在线| 一边摸一边做爽爽视频免费| 99国产精品免费福利视频| 国产精品无大码| 丝袜美足系列| 国产日韩欧美亚洲二区| 伊人亚洲综合成人网| 男女高潮啪啪啪动态图| 伊人久久国产一区二区| 精品少妇一区二区三区视频日本电影 | av.在线天堂| 国产亚洲最大av| 纯流量卡能插随身wifi吗| 男女下面插进去视频免费观看| 亚洲情色 制服丝袜| 新久久久久国产一级毛片| 不卡视频在线观看欧美| 多毛熟女@视频| 黄色视频在线播放观看不卡| 女性生殖器流出的白浆| 日本wwww免费看| 国产精品嫩草影院av在线观看| 天天影视国产精品| av在线观看视频网站免费| 无限看片的www在线观看| 成年动漫av网址| 日本av免费视频播放| 欧美 亚洲 国产 日韩一| 精品亚洲成国产av| 国产黄色免费在线视频| 亚洲综合色网址| 国产精品 欧美亚洲| 国产无遮挡羞羞视频在线观看| 亚洲成人国产一区在线观看 | 又粗又硬又长又爽又黄的视频| 久热爱精品视频在线9| 又大又爽又粗| av网站在线播放免费| 不卡av一区二区三区| 午夜福利在线免费观看网站| 国产成人免费无遮挡视频| 午夜福利免费观看在线| 99久久综合免费| 另类亚洲欧美激情| 在线观看三级黄色| 日韩视频在线欧美| 美女扒开内裤让男人捅视频| 十八禁高潮呻吟视频| 亚洲国产欧美在线一区| 欧美97在线视频| 中国三级夫妇交换| 国产又爽黄色视频| 哪个播放器可以免费观看大片| 韩国高清视频一区二区三区| 国产成人精品久久二区二区91 | 人妻人人澡人人爽人人| 波野结衣二区三区在线| 黑人欧美特级aaaaaa片| 午夜福利一区二区在线看| 国产男人的电影天堂91| 少妇被粗大的猛进出69影院| 青青草视频在线视频观看| 国产成人系列免费观看| 亚洲欧洲日产国产| 嫩草影视91久久| 亚洲国产av新网站| 国产亚洲av高清不卡| 观看av在线不卡| 国产日韩欧美视频二区| 日本av免费视频播放| 亚洲国产欧美日韩在线播放| 久久精品亚洲av国产电影网| 亚洲欧美日韩另类电影网站| 黄片无遮挡物在线观看| 18禁观看日本| 中文精品一卡2卡3卡4更新| 欧美变态另类bdsm刘玥|