免费文献传递   相关文献

Effect of seasonal high temperature and drought on carbon flux of bamboo forest ecosystem in subtropical region.

季节性高温和干旱对亚热带毛竹林碳通量的影响


采用涡度相关技术对安吉亚热带毛竹林生态系统碳通量进行连续观测,选取2011和2013年月尺度净生态系统生产力(NEP)、生态系统呼吸(Re)和总生态系统生产力(GEP)数据,探讨季节性高温、干旱对毛竹林生态系统碳通量的影响.结果表明: 安吉毛竹林年际间NEP有较大差异;2013年7和8月由于水热不同步而造成的高温干旱使其NEP明显下降,相比于2011年同期分别下降了59.9%和80.0%.对2011和2013年月尺度下NEP、Re和GEP与环境因子进行相关分析发现,Re和GEP与温度因子均呈显著相关(P<0.05),但两者对空气和土壤水分的响应方式和程度有所不同,GEP相比于Re更易受到土壤水分降低的影响,而饱和水汽压差的升高会在一定程度上促进Re、同时抑制GEP,这是造成2013年7和8月安吉毛竹林NEP降低的根本原因.

The carbon flux of subtropical bamboo forest ecosystem was continuously measured using eddy covariance technique in Anji County of Zhejiang Province, China. The monthly net ecosystem productivity (NEP), ecosystem respiration (Re) and gross ecosystem productivity (GEP) data from 2011 to 2013 were selected to analyze the impacts of seasonal high temperature and drought on the carbon flux of bamboo forest ecosystem. The results showed that there were big differences among annual NEP of bamboo forest from 2011 to 2013. Because of the asynchronization of precipitation and heat, the seasonal high temperature and drought in July and August of 2013 caused significant decline in NEP by 59.9% and 80.0% when compared with the same months in 2011. Correlation analysis of the NEP, Re, GEP and environmental factors suggested that the atmosphere temperatures were significantly correlated with Re and GEP in 2011 and 2013 (P<0.05). However, to air and soil moisture, Re and GEP had different responses, that was, GEP was more vulnerable by the decrease of the soil moisture compared with Re. Besides, the raising of saturation vapour pressure promoted the Re modestly but inhibited the GEP, which was supposed to be the main reason for NEP decrease of bamboo forest ecosystem in Anji from July to August in 2013.


全 文 :季节性高温和干旱对亚热带毛竹林碳通量的影响
陈晓峰1  江  洪1,2 ∗  牛晓栋1  张金梦1  刘玉莉1  方成圆1
( 1浙江农林大学浙江省森林生态系统碳循环与固碳减排重点实验室, 杭州 311300; 2南京大学国际地球系统科学研究所,
南京 210093)
摘  要  采用涡度相关技术对安吉亚热带毛竹林生态系统碳通量进行连续观测,选取 2011
和 2013年月尺度净生态系统生产力(NEP)、生态系统呼吸(Re)和总生态系统生产力(GEP)
数据,探讨季节性高温、干旱对毛竹林生态系统碳通量的影响.结果表明: 安吉毛竹林年际间
NEP 有较大差异;2013年 7和 8月由于水热不同步而造成的高温干旱使其 NEP 明显下降,相
比于 2011年同期分别下降了 59.9%和 80.0%.对 2011和 2013 年月尺度下 NEP、Re和 GEP 与
环境因子进行相关分析发现,Re和 GEP 与温度因子均呈显著相关(P<0.05),但两者对空气和
土壤水分的响应方式和程度有所不同,GEP 相比于 Re更易受到土壤水分降低的影响,而饱和
水汽压差的升高会在一定程度上促进 Re、同时抑制 GEP,这是造成 2013年 7和 8月安吉毛竹
林 NEP 降低的根本原因.
关键词  涡度相关; 碳通量; 亚热带毛竹林; 干旱胁迫
Effect of seasonal high temperature and drought on carbon flux of bamboo forest ecosystem
in subtropical region. CHEN Xiao⁃feng1, JIANG Hong1,2∗, NIU Xiao⁃dong1, ZHANG Jin⁃meng1,
LIU Yu⁃li1, FANG Cheng⁃yuan1 ( 1Zhejiang Provincial Key Laboratory of Carbon Cycling in Forest
Ecosystems and Carbon Sequestration, Zhejiang Agriculture and Forestry University, Hangzhou
311300, China; 2International Institute for Earth System Science, Nanjing University, Nanjing
210093, China) .
Abstract: The carbon flux of subtropical bamboo forest ecosystem was continuously measured using
eddy covariance technique in Anji County of Zhejiang Province, China. The monthly net ecosystem
productivity (NEP), ecosystem respiration (Re ) and gross ecosystem productivity ( GEP) data
from 2011 to 2013 were selected to analyze the impacts of seasonal high temperature and drought on
the carbon flux of bamboo forest ecosystem. The results showed that there were big differences
among annual NEP of bamboo forest from 2011 to 2013. Because of the asynchronization of precipi⁃
tation and heat, the seasonal high temperature and drought in July and August of 2013 caused sig⁃
nificant decline in NEP by 59.9% and 80.0% when compared with the same months in 2011. Corre⁃
lation analysis of the NEP, Re, GEP and environmental factors suggested that the atmosphere tem⁃
peratures were significantly correlated with Re and GEP in 2011 and 2013 (P<0.05). However, to
air and soil moisture, Re and GEP had different responses, that was, GEP was more vulnerable by
the decrease of the soil moisture compared with Re . Besides, the raising of saturation vapour pres⁃
sure promoted the Re modestly but inhibited the GEP, which was supposed to be the main reason for
NEP decrease of bamboo forest ecosystem in Anji from July to August in 2013.
Key words: eddy covariance; carbon flux; subtropical bamboo forests; drought stress.
本文由国家自然科学重大基金项目(61190114)、国家重点基础研究发展计划项目(2011CB302705)和国家自然科学基金项目(41171324)资助
This work was supported by the Major National Natural Science Foundation Project (61190114), National Key Basic Research Project (2011CB302705)
and Project of National Natural Science Foundation of China (41171324).
2015⁃05⁃04 Received, 2015⁃11⁃20 Accepted.
∗通讯作者 Corresponding author. E⁃mail: jianghong_china@ hotmail.com
应 用 生 态 学 报  2016年 2月  第 27卷  第 2期                                            http: / / www.cjae.net
Chinese Journal of Applied Ecology, Feb. 2016, 27(2): 335-344                    DOI: 10.13287 / j.1001-9332.201602.003
    由于人类活动的影响,大气中 CO2等温室气体
浓度持续升高,地球表面的年均温也在一个多世纪
内上升了 0.6 ℃左右,并有逐渐加快的趋势[1-2] .温
度的升高将会在很大程度上改变全球降水格局[3],
导致高纬度地区冬季降水增加、中纬度内陆地区夏
季降雨减少、干旱和涝灾等极端自然灾害发生的频
率也会相应增加[4] .这一系列气候变化将会直接或
间接影响陆地生态系统的生产力、生物地球化学循
环与水循环[5-6] .森林作为最大的陆地生态系统在深
受气候变化影响的同时,对吸收 CO2等温室气体、减
缓全球气候变化等也具有非常重要的意义[7] .第八
次全国森林资源清查(2008—2014年)结果显示,我
国森林总面积 2.08×108 hm2,森林覆盖率由 2008 年
的 20.4%上升到 21. 6%,森林植被总碳储量也从
2008年的 78.11×108 t 上升到 2014 年的 84.27×108
t,总体表现出巨大的碳汇潜力.
毛竹(Phyllostachys edulis)是我国重要的森林资
源,广泛分布于我国亚热带地区.相比于其他树木,
毛竹具有生长速度快、生物量积累大以及生产周期
短等特点,具有重要生态价值[8] .有研究表明,毛竹
林的碳汇能力明显大于南亚热带常绿阔叶林、中亚
热带针阔混交林与暖温带针阔混交林等森林类
型[9-11] .同时,毛竹对于生长环境的要求也颇为苛
刻,喜温暖湿润气候,对于水热变化十分敏感.我国
南方中亚热带地区由于受到亚热带大陆性季风气候
的影响,年水热资源较为丰富,但在夏季该地区的低
层大气环流受太平洋副热带高压控制,易出现晴热
干燥天气,造成季节性干旱或秋旱.近年来,浙江省
高温干旱的天气有逐年上升趋势,这必然会对毛竹
林生态系统造成很大影响.目前,我国关于毛竹叶片
和土壤层面对气候变化的响应研究已取得一定进
展,如徐升华等[12]通过人工截雨的方法创造极端干
旱环境,得出毛竹林土壤含水量与土壤呼吸关系不
紧密的结论.但毛竹林在生态系统层面上对气候变
化响应的相关研究还十分有限.
本研究以安吉毛竹林生态系统为研究对象,利
用涡度相关技术对其进行长期连续观测,在孙成
等[9]研究基础上,对 2011和 2013年(两年同为毛竹
林生理大年)的观测数据进行分析,初步探讨季节
性高温干旱对毛竹林生态系统碳吸收的影响.本研
究结果可以与江西千烟洲和大岗山观测站、广东鼎
湖山观测站以及湖南会同观测站等亚热带地区 CO2
通量观测站观测数据相结合,完善我国亚热带不同
森林类型的通量观测,以期为预测未来全球气候变
化对我国森林生态系统碳汇的影响以及评价相关模
型提供科学依据.
1  研究地区与研究方法
1􀆰 1  研究区概况
研究区位于浙江省湖州市安吉县山川乡
(30°28′34.5″ N、119°40′25.7″ E).站点海拔 380 m,
下垫面坡度为 2.5° ~ 14.0°,坡向北偏东 8°,研究区
土壤类型主要为黄壤土.观测塔周围植被类型以毛
竹为主,采用粗放式经营,林下只有极少量灌木与草
本植物,研究区竹林面积 2155 hm2,其中,毛竹林占
总面积的 78.6%.毛竹平均高度 15.8 m,胸径为 12 ~
18 cm.区内毛竹有大小年之分,2011与 2013年同为
大年,都表现为 3—5 月为发笋期,5 月底毛竹开始
抽叶,6 月大量展开,7—9 月为主要生长季.该地区
年均气温 16.6 ℃,年降水量 761 ~ 1780 mm,降雨主
要集中在夏季,6月至 7 月初为梅雨季节,7 月中下
旬至 8月受副热带高压影响,易出现季节性干旱,年
均相对湿度约在 70%以上,属典型的亚热带季风
气候.
1􀆰 2  观测方法
观测区内气象观测塔高 40 m,开路涡度相关系
统探头安装在观测塔 38 m 高度上(约 3 倍冠层高
度),包括三维超声风速仪(CAST3, Campbell Inc.,
USA)和开路 CO2 / H2O分析仪(Li⁃7500, LiCor Inc.,
USA),原始采样频率为 10 Hz,通过数据采集器
(CR1000,Campbell Inc. ,USA)存储数据,在线计算
并存储 30 min CO2通量(Fc)、显热通量(HS)和潜热
通量(LE)等相关物理量.观测塔配备 7 层气象梯度
观测系统,利用数据采集器 ( CR1000, Campbell
Inc.,USA)每 30 min 自动记录空气温度与湿度、风
速(1、7、11、17、23、30、38 m)、土壤温度、土壤含水
量(5、50、100 cm)、光合有效辐射和净辐射(38 m)
等气象信息(本研究中空气温度与饱和水汽压差选
用 17 m 高冠层气象数据).仪器设备详情参考孙成
等[9]和刘玉莉等[13]相关研究报道.
1􀆰 3  数据处理
1􀆰 3􀆰 1 CO2通量  采用涡度相关技术(EC)测定植被
与大气间 CO2通量.利用三维超声风速仪测定三维
风速,利用红外线 CO2 / H2O气体分析仪测定 CO2密
度脉动,将所有 10 Hz 的原始数据利用数据采集器
记录并储存.采用 EdiRe软件进行协方差运算,并进
行二次坐标轴旋转、WPL 修正等,最终将 10 Hz 的
原始数据转化为 30 min数据[14] .
633 应  用  生  态  学  报                                      27卷
植被与大气间净生态系统 CO2交换量(net eco⁃
system exchange, NEE,mg CO2·m
-2·s-1)的计算公
式如下:
NEE=Fc+Fs (1)
Fs =(△c / △t)·h (2)
式中:Fc为 CO2湍流通量;Fs为冠层内的碳储存通
量[15];h为观测高度,本文取 38 m;△c 为高度 h 处
前后 两 次 相 邻 时 间 测 定 的 CO2 浓 度 差 值
(mg CO2·m
-3);△t为时间间隔,本文取 1800 s.当
CO2从生态系统排放到大气时,定义 NEE符号为正,
当 CO2从大气进入到生态系统时,定义 NEE为负.净
生态系统生产力(net ecosystem productivity, NEP)
符号的定义与 NEE 正好相反,但绝对值相等,即
NEE=-NEP.
1􀆰 3􀆰 2数据的剔除与插补   在实际通量观测过程
中,由于受到仪器响应、天气状况等问题的影响,会
出现一些异常观测值,因此,需要对数据进行有条件
的剔除.主要包括以下几点:1) 降雨同期的观测数
据;2) 低湍流信号下的数据(U∗ <0.2 m·s-1);3)
观测数据值超过仪器的测量量程或合理范围,即
CO2的浓度超出 500~800 mg CO2·m
-3,CO2通量超
出-2.0~2.0 mg CO2·m
-2·s-1,水汽浓度超出 0~40
g·m-3;4)参与统计的有效样本数不足 15000;5)异
常突出的数据,即某一数值与连续 5 点数值平均值
之差的绝对值大于 5点方差的 2.5倍.按以上方法剔
除后的有效数据占总数据的 61.4%.目前能量平衡
是常用来评价数据质量的方法[16],经过计算本站点
能量闭合度为 85%,根据国际相关研究结果,可以
认定本站点通量数据是可靠的.
由于仪器故障、系统标定等影响,涡度相关数据
存在缺失或异常值,因此需要进行数据插补.本研究
采用查表法[17]对缺失数据进行插补.将温度与光合
有效辐射作为主要环境因子,以相邻两月为一个时
间段建立查找表.在一个时间段内,将光合有效辐射
作为主分隔因子,以 100 μmol·m-2·s-1为间隔,在
每个间隔内以温度作为次分隔因子,以 2 ℃为间隔
划分为若干级别,然后分别计算每个间隔内有效
NEE平均值.不同月份的缺失值再根据缺失时的气
象条件在相应查找表中查找相应有效 NEE来插补.
1􀆰 3􀆰 3生态系统呼吸和生态系统总生产力  利用涡
度相关法仅能直接获取 NEE,为了解生态系统呼吸
(Re)和生态系统总生产力(GEP)的变化,需要将
NEE进行拆分.由于夜间生态系统完全处于 CO2排
放状态,所以夜间生态系统呼吸值(Re nig)等于夜间
NEE值,而白天生态系统呼吸(Re day)则是依据 Re nig
与土壤温度拟合的函数关系外延的方法来估算[18] .
日尺度上的 Re可定义为:
Re = Re day+ Re nig (3)
总生态系统交换量( gross ecosystem exchange,
GEE)可定义为:
GEE=NEE-Re (4)
GEE与 GEP 的绝对值相等但符号正好相反,即
GEP =-GEE.
2  结果与分析
2􀆰 1  安吉毛竹林生态系统气象因子变化
研究区 2011 和 2013 年降雨量分别为 1543. 1
和 1538.5 mm,相差不大,但均高于历史年均降雨量
(1270 mm). 其中,2011年降雨集中于夏季,占全年
降雨量的 68.8%;2013 年降雨较为分散,主要集中
在 2、6、8 和 10 月,在毛竹林生长季的 7、8 月,降雨
量为 214.8 mm,明显低于 2011年的 609.2 mm,与历
史同期降雨量相比下降了 37.8%(图 1).
受太阳高度角和云层变化的影响,两年月平均
光合有效辐射 ( photosynthetically active radiation,
PAR)呈双峰曲线变化,6月出现波谷.由于 6月进入
梅雨季,该月份 PAR 均出现较大程度降低,进入盛
夏后 PAR迅速升高.其中,2013 年 7、8 月由于晴热
干燥天气较多,导致月均 PAR 值分别达到 521.5 和
406. 0 μmol · m-2 · s-1,比 2011 年同期分别高
42􀆰 8%和 29.2%.受太阳辐射的影响,2013 年 7、8 月
的平均温度分别比 2011 年同期高 2. 7 和 3. 0 ℃ .
2013年 7、8 月,日最高气温达到 37 ℃的天数为 39
d,日最高气温达 40 ℃的天数为 9 d.由于温度升高、
降雨量减少,导致研究区 2013 年 7、8 月 5 cm 深土
壤含水量明显低于 2011 年同期值,少数观测日 5
cm深土壤含水量甚至低于 0.15 m3·m-3,50 cm 深
土壤含水量的下降幅度小于 5 cm土壤含水量.饱和
水汽压差(vapor pressure deficit,VPD)是由气温和空
气相对湿度计算得到[19],受水热条件变化的影响,
2013年 7、8月的月平均 VPD高达 1.79和 1.44 kPa.
从 2011 和 2013 年全年数据来看,安吉毛竹林
生态系统水热资源丰富,2011 年雨热同期性较好,
但 2013年夏季雨热不同期,造成了一定程度的土壤
和空气干旱.
2􀆰 2  生态系统碳交换动态
2􀆰 2􀆰 1 总生态系统生产力   2011和2013年安吉毛
7332期                        陈晓峰等: 季节性高温和干旱对亚热带毛竹林碳通量的影响         
图 1  安吉毛竹林生态系统 2011和 2013年气象要素的变化
Fig. 1   Variation of meteorological elements in Anji bamboo
forest ecosystem in 2011 and 2013.
P: 降雨量 Precipitation; SW5: 土壤 5 cm含水量 Soil water content at
the depth of 5 cm; SW50: 土壤 50 cm 含水量 Soil water content at the
depth of 50 cm; Ta: 17 m高空气温度 Air temperature at the height of
17 m; VPD: 17 m 高饱和水汽压差 Vapour pressure deficit at the
height of 17 m; PAR:光合有效辐射 Photosynthetically active radiation.
竹林森林生态系统的总生态系统生产力(GEP)分
别为 1568.6和 1461.4 g C·m-2·a-1 .受生长习性与
环境作用的影响[20-21],月 GEP 总值有较大的变化
波动(图 2).2011 年 1—3月,随着温度和 PAR的升
高,光合作用增强,GEP 逐月增大;4—5月为毛竹林
出笋期,6月毛竹笋爆发式生长完成,但由于进入梅
雨季节,受降雨影响,PAR降低(图 1),导致 GEP 变
化不明显;6月底出梅,PAR 迅速升高,同时伴随着
新竹叶片全部展开,叶面积逐渐稳定,毛竹林生态系
统代谢旺盛,导致 7 月 GEP 迅速升高,达 191􀆰 8
g C·m-2,成为全年最高值;之后,随着 PAR 和温度
的逐渐降低,光合速率降低,GEP 也逐月降低,到 12
月达到全年最低值,为 76.2 g C·m-2 .与 2011 年相
比,2013年 2—3月 GEP 的降低可能与此时段阴雨
天增多导致PAR降低有关,而7—8月GEP的降低
图 2  安吉毛竹林总生态系统生产力(GEP)、生态系统呼吸
(Re)和净生态系统生产力(NEP)的月变化动态
Fig.2  Dynamics of monthly GEP, Re and NEP in Anji bamboo
forest ecosystem.
则应该是期间的高温干旱不利于毛竹光合作用所
致;8月之后,尽管雨热条件也有一定程度差异,但
与 2011年相同月份差异不大.
2􀆰 2􀆰 2生态系统呼吸  2011 和 2013 年生态系统呼
吸(Re)的月变化呈单峰型曲线(图 2).从 2 月开始
随着温度升高,毛竹生长加快,Re逐渐升高,最高值出
现在夏季温度最高的 7月,为 92􀆰 5 g C·m-2,8月以
后,由于气温降低,毛竹生理活动减弱,Re逐渐降低,
12月达到全年最低值,仅 61􀆰 2 g C·m-2 .两年 Re的总
体变化趋势相似,变化幅度超过 5%的月份有 1、3、
4、7、8和 11月.2013年 Re相对较大的月份的空气温
度往往高于 2011 年相对应月份,如 2013 年 3 和 7
月的 Re分别比 2011年对应月份高 14.6%和 6􀆰 4%,
2013年 3 和 7 月的月平均气温分别比 2011 年对应
月份高 2.85 和 2.73 ℃ .但 2013 年 8 月平均空气温
度比 2011年同期高 3.02 ℃,其 Re却比 2011年同期
低 6.5%,这可能与该月长时间的干旱胁迫在一定程
度上抑制了毛竹林的呼吸有关.
833 应  用  生  态  学  报                                      27卷
2􀆰 2􀆰 3净生态系统生产力  2011 和 2013 年安吉毛
竹林森林生态系统的净生态系统生产力(NEP)分
别为 638.8和 532.0 g C·m-2·a-1,各月 NEP 均为
正值且变化幅度较大(图 2). 2011 年 1—3月,随着
温度以及 PAR 的升高,NEP 逐月增大;4、5 月是毛
竹林的出笋高峰期,该阶段生态系统呼吸增加,碳排
放增强,NEP 有所下降;6 月受梅雨季的影响,其
NEP 并未上升,6月底出梅,温度和 PAR 迅速升高,
导致 7月 NEP 迅速升高,达到 99.3 g C·m-2,为全
年最大值;之后;随着 PAR和温度的逐渐降低,NEP
也逐月降低.对比 2011 年,2013 年 NEP 在 3、7 和 8
月出现较大差异,其他月份差异不明显.2013年 3月
NEP 降低值的 47.3%为 GEP 降低造成,52.7%为 Re
升高造成,7、8月降低的 NEP 则几乎全由 GEP 的降
低决定.总体来说,2013 年 NEP 的变化更多取决于
GEP 的影响,而 Re的影响不及前者.
由于 NEP 是 GEP 与 Re 的差值,即 GEP 与 Re
共同决定了 NEP 的大小,本文利用月尺度 Re和
GEP 的耦合关系结合 NEP 值进行对比分析(图 3).
结果表明,2011和 2013年 Re与 GEP 呈极显著正相
关关系(P<0.01),但是 2011 年 3 月、2013 年 7 和 8
月 Re与 GEP 关系出现一定程度偏离现象,这与同期
NEP 出现差异的现象相吻合.表明 Re和 GEP 的耦合
关系很大程度上决定了毛竹林生态系统碳吸收特
征,季节性干旱及其他气候变化将会影响 Re同 GEP
的耦合,从而影响到 NEP.
2􀆰 3  环境要素对安吉毛竹林生态系统碳交换的影响
植物的光合作用主要受到温度、水分和太阳辐
图 3  安吉毛竹林生态系统 2011 和 2013 年生态系统呼吸
(Re)与总生态系统生产力(GEP)的关系
Fig.3  Relationships between Re and GEP in Anji bamboo forest
ecosystem in 2011 and 2013.
实线代表 2011 和 2013 年 Re与 GEP 线性拟合曲线,但不包括 2011
年 3月和 2013年 7和 8月数据 The full line stood for Re and GEP fit⁃
ting curve of 2011 and 2013, and it did not include data in March of
2011, July and August of 2013.
射等环境因子的影响,呼吸作用又受控于土壤水分
和温度.因此,空气和土壤温度、光合有效辐射、饱和
水汽压差和土壤含水量都将影响到安吉毛竹林的碳
交换[20-23] .本研究对毛竹林生态系统 2011 和 2013
年月尺度上的 NEP、GEP 和 Re与各环境因子进行回
归关系分析.
2􀆰 3􀆰 1环境要素对净生态系统生产力的影响  由表
1可以看出,2011 年安吉毛竹林月尺度 NEP 与空气
温度和光合有效辐射呈极显著正相关(P<0.01),与
土壤温度呈显著正相关(P<0.05) .表明空气温度和
光合有效辐射是控制毛竹林生态系统月尺度 NEP
的主要环境要素,但 2013 年月尺度 NEP 与光合有
效辐射和空气温度的相关性不显著(P>0.05).从图
4可以看出,2011 年 NEP 与空气温度和光合有效辐
射呈极显著正相关(P<0.01),NEP 随温度和光合有
效辐射的升高而升高,但在 2013 年 7、8 月高温、强
辐射条件下,NEP 却出现明显降低,表明 2013 年 7、
8月温度与 PAR并不是制约 NEP 的主要因子,这也
是导致 2013 年 NEP 与温度和 PAR 相关系数降低
的主要原因.
    2011年,安吉毛竹林月尺度 NEP 与饱和水汽
压差(VPD)呈显著相关关系(P<0.05),说明 VPD
也是影响生态系统净碳交换的重要因子[24] .随着
VPD的增大,NEP 逐渐升高,当 VPD 大于 0.9 kPa
时,由于植物叶片水势增加,导致气孔关闭,从而降
低了光合作用速率[25],NEP 开始逐渐下降.2013年,
表 1  安吉毛竹林生态系统月净生态系统生产力(NEP)、生
态系统呼吸(Re)、总生态系统生产力(GEP)与环境因子的
相关系数
Table 1   Correlation coefficients of monthly NEP, Re,
GEP with environmental factors in Anji bamboo forest eco⁃
system
指标
Index
年份
Year
PAR Ta Ts VPD SW50 SW5
NEP 2011 0.741∗∗ 0.716∗∗ 0.695∗ 0.747∗ 0.569 0.236
2013 0.205 0.255 0.268 0.567 0.562 0.309
GEP 2011 0.764∗∗ 0.856∗∗ 0.843∗∗ 0.824∗∗ 0.646∗ 0.306
2013 0.502 0.606∗ 0.609∗ 0.539 0.620∗ 0.215
Re 2011 - 0.941∗∗ 0.949∗∗ 0.789∗ 0.642∗ 0.386
2013 - 0.918∗∗ 0.907∗∗ 0.731∗ 0.428 0.081
PAR: 光合有效辐射 Photosynthetically active radiation; Ta: 17 m高空
气温度 Air temperature at the height of 17 m; Ts: 5 cm深处土壤温度
Soil temperature at the depth of 5 cm; VPD:17 m 高饱和水汽压差
Vapour pressure deficit at the height of 17 m; SW5: 土壤 5 cm含水量
Soil water content at the depth of 5 cm; SW50: 土壤 50 cm含水量 Soil
water content at the depth of 50 cm.∗ P<0.05; ∗∗P<0.01. NEP、GEP
和 Re与 VPD呈非线性相关,与其他环境因子呈直线相关 It represen⁃
ted nonlinear correlation between VPD and NEP, GEP, Re, but linear
relationships with other factors.
9332期                        陈晓峰等: 季节性高温和干旱对亚热带毛竹林碳通量的影响         
图 4  安吉毛竹林生态系统 2011和 2013年月净生态系统生产力(NEP)与空气温度(Ta)和光合有效辐射(PAR)的关系
Fig.4  Relationships of monthly NEP with air temperature and photosynthetically active radiation in Anji bamboo forest ecosystem in
2011 and 2013.
◆号代表 2013年 7和 8月 NEP 值,图中实线代表不包括◆号数据的拟合曲线 The symbol “◆” stood for the NEP of July and August 2013. The
full line stood for the fitted curve and it did not include the data of the symbol “◆”.
NEP 与 VPD同样呈二次曲线关系(图 5),虽然其相
关关系不显著(P>0.05),但 VPD和 SW50与 NEP 的
相关系数在各环境因子中最高,分别为 0. 567 和
0􀆰 562,两者相差不大,说明在 2013年 7 和 8 月的季
节性高温干旱情况下,VPD 和 SW50是制约 NEP 的
重要环境因子.
2􀆰 3􀆰 2环境要素对总生态系统生产力的影响  2011
年,安吉毛竹林 GEP 与温度、光合有效辐射和 VPD
均呈极显著相关关系(P<0.01)(表 1),与 SW50呈显
著相关关系(P<0.05),但与表层土壤含水量的相关
图 5  安吉毛竹林月净生态系统生产力(NEP)与饱和水汽
压差(VPD)的关系
Fig.5   Relationships between the monthly NEP and VPD in
Anji bamboo forest ecosystem.
性不显著(P>0.05).这表明,温度和光合有效辐射
是安吉毛竹林生长的主要控制因子,而 VPD 和
SW50也是影响毛竹林 GEP 的重要环境因子. 2013
年,GEP 与 SW50的相关性最高,达到显著水平(P
<0.05),说明在 2013年 7、8 月 SW50是制约 GEP 的
主要因子,同时也导致该年度 PAR、空气和土壤温
度与 GEP 的相关性明显低于 2011 年 (图 6).在
2013年温度最高的 7、8 月,GEP 明显降低,严重偏
离温度与 GEP 的趋势线,但 SW50与 GEP 的趋势线
图中的偏离程度不大.此外,饱和水汽压差是控制气
孔导度和冠层导度的重要因子. 2011 年,当 VPD>
0􀆰 9 kPa(2013 年结果较高,为 1.1 kPa)时,毛竹林
叶片的光饱和点将会减小,从而降低毛竹林的固碳
能力. 2013年 7、8月,在较高的 VPD条件下 GEP 值
降低幅度较小,表明 VPD 虽然对 GEP 有一定程度
影响,但其并不是制约 GEP 的主要因子.
2􀆰 3􀆰 3环境要素对生态系统呼吸的影响  空气和土
壤温度与 Re的相关性达到极显著水平(P<0.01),
VPD和 SW50与 Re呈显著正相关(P<0.05),5 cm 深
土壤含水量与 Re的相关性不显著(P>0.05),表明温
度是控制 Re的主要环境因子,VPD 和 SW50也是影
响 Re的重要环境因子.2013年的季节性高温使得温
度和VPD与Re的相关性与2011年差异不明显,但
043 应  用  生  态  学  报                                      27卷
图 6  安吉毛竹林生态系统 2011和 2013年月总生态系统生
产力(GEP)与空气温度(Ta)、50 cm土壤含水量(SW50)和饱
和水汽压差(VPD)的关系
Fig.6  Relationships of the monthly GEP with air temperature,
soil moisture content at 50 cm and vapor pressure deficit in Anji
bamboo forest ecosystem in 2011 and 2013.
实线代表 2011年线性拟合曲线,虚线代表 2013年线性拟合曲线,但
不包括∗号数据数据 The full line stood for the fitted curves of 2011.
The dotted line stood for the fitted curves of 2013 and it did not include
the data of the symbol “∗”.
SW50与 Re的相关性则明显降低.由图 7 可以看出,
SW50与 Re呈显著正相关(P<0.05),即土壤含水量的
升高将会促进生态系统呼吸,但 2013年 7、8 月土壤
含水量的降低,并未导致 Re值明显降低,说明 2013
年 7、8月较低的土壤含水量并没有影响毛竹林生态
系统呼吸,温度仍是控制 Re的主要环境因子,同时
VPD的升高也在一定程度上促进了生态系统呼吸.
2􀆰 4  水热耦合对于碳吸收的影响
NEP 及其两个组分(GEP 和 Re)都会受到温度
和水分的影响,水分和热量的耦合关系很大程度上
决定了 NEP 值的大小.研究期间,毛竹林 NEP 与
SW50的相关性达到极显著水平(P<0.01),与 VPD
和降雨量的相关性不显著(P>0.05) .因此,本研究
图 7  安吉毛竹林生态系统 2011和 2013年月生态系统呼吸
(Re)与 50 cm土壤含水量(SW50)的关系
Fig.7  Relationships between the monthly Re and soil moisture
content at 50 cm in Anji bamboo forest ecosystem in 2011 and
2013.
∗号代表 2013年 7和 8 月 NEP 值,图中实线代表不包括∗号数据
的拟合曲线 The symbol “∗” stood for the NEP of July and August
2013. The full line stood for the fitted curves and it did not include the
data of the symbol “∗”.
图 8  安吉毛竹林生态系统 2011 和 2013 年空气温度与 50
cm土壤含水量(SW50)的关系
Fig.8  Relationships between air temperature and soil moisture
content at 50 cm in Anji bamboo forest ecosystem in 2011 and
2013.
实线代表 2011和 2013年空气温度与 SW50线性拟合曲线,但不包括
2011年 4和 5月和 2013 年 2、3、7 和 8 月数据 The full line stood for
air temperature and SW50 fitting curves of 2011 and 2013, and it did not
include April and May of 2011 data and February, March, July and Au⁃
gust of 2013 data.
选取 SW50代表水分条件与空气温度进行耦合分析
(图 8),结果表明,2011 和 2013 年温度和土壤含水
量总体呈极显著正相关关系(P<0.01).说明空气温
1432期                        陈晓峰等: 季节性高温和干旱对亚热带毛竹林碳通量的影响         
度和 SW50的耦合关系总体较好,除了少数月份
(2011年的 4、5 月,2013 年的 2、3 月以及出现季节
性高温干旱的 7和 8 月)温度和土壤含水量关系表
现出不同程度的偏离. 2011 年 4、5 月的 NEP 为
2013年相应月份的 76.6%和 85.8%,2013 年 2、3 月
的 NEP 分别为 2011年同期的 75.4%和 73.5%,2013
年 7 和 8 月的 NEP 分别为 2011 年对应月份的
40􀆰 1%和 20.0%(图 2).可见,温度和土壤含水量的
耦合关系很大程度上决定了安吉毛竹林生态系统碳
吸收的特征,当温度和水分同步性较好时,NEP 随
着两者升高而升高,反之,将导致 NEP 降低,且其偏
离程度越高,NEP 降低幅度越大.
3  讨    论
温度和水分是影响生态系统呼吸的主要环境因
素.目前,生态系统呼吸随温度呈指数增加已得到学
术界广泛认同[10,26-27],而水分对生态系统呼吸作用
的影响机理较为复杂.土壤水分通过调控植物、微生
物的生理活动以及土壤通透性等方式来影响生态系
统的呼吸.在土壤水分较低时,土壤微生物活性降
低,凋落物分解减缓,此外,由于受气孔和非气孔等
因素的影响,植物光合作用降低,也会致使向根系输
送的有机物减少,从而降低根系呼吸的底物供应,最
终造成生态系统呼吸降低[28] .本研究区 2013 年 7、8
月出现一定程度的土壤干旱,但主要出现在土壤表
层,深层土壤含水量降低程度不高,不足以抑制植物
体的呼吸作用.另外,土壤含水量降低的同时还伴随
着土壤温度的升高,在一定程度上会促进土壤呼
吸[12],多方面因素综合作用导致生态系统呼吸与土
壤含水量相关性不高.空气饱和水汽压差的升高往
往伴随着温度的升高[29],因此,在一定范围内饱和
水汽压差与生态系统呼吸呈正相关.
光照和温度是植物光合作用的重要影响因素,
一定范围内,两者升高会引起光合速率升高,但超过
最适点后,光合速率则随着两者的升高而减弱.孙成
等[9]和施建敏等[30]研究表明,夏季毛竹林具有明显
的“光合午休”现象.安吉 2013年 7、8月气温和光合
有效辐射明显高出 2011年同期值,从而延长了毛竹
林的“午休”时长,这也是造成 NEP 降低的重要原
因.水分以土壤水分和大气水分的形式来影响生态
系统,两者是森林生态系统总生产力的关键限制因
子.土壤和大气水分的降低所造成的水分胁迫都会
引起植物叶片气孔关闭、光合作用和蒸腾作用降低,
从而使植物生长速度减缓,最终导致总生产力降
低[24] .本研究选取土壤含水量和饱和水汽压差分别
表征土壤水分和大气水分,在饱和水汽压差较高条
件下(月平均 VPD>1. 0 kPa),毛竹林光饱和点降
低,从而抑制总生态系统生产力升高. 2013 年 VPD
与安吉毛竹林生态系统 GEP 的相关系数低于土壤
50 cm含水量与 GEP 的相关系数(表 1),表明土壤
水分是安吉毛竹林 2013 年 7 和 8 月 GEP 下降的主
要影响因素.
安吉毛竹林生态系统 Re和 GEP 都会受到高温
干旱的影响,但两者响应方式和程度的不同是造成
研究区净生态系统生产力降低的根本原因.Van der
Molen等[31]和孙晓敏等[32]研究表明,生态系统碳源
与汇的强度变化取决于 Re和 GEP 对于干旱胁迫的
不同响应方式和程度.季节性高温干旱归根结底是
该地区的水热不同步.但水分过多也会对 NEP 造成
影响,降雨增多,光照不足,也会导致 Re和 GEP 紊
乱.如 2013年 2月研究区降雨量为 2013年同期的 7
倍,PAR仅为 2011年同期的 61.9%,直接导致 GEP
比 2011年同期下降 8.5%,NEP 更是降低了 24.6%.
高温干旱并不是影响森林生态系统碳吸收的唯一因
素.为更准确地了解森林生态系统碳交换过程,还需
要考虑到洪涝、冻害、雪灾、寒流等各种极端环境条
件,并对其发生时期、持续时间和强度等要素进行综
合分析.
4  结    论
本研究利用涡度相关技术对安吉毛竹林生态系
统进行连续观测,对比 2011 和 2013 年碳通量特征
发现,安吉毛竹林净生态系统生产力年际间变化明
显,由于水热不同步造成的季节性高温干旱导致
2013年 7、8月碳吸收能力出现明显下降.对月尺度
NEP、GEP 和 Re与相关环境因子进行相关分析发
现,在高温干旱条件下,SW50和 VPD是制约 NEP 的
主要环境因子,但两者间差异不明显;Re主要受控于
温度和 VPD,VPD的升高一定程度上促进了生态系
统呼吸,土壤含水量对 Re的影响不大;SW50是安吉
毛竹林 2013 年 7、8 月 GEP 降低的主要驱动因子,
VPD虽然也能在一定程度上抑制 GEP,但其影响程
度不及 SW50显著. Re和 GEP 对高温干旱响应方式
与程度的不同是毛竹林碳吸收变化的根本原因.以
SW50表征水分条件与空气温度进行耦合分析,发现
水热不同步会不同程度造成 NEP 的降低,而季节性
高温干旱仅仅是水热不同步的一种主要形式.
243 应  用  生  态  学  报                                      27卷
致谢  实验站点建设工作得到浙江省安吉县林业局的大力
支持,浙江农林大学国际生态中心杨爽、陈云飞、孙成等为本
试验的顺利完成做出巨大贡献,特此感谢!
参考文献
[1]  NOAA. State of the Climate Report. New York: National
Climatic Data Center, 2012
[2]  IPCC. Intergovernmental Panel on Climate Change, Cli⁃
mate change: Synthesis report. Fourth Assessment Re⁃
port of the Intergovernmental Panel on Climate Change.
Cambridge: Cambridge University Press, 2007
[3]  Karl TR, Trenberth KE. Modern global climate change.
Science, 2003, 302: 1719-1723
[4]  Meehl GA, Tebaldi C. More intense, more frequent,
and longer lasting heat waves in the 21st century. Sci⁃
ence, 2004, 305: 994-997
[5]  Ciais P, Reichstein M, Viovy N, et al. Europe⁃wide re⁃
duction in primary productivity caused by the heat and
drought in 2003. Nature, 2005, 437: 529-533
[6]  Zhou J, Zhang ZQ, Sun G. Response of ecosystem car⁃
bon fluxes to drought events in a poplar plantation in
Northern China. Forest Ecology and Management, 2013,
300: 33-42
[7]  Bonan GB. Forests and climate change: Forcings, feed⁃
backs, and the climate benefits of forests. Science,
2008, 320: 1444-1449
[8]  Zhou G⁃M (周国模). Carbon Storage, Fixation and
Distribution in Phyllostachys edulis Stands Ecosystem.
PhD Thesis. Hangzhou: Zhejiang University, 2006 ( in
Chinese)
[9]  Sun C (孙  成), Jiang H (江  洪), Zhou G⁃M (周
国模), et al. Variation characteristics of CO2 flux in
Phyllostachys edulis forest ecosystem in subtropical re⁃
gion of China. Chinese Journal of Applied Ecology (应用
生态学报), 2013, 24(10): 2717-2724 (in Chinese)
[10]  Liu Y⁃F (刘允芬),Yu G⁃R (于贵瑞), Wen X⁃F (温
学发), et al. The season of CO2 flux variation charac⁃
teristics from subtropical plantation ecological system.
Science in China Series D: Earth Sciences(中国科学 D
辑:地球科学), 2006, 36(suppl.1): 91-102 (in Chi⁃
nese)
[11]  Wu J⁃B (吴家兵), Guan D⁃X (关德新), Zhang M
(张  弥), et al. Carbon budget characteristics of the
broadleaved Korean pine forest in Changbai Mountains.
Journal of Beijing Forestry University (北京林业大学学
报), 2007, 29(1): 1-6 (in Chinese)
[12]  Xu S⁃H (徐升华), Zhou B⁃Z (周本智), Li Q (李 
谦), et al. Effects of extreme drought on soil respiration
in the Phyllostachys pubescens forest. Chinese Agricultur⁃
al Science Bulletin (中国农学通报), 2014, 30(1):
58-62 (in Chinese)
[13]  Liu Y⁃L (刘玉莉), Jiang H (江   洪), Zhou G⁃M
(周国模), et al. Water vapor flux variation characteris⁃
tic and the relationship with its environment factors in
Phyllostachys edulis forest in Anji. Acta Ecologica Sinica
(生态学报), 2014, 34(17): 4900 - 4909 ( in Chi⁃
nese)
[14]  Xu Z⁃W (徐自为), Liu S⁃M (刘绍民), Gong L⁃J (宫
丽娟) , et al. A study on the data processing and quali⁃
ty assessment of the eddy covariance system. Advances in
Earth Science (地球科学进展), 2008, 23(4): 357-
370 (in Chinese)
[15]  Wang C⁃L (王春林), Yu G⁃R (于贵瑞), Zhou G⁃Y
(周国逸), et al. Estimation of CO2 flux in a coniferous
and broad⁃leaved mixed forest in Dinghushan. Science in
China Series D: Earth Sciences (中国科学 D 辑:地球
科学), 2006, 36(suppl. 1): 119-129 (in Chinese)
[16]   Schmid HP, Grimmond CSB, Cropley F, et al. Meas⁃
urements of CO2 and energy fluxes over a mixed hard⁃
wood forest in the mid⁃western United States. Agricultur⁃
al and Forest Meteorology, 2000, 103: 357-374
[17]  Falge E, Baldocchi D, Olson R, et al. Gap filling strat⁃
egies for long term energy flux data sets. Agricultural and
Forest Meteorology, 2001, 107: 71-77
[18]  Lloyd J, Taylor JA. On the temperature dependence of
soil respiration. Functional Ecology, 1994, 8: 315-323
[19]  Yu G⁃R (于贵瑞), Wang Q⁃F (王秋凤), Mi N (米
娜). Ecophysiology of Plant Photosynthesis, Transpira⁃
tion, and Water Use. Beijing: Science Press, 2010 ( in
Chinese)
[20]  Cramer W, Kicklighter DW, Bondeau A, et al. Compa⁃
ring global models of terrestrial net primary productivity
(NPP): Overview and key results. Global Change Bio⁃
logy, 1999, 5: 1-15
[21]  Verstraeten WW, Veroustraete F, Feyen J. On tempera⁃
ture and water limitation of net ecosystem productivity:
Implementation in the C⁃Fix model. Ecological Model⁃
ling, 2006, 199: 4-22
[22]  Zhou L⁃Y (周丽艳), Jia B⁃R (贾丙瑞), Zeng W (曾
伟), et al. Carbon exchange of Chinese boreal forest
during its growth season and related regulation mecha⁃
nisms. Chinese Journal of Applied Ecology (应用生态学
报), 2010, 21(10): 2449-2456 (in Chinese)
[23]   Chimner RA. Soil respiration rates of tropical peatlands
in Micronesia and Hawaii. Wetlands, 2004, 24: 51-56
[24]  Wang H, Saigusa N, Yamamoto S, et al. Net ecosystem
CO2 exchange over a larch forest in Hokkaido, Japan.
Atmospheric Environment, 2004, 38: 7021-7032
[25]  Anthoni PM, Law BE, Unsworth MH. Carbon and water
vapor exchange of an open⁃canopied ponderosa pine eco⁃
system. Agricultural and Forest Meteorology, 1999, 95:
151-168
[26]  Falge E, Baldocchi D, Tenhunen J, et al. Seasonality of
ecosystem respiration and gross primary production as
derived from FLUXNET measurements. Agricultural and
Forest Meteorology, 2002, 113: 53-74
3432期                        陈晓峰等: 季节性高温和干旱对亚热带毛竹林碳通量的影响         
[27]  Yu G⁃R (于贵瑞), Wen X⁃F (温学发), Li Q⁃K (李
庆康), et al. Seasonal pattern and environment response
characteristic of typical forest ecosystem respiration in
Chinese subtropics and temperate zone. Science in China
Series D: Earth Sciences (中国科学 D 辑:地球科
学)), 2004, 34(suppl.Ⅱ): 84-94 (in Chinese)
[28]  Boone RD, Nadelhoffer KJ, Canary JD, et al. Roots
exert a strong influence on the temperature sensitivity of
soil respiration. Nature, 1998, 396: 570-572
[29]  Wang HM, Saigusa N, Zu YG, et al. Carbon fluxes and
their response to environmental variables in a Dahurian
larch forest ecosystem in northeast China. Journal of
Forestry Research, 2008, 19: 1-10
[30]  Shi J⁃M (施建敏), Guo Q⁃R (郭起荣), Yang G⁃Y
(杨光耀), et al. Seasonal photosynthetic responses of
Phyllostachys edulis to environmental factors. Guihaia
(广西植物), 2007, 27(6): 923-928 (in Chinese)
[31]  Van der Molen MK, Dolman AJ, Ciais P, et al. Drought
and ecosystem carbon cycling. Agricultural and Forest
Meteorology, 2011, 151: 765-773
[32]  Sun X⁃M (孙晓敏), Wen X⁃F (温学发), Yu G⁃R
(于贵瑞), et al. Seasonal drought effects on carbon
sequestration of a mid⁃subtropical planted forest of sou⁃
thern China. Science in China Series D: Earth Sciences
(中国科学 D 辑:地球科学), 2006, 36( supp. 1):
103-110 (in Chinese)
作者简介  陈晓峰,男,1990年生,硕士研究生. 主要从事森
林生态系统碳循环研究. E⁃mail: 18358107852@ 163.com
责任编辑  杨  弘
封 面 说 明
图片由中国科学院教育部水土保持与生态环境研究中心博士研究生孙龙于 2014年 7 月 26 日
拍摄于子洲县郊区.子洲县是本团队于 2014 年开展样线调查的七个调查地点之一,七个调查地点
由南向北依次是宜君、富县、延安、子长、子洲、榆林、鄂尔多斯.子洲县位于陕西省北部,隶属于陕西
省榆林市,北接米脂,东连绥德,南同清涧、子长接壤.子洲县属中温带大陆季风气候,四季分明,天
气多变,降水量较小,冬季寒冷干燥,雨热同期.子洲县年均气温 9.2 ℃,全年无霜期 164 d,年平均
降水量 499.1 mm,全年 70%降水主要集中在夏季的 6—9月.子洲县属于典型的黄土丘陵区,黄土丘
陵区严重的水土流失导致土壤有机碳流失和土地退化,给区域农业生态和经济发展带了巨大阻碍.
为有效控制水土流失,黄土高原自 20世纪 60 年代开始恢复植被,并于 1999 年进一步开展退耕还
林工程.经过几十年的努力,黄土丘陵区林草覆盖度显著增加,以往的农耕地或者裸地已变成林草
覆盖的绿地,土地利用类型发生了根本变化,因此黄土高原水土流失现象显著降低,黄河年泥沙输
移量也从以往的 16×108 t锐减到 1.5×108 t.植被恢复被认为是提高土壤固碳能力的有效方式.明确
黄土丘陵区植被恢复后土壤有机碳分布特征,对于评价植被恢复固碳和农业可持续发展具有重要
意义.
陈晓峰, 江洪, 牛晓栋, 等. 季节性高温和干旱对亚热带毛竹林碳通量的影响. 应用生态学报, 2016, 27(2): 335-344
Chen X⁃F, Jiang H, Niu X⁃D, et al. Effect of seasonal high temperature and drought on carbon flux of bamboo forest ecosystem in sub⁃
tropical region. Chinese Journal of Applied Ecology, 2016, 27(2): 335-344 (in Chinese)
443 应  用  生  态  学  报                                      27卷