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

    不同GRACE-FO模型反演流域陸地水儲量變化對比

    2022-04-16 10:48:06崔立魯杜安張誠汪曉龍李瓊
    人民長江 2022年3期
    關(guān)鍵詞:重力場時變方差

    崔立魯 杜安 張誠 汪曉龍 李瓊

    摘要:?為了有效比較由美國空間研究中心(Center of Space Research,CSR)、德國地學(xué)中心(Helmholtz-Centre Potsdam-German Research Centre for Geosciences,GFZ)和美國噴氣推進實驗室(Jet Propulsion Laboratory,JPL)機構(gòu)提供的3種GRACE-FO RL06時變重力場模型精度,分別從模型階方差、C? 20 項精度、全球和局部陸地水儲量變化反演等角度對3種模型進行了詳盡的闡述。結(jié)果表明:JPL的RL06模型階方差要大于其他2個機構(gòu)發(fā)布的結(jié)果;而在C? 20 項精度方面,CSR和JPL發(fā)布的數(shù)據(jù)更接近于衛(wèi)星激光測距(Satellite Laser Ranging,SLR);通過對比這3種模型反演的陸地水儲量變化結(jié)果與全球陸面數(shù)據(jù)同化系統(tǒng)(Global Land Data Assimilation system,GLDAS)水文模型的反演結(jié)果,發(fā)現(xiàn)3種模型與GLDAS模型的結(jié)果符合度都比較高,從局部反演結(jié)果來看,在3種模型中CSR與GLDAS模型的符合度較高。

    關(guān) 鍵 詞:陸地水儲量變化; 模型階方差; C? 20 項精度; 相關(guān)系數(shù); GRACE-FO模型

    中圖法分類號: ?P223

    文獻標(biāo)志碼: ?A

    DOI: 10.16232/j.cnki.1001-4179.2022.03.015

    ?0 引 言

    自從2002年3月GRACE(Gravity Recovery and Climate Experiment)衛(wèi)星計劃實施以來,該衛(wèi)星數(shù)據(jù)被廣泛地運用于全球環(huán)境變化監(jiān)測中? [1] ,包括陸地水儲量變化(Terrestrial Water Storage Change,TWSC)? [2-4] 、冰川消融? [5-7] 、海平面上升? [8] 和同震變化? [9-10] 等。該衛(wèi)星任務(wù)于2017年7月結(jié)束,為了延續(xù)該計劃,美國國家航空航天局(National Aeronautics and Space Administration,NASA)和德國宇航中心(Deutsches Zentrum für Luft-und Raumfahrt,DLR)于2018年3月聯(lián)合發(fā)射了GRACE-FO(GRACE Follow-On)衛(wèi)星。目前,美國空間研究中心(Center of Space Research,CSR)、德國地學(xué)中心(Helmholtz-Centre Potsdam-German Research Centre for Geosciences,GFZ)和美國噴氣推進實驗室(Jet Propulsion Laboratory,JPL)三大官方機構(gòu)先后發(fā)布了GRACE-FO時變重力場模型。但是由于這3種模型是由不同機構(gòu)分別解算的,因此判斷3種模型在反演陸地水儲量變化(TWSC)方面的優(yōu)劣性是水文應(yīng)用方面急需解決的問題。同時由于GRACE-FO其中一顆衛(wèi)星上的加速度計出現(xiàn)問題? [11] ,那么采用模型修復(fù)后的GRACE-FO數(shù)據(jù)是否還具備探測TWSC的能力也有待考證。

    李瓊? [12] 從模型階方差、C? 20 項精度和局部TWSC等方面比較驗證了3個機構(gòu)GRACE RL04和RL05模型;郭飛霄等? [13] 從模型階方差、濾波結(jié)果和C? 20 項精度等角度剖析了3個機構(gòu)GRACE RL05和RL06模型。上述研究均取得了一定的成果,但是目前關(guān)于GRACE-FO RL06模型的比較研究較少,同時與水文模型的比較也略顯不足。

    本文針對上述情況從模型階方差、C? 20 項精度、全球TWSC和區(qū)域水儲量變化等4個方面對三大機構(gòu)所發(fā)布的GRACE-FO RL06時變重力場模型進行詳細比較,并引入全球陸面數(shù)據(jù)同化系統(tǒng)(Global Land Data Assimilation system,GLDAS)的水文模型進行驗證。

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

    1.1 實驗數(shù)據(jù)

    本文采用由CSR、GFZ和JPL三大機構(gòu)提供的RL06版本GRACE-FO時變重力場模型數(shù)據(jù)(分別簡稱為CSR RL06,GFZ RL06和JPL RL06),截斷階數(shù)為60階,其時間跨度為2018年6月至2020年3月。3種重力場模型已扣除了由潮汐和非潮汐所引起的大氣和海洋質(zhì)量變化影響,并減去了地球長期靜態(tài)重力場的影響。

    GLDAS水文模型是由NASA與美國海洋和大氣局聯(lián)合發(fā)布的,它包含了全球(除了兩極地區(qū))降雨量、水分蒸發(fā)量、地表徑流、地下徑流、土壤濕度、地表積雪分布等數(shù)據(jù),該數(shù)據(jù)為1°×1°的格網(wǎng)數(shù)據(jù)。本文首先對其進行球諧展開,再截斷至與GRACE時變重力場模型相同階數(shù),采用與GRACE數(shù)據(jù)處理相同的濾波方法得到格網(wǎng)EWH(Equivalent Water Height)值。

    1.2 實驗方法

    采用球諧系數(shù)法反演TWSC的表達式如下所示? [14] :

    Δ h θ,λ = aρ?? ave?? 3ρ? w???SymboleB@ ?l=0??? l ?m=0? P - ??lm?? cos θ ×? 2l+1 1+k l? ?Δ C? lm? cos mλ+ Δ S? lm? sin mλ?? (1)

    式中:Δ h 為TWSC; θ 為地心余緯; λ 為地心經(jīng)度; a 為地球平均半徑,取6 371 km; ρ?? ave 為地球平均密度,即 5 517 ?kg/m 3; ρ? w為水密度,即1 000 kg/m 3; l 為階數(shù); m 為次數(shù); P? - ???lm?? cos θ? 為規(guī)格化締合勒讓德函數(shù); k l 為勒夫數(shù);Δ C??? lm? 與Δ S??? lm? 為重力場球諧系數(shù)變化值,本文中是指每個月的球諧系數(shù)扣除平均重力場(2018年6月至2020年3月的球諧系數(shù)平均值)。

    在反演前,對球諧系數(shù)一階項進行地球質(zhì)心改正? [15] ,并將原始的 C?? 20 項替換為SLR分析得到的結(jié)果? [16] 。同時為了削弱由重力場模型本身誤差和衛(wèi)星軌道及荷載誤差所引起的條帶誤差,本文采用250 km Fan濾波? [17] 加P3M6多項式濾波? [18-19] 的組合算法對球諧系數(shù)進行處理。

    為了評估重力場模型的精度,常常使用由球諧系數(shù)誤差計算得到的階方差作為評價指標(biāo),其表達式如下? [10] :

    δσ l=? 1 2l+1 ??l ?m=0? ?σC 2? lm +σS 2? lm? ??(2)

    式中:σC? lm 和σS? lm 為球諧系數(shù)誤差。

    2 實驗分析

    2.1 模型精度比較

    本文采用2018年6月至2020年3月上述3個機構(gòu)的GRACE-FO模型數(shù)據(jù),計算相應(yīng)時變重力場模型的平均階方差,結(jié)果如圖1所示。由圖1可知,3種不同的GRACE-FO模型平均階方差變化趨勢基本相同,可以分為較為明顯的3個階段:在10階之前表現(xiàn)為減少趨勢,10階到20階時則為緩慢增加趨勢,超過20階以后階方差值的增速在變快。同時JPL RL06的平均階方差值要明顯大于其他2種數(shù)據(jù),CSR RL06的結(jié)果與GFZ RL06較為接近。

    由于衛(wèi)星軌道傾角的原因,GRACE-FO對于球諧系數(shù)中的 C?? 20 項不夠敏感,這導(dǎo)致該系數(shù)存在較大的誤差,因此在GRACE-FO時變重力場數(shù)據(jù)處理過程中,一般采用SLR提供的數(shù)據(jù)進行替換。為了進一步比較3種模型之間的異同,本文將3種模型的 C?? 20 項時間序列與SLR的觀測結(jié)果進行比較(見圖2(a))。CSR和JPL的 C?? 20 項系數(shù)與SLR的結(jié)果基本一致,但是GFZ的 C?? 20 項系數(shù)與其他三者的結(jié)果有著較大的差別,這可能是因為在處理數(shù)據(jù)過程中所采用的靜態(tài)背景場不同所造成的,因此相對于其他2個機構(gòu)的結(jié)果,GFZ存在著系統(tǒng)性的偏差? [9,11] 。

    由于利用球諧系數(shù)法反演TWSC時,采用的是球諧系數(shù)變化量,因此,對 C? ?20 系數(shù)殘差項(即每個月 C?? 20 系數(shù)減去所有月份 C?? 20 系數(shù)平均值)進行比較更具科學(xué)意義(見圖2(b))。從圖2可知:4種不同來源的 C?? 20 殘差時間序列差別不是很大,其中CSR和JPL結(jié)果更接近于SLR,而GFZ結(jié)果要大于其他2個機構(gòu),這與圖2(a)的分析相一致。

    2.2 全球TWSC反演結(jié)果

    本文以GLDAS水文模型結(jié)果作為驗證值,分別與三大機構(gòu)的反演結(jié)果進行比較(見圖3)。由于GLDAS模型在南極大陸和格陵蘭島區(qū)域沒有數(shù)據(jù),因此GRACE-FO結(jié)果扣除了這2個地區(qū)的數(shù)據(jù)。本文中所提到的陸地水儲量包括了地表水(土壤水、植被水、積雪水和徑流等)和地下水兩部分。由圖3可知,4種數(shù)據(jù)的TWSC全球分布基本相同,其中亞馬遜和密西西比河流域TWSC呈現(xiàn)出減少的趨勢,而在剛果河和恒河流域則表現(xiàn)為增加趨勢。但是在局部地區(qū)也存在差別,例如剛果河和恒河流域的增加幅度略有不同,亞馬遜流域的減少幅度也不一樣。這可能是由GRACE-FO和GLDAS本身模型誤差所造成的,并且GRACE-FO的TWSC結(jié)果不僅包括土壤水、積雪水和植被水,還有地下水和人類活動引起的TWSC變化等。

    為了進一步比較分析,本文計算GRACE-FO和GLDAS全球陸地區(qū)域1°×1°TWSC格網(wǎng)數(shù)據(jù)時間序列的相關(guān)系數(shù)(見圖4)。比較結(jié)果表明,3種模型與GLDAS模型的相關(guān)系數(shù)分布大致相同。圖中白色區(qū)域的相關(guān)系數(shù)小于或等于0,可以發(fā)現(xiàn)白色區(qū)域主要集中在沙漠干旱地區(qū)(例如撒哈拉沙漠、塔克拉瑪干沙漠、蒙古高原、阿拉伯半島、澳大利亞沙漠等)。這是因為這些地區(qū)的TWSC非常小,甚至幾乎沒有,因此GRACE-FO衛(wèi)星在這些區(qū)域探測不到相應(yīng)的水文信號。

    2.3 局部TWSC反演結(jié)果

    為了驗證3種GRACE-FO模型反演局部TWSC的結(jié)果,本文選取了世界上較有代表性的六大流域作為實驗區(qū)域,分別為亞馬遜河、密西西比河、長江、奧里諾科河、剛果河和恒河。圖5給出了3種GRACE-FO模型和GLDAS水文模型的TWSC結(jié)果。由圖5可知,6個流域GRACE-FO模型的TWSC反演結(jié)果趨勢與GLDAS模型較為一致,但是GRACE-FO結(jié)果的變化幅度要略大于GLDAS。這是由于GLDAS水文模型只包含土壤水、積雪水和植被水,沒有考慮地下水和人類活動影響。3種GRACE-FO模型反演結(jié)果差別不大。

    本文分別計算了六大流域的TWSC結(jié)果相關(guān)系數(shù),如表1所列。由表1可知,相關(guān)系數(shù)最小的是剛果河,最大的是恒河。所有流域的相關(guān)系數(shù)均大于0.6,這與圖5是相符的。在6個流域中有4個流域的CSR相關(guān)系數(shù)最高,另外2個流域GFZ最高,JPL的結(jié)果較為平衡。

    造成3種GRACE-FO模型反演結(jié)果存在差異的原因主要是:① 在解算過程中采用的模型不同,例如衛(wèi)星攝動模型CSR和GFZ采用的是D430,而JPL采用的D421;② 在計算過程中不可避免出現(xiàn)計算誤差,這種誤差為隨機誤差。從3種模型反演結(jié)果與GLDAS的比較來看,不論在全球范圍內(nèi)還是局部區(qū)域3種模型其實是比較接近的,但是CSR提供的模型精度較其他2個機構(gòu)的要略高一些,因此建議在使用時采用CSR提供的模型。

    3 結(jié) 論

    本文采用CSR、GFZ和JPL三大機構(gòu)發(fā)布的GRACE-FO RL06時變重力場模型數(shù)據(jù)反演了全球和局部地區(qū)TWSC,并比較了3種模型的階方差、C? 20 項精度和TWSC反演結(jié)果。結(jié)果表明:CSR和JPL的GRACE-FO模型階方差和C? 20 項精度均優(yōu)于GFZ。3個模型的TWSC反演結(jié)果與GLDAS結(jié)果基本相同。通過對全球6個流域的TWSC結(jié)果和時間序列的相關(guān)系數(shù)分析可知,GRACE-FO和GLDAS模型的結(jié)果都具有較高的相關(guān)性,在多數(shù)流域中CSR的相關(guān)度較高。

    參考文獻:

    [1]? LONG D,LONGUEVERGNE L,SCANLON B R.Global analysis of approaches for deriving total water storage changes from GRACE satellite[J].Water Resources Research,2015,51(4):2574-2594.

    [2] 崔立魯,宋哲,鄒正波,等.基于GRACE的2009~2010年西南地區(qū)特大干旱分析[J].人民長江,2020,51(7):13-19.

    [3] 崔立魯,杜安,張誠,等.基于重力衛(wèi)星數(shù)據(jù)的華北平原地下水儲量變化分析[J].成都大學(xué)學(xué)報(自然科學(xué)版),2020,39(2):215-219.

    [4] WAHR J,MOLENAAR M,BRYAN F.Time variability of the Earth’s gravity field:hydrological and oceanic effects and their possible detection using GRACE[J].Journal Geophysics Research,1998,103(B12):30205-30229.

    [5] CHEN J L,WILSON C R,BLANKENSHIP D D,et al.Antarctic mass rates from GRACE[J].Geophysics Research Letters,2006,33(11):L11502.

    [6] 崔立魯,張誠,鄒正波,等.2003~2016年南極地區(qū)冰蓋質(zhì)量變化時間序列分析[J].大地測量與地球動力學(xué),2020,40(12):1219-1222,1236.

    [7] 崔立魯,張誠,鄒正波,等.利用重力衛(wèi)星監(jiān)測2003~2016年南極地區(qū)冰蓋質(zhì)量變化[J].科學(xué)技術(shù)與工程,2021,21(1):27-31.

    [8] BAUR O,KUHN M,F(xiàn)EATHERSTONE W E.Continental mass change from GRACE over 2002-2011 and its impact on sea level[J].Journal of Geodesy,2013,87(2):117-125.

    [9] HAN S C,RIVA R,SAUBER J,et al.Source parameter inversion for great earthquakes from a decade-long observation of global gravity fields[J].Journal of Geophysics Research-Solid Earth,2013,118(3):1240-1267.

    [10]? 崔立魯,周甜,張誠,等.利用重力衛(wèi)星監(jiān)測尼泊爾Ms8.1 地震前后重力變化[J].科學(xué)技術(shù)與工程,2020,20(28):11434-11441.

    [11] 冉艷紅,鐘敏,陳威,等.利用GRACE-FO重力衛(wèi)星探測2019年長江中下游極端干旱[J].科學(xué)通報,2021(66):107-117.

    [12] 李瓊.地表物質(zhì)遷移的時變重力場反演方法及其應(yīng)用研究[D].武漢:武漢大學(xué),2014.

    [13] 郭飛霄,孫中苗,任飛虎,等.GRACE RL06與RL05 時變重力場模型數(shù)據(jù)初步比較分析[J].大地測量與地球動力學(xué),2020,40(5):546-551.

    [14] CHENG M,TAPLEY B D.Variations in the Earth’s oblateness during the past 28 years[J].Journal of Geophysics Research-Solid Earth,2004,109(B9):2745-2750.

    [15] SWENSON S,CHAMBERS D,WHAR J.Estimating geocenter variations from a combination of GRACE and ocean model output[J].Journal of Geophysical Research-Solid Earth,2008,113(B8):194-205.

    [16] 崔立魯,宋哲,鄒正波,等.利用重力衛(wèi)星監(jiān)測中國西南地區(qū)陸地水儲量變化[J].科學(xué)技術(shù)與工程,2020,20(30):12313-12317.

    [17] 崔立魯,宋哲,鄒正波,等.重力衛(wèi)星時變重力場位系數(shù)誤差Fan濾波算法分析[J].科學(xué)技術(shù)與工程,2019,19(15):46-51.

    [18] 崔立魯,唐興友,鄒正波,等.重力衛(wèi)星時變重力場位系數(shù)誤差的反向延拓去相關(guān)算法[J].科學(xué)技術(shù)與工程,2019,19(26):23-29.

    [19] 崔立魯,唐興友,鄒正波,等.重力衛(wèi)星時變重力場去相關(guān)濾波算法參數(shù)選擇[J].科學(xué)技術(shù)與工程,2020,20(1):43-47.

    (編輯:謝玲嫻)

    Comparative analysis of terrestrial water storage change inverted by different GRACE-FO models

    CUI Lilu? 1,2 ,DU An 1,ZHANG Cheng 1,WANG Xiaolong 3,LI Qiong 4

    ( 1.School of Architecture and Civil Engineering,Chengdu University,Chengdu 610106,China; 2.School of Geodesy and Geomatics,Wuhan University,Wuhan 430071,China; 3.Nanning Exploration & Survey Geoinformation Institute,Nanning 530022,China; 4.School of Civil Engineering and Geomatics,Southwest Petroleum University,Chengdu 610500,China )

    Abstract:

    In order to effectively compare the accuracy of three GRACE-FO RL06 time-variable gravity field models provided by Center of Space Research (CSR),Helmholtz-Centre Potsdam-German Research Centre for Geosciences (GFZ) and Jet Propulsion Laboratory (JPL),we assessed these three models from the perspectives of model degree variance,C20 accuracy and inversion of global and local terrestrial water storage change.The results showed that the degree variance of JPL RL06 model was greater than those of other two institutions.In term of C20 accuracy,data released by CSR and JPL were closer to Satellite Laser Ranging (SLR).By comparing the results of terrestrial water storage changes inverted by three models and the results of Global Land Data Assimilation system (GLDAS) model,it was found that the results of the three models and GLDAS model were relatively in good agreement.From the local inversion results,the CSR and GLDAS models were in high agreement among the three models.

    Key words:

    terrestrial water storage change;model degree variance;C20 accuracy;correlation coefficient;GRACE-FO model

    猜你喜歡
    重力場時變方差
    方差怎么算
    概率與統(tǒng)計(2)——離散型隨機變量的期望與方差
    計算方差用哪個公式
    基于空間分布的重力場持續(xù)適配能力評估方法
    方差生活秀
    衛(wèi)星測量重力場能力仿真分析
    基于時變Copula的股票市場相關(guān)性分析
    智富時代(2017年4期)2017-04-27 17:08:47
    煙氣輪機復(fù)合故障時變退化特征提取
    基于MEP法的在役橋梁時變可靠度研究
    擾動重力場元無θ奇異性計算公式的推導(dǎo)
    在线免费观看不下载黄p国产| 99热只有精品国产| 国产真实伦视频高清在线观看| 插逼视频在线观看| 三级经典国产精品| 国产精品1区2区在线观看.| 亚洲性久久影院| 国产在线男女| 亚洲在线观看片| 22中文网久久字幕| 99热这里只有是精品在线观看| 一区二区三区免费毛片| av卡一久久| 国模一区二区三区四区视频| 69人妻影院| 青青草视频在线视频观看| 91麻豆精品激情在线观看国产| 亚洲乱码一区二区免费版| 色哟哟哟哟哟哟| 欧美三级亚洲精品| 中文亚洲av片在线观看爽| 在线观看免费视频日本深夜| 人人妻人人澡人人爽人人夜夜 | 老熟妇乱子伦视频在线观看| 成人永久免费在线观看视频| 91精品一卡2卡3卡4卡| 18+在线观看网站| 全区人妻精品视频| 日日撸夜夜添| 嫩草影院新地址| 欧美xxxx性猛交bbbb| 亚洲天堂国产精品一区在线| 一个人看的www免费观看视频| 自拍偷自拍亚洲精品老妇| 波多野结衣高清作品| 99国产精品一区二区蜜桃av| 人妻久久中文字幕网| 人人妻人人看人人澡| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久精品电影| 欧美日韩乱码在线| 成人二区视频| 亚洲最大成人中文| 18禁黄网站禁片免费观看直播| 成人特级av手机在线观看| 国产精品嫩草影院av在线观看| 中文字幕免费在线视频6| 人人妻人人看人人澡| 白带黄色成豆腐渣| 国产探花极品一区二区| 久久这里有精品视频免费| 男女那种视频在线观看| 欧美极品一区二区三区四区| 国内精品宾馆在线| 一级毛片我不卡| 男人和女人高潮做爰伦理| 久99久视频精品免费| 国产中年淑女户外野战色| 美女 人体艺术 gogo| 中文字幕精品亚洲无线码一区| 国产日本99.免费观看| 久久久国产成人精品二区| 中文字幕免费在线视频6| 一本精品99久久精品77| 亚洲三级黄色毛片| 久久久久免费精品人妻一区二区| 婷婷精品国产亚洲av| 日本黄大片高清| 国产老妇女一区| 亚洲乱码一区二区免费版| 精品一区二区三区人妻视频| 2021天堂中文幕一二区在线观| 国产毛片a区久久久久| 日韩一区二区视频免费看| kizo精华| 国产精品一区二区在线观看99 | 啦啦啦韩国在线观看视频| 可以在线观看毛片的网站| 国产激情偷乱视频一区二区| 国产在视频线在精品| 老师上课跳d突然被开到最大视频| 国产精品麻豆人妻色哟哟久久 | a级毛片a级免费在线| 99九九线精品视频在线观看视频| а√天堂www在线а√下载| videossex国产| 国产美女午夜福利| 亚洲精品456在线播放app| av黄色大香蕉| 天堂影院成人在线观看| 三级经典国产精品| 亚州av有码| 日韩高清综合在线| 国产在视频线在精品| av黄色大香蕉| 校园春色视频在线观看| 欧美极品一区二区三区四区| 久久精品国产亚洲av天美| 99热这里只有是精品在线观看| 在线a可以看的网站| 国产成人一区二区在线| 亚洲国产精品合色在线| 看十八女毛片水多多多| 人妻夜夜爽99麻豆av| 欧美三级亚洲精品| 美女脱内裤让男人舔精品视频 | 亚洲av中文字字幕乱码综合| 天堂网av新在线| 在线播放无遮挡| 久久这里有精品视频免费| 看片在线看免费视频| 可以在线观看毛片的网站| 尾随美女入室| 美女脱内裤让男人舔精品视频 | 99久久九九国产精品国产免费| 欧美最黄视频在线播放免费| 国产单亲对白刺激| 国产精华一区二区三区| 日韩成人av中文字幕在线观看| 亚洲18禁久久av| 久久精品久久久久久噜噜老黄 | 少妇的逼水好多| 国产成人一区二区在线| 欧美精品一区二区大全| 少妇丰满av| 亚洲av免费在线观看| 狂野欧美激情性xxxx在线观看| 亚洲七黄色美女视频| 亚洲精品国产av成人精品| 美女 人体艺术 gogo| 中文字幕免费在线视频6| 天堂网av新在线| 12—13女人毛片做爰片一| 床上黄色一级片| 国产人妻一区二区三区在| 精品久久久久久久人妻蜜臀av| 欧美又色又爽又黄视频| 色播亚洲综合网| 天堂影院成人在线观看| 亚洲av中文av极速乱| 一级毛片久久久久久久久女| 日本在线视频免费播放| 日韩欧美 国产精品| 亚洲精品国产成人久久av| 国产成人aa在线观看| 人人妻人人看人人澡| 欧美日韩在线观看h| а√天堂www在线а√下载| 内地一区二区视频在线| 人妻夜夜爽99麻豆av| 成人特级av手机在线观看| 成年免费大片在线观看| 国语自产精品视频在线第100页| 久久久成人免费电影| 久久久久九九精品影院| 午夜免费男女啪啪视频观看| 男人舔女人下体高潮全视频| 国产精品久久久久久精品电影小说 | 特级一级黄色大片| 26uuu在线亚洲综合色| 最近视频中文字幕2019在线8| 真实男女啪啪啪动态图| 内地一区二区视频在线| 老师上课跳d突然被开到最大视频| 国产成人精品一,二区 | 亚洲经典国产精华液单| 亚洲精品自拍成人| 国产高清激情床上av| 久久久久久大精品| 欧美xxxx黑人xx丫x性爽| 神马国产精品三级电影在线观看| av福利片在线观看| 亚洲真实伦在线观看| 亚洲精品色激情综合| 成人无遮挡网站| 啦啦啦啦在线视频资源| 99久久成人亚洲精品观看| 国产探花极品一区二区| 午夜a级毛片| 国产成人91sexporn| а√天堂www在线а√下载| 国产v大片淫在线免费观看| 久久精品国产亚洲av涩爱 | 一个人免费在线观看电影| 蜜桃亚洲精品一区二区三区| 黄色配什么色好看| 97人妻精品一区二区三区麻豆| 欧美性感艳星| 一级黄片播放器| 欧美激情国产日韩精品一区| 久久99蜜桃精品久久| 变态另类丝袜制服| 偷拍熟女少妇极品色| 18禁裸乳无遮挡免费网站照片| 99久国产av精品| 麻豆成人av视频| 久久欧美精品欧美久久欧美| 2022亚洲国产成人精品| 国产高清激情床上av| 亚洲精品日韩av片在线观看| 国产精品免费一区二区三区在线| 女人十人毛片免费观看3o分钟| 成熟少妇高潮喷水视频| 少妇的逼水好多| 亚洲精品亚洲一区二区| 激情 狠狠 欧美| 最近的中文字幕免费完整| 国产精品电影一区二区三区| 日韩一本色道免费dvd| 精品久久久久久久久av| 日本欧美国产在线视频| 高清在线视频一区二区三区 | 久久久久久久久久久丰满| 亚洲成av人片在线播放无| 日本色播在线视频| 亚洲av男天堂| 欧美zozozo另类| 日本黄色片子视频| 国产精品人妻久久久久久| 男的添女的下面高潮视频| 久久精品影院6| 免费看美女性在线毛片视频| 亚洲欧美中文字幕日韩二区| 国产在视频线在精品| 九九在线视频观看精品| 国产一区二区在线观看日韩| 两个人视频免费观看高清| 国产美女午夜福利| 内地一区二区视频在线| 国产亚洲精品av在线| 日本五十路高清| 久久精品综合一区二区三区| videossex国产| 久久精品久久久久久久性| 免费人成在线观看视频色| 国产亚洲精品久久久久久毛片| 亚洲综合色惰| 69av精品久久久久久| 啦啦啦啦在线视频资源| 18+在线观看网站| 黄片无遮挡物在线观看| 亚洲国产欧美在线一区| 成人无遮挡网站| 最近最新中文字幕大全电影3| 熟女电影av网| 国内揄拍国产精品人妻在线| 欧美在线一区亚洲| 91精品国产九色| 免费av观看视频| av在线蜜桃| 99在线人妻在线中文字幕| 男插女下体视频免费在线播放| 国产精品久久电影中文字幕| 国产综合懂色| 婷婷亚洲欧美| av国产免费在线观看| 亚洲国产精品成人综合色| 国产精品蜜桃在线观看 | 在线天堂最新版资源| 美女大奶头视频| 一级二级三级毛片免费看| 亚洲成人久久爱视频| 99热这里只有是精品50| 欧美+亚洲+日韩+国产| 亚洲无线在线观看| 日韩一区二区三区影片| 欧美高清成人免费视频www| 中文字幕av成人在线电影| 亚洲最大成人av| 在线免费十八禁| 中文字幕人妻熟人妻熟丝袜美| 搞女人的毛片| 成人三级黄色视频| 国产成人福利小说| 男插女下体视频免费在线播放| 国产91av在线免费观看| 国产精品人妻久久久久久| 国产成人一区二区在线| 黄色配什么色好看| 国产亚洲精品久久久久久毛片| 日本黄大片高清| 亚洲真实伦在线观看| 亚洲美女视频黄频| 日韩av在线大香蕉| 国产 一区 欧美 日韩| 美女脱内裤让男人舔精品视频 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品.久久久| 亚洲中文字幕一区二区三区有码在线看| 日本三级黄在线观看| 男女那种视频在线观看| 国产白丝娇喘喷水9色精品| www日本黄色视频网| 偷拍熟女少妇极品色| 久久人人爽人人片av| 国产av麻豆久久久久久久| 日本免费一区二区三区高清不卡| 亚洲精品国产成人久久av| 亚洲成a人片在线一区二区| 亚洲国产日韩欧美精品在线观看| 久久99蜜桃精品久久| 69av精品久久久久久| 一级av片app| 大香蕉久久网| 亚洲美女视频黄频| 日韩欧美精品免费久久| 国产亚洲av片在线观看秒播厂 | 久久婷婷人人爽人人干人人爱| 欧美zozozo另类| 久久精品国产亚洲av天美| 三级男女做爰猛烈吃奶摸视频| 国产视频首页在线观看| 国产精品野战在线观看| 国产精品,欧美在线| 欧美成人精品欧美一级黄| 亚洲精品亚洲一区二区| 成人性生交大片免费视频hd| 寂寞人妻少妇视频99o| 国产精品久久久久久av不卡| 久久久久国产网址| 青青草视频在线视频观看| 99久久精品一区二区三区| 97超视频在线观看视频| 乱码一卡2卡4卡精品| 26uuu在线亚洲综合色| 国产精品一二三区在线看| 男的添女的下面高潮视频| 久久久久九九精品影院| 美女脱内裤让男人舔精品视频 | 精品久久久久久久人妻蜜臀av| 爱豆传媒免费全集在线观看| 亚洲久久久久久中文字幕| 久久久a久久爽久久v久久| 人人妻人人澡欧美一区二区| 久久久久性生活片| 亚洲精品乱码久久久久久按摩| 一级毛片我不卡| 国产精品一区二区三区四区免费观看| 69人妻影院| 中国美白少妇内射xxxbb| 麻豆国产97在线/欧美| 国产又黄又爽又无遮挡在线| www.色视频.com| 一卡2卡三卡四卡精品乱码亚洲| 久久99蜜桃精品久久| 狂野欧美白嫩少妇大欣赏| av在线蜜桃| 免费观看的影片在线观看| 国产高清不卡午夜福利| 一边摸一边抽搐一进一小说| 欧美日韩综合久久久久久| 国产黄色小视频在线观看| 日韩欧美在线乱码| 丰满乱子伦码专区| 国产av麻豆久久久久久久| 亚洲国产欧美在线一区| 日韩大尺度精品在线看网址| 久久人人爽人人爽人人片va| 成人美女网站在线观看视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲欧洲国产日韩| 又爽又黄无遮挡网站| 久久久午夜欧美精品| 国产真实乱freesex| 精品熟女少妇av免费看| 在线观看午夜福利视频| 成人亚洲欧美一区二区av| av卡一久久| 欧美不卡视频在线免费观看| av在线播放精品| 欧美日韩国产亚洲二区| 观看美女的网站| 成人午夜精彩视频在线观看| 男女做爰动态图高潮gif福利片| 美女xxoo啪啪120秒动态图| 国内揄拍国产精品人妻在线| 美女cb高潮喷水在线观看| 久久这里有精品视频免费| 村上凉子中文字幕在线| 99久国产av精品国产电影| 国产成人freesex在线| 成年女人永久免费观看视频| 一卡2卡三卡四卡精品乱码亚洲| 91久久精品电影网| 欧美性感艳星| 老女人水多毛片| 欧美bdsm另类| 国产精品嫩草影院av在线观看| 国产精品不卡视频一区二区| 日韩制服骚丝袜av| 好男人在线观看高清免费视频| 麻豆一二三区av精品| 亚洲精品久久久久久婷婷小说 | 夜夜看夜夜爽夜夜摸| 偷拍熟女少妇极品色| 午夜精品一区二区三区免费看| 91aial.com中文字幕在线观看| 日本一二三区视频观看| 亚洲国产日韩欧美精品在线观看| 欧美日韩精品成人综合77777| 免费看美女性在线毛片视频| 婷婷色av中文字幕| 亚洲一区二区三区色噜噜| av卡一久久| 国产亚洲精品av在线| 亚洲av.av天堂| 国产精品一二三区在线看| 国产成人精品久久久久久| 嫩草影院精品99| 日本撒尿小便嘘嘘汇集6| 夫妻性生交免费视频一级片| 亚洲欧美成人综合另类久久久 | 在线免费观看不下载黄p国产| 国产黄片美女视频| 国产精品久久久久久精品电影| 最近2019中文字幕mv第一页| 少妇的逼水好多| 在线免费观看的www视频| 国产69精品久久久久777片| 亚洲第一区二区三区不卡| 我要搜黄色片| 国产精品电影一区二区三区| 日产精品乱码卡一卡2卡三| 能在线免费看毛片的网站| 伦精品一区二区三区| 最近视频中文字幕2019在线8| 啦啦啦韩国在线观看视频| 国产三级中文精品| 天天躁日日操中文字幕| 亚洲在线自拍视频| 婷婷六月久久综合丁香| 国产精品久久久久久久久免| 国产免费一级a男人的天堂| 丝袜美腿在线中文| 啦啦啦观看免费观看视频高清| 边亲边吃奶的免费视频| 天天躁日日操中文字幕| 国产成人午夜福利电影在线观看| 97人妻精品一区二区三区麻豆| 老司机影院成人| 免费av不卡在线播放| 搡老妇女老女人老熟妇| 成人一区二区视频在线观看| av国产免费在线观看| 狠狠狠狠99中文字幕| 人妻少妇偷人精品九色| 三级经典国产精品| 嫩草影院精品99| av在线老鸭窝| 国产真实伦视频高清在线观看| 久久久国产成人免费| 免费看光身美女| 亚洲最大成人av| 国产探花在线观看一区二区| 日韩精品有码人妻一区| 夜夜夜夜夜久久久久| 成人毛片60女人毛片免费| 国模一区二区三区四区视频| 不卡视频在线观看欧美| 国产视频首页在线观看| 国产激情偷乱视频一区二区| 亚州av有码| 内地一区二区视频在线| 亚洲欧美日韩无卡精品| 人妻少妇偷人精品九色| 亚洲欧美精品自产自拍| 国产黄片视频在线免费观看| 亚洲av电影不卡..在线观看| 日韩欧美 国产精品| videossex国产| 久久久久久久久久成人| 波多野结衣高清作品| 中文字幕人妻熟人妻熟丝袜美| 久久国产乱子免费精品| 人妻少妇偷人精品九色| 精品一区二区三区人妻视频| 久久久久九九精品影院| 色综合色国产| 午夜亚洲福利在线播放| 欧洲精品卡2卡3卡4卡5卡区| 12—13女人毛片做爰片一| 男人的好看免费观看在线视频| 亚洲性久久影院| 精品午夜福利在线看| 老女人水多毛片| 免费av毛片视频| 亚洲国产精品久久男人天堂| 一个人看视频在线观看www免费| 亚洲内射少妇av| 秋霞在线观看毛片| 人人妻人人澡人人爽人人夜夜 | 久久精品夜夜夜夜夜久久蜜豆| 最近2019中文字幕mv第一页| 熟女电影av网| 免费搜索国产男女视频| 伊人久久精品亚洲午夜| 欧美性猛交╳xxx乱大交人| 欧美变态另类bdsm刘玥| 国产国拍精品亚洲av在线观看| 国产成人精品一,二区 | 亚洲第一电影网av| 少妇被粗大猛烈的视频| 国产黄色视频一区二区在线观看 | 男女那种视频在线观看| 亚洲自拍偷在线| 三级毛片av免费| 亚洲中文字幕一区二区三区有码在线看| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区在线观看99 | 日本欧美国产在线视频| 亚洲经典国产精华液单| 一级毛片久久久久久久久女| 亚洲欧美成人精品一区二区| 一级毛片电影观看 | 六月丁香七月| 国产精品电影一区二区三区| 日韩欧美三级三区| 国产亚洲av片在线观看秒播厂 | 亚洲av免费高清在线观看| 婷婷精品国产亚洲av| 一级毛片久久久久久久久女| 国产午夜精品论理片| 观看美女的网站| 久久久久久久久久成人| 国产视频首页在线观看| 国产色婷婷99| eeuss影院久久| 免费在线观看成人毛片| 国语自产精品视频在线第100页| 国产成年人精品一区二区| 如何舔出高潮| 自拍偷自拍亚洲精品老妇| 精品欧美国产一区二区三| 三级男女做爰猛烈吃奶摸视频| videossex国产| 天天一区二区日本电影三级| 成人漫画全彩无遮挡| 免费人成在线观看视频色| 99久久精品热视频| 国产v大片淫在线免费观看| 寂寞人妻少妇视频99o| 性色avwww在线观看| 亚洲最大成人中文| 最近视频中文字幕2019在线8| 美女被艹到高潮喷水动态| av免费在线看不卡| 极品教师在线视频| 欧美日韩一区二区视频在线观看视频在线 | 国产淫片久久久久久久久| 午夜福利在线观看免费完整高清在 | 99国产精品一区二区蜜桃av| 中国美白少妇内射xxxbb| 亚洲精品国产成人久久av| 国产精品野战在线观看| 亚洲av成人av| 国产在视频线在精品| 欧美又色又爽又黄视频| 三级毛片av免费| 1000部很黄的大片| 91久久精品国产一区二区成人| 欧美+日韩+精品| 国产av不卡久久| 久久久久久久亚洲中文字幕| 97超视频在线观看视频| 久久久a久久爽久久v久久| 久久中文看片网| 又黄又爽又刺激的免费视频.| 国产伦精品一区二区三区视频9| 又爽又黄a免费视频| 国产精品人妻久久久影院| 联通29元200g的流量卡| 亚洲av免费在线观看| 免费人成在线观看视频色| 国产一区二区激情短视频| 大又大粗又爽又黄少妇毛片口| 亚洲国产欧美在线一区| 精品久久久久久久久久免费视频| 高清日韩中文字幕在线| 久久精品国产鲁丝片午夜精品| 久久这里只有精品中国| 高清毛片免费看| 亚洲美女视频黄频| 日韩高清综合在线| 婷婷亚洲欧美| 99久久无色码亚洲精品果冻| 黄色配什么色好看| 69av精品久久久久久| 变态另类成人亚洲欧美熟女| 亚洲精品久久国产高清桃花| 婷婷精品国产亚洲av| 午夜老司机福利剧场| 色哟哟哟哟哟哟| 久99久视频精品免费| 日本五十路高清| 国产视频内射| 亚洲中文字幕一区二区三区有码在线看| 亚洲七黄色美女视频| 成年版毛片免费区| 亚洲高清免费不卡视频| 免费av不卡在线播放| 在线观看午夜福利视频| 尾随美女入室| 国产成人精品婷婷| 国产精品野战在线观看| 国产高清视频在线观看网站| 极品教师在线视频| 亚洲精品日韩在线中文字幕 | 变态另类成人亚洲欧美熟女| 99riav亚洲国产免费| 欧美三级亚洲精品| 老司机影院成人| 亚洲av中文av极速乱| a级毛色黄片| 久久人人精品亚洲av| 美女高潮的动态|