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

    雷達(dá)氣候研究進(jìn)展及其在城市區(qū)域強(qiáng)天氣臨近預(yù)報中的應(yīng)用

    2014-03-02 06:36:44陳明軒王迎春高峰肖現(xiàn)程叢蘭
    關(guān)鍵詞:風(fēng)暴對流氣候

    陳明軒王迎春高峰肖現(xiàn)程叢蘭

    (1 中國氣象局北京城市氣象研究所,北京 100089;2 北京市氣象局,北京 100089)

    雷達(dá)氣候研究進(jìn)展及其在城市區(qū)域強(qiáng)天氣臨近預(yù)報中的應(yīng)用

    陳明軒1王迎春2高峰1肖現(xiàn)1程叢蘭1

    (1 中國氣象局北京城市氣象研究所,北京 100089;2 北京市氣象局,北京 100089)

    首先介紹了近年來國際上在雷達(dá)氣候研究方面的一些進(jìn)展,特別聚焦于不同地區(qū)對流風(fēng)暴和對流性降水的日變化特征,以及京津冀地區(qū)雷達(dá)氣候研究的一些最新成果。主要揭示了不同地區(qū)對流的形成、加強(qiáng)和日循環(huán)傳播與太陽加熱、地形強(qiáng)迫、盛行風(fēng)及天氣尺度和中尺度系統(tǒng)活動等存在緊密聯(lián)系。其次,對雷達(dá)氣候?qū)W的另一特色研究——采用三維風(fēng)暴屬性追蹤統(tǒng)計的方法也進(jìn)行了介紹。文章也從城市化影響對流風(fēng)暴和對流性降水發(fā)展演變以及從雷達(dá)氣候?qū)W角度研究城市地區(qū)強(qiáng)天氣兩個方面,做了研究進(jìn)展的一些回顧。最后,對雷達(dá)氣候?qū)W在城市地區(qū)強(qiáng)天氣演變機(jī)理研究及精細(xì)化臨近預(yù)報預(yù)警方面的應(yīng)用進(jìn)行了討論。

    雷達(dá)氣候,對流,降水,城市化,臨近預(yù)報

    1 引言

    天氣雷達(dá)可以獲得高時空分辨率的探測資料,是目前對流風(fēng)暴三維特征和強(qiáng)降水最為快速有效的氣象探測工具之一[1]。如果能夠?qū)δ骋坏貐^(qū)長時間的天氣雷達(dá)資料(目前主要是單部雷達(dá)或者雷達(dá)組網(wǎng)觀測得到的反射率因子資料)進(jìn)行統(tǒng)計分析研究,就能夠了解對流風(fēng)暴和對流性降水的精細(xì)時空分布和演變特征,也包括強(qiáng)天氣的出現(xiàn)頻率和風(fēng)暴屬性等,這就是所謂的雷達(dá)氣候研究。自從20世紀(jì)70年代開始,隨著計算機(jī)和數(shù)字化天氣雷達(dá)技術(shù)的發(fā)展,使得利用計算機(jī)處理海量雷達(dá)探測資料成為可能時,就已經(jīng)有學(xué)者開始了雷達(dá)氣候研究。Kuo等[2]利用一個X波段試驗雷達(dá)的4年夏季觀測資料,開展了美國Black Hills地區(qū)的雷達(dá)氣候研究,首次揭示了該地區(qū)對流降水回波的時空分布特征及其在不同盛行風(fēng)形勢下時空分布頻率的差異,可以說是為雷達(dá)氣候研究開創(chuàng)了先河。近幾十年來,隨著天氣雷達(dá)在世界各地

    的廣泛業(yè)務(wù)應(yīng)用和雷達(dá)資料的不斷積累,雷達(dá)氣候研究已經(jīng)成為對流天氣氣候特征研究的主要方向之一,使用的雷達(dá)資料也從單部雷達(dá)的回波觀測轉(zhuǎn)換為使用三維雷達(dá)組網(wǎng)拼圖資料。特別是在美國,自20世紀(jì)80年代后期WSR-88D業(yè)務(wù)雷達(dá)在全美的布網(wǎng)應(yīng)用開始,隨著雷達(dá)資料的積累,已經(jīng)開展了大量的雷達(dá)氣候研究方面的工作。在國際上,也有學(xué)者利用長時間的雷達(dá)估測降水資料進(jìn)行對流性降水的氣候統(tǒng)計研究,這也可以歸屬于雷達(dá)氣候研究范疇。通過雷達(dá)氣候研究,不但能夠為理解強(qiáng)對流天氣的發(fā)展演變機(jī)制提供有效幫助,為對流風(fēng)暴和對流性強(qiáng)降水的臨近預(yù)報和預(yù)警提供背景資料和科學(xué)依據(jù),也能夠為評估和發(fā)展對流尺度數(shù)值預(yù)報、定量降水預(yù)報以及天氣和氣候系統(tǒng)模式提供一個評判基準(zhǔn)。

    隨著全球城市化進(jìn)程的加劇和大型城市的不斷擴(kuò)展,城市和城市群地區(qū)強(qiáng)對流天氣的發(fā)展演變越來越受到人們的關(guān)注。強(qiáng)對流天氣經(jīng)常會給城市地區(qū)帶來重大的災(zāi)害損失,包括人員傷亡、交通阻塞、城市建筑受損等。而城市下墊面的改變和城市熱島效應(yīng)[3-5]等也會對強(qiáng)對流天氣的發(fā)生和發(fā)展產(chǎn)生明顯影響。如果將雷達(dá)氣候與城市氣象特別是城市地區(qū)對流天氣研究相結(jié)合,將有助于理解城市地區(qū)強(qiáng)對流天氣的發(fā)展演變特征以及城市環(huán)境對強(qiáng)對流天氣的影響機(jī)理,對提升城市地區(qū)強(qiáng)天氣的發(fā)展演變和臨近預(yù)警起到重要推動作用。

    本文主要介紹近年來國內(nèi)外特別是國際上關(guān)于雷達(dá)氣候研究方面的一些進(jìn)展,并從城市化影響對流風(fēng)暴和強(qiáng)降水發(fā)展演變以及從雷達(dá)氣候角度研究城市地區(qū)強(qiáng)天氣兩個方面,做了一些研究進(jìn)展的簡介。最后,對雷達(dá)氣候研究及其在城市地區(qū)強(qiáng)天氣臨近預(yù)報方面的應(yīng)用前景進(jìn)行了討論。

    2 雷達(dá)氣候研究進(jìn)展

    2.1 對流日變化氣候特征及盛行風(fēng)影響

    雷達(dá)探測資料具有高時空分辨率的特性,對揭示對流天氣(包括對流風(fēng)暴和對流性降水)時空分布的短周期變化特別是日變化特征具有獨特優(yōu)勢。而對流的日變化特征對強(qiáng)天氣的臨近預(yù)報預(yù)警以及對流尺度數(shù)值天氣預(yù)報物理過程的改進(jìn)等均具有重要意義。因此,對流天氣日變化特征研究成為雷達(dá)氣候研究的主要內(nèi)容之一。世界各地的學(xué)者對此開展了大量的研究,揭示了不同地區(qū)對流風(fēng)暴和對流性降水的日變化特征,以及不同天氣系統(tǒng)、氣候特征和地形強(qiáng)迫對對流日變化的影響機(jī)理。在美國,Carbone等[6-7]和Ahijevych等[8-9]通過對美國大陸多年的WSR-88D雷達(dá)資料的統(tǒng)計研究,分析了美國大陸對流性降水回波出現(xiàn)頻率的時空分布特征,揭示了降水回波存在明顯的日循環(huán)和半日循環(huán)的周期性傳播氣候特征,以及由于落基山等山區(qū)的動力抬升和局地?zé)崃?qiáng)迫作用,導(dǎo)致降水系統(tǒng)存在明顯的山區(qū)新生并向下游(向東)不斷傳播的特征(如圖1和圖2所示)。Carbone等[6]的研究同時也指出,對流的日循環(huán)特征主要是由于太陽加熱作用的日變化而形成的,并且具有不可預(yù)報性,因為數(shù)值預(yù)報包括資料同化系統(tǒng)主要是基于準(zhǔn)平衡氣流和天氣尺度強(qiáng)迫的動力學(xué)基礎(chǔ)而設(shè)計的,其對由于天氣尺度而形成的具有較長周期性變化和生命史特征的對流系統(tǒng)具有較高的可預(yù)報性,而盛夏期間主要以熱力對流系統(tǒng)為主,所以數(shù)值預(yù)報的準(zhǔn)確性要差很多。Hsu等[10]利用覆蓋美國大陸的8年暖季雷達(dá)估測降水資料,通過波譜分解方法,得出了美國大陸對流性降水存在半天、1天、4天等多個周期性變化特征,其中日循環(huán)特征最為明顯,而且小于5天的高頻振蕩表明中尺度對流系統(tǒng)主要以向東傳播為主,與Carbone等[6]得

    出的結(jié)論相似。Lang等[11]和Rowe等[12]利用北美季風(fēng)試驗加強(qiáng)觀測期間的雷達(dá)觀測資料,對2004年7—8月發(fā)生在美國西南部及墨西哥灣北部的強(qiáng)對流個例進(jìn)行了統(tǒng)計分析,發(fā)現(xiàn)地形強(qiáng)迫對對流風(fēng)暴的每日新生和傳播起到關(guān)鍵作用,由于局地加熱對流一般在山區(qū)生成,然后在向下游傳播的過程中得到明顯加強(qiáng)。在歐洲,Weckwerth等[13]利用多年的雷達(dá)觀測資料,對德國西南部到法國北部地區(qū)的對流風(fēng)暴新生和加強(qiáng)的氣候特征進(jìn)行了研究,揭示了這一地區(qū)的對流風(fēng)暴存在顯著的日循環(huán)和頻率峰值特征,每日加熱和山區(qū)地形強(qiáng)迫對風(fēng)暴的新生和加強(qiáng)起到關(guān)鍵作用,對流首先在山區(qū)大量形成,大約2h后到達(dá)平原地區(qū)。Pavai等[14]利用中國研制的WSR-98D雷達(dá)的多年探測資料,研究了羅馬尼亞Bobohalma地區(qū)對流出現(xiàn)頻率的變化特征,發(fā)現(xiàn)這一地區(qū)夏季對流的日變化與地形高度存在緊密關(guān)系,對流頻繁地在較高地形的地區(qū)新生,然后向平原地區(qū)傳播。

    圖1 美國大陸6—8月雷達(dá)回波出現(xiàn)頻率[7](a)2100UTC;(b)0200UTC;(c)0700UTC;(d)1300UTC

    在中國,近年來主要基于地面觀測和衛(wèi)星資料,開展了不同地區(qū)的對流氣候特征研究,特別是對流性降水的日變化特征研究,取得了很多重要的成果[15-24],這里不做詳細(xì)介紹。然而,基于高時空分辨率的雷達(dá)探測資料的對流氣候特征研究工作則相對較少。Lin等[25]利用多年的雷達(dá)探測資料以及閃電、探空和地面資料等,研究了臺灣地區(qū)在弱天氣尺度強(qiáng)迫形勢下午后對流風(fēng)暴的氣候特征,發(fā)現(xiàn)臺灣地區(qū)在弱天氣尺度形勢下的對流風(fēng)暴出現(xiàn)頻率存在顯著的午后最大峰值特征,對流頻繁出現(xiàn)的地區(qū)與臺灣的地形分布一致,且在弱天氣尺度形勢下,對流不存在明顯的傳播特征。Yeung等[26]對36個影響香港地區(qū)的熱帶氣旋個例的雷達(dá)回波和降水特征以及地形作用進(jìn)行了統(tǒng)計分析研究,為熱帶氣旋及其降水臨近預(yù)報提供了背景資料。

    隨著中國大陸新一代天氣雷達(dá)布網(wǎng)的完成和雷達(dá)資料的多年積累,在我國不同地區(qū)開展精細(xì)的雷達(dá)氣候研究的資料條件已經(jīng)逐漸成熟。特別是在京津冀地區(qū),為做好北京奧運會等大型活動的氣象保障服務(wù),已經(jīng)在北京及其周邊地區(qū)建立了日趨完善的新一代天氣雷達(dá)網(wǎng),并且具備了較為可靠的實時雷達(dá)資料質(zhì)量控制體系,也積累了多年的雷達(dá)組網(wǎng)資料[27-28]。這些條件為這一地區(qū)的雷達(dá)氣候研究奠定了良好的基礎(chǔ)。Chen等[29]首次利用2008—2011年暖季5—8月份京津冀地區(qū)6部新一代天氣雷達(dá)逐6min的反射率因子三維拼圖資料,研究了這一地區(qū)對流風(fēng)暴出現(xiàn)頻率的暖季平均日變化特征,以及暖季主要月份(6、7、8月)對流風(fēng)暴日變化的差異。暖季平均日變化特征研究揭示了京津冀西北部山區(qū)是暖季對流活動的主要源地,由于太陽加熱,西北部山區(qū)在午后風(fēng)暴出現(xiàn)頻率顯著增加,且大多數(shù)的對流風(fēng)暴從西北部山區(qū)向山腳甚至東南平原地區(qū)傳播和發(fā)展。到下午接近傍晚時段,在京津冀山腳附近,特別是在北京及其周邊地區(qū),對流風(fēng)暴出現(xiàn)頻率顯著增大,風(fēng)暴的增強(qiáng)和新生與地形對低層偏南暖濕氣流的強(qiáng)迫上升存在密切聯(lián)系(圖3)。在夜間到凌晨時段,對流風(fēng)暴主要在東南部平原地區(qū)發(fā)展,風(fēng)暴頻率隨著時間的推移逐漸減小,風(fēng)暴頻率最低時段發(fā)生在早晨。對比圖2和圖3可以看出,京津冀地區(qū)的對流風(fēng)暴日循環(huán)和傳播特征與Carbone等[7]對美國大陸的對流風(fēng)暴日循環(huán)和傳播的研究結(jié)果存在一定的相似性。上述研究結(jié)論也與Zheng等[16]以及He等[22]利用衛(wèi)星資料開展的這一地區(qū)對流日變化特征的研究結(jié)論類似。從Chen等[29]的研究結(jié)果來看,區(qū)域平均的風(fēng)暴頻率顯示,京津冀地區(qū)的對流風(fēng)暴存在顯著的雙峰值特征:主要峰值是由于地形強(qiáng)迫導(dǎo)致的午后到傍晚時段山腳附近的對流頻發(fā)和增強(qiáng),次要峰值是由于平原地區(qū)有利的觸發(fā)機(jī)制導(dǎo)致的夜對流頻發(fā)(圖4)。He等[22]利用數(shù)值模式資料分析后發(fā)現(xiàn),華北平原地區(qū)有利的夜對流觸發(fā)機(jī)制與山區(qū)—平原的次級環(huán)流上升支以及夜間東南平原地區(qū)的低空西南急流存在密切聯(lián)系。在Chen等[29]的研究中,通過對比暖季各月的對流風(fēng)暴出現(xiàn)頻率特征發(fā)現(xiàn),西北部山區(qū)午后熱

    對流出現(xiàn)頻率最高在6月,最低在8月,而在平原地區(qū)午后對流出現(xiàn)頻率最高在7月。地形強(qiáng)迫效應(yīng)對對流風(fēng)暴的增強(qiáng)和新生的影響,表現(xiàn)最為明顯是在7月,其次是6月。對于平原地區(qū)夜對流而言,出現(xiàn)頻率最高在7月和8月,而6月較低。暖季各月的風(fēng)暴出現(xiàn)頻率的平均日變化趨勢和雙峰值特征與整個暖季的平均趨勢相似,表明6—8月的對流風(fēng)暴控制了整個暖季這一地區(qū)對流風(fēng)暴的日變化特征。平均而言,7月這一地區(qū)的暖季對流風(fēng)暴(包括長生命史的中尺度對流系統(tǒng)和短時生命史的對流單體)的出現(xiàn)頻率最高(圖4)。對流風(fēng)暴在6月和7月存在明顯的自西北向東南的傳播特征,且6月傳播速度最快。

    圖2 美國大陸6—8月雷達(dá)回波出現(xiàn)頻率Hovm?ller圖[7](沿緯度平均的時間—經(jīng)度圖)(a)對應(yīng)的平均地形廓線;(b) 1997—2000年資料統(tǒng)計結(jié)果;(c)1996—2007年資料統(tǒng)計結(jié)果

    上述國內(nèi)外關(guān)于對流日變化氣候特征的研究是目前雷達(dá)氣候研究的主流,其成果主要揭示了不同地區(qū)對流的形成、加強(qiáng)和日循環(huán)傳播與太陽加熱、地形強(qiáng)迫以及天氣尺度和中尺度系統(tǒng)的活動存在緊密聯(lián)系。一般來說,山區(qū)作為主要的熱源,是夏季午后對流的主要源地[9,11-13,25,29-30]。在一些地區(qū),對流具有明顯的向下游移動和發(fā)展傳播的特征,而且在對流下山的過程中,由于地形強(qiáng)迫,使得低層吹來的暖濕空氣產(chǎn)生強(qiáng)烈抬升,從而導(dǎo)致對流(包括對流風(fēng)暴和對流性降水)在山腳附近明顯加強(qiáng),也頻繁出現(xiàn)對流風(fēng)暴的新生[8-13,29];而在另一些地區(qū),特別是以山區(qū)為主的地區(qū)(如我國的臺灣地區(qū)),如果沒有強(qiáng)的天氣尺度系統(tǒng)的影響,對流的傳播特征也不是非常明顯[25,30]。

    早在20世紀(jì)70年代,Kuo等[2]開展的雷達(dá)氣候研究就揭示了暖季對流的發(fā)展演變和時空分布特征與盛行風(fēng)存在密切聯(lián)系。他們發(fā)現(xiàn),對于美國Black Hills地區(qū)來說,因為不同方向的盛行風(fēng)對應(yīng)于不同的天氣系統(tǒng),對流回波在西南、西北向的盛行風(fēng)形勢下出現(xiàn)的頻率最高。近年來,一系列的雷達(dá)氣候研究進(jìn)一步表明,對流風(fēng)暴和對流性降水每日新生、演變和傳播與盛行風(fēng)存在密切關(guān)系[11,26,30-31]。Lang等[11]和Rowe等[12]的研究發(fā)現(xiàn),由于盛行風(fēng)與環(huán)境垂直風(fēng)切變、大氣波動等存在密切聯(lián)系,再加上地形的影響,因此對對流的演變傳播和時空分布特征產(chǎn)生重要影響。Yeung等[26]和Saxen等[30]的雷達(dá)氣候研究也發(fā)現(xiàn),不同的盛行風(fēng)形勢下,對流的出現(xiàn)頻率差異明顯。對流發(fā)展演變與盛行風(fēng)的密切關(guān)系,也被數(shù)值模擬研究所證實[32]。最近,利用京津冀地區(qū)6部雷達(dá)4年暖季的反射率因子資料,計算了500和925hPa

    不用盛行風(fēng)向下的對流風(fēng)暴的出現(xiàn)頻率分布,可以看出其存在顯著的差異:在500hPa,西西南(WSW)和西西北(WNW)風(fēng)向形勢下風(fēng)暴出現(xiàn)頻率顯著偏高;在925hPa,南西南(SSW)和南東南(SSE)風(fēng)向形勢下,風(fēng)暴出現(xiàn)頻率顯著偏高(圖6)。初步分析也表明,在500hPa WSW和WNW風(fēng)向下,對應(yīng)的是偏西氣流,主要與影響京津冀地區(qū)暖季對流風(fēng)暴的三類天氣系統(tǒng)(冷渦、西來槽、西太平洋副熱帶高壓)存在直接聯(lián)系。在925hPa SSW和SSE風(fēng)向下,對應(yīng)的是京津冀南部平原地區(qū)的低層偏南風(fēng),有利于低層偏南暖濕氣流的向北推進(jìn),并在地形強(qiáng)迫抬升作用下,對對流的發(fā)生和加強(qiáng)起到關(guān)鍵作用。

    圖3 京津冀地區(qū)暖季40dBz以上雷達(dá)回波出現(xiàn)頻率Hovm?ller圖[29](a)緯向平均結(jié)果;(b)徑向平均結(jié)果;(c)西南—東北向平均結(jié)果(每幅圖下方的藍(lán)綠色區(qū)域為對應(yīng)的地形平均廓線)

    圖4 暖季平均的逐小時40dBz以上雷達(dá)回波出現(xiàn)的百分率[29](a)京津冀地區(qū)的平均結(jié)果;(b)北京地區(qū)的平均結(jié)果(圖b中的粗黑線表示在每個小時時段過去24h雨量記錄達(dá)到50mm以上的北京地區(qū)自動站站數(shù))

    如雷達(dá)氣候研究結(jié)果所言,不同高度的盛行風(fēng)代表了不同的天氣系統(tǒng)和形勢。一般而言,500~700hPa高度的盛行風(fēng)代表了影響對流移動和傳播的天氣系統(tǒng)或引導(dǎo)氣流,而低層(如925hPa高度)的盛行風(fēng)則代表了低空急流的方向和強(qiáng)度,這些均與對流的每日發(fā)展演變和傳播存在密切關(guān)系。進(jìn)一步地,這種高低空風(fēng)形勢相配合所形成的垂直風(fēng)切變(包括風(fēng)速切變和風(fēng)向切變)以及地形對低空急流帶來的暖濕空氣的強(qiáng)迫抬升作用,對對流日變化也具有顯著的正面影響。這從Parker等[33]針對美國大陸對流日變化以及Lang等[11]和Rowe等[12]針對北美季風(fēng)試驗區(qū)域(墨西哥灣北部地區(qū))對流日變化的雷達(dá)氣候研究中也均得到了很好的驗證。

    另外,基于雷達(dá)觀測資料,也可以開展對流風(fēng)暴和對流性降水的月季、季節(jié)、年季變化等的氣候特征研究。Hocker等[34]利用10年的雷達(dá)反射率因子資料,研究了發(fā)生在美國Oklahoma地區(qū)的颮線系統(tǒng)的出現(xiàn)頻率年季變化、月季變化、日變化以及時空分布和移動路徑等的氣候特征。Murray等[35]利用10年的雷達(dá)資料,研究了美國東北部地區(qū)對流風(fēng)暴的時空分布特征,揭示了這一地區(qū)的對流風(fēng)暴時空分布存在明顯的年季差異,以及與天氣尺度系統(tǒng)的年季活動規(guī)律密切相關(guān)。Carbone等[36]也嘗試?yán)?年暖季的雷達(dá)資料研究揭示美國大陸對流的年季變化,并對北美季風(fēng)環(huán)流和ENSO對對流發(fā)展演變的影響進(jìn)行了初步分析。但是,類似的雷達(dá)氣候研究均需要利用較長時間積累的雷達(dá)資料,這樣從氣候?qū)W意義上來講才具有一定的可靠性,這方面的工作目前開展得還不是很多。

    圖5 京津冀區(qū)域平均的暖季6—8月各月逐小時40dBz以上雷達(dá)回波出現(xiàn)百分率(a),以及出現(xiàn)頻率的月—小時分布圖(b )[29]

    2.2 風(fēng)暴屬性的識別追蹤和統(tǒng)計

    圖6 40dBz以上雷達(dá)反射率因子指示的京津冀區(qū)域平均的逐小時風(fēng)暴出現(xiàn)頻率的風(fēng)向—小時分布圖(a)500hPa;(b)925hPa

    對于雷達(dá)氣候研究,也有一些學(xué)者采用風(fēng)暴屬性追蹤統(tǒng)計的方法進(jìn)行。這種方法首先是基于雷達(dá)三維反射率因子回波,利用風(fēng)暴識別和追蹤分析算法對每個風(fēng)暴單體的屬性進(jìn)行識別、分析和追蹤,然后再對得到的所有風(fēng)暴的有效屬性進(jìn)行統(tǒng)計分析,得出某一地區(qū)對流風(fēng)暴屬性及發(fā)展演變和傳播的氣候特征。目前,比較有代表性的對流風(fēng)暴識別和追蹤分析算法包括風(fēng)暴三維屬性追蹤算法TITAN[37]以及風(fēng)暴質(zhì)心追蹤算法SCIT[38],采用風(fēng)暴屬性追蹤統(tǒng)計的雷達(dá)氣候研

    究也基本上基于這兩種算法。Henry[39]利用TITAN,通過對美國丹佛附近的雷達(dá)觀測資料進(jìn)行風(fēng)暴追蹤和統(tǒng)計研究,揭示了83%的單單體風(fēng)暴的生命史不會超過30min,而超過88%的多單體風(fēng)暴和超級單體風(fēng)暴的生命史超過30min,為丹佛地區(qū)對流風(fēng)暴尺寸和強(qiáng)度的臨近預(yù)報提供了改進(jìn)的科學(xué)基礎(chǔ)。Saxen等[30]利用TITAN算法對美國一個多山地區(qū)的對流風(fēng)暴屬性的追蹤統(tǒng)計也表明,風(fēng)暴的生命史、平均面積、移動傳播速度、強(qiáng)度以及回波頂高等,均表現(xiàn)出顯著的差異和明顯的峰值特征。在歐洲和南美洲,也有學(xué)者嘗試?yán)肨ITAN進(jìn)行風(fēng)暴屬性的識別和追蹤統(tǒng)計,并得出了一些有意義的、定量的風(fēng)暴氣候?qū)W屬性特征,為對流的診斷分析和臨近預(yù)報提供了科學(xué)背景支持[40-41]。Mohee等[42]利用SCIT算法,對美國北Dakota州的對流風(fēng)暴屬性的氣候特征進(jìn)行了研究,發(fā)現(xiàn)風(fēng)暴出現(xiàn)頻率在傍晚到凌晨存在明顯的峰值,風(fēng)暴的平均生命史為24min左右,平均移速為16.4m/s左右,風(fēng)暴向偏東北方向的移動趨勢最為明顯,在暖季每月的變化差異也很明顯,并開展了一系列的對流風(fēng)暴氣候特征屬性的分析。

    在國內(nèi),韓雷等[43]首次利用2003—2007年5—8月天津塘沽的新一代多普勒天氣雷達(dá)三維反射率因子數(shù)據(jù),對天津及鄰近地區(qū)的對流風(fēng)暴的氣候分布特征進(jìn)行了統(tǒng)計分析。通過對風(fēng)暴的識別和追蹤結(jié)果的統(tǒng)計分析,給出了京津地區(qū)對流風(fēng)暴的面積、體積、頂高、最大反射率因子、生命史和運動規(guī)律等定量的強(qiáng)對流風(fēng)暴氣候統(tǒng)計特征。從韓雷等[43]給出的對流風(fēng)暴面積和體積的地理分布可以看出(圖7),北京西北部和西南部地區(qū)風(fēng)暴的面積較小,而東部則明顯較大。根據(jù)Chen等[29]對該地區(qū)對流日變化的雷達(dá)氣候研究,并依據(jù)天氣學(xué)分析,可以推測出西北部山區(qū)的對流風(fēng)暴可能更多為熱力對流,空間與時間尺度都較小,而在東部平原地區(qū)出現(xiàn)具有較好組織性和長生命史特征的對流風(fēng)暴系統(tǒng)的頻率明顯較高。目前,作者正在利用京津冀地區(qū)6部雷達(dá)5年(2008—2012年)暖季5—9月的反射率因子三維拼圖資料和TITAN算法,對這一地區(qū)的對流風(fēng)暴三維屬性進(jìn)行細(xì)致的氣候特征統(tǒng)計研究,相關(guān)的工作這里不再介紹。

    對基于高時空分辨率雷達(dá)三維資料追蹤識別得到的對流風(fēng)暴的三維屬性進(jìn)行定量統(tǒng)計分析,是雷達(dá)氣候研究的一大特色,目前來說,也是其他資料和方法所不能比擬的。與使用常規(guī)觀測資料和靜止衛(wèi)星紅外亮溫的統(tǒng)計結(jié)果相比,使用雷達(dá)資料進(jìn)行這樣的氣候統(tǒng)計不僅可以給出對流風(fēng)暴的三維空間信息的統(tǒng)計結(jié)果(如體積和頂高等),還可以給出風(fēng)暴的生命史、移速、移向等定量氣候分布特征,從而可進(jìn)一步充實對某一地區(qū)對流天氣氣候特征的認(rèn)識。

    2.3 基于雷達(dá)氣候研究的個例分析

    圖7 基于TITAN追蹤識別結(jié)果統(tǒng)計得到的天津及其鄰近地區(qū)對流風(fēng)暴平均面積(a)和體積(b)的地理分布[43]

    這里,基于京津冀地區(qū)雷達(dá)氣候初步研究結(jié)果[29],對三個影響北京地區(qū)的典型對流風(fēng)暴個例進(jìn)行對比分析,分別是2012年“7·21”特大暴雨以及2013年“6·24”和“6·28”強(qiáng)對流天氣個例。圖8是40dBz以上雷達(dá)組合反射率因子回波的時間序列圖。由圖可以看出,不同天氣形勢下的對流風(fēng)暴,均存在明顯的由山區(qū)和山腳向下游平原地區(qū)傳播發(fā)展的特征。而且,從這三個個例來看,在接近傍晚時段,在京津冀山腳附近特別是在北京及其周邊地區(qū),對流風(fēng)暴的面積明顯增大,組織性得到顯著加強(qiáng)。從北京南郊觀象臺的探空來看(圖略),對于這三個個例,低層925hPa環(huán)境風(fēng)均為明顯的偏南風(fēng),這表明對流風(fēng)暴的

    增強(qiáng)和新生與地形對低層偏南暖濕氣流的強(qiáng)迫上升存在密切聯(lián)系。上述個例分析結(jié)果與Chen等[29]的雷達(dá)氣候?qū)W研究結(jié)果一致。對比這三個個例來看,“7·21”特大暴雨個例的雷達(dá)回波傳播速度明顯偏慢,從山腳傳播到渤海灣附近大約耗時9h,而“6·24”和“6·28”強(qiáng)對流個例的雷達(dá)回波傳播速度較快,從山腳傳播到渤海灣僅耗時4h左右。Chen等[29]的雷達(dá)氣候?qū)W研究指出,京津冀地區(qū)的對流風(fēng)暴在6和7月存在明顯的自西北向東南的傳播特征,且6月傳播速度最快,也符合這三個個例的實際觀測結(jié)果。從500hPa探空風(fēng)來看,“7·21”和“6·28”個例500hPa環(huán)境風(fēng)為WSW,風(fēng)暴向下游的傳播方向更偏東,而“6·24”個例的500hPa環(huán)境風(fēng)則為WNW,風(fēng)暴的傳播方向更偏向東南。這也與我們最新得出的雷達(dá)氣候研究結(jié)果一致。

    精細(xì)的雷達(dá)氣候研究結(jié)果是強(qiáng)對流天氣精細(xì)化預(yù)報所必須的科學(xué)背景,可以用來指導(dǎo)預(yù)報人員對實際的強(qiáng)天氣個例進(jìn)行預(yù)報,也可以將雷達(dá)氣候研究結(jié)果集成進(jìn)自動化的強(qiáng)天氣客觀預(yù)報系統(tǒng)中,用來發(fā)展和改進(jìn)強(qiáng)天氣的自動化客觀預(yù)報和預(yù)警[25,28,30]。

    圖8 京津冀地區(qū)雷達(dá)反射率因子≥40dBz的對流風(fēng)暴回波時間序列(a)2012年7月21日;(b)2013年6月24日;(c)2013年6月28日(圖中粗藍(lán)線為200m地形等高線;圖中不同顏色代表不同時間(北京時)的40dBz以上回波覆蓋范圍)

    3 雷達(dá)氣候在城市強(qiáng)天氣臨近預(yù)報中的應(yīng)用

    3.1 城市對局地降水及強(qiáng)天氣的影響

    近年來,國際上開展了大量的城市化和城市區(qū)域特征對局地降水及強(qiáng)天氣影響的研究。對于局地降水,大部分研究都表明,城市會導(dǎo)致局地降水量及強(qiáng)降水頻次的增加,尤其在城市的下風(fēng)方,這種效應(yīng)更加明顯[44-50]。但是,部分學(xué)者認(rèn)為,城市效應(yīng)會導(dǎo)致局地降水有不同程度的減少[51-53]??傊?,目前城市對局地降水的可能影響問題,國際上存在兩種相反的觀點:多數(shù)研究者認(rèn)為,城市的動力、熱力作用使城區(qū)和城市下游地區(qū)降水增加;有一些研究者認(rèn)為,城市大氣污染物的微物理過程會使城市下游地區(qū)的降水減少。我國就城市對局地降水的影響也開展了一些研究,得到的結(jié)論也不盡相同[54-57]。關(guān)于城市對局地強(qiáng)對流天氣的影響,國內(nèi)外研究表明,由于城市規(guī)模的不斷擴(kuò)大,導(dǎo)致城市下墊面變化并引起城市熱島等城市效應(yīng)的不斷加強(qiáng),非常有利于近地面輻合的增強(qiáng)和局地?zé)崃Νh(huán)流的產(chǎn)生,導(dǎo)致暖濕氣流上升強(qiáng)度加大,從而引起城市地區(qū)強(qiáng)對流天氣的新生和強(qiáng)烈發(fā)展[58-68]??偟膩碚f,城市化和城市特征對局地降水和強(qiáng)對流天氣的影響存在較強(qiáng)個性,不同地區(qū)的研究結(jié)論不盡相同。在我國城市化加劇的今天,城市對局地降水和強(qiáng)對流天氣的影響究竟如何,降水及對流的新生、增強(qiáng)或減弱與熱島強(qiáng)度、城市邊界層特征等存在怎樣的一種關(guān)系,尚需進(jìn)一步研究。

    3.2 雷達(dá)氣候與城市強(qiáng)天氣臨近預(yù)報的應(yīng)用結(jié)合

    從上文的介紹來看,一方面,基于高時空分辨率雷達(dá)觀測資料的雷達(dá)氣候研究,能夠揭示對流風(fēng)暴和對流性降水的精細(xì)時空演變特征;另一方面,雖然總體來說城市區(qū)域的尺度相對較小,但城市對對流天氣的影響卻非常顯著。因此,如果利用雷達(dá)氣候研究來揭示城市地區(qū)對流的時空變化特征,就具有一定的優(yōu)越性。但是,就目前來說,國內(nèi)外在這方面開展的工作還不是很多。Bentley等[69]利用10年的雷達(dá)反射率因子拼圖資料,分析了美國Atlanta地區(qū)的對流氣候變化特征,發(fā)現(xiàn)在Atlanta的城市地區(qū)以及城市擴(kuò)展地區(qū),與強(qiáng)的城市熱島效應(yīng)相對應(yīng),對流風(fēng)暴的出現(xiàn)頻率明顯偏高,表明城市對局地強(qiáng)對流的發(fā)生和增強(qiáng)具有明顯正面效應(yīng)。另外,城市下墊面的改變,也能夠引起對流風(fēng)暴的動力和熱動力特征的明顯變化。最近,Niyogi等[70]開展了美國印第安納州最大城市和首府Indianapolis地區(qū)對流風(fēng)暴的雷達(dá)氣候研究,并對城市地區(qū)和周邊非城市地區(qū)的對流特征進(jìn)行了比較,揭示了大城市對區(qū)域?qū)α鞯挠绊?。他們發(fā)現(xiàn),在城市地區(qū)超過60%的對流風(fēng)暴會改變結(jié)構(gòu),而在非城市地區(qū),

    僅為25%。白天的城市環(huán)境對對流風(fēng)暴的影響明顯高于夜晚。雷達(dá)氣候分析指出這一地區(qū)的對流風(fēng)暴在接近城市地區(qū)上風(fēng)方時易出現(xiàn)分裂而在下風(fēng)方時又極易出現(xiàn)合并,使得在城市地區(qū)上風(fēng)方主要以中等尺寸的對流風(fēng)暴為主,而在下風(fēng)方則是小尺寸和大尺寸風(fēng)暴出現(xiàn)的概率較高。同時,Niyogi等[70]開展的敏感性數(shù)值模擬試驗也表明,如果在模式中去除Indianapolis城市地區(qū)邊界層等的影響時,導(dǎo)致模式中的區(qū)域輻合和對流、地面能量平衡、邊界層結(jié)構(gòu)等發(fā)生明顯變化,從而導(dǎo)致實際發(fā)生的對流風(fēng)暴不能夠被很好模擬出來。除了利用雷達(dá)資料外,Ntelekos等[71]也利用閃電資料和地面觀測資料,分析了美國馬里蘭州Baltimore附近的城市群地區(qū)暖季對流風(fēng)暴的時空分布特征,包括季節(jié)變化、日變化和空間頻率等,發(fā)現(xiàn)風(fēng)暴頻率在城市地區(qū)與非城市地區(qū)存在顯著差異,并指出這與城市熱島等城市化效應(yīng)存在密切聯(lián)系。Hand等[31]也利用地面觀測和衛(wèi)星資料,分析了美國Oklahoma城附近的暖季降水氣候特征,發(fā)現(xiàn)在弱天氣尺度強(qiáng)迫形勢下,城市效應(yīng)對降水影響最為明顯,城市不同地區(qū)的濕度分布非常不均勻,城市下風(fēng)方地區(qū)濕度最大,降水也經(jīng)常出現(xiàn)異常。在上述Chen等[29]針對京津冀地區(qū)的雷達(dá)氣候研究中,也透視出北京城市效應(yīng)對對流風(fēng)暴的發(fā)生發(fā)展起到了關(guān)鍵作用。從北京城區(qū)附近風(fēng)暴概率變化的時段來看,20:00—21:00之間觀測的城區(qū)內(nèi)對流風(fēng)暴出現(xiàn)頻率仍然較低,而從21:00開始至00:00,城區(qū)內(nèi)觀測的孤立對流風(fēng)暴單體出現(xiàn)的頻率突然增大(圖9)。上述分析表明此段時間內(nèi)北京城市效應(yīng)對對流風(fēng)暴的加強(qiáng)或者新生的作用最為明顯,這可能與城市熱島及城市陸面不均勻性等存在密切關(guān)系,但還需要進(jìn)一步去研究。上述為數(shù)不多的雷達(dá)氣候?qū)W與城市對流特征相結(jié)合的研究表明,城市局地環(huán)流特征(如城市熱島環(huán)流)對局地強(qiáng)對流的發(fā)生和增強(qiáng)具有明顯的正面效應(yīng)。另外,城市特殊的下墊面和邊界層特征,也對對流風(fēng)暴的結(jié)構(gòu)和熱動力特征產(chǎn)生明顯影響。

    圖9 北京地區(qū)雷達(dá)反射率因子≥40dBz的對流風(fēng)暴回波出現(xiàn)頻率(a)20:00—21:00;(b)21:00—22:00;(c)22:00—23:00;(d)23:00—00:00(北京時;圖中粗黑線為200m地形等高線;修改自文獻(xiàn)[29]的圖2)

    可以看出,如果將雷達(dá)氣候研究與城市地區(qū)對流天氣研究相結(jié)合,不但可以掌握城市地區(qū)和城市化對對流風(fēng)暴和對流性降水的影響,也可以改進(jìn)對城市地區(qū)對流天氣的臨近預(yù)報和預(yù)警。在國際上,雷達(dá)氣候研究結(jié)果作為一個區(qū)域強(qiáng)對流天氣發(fā)生發(fā)展和演變傳播的氣候背景資料,已經(jīng)被應(yīng)用于局地對流天氣的自動臨近預(yù)報試驗中[30,72]。在北京自動臨近預(yù)報系統(tǒng)(BJANC)的參數(shù)調(diào)整優(yōu)化和系統(tǒng)評估中[28],我們也已經(jīng)初步應(yīng)用了一些雷達(dá)氣候研究結(jié)果[43,73],但如何將上述最新的京津冀雷達(dá)氣候研究結(jié)果[29]應(yīng)用于對流天氣的自動臨近預(yù)報中,許多方面還有待于進(jìn)一步深入研究。更為重要的是,如何將雷達(dá)氣候研究專門應(yīng)用于城市地區(qū)對流天氣精細(xì)的臨近預(yù)報和預(yù)警中,目前還未發(fā)現(xiàn)國際上有任何的報道或出版的論文?;谝酝难芯拷?jīng)驗,可以認(rèn)為,對于一個特定的城市或城市群地區(qū),發(fā)展一套集雷達(dá)氣候研究結(jié)果、城市熱島效應(yīng)、城市下墊面影響效應(yīng)、城市冠層和城市水文特征等在內(nèi)的城市地區(qū)對流天氣概念模型,同時在模型中也必須考慮一些精細(xì)的對流熱動力特征[74-80],并將概念模型集成進(jìn)對流天氣臨近預(yù)報的專家系統(tǒng)中(可利用類似于文獻(xiàn)[28]和[72]的模糊邏輯集成方法),再結(jié)合城市地區(qū)精細(xì)的探測資料分析和地理信息系統(tǒng)(GIS),將能夠有效提升對城市地區(qū)對流天氣和對流性降水的監(jiān)測分析和臨近預(yù)報、預(yù)警能力。另外,利用臨近預(yù)報和數(shù)值預(yù)報的融合技術(shù)(Blending),也是提升對流尺度定量降水預(yù)報特別是城市地區(qū)暴洪短時臨近預(yù)報能力的重要手段[81-82]。

    4 結(jié)束語

    通過對國內(nèi)外雷達(dá)氣候研究成果的描述可以看出,雷達(dá)氣候研究對于理解一個區(qū)域內(nèi)對流風(fēng)暴和對流性降水的精細(xì)時空分布和演變特征非常關(guān)鍵,也非常有助于作為中小尺度天氣現(xiàn)象的對流天氣的臨近預(yù)報和預(yù)警能力以及定量降水預(yù)報能力的提升。從本文的介紹可以看出,國際上特別是在美國,已經(jīng)針對不同區(qū)域的對流風(fēng)暴和對流性降水,開展了大量的雷達(dá)氣候研究。目前,美國國家海洋和大氣管理局(NOAA)強(qiáng)風(fēng)暴預(yù)測中心已經(jīng)有專門的網(wǎng)頁提供美國本土幾乎每部雷達(dá)站覆蓋范圍內(nèi)的強(qiáng)天氣氣候特征圖(包括短時大風(fēng)、冰雹、龍卷的出現(xiàn)報告統(tǒng)計圖,強(qiáng)天氣年際變化、月季變化、日變化圖),以及相關(guān)的分析表格,供預(yù)報人員在進(jìn)行強(qiáng)天氣臨近預(yù)報和預(yù)警業(yè)務(wù)時參考(網(wǎng)址:www. spc.noaa.gov/climo/online/rda/)。美國佛羅里達(dá)州立大學(xué)也制作了一套佛羅里達(dá)州和南佐治亞州的雷達(dá)氣候?qū)W數(shù)據(jù)集,為佛羅里達(dá)州環(huán)保部門的地面降水和地下水資源模擬以及NOAA的水文氣象分析和預(yù)報等提供科學(xué)參考(網(wǎng)址:fuelberg.met.fsu.edu/ research/fdep/gquina/climo.html)。本文也介紹了雷達(dá)氣候研究結(jié)果可在強(qiáng)天氣自動臨近預(yù)報專家系統(tǒng)中得到應(yīng)用。另外,雷達(dá)氣候研究結(jié)果可用于雷達(dá)雜波和雷達(dá)觀測偏差分析,從而有助于改進(jìn)雷達(dá)的觀測質(zhì)量,進(jìn)而也能夠改進(jìn)雷達(dá)的定量降水估測效果[83-84]。但是,我國在雷達(dá)氣候研究方面開展的系統(tǒng)性工作還不是很多。目前,我國新一代天氣雷達(dá)已經(jīng)基本完成布網(wǎng),部分地區(qū)也已經(jīng)積累了多年的雷達(dá)觀測資料,開展雷達(dá)氣候研究的資料條件業(yè)已成熟。作者基于上海及其周邊地區(qū)6部新一代天氣雷達(dá)在2010年5月1日—9月14日上海世博會期間的逐6min反射率因子觀測資料,統(tǒng)計得到了上海及其周邊地區(qū)2010年暖季的對流風(fēng)暴日循環(huán)特征(圖10)。從這個初步的雷達(dá)氣候結(jié)果可以看出,至少在2010年暖季,上海及其周邊的對流風(fēng)暴存在明顯的日循環(huán)和傳播特征,其中影響上海地區(qū)的對流峰值主要出現(xiàn)在10:00—12:00之間,但是不同時段影響的區(qū)域完全不同。

    另一方面,近年來,隨著國家經(jīng)濟(jì)發(fā)展和戰(zhàn)略調(diào)整,我國城市化進(jìn)程的速度顯著加快,不但大城市發(fā)展規(guī)模逐漸擴(kuò)大,而且也已經(jīng)涌現(xiàn)出了“長三角”、“珠三角”、“環(huán)渤?!钡却笾谐鞘性萍某鞘腥旱貐^(qū)。國內(nèi)外研究也說明,城市環(huán)境對對流天氣的影響非常明顯。而且更為重要的是,城市地區(qū)對對流天氣災(zāi)害的承載能力更加脆弱,需要更為準(zhǔn)確、及時和精細(xì)化的對流天氣臨近預(yù)報和預(yù)警服務(wù)。因此,如前文所言,這就促使我們應(yīng)該將雷達(dá)氣候研究與城市地區(qū)對流天氣研究有效結(jié)合,從而提升城市地區(qū)對流天氣的精細(xì)臨近預(yù)報和預(yù)警能力。在國際上,也已經(jīng)開始嘗試將雷達(dá)氣候研究與地理信息系統(tǒng)(GIS)相結(jié)合,對對流風(fēng)暴和對流性降水的氣候特征進(jìn)行更為精細(xì)的研究[34,85],這也為城市(群)地區(qū)精細(xì)尺度的對流天氣研究與雷達(dá)氣候研究及城市環(huán)境特征相結(jié)合提供了一種新思路。

    當(dāng)然,雷達(dá)氣候研究也具有一定不確定性。首先,天氣雷達(dá)的探測資料具有一定的局限性,比如雷達(dá)地物雜波(包括城市區(qū)域的高建筑物雜波)、亮帶回波、波束阻擋、不同雷達(dá)之間的標(biāo)定差異等,均影響到天氣雷達(dá)對真實對流回波的探測,進(jìn)而影響到雷達(dá)氣候研究的效果[86]。其次,在雷達(dá)站稀疏的地區(qū),一些區(qū)域成為雷達(dá)探測的“空缺”區(qū),而這些區(qū)域?qū)⒉豢赡塬@得有效的雷達(dá)氣候研究結(jié)果[86]。另外,正如

    Chen等[29]所強(qiáng)調(diào)的一樣,利用幾年的雷達(dá)資料進(jìn)行暖季對流天氣氣候特征的研究,就需要確保影響這一地區(qū)暖季的大氣環(huán)流不出現(xiàn)異常,否則,雷達(dá)氣候研究結(jié)果可能不具有代表性。在如何將雷達(dá)氣候研究有效應(yīng)用于城市地區(qū)對流天氣的臨近預(yù)報和預(yù)警的研究和業(yè)務(wù)中,也面臨許多的課題。

    圖10 上海及其周邊地區(qū)雷達(dá)反射率因子≥35dBz的對流風(fēng)暴回波出現(xiàn)概率(%)(圖中粗黑線表示上海地區(qū))

    [1]俞小鼎, 姚秀萍, 熊廷南, 等. 多普勒天氣雷達(dá)原理與業(yè)務(wù)應(yīng)用.北京: 氣象出版社, 2006.

    [2]Kuo J T, Orville H D. A radar climatology of summertime convective clouds in the Black Hills. J Appl Meteor, 1973, 12(2): 359-368.

    [3]Rao P K. Remote sensing of urban heat islands from an environmental satellite. Bull Amer Meteor Soc, 1972, 53(8):

    647-648.

    [4]Carlson T N, Augustin J A, Boland F E. Potential application of satellite temperatures measurements in the analysis of land use over urban areas. Bull Amer Meteor Soc, 1977, 58(12): 1301-1303.

    [5] MatsonM, McClain E P, McGinnis D F, et al. Satellite detection of urban heat islands. Mon Wea Rev, 1978, 106(12): 1725-1734.

    [6]Carbone R E, Tuttle J D, Ahijevych D A, et al. Inferences of predictability associated with warm season precipitation episodes. J Atmos Sci, 2002, 59(13): 2033-2056.

    [7]Carbone R E, Tuttle J D. Rainfall occurrence in the U.S. warm season: The diurnal cycle. J Climate, 2008, 21(16), 4132-4146.

    [8]Ahijevych D A, Carbone R E, Davis C A. Regional-scale aspects of the diurnal precipitation cycle. 31st Inter Conf on Radar Meteor, Amer Meteor Soc, Seattle, WA, 2003.

    [9]Ahijevych D A, Davis C A, Carbone R E, et al. Initiation of precipitation episodes relative to elevated terrain. J Atmos Sci, 2004, 61(22): 2763-2769.

    [10]Hsu H M, Moncrie ff M W, Tung W W, et al. Multiscale temporal variability of warm-season precipitation over North America: Statistical analysis of radar measurements. J Atmos Sci, 2006, 63(9): 2355-2368.

    [11]Lang T J, Ahijevych D A, Nesbitt S W, et al. Radar-observed characteristics of precipitating systems during NAME 2004. J Climate, 2007, 20(9): 1713-1733.

    [12]Rowe A K, Rutledge S A, Lang T J, et al. Elevation-dependent trends in precipitation observed during NAME. Mon Wea Rev, 2008, 136(12): 4962-4979.

    [13]Weckwerth T M, Wilson J W, Hagen M, et al. Radar climatology of the COPS region. Q J R Meteor Soc, 2011, 137: 31-41.

    [14]Pavai C, Vamos C. Statistical characteristics of re fl ectivity fi eld in areas covered by the Bobohalma WSR-98D radar. Romanian J Meteor, 2006, 8(1-2): 1-20.

    [15]Wang C C, Chen G T J, Carbone R E. A climatology of warmseason cloud patterns over East Asia based on GMS infrared brightness temperature observations. Mon Wea Rev, 2004, 132(7): 1606-1629.

    [16]Zheng Y, Chen J, Chen M, et al. Statistic characteristics and weather significance of infrared TBB during May-August in Beijing and its vicinity. Chin Sci Bull, 2007, 52(24): 3428-3435.

    [17]Zhang H, Zhai P. Temporal and spatial characteristics of extreme hourly precipitation over Eastern China in the warm season. Adv Atmos Sci, 2011, 28(5): 1177-1183.

    [18]Zhai P, Zhang X, Wan H, et al. Trends in total precipitation and frequency of daily precipitation extremes over China. J Climate, 2005, 18(7): 1096-1108.

    [19]Yu R, Zhou T, Xiong A, et al. Diurnal variation of summer precipitation over contiguous China. Geophys Res Lett, 2007, 34: L01704, doi: 10.1029/2006GL028129.

    [20]Yu R, Xu Y, Zhou T, et al. Relation between rainfall duration and diurnal variation in the warm season precipitation over central eastern China. Geophys Res Lett, 2007, 34: L13703, doi: 10.1029/2007GL030315.

    [21]Yin S, Chen D, Xie Y. Diurnal variations of precipitation during the warm season over China. Int J Climatol, 2009, 29: 1154-1170.

    [22]He H, Zhang F. Diurnal variations of warm-season precipitation over Northern China. Mon Wea Rev, 2010, 138(4): 1017-1025.

    [23]Chen G, Sha W, Iwasaki T. Diurnal variation of precipitation over southeastern China: Spatial distribution and its seasonality. J Geophys Res, 2009, 114: D13103, doi: 10.1029/2008JD011103.

    [24]Chen G, Sha W, Iwasaki T, et al. Diurnal variation of rainfall in the Yangtze River Valley during the spring-summer transition from TRMM measurements. J Geophys Res, 2012, 117: D06106, doi:10.1029/2011JD017056.

    [25]Lin P F, Chang P L, Jou B J D, et al. Warm season afternoon thunderstorm characteristics under weak synoptic-scale forcing over Taiwan Island. Wea Forecasting, 2011, 26(1): 44-60.

    [26]Yeung L H, Chan S, Cheng P. Radar climatology of Hong Kong and its application to landfalling tropical cyclone rainfall estimation. 35th Inter Conf on Radar Meteor, Amer Meteor Soc, Pittsburgh, Pa, 2011.

    [27]陳明軒, 高峰. 利用一種自動識別算法移除天氣雷達(dá)反射率因子中的亮帶. 應(yīng)用氣象學(xué)報, 2006, 17(2): 207-214.

    [28]陳明軒, 高峰, 孔榮, 等. 自動臨近預(yù)報系統(tǒng)及其在北京奧運期間的應(yīng)用. 應(yīng)用氣象學(xué)報, 2010, 21 (4): 395-404.

    [29]Chen M, Wang Y, Gao F, et al. Diurnal variations in convective storm activity over contiguous North China during the warmseason based on radar mosaic climatology. J Geophys Res, 2012, 117: D20115, doi:10.1029/2012JD018158.

    [30]Saxen T R, Mueller C K, Warner T T, et al. The operational mesogamma-scale analysis and forecast system of the U.S. Army Test and Evaluation Command. Part IV: The White Sands Missile Range Auto-Nowcast System. J Appl Meteor Climatol, 2008, 47(4): 1123-1139.

    [31]Hand L M, Shepherd J M. An investigation of warm-season spatial rainfall variability in Oklahoma City: Possible linkages to urbanization and prevailing wind. J Appl Meteor Climatol, 2009, 48(2): 251-269.

    [32]Tucker D F, Crook N A. Flow over heated terrain. Part II: Generation of convective precipitation. Mon Wea Rev, 2005, 133(9): 2565-2582.

    [33]Parker M D, Ahijevych D A. Convective episodes in the East-Central United States. Mon Wea Rev, 2007, 135(11): 3707-3727.

    [34]Hocker J E, Basara J B. A 10-year spatial climatology of squall line storms across Oklahoma. Int J Climatol, 2008, 28: 765-775.

    [35]Murray J C, Colle B A. The spatial and temporal variability of convective storms over the Northeast United States during the Warm Season. Mon Wea Rev, 2011, 139(3): 992-1012.

    [36]Carbone R E, Tuttle J D, Ahijevych D A. Inter-annual and semidiurnal variations in summertime precipitation. 31st Inter Conf on Radar Meteor, Amer Meteor Soc, Boston, 2003.

    [37]Dixon M, Wiener G. TITAN: Thunderstorm identification, tracking, analysis, and nowcasting—A radar-based methodology. J Atmos Oceanic Technol, 1993, 10 (6): 785-797.

    [38]Johnson J T, MacKeen P L, Witt A, et al. The storm cell identification and tracking algorithm: An enhanced WSR-88D algorithm. Wea Forecasting, 1998, 13(2): 263-276.

    [39]Henry S G. Analysis of thunderstorm lifetime as a function of size and intensity. 26th Conf on Radar Meteor, Amer Meteor Soc, Norman, OK, 1993.

    [40]Goudenhoofdt E, Delobbe L. Statistical characteristics of convective storms in Belgium derived from volumetric weather radar observations. J Appl Meteor Climatol, 2013, 52(4): 918-934.

    [41]Gomes A M, Held G. Characterization of storm properties during the TroCCiBras Experiment. 32nd Conf on Radar Meteor, Amer Meteor Soc, Albuquerque, 2005.

    [42]Mohee F M, Miller C. Climatology of thunderstorms for North Dakota, 2002-06. J Appl Meteor Climatol, 2010, 49(9): 1881-1890.

    [43]韓雷, 俞小鼎, 鄭永光, 等. 京津及鄰近地區(qū)暖季強(qiáng)對流風(fēng)暴的氣候分布特征. 科學(xué)通報, 2009, 54(11): 1585-1590.

    [44]Takahashi H. Secular variation in the occurrence property of summertime daily rainfall amount in and around the Tokyo Metropolitan area (in Japanese with an English abstract). Tenki, 2003, 50: 31-41.

    [45]Shepherd J M. A review of current investigations of urban-induced rainfall and recommendations for the future. Earth Interactions, 2005, 9: 1-27.

    [46]Diem J E, Brown D P. Anthropogenic impacts on summer precipitation in central Arizona. U S A Prof Geogr, 2003, 55 (3): 343-355.

    [47]Matthew S, Raman S, Suresh R, et al. Urban e ff ects of Chennai on sea breeze induced convection and precipitation. J Earth Sys Sci, 2008, 117(6): 897-909.

    [48]Burian S J, Shepherd J M. E ff ects of urbanization on the diurnal rainfall pattern in Houston: Hydrological processes. Rainfall Hydrol Proc, 2005, 19: 1089-1103.

    [49]Dixon P G, Mote T L. Patterns and causes of Atlanta’s urban heat island-initiated precipitation. J Appl Meteor, 2003, 42: 1273-1284.

    [50]Molders N, Olson M A. Impact of urban e ff ects on precipitation in

    high latitudes. J Hydrometeor, , 2004, 5: 409-429.

    [51]Tayan? M, Karaca M, Yenigün O. Annual and seasonal air temperature trend patterns of climate change and urbanization e ff ects in relation to air pollutants in Turkey. J Geophys Res, 1997, 102 (D2): 1909-1920.

    [52]Rosenfeld D. Suppression of rain and snow by urban and industrial air pollution. Science, 2000, 287: 1793-1796.

    [53]Van den Heever S C, Cotton W R. Urban aerosol impacts on downwind convective storms. J Appl Meteor Climatol, 2007, 46(6): 828-850.

    [54]苗世光, Chen F, 李青春, 等. 北京城市化對夏季大氣邊界層結(jié)構(gòu)及降水的月平均影響. 地球物理學(xué)報, 2010, 53(7): 1580-1593.

    [55]王喜全, 王自發(fā), 齊彥斌, 等. 城市化進(jìn)程對北京地區(qū)冬季降水分布的影響. 中國科學(xué)D輯: 地球科學(xué), 2008, 38(11): 1438-1443.

    [56]吳風(fēng)波, 湯劍平. 城市化對2008年8月25日上海一次特大暴雨的影響. 南京大學(xué)學(xué)報(自然科學(xué)), 2011, 47(1): 71-81.

    [57]Guo X, Fu D, Wang J. Mesoscale convective precipitation system modified by urbanization in Beijing City. Atmos Res, 2006, 82: 112-126.

    [58]Jauregui E, Romales E. Urban e ff ects on convective precipitation in Mexico City. Atmos Environ, 1996, 30: 3383-3389.

    [59]Bornstein R, Lin Q. Urban heat islands and summertime convective thunderstorms in Atlanta: Three cases studies. Atmos Environ, 2000, 34: 507-516.

    [60]Baik J J, Kim Y H, Chun H Y. Dry and moist convection forced by an urban heat island. J Appl Meteor, 2001, 40: 1462-1475.

    [61]Changnon S A, Westcott N E. Heavy rainstorms in Chicago: Increasing frequency, altered impacts, and future implications. J Amer Water Res Assoc, 2002, 38: 1467-1475.

    [62]Craig K, Bornstein R. MM5 simulation of urban induced convective precipitation over Atlanta. 4th Conf on the Urban Environ, Amer Meteor Soc, Norfolk, VA, 2002.

    [63]Rozoff C, Cotton W R, Adegoke J O. Simulation of St. Louis, Missouri, land use impacts on thunderstorms. J Appl Meteor, 2003, 42: 716-738.

    [64]Chen T C, Wang S Y, Yen M C. Enhancement of afternoon thunderstorm activity by urbanization in a valley: Taipei. J Appl Meteor Climatol, 2007, 46(9): 1324-1340.

    [65]Miao S, Chen F. Formation of horizontal convective rolls in urban areas. Atmos Res, 2008, 89: 298-304.

    [66]Shem W, Shepherd M. On the impact of urbanization on summertime thunderstorms in Atlanta: Two numerical model case studies. Atmos Res, 2009, 92: 172-189.

    [67]蒙偉光, 閆敬華, 扈海波. 城市化對珠江三角洲強(qiáng)雷暴天氣的可能影響. 大氣科學(xué), 2007, 31(2): 364-376.

    [68]江曉燕, 劉偉東. 從不同的陸面資料看城市化對北京強(qiáng)降水的影響. 氣象學(xué)報, 2006, 64(4): 527-536.

    [69]Bentley M L, Ashley W S, Stallins J A. Climatological radar delineation of urban convection for Atlanta, Georgia. Int J Climatol, 2009, 30(11): 1589-1594.

    [70]Niyogi D, Pyle P, Lei M, et al. Urban modification of thunderstorms: An observational storm climatology and model case study for the Indianapolis urban region. J Appl Meteor Climatol, 2011, 50(5): 1129-1144.

    [71]Ntelekos A A, Smith J A, Krajewski W F. Climatological analyses of thunderstorms and fl ash fl oods in the Baltimore metropolitan region. J Hydrometeor, 2007, 8(1): 88-101.

    [72]Mueller C, Saxen T, Roberts R, et al. NCAR Auto Nowcast System. Wea Forecasting, 2003, 18(4): 545-561.

    [73]王令, 丁青蘭, 陳明軒, 等. 北京地區(qū)夏季對流風(fēng)暴的雷達(dá)氣候初步研究. 中國氣象學(xué)會雷達(dá)氣象學(xué)與氣象雷達(dá)委員會第二屆學(xué)術(shù)年會文集, 2006.

    [74]孫繼松. 北京地區(qū)夏季邊界層急流的基本特征及形成機(jī)理研究.大氣科學(xué), 2005, 29(3): 445-452.

    [75]王婷婷, 王迎春, 陳明軒, 等. 北京地區(qū)干濕雷暴形成機(jī)制的對比分析. 氣象, 2011, 37(2): 142-155.

    [76]陳雙, 王迎春, 張文龍, 等. 復(fù)雜地形下雷暴增強(qiáng)過程的個例研究. 氣象, 2011, 37(7): 802-813.

    [77]陳明軒, 王迎春. 低層垂直風(fēng)切變和冷池相互作用影響華北地區(qū)一次颮線過程發(fā)展維持的數(shù)值模擬. 氣象學(xué)報, 2012, 70(3): 371-386.

    [78]陳明軒, 王迎春, 肖現(xiàn), 等. 基于雷達(dá)資料四維變分同化和三維云模式對一次超級單體風(fēng)暴發(fā)展維持熱動力機(jī)制的模擬分析. 大氣科學(xué), 2012, 36(5): 929-944.

    [79]孫繼松, 何娜, 郭銳, 等. 多單體雷暴的形變與列車效應(yīng)傳播機(jī)制. 大氣科學(xué), 2013, 37(1): 137-148.

    [80]陳明軒, 王迎春, 肖現(xiàn), 等. 北京7·21暴雨雨團(tuán)的發(fā)生和傳播機(jī)理.氣象學(xué)報, 2013, 71(4): 569-592.

    [81]Kober K, Craig G C, Keil C, et al. Blending a probabilistic nowcasting method with a high-resolution numerical weather prediction ensemble for convective precipitation forecasts. Q J R Meteor Soc, 2012, 138: 755-768.

    [82]程叢蘭, 陳明軒, 王建捷, 等. 基于雷達(dá)外推臨近預(yù)報和中尺度數(shù)值預(yù)報融合技術(shù)的短時定量降水預(yù)報試驗. 氣象學(xué)報, 2013, 71(3): 397-415.

    [83]Chang P L, Lin P F, Jou B J D, et al. An application of re fl ectivity climatology in constructing radar hybrid scans over complex terrain. J Atmos Oceanic Technol, 2009, 26(7): 1315-1327.

    [84]Overeem A, Holleman I, Buishand A. Derivation of a 10-year radar-based climatology of rainfall. J Appl Meteor Climatol, 2009, 48(7): 1448-1463.

    [85]Hocker J E, Basara J B. A geographic information systems-Based analysis of supercells across Oklahoma from 1994 to 2003. J Appl Meteor Climatol, 2008, 47(5): 1518-1538.

    [86]Parker M D, Knievel J C. Do meteorologists suppress thunderstorms? Radar-derived statistics and the behavior of moist convection. Bull Amer Meteor Soc, 2005, 86(3): 341-358.

    An Overview of Progresses in Radar Climatology and Its Prospective Applications in Nowcasting Severe Weather over Urban Regions

    Chen Mingxuan1, Wang Yingchun2, Gao Feng1, Xiao Xian1, Cheng Conglan1
    (1 Institute of Urban Meteorology, China Meteorological Administration, Beijing 100089 2 Beijing Meteorological Service, Beijing 100089)

    International progress in radar climatology in recent years is introduced fi rst in this paper, focusing on diurnal variations of convective storms and precipitation in different regions in the world based on radar climatology, and some research fi ndings of radar climatology for diurnal variations of warm-season convective storms over contiguous North China (Beijing and its vicinity). Studies of radar climatology mainly reveal initiation, enhancement, and diurnal cycle and propagation of convection over different regions, which has a close relationship with solar heating, topographical forcing, prevailing winds, as well as synoptic and mesoscale systems. Statistics of identif i ed and tracked 3D storm properties are another especial study of radar climatology, and are also presented. Studies and progress in effects of urbanization on convective storms and precipitation, and studies on convective weather over urban regions based on radar climatology are also reviewed. Finally, application prospects of radar climatology in analysis on convective mechanism and fi ne-scale nowcasting of convective storms and precipitation over urban regions are discussed in the paper.

    radar climatology, convection, precipitation, urbanization, nowcasting

    10.3969/j.issn.2095-1973.2014.05.003

    2013年5月30日;

    2013年9月5日

    陳明軒(1973—),Email: mxchen@ium.cn

    資助信息:國家自然科學(xué)基金(41075036);公益性行業(yè)(氣象)科研專項(GYHY201306008)

    猜你喜歡
    風(fēng)暴對流氣候
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    瞧,氣候大不同
    氣候變暖會怎樣?
    《風(fēng)暴中奔跑的樹》
    散文詩(2017年12期)2017-08-23 16:07:44
    玩轉(zhuǎn)腦風(fēng)暴
    可怕的大風(fēng)暴
    搖籃(2016年13期)2016-02-28 20:53:17
    基于ANSYS的自然對流換熱系數(shù)計算方法研究
    2015A/W暗黑風(fēng)暴來襲!
    Coco薇(2015年10期)2015-10-19 12:31:52
    二元驅(qū)油水界面Marangoni對流啟動殘余油機(jī)理
    立冬
    中文字幕av成人在线电影| 欧美色视频一区免费| 亚洲av免费高清在线观看| 99久久精品国产国产毛片| 69人妻影院| 搡老妇女老女人老熟妇| 淫秽高清视频在线观看| 久久精品国产亚洲av香蕉五月| 村上凉子中文字幕在线| 国产一区二区三区av在线 | 波多野结衣巨乳人妻| 亚洲人成网站在线观看播放| 国产精品电影一区二区三区| 真人做人爱边吃奶动态| 搡老熟女国产l中国老女人| 久久久久国产网址| 三级经典国产精品| 亚洲第一区二区三区不卡| 大香蕉久久网| 国产精品一区二区三区四区免费观看 | 欧美丝袜亚洲另类| 欧美xxxx性猛交bbbb| 日本一二三区视频观看| 好男人在线观看高清免费视频| 亚洲国产色片| 午夜精品在线福利| 别揉我奶头~嗯~啊~动态视频| 少妇人妻精品综合一区二区 | 亚洲av熟女| 最近中文字幕高清免费大全6| 99国产精品一区二区蜜桃av| 免费观看人在逋| 岛国在线免费视频观看| 男女视频在线观看网站免费| 99国产精品一区二区蜜桃av| 男插女下体视频免费在线播放| 国产精品久久久久久av不卡| 插阴视频在线观看视频| 国内精品宾馆在线| 可以在线观看的亚洲视频| 成人性生交大片免费视频hd| 午夜激情欧美在线| 国产单亲对白刺激| 久久久a久久爽久久v久久| 免费不卡的大黄色大毛片视频在线观看 | 麻豆一二三区av精品| 日韩欧美在线乱码| 搞女人的毛片| 精品不卡国产一区二区三区| 两个人的视频大全免费| 91午夜精品亚洲一区二区三区| 日日啪夜夜撸| 国产不卡一卡二| 久久久久久久久久黄片| 淫秽高清视频在线观看| 99国产精品一区二区蜜桃av| 一级毛片我不卡| 简卡轻食公司| 又爽又黄无遮挡网站| 国内久久婷婷六月综合欲色啪| 青春草视频在线免费观看| 精品少妇黑人巨大在线播放 | 国产真实伦视频高清在线观看| 婷婷精品国产亚洲av在线| 国产精品无大码| 午夜a级毛片| 国产精品一及| 天堂动漫精品| 精品久久久久久久久av| 观看美女的网站| av卡一久久| 成人无遮挡网站| 岛国在线免费视频观看| 69av精品久久久久久| 给我免费播放毛片高清在线观看| 欧美成人精品欧美一级黄| 天堂影院成人在线观看| 国产色爽女视频免费观看| 亚洲人成网站在线观看播放| 国内精品一区二区在线观看| 亚洲第一区二区三区不卡| 2021天堂中文幕一二区在线观| 精品乱码久久久久久99久播| 免费大片18禁| 18禁裸乳无遮挡免费网站照片| 精品无人区乱码1区二区| 高清毛片免费看| 真实男女啪啪啪动态图| 成人欧美大片| 成人精品一区二区免费| 国产精品国产高清国产av| 亚洲中文字幕日韩| av天堂在线播放| 一a级毛片在线观看| 波多野结衣高清作品| 少妇猛男粗大的猛烈进出视频 | 国产男人的电影天堂91| 亚洲最大成人av| 69人妻影院| 好男人在线观看高清免费视频| 欧美成人精品欧美一级黄| 99热精品在线国产| 亚洲最大成人中文| 成人毛片a级毛片在线播放| 色噜噜av男人的天堂激情| 亚洲国产欧洲综合997久久,| 一夜夜www| 久久久国产成人免费| 欧美成人免费av一区二区三区| 免费观看精品视频网站| 欧美不卡视频在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 一区二区三区四区激情视频 | 国产精品伦人一区二区| 久久人妻av系列| 国产亚洲精品综合一区在线观看| 波野结衣二区三区在线| 最好的美女福利视频网| 亚洲av.av天堂| 国产男人的电影天堂91| 国产精品精品国产色婷婷| 插逼视频在线观看| 国产老妇女一区| 亚洲乱码一区二区免费版| 免费在线观看影片大全网站| 国产精品免费一区二区三区在线| 人妻少妇偷人精品九色| 99久久久亚洲精品蜜臀av| 国产高清有码在线观看视频| 亚洲天堂国产精品一区在线| 午夜精品国产一区二区电影 | 国产白丝娇喘喷水9色精品| 中文字幕人妻熟人妻熟丝袜美| 在线免费十八禁| 你懂的网址亚洲精品在线观看 | 精品午夜福利视频在线观看一区| av在线播放精品| 日本爱情动作片www.在线观看 | 日本在线视频免费播放| 国产精品一区二区免费欧美| 老师上课跳d突然被开到最大视频| 亚洲图色成人| 国产乱人偷精品视频| 51国产日韩欧美| 夜夜爽天天搞| 亚洲av成人av| 国产精品一区二区性色av| 亚洲精品一卡2卡三卡4卡5卡| 尾随美女入室| 国产免费男女视频| 国产免费男女视频| 黄色日韩在线| 成人高潮视频无遮挡免费网站| 中文字幕久久专区| 又黄又爽又刺激的免费视频.| 特大巨黑吊av在线直播| 简卡轻食公司| 麻豆av噜噜一区二区三区| 韩国av在线不卡| 可以在线观看的亚洲视频| 69av精品久久久久久| av卡一久久| 校园人妻丝袜中文字幕| 国产精品三级大全| 国产伦精品一区二区三区视频9| 午夜福利18| 精品久久久久久成人av| 美女 人体艺术 gogo| 午夜爱爱视频在线播放| 精品一区二区三区视频在线| 国内精品宾馆在线| 亚洲精品影视一区二区三区av| 偷拍熟女少妇极品色| av在线天堂中文字幕| 91久久精品国产一区二区三区| 热99在线观看视频| 午夜a级毛片| 国产 一区精品| 91av网一区二区| 国产午夜精品久久久久久一区二区三区 | 免费人成视频x8x8入口观看| 国产精品嫩草影院av在线观看| 日韩高清综合在线| 老司机影院成人| 国产v大片淫在线免费观看| 女人被狂操c到高潮| 午夜福利成人在线免费观看| 男人狂女人下面高潮的视频| 干丝袜人妻中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 亚洲图色成人| 国产av麻豆久久久久久久| 国产午夜精品论理片| 别揉我奶头 嗯啊视频| 国产精品电影一区二区三区| 又黄又爽又免费观看的视频| 日日摸夜夜添夜夜添av毛片| 网址你懂的国产日韩在线| .国产精品久久| 少妇熟女aⅴ在线视频| 寂寞人妻少妇视频99o| 午夜a级毛片| 亚洲av免费在线观看| 欧美精品国产亚洲| 欧美精品国产亚洲| 免费人成视频x8x8入口观看| 麻豆成人午夜福利视频| 成人鲁丝片一二三区免费| 国产aⅴ精品一区二区三区波| 麻豆一二三区av精品| 小说图片视频综合网站| av在线蜜桃| 精品一区二区三区人妻视频| 国产精品一区二区性色av| 此物有八面人人有两片| 国产亚洲精品久久久久久毛片| 欧美又色又爽又黄视频| 美女大奶头视频| 亚洲av美国av| 美女黄网站色视频| 成年女人永久免费观看视频| 日韩欧美免费精品| 久久天躁狠狠躁夜夜2o2o| videossex国产| 日韩亚洲欧美综合| 欧美中文日本在线观看视频| 色av中文字幕| 国产伦精品一区二区三区四那| 日韩欧美在线乱码| 亚洲色图av天堂| 欧美三级亚洲精品| 22中文网久久字幕| av天堂在线播放| 日韩精品有码人妻一区| 亚洲内射少妇av| 淫秽高清视频在线观看| 免费观看人在逋| 六月丁香七月| 简卡轻食公司| 级片在线观看| 波多野结衣巨乳人妻| 欧美激情国产日韩精品一区| 成人av一区二区三区在线看| 亚洲国产高清在线一区二区三| 国产v大片淫在线免费观看| 国产伦一二天堂av在线观看| 最近最新中文字幕大全电影3| 国产69精品久久久久777片| 日本撒尿小便嘘嘘汇集6| 亚洲国产色片| 亚洲美女搞黄在线观看 | 久久久久久久久久黄片| 日韩欧美一区二区三区在线观看| 午夜a级毛片| 成年女人永久免费观看视频| 别揉我奶头~嗯~啊~动态视频| 尤物成人国产欧美一区二区三区| 少妇的逼水好多| 亚洲在线观看片| 少妇熟女aⅴ在线视频| 欧美国产日韩亚洲一区| 美女被艹到高潮喷水动态| 看免费成人av毛片| 日韩欧美精品v在线| 久久这里只有精品中国| a级毛片a级免费在线| 99国产精品一区二区蜜桃av| 在现免费观看毛片| 能在线免费观看的黄片| 人妻久久中文字幕网| 精品久久久久久久久久久久久| 亚洲人成网站在线观看播放| 黄色日韩在线| 国产不卡一卡二| 日本精品一区二区三区蜜桃| 国产成人福利小说| 一个人免费在线观看电影| av视频在线观看入口| 成人性生交大片免费视频hd| 国产精品久久久久久久久免| 天堂影院成人在线观看| av在线亚洲专区| 永久网站在线| 搡女人真爽免费视频火全软件 | 成人无遮挡网站| 久久99热6这里只有精品| 亚洲成av人片在线播放无| 久久久久久九九精品二区国产| 丰满人妻一区二区三区视频av| 亚洲精品久久国产高清桃花| 亚洲精品乱码久久久v下载方式| 色综合色国产| 色噜噜av男人的天堂激情| 精品久久久久久久末码| 国产又黄又爽又无遮挡在线| 99视频精品全部免费 在线| 国产高清视频在线观看网站| 午夜激情欧美在线| 中国美白少妇内射xxxbb| 欧美日韩精品成人综合77777| 午夜福利在线观看吧| 十八禁国产超污无遮挡网站| 嫩草影视91久久| 国产精品人妻久久久影院| 美女 人体艺术 gogo| 久久久久国产网址| 菩萨蛮人人尽说江南好唐韦庄 | 久久亚洲国产成人精品v| 18禁黄网站禁片免费观看直播| 日本爱情动作片www.在线观看 | 少妇被粗大猛烈的视频| av在线播放精品| 久久精品国产亚洲网站| 日韩精品中文字幕看吧| 久久久久久伊人网av| 久久精品夜色国产| 国产精品99久久久久久久久| 国产一区二区三区av在线 | 在线观看午夜福利视频| 亚洲人成网站高清观看| 最好的美女福利视频网| 91久久精品国产一区二区成人| 亚洲丝袜综合中文字幕| 日日啪夜夜撸| 最近中文字幕高清免费大全6| 色尼玛亚洲综合影院| 成人av一区二区三区在线看| 亚洲av一区综合| 午夜福利成人在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品一及| 亚洲欧美中文字幕日韩二区| 国产精品人妻久久久久久| 国产精品无大码| 久久综合国产亚洲精品| 欧美成人免费av一区二区三区| 欧美区成人在线视频| 久久久精品欧美日韩精品| 国产精品久久久久久久久免| 在线天堂最新版资源| 露出奶头的视频| 中出人妻视频一区二区| 国产精品久久久久久久电影| 亚洲精品粉嫩美女一区| 非洲黑人性xxxx精品又粗又长| av福利片在线观看| 免费高清视频大片| 免费看av在线观看网站| 波野结衣二区三区在线| 国产私拍福利视频在线观看| 国产黄色视频一区二区在线观看 | 亚洲欧美清纯卡通| 日韩大尺度精品在线看网址| 亚洲激情五月婷婷啪啪| 亚洲五月天丁香| 国产在视频线在精品| 人人妻人人澡欧美一区二区| 成人无遮挡网站| 国产亚洲精品av在线| 久久久久久久久大av| 免费看av在线观看网站| 欧美区成人在线视频| 亚洲最大成人av| 亚洲第一电影网av| 老熟妇仑乱视频hdxx| 黑人高潮一二区| 亚洲欧美精品综合久久99| 国产极品精品免费视频能看的| 国产欧美日韩精品亚洲av| 国产单亲对白刺激| 精品福利观看| 听说在线观看完整版免费高清| 日本三级黄在线观看| 亚洲真实伦在线观看| 99久久精品国产国产毛片| 午夜日韩欧美国产| 搡女人真爽免费视频火全软件 | 国产伦精品一区二区三区四那| 色噜噜av男人的天堂激情| 成人av在线播放网站| 国产精品乱码一区二三区的特点| 悠悠久久av| 麻豆国产97在线/欧美| 老司机影院成人| 人妻制服诱惑在线中文字幕| 天堂动漫精品| 日本撒尿小便嘘嘘汇集6| 在线天堂最新版资源| 69av精品久久久久久| 人人妻人人看人人澡| 美女黄网站色视频| 欧美极品一区二区三区四区| 亚洲中文字幕日韩| 午夜亚洲福利在线播放| 国产高清三级在线| 日本五十路高清| 免费搜索国产男女视频| 又黄又爽又刺激的免费视频.| 男人舔奶头视频| 久久九九热精品免费| 国产探花极品一区二区| 国产黄a三级三级三级人| 色视频www国产| 日本成人三级电影网站| 国产精品三级大全| 国产av一区在线观看免费| 免费看a级黄色片| 亚洲最大成人手机在线| 国产午夜精品论理片| 在线免费观看的www视频| 久久精品国产99精品国产亚洲性色| 亚洲第一区二区三区不卡| 精品免费久久久久久久清纯| a级一级毛片免费在线观看| 久久久久久伊人网av| 午夜福利成人在线免费观看| 国产高潮美女av| 99国产精品一区二区蜜桃av| 日日摸夜夜添夜夜添小说| 一区福利在线观看| 精品一区二区三区视频在线| 国产精品一区二区三区四区免费观看 | 国产成年人精品一区二区| 亚洲欧美日韩高清专用| 国产成人一区二区在线| 亚洲精品粉嫩美女一区| 日本在线视频免费播放| 午夜福利在线在线| 欧美色视频一区免费| 亚洲av五月六月丁香网| 亚洲最大成人av| 乱码一卡2卡4卡精品| 一a级毛片在线观看| 秋霞在线观看毛片| 亚洲成人av在线免费| 欧美一区二区亚洲| 亚洲不卡免费看| 免费看a级黄色片| 日本一二三区视频观看| 少妇高潮的动态图| 国产精品精品国产色婷婷| 亚洲精品乱码久久久v下载方式| 免费看a级黄色片| 男人和女人高潮做爰伦理| 久久热精品热| 日韩欧美免费精品| 亚洲天堂国产精品一区在线| 国产一区二区三区av在线 | 久久中文看片网| 春色校园在线视频观看| 少妇熟女aⅴ在线视频| 日韩强制内射视频| 99精品在免费线老司机午夜| avwww免费| 91精品国产九色| 此物有八面人人有两片| 麻豆成人午夜福利视频| 亚洲第一区二区三区不卡| videossex国产| 国产成人freesex在线 | 啦啦啦观看免费观看视频高清| 99热网站在线观看| 国产精品久久久久久久久免| 婷婷精品国产亚洲av在线| 免费观看的影片在线观看| 亚洲经典国产精华液单| 国产精品爽爽va在线观看网站| 国产探花在线观看一区二区| 国产精品免费一区二区三区在线| 欧美极品一区二区三区四区| 97碰自拍视频| 99在线视频只有这里精品首页| 成人高潮视频无遮挡免费网站| 超碰av人人做人人爽久久| 午夜免费男女啪啪视频观看 | 国产探花在线观看一区二区| 亚州av有码| 观看免费一级毛片| 久久精品夜夜夜夜夜久久蜜豆| 女人被狂操c到高潮| 亚洲精品日韩在线中文字幕 | 秋霞在线观看毛片| 男人舔奶头视频| 最近的中文字幕免费完整| 亚洲第一区二区三区不卡| 日韩欧美精品免费久久| 成人亚洲欧美一区二区av| 色播亚洲综合网| 搡老熟女国产l中国老女人| 精品一区二区免费观看| 午夜福利成人在线免费观看| av中文乱码字幕在线| 午夜激情福利司机影院| 永久网站在线| 禁无遮挡网站| 精品福利观看| 嫩草影院入口| 黄色日韩在线| 免费搜索国产男女视频| 人妻制服诱惑在线中文字幕| 成人漫画全彩无遮挡| 一级毛片aaaaaa免费看小| 国产亚洲91精品色在线| 国产亚洲欧美98| 日韩一本色道免费dvd| 1024手机看黄色片| 欧美区成人在线视频| 日韩欧美精品免费久久| 人人妻人人澡人人爽人人夜夜 | 91久久精品国产一区二区三区| 99热这里只有是精品50| 久久久久久久久大av| 国产色爽女视频免费观看| 国产精品免费一区二区三区在线| 丝袜喷水一区| 免费不卡的大黄色大毛片视频在线观看 | 男人和女人高潮做爰伦理| 亚洲电影在线观看av| 亚洲精品久久国产高清桃花| 三级男女做爰猛烈吃奶摸视频| 女人十人毛片免费观看3o分钟| 色综合色国产| 天天躁日日操中文字幕| 欧美不卡视频在线免费观看| av视频在线观看入口| 久久鲁丝午夜福利片| 一级毛片久久久久久久久女| 特级一级黄色大片| 日日干狠狠操夜夜爽| 亚洲色图av天堂| 精品久久久久久久末码| 大香蕉久久网| 99久久无色码亚洲精品果冻| 亚洲色图av天堂| 一级毛片我不卡| 国产成人精品久久久久久| 午夜爱爱视频在线播放| 日韩一本色道免费dvd| 精品无人区乱码1区二区| 亚洲精品日韩在线中文字幕 | 麻豆精品久久久久久蜜桃| av在线观看视频网站免费| 久久精品国产亚洲网站| 国产三级中文精品| 麻豆精品久久久久久蜜桃| 亚洲精品日韩av片在线观看| 亚洲av电影不卡..在线观看| 日本 av在线| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 国产v大片淫在线免费观看| 久久精品久久久久久噜噜老黄 | 女生性感内裤真人,穿戴方法视频| 99精品在免费线老司机午夜| 精品久久国产蜜桃| 赤兔流量卡办理| 亚洲最大成人av| 欧美高清成人免费视频www| 我要看日韩黄色一级片| 欧美日韩国产亚洲二区| 校园人妻丝袜中文字幕| 成人高潮视频无遮挡免费网站| 女人十人毛片免费观看3o分钟| 亚洲经典国产精华液单| 一本久久中文字幕| 两个人的视频大全免费| 能在线免费观看的黄片| 婷婷精品国产亚洲av在线| 干丝袜人妻中文字幕| 婷婷色综合大香蕉| 日本在线视频免费播放| 97超视频在线观看视频| 免费观看精品视频网站| 露出奶头的视频| 国产成人一区二区在线| 亚洲精品一卡2卡三卡4卡5卡| 午夜日韩欧美国产| 国产精品爽爽va在线观看网站| 久久久久久久久久成人| 日日摸夜夜添夜夜爱| 99视频精品全部免费 在线| 在现免费观看毛片| 1000部很黄的大片| 桃色一区二区三区在线观看| 亚洲国产精品sss在线观看| 精品日产1卡2卡| 国产蜜桃级精品一区二区三区| 最近最新中文字幕大全电影3| 国产黄a三级三级三级人| 午夜视频国产福利| 舔av片在线| 亚洲最大成人中文| 成人无遮挡网站| 91久久精品电影网| 麻豆精品久久久久久蜜桃| 亚洲精品成人久久久久久| 国产av一区在线观看免费| 长腿黑丝高跟| 波多野结衣巨乳人妻| 日本爱情动作片www.在线观看 | 日韩欧美三级三区| 国产色婷婷99| a级毛片a级免费在线| 日韩,欧美,国产一区二区三区 | 亚洲精品亚洲一区二区| 一级黄片播放器| 国产精品无大码| 亚洲图色成人| 亚洲av成人av| 日韩欧美在线乱码| av视频在线观看入口| 成人特级黄色片久久久久久久| 人妻少妇偷人精品九色| 午夜福利18| 又爽又黄a免费视频|