李義宇 郭學良 金蓮姬 李培仁 孫鴻娉 趙德龍 馬新成
1 山西省人工增雨防雷技術中心, 太原 030002
2 中國科學院大氣物理研究所, 北京 100029
3 南京信息工程大學中國氣象局氣溶膠與云降水重點開放實驗室, 南京 210044
4 北京市人工影響天氣辦公室, 北京 100089
氣溶膠是云降水形成的重要組成部分。由于對氣溶膠時空分布、化學組分及其核化特性等方面認識仍然不足,導致氣溶膠的間接氣候效應在輻射強迫中具有很大的不確定性(IPCC, 2007)。氣溶膠的間接氣候效應主要是通過云凝結核(CCN)和冰核(IN)影響云覆蓋、生命史和降水形成和演變過程,從而改變大氣輻射平衡。因此,深入了解氣溶膠的時空分布特征、來源及其與CCN、云滴數濃度的關系,對于認識氣溶膠核化成云及其間接氣候效應具有重要的意義。
在20 世紀50年代,國外科學家就開始關注氣溶膠對云宏、微觀特性的影響(Gunn and Phillips,1957; Squires, 1958; Squires and Twomey, 1960;Warner and Twomey, 1967; Warner, 1968),研究指出人類活動所產生的高濃度氣溶膠、CCN 可以增加云滴數濃度,潛在地降低了降水效率。Terliuc and Gagin(1971)在地面及飛機上進行的CCN觀測研究表明,雨天、干燥天CCN 譜呈典型的大陸性分布,CCN 濃度變化受風向影響顯著,CCN 對降水機制起重要作用。Hudson(1993)的研究表明,人為CCN 會影響層云,最明顯的證據是1987年衛星觀測到了船煙跡云。Rosenfeld(1999)利用TRMM 衛星觀測資料,分析了森林大火后,煙塵對熱帶對流云中暖云降水過程的影響,發現煙塵影響到云的形成,降水受到抑制。Rosenfeld(2000)認為,隨著污染的加重,城市居住區降水將有所減少。Chuang et al.(2000)、Snider and Brenguier(2000)、Ramanathan et al.(2001a, 2001b)通過分析機載CCN、云粒子探測資料發現,CCN 的增加會導致云滴濃度的增加,但CCN 濃度的進一步增加,云滴濃度的增加將趨于停止。
通過飛機觀測可以最直接有效地開展氣溶膠與云物理特性的研究。自20 世紀90年代以來,我國學者相繼在華北、華東和華中等地區開展了一系列氣溶膠與云凝結核的飛機觀測試驗。在氣溶膠特征的飛機觀測方面,Zhang et al.(2006)、范燁等(2007)、陳鵬飛等(2012)、馬新成等(2016)、段婧等(2019)、楊怡曼等(2020)、董曉波等(2021)在我國不同地區開展研究,包括不同天氣背景下氣溶膠分布特征對比,晴天和重污染天氣下氣溶膠垂直分布,氣溶膠垂直分布與大氣邊界層的關系,不同模態氣溶膠空間分布特征,以及海洋與大陸氣溶膠垂直分布特征的對比。范燁等(2007)分析了2004年的8~9月北京地區8架次氣溶膠飛行探測資料,發現風速較小的晴天和陰天都有可能使氣溶膠濃度出現高值,大氣有逆溫層存在時,逆溫層底有氣溶膠粒子的明顯累積,氣溶膠粒子濃度譜為單峰分布。孫霞等(2011)對比分析石家莊地區春季晴、霾天氣溶膠數濃度的垂直分布,發現霾天粗粒子氣溶膠的數濃度比晴天高一個量級。張瑜等(2011)綜合分析2005~2007年河北秋季不同天氣下的氣溶膠垂直分布特征,發現氣溶膠在降水天氣呈雙峰分布,其余天氣以單峰分布為主,陸地上空的氣溶膠比海洋上空高一個量級。在云微物理特征的飛機觀測方面,Deng et al.(2009)、李軍霞等(2014)、龐朝云等(2016)、郝囝等(2019)、蔡兆鑫等(2019)針對不同云系、不同區域和不同過程中的云微物理參數的時空分布開展了相關研究。石立新和段英(2007)發現華北地區的CCN 主要來源于地面,近地層CCN 濃度較高,CCN 濃度隨高度增加而減少;污染地區鄉村上空CCN 濃度比無污染地區鄉村高5 倍以上。盧廣獻和郭學良(2012)研究了2009年春季北京周邊大氣氣溶膠的分布、來源特征及其與云凝結核(CCN)的轉化關系,結果表明,高濃度氣溶膠基本分布在4500 m以下的區域,氣溶膠平均粒子譜呈雙(多)峰分布,而 在4500 m 以 上 基 本 為 單 峰 分 布。Zhang et al.(2011)分析了北京2008年7月的氣溶膠觀測資料,發現云滴數濃度(Nc)和氣溶膠數濃度(Na)的比(Nc/Na)與氣溶膠尺度和化學組分有關,其中Nc/Na與粒徑呈線性關系;尺度大的氣溶膠粒子能促進CCN、云滴的形成。Li et al.(2015)分析了山西地區2013年的氣溶膠觀測資料,發現過飽和度越大,CCN 核化比率越大。然而,由于飛行費用、飛行安全、空域限制等原因,利用飛機開展氣溶膠和云物理空間分布的觀測研究仍然較為缺乏。
云下氣溶膠與CCN 及云中云滴數濃度的關系可以揭示氣溶膠成云的特征,利用飛機對淡積云下氣溶膠、CCN 與云底之上云滴進行原位測量,在國內鮮有,氣溶膠、CCN、云滴三者的擬合關系可以為改進模式提供重要的依據。O’ Dowd et al.(1999)分析了北大西洋和東太平洋云下Na和Nc之間的關系,得到了云滴數和氣溶膠濃度的關系式。本文利用2013~2014年期間在華北中部山西地區開展的9架次飛機觀測資料,分析了氣溶膠垂直分布、譜分布和來源特征,以及與CCN、云中云滴數濃度的轉化關系,為揭示氣溶膠在成云過程中的作用提供科學依據。
觀測資料來源于2013~2014年期間山西省人工降雨防雹辦公室運-12 飛機觀測數據,該飛機觀測 平 臺 搭 載 美 國DMT( Droplet Measurement Technologies)公司的機載云降水粒子測量系統。本文用到的儀器有氣溶膠譜儀PCAPS-100X(Passive Cavity Aerosol Spectrometer Probe),量程為0.1~3 μm,分為15 個不同間隔的通道。雖然PACSPS-100X 探頭測量范圍有限,只能代表積聚模態(0.1~1 μm)和部分粗粒子模態(1~3 μm)的氣溶膠粒子。但PCAPS-100X 探頭是專用的機載設備,其觀測數據相對可靠。蔡兆鑫等(2021),Yang et al.(2020)都用PCAPS-100X 探頭來研究氣溶膠的垂直分布。本文用到0.1~3 μm 的氣溶膠作為CCN,計算所得的CCN 比率較實際值偏大;云 凝 結 核 計 數 器CCN-200(Cloud Condensation Nuclei Counter-200);云 粒 子 譜 儀CDP(Cloud Droplet Probe),量程為3~50 μm,分為30 個不同間隔的通道。儀器采樣頻率都為1 Hz,每年按時返回美國DMT 公司進行標校,并每次觀測前對儀器進行維護,以確保觀測數據質量。
飛行方案是起飛后進行垂直探測獲得氣溶膠和CCN 濃度的垂直分布數據,晴空時在目標區域以300 m 或600 m 為間隔分層盤旋平飛探測,獲得不同高度層的氣溶膠資料。在淡積云天氣時,飛機在云底之下100~300 m 范圍平飛探測,獲得云下氣溶膠和CCN 的分布數據,然后在云底之上100~200 m 水平穿云探測,獲得云中云滴粒子數據。2014年8月11日起飛過程CCN 工作狀態異常,因而沒做這次CCN 數據分析。飛行區域如圖1所示,主要在山西省中部的忻州、文水地區。

圖1 2013~2014年間在華北飛行區域示意圖(矩形區域為飛行區域)Fig. 1 Experimental area for aircraft measurements in northern China from 2013 to 2014. The rectangular area denotes the aircraft flight region
本文依據大氣邊界層逆溫情況,將2013~2014年夏季氣溶膠垂直分布分為兩類,一類是大氣邊界層有明顯的逆溫情況,大氣邊界層頂高度一般在1500 m 左右。另一類無明顯逆溫情況。
3.1.1 不同大氣層結條件下的氣溶膠垂直分布特征
圖2 和圖3 分別給出了邊界層有無逆溫情況的Na與氣溶膠有效直徑(ED)隨高度的變化,同時左邊給出了對應的溫濕度廓線。從圖2 可以看到,在有邊界層逆溫的情況下,由于逆溫層的阻擋作用,使氣溶膠主要集中在邊界層頂(1500 m)以下,氣溶膠在邊界層內出現明顯的累積現象,Na最大值達到3000~4000 cm?3以上。而邊界層頂以上氣溶膠濃度出現劇減現象,可減小2 個量級左右。可見,因逆溫造成的穩定性邊界層對氣溶膠的累積效應,在低能見度霧—霾天氣形成過程中起到至關重要的作用。由于邊界層內強烈的湍流混合作用,使氣溶膠混合均勻,ED 均在0.4 μm 左右。而無邊界層逆溫存在的情況下(圖3),氣溶膠垂直擴散混合作用加強,無明顯累積效應。不同大氣層結條件下的垂直擴散程度不同,圖3 中2013年8月30日顯示的情況表示大氣溫度隨高度的遞減率大,大氣層結不穩定性高,垂直擴散作用強,從地面到3000 m高度的氣溶膠濃度均在500 cm?3,但ED 出現隨高度遞減現象,這是因為上升擴散氣流對不同大小的氣溶膠粒子具有篩選作用,小粒子易被上升氣流帶到高層,而大粒子留在低層。可見此種情況下的對流擴散比較強烈,有利于地面污染氣溶膠向上輸送,對氣溶膠參與云和降水形成過程提供了有利條件。而其他兩個例子2014年8月11日和2014年8月12日的大氣溫度隨高度的遞減率較小,不穩定性較弱,垂直擴散作用偏小,氣溶膠濃度和大小均出現明顯的垂直遞減現象,進入空中的氣溶膠量明顯減小。

圖2 2013年8月3日(綠色)、4日(紅色)和12日(藍色)邊界層逆溫條件下氣溶膠(a)數濃度Na(單位:cm?3)與(b)有效直徑ED(單位:μm)的垂直分布及(c)溫度、濕度廓線(實線為溫度,虛線為相對濕度)Fig. 2 Vertical distributions of (a) number concentration (Na, units: cm?3), (b) effective diameter (ED, units: μm), and (c) profiles of temperature(solid lines) and relative humidity (dashed lines) with atmospheric inversion condition on 3 (green), 4 (red), and 12 (blue) August 2013

圖3 2013年8月30日(青色)、2014年8月11日(藍色)和12日(粉色)邊界層無逆溫條件下(a)Na 與(b)ED 的垂直分布及(c)溫度、濕度廓線Fig. 3 Vertical distributions of (a) number concentration (Na, units: cm?3), (b) effective diameter (ED, units: μm), and (c) profiles of temperature(solid lines) and relative humidity (dashed lines) without atmospheric inversion condition on 30 August 2013 (cyan), 11 August 2014 (blue), and 12 August 2014 (pink)
表1 給出了6 個例子在不同高度的Na統計特征,可以看到,逆溫情況下1500 m 以下Na的平均值為3189 cm?3,而無逆溫情況下邊界層內的Na平均值只有1509 cm?3。逆溫情況下ED 在1500 m 以下平均值為0.36 μm;在1500 m 以上ED 隨高度的增加呈現增大的趨勢,1500~3000 m ED 平均值為0.42 μm,3000 m 以上,ED 平均值為0.53 μm。在逆溫情況下,出現ED 隨高度增大情況應與氣溶膠來源不同有關,在后面詳細討論。沒有邊界層逆溫的情況下,1500 m 以下,ED 平均值為0.45 μm;1500~3000 m ED 平均值為0.35 μm;3000 m 以上,ED 平均值為0.32 μm。ED 隨高度變化很小或出現減小現象,這與垂直擴散氣流對粒子輸送過程的篩選作用有關。

表1 不同高度氣溶膠數濃度Na 統計特征Table 1 Statistical summary of number concentration (Na) at different heights in each flight
3.1.2 不同大氣層結條件下的CCN 濃度垂直分布特征
圖4 是不同大氣層結條件下CCN 數濃度(過飽和度為0.3%)的垂直分布。從圖中可以看出,CCN 數濃度垂直分布與氣溶膠數濃度垂直分布相似,也與大氣穩定度密切相關。在邊界層有明顯逆溫的穩定大氣層結條件下(圖4a),CCN 數濃度在邊界層內也存在明顯累積現象,然后隨高度出現劇減,CCN 數濃度最大值可以達到2000~3000 cm?3,而對應的氣溶膠數濃度最大值為3000~4000 cm?3,說明在過飽和度為0.3%的條件下,邊界層內累積的高濃度氣溶膠具有很高的核化率。在無明顯逆溫的不穩定大氣層結條件下(圖4b),CCN 數濃度在邊界層內無明顯累積現象,而是隨高度呈遞減趨勢。
表2 給出了5 個例子的CCN 數濃度的平均值和最大值。逆溫層結條件下,1500 m 以下CCN 數濃度平均值為1992 cm?3,最大值為2656 cm?3;1500~3000 m 高度,CCN 數濃度平均值為610 cm?3,最大值為1119 cm–3;3000 m 以上,CCN 數濃度平均值僅為475 cm?3,最大值為562 cm?3。無明顯逆溫層結條件下,1500 m 以下的CCN 數濃度平均值為1190 cm?3,最大值為1826 cm?3;1500~3000 m高度層,CCN 數濃度平均值為701 cm?3,最大值為885 cm?3;在3000 m 以上,CCN 數濃度平均值為360 cm?3,平均值為586 cm?3。對比兩種類型的CCN 數濃度垂直分布,可以看出,在邊界層內,逆溫情況下的CCN 數濃度約為無逆溫層結情況下的1.5 倍,而在1500~3000 m、3000 m 以上,兩種類型的CCN 數濃度較為接近,這與兩類氣溶膠數濃度的垂直分布趨勢基本一致。
3.1.3 氣溶膠來源分析
為解釋氣溶膠濃度或尺度隨高度增加,但CCN 濃度不一定增加的現象(如圖2a 和圖4a),了解不同高度氣溶膠的來源具有重要意義。利用HYSPLIT_4模式(HYbrid Single-Particle Lagrangian Integrated Trajectory model)分別對6次個例在不同高度氣溶膠粒子做了后向軌跡分析。圖5 為6次個例在距地面500 m、1200 m、2200 m、2700 m高度的氣團72 小時的后向軌跡。

圖4 (a)有邊界層逆溫和(b)無邊界層逆溫情況下CCN 數濃度(過飽和度為0.3%)垂直分布Fig. 4 Averaged vertical distributions of the CCN number concentration at 0.3% supersaturation (a) with atmospheric inversion condition and (b)without atmospheric inversion condition

圖5 2013年8月(a)3日、(b)4日、(c)12日、(d)30日以及2014年8月(e)11日、(f)12日觀測期間太原不同高度72 小時后向軌跡 紅、藍、綠、青色實線分別代表500 m、1200 m、2200 m、2700 m 高度的后向軌跡。Fig. 5 Backward trajectories of air masses in Taiyuan station on (a) 3 August, (b) 4 August, (c) 12 August, (d) 30 August 2013, (e) 11 August 2014 and (f) 12 August 2014 The red, blue, green, and cyan solid lines represent the backward trajectories at 500 m, 1200 m, 2200 m, and 2700 m height,respectively.
可以看出,2013年8月3日、8月4日和8月12日三次大氣探測過程中,低層氣溶膠主要來自本地源,以城市群人類活動排放源為主,而高層氣溶膠來自西北地區,其中2013年8月3日來自柴達木盆地,2013年8月4日和2013年8月12日來自新疆地區,經過戈壁沙漠之后進入山西。因此,由于這三次探測期間高層氣溶膠來自西北地區遠距離輸送的尺度較大沙塵粒子,而低層氣溶膠主要來自局地城市群人類活動排放的尺度較小的粒子,又因為有邊界層逆溫阻礙垂直輸送,氣溶膠垂直混合受阻,造成氣溶膠有效直徑隨高度增加呈現增大的趨勢。一般情況下,沙塵粒子為非吸濕性氣溶膠粒子,因此,核化成CCN 的濃度較小。2013年8月30日低層氣溶膠來自東北地區輸送,高層氣溶膠來自蒙古地區的輸送。2014年8月11日低層氣溶膠來自蒙古地區而高層氣溶膠來自西北地區。2014年8月12日,低層氣溶膠來自局地的排放,高層氣溶膠來自西北地區。可見2013年8月30日和2014年8月11日低層氣溶膠有效直徑較大,是因為低層氣溶膠有來自其他地方的輸送,而三次觀測的高層氣溶膠均來自西北或蒙古地區的遠距離輸送,而高層氣溶膠有效直徑比較小,這可能是由于云對氣溶膠的濕清除作用,尺度較大的氣溶膠容易形成云滴,因此,氣溶膠有效直徑出現隨高度增加呈現減小或不變的趨勢。
利用2013年7月31日、8月3日、8月4日三次在山西文水地區不同高度層飛機水平盤旋探測數據,分析了不同高度層的氣溶膠粒子譜分布。其中2013年7月31日探測時,空中有一層較薄的云層,2013年8月3日、2013年8月4日為晴空天氣。
這三次飛行探測的方案是,在山西文水地區從最低飛行安全高度1200 m 開始,以300 m 或600 m 為間隔分層水平盤旋飛行探測,在每一個高度層大約飛行5 min。圖6 是三次探測相同高度盤旋平飛時的氣溶膠譜。圖6a 中7月31日、8月4日、8月3日氣溶膠譜均呈單峰分布,不同日期相同尺度的氣溶膠濃度相差較大,圖6b–d 中8月3日和8月4日氣溶膠譜呈雙峰分布,而7月31日氣溶膠譜呈單峰分布,圖6c–d 中7月31日氣溶膠譜與8月3日和8月4日氣溶膠譜型差異較大,直徑小于1.2 μm 的氣溶膠濃度與8月3日和8月4日相比較小,而直徑大于1.2 μm 的氣溶膠濃度較大,這是由于7月31日平飛探測時,飛機在1800~2500 m 位于云中,由于濕增長作業而使得氣溶膠粒子增大;圖6e–f 中三次探測氣溶膠譜型較為一致,都呈單峰分布,說明不同日期3000 m 以上氣溶膠分布差異較小。

圖6 2013年7月31日、8月3日、8月4日不同高度氣溶膠粒子譜分布:(a)1200 m;(b)1500 m;(c)1800 m;(d)2500 m;(e)3100 m;(f)3700 m。圖中N 和D 分別表示粒子的數濃度和直徑Fig. 6 Aerosol particle size distributions at (a) 1200 m, (b) 1500 m, (c) 1800 m, (d) 2500 m, (e) 3100 m, and (f) 3700 m on July 31, August 3, August 4, 2013. N, D represent number concentration and diameter of particle, respectively.
綜上,不同日期氣溶膠譜型在低層差異較大而高層基本一致;因為低層受人為排放的氣溶膠影響較大,造成不同日期氣溶膠譜型相差較大;而不同高度的氣溶膠譜型不同,1200 m 氣溶膠譜呈單峰或雙峰分布,1500~2500 m 氣溶膠譜呈雙峰分布,3100~3700 m,氣溶膠譜呈單峰分布;云中氣溶膠譜呈單峰分布,與晴空相比,直徑小于1.2 μm 的氣溶膠濃度小,而直徑大于1.2 μm 的氣溶膠濃度較大。
圖7 為不同高度的平均氣溶膠譜分布,可以看出,氣溶膠粒子主要集中在小粒子端,大粒子端的Na較小;低層各粒徑段的Na都比較高,隨著高度的增加各粒徑段的Na都有所減小,特別是小粒子端,Na減小非常明顯,從而使高層氣溶膠中大粒子所占的比例有所增加,從后向軌跡來看,觀測期間低層的氣溶膠粒子主要來自近地面的局地排放;而高層的氣溶膠粒子來自西北地區的遠距離輸送,高層氣溶膠大粒子所占比例增加,就是由于來自西北地區沙塵氣溶膠的影響,但整個氣溶膠譜型沒有明顯變化,說明大粒子氣溶膠對本地氣溶膠譜型的影響較小。

圖7 不同高度平均氣溶膠譜分布Fig. 7 Averaged aerosol particle size distributions at different heights
由以上分析可知,氣溶膠數濃度與CCN 數濃度的分布趨勢較為一致,因此我們做了二者的相關性分析。圖8 所示是5架次的飛行探測過程中,氣溶膠數濃度與CCN 數濃度(過飽和度為0.3%,表3)相關性分析。CCN 數濃度與氣溶膠數濃度有較好的相關性,決定系數R2可達到0.8 以上,決定系數R2達到0.9 以上的有兩次,分別是2013年8月3日和2013年8月4日的飛行探測過程,決定系數R2分別為0.91 和0.96。其余三次的R2都在0.8~0.9 之間,為此我們統計了不同日期氣溶膠的平均有效直徑,如表4 所示,可以看出,2013年8月3日和2013年8月4日探測過程中,氣溶膠粒子的平均有效直徑較大,均大于0.4 μm,而其余三日氣溶膠粒子的平均有效直徑都小于0.4 μm,因此,華北中部夏季垂直方向上,氣溶膠數濃度與CCN數濃度有較好的相關性,并且當氣溶膠粒子平均有效直徑大于0.4 μm時,氣溶膠數濃度與CCN 數濃度的相關性更好。

圖8 (a) 2013年8月3日, (b) 2013年8月4日, (c) 2013年8月12日, (d) 2013年8月30日和(e) 2014年8月12日飛機探測過程中CCN 數濃度與氣溶膠數濃度的相關性Fig. 8 Linear regressions between the CCN number concentration and Na in flights of (a) 3 August, 2013, (b) 4 August, 2013, (c) 12 August, 2013 (d)30 August, 2013 and (e) 12 August, 2014

表3 不同高度CCN 數濃度(NCCN;過飽和度為0.3%)統計特征Table 3 Statistical summary of the CCN number concentration (NCCN)at 0.3% supersaturation at different heights in each case

表4 垂直探測過程中不同日期氣溶膠數濃度與CCN 數濃度的決定系數及氣溶膠有效直徑的統計特征Table 4 Statistical properties for the correlation coefficient values between the CCN number concentration, Na, and ED
氣溶膠成為CCN 的能力與尺度和化學組分有很大的關系。在不了解氣溶膠詳細化學組分的情況下,為了進一步了解氣溶膠尺度與氣溶膠活化能力的關系,我們分析了0.3%過飽和度條件下,CCN比率(CCN/CN,CN 為“凝結核”的簡稱)與ED 的關系。圖9 為2013年8月03日飛行探測時,整個探測過程中CCN/CN 與ED 的關系圖。可以看出,CCN/CN 與ED 呈線性關系,隨著ED 的增大,CCN/CN 也增大,當ED 為0.25 μm時,CCN/CN大約為0.25,即約25%的氣溶膠可以作為CCN,當ED 為0.6 μm時,CCN/CN 大約為0.45,即約45%的氣溶膠可以作為CCN,當ED 為0.9 μm時,CCN/CN 大約為0.65,即約65%的氣溶膠可以作為CCN,表明尺度越大的氣溶膠越容易成為CCN。

圖9 2013年8月3日飛行探測過程中CCN 占CN 的比例(CCN/CN)與ED 的關系Fig. 9 Relationship between the CCN to CN ratio (CCN/CN) and ED during the flight on 3 August 2013
為了研究不同高度氣溶膠的活化能力,分析了低層和高層的CCN/CN 與ED 的關系(圖10)。圖10a 為低層(1500 m 以下)CCN/CN 與ED 的關系圖,圖中可以看出,低層氣溶膠ED 值較小,大部分集中在0.5 μm 以下,低層CCN/CN 與ED 也呈線性關系,隨著ED 的增大,CCN/CN 也增大,當ED 為0.25 μm時,CCN/CN 大約為0.2,即約20%的氣溶膠可以作為CCN,當ED 為0.4 μm時,CCN/CN 大約為0.35,即約35%的氣溶膠可以作為CCN;同樣,圖10b 為高層(3000~4000 m)CCN/CN 與ED 的關系圖,圖中可以看出,高層氣溶膠ED 值較大,高層CCN/CN 也與ED 呈線性關系,隨著ED 的增大,CCN/CN 也增大,當ED 為0.4 μm時,CCN/CN 大約為0.4,即約40%的氣溶膠可以作為CCN,當ED 為0.9 μm時,CCN/CN大約為0.65,即約65%的氣溶膠可以作為CCN。

圖10 2013年8月3日飛行探測過程(a)1500m 以下以及(b)3000~4000 m 高度CCN/CN 與ED 的關系Fig. 10 Relationships between the CCN/CN and ED at different heights during the flight on 3 August 2013: (a) 0–1500 m, (b) 3000–4000 m
由以上分析可以得出,CCN/CN 隨著高度增加而增大。CCN 活化率受到氣溶膠粒子尺度和氣溶膠化學成分的影響,Fitzgerald(1973)和Kuwat etal.(2008)等研究表明粒子尺度對活化率的影響高于化學成分的影響。本文中不涉及化學成分的研究,從物理的角度來講,過飽和度較高的情況下,CCN/CN 受到粒子尺度影響很大。故雖然低層和高層粒子來源不同,但由于高層粒子尺度較大,受到粒子尺度的影響,高層CCN/CN 的值較大。
CCN 活化譜是CCN 數濃度隨過飽和度的變化曲線,CCN 數濃度觀測的一個很重要的應用就是CCN 活化譜。通常用經驗公式N=CSK來擬合,其中,C、K為擬合參數,N為過飽和度S下的CCN數濃度,氣溶膠粒子的尺度或化學組分的信息隱含在擬合參數C、K之中。根據擬合參數C、K的值,把CCN 活化譜分為海洋型(C<1000 cm?3,K<1)、過渡型(1000 cm?3
云下氣溶膠粒子能否成為CCN,對云的形成及云的微物理特征有重要的意義。在以往的研究中,很少有在云底之下300 m 范圍內直接對氣溶膠和CCN 進行觀測,所以本次對CCN 活化譜的研究所選取的觀測位置是位于云底之下100~300m 范圍和云外的區域,圖11 中2014年7月15日是在云外做的觀測,高度3700 m 做水平飛行,2014年7月19日是在云底之下100 m(高度為3100 m)做云底之下保持高度水平盤旋飛行,過飽和度分別設置為0.4、0.6、0.8。擬合參數C、K值如表5 所示,2014年7月15日的云外觀測屬于過度型活化譜,過飽和度分別設置為0.2、0.4、0.6。2014年7月19日的云底之下觀測屬于清潔大陸型活化譜,而封秋娟等(2012)對山西省地面做的CCN 活化譜屬于典型大陸型活化譜。

圖11 云外及云下CCN 活化譜Fig. 11 CCN spectra were observed below cloud bases and the environment.

表5 CCN 活化譜擬合參數值Table 5 Parameters of CCN Spectrum
氣溶膠不僅對云的形成及云的宏微觀特征有重要的影響,氣溶膠還可以通過影響云的宏微觀特征而影響氣候,而氣溶膠對云的影響又具有很大的不確定性,本研究想通過對云下氣溶膠粒子與云底之上的云滴粒子進行觀測,確定氣溶膠粒子與云滴粒子的轉化關系。氣溶膠粒子能否活化成為云滴粒子,與氣溶膠粒子的物理化學性質、環境過飽和度以及上升速度等有關系,本研究只做了Na、CCN 濃度和Nc的觀測分析。
利用2014年8月11日的觀測資料,分析了Na與Nc以及CCN 數濃度與Nc之間的關系。本次觀測的飛行方式是在云底之下200 m(高度為3400 m)做保持高度水平盤旋飛行,大約水平盤旋5 min 之后從云底入云,在云底之上100 m(高度3700 m)做保持高度水平盤旋飛行。在云底之下盤旋平飛為了獲取氣溶膠和CCN 的觀測資料,在云底之上盤旋平飛為了獲取云滴粒子的觀測資料。在不了解實際過飽和度、上升速度以及氣溶膠化學組分的情況下,我們希望在氣溶膠和云滴之間找到一種簡單的線性關系,如圖12 所示,擬合方程為y=1.3x?616.3,擬合相關系數為0.96,計算得到氣溶膠轉化為云滴的比率為47%。

圖12 云底之下200 m 高度層Na 與云底之上100 m 高度層云滴數濃度Nc 的關系Fig. 12 Relationship between Na and Nc (cloud droplet number concentration) near the cumulus cloud base
此外,本次研究還做了在固定過飽和度0.3%條件下,CCN 數濃度與Nc之間的關系,如圖13所示,擬合方程為y=1.6x?473.8,擬合相關系數也為0.96,計算得到CCN 轉化為云滴的比率為69%。

圖13 過飽和度0.3%條件下云底之下200 m 高度層CCN 數濃度與云底之上100 m 高度層Nc 的關系Fig. 13 Relationship between the CCN number concentration (0.3%supersaturation) and Nc near the cumulus cloud base
本文分析了山西省9架次的夏季云物理飛機探測資料,得到如下主要結論:
(1)大氣邊界層逆溫層結對大氣氣溶膠粒子的垂直分布有重要的影響。當大氣邊界層有逆溫存在時,邊界層內Na高且ED 在邊界層內比較均勻;當大氣邊界層無逆溫層結時,氣溶膠粒子可隨大氣垂直輸送到較高的高度,邊界層內Na較小。有逆溫時,邊界層內Na大約為無逆溫時的2 倍;兩類Na的垂直分布都呈隨高度增加而減小的趨勢。相同日期里CCN 數濃度的垂直分布趨勢與Na的垂直分布趨勢基本一致,有逆溫時,邊界層內CCN 數濃度大約為無逆溫時1.5 倍。
(2)華北中部夏季高層一般受到來自西北或蒙古等地的氣團影響,造成晴空時高層氣溶膠尺度大,ED 在0.45 μm 以上,有云時,高層尺度較大的氣溶膠活化成為云滴,造成高層氣溶膠尺度小,ED 在0.35 μm 以下;低層氣溶膠以局地排放為主時,氣溶膠尺度小,ED 在0.4 μm 以下。
(3)氣溶膠譜型在低層差異較大而高層基本一致;因為低層受局地排放的氣溶膠影響較大,造成不同日期氣溶膠譜型相差較大;而不同高度的氣溶膠譜型不同,1200 m 氣溶膠譜呈單峰分布,1500~2500m 氣溶膠譜呈雙峰分布,3100~3700 m,氣溶膠譜呈單峰分布;云中氣溶膠譜呈單峰分布。與晴空相比,直徑小于1.2 μm 的氣溶膠濃度小,而直徑大于1.2 μm 的氣溶膠濃度大。低層氣溶膠各粒徑段的粒子濃度均高于高層,高層氣溶膠由于受到來自西北地區沙塵氣溶膠的影響,大粒子所占比例增加,但整個氣溶膠譜型沒有明顯變化,說明大粒子氣溶膠對本地氣溶膠譜型的影響較小。
(4)垂直方向上,氣溶膠數濃度與CCN 數濃度有較好的相關性,相關系數都在0.8 以上,并且當ED 大于0.4 μm時,氣溶膠數濃度與CCN 數濃度的相關性好于氣溶膠粒子平均有效直徑小于0.4 μm 的情況。2013年8月03日整個探測過程中CCN/CN 與ED 呈線性關系,隨著ED 的增大,CCN/CN 也增大,并且低層、高層均表現出這種關系。氣溶膠粒子尺度是影響CCN/CN 的重要因素。
(5)用N=CSK擬合得到的山西地區空中CCN 活化譜參數,云外屬于過度型核譜,云底之下300 m 范圍內屬于清潔大陸型核譜。
(6)積云云下氣溶膠與云滴的線性擬合方程為y=1.3x?616.3,擬合決定系數為0.96,氣溶膠轉化為云滴的比率為47%。過飽和度0.3%條件下,云下CCN 與云滴的線性擬合方程為y=1.6x?473.8,擬合決定系數也為0.96,CCN 轉化為云滴的比率為69%。