The offaxis integrated cavity output spectroscopy technique was used to measure air CO2 concentration, stable carbon (δ13C) and oxygen (δ18O) isotope ratios on the Fourth Ring Road (FRR) and in the green space system of Beijing Institute of Landscape Architecture (BILA) in summer and winter seasons. The variations of CO2 concentration, δ13C value, δ18O value and the differences of them between the FRR and the BILA, which were correlated with traffic volume and meteorological factors, were analyzed at halfhour timescale. The results showed that traffic volume on the FRR was large both in summer and winter with obvious morning and evening rush hours, and more than 150 thousands vehicles were observed everyday during the observation periods. Diurnal variation of the CO2 concentration showed a twopeak curve both on the FRR and in the green space system of the BILA. In contrast, diurnal variation of δ13C value was a twotrough curve while diurnal variation of δ18O value was a singletrough curve. The differences of CO2 concentration, δ13C value and δ18O value between the FRR and the green space system of BILA in summer were greater than those in winter. The carbon isotope partitioning results showed that in summer vehicle exhaust contributed 64.9% to total atmospheric CO2 of the FRR during measurement time, while heterotrophic respiration contributed 56.3% to total atmospheric CO2 of the green space system in BILA. However, in winter atmospheric CO2 from both the FRR and green space system mostly came from vehicle exhaust. Stepwise regression analysis indicated that differences of CO2 concentration between the FRR and green space system were significantly related to vehicle volume and solar radiation at halfhour timescale, while solar radiation and relative humidity were the main meteorological factors causing δ13C and δ18O differences between the FRR and green space system. Plants in the green space system strongly assimilated CO2 from fossil fuel burning by photosynthesis to maintain carbon and oxygen balance of urban area in the growing season, which played an important role in improving urban ecological environment.
全 文 :城市绿地系统内外 CO2浓度、δ
13C、δ18O
差异和来源及影响因素∗
孙守家1,2∗∗ 孟 平1,2 张劲松1,2 舒健骅3 郑 宁1,2
( 1中国林业科学研究院林业研究所 /国家林业局林木培育重点实验室, 北京 100091;2南京林业大学南方现代林业协同创新中
心, 南京 210037; 3北京市园林科学研究院 /绿化植物育种北京市重点实验室, 北京 100102 )
摘 要 采用离轴积分腔输出光谱技术测定夏季和冬季北京市 4 环路和北京园林科学研究
院绿地系统空气中 CO2浓度、δ13C和 δ18O值,在半小时尺度上分析了其变化特征和差异以及
与车流量和气象因子的关联.结果表明: 4环路上车流量较大,夏季和冬季观测期间每日均超
过 15万辆次,有明显的早晚交通高峰.4 环路与绿地系统的空气中 CO2浓度呈双峰曲线日变
化,δ13C值呈双波谷曲线、δ18O值呈单波谷曲线日变化,夏季二者空气中的 CO2浓度、δ13C 和
δ18O差值远大于冬季.同位素定量区分结果显示,夏季观测期间 4环路空气中 CO2主要来源于
机动车尾气,约占 64.9%,而绿地系统空气中 CO2主要来源于自养和异养呼吸,约占 56.3%,冬
季观测期间二者空气中 CO2均主要来源于机动车尾气.逐步回归分析表明,在半小时尺度上车
流量和太阳辐射对绿地系统与 4 环路 CO2浓度差值产生显著影响,太阳辐射和相对湿度则是
影响 δ13C和 δ18O差值的主要气象因子.绿地系统中的植物在生长季节通过光合作用吸收和消
减化石燃料燃烧产生的 CO2来维持城市碳氧平衡,在改善城市生态环境方面发挥重要作用.
关键词 城市绿地系统; 稳定同位素; CO2浓度; 车流量; 来源
∗北京市自然科学基金项目(8132044)和南京林业大学南方现代林业协同创新中心项目资助.
∗∗通讯作者. E⁃mail: ssj1011@ 163.com
2014⁃12⁃30收稿,2015⁃06⁃10接受.
文章编号 1001-9332(2015)10-3000-11 中图分类号 Q149 文献标识码 A
Differences and sources of CO2 concentration, carbon and oxygen stable isotope composition
between inside and outside of a green space system and influencing factors in an urban area.
SUN Shou⁃jia1,2, MENG Ping1,2, ZHANG Jin⁃song1,2, SHU Jian⁃hua3, ZHENG Ning1,2 ( 1Key La⁃
boratory of Tree Breeding and Cultivation of State Forestry Administration, Research Institute of
Forestry, Chinese Academy of Forestry, Beijing 100091, China; 2Collaborative Innovation Center of
Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China; 3 Bei⁃
jing Institute of Landscape Architecture / Beijing Key Laboratory of Greening Plants Breeding, Beijing
100102, China) . ⁃Chin. J. Appl. Ecol., 2015, 26(10): 3000-3010.
Abstract: The off⁃axis integrated cavity output spectroscopy technique was used to measure air CO2
concentration, stable carbon ( δ13C) and oxygen ( δ18O) isotope ratios on the Fourth Ring Road
(FRR) and in the green space system of Beijing Institute of Landscape Architecture (BILA) in
summer and winter seasons. The variations of CO2 concentration, δ13C value, δ18O value and the
differences of them between the FRR and the BILA, which were correlated with traffic volume and
meteorological factors, were analyzed at half⁃hour timescale. The results showed that traffic volume
on the FRR was large both in summer and winter with obvious morning and evening rush hours, and
more than 150 thousands vehicles were observed everyday during the observation periods. Diurnal
variation of the CO2 concentration showed a two⁃peak curve both on the FRR and in the green space
system of the BILA. In contrast, diurnal variation of δ13C value was a two⁃trough curve while diurnal
variation of δ18O value was a single⁃trough curve. The differences of CO2 concentration, δ13C value
and δ18O value between the FRR and the green space system of BILA in summer were greater than
those in winter. The carbon isotope partitioning results showed that in summer vehicle exhaust con⁃
应 用 生 态 学 报 2015年 10月 第 26卷 第 10期
Chinese Journal of Applied Ecology, Oct. 2015, 26(10): 3000-3010
tributed 64.9% to total atmospheric CO2 of the FRR during measurement time, while heterotrophic
respiration contributed 56.3% to total atmospheric CO2 of the green space system in BILA. Howe⁃
ver, in winter atmospheric CO2 from both the FRR and green space system mostly came from vehicle
exhaust. Stepwise regression analysis indicated that differences of CO2 concentration between the
FRR and green space system were significantly related to vehicle volume and solar radiation at half⁃
hour timescale, while solar radiation and relative humidity were the main meteorological factors cau⁃
sing δ13C and δ18O differences between the FRR and green space system. Plants in the green space
system strongly assimilated CO2 from fossil fuel burning by photosynthesis to maintain carbon and
oxygen balance of urban area in the growing season, which played an important role in improving
urban ecological environment.
Key words: urban green space system; stable isotope; CO2 concentration; traffic volume; source.
城市生态系统是由自然、经济和社会所组成的
复合系统,不仅土地利用方式和土地覆被变化强烈,
还是人口和交通的集中区,能源消费和化石燃料燃
烧非常集中.在燃烧过程中,化石燃料的碳大部分以
CO2形式排放,据统计新增 CO2排放量的 80%以上
来自于城区[1],成为大气主要碳源,导致城市生态
系统碳循环不同于自然生态系统.以往对碳循环研
究较少涉及城市碳循环过程,城市在全球碳循环中
的作用并未受到充分重视.近些年来我国城镇化发
展迅速,城市能源消耗持续增加,其碳循环过程对于
区域乃至全球气候变化的影响日益增强[2] .城市绿
地植物具有维持城市碳氧平衡的生态功能,对减缓
城市环境压力具有极其重要作用.因此,对绿地生态
系统影响城市空气中 CO2浓度和碳氧同位素组成进
行研究,对于了解城市绿地生态功能具有重要意义.
CO2中含有碳和氧 2 种元素,其同位素在物理
和化学性质上的差异,导致反应底物和生成产物在
同位素组成上出现差异[3] .因此,碳氧稳定同位素被
认为是理解生态系统碳氧交换的有效工具[4],被广
泛地应用到碳氧分配和转化[5]、生态系统碳转
移[6-7]和气体交换[8]研究中.以往对碳氧循环的研
究多集中于森林[9]、农田[10]、草地[11]、土壤[12]等单
一生态系统,对城市碳氧循环过程研究较少[13-15],
其原因是因为城市碳氧来源较多、不易于区分,且之
前测定碳氧同位素多用质谱仪,只能短暂或离散取
样,时间分辨率较低.近年来,离轴积分腔输出光谱
技术 ( off⁃axis integrated cavity output spectroscopy,
OA⁃ICOS)可以直接和连续观测空气中碳氧同位素
变化,响应时间较短,在野外环境中能用较高时间分
辨率观测和分析空气样品.
作为特大城市的北京,2013 年末全市机动车拥
有量超过 543.7万辆,每年消耗燃油超 700×104 t,化
石燃料排出 CO2数量巨大,但因 CO2无污染性,容易
被人们忽略.CO2是影响全球变暖的主要温室气体,
城市中巨大的排放量对环境的影响不容忽视[16-17] .
绿地系统中植物能够吸收和转化 CO2,其间接减排
无疑是解决城市环境问题的有效方式之一.因此,本
文使用具有高时间分辨率的光谱技术连续测定夏季
和冬季中北京市园林科学研究院绿地系统(Beijing
Institute of Landscape Architecture,BILA)以及相邻
的 4环路(the Fourth Ring Road,FRR)空气中的 CO2
浓度、δ13C 和 δ18O 值,并结合 4 环路上的车流量和
气象因子进行分析,获得夏季和冬季测定期间北京
市绿地系统和 4环路上 CO2浓度、δ13C 和 δ18O 值变
化特征及其差异,定量区分夏季和冬季测定期间 4
环路和绿地系统空气中 CO2的来源,并分析车流量
和气象因子对 CO2浓度、δ13 C 和 δ18O 值差异的影
响,确定绿地系统中植物在吸收和消减 CO2方面作
用,以期为评估城市绿地系统生态功能提供基础
数据.
1 研究地区与研究方法
1 1 试验地概况
本文中的绿地系统于北京市园林科学研究院内
(39°58′34″ N,116°27′48″ E),西南方向紧挨 4 环
路,东南方向 300 m外为京密路和机场高速路,正南
方向 50 m外为重要交通枢纽的四元桥,周边交通繁
忙,车流量大(图 1).园林科学研究院绿地长约 550
m,宽约 300 m,占地面积约 15 hm2,其中绿地面积约
占总面积的 90%,种植了 200 余种植物,以乔木、草
坪和花卉为主.4 环路辅路及绿地系统周边的道路
上均种植国槐(Sophora japonica)行道树,院内围栏
内侧种植银杏 ( Ginkgo biloba )、 楸树 ( Catalpa
bungei)、五角枫(Acer mono)和红叶臭椿(Ailanthus
altissima)等高大乔木作为缓冲林带.
100310期 孙守家等: 城市绿地系统内外 CO2浓度、δ13C、δ18O差异和来源及影响因素
图 1 试验位置及取样点
Fig.1 Location of experimental station and sample plots.
将 DLT⁃CO2气体稳定同位素廓线测量系统
(DLT,Los Gatos Research Inc.,USA)放置在距 4 环
路 200 m的绿地中.测定分 2个地点进行,一个进气
口放置在 4环路上的四元西桥立交桥上,一个放置在
距仪器另一侧 100 m的五角枫林里(图 1),高度均为
4 m,使用外径为 1 / 4英寸的 EATON管连接到测量系
统上,EATON管用橡塑海绵保温管包裹以防温度变
化引起同位素分馏,置入 PVC管后埋入地下 20 cm.
1 2 观测和采样方法
1 2 1空气中 CO2浓度、δ13C 和 δ18O 的测定 试验
在 2013年 8月 1—15日和 2014年 2月 22日—3月
8日进行. DLT⁃CO2气体稳定同位素测量系统使用
前,使用已知 CO2浓度、δ13C 和 δ18O 值的标准气体
(北京理加联合科技有限公司,北京)标定,标定时
将气瓶与仪器进气口段连接,调整减压阀使输出压
力为 0. 1 bar,在仪器桌面点击 setup 图标进入系统
设置,在系统设置中选择 Calibrate 进入到 Calibrate
status界面,然后在 Enter Calibration Data 界面将已
知的 CO2浓度、δ13C 和 δ18O 值填入标定栏,通气后
对仪器进行标定,60 s后标定完成,测试期间每周标
定 1次.测定时分别在 2个进气口端装有 6 μm的不
锈钢过滤器,气体进入分析仪前再经过 4.5 μm的聚
偏氟乙烯(PVDF)滤膜过滤掉空气的水汽.真空泵抽
气的流量为 10 L·min-1,16通道的多路器(Los Ga⁃
tos Research Inc.,USA)控制 2 个进口端的气体采样
时间,其中 4环路端距测定系统约 200 m,设定采样
时间为 9 min,前 4 min对管路进行清洗,后 5 min是
测定时间,绿地内端距测定系统 100 m,设定采样间
隔为 6 min,前 2 min对管路进行清洗,后 4 min是测
定时间,每 15 min完成一个采样循环,采样频率为 1
Hz,其中,13 C 测定精度为 0. 2‰,18 O 测定精度为
2‰,CO2测定精度为 0.1 μmol·mol
-1 .13C和18O的测
定结果分别相对于维也纳 Pee Dee箭石标准(Vienna
Pee Dee Belemnite,VPDB)和维也纳标准平均海水
(Vienna Standard Mean Ocean Water,V⁃SMOW)的千
分率(‰)给出,以 δ13C和 δ18O表示,公式为:
δ13C( 18O)=
Rsam
Rsta
-1
æ
è
ç
ö
ø
÷ ×1000 (1)
式中:Rsam和 Rsta分别表示样品和标准物13C / 12C 或
18O / 16O的摩尔比率.
1 2 2汽车尾气样品的采集和测定 用 FM2002⁃
12V型微型真空泵(成都气海机电制造有限公司,中
国)采集大众、丰田、现代等 9 种品牌 21 款汽车尾
气,外径 8 mm、内径 4 mm乳胶管一端连接 0.45 μm
空气过滤器,另外一端与微型真空泵进气口连接,真
空泵出气口与 1 L 铝塑复合膜气体采样袋(大连海
得科技有限公司,中国)连接.在汽车怠速时,发动机
稳定工作后开始对尾气采样,采样时将空气过滤器
置于汽车排气管内,流量设置为 2 L·min-1,采样袋
装满后用旋塞直通阀密封,带回实验室,抽取 10 mL
注射到测定袋中,用纯氮气稀释至 400 μmol·mol-1
左右,记录稀释倍数,将稀释后的气体样品连接到
DLT⁃CO2同位素测定仪上进行测定.
1 2 3车流量测定 在 4 环路的四元西桥立交桥上
架设视频车辆检测系统(赛道科技有限公司,中国)
测定 4环路上车流量,该系统基于 FPGA、双 DSP 的
图像处理内核,内置视频车辆检测算法,同时测定 3
个车道的车流量,每 1 min 输出 1 组平均值.4 环路
为双向 8 车道,将 3 个车道的车流量平均后乘以 8
获得整个 4环路的车流总量.
2003 应 用 生 态 学 报 26卷
1 2 4气象因子测定 在绿地中设定自动气象站,
采用 LI190SB 光合有效辐射传感器(Li⁃Cor Biosci.
Inc.,美国)、HMP45C 型温湿度传感器 ( Campbell
Sci. Inc.,美国)和 AR⁃100 风速仪(Vector Inc. ,英
国)分别获得太阳辐射(Ra)、温度(Ta)、相对湿度
(HR)和风速(V),数据采集器为 CR10X(Campbell
Sci. Inc.,美国),设定每 1 min 采集 1 次,每 10 min
输出 1组平均值.
1 2 5植物和土壤中 δ13C 和 δ18O 值测定 在夏季
观测期间对银杏、楸树、五角枫和红叶臭椿的叶片及
其生长的土壤进行取样用于测定 δ13C和 δ18O值,取
样时,每种植物分别在东南西北 4 个方位取样后混
合,每种植物 3 个重复,土壤样品 6 个重复.样品在
70 ℃烘箱中烘干 48 h,粉碎过 80 目筛,送中国林科
院稳定同位素比率质谱实验室在 Flash EA1112 HT
元素分析仪(Thermo Scientific公司,美国)中高温氧
化还原为 CO2,用质谱仪(DELTA V Advantage,Ther⁃
mo Scientific公司,美国)检测 CO2中的 δ13C 和 δ18O
值,测定精度分别为 0.1‰和 0.2‰.
1 3 数据处理
生态系统 CO2的碳同位素组成计算采用通量⁃
廓线法,假定13 CO2和12 CO2具有相同的湍流扩散
率[9],用同一高度的 CO2浓度和 δ13C值获得 Keeling
plot方程,其截距即为生态系统 CO2的碳同位素组
成 δE [8,18]:
δb =Ca(δa-δE)(1 / Cb)+δE (2)
式中:δb和 Cb分别是生态系统空气中 δ13C 值和 CO2
浓度;Ca和 δa分别是背景大气 δ13C值和 CO2浓度.
城市系统大气中 CO2主要来源于化石燃料燃烧
和自然生态系统异养自养呼吸.其中,4 环路上化石
燃料燃烧产生的 CO2主要来自于汽车尾气,尾气
CO2碳同位素组成( δV)根据 21 款汽车尾气中 δ13C
值和 CO2浓度利用 Keeling plot 方程获得,自然生态
系统呼吸出 CO2的碳同位素组成根据自养呼吸和异
养呼吸的比例来确定.Bowling 等[19]研究证实,冠层
叶片呼出 CO2的碳同位素比率比叶片碳同位素比率
高约 3‰,土壤微生物或其他生物的异养呼吸过程
中稳定同位素的变化虽不能直接测定,但因线粒体
呼吸过程中不发生碳同位素的分馏[20],故土壤碳同
位素组成可作为土壤异养呼吸出 CO2的碳同位素组
成.典型森林生态系统中自养呼吸约占总呼吸的
0 51,草地生态系统中约占 0.40[11],因此在既有树
木又有草地的城市绿地系统自养呼吸占总呼吸的比
例取其均值为 0.45.绿地系统呼吸出 CO2中碳同位
素组成(δT)根据以下公式获得:
δT = 0.45(δp+3‰)+0.55δs (3)
式中:δp为植物呼吸(自养呼吸)的碳同位素组成;δs
为土壤呼吸(异养呼吸)的碳同位素组成.
绿地系统呼出 CO2的 δT与汽车尾气 CO2的 δV之
间显著差异,利用二元线性混合模型来确定二者对
城市生态系统 δE的贡献率(FV,%) [8]:
FV =
δE-δT
δV-δT
×100% (4)
使用 Excel 2010 进行数据整理并作图,用 IBM
SPSS Statistics 20分析软件进行逐步回归分析.
2 结果与分析
2 1 不同季节 4 环路的车流量特征及夏季气象因
子变化
由图 2 可以看出,夏季观测期间 4 环路上每日
6:00—20:00车流总量约 16.5 万辆,冬季观测期间
约 15.7万辆.夏季和冬季观测期间车流量均有明显
的日变化,但存在较大差异.夏季约在 7:30 形成早
高峰,随后车流量略有下降,10:00 左右达到一天中
的峰值,16:30左右形成晚高峰.冬季早高峰在 8:00
左右,14:00到达最高峰(约每小时 6700 辆),晚高
峰出现在 16:00(约每小时 6200 辆).与夏季相比,
冬季车流量的早高峰晚 0.5 h,晚高峰早 0.5 h.
气象因子影响植物对 CO2的吸收和碳氧稳定同
位素的分馏.夏季测定期间,6:00 太阳辐射开始增
加,8:00 开始迅速升高,12:00 太阳辐射达到峰值
(520.35 W·m-2),此后太阳辐射逐渐减小.温度从
6:00开始升高,12:00—18:00均保持在 30 ℃以上,
其中,14:00到达峰值(31.36 ℃).相对湿度从 5:00
开始下降,13:00达到最低值(54.9%)后逐步回升,
日变化呈单波谷曲线.风速变化比较复杂,2:00—
7:00风速相对较低,随后开始增大,16:00 达到峰值
后逐步下降.
2 2 不同季节绿地系统与 4 环路上 CO2浓度日变
化特征及其差异
夏季观测期间,4 环路上空气中 CO2浓度呈现
双峰曲线的日变化,日均值为 427.77 μmol·mol-1,
第一个 CO2浓度高峰出现在 8:30 左右,峰值为
440 01 μmol·mol-1,比日均值高 16 23 μmol·mol-1,
第二高峰出现在 17: 00 左右,峰值约 438 44
μmol·mol-1,比日均值高 10.67 μmol·mol-1 .绿地系
300310期 孙守家等: 城市绿地系统内外 CO2浓度、δ13C、δ18O差异和来源及影响因素
统中 CO2浓度的日均值为 386.44 μmol·mol
-1,第一
个 CO2浓度高峰出现在 8:30 左右,峰值约 397 69
μmol·mol-1,比日均值高 11.25 μmol·mol-1,第二
峰值出现在 17:00,仅比日均值高 1.66 μmol·mol-1 .
冬季观测期间,4 环路和绿地系统中的 CO2浓度较
高,分别为 481.44和 475.23 μmol·mol-1,其日变化
曲线相似:9:00形成一个 CO2浓度高峰,4 环路上和
绿地系统的 CO2 浓度分别为 495. 91 和 483 16
μmol·mol-1,分 别 比 日 均 值 高 14. 47 和 7 93
μmol·mol-1,17:00形成另外一个高峰,峰值分别为
496.15 和 486. 97 μmol·mol-1,分别比日均值高
14 71和 11.74 μmol·mol-1(图 3).
图 2 夏季和冬季观测期间在半小时尺度上 4环路上车流量以及夏季太阳辐射、温度、相对湿度和风速的变化
Fig.2 Variations of traffic volume on the FRR during measurement time in summer and winter seasons at half⁃hour timescale and vari⁃
ations of solar radiation, temperature, relative humidity and wind speed in summer season (mean±SD).
图 3 夏季(A)和冬季(B)观测期间半小时尺度上 CO2浓度的变化、以及 4环路与北京市园林科学研究院绿地系统中 CO2浓
度差异(C)
Fig.3 Variations of CO2 concentration at half⁃hour timescale during measurement time in summer (A) and winter (B) seasons and
the difference of CO2 concentration between the FRR and the green space system of BILA in summer and winter seasons (C).
1 和 2分别代表 4环路和北京市园林科学研究院绿地系统 Number 1 and 2 represented the FRR and the green space system of BILA, respectively. 下
同 The same below.
4003 应 用 生 态 学 报 26卷
夏季和冬季观测期间的 4 环路与绿地系统中
CO2浓度差值存在明显不同.冬季 CO2浓度差值较
小,日均差值仅为 6.21 μmol·mol-1 .第一个差值高
峰在 9:00,为 13.78 μmol·mol-1,第二个差值高峰
在 16:30,为 12.96 μmol·mol-1 .与冬季相比,夏季观
测期间 4环路上 CO2浓度明显高于绿地系统,日均
差值为 41.43 μmol·mol-1,为冬季的 6.67 倍,其中,
8:00开始差值增大,10:00形成差值第一个高峰,为
53.34 μmol·mol-1,随后开始下降,15:00 又开始升
高, 17: 00 形 成 第 二 个 差 值 高 峰, 为 50 56
μmol·mol-1 .
2 3 不同季节绿地系统与 4 环路上 δ13C 日变化特
征及其差异
夏季观测期间,4 环路上空气中 δ13C 值呈现双
波谷曲线的日变化,日均值为-14.7‰,第一个 δ13C
值波谷出现在 8:00左右,最小值为-15.5‰,比日均
值低 0.9‰,第二个波谷出现在 17:00 左右,峰值约
-13.6‰,比日均值低 1.7‰.绿地系统中 δ13C值也呈
双波谷的日变化,日均值为-11.3‰.δ13C 值从 6:00
开始逐步下降,8:00 到达第一个波谷,为-12.0‰,
比日均值低 0. 7‰,第二个波谷出现在 18:00,为
-12.2‰,比日均值低 1.0‰,随后逐步升高.与夏季
相比,冬季观测期间 4环路和绿地系统中的 δ13C 值
较低,分别低 3. 1‰和 6. 3‰,其日变化曲线相似:
9:00形成第一个 δ13 C 值波谷,4 环路和绿地系统
δ13C值分别为-19.4‰和-19.1‰,分别比日均值低
1 7‰和 1.5‰;在 17:00形成第二个波谷,4 环路上
和绿地系统 δ13C 值分别为-19.1‰和-18.5‰,分别
比日均值低 1.3‰和 1.0‰(图 4).
冬季 4环路与绿地系统中 δ13C 差值较小,日均
差值仅为-0.2‰.夜晚多数时间内 4 环路的 δ13C 值
比绿地系统要高,而在白天 6:00—10:00 和14:00—
18:00时间段内 4 环路的 δ13C 值低于绿地系统,第
一个差值波谷在 8:00,为-0.4‰,第二个差值波谷
在 16:00,为-1.0‰.与冬季相比,夏季观测期间 4 环
路的 δ13 C 值明显低于绿地系统,日均差值为
-3.4‰,为冬季的 16.54 倍.其中 9:00 形成 δ13C 差
值第一个波谷为-4.2‰, 16:30 形成第二个差值波
谷为-4.9‰.
2 4 不同季节绿地系统与 4 环路上 δ18O 日变化特
征及其差异
夏季观测期间,4 环路上空气中 δ18O 的日均值
为 35.0‰,δ18O 值从 9:00 开始逐步下降,12:30 降
至最低点,最小值为 31.4‰,比日均值低 3.7‰,随后
开始升高,最高值出现在 17:30,为 38.0‰.绿地系统
中 δ18O值呈单波谷的日变化,日均值为 30.8‰,比 4
环路日均值低 4.2‰.δ18O值从 9:00 开始逐步下降,
11:00达到波谷,为 20.7‰,比日均值低 2.1‰,随后
逐步升高,17:00达到最高峰 34.7‰.冬季观测期间,
4环路和绿地系统中δ18O日均值分别为36 .4‰和
图 4 夏季(A)和冬季(B)观测期间半小时尺度上 δ13C值的变化、以及 4环路与北京市园林科学研究院绿地系统中 δ13C值的
差异(C)
Fig.4 Variations of δ13C value at half⁃hour in during measurement time in summer (A) and winter (B) seasons and the difference of
δ13C value between the FRR and the green space system of BILA in summer and winter seasons (C).
500310期 孙守家等: 城市绿地系统内外 CO2浓度、δ13C、δ18O差异和来源及影响因素
36.1‰,差别较小,其日变化曲线相似:4环路和绿地
系统中 δ18O值分别在 12:00 和 11:30 降至最低点,
最小值分别为 32. 6‰和 29. 6‰,分别比日均值低
3 8‰和 5.2‰(图 5).
冬季观测期间 4 环路与绿地系统中 δ18O 差值
较小,日均差值仅为 0.3‰.夏季观测期间 4 环路上
δ18O值明显高于绿地系统,日均差值为 4.2‰,为冬
季的 12.67倍.其中 7:00 形成 δ18O差值第一个波峰
(6.4‰),随后开始下降, 15: 30 达到最低值,为
-0.02‰,随后震荡升高至 20:30形成第二个差值波
峰(7.1‰).
2 5 不同季节绿地系统和 4 环路 CO2来源组成及
其差异
将 6:00—18:00 绿地系统和 4 环路以及 21 款
汽车尾气中 CO2浓度和 δ13C 值代入 Keeling plot 方
程.当 CO2浓度趋于无穷大即 1 / CO2趋于 0 时,夏季
观测期间白天绿地系统和 4 环路、冬季观测期间白
天绿地系统和 4环路及汽车尾气的 δ13C和 1 / CO2方
程的截距分别为 - 27. 1‰、 - 25. 7‰、 -27.4‰、
-27.0‰和-30.0‰,方程均达到显著水平,表明夏季
和冬季观测期间 4 环路、绿地系统以及汽车燃烧化
石燃料产生 CO2的 δ13C值分别为-27.1‰、-25.7‰、
-27.4‰、-27.0‰和-30.0‰(表 1).夏季,绿地系统
植物生长旺盛,生态系统呼吸来源于自养和异养呼
吸,故生态系统呼吸出 CO2的 δT = 0.45( δp +3‰) +
0 55δs,计算获得 δT值为-21.8‰,冬季大多数植物
进入休眠,在计算生态系统呼吸时忽略自养呼吸,假
设生态系统呼吸出 CO2的 δT等于异养呼吸 δs,为
-19.0‰.
图 6 结果显示,夏季观测期间白天的 4 环路空
气 CO2主要来源于汽车尾气的排放,约占 64.9%,来
源于植物和土壤呼吸占35.1%,绿地系统空气CO2
图 5 夏季(A)和冬季(B)观测期间半小时尺度上 δ18O值的变化、以及 4环路与北京市园林科学研究院绿地系统中 δ18O 值的
差异(C)
Fig.5 Variations of δ18O at half⁃hour during measurement time in summer (A) and winter (B) seasons and the difference of δ18O
value between the FRR and the green space system of BILA in summer and winter seasons (C).
表 1 汽车尾气(δV)的碳同位素平均值、植物(δp)和土壤(δs)碳同位素值、绿地系统呼吸出 CO2的 δT = 0.45(δp +3‰) +
0 55δs和城市生态系统 CO2(δE)的碳同位素值
Table 1 Average isotopic values of CO2 from exhaust gas (δV), average respiration for plant (δp), soil (δs), the average
isotopic value for green space respiration δT =0.45(δp+3‰)+0.55δs and urban ecosystem (δE)
季节
Season
地点
Location
汽车尾气
δV (‰)
植物叶片
δP (‰)
土壤
δs (‰)
绿地系统呼吸
δT (‰)
城市生态系统
δE (‰)
夏季 4环路 FRR -30.0±1.3 -28.2±1.1 -19.0±0.4 -21.8±0.7 -27.1
Summer 绿地系统 BILA -25.7
冬季 4环路 FRR -30.0±1.3 - -19.0±0.4 -19.0±0.4 -27.4
Winter 绿地系统 BILA -27.0
6003 应 用 生 态 学 报 26卷
图 6 利用碳稳定同位素( δ13C)分割获得夏季(A)和冬季
(B)观测期间城市生态系统中绿地系统呼吸和汽车尾气排
放 CO2所占比例
Fig.6 Fractional contribution of green space respiration CO2 to
total CO2 in the urban ecosystem using isotopic partitioning
(δ13C) during measurement time in summer ( A) and winter
(B) seasons.
则主要来源于植物和土壤的呼吸,占 52.7%,来源于
汽车尾气占 47.3%.绿地系统与 4环路空气中 CO2主
要来源不同,绿地系统中来源于汽车尾气的 CO2比
4环路上少了 17.6%.冬季观测期间白天 4环路和绿
地系统空气CO2均主要来源于汽车尾气的排放,分
别占 76.8%和 73.3%,而来源于土壤呼吸的比例仅
为 23.2%%和 26.7%.绿地系统与 4 环路空气中 CO2
主要来源相同且差异较小,仅差 3.5%.
2 6 不同季节 CO2、δ13C 和 δ18O 差值与车流量的
关系
由图 7 可以看出,夏季观测期间 4 环路与绿地
系统中 CO2浓度差值和 δ13C 差值均与车流量呈极
显著相关关系(P<0.01),δ18O差值与车流量呈显著
相关关系(P<0.05),夏季植物生长旺盛,绿地系统
中高大乔木的阻挡和消减作用,导致绿地系统与 4
环路上的 CO2浓度差值 δ13C 和 δ18O 差值较大.冬季
观测期间 4环路与绿地系统中 CO2浓度、δ13C 和 δ18
O差值均与车流量相关关系不显著(P>0.05),可能
原因是冬季植物落叶,没有林冠的阻挡,CO2扩散较
快,且空气中 CO2主要来源于化石燃料燃烧,来源单
一,故冬季中 4环路与绿地系统中 CO2、δ13C 和 δ18O
差值与车流量无显著相关性.
2 7 夏季 CO2、δ13C 和 δ18O 差值与车流量和环境
因子的关系
夏季中植物生长和车流量均会影响到空气中
CO2变化及其同位素组成.为了评估绿地系统和车流
量对 CO2浓度、δ13C 和 δ18O 值的影响,建立逐步回
归模型来分析车流量和绿地植物生长活动与 CO2浓
度、δ13C和 δ18O差值的关系,选择车流量(x1)、太阳
辐射(x2)、温度(x3)、相对湿度(x4)和风速( x5)作
为影响其 CO2浓度 (y1)、δ13C(y2)和 δ18O(y3)差值
变化的因子,获得逐步回归模型(表 2).
图 7 夏季(A)和冬季(B)观测期间 4环路与北京市园林科学研究院绿地系统中 CO2、δ13C和 δ18O差值与车流量关系
Fig.7 CO2 concentration, δ13C value and δ18O value difference between FRR and the green space system of BILA in relation to traffic
volume during measurement time in summer (A) and winter (B).
700310期 孙守家等: 城市绿地系统内外 CO2浓度、δ13C、δ18O差异和来源及影响因素
表 2 夏季观测期间 4 环路和绿地系统 CO2 ( y1 )、δ13 C
(y2)、δ18O(y3)差值与车流量(x1)和气象因子(太阳辐射
x2、温度 x3、相对湿度 x4和风速 x5)的逐步回归分析
Table 2 Stepwise regression analysis of the difference of
CO2 concentration ( y1), δ13C ( y2) and δ18O ( y3 ) value
between the FRR and green space system in relation to tra⁃
ffic volume (x1) and meteorological factors [( solar radia⁃
tion (x2), temperature (x3), relative humidity (x4) and
wind speed (x5)] during measurement time in summer
逐步回归方程
Stepwise regression equation
相关系数
R
决定系数
R2
P
y1 = 2.110+0.008x1-0.007x2 0.687 0.472 0.000
y2 =-7.585+0.002x2+0.601x4 0.631 0.398 0.001
y3 =-3.403+0.111x4 0.701 0.491 0.000
由表 2 可知,3 个逐步回归方程均达到了极显
著水平,相关系数分别为 0.687、0.631 和 0.701,其
中,车流量和太阳辐射对 CO2浓度差值产生显著影
响,光照是影响植物光合作用的关键因子,除了车辆
直接排放 CO2造成 4环路与绿地系统 CO2浓度差异
外,植物光合作用吸收 CO2也会导致绿地内外的差
异增大.太阳辐射和相对湿度对 δ13C 差值产生显著
影响,则可能是与光合过程中对碳同位素分馏不同
造成的.δ18O差值则受到相对湿度的显著影响,可能
是由于水汽中也含有18O,光合过程中 CO2中的18O
与水汽中的18O 进行交换,从而受到相对湿度的影
响较大.冬季中 CO2浓度、δ13C 和 δ18O 差值较小,与
车流量和气象因子的逐步回归方程不成立(结果未
列出).
3 讨 论
城市生态系统碳氧循环不同于自然生态系统,
受到人为和自然因素的共同影响.本文结果显示,夏
季和冬季观测期间北京市 4环路上的 CO2浓度分别
为 427.77 和 481.44 μmol·mol-1,高于北京市上甸
子区域大气本底监测站监测的 CO2本底浓度,这与
王长科等[21]观测到北京城市大气 CO2浓度明显偏
高的结果一致.4环路上车流量日均超过 15万辆,汽
车在怠速时排出的尾气中 CO2浓度为 69141. 72
μmol·mol-1(结果未列出),为本底数据的 178.18
倍,大量含高浓度 CO2的汽车尾气进入空气中导致
空气中 CO2浓度增加,CO2与城市的交通量有显著
相关关系[13,16] .4环路上 CO2浓度呈双峰曲线变化,
其峰值出现时间与早晚交通高峰时间基本重合,
Song等[22]研究显示,北京的交通状况与 CO2通量变
化趋势基本一致,Hiller等[23]发现,穿过明尼苏达州
草坪的国道上车流量对于观测塔上的 CO2浓度影响
较大.化石燃料燃烧也会导致 CO2中碳氧同位素组
成发生相应改变,北京空气中 δ13C 值呈双波谷曲线
日变化,出现波谷的时间与交通早晚高峰一致,汽车
尾气中 δ13C值较低,导致城市空气 δ13C值低于单一
的生态系统[10,12] .Kèlomé 等[24]研究发现,在科托努
市中交通较少的区域空气中 δ13C 值较高,交通繁忙
地区的 δ13C 值较低,在空间分布上除交通量外,还
与风向有关,顺风和逆风的 Δ13C 与 CO2浓度差值呈
线性关系[17] .同位素定量区分结果显示,4 环路上汽
车尾气排放的 CO2占总量的 64.9%,Wada 等[15]研
究发现,春季名古屋空气中 CO2有 33%来源于汽油
燃烧.夏季北京 4 环路空气中 δ13 C 值比冬季高,
Pataki等[25]对美国盐湖城研究显示,汽油和天然气
燃烧以及生物过程对夜晚城市 δ13C 值的影响存在
明显的季节变化,春季和夏季明显高于冬季.同位素
区分结果显示,冬季化石燃料燃烧产生 CO2占的比
例超过 73%,这可能与冬季供暖增加了 CO2的排放
有关,中纬度的城市如日本东京[26]、波兰罗兹[27]等
也表现出相同趋势,这些研究均表明化石燃料燃烧
对城市环境产生巨大影响.
城市绿地植物具有吸收 CO2维持城市碳平衡的
生态功能,在减缓城市环境压力方面发挥重要作用.
本文结果显示,夏季观测期间绿地系统中 CO2浓度
远低于 4环路,但冬季绿地系统 CO2浓度与 4 环路
相差很小.夏季绿地系统中 CO2减少的其原因可能
有两个:1)夏季,绿地系统中高大乔木形成的景观
带起到了缓解、消减尾气的作用[28],Gu 等[29]研究
认为,相比于裸露街谷,植物的引入减弱了街谷顶部
气体置换速率,而汽车尾气在穿过绿地系统密集多
层次的高大乔木时,由于乔木枝叶的摩擦、阻挡作
用,有效地阻滞气体的水平扩散,4 环路上汽车产生
的 CO2难于水平扩散直接进入绿地系统内部.冬季
由于乔木落叶,4 环路上空气比较容易水平穿过林
带,直接与绿地系统气体进行交换,导致绿地系统与
4环路上空气中的 CO2浓度差别不明显.2)夏季,植
物通过光合作用来同化空气中 CO2,导致绿地系统
中 CO2浓度下降,而冬季植物落叶后无法通过光合
作用来吸收 CO2,因而与 4 环路的浓度差仅为 6.21
μmol·mol-1 .Nordbo等[30]研究发现,当城市区域植
被覆盖率达 80%左右,碳源汇特征表现为中性,而
在伦敦、墨西哥城、东京及墨尔本等城市的研究表
明,城市内年 CO2排放与植被覆盖率呈现一定的负
相关关系[31] .在光合作用过程中,叶片优先吸收12C,
8003 应 用 生 态 学 报 26卷
引起碳稳定同位素分馏[32],夏季绿地系统白天 δ13C
值较高,表明光合作用导致绿地空气中的13C 富集.
同位素区分结果显示,夏季观测期间绿地系统空气
中来源于化石燃料燃烧产生的 CO2比 4 环路上低
17.6%,而冬季植物落叶后仅仅低 3.5%,表明城市
绿地系统在生长季节时降低石化燃料燃烧 CO2浓度
方面发挥了重要作用,这与 Nowak等[33]发现城市树
木在降低 CO2浓度方面比非城市树木的作用更为显
著的结果一致.
夏季,绿地生态系统空气中 CO2浓度、δ13 C 和
δ18O值与道路上存在明显不同,逐步回归分析显示
城市系统碳循环受到自然因素和人类活动的综合影
响.除了车流量直接影响 CO2排放量外,光照通过影
响植物光合作用来同化 CO2也会导致绿地系统与 4
环路差异增大.在半小时尺度上,太阳辐射和相对湿
度是影响绿地系统与 4 环路 δ13C 和 δ18O 差值的主
要环境因子,其通过增强绿地植物碳循环过程来发
挥作用[34] .相似研究显示,亚高山带针叶林冠层内
外 δ13C差值与 VPD有显著的相关关系[35],而 VPD
是由温度和湿度计算获得的.Flanagan 等[36]将微气
象测量与18O 同位素组成相结合,发现光合作用显
著地影响森林生态系统18O变化.城市环境高度异质
性,使得城市系统比结构单一的生态系统更复杂,受
到影响因素更多,同时城市生态系统中石化燃料燃
烧、生物质燃烧、人类呼吸和自然界自养异养呼吸等
多个碳源汇及其易变性使得同位素方法还存在一定
程度的不确定性,需要结合14C、18O 以及 Isosurce 模
型等在进一步的研究中去精确量化和验证.
4 结 论
夏季和冬季观测期间,北京市 4 环路上车流量
较大,有明显的早晚交通高峰.汽车排放出大量
CO2,导致 4环路及其相邻的绿地系统空气中 CO2浓
度呈双峰曲线、δ13C值呈双波谷曲线和 δ18O 值呈单
波谷曲线的日变化,且夏季中二者空气中的 CO2浓
度、δ13C和 δ18O差值较大.夏季 4环路空气中 CO2主
要来源于化石燃料燃烧而绿地系统空气中 CO2主要
来源于生态系统的自养和异养呼吸,冬季二者空气
中 CO2均主要来源于化石燃料燃烧.除了车流量外,
在半小时尺度上太阳辐射和相对湿度是影响 4环路
和绿地系统 CO2浓度、δ13C 和 δ18O 差值的主要环境
因子.绿地系统中的植物在生长季节通过光合作用
同化 CO2来维持城市碳平衡,在改善城市生态环境
方面发挥了重要作用.
参考文献
[1] Churkina G. Modeling the carbon cycle of urban sys⁃
tems. Ecological Modelling, 2008, 216: 107-113
[2] Zhao R⁃Q (赵荣钦), Huang X⁃J (黄贤金). Carbon
cycle of urban system: Characteristics, mechanism and
theoretical framework. Acta Ecologica Sinica (生态学
报), 2013, 33(2): 358-366 (in Chinese)
[3] Lin G⁃H (林光辉). Stable Isotope Ecology. Beijing:
Higher Education Press, 2013 (in Chinese)
[4] Bowling DR, Pataki DE, Richardson JT. Carbon iso⁃
topes in terrestrial ecosystem pools and CO2 fluxes. New
Phytologist, 2008, 178: 24-40
[5] Wang Z⁃P (王智平), Chen Q⁃S (陈全胜). Recently
photosynthesized carbon allocation and turnover: A mi⁃
nor review of the literature. Acta Phytoecologica Sinica
(植物生态学报), 2005, 29(5): 845-850 ( in Chi⁃
nese)
[6] Kodama N, Ferrio JP, Brüggemann N, et al. Short⁃term
dynamics of the carbon isotope composition of CO2 emit⁃
ted from a wheat agroecosystem: Physiological and envi⁃
ronmental controls. Plant Biology, 2011, 13: 115-125
[7] Wingate L, Ogée J, Burlett R, et al. Photosynthetic car⁃
bon isotope discrimination and its relationship to the car⁃
bon isotope signals of stem, soil and ecosystem respira⁃
tion. New Phytologist, 2010, 188: 576-589
[8] Yakir D, Sternberg LDL. The use of stable isotopes to
study ecosystem gas exchange. Oecologia, 2000, 123:
297-311
[9] Lee X, Kim K, Smith R. Temporal variations of the
18O / 16O signal of the whole⁃canopy transpiration in a
temperate forest. Global Biogeochemical Cycles, 2007,
21, Doi: 10.1029 / 2006GB002871
[10] Fassbinder JJ, Griffis TJ, Baker JM. Interannual, sea⁃
sonal, and diel variability in the carbon isotope composi⁃
tion of respiration in a C3 / C4 agricultural ecosystem. Ag⁃
ricultural and Forest Meteorology, 2012, 153: 144-153
[11] Zhu X⁃J (朱先进), Yu G⁃R (于贵瑞), Wang Q⁃F
(王秋凤), et al. The interaction between components of
ecosystem respiration in typical forest and grassland eco⁃
systems. Acta Ecologica Sinica (生态学报), 2013, 33
(21): 6925-6934 (in Chinese)
[12] Goffin S, Aubinet M, Maier M, et al. Characterization
of the soil CO2 production and its carbon isotope compo⁃
sition in forest soil layers using the flux⁃gradient ap⁃
proach. Agricultural and Forest Meteorology, 2014,
188: 45-57
[13] Contini D, Donateo A, Elefante C, et al. Analysis of
particles and carbon dioxide concentrations and fluxes in
an urban area: Correlation with traffic rate and local mi⁃
crometeorology. Atmospheric Environment, 2012, 46:
25-35
[14] Liu W (刘 卫),Wei N⁃N (位楠楠),Wang G⁃H (王
广华), et al. Quantitative estimation source of urban at⁃
mospheric CO2 by carbon isotope composition. Environ⁃
mental Science (环境科学), 2012, 33 ( 4): 1041 -
1049 (in Chinese)
[15] Wada R, Pearce JK, Nakayama T, et al. Observation of
carbon and oxygen isotopic compositions of CO2 at an ur⁃
ban site in Nagoya using Mid⁃IR laser absorption spec⁃
900310期 孙守家等: 城市绿地系统内外 CO2浓度、δ13C、δ18O差异和来源及影响因素
troscopy. Atmospheric Environment, 2011, 45: 1168 -
1174
[16] Velasco E, Pressley S, Grivicke R, et al. Eddy cova⁃
riance flux measurements of pollutant gases in urban
Mexico City. Atmospheric Chemistry and Physics Discus⁃
sions, 2009, 9: 7325-7342
[17] Clark⁃Thorne ST, Yapp CJ. Stable carbon isotope con⁃
straints on mixing and mass balance of CO2 in an urban
atmosphere: Dallas metropolitan area, Texas, USA. Ap⁃
plied Geochemistry, 2003, 18: 75-95
[18] Keeling CD. The concentration and isotopic abundances
of atmospheric carbon dioxide in rural areas. Geochimica
et Cosmochimica Acta, 1958, 13: 322-334
[19] Bowling DR, Pataki DE, Kanderson JT. Carbon isotopes
in terrestrial ecosystem pools and CO2 fluxes: A review.
New Phytologist, 2008, 178: 24-40
[20] Lin GH, Ehleringer JR. Carbon isotopic fractionation
does not occur during dark respiration in C3 and C4
plants. Plant Physiology, 1997, 114: 391-394
[21] Wang C⁃K (王长科), Wang Y⁃S (王跃思), Liu G⁃R
(刘广仁). Characteristics of atmospheric CO2 variations
and some affecting factors in urban area of Beijing. Envi⁃
ronmental Science (环 境 科 学 ), 2003, 24 ( 4 ):
13-17 (in Chinese)
[22] Song T, Wang Y. Carbon dioxide fluxes from an urban
area in Beijing. Atmospheric Research, 2012, 106:
139-149
[23] Hiller R, Zeeman MJ, Eugster W. Eddy⁃covariance flux
measurements in the complex terrain of alpine valley in
Switzerland. Boundary⁃Layer Meteorology, 2008, 127:
449-467
[24] Kèlomé NC, Lévêque J, Andreux F, et al. C4 plant iso⁃
topic composition (δ13C) evidence for urban CO2 pollu⁃
tion in the city of Cotonou, Benin (West Africa). Sci⁃
ence of The Total Environment, 2006, 366: 439-447
[25] Pataki DE, Bowling DR, Ehleringer JR. Seasonal cycle
of carbon dioxide and its isotopic composition in an
urban atmosphere: Anthropogenic and biogenic effects.
Journal of Geophysical Research, 2003, 108, Doi: 10.
1029 / 2003JD003865
[26] Moriwaki R, Kanda M. Seasonal and diurnal fluxes of
radiation, heat, water vapor, and carbon dioxide over a
suburban area. Journal of Applied Meteorology, 2004,
43: 1700-1710
[27] Pawlak W, Fortuniak K, Siedlecki M. Carbon dioxide
flux in the centre of Lodz, Poland: Analysis of a 2⁃year
eddy covariance measurement data set. International
Journal of Climatology, 2011, 31: 232-243
[28] Chen B⁃F (陈步峰), Pan Y⁃J (潘勇军), Xiao Y⁃H
(肖 以 华 ), et al. Characteristics of environmental
effects in the landscape shelter forest along the aero⁃
drome road in Guangzhou. Forest Research (林业科学研
究), 2012, 25(3): 278-284 (in Chinese)
[29] Gu ZL, Zhang YW, Lei KB. Large eddy simulation of
flow in a street canyon with tree planting under various
atmospheric instability conditions. Science China Techno⁃
logical Sciences, 2010, 53: 1928-1937
[30] Nordbo A, Jarvi L, Haapanala S, et al. Fraction of nat⁃
ural area as main predictor of net CO2 emissions from
cities. Geophysical Research Letters, 2012, 39, Doi: 10.
1029 / 2012GL053087
[31] Liu M (刘 敏), Fu Y⁃L (伏玉玲), Yang F (杨
芳). Research progress on urban carbon fluxes based on
eddy covariance technique. Chinese Journal of Applied
Ecology (应用生态学报), 2014, 25(2): 611- 619
(in Chinese)
[32] Zou L⁃L (邹绿柳), Sun G⁃C (孙谷畴), Zhao P (赵
平), et al. Net CO2 exchange and carbon isotope flux in
Acacia mangium plantation. Chinese Journal of Applied
Ecology ( 应 用 生 态 学 报 ), 2009, 20 ( 11 ):
2594-2602 (in Chinese)
[33] Nowak DJ, Crane DE. Carbon storage and sequestration
by urban trees in the USA. Environmental Pollution,
2002, 112: 381-389
[34] Barr JG, Engel V, Smith TJ, et al. Hurricane disturb⁃
ance and recovery of energy balance, CO2 fluxes and
canopy structure in a mangrove forest of the Florida
Everglades. Agricultural and Forest Meteorology, 2012,
153: 54-66
[35] Schaeffer SM, Anderson DE, Burns SP, et al. Canopy
structure and atmospheric flows in relation to the δ13C of
respired CO2 in a subalpine coniferous forest. Agricultu⁃
ral and Forest Meteorology, 2008, 148: 592-605
[36] Flanagan LB, Brooks JR, Varney GT, et al. Discrimina⁃
tion against concentration and stable oxygen isotope ratio
of atmospheric CO2 in boreal forest ecosystems. Global
Biogeochemical Cycles, 1997, 11: 83-98
作者简介 孙守家,男,1975 年生,博士,副研究员. 主要从
事生态学研究. E⁃mail: ssj1011@ 163.com
责任编辑 杨 弘
0103 应 用 生 态 学 报 26卷