珠三角西部一次冬季臭氧污染過程分析

麥健華,于玲玲,鄧雪嬌,歐洪輝,余欣洋

珠三角西部一次冬季臭氧污染過程分析

麥健華1,于玲玲2,鄧雪嬌3*,歐洪輝4,余欣洋1

(1.中山市氣象局,廣東 中山 528400；2.廣東省氣象臺,廣東 廣州 510641；3.中國氣象局廣州熱帶海洋氣象研究所,廣東廣州 510641；4.廣東中山生態環境監測站,廣東 中山 528400)

利用地面氣象和空氣質量觀測資料、高空氣象觀測資料以及再分析數據,對2021年12月8～12日位于珠三角西部的中山市的O3污染過程的特點及成因進行分析.結果表明,該次過程為近年來罕見的冬季O3連續污染過程,期間O3平均峰值濃度是當月非污染日的2.3倍.污染期間受副熱帶高壓和地面大陸高壓脊控制,中山市氣溫、日照和相對濕度偏高,風速偏小,日間風向以北風為主,近地面的下沉氣流、逆溫層和日間相對濕度下降對O3污染起到增幅作用.污染過程存在明顯的外來污染物輸送現象,污染潛在源區主要位于中山市東北方向.在前體物積累階段,大量的NO2從上游地區流入中山市,為O3污染的出現提供了充足的前體物.污染期間NO2在有利的太陽輻射條件下,光解速率比當月非污染日平均偏高23%,導致日間O3濃度快速上升,并在較差的擴散條件下形成污染.

臭氧污染；氣象條件；輸送通量；光解速率

臭氧(O3)是大氣中的氮氧化物(NO)和可揮發性有機物(VOCs)在光化學反應作用下生成的二次空氣污染物[1],近地面高濃度的O3對氣候、植被及人體健康均會產生嚴重的影響[2-4].O3污染事件的發生是光化學反應前體物排放和氣象條件共同作用的結果[5-7],準確的排放源清單為大氣污染事件提供了基礎數據集[8-9],而氣象條件作為污染事件的外因,具有影響要素多、時間變化快的特點,一直是研究O3污染事件的關鍵[10-14].過往研究從宏觀上對容易導致O3污染的天氣類型進行了分類,雖然各個地區由于地理位置差別導致分型有所差異,但我國大部分O3污染天氣類型與副熱帶高壓、臺風外圍及地面高壓脊相關[15-17].氣象條件主導著O3的生消、傳輸和擴散[18],盡管天氣類型有差異,但絕大多數O3污染事件發生在氣溫高、濕度低、輻射強以及一定的風力等大氣條件下[19-22].隨著探測技術的進步,O3激光雷達等垂直探測儀器的出現使大氣環境探測從近地面延展到對流層內[23-24].研究表明,逆溫層、下沉運動等大氣層結分布特征會使O3在邊界層內積聚,對地面污染起到一定的增幅作用[25-26].

經過數年的有效整治,廣東省的顆粒物污染在近年間已經得到有效遏制[27],但隨之而來的O3污染卻呈現出上升趨勢[28],據統計,2006～2019年,珠江三角洲地區O3區域平均值從48μg/m3上升到60μg/m3[29], O3已經取代細顆粒物成為影響廣東省空氣質量最主要的大氣污染物[30].珠三角O3污染具有持續時間長、影響范圍大、污染程度重等特點,因此近年來該地區O3污染的研究已成為熱點,這些研究主要聚焦在O3的污染來源[31]、O3與前體物的關系[32-34]及污染的氣象成因等方面[35-37].O3污染的出現與地區氣候特征密切相關,廣東省的O3污染主要出現在夏、秋兩季[38-39],而冬季由于氣溫低、輻射弱,O3污染出現總體較少.但是2021年12月上中旬,珠江三角洲西部到廣東西部一帶出現了一次罕見的大范圍初冬O3污染過程,污染時段集中在12月8～12日,珠三角重要城市之一的中山市連續5d出現污染,是該次過程污染最為嚴重的城市.本文利用地面空氣質量和氣象觀測數據、探空氣象數據、歐洲中心(ECMWF)的ERA5及美國國家環境預報中心(NCEP)的GDAS再分析資料,對2021年冬季出現在中山市的此次罕見連續O3污染過程進行分析,研究該過程的污染特征及成因,以期為該地區O3污染的預報預警提供參考.

1 數據與方法

1.1 數據來源

本文使用的地面空氣質量觀測數據來源于廣東省生態環境監測中心,包括了2021年12月珠江口周邊城市71個空氣質量監測站逐小時的O3和NO2實況數據,站點分布如圖1所示.氣象數據來源于中山氣象站,站點位于中山市中心的紫馬嶺公園內,包括了2021年12月逐小時的氣溫、相對濕度、降水量、風向風速及太陽輻射數據.在探空氣象數據方面,本文使用了溫廓線雷達及風廓線雷達在垂直方向上的氣溫和垂直速度數據.其中溫廓線雷達數據由廣東中山生態環境監測站提供,雷達同樣位于紫馬嶺公園內,包括了從地面到高空1000m的氣溫數據,垂直分辨率為50m,時間分辨率為1h.垂直速度探空數據采用了位于中山市西邊的江門新會市的風廓線雷達數據,雷達位置位于中山氣象站正西方約35km處,包括了從近地面到高空3000m的數據,垂直分辨率60m,時間分辨率1h.氣象再分析資料采用了ECMWF的ERA5及NCEP的GDAS再分析資料,空間分辨率分別為0.25°×0.25°和1°×1°,時間分辨率分別為3h和6h,垂直方向分別為1000～200hPa共14層及1000～20hPa共23層.

圖1 空氣質量監測站點和中山氣象站位置

NO2是O3的重要前體物,其光解速率是分析大氣光化學污染的重要指標.本文使用了PFS-100光解光譜儀的NO2光解速率觀測資料,并分析其與本次O3污染過程的關系.光解光譜儀通過接收來自各個方向的太陽輻射獲得一定波長范圍內的光譜信息,將光譜信息轉化為光化通量,并計算出多種大氣物質(包括O1D、NO2、NO3、HONO、HCHO、H2O2)的光解速率.儀器觀測地點位于中山市金鐘水庫,與中山氣象站直線距離約4km.

1.2 研究方法

根據《環境空氣質量指數(AQI)技術規定(試行)》[40],O3污染等級根據某天的O3濃度8h滑動平均(O3-8h)劃分為6個級別,分別為優(0～100μg/m3)、良(101～160μg/m3)、輕度污染(161～215μg/m3)、中度污染(216～265μg/m3)、重度污染(266～800μg/m3)和嚴重污染(>800μg/m3),因此當O3-8h超過160μg/m3時,則認為出現了O3污染.中山市的O3濃度采用位于紫馬嶺公園內的空氣質量國控站觀測的O3濃度數據.

利用MeteoInfo軟件和TrajStat插件及NCEP的GDAS再分析數據,計算了中山市O3污染的潛在源區(PSCF).PSCF方法假設后向軌跡經過區域的排放源會對分析區域的污染造成影響,當PSCF值大時,代表該地排放源對分析區域污染的潛在貢獻大.本文把O3濃度污染閾值設置為160μg/m3,即超過閾值的O3濃度對應的軌跡為污染軌跡,進行O3污染軌跡的后向24h聚類分析,并在此基礎上計算PSCF值.

本文通過計算O3污染期間,污染物在中山市范圍內(113.0～113.7°E,22.2～22.8°N)的輸送通量情況,量化了O3和NO2跨區域輸送對中山市污染過程的影響.利用圖1所示的空氣質量監測站點數據,通過Cressman分析方法[41],選取適當的影響半徑,把站點數據插值成格點數據,再結合ERA5的地面10m風場數據,計算污染物的輸送通量.定義

南北向收支=南邊界輸送通量-北邊界輸送通量 (2)

東西向收支=西邊界輸送通量-東邊界輸送通量 (3)

總收支=南北向收支+東西向收支 (4)

式中:南北方向上北風為負值,南風為正值,東西方向上東風為負值,西風為正值,當收支為正時,表示污染物凈流入,收支為負時,表示污染物凈流出.

2 結果與討論

2.1 污染過程概況

2021年12月上、中旬,珠三角西部到粵西出現了一次罕見的大范圍冬季O3污染過程,污染時段集中在12月8～12日,圖2為污染過程期間廣東各地市O3污染天數,可見此次污染過程各城市污染天數從1～5d不等,其中位于珠三角西部的中山市連續5d出現污染,是此次過程污染最為嚴重的地區.

圖2 2021年12月8～12日廣東省O3污染日數

審圖號:GS(2019)1822號

圖3 中山市2013～2021年12月份O3-8h(mg/m3)和污染等級

圖3為2013～2021年間中山市12月份的O3污染日歷,可見12月份出現O3污染的天數較少,近9a僅出現9d,且多是以持續時間為1d的短期污染為主,僅2021年12月8～12日出現了連續5d的長期污染,因此本次污染過程在近年非常罕見.本次過程O3-8h在162～184μg/m3之間,屬輕度污染.圖4為2021年12月中山市O3污染日與非污染日中,各時次平均O3濃度的日變化分布,可見污染日和非污染日中O3濃度均體現出明顯的晝高夜低日變化特征,絕大部分時次污染日的濃度要高于非污染日,在中午到下午O3濃度快速上升并達到峰值時,污染日最高峰值濃度達269μg/m3,平均峰值濃度214μg/m3,為非污染日的2.3倍.傍晚O3濃度開始下降,但污染日濃度仍顯著高于非污染日濃度,部分時次污染日濃度可達非污染日的3倍以上,但在早晨8:00～9:00,污染日濃度與非污染日相當接近,甚至略低于非污染日濃度,可見即使是在連續污染期間,夜間到早晨O3的清除作用仍然十分明顯,出現污染的原因主要在于日間O3濃度的快速增長.

圖4 2021年12月份中山市O3污染日及非污染日平均濃度日變化

2.2 天氣形勢分析

圖5為此次連續O3污染過程的平均500hPa位勢高度和海平面氣壓及其距平,其中距平為相對于過去5a(2016～2020年)同時段的平均500hPa位勢高度和海平面氣壓而言.如圖所示,污染過程期間副熱帶高壓主體位于華南、南海及中南半島一帶,其中588dagpm線位于廣東省中南部,珠江口西岸到粵西沿海等污染主要區域位于副高的主體控制之下,從污染時段與過去5a的平均位勢高度距平來看,廣東地區為正距平,位勢高度偏高4～6dagpm,說明污染期間副高偏強,有利于晴好天氣的出現和光化學反應生成O3[42].從海平面氣壓的分布來看,北高南低的氣壓場分布相當明顯,且氣壓距平為正距平,說明大陸高壓脊較為強盛,地面維持長時間的盛行北風,有利于來自大陸方向的污染物向南輸送,使中山市等下游地區產生污染[43].因此,本次連續污染過程是在500hPa副高及地面冷高壓脊等有利的天氣系統控制之下所產生.

圖5 2021年12月8～12日500hPa平均位勢高度和平均海平面氣壓及距平

2.3 地面氣象要素特征及其與O3污染的關系

本次O3污染過程期間地面氣象要素表現出明顯的特征,圖6和表1分別為2021年12月上旬和中旬各天的O3-8h和氣象要素距平分布及相關系數.

圖6 2021年12月上中旬O3-8h及各氣象要素距平

可以看出污染期間地面氣象要素表現出明顯的氣溫、日照和相對濕度偏高,風速偏小的特征.12月8～12日污染期間5d的日最高氣溫距平均為正,距平在1.0～2.4℃之間.日照時數方面,除12月11日距平為負,其余4d的距平在2.5～2.9h之間.充足的日照以及較高的氣溫有利于光化學反應生成O3,導致污染的出現.8～12日的風速距平在-0.1～ -0.7m/s之間,靜小風不利于污染物的擴散.一般相對濕度較低時有利于O3污染的出現[44],但此次污染期間相對濕度偏高,距平在2%～15%之間,這是因為中山市冬天濕度偏低時,地面一般是受強冷高壓脊控制,此時氣溫低、風速大,并不利于O3污染的出現.從相關系數來看,O3-8h濃度與日最高氣溫及日平均風速的相關性最大,相關系數分別為0.61和-0.65,說明氣溫高、風速小是此次污染過程最直接的氣象影響因素.

表1 2021年12月上中旬O3-8h與各氣象要素距平的相關系數

注:表中相關系數均通過了顯著性水平為0.01的顯著性檢驗.

圖7 2021年12月上中旬日間南北風頻率、日最高氣溫及O3-8h

圖8 2021年12月8～12日中山市后向24h軌跡聚類及潛在源區

審圖號:GS(2019)1822號

除了風速外,本次O3污染過程的出現也跟地面風向密切相關.圖7為2021年12月上中旬日間(09:00～19:00)南、北風風向頻率、日最高氣溫及O3-8h.可以看出受大陸高壓脊影響,12月上中旬中山市大部分日期吹北風為主,僅12月15～16日2d南風頻率達到一半以上.12月8～12日污染期間日間的南風頻率在0～27%之間,且出現南風的時段集中在傍晚時分,對日間O3的快速增長沒有直接影響.中山地處珠三角冬季風的下風向位置,當吹北風時,加強了O3及其前體物向下游地區的輸送,容易使中山出現O3污染,而吹南風時則為清潔的海上氣流,一般不會造成污染.從圖8污染期間到達中山的后向24h軌跡聚類分析及潛在源區分布可知,主要的近地面氣流有兩支,且均為北風氣流,軌跡較長的一支氣流經粵東沿海和深圳到達中山,較短的一支氣流經東莞、廣州到達中山.從O3潛在源區來看,貢獻較大的源區均位于廣東省內,其中中山東北方向的東莞、深圳、惠州一帶潛在貢獻最大,PSCF值達到0.5以上,說明此次污染過程的污染物跨區域傳輸以短距離輸送為主.另外,從圖7可知,12月上中旬中12月16日的日最高氣溫最高,但16日南風頻率較大,且南風出現在中午前后,影響了O3的增長,因此沒有出現O3污染.

2.4 高空氣象要素特征分析

大氣垂直氣象條件對近地面空氣污染物的濃度變化具有重要影響.圖9是12月8～12日污染期間中山市附近垂直速度和相對濕度的垂直分布以及8:00的氣溫垂直廓線.如圖9a所示,由于污染期間受到副熱帶高壓及地面冷高壓脊的控制,中山市附近從地面到3000m高度為大片的下沉氣流區,垂直速度為正值,特別是1000m以下區域基本被下沉氣流控制,除了12月11日外,其余4d日出前在近地面均存在明顯的下沉速度大值區,部分時次和高度下沉速度可達0.6m/s.基于ERA5資料的相對濕度垂直分布顯示(圖9c),從1000hPa～800hPa相對濕度均存在明顯的日變化,夜間濕度較大,日出后濕度明顯下降,在975hPa以下的近地面層中,相對濕度從夜間的80%以上下降至日間的50%以下.日間的云量隨相對濕度的下降而減小,太陽輻射增強,有利于O3生成.在氣溫的垂直變化方面,由于日出后近地面濕度下降,晴空輻射出現,但此時太陽輻射仍較弱,地面熱量散失,有利于近地面逆溫層的形成[45].圖9b為污染期間每天早上08:00從地面到1000m高度的氣溫垂直廓線,可見雖然污染期間各天的氣溫差異明顯,但垂直分布趨勢較為一致,即從地面到400m高度氣溫變化不明顯,而從400～600m高度,5d均存在逆溫現象,兩層氣溫差在0.2～0.3℃之間,在600m高度以上,氣溫隨高度變化迅速下降.因此,本次污染過程大氣層結穩定,近地面存在下沉氣流及逆溫層,減弱了大氣湍流運動,空氣污染物難以向上擴散,同時日間濕度下降,輻射增強,有利于O3生成,加劇了地面污染.

圖9 2021年12月8～12日中山市附近垂直速度垂直分布、08:00氣溫垂直廓線及相對濕度垂直分布

2.5 污染期間O3與NO2的地面輸送通量特征

利用1.2節中污染物輸送通量的計算方法計算了O3污染期間地面O3和NO2的輸送通量特征如表2所示,根據珠三角附近O3及NO2的日變化特征[46],定義前體物積累階段為05:00～08:00,O3增長階段為11:00～17:00,O3下降階段為20:00～次日02:00.在前體物積累階段,由于上游地區的O3和NO2濃度較高,在北風背景下加強了對中山的跨區域傳輸,此時O3和NO2的總輸送通量均為正,說明O3和NO2均為凈流入,由于夜間NO通過化學反應不斷消耗O3[47],因此O3濃度維持在較低水平,但此時NO2的大量輸入為下一階段的光化學反應提供了充足的前體物.從方向上來看,這一階段的O3和NO2輸入以東邊界輸入為主,與圖8中O3污染潛在源區位于中山市東北面相對應.在O3增長階段,O3輸送通量在8日和11日為正,O3凈流入,但在其余3d為負,O3凈流出,即不存在明顯的流入和流出特征.該階段NO2的通量全為負,說明O3增長期間NO2為明顯的凈流出,且以西邊界流出為主,但由于在前體物累積階段中山市范圍內已經累積了大量的前體物,加上氣象條件有利于光化學反應,因此O3濃度快速增長.在O3下降階段,O3通量在12日為負,其余4d為正,說明此階段O3仍以凈流入為主,但由于光化學反應減弱,滴定作用增強,O3濃度快速下降,而這一階段NO2的東西向通量轉為正,但南北風通量有正有負,因此總通量的流入或流出特征不明顯.從以上討論可知,前體物積累階段中大量的NO2從上游地區流入中山市,并在白天有利的氣象條件下發生光化學反應,是此次O3連續污染過程出現的主要原因.

表2 2021年12月8～12日中山市范圍內O3與NO2的輸送通量

圖10 2021年12月8～12日08:00珠江三角洲地面NO2通量及NO2通量散度

圖中風矢為NO2通量(單位:μg/(m2·s)),填色為NO2通量散度(單位:μg/(m3·s)),其中負值表示輻合,正值表示輻散

圖10為2021年12月8～12日8:00珠江三角洲附近的地面NO2輸送通量及輸送通量散度.可以看出O3污染期間,在前體物積累階段地面的NO2輸送通量表現出相似的特征.在持續的東北風影響下,NO2自上游濃度較高的廣州、佛山、東莞等地區向下游地區輸送,因此珠三角北部和東部為NO2通量輻散區,而珠三角南部和西部的中山、江門、珠海等地則是NO2通量輻合區,其中中山市連續5d均存在明顯的NO2通量輻合,大量的NO2輸入為白天時段的光化學反應提供了充足的前體物.

2.6 NO2的光解速率特征

圖11 2021年12月份JNO2與太陽輻射強度及O3濃度變化的關系

NO2在太陽輻射下光解產生NO和O分子,這是對流層形成O3的化學反應的重要一步[48].圖11是2021年12月份NO2光解速率(NO2)與同一時次太陽輻射強度與O3濃度變化(即當前時次O3濃度與上一時次O3濃度的差)的關系,由于NO2在日出前及日落后的量級與日間差距大,因此僅對08:00～17:00數據進行分析.如圖11a所示,隨著太陽輻射強度的增強,NO2的增長相當明顯,二者接近線性關系,相關系數高達0.94.當NO2增強時,從圖10b可知O3濃度總體上呈增加趨勢,但NO2與O3濃度變化的離散度比太陽輻射強度要大,二者相關系數為0.53.從表3可知,當NO2在10-3量級以下時,此時多處于日出或日落前后時段,O3濃度為負增長或0增長,平均濃度變化為-6μg/m3;隨著NO2增大,O3濃度的增長趨勢相當明顯,當NO2從(1～4)′10-3s-1變化到(4～7)′10-3s-1,O3濃度變化的增幅達280%.當NO2繼續變大,O3的濃度變化明顯放緩,結合圖11b,說明NO2對O3濃度變化的影響并非線性,除NO2外,O3濃度變化還受到輸送和擴散等氣象因素的影響.從以上分析可知,當太陽輻射增強,NO2隨之增強,此時O3濃度總體上呈增長趨勢.

表3 不同JNO2范圍下O3濃度平均變化

圖12為2021年12月份O3污染日與非污染日的平均太陽輻射強度和NO2的日變化特征.可見無論是污染日還是非污染日,NO2與太陽輻射強度的日變化趨勢基本一致,在傍晚到第二天早晨接近0,從08:00開始迅速增加,午后達到最大,然后隨著太陽輻射的減弱NO2迅速下降,呈現明顯的單峰結構特征[49].污染日和非污染日的太陽輻射強度均在中午12:00達到峰值,污染日輻射強度相比非污染日偏高-5%～59%,平均偏高22%,在達到峰值的12:00污染日比非污染日偏高110W/m2.與太陽輻射類似,污染日的NO2相比非污染日偏高12%～31%,平均偏高23%,但NO2達到峰值的時間為13:00,相比太陽輻射晚1h,此時污染日NO2相比非污染日偏高22%.因此,12月8～12日在較為有利的太陽輻射條件下,NO2的光解速率上升,導致O3濃度上升,從而產生O3污染.但值得注意的是,較大的NO2并不一定導致O3污染,例如2021年12月中,平均NO2最大為12月1日的5.97×10-3s-1,而在連續污染過程中,最大的日均NO2為12月10日的5.24×10-3s-1,比12月1日偏小.因此,O3污染出現是前體物和氣象條件多種因素綜合作用的結果.

圖12 2021年12月份O3污染日與非污染日的平均太陽輻射強度和JNO2日變化

3 結論

3.1 2021年12月上、中旬珠江三角洲西部到廣東西部出現了一次冬季O3污染過程,中山市12月8～12日連續5d出現污染,為近年來罕見.污染期間日間O3濃度快速增長,平均峰值濃度達214μg/m3,是當月非污染日平均峰值濃度的2.3倍.

3.2 過程期間污染地區主要影響天氣系統為副熱帶高壓和地面大陸高壓脊,500hPa位勢高度及海平面氣壓與歷史同期相比為正距平.污染期間地面氣象要素特征為氣溫、日照和相對濕度偏高,風速偏小,日間風向以北風為主.垂直方向上1000m以下區域為下沉氣流區,大部分污染日在日出前近地面均存在0.6m/s的下沉速度大值區.近地面相對濕度日變化明顯,濕度從夜間的80%以上下降至日間的50%以下.早上08:00近地面層氣溫垂直變化小,從400m～600m高度存在逆溫現象,兩層氣溫差在0.2～0.3℃之間.有利的天氣形勢、地面及高空氣象條件是此次O3污染過程出現的主要原因.

3.3 此次污染過程存在明顯的外來污染物輸送現象,跨區域傳輸以短距離輸送為主,潛在源區主要位于中山市東北方向.在前體物積累階段,大量的NO2從上游地區凈流入中山市,為O3污染的出現提供了充足的光化學反應前體物.

3.4NO2日變化趨勢呈單峰分布且隨著太陽輻射的增強而增強,當NO2增強時,O3濃度總體上呈增長趨勢.污染日的平均太陽輻射強度相比當月非污染日平均偏高22%,峰值比非污染日偏高110W/m2.污染日的NO2相比非污染日平均偏高23%,峰值偏高22%.污染期間在有利的太陽輻射條件下,NO2的光解速率上升,導致O3濃度快速上升.

[1] QX/T 240-2014 光化學煙霧判識[S].QX/T 240-2014 Identification of photochemical smog [S].

[2] 趙 楠,盧毅敏.中國地表臭氧濃度估算及健康影響評估[J]. 環境科學, 2022,43(3):1235-1245.

Zhao N, Lu Y M. Estimation of surface ozone concentration and health impact assessment in China [J]. Environmental Science, 2022,43(3): 1235-1245.

[3] Pleijel H, Broberg M C, Uddling J. Ozone impact on wheat in Europe, Asia and North America-A comparison [J]. Science of The Total Environment, 2019,664:908-914.

[4] Pu X, Wang T J, Huang X, et al. Enhanced surface ozone during the heat wave of 2013 in Yangtze River Delta region, China [J]. Science of The Total Environment, 2017,603-604:807-816.

[5] 徐芯蓓,劉濤濤,徐玲玲,等.2017年廈門金磚會晤期間人為減排和氣象條件變化對臭氧污染特征的影響[J]. 環境科學學報, 2020, 40(12):4380-4389.

Xu X B, Liu T T, Xu L L, et al. Effects of emission control and changes in meteorological conditions on the characteristics of ozone pollution during the 2017 BRICS Summit in Xiamen [J]. Acta Scientiae Circumstantiae, 2020,40(12):4380-4389.

[6] Cheng N L, Chen Z Y, Sun F, et al. Ground ozone concentrations over Beijing from 2004 to 2015: Variation patterns, indicative precursors and effects of emission-reduction [J]. Environmental Pollution, 2018, 237:262-274.

[7] 黃小剛,邵天杰,趙景波,等.氣象因素和前體物對中國東部O3濃度分布的影響[J]. 中國環境科學, 2019,39(6):2273-2282.

Huang X G, Shao T J, Zhao J B, et al. Impact of meteorological factors and precursors on spatial distribution of ozone concentration in Eastern China [J]. China Environmental Science, 2019,39(6):2273- 2282.

[8] Ohara T, Akimoto H, Kurokawa J, et al. An Asian emission inventory of anthropogenic emission sources for the period 1980～2020 [J]. Atmospheric Chemistry and Physics, 2007,7:4419-4444.

[9] Li M, Zhang Q, Zheng B, et al. Persistent growth of anthropogenic non-methane volatile organic compound (NMVOC) emissions in China during 1990～2017: drivers, speciation and ozone formation potential [J]. Atmospheric Chemistry and Physics, 2019,19(13):8897- 8913.

[10] Deng T, Wang T J, Wang S Q, et al. Impact of typhoon periphery on high ozone and high aerosol pollution in the Pearl River Delta region [J]. Science of the Total Environment, 2019,668:617-630.

[11] 李莉莉,王 隆,劉喜平,等.哈爾濱市臭氧時空分布特征及氣象要素的關系[J]. 中國環境科學, 2020,40(5):1991-1999.

Li L L, Wang L, Liu X P, et al. Temporal and spatial distribution characteristics of ozone and its relationship with meteorological factors in Harbin [J]. China Environmental Science, 2020,40(5):1991- 1999.

[12] 王 玫,鄭有飛,柳艷菊,等.京津冀臭氧變化特征及與氣象要素的關系[J]. 中國環境科學, 2019,39(7):2689-2698.

Wang M, Zheng Y F, Liu Y J, et al. Characteristics of ozone and its relationship with meteorological factors in Beijing-Tianjin-Hebei region [J]. China Environmental Science, 2019,39(7):2689-2698.

[13] 劉 超,張恒德,張天航,等.青島“上合峰會”期間夜間臭氧增長成因分析[J]. 中國環境科學, 2020,40(8):3332-3341.

Liu C, Zhang H D, Zhang T H, et al. The causes of ozone concentration growth in the night during the “Shanghai Cooperation Organization Summit” in Qingdao [J]. China Environmental Science, 2020,40(8):3332-3341.

[14] Chen Z Y, Zhuang Y, Xie X M, et al. Understanding long-term variations of meteorological influences on ground ozone concentrations in Beijing during 2006～2016 [J]. Environmental Pollution, 2019,245: 29-37.

[15] 李婷苑,陳靖揚,翁佳烽,等.廣東省臭氧污染天氣型及其變化特征[J]. 中國環境科學, 2022,42(5):2015-2024.

Li T Y, Chen J Y, Weng J F, et al. Ozone pollution synoptic patterns and their variation characteristics in Guangdong Province [J]. China Environmental Science, 2022,42(5):2015-2024.

[16] 蒲 茜,李振亮,張 悅,等.重慶市O3污染日的大氣環流分型與傳輸特征[J]. 中國環境科學, 2021,41(1):18-27.

Pu X, Li Z L, Zhang Y, et al. Characterization of atmospheric circulation and transmission in Chongqing City during ozone polluted days [J]. China Environmental Science, 2021,41(1):18-27.

[17] 常爐予,許建明,瞿元昊,等.上海市臭氧污染的大氣環流客觀分型研究[J]. 環境科學學報, 2019,39(1):169-179.

Chang L Y, Xu J M, Qu Y H, et al. Study on objective synoptic classification on ozone pollution in Shanghai [J]. Acta Scientiae Circumstantiae, 2019,39(1):169-179.

[18] Tie X X, Geng F H, 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.

[19] 楊 健,尹沙沙,于世杰,等.安陽市近地面臭氧污染特征及氣象影響因素分析[J]. 環境科學, 2020,41(1):115-124.

Yang J, Yin S S, Yu S J, et al. Characteristic of surface ozone and meteorological parameters analysis in Anyang City [J]. Environmental Science, 2020,41(1):115-124.

[20] Alghamdi M A, Khoder M, Harrison R M, et al. Temporal variations of O3and NOin the urban background atmosphere of the coastal city Jeddah, Saudi Arabia [J]. Atmospheric Environment, 2014,94:205- 214.

[21] 陳志青,邵天杰,趙景波,等.東北地區臭氧濃度空間格局演變規律及其影響因素[J]. 環境科學學報, 2020,40(9):3071-3080.

Chen Z Q, Shao T J, Zhao J B, et al. Evolution and influencing factors of ozone concentration spatial distribution in Northeastern region [J]. Acta Scientiae Circumstantiae, 2020,40(9):3071-3080.

[22] Ejimofor C S, Okoro E C, Sivla W T. Effects of elevated humidity on stratospheric ozone content in the tropics [J]. International Journal of Physical Sciences, 2021,15(4):182-193.

[23] Klein A, Ancellet G, Ravetta F, et al. Characterizing the seasonal cycle and vertical structure of ozone in Paris, France using four years of ground based LIDAR measurements in the lowermost troposphere [J]. Atmospheric Environment, 2017,167:603-615.

[24] Su W J, Liu C, Hu Q H, et al. Characterization of ozone in the lower troposphere during the 2016G20 conference in Hangzhou [J]. Scientific Reports, 2017,7(1):1-11.

[25] He G W, Deng T, Wu Dui, et al. Characteristics of boundary layer ozone and its effect on surface ozone concentration in Shenzhen, China: A case study [J]. Science of the Total Environment, 2021,791: 148044.

[26] Zhu X W, Ma Z Q, Qiu Y L, et al. An evaluation of the interaction of morning residual layer ozone and mixing layer ozone in rural areas of the North China Plain [J]. Atmospheric Research, 2020,236:104788.

[27] 趙卉卉,王明旭,張永波.廣東省“十三五”環境保護戰略思考[J]. 環境保護科學, 2016,42(1):28-32.

Zhao H H, Wang M X, Zhang Y B, et al. Strategic ideas of the 13th five-year environment protection plan in Guangdong Province [J]. Environmental Protection Science, 2016,42(1):28-32.

[28] 沈 勁,楊土士,晏平仲,等.廣東省臭氧污染特征及其成因分析[J]. 環境科學與技術, 2020,43(12):90-95.

Shen J, Yang T S, Yan P Z, et al. Characteristics and causes of ozone pollution in Guangdong Province [J]. Environmental Science & Technology, 2020,43(12):90-95.

[29] 趙 偉,高 博,盧 清,等.2006～2019年珠三角地區臭氧污染趨勢[J]. 環境科學, 2021,42(1):97-105.

Zhao W, Gao B, Lu Q, et al. Ozone pollution trend in the Pearl River Delta Region during 2006～2019 [J]. Environmental Science, 2021,42 (1):97-105.

[30] 廣東省生態環境廳.2021廣東省生態環境狀況公報[Z].http://gdee. gd.gov.cn/attachment/0/493/493184/3927093.pdf,2022-04-26.

Guangdong Provincial Department of ecology and environment. 2021Report on the state of Guangdong provincial ecology and environment [Z]. http://gdee.gd.gov.cn/attachment/0/493/493184/ 3927093.pdf,2022-04-26.

[31] 裴成磊,謝雨彤,陳 希,等.廣州市冬季一次典型臭氧污染過程分析[J]. 環境科學, 2022,43(10):4305-4315.

Pei C l, Xie Y T, Chen X, et al. Analysis of a typical ozone pollution process in Guangzhou in winter [J]. Environmental Science, 2022, 43(10):4305-4315.

[32] 王 川,夏士勇,曹禮明,等.深圳西部城區大氣O3污染特征及超標成因[J]. 中國環境科學, 2020,40(4):1414-1420.

Wang C, Xia S Y, Cao L M, et al. Study on the characteristics and the cause of atmospheric O3pollution in western urban of Shenzhen [J]. China Environmental Science, 2020,40(4):1414-1420.

[33] 馮 凝,唐夢雪,李孟林,等.深圳市城區VOCs對PM2.5和O3耦合生成影響研究[J]. 中國環境科學, 2021,41(1):11-17.

Feng N, Tang M X, Li M L, et al. Research on the influence of VOCs on the coupling generation of PM2.5and O3in Shenzhen [J]. China Environmental Science, 2021,41(1):11-17.

[34] 高 戈,李海萍,魏 冰,等.珠江三角洲干線公路夏季臭氧污染特征[J]. 中國環境科學, 2020,40(3):929-937.

Gao G, Li H P, Wei B, et al. Ozone pollution characteristics of nation trunk highway in Pearl River Delta during summer [J]. China Environmental Science, 2020,40(3):929-937.

[35] 黃 俊,廖碧婷,吳 兌,等.廣州近地面臭氧濃度特征及氣象影響分析[J]. 環境科學學報, 2018,38(1):23-31.

Huang J, Liao B T, Wu Dui, et al. Guangzhou ground level ozone concentration characteristics and associated meteorological factors [J]. Acta Scientiae Circumstantiae, 2018,38(1):23-31.

[36] 趙 偉,高 博,劉 明,等.氣象因素對香港地區臭氧污染的影響[J]. 環境科學, 2019,40(1):55-66.

Zhao W, Gao B, Liu M, et al. Impact of meteorological factors on the ozone pollution in Hong Kong [J]. Environment Science, 2019,40(1): 55-66.

[37] 梁碧玲,張 麗,賴 鑫,等.深圳市臭氧污染特征及其與氣象條件的關系[J]. 氣象與環境學報, 2017,33(1):66-71.

Liang B L, Zhang L, Lai X, et al. Analysis of the characteristics of ozone pollution and its relationship with meteorological conditions in Shenzhen [J]. Journal of Meteorology and Environment, 2017,33(1): 66-71.

[38] 沈 勁,黃曉波,汪 宇,等.廣東省臭氧污染特征及其來源解析研究[J]. 環境科學學報, 2017,37(12):4449-4457.

Shen J, Huang X B, Wang Y, et al. Study on ozone pollution characteristics and source apportionment in Guangdong Province [J]. Acta Scientiae Circumstantiae, 2017,37(12):4449-4457.

[39] 李嫣婷,孫天樂,何 龍,等.深圳市秋季大氣臭氧立體分布特征[J]. 中國環境科學, 2020,40(5):1975-1981.

Li Y T, Sun T L, He L, et al. Vertical distribution characteristics of ozone pollution in Shenzhen in autumn [J]. China Environmental Science, 2020,40(5):1975-1981.

[40] HJ633-2012 環境空氣質量指數(AQI)技術規定(試行) [S]. HJ633-2012 Technical regulation on ambient air quality index(on trial) [S].

[41] Cressman G P. An operational objective analysis system [J]. Monthly Weather Review, 1959,87(10):367-374.

[42] Shen J, Zhang Y H, Wang X S, et al. An ozone episode over the Pearl River Delta in October 2008 [J]. Atmospheric Environment, 2015, 122:852-863.

[43] 周學思,廖志恒,王 萌,等.2013～2016年珠海地區臭氧濃度特征及其與氣象因素的關系[J]. 環境科學學報, 2019,39(1):143-153.

Zhou X S, Liao Z H, Wang M, et al. Characteristics of ozone concentration and its relationship with meteorological factors in Zhuhai during 2013～2016 [J]. Acta Scientiae Circumstantiae, 2019, 39(1):143-153.

[44] 王 磊,劉端陽,韓桂榮,等.南京地區近地面臭氧濃度與氣象條件關系研究[J]. 環境科學學報, 2018,38(4):1285-1296.

Wang L, Liu D Y, Han G R, et al. Study on the relationship between surface ozone concentrations and meteorological conditions in Nanjing, China [J]. Acta Scientiae Circumstantiae, 2018,38(4):1285-1296.

[45] 廣東省氣象局.廣東省天氣預報技術手冊[M]. 北京,氣象出版社, 2006:388-389.

Guangdong Meteorological Service. Technical manual of weather forecast for Guangdong Province [M]. Beijing, China Meteorological Press, 2006:388-389.

[46] 劉 建,吳 兌,范紹佳,等.前體物和氣象因子對珠江三角洲臭氧污染的影響[J]. 中國環境科學, 2017,37(3):813-820.

Liu J, Wu D, Fan S J, et al. Impacts of precursors and meteorological factors on ozone pollution in Pearl River Delta [J]. China Environmental Science, 2017,37(3):813-820.

[47] 王占山,李云婷,陳 添,等.北京城區臭氧日變化特征及與前體物的相關性分析[J]. 中國環境科學, 2014,34(12):3001-3008.

Wang Z S, Li Y T, Chen T, et al. Analysis on diurnal variation characteristics of ozone and correlations with its precursors in urban atmosphere of Beijing [J]. China Environmental Science, 2014, 34(12):3001-3008.

[48] 肖鐘湧,謝先全,陳穎鋒,等.粵港澳大灣區NO2污染的時空特征及影響因素分析[J]. 中國環境科學, 2020,40(5):2010-2017.

Xiao Z Y, Xie X Q, Chen Y F, et al. Temporal and spatial characteristics and influencing factors of NO2pollution over Guangdong-Hong Kong-Macao Greater Bay Area, China [J]. China Environmental Science, 2020,40(5):2010-2017.

[49] 陶麗萍,鄧 濤,吳 兌,等.廣州旱季雙高污染及消光系數垂直分布特征[J]. 中國環境科學, 2022,42(2):497-508.

Tao L P, Deng T, Wu D, et al. High aerosol and high ozone pollution and vertical distribution of extinction coefficients in Guangzhou during the dry season [J]. China Environmental Science, 2022,42(2): 497-508.

Analysis of a winter ozone pollution process in the western Pearl River Delta.

MAI Jian-hua1, YU Ling-ling2, DENG Xue-jiao3*, OU Hong-hui4, YU Xin-yang1

(1.Zhongshan Meteorological Service, Zhongshan 528400, China；2.Guangdong Meteorological Observatory, Guangzhou 510641, China；3.Guangzhou Institute of Tropical and Marine Meteorology, China Meteorological Administration, Guangzhou 510641, China；4.Zhongshan Ecology and Environment Monitoring Station of Guangdong Province, Zhongshan 528400, China)., 2023,43(12)：6235～6245

Using surface meteorological and air quality observational data, the vertical sounding meteorological data and the reanalysis data, the characteristics and cause of an ozone pollution process in Zhongshan City, located in the western part of Pearl River Delta from December 8 to 12, 2021 were analyzed. The results showed that this process was a rare winter continuous O3pollution process in recent years, during which the average O3peak concentration in pollution days was 2.3 times higher than that of non-pollution days of the month. Under the control of subtropical high and the surface cold high ridge during the pollution days, the temperature, sunshine hour and humidity in Zhongshan were relatively high, while the wind speed was relatively low, and the daytime wind direction was dominated by north wind. The downdraft, the inversion layer and the descending relative humidity in the daytime near surface enhanced the ozone pollution. There was a clear phenomenon of external pollutant transport during the pollution process, and the potential sources areas of pollution were mainly located in the northeast direction of Zhongshan. A large amount of NO2flowed into Zhongshan from upstream areas during the accumulation stage of precursors, providing sufficient precursors for the occurrence of O3pollution. During the pollution days, under favorable solar radiation conditions, the average photolysis rate of NO2was 23% higher than that of non-pollution days of that month, leading to a rapid increase in daytime O3concentration and the formation of pollution under poor diffusion conditions.

ozone pollution；meteorological conditions；transport flux；photolysis rate

X511

A

1000-6923(2023)12-6235-11

麥健華,于玲玲,鄧雪嬌,等.珠三角西部一次冬季臭氧污染過程分析 [J]. 中國環境科學, 2023,43(12):6235-6245.

Mai J H, Yu L L, Deng X J, et al. Analysis of a winter ozone pollution process in the western Pearl River Delta [J]. China Environmental Science, 2023,43(12):6235-6245.

2023-04-03

廣東省重點領域研發計劃項目(2020B1111360003);國家自然科學基金資助面上項目(42275123);廣東省氣象局科技創新團隊計劃項目(GRMCTD202003);中山市氣象局科學技術研究項目(JKT201913)

* 責任作者, 研究員, dxj@gd121.cn

麥健華(1985-),男,廣東中山人,高級工程師,碩士,主要從事環境氣象研究工作.發表論文10余篇.maigua@163.com.