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

1中國(guó)氣象科學(xué)研究院中國(guó)氣象局人工影響天氣中心，北京100081

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

1 引言

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

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

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

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

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

2 降水過(guò)程實(shí)況分析

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

2.1 天氣背景分析

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

2.2 云系的演變特征

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

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

圖1 2007年10月5日20時(shí)～6日20時(shí)24 h 地面降水量（單位：mm）：（a）實(shí)況；（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時(shí)的（a）500 hPa高空天氣形勢(shì)圖和（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時(shí)和（b）6日05時(shí)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 積層混合云降水過(guò)程的數(shù)值模擬

3.1 模式及模擬方案介紹

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

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

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

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

圖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 模擬（上）和實(shí)測(cè)（下）的6小時(shí)地面降水量（單位：mm）：（a、d）2007年10月5日20時(shí)～6日02時(shí)；（b、e）2007年10月6日02時(shí)～08時(shí)（圖5b中矩形代表北京區(qū)域）；（c、f）2007年10月6日08時(shí)～14時(shí)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

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

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

4.1 云系的水平分布

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

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

由此可見(jiàn)，積層混合云系是導(dǎo)致此次北京降水的主要云系。隨著積層混合云中對(duì)流云團(tuán)移入北京并逐漸加強(qiáng)，北京地區(qū)的降水增強(qiáng)；隨后，積層混合云中對(duì)流云團(tuán)東移出北京，云系出現(xiàn)減弱的趨勢(shì)，北京區(qū)域降水減弱。

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

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

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

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

之后，強(qiáng)降水云東移并逐漸變?nèi)酰绷Ｗ雍繙p小，大值區(qū)東移出北京，冷云降水作用減弱，該地區(qū)降水也隨之減弱。

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

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

雨水的形成貢獻(xiàn)也較小。相對(duì)于西部和東部，北京中部云區(qū)中霰的融化過(guò)程對(duì)雨水的形成貢獻(xiàn)較大，冷云降水過(guò)程更強(qiáng)，其地面雨水含量最大。

圖6 2007年10月6日三個(gè)時(shí)刻模擬的云系總水凝物含水量厚度（單位：mm）的分布：（a）01時(shí)；（b）05時(shí)；（c）15時(shí)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時(shí)，通過(guò)北京站（40°N）含水量（單位：g kg?1）的經(jīng)向—垂直剖面：（a）云水（彩色陰影）、冰晶（紅色實(shí)線）含水量；（b）雪（彩色陰影）、霰（紅色實(shí)線）和雨水（綠色實(shí)線）含水量。黑線：等溫線（°C）；箭頭：風(fēng)場(chǎng)（單位：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日三個(gè)時(shí)刻水凝物含水量區(qū)域累加量（單位：g kg?1）的垂直廓線：（a）01 時(shí)；（b）05 時(shí)；（c）11 時(shí)。（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時(shí)北京區(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)

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

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

5.1 水汽的源地及輸送通道

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

圖10 2007年10月6日三個(gè)時(shí)刻模擬的700 hPa上水汽通量散度（彩色陰影，單位：g s?1 cm?1 hPa?1；負(fù)值表示水汽的輻合）和風(fēng)場(chǎng) (單位：m s?1) 的水平分布：（a）01 時(shí)；（b）05 時(shí)；（c）15 時(shí)。圖中矩形區(qū)域?yàn)楸本┦蟹秶鶩ig.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è)，其中水平范圍較小的一個(gè)水汽輻合中心的前沿開(kāi)始移入北京西南部，此時(shí)北京區(qū)域并未出現(xiàn)明顯降水。隨著天氣系統(tǒng)的發(fā)展移動(dòng)，水汽輻合大值區(qū)隨之向東北伸展加強(qiáng)，在降水強(qiáng)盛期（圖10b），水汽輻合高值區(qū)基本上覆蓋了北京大部分地區(qū)，其輻合中心的最大值由降水初期的 5 g s?1cm?1hPa?1增加到 10 g s?1cm?1hPa?1，可見(jiàn)充足的水汽輸送，為北京地區(qū)的強(qiáng)降水創(chuàng)造了條件。之后，隨著系統(tǒng)的進(jìn)一步東移，西南向水汽輸送逐步減弱后消失，水汽通量的大值區(qū)也隨之移出北京區(qū)域，北京區(qū)域降水逐漸結(jié)束（圖10c）。可見(jiàn)，水汽強(qiáng)輻合區(qū)與強(qiáng)降水區(qū)相對(duì)應(yīng)，隨著北京地區(qū)水汽輻合量的增加，北京地區(qū)降水量開(kāi)始加大。

圖11 2007年10月6日02時(shí)700 hPa上實(shí)測(cè)的水平流場(chǎng)Fig.11 The observed flow field of wind at 700 hPa at 0200 BJT 6 October 2007

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

為了了解各邊界不同高度上水物質(zhì)輸送特點(diǎn)，對(duì)北京區(qū)域（115.3°N～117.5°N，39.3°E～41.5°E）各邊界不同高度上的的水物質(zhì)通量和輻合量。輻合量＝（東邊界水物質(zhì)通量－西邊界水物質(zhì)通量）＋（北邊界水物質(zhì)通量－南邊界水物質(zhì)通量），若大于 0，則水物質(zhì)通量是輻散（凈流出），若小于0，則水物質(zhì)通量是輻合（凈流入）。

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

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

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

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

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

圖13 穿過(guò)北京區(qū)域各邊界截面積的（a）水汽通量（qv）和輻合量的垂直積分和（b）總水凝物（云水、雨水、霰、雪和冰晶）通量和輻合量的垂直積分隨時(shí)間變化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 水分平衡等式中各項(xiàng)值（單位：kg）Table 1 The computational values of each term in the equilibrium equations of water substances (units: kg)

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

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

圖14 北京區(qū)域每小時(shí)總水物質(zhì)降水效率、凝結(jié)率、凝華率、水凝物降水效率和地面平均小時(shí)降水量（單位：mm h?1）的時(shí)間分布（注：每小時(shí)的總水物質(zhì)降水效率、凝結(jié)率和凝華率、水凝物降水效率對(duì)應(yīng)左邊y軸，地面平均小時(shí)降水量對(duì)應(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

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

6 結(jié)論與討論

（1）本次降水過(guò)程是受高空槽影響的積層混合云降水。降水分布不均勻，雨帶中存在多個(gè)強(qiáng)降水中心，其中位于北京地區(qū)的強(qiáng)降水為中雨量級(jí)。

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

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

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

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

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

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

國(guó)際上主要針對(duì)熱帶海洋性對(duì)流云（如：臺(tái)風(fēng)云系等）開(kāi)展降水效率的評(píng)估研究，利用不同降水效率的定義，對(duì)比分析其雨水微物理降水效率（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）。國(guó)外多數(shù)研究中利用CMPE定義評(píng)估降水效率的對(duì)象主要是海洋性對(duì)流云，而本文評(píng)估降水效率的對(duì)象是大陸性積層混合云。由于不同性質(zhì)云其降水強(qiáng)度不同，導(dǎo)致不同類型云的降水效率存在較大差異。與海洋性對(duì)流云相比，大陸性積層混合云的降水強(qiáng)度較小，其降水效率也較小。因此，利用CMPE定義評(píng)估不同類型云的降水效率，其結(jié)果不具備可比性。

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

鑒于不同降水效率定義所得結(jié)果存在差異以及不同類型云的云微物理降水效率計(jì)算結(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.

丁一匯, 胡國(guó)權(quán).2003.1998 年中國(guó)大洪水時(shí)期的水汽收支研究 [J].氣象學(xué)報(bào), 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.層狀云系人工增雨潛力評(píng)估研究 [J].大氣科學(xué),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.

胡國(guó)權(quán), 丁一匯.2003.1991 年江淮暴雨時(shí)期的能量和水汽循環(huán)研究[J].氣象學(xué)報(bào), 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.積雨云微物理過(guò)程的數(shù)值模擬(一)微物理模式[J].氣象學(xué)報(bào), 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].大氣科學(xué), 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.

連志鸞, 李國(guó)翠.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)國(guó), 劉奇駿.2007a.祁連山夏季地形云結(jié)構(gòu)和云微物理過(guò)程的模擬研究 (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)國(guó), 劉奇駿.2007b.祁連山夏季地形云結(jié)構(gòu)和云微物理過(guò)程的模擬研究 (II): 云微物理過(guò)程和地形影響 [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 模式的新顯式云物理方案的建立和耦合及原微物理方案的對(duì)比分析 [D].北京大學(xué)博士學(xué)位論文, 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.

馬開(kāi)玉, 高國(guó)棟.1992.長(zhǎng)江流域典型旱/澇夏季大氣中的水汽輸送 [J].氣象科學(xué), 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.熱帶對(duì)流熱量與水汽收支的數(shù)值模擬研究 [J].地球物理學(xué)報(bào), 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.祁連山冬季降雪個(gè)例模擬分析 (I): 降雪過(guò)程和地形影響 [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.祁連山冬季降雪個(gè)例模擬分析 (II):人工催化試驗(yàn) [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].大氣科學(xué), 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í)例的數(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.

張佃國(guó), 姚展予, 龔佃利, 等.2011.環(huán)北京地區(qū)積層混合云微物理結(jié)構(gòu)飛機(jī)聯(lián)合探測(cè)研究 [J].大氣科學(xué)學(xué)報(bào), 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.一次層狀云系水分收支和降水機(jī)制的數(shù)值研究 [J].氣象學(xué)報(bào), 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.一次典型積層混合云降水過(guò)程雨滴譜特征 [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)國(guó).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.