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

    基于溫度植被干旱指數(shù)的黑龍江省20年干旱時空特征研究

    2022-08-25 08:59:50解文歡張有智宋麗娟呂志群
    水土保持研究 2022年5期
    關(guān)鍵詞:黑龍江省植被面積

    吳 黎, 解文歡, 張有智, 宋麗娟, 李 巖, 呂志群

    (黑龍江省農(nóng)業(yè)科學(xué)院 農(nóng)業(yè)遙感與信息研究所, 哈爾濱 150086)

    農(nóng)業(yè)干旱是作物生長過程中因供水不足,阻礙作物正常生長而發(fā)生的水量供應(yīng)不平衡的現(xiàn)象[1]。黑龍江省是典型的旱作農(nóng)業(yè)區(qū),近20 a來黑龍江省旱災(zāi)占總農(nóng)業(yè)氣象災(zāi)害總面積的63%,占據(jù)省內(nèi)三大氣象災(zāi)害首位。遙感監(jiān)測技術(shù)以其宏觀、快速、動態(tài)、經(jīng)濟的優(yōu)勢,成為旱災(zāi)監(jiān)測研究的有效重要手段。眾多學(xué)者研究表明利用地面溫度和植被指數(shù)特征空間耦合而成的溫度植被干旱指數(shù)(TVDI)監(jiān)測旱情方法優(yōu)越、可行性高[2]。利用TVDI模型監(jiān)測地表旱情時,其受傳感器類型和數(shù)據(jù)空間分辨率的影響較小,因此被廣泛應(yīng)用到干旱監(jiān)測中[3-8]。

    眾多學(xué)者對黑龍江省的干旱研究主要歸納為兩類:一是基于長時間的氣象站點數(shù)據(jù)采用不同干旱方法進行時空特征、演變規(guī)律、變化趨勢、災(zāi)害評價等研究,二是基于遙感數(shù)據(jù)進行干旱監(jiān)測及預(yù)測模型分析。分析干旱時空特征研究基本都是以氣象數(shù)據(jù)為基礎(chǔ)的,其原因是氣象數(shù)據(jù)歷史時間更久,從1951年開始我國便已有氣象站觀測數(shù)據(jù),時間越長對規(guī)律和趨勢分析就越有利,這是氣象數(shù)據(jù)的優(yōu)點。但氣象站屬于點數(shù)據(jù),不能代表面上干旱情況,學(xué)者需采用一系列方法進行面數(shù)據(jù)的轉(zhuǎn)化,這樣在監(jiān)測精度上有所損失。遙感數(shù)據(jù)屬于面數(shù)據(jù),以像元為單元,可監(jiān)測每一個像元的干旱程度彌補氣象數(shù)據(jù)的不足,且MODIS數(shù)據(jù)具有時效性好、幅寬大、波段范圍廣,數(shù)據(jù)開放的優(yōu)勢,可快速實現(xiàn)全省范圍內(nèi)土壤干旱監(jiān)測,在防治干旱,指導(dǎo)精準(zhǔn)灌溉,節(jié)約水資源,降低農(nóng)民損失方面均具有重要意義。目前對黑龍江省利用遙感數(shù)據(jù)進行20 a的時空變化特征研究特別缺少,學(xué)者主要集中在對單一時刻或某一時段的干旱監(jiān)測上。針對如上問題本文基于MODIS時序數(shù)據(jù)計算逐月TVDI,累計20 a間黑龍江省干旱發(fā)生頻率及干旱強度,分析黑龍江省干旱時間特征,找到規(guī)律并同時開展空間變化特征研究,為后續(xù)作者進行黑龍江省干旱預(yù)警研究奠定基礎(chǔ),為政府抗旱救災(zāi)、制定農(nóng)業(yè)政策和管控農(nóng)業(yè)項目提供科學(xué)數(shù)據(jù)支撐。

    1 研究區(qū)概況

    黑龍江省地理位置介于121°11′—135°05′E和43°25′—53°33′N。年平均氣溫在-6~4℃,無霜期多在100~140 d,年降雨量介于400~650 mm,中部最多,東部次之,北部和西部最少,黑龍江省地勢西北高東南低,農(nóng)業(yè)地域特征是東西兩大平原、南北兩大山地。

    2 數(shù)據(jù)來源與研究方法

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

    中國科學(xué)院計算機網(wǎng)絡(luò)信息中心(http:∥www.nsdata.cn/)獲取MODIS標(biāo)準(zhǔn)產(chǎn)品中16 d合成的植被指數(shù)MOD13A2,8 d合成的地表溫度MOD11A2數(shù)據(jù)。構(gòu)建溫度植被干旱指數(shù)(TVDI)需將植被指數(shù)(MOD13A2)及地表溫度(MOD11A2)數(shù)據(jù)統(tǒng)一到相同的時間分辨率和相同的空間分辨率上,將臨近的兩個8 d的MOD11A2以最大值法(同一點兩景影像中取最大值作為新影像的相元值·)合成16 d的數(shù)據(jù)。兩種數(shù)據(jù)的空間分辨率均為1 km。數(shù)據(jù)時間從2000—2019年,每年的5月初至9月末,涵蓋農(nóng)作物的整個生育期。從(https:∥ladsweb.nascom.nasa.gov/data/search.html)網(wǎng)站獲取土地覆蓋類型產(chǎn)品數(shù)據(jù)(MCD12Q1),獲取黑龍江省耕地區(qū)域界限范圍。根據(jù)國際地圈生物圈計劃(IGBP)的全球植被分類方案,選取DN值為12的像元為耕地。統(tǒng)一全部圖幅投影,后進行圖幅拼接,再利用黑龍江省行政區(qū)劃圖進行黑龍江省遙感影像的裁剪,得到黑龍江省歸一化植被指數(shù)和地表溫度數(shù)據(jù);最后利用MCD12Q1數(shù)據(jù)提取黑龍江省耕地,對MODIS所有時序數(shù)據(jù)進行掩膜處理,獲得最終的農(nóng)作物研究區(qū)的歸一化植被指數(shù)和地表溫度時序數(shù)據(jù)集。

    從中國氣象數(shù)據(jù)共享網(wǎng)(http:∥cdc.nmic.cn/)獲取黑龍江省40個氣象站2000—2019年氣象觀測降雨量數(shù)據(jù)用于計算標(biāo)準(zhǔn)化降水指數(shù)(SPI)[9]。

    2.2 研究方法

    (1) TVDI模型法。國內(nèi)外眾多學(xué)者研究不同時間分辨率和空間尺度的植被指數(shù)(NDVI)和地表溫度(Ts)的關(guān)系表明,二者存在明顯的負(fù)相關(guān)關(guān)系。Sandholt等[10]在研究NDVI-Ts特征空間時發(fā)現(xiàn)了很多條直線,據(jù)此提出了溫度植被干旱指數(shù)(TVDI)監(jiān)測地表相對含水量狀況,計算公式為:

    (1)

    Tsmax=a1+b1×NDVI

    (2)

    Tsmin=a2+b2×NDVI

    (3)

    將公式(2—3)代入公式(1)

    (4)

    式中:Ts為地表溫度;Tsmin為相同NDVI條件下的最小地表溫度;Tsmax為相同NDVI條件下的最大地表溫度;a1,a2,b1,b2為擬合方程系數(shù)。

    TVDI值為0~1,引用作者較前TVDI等級劃分結(jié)果[11],共分5個強度其中無旱對應(yīng)的強度等級I(Intensity)為1,輕旱為2,中旱為3,重旱為4,特旱為5,相對應(yīng)的TVDI值范圍見表1;20 a總干旱強度為I20,定義I20<60為干旱頻率低強度弱,60≤I20<80為干旱發(fā)生頻率中強度中,80≤I20為干旱頻率高強度大。

    表1 溫度植被干旱指數(shù)干旱強度

    為了研究黑龍江省干旱特征,分析年平均干旱強度和年內(nèi)極端干旱(年內(nèi)月平均TVDI最大值)情況,本文將2000—2019年每年5—9月月溫度植被干旱指數(shù)數(shù)據(jù)求平均獲得年平均溫度植被干旱指數(shù)(TVDImean),其代表的是作物生育期內(nèi)平均受旱程度;一年中5—9月的每月TVDI值中最大溫度植被干旱指數(shù)為TVDImax,其代表的是作物生育期內(nèi)極端受旱情況。采用一元線性回歸分析方法研究時序TVDI數(shù)據(jù)的趨勢特征。

    (2) 標(biāo)準(zhǔn)化降水指數(shù)。降雨分布是一種偏態(tài)分布,因此直接用降雨量很難在不同時間、不同地區(qū)上進行相互比較,標(biāo)準(zhǔn)化降水指數(shù)(SPI)是表示某時段降水量出現(xiàn)概率的指標(biāo)之一,因其具有概率屬性,其能很好的比較不同區(qū)域的干旱狀況[12]。通過計數(shù)出累積降水的概率密度函數(shù)Γ,進行正態(tài)標(biāo)準(zhǔn)化處理,其計算方法是在特定時間段內(nèi)產(chǎn)生的降水量與平均值之差除以標(biāo)準(zhǔn)差,其中降水量平均值和標(biāo)準(zhǔn)差是根據(jù)過去的歷年氣象觀測記錄來確定的[13]。公式如下:

    式中:x為降水樣本值;H(x)為與Γ函數(shù)相關(guān)的降水分布累積概率;c0=2.515 517;c1=0.802 853;c2=0.010 328;d1=1.432 788;d2=0.189 269;d3=0.001 308。

    本文計算的是1個月時間尺度SPI下黑龍江各站點的SPI指數(shù)值。

    3 結(jié)果與分析

    3.1 精度驗證

    標(biāo)準(zhǔn)化降水指數(shù)(SPI)是氣象干旱監(jiān)測中被廣泛應(yīng)用的指標(biāo),本研究利用SPI與TVDI進行相互精度驗證,可明確遙感TVDI干旱指標(biāo)的指導(dǎo)意義[14-16]。SPI選取2000—2019年黑龍江省40個氣象站的5—9月的月降水量數(shù)據(jù),計算1個月尺度SPI值,求40個站點平均值代表全省1個月尺度SPI情況,最后將1個月SPI值累積求平均獲得年SPI值。同時獲取與SPI時間相匹配的月TVDI值,將月TVDI值累積求平均獲得年平均TVDI(TVDImean)。對年SPI與TVDImean進行關(guān)系分析,可見年平均SPI與TVDImean有負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)為-0.67,公式為:SPI=-3.2737TVDI+1.9136,p值為0.002,表明年平均SPI與TVDImean存在極顯著負(fù)相關(guān)性。說明TVDImean可以很好的體現(xiàn)年干旱情況,即TVDImean越高,年平均標(biāo)準(zhǔn)化降水指數(shù)越小,越干旱,且氣象干旱也是造成農(nóng)業(yè)干旱的主要原因。

    3.2 黑龍江省干旱時間特征分析

    研究全省20 a來干旱時間特征變化,TVDImean時序數(shù)據(jù)(圖1)顯示,黑龍江省2006年是省內(nèi)干旱趨勢變化轉(zhuǎn)折年,干旱發(fā)生較重的在2000—2006年,此時間段內(nèi)TVDImean均值為0.63;到2013年干旱水平降到最低,2007—2013年此時間段內(nèi)TVDImean均值為0.57;2015年TVDImean值雖有所升高,但2014—2019年時間段內(nèi)TVDImean均值依然有所降低,為0.56。3個時段內(nèi)TVDImean均值呈下降趨勢。對2000—2019年TVDImean值與時間進行回歸分析,表明TVDImean值20 a間總體呈減弱的趨勢,經(jīng)過顯著性分析得出p<0.05,表明黑龍江省年平均干旱有顯著減弱趨勢。

    TVDImax監(jiān)測(圖1)出黑龍江省2006年、2009年和2015年極端干旱情況最重(TVDImax最大),且分布在3個不同干旱時間內(nèi),TVDImean與TVDImax關(guān)系不大,從變化趨勢來看黑龍江省近20 a的TVDImax有減弱的趨勢,但未通過顯著性分析。

    圖1 2000-2019年TVDI變化特征

    分析2000—2019年TVDImean中干旱面積占比情況,由圖2可以看出,黑龍江省年干旱類型中,無旱與特旱面積占比較少,主要集中在輕旱、中旱和重上,其中中旱面積占比最大,表明黑龍江省年干旱主要以中旱為主。輕旱和無旱面積占比呈增加趨勢,其他干旱等級面積占比呈減少趨勢,該趨勢與TVDImean呈減弱趨勢相一致。

    圖2 2000-2019年黑龍江省不同等級干旱面積占比

    分析5—9月20 a平均干旱面積占比的月動態(tài)情況,由圖3可知,黑龍江省無旱面積占比呈先增加后降低現(xiàn)象,在7月無旱面積占比達到峰值,8月無旱面積占比劇降其比例小于5月和6月,9月無旱面積占比降到最低;輕旱面積占比也呈先增加后降低現(xiàn)象,7月輕旱面積占比達到峰值,8月輕旱面積占比較7月減少但高于6月,9月輕旱面積占比降到最低;中旱面積占比呈先降低后升高現(xiàn)象,7月中旱面積占比降到最低,8月中旱面積占比逆轉(zhuǎn)升高比例高于5月和6月,9月中旱面積占比達到峰值;重旱面積占比6月份最高,5月重旱面積占比僅次于6月居第二,9月重旱面積占比小于6月和5月,7月和8月重旱面積占比特??;特旱面積占比呈降低現(xiàn)象5月最大,7月和8月特旱面積占比最小,幾乎接近0值,9月特旱面積占比有所增加但仍小于6月。中旱、重旱、特旱3個等級干旱面積占比總和由大到?。?月>5月>6月>8月>7月,其中雖然9月3個等級干旱面積占比總量最大,但以中旱為主,特旱特少,重旱占比也比5月和6月少,中旱本身級別不高,9月又屬秋收季節(jié),因此對作物影響不是很大;5月雖不是3個等級干旱面積占比總量最大,但其重旱、特旱面積占比均較大。因此分析表明5月是黑龍江省干旱最嚴(yán)重期,6月是干旱期,7月是黑龍江省干旱最弱時期,8月是干旱較弱時期,9月是干旱期以中旱為主。

    圖3 2000-2019年黑龍江省5-9月干旱面積占比

    3.3 黑龍江省干旱空間特征分析

    根據(jù)定義TVDImean值中各種干旱類型的干旱強度:無旱干旱強度為1,輕旱干旱強度為2,中旱干旱強度為3,重旱干旱強度為4,特旱干旱強度為5,利用公式(5)計算I20,本文假設(shè)20 a同一像元干旱強度都為4,那80就是重旱和中旱的臨界值,因此認(rèn)為80≤I20的像元點為干旱頻率高且干旱強度大的區(qū)域稱為重風(fēng)險區(qū),主要分布在訥河、甘南縣、市齊齊哈爾市、龍江縣、林甸縣、泰來縣、大慶市、肇源縣和杜爾伯特縣。假設(shè)20 a同一點干旱強度都為3,那60就是中旱和輕旱的臨界值,因此認(rèn)為60≤I20<80的像元點為干旱頻率中干旱強度中區(qū)域稱為中風(fēng)險區(qū),主要分布在西部地區(qū)、寶清縣、富錦市、友誼縣、集賢縣、勃利縣和寧安市干旱頻率中強度中。I20<60的像元點為干旱頻率輕干旱強度輕區(qū)域稱為低風(fēng)險區(qū),主要分布在北部、中部和東北部地區(qū)(圖4)。

    分析5—9月逐月干旱空間特征情況,5月重風(fēng)險區(qū)內(nèi)含縣市較多,主要分布在西部地區(qū)、寧安市、勃利縣、雞東縣、友誼縣、寶清縣和集賢縣,表明該地區(qū)每年5月份發(fā)生重旱、特旱的概率大;中風(fēng)險區(qū)主要分布在北部、中部;干旱頻率輕干旱強度輕區(qū)域主要分布在東部地區(qū)。6月重風(fēng)險區(qū)明顯減少,主要分布在西部部分地區(qū);中風(fēng)險區(qū)域增加,增加主要為5月份的重風(fēng)險區(qū),也說明6月份較5月份干旱有明顯減弱;低風(fēng)險區(qū)分布在北部和東部地區(qū)。7月、8月份重風(fēng)險區(qū)幾乎沒有,此結(jié)果與上文時間特征結(jié)果相吻合;中風(fēng)險區(qū)零星分布在西部地區(qū),該區(qū)面積8月份大于7月份;低風(fēng)險區(qū)7月份幾乎覆蓋全省(除齊齊哈爾市和大慶市附近外),8月份除西南部外的全部區(qū)域,表明7月、8月較少發(fā)生干旱。9月重風(fēng)險區(qū)主要分布在大慶市周邊及泰來縣、龍江縣、甘南縣、齊齊哈爾市、訥河市等地;中風(fēng)險區(qū)分布在西部地區(qū)和東南部地區(qū);低風(fēng)險區(qū)主要分布在北部和東北部(圖4)。

    圖4 黑龍江省20 a及逐月干旱空間特征分布

    4 討 論

    作者較前研究[11]已對溫度植被干旱指數(shù)干旱監(jiān)測指標(biāo)的干旱強度進行了劃分,且該劃分結(jié)果也在黑龍江省進行應(yīng)用與驗證[2,17]證明了劃分的TVDI干旱強度對黑龍江省土壤干旱情況的反演具有較好的適用性。本文就是基于此對黑龍江省開展2000—2019年的干旱時空變化特征研究,張劍俠等[18]以氣象數(shù)據(jù)和統(tǒng)計數(shù)據(jù)計算月干旱指數(shù),統(tǒng)計30 a(1975—2005)內(nèi)月干旱發(fā)生的次數(shù),區(qū)域劃分干旱多發(fā)區(qū)、干旱次發(fā)區(qū)和較少發(fā)生干旱區(qū),區(qū)域劃分結(jié)果與本研究結(jié)果具有一致性?;赥VDI的干旱空間特征分析是以像元為基礎(chǔ)單元,屬于面狀數(shù)據(jù),對于全省研究更具有代表性,而傳統(tǒng)的氣象數(shù)據(jù)和統(tǒng)計數(shù)據(jù)均是以氣象監(jiān)測站和縣為基礎(chǔ)單元,屬于點狀數(shù)據(jù),對全省研究屬以點概面,因此TVDI監(jiān)測更具優(yōu)勢。

    TVDI的構(gòu)建是基于溫度—植被的理想假設(shè),即在相同植被覆蓋度的前提下,地表溫度越高,土壤水分越低,這一區(qū)域就越干旱,反之就越濕潤。本文研究的黑龍江省耕地范圍生長季內(nèi)的干旱情況,僅從區(qū)域角度和時間跨度上進行分析,并未考慮到地面附著植被類型,由于植被類型的不同,植株高矮存在一定的差異,對TVDI的計算會有一定的影響,因此后續(xù)作者會從植被覆蓋類型出發(fā),研究不同植被類型、不同時期、不同區(qū)域?qū)Ω珊档捻憫?yīng)程度的研究將更有意義。

    5 結(jié) 論

    (1) 從時間尺度上來分析,2000—2019年黑龍江省20 a間2000—2006年干旱嚴(yán)重、2007—2013年干旱減弱、2014—2019年干旱震蕩減弱;TVDImax監(jiān)測出黑龍江省2006年、2009年、2015年三年中極端干旱情況最重;20 a間TVDImean有顯著減弱趨勢而TVDImax未通過顯著性檢驗。分析20 a間TVDImean中干旱面積占比可知,全省中旱現(xiàn)象較為普遍。分析月內(nèi)干旱情況可知,5月是黑龍江省干旱發(fā)生的最嚴(yán)重時期,6月是干旱期,7月是黑龍江省干旱最弱時期,8月是干旱較弱時期,9月也是干旱期以中旱為主。

    (2) 從空間尺度上來分析,20 a重風(fēng)險區(qū)主要分布在訥河市、甘南縣、齊齊哈爾市、龍江縣、林甸縣、泰來縣、大慶市、肇源縣和杜爾伯特縣。中風(fēng)險區(qū)主要分布在西部地區(qū)、寶清縣、富錦市、友誼縣、集賢縣、勃利縣和寧安市。低風(fēng)險區(qū)主要分布在北部、中部和東北部地區(qū)。分析逐月干旱空間特征:西部地區(qū)、寧安市、勃利縣、雞東縣、友誼縣、寶清縣和集賢縣等地區(qū)5月干旱發(fā)生頻率高強度大,西部部分地區(qū)6月持續(xù)高頻率高強度干旱,大慶市周邊及泰來縣、龍江縣、甘南縣、齊齊哈爾市、訥河市等地區(qū)9月干旱發(fā)生頻率高強度大。

    猜你喜歡
    黑龍江省植被面積
    基于植被復(fù)綠技術(shù)的孔植試驗及應(yīng)用
    怎樣圍面積最大
    最大的面積
    巧用面積法解幾何題
    黑龍江省節(jié)能監(jiān)測中心
    巧用面積求坐標(biāo)
    綠色植被在溯溪旅游中的應(yīng)用
    黑龍江省土壤污染防治實施方案
    黑龍江省人民政府令
    黑龍江省人民政府令
    一本一本综合久久| 国产午夜精品论理片| 国产高清videossex| 黑人欧美特级aaaaaa片| 国产1区2区3区精品| 欧洲精品卡2卡3卡4卡5卡区| 国产精华一区二区三区| 色播亚洲综合网| 国产精品av久久久久免费| 99在线人妻在线中文字幕| www日本黄色视频网| 午夜免费成人在线视频| 免费在线观看影片大全网站| 久久欧美精品欧美久久欧美| 欧美另类亚洲清纯唯美| 国产视频内射| 亚洲午夜理论影院| 亚洲真实伦在线观看| 十八禁网站免费在线| 欧美黑人巨大hd| 又爽又黄无遮挡网站| 国产激情欧美一区二区| 日韩欧美 国产精品| 香蕉av资源在线| 欧美激情在线99| 首页视频小说图片口味搜索| 国产精品99久久久久久久久| 久久国产乱子伦精品免费另类| 床上黄色一级片| 俺也久久电影网| 亚洲avbb在线观看| 老熟妇乱子伦视频在线观看| 久久这里只有精品中国| 中亚洲国语对白在线视频| 免费在线观看成人毛片| 久久人妻av系列| 国产乱人伦免费视频| 美女大奶头视频| 欧美极品一区二区三区四区| 91在线观看av| 日本五十路高清| 亚洲中文av在线| 精品欧美国产一区二区三| 国产av一区在线观看免费| 又黄又爽又免费观看的视频| 最新中文字幕久久久久 | 国产精华一区二区三区| 天堂影院成人在线观看| 白带黄色成豆腐渣| 国产高清激情床上av| 51午夜福利影视在线观看| 日本撒尿小便嘘嘘汇集6| 国产激情偷乱视频一区二区| 大型黄色视频在线免费观看| 国产精品爽爽va在线观看网站| 欧美色视频一区免费| 国产aⅴ精品一区二区三区波| 一本综合久久免费| ponron亚洲| 国产成年人精品一区二区| 亚洲精华国产精华精| 一卡2卡三卡四卡精品乱码亚洲| 91九色精品人成在线观看| 又粗又爽又猛毛片免费看| 哪里可以看免费的av片| 99在线视频只有这里精品首页| av黄色大香蕉| 欧美三级亚洲精品| 免费电影在线观看免费观看| 国产亚洲精品久久久com| 高清在线国产一区| 亚洲人成网站在线播放欧美日韩| 免费av毛片视频| 久久亚洲真实| 久久国产精品人妻蜜桃| 麻豆成人av在线观看| 亚洲片人在线观看| 嫩草影院入口| 两个人看的免费小视频| 午夜福利高清视频| 琪琪午夜伦伦电影理论片6080| 真人一进一出gif抽搐免费| 制服人妻中文乱码| 三级国产精品欧美在线观看 | 看黄色毛片网站| 两个人的视频大全免费| 五月玫瑰六月丁香| www.www免费av| 国内久久婷婷六月综合欲色啪| 老熟妇乱子伦视频在线观看| 97超视频在线观看视频| 午夜精品久久久久久毛片777| 国产又黄又爽又无遮挡在线| 每晚都被弄得嗷嗷叫到高潮| 成人特级av手机在线观看| 国产麻豆成人av免费视频| 国产激情偷乱视频一区二区| 午夜影院日韩av| 亚洲性夜色夜夜综合| xxxwww97欧美| 黄色日韩在线| 中文字幕人妻丝袜一区二区| 99久久国产精品久久久| 国产激情偷乱视频一区二区| 亚洲国产精品sss在线观看| 18禁美女被吸乳视频| 欧美黄色淫秽网站| 国产精品98久久久久久宅男小说| 五月玫瑰六月丁香| xxxwww97欧美| 在线观看舔阴道视频| 国产男靠女视频免费网站| 亚洲无线在线观看| 精品一区二区三区四区五区乱码| 好看av亚洲va欧美ⅴa在| 亚洲精品在线观看二区| 日韩欧美免费精品| 午夜福利成人在线免费观看| h日本视频在线播放| 久久精品国产亚洲av香蕉五月| 老鸭窝网址在线观看| 久久伊人香网站| 欧美日韩综合久久久久久 | 国产97色在线日韩免费| 国产乱人伦免费视频| 国产伦精品一区二区三区视频9 | 亚洲成人久久爱视频| 亚洲欧美日韩高清在线视频| 国产精品影院久久| 啪啪无遮挡十八禁网站| 啦啦啦观看免费观看视频高清| 两个人的视频大全免费| 国产免费av片在线观看野外av| 制服丝袜大香蕉在线| 国产97色在线日韩免费| 曰老女人黄片| 色综合亚洲欧美另类图片| 国产亚洲欧美在线一区二区| 99久久国产精品久久久| xxx96com| 亚洲精品在线美女| 国产精品爽爽va在线观看网站| 成人鲁丝片一二三区免费| 在线十欧美十亚洲十日本专区| 国产主播在线观看一区二区| 国产淫片久久久久久久久 | 国产麻豆成人av免费视频| 搡老熟女国产l中国老女人| 美女大奶头视频| 久久亚洲真实| 精品久久久久久久人妻蜜臀av| 婷婷精品国产亚洲av在线| 午夜亚洲福利在线播放| 亚洲av成人精品一区久久| ponron亚洲| 叶爱在线成人免费视频播放| 国产高清videossex| 身体一侧抽搐| 51午夜福利影视在线观看| 高清在线国产一区| 人人妻,人人澡人人爽秒播| 51午夜福利影视在线观看| 成人18禁在线播放| 日韩人妻高清精品专区| 免费搜索国产男女视频| 小蜜桃在线观看免费完整版高清| 亚洲真实伦在线观看| 一夜夜www| 成熟少妇高潮喷水视频| 麻豆成人av在线观看| 九色成人免费人妻av| 日本撒尿小便嘘嘘汇集6| 又黄又爽又免费观看的视频| 哪里可以看免费的av片| 精品99又大又爽又粗少妇毛片 | 1024手机看黄色片| 少妇熟女aⅴ在线视频| 国产伦精品一区二区三区视频9 | 国产亚洲欧美98| 国产精品一区二区三区四区免费观看 | 手机成人av网站| 亚洲美女视频黄频| 亚洲成人久久爱视频| 亚洲成av人片在线播放无| 丁香欧美五月| av女优亚洲男人天堂 | 欧美zozozo另类| 一本精品99久久精品77| 国产男靠女视频免费网站| 久久亚洲真实| 91在线观看av| 亚洲av美国av| 99久国产av精品| 99视频精品全部免费 在线 | 曰老女人黄片| 90打野战视频偷拍视频| av黄色大香蕉| 成人三级黄色视频| 成人性生交大片免费视频hd| 99国产精品一区二区三区| 国产私拍福利视频在线观看| 国产69精品久久久久777片 | 国产97色在线日韩免费| 九九在线视频观看精品| 日本精品一区二区三区蜜桃| 老鸭窝网址在线观看| 国产亚洲欧美98| 国产欧美日韩精品一区二区| 在线a可以看的网站| 久久国产乱子伦精品免费另类| 国产精品九九99| 性色avwww在线观看| 色综合欧美亚洲国产小说| 国产激情偷乱视频一区二区| 男人舔女人的私密视频| 久久久久久久午夜电影| 在线观看免费午夜福利视频| 美女被艹到高潮喷水动态| 欧美日韩综合久久久久久 | 久久人人精品亚洲av| 亚洲精品美女久久久久99蜜臀| 色精品久久人妻99蜜桃| 亚洲第一电影网av| 女同久久另类99精品国产91| 18禁国产床啪视频网站| 一个人观看的视频www高清免费观看 | 亚洲欧美激情综合另类| 搡老熟女国产l中国老女人| 久久久久久人人人人人| 亚洲av成人av| 一级作爱视频免费观看| 国产av麻豆久久久久久久| 国产真人三级小视频在线观看| 久久九九热精品免费| 黑人操中国人逼视频| 一进一出好大好爽视频| 国内精品久久久久久久电影| 亚洲精品在线美女| 99久久无色码亚洲精品果冻| 亚洲国产欧美网| 免费av毛片视频| 在线免费观看不下载黄p国产 | 中国美女看黄片| 中文字幕av在线有码专区| 男女做爰动态图高潮gif福利片| 欧美色视频一区免费| 日本熟妇午夜| 天堂动漫精品| 国产成人啪精品午夜网站| 一区二区三区高清视频在线| 色精品久久人妻99蜜桃| 99国产精品一区二区三区| 亚洲18禁久久av| 国产单亲对白刺激| 亚洲 国产 在线| 午夜福利在线观看吧| 久久久久久久久中文| 国产高清激情床上av| 热99re8久久精品国产| 精品福利观看| 国产精品女同一区二区软件 | 中文字幕人妻丝袜一区二区| 国产久久久一区二区三区| 男人舔女人的私密视频| 成人无遮挡网站| 在线观看舔阴道视频| 97超视频在线观看视频| www.精华液| 大型黄色视频在线免费观看| 成年人黄色毛片网站| 国产精品乱码一区二三区的特点| 亚洲av电影在线进入| 国产精品永久免费网站| 丁香欧美五月| 97超级碰碰碰精品色视频在线观看| 熟女电影av网| 国产精品,欧美在线| 免费看美女性在线毛片视频| 窝窝影院91人妻| 久久99热这里只有精品18| 精品一区二区三区四区五区乱码| 不卡一级毛片| 这个男人来自地球电影免费观看| 婷婷亚洲欧美| 精品国产乱码久久久久久男人| 亚洲人成网站高清观看| 免费av毛片视频| 成年免费大片在线观看| 男人舔女人的私密视频| 一本一本综合久久| 国产亚洲av高清不卡| 少妇裸体淫交视频免费看高清| 精品国产乱子伦一区二区三区| 色噜噜av男人的天堂激情| 黄色片一级片一级黄色片| av福利片在线观看| 最近最新中文字幕大全免费视频| 小蜜桃在线观看免费完整版高清| 亚洲精品美女久久av网站| 熟女少妇亚洲综合色aaa.| 好男人电影高清在线观看| 18禁裸乳无遮挡免费网站照片| 国产精品99久久久久久久久| 嫩草影院精品99| 免费av毛片视频| 久久精品国产99精品国产亚洲性色| 免费在线观看影片大全网站| 国产亚洲精品久久久com| 成人国产综合亚洲| 亚洲七黄色美女视频| 国产激情偷乱视频一区二区| 长腿黑丝高跟| 最新美女视频免费是黄的| 五月伊人婷婷丁香| 亚洲一区高清亚洲精品| 久久欧美精品欧美久久欧美| 国产精品美女特级片免费视频播放器 | 亚洲精品国产精品久久久不卡| 一卡2卡三卡四卡精品乱码亚洲| 欧美乱色亚洲激情| 日韩成人在线观看一区二区三区| 欧美zozozo另类| 色综合婷婷激情| 欧美黄色片欧美黄色片| 日韩欧美在线二视频| 身体一侧抽搐| 国产亚洲av嫩草精品影院| 欧美激情久久久久久爽电影| 在线观看美女被高潮喷水网站 | 亚洲av美国av| 黑人欧美特级aaaaaa片| 97超视频在线观看视频| 久久人妻av系列| 九九久久精品国产亚洲av麻豆 | 国产高清视频在线播放一区| av女优亚洲男人天堂 | 99热6这里只有精品| 日韩欧美 国产精品| 桃红色精品国产亚洲av| 国产亚洲精品久久久久久毛片| 少妇裸体淫交视频免费看高清| 国产亚洲精品综合一区在线观看| 久久久久久大精品| 一区二区三区国产精品乱码| 90打野战视频偷拍视频| 国产亚洲精品综合一区在线观看| 嫩草影院精品99| 一区二区三区国产精品乱码| a级毛片a级免费在线| 丁香六月欧美| 视频区欧美日本亚洲| 国产又黄又爽又无遮挡在线| 欧美黄色片欧美黄色片| 亚洲在线自拍视频| 波多野结衣高清作品| 亚洲精品乱码久久久v下载方式 | 中国美女看黄片| 国语自产精品视频在线第100页| 日日干狠狠操夜夜爽| 香蕉久久夜色| 51午夜福利影视在线观看| 亚洲无线观看免费| 午夜亚洲福利在线播放| www日本在线高清视频| 国产黄色小视频在线观看| 91九色精品人成在线观看| 露出奶头的视频| 国产激情偷乱视频一区二区| 国产精品爽爽va在线观看网站| 真人做人爱边吃奶动态| 国产一区二区三区在线臀色熟女| 亚洲午夜理论影院| 在线观看美女被高潮喷水网站 | 日本撒尿小便嘘嘘汇集6| 97碰自拍视频| 免费在线观看成人毛片| 久久国产乱子伦精品免费另类| 国产97色在线日韩免费| 18禁黄网站禁片免费观看直播| 欧美黄色淫秽网站| www日本黄色视频网| 99久久综合精品五月天人人| 久久精品夜夜夜夜夜久久蜜豆| 欧美一级毛片孕妇| 色噜噜av男人的天堂激情| 精品日产1卡2卡| 成人国产一区最新在线观看| 叶爱在线成人免费视频播放| 亚洲人成电影免费在线| 欧美绝顶高潮抽搐喷水| 长腿黑丝高跟| 国产人伦9x9x在线观看| 波多野结衣高清无吗| 国产成人精品久久二区二区91| 久久人妻av系列| 五月伊人婷婷丁香| 男人和女人高潮做爰伦理| 欧美乱色亚洲激情| 中文字幕熟女人妻在线| 观看免费一级毛片| 女人高潮潮喷娇喘18禁视频| www.www免费av| 亚洲av成人一区二区三| 久久婷婷人人爽人人干人人爱| 欧美日韩福利视频一区二区| 老司机在亚洲福利影院| 90打野战视频偷拍视频| av黄色大香蕉| 亚洲国产精品sss在线观看| 欧美日韩一级在线毛片| h日本视频在线播放| 亚洲性夜色夜夜综合| 欧美一级a爱片免费观看看| 国产不卡一卡二| 国产精品美女特级片免费视频播放器 | 色哟哟哟哟哟哟| 欧美激情在线99| 一区二区三区国产精品乱码| 99热只有精品国产| 校园春色视频在线观看| 国产高潮美女av| 日本免费a在线| 午夜精品一区二区三区免费看| 国产一区二区在线观看日韩 | 欧美zozozo另类| 99精品在免费线老司机午夜| 啦啦啦免费观看视频1| 欧美日韩福利视频一区二区| 亚洲专区国产一区二区| 91av网站免费观看| 窝窝影院91人妻| 亚洲激情在线av| 老司机午夜福利在线观看视频| 午夜福利在线观看免费完整高清在 | 亚洲成av人片在线播放无| 国产精品av久久久久免费| 久久精品国产综合久久久| or卡值多少钱| 久久久久国产精品人妻aⅴ院| 男女视频在线观看网站免费| 午夜福利在线观看免费完整高清在 | 亚洲片人在线观看| 欧美大码av| АⅤ资源中文在线天堂| 亚洲成av人片在线播放无| 一个人免费在线观看电影 | 性色avwww在线观看| 亚洲成人免费电影在线观看| 天天一区二区日本电影三级| 亚洲国产欧美人成| 亚洲自偷自拍图片 自拍| 国产精品久久久久久人妻精品电影| 欧美激情久久久久久爽电影| 俺也久久电影网| 成人av一区二区三区在线看| 日本撒尿小便嘘嘘汇集6| www日本黄色视频网| 久久久久精品国产欧美久久久| 五月玫瑰六月丁香| 国产亚洲精品av在线| 极品教师在线免费播放| 国语自产精品视频在线第100页| 嫁个100分男人电影在线观看| 亚洲一区二区三区色噜噜| 麻豆成人午夜福利视频| 国产三级在线视频| 蜜桃久久精品国产亚洲av| 丰满人妻一区二区三区视频av | 一个人看视频在线观看www免费 | 制服人妻中文乱码| 久久精品综合一区二区三区| 免费看光身美女| 日韩高清综合在线| 日韩人妻高清精品专区| 91麻豆av在线| 国产精品一及| 色综合婷婷激情| 国产精品自产拍在线观看55亚洲| 成年版毛片免费区| 麻豆国产av国片精品| 99国产极品粉嫩在线观看| 色av中文字幕| 国产精品一区二区三区四区久久| 亚洲av成人一区二区三| 母亲3免费完整高清在线观看| 日韩中文字幕欧美一区二区| 午夜激情欧美在线| 国产一区在线观看成人免费| 国产高清三级在线| 听说在线观看完整版免费高清| 欧美精品啪啪一区二区三区| 嫁个100分男人电影在线观看| 好男人电影高清在线观看| 亚洲第一欧美日韩一区二区三区| 在线免费观看的www视频| 99国产精品99久久久久| 国产成人精品无人区| 熟女电影av网| 日本撒尿小便嘘嘘汇集6| 久久精品夜夜夜夜夜久久蜜豆| 欧美一级a爱片免费观看看| 国产一区二区三区在线臀色熟女| 国产精品99久久99久久久不卡| 麻豆国产av国片精品| 午夜福利18| 国产精品永久免费网站| 久久精品综合一区二区三区| 成人国产综合亚洲| 99在线视频只有这里精品首页| 久久人人精品亚洲av| 国产高清videossex| 琪琪午夜伦伦电影理论片6080| 天天一区二区日本电影三级| 亚洲成av人片免费观看| 国产男靠女视频免费网站| 国产精品永久免费网站| 欧美日韩中文字幕国产精品一区二区三区| 757午夜福利合集在线观看| 青草久久国产| 一级毛片高清免费大全| 欧美在线一区亚洲| 国产精品综合久久久久久久免费| 色综合欧美亚洲国产小说| 99精品欧美一区二区三区四区| 午夜日韩欧美国产| a级毛片在线看网站| 精品电影一区二区在线| 精品日产1卡2卡| 亚洲精品美女久久久久99蜜臀| 国产伦在线观看视频一区| 叶爱在线成人免费视频播放| 亚洲午夜精品一区,二区,三区| 成年人黄色毛片网站| 亚洲乱码一区二区免费版| 亚洲第一电影网av| 久久久久久大精品| 婷婷精品国产亚洲av| 别揉我奶头~嗯~啊~动态视频| 久久国产乱子伦精品免费另类| 级片在线观看| 久久久国产精品麻豆| 精品福利观看| 国产成人aa在线观看| 99re在线观看精品视频| 欧美一级a爱片免费观看看| 免费大片18禁| 亚洲成人免费电影在线观看| 国产视频内射| 香蕉av资源在线| 国产亚洲精品综合一区在线观看| 色精品久久人妻99蜜桃| 男人舔奶头视频| 欧美国产日韩亚洲一区| 久久久久久大精品| 一二三四在线观看免费中文在| 久久久国产成人免费| 少妇丰满av| 丁香六月欧美| a级毛片在线看网站| 日韩欧美国产一区二区入口| 国产高清videossex| 听说在线观看完整版免费高清| 熟女电影av网| www国产在线视频色| 亚洲av免费在线观看| 国产欧美日韩精品一区二区| 欧美一区二区国产精品久久精品| 9191精品国产免费久久| av视频在线观看入口| 日韩欧美 国产精品| 亚洲国产中文字幕在线视频| 90打野战视频偷拍视频| 午夜两性在线视频| 国产单亲对白刺激| 变态另类丝袜制服| 国产精品一区二区精品视频观看| 亚洲人成伊人成综合网2020| 国产精品久久久久久精品电影| 欧美日本视频| 香蕉国产在线看| 中文字幕人成人乱码亚洲影| 亚洲av中文字字幕乱码综合| a级毛片在线看网站| 婷婷精品国产亚洲av在线| 最近最新中文字幕大全电影3| 俄罗斯特黄特色一大片| 国产午夜福利久久久久久| 欧美成人免费av一区二区三区| 国产一区二区激情短视频| 一级黄色大片毛片| 性欧美人与动物交配| 亚洲精品美女久久久久99蜜臀| 欧美不卡视频在线免费观看| 黄色片一级片一级黄色片| 亚洲激情在线av| 一个人看视频在线观看www免费 | 99精品在免费线老司机午夜| 欧美黑人巨大hd| 人人妻人人澡欧美一区二区| 国内精品久久久久精免费| 色播亚洲综合网| 国产av在哪里看| 亚洲欧美日韩高清在线视频| 久久亚洲真实| 欧美日本亚洲视频在线播放| 免费观看精品视频网站| 国产精品 国内视频| 日本与韩国留学比较| 久久久精品大字幕| 国产精品久久久久久人妻精品电影| 亚洲精品在线观看二区| 亚洲在线观看片| 99久久无色码亚洲精品果冻| 亚洲精品一区av在线观看| 亚洲在线自拍视频|