采用离轴积分腔输出光谱技术测定华北低丘山区栓皮栎人工林冠层上缘(11 m)和下部(6 m)大气CO2浓度和δ13C值,在小时尺度上分析冠层CO2浓度和δ13C变化及其影响因素.结果表明: 冠层CO2浓度呈先高后低再升高的日变化趋势,而δ13C值没有明显一致的日变化规律.白天大气不稳定状态出现的频率为70.2%,在光合作用和林内湍流的共同作用下,栓皮栎冠层下部CO2浓度高于冠层上缘约1.70 μmol·mol-1,而δ13C值低于冠层上缘约0.81‰.晚上大气稳定状态出现的频率为76.2%,湍流弱,冠层叶片呼出的CO2不易流动,导致冠层下部CO2浓度高于上缘约1.24 μmol·mol-1,δ13C值低于冠层上缘约0.58‰.白天和晚上冠层上下缘的CO2浓度差值与δ13C差值均呈显著的相关关系.逐步回归分析表明,白天太阳辐射和相对湿度是影响冠层CO2浓度和δ13C值差异的主要环境因子,晚上温度显著影响冠层下部与上缘δ13C值的变化,这些环境因子通过增强或减弱光合和呼吸作用来影响冠层大气中CO2浓度和δ13C值的变化.
The offaxis integrated cavity output spectroscopy technique was used to measure air CO2 concentration and stable carbon isotope ratio (δ13C) above (11 m) and at the bottom (6 m) of canopy of a Quercus variabilis plantation in a low hilly area of North China. The variations of CO2 concentration and δ13C value in Q. variabilis plantation canopy and the influencing factors were analyzed at hourly timescale. The results showed that diurnal variation in the CO2 concentration had a trend, while there was no obvious similar tendency in the diurnal change of δ13C value. Daytime atmosphere stability frequency during the experiment time was 70.2%. With the combined effects of photosynthesis and turbulent in the canopy, CO2 concentration at the bottom of canopy was 1.70 μmol·mol-1 higher than that above the canopy, while the δ13C value was 0.81‰ lower than that above the canopy. Atmosphere stability frequency was 76.2% at night. The CO2 from leaf was not easy to move because of weak turbulent. Thus, CO2 concentration at the bottom of canopy was 1.24 μmol·mol-1 higher than that above canopy, while the δ13C value was 0.58‰ lower than that above canopy. The difference of CO2 concentration between above and at the bottom of the canopy was strongly correlated with their δ13C difference both in daytime and at nighttime. Stepwise regression analysis indicated that solar radiation and relative humidity in daytime were the main environmental factors causing CO2 concentration and δ13C difference between above and at the bottom of the canopy, whereas at nighttime temperature was a key environmental factor influencing δ13C value. The above environmental factors strongly influenced CO2 concentration and δ13C value in air above and at the bottom of Q. variabilis plantation canopy by increasing or decreasing photosynthesis and respiration.
全 文 :华北低丘山区栓皮栎人工林冠层 CO2浓度
和 δ13C变化及其影响因素∗
孙守家1,2∗∗ 孟 平1,2 张劲松1,2 郑 宁1,2 何春霞1,2 李岩泉1,2
( 1中国林业科学研究院林业研究所国家林业局林木培育重点实验室, 北京 100091; 2南京林业大学南方现代林业协同创新中
心, 南京 210037)
摘 要 采用离轴积分腔输出光谱技术测定华北低丘山区栓皮栎人工林冠层上缘(11 m)和
下部(6 m)大气 CO2浓度和 δ13C值,在小时尺度上分析冠层 CO2浓度和 δ13C变化及其影响因
素.结果表明: 冠层 CO2浓度呈先高后低再升高的日变化趋势,而 δ13C值没有明显一致的日变
化规律.白天大气不稳定状态出现的频率为 70.2%,在光合作用和林内湍流的共同作用下,栓
皮栎冠层下部 CO2浓度高于冠层上缘约 1.70 μmol·mol
-1,而 δ13C值低于冠层上缘约 0.81‰.
晚上大气稳定状态出现的频率为 76.2%,湍流弱,冠层叶片呼出的 CO2不易流动,导致冠层下
部 CO2浓度高于上缘约 1.24 μmol·mol
-1,δ13C 值低于冠层上缘约 0.58‰.白天和晚上冠层上
下缘的 CO2浓度差值与 δ13C差值均呈显著的相关关系.逐步回归分析表明,白天太阳辐射和
相对湿度是影响冠层 CO2浓度和 δ13C值差异的主要环境因子,晚上温度显著影响冠层下部与
上缘 δ13C值的变化,这些环境因子通过增强或减弱光合和呼吸作用来影响冠层大气中 CO2浓
度和 δ13C值的变化.
关键词 栓皮栎; 稳定碳同位素; CO2浓度; 冠层; 光合作用
∗林业公益性行业科研专项经费项目(201404206)和南京林业大学南方现代林业协同创新中心项目资助.
∗∗通讯作者. E⁃mail: ssj1011@ 163.com
2014⁃06⁃18收稿,2014⁃12⁃09接受.
文章编号 1001-9332(2015)02-0370-09 中图分类号 Q948.1 文献标识码 A
Variations of CO2 concentration and δ13C and influencing factors of Quercus variabilis plan⁃
tation in low hilly area of North China. SUN Shou⁃jia1,2, MENG Ping1,2, ZHANG Jin⁃song1,2,
ZHENG Ning1,2, HE Chun⁃xia1,2, LI Yan⁃quan1,2 ( 1Key Laboratory of Tree Breeding and Cultiva⁃
tion, State Forestry Administration, Research Institute of Forestry, Chinese Academy of Forestry, Bei⁃
jing 100091, China; 2Collaborative Innovation Center of Sustainable Forestry in Southern China,
Nanjing Forestry University, Nanjing 210037) . ⁃Chin. J. Appl. Ecol., 2015, 26(2): 370-378.
Abstract: The off⁃axis integrated cavity output spectroscopy technique was used to measure air CO2
concentration and stable carbon isotope ratio (δ13C) above (11 m) and at the bottom (6 m) of
canopy of a Quercus variabilis plantation in a low hilly area of North China. The variations of CO2
concentration and δ13C value in Q. variabilis plantation canopy and the influencing factors were ana⁃
lyzed at hourly timescale. The results showed that diurnal variation in the CO2 concentration had a
trend, while there was no obvious similar tendency in the diurnal change of δ13C value. Daytime at⁃
mosphere stability frequency during the experiment time was 70.2%. With the combined effects of
photosynthesis and turbulent in the canopy, CO2 concentration at the bottom of canopy was 1.70
μmol·mol-1 higher than that above the canopy, while the δ13C value was 0.81‰ lower than that
above the canopy. Atmosphere stability frequency was 76.2% at night. The CO2 from leaf was not
easy to move because of weak turbulent. Thus, CO2 concentration at the bottom of canopy was 1.24
μmol·mol-1 higher than that above canopy, while the δ13C value was 0.58‰ lower than that above
canopy. The difference of CO2 concentration between above and at the bottom of the canopy was
strongly correlated with their δ13C difference both in daytime and at nighttime. Stepwise regression
analysis indicated that solar radiation and relative humidity in daytime were the main environmental
factors causing CO2 concentration and δ13C difference between above and at the bottom of the canopy,
应 用 生 态 学 报 2015年 2月 第 26卷 第 2期
Chinese Journal of Applied Ecology, Feb. 2015, 26(2): 370-378
whereas at nighttime temperature was a key environmental factor influencing δ13C value. The above
environmental factors strongly influenced CO2 concentration and δ13C value in air above and at the
bottom of Q. variabilis plantation canopy by increasing or decreasing photosynthesis and respiration.
Key words: Quercus variabilis; stable carbon isotope; CO2 concentration; canopy; photosynthesis.
森林生态系统每年固定碳约占整个陆地生态系
统的三分之二[1],在调节全球碳循环过程中起着巨
大的作用[2] .植物通过光合作用固定碳,经呼吸作用
再次释放 CO2,其碳循环能否平衡取决于 CO2 的流
向[3] .CO2 的流向与植物新陈代谢及其环境因子密
切相关[4],所以在较高时空尺度上了解植物生理活
动对 CO2 交换过程的影响对于研究植物碳循环非
常重要.尽管 Keeling 曲线[5-6]、流量梯度[7]、涡度相
关[8-9]和箱体法[10]被广泛用于碳循环研究,并成功
应用到生态系统和全球尺度上[11-12],但在呼吸过
程[13]、生态系统 CO2 来源[14]、光合和呼吸作用 CO2
流向[15]等方面仍然存在分歧,其主要原因之一是在
CO2 固定和释放过程中发生诸多生化反应.碳稳定
同位素被认为是理解叶片到生态系统碳交换过程及
其分析光合和呼吸途径中生理生化反应的有效工
具[16],已经被应用到光合过程、生态系统碳转
移[4,17]和其他生态学[18-19]研究中.不过,早期研究多
假定大气碳同位素比率或生态系统中植物组织碳同
位素比率不随时间变化而改变[20],事实上13C 明显
受到环境和生物因素的调节[21],在较短时间尺度上
存在很大差异,假定其保持恒定易掩盖生态系统与
大气之间的碳同位素失衡现象[22] .冠层是林木进行
生理活动的主要层次,目前,对冠层碳循环的研究多
集中于林冠层 CO2 交换特征[23]及其与大气交换过
程的模拟[24],而对于冠层内外13C 的时空变化研究
较少[9,12],主要是因为先前的研究中测定空气中碳
同位素多用质谱仪,对空气中碳同位素只能短暂或
离散取样,时间分辨率较低.在 CO2 纯化过程中,需
要使用液氮和 3级串联式冷阱,进行 2次 CO2 捕集,
操作成本较高,纯化过程中容易产生不完全捕集,易
影响测定结果.近年来,离轴积分腔输出光谱技术
( off⁃axis integrated cavity output spectroscopy, OA⁃
ICOS)的发展使得对空气中碳同位素连续观测成为
可能,其响应时间短,能在野外较高时间分辨率连续
观测和进行大量样品分析.
栓皮栎(Quercus variabilis)人工林是华北低丘
山区主要林型,对当地生态环境具有重要的调节作
用.在较短时间尺度上研究其冠层的碳循环,有助于
揭示栓皮栎人工林的生态功能.因此,本文使用具有
高时间分辨率的光谱技术连续测定栓皮栎冠层内外
的 CO2 和 δ13C 值,并结合复杂的大气环境进行分
析,主要研究目标是:1)获得冠层内外 CO2 浓度和
δ13C的差异,2)分析复杂大气环境下栓皮栎光合和
呼吸过程对冠层 δ13C的影响,3) 确定影响冠层 CO2
和 δ13C 差异的主要环境因子,为进一步研究华北低
丘山区人工林生态系统的生态服务功能提供基础.
1 研究地区与研究方法
1 1 研究区概况
试验于 2011年 5月 15—31日在中国林业科学
研究院黄河小浪底森林生态站进行,该站隶属于中
国森林生态系统定位研究网络(CFERN),地处河南
省济源市境内太行山南段与黄河流域的交界处,定
位站中心地区海拔 410 m.属暖温带大陆性季风气
候,年均温 12.4 ~ 14.3 ℃,年日照时数 2367.7 h,年
日照率 54%,年均降水量 641.7 mm,其中,6—9月为
438.0 mm,占全年的 68.3%.同位素廓线测量系统和
通量观测塔位于东西走向的低山山脊,顶部地势平
坦,核心区周围 1 km2范围内植被覆盖率约为 96%,
其中林木盖度约占总盖度的 71.7%,林木中以栓皮
栎为主,占 94.2%,其余为侧柏(Platycladus orienta⁃
lis)和刺槐(Robinia pseudoacacia),栓皮栎林分密度
约 1850 株·hm-2,叶面积指数 ( LAI)为 ( 2. 96 ±
0 02),郁闭度为 0. 8,林龄为 35 a,平均胸径为
(11.5±1.9) cm,平均树高为(10.7±1.1) m.
1 2 研究方法
1 2 1 CO2 浓度和同位素组成 δ13 C 测定 利用
DLT CO2 气体稳定同位素廓线测量系统(DLT,Los
Gatos Research Inc.,USA)连续观测 CO2 浓度及其同
位素组成( 13C / 12C)的变化,该系统基于 OA⁃ICOS技
术,测定精度为 0.1‰,使用前用已知 CO2 浓度和
δ13C值的标准气体(北京理加联合科技有限公司,北
京)进行标定.测定分为 2 层进行,高度分别为 6 m
(冠层下部,CB)和 11 m(冠层上缘,CT),在气管的
进口端装有 6 μm 的不锈钢过滤器,进入分析仪前
再经过 4.5 μm的过滤器过滤.真空泵抽气,16 通道
的多路器(Los Gatos Research Inc.,USA)控制不同
高度气体采样时间,每层取样 3 min,前 1 min 对管
1732期 孙守家等: 华北低丘山区栓皮栎人工林冠层 CO2浓度和 δ13C变化及其影响因素
路进行清洗,后 2 min 是测定时间,每 6 min 完成一
个采样循环,采样频率为 1 Hz. δ13 C 值的输出以
VPDB(Vienna Pee Dee Belemnite)为标准.
δ13C=(Rsam / Rsta-1)×1000 (1)
式中:Rsam和 Rsta分别表示样品和标准物13C / 12C 的
摩尔比率.
1 2 2环境因子测定 在同一观测点装有一套涡度
相关系统,HMP45C 型温湿度传感器(Vaisala,Hel⁃
sinki,Finland)、CSAT3 三维超声风速仪( Campbell
Sci. Inc., USA)和 AR⁃100 型 3 杯风速计 ( Vector
Instruments,UK)可以获得距地表 6、11和 26 m高度
处的空气相对湿度、温度、风速和风向量,并由
CR5000数据采集器(Campbell Sci. Inc., USA)记录
和保存,采样频率为 10 Hz.
1 3 数据处理
栓皮栎人工林及其生理活动与环境之间的相互
作用发生在冠层,其 CO2 输送与冠层大气的稳定状
态有直接关系,对大气稳定度进行分析有助于确定
冠层内 CO2 的流动方向,稳定度参数用 Monin⁃
Obkhov法获得:
Ç=( z-d0) / L (2)
式中:无量纲参数 Ç为大气稳定度参数;d0为零平面
位移高度;d0 = 0.63 h;h 为树平均高度 11 m;L 是
Monin⁃Obukhov长度.Ç<-1为不稳定状态,Ç>1为稳
定状态,-1≤Ç≤1为中性.
使用 Excel 2010 作图,用 PASW Statistics 18 进
行相关性和逐步回归分析(α= 0.05).
2 结果与分析
2 1 测定期间的环境因子和大气稳定度频率变化
由图 1可以看出,5月 15—31日,除了 5 月 19、
21、28、29和30日,其他天气均为晴天,日际间的太
图 1 试验期间太阳辐射、温度、相对湿度、大气压、风速及大气稳定度频率变化
Fig.1 Variations of solar radiation, temperature, relative humidity, atmospheric pressure, wind speed and atmosphere stability
frequency during the experiment.
Ⅰ: 稳定 Stable; Ⅱ: 中性 Neutral; Ⅲ: 不稳定 Unstable.
273 应 用 生 态 学 报 26卷
阳辐射最大值为 600~900 W·m-2 .日际间气温相差
比较大,5月 17日日均温度最高,为 30.4 ℃,5月 20
日日均温度降至最低值 14.8 ℃,随后上升至 5月 31
日的 27.5 ℃,总体呈现高低高的变化趋势.相对湿
度变幅较大,5 月 17 日相对湿度最低,为 21.7%,5
月 20日相对湿度最高,为 76.5%,总体呈现低高低
的变化趋势.大气压呈现先降后升的变化趋势,5 月
16—18日大气压明显低于其他日期.日际间的风速
差别较大,5月 18、25、26和 30日风速较大.总体上,
测定期间的环境因子差异大,环境条件相对复杂.
从 8:00开始栓皮栎人工林林冠上大气不稳定
状态频率开始迅速增加,在 13:00—14:00 达到峰
值,随后平缓下降,到 17:00 迅速下降,19:00 后大
气不稳定态频率趋于稳定.与不稳定状态变化相反,
林冠上的大气稳定状态频率峰值出现在凌晨 3:00
左右,为 88.2%,从 6:00 开始迅速下降,在 9:00—
16:00保持较低水平,均值为 5.5%,18:00 以后频率
迅速增加至 70%以上.在 7:30—8:30 和 17:30—
18:00期间,栓皮栎林冠上方大气稳定度变化最为复
杂,稳定、中性和不稳定状态出现的频率均在 30%
左右.统计显示,9:00—15:00 期间,大气不稳定状
态、中性和稳定状态的频率均值分别为 70. 2%、
23 9%和 5.9%,19:00 至次日 6:00,大气不稳定状
态、中性和稳定状态的频率均值分别为 11. 3%、
12 6%和 76.2%,由此可见,昼间大气不稳定,气体
交换活动较强,夜晚大气较稳定,湍流弱,空气流动
性较差.
2 2 栓皮栎人工林冠层 CO2 浓度和 δ13C的变化
由图 2可以看出,测定期间栓皮栎人工林 CO2
浓度均值为 396.08 μmol·mol-1,最高值为 439.79
μmol·mol-1,最低值为 364.52 μmol·mol-1 .CO2 浓
度日变化呈现高低高的趋势,一天内最大值多数出
现在 2:00—6:00,占总测定天数的 70.6%,最小值多
数出现在 12:00—16:00,占总测定天数的 76.5%.
δ13C均值为-14.01‰,最大值为-2.20‰,最小值为
-30.67‰,变幅高达 28.47‰,但 δ13C 值的日变化没
有明显的规律性,一天内最大值出现在 5:00—8:00
的天数占总测定天数的 52. 9%,出现在 10:00—
14:00的天数占 35. 3%,而最小值出现在 12:00—
16:00的天数占总测定天数的 58. 8%,出现在
23:00—次日 3:00的天数占 29.4 %.冠层下部和上
缘的 CO2 浓度和 δ13C值的变化趋势一致,曲线几乎
重合,表明这种变化趋势主要是由于 CO2 和 δ13C 本
底变化造成的 .冠层下部和上缘CO2浓度均值分别
图 2 小时尺度上栓皮栎人工林冠层下部与上缘的 CO2 浓
度和 δ13C的变化
Fig.2 Variations of CO2 concentration and δ13C value between
the canopy bottom and top of Quercus variabilis plantation at
hourly timescale.
Ⅰ: 冠层下部 Bottom of canopy; Ⅱ: 冠层上缘 Top of canopy.
为 396.78 和 395.30 μmol·mol-1,δ13C 均值分别为
-14.35‰和-13.65‰.虽然冠层下部和上缘 CO2 浓
度和 δ13C均值仅相差 1.48 μmol·mol-1和 0.70‰,
但冠层下部 CO2 浓度在曲线的波峰处要明显高于
冠层上缘,而冠层下部 δ13C 值在曲线的波谷处要明
显低于冠层上缘,因而需要结合白天光合作用和夜
晚呼吸作用 2个生理过程在较短的时间尺度(h)上
开展分析,才能更准确地确定造成差异的原因.
2 3 白天栓皮栎人工林冠层 CO2 浓度和 δ13C的差异
由图 3可以看出,在 6:00—18:00,栓皮栎人工
林冠层上缘 CO2 浓度低于冠层下部的时间占到
93 1%,其浓度相差-1.07 ~ 5.72 μmol·mol-1,平均
相差 1.70 μmol·mol-1 .浓度差曲线变化较复杂,规
律性较差,最大浓度差出现时间多在午后,最大浓度
差在午后的日数占全部测定日数的 70.6%,这可能
与栓皮栎午后较高的光合速率有关.在光合作用过
程中,叶片优先吸收12C,引起碳稳定同位素分馏,导
致冠层上缘13C富集,统计显示 92.2%的时间内冠层
上缘 δ13C 值高于冠层下部,δ13C 值相差-3.21‰~
0 65‰,平均差值为-0.81‰.冠层下部与上缘的δ13C
差值曲线变化与 CO2 浓度差值曲线变化趋势相反,
最大 δ13C 差值在午后的天数占全部测定天数的
64 7%.
3732期 孙守家等: 华北低丘山区栓皮栎人工林冠层 CO2浓度和 δ13C变化及其影响因素
图 3 白天(A)和晚上(B)栓皮栎人工林冠层下部与上缘的 CO2 浓度差和 δ13C差
Fig.3 Differences of CO2 concentration and δ13C value between the canopy bottom and top of Quercus variabilis plantation at hourly
timescale in the daytime (A)and nighttime (B).
由图 4可以看出,白天栓皮栎人工林冠层上缘
和下部的 CO2 浓度差值与 δ13C 差值呈显著负相关
(n= 204),表明冠层上部叶片光合作用吸收 CO2 越
多,冠层上缘大气中13C 越富集,冠层下部与上缘的
δ13C差值绝对值越大.
栓皮栎冠层 CO2 浓度和 δ13C 的变化与大气稳
定度和植物生理活动强弱有关.9:00—15:00,大气
不稳定状态占 70.2%,此时冠层内大气由于热力湍
流发展旺盛产生向上的运动,冠层下部的 CO2 随气
流透过冠层向林冠上缘流动.由图 5 可以看出,
92 3%的大气不稳定态的林冠下部 CO2 浓度高于林
冠上缘,97.8%大气不稳定态的林冠下部 δ13C 值低
于林冠上缘.林冠下部 CO2 随气流透过冠层向上流
动时,冠层上部叶片光合作用过程中吸收了部分
CO2,同时冠层下部被遮阴,呼吸释放的 CO2 多于光
合所吸收的 CO2,从而导致冠层下部与上缘的 CO2
浓度和 δ13C值存在差别.
2 4 晚上栓皮栎人工林冠层 CO2 浓度和 δ13C 的
差异
由图 3可以看出,18:00—次日 6:00,栓皮栎人
工林冠层下部 CO2 浓度高于冠层上缘的时间高达
87.3%,浓度相差-2.30~5.92 μmol·mol-1,平均相
图 4 白天(A)和晚上(B)在小时尺度上栓皮栎人工林冠层下部和上缘的 CO2 浓度差与 δ13C差的关系
Fig.4 CO2 concentration difference in relation to δ13C value difference between the canopy bottom and top of Quercus variabilis planta⁃
tion at hourly timescale in the daytime (A) and nighttime (B).
∗∗P<0.01. 下同 The same below.
473 应 用 生 态 学 报 26卷
图 5 稳定 Ç>1和不稳定状态下(Ç<-1)栓皮栎人工林冠层下部与上缘的 CO2 浓度和 δ13C的差异
Fig.5 Differences of CO2 concentration and δ13C value between the canopy bottom and top of Quercus variabilis plantation under stable
(Ç>1) unstable and (Ç<-1) atmospheric conditions.
差 1.24 μmol·mol-1 .浓度差曲线变化的规律性较
差,最大浓度差在凌晨的天数占全部测定天数的
76 5%.叶片呼吸作用释放出的 CO2 中的13C 较少,
导致冠层下部的13C 贫化,冠层下部的 δ13C 值低于
冠层上缘时间高达 91.2%,δ13 C 值相差-2.25‰ ~
0 56‰,平均相差-0.58‰.冠层下部与上缘的 δ13C
差值曲线变化与 CO2 浓度差值曲线变化趋势相反,
最大 δ13C 差值在凌晨的天数占全部测定天数的
58 8%.
夜晚栓皮栎冠层下部和上缘 CO2 浓度差值与
δ13C差值呈显著负相关(n = 204),表明呼吸产生的
CO2 显著地影响了冠层内 δ13C 值变化,呼出的 CO2
越多,冠层下缘空气中12C越贫化(图 4).
19:00—次日 6:00,大气稳定状态高达 76.2%,
由于夜间树冠的辐射作用,使冠层温度低于林冠上
方大气,易形成逆温层,大气层结构稳定.由图 5 可
以看出,84.8 %大气稳定态的林冠下部 CO2 浓度高
于林冠上缘,91.1%大气稳定态的林冠下部 δ13C 低
于林冠上缘.冠层叶片在呼吸过程中释放出 CO2,稳
定大气使得林冠内空气流动性较差,CO2 容易下沉,
在冠层下部聚集,造成 CO2 浓度升高和 δ13 C 值
降低.
2 5 冠层 CO2 和 δ13C差值与环境因子的关联
选择太阳辐射(x1)、温度(x2)、相对湿度(x3)、
大气压(x4)和风速(x5)5 个因子作为影响栓皮栎冠
层 CO2 浓度差值(y1)和 δ13C 差值(y2)变化的环境
因子,采用逐步回归分析法建立白天和晚上的线性
回归模型(表 1).
由表 1可以看出,白天 CO2 浓度差值与环境因
子的线性回归方程达到了显著水平( r = 0.511),太
阳辐射和相对湿度能显著影响冠层下部和上缘
的 CO2浓度,白天上升气流促使冠层内CO2向上流
表 1 栓皮栎人工林冠层 CO2 浓度和 δ13C差值与环境因子间的逐步回归分析
Table 1 Stepwise regression analysis of CO2 concentration and δ13C value differences between the canopy bottom and top of
Quercus variabilis plantation in relation to environmental factors (n=204)
项目
Item
时间
time
逐步回归方程
Stepwise regression equation
r R2 显著性
Sig.
CO2 差值 白天 Daytime 6:00—18:00 y1 = 0.001 x1-0.023 x3+2.55 0.511 0.261 0.000
ΔCO2 晚上 Nighttime 18:00—次日 6:00 -
δ13C差值 白天 Daytime 6:00—18:00 y2 =-0.001x1+0.005x3-0.760 0.525 0.276 0.000
Δ13C 晚上 Nighttime 18:00—次日 6:00 y2 =-0.044x2+0.311 0.476 0.190 0.000
- 未输入任何变量到方程 No variables were involved to the stepwise equation.
5732期 孙守家等: 华北低丘山区栓皮栎人工林冠层 CO2浓度和 δ13C变化及其影响因素
动,经过冠层时部分 CO2 被叶片的光合作用吸收,
从而造成冠层上缘 CO2 浓度下降.但是,夜晚 CO2
浓度差值与环境因子的回归方程不成立.白天和晚
上 δ13C差值与环境因子的线性回归方程均达到显
著水平,2个方程的相关系数分别是 0.525 和 0.476.
白天,太阳辐射和相对湿度影响冠层叶片光合作用,
叶片优先吸收12C,从而导致冠层上缘与下部的 δ13C
存在差异.晚上温度显著影响冠层 δ13C 值变化,温
度越高呼出的 CO2 越多,同时由于夜晚大气稳定性
高,冠层内空气流动性差,呼出的 CO2 容易下沉,导
致冠层下部与上缘的 δ13C存在差异.
3 讨 论
CO2 浓度和 δ13C 值作为植物功能特征的潜在
指标[25],可以用于研究植物生理与生态环境之间的
关系.本研究中,冠层 CO2 浓度呈现高低高的日变化
规律,波峰多出现在 2:00—6:00,波谷多出现在午
后 12:00—16:00,但 δ13C值的日变化没有明显一致
的规律性,主要是 δ13C 值与植物光合作用过程中气
孔的传导和 CO2 的固定有关[26],许多环境因子在不
同程度上影响 CO2 传输和固定,导致 δ13C值变化复
杂.Keeling等[27]和 Pataki等[28]认为,在较长时间尺
度上冠层大气中 CO2 浓度与 δ13C 表现出负相关关
系,许多研究[3,12]也证实了这一结果;但本研究在小
时尺度上的研究结果显示,栓皮栎冠层中 CO2 浓度
与 δ13C相关性不显著,除了受环境因子影响外,栓
皮栎光合和呼吸过程中碳同位素的分馏效应不
同[12]和时间尺度较短也是重要原因.
冠层大气 CO2 的碳同位素组成不是一个常数,
在小时尺度上有复杂的时空变化[8] .叶片光合和呼
吸过程中吸收或释放 CO2 导致冠层 CO2 浓度下降
或上升,同时叶片光合和呼吸过程优先吸收或释
放12CO2,造成冠层大气中13C 富集或贫化[22],导致
冠层内外的 δ13C 值不一样.本研究表明,栓皮栎冠
层下部和上缘 CO2 浓度和 δ13C 值在小时尺度上存
在差异,其中,>87%的时间冠层上缘 CO2 浓度低于
冠层下部,>91%的时间冠层上缘 δ13C 值高于冠层
下部,这与 Broadmeadow等[29]观测到的热带森林不
同冠层 CO2 浓度和 δ13C结果相似,不同高度冠层的
CO2 浓度和 δ13C成层分布现象明显.
CO2 浓度和 δ13C 成层分布现象容易受到大气
稳定状态的影响.白天,热力湍流发展旺盛产生向上
的气流,带动冠层下部的 CO2 向林冠上缘流动,冠
层叶片进行光合作用吸收 CO2,从而导致栓皮栎冠
层上缘 CO2 浓度平均值比冠层下部低 1.70 μmol·
mol-1,冠层上缘的 δ13C 值比冠层下部高约 0.81‰.
Francey等[30]发现,在松树冠层的顶部和底部之间
δ13C 差别达 0.8‰,造成梯度变化的原因可能是树
冠光合、土壤呼吸及落叶分解共同作用的结果.由于
中午光合作用强烈,栓皮栎冠层内外 CO2 和 δ13C 差
值最大时刻多在中午,分别占到总天数的 70.6%和
64.7%.Buchmann等[31]发现,C3和 C4植物中午时刻
的冠层光合汇强度最大,使得冠层 CO2 浓度下降.
Ogée等[15]发现,松树冠层中碳失衡最大差异是在
中午,一般<3‰,主要是由于午后冠层光合汇强度
较大造成的.晚上,大气多处于稳定状态,气流上下
流动性较弱,造成冠层下部和上缘 CO2 和δ13C差异
的原因主要是冠层叶片的呼吸,栓皮栎在呼吸过程
中释放含有较低13C 的 CO2 [32],因此造成冠层下部
CO2 浓度比上缘高 1.24 μmol·mol
-1,冠层下部 δ13C
值低于冠层上缘 0.58‰,这与 Zobitz 等[33]研究冷杉
林气体交换时发现的冠层碳失衡值在±2‰之间的
结果相似.Flanagan等[34]将森林生态系统 CO2 的碳
同位素组成与微气象测量相结合进行分析,发现夜
晚呼吸释放 CO2 通常导致冠层 CO2 中碳同位素贫
化.Schaeffer等[9]发现,在夜晚稳定状态下亚高山带
针叶林冠层内外 δ13C 值相差 0.64‰~0 81‰,大气
稳定性的增加加剧了冠层内外 CO2 和 δ13C差异
程度.
除了植物生理活动强弱和大气稳定性外,降水、
温度、光照和土壤湿度等因素都有可能对冠层 CO2
浓度和 δ13C 产生影响[32-33] .逐步回归分析显示,在
小时尺度上太阳辐射、相对湿度和温度是影响栓皮
栎人工林冠层 CO2 浓度和 δ13C 差异的主要环境因
子.对亚高山带针叶林的研究显示,冠层内外 δ13C
差值与饱和水气压(VPD)有显著的相关关系[9],而
VPD是由温度和湿度计算获得的.因此,光照、温度
和湿度是在较短时间尺度上影响栓皮栎冠层大气
中13C值的重要环境因子,其通过增强栓皮栎的碳
循环过程(光合和呼吸作用)来发挥作用[22,35] .但也
有研究表明,在季节尺度上,马占相思人工林冠层
δ13C差异随气温和 VPD的升高而增高,但没有显著
的正相关关系[22] .
4 结 论
栓皮栎人工林冠层 CO2 浓度有相对一致的日
673 应 用 生 态 学 报 26卷
变化,而 δ13C值的日变化规律不明显.白天,冠层叶
片光合作用吸收 CO2 和大气的不稳定状态共同导
致冠层下部 CO2 浓度高于冠层上缘,而 δ13C值低于
冠层上缘.晚上,稳定的大气状态使得冠层叶片呼吸
释放出 CO2 不易流动,导致冠层下部 CO2 浓度高于
冠层上缘而 δ13C 低于冠层上缘.白天太阳辐射和相
对湿度是影响冠层 CO2 浓度和 δ13C 差异的主要环
境因子,晚上温度能显著地影响冠层 δ13C 值变化,
这些环境因子通过增强光合和呼吸作用来实现栓皮
栎对生态环境的调节作用.
参考文献
[1] Kramer PJ. Carbon dioxide concentration, photosynthe⁃
sis, and dry matter production. Bioscience, 1981, 31:
29-33
[2] Canadell JG, Mooney HA, Baldocchi DD, et al. Carbon
metabolism of the terrestrial biosphere: A multi⁃tech⁃
nique approach for improved understanding. Ecosystems,
2000, 3: 115-130
[3] Yakir D, Wang XF. Fluxes of CO2 and water between
terrestrial vegetation and the atmosphere estimated from
isotope measurements. Nature, 1996, 380: 515-517
[4] 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
[5] Sun W (孙 伟), Lin G⁃H (林光辉), Chen S⁃P (陈
世苹), et al. Applications of stable isotope techniques
and Keeling Plot approach to carbon and water exchange
studies of terrestrial ecosystems. Acta Phytoecologica
Sinica (植物生态学报), 2005, 29(5): 851-862 ( in
Chinese)
[6] Yuan G⁃F (袁国富), Zhang N (张 娜), Sun X⁃M
(孙晓敏), et al. Partitioning wheat field evapotranspi⁃
ration using Keeling Plot method and continuous atmos⁃
pheric vapor δ18O data. Chinese Journal of Plant Ecology
(植物生态学报), 2010, 34(2): 170-178 ( in Chi⁃
nese)
[7] Zhang M (张 弥), Wen X⁃F (温学发), Yu G⁃R
(于贵瑞), et al. Effects of CO2 storage flux on carbon
budget of forest ecosystem. Chinese Journal of Applied
Ecology (应用生态学报), 2010, 21(5): 1201-1209
(in Chinese)
[8] Lee XH, 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, GB3013, doi: 10.1029 / 2006GB002871
[9] 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. Agricul⁃
tural and Forest Meteorology, 2008, 148: 592-605
[10] Griffith D, Deutscher N, Caldow C, et al. A Fourier
transform infrared trace gas and isotope analyzer for
atmospheric applications. Atmospheric Measurement
Techniques, 2012, 5: 2481-2498
[11] Ogée J, Brunet Y, Loustau D, et al. MuSICA, a CO2,
water and energy multilayer, multileaf pine forest mod⁃
el: Evaluation from hourly to yearly time scales and sen⁃
sitivity analysis. Global Change Biology, 2003, 9: 697-
717
[12] Lai CT, Ehleringer JR, Schauer AJ, et al. Canopy⁃scale
δ13C of photosynthetic and respiratory CO2 fluxes: Ob⁃
servations in forest biomes across the United States.
Global Change Biology, 2005, 11: 633-643
[13] Atkin OK, Millar H, Turnbull MH. Plant respiration in
a changing world. New Phytologist, 2010, 187: 268-
272
[14] Griffis TG. Tracing the flow of carbon dioxide and water
vapor between the biosphere and atmosphere: A review
of optical isotope techniques and their application. Agri⁃
cultural and Forest Meteorology, 2013, 174 / 175: 85-
109
[15] Ogée J, Peylin P, Ciais P, et al. Partitioning net eco⁃
system carbon exchange into net assimilation and respi⁃
ration with canopy⁃scale isotopic measurements: An
error propagation analysis with 13CO2 and CO18 O data.
Global Biogeochemical Cycles, 2004, 18, GB2019, doi:
10. 1029 / 2003GB002166
[16] Bowling DR, Pataki DE, Richardson JT. Carbon iso⁃
topes in terrestrial ecosystem pools and CO2 fluxes. New
Phytologist, 2008, 178: 24-40
[17] 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
[18] Ponton S, Flanagan LB, Alstad KP, et al. Comparison
of ecosystem water⁃use efficiency among Douglas⁃fir for⁃
est, aspen forest and grassland using eddy covariance
and carbon isotope techniques. Global Change Biology,
2006, 12: 294-310
[19] He C⁃X (何春霞), Meng P (孟 平), Zhang J⁃S (张
劲松), et al. A study on water use of two fruit tree⁃
wheat intercropping systems in the rocky hilly region of
north china with stable carbon isotope technique. Scien⁃
tia Silvae Sinicae (林业科学), 2012, 48(5): 1-7 (in
Chinese)
[20] Randerson JT, Collatz GJ, Fessenden JE, et al. A pos⁃
sible global covariance between terrestrial gross primary
production and 13C discrimination: Consequence for the
atmospheric 13C budget and its response to ENSO.
Global Biogeochemical Cycles, 2002, 16: 83⁃1-83⁃16
[21] Bowling DR, McDowell NG, Bond BJ. 13C content of
ecosystem respiration is linked to precipitation and vapor
pressure deficit. Oecologia, 2002, 131: 113-124
[22] Zhou 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)
[23] Song Q⁃H (宋清海), Zhang Y⁃P (张一平), Yu G⁃R
7732期 孙守家等: 华北低丘山区栓皮栎人工林冠层 CO2浓度和 δ13C变化及其影响因素
(于贵瑞), et al. Characteristics of CO2 exchange at
leaf and canopy levels in tropical seasonal rain forest of
Xishuangbanna. Chinese Journal of Applied Ecology (应
用生态学报), 2008, 19(4): 723-728 (in Chinese)
[24] Diao Y⁃W (刁一伟), Wang A⁃Z (王安志), Jin C⁃J
(金昌杰), et al. Simulation of CO2 exchange between
forest canopy and atmosphere. Chinese Journal of Applied
Ecology (应用生态学报), 2006, 17 ( 12): 2261 -
2265 (in Chinese)
[25] Tan W⁃B (檀文炳), Wang G⁃A (王国安), Han J⁃M
(韩家懋), et al. δ13C and water⁃use efficiency indicated
by δ13C of different plant functional groups on Changbai
Mountains, Northeast China. Chinese Science Bulletin
(科学通报), 2009, 54(13): 1759 - 1764 ( in Chi⁃
nese)
[26] Yin S⁃P (殷树鹏), Zhang C⁃J (张成君), Guo F⁃Q
(郭方琴), et al. Effect of environmental factors on sta⁃
ble carbon isotopic composition of plants and application
in water use efficiency. Journal of Isotopes (同位素),
2008, 21(1): 46-53 (in Chinese)
[27] Keeling CD, Mook WG, Tans PP. Recent trends in the
13C / 12 C ratio of atmospheric carbon⁃dioxide. Nature,
1979, 277: 121-123
[28] Pataki DE, Ehleringer JR, Flanagan LB, et al. The
application and interpretation of Keeling plots in terres⁃
trial carbon cycle research. Global Biogeochemical Cy⁃
cles, 2003, 17: 1022, doi: 10.1029 / 2001GB001850
[29] Broadmeadow MSJ, Griffiths H, Maxwell C, et al. The
carbon isotope ratio of plant organic material reflects
temporal and spatial variations in CO2 within tropical for⁃
est formations in Trinidad. Oecologia, 1992, 89: 435-
441
[30] Francey RJ, Farquhar GD. An explanation of 13C / 12 C
variations in tree rings. Nature, 1982, 297: 28-31
[31] Buchmann N, Kao WY, Ehleringer JR. Carbon dioxide
concentrations within forest canopies: Variation with
time, stand structure, and vegetation type. Global
Change Biology, 1996, 2: 421-432
[32] Barbour MM, Hunt JE, Kodama N, et al. Rapid chan⁃
ges in δ13C of ecosystem⁃respired CO2 after sunset are
consistent with transient 13C enrichment of leaf respired
CO2 . New Phytologist, 2011, 190: 990-1002
[33] Zobitz JM, Burns SP, Ogée J, et al. Partitioning net
ecosystem exchange of CO2: A comparison of a Baye⁃
sian / isotope approach to environmental regression meth⁃
ods. Journal of Geophysical Research, 2007, 112,
G03013, doi: 10.1029 / 2006JG000282
[34] Flanagan L, Cai T, Black TA, et al. Measuring and
modeling ecosystem photosynthesis and the carbon iso⁃
tope composition of ecosystem⁃respired CO2 in three bo⁃
real coniferous forests. Agricultural and Forest Meteorolo⁃
gy, 2012, 153: 165-176
[35] 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
作者简介 孙守家,男,1975年生,博士,副研究员.主要从事
植物生理生态学研究,发表论文 20 余篇. E⁃mail: ssj1011@
163.com
责任编辑 孙 菊
873 应 用 生 态 学 报 26卷