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

    1981—2017年黃河流域內(nèi)蒙古段氣候及植被生態(tài)變化狀況

    2021-07-02 00:56:36桑婧
    內(nèi)蒙古氣象 2021年1期
    關(guān)鍵詞:趨勢區(qū)域

    桑婧

    (內(nèi)蒙古生態(tài)與農(nóng)業(yè)氣象中心,內(nèi)蒙古 呼和浩特 010051)

    關(guān)鍵字 黃河流域內(nèi)蒙古段;NDVI;氣候變化;Mann-Kendall檢驗

    引言

    植被泛指地球表面某一區(qū)域所覆蓋的植物群落,是生態(tài)系統(tǒng)的重要組成部分之一,連接了大氣、土壤和水,在圈層間物質(zhì)循環(huán)和能量、信息傳遞中起到了溝通的作用[1]。在全球范圍內(nèi),植被在保持土壤水分、調(diào)節(jié)大氣氣候、維穩(wěn)生態(tài)環(huán)境等方面具有不可忽視的作用。植被歸一化指數(shù)(NDVI)與植物的生物量和葉面積指數(shù)有相對較好的相關(guān)關(guān)系,隨著遙感技術(shù)的發(fā)展,目前NDVI被廣泛地應(yīng)用于水土保持和生態(tài)環(huán)境的研究中,可以代表地面植被覆蓋程度的變化[2]。NOAA/AVHHRR和SPOT數(shù)據(jù)具有覆蓋范圍廣,分辨率高且獲取方便的優(yōu)勢,因此適合作為數(shù)據(jù)源研究長時間序列NDVI變化[3]。

    黃河流域地理位置和氣候條件較為特殊,加之頻繁的人類活動,流域內(nèi)水土流失嚴(yán)重導(dǎo)致黃河流域生態(tài)環(huán)境十分脆弱,研究黃河流域內(nèi)植被覆蓋的變化對保護(hù)黃河流域生態(tài)環(huán)境具有重要意義[4-5]。新中國以來,我國在黃河流域植被覆蓋變化領(lǐng)域做了大量的研究工作,得到了許多有意義的結(jié)論。袁麗華等[3]利用2000—2010 年的MODI3QI數(shù)據(jù)研究了黃河流域植被覆蓋區(qū)域的NDVI時空變化特征,結(jié)果發(fā)現(xiàn),2000—2010 年間黃河流域NDVI年值呈現(xiàn)出西部和東南部高,均值在0.3~0.4之間,黃河流域植被改善的區(qū)域遠(yuǎn)遠(yuǎn)大于退化的區(qū)域。賀振等[2]利用GIMMS-NDVI數(shù)據(jù)分析了1982—2013黃河流域植被覆蓋時空演化趨勢,發(fā)現(xiàn),黃河流域植被覆蓋增長較為緩慢,增速為每10 年增長0.018,32 年間黃河流域植被覆蓋整體得到改善,得到改善的區(qū)域約為59.49%。田智慧等[6]基于MOD17A3 NPP數(shù)據(jù)及氣象和土地利用數(shù)據(jù)分析了2000—2015年植被NPP變化,結(jié)果表明,上、中、下游植被NPP年均值呈明顯的梯度分布,即上游<中游<下游,說明中上游區(qū)域生態(tài)環(huán)境比較脆弱,且中上游植被NPP總量占整個流域的96%,可見中上游生態(tài)環(huán)境對整個黃河流域生態(tài)環(huán)境的影響較大,故加強(qiáng)對中上游流域生態(tài)環(huán)境保護(hù)和改善至關(guān)重要。郭帥等[7]基于GIMMS NDVI分析了黃河流域NDVI的時空變化,發(fā)現(xiàn),1982—2015年黃河流域及上游、中游NDVI均呈現(xiàn)顯著線性增加的趨勢,中游的增速最快,顯著增加的區(qū)域面積也最大。

    綜上所述,雖然前人對黃河流域植被覆蓋變化趨勢做了較多的研究,但鮮見長時間序列下只針對黃河流域內(nèi)蒙古段的研究。本研究融合多源遙感數(shù)據(jù),通過預(yù)處理得到了1981—2017 年黃河流域年最大NDVI數(shù)據(jù),采用了目前被廣泛認(rèn)可的Mann-Kendall趨勢檢驗研究黃河流域內(nèi)蒙古段植被覆蓋空間分布特征和變化規(guī)律,以期為黃河流域內(nèi)蒙古段生態(tài)環(huán)境保護(hù)和治理提供科學(xué)依據(jù)。

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

    1.1 研究區(qū)概況

    黃河發(fā)源于巴顏喀拉山北麓各姿各雅山下的卡日曲河谷,海拔高度4675 m,平均流量1774.5 m3·s-1,一路歷經(jīng)5464 km,最后在山東省注入渤海。內(nèi)蒙古自治區(qū)范圍內(nèi)黃河干流830 km,約占黃河中上游流域的17.7%,流經(jīng)內(nèi)蒙古7 個盟市。黃河流域內(nèi)蒙古段地貌多樣,屬于干旱半干旱氣候區(qū),常年氣候干旱寒冷,流經(jīng)區(qū)域多為荒漠和半荒漠,氣候干燥再加上人類活動的影響使黃河生態(tài)環(huán)境較為脆弱。

    1.2 數(shù)據(jù)來源及處理

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

    1981—2001 年的NOAA/AVHHRR逐旬NDVI數(shù)據(jù)(分辨率為2.2 km)來源于中國農(nóng)業(yè)科學(xué)研究院,1998—2017 年的SPOT-Vegetation NDVI(分辨率為1.0 km)數(shù)據(jù)來源于美國國家海洋和大氣管理局官網(wǎng)。黃河流域1:250000二級流域分級數(shù)據(jù)集(2002 年)來源于國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心。黃河流域內(nèi)蒙古段及周邊內(nèi)蒙古境內(nèi)44 個氣象站點對應(yīng)年份的年平均氣溫、年降水量和相對濕度數(shù)據(jù)(黃河流域內(nèi)蒙古段區(qū)域及氣象站點位置見圖1)。

    圖1 黃河流域內(nèi)蒙古段位置及氣象站點分布情況

    1.2.2 非同源遙感影像處理

    NOAA/AVHRR NDVI(1981—2001)、SPOTVegetation(1998—2017)數(shù)據(jù)集在時間長度和空間分辨率方面各具優(yōu)勢,通過重采樣(ArcTool box完成)、相關(guān)分析和最小二乘擬合等方法,依據(jù)重合時段(1998—2001年)的40 期月值數(shù)據(jù)的對比分析,將以上非同源數(shù)據(jù)序列進(jìn)行融合,構(gòu)建1981—2017年長時間序列的植被遙感數(shù)據(jù)集。而后將1981—2017 年長時間序列的NDVI數(shù)據(jù)進(jìn)行年最大值合成。

    最小二乘法在估計誤差、擬合多源數(shù)據(jù)、預(yù)報預(yù)測等許多學(xué)科領(lǐng)域得到廣泛應(yīng)用[10]。最小二乘法(又稱最小平方法)屬于一種數(shù)學(xué)優(yōu)化技術(shù),它通過最小化的誤差的平方和已達(dá)到尋找數(shù)據(jù)的最佳函數(shù)匹配。利用最小二乘法可以簡便地通過已知數(shù)據(jù)求得未知的數(shù)據(jù),并使得求得的數(shù)據(jù)與實際數(shù)據(jù)之間誤差的平方和為最小。根據(jù)最小二乘原理,本研究通過重合時段的數(shù)據(jù)求得a1和a0,代入到原公式Y(jié)=a1X+a0中,其中Y為校準(zhǔn)后的遙感數(shù)據(jù),X為需校準(zhǔn)的遙感數(shù)據(jù):

    式中:n= 40。xi、yj分別為重合時段對應(yīng)NOAA和SPOT遙感數(shù)據(jù)。

    1.3 MK檢驗計算方法

    Mann-Kendall趨勢檢驗是目前被廣泛認(rèn)可的可以使用在氣象和生態(tài)研究的領(lǐng)域的趨勢檢驗方法。Mann-Kendall趨勢檢驗適用于趨勢可能是單調(diào)的,因此數(shù)據(jù)中不存在周期性的情況,該檢驗允許缺少值,并且數(shù)據(jù)不需要符合正態(tài)分布[11]。計算公式如下:

    首先,S的方差由公式(3)計算,q是樣本的數(shù)量,tp是樣本中第p個樣本的值。

    S和VAR(S)的值用于計算測試統(tǒng)計量MK:

    MK值評估在統(tǒng)計上顯示趨勢,MK的正值表示上升趨勢,負(fù)值表示下降趨勢。且統(tǒng)計量MK具有正態(tài)分布,取值范圍為(-∞,+∞)。為了測試在顯著性水平為a的向上或向下單調(diào)趨勢(雙尾檢驗),在給定顯著性水平α下,當(dāng)|MK|>u1-α/2時,表示研究序列在α水平上存在顯著的變化。在Mann-Kendall檢驗的顯著性水平α為0.001、0.01、0.05和0.1,常用的是0.1和0.05,即90%和95%置信水平。本研究取α=0.05,本文判斷在95%置信水平上NDVI、年降水量、年平均氣溫和相對濕度時間序列變化趨勢的顯著性。

    對于研究區(qū)內(nèi)氣象站點的氣象數(shù)據(jù)的Mann-Kendall趨勢檢驗結(jié)果,我們采用反距離權(quán)重法進(jìn)行插值形成格點數(shù)據(jù)(分辨率為2.2 km)。對于NDVI數(shù)據(jù),使用Matlab等工具,根據(jù)Mann-Kendall檢驗公式進(jìn)行編程,將柵格數(shù)據(jù)轉(zhuǎn)化成ASCII碼進(jìn)行計算,然后格點結(jié)果在Arcgis軟件上進(jìn)行繪制。

    2 結(jié)果與分析

    2.1 黃河流域內(nèi)蒙古段氣候條件變化趨勢分析

    黃河流域內(nèi)蒙古段和周邊內(nèi)蒙古境內(nèi)氣象站點年平均氣溫、年降水量和年平均相對濕度Mann-Kendall檢驗結(jié)果(圖2),黃河流域內(nèi)蒙段區(qū)域年平均氣溫和年降水量變化(圖3)。從圖2a和圖3a可見,年降水量在空間上基本沒有明顯的上升或者下降趨勢,區(qū)域平均年降水量在時間上略有上升趨勢。區(qū)域年降水量最高值是4017.5 mm(2012 年),最低值是1925.9 mm(2005 年),最高值和最低值相差一倍多。從圖2b和圖3b可見,年平均氣溫除了個別氣象站點的變化趨勢不明顯之外,其他的站點均呈現(xiàn)顯著增加的趨勢,且區(qū)域年平均氣溫也呈現(xiàn)明顯的上升趨勢,氣候傾向斜率為0.048 ℃/a,通過了95%置信水平檢驗。從時間上看,區(qū)域年平均氣溫最高值為8.1 ℃(1998年),最低值為5.1 ℃(1984 年),1981—2017 年氣溫總體呈現(xiàn)波動上升的趨勢。

    圖2 黃河流域內(nèi)蒙古段1981—2017 年年降水量(a)、年平均氣溫(b)和年平均相對濕度(c)變化趨勢

    從圖2c可以看出,總體上,年平均相對濕度在空間上呈現(xiàn)北部明顯減少,中部和南部基本沒有明顯變化。在呼和浩特市、包頭市和鄂爾多斯市三市交界部分、河套平原和庫布齊沙漠的東部、烏蘭布和沙漠的東南部年相對濕度均呈現(xiàn)顯著的降低。從時間上看(圖3c),區(qū)域的年平均相對濕度總體呈現(xiàn)波動下降的趨勢,但趨勢并不明顯。最高值為57%(2003年),最低值為46%(2013年),近10 年的相對濕度平均水平較低,為48%,低于1981—2017年的平均值50%。

    圖3 黃河流域內(nèi)蒙古段1981—2017 年區(qū)域年降水量(a)、年平均氣溫(b)和年平均相對濕度(c)變化

    綜上所述,黃河流域內(nèi)蒙古段和周邊內(nèi)蒙古境內(nèi)氣象站點的年平均氣溫?zé)o論從空間還是時間上均呈現(xiàn)明顯的上升趨勢,氣候傾向斜率為0.048 ℃/a;年降水量在時間上和空間上均沒有明顯的上升和下降的趨勢,但降水分布不均,區(qū)域年降水量最高值和最低值相差一倍多;年平均相對濕度在空間上主要表現(xiàn)為北部明顯減少,中部和南部無明顯變化趨勢,時間上呈現(xiàn)波動下降的趨勢。

    2.2 黃河流域內(nèi)蒙古段植被覆蓋及其空間分布

    利用1981—2017 年年最大值合成的NDVI數(shù)據(jù)進(jìn)行平均值計算,得到了37 年黃河流域內(nèi)蒙古段的植被歸一化指數(shù)空間分布(圖4)。從圖4中可以得出,總體上,黃河流域內(nèi)蒙古段的NDVI值以0.1~0.4的低值為主,總體呈現(xiàn)從北到南,低到高再到低的變化過程,即西北和南部較低,大部分在0.4以內(nèi);中部偏北的地區(qū)較高,大部在0.4以上。黃河流域內(nèi)蒙古段的高值區(qū)主要位于黃河干流以北,陰山山脈以南的區(qū)域,此地的地形較為平坦,土壤較為肥沃,有利于植被的生長。

    圖4 黃河流域內(nèi)蒙古段1981—2017年植被歸一化指數(shù)平均年最大值空間分布

    黃河流域內(nèi)蒙古段所流經(jīng)的土地主要以草地和耕地為主,有少部分荒地和林地[10]。其中耕地和林地的植被類型較為豐富,生態(tài)環(huán)境較好,所以耕地和林地的NDVI值較高,大部分在0.4以上,小部分地區(qū)甚至可以達(dá)到0.7以上。草地區(qū)域的NDVI值較低,主要在0.1~0.3之間,由于草地和荒地是穿插分布,所以部分NDVI<0.1的區(qū)域基本與黃河流域內(nèi)蒙古段荒地出現(xiàn)的范圍吻合,荒地的生態(tài)環(huán)境脆弱植被稀少,有些是裸露的沙地甚至是荒漠。

    對1981—2017年NDVI年最大值平均的像元分級統(tǒng)計結(jié)果顯示:NDVI<0.1的無植被覆蓋區(qū)面積為2910 km2,占黃河流域內(nèi)蒙古段的2%,NDVI>0.1的有植被覆蓋區(qū)域面積占98%。其中NDVI在0.1~0.4的低值區(qū)面積為108295 km2,約占黃河流域內(nèi)蒙古段的74.4%;NDVI>0.4的高值區(qū)面積為34270 km2,約占該流域面積的23.6%,其中NDVI>0.6的高值區(qū)面積為3470 km2,約占該流域面積的2.39%。NDVI值像元分級面積(表1)。

    表1 NDVI各級面積和所占該流域面積比例

    2.3 黃河流域內(nèi)蒙古段植被覆蓋變化趨勢分析

    黃河流域內(nèi)蒙古段植被歸一化指數(shù)(NDVI)Mann-Kendall檢驗結(jié)果如表2和圖5所示,將Mann-Kendall趨勢檢驗按照95%置信水平劃分為顯著變化和無明顯變化,即MK>1.96顯著增加,-1.96≤MK≤1.96無明顯變化,MK<-1.96顯著減少,NDVI變化趨勢可以明顯看出像元尺度上,黃河流域內(nèi)蒙古段的NDVI值顯著增加的區(qū)域占該流域總面積的43.39%,顯著減少區(qū)域面積占該流域面積的38.70%,無明顯變化區(qū)域面積占該流域總面積的17.91%,顯著增加區(qū)域面積占比最高,無明顯變化其次,顯著減少最少,總體上來說黃河流域內(nèi)蒙古段植被歸一化指數(shù)是增加趨勢。

    表2 NDVI變化趨勢統(tǒng)計

    從圖5可見,顯著增加的區(qū)域主要分布在河套平原、呼和浩特市、包頭市和鄂爾多斯市三市交界處;無明顯變化的區(qū)域主要分布在陰山山脈、鄂爾多斯高原;顯著減少區(qū)域主要分布在大青山、庫布齊沙漠。在毛烏素沙地零散分布著顯著增加和顯著減少的區(qū)域,可能是由于毛烏素沙地處于幾個自然地帶的交接處,植被和土壤都呈現(xiàn)過渡性特點,分布交錯類型較為復(fù)雜。

    圖5 黃河流域內(nèi)蒙古段1981—2017年植被歸一化指數(shù)年最大值 變化趨勢

    結(jié)合土地利用類型來看,黃河流域內(nèi)蒙古段的耕地NDVI值以顯著增加為主,包括有“塞外江南”和“塞上糧倉”之稱的河套平原、以及大青山南麓;林地和草地的NDVI值主要以無明顯變化為主;荒地和沙地主要以顯著減少為主。

    結(jié)合年平均值和變化趨勢圖來看,在NDVI年平均值較高的區(qū)域,變化趨勢基本以顯著增加為主,從河套平原和大青山南麓可以明顯的看出這一結(jié)果。在NDVI值較低的地區(qū),特別是鄂爾多斯高原的西北、河套平原的西北,NDVI的變化趨勢基本是顯著減少為主。值得關(guān)注的是在準(zhǔn)格爾旗和呼和浩特市的交界處,NDVI平均值是較低的,在0.2~0.3之間為主,但NDVI的趨勢卻是以顯著增加為主,這說明該地的植被覆蓋程度雖然較差,卻在逐步改善中。

    3 結(jié)論與討論

    3.1 結(jié)論

    (1)黃河流域內(nèi)蒙古段和周邊內(nèi)蒙古境內(nèi)氣象站點的年平均氣溫呈現(xiàn)明顯的上升趨勢,氣候傾向斜率為0.048 ℃/a;年降水量在時間上和空間上均沒有明顯的上升和下降的趨勢;年平均相對濕度在空間上主要表現(xiàn)為北部明顯減少,中部和南部無明顯變化趨勢,時間上呈現(xiàn)波動下降的趨勢。

    (2)從空間上看,黃河流域內(nèi)蒙古段總體呈現(xiàn)從北到南,NDVI值由低到高再變低的變化過程,即西北和南部較低,以0.1~0.4低值為主;中部偏北的地區(qū)較高,大部在0.4以上。黃河流域內(nèi)蒙古段的NDVI值顯著增加的區(qū)域面積占整個流域的43.39%多于顯著減少區(qū)域(面積占比僅為17.91%)。

    (3)耕地的NDVI值較高,大部分在0.4以上,且以NDVI值從1981—2017年顯著增加為主;林地的NDVI值大部分也在0.4以上,但從1981—2017年無明顯的變化趨勢;草地和荒地的NDVI值較低,在0.1~0.3之間,以顯著減少為主。

    (4)NDVI年平均值較高的區(qū)域,變化趨勢基本也以顯著增加為主,在NDVI值較低的地區(qū),NDVI的變化趨勢基本是顯著減少為主,但準(zhǔn)格爾旗和呼和浩特市的交界的植被覆蓋程度雖然較差,改善的趨勢卻十分明顯。

    3.2 討論

    本研究使用NOAA/AVHHRR和SPOT-Vegetation數(shù)據(jù)制作成長時間序列的NDVI數(shù)據(jù),利用Mann-Kendall趨勢檢驗對黃河流域內(nèi)蒙古段植被變化趨勢特征進(jìn)行分析。由于時間序列較長,可以更有效地展現(xiàn)黃河流域內(nèi)蒙古段植被時間變化的趨勢特征,具有一定的說服力。同時本研究也有一些局限性,比如,沒有針對不用的植被類型進(jìn)行分析;沒有進(jìn)行實地考察,將考察結(jié)果與Mann-Kendall趨勢檢驗結(jié)果進(jìn)行對比;對于由于毛烏素沙地土地類型較為復(fù)雜,Mann-Kendall趨勢檢驗的結(jié)果也出現(xiàn)了增加和減少交錯分布的情況,對此情況還需調(diào)查后給出合理的解釋。

    猜你喜歡
    趨勢區(qū)域
    永久基本農(nóng)田集中區(qū)域“禁廢”
    趨勢
    分割區(qū)域
    初秋唇妝趨勢
    Coco薇(2017年9期)2017-09-07 21:23:49
    SPINEXPO?2017春夏流行趨勢
    關(guān)于四色猜想
    分區(qū)域
    “去編”大趨勢
    基于嚴(yán)重區(qū)域的多PCC點暫降頻次估計
    電測與儀表(2015年5期)2015-04-09 11:30:52
    趨勢
    汽車科技(2015年1期)2015-02-28 12:14:44
    大型黄色视频在线免费观看| 99热这里只有精品一区| 日本黄色片子视频| 婷婷色av中文字幕| 热99re8久久精品国产| 人妻系列 视频| 国产v大片淫在线免费观看| 桃色一区二区三区在线观看| 国产老妇伦熟女老妇高清| 国产成人影院久久av| 卡戴珊不雅视频在线播放| 91久久精品国产一区二区成人| 69av精品久久久久久| 毛片一级片免费看久久久久| 18禁在线播放成人免费| 久久久精品94久久精品| 精品一区二区免费观看| 久久国产乱子免费精品| 少妇丰满av| 亚洲精品成人久久久久久| 91av网一区二区| av免费观看日本| 最新中文字幕久久久久| 午夜a级毛片| av免费在线看不卡| 国产一区亚洲一区在线观看| 哪里可以看免费的av片| 中文欧美无线码| 久久久成人免费电影| 精华霜和精华液先用哪个| 亚洲一区二区三区色噜噜| 成人高潮视频无遮挡免费网站| 久久精品国产亚洲av涩爱 | 日本免费一区二区三区高清不卡| 国产高潮美女av| 久久热精品热| 神马国产精品三级电影在线观看| av天堂在线播放| 亚洲综合色惰| 国产一区二区激情短视频| 日本av手机在线免费观看| 午夜免费激情av| 一级黄片播放器| 欧美成人精品欧美一级黄| 高清在线视频一区二区三区 | 99久久精品国产国产毛片| 亚洲欧美清纯卡通| 国产麻豆成人av免费视频| 国产精品蜜桃在线观看 | 久久久欧美国产精品| 国产在线男女| 国产黄片美女视频| 全区人妻精品视频| 蜜臀久久99精品久久宅男| 亚洲欧洲国产日韩| 亚洲精品色激情综合| 国产精品精品国产色婷婷| 欧美xxxx黑人xx丫x性爽| 久久人人爽人人片av| 久久久久久久久久久丰满| a级一级毛片免费在线观看| 大香蕉久久网| av在线蜜桃| www日本黄色视频网| 国产亚洲91精品色在线| 国产黄色小视频在线观看| 亚洲精品国产av成人精品| 99热全是精品| 国产大屁股一区二区在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 两个人的视频大全免费| 一区福利在线观看| 亚洲国产欧洲综合997久久,| 欧美一区二区亚洲| 国产成人a∨麻豆精品| 一本一本综合久久| 亚洲精品乱码久久久v下载方式| 在线播放国产精品三级| 亚洲国产精品国产精品| 大又大粗又爽又黄少妇毛片口| 一个人看视频在线观看www免费| 在线天堂最新版资源| 亚洲精品亚洲一区二区| 国内揄拍国产精品人妻在线| 亚洲国产欧洲综合997久久,| 美女xxoo啪啪120秒动态图| 一级毛片我不卡| 色视频www国产| 亚洲国产精品国产精品| 亚洲在久久综合| 亚洲电影在线观看av| 2021天堂中文幕一二区在线观| av又黄又爽大尺度在线免费看 | 久久婷婷人人爽人人干人人爱| 我要看日韩黄色一级片| 国内精品宾馆在线| 99热这里只有精品一区| 97超碰精品成人国产| 插逼视频在线观看| 久久久久久九九精品二区国产| eeuss影院久久| 床上黄色一级片| 国产高清有码在线观看视频| 最近中文字幕高清免费大全6| 能在线免费看毛片的网站| 日韩欧美精品免费久久| 99热这里只有是精品在线观看| 免费观看精品视频网站| 免费看光身美女| 一区二区三区高清视频在线| 日韩欧美 国产精品| 菩萨蛮人人尽说江南好唐韦庄 | 老熟妇乱子伦视频在线观看| 变态另类成人亚洲欧美熟女| 一个人看的www免费观看视频| 99视频精品全部免费 在线| 亚洲成人久久性| 波多野结衣高清无吗| 国内揄拍国产精品人妻在线| 久久精品夜夜夜夜夜久久蜜豆| 国产成年人精品一区二区| АⅤ资源中文在线天堂| 日日干狠狠操夜夜爽| 免费观看的影片在线观看| 一进一出抽搐动态| 免费观看的影片在线观看| 麻豆成人午夜福利视频| 熟女人妻精品中文字幕| 精品一区二区三区视频在线| 亚洲国产精品成人综合色| 能在线免费看毛片的网站| 男插女下体视频免费在线播放| 日韩一区二区三区影片| 麻豆成人午夜福利视频| 亚洲婷婷狠狠爱综合网| 免费观看人在逋| 狂野欧美激情性xxxx在线观看| 亚洲国产欧美在线一区| 九九久久精品国产亚洲av麻豆| 久久人人精品亚洲av| www日本黄色视频网| 黄片无遮挡物在线观看| 97人妻精品一区二区三区麻豆| 日韩中字成人| 久久久久久久久大av| 噜噜噜噜噜久久久久久91| 国产精品三级大全| 少妇被粗大猛烈的视频| 国产三级在线视频| 天天躁夜夜躁狠狠久久av| 国产熟女欧美一区二区| 不卡视频在线观看欧美| av福利片在线观看| 日本黄大片高清| 只有这里有精品99| 久久精品综合一区二区三区| 免费观看人在逋| 久久韩国三级中文字幕| 国产亚洲精品久久久com| 国产精品伦人一区二区| 直男gayav资源| 午夜久久久久精精品| 欧美最黄视频在线播放免费| 91精品一卡2卡3卡4卡| 日本熟妇午夜| 你懂的网址亚洲精品在线观看 | 欧美变态另类bdsm刘玥| 欧美高清成人免费视频www| 日韩国内少妇激情av| 国产免费男女视频| 少妇熟女aⅴ在线视频| 亚洲在久久综合| 国语自产精品视频在线第100页| 亚洲国产精品合色在线| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧美日韩东京热| 日本成人三级电影网站| 高清午夜精品一区二区三区 | 欧美性感艳星| 欧美性猛交╳xxx乱大交人| 午夜亚洲福利在线播放| 精品久久久久久久末码| av黄色大香蕉| 天堂网av新在线| 日韩一区二区视频免费看| 亚洲丝袜综合中文字幕| 亚洲国产日韩欧美精品在线观看| 91狼人影院| 国产白丝娇喘喷水9色精品| 一个人看视频在线观看www免费| 久久鲁丝午夜福利片| 国产毛片a区久久久久| 夜夜看夜夜爽夜夜摸| 1000部很黄的大片| 欧美日韩一区二区视频在线观看视频在线 | 寂寞人妻少妇视频99o| 在线观看av片永久免费下载| 亚洲真实伦在线观看| av国产免费在线观看| 日韩制服骚丝袜av| 国产视频首页在线观看| 六月丁香七月| 亚洲欧美中文字幕日韩二区| 性欧美人与动物交配| 精品久久久久久久人妻蜜臀av| 亚洲色图av天堂| 亚洲激情五月婷婷啪啪| 国产高清激情床上av| 少妇裸体淫交视频免费看高清| 欧美色视频一区免费| 高清毛片免费看| 欧美最黄视频在线播放免费| 97超碰精品成人国产| 日韩,欧美,国产一区二区三区 | 亚洲成av人片在线播放无| 中文资源天堂在线| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩无卡精品| 成年女人看的毛片在线观看| 日韩制服骚丝袜av| 欧美日韩乱码在线| 人妻少妇偷人精品九色| 最近最新中文字幕大全电影3| 日韩大尺度精品在线看网址| .国产精品久久| 国产精品综合久久久久久久免费| 黄片wwwwww| 久久久精品欧美日韩精品| 一个人观看的视频www高清免费观看| 日韩亚洲欧美综合| 极品教师在线视频| 黄色欧美视频在线观看| 一个人看视频在线观看www免费| 亚洲一区二区三区色噜噜| 村上凉子中文字幕在线| 91麻豆精品激情在线观看国产| or卡值多少钱| 能在线免费看毛片的网站| 国产中年淑女户外野战色| 国内久久婷婷六月综合欲色啪| 久久这里有精品视频免费| 国产探花极品一区二区| 两个人的视频大全免费| 日韩欧美精品v在线| 国产成人91sexporn| 国产一区二区激情短视频| 狠狠狠狠99中文字幕| 国产美女午夜福利| 国产久久久一区二区三区| av国产免费在线观看| 黑人高潮一二区| 91久久精品国产一区二区三区| 午夜福利在线观看吧| 少妇人妻精品综合一区二区 | 中文字幕av在线有码专区| 久久久久国产网址| 国产精品一区二区性色av| 男女啪啪激烈高潮av片| 国产成人91sexporn| 久久久精品欧美日韩精品| 国产精品久久视频播放| 成人二区视频| 精品久久久噜噜| 免费看光身美女| 亚洲av中文av极速乱| 97超视频在线观看视频| 69人妻影院| 给我免费播放毛片高清在线观看| 夜夜夜夜夜久久久久| 欧美潮喷喷水| 少妇高潮的动态图| 亚洲电影在线观看av| avwww免费| 国产成人福利小说| 中国美女看黄片| 麻豆精品久久久久久蜜桃| 乱码一卡2卡4卡精品| 中国美白少妇内射xxxbb| 欧美区成人在线视频| 日韩欧美精品v在线| 久久久精品94久久精品| 亚洲精品国产av成人精品| 成人鲁丝片一二三区免费| 久久久精品大字幕| 成人高潮视频无遮挡免费网站| 在线观看av片永久免费下载| 久久精品国产亚洲av香蕉五月| 欧美一区二区亚洲| a级毛色黄片| 国产午夜精品一二区理论片| 中文亚洲av片在线观看爽| 色噜噜av男人的天堂激情| 91久久精品国产一区二区三区| av女优亚洲男人天堂| 亚洲天堂国产精品一区在线| 中文亚洲av片在线观看爽| 欧美精品国产亚洲| 亚洲欧美中文字幕日韩二区| 99热精品在线国产| 99久久久亚洲精品蜜臀av| 久久精品综合一区二区三区| 国产69精品久久久久777片| 两性午夜刺激爽爽歪歪视频在线观看| 精品人妻偷拍中文字幕| 精品午夜福利在线看| 最近2019中文字幕mv第一页| 色5月婷婷丁香| 97人妻精品一区二区三区麻豆| 毛片女人毛片| 欧美日本视频| 精品国产三级普通话版| 午夜精品在线福利| 中文字幕av在线有码专区| 国产一区二区在线观看日韩| 给我免费播放毛片高清在线观看| 赤兔流量卡办理| 欧美日韩精品成人综合77777| 久久午夜亚洲精品久久| 最近最新中文字幕大全电影3| 成人毛片60女人毛片免费| 舔av片在线| 麻豆av噜噜一区二区三区| 小说图片视频综合网站| 亚洲国产欧美在线一区| 少妇熟女欧美另类| 变态另类成人亚洲欧美熟女| 床上黄色一级片| 欧美极品一区二区三区四区| 12—13女人毛片做爰片一| 亚洲七黄色美女视频| 国产高清三级在线| 国产精品福利在线免费观看| 美女cb高潮喷水在线观看| 久久久久国产网址| 色哟哟·www| 人妻久久中文字幕网| 麻豆久久精品国产亚洲av| 欧美性猛交黑人性爽| 高清在线视频一区二区三区 | 级片在线观看| 欧美日韩综合久久久久久| 国产黄a三级三级三级人| 狂野欧美白嫩少妇大欣赏| 久久鲁丝午夜福利片| 丰满人妻一区二区三区视频av| 桃色一区二区三区在线观看| 高清午夜精品一区二区三区 | av福利片在线观看| 哪里可以看免费的av片| 99热这里只有精品一区| 日韩人妻高清精品专区| 一个人免费在线观看电影| 中文字幕久久专区| 日韩欧美一区二区三区在线观看| 性欧美人与动物交配| 国产一级毛片七仙女欲春2| 嫩草影院精品99| 亚洲图色成人| 三级国产精品欧美在线观看| 岛国在线免费视频观看| 日韩视频在线欧美| 久久精品久久久久久噜噜老黄 | 国产成人aa在线观看| 国语自产精品视频在线第100页| 色综合站精品国产| 亚洲第一电影网av| 91午夜精品亚洲一区二区三区| 少妇猛男粗大的猛烈进出视频 | 22中文网久久字幕| 日日摸夜夜添夜夜爱| 欧美潮喷喷水| 搡女人真爽免费视频火全软件| a级毛色黄片| 日韩欧美 国产精品| 18禁在线无遮挡免费观看视频| 婷婷亚洲欧美| 亚洲五月天丁香| .国产精品久久| 欧美成人一区二区免费高清观看| 亚洲精品成人久久久久久| 又爽又黄无遮挡网站| 久久99蜜桃精品久久| 99热这里只有是精品在线观看| 自拍偷自拍亚洲精品老妇| 人妻系列 视频| 人人妻人人澡欧美一区二区| 男人狂女人下面高潮的视频| 亚洲欧洲国产日韩| 中文字幕人妻熟人妻熟丝袜美| 免费电影在线观看免费观看| 久久久精品大字幕| 久久久精品94久久精品| 亚洲性久久影院| 精品免费久久久久久久清纯| 久久久久久久久久久免费av| 午夜老司机福利剧场| 高清午夜精品一区二区三区 | 黄片无遮挡物在线观看| 中文字幕免费在线视频6| 亚洲美女视频黄频| 亚洲人成网站高清观看| 色哟哟·www| 久久久久久久久中文| 一级二级三级毛片免费看| 国产亚洲精品av在线| 亚洲精品粉嫩美女一区| av在线亚洲专区| 国产高清有码在线观看视频| 日本成人三级电影网站| 麻豆成人午夜福利视频| 只有这里有精品99| 亚洲人成网站在线播| 国产伦精品一区二区三区视频9| 日韩欧美国产在线观看| 国产亚洲91精品色在线| 人人妻人人看人人澡| 欧美性感艳星| 看十八女毛片水多多多| 亚洲最大成人手机在线| 边亲边吃奶的免费视频| 麻豆乱淫一区二区| 国产一区二区在线观看日韩| 观看免费一级毛片| 丰满乱子伦码专区| 久久精品国产亚洲av涩爱 | 成人毛片60女人毛片免费| 欧美xxxx性猛交bbbb| 日产精品乱码卡一卡2卡三| 男人狂女人下面高潮的视频| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲,欧美,日韩| 99国产精品一区二区蜜桃av| 在线观看免费视频日本深夜| 国产人妻一区二区三区在| 午夜福利在线在线| 欧美成人免费av一区二区三区| 婷婷六月久久综合丁香| 日韩大尺度精品在线看网址| 欧美性猛交╳xxx乱大交人| 国内精品一区二区在线观看| 亚洲欧美清纯卡通| 日韩一本色道免费dvd| 麻豆av噜噜一区二区三区| 久久精品久久久久久噜噜老黄 | 亚洲一区二区三区色噜噜| 热99re8久久精品国产| 国产 一区精品| 国产精品永久免费网站| 国产亚洲91精品色在线| 欧美日韩国产亚洲二区| 久久午夜福利片| 麻豆久久精品国产亚洲av| 精品无人区乱码1区二区| 精品人妻视频免费看| 内地一区二区视频在线| 国产精品一及| 免费一级毛片在线播放高清视频| 免费大片18禁| 三级男女做爰猛烈吃奶摸视频| 1024手机看黄色片| 美女xxoo啪啪120秒动态图| 午夜免费男女啪啪视频观看| 两个人视频免费观看高清| 人体艺术视频欧美日本| 热99re8久久精品国产| av国产免费在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲内射少妇av| 成人永久免费在线观看视频| 午夜老司机福利剧场| 日本免费一区二区三区高清不卡| 99久国产av精品| 我的女老师完整版在线观看| 日韩欧美在线乱码| 卡戴珊不雅视频在线播放| 国产极品精品免费视频能看的| 能在线免费观看的黄片| 日韩中字成人| 欧美一级a爱片免费观看看| 黄色视频,在线免费观看| 网址你懂的国产日韩在线| 国产亚洲精品久久久com| 日日啪夜夜撸| 成年av动漫网址| 国产精品一区二区三区四区久久| 国产伦一二天堂av在线观看| 国模一区二区三区四区视频| 大型黄色视频在线免费观看| 午夜福利成人在线免费观看| 国产亚洲91精品色在线| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利视频1000在线观看| 免费看a级黄色片| 99热这里只有精品一区| 国产黄a三级三级三级人| 两性午夜刺激爽爽歪歪视频在线观看| 99久久精品国产国产毛片| 日本五十路高清| 直男gayav资源| 久久久久久国产a免费观看| 熟妇人妻久久中文字幕3abv| 午夜福利视频1000在线观看| 最好的美女福利视频网| 亚洲第一电影网av| 国产视频首页在线观看| 在线观看66精品国产| 国产黄色视频一区二区在线观看 | а√天堂www在线а√下载| 久久精品91蜜桃| 成年免费大片在线观看| 美女黄网站色视频| 亚洲国产日韩欧美精品在线观看| 欧美xxxx黑人xx丫x性爽| 男女啪啪激烈高潮av片| 久久精品国产亚洲av天美| 男女边吃奶边做爰视频| 白带黄色成豆腐渣| 一进一出抽搐动态| 久久这里只有精品中国| 久久久久久久久中文| 六月丁香七月| 久久欧美精品欧美久久欧美| 久久精品国产亚洲av天美| 色综合站精品国产| 免费av不卡在线播放| 在线免费观看不下载黄p国产| 成人国产麻豆网| 国产激情偷乱视频一区二区| 卡戴珊不雅视频在线播放| 久久热精品热| 蜜桃久久精品国产亚洲av| 国国产精品蜜臀av免费| 免费av不卡在线播放| 三级毛片av免费| 久久亚洲精品不卡| 久久久久久久久中文| 一区二区三区免费毛片| 此物有八面人人有两片| 日本三级黄在线观看| 久久久久久大精品| 欧美区成人在线视频| 国产精品一二三区在线看| 久久久久久久久大av| 日日摸夜夜添夜夜添av毛片| 亚洲va在线va天堂va国产| 国产三级在线视频| 亚洲在线自拍视频| 欧美zozozo另类| 在线观看av片永久免费下载| 日本在线视频免费播放| 狂野欧美激情性xxxx在线观看| 欧美日韩综合久久久久久| 国产精品一二三区在线看| 99热只有精品国产| 日本熟妇午夜| 成人性生交大片免费视频hd| 国内揄拍国产精品人妻在线| 在线国产一区二区在线| 亚洲欧美清纯卡通| av免费观看日本| 日韩大尺度精品在线看网址| 卡戴珊不雅视频在线播放| 国产三级中文精品| 久久人妻av系列| 精品免费久久久久久久清纯| 成人毛片a级毛片在线播放| 亚洲美女搞黄在线观看| 舔av片在线| 亚洲国产高清在线一区二区三| 国产一级毛片在线| 亚洲欧洲日产国产| 欧美日韩综合久久久久久| 成人鲁丝片一二三区免费| 久久精品影院6| 国产精品乱码一区二三区的特点| 久久午夜福利片| 一进一出抽搐gif免费好疼| 亚洲成a人片在线一区二区| 最新中文字幕久久久久| 蜜桃久久精品国产亚洲av| 日韩在线高清观看一区二区三区| 中文亚洲av片在线观看爽| 精品无人区乱码1区二区| 婷婷六月久久综合丁香| 人人妻人人澡欧美一区二区| 啦啦啦啦在线视频资源| 热99在线观看视频| 美女被艹到高潮喷水动态| 我要搜黄色片| av福利片在线观看| 日韩 亚洲 欧美在线| 美女大奶头视频| 十八禁国产超污无遮挡网站| 青春草亚洲视频在线观看| 在线观看66精品国产| 国产成人aa在线观看| 少妇丰满av| 十八禁国产超污无遮挡网站| a级毛片a级免费在线| 免费看光身美女| 一本久久中文字幕| 亚洲精品日韩在线中文字幕 | 久久99蜜桃精品久久| 老师上课跳d突然被开到最大视频| a级毛片a级免费在线| 欧美最新免费一区二区三区| 国产精品美女特级片免费视频播放器| 欧美日韩精品成人综合77777| 人妻久久中文字幕网| 亚洲国产欧洲综合997久久,| 久久精品国产亚洲av涩爱 | 91久久精品电影网| 老司机福利观看| 亚洲精品粉嫩美女一区| 白带黄色成豆腐渣|