• <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网| 亚洲,欧美精品.| 一个人免费在线观看电影 | 国产三级在线视频| 国产成人福利小说| 精品欧美国产一区二区三| 亚洲欧洲精品一区二区精品久久久| 9191精品国产免费久久| 久久欧美精品欧美久久欧美| 日韩欧美一区二区三区在线观看| 十八禁人妻一区二区| 国内精品久久久久久久电影| 国产1区2区3区精品| 国产极品精品免费视频能看的| 五月伊人婷婷丁香| 亚洲国产精品成人综合色| 成年女人毛片免费观看观看9| 看免费av毛片| 1000部很黄的大片| 美女高潮喷水抽搐中文字幕| 后天国语完整版免费观看| 国产综合懂色| 小说图片视频综合网站| 久久精品夜夜夜夜夜久久蜜豆| 99国产综合亚洲精品| bbb黄色大片| 日韩大尺度精品在线看网址| 日韩欧美 国产精品| 亚洲18禁久久av| 性色avwww在线观看| 法律面前人人平等表现在哪些方面| 国产成+人综合+亚洲专区| 97碰自拍视频| 成在线人永久免费视频| 亚洲精品456在线播放app | av视频在线观看入口| 老司机午夜十八禁免费视频| 中文资源天堂在线| 国产av在哪里看| 色综合站精品国产| av片东京热男人的天堂| 日韩国内少妇激情av| 村上凉子中文字幕在线| 黄色视频,在线免费观看| bbb黄色大片| 麻豆一二三区av精品| 我的老师免费观看完整版| 欧美黄色片欧美黄色片| 女警被强在线播放| 91九色精品人成在线观看| 此物有八面人人有两片| 我要搜黄色片| 精品久久久久久成人av| 国产亚洲av高清不卡| 国产精品久久久久久精品电影| 亚洲国产欧美一区二区综合| 亚洲在线观看片| 精品一区二区三区视频在线观看免费| 丰满人妻一区二区三区视频av | 国产日本99.免费观看| 亚洲av片天天在线观看| 日韩欧美精品v在线| 国产精品女同一区二区软件 | 婷婷六月久久综合丁香| a级毛片在线看网站| 国产成人啪精品午夜网站| 亚洲色图av天堂| 男女下面进入的视频免费午夜| 51午夜福利影视在线观看| www国产在线视频色| 亚洲自偷自拍图片 自拍| 精品久久久久久,| 偷拍熟女少妇极品色| 欧美丝袜亚洲另类 | 久久精品91无色码中文字幕| xxx96com| 一区二区三区高清视频在线| 琪琪午夜伦伦电影理论片6080| 久久久色成人| 中文字幕高清在线视频| 午夜福利欧美成人| 美女cb高潮喷水在线观看 | 成在线人永久免费视频| 美女午夜性视频免费| 精品一区二区三区四区五区乱码| 在线国产一区二区在线| 无遮挡黄片免费观看| 一个人免费在线观看电影 | 村上凉子中文字幕在线| 成人无遮挡网站| 天堂影院成人在线观看| 国产爱豆传媒在线观看| 欧美一区二区国产精品久久精品| 免费在线观看成人毛片| 国产乱人伦免费视频| 国产私拍福利视频在线观看| 香蕉国产在线看| 国产精品野战在线观看| 国产精品久久久av美女十八| www.www免费av| 亚洲av成人精品一区久久| 国产精品综合久久久久久久免费| av在线天堂中文字幕| 99久久精品国产亚洲精品| 制服人妻中文乱码| 婷婷精品国产亚洲av| 亚洲国产精品999在线| 深夜精品福利| 精品一区二区三区视频在线 | 午夜影院日韩av| 久久中文字幕一级| 成熟少妇高潮喷水视频| 小说图片视频综合网站| 国产一区二区三区在线臀色熟女| 午夜成年电影在线免费观看| 两个人的视频大全免费| 国产黄色小视频在线观看| 老司机福利观看| 国产精品98久久久久久宅男小说| 亚洲国产欧美人成| 成人av在线播放网站| 99久久精品一区二区三区| 国产爱豆传媒在线观看| 免费人成视频x8x8入口观看| 成人一区二区视频在线观看| 男人舔奶头视频| 免费观看人在逋| 精品久久久久久成人av| 999久久久国产精品视频| 舔av片在线| 91麻豆av在线| 两个人的视频大全免费| 久久久久久久精品吃奶| 久久热在线av| 亚洲自拍偷在线| 亚洲aⅴ乱码一区二区在线播放| 国产高清videossex| 国产一级毛片七仙女欲春2| 国产av不卡久久| 亚洲人成伊人成综合网2020| 老司机深夜福利视频在线观看| 国产精品久久久av美女十八| 麻豆国产av国片精品| 黑人操中国人逼视频| 亚洲成人中文字幕在线播放| 真人一进一出gif抽搐免费| 国产伦精品一区二区三区四那| 亚洲精品一区av在线观看| 国产高清有码在线观看视频| 免费观看精品视频网站| 少妇丰满av| 俺也久久电影网| 成人一区二区视频在线观看| www日本黄色视频网| 国产精品久久视频播放| 韩国av一区二区三区四区| www日本黄色视频网| 日本一二三区视频观看| 又黄又粗又硬又大视频| 国产蜜桃级精品一区二区三区| 麻豆国产av国片精品| 香蕉久久夜色| 亚洲熟妇中文字幕五十中出| 成人亚洲精品av一区二区| 中文亚洲av片在线观看爽| 国产精品九九99| 精品久久蜜臀av无| 国产乱人视频| 99国产精品一区二区蜜桃av| 99国产精品99久久久久| 九九久久精品国产亚洲av麻豆 | 精品不卡国产一区二区三区| 黄色视频,在线免费观看| 日韩欧美在线二视频| 亚洲国产高清在线一区二区三| 露出奶头的视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲熟女毛片儿| 琪琪午夜伦伦电影理论片6080| 老司机福利观看| 不卡一级毛片| 免费大片18禁| 怎么达到女性高潮| 国产蜜桃级精品一区二区三区| 熟女人妻精品中文字幕| 黑人巨大精品欧美一区二区mp4| 成人特级av手机在线观看| 亚洲国产中文字幕在线视频| 日本五十路高清| 久久亚洲真实| 欧美最黄视频在线播放免费| 久久这里只有精品中国| 99久久精品国产亚洲精品| h日本视频在线播放| 国产精品九九99| 亚洲第一欧美日韩一区二区三区| 久久久久久九九精品二区国产| 亚洲七黄色美女视频| 亚洲人成网站高清观看| 亚洲自拍偷在线| 国产精华一区二区三区| 又黄又爽又免费观看的视频| 国产一区二区三区在线臀色熟女| 怎么达到女性高潮| 久久天堂一区二区三区四区| 精品日产1卡2卡| 每晚都被弄得嗷嗷叫到高潮| 亚洲精华国产精华精| 日日摸夜夜添夜夜添小说| 美女免费视频网站| 中文字幕精品亚洲无线码一区| 亚洲av免费在线观看| 国产真实乱freesex| 国产精品乱码一区二三区的特点| 88av欧美| 国产成人啪精品午夜网站| 久久精品人妻少妇| 成年人黄色毛片网站| 国产伦在线观看视频一区| 欧美午夜高清在线| 99久国产av精品| 午夜视频精品福利| 一本久久中文字幕| 成人国产综合亚洲| 国产精品女同一区二区软件 | 免费看美女性在线毛片视频| 国产av麻豆久久久久久久| 色噜噜av男人的天堂激情| 国产真实乱freesex| 亚洲一区二区三区不卡视频| 亚洲国产精品合色在线| av在线天堂中文字幕| 日本一本二区三区精品| 搡老岳熟女国产| 精品一区二区三区视频在线 | 五月伊人婷婷丁香| 老司机午夜十八禁免费视频| 成人av一区二区三区在线看| 桃红色精品国产亚洲av| 国产精品98久久久久久宅男小说| 男女之事视频高清在线观看| 中文字幕av在线有码专区| 婷婷六月久久综合丁香| 欧美成人免费av一区二区三区| 日韩免费av在线播放| 麻豆av在线久日| 99精品在免费线老司机午夜| 一a级毛片在线观看| 可以在线观看的亚洲视频| 国产精品99久久99久久久不卡| 视频区欧美日本亚洲| 伦理电影免费视频| 男女午夜视频在线观看| 热99re8久久精品国产| 亚洲专区国产一区二区| 69av精品久久久久久| 国产亚洲精品综合一区在线观看| 美女黄网站色视频| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利高清视频| 亚洲精品久久国产高清桃花| 网址你懂的国产日韩在线| 露出奶头的视频| 亚洲人成网站高清观看| 国产私拍福利视频在线观看| 国产激情久久老熟女| 一个人看视频在线观看www免费 | 免费搜索国产男女视频| 国产熟女xx| 岛国视频午夜一区免费看| x7x7x7水蜜桃| 一进一出抽搐gif免费好疼| 免费人成视频x8x8入口观看| 熟妇人妻久久中文字幕3abv| 露出奶头的视频| 天堂影院成人在线观看| 一区福利在线观看| 欧美一级毛片孕妇| 国产精品久久久久久亚洲av鲁大| 麻豆国产av国片精品| 亚洲国产欧美网| 日韩欧美三级三区| 精品久久久久久久久久久久久| 大型黄色视频在线免费观看| 国产精华一区二区三区| 亚洲中文字幕日韩| 午夜福利在线观看吧| 又黄又粗又硬又大视频| 宅男免费午夜| 午夜福利成人在线免费观看| 99视频精品全部免费 在线 | 99视频精品全部免费 在线 | 久久精品人妻少妇| 亚洲成av人片免费观看| 国产乱人伦免费视频| 国产精品av视频在线免费观看| 久久这里只有精品19| 黑人巨大精品欧美一区二区mp4| 性色av乱码一区二区三区2| 丁香六月欧美| 国产精品一及| 日韩欧美国产在线观看| 九九久久精品国产亚洲av麻豆 | 成人国产一区最新在线观看| 亚洲国产色片| 亚洲电影在线观看av| 久久亚洲真实| 欧美不卡视频在线免费观看| 一个人免费在线观看的高清视频| 精品国产乱子伦一区二区三区| 欧美午夜高清在线| 中文字幕最新亚洲高清| 宅男免费午夜| 精品久久久久久久久久久久久| 亚洲国产看品久久| 少妇的丰满在线观看| 成年女人看的毛片在线观看| 搡老妇女老女人老熟妇| 日本熟妇午夜| 亚洲欧美一区二区三区黑人| 黑人操中国人逼视频| 琪琪午夜伦伦电影理论片6080| 国产精品一区二区免费欧美| 五月玫瑰六月丁香| 亚洲人成网站高清观看| 男女床上黄色一级片免费看| 99热这里只有是精品50| 国产蜜桃级精品一区二区三区| 久久草成人影院| 亚洲美女黄片视频| 国产亚洲精品久久久com| 黄色成人免费大全| а√天堂www在线а√下载| 一本综合久久免费| 一夜夜www| 人人妻人人看人人澡| 琪琪午夜伦伦电影理论片6080| 脱女人内裤的视频| 国产三级黄色录像| 午夜日韩欧美国产| 噜噜噜噜噜久久久久久91| 看片在线看免费视频| 亚洲精品美女久久av网站| 亚洲人成网站高清观看| 久久精品影院6| 欧美一级a爱片免费观看看| 国产一区在线观看成人免费| av片东京热男人的天堂| 欧美av亚洲av综合av国产av| 亚洲欧美精品综合久久99| 国产成人系列免费观看| 99久久综合精品五月天人人| 午夜a级毛片| 女人被狂操c到高潮| 久久久国产欧美日韩av| 99久久综合精品五月天人人| 精品国产三级普通话版| 欧美中文综合在线视频| 午夜福利18| 桃色一区二区三区在线观看| 亚洲欧美日韩无卡精品| 99久久综合精品五月天人人| 精品国产三级普通话版| 级片在线观看| 在线免费观看的www视频| 中文字幕最新亚洲高清| www.999成人在线观看| 亚洲av成人精品一区久久| 欧美黄色片欧美黄色片| 99精品在免费线老司机午夜| 综合色av麻豆| 成人永久免费在线观看视频| 中文字幕高清在线视频| 国产精品女同一区二区软件 | 免费av毛片视频| 久久久久精品国产欧美久久久| 热99在线观看视频| 亚洲国产中文字幕在线视频| 1024香蕉在线观看| 两个人看的免费小视频| 精品国产美女av久久久久小说| 国产人伦9x9x在线观看| 97超视频在线观看视频| 可以在线观看毛片的网站| 哪里可以看免费的av片| 午夜日韩欧美国产| 日本a在线网址| 国产精品一区二区三区四区免费观看 | 中文字幕人妻丝袜一区二区| 久久久精品大字幕| 国产三级黄色录像| 国产伦一二天堂av在线观看| 亚洲av第一区精品v没综合| 国内精品一区二区在线观看| 亚洲av第一区精品v没综合| 国产毛片a区久久久久| 国产精品 国内视频| 在线观看66精品国产| 在线观看午夜福利视频| 国产极品精品免费视频能看的| 最近在线观看免费完整版| 悠悠久久av| 国产成人av教育| 18禁观看日本| 国产69精品久久久久777片 | 国产成人啪精品午夜网站| 日韩欧美一区二区三区在线观看| 精品一区二区三区av网在线观看| 波多野结衣高清作品| 一级毛片精品| 亚洲18禁久久av| 国产av麻豆久久久久久久| 亚洲av电影在线进入| 国产精品99久久久久久久久| 久久亚洲真实| 不卡一级毛片| 日本撒尿小便嘘嘘汇集6| 人人妻,人人澡人人爽秒播| 亚洲欧美激情综合另类| 日韩 欧美 亚洲 中文字幕| 欧美国产日韩亚洲一区| 18禁国产床啪视频网站| 高潮久久久久久久久久久不卡| 国产精品日韩av在线免费观看| 亚洲中文日韩欧美视频| 国产高清有码在线观看视频| 国产精品亚洲一级av第二区| 亚洲天堂国产精品一区在线| 999久久久国产精品视频| 午夜激情福利司机影院| av国产免费在线观看| 国产精品久久电影中文字幕| 在线免费观看的www视频| 国产精品电影一区二区三区| 国产精品自产拍在线观看55亚洲| 国产毛片a区久久久久| 又大又爽又粗| 日本黄大片高清| 国产成人aa在线观看| 亚洲精品在线观看二区| 校园春色视频在线观看| 18禁国产床啪视频网站| 成人精品一区二区免费| 久久久精品欧美日韩精品| 成人高潮视频无遮挡免费网站| 丁香六月欧美| 女生性感内裤真人,穿戴方法视频| 黑人操中国人逼视频| 亚洲七黄色美女视频| 在线观看66精品国产| 欧美色视频一区免费| 大型黄色视频在线免费观看| 精品午夜福利视频在线观看一区| 国产私拍福利视频在线观看| 国产激情欧美一区二区| 18禁美女被吸乳视频| 国产精品98久久久久久宅男小说| 久9热在线精品视频| 国产三级中文精品| 99久国产av精品| 午夜免费观看网址| 天堂av国产一区二区熟女人妻| 91av网站免费观看| 国产亚洲av嫩草精品影院| 91九色精品人成在线观看| 国产成人精品无人区| 国产乱人伦免费视频| 欧美成狂野欧美在线观看| 欧美日韩黄片免| 在线免费观看不下载黄p国产 | 亚洲精品乱码久久久v下载方式 | 午夜激情欧美在线| 丁香欧美五月| 中文在线观看免费www的网站| 久久精品影院6| 国产伦一二天堂av在线观看| 免费看十八禁软件| 欧美+亚洲+日韩+国产| 精品午夜福利视频在线观看一区| 国产成年人精品一区二区| 99久久99久久久精品蜜桃| 日韩欧美在线二视频| 午夜a级毛片| 怎么达到女性高潮| 搡老妇女老女人老熟妇| 日韩中文字幕欧美一区二区| 日本撒尿小便嘘嘘汇集6| 久久精品aⅴ一区二区三区四区| 日韩大尺度精品在线看网址| 国产伦人伦偷精品视频| h日本视频在线播放| 别揉我奶头~嗯~啊~动态视频| 欧美日本亚洲视频在线播放| 亚洲天堂国产精品一区在线| a级毛片a级免费在线| 亚洲色图av天堂| 国产淫片久久久久久久久 | 又黄又粗又硬又大视频| 男人的好看免费观看在线视频| 国内少妇人妻偷人精品xxx网站 | 日本精品一区二区三区蜜桃| 99久久99久久久精品蜜桃| 日本 欧美在线| ponron亚洲| 亚洲精品456在线播放app | 国产熟女xx| 欧美成人一区二区免费高清观看 | 中出人妻视频一区二区| 亚洲18禁久久av| 欧美乱色亚洲激情| 国产亚洲精品久久久com| 人人妻人人澡欧美一区二区| 亚洲天堂国产精品一区在线| 听说在线观看完整版免费高清| 日本熟妇午夜| 最近在线观看免费完整版| 天堂动漫精品| 国产 一区 欧美 日韩| 一个人免费在线观看的高清视频| 黄色成人免费大全| 精品一区二区三区四区五区乱码| 久久99热这里只有精品18| 国产激情欧美一区二区| 免费人成视频x8x8入口观看| 精品电影一区二区在线| 精品久久久久久成人av| 国产亚洲欧美在线一区二区| 激情在线观看视频在线高清| 国产淫片久久久久久久久 | 成人国产综合亚洲| 成人av在线播放网站| 国产又色又爽无遮挡免费看| 国产欧美日韩精品一区二区| 一个人观看的视频www高清免费观看 | 88av欧美| 日本熟妇午夜| 欧美日韩一级在线毛片| 一a级毛片在线观看| 日日摸夜夜添夜夜添小说| 女人高潮潮喷娇喘18禁视频| svipshipincom国产片| 亚洲中文av在线| 老熟妇乱子伦视频在线观看| 国产麻豆成人av免费视频| 婷婷六月久久综合丁香| 亚洲av成人一区二区三| 丝袜人妻中文字幕| 久久久国产欧美日韩av| 熟女电影av网| 午夜亚洲福利在线播放| 精品国产超薄肉色丝袜足j| 久久精品国产99精品国产亚洲性色| 美女黄网站色视频| 老司机深夜福利视频在线观看| 后天国语完整版免费观看| 无遮挡黄片免费观看| 国产精品一区二区免费欧美| 日本a在线网址| 欧美一区二区精品小视频在线| 成熟少妇高潮喷水视频| 色综合站精品国产| 欧美日韩综合久久久久久 | ponron亚洲| 亚洲国产精品久久男人天堂| 视频区欧美日本亚洲| 精品不卡国产一区二区三区| 欧美黄色淫秽网站| 老汉色av国产亚洲站长工具| 久久久久国产精品人妻aⅴ院| a在线观看视频网站| 人人妻人人澡欧美一区二区| 久久久精品大字幕| 黄色片一级片一级黄色片| 久久天堂一区二区三区四区| 欧美乱色亚洲激情| 别揉我奶头~嗯~啊~动态视频| 国产精品av视频在线免费观看| 亚洲成av人片在线播放无| 免费观看的影片在线观看| 99国产极品粉嫩在线观看| 午夜久久久久精精品| 亚洲国产欧洲综合997久久,| 成年版毛片免费区| 非洲黑人性xxxx精品又粗又长| 黄色女人牲交| 久久久久久久久中文| 日本一二三区视频观看| 中文字幕人成人乱码亚洲影| 国产精品野战在线观看| 丝袜人妻中文字幕| 精华霜和精华液先用哪个| 久久人妻av系列| 51午夜福利影视在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久伊人香网站| 免费大片18禁| 国内少妇人妻偷人精品xxx网站 | 淫秽高清视频在线观看| 动漫黄色视频在线观看| 美女扒开内裤让男人捅视频| 亚洲性夜色夜夜综合| 身体一侧抽搐| 成人18禁在线播放| 观看美女的网站| 老汉色∧v一级毛片| 国产精品永久免费网站| 观看美女的网站| 亚洲 国产 在线| 91老司机精品| a级毛片在线看网站| 日韩成人在线观看一区二区三区| 精品国内亚洲2022精品成人| 一边摸一边抽搐一进一小说|