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

    基于Maxent模型的西南地區(qū)香果樹地理分布預(yù)測研究

    2021-09-24 01:07:56魏俊朱坤陳文德盧雨田蔡艷琨吳群黃仲宣
    中國野生植物資源 2021年8期
    關(guān)鍵詞:適生區(qū)西南地區(qū)表層

    魏俊,朱坤,2,陳文德,2*,盧雨田,蔡艷琨,吳群,黃仲宣

    (1.成都理工大學(xué)旅游與城鄉(xiāng)規(guī)劃學(xué)院,四川成都610059;2.四川省社會科學(xué)重點(diǎn)研究基地(擴(kuò)展)國家公園研究中心,四川成都610059)

    香果樹(Emmenopterys henryiOliv.)隸屬茜草科(Rubiaceae)、香果樹屬(EmmenopterysOliv.)落葉大喬木,是古老孑遺植物。在《中華本草》和《新華本草綱要》中記載香果樹具有濕中和胃,降逆止嘔之功效。香果樹種子萌發(fā)力較低,主要靠根系或樹樁萌條更新,繁殖比較困難,尚未發(fā)現(xiàn)純林,只見其星散分布,野外資源量少,1991 年被認(rèn)定為近危保護(hù)植物納入《中國植物紅皮書》,1999 年被列為第一批中國特有二級珍稀瀕危保護(hù)植物,有研究認(rèn)為香果樹是香果樹屬中己知的唯一現(xiàn)存種[1],具有較高的經(jīng)濟(jì)和科研價(jià)值[2]。

    目前,關(guān)于香果樹主要研究具體保護(hù)區(qū)等小區(qū)域的種群特征、群落結(jié)構(gòu)、生長動態(tài)、基因組特征、繁殖能力和病害等方面且以中東部省區(qū)為主[3-8],對西南地區(qū)香果樹的地理空間范圍分布情況的研究較少,其中張永華對全國香果樹的親緣地理學(xué)和景觀遺傳學(xué)及其生態(tài)位分化進(jìn)行了研究[1];周文娟等對貴州省境內(nèi)香果樹的分布現(xiàn)狀研究發(fā)現(xiàn)其在貴州屬廣布種,但多為零星分布[9]。

    最大熵模型(Maximum entropy model,Maxent)是近年來用作預(yù)測物種在地理尺度上空間分布的重要方法之一。王娟娟等利用Maxent 模型對川貝母的潛在分布和適應(yīng)性進(jìn)行了研究[10];崔晉亮等利用該模型對藍(lán)莓的潛在分布區(qū)進(jìn)行了預(yù)測研究[11];張東方等采用最大熵模型研究了當(dāng)歸在全球的生態(tài)適宜區(qū)和生態(tài)特征[12];王國崢等采用4 種基于不同算法的模型(GARP、Bioclim、Domain 和Maxent)研究金錢松的潛在分布區(qū),結(jié)果表明Maxent 模型對植物潛在地理分布預(yù)測結(jié)果比其他模型更為精確[13]。因此,本研究采用最大熵模型,以對西南地區(qū)涼山山系樣地調(diào)研為基礎(chǔ),結(jié)合香果樹標(biāo)本數(shù)據(jù),對西南地區(qū)(川、云、貴、藏、渝)香果樹在歷史時(shí)期、當(dāng)前時(shí)期和未來時(shí)期適生區(qū)的空間分布情況進(jìn)行預(yù)測分析研究。

    1 數(shù)據(jù)來源與研究方法

    1.1 樣本數(shù)據(jù)來源

    采用樣地調(diào)查法對涼山山系香果樹的分布情況進(jìn)行實(shí)地調(diào)研,獲得樣本數(shù)據(jù)點(diǎn)44 個(gè),在全球生物多樣性信息數(shù)據(jù)庫網(wǎng)絡(luò)(GBIF,http://www.gbif.org/)、中國數(shù)字植物標(biāo)本館(CVH,http://www.cvh.ac.cn/)和中國自然保護(hù)區(qū)資源平臺(http://www.papc.cn/)中共查詢到西南地區(qū)香果樹樣本點(diǎn)117 個(gè),其中貴州 19 個(gè)、云南 5 個(gè)、重慶 7 個(gè)、四川 86 個(gè),數(shù)字植物標(biāo)本館中顯示西藏地區(qū)有1 份香果樹標(biāo)本,但是沒有提供坐標(biāo)信息。加上實(shí)地調(diào)研的數(shù)據(jù)一共得到樣本點(diǎn)161 個(gè),統(tǒng)計(jì)發(fā)現(xiàn)有5 個(gè)重復(fù)樣本點(diǎn),去除后共有156個(gè)樣本點(diǎn)(圖1)。

    圖1 西南地區(qū)香果樹分布情況Fig.1 Distribution of E.henryiin Southwest China

    1.2 生態(tài)因子數(shù)據(jù)來源

    生態(tài)因子共36 個(gè),其中環(huán)境因子包含19 個(gè)生物氣候因子、7 類氣候因子和1 個(gè)地形因子,環(huán)境因子 數(shù) 據(jù) 來 源 于 Worldclim(http://www. worldclim.org/),此數(shù)據(jù)集為當(dāng)前條件(2000~2020 年)且分辨率為30S(約1km2)的生物氣候變量,這些變量主要反映溫度和降水的特點(diǎn)及其季節(jié)性變化特征,地形因子來源于地理空間數(shù)據(jù)云(http://www.gscloud.cn/),下載分辨率為90 m 的數(shù)字高程地圖(DEM),從中提取海拔,然后在ArcGIS10.2 中重采樣生成分辨率為30″的柵格數(shù)據(jù)圖層[14]。土壤數(shù)據(jù)來自聯(lián)合國糧食及農(nóng)業(yè)組織網(wǎng)站中的世界土壤數(shù)據(jù)庫(Har?monized world soil database v1.2|FAO SOILS POR?TAL),坐標(biāo)系為WGS-1984,柵格大小約為1km2,共選取了9個(gè)土壤因子。

    1.3 研究方法

    1.3.1 數(shù)據(jù)預(yù)處理

    表1 生態(tài)因子縮寫及含義Table 1 Abbreviations and meanings of ecological factors

    從國家基礎(chǔ)地理信息系統(tǒng)(http:nfgis.nsdi.gov.cn/)下載獲得的1∶400 000 0 的西南地區(qū)(四川、重慶、云南、貴州、西藏)標(biāo)準(zhǔn)地圖作為分析底圖,再利用ArcGIS 10.2 軟件將西南地區(qū)的生物變量數(shù)據(jù)掩膜出來,將圖層的坐標(biāo)系定義為GCS-WGS-1984,并導(dǎo)出為ASCII格式。將香果樹樣本點(diǎn)數(shù)據(jù)文件保存為“.csv”格式。

    1.3.2 生態(tài)因子篩選

    由于環(huán)境變量之間具有多重共線性,會導(dǎo)致預(yù)測分布過度擬合[13],用 ArcGIS 10.2 中 ArcTool?box 下多元分析的波段集統(tǒng)計(jì)功能計(jì)算生態(tài)因子的協(xié)方差和相關(guān)矩陣進(jìn)行相關(guān)性檢驗(yàn),篩選出相關(guān)系數(shù)絕對值小于0.75 的生態(tài)因子,并依據(jù)香果樹的生長環(huán)境再次篩選,最終得到4 個(gè)土壤因子:表層土礫石含量、表層土粘土含量、表層土壤有機(jī)碳含量、表層土基本飽和度,13 個(gè)環(huán)境因子:海拔、1 月和 7 月平均降水量、1 月太陽輻射、7 月最高溫、1 月最低溫、年均溫、等溫性、溫度季節(jié)性變化的標(biāo)準(zhǔn)差、年均降水量、最干月降水量、降水量變異系數(shù)、最冷季度降水量,共篩選出17 個(gè)生態(tài)因子。

    圖2 生物氣候因子相關(guān)系數(shù)圖Fig.2 Correlation coefficient of bioclimatic factors

    1.3.3 最大熵模型

    最大熵原理(The principle of maximum entropy)起源于信息論和統(tǒng)計(jì)力學(xué),是基于有限的已知信息對未知分布進(jìn)行無偏推斷的一種數(shù)學(xué)方法[15]。Phil?lips 團(tuán)隊(duì)基于生態(tài)位理論,考慮氣候、海拔、植被等環(huán)境因子,用最大熵原理作為統(tǒng)計(jì)推斷工具,構(gòu)建了物種地理尺度上空間分布的最大熵模型[16-17],編寫了一個(gè)可以免費(fèi)獲取的計(jì)算軟件(http://www.cs.princeton.edu/~schapire/maxent/)[15]。該模型是借助機(jī)器學(xué)習(xí)方法建立的對最大熵原理的應(yīng)用[16],Elith等[18]在2011年對其做了統(tǒng)計(jì)上學(xué)的解釋,使這個(gè)模型在生態(tài)學(xué)領(lǐng)域得到了更加廣泛的應(yīng)用。

    最大熵模型軟件的具體操作為:將已知的物種分布數(shù)據(jù)的“.csv”格式文件和篩選出來的生態(tài)因子導(dǎo)入Maxent 軟件,定義結(jié)果輸出的位置以及生態(tài)因子圖層位置,勾選創(chuàng)建響應(yīng)曲線和刀切法復(fù)選框?qū)ψ兞恐匾赃M(jìn)行測試,設(shè)置采取25%的樣本點(diǎn)作為測試集,剩下75%的樣本點(diǎn)作為訓(xùn)練集,其余選項(xiàng)采用模型的默認(rèn)設(shè)定,重復(fù)運(yùn)行10次進(jìn)行建模。受試者工作特征曲線(Receiver operating characteristic?curve,ROC)分析法在物種潛在分布預(yù)測模型評價(jià)中應(yīng)用比較普遍,這一診斷試驗(yàn)評價(jià)指標(biāo)的認(rèn)可度較高[12,19]。計(jì)算 ROC 曲線下方的面積(Area under curve,AUC),其中AUC 的數(shù)值范圍為 0.6~0.7 時(shí)預(yù)測精度較低,0.7~0.8 時(shí)預(yù)測精度中等,0.8~0.9 時(shí)預(yù)測精度良好,大于0.9 時(shí)預(yù)測精度優(yōu)秀;并利用Jackknife刀切圖來評價(jià)各環(huán)境因子的權(quán)重[14]。

    2 結(jié)果分析

    2.1 預(yù)測精度分析

    將篩選出來的17 個(gè)生態(tài)因子以及香果樹樣本點(diǎn)數(shù)據(jù)分別導(dǎo)入Maxent3.3.3軟件中,重復(fù)運(yùn)行10次上述方法建立最大熵模型,得到生態(tài)因子Maxent 模型平均訓(xùn)練集的AUC 值為0.956(圖3),表明模型預(yù)測結(jié)果為優(yōu)秀,滿足進(jìn)行香果樹地理分布和生態(tài)因子關(guān)系模擬研究的精度要求。

    圖3 生態(tài)因子預(yù)測ROC曲線Fig.3 Receiver operating characteristic curve of ecological factor prediction

    2.2 潛在分布區(qū)及其適生區(qū)等級劃分

    為了對各區(qū)域的適生程度加以區(qū)分,按照政府間氣候變化專門委員會(Intergovernmental panel on climate change,IPCC)關(guān)于可能性的劃分標(biāo)準(zhǔn)并結(jié)合文獻(xiàn)研究[10],將香果樹潛在分布預(yù)測區(qū)劃分為4個(gè)等級:P < 0.1 為不適宜區(qū);0.1 ≤ P < 0.4 低適宜區(qū);0.4 ≤ P < 0.6 為較適宜區(qū);P ≥ 0.6 為高適宜區(qū)。將最大熵模型模擬導(dǎo)出的“Oliv._avg.asc”文件導(dǎo)入ArcGIS 10.2,利用工具欄將“.asc”文件轉(zhuǎn)為柵格文件并進(jìn)行重分類,得到香果樹適生區(qū)分布圖(圖4)。

    圖4 生態(tài)因子模擬香果樹潛在適生區(qū)Fig.4 Ecological factors simulating the potential adaptive area of E.henryi Oliv

    利用ArcGIS10.2 統(tǒng)計(jì)分析工具,將模擬結(jié)果的柵格數(shù)據(jù)轉(zhuǎn)換為矢量數(shù)據(jù)并進(jìn)行面積計(jì)算,得到香果樹當(dāng)前潛在分布區(qū)及其適生區(qū)面積。生態(tài)因子模擬結(jié)果中香果樹分布的高適宜區(qū)面積約為3.27×104km2,中適宜區(qū)面積約為 16.78 × 104km2,低適宜區(qū)面積約為24.13 × 104km2,不適宜區(qū)面積約為198.82 × 104km2,其中適宜區(qū)總面積僅占西南地區(qū)總面積的18.19%。

    2.3 主導(dǎo)生態(tài)因子特征分析

    根據(jù)參與模型建立的生態(tài)因子對最大熵模型的貢獻(xiàn)率(表2),可以判斷影響物種分布的主要生態(tài)因子。生態(tài)因子貢獻(xiàn)率結(jié)果顯示共有8個(gè)環(huán)境因子對模擬結(jié)果的貢獻(xiàn)率大于4%,由高到低依次為表層粘土含量(25.25%)、最干月降水量(14.15%)、海拔(13.51%)、表層土礫石含量(13.1%)、1 月太陽輻射強(qiáng)度(11.69%)、表層土基本飽和度(5%)、1 月最低溫(4.38%)、表層土壤有機(jī)碳含量(4.1%),累積貢獻(xiàn)率高達(dá)91.12%。

    表2 生態(tài)因子貢獻(xiàn)率Table 2 The percent contribution of ecological factors

    2.3.1 香果樹分布區(qū)環(huán)境特征

    基于各氣候因子響應(yīng)曲線(圖5)可知,當(dāng)潛在分布概率大于0.5 時(shí),其對應(yīng)的生態(tài)因子值比較適合香果樹的生存[12],曲線中藍(lán)色部分表示10次重復(fù)運(yùn)算的波動范圍,紅線表示10次運(yùn)算的平均值。最干月降水量小于8 mm 時(shí)香果樹適生概率小于0.5;海拔在950~2 200 m 時(shí)其適生概率大于0.5;1 月太陽輻射在6 900~8 100 kJ/m2·day 時(shí)其適生概率大于0.5;1 月最低溫在-4~1℃時(shí)其適生概率大于0.5,香果樹在極端最低溫-15℃時(shí)都可能生存下來;7 月最高溫在21~32℃時(shí)其適生概率大于0.5,在35℃的高溫下仍具有9%的生存概率;降水量變異系數(shù)的范圍為60~95 時(shí),適生概率大于0.5;等溫性的適宜范圍在25~34。

    圖5 香果樹的潛在分布概率與環(huán)境主導(dǎo)影響因子的關(guān)系Fig.5 Relationship between the potential distribution probability of E.henryi Oliva environmental dominant impact factors

    2.3.2 香果樹分布區(qū)土壤特征

    基于各土壤因子響應(yīng)曲線(圖6)可知適合香果樹生存的土壤因子范圍。表層土粘土含量在21%~40%時(shí)香果樹適生概率大于0.5;表層土礫石含量在4%的時(shí)候其適生概率達(dá)到了0.5,隨著表層土礫石含量升高,適生概率短暫下降,礫石含量達(dá)到5.5%之后其適生概率開始增加,當(dāng)?shù)[石含量增加到8%時(shí)適生概率最大,礫石含量適宜范圍為4%~10%;表層土基本飽和度在15~22之間時(shí)香果樹適生概率逐漸下降,之后香果樹存在概率逐漸增加,當(dāng)表層土基本飽和度達(dá)到33時(shí)香果樹存在概率最大,表層土基本飽和度適宜范圍為10~40;香果樹適生概率隨著表層土壤有機(jī)碳含量的增加逐漸增大,表層土壤有機(jī)碳含量達(dá)到1.5%的時(shí)候香果樹存在概率最大,表層土壤有機(jī)碳含量適宜范圍為0.5%~1.5%。

    圖6 香果樹的潛在分布概率與土壤主導(dǎo)影響因子的關(guān)系Fig.6 Relationship between the potential distribution probability of E.henryi Oliva soil dominant impact factors

    2.4 未來時(shí)空分布預(yù)測

    由于香果樹屬于落葉大喬木,其群落演替時(shí)間較長,因此基于歷史氣候數(shù)據(jù)和未來氣候數(shù)據(jù)對其近100 年的時(shí)空分布變化情況進(jìn)行預(yù)測。氣候數(shù)據(jù)來源于 Worldclim(http://www.worldclim.org/),其中歷史氣候數(shù)據(jù)是該網(wǎng)站發(fā)布的1 970~2 000 年世界氣候2.1 版氣候數(shù)據(jù);未來氣候數(shù)據(jù)(2030s、2050s、2090s)來自于該網(wǎng)站 BCC-CSM2-MR 模型ssp585 共享社會經(jīng)濟(jì)路徑,目前提供的未來氣候數(shù)據(jù)的空間分辨率最高為2.5 min,提供的氣候數(shù)據(jù)類型為月最低溫、月最高溫、月平均降水量和19 個(gè)生物氣候因子。篩選出1 月和7 月平均降水量、1月和7 月最低溫、1 月和7 月最高溫、年均溫、等溫性、溫度季節(jié)性變化的標(biāo)準(zhǔn)差、年均降水量、最干月降水量、降水量變異系數(shù)、最冷季度降水量,假設(shè)海拔因子沒有變化,沿用前面的DEM 數(shù)據(jù),共得到14項(xiàng)環(huán)境因子。

    將4 個(gè)時(shí)期的數(shù)據(jù)按照前面的方式進(jìn)行處理,并分別用Maxent3.3.3 軟件進(jìn)行模擬預(yù)測,平均訓(xùn)練集的 AUC 值分別為 0.942、0.939、0.933、0.930,預(yù)測結(jié)果達(dá)到了優(yōu)秀水平。利用ArcGIS 10.2 進(jìn)行加工處理,得到 20 世紀(jì)末、2030s、2050s和2090s 各時(shí)期香果樹適生區(qū)分布情況(圖7),同樣利用統(tǒng)計(jì)分析工具計(jì)算4 個(gè)時(shí)期香果樹各等級適生區(qū)面積(圖8)、Maxent3.3.3 軟件輸出結(jié)果中包含環(huán)境變量對最大熵模型的相對貢獻(xiàn)的估計(jì)值(圖9)。

    IPCC 第五次評估報(bào)告利用全球氣候耦合模式,模擬發(fā)布了在溫室氣體排放濃度最低情景下,到本世紀(jì)中葉全球平均地表溫度將上升0.4~1.6℃[20]。根據(jù)圖 8 可以看出,在未來的 3 個(gè)時(shí)期,隨著氣候的變化,香果樹在西南地區(qū)的不適宜分布區(qū)面積將會減少,低適宜區(qū)和中適宜區(qū)面積將會有所增加,高適宜區(qū)面積在2030s 和2050s 將會有所增加,但是在2090s 時(shí)期內(nèi)相比于前一時(shí)期將會有所減少。由圖7 可以看出,未來的三個(gè)時(shí)期內(nèi)低適宜區(qū)分布范圍在西藏的東南邊緣和云南的東北部將會有所增加,中適宜區(qū)分布范圍將主要在貴州省境內(nèi)有所增加,2030s 和2050s 高適宜區(qū)分布范圍將在四川的東北邊緣、重慶的東南部和貴州境內(nèi)有所增加,2090s 高適宜區(qū)分布范圍在貴州中西部地區(qū)將會有所減少。

    圖7 不同時(shí)期香果樹適生區(qū)的預(yù)測分布圖Fig.7 Distribution map for predicting the suitable areas of E.henryi Oliv in different periods

    圖8 不同時(shí)期預(yù)測適生區(qū)分布面積Fig.8 Distribution acreage of adaptive area predicted in different periods

    根據(jù)圖9可以看出最干月降水量因子在各時(shí)期貢獻(xiàn)率最大,尤其是2030s、2050s 這兩個(gè)時(shí)期內(nèi)其貢獻(xiàn)率超過了35%;其次是海拔和等溫性因子,在四個(gè)時(shí)期內(nèi)貢獻(xiàn)率都超過了10%;1 月平均降水量、溫度季節(jié)性變化的標(biāo)準(zhǔn)差、7 月最高溫、7 月平均降水量這幾個(gè)因子的貢獻(xiàn)率最小,不超過1%。由此可知,影響香果樹適生區(qū)在西南地區(qū)分布范圍的主要環(huán)境因子是最干月降水量、海拔和等溫性,表明香果樹具有喜濕潤的特點(diǎn)。

    圖9 環(huán)境因子在各時(shí)期貢獻(xiàn)率百分比Fig.9 The percent contribution of environmental factors in each period

    3 結(jié)論

    將 ArcGIS 10.2 與 Maxent3.3.3 結(jié)合,通過 17 個(gè)生態(tài)因子對西南地區(qū)香果樹當(dāng)前的潛在分布區(qū)進(jìn)行預(yù)測,預(yù)測精度達(dá)到優(yōu)秀水平(AUC=0.956),利用ArcGIS 10.2 統(tǒng)計(jì)分析工具計(jì)算得到西南地區(qū)香果樹低適宜區(qū)面積約為21.54×104km2,占西南地區(qū)總面積的9.93%;中適宜區(qū)面積約為16.78 × 104km2,占西南地區(qū)總面積的6.91%;高適宜區(qū)面積約為3.28 × 104km2,占西南地區(qū)總面積的1.35%。香果樹分布中等適宜以上區(qū)域面積占比為8.26%。分析結(jié)果表明西南地區(qū)香果樹的適生區(qū)范圍很小,中高適宜區(qū)主要分布在四川南部、云南東北部、貴州省和重慶市全域范圍內(nèi),分布比較零散,西藏的東南邊緣的唐古拉山區(qū)分布有少量中低適宜區(qū)。

    在影響香果樹分布的主導(dǎo)環(huán)境因子中,最干月降水量、海拔、1 月太陽輻射這三項(xiàng)的貢獻(xiàn)率都超過了10%。適宜區(qū)最干月降水量的閾值為5~43 mm,海拔的閾值為133~3 000 m,1 月太陽輻射的閾值為5 030~9 700 kJ/m2·day,1 月最低溫的閾值為-14℃~7℃。在影響香果樹分布的主導(dǎo)土壤因子中,表層粘土含量、表層土礫石含量的貢獻(xiàn)率超過了10%,其中表層土粘土含量的貢獻(xiàn)率達(dá)到了25.25%;適宜區(qū)表層粘土含量的閾值為15%~56%,表層土礫石含量的閾值為0~15%,表層土基本飽和度的閾值為10%~100%,高適宜區(qū)表層土壤有機(jī)碳含量的閾值為0.2%~2.5%。

    對香果樹在西南地區(qū)近100 年的適生區(qū)的預(yù)測結(jié)果表明,在未來的幾十年內(nèi)香果樹的適生區(qū)面積整體會有小幅度增加,但主要是中低適宜區(qū)取代小部分不適宜區(qū),高適宜區(qū)面積將僅在2050s 有所增加,到2090s 則會減少至現(xiàn)有水平。其中最干月降水量、海拔和等溫性等環(huán)境因子是主要影響因子。

    通過對西南地區(qū)香果樹的潛在分布及適生區(qū)分布的預(yù)測研究,為西南地區(qū)的中國二級珍稀瀕危保護(hù)物種香果樹的樹種分布調(diào)查和保護(hù)提供了較為詳細(xì)的依據(jù);對西南地區(qū)香果樹分布主導(dǎo)生態(tài)影響因子的分析,得出了香果樹生長的生境條件;近100 年的預(yù)測結(jié)果展示了只在氣候環(huán)境變化的影響下,西南地區(qū)香果樹適生區(qū)范圍的變化趨勢,為香果樹的植育保護(hù)提供了參考。

    猜你喜歡
    適生區(qū)西南地區(qū)表層
    氣候變化下中國蒟蒻薯科箭根薯的地理分布格局預(yù)測
    未來氣候條件下當(dāng)歸適生區(qū)預(yù)測及時(shí)空變化分析
    基于MCI指數(shù)的西南地區(qū)近60年夏季干旱特征
    西南地區(qū)干濕演變特征及其對ENSO事件的響應(yīng)
    半潛式平臺表層卡套管處理與認(rèn)識
    海洋石油(2021年3期)2021-11-05 07:43:10
    水體表層沉積物對磷的吸收及釋放研究進(jìn)展
    氣候變化下瀕危植物半日花在中國的潛在分布
    山西隊(duì)?wèi)蛳蛭髂系貐^(qū)的傳播
    中華戲曲(2019年2期)2019-02-06 06:54:16
    巴拉圭瓜多竹適生區(qū)分布研究
    Завершено строительство крупнейшего комплексного транспортного узла на юго
    --западе Китая
    中亞信息(2016年7期)2016-10-20 01:41:30
    欧美成人精品欧美一级黄| av在线天堂中文字幕| avwww免费| 一区福利在线观看| 国内揄拍国产精品人妻在线| 在线观看一区二区三区| 久久久久久伊人网av| a级毛片a级免费在线| 一级毛片我不卡| aaaaa片日本免费| 亚洲一区高清亚洲精品| 国产在视频线在精品| 日日摸夜夜添夜夜添av毛片| 波野结衣二区三区在线| 夜夜看夜夜爽夜夜摸| 成年女人看的毛片在线观看| 日本五十路高清| 天堂动漫精品| 国产国拍精品亚洲av在线观看| 99久久精品热视频| 午夜亚洲福利在线播放| 男女边吃奶边做爰视频| 日韩强制内射视频| 天美传媒精品一区二区| 高清日韩中文字幕在线| 听说在线观看完整版免费高清| 深爱激情五月婷婷| 国产精品久久视频播放| 国产美女午夜福利| 简卡轻食公司| 极品教师在线视频| 天堂影院成人在线观看| 免费一级毛片在线播放高清视频| 久久亚洲国产成人精品v| 欧美绝顶高潮抽搐喷水| 91在线观看av| 啦啦啦韩国在线观看视频| 赤兔流量卡办理| 中文资源天堂在线| 小蜜桃在线观看免费完整版高清| 国产精品99久久久久久久久| 夜夜看夜夜爽夜夜摸| 波多野结衣高清无吗| 午夜影院日韩av| 亚洲欧美成人综合另类久久久 | 国产av麻豆久久久久久久| 久久久久久久久久成人| 中文亚洲av片在线观看爽| 国产真实伦视频高清在线观看| 日日撸夜夜添| 长腿黑丝高跟| 精品午夜福利在线看| 精品人妻偷拍中文字幕| av中文乱码字幕在线| 国产精品,欧美在线| 亚洲精品国产成人久久av| 亚洲精品粉嫩美女一区| 搡老岳熟女国产| 男人舔女人下体高潮全视频| 丰满人妻一区二区三区视频av| 女生性感内裤真人,穿戴方法视频| 俄罗斯特黄特色一大片| 久久综合国产亚洲精品| 免费人成视频x8x8入口观看| 亚洲高清免费不卡视频| 1024手机看黄色片| 夜夜夜夜夜久久久久| 久久久久久久久久久丰满| 国产激情偷乱视频一区二区| 老熟妇乱子伦视频在线观看| 在线观看av片永久免费下载| 久久精品国产自在天天线| 欧美精品国产亚洲| 一区二区三区四区激情视频 | 99热6这里只有精品| 性色avwww在线观看| 91午夜精品亚洲一区二区三区| 香蕉av资源在线| 中国国产av一级| 麻豆国产av国片精品| 免费看日本二区| 中文在线观看免费www的网站| 三级经典国产精品| 国产精品久久久久久亚洲av鲁大| 国产高潮美女av| 国内揄拍国产精品人妻在线| 国产男人的电影天堂91| 中文在线观看免费www的网站| 国产大屁股一区二区在线视频| 亚洲不卡免费看| 亚洲丝袜综合中文字幕| 18+在线观看网站| 国产成人精品久久久久久| 欧美日韩综合久久久久久| 人妻夜夜爽99麻豆av| 欧美日本亚洲视频在线播放| 国产精品三级大全| 色播亚洲综合网| 亚洲欧美日韩高清专用| 能在线免费观看的黄片| 国产精品永久免费网站| 欧美高清成人免费视频www| 在线国产一区二区在线| 黄色配什么色好看| 国产精品一区二区三区四区免费观看 | 国内精品宾馆在线| 久久人人爽人人片av| 极品教师在线视频| 国产精品野战在线观看| 国产在线精品亚洲第一网站| 亚洲性久久影院| 婷婷精品国产亚洲av| 久久久午夜欧美精品| 精品一区二区免费观看| 俄罗斯特黄特色一大片| 男插女下体视频免费在线播放| 国产69精品久久久久777片| 一本精品99久久精品77| 国产午夜精品论理片| 国产熟女欧美一区二区| 国产三级在线视频| 亚洲欧美日韩无卡精品| 国产一区二区激情短视频| 18禁裸乳无遮挡免费网站照片| 啦啦啦啦在线视频资源| 中文字幕久久专区| 亚洲成人av在线免费| 亚洲成人中文字幕在线播放| 一进一出好大好爽视频| 嫩草影院新地址| 1024手机看黄色片| 午夜精品一区二区三区免费看| 在线播放无遮挡| 夜夜看夜夜爽夜夜摸| 国产三级在线视频| АⅤ资源中文在线天堂| 给我免费播放毛片高清在线观看| 亚洲久久久久久中文字幕| 国产av一区在线观看免费| 日本一本二区三区精品| .国产精品久久| 国产成人a区在线观看| aaaaa片日本免费| 欧美成人免费av一区二区三区| av在线老鸭窝| 国产白丝娇喘喷水9色精品| 我的老师免费观看完整版| 久久欧美精品欧美久久欧美| 日韩精品有码人妻一区| 深夜精品福利| 91久久精品国产一区二区三区| 亚洲欧美清纯卡通| 欧美另类亚洲清纯唯美| 观看免费一级毛片| 老女人水多毛片| 中国美白少妇内射xxxbb| 熟女电影av网| 国产视频一区二区在线看| 人妻夜夜爽99麻豆av| 99国产极品粉嫩在线观看| 国产真实伦视频高清在线观看| 日韩欧美免费精品| 成人二区视频| 看黄色毛片网站| 国产极品精品免费视频能看的| 国产成年人精品一区二区| 国产人妻一区二区三区在| 日本三级黄在线观看| 人妻制服诱惑在线中文字幕| 亚洲电影在线观看av| 天堂动漫精品| 男人和女人高潮做爰伦理| 天堂影院成人在线观看| 特级一级黄色大片| 男人舔奶头视频| 国产精品一区二区三区四区久久| 成人无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 午夜日韩欧美国产| 男人舔奶头视频| 天堂动漫精品| 欧美zozozo另类| 麻豆av噜噜一区二区三区| 欧美潮喷喷水| 日本黄色视频三级网站网址| 精品一区二区免费观看| 99九九线精品视频在线观看视频| 国产精品久久久久久av不卡| 久久99热6这里只有精品| 大又大粗又爽又黄少妇毛片口| av天堂中文字幕网| 日产精品乱码卡一卡2卡三| 久久精品影院6| 99riav亚洲国产免费| 晚上一个人看的免费电影| 人人妻人人看人人澡| 日本-黄色视频高清免费观看| 久久久成人免费电影| 我要搜黄色片| 日本黄色片子视频| 欧美高清成人免费视频www| 久久热精品热| 国产黄a三级三级三级人| 国产日本99.免费观看| 国产精品一区二区三区四区免费观看 | 日韩一区二区视频免费看| 亚洲国产精品成人综合色| 久久久久国内视频| 天堂影院成人在线观看| 日本黄大片高清| 哪里可以看免费的av片| 男插女下体视频免费在线播放| 国产成人a区在线观看| 国产视频一区二区在线看| 亚洲av免费在线观看| 男插女下体视频免费在线播放| 国产色爽女视频免费观看| 午夜福利在线观看免费完整高清在 | 欧美人与善性xxx| 精品国产三级普通话版| 国产黄色小视频在线观看| 日本黄色片子视频| 中国美白少妇内射xxxbb| 18禁黄网站禁片免费观看直播| 综合色丁香网| 精品一区二区三区视频在线| 欧美成人免费av一区二区三区| 午夜视频国产福利| 大型黄色视频在线免费观看| 日韩欧美免费精品| 美女高潮的动态| 亚洲中文日韩欧美视频| 麻豆久久精品国产亚洲av| 亚洲av美国av| 免费搜索国产男女视频| 老司机影院成人| 蜜桃久久精品国产亚洲av| 国产白丝娇喘喷水9色精品| 男人舔女人下体高潮全视频| 少妇猛男粗大的猛烈进出视频 | 国产精品伦人一区二区| 久久热精品热| 亚洲va在线va天堂va国产| 日韩在线高清观看一区二区三区| 十八禁网站免费在线| 联通29元200g的流量卡| 五月伊人婷婷丁香| 日本撒尿小便嘘嘘汇集6| 熟女电影av网| 天天躁日日操中文字幕| 伦精品一区二区三区| 亚洲自偷自拍三级| 午夜福利成人在线免费观看| 色在线成人网| 99热6这里只有精品| 亚洲专区国产一区二区| 亚洲av熟女| 精品久久久久久久人妻蜜臀av| 一进一出好大好爽视频| 色av中文字幕| 国产精品福利在线免费观看| 好男人在线观看高清免费视频| 久久久久国产精品人妻aⅴ院| 日本精品一区二区三区蜜桃| 中文字幕精品亚洲无线码一区| 天天一区二区日本电影三级| 久久综合国产亚洲精品| 啦啦啦韩国在线观看视频| avwww免费| 一级黄片播放器| 熟女电影av网| 欧美xxxx黑人xx丫x性爽| 91在线观看av| 女生性感内裤真人,穿戴方法视频| 亚洲国产精品久久男人天堂| 久久欧美精品欧美久久欧美| av视频在线观看入口| 欧美国产日韩亚洲一区| 免费av观看视频| 白带黄色成豆腐渣| 亚洲av成人av| 欧美色视频一区免费| 国产在线男女| 国产成人91sexporn| 国产人妻一区二区三区在| 大香蕉久久网| 天堂动漫精品| 成年免费大片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区亚洲精品在线观看| 成年女人毛片免费观看观看9| 身体一侧抽搐| 国产黄色小视频在线观看| 成人美女网站在线观看视频| 日产精品乱码卡一卡2卡三| 日日摸夜夜添夜夜爱| 又黄又爽又刺激的免费视频.| 在线天堂最新版资源| 久久久久久大精品| 欧美极品一区二区三区四区| 中文亚洲av片在线观看爽| 国产精品精品国产色婷婷| 精品99又大又爽又粗少妇毛片| 色哟哟哟哟哟哟| 观看美女的网站| 十八禁国产超污无遮挡网站| 日韩欧美精品v在线| 丝袜美腿在线中文| 又爽又黄a免费视频| 两个人的视频大全免费| 欧美一区二区亚洲| 天堂影院成人在线观看| 国产三级中文精品| 一区二区三区四区激情视频 | 美女 人体艺术 gogo| 黄色欧美视频在线观看| 亚洲人成网站在线播放欧美日韩| 久久久久久伊人网av| 能在线免费观看的黄片| 看黄色毛片网站| 少妇丰满av| 最近手机中文字幕大全| av在线亚洲专区| 国产一区二区激情短视频| 国产 一区 欧美 日韩| 国产精品av视频在线免费观看| 免费看av在线观看网站| 午夜福利18| 99热全是精品| 国产高潮美女av| 日本欧美国产在线视频| 亚洲一区二区三区色噜噜| 亚洲精品久久国产高清桃花| .国产精品久久| 欧美zozozo另类| 亚洲欧美成人精品一区二区| 黄色视频,在线免费观看| 一级毛片aaaaaa免费看小| 日韩欧美三级三区| 久久久a久久爽久久v久久| 亚洲欧美成人综合另类久久久 | 狂野欧美白嫩少妇大欣赏| 免费在线观看影片大全网站| 午夜福利视频1000在线观看| 国产一级毛片七仙女欲春2| 日韩欧美精品免费久久| 熟女电影av网| 给我免费播放毛片高清在线观看| 成人特级黄色片久久久久久久| 国产精品久久视频播放| 国产精品亚洲一级av第二区| 国产精华一区二区三区| 永久网站在线| 99在线视频只有这里精品首页| 日本-黄色视频高清免费观看| 18禁裸乳无遮挡免费网站照片| av天堂中文字幕网| videossex国产| 国产精品1区2区在线观看.| 一夜夜www| 精品少妇黑人巨大在线播放 | 99在线人妻在线中文字幕| 中国国产av一级| 国产在线精品亚洲第一网站| 国产精品,欧美在线| 国产精品永久免费网站| 国产精品一及| 精品久久久噜噜| 亚洲18禁久久av| 久久久久久久亚洲中文字幕| 亚洲熟妇中文字幕五十中出| a级一级毛片免费在线观看| АⅤ资源中文在线天堂| 亚洲丝袜综合中文字幕| or卡值多少钱| 一个人看视频在线观看www免费| av在线天堂中文字幕| 国产色爽女视频免费观看| av国产免费在线观看| 男女视频在线观看网站免费| 亚州av有码| 精品久久久久久成人av| 色在线成人网| 99热全是精品| 成人av一区二区三区在线看| 99riav亚洲国产免费| 精品福利观看| 亚洲av第一区精品v没综合| 日韩在线高清观看一区二区三区| 一进一出好大好爽视频| 国产精品一二三区在线看| 九九在线视频观看精品| 长腿黑丝高跟| 久久久久精品国产欧美久久久| 色在线成人网| 99热6这里只有精品| 日本免费一区二区三区高清不卡| 男人和女人高潮做爰伦理| .国产精品久久| 日本爱情动作片www.在线观看 | 久久久久国产网址| 日本与韩国留学比较| 18禁裸乳无遮挡免费网站照片| 黄色一级大片看看| 日韩欧美三级三区| 老师上课跳d突然被开到最大视频| 国产不卡一卡二| 免费观看的影片在线观看| 国语自产精品视频在线第100页| 久久综合国产亚洲精品| 一级毛片aaaaaa免费看小| 色综合亚洲欧美另类图片| 久久久久性生活片| 亚洲人成网站在线播放欧美日韩| 日本免费a在线| a级一级毛片免费在线观看| 国产精品一及| 女的被弄到高潮叫床怎么办| 国内少妇人妻偷人精品xxx网站| 毛片一级片免费看久久久久| 亚洲成av人片在线播放无| 亚洲精品色激情综合| 欧美另类亚洲清纯唯美| 欧美又色又爽又黄视频| 99久国产av精品国产电影| 黄色配什么色好看| 中文字幕熟女人妻在线| 精品久久久噜噜| 香蕉av资源在线| 精品午夜福利视频在线观看一区| 亚洲精品亚洲一区二区| 好男人在线观看高清免费视频| 精品一区二区三区人妻视频| 最近的中文字幕免费完整| 免费看a级黄色片| 国产精品久久久久久久久免| 草草在线视频免费看| 国产乱人偷精品视频| 日韩三级伦理在线观看| 真人做人爱边吃奶动态| 欧美一区二区亚洲| 婷婷精品国产亚洲av| 亚洲va在线va天堂va国产| 国内精品美女久久久久久| 国产精品一区www在线观看| 免费搜索国产男女视频| 日韩欧美在线乱码| 国产高清有码在线观看视频| 女同久久另类99精品国产91| 国产精品嫩草影院av在线观看| 国产免费一级a男人的天堂| 最近2019中文字幕mv第一页| 亚洲国产精品sss在线观看| 午夜精品国产一区二区电影 | 日韩制服骚丝袜av| 午夜久久久久精精品| 亚洲av成人av| 亚洲成av人片在线播放无| 欧美一区二区亚洲| .国产精品久久| 男人和女人高潮做爰伦理| 亚洲激情五月婷婷啪啪| 国产v大片淫在线免费观看| 丰满的人妻完整版| 日日摸夜夜添夜夜添小说| 我要看日韩黄色一级片| 又黄又爽又免费观看的视频| 午夜免费激情av| 国产精品综合久久久久久久免费| 99热这里只有是精品在线观看| 亚洲精品成人久久久久久| 久久精品国产亚洲av天美| 久久久久久久久久黄片| 国产亚洲精品综合一区在线观看| 亚洲国产精品成人综合色| 国产视频一区二区在线看| 好男人在线观看高清免费视频| 午夜福利在线在线| 国产91av在线免费观看| 亚洲精品色激情综合| 看黄色毛片网站| 岛国在线免费视频观看| 夜夜看夜夜爽夜夜摸| 亚洲人与动物交配视频| 婷婷色综合大香蕉| 老女人水多毛片| 床上黄色一级片| 国产精品99久久久久久久久| 成人无遮挡网站| 99久久久亚洲精品蜜臀av| 夜夜夜夜夜久久久久| 久久亚洲精品不卡| 中文字幕久久专区| 俄罗斯特黄特色一大片| 成人亚洲欧美一区二区av| 成人性生交大片免费视频hd| 国产 一区精品| 亚洲第一电影网av| 男女边吃奶边做爰视频| 久久久久国产精品人妻aⅴ院| 一级毛片我不卡| 大型黄色视频在线免费观看| 午夜精品国产一区二区电影 | 国产黄色小视频在线观看| 亚洲人成网站在线播放欧美日韩| 秋霞在线观看毛片| 久久精品久久久久久噜噜老黄 | 亚洲第一电影网av| 欧美一区二区国产精品久久精品| 国产毛片a区久久久久| 日韩欧美精品免费久久| 国产精品久久久久久亚洲av鲁大| 亚洲精品乱码久久久v下载方式| 国产私拍福利视频在线观看| 欧美日本视频| 一区二区三区免费毛片| 日韩欧美国产在线观看| 亚洲18禁久久av| 亚洲性久久影院| 97热精品久久久久久| ponron亚洲| 中文字幕免费在线视频6| 国产国拍精品亚洲av在线观看| 精品少妇黑人巨大在线播放 | 日日摸夜夜添夜夜爱| 精品一区二区三区人妻视频| 国产v大片淫在线免费观看| 国产高清视频在线观看网站| 免费看光身美女| 男人和女人高潮做爰伦理| 日韩欧美精品v在线| 国产高清三级在线| 亚洲av一区综合| 在线天堂最新版资源| 99热这里只有是精品在线观看| 国产精品久久久久久精品电影| 国产伦精品一区二区三区四那| 一区二区三区四区激情视频 | 欧美一区二区国产精品久久精品| h日本视频在线播放| 熟女人妻精品中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 一级毛片aaaaaa免费看小| 欧美zozozo另类| 欧美一区二区精品小视频在线| 中文字幕免费在线视频6| 国产成人福利小说| 精品福利观看| 精品无人区乱码1区二区| 亚洲美女视频黄频| 中文字幕av在线有码专区| 天天躁日日操中文字幕| 人人妻,人人澡人人爽秒播| 中文在线观看免费www的网站| 久久久精品欧美日韩精品| 插逼视频在线观看| 干丝袜人妻中文字幕| 久久久久久久久久黄片| 日韩国内少妇激情av| 午夜精品在线福利| 久久精品国产自在天天线| 国产激情偷乱视频一区二区| 夜夜爽天天搞| 91午夜精品亚洲一区二区三区| 国产久久久一区二区三区| 一级毛片电影观看 | 国产大屁股一区二区在线视频| 久久久精品欧美日韩精品| 你懂的网址亚洲精品在线观看 | 在线免费观看不下载黄p国产| 国产久久久一区二区三区| 亚洲精品国产成人久久av| 在线a可以看的网站| 少妇熟女欧美另类| 91麻豆精品激情在线观看国产| 亚洲性久久影院| 欧美色视频一区免费| 男女下面进入的视频免费午夜| 成人精品一区二区免费| 俺也久久电影网| 亚洲精品日韩av片在线观看| 国产探花在线观看一区二区| 毛片一级片免费看久久久久| 人人妻,人人澡人人爽秒播| 又黄又爽又刺激的免费视频.| 91在线精品国自产拍蜜月| 老熟妇仑乱视频hdxx| 99久久成人亚洲精品观看| 在线观看美女被高潮喷水网站| 免费在线观看影片大全网站| 亚洲不卡免费看| 99九九线精品视频在线观看视频| 亚洲国产高清在线一区二区三| 国产精品伦人一区二区| 女生性感内裤真人,穿戴方法视频| 变态另类成人亚洲欧美熟女| 麻豆av噜噜一区二区三区| 国产精品人妻久久久影院| 又黄又爽又免费观看的视频| 一级黄片播放器| 精品乱码久久久久久99久播| 亚洲欧美日韩东京热| 三级毛片av免费| 亚洲精品色激情综合| 老熟妇乱子伦视频在线观看| 在线免费观看不下载黄p国产| 久久久久久久久久黄片| 国产精品一区二区三区四区免费观看 | 美女高潮的动态| 国产黄a三级三级三级人| av天堂中文字幕网| 久久久国产成人精品二区| 色哟哟·www| 夜夜看夜夜爽夜夜摸| 精品久久久久久成人av|