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

    基于改進作物散射模型的陜西楊凌區(qū)麥田土壤水分反演研究①

    2020-10-05 07:30:58蔡慶空李二俊陶亮亮潘潔晨
    土壤 2020年4期
    關鍵詞:后向散射系數(shù)冠層

    蔡慶空,李二俊,陶亮亮,潘潔晨,陳 超,王 果

    (1 河南工程學院土木工程學院,鄭州 451191;2 河南工程學院人文社會科學學院,鄭州 451191;3 南京信息工程大學地理科學學院,南京 210044)

    遙感方法作為目前最常用的地表土壤水分探測手段,為獲取大面積土壤水分時空信息提供了可能,可以快速獲取區(qū)域地表土壤水分,有效地監(jiān)測土壤水分在空間和時間上的變化,對農(nóng)情預測、田間灌溉等舉措有著較好的指導作用[1]。目前地表土壤水分的監(jiān)測方法根據(jù)波段類型分為光學遙感、微波遙感兩種反演方法。光學遙感進行土壤水分監(jiān)測主要根據(jù)土壤的光譜反射率以及植被葉片的光譜特性變化來間接獲得土壤水分值,常常會受到觀測條件的限制,云、霧等天氣現(xiàn)象對光學遙感影響較大。而微波遙感由于其波長較長,受到天氣影響較小,能夠穿透云層,實現(xiàn)全天候、全天時對地面進行觀測,同時微波能夠穿透一定厚度的地物,對土壤介電常數(shù)非常敏感,能夠獲取地表土壤水分更加全面的信息[2-3]。

    合成孔徑雷達(synthetic aperture radar,SAR)作為微波遙感最主要的反演手段之一,廣泛應用于不同地表覆蓋類型的土壤水分反演。當利用 SAR反演裸土區(qū)的土壤水分含量時,雷達信號主要受到土壤介電常數(shù)和地表粗糙度的影響[4]。而針對植被覆蓋區(qū)的土壤水分反演,還受到地表植被的影響。目前常用的植被散射模型有MIMICS 模型、水云模型等,該類模型的本質(zhì)是研究如何有效地消除植被對雷達信號的影響,然后將其轉(zhuǎn)化為裸土區(qū)土壤水分反演問題,進而估算出植被覆蓋下的土壤含水量[5-6]。

    目前對于在植被覆蓋下的土壤水分反演多是將光學和雷達遙感結合起來進行研究,發(fā)展出以雷達反演模型為主,光學數(shù)據(jù)為輔的協(xié)同模式來反演地表土壤水分。Paris[7]利用水云模型研究玉米葉片大小對其后向散射系數(shù)的影響,假設植被的生物物理特性和散射特性與單一散射體相關,研究表明其后向散射系數(shù)與地表土壤水分直接存在線性關系。Lievens 和Verhoest[8]利用ALOS/PALSAR 雷達HH 極化數(shù)據(jù),結合水云模型和有效粗糙度參數(shù)對植被覆蓋區(qū)土壤水分進行研究,最低的土壤水分反演誤差接近5.5%。Bindlish 和Barros[9]在水云模型的基礎上通過引入植被相關長度信息建立了一個半經(jīng)驗的植被散射參數(shù)化土壤水分反演算法,消除植被對后向散射系數(shù)的影響,提高反演精度。拉巴等[10]利用MODIS 衛(wèi)星數(shù)據(jù)第7 波段對水分變化較敏感的特點,構建了一種簡單實用的土壤水分反演回歸模型,為高原土壤水分遙感監(jiān)測提供了可行的方法。蔣金豹等[11]利用MIMICS模型以及RADARSAT-2 和Landsat8 遙感數(shù)據(jù)對麥田地表土壤水分進行反演,取得了較好的反演精度。劉萬俠等[12]在分析植被對微波信號影響的基礎上,結合ASAR 雷達數(shù)據(jù),利用水云模型去除植被覆蓋影響,對后向散射系數(shù)與實測土壤水分關系進行研究,表明VV 極化的擬合效果比HH 極化要好。張友靜等[13]利用ASAR-APP 和光學影像數(shù)據(jù),結合水云模型對小麥覆蓋下地表土壤水分進行反演研究,結果表明土壤水分反演R2可達90%,均方根誤差為3.83%。趙天杰等[14]在ASAR 以及PALSAR 數(shù)據(jù)的基礎上,采用MIMICS 模型對農(nóng)作物的后向散射特性進行模擬分析,結果表明同極化的模擬效果比交叉極化好。Wang 等[15]結合水云模型,利用光學和雷達數(shù)據(jù)對植被覆蓋區(qū)進行研究,對該協(xié)同反演模型的適用范圍進行了拓展。余凡和趙英時[4]將光學與雷達數(shù)據(jù)協(xié)同起來建立半經(jīng)驗耦合模型,其反演精度顯著高于MIMICS 模型單獨反演得到的結果。趙昕等[16]利用RADARSAT-2 和Landsat-8 數(shù)據(jù)結合水云模型構建植被覆蓋區(qū)土壤水分反演的半經(jīng)驗耦合模型,反演精度有了較大提高。

    本研究將一種改進的作物散射雷達模型與PROSAIL 光學模型結合起來反演麥田土壤水分,充分考慮研究區(qū)內(nèi)冬小麥的生長狀況及空間分布,利用RADARSAT-2 雷達數(shù)據(jù)及陜西楊凌區(qū)地面實測數(shù)據(jù)對模型進行構建與驗證,估算出麥田土壤水分值。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于楊凌農(nóng)業(yè)高新技術產(chǎn)業(yè)示范區(qū),介于107.8° ~ 108.3°E 和34.1° ~ 34.5°N,地理位置如圖1所示。屬東亞暖溫帶半濕潤半干旱氣候區(qū),具有春暖多風、夏熱多雨、秋涼清爽、冬寒干燥等明顯的大陸性季風氣候特征。年均氣溫12.9 ℃,無霜期211 d。年均日照時數(shù) 2 163.8 h,年總太陽輻射量114.86 kcal/cm2,年均降水量635.1 mm。

    圖1 研究區(qū)概況(A:扶風;B:杏林;C:楊凌)Fig. 1 Overview of study area(A: Fufeng; B: Xinglin; C: Yangling)

    1.2 試驗設計

    本研究主要以冬小麥為研究對象,研究區(qū)共有3個小麥試驗站點,包括楊凌、杏林及扶風巨良農(nóng)場小麥種植示范區(qū),每個觀測站點根據(jù)小麥種植面積的差異設立不同的固定觀測點,根據(jù)試驗時間的不同選取不同的固定樣點進行觀測,由于遙感影像的限制,樣點大小為30 m×30 m,在每個樣點中測得的地面實測數(shù)據(jù)包括地表土壤水分、地表粗糙度、植株高度、葉面積指數(shù)(LAI)等。LAI 主要采用LI. COR 公司生產(chǎn)的 LAI2000 植被冠層分析儀測量,測量時間主要選在當日上午7:00—9:30 以及下午16:00—18:00,避免太陽光線直射從而帶入測量誤差,每個采樣點共測10 次,取10 次的平均值為每個樣點的LAI 值。土壤水分主要采用TDR300 型土壤水分測試儀測得,每一樣點共選取并測得5 個均勻的測量值,并將它們的平均值作為該樣點最終的土壤水分值。植株高度由卷尺測量得到,測量時選取10 株左右進行測量取平均,以地表到植株自由冠層之間的高度為準。地表粗糙度的主要測量方法為剖面板法,測量中利用相機對粗糙度板進行拍照,將所測得的照片按照剖面板輪廓進行數(shù)字化處理,根據(jù)描繪的地表起伏情況計算得到地表粗糙度參數(shù)。

    1.3 衛(wèi)星數(shù)據(jù)與預處理

    Radarsat-2 雷達衛(wèi)星一顆C 波段多極化高分辨率衛(wèi)星,可以獲取多種分辨率影像,廣泛應用于海冰海洋監(jiān)測、災害管理、農(nóng)業(yè)、水文等領域。本研究選取SLC(single look complex)格式的Radarsat-2 雷達影像,過境時間為2014 年3 月29 日,中心頻率為5.405 GHz,影像分辨率約為8 m,入射角為27.8°。對Radarsat-2 遙感影像的預處理主要在Next ESA SAR Toolbox (NEST, version 5.0.16)和ENVI5.0 軟件中進行,包括輻射校正、斑噪去除和幾何校正。

    Landsat-8 光學數(shù)據(jù)是美國國家航空航天局(NASA)于2013 年2 月11 日在加州范登堡空軍基地發(fā)射的第八顆陸地衛(wèi)星,主要為地球資源、水、森林、環(huán)境以及城市規(guī)劃等領域提供數(shù)據(jù)支撐。Landsat-8衛(wèi)星上攜帶有兩個主要載荷運行陸地成像儀(operational land imager,OLI)和熱紅外傳感器(thermal infrared sensor,TIRS),其中OLI 陸地成像儀包括9 個波段,空間分辨率為30 m,成像寬幅可達185 km×185 km。本文選取與Radarsat-2 雷達影像相對應的過境時間為2014 年3 月24 日的光學影像數(shù)據(jù),空間分辨率為30 m。Landsat-8 數(shù)據(jù)的預處理主要包括輻射定標、大氣校正、幾何校正等,主要在ENVI5.0 軟件中進行。

    1.4 改進的作物散射模型

    植被覆蓋層的雷達后向散射系數(shù)主要由與植被和土壤結構相關的幾何參數(shù)以及植被和下行土壤含水量相關的介電參數(shù)所決定,因此雷達數(shù)據(jù)所接收到的冠層水分輻射主要由植被結構和與水分相關的信息組成[17]。目前對于從雷達數(shù)據(jù)中估計植被水分的研究主要以回歸模型為主,而植被結構和水分參數(shù)會受到雷達數(shù)據(jù)極化、頻率以及入射角度的影響,因此目前的這些算法往往會局限于特定的研究區(qū)域而對于不同的環(huán)境條件將不再適用[18]。

    Saatchi 和Moghaddam[19]提出一種通過估算樹木冠層水分來反演其生物量的半經(jīng)驗模型,該模型簡化成兩層森林后向散射模型并將冠層與枝干層的結構和介電參數(shù)分離開來,從而適用于多種估算樹木冠層水分和生物量的森林類型。該模型假設在估算森林冠層參數(shù)時,總的后向散射系數(shù)主要由3 個主要散射項組成:冠層直接體散射部分、冠層與地表相互耦合散射部分以及枝干與地表相互耦合散射部分。每一個散射項都與其相關的組成部分(葉、枝干、莖)的影響參數(shù)(大小、角分布、介電常數(shù)等)有著密切的關系,其總的后向散射系數(shù)可以表示為:

    從以上的表達式可以看出,該模型假設森林冠層輻射主要受到樹冠、枝干與地表之間的交叉影響,而在致密森林條件下忽略了土壤表面的直接散射部分。上式中每一項可以具體表示為:

    式中: k0表示自由波數(shù);θ 表示入射角;是水的介電常數(shù);s 表示均方根高度;Γpq表示菲涅爾反射率;cW 表示冠層含水量;tW 表示樹干含水量;βpq,c和βpq,t分別表示冠層和樹干平均衰減系數(shù);γpq,c,γpq,cg和γpq,tg分別表示植被冠層內(nèi)整體散射截面系數(shù),依賴于冠層的幾何特性,與雷達頻率及植被含水量無關。

    為了將Saatchi和Moghaddam提出的森林散射模型簡化成適用于農(nóng)作物等低矮植被,本研究根據(jù)農(nóng)作物植被冠層和莖干沒有明顯區(qū)別的特征,將植被冠層和莖干作為一層來處理,即去除森林散射模型的莖干層,只保留其冠層直接散射部分和植被冠層與地表相互耦合部分的信息,建立低矮植被后向散射模型,即作物散射模型(crop scattering model,CSM),其植被層后向散射系數(shù)可以表示為:

    在稠密森林植被覆蓋的情況下,下墊面地表土壤的后向散射部分可以忽略不計,但是本研究的主要對象是低矮植被,其受到土壤下墊面的散射往往比較強烈,其對雷達信號的影響不能忽略,因此在構建低矮植被散射模型時,必須考慮土壤表面直接的散射信息,并將其作為一個重要的參量信息,另外根據(jù)研究區(qū)實際地表覆蓋狀況,在作物散射模型中引入經(jīng)驗參數(shù),使得地表覆蓋在像元層按照一定比例完全分離出植被層和裸土層的散射貢獻,有利于在反演土壤水分時在總的后向散射系數(shù)中剝離植被層的散射貢獻信息,構建出稀疏植被覆蓋下各地物對雷達信號散射貢獻的生成模型。由此可以得出,低矮植被覆蓋區(qū)總的后向散射系數(shù)分為植被冠層的散射項與土壤表面的直接散射項,即建立改進的作物散射模型(modified CSM,MCSM),表達式為:

    將式(5)代入(7),總的后向散射系數(shù)可以改寫為:

    對于改進模型中裸土的后向散射項,選擇AIEM進行研究,因為AIEM 能夠有效地模擬隨機粗糙度地表的后向散射系數(shù),并且該模型適用范圍較廣,模擬精度較高。

    2 結果與討論

    本研究遙感影像選擇 2014 年 3 月 29 日RADARSAT-2雷達數(shù)據(jù)和2014年3月24日Landsat-8光學數(shù)據(jù),此時冬小麥對應的生長期為拔節(jié)后期。利用實測的驗證數(shù)據(jù),從定量的角度對改進模型的模擬結果進行進一步驗證,與裸土散射模型AIEM 共同模擬植被冠層和裸土層的散射貢獻,得到總的后向散射系數(shù)模擬值,將其與雷達數(shù)據(jù)實測值進行比較,同時本研究利用未改進的作物散射模型模擬出植被總的后向散射系數(shù),與改進作物散射模型進行比較,進一步驗證改進模型的模擬精度。

    首先將 2/3 的采樣點實測數(shù)據(jù)代入公式(2)、(6),分別計算出利用散射模型AIEM計算出土壤層直接后向散射系數(shù)然后通過式(8)計算得到模型總的后向散射系數(shù)并對模型進行訓練,采用Levenberg-Marquardt 非線性最小二乘法獲取模型在HH 和VV 極化下的經(jīng)驗參數(shù),如表1所示。獲取模型經(jīng)驗參數(shù)后,本研究將剩余的1/3 采樣點驗證數(shù)據(jù)代入式(8)計算得到作物散射模型與改進作物散射模型模擬的后向散射系數(shù)值,同時和雷達后向散射系數(shù)實測值構成如圖2 所示的散點圖。

    表1 改進的作物散射模型HH 和VV 極化結構參數(shù)的估計值Table 1 Structural parameter values calculated from modified crop scattering model in HH and VV polarizations

    從圖2 可以看出,未改進的作物散射模型在HH和VV 極化下模擬得到的后向散射系數(shù)與雷達實測值之間呈線性關系,而在通過引入經(jīng)驗權重和土壤直接背景散射信息對模型進行改進,改進的作物散射模型模擬精度有顯著的提高,其線性關系R2在HH 和VV 極化下都達到80% 以上,主要是由于采樣點數(shù)據(jù)充足,同時雷達數(shù)據(jù)與光學影像的獲取時間最近,匹配程度最好。另外可以看出,HH 極化下改進模型的模擬精度都比VV 極化的要高,而研究也表明冬小麥的冠層和土壤對雷達信號的影響不同,雷達波在C波段下,HH 極化的后向散射系數(shù)要比VV 極化對植被覆蓋情況更加敏感[20-21]。

    圖2 作物散射模型與改進模型模擬得到的后向散射系數(shù)與Radarsat-2 雷達影像HH 和VV 極化實測值散點圖Fig. 2 Scatterplots of backscatter coefficients calculated from crop scattering model and modified model vs observations in HH (A, C) and VV(B, D) polarizations

    為進一步驗證該改進模型的有效性,本研究結合地面實測數(shù)據(jù)以及遙感影像數(shù)據(jù),將改進的作物散射模型土壤水分反演結果與光學指數(shù)模型TVDI、簡化的MIMICS 模型進行比較分析。目前光學遙感反演地表土壤水分最廣泛應用的是溫度植被干旱指數(shù)(TVDI)[22],該指數(shù)將與土壤水分最密切相關的兩個重要參數(shù)植被指數(shù)和地表溫度聯(lián)合起來,使得聯(lián)合后的結果地表信息更加豐富,更有利于地表土壤水分的反演。根據(jù)TVDI 與土壤含水量的對應關系,只需要利用地面實測數(shù)據(jù)計算出TVDI 值與兩者的函數(shù)關系式,即可實現(xiàn)土壤水分的反演。而MIMICS 模型[5]可以用于高大植被后向散射系數(shù)的精確模擬,但是其輸入?yún)?shù)較多,無法適用于農(nóng)作物等低矮植被后向散射系數(shù)的模擬,因此本研究在MIMICS 模型的基礎上,將植被層作為一層處理,將模型中與莖干有關的散射項消除掉,使得總的后向散射系數(shù)只包括4 部分,由此可以得到簡化的MIMICS 模型。根據(jù)地面實測數(shù)據(jù)以及研究區(qū)地理特征和地物分布情況,本研究確定MIMICS 模型的輸入?yún)?shù)為:入射角28°,頻率5.405 GHz,冠層及地表溫度為15 ℃,冠層高度為1.2 m,地表土壤體積含水量為0.5 cm3/cm3,均方根高度為2 cm,相關長度為50 cm,土壤質(zhì)地主要為沙土占 57.42%,泥沙占 20.99%,黏土占20.59%,葉片長度9 cm,寬度1 cm,植被冠層含水量為0.928 kg/m2。在簡化的MIMICS 模型中代入地面實測數(shù)據(jù)以及遙感影像數(shù)據(jù)即可反演得到地表土壤水分。

    本研究將改進的作物散射模型與TVDI 光學模型以及簡化的MIMICS 雷達模型分別對植被覆蓋下地表土壤水分進行反演,反演結果如圖3 所示。

    圖3 三種模型反演得到的土壤水分結果驗證Fig. 3 Validation of soil moistures inverted by three models

    由圖3 可以得到,改進的作物散射模型要比TVDI 和簡化的MIMICS 模型好,R2達到84.3%,均方根誤差為0.028 cm3/cm3;而TVDI 的反演結果誤差較大,R2為58.7%,均方根誤差為0.045 cm3/cm3,可能是由于光學模型對于土壤水分不是很敏感,而且光學遙感的時間差異也帶來一定的誤差。簡化的MIMICS 模型反演結果比TVDI 要好一點,但是精度不高,R2為66.9%,均方根誤差為0.043 cm3/cm3,可能是由于該簡化模型比較復雜,參數(shù)比較多,同時忽略土壤層直接的散射貢獻信息,使得反演結果受到較大的影響。

    3 結論

    本文通過對應用于森林等高大植被的散射模型進行改進,提出一種適用于農(nóng)作物、草地等低矮植被的作物散射模型,結合RADARSAT-2 雷達數(shù)據(jù)與Landsat-8 光學數(shù)據(jù),反演得到植被覆蓋區(qū)地表土壤水分。由反演結果可以得出:改進的作物散射模型反演結果比TVDI 和簡化的MIMICS 模型精度高,R2達到0.843,均方根誤差為0.028 cm3/cm3,同時該模型對土壤水分比較敏感,所需參數(shù)較少,易于使用,另外該散射模型考慮土壤自身背景直接散射信息,對植被和土壤對雷達信號的影響描述較全面,適合植被覆蓋區(qū)地表土壤水分的反演。此外該改進模型可以有效地分離出植被層和裸土層對雷達信號的散射貢獻信息,可以應用于從稀疏植被覆蓋到全植被覆蓋的多種地表覆蓋類型中,同時也可以應用于多種植被覆蓋下的復雜地表土壤水分的反演,下一步將對該改進模型的適用性作進一步的研究。

    猜你喜歡
    后向散射系數(shù)冠層
    等離子體層嘶聲波對輻射帶電子投擲角散射系數(shù)的多維建模*
    物理學報(2022年22期)2022-12-05 11:16:04
    隨機波動格點方程的后向緊隨機吸引子①
    隨機Kuramoto-Sivashinsky格點方程的后向緊隨機吸引子①
    基于低空遙感的果樹冠層信息提取方法研究
    基于激光雷達的樹形靶標冠層葉面積探測模型研究
    安徽省淮南森林冠層輻射傳輸過程的特征
    大氣科學(2021年1期)2021-04-16 07:34:18
    北部灣后向散射系數(shù)的時空分布與變化分析
    無界域上非自治Navier-Stokes方程的后向緊動力學
    施氮水平對冬小麥冠層氨揮發(fā)的影響
    一類具有飽和發(fā)生率和治療的SIS傳染病模型的后向分支及動力學行為
    日韩亚洲欧美综合| 国产综合懂色| 免费看不卡的av| 欧美区成人在线视频| 2022亚洲国产成人精品| 国产精品一及| 嫩草影院新地址| 高清欧美精品videossex| 人妻系列 视频| 乱码一卡2卡4卡精品| 91狼人影院| 老司机影院毛片| 日本与韩国留学比较| 国产精品人妻久久久影院| 听说在线观看完整版免费高清| 亚洲精品国产av成人精品| 啦啦啦啦在线视频资源| 久热这里只有精品99| 久久精品国产亚洲av涩爱| 91精品一卡2卡3卡4卡| 99热这里只有是精品50| 日日啪夜夜爽| 国产v大片淫在线免费观看| 热99国产精品久久久久久7| 69av精品久久久久久| 欧美性猛交╳xxx乱大交人| 2018国产大陆天天弄谢| 91aial.com中文字幕在线观看| 国产淫语在线视频| 亚洲欧美清纯卡通| 在线观看一区二区三区| 五月开心婷婷网| 人妻系列 视频| 国产综合精华液| 亚洲精品国产色婷婷电影| 一边亲一边摸免费视频| 色哟哟·www| 精品午夜福利在线看| 七月丁香在线播放| 精品午夜福利在线看| 亚洲av男天堂| 爱豆传媒免费全集在线观看| 午夜激情久久久久久久| 国产亚洲最大av| 有码 亚洲区| 午夜激情久久久久久久| 久久久精品欧美日韩精品| 久久97久久精品| 久久精品久久久久久久性| 国产伦精品一区二区三区四那| 夜夜看夜夜爽夜夜摸| 亚洲av成人精品一二三区| 少妇裸体淫交视频免费看高清| 秋霞在线观看毛片| 久久久精品免费免费高清| 大陆偷拍与自拍| 国产一区二区在线观看日韩| 天美传媒精品一区二区| 色视频www国产| 国产午夜精品一二区理论片| 中文字幕av成人在线电影| 精品午夜福利在线看| 男插女下体视频免费在线播放| 亚洲不卡免费看| 成人高潮视频无遮挡免费网站| 2021天堂中文幕一二区在线观| 亚洲精品国产成人久久av| 水蜜桃什么品种好| 日韩一区二区视频免费看| 国产成人午夜福利电影在线观看| 成年女人看的毛片在线观看| 国产精品人妻久久久影院| 亚洲欧美清纯卡通| 国产色婷婷99| 欧美精品一区二区大全| 99久久九九国产精品国产免费| 伊人久久国产一区二区| 国国产精品蜜臀av免费| 欧美另类一区| 成人漫画全彩无遮挡| 最近中文字幕高清免费大全6| 国产在视频线精品| 婷婷色av中文字幕| 国产v大片淫在线免费观看| videos熟女内射| 香蕉精品网在线| 热99国产精品久久久久久7| 日韩成人av中文字幕在线观看| 大片电影免费在线观看免费| 哪个播放器可以免费观看大片| 免费观看av网站的网址| 亚洲成人精品中文字幕电影| 免费电影在线观看免费观看| 国产亚洲午夜精品一区二区久久 | av女优亚洲男人天堂| 国产精品伦人一区二区| 免费少妇av软件| 国产欧美亚洲国产| 国产午夜福利久久久久久| 精品少妇黑人巨大在线播放| 国产在线男女| 日韩,欧美,国产一区二区三区| 日韩视频在线欧美| 国产精品人妻久久久久久| 夜夜爽夜夜爽视频| 国产探花极品一区二区| 99精国产麻豆久久婷婷| 久久久久久九九精品二区国产| 少妇丰满av| 亚洲欧洲日产国产| 好男人视频免费观看在线| 一级黄片播放器| tube8黄色片| 日日摸夜夜添夜夜添av毛片| 欧美日韩一区二区视频在线观看视频在线 | 高清视频免费观看一区二区| 久久久久久久久久久免费av| 人妻夜夜爽99麻豆av| 少妇被粗大猛烈的视频| 国产视频首页在线观看| 一区二区三区乱码不卡18| 成人鲁丝片一二三区免费| 亚洲在线观看片| 五月开心婷婷网| 91久久精品电影网| 国产精品国产三级专区第一集| 汤姆久久久久久久影院中文字幕| 色吧在线观看| 久久精品综合一区二区三区| 亚洲高清免费不卡视频| 内地一区二区视频在线| 麻豆成人av视频| 国产极品天堂在线| 国产精品一及| 69av精品久久久久久| 久久人人爽人人片av| 日韩电影二区| 两个人的视频大全免费| 七月丁香在线播放| 亚洲精品影视一区二区三区av| 青春草国产在线视频| 国产一区有黄有色的免费视频| 久久久久久久久久久免费av| av国产免费在线观看| 国产男人的电影天堂91| 91久久精品国产一区二区三区| 嫩草影院新地址| 亚洲欧美中文字幕日韩二区| 亚洲精品456在线播放app| 中文字幕人妻熟人妻熟丝袜美| 精品视频人人做人人爽| 日韩成人伦理影院| 天堂中文最新版在线下载 | 国产精品女同一区二区软件| 国产精品成人在线| 午夜福利网站1000一区二区三区| 黄片无遮挡物在线观看| 日韩三级伦理在线观看| av在线播放精品| 亚洲精品日本国产第一区| 免费看av在线观看网站| 久久精品国产a三级三级三级| 欧美变态另类bdsm刘玥| 特级一级黄色大片| 亚洲电影在线观看av| 亚洲精品自拍成人| 日韩精品有码人妻一区| 一级毛片电影观看| 精品久久久精品久久久| 日本黄色片子视频| 日本一本二区三区精品| 一个人看的www免费观看视频| 亚洲精品亚洲一区二区| 国产一区亚洲一区在线观看| 一二三四中文在线观看免费高清| 欧美精品人与动牲交sv欧美| 校园人妻丝袜中文字幕| 国产亚洲av片在线观看秒播厂| 两个人的视频大全免费| 在线观看国产h片| 国产精品伦人一区二区| 日本-黄色视频高清免费观看| 亚洲欧洲日产国产| 男人狂女人下面高潮的视频| 插逼视频在线观看| 国产成人福利小说| 欧美日韩亚洲高清精品| 日本黄大片高清| 中国国产av一级| 国产爽快片一区二区三区| 国产熟女欧美一区二区| 搡老乐熟女国产| 亚洲av不卡在线观看| 99精国产麻豆久久婷婷| 国产精品精品国产色婷婷| 午夜精品国产一区二区电影 | 免费av毛片视频| 五月天丁香电影| 免费av观看视频| 少妇被粗大猛烈的视频| 下体分泌物呈黄色| 不卡视频在线观看欧美| xxx大片免费视频| 天堂俺去俺来也www色官网| 色综合色国产| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美日韩国产mv在线观看视频 | 国产精品一区二区性色av| 国产精品精品国产色婷婷| 欧美成人午夜免费资源| 神马国产精品三级电影在线观看| 国产免费一区二区三区四区乱码| 国产精品人妻久久久久久| 久久精品久久久久久噜噜老黄| 禁无遮挡网站| 日韩在线高清观看一区二区三区| av免费在线看不卡| 亚洲成人一二三区av| 波野结衣二区三区在线| 狂野欧美激情性xxxx在线观看| 亚洲国产成人一精品久久久| 亚洲第一区二区三区不卡| 中国国产av一级| 久久久久久久久久久丰满| 欧美精品人与动牲交sv欧美| 精品久久久久久久末码| 亚洲欧美成人精品一区二区| 国产成年人精品一区二区| 99热这里只有是精品在线观看| 一个人看的www免费观看视频| 精品久久久久久久久av| 日本三级黄在线观看| 日韩欧美 国产精品| 午夜亚洲福利在线播放| 精品视频人人做人人爽| 99热国产这里只有精品6| 日本与韩国留学比较| 亚洲av电影在线观看一区二区三区 | 你懂的网址亚洲精品在线观看| 久久久精品94久久精品| videos熟女内射| 只有这里有精品99| 2021少妇久久久久久久久久久| 国产毛片在线视频| av免费在线看不卡| 少妇猛男粗大的猛烈进出视频 | 久久99精品国语久久久| 亚洲国产欧美在线一区| 一二三四中文在线观看免费高清| 色视频www国产| 在线播放无遮挡| 边亲边吃奶的免费视频| 久久99热这里只有精品18| av在线天堂中文字幕| 午夜激情福利司机影院| 老师上课跳d突然被开到最大视频| 精品少妇黑人巨大在线播放| 国产精品成人在线| 国国产精品蜜臀av免费| 春色校园在线视频观看| 中文字幕亚洲精品专区| 高清视频免费观看一区二区| 亚洲精品,欧美精品| 99九九线精品视频在线观看视频| 亚洲天堂国产精品一区在线| 97人妻精品一区二区三区麻豆| 国产精品嫩草影院av在线观看| 欧美激情国产日韩精品一区| 亚洲国产欧美在线一区| 丝瓜视频免费看黄片| 国内少妇人妻偷人精品xxx网站| 丝袜美腿在线中文| 91狼人影院| 精品午夜福利在线看| 国产v大片淫在线免费观看| 人妻夜夜爽99麻豆av| 一二三四中文在线观看免费高清| 亚洲国产精品999| 97超视频在线观看视频| 九九爱精品视频在线观看| 久久精品人妻少妇| 国内精品宾馆在线| 亚洲最大成人中文| 91aial.com中文字幕在线观看| 久久影院123| 亚洲成人一二三区av| 三级经典国产精品| 亚洲伊人久久精品综合| 69av精品久久久久久| 在线亚洲精品国产二区图片欧美 | www.av在线官网国产| 亚洲怡红院男人天堂| 久久久久九九精品影院| 在线观看一区二区三区激情| 纵有疾风起免费观看全集完整版| 日本猛色少妇xxxxx猛交久久| 亚洲精品成人av观看孕妇| 内射极品少妇av片p| 国产一区二区在线观看日韩| 好男人视频免费观看在线| 欧美一级a爱片免费观看看| 91午夜精品亚洲一区二区三区| 青青草视频在线视频观看| 久久久久久久大尺度免费视频| 成年版毛片免费区| 国产精品一及| 最近的中文字幕免费完整| 国产毛片在线视频| 亚洲欧洲国产日韩| 性色av一级| 边亲边吃奶的免费视频| 日韩亚洲欧美综合| 日日撸夜夜添| 在线亚洲精品国产二区图片欧美 | 中文字幕制服av| 汤姆久久久久久久影院中文字幕| 亚洲人与动物交配视频| 国产永久视频网站| 精品一区二区三区视频在线| 亚洲欧美一区二区三区黑人 | 我要看日韩黄色一级片| 人人妻人人看人人澡| 在线观看三级黄色| 亚洲在线观看片| 99热网站在线观看| 97在线人人人人妻| 美女国产视频在线观看| 久久精品国产鲁丝片午夜精品| 亚洲天堂av无毛| 久久久久性生活片| 亚洲真实伦在线观看| 久久久久精品性色| 欧美精品人与动牲交sv欧美| 简卡轻食公司| 欧美精品一区二区大全| 在线免费十八禁| 全区人妻精品视频| 亚洲精品国产av成人精品| 一个人看的www免费观看视频| 69人妻影院| 色5月婷婷丁香| 久久久久久国产a免费观看| 欧美一级a爱片免费观看看| 高清日韩中文字幕在线| 国产日韩欧美亚洲二区| 最近最新中文字幕免费大全7| 在线亚洲精品国产二区图片欧美 | 久久精品国产亚洲av天美| 亚洲熟女精品中文字幕| 国产精品伦人一区二区| 精华霜和精华液先用哪个| 天堂俺去俺来也www色官网| 亚洲精品视频女| 在线免费十八禁| 在线看a的网站| 禁无遮挡网站| 国产成人精品久久久久久| 青青草视频在线视频观看| 精品亚洲乱码少妇综合久久| 又大又黄又爽视频免费| 国产久久久一区二区三区| 哪个播放器可以免费观看大片| 99热6这里只有精品| 色综合色国产| 听说在线观看完整版免费高清| 少妇的逼水好多| 精品国产乱码久久久久久小说| 毛片女人毛片| 精品熟女少妇av免费看| 五月天丁香电影| 永久网站在线| 亚洲av免费高清在线观看| 免费观看在线日韩| 精品一区二区三区视频在线| 少妇高潮的动态图| 色哟哟·www| 国产综合精华液| 日日撸夜夜添| 男人狂女人下面高潮的视频| .国产精品久久| 丰满乱子伦码专区| 观看免费一级毛片| 国产精品成人在线| 中文乱码字字幕精品一区二区三区| 成人高潮视频无遮挡免费网站| 99热网站在线观看| 熟女人妻精品中文字幕| 一个人看视频在线观看www免费| xxx大片免费视频| 久久久久性生活片| 国产黄片美女视频| 国产精品三级大全| 亚洲久久久久久中文字幕| 22中文网久久字幕| av卡一久久| 97超碰精品成人国产| 色吧在线观看| 美女高潮的动态| 国产成人精品久久久久久| 69人妻影院| 成人欧美大片| 亚洲欧美成人综合另类久久久| 亚洲人与动物交配视频| 最近中文字幕2019免费版| 欧美一区二区亚洲| 又黄又爽又刺激的免费视频.| av播播在线观看一区| 高清欧美精品videossex| av福利片在线观看| 亚洲,一卡二卡三卡| 波多野结衣巨乳人妻| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 看免费成人av毛片| av卡一久久| 99精国产麻豆久久婷婷| 国产精品熟女久久久久浪| 国产午夜福利久久久久久| 日韩亚洲欧美综合| 午夜日本视频在线| 嘟嘟电影网在线观看| 久久精品国产自在天天线| 中文字幕免费在线视频6| 日韩精品有码人妻一区| 丝袜喷水一区| 一级a做视频免费观看| 精品久久久久久久人妻蜜臀av| 噜噜噜噜噜久久久久久91| 人妻一区二区av| 天美传媒精品一区二区| 久久女婷五月综合色啪小说 | 伊人久久精品亚洲午夜| 天天躁日日操中文字幕| 熟女电影av网| 国产人妻一区二区三区在| 日韩人妻高清精品专区| 美女高潮的动态| 看非洲黑人一级黄片| 欧美高清性xxxxhd video| 久久久成人免费电影| 新久久久久国产一级毛片| 免费播放大片免费观看视频在线观看| 久久精品国产亚洲网站| 国产精品久久久久久久电影| 国产成人一区二区在线| 最近最新中文字幕免费大全7| 少妇 在线观看| 色网站视频免费| a级毛片免费高清观看在线播放| 国产大屁股一区二区在线视频| 成人黄色视频免费在线看| 欧美xxxx性猛交bbbb| 精品国产露脸久久av麻豆| 校园人妻丝袜中文字幕| a级一级毛片免费在线观看| 久久久久久伊人网av| 亚洲国产成人一精品久久久| 亚洲欧美精品专区久久| 黄色视频在线播放观看不卡| 全区人妻精品视频| 毛片一级片免费看久久久久| 天美传媒精品一区二区| 特级一级黄色大片| 亚洲精品影视一区二区三区av| 国产精品久久久久久久电影| 视频中文字幕在线观看| 成年版毛片免费区| 禁无遮挡网站| 亚洲精品乱码久久久久久按摩| 性色av一级| 日韩不卡一区二区三区视频在线| 国产色爽女视频免费观看| 亚洲欧洲日产国产| 久久久久久久久久成人| 久久韩国三级中文字幕| 永久网站在线| 国产片特级美女逼逼视频| 午夜福利在线在线| 18禁动态无遮挡网站| 高清日韩中文字幕在线| av专区在线播放| 在线观看美女被高潮喷水网站| 伦理电影大哥的女人| 高清欧美精品videossex| 天天躁夜夜躁狠狠久久av| 欧美丝袜亚洲另类| 亚洲天堂av无毛| 日日啪夜夜爽| 有码 亚洲区| 在线精品无人区一区二区三 | 国产精品一区www在线观看| 国产精品国产三级专区第一集| 欧美日韩亚洲高清精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美日韩无卡精品| 美女视频免费永久观看网站| 国产黄频视频在线观看| 久久久久国产精品人妻一区二区| 麻豆久久精品国产亚洲av| 亚洲精品国产av成人精品| 在线播放无遮挡| 亚洲国产欧美人成| 丝袜喷水一区| 99久久中文字幕三级久久日本| 国产久久久一区二区三区| 国产 一区精品| 99热这里只有是精品在线观看| 麻豆久久精品国产亚洲av| 卡戴珊不雅视频在线播放| 99re6热这里在线精品视频| 九色成人免费人妻av| 在线观看一区二区三区激情| 视频区图区小说| 蜜桃亚洲精品一区二区三区| 永久免费av网站大全| 午夜老司机福利剧场| 网址你懂的国产日韩在线| 女人十人毛片免费观看3o分钟| 尾随美女入室| 久久久久久国产a免费观看| 亚洲三级黄色毛片| 2018国产大陆天天弄谢| 免费少妇av软件| 免费不卡的大黄色大毛片视频在线观看| 国产精品蜜桃在线观看| 美女主播在线视频| 中文字幕免费在线视频6| xxx大片免费视频| 人妻一区二区av| 18禁裸乳无遮挡免费网站照片| 亚洲一区二区三区欧美精品 | 18禁裸乳无遮挡动漫免费视频 | 在线观看免费高清a一片| 制服丝袜香蕉在线| 热re99久久精品国产66热6| 精品一区二区免费观看| 亚洲人成网站在线观看播放| 天天躁日日操中文字幕| 男女边摸边吃奶| 欧美高清性xxxxhd video| 精品久久国产蜜桃| 国产精品.久久久| 男人爽女人下面视频在线观看| 免费观看a级毛片全部| 男女那种视频在线观看| 在线观看国产h片| 99精国产麻豆久久婷婷| 国产精品国产三级国产专区5o| 少妇人妻久久综合中文| 国产午夜精品一二区理论片| 亚洲高清免费不卡视频| 国产黄色免费在线视频| 国产乱人偷精品视频| av一本久久久久| 三级男女做爰猛烈吃奶摸视频| 亚洲精品成人av观看孕妇| 少妇人妻久久综合中文| 久久亚洲国产成人精品v| videossex国产| 日产精品乱码卡一卡2卡三| 新久久久久国产一级毛片| 中国国产av一级| 高清午夜精品一区二区三区| 久久精品国产亚洲网站| 最近最新中文字幕免费大全7| 人妻一区二区av| 大香蕉97超碰在线| 亚洲成人中文字幕在线播放| 国产淫语在线视频| av卡一久久| 日日啪夜夜撸| 老女人水多毛片| 深爱激情五月婷婷| 大话2 男鬼变身卡| 狠狠精品人妻久久久久久综合| 男女边吃奶边做爰视频| 成人免费观看视频高清| 亚洲欧美日韩另类电影网站 | 日韩av免费高清视频| 国产毛片a区久久久久| 亚洲成人久久爱视频| 久久99热这里只有精品18| 久久亚洲国产成人精品v| 国产亚洲av嫩草精品影院| 亚洲成人中文字幕在线播放| 秋霞伦理黄片| 国模一区二区三区四区视频| 中文欧美无线码| 精品一区二区三卡| 毛片一级片免费看久久久久| 久久久午夜欧美精品| av在线老鸭窝| 亚洲色图综合在线观看| 天天躁夜夜躁狠狠久久av| 亚洲精品一区蜜桃| av在线天堂中文字幕| 亚洲欧美一区二区三区黑人 | 麻豆成人午夜福利视频| eeuss影院久久| 精品久久久噜噜| 国产高清国产精品国产三级 | 成人亚洲欧美一区二区av| 成年人午夜在线观看视频| 免费看a级黄色片| 欧美成人午夜免费资源| 国产男人的电影天堂91| 日日摸夜夜添夜夜添av毛片| 色播亚洲综合网| 亚洲av日韩在线播放| 三级国产精品片| 又大又黄又爽视频免费| 国产成人精品一,二区| 国产高清国产精品国产三级 | 夜夜看夜夜爽夜夜摸| 97精品久久久久久久久久精品| 亚洲精品成人久久久久久|