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

    松遼平原西部地表土壤水分時空變化與影響因素分析

    2017-12-19 02:52:28閆利爽
    關(guān)鍵詞:方根柵格土壤水分

    閆利爽,黃 方,常 帥,李 博

    (東北師范大學(xué)地理科學(xué)學(xué)院,吉林 長春 130024)

    松遼平原西部地表土壤水分時空變化與影響因素分析

    閆利爽,黃 方,常 帥,李 博

    (東北師范大學(xué)地理科學(xué)學(xué)院,吉林 長春 130024)

    基于土壤濕度和海洋鹽度衛(wèi)星(SMOS)土壤水分數(shù)據(jù),采用時空均方根誤差、趨勢分析、重心遷移模型和相關(guān)分析方法,分析了松遼平原西部2010—2014年生長季地表土壤水分時空變化特征及影響因素.結(jié)果表明:研究區(qū)地表土壤水分呈東北高、西南低的空間分布格局,5年間平均地表土壤水分含量從0.110 m3/m3增加到0.120 m3/m3.扶余、大安、通遼、奈曼旗等地表層土壤水分空間變化顯著.在不同年份,地表土壤水分隨季節(jié)變化明顯的地區(qū)不同.土壤水分重心向東北方向移動,遷移距離達30.809 km.降水、植被覆蓋因素均與地表土壤水分呈顯著正相關(guān),降水影響更強.在海拔小于200 m或坡度為10°~20°的地區(qū)地表土壤水分相對較高,平均值分別為0.140,0.116 m3/m3.黑土表層平均土壤水分更高,且具有明顯的增長趨勢.

    土壤水分;SMOS數(shù)據(jù);TRMM降水?dāng)?shù)據(jù);NDVI;松遼平原西部

    土壤水分是控制陸地和大氣間物質(zhì)、能量交換的關(guān)鍵參數(shù),對陸地表面蒸散發(fā)、水分運移和碳收支有很強的控制作用,為植被和土壤生物生存提供物質(zhì)基礎(chǔ).[1-3]表層土壤水分可能通過陸地-地表大氣的反饋機制來調(diào)控地表水分布過程,對土壤降雨入滲、產(chǎn)流和抗侵蝕能力產(chǎn)生影響,進而對深層土壤水分的動態(tài)變化起決定作用.[4]因此,及時準確地獲取與評估區(qū)域地表土壤水分時空分布規(guī)律,對于水循環(huán)和水量平衡、生態(tài)系統(tǒng)能量與物質(zhì)交換以及氣候變化等研究具有重要意義.

    實地土壤濕度測量耗時費力,時間間隔長、稀少且不連續(xù),難以建立大范圍多層深度的土壤水分監(jiān)測網(wǎng).[5]衛(wèi)星遙感為大范圍土壤水分動態(tài)監(jiān)測提供了有效途徑.土壤濕度和海洋鹽度衛(wèi)星(soil moisture and ocean salinity,SMOS)是首顆以監(jiān)測全球表面土壤水分含量和海洋鹽分為目標設(shè)計的對地觀測科學(xué)試驗衛(wèi)星,[6-7]搭載L波段(1.4 GHz)的合成孔徑微波成像輻射計(microwave imaging radiometer with aperture synthesis,MIRAS),從不同入射角度(0°~55°),以全極化、二維綜合孔徑干涉方式進行探測.[8-10]SMOS衛(wèi)星空間分辨率35~50 km,時間分辨率平均為3 d,對地表(約5 cm深度)土壤水分的探測精度可達到0.04 m3/m3,適于表現(xiàn)大范圍內(nèi)土壤水分的時空變化[4].一些研究[11-13]已對SMOS數(shù)據(jù)精度及其與實測數(shù)據(jù)的相關(guān)關(guān)系進行了分析和驗證.松遼平原西部是我國重要的農(nóng)牧業(yè)生產(chǎn)基地,但區(qū)域水分空間分配不均、水資源短缺,干旱頻發(fā),生態(tài)環(huán)境脆弱.境內(nèi)農(nóng)業(yè)氣象站點少,傳統(tǒng)技術(shù)手段很難實時獲取區(qū)域土壤干濕狀況及變化態(tài)勢.以往在本區(qū)及相鄰區(qū)域開展的大尺度研究[14-15]主要是年內(nèi)尺度的土壤水分變化規(guī)律分析,本文以SMOS土壤水分數(shù)據(jù)為基礎(chǔ),分析了松遼平原西部地表土壤水分多年的時空分布變化特征,揭示了其與降水、植被、地形及土壤類型的關(guān)系,旨在為農(nóng)業(yè)干旱監(jiān)測、農(nóng)作物估產(chǎn)、農(nóng)業(yè)資源區(qū)劃研究等提供數(shù)據(jù)支持和決策參考.

    1 研究區(qū)概況

    本文的研究區(qū)地處松嫩平原和遼河平原的交界處,介于116.271°~126.092°E,41.282°~46.275°N之間,總面積約為23.8萬km2.地勢西南高東北低,西部多山、丘陵,中東部地區(qū)為沖積平原.研究區(qū)屬半濕潤、半干旱的溫帶大陸性季風(fēng)氣候,年平均降水量300 mm,年均蒸發(fā)量是降水量的5~6倍.植被包括羊草、大針茅、隱子草、小黃柳、冰草、拂子茅和馬藺等.主要土壤類型為森林土、針葉林土、栗鈣土、沼澤土、黑鈣土等.

    2 研究方法

    2.1 數(shù)據(jù)及預(yù)處理

    本文以SMOS衛(wèi)星Level-3逐月土壤水分數(shù)據(jù)為基礎(chǔ),空間分辨率0.26°×0.26°,投影為WGS_1984.時相為2010—2014年生長季(5—10月),共30幅.Level-3數(shù)據(jù)提供近地表0~5 cm的土壤水分信息,數(shù)據(jù)格式為NetCDF[9].降水?dāng)?shù)據(jù)TRMM 3B43產(chǎn)品來自美國NASA官方網(wǎng)站(http://trmm.gsfc.nasa.gov/),空間分辨率0.25°×0.25°,時間分辨率3 h,投影與SMOS數(shù)據(jù)相同.TRMM數(shù)據(jù)在ArcGIS中用克呂格插值方法進行了重采樣,使其與SMOS數(shù)據(jù)空間分辨率一致,之后再轉(zhuǎn)為月數(shù)據(jù).其他數(shù)據(jù)包括2010—2013年5—10月的10 d最大值合成SPOT-VGT NDVI產(chǎn)品數(shù)據(jù)(VGT-S10),空間分辨率為1 km×1 km,來自VITO影像中心(http://www.vito-eodata.be),采用最大值合成法[10]獲得24幅研究區(qū)生長季各月NDVI影像.DEM數(shù)據(jù)來自地理空間數(shù)據(jù)云(http://www.gscloud.cn),利用ArcGIS空間分析模塊計算得到高程、坡度和坡向.土壤類型分布從中國科學(xué)院林業(yè)土壤研究所編制的《中國東北地區(qū)土壤圖》(1:50 000)矢量化獲得.

    2.2 時空均方根誤差

    地表土壤水分的時空變化可以通過研究時段內(nèi)每個土壤水分柵格時間、空間平均相對差值的均方根誤差(root mean square error,RMSE)予以反映[4].本研究中,首先計算土壤水分(soil moisture,SM)的空間平均值和時間平均值,公式如下:

    (1)

    (2)

    (3)

    (3)

    (5)

    (6)

    最后,土壤水分的空間和時間均方根誤差RS和RT可計算如下:

    (7)

    (8)

    RS反映土壤水分變化與區(qū)域平均狀況的差異,數(shù)值高表示土壤水分空間分布差異大,反之說明對應(yīng)柵格土壤水分變化與平均狀況一致性較強.土壤水分隨季節(jié)的變化則通過時間均方根誤差來表現(xiàn),RT高的柵格表示土壤水分與該時段的柵格平均值差異大,季節(jié)變化明顯,反之說明土壤水分值接近于該位置的時間平均值,隨季節(jié)變化小.

    2.3 趨勢分析

    通過對隨時間變化的變量與時間進行一元回歸分析,得到變化斜率,進而模擬該變量在研究時段內(nèi)的變化趨勢.[16]土壤水分的變化斜率計算如下:

    (9)

    式中:n為年數(shù),i為年份序列,以2010年為1,依次類推至2014年;Mi表示第i年SMOS影像各柵格的地表土壤水分值;SM為各柵格土壤水分趨勢線的斜率,SM為正值,說明研究期間土壤水分含量總體趨勢增加,反之則趨于減少.

    2.4 重心遷移模型

    以2010—2014年研究區(qū)地表土壤水分平均值作為重心的加權(quán)系數(shù),計算各年地表土壤水分的分布重心,得到土壤水分的空間變化軌跡.土壤水分分布重心表示為

    (10)

    式中:(XM,YM)為某年土壤水分分布重心坐標;Xi,Yi為該期SMOS影像中第i個土壤水分柵格中心的經(jīng)度和緯度;Mi為該柵格土壤水分平均值;n為柵格數(shù).重心遷移距離和角度求算公式為:

    (11)

    (12)

    其中:Dt+1,t和θt+1,t為相鄰年份土壤水分重心的遷移距離和角度;(XMt,YMt)和(XMt+1,YMt+1)分別表示為t年和ti+1年區(qū)域平均地表土壤水分重心坐標.

    2.5 相關(guān)分析

    利用Pearson相關(guān)分析方法分析地表土壤水分變化對降水、植被覆蓋的響應(yīng)關(guān)系,逐柵格分別提取地表土壤水分、NDVI值和TRMM降水值,利用SPSS統(tǒng)計軟件,進行雙變量相關(guān)分析,顯著性檢驗采用t檢驗法.

    3 地表土壤水分的時空變化特征

    3.1 空間分布格局

    2010—2014年生長季期間,松遼平原西部地表土壤水分平均值在0.059~0.205 m3/m3之間變化,空間分布差異較大,表現(xiàn)為東北部高、中西部低.近5年,地表土壤水分含量呈增長態(tài)勢.2010年區(qū)域地表土壤水分含量平均值為0.110 m3/m3,2011年土壤水分下降,僅0.095 m3/m3,是研究時段內(nèi)的最低值.此后表層水分持續(xù)上升,到2013年達到最大值0.130 m3/m3.2014年土壤水分含量又下降,平均值為0.120 m3/m3.從各月變化情況看,5—7月區(qū)域平均地表土壤水分逐漸增加,7月達最大值(0.130 m3/m3).9月地表土壤水分僅為0.100 m3/m3,降至最小值,10月平均地表土壤水分略有增加(0.110 m3/m3).

    3.2 空間變化

    2010年松遼平原西部地表土壤水分空間均方根誤差RS在0.056~0.673之間變化,2011—2013年RS的最大值均在0.9以上,土壤水分變化與區(qū)域平均狀況的差異更明顯,2014年空間均方根誤差RS在0.047~0.703之間(見圖1).各年中扶余、大安、鎮(zhèn)賚、奈曼旗等地的空間均方根誤差值均較高,說明在這些地區(qū)土壤水分含量的變化與全區(qū)平均狀況存在較大差異.在林西縣、翁牛特旗、巴林右旗、科爾沁右翼中旗、赤峰和巴林左旗,各年份RS較小,說明上述地區(qū)在同一時期土壤水分與研究區(qū)平均表層土壤水分狀況一致性較強.

    圖1 松遼平原西部地表土壤水分含量的空間均方根誤差(左)和時間均方根誤差(右)分布(2014年)

    3.3 季節(jié)變化

    松遼平原西部各年土壤水分的時間均方根誤差值均低于空間均方根誤差值(見圖1),說明區(qū)域土壤水分的季節(jié)波動弱于空間變化.2011年各地土壤水分時間均方根誤差RT值差別最大,變化范圍在0.090~0.632之間;2000年RT值次之,在0.051~0.510變化;2012年各地RT差異最小.松原、扶余等地RT相對較小,表明地表土壤水分隨季節(jié)變化不明顯.若RT較大,說明在時間序列上波動明顯.各年中,地表土壤水分顯著季節(jié)性變化分別出現(xiàn)在不同地點,如巴林左旗和霍林郭勒市(2010年)、扎魯特旗和庫倫旗(2011年)、霍林郭勒市和敖漢旗等(2012年)、通遼市和庫倫旗(2013年)、奈曼旗和庫倫旗(2014年).

    3.4 空間變化趨勢和重心遷移

    地表土壤水分變化斜率SM在0.021~0.117之間變化,說明研究區(qū)土壤水分含量呈上升趨勢,并在空間上自南向北逐漸加快.土壤水分含量分布重心遷移方向是東北—東北—東北—西南,總體向東北移動,遷移距離達30.809 km(見表1).2012—2013年,土壤水分含量分布重心雖然曾移到科爾沁左翼中旗境內(nèi),但主要分布于扎魯特旗.

    表1 地表土壤水分分布重心的變化

    4 地表土壤水分變化的影響因素

    4.1 降水

    2010—2014年生長季,研究區(qū)的平均降水量呈上升趨勢,與地表土壤水分含量增加態(tài)勢基本一致.其中:2011年僅為19 mm,同年的地表土壤水分均值也最低;2012年區(qū)域降水量最高(26 mm),但之后兩年降水持續(xù)下降.由表2可見,生長季降水量與地表土壤水分呈顯著正相關(guān).2012年,土壤水分對降水的響應(yīng)明顯,相關(guān)系數(shù)達到0.513.月尺度上,7月份的相關(guān)系數(shù)更高,說明土壤水分含量對當(dāng)月降水的響應(yīng)最為明顯,由于雨季降水豐沛,對植被生長所消耗的土壤水分有所補充.

    表2 生長季降水與地表土壤水分的相關(guān)系數(shù)

    注:**在0.01水平顯著相關(guān),*在0.05水平顯著相關(guān),下同.

    4.2 植被覆蓋

    2010—2013年,生長季平均NDVI值呈增加趨勢,植被覆蓋狀況有所改善.2010年NDVI值最低,僅為0.45,隨后持續(xù)上升,到2012年達到最大值0.53;2013年NDVI值略下降,減至0.51.由表3可見,生長季地表土壤水分與NDVI值呈顯著正相關(guān),6月和7月相關(guān)性最強.除2012年,其余年份在8月份與植被覆蓋相關(guān)性均顯著.夏季區(qū)域植被覆蓋狀況最好,保水能力強.隨著植被長勢減弱,需水量開始減少,對土壤水分的影響降低,9—10月NDVI與地表土壤水分則呈負相關(guān)關(guān)系.

    表3 NDVI與地表土壤水分的相關(guān)系數(shù)

    4.3 地形

    海拔低于200 m的地區(qū)土壤水分含量平均值為0.140 m3/m3,而200~500 m的區(qū)域平均地表土壤水分值最低(0.100 m3/m3),在海拔500~1 000 m和大于1 000 m的區(qū)域,土壤水分含量平均為0.110 m3/m3.坡度10°~20°的地區(qū)平均地表土壤水分含量最高(0.116 m3/m3),小于10°區(qū)域為0.113 m3/m3,大于20°的地區(qū)為0.111 m3/m3.從地形坡向看,南坡地表土壤水分值平均為0.121 m3/m3,向北逐漸降低,東/東南坡、西/西南坡的平均土壤水分含量分別為0.113和0.115 m3/m3,北坡最低(0.108 m3/m3),這與南坡植被生長茂盛,保水能力強有一定關(guān)聯(lián).

    4.4 土壤類型

    不同土壤因質(zhì)地不同對水分的涵養(yǎng)能力存在差別.將土壤類型圖分別與地表土壤水分平均值及變化趨勢圖進行空間疊加,提取各土壤類型區(qū)的土壤水分含量值.結(jié)果表明,黑土表層平均土壤水分含量最高,為0.150 m3/m3;水稻土次之,平均土壤水分為0.137 m3/m3;栗鈣土土壤水分含量最低(0.096 m3/m3).近5年,黑土和水稻土表層土壤水分的變化斜率分別為0.073和為0.072,增長趨勢明顯.褐土的表層水分變化趨勢最小,平均SM值僅為0.035.

    5 結(jié)論

    (1) 2010—2014年生長季,松遼平原西部表層土壤水分空間分布上呈現(xiàn)東北高、西南低的特征.年際尺度上,區(qū)域平均地表土壤水分含量值略有增加.

    (2) 近5年,松遼平原西部地區(qū)地表土壤水分的季節(jié)波動弱于空間變化.空間變化顯著的地區(qū)主要包括扶余、大安、通遼、奈曼旗等縣市.松原、扶余等地時間均方根誤差相對較小,隨季節(jié)變化不顯著.在不同年份,土壤水分隨季節(jié)變化強烈的地區(qū)不盡相同.研究區(qū)地表土壤水分含量分布重心總體上向東北方向移動.

    (3) 松遼平原西部地區(qū)降水、植被覆蓋與地表土壤水分之間存在顯著的正相關(guān)關(guān)系,各年中降水量與地表土壤水分含量的相關(guān)系數(shù)相對于植被指數(shù)要更高,說明在半濕潤-半干旱地區(qū)降水對地表土壤水分的影響更大.隨著海拔的變化,土壤水分含量波動變化.地表土壤水分自南坡向北坡逐漸降低.黑土區(qū)表層平均土壤水分含量最高,且呈更明顯的增長趨勢.

    [1] 楊濤,宮輝力,李小娟,等.土壤水分遙感監(jiān)測研究進展[J].生態(tài)學(xué)報,2010,30(22):6264-6277.

    [2] 陳亮.SMOS土壤水分反演算法研究[D].北京:中國科學(xué)院遙感應(yīng)用研究所,2009.

    [3] 王樹果.微波遙感反演地表土壤水分研究——以黑河試驗數(shù)據(jù)為例[D].北京:中國科學(xué)院研究生院,2010.

    [4] 李小英,段爭虎.基于SMOS的黃土高原區(qū)域尺度表層土壤水分時空變化[J].中國沙漠,2014,34(1):133-139.

    [5] 蘭鑫宇,郭子祺,田野,等.土壤濕度遙感估算同化研究綜述[J].地球科學(xué)進展,2015,30(6):668-679.

    [6] FONT J,CAMPS A,BORGES A,et al.SMOS:the challenging sea surface salinity measurement from space[J].Proceedings of the IEEE:Satellite Remote Sensing Missions for Monitoring Water,Carbon and Global Climate Change,2010,98(5):649-665.

    [7] SILVESTRIN P,BERGER M,KERR Y H,et al.ESA’s second earth explorer opportunity mission:the soil moisture and ocean salinity mission-SMOS[J].IEEE Transactions on Geoscience and Remote Sensing Newsletter,2001,118:11-14.

    [8] CAMPS A,F(xiàn)ONT J,VALL-LLOSSERA M,et al.The WISE 2000 and 2001 field experiments in support of the SMOS mission:sea surface L-band brightness temperature observations and their application to sea surface salinity retrieval[J].IEEE Transactions on Geoscience and Remote Sensing,2004,42(4):804-823.

    [9] 李一楠,李浩,呂容川,等.SMOS在軌定標概述[J].空間電子技術(shù),2012,47(2):20-25,36.

    [10] 唐治華.國外海洋鹽度與土壤濕度探測衛(wèi)星的發(fā)展[J].航天器工程,2013,22(3):83-89.

    [11] 陳建,張韌,王輝贊,等.Matlab對基于二進制-XML混合格式的SMOS鹽度數(shù)據(jù)的提取方法與實現(xiàn)[J].海洋通報,2011,30(6):692-696.

    [12] DENTE L,SU Z B,WEN J.Validation of SMOS soil moisture products over the Maqu and Twente regions[J].Sensors,2012,12(8):9965-9986.

    [13] DALL’AMICO J T,SCHLENZ F,LOEW A,et al.The SMOS validation campaign 2010 in the upper Danube catchment:a data set for studies of soil moisture,brightness temperature,and their spatial variability over a heterogeneous land surface[J].IEEE Transactions on Geoscience and Remote Sensing,2013,51(1):364-377.

    [14] 楊玲,楊艷昭.基于TVDI 的西遼河流域土壤濕度時空格局及其影響因素[J].干旱區(qū)資源與環(huán)境,2016,30(2):76-81.

    [15] 魏寶成,銀山,賈旭,等.蒙古高原植物生長期土壤水分時空變化特征[J].干旱區(qū)研究,2016,33(3):467-475.

    [16] HOPE A,BOYNTON W,STOW D et al.Inter-annual growth dynamics of vegetation in the Kuparuk River watershed based on the normalized difference vegetation index[J].International Journal of Remote Sensing,2003,24(17):3413-3425.

    Spatio-temporalchangesofsurfacesoilmoistureandimpactfactorsintheWesternSongliaoPlain

    YAN Li-shuang,HUANG Fang,CHANG Shuai,LI Bo

    (School of Geographical Sciences,Northeast Normal University,Changchun 130024,China)

    Integrating soil moisture data from soil moisture and ocean salinity satellite (SMOS),TRMM precipitation data,SPOT-VGT NDVI data,DEM and soil map,changes of surface soil moisture and the impact factors in the Western Songliao Plain during 2010—2014 were analyzed.The spatial and temporal root mean square error (RMSE),trend analysis and distribution center migration model were used.The results indicated that surface soil moisture decreased from the northeast to the southwest spatially.The regional average surface soil moisture increased slightly by 0.01 m3/m3.The increasing trend accelerated from the south to the north.The distribution center moved towards the northeast with 30.809 km.The greater spatial RMSE of soil moisture were found in Fuyu County,Tongliao City,Daan County and Naiman Banner etc.In different years,according to the values of temporal RMSE,the significant seasonal changes of soil moisture were found in different places.The positive correlation coefficients of surface soil moisture and precipitation were greater than those of NDVI.In the regions with the elevation lower than 200 m or 10°~20° slope,the average soil moisture was 0.140 m3/m3and 0.116 m3/m3,respectively.Averaging surface soil moisture of black soil was highest with a significantly increasing trend.

    soil moisture;SMOS;TRMM;NDVI;the Western Songliao Plain

    1000-1832(2017)04-0128-06

    10.16163/j.cnki.22-1123/n.2017.04.024

    2016-09-18

    國家自然科學(xué)基金資助項目(41571405;41671379).

    閆利爽(1989—),女,碩士研究生;通訊作者:黃方(1971—),女,博士,副教授,主要從事遙感信息分析與地表參數(shù)反演研究.

    X 87學(xué)科代碼420·20

    A

    (責(zé)任編輯:方林)

    猜你喜歡
    方根柵格土壤水分
    方根拓展探究
    基于鄰域柵格篩選的點云邊緣點提取方法*
    均方根嵌入式容積粒子PHD 多目標跟蹤方法
    西藏高原土壤水分遙感監(jiān)測方法研究
    揭開心算方根之謎
    不同剖面形狀的柵格壁對柵格翼氣動特性的影響
    數(shù)學(xué)魔術(shù)
    基于CVT排布的非周期柵格密度加權(quán)陣設(shè)計
    不同覆蓋措施對棗園土壤水分和溫度的影響
    植被覆蓋區(qū)土壤水分反演研究——以北京市為例
    久久99蜜桃精品久久| 香蕉精品网在线| 看免费av毛片| 亚洲欧洲精品一区二区精品久久久 | 免费看av在线观看网站| 午夜老司机福利剧场| 久久影院123| 性少妇av在线| 亚洲成av片中文字幕在线观看 | 亚洲中文av在线| 伊人久久国产一区二区| 黄网站色视频无遮挡免费观看| 伦理电影免费视频| 午夜老司机福利剧场| 综合色丁香网| 在线观看美女被高潮喷水网站| 99热网站在线观看| 男人操女人黄网站| 亚洲激情五月婷婷啪啪| 男女边摸边吃奶| 日产精品乱码卡一卡2卡三| 一级毛片黄色毛片免费观看视频| 丝瓜视频免费看黄片| 午夜福利影视在线免费观看| 亚洲 欧美一区二区三区| www.精华液| 在线观看国产h片| 国产精品 欧美亚洲| 日本vs欧美在线观看视频| 亚洲欧洲国产日韩| 久久久精品国产亚洲av高清涩受| 精品久久蜜臀av无| 国产精品一区二区在线不卡| 日本av手机在线免费观看| 新久久久久国产一级毛片| 亚洲av综合色区一区| 日韩一区二区三区影片| 自线自在国产av| 国产av一区二区精品久久| 国产伦理片在线播放av一区| 国产在线免费精品| 亚洲人成77777在线视频| 精品国产露脸久久av麻豆| 国产av码专区亚洲av| 国产熟女欧美一区二区| 永久网站在线| 欧美xxⅹ黑人| 国产有黄有色有爽视频| 一级毛片黄色毛片免费观看视频| 五月开心婷婷网| 九草在线视频观看| 90打野战视频偷拍视频| 国产熟女欧美一区二区| 国产熟女欧美一区二区| 精品国产露脸久久av麻豆| 人人妻人人澡人人爽人人夜夜| 亚洲国产日韩一区二区| 国产又色又爽无遮挡免| 亚洲三级黄色毛片| 天天影视国产精品| 晚上一个人看的免费电影| 热re99久久精品国产66热6| 下体分泌物呈黄色| 日韩免费高清中文字幕av| 卡戴珊不雅视频在线播放| 黄色怎么调成土黄色| 欧美国产精品一级二级三级| 久久久久国产一级毛片高清牌| 美女国产视频在线观看| 99热国产这里只有精品6| 欧美日韩成人在线一区二区| 交换朋友夫妻互换小说| 亚洲av日韩在线播放| 男女边摸边吃奶| 香蕉丝袜av| 香蕉丝袜av| 欧美精品一区二区大全| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 狠狠精品人妻久久久久久综合| 少妇人妻精品综合一区二区| 亚洲一区中文字幕在线| 亚洲国产精品成人久久小说| 一级毛片我不卡| 国产免费一区二区三区四区乱码| 天美传媒精品一区二区| 亚洲国产欧美在线一区| 亚洲精品国产一区二区精华液| 男女高潮啪啪啪动态图| 男女边吃奶边做爰视频| 久久久久国产一级毛片高清牌| 91久久精品国产一区二区三区| 777米奇影视久久| av片东京热男人的天堂| 两性夫妻黄色片| 只有这里有精品99| 亚洲激情五月婷婷啪啪| 久久免费观看电影| 日本免费在线观看一区| 国产成人精品久久二区二区91 | 九九爱精品视频在线观看| 九九爱精品视频在线观看| 亚洲美女搞黄在线观看| 亚洲av在线观看美女高潮| 黑人巨大精品欧美一区二区蜜桃| 国产在线一区二区三区精| 欧美国产精品一级二级三级| 亚洲精品国产一区二区精华液| 国产激情久久老熟女| 国产高清国产精品国产三级| 有码 亚洲区| 美女xxoo啪啪120秒动态图| av天堂久久9| 成人影院久久| 欧美少妇被猛烈插入视频| 777米奇影视久久| 丰满乱子伦码专区| 制服人妻中文乱码| 国产成人a∨麻豆精品| 一本久久精品| av国产久精品久网站免费入址| 国产精品熟女久久久久浪| 丰满饥渴人妻一区二区三| 欧美成人午夜精品| 另类精品久久| 国产在线一区二区三区精| 国产淫语在线视频| 日韩在线高清观看一区二区三区| 在线精品无人区一区二区三| 十八禁网站网址无遮挡| 国产免费福利视频在线观看| 男女啪啪激烈高潮av片| 大片电影免费在线观看免费| 久久精品国产综合久久久| 春色校园在线视频观看| 免费日韩欧美在线观看| 久久久久国产精品人妻一区二区| a级毛片黄视频| 亚洲精品第二区| 精品亚洲乱码少妇综合久久| 国产又爽黄色视频| 18+在线观看网站| 成人午夜精彩视频在线观看| 18+在线观看网站| 成人免费观看视频高清| 国产高清不卡午夜福利| 男女下面插进去视频免费观看| 日韩人妻精品一区2区三区| 亚洲婷婷狠狠爱综合网| 国产极品粉嫩免费观看在线| 婷婷色综合大香蕉| 国产一区二区在线观看av| 久久久国产一区二区| 两个人看的免费小视频| 最近中文字幕高清免费大全6| 免费观看无遮挡的男女| 一本久久精品| 9热在线视频观看99| 激情视频va一区二区三区| 欧美成人精品欧美一级黄| 纯流量卡能插随身wifi吗| 国产亚洲av片在线观看秒播厂| 欧美97在线视频| 亚洲人成网站在线观看播放| 看免费成人av毛片| 免费观看在线日韩| 亚洲av欧美aⅴ国产| av.在线天堂| 国产精品香港三级国产av潘金莲 | 曰老女人黄片| 多毛熟女@视频| 极品人妻少妇av视频| 欧美精品一区二区大全| 在线天堂中文资源库| 80岁老熟妇乱子伦牲交| 香蕉精品网在线| 激情视频va一区二区三区| 少妇精品久久久久久久| 欧美变态另类bdsm刘玥| 欧美xxⅹ黑人| 久久久精品国产亚洲av高清涩受| √禁漫天堂资源中文www| 国产成人av激情在线播放| 亚洲婷婷狠狠爱综合网| 亚洲三区欧美一区| 久久久久久久大尺度免费视频| 飞空精品影院首页| 午夜影院在线不卡| 九草在线视频观看| 老女人水多毛片| 亚洲视频免费观看视频| 精品酒店卫生间| 久久精品国产亚洲av高清一级| 久久精品国产自在天天线| 国产一区亚洲一区在线观看| 欧美少妇被猛烈插入视频| 麻豆精品久久久久久蜜桃| 久久99热这里只频精品6学生| 亚洲欧洲精品一区二区精品久久久 | 王馨瑶露胸无遮挡在线观看| 人妻人人澡人人爽人人| 999精品在线视频| 免费少妇av软件| 女性生殖器流出的白浆| 午夜日本视频在线| 亚洲成国产人片在线观看| 青春草亚洲视频在线观看| 不卡av一区二区三区| 18禁动态无遮挡网站| 久久国产精品男人的天堂亚洲| 欧美+日韩+精品| 边亲边吃奶的免费视频| 亚洲综合色惰| 观看美女的网站| 国产老妇伦熟女老妇高清| 蜜桃国产av成人99| 亚洲精品美女久久久久99蜜臀 | 久久久久久久久免费视频了| 综合色丁香网| 午夜日韩欧美国产| 中文乱码字字幕精品一区二区三区| 亚洲五月色婷婷综合| 国产爽快片一区二区三区| 日韩电影二区| 大香蕉久久网| 制服丝袜香蕉在线| 大香蕉久久成人网| 亚洲国产av影院在线观看| 成年女人毛片免费观看观看9 | 18禁国产床啪视频网站| 在线天堂中文资源库| 亚洲第一区二区三区不卡| av在线观看视频网站免费| 一区二区av电影网| av一本久久久久| 大片免费播放器 马上看| a 毛片基地| 三上悠亚av全集在线观看| 美女国产高潮福利片在线看| 一区二区三区乱码不卡18| 一二三四在线观看免费中文在| 伊人久久大香线蕉亚洲五| 两个人免费观看高清视频| 成人免费观看视频高清| 日韩av不卡免费在线播放| 有码 亚洲区| 91精品国产国语对白视频| 曰老女人黄片| 天天躁夜夜躁狠狠躁躁| 少妇的逼水好多| 看十八女毛片水多多多| a级毛片在线看网站| 搡女人真爽免费视频火全软件| 多毛熟女@视频| 嫩草影院入口| 黄网站色视频无遮挡免费观看| 国产乱人偷精品视频| 亚洲欧美成人精品一区二区| 欧美日韩一区二区视频在线观看视频在线| 黄片无遮挡物在线观看| 最近最新中文字幕免费大全7| av一本久久久久| 国产在视频线精品| 成人二区视频| 99久久人妻综合| 狠狠精品人妻久久久久久综合| 丰满饥渴人妻一区二区三| 久久人人爽人人片av| 精品亚洲乱码少妇综合久久| 成人午夜精彩视频在线观看| 国产一区二区三区综合在线观看| 中文字幕最新亚洲高清| 亚洲欧美一区二区三区国产| 亚洲三级黄色毛片| 老鸭窝网址在线观看| 男人操女人黄网站| www.熟女人妻精品国产| 99九九在线精品视频| 精品国产超薄肉色丝袜足j| av在线app专区| 精品人妻熟女毛片av久久网站| 国产精品 国内视频| av片东京热男人的天堂| 亚洲精品av麻豆狂野| 一二三四在线观看免费中文在| 国产野战对白在线观看| 久久久久久久久久久免费av| 国产一区二区三区av在线| 国产亚洲欧美精品永久| 美国免费a级毛片| 亚洲国产成人一精品久久久| a级片在线免费高清观看视频| 午夜福利在线观看免费完整高清在| 十分钟在线观看高清视频www| 国产精品人妻久久久影院| 色播在线永久视频| 婷婷色av中文字幕| 成年女人在线观看亚洲视频| 香蕉丝袜av| 免费在线观看黄色视频的| 两个人看的免费小视频| 免费日韩欧美在线观看| 久久久a久久爽久久v久久| 久久久久久久久久久免费av| av在线观看视频网站免费| 亚洲精品国产一区二区精华液| 国产精品久久久久成人av| 精品亚洲成a人片在线观看| 国产一区二区三区综合在线观看| 最新的欧美精品一区二区| 中文字幕色久视频| 日日摸夜夜添夜夜爱| 国产日韩欧美视频二区| 久久精品aⅴ一区二区三区四区 | 波多野结衣一区麻豆| 久久久久久伊人网av| 毛片一级片免费看久久久久| 久久久久久免费高清国产稀缺| 巨乳人妻的诱惑在线观看| 精品国产一区二区三区久久久樱花| 亚洲天堂av无毛| 亚洲国产精品一区三区| 最近的中文字幕免费完整| 日韩一区二区视频免费看| 精品福利永久在线观看| 亚洲av男天堂| 卡戴珊不雅视频在线播放| 超碰97精品在线观看| 国产精品不卡视频一区二区| 这个男人来自地球电影免费观看 | 电影成人av| www.精华液| 侵犯人妻中文字幕一二三四区| 国产又爽黄色视频| 肉色欧美久久久久久久蜜桃| 中文字幕色久视频| 我要看黄色一级片免费的| 老女人水多毛片| 午夜福利在线观看免费完整高清在| 免费不卡的大黄色大毛片视频在线观看| 日韩免费高清中文字幕av| 老女人水多毛片| 免费高清在线观看日韩| 少妇人妻久久综合中文| 午夜免费鲁丝| 热re99久久精品国产66热6| 久久久精品免费免费高清| 纵有疾风起免费观看全集完整版| 日韩av免费高清视频| 日本欧美国产在线视频| 国产xxxxx性猛交| 精品国产一区二区三区久久久樱花| 日日撸夜夜添| 一本大道久久a久久精品| 国产成人精品一,二区| 亚洲经典国产精华液单| 26uuu在线亚洲综合色| 少妇猛男粗大的猛烈进出视频| 日韩制服丝袜自拍偷拍| 精品国产一区二区三区久久久樱花| 极品人妻少妇av视频| 9191精品国产免费久久| 成年av动漫网址| 丝袜美腿诱惑在线| 成年女人毛片免费观看观看9 | 亚洲成国产人片在线观看| 午夜91福利影院| 久久精品国产鲁丝片午夜精品| 色94色欧美一区二区| 国产精品秋霞免费鲁丝片| 精品卡一卡二卡四卡免费| 高清在线视频一区二区三区| 亚洲中文av在线| 国产欧美日韩综合在线一区二区| 久久青草综合色| 国产精品国产三级国产专区5o| 亚洲第一青青草原| 成人手机av| 亚洲情色 制服丝袜| 久久人人爽人人片av| 午夜福利视频在线观看免费| 人妻 亚洲 视频| 精品国产乱码久久久久久男人| 伊人亚洲综合成人网| 成人影院久久| 久久久久久免费高清国产稀缺| 人妻一区二区av| 91国产中文字幕| 精品少妇久久久久久888优播| 精品国产超薄肉色丝袜足j| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | a级毛片在线看网站| 青青草视频在线视频观看| 少妇 在线观看| 日韩精品免费视频一区二区三区| 久久99一区二区三区| 国产一区亚洲一区在线观看| 99久国产av精品国产电影| 精品亚洲成a人片在线观看| www.av在线官网国产| 日韩制服骚丝袜av| 中文字幕av电影在线播放| 亚洲国产欧美日韩在线播放| 欧美激情极品国产一区二区三区| 国产精品秋霞免费鲁丝片| 99久久综合免费| 久久久久久人妻| 久久国产精品大桥未久av| 在线观看国产h片| 亚洲美女搞黄在线观看| 国产av码专区亚洲av| av一本久久久久| 精品99又大又爽又粗少妇毛片| 亚洲欧美一区二区三区黑人 | 美女中出高潮动态图| a 毛片基地| 午夜福利乱码中文字幕| 日韩中字成人| 又黄又粗又硬又大视频| kizo精华| 国产精品秋霞免费鲁丝片| 97人妻天天添夜夜摸| 国产成人91sexporn| 免费在线观看视频国产中文字幕亚洲 | 少妇猛男粗大的猛烈进出视频| 国产精品麻豆人妻色哟哟久久| 汤姆久久久久久久影院中文字幕| 91精品伊人久久大香线蕉| 日韩精品免费视频一区二区三区| 黄片无遮挡物在线观看| 丰满少妇做爰视频| 国产又色又爽无遮挡免| 国产一区二区在线观看av| a级毛片黄视频| 欧美日本中文国产一区发布| 99久久中文字幕三级久久日本| 99热全是精品| 国产av精品麻豆| 国产精品一区二区在线不卡| 99久久人妻综合| 国产日韩欧美在线精品| 90打野战视频偷拍视频| 桃花免费在线播放| 三级国产精品片| 亚洲人成77777在线视频| 久久久精品94久久精品| 久久久久久久大尺度免费视频| 肉色欧美久久久久久久蜜桃| 国产精品女同一区二区软件| 亚洲av在线观看美女高潮| 男女下面插进去视频免费观看| 黄色怎么调成土黄色| 日本91视频免费播放| 久久精品亚洲av国产电影网| 好男人视频免费观看在线| 少妇的逼水好多| 极品人妻少妇av视频| 久久精品夜色国产| 国产男女超爽视频在线观看| 激情五月婷婷亚洲| 久久久久久人妻| 精品人妻偷拍中文字幕| 9色porny在线观看| 久久久久久久亚洲中文字幕| 亚洲熟女精品中文字幕| 老司机影院成人| 中文字幕人妻丝袜一区二区 | av片东京热男人的天堂| 亚洲精品第二区| 成人亚洲欧美一区二区av| 黑人欧美特级aaaaaa片| 国产成人精品一,二区| 香蕉国产在线看| 成年人午夜在线观看视频| 成年人免费黄色播放视频| 亚洲av福利一区| 男人舔女人的私密视频| 视频在线观看一区二区三区| 亚洲国产最新在线播放| 人妻一区二区av| 亚洲av福利一区| 午夜精品国产一区二区电影| 一级片免费观看大全| 亚洲国产欧美网| 伦精品一区二区三区| 成年女人毛片免费观看观看9 | 午夜老司机福利剧场| 亚洲国产看品久久| 国产麻豆69| 精品国产超薄肉色丝袜足j| 在线看a的网站| 精品少妇内射三级| 国产午夜精品一二区理论片| 成人国语在线视频| 在线亚洲精品国产二区图片欧美| 一级片免费观看大全| 美国免费a级毛片| 捣出白浆h1v1| 91久久精品国产一区二区三区| 男女午夜视频在线观看| 欧美日韩精品成人综合77777| 亚洲欧美精品综合一区二区三区 | 亚洲天堂av无毛| 亚洲四区av| 岛国毛片在线播放| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美色中文字幕在线| 日韩伦理黄色片| 国产亚洲一区二区精品| 人成视频在线观看免费观看| 午夜91福利影院| 91精品三级在线观看| 亚洲精品国产色婷婷电影| 精品一区二区三卡| xxx大片免费视频| 欧美av亚洲av综合av国产av | 国产熟女午夜一区二区三区| 丝袜脚勾引网站| 成年美女黄网站色视频大全免费| 亚洲色图综合在线观看| a级片在线免费高清观看视频| 国产男女内射视频| 丝袜喷水一区| 国产麻豆69| 999精品在线视频| 久久精品亚洲av国产电影网| 国产精品99久久99久久久不卡 | 亚洲成人av在线免费| 观看美女的网站| 亚洲精华国产精华液的使用体验| 一本色道久久久久久精品综合| 九九爱精品视频在线观看| 亚洲天堂av无毛| 丰满饥渴人妻一区二区三| 久久久久久久久久久久大奶| 美国免费a级毛片| 国产av精品麻豆| 精品少妇内射三级| 九草在线视频观看| 精品少妇久久久久久888优播| 亚洲欧美成人精品一区二区| 亚洲精品国产一区二区精华液| 亚洲人成电影观看| 中文乱码字字幕精品一区二区三区| 人体艺术视频欧美日本| 国精品久久久久久国模美| 久久亚洲国产成人精品v| 日本猛色少妇xxxxx猛交久久| 精品亚洲成国产av| 最黄视频免费看| 少妇人妻精品综合一区二区| 亚洲经典国产精华液单| 欧美精品国产亚洲| 国产亚洲一区二区精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品国产三级专区第一集| 黄色怎么调成土黄色| 久久久欧美国产精品| 97在线视频观看| 日韩一区二区视频免费看| 国产精品免费大片| 午夜福利乱码中文字幕| 日韩av免费高清视频| 日韩视频在线欧美| 亚洲精品乱久久久久久| 欧美精品高潮呻吟av久久| 七月丁香在线播放| 久久精品国产鲁丝片午夜精品| 视频区图区小说| 婷婷色av中文字幕| 午夜av观看不卡| 国产一级毛片在线| 久久精品亚洲av国产电影网| 日韩中文字幕欧美一区二区 | 欧美成人午夜精品| 超色免费av| 麻豆精品久久久久久蜜桃| 最新的欧美精品一区二区| 亚洲欧美日韩另类电影网站| 男女啪啪激烈高潮av片| 少妇人妻精品综合一区二区| 国产黄色免费在线视频| 99热全是精品| 最新的欧美精品一区二区| 国产精品国产三级国产专区5o| 最新的欧美精品一区二区| 麻豆乱淫一区二区| 国产精品偷伦视频观看了| 婷婷成人精品国产| 国产精品无大码| 久热久热在线精品观看| 最新中文字幕久久久久| 亚洲图色成人| 亚洲美女视频黄频| 精品人妻一区二区三区麻豆| 亚洲av.av天堂| 国产精品久久久av美女十八| av网站在线播放免费| 日韩精品免费视频一区二区三区| 日本猛色少妇xxxxx猛交久久| av视频免费观看在线观看| 国产亚洲最大av| 亚洲欧美成人精品一区二区| 黄片播放在线免费| 考比视频在线观看| 一级毛片 在线播放| 日韩制服丝袜自拍偷拍| 久久久久精品久久久久真实原创| 侵犯人妻中文字幕一二三四区| 久久人妻熟女aⅴ| 9热在线视频观看99| 这个男人来自地球电影免费观看 | 国产黄色视频一区二区在线观看| 成人毛片a级毛片在线播放| 日韩一区二区三区影片| 欧美日韩视频精品一区|