海陸風環流中海鹽氣溶膠對大氣影響的模擬

沈 傲,田春艷,劉一鳴,申 沖,陳曉陽,陳訓來,王明潔,藍 靜,樊 琦*

?

海陸風環流中海鹽氣溶膠對大氣影響的模擬

沈 傲1,田春艷1,劉一鳴1,申 沖1,陳曉陽1,陳訓來2,王明潔2,藍 靜3,樊 琦1*

(1.中山大學大氣科學學院,廣東省氣候變化與自然災害研究重點實驗室,廣東 廣州 510275；2.深圳市國家氣候觀象臺,廣東 深圳 518040；3.廣州市氣候與農業氣象中心,廣東 廣州 511430)

采用WRF-CMAQ模式對2013年9月9~13日珠江口一次典型海陸風過程進行數值模擬,并進一步探討了在此類型環流系統下，海鹽氣溶膠及其對大氣環境的影響.結果表明:珠江口附近海風環流一般在14:00左右開始形成,17:00左右發展最旺盛,且在珠江口沿岸形成扇形風系.在海陸風環流影響下,珠江三角洲地區陸地上Na+濃度從17:00開始增加,次日02:00左右達到最大,10:00以后開始下降.海鹽氣溶膠存在明顯的氯虧損現象,其中白天的氯虧損情況遠比晚上嚴重,細粒徑的氯虧損情況比粗粒徑嚴重.設置有無海鹽氣溶膠排放的敏感性試驗進行對比研究發現:考慮海鹽排放后,珠江口地區HCl,SO42-,NO3-的濃度均有增加,其中,HCl和NO3-的增加較多,SO42-較少.HCl濃度增加的高值在Cl-高值之前;SO42-濃度變化主要受海鹽氣溶膠排放中的硫酸鹽組分影響;NO3-濃度變化則主要受氯虧損機制和NO3-干沉降過程共同影響.

珠江口；海陸風環流；海鹽氣溶膠；氯虧損；數值模擬

海陸風是發生在海岸附近的由于海陸熱力差異導致的一種局地環流,具有明顯的日變化[1].國內外學者對此現象進行了大量的研究[2-3].近些年來,隨著數值模擬技術的不斷提高,MM5、WRF等中尺度氣象模式得到廣泛運用,為海陸風的研究提供了良好的工具.部分學者對國內外各地海陸風環流進行數值模擬和檢驗分析,均取得了較好的研究成果[4-7].研究表明,對于珠江三角洲地區,海陸風發生頻率夏季大于冬季,并且海風開始和結束的時間都較晚[8].張立鳳等[9]由Mass模式的診斷結果得到了珠江口地區海陸風風系分布的概念模型:吹海風時,珠江口沿岸存在扇形風系和扇形輻合帶;吹陸風時,珠江口海面上存在東西兩岸陸風的輻合. Ding等[10]和Lu等[11]分別利用MM5模式對2001年9月和2003年整年的珠江三角洲地區海陸風變化進行模擬,均發現,夏秋季珠江口的海風形成時間較晚,一般為14:00以后,且珠江口會有較明顯的扇形輻散風系.

海陸風環流會對沿海地區大氣中顆粒物的輸送和反應有較大的影響[12-13],其中海鹽氣溶膠(SSA)是大氣中最主要的自然源氣溶膠之一[14],尤其是在海區和沿海大氣,也是對流層中最大的顆粒物來源[15].據估計,海鹽氣溶膠的全球通量每年為1′109~ 1′1010g[16].根據IPCC在2001年的報告估計,每年由于作用于海洋表面的風應力,大概有3.3′109g的海鹽氣溶膠進入海洋大氣邊界層[17].

海鹽氣溶膠主要來源于海水飛沫,其主要化學組分也與海水相似,主要為Cl-,Na+,Mg2+,Ca2+, K+,S及Br-,且其關系比值也與海水中的相應比值接近.海鹽氣溶膠除了本身的排放外,還會因為氯化學反應過程增加沿海地區顆粒物的濃度,改變顆粒物化學組分.其中Cl-顆粒物易與H2SO4,HNO3等反應生成HCl氣體而發生氯虧損現象:

式中:g代表氣態,cd代表凝聚態.

國內外的研究中,學者們利用觀測資料、實驗室研究和數值模擬等方法,均發現海鹽氣溶膠會對空氣中的部分物質產生重要的影響[18-21].研究表明,海鹽氣溶膠會使海洋和沿海地區的硫酸鹽、硝酸鹽、氨氣和PM10濃度增加,而硫酸氣體、硝酸氣體、銨鹽濃度減小,其中硫酸鹽顆粒物的一次排放引起的,硝酸鹽濃度的增加主要是由海鹽氣溶膠的氯化學反應引起的,且硝酸鹽和硝酸氣體的改變最為明顯.

目前針對沿海地區特有的海陸風環流系統對海鹽氣溶膠的輸送及其影響方面的研究較少.本文選取2013年9月9~13日珠江口(PRE)地區的海陸風過程,利用三維空氣質量模式CMAQ對此過程進行數值模擬,分析海陸風過程中珠江三角洲地區的氣象場特征、海鹽氣溶膠的輸送及其對沿海地區大氣環境的影響.

1 研究方法

1.1 模式及其設置

本文采用WRF-CMAQ模式系統:由中尺度模式WRF[22]提供氣象場數據,并由MCIP[23]轉換程序把WRF輸出結果轉化成為空氣質量模式CMAQ[24]能讀取的文件格式.排放源模塊同時考慮了人為排放源和生物排放源.CMAQ利用氣象模塊和排放源模塊的輸入結果,對多種污染物的濃度進行計算.其中CMAQ的海鹽氣溶膠模塊的優勢在于,一方面考慮了海溫的影響,另一方面更加注重海鹽細顆粒物的模擬效果[25].CMAQ中海鹽排放按來源區域分為廣闊海洋面排放和沿海碎浪帶排放.沿海碎浪帶的海鹽排放要明顯高于廣闊海洋面上的排放.對于廣闊海洋面,海鹽排放量計算公式參考文獻[26],而對于沿海碎浪帶,海鹽排放量的計算參考文獻[27],即在文獻[26]公式的基礎上假設沿海碎浪帶的白浪占比為100%,且碎浪帶的寬度為25m.在計算海鹽排放時,需要對模擬區域中每個網格的廣闊海洋面和沿海碎浪帶的比例進行計算.本文利用ArcGIS軟件從海岸線信息中獲取這些數據.

對于WRF-CMAQ系統模擬能力的評估已得到了驗證[21],說明WRF-CMAQ系統總體上能較好地反映實際大氣情況,為之后的分析提供較可靠的氣象場和污染物濃度結果.

本文使用的WRF版本為v3.6.1,模式采用蘭伯特投影,中心位于(110°E,30°N).模擬區域采用三重網格嵌套(圖1),水平格距分別為27,9, 3km,格點數為240×240,148×160,199×175.垂直方向分為27層,在大氣邊界層以下劃分更細致一些.地形場資料采用的是歐盟聯合研究中心(JRC)制作的GLC2009[28],分辨率為30s.氣象初始場和邊界場資料采用美國國家環境預報中心(NCEP)/美國國家大氣研究中心(NCAR)的再分析資料,空間分辨率為1°×1°,時間間隔為6h.海溫數據來自于NCEP根據衛星數據反演得到的全球實時(RTG)海溫數據,空間分辨率為0.5°×0.5°.WRF模式采用的物理參數化方案包括:Lin微物理方案,Goddard短波輻射方案,RRTM長波輻射方案,Monin-Obukhov近地面層方案,Noah陸面過程方案,YSU邊界層方案,Kain-Fritsch積云參數化方案,其中第三重網格不采用積云參數化方案,并采用了單層城市冠層模式(UCM).模擬時間為2013年9月7日08:00~9月13日08:00,其中前48h為啟動時間.

本文使用的CMAQ版本為v5.1,模擬區域采用三重網格嵌套模式(圖1),水平格距與WRF模式設置相同,格點數分別為119′129,102′81,153′135. 垂直方向分為27層,在大氣邊界層以下劃分更細致一些,最底層高度為22m.內重網格的化學邊界條件由外重網格的模擬輸出結果提供.模式模擬第1d的化學初始條件由模式默認的結果提供,第2d的化學初始條件由第1d的模擬結果提供.CMAQ模式采用的化學參數化方案包括:SAPRC07TIC氣相化學機制,AERO6I氣溶膠機制,AQKMTI液相化學機制,其中光化學反應速率根據模式模擬的氣溶膠和臭氧濃度在線計算. 模擬時間為2013年9月8日08:00~ 9月13日08:00,其中前24h為啟動時間.

圖1 WRF-CMAQ模式網格范圍示意

排放源部分,人為源排放采用清華大學制作的MIX2010排放源清單[29],生物排放源由MEGAN[30]模式計算得到.

1.2 敏感性試驗設置

為具體研究海鹽氣溶膠對大氣環境的影響,本文采用設置敏感性試驗的方法定量得到海鹽氣溶膠對大氣環境的貢獻.基準試驗(SACL)考慮了海鹽的排放量,而控制試驗(NOSA)不考慮海鹽的排放量,對比2個數值試驗的模擬結果即可得到海鹽氣溶膠對大氣環境的影響.

2 結果分析

2.1 海陸風特征的模擬分析

沿海地區海陸風會受到多重因素的影響.本研究綜合不同學者對海陸風的研究以及研究區域的位置特點,對于海陸風的識別方法采用:風分量為正值時,為向岸流即海風;為負值時,為離岸流即陸風.在24h內,風分量發生變化且持續時間不少于4h,則認為有海風或陸風發生.風分量持續時間為陸風或海風持續時間[31].

取113°45′E,22°30′N(珠江口海上)做垂直層次上風沿時間變化的剖面圖(圖2),可以發現,風具有明顯的日變化特征.這4d內,海風環流一般在午后14:00左右開始形成,17:00~18:00左右發展最旺盛,風風速最大可達8m/s.海風一般可維持至次日凌晨02:00左右,其中10/11日海風維持到了次日06:00左右,且11日的風強度在這4d中達到最大.相比海風,陸風的維持時間較短,一般僅能維持6h左右,強度也較小,風風速一般為2~3m/s.垂直方向上,海風影響范圍可達1000m,并且強度隨高度的增加而減小.

圖2 模式模擬的2013年9月9~13日113°45′E 22°30′N 垂直方向上V風隨時間變化

以9月11~12日海陸風過程為例,選取典型時刻(9月11日10:00、14:00、17:00,9月12日00:00、06:00、10:00)做珠江口10m風場與2m溫度場疊加圖(圖3).可以發現,11日10:00,珠江口為較弱的離岸陸風,此時珠江口風速接近于0,陸地上風速也較小,約為2~3m/s,且無主導風向.但陸上溫度已開始高于海上溫度,溫差約為2~3℃.下午14:00,珠江口開始形成向岸海風,海上風速最大達6m/s,海陸溫差達8℃.由于珠江口特殊的地形,海風在珠江口處形成扇形輻散風系,這與張立鳳等[9]的研究結果一致,且由于右岸溫度略高于左岸溫度,風系的西風分量要略強于東風分量.至17:00,海風發展最為旺盛,海上風速最大達8m/s,陸上風速最大達6m/s.此時珠江口東西兩岸溫度一致,風系的東西風分量也一致.隨后,海風一直維持至12日凌晨6:00,此時由于陸上溫度已普遍低于海上溫度3℃左右,海風逐漸向陸風轉變,并且由于海上的背景主導風,珠江口出現較強的東風分量.至10:00,海風已完全轉變為陸風.其余幾日與此類似,故不再贅述.

圖4 模式模擬的2013年9月11日17:00沿113°45′E V風垂直剖面

選取11~12日海風強度最大時(11日17:00),沿113°45′E做風的垂直剖面圖(圖4),可以發現,地面和高空存在明顯的反向氣流,且海風即將結束前風場分量值有方向相反對峙區域,說明了海風對內陸的滲透距離.由圖4可知,海風強度發展至最大時,海風可推進至23°5′N左右.由于珠江口位于22°45′N,可得知海風向陸地推進了20′,即為向陸地推進將近37km,可影響到廣州中南部地區.

2.2 海陸風環流下海鹽氣溶膠的模擬分析

在CMAQ模式中,海鹽氣溶膠各組分的占比如表1所示,其中氯鹽和鈉鹽是海鹽中的主要成分,各占55%和31%.而Na+化學性質不活躍,故常用于表征海鹽氣溶膠.模式中,海鹽排放在海面上完成,物質的平流過程依賴于連續方程的質量守恒特性,而在垂直擴散上,模式采用的是非對稱對流法ACM[32]的更新版本ACM2[33].在對流條件下(表面變暖時),被加熱的空氣通過浮力垂直運輸,并與表面上方每個層次的環境空氣混合,直到上升空氣的溫度等于周圍環境的溫度.由于該過程是由快速上升的窄氣流和緩慢下降的寬氣流所導致,所以垂直擴散是非對稱的,在CMAQ計算期間通過使用ACM2,在線垂直分配生物源和點源排放.在非對流條件下(表面冷卻時),垂直擴散由渦流擴散方法表示.渦流擴散是一種局部混合方案,并且使用與區域酸沉積模型[34]中相同的基于行星邊界層(PBL)相似性的算法估算.

表1 CMAQ模式中海鹽氣溶膠干顆粒組分占比(%)

2.2.1 基準試驗中海鹽氣溶膠的分布 選取113°35′E 23°N(廣州南部)做垂直方向上環境中Na+質量濃度隨時間的變化(圖5).由于在CMAQ模式中Na+根據粒徑大小不同分為核膜態、積聚態、粗粒態3種模態,其中粗粒態占比可達90%以上,故僅研究3種模態下Na+質量濃度的總和.可以發現,陸地上Na+本底濃度極低,一般不超過1μg/m3.在17:00左右,也就是海風強度最大時,該地Na+濃度開始增加,一般在第二日02:00前后達到峰值,最大可達7μg/m3.在陸風形成后,該地Na+濃度逐漸降低.在垂直方向上,Na+濃度隨高度的增加而減小,最高可超過2500m,濃度相對較大區(33μg/m3)一般分布在500m以下高度處.

選取典型時刻(9月11日10:00,14:00,17:00;9月12日00:00,06:00,10:00)對模式最底層的珠江三角洲(PRD)風向與Na+質量濃度做疊加圖(圖6).可以發現,在11日10:00,由于盛行東風,珠江三角洲西部Na+質量濃度比東部高.14:00,海風剛開始生成,由于之前持續的陸風,陸上Na+濃度較低,此時海上Na+的濃度可達10μg/m3.17:00,海風強度達到最大,此時沿海陸地Na+濃度開始增加.此后,Na+一直向陸地推進.在12日00:00,整個海岸線地區Na+濃度達到最大,普遍可達10μg/m3以上,廣州地區Na+可推進至24°45′左右,南部地區濃度較高.直至12日06:00,Na+推進至最遠,中心濃度較大區數值可超過10μg/m3.并且由于此時的背景風為東風或偏東風,珠江口西部陸地Na+濃度普遍高于東部,廣州南部Na+濃度也較高.12日10:00,珠江三角洲地區已全被陸風控制,陸地上Na+濃度也逐漸降至1μg/m3左右.

圖5 模式模擬的2013年9月9~13日113°35′E 23°N 垂直方向上Na+質量濃度隨時間變化

2.2.2 基準試驗中海鹽氣溶膠的氯虧損現象 選取113°35′E 23°N(廣州南部)做3種模態下Na+與Cl-在模式最底層的質量濃度總和隨時間的變化(圖7).可以發現,Cl-和Na+的變化基本同步,說明該地Cl-主要來源于海鹽.除了每日00:00前后Cl-的質量濃度要高于Na+,其余時間均會低于Na+的質量濃度.并且在每日00:00前后,Cl-濃度達到峰值,16:00前后達到谷值.而排放初期的海鹽氣溶膠和海水成分一致,Cl-和Na+的摩爾比為常數1.16.對9月9~13日Cl-,Na+,Ca2+分別做摩爾濃度的比值,其中[Cl-]/[Na+]若小于1.16,說明出現了氯虧損狀況;[Na+]/[Ca2+]若大于6,說明海洋源占優勢[35].可以發現,9月9~13日,絕大多數時間海洋源氣溶膠占優勢,并且每日00:00前后,海洋氣溶膠的濃度達到最大.此地一直出現較為明顯的氯虧損現象,并且白天的氯虧損狀況([Cl-]/ [Na+]約為0.1左右)要遠比晚上([Cl-]/[Na+]約為0.85左右)嚴重,這與Cl-濃度變化曲線趨勢一致.

由9月9~13日3個模態(核膜態、積聚態、粗模態)下的摩爾濃度比值[Cl-]/[Na+](圖8)可以看出,由于粗粒態下顆粒物的占比可達總顆粒物的90%以上,所以該模態下[Cl-]/[Na+]比值變化與3模態總和下的變化十分接近,但在核膜態和積聚態下,[Cl-]/ [Na+]比值要遠小于粗粒態下的比值,且其變化趨勢也與粗粒態下的相差甚遠,說明細粒徑下的氯虧損情況要比粗粒徑下的嚴重.究其原因,由于細粒子在大氣中停留時間長,比表面積大,發生反應的幾率大,因而氯虧損可能也較多[36].

圖8 2013年9月9~13日113°35′E 23°N模式最底層摩爾濃度比值[Cl-]/[Na+]

2.3 海鹽氣溶膠對大氣環境的影響

由式(1)、(2)可知，Cl-會與酸性氣體反應生成二次無機鹽和HCl氣體釋放到大氣當中,結果就是Cl-濃度減小,HCl濃度增加.

選取113°35′E 23°N(廣州南部)做Cl-質量濃度差值和HCl體積濃度差值以及NOSA試驗模擬下的HCl體積濃度(圖9(a)).可以發現,大氣中HCl的本底濃度極低,一般僅能維持在0.2′10-9左右.考慮海鹽排放后,該地HCl濃度最多可增加2.5′10-9,能占大氣中HCl總和的92.3% (12日18:00).Cl-濃度增加的高值分別在9日23:00,11日01:00,12日00:00,13日00:00.而HCl濃度增加的高值分別在9日19:00,10日18:00,11日18:00,12日18:00.再選取113°45′E 22°30′N(珠江口海上)同樣做Cl-質量濃度差值和HCl體積濃度差值和NOSA試驗模擬下的HCl體積濃度(圖9(b)),可以發現,海上HCl濃度增加的高值分別在9日15:00,10日16:00,11日14:00,12日14:00,均處在陸風轉變為海風的交界時.可見,在陸風盛行時,陸地上許多污染物輸送至珠江口海面上.當陸風轉變為海風時,離陸地較遠的新鮮海鹽氣溶膠輸送至珠江口海面處,并在此通過氯化學反應生成HCl.18:00前后,由于此時珠江口海風強度達到最大,并且HCl作為氣體比海鹽顆粒輕,所以海上堆積的HCl又重被輸送至陸地.到第二日00:00前后, Cl-輸送至內陸,HCl濃度增加也一直維持在1′10-9左右.

選取113°35′E 23°N(廣州南部)分別做Cl-和SO42-質量濃度差值以及NOSA試驗模擬下的SO42-質量濃度(圖10).可以發現,大氣中SO42-本底濃度較高,而考慮海鹽排放后,SO42-濃度的增加量僅占大氣中總SO42-的8.5%左右.由圖也可看出,SO42-質量濃度差值的變化趨勢與Cl-大致相似,兩者的相關系數達0.82.由表1和公式(2)可知, SO42-濃度的增加量由2部分組成,一是海鹽氣溶膠中硫酸鹽組分的直接排放,二是海鹽氣溶膠中的Cl-與大陸污染氣團中的硫酸氣體發生反應生成SO42-,但二者總體增加的量很少.SO42-濃度的增加中氯化學反應占的比重很小,主要是硫酸鹽組分的直接排放.

圖10 2013年9月9~13日113°35′E 23°N模式最底層SACL和NOSA試驗模擬的Cl-與SO42-質量濃度差值和NOSA試驗模擬的SO42-質量濃度

選取113°35′E 23°N(廣州南部)分別做Cl-和NO3-質量濃度差值以及NOSA試驗模擬下的NO3-質量濃度(圖11(a)).可以發現,大氣中NO3-本底濃度大多數時間極低,考慮海鹽排放后NO3-濃度的增加量占大氣中總NO3-的比值最大可達100%(9日16:00,10日22:00,11日15:00),在數值上NO3-濃度最多可增加10.36μg/m3(13日07:00),NO3-濃度增加量占總NO3-的比值最小也能達46.1%(10日07:00),說明大氣中NO3-的變化絕大多數是由海鹽氣溶膠的氯虧損機制造成的.對比圖10,NO3-質量濃度的增加量明顯大于SO42-的增加量,可見,在氯虧損機制中, HNO3的置換反應要強于H2SO4的置換反應.在變化趨勢上,NO3-濃度增加量的高值在Cl-高值前后均有出現.再選取113°45′E 22°30′N(珠江口海上)同樣做Cl-和NO3-質量濃度差值和NOSA試驗模擬下的NO3-質量濃度(圖11(b)),可以發現,同HCl一樣,珠江口海面上出現在Cl-高值前的NO3-高值也是出現在陸風轉變成海風的時段,類似圖9的分析,可以認為盛行海風將珠江口堆積的NO3-輸送至陸地.出現在Cl-高值后的NO3-高值是由于盛行東風將臺灣地區的NO3-輸送至珠江口,隨后通過海風輸送至陸地,在此不作討論.在高值后NO3-質量濃度的增加趨勢隨后會顯著減小,是由于海鹽氣溶膠粒徑較大,NO3-與其混合后受重力影響沉降到地面.從圖11(c)可以看出,NO3-干沉降量的增加與NO3-質量濃度的減小有較好的一致性.

3 結論

3.1 珠江口海風環流一般在午后14:00左右開始形成,17:00~18:00左右發展最旺盛,強度較大,一般可維持至第二日凌晨2:00左右.陸風的維持時間較短,一般為維持6h左右,強度也較小.垂直方向上,海風影響范圍可達1000m.海風發生時,珠江口處形成扇形輻散風系.

3.2 用Na+表征海鹽氣溶膠來看,在17:00左右,沿海地區Na+濃度開始增加,第二日2:00前后達到峰值.在陸風形成后,Na+濃度逐漸降低,垂直方向上, Na+最高可影響到2500m以上高度,濃度相對較大區(33μg/m3)分布在500m高度以下.內陸存在氯虧損現象,白天的氯虧損情況遠比晚上嚴重,細粒徑的氯虧損情況比粗粒徑的嚴重.

3.3 考慮海鹽排放后,珠江三角洲地區HCl,SO42-, NO3-的濃度均有增加,其中,HCl和NO3-增加的比例較多,SO42-則較少.在時間上,HCl濃度增加的高值在Cl-高值之前,SO42-濃度增加量主要受海鹽氣溶膠的硫酸鹽組分影響,但貢獻較少,NO3-濃度變化主要受氯虧損機制和NO3-干沉降過程影響,且其濃度增加的高值在Cl-高值前后均有出現.

[1] Brook J R, Makar P A, Sills D M L, et al. Exploring the nature of air quality over southwestern ontario: Main findings from the border air quality and meteorology study [J]. Atmospheric Chemistry and Physics, 2013,13(20):10461-10482.

[2] Jeffreys H. On the dynamics of wind [J]. Quarterly Journal of the Royal Meteorological Society, 1922,48(201):29-48.

[3] Haurwitz B. Comments on the sea-breeze circulation [J]. Journal of Meteorology, 1947,4(1):1-8.

[4] 孫 貞,高榮珍,張 進,等.青島地區8月一次海風環流實例分析和WRF模擬[J]. 氣象, 2009.35(8):76-84. Sun Z, Gao R Z, Zhang J, et al. An analysis and WRF simulation of a case of sea-breeze circulation in August at Qingdao [J]. Meteorological Monthly, 2009,35(8):76-84.

[5] Papanastasiou D K, Melas D, Lissaridis I. Study of wind field under sea breeze conditions; an application of WRF model [J]. Atmospheric Research, 2010,98(1):102-117.

[6] 柯史釗,黃 健.華南海陸風的數值模擬 [J]. 熱帶氣象學報, 1993, 9(2):169-176. Ke S Z, Huang J. Numerical simulation of sea-land breeze in South China [J]. Journal of Tropical Meteorology, 1993,9(2):169-176.

[7] 文偉俊,沈桐立,丁治英,等.城市化對廣西夏季海陸風影響的數值試驗[J]. 熱帶氣象學報, 2009,25(3):350-356. Wen W J, Shen T L, Ding Z Y, et al. Numerical experiment on the influence of urbanization on sea-land breeze in Guangxi in summer [J]. Journal of Tropical Meteorology, 2009,25(3):350-356.

[8] 邱曉暖,范紹佳.海陸風研究進展與我國沿海三地海陸風主要特征[J]. 氣象, 2013,39(2):186-193. Qiu X N, Fan S J. Progress of sea-land breeze study and the characteristics of sea-land breeze in three coastal areas in China [J]. Meteorological Monthly, 2013,39(2):186-193.

[9] 張立鳳,張 銘,林宏源.珠江口地區海陸風系的研究 [J]. 大氣科學, 1999,33(5):581-589. Zhang L F, Zhang M, Lin H Y. Study on the sea-land-breeze system in the Pearl River Estuary region [J]. Chinese Journal of Atmospheric Sciences, 1999,33(5):581-589.

[10] Ding A J, Wang T, Zhao M, et al. Simulation of sea-land breezes and a discussion of their implications on the transport of air pollution during a multi-day ozone episode in the Pearl River Delta of China [J]. Atmospheric Environment, 2004,38(39):6737-6750.

[11] Lu X, Chow K C, Yao T, et al. Seasonal variation of the land-sea breeze circulation in the Pearl River Delta region [J]. Journal of Geophysical Research, 2009,114, D17112, doi:10.1029/2009JD011764.

[12] 吳 蒙,羅 云,吳 兌,等.珠三角干季海陸風特征及其對空氣質量影響的觀測 [J]. 中國環境科學, 2016,36(11):3263-3272. Wu M, Luo Y, Wu D, et al. Observation on the characteristics of sea-land breezes and its influence to air quality over Pearl River Delta region during dry season [J]. China Environmental Science, 2016, 36(11):3263-3272.

[13] 莊延娟,王伯光,劉 燦.珠江口西岸冬季海陸風背景下羰基化合物的初步研究 [J]. 中國環境科學, 2011,31(4):568-575. Zhuang Y J, Wang B G, Liu C. Preliminary research on the carbonyl compounds under of sea-land breeze at the west coast of the Pearl River Estuary [J]. China Environmental Science, 2011,31(4):568-575.

[14] Roth B, Okada K. On the modification of sea-salt particles in the coastal atmosphere [J]. Atmospheric Environment, 1998,32:1555-1569.

[15] Weis D D, Ewing G E. The reaction of nitrogen dioxide with sea salt aerosol [J]. The Journal of Physical Chemistry A, 1999,103:4865-4873.

[16] Jaenicke R. Tropospheric aerosols. In: Aerosol Cloud Climate Interact [M]. Ed. San Diego: Academic Press, 1993:1-31.

[17] Houghton J T, Ding Y, Griggs D J, et al. IPCC. Climate Change: The Scientific Basis [R]. Ed. Cambridge: Cambridge University Press, 2001.

[18] 劉 建,吳 兌,范紹佳,等.臺風對沿海城市細粒子中海鹽氣溶膠的影響 [J]. 環境科學學報, 2017,37(9):3255-3261. Liu J, Wu D, Fan S J, et al. Impact of Typhoon heavy precipitation on sea-salt aerosol in coastal city [J]. Journal of Environmental Science, 2017,37(9):3255-3261.

[19] Im U. Impact of sea-salt emissions on the model performance and aerosol chemical composition and deposition in the east Mediterranean coastal regions [J]. Atmospheric Environment, 2013,75:329-340.

[20] 劉 倩,王體健,李 樹,等.海鹽氣溶膠影響酸堿氣體及無機鹽氣溶膠的敏感性試驗 [J]. 氣候與環境研究, 2008,13(5):598-607. Liu Q, Wang T J, Li S, et al. Sensitivity tests of the impacts of sea salt aerosol on the acid-base gases and the inorganic salts [J]. Climatic and Environmental Research, 2008,13(5):598-607.

[21] Liu Y M, Zhang S T, Fan Q, et al. Accessing the impact of sea-salt emissions on aerosol chemical formation and deposition over Pearl River Delta, China [J]. Aerosol and Air Quality Research, 2015,15: 2232–2245.

[22] Skamarock W C, Klemp J B, Dudhia J, et al. A description of the Advanced Research WRF version 3 [Z]. NCAR Technical Note: NCAR/TN-475+STR, 2008.

[23] Otte T L, Pleim J E. The meteorology-chemistry interface processor (mcip) for the cmaq modeling system: Updates through mcipv3.4.1 [J]. Geoscientific Model Development, 2010,3(1):243-256.

[24] Byun D, Schere K L. Review of the governing equations, computational algorithms, and other components of the Models- 3Community Multiscale Air Quality (CMAQ) modeling system [J]. Applied Mechanics Reviews, 2006,59(2):51-77.

[25] Gantt B, Kelly J T, Bash J O. Updating sea spray aerosol emissions in the Community Multiscale Air Quality (CMAQ) model version 5.0.2 [J]. Geoscientific Model Development, 2015,8:3733–3746.

[26] Gong S L. A parameterization of sea-salt aerosol source function for sub- and super-micron particles [J]. Global Biogeochemical Cycles, 2003,17(1097).

[27] Leeuw G D, Neele F P, Hill M, et al. Production of sea spray aerosol in the surf zone [J]. Journal of Geophysical Research-Atmosphere, 2000,105:29397–29409.

[28] 常 鳴,樊少芬,王雪梅.珠三角土地覆被資料優選及在WRF模式中的初步應用 [J]. 環境科學學報, 2014,34(8):1922-1933. Chang M, Fan S F, Wang X M. Impact of refined land-cover data on WRF performance over the Pearl River Delta region, China [J]. Acta Scientiae Circumstantiae, 2014,34(8):1922-1933.

[29] Li M, Zhang Q, Kurokawa J I, et al. MIX: a mosaic Asian anthropogenic emission inventory under the international collaboration framework of the MICS-Asia and HTAP [J]. Atmospheric Chemistry and Physics, 2017,17:935-963.

[30] Guenther A, Karl T, Harley P, et al. Estimates of global terrestrial isoprene emissions using MEGAN (model of emissions of gases and aerosols from nature) [J]. Atmospheric Chemistry and Physics, 2006,6:3181-3210.

[31] 趙文龍.廣州地區海陸風條件下臭氧污染的數值模擬研究 [D]. 廣州:暨南大學, 2017. Zhao W L. Numerical simulation of ozone pollution in Guangzhou Area under sea-land breeze condition [D]. Guangzhou: Jinan University, 2017.

[32] Pleim J E, Chang J S. A non local closure model in the convective boundary layer [J]. Atmospheric Environment, 1992,26A:965-981.

[33] Pleim J E. A combined local and nonlocal closure model for the atmospheric boundary layer. Part I: model description and testing [J]. Journal of Applied Meteorology and Climatology, 2007,46:1383-1395.

[34] Chang J S, Middleton P B, Stockwell W R, et al. Acidic deposition: State of science and technology. Report 4. The regional acid deposition model and engineering model. Final report [J]. Journal of Materials Science Letters, 1990,9(7):772-773.

[35] 姚小紅,黃美元,高會旺,等.海鹽氯循環對輸送污染物致酸的緩沖作用[J]. 大氣科學, 1998,22(6):913-918. Yao X H, Huang M Y, Gao H W, et al. The buffering effect of sea salt chlorine cycle on acid transport of pollutants [J]. Chinese Journal of Atmospheric Sciences, 1998,22(6):913-918.

[36] McInnes L M, Covert D S, Quinn P K, et al. Measurements of chloride depletion and sulfur enrichment in individual sea-salt particles colllected from the remote marine boundary layer [J]. Journal of Geophysical Research, 1994,99(D4):8257-8268.

The influence of sea salt aerosols on atmospheric environment under sea-land-breeze circulation: a numerical simulation.

SHEN Ao1, TIAN Chun-yan1, LIU Yi-ming1, SHEN Chong1, CHEN Xiao-yang1, CHEN Xun-lai2, WANG Ming-jie2, LAN Jing3, FAN Qi1*

(1.School of Atmospheric Sciences, Guangdong Province Key Laboratory for Climate Change and Natural Disaster Studies, Sun Yat-sen University, Guangzhou 510275, China；2.Shenzhen National Climate Observatory, Shenzhen 518040, China；3.Guangzhou Climate and Agrometeorology Center, Guangzhou 511430, China)., 2019,39(4)：1427~1435

The influence of sea salt aerosols on atmospheric environment in a typical sea-land-breeze circulation over Pearl River Estuary (PRE) during 9~13September 2013 was investigated with the simulation based on the WRF-CMAQ modeling. The results indicated that the sea-breeze circulation was initiated at around 14:00 LST (local standard time) and subsequently reached its peak at around 17:00 LST, and formed a fan wind circulation along PRE. The mass concentration of Na+in the Pearl River Delta (PRD) increased from 17:00 LST and reached the maximum value at around 02:00 LST the next day, then gradually decreased after 10:00 LST. The substantial chloride depletion was observed for the sea salt aerosols with more depletion being seen during daytime than nighttime, and more depletion for fine particles than coarse particles. Sensitivity experiments were respectively conducted with and without the emissions of sea salt aerosols in the model system. The results showed that the mass concentrations of HCl, SO42-, and NO3-all increased with the emissions and more increases were seen for HCl and NO3-than for SO42-. The peak concentration of HCl arrived earlier than that of Cl-. The SO42-concentration was mainly affected by the sulfate component in the emitted sea salt aerosols, while the NO3-concentration was mainly affected by both the chloride depletion and its dry deposition.

Pearl River Estuary；sea-land-breeze circulation；sea salt aerosols；chloride depletion；numerical simulation

X131.1

A

1000-6923(2019)04-1427-09

2018-09-20

國家重點研發計劃大氣專項課題(2017YFC0210105,2016YFC0203305, 2016YFC0203600);深圳市科技計劃項目(JCYJ20170306150333250);國家自然科學基金資助項目(91544102,41630422);氣象行業專項項目(GYHY201406031)

*責任作者, 教授, eesfq@mail.sysu.edu.cn

沈 傲(1996-),女,安徽合肥人,中山大學博士研究生,主要從事大氣環境數值模擬研究.發表論文1篇.