基于樹干液流技術(shù)的北京市刺槐冠層吸收臭氧特征研究

王 華，歐陽志云，任玉芬，張紅星，王效科，郜世奇，關(guān)云飛，高付元

(1.中國科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室，北京 100085;2.北京市農(nóng)林科學(xué)院林業(yè)果樹研究所，北京 100093;3.北京教學(xué)植物園，北京 100061)

近地層臭氧(O3)作為影響全球氣候變化的重要因素和危害人類健康、動(dòng)植物生長的二次污染物，受到廣泛關(guān)注［1－3］。隨著城市化進(jìn)程的推進(jìn)，近地層臭氧濃度迅速增加［4］。樹木是城市生態(tài)系統(tǒng)的重要組成部分，可通過氣孔途徑和非氣孔途徑有效去除大氣臭氧，在凈化空氣方面發(fā)揮著重要作用［5］。其中，通過氣孔途徑吸收是林木去除臭氧的主要途徑，并與臭氧引起林木傷害密切相關(guān)［6］。因此，城市樹木吸收臭氧方面的研究已成為目前植物生理生態(tài)學(xué)和城市生態(tài)學(xué)研究的熱點(diǎn)。然而，現(xiàn)有的大部分研究主要集中在天然林吸收臭氧的影響機(jī)制［7－10］，城市樹木吸收臭氧的環(huán)境影響機(jī)制研究較少。

葉片尺度和生態(tài)系統(tǒng)尺度上森林去除臭氧的研究均有報(bào)道，但是冠層尺度上森林去除臭氧的研究很少。樹干液流技術(shù)已經(jīng)廣泛應(yīng)用于森林水分利用研究，由于水分蒸騰和吸收臭氧通過氣孔行為耦合，故基于樹干液流技術(shù)能夠衡量森林冠層吸收臭氧量，這種方法考慮了邊界層的影響，適用于異質(zhì)性景觀和山地景觀［8，11］。因此，結(jié)合城市環(huán)境的異質(zhì)化特征，本研究采用基于樹干液流的方法衡量北京市典型綠化樹種刺槐冠層吸收臭氧規(guī)律。

我國首都北京是世界上面積最大、歷史最悠久的城市之一，是現(xiàn)代城市化的典型代表。由于大范圍建設(shè)和機(jī)動(dòng)車數(shù)量迅速增加，臭氧污染問題日趨嚴(yán)峻［12］。改善空氣質(zhì)量一直是該城市面臨的重大難題。本文以北京市典型綠化樹種刺槐為研究對(duì)象，通過測(cè)定樹干液流密度，結(jié)合環(huán)境因子監(jiān)測(cè)，探索冠層尺度上刺槐吸收臭氧特征及其與環(huán)境因子之間的關(guān)系。研究結(jié)果和方法可為明確環(huán)境如何影響樹木吸收臭氧功能、空氣質(zhì)量提高、以及城市規(guī)劃提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究地點(diǎn)和研究對(duì)象

研究地點(diǎn)位于北京市南二環(huán)的北京教學(xué)植物園(116°25'37″— 116°25'50″E，39°52'20″— 39°52'28″N)，占地面積116，500m2，周邊商業(yè)區(qū)和住宅區(qū)密集，行人和機(jī)動(dòng)車流量較大。該地區(qū)屬于溫帶半濕潤大陸性季風(fēng)氣候，年均氣溫11—12℃，年均降水量約為586 mm，超過70%發(fā)生在6—8月。

研究對(duì)象為北京市典型綠化樹種刺槐，選擇健康、生長環(huán)境一致的單株作為研究對(duì)象，2009年樣樹樹形測(cè)定結(jié)果為:平均胸徑為(38.43±2.43)cm，平均樹高(12.87±0.20)m，平均邊材面積為(129.12±12.11)cm2(表 1)。

表1 測(cè)定樹干液流樣樹的樹形特征Table 1 Characteristics of the sampled trees for sap flow measurements

1.2 樹形測(cè)定

樹冠垂直投影面積采用測(cè)定東西和南北冠幅，結(jié)合橢圓面積計(jì)算公式確定。為避免傷害樣樹，于北京教學(xué)植物園和北京林業(yè)大學(xué)鷲峰教學(xué)林場(chǎng)選取各個(gè)徑級(jí)的23株樹木，采用生長錐鉆取木芯，再用直尺測(cè)定邊材厚度，建立邊材面積與胸徑之間的關(guān)系式，進(jìn)而反推測(cè)定樣樹邊材面積(圖1)。

1.3 樹干液流密度測(cè)定

2009年7月1日到2009年10月31日采用熱消散傳感器(Dynamax，Houston，TX，USA)連續(xù)測(cè)定單位邊材面積上的瞬時(shí)樹干液流密度(Js，g H2O cm－2s－1)。傳感器的工作原理、安裝、液流密度計(jì)算以及整樹蒸騰計(jì)算方法參見文獻(xiàn)［13－14］。將傳感器饋線與數(shù)據(jù)采集器(CR1000 Dynamax，USA)連接，組成完整的樹干液流測(cè)定系統(tǒng)，每10 min進(jìn)行平均并儲(chǔ)存數(shù)據(jù)，將筆記本電腦與數(shù)據(jù)采集器連接，定期采集數(shù)據(jù)。

在樹干液流密度點(diǎn)的測(cè)定值尺度推廣到整樹尺度的過程中，忽視樹干液流密度的軸向和徑向變異，將引起結(jié)果的誤差［15－16］。為明確樹干液流密度的軸向變異，Herbst等建議小樹采用1個(gè)或者2個(gè)傳感器，大樹采用4個(gè)傳感器［17］，因而，在每株刺槐樣樹東南西北面各安裝一個(gè)傳感器。由于刺槐樣樹邊材厚度為2 cm左右，故樹干液流密度徑向變異不明顯，因而，采用長度為30mm的傳感器。

1.4 環(huán)境因素監(jiān)測(cè)

圖1 刺槐邊材面積與胸徑之間的關(guān)系Fig.1 The relationship between sapwood area(A s)and diameter at breast height(DBH)for R.pseudoacacia

氣象站位于北京教學(xué)植物園中較為開闊的一塊空地上，避免樹木、建筑物和其它障礙物的影響，同時(shí)距離研究對(duì)象近?？諝鉁囟葷穸葌鞲衅?HMP45C，Vaisala Inc.，Helsinki，F(xiàn)inland)，風(fēng)速風(fēng)向傳感器(034B，Met One Instruments，Grants Pass，Oregon，USA)安裝在高度為10 m的標(biāo)準(zhǔn)桅桿上。日照強(qiáng)度計(jì)(CMP－11，Kipp and Zonen，Delft，Netherlands)安裝在高度為1.5 m的標(biāo)準(zhǔn)桅桿上測(cè)定太陽輻射。降雨由雨量計(jì)(TE525MM，Campbell Scientific Inc.，Logan，UT)在地面以上2 m以下高度測(cè)定。土壤溫度傳感器(109，Campbell Scientific Inc.，Logan，UT)埋設(shè)深度為10、30、50 cm和 80 cm。土壤含水量傳感器(ECH2O，Decagon Devices Inc.，Pullman，WA，USA)埋設(shè)位于研究樹木之間，深度為30 cm。水汽壓虧缺是由大氣溫度和空氣相對(duì)濕度計(jì)算得到［18］。上述氣象和土壤參數(shù)監(jiān)測(cè)與樹干液流測(cè)定同步。

大氣臭氧濃度采用 TEI Model 49i氣體分析儀(Thermo Environmental Instruments Inc.，F(xiàn)ranklin，MA，USA)監(jiān)測(cè)。臭氧監(jiān)測(cè)儀器距離研究對(duì)象近，每小時(shí)記錄一個(gè)濃度值。

1.5 冠層導(dǎo)度計(jì)算

城市樹木冠層開放，通風(fēng)好，空氣動(dòng)力學(xué)粗糙，表明冠層表面與大氣之間耦合度強(qiáng)［19］。因此，研究樹種刺槐的整樹冠層導(dǎo)度(Gc，mmol H2O·m－2·s－1)計(jì)算公式如下［20］:

式中，Ec是冠層蒸騰，將樹干液流密度乘以邊材面積再除以冠層投影面積得到［13］;ρ是水的密度(998 kg/m3);GV是水蒸氣的通用氣體常數(shù)(0.462 m3kPa K－1kg－1);Ta是大氣溫度;D是水汽壓虧缺(kPa)。參考Granier等人的建議［21］，基于如下原則選擇樹干液流密度數(shù)據(jù)計(jì)算冠層導(dǎo)度:1)去除降雨或者降雨以后2h的數(shù)據(jù)，從而避免蒸散和樹木蒸騰之間不一致;2)去除總輻射、水汽壓虧缺、或者林段蒸騰非常低的數(shù)據(jù)(小于最大值的5%)，因?yàn)樵谶@些條件下計(jì)算冠層導(dǎo)度的不確定性較大。

Gc對(duì)水汽壓虧缺(D)和總輻射(Rs)的響應(yīng)規(guī)律通過邊界線分析得到［22－23］，采用不同D(步長0.2 kPa)或者Rs(步長50 W/m2)等級(jí)下3株樣樹的最大冠層導(dǎo)度。采用如下公式描述Gc響應(yīng)D變化［19］:

式中，參數(shù)－m為回歸方程斜率，量化了氣孔導(dǎo)度響應(yīng)D變化的敏感度;參數(shù)b為回歸方程截距，為D=1 kPa的參比冠層導(dǎo)度。采用Sigma－Plot 10.0軟件(Systat Software Inc.，San Jose，California)的最小二乘回歸確定m和b。上述Gc和ln(D)之間的關(guān)系式為各樣樹之間的比較提供了簡便的基準(zhǔn)。

1.6 整樹冠層吸收臭氧速率計(jì)算

葉肉中臭氧濃度近乎為0［24］，因此，可以根據(jù)如下通量公式計(jì)算冠層吸收臭氧速率:

式中，F(xiàn)O3是單位冠層投影面積的冠層吸收臭氧速率，［O3］是大氣臭氧濃度，Gc是冠層對(duì)水蒸氣的導(dǎo)度，0.613是轉(zhuǎn)換系數(shù)，代表大氣中臭氧和水蒸氣擴(kuò)散系數(shù)比［7，25］。

2 結(jié)果

2.1 刺槐吸收臭氧的日變化規(guī)律

夏季晴天天氣條件下，大氣臭氧濃度隨著總輻射增大逐漸上升，于15:00左右達(dá)到峰值，之后開始下降;刺槐冠層對(duì)臭氧的導(dǎo)度上午達(dá)到峰值，之后迅速下降;刺槐吸收臭氧速率呈單峰曲線，于15:00左右達(dá)到最大值，峰值范圍較窄，盡管上午冠層對(duì)臭氧的導(dǎo)度最高，但由于臭氧濃度很低，故上午刺槐吸收臭氧速率并沒有達(dá)到最大值;刺槐吸收臭氧累積量不斷增加，中午前后增加最明顯，夜間隨著臭氧濃度降低和氣孔關(guān)閉，其吸收臭氧累積量變化不大(圖2)。

秋季晴天天氣條件下，總輻射和水汽壓虧缺明顯減少，進(jìn)而大氣臭氧濃度也降低;刺槐冠層對(duì)臭氧的導(dǎo)度明顯降低，且下午刺槐冠層對(duì)臭氧的導(dǎo)度響應(yīng)水汽壓虧缺下降而減少的幅度較小;刺槐吸收臭氧速率呈單峰曲線，峰值范圍較寬;刺槐冠層吸收臭氧累積量明顯減少，中午前后增加最明顯，且從傍晚開始保持較為穩(wěn)定的狀態(tài)(圖2)。

2.2 刺槐吸收臭氧的季節(jié)變化規(guī)律

圖2 夏(2009年7月18日)秋(2009年10月4日)季節(jié)晴天天氣下總輻射(R s)、水汽壓虧缺(D)、大氣臭氧濃度(O3)、刺槐冠層對(duì)臭氧的導(dǎo)度(G O3)、吸收臭氧速率的日變化(F O3)以及吸收臭氧累積量(AFst)Fig.2 Diurnal course of total radiation(R s)，vapour pressure deficit(D);ambient O3 concentration(middle)and canopy conductance for ozone(G O3);canopy O3 uptake rate(F O3)and cumulative O 3 uptake(AFst)in R.pseudoacacia on clear days in the summer(18th July 2009，left)and the autumn(4th October 2009，right)

大氣臭氧濃度、刺槐冠層對(duì)臭氧的導(dǎo)度、吸收臭氧速率日均值以及吸收臭氧累積量的季節(jié)變化規(guī)律如圖3所示。臭氧濃度日均值季節(jié)變化明顯，夏季出現(xiàn)最大值146.4μg/m3，秋季出現(xiàn)最小值1.5μg/m3。夏秋季節(jié)刺槐冠層對(duì)臭氧的導(dǎo)度和吸收臭氧速率日均值波動(dòng)較大，陰雨天氣條件下冠層對(duì)臭氧的導(dǎo)度和吸收臭氧速率較低，而晴天天氣條件下冠層對(duì)臭氧的導(dǎo)度和吸收臭氧速率較高;然而，從總體趨勢(shì)上看，刺槐冠層對(duì)臭氧的導(dǎo)度和吸收臭氧速率從夏季到秋季逐漸降低，可見，刺槐吸收臭氧速率的日變化與季節(jié)變化規(guī)律交疊。刺槐吸收臭氧累積量夏季顯著增加，秋季略有增加。

2.3 總輻射和水汽壓虧缺對(duì)刺槐吸收臭氧的影響

刺槐冠層吸收臭氧速率日變化規(guī)律與大氣臭氧濃度日變化規(guī)律相似(圖2)。在一定的臭氧濃度下，刺槐冠層吸收臭氧速率的變化主要由冠層對(duì)臭氧的導(dǎo)度調(diào)控，進(jìn)而受水汽壓虧缺和總輻射影響(圖4)。3株刺槐樣樹冠層對(duì)臭氧的導(dǎo)度接近，且對(duì)水汽壓虧缺和總輻射的響應(yīng)方式相似(圖4)。

水汽壓虧缺較高的條件下，臭氧濃度較高，但刺槐冠層對(duì)臭氧的導(dǎo)度降低，進(jìn)而導(dǎo)致冠層吸收臭氧速率降低;水汽壓虧缺較低的條件下，刺槐冠層對(duì)臭氧的導(dǎo)度很高，故在臭氧濃度較低的條件下，冠層吸收臭氧速率較高;水汽壓虧缺很小的條件下，刺槐冠層吸收臭氧速率略降，可能是由于早晨光合作用較弱，氣孔開度較小(圖4)。

隨著總輻射上升，臭氧濃度增加，當(dāng)總輻射大于800 W/m2，臭氧濃度顯著降低。當(dāng)總輻射小于600 W/m2，刺槐冠層對(duì)臭氧的導(dǎo)度變化很小;當(dāng)總輻射大于600 W/m2，刺槐冠層對(duì)臭氧的導(dǎo)度迅速下降，這個(gè)時(shí)候氣孔可能達(dá)到了光飽和。因此，受臭氧濃度和冠層對(duì)臭氧的導(dǎo)度的雙重影響，刺槐冠層吸收臭氧速率呈不對(duì)稱單峰格局:當(dāng)總輻射小于800 W/m2，隨總輻射增強(qiáng)，刺槐冠層吸收臭氧速率略有增加;當(dāng)總輻射在800 W/m2左右時(shí)，冠層吸收臭氧速率達(dá)到峰值，之后迅速下降(圖4)。

3 討論

3.1 刺槐整樹冠層吸收臭氧特征

北京市典型綠化樹種刺槐整樹冠層吸收臭氧速率夏季高于秋季，與夏季大氣溫度、輻射、臭氧濃度等較高有關(guān)，這與高海拔地區(qū)成年挪威云杉(Picea abies)、瑞士五葉松(Pinus cembra)和歐洲落葉松(Larix decidua)吸收臭氧速率的季節(jié)變化規(guī)律一致［8］。與上述研究不一致，加那利群島松林吸收臭氧速率季節(jié)變化不大，可能是如下原因?qū)е?首先，當(dāng)空氣和土壤干旱時(shí)，在濃霧對(duì)空氣濕度的影響和松林根系能獲取土壤深層水分的聯(lián)合作用下，松林吸收臭氧速率降低不顯著［26，27］。其次，由于地處低緯度和冷洋流中，該群島氣候季節(jié)變化不大。最后，群島的經(jīng)濟(jì)長期以農(nóng)業(yè)為主，大城市很少，故污染源少。

圖4 2009年7月1日到2009年10月31日大氣臭氧濃度(O3)、刺槐冠層對(duì)臭氧的導(dǎo)度(G O3)、吸收臭氧速率(F O3)對(duì)水汽壓虧缺(D)和總輻射(R s)的響應(yīng)規(guī)律Fig.4 Means of ambient air ozone concentration(O3)，canopy conductance for ozone(G O3)，and canopy O3 uptake rate(F O3)in relation to vapour pressure deficit(D)and total radiation(R s)from 1st July，2009 to 31th October，2009

北京市典型綠化樹種刺槐單位冠層投影面積上年吸收臭氧量約為0.16 g/m2(例如，3株樣樹單位冠層投影面積上月均吸收臭氧量0.57 mmol·m－2·月－1×6個(gè)月×48 g/mol÷1000=0.16 g/m2)，明顯低于模型得到的結(jié)果［5，28］。采用城市森林影響模型，Nowak和 Dwyer［5］總結(jié)美國城市森林年均去除污染物量為 9.3 g/m2，變化范圍為6.6—12.0 g/m2，臭氧占污染物比率為31.8%，因此其去除臭氧量約為2.96 g/m2(例如，9.3 g/m2×31.8%=2.96 g/m2)。采用同樣的模型，北京城市森林去除污染物量較高(27.5 g/m2)，臭氧占污染物比率為20%，因此其去除臭氧量約為5.50 g/m2［28］。野外觀測(cè)數(shù)據(jù)得到的結(jié)果和模型方法得到結(jié)果區(qū)別較大，可能是如下原因?qū)е隆Ｊ紫?，基于樹干液流的方法確定的是氣孔吸收臭氧量，而模型方法確定的是干沉降值，包括氣孔吸收量和樹木非蒸騰表面吸附量。其次，模型方法得到的結(jié)果為森林去除污染物量的粗略估計(jì)，需要在如下幾個(gè)方面調(diào)整參數(shù)以提高模型的準(zhǔn)確度:模型產(chǎn)生結(jié)果的不確定性分析;測(cè)定樹木個(gè)體去除污染物以改進(jìn)模型的假設(shè)－樹木處于相同層次;準(zhǔn)確模擬大氣污染物濃度及其與地點(diǎn)、城市森林的面積和樹種組成之間的關(guān)系［28］。再次，老樹較低的水力導(dǎo)度限制其光合速率［29］，研究對(duì)象刺槐樹齡超過100 a，因而推斷其氣孔導(dǎo)度明顯降低，進(jìn)而導(dǎo)致冠層吸收臭氧速率較低。最后，近年來北京市空氣質(zhì)量逐漸改善，大氣臭氧濃度有所降低。兩種方法得到結(jié)果的區(qū)別表明:評(píng)估森林受臭氧危害的風(fēng)險(xiǎn)應(yīng)考慮各樹種冠層臭氧通量。

3.2 刺槐整樹冠層吸收臭氧的環(huán)境調(diào)控

城市典型綠化樹種刺槐整樹冠層吸收臭氧特征受臭氧濃度和冠層導(dǎo)度的共同影響。本研究地點(diǎn)臭氧濃度季節(jié)變化顯著，夏季較高，秋季較低。在一定程度上，這導(dǎo)致了研究樹種刺槐冠層吸收臭氧速率的季節(jié)變化規(guī)律與臭氧濃度的變化規(guī)律相似(圖3)。在一定的臭氧濃度下，氣孔導(dǎo)度是整樹吸收臭氧速率的關(guān)鍵限制因子，進(jìn)而受微氣候、樹齡、海拔等因素影響［7，25］。與上述研究一致，刺槐整樹冠層對(duì)臭氧的導(dǎo)度主要受水汽壓虧缺和總輻射影響(圖4)。低輻射條件下刺槐保持一定的吸收臭氧速率，可能由于夜間、凌晨或者陰雨天氣條件下氣孔保持部分開放。很多研究表明某些植物夜間氣孔開放，保持一定速率的蒸騰［30］。城市環(huán)境中較高的空氣溫度與較低的平均風(fēng)速，可能誘導(dǎo)植物氣孔增大開度。因此，城市環(huán)境下刺槐吸收臭氧速率可能高于自然環(huán)境下刺槐吸收臭氧速率。

3.3 冠層尺度上城市樹木吸收臭氧研究方法

量化樹木吸收臭氧的方法包括葉片水平上的氣體交換法和冠層尺度上的渦度相關(guān)法、模型、基于樹干液流的方法。氣體交換法測(cè)定葉片吸收臭氧量，但是忽視了邊界層的影響［31］。渦度相關(guān)技術(shù)測(cè)定冠層臭氧沉降量，但是無法區(qū)分氣孔吸收臭氧與樹木表面吸附臭氧［32］。模型則需要根據(jù)當(dāng)?shù)丨h(huán)境以及植被的特點(diǎn)對(duì)模型參數(shù)化［33］?；跇涓梢毫鞯姆椒紤]了邊界層的影響，適用于異質(zhì)性景觀和山地景觀［7］。

由于城市環(huán)境突出的異質(zhì)性特征，如何選擇適合的方法研究城市樹木吸收臭氧的一個(gè)難點(diǎn)。目前，城市林木吸收臭氧速率主要采用整合植被、氣象和大氣臭氧濃度信息的模型模擬的方法［28］。此類研究能反映區(qū)域尺度上城市樹木去除臭氧速率，但是不能反映樹木直接吸收臭氧量，也無法進(jìn)行種間比較。為解決上述問題，根據(jù)樹木吸收臭氧和蒸騰作用通過氣孔行為耦合，提出基于樹干液流技術(shù)和大氣臭氧濃度監(jiān)測(cè)確定城市樹木整樹吸收臭氧速率方法，為衡量城市林木吸收臭氧及其它痕量氣體，以及評(píng)價(jià)城市樹木受臭氧危害的風(fēng)險(xiǎn)等提供了準(zhǔn)確可行的方法。

［1］Karnosky D F，Skelly JM，Percy K E，Chappelka A H.Perspectives regarding50 years of research on effects of tropospheric ozone air pollution on USforests.Environmental Pollution，2007，147(3):489－506.

［2］Sitch S，Cox P M，Collins WJ，Huntingford C.Indirect radiative forcing of climate change through ozone effects on the land－carbon sink.Nature，2007，448(7155):791－794.

［3］Selin N E，Wu S，Nam K M，Reilly J M，Paltsev S，Prinn R G，Webster M D.Global health and economic impacts of future ozone pollution.Environmental Research Letters，2009，4(4):044014.

［4］Vingarzan R.A review of surface ozone background levels and trends.Atmospheric Environment，2004，38(21):3431－3442.

［5］Nowak D J，Dwyer JF.Understanding the benefits and costs of urban forest ecosystem ∥ Kuser J E，ed.Urban and Community Forestry in the Northeast.New York:Springer Science and Business Media，2007:25－46.

［6］Landolt W，Keller T.Uptake and effects of air pollutants on woody plants.Experientia，1985，41(3):301－309.

［7］Wieser G，H?sler R，G?tz B，Koch W，Havranek WM.Role of climate，crown position，tree age and altitude in calculated ozone flux into needles of Picea abies and Pinus cembra:a synthesis.Environmental Pollution，2000，109(3):415－422.

［8］Wieser G，Matyssek R，K?stner B，Oberhuber W.Quantifying ozone uptake at the canopy level of spruce，pine and larch trees at the alpine timberline:an approach based on sap flow measurement.Environmental Pollution，2003，126(1):5－8.

［9］Nunn A J，Wieser G，Metzger U，L?w M，Wipfler P，H?berle K，Matyssek R.Exemplifying whole－plant ozone uptake in adult forest trees of contrasting species and site conditions.Environmental Pollution，2007，146(3):629－639.

［10］Braun S，Schindler C，Leuzinger S.Use of sap flow measurements tovalidatestomatal functions for mature beech(Fagussylvatica)in view of ozone uptake calculations.Environmental Pollution，2010，158(9):2954－2963.

［11］K?stner B，Matyssek R，Heilmeier H，Clausnitzer F，Nunn A J，Wieser G.Sap flow measurements as a basis for assessing trace－gas exchange of trees.Flora－Morphology，Distribution，F(xiàn)unctional Ecology of Plants，2008，203(1):14－33.

［12］Chan C K，Yao X H.Air pollution in mega cities in China.Atmospheric Environment，2008，42(1):1－42.

［13］Granier A.Evaluation of transpiration in a Douglas－fir stand by means of sap flow measurements.Tree Physiology，1987，3(4):309－320.

［14］Zhao P，Rao X Q，Ma L，Cai X A，Zeng X P.Application of Granier’s sap flow system in water use of Acacia mangium forest.Journal of Tropical and Subtropical Botany，2005，13(6):457－468.

［15］Delzon S，Sartore M，Granier A，Loustau D.Radial profiles of sap flow with increasing tree size in maritime pine.Tree Physiology，2004，24(11):1285－1293.

［16］Lu P，Müller W J，Chacko E K.Spatial variations in xylem sap flux density in the trunk of orchard－grown，mature mango trees under changing soil water conditions.Tree Physiology，2000，20(10):683－692.

［17］Herbst M，Roberts JM，Rosier P T W，Taylor M E，Gowing DJ.Edge effects and forest water use:a field study in a mixed deciduous woodland.Forest Ecology and Management，2007，250(3):176－186.

［18］Campbell G S，Norman JM.An Introduction to Environmental Biophysics.2nd ed.New York:Springer－Verlag，1998.

［19］Oren R，Sperry JS，Katul G G，Pataki D E，Ewers B E，Phillips N，Sch?fer K V R.Survey and synthesis of intra－and interspecific variation in stomatal sensitivity to vapour pressure deficit.Plant，Cell and Environment，1999，22(12):1515－1526.

［20］K?stner B M M，Schulze E D，Kelliher F M，Hollinger D Y，Byers J N，Hunt J E，McSeveny T M，Meserth R，Weir P L.Transpiration and canopy conductance in a pristine broad－leaved forest of Nothofagus:an analysis of xylem sap flow and eddy correlation measurements.Oecologia，1992，91(3):350－359.

［21］Granier A，Loustau D，Bréda N.A generic model of forest canopy conductance dependent on climate，soil water availability and leaf area index.Annals of Forest Science，2000，57(8):755－765.

［22］Jarvis P G.The interpretation of the variations in leaf water potential and stomatal conductance found in canopies in the field.Philosophical Transactions of the Royal Society B:Biological Sciences，1976，273(927):593－610.

［23］Dye P J，Olbrich B W.Estimating transpiration from 6－year－old Eucalyptus grandis trees:development of a canopy conductance model and comparison with independent sap flux measurements.Plant，Cell and Environment，1993，16(1):45－53.

［24］Laisk A，Kull O，Moldau H.Ozone concentration in leaf intercellular air spaces is close to zero.Plant Physiology，1989，90(3):1163－1167.

［25］Wieser G，Havranek WM.Environmental control of ozone uptake in Larix decidua Mill.:a comparison between different altitudes.Tree Physiology，1995，15(4):253－258.

［26］Luis V C，Jiménez M S，Morales D，Kucera J，Wieser G.Canopy transpiration of a Canary Islands pine forest.Agricultural and Forest Meteorology，2005，135(1/4):117－123.

［27］Wieser G，Luis V C，Cuevas E.Quantification of ozone uptake at the stand level in a Pinus canariensis forest in Tenerife，Canary Islands:An approach based on sap flow measurements.Environmental Pollution，2006，140(3):383－386.

［28］Yang J，McBride J，Zhou J X，Sun Z Y.The urban forest in Beijing and its role in air pollution reduction.Urban Forestry and Urban Greening，2005，3(2):65－78.

［29］Hubbard R M，Bond B J，Ryan M G.Evidence that hydraulic conductance limits photosynthesis in old Pinus ponderosa trees.Tree Physiology，1999，19(3):165－172.

［30］Daley M J，Phillips N G.Interspecific variation in nighttime transpiration and stomatal conductance in a mixed New England deciduous forest.Tree Physiology，2006，26(4):411－419.

［31］Wieser G，Havranek W M.Ozone uptake in the sun and shade crown of spruce:quantifying the physiological effects of ozone exposure.Trees，1993，7(4):227－232.

［32］Zeller K F，Nikolov N T.Quantifying simultaneous fluxes of ozone，carbon dioxide and water vapor above a subalpine forest ecosystem.Environmental Pollution，2000，107(1):1－20.

［33］Wieser G，Emberson L D.Evaluation of the stomatal conductance formulation in the EMEP ozone deposition model for Picea abies.Atmospheric Environment，2004，38(15):2339－2348.

參考文獻(xiàn):

［14］趙平，饒興權(quán)，馬玲，蔡錫安，曾小平.Granier樹干液流測(cè)定系統(tǒng)在馬占相思的水分利用研究中的應(yīng)用.熱帶亞熱帶植物學(xué)報(bào)，2005，13(6):457－468.