李 品,馮兆忠,*,尚 博,袁相洋,代碌碌,徐彥森

1 中國(guó)科學(xué)院生態(tài)環(huán)境研究中心 城市與區(qū)域生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100085 2 中國(guó)科學(xué)院大學(xué)資源環(huán)境學(xué)院,北京 101408

地表臭氧(O3)主要由氮氧化物(NOx)和揮發(fā)性有機(jī)化合物(VOCs)等前體物經(jīng)光化學(xué)反應(yīng)生成,是一種具有極強(qiáng)氧化毒性作用的二次污染物和重要的溫室氣體[1]。隨著我國(guó)工業(yè)化與城市化進(jìn)程的加快,NOx和VOCs等前體物顯著增加[2],地表O3濃度也不斷升高。根據(jù)環(huán)保部網(wǎng)站監(jiān)測(cè)數(shù)據(jù)顯示,O3正在成為我國(guó)主要城市夏季的首要空氣污染物,其中夏季小時(shí)平均O3濃度高達(dá)300nmol/mol以上,年均最大O3濃度達(dá)到60nmol/mol[3- 4],遠(yuǎn)超過敏感樹種的受害臨界負(fù)荷[5],對(duì)樹木生長(zhǎng)和森林生態(tài)系統(tǒng)服務(wù)功能產(chǎn)生嚴(yán)重威脅[6]。

城市樹木在維持城市生態(tài)平衡,改善和緩解城市生態(tài)環(huán)境問題等方面起著無可替代的重要作用。但隨著O3濃度的升高,研究發(fā)現(xiàn)暴露在O3下的樹木出現(xiàn)一系列的傷害表征,如葉片出現(xiàn)缺綠斑、變黃、早衰、脫落等可見傷害癥狀[7- 9],光合速率降低[10- 12],氣孔反應(yīng)滯后甚至失靈[13- 14],水分消耗量增加[14- 15],樹木生物量降低[16],改變生物量分配格局[11]和樹木對(duì)病蟲害以及環(huán)境脅迫的抵御能力[17],從而降低城市樹木作為“碳匯”和“清潔工”改善生態(tài)環(huán)境問題的潛力[18- 19]。

植物氣孔是進(jìn)行光合作用和蒸騰作用等氣體交換的重要調(diào)控器官,其氣孔特征(數(shù)量、密度、開度、大小等)的變化體現(xiàn)了對(duì)環(huán)境變化及逆境脅迫響應(yīng)與適應(yīng)的結(jié)果。O3作為逆境脅迫因子,通過氣孔進(jìn)入植物體內(nèi),氣孔特征決定了O3的進(jìn)入量[20],也是評(píng)價(jià)O3對(duì)植物傷害的關(guān)鍵因素[21]。目前,高濃度O3對(duì)城市綠化樹種的葉片氣孔數(shù)量特征影響的研究有一些報(bào)道[22- 25],得出有差異的結(jié)論：高濃度O3處理顯著降低葉片氣孔開度和大小,但對(duì)氣孔密度沒有顯著影響[22, 24]或者增加了氣孔密度[23, 25]。以上研究都是基于單個(gè)樹種單一O3升高濃度的熏蒸結(jié)果,由于基于不同的樹種、O3濃度和實(shí)驗(yàn)持續(xù)期,很難評(píng)價(jià)O3對(duì)植物氣孔特征的影響以及對(duì)植物產(chǎn)生的定量傷害。目前,根據(jù)不同濃度梯度的O3與植物生理指標(biāo)、生物量、產(chǎn)量、可見傷害癥狀等建立劑量響應(yīng)關(guān)系[16, 26],從而定量評(píng)價(jià)O3對(duì)植物的傷害,是有效可行的方法?？紤]到樹木生產(chǎn)緩慢,若干年后收獲生物量來評(píng)估在實(shí)際過程中很難實(shí)現(xiàn)。而根據(jù)生理指標(biāo),特別是敏感性的指標(biāo)——?dú)饪滋卣?進(jìn)行相應(yīng)的O3劑量關(guān)系研究,尚未見報(bào)道。本實(shí)驗(yàn)選取的6種樹種,作為城市行道木、風(fēng)景園林綠化、防護(hù)林體系構(gòu)成在北方地區(qū)廣泛栽種。通過研究不同O3濃度脅迫對(duì)這6種城市綠化樹種葉片氣孔的影響效應(yīng),并根據(jù)O3濃度梯度下氣孔特征指標(biāo)的響應(yīng),建立O3濃度與氣孔指標(biāo)的劑量關(guān)系,以期為地表O3污染下篩選抗污染的樹種及保護(hù)城市綠化樹木少受O3傷害提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗(yàn)區(qū)概況

實(shí)驗(yàn)區(qū)位于北京市延慶區(qū)延慶鎮(zhèn)唐家堡村(40°45′N, 115°97′E),在北京市區(qū)的東北方向,屬大陸性季風(fēng)氣候,是溫帶與中溫帶、半干旱與半濕潤(rùn)帶的過渡連帶。氣候冬冷夏涼,年平均氣溫10.2℃,夏季平均氣溫22—25℃,年降水量567mm,年平均日照2800h。

1.2 實(shí)驗(yàn)設(shè)計(jì)

實(shí)驗(yàn)選用6種城市綠化樹種的一年生幼苗：2種不同基因型的楊樹‘546’(Populusdeltoidescv. ‘55/56’ ×P.deltoidescv. ‘Imperial’)和‘107’(P.euramericanacv. ‘74/76’)、白蠟(Fraxinuschinensis)、法桐(Platanusorientalis)、洋槐(Robiniapseudoacacia)、國(guó)槐(Sophorajaponica)。2016年4月10日將長(zhǎng)勢(shì)較為一致的幼苗統(tǒng)一移入20L的圓形花盆中,用當(dāng)?shù)赝寥琅c購(gòu)買的草炭土按體積比3∶1混勻后裝盆作為栽培基質(zhì)。待移植兩周后,選擇樹高和基徑相差不大的幼苗,移入開頂式氣室(OTCs)中,每間隔1—2d澆水,以適應(yīng)氣室內(nèi)環(huán)境。10d后進(jìn)行O3熏蒸處理。

本實(shí)驗(yàn)在15個(gè)OTC中進(jìn)行O3熏氣,OTC為高3.0m,邊長(zhǎng)1.0m的正八面柱體。實(shí)驗(yàn)共設(shè)5個(gè)處理：CF(過濾大氣,charcoal filtered ambient air,[O3]<40nmol/mol),NF(環(huán)境大氣[O3])、NF20(NF+20nmol/mol[O3])、NF40(NF+40nmol/mol[O3])、NF60(NF+60nmol/mol [O3])。每個(gè)處理有3個(gè)OTC重復(fù),每個(gè)OTC內(nèi)每種樹5—6盆,并隨機(jī)分布。實(shí)驗(yàn)期間各個(gè)罩子的O3濃度由美國(guó)熱電公司49i臭氧分析儀監(jiān)測(cè)(Model 49i; Thermo Scientific, Franklin, MA, USA),熏氣時(shí)間從2016年6月26日到9月30日,持續(xù)96d,每天熏氣時(shí)間為10h(08:00—18:00),陰雨天停止熏氣。每個(gè)OTC外面和內(nèi)部的溫度和濕度由CR1000數(shù)據(jù)記錄器監(jiān)測(cè)(Campbell Scientific, North Logan, Utah, USA),每5min記錄一次。在整個(gè)熏氣期間,OTC內(nèi)部的溫度比外面環(huán)境平均高1.7℃,濕度平均低4.5%。分別于8月中上旬(生長(zhǎng)旺期)和9月中下旬(生長(zhǎng)末期),分兩次進(jìn)行取樣測(cè)定指標(biāo)。所有樹種在取樣時(shí)的實(shí)驗(yàn)期間小時(shí)O3濃度超過40nmol/mol的累計(jì)值(AOT40,the accumulated hourly mean O3concentration over a threshold of 40nmol/mol)詳見表1。

表1 6種樹種2次取樣時(shí)的AOT40值

1.3 指標(biāo)測(cè)定

在8月中上旬(生長(zhǎng)旺期)和9月中下旬(生長(zhǎng)末期)兩個(gè)時(shí)間段,分別在每個(gè)OTC中選擇每種樹苗3盆,于9:00—11:00間采集中上部(從頂端數(shù)第5—8片)陽生葉片,用指甲油印跡法進(jìn)行取樣：輕拭葉片表面灰塵,選擇葉片遠(yuǎn)軸面中部靠近葉脈部位,均勻涂上無色指甲油,將其自然晾干,用鑷子撕下已干的指甲油薄膜層,制作裝片進(jìn)行顯微鏡觀察和測(cè)量[27- 28]。

用Dino-Lite AM7013MZT4(R4)顯微鏡照相儀拍照,用DinoCapture 2.0軟件計(jì)算每幅圖片上選定的1mm2視野面積中的氣孔個(gè)數(shù)(單位面積植物葉片包含的氣孔數(shù)目),即為氣孔密度(氣孔數(shù)/mm2)。隨機(jī)選取每幅圖片上完整清晰的10個(gè)氣孔,用測(cè)微尺測(cè)量氣孔長(zhǎng)度、寬度,氣孔長(zhǎng)度為平行于氣孔器的最長(zhǎng)值,氣孔寬度為垂直于氣孔器的最寬值。計(jì)算氣孔開度(氣孔開度=πab,其中a表示1/2氣孔長(zhǎng)度,b表示1/2氣孔寬度)和氣孔大小(氣孔大小=氣孔長(zhǎng)度×氣孔寬度),然后取10個(gè)氣孔的平均值作為這幅圖片上氣孔特征的值。

1.4 數(shù)據(jù)分析

實(shí)驗(yàn)數(shù)據(jù)按OTC為重復(fù)單元(n=3)。運(yùn)用SPSS 18.0軟件進(jìn)行統(tǒng)計(jì)分析,選擇Tukey′s Honestly Significant Difference (HSD) 檢驗(yàn)法分析被測(cè)指標(biāo)在不同O3濃度之間的差異,并選擇一般線性模型(General Linear Model)將O3處理、采樣時(shí)間與樹種進(jìn)行多因素方差分析。文中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)差(Mean±SD)。運(yùn)用Sigmaplot 10.0繪制氣孔特征指標(biāo)與AOT40的線性關(guān)系圖(普通最小二乘法),并用Spearman′s相關(guān)關(guān)系檢驗(yàn)線性關(guān)系的顯著性。圖中實(shí)線代表線性擬合關(guān)系顯著(P<0.05),虛線代表不顯著(P>0.05)。

2 結(jié)果與分析

2.1 O3濃度升高對(duì)氣孔密度的影響

三因素方差分析結(jié)果表明(表2),O3處理、樹種之間以及O3和樹種的交互作用對(duì)氣孔密度的影響均有顯著性差異；采樣時(shí)期對(duì)氣孔密度的影響不顯著(P=0.1498),但采樣時(shí)期與樹種的交互作用對(duì)氣孔密度的影響是顯著的。具體來講,隨著O3濃度的升高,6種樹種的氣孔密度在兩個(gè)時(shí)期均呈顯著的線性下降趨勢(shì)(P<0.05),且8月份的下降斜率均大于9月份(表3和圖1)。與CF相比,O3濃度升高至NF20已引起6種樹的氣孔密度普遍降低(圖2)；NF60處理下楊樹546、楊樹107、白蠟、法桐、洋槐和國(guó)槐的氣孔密度與CF相比在生長(zhǎng)旺期(8月份)分別降低26%、41%、39%、55%、43%和22%,在生長(zhǎng)末期(9月份)分別降低14%、31%、35%、50%、40%和14%(圖2)。其中,楊樹546和法桐的氣孔密度在兩個(gè)采樣時(shí)期呈現(xiàn)顯著差異(P<0.001),在未受O3熏蒸的CF中,楊樹546和法桐的氣孔密度隨著生長(zhǎng)進(jìn)入末期分別下降17%和14%(圖2)；并且楊樹546出現(xiàn)O3與采樣時(shí)期的交互作用。

表2O3處理、采樣時(shí)間與樹種的交互作用對(duì)氣孔密度、氣孔開度和氣孔大小的影響的方差分析

Table2ANOVAresults(Pvalues)formaineffectsandinteractionsofO3treatment(CF, NF, NF20, NF40, and NF60),samplingtime(August and September)andsixtreespeciesonstomataldensity,stomatalapertureandstomatalsize

變量 Parameters氣孔密度Stomatal density氣孔開度Stomatal aperture氣孔大小Stomatal sizeO3 Ozone<0.0001<0.0001<0.0001采樣時(shí)期 Sampling time0.1498<0.0001<0.0001樹種Tree species<0.0001<0.0001<0.0001O3×采樣時(shí)期Ozone×sampling time0.20390.06600.0070O3×樹種Ozone×tree species<0.0001<0.0001<0.0001采樣時(shí)期×樹種Sampling time×tree species<0.0001<0.0001<0.0001O3×采樣時(shí)期×樹種Ozone×sampling time×tree species0.40660.0519<0.0001

2.2 O3濃度升高對(duì)氣孔開度的影響

方差分析結(jié)果表明,O3處理、樹種之間、采樣時(shí)期均顯著影響了氣孔開度,其中O3處理與樹種之間、采樣時(shí)期與樹種之間表現(xiàn)出顯著的交互作用(表2)。隨著O3濃度的升高,6種樹種的氣孔開度在兩個(gè)時(shí)期均呈現(xiàn)顯著的線性下降趨勢(shì)(P<0.05),且8月份的下降斜率大于9月份的(表3和圖3)。具體來看(圖4),O3處理顯著影響了每種樹種的氣孔開度,但采樣時(shí)期僅對(duì)楊樹546和白蠟呈顯著影響；并且在白蠟、法桐和國(guó)槐中出現(xiàn)O3與采樣時(shí)期的交互作用。與CF相比,O3濃度為環(huán)境濃度(即NF,AOT40在8.3—11.4μmol/mol之間,表1)時(shí),白蠟、法桐、洋槐和國(guó)槐的氣孔開度已顯著降低(圖2)。

2.3 O3濃度升高對(duì)氣孔大小的影響

三因素方差分析結(jié)果表明(表2),O3處理、樹種之間、采樣時(shí)期以及兩兩交互作用、三個(gè)因素交互作用對(duì)6種樹種的氣孔大小均有極顯著的影響(P<0.01,表2)。從圖5中可知,隨著O3濃度升高,所有6種樹種在2個(gè)時(shí)期的氣孔大小均呈現(xiàn)顯著線性下降(除了楊樹兩個(gè)基因型品種546和107在9月份的趨勢(shì)不顯著外),同樣8月份的下降斜率均大于9月份(表3和圖5)。O3處理、采樣時(shí)期以及二者的交互作用對(duì)楊樹546、白蠟、法桐和國(guó)槐的氣孔大小有顯著影響,但楊樹107和洋槐的氣孔大小僅受O3處理的影響(圖6)。

表3 不同O3處理下6種樹種2個(gè)時(shí)期的氣孔密度、氣孔開度、氣孔大小與AOT40的線性關(guān)系

圖1 不同O3處理下6種樹種在2個(gè)采樣時(shí)期的氣孔密度的變化Fig.1 Ozone exposure-response relationships with stomatal density for six tree species at two sampling times (August and September)

圖2 不同O3處理下6種樹種2個(gè)采樣時(shí)期的氣孔密度差異比較Fig.2 The six tree species comparison of stomatal density (mean±S.D., n=3) across different O3 treatments at two sampling times (August and September)不同小寫字母間表示差異顯著

圖3 不同O3處理下6種樹種在2個(gè)采樣時(shí)期的氣孔開度的變化Fig.3 Ozone exposure-response relationships with stomatal aperture for six tree species at two sampling times (August and September)

圖4 不同O3處理下6種樹種2個(gè)采樣時(shí)期的氣孔開度差異比較Fig.4 The six tree species comparison of stomatal aperture (mean±S.D., n=3) across different O3 treatments at two sampling times (August and September). Different lowercase letters denote significant difference among the bars in each panel

圖5 不同O3處理下6種樹種在2個(gè)采樣時(shí)期氣孔大小的變化Fig.5 Ozone exposure-response relationships with stomatal size for six tree species at two sampling times (August and September)

圖6 不同O3處理下6種樹種2個(gè)采樣時(shí)期氣孔大小的差異比較Fig.6 The six tree species comparison of stomatal size (mean±S.D., n=3) across different O3 treatments at two sampling times (August and September). Different lowercase letters denote significant difference among the bars in each panel

3 討論

植物葉片上的氣孔數(shù)量、大小和功能等屬性對(duì)環(huán)境變化十分敏感[29- 30]。O3對(duì)植物的影響是一個(gè)累積的過程,隨著時(shí)間的延長(zhǎng),O3脅迫對(duì)葉片氣孔屬性的損害會(huì)不斷加深[31- 32]。本研究發(fā)現(xiàn),隨著O3濃度的不斷增加和熏蒸時(shí)間的持續(xù),6種樹種葉片氣孔密度、開度和大小都不同程度的顯著減小。這些都表明O3脅迫對(duì)這些樹種葉片氣孔影響具有累積效應(yīng)。當(dāng)O3累積吸收量超過自身解毒和修復(fù)能力時(shí),氣孔結(jié)構(gòu)受損,發(fā)育受到影響[25],導(dǎo)致氣孔密度、開度和大小呈現(xiàn)顯著變化。氣孔數(shù)變少、氣孔變小,開張度變小,很可能是植物為了適應(yīng)不斷升高的O3環(huán)境而迅速做出的反應(yīng)(馴化作用)。雖然這些變化不利于植物吸收CO2進(jìn)行光合作用固定生物量,但這可能是植物在脅迫環(huán)境下把更多資源分配給維持性存活而不是生長(zhǎng)性存活的一種權(quán)衡(trade-off)策略,也即在高濃度O3脅迫下抑制了氣孔的生長(zhǎng)發(fā)育。

本研究中O3脅迫減少氣孔密度,與前人的研究結(jié)果不一致。徐文鐸等[22]對(duì)4年生銀杏葉片研究發(fā)現(xiàn),高濃度O3使氣孔密度下降,但未達(dá)到顯著水平(P>0.05)；李梅等[24]發(fā)現(xiàn)高濃度O3熏蒸下1年生水蠟葉片的氣孔密度變化不顯著；文志等[23]發(fā)現(xiàn)O3脅迫62天和110天,1年生元寶楓葉片的氣孔密度都顯著增加。付偉等[25]發(fā)現(xiàn)高濃度O3熏蒸顯著增加了4年生黃檗葉片的氣孔密度。相關(guān)研究結(jié)果存在較大差異,可能與受O3熏蒸的時(shí)間和濃度影響、樹種、以及葉發(fā)育階段有關(guān)。在葉片形態(tài)學(xué)上,氣孔產(chǎn)生和形成是在葉片發(fā)育早期[33- 34],O3濃度變化對(duì)處于發(fā)育早起葉片的氣孔數(shù)量影響較大,成年葉片氣孔數(shù)目除受環(huán)境因素影響外,還受遺傳基因的控制。目前以上研究都是基于人工控制OTC中的短期觀測(cè)O3增加對(duì)氣孔參數(shù)的影響,這不一定反應(yīng)植物在長(zhǎng)期進(jìn)化過程中對(duì)環(huán)境脅迫的適應(yīng)。對(duì)于4年生的銀杏成年葉片來說,即使成年葉片氣孔數(shù)量基本形成,O3脅迫仍在一定程度上使氣孔數(shù)量和密度有所下降[22],這可能是由于高濃度O3傷害了氣孔器細(xì)胞,也破壞表皮細(xì)胞的發(fā)生,因此氣孔數(shù)量和密度有所降低[22]。因此,開展長(zhǎng)期的實(shí)驗(yàn)有助于研究葉片生長(zhǎng)發(fā)育階段下氣孔參數(shù)對(duì)O3濃度增加的響應(yīng)這一較復(fù)雜的問題。

氣孔大小、氣孔開度決定了進(jìn)入植物體內(nèi)O3的量。氣孔越大,開度越大,O3進(jìn)入量越多。本研究中O3脅迫減少了氣孔開度和氣孔大小,與前人的研究結(jié)果一致[11, 22- 25, 35]。但也有結(jié)果發(fā)現(xiàn),由于O3導(dǎo)致氣孔反應(yīng)滯后(sluggishness)甚至失靈[21, 35- 36],使氣孔面對(duì)脅迫時(shí)無法快速做出應(yīng)激性的關(guān)閉反應(yīng)來減少O3攝入量,從而降低O3對(duì)自身的傷害。這與不同樹種對(duì)O3敏感性差異有關(guān),相對(duì)敏感的樹種很可能在相同的O3濃度下自身的解毒防御系統(tǒng)已經(jīng)崩潰,保衛(wèi)細(xì)胞作用失效從而使氣孔無法關(guān)閉,增大了氣孔開度。本研究的6種樹種在持續(xù)的不同O3濃度熏蒸下氣孔開度都減小,并隨著O3劑量的增加,減小增加,這可能是氣孔對(duì)環(huán)境做出的一種適應(yīng)調(diào)節(jié)反應(yīng)。

本研究發(fā)現(xiàn)O3處理與樹種之間、采樣時(shí)期與樹種之間在氣孔密度、開度和大小屬性間都存在顯著交互效應(yīng)(表2)。由于每種樹種固有的氣孔屬性本身有很大差別(見圖2、4、6中CF的對(duì)照數(shù)值)；并且隨著生長(zhǎng)期的持續(xù),固有氣孔屬性會(huì)隨著環(huán)境變化和內(nèi)在物候期而發(fā)生變化[29- 30],因此表現(xiàn)出樹種與采樣時(shí)期之間存在顯著交互效應(yīng)。另外,樹種對(duì)O3的敏感性具有種屬特異性[7, 9, 26],越敏感的樹種會(huì)越早做出生理生態(tài)(光合生理、水分生理等)的響應(yīng)變化[26]；本研究設(shè)置了5種O3濃度梯度,從每種樹種的AOT40值來看(表1),這6種樹種的敏感性排序從高到低為：楊樹546>楊樹107>法桐>白蠟>洋槐>國(guó)槐,這與其他研究中的敏感性結(jié)果一致[9, 12, 16]。因此,隨著O3濃度的升高,不同敏感性樹種的氣孔屬性做出響應(yīng)的時(shí)間和O3劑量不同,表現(xiàn)出樹種與O3處理(不同濃度)之間存在顯著交互效應(yīng)。

本研究?jī)H對(duì)氣孔特征指標(biāo)與O3劑量之間的線性關(guān)系進(jìn)行了分析, 進(jìn)一步對(duì)其他生理指標(biāo)(如光合作用特征和抗氧化系統(tǒng)等)進(jìn)行定量化分析將是今后研究的重點(diǎn)[26]。植物既可以短期調(diào)節(jié)氣孔開閉又可以長(zhǎng)期控制氣孔發(fā)育來適應(yīng)環(huán)境的變化。O3對(duì)植物的累積性以及植物對(duì)O3脅迫的適應(yīng)性,使得評(píng)價(jià)O3脅迫對(duì)不同敏感性植物以及植物不同發(fā)育階段的氣孔運(yùn)動(dòng),從而最終評(píng)價(jià)O3對(duì)不同植物的影響和危害變得具有挑戰(zhàn)性,也將是今后開展多年長(zhǎng)期實(shí)驗(yàn)研究的必要性所在。

4 結(jié)論

(1)6種城市綠化樹木的氣孔特征指標(biāo)對(duì)O3濃度升高極為敏感,隨著O3濃度的升高,6種樹種的氣孔密度、開度和大小出現(xiàn)不同程度的顯著降低。

(2)樹種與O3處理之間,以及樹種與采樣時(shí)間之間產(chǎn)生顯著性的交互影響,可能與不同樹種之間的氣孔特征存在顯著差異,以及生長(zhǎng)旺期與末期葉片不同的氣孔發(fā)育階段有關(guān)。

(3)本研究為深入研究環(huán)境變化背景下城市植物的葉組織結(jié)構(gòu)和功能在適應(yīng)地表O3污染方面提供理論證據(jù)。

參考文獻(xiàn)(References):

[1] Wang T, Xue L K, Brimblecombe P, Lam Y F, Li L, Zhang L. Ozone pollution in China: a review of concentrations, meteorological influences, chemical precursors, and effects. Science of the Total Environment, 2017, 575: 1582- 1596.

[2] Liu X J, Zhang Y, Han W X, Tang A H, Shen J L, Cui Z L, Vitousek P, Erisman J W, Goulding K, Christie P, Fangmeier A, Zhang F S. Enhanced nitrogen deposition over China. Nature, 2013, 494(7438): 459- 462.

[3] Chen W, Tang H Z, Zhao H M. Diurnal, weekly and monthly spatial variations of air pollutants and air quality of Beijing. Atmospheric Environment, 2015, 119: 21- 34.

[4] Feng Z Z, Hu E Z, Wang X K, Jiang L J, Liu X J. Ground-level O3pollution and its impacts on food crops in China: a review. Environmental Pollution, 2015, 199: 42- 48.

[5] Novak K, Skelly J M, Schaub M, Kr?uchi N, Hug C, Landolt W, Bleuler P. Ozone air pollution and foliar injury development on native plants of Switzerland. Environmental Pollution, 2003, 125(1): 41- 52.

[6] Ren W, Tian H Q, Tao B, Chappelka A, Sun G, Lu C Q, Liu M L, Chen G S, Xu X F. Impacts of tropospheric ozone and climate change on net primary productivity and net carbon exchange of China′s forest ecosystems. Global Ecology and Biogeography, 2011, 20(3): 391- 406.

[7] Zhang W W, Feng Z Z, Wang X K, Niu J F. Responses of native broadleaved woody species to elevated ozone in subtropical China. Environmental Pollution, 2012, 163: 149- 157.

[8] Feng Z Z, Sun J S, Wan W X, Hu E Z, Calatayud V. Evidence of widespread ozone-induced visible injury on plants in Beijing, China. Environmental Pollution, 2014, 193: 296- 301.

[9] Li P, Calatayud V, Gao F, Uddling J, Feng Z Z. Differences in ozone sensitivity among woody species are related to leaf morphology and antioxidant levels. Tree Physiology, 2016, 36(9): 1105- 1116.

[10] Feng Z Z, Zeng H Q, Wang X K, Zheng Q W, Feng Z W. Sensitivity ofMetasequoiaglyptostroboidesto ozone stress. Photosynthetica, 2008, 46(3): 463- 465.

[11] Wittig V E, Ainsworth E A, Naidu S L, Karnosky D F, Long S P. Quantifying the impact of current and future tropospheric ozone on tree biomass, growth, physiology and biochemistry: a quantitative meta-analysis. Global Change Biology, 2009, 15(2): 396- 424.

[12] Gao F, Calatayud V, García-Breijo F, Reig-Armiana J, Feng Z Z. Effects of elevated ozone on physiological, anatomical and ultrastructural characteristics of four common urban tree species in China. Ecological Indicators, 2016, 67: 367- 379.

[13] Nikula S, Percy K, Oksanen E, Holopainen T, Manninen S. Effects of elevated ozone on growth and foliar traits of European and hybrid aspen. Boreal Environment Research, 2009, 14(SA): 29- 47.

[14] Sun G E, McLaughlin S B, Porter J H, Uddling J, Mulholland P J, Adams M B, Pederson N. Interactive influences of ozone and climate on streamflow of forested watersheds. Global Change Biology, 2012, 18(11): 3395- 3409.

[15] Uddling J, Teclaw R M, Pregitzer K S, Ellsworth D S. Leaf and canopy conductance in aspen and aspen-birch forests under free-air enrichment of carbon dioxide and ozone. Tree Physiology, 2009, 29(11): 1367- 1380.

[16] Hu E Z, Gao F, Xin Y, Jia H X, Li K H, Hu J J, Feng Z Z. Concentration-and flux-based ozone dose-response relationships for five poplar clones grown in North China. Environmental Pollution, 2015, 207: 21- 30.

[17] Krupa S, McGrath M T, Andersen C P, Booker F L, Burkey K O, Chappelka A, Chevone B I, Pell E J, Zilinskas B A. Ambient ozone and plant health. Plant Disease, 2001, 85(1): 4- 12.

[18] Sitch S, Cox PM, Collins WJ, Huntingford C. Indirect radiative forcing of climate change through ozone effects on the land-carbon sink. Nature, 2007, 448(7155): 791- 794.

[19] Ainsworth E A, Yendrek C R, Sitch S, Collins W J, Emberson L D. The effects of tropospheric ozone on net primary productivity and implications for climate change. Annual Review of Plant Biology, 2013, 63: 637- 661.

[20] 任巍, 田漢勤. 臭氧污染與陸地生態(tài)系統(tǒng)生產(chǎn)力. 植物生態(tài)學(xué)報(bào), 2007, 31(2): 219- 230.

[21] Hoshika Y, Omasa K, Paoletti E. Both ozone exposure and soil water stress are able to induce stomatal sluggishness. Environmental and Experimental Botany, 2013, 88: 19- 23.

[22] 徐文鐸, 齊淑艷, 何興元, 陳瑋, 趙桂玲, 周園. 大氣中CO2、O3濃度升高對(duì)銀杏成年葉片氣孔數(shù)量特征的影響. 生態(tài)學(xué)雜志, 2008, 27(7): 1059- 1063.

[23] 文志, 王麗, 王效科, 李麗, 崔健. O3和干旱脅迫對(duì)元寶楓葉片氣孔特征的復(fù)合影響. 生態(tài)學(xué)雜志, 2014, 33(3): 560- 566.

[24] 李梅, 徐勝, 張恒慶, 何興元, 陳瑋, 趙詣, 李波. 高濃度O3對(duì)水蠟葉表皮氣孔及其葉組織結(jié)構(gòu)特征的影響. 生態(tài)學(xué)雜志, 2014, 33(1): 53- 58.

[25] 付偉, 鄧?yán)蛱m, 徐勝, 何興元, 陳瑋, 趙詣, 張娜, 蘇麗麗. 臭氧對(duì)黃檗幼苗葉片可見傷害及氣孔特征的影響. 東北林業(yè)大學(xué)學(xué)報(bào), 2015, 43(2): 14- 18.

[26] 辛月, 高峰, 馮兆忠. 不同基因型楊樹的光合特征與臭氧劑量的響應(yīng)關(guān)系. 環(huán)境科學(xué), 2016, 37(6): 2359- 2367.

[27] 劉明智, 努爾巴衣·阿布都沙力克, 潘曉玲. 指甲油涂抹撕取法制取植物葉氣孔裝片. 生物學(xué)通報(bào), 2005, 40(10): 44- 44.

[28] Tricker P J, Trewin H, Kull O, Clarkson G J J, Eensalu E, Tallis M J, Colella A, Doncaster C P, Sabatti M, Taylor G. Stomatal conductance and not stomatal density determines the long-term reduction in leaf transpiration of poplar in elevated CO2. Oecologia, 2005, 143(4): 652- 660.

[29] 鄭淑霞, 上官周平. 近70年黃土高原3種植物葉片氣孔特征參數(shù)比較. 植物資源與環(huán)境學(xué)報(bào), 2005, 14(1): 1- 5.

[30] 鄭淑霞, 上官周平. 遼東櫟葉片氣孔密度及δ13C值的時(shí)空變異. 林業(yè)科學(xué), 2005, 41(2): 30- 36.

[31] Asensi-Fabado A, García-Breijo F J, Reig-Armiana J. Ozone-induced reductions in below-ground biomass: an anatomical approach in potato. Plant, Cell & Environment, 2010, 33(7): 1070- 1083.

[32] Giacomo B, Forino L M C, Tagliasacchi A M, Bernardi R, Durante M. Ozone damage and tolerance in leaves of two poplar genotypes. Caryologia, 2010, 63(4): 422- 434.

[33] Gay A P, Hurd R G. The influence of light on stomatal density in the tomato. New Phytologist, 1975, 75(1): 37- 46.

[34] Schoch P G, Zinsou C, Sibi M. Dependence of the stomatal index on environmental factors during stomatal differentiation in leaves ofVignasinensisL. I. Effect of light intensity. Journal of Experimental Botany, 1980, 31(5): 1211- 1216.

[35] Hoshika Y, Carriero G, Feng Z Z, Zhang Y L, Paoletti E. Determinants of stomatal sluggishness in ozone-exposed deciduous tree species. Science of the Total Environment, 2014, 481: 453- 458.

[36] Paoletti E, Grulke N E. Ozone exposure and stomatal sluggishness in different plant physiognomic classes. Environmental Pollution, 2010, 158: 2664- 2671.