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

    北京一次積層混合云系結(jié)構(gòu)和水分收支的數(shù)值模擬分析

    2015-09-22 05:33:58陶玥李軍霞黨娟李宏宇孫晶
    大氣科學 2015年3期
    關(guān)鍵詞:云系水汽含水量

    陶玥 李軍霞 黨娟 李宏宇 孫晶

    1中國氣象科學研究院中國氣象局人工影響天氣中心,北京100081

    2山西省人工降雨防雹辦公室,太原030002

    1 引言

    積層混合云是我國主要降水云型,同時是一種非常有效地降水系統(tǒng),在冷鋒、梅雨鋒、冷渦、地形云等多種天氣系統(tǒng)中都會出現(xiàn)。積層混合云也是環(huán)北京地區(qū)一種主要的降水云系(王昂生等,1987;連志鸞和李國翠,2005),是我國北方實施人工增雨作業(yè)的主要對象。結(jié)合觀測和數(shù)值模擬結(jié)果深入分析積層混合云系宏微觀結(jié)構(gòu)、水分收支、降水效率,可為人工影響天氣工作提供理論支持,對環(huán)北京地區(qū)的人工增雨工作有重要意義。

    過去對積層混合云降水的研究多集中在梅雨鋒、東北冷渦等暴雨天氣,對積層混合云降水發(fā)生發(fā)展的天氣學過程和降水機理研究較多,對其云微物理結(jié)構(gòu)、云系中各種水分收支狀況和降水效率研究較少。少數(shù)觀測資料分析和數(shù)值模擬研究指出積層混合云系中微物理量在水平和垂直方向上分布都不均勻,積云中的液態(tài)水垂直積分含量大大高于層云中含量;在降水物理過程中,積層混合云不僅具有積云和層云形態(tài)混合的特征,還具有冷云過程和暖云過程共存的相態(tài)混合特征(鄒倩等,2008;于翡和姚展予,2009;周黎明等,2010;張佃國等,2011)。

    目前,國內(nèi)外許多研究人員對大氣水分收支方面開展了大量研究。水汽輸送狀況的改變可造成區(qū)域的水汽含量變化及區(qū)域水汽收支平衡的變化,并最終導致降水的變化。大氣中水汽輸送的重要性已為廣大學者所公認,丁一匯和胡國權(quán)(2003)、胡國權(quán)等(1991)對 1998、1991 年我國江淮等地持續(xù)暴雨期的水汽收支、循環(huán)特征進行了定量計算和分析。丁一匯等指出中國東部暴雨區(qū)的水汽主要是從南邊界和西邊界流入,東邊界和北邊界流出,且水汽的流入流出主要集中在中低層。馬開玉和高國棟(1992)也指出,長江流域典型旱澇年夏季,我國大氣中水汽輸送的特征明顯不同,水汽輸送通量的分布在大范圍內(nèi)與降水的分布和降水的特點是一致的。平凡和羅哲賢(2007)診斷分析了熱帶地區(qū)對流活動的水汽收支,指出深對流與淺對流在熱量及水汽循環(huán)中存在較大差異,深對流中水汽的凝結(jié)及潛熱釋放起著主導作用。Tao et al.(1983)分析了颮線系統(tǒng)的水分收支,指出對流區(qū)和層狀區(qū)對于地面降水都有重要作用,對流區(qū)在產(chǎn)生層狀區(qū)降水方面起主要作用,赤道與中緯度的颮線系統(tǒng)的水分收支間存在差別。降水效率是研究水分收支常用的物理量。

    降水效率的定義有多種。Braham(1952)首次把降水效率(PE)定義為降水量與有假絕熱上升凝結(jié)的水量之比。目前,國際上降水效率的定義主要分為三類:雨水微物理降水效率(RMPE)、云微物理降水效率(CMPE)和大尺度降水效率(LSPE)(Gao and Li,2011;Shen et al.,2013)。RMPE 定義為雨強與雨水微物理收支中的降雨源項之比(Auer and Marwitz,1968;Heymsfield and Schotz,1985;Chong and Hauser,1989;Doswell et al., 1996;Ferrier et al., 1996;Li et al., 2002;Tao et al., 2004;Sui et al., 2005)。由于數(shù)值模式中云微物理參數(shù)化方案廣泛用于云微物理量的預報,降水效率也可以通過云微物理收支來定義,即,CMPE為雨強與云微物理收支中的降雨源項之比,云微物理收支源于云中水凝物(云水、冰晶、雪、霰)的微物理過程(Weisman and Klemp,1982;Lipps and Hemler,1986;Ferrier et al., 1996;Li et al., 2002;Sui et al.,2005)。LSPE定義為雨強與地面降雨收支中的降雨源項之比(Sui et al., 2007)。由于降雨源項在不同收支中不同,導致相同雨強下計算的降水效率也不同(Gao and Li,2011;Shen et al., 2013)。Colle et al.(2005)研究了迎風坡的微物理收支狀況和各云微物理過程的重要性,指出迎風坡?lián)p耗水汽最多,這與水汽凝結(jié)核雪凝華過程有關(guān)。在高層產(chǎn)生的雪有2/3流入到背風坡,導致迎風坡降水效率只有50%。Hobbs et al.(1980)把降水效率定義為垂直于雨帶1 m寬距離的總雨強與總凝結(jié)率之比,通過對冷鋒云系的研究得出,暖區(qū)和窄冷鋒雨帶的降水效率分別在40%~50%和30%~50%,高層淺對流泡里的冰粒子對寬冷鋒雨帶和一些暖區(qū)雨帶降水的發(fā)生產(chǎn)生了重要的作用。李宏宇等(2006)提出了“擬降水效率”的定義,擬降水效率定義為地面降水量與總水汽量之比,認為總水汽量由懸滯水汽總量和總濕物質(zhì)量及地面累積降水組成;分析了河南省春季低槽冷鋒降水過程,認為降水主要集中在冷鋒面附近,鋒前鋒后擬降水效率通常介于10%~30%,鋒區(qū)附近達40%~60%,積層混合區(qū)域高達70%以上。洪延超和周非非(2006)利用數(shù)值模式,分別計算了層狀云中凝結(jié)水和凝華水的降水效率,即凝結(jié)水和凝華水的降水量與雨水中的凝結(jié)水和凝華水的比率。周非非等(2010)利用數(shù)值模式,在水物質(zhì)基本達到收支平衡的條件下,估算的河南省域總水物質(zhì)降水效率、凝結(jié)率、凝華率和水凝物降水效率及水汽降水效率。

    降水云系的水資源分布特征及微物理轉(zhuǎn)化機制、降水云系區(qū)域大氣水分的收支、循環(huán)規(guī)律,對于開發(fā)利用空中云水資源具有重要意義,是研究人工增雨原理、催化作業(yè)條件和催化技術(shù)方法的重要基礎(chǔ)。人工增雨的先決條件是降水云系的存在,而水汽是成云致雨的必要條件,云水形成、冰粒子的增長都需要水汽。同時,輸入云中水汽如何轉(zhuǎn)化為地面降水,在降水形成過程中水汽如何分配和轉(zhuǎn)化?其轉(zhuǎn)化效率和降水效率如何?這是人工增雨潛力研究的重要基礎(chǔ),有助于深入了解云和降水形成的云物理過程。弄清大氣水分收支是實施人工增雨的首要問題,也是積層混合云系人工增雨潛力研究的需要。因此,本文選取2007年10月5日低槽冷鋒影響北京產(chǎn)生中雨的積層混合云系作為研究對象,利用中國氣象科學研究院(CAMS)中尺度云分辨模式,結(jié)合實測資料,分析了此次積層混合云系宏微觀結(jié)構(gòu)和降水特征,在此基礎(chǔ)上,重點分析了北京區(qū)域內(nèi)云系中各種水分收支狀況與降水效率。

    2 降水過程實況分析

    2007年10月5日20時(北京時,下同)至6日 20時,陜西北部、山西中北部、河北中部和東北部、京津的部分地區(qū)、遼寧西部等地區(qū)出現(xiàn)了一次大范圍的降水過程,部分地區(qū)出現(xiàn)中到大雨。24小時雨帶呈西南—東北走向分布,降水分布不均勻,有三個強降水中心,其中位于北京的降水中心雨量約20 mm左右,為中雨量級(圖1a)。此次降水主要出現(xiàn)在晚上至次日下午,其中山西至河北大范圍雨帶出現(xiàn)在晚上,北京地區(qū)降雨主要集中在 6日02~14 時。

    2.1 天氣背景分析

    從天氣形勢看,5日20時,500 hPa為西風槽環(huán)流形勢(圖 2a),并不斷東移,北京位于高空槽的東南部。700 hPa有槽線(圖略),在槽前為冷平流,槽后為暖平流,有西南氣流提供水汽輸送。5日20時,850 hPa上,河北中南部產(chǎn)生氣旋性環(huán)流,中心位于(38°N,112°E)。地面,東北東部為阻塞高壓,西北為貝加爾湖冷高壓,兩高壓之間為低壓輻合區(qū),有蒙古氣旋冷鋒生成(圖 2b),自西向東移動經(jīng)過北京??傊?,低槽冷鋒是此次降水過程的主要影響系統(tǒng)。

    2.2 云系的演變特征

    從衛(wèi)星亮溫TBB(Black-Body Temperature)圖(圖 3)上看,低槽冷鋒云系水平范圍較大,亮溫較高且相對均勻的云是層狀云;云系中存在多個低溫中心,云頂亮溫最低達-50°C到-60°C,說明在大范圍的層狀云中存在多個對流中心,即這次降水云系有明顯的積層混合云特征。在水平方向表現(xiàn)為不均勻分布,即,大片的層狀云中有多個絮狀對流云團。5日22時,對應(yīng)于高空槽前,陜西、山西北部至內(nèi)蒙中部上空有大片云系,呈西南—東北走向分布(圖3a),受高空槽影響,云系東移發(fā)展加強。降水云系6日01時移入北京境內(nèi),05時,降水云系中的一個對流中心覆蓋北京上空,其云頂溫度較低,云頂亮溫達-40°C(圖 3b),說明云內(nèi)有弱對流的積云,積云的云頂較高,云層較厚。隨著高空槽繼續(xù)東移,北京和河北地區(qū)云系逐漸東移減弱。結(jié)合降水演變可見,受高空槽影響,云系自西南向東北逐漸向北京推進;隨著積層混合云中對流云團移入北京并逐漸加強,北京地區(qū)的降水增強;隨后,積層混合云中對流云團東移出北京,北京區(qū)域降水減弱,云系出現(xiàn)減弱的趨勢。隨著北京北部云系消散,北京區(qū)域降水逐漸結(jié)束。對應(yīng)24 h地面降水(圖 1a),可見,此次北京地區(qū)的降水過程主

    要由高空槽云系的東移發(fā)展造成,積層混合云是產(chǎn)生此次降水的主要云系。

    圖1 2007年10月5日20時~6日20時24 h 地面降水量(單位:mm):(a)實況;(b)模擬Fig.1 (a) Observed and (b) simulated 24-hour cumulative rainfall (units: mm) from 2000 BJT (Beijing time) 5 to 2000 BJT 6 October 2007

    圖2 2007年10月5日20時的(a)500 hPa高空天氣形勢圖和(b)地面天氣圖Fig.2 (a) 500-hPa synoptic situation chart and (b) the surface synoptic chart at 2000 BJT 5 October 2007

    圖3 2007年(a)10月5日22時和(b)6日05時FY-2C衛(wèi)星的TBB(°C)分布Fig.3 Black-Body Temperature (TBB) from satellite FY-2C (units: °C) at (a) 2200 BJT 5 October 2007 and (b) 0500 BJT 6 October 2007

    3 積層混合云降水過程的數(shù)值模擬

    3.1 模式及模擬方案介紹

    本文選用中國氣象科學研究院(CAMS)中尺度云分辨模式對此次積層混合云過程進行數(shù)值模擬。該模式耦合了PSU/NCAR( Penn State University/National Center for Atmospheric Research )的非靜力平衡模式 MM5v3(Mesoscale Model)和 CAMS(Chinese Academy of Meteorological Sciences)復雜云微物理方案,此方案由樓小鳳等(2002)以胡志晉和何觀芳(1987)早期的對流云模式和層狀云模式為基礎(chǔ)發(fā)展而來,是一套準隱式格式的混合相雙參數(shù)微物理雪晶方案。方案中水物質(zhì)分別是水汽、云水、雨滴、冰晶、雪、霰,云物理預報量有11個,包括各種水成物的比質(zhì)量,雨滴、冰晶、雪、霰的比濃度以及云滴譜拓寬度。模式考慮了 31種云物理過程,包含 56個方程式。模式采用準隱式計算格式,保證了計算的正定性、穩(wěn)定性和水物質(zhì)守恒。與 MM5原有的顯式云降水方案相比, 預報量增加了雨、雪和霰的粒子比濃度,完善了多個微物理過程的描述,如云雨自動轉(zhuǎn)化、冰晶核化、冰雪和雪霰的自動轉(zhuǎn)化、冰晶繁生等過程。通過對降雪、臺風降水、華南暴雨、長江梅雨、地形云降水等實例的模擬,并與MM5原有的方案進行對比分析表明, 該方案能合理地模擬各種云降水過程,模擬的云微物理結(jié)果合理(孫晶等,2009a,2009b;劉衛(wèi)國和劉奇駿,2007a,2007b)。

    對此次低槽冷鋒云系降水過程,模擬使用了NCEP (National Centers for Environmental Prediction)/NCAR每6小時1次的1°(經(jīng)度)×1°(緯度)再分析資料作為初始場。采用三重嵌套結(jié)構(gòu),D1、D2和D3模擬域的分辨率分別為30 km、10 km、3.3 km(圖4)。D1、D2采用KF(Kain-Fritsch)對流參數(shù)化方案和顯式方案(CAMS雙參數(shù)冰相微物理方案)相結(jié)合的方法,D3不采用對流參數(shù)化方案,僅采用顯式方案。本文分析主要針對北京區(qū)域(115.3°N–117.5°N,39.3°E–41.5°E,見圖 5 中矩形區(qū))的模擬結(jié)果分析云系的水分收支。

    3.2 模擬結(jié)果的驗證

    為了驗證模擬結(jié)果的可靠性,對觀測和模擬的降水量、形勢場分布進行了對比分析。對比模擬的24小時地面降水量(圖1b)和實況(圖1a)可以看出,模擬的雨帶走向、形狀、降水落區(qū)的分布情況與實況降水基本一致,但模擬的雨帶略窄,尤其是模擬的山西北部強降水區(qū)的范圍和雨量比實況偏弱較多(實況最大雨量達到60 mm,模擬最大雨量達到20 mm)。由于數(shù)值模擬未采用雷達或衛(wèi)星的同化技術(shù),模式的初始場中不含有云方面的信息,在模擬初期6小時內(nèi)(10月05日20時~6日02時),模擬雨帶的強度偏弱、其移動速度比實況略慢,且云帶和雨帶位置比實況偏西,表現(xiàn)在6日02時6小時累積地面降水圖5a1、5b1中,觀測雨區(qū)分布在山西西北、河北中東部,山西西北大部雨量為10 mm左右,局地達30 mm左右;模擬雨區(qū)位置與實測接近,但強度偏弱,尤其是山西西北部強降水中心偏弱,最大達10 mm。從模擬和觀測每6小時降水演變來看(圖5),模式運行6小時后,模擬雨帶和實況差異減小,強降水中心的雨量和實況差異比初期6小時明顯減小。6日02時后,雨帶自西北向東南移動,降水開始移入北京,06日08時~14時,觀測和模擬的西南—東北走向雨帶形成,分布于山西北部—北京—河北中北部地區(qū),觀測中大范圍雨帶上分布2個10 mm左右降水中心,北京位于降水中心之一,而模擬降水的10 mm降水中心合并為一個,且北京 10 mm降水范圍比實況的略偏南,強度和實況相當。6日08時~6日14時,模擬和觀測的雨帶向東南移動,降水中心移出北京??梢?,模擬 24小時地面降水中山西北部強降水區(qū)的范圍和雨量偏弱的主要原因在于模式采用的初始場中沒有云方面的信息,模擬初期,模擬云帶和雨帶的強度偏弱所致。本文研究的重點區(qū)域是北京地區(qū),北京降水的主要時段內(nèi)(06日02時~14時),觀測和模擬的北京地區(qū)的降水強度和實況相當。

    圖4 三重嵌套的模擬區(qū)域(D1、D2、D3分別代表第1、2、3層模擬域)Fig.4 The triple nested simulation domains used in Chinese Academy of Meteorological Sciences (CAMS) model

    圖5 模擬(上)和實測(下)的6小時地面降水量(單位:mm):(a、d)2007年10月5日20時~6日02時;(b、e)2007年10月6日02時~08時(圖5b中矩形代表北京區(qū)域);(c、f)2007年10月6日08時~14時Fig.5 Simulated (upper) and observed (lower) 6-hour cumulative rainfall (shaded): (a, d) From 2000 BJT 5 to 0200 BJT 6 October 2007; (b, e) from 0200 BJT to 0800 BJT 6 October 2007 (the rectangle area represents the Beijing region in Fig.5b); (c, f) from 0800 BJT to 1400 BJT 6 October 2007

    總的來看,模擬的雨帶范圍和演變過程與實況基本一致,反映了雨帶向東移動這一特點。模式模擬的各時次高、中、低層位勢高度場、風場的走向和位置與同時次的實況客觀分析場(圖略)也基本吻合,說明模擬結(jié)果能較好的反映實況,模擬結(jié)果再現(xiàn)了大氣環(huán)流形勢和地面降水分布。

    4 模擬云系的結(jié)構(gòu)分布

    4.1 云系的水平分布

    為便于后面分析,首先給出模擬的北京區(qū)域平均雨強隨時間變化(圖略),可見模擬的北京區(qū)域降水主要集中在6日02時至14時,與觀測一致。降雨旺盛時段為6日04時~07時,05時的區(qū)域平均最大雨強達2.3 mm h?1。下面分析北京區(qū)域不同降水階段模擬云系的演變情況。通過對比FY-2C衛(wèi)星的TBB分布(圖3)和模擬云系,模擬云系的走向、落區(qū)和發(fā)展演變與實況云系基本一致,模擬出了云系受高空槽影響自西南向東北向北京移動,并逐漸加強,移出北京后云系逐漸減弱。模擬云系的水平結(jié)構(gòu)呈不均勻分布,即,大片的層狀云中有多個絮狀對流云團,具有明顯的積層混合云特征。

    積層混合云系中不同強度的云區(qū)移動經(jīng)過北京,引起北京地區(qū)降水強度的變化。根據(jù)北京區(qū)域平均雨強,將北京地區(qū)的降水劃分為降水初期、降水強盛期、降水末期。圖6給出了代表北京區(qū)域降水初期、降水強盛期、降水末期三個時刻模擬的云系總水凝物含水量厚度分布,分析移入北京的云系情況及其降水特征。6日01時,云系范圍較大,大部分地區(qū)的總水凝物含水量厚度分布相對比較均勻,一般在0.1~1 mm之間,大值區(qū)在山西北部,最大值為3 mm,云系前沿已移入北京西部,其上空的云層較薄,含水量較?。▓D 6a),北京西部出現(xiàn)弱降水。云系繼續(xù)向東北移動并逐漸加強,05時,云系整體范圍加大,層狀云中有多個積云云團,整個云系呈現(xiàn)積層混合云特征,表現(xiàn)為云系中的總含水量水平分布不均勻,有多個高含水量中心,其中水平范圍最大的含水量大值區(qū)覆蓋北京大部,總水凝物含水量厚度最大值達3.5 mm(圖6b),大值區(qū)是層狀云區(qū)的3倍左右,此時,北京區(qū)域降水量達到最大。隨后,云系東移,高含水量中心移出北京,北京區(qū)域降水減小,云系出現(xiàn)減弱的趨勢。14時之后,云系明顯收縮,北京北部云系消散,北京區(qū)域降水逐漸結(jié)束(圖6c)。

    由此可見,積層混合云系是導致此次北京降水的主要云系。隨著積層混合云中對流云團移入北京并逐漸加強,北京地區(qū)的降水增強;隨后,積層混合云中對流云團東移出北京,云系出現(xiàn)減弱的趨勢,北京區(qū)域降水減弱。

    4.2 云系的垂直結(jié)構(gòu)及降水機制

    利用模式細網(wǎng)格區(qū)域輸出的云水、雨水、雪、霰和冰晶比含水量,分析此次積層混合云降水過程及北京區(qū)域上空云的微物理結(jié)構(gòu)特征。

    降水初期(6日01時),沿40°N的云系垂直剖面(圖略)可見,有大范圍積層混合云系東移進入北京西部山區(qū)。在5 km高度上,有層云存在,云水含量較小,115.5°E上空有對流泡體存在,雪含量豐富,大值區(qū)在5.5 km高度上(0.3 g kg?1),雪下落到暖區(qū)融化成雨水,此時西部山區(qū)開始降雨。隨著高空槽移動,云系東移發(fā)展加強,高空雪含量增加,大值區(qū)東移,北京區(qū)域地面雨水增加。

    6日05時為北京區(qū)域降水強盛期,圖7是此時刻通過北京站(40°N)含水量的經(jīng)向—垂直剖面,北京地區(qū)(115.3°E~117.5°E)上空云層較為深厚,層狀云中有多個對流泡體存在,云頂高度達到-40°C層以上,含有冰晶和雪的云層較厚,云水層范圍大,含水量高,尤其是北京東部云水深厚且過冷云水豐富。高空雪覆蓋范圍大,比含水量也最大,存在多個雪的大值區(qū),大值(0.9 g kg?1)區(qū)在5.5 km高度上。在同時含有冰晶、雪和過冷云水的區(qū)域,形成了霰,霰的含水量中心與冰晶和雪的高含水量區(qū)對應(yīng)。冰晶雪下落撞凍過冷水,在零度層上方形成了較多的霰。雨水主要存在于云的暖區(qū),雪、霰與雨水在空間上完全匹配,雨水含水量中心的上方就是霰和雪的含水量中心,看來是雪霰粒子融化形成雨水,0°C層以下區(qū)域仍有一些未融化的雪和霰粒子。說明降水過程有冰相粒子參與,冷云過程對降水形成起到重要的作用(圖7b)。

    之后,強降水云東移并逐漸變?nèi)?,雪霰粒子含量減小,大值區(qū)東移出北京,冷云降水作用減弱,該地區(qū)降水也隨之減弱。

    為了分析云系的微觀結(jié)構(gòu)和降水機制,圖8給出了北京區(qū)域(39.3°N~41.5°N,115.3°E~117.5°E)內(nèi)水凝物含水量區(qū)域累加量(區(qū)域內(nèi)每個格點上的水凝物含水量進行累加)垂直廓線的演變。由圖可見,總體而言,按水凝物粒子含水量最大值出現(xiàn)的高度云中由高到低依次為冰晶、雪、霰、云水和雨水。6日01時,云系開始移入北京區(qū)域,云中水凝物含量都較小,但云系中雪的含量最大,其次是云水和雨水,霰和冰晶最小(圖 8a)。隨著云系移入北京,雪的含量不斷增多,隨之雨水不斷增多,到6日05時云中水凝物含量都增加,其中雪含量增加的最多。雪的含水量在500 hPa達到最大,由01時的0.8 g kg?1增加到2.55 g kg?1,同時地面雨水的含水量也達到最大(0.4 g kg?1)。在雪完全融化的高度,雨水的含水量出現(xiàn)大值,說明雪的融化對雨的形成貢獻大,即雨水主要是雪降落到 0°C層(600 hPa)以下融化形成的。這表明,冷云降水過程是北京區(qū)域雨水形成的主要過程。由于0°C層以下還存在較多的云水,0°C層附近存在少量的霰(霰的含水量比云水?。?,所以,暖云過程(凝結(jié)碰并過程)和霰的融化對雨水的形成也有貢獻,其中霰的融化對雨水的形成貢獻最?。▓D 8b)。隨后,云系變?nèi)鯑|移,雪和霰粒子含量減小,大值區(qū)東移,地面降水減弱,如:11時(圖8c),云中雪和云水含水量極大值差不多大,冰晶和霰含量較小,可見地面降水主要由冷云過程和暖云過程共同導致。與降水強盛期時(圖8b)相比,可見隨著冷云降水作用減弱,該地區(qū)降水也隨之減弱。

    結(jié)合云系的垂直剖面(圖7),圖9給出了北京區(qū)域降水強盛期(6日05時)云系中不同云區(qū)(代表北京西部、中部、東部)的云微物理結(jié)構(gòu)。從圖9可見,云系中不同云區(qū)的微物理結(jié)構(gòu)相同,均為混合相云,中高層為冰晶、雪,中層為霰,低層為云水和雨水,但水凝物粒子含水量最大值的排序不同;不同云區(qū)中雪的含量均為最大;云系中不同云區(qū)形成雨水的主要微物理過程均為雪的融化過程、霰的融化過程和凝結(jié)碰并過程,均以冷云降水過程為主,但這三個過程對形成雨水的貢獻大小不同,對應(yīng)地面雨水含量不同。如:北京西部(圖9a),在雪完全融化的高度,雨水的含水量出現(xiàn)大值;霰和云水含量較小,其中云水含量最?。徽f明雪的融化過程對雨水形成貢獻最大。北京中部(圖9b)霰的含量明顯多于北京西部和東部,中部云區(qū)中雪和霰的含量相當,其中雪含量最大值為0.8 g kg?1,霰含量最大值為0.73 g kg?1,其次是云水,云水含量最大值為0.45 g kg?1;在霰完全融化的高度,雨水的含水量出現(xiàn)大值;零度層以下幾乎不存在雪,而霰的含水量超過了0.7 g m?3,霰的融化對雨水形成有較大貢獻;此外,過冷水含量大,對雪和霰的撞凍增長有利??梢姡凑招纬捎晁呢暙I大小,冰相粒子的融化過程最大,凝結(jié)碰并過程最小,冷云過程對降水形成起到重要的作用。北京東部(圖9c)含水量由大到小依次為雪、云水、霰、雨水、冰晶;霰的含量的最大值略小于云水;說明雪的融化過程最大,凝結(jié)碰并過程次之,霰的融化過程最小。北京西部和東部地區(qū)云中霰的含量較小,霰的融化對

    雨水的形成貢獻也較小。相對于西部和東部,北京中部云區(qū)中霰的融化過程對雨水的形成貢獻較大,冷云降水過程更強,其地面雨水含量最大。

    圖6 2007年10月6日三個時刻模擬的云系總水凝物含水量厚度(單位:mm)的分布:(a)01時;(b)05時;(c)15時Fig.6 The water content (mm) of the total hydrometeors for the simulated cloud system at (a) 0100 BJT, (b) 0500 BJT, and (c) 1500 BJT on 6 October 2007

    圖7 2007年10月6日05時,通過北京站(40°N)含水量(單位:g kg?1)的經(jīng)向—垂直剖面:(a)云水(彩色陰影)、冰晶(紅色實線)含水量;(b)雪(彩色陰影)、霰(紅色實線)和雨水(綠色實線)含水量。黑線:等溫線(°C);箭頭:風場(單位:m s?1)Fig.7 Vertical sections of water content (units: g kg?1) of hydrometeors in simulated cloud along Beijing station (40°N) at 0530 BJT 6 October 2007:(a) Cloud water (shaded) and ice crystal (red solid line); (b) snow (shaded),graupel (red solid line), and rain (green solid line).Black line: isotherm(°C); arrow: wind field

    圖8 北京區(qū)域(39.3°N~41.5°N,115.3°E~117.5°E)內(nèi),2007 年 10月6日三個時刻水凝物含水量區(qū)域累加量(單位:g kg?1)的垂直廓線:(a)01 時;(b)05 時;(c)11 時。(qc:云水,qr:雨水,qg:霰,qs:雪,qi:冰晶)Fig.8 Vertical profiles of area accumulation of total water content of hydrometeors over Beijing region(39.3°N?41.5°N,115.3°E?117.5°E)for the simulated cloud system at (a) 0100 BJT, (b) 0500 BJT, and (c) 1100 BJT on 6 October 2007.(qc: cloud water, qr: rain water, qg: graupel, qs:snow, qi: ice crystal)

    圖9 2007年10月6日05時北京區(qū)域云系不同部位的水凝物含水量(單位:g kg?1)的垂直廓線:(a)(40°N, 115.5°E);(b)(40°N, 116.1°E);(c)(40°N, 116.9°E)Fig.9 Vertical profiles of total water content of hydrometeors in different places of Beijing cloud system at 0500 BJT 6 October 2007: (a) (40°N,115.5°E); (b) (40°N, 116.1°E); (c) ( 40°N, 116.9°E)

    由上可見,北京地區(qū)上空積層混合云系的微物理結(jié)構(gòu)具有混合相云特征。高層為冰晶、雪,中層為霰,低層為云水和雨水;其中,雪的含水量最高;冷云降水過程(雪和霰的融化)是雨水形成的主要過程,暖云降水過程對雨水的形成也有貢獻,其中冰相降水過程占主導地位。北京地區(qū)云系中不同云區(qū)的降水形成機制相同,均以冷云過程為降水形成的主要過程,冰相過程對降水形成起到重要的作用,但微物理過程對形成雨水的貢獻大小不同。

    5 模擬云系的水分收支分析

    5.1 水汽的源地及輸送通道

    通過對比低層等壓面上北京區(qū)域的水汽輸送,發(fā)現(xiàn)在700 hPa北京區(qū)域水汽輻合量最大,說明影響此次降水過程的水汽輸送主要集中在對流層低層。圖 10給出了代表北京區(qū)域降水初期、降水強盛期、降水末期三個時刻模擬的700 hPa水汽通量散度及風場分布圖。利用6小時1°×1° NCEP/NCAR再分析資料,圖11給出了10月6日02時700 hPa上實測的水平風場。從模擬的水汽輻合區(qū)、風場水平分布結(jié)合實測風場,可以看出,本次北京區(qū)域降水過程的主要水汽源地為黃海海面及蒙古國。一支水汽從黃海海面經(jīng)過江蘇、安徽、湖南,沿四川和陜西交界向東北輸送,另一支水汽由蒙古氣旋西側(cè)的西北氣流向東南輸送,兩支氣流在陜西北部匯合后的西南氣流將水汽輸送到華北地區(qū)。北京區(qū)域的水平水汽輸送以西、南入,東、北出為主。

    圖10 2007年10月6日三個時刻模擬的700 hPa上水汽通量散度(彩色陰影,單位:g s?1 cm?1 hPa?1;負值表示水汽的輻合)和風場 (單位:m s?1) 的水平分布:(a)01 時;(b)05 時;(c)15 時。圖中矩形區(qū)域為北京市范圍Fig.10 The water vapor flux divergence (shading, units: g s?1 cm?1 hPa?1;negative value: convergence ) and wind at 700 hPa for the simulated cloud system at (a) 0100 BJT, (b) 0500 BJT, and (c) 1500 BJT on 6 October 2007(the rectangle area represents the Beijing region)

    另外,從圖中可看出,在降水初期(圖10a),水汽輻合大值區(qū)主體位于北京區(qū)域的西側(cè),其中水平范圍較小的一個水汽輻合中心的前沿開始移入北京西南部,此時北京區(qū)域并未出現(xiàn)明顯降水。隨著天氣系統(tǒng)的發(fā)展移動,水汽輻合大值區(qū)隨之向東北伸展加強,在降水強盛期(圖10b),水汽輻合高值區(qū)基本上覆蓋了北京大部分地區(qū),其輻合中心的最大值由降水初期的 5 g s?1cm?1hPa?1增加到 10 g s?1cm?1hPa?1,可見充足的水汽輸送,為北京地區(qū)的強降水創(chuàng)造了條件。之后,隨著系統(tǒng)的進一步東移,西南向水汽輸送逐步減弱后消失,水汽通量的大值區(qū)也隨之移出北京區(qū)域,北京區(qū)域降水逐漸結(jié)束(圖10c)??梢姡麖娸椇蠀^(qū)與強降水區(qū)相對應(yīng),隨著北京地區(qū)水汽輻合量的增加,北京地區(qū)降水量開始加大。

    圖11 2007年10月6日02時700 hPa上實測的水平流場Fig.11 The observed flow field of wind at 700 hPa at 0200 BJT 6 October 2007

    5.2 水汽和水凝物的邊界輸送特征

    為了了解各邊界不同高度上水物質(zhì)輸送特點,對北京區(qū)域(115.3°N~117.5°N,39.3°E~41.5°E)各邊界不同高度上的的水物質(zhì)通量和輻合量。輻合量=(東邊界水物質(zhì)通量-西邊界水物質(zhì)通量)+(北邊界水物質(zhì)通量-南邊界水物質(zhì)通量),若大于 0,則水物質(zhì)通量是輻散(凈流出),若小于0,則水物質(zhì)通量是輻合(凈流入)。

    6日05時是北京區(qū)域降水強度較大的時次,從水汽和水凝物通量和輻合量隨高度分布可以看出(圖12),這一區(qū)域水汽輻合主要發(fā)生在800~450 hPa,最大水汽輻合量位于約700 hPa高度,水汽從西邊界、南邊界流入,東邊界、北邊界流出,水汽輻散主要在800 hPa以下,最大水汽輻散量在900 hPa,且輻散層較薄,水汽從東邊界流入,西邊界、北邊界流出。從各水凝物的通量來看,雪的通量最大,其次是云水、雨水、霰和冰晶。雪在620~650 hPa上為輻合,從西邊界、南邊界流入,從東邊界、北邊界流出,420 hPa以上為輻散,從東邊界、北邊界流出。云水在800 hPa以下為輻合,從南邊界、西邊界流入,東邊界流出。雨水輻合主要在 600~800 hPa,從南邊界、西邊界流入,從東邊界、北邊界流出,800 hPa以下為輻散,從東流入,從南邊界、北邊界、西邊界流出。冰晶在420 hPa以上的高層為輻散,從東邊界、北邊界流出,500~420 hPa為輻合,從西邊界、南邊界流入,從東邊界、北邊界流出。霰在570 hPa以下為輻合,從西邊界和南邊界流入,東邊界、北邊界流出,在550 hPa以上為輻散,主要從東邊界流出。以上分析表明,在降水較大時段,北京區(qū)域以外,水汽和水凝物主要從西邊界和南邊界輸送到域內(nèi),具有西向風分量和南向分量的暖濕氣流是云系的水分輸送帶。雪、云水存在明顯的輸送,其中雪的輻合量最大,說明平流輸送從區(qū)域外凈流入的雪量對區(qū)域內(nèi)降水的增加有一定影響。主要受高層氣流影響, 在高層冰晶、雪、霰存在明顯的輻散。

    將圖 12中各氣壓層的水物質(zhì)通量對高度作積分,可得穿過高為從地表到模式頂、長為各邊界長度的截面積的通量如圖 13所示,其正(負)值表示水物質(zhì)從北京區(qū)域水平輸送出去(進來)。5日20時至6日11時,區(qū)域整層水汽通量輻合,最大值出現(xiàn)在 6日03~05時,從各邊界通量來看,在云系發(fā)展過程中,主要從西、南邊界向區(qū)域內(nèi)輸送水汽,從東、北邊界輸出水汽。6日11時至14時,區(qū)域整層水汽通量輻散,主要從西、南邊界輸入水汽,從東、北邊界輸出水汽。結(jié)合區(qū)域平均降水量變化,北京降水區(qū)域的總體水汽收支和平均面降雨量的變化趨勢相吻合,隨著水汽收入的增加,降水量開始加大。而從水凝物各邊界通量的垂直積分隨時間演變可以看出(圖13b),這一區(qū)域總水凝物在西、南邊界上流入,到6日04時輻合量達最大。隨著東邊界和北邊界水凝物流出增大,總水凝物開始水平輻散,到6日06~08時,水凝物的輻散達到最大。北京區(qū)域降水主要時段內(nèi),水汽及水凝物通量在水平方向上為凈流入。

    從以上區(qū)域邊界輸送特征看,西邊界和南邊界是把北京區(qū)域外水汽和水凝物輸送到域內(nèi)的主要邊界,西、南邊界上水汽的主要輸入途徑分別位于高層和中低層,因而,具有西向風分量和南向分量的暖濕氣流是云系的水分輸送帶,雪、云水存在明顯的輸送,其中雪的輻合量最大,說明平流輸送導致從區(qū)域外凈流入的雪量對區(qū)域內(nèi)降水的增加有一定影響,冰晶、雪和霰受高層氣流影響存在明顯的輻散,總的水凝物表現(xiàn)為輻散,但與水汽輸送相比要小1個量級。

    5.3 水分平衡和轉(zhuǎn)化效率分析

    利用周非非等(2010)提出的水分平衡和轉(zhuǎn)化效率分析方法,由 CAMS顯式云微物理方案中水汽、云水、雨水、冰晶、雪和霰的比含水量方程可以得到水汽、總水凝物和總水物質(zhì)的平衡等式,分別如(1)~(3)式:

    圖13 穿過北京區(qū)域各邊界截面積的(a)水汽通量(qv)和輻合量的垂直積分和(b)總水凝物(云水、雨水、霰、雪和冰晶)通量和輻合量的垂直積分隨時間變化Fig.13 Time series of the vertically integrated fluxes (solid lines) and flux convergences (dashed line) of (a) water vapor (qv) and (b) all hydrometeors across each of the boundaries of the Beijing region

    表1 水分平衡等式中各項值(單位:kg)Table 1 The computational values of each term in the equilibrium equations of water substances (units: kg)

    對式(1)~(3)中各項作從10月5日20時至6日14時北京區(qū)域的時間和空間的積分,估算的值列在表 1中。其中,凝華量=冰晶凝華量+雪凝華量+霰凝華量。最終計算出式(1)等號左邊所有項的代數(shù)和為-7.7×108kg,雖然不等于0,但僅占式(1)等號左邊所有正值項之和的0.01%,基本上達到平衡。式(2)等號左邊約等于-0.72×1011kg,占式(1)等號左邊所有正值項之和的6.4%。式(3)等號左邊約等于-7.12×109kg,占式(1)等號左邊所有正值項之和的0.94%。由于這是對北京區(qū)域在整個模擬時段的時空積分的計算結(jié)果,計算的水汽和水凝物收支出現(xiàn)相對較小的偏差,是在可以接受的范圍內(nèi),故本文研究的水汽、水凝物和總水物質(zhì)的收支均基本上達到了平衡。

    在估算出水物質(zhì)基本達到平衡的前提下,根據(jù)式(1)~(3)能夠計算總水物質(zhì)降水效率、水汽凝結(jié)(華)率和水凝物降水效率:

    圖14 北京區(qū)域每小時總水物質(zhì)降水效率、凝結(jié)率、凝華率、水凝物降水效率和地面平均小時降水量(單位:mm h?1)的時間分布(注:每小時的總水物質(zhì)降水效率、凝結(jié)率和凝華率、水凝物降水效率對應(yīng)左邊y軸,地面平均小時降水量對應(yīng)右邊y軸)Fig.14 Time series of the hourly precipitation efficiency of all water substances, hourly condensation efficiency, hourly deposition efficiency,hourly precipitation efficiency of hydrometeors (the left y axis), and average hourly surface precipitation (the right y axis, units: mm h?1) in the Beijing region

    通過分別估算(4)~(6)式中各項對北京區(qū)域每小時的累積量,求得的轉(zhuǎn)化效率即為每小時的轉(zhuǎn)化效率,其中水汽(水凝物、總水物質(zhì))的初始內(nèi)留量指的是每小時初始時刻對應(yīng)的內(nèi)留量。圖14給出了每小時的總水物質(zhì)降水效率、凝結(jié)率和凝華率以及水凝物降水效率隨時間的演變。可以看出,水凝物降水效率明顯高于降水效率和凝結(jié)(華)率,前者在在6日01時至07時一直大于10%,6日05時達到最大,而后者最大也不足5%。降水效率、凝結(jié)率、凝華率及水凝物降水效率的最大值分別為 4.02%、3.01%、2.14%、26.5%。平均小時降水量的時間分布曲線與降水效率、凝結(jié)率和凝華率的曲線較相似,平均雨強的最大值出現(xiàn)在 6日 05時。根據(jù)(4)~(6)式估算在 5日 20時~6日14時北京區(qū)域的時空累積量,則總體上算得的總水物質(zhì)降水效率、凝結(jié)率、凝華率及總水凝物降水效率,其值分別為5.6%、4.77%、4.19%、44.9%。

    6 結(jié)論與討論

    (1)本次降水過程是受高空槽影響的積層混合云降水。降水分布不均勻,雨帶中存在多個強降水中心,其中位于北京地區(qū)的強降水為中雨量級。

    (2)積層混合云系是導致此次北京降水的主要云系。受高空槽影響,云系自西南向東北逐漸向北京推進;隨著積層混合云中對流云團移入北京并逐漸加強,北京地區(qū)的降水增強;隨后,積層混合云中對流云團東移出北京,北京區(qū)域降水減弱,云系出現(xiàn)減弱的趨勢。隨著北京北部云系消散,北京區(qū)域降水逐漸結(jié)束。積層混合云在衛(wèi)星 TBB分布圖上表現(xiàn)為大片層狀云中有多個絮狀對流云團,其云頂亮溫達-40°C。積層混合云系在水平和垂直方向都表現(xiàn)為不均勻分布。云系中的含水量水平分布不均勻,有多個高含水量中心,含水量大值區(qū)是層狀云區(qū)的3倍左右;垂直方向上看,有層云存在,云水含量較小,其中有多個積云單體存在,積層混合云中積云區(qū)對應(yīng)的降水較大。隨著積層混合云的東移加強,積云區(qū)加強且范圍增大,積云區(qū)的高空雪含量增大,對應(yīng)地面雨水增加。

    (3)北京地區(qū)上空積層混合云系的微物理結(jié)構(gòu)具有混合相云特征。高層為冰晶、雪,中層為霰,低層為云水和雨水;其中,雪的含水量最高。冷云降水過程(雪和霰的融化)是雨水形成的主要過程,暖云降水過程對雨水的形成也有貢獻;其中冷云降水過程占主導地位,雪的融化對雨水的形成貢獻最大。北京地區(qū)云系中不同云區(qū)的降水形成機制相同,均以冷云過程為降水形成的主要過程,冰相過程對降水形成起到重要的作用,但微物理過程對形成雨水的貢獻大小不同。

    (4)北京區(qū)域降水過程的主要水汽源地為黃海海面及蒙古國。一支水汽從黃海海面經(jīng)過江蘇、安徽、湖南,沿四川西部向東北輸送,另一支水汽由蒙古氣旋西側(cè)的西北氣流向東南輸送,兩支氣流在陜西北部匯合后的西南氣流將水汽輸送到華北地區(qū)。北京區(qū)域的水平水汽輸送以西、南入,東、北出為主。水汽強輻合區(qū)與強降水區(qū)相對應(yīng),降水發(fā)生區(qū)域的總體水汽收支和平均面降雨量的變化趨勢相吻合,隨著水汽收入的增加,降水量開始加大。

    (5)北京市區(qū)域以外,水汽和水凝物主要從西邊界和南邊界輸送到域內(nèi),具有西向風分量和南向分量的暖濕氣流是云系的水分輸送帶。雪、云水存在明顯的輸送,其中雪的輻合量最大,說明平流輸送從區(qū)域外凈流入的雪量對區(qū)域內(nèi)降水的增加有一定影響。在高層冰晶、雪、霰存在明顯的輻散,主要受高層氣流影響。北京區(qū)域降水主要時段內(nèi),水物質(zhì)通量在水平方向上為凈流入。

    (6)通過對北京區(qū)域水汽、水凝物和總水物質(zhì)的水分收支各項的估算表明,水物質(zhì)基本達到平衡。在此基礎(chǔ)上,估算了總水物質(zhì)降水效率、凝結(jié)(華)率和水凝物降水效率,北京市區(qū)域內(nèi)每小時的降水效率和凝結(jié)(華)率的演變分布與雨強的分布相似,且極值出現(xiàn)的時刻相近。北京市區(qū)域從 5日20時至6日14時,總體上總水物質(zhì)降水效率、凝結(jié)率、凝華率及水凝物降水效率分別為 5.6%、4.77%、4.19%、44.9%。

    降水效率有多種不同的定義。目前,國際上降水效率的定義主要分為三類。本文中的降水效率的定義基本屬于國際中云微物理降水效率(CMPE)這一類,它是總水物質(zhì)(總水凝物)的降水效率,全面考慮了水汽和各種水凝物(總水凝物)對降水的貢獻。本文中總水物質(zhì)(總水凝物)的降水效率計算,首先從中尺度數(shù)值模式預報水汽、云水、雨水、冰晶、雪和霰的方程得到的水汽、總水凝物和總水物質(zhì)的平衡等式出發(fā),在計算出水汽、水凝物和總水物質(zhì)收支基本達到平衡的前提下,根據(jù)平衡等式得出計算總水物質(zhì)(總水凝物)的降水效率公式。

    國際上主要針對熱帶海洋性對流云(如:臺風云系等)開展降水效率的評估研究,利用不同降水效率的定義,對比分析其雨水微物理降水效率(RMPE)、云微物理降水效率(CMPE)和大尺度降水效率(LSPE)的差異和原因(Weisman and Klemp,1982;Lipps and Hemler,1986;Ferrier et al.,1996;Li et al., 2002;Sui et al., 2005,2007;Gao et al., 2011;Shen et al.,2013)。國外多數(shù)研究中利用CMPE定義評估降水效率的對象主要是海洋性對流云,而本文評估降水效率的對象是大陸性積層混合云。由于不同性質(zhì)云其降水強度不同,導致不同類型云的降水效率存在較大差異。與海洋性對流云相比,大陸性積層混合云的降水強度較小,其降水效率也較小。因此,利用CMPE定義評估不同類型云的降水效率,其結(jié)果不具備可比性。

    另外,模式中不同云微物理參數(shù)化方案中云中水凝物(云水、冰晶、雪、霰等)的微物理過程不同,也能導致計算出的云微物理降水效率(CMPE)不同。本文的云微物理參數(shù)化方案采用中國氣象科學研究院CAMS復雜云微物理方案,此方案與中尺度模式中的顯示云微物理參數(shù)化方案不同,則云中水凝物(云水、冰晶、雪、霰等)的微物理過程不同,導致計算出的總水物質(zhì)(總水凝物)的降水效率不同。

    鑒于不同降水效率定義所得結(jié)果存在差異以及不同類型云的云微物理降水效率計算結(jié)果存在較大差異,本文中降水效率的結(jié)果與以往降水效率的研究結(jié)果不具有可比性。

    (References)

    Auer A H Jr, Marwitz J D.1968.Estimates of air and moisture flux into hailstorms on the high plains [J].J.Appl.Meteor., 7: 196–198.

    Braham R R Jr.1952.The water and energy budgets of the thunderstorm and their relation to thunderstorm development [J].J.Meteor., 9 (4):227–242.

    Chong M, Hauser D.1989.A tropical squall line observed during the CORT 81 experiment in West Africa.Part II: Water budget [J].Mon.Wea.Rev.,117: 728–744.

    Colle B A, Garvert M F, Wolfe J B, et al.2005.The 13–14 December 2001 IMPROVE-2 event.Part III: Simulated microphysical budgets and sensitivity studies [J].J.Atmos.Sci., 62 (10): 3535–3558.

    丁一匯, 胡國權(quán).2003.1998 年中國大洪水時期的水汽收支研究 [J].氣象學報, 61 (2): 129–144.Ding Yihui, Hu Guoquan.2003.A study on water vapor budget over China during the 1998 severe flood periods [J].Acta Meteor.Sinica (in Chinese), 61 (2): 129–144.

    Doswell C A III, Brooks H E, Maddox R A.1996.Flash flood forecasting:An ingredients-based methodology.Wea.Forecasting, 11 (4): 560–581.

    Ferrier B S, Simpson J, Tao W K.1996.Factors responsible for precipitation efficiencies in midlatitude and tropical squall simulations [J].Mon.Wea.Rev., 124 (10): 2100–2125.

    Gao S T, Li X F.2011.Can water vapour process data be used to estimate precipitation efficiency [J].Quart.J.Roy.Meteor.Soc., 137 (657): 969–978.

    Heymsfield G M, Schotz S.1985.Structure and evolution of a severe squall line over Oklahoma [J].Mon.Wea.Rev., 113 (9): 1563–1589.

    Hobbs P V, Matejka T J, Herzegh P H, et al.1980.The mesoscale and microscale structure and organization of clouds and precipitation in midlatitude cyclones.I: A case study of a cold front [J].J.Atmos.Sci., 37(3): 568–596.

    洪延超, 周非非.2006.層狀云系人工增雨潛力評估研究 [J].大氣科學,30 (5): 913–926.Hong Yanchao, Zhou Feifei.2006.The study of evaluation of potential of artificial precipitation enhancement in stratiform cloud system [J].Chinese Journal of Atmospheric Sciences (in Chinese), 30 (5): 913–926.

    胡國權(quán), 丁一匯.2003.1991 年江淮暴雨時期的能量和水汽循環(huán)研究[J].氣象學報, 61 (2): 146–163.Hu Guoquan, Ding Yihui.2003.A study on the energy and water cycles over Changjiang–Huaihe River basins during the 1991 heavy rain periods [J].Acta Meteor.Sinica (in Chinese), 61 (2): 146–163.

    胡志晉, 何觀芳.1987.積雨云微物理過程的數(shù)值模擬(一)微物理模式[J].氣象學報, 45 (4): 467–484.Hu zhijin, He Guanfang.1987.Numerical simulation of microprocesses in cumulonimbus clouds (I)microphysical model [J].Acta Meteor.Sinica (in Chinese), 45 (4):467–484.

    李宏宇, 王華, 洪延超.2006.鋒面云系降水中的增雨潛力數(shù)值研究 [J].大氣科學, 30 (2): 341–350.Li Hongyu, Wang Hua, Hong Yanchao.2006.A numerical study of precipitation enhancement potential in frontal cloud system [J].Chinese Journal of Atmospheric Sciences (in Chinese),30 (2): 341–350.

    Li X F, Sui C H, Lau K M.2002.Precipitation efficiency in the tropical deep convective regime: A 2-D cloud resolving modeling study [J].J.Meteor.Soc.Japan, 80 (2): 205–212.

    連志鸞, 李國翠.2005.石家莊的云、降水和水汽特征 [J].氣象科技, 33(S1): 21–26.Lian Zhiluan, Li Guocui.2005.Analysis of characteristics of clouds precipitation and water vapor over Shijiazhuang [J].Meteorological Science and Technology (in Chinese), 2005, 33(S1):21–26.

    Lipps F B, Hemler R S.1986.Numerical simulation of deep tropical convection associated with large-scale convergence [J].J.Atmos.Sci., 43(17): 1796–1816.

    劉衛(wèi)國, 劉奇駿.2007a.祁連山夏季地形云結(jié)構(gòu)和云微物理過程的模擬研究 (I): 模式云物理方案和地形云結(jié)構(gòu) [J].高原氣象, 26 (1): 1–15.Liu Weiguo, Liu Qijun.2007a.The numerical simulation of orographic cloud structure and cloud microphysical processes in Qilian Mountains in summer.Part (I): Cloud microphysical scheme and orographic cloud structure [J].Plateau Meteorology (in Chinese), 26 (1): 1–15.

    劉衛(wèi)國, 劉奇駿.2007b.祁連山夏季地形云結(jié)構(gòu)和云微物理過程的模擬研究 (II): 云微物理過程和地形影響 [J].高原氣象, 26 (1): 16–29.Liu Weiguo, Liu Qijun.2007b.The numerical simulation of orographic cloud structure and cloud microphysical processes in Qilian Mountains in summer.Part (II): Cloud microphysical processes and orographic influence [J].Plateau Meteorology (in Chinese), 26 (1): 16–29.

    樓小鳳.2002.MM5 模式的新顯式云物理方案的建立和耦合及原微物理方案的對比分析 [D].北京大學博士學位論文, 140pp.Lou Xiaofeng.2002.Development and implementation of a new explicit microphysical scheme and comparisons of original schemes of MM5 [D].Ph.D.dissertation (in Chinese), Peking University, 140pp.

    馬開玉, 高國棟.1992.長江流域典型旱/澇夏季大氣中的水汽輸送 [J].氣象科學, 12 (1): 48–56.Ma Kaiyu, Gao Guodong.1992.Analyses of water vapor transfer characteristics over China in typical drought/flood summer in the Yangtze River valley [J].Scientia Meteorologica Sinica (in Chinese), 12 (1): 48–56.

    平凡, 羅哲賢.2007.熱帶對流熱量與水汽收支的數(shù)值模擬研究 [J].地球物理學報, 50 (5): 1351–1361.Ping Fan, Luo Zhexian.2007.The numeral simulated study of convective heat and moisture budget in the tropical [J].Chinese Journal of Geophysics (in Chinese), 50 (5): 1351–1361.

    Shen X Y, Qin T, Li X F.2013.Effects of clouds, sea surface temperature,and its diurnal variation on precipitation efficiency [J].Chinese Physics B,22 (9): 1– 9.

    Sui C H, Li X, Yang M J, et al.2005.Estimation of oceanic precipitation efficiency in cloud models [J].J.Atmos.Sci., 62(12): 4358–4370.

    Sui C H, Li X, Yang M J.2007.On the definition of precipitation efficiency[J].J.Atmos.Sci., 64 (12): 4506–4513.

    孫晶, 樓小鳳, 胡志晉.2009a.祁連山冬季降雪個例模擬分析 (I): 降雪過程和地形影響 [J].高原氣象, 28 (3): 485–495.Sun Jing, Lou Xiaofeng, Hu Zhijin.2009a.Numerical simulation of snowfall in winter of Qilian Mountains.Part (I): Snowfall process and orographic influence[J].Plateau Meteorology (in Chinese), 28 (3): 485–495.

    孫晶, 樓小鳳, 史月琴, 等.2009b.祁連山冬季降雪個例模擬分析 (II):人工催化試驗 [J].高原氣象, 28 (3): 496–506.Sun Jing, Lou Xiaofeng, Shi Yueqing, et al.2009b.Numerical simulation of snowfall in winter of Qilian Mountains.Part (II): Seeding test [J].Plateau Meteorology (in Chinese), 28 (3): 496–506.

    Tao W K, Simpson J, Sui C H, et al.1983.Heating, moisture, and water budgets of tropical and midlatitude squall lines: Comparisons and sensitivity to longwave radiation [J].J.Atmos.Sci., 50 (5): 673–690.

    Tao W-K, Johnson D, Shie C-L, et al.2004.The atmospheric energy budget and large-scale precipitation efficiency of convective systems during TOGA COARE, GATE, SCSMEX, and ARM: Cloud-resolving model simulations [J].J.Atmos.Sci., 61 (20): 2405–2423.

    王昂生, 徐乃璋, 梁飛燕.1987.北京1980~1983年夏季云和降水的分析研究 [J].大氣科學, 11 (3): 297–303.Wang Angsheng, Xu Naizhang,Liang Feiyan.1987.The analysis and research of the cloud and the precipitation in the summers from 1980 to 1983 over the Beijing area [J].Chinese Journal of Atmospheric Sciences (in Chinese), 11 (3): 297–303.

    Weisman M L, Klemp J B.1982.The dependence of numerically simulated convective storms on vertical wind shear and buoyancy [J].Mon.Wea.Rev., 110 (6): 504–520.

    于翡, 姚展予.2009.一次積層混合云降水實例的數(shù)值模擬分析 [J].氣象, 35 (12): 3–11.Yu Fei, Yao Zhanyu.2009.Numerical study on the complex of the stratiform and embedded convective cloud precipitation:A case study [J].Meteorological Monthly (in Chinese), 35 (12): 3–11.

    張佃國, 姚展予, 龔佃利, 等.2011.環(huán)北京地區(qū)積層混合云微物理結(jié)構(gòu)飛機聯(lián)合探測研究 [J].大氣科學學報, 34 (1): 109–121.Zhang Dianguo, Yao Zhanyu, Gong Dianli, et al.2011.Microphysical structures of stratocumulus mixed cloud detected by aircraft around Beijing area [J].Transactions of Atmospheric Sciences (in Chinese), 34 (1): 109–121.

    周非非, 洪延超, 趙震.2010.一次層狀云系水分收支和降水機制的數(shù)值研究 [J].氣象學報, 68 (2): 182–194.Zhou Feifei, Hong Yanchao,Zhao Zhen.2010.A numerical study of the moisture budget and the mechanism for precipitation for a stratiform cloud system [J].Acta Meteor.Sinica (in Chinese), 68 (2): 182–194.

    周黎明, 張洪生, 王俊, 等.2010.一次典型積層混合云降水過程雨滴譜特征 [J].氣象科技, 38 (S1): 73–77.Zhou Liming, Zhang Hongsheng,Wang Jun, et al.2010.Raindrop spectral characteristics of mixed-cloud precipitation in Shandong Province [J].Meteorological Science and Technology (in Chinese), 38 (S1): 73–77.

    鄒倩, 劉奇駿, 劉衛(wèi)國.2008.積層混合云結(jié)構(gòu)和云微物理的數(shù)值模擬[J].干旱氣象, 26 (2): 74–83.Zou Qian, Liu Qijun, Liu Weiguo.2008.Numerical simulation of stratocumulus cloud structure and cloud microphysical processes [J].Arid Meteorology (in Chinese), 26 (2): 74–83.

    猜你喜歡
    云系水汽含水量
    青藏高原上空平流層水汽的時空演變特征
    2020年江西汛期大暴雨衛(wèi)星云圖特征分析
    2019年5月26日朝陽飛機人工增雨作業(yè)分析
    結(jié)合Sentinel- 1B和Landsat8數(shù)據(jù)的針葉林葉片含水量反演研究
    森林工程(2018年4期)2018-08-04 03:23:16
    一次冰雹過程成雹機理的數(shù)值模擬
    廣西11—12月人工增雨天氣研究
    1979~2011年間平流層溫度及平流層水汽的演變趨勢
    基于兩相混合流理論PEMFC含水量特性分析
    深圳“5·11”特大暴雨過程的水汽輸送特征分析
    液氯中含水量測定注意事項及改進
    中國氯堿(2014年10期)2014-02-28 01:05:00
    久久久久久久久久人人人人人人| 人妻一区二区av| 亚洲欧美精品自产自拍| 欧美日本中文国产一区发布| 一级av片app| 熟妇人妻不卡中文字幕| 麻豆成人午夜福利视频| 能在线免费看毛片的网站| 国产69精品久久久久777片| 国产欧美日韩精品一区二区| 亚洲欧美成人精品一区二区| 国产午夜精品久久久久久一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲成人一二三区av| 成人美女网站在线观看视频| 99国产精品免费福利视频| 日韩av不卡免费在线播放| 日韩中文字幕视频在线看片| 春色校园在线视频观看| 在线观看人妻少妇| 成年美女黄网站色视频大全免费 | 美女内射精品一级片tv| 国产精品蜜桃在线观看| 欧美人与善性xxx| 亚洲精品自拍成人| 亚洲国产欧美在线一区| 中文字幕av电影在线播放| 乱系列少妇在线播放| 曰老女人黄片| 午夜免费鲁丝| 成人无遮挡网站| 亚洲高清免费不卡视频| 夜夜骑夜夜射夜夜干| 亚洲国产色片| 韩国高清视频一区二区三区| 精品99又大又爽又粗少妇毛片| 久久精品国产鲁丝片午夜精品| av视频免费观看在线观看| 国产精品蜜桃在线观看| 我要看日韩黄色一级片| 欧美日韩视频精品一区| 777米奇影视久久| 青春草亚洲视频在线观看| 婷婷色综合www| 欧美3d第一页| 又黄又爽又刺激的免费视频.| 99re6热这里在线精品视频| 亚洲欧美精品自产自拍| 精品久久久久久久久亚洲| 丰满人妻一区二区三区视频av| 日韩精品有码人妻一区| 狂野欧美激情性xxxx在线观看| 欧美国产精品一级二级三级 | 亚洲国产欧美日韩在线播放 | 九色成人免费人妻av| 啦啦啦中文免费视频观看日本| 最近的中文字幕免费完整| 免费黄频网站在线观看国产| 欧美老熟妇乱子伦牲交| 国产成人精品无人区| 在线看a的网站| 亚洲av国产av综合av卡| 五月伊人婷婷丁香| 亚洲伊人久久精品综合| 少妇的逼好多水| 亚洲内射少妇av| 亚洲精品乱久久久久久| 国产综合精华液| 亚洲精品第二区| 91成人精品电影| 欧美激情国产日韩精品一区| 女性生殖器流出的白浆| 亚洲精品一二三| 亚洲久久久国产精品| 欧美高清成人免费视频www| 看免费成人av毛片| 91aial.com中文字幕在线观看| 久久久久久久久久成人| 在线播放无遮挡| 亚洲精品一二三| 色5月婷婷丁香| 久久久久久久久久人人人人人人| 久久久国产精品麻豆| av国产久精品久网站免费入址| 9色porny在线观看| 精品午夜福利在线看| 中国三级夫妇交换| 狂野欧美激情性xxxx在线观看| 日韩中字成人| 日本-黄色视频高清免费观看| 欧美日韩视频精品一区| 欧美日韩精品成人综合77777| 亚洲精品乱码久久久久久按摩| 成人毛片a级毛片在线播放| 国产免费一级a男人的天堂| 亚洲欧美精品自产自拍| 日本免费在线观看一区| av天堂中文字幕网| 亚洲精品久久久久久婷婷小说| 日本黄大片高清| 在线观看av片永久免费下载| 亚洲av欧美aⅴ国产| 男女啪啪激烈高潮av片| 91久久精品电影网| 成年美女黄网站色视频大全免费 | 亚洲精品乱码久久久久久按摩| 久久99一区二区三区| 国产有黄有色有爽视频| 久久人妻熟女aⅴ| 自拍偷自拍亚洲精品老妇| a级毛片在线看网站| 制服丝袜香蕉在线| 岛国毛片在线播放| 蜜桃在线观看..| 黄色日韩在线| 26uuu在线亚洲综合色| 国产精品国产三级国产av玫瑰| 18禁裸乳无遮挡动漫免费视频| 亚洲美女视频黄频| 国产欧美日韩综合在线一区二区 | a级毛片免费高清观看在线播放| 夫妻午夜视频| 亚洲精品国产av成人精品| 国产91av在线免费观看| 精品99又大又爽又粗少妇毛片| 亚洲国产成人一精品久久久| 美女视频免费永久观看网站| 亚洲精品国产av蜜桃| 99热这里只有精品一区| 亚洲人成网站在线播| 91精品一卡2卡3卡4卡| 99久久精品一区二区三区| 观看免费一级毛片| 丝袜喷水一区| 在线观看一区二区三区激情| 一个人免费看片子| 国产一区有黄有色的免费视频| 日韩视频在线欧美| 麻豆乱淫一区二区| 天美传媒精品一区二区| 校园人妻丝袜中文字幕| 少妇猛男粗大的猛烈进出视频| 99九九在线精品视频 | 久久国产精品大桥未久av | 久久午夜福利片| 精品久久久精品久久久| 亚洲精品国产av蜜桃| 秋霞伦理黄片| 午夜福利在线观看免费完整高清在| 91久久精品电影网| 国产一区二区三区av在线| 亚洲av成人精品一二三区| 国产亚洲91精品色在线| 91久久精品国产一区二区三区| 91在线精品国自产拍蜜月| 精品少妇内射三级| 99久久精品国产国产毛片| 熟女人妻精品中文字幕| 精品久久久噜噜| 老司机影院毛片| 日本色播在线视频| 亚洲欧洲精品一区二区精品久久久 | 国产成人免费观看mmmm| 少妇的逼好多水| 精品久久久久久久久av| 国产亚洲91精品色在线| 少妇裸体淫交视频免费看高清| 国产精品女同一区二区软件| 精品少妇内射三级| 熟女人妻精品中文字幕| 91精品国产九色| 老女人水多毛片| 黄色欧美视频在线观看| 亚洲欧美一区二区三区黑人 | 七月丁香在线播放| 精品一区在线观看国产| 国产淫片久久久久久久久| 成年av动漫网址| 我要看日韩黄色一级片| 我要看黄色一级片免费的| 中文字幕亚洲精品专区| 国产av一区二区精品久久| 大话2 男鬼变身卡| 中国三级夫妇交换| 国产精品嫩草影院av在线观看| 中文在线观看免费www的网站| 免费看光身美女| 精华霜和精华液先用哪个| 久久6这里有精品| 少妇 在线观看| 黄色欧美视频在线观看| 亚洲欧美一区二区三区黑人 | 精品一区二区免费观看| 精品午夜福利在线看| 噜噜噜噜噜久久久久久91| 在线观看美女被高潮喷水网站| 91久久精品国产一区二区成人| 国模一区二区三区四区视频| 99久国产av精品国产电影| 亚洲精品第二区| 中文字幕亚洲精品专区| 久久久a久久爽久久v久久| 免费看不卡的av| 免费看不卡的av| 免费大片黄手机在线观看| 国产黄片美女视频| 三级国产精品欧美在线观看| 最新的欧美精品一区二区| 国产欧美日韩综合在线一区二区 | 六月丁香七月| 一本大道久久a久久精品| 亚洲精品第二区| 欧美日本中文国产一区发布| 久久女婷五月综合色啪小说| 极品人妻少妇av视频| av一本久久久久| 春色校园在线视频观看| 午夜影院在线不卡| 纵有疾风起免费观看全集完整版| 一级毛片 在线播放| 热99国产精品久久久久久7| 成人综合一区亚洲| 欧美激情国产日韩精品一区| 久久人妻熟女aⅴ| 成人国产麻豆网| 在线观看三级黄色| 中国国产av一级| 一二三四中文在线观看免费高清| 男男h啪啪无遮挡| 国产精品久久久久成人av| av女优亚洲男人天堂| 久久99热6这里只有精品| 久久精品国产鲁丝片午夜精品| 在线观看一区二区三区激情| 97超碰精品成人国产| 亚洲在久久综合| 国产男女超爽视频在线观看| 国产成人精品久久久久久| 成人毛片60女人毛片免费| 国产精品国产三级专区第一集| 亚洲精品日韩av片在线观看| 99久久精品热视频| 黑人高潮一二区| 国产精品三级大全| 伦理电影免费视频| 人人妻人人添人人爽欧美一区卜| 久久6这里有精品| 国模一区二区三区四区视频| 免费黄频网站在线观看国产| 成人毛片a级毛片在线播放| 18禁在线播放成人免费| 蜜桃久久精品国产亚洲av| 人人妻人人澡人人爽人人夜夜| 国产av码专区亚洲av| 在线观看免费日韩欧美大片 | 国产 一区精品| 国产精品国产三级国产av玫瑰| 天堂8中文在线网| 男人爽女人下面视频在线观看| 欧美成人精品欧美一级黄| 国产一区有黄有色的免费视频| 国产在线视频一区二区| 91久久精品电影网| 亚洲av成人精品一二三区| 男女国产视频网站| 高清欧美精品videossex| 热re99久久精品国产66热6| 99热这里只有是精品50| 精品人妻熟女毛片av久久网站| 肉色欧美久久久久久久蜜桃| 免费人妻精品一区二区三区视频| 十分钟在线观看高清视频www | 女性被躁到高潮视频| 熟妇人妻不卡中文字幕| 午夜视频国产福利| 黄色日韩在线| 国产精品国产三级专区第一集| 99久久中文字幕三级久久日本| 成人18禁高潮啪啪吃奶动态图 | 国产男人的电影天堂91| 美女福利国产在线| 国产亚洲一区二区精品| 新久久久久国产一级毛片| 久热这里只有精品99| 在线天堂最新版资源| 成人亚洲欧美一区二区av| 最近中文字幕2019免费版| 26uuu在线亚洲综合色| 久久97久久精品| 在现免费观看毛片| 麻豆乱淫一区二区| 国产高清三级在线| 在线精品无人区一区二区三| 亚洲欧美精品自产自拍| 欧美bdsm另类| 精品久久国产蜜桃| 美女福利国产在线| 国产成人aa在线观看| 亚洲国产精品国产精品| 国产伦理片在线播放av一区| 成人亚洲欧美一区二区av| 一级,二级,三级黄色视频| 国产欧美日韩综合在线一区二区 | 色视频在线一区二区三区| 曰老女人黄片| 久久久久久久亚洲中文字幕| 亚洲成色77777| 成人国产av品久久久| 久久综合国产亚洲精品| 国产精品99久久99久久久不卡 | 色视频在线一区二区三区| 最近的中文字幕免费完整| 久久国产精品男人的天堂亚洲 | 久久精品久久久久久久性| 国产成人一区二区在线| 久久久久久久久久久久大奶| 亚洲精品乱久久久久久| 中文乱码字字幕精品一区二区三区| 日韩精品免费视频一区二区三区 | 秋霞在线观看毛片| 国产白丝娇喘喷水9色精品| 国产一区亚洲一区在线观看| 欧美成人精品欧美一级黄| 人妻一区二区av| 欧美日韩在线观看h| 亚洲色图综合在线观看| 国产亚洲av片在线观看秒播厂| 老司机影院毛片| 国产精品人妻久久久影院| 亚洲欧美日韩另类电影网站| 黄色日韩在线| 多毛熟女@视频| 性色avwww在线观看| 欧美成人精品欧美一级黄| kizo精华| 天天躁夜夜躁狠狠久久av| 久久久久视频综合| videossex国产| 成年av动漫网址| 免费看光身美女| 日本与韩国留学比较| 欧美国产精品一级二级三级 | 人人妻人人澡人人爽人人夜夜| 欧美日韩视频精品一区| 日韩强制内射视频| a 毛片基地| 成人综合一区亚洲| 99精国产麻豆久久婷婷| 亚洲伊人久久精品综合| 久久精品夜色国产| 久久久久国产网址| 亚洲精品乱久久久久久| 亚洲av综合色区一区| 午夜福利网站1000一区二区三区| 午夜福利视频精品| 亚洲人成网站在线观看播放| 亚洲精品色激情综合| 免费黄网站久久成人精品| 亚洲精品亚洲一区二区| 亚洲高清免费不卡视频| 美女中出高潮动态图| 人妻夜夜爽99麻豆av| 天堂8中文在线网| 免费黄网站久久成人精品| 99精国产麻豆久久婷婷| 黑人高潮一二区| 国精品久久久久久国模美| 婷婷色av中文字幕| 在线精品无人区一区二区三| 国产高清三级在线| 女人久久www免费人成看片| 亚洲欧洲国产日韩| 国产精品熟女久久久久浪| 99热这里只有是精品50| 成人毛片60女人毛片免费| 偷拍熟女少妇极品色| 美女主播在线视频| 国产成人免费无遮挡视频| 九九爱精品视频在线观看| 欧美精品人与动牲交sv欧美| 欧美日韩视频高清一区二区三区二| videos熟女内射| 波野结衣二区三区在线| 草草在线视频免费看| 久久午夜福利片| 国产日韩欧美亚洲二区| 国产av码专区亚洲av| 一本色道久久久久久精品综合| 九九久久精品国产亚洲av麻豆| 日韩亚洲欧美综合| 高清视频免费观看一区二区| 18禁裸乳无遮挡动漫免费视频| 国产无遮挡羞羞视频在线观看| 99re6热这里在线精品视频| 日本色播在线视频| 国产av国产精品国产| 乱码一卡2卡4卡精品| 国产精品一区二区性色av| 亚洲情色 制服丝袜| 久久精品国产亚洲av天美| 国产午夜精品久久久久久一区二区三区| 亚洲国产精品一区三区| 国产无遮挡羞羞视频在线观看| 亚洲在久久综合| 26uuu在线亚洲综合色| 久久韩国三级中文字幕| 精品亚洲乱码少妇综合久久| 99视频精品全部免费 在线| 多毛熟女@视频| 国产视频内射| 久久久久国产网址| 日韩制服骚丝袜av| 一本大道久久a久久精品| 在线观看免费视频网站a站| 91午夜精品亚洲一区二区三区| 亚洲精品国产色婷婷电影| 妹子高潮喷水视频| 成年人免费黄色播放视频 | 99热网站在线观看| 国产老妇伦熟女老妇高清| 80岁老熟妇乱子伦牲交| 亚洲真实伦在线观看| 91在线精品国自产拍蜜月| 久久精品国产a三级三级三级| 精品国产一区二区三区久久久樱花| av又黄又爽大尺度在线免费看| 亚洲欧美一区二区三区国产| 午夜福利,免费看| av福利片在线| 成人无遮挡网站| 国产欧美日韩精品一区二区| 自拍偷自拍亚洲精品老妇| 黄色毛片三级朝国网站 | 国产伦精品一区二区三区视频9| 一区二区三区免费毛片| 国产免费又黄又爽又色| 亚洲一区二区三区欧美精品| 嫩草影院新地址| 伊人久久国产一区二区| 两个人的视频大全免费| 亚洲av电影在线观看一区二区三区| 十分钟在线观看高清视频www | 黑人巨大精品欧美一区二区蜜桃 | 国语对白做爰xxxⅹ性视频网站| 欧美一级a爱片免费观看看| 亚洲av中文av极速乱| 久久99热这里只频精品6学生| 亚洲av综合色区一区| 看十八女毛片水多多多| 国产精品人妻久久久影院| 少妇熟女欧美另类| 日本黄色日本黄色录像| 99热网站在线观看| 一本久久精品| 中文字幕亚洲精品专区| 全区人妻精品视频| 又大又黄又爽视频免费| 亚洲,一卡二卡三卡| 亚洲av成人精品一区久久| 国产精品国产三级专区第一集| 久久国产亚洲av麻豆专区| 男女无遮挡免费网站观看| 国产男人的电影天堂91| 黄色欧美视频在线观看| 久久久久视频综合| 精品人妻一区二区三区麻豆| 3wmmmm亚洲av在线观看| 国产中年淑女户外野战色| 久久久久久久久久久久大奶| 校园人妻丝袜中文字幕| 日本91视频免费播放| 亚洲av国产av综合av卡| 国产精品久久久久久av不卡| 日韩一区二区视频免费看| 成人黄色视频免费在线看| 欧美 亚洲 国产 日韩一| 在线播放无遮挡| 精品一区二区免费观看| 久久精品夜色国产| 国产精品一区二区在线观看99| 国产精品三级大全| 蜜桃久久精品国产亚洲av| 22中文网久久字幕| 亚洲成人手机| 赤兔流量卡办理| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品成人久久小说| 国产在线免费精品| 日韩av免费高清视频| 精品99又大又爽又粗少妇毛片| 天天操日日干夜夜撸| 国产免费视频播放在线视频| 成人综合一区亚洲| av福利片在线| 国产黄色视频一区二区在线观看| 性高湖久久久久久久久免费观看| 丰满人妻一区二区三区视频av| 最近的中文字幕免费完整| 亚洲中文av在线| av一本久久久久| videos熟女内射| 如日韩欧美国产精品一区二区三区 | 欧美高清成人免费视频www| 桃花免费在线播放| 国产精品伦人一区二区| 边亲边吃奶的免费视频| av卡一久久| 国产成人午夜福利电影在线观看| 亚洲av男天堂| 亚洲精品456在线播放app| 纵有疾风起免费观看全集完整版| 少妇猛男粗大的猛烈进出视频| 人人妻人人看人人澡| 爱豆传媒免费全集在线观看| 欧美 日韩 精品 国产| 欧美精品一区二区大全| 卡戴珊不雅视频在线播放| 国产在线男女| 丰满人妻一区二区三区视频av| 国产精品久久久久久久久免| 22中文网久久字幕| 一级片'在线观看视频| 国产欧美日韩一区二区三区在线 | 欧美激情国产日韩精品一区| 观看av在线不卡| 亚洲欧美日韩另类电影网站| 久久毛片免费看一区二区三区| 男女啪啪激烈高潮av片| 夜夜骑夜夜射夜夜干| av天堂中文字幕网| 美女内射精品一级片tv| 看免费成人av毛片| 噜噜噜噜噜久久久久久91| 精品久久久噜噜| 99久久精品热视频| 国产精品久久久久久久久免| 亚洲精品视频女| 老司机亚洲免费影院| 精品酒店卫生间| 黑人猛操日本美女一级片| 国产日韩欧美亚洲二区| 国产一区二区在线观看日韩| 少妇人妻精品综合一区二区| 欧美xxⅹ黑人| 能在线免费看毛片的网站| www.av在线官网国产| 国产高清国产精品国产三级| 青春草视频在线免费观看| 亚洲综合色惰| 熟妇人妻不卡中文字幕| 成人黄色视频免费在线看| 免费在线观看成人毛片| 男的添女的下面高潮视频| 欧美 亚洲 国产 日韩一| 波野结衣二区三区在线| av国产精品久久久久影院| 亚洲国产精品999| 男女国产视频网站| 免费人成在线观看视频色| 亚州av有码| 91成人精品电影| 色婷婷av一区二区三区视频| 99热这里只有精品一区| 观看av在线不卡| 亚洲成人一二三区av| 国产精品国产三级专区第一集| 精品人妻熟女av久视频| 国产精品一区二区性色av| 一本久久精品| 精品国产一区二区三区久久久樱花| 国产免费视频播放在线视频| 国产在线视频一区二区| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华液的使用体验| 精品久久国产蜜桃| 亚洲丝袜综合中文字幕| 国产男人的电影天堂91| 国产精品麻豆人妻色哟哟久久| 七月丁香在线播放| 乱码一卡2卡4卡精品| 少妇人妻久久综合中文| 国产91av在线免费观看| 国产一区二区在线观看日韩| 香蕉精品网在线| 五月伊人婷婷丁香| 晚上一个人看的免费电影| av免费观看日本| 晚上一个人看的免费电影| 观看av在线不卡| 国产欧美日韩一区二区三区在线 | 亚洲真实伦在线观看| 日韩一区二区三区影片| 欧美亚洲 丝袜 人妻 在线| 一级黄片播放器| 日本午夜av视频| 亚洲精品乱码久久久v下载方式| 久久久久久伊人网av| 亚洲美女搞黄在线观看| 精品国产一区二区久久| 老熟女久久久| 国产成人免费观看mmmm| 丝袜喷水一区| 国产亚洲91精品色在线| 男女无遮挡免费网站观看| 精品久久国产蜜桃| 少妇的逼水好多| 青春草亚洲视频在线观看| 国产毛片在线视频| 日韩成人av中文字幕在线观看| 99re6热这里在线精品视频| 日本爱情动作片www.在线观看| 亚洲在久久综合| 亚洲欧美日韩卡通动漫| 成人免费观看视频高清| 国产成人a∨麻豆精品|