韓力慧,蘭 童,程水源,王慎澳,田 健,齊超楠,肖 茜,王海燕,韓登越,王迎澳(.北京工業(yè)大學(xué)環(huán)境與生命學(xué)部,北京 0024；2.區(qū)域大氣復(fù)合污染防治北京市重點(diǎn)實(shí)驗(yàn)室,北京 0024)

唐山市是一座典型的重工業(yè)城市,顆粒物污染現(xiàn)象常在春、冬季頻繁發(fā)生,這可能與燃煤、不利的氣象條件、氣象傳輸?shù)纫蛩赜嘘P(guān).近年來隨著顆粒物治理工作的不斷開展,唐山市PM2.5和PM10濃度逐年下降,而O3問題逐漸凸顯,顆粒物和O3復(fù)合污染的治理已經(jīng)成為改善當(dāng)?shù)乜諝赓|(zhì)量的關(guān)鍵.

顆粒物和O3濃度存在多個(gè)時(shí)間尺度的波動(dòng),短期排放以及氣候的變化均會(huì)導(dǎo)致污染物濃度的變化.因此,將污染物序列進(jìn)行多時(shí)間尺度的劃分可以更加深入的探究污染物的變化規(guī)律,識(shí)別造成其濃度波動(dòng)的原因.KZ 濾波法通過濾波器對(duì)污染物原始序列進(jìn)行劃分,分離出原始時(shí)間序列中的高頻和低頻信號(hào)[1],將KZ 濾波結(jié)合逐步回歸模型可以去除氣象因素的干擾[2-3],識(shí)別源排放及氣象因素對(duì)于污染物濃度變化的貢獻(xiàn).KZ 濾波法可用于評(píng)估O3的控制效果[4],及顆粒物的研究[5].

污染物之間或污染物與氣象要素間的相互作用通常用相關(guān)性來表示[6-8],然而污染物間或污染物與氣象要素間的相關(guān)性與選取的時(shí)間尺度有關(guān),隨著時(shí)間尺度的增大相關(guān)系數(shù)也會(huì)變化.不同時(shí)間尺度下,顆粒物和O3間存在不同的相互作用,兩者間不僅存在共同的前體物,而且在各種轉(zhuǎn)化途徑上相互影響[9-10],同時(shí),氣象要素也會(huì)影響顆粒物與O3的生成過程[11-12].小波相干性分析的方法可以揭示不同時(shí)間尺度下污染物之間或污染物與氣象要素間的相干性[13].小波相干性分析已經(jīng)用于探討武漢、山西、長三角等地PM2.5、PM10、O3受其他污染物(CO、NO2、SO2等)或氣象要素中溫度(T)、相對(duì)濕度(RH)、地表壓力(P)等的影響,其研究結(jié)論也證明了小時(shí)間尺度和大時(shí)間尺度下彼此間存在著不同的作用機(jī)制[14-15].

因此,本文基于唐山市2015～2022 年的PM2.5、PM10以及O3濃度數(shù)據(jù),使用KZ 濾波法探究顆粒物和O3的多時(shí)間尺度變化特征,同時(shí)將KZ 濾波法與多元逐步回歸法相結(jié)合,定量識(shí)別氣象因素以及源排放作用對(duì)于PM2.5、PM10以及O3濃度變化的影響,最后使用小波相干性分析法探究PM2.5、PM10和O3間多時(shí)間尺度的相互作用機(jī)制以及氣象因素對(duì)于各污染物生成的影響,以期為唐山市顆粒物和O3污染防治工作提供參考.

1 材料與方法

1.1 數(shù)據(jù)來源

通過空氣質(zhì)量在線監(jiān)測平臺(tái)(https://www.aqistudy.cn/historydata/)下載2015 年1 月1 日～ 2022年12 月31 日,唐山市PM2.5、PM10日均濃度以及O38h 日均濃度數(shù)據(jù).KZ 濾波分析和多元線性回歸擬合時(shí),所需氣象數(shù)據(jù)均來源于歐洲中期天氣預(yù)報(bào)中心(ECMWF)的再分析資料ERA- Interim 中的地面氣象要素參數(shù),包括每日的平均溫度T( ℃)、相對(duì)濕度RH(%)、氣壓P(kPa)、10m 水平風(fēng)速U(m/s)、10m 垂直風(fēng)速V(m/s)、太陽短波輻射SR(J/m2)、太陽長波輻射LR(J/m2)以及單日累積降水量Pre(mm).

1.2 KZ 濾波法

KZ 濾波法是一種時(shí)間序列分析的方法,它通過滑動(dòng)平均的原理對(duì)原始序列進(jìn)行濾波.KZ(m,p)濾波器對(duì)原始序列的處理過程如式(1)所示:

式中: X 為污染物濃度原始時(shí)間序列; Xi表示原始時(shí)間序列中第i 日污染物濃度; j 為滑動(dòng)窗口變量; k 為對(duì)原始序列X 進(jìn)行濾波時(shí)的單邊滑動(dòng)窗口長度; m為滑動(dòng)窗口長度,且m=2k+1.Yi為原始序列經(jīng)過一次濾波后第i 日污染物濃度.使用KZ(m,p)濾波器進(jìn)行濾波時(shí),需要通過公式(1)進(jìn)行多次濾波,其中將原始序列第一次濾波后輸出的序列作為第二次濾波的輸入序列,進(jìn)行第二次濾波的過程稱為1 次迭代,重復(fù)進(jìn)行p 次迭代后,即可獲得原始序列的不同尺度分量.m 和p 的取值與劃分的時(shí)間尺度有關(guān),例如,對(duì)于長期分量,通常m 取365d,p 為3 次.當(dāng)滿足式(2)時(shí),可將天數(shù)小于等于N 的高頻濃度濾除[16]:

在濾波過程中,假定污染物的原始濃度序列可被劃分為短期分量、季節(jié)分量以及長期分量,即:

式中:X(t)為原始時(shí)間序列;E(t)、S(t)和W(t)依次為長期分量、季節(jié)分量和短期分量時(shí)間序列.

長期分量E(t)可通過m 為365d,p 為3 次的濾波器KZ(365,3)對(duì)原始序列進(jìn)行濾波得到,如公式(4)所示.為了獲得季節(jié)分量和短期分量,引入了基線分量XB(t),定義為長期分量與季節(jié)分量之和.該基線分量可通過m 為15d,p 為5 次的濾波器KZ(15,5)對(duì)原始序列進(jìn)行濾波得到.因此,季節(jié)分量S(t)可通過基線分量XB(t)與長期分量E(t)的差值來計(jì)算,如公式(6)所示,短期分量可由原始序列與基線分量的差值計(jì)算得到,如公式(7)所示.由此可見,使用KZ(365,3)和KZ(15,5)對(duì)原始序列進(jìn)行濾波后,可分別根據(jù)以下公式計(jì)算長期分量、季節(jié)分量和短期分量[17].

當(dāng)長期分量、季節(jié)分量以及短期分量3 個(gè)分量之間的關(guān)系相互獨(dú)立時(shí),其方差關(guān)系應(yīng)滿足:

將長期分量、季節(jié)分量以及短期分量的方差和與原始序列進(jìn)行比較,并將兩者的比值(即可解釋方差)作為判斷原始序列劃分可靠性的主要依據(jù).

1.3 多元逐步回歸

KZ 濾波后污染物的各分量均包含了氣象因素和源排放因素的影響,針對(duì)氣象因素,采用多元逐步回歸的方法[18],即以污染物濃度作為因變量,各氣象因素作為自變量進(jìn)行逐步回歸,剔除對(duì)于污染物濃度影響較弱的氣象因子,得到基于多個(gè)重要?dú)庀笠蛩氐奈廴疚餄舛茸顑?yōu)擬合模型.此時(shí),該污染物濃度分量與其基于氣象因素的最優(yōu)擬合模型的殘差可近似看作主要受源排放影響的部分,氣象因素影響極小,可以忽略,如公式(9)和(10)所示.本研究主要探究污染物的長期分量中氣象因素和源排放因素的影響,氣象因素的剝離過程如下:

式(9)和(10)中:ωi(t)分別為短期回歸分量以及基線回歸分量中的第i 個(gè)氣象因子; ai、bi為第i 個(gè)氣象因子的偏回歸系數(shù); a0、b0為回歸方程的截距項(xiàng),多元逐步回歸結(jié)果為各污染物在95%置信水平上的統(tǒng)計(jì)模型.εw(t)和εB(t)分別是短期回歸序列和基線回歸序列的殘差序列,式(11)中兩者相加即為總殘差序列ε(t)使用KZ(365,3)濾波后可獲得總殘差序列的長期分量εLT(t)式(12)中將εLT(t)與長期序列的均值相加即為調(diào)整后的長期序列XLT(t),該序列僅與源排放作用有關(guān).本研究使用逐步回歸法進(jìn)行分析時(shí),初始采用的氣象因子包括T、RH、P、U、V、SR、LR 以及Pre.

1.4 小波相干性分析

KZ 濾波法與多元逐步回歸法相結(jié)合的最大特點(diǎn)在于可以分別定量估算源排放和氣象因素對(duì)污染物長期分量的貢獻(xiàn),同時(shí),也可以定性分析不同時(shí)間尺度、污染物與氣象因素之間的關(guān)系,但在時(shí)間尺度上有一定的局限性.小波相干性分析以1d 為時(shí)間單位,可以實(shí)現(xiàn)在任何時(shí)間尺度上分析污染物與氣象因素之間的影響關(guān)系,還能揭示污染物與污染物之間的相互影響關(guān)系,相比于多元逐步回歸的方法,凸顯出較全面的優(yōu)勢.因此,本研究在KZ 濾波與多元逐步回歸分析的基礎(chǔ)上,結(jié)合小波相干性分析進(jìn)一步探究不同時(shí)間尺度上污染物PM2.5、PM10及O3之間,污染物與氣象因素之間的相互作用機(jī)制.

小波變換繼承和發(fā)展了傅里葉局部化變換的思想,克服了窗口大小不隨頻率變化等缺點(diǎn),是一種進(jìn)行信號(hào)時(shí)頻分析和處理的理想工具.污染物的濃度序列實(shí)質(zhì)上也是一種信號(hào),因此,可用小波變換來探究不同時(shí)間尺度污染物濃度序列特征.小波變換如公式(13)所示:

式中: WX(a,τ)為小波系數(shù); X(t)為原始時(shí)間序列; ψ(t)為基小波;a 為時(shí)間尺度因子;τ 為平移因子;t 為時(shí)間.

小波相干系數(shù)是將兩個(gè)相同維度的時(shí)間序列進(jìn)行小波變換后在頻域內(nèi)的相關(guān)系數(shù),它能清晰的反映兩個(gè)等長時(shí)間序列的數(shù)據(jù)在時(shí)間和頻率上的線性相關(guān)或非線性相關(guān)性.小波相干系數(shù)定義如下[19]:

式中:R(a,τ)為小波相干系數(shù);a 表示時(shí)間尺度;τ 代表平移因子,同式(13)所示;S 表示平滑算子;WX(a,τ)和WY(a,τ)分別表示X 和Y 的小波系數(shù); WXY(a,τ)是X 序列和Y 序列的交叉小波譜,其計(jì)算公式如下:

式中: W※Y(a,τ)表示W(wǎng)Y(a,τ)的復(fù)共軛.通過改變a 和τ的大小可以做出不同時(shí)間尺度的相干性譜圖.譜圖中灰白色區(qū)域以上的點(diǎn)通過了90%的可信度檢驗(yàn).此外,譜圖中箭頭與x 軸正方向的夾角β 反映兩污染物的相關(guān)性及相位差:若β∈(-π/2,π/2)則兩污染物之間正相關(guān),若β∈(π/2,3π/2)則兩污染物間負(fù)相關(guān);且當(dāng)β=π/2 時(shí),前者比后者提前變化1/4 個(gè)周期,當(dāng)β=π時(shí),前者比后者提前變化1/2 個(gè)周期,其余相位關(guān)系可類比計(jì)算.

2 結(jié)果與討論

2.1 污染物不同時(shí)間尺度的變化特征

2015～2022 年期間,PM2.5、PM10濃度原始序列均呈下降的趨勢,如圖1 所示,其中PM2.5從2015 年的84.17μg/m3下降到2022 年的36.88μg/m3, PM10從 2015 年的 140.53μg/m3下降到 2022 年的69.52μg/m3,兩者降幅分別為56%和50%.而O3濃度原始序列呈周期性波動(dòng),平均濃度為97.64μg/m3,年際整體變化幅度小于10μg/m3.

圖1 2015～2022 年期間唐山市PM2.5,PM10 和O3 不同時(shí)間尺度濃度變化Fig.1 Concentration variations of PM2.5, PM10 and O3 at different time scales in Tangshan City during 2015～2022

污染物原始序列經(jīng)過KZ 濾波后,PM2.5和PM10長期分量序列與原始序列變化趨勢一致,均呈現(xiàn)顯著下降趨勢,這可能與2013 年9 月《大氣污染防治行動(dòng)計(jì)劃》頒布以來,唐山市顆粒物的治理效果逐年提高有關(guān).例如,2021 年唐山市政府對(duì)鋼鐵、焦化、水泥等12 個(gè)行業(yè)企業(yè)實(shí)施月度總量減排,累計(jì)減少顆粒物4.40t、NOx2.67 萬t、CO 96.13 萬t,推進(jìn)PM2.5與O3、NOx與VOCs 協(xié)同控制[20].2022 年唐山市政府完成204 項(xiàng)超低排放改造工程和175 項(xiàng)CO、125 項(xiàng)NOx、335項(xiàng)VOCs治理工程,減排各項(xiàng)污染物101.12萬t[21],這些政策的施行使大氣污染物濃度大大降低.說明長期分量的變化可能主要受污染物減排的影響,這與孟曉艷等[22]的研究結(jié)果相一致.PM2.5和PM10的季節(jié)分量和短期分量均在每年的12 月～第二年2 月呈現(xiàn)強(qiáng)烈的波動(dòng),如圖1 中PM2.5和PM10的季節(jié)分量和短期分量所示,這可能與冬季供暖期大量燃煤排放致使唐山市冬季大氣顆粒物的濃度顯著升高,且高于其他季節(jié),同時(shí),冬季不利的氣象條件如逆溫層的頻繁出現(xiàn),以及對(duì)流層較低導(dǎo)致顆粒物濃度驟升[23]有很大關(guān)系,說明顆粒物濃度的季節(jié)分量和短期分量受污染物的排放和氣象因素的雙重影響較大.

O3的長期分量序列不同于原始序列的變化,其變化幅度較小,為11μg/m3(為了清晰顯示其變化特征,在其長期分量圖中將縱坐標(biāo)的單位刻度值選為3μg/m3,遠(yuǎn)小于原始序列中縱坐標(biāo)的單位刻度值75μg/m3),遠(yuǎn)小于原始序列的變化幅度277μg/m3,其分別在2015、2017 和2022 年出現(xiàn)峰值,自2017 年出現(xiàn)上升后,從2018～2021 年逐年緩步下降,但2022年濃度再次出現(xiàn)反彈,這可能與2022 年5～10 月,降水和云量普遍偏少、而日照較強(qiáng)[24]有利于O3的生成有關(guān).O3季節(jié)分量和短期分量每年在5～7 月間出現(xiàn)強(qiáng)烈波動(dòng),如圖1 所示,季節(jié)分量和短期分量濃度與其前體物VOCx和NOx的變化有關(guān)[25],而溫度和紫外線強(qiáng)度間接影響了O3的生成過程[26].由于唐山市日均溫度和太陽短波輻射強(qiáng)度在每年均呈現(xiàn)周期性變化,如圖2所示,春夏季呈上升趨勢,秋冬季呈下降趨勢,因此,O3濃度在每年中也呈現(xiàn)周期性的變化,即春夏季呈上升趨勢,這可能與春夏季O3前體物VOCx和NOx在較強(qiáng)太陽光作用下促進(jìn)了O3的生成;而秋冬季一方面由于顆粒物濃度的升高影響了太陽輻射的透過率,另一方面較低的溫度和光強(qiáng)不利于O3的生成,因此每年秋冬時(shí)段O3濃度呈下降趨勢.

圖2 唐山市2015～2022 年氣象要素T、RH、WS、V、Pre 和SR 的變化趨勢Fig.2 Variation trends of meteorological elements as T, RH, WS, V, Pre and SR in Tangshan City from 2015 to 2022

圖3 污染物的長期序列及經(jīng)氣象調(diào)整后的長期序列Fig.3 Long-term component and the adjusted long-term component of pollutants

為了識(shí)別影響污染物濃度變化的主要因素,本研究針對(duì)污染物原始序列的3 個(gè)分量進(jìn)行了方差分析,其占總方差的百分?jǐn)?shù)如表1 所示.這3 種污染物的長期、季節(jié)以及短期分量的方差占總方差的百分?jǐn)?shù)之和均接近100%,說明這3 個(gè)分量可以很好的解釋原始序列.PM2.5和PM10的濃度分量中短期分量占主要地位,其占比均超過65%,表明不利的短期氣象條件(高相對(duì)濕度、低風(fēng)速、逆溫層等)[27]及污染物短期排放(如早、晚高峰)是顆粒物濃度變化的主因.O3的濃度分量中季節(jié)分量的占比最大,氣象條件中溫度和光照的季節(jié)性變化對(duì)O3的生成影響較大.

表1 2015～2022 年唐山市顆粒物和O3各分量方差的貢獻(xiàn)(%)Table 1 Contribution of variances of particulate matter and O3 components in Tangshan City from 2015 to 2022(%)

2.2 氣象因素和源排放對(duì)污染物的貢獻(xiàn)

使用多元逐步回歸法對(duì)顆粒物和O3的長期分量序列進(jìn)行調(diào)整后,得到僅由排放源過程控制的污染物的長期序列.在進(jìn)行氣象調(diào)整的過程中,短期分量和基線分量基于氣象因素的逐步回歸方程及擬合度(R)如表2 所示,所有回歸方程的相關(guān)系數(shù)R 均超過0.3,說明污染物和各氣象要素間存在不同程度的相關(guān)性.

表2 2015～2022 年唐山市顆粒物和O3 短期、基線分量的多氣象因子逐步回歸方程Table 2 stepwise multi-meteorological linear regression equation of short-term and baseline components of particulate matter and O3 in Tangshan City from 2015 to 2022

PM2.5和PM10調(diào)整后的序列濃度在2018 年1月之前均低于長期序列濃度,說明從2015 年1 月到2018 年1 月,氣象因素有利于污染物的積累.此外,兩者在2022年調(diào)整后序列濃度均略高于長期序列,說明氣象因素有利于顆粒物的擴(kuò)散.O3調(diào)整后的序列濃度在2015 年、2019 年和2022 年大多低于長期序列,說明此期間氣象因素的變化主要有利于O3的生成,特別是短波輻射強(qiáng)度的變化,如圖 2 所示,2015、2019、2022 3a 夏季的短波輻射平均強(qiáng)度高于其余年份,更加有利于O3的生成.其余時(shí)期氣象因素則有利于O3的擴(kuò)散,而不利于O3的形成,可能主要與氣壓和相對(duì)濕度的變化有關(guān).張淼等[28]指出低壓系統(tǒng)下污染物的垂直運(yùn)動(dòng)增加,對(duì)O3的擴(kuò)散有一定的積極作用;而當(dāng)相對(duì)濕度升高時(shí),成云量也會(huì)增加,太陽輻射強(qiáng)度減弱,影響光化學(xué)反應(yīng)的速率,抑制O3的生成[29].

由于污染物氣象調(diào)整后的序列濃度僅代表源排放的貢獻(xiàn),污染物長期序列濃度與相應(yīng)氣象調(diào)整后序列濃度之差值就代表氣象因素對(duì)其長期分量的貢獻(xiàn).因此,可估算出 2015～2022 年間污染物PM2.5、PM10和O3源排放引起的濃度變化,氣象因素引起的濃度變化,以及這些變化占總變化的百分?jǐn)?shù),如表3 所示.對(duì)于PM2.5和PM10,源排放作用對(duì)于顆粒物濃度的貢獻(xiàn)約為 36.22～80.86,69.29～137.81μg/m3,源排放強(qiáng)度的減弱對(duì)于污染物濃度減少的貢獻(xiàn)超過90%,進(jìn)一步說明長期分量的變化受減排政策的影響較大,而氣象因素對(duì)于污染物濃度的貢獻(xiàn)以及污染物濃度變化的貢獻(xiàn)較小,因此對(duì)于污染源的控制是治理唐山市顆粒物污染的首要舉措.

表3 2015～2022 年唐山市顆粒物和O3 源排放及氣象因素的貢獻(xiàn)Table 3 Contribution of source emissions and meteorological factors of particulate matter and O3 in Tangshan City from 2015 to 2022

源排放作用對(duì)于O3濃度的貢獻(xiàn)約為96.57～102.86μg/m3,顯著高于氣象因素的貢獻(xiàn).氣象因素的變化對(duì)O3濃度改變的貢獻(xiàn)約為58.45%,超過了源的貢獻(xiàn),且該結(jié)果與Gao 等[30]對(duì)于河北省的研究結(jié)果相近.O3的生成過程主要是大氣中的NOx和VOCs在太陽光照射下發(fā)生光化學(xué)反應(yīng)的結(jié)果.溫度升高,輻射強(qiáng)度增大,O3的形成速率提升[31],因此每年夏季5～7 月是O3污染出現(xiàn)最多的時(shí)段,在治理唐山市O3污染時(shí),既要監(jiān)測O3前體物NOx和VOCs的變化情況,同時(shí)也要關(guān)注夏季光照、溫度等因素的變化對(duì)于光化學(xué)反應(yīng)速率的影響.

2.3 影響顆粒物與O3的重要因素

2.3.1 氣象因素對(duì)污染物的影響 T、RH、V 和Pre是影響PM2.5濃度的重要因素;而T、SR、RH 和Pre是影響O3生成的重要因素[32-34],如表2 所示.由于氣象因素在不同時(shí)間尺度上對(duì)PM2.5和O3的影響可能會(huì)有所不同,因此,為了進(jìn)一步探究氣象因素在不同時(shí)間尺度上對(duì)污染物PM2.5和O3的影響,本研究采用小波相干性分析法估算了2015～2022 年間不同時(shí)間尺度污染物PM2.5和O3分別與氣象要素的相干系數(shù)R2,探討了不同時(shí)間尺度氣象因素對(duì)PM2.5和O3的影響.

如圖4 所示,T 和RH 在小于32d 的時(shí)間尺度上與PM2.5濃度的R2大多在0.7 以上,如圖4 中圖(a)和(b)所示,箭頭向右,說明兩者均在該時(shí)間尺度上與PM2.5濃度呈正相關(guān)關(guān)系.且T 與PM2.5濃度的相位差不恒定,兩者間存在約1/8 周期的滯后關(guān)系,RH 與PM2.5濃度的相位差較小,特別是在4～16d 的時(shí)間尺度上兩者幾乎不存在相位差,這可能與小時(shí)間尺度上T 的影響不僅代表溫度對(duì)于PM2.5的影響,還反映了排放源的排放規(guī)律,而相對(duì)濕度的增加有利于空氣中細(xì)顆粒物的吸附作用,促進(jìn)二次組分的生成有關(guān)[35].256～512d 的時(shí)間尺度上兩者與PM2.5濃度的R2均大于0.9,箭頭向左,且有不同的相位差,說明T和RH 在該時(shí)間尺度上與PM2.5濃度呈負(fù)相關(guān)關(guān)系.地表垂直風(fēng)速V 與PM2.5濃度在4～190d 時(shí)間尺度上為R2在0.6～0.8,箭頭向右,且有不同的相位差,表明V在該時(shí)間尺度上與PM2.5濃度呈正相關(guān)關(guān)系,可能是因?yàn)轱L(fēng)速的增加有利于地面揚(yáng)塵進(jìn)入大氣中,造成大氣中顆粒物濃度升高.在256d～512d 的時(shí)間尺度上,V 與PM2.5濃度表現(xiàn)出負(fù)的相關(guān)性,如圖4(c)所示,說明大時(shí)間尺度上V 的增加有利于PM2.5濃度擴(kuò)散.單日累積降水量Pre 在256～512d 的時(shí)間尺度上與PM2.5濃度存在負(fù)相關(guān)關(guān)系,這可能與降水對(duì)PM2.5具有清除作用有關(guān)[36-37].

圖4 2015～2022 年不同時(shí)間尺度下PM2.5 與氣象因素的相關(guān)系數(shù)變化Fig.4 Variations of the coherence coefficients between PM2.5 and meteorological factors at different time scales from 2015 to 2022

如圖4 所示,T 和SR 與O3的R2在各時(shí)間尺度上幾乎均在0.65 以上,且箭頭向右,如圖4(e)和(g)所示,呈顯著的正相關(guān)關(guān)系,這可能是因?yàn)門 和SR的增加可以促進(jìn)O3的生成.但T 在64～192d 的尺度上、SR 在128～256d 的尺度上與O3濃度基本沒有明顯的相關(guān)性,這可能與兩種氣象要素的季節(jié)性變化有關(guān).RH 在小于32d 的時(shí)間尺度上與O3濃度的R2在0.6 以上,箭頭向左,呈現(xiàn)負(fù)相關(guān)關(guān)系,然而在256～512d 的時(shí)間尺度上與O3濃度的R2在0.9 以上,箭頭向右,呈現(xiàn)顯著的正相關(guān)關(guān)系,這可能與唐山市小時(shí)間尺度上RH 主要在55%～65%的范圍內(nèi),有利于O3濃度的增加,大時(shí)間尺度上RH 會(huì)超過80%,不利于O3的生成有關(guān)[44]. Pre 與O3濃度在0～128d的尺度上R2在0.4～0.8 之間,箭頭向左,呈現(xiàn)明顯負(fù)相關(guān)關(guān)系,在256～512d 尺度上在0.9 以上,箭頭向右,呈現(xiàn)顯著正相關(guān)關(guān)系,這可能是因?yàn)樘粕绞蠴3的高值時(shí)段集中于夏季,該時(shí)段內(nèi)降水頻繁且降水量較大,因此對(duì)O3的清除作用較為明顯,然而,大尺度上降水量因降水速率不同而對(duì)O3的作用有所不同,致使降水量與O3濃度呈正相關(guān)關(guān)系.

2.3.2 顆粒物與O3的相互影響 如圖5 所示,PM2.5和PM10之間在不同時(shí)間尺度上均呈現(xiàn)較高的相關(guān)性數(shù),且箭頭向右,在0～128d 的尺度上幾乎不存在相位差.但在128～256d 的時(shí)間尺度上,部分時(shí)段 PM2.5和PM10出現(xiàn)約1/8 周期的相位差,即兩者間存在16～32d 的滯后關(guān)系,而在256～512d 的時(shí)間尺度上相位差消失,說明PM2.5和PM10之間始終存在較強(qiáng)的正相關(guān)關(guān)系.PM2.5和O3在4～32d 的尺度上R2在0.4～0.9之間,大部分箭頭向右,主要呈現(xiàn)正相關(guān)關(guān)系.但在256～512d 的時(shí)間尺度上,PM2.5和O3相干系數(shù)大于0.9,箭頭向左,出現(xiàn)約1/2 周期的相位差,兩者間存在128～256d 的滯后關(guān)系,呈現(xiàn)負(fù)相關(guān)關(guān)系.小時(shí)間尺度上兩者的正相關(guān)性,一方面源于O3濃度的增加提升了大氣氧化能力,有助于二次污染物的生成;另一方面當(dāng)源排放作用增強(qiáng)時(shí),PM2.5濃度增加, VOCx和NOx等物質(zhì)作為O3的重要前體物也會(huì)促進(jìn)O3的生成.大時(shí)間尺度上兩者呈負(fù)相關(guān)性,可能是因?yàn)轭w粒物濃度增加會(huì)削弱太陽輻射強(qiáng)度,抑制O3的生成,再者,顆粒物對(duì)自由基還存在猝滅作用[39],顆粒物的增加會(huì)降低大氣中·OH 的濃度,從而減弱O3的生成.PM10和O3與PM2.5和O3的相關(guān)性圖譜非常相似,均呈現(xiàn)小時(shí)間尺度的正相關(guān)及大時(shí)間尺度的負(fù)相關(guān)關(guān)系.

圖5 2015～2022 年不同時(shí)間尺度下O3 與氣象因素的相關(guān)系數(shù)變化Fig.5 Variations of the coherence coefficients between O3 and meteorological factors at different time scales from 2015 to 2022

3 結(jié)論

3.1 2015～2022 年間,唐山市大氣污染物PM2.5、PM10和O3濃度的長期分量、季節(jié)分量和短期分量能很好的解釋其原始序列濃度變化.顆粒物PM2.5和PM10濃度長期分量均呈現(xiàn)顯著下降趨勢,且主要受污染源排放的影響,而短期分量和季節(jié)分量均在12～次年2 月之間呈現(xiàn)不同程度的周期波動(dòng),主要受源排放和氣象因素的影響.O3濃度長期分量變化幅度較小,且除受到污染源排放的影響外,還受到氣象因素的影響,其短期分量和季節(jié)分量均在5～7 月之間出現(xiàn)不同程度的周期波動(dòng),主要受氣象因素的影響.

3.2 KZ 濾波法與多元逐步回歸法相結(jié)合是一種有效區(qū)分和定量估算源排放和氣象因素對(duì)大氣污染物濃度長期分量變化貢獻(xiàn)的技術(shù).2015～2022 年間,唐山市大氣顆粒物PM2.5和PM10濃度的長期分量變化中,源排放貢獻(xiàn)均占90%以上,而氣象因素的貢獻(xiàn)則占10%以下.O3濃度的長期分量變化,源排放的貢獻(xiàn)占41.55%,氣象因素的貢獻(xiàn)占58.45%.因此,對(duì)于O3污染的治理不僅要控制源排放強(qiáng)度,也需要關(guān)注氣象因素對(duì)其前體物VOCx和NOx向O3轉(zhuǎn)化過程的影響.

3.3 2015～2022 年間,唐山市氣象因素溫度、相對(duì)濕度、地面垂直風(fēng)速和降水量對(duì)PM2.5的影響主要表現(xiàn)為小時(shí)間尺度的正向作用和大時(shí)間尺度的負(fù)向作用.而氣象因素溫度和短波輻射強(qiáng)度對(duì)O3的影響在不同時(shí)間尺度大多呈正向作用,相對(duì)濕度和降水量對(duì)O3的影響則表現(xiàn)為小時(shí)間尺度的負(fù)向影響和大時(shí)間尺度的正向影響.顆粒物PM2.5、PM10和O3之間存在小時(shí)間尺度的正向作用以及大時(shí)間尺度的負(fù)向影響.