李 嬡 魯紹偉 李少寧 趙 娜 徐曉天

（1. 沈陽農(nóng)業(yè)大學(xué)，遼寧 沈陽 110866；2. 北京市林業(yè)果樹科學(xué)研究院，北京燕山森林生態(tài)系統(tǒng)定位觀測研究站，北京 100093）

空氣負(fù)離子（NAI）是帶負(fù)電荷的單個氣體分子和輕離子團(tuán)的總稱[1]。Luts[2]指出負(fù)離子在穩(wěn)定的狀態(tài)下，以O(shè)2·-占據(jù)主導(dǎo)地位，因此NAI又稱負(fù)氧離子。NAI作為衡量和評價空氣質(zhì)量高低的重要指標(biāo)之一，對人們的情緒、記憶以及行為有影響，有助于使人在運動后恢復(fù)生理反應(yīng)[3]，還能減少霧氣、煙霧、車輛尾氣對人體的傷害等[4]，具有殺菌、降塵、清潔空氣的作用[5]。我國已于2008年將NAI濃度列入森林生態(tài)系統(tǒng)服務(wù)功能的一項評價指標(biāo)。這些研究表明，NAI對人類的健康問題以及生存環(huán)境有著重要的作用。

臭氧（O3）是地球大氣中一種微量氣體，它是由于大氣中氧分子受太陽輻射分解成氧原子后，氧原子又與周圍的氧分子結(jié)合而形成的，含有3個氧原子。O3污染對大氣能見度、人類健康和氣候變化造成了負(fù)面影響，尤其是在人類健康方面，許多的流行病學(xué)研究發(fā)現(xiàn)O3與早死有關(guān)聯(lián)[6-7]，最近一項研究報告指出，在中國經(jīng)濟(jì)發(fā)達(dá)地區(qū)（如京津冀地區(qū)），O3已取代PM2.5成為主要污染物[8]。近年來對于NAI濃度的研究大多數(shù)都密切關(guān)注其自身的時空變化特征[9-11]，也有研究發(fā)現(xiàn)NAI濃度與溫度、濕度、環(huán)境污染、風(fēng)效應(yīng)指數(shù)密切相關(guān)[12-15]，與大氣污染物（如顆粒物、二氧化氮、二氧化硫等）有一定的相關(guān)性[16-17]，而NAI與O3都來自于空氣中的氧分子，卻很少有人關(guān)注于二者的相關(guān)關(guān)系。

城市森林是城市中有生命的基礎(chǔ)設(shè)施，與城市居民的聯(lián)系最為密切，具有一系列的生態(tài)保健功能，在保護(hù)人類生存和健康方面發(fā)揮著越來越重要的作用，其生態(tài)功能和生態(tài)效益日益受到環(huán)保、園林部門的重視以及城市居民的關(guān)注[18]。城市森林環(huán)境中的空氣負(fù)氧離子能促進(jìn)和維護(hù)城市居民健康[19]，而且森林環(huán)境中樹木葉片的存在更有利于植物清除包含O3等在內(nèi)空氣污染物[20]。因此，本研究致力于探討城市森林公園中NAI濃度和O3濃度的時間和空間變化特征以及二者的相互作用，為城市公園綠地規(guī)劃和管理提供參考依據(jù)，對建設(shè)高質(zhì)量的城市公園綠地具有實際應(yīng)用價值，對指導(dǎo)居民進(jìn)行合理的綠地游憩活動具有借鑒意義。

1 研究方法

1.1 研究地概況

研究地位于北京市（東經(jīng)115.7°～117.4°，北緯39.4°～41.6°），屬于北溫帶半濕潤大陸性季風(fēng)氣候，夏季高溫多雨，冬季寒冷干燥，春、秋短促。全年無霜期180～200 d，西部山區(qū)較短。降水季節(jié)分配很不均勻，全年降水的80%集中在夏季6—8月，7—8月有大雨。

研究選取北京市的4個城市森林公園：北京朝陽公園（東經(jīng)116.45°，北緯39.92°）、南海子郊野公園（東經(jīng)116.47°，北緯39.77°）、北京西山森林公園（東經(jīng)116.19°，北緯39.97°）和北京松山國家級自然保護(hù)區(qū)（東經(jīng)115.82°，北緯40.50°），園內(nèi)均設(shè)有城市森林生態(tài)環(huán)境監(jiān)測站，4個城市森林公園的環(huán)境特征見表1。監(jiān)測點主要綠地類型為針闊、喬灌草混交林，森林樹種主要有側(cè)柏（Platycladus orientalis）、油松（Pinus tabulaeformis）、刺 槐（Robinia pseudoacacia）、蒙古櫟（Quercus mongolica）等基本代表了北京城市森林公園的植被特征。

表 1 城市森林公園的具體環(huán)境特征Table 1 Environmental characteristics of urban forest parks

1.2 數(shù)據(jù)采集

4個城市森林生態(tài)環(huán)境監(jiān)測站分別設(shè)置49i型雙光室紫外光度法O3分析儀（Thermo Fisher Scientific，美國）和Yconetest-500負(fù)離子自動監(jiān)測系統(tǒng)（廣東佛山云創(chuàng)氣象服務(wù)技術(shù)有限公司，中國），對O3濃度及NAI濃度實行野外自動在線監(jiān)測，頻次為1 h測定1次，全天24 h不間斷采樣。根據(jù)氣候統(tǒng)計法，選取2018年夏季（2018年6—8月）和2018年冬季（2018年12月至2019年2月）的數(shù)據(jù)進(jìn)行分析。

1.3 數(shù)據(jù)處理與分析

通過計算NAI濃度和O3濃度的日平均濃度來研究其季節(jié)變化；同時為降低誤差，選取了2個季節(jié)晴朗無風(fēng)連續(xù)3 d的數(shù)據(jù)，平均了所有時間點的濃度來研究NAI濃度和O3濃度的日變化。世界衛(wèi)生組織規(guī)定的清新空氣負(fù)氧離子標(biāo)準(zhǔn)濃度為1 000～1 500 ion/cm3，一級空氣質(zhì)量標(biāo)準(zhǔn)濃度為2 000 ion/cm3以上[21]。每日最大1 h平均濃度是指任何1 h污染物濃度的算術(shù)平均值，每日最大8 h平均濃度指連續(xù)8 h污染物的算術(shù)平均值，也稱8 h滑動平均。 GB 3095—2012[22]規(guī)定O3每日最大1 h濃度（O31 h max）和每日最大8 h平均濃度（O38 h max）一級標(biāo)準(zhǔn)限值分別為160 μg/m3和100 μg/m3。通過研究污染比較嚴(yán)重的公園來確定北京市森林公園的NAI濃度是否達(dá)到標(biāo)準(zhǔn)，由于南海子公園是北京市的近郊開發(fā)區(qū)，污染相對比較嚴(yán)重，且植被覆蓋率在4個公園中最低，生成負(fù)離子濃度條件不高，所以本文只研究了南海子公園的季節(jié)變化。

根據(jù)各站點的地理位置劃分北京市各區(qū)域，統(tǒng)計不同區(qū)域NAI濃度和O3濃度，比較空間變化。為了分析NAI濃度與O3的相關(guān)性，采用雙尾檢驗進(jìn)行了統(tǒng)計分析，結(jié)果表明這些數(shù)據(jù)不服從正態(tài)分布。然后，計算Pearson相關(guān)系數(shù)量化NAI濃度和O3之間的相關(guān)性。所有P值均為雙尾，P＜0.05具有統(tǒng)計學(xué)意義。所有的分析均使用Origin 2018軟件。

2 結(jié)果與分析

2.1 NAI濃度與O3濃度的時間變化特征

2.1.1 NAI濃度與O3濃度的季節(jié)變化特征

北京近郊開發(fā)區(qū)南海子公園夏季NAI濃度和O3日均濃度明顯高于冬季，見圖1。由圖1可知，夏冬季NAI日均濃度分別為（1 761±57）ion/cm3和（1 481±78）ion/cm3，夏季NAI日均濃度比冬季高出300 ion/cm3左右；O3日均濃度分別為（36±10）μg/m3和（21±10）μg/m3，夏季O3平均濃度約為冬季的1.7倍。夏季剛開始NAI濃度變化趨勢比較穩(wěn)定，7月第2周開始NAI濃度變化幅度較大，最大值與最小值相差353 ion/cm3；冬季NAI濃度在1 500 ion/cm3上下浮動，在2月的時候突然下降到1 500 ion/cm3以下。夏季O3濃度在7月中旬之前O3日均濃度變化浮動較大，最大值與最小值相差45 μg/m3，中旬之后O3濃度降低，變化趨于穩(wěn)定；冬季O3濃度的變化趨勢呈逐漸遞增的趨勢，12月出現(xiàn)濃度的最低值6 μg/m3，次年2月濃度達(dá)到了最大值為46 μg/m3。因此，夏季NAI濃度和O3濃度的變化趨勢相反，冬季二者變化趨勢不明顯。

圖 1 南海子公園NAI和O3日均濃度Fig. 1 Diurnal concentration of NAI and O3 in Nanhaizi Park

2.1.2 NAI濃度與O3濃度的日變化特征

北京城市森林夏季和冬季NAI濃度呈現(xiàn)出不同的變化趨勢，見圖2。由圖2可知，夏季為“單峰單谷”曲線，冬季NAI濃度呈現(xiàn)“單峰”曲線；2個季節(jié)的O3濃度日變化趨勢相似，均呈現(xiàn)“單峰單谷”曲線。夏季NAI濃度日變化的總體規(guī)律為清晨（5：00—9：00）NAI濃度為高值區(qū)，出現(xiàn)當(dāng)日的最高峰值；中午（11：00—15：00）為低值區(qū)，出現(xiàn)當(dāng)日的最低峰值，傍晚之后NAI濃度開始回升。如圖2a所示，3個區(qū)域的NAI濃度日變化均表現(xiàn)出相似的特征，分別在7：00—9：00出現(xiàn)峰值，在15：00左右出現(xiàn)谷值。市中心、近郊淺山林區(qū)和遠(yuǎn)郊山地林區(qū)NAI濃度最高分別為2 348、1 883 ion/cm3和2 465 ion/cm3，最低分 別 為2 016、1 590 ion/cm3和2 037 ion/cm3。冬季NAI濃度則呈現(xiàn)出與夏季不同的模式，3個區(qū)域NAI濃度變化趨勢相似（如圖2b所示），呈現(xiàn)“單峰”曲線，NAI濃度在0：00—7：00為低值區(qū)，7：00—8：00開始逐漸上升，在10：00—11：00左右出現(xiàn)一天當(dāng)中的最高峰值，在2 000～2 300 ion/cm3之間，12：00—13：00開始下降，19：00—23：00回到低值區(qū)，隨后趨于穩(wěn)定。

夏季和冬季不同區(qū)域的O3濃度差異比較明顯，且遠(yuǎn)郊山地林區(qū)O3濃度始終低于其他2個區(qū)域。夏季O3濃度最低出現(xiàn)在早上7：00—8：00，遠(yuǎn)郊區(qū)濃度最低，不超過20 μg/m3，而其他2個區(qū)域最低濃度為20～50 μg/m3；3個區(qū)域的O3濃度從8：00開始增加，直到15：00—17：00左右到達(dá)峰值，市中心O3濃度最高，達(dá)到200 μg/m3以上，近郊區(qū)次之，125 μg/m3，而遠(yuǎn)郊區(qū)未到50 μg/m3；隨后呈現(xiàn)下降的趨勢，直到第2天清晨。冬季O3日變化規(guī)律與夏季相似，但波動幅度較小。O3濃度最低出現(xiàn)在8：00（遠(yuǎn)郊區(qū)除外），15：00時上升到最高值，然后逐漸下降，直到次日早上。

圖 2 夏季和冬季NAI和O3平均日變化Fig. 2 Diurnal variation of NAI and O3 in summer and winter

2.2 NAI濃度與O3濃度的空間變化特征

由圖3a可知，夏季市中心公園O3濃度最大為（84±64）μg/m3，近郊淺山林區(qū)濃度次之（58±45）μg/m3，然后是近郊開發(fā)區(qū)（36±10）μg/m3，遠(yuǎn) 郊 區(qū)O3濃 度 最 ?。?2±13）μg/m3，北 京 市中心污染程度比較嚴(yán)重，是近郊區(qū)O3濃度的1.5～2.5倍，是遠(yuǎn)郊區(qū)的4倍左右。北京冬季O3濃度空間變化與夏季相似，從市中心到郊區(qū)呈遞減的趨勢，但整體變化趨勢不明顯，濃度變化范圍在15～30 μg/m3之間。

圖 3 夏冬季北京各區(qū)域NAI和O3的平均濃度Fig. 3 Average concentration of NAI and O3 in various regions of Beijing in summer and winter

由圖3b可知，NAI濃度空間變化特征夏季表現(xiàn)為市中心公園＞遠(yuǎn)郊區(qū)公園＞近郊淺山林區(qū)＞近郊開發(fā)區(qū)，市中心公園NAI濃度最大，比遠(yuǎn)郊區(qū)、近郊區(qū)淺山林區(qū)、近郊開發(fā)區(qū)濃度分別上升了7%、14%、19%；冬季NAI濃度變化趨勢與之類似，NAI濃度為1 500～2 000 ion/cm3，市中心NAI濃 度 最 大（2 002±90）ion/cm3，比 遠(yuǎn) 郊 區(qū)、近郊淺山林區(qū)、近郊開發(fā)區(qū)分別升高7%、9%、26%。夏冬季遠(yuǎn)郊區(qū)公園的NAI濃度次于市中心公園，說明NAI濃度變化與北京地區(qū)污染程度總趨勢不一致。

2.3 NAI濃度與O3濃度的相互作用

由圖4可知，夏季NAI濃度與O3濃度呈負(fù)相關(guān)關(guān)系（P＜0.01，r=-0.303），冬季NAI濃度與O3濃度呈正相關(guān)關(guān)系（P＜0.01，r= 0.388）。

圖 4 南海子公園夏季和冬季NAI濃度和O3濃度的散點圖Fig. 4 Scattered points of NAI and O3 concentrations in Nanhaizi Park in summer and winter

由圖5a可知，研究了夏季O3濃度對NAI濃度的影響，在4個O3max（日最大O3濃度）水平下，NAI濃度的日變化規(guī)律，發(fā)現(xiàn)越高的O3max濃度水平下，NAI濃度越低。當(dāng)O3max＞50 μg/m3時，NAI濃度最低，且日變化顯著，最低值出現(xiàn)在14：00；當(dāng)20 μg/m3＜O3max＜50 μg/m3時，NAI濃度較低，最低值出現(xiàn)在18：00；當(dāng)O3max＜20 μg/m3時，NAI濃 度 最 高，日 變 化 不 顯 著?？?知，當(dāng)O3max＞50 μg/m3時，會 影 響NAI濃 度變化。

圖 5 南海子公園不同濃度水平下NAI濃度與O3濃度的日變化Fig. 5 Diurnal changes of NAI and O3 concentrations in Nanhaizi Park

由圖5b可知，在冬季3個NAI濃度水平下，O3濃度呈現(xiàn)出相似的變化曲線。當(dāng)NAI濃度＞1 500 ion/cm3時，O3濃度最高，每小時變化較大；當(dāng)1 200 ion/cm3＜NAI濃度＜1 500 ion/cm3時，O3濃度與上一個NAI濃度水平相比，進(jìn)一步降低；當(dāng)NAI濃度＜1 200 ion/cm3時，O3濃度最低，每小時變化較小?？傮w而言，冬季NAI濃度水平越高，O3濃度越高。

3 結(jié)論與討論

北京近郊南海子公園夏季和冬季NAI日均濃度分別為（1 761±57）ion/cm3和（1 481±78）ion/cm3，達(dá)到了世界衛(wèi)生組織規(guī)定的清新空氣負(fù)氧離子標(biāo)準(zhǔn)濃度。夏季和冬季O3日均濃度分別為（36±10）μg/m3和（21±10）μg/m3，且O31 h max和O38 h max均小于100 μg/m3，達(dá)到了GB 3095—2012[22]規(guī)定標(biāo)準(zhǔn)限值，說明北京城市森林公園符合我國風(fēng)景區(qū)等空氣清潔地區(qū)的水平。

研究發(fā)現(xiàn)北京森林公園夏季NAI濃度和O3濃度均比冬季高。夏季植物受太陽日照時間長、紫外線強等因素，光合作用旺盛，其生長茂盛，降水量豐富，空氣濕度較大，為NAI的生產(chǎn)提供了極為有利的條件[23-24]，因此，夏季NAI濃度較高，這與國內(nèi)其他森林公園的觀測結(jié)果一致[25]。O3的形成與其前體物和氣象因素有關(guān)[26]，夏季高溫和強的太陽輻射會導(dǎo)致光化學(xué)反應(yīng)速率的提高，從而促進(jìn)O3的生成[27]，所以夏季O3濃度高于冬季。

已有研究表明，NAI濃度的日變化峰值出現(xiàn)在日出前后，隨后開始下降，最低值出現(xiàn)在下午，然后又逐步升高[24-25,28]，但這些研究偏向于夏季和秋季，對于冬季NAI濃度的日變化幾乎沒有。本研究中夏季NAI濃度日變化趨勢和前人的研究一致，呈“單峰單谷”型，冬季則出現(xiàn)不同的變化趨勢，呈“單峰”型，且在正午11：00—12：00出現(xiàn)峰值，這是因為冬季正午的太陽輻射是一天當(dāng)中最強的時候，空氣會電離產(chǎn)生更多的NAI。O3濃度在2個季節(jié)的日變化基本相同，夏季17：00和冬季15：00出現(xiàn)單峰。在這兩個季節(jié)，O3濃度峰值存在差別是由于夏季一直持續(xù)更高的太陽輻射和溫度而促進(jìn)光化學(xué)反應(yīng)形成O3[29]。

根據(jù)離城市中心距離可知，4個城市森林污染程度分別為朝陽公園＞南海子公園＞西山公園＞松山自然保護(hù)區(qū)，可以看出NAI濃度變化與這4個地區(qū)污染程度總趨勢不一致，說明與其他影響因素相比，城市污染梯度對城市森林中NAI濃度影響較小。邵海榮等研究表明北京地區(qū)NAI濃度從市中心向近郊、遠(yuǎn)郊逐漸增大[11]。Liang等[30]研究發(fā)現(xiàn)城市化可以引起NAI濃度沿城鄉(xiāng)梯度的變化，但不會改變其時間變化規(guī)律；城市化對可以誘發(fā)NAI濃度沿城鄉(xiāng)梯度變化，但僅限于非植被或植被較少的地區(qū)（如開放空間），而在植被豐富的地區(qū)（如城市森林）則不明顯，是因為森林有更多的能夠促進(jìn)電解的葉片生物量，其次森林釋放的芳香物質(zhì)也可以促進(jìn)空氣電離，增加森林中的NAI濃度[31]，這也是本研究中NAI濃度沒有隨著污染程度嚴(yán)重而呈現(xiàn)下降的趨勢的原因。

O3濃度在2個季節(jié)中由市中心公園向郊區(qū)公園呈下降趨勢，這種模式在夏季更為明顯，這與四個城市森林公園的污染程度有極大的關(guān)系。張紅星等[32]采用移動監(jiān)測車對北京城市道路進(jìn)行多點位連續(xù)測定地表O3濃度，其采樣管距離地面約3.3 m，研究發(fā)現(xiàn)北京O3濃度從市中心城區(qū)向郊區(qū)公園呈現(xiàn)上升趨勢。王占山等[33]通過對北京市35個自動空氣監(jiān)測子站的O3濃度進(jìn)行分析，得出北京市城區(qū)O3濃度相對較低，周邊區(qū)縣相對較高，生態(tài)植被優(yōu)良的東北地區(qū)最高。本研究與前人的研究內(nèi)容同為北京市O3的空間變化，不同的是，本研究對象是不同城市森林公園O3濃度的變化。由此可以看出，在大尺度上O3濃度從城區(qū)到遠(yuǎn)郊逐漸增加，而城市公園的存在改變了這一規(guī)律，使得城區(qū)公園內(nèi)的O3濃度比郊區(qū)公園高，原因是城區(qū)公園四周為高樓大廈，風(fēng)速較低，不利于O3的擴(kuò)散，且城區(qū)內(nèi)溫度較高，光化學(xué)反應(yīng)劇烈，有利于O3產(chǎn)生。

研究從不同季節(jié)的角度分析NAI濃度和O3濃度的相關(guān)性，發(fā)現(xiàn)NAI濃度和O3濃度在夏季呈顯著負(fù)相關(guān)關(guān)系，相關(guān)系數(shù)為-0.303，在冬季呈顯著正相關(guān)關(guān)系相關(guān)系數(shù)為0.388。趙艷佩[14]研究發(fā)現(xiàn)NAI濃度和O3濃度在1—5月呈顯著正相關(guān)，相關(guān)系數(shù)為0.265，同時也驗證了O3與其他污染物均呈顯著負(fù)相關(guān)，和NAI與其他污染物的相關(guān)關(guān)系一致，從側(cè)面驗證了O3與NAI之間的正相關(guān)性；曹建新等[9]研究發(fā)現(xiàn)NAI濃度和O3濃度全年呈顯著正相關(guān)，相關(guān)系數(shù)為0.217，NAI濃度和O3濃度呈正相關(guān)，這是因為空氣電離產(chǎn)生NAI過程中釋放的氧原子，會與空氣中的氧氣結(jié)合生成O3。

夏季NAI濃度和O3濃度呈負(fù)相關(guān)關(guān)系，可能與夏季植物生長茂盛有關(guān)。前人的研究不是在城市森林或公園中或者沒有單獨討論夏季二者的相關(guān)性，因此沒有植物的作用。有研究表明日間植物通過氣孔途徑吸收生態(tài)系統(tǒng)中80%的O3[34]；也有研究發(fā)現(xiàn)有些植物可以吸收空氣中的O3，并利用其體內(nèi)的一系列酶如超氧化物歧化酶等進(jìn)行轉(zhuǎn)化清除[35]；且植物自身釋放的一氧化氮和萜類等氣體能中和大氣中的O3等一系列研究表明植物降低空氣中的O3濃度[36]。相反，某些植物的葉片呈針狀等曲率半徑較小，具有尖端放電的功能，使空氣發(fā)生電離，NAI濃度增加；植物的葉表面在短波紫外線的作用下，也能夠發(fā)生光電效應(yīng)，使NAI濃度升高[37]，李少寧等研究發(fā)現(xiàn)樹木生長季能提供更多的NAI[38]，所以說植物會影響NAI濃度和O3濃度相關(guān)性。

研究選取夏季和冬季北京市中心到遠(yuǎn)郊區(qū)四個不同區(qū)域NAI濃度和O3濃度的數(shù)據(jù)進(jìn)行分析，探討NAI濃度和O3濃度的時空變化特征及其相互作用，得到如下結(jié)論。時間變化上，北京城市森林夏季NAI濃度和O3平均濃度均高于冬季；夏季和冬季NAI濃度日變化分別呈現(xiàn)“單峰單谷”和“單峰”曲線；兩季節(jié)O3日變化趨勢相似，均呈“單峰單谷”曲線。夏季NAI濃度分別在7：00—9：00出現(xiàn)峰值和15：00左右出現(xiàn)谷值；冬季NAI濃度峰值在12：00左右；夏季和冬季O3濃度谷值在7：00—9：00，冬季峰值在15：00，比夏季早兩個小時左右。二者在夏季和冬季均具有明顯的日變化特征?？臻g變化上，NAI濃度變化表現(xiàn)為：市中心公園＞遠(yuǎn)郊區(qū)公園＞近郊淺山林區(qū)＞近郊開發(fā)區(qū)；O3濃度變化表現(xiàn)為：市中心公園＞近郊淺山林區(qū)＞近郊開發(fā)區(qū)＞遠(yuǎn)郊區(qū)?？梢姵鞘形廴咎荻葘Τ鞘猩止珗@NAI濃度影響小，而O3濃度空間變化符合城市污染梯度。夏季NAI濃度與O3濃度呈負(fù)相關(guān)關(guān)系（P＜0.01，r=-0.303），冬季NAI濃度與O3濃度呈正相關(guān)關(guān)系（P＜0.01，r=0.388），說明植物的生長季與NAI濃度和O3濃度的相關(guān)性有關(guān)。