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

    GPM衛(wèi)星降水?dāng)?shù)據(jù)產(chǎn)品適用性評(píng)價(jià)
    ——以雅魯藏布江為例

    2021-07-03 06:24:24劉艷麗吳永祥王高旭金君良賀瑞敏劉翠善
    中國農(nóng)村水利水電 2021年6期
    關(guān)鍵詞:雅江尺度降水

    姜 曦,劉艷麗,吳永祥,王高旭,吳 巍,金君良,賀瑞敏,劉翠善

    (1.南京水利科學(xué)研究院水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,南京210029;2.河海大學(xué)水文水資源學(xué)院,南京210098;3.長江保護(hù)與綠色發(fā)展研究院,南京210098;4.水利部應(yīng)對氣候變化研究中心,南京210029)

    降水是水文循環(huán)的重要環(huán)節(jié)[1],是水文分析的重要指標(biāo),同時(shí)也是水文、氣象、環(huán)境等學(xué)科主要研究對象[2]。降水資料的準(zhǔn)確性、完整性和客觀性對工程實(shí)踐與科學(xué)研究具有至關(guān)重要的作用[3]。我國主要的降水信息來源于水文站、氣象站等地面雨量計(jì),其能夠提供長系列、高精度的降水信息,但對于地面觀測站點(diǎn)稀疏地區(qū),地面觀測站的降水資料代表性較差,不能準(zhǔn)確的反映該區(qū)域降水的時(shí)空分布特征,這對水文分析計(jì)算與水資源優(yōu)化配置產(chǎn)生較大影響[4-6]。隨著遙感衛(wèi)星應(yīng)用的發(fā)展,多源、多時(shí)相、多尺度以及不同類型遙感衛(wèi)星的應(yīng)用,衛(wèi)星降水產(chǎn)品的開發(fā)利用發(fā)展為降水研究的重要趨勢[7],其覆蓋范圍與時(shí)空連續(xù)性均能作為地面站點(diǎn)的有效補(bǔ)充。這為地面站點(diǎn)分布稀疏、降水資料缺失區(qū)域的降水研究提供了有力的數(shù)據(jù)支撐[8]。雅魯藏布江(雅江)流域地處高寒區(qū)域,地面雨量計(jì)分布稀疏、不均勻,地面資料代表性不足,是我國典型的缺資料地區(qū),研究衛(wèi)星降水產(chǎn)品的適用性具有重要現(xiàn)實(shí)意義[9]。

    全球范圍內(nèi)應(yīng)用較為廣泛的降水?dāng)?shù)據(jù)產(chǎn)品有TRMM(Tropical Rainfall Measuring Mission)、GPM-IMERG(Integrated Multisatellite Retrievals for GPM)、GSMap(Global Satellite Mapping of Precipitation)、PERSIANN(Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks)、APHRODITE(Asian Precipitation-Highly-Resolved Observational Data Integration Towards Evaluation)、CMORPH(Climate Prediction Center Morphing Technique)、CRU(Climatic Research Unit)、GLDAS(Global Land Data Assimilation System)等(表1)[1-6],這些數(shù)據(jù)產(chǎn)品已經(jīng)應(yīng)用在水文、氣象、生態(tài)等研究領(lǐng)域[10]。TRMM 降水?dāng)?shù)據(jù)已被驗(yàn)證在亞洲、美洲等多地具有較好的精度和適用性[4,11,12],GPM 衛(wèi)星降水產(chǎn)品為NASA 與JAXA 聯(lián)合發(fā)布的TRMM 降水產(chǎn)品的新一代產(chǎn)品,GPM-IMERG 與GSMap 為GPM(Global Precipitation Measurement)系列數(shù)據(jù)產(chǎn)品中的兩個(gè)代表,其空間分辨率由TRMM 衛(wèi)星降水產(chǎn)品的0.25°提升至0.1°,GSMap 時(shí)間分辨率由TRMM 衛(wèi)星降水產(chǎn)品的3 h 提升至1 h,GPM IMERG 提升至0.5 h[12]。其時(shí)空精度均高于TRMM 降水產(chǎn)品,更能滿足水文模擬、預(yù)報(bào)等研究對降水?dāng)?shù)據(jù)的需求,更加具有應(yīng)用潛力。

    表1 常見數(shù)據(jù)產(chǎn)品時(shí)空精度信息Tab.1 Spatio-temporal accuracy information of common data products

    衛(wèi)星數(shù)據(jù)產(chǎn)品的精度和準(zhǔn)確性對后續(xù)的研究有著至關(guān)重要的作用。目前,常用的降水?dāng)?shù)據(jù)產(chǎn)品質(zhì)量評(píng)估方法有定量分析與分類分析等,定量分析包括相關(guān)性分析,偏差分析,均方誤差分析、均方根誤差分析等;分類分析包括探照率分析,錯(cuò)報(bào)率分析等[10,13]。陳曉宏[14]等通過相關(guān)性、錯(cuò)報(bào)率等指標(biāo)驗(yàn)證GPM IMERG降水產(chǎn)品在北江流域具有良好的適用性;XU R[8]等按氣象站點(diǎn)高程進(jìn)行劃分,發(fā)現(xiàn)在青藏高原高海拔地區(qū),GPM IMERG 與TRMM 3B42 V7 對降水的檢測能力較低;M Yang[15]等按降雨強(qiáng)度對降雨事件進(jìn)行劃分,通過雨量站降水?dāng)?shù)據(jù)與IMERG 降水?dāng)?shù)據(jù)進(jìn)行相關(guān)性、均方根誤差等的分析,檢驗(yàn)IMERG 衛(wèi)星降水產(chǎn)品在四川省的精度。FRANCHITO S H[16]等通過不同季節(jié)和不同區(qū)域的劃分,在巴西境內(nèi)進(jìn)行TRMM 衛(wèi)星降水產(chǎn)品評(píng)估。

    在雅江流域,黃浠[17]等在年、月尺度上對中國地面降水網(wǎng)格數(shù)據(jù)、CRU 和GLDAS 的精度進(jìn)行評(píng)估,發(fā)現(xiàn)在月尺度和年尺度的網(wǎng)格數(shù)據(jù)與實(shí)測降水?dāng)?shù)值最接近,而CRU 和GLDAS 與實(shí)測降水量相差較大;孫赫[9]等對選取CMA和APHRODITE、PERSIANN 和GPM、GLDAS、HAR(High Asia Refined analysis)與實(shí)測數(shù)據(jù)對比,結(jié)果表明APHRODITE、GPM 和HAR 降水低估(10%~30%),PERSIANN 和GLDAS 整體高估上游流域站點(diǎn)降水(28%~60%),但低估下游流域站點(diǎn)降水(11%~21%);呂洋[18]等選取相關(guān)系數(shù)、相對偏差對TRMM 3B42日降水?dāng)?shù)據(jù)與TRMM 3B43月降水?dāng)?shù)據(jù)在雅江流域16 個(gè)地面站進(jìn)行站點(diǎn)精度評(píng)估,表明TRMM月降水產(chǎn)品質(zhì)量很高,日尺度產(chǎn)品相關(guān)性較差;劉江濤[3]等對TRMM 3B42 與PERSIANN 在降水產(chǎn)品精度進(jìn)行的定量指標(biāo)和分類指標(biāo)評(píng)估,結(jié)果顯示兩種降水?dāng)?shù)據(jù)對弱降水的估值偏高,而對于強(qiáng)降水的估值偏低。目前,關(guān)于雅江流域的衛(wèi)星降水?dāng)?shù)據(jù)評(píng)估已有一些成果,但主要以TRMM 衛(wèi)星產(chǎn)品為主,GPM 降水產(chǎn)品空間范圍覆蓋更廣,時(shí)空精度更高,更符合水文模擬與預(yù)報(bào)的數(shù)據(jù)需求,本文將通過定量分析與分類分析結(jié)合的方式,對GPM 系列降水產(chǎn)品(GPM IMERG、GSMap)進(jìn)行時(shí)空尺度的適用性分析。

    1 研究區(qū)及降水?dāng)?shù)據(jù)概況

    1.1 研究區(qū)域概況

    雅魯藏布江起源于我國青藏高原寒區(qū),是我國最長的高原河流,水資源量豐富,其流域內(nèi)冰川廣泛分布,流域面積24.04萬km2,河長2 057 km。與瀾滄江、怒江每年出境水量高達(dá)5 000 多億m3,接近我國年用水總量6 000 億m3左右。相對于中國東部地區(qū),雅江流域獨(dú)特的地理特征與氣候環(huán)境為水文要素的監(jiān)測增加困難,地面監(jiān)測站點(diǎn)稀少,氣象站僅有13 個(gè)(圖1),對流域時(shí)空特征的評(píng)估帶來較大困難,對流域進(jìn)行基本水文模擬造成一定程度的影響。衛(wèi)星降水產(chǎn)品在很大程度上可以彌補(bǔ)雅江流域地面站點(diǎn)稀少、數(shù)據(jù)代表性差的不足,為雅江流域的水文氣象等分析提供數(shù)據(jù)基礎(chǔ)。

    圖1 雅魯藏布江流域概況圖Fig.1 Overview of the Yarlung Zangbo River Basin

    1.2 降水資料選取

    1.2.1 實(shí)測數(shù)據(jù)的選取

    對于單一站點(diǎn)數(shù)據(jù)精度的分析,考慮降水資料的完整性與代表性,本文選取雅江流域內(nèi)部波密站、林芝站、江孜站、拉孜站、當(dāng)雄站、普蘭站共7個(gè)地面觀測站實(shí)測值作為評(píng)估依據(jù)。時(shí)間序列為2001.01.01-2019.12.31年的日降水序列與月降水序列,數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn)。

    由于衛(wèi)星降水?dāng)?shù)據(jù)產(chǎn)品多為柵格數(shù)據(jù),本文在衛(wèi)星數(shù)據(jù)產(chǎn)品時(shí)空精度分析過程中,選取中國地面降水格點(diǎn)數(shù)據(jù)集(V2.0)作為實(shí)測降水?dāng)?shù)據(jù)與衛(wèi)星降水產(chǎn)品進(jìn)行時(shí)空精度分析,空間分別率為0.5°×0.5°,時(shí)間序列為2001.01.01-2019.12.31年日降水序列與2014.03-2018.08月降水序列,數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn)。

    1.2.2 衛(wèi)星數(shù)據(jù)的選取

    衛(wèi)星降水?dāng)?shù)據(jù)選取為GPM 系列降水產(chǎn)品,GPM-IMERG(https://gpm.nasa.gov/data-access/downloads/gpm)與 GSMap(https://sharaku.eorc.jaxa.jp/GsMap/index.htm)的日尺度降水與月尺度降水,空間分辨率為0.1°×0.1°。GPM-IMERG 時(shí)間序列選擇為2014.03.01-2019.12.31 的GPM-IMERG_Final日尺度數(shù)據(jù)產(chǎn)品與月尺度數(shù)據(jù)產(chǎn)品,GSMap日尺度衛(wèi)星降水產(chǎn)品時(shí)間序列選取為2001.01.01-2019.12.31,其中2001.01-2010.12 為GSMap_MVK V5 數(shù)據(jù)集,2011.01-2014.02 為GSMaP_RNL 數(shù)據(jù)集,2014.03-2019.12 為GSMap_MVK V7 數(shù)據(jù)集。由于GSMap MVK月尺度衛(wèi)星降水產(chǎn)品起始時(shí)間為2014年3月,本文月尺度適用性分析的時(shí)間序列選取為2014.03-2018.08。

    2 數(shù)據(jù)產(chǎn)品質(zhì)量評(píng)估方法

    由于中國地面降水格點(diǎn)數(shù)據(jù)集(V2.0)空間精度為0.5°的格網(wǎng)數(shù)據(jù),GPM 系列降水?dāng)?shù)據(jù)集空間精度為0.1°,本文選取最鄰近邊界插值法對地面格網(wǎng)數(shù)據(jù)進(jìn)行空間插值,從而實(shí)現(xiàn)數(shù)據(jù)集空間精度的統(tǒng)一。

    質(zhì)量評(píng)估采用定量分析和分類分析結(jié)合的方式,定量分析指標(biāo)選擇Person 相關(guān)系數(shù)(CC)、偏差(Bias)和均方根誤差(Rmse)。分類分析選取指標(biāo)為探照率(POD)、錯(cuò)報(bào)率(FAR)。Person 相關(guān)系數(shù)(CC)用來描述衛(wèi)星觀測數(shù)據(jù)與地面實(shí)測數(shù)據(jù)的相關(guān)關(guān)系強(qiáng)弱,偏差(Bias)與均方根誤差(Rmse)用來描述衛(wèi)星觀測數(shù)據(jù)與地面實(shí)測數(shù)據(jù)的偏離程度,代表數(shù)據(jù)的準(zhǔn)確程度。探照率(POD)與錯(cuò)報(bào)率(FAR)用來描述衛(wèi)星監(jiān)測能力與監(jiān)測的準(zhǔn)確程度。各指標(biāo)具體計(jì)算方法如下:

    式中:n代表參與計(jì)算的降水總天數(shù);i代表第i個(gè)衛(wèi)星降水產(chǎn)品降水值或地面觀測站點(diǎn)實(shí)測值;A表示與對應(yīng)于地面觀測站的衛(wèi)星數(shù)據(jù)產(chǎn)品降水值,mm;表示衛(wèi)星數(shù)據(jù)產(chǎn)品降水值的平均值,mm;B表示地面雨量站測得降水值,mm;表示地面雨量站測得降水值的平均值,mm;T代表降水衛(wèi)星正確監(jiān)測降水事件的次數(shù);F代表實(shí)際未發(fā)生降雨事件,降水衛(wèi)星誤判為發(fā)生降水事件的次數(shù);M代表實(shí)際發(fā)生降水事件,而降水衛(wèi)星未能探測到降水的次數(shù)(表2)。

    表2 衛(wèi)星降水與地面觀測站降水觀測情況表Tab.2 Table of precipitation observation about satellite precipitation and ground observation station

    3 結(jié)果與討論

    3.1 流域整體精度評(píng)定

    3.1.1 時(shí)間尺度

    地面實(shí)測降水量與GPM 降水系列兩種數(shù)據(jù)產(chǎn)品月、日尺度數(shù)據(jù)集的定量分析結(jié)果見圖2、圖3,可以看出,對于月尺度降水產(chǎn)品,兩種衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測數(shù)據(jù)一致性表現(xiàn)較好,均達(dá)到顯著水平,GSMap 與GPM-IMERG 衛(wèi)星降水?dāng)?shù)據(jù)與實(shí)測降水?dāng)?shù)據(jù)的CC值分別達(dá)到0.927、0.949,對于偏離程度,GSMap 衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測數(shù)據(jù)相比具有一定程度的低估Bias=-0.124,Rmse=16.989,而GPM-IMERG 衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測數(shù)據(jù)相比則有一定程度的高估Bias=0.478,Rmse=42.144。兩種衛(wèi)星降水產(chǎn)品日尺度數(shù)據(jù)集與實(shí)測數(shù)據(jù)相關(guān)性不顯著,相關(guān)系數(shù)均在0.5 左右,但在偏離程度方面,GSMap 與GPM-IMERG衛(wèi)星降水?dāng)?shù)據(jù)與實(shí)測數(shù)據(jù)相比均有一定程度的高估,Bias分別為0.033 與0.378,Rmse分別為2.287 與2.681。兩種衛(wèi)星降水的偏離程度(Bias)具有較大的差異,可能是由于兩種衛(wèi)星降水產(chǎn)品的反演算法不同的原因造成,GSMap_MVK 采用雙向云移動(dòng)矢量法對源數(shù)據(jù)進(jìn)行處理,GPM-IMERG_Final是全球降水氣候中心站點(diǎn)數(shù)據(jù)校正后的數(shù)據(jù)集。兩種衛(wèi)星降水日尺度數(shù)據(jù)與地面實(shí)測數(shù)據(jù)一致性表現(xiàn)較差主要由于衛(wèi)星觀測降水事件受地形、氣候等多種原因的影響顯著,由于雅江位于高寒地區(qū),地形復(fù)雜,對衛(wèi)星監(jiān)測降水增加了難度。

    圖2 月尺度GSMap與GPM-IMERG 降水產(chǎn)品與實(shí)測數(shù)據(jù)散點(diǎn)密度圖Fig.2 Monthly scale GSMap,GPM-IMERG precipitation products and measured data scatter density map

    圖3 日尺度GSMap與GPM-IMERG 降水產(chǎn)品與實(shí)測數(shù)據(jù)散點(diǎn)密度圖Fig.3 Daily scale GSMap,GPM-IMERG precipitation products and scatter density map

    雅江流域降水主要集中在6-9月,夏季、秋季降水量較大,降雨次數(shù)較多,而春季和冬季降水量很少,對夏、秋兩季衛(wèi)星降水產(chǎn)品日尺度數(shù)據(jù)進(jìn)行定量指標(biāo)分析(圖4),發(fā)現(xiàn),夏季GPM IMERG 衛(wèi)星降水與地面實(shí)測降水的相關(guān)性與精確度優(yōu)于GSMap 衛(wèi)星降水產(chǎn)品,而秋季GSMap 衛(wèi)星降水與地面實(shí)測降水的相關(guān)性與精確度優(yōu)于GPM IMERG。參考夏季、秋季在對雅江流域進(jìn)行降水分析結(jié)果,在以后進(jìn)行降水分析時(shí)可根據(jù)不同季節(jié)選取不同衛(wèi)星降水產(chǎn)品。

    圖4 GSMap與GPM IMERG衛(wèi)星降水產(chǎn)品在春夏兩季與實(shí)測降水?dāng)?shù)據(jù)散點(diǎn)圖Fig.4 Scatter plot of precipitation data of GSMap,GPM IMERG satellite precipitation products and measured data in spring and summer

    黃浠[17]曾對CRU、GLADS 兩種降水產(chǎn)品進(jìn)行評(píng)估,兩種降水產(chǎn)品所測數(shù)據(jù)與地面實(shí)測數(shù)據(jù)相關(guān)性分別為CC=0.24,CC=0.35,劉江濤[3]曾對TRMM 3B42與PERSIAN 在雅江流域進(jìn)行定量評(píng)估,計(jì)算得出其與實(shí)測數(shù)據(jù)偏離程度分別為Bias=0.579、Bias=0.845,呂洋[18]對TRMM 3B42 與3B43 進(jìn)行全流域一致性評(píng)估,計(jì)算得出月尺度衛(wèi)星降水?dāng)?shù)據(jù)與實(shí)測數(shù)據(jù)CC=0.902,日尺度衛(wèi)星降水?dāng)?shù)據(jù)與實(shí)測數(shù)據(jù)CC=0.465。孫赫[9]等對APHRODITE、GPM、PERSIANN-CDR、GLDAS、HAR 五種降水產(chǎn)品日尺度數(shù)據(jù)集在雅江子流域的適用性進(jìn)行分析,GPM 降水產(chǎn)品相關(guān)性與偏離程度表現(xiàn)明顯優(yōu)于PERSIANN-CDR、GLDAS、HAR,與APHRODITE 相比,在不同子流域表現(xiàn)不同。GPM 衛(wèi)星降水產(chǎn)品作為TRMM 的升級(jí)版本,結(jié)合本文分析結(jié)果,可以看出,GSMap 與GPM-IMERG 衛(wèi)星降水產(chǎn)品在精確度方面明顯優(yōu)于TRMM(Bias<0.579),兩種衛(wèi)星降水產(chǎn)品的數(shù)據(jù)與實(shí)測數(shù)據(jù)相關(guān)性優(yōu)于TRMM(月尺度CC>0.902,日尺度CC>0.465)。與TRMM衛(wèi)星降水?dāng)?shù)據(jù)相比,GSMap 與GPM-IMERG 均能更加準(zhǔn)確的對實(shí)際降水情況進(jìn)行估計(jì),更好地滿足水文、氣象等領(lǐng)域?qū)τ跀?shù)據(jù)精度的要求。

    3.1.2 空間尺度

    兩種衛(wèi)星降水產(chǎn)品與地面實(shí)測柵格數(shù)據(jù)集的相關(guān)性在空間上分布趨勢一致(圖5),均為在雅江中游相關(guān)性顯著,而上游與下游表現(xiàn)相對較差。對于兩種產(chǎn)品月尺度降水?dāng)?shù)據(jù),GSMap衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測數(shù)據(jù)的相關(guān)系數(shù)范圍為0.177~0.954,GPM-IMERG 衛(wèi)星降水與地面實(shí)測數(shù)據(jù)的相關(guān)系數(shù)范圍為0.293~0.993,GSMap 衛(wèi)星降水與IMERG 衛(wèi)星降水日尺度數(shù)據(jù)與地面實(shí)測數(shù)據(jù)相關(guān)系數(shù)分布范圍分別為0.001~0.420、0.003~0.553。分析認(rèn)為雅江流域上游與下游衛(wèi)星降水產(chǎn)品與地面實(shí)測柵格數(shù)據(jù)集的相關(guān)性較差,主要是因?yàn)閷?shí)測柵格數(shù)據(jù)集由地面測站監(jiān)測降水?dāng)?shù)據(jù)通過樣條插值得到,而上下游區(qū)域地面觀測站點(diǎn)設(shè)置稀疏,在插值過程中產(chǎn)生一定程度的誤差。

    圖5 雅江流域衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測數(shù)據(jù)相關(guān)系數(shù)空間分布圖Fig.5 The spatial distribution of correlation coefficients between satellite precipitation data and ground-measured data in the Yarlung Zangbo River Basin

    3.2 單一站點(diǎn)數(shù)據(jù)精度評(píng)定

    通過對單一站點(diǎn)相關(guān)性與偏差的分析,可以進(jìn)一步驗(yàn)證時(shí)空精度評(píng)估的結(jié)果(表3、4),7 個(gè)地面觀測站對應(yīng)的GSMap 與GPM-IMERG 兩種衛(wèi)星降水產(chǎn)品月尺度數(shù)據(jù)與地面實(shí)測柵格數(shù)據(jù)相關(guān)性表現(xiàn)較好,在雅江中游區(qū)域,各站點(diǎn)衛(wèi)星降水量數(shù)據(jù)與地面實(shí)測數(shù)據(jù)相關(guān)系數(shù)CC均大于0.8(表3);相反,日尺度的衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測柵格數(shù)據(jù)的相關(guān)性較差(表4),除波密站外,CC值均小于0.5。由7 個(gè)實(shí)測站的相對偏差Bias的范圍為0.007 至0.993 不等,可以看出Bias隨地理位置的變化,表現(xiàn)出明顯的差異性。

    表3 月尺度兩種衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測數(shù)據(jù)相關(guān)性、相對偏差統(tǒng)計(jì)表Tab.3 Statistical table of correlation and relative deviation between two kinds of satellite precipitation data and ground-measured data at monthly scale

    表4 日尺度兩種衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測數(shù)據(jù)相關(guān)性、相對偏差統(tǒng)計(jì)表Tab.4 Statistical table of correlation and relative deviation between two kinds of satellite precipitation data and ground-measured data at daily scale

    3.3 分類指標(biāo)評(píng)定

    在定量分析的基礎(chǔ)上進(jìn)行衛(wèi)星降水?dāng)?shù)據(jù)產(chǎn)品的分類指標(biāo)分析,可以更好的反應(yīng)衛(wèi)星降水產(chǎn)品對降水時(shí)間的監(jiān)測捕捉能力。本文對分類指標(biāo)中的探照率(POD)與錯(cuò)報(bào)率(FAR)進(jìn)行分析,GSMap 衛(wèi)星降水產(chǎn)品的POD為0.958,GPM_IMERG 衛(wèi)星降水產(chǎn)品的POD為0.924,兩種衛(wèi)星降水產(chǎn)品在探照率上表現(xiàn)差異不大,說明兩種衛(wèi)星降水產(chǎn)品均能夠探測到大部分的降水事件;GSMap與GPM-IMERG 的FAR分別為0.100與0.172,說明兩種衛(wèi)星降水產(chǎn)品對降水事件的錯(cuò)誤判斷情況相似,僅有小部分降水事件被錯(cuò)誤判斷。

    表5 分類分析結(jié)果表Tab.5 Table of classification analysis results table

    4 結(jié) 論

    選取定量與定性分析相結(jié)合的方式構(gòu)建適用性評(píng)估體系,分別在日尺度、月尺度上對GPM 系列降水產(chǎn)品中GSMap 與GPM IMERG 的精度進(jìn)行評(píng)估,并對兩種衛(wèi)星降水產(chǎn)品在雅江流域的適用性進(jìn)行分析,主要得出以下結(jié)論。

    (1)從全流域角度來看,GPM-IMERG 衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測數(shù)據(jù)的相關(guān)性略優(yōu)于GsMap,而GSMap 衛(wèi)星降水?dāng)?shù)據(jù)在偏離程度上表現(xiàn)優(yōu)于GPM-IMERG。GSMap衛(wèi)星降水產(chǎn)品與PMIMERG 衛(wèi)星降水產(chǎn)品的月尺度降水?dāng)?shù)據(jù)與實(shí)際觀測降水?dāng)?shù)據(jù)具有較高的相關(guān)性,但日尺度降水?dāng)?shù)據(jù)與地面觀測站點(diǎn)所測降水?dāng)?shù)據(jù)相關(guān)性較差,建議通過與地面觀測數(shù)據(jù)融合后使用;兩種數(shù)據(jù)產(chǎn)品精度存在季節(jié)性差異,夏季推薦采用GPM IMERG衛(wèi)星降水產(chǎn)品,而秋季推薦采用GSMap衛(wèi)星降水產(chǎn)品。

    (2)兩種衛(wèi)星降水的日尺度、月尺度數(shù)據(jù)集均表明雅江流域中游衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測網(wǎng)格數(shù)據(jù)相關(guān)性表現(xiàn)較好,而上游與下游衛(wèi)星降水?dāng)?shù)據(jù)與地面實(shí)測網(wǎng)格數(shù)據(jù)相關(guān)性不顯著。

    (3)通過選擇的7 個(gè)地面實(shí)測站點(diǎn)進(jìn)行相關(guān)性與相對偏差分析,進(jìn)一步驗(yàn)證了兩種衛(wèi)星降水產(chǎn)品月尺度數(shù)據(jù)集與地面實(shí)測站點(diǎn)相關(guān)性顯著,而日尺度數(shù)據(jù)集相關(guān)性較差的結(jié)論,同時(shí)發(fā)現(xiàn)兩種衛(wèi)星降水的Bias隨地理位置的變化,表現(xiàn)出明顯的差異性。

    (4)GSMap 與GPM-IMERG 兩種衛(wèi)星降水產(chǎn)品的探照率和錯(cuò)報(bào)率均表現(xiàn)較好,均能準(zhǔn)確的探測到大部分降水事件,且僅錯(cuò)誤判斷少部分降水事件。

    (5)在雅魯藏布江流域,GPM 衛(wèi)星降水?dāng)?shù)據(jù)與TRMM 衛(wèi)星降水?dāng)?shù)據(jù)相比,具有更好的適用性,更能滿足水資源開發(fā)利用對數(shù)據(jù)的需求。

    綜上,GPM 系列降水產(chǎn)品延續(xù)了TRMM 在時(shí)空精度上的優(yōu)勢,并具有更好的適用性,但在雅江的水文模擬及預(yù)報(bào)研究中,還要考慮季節(jié)和上下游差異,通過與地面觀測數(shù)據(jù)的融合分析,為水文模擬與預(yù)報(bào)提供更為可靠的降水?dāng)?shù)據(jù)支撐。

    猜你喜歡
    雅江尺度降水
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    雅江光電舉辦“光質(zhì)造 創(chuàng)未來”2018銷售年會(huì)
    演藝科技(2018年2期)2018-03-26 12:30:58
    雅江亮相2017年中國文化館年會(huì)·文化藝術(shù)博覽會(huì)
    演藝科技(2017年12期)2018-01-27 00:10:55
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    蓮花盛放雅江點(diǎn)亮杭州奧體中心主體育場
    中國照明(2016年5期)2016-06-15 20:30:13
    降水現(xiàn)象儀模擬軟件設(shè)計(jì)與實(shí)現(xiàn)
    “雅江”在廣州展示LED產(chǎn)品
    演藝科技(2016年6期)2016-05-14 12:10:33
    ESSENTIAL NORMS OF PRODUCTS OF WEIGHTED COMPOSITION OPERATORS AND DIFFERENTIATION OPERATORS BETWEEN BANACH SPACES OF ANALYTIC FUNCTIONS?
    日韩精品有码人妻一区| 噜噜噜噜噜久久久久久91| 久久99热这里只频精品6学生| videossex国产| 干丝袜人妻中文字幕| 成人综合一区亚洲| 中文天堂在线官网| 啦啦啦在线观看免费高清www| 久久午夜福利片| 男人狂女人下面高潮的视频| 色5月婷婷丁香| 交换朋友夫妻互换小说| 自拍欧美九色日韩亚洲蝌蚪91 | 看非洲黑人一级黄片| 在线观看国产h片| 国产精品一区www在线观看| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩东京热| 欧美亚洲 丝袜 人妻 在线| 国产一区二区亚洲精品在线观看| 性色avwww在线观看| 一二三四中文在线观看免费高清| 天天躁日日操中文字幕| 毛片一级片免费看久久久久| 在线免费十八禁| 麻豆成人午夜福利视频| 狂野欧美白嫩少妇大欣赏| 毛片女人毛片| 下体分泌物呈黄色| 日本黄色片子视频| 伊人久久国产一区二区| 精品久久久噜噜| 老女人水多毛片| 欧美潮喷喷水| 成人欧美大片| 国产精品99久久久久久久久| 伦理电影大哥的女人| 永久免费av网站大全| 亚洲av男天堂| 国产亚洲精品久久久com| 高清在线视频一区二区三区| 免费av不卡在线播放| 亚洲欧美成人综合另类久久久| av国产精品久久久久影院| 久久久久性生活片| 岛国毛片在线播放| 免费电影在线观看免费观看| 热99国产精品久久久久久7| 乱码一卡2卡4卡精品| 在线a可以看的网站| 精品一区二区免费观看| 插阴视频在线观看视频| 精品视频人人做人人爽| 免费看a级黄色片| 青春草亚洲视频在线观看| 国产黄色视频一区二区在线观看| 国产精品99久久久久久久久| 26uuu在线亚洲综合色| 久久6这里有精品| 在线亚洲精品国产二区图片欧美 | a级一级毛片免费在线观看| 国产成人精品久久久久久| 黑人高潮一二区| 2021少妇久久久久久久久久久| 久久久a久久爽久久v久久| 国产精品久久久久久久电影| 大片电影免费在线观看免费| 国产乱来视频区| 亚洲精品一二三| 国产极品天堂在线| 国产成人免费无遮挡视频| 九九爱精品视频在线观看| 伊人久久精品亚洲午夜| 在线观看一区二区三区| 热99国产精品久久久久久7| 亚洲在线观看片| 麻豆久久精品国产亚洲av| 大片电影免费在线观看免费| 97精品久久久久久久久久精品| 99热这里只有精品一区| 一级爰片在线观看| 成人午夜精彩视频在线观看| 欧美人与善性xxx| 国产成人精品福利久久| 国产有黄有色有爽视频| 成人毛片60女人毛片免费| 精品久久久久久久久亚洲| 一边亲一边摸免费视频| 亚洲色图综合在线观看| 麻豆成人午夜福利视频| 97精品久久久久久久久久精品| 如何舔出高潮| 成人综合一区亚洲| 色5月婷婷丁香| 免费观看在线日韩| 婷婷色综合大香蕉| 777米奇影视久久| 五月玫瑰六月丁香| 国产精品爽爽va在线观看网站| 好男人视频免费观看在线| 直男gayav资源| 在线天堂最新版资源| eeuss影院久久| 日本与韩国留学比较| 亚洲av免费在线观看| av在线播放精品| 在现免费观看毛片| 国产精品一二三区在线看| 欧美潮喷喷水| 中文资源天堂在线| 内射极品少妇av片p| 精华霜和精华液先用哪个| 久久6这里有精品| 精品人妻一区二区三区麻豆| 日韩av在线免费看完整版不卡| 久久影院123| 女人被狂操c到高潮| 亚洲欧美日韩东京热| 搡女人真爽免费视频火全软件| 日韩免费高清中文字幕av| 日韩电影二区| 久久国产乱子免费精品| 毛片女人毛片| 国产精品99久久久久久久久| 亚洲精品乱久久久久久| 在线免费观看不下载黄p国产| 尾随美女入室| 七月丁香在线播放| 国产中年淑女户外野战色| 丝袜喷水一区| 日韩欧美一区视频在线观看 | 国产 精品1| 久久久a久久爽久久v久久| 国产欧美亚洲国产| 少妇人妻一区二区三区视频| 18禁裸乳无遮挡动漫免费视频 | 国产伦精品一区二区三区四那| 我的老师免费观看完整版| 91久久精品国产一区二区三区| 啦啦啦啦在线视频资源| 肉色欧美久久久久久久蜜桃 | 黄色一级大片看看| 亚洲天堂av无毛| 好男人视频免费观看在线| 欧美97在线视频| 一区二区三区精品91| 亚洲aⅴ乱码一区二区在线播放| videos熟女内射| 在线观看av片永久免费下载| 亚洲aⅴ乱码一区二区在线播放| tube8黄色片| 日韩人妻高清精品专区| 嫩草影院入口| 久久久精品免费免费高清| av网站免费在线观看视频| 丰满人妻一区二区三区视频av| 日韩成人av中文字幕在线观看| 久久久久性生活片| 成人欧美大片| 国产女主播在线喷水免费视频网站| 最新中文字幕久久久久| 国产男女内射视频| 日本爱情动作片www.在线观看| 亚洲国产色片| 真实男女啪啪啪动态图| 国产毛片a区久久久久| 亚洲精品亚洲一区二区| av免费在线看不卡| 欧美zozozo另类| 午夜激情福利司机影院| 夜夜爽夜夜爽视频| .国产精品久久| 亚洲天堂av无毛| 美女xxoo啪啪120秒动态图| 性色avwww在线观看| 在线亚洲精品国产二区图片欧美 | 日日撸夜夜添| 亚洲精品456在线播放app| 99视频精品全部免费 在线| 五月玫瑰六月丁香| 嫩草影院新地址| 亚洲精品色激情综合| 免费观看无遮挡的男女| 美女高潮的动态| 青春草视频在线免费观看| 免费av观看视频| 亚洲电影在线观看av| 日韩精品有码人妻一区| 欧美成人一区二区免费高清观看| 中国美白少妇内射xxxbb| 草草在线视频免费看| 国产精品久久久久久精品电影小说 | 欧美3d第一页| 交换朋友夫妻互换小说| 最后的刺客免费高清国语| 联通29元200g的流量卡| 国产一区二区在线观看日韩| 少妇 在线观看| 国产精品一区二区在线观看99| 国产亚洲午夜精品一区二区久久 | 精品午夜福利在线看| 亚洲成人一二三区av| 日本-黄色视频高清免费观看| 欧美性猛交╳xxx乱大交人| 97在线人人人人妻| 亚洲av中文av极速乱| 六月丁香七月| 99re6热这里在线精品视频| 国产精品国产三级国产av玫瑰| 秋霞在线观看毛片| 亚洲国产高清在线一区二区三| av在线天堂中文字幕| 国产 一区 欧美 日韩| 免费看a级黄色片| 美女cb高潮喷水在线观看| 91狼人影院| 一本色道久久久久久精品综合| 啦啦啦啦在线视频资源| 欧美精品人与动牲交sv欧美| 最近手机中文字幕大全| 伦精品一区二区三区| 日韩三级伦理在线观看| 亚洲美女搞黄在线观看| 97人妻精品一区二区三区麻豆| 日韩亚洲欧美综合| av福利片在线观看| 久久6这里有精品| av专区在线播放| 国内揄拍国产精品人妻在线| 美女高潮的动态| 身体一侧抽搐| 插阴视频在线观看视频| 久久久久久久久久久免费av| 日韩成人av中文字幕在线观看| 成人综合一区亚洲| 免费电影在线观看免费观看| 丝瓜视频免费看黄片| 国产高潮美女av| 我要看日韩黄色一级片| 麻豆乱淫一区二区| 赤兔流量卡办理| 99视频精品全部免费 在线| 亚州av有码| 男插女下体视频免费在线播放| 国产精品蜜桃在线观看| 天天一区二区日本电影三级| 国产亚洲91精品色在线| 久久精品夜色国产| 免费看不卡的av| 国产亚洲5aaaaa淫片| 国产av不卡久久| 高清欧美精品videossex| 精品亚洲乱码少妇综合久久| 国产老妇女一区| 亚洲内射少妇av| 欧美日韩亚洲高清精品| 欧美三级亚洲精品| 校园人妻丝袜中文字幕| 久久久精品免费免费高清| 日日撸夜夜添| 久久精品国产鲁丝片午夜精品| 自拍偷自拍亚洲精品老妇| 五月天丁香电影| 黄色日韩在线| 久久6这里有精品| 波多野结衣巨乳人妻| 日韩大片免费观看网站| 欧美高清性xxxxhd video| 久久99精品国语久久久| 亚洲精品第二区| 日本免费在线观看一区| 一区二区三区免费毛片| 三级经典国产精品| 噜噜噜噜噜久久久久久91| 夜夜爽夜夜爽视频| 日韩大片免费观看网站| 欧美极品一区二区三区四区| 亚洲成人中文字幕在线播放| 成年免费大片在线观看| 久久久亚洲精品成人影院| 久久久久九九精品影院| 99精国产麻豆久久婷婷| 免费观看无遮挡的男女| 亚洲欧美日韩东京热| 亚洲最大成人av| 欧美一级a爱片免费观看看| 亚洲怡红院男人天堂| 成人二区视频| 少妇人妻一区二区三区视频| 少妇 在线观看| 亚洲av日韩在线播放| 久久午夜福利片| 91狼人影院| 亚洲综合色惰| 91精品伊人久久大香线蕉| 99精国产麻豆久久婷婷| 精品一区在线观看国产| av线在线观看网站| 亚洲av欧美aⅴ国产| 丝袜美腿在线中文| 免费看不卡的av| 国产精品一及| 三级国产精品片| 国产精品一区二区在线观看99| 国产男人的电影天堂91| 三级国产精品片| 制服丝袜香蕉在线| 午夜福利视频精品| 午夜免费观看性视频| 免费av毛片视频| 国产精品久久久久久精品古装| 免费观看在线日韩| 亚洲色图综合在线观看| 国产在线男女| 18禁裸乳无遮挡动漫免费视频 | 九九在线视频观看精品| 久久精品综合一区二区三区| 成年av动漫网址| 麻豆久久精品国产亚洲av| 一二三四中文在线观看免费高清| 日韩欧美一区视频在线观看 | 国产伦精品一区二区三区四那| 亚洲精品日韩av片在线观看| 高清欧美精品videossex| 免费看光身美女| 国产黄a三级三级三级人| 九九爱精品视频在线观看| 日日啪夜夜撸| 69av精品久久久久久| 久久久久久久亚洲中文字幕| 亚洲熟女精品中文字幕| 亚洲精品456在线播放app| 久久99蜜桃精品久久| 久久人人爽人人片av| 久久久亚洲精品成人影院| 亚洲欧美中文字幕日韩二区| 99视频精品全部免费 在线| 丝袜美腿在线中文| 网址你懂的国产日韩在线| 免费观看性生交大片5| 色5月婷婷丁香| 嫩草影院精品99| 午夜免费男女啪啪视频观看| 免费av毛片视频| 久久国内精品自在自线图片| 午夜亚洲福利在线播放| 亚洲av成人精品一区久久| 免费观看a级毛片全部| 国产黄色视频一区二区在线观看| 成人高潮视频无遮挡免费网站| 久久久久精品久久久久真实原创| 天天躁日日操中文字幕| 国产亚洲av片在线观看秒播厂| 身体一侧抽搐| 视频中文字幕在线观看| 国产大屁股一区二区在线视频| 久久99蜜桃精品久久| a级毛色黄片| 午夜福利在线在线| 成年免费大片在线观看| 日本午夜av视频| 欧美 日韩 精品 国产| 哪个播放器可以免费观看大片| 美女高潮的动态| 黄片无遮挡物在线观看| 在线亚洲精品国产二区图片欧美 | 国产美女午夜福利| 日日啪夜夜撸| 国产一级毛片在线| 日韩不卡一区二区三区视频在线| 伊人久久国产一区二区| 亚洲自偷自拍三级| 人妻 亚洲 视频| 一边亲一边摸免费视频| 国产精品国产av在线观看| 97在线视频观看| 国产探花在线观看一区二区| 看非洲黑人一级黄片| 久久综合国产亚洲精品| 亚洲成人精品中文字幕电影| videossex国产| 女人久久www免费人成看片| 99热网站在线观看| 一级黄片播放器| 成人毛片60女人毛片免费| av在线播放精品| 观看免费一级毛片| 国产黄频视频在线观看| 大香蕉97超碰在线| 日本熟妇午夜| 日韩av免费高清视频| 热99国产精品久久久久久7| 白带黄色成豆腐渣| 18禁在线播放成人免费| 国产色婷婷99| 久久久久精品久久久久真实原创| 97精品久久久久久久久久精品| 亚洲av欧美aⅴ国产| 大片免费播放器 马上看| 国产成人一区二区在线| 日韩av免费高清视频| 国产免费又黄又爽又色| 国产伦精品一区二区三区视频9| 天堂俺去俺来也www色官网| 午夜福利高清视频| 精品久久久久久久末码| 国产毛片a区久久久久| 精品少妇久久久久久888优播| 熟女av电影| 一级毛片电影观看| 老女人水多毛片| 欧美bdsm另类| 欧美成人精品欧美一级黄| 国产真实伦视频高清在线观看| 午夜激情福利司机影院| 身体一侧抽搐| 六月丁香七月| 色综合色国产| 丝袜美腿在线中文| 99re6热这里在线精品视频| 国产成年人精品一区二区| 观看免费一级毛片| 国内揄拍国产精品人妻在线| 成年免费大片在线观看| 一区二区三区精品91| 天堂俺去俺来也www色官网| 成人毛片60女人毛片免费| 嫩草影院入口| 国产欧美亚洲国产| 亚洲,一卡二卡三卡| 国产一区二区三区综合在线观看 | 在线观看av片永久免费下载| 免费观看的影片在线观看| 欧美日韩综合久久久久久| 国产亚洲最大av| 成人午夜精彩视频在线观看| 男女国产视频网站| 成年女人在线观看亚洲视频 | 成年av动漫网址| 97人妻精品一区二区三区麻豆| 免费看日本二区| 国产精品成人在线| 男插女下体视频免费在线播放| 你懂的网址亚洲精品在线观看| 亚洲最大成人av| 看黄色毛片网站| 精品一区在线观看国产| 小蜜桃在线观看免费完整版高清| 一边亲一边摸免费视频| 午夜免费观看性视频| a级毛片免费高清观看在线播放| 国产伦精品一区二区三区四那| 一级毛片电影观看| a级一级毛片免费在线观看| 又粗又硬又长又爽又黄的视频| 又爽又黄a免费视频| 我的女老师完整版在线观看| 涩涩av久久男人的天堂| 国产成人免费观看mmmm| 观看免费一级毛片| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久亚洲中文字幕| av免费在线看不卡| 中文在线观看免费www的网站| 五月伊人婷婷丁香| 18禁动态无遮挡网站| 亚洲欧美日韩无卡精品| 国产视频内射| 亚洲一区二区三区欧美精品 | 国产亚洲5aaaaa淫片| 男人狂女人下面高潮的视频| 97人妻精品一区二区三区麻豆| 在线亚洲精品国产二区图片欧美 | 麻豆成人午夜福利视频| 免费av不卡在线播放| 麻豆国产97在线/欧美| 人人妻人人爽人人添夜夜欢视频 | av免费在线看不卡| 80岁老熟妇乱子伦牲交| 久久97久久精品| 夜夜看夜夜爽夜夜摸| 狂野欧美激情性bbbbbb| 最新中文字幕久久久久| 精品视频人人做人人爽| 最近中文字幕高清免费大全6| 人体艺术视频欧美日本| 秋霞在线观看毛片| 久久99热这里只有精品18| 精品视频人人做人人爽| 91在线精品国自产拍蜜月| 亚洲国产高清在线一区二区三| 午夜激情福利司机影院| 免费高清在线观看视频在线观看| 亚洲人成网站在线观看播放| xxx大片免费视频| 又黄又爽又刺激的免费视频.| 午夜激情久久久久久久| 成人欧美大片| 在线观看国产h片| 99re6热这里在线精品视频| 一个人观看的视频www高清免费观看| h日本视频在线播放| 在线观看免费高清a一片| 夫妻性生交免费视频一级片| 18禁在线无遮挡免费观看视频| 亚洲成色77777| 18禁裸乳无遮挡免费网站照片| 在线看a的网站| 男人爽女人下面视频在线观看| 久久久午夜欧美精品| 国产一级毛片在线| 视频区图区小说| 最近最新中文字幕大全电影3| 国产欧美另类精品又又久久亚洲欧美| 日韩一区二区视频免费看| 舔av片在线| 日韩国内少妇激情av| 亚洲,一卡二卡三卡| 国产高清不卡午夜福利| 亚洲va在线va天堂va国产| 国产有黄有色有爽视频| 久久久色成人| 男女边摸边吃奶| 99热6这里只有精品| 街头女战士在线观看网站| 成人鲁丝片一二三区免费| 97人妻精品一区二区三区麻豆| 日本三级黄在线观看| 国产一区亚洲一区在线观看| 午夜免费观看性视频| 搡老乐熟女国产| 搞女人的毛片| 亚洲国产精品成人综合色| 一个人看的www免费观看视频| 精品久久久久久久末码| 波野结衣二区三区在线| 亚洲国产色片| 精品国产三级普通话版| 国产午夜精品一二区理论片| 97热精品久久久久久| 欧美亚洲 丝袜 人妻 在线| 乱码一卡2卡4卡精品| 欧美一区二区亚洲| 欧美日韩一区二区视频在线观看视频在线 | 午夜福利视频1000在线观看| 国产午夜精品久久久久久一区二区三区| 97在线人人人人妻| 一级二级三级毛片免费看| 夜夜看夜夜爽夜夜摸| 成人鲁丝片一二三区免费| 大码成人一级视频| www.av在线官网国产| 久久久久精品性色| 国产精品av视频在线免费观看| 亚洲人成网站在线观看播放| 日韩欧美精品免费久久| 黄片无遮挡物在线观看| 成人综合一区亚洲| 久久韩国三级中文字幕| 亚洲精品久久久久久婷婷小说| av黄色大香蕉| 欧美性感艳星| 大话2 男鬼变身卡| 99热这里只有精品一区| 国产精品爽爽va在线观看网站| 99视频精品全部免费 在线| 啦啦啦啦在线视频资源| 久久久国产一区二区| 国产黄片美女视频| 免费在线观看成人毛片| 国产欧美另类精品又又久久亚洲欧美| 久久鲁丝午夜福利片| 久久久久精品性色| 亚洲伊人久久精品综合| 国产精品99久久久久久久久| 国产精品久久久久久精品电影小说 | 亚洲自偷自拍三级| 黄色一级大片看看| 久久人人爽人人片av| 免费不卡的大黄色大毛片视频在线观看| 久久精品国产亚洲av涩爱| 五月天丁香电影| 日韩不卡一区二区三区视频在线| 日韩在线高清观看一区二区三区| 禁无遮挡网站| 99热国产这里只有精品6| 18+在线观看网站| 春色校园在线视频观看| 婷婷色综合大香蕉| 看黄色毛片网站| 久久国内精品自在自线图片| 99热国产这里只有精品6| kizo精华| 亚洲精品日本国产第一区| 草草在线视频免费看| av在线亚洲专区| 五月天丁香电影| 99久久精品热视频| 99re6热这里在线精品视频| 黄色视频在线播放观看不卡| 七月丁香在线播放| 亚洲国产日韩一区二区| 亚洲精品第二区| 大码成人一级视频| 一个人看的www免费观看视频| 99热这里只有是精品在线观看| 大片电影免费在线观看免费| 三级男女做爰猛烈吃奶摸视频| 成人二区视频| 丰满乱子伦码专区| 亚洲第一区二区三区不卡| 熟妇人妻不卡中文字幕| av在线天堂中文字幕| 91久久精品电影网| 国产综合懂色|