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

    IMERG系列產(chǎn)品在大汶河流域的降水監(jiān)測能力評價

    2023-08-08 14:59:36段震王超邸燕劉強王剛
    人民長江 2023年7期
    關(guān)鍵詞:大汶河降水流域

    段震 王超 邸燕 劉強 王剛

    摘要:

    為探究高分辨率多衛(wèi)星聯(lián)合反演IMERG(Integrated Multi-satellitE Retrievals for GPM)衛(wèi)星降水產(chǎn)品在大汶河流域的監(jiān)測精度與適用性,以2001~2019年較長系列地面雨量站實測數(shù)據(jù)為基準,采用探測率、虛報率等11個指標,對IMERG系列Early、Late、Final產(chǎn)品在研究區(qū)的降水監(jiān)測能力進行評價。結(jié)果表明:① 日、月、年尺度下,F(xiàn)inal產(chǎn)品的精度最好,Early、Late兩個產(chǎn)品的精度基本一致,在時間上,IMERG系列產(chǎn)品均能較好反映流域降水的年際變化與年內(nèi)分配;在空間上,表現(xiàn)為整體高估研究區(qū)降水量,低估泰山區(qū)域降水量(相對偏差小于-20%)。② 在不同強度日降水探測能力方面,IMERG系列產(chǎn)品間差異不大,僅對小雨有較好的命中率(探測率>0.6)。各產(chǎn)品探測能力一般,整體表現(xiàn)為誤報率較高(誤報率>0.7)、探測率與關(guān)鍵成功率較低,同時,隨著降水強度的增大,各產(chǎn)品探測能力呈下降趨勢。③ IMERG系列產(chǎn)品均具有一定的極端降水監(jiān)測能力,F(xiàn)inal產(chǎn)品評價結(jié)果更優(yōu)、表現(xiàn)更為穩(wěn)定,在參與評價的5個指數(shù)中,F(xiàn)inal產(chǎn)品有4個指數(shù)結(jié)果最優(yōu),Late產(chǎn)品有1個最優(yōu),Early產(chǎn)品則有4個指數(shù)結(jié)果優(yōu)于Late產(chǎn)品??傮w上,IMERG系列數(shù)據(jù)均能較好反映大汶河流域(除泰山區(qū)域)降水的時空分布與極端降水情況,F(xiàn)inal產(chǎn)品整體精度最高,Early、Late產(chǎn)品精度較為接近,其中,Early產(chǎn)品在極端降水監(jiān)測中的表現(xiàn)好于Late產(chǎn)品,而對不同強度日降水事件的監(jiān)測能力差于Late產(chǎn)品。

    關(guān) 鍵 詞:

    衛(wèi)星遙感降水反演; IMERG; 精度評估; 極端降水; 大汶河流域

    中圖法分類號: P426

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2023.07.013

    0 引 言

    降水是水文循環(huán)的關(guān)鍵環(huán)節(jié),是地球上各類水體的直接或間接補給源,準確、可靠的降水監(jiān)測對于區(qū)域水資源評價與管理、洪水預(yù)報、極端降水研究等工作具有重要意義[1]。目前,主要的降水觀測手段有雨量站監(jiān)測、地面雷達監(jiān)測與衛(wèi)星傳感器監(jiān)測[2]。其中,雨量站與地基雷達監(jiān)測均依賴于地面布設(shè)的儀器設(shè)備,雖然整體監(jiān)測精度較高,但受維護成本、布設(shè)密度、地形等因素的影響,往往難以準確反映較大尺度降水的空間分布情況[3-4]。相較于地面觀測,衛(wèi)星降水監(jiān)測受人類活動影響更小、覆蓋范圍更廣、數(shù)據(jù)連續(xù)性更強,隨著遙感技術(shù)與反演算法的不斷完善,一系列衛(wèi)星遙感降水產(chǎn)品為氣候變化、水文研究等工作提供了新的數(shù)據(jù)支撐[5]。

    全球降水觀測(Global Precipitation Measurement,GPM)計劃是目前全球最為先進的衛(wèi)星降水觀測計劃,其三級降水產(chǎn)品IMERG 的發(fā)布,則被認為是衛(wèi)星遙感降水反演技術(shù)進入新時代的標志[6-9]。自IMERG產(chǎn)品2014年3月正式發(fā)布以來,國內(nèi)外學(xué)者在全球眾多國家、地區(qū)對IMERG產(chǎn)品精度進行了檢驗,并與其他主流衛(wèi)星降水產(chǎn)品進行對比,證明了該產(chǎn)品在降水監(jiān)測方面具有一定的能力[10-12]。隨著V06B算法版本的發(fā)布,IMERG產(chǎn)品的時間序列延長到了2000年6月,使得該產(chǎn)品在氣象研究、水文應(yīng)用中具備了更大的潛力。

    近年來,國內(nèi)對于該產(chǎn)品的研究逐漸增多??子睿?3]、張茹等[14]分析發(fā)現(xiàn)IMERG產(chǎn)品能較好反映中國大陸區(qū)域的降水分布情況;劉若蘭[15]、Fang[16]等則證明該產(chǎn)品對極端降水事件有一定的監(jiān)測能力,可以較好地擬合其年際間的變化趨勢;但不同區(qū)域的研究成果[17-19]表明,IMERG產(chǎn)品受地形、天氣等因素的影響,存在著區(qū)域差異性與不確定性等特點,因此在應(yīng)用時需要進行精度檢驗,必要時還需對產(chǎn)品進行校正。雖然目前對于IMERG產(chǎn)品的研究已取得了豐碩成果,但現(xiàn)有研究仍存在一定不足。首先,IMERG依據(jù)算法的不同分為了Early、Late、Final共3套產(chǎn)品,因Final產(chǎn)品經(jīng)過GPCC(Global Precipitation Climatology Center)月數(shù)據(jù)校正,被GPM推薦為研究首選,所以現(xiàn)有的研究中主要采用Final數(shù)據(jù),然而,該數(shù)據(jù)從監(jiān)測到發(fā)布時間超過3個月,不能滿足水文預(yù)報等研究的時效性要求,因此,對發(fā)布時間不同的3套產(chǎn)品適用情況進行對比檢驗,評價各產(chǎn)品的精度情況,對于不同時間需求的水文研究具有重要意義;其次,目前的研究中,采用數(shù)據(jù)的時間系列普遍較短[20-21],大多小于5 a,檢驗樣本偏少,難以反映IMERG產(chǎn)品真實的降水監(jiān)測能力。

    大汶河是黃河下游最大一級支流,是保障和實現(xiàn)黃河流域生態(tài)保護和高質(zhì)量發(fā)展的重要組成區(qū)域。其流域內(nèi)地形起伏跌宕,既有海拔高達1 545 m連綿的中低山(泰山)山脈,又有平均海拔約為110 m較為平坦的山前平原,相對高差達1 400余米。受地形、人類活動等因素的影響,流域內(nèi)水文站(網(wǎng))有限,且分布不均。近幾年受極端氣候影響,流域內(nèi)小概率的洪水事件頻發(fā),對流域內(nèi)水文信息采集與水文預(yù)報工作提出了更高的要求。本文基于流域內(nèi)2001~2019年19 a較長系列的實測水文資料,選取探測率、虛報率、關(guān)鍵成功率等11個參數(shù)指標,從不同方面定量評價了IMERG系列不同產(chǎn)品(Early、Late、Final)在該流域內(nèi)的不同強度與不同時空條件下的降水監(jiān)測能力。

    1 研究區(qū)域與方法

    1.1 研究區(qū)概況

    大汶河位于山東省中西部,地理坐標為35°59′N~36°7′N、116°12′E~117°56′E,自東向西流經(jīng)濟南、泰安、濟寧的8個縣市區(qū),于東平縣馬口村入東平湖,再由東平湖入黃河。流域?qū)偃A北暖溫帶半濕潤季風(fēng)性大陸氣候區(qū),多年平均降水量708.9 mm,其中約70%的降水發(fā)生在汛期。流域地形東寬西窄、地勢東高西低,上游為魯中山區(qū),中下游為丘陵平原,山區(qū)與丘陵約占總流域面積的68%。受降水、地形等因素的影響,流域內(nèi)水旱災(zāi)害頻發(fā),由于雨量站主要布設(shè)于人類活動的城市、農(nóng)村等區(qū)域,站點分布不均,局部實測數(shù)據(jù)往往難以準確反映區(qū)域空間降水情況。綜合考慮流域的地形、降水、雨量站分布等情況,本次選取大汶河戴村壩以上流域為研究區(qū)(見圖1),區(qū)域面積8 264 km2,約占總流域的90%。

    1.2 數(shù)據(jù)來源

    本次衛(wèi)星遙感降水?dāng)?shù)據(jù)采用IMERG系列提供的Early、Late兩個近實時的降水產(chǎn)品以及經(jīng)GPCC月數(shù)據(jù)校正的Final降水產(chǎn)品。3套產(chǎn)品發(fā)布延遲時間分別為4 h,12 h和3.5個月。產(chǎn)品數(shù)據(jù)下載自美國宇航局GPM計劃網(wǎng)站(https:∥gpm.nasa.gov/data/),數(shù)據(jù)算法版本采用最新的V06B,系列為2001~2019年,時間分辨率為1 d,空間分辨率為0.1°×0.1°,變量為PrecipitationCal(經(jīng)校準的降水?dāng)?shù)據(jù))。

    地面資料采用流域內(nèi)18個雨量站2001~2019年的實測日降水?dāng)?shù)據(jù),站點位置分布情況見圖1。

    1.3 研究方法

    本次研究以地面雨量站觀測的降水?dāng)?shù)據(jù)為基準,選用11個參數(shù)指標(見表1)對IMERG系列3套產(chǎn)品在大汶河流域的降水監(jiān)測能力進行分析評估。具體評估內(nèi)容分為以下3部分:① 采用皮爾遜相關(guān)系數(shù)(Correlation Coefficient,CC)、相對偏差(Relative Bias,RB)、均方根誤差(Root Mean-Squared Error,RMSE)對降水產(chǎn)品日、月、年尺度下的精度進行定量評價,其中CC用于評價 IMERG 產(chǎn)品與地面數(shù)據(jù)在時間序列上的相關(guān)程度,RB用于衡量 IMERG 產(chǎn)品較地面數(shù)據(jù)的偏離程度,負值為低估,正值為高估,RMSE用于估計IMERG產(chǎn)品與地面數(shù)據(jù)在時間序列上的平均誤差幅度;② 采用探測率(Probability of Detection,POD)、虛報率(False Alarm Ratio,F(xiàn)AR)、關(guān)鍵成功率(Critical Success-Index,CSI)[22]對IMERG產(chǎn)品的不同強度日降水探測能力進行評價;③ 采用持續(xù)干燥日數(shù)CDD、持續(xù)濕潤日數(shù)CWD、最大1 d降水量Rx1day、最大5 d降水量Rx5day、強降水量R95p[23]共5個指標對IMERG產(chǎn)品的極端降水監(jiān)測能力進行評估。本文采用泰森多邊形法計算站點面雨量,采用柵格面積加權(quán)計算IMERG產(chǎn)品面雨量。

    2 結(jié)果與分析

    2.1 不同時空尺度的精度評價

    2.1.1 日尺度空間精度評價

    圖2給出了IMERG 3套產(chǎn)品數(shù)據(jù)日尺度下的空間精度分布情況,從該圖可以看出:

    (1) 3套產(chǎn)品與地面站點數(shù)據(jù)在時間序列上的相關(guān)程度較好,Early產(chǎn)品的CC值為0.60~0.76,Late產(chǎn)品為0.64~0.78,F(xiàn)inal產(chǎn)品為0.70~0.80。從空間上看,Early、Late產(chǎn)品的CC值在流域各區(qū)域的差異較大,F(xiàn)inal產(chǎn)品的表現(xiàn)則更為穩(wěn)定,同時各產(chǎn)品均在流域東南部有最高的相關(guān)程度,在中北部泰山區(qū)域相關(guān)程度最低。

    (2) 從RMSE值的空間分布情況看,3套產(chǎn)品在泰山區(qū)域的誤差幅度明顯高于其余區(qū)域,泰山站處的RMSE值為8.4~9.5 mm,其余區(qū)域大多為5.8~7.0 mm;

    (3) 從RB值的空間分布情況看,3套產(chǎn)品均低估了泰山區(qū)域的降水情況,其中Final產(chǎn)品低估最為明顯,RB值達到了-30.1%,在其他區(qū)域,各產(chǎn)品又表現(xiàn)出明顯高估降水的情況,Early與Late兩個產(chǎn)品在大多數(shù)區(qū)域的RB值超過了15%,而Final產(chǎn)品約2/3區(qū)域RB值介于0~10%之間,表現(xiàn)較好。

    綜合來看,F(xiàn)inal產(chǎn)品的日尺度精度最好,Early與Late產(chǎn)品表現(xiàn)較為接近,3套產(chǎn)品均整體上高估了(除泰山區(qū)域外)地面的降水量,且在泰山區(qū)域形成一個明顯的降水精度低值區(qū)。該區(qū)域的精度主要是由參與評價的泰山站反映,由于該站點位于泰山山頂(海拔大于1 500 m),有獨特的小氣候,是流域內(nèi)的一個海拔高點與降水量極大值點,而與之對應(yīng)的IMERG柵格覆蓋范圍內(nèi)既有泰山主峰又有海拔較低的城區(qū),最大高差超過了1 300 m,因此,以IMERG柵格的面雨量作為站點對應(yīng)的點雨量值進行評價,可能是造成該處精度較低的原因之一。

    2.1.2 月尺度面雨量精度評價

    圖3給出了IMERG 3套產(chǎn)品各月份數(shù)據(jù)的定量指標結(jié)果,從該圖可以看出,F(xiàn)inal產(chǎn)品與地面站點數(shù)據(jù)的相關(guān)性最好,各月份的CC值均大于0.90,Early、Late產(chǎn)品CC值的年內(nèi)差異較大,最差結(jié)果均出現(xiàn)在2月,分別為0.17和0.23,最好結(jié)果均出現(xiàn)在11月份,分別為 0.94和0.93;RMSE值的大小主要受月降水量的影響,月降水量越大均方根誤差越大,Early、Late 兩個產(chǎn)品的RMSE明顯大于Final產(chǎn)品,2個產(chǎn)品各有7個月的RMSE超過了20 mm,F(xiàn)inal 產(chǎn)品僅有2個月超過20 mm且有7個月未超過10 mm;從相對偏差結(jié)果來看,3套產(chǎn)品均有11個月的結(jié)果大于0,同時,Early產(chǎn)品在6~8月的RB為-6.6%~3.6%,其余月份為38.0%~170.3%,Late 產(chǎn)品在6~8月的RB為-1.9%~5.2%,其余月份為25.7%~138.1%,F(xiàn)inal產(chǎn)品在6~9月的RB為-1.2%~4.5%,其余月份為17.3%~60.0%,表明3套產(chǎn)品汛期(山東省汛期為6~9月)觀測結(jié)果較地面數(shù)據(jù)的偏離程度好于其余月份。綜合圖3結(jié)果,F(xiàn)inal產(chǎn)品在各月份的精度均高于其余2個產(chǎn)品,Early、Late產(chǎn)品在4,6,7,8,9,11月中有較好的表現(xiàn),特別是7,8月份。因此,F(xiàn)inal產(chǎn)品應(yīng)作為IMERG系列資料應(yīng)用的首選,而對于時效性要求較高的研究,若研究期為上述月份時可根據(jù)需求選用Early、Late數(shù)據(jù),若為其他月份則應(yīng)謹慎選用。

    圖4給出了IMERG 3套產(chǎn)品以及地面站點數(shù)據(jù)的多年平均月降水量結(jié)果。從該圖可以看出,3套產(chǎn)品均能較好地反映研究區(qū)年內(nèi)降水的分布情況,全年各月中,Early產(chǎn)品在4,8兩個月份的數(shù)據(jù)最接近地面數(shù)據(jù),Late產(chǎn)品在6月份的數(shù)據(jù)最接近地面數(shù)據(jù),F(xiàn)inal產(chǎn)品在剩余的9個月份中表現(xiàn)最好。

    2.1.3 年尺度面雨量精度評價

    圖5給出了IMERG 3套產(chǎn)品以及地面站點數(shù)據(jù)的年降水量年際變化過程。從該圖可以看出,IMERG 3套降水產(chǎn)品均能較好地反映流域年降水量的年際變化情況,其中,F(xiàn)inal產(chǎn)品與站點資料更為接近,Early、Late產(chǎn)品表現(xiàn)較為一致。

    為更好地分析降水產(chǎn)品的空間精度,繪制了流域多年平均年降水量空間分布圖(見圖6)。結(jié)果表明,3套產(chǎn)品雖然模擬出了大汶河流域年降水量東高西低的趨勢,但均未正確模擬出在泰山處的降水極大值區(qū)域,對該區(qū)域的降水量存在明顯低估,而在其他區(qū)域則表現(xiàn)出高估地面降水的情況,這與日尺度下的檢驗結(jié)果基本一致。

    2.2 不同降水強度的探測能力評價

    2.2.1 不同降水量等級評價

    圖7為IMERG 3套產(chǎn)品在不同日降水等級(等級劃分見表2)下POD、FAR、CSI指標的箱線圖。從圖7(a)可以看出,3套產(chǎn)品對小雨事件的探測率整體在0.6以上,顯著高于其他降水事件,這表明IMERG產(chǎn)品對小雨等級日降水事件命中率更高;圖7(b)、圖7(c)中,F(xiàn)AR值主要分布在0.6~0.8之間,CSI值主要分布在0.1~0.3之間,說明3套產(chǎn)品對各降水事件都存在較高的誤報率與較低的關(guān)鍵成功率。從圖中箱體及范圍的變化情況可以看出,隨著日降水等級的提高,指標箱體變長,范圍變大,表明IMERG產(chǎn)品在各站點處的精度差異逐漸加大。從整體結(jié)果看,IMERG產(chǎn)品對于準確監(jiān)測研究區(qū)日降水等級的能力較差,各產(chǎn)品間的探測效果接近,其中,Late產(chǎn)品對于大暴雨事件整體監(jiān)測效果最好,F(xiàn)inal產(chǎn)品在其余降水事件中效果最好。

    圖8給出了各降水?dāng)?shù)據(jù)的不同強度面降水事件發(fā)生頻率。從該圖可以看出,IMERG數(shù)據(jù)明顯低估無雨事件的發(fā)生頻率,3套產(chǎn)品較地面數(shù)據(jù)分別低22.25%,18.06%,16.28%;明顯高估小雨事件發(fā)生頻率,各產(chǎn)品分別高21.43%,17.15%,16.10%,即IMERG產(chǎn)品在大量無雨(日降水<0.1 mm)日內(nèi)監(jiān)測到了0.1~9.9 mm的降水量。結(jié)果表明,IMERG產(chǎn)品對小于0.1 mm的微量降水存在明顯高估的情況。而其他強度降水事件發(fā)生頻率較低,各數(shù)據(jù)間差異不大。

    2.2.2 不同降水量閾值評價

    通過統(tǒng)計分析流域內(nèi)各站點的日降水事件發(fā)現(xiàn),流域內(nèi)發(fā)生大雨以下事件的頻率達到了90%以上。為了更好地分析IMERG產(chǎn)品對日降水的探測能力,本文分別以0.1,1,5,10,20,30,40,50,100 mm為閾值,對IMERG產(chǎn)品在各雨量站點處的POD、FAR、CSI指標進行計算。

    圖9給出了IMERG 3套產(chǎn)品在不同降水量閾值下POD、FAR、CSI指標的箱線圖。從該圖可以看出,3套產(chǎn)品的表現(xiàn)基本一致,在探測能力上,隨著閾值的增大,IMERG產(chǎn)品的POD呈現(xiàn)下降趨勢;從誤報率上看,各產(chǎn)品整體在0.3以上,隨著閾值的增大,F(xiàn)AR呈現(xiàn)出先減小后增大的趨勢,當(dāng)閾值為5 mm時,F(xiàn)AR值最??;從關(guān)鍵成功率上看,各產(chǎn)品整體在0.6以下,隨著閾值的增大,CSI呈現(xiàn)出先增大后減小的趨勢,當(dāng)閾值為1 mm和5 mm時,CSI表現(xiàn)最好。從FAR、CSI兩個指標可以看出,當(dāng)閾值由0.1 mm增大為1 mm后,IMERG產(chǎn)品的監(jiān)測精度得到了明顯改善,這可能是由于GPM衛(wèi)星傳感器對于弱降水的監(jiān)測更為敏感,以及部分月份中地面雨量計的分辨率為0.5 mm不能滿足0.1 mm的閾值要求引起的。綜合3個指標成果,IMERG產(chǎn)品對于0.1,1,5,10 mm等4個閾值強度的降水事件探測能力更強,當(dāng)閾值為100 mm時,F(xiàn)inal產(chǎn)品的探測能力最差,而在其他閾值下Final綜合監(jiān)測效果最好。綜合圖7與圖9指標結(jié)果來看,當(dāng)降水等級不超過暴雨時,F(xiàn)inal產(chǎn)品仍是研究首選,其次為Late產(chǎn)品;對于大暴雨事件,Late產(chǎn)品的監(jiān)測效果最好,其次為Early產(chǎn)品。

    2.3 極端降水監(jiān)測能力評估

    圖10給出了各降水?dāng)?shù)據(jù)極端降水指數(shù)的年際變化過程。由該圖可知,IMERG 3套產(chǎn)品基本能反映5個極端降水指數(shù)的年際變化趨勢,但部分年份的指數(shù)誤差較大,如圖10(a)中,基于地面數(shù)據(jù)計算的2003年持續(xù)濕潤日數(shù)為5 d,IMERG產(chǎn)品的數(shù)值均為13 d,該年份誤差超過了100%。

    表3給出了各指數(shù)多年平均成果,表4給出了IMERG 3套產(chǎn)品極端降水指數(shù)的精度結(jié)果。結(jié)合表3、表4具體指標結(jié)果來看:① 大汶河流域多年平均持續(xù)濕潤日數(shù)為6.5 d,Late產(chǎn)品在該指標上表現(xiàn)最好,多年平均CWD為6.5,CC值為0.48,RB值僅為-0.8%。② 流域多年平均持續(xù)干燥日數(shù)為43.2 d,Early、Late兩種產(chǎn)品的結(jié)果明顯偏小,RB值分別為-20.6%與-22.9%,F(xiàn)inal產(chǎn)品表現(xiàn)稍好,多年平均CDD值為41.5,CC值為0.53,但同樣低估了流域持續(xù)干燥日數(shù);IMERG產(chǎn)品整體上高估CWD、低估CDD的情況可能是由于IMERG產(chǎn)品對于弱降水的高估效應(yīng)引起的。③ 對于R95p指數(shù),3套產(chǎn)品均有較好的CC值成果(CC>0.8),F(xiàn)inal產(chǎn)品的RB值為-1.5%,成果最優(yōu),Late產(chǎn)品RB值為11.7%,結(jié)果最差。④ 流域多年平均最大1 d降水量為71.0 mm,該指標下各產(chǎn)品的表現(xiàn)與R95p相似,F(xiàn)inal產(chǎn)品精度最高,Rx1day值為70.3。⑤ 對于Rx5day指數(shù),雖然Early產(chǎn)品的誤差?。≧B=-0.3%),但相關(guān)系數(shù)(0.47)明顯小于Final 產(chǎn)品值(0.78),綜合來看,F(xiàn)inal產(chǎn)品對于流域最大5 d降水量的監(jiān)測能力更好。

    由上述分析可知,IMERG產(chǎn)品對大汶河流域極端降水具有一定的監(jiān)測能力,其中,F(xiàn)inal產(chǎn)品的綜合監(jiān)測能力最強,其他2種產(chǎn)品對不同指標的表現(xiàn)各有優(yōu)劣。具體來看,Late產(chǎn)品對CWD的監(jiān)測效果最好,F(xiàn)inal 產(chǎn)品對CDD、R95p、Rx1day、Rx5day共4個指標的監(jiān)測效果最好,Early產(chǎn)品對除CWD外的其余4個指標整體監(jiān)測效果優(yōu)于Late產(chǎn)品。

    3 討論與結(jié)論

    3.1 討 論

    衛(wèi)星降水產(chǎn)品的發(fā)布,是由美國宇航局GPM計劃網(wǎng)站按照同一計算方法與標準發(fā)布的。其不同產(chǎn)品在不同區(qū)域、流域中的應(yīng)用受其區(qū)域、流域下墊面條件、地形氣候等不同因素的影響,在對某一水文現(xiàn)象特征的識別與響應(yīng)中,既有宏觀上的“統(tǒng)一性”特征,又有其區(qū)域、流域上的“特殊”性特征。從日尺度下空間精度分析來看,IMERG各產(chǎn)品對流域東部山丘區(qū)的降水監(jiān)測精度好于西部平原,這與李芳等[24]在大汶河相鄰的魯南區(qū)域研究結(jié)果是相近的。同時,各產(chǎn)品在空間上的精度變化趨勢基本一致,并在流域地形急劇變化且海拔最高的泰山山頂區(qū)域出現(xiàn)了極差的檢驗結(jié)果,這與王思夢等[25]得出的地形影響產(chǎn)品精度的結(jié)論是一致的,但其“隨海拔的升高,相關(guān)系數(shù)呈現(xiàn)出增加-減少-增加”的結(jié)論又和本文不同,這可能是由流域大小、位置、氣候以及地形等因素引起的,說明IMERG產(chǎn)品精度存在一定的區(qū)域性差異。因此,在不同區(qū)域、流域中應(yīng)用時,需要事先從各方面對其適用性進行對比分析和探討。

    3.2 結(jié) 論

    本文以地面雨量站資料為基準,采用11個參數(shù)指標,從數(shù)據(jù)精度定量評價、不同強度日降水探測能力以及極端降水監(jiān)測能力共3個方面對IMERG系列3套產(chǎn)品在大汶河流域的降水監(jiān)測精度進行評價,主要結(jié)論如下。

    (1) IMERG系列產(chǎn)品均能較好反映流域內(nèi)降水的空間分布、年內(nèi)分配與年際變化情況。各產(chǎn)品中,F(xiàn)inal數(shù)據(jù)整體精度最高,Early、Late表現(xiàn)基本一致;空間上,各產(chǎn)品均存在高估流域(除泰山區(qū)域)、低估泰山區(qū)域降水量的情況;在時間上,F(xiàn)inal產(chǎn)品對各月降水量均有較好的監(jiān)測效果,Early、Late產(chǎn)品各月份差異較大,在4,6,7,8,9,11月中有較高的精度。

    (2) IMERG系列產(chǎn)品對于不同強度日降水的探測能力較為接近且表現(xiàn)一般,均存在較高的誤報率與較低的探測率、關(guān)鍵成功率。隨著日降水強度的增大,各產(chǎn)品的探測能力呈下降趨勢,同時各產(chǎn)品在流域不同區(qū)域的精度差異也逐漸增大。降水強度不超過暴雨時,產(chǎn)品精度為Final>Late>Early;降水強度為大暴雨時,產(chǎn)品精度為Late>Early>Final。

    (3) IMERG系列產(chǎn)品在大汶河流域具有一定的極端降水監(jiān)測能力,從選取的5個極端降水指數(shù)精度結(jié)果上看,各產(chǎn)品的監(jiān)測能力為Final>Early>Late,特別是在R95p、Rx1day、Rx5day等3個指數(shù)上,Early產(chǎn)品表現(xiàn)出與Final產(chǎn)品較為接近的精度結(jié)果。

    綜合上述結(jié)論,F(xiàn)inal數(shù)據(jù)應(yīng)作為IMERG系列產(chǎn)品在大汶河流域應(yīng)用的首選,若所需的Final數(shù)據(jù)尚未發(fā)布,則可根據(jù)以下情況選用Early、Late數(shù)據(jù):① 當(dāng)研究期為4,6,7,8,9,11月時,可根據(jù)研究的時效性需求選用Early、Late數(shù)據(jù),其他月份應(yīng)謹慎采用;② 當(dāng)研究中涉及不同等級日降水時,則應(yīng)先選Late數(shù)據(jù);③ 當(dāng)研究涉及極端降水時,則應(yīng)優(yōu)先選擇Early數(shù)據(jù)。若研究區(qū)域為泰山區(qū)域時,IMERG系列3套產(chǎn)品均應(yīng)謹慎選擇。

    參考文獻:

    [1] 唐國強,萬瑋,曾子悅,等.全球降水測量(GPM)計劃及其最新進展綜述[J].遙感技術(shù)與應(yīng)用,2015,30(4):607-615.

    [2] 唐國強.衛(wèi)星遙感降水在全球及典型區(qū)域的檢驗、應(yīng)用和改進[D].北京:清華大學(xué),2019.

    [3] 劉俊峰,陳仁升,韓春壇,等.多衛(wèi)星遙感降水?dāng)?shù)據(jù)精度評價[J].水科學(xué)進展,2010,21(3),343-348.

    [4] KIDD C,HUFFMAN G.Global precipitation measurement[J].Meteorological Applications,2011,18(3):334-353.

    [5] 吳一凡.GPM IMERG v5和TRMM 3B42 v7衛(wèi)星降水產(chǎn)品在長江流域的精度評價與對比[D].南京:南京林業(yè)大學(xué),2019.

    [6] HOU A Y,KAKAR R K,NEECK S,et al.The global precipitation measurement mission[J].Bulletin of the American Meteorological Society,2014,95(5):701-722.

    [7] KIRSCHBAUM D B,HUFFMAN G J,ADLER R F,et al.NASA′s remotely sensed precipitation:a reservoir for applications users[J].Bulletin of the American Meteorological Society,2017,98(6):1169-1184.

    [8] HUFFMAN G J,BOLVIN D T,BRAITHWAITE D,et al.Integrated multi-satellite retrievals for the global precipitation measurement(GPM)mission(IMERG)[M].Berlin:Springer,2020:343-353.

    [9] 許金濤.中國大陸IMERG降水產(chǎn)品精度評價與融合校正研究[D].杭州:浙江大學(xué),2020.

    [10] PERERA H,F(xiàn)ERNANDO S,GUNATHILAKE M B,et al.Evaluation of satellite rainfall products over the Mahaweli River Basin in Sri Lanka[J].Advances in Meteorology,2022:1926854.

    [11] TANG G,CLARK M P,PAPALEXIOU S M,et al.Have satellite precipitation products improved over last two decades? A comprehensive comparison of GPM IMERG with nine satellite and reanalysis datasets[J].Remote Sensing of Environment,2020,240:111697.

    [12] JIANG S,WEI L,REN L,et al.Utility of integrated IMERG precipitation and GLEAM potential evapotranspiration products for drought monitoring over mainland China[J].Atmospheric Research,2021,247:105141.

    [13] 孔宇.中國大陸GPM/IMERG產(chǎn)品的精度評估[D].南京:南京信息工程大學(xué),2017.

    [14] 張茹,雍斌,曾歲康.GPM衛(wèi)星降水產(chǎn)品在中國大陸的精度評估[J].人民長江,2021,52(5):50-59.

    [15] 劉若蘭,江善虎,任立良,等.全球降水觀測計劃IMERG降水產(chǎn)品對中國大陸極端降雨監(jiān)測能力評估[J].中國農(nóng)村水利水電,2021,(4):57-63.

    [16] FANG J,YANG W,LUAN Y,et al.Evaluation of the TRMM 3B42 and GPM IMERG products for extreme precipitation analysis over China[J].Atmospheric Research,2019,223:24-38.

    [17] 李媛媛,寧少尉,丁偉,等.最新GPM降水?dāng)?shù)據(jù)在黃河流域的精度評估[J].國土資源遙感,2019,31(1):164-170.

    [18] 楊震宇,寧少尉,金菊良.最新GSMap-gauged與GPM-IMERG衛(wèi)星降水產(chǎn)品性能評估及其在長江流域中的應(yīng)用[J].水電能源科學(xué),2018,36(11):5-8,13.

    [19] CHEN C,CHEN Q,DUAN Z,et al.Multiscale comparative evaluation of the GPM IMERG v5 and TRMM 3B42 v7 precipitation products from 2015 to 2017 over a climate transition area of China[J].Remote Sensing,2018,10(6):944.

    [20] 慶丹丹,宋凡,寧少尉,等.最新GPM遙感數(shù)據(jù)在淮河和海河流域的誤差分析[J].水文,2018,38(5):52-58.

    [21] 曾歲康,雍斌.全球降水計劃IMERG和GSMaP反演降水在四川地區(qū)的精度評估[J].地理學(xué)報,2019,74(7):1305-1318.

    [22] 覃俊凱.IMERG衛(wèi)星降水?dāng)?shù)據(jù)在巖溶區(qū)徑流模擬中的適用性研究[D].南寧:廣西大學(xué),2019.

    [23] 劉玄,唐培軍,吳同帥,等.山東省極端氣候指數(shù)變化特征研究[J].水利水運工程學(xué)報,2022(2):40-50.

    [24] 李芳,孔宇,高謙.GPM/IMERG產(chǎn)品在魯南地區(qū)的精度評估[J].氣象科技,2020,48(4):474-481.

    [25] 王思夢,王大釗,黃昌.GPM衛(wèi)星降水?dāng)?shù)據(jù)在黑河流域的適用性評價[J].自然資源學(xué)報,2018,33(10):1847-1860.

    (編輯:謝玲嫻)

    Evaluation of precipitation monitoring capability of IMERG series products in Dawen River basin

    DUAN Zhen1,2,WANG Chao2,DI Yan3,LIU Qiang2,WANG Gang2

    (1.Tai′an Hydrology Center,Tai′an 271000,China, 2.College of Water Conservation & Civil Engineering,Shandong Agricultural University,Tai′an 271000,China; 3.Feicheng Water Resources Protection Center,Tai′an 271600,China)

    Abstract:

    In order to explore the monitoring accuracy and applicability of high-resolution multi-satellite joint inversion IMERG(Integrated Multi-satellite Retrievals for GPM ) satellite precipitation products in the Dawen River basin,based on a relatively long series of measured data of ground rainfall stations from 2001 to 2019,this paper used 11 indicators such as detection rate and false alarm rate to evaluate the precipitation monitoring capabilities of Early,Late and Final products of IMERG series in the study area.The results show that:①At the daily,monthly and annual scales,the accuracy of Final products is the best,and the accuracy of Early and Late products is basically same.In terms of time,the IMERG series products can better reflect the inter-annual variation and intra-annual distribution of precipitation in the basin.In terms of space,the precipitation in the study area was overestimated and the precipitation in the Taishan mountain area was underestimated(relative deviation less than -20%).②In view of the ability to detect daily precipitation of different intensities,the IMERG series products have minor differences and only have a good hit rate for light rain(detection rate > 0.6).The detection ability of each product is general,and the overall performance is a high false alarm rate(false alarm rate > 0.7),low detection rate and key success rate.At the same time,with the increase of precipitation intensity,the detection ability of each product shows a downward trend.③The IMERG series products have certain extreme precipitation monitoring capabilities.The Final product has better evaluation results and more stable performance.Because among the five indexes participating in the evaluation,the Final product has four index gaining optimal results,the Late product has one optimal,and the Early product has four index that are better than the Late product.In general,the IMERG series data can better reflect the spatial and temporal distribution of precipitation and extreme precipitation in the Dawen River basin(except Taishan mountain area).The overall accuracy of the Final product is the highest,and the accuracies of the Early and the Late products are relatively close.Among them,the performance of the Early products for extreme precipitation monitoring is better than that of the Late product,while the monitoring ability of daily precipitation events with different intensities is worse than that of the Late product.

    Key words:

    satellite remote sensing precipitation inversion;IMERG;accuracy evaluation;extreme precipitation;Dawen River basin

    猜你喜歡
    大汶河降水流域
    由大汶河到泰山
    壓油溝小流域
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    堡子溝流域綜合治理
    羅堰小流域
    大汶河流域水資源質(zhì)量狀況及變化趨勢分析
    中國水運(2020年11期)2020-12-23 04:22:54
    大汶河梯級攔蓄水庫富營養(yǎng)化評估
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    打造智慧流域的思路及構(gòu)想——以討賴河流域為例
    降水現(xiàn)象儀模擬軟件設(shè)計與實現(xiàn)
    亚洲精品自拍成人| 视频区图区小说| 欧美日韩福利视频一区二区| 日韩中文字幕欧美一区二区 | 亚洲国产精品一区二区三区在线| 国产深夜福利视频在线观看| 亚洲精品国产av成人精品| 国产99久久九九免费精品| 啦啦啦中文免费视频观看日本| 欧美精品一区二区免费开放| 19禁男女啪啪无遮挡网站| h视频一区二区三区| 午夜影院在线不卡| 免费少妇av软件| 女人高潮潮喷娇喘18禁视频| 国产欧美日韩一区二区三区在线| 欧美av亚洲av综合av国产av | 别揉我奶头~嗯~啊~动态视频 | 美女扒开内裤让男人捅视频| 亚洲天堂av无毛| 亚洲免费av在线视频| 国产精品偷伦视频观看了| 国产精品 国内视频| 午夜日本视频在线| 哪个播放器可以免费观看大片| 免费观看a级毛片全部| 女人精品久久久久毛片| 黑人猛操日本美女一级片| 亚洲 欧美一区二区三区| 韩国av在线不卡| 女人高潮潮喷娇喘18禁视频| 秋霞伦理黄片| 最近中文字幕2019免费版| 国产精品av久久久久免费| 久久热在线av| 国产精品 欧美亚洲| 黄片播放在线免费| 少妇猛男粗大的猛烈进出视频| 国语对白做爰xxxⅹ性视频网站| 777久久人妻少妇嫩草av网站| 婷婷色av中文字幕| 国产野战对白在线观看| 亚洲人成网站在线观看播放| 欧美激情极品国产一区二区三区| 老司机影院毛片| 男女边摸边吃奶| 在线观看免费午夜福利视频| 日韩中文字幕视频在线看片| 久久热在线av| 亚洲欧美清纯卡通| 一边摸一边做爽爽视频免费| 亚洲美女搞黄在线观看| 久久婷婷青草| 五月开心婷婷网| 国产成人精品久久二区二区91 | 亚洲精品国产一区二区精华液| 青青草视频在线视频观看| 一本久久精品| 久久毛片免费看一区二区三区| 一个人免费看片子| 亚洲,欧美精品.| 大陆偷拍与自拍| 亚洲国产精品一区二区三区在线| 国产精品一国产av| 一本色道久久久久久精品综合| 欧美少妇被猛烈插入视频| 在线观看免费高清a一片| 欧美 亚洲 国产 日韩一| 免费黄色在线免费观看| 亚洲欧美一区二区三区久久| 午夜激情av网站| 国产爽快片一区二区三区| 亚洲国产看品久久| 亚洲久久久国产精品| 男女之事视频高清在线观看 | 亚洲第一区二区三区不卡| 亚洲第一区二区三区不卡| 日韩大码丰满熟妇| 国产精品秋霞免费鲁丝片| 亚洲精品国产一区二区精华液| 啦啦啦视频在线资源免费观看| 久久久久精品性色| 欧美日韩亚洲高清精品| 精品国产一区二区三区四区第35| 亚洲成av片中文字幕在线观看| 丝袜美腿诱惑在线| 汤姆久久久久久久影院中文字幕| 国产一区二区在线观看av| 一级毛片电影观看| 看免费av毛片| 五月天丁香电影| 如何舔出高潮| 国产日韩欧美在线精品| 啦啦啦 在线观看视频| 交换朋友夫妻互换小说| 一区二区av电影网| 激情五月婷婷亚洲| 色播在线永久视频| 一级黄片播放器| 十八禁高潮呻吟视频| 国产成人午夜福利电影在线观看| 两个人免费观看高清视频| 亚洲国产精品成人久久小说| 妹子高潮喷水视频| 亚洲av成人不卡在线观看播放网 | 波多野结衣一区麻豆| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产欧美在线一区| 天堂8中文在线网| 亚洲四区av| 国产精品一二三区在线看| 国产乱人偷精品视频| 久久精品熟女亚洲av麻豆精品| 久久鲁丝午夜福利片| 亚洲国产精品一区二区三区在线| 99香蕉大伊视频| 可以免费在线观看a视频的电影网站 | 天堂中文最新版在线下载| 国产国语露脸激情在线看| 亚洲精品国产色婷婷电影| 一区二区三区激情视频| 日本黄色日本黄色录像| 中文字幕色久视频| 亚洲精品在线美女| av.在线天堂| 色综合欧美亚洲国产小说| 丝袜脚勾引网站| 亚洲精华国产精华液的使用体验| 精品国产乱码久久久久久小说| 精品卡一卡二卡四卡免费| 精品亚洲成国产av| 男的添女的下面高潮视频| av天堂久久9| 在线亚洲精品国产二区图片欧美| 制服丝袜香蕉在线| 校园人妻丝袜中文字幕| 青春草视频在线免费观看| av在线观看视频网站免费| 卡戴珊不雅视频在线播放| 男女午夜视频在线观看| 在线观看免费日韩欧美大片| 女人爽到高潮嗷嗷叫在线视频| 国产免费又黄又爽又色| 亚洲色图 男人天堂 中文字幕| 色播在线永久视频| 久久天堂一区二区三区四区| 欧美日韩福利视频一区二区| 中文乱码字字幕精品一区二区三区| 亚洲精品一二三| 亚洲天堂av无毛| 咕卡用的链子| 亚洲视频免费观看视频| 精品一区二区免费观看| 精品少妇一区二区三区视频日本电影 | 天天操日日干夜夜撸| 久久久久精品国产欧美久久久 | 国产1区2区3区精品| 亚洲少妇的诱惑av| 人成视频在线观看免费观看| 91精品伊人久久大香线蕉| 亚洲欧美色中文字幕在线| 18在线观看网站| 啦啦啦在线免费观看视频4| 亚洲成人免费av在线播放| 一级,二级,三级黄色视频| av在线播放精品| 国产一区有黄有色的免费视频| 一级片免费观看大全| 国产成人av激情在线播放| 亚洲图色成人| 国产成人精品久久久久久| 成人影院久久| 一个人免费看片子| 九色亚洲精品在线播放| 国产成人一区二区在线| 中文字幕色久视频| 中文字幕亚洲精品专区| 午夜日本视频在线| 中文字幕精品免费在线观看视频| 国产一区二区三区综合在线观看| 亚洲一级一片aⅴ在线观看| 韩国av在线不卡| 久久久久久久久久久免费av| 久久ye,这里只有精品| 狠狠精品人妻久久久久久综合| 欧美激情高清一区二区三区 | 成人亚洲精品一区在线观看| 国产亚洲欧美精品永久| 国产日韩欧美视频二区| 免费高清在线观看日韩| 9色porny在线观看| 久久国产精品男人的天堂亚洲| 十八禁人妻一区二区| 91精品伊人久久大香线蕉| 亚洲精品久久成人aⅴ小说| av免费观看日本| 天天躁日日躁夜夜躁夜夜| 国产成人午夜福利电影在线观看| 成年女人毛片免费观看观看9 | 91成人精品电影| 精品福利永久在线观看| av电影中文网址| 亚洲国产精品999| 成人国产av品久久久| 男的添女的下面高潮视频| 亚洲精品av麻豆狂野| 日韩电影二区| 女人精品久久久久毛片| 18禁动态无遮挡网站| 午夜福利网站1000一区二区三区| 99久久99久久久精品蜜桃| 久久精品aⅴ一区二区三区四区| 欧美日韩一级在线毛片| 男人爽女人下面视频在线观看| 精品少妇黑人巨大在线播放| www.熟女人妻精品国产| 亚洲国产毛片av蜜桃av| 午夜老司机福利片| 亚洲av综合色区一区| av.在线天堂| 天美传媒精品一区二区| 巨乳人妻的诱惑在线观看| 欧美人与性动交α欧美精品济南到| 婷婷色综合www| 男女床上黄色一级片免费看| av又黄又爽大尺度在线免费看| 亚洲精品国产区一区二| 久久久久久久久久久久大奶| 亚洲成人手机| 在线精品无人区一区二区三| 精品一区在线观看国产| 国产xxxxx性猛交| 免费观看a级毛片全部| av女优亚洲男人天堂| 九草在线视频观看| 亚洲国产精品一区三区| 精品一品国产午夜福利视频| 最近最新中文字幕免费大全7| 亚洲三区欧美一区| av网站免费在线观看视频| 少妇猛男粗大的猛烈进出视频| 国产av一区二区精品久久| 少妇人妻 视频| 精品少妇一区二区三区视频日本电影 | 别揉我奶头~嗯~啊~动态视频 | 我要看黄色一级片免费的| 国产精品二区激情视频| 熟女少妇亚洲综合色aaa.| 三上悠亚av全集在线观看| 久久天躁狠狠躁夜夜2o2o | 男女午夜视频在线观看| 亚洲欧洲国产日韩| 国产精品久久久久久久久免| 1024视频免费在线观看| 我要看黄色一级片免费的| www.av在线官网国产| √禁漫天堂资源中文www| 日韩制服丝袜自拍偷拍| 最黄视频免费看| 欧美精品亚洲一区二区| 国产免费福利视频在线观看| 国产av国产精品国产| 考比视频在线观看| 大片电影免费在线观看免费| 涩涩av久久男人的天堂| 精品一区二区三区四区五区乱码 | 无遮挡黄片免费观看| 免费观看a级毛片全部| 九色亚洲精品在线播放| 亚洲国产精品一区二区三区在线| 在线观看免费高清a一片| 亚洲精华国产精华液的使用体验| 黑人猛操日本美女一级片| 91成人精品电影| 欧美日韩一级在线毛片| 韩国av在线不卡| 高清黄色对白视频在线免费看| 啦啦啦视频在线资源免费观看| 亚洲成人国产一区在线观看 | 国产黄色视频一区二区在线观看| 满18在线观看网站| 我要看黄色一级片免费的| 欧美激情极品国产一区二区三区| 大码成人一级视频| 国产日韩欧美视频二区| 国产福利在线免费观看视频| 我要看黄色一级片免费的| 国产在线视频一区二区| 色视频在线一区二区三区| 天天躁夜夜躁狠狠躁躁| 男女下面插进去视频免费观看| 午夜日韩欧美国产| 亚洲欧美成人精品一区二区| 午夜老司机福利片| 国产精品蜜桃在线观看| 91aial.com中文字幕在线观看| 97人妻天天添夜夜摸| 99精国产麻豆久久婷婷| 亚洲国产欧美网| 十八禁网站网址无遮挡| 亚洲情色 制服丝袜| 午夜免费观看性视频| 最新在线观看一区二区三区 | 在线观看一区二区三区激情| 最新在线观看一区二区三区 | 在线观看免费视频网站a站| 国产97色在线日韩免费| 精品亚洲乱码少妇综合久久| 人人妻人人澡人人看| 热re99久久精品国产66热6| 亚洲欧美精品自产自拍| 亚洲精品国产av成人精品| 女性生殖器流出的白浆| 久久97久久精品| 国产成人一区二区在线| 人人澡人人妻人| 国产av一区二区精品久久| 自拍欧美九色日韩亚洲蝌蚪91| 涩涩av久久男人的天堂| 色吧在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产97色在线日韩免费| 波多野结衣av一区二区av| 国产国语露脸激情在线看| 如何舔出高潮| 操出白浆在线播放| 久久久久久久大尺度免费视频| 国产欧美日韩一区二区三区在线| 欧美人与性动交α欧美精品济南到| 国产精品一二三区在线看| 日韩av不卡免费在线播放| a级毛片黄视频| 51午夜福利影视在线观看| 一级毛片我不卡| 欧美黄色片欧美黄色片| 青青草视频在线视频观看| kizo精华| 久久午夜综合久久蜜桃| 1024香蕉在线观看| 久久ye,这里只有精品| 两个人免费观看高清视频| 日韩精品有码人妻一区| 久久久久久久国产电影| 日韩伦理黄色片| 另类精品久久| 欧美人与性动交α欧美精品济南到| 啦啦啦啦在线视频资源| 免费观看a级毛片全部| 亚洲av综合色区一区| 一区二区三区精品91| 少妇 在线观看| 男的添女的下面高潮视频| 成年人午夜在线观看视频| 你懂的网址亚洲精品在线观看| 99九九在线精品视频| 精品酒店卫生间| 啦啦啦在线免费观看视频4| 老汉色av国产亚洲站长工具| 久久综合国产亚洲精品| 亚洲成av片中文字幕在线观看| 亚洲国产精品一区三区| 精品少妇黑人巨大在线播放| 中国国产av一级| 嫩草影视91久久| 色94色欧美一区二区| 国产精品欧美亚洲77777| 女人高潮潮喷娇喘18禁视频| 中国三级夫妇交换| 午夜福利视频在线观看免费| 亚洲免费av在线视频| 男女午夜视频在线观看| 久久久久国产一级毛片高清牌| 午夜影院在线不卡| 日本猛色少妇xxxxx猛交久久| 成年动漫av网址| 啦啦啦在线观看免费高清www| 久久国产精品男人的天堂亚洲| 在线观看三级黄色| h视频一区二区三区| 男女边吃奶边做爰视频| 国产一卡二卡三卡精品 | 欧美精品高潮呻吟av久久| 亚洲一级一片aⅴ在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品国产一区二区精华液| 成年动漫av网址| √禁漫天堂资源中文www| 国产av一区二区精品久久| 国产视频首页在线观看| 99国产精品免费福利视频| 国产人伦9x9x在线观看| 深夜精品福利| av一本久久久久| 高清在线视频一区二区三区| av片东京热男人的天堂| 久久久国产欧美日韩av| 少妇被粗大猛烈的视频| 日日啪夜夜爽| 中文字幕色久视频| 丝瓜视频免费看黄片| 国产有黄有色有爽视频| 五月天丁香电影| 久久精品国产亚洲av高清一级| 国精品久久久久久国模美| 国产精品久久久人人做人人爽| 亚洲精品乱久久久久久| 国产又爽黄色视频| 免费久久久久久久精品成人欧美视频| 午夜免费男女啪啪视频观看| 天堂8中文在线网| 久久鲁丝午夜福利片| 久久国产亚洲av麻豆专区| 国产在线免费精品| 王馨瑶露胸无遮挡在线观看| 51午夜福利影视在线观看| 国产精品久久久久久精品电影小说| 日韩人妻精品一区2区三区| 97精品久久久久久久久久精品| 亚洲av中文av极速乱| 在线观看www视频免费| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻,人人澡人人爽秒播 | 色视频在线一区二区三区| 九草在线视频观看| 人妻人人澡人人爽人人| 男的添女的下面高潮视频| 精品少妇一区二区三区视频日本电影 | 黑人猛操日本美女一级片| 1024香蕉在线观看| 久久99一区二区三区| 成人漫画全彩无遮挡| 国产毛片在线视频| 日日摸夜夜添夜夜爱| 伦理电影大哥的女人| 天天躁日日躁夜夜躁夜夜| av有码第一页| 免费在线观看视频国产中文字幕亚洲 | 久久久国产欧美日韩av| 国产精品三级大全| 国产伦理片在线播放av一区| 高清不卡的av网站| 香蕉丝袜av| 久久久久网色| 一区二区三区精品91| 久久久久久人人人人人| 色婷婷av一区二区三区视频| 久久婷婷青草| 色视频在线一区二区三区| 看十八女毛片水多多多| 欧美国产精品va在线观看不卡| 国产成人精品久久久久久| 国产黄色视频一区二区在线观看| 波多野结衣一区麻豆| 亚洲伊人色综图| 少妇的丰满在线观看| 午夜福利乱码中文字幕| 黄色怎么调成土黄色| 久久精品人人爽人人爽视色| 亚洲精品日本国产第一区| 18在线观看网站| 国产精品蜜桃在线观看| 亚洲成国产人片在线观看| 午夜免费观看性视频| 青青草视频在线视频观看| 色94色欧美一区二区| 免费日韩欧美在线观看| 久久久久久久久久久久大奶| 老司机靠b影院| 日韩伦理黄色片| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕最新亚洲高清| 女性生殖器流出的白浆| www日本在线高清视频| 亚洲精品自拍成人| 晚上一个人看的免费电影| av片东京热男人的天堂| 国产av一区二区精品久久| 久久久久久久精品精品| 亚洲成人av在线免费| 老司机亚洲免费影院| 国产黄色免费在线视频| 久久久久久人人人人人| 在线观看三级黄色| 亚洲,欧美,日韩| 精品少妇内射三级| 国产有黄有色有爽视频| 精品少妇一区二区三区视频日本电影 | 日本av手机在线免费观看| 久久久久久久久久久免费av| 一边摸一边抽搐一进一出视频| 街头女战士在线观看网站| 熟妇人妻不卡中文字幕| 一边亲一边摸免费视频| 日本黄色日本黄色录像| 无限看片的www在线观看| 麻豆av在线久日| 欧美成人精品欧美一级黄| 一本—道久久a久久精品蜜桃钙片| 日韩欧美一区视频在线观看| 久久精品久久精品一区二区三区| 蜜桃在线观看..| 亚洲 欧美一区二区三区| 在线天堂最新版资源| 精品少妇一区二区三区视频日本电影 | av在线播放精品| 天堂中文最新版在线下载| 亚洲熟女毛片儿| 菩萨蛮人人尽说江南好唐韦庄| 午夜免费观看性视频| 国产av国产精品国产| 满18在线观看网站| 国产精品久久久久久久久免| xxx大片免费视频| 久久亚洲国产成人精品v| 色视频在线一区二区三区| 大陆偷拍与自拍| 侵犯人妻中文字幕一二三四区| 美女主播在线视频| 精品国产一区二区三区四区第35| 久久青草综合色| 欧美日韩一区二区视频在线观看视频在线| 麻豆精品久久久久久蜜桃| 日韩制服骚丝袜av| 搡老岳熟女国产| 综合色丁香网| 免费观看人在逋| 韩国av在线不卡| 色婷婷av一区二区三区视频| 18禁国产床啪视频网站| 亚洲成人手机| 又大又爽又粗| 涩涩av久久男人的天堂| 久久热在线av| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩欧美视频二区| 最新在线观看一区二区三区 | 国产高清不卡午夜福利| 亚洲精品第二区| 菩萨蛮人人尽说江南好唐韦庄| 久久毛片免费看一区二区三区| 午夜日韩欧美国产| 国产黄色免费在线视频| 男女午夜视频在线观看| 黑人欧美特级aaaaaa片| 国产伦人伦偷精品视频| 热re99久久精品国产66热6| 各种免费的搞黄视频| 一本一本久久a久久精品综合妖精| 欧美日韩国产mv在线观看视频| 亚洲激情五月婷婷啪啪| 亚洲少妇的诱惑av| 久久热在线av| 制服丝袜香蕉在线| 亚洲精品中文字幕在线视频| 久久久久国产一级毛片高清牌| 不卡视频在线观看欧美| 亚洲精品日本国产第一区| 少妇被粗大猛烈的视频| 久热这里只有精品99| 99久久综合免费| 亚洲av国产av综合av卡| 亚洲av在线观看美女高潮| 久久久精品94久久精品| 人人妻人人澡人人爽人人夜夜| 亚洲精品成人av观看孕妇| 免费不卡黄色视频| 亚洲男人天堂网一区| 免费在线观看视频国产中文字幕亚洲 | 一本一本久久a久久精品综合妖精| 我要看黄色一级片免费的| 日本色播在线视频| 久久狼人影院| 丝袜美足系列| 免费人妻精品一区二区三区视频| tube8黄色片| 人人妻人人添人人爽欧美一区卜| 最近最新中文字幕大全免费视频 | 久久久精品免费免费高清| 丰满迷人的少妇在线观看| 欧美激情极品国产一区二区三区| 人妻人人澡人人爽人人| 人成视频在线观看免费观看| 久久人人爽人人片av| 日本一区二区免费在线视频| 亚洲成人免费av在线播放| 免费av中文字幕在线| 国产精品久久久人人做人人爽| 天堂8中文在线网| 亚洲第一青青草原| 久久人人爽av亚洲精品天堂| 精品人妻熟女毛片av久久网站| 少妇猛男粗大的猛烈进出视频| 中文字幕人妻丝袜制服| av女优亚洲男人天堂| 欧美日韩成人在线一区二区| 如何舔出高潮| 色婷婷av一区二区三区视频| 精品视频人人做人人爽| 亚洲情色 制服丝袜| 久久久久久久久久久免费av| 欧美 亚洲 国产 日韩一| 亚洲精品视频女| 国产一级毛片在线| 色吧在线观看| 国产精品国产av在线观看| 啦啦啦中文免费视频观看日本| 最近中文字幕2019免费版| 久久久久精品性色| 国产一区二区三区av在线| 国产乱人偷精品视频| 国产激情久久老熟女| 99re6热这里在线精品视频| 午夜91福利影院| 欧美精品一区二区大全| 国产精品一二三区在线看|