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

    新疆典型蝗蟲適生區(qū)分布預(yù)測

    2022-11-26 08:17:44巨喜鋒吳建國陳吉軍關(guān)靖云李漠巖鄭江華
    生態(tài)學(xué)報 2022年21期
    關(guān)鍵詞:物種新疆模型

    巨喜鋒,林 峻,吳建國,陳吉軍,關(guān)靖云,李漠巖,鄭江華,,*

    1 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院, 烏魯木齊 830046

    2 新疆治蝗滅鼠指揮部辦公室, 烏魯木齊 830001

    3 綠洲生態(tài)教育部重點實驗室, 烏魯木齊 830046

    草地是自然生態(tài)系統(tǒng)的重要組成部分,不僅在水土保持、生態(tài)平衡、碳固定方面具有重要生態(tài)價值,而且在政治、經(jīng)濟、文化、生態(tài)環(huán)境建設(shè)方面扮演著重要角色[1—2]。新疆草原草地面積為57.2588萬km2,在全國草地面積中排名第三[3],其健康受包括蝗災(zāi)在內(nèi)的各類草原生物災(zāi)害嚴重影響,其中蝗災(zāi)與水災(zāi)、旱災(zāi)并稱為三大天災(zāi),新疆蝗災(zāi)發(fā)生最早可追溯到東漢時期[4],1986—1991年發(fā)生面積近4萬km2,其中嚴重危害區(qū)面積近2.8萬km2,發(fā)生面積占可利用面積的8.5%[5]。新疆草原蝗蟲具有種類多,分布廣,數(shù)量大,災(zāi)害嚴重的特點。20世紀陳永林[6]調(diào)查新疆蝗蟲有62屬157種;近期于冰潔等[7]指出新疆蝗蟲種類達65屬177種;而徐超民等[8]認為新疆草原蝗蟲發(fā)生區(qū)蝗蟲記載有8科16屬171種,其中優(yōu)勢種包括亞洲飛蝗(Locustamigratoriamigratoria)、意大利蝗(CalliptamusItalicus)、黑腿星翅蝗(C.barbarus)、西伯利亞蝗(GomphocerusSibirius)、寬須蟻蝗(Myrmeleotettixpalpalis)、黑條小車蝗(Oedaleusdecorus)、紅脛戟紋蝗(Dociostauruskraussikraussi)等十余種。近年新疆極端氣候事件日漸頻繁,氣候變暖和極端降水增加與新疆典型蝗蟲持續(xù)爆發(fā)關(guān)系緊密[9],此外肆虐亞非的沙漠蝗(Schistocercagregaria),使得全球糧食安全和生態(tài)安全備受關(guān)注,文獻也應(yīng)用多種技術(shù)方法研究了沙漠蝗的遷飛通道以及入境中國的風(fēng)險,此外Guan等還通過MaxEnt模型研究了在未來氣候條件下沙漠蝗的適生區(qū)分布情況[10—11]。在常態(tài)化疫情防控條件下,切實做好新疆本地典型蝗蟲適生區(qū)預(yù)測與境外蝗蟲遷飛入境監(jiān)測工作異常重要。新疆草原蝗蟲在日常監(jiān)測、防治方面面臨防治范圍廣、防治監(jiān)測蝗蟲類型多的嚴峻挑戰(zhàn)。

    物種分布模型(Species Distribution Models, SDMs)以已知物種分布點及環(huán)境數(shù)據(jù)作為輸入,利用數(shù)學(xué)統(tǒng)計方法或機器學(xué)習(xí)理論探究物種在生態(tài)位空間實際狀況[12],在生物多樣性保護研究、物種入侵研究及物種適生區(qū)預(yù)測等方面得到廣泛應(yīng)用[13]。在模型發(fā)展中形成了基于生態(tài)位理論、數(shù)理統(tǒng)計和機器學(xué)習(xí)物種分布模型[14],其中BIOCLIM模型、HABITAT模型及領(lǐng)域模型(DOMAIN)等都是典型的基于生態(tài)位量化理論而形成的模型,該類模型基于典型已知物種分布點及環(huán)境空間中多變量極值組成的包絡(luò)體外邊緣,預(yù)測出物種分布的潛在區(qū)域[15]。在此基礎(chǔ)上,隨著數(shù)學(xué)方法的發(fā)展,出現(xiàn)了應(yīng)用數(shù)理統(tǒng)計方法探索物種分布點與預(yù)測區(qū)統(tǒng)計學(xué)關(guān)系以預(yù)測物種潛在分布區(qū)的模型算法,該類模型算法包括廣義線性模型(GLM)、隨機森林模型(RF)和馬氏距離(MAHAL)等[16]。近些年伴隨著計算機科學(xué)的發(fā)展,機器學(xué)習(xí)的理論思想逐步引用到物種分布模型預(yù)測中,其中應(yīng)用相對較為廣泛的有人工神經(jīng)網(wǎng)絡(luò)(ANN)、支持向量機(SVM)、提升回歸樹(BRT)以及最大熵模型(MaxEnt)等一大批優(yōu)秀模型[17]。但基于不同模型的物種分布區(qū)及潛在適生區(qū)預(yù)測結(jié)果具有很強的不確定性,從而導(dǎo)致單個最優(yōu)模型的選擇具有一定困難,因此將多個不同模型預(yù)測輸出組合為一個集成模型變得尤為重要[18]。集成建模的基本思想是將各單一模型中捕獲到的真實信號予以保留,對數(shù)據(jù)和模型中產(chǎn)生的錯誤和不確定更好地分離和降低,從而使得建模的結(jié)果更為可靠[19—20]。

    在新疆典型蝗蟲發(fā)生預(yù)測相關(guān)研究中,楊洪升等[21]探究了新疆蝗蟲與植被長勢、氣候因子間的關(guān)系,指出蝗蟲發(fā)生與植被長勢息息相關(guān),且溫度是影響植被長勢的主要因子。Zheng等[22]使用基于GIS多準則分析方法預(yù)測了新疆蝗蟲可能爆發(fā)的區(qū)域,得到較高的預(yù)測精度。楊會楓等[23]使用MaxEnt模型預(yù)測了西伯利亞蝗在新疆的潛在分布,確定了影響西伯利亞蝗分布的重要環(huán)境因子,得出模型精度較高。李培先等[24]基于MaxEnt模型預(yù)測了未來不同氣候變化情景下新疆意大利蝗潛在分布,指出極端水分條件和水熱條件對其潛在分布影響巨大。因不同物種分布模型自身差異,各模型算法對模型不確定性的影響會進一步影響到模型預(yù)測結(jié)果的準確度[25]?;诩山K枷氚l(fā)展起來的利用多個模型信息的集合模型可以有效降低單一模型預(yù)測的不確定性,從而提升模型的穩(wěn)定性及預(yù)測結(jié)果的可靠性[26]。當(dāng)前也有大量的文獻應(yīng)用多個模型集成對不同生物物種的分布做了預(yù)測,如韓欣嬈等[27]通過多個物種分布模型的集合對草地貪夜蛾的潛在分布區(qū)做了預(yù)測,Lippi等[28]應(yīng)用多個物種分布模型對美國犬蜱的分布區(qū)做了預(yù)測。

    前人的研究集中于單一種類蝗蟲在不同氣候條件下的潛在適生區(qū)展開,但各類主要蝗蟲對地理環(huán)境、草地類型、氣候狀況、土壤條件等環(huán)境的適宜性存在一致性,因此文獻[29]將新疆草原蝗蟲發(fā)生區(qū)分為八大亞區(qū);此外前人針對草原的相關(guān)研究主要通過單一物種分布模型開展,但單一模型預(yù)測的不確定性較高,而通過多模型集成預(yù)測的方法則受到了很大的關(guān)注,相關(guān)研究成果也逐步增多;且相聚焦于與蝗蟲生長密切相關(guān)的氣候因素、植被等因素展開,影響蝗蟲產(chǎn)卵、繁殖關(guān)系密切的土壤環(huán)境研究有待深入。基于此,本文綜合考慮了新疆草原優(yōu)勢蝗種,而非某單一類型蝗蟲;而且集成建??梢杂行Ы档皖A(yù)測不確定性、提升模型的可靠性[30];研究表明土壤因素對蝗蟲的繁殖、蟲卵孵化過程具有顯著影響[31—32]。

    本文以與蝗蟲活動緊密相關(guān)的環(huán)境數(shù)據(jù)為基礎(chǔ),應(yīng)用經(jīng)典物種分布模型對新疆典型蝗蟲的潛在適生區(qū)進行研究,探究模型在新疆典型蝗蟲預(yù)測研究中的性能表現(xiàn);選擇高性能模型構(gòu)建草原蝗蟲預(yù)測模型集合對新疆典型蝗蟲的適生區(qū)展開預(yù)測,以期為一線蝗蟲監(jiān)測工作提供可靠的技術(shù)支持。

    1 研究區(qū)概況

    新疆位于亞歐大陸腹地,處于我國北方干旱半干旱區(qū),地理位置在東經(jīng)73°40′—96°23′,北緯34°25′—49°10′,面積166萬km2,如圖1。其地形特征是“三山夾兩盆”——自北向南依次為阿爾泰山、準噶爾盆地、天山、塔里木盆地和昆侖山。因深居內(nèi)陸,遠離海洋而形成典型的干旱氣候特征,年均氣溫10℃,年均降水量167.1 mm,降水空間差異大,北部多于南部,由此形成多種草地類型。草場面積占自治區(qū)國土面積一半以上,是當(dāng)?shù)剞r(nóng)牧業(yè)重要的生產(chǎn)資料[33]。

    圖1 研究區(qū)位置及新疆典型蝗蟲發(fā)生監(jiān)測點

    2 數(shù)據(jù)與方法

    2.1 數(shù)據(jù)來源

    2.1.1 蝗蟲發(fā)生數(shù)據(jù)

    使用由新疆維吾爾自治區(qū)蝗蟲鼠害預(yù)測預(yù)報防治中心站提供的全疆范圍內(nèi)各類蝗蟲發(fā)生歷史記錄,并通過收集文獻[23—24,33—36]及全球生物多樣性信息網(wǎng)絡(luò)數(shù)據(jù)庫(GBIF)中新疆典型蝗蟲數(shù)據(jù)記錄,此外還有研究團隊長期在野外調(diào)查中采集記錄的蝗蟲發(fā)生區(qū)GPS數(shù)據(jù)。原始數(shù)據(jù)記錄共計3800條,通過篩選及預(yù)處理得到3714條數(shù)據(jù)參與后期處理與分析。

    2.1.2 氣象數(shù)據(jù)

    氣象數(shù)據(jù)包括與新疆典型草原蝗蟲繁殖、生長密切相關(guān)的氣溫降水?dāng)?shù)據(jù),其中氣溫數(shù)據(jù)為月最高氣溫和月最低氣溫。氣象數(shù)據(jù)空間分辨率為2.5′,該數(shù)據(jù)時間范圍基本與所收集的蝗蟲發(fā)生記錄相一致,為2000—2018年。該數(shù)據(jù)集由WorldClim(https://www.worldclim.org/data/monthlywth.html)提供,數(shù)據(jù)包括平均最低溫度、平均最高溫度和總降水量3個變量。通過處理與合成生成多年平均最低溫度、多年平均最高溫度和多年總降水量數(shù)據(jù)。

    2.1.3 土壤數(shù)據(jù)

    土壤濕度數(shù)據(jù)來自全球土地數(shù)據(jù)同化系統(tǒng)(https://disc.gsfc.nasa.gov/),使用空間分辨率為0.25°的表層土壤濕度(0—10 cm)數(shù)據(jù),該數(shù)據(jù)時間范圍與氣象數(shù)據(jù)一致。采用表層土壤酸堿度數(shù)據(jù)(0—30 cm),獲取自中國科學(xué)院南京土壤研究所(http://www.soilinfo.cn/),其空間分辨率為1 km。

    2.1.4 植被數(shù)據(jù)

    使用MODIS NDVI月度歸一化植被指數(shù)(MOD13A3)表征研究區(qū)植被生長狀況,數(shù)據(jù)時間范圍與氣象數(shù)據(jù)保持一致,空間分辨率為1km,通過最大值合成方法,形成多年植被指數(shù)平均值。

    2.1.5 地形數(shù)據(jù)

    本研究使用的地形數(shù)據(jù)為數(shù)字高程模型(DEM),數(shù)據(jù)獲取自中國科學(xué)院地理科學(xué)與資源研究所資源環(huán)境科學(xué)與數(shù)據(jù)中心(http://www.resdc.cn/),數(shù)據(jù)空間分辨率為500 m,使用ArcGIS軟件生成研究區(qū)坡度、坡向數(shù)據(jù)。

    根據(jù)實際狀況,在開展數(shù)據(jù)分析和模型構(gòu)建前,使用ArcGIS軟件對所有環(huán)境因子數(shù)據(jù)重投影為WGS—84坐標,根據(jù)新疆行政區(qū)劃數(shù)據(jù)(由國家基礎(chǔ)地理信息中心發(fā)布的國家基礎(chǔ)地理信息數(shù)據(jù)庫提供)裁剪出研究區(qū)數(shù)據(jù),將所有數(shù)據(jù)重采樣至1 km,并將處理結(jié)果統(tǒng)一轉(zhuǎn)為Geo TIFF格式。

    2.2 研究方法

    為避免環(huán)境因子間共線性問題導(dǎo)致模型穩(wěn)健性降低,采用相關(guān)性分析方法和方差膨脹因子進行篩選。Green團隊推薦相關(guān)系數(shù)(r=0.7)為最佳閾值[37]。由于相關(guān)性分析有時存在非顯性的相關(guān)結(jié)構(gòu),學(xué)者[38]提出使用方差膨脹系數(shù)度量多重共線性影響的嚴重程度,一般認為當(dāng)方差膨脹系數(shù)高于10時,則模型會存在較為嚴重的多重共線性問題[39]。移除方差膨脹系數(shù)最大且相關(guān)性高于0.7的環(huán)境變量,使用其他變量構(gòu)建模型。

    物種分布模型在各領(lǐng)域獲得了很高的評價,如BIOCLIM、DOMAIN及最大熵(MaxEnt)等經(jīng)典模型[40]。但由于數(shù)據(jù)量的增大及各類模型應(yīng)用場景差異,通過集成建??山档湍P瓦x擇中的誤差,同時提供所有候選模型中各變量的相對重要性[41]。集成建模分簡單平均和加權(quán)平均兩種方法[20],其中加權(quán)平均基于單個模型的精度系數(shù)(AUC)分配權(quán)重[42],如式1:

    (1)

    式中,WApi為i的集成模型預(yù)測,AUCj為單一模型的AUC值;pij為各單一模型對i的預(yù)測值。

    本文使用R語言dismo(https://cran.r—project.org/web/packages/dismo/index.html)包實現(xiàn)各典型模型的預(yù)測,之后使用加權(quán)平均方法實現(xiàn)集成模型對新疆典型蝗蟲適生區(qū)預(yù)測。

    眾多模型評價方法中接受者操作特性曲線(ROC)曲線下面積(AUC)以及真實技巧統(tǒng)計(TSS)得到了廣泛的應(yīng)用,AUC將所有閾值范圍的值集成到一個單一最終指標中,而非在所有閾值上尋找度量的最大值[43];TSS是對Kappa系數(shù)優(yōu)點的繼承和對物種分布模型響應(yīng)曲線不足的克服,通過敏感性和特異性度量[44]。Araujo等[45]認為AUC低于0.5為反預(yù)測(即觀測值與預(yù)測值間存在負相關(guān)),AUC低于0.6為失敗,低于0.7為差,低于0.8為可接受,低于0.9為良好,高于0.9為優(yōu)秀;張雷等[46]認為TSS低于0.4為失敗,低于0.55為一般,低于0.70為好,低于0.85為很好,高于0.85為極好。

    3 結(jié)果與分析

    3.1 環(huán)境變量篩選

    通過相關(guān)分析得出變量間存在不同相關(guān)性,其中土壤濕度與降水變量間相關(guān)性最高(0.5027),而坡度與坡向間相關(guān)性最小(0.0021),相關(guān)性最高的變量未達到最佳閾值。將方差膨脹系數(shù)的相關(guān)性閾值設(shè)定為0.7,再次對所有變量進行篩選,得出月最低氣溫變量及地形變量與其他變量間共線性問題更為突出,如表1,因此予以剔除。在新疆典型蝗蟲適生區(qū)預(yù)測中使用除月最低氣溫數(shù)據(jù)及地形外的其余七個環(huán)境因子參與建模。

    表1 各環(huán)境變量方差膨脹系數(shù)

    3.2 模型精度評價

    將發(fā)生監(jiān)測點數(shù)據(jù)分為訓(xùn)練集和驗證集兩部分,其中驗證樣本集通過隨機選取20%樣點得到。利用AUC驗證各預(yù)測模型的性能得出均在可接受的范圍以內(nèi),如圖2和表2所示,其中DOMAIN模型的AUC值和TSS值均最低,分別為0.688和0.301;提升回歸樹模型最高,分別為0.920和0.910。

    圖2 各模型接受者操作特性曲線

    表2 各模型真實技巧統(tǒng)計值

    3.3 模型預(yù)測結(jié)果

    首先使用單一模型對新疆典型蝗蟲的適生區(qū)進行預(yù)測,得出高適宜性區(qū)域集中于天山以北區(qū)域,但預(yù)測結(jié)果存在一定差異,其中應(yīng)用BIOCLIM模型得出適生區(qū)面積最小,DOMAIN模型得出全疆大范圍內(nèi)均為高度適生區(qū),如圖3。地理分布上,高度適生區(qū)集中于北疆各大草原及天山山麓、阿爾泰山附近;行政區(qū)劃上,則集中于阿勒泰地區(qū)、塔城地區(qū)、博樂地區(qū)、石河子市、昌吉州、烏魯木齊市一帶,此外伊犁州、巴音郭楞蒙古自治州、克孜勒蘇柯爾克孜自治州及哈密部分區(qū)域也有分布。

    圖3 各模型對典型蝗蟲預(yù)測結(jié)果

    由于各模型性能存在差異,為提高預(yù)測模型可靠性和穩(wěn)健度,本研究將模型AUC和TSS值均更高的提升回歸樹、支持向量機和最大熵等3個高精度模型進行集成,以3個模型的加權(quán)平均最終預(yù)測新疆典型蝗蟲適生區(qū),如圖4。通過集成模型預(yù)測結(jié)果與單一模型相似,但通過集成模型預(yù)測的結(jié)果結(jié)合了三種模型的優(yōu)勢,綜合各模型在不同區(qū)域的預(yù)測結(jié)果,從而使得預(yù)測的區(qū)域達到更合理的范圍。在對三種典型蝗蟲統(tǒng)一預(yù)測的基礎(chǔ)上,運用集成模型分別針對三種典型蝗蟲進行適生區(qū)的預(yù)測,有利于認識到不同類型蝗蟲在全疆各地的地理分布特征和規(guī)律,以便在后期工作中對單一蝗蟲采取相對應(yīng)的監(jiān)測與防控措施。作為新疆典型蝗蟲適生區(qū)集成預(yù)測的重要基礎(chǔ),也分別對亞洲飛蝗、意大利蝗和西伯利亞蝗在全疆的適生區(qū)進行了模擬預(yù)測。三類蝗蟲的適生區(qū)具有較大比例的重合區(qū),從圖中看,蝗蟲的高度適生區(qū)基本都集中在新疆北部天山段及伊犁河谷等區(qū)域,對三類蝗蟲適生區(qū)的分別預(yù)測將會對蝗蟲的監(jiān)測、防治等各項工作的開展提供便利和支持。

    圖4 集成模型對不同蝗蟲類型預(yù)測結(jié)果

    使用ArcGIS自然間斷法,將集成模型預(yù)測的新疆典型蝗蟲適生區(qū)劃分為5個等級[24],如表3。根據(jù)劃分標準對集成模型預(yù)測結(jié)果分類,無論是新疆典型草原蝗蟲還是亞洲飛蝗、意大利蝗和西伯利亞蝗單一類型蝗蟲優(yōu)勢種而言,北疆各地州高適生區(qū)面積占比較大,南疆各地州則相對較少。而在北疆的地理分布則主要集中于天山一帶、阿爾泰山一帶、東部的吐哈盆地和西部的伊犁河谷一帶,南疆的適生區(qū)則呈帶狀分布于塔里木盆地南緣的綠洲一帶,如圖5。

    表3 新疆典型蝗蟲適生區(qū)劃分

    圖5 新疆典型蝗蟲適生區(qū)分級

    統(tǒng)計新疆典型蝗蟲不同適生區(qū)分區(qū)面積和比例,如表4,各類適生區(qū)(低適生區(qū)、邊緣適生區(qū)、適生區(qū)、高危區(qū))面積共約59.844萬km2,占新疆所有面積的36%。其中低適生區(qū)面積為28.284萬km2,所占比例最大,占所有適生區(qū)面積的近二分之一(47%),占所有總面積的17%;適生區(qū)面積1.735萬km2,所占比例最小,僅為不到百分之一(0.87%);邊緣適生區(qū)所占面積較低適生區(qū)低,但也高于其他適生區(qū)分區(qū);而高危區(qū)面積為16.568萬km2,所占面積僅次于低適生區(qū),所占比例也高達10%。

    表4 不同適生區(qū)分區(qū)面積及占比

    3.4 環(huán)境因子重要性

    新疆各典型蝗蟲對不同環(huán)境因子具有不同耐受性,本文綜合考慮影響意大利蝗、西伯利亞蝗和亞洲飛蝗生存與繁殖等活動的重要生態(tài)環(huán)境變量,各環(huán)境變量對預(yù)測新疆蝗蟲適生區(qū)分布的貢獻率如表5。綜合而言,在新疆三類典型蝗蟲適生區(qū)預(yù)測的建模中,不同環(huán)境變量在各主要預(yù)測模型中的貢獻率存在較大差異性,但是各個環(huán)境變量貢獻率的排名整體上較為一致。具體而言,從3個預(yù)測模型的表現(xiàn)上看,在對三類蝗蟲預(yù)測中,降水、歸一化植被指數(shù)及月最高氣溫3個變量始終是對蝗蟲適生區(qū)預(yù)測貢獻最大的幾類,其中在提升回歸樹模型中,對亞洲飛蝗、意大利蝗和西伯利亞蝗適生區(qū)預(yù)測貢獻率最高的變量均為降水,其表現(xiàn)分別為20.8%、53.7%和48.2%;支持向量機模型中,對亞洲飛蝗、意大利蝗和西伯利亞蝗適生區(qū)預(yù)測貢獻率最高的變量分別為土壤酸堿度、月最高氣溫和降水,其表現(xiàn)分別為20.13%、26.28%和23.55%;最大熵模型中,對亞洲飛蝗、意大利蝗和西伯利亞蝗適生區(qū)預(yù)測貢獻率最高的變量均為降水,其表現(xiàn)分別為31%、62.2%和67.3%。

    表5 環(huán)境變量貢獻率

    土壤濕度和土壤酸堿度及地形因子也同樣對新疆草原典型蝗蟲的適生區(qū)分布有較大影響,而坡向則在新疆典型蝗蟲發(fā)生和分布中的影響較小。環(huán)境變量在各預(yù)測模型分別針對三類典型蝗蟲適生區(qū)預(yù)測中的貢獻率存在一定差別,但根據(jù)3個模型的綜合表現(xiàn)月最高溫度、降水和植被3個變量始終是蝗蟲生長和繁育的重要且直接的環(huán)境因素,對于其生境的選擇具有很大的影響。

    4 討論與結(jié)論

    4.1 討論

    物種分布模型準確預(yù)測基礎(chǔ)之一是環(huán)境變量篩選,研究表明[47]環(huán)境變量會對模型預(yù)測精度造成較為顯著的影響。當(dāng)各環(huán)境變量間存在較高相關(guān)性時,會難以避免地帶來自相關(guān)和多重共線性等問題,從而在建模中造成不必要的信息冗余,影響模型的穩(wěn)健性。本研究采用相關(guān)性及方差膨脹系數(shù)移除月最低氣溫變量構(gòu)建預(yù)測模型,這與眾多學(xué)者在相關(guān)研究中采用的變量篩選方法基本一致,如學(xué)者Pradhan[48]使用皮爾遜相關(guān)系數(shù)對環(huán)境變量進行了相關(guān)性分析,其后計算方差膨脹系數(shù),對南亞和東南亞的生物氣候變量進行了多重共線性和冗余的篩選。

    對模型性能度量的重要標準之一是接受者操作特性曲線(即AUC值),AUC受閾值影響較低,可對所有閾值范圍內(nèi)的性能評價結(jié)果予以顯示,因此也受到生態(tài)環(huán)境領(lǐng)域?qū)W者重視[49]。AUC值越接近1即表明模型的性能越好,本研究中各模型的AUC基本都保持在0.7以上,說明各模型的性能可較好地預(yù)測新疆典型蝗蟲的適生區(qū),但模型間存在一定差異,其中DOMAIN性能相對最差,對新疆典型蝗蟲的預(yù)測能力最低;而提升回歸樹、支持向量機、最大熵3個模型的AUC均大于0.9,其性能在新疆典型蝗蟲的預(yù)測中高于其他各模型,這也與相關(guān)研究中各模型在不同物種分布預(yù)測模型中的性能表現(xiàn)基本一致[50]。具體如Duan等[51]使用BIOCLIM、DOMAIN、MAHAL、RF、MaxEnt和SVM等六個模型預(yù)測了五種樹種的潛在分布區(qū),以標準差、變異系數(shù)、Kappa值和AUC值評估模型的穩(wěn)定性,指出MAHAL、RF、MaxEnt和SVM具有更高預(yù)測精度。

    物種分布模型模擬出物種生存的可能性空間依賴于基本生態(tài)位的穩(wěn)定性及空間分布的平衡[52]。通過物種分布模型預(yù)測的物種潛在適生區(qū)與真實適生區(qū)之間存在一定差異,但預(yù)測的潛在適生區(qū)基本覆蓋了其實際的分布狀況[53]。本文參考陳健等[54]關(guān)于蝗蟲生境研究應(yīng)用進展中不同生境因子對蝗蟲影響的研究成果,結(jié)合新疆草原典型蝗蟲具體實際,綜合考慮了對蝗蟲有顯著影響的氣溫、地形、植被和土壤四個方面的因素,使預(yù)測的潛在適生區(qū)較真實蝗蟲適生區(qū)范圍最大化[55]。在新疆典型蝗蟲適生區(qū)預(yù)測中各生境因子對預(yù)測模型的貢獻率表現(xiàn)差異顯著,其中降水量對蝗蟲適生區(qū)的影響尤為顯著,這與王晗等[56]在新疆意大利蝗適生區(qū)與降水量、最高氣溫和最低氣溫的變化相關(guān)研究中相關(guān)的結(jié)論有一定的關(guān)聯(lián)性。因本文選取的環(huán)境變量并非在不同氣候情境下的生物氣候變量,預(yù)測變量的貢獻情況與李培先[24]等關(guān)于意大利蝗潛在分布中各變量貢獻率和楊會楓[23]關(guān)于西伯利亞蝗在新疆潛在分布中各變量貢獻率存在一定差異。

    本文對新疆各類典型蝗蟲的適生區(qū)預(yù)測結(jié)果顯示,典型蝗蟲的適生區(qū)具有顯著的區(qū)域差異,各類適生區(qū)面積近60萬km2,適生區(qū)地理分布主要集中在新疆天山一帶、北部阿爾泰山及伊犁河谷一帶,此外在東部的吐哈盆地和南疆塔里木盆地邊緣也有適生區(qū),這些基本與文獻有關(guān)新疆蝗蟲地理分布的表述基本一致。同時本文預(yù)測的適生區(qū)地理分布及其面積也與鄭淑丹對2012年新疆草地蝗蟲發(fā)生分布擬合結(jié)果一致,但本文預(yù)測結(jié)果無論是從預(yù)測方法上還是預(yù)測結(jié)果上均有了較大的提升和優(yōu)化;更近一步,與楊會楓等人使用MaxEnt模型對西伯利亞蝗潛在分布區(qū)的預(yù)測結(jié)果同樣具有相似性,這也印證了本文構(gòu)建的新疆典型蝗蟲預(yù)測模型及其預(yù)測結(jié)果對新疆典型蝗蟲適生區(qū)預(yù)測存在一定合理性和可靠性。

    本文基于新疆典型蝗蟲發(fā)生數(shù)據(jù)和相關(guān)環(huán)境變量,使用經(jīng)典物種分布模型對新疆典型蝗蟲適生區(qū)進行了預(yù)測,但還存在一些不足:(1)在環(huán)境變量選擇方面未全面考慮人類活動等非自然因素影響,因此可能對預(yù)測結(jié)果造成一定影響;(2)僅選用部分經(jīng)典物種分布模型進行了預(yù)測,暫未將如生態(tài)位因素分析模型(ENFA)、廣義加性模型(GAM)、人工神經(jīng)網(wǎng)絡(luò)(ANNs)等其他具有同樣代表性的模型考慮在內(nèi),可能也會對更精準預(yù)測結(jié)果的展示產(chǎn)生一定影響。(3)僅對新疆全區(qū)尺度典型蝗蟲適生區(qū)進行了預(yù)測,未對地州甚至區(qū)縣尺度狀況展開,另一方面也未對新疆周邊省區(qū)以及國外地區(qū)展開研究,在未來條件允許的情況下也可考慮在不同尺度下推進。總之,本研究對新疆典型蝗蟲適生區(qū)預(yù)測一方面于區(qū)域相關(guān)研究深入開展具有一定探索意義,另一方面于基層實際蝗蟲監(jiān)測、防治工作的推進亦可提供一定建議和支持。

    4.2 結(jié)論

    基于物種分布模型對新疆典型蝗蟲適生區(qū)預(yù)測,得出:

    (1)新疆典型蝗蟲適生區(qū)受多種環(huán)境因素共同作用,各因素相互影響對新疆蝗蟲適生區(qū)影響存在差異。本研究選取的眾多影響蝗蟲繁殖、生長等不同生命活動階段的環(huán)境因子中,月最低氣溫變量及地形相較于其他變量對蝗蟲適生區(qū)的影響最小。降水、氣溫和植被對物種分布預(yù)測模型的貢獻率最高,此外土壤濕度和土壤酸堿度對新疆典型蝗蟲適生區(qū)預(yù)測也有較高貢獻率。

    (2)各預(yù)測模型性能存在顯著差異,其中DOMAIN模型在新疆典型蝗蟲適生區(qū)預(yù)測中性能最低(AUC=0.688,TSS=0.301),而提升回歸樹模型的預(yù)測性能最佳(AUC=0.920,TSS=0.910);綜合預(yù)測性能最優(yōu)的3個經(jīng)典模型集成最終預(yù)測模型對新疆典型蝗蟲潛在適生區(qū)的預(yù)測具有更高的準確性與合理性。

    (3)通過集成預(yù)測模型的預(yù)測,新疆典型蝗蟲的潛在適生區(qū)分布具有區(qū)域性,其中天山以北的廣大草原地區(qū)存在各類蝗蟲適生區(qū),此外在東部的哈密地區(qū)、南疆各綠洲邊緣地帶亦有零星分布;各類適生區(qū)面積56.844萬 km2,其中高危區(qū)面積16.568萬 km2。

    (4)新疆草原典型蝗蟲適生區(qū)存在高度集中、區(qū)域差異顯著、地理分布結(jié)構(gòu)復(fù)雜等多種特點,不同區(qū)域、不同地州草原部門對蝗蟲的日常監(jiān)測和防治工作需緊密結(jié)合當(dāng)?shù)毓ぷ鲗嶋H,分區(qū)、分時制定嚴密的蝗蟲日常監(jiān)測和防治計劃,在做好本地蝗蟲監(jiān)測防治工作的同時,各邊境地區(qū)草原部門還需做好跨境蝗蟲的遷飛入境的監(jiān)測、預(yù)防工作。

    致謝:自治區(qū)級產(chǎn)學(xué)研聯(lián)合培養(yǎng)研究生基地——新疆維吾爾自治區(qū)治蝗滅鼠指揮部辦公室提供蝗蟲發(fā)生點數(shù)據(jù),國家基礎(chǔ)地理信息數(shù)據(jù)庫提供本文所有基礎(chǔ)底圖數(shù)據(jù)。

    猜你喜歡
    物種新疆模型
    一半模型
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    在新疆(四首)
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    回首2018,這些新物種值得關(guān)注
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    3D打印中的模型分割與打包
    瘋狂的外來入侵物種
    新疆多怪
    絲綢之路(2014年9期)2015-01-22 04:24:46
    在线观看美女被高潮喷水网站 | 99国产精品一区二区三区| 久久国产精品影院| 一级a爱片免费观看的视频| 亚洲国产欧美人成| 国模一区二区三区四区视频| 欧美国产日韩亚洲一区| 日韩欧美在线乱码| 欧美日韩黄片免| 亚洲美女视频黄频| 天天躁日日操中文字幕| 免费在线观看亚洲国产| 午夜免费男女啪啪视频观看 | 精品福利观看| 草草在线视频免费看| 综合色av麻豆| 欧美日韩福利视频一区二区| 啦啦啦观看免费观看视频高清| .国产精品久久| 村上凉子中文字幕在线| 国产精品爽爽va在线观看网站| 久9热在线精品视频| 乱码一卡2卡4卡精品| 人妻久久中文字幕网| 波多野结衣巨乳人妻| 亚洲一区二区三区色噜噜| 亚洲av.av天堂| 啪啪无遮挡十八禁网站| 欧美极品一区二区三区四区| 激情在线观看视频在线高清| www日本黄色视频网| 少妇裸体淫交视频免费看高清| 三级毛片av免费| 精品国内亚洲2022精品成人| 欧美三级亚洲精品| 免费高清视频大片| 精品一区二区三区视频在线| 国产精品电影一区二区三区| 91字幕亚洲| 有码 亚洲区| 欧美成狂野欧美在线观看| 直男gayav资源| 国产毛片a区久久久久| 日韩成人在线观看一区二区三区| 欧美激情在线99| 精品午夜福利视频在线观看一区| 欧美激情国产日韩精品一区| 中文字幕人成人乱码亚洲影| 亚洲黑人精品在线| 久久精品国产亚洲av天美| 黄色一级大片看看| 热99在线观看视频| 久久久久九九精品影院| 怎么达到女性高潮| 亚洲专区中文字幕在线| 99在线人妻在线中文字幕| 波多野结衣高清作品| 色噜噜av男人的天堂激情| 婷婷六月久久综合丁香| 久久精品国产清高在天天线| 亚洲 欧美 日韩 在线 免费| 亚洲va日本ⅴa欧美va伊人久久| 国产v大片淫在线免费观看| 老司机午夜十八禁免费视频| 精品久久久久久成人av| 日日夜夜操网爽| 久久久久免费精品人妻一区二区| 亚洲国产精品合色在线| 成人特级黄色片久久久久久久| 久久精品人妻少妇| 精华霜和精华液先用哪个| 日韩欧美一区二区三区在线观看| 亚洲自偷自拍三级| 最近中文字幕高清免费大全6 | 亚洲国产精品999在线| 欧美午夜高清在线| 九九热线精品视视频播放| 嫩草影院入口| xxxwww97欧美| 中文在线观看免费www的网站| 午夜免费男女啪啪视频观看 | 99久国产av精品| 好男人在线观看高清免费视频| 又爽又黄无遮挡网站| 99热这里只有精品一区| 免费看日本二区| 性欧美人与动物交配| 国产 一区 欧美 日韩| 91九色精品人成在线观看| 三级男女做爰猛烈吃奶摸视频| 人人妻,人人澡人人爽秒播| 一级a爱片免费观看的视频| 国产伦人伦偷精品视频| 精品99又大又爽又粗少妇毛片 | 久久久成人免费电影| 亚洲五月天丁香| 高清在线国产一区| 国产一区二区三区视频了| 亚洲熟妇熟女久久| 两人在一起打扑克的视频| 在线播放无遮挡| 亚洲成人精品中文字幕电影| 天天躁日日操中文字幕| 99热精品在线国产| 久久人人精品亚洲av| 免费观看人在逋| 校园春色视频在线观看| 午夜日韩欧美国产| 欧美最新免费一区二区三区 | 精品乱码久久久久久99久播| 亚洲av一区综合| 欧美色视频一区免费| 老司机福利观看| 9191精品国产免费久久| 国产一级毛片七仙女欲春2| 婷婷精品国产亚洲av在线| 中国美女看黄片| 人妻久久中文字幕网| 色综合婷婷激情| 可以在线观看毛片的网站| h日本视频在线播放| 99国产极品粉嫩在线观看| 一a级毛片在线观看| 五月玫瑰六月丁香| 精品国产亚洲在线| 亚洲国产欧洲综合997久久,| 久久午夜亚洲精品久久| 又爽又黄无遮挡网站| 精品久久久久久久人妻蜜臀av| 国产精品98久久久久久宅男小说| 精品一区二区三区视频在线观看免费| 久久精品91蜜桃| 男人舔女人下体高潮全视频| 一个人免费在线观看的高清视频| 久久中文看片网| av视频在线观看入口| 亚洲第一电影网av| 国产免费一级a男人的天堂| 九九热线精品视视频播放| 国产成人欧美在线观看| 亚洲 欧美 日韩 在线 免费| 国产精品女同一区二区软件 | 国产成人啪精品午夜网站| 国产熟女xx| 亚洲国产高清在线一区二区三| 色在线成人网| 久久精品久久久久久噜噜老黄 | 国产激情偷乱视频一区二区| 国产乱人伦免费视频| 欧美+亚洲+日韩+国产| 国产高清视频在线观看网站| 真人一进一出gif抽搐免费| 日韩欧美在线乱码| 国产精品一及| 免费在线观看成人毛片| 亚洲专区国产一区二区| 亚洲av电影在线进入| 精华霜和精华液先用哪个| 狂野欧美白嫩少妇大欣赏| 色尼玛亚洲综合影院| 国产在视频线在精品| 真人做人爱边吃奶动态| 悠悠久久av| 亚洲av成人av| 亚洲五月天丁香| 成年免费大片在线观看| 精品无人区乱码1区二区| 国内精品久久久久久久电影| 免费观看精品视频网站| 日韩欧美 国产精品| 久久精品影院6| 美女cb高潮喷水在线观看| 99热这里只有是精品在线观看 | 久久人人精品亚洲av| 亚洲国产欧美人成| 能在线免费观看的黄片| 日日夜夜操网爽| 国内久久婷婷六月综合欲色啪| 久久久久久久久中文| 国产黄片美女视频| 日韩中文字幕欧美一区二区| 色哟哟哟哟哟哟| 宅男免费午夜| 久久午夜福利片| 如何舔出高潮| 天天一区二区日本电影三级| 久久久久久九九精品二区国产| 亚洲 国产 在线| 色尼玛亚洲综合影院| 日韩精品中文字幕看吧| 亚洲国产精品999在线| 乱人视频在线观看| 久久午夜福利片| 俺也久久电影网| 久久人人精品亚洲av| 黄片小视频在线播放| 日韩av在线大香蕉| 婷婷丁香在线五月| 日韩国内少妇激情av| 精品人妻视频免费看| 精品免费久久久久久久清纯| 亚洲人成伊人成综合网2020| 国内精品久久久久精免费| 成人国产综合亚洲| 午夜影院日韩av| 九色成人免费人妻av| 青草久久国产| 午夜久久久久精精品| 在线观看免费视频日本深夜| 日本 av在线| 欧美成人一区二区免费高清观看| 免费看日本二区| 成人永久免费在线观看视频| 免费av不卡在线播放| 亚洲欧美日韩无卡精品| АⅤ资源中文在线天堂| 国产欧美日韩一区二区精品| 日韩大尺度精品在线看网址| 久久香蕉精品热| 国产精品免费一区二区三区在线| 免费观看的影片在线观看| 搡老熟女国产l中国老女人| 黄片小视频在线播放| 夜夜夜夜夜久久久久| 在线a可以看的网站| 亚洲狠狠婷婷综合久久图片| 亚洲三级黄色毛片| 99久久精品热视频| 国产成人欧美在线观看| 97碰自拍视频| 男人和女人高潮做爰伦理| 亚洲最大成人中文| 欧美性猛交╳xxx乱大交人| 精品人妻一区二区三区麻豆 | 男人狂女人下面高潮的视频| 动漫黄色视频在线观看| 亚洲国产精品sss在线观看| 亚洲精品粉嫩美女一区| 桃红色精品国产亚洲av| 国产精品亚洲av一区麻豆| 伦理电影大哥的女人| 国产高清激情床上av| 男女视频在线观看网站免费| 亚洲欧美精品综合久久99| 99久久99久久久精品蜜桃| 亚洲精品在线美女| 在线免费观看的www视频| 啪啪无遮挡十八禁网站| 中文字幕免费在线视频6| 网址你懂的国产日韩在线| 草草在线视频免费看| 欧美日本亚洲视频在线播放| 亚洲精品日韩av片在线观看| 一个人免费在线观看电影| 99热这里只有是精品在线观看 | 久久婷婷人人爽人人干人人爱| 午夜老司机福利剧场| 欧美日本视频| 18禁黄网站禁片免费观看直播| 一区二区三区四区激情视频 | 女人十人毛片免费观看3o分钟| 国产av麻豆久久久久久久| 99久久精品热视频| 中文字幕免费在线视频6| 国产精品久久久久久亚洲av鲁大| 欧美色视频一区免费| 国产白丝娇喘喷水9色精品| 露出奶头的视频| 最近最新中文字幕大全电影3| 日韩欧美国产一区二区入口| 深爱激情五月婷婷| 可以在线观看的亚洲视频| 看片在线看免费视频| 成年女人看的毛片在线观看| 国产av不卡久久| 国产美女午夜福利| 9191精品国产免费久久| 午夜福利18| 国产精品女同一区二区软件 | 夜夜爽天天搞| 可以在线观看的亚洲视频| 日韩欧美在线乱码| 亚洲精品日韩av片在线观看| 99久久无色码亚洲精品果冻| 高清在线国产一区| 两个人的视频大全免费| 免费搜索国产男女视频| 欧美激情久久久久久爽电影| 亚洲av一区综合| 天堂av国产一区二区熟女人妻| www.色视频.com| 蜜桃亚洲精品一区二区三区| 亚洲成人精品中文字幕电影| 国产视频内射| 日韩欧美国产一区二区入口| bbb黄色大片| 一区二区三区激情视频| a在线观看视频网站| 国内精品一区二区在线观看| 男女之事视频高清在线观看| 欧美+日韩+精品| 亚洲国产欧美人成| 国产伦精品一区二区三区四那| 91麻豆精品激情在线观看国产| 色av中文字幕| eeuss影院久久| 免费在线观看亚洲国产| 在线观看av片永久免费下载| 精品熟女少妇八av免费久了| 国产精品一区二区性色av| 日韩欧美精品免费久久 | 男人和女人高潮做爰伦理| av福利片在线观看| 丝袜美腿在线中文| 欧美不卡视频在线免费观看| 日日干狠狠操夜夜爽| 在线播放无遮挡| 女生性感内裤真人,穿戴方法视频| 九色国产91popny在线| 在线观看一区二区三区| 亚洲av电影不卡..在线观看| 18禁黄网站禁片免费观看直播| 在线看三级毛片| 99久久久亚洲精品蜜臀av| 男插女下体视频免费在线播放| 国产精品1区2区在线观看.| 天美传媒精品一区二区| 欧美乱色亚洲激情| 日本免费a在线| 亚洲午夜理论影院| 90打野战视频偷拍视频| 人人妻人人看人人澡| 免费电影在线观看免费观看| 国产视频内射| 亚洲精品在线美女| 亚洲av一区综合| 又爽又黄a免费视频| 国产高清视频在线播放一区| 国产高清有码在线观看视频| 欧美最黄视频在线播放免费| 日本 av在线| 国产淫片久久久久久久久 | 18禁黄网站禁片免费观看直播| 亚洲五月天丁香| 99久久精品一区二区三区| 丁香欧美五月| 能在线免费观看的黄片| 亚洲专区国产一区二区| 乱人视频在线观看| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| 国产一区二区亚洲精品在线观看| 少妇裸体淫交视频免费看高清| 亚洲最大成人手机在线| 99热精品在线国产| 中文字幕人成人乱码亚洲影| 国产淫片久久久久久久久 | 久久欧美精品欧美久久欧美| 夜夜爽天天搞| 国产在线男女| 宅男免费午夜| 麻豆av噜噜一区二区三区| h日本视频在线播放| 亚洲av免费在线观看| 欧美3d第一页| 深夜a级毛片| 国内精品久久久久久久电影| 日本与韩国留学比较| xxxwww97欧美| 国产午夜精品久久久久久一区二区三区 | 欧美日韩国产亚洲二区| or卡值多少钱| 久久久久久久久久成人| 两个人视频免费观看高清| 特级一级黄色大片| 久久6这里有精品| 757午夜福利合集在线观看| 最近中文字幕高清免费大全6 | 欧美激情久久久久久爽电影| 亚洲av电影在线进入| 精品久久久久久成人av| 亚洲第一区二区三区不卡| 亚洲精品粉嫩美女一区| 又爽又黄无遮挡网站| 午夜两性在线视频| 波多野结衣高清作品| 成年女人毛片免费观看观看9| 欧美不卡视频在线免费观看| 有码 亚洲区| 亚洲久久久久久中文字幕| 精品人妻偷拍中文字幕| 欧美激情在线99| 麻豆一二三区av精品| 国产成人福利小说| 少妇的逼水好多| 欧美一区二区国产精品久久精品| 真人做人爱边吃奶动态| 美女xxoo啪啪120秒动态图 | 色播亚洲综合网| 少妇高潮的动态图| 久久久久国内视频| 亚洲中文字幕一区二区三区有码在线看| 成年人黄色毛片网站| 麻豆久久精品国产亚洲av| 1024手机看黄色片| 欧美一级a爱片免费观看看| 亚洲av免费在线观看| 淫妇啪啪啪对白视频| 国产成年人精品一区二区| 国产成人a区在线观看| 长腿黑丝高跟| 99久久99久久久精品蜜桃| av在线天堂中文字幕| 国产亚洲精品av在线| 超碰av人人做人人爽久久| 欧美绝顶高潮抽搐喷水| 性插视频无遮挡在线免费观看| 三级男女做爰猛烈吃奶摸视频| www.www免费av| 成人特级黄色片久久久久久久| 成年人黄色毛片网站| 日韩欧美 国产精品| 无人区码免费观看不卡| 国产午夜精品久久久久久一区二区三区 | 日韩欧美国产一区二区入口| 禁无遮挡网站| 免费无遮挡裸体视频| 黄色视频,在线免费观看| 国产成人a区在线观看| 色在线成人网| 国产欧美日韩一区二区精品| 精品福利观看| 美女xxoo啪啪120秒动态图 | 亚洲av美国av| 校园春色视频在线观看| 国产av在哪里看| 日韩高清综合在线| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久国产蜜桃| 精品久久久久久久久久久久久| 亚洲美女黄片视频| 国产人妻一区二区三区在| 午夜两性在线视频| 一边摸一边抽搐一进一小说| 999久久久精品免费观看国产| av女优亚洲男人天堂| 国产三级中文精品| 特大巨黑吊av在线直播| 久久久成人免费电影| 国产大屁股一区二区在线视频| 亚洲av免费在线观看| 国产精品野战在线观看| 免费观看人在逋| 网址你懂的国产日韩在线| 俺也久久电影网| 国产高潮美女av| xxxwww97欧美| 午夜a级毛片| 久久人妻av系列| 精品久久久久久久久av| 欧美激情国产日韩精品一区| 最近最新中文字幕大全电影3| 看片在线看免费视频| 亚洲国产欧洲综合997久久,| 久久人妻av系列| 黄色女人牲交| 在线播放国产精品三级| www.色视频.com| 欧美高清性xxxxhd video| 午夜影院日韩av| 男人狂女人下面高潮的视频| 美女免费视频网站| 老司机午夜福利在线观看视频| 色在线成人网| 首页视频小说图片口味搜索| 久久国产乱子免费精品| 日本黄色片子视频| 少妇人妻一区二区三区视频| 国产成人欧美在线观看| 亚洲欧美日韩无卡精品| 久久草成人影院| 日韩中字成人| 人妻久久中文字幕网| 美女免费视频网站| 人妻久久中文字幕网| 国产不卡一卡二| 免费黄网站久久成人精品 | 久久久久九九精品影院| 亚洲欧美日韩东京热| 丰满乱子伦码专区| 国产探花在线观看一区二区| 日本精品一区二区三区蜜桃| 级片在线观看| 午夜a级毛片| 国产69精品久久久久777片| 久久欧美精品欧美久久欧美| 成人av在线播放网站| 日本黄大片高清| 日本一二三区视频观看| 国产精品亚洲一级av第二区| 一区二区三区四区激情视频 | 村上凉子中文字幕在线| 97超级碰碰碰精品色视频在线观看| 亚洲精品在线观看二区| 亚洲人成伊人成综合网2020| 亚洲乱码一区二区免费版| 99久国产av精品| avwww免费| 亚洲男人的天堂狠狠| www.色视频.com| 在线播放国产精品三级| 日韩欧美精品v在线| 一级黄片播放器| 亚洲av成人不卡在线观看播放网| 简卡轻食公司| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美日韩高清专用| 国产高潮美女av| 国产美女午夜福利| 91字幕亚洲| 亚洲人与动物交配视频| 久久久久亚洲av毛片大全| 国产单亲对白刺激| 亚洲国产精品sss在线观看| 色综合欧美亚洲国产小说| 亚洲综合色惰| 久久久精品大字幕| 麻豆国产97在线/欧美| 午夜日韩欧美国产| 最近视频中文字幕2019在线8| 中文资源天堂在线| 天美传媒精品一区二区| h日本视频在线播放| 久久午夜福利片| 亚洲 欧美 日韩 在线 免费| 欧美色视频一区免费| 日韩欧美国产一区二区入口| 桃色一区二区三区在线观看| 深夜a级毛片| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产清高在天天线| 欧美国产日韩亚洲一区| 一级a爱片免费观看的视频| 蜜桃亚洲精品一区二区三区| 男女下面进入的视频免费午夜| 亚洲第一欧美日韩一区二区三区| 一区福利在线观看| 内地一区二区视频在线| 国产三级在线视频| 女人被狂操c到高潮| 一a级毛片在线观看| 亚洲男人的天堂狠狠| 亚洲精品色激情综合| 欧美潮喷喷水| 亚洲av日韩精品久久久久久密| 国产成+人综合+亚洲专区| 国产黄色小视频在线观看| 黄色配什么色好看| 久久久久久久精品吃奶| 波野结衣二区三区在线| 日韩欧美精品免费久久 | 国产蜜桃级精品一区二区三区| 久久国产精品人妻蜜桃| 亚洲熟妇中文字幕五十中出| 亚洲三级黄色毛片| 国产av不卡久久| 免费一级毛片在线播放高清视频| 久久久久久国产a免费观看| 欧美日韩国产亚洲二区| 最近中文字幕高清免费大全6 | 免费观看精品视频网站| 欧美日韩中文字幕国产精品一区二区三区| netflix在线观看网站| 国模一区二区三区四区视频| av天堂中文字幕网| 成年人黄色毛片网站| 亚洲 国产 在线| 久久精品国产99精品国产亚洲性色| 久久精品夜夜夜夜夜久久蜜豆| 国内久久婷婷六月综合欲色啪| 哪里可以看免费的av片| 欧美一级a爱片免费观看看| 久久精品综合一区二区三区| 日韩欧美在线乱码| 亚洲七黄色美女视频| 国产白丝娇喘喷水9色精品| 久久性视频一级片| 午夜福利高清视频| 欧美3d第一页| 久久精品国产亚洲av香蕉五月| 中文字幕人妻熟人妻熟丝袜美| www.999成人在线观看| 99久久99久久久精品蜜桃| 成人午夜高清在线视频| 女同久久另类99精品国产91| 一本一本综合久久| 久久欧美精品欧美久久欧美| 91麻豆av在线| 国产欧美日韩一区二区三| 99热精品在线国产| 欧美一区二区亚洲| 一个人看的www免费观看视频| 美女cb高潮喷水在线观看| 午夜福利在线在线| 免费黄网站久久成人精品 | av黄色大香蕉| 亚洲熟妇熟女久久| 亚洲中文日韩欧美视频| 男女之事视频高清在线观看| 亚洲av二区三区四区| 欧美色视频一区免费| 自拍偷自拍亚洲精品老妇| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美三级三区| 9191精品国产免费久久|