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

    基于擾動(dòng)指數(shù)的岷江上游森林?jǐn)_動(dòng)時(shí)間序列研究
    ——以理縣為例

    2017-08-09 12:17:22張宇欣楊存建羅銀建馬洋洋
    浙江林業(yè)科技 2017年3期
    關(guān)鍵詞:理縣擾動(dòng)植被

    張宇欣,倪 靜,楊存建,羅銀建,馬洋洋

    (四川師范大學(xué) 西南土地資源評(píng)價(jià)與監(jiān)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室,四川 成都 610068)

    基于擾動(dòng)指數(shù)的岷江上游森林?jǐn)_動(dòng)時(shí)間序列研究
    ——以理縣為例

    張宇欣,倪 靜,楊存建,羅銀建,馬洋洋

    (四川師范大學(xué) 西南土地資源評(píng)價(jià)與監(jiān)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室,四川 成都 610068)

    利用擾動(dòng)指數(shù)算法對(duì)四川省阿壩州理縣的森林?jǐn)_動(dòng)進(jìn)行時(shí)間序列監(jiān)測(cè)研究?;赥M數(shù)據(jù)對(duì)1994-2001,2001-2007,2007-2008,2008-2011年4個(gè)時(shí)間段的理縣森林?jǐn)_動(dòng)進(jìn)行監(jiān)測(cè)。結(jié)果表明,4個(gè)時(shí)間段理縣受擾動(dòng)的森林面積分別為437.94,278.46,260.46,184.14 hm2,主要原因依次為人為活動(dòng),洪水及次生災(zāi)害影響,汶川地震破壞,滑坡、塌方和堰塞湖等地震次生災(zāi)害。經(jīng)過(guò)驗(yàn)證kappa系數(shù)達(dá)到0.84。數(shù)據(jù)顯示,2001年前該縣森林?jǐn)_動(dòng)較為嚴(yán)重,2001年后該縣受擾動(dòng)森林面積逐步減少,說(shuō)明森林狀況逐漸轉(zhuǎn)好。

    理縣;森林?jǐn)_動(dòng);擾動(dòng)指數(shù);時(shí)間序列;監(jiān)測(cè)

    森林是地表土壤、水體與大氣之間進(jìn)行物質(zhì)、能量交換的關(guān)鍵要素,在地表能量循環(huán)、水循環(huán)和生物化學(xué)循環(huán)中發(fā)揮著重要的作用[1]。森林?jǐn)_動(dòng)是區(qū)域生態(tài)系統(tǒng)環(huán)境變化的重要指示,更是全球生態(tài)環(huán)境的研究熱點(diǎn),擾動(dòng)類型主要為森林采伐和自然災(zāi)害。國(guó)外學(xué)者對(duì)于森林?jǐn)_動(dòng)研究較多[2]。Masek J G等[3]采用Landsat數(shù)據(jù)檢測(cè)北美森林的擾動(dòng)。Wylie B等[4]研究了氣候變化對(duì)于生態(tài)系統(tǒng)的干擾和影響及其在加拿大育空河流域北部森林的表現(xiàn)。Frantz D等[5]采用Landsat和Modis數(shù)據(jù)對(duì)昆士蘭的森林?jǐn)_動(dòng)進(jìn)行了時(shí)間序列的研究。Misurec J等[6]采用DI算法研究挪威云杉的森林?jǐn)_動(dòng)時(shí)空變化。楊存建等[7]利用GIS和RS技術(shù)揭示出四川西部2002-2008年植被變化的時(shí)空特征。對(duì)四川省阿壩藏族羌族自治州理縣進(jìn)行了區(qū)域尺度的森林?jǐn)_動(dòng)研究。

    1 材料與方法

    1.1 研究區(qū)概況

    理縣隸屬于四川省阿壩藏族羌族自治州,102°32′46″ ~ 103°30′30″ E,30°54′43″ ~ 31°12′12″ N。全縣面積為4 318 km2。地質(zhì)結(jié)構(gòu)屬龍門山斷裂帶中斷,地形呈蜿蜒起伏的立體單元,地表由西北向東南傾斜。地貌類型為低中山-中山-高山-極高山,是典型的中高山峽谷區(qū),海拔1 422 ~ 5 922 m,平均海拔2 700 m,植被垂直帶譜明顯。氣候受西伯利亞西風(fēng)氣流、印度洋暖流和太平洋東南季風(fēng)三個(gè)環(huán)流的影響,形成季風(fēng)氣候,有干濕季之分。因海拔高差懸殊,地形復(fù)雜,氣候差異顯著,具有山地立體型氣候特征。冬季受青藏高原的北方冷氣流的影響,降水稀少,日照強(qiáng)烈,晴朗多大風(fēng),五、六月份西南季風(fēng)加強(qiáng),溫暖濕潤(rùn),降水增多,形成雨季。七、八月份青藏高壓穩(wěn)定,副熱帶高壓西伸,降水減少,形成伏旱。九、十月份雨量增加,形成低溫降雨季節(jié)。年降水量650 ~ 1 000 mm,為植被的生長(zhǎng)發(fā)育提供有利的雨水條件。

    理縣位于青藏高原東部與四川盆地的過(guò)渡地帶,是岷江上游生態(tài)屏障的重要組成。區(qū)域內(nèi)森林資源豐富,研究時(shí)序內(nèi)受到擾動(dòng)變化較大,1994-2000年受人為采伐影響,自2000年開展退耕還林工程,2008年受汶川地震破壞。研究其森林?jǐn)_動(dòng)對(duì)維護(hù)和提升該縣生態(tài)系統(tǒng)服務(wù)功能具有重要意義,是提高林地監(jiān)管能力、開展生態(tài)空間規(guī)劃、兌現(xiàn)生態(tài)效益補(bǔ)償,加強(qiáng)林地保護(hù)利用管理和科學(xué)預(yù)測(cè)植被長(zhǎng)勢(shì)的重要基礎(chǔ)和科學(xué)支撐,為該縣區(qū)域經(jīng)濟(jì)社會(huì)發(fā)展、城鄉(xiāng)建設(shè)、土地利用及生態(tài)環(huán)境保護(hù)等奠定基礎(chǔ)。研究區(qū)如圖1所示。

    1.2 研究數(shù)據(jù)

    1.2.1 遙感影像及預(yù)處理 共選用5期TM影像Path 130/Row 38(WRS-2),獲取日期分別為1994年6月26日,2001年6月13日,2007年9月18日,2008年10月6日和2011年8月12日,云量皆小于10%,數(shù)據(jù)級(jí)別皆為L(zhǎng)1。數(shù)據(jù)均來(lái)源于美國(guó)地質(zhì)勘探局(USGS)并且均選取植被生長(zhǎng)季節(jié)內(nèi)的影像(6-10月上旬)。通過(guò)ENVI對(duì)遙感數(shù)據(jù)進(jìn)行輻射定標(biāo)、大氣校正后再對(duì)研究區(qū)內(nèi)的云及其陰影進(jìn)行目視識(shí)別后建立掩膜。

    1.2.2 地面調(diào)查資料 調(diào)查資料包括矢量數(shù)據(jù)和柵格數(shù)據(jù)。矢量數(shù)據(jù)為研究區(qū)2000年、2010年土地利用數(shù)據(jù)。通過(guò)土地利用數(shù)據(jù)提取森林植被的覆蓋范圍,再對(duì)區(qū)域內(nèi)的影像進(jìn)行森林植被的掩膜提取;柵格數(shù)據(jù)為研究區(qū)30 mDEM數(shù)據(jù)。

    1.3 研究方法

    1.3.1 森林?jǐn)_動(dòng)指數(shù)時(shí)間序列分析 纓帽變換(Tasseled Cap Transformation,TCT)將Landsat的6個(gè)原始波段(1 ~ 5和7)折合成亮度、綠度和濕度3個(gè)地表特征的分量[8]。擾動(dòng)指數(shù)(DI)是一種基于Landsat TM/ETM+的纓帽變換3個(gè)地表特征分量的指數(shù)[9]。具體公式見公式1。

    圖1 研究區(qū)域Figure 1 Map of the study area

    式中,B’,G’和W’代表TCT中亮度、綠度和濕度的標(biāo)準(zhǔn)化指數(shù)。其標(biāo)準(zhǔn)化指數(shù)的計(jì)算公式見公式2。

    式中,μB表示森林平均亮度;σB表示森林亮度的標(biāo)準(zhǔn)差。

    經(jīng)過(guò)計(jì)算,受到劇烈擾動(dòng)的森林將表現(xiàn)出較高的B’、較低的G’和較低的W’值,此時(shí)的DI值較高。相對(duì)的,未擾動(dòng)和輕微擾動(dòng)區(qū)域則表現(xiàn)出較低的B’,較高的G’和較高的W’值,此時(shí)的DI值較低[3]。

    1.3.2 降噪處理 考慮到纓帽變換后的影像中仍可能存在噪聲等干擾信息,選用3×3窗口作為模版對(duì)結(jié)果進(jìn)行離散化窗口滑窗卷積,使中間點(diǎn)的DI值等于八連通區(qū)域的加權(quán)DI均值,其具體公式表示為:

    式中,x和h分別為兩個(gè)序列,i與n皆為自變量。

    1.3.3 擾動(dòng)區(qū)域提取 擾動(dòng)區(qū)域通過(guò)構(gòu)建年間差值指數(shù)ΔDI進(jìn)行研究,經(jīng)過(guò)多次試驗(yàn),確定合適的閾值進(jìn)行擾動(dòng)區(qū)域提取,具體公式見公式4。

    式中,n為后一年份,l為前一年份。

    2 結(jié)果與分析

    2.1 精度驗(yàn)證

    驗(yàn)證資料:(1)ΔNDVI:通過(guò)5期影像計(jì)算NDVI,構(gòu)建差值指數(shù)計(jì)算NDVI;(2)2000,2010年理縣土地利用數(shù)據(jù);(3)高分遙感影像(Google E arth)進(jìn)行輔助驗(yàn)證。驗(yàn)證方法:本小節(jié)采用分層抽樣的方法,通過(guò)對(duì)各區(qū)間擾動(dòng)區(qū)域及未擾動(dòng)區(qū)域隨機(jī)抽取30個(gè)采樣點(diǎn),總計(jì)240個(gè)采樣點(diǎn)。利用混淆矩陣建立精度評(píng)價(jià)表,最終計(jì)算kappa系數(shù)為0.84,詳見表1。

    表1 混淆矩陣Table 1 Confusion matrix

    2.2 結(jié)果及影響因子分析

    2.2.1 趨勢(shì)線分析 通過(guò)趨勢(shì)線分析法模擬1994-2011年理縣森林?jǐn)_動(dòng)面積變化趨勢(shì),其計(jì)算公式為:

    從圖2可以看出,森林受擾動(dòng)區(qū)域面積呈逐步下降趨勢(shì),說(shuō)明森林結(jié)構(gòu)不斷轉(zhuǎn)好,森林保護(hù)意識(shí)及力度逐步加強(qiáng)。

    圖2 1994-2011年森林?jǐn)_動(dòng)面積Figure 3 Forest disturbance area from 1994 to 2011

    2.2.2 影響因子分析 依據(jù)自然間斷點(diǎn)分級(jí)法,將研究區(qū)的30 m DEM數(shù)據(jù)進(jìn)行分級(jí),將高程1 396 ~ 2 768 m,2 768 ~ 3 520 m,3 520 ~ 4 173 m,4 173 ~ 5 854 m分別劃分為1 ~ 4級(jí)(圖3)。

    圖3 研究區(qū)DEM(a)與重分類(b)Figure 4 DEM in study area and reclassification

    在Arcgis10.2軟件中,利用面積制表功能分時(shí)間段、分高程統(tǒng)計(jì)落在擾動(dòng)區(qū)域在不同海拔等級(jí)中的分布情況。詳見表2。

    由表2可知,1994-2011年的擾動(dòng)總面積達(dá)1 160.91 hm2。1994-2011,2001-2007,2007-2008,2008-2011年的森林?jǐn)_動(dòng)區(qū)域分別主要分布在第三、第四等級(jí),分別占總擾動(dòng)面積的35.51%,22.96%,7.60%,6.85%。由此可知,較高海拔即3 520 ~ 5 854 m范圍內(nèi)的森林較易受到擾動(dòng)。

    表2 擾動(dòng)面積海拔分布Table 2 Disturbance area and their distribution of elevation

    DEM生成研究區(qū)坡度和坡向數(shù)據(jù)。依據(jù)《中國(guó)土壤侵蝕分類標(biāo)準(zhǔn)(SL190)》中坡度的劃分標(biāo)準(zhǔn),將坡度進(jìn)行重新分級(jí)。具體劃分標(biāo)準(zhǔn)見表3。

    在Arcgis10.2軟件中,利用面積制表功能分時(shí)間段、分坡度統(tǒng)計(jì)落在擾動(dòng)區(qū)域在不同坡度等級(jí)中的分布情況。詳見表4。

    由表4可知,1994-2011,2001-2007,2007-2008,2008-2011年的森林?jǐn)_動(dòng)區(qū)域分別主要分布在緩的、陡的、極陡的三個(gè)坡度級(jí)上,擾動(dòng)面積分別為148.05,91.89,119.79,64.71 hm2,占總擾動(dòng)面積的12.75%,7.92%,10.32%,5.57%。由此可知較高的坡度即15°以上的森林植被更加容易受到擾動(dòng)影響。

    依據(jù)全國(guó)數(shù)字地貌制圖中的坡向劃分標(biāo)準(zhǔn),對(duì)坡向進(jìn)行重新分級(jí)。劃分結(jié)果見表5。

    表3 坡度劃分標(biāo)準(zhǔn)Table 3 Grading of slope

    表4 擾動(dòng)面積坡度分布Table 4 Disturbance area and their distribution at slope

    在Arcgis10.2軟件中,利用面積制表功能分時(shí)間段、分坡度統(tǒng)計(jì)落在擾動(dòng)區(qū)域在不同坡向中的分布情況,見表6。

    由表6可知,1994-2011,2001-2007,2007-2008,2008-2011年的森林?jǐn)_動(dòng)區(qū)域分別主要分布在NE,S,NW三個(gè)坡向上,擾動(dòng)面積分別為100.44,89.37,55.98,47.79 hm2,占總擾動(dòng)面積的比例分別為8.65%,7.70%,4.82%,4.12%。由此可知在22.5 ~ 67.5°,157.5 ~ 202.5°,292.5 ~ 337.5°三個(gè)范圍內(nèi)森林植被較易受到擾動(dòng)。

    表5 研究區(qū)坡向劃分標(biāo)準(zhǔn)Table 5 Classification of aspect in study area

    表6 擾動(dòng)面積坡向分布Table 6 Disturbance area and their distribution at aspect

    3 結(jié)論

    綜上所述,(1)1994-2001年理縣受擾動(dòng)的森林面積為437.94 hm2,分布區(qū)域主要為森林邊緣、鄰近建筑用地,人為活動(dòng)是這期間森林?jǐn)_動(dòng)的主要原因,開墾、亂砍亂伐等對(duì)森林造成破壞;(2)2001-2007年森林?jǐn)_動(dòng)面積為278.46 hm2,擾動(dòng)區(qū)域主要為海拔較高的山頂附近。由于山高坡陡,地表土壤松散,且研究時(shí)期為多雨季節(jié),擾動(dòng)主要原因?yàn)楹樗按紊鸀?zāi)害影響[10];(3)2007-2008年森林受擾動(dòng)面積為260.46 hm2,主要分布于山坡、山谷、河谷附近,由于理縣為受汶川地震影響最嚴(yán)重的區(qū)域之一[11-12],因此森林區(qū)域受地震、滑坡、泥石流等災(zāi)害破壞嚴(yán)重;(4)2008-2011年森林?jǐn)_動(dòng)面積為184.14 hm2,滑坡、塌方和堰塞湖等地震次生災(zāi)害是這一階段導(dǎo)致森林受到破壞的主要原因。

    從擾動(dòng)因子分析:(1)2001年前森林受擾動(dòng)影響較為劇烈,主要原因是人為采伐;(2)2001年后森林?jǐn)_動(dòng)面積逐漸減少,說(shuō)明森林狀況逐漸轉(zhuǎn)好,森林保護(hù)意識(shí)不斷提高;(3)2008年由于受到汶川地震影響,導(dǎo)致森林?jǐn)_動(dòng)區(qū)域再次升高。

    1994-2011年理縣森林受擾動(dòng)區(qū)域面積呈逐步下降趨勢(shì),說(shuō)明該縣森林狀況不斷改善,證明該區(qū)域正逐步加強(qiáng)林地保護(hù)和林地監(jiān)管能力。研究采用DI擾動(dòng)指數(shù)是基于TM/ETM+纓帽變換的亮度、綠度和濕度分量的算法;研究中對(duì)于閾值的選取經(jīng)過(guò)多次試驗(yàn)后選取最佳閾值進(jìn)行提取,因此閾值的選取對(duì)于研究結(jié)果有一定的影響。

    [1] 陳效逑,王恒. 1982-2003年內(nèi)蒙古植被帶和植被覆蓋度的時(shí)空變化[J]. 地理學(xué)報(bào),2009,64(1):84-94.

    [2] 楊辰,沈潤(rùn)平. 森林?jǐn)_動(dòng)遙感監(jiān)測(cè)研究進(jìn)展[J]. 國(guó)土資源遙感,2015,27(1):1-8.

    [3] Masek J G,Huang C,Wolfe R,et al. North American forest disturbance mapped from a decadal Landsat record[J]. Remote Sens Environ,2008,112(6):2914-2926.

    [4] Wylie B,Rigge M,Brisco B,et al. Effects of Disturbance and Climate Change on Ecosystem Performance in the Yukon River Basin Boreal Forest[J]. Remote Sens,2014,6(10):9145-9169.

    [5] Frantz D,R?der A,Udelhoven T,et a l. Forest Disturbance M apping U sing D ense S ynthetic Landsat/MODIS Time-Series a nd Permutation-Based Disturbance Index Detection[J]. Remote Sens,2016,8(4):277.

    [6] Misurec J,Kopa?ková V,Lhotáková Z,et al. Detection of spatio-temporal changes of Norway spruce forest stands in Ore Mountains using airborne hyperspectral imagery[C]. International Electronic Conference on Remote Sensing. 2016:1-20.

    [7] 楊存建,趙梓健,任小蘭,等.基于遙感和GIS的川西綠被時(shí)空變化研究[J]. 生態(tài)學(xué)報(bào),2012,32(2):632-640.

    [8] Crist E P. A TM Tasseled Cap equivalent transformation for reflectance factor data[J]. Remote Sens Environ,1985,17(3):301-306.

    [9] Healey S P,Cohen W B,Yang Z,et al. Comparison of Tasseled Cap-based Landsat data structures for use in forest disturbance detection[J]. Remote Sens Environ,2005,97(3):301-310.

    [10] 郭兵,陶和平,劉斌濤,等. 基于GIS和USLE的汶川地震后理縣土壤侵蝕特征及分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(14):118-126.

    [11] 甘建軍,吳晗,唐春,等. 四川理縣塔斯溝多發(fā)性特大泥石流成災(zāi)機(jī)理與特征[J]. 災(zāi)害學(xué),2015(4):59-63.

    [12] 趙旦,張淼,于名召,等. 汶川地震災(zāi)后農(nóng)田和森林植被恢復(fù)遙感監(jiān)測(cè)[J]. 遙感學(xué)報(bào),2014,8(4):958-970.

    Study on Spatio-temporal Variation of Forest Disturbance in Lixian by Disturbance Index

    ZHANG Yu-xin,NI Jing,YANG Cun-jian,LUO Yin-jian,MA Yang-yang
    (Key Laboratory of the Evaluation and Monitoring of Southwest Land Resources of Ministry of Education, Sichuan Normal University, Chengdu 610068, China)

    Research was made on forest disturbance in time series in Lixian county, Sichuan province by disturbance index. Based on TM remote sensing data, determinations were conducted on forest disturbance during four periods, like 1994-2001, 2001-2007, 2007- 2008 and 2008-2011. The result showed that during the 4 periods, forest disturbance area in Lixian was 437.94ha, 278.46ha, 260.46ha and184.14ha. The main causes were human activities, flood and secondary disasters, Wenchuan Earthquake and its secondary disasters of the earthquake such as landslide and dammed lake. Kappa coefficient was 0.84. Analysis demonstrated that forest disturbance was serious before 2001, but decreased then, indicating healthier forest.

    Lixian county; forest disturbance; disturbance index; time series; monitoring

    S757.2;TP73

    :A

    :1001-3776(2017)03-0048-06

    10.3969/j.issn.1001-3776.2017.03.009

    2016-12-09 ;

    2017-03-17

    四川師范大學(xué)研究生優(yōu)秀學(xué)位論文培育基金資助(校研字(2016)4-42)

    張宇欣,碩士研究生,從事遙感與GIS理論技術(shù)及應(yīng)用研究;E-mail:zhangyuxin11223@163.com。通信作者:楊存建,教授,從事遙感和地理信息系統(tǒng)應(yīng)用研究;E-mail:yangcj2008@126.com。

    猜你喜歡
    理縣擾動(dòng)植被
    基于植被復(fù)綠技術(shù)的孔植試驗(yàn)及應(yīng)用
    Bernoulli泛函上典則酉對(duì)合的擾動(dòng)
    基于“雙評(píng)價(jià)”的水土流失敏感性評(píng)價(jià)——以四川省理縣為例
    (h)性質(zhì)及其擾動(dòng)
    羌族歌舞探究——以四川理縣蒲溪鄉(xiāng)為例
    理縣蒲溪村羌族釋比戲文化初探
    綠色植被在溯溪旅游中的應(yīng)用
    小噪聲擾動(dòng)的二維擴(kuò)散的極大似然估計(jì)
    基于原生植被的長(zhǎng)山群島植被退化分析
    用于光伏MPPT中的模糊控制占空比擾動(dòng)法
    最近最新中文字幕大全免费视频| 啦啦啦 在线观看视频| 久久久久国内视频| а√天堂www在线а√下载| 嫩草影视91久久| 亚洲成人精品中文字幕电影| 色综合亚洲欧美另类图片| 久久人人97超碰香蕉20202| 麻豆成人av在线观看| 9191精品国产免费久久| 欧美成人午夜精品| 久久国产精品男人的天堂亚洲| 亚洲欧美日韩另类电影网站| 999久久久精品免费观看国产| 欧美一区二区精品小视频在线| 日韩欧美一区视频在线观看| 大型av网站在线播放| 久久国产精品影院| 一本综合久久免费| 在线观看www视频免费| 日韩精品中文字幕看吧| 亚洲精品国产色婷婷电影| 久久精品国产亚洲av香蕉五月| 日本撒尿小便嘘嘘汇集6| www.自偷自拍.com| 午夜福利在线观看吧| 淫妇啪啪啪对白视频| 亚洲情色 制服丝袜| 亚洲精品一卡2卡三卡4卡5卡| 亚洲久久久国产精品| aaaaa片日本免费| 一边摸一边抽搐一进一出视频| 国产欧美日韩精品亚洲av| 亚洲成a人片在线一区二区| 欧美精品亚洲一区二区| 美女大奶头视频| 好男人在线观看高清免费视频 | 精品免费久久久久久久清纯| svipshipincom国产片| www.自偷自拍.com| 男男h啪啪无遮挡| 精品久久久久久,| 日韩av在线大香蕉| 国产成人av教育| 黑人巨大精品欧美一区二区mp4| 88av欧美| 少妇的丰满在线观看| 久久精品成人免费网站| 亚洲久久久国产精品| 成人精品一区二区免费| 亚洲av成人av| 亚洲精品一卡2卡三卡4卡5卡| 18禁国产床啪视频网站| 在线av久久热| 亚洲一区二区三区色噜噜| av在线天堂中文字幕| 高清在线国产一区| 美女大奶头视频| av福利片在线| 精品福利观看| 69av精品久久久久久| 亚洲成人精品中文字幕电影| √禁漫天堂资源中文www| 亚洲国产高清在线一区二区三 | 人人妻人人澡欧美一区二区 | 黄色视频,在线免费观看| 欧美日韩福利视频一区二区| 欧美不卡视频在线免费观看 | 97超级碰碰碰精品色视频在线观看| av免费在线观看网站| 久久人妻熟女aⅴ| 黄色毛片三级朝国网站| 91麻豆精品激情在线观看国产| 多毛熟女@视频| 久久香蕉国产精品| 午夜激情av网站| 日韩有码中文字幕| 一a级毛片在线观看| 亚洲精品av麻豆狂野| 免费人成视频x8x8入口观看| 亚洲成人久久性| 日本三级黄在线观看| 亚洲人成77777在线视频| 国产视频一区二区在线看| e午夜精品久久久久久久| 他把我摸到了高潮在线观看| 麻豆一二三区av精品| 这个男人来自地球电影免费观看| 亚洲 国产 在线| 精品久久久久久成人av| 久久久久久大精品| 国产亚洲精品久久久久久毛片| 日本vs欧美在线观看视频| 久久亚洲精品不卡| 精品久久久久久久人妻蜜臀av | 香蕉丝袜av| 精品少妇一区二区三区视频日本电影| 久久久精品国产亚洲av高清涩受| 他把我摸到了高潮在线观看| 香蕉久久夜色| 99在线视频只有这里精品首页| av欧美777| 99riav亚洲国产免费| 好男人电影高清在线观看| 国产片内射在线| 久久精品人人爽人人爽视色| 国产亚洲精品久久久久久毛片| 精品午夜福利视频在线观看一区| 欧美激情 高清一区二区三区| 午夜日韩欧美国产| 免费看十八禁软件| 巨乳人妻的诱惑在线观看| 色综合婷婷激情| 波多野结衣一区麻豆| 色综合亚洲欧美另类图片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产av一区二区精品久久| 免费在线观看视频国产中文字幕亚洲| 欧美日韩福利视频一区二区| 18禁国产床啪视频网站| 激情在线观看视频在线高清| 国产成人一区二区三区免费视频网站| 免费av毛片视频| 激情在线观看视频在线高清| 女人被狂操c到高潮| 精品久久久久久久毛片微露脸| 身体一侧抽搐| 精品第一国产精品| 最近最新中文字幕大全免费视频| 人人妻,人人澡人人爽秒播| 脱女人内裤的视频| 在线观看免费日韩欧美大片| 丁香欧美五月| 亚洲 国产 在线| 脱女人内裤的视频| 青草久久国产| 黄色女人牲交| 99国产综合亚洲精品| 久久精品国产亚洲av香蕉五月| 嫩草影视91久久| 亚洲专区字幕在线| xxx96com| 高清黄色对白视频在线免费看| 亚洲少妇的诱惑av| 在线观看舔阴道视频| 免费不卡黄色视频| or卡值多少钱| 欧美日本亚洲视频在线播放| 精品久久久精品久久久| 操美女的视频在线观看| 大陆偷拍与自拍| 午夜福利18| 97人妻天天添夜夜摸| 91老司机精品| 欧美久久黑人一区二区| 最新美女视频免费是黄的| 午夜久久久久精精品| 正在播放国产对白刺激| 女人被狂操c到高潮| av视频在线观看入口| 两个人免费观看高清视频| 亚洲专区国产一区二区| 精品第一国产精品| 国产熟女午夜一区二区三区| 女生性感内裤真人,穿戴方法视频| 成人av一区二区三区在线看| 亚洲 欧美一区二区三区| 国产欧美日韩综合在线一区二区| videosex国产| 亚洲少妇的诱惑av| 99国产精品一区二区三区| 欧美成人免费av一区二区三区| 欧美在线一区亚洲| 色综合亚洲欧美另类图片| 中文亚洲av片在线观看爽| 日韩大尺度精品在线看网址 | 精品高清国产在线一区| 97超级碰碰碰精品色视频在线观看| 国内毛片毛片毛片毛片毛片| 婷婷六月久久综合丁香| 亚洲精品国产一区二区精华液| 国产成人精品久久二区二区91| 欧美黑人精品巨大| 日韩中文字幕欧美一区二区| 黑人欧美特级aaaaaa片| 别揉我奶头~嗯~啊~动态视频| 亚洲av日韩精品久久久久久密| 乱人伦中国视频| 天天一区二区日本电影三级 | 亚洲人成网站在线播放欧美日韩| 高潮久久久久久久久久久不卡| 69精品国产乱码久久久| www.熟女人妻精品国产| 后天国语完整版免费观看| 黄色视频不卡| av中文乱码字幕在线| 欧美激情久久久久久爽电影 | 国产av在哪里看| 老熟妇乱子伦视频在线观看| 日韩欧美国产在线观看| 成年人黄色毛片网站| 国产成人精品在线电影| 免费少妇av软件| 精品不卡国产一区二区三区| 人人澡人人妻人| 性少妇av在线| cao死你这个sao货| 88av欧美| 涩涩av久久男人的天堂| 久热爱精品视频在线9| 91老司机精品| 夜夜夜夜夜久久久久| 欧美日韩中文字幕国产精品一区二区三区 | 一卡2卡三卡四卡精品乱码亚洲| 国产成人影院久久av| 亚洲中文av在线| 国产xxxxx性猛交| 久久精品91蜜桃| 好看av亚洲va欧美ⅴa在| 丰满人妻熟妇乱又伦精品不卡| 99香蕉大伊视频| 国产精品精品国产色婷婷| 亚洲精品一区av在线观看| 亚洲av五月六月丁香网| 成人欧美大片| 国产亚洲欧美精品永久| 精品欧美国产一区二区三| 国产伦人伦偷精品视频| 国产主播在线观看一区二区| 12—13女人毛片做爰片一| 欧美国产精品va在线观看不卡| 精品久久久久久久久久免费视频| 夜夜夜夜夜久久久久| 高清毛片免费观看视频网站| 亚洲av电影在线进入| 国产av一区二区精品久久| 一区二区三区高清视频在线| 日韩精品中文字幕看吧| 国内精品久久久久精免费| 亚洲欧美精品综合久久99| 一区在线观看完整版| 黄频高清免费视频| 精品国产一区二区久久| 我的亚洲天堂| 午夜福利高清视频| 久久婷婷成人综合色麻豆| 久久久久国产精品人妻aⅴ院| 国产精品电影一区二区三区| 欧美一级毛片孕妇| 超碰成人久久| 成人av一区二区三区在线看| 大型黄色视频在线免费观看| 男女下面进入的视频免费午夜 | 免费无遮挡裸体视频| 色播在线永久视频| 精品国产乱子伦一区二区三区| 亚洲国产中文字幕在线视频| 国产成人av教育| 如日韩欧美国产精品一区二区三区| 国产高清有码在线观看视频 | 亚洲专区字幕在线| 国产精品电影一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 成人国语在线视频| 中文字幕av电影在线播放| 91成人精品电影| 亚洲五月天丁香| 丁香六月欧美| 国产又爽黄色视频| a级毛片在线看网站| 久9热在线精品视频| 香蕉久久夜色| 婷婷六月久久综合丁香| 一进一出抽搐动态| 成人国产综合亚洲| 亚洲第一青青草原| 男女之事视频高清在线观看| 亚洲精品国产一区二区精华液| 亚洲精品在线美女| 亚洲久久久国产精品| 国产国语露脸激情在线看| 亚洲中文日韩欧美视频| 大型黄色视频在线免费观看| 搡老妇女老女人老熟妇| 两个人看的免费小视频| 伊人久久大香线蕉亚洲五| 精品免费久久久久久久清纯| 亚洲免费av在线视频| 亚洲国产日韩欧美精品在线观看 | 12—13女人毛片做爰片一| 久久精品人人爽人人爽视色| 女人被狂操c到高潮| 在线观看66精品国产| 久久精品影院6| 美女扒开内裤让男人捅视频| 精品国产超薄肉色丝袜足j| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩一区二区三区精品不卡| 无遮挡黄片免费观看| 日韩av在线大香蕉| 国产1区2区3区精品| 亚洲片人在线观看| 久久久精品国产亚洲av高清涩受| 九色亚洲精品在线播放| 精品久久久久久久毛片微露脸| 制服人妻中文乱码| 97人妻精品一区二区三区麻豆 | 1024视频免费在线观看| 黑人操中国人逼视频| 国产91精品成人一区二区三区| 久久久精品国产亚洲av高清涩受| 女性被躁到高潮视频| 国产欧美日韩精品亚洲av| 校园春色视频在线观看| 嫩草影视91久久| 99riav亚洲国产免费| 欧美乱码精品一区二区三区| 欧美激情 高清一区二区三区| 亚洲第一欧美日韩一区二区三区| 亚洲中文字幕日韩| 搡老妇女老女人老熟妇| 高清在线国产一区| 久久婷婷成人综合色麻豆| 大陆偷拍与自拍| 两个人看的免费小视频| 亚洲国产精品sss在线观看| 亚洲成av人片免费观看| 99久久99久久久精品蜜桃| 久久中文看片网| 久久久久国产一级毛片高清牌| 色哟哟哟哟哟哟| 狠狠狠狠99中文字幕| 成人欧美大片| 久久狼人影院| 午夜精品在线福利| 天天一区二区日本电影三级 | 国产人伦9x9x在线观看| 亚洲成人久久性| 一夜夜www| 怎么达到女性高潮| 一区福利在线观看| 国产精华一区二区三区| 亚洲国产欧美网| 国产欧美日韩一区二区三| 波多野结衣一区麻豆| 又黄又粗又硬又大视频| 精品国产美女av久久久久小说| 亚洲五月婷婷丁香| 精品久久久久久久久久免费视频| 午夜影院日韩av| 亚洲色图 男人天堂 中文字幕| 又黄又粗又硬又大视频| 精品国产美女av久久久久小说| 嫁个100分男人电影在线观看| 国内精品久久久久精免费| 亚洲国产精品999在线| 两性夫妻黄色片| 好看av亚洲va欧美ⅴa在| 一a级毛片在线观看| 欧美色视频一区免费| 亚洲午夜精品一区,二区,三区| 嫩草影院精品99| 欧美日韩亚洲国产一区二区在线观看| 在线观看66精品国产| 美女高潮到喷水免费观看| 欧美国产日韩亚洲一区| 精品人妻1区二区| 在线观看免费日韩欧美大片| 久久久久国产精品人妻aⅴ院| 国产私拍福利视频在线观看| 最新在线观看一区二区三区| 一区二区三区激情视频| 久久香蕉国产精品| 成年女人毛片免费观看观看9| 亚洲国产高清在线一区二区三 | av免费在线观看网站| 90打野战视频偷拍视频| 最近最新中文字幕大全电影3 | 一本大道久久a久久精品| 亚洲男人天堂网一区| 亚洲欧美精品综合久久99| 好男人电影高清在线观看| 日韩大码丰满熟妇| 黑人操中国人逼视频| 88av欧美| 91麻豆av在线| 激情视频va一区二区三区| 午夜亚洲福利在线播放| 亚洲欧美日韩高清在线视频| 国产av在哪里看| 国产成人影院久久av| 国产熟女xx| 亚洲第一青青草原| 18禁黄网站禁片午夜丰满| 一级黄色大片毛片| 国产一区二区三区在线臀色熟女| 91大片在线观看| 午夜福利18| 操出白浆在线播放| 九色亚洲精品在线播放| 久久久精品欧美日韩精品| xxx96com| 午夜福利欧美成人| 欧美日韩福利视频一区二区| 国产成人精品久久二区二区免费| 村上凉子中文字幕在线| 电影成人av| 50天的宝宝边吃奶边哭怎么回事| 久久国产精品男人的天堂亚洲| 很黄的视频免费| 久久久久亚洲av毛片大全| 极品人妻少妇av视频| 在线视频色国产色| 在线av久久热| 久热这里只有精品99| 国产日韩一区二区三区精品不卡| 亚洲精品美女久久av网站| 亚洲中文日韩欧美视频| 成人免费观看视频高清| netflix在线观看网站| 国产精品av久久久久免费| 色婷婷久久久亚洲欧美| 99riav亚洲国产免费| 身体一侧抽搐| av片东京热男人的天堂| 亚洲片人在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产精品爽爽va在线观看网站 | 又黄又粗又硬又大视频| 母亲3免费完整高清在线观看| 中文亚洲av片在线观看爽| 久9热在线精品视频| 久久欧美精品欧美久久欧美| 中出人妻视频一区二区| 日韩欧美一区二区三区在线观看| 91老司机精品| 日韩视频一区二区在线观看| 手机成人av网站| 嫁个100分男人电影在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成伊人成综合网2020| 国产不卡一卡二| 国产aⅴ精品一区二区三区波| 亚洲熟妇熟女久久| 一本大道久久a久久精品| 久久久久久久久久久久大奶| 色哟哟哟哟哟哟| 可以在线观看的亚洲视频| 久久人人精品亚洲av| 欧洲精品卡2卡3卡4卡5卡区| 久久精品影院6| 咕卡用的链子| 久久精品91无色码中文字幕| 不卡一级毛片| 国产97色在线日韩免费| 国产蜜桃级精品一区二区三区| 精品第一国产精品| 国产一区二区三区在线臀色熟女| 一区二区三区国产精品乱码| 91精品三级在线观看| 国产亚洲精品综合一区在线观看 | 美女高潮到喷水免费观看| av电影中文网址| a级毛片在线看网站| 国产精品1区2区在线观看.| 91国产中文字幕| 在线观看66精品国产| 久久久国产成人免费| 最好的美女福利视频网| 高清在线国产一区| av有码第一页| 美女国产高潮福利片在线看| 国产欧美日韩一区二区精品| 天天躁狠狠躁夜夜躁狠狠躁| 最新美女视频免费是黄的| 久久久国产精品麻豆| 女性被躁到高潮视频| 91国产中文字幕| 亚洲伊人色综图| 免费久久久久久久精品成人欧美视频| 十八禁网站免费在线| 好看av亚洲va欧美ⅴa在| 亚洲精华国产精华精| 国产高清视频在线播放一区| 91麻豆av在线| 天堂动漫精品| 老司机深夜福利视频在线观看| 日韩大码丰满熟妇| 人人妻,人人澡人人爽秒播| 国产精品久久电影中文字幕| 国产黄a三级三级三级人| 欧美绝顶高潮抽搐喷水| 两性夫妻黄色片| 国产三级在线视频| 久久午夜综合久久蜜桃| 国产精品自产拍在线观看55亚洲| 少妇 在线观看| 国产麻豆成人av免费视频| 一a级毛片在线观看| 人妻丰满熟妇av一区二区三区| 亚洲人成77777在线视频| 国产精品99久久99久久久不卡| 午夜精品久久久久久毛片777| 亚洲国产精品999在线| 九色国产91popny在线| 国产精品自产拍在线观看55亚洲| 在线观看66精品国产| 无人区码免费观看不卡| 男女做爰动态图高潮gif福利片 | 国产91精品成人一区二区三区| 正在播放国产对白刺激| 欧美黄色片欧美黄色片| 99国产精品99久久久久| 亚洲aⅴ乱码一区二区在线播放 | 久久国产精品人妻蜜桃| or卡值多少钱| 亚洲自偷自拍图片 自拍| 亚洲精品中文字幕一二三四区| 成人18禁高潮啪啪吃奶动态图| 热re99久久国产66热| АⅤ资源中文在线天堂| 人妻丰满熟妇av一区二区三区| 亚洲精品国产精品久久久不卡| 中文字幕最新亚洲高清| 美女免费视频网站| 性欧美人与动物交配| 在线视频色国产色| 中文字幕色久视频| 丁香欧美五月| 亚洲欧美精品综合久久99| 国产视频一区二区在线看| 成人亚洲精品一区在线观看| 黄色a级毛片大全视频| 亚洲av美国av| 正在播放国产对白刺激| 一进一出好大好爽视频| 精品久久蜜臀av无| 99re在线观看精品视频| 久久久国产成人精品二区| 日韩精品青青久久久久久| 亚洲人成电影免费在线| 999久久久国产精品视频| 久9热在线精品视频| 黑人巨大精品欧美一区二区蜜桃| 午夜老司机福利片| 国产野战对白在线观看| 精品久久蜜臀av无| 欧美黑人精品巨大| 九色国产91popny在线| 极品教师在线免费播放| 在线观看www视频免费| 老司机午夜福利在线观看视频| 老汉色av国产亚洲站长工具| 少妇熟女aⅴ在线视频| 亚洲人成电影免费在线| av视频免费观看在线观看| 国产亚洲精品av在线| 在线观看免费日韩欧美大片| 日本欧美视频一区| 啦啦啦韩国在线观看视频| 久久亚洲真实| 精品人妻1区二区| 国产精品一区二区免费欧美| 狠狠狠狠99中文字幕| 日韩免费av在线播放| 高清在线国产一区| 91av网站免费观看| 少妇粗大呻吟视频| www.熟女人妻精品国产| 久久人妻福利社区极品人妻图片| 欧美亚洲日本最大视频资源| 国产成人啪精品午夜网站| 亚洲狠狠婷婷综合久久图片| 成人特级黄色片久久久久久久| 精品久久蜜臀av无| 琪琪午夜伦伦电影理论片6080| 在线观看66精品国产| 欧美黄色淫秽网站| 啦啦啦 在线观看视频| 黄色成人免费大全| 亚洲成国产人片在线观看| 啦啦啦 在线观看视频| 国产精品免费视频内射| 欧美黄色淫秽网站| 国产黄a三级三级三级人| 亚洲专区国产一区二区| 欧美中文日本在线观看视频| 国产一区二区在线av高清观看| 一级黄色大片毛片| 久久精品国产99精品国产亚洲性色 | 亚洲成av片中文字幕在线观看| 国产精品国产高清国产av| 国产精品美女特级片免费视频播放器 | 久久久久久亚洲精品国产蜜桃av| 亚洲精品美女久久久久99蜜臀| 久久久久国产精品人妻aⅴ院| 精品国产一区二区久久| 美女高潮喷水抽搐中文字幕| 老熟妇乱子伦视频在线观看| 午夜老司机福利片| 99精品欧美一区二区三区四区| 9热在线视频观看99| 桃色一区二区三区在线观看| 亚洲精品美女久久av网站| 变态另类成人亚洲欧美熟女 | 在线免费观看的www视频| 多毛熟女@视频| 国产精品av久久久久免费| 亚洲国产欧美日韩在线播放| 高清在线国产一区| 亚洲av第一区精品v没综合| 自线自在国产av| 久久国产精品人妻蜜桃| 日韩av在线大香蕉| 久久人人精品亚洲av|