• <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?
    能在线免费看毛片的网站| 精品久久久久久久久av| 精品久久国产蜜桃| 国产伦理片在线播放av一区 | 久久久久久久久大av| 99久久中文字幕三级久久日本| 免费观看人在逋| 狂野欧美白嫩少妇大欣赏| 干丝袜人妻中文字幕| 99久久精品一区二区三区| 欧美日本视频| 欧美日韩在线观看h| 99久久中文字幕三级久久日本| 久久韩国三级中文字幕| 精品久久久久久成人av| 网址你懂的国产日韩在线| 91aial.com中文字幕在线观看| 成年女人永久免费观看视频| 国产黄色视频一区二区在线观看 | 日本黄色视频三级网站网址| 久久亚洲精品不卡| 一边摸一边抽搐一进一小说| 国产又黄又爽又无遮挡在线| 婷婷精品国产亚洲av| 大型黄色视频在线免费观看| 日韩强制内射视频| 国产精品一二三区在线看| 午夜视频国产福利| 午夜福利在线在线| 亚洲无线在线观看| 麻豆av噜噜一区二区三区| 麻豆一二三区av精品| 国产精品一区二区性色av| 我的女老师完整版在线观看| 国产亚洲精品久久久com| 久久人人爽人人片av| 国产精品野战在线观看| 国产乱人视频| 麻豆国产97在线/欧美| 99久国产av精品国产电影| 午夜久久久久精精品| av天堂中文字幕网| 亚州av有码| 日韩高清综合在线| 久久久精品大字幕| 熟女电影av网| 日韩三级伦理在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站在线播| 久久鲁丝午夜福利片| 成人三级黄色视频| 夜夜爽天天搞| 中文在线观看免费www的网站| 亚洲经典国产精华液单| 美女黄网站色视频| 亚洲国产色片| 日韩在线高清观看一区二区三区| 中文欧美无线码| 婷婷六月久久综合丁香| 少妇裸体淫交视频免费看高清| 久久99热6这里只有精品| 亚洲欧美成人精品一区二区| 国产一区二区三区av在线 | 国产在线男女| 婷婷亚洲欧美| 伦精品一区二区三区| 男女边吃奶边做爰视频| 午夜福利在线在线| 国产成人精品久久久久久| 99久久精品一区二区三区| 综合色av麻豆| 一级黄片播放器| 麻豆一二三区av精品| 精品人妻偷拍中文字幕| 偷拍熟女少妇极品色| 欧美色视频一区免费| 又黄又爽又刺激的免费视频.| 美女脱内裤让男人舔精品视频 | 青春草亚洲视频在线观看| 午夜福利视频1000在线观看| 美女被艹到高潮喷水动态| 国产女主播在线喷水免费视频网站 | 免费不卡的大黄色大毛片视频在线观看 | 男女下面进入的视频免费午夜| 99九九线精品视频在线观看视频| 热99re8久久精品国产| 赤兔流量卡办理| 精品人妻偷拍中文字幕| 日本免费a在线| 欧美丝袜亚洲另类| 午夜视频国产福利| 免费观看a级毛片全部| 校园人妻丝袜中文字幕| 亚洲美女搞黄在线观看| 国产午夜福利久久久久久| 国产精品国产三级国产av玫瑰| 午夜免费激情av| 国产亚洲av片在线观看秒播厂 | 一进一出抽搐动态| 久久久久久久久久成人| 国产精品av视频在线免费观看| 久久久久久国产a免费观看| 乱系列少妇在线播放| 成人一区二区视频在线观看| 大型黄色视频在线免费观看| av视频在线观看入口| 午夜精品国产一区二区电影 | 日本爱情动作片www.在线观看| 一级av片app| 麻豆成人av视频| 国产午夜精品久久久久久一区二区三区| 色哟哟·www| 欧美xxxx性猛交bbbb| 在线国产一区二区在线| 黄片wwwwww| 99久久无色码亚洲精品果冻| 啦啦啦韩国在线观看视频| 久久久久久九九精品二区国产| a级毛片a级免费在线| 久久人人爽人人片av| 欧美+日韩+精品| av免费观看日本| 欧美激情在线99| 国产久久久一区二区三区| 成人亚洲欧美一区二区av| 此物有八面人人有两片| av在线天堂中文字幕| 女人被狂操c到高潮| 成人高潮视频无遮挡免费网站| 一区福利在线观看| 好男人在线观看高清免费视频| 啦啦啦韩国在线观看视频| 亚洲精品日韩在线中文字幕 | 成人特级黄色片久久久久久久| 99在线人妻在线中文字幕| 三级毛片av免费| av在线亚洲专区| 婷婷色综合大香蕉| 99久国产av精品国产电影| 免费一级毛片在线播放高清视频| 亚洲国产精品sss在线观看| 亚洲欧美日韩卡通动漫| 国产成人a区在线观看| 一区二区三区高清视频在线| 亚洲国产精品sss在线观看| 六月丁香七月| 天堂av国产一区二区熟女人妻| 12—13女人毛片做爰片一| 噜噜噜噜噜久久久久久91| 性欧美人与动物交配| 少妇高潮的动态图| 免费av毛片视频| 亚洲av免费高清在线观看| 小蜜桃在线观看免费完整版高清| 欧美bdsm另类| 毛片女人毛片| 黄色欧美视频在线观看| 国产 一区精品| 嫩草影院入口| 国产在线男女| 亚洲综合色惰| 日本撒尿小便嘘嘘汇集6| 九草在线视频观看| 亚洲18禁久久av| 非洲黑人性xxxx精品又粗又长| 亚洲人成网站在线观看播放| 久久99热这里只有精品18| 久久九九热精品免费| 亚洲三级黄色毛片| 3wmmmm亚洲av在线观看| 亚洲va在线va天堂va国产| 国产av麻豆久久久久久久| 精品人妻视频免费看| 欧美成人一区二区免费高清观看| kizo精华| 久久中文看片网| 亚洲欧美中文字幕日韩二区| 少妇高潮的动态图| 亚洲精品国产成人久久av| 亚洲av电影不卡..在线观看| 成人毛片a级毛片在线播放| a级毛片a级免费在线| 一卡2卡三卡四卡精品乱码亚洲| 国产白丝娇喘喷水9色精品| 国产精品无大码| 听说在线观看完整版免费高清| 中国美女看黄片| 日韩欧美国产在线观看| 天堂中文最新版在线下载 | 欧美一区二区国产精品久久精品| 亚洲四区av| 亚洲人成网站高清观看| 两个人的视频大全免费| 99热网站在线观看| 在线免费十八禁| 级片在线观看| 九九在线视频观看精品| 国产精品久久久久久精品电影| 狠狠狠狠99中文字幕| 99久国产av精品| 婷婷亚洲欧美| 亚洲自偷自拍三级| 能在线免费观看的黄片| 午夜视频国产福利| 91aial.com中文字幕在线观看| 岛国在线免费视频观看| 国产91av在线免费观看| 国产精品国产高清国产av| 草草在线视频免费看| 男女视频在线观看网站免费| 成熟少妇高潮喷水视频| 在线国产一区二区在线| 成年女人看的毛片在线观看| 国产精品av视频在线免费观看| 亚洲欧洲国产日韩| av视频在线观看入口| 51国产日韩欧美| 成人亚洲精品av一区二区| 欧美一级a爱片免费观看看| 夫妻性生交免费视频一级片| www.色视频.com| 18禁裸乳无遮挡免费网站照片| 色尼玛亚洲综合影院| 国产熟女欧美一区二区| 久久久久久伊人网av| 日韩亚洲欧美综合| 中文字幕久久专区| 日产精品乱码卡一卡2卡三| 午夜老司机福利剧场| 中文字幕av在线有码专区| 九九久久精品国产亚洲av麻豆| 最近2019中文字幕mv第一页| .国产精品久久| 国产精品爽爽va在线观看网站| 色尼玛亚洲综合影院| 99国产极品粉嫩在线观看| 一级毛片我不卡| 成人永久免费在线观看视频| 亚洲欧洲日产国产| 又黄又爽又刺激的免费视频.| 国产 一区精品| av女优亚洲男人天堂| www.av在线官网国产| 一边亲一边摸免费视频| 国产 一区 欧美 日韩| 男的添女的下面高潮视频| 日韩av在线大香蕉| 看片在线看免费视频| 2021天堂中文幕一二区在线观| 久久久久九九精品影院| 伦理电影大哥的女人| 久久久久久国产a免费观看| 国产av一区在线观看免费| av天堂在线播放| 亚洲欧美清纯卡通| 国产黄片美女视频| 搡老妇女老女人老熟妇| 欧美高清成人免费视频www| 亚洲av男天堂| 一卡2卡三卡四卡精品乱码亚洲| 日韩av不卡免费在线播放| 97热精品久久久久久| 色5月婷婷丁香| 久久人人爽人人爽人人片va| 亚洲精品乱码久久久v下载方式| 午夜福利在线在线| 国产精品久久久久久精品电影小说 | 一夜夜www| 成人亚洲精品av一区二区| 搞女人的毛片| 国产三级在线视频| 日本免费a在线| 国产在视频线在精品| 国产精品一区二区性色av| 中国国产av一级| 亚洲欧美精品综合久久99| h日本视频在线播放| 国产成年人精品一区二区| 国产人妻一区二区三区在| 久久久久免费精品人妻一区二区| 久久婷婷人人爽人人干人人爱| 一本久久中文字幕| 小蜜桃在线观看免费完整版高清| 精品人妻视频免费看| 欧美性感艳星| 精品熟女少妇av免费看| 久久这里只有精品中国| 成人av在线播放网站| 欧美高清成人免费视频www| www.色视频.com| 最近的中文字幕免费完整| 一个人看的www免费观看视频| 日日撸夜夜添| 久久人人爽人人片av| 成年版毛片免费区| 欧美又色又爽又黄视频| 亚洲欧美日韩卡通动漫| 亚洲av成人精品一区久久| 精品人妻偷拍中文字幕| 搡老妇女老女人老熟妇| 日韩欧美国产在线观看| 亚洲精品乱码久久久久久按摩| 18禁在线无遮挡免费观看视频| 看黄色毛片网站| 亚洲欧洲国产日韩| 亚洲欧美清纯卡通| 久久久久久久久中文| 九九久久精品国产亚洲av麻豆| 国产成人aa在线观看| 欧美zozozo另类| 中国美女看黄片| 亚洲欧美成人精品一区二区| 久久久久久久久大av| 一级黄片播放器| 国产午夜福利久久久久久| 两个人视频免费观看高清| 日本一本二区三区精品| 97人妻精品一区二区三区麻豆| 99在线视频只有这里精品首页| 麻豆成人av视频| 99久久精品热视频| 人妻系列 视频| 成年女人看的毛片在线观看| 中文精品一卡2卡3卡4更新| 女的被弄到高潮叫床怎么办| 全区人妻精品视频| 一个人免费在线观看电影| 欧美另类亚洲清纯唯美| 人体艺术视频欧美日本| 日本爱情动作片www.在线观看| 亚洲av中文字字幕乱码综合| 亚洲成人精品中文字幕电影| 国产片特级美女逼逼视频| 成人亚洲欧美一区二区av| 日韩成人av中文字幕在线观看| 欧美最黄视频在线播放免费| 男女下面进入的视频免费午夜| 日韩 亚洲 欧美在线| 日韩一本色道免费dvd| 免费人成视频x8x8入口观看| 中国美白少妇内射xxxbb| 久久人人爽人人爽人人片va| 国产精品久久久久久久久免| 国产精品久久久久久亚洲av鲁大| 亚洲av免费高清在线观看| 天天一区二区日本电影三级| 中文精品一卡2卡3卡4更新| 午夜福利成人在线免费观看| 久久久久久九九精品二区国产| 少妇猛男粗大的猛烈进出视频 | 22中文网久久字幕| 亚洲国产精品成人综合色| 中文资源天堂在线| 国产成年人精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 久久人人精品亚洲av| 欧美性猛交黑人性爽| 真实男女啪啪啪动态图| 国产精品av视频在线免费观看| 日韩 亚洲 欧美在线| 国产黄片美女视频| 亚洲av第一区精品v没综合| 欧美变态另类bdsm刘玥| 国产精品国产三级国产av玫瑰| 亚洲人成网站在线播| 久久99精品国语久久久| 最近手机中文字幕大全| 国产毛片a区久久久久| 一级毛片电影观看 | 国产亚洲欧美98| 国产午夜福利久久久久久| av视频在线观看入口| 啦啦啦观看免费观看视频高清| 一进一出抽搐动态| 午夜视频国产福利| 久久国产乱子免费精品| av卡一久久| 亚洲在线自拍视频| 国产精品一区二区性色av| 久久6这里有精品| 久久久国产成人免费| 午夜激情福利司机影院| 久久精品国产鲁丝片午夜精品| 蜜臀久久99精品久久宅男| av在线播放精品| 大又大粗又爽又黄少妇毛片口| 国产不卡一卡二| 国产精品一区www在线观看| 日本成人三级电影网站| 搡老妇女老女人老熟妇| av天堂在线播放| 国产精品伦人一区二区| 国产极品天堂在线| 免费搜索国产男女视频| 麻豆一二三区av精品| 九九久久精品国产亚洲av麻豆| 狂野欧美激情性xxxx在线观看| 日本av手机在线免费观看| 国产午夜精品一二区理论片| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 1000部很黄的大片| videossex国产| 国产白丝娇喘喷水9色精品| 岛国在线免费视频观看| 91精品一卡2卡3卡4卡| 日韩三级伦理在线观看| 日本三级黄在线观看| 日韩中字成人| 日韩三级伦理在线观看| 97超视频在线观看视频| 精品欧美国产一区二区三| 看片在线看免费视频| 日韩欧美三级三区| 亚洲国产精品国产精品| 国产不卡一卡二| 91精品国产九色| 日本与韩国留学比较| 1024手机看黄色片| 欧美区成人在线视频| 97在线视频观看| 蜜桃亚洲精品一区二区三区| 亚洲成人久久性| 一本久久中文字幕| 欧美不卡视频在线免费观看| 日韩强制内射视频| 免费观看精品视频网站| 亚洲国产欧美在线一区| 男女那种视频在线观看| 少妇被粗大猛烈的视频| 日本黄色片子视频| 日韩精品青青久久久久久| 国产精品一区www在线观看| 99久久久亚洲精品蜜臀av| 好男人在线观看高清免费视频| 中文亚洲av片在线观看爽| 超碰av人人做人人爽久久| 欧美高清成人免费视频www| 国产亚洲91精品色在线| 人妻系列 视频| 亚洲精品日韩在线中文字幕 | 嫩草影院入口| 日日摸夜夜添夜夜添av毛片| 有码 亚洲区| 成人欧美大片| 亚洲自偷自拍三级| 成人美女网站在线观看视频| 亚洲综合色惰| 中文字幕制服av| 91av网一区二区| 亚洲成a人片在线一区二区| 亚洲成人av在线免费| av又黄又爽大尺度在线免费看 | 又粗又硬又长又爽又黄的视频 | 成人欧美大片| 日产精品乱码卡一卡2卡三| 亚洲成人久久爱视频| 你懂的网址亚洲精品在线观看 | 禁无遮挡网站| 九九在线视频观看精品| 天天躁日日操中文字幕| .国产精品久久| 精品久久久久久久末码| 亚洲经典国产精华液单| 国产成人a∨麻豆精品| 免费av不卡在线播放| 99久久无色码亚洲精品果冻| 欧美激情国产日韩精品一区| 深爱激情五月婷婷| 在线播放无遮挡| 精品久久国产蜜桃| 国产成人精品一,二区 | 99久久人妻综合| 精华霜和精华液先用哪个| 欧美丝袜亚洲另类| 日日啪夜夜撸| 国产中年淑女户外野战色| 亚洲中文字幕一区二区三区有码在线看| 亚洲婷婷狠狠爱综合网| 亚洲最大成人中文| 午夜亚洲福利在线播放| 久久久久久久久久成人| 99九九线精品视频在线观看视频| 亚洲精华国产精华液的使用体验 | 熟女电影av网| 51国产日韩欧美| 婷婷色综合大香蕉| 国产日韩欧美在线精品| 超碰av人人做人人爽久久| 欧美最新免费一区二区三区| 国产精品久久久久久精品电影小说 | 三级国产精品欧美在线观看| 久久精品夜色国产| 精品国产三级普通话版| 啦啦啦观看免费观看视频高清| 美女xxoo啪啪120秒动态图| 国语自产精品视频在线第100页| 国内精品宾馆在线| 日韩欧美精品v在线| 99久久成人亚洲精品观看| 国产一区二区激情短视频| 搞女人的毛片| 免费观看的影片在线观看| av天堂在线播放| 亚洲图色成人| 成人特级av手机在线观看| 亚洲欧美日韩东京热| 高清毛片免费观看视频网站| 亚洲av中文av极速乱| 中文字幕制服av| 国产男人的电影天堂91| 久久久久免费精品人妻一区二区| 性色avwww在线观看| 男女啪啪激烈高潮av片| 精品99又大又爽又粗少妇毛片| 精品欧美国产一区二区三| 三级国产精品欧美在线观看| 毛片一级片免费看久久久久| 国产成人影院久久av| av天堂中文字幕网| 嫩草影院入口| 熟女电影av网| 在线免费十八禁| 22中文网久久字幕| 99热全是精品| 国产精品麻豆人妻色哟哟久久 | 女人被狂操c到高潮| 亚洲自拍偷在线| 亚洲一区高清亚洲精品| 高清日韩中文字幕在线| 久久久欧美国产精品| 久久午夜福利片| 两个人的视频大全免费| 特大巨黑吊av在线直播| 久久中文看片网| 亚洲国产日韩欧美精品在线观看| 高清毛片免费观看视频网站| 国产午夜精品久久久久久一区二区三区| 国产亚洲av片在线观看秒播厂 | 18禁黄网站禁片免费观看直播| 丝袜喷水一区| 深爱激情五月婷婷| 天堂中文最新版在线下载 | 蜜臀久久99精品久久宅男| 级片在线观看| 国产精品女同一区二区软件| 国产午夜福利久久久久久| 免费看日本二区| 又粗又硬又长又爽又黄的视频 | 久久综合国产亚洲精品| 永久网站在线| 嘟嘟电影网在线观看| 国产一区亚洲一区在线观看| 成人高潮视频无遮挡免费网站| 精品一区二区三区视频在线| 久久精品夜夜夜夜夜久久蜜豆| 午夜视频国产福利| 能在线免费看毛片的网站| 人人妻人人澡欧美一区二区| 十八禁国产超污无遮挡网站| 国产黄片视频在线免费观看| 久久久久久九九精品二区国产| 亚洲精品乱码久久久v下载方式| 淫秽高清视频在线观看| 在线天堂最新版资源| 国产精品美女特级片免费视频播放器| 国产精品福利在线免费观看| 亚洲国产精品久久男人天堂| 爱豆传媒免费全集在线观看| 亚洲精品影视一区二区三区av| 久久人人爽人人片av| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲网站| 亚洲精品自拍成人| 国产精品伦人一区二区| 久久精品国产亚洲av香蕉五月| 国产乱人偷精品视频| 国产高清视频在线观看网站| 成人一区二区视频在线观看| 丰满乱子伦码专区| 国产成人午夜福利电影在线观看| 久久6这里有精品| 国产乱人视频| 老师上课跳d突然被开到最大视频| 一级二级三级毛片免费看| 国产一级毛片七仙女欲春2| 校园春色视频在线观看| 久久久精品欧美日韩精品| 亚洲电影在线观看av| 成人亚洲精品av一区二区| 国产精品爽爽va在线观看网站| 成人特级黄色片久久久久久久| 99久国产av精品| 久久精品综合一区二区三区| 在线a可以看的网站| 简卡轻食公司| 精品一区二区免费观看| 特级一级黄色大片| 久99久视频精品免费| 精品一区二区免费观看| 国产av麻豆久久久久久久| 精品免费久久久久久久清纯| 国产精品久久久久久久电影| 久久草成人影院| 欧美性猛交黑人性爽| 亚洲色图av天堂| 在线免费十八禁| 97人妻精品一区二区三区麻豆| 亚洲三级黄色毛片| 一区二区三区免费毛片| 亚洲自拍偷在线| 亚洲av中文字字幕乱码综合| 国产精华一区二区三区| 久久精品久久久久久噜噜老黄 | av天堂在线播放|