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

    珍稀瀕危樹種長序榆潛在適生區(qū)及其動態(tài)變化預測

    2021-11-28 11:10:14羅喻才陳禹錦余澤平楊光耀張文根
    南方林業(yè)科學 2021年5期
    關(guān)鍵詞:環(huán)境變量適生區(qū)貢獻率

    羅喻才,陳禹錦,余澤平,楊光耀,張文根★

    (1.中國科學院華南植物園,廣東 廣州510650;2.江西農(nóng)業(yè)大學·野生動植物保護與利用研究中心,江西 南昌330045;3.江西官山國家級自然保護區(qū)管理局,江西 宜春330063;4.中國科學院大學,北京100049)

    為適應(yīng)不斷變化的環(huán)境,物種或表現(xiàn)為分布格局的改變,或發(fā)生適應(yīng)性進化,更甚者則表現(xiàn)為物種瀕危和滅絕[1]。研究環(huán)境變化引起的物種分布格局變遷,了解物種分布與生態(tài)因子之間的關(guān)系,不僅能夠掌握瀕危物種的分布現(xiàn)狀,同時還能為物種的保護、引種以及馴化提供針對性策略[2-3]?;谖锓N分布模型和地理信息系統(tǒng)探討氣候變化對物種潛在適生區(qū)的影響已成為近些年生物多樣性和保護生物學研究熱點之一。在眾多物種分布模型中,MaxEnt利用分布數(shù)據(jù)與環(huán)境變量進行分析,探尋與此環(huán)境因子相似的環(huán)境像元,以此對目標物種的潛在地理分布進行無偏判斷,即使在樣本數(shù)量較少的情況下依舊能獲得較好的預測結(jié)果,為物種分布預測研究的熱點工具[4-7]。

    長序榆(Ulmus elongataL.K.Fu & C.S.Ding),是榆科(Ulmaceae)榆屬(Ulmus)在東亞的唯一代表種,主要分布在海拔400~900 m的東南季風氣候地區(qū),在研究北美和東亞植物地理學及系統(tǒng)演化中具有重要意義[8-9]。因受到長期的自然、人為因素干擾,長序榆所處生境破壞嚴重,野外種群數(shù)量較少,現(xiàn)已列入中國被子植物關(guān)鍵類群中高度瀕危種類和中國生物多樣性保護行動計劃中優(yōu)先保護的物種,同時也是國家Ⅱ級保護植物和全國極小種群野生保護植物(PSESP)[10-14]。目前,針對長序榆的保護生物學研究主要集中在其天然更新、群落特征分析、繁殖特性等方面[15-17],而關(guān)于環(huán)境變化引起長序榆適生區(qū)動態(tài)改變的相關(guān)研究尚未見報道。

    本研究基于長序榆的地理分布信息,結(jié)合氣候、地形、土壤等環(huán)境變量,利用MaxEnt模型和ArcGIS對長序榆當前和未來時期(2020—2040年)的潛在適生區(qū)進行了預測,旨在探討生態(tài)因子對長序榆分布區(qū)的影響以及在氣候變化下其適生區(qū)的變化規(guī)律,為長序榆的管理保護、引種栽培和開發(fā)利用提供參考。

    1 材料與方法

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

    基于國家標本資源共享平臺(http://www.nsii.org.cn/2017/home.php)、中國數(shù)字植物標本館(http://www.cvh.ac.cn/)、教學標本資源共享平臺(http://mnh.scu.edu.cn)和文獻資料,對所獲得的分布數(shù)據(jù)進行篩選,剔除重復、模糊的分布地點,并利用Google Earth定位系統(tǒng)提取坐標信息,共獲得長序榆分布點37份(表1)。

    表1 長序榆的分布信息Tab.1 The distribution imformation of U.elongata

    當前時期氣候以及海拔數(shù)據(jù)來自世界氣候數(shù)據(jù)庫(http://worldclim.org)。未來時期氣候數(shù)據(jù)采用聯(lián)合國政府間氣候變化專門委員會(IPCC)第6次氣候變化評估報告中公布的在未來2040年共享社會途徑(SSPs)下SSPs126和SSPs585兩種氣候情景模型,其中,在SSPs126情境下,全球溫度上升將限制在2℃以內(nèi);在SSPs585情境下,全球氣溫將平均升高5℃。氣候數(shù)據(jù)19個,即Bio_01至Bio_19,依次為年平均氣溫、晝夜溫差月均值、等溫性、溫度季節(jié)變化、最暖月最高溫度、最冷月最低溫度、年均溫變化范圍、最濕季度平均溫度、最干季度平均溫度、最暖季度平均溫度、最冷季度平均溫度、年均降水量、最濕月降水量、最干月降水量、降水量變異系數(shù)、最濕季度降水量、最干季度降水量、最暖季度降水量和最冷季度降水量。

    土壤數(shù)據(jù)來自世界土壤數(shù)據(jù)庫(HWSD)中的中國土壤數(shù)據(jù)集(http://data.casnw.net/portal/)。土壤因子16個:碎石體積百分比、沙含量、淤泥含量、粘土含量、USDA土壤質(zhì)地分類、土壤容重、有機碳含量、酸堿度、粘性層土壤的陽離子交換能力、土壤的陽離子交換能力、基本飽和度、交換性鹽基、碳酸鹽或石灰含量、硫酸鹽含量、可交換鈉鹽以及電導率。

    1.2 環(huán)境變量篩選及預處理

    為避免因環(huán)境變量間共線性導致的過度擬合,先將長序榆分布數(shù)據(jù)和環(huán)境變量載入MaxEnt模型,進行預模擬試驗。再利用ArcGIS將分布點與環(huán)境變量進行數(shù)值耦合,對耦合后的變量進行“Pearson”相關(guān)性分析?;陬A模擬試驗結(jié)果和變量之間的相關(guān)性分析,去除貢獻率為0以及相關(guān)系數(shù)|r|>0.8的2個變量中貢獻率較小的變量,最終獲得12個環(huán)境變量用于長序榆的適生區(qū)預測(表2)。

    1.3 適生區(qū)模擬及精度檢驗

    將長序榆分布數(shù)據(jù)以及環(huán)境變量載入MaxEnt模型,隨機選擇75%的長序榆分布點作為訓練集,另外25%作為測試集,選擇“刀切法(Jackknife)”計算各環(huán)境因子的貢獻率,并繪制相應(yīng)的環(huán)境響應(yīng)曲線,重復運算20次,得到20個隨機預測模型。利用受試者工作特征曲線(ROC)和AUC值(受試者工作特征曲線下的面積值)對模型預測的準確性進行檢測,AUC值越接近1,則表明預測結(jié)果越精確。

    1.4 適宜性等級劃分

    將MaxEnt模型輸出的結(jié)果導入ArcGIS,利用ArcGIS的ArcToolbox工具箱將其轉(zhuǎn)換成Raster格式。然后利用ArcToolbox空間分析工具中重分類命令的自然間斷點分級法將圖層的適生等級分為4類,分別為非適生區(qū)(適生概率P<0.1)、低適生區(qū)(適生概率0.1≤P<0.3)、一般適生區(qū)(適生概率0.3≤P<0.5)、中度適生區(qū)(適生概率0.5≤P<0.7)和高度適生區(qū)(適生概率P≥0.7),從而獲得長序榆的潛在適生分布圖。

    1.5 未來不同氣候背景下的適宜性動態(tài)變化

    在該研究中,分布概率小于0.1為非適生區(qū)。故以0.1為閾值,利用ArcGIS 10.2對不同時期的分布結(jié)果進行適宜性劃分。概率P≥0.1作為適生區(qū),以“1”表示,而概率P<0.1作為不適生區(qū),以“0”表示,通過適宜性劃分得到不同時期物種分布的二進制分布圖。將未來2040年的二進制分布圖和當前時期的二進制分布圖導入ArcGIS進行疊加運算,從而獲得當前至未來時期兩種情景下長序榆適生區(qū)動態(tài)變化圖。

    2 結(jié)果

    2.1 模型準確性檢驗

    在建模的過程中,利用ROC曲線以及AUC值來檢驗長序榆MaxEnt模型預測的準確性。一般認為,AUC值越接近1,模擬效果越好。其中,當0.5≤AUC<0.6時,模擬效果失??;0.6≤AUC<0.7,模擬效果較差;0.7≤AUC<0.8,模擬效果一般;0.8≤AUC<0.9,模擬效果良好;0.9≤AUC<1,模擬效果極好[19]。如圖1所示,在20次重復運算后,模型的平均AUC值為0.988,表明本研究中長序榆MaxEnt模型預測結(jié)果極好。

    圖1 長序榆MaxEnt模型的ROC曲線Fig.1 ROC curve of U.elongata MaxEnt model

    2.2 影響長序榆分布的環(huán)境變量

    在參與建模的12個環(huán)境變量中(表2),對長序榆潛在適生區(qū)分布影響最大的變量是水分,總貢獻率高達77.6%,其中較高的有最冷季度降水量(55.5%)、最干季度降水量(18.5%)以及最濕月降水量(1.9%)。其次是海拔因子,貢獻率為16.5%。溫度因子的總貢獻率為3.8%(溫度季節(jié)變化貢獻率2.2%,年平均氣溫貢獻率1.6%)。而其他環(huán)境因子,如粘性層土壤的陽離子交換能力、有機碳含量等對長序榆潛在適生區(qū)分布的貢獻率均較小。

    表2 環(huán)境變量對長序榆預測分布的相對貢獻Tab.2 Relative contribution of environmental variables to the predicted distribution of U.elongata

    長序榆的適生概率隨著每個環(huán)境因子的增加呈現(xiàn)先增后減的趨勢(圖2)。中度適生等級以上的分布區(qū)域(分布概率P>0.7),主要環(huán)境變量最冷季度降水量的適宜值范圍為192.07~535.7 mm,在243.5 mm時分布概率最大;最干季度降水量為158.84~359.72 mm,在209.064 mm處分布概率最大;海拔為454.4~980.72 m,在624.68 m時分布概率最大;溫度季節(jié)變化為687.65~814.22,在754.45時分布概率最大;最濕月降水量為275.79~483.05 mm,在366.89 mm時適生度最高;年平均溫為12.75~17.38℃,在14.62℃時適生度最高。

    圖2 主要環(huán)境因子的響應(yīng)曲線Fig.2 Response curve of major environmental factors

    2.3 長序榆的潛在適生區(qū)分布

    從圖3A和圖3B來看,長序榆現(xiàn)存分布點均在所預測的范圍內(nèi),表明模型預測得到的適生區(qū)與實際的適生區(qū)較為吻合,MaxEnt模型預測的可信度高。在當前適生區(qū)的4個適生等級中,高適生區(qū)的面積為15 952.5 km2,其核心分布在浙江省的西南部(慶元、遂昌、松陽等縣),呈塊狀分布,與江西、安徽、福建三省交界的區(qū)域(包括武夷山脈、白際山脈等)呈帶狀分布。中度適生區(qū)的面積為42 570 km2,除分布在高等適生區(qū)周邊外,廣西東北部、湘贛邊界的羅霄山脈、閩贛邊界的武夷山脈以及福建東北部等地區(qū)均有零星分布。其它等級適生區(qū)的面積較大,在廣西、浙江、福建、臺灣等省區(qū)呈間斷分布(圖3B)。

    圖3 長序榆分布、潛在適生區(qū)及未來適生區(qū)預測Fig.3 Potentially distribution of U.elongata in current and future

    2.4 未來長序榆的潛在適生區(qū)及其動態(tài)變化

    在未來2040年的不同氣候模式下,各等級適生區(qū)的面積均有不同程度的縮減(圖3和圖4)。其中,在SSPs126情景下,高適生區(qū)的面積為1 215 km2,相比于當前時期有大幅度縮減,當前時期的高度適生區(qū)降為中度適生區(qū)(共26 662.5 km2)。在SSPs585情景下,高適生區(qū)的面積為5 535 km2,中度適生區(qū)的面積為41 715 km2。其他等級的適生區(qū)面積相比于當前時期均有所下降。從整體來看,長序榆的潛在適生區(qū)分布面積在未來兩種情景下均呈現(xiàn)減少且向高緯度、高海拔地區(qū)移動的趨勢(圖5 A、B)。在SSPs126模式下,長序榆潛在適生面縮減了126 637.48 km2,其中,在江西、福建、湖南等地有明顯的減少,而在高緯度地區(qū)有小幅度擴張(2 713.99 km2)。在SSPs585模式下,長序榆潛在適生面積縮減了66 481.14 km2,并向內(nèi)陸高海拔地區(qū)擴張了29 641.57 km2(表3)。

    圖4 不同時期長序榆不同等級適宜分布區(qū)面積的比較Fig.4 Comparison of different grades of suitable areas of U.elongata under different climate scenarios

    圖5 未來時期長序榆適生區(qū)動態(tài)變化Fig.5 Possible change of potentially distribution of U.elongata in future

    表3 當前至2040年兩種氣候模式下長序榆的適生區(qū)面積變化Tab.3 Change of U.elongata suitable areas under two climate models from current to 2040

    3 討論

    本研究利用MaxEnt模型對瀕危樹種長序榆不同時期的潛在適生區(qū)進行了預測,分析了各環(huán)境變量對長序榆的重要性。從各環(huán)境因子的貢獻率來看,長序榆的生長主要受水分因素的影響,其貢獻率在參與建模的環(huán)境變量中占77.6%。研究表明,長序榆幼苗的光合作用對水分的缺失較為敏感,因此最冷季度、干旱季度降水量是影響長序榆幼苗生長、居群更新的限制因素[20]。除水分因素外,海拔和溫度也是影響長序榆生長的主要因素。長序榆主要分布在亞熱帶地區(qū)海拔600~900 m的山坡上,這些區(qū)域年均溫在16.2~19.5℃,與模型運算的結(jié)果大致相符[21]。土壤變量的貢獻率較低,但土壤對植物的生長至關(guān)重要,在對長序榆群落進行保護時應(yīng)綜合考慮各環(huán)境變量的關(guān)系,適當提高土壤肥力,有利于為長序榆成為群落優(yōu)勢種提供幫助[22]。

    從長序榆分布現(xiàn)狀和當前時期潛在適生區(qū)分布可以看出,其核心適生區(qū)(分布概率P>0.7)主要分布在浙江省杭州市清涼峰國家級自然保護區(qū)、浙江省麗水市九龍山國家級自然保護區(qū)、安徽省黃山市牯牛降國家級自然保護區(qū)以及江西、安徽、福建三省交界的地區(qū)。這些核心分布區(qū)呈狹窄、間斷的帶狀分布,居群分布碎片化嚴重。此外,影響物種分布和生境選擇的氣候因子對于珍稀瀕危植物的研究與管理至關(guān)重要,氣候的急劇變化往往導致許多野生植物種群數(shù)量不斷減少,成為珍稀瀕危植物[23-24]。本研究中,長序榆適生區(qū)面積在未來兩種不同氣候變化情景下均有不同程度的縮減。與當前時期相比,隨著未來氣溫的逐漸升高,長序榆在低緯度的適生區(qū)逐漸消失,整體向北擴張;同海拔地區(qū)的適生度也隨著氣溫升高而降低,高適生區(qū)面積大幅下降。

    因此,基于本文研究結(jié)果,對珍稀瀕危植物長序榆的保護提出如下建議:1)安徽南部、浙江西南及西北部、福建北部、江西的武夷山脈等地區(qū)為長序榆的核心分布區(qū),因此,除了對長序榆種群做就地保護和近地保護外,還需著重加強對這些區(qū)域的養(yǎng)護和管理,完善野生長序榆居群的本底摸查;2)目前雖然已對長序榆進行了一系列的引種栽培以及苗木繁育試驗,但引種區(qū)僅局限于現(xiàn)有的分布區(qū)內(nèi)[25],應(yīng)增加在其他適生等級較高的區(qū)域的引種栽培試驗;3)加強長序榆生長影響的生態(tài)因子研究,種子繁育及無性繁殖技術(shù)研究,探索其瀕危的生理、生態(tài)和分子機制;4)加大瀕危植物的科普宣傳力度,禁止亂砍濫伐,維護好現(xiàn)有長序榆生長環(huán)境。

    致謝:江西農(nóng)業(yè)大學陳伏生教授、葉學敏博士對本文提出建設(shè)性意見和建議。

    猜你喜歡
    環(huán)境變量適生區(qū)貢獻率
    基于最大熵模型的云南思茅松潛在分布區(qū)
    氣候變化下中國蒟蒻薯科箭根薯的地理分布格局預測
    未來氣候條件下當歸適生區(qū)預測及時空變化分析
    一種通用的裝備體系貢獻率評估框架
    從桌面右鍵菜單調(diào)用環(huán)境變量選項
    氣候變化下瀕危植物半日花在中國的潛在分布
    徹底弄懂Windows 10環(huán)境變量
    巴拉圭瓜多竹適生區(qū)分布研究
    竹子學報(2018年2期)2018-10-25 02:02:12
    關(guān)于裝備體系貢獻率研究的幾點思考
    В первой половине 2016 года вклад потребления в рост китайской экономики достиг 73,4 процента
    中亞信息(2016年10期)2016-02-13 02:32:45
    日本色播在线视频| 偷拍熟女少妇极品色| 熟女av电影| 国产精品伦人一区二区| 中文字幕精品免费在线观看视频 | 老司机影院毛片| 亚洲丝袜综合中文字幕| 国产91av在线免费观看| 人妻系列 视频| 久久精品夜色国产| 久久97久久精品| 一级毛片久久久久久久久女| 麻豆国产97在线/欧美| 久久国产乱子免费精品| 免费人妻精品一区二区三区视频| 亚洲激情五月婷婷啪啪| 亚洲一区二区三区欧美精品| 亚洲精品自拍成人| 国产高潮美女av| 亚洲欧美精品自产自拍| 欧美日韩综合久久久久久| 成人高潮视频无遮挡免费网站| 精品久久久噜噜| 亚洲av在线观看美女高潮| 亚洲成人一二三区av| 日本黄色片子视频| 18禁动态无遮挡网站| 国产黄色视频一区二区在线观看| 大香蕉97超碰在线| 成人无遮挡网站| 久久 成人 亚洲| 丰满迷人的少妇在线观看| 各种免费的搞黄视频| 少妇人妻 视频| 26uuu在线亚洲综合色| 精品国产乱码久久久久久小说| 久久99蜜桃精品久久| 欧美精品一区二区大全| 22中文网久久字幕| 人妻夜夜爽99麻豆av| 国产有黄有色有爽视频| 亚洲国产高清在线一区二区三| 国产有黄有色有爽视频| 三级经典国产精品| 99国产精品免费福利视频| 伦精品一区二区三区| 99re6热这里在线精品视频| 伦精品一区二区三区| 精品酒店卫生间| 亚洲成色77777| av天堂中文字幕网| 高清午夜精品一区二区三区| 国产av精品麻豆| 日本免费在线观看一区| 男的添女的下面高潮视频| 2022亚洲国产成人精品| 亚洲一区二区三区欧美精品| 深夜a级毛片| 国产精品一及| 又粗又硬又长又爽又黄的视频| 亚洲内射少妇av| 爱豆传媒免费全集在线观看| 啦啦啦在线观看免费高清www| 精品一区在线观看国产| 久久久国产一区二区| 高清av免费在线| 少妇裸体淫交视频免费看高清| 色视频在线一区二区三区| 国产成人a∨麻豆精品| 能在线免费看毛片的网站| 丰满少妇做爰视频| 在线免费十八禁| 最近的中文字幕免费完整| 少妇人妻精品综合一区二区| 精品人妻熟女av久视频| 免费观看的影片在线观看| 亚洲精品日本国产第一区| 国产老妇伦熟女老妇高清| 日韩 亚洲 欧美在线| 狠狠精品人妻久久久久久综合| 精品亚洲乱码少妇综合久久| 一个人免费看片子| 国产高清三级在线| 人妻 亚洲 视频| 国产成人freesex在线| av天堂中文字幕网| 男女边摸边吃奶| 哪个播放器可以免费观看大片| 黄片无遮挡物在线观看| 亚洲怡红院男人天堂| 亚洲欧洲日产国产| 狂野欧美激情性xxxx在线观看| 少妇的逼好多水| 18禁在线播放成人免费| 51国产日韩欧美| 久久精品久久久久久久性| 欧美日韩亚洲高清精品| av在线老鸭窝| 老女人水多毛片| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品一区蜜桃| 精品一品国产午夜福利视频| 久久这里有精品视频免费| 九草在线视频观看| 午夜福利高清视频| 久久久成人免费电影| 国产视频首页在线观看| 插阴视频在线观看视频| 最近2019中文字幕mv第一页| 免费看日本二区| 国产高清三级在线| 日韩在线高清观看一区二区三区| 亚洲欧美成人精品一区二区| 伦理电影大哥的女人| 成人毛片60女人毛片免费| 亚洲熟女精品中文字幕| av国产免费在线观看| 久久久成人免费电影| 九色成人免费人妻av| 一级毛片久久久久久久久女| 99久久综合免费| 一本一本综合久久| 美女视频免费永久观看网站| 午夜精品国产一区二区电影| 伊人久久国产一区二区| av视频免费观看在线观看| 网址你懂的国产日韩在线| 久久人妻熟女aⅴ| 蜜臀久久99精品久久宅男| 黄片无遮挡物在线观看| 亚洲综合色惰| 九色成人免费人妻av| av在线播放精品| 最后的刺客免费高清国语| 青春草国产在线视频| 狠狠精品人妻久久久久久综合| 成人综合一区亚洲| 国产探花极品一区二区| 国产免费福利视频在线观看| 国产黄频视频在线观看| 日日啪夜夜撸| 日本vs欧美在线观看视频 | 国产黄片美女视频| 22中文网久久字幕| 日日啪夜夜撸| 亚洲国产精品专区欧美| 成人特级av手机在线观看| 在线观看免费高清a一片| 毛片女人毛片| 欧美精品亚洲一区二区| 秋霞在线观看毛片| 尤物成人国产欧美一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 人人妻人人看人人澡| 黄色配什么色好看| 男女免费视频国产| 欧美精品一区二区免费开放| 国产精品av视频在线免费观看| 国产亚洲欧美精品永久| 久久6这里有精品| 亚洲欧美清纯卡通| 免费看av在线观看网站| 老女人水多毛片| 不卡视频在线观看欧美| 日本免费在线观看一区| 亚洲av欧美aⅴ国产| 少妇人妻久久综合中文| 成人漫画全彩无遮挡| 久久 成人 亚洲| 国产精品一二三区在线看| 在线观看一区二区三区激情| 亚洲国产精品999| 大话2 男鬼变身卡| 观看av在线不卡| 精品一区在线观看国产| 亚洲va在线va天堂va国产| 日本黄色片子视频| 久久婷婷青草| 国产白丝娇喘喷水9色精品| 精品人妻偷拍中文字幕| 国产一区二区在线观看日韩| 啦啦啦在线观看免费高清www| 日韩成人av中文字幕在线观看| 中文字幕精品免费在线观看视频 | videossex国产| 干丝袜人妻中文字幕| 精品一品国产午夜福利视频| av福利片在线观看| 美女主播在线视频| 亚洲av电影在线观看一区二区三区| av不卡在线播放| 亚洲精品国产av蜜桃| 日韩视频在线欧美| 99精国产麻豆久久婷婷| 亚洲欧美日韩无卡精品| 少妇的逼好多水| 91精品国产国语对白视频| 妹子高潮喷水视频| 80岁老熟妇乱子伦牲交| 国产精品久久久久久久电影| 久久99蜜桃精品久久| 亚洲精品色激情综合| 久久6这里有精品| 国产精品一二三区在线看| 精品一品国产午夜福利视频| 久久人妻熟女aⅴ| 亚洲国产av新网站| 激情五月婷婷亚洲| 女性被躁到高潮视频| 亚洲不卡免费看| 亚洲欧美日韩东京热| 国产一区二区三区av在线| 搡老乐熟女国产| 国产精品99久久99久久久不卡 | 日日摸夜夜添夜夜添av毛片| 国产av精品麻豆| 欧美bdsm另类| 免费观看无遮挡的男女| 国产一级毛片在线| 国产成人午夜福利电影在线观看| 我的老师免费观看完整版| 中文精品一卡2卡3卡4更新| 色视频在线一区二区三区| 久久精品国产亚洲网站| 久热这里只有精品99| 我的老师免费观看完整版| 免费久久久久久久精品成人欧美视频 | 91精品国产九色| 成人综合一区亚洲| 久久人妻熟女aⅴ| 亚洲综合色惰| 性色avwww在线观看| 欧美成人午夜免费资源| 欧美激情极品国产一区二区三区 | av免费在线看不卡| av国产久精品久网站免费入址| 热re99久久精品国产66热6| 久久这里有精品视频免费| 人妻夜夜爽99麻豆av| 亚洲av福利一区| 性色avwww在线观看| 免费黄网站久久成人精品| 国产在线视频一区二区| 在线观看免费日韩欧美大片 | av专区在线播放| 插逼视频在线观看| 草草在线视频免费看| 91久久精品国产一区二区成人| 久久久久性生活片| av视频免费观看在线观看| 欧美日韩综合久久久久久| 日韩不卡一区二区三区视频在线| 人人妻人人添人人爽欧美一区卜 | 亚洲精品乱码久久久v下载方式| 久久亚洲国产成人精品v| 亚洲中文av在线| 丝袜喷水一区| 日本av免费视频播放| 亚洲色图综合在线观看| 精品视频人人做人人爽| 偷拍熟女少妇极品色| 大码成人一级视频| 少妇裸体淫交视频免费看高清| 国产片特级美女逼逼视频| 国产 一区 欧美 日韩| 夜夜看夜夜爽夜夜摸| 久久久久久久久大av| 一区二区三区免费毛片| 亚洲美女搞黄在线观看| 国产精品99久久久久久久久| 在线观看免费日韩欧美大片 | 亚洲人成网站在线播| 国产成人午夜福利电影在线观看| 男女免费视频国产| 国产男女内射视频| 我要看黄色一级片免费的| 精品亚洲成国产av| 日本黄大片高清| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美日韩卡通动漫| 男人狂女人下面高潮的视频| 欧美激情国产日韩精品一区| av国产免费在线观看| 久久精品久久精品一区二区三区| 人妻系列 视频| 国产精品99久久久久久久久| 国产一区二区三区av在线| 搡老乐熟女国产| 亚洲成人一二三区av| 九草在线视频观看| 亚洲美女视频黄频| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久久精品古装| 午夜视频国产福利| 亚洲激情五月婷婷啪啪| 中国国产av一级| 熟女人妻精品中文字幕| 赤兔流量卡办理| 永久免费av网站大全| 777米奇影视久久| 亚洲四区av| 日韩一区二区视频免费看| 日韩一本色道免费dvd| 一级毛片电影观看| 最近的中文字幕免费完整| 人人妻人人看人人澡| 少妇人妻一区二区三区视频| 亚洲aⅴ乱码一区二区在线播放| 激情五月婷婷亚洲| 日韩亚洲欧美综合| 国产亚洲av片在线观看秒播厂| 九九在线视频观看精品| 高清日韩中文字幕在线| 我要看日韩黄色一级片| 最近2019中文字幕mv第一页| 人妻一区二区av| 观看av在线不卡| 一本久久精品| 亚洲丝袜综合中文字幕| 黄片无遮挡物在线观看| 91在线精品国自产拍蜜月| 狠狠精品人妻久久久久久综合| 亚洲久久久国产精品| 哪个播放器可以免费观看大片| 91久久精品国产一区二区成人| 伊人久久国产一区二区| 熟女电影av网| 欧美国产精品一级二级三级 | 国产老妇伦熟女老妇高清| 蜜桃久久精品国产亚洲av| 精品熟女少妇av免费看| 日韩av免费高清视频| 97精品久久久久久久久久精品| 久久av网站| 99久久精品热视频| a 毛片基地| 欧美人与善性xxx| 搡女人真爽免费视频火全软件| 午夜老司机福利剧场| 精品亚洲成a人片在线观看 | 啦啦啦视频在线资源免费观看| 亚洲国产精品国产精品| 寂寞人妻少妇视频99o| 激情 狠狠 欧美| 99热这里只有精品一区| 成人亚洲欧美一区二区av| 久久人人爽人人爽人人片va| 免费黄色在线免费观看| 午夜免费观看性视频| av播播在线观看一区| 国产精品久久久久久av不卡| 欧美97在线视频| 亚洲成人中文字幕在线播放| 午夜福利视频精品| 国产色婷婷99| 欧美日韩精品成人综合77777| 国语对白做爰xxxⅹ性视频网站| 成人美女网站在线观看视频| 国产精品偷伦视频观看了| 毛片一级片免费看久久久久| 久久精品国产鲁丝片午夜精品| 寂寞人妻少妇视频99o| 免费大片18禁| 国产深夜福利视频在线观看| 精品人妻偷拍中文字幕| 亚洲色图综合在线观看| 老女人水多毛片| 日本与韩国留学比较| 亚洲综合精品二区| 国产精品不卡视频一区二区| 九九久久精品国产亚洲av麻豆| 婷婷色综合www| 伦精品一区二区三区| 色吧在线观看| 韩国av在线不卡| 成人亚洲欧美一区二区av| 久久久久人妻精品一区果冻| 女的被弄到高潮叫床怎么办| 99视频精品全部免费 在线| 久久久精品94久久精品| 特大巨黑吊av在线直播| 1000部很黄的大片| 嫩草影院新地址| 黄色欧美视频在线观看| 成人黄色视频免费在线看| 男人和女人高潮做爰伦理| 国产午夜精品一二区理论片| 国产伦精品一区二区三区视频9| 水蜜桃什么品种好| 国产欧美亚洲国产| 亚洲av福利一区| 中文字幕制服av| 边亲边吃奶的免费视频| 黄色怎么调成土黄色| 国产成人精品一,二区| 少妇人妻一区二区三区视频| 日本欧美国产在线视频| 国产探花极品一区二区| 777米奇影视久久| 中文字幕免费在线视频6| 在线观看av片永久免费下载| 久久6这里有精品| 久久婷婷青草| 26uuu在线亚洲综合色| 欧美xxxx性猛交bbbb| 五月开心婷婷网| 夫妻午夜视频| 国产深夜福利视频在线观看| 午夜日本视频在线| 日韩欧美精品免费久久| 欧美+日韩+精品| 一边亲一边摸免费视频| 国产爽快片一区二区三区| 男女啪啪激烈高潮av片| 六月丁香七月| 国产伦理片在线播放av一区| 日韩大片免费观看网站| 香蕉精品网在线| 在线观看av片永久免费下载| 一级毛片我不卡| 亚洲欧美日韩东京热| 多毛熟女@视频| 最近最新中文字幕大全电影3| 日韩av不卡免费在线播放| 王馨瑶露胸无遮挡在线观看| 午夜免费观看性视频| 久久久久久久久大av| 亚洲精品456在线播放app| 欧美xxⅹ黑人| 日日摸夜夜添夜夜爱| 另类亚洲欧美激情| 国产精品一及| 亚洲成人一二三区av| 日韩欧美 国产精品| 欧美成人a在线观看| 日日摸夜夜添夜夜添av毛片| 亚洲精品色激情综合| 亚洲成人一二三区av| 午夜福利网站1000一区二区三区| 久久国产亚洲av麻豆专区| 亚洲欧美成人精品一区二区| 一区二区av电影网| 热re99久久精品国产66热6| 两个人的视频大全免费| 91精品伊人久久大香线蕉| 国产毛片在线视频| 亚洲精品色激情综合| 亚洲国产成人一精品久久久| 日韩av在线免费看完整版不卡| 欧美xxxx黑人xx丫x性爽| 亚洲久久久国产精品| 99热国产这里只有精品6| 春色校园在线视频观看| 人妻制服诱惑在线中文字幕| 观看av在线不卡| 精品久久久精品久久久| 麻豆成人av视频| 久久久国产一区二区| 久久久久久久久久久丰满| 男女无遮挡免费网站观看| 免费播放大片免费观看视频在线观看| 99久久精品国产国产毛片| 成人免费观看视频高清| 深爱激情五月婷婷| 免费人成在线观看视频色| 亚洲av.av天堂| 久久久亚洲精品成人影院| 久久婷婷青草| 免费播放大片免费观看视频在线观看| 精品人妻熟女av久视频| 国产白丝娇喘喷水9色精品| 激情 狠狠 欧美| 午夜免费观看性视频| 97超视频在线观看视频| 九九爱精品视频在线观看| 自拍偷自拍亚洲精品老妇| 亚洲欧美一区二区三区国产| 欧美日韩综合久久久久久| 最近中文字幕高清免费大全6| 国产成人a区在线观看| www.色视频.com| 纯流量卡能插随身wifi吗| 亚洲欧美日韩无卡精品| 国产中年淑女户外野战色| 永久免费av网站大全| 国产在视频线精品| 亚洲av中文字字幕乱码综合| 亚洲一级一片aⅴ在线观看| 免费黄网站久久成人精品| 自拍偷自拍亚洲精品老妇| 少妇 在线观看| 国产精品一区二区性色av| 亚洲精品国产av成人精品| 国产在线视频一区二区| 成人18禁高潮啪啪吃奶动态图 | 人妻制服诱惑在线中文字幕| 久久久久人妻精品一区果冻| 国产欧美日韩精品一区二区| 97超视频在线观看视频| 又爽又黄a免费视频| 久久精品久久久久久噜噜老黄| 国产无遮挡羞羞视频在线观看| 亚洲一区二区三区欧美精品| 亚洲欧美中文字幕日韩二区| 欧美三级亚洲精品| 成人无遮挡网站| 女性生殖器流出的白浆| 色视频在线一区二区三区| 欧美人与善性xxx| 偷拍熟女少妇极品色| 中文资源天堂在线| 欧美国产精品一级二级三级 | 日韩电影二区| 欧美精品国产亚洲| 亚洲av成人精品一区久久| 日日撸夜夜添| 在线亚洲精品国产二区图片欧美 | 免费黄网站久久成人精品| 免费人妻精品一区二区三区视频| 97精品久久久久久久久久精品| 身体一侧抽搐| 亚洲四区av| 免费观看av网站的网址| 一本久久精品| 高清午夜精品一区二区三区| 能在线免费看毛片的网站| 成人亚洲精品一区在线观看 | 久久精品国产亚洲网站| 亚洲av中文av极速乱| 国产中年淑女户外野战色| 中文字幕av成人在线电影| 简卡轻食公司| 男的添女的下面高潮视频| 亚洲精品日本国产第一区| 身体一侧抽搐| 麻豆成人午夜福利视频| av网站免费在线观看视频| 色婷婷久久久亚洲欧美| 午夜福利在线观看免费完整高清在| av国产久精品久网站免费入址| 人人妻人人澡人人爽人人夜夜| 国产成人免费无遮挡视频| 久久久久久久久久久免费av| 大片免费播放器 马上看| 亚洲成人av在线免费| 欧美成人一区二区免费高清观看| 国产精品一区二区在线观看99| 91久久精品国产一区二区成人| 岛国毛片在线播放| 少妇猛男粗大的猛烈进出视频| 亚洲色图av天堂| 国产精品久久久久久av不卡| 亚洲av二区三区四区| 亚洲婷婷狠狠爱综合网| 在线天堂最新版资源| 国语对白做爰xxxⅹ性视频网站| 欧美最新免费一区二区三区| 久久6这里有精品| 久久女婷五月综合色啪小说| 精品久久久精品久久久| 久久人妻熟女aⅴ| 久久久久精品性色| 精品人妻偷拍中文字幕| 亚洲欧美日韩另类电影网站 | 丰满乱子伦码专区| 日韩在线高清观看一区二区三区| 在线观看国产h片| 免费看光身美女| 91精品国产国语对白视频| 波野结衣二区三区在线| 一区二区三区乱码不卡18| 人体艺术视频欧美日本| 一级毛片我不卡| 三级国产精品欧美在线观看| 久久久久精品性色| 亚洲av中文av极速乱| 妹子高潮喷水视频| 亚洲成人中文字幕在线播放| 九九爱精品视频在线观看| 99热网站在线观看| 99热这里只有是精品50| 欧美国产精品一级二级三级 | 18+在线观看网站| 一级毛片久久久久久久久女| 看免费成人av毛片| 欧美日韩国产mv在线观看视频 | 一区二区三区免费毛片| 亚洲欧美一区二区三区黑人 | 免费播放大片免费观看视频在线观看| 超碰97精品在线观看| 国内揄拍国产精品人妻在线| 中文在线观看免费www的网站| 大码成人一级视频| a级一级毛片免费在线观看| 国产精品一区二区在线不卡| 日韩,欧美,国产一区二区三区| 少妇猛男粗大的猛烈进出视频| 精品国产一区二区三区久久久樱花 | .国产精品久久| 女性被躁到高潮视频| 汤姆久久久久久久影院中文字幕| 欧美3d第一页| 亚洲美女视频黄频| 国产日韩欧美在线精品| 舔av片在线| 五月开心婷婷网| 亚洲av免费高清在线观看| 亚洲av中文av极速乱| 狠狠精品人妻久久久久久综合| 色婷婷av一区二区三区视频| 免费少妇av软件| 狂野欧美激情性xxxx在线观看| 国产成人freesex在线| 少妇猛男粗大的猛烈进出视频| 国产有黄有色有爽视频|