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

    基于MaxEnt模型預測湖南省榧樹適宜分布區(qū)

    2024-09-13 00:00:00易烜楊振宇朱晉梅謝旭黃靚劉文劍王琢玙吳磊
    中南林業(yè)科技大學學報 2024年5期

    摘 要:【目的】榧樹Torreya grandis是我國特有的珍稀經濟樹種,預測榧樹在湖南省的適宜分布區(qū),探究影響其分布的主導環(huán)境因子,有助于其保護策略的制定和開發(fā)利用適宜地的選擇?!痉椒ā恳院鲜?93個野外調查榧樹分布點為基礎數(shù)據源,選取了22個影響其分布的環(huán)境變量,運用MaxEnt最大熵模型和ArcGIS 10.7軟件,分析影響其分布的主導環(huán)境因子,預測其在湖南省的適宜分布區(qū)?!窘Y果】1)ROC曲線下訓練數(shù)據集10次交叉驗證擬合的平均AUC值為0.914,測試數(shù)據集10次交叉驗證擬合的平均AUC值為0.867,標準差為0.074,模型預測精度可靠。2)限制榧樹當前潛在適生區(qū)分布的4個關鍵因素依次是海拔、坡度、最干季平均溫度以及最熱季降水量,其百分比貢獻率依次為29.40、16.90、8.60、7.30;置換重要性依次為53.30、14.60、5.70、6.50。3)榧樹在湖南省的適生區(qū)總面積為9.114×104 km2,高適生區(qū)、中適生區(qū)、低適生區(qū)面積分別為0.872×104、2.659×104、5.583×104 km2;榧樹適生區(qū)主要集中在湘中、西、東南地區(qū),高、中適生區(qū)集中分布在邵陽北部、南部、西南部,益陽西南部,永州中南部,婁底西北部以及懷化中南部?!窘Y論】地形、氣候等因子綜合影響著湖南省榧樹的地理分布格局,邵陽、益陽、永州、婁底、懷化等地區(qū)中,海拔較低、地勢相對平緩、氣候涼爽的溪谷地區(qū)可能適合榧樹的引種栽培。

    關鍵詞:最大熵MaxEnt模型;榧樹;適宜分布區(qū);影響因子

    中圖分類號:S718.54+5 文獻標志碼:A 文章編號:1673-923X(2024)05-0046-10

    基金項目:國家自然科學基金項目(31600160);中央財政林業(yè)科技推廣示范項目([2023]XTO7號);中央財政野生動植物保護項目(65235024);湖南省自然科學基金項目(2021JJ31152)。

    Identification of suitable distribution areas for Torreya grandis in Hunan province based on the MaxEnt model

    YI Xuan1, YANG Zhenyu2, ZHU Jinmei3, XIE Xu4, HUANG Liang1, LIU Wenjian1, WANG Zhuoyu1, WU Lei2

    (1. Qingyanghu State Owned Forest Farm, Ningxiang 410627, Hunan, China; 2. College of Forestry, Central South University of Forestry Technology, Changsha 410004, Hunan, China; 3. The Forestry Department of Hunan Province, Changsha 410004, Hunan, China; 4. Hunan Prospecting Designing Research Institute for Agriculture Forestry Industry, Changsha 410007, Hunan, China)

    Abstract:【Objective】Torreya grandis is a rare economic tree species endemic to China. Predicting the suitable distribution area of T. grandis in Hunan province and exploring the dominant environmental factors affecting the distribution of T. grandis will help to formulate the conservation strategy for the species and select the suitable sites for development and utilization.【Method】193 field survey points of T. grandis in Hunan province were taken, 22 environmental variables affecting the distribution of T. grandis were selected, and the MaxEnt model was used in conjunction with the ArcGIS 10.7 software to analyze the dominant environmental factors affecting the distribution of T. grandis and predict the suitable distribution areas of the species in Hunan province.【Result】1) the average AUC value of 10 cross-validated fits under the ROC curve for the training dataset was 0.914, and the average AUC value of 10 cross-validated fits for the test dataset was 0.867, with a standard deviation of 0.074, which made the model prediction results accurate; 2) Altitude, slope, mean temperature in the driest quarter, and precipitation of warmest quarter are the four key factors limiting T. grandis’s current potential suitable distribution, with percentage contributions of 29.40, 16.90, 8.60, and 7.30, respectively; and replacement importance of 53.30, 14.60, 5.70, and 6.50, respectively. 3) The total area of T. grandis suitable area in Hunan province was 9.114×104 km2, and the areas of high, medium and low suitable distribution areas were 0.872×104 km2, 2.659×104 km2 and 5.583×104 km2, respectively; the suitable distribution areas were mainly concentrated in the central, western and southeastern parts of Hunan province, and the high and medium suitable distribution areas were concentrated in the northern, southern and southwestern part of Shaoyang, southwestern Yiyang, central and southern Yongzhou, and the central and southern part of Yongzhou, and the southwestern part of Yiyang, and the southwestern part of Yongzhou. The suitable distribution areas were mainly concentrated in the central and southeastern areas of Hunan, with high and medium suitable distribution areas concentrated in the northern, southern and southwestern parts of Shaoyang, Yiyang, the central and southern parts of Yongzhou, the northwestern part of Loudi and the central and southern parts of Huaihua.【Conclusion】Topography, climate and other factors affect the geographical distribution pattern of Torreya grandis in Hunan Province. The valley areas of Shaoyang, Yiyang, Yongzhou, Loudi and Huaihua with low altitude, relatively flat terrain and cool climate may be suitable for the introduction and cultivation of Torreya grandis.

    Keywords: MaxEnt model; Torreya grandis; suitable distribution area; impact factor

    榧樹Torreya grandis,俗稱香榧,隸屬于紅豆杉科Taxaceae榧樹屬Torreya,是中國特有的珍稀經濟樹種,屬于國家二級保護野生植物。榧樹用途廣泛,除了觀賞與建筑,其果實也用途頗多,具有制作食用油、堅果和入藥等經濟價值,市場前景十分廣闊,適合在山區(qū)發(fā)展香榧產業(yè),能夠有效增加山區(qū)農民收入,促進當?shù)剞r林經濟發(fā)展,助力鄉(xiāng)村振興[1]。湖南是榧樹的主要分布區(qū)之一,現(xiàn)有的榧樹林多為天然林,如何合理布局其栽種區(qū)域,發(fā)揮其最大的經濟效益、生態(tài)效益,適地適樹是關鍵[2]。科學預測和評價榧樹的潛在適生區(qū),在生態(tài)經濟林建設中因地制宜,合理種植榧樹這一珍稀經濟樹種具有極其重要的理論與實踐意義。

    影響植物生長與地理分布的因素有很多,除了植物本身特性、人為干擾外,地形地貌、土壤因子、氣候因子等都是影響植物生長與地理分布的關鍵因素[3–7]。許多生態(tài)學與林學研究者通過構建生態(tài)位模型,分析物種分布點和影響其分布的環(huán)境因子之間的關系,推斷物種的生態(tài)需求,從而預測其潛在的分布區(qū)域[8-9]。如今,很多模型已經被用于預測物種分布,包括BRT、RF、BN和MaxEnt等模型[10–14]。其中,基于因子相似性原理的最大熵模型(maximum entropy modeling,MaxEnt)具有預測結果精度較高、適用于小樣本量、可操作性強(MaxEnt 軟件是免費的并且有友好的用戶界面)等特點,已被廣泛應用于許多物種適宜性生態(tài)分布區(qū)的擬合[4-5,15-16]。

    本研究以湖南省為研究區(qū),以珍稀經濟樹種榧樹為研究對象,構建MaxEnt模型,模擬湖南省榧樹潛在空間分布格局,劃分榧樹適生區(qū)等級,進一步確定影響湖南省榧樹分布的關鍵環(huán)境因子,為榧樹在湖南省種植與經營管理提供適地適樹的科學依據。

    1 材料與方法

    1.1 數(shù)據來源與處理

    除了氣候、土壤、地形、水文、人類活動和生物等限制條件外,擬合MaxEnt模型等物種分布模型時,是否對物種分布數(shù)據以及環(huán)境變量進行預處理對于物種生境適宜性評價結果的可靠性起著重要的作用[17-19]。物種地理分布坐標過于密集以及環(huán)境變量之間存在共線性等因素對模型產生了復雜的影響,從而降低了模型的適用性,并且導致了模擬結果的偏差[20]。本研究對193個物種分布數(shù)據、22個環(huán)境數(shù)據均做了相應的預處理。

    1.1.1 物種分布數(shù)據及預處理

    榧樹的實際分布數(shù)據來源于湖南省野外調查獲取的193個榧樹分布點(圖1a)。將榧樹分布數(shù)據整理為WGS84坐標系,并轉成便于模型計算的CSV格式??紤]到榧樹采樣分布點的空間自相關性[21],以30″×30″為一個樣本單元對原采樣點進行空間稀釋,采用ENMtools軟件進行自相關篩選后得到有效的榧樹分布點為41個(圖1b)。

    1.1.2 環(huán)境數(shù)據及預處理

    本研究選取19個生物氣候變量(bio1~bio19)以及3個地形數(shù)據(海拔Altitude、坡度Slope、坡向Aspect),這些數(shù)據源自世界氣候數(shù)據庫WorldClim(http: //www.worldclim.org/)網站(表1)。這22個環(huán)境變量圖層均為TIF格式(柵格數(shù)據),均為30″(約1 km)的空間分辨率,地理坐標系均為WGS84[22]。環(huán)境因子間往往存在一定的相關性,為避免MaxEnt模型過擬合[15],通過ArcGIS 10.7軟件將22個柵格格式環(huán)境變量轉換為ASCII格式;使用Maxent 3.4.4軟件將41個榧樹分布點與22個環(huán)境變量進行預擬合,計算環(huán)境因子的百分比貢獻率,取百分比累計貢獻率≥90%的環(huán)境變量;利用ENMtools軟件對22個環(huán)境因子進行相關性分析,對于相關系數(shù)|r|≥0.8的相關環(huán)境變量,只保留一個貢獻率高且有代表性的變量作為最后擬合模型的變量。根據相關分析、貢獻率的計算結果(圖2),最終保留8個環(huán)境變量用于MaxEnt模型構建,包括3個地形變量(海拔、坡度、坡向)和5個氣候變量(晝夜溫差與年溫差比值bio_3、最干季度平均溫度bio_9、降水量季節(jié)性變異性系數(shù)bio_15、最暖季度平均降水量bio_18、最冷季度平均降水量bio_19)。

    1.2 MaxEnt模型構建

    將榧樹的地理分布數(shù)據和環(huán)境變量數(shù)據導入Maxent 3.4.4軟件進行建模運算。采用Logistic方程揭示榧樹在湖南的潛在分布區(qū)。將75%的分布點作為訓練集,25%的分布點作為測試集。使用交叉驗證方法(Crossvalidate)處理物種分布數(shù)據,該方法將數(shù)據隨機平均分成10個子集,每次將其中1個子集作為測試集,其余9個子集作為訓練集,這樣重復運行10次,這種方式能夠充分利用數(shù)據。本研究設置的最大的迭代次數(shù)為5 000次,最大背景點數(shù)量為10 000個,其他參數(shù)采用默認設置。最終,將10次運行的結果取平均值,并將其保存為ASCII格式文件[23]。將結果文件導入到AreGIS10.7軟件中,采用自然間斷點分級法(Jenks’ natural breaks)將適宜度分為4個等級[24]:非適生區(qū)(0~0.08]、低適生區(qū)(0.08~0.26]、高適生區(qū)(0.26~0.51]和最適生區(qū)(0.51~1.00],從而得到香榧在湖南的潛在地理分布區(qū)[25],并通過堆積柱狀圖、三線表等將湖南省及各市級行政單位各個適生等級的面積進行分類統(tǒng)計。

    1.3 模型精度評價

    分類模型常使用接受者操作特征曲線(receiver operating characteristic curve,ROC)對模型擬合預測效果進行檢驗。在該曲線中,每個可能的判斷閾值都被視為預測結果的值,通過計算得到相應的靈敏度和特異度。模型預測準確度采用曲線下的面積值(Area under curve,AUC)來衡量。AUC的取值為0~1,值越大表示與隨機分布相距越遠,預測效果越好[26]。AUC值的優(yōu)點在于不受判斷閾值的影響,因而可以用來比較不同模型的預測準確度。根據AUC值的大小,可以對模型進行評估:0.7

    1.4 環(huán)境變量重要性評估

    MaxEnt模型中,百分比貢獻率、置換重要值、Jackknife檢驗3個指標或方法常用于評估環(huán)境因子對物種分布的重要性[9]。這3個指標或方法均能在Maxent 3.4.4軟件中實現(xiàn)。

    2 結果與分析

    2.1 模型精度評價

    ROC曲線下訓練數(shù)據集10次交叉驗證擬合的平均AUC值為0.914,測試數(shù)據集10次交叉驗證擬合的平均AUC值為0.867,標準差為0.074(圖3)?;谟柧?、測試數(shù)據集擬合AUC值均大于0.850,顯著高于隨機預測值(0.500)。預測結果很準確,這表示MaxEnt模型對榧樹在湖南省的潛在適生區(qū)的預測結果是可靠的。

    2.2 影響榧樹地理分布的主要環(huán)境變量和適宜性范圍

    由Jackknife(刀切法)檢驗結果(圖4)可知:地形因子中的坡向、氣候因子中的降水量季節(jié)性變異性系數(shù)(bio_15)的正規(guī)化訓練增益、測試增益均接近于0,AUC值均小于0.6,表明這兩個變量對當前湖南省榧樹的潛在適生區(qū)預測重要性較小。另外,等溫性、最冷季降水量(bio_19)的正規(guī)化訓練增益、測試增益表現(xiàn)不佳,均小于0.3。相反地,正規(guī)化訓練增益最大值為海拔,從大到小依次排序為海拔、坡度、最干季平均溫度、最熱季降水量;測試增益最大值也為海拔,最干季平均溫度次之,接著是坡度、最熱季降水量;AUC值最大的為最干季平均溫度,海拔次之,之后為最熱季降水量、坡度,這些環(huán)境變量展示了對數(shù)據的良好匹配,并且包含了更多有效的信息[18]。由表2可得,貢獻率、置換重要性前3名的環(huán)境變量依次為海拔、坡度、等溫性,累計貢獻率為62.5%,累計置換重要性為79.6%。綜上,當前影響榧樹在湖南省地理分布的主要氣候因子為最干季平均溫度、最熱季降水量;主要地形因子為海拔、坡度。

    圖5顯示,海拔、坡度、最干季平均溫度以及最熱季降水量是影響榧樹地理分布的主要環(huán)境因子,均具有一個峰值。在一定范圍內,榧樹的分布概率隨海拔、坡度以及最熱季降水量的升高而升高,到達一定峰值后,分布概率隨著環(huán)境因子的升高而降低。另外,榧樹分布概率先隨最干季平均溫度的升高而有所緩慢下降,再較快地升高,達到一定峰值后,呈較為迅速的下降趨勢。榧樹對4個主要環(huán)境變量具有一定的適應性,閾值分別為:19~1 864 m(海拔)、0~45.151°(坡度)、0.700~11.483 ℃(最干季平均溫度)、 381~698 mm(最熱季降水量);榧樹分布概率大于0.5的閾值分別為:397~884 m(海拔)、2.149°~19.758°(坡度)、6.687~8.253 ℃(最干季平均溫度)、448~500 mm(最熱季降水量)。

    2.3 榧樹適宜生長區(qū)的劃分

    湖南省榧樹適宜生長區(qū)劃分結果如圖6、7以及表3所示。適生區(qū)總面積為9.114×104 km2,占湖南省總面積的43.02%,其中邵陽市、張家界市、懷化市、婁底市以及湘西土家族苗族自治州的適生區(qū)面積均超過該行政區(qū)總面積的50%。高適生區(qū)主要集中在邵陽市北部和南部、永州市中部、益陽市西南部以及婁底市西北部,面積為0.872×104 km2。中適生區(qū)總面積主要集中在懷化市中南部、邵陽市西南部、益陽市西部、永州市中南部、婁底市西北部,面積為2.659×104 km2。低適生區(qū)主要集中在懷化市北部、湘西土家族苗族自治州和張家界市大部分地區(qū)、邵陽市東部、郴州市東南部,面積為5.583×104 km2。

    3 結論與討論

    3.1 討 論

    3.1.1 影響湖南省榧樹分布的主要環(huán)境變量

    研究區(qū)尺度不同,主導植被地理分布的環(huán)境變量也不盡相同。在大尺度(如全球范圍)下,水熱條件是決定植物分布的決定因素,而在較小尺度(景觀或更小尺度)下,區(qū)域環(huán)境因子是影響植被分布格局的主導因素[28]。本研究結果顯示,限制榧樹當前潛在適生區(qū)分布的4個最重要因素依次是海拔、坡度、最干季平均溫度、最熱季降水量。本研究區(qū)域范圍是湖南省,屬于上述研究中的小尺度區(qū)域。本研究區(qū)域海拔對榧樹的地理分布影響最大,這是由于植物種群的空間分布格局受海拔梯度的變化產生的溫度、濕度和光照等多種環(huán)境因子共同影響[29]。王義貴等[21]對我國特有樹種大果圓柏生境適宜性進行了評價,也得出了類似的結果,指出海拔是限制大果圓柏地理分布的最主要因子,認為海拔為主的地形因素協(xié)同降水、溫度、土壤因素綜合影響著大果圓柏的地理分布格局。另外,本研究發(fā)現(xiàn)湖南省榧樹分布概率大于0.5的閾值分別為:397~884 m(海拔)、2.149°~19.758°(坡度)、6.687~8.253 ℃(最干季平均溫度)、448~500 mm(最熱季降水量)。這與榧樹的生態(tài)習性一致:該種喜溫暖潮濕、光照較弱涼爽的氣候條件,最適生長于溪流兩旁和山腰谷地等朝夕多霧、直射光較少但散射光較多的區(qū)域。

    3.1.2 湖南省榧樹當前潛在適生區(qū)分析

    湖南省榧樹的適宜分布區(qū)主要集中在湘中、湘西、湘東南地區(qū),高、中適生區(qū)集中分布在邵陽(北部、南部、西南部)、益陽西南部、永州中南部、婁底西北部以及懷化中南部。在區(qū)域尺度下,地形地貌是主導植被空間分布格局的關鍵因素,降水和溫度等小氣候也受其影響[9]。湖南省的地形地貌呈朝東北開口的不對稱馬蹄形,東、南、西三面環(huán)山,中部丘陵起伏,北部平原展開[30]。湖南省榧樹分布格局與其地形地貌的分布密切相關,高、中適生區(qū)主要分布在湘中丘陵地區(qū),該地區(qū)平均海拔在500 m以下,溫度適宜且降水量充足,適合榧樹的生長。低適生區(qū)主要分布在湘西、湘西北地區(qū),這些地區(qū)主要由山地、山原組成,海拔多在1 000 m以上,冬季氣溫普遍偏低,低溫可能是制約榧樹分布的關鍵因子。而非適生區(qū)主要分布在洞庭湖平原沿湘、資、沅、澧等水系方向的區(qū)域,雖然這些地方地勢平坦,水分充足,但其日照較強,影響了喜陰的幼年榧樹正常生長;此外,人類活動也是影響植被分布的主要因素,其多屬于人類居住和活動相對集中的位置,人為干擾較為嚴重,不利于榧樹的存活[31]。

    MaxEnt模型只需要依據物種存在數(shù)據來構建,而缺少物種不存在數(shù)據來約束,很容易導致物種分布區(qū)的預測范圍偏大[32-33]。本研究榧樹的分布數(shù)據在多行政單位沒有出現(xiàn),例如張家界市、常德市等,模型擬合預測的榧樹分布范圍可能比實際范圍偏大。另外,土壤環(huán)境對植被的生長和演替有重要影響,不同的土壤理化性質(土壤養(yǎng)分、含水量等)作用于群落內的眾多生態(tài)過程,其中全磷、全鉀以及土壤有機質是制約天然林植物物種組成和分布的關鍵因子[34];相反,植被格局也影響著土壤養(yǎng)分、水分等輸送過程[35]。因此,在今后的研究中,還應考慮土壤理化性質對湖南榧樹分布與生長的影響。

    3.2 結 論

    基于湖南省榧樹193個野外調查的榧樹分布數(shù)據和19個生物氣候變量以及3個地形數(shù)據,采用MaxNnt物種分布模型,對湖南省榧樹的潛在適生區(qū)進行了可靠預測。總體而言,湖南省榧樹地理分布受地形因子影響比受氣候因子影響大,認為海拔、坡度是影響湖南省榧樹地理分布的主要地形因子,最干季平均溫度、最熱季降水量是影響湖南省榧樹地理分布的主要氣候因子。湖南省榧樹高、中適生區(qū)集中分布在邵陽(北部、南部、西南部)、益陽西南部、永州中南部、婁底西北部以及懷化中南部,這些地方可能適合湖南省榧樹的引種與種植。研究結果可為湖南省榧樹的合理管理和種植提供理論依據,并為今后湖南省榧樹森林資源的可持續(xù)發(fā)展提供理論參考。

    參考文獻:

    [1] 羅毅,張雪,佘新松,等.香榧林下復合種植模式對皖南山區(qū)土壤養(yǎng)分保持的影響[J].中國農學通報,2023,39(20):51-59. LUO Y, ZHANG X, SHE X S, et al. Effects of compounds plantations of Torreya grandis on soil nutrient conservation in mountainous areas of southern Anhui province[J]. Chinese Agricultural Science Bulletin,2023,39(20):51-59.

    [2] 李莎,史明昌,劉小宇,等.湖北省馬尾松精準適宜性空間分布[J].生態(tài)學報,2019,39(6):1960-1968. LI S, SHI M C, LIU X Y, et al. Precise spatial distribution of suitability of masson pine in Hubei province[J]. Acta Ecologica Sinica,2019,39(6):1960-1968.

    [3] 張曉瑋,蔣玉梅,畢陽,等.基于MaxEnt模型的中國沙棘潛在適宜分布區(qū)分析[J].生態(tài)學報,2022,42(4):1420-1428. ZHANG X W, JIANG Y M, BI Y, et al. Identification of potential distribution area for Hippophae rhamnoides subsp. sinensis by the MaxEnt model[J]. Acta Ecologica Sinica,2022,42(4):1420-1428.

    [4] 李燦,劉賢安,王娟,等.基于MaxEnt的四川省紅豆杉潛在分布區(qū)分析及適宜性評價[J].四川林業(yè)科技,2017,38(5): 1-7,32. LI C, LIU X A, WANG J, et al. Predictive distribution and habitat suitability assessment of Taxus chinensis based on MaxEnt in Sichuan province[J]. Journal of Sichuan Forestry Science and Technology,2017,38(5):1-7,32.

    [5] 付小勇,澤桑梓,周曉,等.基于MaxEnt的云南省薇甘菊分布預測及評價[J].廣東農業(yè)科學,2015,42(12):159-162. FU X Y, ZE S Z, ZHOU X, et al. Distribution prediction and assessment of Mikania micrantha in Yunnan province based on MaxEnt model[J]. Guangdong Agricultural Sciences,2015,42(12): 159-162.

    [6] 萬勝,孫佳,王磊,等.新疆不同產區(qū)棗果品質及其對生境因子的響應[J].經濟林研究,2023,41(1):97-105. WAN S, SUN J, WANG L, et al. Jujube fruit quality and its response to environment factors in Xinjiang different plantation areas[J]. Non-wood Forest Research,2023,41(1):97-105.

    [7] 陶紅,韓宏,馮廷敏,等.基于主成分分析的不同碧桃品種在甘肅地區(qū)的適應性評價[J].經濟林研究,2023,41(1):236-244. TAO H, HAN H, FENG T M, et al. Adaptability evaluation of Prunus persica cultivars in Gansu province based on principal component analysis[J]. Non-wood Forest Research,2023,41(1): 236-244.

    [8] 朱耿平,劉國卿,卜文俊,等.生態(tài)位模型的基本原理及其在生物多樣性保護中的應用[J].生物多樣性,2013,21(1):90-98. ZHU G P, LIU G Q, BU W J, et al. Ecological niche modeling and its applications in biodiversity conservation[J]. Biodiversity Science,2013,21(1):90-98.

    [9] 李璇,李垚,方炎明.基于優(yōu)化的Maxent模型預測白櫟在中國的潛在分布區(qū)[J].林業(yè)科學,2018,54(8):153-164. LI X, LI Y, FANG Y M. Prediction of potential suitable distribution areas of Quercus fabri in China based on an optimized Maxent model[J]. Scientia Silvae Sinicae,2018,54(8): 153-164.

    [10] LI X H, WANG Y. Applying various algorithms for species distribution modelling[J]. Integrative Zoology,2013,8(2): 124-135.

    [11] TANTIPISANUH N, GALE G A, POLLINO C. Bayesian networks for habitat suitability modeling: a potential tool for conservation planning with scarce resources[J]. Ecological Applications,2014,24(7):1705-1718.

    [12] 許仲林,彭煥華,彭守璋.物種分布模型的發(fā)展及評價方法[J].生態(tài)學報,2015,35(2):557-567. XU Z L, PENG H H, PENG S Z. The development and evaluation of species distribution models[J]. Acta Ecologica Sinica, 2015,35(2):557-567.

    [13] MI C R, HUETTMANN F, GUO Y M, et al. Why choose random forest to predict rare species distribution with few samples in large undersampled areas? Three Asian crane species models provide supporting evidence[J]. PeerJ,2017,5:e2849.

    [14] 郭彥龍,趙澤芳,喬慧捷,等.物種分布模型面臨的挑戰(zhàn)與發(fā)展趨勢[J].地球科學進展,2020,35(12):1292-1305. GUO Y L, ZHAO Z F, QIAO H J, et al. Challenges and development trend of species distribution model[J]. Advances in Earth Science,2020,35(12):1292-1305.

    [15] ZHOU Y C, ZHANG Z X, ZHU B, et al. MaxEnt modeling based on CMIP6 models to project potential suitable zones for Cunninghamia lanceolata in China[J]. Forests,2021,12(6):752.

    [16] HEBBAR K B, ABHIN P S, SANJO JOSE V, et al. Predicting the potential suitable climate for coconut (Cocos nucifera L.) cultivation in India under climate change scenarios using the MaxEnt model[J]. Plants,2022,11(6):731.

    [17] POUTEAU R, MEYER J-Y, LARRUE S. Using range filling rather than prevalence of invasive plant species for management prioritisation: the case of Spathodea campanulata in the Society Islands (South Pacific)[J]. Ecological Indicators,2015,54:87-95.

    [18] SILLERO N. What does ecological modelling model? A proposed classification of ecological niche models based on their underlying methods[J]. Ecological Modelling,2011,222(8): 1343-1346.

    [19] GUO Y L, LI X, ZHAO Z F, et al. Modeling the distribution of Populus euphratica in the Heihe River Basin, an inland river basin in an arid region of China[J]. Science China: Earth Sciences, 2018,61(11):1669-1684.

    [20] VERBRUGGEN H, TYBERGHEIN L, BELTON G S, et al. Improving transferability of introduced species’ distribution models: new tools to forecast the spread of a highly invasive seaweed[J]. PLoS ONE,2013,8(6):e68337.

    [21] 王義貴,何學高,胡云云,等.基于MaxEnt模型的大果圓柏生境適宜性評價[J].中南林業(yè)科技大學學報,2023,43(7):41-51. WANG Y G, HE X G, HU Y Y, et al. Prediction of potentially suitable distribution areas for Juniperus tibetica based on MaxEnt model[J]. Journal of Central South University of Forestry Technology,2023,43(7):41-51.

    [22] FICK S E, HIJMANS R J. WorldClim 2: new 1‐km spatial resolution climate surfaces for global land areas[J]. International Journal of Climatology,2017,37(12):4302-4315.

    [23] 李垚,張興旺,方炎明.小葉櫟分布格局對末次盛冰期以來氣候變化的響應[J].植物生態(tài)學報,2016,40(11):1164-1178. LI Y, ZHANG X W, FANG Y M. Responses of the distribution pattern of Quercus chenii to climate change following the Last Glacial Maximum[J]. Chinese Journal of Plant Ecology, 2016,40(11):1164-1178.

    [24] 寧瑤,雷金睿,宋希強,等.石灰?guī)r特有植物海南鳳仙花潛在適宜生境分布模擬[J].植物生態(tài)學報,2018,42(9):946-954. NING Y, LEI J R, SONG X Q, et al. Modeling the potential suitable habitat of Impatiens hainanensis, a limestone-endemic plant[J]. Chinese Journal of Plant Ecology,2018,42(9):946-954.

    [25] 張華,趙浩翔,王浩.基于Maxent模型的未來氣候變化情景下胡楊在中國的潛在地理分布[J].生態(tài)學報,2020,40(18): 6552-6563. ZHANG H, ZHAO H X, WANG H. Potential geographical distribution of Populus euphratica in China under future climate change scenarios based on Maxent model[J]. Acta Ecologica Sinica, 2020,40(18):6552-6563.

    [26] SHEKHAR C, GINWAL H S, MEENA R K, et al. Spatiotemporal distribution of broad-leaved Quercus semecarpifolia indicates altitudinal shift in northwestern Himalayas[J]. Plant Ecology,2022,223(6):671-697.

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

    [28] BURKE A. Classification and ordination of plant communities of the Naukluft Mountains, Namibia[J]. Journal of Vegetation Science,2001,12(1):53-60.

    [29] 唐志紅,尉秋實,劉虎俊,等.祁連山東段高寒植被群落特征及其與地形氣候因子關系研究[J].生態(tài)學報,2020,40(1): 223-232. TANG Z H, YU Q S, LIU H J, et al. Characteristics of alpine vegetation community and its relationship to topographic climate factors in the eastern Qilian Mountain[J]. Acta Ecologica Sinica, 2020,40(1):223-232.

    [30] 黃國林,曾斌,李衛(wèi)東,等.基于地貌形態(tài)的湖南星級農莊空間分布[J].江蘇農業(yè)科學,2019 47(21):28-32. HUANG G L, ZENG B, LI W D, et al. Spatial distribution of star farms in Hunan based on geomorphological morphology[J]. Jiangsu Agricultural Sciences,2019,47(21):28-32.

    [31] 周炳江,王玉潔,馬長樂,等.基于MaxEnt與ArcGIS的云南榧樹潛在生境分析[J].生態(tài)學報,2022,42(11):4485-4493. ZHOU B J, WANG Y J, MA C L, et al. Analysis of potential habitat of Torreya yunnanensis based on MaxEnt and ArcGIS[J]. Acta Ecologica Sinica,2022,42(11):4485-4493.

    [32] CHEFAOUI R M, LOBO J M. Assessing the effects of pseudoabsences on predictive distribution model performance[J]. Ecological Modelling,2008,210(4):478-486.

    [33] 陳禹光,樂新貴,陳宇涵,等.基于MaxEnt模型預測氣候變化下杉木在中國的潛在地理分布[J].應用生態(tài)學報,2022,33(5): 1207-1214. CHEN Y G, LE X G. CHEN Y H, et al. Identification of the potential distribution area of Cunninghamia lanceolata in China under climate change based on the MaxEnt model[J]. Chinese Journal of Applied Ecology,2022,33(5):1207-1214.

    [34] 徐明鋒,胡硯秋,李文斌,等.土壤養(yǎng)分對亞熱帶天然林物種分布的影響[J].中南林業(yè)科技大學學報,2014,34(9):91-97. XU M F, HU Y Q, LI W B, et al. Effects of soil nutrients on species distribution in a subtropical natural forest[J]. Journal of Central South University of Forestry Technology,2014,34(9): 91-97.

    [35] 馬宗文,謝正磊,段曉峰,等.黃河三角洲自然保護區(qū)植物與土壤因子關系及生態(tài)位分析[J].北京大學學報(自然科學版), 2012,48(5):801-811. MA Z W, XIE Z L, DUAN X F, et al. Plant-soil relationship and plant niche in the Yellow River Delta National Natural Reserve, China[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2012,48(5):801-811.

    [本文編校:謝榮秀]

    欧美人与善性xxx| 国产一区二区在线观看日韩| 亚洲美女黄色视频免费看| 日产精品乱码卡一卡2卡三| 亚洲美女黄色视频免费看| 亚洲怡红院男人天堂| 国产精品一及| 国产成人精品婷婷| 亚洲精品日本国产第一区| 最近手机中文字幕大全| 亚洲精品日本国产第一区| 99热6这里只有精品| 欧美成人精品欧美一级黄| 18禁动态无遮挡网站| 亚洲精品乱码久久久v下载方式| 亚洲精品国产av蜜桃| 久久久久久久国产电影| 老女人水多毛片| 一级毛片aaaaaa免费看小| 色5月婷婷丁香| 国产一区二区三区综合在线观看 | 久久女婷五月综合色啪小说| 男女免费视频国产| 欧美精品一区二区免费开放| 少妇被粗大猛烈的视频| 最近最新中文字幕免费大全7| 色婷婷av一区二区三区视频| 精品久久国产蜜桃| 午夜免费鲁丝| 国模一区二区三区四区视频| 成年人午夜在线观看视频| 久久久久性生活片| 中文字幕制服av| 国产大屁股一区二区在线视频| 久久精品熟女亚洲av麻豆精品| 亚洲性久久影院| 少妇高潮的动态图| 国产免费福利视频在线观看| 亚洲av国产av综合av卡| 婷婷色综合大香蕉| 99热网站在线观看| 国产亚洲午夜精品一区二区久久| 一级二级三级毛片免费看| 少妇猛男粗大的猛烈进出视频| 精品久久久精品久久久| 日韩 亚洲 欧美在线| 在线 av 中文字幕| 亚洲成人一二三区av| av黄色大香蕉| 久久精品熟女亚洲av麻豆精品| 日韩欧美一区视频在线观看 | 直男gayav资源| 又爽又黄a免费视频| a级一级毛片免费在线观看| 中文资源天堂在线| 国产av精品麻豆| 久久热精品热| 精品国产乱码久久久久久小说| 五月玫瑰六月丁香| 一级毛片久久久久久久久女| 一级毛片我不卡| 久久精品国产亚洲网站| 免费看av在线观看网站| 国产亚洲av片在线观看秒播厂| 国产av一区二区精品久久 | 精华霜和精华液先用哪个| 日日撸夜夜添| 久久久精品94久久精品| 国产av国产精品国产| 老司机影院成人| 极品教师在线视频| 联通29元200g的流量卡| 亚洲欧美日韩东京热| 欧美日韩一区二区视频在线观看视频在线| 好男人视频免费观看在线| 亚洲aⅴ乱码一区二区在线播放| 国产精品熟女久久久久浪| 国产精品久久久久久精品古装| av在线老鸭窝| av专区在线播放| av天堂中文字幕网| 婷婷色麻豆天堂久久| www.色视频.com| 在线观看国产h片| 日韩欧美精品免费久久| 麻豆乱淫一区二区| 又爽又黄a免费视频| 亚洲精华国产精华液的使用体验| .国产精品久久| 亚洲精品日韩在线中文字幕| 九九久久精品国产亚洲av麻豆| 噜噜噜噜噜久久久久久91| 欧美日韩国产mv在线观看视频 | 婷婷色麻豆天堂久久| 最近最新中文字幕免费大全7| 久久影院123| 欧美成人一区二区免费高清观看| 免费观看的影片在线观看| 亚洲va在线va天堂va国产| 亚洲,一卡二卡三卡| 91精品国产国语对白视频| 国产综合精华液| 热99国产精品久久久久久7| 国产中年淑女户外野战色| 老熟女久久久| 蜜桃在线观看..| 在线天堂最新版资源| 少妇人妻一区二区三区视频| 亚洲精品色激情综合| 哪个播放器可以免费观看大片| 日韩不卡一区二区三区视频在线| 欧美国产精品一级二级三级 | 免费观看的影片在线观看| 亚洲精品中文字幕在线视频 | 亚洲欧美中文字幕日韩二区| 精品99又大又爽又粗少妇毛片| 免费看不卡的av| 免费观看av网站的网址| 成人黄色视频免费在线看| www.色视频.com| 黄色怎么调成土黄色| av女优亚洲男人天堂| 亚洲美女搞黄在线观看| 中文资源天堂在线| 久久人妻熟女aⅴ| 最新中文字幕久久久久| 国产视频首页在线观看| av不卡在线播放| 在线亚洲精品国产二区图片欧美 | 插阴视频在线观看视频| 亚洲精品456在线播放app| 国产伦精品一区二区三区四那| av天堂中文字幕网| 精华霜和精华液先用哪个| 成人特级av手机在线观看| 色婷婷久久久亚洲欧美| 国产精品蜜桃在线观看| h视频一区二区三区| 精品国产一区二区三区久久久樱花 | 女人久久www免费人成看片| 亚洲真实伦在线观看| 一级片'在线观看视频| 国语对白做爰xxxⅹ性视频网站| 韩国高清视频一区二区三区| 国产日韩欧美在线精品| 国产成人精品福利久久| 自拍偷自拍亚洲精品老妇| 午夜激情久久久久久久| 一区二区三区免费毛片| 香蕉精品网在线| 国产午夜精品一二区理论片| 成人毛片60女人毛片免费| 美女国产视频在线观看| 国产av国产精品国产| 久久精品国产亚洲网站| 一级毛片aaaaaa免费看小| 亚洲电影在线观看av| 久久精品久久久久久噜噜老黄| 成人漫画全彩无遮挡| 国产精品嫩草影院av在线观看| 亚洲精品乱码久久久v下载方式| 亚洲久久久国产精品| 91午夜精品亚洲一区二区三区| 在线观看三级黄色| 亚洲高清免费不卡视频| 亚洲精品456在线播放app| 在线播放无遮挡| 久久热精品热| 日本猛色少妇xxxxx猛交久久| 亚洲欧美清纯卡通| 又爽又黄a免费视频| 国产精品一区二区三区四区免费观看| freevideosex欧美| 国产精品三级大全| 男人爽女人下面视频在线观看| 精品亚洲成国产av| av福利片在线观看| 美女主播在线视频| 成人18禁高潮啪啪吃奶动态图 | 国产熟女欧美一区二区| 26uuu在线亚洲综合色| 日韩成人av中文字幕在线观看| 人妻系列 视频| 成人午夜精彩视频在线观看| 纵有疾风起免费观看全集完整版| 日韩三级伦理在线观看| 久久午夜福利片| 晚上一个人看的免费电影| 欧美人与善性xxx| 街头女战士在线观看网站| 国产视频首页在线观看| 一级黄片播放器| 蜜桃在线观看..| 这个男人来自地球电影免费观看 | 丝袜喷水一区| 欧美区成人在线视频| 五月伊人婷婷丁香| 久久精品久久精品一区二区三区| 妹子高潮喷水视频| 亚洲,欧美,日韩| 国产精品福利在线免费观看| 欧美一区二区亚洲| 色视频www国产| 人妻一区二区av| 国产精品成人在线| 久久久久人妻精品一区果冻| 国产中年淑女户外野战色| 精品视频人人做人人爽| av在线蜜桃| 国产亚洲欧美精品永久| 国内精品宾馆在线| 久久国产精品男人的天堂亚洲 | 九色成人免费人妻av| 成年人午夜在线观看视频| 美女中出高潮动态图| 91久久精品电影网| 多毛熟女@视频| 亚洲人与动物交配视频| 欧美变态另类bdsm刘玥| 交换朋友夫妻互换小说| 亚洲av福利一区| 在线 av 中文字幕| 精品久久久久久久久av| 一本一本综合久久| 黄色欧美视频在线观看| 边亲边吃奶的免费视频| 一边亲一边摸免费视频| 色婷婷久久久亚洲欧美| 亚洲av中文字字幕乱码综合| 亚洲精品亚洲一区二区| 夫妻午夜视频| 国产精品久久久久久精品电影小说 | 亚洲性久久影院| 国产精品福利在线免费观看| 亚洲国产欧美人成| 久久久久精品性色| 久久婷婷青草| 日日啪夜夜撸| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品古装| 不卡视频在线观看欧美| 日韩中字成人| 毛片一级片免费看久久久久| 直男gayav资源| 欧美三级亚洲精品| 日韩av免费高清视频| 国产精品久久久久久精品电影小说 | 1000部很黄的大片| 有码 亚洲区| 少妇熟女欧美另类| 久久精品国产亚洲av涩爱| av免费观看日本| 91久久精品国产一区二区成人| 国产久久久一区二区三区| 人人妻人人澡人人爽人人夜夜| 我的女老师完整版在线观看| 五月天丁香电影| 99久久精品一区二区三区| freevideosex欧美| 天美传媒精品一区二区| 嘟嘟电影网在线观看| 亚洲最大成人中文| 久久精品国产自在天天线| 久久韩国三级中文字幕| 新久久久久国产一级毛片| 免费观看av网站的网址| 亚洲国产欧美在线一区| 国产精品一区www在线观看| 欧美高清成人免费视频www| 街头女战士在线观看网站| 免费观看无遮挡的男女| 久久精品国产a三级三级三级| 日本免费在线观看一区| 国产亚洲欧美精品永久| 亚洲精品成人av观看孕妇| 色5月婷婷丁香| 精品国产三级普通话版| 日韩免费高清中文字幕av| 午夜免费观看性视频| 亚洲欧美精品专区久久| 毛片一级片免费看久久久久| 99热这里只有是精品50| 亚洲精品乱码久久久久久按摩| 国产综合精华液| 久久精品国产亚洲av天美| 你懂的网址亚洲精品在线观看| 日韩av免费高清视频| 色视频www国产| 人人妻人人澡人人爽人人夜夜| xxx大片免费视频| a级毛色黄片| 看十八女毛片水多多多| 久久久久视频综合| 最近的中文字幕免费完整| 亚洲在久久综合| 七月丁香在线播放| 亚洲人与动物交配视频| 日韩国内少妇激情av| 一级毛片 在线播放| 久久精品夜色国产| 黄色一级大片看看| 十分钟在线观看高清视频www | 99久久综合免费| 赤兔流量卡办理| av国产免费在线观看| 精品久久久精品久久久| 亚洲av中文字字幕乱码综合| 日本午夜av视频| 丰满少妇做爰视频| 人妻 亚洲 视频| 亚洲精品日韩在线中文字幕| 午夜免费鲁丝| 国产 一区 欧美 日韩| 欧美精品人与动牲交sv欧美| 亚洲国产精品一区三区| 亚洲图色成人| 纵有疾风起免费观看全集完整版| 国产高清国产精品国产三级 | 最新中文字幕久久久久| av不卡在线播放| 中国美白少妇内射xxxbb| 亚洲欧美成人综合另类久久久| 亚洲欧美清纯卡通| 久久久欧美国产精品| 少妇丰满av| 亚洲av免费高清在线观看| 国产精品国产av在线观看| 观看av在线不卡| 亚洲婷婷狠狠爱综合网| 联通29元200g的流量卡| 国产 一区精品| 性色avwww在线观看| 午夜视频国产福利| 亚洲美女视频黄频| 在线免费十八禁| 伦精品一区二区三区| 十分钟在线观看高清视频www | 最近手机中文字幕大全| 97超碰精品成人国产| 亚洲欧洲国产日韩| av女优亚洲男人天堂| 国产在线视频一区二区| 大码成人一级视频| 日本黄色片子视频| 欧美日韩综合久久久久久| 多毛熟女@视频| 成年av动漫网址| 交换朋友夫妻互换小说| 国产高清有码在线观看视频| 久久久久久久久久久丰满| 午夜福利高清视频| 一区二区三区四区激情视频| 亚洲人成网站高清观看| 街头女战士在线观看网站| 最近手机中文字幕大全| 久久毛片免费看一区二区三区| 少妇人妻 视频| 国产黄色视频一区二区在线观看| 国产黄片美女视频| 亚洲国产欧美人成| 亚洲不卡免费看| 久热久热在线精品观看| 国产av精品麻豆| 国产精品一区二区在线不卡| 亚洲图色成人| av在线观看视频网站免费| 一区二区三区乱码不卡18| 亚洲精华国产精华液的使用体验| 一区二区av电影网| 免费黄色在线免费观看| 国产精品99久久久久久久久| 久久久成人免费电影| 亚洲人与动物交配视频| 免费av不卡在线播放| 欧美少妇被猛烈插入视频| 老熟女久久久| 欧美亚洲 丝袜 人妻 在线| 王馨瑶露胸无遮挡在线观看| 少妇的逼水好多| 成人国产av品久久久| 在线亚洲精品国产二区图片欧美 | 精品久久久噜噜| 国产一区二区在线观看日韩| 狂野欧美白嫩少妇大欣赏| 99精国产麻豆久久婷婷| 黄色怎么调成土黄色| 成人综合一区亚洲| 免费播放大片免费观看视频在线观看| 熟女人妻精品中文字幕| 国产黄片美女视频| 国产精品女同一区二区软件| 黄色欧美视频在线观看| 搡女人真爽免费视频火全软件| 少妇熟女欧美另类| 成年美女黄网站色视频大全免费 | 汤姆久久久久久久影院中文字幕| 97热精品久久久久久| 亚洲精华国产精华液的使用体验| 最近最新中文字幕免费大全7| 国产v大片淫在线免费观看| 久久精品人妻少妇| 能在线免费看毛片的网站| 丰满人妻一区二区三区视频av| 老司机影院成人| 国产精品一区二区在线观看99| 18禁在线播放成人免费| 亚洲第一av免费看| 搡女人真爽免费视频火全软件| 建设人人有责人人尽责人人享有的 | 亚洲欧美一区二区三区黑人 | 国产精品免费大片| 又爽又黄a免费视频| 久久国产乱子免费精品| 能在线免费看毛片的网站| 国产成人一区二区在线| 九色成人免费人妻av| 精品亚洲成a人片在线观看 | 性色av一级| 欧美精品人与动牲交sv欧美| 好男人视频免费观看在线| kizo精华| 18禁裸乳无遮挡动漫免费视频| 欧美3d第一页| 亚洲av欧美aⅴ国产| 97热精品久久久久久| 亚洲最大成人中文| 日韩强制内射视频| 国产免费又黄又爽又色| 亚洲av.av天堂| 国产无遮挡羞羞视频在线观看| 五月伊人婷婷丁香| 亚洲精品一区蜜桃| 国产国拍精品亚洲av在线观看| 久久这里有精品视频免费| 丰满迷人的少妇在线观看| 看免费成人av毛片| 一区二区三区乱码不卡18| 三级国产精品欧美在线观看| 韩国av在线不卡| 亚洲精品日韩av片在线观看| 亚洲精品456在线播放app| 亚洲真实伦在线观看| 乱系列少妇在线播放| 欧美日韩综合久久久久久| 噜噜噜噜噜久久久久久91| 欧美高清性xxxxhd video| 国产黄频视频在线观看| 能在线免费看毛片的网站| 少妇高潮的动态图| 视频中文字幕在线观看| 成人午夜精彩视频在线观看| 久久久久视频综合| av在线app专区| 视频区图区小说| 精品99又大又爽又粗少妇毛片| 我的女老师完整版在线观看| av视频免费观看在线观看| 午夜免费鲁丝| 国产有黄有色有爽视频| 国产中年淑女户外野战色| av国产免费在线观看| av黄色大香蕉| 黑人猛操日本美女一级片| 18禁在线无遮挡免费观看视频| 七月丁香在线播放| 建设人人有责人人尽责人人享有的 | 国产成人精品福利久久| 色视频在线一区二区三区| 大又大粗又爽又黄少妇毛片口| 久久久久久久久大av| 九九爱精品视频在线观看| 日韩人妻高清精品专区| 在现免费观看毛片| 制服丝袜香蕉在线| 2021少妇久久久久久久久久久| 99精国产麻豆久久婷婷| 久久韩国三级中文字幕| 国产成人午夜福利电影在线观看| 国产精品三级大全| 亚洲第一av免费看| 国产男人的电影天堂91| 老师上课跳d突然被开到最大视频| 日本av免费视频播放| 一本一本综合久久| 老司机影院毛片| 日本爱情动作片www.在线观看| 亚洲无线观看免费| 国产成人精品婷婷| 老师上课跳d突然被开到最大视频| 亚洲伊人久久精品综合| 极品教师在线视频| 国产午夜精品久久久久久一区二区三区| 国产黄色视频一区二区在线观看| 交换朋友夫妻互换小说| av不卡在线播放| 日本-黄色视频高清免费观看| 一本久久精品| 欧美日韩综合久久久久久| 亚洲婷婷狠狠爱综合网| 国产精品久久久久久久电影| 麻豆精品久久久久久蜜桃| 日韩精品有码人妻一区| 久久精品国产自在天天线| 成人美女网站在线观看视频| 天堂8中文在线网| 99精国产麻豆久久婷婷| 久久精品国产亚洲网站| 亚洲综合色惰| 国产精品不卡视频一区二区| 天天躁日日操中文字幕| 交换朋友夫妻互换小说| 国产欧美日韩一区二区三区在线 | 国产精品福利在线免费观看| 老女人水多毛片| 成人18禁高潮啪啪吃奶动态图 | 国产大屁股一区二区在线视频| a级毛片免费高清观看在线播放| www.av在线官网国产| 老司机影院成人| 日本色播在线视频| 丝袜喷水一区| 欧美日韩视频高清一区二区三区二| 国产av码专区亚洲av| 三级国产精品片| 亚洲av.av天堂| 亚洲精品一二三| 天天躁日日操中文字幕| 亚洲精品亚洲一区二区| 久久久久久久久久久丰满| 亚洲精品国产成人久久av| 性高湖久久久久久久久免费观看| 水蜜桃什么品种好| 日韩不卡一区二区三区视频在线| 中文字幕免费在线视频6| 99久久精品热视频| 一级爰片在线观看| 久久精品人妻少妇| 欧美精品人与动牲交sv欧美| 国产精品久久久久久av不卡| 日日摸夜夜添夜夜添av毛片| 国产精品av视频在线免费观看| 亚洲精品456在线播放app| 精品人妻视频免费看| 免费看日本二区| 少妇被粗大猛烈的视频| 亚洲不卡免费看| 最近最新中文字幕免费大全7| 亚洲av成人精品一区久久| 国产亚洲欧美精品永久| a级一级毛片免费在线观看| 有码 亚洲区| 久久久久国产网址| 亚洲欧美中文字幕日韩二区| 一区二区av电影网| 国产精品不卡视频一区二区| 我要看黄色一级片免费的| 欧美日韩在线观看h| 国产高清三级在线| 啦啦啦中文免费视频观看日本| 午夜福利网站1000一区二区三区| 狠狠精品人妻久久久久久综合| 多毛熟女@视频| 黑丝袜美女国产一区| 韩国高清视频一区二区三区| 丰满人妻一区二区三区视频av| 99国产精品免费福利视频| 女性被躁到高潮视频| 久久精品国产自在天天线| 欧美变态另类bdsm刘玥| 插阴视频在线观看视频| 少妇人妻一区二区三区视频| 免费人成在线观看视频色| 国产精品久久久久久av不卡| 日产精品乱码卡一卡2卡三| 又黄又爽又刺激的免费视频.| 免费观看性生交大片5| 国产永久视频网站| av在线app专区| 欧美性感艳星| 人妻夜夜爽99麻豆av| 久久久久久人妻| 亚洲精品国产成人久久av| 熟女电影av网| 特大巨黑吊av在线直播| 一级a做视频免费观看| 老司机影院毛片| 美女高潮的动态| 2018国产大陆天天弄谢| av天堂中文字幕网| 国产精品久久久久久精品古装| 日韩av免费高清视频| 欧美最新免费一区二区三区| 国产精品一二三区在线看| 欧美精品国产亚洲| 精品酒店卫生间| 国产乱来视频区| 中文精品一卡2卡3卡4更新| h日本视频在线播放| 国产精品国产三级国产专区5o| 在线观看三级黄色| 国产黄片美女视频| 啦啦啦啦在线视频资源| 久久久国产一区二区| 久久久久久久久久人人人人人人| 成人二区视频| 国产精品三级大全| 一级毛片黄色毛片免费观看视频| 欧美成人精品欧美一级黄| 伦理电影免费视频| 欧美三级亚洲精品| 久久久久久人妻| 一级毛片久久久久久久久女| 午夜福利高清视频| 一个人免费看片子| 美女福利国产在线 |