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

    2000年以來我國長江中游區(qū)域暴雨研究進展

    2014-03-02 01:48:29崔春光林春澤王曉芳李俊彭濤張文剛
    關鍵詞:系統(tǒng)

    崔春光 林春澤 王曉芳 李俊 彭濤 張文剛

    (中國氣象局武漢暴雨研究所 暴雨監(jiān)測預警湖北省重點實驗室,武漢 430074)

    2000年以來我國長江中游區(qū)域暴雨研究進展

    崔春光 林春澤 王曉芳 李俊 彭濤 張文剛

    (中國氣象局武漢暴雨研究所 暴雨監(jiān)測預警湖北省重點實驗室,武漢 430074)

    長江中游地區(qū)是暴雨多發(fā)區(qū)和重災區(qū),深入認識該區(qū)域暴雨形成的中、小尺度天氣系統(tǒng)、分析暴雨發(fā)生機理與發(fā)展暴雨、洪水預報技術,對于我國防洪減災具有重要意義。針對長江中游區(qū)域暴雨從觀測到資料分析、機理研究,再到數值預報與水文應用研究,重點從暴雨外場觀測試驗及資料反演和中尺度資料融合技術、長江中游暴雨MCS(中尺度對流系統(tǒng))觀測分析和理論研究、中尺度暴雨數值模式研究和長江中游水文氣象耦合的暴雨洪水預報系統(tǒng)四個方面介紹了2000年以來我國科研業(yè)務單位的主要研究進展,并對制約長江中游區(qū)域暴雨研究深入與預報水平提高的主要問題進行了簡要討論,對未來的研究趨勢進行了展望。

    資料融合,中尺度對流系統(tǒng),水文氣象耦合

    1 引言

    長江中游是我國暴雨頻發(fā)區(qū)域,受中高緯西風帶環(huán)流系統(tǒng)、亞洲季風系統(tǒng)以及高原東移低值系統(tǒng)的共同作用,每年5到9月,暴雨過程頻繁,且持續(xù)時間長,時空分布不均,既有梅雨期的大范圍持續(xù)性降水,也有盛夏局地突發(fā)性強降水,所造成的災害重、損失大。對長江中游影響最大主要降水集中在每年6月中旬到7月中旬,即長江梅雨期,該時期長江流域是極地冷空氣和熱帶暖濕空氣頻繁交綏的地帶,梅雨鋒上的中尺度天氣系統(tǒng)活躍,它們的活動不僅維持了梅雨期持續(xù)性的降水,而且給暴雨的產生創(chuàng)造了十分有利的條件[1]。暴雨是由中小尺度天氣系統(tǒng)直接造成的,但它的發(fā)生發(fā)展是多尺度環(huán)流和天氣系統(tǒng)相互作用的結果,因而暴雨形成是一個十分復雜的問題[2]。長江中游區(qū)域長江、漢江東西橫貫,洞庭湖、鄱陽湖水系密布,大小水庫眾多,是我國防汛的關鍵區(qū)域。因此長江中游區(qū)域暴雨尤其是梅雨鋒暴雨的研究受到國內外氣象學者和我國科研、業(yè)務機構的廣泛關注和高度重視,在國家和地方各類科研項目的推動下,研究工作取得了明顯進步。本文對2000年以來我國長江中游區(qū)域暴雨研究進展做一些概括性歸納和總結,重點從暴雨外場觀測試驗及資料反演和中尺度資料融合技術、長江中游暴雨MCS觀測分析和理論研究、中尺

    度暴雨數值模式研究和長江中游水文氣象耦合的暴雨洪水預報系統(tǒng)四個方面進行介紹,并對制約長江中游區(qū)域暴雨研究深入與預報水平提高的主要因素進行討論,對未來開展的重點工作與研究趨勢進行了展望。

    2 暴雨外場觀測試驗及資料反演和中尺度資料融合技術

    2.1 暴雨外場觀測試驗

    由于引發(fā)暴雨的中尺度對流系統(tǒng)具有一定的頻發(fā)性、突發(fā)性以及較小的時空分辨率,使得常規(guī)氣象觀測難以捕獲,因此在重點區(qū)域開展中尺度暴雨研究的外場觀測試驗以獲取中尺度暴雨三維結構信息對于深入研究暴雨形成機理與提高暴雨預報水平尤為重要。1998年長江流域大洪水以來,國內開展了多次暴雨野外觀測試驗。如國家973“我國重大天氣災害形成機理與預測理論研究”2001—2002年汛期開展的長江中下游暴雨外場觀測試驗以梅雨鋒內中尺度暴雨為主要觀測對象,獲取了大量暴雨過程實時加密資料, 極大地推動了長江梅雨鋒暴雨中尺度研究和預報[3]。國家973“我國南方致洪暴雨監(jiān)測與預測理論和方法研究”于2008 年、2009 年5 月1 日—7 月20 日實施了中國南方暴雨野外科學試驗[4-5],該觀測試驗是迄今為止國內在單一項目支持下規(guī)模最大、裝備最為先進的一次暴雨外場觀測試驗,對我國汛期雨帶的主要影響區(qū)域及其北抬的完整過程實施加密觀測,涉及四個觀測區(qū)和十三省兩市,幾乎覆蓋整個長江以南地區(qū)。投入的設備有55部雷達,48部探空,7000個地面觀測站,以及特種設備如C波段和X波段雙偏振多普勒雷達與8毫米波云雷達,組織的多雷達組網觀測獲取的資料為分析研究暴雨云內水凝物分布、發(fā)展對流云微物理過程提供了良好基礎,并首次在外場觀測試驗中實施了機載下投式探空觀測,獲取了近40條氣象廓線的觀測資料。

    在前述兩個973項目支持下開展的暴雨外場觀測試驗,是國內無論從觀測、分析到應用都較為先進的β中尺度暴雨加密觀測試驗,尤其是在多遙感途徑聯合觀測、新型裝備研制與中試、高時空分辨率的中尺度氣象分析場生成等方面進行了有力探索,武漢暴雨研究所作為主要參與單位之一,在此過程中逐步打下了堅實的設備與技術基礎。2002年至今,武漢暴雨研究所依托氣象業(yè)務觀測網,在氣象災害頻發(fā)的長江中游地區(qū)建設具備高密度水汽監(jiān)測、高空觀測、云觀測、邊界層觀測、定量降水估算校準觀測等功能的β中尺度暴雨監(jiān)測外場試驗基地(見圖1),目前已建成了咸寧、武漢、荊州三個探測基地,與黃石、孝感、荊門、仙桃、神農架、黃岡等多個觀測點,在長江中游初步形成了點面相結合的觀測布局?;靥綔y設備包括不同類型和功能的車載雷達、邊界層風廓線雷達、GPS移動探空系統(tǒng)、微波輻射計、地基GPS/ MET站、激光雨滴譜儀、神農架高山梯度氣象觀測網等。此外,各個基地還設有業(yè)務天氣雷達、業(yè)務高空探空站、自動氣象站、雨量站等多種常規(guī)觀測設備,可對長江中游重要暴雨天氣系統(tǒng)的風場、水汽場、降水場、云和降水粒子相態(tài)分布進行固定的實時綜合觀測。2010—2013年汛期,武漢暴雨所在基地內組織多次長江中游梅雨鋒暴雨多遙感途徑聯合觀測試驗,成

    功捕捉了多次梅雨期重要天氣過程,積累了大量多途徑高精度同步觀測資料,有效推動了暴雨監(jiān)測預警技術、暴雨機理、暴雨預報等方面研究工作的開展。

    圖1 長江中游暴雨監(jiān)測外場試驗基地區(qū)域示意圖

    2.2 資料反演技術

    氣象衛(wèi)星和雷達是目前遙感探測暴雨的兩種主要手段,在暴雨的監(jiān)測預警中發(fā)揮了重要作用,但其探測得到的基本信息遠不能滿足暴雨分析的需求,為了得到更多有效的信息,資料反演技術應運而生??茖W有效的反演技術能夠準確地解譯遙感資料中包含的暴雨相關信息,對分析暴雨的發(fā)生及發(fā)展機制、改進和提高暴雨的預報水平有非常重要的意義。

    目前對于暴雨的分析預報具有重要指示作用的參數包括風場、溫度場、云特性等,能夠詳細描述暴雨云團的動力及熱力形成條件、云團特征等信息?;趩味嗥绽绽走_和靜止衛(wèi)星云圖的風場反演已有數十年歷史,其反演技術已比較成熟,隨著遙感技術和數學方法的不斷發(fā)展,數值微分方法和正則化方法被用來提高云導風的反演精度[6],取得了理想的效果?;诶走_風場反演的基本算法,國內外許多學者對其進行了改進[7-9],一定程度上提高了反演精度,但距離業(yè)務使用還有一定差距,因此變分法和多雷達風場反演方法被提出。通過比較驗證發(fā)現,變分法對于風場反演具有較高精度,四維變分法能夠反演得到低層風場信息[10],劉黎平[11]采用三維變分法實現了雙多基地多普勒雷達系統(tǒng)三維風場的有效反演。利用雙多普勒雷達同步觀測,得到兩個方向的徑向速度資料,再與其他約束條件結合就可以得到風場的三維結構[12],羅昌榮等[13]利用該方法有效地分析了臺風外圍強帶狀回波風場結構特征。彭亮等[14]將小粒子示蹤法用于3mm多普勒云雷達云內空氣垂直速度的反演。Rosenfeld提出的衛(wèi)星微物理反演分析方法,對于云的微物理特征的獲取提供了非常有效的方法,該方法已在NOAA[15-16]、FY[17-19]、MODIS[20]等衛(wèi)星資料分析中應用,得到了云的微物理及垂直結構特征。根據Z-qr關系及定常暖云方案,利用雷達反射率因子可以對云的微物理量進行反演[21]。劉黎平等[22]利用地基Ka波段雷達探測的回波強度、粒子下落速度和速度譜寬數據,在忽略空氣上升速度和湍流對雷達觀測的速度譜寬的影響前提下,對云液態(tài)水含量及滴譜參數進行了反演,并與機載探測數據進行了對比,結果顯示兩者結果比較一致。大氣溫度及濕度廓線的探測主要利用星載大氣探測儀,其中ATOVS是具有代表性的一種,采用一維變分技術,同時對輻射傳輸方程做擾動處理,對大氣廓線的反演效果較好[23-24]。最小方差算法被用于NOAA/ AMSU資料海面大氣溫度垂直廓線的反演,并明顯提高了對流層大部分區(qū)域的反演精度[25],而利用特征向量統(tǒng)計反演算法由MODIS資料同樣可以反演得到大氣溫濕廓線[26]。王曦等[27]利用FY-3A微波輻射計亮溫數據,通過建立適用于微波波段的物理迭代算法,反演得到了晴空海域的大氣濕度廓線。王祥等[28]利用FY-3A衛(wèi)星近紅外波段資料,分別利用雙通道比值法與三通道比值法反演了大氣水汽含量,并與實測數據進行對比,發(fā)現三通道比值法精度更高。經過長時間的積累和不斷改進,資料反演技術較之前取得了明顯提高,但其中仍存在一些問題,例如單雷達風場反演算法主要以科研為主,目前還無法達到普適性,難以業(yè)務化應用,衛(wèi)星資料反演主要針對海洋區(qū)域,對受復雜地形影響的陸地區(qū)域難以開展有效的反演等。隨著計算機水平及傳感器探測精度的提高,資料反演技術必將會越來越高效,探測資料中更多的有效信息將會被解譯,為暴雨的分析研究提供幫助。

    2.3 中尺度資料融合技術

    隨著我國中尺度觀測系統(tǒng)的建設與完善,衛(wèi)星、雷達、微波輻射計、GPS、風廓線儀等非常規(guī)觀測手段日益豐富,所獲取的高時空分辨率的探測資料為暴雨中尺度系統(tǒng)的深入分析和預報能力的提高提供了可靠保證,然而這些非常規(guī)資料各自獨立,且能提供的直接產品大多采用形態(tài)學(如回波形狀、云圖形態(tài))來做定性分析,對中尺度暴雨系統(tǒng)結構特征分析十分有限[5]。中尺度資料融合技術可以將分布在全國的地面站觀測要素,連同來自衛(wèi)星、雷達、微波輻射計、GPS、風廓線儀等的各種資料融合成整體,形成包含豐富中尺度信息的中尺度分析場,并可為中尺度數值模式提供高質量的初始場。這為解決多源非常規(guī)觀測資料在暴雨中尺度天氣分析與預報中有效應用的難題提供了一個有效手段。

    20世紀90年代以來,國際上普遍采用資料同化技術綜合處理各類氣象觀測數據來重建長期而連續(xù)的格點再分析資料,并逐漸使用變分同化取代了最優(yōu)插值技術,以四維變分同化為基礎的集合—變分同化和混合資料同化技術更是成為了未來資料同化的主流發(fā)展方向[29]。近10年,我國資料同化技術取得了長足進步,但依然缺少系統(tǒng)的中尺度資料融合分析工具。武漢暴雨研究所于2005年與美國NOAA下屬的地球系統(tǒng)實驗室合作引進了局地分析系統(tǒng)(LAPS),并成功進行了本地化開發(fā),完成了LAPS 系統(tǒng)對我國探空資料、地面資料、多普勒雷達資料等的融合同化[30]。LAPS系統(tǒng)的基礎算法是在背景場基礎上采用

    距離權重插值得到網格點值,然后對氣溫、氣壓、風的關系采用三維變分進行約束,對垂直水汽分布采用一維變分進行約束。早在1996年LAPS系統(tǒng)已為夏季亞特蘭大奧運會提供服務[31]。GOES衛(wèi)星輻射資料、GPS水汽等資料同化進入LAPS系統(tǒng)后提高了LAPS系統(tǒng)的水汽分析能力[32]。武漢暴雨所結合暴雨過程對雷達資料、GPS水汽資料的應用開展了相關研究[33-34]。國家973“我國南方致洪暴雨監(jiān)測與預測理論和方法研究”利用由武漢暴雨所引進的LAPS 分析系統(tǒng),融合了2008—2009年5月1日—7月20日我國南方暴雨野外科學試驗期間所獲取的華南、華中、江淮和長三角觀測區(qū)加密觀測探空資料、多普勒雷達資料和地面資料,生成了水平分辨為5km,時間分辨為3h的LAPS中尺度分析資料,批量檢驗結果表明,LAPS 生成的中尺度分析系統(tǒng)能抓住80%以上引發(fā)暴雨過程的中尺度影響系統(tǒng),尤其對α中尺度或β中尺度系統(tǒng)具有較好的中尺度再現能力,對于觀測區(qū)域內由β中尺度低渦東移及切變線引成的暴雨過程,LAPS中尺度分析場均能很好地再現其中尺度影響系統(tǒng)的發(fā)展演變過程,LAPS中尺度分析場的位勢高度、溫度、相對濕度、風向和風速的分析誤差均在觀測誤差范圍內[5]。目前在武漢暴雨所業(yè)務運行的LAPS系統(tǒng),融合了長江流域30多部雷達資料及衛(wèi)星、GPS、自動氣象站等資料,實現大范圍雷達資料的綜合應用,提供3h更新5km分辨率垂直21層近40種物理量的中尺度氣象分析場(見圖2)。NOAA新開發(fā)的變分LAPS系統(tǒng)(VARLAPS),是基于多重網格法上的變分同化系統(tǒng),具有傳統(tǒng)LAPS系統(tǒng)計算速度快等優(yōu)點,并可靈活地在不同的網格上使用不同的約束條件以達到最優(yōu)的效果[35]。VAR-LAPS利用多尺度變分技術解決了舊版本在客觀分析中平衡與約束不能同時應用、遙感觀測數據沒有最佳融合、交叉變量協(xié)方差不能直接應用等問題。Li等[36-37]利用該系統(tǒng)開展多普勒雷達資料同化對中國海域海溫的預報。武漢暴雨所正在開展VARLAPS對新型探測資料如風廓線、GPS水汽等的融合同化試驗研究。VAR-LAPS作為中尺度資料融合支撐系統(tǒng)將在建立更優(yōu)暴雨中尺度分析場數據集與提高數值預報模式初始場質量等方面有著廣泛的應用前景。

    圖2 武漢暴雨研究所LAPS系統(tǒng)流程圖

    3 長江中游暴雨MCS觀測分析和理論研究

    3.1 長江中游MCS的觀測特征

    產生暴雨的MCS雷達回波上呈現了多種多樣的組織結構形態(tài),王曉芳等[38]將MCS分為線狀和非線狀,其中線狀又為8個確定的類型(TS、LS、PS、BL四類MCS為移動型,而BB、TL/AS、EL、LL四類MCS是準靜止的),而非線狀的MCS組織結構特征不顯著,主要是大片回波中嵌有或多或少的無組織強對流回波。Zheng等[39]通過研究我國中東部地區(qū)的MCS,考慮到線狀弓形回波對流和層狀對流易造成大風等劇烈天氣,將弓狀回波(BE)和無層狀降水的隨流線(NS)單獨分類。颮線是產生劇烈天氣的中尺度對流系統(tǒng),最近孟智勇[40]對我國東部的颮線特征也進行了研究。

    MCS常形成在有利的天氣尺度條件下,如副熱帶高壓、熱帶氣旋、中高緯度槽前、冷渦、槽后等, Meng等[40]統(tǒng)計發(fā)現有三分之一的颮線形成于短波槽天氣形勢下。當副熱帶高壓、中高緯度冷空氣活動、高原中尺度系統(tǒng)東移、南海季風涌等多尺度天氣系統(tǒng)的協(xié)同作用在這些系統(tǒng)處于活躍期時,梅雨鋒上出現暴雨,但并不是整個梅雨鋒帶上普降強降水,降水還存在不均勻性,這說明降水強弱還與局地的氣象條件特別是垂直風切變密切相關[41-43]。根據實際個例分析和數值模擬,總結幾類線狀MCS的結構模型,并探討了它們發(fā)生發(fā)展的可能機制[44-45]。

    長江流域線狀MCS多形成在長江兩岸附近,在宜昌至武漢沿江兩岸、湖南東北部是MCS的高發(fā)地,其次為重慶沿江北部、重慶南部至鄂西、大別山兩側、湖南和江西交界山區(qū),沿地形的低洼地和河谷平地是MCS的多發(fā)地。MCS移動路徑為東、東偏北、東偏南、南四種,這與環(huán)境場的引導氣流有關,MCS日變化活動特征表現為多鋒型[38,46],MCS生命史較長,水平尺度大,日落后減弱,午夜再度發(fā)展,即濕對流。陳炯等[47]分析中國暖季短時強降水特征,發(fā)現它與MCS的日變化特征一致。祁秀香等[48]統(tǒng)計得到2007年夏季我國大陸深對流活動主要集中在4個區(qū)域,除華南沿海地區(qū)外,其余三個:青藏高原,云貴高原東部及四川、重慶,江淮流域等都是中緯度帶區(qū),江淮流域深對流日變化也呈多峰型且明顯東傳。

    3.2 MCS的理論研究

    胡伯威[49]用一個包含簡單積云對流參數化的線性模式討論在長江流域相當正壓性質的梅雨鋒上深厚慣性重力波的發(fā)生、傳播和頻散性質,發(fā)現與積云對流加熱場有最密切關系的低層濕度條件及其水平分布不均勻性質都對擾動的波速和穩(wěn)定度有重要影響,最有利出現不穩(wěn)定發(fā)展的地方是在“低層濕度鋒區(qū)”南界附近,因此這里是中尺度對流系統(tǒng)(MCS)的活躍地帶。關于頻散性質的分析表明,能允許群速度大于相速度的低層濕度條件及其水平分布的特征值范圍比較寬。而在上述穩(wěn)定度接近中性的條件下最有利于出現沿“濕度鋒區(qū)”軸向的上、下游效應。這個結果也許能很好地解釋:在梅雨鋒上,恰恰是當一個旺盛的MCS達到頂點而開始趨向衰減時,在它的上、下游很快發(fā)展起新的MCS,新、舊MCS一起形成“云團波串”這樣一個重要而有趣的觀測事實。同樣在虛擬高度坐標系中,水平風切變強度不均勻分布對CISK慣性重力波不穩(wěn)定有重要作用。在一般干的層結大氣中,實際可能出現再強的水平風切變的影響也難以使慣性重力波變得不穩(wěn)定;只有在積云對流潛熱參與,原為弱穩(wěn)定條件下,水平風切變強度不均勻能促使低空急流北側不穩(wěn)定擾動的發(fā)生發(fā)展。而水平風切變強度不均勻對不穩(wěn)定貢獻最大的區(qū)域是梅雨鋒南側的急流軸附近[50]。

    圖3 濕度鋒上云團波串發(fā)生發(fā)展示意圖[51]

    3.3 中尺度地形對MCS的影響

    長江中游地區(qū)多中尺度地形,長江中游暴雨尤其是(特)大暴雨的發(fā)生,與這些地形密不可分。在估計山脈迎風坡降水方面,早期Smith[57]提出的簡單線性模型降水率的大小由水汽混合比、正交于山脈的風速大小和地形梯度決定。之后,考慮了云微物理過程中云水—雨水的轉化時間、雨水下落時間對地形降水的影響[58-59]。當濕大氣地形費勞德數Fr>>1或者Fw>>1時,氣流爬過山脈產生的水平擾動小,線性理論對地形降水能做出較好的解釋;當Fr<<1或者Fw<<1時,氣流發(fā)生水平擾動大,線性理論很難再適用[60]。地形阻滯效應限制了該模型的應用[61-62]。

    趙玉春等[63]根據長江中下游中小尺度地形特征和梅雨鋒上MCS發(fā)展的典型探空曲線,設計理想數值試驗來研究中尺度地形對流降水模態(tài),探討了不同地形

    高度、尺度、形狀(水平形態(tài)比)、環(huán)境氣流以及旋轉效應對降水模態(tài)和分布的影響,發(fā)現地形重力波破碎觸發(fā)對流是地形對流降水的主要機制之一,地形重力波的形成、特征和傳播受到環(huán)境條件制約,也受到地形自身參數的影響。單一的地形動力參數不能完全決定地形對流降水的模態(tài)。同時,在天氣系統(tǒng)發(fā)展的不同階段,地形相對于對流暴雨系統(tǒng)的位置不同,地形對暴雨對流系統(tǒng)的影響存在著差異[64]。廖移山等[65]通過對一次特大暴雨過程的天氣學機理分析,探討了冷空氣與復雜地形的交叉阻擋對暖濕氣流的強迫抬升作用。地形對暴雨對流系統(tǒng)的影響還值得進一步研究。

    4 中尺度暴雨數值模式研究

    4.1 中尺度暴雨數值模式發(fā)展概述

    武漢暴雨所自20世紀80年代開始發(fā)展降水模式的研究和業(yè)務工作,先后發(fā)展了華中暴雨數值預報準業(yè)務系統(tǒng)(簡稱CRONWS)、暴雨中尺度分析預報系統(tǒng)(簡稱MAPS)等,1999年開始,與中科院大氣所合作,在國家“973”項目的支持下,基于LASG“REM”模式,發(fā)展了新一代暴雨數值模式系統(tǒng)AREM[66],采用η坐標,抓住了東亞區(qū)域的地形特點[67],充分考慮暴雨模式中水汽等要素的強梯度分布特征,設計了兩步正定保形水汽平流方案[68],設計了獨特的E網格“半格距差分方案”,保證了模式運行效率、精度和穩(wěn)定度等特點[69]。2007年,在行業(yè)專項的支持下,針對制約暴雨模式發(fā)展的多源資料利用、資料同化、物理過程等關鍵技術開展了研究,在原有簡單初值系統(tǒng)(客觀分析)的基礎上,發(fā)展了先進的多源資料融合技術、四維變分同化技術等,發(fā)展了冷云微物理方案,改進和完善了模式邊界層、陸面物理過程等方案,同時建立了快速循環(huán)預報系統(tǒng)和集合預報系統(tǒng)等,暴雨的定量預報水平得到了顯著提高。

    4.2 資料同化技術

    李紅莉等[30,33-34]開展了多普勒雷達、風云二號衛(wèi)星云圖、GPS水汽、風廓線雷達、微波輻射計、自動站、地面降水等多源觀測資料的同化應用研究,并實現了上述多源觀測資料在同一個平臺進行同化,并開展了相應的數值預報試驗。結果表明,融合同化雷達徑向風資料能顯著改善中尺度風場結構,GPS水汽資料能調整分析區(qū)的濕度場,風廓線雷達能調整風場結構、微波輻射計主要通過反演相對濕度廓線調整模式的濕度場。地面資料是資料同化中的難點,在解決稠密地面資料在模式中的應用方面,武漢暴雨研究所發(fā)展了一套降水資料的一維變分同化技術,通過反演大氣溫、濕廓線調整模式的初始場[70]。以上多源資料的融合均明顯增加了模式初始場的中尺度信息,進而提高了模式對強降水的預報能力。

    此外,武漢暴雨研究所與地球流體力學數值模擬國家重點實驗室合作,發(fā)展了一種新的快速有效的四維變分資料同化方法——基于歷史擬合投影技術的四維變分同化方法(HSP-4Dvar),該方法解決了國際上經典四維變分中B矩陣低估和奇異處理困難的難題,同時克服了經典四維變分依靠伴隨模式和計算量巨大問題,不僅提高了四維變分同化的精度,同時計算時效節(jié)省20倍[71-72]。

    近年來,適應性觀測成為數值預報領域發(fā)展的熱點,武漢暴雨研究所利用WRF模式的伴隨同化系統(tǒng),開展了短期預報誤差對觀測資料的敏感性研究,該研究可為提高資料同化效率、優(yōu)化觀測網絡以及開展目標觀測試驗提供科學依據。

    4.3 中尺度模式關鍵物理過程的改進

    在物理過程改進中,重點考慮了那些對中尺度暴雨過程短期預報有重要影響的物理過程,主要包括顯、隱式網格尺度云和降水過程、積云對流過程、邊界層過程以及地-氣交換過程等。在新發(fā)展的冷云微物理過程中,考慮了包括水汽、云水、云冰、雨水、雪5種不同相態(tài)水物質,以及各種凝結水物質之間凝結/蒸發(fā)、凝華/升華、凍結/融化、集聚和碰并等25種復雜的云微物理過程,更好地描述了暴雨成云至雨的變化過程。通過網格變換和變量變換,解決了MYJ方案與AREM模式的匹配問題,通過對邊界層觀測資料的研究,用更能合理反映地表濕度的變化的通量方案替代了蒸發(fā)方案[73]。在引進Wang云微物理過程的基礎上,修改了冰核的數密度以及雹的截獲參數,并對時間尺度的公式進行了修改,改進了對臺風和熱帶、副熱帶降水的描述能力。此外,通過陸面物理過程引入,改進了模式對下墊面強迫(包括感熱,潛熱)的估計,從而改進模式預報效果[74-75]。

    4.4 快速更新循環(huán)預報技術

    在LAPS云分析技術的基礎上,武漢暴雨研究所發(fā)展了一套3h快速更新同化預報系統(tǒng)(AREMRUC)[76]。該系統(tǒng)于2009年5月下旬投入準業(yè)務運行,AREM-RUC每3h更新同化一次,實現了多普勒雷達、加密自動站、GPS水汽等局地高頻資料在模式中的應用,這對及時捕捉大氣變化和中小尺度天氣系統(tǒng)、提高降水預報準確率有重要意義;同時熱啟動技術的加入[77],有效地抑制了模式的SPIN-UP現象,使模式在短時預報上的應用成為可能。AREM-RUC成

    功預報出湖北汛期多次強降水過程,TS評分表明,該系統(tǒng)明顯高于原有的模式系統(tǒng)[78],特別是在大雨、暴雨、大暴雨等強降水等級。

    4.5 暴雨集合預報技術

    暴雨對模式的初值和物理過程都非常敏感,初值擾動的結構比擾動的大小更重要[79],精細的小尺度擾動可以改進對強降水的集合預報[80]。在對物理過程擾動的研究中,初步的結果表明隨機物理過程方法更適合于接近正壓的暖濕環(huán)境中的梅雨鋒暴雨。此外還開展了基于集合敏感性的目標觀測以改進暴雨預報的研究。圖4給出了針對一次西南低渦東移過程的集合敏感性試驗,通過改善敏感區(qū)中敏感因子(水汽和風場)的觀測可以改進對西南低渦的路徑預報,結果表明,集合預報在認識暴雨機制、開展目標觀測方面的潛在用途[81]。基于Breeding方法,2009年建立了11個成員的中尺度集合預報系統(tǒng),提供降水概率等級和預報產品,檢驗結果表明,由于能提供更多預報信息,集合預報系統(tǒng)優(yōu)于單一的確定性預報,能有效提高對暴雨預報決策的科學性[82]。

    4.6 暴雨數值預報業(yè)務支撐系統(tǒng)

    基于以上研究,武漢暴雨研究所建立起了相對完善的區(qū)域數值預報體系,其中包括資料融合系統(tǒng)、確定性預報系統(tǒng)、快速更新循環(huán)和集合預報系統(tǒng)等。此外,針對精細化預報的需求,建立了基于WRF3.4模式的精細化預報系統(tǒng),采用WRF3Dvar,目前以上系統(tǒng)均已實現業(yè)務化運行,在預報業(yè)務中發(fā)揮著技術支撐作用。

    5 長江中游水文氣象耦合的暴雨洪水預報系統(tǒng)

    長江中游河流、水庫、湖泊眾多,汛期防汛形勢緊張,迫切需求水文氣象耦合的暴雨洪水預報系統(tǒng)來開展預報與防汛決策服務。提高洪水預報精度和延長預見期,必須從提高降雨估算精度開始,結合降雨預報,采取流域降雨徑流模型和河道洪水演算模型的途徑來實現。早在30多年前,陳金榮等[83]從長江流域防洪的豐富實踐中對水文氣象預報已做了精辟而深刻的論述,指出氣象預報同水文預報的結合為防洪和水庫調度等工作,帶來巨大的優(yōu)越性和必要性,是水文預報發(fā)展的必然趨勢。文中闡述的基本思想在水文氣象預報技術已有許多新的發(fā)展和進步的今天,仍然具有指導意義。當前,氣象模式和水文模型分別是氣象和水文領域里最先進的工具,水文模型與氣象模式耦合開展洪水預報必然成為未來發(fā)展的方向。

    5.1 降水降尺度技術

    洪水預報中充分有效地利用天氣預報及氣象信息,其首要任務就是解決降水信息場與水文模型之間的時間和空間尺度上匹配的問題,縮小兩者的尺度差異,尋找水文氣象結合的契機。為了有效地提高定量降水估算預報產品在水文預報中的應用,國內外學者從雷達資料同化、模式嵌套、高分辨的區(qū)域數值模式、動力釋用技術等多方面系統(tǒng)地開展了降尺度方法研究。

    5.1.1 觀測降水降尺度

    基于新一代多普勒雷達探測技術和雨量計實測資料,采用卡爾曼濾波、最優(yōu)插值、變分、統(tǒng)計權重等多種校正方法,研發(fā)高分辨的定量降水估算校正技術[84-85]。萬玉發(fā)等提出了一種基于準同雨團樣本概念的雷達與雨量計實時同步結合的降雨估算新技術[86]。

    圖4 改進敏感區(qū)觀測后的低渦路徑預報試驗(標識1和2的點線分別為改進水汽場和風場的低渦預報路徑)

    5.1.2 預報降水降尺度

    基于中尺度數值模式通過資料同化、模式嵌套、動力釋用、空間統(tǒng)計等技術提高定量降水預報精度,提高降水分辨率,縮小了天氣模式與水文模型網格尺度的差異。例如,中國氣象局武漢暴雨研究所利用數值模式動力釋用降尺度技術將強降水構成分解為大尺度降水、不穩(wěn)定降水和地形影響降水,基于AREM模式利用高分辨率(2km×2km)的地形因子(海拔高度、坡度、坡向等)診斷數值模式預報降水,建立精簡實用的動力釋用降尺度模型,開展試驗流域降水預報降尺度試驗,獲取高時空分辨率(1h、2km×2km)的數值模式降水預報信息[87]。

    5.2 水文氣象耦合技術與實踐

    水文氣象耦合技術主要是利用氣象業(yè)務中雷達估算降雨及數值預報降雨技術來刻畫水文模型的輸入場,以改變目前水文模型中降雨量估算的落后方式,進一步提高降雨場估算和預報的時空精度,從而提高洪水模擬預報水平。

    5.2.1 雷達定量降水估算與水文模型的耦合

    雷達估算降水有時空分辨率高的優(yōu)點,可以比較客觀地反映降水量相對大小的分布趨勢。許繼軍等[88]、彭濤等[89]利用武漢、宜昌天氣雷達資料結合水文模型對三峽庫區(qū)清港河、香溪河、白蓮河等流域進行了水文預報試驗,結果表明:利用雷達估測降雨,與水文模型相匹配,較好地發(fā)揮了二者在表現降雨空間分布變化方面的優(yōu)勢,能夠有效地提高洪水預報精度。

    5.2.2 定量降水預報與水文模型的耦合

    目前,定量降水預報(QPF)應用于水文預報模型是被水文氣象學界普遍認同的發(fā)展方向之一。為此長江水利委員會、武漢大學、武漢暴雨所已開展過相關試驗:楊文發(fā)等[90]考慮預見期內的降水預報信息,采用降水預報與洪水預報耦合的方法,對1998年汛期發(fā)生在長江上游三峽區(qū)間的一次暴雨洪水過程進行了預報試驗;李超群、郭靖等[91-92]基于數值模式定量降水預報產品構建了清江隔河巖、漢江丹江口流域的水文氣象耦合模型,開展洪水預報的研究;崔春光、彭濤等[93-94]基于中尺度數值模式(AREM)將定量降水預報及其降水集合預報產品引入洪水預報領域,延長了洪水預見期,提高了洪水預報精度,拓展了洪水不確定性預報的新方法。

    5.3 水文氣象耦合的暴雨洪水預報系統(tǒng)

    洪水預報系統(tǒng)是在計算機上實現洪水預報聯機作業(yè)的運行系統(tǒng),它靠快速、準確地收集、存儲和處理水情、雨情,通過各種專業(yè)數學模型進行洪水預報,從而及時、準確地做出洪水流量過程的預報,提高了洪水預報的時效性和精確度。

    5.3.1 中國洪水預報系統(tǒng)

    在我國,洪水預報系統(tǒng)研制始于1985年,至20世紀90年代中期,我國從無到有基本上普及了洪水預報系統(tǒng)。當前,洪水預報系統(tǒng)研制進展較快,已開發(fā)了具有中國特色的專家交互式洪水預報系統(tǒng)。水利部水文局開發(fā)的“中國洪水預報系統(tǒng)”(CNFFS)和長江委水文局開發(fā)的“WIS水文預報平臺”(WISHFS)便是代表性的成果。中國洪水預報系統(tǒng)采用C/S結構,以全國統(tǒng)一的實時水情數據庫為依托,以地理信息系統(tǒng)為平臺,能方便地構建五類洪水預報方案,具有標準的、通用的預報模型方法庫,可任意選擇多模型、多方法制定預報方案,具有人工和自動優(yōu)選結合的模型率定功能,具有定時預報和人機交互預報功能,可干預任何信息源和預報過程,具有全面完善的系統(tǒng)管理功能等,已在25個流域機構和省(直轄市、自治區(qū))的水文部門推廣應用[95]。

    5.3.2 流域水文氣象實時預報系統(tǒng)

    近幾年,中國氣象局武漢暴雨研究所在GIS技術以及高分辨的DEM數據的支持下,基于定量降雨估算(QPE)、定量降雨預報(QPF)、水文預報模型和實時洪水預報技術,通過雷達測雨技術、中尺度暴雨預報模式技術獲取高時空分辨率的水文資料(降雨等)場輸入水文模型來進行實時水文氣象預報(見圖5),基于Asp.net技術自主設計研發(fā)了web形式的流域水文氣象實時預報系統(tǒng),并實現了與氣象信息綜合分析處理系統(tǒng)(MICAPS)有機融合,研發(fā)了基于MICAPS系統(tǒng)平臺的暴雨洪澇預報預警模塊。目前系統(tǒng)成功應用于三峽庫區(qū)、淮河王家壩、清江水布埡等12個流域[96]。

    6 討論與展望

    對近年來我國長江中游區(qū)域暴雨研究進展進行了回顧與總結,得出制約該區(qū)域暴雨研究深入與預報水平提高的因素主要體現在以下幾個方面。

    (1)已有的外場加密觀測試驗沒有得到能清晰反映暴雨β中尺度對流系統(tǒng)三維動力熱力結構、云微物理結構以及用于診斷分析其發(fā)生發(fā)展物理機制的加密觀測資料,更沒有針對暴雨數值預報模式物理過程方案研發(fā)和改進進行有針對性的觀測設計和野外觀測試驗。

    (2)長江中游處于復雜的中小尺度地形環(huán)境中,且引發(fā)暴雨的中尺度對流系統(tǒng)存在著多種組織類型。不同復雜地形環(huán)境下、不同組織類型的暴雨中尺度對流系統(tǒng)啟動、組織、發(fā)展、維持和消亡的主要物理機制是什么、之間存在著何種差異并不清楚。

    圖5 流域水文氣象實時預報結構及數據流程

    (3)基于高時空分辨率的衛(wèi)星、雷達、GPS水汽、風廓線儀、微波輻射計以及地面加密自動站等多源觀測,眾多資料亟待質量控制和最大效益的綜合應

    用,尤其是衛(wèi)星、雷達等高分辨率的遙感資料急需得到有效的融合和同化。區(qū)域暴雨數值預報模式中缺乏針對東亞天氣特點和青藏高原大地形影響下模式物理過程的描述,對暴雨預報的不確定性研究不夠充分。

    (4)水文預報模式和氣象預報模式耦合技術研發(fā)進展緩慢,時空分辨率的差異導致氣象模式與流域水文模型耦合效果不佳,定量降水估算產品在流域洪水預報中未能得到有效合理的應用。

    針對上述問題,未來我國長江中游區(qū)域暴雨研究可重點開展以下幾點工作。

    (1)開展青藏高原東側及其下游關鍵區(qū)大氣邊界層科學試驗、沿梅雨鋒切變線、低空急流帶、水汽輸送帶剖面觀測、深入內陸的臺風暴雨系統(tǒng)野外科學觀測,中尺度觀測資料融合與反演技術。

    (2)重點研究水汽輸送、對流擾動、大氣日變化、地形以及慣性重力波等因子對暴雨中尺度對流系統(tǒng)發(fā)生發(fā)展的影響機制,揭示長江中游暴雨天氣系統(tǒng)結構和形成機理。

    (3)開展多源探測資料的同化和云分析技術研究,發(fā)展高分辨率非靜力暴雨數值模式系統(tǒng)與適合我國區(qū)域特色的模式邊界層和云微物理過程方案。

    (4)突破以暴雨洪水為對象的水文氣象耦合關鍵技術,建立長江中游干、支流水文氣象預報系統(tǒng)。

    [1]陶詩言, 等. 中國之暴雨. 北京: 科學出版社, 1980.

    [2]丁一匯, 張建云, 等. 暴雨洪澇. 北京: 氣象出版社, 2009.

    [3]倪允琪, 劉黎平, 高梅, 等. 長江中下游梅雨鋒暴雨野外科學試驗.北京: 氣象出版社, 2004.

    [4]倪允琪, 周秀驥, 張人禾, 等. 我國南方暴雨的試驗與研究. 應用氣象學報, 2006, 17(6): 690-704.

    [5]崔春光, 倪允琪, 李紅莉, 等. 中國南方暴雨野外試驗(SCHeREX)中尺度氣象分析場的建立及其質量評估. 氣象學報, 2011, 69(1):26-40.

    [6]龍智勇, 石漢青, 黃思訓. 利用衛(wèi)星云圖反演云導風的新思路. 物理學報, 2011, 60(5): 834-839.

    [7]周振波, 閔錦忠, 彭霞云, 等. 單多普勒雷達風場反演的擴展方法(1): 方法與對比試驗. 高原氣象, 2006, 25(3): 516-524.

    [8]李華宏, 薛紀善, 王曼, 等. 多普勒雷達風廓線的反演及變分同化試驗. 應用氣象學報, 2007, 18(1): 50-57.

    [9]方德賢, 劉國慶, 董新寧, 等. 單多普勒雷達二維風場反演——Extended VPP方法. 氣象學報, 2007, 65(2): 231-240.

    [10]牟蓉, 劉黎平, 許小永, 等. 四維變分方法反演低層風場能力研究. 氣象, 2007, 33(1): 11-20.

    [11]劉黎平, 莫月琴, 沙雪松, 等. C波段雙多基地多普勒雷達資料處理和三維變分風場反演方法研究. 大氣科學, 2005, 29(6): 986-996.

    [12]劉黎平. 用雙多普勒雷達反演降水系統(tǒng)三維風場試驗研究. 應用氣象學報, 2003, 14(4): 502-504.

    [13]羅昌榮, 池艷珍, 周海光. 雙雷達反演臺風外圍強帶狀回波風場結構特征研究. 大氣科學, 2012, 36(2): 247-258.

    [14]彭亮, 陳洪濱, 李柏. 3mm多普勒云雷達測量反演云內空氣垂直速度的研究. 大氣科學, 2012, 36(1): 1-10.

    [15]劉貴華, 余興, 戴進. 不同積云微物理特征的衛(wèi)星反演分析. 南京氣象學院學報, 2009, 32(1): 17-23.

    [16] 徐小紅, 余興, 朱延年, 等. 一次強颮線云結構特征的衛(wèi)星反演分析. 高原氣象, 2012, 31(1): 258-268.

    [17]戴進, 余興, 劉貴華, 等. 一次暴雨過程中云微物理特征的衛(wèi)星反演分析. 氣象學報, 2010, 68(3): 387-397.

    [18]劉貴華, 余興, 師春香, 等. FY-3A/VIRR反演云微物理特征及與TERRA/MODIS反演結果的比較. 高原氣象, 2011, 30(2): 461-470.

    [19]陳英英. 利用FY-3A衛(wèi)星云圖對一次暴雨過程的特征云參數分析. 暴雨災害, 32(1): 24-31.

    [20]Lensky I M, Rosenfeld D. Clouds-aerosols-precipitation satellite analysis tool(CAPSAT). Atmos Chem Phys Discuss, 2008, 8: 4765-4890.

    [21]劉紅亞, 徐海明, 胡志晉, 等. 雷達反射率因子在中尺度云分辨模式初始化中的應用Ⅰ: 云微物理量和垂直速度的反演. 氣象學報, 2007, 65(6): 896-905.

    [22]劉黎平, 宗蓉, 齊彥斌, 等. 云雷達反演層狀云微物理參數及其與飛機觀測數據的對比. 中國工程科學, 2012, 14(9): 64-71.

    [23]王穎, 黃勇, 黃思源. 大氣溫濕廓線反演問題的研究. 國土資源遙感, 2008(1): 23-26.

    [24]潘寧, 郁凡. 用一維變分法校正衛(wèi)星資料反演地球大氣濕度場的試驗研究. 地球物理學報, 2009, 52(12): 2973-2986.

    [25]王靜, 趙朝方. 最小方差算法反演NOAA/AMSU海面大氣溫度垂直廓線. 海洋湖沼通報, 2008(2): 16-23.

    [26]曲思邈, 李國春. 利用MODIS紅外資料反演大氣溫濕度廓線的研究. 氣象與環(huán)境學報, 2012, 28(3): 2l-24.

    [27]王曦, 李萬彪. 應用FY-3A/MWHS資料反演太平洋海域晴空大氣濕度廓線. 熱帶氣象學報, 2013, 29(1): 47-54.

    [28]王祥, 趙冬至, 蘇岫, 等. 基于岸基實測數據的近紅外通道海洋大氣水汽反演. 紅外與毫米波學報, 2012, 31(6): 550-555.

    [29]龔建東. 同化技術: 數值天氣預報突破的關鍵——以歐洲中期天氣預報中心同化技術演進為例. 氣象科技進展, 2013, 3( 3): 6-13.

    [30]李紅莉, 崔春光, 王志斌. LAPS 的設計原理、模塊功能與產品應用. 暴雨災害, 2009, 28(1): 64-70.

    [31]Snook J S. Local domain forecasting support to the 1996 Atlanta Olympic Games. 12th Int Conf on Interactive and Processing Systems for Meteorology, Oceanography, and Hydrology. Atlanta: Amer Meteor Soc, 1996.

    [32]Birkenheuer D. Progress in applying GOES-derived data in local data assimilation. 10th Conf on Satellite Meteorology and Oceanography. Long Beach: Amer Meteor Soc, 2000.

    [33]李紅莉, 萬蓉, 謝有才. 利用LAPS 系統(tǒng)同化地基GPS 水汽資料的應用研究. 熱帶氣象學報, 2010, 26(6): 702-709.

    [34]李紅莉, 崔春光, 王志斌, 等. 中尺度分析系統(tǒng)LAPS 應用雷達資料的個例研究. 高原氣象, 2009, 28(6): 1443-1452.

    [35]Xie Y, Koch S, McGinley J, et al. A space-time multiscale analysis system: A sequential variational analysis approach. Mon Wea Rev, 2011, 139: 1224-1240.

    [36]Li W, Xie Y F, Deng S M, et al. Application of the multigrid method to the two-dimensional doppler radar radial velocity data assimilation. J Atmos Oceanic Technol, 2010, 27: 319-332.

    [37]Li W, Xie Y F, He Z J, et al. Application of the multigrid data assimilation scheme to the China seas’ temperature forecast. J Atmos Oceanic Technol, 2008, 25: 2106-2116.

    [38]王曉芳, 崔春光. 長江中下游地區(qū)梅雨期線狀中尺度對流系統(tǒng)分析Ⅰ: 組織類型特征. 氣象學報, 2012, 70(5): 909-923.

    [39]Zheng L L, Sun J H, Zhang X L, et al. Organizational modes of mesoscale convective systems over central East China. Weather and Forecasting, 2013, 28(5): 1081-1098.

    [40]Meng Z Y, Yan D C, Zhang Y J. General features of squall lines in East China. Monthly Weather Review, 2013, 141(5): 1629-1647.

    [41]王曉芳. 長江中下游地區(qū)梅雨期線狀中尺度對流系統(tǒng)分析Ⅱ:環(huán)境特征. 氣象學報, 2012, 70(5): 924-935.

    [42]鄭淋淋, 孫建華. 干、濕環(huán)境下中尺度對流系統(tǒng)發(fā)生的環(huán)流背景和地面特征分析. 大氣科學, 2013, 37(4): 891-904.

    [43]崔春光, 王曉芳, 付志康. 多源探測資料在一次非線狀MCS分析中的綜合應用. 氣象, 2013, 9(5): 556-566.

    [44]王曉芳, 胡伯威, 李燦. 湖北一次颮線過程的觀測分析及數值模擬. 高原氣象, 2010, 29(2): 471-485.

    [45]王曉芳, 胡伯威, 李紅莉, 等. 梅雨期一個伴有前導層狀降水對流線的結構特征. 高原氣象, 2011, 30(4): 1052-1066.

    [46]鄭永光, 陳炯, 朱佩君. 中國及周邊地區(qū)夏季中尺度對流系統(tǒng)分布及其日變化特征. 科學通報, 2008, 53(4): 471-481.

    [47]陳炯, 鄭永光, 張小玲, 等. 中國暖季短時強降水分布和日變化特征及其與中尺度對流系統(tǒng)日變化關系分析. 氣象學報, 2013, 71(3): 367-382.

    [48] 祁秀香, 鄭永光. 2007年夏季我國深對流活動時空分布特征. 應用氣象學報, 2009, 20(3): 286-294.

    [49]胡伯威. 梅雨鋒上MCS的發(fā)展、傳播以及與低層“濕度鋒”相關聯的CISK慣性重力波. 大氣科學, 2005, 2(6): 845-853.

    [50]王曉芳, 崔春光, 胡伯威. 與水平風切變強度不均勻相聯系的CISK慣性重力波. 應用氣象學報, 2007, 18(6): 760-768.

    [51]許小峰, 高麗, Jin F F. 梅雨濕度鋒上中尺度對流云團生成和發(fā)展的波捕獲機制. 自然科學進展,2009, 19(11): 1203-1208.

    [52]陳忠明. 非絕熱加熱、濕斜壓熱動力耦合強迫激發(fā)輻合增長的一種機制. 熱帶氣象學報, 2010, 26(6): 710-715.

    [53]Rotunno R, Klemp J B, Weisman M L. A theory for strong, longlived squall lines. J Atmos Sci, 1988, 45: 463-485.

    [54]Parker M D, Johnson R H. Simulated convective lines with leading precipitation. Part II: Evolution and Maintenance. J Atmos Sci, 2004, 61: 1656-1672.

    [55]Parker M D. Simulated convective lines with parallel stratiform precipitation. Part I: An archetype for convection in along-line shear. J Atmos Sci, 2007, 64: 267-288.

    [56]Bao X W, Tan Z M. A new mechanism of convective cell regeneration and development within a two-dimensional multicell storm. Acta Meteorological Sinica, 2009, 23(4): 414-427.

    [57]Smith R B. The influence of mountains on the atmosphere. Advances in Geophysics, 1979, 21: 87-230.

    [58]Smith R B. A linear upslope-time-delay model for orographic precipitation. J Hydro, 2003, 282: 2-9.

    [59]Smith R B, Barstad I. A linear theory of orographic precipitation. J Atmos Sci, 2004, 61 (12): 1377-1391.

    [60]趙玉春, 許小峰, 崔春光. 川西高原東坡地形對流暴雨的研究. 氣候與環(huán)境研究, 2012, 17(5): 607-616.

    [61]Hughes M M, HallA, Fovell R G. Blocking in areas of complex topography, and its in fl uence on rainfall distribution. J Atmos Sci, 2008,66(2): 508-518.

    [62]Zhao Y C. Numerical investigation of a localized extremely heavy rainfall event in complex topographic area during midsummer. Atmospheric Research, 2012, 113: 22-39.

    [63]趙玉春, 王葉紅, 崔春光. 大對流有效位能和條件不穩(wěn)定下地形對流降水的三維理想數值研究. 氣象學報, 2011, 69(5): 782-798.

    [64]趙玉春, 許小峰, 崔春光. 中尺度地形對梅雨鋒暴雨影響的個例研究. 高原氣象, 2012, 31(5): 1268-1282.

    [65]廖移山, 馮新, 石燕, 等. 2008年“7·22”襄樊特大暴雨的天氣學機理分析及地形的影響. 氣象學報, 2011, 69(6): 945-955.

    [66]宇如聰, 薛紀善, 徐幼平. AREMS中尺度暴雨數值預報模式系統(tǒng). 北京: 氣象出版社, 2004.

    [67]Mesinger F. A blocking technique for representation of mountains in atmospheric models. Rev Meteor Aeronaut, 1984, 44: 195-202.

    [68]Yu R C.Application of shape-preserving advection scheme to the moisture equation in an E-grid regional forecast model. Advances of Atmospheric Sciences, 1995, 12(1): 13-19.

    [69]宇如聰. E網格變量分布下差分格式的性質. 大氣科學, 1994, 18(2): 152-162.

    [70]王葉紅, 賴安偉, 趙玉春, 等. 降水資料同化在GRAPES-MESO模式中應用試驗研究. 大氣科學, 2013, 37(3): 645-667.

    [71]趙海貝, 王斌, 戴永久. 基于歷史樣本投影的四維變分陸面數據同化方法及其初步應用. 氣候與環(huán)境研究, 2009, 14(4): 383-389.

    [72]趙明, 王斌, 成巍, 等. 一種基于歷史樣本投影的三維變分同化方法. 氣候與環(huán)境研究, 2010, 15(2): 199-209.

    [73]蔣興文, 何光碧, 李躍清. AREM模式中引入MYJ邊界層方案介紹及其試驗. 高原山地氣象研究, 2010, 30(4): 1-7.

    [74]楊磊, 徐幼平, 程銳, 等. AREM模式中陸面過程的改進與試驗.暴雨災害, 2008, 27(1): 1-8.

    [75]李昀英, 葉成志, 鐘中. 陸面參數化方案對兩例不同類型暴雨可預報性的影響. 大氣科學, 2010, 34(2): 407-417.

    [76]王葉紅, 彭菊香, 公穎, 等. AREM-RUC 3h快速更新同化預報系統(tǒng)的建立與實時預報對比檢驗. 暴雨災害, 2011, 30(4): 296-304.

    [77]王葉紅, 趙玉春, 李紅莉, 等. AREM模式的熱啟動數值模擬——以2007年7月13日暴雨過程為例. 高原氣象, 2011, 30(6): 1488-1504.

    [78]王葉紅, 彭菊香, 趙玉春. 兩種初值方案對2008年中國汛期降水影響的數值試驗. 熱帶氣象學報, 2012, 28(6): 838-852.

    [79]李俊, 杜鈞, 王明歡, 等. 中尺度暴雨集合預報系統(tǒng)研發(fā)中的初值擾動試驗. 高原氣象, 2009, 28(6): 1365-1375.

    [80]李俊, 杜鈞, 王明歡, 等. AREM模式兩種初值擾動方案的集合降水預報試驗及檢驗. 熱帶氣象學報, 2010, 26(6): 217-228.

    [81]Li J, Du J, ZhangD L, et al. Ensemble-based analysis and sensitivity of mesoscale forecasts of a vortex over southwest China. Q J R, 2013, doi: 10.1002/qj.2200.

    [82]李俊, 王明歡, 公穎, 等. AREM短期集合預報系統(tǒng)及其降水預報檢驗. 暴雨災害, 2010, 29(1): 30-37.

    [83]陳金榮, 羅伯昆. 對水文氣象預報的認識和實踐. 人民長江, 1982(6): 52-59.

    [84]李建通, 楊維生. 提高最優(yōu)插值法測量區(qū)域降水量精度的探討.大氣科學, 2000, 24(2): 263-270.

    [85] 尹忠海, 張沛源.利用卡爾曼濾波校準方法估算區(qū)域降水量. 應用氣象學報, 2005, 16(2): 213-219.

    [86]萬玉發(fā), 吳翠紅, 金鴻祥. 基于準同雨團樣本概念雷達和雨量計的實時同步結合方法. 氣象學報, 2008, 66(1): 262-273.

    [87]殷志遠, 賴安偉, 公穎, 等. 氣象水文耦合中的降尺度方法研究進展. 暴雨災害, 2010, 29(1): 89-95.

    [88]許繼軍, 楊大文, 蔡治國. 分布式水文模型結合雷達測雨用于三峽區(qū)間的洪水預報. 長江科學院院報, 2007, 24(6): 42-48.

    [89]彭濤, 宋星原, 殷志遠, 等. 雷達定量估算降水在水文模式汛期洪水預報中的應用試驗. 氣象, 2010, 36(12): 50-55.

    [90]楊文發(fā), 等. 降水預報與洪水預報耦合應用初探. 水資源調查, 2003, 24(1): 38-40.

    [91]李超群, 郭生練, 張洪剛. 基于短期定量降水預報的隔河巖洪水預報研究. 水電能源科學, 2006, 24(4): 31-35.

    [92] 郭靖, 郭生練, 張俊, 等. 漢江流域未來降水徑流預測分析研究.水文, 2009(5): 18-22.

    [93]崔春光, 彭濤, 沈鐵元, 等. 定量降水預報與水文模型耦合的中小流域汛期洪水預報試驗. 氣象, 2010, 36(12): 56-61.

    [94]彭濤, 李俊, 殷志遠, 等. 基于集合降水預報產品的汛期洪水預報試驗. 暴雨災害, 2010, 29(3): 274-278.

    [95]劉志雨. 我國洪水預報技術研究進展與展望. 中國防汛抗旱, 2009(5): 13-16.

    [96] PengT, Wei C Z, Ye J T, et al. Basin hydrometeorology real-time forecast system. Information Computing and Applications, 2012, CCIS307: 665-671.

    Review of the Study on Heavy Rain in the Middle Reaches of the Yangtze River Since 2000

    Cui Chunguang, Lin Chunze, Wang Xiaofang, Li Jun, Peng Tao, Zhang Wengang
    (Hubei Key Laboratory for Heavy Rain Monitoring and Warning Research, Institute of Heavy Rain, China Meteorological Administration, Wuhan 430074)

    The middle reaches of the Yangtze River are frequently-occurring heavy rain areas. It is important to have the knowledge of heavy rain mesoscale and small-scale weather systems in these areas, their mechanism, and develop the heavy rain and f l ood forecasting technology for the f l ood control and hazard mitigation. From the heavy rain f i eld observation experiment, data inversion and fusion technology, observation analysis and theoretical research on heavy rain mesoscale convective systems, the mesoscale numerical model, and the hydrometeorological coupling rainstorm and f l ood forecasting system which are the most crucial aspects of observation data analysis, mechanism, numerical prediction, and hydrological application of the study on heavy rain in the middle reaches of the Yangtze River, the main research progress has been brief l y reviewed in China since 2000. The problems which obstruct the research depth and forecast quality were also discussed. Finally, the research development tendency was prospected.

    data fusion; mesoscale convective systems; hydrometeorological coupling

    10.3969/j.issn.2095-1973.2014.02.001

    2013年9月30日;

    2014年1月20日

    崔春光(1964—),Email: Cgcui@whihr.com.cn

    資助信息:公益性行業(yè)(氣象)科研專項(GYHY201106006、GYHY201206028),國家自然科學基金(40930951)

    猜你喜歡
    系統(tǒng)
    Smartflower POP 一體式光伏系統(tǒng)
    WJ-700無人機系統(tǒng)
    ZC系列無人機遙感系統(tǒng)
    北京測繪(2020年12期)2020-12-29 01:33:58
    基于PowerPC+FPGA顯示系統(tǒng)
    基于UG的發(fā)射箱自動化虛擬裝配系統(tǒng)開發(fā)
    半沸制皂系統(tǒng)(下)
    FAO系統(tǒng)特有功能分析及互聯互通探討
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    一德系統(tǒng) 德行天下
    PLC在多段調速系統(tǒng)中的應用
    18禁在线播放成人免费| 我的老师免费观看完整版| 白带黄色成豆腐渣| 国产白丝娇喘喷水9色精品| 九九在线视频观看精品| 日韩欧美 国产精品| 久久久亚洲精品成人影院| 亚洲激情五月婷婷啪啪| 国产午夜精品久久久久久一区二区三区| 亚洲精品日韩av片在线观看| 只有这里有精品99| 免费播放大片免费观看视频在线观看| 国产精品久久久久久精品电影| 免费看a级黄色片| 国产精品偷伦视频观看了| 黄片无遮挡物在线观看| 搡女人真爽免费视频火全软件| 国产精品人妻久久久影院| 国产男人的电影天堂91| 午夜激情久久久久久久| 国产一区亚洲一区在线观看| 九色成人免费人妻av| 久久影院123| 女人十人毛片免费观看3o分钟| 亚洲综合精品二区| 亚洲精品日本国产第一区| 国产中年淑女户外野战色| 91精品国产九色| 亚洲人成网站高清观看| 久久久久久久久久成人| 熟女电影av网| 97超视频在线观看视频| 69av精品久久久久久| 一级毛片我不卡| 国模一区二区三区四区视频| 99久国产av精品国产电影| 久久精品国产亚洲网站| 国产久久久一区二区三区| 内地一区二区视频在线| 一级爰片在线观看| 男女无遮挡免费网站观看| 亚洲精华国产精华液的使用体验| 国产欧美亚洲国产| 国产午夜福利久久久久久| 2021天堂中文幕一二区在线观| 亚洲欧美成人综合另类久久久| 大码成人一级视频| 在线观看国产h片| 少妇猛男粗大的猛烈进出视频 | 国产午夜精品一二区理论片| h日本视频在线播放| 汤姆久久久久久久影院中文字幕| 激情 狠狠 欧美| 国产精品99久久99久久久不卡 | 美女视频免费永久观看网站| 岛国毛片在线播放| 搡女人真爽免费视频火全软件| 午夜福利视频精品| 三级国产精品片| 老女人水多毛片| a级一级毛片免费在线观看| 99视频精品全部免费 在线| 欧美精品一区二区大全| av一本久久久久| 亚洲人与动物交配视频| av专区在线播放| 国产女主播在线喷水免费视频网站| 在线观看三级黄色| 欧美xxxx性猛交bbbb| 久久久精品免费免费高清| 91久久精品电影网| 色视频在线一区二区三区| 成人鲁丝片一二三区免费| 在线精品无人区一区二区三 | 一级毛片久久久久久久久女| 中文字幕制服av| 白带黄色成豆腐渣| 五月玫瑰六月丁香| 秋霞在线观看毛片| 久久鲁丝午夜福利片| 五月伊人婷婷丁香| 好男人视频免费观看在线| 日本黄色片子视频| 国产亚洲av片在线观看秒播厂| 午夜免费观看性视频| 女人久久www免费人成看片| 精品国产三级普通话版| 日本一二三区视频观看| 小蜜桃在线观看免费完整版高清| 国产一区二区三区av在线| 亚洲精品久久午夜乱码| 国产日韩欧美亚洲二区| 91久久精品电影网| 超碰av人人做人人爽久久| 一级av片app| 看免费成人av毛片| 99热网站在线观看| 久久久久网色| 80岁老熟妇乱子伦牲交| 五月伊人婷婷丁香| 最后的刺客免费高清国语| 亚洲av成人精品一区久久| 交换朋友夫妻互换小说| 国产黄片视频在线免费观看| 国产大屁股一区二区在线视频| 你懂的网址亚洲精品在线观看| 国产精品人妻久久久影院| 国产人妻一区二区三区在| 小蜜桃在线观看免费完整版高清| 看非洲黑人一级黄片| 亚洲婷婷狠狠爱综合网| 一级毛片aaaaaa免费看小| 91午夜精品亚洲一区二区三区| av在线天堂中文字幕| 国产精品伦人一区二区| 国产熟女欧美一区二区| 舔av片在线| 婷婷色av中文字幕| 精品国产一区二区三区久久久樱花 | 免费人成在线观看视频色| 神马国产精品三级电影在线观看| 91久久精品国产一区二区成人| 高清在线视频一区二区三区| 午夜福利网站1000一区二区三区| 亚洲av一区综合| 99热这里只有精品一区| 久久久久久九九精品二区国产| a级毛片免费高清观看在线播放| 天堂俺去俺来也www色官网| 欧美+日韩+精品| 欧美三级亚洲精品| 夫妻午夜视频| 99热6这里只有精品| 直男gayav资源| 欧美xxⅹ黑人| 视频中文字幕在线观看| 国产 一区 欧美 日韩| av网站免费在线观看视频| av在线天堂中文字幕| 国产91av在线免费观看| 永久免费av网站大全| 久久精品综合一区二区三区| 欧美成人午夜免费资源| 久久久久久九九精品二区国产| av在线天堂中文字幕| 亚洲人成网站在线播| 美女脱内裤让男人舔精品视频| 国产男人的电影天堂91| 国产精品一二三区在线看| 久久久久国产网址| 国产成人精品婷婷| 国产精品久久久久久久久免| 国产免费福利视频在线观看| 视频中文字幕在线观看| 日韩视频在线欧美| 成人黄色视频免费在线看| 大片免费播放器 马上看| 男人添女人高潮全过程视频| 国产一区二区在线观看日韩| 精品人妻熟女av久视频| 亚洲av二区三区四区| 一级片'在线观看视频| 成年免费大片在线观看| av.在线天堂| 久久久久久国产a免费观看| 99热这里只有是精品在线观看| 白带黄色成豆腐渣| 日韩强制内射视频| 男人狂女人下面高潮的视频| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 久久6这里有精品| 国产精品伦人一区二区| 亚洲无线观看免费| 久久热精品热| 各种免费的搞黄视频| 性插视频无遮挡在线免费观看| 男人狂女人下面高潮的视频| 日日摸夜夜添夜夜添av毛片| 寂寞人妻少妇视频99o| 各种免费的搞黄视频| 丰满少妇做爰视频| 高清毛片免费看| 亚洲国产日韩一区二区| 自拍偷自拍亚洲精品老妇| 国产真实伦视频高清在线观看| 一级二级三级毛片免费看| 免费在线观看成人毛片| 99久久精品一区二区三区| 在线看a的网站| 国产男人的电影天堂91| 免费黄网站久久成人精品| 国产精品一区二区性色av| 欧美性感艳星| 国产淫语在线视频| www.色视频.com| 秋霞在线观看毛片| 波野结衣二区三区在线| 在线观看美女被高潮喷水网站| 新久久久久国产一级毛片| 日韩av不卡免费在线播放| 欧美激情在线99| 一级毛片aaaaaa免费看小| eeuss影院久久| 国产淫片久久久久久久久| 特级一级黄色大片| 日本猛色少妇xxxxx猛交久久| 亚洲精品456在线播放app| 老司机影院毛片| 噜噜噜噜噜久久久久久91| 国产久久久一区二区三区| 男女边摸边吃奶| 午夜福利视频精品| 特级一级黄色大片| 丰满乱子伦码专区| 亚洲av成人精品一区久久| 少妇人妻一区二区三区视频| 人妻少妇偷人精品九色| 在线a可以看的网站| 欧美性猛交╳xxx乱大交人| 老女人水多毛片| 亚洲成人精品中文字幕电影| 日韩 亚洲 欧美在线| 网址你懂的国产日韩在线| 精品一区在线观看国产| 国产大屁股一区二区在线视频| 只有这里有精品99| 久久国内精品自在自线图片| 亚洲在久久综合| 午夜激情久久久久久久| 在线观看人妻少妇| 国产精品av视频在线免费观看| 一级黄片播放器| 人体艺术视频欧美日本| 日韩 亚洲 欧美在线| 成人欧美大片| 91午夜精品亚洲一区二区三区| 久久久精品欧美日韩精品| 五月玫瑰六月丁香| 欧美少妇被猛烈插入视频| 久久久久性生活片| 日韩一区二区视频免费看| 亚洲真实伦在线观看| 欧美精品国产亚洲| 好男人在线观看高清免费视频| 午夜老司机福利剧场| 久久久久久久国产电影| 国产欧美亚洲国产| 国产一区二区三区av在线| 免费观看性生交大片5| 亚洲国产精品国产精品| 国产永久视频网站| 一级二级三级毛片免费看| 51国产日韩欧美| 亚洲av中文字字幕乱码综合| 亚洲精品日韩在线中文字幕| 国产高清国产精品国产三级 | 男插女下体视频免费在线播放| 国产高清三级在线| 久久久久久九九精品二区国产| 亚洲最大成人av| 国产精品久久久久久精品古装| 九草在线视频观看| 国产伦精品一区二区三区视频9| 王馨瑶露胸无遮挡在线观看| 国产大屁股一区二区在线视频| 永久免费av网站大全| 我要看日韩黄色一级片| 亚洲国产成人一精品久久久| 黄色日韩在线| 国产精品三级大全| 国产日韩欧美在线精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜福利网站1000一区二区三区| 日本一二三区视频观看| 最近中文字幕2019免费版| 国产黄片视频在线免费观看| 伦精品一区二区三区| 久久久久久国产a免费观看| 老师上课跳d突然被开到最大视频| 91aial.com中文字幕在线观看| 国产精品一区www在线观看| 美女脱内裤让男人舔精品视频| 久久精品久久精品一区二区三区| 欧美激情国产日韩精品一区| 欧美一区二区亚洲| 国产探花在线观看一区二区| 中文精品一卡2卡3卡4更新| 别揉我奶头 嗯啊视频| 久久99蜜桃精品久久| 五月天丁香电影| 久久久久久久久久人人人人人人| 国产av国产精品国产| 人人妻人人看人人澡| 国产亚洲91精品色在线| av在线老鸭窝| 国产精品国产三级国产av玫瑰| 中文字幕免费在线视频6| 日产精品乱码卡一卡2卡三| 在线 av 中文字幕| 成年女人看的毛片在线观看| 能在线免费看毛片的网站| 国产欧美亚洲国产| 亚洲欧美日韩卡通动漫| 岛国毛片在线播放| 亚洲,一卡二卡三卡| 97超碰精品成人国产| 亚洲av免费在线观看| 丝瓜视频免费看黄片| 97人妻精品一区二区三区麻豆| 99热全是精品| freevideosex欧美| 韩国av在线不卡| 中文天堂在线官网| 精品一区二区三区视频在线| 国产成人精品久久久久久| 性色avwww在线观看| h日本视频在线播放| 69人妻影院| 尾随美女入室| 超碰av人人做人人爽久久| 国产成人午夜福利电影在线观看| 精品一区在线观看国产| 免费人成在线观看视频色| 免费不卡的大黄色大毛片视频在线观看| 免费黄网站久久成人精品| 街头女战士在线观看网站| 国产黄色视频一区二区在线观看| 中文天堂在线官网| 亚洲欧美一区二区三区国产| 亚洲精品视频女| 国产片特级美女逼逼视频| 国产有黄有色有爽视频| 一本久久精品| 你懂的网址亚洲精品在线观看| 18禁在线无遮挡免费观看视频| 国产毛片a区久久久久| 久久人人爽av亚洲精品天堂 | 自拍偷自拍亚洲精品老妇| 97精品久久久久久久久久精品| 狠狠精品人妻久久久久久综合| 深夜a级毛片| 国产亚洲一区二区精品| 日韩欧美精品v在线| 久久久久性生活片| 中文欧美无线码| 建设人人有责人人尽责人人享有的 | 99久久九九国产精品国产免费| 超碰av人人做人人爽久久| 熟女av电影| 国产精品爽爽va在线观看网站| 免费人成在线观看视频色| 日韩不卡一区二区三区视频在线| 国产精品久久久久久精品电影小说 | av福利片在线观看| 国产一区二区在线观看日韩| 18禁裸乳无遮挡免费网站照片| 成人美女网站在线观看视频| 22中文网久久字幕| 卡戴珊不雅视频在线播放| 欧美日韩亚洲高清精品| 九草在线视频观看| 少妇的逼好多水| 国产在线一区二区三区精| 亚洲精品自拍成人| 搡老乐熟女国产| 欧美成人午夜免费资源| 人妻夜夜爽99麻豆av| 午夜免费鲁丝| 国产高清三级在线| 婷婷色麻豆天堂久久| 久久久久九九精品影院| 另类亚洲欧美激情| 成人国产麻豆网| 在线免费十八禁| 大香蕉久久网| 夫妻午夜视频| 国产69精品久久久久777片| 亚洲成色77777| 激情五月婷婷亚洲| 欧美97在线视频| 久久人人爽人人爽人人片va| 欧美变态另类bdsm刘玥| 偷拍熟女少妇极品色| 人妻 亚洲 视频| av一本久久久久| kizo精华| 免费观看性生交大片5| 看免费成人av毛片| 97热精品久久久久久| 三级男女做爰猛烈吃奶摸视频| 嫩草影院入口| 国产有黄有色有爽视频| 中文乱码字字幕精品一区二区三区| 日本黄大片高清| 中文乱码字字幕精品一区二区三区| 亚洲精品aⅴ在线观看| 中文乱码字字幕精品一区二区三区| 亚洲av欧美aⅴ国产| 亚洲人与动物交配视频| 中文欧美无线码| 久热久热在线精品观看| 精品国产三级普通话版| 国产高清三级在线| 欧美国产精品一级二级三级 | 在线亚洲精品国产二区图片欧美 | 久热这里只有精品99| 亚洲精品国产成人久久av| 亚洲精品一二三| av在线天堂中文字幕| 在线精品无人区一区二区三 | av线在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 精品视频人人做人人爽| 舔av片在线| 搞女人的毛片| 一个人观看的视频www高清免费观看| 国产精品一区二区在线观看99| 性色avwww在线观看| 日本黄大片高清| 国产精品99久久99久久久不卡 | 乱码一卡2卡4卡精品| 一本色道久久久久久精品综合| 欧美97在线视频| 久久久久久久午夜电影| 伦精品一区二区三区| 欧美性感艳星| 80岁老熟妇乱子伦牲交| 亚洲精品乱久久久久久| 精品国产三级普通话版| 51国产日韩欧美| 国产黄片美女视频| 听说在线观看完整版免费高清| 日韩视频在线欧美| 亚洲av.av天堂| 欧美丝袜亚洲另类| 精品酒店卫生间| 欧美性猛交╳xxx乱大交人| 亚洲欧洲日产国产| 亚洲久久久久久中文字幕| 亚洲精品日韩在线中文字幕| 国产亚洲av嫩草精品影院| 白带黄色成豆腐渣| 麻豆成人av视频| 亚洲精品日韩av片在线观看| 一本久久精品| 欧美最新免费一区二区三区| 久久97久久精品| 91狼人影院| 丰满乱子伦码专区| 成人综合一区亚洲| 一级毛片 在线播放| 熟妇人妻不卡中文字幕| 国产真实伦视频高清在线观看| 狂野欧美激情性bbbbbb| 夜夜看夜夜爽夜夜摸| 性插视频无遮挡在线免费观看| 内射极品少妇av片p| a级毛片免费高清观看在线播放| 成人美女网站在线观看视频| 美女国产视频在线观看| 精品99又大又爽又粗少妇毛片| 欧美日韩在线观看h| 午夜激情福利司机影院| 交换朋友夫妻互换小说| 日韩亚洲欧美综合| 欧美三级亚洲精品| 99久国产av精品国产电影| 永久免费av网站大全| 久久久色成人| 午夜福利在线观看免费完整高清在| 亚洲精品国产色婷婷电影| 久久久精品94久久精品| tube8黄色片| 免费播放大片免费观看视频在线观看| a级毛色黄片| 成人毛片a级毛片在线播放| 久久久久久久精品精品| 18禁在线无遮挡免费观看视频| 少妇的逼水好多| 少妇的逼好多水| 天堂网av新在线| 国产真实伦视频高清在线观看| 久久久久精品性色| 一级毛片黄色毛片免费观看视频| 国产爽快片一区二区三区| 蜜臀久久99精品久久宅男| av女优亚洲男人天堂| 国产精品久久久久久久电影| 毛片一级片免费看久久久久| 十八禁网站网址无遮挡 | 欧美+日韩+精品| 狂野欧美激情性xxxx在线观看| a级一级毛片免费在线观看| 777米奇影视久久| 国产高清三级在线| 国产精品av视频在线免费观看| 亚洲人成网站在线播| 精华霜和精华液先用哪个| 亚洲成人一二三区av| 亚洲精品日韩av片在线观看| 又爽又黄无遮挡网站| 亚洲成人精品中文字幕电影| 热re99久久精品国产66热6| 国产精品成人在线| 91精品国产九色| 黄片wwwwww| 亚洲人成网站在线播| 男女下面进入的视频免费午夜| 97精品久久久久久久久久精品| 久久久久网色| 国产高潮美女av| 老司机影院成人| 精品熟女少妇av免费看| 好男人视频免费观看在线| 在线天堂最新版资源| 午夜老司机福利剧场| 成人亚洲精品一区在线观看 | 亚洲精品自拍成人| 黄色视频在线播放观看不卡| 国产永久视频网站| 欧美xxxx性猛交bbbb| freevideosex欧美| 三级国产精品片| av播播在线观看一区| 免费看日本二区| 美女主播在线视频| 欧美 日韩 精品 国产| 日韩伦理黄色片| 亚洲av在线观看美女高潮| 少妇人妻久久综合中文| 内地一区二区视频在线| 久久久成人免费电影| av在线老鸭窝| 午夜激情久久久久久久| 亚洲丝袜综合中文字幕| av一本久久久久| 色5月婷婷丁香| 日本av手机在线免费观看| 久久鲁丝午夜福利片| 国产成人a∨麻豆精品| av天堂中文字幕网| 日韩免费高清中文字幕av| 永久网站在线| 五月天丁香电影| 只有这里有精品99| 国产色爽女视频免费观看| 欧美日韩视频高清一区二区三区二| 亚洲第一区二区三区不卡| 国产69精品久久久久777片| 日韩三级伦理在线观看| 国产一区二区三区综合在线观看 | 国产在线男女| 国模一区二区三区四区视频| 亚洲精品久久久久久婷婷小说| 国产日韩欧美亚洲二区| 一级毛片久久久久久久久女| 三级国产精品欧美在线观看| 亚洲av不卡在线观看| 国产亚洲午夜精品一区二区久久 | 国产成人aa在线观看| 日本wwww免费看| 亚洲人成网站在线观看播放| 国产成人a∨麻豆精品| 国产黄色免费在线视频| 国产 一区 欧美 日韩| 午夜亚洲福利在线播放| 日本黄大片高清| 欧美成人精品欧美一级黄| 有码 亚洲区| av天堂中文字幕网| 色视频www国产| 一级毛片 在线播放| 精品少妇黑人巨大在线播放| 好男人在线观看高清免费视频| 久久精品久久久久久噜噜老黄| 小蜜桃在线观看免费完整版高清| 久久精品久久精品一区二区三区| h日本视频在线播放| 91在线精品国自产拍蜜月| 亚洲精品日本国产第一区| 亚洲av日韩在线播放| 黄色视频在线播放观看不卡| 欧美日韩综合久久久久久| 久久久久九九精品影院| 麻豆久久精品国产亚洲av| 亚洲欧美日韩另类电影网站 | 亚洲丝袜综合中文字幕| 国产精品爽爽va在线观看网站| 精品少妇黑人巨大在线播放| 性插视频无遮挡在线免费观看| 日本-黄色视频高清免费观看| 少妇 在线观看| 哪个播放器可以免费观看大片| 精品亚洲乱码少妇综合久久| 久久人人爽人人爽人人片va| 看黄色毛片网站| 黄色怎么调成土黄色| 国产老妇女一区| 国产 精品1| 久久国产乱子免费精品| 精品人妻一区二区三区麻豆| 一级片'在线观看视频| 91狼人影院| 欧美日韩国产mv在线观看视频 | 街头女战士在线观看网站| 亚洲精品乱码久久久v下载方式| 国产男人的电影天堂91| 尾随美女入室| 免费在线观看成人毛片| 亚洲精品成人久久久久久| 又爽又黄a免费视频| 一边亲一边摸免费视频| 特级一级黄色大片| 国模一区二区三区四区视频|