楊清健,趙天良*,鄭小波,史建武,張朝能,常嘉成,張 凱,鐘曜謙,于 超

亞洲季風(fēng)強(qiáng)弱年蒙自市大氣環(huán)境容量差異估算

楊清健1,趙天良1*,鄭小波2,史建武3,張朝能3,常嘉成1,張 凱1,鐘曜謙3,于 超1

(1.南京信息工程大學(xué),中國(guó)氣象局氣溶膠-云-降水重點(diǎn)開放實(shí)驗(yàn)室,氣候與環(huán)境變化國(guó)際合作聯(lián)合實(shí)驗(yàn)室,江蘇 南京 210044；2.貴州省山地環(huán)境氣候研究所,貴州 貴陽 550002；3.昆明理工大學(xué),環(huán)境科學(xué)與工程學(xué)院,云南 昆明 650500)

利用區(qū)域空氣質(zhì)量模式WRF-Chem,對(duì)亞洲季風(fēng)氣候變化背景下云南省蒙自市大氣環(huán)境容量進(jìn)行模擬評(píng)估.根據(jù)標(biāo)準(zhǔn)化南亞夏季風(fēng)指數(shù)分別選取2005年和2015年為強(qiáng)、弱季風(fēng)年.對(duì)2015年四季 (以1月、4月、7月和10月為代表月)和2005年夏季(7月為代表月)的主要大氣污染物濃度進(jìn)行模擬.結(jié)果表明蒙自市2015年全年CO、NO2、SO2、PM2.5、PM10的大氣環(huán)境容量分別為120.31、1.127、1.875、1.267、1.688(×104t/a),其中各污染物冬季大氣環(huán)境容量最小,春季的最大(PM10除外),且PM2.5在冬季排放量已飽和.強(qiáng)季風(fēng)年相對(duì)弱季風(fēng)年夏季CO、NO2、SO2、PM2.5、PM10的大氣環(huán)境容量分別提升4.81%、3.86%、12.6%、18.4%、8.7%,其中 PM2.5的容量提升最高.亞洲季風(fēng)年際變化對(duì)云南高原空氣質(zhì)量及大氣環(huán)境容量具有重要的調(diào)制作用.

大氣環(huán)境容量；WRF-Chem；云南高原；季風(fēng)氣候變化

隨著近年來經(jīng)濟(jì)的發(fā)展,云貴高原(以下簡(jiǎn)稱高原)地區(qū)空氣質(zhì)量開始出現(xiàn)較大變化,近30多年來,整個(gè)高原的平均干能見度均逐漸降低,四季的平均霾日數(shù)均逐漸升高[1];鄭小波等[2]研究表明,從20世紀(jì)60年代到21世紀(jì)初,高原地區(qū)的年平均干能見度從34km降低到了27km,年平均干消光系數(shù)從20世紀(jì)80年代的0.176km-1升至2005年的0.190km-1;氣溶膠的增多導(dǎo)致日照時(shí)長(zhǎng)和地表溫度改變[3];且在蒙自區(qū)域,能見度的降低導(dǎo)致到達(dá)地表的總輻射下降[4].高原出現(xiàn)的空氣質(zhì)量不利變化,其原因和后果需要深入研究.

紅河州地處云南省的東南部,與越南接壤,是典型的季風(fēng)氣候區(qū),其夏季受南亞季風(fēng)影響較大[5].紅河州首府蒙自市位于滇南中心城市核心區(qū).近年來隨著紅河工業(yè)園和蒙自經(jīng)濟(jì)開發(fā)區(qū)的建立,本地排放增長(zhǎng)迅速,空氣質(zhì)量也有所下降.研究蒙自地區(qū)的大氣環(huán)境容量對(duì)探究云南邊陲地區(qū)空氣質(zhì)量現(xiàn)狀和潛在的大氣污染趨勢(shì)有重要意義.

大氣環(huán)境容量是指在滿足大氣環(huán)境目標(biāo)值的條件下,某區(qū)域大氣環(huán)境所能承納污染物的最大能力或所能允許排放的污染物的總量,是支撐空氣質(zhì)量管理決策和大氣污染物總量控制的重要參照.大氣環(huán)境容量主要取決于環(huán)境對(duì)污染物的自凈能力與自凈空間.若超過了容量的閾值,大氣環(huán)境就不能發(fā)揮其正常的功能,進(jìn)而使生態(tài)環(huán)境、人群健康及物質(zhì)財(cái)產(chǎn)受到損害[6].氣象條件是影響大氣環(huán)境容量最重要的因素之一,它影響著區(qū)域的大氣擴(kuò)散、稀釋能力,從而影響大氣環(huán)境容量[7].計(jì)算蒙自各個(gè)季節(jié)和不同大氣環(huán)流背景下的大氣環(huán)境容量,可為蒙自地區(qū)大氣環(huán)境質(zhì)量的管理控制提供科學(xué)依據(jù).

大氣環(huán)境容量的估算已經(jīng)有了不少研究和發(fā)展,目前主要有3種方法,分別為①A-P值法:如王涵瑾等[8]通過建立修正A值法得出成都市各污染物均為夏季環(huán)境容量最大,冬季的最小;②線性優(yōu)化法:如肖楊等[9]基于線性規(guī)劃模型,估算北京市通州區(qū)的二氧化硫最大允許排放量為41.311t/a;③模式模擬法:如郝吉明等[7]利用GEOS—Chem模擬計(jì)算出2013年京津冀和西北五省(自治區(qū))5種大氣污染物的排放量超出大氣環(huán)境容量.但是目前對(duì)于大氣環(huán)境容量的研究多集中于大氣環(huán)境的計(jì)算和計(jì)算方法、大氣環(huán)境容量與區(qū)域環(huán)境之間的關(guān)系等[10-12],缺少氣象條件變化特別是本世紀(jì)以來氣候變化導(dǎo)致季風(fēng)減弱的趨勢(shì)對(duì)大氣環(huán)境影響的研究.本文不僅意在厘清蒙自市大氣環(huán)境容量與空氣質(zhì)量變化的關(guān)系,而且研究不同季風(fēng)條件下主要污染物大氣環(huán)境容量的變化特點(diǎn).

由于A-P 值法和線性優(yōu)化法不能很好考慮如不同季風(fēng)條件對(duì)污染物的輸送與清除作用、大氣污染物的非均相化學(xué)轉(zhuǎn)化等因素,本研究采用以WRF-Chem[13]為基礎(chǔ)的模式模擬法計(jì)算大氣環(huán)境容量[14].該模式的模擬效果已經(jīng)得到廣泛驗(yàn)證.如周廣強(qiáng)等[15]利用WRF-Chem建立空氣質(zhì)量預(yù)報(bào)系統(tǒng),并有很好的預(yù)報(bào)效;Tie等[16]模擬上海市臭氧的變化情況,模擬結(jié)果與觀測(cè)有很高的一致性.

本文基于WRF-Chem,以污染物濃度達(dá)到二級(jí)年均(或日均)標(biāo)準(zhǔn)(GB3095-2012)[17]為環(huán)境目標(biāo),建立蒙自市主要大氣污染物環(huán)境容量迭代計(jì)算方法.首先以2015年蒙自市4個(gè)季節(jié)代表月(1月、4月、7月和10月)5種污染物(CO、NO2、SO2、PM2.5、PM10)達(dá)到二級(jí)年均(其中CO為日均)標(biāo)準(zhǔn)(GB3095-2012)為約束目標(biāo),計(jì)算了蒙自市的最大允許排放量(大氣環(huán)境容量);然后,以典型的強(qiáng)和弱夏季風(fēng)年(2005年和2015年)為例,對(duì)比分析蒙自市夏季大氣環(huán)境容量的變化特征,研究季風(fēng)變化可能對(duì)大氣環(huán)境容量產(chǎn)生哪些重要影響.

1 數(shù)據(jù)和方法

1.1 算法設(shè)計(jì)

本文用WRF-Chem模式估算大氣環(huán)境容量的技術(shù)路線如圖1所示:

圖1 大氣環(huán)境容量計(jì)算流程

1.1.1 基準(zhǔn)情形大氣組分濃度模擬 利用WRF- Chem模擬出2005年夏季和2015年四季蒙自市5種污染物的平均濃度,并將2015年四季5種主要污染物模擬濃度日均值與蒙自市環(huán)境部門3個(gè)國(guó)控站點(diǎn)(以下簡(jiǎn)稱國(guó)控站)5種主要污染物監(jiān)測(cè)數(shù)據(jù)日均值做對(duì)比,驗(yàn)證模式的合理性.

1.1.2 目標(biāo)值 以《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB3095- 2012)規(guī)定二級(jí)標(biāo)準(zhǔn)下CO的日平均濃度值、NO2、SO2、PM2.5、PM10的年平均濃度值4mg/cm3、40μg/cm3、60μg/cm3、35μg/cm3、70μg/cm3為目標(biāo)限值.計(jì)算各代表月份各種污染物目標(biāo)限值濃度與模式模擬濃度(提取3個(gè)國(guó)控站所在經(jīng)緯度的模擬數(shù)據(jù)并取平均)的比值.

1.1.3 調(diào)整排放源清單 基于污染物濃度與排放量呈線性關(guān)系的假定,將目標(biāo)區(qū)域排放源數(shù)據(jù)改為原來的倍,并對(duì)排放源清單進(jìn)行調(diào)整.

1.1.4 大氣環(huán)境容量估算 利用新的排放清單,再次模擬計(jì)算各季節(jié)各污染物平均濃度,重復(fù)(2),(3)步驟,直至各污染物均接近達(dá)標(biāo)限制,由此估算二級(jí)年均(或日均)標(biāo)準(zhǔn)下的大氣環(huán)境容量.

在迭代過程中,氣態(tài)污染物迭代首先達(dá)到最大允許排放量,由此可基本確定由氣態(tài)污染物生成的二次PM2.5(硫酸鹽、硝酸鹽等)的“最大允許排放量”.然后維持氣態(tài)污染物排放量不變,繼續(xù)迭代計(jì)算一次PM2.5和一次PM10的排放量,最終確定一次PM2.5和一次PM10的最大允許排放量.

1.2 模式設(shè)置

利用WRF-Chem模式的3.8.1版本,采用3層嵌套網(wǎng)格,如圖2所示,d01層(最外層)覆蓋了中國(guó)大部分地區(qū)及部分中亞、南亞地區(qū),d02層覆蓋了云南省及其周邊地區(qū),d03層覆蓋了紅河州及部分周邊地區(qū).最外層網(wǎng)格中心位于23.71°N、103.25°E, 3層網(wǎng)格的格點(diǎn)數(shù)分別是99×99、106×109、88×91,對(duì)應(yīng)的分辨率分別是27km、9km和3km,積分時(shí)間步長(zhǎng)為135s.模式垂直方向上采用地形跟隨坐標(biāo)系,分為27層.模擬層頂氣壓為100hPa.采用物理和化學(xué)過程參數(shù)化方案設(shè)置如表1所示.氣象初始場(chǎng)為FNL再分析資料.排放源為清華大學(xué)MIX 2010 年排放清單,該清單提供了包括SO2, NO2,CO, NH3, NMVOC(揮發(fā)性有機(jī)物), PM10, PM2.5, BC, OC, CO2等10種主要大氣化學(xué)成分,以及CB05和SAPRC-99兩種大氣化學(xué)機(jī)制的分組分NMVOC排放數(shù)據(jù)[18].模擬時(shí)段為2005年7月和2015年1、4、7、10共5個(gè)月份,模式輸出時(shí)間間隔為1h.模式氣象場(chǎng)spin-up時(shí)間為18h,化學(xué)場(chǎng)spin-up時(shí)間為3d.

圖2 WRF-Chem模擬區(qū)域3層網(wǎng)格設(shè)置和地形海拔高度(m)

表1 WRF-Chem模擬的物理化學(xué)參數(shù)化方案選擇

2 結(jié)果與分析

2.1 模擬驗(yàn)證

將模擬的四季各主要污染物模擬濃度日均值與國(guó)控站各主要污染物監(jiān)測(cè)數(shù)據(jù)日均值做對(duì)比,評(píng)估化學(xué)場(chǎng)模擬效果;并將模擬的地表風(fēng)速和地面2m溫度日均值與中國(guó)氣象局MICAPS(氣象信息綜合分析處理系統(tǒng))中蒙自站觀測(cè)資料日均值對(duì)比,評(píng)估氣象場(chǎng)模擬效果.其中為相關(guān)系數(shù),為顯著性水平.RMSE為均方根誤差.IOA為一致性指數(shù),計(jì)算式見式(1),IOA的變化范圍為0~1之間,且IOA指數(shù)越接近于1代表模擬效果越好.MAD為平均絕對(duì)偏差.

利用MICAPS蒙自站觀測(cè)資料(站號(hào)56985)分別對(duì)模擬的強(qiáng)和弱季風(fēng)年夏季的地面溫度場(chǎng)和風(fēng)場(chǎng)進(jìn)行驗(yàn)證.結(jié)果如圖3和表2所示.

(a)弱季風(fēng)年夏季溫度場(chǎng),(b)弱季風(fēng)年夏季風(fēng)場(chǎng),(c)強(qiáng)季風(fēng)年夏季溫度場(chǎng),(d)強(qiáng)季風(fēng)年夏季風(fēng)場(chǎng)

從圖3可知,不論是強(qiáng)還是弱季風(fēng)年,溫度場(chǎng)模擬值和觀測(cè)值變化趨勢(shì)的一致性較好,相關(guān)系數(shù)均通過99%的顯著性檢驗(yàn).模擬風(fēng)速變化趨勢(shì)與觀測(cè)基本一致,但在一些時(shí)段出現(xiàn)模擬較觀測(cè)值略微偏低現(xiàn)象.分析誤差原因(1)可能是地處高原山地的蒙自市周邊地形復(fù)雜,氣象站所在盆地與城市周邊山地的高差較大,模式難以精確模擬地表風(fēng)速在此下墊面上隨時(shí)間的變化;(2)模式氣象初始場(chǎng)為再分析資料驅(qū)動(dòng),與觀測(cè)資料有一定差異,且再分析資料分辨率為1°×1°,難以提供精準(zhǔn)的初始邊界條件;(3)模式采用的參數(shù)化方案并不能準(zhǔn)確反映某一物理過程的發(fā)生機(jī)制,也不一定適用與各種復(fù)雜地形和不同外部條件.總的來說,可以認(rèn)為氣象場(chǎng)模擬效果較好.從圖3(b)和(d)還可知,強(qiáng)季風(fēng)年的風(fēng)速明顯大于弱季風(fēng)年.

表2 WRF-Chem氣象場(chǎng)模擬評(píng)估統(tǒng)計(jì)參量

注:V05:強(qiáng)季風(fēng)年地表風(fēng)速; V15:弱季風(fēng)年地表風(fēng)速; T05:強(qiáng)季風(fēng)年地面2m溫度; T15:弱季風(fēng)年地面2m溫度.

將模擬的2015年4季5種主要污染物日均值與3個(gè)國(guó)控站點(diǎn)同期監(jiān)測(cè)日均值的平均值進(jìn)行比較,評(píng)估其模擬效果(2005年缺少環(huán)境觀測(cè)數(shù)據(jù),以2015年數(shù)據(jù)進(jìn)行檢驗(yàn)).檢驗(yàn)圖如圖4所示

表3 WRF- Chem化學(xué)場(chǎng)模擬效果統(tǒng)計(jì)參量

由圖4和表3可知,5種污染物的相關(guān)系數(shù)均通過了99%的顯著性檢驗(yàn).且PM2.5,PM10,CO, NO24種污染物90%以上的點(diǎn)均分布在2倍誤差線以內(nèi),同時(shí)MAD和RMSE均較低,IOA除了SO2都在0.5以上;出現(xiàn)部分模擬偏差主要原因是:(1)排放源的制作年份為2010年,未及時(shí)更新,且在云南偏遠(yuǎn)地區(qū)排放源存在誤差,而污染物濃度對(duì)排放源網(wǎng)格數(shù)據(jù)的敏感性非常強(qiáng),提升排放源數(shù)據(jù)的準(zhǔn)確度能提高模擬精度[19];(2)模式對(duì)污染物的一些反應(yīng)機(jī)理尤其是二次有機(jī)氣溶膠的模擬機(jī)理還不是很完善,對(duì)二次PM2.5的模擬和一些顆粒物表面的異相模擬未進(jìn)行考慮[20].

綜合氣象場(chǎng)驗(yàn)證、模擬污染物濃度與觀測(cè)值驗(yàn)證,并考慮到蒙自地區(qū)地形和排放源的影響,可以認(rèn)為WRF-Chem模式能夠合理模擬蒙自市主要大氣污染物,并適用于蒙自市大氣環(huán)境容量的估算和對(duì)比分析.

2.2 蒙自市大氣環(huán)境容量估算

各季節(jié),各污染物在蒙自市的大氣環(huán)境容量如圖5所示,其中已知容量為MIX排放源中整個(gè)蒙自市所在格點(diǎn)的排放量之和,剩余容量為大氣環(huán)境容量減去已知容量.全市的全年CO、NO2、SO2、PM2.5和PM10的大氣環(huán)境容量分別為120.31′104t/a、1.127′104t/a、1.875′104t/a、1.267′104t/a和1.688′104t/a.

由于不同季節(jié)大氣環(huán)境條件不同,造成了四季有不同的大氣環(huán)境容量.對(duì)于PM2.5,大氣環(huán)境容量在冬季最小,秋季次之,春季和夏季相當(dāng).PM2.5剩余容量在冬季為0,即冬季排放已飽和;在夏季最大,其次是春季和秋季.對(duì)于PM10,大氣環(huán)境容量從小到大按季節(jié)排列為冬季<秋季<春季<夏季,且剩余容量有相同趨勢(shì).對(duì)于CO,其大氣環(huán)境容量表現(xiàn)為冬季<秋季<夏季<春季;而NO2和SO22種污染物則表現(xiàn)為冬季<夏季<秋季<春季;CO、NO2、SO23種污染物的剩余容量均與大氣環(huán)境容量呈相同趨勢(shì).

從全年平均剩余容量看,CO剩余容量最大,其所占比例達(dá)到了94.1%,而PM2.5剩余容量最小,僅為23.4%,PM10也僅有25.5%.

綜上所述,蒙自地區(qū)顆粒物排放如果不加控制,則今后成為環(huán)境污染主因的可能性較大.

鑒于此,可針對(duì)不同污染物、不同季節(jié)采取差異化管理.特別是冬季PM2.5排放現(xiàn)已經(jīng)飽和的狀況下,建議在工廠較多區(qū)域,對(duì)一月份的顆粒物排放采取管控措施,同時(shí)可將排放較多的生產(chǎn)活動(dòng)盡量安排在4~7月,以充分利用大氣自凈能力.

2.3 強(qiáng)弱夏季風(fēng)年蒙自市夏季大氣環(huán)境容量分析

地處低緯度的蒙自市受季風(fēng)影響較大,形成冬干夏雨的干濕分明季風(fēng)氣候[21].由于蒙自在夏季主要受到南亞季風(fēng)影響[5],根據(jù)本世紀(jì)以來標(biāo)準(zhǔn)化的南亞夏季風(fēng)指數(shù)[22],如圖6所示, 2005年是一個(gè)相對(duì)強(qiáng)的夏季風(fēng)年,2015年則是絕對(duì)弱的夏季風(fēng)年,結(jié)合考慮蒙自的地形、排放源等具有時(shí)空變化特點(diǎn)的因素,選擇相距10a的2005年和2015年為強(qiáng)弱夏季風(fēng)代表年,并取相同的排放源,對(duì)其夏季(7月)的大氣環(huán)境容量進(jìn)行模擬估算.

由表4可知,強(qiáng)夏季風(fēng)年P(guān)M2.5、PM10、NO2、SO2、CO 5種污染物的大氣環(huán)境容量相對(duì)于弱夏季風(fēng)年分別有18.4%、8.7%、3.86%、12.6%、4.81%的提升;其中,PM2.5的提升最大,SO2和PM10次之,強(qiáng)季風(fēng)對(duì)細(xì)顆粒物的擴(kuò)散作用對(duì)于當(dāng)?shù)叵募綪M2.5大氣環(huán)境容量的提升影響最為顯著;而對(duì)NO2的提升最小,說明在排放量等其他影響條件一致時(shí),強(qiáng)弱季風(fēng)對(duì)NO2的污染物濃度影響程度較小.

在不斷調(diào)整增減排放源清單,迭代計(jì)算的過程中,可得到污染物濃度隨排放量變化的關(guān)系,如圖7所示.其中PM2.5,PM10,SO23種污染物在相同排放量強(qiáng)度下,強(qiáng)季風(fēng)年對(duì)應(yīng)的污染物濃度均比弱季風(fēng)年對(duì)應(yīng)的低,而CO和NO2變化不大.對(duì)于顆粒物,其總體變化率(提升相同的污染物濃度需要增加的污染物排放量)均是強(qiáng)季風(fēng)年大于弱季風(fēng)年,這反映了強(qiáng)季風(fēng)條件(一般表現(xiàn)為較大的風(fēng)速和更多的降水)對(duì)顆粒物的良好清除效應(yīng);且隨著污染物排放量的提升,顆粒物強(qiáng)弱季風(fēng)年的污染物濃度隨排放量增加的變化率均有所減小,說明當(dāng)污染物超過一定量時(shí),大氣的自凈能力也有所下降.

圖6 標(biāo)準(zhǔn)化南亞夏季風(fēng)指數(shù)年際變化

表4 強(qiáng)弱夏季風(fēng)年夏季大氣環(huán)境容量比較

在不斷調(diào)整增減排放源清單,迭代計(jì)算的過程中,也可得到以二級(jí)年均濃度為達(dá)標(biāo)值(CO為二級(jí)日均濃度)的達(dá)標(biāo)率隨排放量變化的關(guān)系(達(dá)標(biāo)率為污染物日均濃度小于二級(jí)年均濃度標(biāo)準(zhǔn)值的天數(shù)與總天數(shù)的比值),如圖8所示.因?yàn)榈恍枭倭看螖?shù),所以該關(guān)系圖僅以離散化數(shù)據(jù)作圖,可能會(huì)存在一些誤差,但仍能看出明顯變化趨勢(shì)和差異.

對(duì)于PM2.5,取相同的達(dá)標(biāo)率(固定縱坐標(biāo)),可以看到,強(qiáng)季風(fēng)年比弱季風(fēng)年對(duì)應(yīng)的排放量(橫坐標(biāo))更大,且當(dāng)排放量增大到10000t/a以上時(shí),這種差異趨勢(shì)變得越來越大.SO2、CO和NO2也有此差異化特點(diǎn),但沒有PM2.5那么顯著.對(duì)于PM10,總體從斜率上看,相同達(dá)標(biāo)率情況下,強(qiáng)季風(fēng)年可排放更多;但中間有一部分在相同達(dá)標(biāo)率情況下,強(qiáng)季風(fēng)年可排放反而低于弱季風(fēng)年,這可能與強(qiáng)季風(fēng)將外沿地區(qū)PM10輸送至本地區(qū),導(dǎo)致蒙自PM10濃度上升有關(guān)[23];但隨著排放量的持續(xù)加大,強(qiáng)季風(fēng)年的清除效應(yīng)更強(qiáng),所以弱季風(fēng)年的PM10達(dá)標(biāo)率急劇下降,而強(qiáng)季風(fēng)年則緩慢下降.

3 討論

基于目前大氣物理化學(xué)過程機(jī)制的認(rèn)識(shí),WRF- Chem主要考慮了污染物排放、傳輸、化學(xué)轉(zhuǎn)化、干濕沉降等過程,能反映出大氣污染物濃度與一次排放源和其前體物排放量之間的非線性關(guān)系.但考慮到實(shí)際計(jì)算大氣環(huán)境容量的可操作性,本文對(duì)排放量和模擬污染物濃度之間的復(fù)雜非線性關(guān)系作了線性近似.同時(shí),本文在迭代計(jì)算的過程中,首先使氣態(tài)污染物迭代達(dá)到最大允許排放量,由此可基本確定由氣態(tài)污染物生成的二次PM2.5(硫酸鹽、硝酸鹽等)的最大允許排放量.然后維持氣態(tài)污染物排放量不變,繼續(xù)迭代計(jì)算一次PM2.5和一次PM10的排放量,并最終確定一次PM2.5和一次PM10的最大允許排放量.盡管這種方法有一定的局限性,但是一些研究者已經(jīng)將該方法用于大氣環(huán)境容量的計(jì)算[11,14,24]后續(xù)工作也將進(jìn)一步研究并完善計(jì)算方法,充分考慮二次污染物對(duì)大氣環(huán)境容量的影響.

除了氣象條件外,大氣污染源區(qū)域空間分布及其排放強(qiáng)度變化等也是影響大氣環(huán)境容量的重要因素.本文主要針對(duì)局地大氣污染物控制探究氣象條件對(duì)大氣環(huán)境容量影響程度,故所有模擬均采用同樣的排放源空間分布.后續(xù)工作也會(huì)通過改變排放源空間分布的敏感性試驗(yàn)來具體探究不同排放源分布對(duì)大氣環(huán)境容量的影響.

此外,目前我國(guó)污染物排放源清單在云貴高原地區(qū)缺少完整實(shí)測(cè)數(shù)據(jù),因此,與中東部地區(qū)相比具有更大的不確定性;同時(shí),目前有關(guān)污染物的化學(xué)反應(yīng)機(jī)制仍不完善.模式并不能精確捕捉完整的反應(yīng)機(jī)制.所以,若想提高大氣環(huán)境容量的計(jì)算精度,就必須提升排放源的精度和完善污染物生成相關(guān)的物理化學(xué)反應(yīng)機(jī)制.

4 結(jié)論

4.1 利用在線耦合氣象與化學(xué)模式的區(qū)域空氣質(zhì)量模式WRF-Chem,對(duì)云南省蒙自市進(jìn)行大氣環(huán)境容量精細(xì)估算,得出了2015年蒙自市CO、NO2、SO2、PM2.5和PM10的大氣環(huán)境容量為120.31,1.127, 1.875,1.267,1.688(×104t/a),其中各污染物冬季大氣環(huán)境容量最小,除PM10外的污染物春季大氣環(huán)境容量最大.

4.2 根據(jù)標(biāo)準(zhǔn)化南亞夏季風(fēng)指數(shù)選取2005年為強(qiáng)季風(fēng)年,2015年為弱季風(fēng)年,計(jì)算得到強(qiáng)季風(fēng)年夏季相對(duì)弱季風(fēng)年夏季CO、NO2、SO2、PM2.5和PM10的大氣環(huán)境容量分別提升4.81%、3.86%、12.6%、18.4%和8.7%,其中,PM2.5的提升最高,并得到不同季風(fēng)條件下污染物濃度隨排放量變化的關(guān)系,及達(dá)標(biāo)率隨排放量變化的關(guān)系.分析結(jié)果可為空氣質(zhì)量預(yù)報(bào)、污染控制的空間區(qū)域差異化管理、不同季風(fēng)強(qiáng)弱導(dǎo)致的污染物變化分析、不同季節(jié)污染物變化分析等提供參考.

4.3 利用在線耦合氣象與化學(xué)模式的區(qū)域空氣質(zhì)量模式WRF-Chem能較好模擬各個(gè)空間尺度上大氣物理和化學(xué)過程及不同氣象條件對(duì)污染物排放、擴(kuò)散、傳輸?shù)挠绊?可反映不同氣象條件下污染物排放量的差異與具有時(shí)空變化特征的大氣環(huán)境容量情況.

[1] 楊樂心.云貴高原區(qū)域大氣環(huán)境及氣候資源時(shí)空特征分析 [D]. 南京:南京信息工程大學(xué), 2015. Yang L X. Tempor-spatial characterization of atmospheric environment and climatic resources over the Yunnan-Guizhou Plateau [D]. Nanjing: Nanjing University of Information Science & Technology, 2015.

[2] 鄭小波,王學(xué)鋒,羅宇翔,等.云貴高原1961~2006年大氣能見度和消光因素變化趨勢(shì)及原因[J]. 生態(tài)環(huán)境學(xué)報(bào), 2010,19(2):314-319. Zheng X B, Wang X F, Luo Y X, et al. Long-term trends in visibility and atmospheric extinction coefficient over Yunnan-Guizhou Plateau in southwest China for 1961~2006 [J]. Ecology and Environmental Sciences, 2010,19(2):314-319.

[3] Zheng X B, Kang W M, Zhao T L, et al. Long-term trends in sunshine duration over Yunnan-Guizhou Plateau in Southwest China for 1961~ 2005 [J]. Geophysics Research Letter, 2008,35,L15707,doi:10.1029/ 2008GL034482.

[4] 鄭小波,王學(xué)鋒,羅宇翔,等.1961~2005年云貴高原太陽輻射變化特征及其影響因子 [J]. 氣候與環(huán)境研究, 2011,16(5):657-664. Zheng X B, Wang X F, Luo Y X, et al. Variation characteristics of global radiation and the associated climatic factors over Yunnan- Guizhou Plateau in the southwestern part of China for 1961~2005 [J]. Climatic and Environmental Research, 2011,16(5):657-664.

[5] 晏紅明,杞明輝,肖子牛.云南5月雨量與熱帶海溫異常及亞洲季風(fēng)變化的關(guān)系 [J]. 應(yīng)用氣象學(xué)報(bào), 2001,12(3):368-376. Yan H M, Qi M H, Xiao Z N.Relationship between precipitation in May over Yunnan and the changes of tropical sea surface temperature and Asia monsoon [J]. Quarterly Journal of applied meteorology, 2001, 12(3):368-376.

[6] 胡 毅,李 萍,楊建功,等.應(yīng)用氣象學(xué) [M]. 北京:氣象出版社, 2007. Hu Y, Li P, Yang J G, et al. Applied meteorology [M]. Beijing: China Meteorological Press, 2017.

[7] 郝吉明,許嘉鈺,吳 劍,等.我國(guó)京津冀和西北五省(自治區(qū))大氣環(huán)境容量研究 [J]. 中國(guó)工程科學(xué), 2017,19(4):13-19. Hao J M, Xu J Y, Wu J, et al. A Study of the atmospheric environmental capacity of Jingjinji and of the five northwestern provinces and autonomous regions in China [J].Engineering Sciences, 2017,19(4):13-19.

[8] 王涵瑾,王源程,倪長(zhǎng)健.基于修正A值法核算成都市季節(jié)大氣環(huán)境容量 [J]. 環(huán)境與可持續(xù)發(fā)展, 2015,40(3):71-74. Wang H J, Wang Y C, Ni C J.Calculating the seasonal atmospheric environmental capacity in Chengdu based on modified A value method [J]. Environment and Sustainable Development, 2015,40(3): 71-74.

[9] 肖 楊,毛顯強(qiáng),馬根慧,等.基于ADMS和線性規(guī)劃的區(qū)域大氣環(huán)境容量測(cè)算 [J]. 環(huán)境科學(xué)研究, 2008,21(3):13-16. Xiao Y, Mao X Q, Ma G H, et al. Atmospheric environmental capacity study based on ADMS Model and linear programming [J]. Research of Environmental Sciences, 2008,21(3):13-16.

[10] 李海晶.大氣環(huán)境容量估算及總量控制方法的研究進(jìn)展 [J]. 四川環(huán)境, 2007,(1):67-71. Li H J. Research progress on the methods of the atmospheric environmental capacity calculation and total quantity control [J]. Sichuan Environment, 2007,(1):67-71.

[11] 薛文博,付 飛,王金南,等.基于全國(guó)城市PM2.5達(dá)標(biāo)約束的大氣環(huán)境容量模擬 [J]. 中國(guó)環(huán)境科學(xué), 2014,34(10):2490-2496. Xue W B, Fu F, Wang J N, et al. Modeling study on atmospheric environmental capacity of major pollutants constrained by PM2.5compliance of Chinese cities [J]. China Environmental Science, 2014, 34(10):2490-2496.

[12] 錢躍東,王勤耕.針對(duì)大尺度區(qū)域的大氣環(huán)境容量綜合估算方法 [J]. 中國(guó)環(huán)境科學(xué), 2011,31(3):504-509. Qian Y D, Wang Q G.An integrated method of atmospheric environmental capacity estimation for large-scale region [J]. China Environmental Science, 2011,31(3):504-509.

[13] Grell G A, Schmitz P R, Mckeen S A, et al. Fully coupled 'online' chemistry within the WRF model [J]. Atmospheric Environment, 2005, 39(37):6957-6975.

[14] 常嘉成,趙天良,譚成好,等.基于WRF-Chem模擬的玉溪市大氣環(huán)境容量精細(xì)估算 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2017,37(10):3876-3884. Chang J C, Zhao T L, Tan C H, et al. An elaborative assessment of atmospheric environmental capacity in Yuxi based on WRF-Chem modeling [J]. Acta Scientiae Circumstantiae, 2017,37(10):3876-3884.

[15] 周廣強(qiáng),謝 英,吳劍斌,等.基于WRF-Chem模式的華東區(qū)域PM2.5預(yù)報(bào)及偏差原因 [J]. 中國(guó)環(huán)境科學(xué), 2016,36(8):2251-2259. Zhou G Q, Xie Y, Wu J B, et al. WRF-Chem based PM2.5forecast and bias analysis over the East China Region [J]. China Environment Science, 2016,36(8):2251-2259.

[16] Tie X, Geng F, Peng L, et al. Measurement and modeling of O3variability in Shanghai, China: Application of the WRF-Chem model [J]. Atmospheric Environment, 2009,43(28):4289-4302.

[17] GB 3095-2012 環(huán)境空氣質(zhì)量標(biāo)準(zhǔn) [S]. GB 3095-2012 Ambient air quality standard [S].

[18] 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(2):935-963.

[19] Theobald M R, David S, Massimo V. Improving the spatial resolution of air-quality modelling at a European scale – development and evaluation of the Air Quality Re-gridder Model (AQR v1.1) [J]. Geoscientific Model Development, 2016,9(12):4475-4489.

[20] Zheng B, Zhang Q, Zhang Y, et al. Heterogeneous chemistry: a mechanism missing in current models to explain secondary inorganic aerosol formation during the January 2013 haze episode in North China [J]. Atmospheric Chemistry and Physics, 2015,15(4):2031- 2049.

[21] 陳皖滇.蒙自市環(huán)境空氣質(zhì)量變化特征[J]. 環(huán)境科學(xué)導(dǎo)刊, 2016, 35(S1):82-86. Chen W D. Changing characteristics of air quality in Mengzi City [J]. Environmental Science Survey, 2016,35(S1):82-86.

[22] Li J P, Zeng Q C. A new monsoon index and the geographical distribution of the global monsoons [J]. Advances in Atmospheric Sciences, 2003,20(2):299-302.

[23] 程敘耕,何金海,車慧正,等.1980~2010年中國(guó)區(qū)域地面風(fēng)速對(duì)能見度影響的地理分布特征 [J]. 中國(guó)沙漠, 2013,33(6):1832-1839. Cheng X G, He J H, Che H Z, et al. Impact of surface wind speed on atmospheric visibility and its geographic pattern over China in 1980~2010 [J]. Journal of Desert Research, 2013,33(6):1832-1839.

[24] 李 莉,程水源,陳東升,等.基于CMAQ的大氣環(huán)境容量計(jì)算方法及控制策略 [J]. 環(huán)境科學(xué)與技術(shù), 2010,33(8):162-166. Li L, Cheng S Y, Chen D S, et al. A Calculated methodology of atmospheric environmental capacity based on CMAQ and control strategy [J]. Environmental Science & Technology, 2010,33(8):162- 166.

致謝：感謝中國(guó)環(huán)境監(jiān)測(cè)總站提供的空氣質(zhì)量數(shù)據(jù).

Assessments of the differences of atmospheric environmental capacity between strong and weak Asian monsoon years in Mengzi.

YANG Qing-jian1, ZHAO Tian-liang1*, ZHENG Xiao-bo2, SHI Jian-wu3, ZHANG Chao-neng3, CHANG Jia-cheng1, ZHANG Kai1, ZHONG Yao-qian3, YU Chao1

(1.Joint International Research Laboratory of Climate and Environment Change, Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing 210044, China；2.Guizhou Institute of Mountainous Environment and Climate, Guiyang 550002, China；3.Faculty of Environment Science and Engineering, Kunming University of Science and Technology, Kunming 650500, China)., 2019,39(10)：4054~4064

Atmospheric environmental capacity is an important index to support air quality management and the total quantity control of air pollutant emissions. The air quality model WRF-Chem was used to simulate and estimate the atmospheric environmental capacity under the background of Asian monsoon climate change in Mengzi over the Yunnan Plateau. 2005 and 2015 respectively was chosen as the strong and weak monsoon years according to the normalized South-Asian Monsoon Index. Concentrations of major atmospheric pollutants at all seasons (represented by January, April, July and October respectively) of 2015 and summer (represented by July) of 2005 were simulated, the atmospheric environmental capacities of CO, NO2, SO2, PM2.5and PM10in Mengzi were estimated with 120.31、1.127、1.875、1.267、1.688(×104t/a), respectively. The atmospheric environmental capacity of major air pollutants in winter was the smallest; and in spring was the largest excepting PM10, and the PM2.5emissions in winter reached saturation. Compared with the 2015 summer with weak monsoon, the atmospheric environmental capacities of CO, NO2, SO2, PM2.5and PM10increased respectively by 4.81%, 3.86%, 12.6%, 18.4% and 8.7% in 2005 summer with strong monsoon, and the PM2.5increased most. The interannual variation of the Asian monsoon plays an important role in regulating the air quality and atmospheric environmental capacity of the Yunnan plateau.

atmospheric environmental capacity；WRF-Chem；Yunnan plateau；monsoon climate change

X51

A

1000-6923(2019)10-4054-11

楊清健(1994-),男,河南鄭州人,南京信息工程大學(xué)碩士研究生,主要研究大氣環(huán)境數(shù)值模擬及觀測(cè)分析.

2019-04-01

國(guó)家自然科學(xué)基金資助項(xiàng)目(91744209,21667014);江蘇省研究生科研與實(shí)踐創(chuàng)新計(jì)劃項(xiàng)目(KYCX18_1027)

* 責(zé)任作者, 教授, tlzhao@nuist.edu.cn