吳 波,劉春瓊,張 嬌,李彥輝,陳郁兵,文 燁,鮑冰逸,杜 娟,史 凱*

COVID-19期間區(qū)域大氣高污染發(fā)生的非線性動(dòng)力機(jī)制

吳 波1,劉春瓊2,張 嬌1,李彥輝1,陳郁兵1,文 燁1,鮑冰逸1,杜 娟1,史 凱2*

(1.吉首大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,湖南 吉首 416000；2.吉首大學(xué)生物資源與環(huán)境科學(xué)學(xué)院,湖南 吉首 416000)

以長(zhǎng)株潭城市群3個(gè)城市(長(zhǎng)沙、株洲和湘潭)疫情期間(2020年1月24日～2020年5月31日)大氣PM2.5和O3小時(shí)平均濃度時(shí)間序列數(shù)據(jù)為研究對(duì)象,對(duì)污染物日變化規(guī)律、長(zhǎng)期持續(xù)性、多重分形性及自組織演化動(dòng)力學(xué)特性進(jìn)行研究.以期闡釋疫情期間高污染事件發(fā)生及演化的內(nèi)在動(dòng)力機(jī)制.首先,對(duì)3個(gè)城市PM2.5和O3質(zhì)量濃度的日變化規(guī)律進(jìn)行分析,發(fā)現(xiàn)O3呈現(xiàn)出晝高夜低的單峰型,而PM2.5日變化規(guī)律表現(xiàn)出晝低夜高單峰型,與非疫情期間的特征有所差異.進(jìn)一步,應(yīng)用消除趨勢(shì)波動(dòng)分析法(DFA)、多重分形消除趨勢(shì)波動(dòng)分析法(MFDFA)和概率統(tǒng)計(jì)分析,研究了大氣復(fù)合系統(tǒng)中PM2.5和O3質(zhì)量濃度序列的長(zhǎng)期持續(xù)性和多重分形結(jié)構(gòu).結(jié)果發(fā)現(xiàn)3個(gè)城市PM2.5和O3濃度序列均具有顯著的長(zhǎng)期持續(xù)性特征和較強(qiáng)的多重分形結(jié)構(gòu),同時(shí),應(yīng)用去趨勢(shì)互相關(guān)分析法(DCCA)方法和多重分形去趨勢(shì)互相關(guān)分析法(MFDCCA)對(duì)PM2.5和O3兩者之間的互相關(guān)性進(jìn)行了分析,發(fā)現(xiàn)PM2.5-O3之間的互相關(guān)也存在顯著的長(zhǎng)期持續(xù)性特征以及在不同時(shí)間尺度存在多重分形特征.進(jìn)一步將疫情期間得到的非線性分析結(jié)果與2018年和2019年同期非疫情期間的結(jié)果進(jìn)行了對(duì)比分析.最后,基于自組織臨界理論(SOC)探討了大氣PM2.5和O3濃度時(shí)空演化的內(nèi)在動(dòng)力規(guī)律,并結(jié)合典型區(qū)域氣象特征,闡明了SOC內(nèi)稟動(dòng)力機(jī)制可能是COVID-19疫情期間大氣高污染形成的主導(dǎo)機(jī)制之一.疫情期間大氣PM2.5和O3濃度并非分別獨(dú)自演化,而是依然保持復(fù)雜的相互作用.靜穩(wěn)氣象條件下,大氣復(fù)合污染內(nèi)部的非線性耦合作用可能達(dá)到動(dòng)力學(xué)臨界狀態(tài),將導(dǎo)致長(zhǎng)株潭城市群在疫情期間仍有發(fā)生大氣高污染的風(fēng)險(xiǎn).

日變化規(guī)律；長(zhǎng)期持續(xù)性；多重分形；冪律統(tǒng)計(jì)；自組織臨界；PM2.5；O3

2020年,為阻斷COVID-19疫情的蔓延,我國(guó)政府在全國(guó)范圍內(nèi)實(shí)施了遏制病毒傳播的封鎖措施,例如:嚴(yán)格交通管制,限制居民外出,停止集市集會(huì),停工停業(yè)停課等[1].從環(huán)境學(xué)角度來(lái)看,COVID-19疫情爆發(fā)可以視為一場(chǎng)代價(jià)巨大的環(huán)境學(xué)實(shí)驗(yàn),為研究疫情期間大氣污染物的演化規(guī)律提供了一個(gè)獨(dú)特的機(jī)會(huì).

在疫情防控期間,工廠和企業(yè)停產(chǎn)停工、城市交通量銳減導(dǎo)致一次污染物排放量大幅度減少[2].但COVID-19疫情期間環(huán)境狀況并未得到顯著改善,部分地區(qū)高污染事件還時(shí)有發(fā)生.例如,Sicard等[3]定量分析武漢市COVID-19疫情期間人為排放污染物的減少對(duì)空氣污染的影響.結(jié)果發(fā)現(xiàn)與2017～2019年同期相比,2020年疫情期間武漢市O3日平均濃度增加了36%.Huang等[4]發(fā)現(xiàn)2020年COVID-19疫情期間,雖然一次污染物的排放量大量減低,但中國(guó)東部地區(qū)嚴(yán)重霧霾污染事件并未能得到有效減少.Chang等[5]通過(guò)模擬分析上海市2020年與2019年春節(jié)期間氣溶膠變化趨勢(shì),發(fā)現(xiàn)2020年春節(jié)期間(疫情期間)PM2.5中的二次氣溶膠分?jǐn)?shù)比2019年春節(jié)期間的上升了16%.相關(guān)研究表明在工業(yè)生產(chǎn)和交通排放銳減的情況下,疫情期間高污染事件的發(fā)生顯然不是人為排放因素導(dǎo)致的,而極有可能是由其內(nèi)在演化的動(dòng)力學(xué)機(jī)制所主導(dǎo).

近年來(lái),長(zhǎng)株潭城市群高濃度、強(qiáng)氧化性的大氣復(fù)合污染態(tài)勢(shì)明顯[6-7].在大氣復(fù)合污染中,PM2.5與O3之間的相互作用關(guān)系及其形成機(jī)理存在著復(fù)雜的關(guān)聯(lián)性,二者不僅具有共同的前體物(如NO和揮發(fā)性有機(jī)物(VOCs)),而且還能通過(guò)多種途徑相互作用影響(如非均相反應(yīng)等).前人研究發(fā)現(xiàn),二次氣溶膠形成的主要化學(xué)途徑包括氣相光化學(xué)氧化反應(yīng)、顆粒物表面的非均相反應(yīng)以及顆粒物內(nèi)部的液相氧化反應(yīng)[8].在高濃度O3的催化氧化作用下,SO2和NO通過(guò)均相和非均相反應(yīng)生成硫酸鹽和硝酸鹽[9-12].同時(shí),VOCs與OH、NO3自由基及O3等氧化物通過(guò)光氧化反應(yīng)也能生成二次氣溶膠[13].當(dāng)細(xì)顆粒物濃度上升時(shí),細(xì)顆粒物能通過(guò)散射、吸收紫外光輻射及改變大氣紫外光輻射強(qiáng)度直接影響大氣中O3生成.此外,大氣中細(xì)顆粒物還可以通過(guò)影響云的光學(xué)厚度、云滴的有效半徑及云滴數(shù)濃度等間接影響O3生成[14].大氣復(fù)合系統(tǒng)中PM2.5和O3之間的演化,不僅受到微觀物理化學(xué)等機(jī)制的作用,同時(shí)還受氣象、地形等因素的影響.由于PM2.5和O3之間的演化涉及眾多復(fù)雜因素(如太陽(yáng)光輻射、氣象等)的控制,所以其二者濃度波動(dòng)演化具有非線性、非平穩(wěn)、非周期等特征[15].傳統(tǒng)統(tǒng)計(jì)學(xué)方法不能有效地消除非平穩(wěn)、復(fù)雜性等特征所導(dǎo)致的偽相關(guān)現(xiàn)象,因此很難通過(guò)傳統(tǒng)統(tǒng)計(jì)學(xué)方法準(zhǔn)確地刻畫(huà)出具有非線性特征的污染物內(nèi)在演化規(guī)律特征.分形方法能夠從宏觀、整體上分析大氣污染物濃度序列的復(fù)雜非線性特征,及大氣污染物隨時(shí)間演化的內(nèi)在標(biāo)度規(guī)律.因此,應(yīng)用非線性分形方法研究COVID-19疫情期間大氣復(fù)合系統(tǒng)中PM2.5和O3濃度時(shí)空演化非線性規(guī)律及二者相互作用的時(shí)間尺度特征具有重要的科學(xué)意義.

本研究重點(diǎn)以長(zhǎng)株潭城市群長(zhǎng)沙、株洲和湘潭3個(gè)城市疫情期間(2020年1月24日～2020年5月31日)大氣PM2.5和O3小時(shí)平均濃度時(shí)間序列為研究對(duì)象.首先應(yīng)用消除趨勢(shì)波動(dòng)分析法(DFA)和去趨勢(shì)互相關(guān)分析法(DCCA)研究長(zhǎng)沙、株洲和湘潭3個(gè)城市疫情期間PM2.5和O3濃度序列的時(shí)間演化規(guī)律及內(nèi)在動(dòng)力學(xué)機(jī)制.其次,應(yīng)用多重分形消除趨勢(shì)波動(dòng)分析法(MFDFA)和多重分形去趨勢(shì)互相關(guān)分析法(MFDCCA)分析3個(gè)城市疫情期間PM2.5和O3濃度序列多重分形性及二者之間的互相關(guān)性的多重分形特征,并與2018年和2019年同期數(shù)據(jù)結(jié)果進(jìn)行對(duì)比分析.最后,根據(jù)自組織臨界理論(SOC)闡明疫情期間長(zhǎng)株潭城市群大氣復(fù)合污染高污染發(fā)生的演化特征和內(nèi)在機(jī)制,以期為城市大氣高污染事件的預(yù)測(cè)、預(yù)警和風(fēng)險(xiǎn)評(píng)價(jià)提供科學(xué)參考.

1 材料與方法

1.1 材料

長(zhǎng)株潭城市群位于湖南省中東部(26°03￠～ 28o40￠N,111o54′～114o15′E),包括長(zhǎng)沙、株洲和湘潭3個(gè)主要城市.該城市群位于雪峰山脈、羅霄山脈之間的湘江谷地,中心地勢(shì)低,周邊以山地丘陵地貌為主.城市群屬于亞熱帶季風(fēng)氣候區(qū),四季分明,夏熱期長(zhǎng),降水豐沛,年平均降雨量1400mm左右.該城市群作為湖南省工業(yè)化、城鎮(zhèn)化發(fā)展的核心區(qū)域,其大氣污染物排放量高,且地形條件不利于污染物的擴(kuò)散[16].因此,該城市群內(nèi)很容易發(fā)生高污染事件.本文選擇長(zhǎng)沙、株洲和湘潭3個(gè)城市2020年1月24日～2020年5月31日大氣PM2.5、O3小時(shí)平均濃度作為主要研究對(duì)象(圖1).實(shí)時(shí)監(jiān)測(cè)數(shù)據(jù)來(lái)源于全國(guó)城市空氣質(zhì)量實(shí)時(shí)發(fā)布平臺(tái).由于儀器校準(zhǔn)、維護(hù)、停電和故障等原因,可能會(huì)造成短時(shí)間內(nèi)各城市監(jiān)測(cè)子站內(nèi)的O3、PM2.5質(zhì)量濃度數(shù)據(jù)有不同程度的缺失,但數(shù)據(jù)的缺失比率均不超過(guò)0.9%,且缺失數(shù)據(jù)在研究時(shí)段幾乎均勻分布,不存在連續(xù)缺失的現(xiàn)象.并且利用缺失數(shù)據(jù)前后時(shí)刻的濃度平均值對(duì)缺失數(shù)據(jù)進(jìn)行填補(bǔ).

圖1 2020年1月24日～2020年5月31日長(zhǎng)沙、株洲和湘潭3個(gè)城市O3和PM2.5原始數(shù)據(jù)序列圖

1.2 研究方法

1.2.1 去趨勢(shì)互相關(guān)分析法(DCCA) DCCA能有效地避免非平穩(wěn)序列所引起的偽相關(guān)的影響,現(xiàn)已成為研究?jī)山M非平穩(wěn)序列之間相關(guān)性的最科學(xué)有效方法之一.具體算法如下.

首先,對(duì)兩類污染物濃度時(shí)間序列()和(),=1,2,…,進(jìn)行積分.

然后,對(duì)給定的時(shí)間尺度,用積分序列減趨勢(shì)信號(hào),得到兩組數(shù)據(jù)殘余信號(hào)的協(xié)方差.

式中:()為時(shí)間尺度下()和()污染物殘余信號(hào)的協(xié)方差.

最后,取不同的時(shí)間尺度,重復(fù)上述過(guò)程,得到不同時(shí)間尺度下的().若在雙對(duì)數(shù)坐標(biāo)下ln～ ln()滿足線性關(guān)系,則存在冪律關(guān)系

式中:表示DCCA指數(shù), 即自相似性參數(shù).

其中,值得注意的是,如果研究序列()與()為同一時(shí)間序列時(shí),則DCCA算法等效于DFA. DCCA指數(shù)反映了兩組非平穩(wěn)時(shí)間序列之間的長(zhǎng)期互相關(guān)性.≈0.5時(shí),表明分析的兩組時(shí)間序列中不存在長(zhǎng)期互相關(guān)性,即兩組時(shí)間序列任何時(shí)刻均是完全相互獨(dú)立;≠0.5時(shí),表示兩組時(shí)間序列中存在長(zhǎng)期互相關(guān)性,這意味著兩組序列之間的數(shù)據(jù)值之間不是完全獨(dú)立的;當(dāng)>0.5時(shí),表示兩組時(shí)間序列之間存在長(zhǎng)期正相關(guān)性,即如果過(guò)去一段時(shí)間內(nèi)一類污染物呈增加(減小)的趨勢(shì),則會(huì)導(dǎo)致另一組污染物呈現(xiàn)出相同的趨勢(shì),且這種影響在一定時(shí)間尺度上表現(xiàn)出無(wú)標(biāo)度的冪律衰減形式.當(dāng)<0.5則表示兩組污染物濃度序列具有反持續(xù)性的長(zhǎng)期互相關(guān)性.

1.2.2 多重分形去趨勢(shì)互相關(guān)分析(MFDCCA) MFDCCA是由Zhou[17]提出的一種可以識(shí)別2個(gè)非平穩(wěn)時(shí)間序列之間的相互關(guān)系和多重分形特征的新統(tǒng)計(jì)方法[23].其算法具體如下.

首先確定2個(gè)污染濃度時(shí)間序列()和(),=1,2,…,,其中是時(shí)間序列的長(zhǎng)度.其序列()和()的累積離差序列如下:

其次,將累積離差序列()和()劃分為具有相同長(zhǎng)度N=/的互不相交的區(qū)間,為時(shí)間尺度.由于時(shí)間序列長(zhǎng)度通常不是時(shí)間尺度的整倍數(shù),為了避免數(shù)據(jù)遺漏,將其序列從尾部向前重復(fù)上述相同的過(guò)程.因此,得到2N個(gè)區(qū)間.

和

每個(gè)區(qū)間,=+1,+2,…,2.然后,對(duì)所有分段進(jìn)行平均,得到(≠0)階波動(dòng)函數(shù)().

最后,計(jì)算不同時(shí)間尺度下的波動(dòng)函數(shù)().如果2個(gè)時(shí)間序列之間存在長(zhǎng)期冪律互相關(guān)關(guān)系,則波動(dòng)函數(shù)與時(shí)間尺度將存在冪律關(guān)系

式中:()為廣義赫斯特指數(shù).

特別需要注意的是,如果系列()與()相同,則MFDCCA算法等效于MFDFA.當(dāng)=2時(shí),廣義互相關(guān)指數(shù)(2)可以表征赫斯特指數(shù)所具備的性質(zhì).(2)=0.5表示兩序列之間不存在相互關(guān)系.(2)>0.5表示兩序列之間呈現(xiàn)出長(zhǎng)期冪律互相關(guān),且以冪律形式呈現(xiàn)出正相關(guān)關(guān)系,這意味某一個(gè)污染物增大可能會(huì)在未來(lái)一定時(shí)間尺度內(nèi)造成另一個(gè)污染物濃度增大.這也表明兩個(gè)污染物之間的相關(guān)性具有長(zhǎng)期持續(xù)性特征.而(2)<0.5則與此相反.

如果()不依賴的變化而變化,則兩個(gè)序列之間的相互關(guān)系為單一分形.如果()表現(xiàn)出隨增加而單調(diào)遞減形式,則兩個(gè)序列之間的互相關(guān)性具有典型的多重分形特征.在多重分形模式下,()隨著的變化表現(xiàn)為在不同尺度下,兩個(gè)序列之間的小波動(dòng)和大波動(dòng)互相關(guān)的長(zhǎng)期持續(xù)性特征具有異質(zhì)性行為.多重分形的強(qiáng)度可通過(guò)()的范圍來(lái)計(jì)算.

當(dāng)?越大, 多重分形性越強(qiáng).

1.2.3 污染濃度概率統(tǒng)計(jì)分析 過(guò)去研究中,前人已經(jīng)證明了大氣污染濃度波動(dòng)并非完全隨機(jī)而是在統(tǒng)計(jì)上遵從典型的冪律分布[18-19].如果某一次污染事件的發(fā)生規(guī)模在統(tǒng)計(jì)上呈現(xiàn)出冪律關(guān)系,則可認(rèn)為污染系統(tǒng)出現(xiàn)自組織臨界(SOC)行為的一個(gè)顯著標(biāo)致.通常表現(xiàn)為污染事件發(fā)生的頻率隨著事件規(guī)模大小呈現(xiàn)冪指數(shù)下降.如果污染濃度序列波動(dòng)具有相類似的特征, 則需滿足下列關(guān)系式.

式中:表示某一污染物濃度值,μg/m3;為標(biāo)度指數(shù),為大于某一污染物濃度值(0)的累計(jì)次數(shù).

1.2.4 自組織臨界理論 Bak等[20]于1987年提出的自組織臨界(SOC)這一概念,從整體、宏觀角度闡釋了復(fù)雜系統(tǒng)中長(zhǎng)期持續(xù)性特征和冪律統(tǒng)計(jì)分布的產(chǎn)生動(dòng)力學(xué)機(jī)制.該類SOC系統(tǒng)是由數(shù)目龐大發(fā)生短程近鄰相互作用的組元所構(gòu)成.在系統(tǒng)內(nèi)部和外部因素共同作用下,該系統(tǒng)將自發(fā)地向臨界狀態(tài)演化并最終穩(wěn)定在該臨界狀態(tài).當(dāng)系統(tǒng)穩(wěn)定在臨界狀態(tài)時(shí),外界任何一個(gè)小事件(物質(zhì)或能量輸入)都將可能引起影響眾多組元的連鎖反應(yīng),最終造成大規(guī)模事件發(fā)生.在復(fù)雜SOC系統(tǒng)中大規(guī)模事件發(fā)生的概率往往比小規(guī)模事件發(fā)生的概率小,其大規(guī)模事件發(fā)生主要是由系統(tǒng)內(nèi)眾多大大小小規(guī)模的連鎖反應(yīng)所導(dǎo)致,具有標(biāo)度不變性.因此,系統(tǒng)內(nèi)的時(shí)空關(guān)聯(lián)函數(shù)在宏觀上涌現(xiàn)出冪律標(biāo)度特征及長(zhǎng)期持續(xù)性等特征.這樣,冪律分布統(tǒng)計(jì)特征及長(zhǎng)期持續(xù)性特征可以作為復(fù)雜系統(tǒng)內(nèi)存在SOC特性的重要依據(jù)[21].沙堆模型是理解SOC理論的關(guān)鍵性模型.緩慢的向沙堆中投加沙粒時(shí),當(dāng)沙粒累積到一定程度時(shí)沙堆就會(huì)產(chǎn)生崩塌,這類崩塌屬于一種非線性的臨界行為.從宏觀上來(lái)看,這種臨界行為產(chǎn)生的宏觀效應(yīng)具有冪律統(tǒng)計(jì)規(guī)律.

2 結(jié)果與分析

2.1 PM2.5和O3濃度日變化規(guī)律

從圖2中可以發(fā)現(xiàn),3個(gè)城市PM2.5和O3質(zhì)量濃度的日變化規(guī)律曲線基本一致.3個(gè)城市O3質(zhì)量濃度日變化規(guī)律與近地面大氣光化學(xué)過(guò)程密切相關(guān),均呈現(xiàn)出明顯的晝高夜低單峰型變化規(guī)律.3個(gè)城市O3質(zhì)量濃度均于下午16:00左右達(dá)到一天中的最高值,隨后濃度慢慢下降(圖2),直到次日日出時(shí)分(07:00左右)達(dá)到一天最小濃度值.O3質(zhì)量濃度于下午16:00點(diǎn)左右達(dá)到最高值,這主要是由于O3對(duì)太陽(yáng)輻射及溫度的響應(yīng)具有明顯延遲性,即前體物通過(guò)復(fù)雜的光化學(xué)反應(yīng)生成臭氧需要一定時(shí)間[22].但16:00左右以后太陽(yáng)輻射強(qiáng)度降低,光化學(xué)反應(yīng)過(guò)程中O3的生成量低于氧化反應(yīng)中(例如NO氧化成NO2)O3的消耗量.因此,O3質(zhì)量濃度在16:00點(diǎn)以后會(huì)逐漸降低,直至次日07:00點(diǎn)左右達(dá)到一天中最低濃度值[23-24].此外,3個(gè)城市PM2.5質(zhì)量濃度日變化規(guī)律特征均表現(xiàn)出凌晨(03:00點(diǎn)左右)最高,隨后濃度逐漸下降直到下午(16:00點(diǎn)左右)達(dá)到一天中的最低點(diǎn).這是因?yàn)镻M2.5日變化濃度主要是受太陽(yáng)輻射和其重力作用影響.PM2.5等顆粒物經(jīng)過(guò)沉降作用,晚上近地面的PM2.5濃度不斷積累,于凌晨03:00點(diǎn)左右達(dá)到一天中的最大值.但隨著太陽(yáng)輻射的出現(xiàn)和增強(qiáng),近地面溫度上升,近地面附近將形成暖氣團(tuán),在暖氣團(tuán)的攜帶下PM2.5將會(huì)一起上升.因此,PM2.5質(zhì)量濃度會(huì)隨著太陽(yáng)輻射增強(qiáng)逐漸減少.PM2.5污染物經(jīng)過(guò)一天的輻射、升溫、上升,直至下午16:00點(diǎn)左右達(dá)到一天的最低值[25].上述結(jié)果與以前的研究存在一定的差異,前人研究發(fā)現(xiàn)非疫情期間PM2.5濃度呈現(xiàn)出雙峰型日變化規(guī)律,其主要原因可能是由于上下班高峰期的機(jī)動(dòng)車(chē)尾氣的排放[26-27].而COVID-19疫情期間PM2.5呈現(xiàn)出單峰型日變化規(guī)律主要可能是由于人們出行受限機(jī)動(dòng)車(chē)尾氣排放顯著減少的原因.

圖2 2020年1月24日～2020年5月31日長(zhǎng)株潭城市群3個(gè)城市PM2.5和O3濃度日變化規(guī)律

2.2 PM2.5和O3長(zhǎng)期持續(xù)性及長(zhǎng)期相關(guān)性特征

應(yīng)用消除趨勢(shì)波動(dòng)分析法(DFA)分別對(duì)長(zhǎng)沙、株洲和湘潭3個(gè)城市的PM2.5和O3濃度時(shí)間序列進(jìn)行分析,結(jié)果如圖3所示.COVID-19疫情期間長(zhǎng)株潭城市群長(zhǎng)沙、株洲和湘潭3個(gè)城市PM2.5和O3的ln()～ln()關(guān)系呈現(xiàn)出顯著的線性關(guān)系.并利用最小二乘法對(duì)3個(gè)城市PM2.5和O3濃度時(shí)間序列的DFA指數(shù)進(jìn)行線性估計(jì),其指數(shù)分別為長(zhǎng)沙((PM2.5)=0.968,(O3)=0.771),株洲((PM2.5)= 0.997,(O3)=0.812),湘潭((PM2.5)=1.000,(O3)= 0.768).污染物濃度序列的長(zhǎng)期持續(xù)性特征能被DFA指數(shù)直接描述,且3個(gè)城市PM2.5和O3的指數(shù)均高于0.5,表示其污染物具有較強(qiáng)的長(zhǎng)期持續(xù)性特征.長(zhǎng)期持續(xù)性反映了COVID-19疫情期間長(zhǎng)沙、株洲和湘潭3個(gè)城市的PM2.5和O3濃度時(shí)間序列變化前后的相關(guān)性規(guī)律特征.換言之,PM2.5和O3濃度時(shí)間序列在一定時(shí)間尺度范圍內(nèi)其相關(guān)性變化并不嚴(yán)格遵循馬爾科夫過(guò)程,而是遵循冪律衰減規(guī)律.其指數(shù)越大,污染物的長(zhǎng)期持續(xù)性越強(qiáng),則意味著如果過(guò)去一段時(shí)間內(nèi)污染源排放的PM2.5(O3)增加將會(huì)導(dǎo)致未來(lái)一段時(shí)間內(nèi)PM2.5(O3)濃度持續(xù)增加.即過(guò)去時(shí)刻的PM2.5(O3)濃度會(huì)對(duì)未來(lái)一定時(shí)間尺度內(nèi)PM2.5(O3)質(zhì)量濃度積累或消散產(chǎn)生持續(xù)性的影響.

圖3 疫情期間長(zhǎng)沙、株洲和湘潭3個(gè)城市PM2.5和O3的DFA以及PM2.5-O3的DCCA圖

為了定量分析PM2.5和O3之間相互作用對(duì)二者濃度演化過(guò)程中的影響,應(yīng)用去趨勢(shì)互相關(guān)分析(DCCA)方法對(duì)長(zhǎng)株潭城市群3個(gè)城市的PM2.5和O3濃度相互作用進(jìn)行分析,其結(jié)果如圖3所示.長(zhǎng)沙、株洲和湘潭3個(gè)城市在研究時(shí)段內(nèi)均表現(xiàn)出相似演化特征,呈現(xiàn)出較好的一致性.該結(jié)果與PM2.5和O3的DFA分析結(jié)果類似,研究期間內(nèi)PM2.5和O3之間互相關(guān)的ln()～ln()也呈現(xiàn)出顯著的線性關(guān)系.其中長(zhǎng)沙、株洲和湘潭DCCA標(biāo)度指數(shù)分別為0.893、0.910、0.900.長(zhǎng)株潭城市群3個(gè)城市的分析結(jié)果表明,在COVID-19疫情期間內(nèi)PM2.5和O3濃度的相關(guān)性均表現(xiàn)出較強(qiáng)的長(zhǎng)期持續(xù)性特征.這意味著在過(guò)去的一段時(shí)間內(nèi)3個(gè)城市PM2.5(O3)污染物排放增加將導(dǎo)致未來(lái)一段時(shí)間內(nèi)3個(gè)城市O3(PM2.5)濃度持續(xù)增加.

由上述結(jié)果可知,長(zhǎng)株潭城市群3個(gè)城市PM2.5和O3濃度及PM2.5-O3互相關(guān)性都有顯著的長(zhǎng)期持續(xù)性特征.這表明在各個(gè)城市大氣系統(tǒng)中PM2.5和O3濃度演化不僅受到自身濃度波動(dòng)的影響,同時(shí)還受到二者之間的相互作用的影響.

2.3 O3和PM2.5濃度序列多重分形特征

為了進(jìn)一步研究局域時(shí)間尺度下O3濃度和PM2.5濃度序列的非線性特征及二者之間非線性互相關(guān)關(guān)系,對(duì)COVID-19疫情期間長(zhǎng)株潭城市群3個(gè)城市O3濃度,PM2.5濃度序列和O3-PM2.5互相關(guān)性進(jìn)行多重分形消除趨勢(shì)波動(dòng)分析法(MFDFA)和多重分形去趨勢(shì)互相關(guān)分析法(MFDCCA)分析.

從圖4、5、6可知,在整個(gè)時(shí)間尺度上,長(zhǎng)沙、株洲和湘潭3個(gè)城市PM2.5、O3濃度和PM2.5-O3互相關(guān)性的雙對(duì)數(shù)坐標(biāo)圖中不同階下ln()均能很好地?cái)M合去趨勢(shì)波動(dòng)函數(shù)lnF().其中,不同(=0,±5,±10,±20)階的ln()～ln()均表現(xiàn)出良好的線性關(guān)系且回歸線收斂.這表示3個(gè)城市PM2.5、O3濃度波動(dòng)和PM2.5-O3之間的互相關(guān)演化在不同尺度上均具有不同冪律形式的長(zhǎng)期持續(xù)性特征.這也說(shuō)明在一定時(shí)間尺度上,PM2.5、O3濃度和PM2.5-O3互相關(guān)性不服從經(jīng)典的馬爾可夫過(guò)程,而是以冪律形式隨時(shí)間緩慢衰減.從圖4～6中階廣義Hurst指數(shù)曲線變化來(lái)看,長(zhǎng)沙、株洲和湘潭3個(gè)城市的O3、PM2.5濃度和PM2.5-O3互相關(guān)性的()均表現(xiàn)出隨著值的增大而減小的規(guī)律,這證明3個(gè)城市中PM2.5、O3濃度和PM2.5-O3互相關(guān)性的長(zhǎng)期持續(xù)特征均具有多重分形標(biāo)度不變性的結(jié)構(gòu).

進(jìn)一步依據(jù)(10)式可計(jì)算出多重分形的強(qiáng)度?.在多重分形譜中,?值越高,則表示序列的多重分形強(qiáng)度越高,即污染物濃度時(shí)間序列分布越不均勻,其分形特征越明顯.計(jì)算得到的疫情期間各城市O3、PM2.5濃度和PM2.5-O3互相關(guān)性的?h值見(jiàn)表1.從整體上看,各城市PM2.5濃度序列的多重分形強(qiáng)度最大.這種情況的原因可能是,本研究探討的疫情時(shí)段主要涉及冬季和春季兩個(gè)季節(jié),在靜風(fēng)頻率和逆溫作用下PM2.5濃度演化波動(dòng)不穩(wěn)定,所以PM2.5表現(xiàn)出更強(qiáng)的多重分形特征.此外,3個(gè)城市的各非線性特征分析結(jié)果具有顯著的一致性,這充分表明長(zhǎng)沙、株洲和湘潭3個(gè)城市污染具有區(qū)域性特征.

根據(jù)上述結(jié)果可知,長(zhǎng)沙、株洲和湘潭3個(gè)城市在COVID-19疫情期間PM2.5、O3濃度和PM2.5-O3互相關(guān)性的多重分形特征呈現(xiàn)出顯著的動(dòng)力學(xué)同步性.這表示了3個(gè)城市中PM2.5和O3之間可能存在相似的非線性動(dòng)力學(xué)演化機(jī)制.此外,3個(gè)城市PM2.5、O3濃度均具有較強(qiáng)的多重分形特征,且PM2.5-O3之間的互相關(guān)也存在顯著的多重分形特征.這表明了疫情期間大氣系統(tǒng)中PM2.5和O3濃度并非分別獨(dú)自演化,而是二者演化過(guò)程中還存在著復(fù)雜的相互作用關(guān)系.

圖5 疫情期間株洲PM2.5、O3和O3-PM2.5雙對(duì)數(shù)波動(dòng)曲線及q階廣義Hurst指數(shù)圖

圖6 疫情期間湘潭PM2.5、O3和O3-PM2.5雙對(duì)數(shù)波動(dòng)曲線及q階廣義Hurst指數(shù)圖

2.4 O3和PM2.5污染濃度概率統(tǒng)計(jì)分布

為了進(jìn)一步驗(yàn)證COVID-19疫情期間PM2.5、O3污染演化內(nèi)在動(dòng)力學(xué)機(jī)制,根據(jù)公式(11)計(jì)算了長(zhǎng)沙、株洲和湘潭3個(gè)城市PM2.5、O3小時(shí)平均濃度的累計(jì)頻度分布(圖7).利用最小二乘法對(duì)O3和PM2.5污染濃度數(shù)據(jù)點(diǎn)進(jìn)行線性回歸,結(jié)果發(fā)現(xiàn)在一定區(qū)間內(nèi),長(zhǎng)沙、株洲和湘潭3個(gè)城市均呈現(xiàn)出顯著的冪律分形結(jié)構(gòu).3個(gè)城市冪律分形結(jié)構(gòu)的標(biāo)度指數(shù)分別為O3(長(zhǎng)沙=1.651,株洲=2.531,湘潭=1.609), PM2.5(長(zhǎng)沙=3.071,株洲=2.809,湘潭=3.574).各城市兩類污染物的標(biāo)度指數(shù)有所差異,主要原因可能是雖然長(zhǎng)株潭城市群具有區(qū)域污染特征,但各個(gè)城市地理位置,大氣環(huán)境等還是存在一定的差異性,導(dǎo)致不同城市同一類污染物的冪律分形結(jié)構(gòu)有所不同.O3和PM2.5濃度波動(dòng)呈現(xiàn)出的冪律分形結(jié)構(gòu),表明COVID-19疫情期間長(zhǎng)沙、株洲和湘潭3個(gè)城市O3和PM2.5濃度演化不存在典型的特征濃度,而是具有標(biāo)度不變性特征.在復(fù)雜系統(tǒng)中,穩(wěn)定的冪律分布結(jié)構(gòu)可以作為復(fù)雜系統(tǒng)演化到自組織臨界(SOC)狀態(tài)的一個(gè)顯著標(biāo)志.因此,COVID-19疫情期間PM2.5、O3污染時(shí)空演化可能具有SOC特性.

圖7 長(zhǎng)株潭城市群O3和PM2.5污染小時(shí)平均濃度的概率統(tǒng)計(jì)分布

2.5 歷史同期數(shù)據(jù)對(duì)比分析結(jié)果

為了進(jìn)一步說(shuō)明疫情期間大氣污染的非線性演化特征,本文將非疫情期間(2018年1月24日～2018年5月31日和2019年1月24日～2019年5月31日)與疫情期間(2020年1月24日～2020年5月31日)O3和PM2.5污染的分形結(jié)果進(jìn)行對(duì)比分析,其結(jié)果如表1所示.對(duì)長(zhǎng)株潭城市群3個(gè)城市來(lái)說(shuō),2018年和2019年的分形參數(shù)數(shù)值大小均表現(xiàn)出類似的結(jié)果,而2020年的結(jié)果有較大的變化.主要特征差異表現(xiàn)為,各城市2020年O3濃度演化的長(zhǎng)期持續(xù)性指數(shù)顯著大于2018年和2019年,而2020年大氣污染演化的其余分形參數(shù)均小于2018年和2019年.

表1 長(zhǎng)株潭城市群各城市3a同期分形參數(shù)結(jié)果

疫情期間,人類活動(dòng)排放的劇減也導(dǎo)致消耗O3的前體物質(zhì)(如NO等)的排放也相應(yīng)減少,這也使得大氣環(huán)境中O3污染濃度的累積得到增強(qiáng),表現(xiàn)為疫情期間O3污染的長(zhǎng)期持續(xù)性較強(qiáng).而在非疫情期間,大氣環(huán)境中PM2.5和O3污染除了自然演化之外,還受到人類活動(dòng)排放(如工業(yè)、生活、交通等)的影響,這也就增加了其波動(dòng)演化的非線性復(fù)雜性,表現(xiàn)為非疫情期間反映大氣污染演化非線性特征的分形參數(shù)數(shù)值均較高.

2.6 討論

自組織臨界態(tài)(SOC)理論能從宏觀整體性的角度為復(fù)雜系統(tǒng)演化中產(chǎn)生的長(zhǎng)期持續(xù)性和冪律分形結(jié)構(gòu)提供內(nèi)在動(dòng)力機(jī)制解釋.前人研究中,Shi等[28]基于自組織臨界數(shù)值模型對(duì)PM10、SO2和NO2濃度進(jìn)行分析,清楚的闡釋了3類污染物濃度演化過(guò)程中呈現(xiàn)出分形特征的內(nèi)在動(dòng)力學(xué)機(jī)制.蘇蓉等[29]應(yīng)用R/S方法、功率譜、強(qiáng)度-頻度法及SOC理論對(duì)上海市PM10污染進(jìn)行分析,發(fā)現(xiàn)SOC內(nèi)稟機(jī)制是導(dǎo)致其濃度演化呈現(xiàn)出長(zhǎng)期持續(xù)性的主要原因.Shi等[30]提出了一種新的PM2.5自組織臨界演化模型,并進(jìn)一步驗(yàn)證了SOC內(nèi)稟機(jī)制是PM2.5演化中長(zhǎng)期持續(xù)性和冪律統(tǒng)計(jì)分布的產(chǎn)生根源.相關(guān)研究已表明城市空氣污染演化的內(nèi)在動(dòng)力學(xué)機(jī)制可能主要受SOC行為主導(dǎo).

本文進(jìn)一步研究發(fā)現(xiàn)COVID-19疫情期間長(zhǎng)株潭城市群大氣復(fù)合污染中PM2.5和O3動(dòng)力學(xué)過(guò)程,其特征與SOC機(jī)制特征高度相似,具體闡述如下.

在大氣復(fù)合污染系統(tǒng)中存在眾多發(fā)生短程相互作用的組元(O3、PM2.5和天氣氣象要素等),其中短程相互作用主要包括兩個(gè)方面,其一,O3可以通過(guò)對(duì) OH、H2O2、RCHO等氧化劑濃度的控制來(lái)影響二次氣溶膠的生成.其二,二次氣溶膠能通過(guò)影響大氣動(dòng)力學(xué)、光解速率、云光學(xué)厚度和非均相反應(yīng)過(guò)程來(lái)直接或間接影響O3生成.

大氣復(fù)合污染系統(tǒng)作為一個(gè)開(kāi)放性復(fù)雜巨系統(tǒng),外界能持續(xù)性地向系統(tǒng)內(nèi)輸入物質(zhì)或能量(大氣中揮發(fā)性有機(jī)物(VOCs)和氮氧化物(NO)將通過(guò)氧化反應(yīng)生成O3或某些大氣活性氣體經(jīng)過(guò)復(fù)雜的化學(xué)反應(yīng)生成二次氣溶膠等).同時(shí),該系統(tǒng)也會(huì)以某種方式、途徑向系統(tǒng)外耗散其能量或物質(zhì)(O3濃度上升將增強(qiáng)大氣氧化能力,加劇光化學(xué)反應(yīng),該過(guò)程O(píng)3迅速消耗.二次氣溶膠可通過(guò)降雨洗刷作用脫離系統(tǒng)或在大風(fēng)中碰并、凝結(jié)等過(guò)程形成更大粒徑的顆粒物脫離大氣系統(tǒng)).最終,大氣系統(tǒng)將在近地面內(nèi)外界輸入的能量或物質(zhì)與系統(tǒng)內(nèi)部耗散的能量或物質(zhì)形成一個(gè)穩(wěn)定的耗散結(jié)構(gòu)體系.

靜穩(wěn)天氣條件下,大氣復(fù)合污染系統(tǒng)會(huì)在其短程近鄰相互作用下系統(tǒng)將會(huì)自發(fā)地向臨界狀態(tài)演化.當(dāng)系統(tǒng)達(dá)到臨界狀態(tài)時(shí),大氣復(fù)合系統(tǒng)的時(shí)間演化過(guò)程中系統(tǒng)將鎖定并維持穩(wěn)定狀態(tài).

當(dāng)大氣復(fù)合污染系統(tǒng)自發(fā)演化到臨界狀態(tài)并維持著臨界態(tài)時(shí),系統(tǒng)外界任何微小的擾動(dòng)(天氣條件等)均可能造成該系統(tǒng)發(fā)生全局性連鎖反應(yīng).因此,系統(tǒng)局部區(qū)域?qū)Τ跏紬l件(天氣氣象、地形等)極其敏感.同時(shí),在系統(tǒng)向臨界狀態(tài)演化時(shí),系統(tǒng)內(nèi)眾多短程近鄰相互作用(O3通過(guò)控制氧化性間接影響二次氣溶膠的生成.而二次氣溶膠可通過(guò)改變光化學(xué)通量來(lái)影響O3生成)組元(O3、PM2.5天氣氣象要素等)將演化出長(zhǎng)期持續(xù)性特征.即一次污染事件發(fā)生的規(guī)模與其過(guò)去一段時(shí)間內(nèi)污染事件發(fā)生的規(guī)模有著密切關(guān)系.最終導(dǎo)致大氣復(fù)合系統(tǒng)耗散外界輸入和輸出物質(zhì)或能量的時(shí)空關(guān)聯(lián)函數(shù)呈現(xiàn)出冪律統(tǒng)計(jì)分布或長(zhǎng)期持續(xù)性特征.

此外,根據(jù)同期長(zhǎng)株潭城市群地面觀測(cè)數(shù)據(jù)統(tǒng)計(jì)發(fā)現(xiàn),COVID-19疫情期間長(zhǎng)株潭城市群94%以上天氣處于靜小風(fēng)環(huán)境,近地面最大瞬時(shí)風(fēng)速不高于3.3m/s.疫情期間該城市群地區(qū)水平輸送和垂直擴(kuò)散條件較差.同時(shí),長(zhǎng)株潭城市群位于雪峰山脈,羅霄山脈之間的湘江谷地,中心地勢(shì)低,周邊以山地丘陵地貌為主,地形條件不利于污染物的擴(kuò)散.區(qū)域內(nèi)不利于污染物自身稀釋和擴(kuò)散的靜穩(wěn)天氣為大氣復(fù)合污染演化的SOC機(jī)制建立了穩(wěn)定的外在環(huán)境基礎(chǔ).

Le等[31]在研究北京疫情期間重度灰霾發(fā)生的成因時(shí),發(fā)現(xiàn)在COVID-19疫情期間一次污染物排放量下降的情況下,PM2.5和臭氧含量增加的原因主要可能是二次氣溶膠和臭氧的復(fù)雜化學(xué)作用. PM2.5主要來(lái)源于其前體物通過(guò)復(fù)雜的化學(xué)反應(yīng)生成.其一,大氣中SO2經(jīng)過(guò)氣相氧化或液相氧化反應(yīng)生成硫酸鹽.其二,大氣中NO2被OH氧化生成HNO3,而HNO3與堿性物質(zhì)反應(yīng)生成硝酸鹽,或NO2與O3反應(yīng)產(chǎn)生的N2O5在顆粒物表面非均相水解生成硝酸鹽.其三,揮發(fā)性有機(jī)物(VOCs)能通過(guò)大氣光氧化過(guò)程、成核過(guò)程、凝結(jié)和氣粒分配過(guò)程及非均相反應(yīng)等化學(xué)過(guò)程生成二次氣溶膠[31-32].同理,O3污染主要是大氣環(huán)境中NO、CO、VOCs通過(guò)大氣光化反應(yīng)生成[33].上述過(guò)程中,生成的二次氣溶膠和O3進(jìn)入大氣系統(tǒng),這一過(guò)程為大氣復(fù)合系統(tǒng)演化至SOC狀態(tài),提供了緩慢的物質(zhì)或能量的輸入來(lái)源.在城市局域大氣環(huán)境中,當(dāng)O3和二次氣溶膠的前體物聚集時(shí),高氧化污染物O3能將大氣中SO2和NO通過(guò)復(fù)雜的化學(xué)反應(yīng)氧化成硫酸鹽和硝酸鹽.二次氣溶膠質(zhì)量濃度能通過(guò)增強(qiáng)對(duì)紫外光輻射的散射,吸收等能力減低大氣中生成O3的光化學(xué)反應(yīng)速率.在大氣復(fù)合污染系統(tǒng)中組元間短程近鄰相互作用下,O3和二次氣溶膠非線性演化將自發(fā)地向SOC狀態(tài)演化并維持著該臨界狀態(tài).當(dāng)大氣復(fù)合系統(tǒng)穩(wěn)定在某一臨界狀態(tài)時(shí),外界任何一點(diǎn)污染物的輸入都可能導(dǎo)致大氣復(fù)合系統(tǒng)發(fā)生鏈?zhǔn)椒磻?yīng),在城市局部空間內(nèi)形成較嚴(yán)重的PM2.5(或O3)污染事件,這現(xiàn)象類似于沙堆模型中的崩塌事件.此外,研究還發(fā)現(xiàn)長(zhǎng)株潭城市群3個(gè)城市PM2.5和O3小時(shí)平均濃度序列具有顯著的長(zhǎng)期持續(xù)性和冪律分形結(jié)構(gòu)特征,這兩類特征是系統(tǒng)存在SOC 特性的重要判斷依據(jù).盡管疫情期間人類活動(dòng)產(chǎn)生的污染物銳減,但這并沒(méi)有改變大氣污染系統(tǒng)內(nèi)在的非線性演化機(jī)制.因此,可以認(rèn)為大氣復(fù)合污染的內(nèi)在SOC內(nèi)稟機(jī)制是使得COVID-19疫情期間長(zhǎng)株潭城市群發(fā)生高污染的主導(dǎo)動(dòng)力學(xué)機(jī)制.

3 結(jié)論

3.1 長(zhǎng)株潭城市群3個(gè)城市PM2.5和O3日內(nèi)平均濃度規(guī)律呈現(xiàn)出顯著的相似性. 3個(gè)城市O3呈現(xiàn)出晝高夜低的單峰型,而PM2.5日變化規(guī)律表現(xiàn)出晝低夜高單峰型,與非疫情期間的特征有所差異.從實(shí)證研究結(jié)果看,O3晝高夜低可能由太陽(yáng)光輻射強(qiáng)度變化引起的,而PM2.5晝低夜高的主要原因可能是太陽(yáng)輻射和本身的沉降作用共同作用的.

3.2 長(zhǎng)株潭城市群3個(gè)城市PM2.5和O3小時(shí)平均濃度序列呈現(xiàn)出的長(zhǎng)期持續(xù)性特征,且PM2.5-O3互相關(guān)性也呈現(xiàn)出顯著的長(zhǎng)期持續(xù)性特征.該結(jié)果表明大氣復(fù)合系統(tǒng)中PM2.5和O3濃度演化并非相互獨(dú)立的,二者之間的相互作用也是不可忽視的.

3.3 長(zhǎng)株潭城市群3個(gè)城市PM2.5,O3小時(shí)平均濃度序列以及PM2.5-O3互相關(guān)性均存在著顯著的多重分形特征.這個(gè)結(jié)果進(jìn)一步驗(yàn)證了疫情期間大氣系統(tǒng)中PM2.5和O3濃度并非分別獨(dú)自演化,而是依然存在復(fù)雜的相互作用.

3.4 長(zhǎng)株潭城市群PM2.5和O3濃度演化可能不僅受到污染物本身的長(zhǎng)期持續(xù)性影響,同時(shí)還受到二者之間互相關(guān)的長(zhǎng)期持續(xù)性影響.此外,長(zhǎng)沙、株洲和湘潭3個(gè)城市PM2.5和O3濃度還呈現(xiàn)出顯著的冪律分形結(jié)構(gòu).長(zhǎng)期持續(xù)性和冪律分形結(jié)構(gòu)都是SOC內(nèi)稟機(jī)制的宏觀涌現(xiàn).因此,該研究認(rèn)為這種SOC內(nèi)稟機(jī)制可能是使得COVID-19疫情期間高污染事件發(fā)生的主導(dǎo)機(jī)制.

[1] 樂(lè) 旭,雷亞棟,周 浩,等.新冠肺炎疫情期間中國(guó)人為碳排放和大氣污染物的變化 [J]. 大氣科學(xué)學(xué)報(bào), 2020,43(2):265-274. Le X, Lei Y D, Zhou H, et al. Changes of anthropogenic carbon emissions and air pollutants during the COVID-19 epidemic in China [J]. Transactions of Atmospheric Sciences, 2020,43(2):265-274.

[2] Fan C, Li Y, Guang J, et al. The impact of the control measures during the COVID-19 outbreak on air pollution in China [J]. Remote Sensing, 2020,12(10):1613.

[3] Sicard P, Marco A D, Agathokleous E, et al. Amplified ozone pollution in cities during the COVID-19lockdown [J]. Science of The Total Environment, 2020,735:3382-3393.

[4] Huang X, Ding A J, Gao J, et al. Enhanced secondary pollution offset reduction of primary emissions during COVID-19 lockdown in China [J]. National Science Review, 2020,0:1-9.

[5] Chang Y H, Huang R J, Ge X L, et al. Puzzling haze events in China during the Coronavirus (COVID-19) Shutdown [J]. Geophysical Research Letters, 2020,47(12):e2020GL088533.

[6] 羅岳平,劉孟佳,甘 杰,等.長(zhǎng)株潭城市環(huán)境空氣中PM2.5和O3質(zhì)量濃度的相關(guān)性研究 [J]. 安全與環(huán)境學(xué)報(bào), 2015,15(4):313-317. Luo Y P, Liu M J, Gan J, et al. Correlation study on PM2.5and O3mass concentrations in ambient air by taking urban cluster of Changsha, Zhuzhou and Xiangtan as an example [J]. Journal of Safety and Environment, 2015,15(4):313-317.

[7] 賴安琪,陳曉陽(yáng),劉一鳴,等.珠江三角洲高質(zhì)量濃度PM2.5和O3復(fù)合污染特征 [J]. 中山大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018,57(4):30-36. Lai A Q, Chen X Y, Liu Y M, et al. Characteristics of complex pollution with high concentrations of PM2.5and O3over the Pearl River Delta, China [J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2018,57(4):30-36.

[8] Wang G H, Zhang R Y, Gomez M E, et al. Persistent sulfate formation from London Fog to Chinese haze [J]. Proceedings of the National Academy of sciences of the United States of America, 2016,113(48):13630-13635.

[9] Chen L, Liu R M, Huang Q, et al. An integrated simulation- monitoring framework for nitrogen assessment: A case study in the Baixi watershed, China [J]. Procedia Environmental Sciences, 2012, 13(4):1076-1090.

[10] Husain L, Rattigan O V, Dutkiewicz V, et al. Case studies of the SO2+ H2O2reaction in clouds [J]. Journal of Geophysical Research Atmospheres, 2000,105(D8):9831-9841.

[11] Sarwar G, Fahey K, Kwok R, et al. Potential impacts of two SO2oxidation pathways on regional sulfate concentrations: Aqueous-phase oxidation by NO2and gas-phase oxidation by Stabilized Criegee Intermediates [J]. Atmospheric Environment, 2013,68(4):186-197.

[12] Griffith S M, Huang X H H, Louie P K K, et al. Characterizing the thermodynamic and chemical composition factors controlling PM2.5nitrate: Insights gained from two years of online measurements in Hong Kong [J]. Atmospheric Environment, 2015,122(12):864-875.

[13] 李 紅,彭 良,畢 方,等.我國(guó)PM2.5與臭氧污染協(xié)同控制策略研究 [J]. 環(huán)境科學(xué)研究, 2019,32(10):1763-1778. Li H, Peng L, Bi F, et al. Strategy of coordinated control of PM2.5and Ozone in China [J]. Research of Environmental Sciences, 2019,32(10): 1763-1778.

[14] 張宇靜,趙天良,殷翀之,等.徐州市大氣PM2.5與O3作用關(guān)系的季節(jié)變化 [J]. 中國(guó)環(huán)境科學(xué), 2019,39(6):2267-2272. Zhang Y J, Zhao T L, Yin C Z, et al. Seasonal variation of the relationship between surface PM2.5and O3concentrations in Xuzhou [J]. China Environmental Science, 2019,39(6):2267-2272.

[15] 杜 娟,劉春瓊,黃紅良,等.成都冬季PM2.5演化的EEMD分解及自組織臨界態(tài) [J]. 環(huán)境科學(xué)與技術(shù), 2019,42(8):133-141.Du J, Liu C Q, Huang H L, et al. Multiscale EEMD decomposition of PM2.5in Chengdu during winter and self-organized criticality of air pollution [J]. Environmental Science & Technology, 2019,42(8):133-141.

[16] 胡晨霞,鄒 濱,李沈鑫,等.城市微環(huán)境PM2.5濃度空間分異特征分析 [J]. 中國(guó)環(huán)境科學(xué), 2018,38(3):910-916. Hu C X, Zou B, Li S X, et al. Spatial heterogeneity analysis of PM2.5concentrations in intra-urban microenvironments [J]. China Environmental Science, 2018,38(3):910-916.

[17] Zhou W X. Multifractal detrended cross-correlation analysis for two nonstationary signals [J]. Physical Review E Statal Nonlinear & Soft Matter Physics, 2008,77(6):066211.

[18] 史 凱.成都市一次重度灰霾期間大氣PM2.5的自組織臨界特性 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2014,34(10):2645-2653. Shi K. Self-organized criticality of PM2.5during a typical haze period in Chengdu, China [J]. Acta Scientiae Circumstantiae, 2014,34(10): 2645-2653.

[19] 史 凱,劉春瓊,吳生虎.大氣PM10跨界輸送的自組織動(dòng)力機(jī)制——以舟山市為例 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2014,34(5):1125-1132. Shi K, Liu C Q, Wu S H. Self-organized evolution of trans-boundary PM10pollution in Zhoushan City [J]. Acta Scientiae Circumstantiae, 2014,34(5):1125-1132.

[20] Bak P, Tang C, Wiesenfeld K. Self-organized criticality: an explanation of the noise [J]. Physics Review Letters, 1987,59(4):381-384.

[21] Turcotte D L, Malamud B D. Landslides, forest fires, and earthquakes: examples of self-organized critical behavior [J]. Statistical Mechanics and its Applications, 2004,340(4):580-589.

[22] 蘇志華,韓會(huì)慶,李 莉,等.貴陽(yáng)市臭氧的時(shí)空分布,氣象作用及其與前體物的關(guān)系 [J]. 中山大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020,59(5):102-112. Su Z H, Han H Q, Li L, et al. Temporal and spatial distribution of ozone, meteorological effect and relation with its precursors in Guiyang city [J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2020,59(5):102-112.

[23] 殷永泉,單文坡,紀(jì) 霞,等.濟(jì)南大氣臭氧濃度變化規(guī)律 [J]. 環(huán)境科學(xué), 2006,27(11):2299-2302.Yin Y Q, Shan W P, Ji X, et al.Characteristics of atmospheric ozone in the urban area of Ji'nan [J]. Environmental Science, 2006,27(11):2299-2302.

[24] 楊關(guān)盈,樊曙先,湯 潔,等.臨安近地面臭氧變化特征分析 [J]. 環(huán)境科學(xué)研究, 2008,21(3):31-35. Yang G Y, fan S X, Tang J, et al. Characteristic of surface ozone concentrations at Lin'an [J]. Research of Environmental Sciences, 2008,21(3):31-35.

[25] 郭含文,丁國(guó)棟,趙媛媛,等.城市不同綠地PM2.5質(zhì)量濃度日變化規(guī)律 [J]. 中國(guó)水土保持科學(xué), 2013,11(4):99-103. Guo H W, Ding G D, Zhao Y Y, et al.Diurnal variations in the mass concentration of suspended Particulate Matter 2.5(PM2.5) of different urban green space [J].Science of Soil and Water Conservation, 2013, 11(4):99-103.

[26] 張冬健,薛惠文.2015年北京,上海和拉薩PM2.5濃度日變化和日際變化分析 [J]. 北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018,54(4):705-712. Zhang D J, Xue H W. Interdiurnal and diurnal variation of PM2.5concentration in Beijing, Shanghai, and Lhasa in 2015 [J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2018,54(4):705-712.

[27] 于 洲,劉壽東,王詠薇,等.杭州市2014年城區(qū)大氣污染物濃度變化特征觀測(cè)分析 [J]. 科學(xué)技術(shù)與工程, 2016,16(16):95-104. Yu Z, Liu S D, Wang Y W, et al. Observation of atmospheric pollutants in the urban area of Hangzhou [J].Science Technology and Engineering, 2016,16(16):95-104.

[28] Shi K, Liu C Q. Self-organized criticality of air pollution [J]. Atmospheric Environment, 2009,43(21):3301-3304.

[29] 蘇 蓉,史 凱,黃正文,等.基于自組織臨界的上海PM10污染預(yù)警與風(fēng)險(xiǎn)評(píng)價(jià) [J]. 中國(guó)環(huán)境科學(xué), 2010,30(7):888-892. Su R, Shi K, Huang Z W, et al. Risk assessment and for warning of PM10pollution in Shanghai based on self-organized criticality theory [J]. China Environmental Science, 2010,30(7):888-892.

[30] Shi K, Liu C, Huang Y. Multifractal processes and self-organized criticality of PM2.5during a typical haze period in Chengdu, China [J]. Aerosol and Air Quality Research, 2015,15(3):926-934.

[31] Le T H, Wang Y, Liu L, et al. Unexpected air pollution with marked emission reductions during the COVID-19 outbreak in China [J]. Science, 2020,369(6504):702-706.

[32] 葉興南,陳建民.大氣二次細(xì)顆粒物形成機(jī)理的前沿研究 [J]. 化學(xué)進(jìn)展, 2009,21(Z1):288-296. Ye X G, Chen J M. Advances in the mechanism of secondary fine particulate matters formation [J]. Progress in Chemistry, 2009, 21(Z1):288-296.

[33] 賴安琪,陳曉陽(yáng),劉一鳴,等.珠江三角洲PM2.5和O3復(fù)合污染過(guò)程的數(shù)值模擬 [J]. 中國(guó)環(huán)境科學(xué), 2017,37(11):4022-4031. Lai A Q, Chen X Y, Liu Y M, et al. Numerical simulation of a complex pollution episode with high concentrations of PM2.5and O3over the Pearl River Delta region, China [J]. China Environmental Science, 2017,37(11):4022-4031.

Nonlinear dynamic mechanism of regional air pollution during the COVID-19 lockdown.

WU Bo1, LIU Chun-qiong2, ZHANG Jiao1, LI Yan-hui1, CHEN Yu-bing1,WEN Ye1, BAO Bing-yi1, DU Juan1, SHI Kai2*

(1.College of Mathematics and Statistics, Jishou University, Jishou 416000, China；2.College of Biology and Environmental Science, Jishou University, Jishou 416000, China)., 2021,41(5)：2028～2039

Based on hourly concentration of PM2.5and O3during the epidemic period (January 24, 2020 to May 31, 2020) in Changsha, Zhuzhou and Xiangtan, the diurnal patterns, long-term persistence, multifractality and self-organization evolution dynamics of these two pollutants were studied to reveal the internal dynamic mechanism of the occurrence and evolution of heavy pollution events during the epidemic period. Firstly, the diurnal patterns of PM2.5and O3concentrations were investigated. It showed that O3showed a single peak of high concentration in the daytime and low in the night, while PM2.5showed a single lowest peak concentration in the day and high in the night, which was different from the pattern in non-epidemic periods. Furthermore, detrended fluctuation analysis (DFA), the multifractal detrended fluctuation analysis (MFDFA) and probability statistical analysis were applied to study the long-term persistence, multi-fractal structure of PM2.5and O3series. The results showed that PM2.5and O3series had significant long-term persistence characteristics and strong multi-fractal structures for the three cities. Meanwhile, detrended cross-correlation analysis (DCCA) and multifractal detrended cross-correlation analysis (MFDCCA) were conducted to estimate the cross-correlations between PM2.5and O3series. Long-term persistence as well as multifractal features at different time scales was also observed in PM2.5-O3cross-correlations. Next, nonlinear analysis results obtained during epidemic period were compared with those obtained in the same periods of non-epidemic years of 2019 and 2018. Finally, based on the self-organized criticality (SOC) theory, the internal dynamic law of spatial and temporal evolution of PM2.5and O3series was discussed. Combined with the typical regional meteorological characteristics, it was found that the intrinsic dynamic mechanism of SOC may be one of the leading mechanisms of heavy air pollution episodes during the COVID-19 lockdown period. During the epidemic period, PM2.5and O3concentrations did not evolve independently but remained complex interactions. Under the stable meteorological conditions, the nonlinear coupling effect inside the air combined pollution might reach the dynamic critical state, thus, lead to the risk of heavy air pollution in Greater Changsha Metropolitan Region during the epidemic period.

diurnal patterns；long-term persistence characteristics；multifractal characteristics；power-law；self-organized criticality；PM2.5；O3

X513

A

1000-6923(2021)05-2028-12

吳 波(1995-),男,湖南吉首人,吉首大學(xué)碩士研究生,主要從事應(yīng)用統(tǒng)計(jì)分析.發(fā)表論文1篇.

2020-09-27

湖南省自然科學(xué)基金資助項(xiàng)目(2020JJ4504);吉首大學(xué)校級(jí)科研項(xiàng)目資助(Jdy20062,Jdx02)

* 責(zé)任作者, 教授, einboplure@163.com