999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于ECMWF集合預(yù)報(bào)的淮河流域面雨量預(yù)報(bào)性能評(píng)估

2019-09-02 12:06:18葉金印
中低緯山地氣象 2019年4期

羅 靜,葉金印,王 喜,劉 靜

(1.安徽省氣象臺(tái),安徽 合肥 230031;2.中國(guó)氣象局氣象干部培訓(xùn)學(xué)院安徽分院,安徽 合肥 230031;3.江蘇省泰州市氣象局,江蘇 泰州 225300;4.安徽省蚌埠市氣象局,安徽 蚌埠 233040)

0 引言

面雨量在洪水模擬和預(yù)報(bào)中扮演極其重要的角色,無(wú)論是流域防汛抗洪,還是水利調(diào)度都離不開(kāi)準(zhǔn)確及時(shí)的面雨量資料。業(yè)務(wù)中常用的面雨量監(jiān)測(cè)產(chǎn)品是基于常規(guī)觀測(cè)資料計(jì)算出來(lái)的,而常規(guī)觀測(cè)站因站點(diǎn)分布不均勻、密度低易導(dǎo)致觀測(cè)數(shù)據(jù)空間代表性差,對(duì)面雨量計(jì)算結(jié)果有較大影響。近年來(lái),隨著雷達(dá)、衛(wèi)星資料的融入,無(wú)資料或缺資料流域亦可獲取高分辨率的降水場(chǎng)。CMORPH融合降水產(chǎn)品就是基于全國(guó)自動(dòng)站逐小時(shí)觀測(cè)數(shù)據(jù)和衛(wèi)星反演降水產(chǎn)品,利用概率密度匹配法、最優(yōu)插值法生成的融合降水產(chǎn)品[1]。該產(chǎn)品充分結(jié)合了地面觀測(cè)資料和衛(wèi)星反演降水各自的優(yōu)勢(shì),可提供高分辨率的實(shí)況降水格點(diǎn)場(chǎng),為開(kāi)展淮河流域精細(xì)化面雨量監(jiān)測(cè)提供一種新的思路。

近年來(lái),面雨量預(yù)報(bào)方法研究多側(cè)重于數(shù)理統(tǒng)計(jì)方法、天氣學(xué)診斷與數(shù)值預(yù)報(bào)相結(jié)合方法、多模式集成方法等[2-3],而對(duì)面雨量集合預(yù)報(bào)方法和檢驗(yàn)評(píng)估研究較少。ECMWF集合預(yù)報(bào)產(chǎn)品自2012年下半年開(kāi)始由國(guó)家氣象信息中心業(yè)務(wù)下發(fā),時(shí)空分辨率高、預(yù)報(bào)時(shí)效長(zhǎng),為淮河流域開(kāi)展面雨量集合預(yù)報(bào)業(yè)務(wù)提供了有利的條件。為更好地發(fā)揮ECMWF集合預(yù)報(bào)在淮河流域洪水風(fēng)險(xiǎn)管理和水資源調(diào)度工作中的決策支撐作用,對(duì)其預(yù)報(bào)性能進(jìn)行評(píng)估十分必要。

本文在對(duì)兩種面雨量監(jiān)測(cè)產(chǎn)品進(jìn)行對(duì)比分析的基礎(chǔ)上,分山區(qū)和平原地區(qū)進(jìn)一步檢驗(yàn)評(píng)估基于ECMWF集合預(yù)報(bào)產(chǎn)品的淮河流域面雨量預(yù)報(bào)性能。

1 流域子單元?jiǎng)澐趾唾Y料

1.1 淮河流域子單元?jiǎng)澐?/h3>

淮河流域橫跨山東、河南、安徽、江蘇四省,包含淮河和沂沭泗兩大水系?;春铀蛋l(fā)源于河南省桐柏山區(qū),全長(zhǎng)超過(guò)1 000 km,總落差200 m,其中上游落差高達(dá)178 m,中下游落差僅22 m,如遇洪水下泄不暢,極易導(dǎo)致洪澇災(zāi)害;沂沭泗水系發(fā)源于山東省沂蒙山區(qū),主要流經(jīng)山東、江蘇兩省,最終匯入江蘇境內(nèi)的洪澤湖。根據(jù)洪水匯流特征并結(jié)合水文氣象業(yè)務(wù)發(fā)展需求,淮河流域被劃分成15個(gè)子單元(如圖1)。各子單元的地理及氣候特征如表1所示,流域西部自北向南依次為伏牛山區(qū)(子單元5)、桐柏山區(qū)(子單元1、2、3)和大別山區(qū)(子單元4、7),東北部為沂蒙山區(qū)(子單元12、13),其余為廣闊的平原(子單元6、8、9、10、11、14、15)。

圖1 淮河流域15個(gè)子單元?jiǎng)澐諪ig.1 Division of 15 subunits in Huaihe River Basin

1.2 資料來(lái)源

CMORPH融合降水產(chǎn)品為國(guó)家氣象信息中心提供的時(shí)間分辨率為1 h、空間分辨率為0.1°×0.1°的格點(diǎn)降水資料,該資料始于2008年1月1日(北京時(shí),下同),本文所用資料時(shí)段為2008年1月1日—2017年12月31日。地面降水觀測(cè)資料采用國(guó)家氣象信息中心提供的淮河流域172個(gè)常規(guī)氣象站逐日降水(20—20時(shí))觀測(cè)整編資料,時(shí)段與CMORPH融合降水產(chǎn)品保持一致。ECMWF集合預(yù)報(bào)資料來(lái)源于國(guó)家氣象信息中心業(yè)務(wù)廣播系統(tǒng),所用資料時(shí)段為2013年1月1日—2017年12月31日,包含1個(gè)控制預(yù)報(bào)(確定性預(yù)報(bào))和50個(gè)擾動(dòng)預(yù)報(bào)(集合預(yù)報(bào)),預(yù)報(bào)時(shí)效360 h,0~72 h為3 h間隔,72~240 h為6 h間隔,240~360 h為12 h間隔,空間分辨率為0.5°×0.5°。

表1 流域各子單元地理、氣候特征Tab.1 The geographical and climatic characteristics of each subunit in the basin

1.3 面雨量計(jì)算方法

一次降水過(guò)程中整個(gè)流域的降水情況可用面雨量即流域面上的平均降水量來(lái)描述。算術(shù)平均法和泰森多邊形法是實(shí)際業(yè)務(wù)中最常用的兩種面雨量計(jì)算方法[4]。本文計(jì)算面雨量對(duì)于站點(diǎn)資料采用泰森多邊形法,對(duì)于格點(diǎn)資料采用算術(shù)平均法。

2 兩種雨量資料對(duì)比分析

分別計(jì)算基于兩種雨量資料的淮河流域所有子單元逐日面雨量,對(duì)比分析CMORPH融合降水資料在淮河流域面雨量監(jiān)測(cè)中的可用性。

2.1 相關(guān)分析

兩種面雨量計(jì)算結(jié)果分別用C(COMRPH融合產(chǎn)品)和G(地面觀測(cè)雨量資料)來(lái)表示,下同。對(duì)于所有樣本,C平均值(1.9 mm)小于G平均值(2.4 mm),整體表現(xiàn)為低估;從兩種計(jì)算結(jié)果的線性擬合(如圖2)來(lái)看,C和G存在顯著的線性相關(guān),總體相關(guān)系數(shù)高達(dá)0.956。分子單元統(tǒng)計(jì),子單元10相關(guān)最好為0.983,子單元12相關(guān)最差為0.916。進(jìn)一步統(tǒng)計(jì)相關(guān)系數(shù)與各子單元面積、平均海拔高度的關(guān)系(如圖3)可知,面積相對(duì)較大的平原地區(qū)相關(guān)較好,面積相對(duì)較小的山區(qū)相關(guān)較差。

圖2 兩種面雨量計(jì)算結(jié)果的一元線性回歸Fig.2 Linear regression of two kinds of areal precipitation

圖3 相關(guān)系數(shù)隨子流域面積、平均海拔高度的分布Fig.3 Distributions of correlation coefficient with area and average elevation of the subbasin

受特殊的地理氣候特征影響,淮河流域降水多集中在每年的5—9月。將所有樣本分汛期(5—9月)和非汛期(10月—次年4月)統(tǒng)計(jì)可知:對(duì)于整個(gè)流域,汛期C平均值為3.8 mm,G平均值為4.4 mm,相關(guān)系數(shù)為0.961;非汛期C平均值為1.0 mm,G平均值為1.3 mm,相關(guān)系數(shù)為0.938。對(duì)于各子單元(如表2),汛期相關(guān)總體好于非汛期;無(wú)論汛期還是非汛期,相關(guān)都是在山區(qū)較差,而在平原地區(qū)相對(duì)較好。

表2 流域各子單元相關(guān)系數(shù)Tab.2 Correlation coefficient of each subunit in the basin

2.2 量級(jí)偏差分析

根據(jù)我國(guó)24 h江河面雨量等級(jí)劃分標(biāo)準(zhǔn)[5],分小雨(0,6)、中雨[6,15)、大雨[15,30)、暴雨[30,60)、大暴雨[60,150)、大暴雨以上[150,+∞)6個(gè)量級(jí)對(duì)C和G進(jìn)行統(tǒng)計(jì)。如表3所示,除了無(wú)雨情況,對(duì)于各量級(jí)的G降水,C樣本集中在與G樣本同一量級(jí)和比其小一個(gè)量級(jí)的范圍內(nèi),這與2.1得出的C對(duì)G存在低估的結(jié)論相一致。各量級(jí)C與G一一對(duì)應(yīng)的比例分別為:小雨89.81%;中雨59.07%;大雨59.25%;暴雨57.08%;大暴雨54.4%;大暴雨以上0。分汛期和非汛期統(tǒng)計(jì)得出的結(jié)論與整體統(tǒng)計(jì)結(jié)果相一致(圖略)。

表3 兩種面雨量計(jì)算結(jié)果分量級(jí)統(tǒng)計(jì)Tab.3 Statistics of two kinds of areal precipitation on different grades

3 基于ECMWF的面雨量確定性預(yù)報(bào)性能評(píng)估

采用晴雨預(yù)報(bào)正確率、TS評(píng)分、漏報(bào)率和空?qǐng)?bào)率等指標(biāo)檢驗(yàn)評(píng)估基于ECMWF控制預(yù)報(bào)的淮河流域面雨量確定性預(yù)報(bào)性能。由于樣本數(shù)量有限,只分山區(qū)和平原地區(qū)分別對(duì)汛期和非汛期面雨量預(yù)報(bào)性能進(jìn)行評(píng)估分析。

3.1 正確率

正確率(PC)可以對(duì)“有”、“無(wú)”降水進(jìn)行預(yù)報(bào)檢驗(yàn),表達(dá)式為:

(1)

式中:NA為有雨預(yù)報(bào)正確次數(shù);NB為無(wú)雨空?qǐng)?bào)次數(shù);NC為有雨漏報(bào)次數(shù);ND為無(wú)雨預(yù)報(bào)正確次數(shù)。如圖4,無(wú)論山區(qū)還是平原地區(qū),汛期和非汛期24 h預(yù)報(bào)正確率均最高,接近80%,且正確率隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸降低。對(duì)于山區(qū),0~96 h正確率汛期高于非汛期,120~360 h反之;對(duì)于平原地區(qū),0~72 h正確率汛期高于非汛期,96~288 h非汛期高于汛期,312~360 h又是汛期高于非汛期;由此說(shuō)明,隨著預(yù)報(bào)時(shí)效的延長(zhǎng),對(duì)于晴雨預(yù)報(bào)的正確率汛期比非汛期衰減的速度更快,這種現(xiàn)象在山區(qū)表現(xiàn)更為明顯。

圖4 山區(qū)和平原不同預(yù)報(bào)時(shí)效降水預(yù)報(bào)的正確率Fig.4 Accuracy of precipitation forecast of mountain area and plain in different forecast periods

3.2 TS評(píng)分、漏報(bào)率、空?qǐng)?bào)率

TS評(píng)分(TS)、漏報(bào)率(PO)、空?qǐng)?bào)率 (FAR)可用來(lái)檢驗(yàn)?zāi)J綄?duì)不同量級(jí)降水的預(yù)報(bào)效果,表達(dá)式如下:

(2)

(3)

(4)

式(2)—(4)中k代表降水量級(jí);NAk、NBk、NCk分別為k量級(jí)降水預(yù)報(bào)正確(預(yù)報(bào)與實(shí)況量級(jí)相同)、空?qǐng)?bào)(預(yù)報(bào)量級(jí)大于實(shí)況量級(jí))及漏報(bào)(預(yù)報(bào)量級(jí)小于實(shí)況量級(jí))的樣本數(shù)。由于實(shí)況樣本和預(yù)報(bào)樣本均無(wú)大暴雨以上量級(jí)降水出現(xiàn),此節(jié)只分小雨、中雨、大雨、暴雨、大暴雨5個(gè)量級(jí)進(jìn)行統(tǒng)計(jì)。

結(jié)果如圖5、圖6所示:①對(duì)于汛期,無(wú)論是山區(qū)還是平原地區(qū),模式對(duì)各量級(jí)降水預(yù)報(bào)的性能隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸下降,且平原地區(qū)的預(yù)報(bào)性能優(yōu)于山區(qū);小雨、中雨量級(jí)降水預(yù)報(bào)的空?qǐng)?bào)率明顯大于漏報(bào)率,大雨及以上量級(jí)差異不明顯,說(shuō)明汛期模式對(duì)小雨、中雨量級(jí)降水預(yù)報(bào)的范圍和量級(jí)偏大,這與劉靜等[6]的研究結(jié)果較為一致。②對(duì)于非汛期,無(wú)論是山區(qū)還是平原地區(qū),對(duì)大暴雨基本無(wú)預(yù)報(bào)能力,其它量級(jí)降水預(yù)報(bào)性能隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸下降,且平原地區(qū)的預(yù)報(bào)性能優(yōu)于山區(qū);暴雨及其以下量級(jí)降水預(yù)報(bào)同樣存在高估現(xiàn)象。③無(wú)論是山區(qū)還是平原地區(qū),模式對(duì)小雨和極端降水的預(yù)報(bào)能力汛期好于非汛期,對(duì)中間量級(jí)降水預(yù)報(bào)則是非汛期表現(xiàn)更好。

圖5 汛期各量級(jí)降水預(yù)報(bào)的TS評(píng)分、漏報(bào)率、空?qǐng)?bào)率Fig.5 Treat Score, leak forecast quotiety and absent forecast quotiety of precipitation in flood season on different grades

圖6 非汛期各量級(jí)降水預(yù)報(bào)的TS評(píng)分、漏報(bào)率、空?qǐng)?bào)率Fig.6 Treat Score, leak forecast quotiety and absent forecast quotiety of precipitation in non flood season on different grades

4 基于ECMWF的面雨量集合預(yù)報(bào)性能評(píng)估

近年來(lái),集合預(yù)報(bào)技術(shù)廣泛應(yīng)用于天氣預(yù)報(bào)領(lǐng)域,將單一的確定性預(yù)報(bào)轉(zhuǎn)換為概率預(yù)報(bào),更好地服務(wù)于洪水預(yù)報(bào)預(yù)警和防災(zāi)減災(zāi)決策。集合預(yù)報(bào)性能評(píng)估的方法很多[7-10],本文采用CRPS(連續(xù)排序概率評(píng)分)及CRPS的兩個(gè)變體(RCRPS、CRPSS)來(lái)評(píng)估基于ECMWF擾動(dòng)預(yù)報(bào)的淮河流域面雨量集合預(yù)報(bào)性能。表達(dá)式為:

(5)

(6)

(7)

CRPS可用來(lái)評(píng)估具有一定概率密度分布的集合預(yù)報(bào)量與單一觀測(cè)值的一致程度,式(5)中x是預(yù)報(bào)量,xa是觀測(cè)值,P(x)是x的累積分布函數(shù),H(x-xa)是赫維賽德函數(shù)。一般意義上講,CRPS的大小與集合預(yù)報(bào)的性能成反比。但值得注意的是,CRPS的大小跟x的大小密切相關(guān),一個(gè)特定的地區(qū)或者季節(jié),一個(gè)較低的CRPS相比于其他區(qū)域或者季節(jié)并不一定等同于一個(gè)更好的預(yù)報(bào)結(jié)果。

可通過(guò)標(biāo)準(zhǔn)化處理來(lái)消除x對(duì)CRPS的影響,RCRPS就是其中一種標(biāo)準(zhǔn)化形式,式(6)中σa是某區(qū)域和某些研究樣本中所有xa的標(biāo)準(zhǔn)差。

CRPSS是CRPS的另一種標(biāo)準(zhǔn)化形式,式(7)中CRPSF代表預(yù)報(bào)值的CRPS評(píng)分,CRPSR代表同一預(yù)報(bào)變量的氣候預(yù)報(bào)值的CRPS評(píng)分。CRPSS不僅可以消除降水等級(jí)的影響,還能用來(lái)衡量集合預(yù)報(bào)系統(tǒng)相對(duì)于氣候預(yù)報(bào)的改進(jìn)程度。

4.1 面雨量集合預(yù)報(bào)整體性能分析

利用逐日累積降水量計(jì)算出的伏牛山區(qū)、桐柏山區(qū)、大別山區(qū)、沂蒙山區(qū)及平原地區(qū)5個(gè)汛期和4個(gè)非汛期平均的CRPSS(圖7)來(lái)評(píng)估淮河流域ECMWF面雨量集合預(yù)報(bào)相對(duì)于氣候預(yù)報(bào)的整體性能??梢?jiàn):無(wú)論是山區(qū)還是平原,汛期集合預(yù)報(bào)性能均隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸降低但仍優(yōu)于氣候均值預(yù)報(bào)。非汛期伏牛山區(qū)、大別山區(qū)和平原地區(qū)集合預(yù)報(bào)性能優(yōu)于氣候均值預(yù)報(bào),桐柏山區(qū)的優(yōu)勢(shì)只體現(xiàn)在0~96 h,沂蒙山區(qū)288 h以后CRPSS才大于0;沂蒙山區(qū)出現(xiàn)集合預(yù)報(bào)整體性能隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸增加的異?,F(xiàn)象,伏牛山區(qū)、桐柏山區(qū)、大別山區(qū)及平原地區(qū)240 h以后同樣有增加趨勢(shì)。

圖7 汛期(空心圓)和非汛期(實(shí)心圓)不同預(yù)報(bào)時(shí)效的平均CRPSSFig.7 Average CRPSS in flood season (hollow circle) and non flood season (solid circle) on different forecast periods

為了進(jìn)一步考察集合預(yù)報(bào)性能是否依賴于預(yù)報(bào)時(shí)效和預(yù)報(bào)累積時(shí)段,分別計(jì)算3 h、6 h、12 h、24 h所有可用預(yù)報(bào)時(shí)效的汛期和非汛期CRPS、RCRPS和CRPSS的平均值(圖8、圖9)。

4.2 不同預(yù)報(bào)時(shí)效和不同預(yù)報(bào)累積時(shí)段的集合預(yù)報(bào)性能分析

從CRPS來(lái)看:汛期任意預(yù)報(bào)累積時(shí)段的CRPS均隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸增加。3 h累積時(shí)段的CRPS隨著預(yù)報(bào)時(shí)效的延長(zhǎng)出現(xiàn)明顯的波動(dòng)狀態(tài),3~12 hCRPS逐漸減小,12~15 hCRPS迅速增大,15~24 hCRPS又開(kāi)始減小,如此往復(fù);隨著預(yù)報(bào)累積時(shí)段的增加,這種波動(dòng)現(xiàn)象越來(lái)越不明顯,12 h、24 h累積時(shí)段CRPS基本無(wú)波動(dòng)。非汛期CRPS的變化趨勢(shì)與汛期類似。各累積時(shí)段CRPS的值汛期明顯大于非汛期,且伏牛山區(qū)及沂蒙山區(qū)CRPS值明顯小于其他地區(qū),這與流域降水主要集中在汛期以及上述山區(qū)年平均降水量(表1)相對(duì)較小有關(guān)。

RCRPS是一個(gè)與預(yù)報(bào)變量量級(jí)無(wú)關(guān)的標(biāo)準(zhǔn)化評(píng)價(jià)指標(biāo),可用于對(duì)比評(píng)估不同地區(qū)、不同季節(jié)集合預(yù)報(bào)的性能。從RCRPS來(lái)看:RCRPS隨著預(yù)報(bào)時(shí)效及預(yù)報(bào)累積時(shí)段的變化趨勢(shì)與CRPS類似;汛期各地區(qū)預(yù)報(bào)性能差別相對(duì)較小,且山區(qū)整體預(yù)報(bào)性能好于平原地區(qū);非汛期各地區(qū)預(yù)報(bào)性能差別相對(duì)較大,平原地區(qū)預(yù)報(bào)性能表現(xiàn)更好,沂蒙山區(qū)出現(xiàn)預(yù)報(bào)性能隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸增加的異?,F(xiàn)象。

從CRPSS來(lái)看:汛期任意一個(gè)預(yù)報(bào)累計(jì)時(shí)段CRPSS的值基本都大于零,且隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸減小,說(shuō)明集合預(yù)報(bào)相對(duì)于氣候均值預(yù)報(bào)的優(yōu)勢(shì)逐漸降低,并且這種下降趨勢(shì)在0~96 h最為明顯;非汛期6、12、24 h累積時(shí)段沂蒙山區(qū)CRPSS隨著預(yù)報(bào)時(shí)效的增加逐漸增大,其它地區(qū)0~240 hCRPSS隨著預(yù)報(bào)時(shí)效的延長(zhǎng)逐漸減小,240 h以后同樣有增大趨勢(shì)。

圖8 汛期3 h、6 h、12 h、24 h預(yù)報(bào)累積時(shí)段CRPS、RCRPS、CRPSS平均值Fig.8 Average values of CRPS, RCRPS and CRPSS in cumulative periods of 3, 6, 12 and 24 hours in flood season

圖9 非汛期(b)3 h、6 h、12 h、24 h預(yù)報(bào)累積時(shí)段CRPS、RCRPS、CRPSS平均值Fig.9 Average values of CRPS, RCRPS and CRPSS in cumulative periods of 3, 6, 12 and 24 hours in non flood season

5 結(jié)語(yǔ)

基于CMORPH融合降水產(chǎn)品的淮河流域面雨量計(jì)算結(jié)果對(duì)基于地面降水觀測(cè)資料的面雨量計(jì)算結(jié)果存在低估現(xiàn)象,但二者總體相關(guān)系數(shù)高達(dá)0.956,COMRPH融合降水產(chǎn)品可應(yīng)用于淮河流域精細(xì)化面雨量監(jiān)測(cè)業(yè)務(wù)中。

基于ECMWF的確定性面雨量預(yù)報(bào)與實(shí)況量級(jí)整體相差不大,但存在小量級(jí)降水過(guò)度預(yù)報(bào)現(xiàn)象且對(duì)強(qiáng)降水量級(jí)預(yù)報(bào)不足。當(dāng)確定性預(yù)報(bào)有弱降水過(guò)程時(shí),要考慮降水量級(jí)預(yù)報(bào)是否偏大;預(yù)報(bào)有大范圍降水過(guò)程時(shí),特別是非汛期需要防范大范圍降水過(guò)程中局地強(qiáng)降水可能造成的山洪、局部?jī)?nèi)澇等災(zāi)害。

汛期平原地區(qū)和山區(qū)集合預(yù)報(bào)相對(duì)于氣候均值預(yù)報(bào)的優(yōu)勢(shì)都很明顯,且山區(qū)集合預(yù)報(bào)整體性能優(yōu)于平原地區(qū);非汛期除了沂蒙山區(qū),集合預(yù)報(bào)性能整體優(yōu)于氣候均值預(yù)報(bào),且平原地區(qū)表現(xiàn)更好。

面雨量預(yù)報(bào)可為水庫(kù)洪水預(yù)報(bào)、調(diào)度和管理提供技術(shù)支持手段,在流域防汛抗旱工作中具有十分重要的作用。防汛指揮部門(mén)在應(yīng)用面雨量預(yù)報(bào)時(shí),要注意長(zhǎng)、中、短期預(yù)報(bào)相結(jié)合,并根據(jù)最新降水預(yù)報(bào)及時(shí)調(diào)整防汛抗旱調(diào)度方案。

主站蜘蛛池模板: 国产女同自拍视频| 亚洲成人一区在线| 啪啪免费视频一区二区| 无码粉嫩虎白一线天在线观看| 超清无码一区二区三区| 国产视频a| 97青草最新免费精品视频| 国产午夜福利亚洲第一| 一区二区影院| 鲁鲁鲁爽爽爽在线视频观看 | 国产黄在线观看| 国产第一福利影院| 色视频久久| 国产资源免费观看| 久久亚洲中文字幕精品一区| 乱人伦99久久| 亚洲天堂久久新| 欧美乱妇高清无乱码免费| 亚洲精品久综合蜜| 免费人成视网站在线不卡| 91免费在线看| 国产欧美日韩资源在线观看| 国产精品无码AV片在线观看播放| 影音先锋丝袜制服| 欧美一级特黄aaaaaa在线看片| 亚洲国产av无码综合原创国产| 五月婷婷亚洲综合| 国产精品99在线观看| 99久久99视频| 女人18毛片一级毛片在线 | 久久久久无码精品| 91人妻在线视频| a色毛片免费视频| 色播五月婷婷| 国产日本欧美在线观看| www.亚洲国产| 日韩最新中文字幕| 精品丝袜美腿国产一区| 欧美日韩成人在线观看| 国产伦精品一区二区三区视频优播 | 国产福利小视频在线播放观看| 国产激情国语对白普通话| 欧美69视频在线| 鲁鲁鲁爽爽爽在线视频观看| 欧美a级完整在线观看| 国产精品亚洲片在线va| 久草视频中文| 99re热精品视频中文字幕不卡| 欧美午夜久久| 日韩精品无码一级毛片免费| 免费观看男人免费桶女人视频| 亚洲最猛黑人xxxx黑人猛交| 欧洲极品无码一区二区三区| 久久77777| 国产精品永久不卡免费视频| 婷婷99视频精品全部在线观看| 91精品在线视频观看| av在线无码浏览| 免费在线视频a| 国产精品手机视频一区二区| 无码国产偷倩在线播放老年人 | 欧美.成人.综合在线| 免费a级毛片18以上观看精品| 99久久精品久久久久久婷婷| 国产免费久久精品99re不卡 | 久久精品66| a国产精品| 99热这里只有免费国产精品 | 亚洲美女操| 毛片在线区| 亚洲天堂福利视频| 在线国产91| 国产精品亚欧美一区二区| 国产高清国内精品福利| 国产99在线观看| 日韩麻豆小视频| 毛片手机在线看| 中文一区二区视频| 最新国产成人剧情在线播放| 91小视频在线播放| 91免费观看视频| a级毛片在线免费观看|