• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2021-05-29 03:55:02劉春瓊李彥輝陳郁兵鮑冰逸
    中國環(huán)境科學 2021年5期
    關(guān)鍵詞:大氣特征疫情

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

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

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

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

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

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

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

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

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

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

    1 材料與方法

    1.1 材料

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

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

    1.2 研究方法

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    當?越大, 多重分形性越強.

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

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

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

    2 結(jié)果與分析

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

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

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

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

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

    圖3 疫情期間長沙、株洲和湘潭3個城市PM2.5和O3的DFA以及PM2.5-O3的DCCA圖

    為了定量分析PM2.5和O3之間相互作用對二者濃度演化過程中的影響,應(yīng)用去趨勢互相關(guān)分析(DCCA)方法對長株潭城市群3個城市的PM2.5和O3濃度相互作用進行分析,其結(jié)果如圖3所示.長沙、株洲和湘潭3個城市在研究時段內(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)系.其中長沙、株洲和湘潭DCCA標度指數(shù)分別為0.893、0.910、0.900.長株潭城市群3個城市的分析結(jié)果表明,在COVID-19疫情期間內(nèi)PM2.5和O3濃度的相關(guān)性均表現(xiàn)出較強的長期持續(xù)性特征.這意味著在過去的一段時間內(nèi)3個城市PM2.5(O3)污染物排放增加將導致未來一段時間內(nèi)3個城市O3(PM2.5)濃度持續(xù)增加.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    2.6 討論

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

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

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

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

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

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

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

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

    3 結(jié)論

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

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

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

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

    [1] 樂 旭,雷亞棟,周 浩,等.新冠肺炎疫情期間中國人為碳排放和大氣污染物的變化 [J]. 大氣科學學報, 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] 羅岳平,劉孟佳,甘 杰,等.長株潭城市環(huán)境空氣中PM2.5和O3質(zhì)量濃度的相關(guān)性研究 [J]. 安全與環(huán)境學報, 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] 賴安琪,陳曉陽,劉一鳴,等.珠江三角洲高質(zhì)量濃度PM2.5和O3復合污染特征 [J]. 中山大學學報(自然科學版), 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] 李 紅,彭 良,畢 方,等.我國PM2.5與臭氧污染協(xié)同控制策略研究 [J]. 環(huán)境科學研究, 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]. 中國環(huán)境科學, 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)境科學與技術(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]. 中國環(huán)境科學, 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)境科學學報, 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跨界輸送的自組織動力機制——以舟山市為例 [J]. 環(huán)境科學學報, 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] 蘇志華,韓會慶,李 莉,等.貴陽市臭氧的時空分布,氣象作用及其與前體物的關(guān)系 [J]. 中山大學學報(自然科學版), 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] 殷永泉,單文坡,紀 霞,等.濟南大氣臭氧濃度變化規(guī)律 [J]. 環(huán)境科學, 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)境科學研究, 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] 郭含文,丁國棟,趙媛媛,等.城市不同綠地PM2.5質(zhì)量濃度日變化規(guī)律 [J]. 中國水土保持科學, 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]. 北京大學學報(自然科學版), 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ū)大氣污染物濃度變化特征觀測分析 [J]. 科學技術(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污染預警與風險評價 [J]. 中國環(huán)境科學, 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] 葉興南,陳建民.大氣二次細顆粒物形成機理的前沿研究 [J]. 化學進展, 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] 賴安琪,陳曉陽,劉一鳴,等.珠江三角洲PM2.5和O3復合污染過程的數(shù)值模擬 [J]. 中國環(huán)境科學, 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-),男,湖南吉首人,吉首大學碩士研究生,主要從事應(yīng)用統(tǒng)計分析.發(fā)表論文1篇.

    2020-09-27

    湖南省自然科學基金資助項目(2020JJ4504);吉首大學校級科研項目資助(Jdy20062,Jdx02)

    * 責任作者, 教授, einboplure@163.com

    猜你喜歡
    大氣特征疫情
    大氣的呵護
    軍事文摘(2023年10期)2023-06-09 09:15:06
    戰(zhàn)疫情
    抗疫情 顯擔當
    疫情中的我
    37°女人(2020年5期)2020-05-11 05:58:52
    如何表達“特征”
    不忠誠的四個特征
    當代陜西(2019年10期)2019-06-03 10:12:04
    抓住特征巧觀察
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    寂寞人妻少妇视频99o| 丝袜美足系列| 91精品国产国语对白视频| 高清不卡的av网站| 国内精品宾馆在线| 日韩一区二区三区影片| 纯流量卡能插随身wifi吗| 久久亚洲国产成人精品v| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 热99国产精品久久久久久7| 久久久久久久大尺度免费视频| 国产精品不卡视频一区二区| 欧美日韩综合久久久久久| 久久精品久久久久久噜噜老黄| 久久午夜综合久久蜜桃| 人人妻人人添人人爽欧美一区卜| 亚洲国产欧美日韩在线播放| 日韩成人伦理影院| 黄片无遮挡物在线观看| 亚洲色图综合在线观看| 男人舔女人的私密视频| 老司机影院毛片| 观看av在线不卡| 国产综合精华液| 亚洲精品第二区| h视频一区二区三区| 日本-黄色视频高清免费观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品av麻豆狂野| 热re99久久精品国产66热6| av视频免费观看在线观看| 国产成人精品无人区| 日韩大片免费观看网站| 2021少妇久久久久久久久久久| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人 | videossex国产| 熟妇人妻不卡中文字幕| 色网站视频免费| 日韩制服骚丝袜av| 欧美 亚洲 国产 日韩一| 丁香六月天网| 亚洲人成77777在线视频| 欧美bdsm另类| 看免费av毛片| 欧美亚洲日本最大视频资源| 亚洲久久久国产精品| 亚洲激情五月婷婷啪啪| 国产成人免费观看mmmm| 精品久久久精品久久久| 少妇精品久久久久久久| 成年美女黄网站色视频大全免费| 日韩av免费高清视频| 女性被躁到高潮视频| 久久久久久久亚洲中文字幕| 美女大奶头黄色视频| 亚洲综合色惰| 色哟哟·www| 精品少妇久久久久久888优播| 成人无遮挡网站| 日韩欧美精品免费久久| 免费日韩欧美在线观看| 精品一区二区三区视频在线| 建设人人有责人人尽责人人享有的| 韩国高清视频一区二区三区| 国产精品蜜桃在线观看| 97在线人人人人妻| 久久精品熟女亚洲av麻豆精品| 亚洲人成网站在线观看播放| 久久鲁丝午夜福利片| 大陆偷拍与自拍| 极品人妻少妇av视频| 国产成人免费无遮挡视频| 另类精品久久| kizo精华| 人人妻人人澡人人看| 美女福利国产在线| 国产亚洲精品第一综合不卡 | 在线 av 中文字幕| 80岁老熟妇乱子伦牲交| 免费人妻精品一区二区三区视频| xxxhd国产人妻xxx| av在线app专区| 精品国产一区二区三区久久久樱花| 久久久久久伊人网av| 嫩草影院入口| 乱码一卡2卡4卡精品| 国产有黄有色有爽视频| 一级片'在线观看视频| av女优亚洲男人天堂| 日韩欧美一区视频在线观看| 国产精品偷伦视频观看了| 精品第一国产精品| 18禁观看日本| 国产熟女欧美一区二区| 九色成人免费人妻av| a 毛片基地| 考比视频在线观看| 中文字幕人妻熟女乱码| 高清欧美精品videossex| 国产成人精品一,二区| 纵有疾风起免费观看全集完整版| 中国国产av一级| 国产精品久久久久久精品古装| 十八禁网站网址无遮挡| 精品一区二区三区四区五区乱码 | 久久久久久久久久成人| 日韩伦理黄色片| 国产精品久久久久久av不卡| 国产色爽女视频免费观看| 男女免费视频国产| 丰满迷人的少妇在线观看| 国产成人欧美| 国产熟女欧美一区二区| 中国国产av一级| 黄色 视频免费看| 亚洲中文av在线| 色婷婷av一区二区三区视频| 51国产日韩欧美| 国产精品不卡视频一区二区| av福利片在线| 大香蕉97超碰在线| 人人妻人人爽人人添夜夜欢视频| 亚洲伊人久久精品综合| 中文乱码字字幕精品一区二区三区| 欧美变态另类bdsm刘玥| 成人黄色视频免费在线看| 日日爽夜夜爽网站| 国产又爽黄色视频| 国产精品无大码| 91精品伊人久久大香线蕉| 一个人免费看片子| 一区二区三区精品91| 中文字幕av电影在线播放| av线在线观看网站| 成年美女黄网站色视频大全免费| 天堂8中文在线网| 免费看光身美女| 免费观看a级毛片全部| 在线观看免费高清a一片| 老女人水多毛片| 天天操日日干夜夜撸| 国产日韩一区二区三区精品不卡| 欧美成人午夜免费资源| 天堂中文最新版在线下载| 色哟哟·www| 国产无遮挡羞羞视频在线观看| 免费看光身美女| av片东京热男人的天堂| 99热全是精品| 99国产综合亚洲精品| 一级a做视频免费观看| 99热这里只有是精品在线观看| 伦理电影大哥的女人| 黑人高潮一二区| 少妇熟女欧美另类| 中文字幕人妻丝袜制服| av天堂久久9| 成年av动漫网址| 狠狠婷婷综合久久久久久88av| 国产精品三级大全| 精品99又大又爽又粗少妇毛片| 国产熟女欧美一区二区| 伦理电影免费视频| 久久精品国产鲁丝片午夜精品| 欧美xxⅹ黑人| 日韩成人av中文字幕在线观看| 亚洲精品一区蜜桃| 天天躁夜夜躁狠狠久久av| xxxhd国产人妻xxx| 日韩欧美一区视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产色片| 99国产精品免费福利视频| 久久精品久久久久久噜噜老黄| 在现免费观看毛片| 91精品三级在线观看| 一级毛片 在线播放| 51国产日韩欧美| 黄色一级大片看看| 多毛熟女@视频| 中文字幕亚洲精品专区| 天天躁夜夜躁狠狠久久av| xxxhd国产人妻xxx| 国产色婷婷99| 精品国产乱码久久久久久小说| 国内精品宾馆在线| 伦理电影大哥的女人| 国产国拍精品亚洲av在线观看| 国产亚洲一区二区精品| 制服丝袜香蕉在线| 91精品伊人久久大香线蕉| 欧美97在线视频| 国产黄频视频在线观看| 高清不卡的av网站| 秋霞伦理黄片| 国产福利在线免费观看视频| 国产成人91sexporn| 久久久精品94久久精品| 五月天丁香电影| 亚洲精品一区蜜桃| 少妇高潮的动态图| 91国产中文字幕| 欧美国产精品一级二级三级| 最后的刺客免费高清国语| 晚上一个人看的免费电影| 免费在线观看完整版高清| 边亲边吃奶的免费视频| 国产成人精品婷婷| 精品久久久精品久久久| 欧美3d第一页| 国产不卡av网站在线观看| 久久免费观看电影| xxx大片免费视频| 午夜免费鲁丝| 色哟哟·www| 大香蕉97超碰在线| 巨乳人妻的诱惑在线观看| 一级毛片黄色毛片免费观看视频| 日日撸夜夜添| 又黄又粗又硬又大视频| 丁香六月天网| 纵有疾风起免费观看全集完整版| √禁漫天堂资源中文www| 又黄又粗又硬又大视频| 国产爽快片一区二区三区| 国产伦理片在线播放av一区| 日韩av不卡免费在线播放| 国产精品三级大全| 亚洲国产精品成人久久小说| 男人添女人高潮全过程视频| 在线看a的网站| 国产亚洲欧美精品永久| 日本欧美国产在线视频| 国产成人精品一,二区| 免费观看在线日韩| 少妇的丰满在线观看| 好男人视频免费观看在线| 熟妇人妻不卡中文字幕| videosex国产| a级毛片在线看网站| 99re6热这里在线精品视频| 最新的欧美精品一区二区| 97在线视频观看| 建设人人有责人人尽责人人享有的| 成年人免费黄色播放视频| 免费在线观看黄色视频的| 美国免费a级毛片| 国产xxxxx性猛交| 婷婷色麻豆天堂久久| 黑人欧美特级aaaaaa片| 黄色 视频免费看| 国产探花极品一区二区| 欧美亚洲 丝袜 人妻 在线| www.熟女人妻精品国产 | 欧美日本中文国产一区发布| 男人操女人黄网站| 亚洲欧美日韩卡通动漫| 纯流量卡能插随身wifi吗| 国产成人精品福利久久| 日本午夜av视频| 亚洲色图综合在线观看| 高清毛片免费看| 久久久久久久久久人人人人人人| 色婷婷av一区二区三区视频| 中文字幕精品免费在线观看视频 | 大片电影免费在线观看免费| 精品亚洲成国产av| 精品一区二区三区四区五区乱码 | 18禁国产床啪视频网站| 国产免费一级a男人的天堂| 亚洲精品美女久久久久99蜜臀 | 黄色毛片三级朝国网站| 十八禁高潮呻吟视频| 国产精品 国内视频| 亚洲色图综合在线观看| 久久久久久久亚洲中文字幕| 日本av免费视频播放| 国语对白做爰xxxⅹ性视频网站| 十八禁高潮呻吟视频| 亚洲av国产av综合av卡| 卡戴珊不雅视频在线播放| 男人添女人高潮全过程视频| 亚洲熟女精品中文字幕| 亚洲色图综合在线观看| 国产成人精品一,二区| 国产黄色免费在线视频| 一二三四中文在线观看免费高清| 天天躁夜夜躁狠狠躁躁| 亚洲中文av在线| 观看av在线不卡| 丝袜美足系列| 日本91视频免费播放| 人成视频在线观看免费观看| 久久ye,这里只有精品| 韩国高清视频一区二区三区| 岛国毛片在线播放| 哪个播放器可以免费观看大片| 国产精品久久久久久av不卡| 国产亚洲av片在线观看秒播厂| 亚洲av日韩在线播放| 国产精品三级大全| 黄网站色视频无遮挡免费观看| 精品人妻熟女毛片av久久网站| 亚洲美女视频黄频| 99久国产av精品国产电影| 美国免费a级毛片| 好男人视频免费观看在线| 乱码一卡2卡4卡精品| 久久人妻熟女aⅴ| 又粗又硬又长又爽又黄的视频| 97人妻天天添夜夜摸| 成人午夜精彩视频在线观看| 国产永久视频网站| 成人影院久久| 久久久久久人人人人人| 色网站视频免费| xxx大片免费视频| 三级国产精品片| 岛国毛片在线播放| 爱豆传媒免费全集在线观看| 亚洲国产色片| 伦精品一区二区三区| 亚洲av成人精品一二三区| 最后的刺客免费高清国语| 观看美女的网站| 下体分泌物呈黄色| 亚洲精品一二三| 91国产中文字幕| 成人18禁高潮啪啪吃奶动态图| 久久久久久人人人人人| av.在线天堂| 久久久久人妻精品一区果冻| 三上悠亚av全集在线观看| 99视频精品全部免费 在线| 精品午夜福利在线看| 国产 精品1| 国产探花极品一区二区| 色婷婷av一区二区三区视频| 美女国产高潮福利片在线看| 18禁国产床啪视频网站| 国产麻豆69| 国产老妇伦熟女老妇高清| 人妻一区二区av| 欧美少妇被猛烈插入视频| 日韩视频在线欧美| 一级片'在线观看视频| 在线观看一区二区三区激情| 精品久久蜜臀av无| 国产精品久久久av美女十八| 精品久久蜜臀av无| 久久国产精品男人的天堂亚洲 | 99久久人妻综合| 国产一区二区在线观看日韩| 男人舔女人的私密视频| 菩萨蛮人人尽说江南好唐韦庄| 男人舔女人的私密视频| 精品一区二区免费观看| 乱码一卡2卡4卡精品| 七月丁香在线播放| 日本91视频免费播放| 十分钟在线观看高清视频www| 日日啪夜夜爽| av免费在线看不卡| 日本91视频免费播放| 成人国语在线视频| 久久久久久久亚洲中文字幕| 2018国产大陆天天弄谢| 久久久精品免费免费高清| 男女边吃奶边做爰视频| 18+在线观看网站| av天堂久久9| 亚洲精华国产精华液的使用体验| 成人18禁高潮啪啪吃奶动态图| 日韩在线高清观看一区二区三区| 天堂8中文在线网| 色5月婷婷丁香| 伦理电影免费视频| 伦精品一区二区三区| 天堂8中文在线网| 99国产精品免费福利视频| 日本av手机在线免费观看| 成年人免费黄色播放视频| 一区二区三区四区激情视频| 成人免费观看视频高清| 人妻人人澡人人爽人人| 国产免费一级a男人的天堂| 美女国产高潮福利片在线看| 国产黄色免费在线视频| 999精品在线视频| 午夜av观看不卡| 国产成人欧美| 国产一区二区在线观看av| 新久久久久国产一级毛片| 亚洲一级一片aⅴ在线观看| 亚洲国产毛片av蜜桃av| 日韩三级伦理在线观看| 亚洲精品久久成人aⅴ小说| 一级毛片电影观看| 亚洲久久久国产精品| 亚洲 欧美一区二区三区| 久久久久国产网址| 咕卡用的链子| 国产女主播在线喷水免费视频网站| 国产片特级美女逼逼视频| 边亲边吃奶的免费视频| 国产乱来视频区| 欧美精品一区二区免费开放| 久久人人97超碰香蕉20202| 国产成人91sexporn| 日本91视频免费播放| 国产成人精品婷婷| 亚洲精品一二三| 视频中文字幕在线观看| 少妇高潮的动态图| 夜夜骑夜夜射夜夜干| 国产福利在线免费观看视频| 国产日韩欧美视频二区| 欧美人与性动交α欧美精品济南到 | 一本色道久久久久久精品综合| 少妇精品久久久久久久| 欧美国产精品一级二级三级| 亚洲,欧美精品.| 亚洲国产精品一区三区| 女人久久www免费人成看片| 成人无遮挡网站| 国产免费福利视频在线观看| 欧美xxxx性猛交bbbb| 大香蕉久久成人网| 高清毛片免费看| 国产精品国产三级国产av玫瑰| 男人添女人高潮全过程视频| 丰满迷人的少妇在线观看| 成年人免费黄色播放视频| 丝袜美足系列| 男女下面插进去视频免费观看 | 久久久久精品久久久久真实原创| 国产日韩欧美亚洲二区| 大片电影免费在线观看免费| 久久av网站| 中文字幕av电影在线播放| 91aial.com中文字幕在线观看| 女的被弄到高潮叫床怎么办| av免费在线看不卡| 国产av一区二区精品久久| 在线观看美女被高潮喷水网站| 国产在线免费精品| 欧美 亚洲 国产 日韩一| a级毛片在线看网站| 另类精品久久| 99热网站在线观看| 男女无遮挡免费网站观看| 在线观看国产h片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲,一卡二卡三卡| 亚洲成色77777| 国产精品熟女久久久久浪| 亚洲一区二区三区欧美精品| 久久久久久人人人人人| 亚洲精品色激情综合| 丰满迷人的少妇在线观看| 久久精品国产亚洲av涩爱| 亚洲欧美日韩卡通动漫| 午夜福利在线观看免费完整高清在| 亚洲av免费高清在线观看| 亚洲精品自拍成人| 亚洲国产成人一精品久久久| 黑人高潮一二区| 青春草视频在线免费观看| 久久久久视频综合| 激情视频va一区二区三区| 国产成人欧美| 1024视频免费在线观看| 国产片特级美女逼逼视频| 欧美日韩综合久久久久久| 久久精品久久精品一区二区三区| 欧美变态另类bdsm刘玥| 免费久久久久久久精品成人欧美视频 | 国产精品久久久久久精品电影小说| 国产有黄有色有爽视频| 99九九在线精品视频| videosex国产| 在线观看免费高清a一片| www.av在线官网国产| xxxhd国产人妻xxx| 在线精品无人区一区二区三| 青春草亚洲视频在线观看| 国产成人午夜福利电影在线观看| 欧美另类一区| 80岁老熟妇乱子伦牲交| √禁漫天堂资源中文www| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美精品永久| av在线老鸭窝| 美女主播在线视频| 丰满饥渴人妻一区二区三| 国产黄频视频在线观看| 国产精品国产三级国产专区5o| 国产精品偷伦视频观看了| 啦啦啦啦在线视频资源| 18在线观看网站| 九色亚洲精品在线播放| 成人毛片60女人毛片免费| 国产精品嫩草影院av在线观看| 日本黄大片高清| 国产精品偷伦视频观看了| 国内精品宾馆在线| 欧美xxⅹ黑人| 黄色毛片三级朝国网站| 国产av精品麻豆| 日韩一本色道免费dvd| 欧美日韩av久久| 亚洲激情五月婷婷啪啪| 日韩一区二区三区影片| 蜜臀久久99精品久久宅男| 婷婷成人精品国产| 久久久久精品久久久久真实原创| 国产亚洲精品第一综合不卡 | 久久人人97超碰香蕉20202| 最新中文字幕久久久久| 久久精品国产亚洲av天美| 全区人妻精品视频| 国语对白做爰xxxⅹ性视频网站| 熟女电影av网| 99九九在线精品视频| 亚洲成人av在线免费| 一本久久精品| 成年美女黄网站色视频大全免费| 免费日韩欧美在线观看| 男的添女的下面高潮视频| 下体分泌物呈黄色| 亚洲av欧美aⅴ国产| 欧美精品亚洲一区二区| 国产精品久久久久久久电影| 国产精品久久久久成人av| 久久99一区二区三区| 老司机影院成人| 天天躁夜夜躁狠狠久久av| 秋霞在线观看毛片| av在线观看视频网站免费| 黄网站色视频无遮挡免费观看| 国产精品偷伦视频观看了| 蜜桃国产av成人99| 精品国产乱码久久久久久小说| 韩国精品一区二区三区 | 黑人巨大精品欧美一区二区蜜桃 | 精品亚洲成国产av| 亚洲av福利一区| 国产av一区二区精品久久| 亚洲av电影在线观看一区二区三区| 精品少妇黑人巨大在线播放| 卡戴珊不雅视频在线播放| 99精国产麻豆久久婷婷| 夜夜骑夜夜射夜夜干| 国产精品国产三级专区第一集| 男女午夜视频在线观看 | videosex国产| 久久人人爽人人片av| 我的女老师完整版在线观看| 欧美人与善性xxx| 老司机亚洲免费影院| 中文欧美无线码| 国产又爽黄色视频| 啦啦啦视频在线资源免费观看| 国产女主播在线喷水免费视频网站| 中文字幕免费在线视频6| 一级毛片电影观看| 美女脱内裤让男人舔精品视频| 久久久精品区二区三区| 精品久久蜜臀av无| 国产在线免费精品| 亚洲五月色婷婷综合| 日韩大片免费观看网站| 捣出白浆h1v1| 水蜜桃什么品种好| 日日撸夜夜添| 十八禁网站网址无遮挡| 人妻少妇偷人精品九色| www.色视频.com| 国产成人一区二区在线| 国产片特级美女逼逼视频| 亚洲欧美成人综合另类久久久| 国产不卡av网站在线观看| 欧美人与性动交α欧美软件 | 丝袜喷水一区| av在线播放精品| 男女啪啪激烈高潮av片| 丁香六月天网| 久久青草综合色| 日韩一本色道免费dvd| 丁香六月天网| av在线播放精品| 视频在线观看一区二区三区| 免费黄频网站在线观看国产| 久久精品久久精品一区二区三区| 日韩av不卡免费在线播放| 在线观看人妻少妇| 亚洲,一卡二卡三卡| 亚洲av中文av极速乱| 国产高清三级在线| 久久热在线av| 免费高清在线观看日韩| 欧美日韩av久久| 欧美精品一区二区免费开放| 男人舔女人的私密视频| 99热国产这里只有精品6| 少妇精品久久久久久久| 久久影院123| 99视频精品全部免费 在线| 99热6这里只有精品| 人人妻人人澡人人看| 亚洲国产看品久久| 视频中文字幕在线观看| 国产有黄有色有爽视频| 青春草国产在线视频|