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

    多種降水再分析數(shù)據(jù)在長江流域的適用性對比

    2022-10-06 03:12:26旭,郭力,,3,湯陽,張榮,郭
    人民長江 2022年9期
    關(guān)鍵詞:長江流域氣象站降水量

    丁 光 旭,郭 家 力,,3,湯 正 陽,張 海 榮,郭 豪 豪

    (1.三峽大學(xué) 水利與環(huán)境學(xué)院,湖北 宜昌 443002; 2.中國長江電力股份有限公司 智慧長江與水電科學(xué)湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 宜昌 443000; 3.三峽大學(xué) 三峽庫區(qū)生態(tài)環(huán)境教育部工程研究中心,湖北 宜昌 443002)

    0 引 言

    長江是世界第三大河流,長江流域人口密集,農(nóng)業(yè)發(fā)達(dá),水利工程眾多(如三峽工程、南水北調(diào)中線工程等),但流域內(nèi)地形和氣候復(fù)雜、自然災(zāi)害頻發(fā)[1]。獲取準(zhǔn)確的降水信息對長江流域水資源管理、農(nóng)業(yè)生產(chǎn)、旱澇災(zāi)害防控等領(lǐng)域都起著至關(guān)重要的作用。

    地面站點(diǎn)雨量計(jì)觀測是傳統(tǒng)獲取降水資料最直接有效的辦法,這種方法獲得的降水資料精度較高,然而流域內(nèi)降水過程一般較復(fù)雜且具有高度時空變異性,再加上流域內(nèi)雨量站分布不均,使得很難準(zhǔn)確獲取地理位置偏遠(yuǎn)且地形復(fù)雜地區(qū)的降水資料。隨著衛(wèi)星遙感技術(shù)和數(shù)據(jù)反演算法的發(fā)展,同化了衛(wèi)星及地面、高空等常規(guī)觀測資料的氣象再分析資料得到深入應(yīng)用,再分析資料已成為新的降水?dāng)?shù)據(jù)來源。這些降水再分析數(shù)據(jù)可在一定程度上彌補(bǔ)傳統(tǒng)站點(diǎn)降水觀測的不足[2]。因此分析不同來源降水再分析數(shù)據(jù)在長江流域的誤差及其適用性,對長江流域降水?dāng)?shù)據(jù)融合、水文預(yù)報(bào)、防汛抗旱等具有重要的數(shù)據(jù)支撐作用。

    近幾年來,一些具有較高時空分辨率的降水再分析數(shù)據(jù)相繼發(fā)布,如中國區(qū)域高時空分辨率地面氣象要素驅(qū)動數(shù)據(jù)集(China Meteorological Forcing Dataset,CMFD)、歐洲中期天氣預(yù)報(bào)中心全球氣候第五代產(chǎn)品ERA5以及多源加權(quán)集合降水?dāng)?shù)據(jù)(Multi-Source Weighted-Ensemble Precipitation,MSWEP)等。目前這3種降水?dāng)?shù)據(jù)已有針對國內(nèi)一些區(qū)域的應(yīng)用研究,鄧越等[3]基于1979~2015年中國大陸824個氣象站點(diǎn)觀測數(shù)據(jù),評估MSWEP在中國大陸區(qū)域的精度,發(fā)現(xiàn)MSWEP在中國大陸區(qū)域降水整體上存在高估現(xiàn)象。黃琦等[4]對MSWEP在祁連山地區(qū)的適用性進(jìn)行了分析,結(jié)果表明訂正后的MSWEP與地表降水更為吻合。羅瀟等[5]對比分析了CMFD與CRA-40在湖南省的適用性,研究表明兩套資料均能較好地反映湖南省降水空間分布特征。黃琦等[6]在雅礱江流域?qū)?種降水產(chǎn)品進(jìn)行適用性分析,發(fā)現(xiàn)CMFD在各時間尺度表現(xiàn)均優(yōu)于IMERG和MSWEP。崔豪等[7]對CMFD在三峽庫區(qū)的適應(yīng)性進(jìn)行分析,發(fā)現(xiàn)CMFD能較好地反映三峽庫區(qū)內(nèi)降水的空間分布特征;Hu等[8]研究發(fā)現(xiàn),ERA5在青藏高原雖然能夠反映降水總體時空特征,但仍存在明顯的偏差。Jiang等[9]分析評估了ERA5在中國大陸地區(qū)的適用性,發(fā)現(xiàn)ERA5對中國大陸地區(qū)的降水有較大的偏差。Xin等[10]對ERA5在粵港澳大灣區(qū)的適用性進(jìn)行了分析,發(fā)現(xiàn)ERA5降水?dāng)?shù)據(jù)在城鎮(zhèn)化水平較高的地區(qū)精度不高,還需進(jìn)一步改進(jìn)。

    對于降水再分析數(shù)據(jù)產(chǎn)品的適應(yīng)性分析,不少學(xué)者僅針對單一降水產(chǎn)品進(jìn)行研究,而對于同一流域或區(qū)域的不同種降水再分析數(shù)據(jù),由于選用不同模式運(yùn)算、觀測資料及同化方法等帶來的差異,使得不同數(shù)據(jù)之間的可信度也存在差異[11]。因此有必要對多種降水再分析數(shù)據(jù)之間的差異性進(jìn)行對比研究。此外,對于CMFD、MSWEP降水再分析數(shù)據(jù),雖然有對其在國內(nèi)一些流域或區(qū)域進(jìn)行適應(yīng)性研究[3-7],但針對整個長江流域的研究較少。特別是對于ERA5,其包括溫度、濕度、壓強(qiáng)等多種要素,而對降水要素的適應(yīng)性分析還只是針對青藏高原[8]、中國大陸區(qū)域[9]、粵港澳大灣區(qū)[10]等,對于整個長江流域的適應(yīng)性研究較少。因此,本文擬采用1979~2016年長江流域226個地面氣象站點(diǎn)資料作為基準(zhǔn)資料,對比分析CMFD、MSWEP、ERA5這3種降水再分析產(chǎn)品在不同時間尺度以及不同量級降水強(qiáng)度下的精度,分析3種降水?dāng)?shù)據(jù)對比地面站點(diǎn)降水精度在長江流域的分布情況,并探索誤差的時空分布規(guī)律,為長江流域降水及水文模擬等相關(guān)研究與應(yīng)用選擇可靠的降水再分析產(chǎn)品提供依據(jù)。

    1 研究區(qū)域、數(shù)據(jù)和方法

    1.1 研究區(qū)域概況

    長江流域位于24°27′N~35°54′N,90°33′E~122°19′E之間,流域面積180萬 km2,占全國總面積的18.75%,流域橫跨青藏高原、四川盆地和長江中下游平原三大階梯地貌,流域地形復(fù)雜,夏季從高原季風(fēng)氣候過渡到東亞季風(fēng)氣候,雨量充沛,極易發(fā)生洪澇災(zāi)害[12]。本次研究一共選取流域內(nèi)226個地面氣象站點(diǎn),地面氣象站點(diǎn)分布如圖1所示。

    1.2 研究數(shù)據(jù)

    研究所采用地面氣象站點(diǎn)資料來自國家氣象信息中心的中國地面氣候資料日值數(shù)據(jù)集(V3.0),觀測數(shù)據(jù)經(jīng)過氣候界限值檢查、臺站極值檢查、時空一致性檢查等一系列質(zhì)量控制。降水再分析資料包括中國區(qū)域高時空分辨率地面氣象要素驅(qū)動數(shù)據(jù)集CMFD、多源加權(quán)集合降水?dāng)?shù)據(jù)MSWEP以及歐洲中期天氣預(yù)報(bào)中心ERA5降水?dāng)?shù)據(jù)。

    多源加權(quán)集合降水?dāng)?shù)據(jù)集(MSWEP)[13]結(jié)合站點(diǎn)數(shù)據(jù)、衛(wèi)星觀測數(shù)據(jù)以及模型模擬數(shù)據(jù),具有時間尺度長、空間分辨率較高等優(yōu)點(diǎn)。中國區(qū)域高時空分辨率地面氣象要素驅(qū)動數(shù)據(jù)集(CMFD)[14]是中國科學(xué)院青藏高原研究所開發(fā)的一套近地面氣象與環(huán)境要素再分析數(shù)據(jù)集,包含降水在內(nèi)共7個要素,是以國際上現(xiàn)有的Princeton再分析資料、GLDAS資料、GEWEX-SRB輻射資料,以及TRMM降水資料為背景場,融合了中國氣象局常規(guī)氣象觀測數(shù)據(jù)制作而成。ERA5是歐洲中期天氣預(yù)報(bào)中心(ECMWF)發(fā)布的全球氣候第五代產(chǎn)品,其性能與上一代產(chǎn)品ERA-Interim相比有了很大提升[15-16]。數(shù)據(jù)詳細(xì)信息如表1所列。

    表1 研究數(shù)據(jù)概況Tab.1 Overview of the research data

    考慮到部分氣象站點(diǎn)資料僅更新到2017年,因此選取1979~2016年作為研究時段,其中ERA5數(shù)據(jù)空間分辨率為0.25°×0.25°,為了排除空間分辨率對降水產(chǎn)品精度評估的影響,本文將ERA5數(shù)據(jù)經(jīng)反距離加權(quán)(IDW)法插值到與MSWEP以及CMFD相同的空間分辨率0.1°×0.1°。CMFD與MSWEP時間分辨率為3 h,而ERA5時間分辨率為1 h,因此計(jì)算時分別將3 h和1 h降水量累加到日、月、年值降水量序列。

    1.3 研究方法

    1.3.1數(shù)據(jù)關(guān)聯(lián)方法

    由于將氣象站點(diǎn)數(shù)據(jù)直接插值成空間數(shù)據(jù)有一定的誤差[17],因此本文采用點(diǎn)對點(diǎn)方法進(jìn)行數(shù)據(jù)之間的關(guān)聯(lián),即利用站點(diǎn)坐標(biāo)提取其所在格網(wǎng)的降水?dāng)?shù)據(jù)值,以氣象站點(diǎn)資料為基準(zhǔn),對比分析CMFD、MSWEP、ERA5這3種降水?dāng)?shù)據(jù)的質(zhì)量精度。

    1.3.2評價指標(biāo)選取

    采用皮爾遜相關(guān)系數(shù)(R)、均方根誤差(RMSE)、平均絕對誤差(MAE)和相對偏差(BIAS)4種定量評價指標(biāo)評價降水再分析產(chǎn)品精度[18]。其中,R表示降水再分析數(shù)據(jù)與氣象站點(diǎn)觀測數(shù)據(jù)之間的線性相關(guān)程度,相關(guān)度越強(qiáng),相關(guān)系數(shù)R越接近于1或-1;RMSE表征降水再分析數(shù)據(jù)與氣象站點(diǎn)觀測數(shù)據(jù)之間的離散程度,對特大特小誤差反映非常敏感,值越大表示兩組數(shù)據(jù)越離散,最優(yōu)值為0;MAE反映降水再分析數(shù)據(jù)與氣象站點(diǎn)觀測數(shù)據(jù)之間的平均絕對偏差程度,可以避免誤差相互抵消的問題,最優(yōu)值為0;BIAS直觀反映出降水再分析數(shù)據(jù)存在的系統(tǒng)誤差大小及相對偏差程度,最優(yōu)值為0。各指標(biāo)計(jì)算公式如下:

    (1)

    (2)

    (3)

    (4)

    此外,采用探測率(Probability of Detection,POD)、空報(bào)率(False Alarm Ratio,F(xiàn)AR)、偏差評分(B)、公平技巧評分(Equitable Threat Score,ETS)4個分類統(tǒng)計(jì)指標(biāo)評估3種降水再分析數(shù)據(jù)對不同量級降水事件的捕捉能力。其中,POD表示降水再分析數(shù)據(jù)準(zhǔn)確捕捉實(shí)際降水事件的能力;FAR表示降水再分析數(shù)據(jù)對降水事件的空報(bào)情況;B表示降水再分析數(shù)據(jù)是否高估(B>1)或者低估(B<1)了降水事件;ETS越大表示降水再分析數(shù)據(jù)對降水事件的綜合預(yù)測能力越強(qiáng)。其計(jì)算公式如下:

    (5)

    (6)

    (7)

    (8)

    (9)

    式中:H代表命中降水?dāng)?shù)量;F代表空報(bào)降水?dāng)?shù)量;M代表漏報(bào)降水?dāng)?shù)量。POD和FAR的取值介于[0,1],B取值介于[0,+∞),ETS取值介于[-1/3,1];最優(yōu)值POD=1,F(xiàn)AR=0,B=1,ETS=1。

    依據(jù)國家氣象局頒布的降水強(qiáng)度等級劃分標(biāo)準(zhǔn),選取降雨強(qiáng)度0.1~9.9,10.0~24.9,25.0~49.9,50.0~99.9,>100.0 mm/d作為發(fā)生小雨、中雨、大雨、暴雨及以上量級降水的判定標(biāo)準(zhǔn)。其中,以0.1 mm/d作為降水事件發(fā)生與否的閾值。

    2 結(jié)果分析

    2.1 多源降水產(chǎn)品年系列變化分析

    圖2為流域內(nèi)226個氣象站點(diǎn)與各降水再分析數(shù)據(jù)1979~2016年多年平均降水量散點(diǎn)圖,其中每個點(diǎn)即代表一個站點(diǎn),黑色虛線為理想狀態(tài)分布線。可以看出ERA5多年平均降水量在絕大部分站點(diǎn)上均高估站點(diǎn)降水量,而MSWEP以及CMFD則能較密集且均勻地分布在理想分布線兩側(cè),兩者相關(guān)系數(shù)分別為0.93 與0.98,高于ERA5的0.82。

    為了更直觀地表示不同來源降水?dāng)?shù)據(jù)的空間差異,分別基于長江流域1979~2016年地面氣象站點(diǎn)數(shù)據(jù)以及同期的CMFD、MSWEP和ERA5降水?dāng)?shù)據(jù)繪制多年平均降水量空間分布圖,結(jié)果如圖3所示。從圖3中可以看出,長江流域多年平均降水量總體上呈現(xiàn)東高西低的趨勢,降水的時空分布不均勻,東部地區(qū)多年平均降水量大于2 000 mm,西部地區(qū)則不足400 mm。ERA5在長江流域絕大部分地區(qū)都呈現(xiàn)高估的現(xiàn)象,部分站點(diǎn)高估程度甚至超過100%,高估較為嚴(yán)重。MSWEP在橫斷山脈、四川盆地以及長江下游表現(xiàn)為低估現(xiàn)象。CMFD在總體上則能較好地把握長江流域的降水空間分布。

    表2為長江流域226個站點(diǎn)多年平均降水量在各相對偏差區(qū)間所占比重分析,可以看出ERA5對長江流域降水總體表現(xiàn)為高估現(xiàn)象,這與Jiang等[9]在中國大陸地區(qū)的研究結(jié)論相似。其中僅有3個站點(diǎn),占站點(diǎn)總數(shù)的1.32%,相對偏差處于-20%~20%(可接受的范圍)之間,其余站點(diǎn)相對偏差均大于20%,其中有48.23%的站點(diǎn)降水值高估程度處于50%~100%之間,有15.04%的站點(diǎn)降水值高估程度達(dá)到100%以上。對于MSWEP,其絕大部分站點(diǎn)(92.03%)的相對偏差處于-20%~20%之間,其中相對偏差低于0的站點(diǎn)占到總數(shù)的59.19%,說明MSWEP對大部分站點(diǎn)表現(xiàn)為低估的現(xiàn)象。同樣對于CMFD,有99.12%的站點(diǎn)相對偏差處于-20%~20%之間,其中相對偏差處于0~20%之間的站點(diǎn)占到總數(shù)的71.24%,說明其整體上表現(xiàn)為略微高估的現(xiàn)象。相較于ERA5,CMFD與MSWEP這兩種數(shù)據(jù)對長江流域降水值估算較好,精度較高,其中CMFD優(yōu)于MSWEP。

    表2 站點(diǎn)在各相對偏差區(qū)間所占比重Tab.2 The proportion of the stations in each BIAS interval %

    2.2 多源降水產(chǎn)品月系列比較分析

    根據(jù)長江流域226個站點(diǎn)算術(shù)面平均逐月降水量及其對應(yīng)格點(diǎn)的CMFD、MSWEP、ERA5算術(shù)面平均逐月降水量,繪制多年月平均降水年內(nèi)分布圖,如圖4所示??梢钥闯?,3種降水產(chǎn)品中CMFD與MSWEP能較好地捕捉月平均降水年內(nèi)分布特征,基本與站點(diǎn)月平均降水過程線重合,精度較高,其月平均降水在6月份達(dá)到最大,12月份達(dá)到最低。而對于ERA5降水?dāng)?shù)據(jù),其各月平均降水均高于站點(diǎn)降水,雖然ERA5在長江流域月尺度上精度不高,但其年內(nèi)各月降水分布趨勢與地面站點(diǎn)降水一致。

    依據(jù)長江流域226個站點(diǎn)算術(shù)面平均逐月降水量與其對應(yīng)格點(diǎn)的CMFD、MSWEP、ERA5逐月降水量,計(jì)算各定量評價指標(biāo),結(jié)果如圖5所示,可以明顯看出ERA5的4種評價指標(biāo)表現(xiàn)均較差。就相關(guān)系數(shù)R而言,CMFD與MSWEP表現(xiàn)均較好,其中CMFD的相關(guān)系數(shù)各月均在0.9以上,平均值達(dá)到0.95,MSWEP的平均值為0.92,相比之下,ERA5相關(guān)系數(shù)平均值只有0.61。就均方根誤差RMSE而言,CMFD最大值為7月的46.42 mm,最小值為12月的10.62 mm,MSWEP最大值為7月的58.60 mm,最小值為12月的10.71 mm,ERA5最大值為8月127.70 mm,最小值為2月的55.41 mm。同樣對于平均絕對誤差MAE,3種降水再分析產(chǎn)品最大最小值發(fā)生月份與均方根誤差相同,表明3種降水再分析產(chǎn)品誤差在夏季月份要大于冬季月份。而就相對偏差BIAS而言,ERA5在11月至次年1月偏差最大,主要原因是在冬季月份降水量較少,在誤差相同的情況下會導(dǎo)致偏差程度較其他月份高。綜上可以看出,CMFD在月尺度上性能表現(xiàn)最好,其與各月份站點(diǎn)降水量擬合程度較高,相關(guān)系數(shù)最大,均方根誤差以及平均絕對誤差最小,MSWEP與ERA5次之。

    從圖5可以看出,汛期的月份和非汛期的月份3種降水產(chǎn)品的差異較大。因此,本文以5~10月作為汛期,11月至次年4月作為非汛期[19],進(jìn)一步評估了3種降水再分析數(shù)據(jù)在汛期與非汛期的適用性。表3列出了3種降水產(chǎn)品在汛期與非汛期的精度對比,可以看出,3種降水產(chǎn)品在汛期與非汛期的精度均表現(xiàn)為CMFD>MSWEP>ERA5,其中CMFD與MSWEP的各指標(biāo)值明顯優(yōu)于ERA5。就相對偏差BIAS而言,CMFD在汛期與非汛期均表現(xiàn)為輕微高估,MSWEP則表現(xiàn)為輕微低估。而ERA5均高估了汛期與非汛期的降水量,其BIAS值在汛期高估程度為47.74%,在非汛期高估程度達(dá)到88.82%。3種降水產(chǎn)品除相對偏差BIAS外,其余3個指標(biāo)均是汛期表現(xiàn)要劣于非汛期。

    表3 汛期與非汛期3種降水產(chǎn)品精度指標(biāo)Tab.3 The accuracy indices for three precipitation products during flood season and non-flood season

    2.3 多源降水產(chǎn)品日系列比較分析

    基于日尺度序列計(jì)算3種降水產(chǎn)品相對226個站點(diǎn)觀測數(shù)據(jù)的相關(guān)系數(shù)(R)、均方根誤差(RMSE)、平均絕對誤差(MAE)、相對偏差(BIAS)等指標(biāo),各指標(biāo)空間分布如圖6所示。圖7為日尺度上各評價指標(biāo)箱形圖。通過對比發(fā)現(xiàn),CMFD在長江流域表現(xiàn)出較高的相關(guān)性,其相關(guān)系數(shù)變化范圍在0.65~0.97之間,平均值為0.87,遠(yuǎn)大于MSWEP的0.51與ERA5的0.49,原因可能是CMFD為中國區(qū)域高時空分辨率數(shù)據(jù)集,相較于MSWEP與ERA5,其結(jié)合了更多的中國地面氣象站點(diǎn),對降水的把握更精準(zhǔn)[20]。對于MSWEP來說,其相關(guān)系數(shù)在金沙江下游及烏江流域較低,其值約在0.20~0.40之間。該地區(qū)地形相對復(fù)雜,可能導(dǎo)致該地區(qū)誤差較大,該發(fā)現(xiàn)與鄧越等[3]在中國大陸區(qū)域的研究結(jié)果相似。而ERA5則表現(xiàn)出相關(guān)系數(shù)自西向東逐漸增大的趨勢。從均方根誤差以及平均絕對誤差來分析,3種產(chǎn)品都大致表現(xiàn)為自西向東逐漸增大的趨勢,其中MSWEP與ERA5對這一趨勢表現(xiàn)較為明顯,這一現(xiàn)象與降水強(qiáng)度分布密切相關(guān),長江下游降水強(qiáng)度大,降水量變化也比較大,最終使得誤差也比較大。而長江上游流域地勢較高,全年降水量較少,也就使得誤差較小[21]。而CMFD無論是均方根誤差還是平均絕對誤差,其值在整個流域上都明顯低于另外兩種降水產(chǎn)品,對比圖7(b)與圖7(c)也能印證這一點(diǎn)。就相對偏差而言,CMFD與MSWEP都有較高的精確度,兩者絕大部分站點(diǎn)的相對偏差都集中分布在-20%~20%之間,其平均值分別為1.09%和-1.18%。從圖7(d)也能看出CMFD與MSWEP相對偏差值分布在0的兩側(cè),因此平均值接近于0。從圖6(l)可以看出,ERA5總體上高估日降水量,其各站點(diǎn)相對偏差平均值達(dá)到68%,高估較為嚴(yán)重。從空間上看,高估嚴(yán)重(BIAS>100%)的站點(diǎn)大多位于金沙江流域,該地區(qū)地處橫斷山脈,復(fù)雜的地形限制了衛(wèi)星對降水的反演,阻礙了再分析資料的優(yōu)化處理,趙平偉等[22]的研究結(jié)果也證實(shí)了這一點(diǎn)。

    圖8給出了3種降水產(chǎn)品在0.1 mm/d閾值下各分類統(tǒng)計(jì)指標(biāo)的空間分布圖,探測率(POD)、空報(bào)率(FAR)、偏差評分(B)以及公平技巧評分(ETS)表征3種降水再分析數(shù)據(jù)對降水/非降水事件的捕捉能力。對比圖8(a)以及圖8(c)可以發(fā)現(xiàn),CMFD以及ERA5的探測率(POD)較高,絕大部分站點(diǎn)POD值都在0.9以上,表明CMFD與ERA5有較高的探測能力,其中ERA5的探測率(POD)又要明顯高于CMFD,主要原因在于ERA5總體上高估日降水量,使得對0.1 mm/d 閾值下的降水事件有較高的探測能力。而MSWEP的POD值則相對較低,其值主要集中在0.75~0.90之間,從空間上看,MSWEP在金沙江流域的探測率最優(yōu),主要集中在0.85~0.95之間。

    就空報(bào)率(FAR)而言,CMFD與MSWEP具有相似的空間分布格局,兩者空報(bào)率(FAR)都集中在0.16~0.50之間,平均值分別為0.33與0.35,而ERA5則表現(xiàn)為長江下游流域的FAR值小于上游流域。

    就偏差評分(B)而言,CMFD與MSWEP同樣優(yōu)于ERA5,兩者主要集中在1.0~1.8之間,表明對降水的高估。ERA5則集中在1.2~3.0之間,從空間上看呈現(xiàn)自東向西逐漸增加的趨勢。

    從圖8(j)可以看出,CMFD在烏江以及四川盆地的ETS評分要低于其他區(qū)域,其值主要集中在0.25~0.35之間,而MSWEP同樣如此,ETS評分集中在0.15~0.35之間。ERA5則表現(xiàn)為下游ETS評分要明顯高于上游流域,上游流域主要集中在0.10~0.30之間,下游流域主要集中在0.30~0.45之間。綜上,CMFD對于各評價指標(biāo)的空間分布要優(yōu)于MSWEP與ERA5,表明CMFD對于降水事件發(fā)生與否的捕捉能力最強(qiáng)。

    表4計(jì)算了CMFD、MSWEP和ERA5在不同量級降水上的POD、FAR、B以及ETS,從表中可以看出,3種降水產(chǎn)品的探測率POD均表現(xiàn)出一致的變化趨勢,隨著降水量級的增加,3種降水產(chǎn)品對降水的探測能力均減弱,但在不同量級降水上均表現(xiàn)為CMFD>ERA5>MSWEP的規(guī)律,表明CMFD對不同量級降水的探測能力最強(qiáng)。對于空報(bào)率FAR而言,CMFD的空報(bào)情況隨著降水量級的增加有所緩解,ERA5則相反,MSWEP則對大雨空報(bào)最嚴(yán)重,對小雨空報(bào)最優(yōu)。就偏差評分B而言,CMFD與MSWEP都存在小雨及中雨高估,大雨及以上量級降水低估的問題,這與武海喆等[23]在黃土高原的研究結(jié)果相似,而ERA5對不同量級降水均表現(xiàn)為高估。就ETS而言,CMFD在不同量級降水上均大于0.3,MSWEP以及ERA5都小于0.2,其中在小雨、大暴雨的評分上,MSWEP要優(yōu)于ERA5,而在中雨、大雨、暴雨上,ERA5要優(yōu)于MSWEP??傮w而言,CMFD在4種指標(biāo)上的表現(xiàn)均要優(yōu)于MSWEP與ERA5,表明CMFD對不同量級降水的捕捉能力要強(qiáng)于另外兩種產(chǎn)品。

    表4 不同量級雨量探測能力Fig.4 The detecting capacities of rainfall at different magnitudes

    3 結(jié) 論

    本文基于長江流域1979~2016年226個地面氣象站點(diǎn)以及同期的CMFD、MSWEP和ERA5這3種降水再分析數(shù)據(jù),利用定量評價指標(biāo)以及分類統(tǒng)計(jì)指標(biāo)分析研究3種降水再分析產(chǎn)品在長江流域的適用性,得到以下結(jié)論:

    (1) 3種降水?dāng)?shù)據(jù)中,ERA5在長江流域不同時間尺度上均呈現(xiàn)出嚴(yán)重高估現(xiàn)象,在使用時可考慮對其進(jìn)行校正。

    (2) 日尺度上,CMFD與站點(diǎn)觀測資料一致性最高。相關(guān)系數(shù)各站點(diǎn)平均達(dá)到0.87。3種產(chǎn)品自上游向下游均不同程度表現(xiàn)出誤差逐漸增大的趨勢。對于不同量級降水的探測能力,總體上表現(xiàn)為CMFD>ERA5>MSWEP。

    (3) 月尺度上,CMFD和MSWEP與地面站點(diǎn)降水具有較高的一致性,其相關(guān)系數(shù)明顯高于ERA5,CMFD與MSWEP在汛期與非汛期的表現(xiàn)均優(yōu)于ERA5。

    (4) CMFD與MSWEP均能較好地反映長江流域多年平均降水量的空間分布特征。

    綜上,CMFD在長江流域不同時間尺度上均具有較高的精度,能在一定程度上克服地面雨量站點(diǎn)稀缺和空間分布不均的局限性。長江作為中國第一大河,流域內(nèi)地形地勢復(fù)雜、氣候差異較大,因此在未來的研究中應(yīng)該進(jìn)一步精細(xì)化評估不同氣候地形條件對降水產(chǎn)品造成的影響。此外,針對降水產(chǎn)品的精度評估,應(yīng)當(dāng)將其延長至水文驗(yàn)證端,即將降水再分析數(shù)據(jù)引入水文模型進(jìn)行驗(yàn)證,增加降水再分析數(shù)據(jù)的實(shí)際應(yīng)用價值。

    猜你喜歡
    長江流域氣象站降水量
    繪制和閱讀降水量柱狀圖
    珠峰上架起世界最高氣象站
    走遍長江流域的英國小伙
    降水量是怎么算出來的
    啟蒙(3-7歲)(2019年8期)2019-09-10 03:09:08
    心靈氣象站
    趣味(語文)(2019年3期)2019-06-12 08:50:14
    1988—2017年呼和浩特市降水演變特征分析
    長江流域園區(qū)的府際合作研究
    基于小波變換的三江平原旬降水量主周期識別
    長江流域徑流演變規(guī)律研究
    自動氣象站應(yīng)該注意的一些防雷問題
    videossex国产| 最近中文字幕高清免费大全6| 日韩人妻高清精品专区| 尤物成人国产欧美一区二区三区| 国产一区二区亚洲精品在线观看| 纵有疾风起免费观看全集完整版 | 黄色一级大片看看| 午夜激情福利司机影院| 精品国产三级普通话版| 亚洲美女搞黄在线观看| 久久久精品欧美日韩精品| 欧美三级亚洲精品| 国产伦一二天堂av在线观看| 色网站视频免费| 精品国产一区二区三区久久久樱花 | 色播亚洲综合网| 久久久久久久午夜电影| 国产成人一区二区在线| 色播亚洲综合网| av国产免费在线观看| 麻豆一二三区av精品| 国产老妇伦熟女老妇高清| 搡女人真爽免费视频火全软件| 亚洲最大成人中文| 桃色一区二区三区在线观看| 最近视频中文字幕2019在线8| 日韩欧美 国产精品| 99久久人妻综合| 亚洲高清免费不卡视频| 男人舔女人下体高潮全视频| 国内精品一区二区在线观看| 51国产日韩欧美| 久久久精品94久久精品| 日韩欧美 国产精品| 久久久a久久爽久久v久久| 亚洲av电影在线观看一区二区三区 | 欧美日韩一区二区视频在线观看视频在线 | 床上黄色一级片| 国产精品一区二区三区四区免费观看| 国产精品美女特级片免费视频播放器| 精品国产三级普通话版| 亚洲在久久综合| 亚洲欧美日韩高清专用| 久久婷婷人人爽人人干人人爱| 国产毛片a区久久久久| 非洲黑人性xxxx精品又粗又长| 在线观看一区二区三区| 亚洲国产高清在线一区二区三| 天堂影院成人在线观看| 99久久精品国产国产毛片| 激情 狠狠 欧美| 精品一区二区免费观看| 成人毛片a级毛片在线播放| 搞女人的毛片| 国产日韩欧美在线精品| 秋霞在线观看毛片| 九色成人免费人妻av| 欧美日本亚洲视频在线播放| 国内精品一区二区在线观看| 亚洲在线自拍视频| 一级av片app| 午夜福利视频1000在线观看| 国产精品人妻久久久影院| 亚洲欧美精品自产自拍| 成人高潮视频无遮挡免费网站| 丝袜喷水一区| 99热这里只有是精品在线观看| 91精品伊人久久大香线蕉| 长腿黑丝高跟| 亚洲欧洲国产日韩| 天天躁夜夜躁狠狠久久av| 国产黄片美女视频| 久久亚洲国产成人精品v| 国产成人精品一,二区| 99热这里只有是精品在线观看| 日韩欧美精品免费久久| 免费看美女性在线毛片视频| 床上黄色一级片| 99久久精品热视频| 免费看a级黄色片| 三级男女做爰猛烈吃奶摸视频| 一级爰片在线观看| 麻豆乱淫一区二区| 日本爱情动作片www.在线观看| 99在线人妻在线中文字幕| 亚洲乱码一区二区免费版| 亚洲色图av天堂| 成人午夜精彩视频在线观看| 亚洲av电影不卡..在线观看| 免费在线观看成人毛片| 狠狠狠狠99中文字幕| 99在线人妻在线中文字幕| 日韩制服骚丝袜av| 国产av在哪里看| 午夜福利高清视频| 亚洲国产精品合色在线| 国产日韩欧美在线精品| 欧美成人免费av一区二区三区| 国内精品宾馆在线| 亚洲av电影在线观看一区二区三区 | 日本五十路高清| 麻豆av噜噜一区二区三区| 精华霜和精华液先用哪个| 国产精品乱码一区二三区的特点| 麻豆一二三区av精品| 免费黄网站久久成人精品| 免费在线观看成人毛片| 能在线免费看毛片的网站| 欧美性感艳星| 久久久久久久午夜电影| 国产一区亚洲一区在线观看| 中文欧美无线码| 国内少妇人妻偷人精品xxx网站| 精品国产一区二区三区久久久樱花 | 午夜福利在线在线| 日韩,欧美,国产一区二区三区 | 久久精品国产自在天天线| 亚洲av日韩在线播放| 精品不卡国产一区二区三区| 成人高潮视频无遮挡免费网站| 欧美高清成人免费视频www| 国产免费又黄又爽又色| 欧美3d第一页| 99热全是精品| 啦啦啦啦在线视频资源| 青春草国产在线视频| 最近2019中文字幕mv第一页| 水蜜桃什么品种好| 在线免费观看不下载黄p国产| 亚洲成人久久爱视频| 国产淫片久久久久久久久| 好男人在线观看高清免费视频| 日韩高清综合在线| 国模一区二区三区四区视频| 免费观看性生交大片5| or卡值多少钱| 成年女人永久免费观看视频| 在线免费观看的www视频| 亚洲av男天堂| 久久久久久伊人网av| 好男人在线观看高清免费视频| 欧美成人a在线观看| 日本免费a在线| 久久久久久大精品| 欧美一区二区国产精品久久精品| 99国产精品一区二区蜜桃av| 国产黄色视频一区二区在线观看 | 波多野结衣巨乳人妻| 国产大屁股一区二区在线视频| 国产日韩欧美在线精品| 成年av动漫网址| 美女高潮的动态| 在线免费十八禁| 性色avwww在线观看| 69人妻影院| 免费观看性生交大片5| 中国国产av一级| 国产视频首页在线观看| 国产一区亚洲一区在线观看| 亚洲一级一片aⅴ在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩精品有码人妻一区| 免费观看a级毛片全部| 久久久久久久午夜电影| 99久久精品国产国产毛片| 最近手机中文字幕大全| 国产亚洲av嫩草精品影院| 亚洲综合精品二区| 男人的好看免费观看在线视频| 亚洲国产色片| 日本一二三区视频观看| 亚洲四区av| 精品欧美国产一区二区三| 97热精品久久久久久| 毛片女人毛片| 能在线免费看毛片的网站| 久久久久九九精品影院| 建设人人有责人人尽责人人享有的 | 久久久久久久久久久免费av| 99久久九九国产精品国产免费| 亚洲欧美成人精品一区二区| 亚洲av成人精品一区久久| 国产 一区 欧美 日韩| 亚洲国产精品久久男人天堂| 在线播放国产精品三级| 国产亚洲av片在线观看秒播厂 | 精品久久久久久久久av| 黑人高潮一二区| 日韩三级伦理在线观看| 国产一区二区在线av高清观看| 一级毛片电影观看 | 日韩高清综合在线| 亚洲精品影视一区二区三区av| 18禁在线无遮挡免费观看视频| 午夜福利在线观看吧| 美女内射精品一级片tv| 97热精品久久久久久| 国产精品久久视频播放| 亚洲婷婷狠狠爱综合网| 国产在视频线精品| 亚洲欧美一区二区三区国产| 国产69精品久久久久777片| 亚洲久久久久久中文字幕| 大又大粗又爽又黄少妇毛片口| 特级一级黄色大片| h日本视频在线播放| .国产精品久久| 永久免费av网站大全| 午夜福利在线在线| 精品国产三级普通话版| 一边亲一边摸免费视频| 最近视频中文字幕2019在线8| 一区二区三区免费毛片| 国产麻豆成人av免费视频| 免费黄色在线免费观看| 精品久久久久久成人av| 日韩av在线免费看完整版不卡| 麻豆国产97在线/欧美| 久久久久久九九精品二区国产| 午夜日本视频在线| 欧美又色又爽又黄视频| av播播在线观看一区| 神马国产精品三级电影在线观看| 国产午夜精品久久久久久一区二区三区| 91av网一区二区| 男女下面进入的视频免费午夜| 久久婷婷人人爽人人干人人爱| 亚洲熟妇中文字幕五十中出| 一级毛片电影观看 | 在线免费观看的www视频| 在线免费观看的www视频| 噜噜噜噜噜久久久久久91| 日日干狠狠操夜夜爽| 激情 狠狠 欧美| 日本免费a在线| 精品免费久久久久久久清纯| 好男人视频免费观看在线| 一本一本综合久久| 狂野欧美白嫩少妇大欣赏| 三级国产精品片| 韩国高清视频一区二区三区| 日韩欧美国产在线观看| 午夜视频国产福利| 中国美白少妇内射xxxbb| 免费av不卡在线播放| 日韩一本色道免费dvd| 亚洲成人精品中文字幕电影| www.色视频.com| 91久久精品国产一区二区成人| 亚洲人成网站在线播| 小蜜桃在线观看免费完整版高清| 久久精品国产亚洲av涩爱| 亚洲中文字幕一区二区三区有码在线看| 国产午夜精品久久久久久一区二区三区| 美女国产视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | av国产久精品久网站免费入址| 一二三四中文在线观看免费高清| 久久久久网色| 亚洲人成网站在线播| 亚洲丝袜综合中文字幕| 看十八女毛片水多多多| 国内少妇人妻偷人精品xxx网站| 男人舔女人下体高潮全视频| 色综合站精品国产| 少妇猛男粗大的猛烈进出视频 | h日本视频在线播放| 中文字幕av成人在线电影| 亚洲欧美成人综合另类久久久 | 免费在线观看成人毛片| 最新中文字幕久久久久| 日韩一区二区三区影片| 青青草视频在线视频观看| 日韩国内少妇激情av| 久久精品国产亚洲av涩爱| 欧美区成人在线视频| 日韩精品青青久久久久久| 在线免费观看的www视频| 成人午夜高清在线视频| 大香蕉97超碰在线| 国产精品久久久久久久久免| 日韩一区二区视频免费看| 在线a可以看的网站| 国产一区二区三区av在线| 99在线人妻在线中文字幕| 日日摸夜夜添夜夜爱| 国产中年淑女户外野战色| 中文在线观看免费www的网站| 欧美日韩国产亚洲二区| 久久久午夜欧美精品| 婷婷六月久久综合丁香| 亚洲欧洲日产国产| 丰满少妇做爰视频| 三级国产精品片| 久久精品熟女亚洲av麻豆精品 | 最近2019中文字幕mv第一页| 欧美高清性xxxxhd video| 国产精品人妻久久久久久| 国产成人精品久久久久久| 久久亚洲精品不卡| 亚洲av一区综合| 国产精品美女特级片免费视频播放器| 九九爱精品视频在线观看| 国产成人一区二区在线| 1000部很黄的大片| av免费观看日本| 在线观看美女被高潮喷水网站| 亚洲av免费高清在线观看| 一卡2卡三卡四卡精品乱码亚洲| 日韩在线高清观看一区二区三区| 精品不卡国产一区二区三区| 村上凉子中文字幕在线| 亚洲国产欧洲综合997久久,| 亚洲国产精品成人综合色| 国产黄色视频一区二区在线观看 | 久久这里只有精品中国| 成人漫画全彩无遮挡| 好男人在线观看高清免费视频| 中文字幕人妻熟人妻熟丝袜美| 不卡视频在线观看欧美| 神马国产精品三级电影在线观看| 亚洲自偷自拍三级| 丝袜美腿在线中文| 最近的中文字幕免费完整| 亚洲天堂国产精品一区在线| 欧美不卡视频在线免费观看| 亚洲精品久久久久久婷婷小说 | 男女国产视频网站| 精品一区二区三区视频在线| 中文在线观看免费www的网站| 国产免费男女视频| 又粗又爽又猛毛片免费看| 国产成人福利小说| 国产精品嫩草影院av在线观看| 亚洲国产精品久久男人天堂| 日韩国内少妇激情av| 日韩欧美精品v在线| 国产精品精品国产色婷婷| 国产精品一二三区在线看| 人人妻人人看人人澡| 欧美又色又爽又黄视频| 激情 狠狠 欧美| 非洲黑人性xxxx精品又粗又长| 成人一区二区视频在线观看| 日本黄色视频三级网站网址| www.色视频.com| 国模一区二区三区四区视频| 天堂影院成人在线观看| 小说图片视频综合网站| 亚洲伊人久久精品综合 | 欧美丝袜亚洲另类| videos熟女内射| 色网站视频免费| 免费av毛片视频| 国产精品久久视频播放| 亚洲人成网站在线播| 99久久九九国产精品国产免费| 久久精品久久久久久噜噜老黄 | 18+在线观看网站| 日日啪夜夜撸| 国产国拍精品亚洲av在线观看| 欧美色视频一区免费| 国产淫语在线视频| 日韩av不卡免费在线播放| 色综合站精品国产| 水蜜桃什么品种好| 看黄色毛片网站| 亚洲精品久久久久久婷婷小说 | 51国产日韩欧美| 波多野结衣巨乳人妻| 亚洲精品日韩在线中文字幕| 亚洲精品国产成人久久av| 麻豆av噜噜一区二区三区| 欧美+日韩+精品| 精品久久久久久久久亚洲| 免费观看性生交大片5| 亚洲欧美精品综合久久99| 国产精品99久久久久久久久| 男女啪啪激烈高潮av片| 中文字幕熟女人妻在线| 婷婷色av中文字幕| 亚洲国产精品合色在线| 在线观看一区二区三区| 少妇猛男粗大的猛烈进出视频 | 99久久精品热视频| 最后的刺客免费高清国语| 国产一区有黄有色的免费视频 | 十八禁国产超污无遮挡网站| 精品人妻偷拍中文字幕| 一级毛片久久久久久久久女| 晚上一个人看的免费电影| 久久鲁丝午夜福利片| 亚洲美女视频黄频| 久久久久九九精品影院| 国产欧美日韩精品一区二区| 蜜臀久久99精品久久宅男| 久久99热这里只频精品6学生 | 国产精品乱码一区二三区的特点| 少妇丰满av| 日韩高清综合在线| 中文字幕久久专区| 日产精品乱码卡一卡2卡三| 嫩草影院新地址| 尤物成人国产欧美一区二区三区| av.在线天堂| 噜噜噜噜噜久久久久久91| 色视频www国产| 色综合色国产| 国产成人精品婷婷| 亚洲一区高清亚洲精品| 久久精品国产亚洲av涩爱| 草草在线视频免费看| 午夜福利视频1000在线观看| 夜夜爽夜夜爽视频| 国产av在哪里看| 日韩av在线大香蕉| 国产日韩欧美在线精品| 国产乱人视频| 免费观看性生交大片5| 午夜久久久久精精品| 男的添女的下面高潮视频| 免费看美女性在线毛片视频| 国产av一区在线观看免费| 三级国产精品欧美在线观看| or卡值多少钱| 蜜桃久久精品国产亚洲av| 大话2 男鬼变身卡| 听说在线观看完整版免费高清| 国产久久久一区二区三区| 在线a可以看的网站| 亚洲国产欧美人成| av又黄又爽大尺度在线免费看 | 18禁动态无遮挡网站| 国产熟女欧美一区二区| 欧美性猛交黑人性爽| 欧美成人精品欧美一级黄| 亚洲国产高清在线一区二区三| 欧美精品国产亚洲| 亚洲精品456在线播放app| 精品熟女少妇av免费看| 我要搜黄色片| 亚洲人成网站高清观看| 村上凉子中文字幕在线| 国产在视频线精品| 久久亚洲国产成人精品v| 草草在线视频免费看| 国产麻豆成人av免费视频| 99久久九九国产精品国产免费| 欧美精品国产亚洲| 久热久热在线精品观看| 身体一侧抽搐| 亚洲av熟女| 午夜精品一区二区三区免费看| 日产精品乱码卡一卡2卡三| 久久久久九九精品影院| 尾随美女入室| 日韩欧美三级三区| 一级av片app| 我的老师免费观看完整版| 久久精品久久精品一区二区三区| 亚洲人成网站高清观看| 亚洲国产精品sss在线观看| 免费观看的影片在线观看| 久久亚洲国产成人精品v| 最近最新中文字幕免费大全7| 久久精品久久久久久久性| 三级国产精品欧美在线观看| 观看美女的网站| 99热这里只有是精品50| av免费观看日本| 国产精品久久久久久精品电影小说 | 麻豆成人午夜福利视频| 国产精品乱码一区二三区的特点| 少妇的逼好多水| 九九爱精品视频在线观看| 亚洲精品乱码久久久久久按摩| 性色avwww在线观看| 国产成人一区二区在线| 99久久中文字幕三级久久日本| 成人鲁丝片一二三区免费| 亚洲欧美精品自产自拍| 91久久精品国产一区二区成人| 波多野结衣巨乳人妻| 免费不卡的大黄色大毛片视频在线观看 | 黑人高潮一二区| 亚洲av电影不卡..在线观看| 男女啪啪激烈高潮av片| 国产高清有码在线观看视频| 天堂av国产一区二区熟女人妻| 久久久成人免费电影| 欧美丝袜亚洲另类| 1000部很黄的大片| 天天躁日日操中文字幕| 亚洲无线观看免费| 少妇的逼好多水| 精品久久久久久久久久久久久| 亚洲自拍偷在线| 成人av在线播放网站| 欧美最新免费一区二区三区| 精品免费久久久久久久清纯| 成人欧美大片| 国产精品,欧美在线| 国产一级毛片在线| 亚洲精品乱码久久久久久按摩| 国产在视频线精品| 亚洲经典国产精华液单| 在线观看美女被高潮喷水网站| 日韩一区二区视频免费看| 最近最新中文字幕免费大全7| 欧美成人免费av一区二区三区| 神马国产精品三级电影在线观看| 好男人在线观看高清免费视频| 免费无遮挡裸体视频| 超碰av人人做人人爽久久| 国产在线男女| 亚洲av男天堂| 国产精品久久久久久久久免| 日日啪夜夜撸| 最近2019中文字幕mv第一页| 全区人妻精品视频| 免费搜索国产男女视频| 国产精品不卡视频一区二区| 变态另类丝袜制服| 久久人妻av系列| 日韩精品有码人妻一区| 又爽又黄无遮挡网站| 国产精品国产三级专区第一集| 久久久久网色| 亚洲欧美精品自产自拍| 一区二区三区免费毛片| 黄色配什么色好看| 一级毛片久久久久久久久女| 日韩一区二区视频免费看| 亚洲精品乱久久久久久| 中文字幕免费在线视频6| 亚洲成人中文字幕在线播放| 久久久久性生活片| 精品人妻偷拍中文字幕| 亚洲欧美日韩高清专用| 久久精品91蜜桃| 亚洲国产最新在线播放| 看十八女毛片水多多多| 午夜亚洲福利在线播放| 看片在线看免费视频| 精品午夜福利在线看| 天堂影院成人在线观看| 直男gayav资源| 亚洲精品自拍成人| 午夜激情欧美在线| 午夜福利网站1000一区二区三区| 1024手机看黄色片| 麻豆久久精品国产亚洲av| 亚洲精品国产av成人精品| 国产高清不卡午夜福利| 日日摸夜夜添夜夜爱| 搞女人的毛片| av在线天堂中文字幕| av专区在线播放| 久久精品影院6| 国产亚洲av片在线观看秒播厂 | 国产成人一区二区在线| 床上黄色一级片| 国产精品久久电影中文字幕| 1000部很黄的大片| 嫩草影院新地址| 国产在线一区二区三区精 | 午夜福利视频1000在线观看| 欧美xxxx性猛交bbbb| 99国产精品一区二区蜜桃av| 精品人妻一区二区三区麻豆| 五月伊人婷婷丁香| 亚洲一区高清亚洲精品| 国产精品福利在线免费观看| 国产视频内射| 欧美日本亚洲视频在线播放| 51国产日韩欧美| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区国产| 国产免费一级a男人的天堂| 免费观看的影片在线观看| 免费av不卡在线播放| 好男人视频免费观看在线| 亚洲色图av天堂| 久久韩国三级中文字幕| 亚洲经典国产精华液单| 久久6这里有精品| 国内精品一区二区在线观看| 69人妻影院| 人妻制服诱惑在线中文字幕| 69人妻影院| 伦理电影大哥的女人| 一夜夜www| 内射极品少妇av片p| 一夜夜www| 欧美激情久久久久久爽电影| 日韩av在线大香蕉| h日本视频在线播放| 一区二区三区免费毛片| 亚洲成人av在线免费| 亚洲精品影视一区二区三区av| 三级经典国产精品| 男女视频在线观看网站免费| av天堂中文字幕网| 一区二区三区乱码不卡18| 精品国产露脸久久av麻豆 | 国产极品天堂在线| 熟女人妻精品中文字幕| 日韩一本色道免费dvd| 国产黄色视频一区二区在线观看 | 卡戴珊不雅视频在线播放| 国产毛片a区久久久久| 男女视频在线观看网站免费| 国语自产精品视频在线第100页| 六月丁香七月| 插逼视频在线观看|