北京奧林匹克森林公園綠地碳交換動(dòng)態(tài)及其環(huán)境控制因子

陳文婧，李春義，何桂梅，王小平，查天山, 3，賈 昕,*

(1. 北京林業(yè)大學(xué)水土保持學(xué)院，北京 100083；2. 北京市林業(yè)碳匯工作辦公室，北京 100013；3. 中關(guān)村綠色碳匯研究院，北京 100083)

北京奧林匹克森林公園綠地碳交換動(dòng)態(tài)及其環(huán)境控制因子

陳文婧1，李春義1，何桂梅2,3，王小平2,3，查天山1, 3，賈 昕1,*

(1. 北京林業(yè)大學(xué)水土保持學(xué)院，北京 100083；2. 北京市林業(yè)碳匯工作辦公室，北京 100013；3. 中關(guān)村綠色碳匯研究院，北京 100083)

隨著城市化進(jìn)程的推進(jìn)，城市公園綠地的面積也在不斷地增加。在碳循環(huán)與氣候變化研究中，以人工植被為主要存在形態(tài)的城市綠地生態(tài)系統(tǒng)，其潛在的碳匯功能亦不容忽視。基于渦度相關(guān)技術(shù)，于2011年12月1日至2012年11月30日對(duì)北京奧林匹克森林公園城市綠地生態(tài)系統(tǒng)進(jìn)行了碳通量觀測(cè)，以探討城市綠地生態(tài)系統(tǒng)碳交換及其與環(huán)境因子的關(guān)系及其源/匯屬性和強(qiáng)度。研究發(fā)現(xiàn)：奧林匹克森林公園綠地年總生態(tài)系統(tǒng)生產(chǎn)力(GEP)、生態(tài)系統(tǒng)呼吸(Re)、生態(tài)系統(tǒng)凈生產(chǎn)力(NEP)具有明顯的季節(jié)變化，生長(zhǎng)季(4月—11月)以吸收二氧化碳(CO2)為主，非生長(zhǎng)季以釋放CO2為主。Re隨空氣溫度(Ta)呈指數(shù)增加，溫度敏感性系數(shù)(Q10)為2.5；GEP也隨Ta的升高而增加；GEP與Re對(duì)Ta的響應(yīng)差異決定著NEP與Ta的關(guān)系：當(dāng)Talt; 10.0 ℃時(shí)，NEP隨Ta升高而下降；當(dāng)Tagt; 10.0 ℃時(shí)， NEP隨Ta升高而增加。在生長(zhǎng)季各月，日總GEP隨日光合有效輻射(PAR)的升高而增加，生態(tài)系統(tǒng)光合作用表觀光量子效率(α)和平均最大光合速率(Amax)也表現(xiàn)出明顯的季節(jié)變化，最大值出現(xiàn)在7月，分別為0.083 μmol CO2/μmol PAR 和29.46 μmol·m-2·s-1，最小值出現(xiàn)在11月，分別為0.017 μmol CO2/μmol PAR和4.16 μmol·m-2·s-1。奧林匹克森林公園綠地全年GEP 、Re、NEP的年總量分別為1192、1028、164 g C/m2。該研究結(jié)果可用于估算、模擬預(yù)測(cè)相似城市生態(tài)系統(tǒng)在氣候變化背景下生態(tài)系統(tǒng)凈碳交換，可作為城市綠地生態(tài)系統(tǒng)管理與應(yīng)對(duì)氣候變化的重要理論基礎(chǔ)。

城市綠地；碳通量；渦度協(xié)方差；生態(tài)系統(tǒng)生產(chǎn)力；生態(tài)系統(tǒng)呼吸

近20年來(lái)，我國(guó)森林碳(C)儲(chǔ)量增加了約0. 4 Pg C，年增加0. 011 — 0. 035 Pg C[1]。人工造林面積的不斷擴(kuò)大對(duì)森林碳儲(chǔ)量的增加起著重要貢獻(xiàn)。目前，我國(guó)人工林占全球人工林面積的近1/3；年增量占全球的53.2%[2]。隨著城市化進(jìn)程的推進(jìn)，城市森林公園的面積也在不斷增加，城市綠地/人工林作為城市生態(tài)的理想“調(diào)節(jié)器”，深刻影響著城市的生態(tài)環(huán)境[3]。在碳循環(huán)與氣候變化研究中，以人工林植被為主要存在形態(tài)的城市綠地作為重要的碳平衡組分，對(duì)減輕氣候變化的影響發(fā)揮著重要作用，其潛在的碳匯功能(或碳固持)不可低估[4]。

城市由于受特殊的下墊面、建筑物及人類活動(dòng)的影響，表現(xiàn)出高溫、高大氣二氧化碳(CO2)濃度、高氮沉降以及高氣溶膠濃度等環(huán)境特征[5]，城市中植被的性質(zhì)和特征也與自然植被有很大不同[6]。作為城市綠化美化、改善生態(tài)環(huán)境的重要載體的城市綠地生態(tài)系統(tǒng)，其豐富的林種因人工管理、熱島效應(yīng)和相對(duì)較高的氮沉降等可能使單株樹木的生長(zhǎng)速度高于野外[7]；同時(shí)，城市中較高的氣溶膠濃度，削弱了地面接收的太陽(yáng)輻射，但增加了散射輻射的比例，而大氣懸浮氣溶膠粒子直接吸收短波輻射使低層大氣中上部增溫[8]，使不同植物凈光合速率、蒸騰速率和氣孔導(dǎo)度發(fā)生不同程度的增加或減少，造成城市單位面積碳固持能力增加或降低[9]。近年的一些研究認(rèn)為，城市的特征性環(huán)境可作為研究生態(tài)系統(tǒng)對(duì)全球氣候變化的響應(yīng)和預(yù)測(cè)未來(lái)生態(tài)系統(tǒng)的變化趨勢(shì)的“自然大實(shí)驗(yàn)場(chǎng)”[10- 11]。然而，這一方法仍處于前期的討論階段，相關(guān)的研究也僅限于一些有關(guān)城市熱島和城市污染等碳源方面的報(bào)道，對(duì)城市綠地生態(tài)系統(tǒng)碳匯功能及與全球變化相關(guān)的研究仍相對(duì)薄弱。通過(guò)對(duì)城市綠地生態(tài)系統(tǒng)碳匯功能的量化研究，可以評(píng)估并預(yù)測(cè)其在減少大氣CO2中的潛在作用[12]。歐美國(guó)家對(duì)城市綠地生態(tài)系統(tǒng)碳固持功能已有大量研究[13- 15]，近年來(lái)許多發(fā)展中國(guó)家對(duì)城市綠地生態(tài)系統(tǒng)碳固持功能也越來(lái)越重視[16- 18]。對(duì)城市綠地生態(tài)系統(tǒng)碳交換規(guī)律及其環(huán)境控制因素的理解和認(rèn)識(shí)對(duì)我國(guó)陸地生態(tài)系統(tǒng)碳循環(huán)和氣候變化研究具有重要的理論和實(shí)踐意義[19]。

渦度協(xié)方差法已成為目前直接測(cè)定植被與大氣間的碳水交換，認(rèn)識(shí)和理解生態(tài)系統(tǒng)碳源/匯屬性動(dòng)態(tài)的最可靠方法之一[20]，廣泛應(yīng)用于估算從局域到景觀尺度的物質(zhì)與能量收支[21]。渦度相關(guān)系統(tǒng)通過(guò)對(duì)風(fēng)速、溫度、水汽(H2O)和CO2脈動(dòng)值的觀測(cè)，可實(shí)現(xiàn)對(duì)城市綠地生態(tài)系統(tǒng)的碳-水-能量通量的監(jiān)測(cè)，不但有助于分析城市綠地與大氣間CO2交換的規(guī)律及其環(huán)境控制因素，而且可以評(píng)價(jià)渦度協(xié)方差法在城市綠地生態(tài)系統(tǒng)研究中的適用性。然而已有的渦度通量研究主要集中在森林、農(nóng)田和草地等生態(tài)系統(tǒng)，利用該技術(shù)對(duì)城市綠地生態(tài)系統(tǒng)物質(zhì)與能量交換的連續(xù)觀測(cè)仍不多見[22]。通過(guò)探索城市綠地生態(tài)系統(tǒng)對(duì)氣候變化的響應(yīng)，利用碳通量數(shù)據(jù)，可預(yù)測(cè)氣候變化背景下城市綠地生態(tài)系統(tǒng)碳匯功能，為反演生態(tài)系統(tǒng)對(duì)氣候和環(huán)境變化的響應(yīng)關(guān)鍵碳循環(huán)參數(shù)提供科學(xué)依據(jù)。本研究運(yùn)用渦度協(xié)方差法對(duì)北京奧林匹克公園典型綠地生態(tài)系統(tǒng)進(jìn)行了全年連續(xù)監(jiān)測(cè)，分析了(1)總生態(tài)系統(tǒng)生產(chǎn)力(GEP)，生態(tài)系統(tǒng)呼吸(Re)和生態(tài)系統(tǒng)凈生產(chǎn)力(NEP)的季節(jié)動(dòng)態(tài)；(2)環(huán)境因素對(duì)碳交換各組分的控制作用；(3)全年碳源/匯屬性及強(qiáng)度。

1 研究地區(qū)與研究方法

1.1 研究地區(qū)概況

研究區(qū)位于北京奧林匹克森林公園北園(40°01′N，116°23′E，海拔51 m)，占地680 hm2，植被覆蓋率達(dá)90%以上，是亞洲最大的城市森林公園。該區(qū)域?yàn)榈湫偷呐瘻貛О霛駶?rùn)大陸性季風(fēng)氣候，年均溫14.0 ℃(北京市氣象局，2007年)，1月均溫 -5.8 ℃，極端最低-27.4 ℃，7月均溫25.6 ℃，極端最高42.0 ℃，年均降雨量約600 mm，全年無(wú)霜期180—200d。降水季節(jié)分配不均，全年降水的80%集中在6、7、8三個(gè)月，7、8月暴雨和強(qiáng)對(duì)流天氣較頻繁。土壤類型主要為潮褐土，土壤容重為1.62 g/cm3，土壤孔隙度約為40.30%。通量塔足跡內(nèi)主要植被為人工營(yíng)造的喬灌草復(fù)層景觀林，喬木主要為油松(Pinustabulaeformis)，側(cè)柏(Platycladusorientalis)，國(guó)槐(Sophorajaponica)，白蠟(Fraxinuschinensis)和銀杏(Ginkgobiloba)，林分密度約為210 株/hm2，平均樹高7.7 m，胸徑20.6 cm；灌木樹種主要為山桃(Prunusdavidiana)和榆葉梅(Amygdalustriloba)，叢生灌木主要為紅瑞木(SwidaalbaOpiz)，丁香(Syzygiumaromaticum)，林分密度約為351株/hm2，平均樹高2.8 m，平均基徑72.0 cm；地被植物主要為鳶尾(Lristectorum)和石竹(Dianthuschinensis)。

1.2 研究方法

渦度通量觀測(cè)儀器安裝高度為11.5 m，包括三維超聲風(fēng)速儀(CSAT3, Campbell Scientific Ltd, USA)和CO2/H2O紅外氣體分析儀(EC155, Campbell Scientific Ltd, USA)。凈輻射儀(CNR- 4, Kipp amp; Zonen Inc., Netherlands)及光量子傳感器(PARLITE, Kipp amp; Zonen Inc., Netherlands)均安裝在12 m處；兩個(gè)翻斗式雨量筒(TE525，Campbell Scientific Ltd, USA)分別安裝在距離通量塔基礎(chǔ)2 m和5 m處迎風(fēng)向的位置；5個(gè)空氣溫濕度傳感器(HMP45C，Campbell Scientific Ltd, USA)分別安裝在在1.5、5、8、10 m和12 m高度；土壤溫度傳感器(Campbell- 109, Campbell Scientific Ltd, USA)及土壤含水量傳感器(CS616, Campbell Scientific Ltd, USA)分5個(gè)方向布設(shè)在塔周邊10—15 m處，其中4個(gè)點(diǎn)安裝在10 cm深度，一個(gè)點(diǎn)安裝在30 cm深度；土壤熱通量板(HFP01, Campbell Scientific Ltd, USA)分4個(gè)方向安裝在10 cm深度。數(shù)據(jù)采集器(CR3000, Campbell Scientific Ltd, USA)以10Hz頻率記錄渦度協(xié)方差系統(tǒng)的數(shù)據(jù)，對(duì)于微氣象系統(tǒng)記錄30 min平均值(掃描頻率0.1 Hz)。

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

本研究所用數(shù)據(jù)從2011年12月1日至2012年11月30日，期間由于儀器故障、停電等造成16.7%渦度數(shù)據(jù)缺失。對(duì)10 Hz數(shù)據(jù)進(jìn)行了校正和質(zhì)量控制，包括剔除野點(diǎn)[23]、二次坐標(biāo)軸旋轉(zhuǎn)[24]、去除傳感器延時(shí)影響[25]、頻率響應(yīng)校正[26]和WPL校正[27]等。用校正過(guò)的10 Hz數(shù)據(jù)計(jì)算半小時(shí)通量，然后利用3倍標(biāo)準(zhǔn)差法剔除野點(diǎn)。生態(tài)系統(tǒng)凈碳交換(NEE)被定義為公式(1)[28]：

NEE=Fc+Fs+Vc

(1)

式中，F(xiàn)c是通量觀測(cè)塔在植被上部所觀測(cè)到的CO2通量，F(xiàn)s是渦度相關(guān)系統(tǒng)安裝高度下冠層內(nèi)儲(chǔ)存通量，Vc是垂直和水平平流效應(yīng)。在地勢(shì)平坦，植被類型空間分布均勻的植被下墊面，Vc可忽略不計(jì)。對(duì)于較高的植被冠層而言，觀測(cè)高度以下的Fs變化常常不為零，因此，需要對(duì)Fs進(jìn)行計(jì)算[29]：

(2)

式中,ΔC(Z) 為觀測(cè)高度處的CO2濃度變化(μmol/m3)，Δt為采樣周期(s)，ΔZ為觀測(cè)高度(m)。

對(duì)夜間NEE(NEEnight)數(shù)據(jù)進(jìn)行摩擦風(fēng)速(u*)篩選，剔除u*lt;0.15 m/s的數(shù)據(jù)。篩選后得到有效數(shù)據(jù)占總數(shù)據(jù)比例約為53%。通過(guò)通量貢獻(xiàn)分析(足跡分析)[30]，得出90%通量足跡為上風(fēng)向150 m范圍內(nèi)，在100 m的footprint貢獻(xiàn)為70%，而通量貢獻(xiàn)的峰值則出現(xiàn)在50 m處。

對(duì)由于儀器故障、停電和數(shù)據(jù)質(zhì)量控制等造成的缺失數(shù)據(jù)用如下方法進(jìn)行插補(bǔ)：(1) 小于2 h的缺失數(shù)據(jù)用線性內(nèi)插法插補(bǔ)；(2) 大于2 h 的缺失數(shù)據(jù)采用非線性回歸法進(jìn)行插補(bǔ)[31]。對(duì)于生長(zhǎng)季白天光合有效輻射(PAR)大于3 μmol·m-2·s-1[32]的渦度數(shù)據(jù)插補(bǔ)采用如下光響應(yīng)模型[33]：

(3)

式中，NEEday是白天的生態(tài)系統(tǒng)凈碳交換量(μmol CO2·m-2·s-1)，α為生態(tài)系統(tǒng)表觀光量子效率(μmol CO2/μmol光量子)，PAR為光合有效輻射(μmol·m-2·s-1)，Amax為生態(tài)系統(tǒng)最大光合速率(μmol CO2·m-2·s-1)，Rd為白天平均生態(tài)系統(tǒng)呼吸速率(μmol CO2·m-2·s-1)?？紤]到植物的物候和光合參數(shù)(α,Amax)的季節(jié)變化，采用公式(3)按月對(duì)數(shù)據(jù)進(jìn)行擬合和插補(bǔ)。

預(yù)分析發(fā)現(xiàn)：冬季和夏季溫度變化范圍較窄，NEEnight對(duì)空氣溫度(Ta)的回歸效果不佳(R2lt; 0.2，Pgt; 0.05)，而春季和秋季NEEnight與Ta的反應(yīng)相似。因此，在插補(bǔ)夜間數(shù)據(jù)時(shí)全年采用單一函數(shù)(公式4)對(duì)夜間日平均Ta與NEEnight進(jìn)行擬合[34]。

(4)

式中，R0為當(dāng)Ta為0 ℃時(shí)的呼吸速率；Q10為溫度增加10 ℃時(shí)呼吸速率的相對(duì)增加量(i.e., 溫度敏感性)。再將公式(4)獲得的呼吸參數(shù)(R0,Q10)外推到白天獲得白天生態(tài)系統(tǒng)呼吸的日變化。

因此凈生態(tài)系統(tǒng)生產(chǎn)力(NEP)為：

NEP=-NEE=GEP-Re

(5)

本文中的NEE和Re正值表示生態(tài)系統(tǒng)向大氣中放出CO2，NEP和GEP正值則表示生態(tài)系統(tǒng)從大氣中吸收CO2。通過(guò)回歸分析生態(tài)系統(tǒng)碳交換及其組分與環(huán)境變量之間的關(guān)系。數(shù)據(jù)統(tǒng)計(jì)分析與作圖均采用Matlab(Version 7.12.0., The MathWorks, Natick, MA, USA)完成。

圖1 奧林匹克森林公園日平均氣溫(Ta)、日累積光和有效輻射(PAR)、日降雨量(precipitation)和飽和水汽壓差 (VPD) 季節(jié)動(dòng)態(tài)Fig.1 Mean daily air temperature (Ta), integrated daily photosynthetically active radiation (PAR), precipitation and daily mean vapour pressure deficit (VPD)

2 結(jié)果與分析

2.1 環(huán)境因子的季節(jié)變化

生態(tài)系統(tǒng)碳交換通常受溫度、光合有效輻射、濕度、風(fēng)速和下墊面特征等因素的影響。在本研究數(shù)據(jù)時(shí)間段內(nèi)北京奧林匹克森林公園最高氣溫(Ta)發(fā)生在7月2日(29.6 ℃)，最低氣溫發(fā)生在1月22日(-7.9 ℃)，月平均氣溫從3月到11月均大于0 ℃ (圖1)。光合有效輻射(PAR)的季節(jié)趨勢(shì)為單峰型，最大日總量出現(xiàn)在6月，為569.19 mol·m-2·d-1(圖1)。年降雨總量716.00 mm，其中2012年7月21日發(fā)生61年一遇的特大暴雨，日降雨量達(dá)到176.85 mm (圖1)。飽和水汽壓差(VPD)變化范圍為0.15 — 2.81 Kpa (圖1)，呈明顯季節(jié)變化，總體表現(xiàn)為夏季高、冬季低。

2.2 生態(tài)系統(tǒng)的GEP、Re及NEP的年總量和季節(jié)動(dòng)態(tài)

圖2 2011年12月—2012年11月，日總生態(tài)系統(tǒng)生產(chǎn)(GEP)、日生態(tài)系統(tǒng)呼吸(Re)和日生態(tài)系統(tǒng)凈生產(chǎn)力(NEP)Fig.2 Daily gross ecosystem productivity (GEP)，daily ecosystem respiration (Re) and daily net ecosystem productivity (NEP) from December, 2011 to November, 2012

2012年觀測(cè)點(diǎn)年總生態(tài)系統(tǒng)生產(chǎn)(GEP)、生態(tài)系統(tǒng)呼吸(Re) 、生態(tài)系統(tǒng)凈生產(chǎn)力(NEP)的年總量分別為1192、1028、164 g C/m2。在城市綠地植被中，喬木的儲(chǔ)碳能力是整個(gè)城市綠地中儲(chǔ)碳的主要部分，且相同數(shù)量的老樹要比幼樹對(duì)碳儲(chǔ)存的貢獻(xiàn)高很多[35]。與喬木相比，草本植物固定下來(lái)的大部分碳并不能長(zhǎng)期積累；灌叢由于修剪受到嚴(yán)格的形體控制，其徑向增粗并不明顯，因而其固碳量隨時(shí)間的增長(zhǎng)也不大。在奧林匹克森林公園中，由于多數(shù)植被為人工建設(shè)，歷史遺留較少，因而大樹的比例很低。因此，奧林匹克森林公園這個(gè)年輕公園的植被年凈碳吸收量6.04 tCO2/hm2小于北京海淀公園綠地8.76 tCO2/hm2[3]。而Jo[36]對(duì)韓國(guó)Chuncheon, Kangleung和Seoul 3個(gè)城市的綠地碳通量進(jìn)行了調(diào)查結(jié)果表明，在不考慮草地與土壤的碳通量的情況下，木本植物年凈碳吸收量為5.87—14.33 tCO2/hm2。與其他林地碳通量觀測(cè)結(jié)果相比較，美國(guó)密執(zhí)安大學(xué)生物站的年NEP在2.93—6.23 tCO2/hm2之間[ 37]，美國(guó)哈佛林地的年NEP是7. 33 tCO2/hm2[38]。由于引起碳吸收能力不同的原因仍存在很多的不確定性，不同年際之間環(huán)境因素對(duì)碳交換各組分的控制作用亦存在差異，導(dǎo)致本研究結(jié)果與以上研究結(jié)果相比相當(dāng)，但顯著小于野外年輕的森林44.6 tCO2/hm2[39]。因此，必須進(jìn)一步發(fā)展針對(duì)城市植被的碳循環(huán)研究方法以更準(zhǔn)確地估算城市地區(qū)植被的碳吸收和碳固持[40]。

2.3 溫度對(duì)生態(tài)系統(tǒng)碳交換的影響

日平均NEEnight隨Ta的增加呈指數(shù)增長(zhǎng)(圖3)，Q10為2.5，R0為0.55 μmol CO2·m-2·s-1，Ta可解釋73%的NEEnight變異。本研究得到的Q10值與其它溫帶地區(qū)的研究結(jié)果接近，比如鄭澤梅等[41]對(duì)長(zhǎng)白山溫帶混交林的研究發(fā)現(xiàn)林冠下層夜間CO2通量隨溫度呈指數(shù)增長(zhǎng)，Q10為2.7，高于在模擬呼吸時(shí)通常作為默認(rèn)值的2.0[42]。

圖3 夜間生態(tài)系統(tǒng)凈交換(NEEnight)、日平均生態(tài)系統(tǒng)生產(chǎn)(GEP)、凈生態(tài)系統(tǒng)碳交換(NEE)分別與同期空氣溫度的關(guān)系Fig.3 Relationship between nighttime net ecosystem exchange (NEEnight), daily mean gross ecosystem productivity (GEP), net ecosystem exchange (NEE) and the corresponding period air temperature, respectively

溫度除了影響生態(tài)系統(tǒng)呼吸，還通過(guò)影響光合作用來(lái)影響生態(tài)系統(tǒng)的CO2交換(圖3)。隨著季節(jié)的變化，春季GEP隨著Ta的升高逐漸增加，隨后在夏季達(dá)到最高，隨著秋季Ta的降低，GEP也開始降低。許多研究發(fā)現(xiàn)，GEP隨溫度呈拋物線變化趨勢(shì)[43]，在溫度上升到一個(gè)最適溫度(GEP達(dá)到最大值時(shí)的溫度)后，GEP不再繼續(xù)升高，而是隨著Ta的升高而降低，說(shuō)明當(dāng)Ta升高到一定程度時(shí)(即超過(guò)臨界的最適溫度)，對(duì)GEP會(huì)產(chǎn)生抑制作用。但在本研究中，Ta對(duì)GEP的影響不存在最適溫度，即當(dāng)Ta達(dá)到最高溫度時(shí)，也未對(duì)GEP產(chǎn)生抑制作用，植物對(duì)Ta具有一定的適應(yīng)性。

奧林匹克森林公園綠地凈生態(tài)系統(tǒng)交換在氣溫小于10 ℃時(shí)，整個(gè)生態(tài)系統(tǒng)以呼吸作用為主，此時(shí)公園綠地處于非生長(zhǎng)季，Re隨溫度的升高而增強(qiáng), 導(dǎo)致日總NEE隨溫度增大(圖3)。氣溫大于10 ℃植被進(jìn)入生長(zhǎng)季，光合作用隨溫度增加比呼吸快，日總NEE隨著溫度的增加發(fā)生由正到負(fù)的變化，表現(xiàn)為凈吸收大氣CO2，30 ℃時(shí)表現(xiàn)為最大吸收(NEE =-4.01 g C·m-2·d-1)。該結(jié)果與Niu[43]對(duì)20個(gè)中高緯度通量站點(diǎn)的研究結(jié)果不盡相同，本研究結(jié)果僅與該研究中的兩個(gè)落葉闊葉林站點(diǎn)(US-UMB站及DE-Hai站)的NEE與Ta的關(guān)系一致，而其它18個(gè)站點(diǎn)的日總NEE均先隨Ta的增加而增加，當(dāng)Ta達(dá)到最適溫度(即日總NEE最大吸收的峰值時(shí))后，隨著Ta的繼續(xù)增加，日總NEE開始減少與降低，日總NEE隨Ta呈拋物線變化趨勢(shì)。

表1白天生態(tài)系統(tǒng)凈碳交換(NEE)與光合有效輻射(PAR)之間的相關(guān)參數(shù)比較

Table1ComparisonofparametersfortherelationshipbetweendaytimenetecosystemCO2exchange(NEE)andphotosyntheticallyactiveradiation(PAR)

月份Month/月光合作用表觀光量子效率Ecosystemquantumyield(α)/(μmolCO2/μmolPAR)最大光合速率Maximumphotosynthesis(Amax)/(μmol·m-2·s-1)呼吸速率RespirationrateRd/(μmol·m-2·s-1)40．0129．912．7750．04816．764．1860．06918．335．1470．08329．468．1680．07527．856．8690．03524．213．04100．01924．091．90110．0114．161．63

2.4 光合有效輻射與白天凈生態(tài)系統(tǒng)交換的關(guān)系

在生長(zhǎng)季，白天凈生態(tài)系統(tǒng)碳交換主要受PAR調(diào)控。生長(zhǎng)季各月PAR與對(duì)應(yīng)的白天凈生態(tài)系統(tǒng)碳交換(NEEday)關(guān)系顯著(圖4，表1)。光合作用隨著PAR的增大而逐漸增強(qiáng)，NEE負(fù)值也隨之增大，隨著生態(tài)系統(tǒng)對(duì)CO2的吸收量增大的趨勢(shì), 系統(tǒng)的碳匯能力也逐漸增大，白天最大光合速率(Amax)和平均白天生態(tài)系統(tǒng)呼吸速率(Rd)均在7月最大，與光和有效輻射表現(xiàn)出相同的季節(jié)變化趨勢(shì)(表1)。6—8月的生態(tài)系統(tǒng)光合作用表觀光量子效率α為0.068 — 0.083，平均值為0.075，該結(jié)果比Ruimy[44]報(bào)道的一般針葉林的0.024以及闊葉森林的0.037高，比Havard中緯度森林的0.04[45]、丹麥山毛櫸森林的0.06[46]以及華北平原冬小麥農(nóng)田生態(tài)系統(tǒng)的0.068也高[47]。由于城市中較高的氣溶膠濃度削弱了地面接收的太陽(yáng)輻射，減弱了植物的光合作用，但其散射輻射的比例很大，本研究結(jié)果與很多研究者認(rèn)為的光合利用效率的強(qiáng)度在散射輻射條件中是直射輻射中的兩倍，碳收支對(duì)散射輻射變化的敏感度高于對(duì)直射輻射變化的敏感度的結(jié)論相吻合[48]。

圖4 2012年生長(zhǎng)季各月運(yùn)用直角雙曲線模型(公式3)擬合的白天生態(tài)系統(tǒng)碳交換(NEEday)和光合有效輻射(PAR)的關(guān)系Fig.4 Relationship between daytime NEE and photosynthetically active radiation (PAR) using rectangular hyperbolic model(Formula 3)for each month during the growing season in 2012

3 結(jié)論

生態(tài)系統(tǒng)的碳循環(huán)已經(jīng)成為全球研究的熱點(diǎn), 而城市地區(qū)碳循環(huán)數(shù)據(jù)的缺乏則成為研究區(qū)域乃至全球碳循環(huán)的瓶頸[49]。對(duì)城市公園綠地CO2通量的長(zhǎng)期連續(xù)觀測(cè)，有助于分析城市綠地生態(tài)系統(tǒng)光合生產(chǎn)和呼吸的動(dòng)態(tài)變化，了解綠地的碳收支狀況及其固碳能力。本文通過(guò)對(duì)奧林匹克森林公園綠地2011年12月2012年11月的通量觀測(cè)數(shù)據(jù)進(jìn)行處理與缺失值插補(bǔ)，并通過(guò)與同期空氣溫度、太陽(yáng)輻射等氣象數(shù)據(jù)的相關(guān)性分析，定量研究了奧林匹克森林公園綠地CO2通量的季節(jié)變化以及影響因子，評(píng)價(jià)了綠地的固碳能力。

研究發(fā)現(xiàn)：(1)奧林匹克森林公園綠地年總生態(tài)系統(tǒng)生產(chǎn)(GEP)、生態(tài)系統(tǒng)呼吸(Re) 、年生態(tài)系統(tǒng)凈生產(chǎn)力(NEP)在2012年內(nèi)具有明顯的季節(jié)變化，植物生長(zhǎng)季4月—11月份以吸收CO2為主，非生長(zhǎng)季以釋放CO2為主。(2)Re隨溫度(Ta)呈指數(shù)增加，溫度敏感系數(shù)(Q10)為2.5; GEP也隨Ta的升高而增加；GEP與Re對(duì)溫度的響應(yīng)差異決定著NEP與Ta的關(guān)系：當(dāng)Talt; 10 ℃時(shí)，NEP隨Ta升高而下降；當(dāng)Tagt; 10 ℃時(shí)，NEP隨Ta升高而增加。(3)在生長(zhǎng)季各月，日總GEP隨日光合有效輻射(PAR)的升高而增加，生態(tài)系統(tǒng)光合作用表觀光量子效率(α)和平均最大光合速率(Amax)也表現(xiàn)出明顯的季節(jié)變化，最大值出現(xiàn)在7月，分別為0.083 μmol CO2/μmol PAR 和29.46 μmol·m-2·s-1，最小值出現(xiàn)在11月，分別為0.017 μmol CO2/μmol PAR和4.16 μmol·m-2·s-1。(4)奧林匹克森林公園綠地全年GEP 、Re、NEP的年總量分別為1192、1028、164 g C/m2。由于城市綠地生態(tài)系統(tǒng)在群落組成結(jié)構(gòu)、氣候條件和干擾歷史等方面的差異，本研究結(jié)果在區(qū)域尺度的空間代表性有限。未來(lái)研究需要結(jié)合多個(gè)監(jiān)測(cè)站點(diǎn)的數(shù)據(jù)，通過(guò)數(shù)據(jù)整合與模型參數(shù)化實(shí)現(xiàn)從生態(tài)系統(tǒng)水平到區(qū)域的擴(kuò)展，從而實(shí)現(xiàn)未來(lái)氣候變化情境下精確預(yù)測(cè)城市綠地對(duì)區(qū)域碳收支的貢獻(xiàn)。

致謝:唐祥和王奔對(duì)渦度通量監(jiān)測(cè)設(shè)備維護(hù)與運(yùn)行、數(shù)據(jù)整理與分析提供幫助。

[1] Fang J Y, Chen A P, Peng C H, Zhao S, Ci L. Changes in forest biomass carbon storage in China between 1949 and 1998. Science, 2001, 291(5525): 2320- 2322.

[2] Tong X J, Zhang J S, Meng P, Huang H, Guo L, Yin C J, Cao J. Characteristics of net carbon exchange over a mixed plantation in a hilly area of the north China. Scientia Silvae Sinicae, 2010, 46(3): 37- 43.

[3] Li X, Sun R, Li Y, Wang X X, Xie G H, Yan X D, Zhu Q J. Carbon dioxide fluxes on green space in Haidian Park, Beijing. Acta Ecologica Sinica, 2010, 30(24): 6715- 6725.

[4] Millward A A, Sabir S. Benefits of a forested urban park: What is the value of Allan Gardens to the city of Toronto, Canada? Landscape and Urban Planning, 2011, 100(3): 177- 188.

[5] George K, Ziska L H, Bunce J A, Quebedeaux B. Elevated atmospheric CO2concentration and temperature across an urban-rural transect. Atmospheric Environment, 2007, 41(35): 7654- 7665.

[6] Wen J S, Ge Y, Jiao L, Deng Z P, Peng C H, Chang J. Does urban land use decrease carbon sequestration? -A case study in Taizhou, China. Chinese Journal of Plant Ecology, 2010, 34(6): 651- 660.

[7] Green D M, Olesyszyn M. Enzyme activities and carbon dioxide flux in a Sonoran Desert urban ecosystem. Soil Science Society of America Journal, 2002, 66(6): 2002- 2008.

[8] Wang H X, Huang J G, Chen C H. The influences of urban aerosols on solar radiation and their reflection in urban boundary layer temperature variation. Acta Meteorologica Sinica, 1993, 51(4): 457- 464.

[9] Shi F S, Wu N, Wu Y, Wang Q. Effect of simulated temperature enhancement on growth and photosynthesis of Deschampsia caespitosa and Thlaspi arvense in Northwestern Sichuan, China. Chinese Journal of Applied amp; Environmental Biology, 2009, 15(6): 750- 755.

[10] Shen W J, Wu J G, Grimm N B, Hope D. Effects of urbanization-induced environmental changes on ecosystem functioning in the Phoenix metropolitan region, USA. Ecosystems, 2008, 11(1): 138- 155.

[11] Zhao D H, Ouyang Y, Qi J G, Fan Y M, An S Q. Near-surface atmospheric properties along an urban center-suburban-urban forest gradient in summer in Nanjing City. Acta Ecologica Sinica, 2009, 29(12): 6654- 6663.

[12] Nowak D J, Crane D E. Carbon storage and sequestration by urban trees in the USA. Environmental Pollution, 2002, 116(3): 381- 389.

[13] McHale M R, Mcpherson E G, Burke I C. The potential of urban tree plantings to be cost effective in carbon credit markets. Urban Forestry and Urban Greening, 2007, 6(1): 49- 60.

[14] McHale M R, Burke I C, Lefsky M A, Peper P J, McPherson E G. Urban forest biomass estimates: is it important to use allometric relationships developed specifically for urban trees? Urban Ecosystems, 2009, 12(1): 95- 113.

[15] Escobedo F, Varela F, Zhao M, Wagner J E, Zipperer W. Analyzing the efficacy of subtropical urban forests in offsetting carbon emissions from cities. Environmental Science and Policy, 2010, 13(15): 362- 372.

[16] Yang J, Mcbride J, Zhou J J, Shen Z Y. The urban forest in Beijing and its role in air pollution reduction. Urban Forestry and Urban Greening, 2005, 3(2): 65- 78.

[17] Stoffberg G H, van Rooyen M W, van der Linde M J, Groeneveld H T. Carbon sequestration estimates of indigenous street trees in the City of Tshwane South Africa. Urban Forestry and Urban Greening, 2010, 9(1): 9- 14.

[18] Zhao M, Kong Z H, Escobedo F J, Gao J. Impacts of urban forests on offsetting carbon emissions from industrial energy use in Hangzhou, China. Journal of Environmental Management, 2010, 91(4): 807- 813.

[19] Liu C F, Li X. Carbon storage and sequestration by urban forests in Shenyang, China. Urban Forestry and Urban Greening, 2011, 11(2): 121- 128.

[20] Baldocchi D D. Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: past, present and future. Global Change Biology, 2003, 9(4): 479- 492.

[21] Papale D, Reichstein M, Aubinet M, Canfora E, Bernhofer C, Kutsch W, Longdoz B, Rambal S, Valentini R, Vesala T, Yakir D. Towards a standardized processing of Net Ecosystem Exchange measured with eddy covariance technique: algorithms and uncertainty estimation. Biogeosciences, 2006, 3(4): 571- 583.

[22] Wang X X, Zhu Q J, Chen S H, Liu X, Hu Y M. Remote sensing retrieval of water and heat fluxes over urban green space and experimental validation. Areal Research and Development, 2010, 29(3): 63- 66.

[23] Vickers D, Mahrt L. Quality control and flux sampling problems for tower and aircraft data. Journal of Atmospheric and Oceanic Technology, 1997, 14(3): 512- 526.

[24] Nakai M, Kashiwazaki N, Takizawa A, Maedomari N, Ozawa M, Noguchi J, Kaneko H, Shino M, Kikuchi K. Morphologic changes in boar sperm nuclei with reduced disulfide bonds in electro stimulated porcine oocytes. Reproduction, 2006, 131(3): 603- 611.

[25] Fan S M, Wofsy S C, Bakwin P S, Jacob D J, Fitzjarrald D R. Atmosphere-biosphere exchange of CO2and O3in the central Amazon forest. Journal of Geophysical Research, 1990, 95(16): 851- 864.

[26] Moncrieff J B, Massheder J M, De Bruin H, Elbersc J, Friborg T, Heusinkveld B, Kabat P, Scott S, Soegaard H, Verhoef A. A system to measure surface fluxes of momentum, sensible heat, water vapour and carbon dioxide. Journal of Hydrology, 1997, 188- 189: 589- 611.

[27] Webb E K, Pearman G I, Leuning R. Correction of flux measurements for density effects due to heat and water vapour transfer. Quarterly Journal of the Royal Meteorological Society, 1980, 106(447): 85- 100.

[28] Aubinet M, Grelle A, Ibrom A, Rannik ü, Moncrieff J, Foken T, Kowalski A S, Martin P H, Berbigier P, Bernhofer C, Clement R, Elbers J, Granier A, Grünwald T, Morgenstern K, Pilegaard K, Rebmann C, Snijders W, Valentini R, Vesala T. Estimates of the annual net carbon and water exchange of European forests: the EUROFLUX methodology. Advances in Ecological Research, 2000, 30: 113- 175.

[29] Li C, He H L, Liu M, Su W, Fu Y L, Zhang L M, Wen X F, Yu G R. The design and application of CO2flux data processing system at ChinaFLUX. Gen-Information Science, 2008, 18(5): 557- 565.

[30] Schmid H P. Footprint modeling for vegetation atmosphere exchange studies: a review and perspective. Agricultural and Forest Meteorology, 2002, 113(1/4): 159- 183.

[31] Fange E, Baldocchi D, Olson R, Anthoni P, Aubinet M, Bernhofer C, Burba G, Ceulemans R, Clement R, Dolman H, Granier A, Gross P, Grünwald T, Hollinger D, Jensen N O, Katul G, Keronen P, Kowalski A, Lai C T, Law B E, Meyers T, Moncrieff J, Moors E, Munger J W, Pilegaard K, Rannik ü, Rebmann C, Suyker A, Tenhunen J, Tu K, Verma S, Vesala T, Wilson K, Wofsy S. Gap filling strategies for defensible annual sums of net ecosystem exchange. Agricultural and Forest Meteorology, 2001, 107(1): 43- 69.

[32] Zha T S, Xing Z S, Wang K Y, Kellom?ki S, Barr A G. Total and component carbon fluxes of a scots pine ecosystem from chamber measurements and eddy covariance. Annals of Botany, 2007, 99(2): 345- 353.

[33] Moren A S, Lindroth A. CO2exchange at the floor of a boreal forest. Agricultural and Forest Meteorology, 2000, 101(1): 1- 14.

[34] Zha T S, Kellom?ki S, Wang K Y, Rouvinen I. Carbon sequestration and ecosystem respiration for 4 years in a Scots pine forest. Global Change Biology, 2004, 10(9): 1492- 1503.

[35] Giardina C P, Harmon M E, Heath L S, Houghton R A, Jackson R B, Morrison J F, Murray B C, Pataki D E, Skog K E. A synthesis of current knowledge on forests and carbon storage in the United States. Ecological Applications, 2011, 21(6): 1902- 1924.

[36] Jo H K. Impacts of urban greenspace on offsetting carbon emissions for middle Korea. Journal of Environmental Management, 2002, 64(2): 115- 126.

[37] Schmid H P, Su H B, Vogel C S, Curtis P S. Ecosystem atmosphere exchange of carbon dioxide over a mixed hardwood forest in northern lower Michigan. Journal of Geophysical Research, 2003, 108(D14): 4417.

[38] Barford C C, Wofsy S C, Goulden M L, Munger J W, Pyle E H, Urbanski S P, Hutyra L, Saleska S R, Fitzjarrald D, Moore K. Factors controlling long-and short-term sequestration of atmospheric CO2in a mid-latitude forest. Science, 2001, 294(5547): 1688- 1691.

[39] Zhang J, Ge Y, Chang J, Jiang B, Jiang H, Peng C H, Zhu J R, Yuan W G, Qi L Z, Yu S Q. Carbon storage by ecological service forests in Zhejiang Province, subtropical China. Forest Ecology and Management, 2007, 245(1/3): 64- 75.

[40] Jo H K, McPherson E G. Carbon Storage and Flux in Urban Residential Greenspace. Journal of Environmental Management, 1995, 45(2): 109- 133.

[41] Zheng Z M, Zhang M, Wen X F, Sun X M, Yu G R, Zhang L M, Han S J, Wu J B. Contributions of below-canopy CO2fluxes to ecosystem carbon budget of a temperate mixed forest in Changbai Mountain. Acta Ecologica Sinica, 2009, 29(1): 1- 8.

[42] Goulden M L, Daube B C, Fan S M, Sutton D J, Bazzaz A, Munger J W, Wofsy S C. Gross CO2uptake by a black spruce forest. Journal of Geophysical Research, 1997, 102: 28987- 28996.

[43] Niu S L, Luo Q, Fei S F, Montagnani L, Bohrer G, Janssens I A, Gielen B, Rambal S, Moors E, Matteucci G. Seasonal hysteresis of net ecosystem exchange in response to temperature change: patterns and causes. Global Change Biology, 2011, 17(10): 3102- 3114.

[44] Ruim A, Jarvis P G, Baldocchi D D, Saugier B. CO2fluxes over plant canopies and solar radiation, A review. Advances in Ecological Research, 1996, 26: 1- 68.

[45] Wofsy S C, Goulden M L, Munger J W, Fan S M, Bakwin P S, Daube B C, Bassow S L, Bazzaz F A. Net exchange of CO2in a mid-latitude forest. Science, 1993, 260(5112): 1314- 1317.

[46] Pilegaard K, Hummelsho J, Jensen N O, Chen Z. Two years of continuous CO2eddy-flux measurements over a Danish beech forest. Agricultural and Forest Meteorology, 2001, 107(1): 29- 41.

[47] Wang J L, Wen X F, Sun X M, Wang Y Y. Asymmetry characteristic on the diurnal changes of CO2and H2O fluxes at full heading time of winter-wheat in north China. Acta Agriculturae Boreali-Sinica, 2009, 24(5): 159- 163.

[48] Gu L H, Baldocchi D, Verma S B, Black T A, Vesala T, Falge E M, Dowty P R. Advantages of diffuse radiation for terrestrial ecosystem productivity. Journal of Geophysical Research, 2002, 107(D6): ACL 2- 1-ACL 2- 23.

[49] Piao S L, Fang J Y, Ciais P, Peylin P, Huang Y, Sitch S, Wang T. The carbon balance of terrestrial ecosystems in China. Nature, 2009, 458(7241): 1009- 1014.

參考文獻(xiàn):

[2] 同小娟, 張勁松, 孟平, 黃輝, 國(guó)琳, 尹昌君, 高峻. 華北低丘山地人工混交林凈生態(tài)系統(tǒng)碳交換的變化特征. 林業(yè)科學(xué), 2010, 46(3): 37- 43.

[3] 李霞, 孫睿, 李遠(yuǎn), 王修信, 謝東輝, 嚴(yán)曉丹, 朱啟疆. 北京海淀公園綠地二氧化碳通量. 生態(tài)學(xué)報(bào), 2010, 30(24): 6715- 6725.

[6] 溫家石, 葛瀅, 焦荔, 鄧志平, 彭長(zhǎng)輝, 常杰. 城市土地利用是否會(huì)降低區(qū)域碳吸收能力? —臺(tái)州市案例研究. 植物生態(tài)學(xué)報(bào), 2010, 34(6): 651- 660.

[8] 王海嘯, 黃建國(guó), 陳長(zhǎng)和. 城市氣溶膠對(duì)太陽(yáng)輻射的影響及其在邊界層溫度變化中的反映. 氣象學(xué)報(bào), 1993, 51(4): 457- 464.

[9] 石福孫, 吳寧, 吳彥, 王乾. 模擬增溫對(duì)川西北高寒草甸兩種典型植物生長(zhǎng)和光合特征的影響. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2009, 15(6): 750- 755.

[11] 趙德華, 歐陽(yáng)琰, 齊家國(guó), 范亞民, 安樹青. 夏季南京市中心- 郊區(qū)- 城市森林梯度上的近地層大氣特征. 生態(tài)學(xué)報(bào), 2009, 29(12): 6654- 6663.

[22] 王修信, 朱啟疆, 陳聲海, 劉馨, 胡玉梅. 城市綠地水熱通量的遙感反演與實(shí)驗(yàn)驗(yàn)證. 地域研究與開發(fā), 2010, 29(3): 63- 66.

[29] 李春, 何洪林, 劉敏, 蘇文, 伏玉玲, 張雷明, 溫學(xué)發(fā), 于貴瑞. ChinaFLUX CO2通量數(shù)據(jù)處理系統(tǒng)與應(yīng)用. 地球信息科學(xué), 2008, 10(5): 557- 565.

[41] 鄭澤梅, 張彌, 溫學(xué)發(fā), 孫曉敏, 于貴瑞, 張雷明, 韓士杰, 吳家兵. 長(zhǎng)白山溫帶混交林林冠下層CO2通量對(duì)生態(tài)系統(tǒng)碳收支的貢獻(xiàn). 生態(tài)學(xué)報(bào), 2009, 29(1): 1- 8.

[47] 王建林, 溫學(xué)發(fā), 孫曉敏, 王玉英. 華北平原冬小麥生態(tài)系統(tǒng)齊穗期水碳通量日變化的非對(duì)稱響應(yīng). 華北農(nóng)學(xué)報(bào), 2009, 24(5): 159- 163.

DynamicsofCO2exchangeanditsenvironmentalcontrolsinanurbangreen-landecosysteminBeijingOlympicForestPark

CHEN Wenjing1, LI Chunyi1, HE Guimei2, 3, WANG Xiaoping2, 3, ZHA Tianshan1,3, JIA Xin1,*

1SchoolofSoilandWaterConservation,BeijingForestryUniversity,Beijing100083,China2BeijingForestryCarbonAdministration,Beijing100013,China3ChineseInstituteofGreenCarbon,Beijing100083,China

The area of urban forests and green-land is expanding dramatically across China in order to face rapid urbanization. Urban green-land ecosystems with plantations as their main vegetation type, have the great potential to sequestrate atmospheric carbon. Continuous measurements of CO2flux were made using eddy covariance technique from December 2011 to November 2012 in a mixed forest in Beijing Olympic Forest Park to quantify the seasonal dynamics of net ecosystem CO2exchange (NEE) and its responses to environmental factors. Gross ecosystem productivity (GEP), ecosystem respiration (Re), and net ecosystem productivity (NEP = -NEE) showed strong seasonal pattern, with CO2uptake dominating during the growing season from April to November, and a respiratory release of CO2dominating during the non-growing season. The carbon flux was influenced by photosynthetically active radiation (PAR), water vapor pressure deficit (VPD) and air temperature (Ta). In growing season, daytime net ecosystem carbon exchange (NEEday) increased with increasing PAR. The ecosystem quantum yield (α) and maximum photosynthesis (Amax) showed an apparent seasonal pattern, both peaking in July. VPD also affected NEE through its direct effect on photosynthesis. NEE increased with the increasing PAR up to a threshold of 1200 μmol·m-2·s-1, then decreased with increasing PAR above this threshold. GEP,Reand NEP were all influenced byTa, but responded differently.Reincreased exponentially with air temperature (Ta), with the temperature sensitivity (Q10) being 2.5. GEP also increased withTa. This differential response of GEP andRedetermined the relationship between NEP andTa. NEP decreased with increasingTawhenTalt; 10.0 ℃, but increased whenTagt; 10.0 ℃. NEEdayincreased with PAR. The ecosystem quantum yield (α) and maximum photosynthesis (Amax) showed an apparent seasonal pattern, both peaking in July with the value of 0.083 μmol CO2/μmol PAR and 29.46 μmol·m-2·s-1, respectively, and reaching a minimum in November with the value of 0.017 μmol CO2/μmol PAR and 4.16 μmol·m-2·s-1. The predicted annual totals of GEP,Reand NEP were 1192, 1028 and 164 g C/m2, respectively. The present results could contribute to the carbon budget of urban ecosystems, and help make carbon-oriented management strategies for sustainable urban development under global climate change.

urban green-land; CO2flux; eddy covariance; ecosystem productivity; ecosystem respiration

國(guó)家“十一五”科技支撐計(jì)劃項(xiàng)目(2008BAD95B07)

2013- 01- 07;

2013- 08- 07

*通訊作者Corresponding author.E-mail: xinjia@bjfu.edu.cn

10.5846/stxb201301070053

陳文婧，李春義，何桂梅，王小平，查天山，賈昕.北京奧林匹克森林公園綠地碳交換動(dòng)態(tài)及其環(huán)境控制因子.生態(tài)學(xué)報(bào),2013,33(20):6712- 6720.

Chen W J, Li C Y, He G M, Wang X P, Zha T S, Jia X.Dynamics of CO2exchange and its environmental controls in an urban green-land ecosystem in Beijing Olympic Forest Park.Acta Ecologica Sinica,2013,33(20):6712- 6720.