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

    植被指數(shù)對藏北高寒草甸干旱的敏感性分析

    2018-03-05 08:51:30周艷蓮居為民戴聲佩包永康
    江蘇農(nóng)業(yè)科學 2018年1期
    關鍵詞:時間尺度草甸植被指數(shù)

    李 晗,周艷蓮,居為民,戴聲佩,宋 練,包永康

    (1.南京大學國際地球系統(tǒng)科學研究所,江蘇南京 210046;2.南京大學地理與海洋科學學院,江蘇南京 210046;3.江蘇省地理信息技術(shù)重點實驗室,江蘇南京 210046)

    當前人類活動的不斷加劇會對全球大氣循環(huán)的格局造成越來越嚴重的影響。許多大氣循環(huán)模型預測未來全球氣候變化的特征,結(jié)果是極端降水事件和極端干旱事件發(fā)生的頻率會顯著增加[1-4]。植被生長季干旱嚴重會影響植被指數(shù)、地表溫度、土壤濕度、反射率、蒸散等地表生物物理參數(shù)[5]。

    干旱指數(shù)是反映氣候干旱程度的指標,許多學者利用氣象資料、水文資料、土壤資料或衛(wèi)星遙感資料等來建立不同的干旱指數(shù),從而進行干旱監(jiān)測預警。資料顯示常用的干旱指數(shù)有50多種,例如標準化降水指數(shù)(SPI)[6]、Palmer干早指數(shù)(PDSI)[7]、Z指數(shù)[8]、綜合氣象干旱指數(shù)(CI)[9]等。這些指數(shù)在一定程度上為全球或者地區(qū)的干旱監(jiān)測提供了幫助,但都不可避免地需要降水數(shù)據(jù)的地面采集[10-11]。

    與常規(guī)地面干旱監(jiān)測相比,利用遙感手段進行干旱監(jiān)測具有時空分辨率相對較高、傳感器覆蓋面廣、數(shù)據(jù)獲取較為容易、消耗較少人力物力等特點。隨著遙感對地觀測能力的提升,大量研究者利用遙感數(shù)據(jù)不同波段組合獲取植被指數(shù),通過簡化復雜的干旱現(xiàn)象評估干旱的嚴重程度。遙感數(shù)據(jù)可以獲取與植物生長有關的植被指數(shù),如歸一化差值植被指數(shù)(NDVI)和增強型植被指數(shù)(EVI)等。NDVI可以反映年尺度和季節(jié)尺度上植被生長活動的變化,由于它進行了比值和歸一化處理,在一定程度上可以消除如太陽高度角、衛(wèi)星觀測角、地形、云、陰影和大氣條件有關的輻照度變化等的部分影響[12]。EVI利用藍光波段和紅光波段通過大氣氣溶膠時的差別,補償了殘留氣溶膠對紅光的吸收,綜合采用抗大氣植被指數(shù)和抗土壤植被指數(shù),克服了土壤背景的影響和NDVI在植被高覆蓋區(qū)易飽和、植被低覆蓋區(qū)受土壤植被影響較大、對大氣衰減去除不徹底等缺點[13]。Song等利用中國地區(qū)1982—2005年GIMMS NDVI數(shù)據(jù)集和降水數(shù)據(jù)分析認為,在所有氣象要素中,每月NDVI的最大值與降水的相關性最強,尤其是在干旱與半干旱地區(qū)其相關性表現(xiàn)最為明顯[14]。趙偉等采用NDVI時間序列諧波分析法(HANTS)對中分辨率成像光譜儀MODIS產(chǎn)品的EVI數(shù)據(jù)進行去云處理,重構(gòu)EVI數(shù)據(jù)并結(jié)合地面的氣象數(shù)據(jù),分析表明處理后的EVI時間序列數(shù)據(jù)能夠較好地反映干旱對地表植被的影響[12]。但是,NDVI、EVI等均是比較保守的指數(shù),在葉片受到外界環(huán)境脅迫而導致水分含量發(fā)生變化時,NDVI與EVI自身不一定能及時地響應[15-16]。因此,近年來一些與水分相關的植被指數(shù)被越來越多地用于植被水分含量的估計與干旱監(jiān)測。無論是物理輻射傳輸模型還是試驗統(tǒng)計均表明,植物的水分含量在很大程度上影響著葉片在短波紅外波段下的反射率,由于這個區(qū)域是葉片水分光譜的吸收范圍,所以植被在短波紅外的反射率與它們?nèi)~片的水分含量是呈負相關的[16]。研究表明,將短波紅外波段反射率和近紅外波段反射率一起組合建立如歸一化差值水體指數(shù)(NDWI)和地表水分指數(shù)(LSWI)等植被指數(shù)可以對植被水分含量進行監(jiān)測[17]。一些國外學者研究發(fā)現(xiàn),與植被水分有關的植被指數(shù)(如LSWI)在植被干旱監(jiān)測上的敏感性優(yōu)于與植被生長狀況有關的植被指數(shù)[18-19]。Bajgain等以美國季風區(qū)高桿草原為研究區(qū),分析認為當干旱發(fā)生時,LSWI在年際尺度上的下降幅度優(yōu)于NDVI與EVI,當土壤水分含量發(fā)生變化時,LSWI更能對變化做出及時的響應[20]。但目前在國內(nèi),尤其是關于半干旱區(qū)高原草甸不同類型植被指數(shù)對干旱的敏感性研究尚少,不同時間尺度、不同類型植被指數(shù)對干旱追蹤能力的分析也較欠缺,因此利用不同類型植被指數(shù)對半干旱區(qū)高原草甸干旱的敏感性進行系統(tǒng)的分析,就顯得很有必要。

    高寒草甸是青藏高原廣泛分布的典型植被類型之一,面積約為1.2×106km2,是亞洲中部高寒地區(qū)乃至世界高寒地區(qū)極具代表性的植被[21-22]。水分作為高寒草甸這種干旱、半干旱區(qū)草原植物生長發(fā)育的限制性資源,其缺失會直接導致區(qū)域內(nèi)地生物量的減少,而草原生態(tài)系統(tǒng)是陸地生態(tài)系統(tǒng)中對干旱非常敏感的系統(tǒng)[1]。因此,對高寒草甸的干旱進行分析與評估是至關重要的。目前由于人類活動的加劇,藏北牧區(qū)超載放牧現(xiàn)象十分嚴重,極端天氣事件時有發(fā)生,草地利用壓力很大,因此對類似于當雄草原這種典型高寒草甸植被的干旱監(jiān)測研究顯得很有必要[21]。本研究以高寒草甸干旱的監(jiān)測與評估為目的,以MODIS的反射率數(shù)據(jù)為遙感數(shù)據(jù)源,計算NDVI、EVI、LSWI等植被指數(shù),結(jié)合同時相內(nèi)觀測的降水與土壤濕度數(shù)據(jù),分析在不同時間尺度上,不同類型植被指數(shù)對干旱的敏感性,篩選出適宜高寒草甸干旱監(jiān)測與評估的最佳遙感指數(shù)。

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

    1.1 研究地概況

    西藏自治區(qū)拉薩市當雄縣草原站,地理坐標為 30°29′50.21″N,91°03′58.90″E,地處念青唐古拉山的南緣。當雄草原植被屬于典型的高寒嵩草草甸植被,建群種有高山嵩草、矮嵩草、絲穎針茅,伴生種有多種密叢生嵩草,本類草甸由于超載過牧而退化,針茅類入侵較多,常形成草原草甸,群落覆蓋度約80%。氣候為高原季風氣候,年平均氣溫1.3 ℃,最冷月1月最低平均氣溫-10.4 ℃,最熱月7月10.7 ℃。多年平均降水量450 mm,其中85%集中在6—8月。冰凍期較長,持續(xù)時間為11月至次年1月。年太陽總輻射 7 527.6 MJ/m2,光合有效輻射3 213.3 MJ/m2。土壤類型為高寒草甸土,土壤厚度為30~50 cm。

    1.2 降水與土壤濕度數(shù)據(jù)

    采用2個不同深度(5、10 cm)的土壤水分含量(SWC)來表征土壤的濕潤程度。同時,為了與MODIS地表反射率數(shù)據(jù)匹配,對站點獲取的日平均降水和土壤水分含量數(shù)據(jù)分別作求和與平均處理,獲得時間尺度為8 d的數(shù)據(jù),進一步計算得到月平均降水、年平均降水和土壤水分含量。

    1.3 遙感數(shù)據(jù)源及數(shù)據(jù)處理

    MODIS數(shù)據(jù)可以同時獲得可見光至近紅外的36個波段,常被用于植被遙感研究。本研究通過美國航空航天宇航局的戈達德航天中心數(shù)據(jù)網(wǎng)站(LAADS,https://ladsweb.nascom.nasa.gov/)下載覆蓋當雄通量觀測站點每8 d的合成地表反射率產(chǎn)品MOD09A1,時間序列為2004年至2012年,空間分辨率為500 m,該產(chǎn)品主要包括MODIS傳感器1~7波段的地表反射率數(shù)據(jù)和各個波段的質(zhì)量控制數(shù)據(jù),同時還包括反射率數(shù)據(jù)狀態(tài)質(zhì)量保證數(shù)據(jù)。利用MODIS數(shù)據(jù)的處理工具(MRT)將投影方式由Sinusoidal投影轉(zhuǎn)換為UTM/WGS 84投影,并通過最近鄰法重采樣,獲取站點中心單像元紅光波段(620~670 nm)、近紅外波段1(841~876 nm)、藍光波段(459~479 nm)以及短波紅外波段1(1 628~1 652 nm)反射率數(shù)據(jù)。根據(jù)公式(1~3)分別計算不同時間尺度的NDVI、EVI和LSWI。

    (1)

    (2)

    (3)

    式中:ρX為MOIDS數(shù)據(jù)X波段的反射率;公式(2)中C1和C2是氣溶膠等大氣影響物的校正系數(shù),L為背景(土壤)調(diào)整系數(shù),G為常量,取值分別為C1=6,C2=7.5,L=1,G=1。

    由于LSWI與天氣條件有關,無法像葉面積指數(shù)(LAI)和光合作用有效輻射比率(FPAR)一樣進行數(shù)據(jù)平滑處理[21]。本研究結(jié)合質(zhì)量控制數(shù)據(jù),采用兩步差補方法處理時間序列中受云影響、質(zhì)量較差的植被指數(shù)數(shù)據(jù)[22]。NDVI、EVI和LSWI均未作平滑處理。

    2 結(jié)果與分析

    2.1 年時間尺度敏感性分析

    圖1表明,2004—2012年,拉薩當雄高原草甸站點植被生長季的年平均降水量為428 mm(±126 mm)。年平均降水量大于2004—2012年平均降水量與標準差之和(554 mm)的年份被認為是濕潤年,而年平均降水量小于2004—2012年平均降水量與標準差之差(302 mm)的年份被認為是干旱年。拉薩當雄高原草甸站點2006、2009年的降水量分別只有241、291 mm,因此為干旱年。而2008年的降水量為648 mm,因此為濕潤年。干旱年份(2006、2009年)植被生長季(5—9月)的不同深度(5、10 cm)土壤水分含量均<9%,而濕潤年份(2008年)土壤水分含量為12.6%。

    圖2表明,干旱年份(2006、2009年)植被生長季平均NDVI(NDVIavg)、EVI(EVIavg)和LSWI(LSWIavg)在分析的時間序列中均較低(NDVIavg<0.28、EVIavg<0.18、LSWIavg<0),而在濕潤年份(2008年)植被生長季NDVI、EVI和LSWI均達到最大值(NDVIavg=0.41、EVIavg=0.28、LSWIavg=0.07)。

    干旱年份(2006、2009年)植被生長季植被指數(shù)最大值NDVI(NDVImax)、EVI(EVImax)和LSWI(LSWImax)均取得最小值(NDVImax<0.44,EVImax<0.29,LSWImax<0.07),在濕潤年份(2008年)植被生長季植被指數(shù)只有最大值NDVI和LSWI達到峰值(NDVImax=0.64,LSWImax=0.20),而EVImax僅為0.48,遠遠小于2011年(EVImax=0.67)。

    植被指數(shù)均值和最大值的變化幅度可以用來表征植被指數(shù)對干旱的相對敏感性。變化幅度用當年植被指數(shù)的均值或最大值減去多年相應的平均值表示。變化幅度百分比用當年植被指數(shù)的均值或最大值減去多年相應的平均值后的結(jié)果再除以多年相應的平均值表示。如圖3所示,2006年LSWIavg與2004—2012年平均值相比,減少了0.052(316%),減少幅度最大。2009年的結(jié)果與2006年相似,LSWIavg、EVIavg、NDVIavg

    分別比2004—2012年平均值減少0.037(227%)、0.068(16%)、0.054(29%)。植被生長季LSWImax與NDVImax、EVImax相比也表現(xiàn)出同樣的敏感性,2006年LSWImax減少0.13(96%),均大于EVImax減少的0.14(35%)與NDVImax減少的0.11(22%)。同樣的,2009年LSWImax減少0.07(52%),均大于EVImax減少的0.13(32%)與NDVImax減少的0.08(16%)。

    濕潤年份(2008年)LSWI持續(xù)大于0的時間最大,而2006年僅有1個LSWI大于0。由于氣溫低,海拔高,4月中旬生長季才剛剛開始。經(jīng)過春季生長之后,植被生長的速率很大一部分取決于土壤水分含量是否充足。若土壤水分含量不充足,植被的生長就會受到抑制。相反,干旱年份LSWI小于0的天數(shù)多于濕潤年份與正常年份,而濕潤年份LSWI大于0.1的時間多于干旱年份與正常年份。如圖4所示,LSWI小于0的天數(shù)隨降水量的增加整體是下降的。2006年LSWI小于0的天數(shù)為120 d,大于2007年的72 d和2008年的 40 d,而2006年LSWI大于0.1的天數(shù)僅為8 d,而2008年為88 d,大于2007年的56 d。

    2.2 月時間尺度敏感性分析

    在分析以植被指數(shù)為基礎的干旱指標對降水的敏感性時,易忽略植被指數(shù)和降水的同時相季節(jié)變化,可能對其相關性產(chǎn)生影響。本研究為消除季節(jié)變化引起的降水與植被指數(shù)的正弦波動,對月時間尺度的降水數(shù)據(jù)和植被指數(shù)數(shù)據(jù)均作了距平處理。具體做法是用當月降水的總累積量和植被指數(shù)的平均值減去多年的月平均值。圖5顯示,當前月份LSWI與降水距平的相關性最高,r2達到0.32(P<0.001),其次是NDVI,r2為0.16(P<0.1),EVI的r2最低,僅為0.05。而在1個月滯后的情況下,NDVI與EVI的r2分別為0.20(P<0.01)和0.24(P<0.01),與當前月份降水與植被指數(shù)距平相關性相比有一定的提高,而LSWI和降水距平的r2僅為0.17(P<0.1),有一定的下降。

    3 結(jié)論與討論

    相比利用氣象、水文、土壤等資料建立干旱指數(shù)用于干旱監(jiān)測等方法而言,以植被指數(shù)為核心的遙感干旱監(jiān)測指數(shù)具有很多優(yōu)勢。其成本低、時空分辨率較高、對降水能夠及時響應,而且它對地表植被的葉綠素含量、水分含量均比較敏感,優(yōu)勢明顯。植被在不同波段的光譜特征能反映指數(shù)的不同要素或者某種特征的差異,植被指數(shù)利用多個光譜波段的反射率經(jīng)分析運算(線性或非線性組合方式),產(chǎn)生對植被長勢、水分含量等有一定指示意義的數(shù)值。這是利用植被指數(shù)進行干旱監(jiān)測與評估的理論基礎[20]。Chandrasekar等以印度南岸的城市——安得拉邦為研究區(qū),以整個地區(qū)MODIS反射率數(shù)據(jù)為基礎,結(jié)合氣象數(shù)據(jù)分析認為,LSWI在低降水地區(qū)與降水的相關性更好,優(yōu)于NDVI;而在高降水地區(qū),LSWI與降水的相關性要差于NDVI,且都存在滯后現(xiàn)象[23]。但是該研究沒有考慮因季節(jié)變化引起的降水和植被指數(shù)在年內(nèi)的正弦波動,這種季節(jié)特性在一定程度上影響了結(jié)果的可靠性。本研究利用NDVI、EVI、LSWI這3種指數(shù)對高寒草甸進行干旱監(jiān)測,同時為消除季節(jié)變化引起的降水與植被指數(shù)的正弦波動,對降水和植被指數(shù)在月尺度上進行距平處理,結(jié)果顯示月時間尺度上LSWI對降水的敏感性優(yōu)于NDVI與EVI,且LSWI不存在滯后現(xiàn)象。這與草原生態(tài)系統(tǒng)下墊面比較均一、受背景反射率影響較小等原因有關。2006年夏季站點周圍發(fā)生了嚴重的干旱現(xiàn)象,LSWI在8、9月已經(jīng)顯示出明顯的下降趨勢,而NDVI和EVI直到10月才開始下降。進一步說明LSWI對降水的敏感性在月尺度上較NDVI、EVI更好,且不存在滯后現(xiàn)象(圖6)。

    Bajgain以美國季風區(qū)的2個高桿草原站點為研究區(qū)、以遙感獲取的LSWI為基礎對干旱進行的分級(LSWI>0.1;00.1(0)表示正常。海北高寒草甸位于半干旱區(qū),正常年份均值處于0

    以植被指數(shù)為核心的遙感干旱監(jiān)測指數(shù)雖然具有很多優(yōu)勢,但在其大范圍應用之前應該對不同的植被類型進行遙感干旱檢測研究,以確定合適的研究指標和研究手段。本研究的目的就是探討藏北典型高寒草甸干旱檢測的最佳遙感指數(shù),以期更精確地對高寒草甸干旱事件的發(fā)生進行有效地監(jiān)測與評估。當然利用植被指數(shù)進行干旱監(jiān)測與評估也有諸多限制,例如遙感觀測資料的時間尺度,在天時間尺度或者更高時間尺度上的實時監(jiān)測目前來講仍存在極大挑戰(zhàn),不過隨著區(qū)域性小衛(wèi)星的陸續(xù)發(fā)射升空,也為干旱的實時遙感監(jiān)測提供了可能。

    對干旱指數(shù)的評價應該從多方面進行,綜合考慮。本研究評價干旱指數(shù)的效果時,采用被廣泛使用的相關系數(shù)r2和顯著性水平P值作為主要評價指標,也定性地分析了在干旱發(fā)生當月,最早作出反應的植被指數(shù),并驗證了以LSWI為基礎的干旱指標在半干旱區(qū)的適應性問題,進一步說明LSWI對草原干旱事件的響應優(yōu)于NDVI與EVI。

    本研究利用MODIS反射率數(shù)據(jù)計算NDVI、EVI和LSWI,結(jié)合降水與土壤濕度數(shù)據(jù)對藏北高寒草甸干旱進行監(jiān)測,干旱事件發(fā)生時植被生長受到了很大的限制,3種植被指數(shù)在干旱年份(2006年)減小幅度也基本相同,LSWI略為敏感。月尺度上LSWI與降水距平的相關性最高,r2達到0.32,其次是NDVI與EVI,且LSWI不存在滯后現(xiàn)象。LSWI<0的天數(shù)隨降水量的增加呈明顯的減少趨勢。利用LSWI對干旱進行分級(LSWI>0.1;0

    [1]張 彬,朱建軍,劉華民,等.極端降水和極端干旱事件對草原生態(tài)系統(tǒng)的影響[J].植物生態(tài)學報,2014,38(9):1008-1018.

    [2]International Panel on Climate Change.Climate change 2007:impacts,adaption and vulnerability[R].New York:Cambridge University Press,2007.

    [3]Easterling D R,Meehl G A,Parmesan C,et al.Climate extremes:observations,modeling,and impacts[J].Science,2000,289(5487):2068-2074.

    [4]Meehl G A,Tebaldi C.More intense,more frequent,and longer lasting heat waves in the 21st century[J].Science,2004,305(5686):994-997.

    [5]Ghulam A,Li Z L,Qin Q,et al.Exploration of the spectral space basedon vegetation index and lbedo for surface drought estimation[J].Journal of Applied Remote Sensing,2007,1(1):13529-13529.

    [6]王勁松,李耀輝,王潤元,等.我國氣象干旱研究進展評述[J].干旱氣象,2012,30(4):497-508.

    [7]楊小利,楊興國,馬鵬里,等.PDSI在甘肅中東部地區(qū)的修正和應用[J].地球科學進展,2005,20(9):1022-1028.

    [8]王勁松,黃玉霞,馮建英,等.徑流量Z指數(shù)與Palmer指數(shù)對河西干旱的監(jiān)測[J].應用氣象學報,2009,20(4):471-477.

    [9]李樹巖,劉榮花,師麗魁,等.基于CI指數(shù)的河南省近40年干旱特征分析[J].干旱氣象,2009,27(2):97-102.

    [10]王春林,郭 晶,薛麗芳,等.改進的綜合氣象干旱指數(shù)CInew及其適用性分析[J].中國農(nóng)業(yè)氣象,2011,32(4):621-626.

    [11]張紅衛(wèi),陳懷亮,申雙和,等.基于表層水分含量指數(shù)(SWCI)的土壤干旱遙感監(jiān)測[J].遙感技術(shù)與應用,2008,23(6):624-628.

    [12]趙 偉,李召良.利用MODIS/EVI時間序列數(shù)據(jù)分析干旱對植被的影響[J].地理科學進展,2007,26(6):40-46.

    [13]彭代亮,黃敬峰,王秀珍.基于MODIS-EVI區(qū)域植被季節(jié)變化與氣象因子的關系[J].應用生態(tài)學報,2007,18(5):983-989.

    [14]Song Y,Ma M G.A statistical analysis of the relationship between climatic factors and the normalized difference vegetation index in China[J].International Journal of Remote Sensing,2011,32(14):3947-3965.

    [15]劉小磊,覃志豪.NDWI與NDVI指數(shù)在區(qū)域干旱監(jiān)測中的比較分析——以2003年江西夏季干旱為例[J].遙感技術(shù)與應用,2007,22(5):608-612.

    [16]劉立文,張吳平,段永紅,等.TVDI模型的農(nóng)業(yè)旱情時空變化遙感應用[J].生態(tài)學報,2014,34(13):3704-3711.

    [17]Chen D Y,Huang J F,Jackson T J.Vegetation water content estimation for corn and soybeans using spectral indices derived from MODIS near- and short-wave infrared bands[J].Remote Sensing of Enviroment,2005,98(2):225-236.

    [18]張友水,謝元禮.MODIS影像的NDVI和LSWI植被水分含量估算[J].地理科學,2008,28(1):72-76.

    [19]劉元亮,李 艷,吳劍亮.基于LSWI和NDVI時間序列的水田信息提取研究[J].地理與地理信息科學,2015,31(3):32-37.

    [20]Bajgain R,Xiao X M,Wagle P,et al.Sensitivity analysis of vegetation indices to drought over two tallgrass prairie sites[J].ISPRS Journal of Photogrammetry and Remote Sensing,2015(108):151-160.

    [21]周宇庭,付 剛,沈振西,等.藏北典型高寒草甸地上生物量的遙感估算模型[J].草業(yè)學報,2013,22(1):120-129.

    [22]Zhang Y,Song C H,Sun G,et al.Understanding moisture stress on light use efficiency across terrestrial ecosystems based on global flux and remote sensing data[J].Journal of Geophysical Research:Biogeosciences,2015,120(10):2053-2066.

    [23]Chandrasekar K,Sesha S,Roy P S,et al.Land surface water index (LSWI) response to rainfall and NDVI using the MODIS vegetation index product[J].International Journal of Remote Sensing,2010,31(15):3987-4005.

    猜你喜歡
    時間尺度草甸植被指數(shù)
    高山草甸
    時間尺度上非完整系統(tǒng)的Noether準對稱性與守恒量
    The world's narrowest river
    山地草甸
    時間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    力學學報(2021年10期)2021-12-02 02:32:04
    交直流混合微電網(wǎng)多時間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    大連市暴雨多時間尺度研究分析
    久久精品夜夜夜夜夜久久蜜豆| 色噜噜av男人的天堂激情| 国内揄拍国产精品人妻在线| 男女之事视频高清在线观看| 日韩精品中文字幕看吧| 91九色精品人成在线观看| 久久精品91蜜桃| 午夜福利18| 长腿黑丝高跟| 一本久久中文字幕| 黄色丝袜av网址大全| 伊人久久精品亚洲午夜| 亚洲精品一区av在线观看| 九色成人免费人妻av| 亚洲精品影视一区二区三区av| 国产视频内射| 亚洲,欧美,日韩| 97碰自拍视频| 深爱激情五月婷婷| 精品无人区乱码1区二区| 免费人成在线观看视频色| 夜夜看夜夜爽夜夜摸| 国产伦精品一区二区三区视频9| 最近中文字幕高清免费大全6 | 国产中年淑女户外野战色| 国产探花在线观看一区二区| 久久精品人妻少妇| 亚洲性夜色夜夜综合| 国产精品美女特级片免费视频播放器| 中文字幕av在线有码专区| 两个人的视频大全免费| 久久中文看片网| 国产精品久久久久久人妻精品电影| 亚洲成a人片在线一区二区| 18+在线观看网站| 久久久久久九九精品二区国产| 男女床上黄色一级片免费看| 一级作爱视频免费观看| 国内精品久久久久精免费| 全区人妻精品视频| 在线免费观看的www视频| 一边摸一边抽搐一进一小说| 精品久久久久久久久久免费视频| 老熟妇仑乱视频hdxx| 亚洲av美国av| 最后的刺客免费高清国语| 欧美国产日韩亚洲一区| 亚洲最大成人中文| avwww免费| 女人十人毛片免费观看3o分钟| 老司机午夜福利在线观看视频| 亚洲av中文字字幕乱码综合| 无遮挡黄片免费观看| 欧美性感艳星| 中文字幕免费在线视频6| 亚洲最大成人中文| 一进一出好大好爽视频| 午夜精品在线福利| 波多野结衣巨乳人妻| 日韩人妻高清精品专区| 国产精品爽爽va在线观看网站| 日韩 亚洲 欧美在线| 天堂网av新在线| 欧美三级亚洲精品| 99热只有精品国产| 国产成人aa在线观看| 国产在线男女| 91av网一区二区| 成人午夜高清在线视频| 蜜桃久久精品国产亚洲av| 国产精品一区二区三区四区免费观看 | 精品一区二区三区视频在线| 国产精品人妻久久久久久| 嫩草影视91久久| 九九热线精品视视频播放| 亚洲人与动物交配视频| 亚洲,欧美,日韩| av国产免费在线观看| 男女之事视频高清在线观看| 精品久久久久久久久久久久久| 熟女电影av网| 国产精品日韩av在线免费观看| 熟女人妻精品中文字幕| 特级一级黄色大片| 精品久久久久久久久亚洲 | 中文字幕免费在线视频6| 国产乱人伦免费视频| 婷婷六月久久综合丁香| 久久午夜亚洲精品久久| 人妻丰满熟妇av一区二区三区| 国产一区二区三区视频了| 久久人人爽人人爽人人片va | 熟女电影av网| 亚洲熟妇中文字幕五十中出| 亚洲人成电影免费在线| 搡老熟女国产l中国老女人| 久久国产乱子免费精品| 成人午夜高清在线视频| 男女那种视频在线观看| 国产精品久久久久久亚洲av鲁大| av视频在线观看入口| 亚洲人成网站在线播| 国产白丝娇喘喷水9色精品| 日韩欧美一区二区三区在线观看| 午夜免费激情av| 精品久久久久久久久久免费视频| 国产精品综合久久久久久久免费| 日韩欧美免费精品| 别揉我奶头 嗯啊视频| 精品不卡国产一区二区三区| 亚洲av日韩精品久久久久久密| 两人在一起打扑克的视频| 搡老熟女国产l中国老女人| 美女免费视频网站| 日本 欧美在线| 真实男女啪啪啪动态图| ponron亚洲| or卡值多少钱| 国产精品亚洲av一区麻豆| 免费无遮挡裸体视频| 免费av毛片视频| 午夜精品在线福利| 欧美+日韩+精品| 久久亚洲精品不卡| 国产人妻一区二区三区在| 韩国av一区二区三区四区| 少妇人妻一区二区三区视频| 婷婷精品国产亚洲av在线| 国产一区二区激情短视频| 很黄的视频免费| 亚洲最大成人中文| 久久人人爽人人爽人人片va | 免费在线观看影片大全网站| 国产黄色小视频在线观看| 色av中文字幕| 国产精品综合久久久久久久免费| 国产黄a三级三级三级人| 在线看三级毛片| 亚洲国产精品合色在线| 88av欧美| 麻豆av噜噜一区二区三区| 欧美精品啪啪一区二区三区| 亚洲18禁久久av| bbb黄色大片| 色5月婷婷丁香| 精品人妻1区二区| 亚洲av第一区精品v没综合| 99久久成人亚洲精品观看| 99久久久亚洲精品蜜臀av| 国产一区二区三区在线臀色熟女| 别揉我奶头 嗯啊视频| 1000部很黄的大片| 亚洲国产精品成人综合色| 中文在线观看免费www的网站| 三级国产精品欧美在线观看| 精品免费久久久久久久清纯| 精品久久久久久久久久免费视频| 97人妻精品一区二区三区麻豆| 性插视频无遮挡在线免费观看| 观看免费一级毛片| 国产三级中文精品| 亚洲av熟女| 亚洲人成网站高清观看| 91av网一区二区| 久久久久九九精品影院| 久久精品国产亚洲av涩爱 | 午夜免费男女啪啪视频观看 | 精品人妻1区二区| 一个人观看的视频www高清免费观看| 亚洲av成人精品一区久久| 嫩草影视91久久| 国产黄a三级三级三级人| 首页视频小说图片口味搜索| 在线观看66精品国产| 天堂动漫精品| 色视频www国产| 精品福利观看| 久久久精品欧美日韩精品| 国产成人啪精品午夜网站| 精品一区二区三区视频在线观看免费| 一区二区三区四区激情视频 | www日本黄色视频网| 久久婷婷人人爽人人干人人爱| 人人妻人人看人人澡| 国产探花在线观看一区二区| 亚洲午夜理论影院| 亚洲国产精品合色在线| 欧美+亚洲+日韩+国产| 波多野结衣高清无吗| 嫩草影视91久久| 欧美又色又爽又黄视频| 亚洲精品亚洲一区二区| 亚洲精品一区av在线观看| 身体一侧抽搐| 波多野结衣高清无吗| 久久婷婷人人爽人人干人人爱| 欧美日韩中文字幕国产精品一区二区三区| 美女xxoo啪啪120秒动态图 | 小蜜桃在线观看免费完整版高清| 国产精品一区二区三区四区免费观看 | 校园春色视频在线观看| 老熟妇仑乱视频hdxx| 99国产精品一区二区蜜桃av| 午夜视频国产福利| 美女黄网站色视频| 内地一区二区视频在线| 99国产极品粉嫩在线观看| 一边摸一边抽搐一进一小说| 成人高潮视频无遮挡免费网站| 丰满的人妻完整版| 色综合欧美亚洲国产小说| 欧美中文日本在线观看视频| 亚洲avbb在线观看| 色5月婷婷丁香| 亚洲自拍偷在线| 亚洲五月天丁香| 国产又黄又爽又无遮挡在线| 亚洲国产高清在线一区二区三| 十八禁网站免费在线| 能在线免费观看的黄片| 免费黄网站久久成人精品 | 国产一区二区三区视频了| 久久人人爽人人爽人人片va | 国产探花在线观看一区二区| 午夜两性在线视频| 免费看a级黄色片| 欧美性感艳星| 亚洲国产精品久久男人天堂| 国产乱人伦免费视频| 久久精品国产清高在天天线| 久久九九热精品免费| 免费在线观看日本一区| 国产精品影院久久| 欧美一区二区亚洲| 亚洲人成伊人成综合网2020| 国产伦一二天堂av在线观看| 日韩欧美国产一区二区入口| 国产老妇女一区| 日韩精品青青久久久久久| 国语自产精品视频在线第100页| av欧美777| 长腿黑丝高跟| 国产伦人伦偷精品视频| 日本黄色片子视频| 亚洲va日本ⅴa欧美va伊人久久| 色哟哟哟哟哟哟| 欧美精品国产亚洲| 精品欧美国产一区二区三| 精品午夜福利视频在线观看一区| 国产探花极品一区二区| h日本视频在线播放| 我的老师免费观看完整版| 亚洲,欧美精品.| 久久九九热精品免费| 啦啦啦韩国在线观看视频| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久精品电影| 午夜精品久久久久久毛片777| 成人高潮视频无遮挡免费网站| 国产精品三级大全| 亚洲第一电影网av| 亚洲成人精品中文字幕电影| 99riav亚洲国产免费| www.熟女人妻精品国产| 搞女人的毛片| 亚洲国产高清在线一区二区三| 精品一区二区免费观看| 在线观看美女被高潮喷水网站 | 日本一本二区三区精品| 久久香蕉精品热| av在线观看视频网站免费| 久久99热这里只有精品18| 嫩草影院新地址| 51国产日韩欧美| 久久久精品大字幕| 国产精品1区2区在线观看.| 国产蜜桃级精品一区二区三区| 亚洲成人久久爱视频| 精品不卡国产一区二区三区| 一区二区三区四区激情视频 | 久久天躁狠狠躁夜夜2o2o| 在线观看舔阴道视频| 日本精品一区二区三区蜜桃| 精品人妻一区二区三区麻豆 | 99热这里只有精品一区| 精华霜和精华液先用哪个| 亚洲成人久久爱视频| 日本黄色视频三级网站网址| 18美女黄网站色大片免费观看| 成人永久免费在线观看视频| 亚洲狠狠婷婷综合久久图片| 国产av麻豆久久久久久久| 久久中文看片网| 内地一区二区视频在线| 日韩欧美精品v在线| 亚洲人成电影免费在线| 成人毛片a级毛片在线播放| 亚洲欧美清纯卡通| 午夜两性在线视频| 性插视频无遮挡在线免费观看| 久久久久久久久久成人| 国产精品精品国产色婷婷| 丁香六月欧美| 波多野结衣高清作品| 别揉我奶头 嗯啊视频| 亚洲成人中文字幕在线播放| 中出人妻视频一区二区| 美女免费视频网站| 国产精品亚洲av一区麻豆| 又粗又爽又猛毛片免费看| 国产亚洲av嫩草精品影院| 国产白丝娇喘喷水9色精品| 国产黄色小视频在线观看| 91在线观看av| 国产一区二区三区在线臀色熟女| 99久久99久久久精品蜜桃| 亚洲熟妇中文字幕五十中出| 啦啦啦观看免费观看视频高清| 国产免费av片在线观看野外av| 老司机午夜福利在线观看视频| 欧美又色又爽又黄视频| 欧美色欧美亚洲另类二区| 大型黄色视频在线免费观看| 亚洲最大成人手机在线| 身体一侧抽搐| 简卡轻食公司| 亚洲一区二区三区色噜噜| 俺也久久电影网| 在线a可以看的网站| av视频在线观看入口| 国产亚洲精品综合一区在线观看| 欧美bdsm另类| 18禁裸乳无遮挡免费网站照片| 99精品在免费线老司机午夜| 国产极品精品免费视频能看的| 最新中文字幕久久久久| www.www免费av| 看片在线看免费视频| 久久精品国产亚洲av涩爱 | 日本精品一区二区三区蜜桃| 青草久久国产| 精品一区二区三区av网在线观看| 国产高清三级在线| 动漫黄色视频在线观看| 波多野结衣高清无吗| 变态另类丝袜制服| 很黄的视频免费| 三级男女做爰猛烈吃奶摸视频| 国产一区二区激情短视频| 蜜桃亚洲精品一区二区三区| 亚洲第一区二区三区不卡| 噜噜噜噜噜久久久久久91| 99久国产av精品| 国产午夜福利久久久久久| 长腿黑丝高跟| 一进一出抽搐动态| 麻豆一二三区av精品| 国产真实伦视频高清在线观看 | 俺也久久电影网| 精品国内亚洲2022精品成人| 少妇被粗大猛烈的视频| 嫁个100分男人电影在线观看| 国产亚洲精品av在线| 欧美三级亚洲精品| 亚洲 国产 在线| 99久久久亚洲精品蜜臀av| 国产黄a三级三级三级人| 久久99热6这里只有精品| 伦理电影大哥的女人| 色吧在线观看| 18禁黄网站禁片午夜丰满| 国产精品电影一区二区三区| 大型黄色视频在线免费观看| 俄罗斯特黄特色一大片| 黄色丝袜av网址大全| 国内少妇人妻偷人精品xxx网站| 听说在线观看完整版免费高清| 少妇人妻精品综合一区二区 | 99久久精品一区二区三区| 国产真实伦视频高清在线观看 | 亚洲最大成人手机在线| 高潮久久久久久久久久久不卡| 亚洲国产精品sss在线观看| 国产视频内射| eeuss影院久久| 99热这里只有精品一区| 宅男免费午夜| 欧美潮喷喷水| 欧美性感艳星| 国产精品一区二区三区四区免费观看 | 午夜福利高清视频| 午夜激情欧美在线| 亚洲国产精品999在线| 久久久国产成人精品二区| 91av网一区二区| 成人亚洲精品av一区二区| 久久精品夜夜夜夜夜久久蜜豆| 国产精品永久免费网站| 精品午夜福利在线看| 欧美一级a爱片免费观看看| 国产高清视频在线观看网站| 少妇的逼水好多| 俄罗斯特黄特色一大片| 久久久久国产精品人妻aⅴ院| 国产高清有码在线观看视频| 小说图片视频综合网站| 好男人在线观看高清免费视频| 老熟妇仑乱视频hdxx| 中文字幕精品亚洲无线码一区| 久久久久亚洲av毛片大全| 99热这里只有是精品在线观看 | av在线老鸭窝| 久久精品综合一区二区三区| 日日摸夜夜添夜夜添av毛片 | 婷婷丁香在线五月| 香蕉av资源在线| 国内精品久久久久精免费| 国产精品日韩av在线免费观看| 国产亚洲精品av在线| 搡老妇女老女人老熟妇| 免费av不卡在线播放| 久久久久国内视频| 欧美一区二区精品小视频在线| 国内精品久久久久久久电影| 91av网一区二区| 久久精品国产亚洲av天美| 国产高清激情床上av| 99久久成人亚洲精品观看| 日本精品一区二区三区蜜桃| 人人妻人人看人人澡| 观看美女的网站| 日韩精品中文字幕看吧| 美女cb高潮喷水在线观看| 99久久久亚洲精品蜜臀av| 很黄的视频免费| 亚洲成人中文字幕在线播放| 精品人妻一区二区三区麻豆 | 成人国产一区最新在线观看| 成人鲁丝片一二三区免费| avwww免费| 久久国产精品人妻蜜桃| 欧美性猛交黑人性爽| 老司机深夜福利视频在线观看| 亚洲五月天丁香| 亚洲国产高清在线一区二区三| 日韩av在线大香蕉| 国产精品亚洲一级av第二区| 欧美日韩福利视频一区二区| 久久国产乱子伦精品免费另类| 蜜桃久久精品国产亚洲av| 在线免费观看的www视频| av专区在线播放| 日韩欧美免费精品| 黄色日韩在线| 简卡轻食公司| 亚洲精品成人久久久久久| 久久伊人香网站| 长腿黑丝高跟| 熟妇人妻久久中文字幕3abv| 性插视频无遮挡在线免费观看| 亚洲第一区二区三区不卡| 热99re8久久精品国产| 成人三级黄色视频| 国产亚洲欧美98| 国产精品av视频在线免费观看| 天堂网av新在线| 最近视频中文字幕2019在线8| 国内精品一区二区在线观看| 久久九九热精品免费| 亚洲av成人不卡在线观看播放网| 中出人妻视频一区二区| 亚洲av免费在线观看| 亚洲av中文字字幕乱码综合| 亚洲欧美日韩卡通动漫| 亚洲中文字幕一区二区三区有码在线看| 美女 人体艺术 gogo| 亚洲国产精品sss在线观看| 精品人妻一区二区三区麻豆 | 全区人妻精品视频| 一个人观看的视频www高清免费观看| 成人国产综合亚洲| 亚洲狠狠婷婷综合久久图片| 别揉我奶头~嗯~啊~动态视频| 中文亚洲av片在线观看爽| 国模一区二区三区四区视频| 国产精品久久久久久久电影| 国产亚洲av嫩草精品影院| av在线蜜桃| 在线免费观看不下载黄p国产 | 韩国av一区二区三区四区| 51午夜福利影视在线观看| 国产精品国产高清国产av| 嫩草影院新地址| 99国产精品一区二区三区| 69人妻影院| 黄色日韩在线| 午夜老司机福利剧场| 亚洲成人免费电影在线观看| av视频在线观看入口| 精品无人区乱码1区二区| 久久国产乱子伦精品免费另类| 中文字幕免费在线视频6| 亚洲激情在线av| 国内精品久久久久精免费| 综合色av麻豆| 免费看光身美女| 亚洲最大成人av| 国产伦在线观看视频一区| а√天堂www在线а√下载| 欧美乱色亚洲激情| 极品教师在线免费播放| 国产伦在线观看视频一区| 一夜夜www| 国产精品一区二区三区四区免费观看 | 在线国产一区二区在线| 嫩草影院精品99| 小蜜桃在线观看免费完整版高清| a级毛片a级免费在线| 亚洲无线观看免费| 国产伦人伦偷精品视频| 日韩精品中文字幕看吧| 毛片一级片免费看久久久久 | 久久久久久久久久成人| 日本a在线网址| 白带黄色成豆腐渣| 国模一区二区三区四区视频| 亚洲中文字幕日韩| 久久午夜亚洲精品久久| 久久久久久九九精品二区国产| 午夜激情欧美在线| 2021天堂中文幕一二区在线观| 麻豆av噜噜一区二区三区| 久久国产乱子伦精品免费另类| 国产v大片淫在线免费观看| 中文字幕人成人乱码亚洲影| 别揉我奶头 嗯啊视频| 亚洲激情在线av| 国产在线精品亚洲第一网站| 国产av不卡久久| 国产av麻豆久久久久久久| 亚洲av电影在线进入| 久久久久免费精品人妻一区二区| 国产黄色小视频在线观看| 国产白丝娇喘喷水9色精品| 久久久久国内视频| 亚洲一区二区三区色噜噜| 亚洲成人久久爱视频| 人人妻人人澡欧美一区二区| 最新在线观看一区二区三区| 婷婷色综合大香蕉| 99久久成人亚洲精品观看| 亚洲自偷自拍三级| 国产激情偷乱视频一区二区| 久久久精品欧美日韩精品| 日韩国内少妇激情av| 中亚洲国语对白在线视频| 亚洲第一区二区三区不卡| 日韩 亚洲 欧美在线| 亚洲成人精品中文字幕电影| 色av中文字幕| 色吧在线观看| 啦啦啦观看免费观看视频高清| 国产三级黄色录像| 91在线观看av| 欧美激情在线99| 在线天堂最新版资源| 国产欧美日韩一区二区精品| 午夜福利在线在线| 国产v大片淫在线免费观看| 日韩欧美一区二区三区在线观看| 亚洲在线自拍视频| or卡值多少钱| 国产毛片a区久久久久| 午夜免费激情av| 中文字幕av成人在线电影| 两人在一起打扑克的视频| 美女免费视频网站| 尤物成人国产欧美一区二区三区| 日本精品一区二区三区蜜桃| 日日摸夜夜添夜夜添小说| 亚洲av成人av| 可以在线观看毛片的网站| 舔av片在线| 中文资源天堂在线| 亚洲黑人精品在线| 99视频精品全部免费 在线| 免费观看的影片在线观看| 欧美一区二区精品小视频在线| 最新在线观看一区二区三区| 热99在线观看视频| 久久久精品大字幕| 乱人视频在线观看| 成年版毛片免费区| 真人一进一出gif抽搐免费| 国产爱豆传媒在线观看| 最新在线观看一区二区三区| 草草在线视频免费看| 精品久久久久久久久久免费视频| 亚洲最大成人av| 亚洲av五月六月丁香网| 精品一区二区三区视频在线| 国产野战对白在线观看| 国产av一区在线观看免费| 国产三级黄色录像| 国内精品一区二区在线观看| 久久午夜亚洲精品久久| 国产高清激情床上av| 国产精品亚洲美女久久久| 欧美日韩福利视频一区二区| 国产精品一及| 欧美绝顶高潮抽搐喷水| 热99re8久久精品国产| 亚洲自偷自拍三级| 亚洲精品一卡2卡三卡4卡5卡| 桃色一区二区三区在线观看|