基于涡动相关技术的森林生态系统二氧化碳通量研究进展

龚 元，纪小芳，花雨婷，张银龙，李 楠

（1.南京林业大学 生物与环境学院，江苏 南京 210037；2.南京林业大学 江苏省南方现代林业协同创新中心，江苏 南京 210037；3.阿拉巴马大学 生物科学系，阿拉巴马 塔斯卡卢萨 AL35487；4.南京林业大学 林学院，江苏 南京 210037；5.南京林业大学 风景园林学院，江苏 南京 210037）

18世纪60年代工业革命以来，由于全球土地利用变化、工业活动和化石燃料的使用，大量温室气体排放至大气中[1]。甲烷(CH4)和二氧化碳(CO2)等温室气体的大量排放，导致全球增温和气候变化，人类开始关注陆地生态系统碳循环过程和消失的碳汇(missing carbon sink)[1]。了解植物(草地、森林生态系统)对大气CO2吸收和排放的调节作用是人类应对全球气候变化的基础[1-2]。森林生态系统是陆地生态系统物质循环和能量流动的重要组成部分，开展不同生态尺度的森林生态系统碳循环研究对更好地了解陆地生态系统与大气间的碳交换过程，制定全球碳排放和碳交易政策(如欧洲碳排放交易体系等)，应对全球气候变化具有参考和服务作用[2-4]。涡动相关系统是目前应用较多的温室气体观测技术。本文主要依据涡动相关观测系统对森林生态系统CO2通量观测值、微气象数据等相关研究成果，讨论森林生态系统CO2通量观测原理和数据预处理、CO2通量动态特征、CO2通量环境影响因素、CO2通量足迹估算等方面的研究成果，为其他森林生态系统碳循环研究提供服务和方法参考。

1 CO2通量观测方法

研究森林生态系统碳循环过程就必须对其和大气间碳交换过程进行不同时间尺度/间隔的观测，其中对森林生态系统与大气间垂直碳交换即CO2通量观测是研究森林生态系统碳循环过程的重要研究内容[5-6]。CO2通量是指在一定生态尺度下单位时间单位面积内CO2流通的量，一般以μmol·m-2·s-1为单位，亦可依据不同的科学问题的需要对其质量单位、面积单位和时间单位进行换算[3]。对森林生态系统CO2通量进行观测的方式也较多，如自下而上(bottom-up)的涡动协方差/涡动相关技术和自上而下(topdown)的参数化遥感产品等[3]。20世纪90年代以来，涡动相关系统开始广泛应用于高大植被较多的森林生态系统CO2通量观测[7-8]。如今随着涡动相关系统的推广，为了数据共享和数据处理标准化形成了国际通量网(FLUXNET，https://fluxnet.fluxdata.org/)，中国通量网(ChinaFLUX，http://www.chinaflux.org/)等全球-区域的通量和微气象数据共享网络平台[9-10]。依据国际通量网2017年2月的统计结果，目前在其注册的站点全球共914个，其中较多站点为森林生态系统CO2通量观测站点，且多分布于温带地区。基于涡动相关系统森林生态系统CO2通量观测的研究也较多，如森林生态系统的碳源/汇估算、CO2通量观测源区/足迹的计算、CO2通量动态特征的提取、CO2通量预测/建模、地面观测值与遥感观测值的验证等[5-8,11-18]。开展基于涡动相关系统的森林生态系统CO2通量研究，可以为了解陆地生态系统碳循环过程，应对全球气候变化和评估全球碳平衡提供服务和参考[17]。

1.1 涡动相关系统工作原理

大气湍流是下垫面与大气间进行物质交换和能量流动的主要方式，为记录其过程中物质和能量流通的量，涡动相关系统便因此诞生[8]。该系统一般安装在地面边界层内，在该层中通量随高度变化的影响较小，在该层开展下垫面与大气的碳交换具有较高的空间代表性[8]。涡动相关系统一般以10～20 Hz的采样频率来记录一定观测高度(zm，m)上的微气象信息(三维风速/风向、超声虚温等)和某物质(CO2、H2O和CH4等)流通的量，后一般取一定时间间隔内(30 min)的平均值来记录微气象数据和通量数据，其中已经假设垂直风速的时间平均值为0 m·s-1，CO2通量可依据下式计算[8]：

式(1)中：FC为CO2通量(μmol·m-2·s-1或g·m-2·s-1)。w为垂直风速(m·s-1)，w'为垂直风速脉动量(m·s-1)，代表垂直风速瞬时值与平均值的偏差。ρC为CO2密度(kg·m-3)，ρC'为CO2密度脉动量(kg·m-3)，代表CO2密度瞬时值与平均值的偏差。公式横上线代表时间平均[8]。当FC＜0时代表所观测的区域/下垫面为碳汇(carbon sink)，当FC＞0时代表所观测区域/下垫面为碳源(carbon source)[8]。

1.2 通量观测源区与足迹

通量源区(source area)即代表涡动相关系统所测量CO2通量的下垫面来源。因涡动相关系统的观测受到下垫面、大气边界层环境等因素的影响，观测到的CO2通量只能代表传感器上风向一定区域内的状况，因此在使用这些CO2通量观测值时需要对CO2通量的观测源区的空间代表性进行评估[14]。

通量源区的足迹函数(footprint)代表上风向下垫面中若干个点源所形成的源区对在zm高度上所观测CO2通量F(0, 0,zm)的贡献强度。可由下式计算[18]：

式(2)中：通量观测点为原点(0, 0)，x轴指向来风方向，F0(x,y,0)为上风向源区中某一个点源(x,y)的通量贡献强度[源强，量纲与F(0, 0,zm)一致]。φ(x,y,zm)即为通量足迹函数(flux footprint predictions，FFP)[16]，足迹函数的量纲为m-2。通量足迹函数的计算结果主要包括上风向通量足迹贡献峰值的位置(m)，以及不同通量足迹贡献百分比(10%～95%)在上风向的位置(m)等信息[16]。由于当通量足迹贡献率为1时，通量源区为无限大，因此在进行通量源区的空间代表性的评估时，一般将通量足迹贡献率设定为80%～90%[16]。近年来，随着通量足迹函数的发展，出现了较多优秀通量足迹计算模型，依据其计算模式可分为解析式模型、拉格朗日随机扩散模型、大涡模拟模型、闭合模型等，具体为FSAM模型[19]、KM01模型[18]、KLJUN模型[16]等。KLJUN模型是目前较新颖的基于尺度(量纲)的通量足迹模型，也是目前应用至森林生态系统通量足迹计算较多的模型之一[14]，KLJUN等[16]提供了该通量足迹计算模型的在线计算平台(http://footprint.kljun.net/index.php)，并且还提供了该通量足迹计算模型在Matlab和R语言软件平台下的开源代码，可供研究人员进行下载和使用。

1.3 CO2通量预处理与数据质量控制

在使用和分析涡动相关系统对森林生态系统的CO2通量观测值时，需要对CO2通量数据进行质量控制和数据插补。CO2通量数据预处理过程如下表1所示[8]。为了方便和简化CO2通量数据的预处理过程，美国LI-COR公司(https://www.licor.com/)发展了几套基于涡动相关系统的配套通量数据计算软件，包括开源式涡动相关通量处理软件EddyPro。EddyPro集成了包括通量数据修正、质量控制、野点去除、通量足迹计算等预处理过程，使用方便。另外较新的通量处理软件Tovi，除了通量数据预处理功能之外还集成了通量源区绘制、数据插补等数据可视化功能。除了美国LI-COR公司之外，由德国马克斯普朗克生物地球化学研究所发展的通量数据在线预处理，插补和制图工具(https://www.bgc-jena.mpg.de/bgi/index.php/Main/HomePage)也为CO2通量数据和微气象数据的使用和分析提供了方便。

表1 CO2通量数据预处理方法Table 1 CO2 flux data preprocessing method

在CO2通量数据插补方面一般采用平均日变化、查表法和非线性回归法等[18]进行。其中非线性回归法是基于CO2通量与环境因子关系的经验模型，主要包括基于白天生态系统净交换量和光合有效辐射量对Michaelis-Menten光响应曲线模型和采用夜间生态系统呼吸的Arrhenius模型等[5-6]。一般在森林生态系统CO2通量插补时依据应所缺失通量数据的时间长度来确定插补方法，对缺失时间长度较短(≤14 d)使用平均日变化法和查表法[8]，对缺失时间长度较长的(>14 d)则使用非线性回归法，如使用基于光合有效辐射的Landsberg模型[20]和基于土壤温度(Ts)的Arrhenius模型等[21]。

2 涡动相关系统在森林生态系统的应用

在完成通量数据质量控制、通量数据插补、通量源区空间代表性检验后便可开始分析和使用所观测到的局域尺度(local-scale)森林生态系统CO2通量数据，亦可称之为森林生态系统净CO2交换，也可将森林生态系统净CO2交换数据拆分为生态系统呼吸和总初级生产力应用于森林生态系统碳循环过程研究。目前，涡动相关系统在森林生态系统的应用也有较多，依据国际通量网的统计结果，全球大部分涡动相关通量观测站点位于温带地区，且多为森林通量站点[14]。目前，基于涡动相关系统的森林生态系统碳循环研究的数据来源主要分2类，即局域测量和国际通量网的共享通量数据。基于这些通量和气象数据的分析主要包括：CO2通量数据标准化、CO2通量特征分析、CO2通量建模、通量足迹分析和多站点大数据分析等[18-30]。由于树龄、地理和气候等因素的影响，不同森林生态系统的碳源/汇状态不同。一般认为，森林生态系统是陆地生态系统的重要碳汇[25]。由于通量塔架设和维护成本较高，基于涡动相关技术森林生态系统的碳循环研究多依托单个通量塔的观测数据，多站点、长时间序列的数据整合分析较少。涡动相关系统对地形、下垫面完整度、微气象环境等因素要求较高，通量数据的有效率研究依旧是重要研究内容。由于通量数据插补方法的差异，同一时间序列的通量数据会由于不同的数据插补方法对生态系统与大气间碳交换估算产生影响，即影响生态系统碳源/汇功能的判定，因此开展更多通量数据插补方法研究是通量数据质量控制和CO2通量研究的保障和前提。

3 森林生态系统CO2通量动态特征及其影响因子

基于涡动相关系统对森林生态系统CO2通量研究包括CO2通量特征和通量源区的研究，其中CO2通量源区的计算包括足迹模型验证、足迹气候态和不同景观类型通量贡献差异等研究。CO2通量动态特征的讨论基于不同时间尺度、植物群落类型、气候条件等背景下森林生态系统与大气间碳交换的特征。

3.1 森林生态系统CO2通量足迹分析

涡动相关系统对CO2通量观测受到多种环境因素的影响，例如：观测高度、下垫面、大气稳定度、大气边界层高度等。因此，在使用CO2通量数据前需要对通量源区进行检验，即所观测的CO2通量值是否来源于感兴趣区域，即通量足迹分析[14]。

BALDOCCHI等[29]解释了拉格朗日随机模型(Lagrangian model)在高大植被冠层中通量足迹的应用，认为拉格朗日随机模型可应用于森林冠层上下2个部分的通量足迹分析，并强调在应用拉格朗日随机模型进行高大植被林冠层通量足迹分析时需要服从高斯力作用(Gaussian forcing)。KLJUN等[16]基于拉格朗日随机扩散模型所发展的三维拉格朗日足迹模型是目前应用较多的参数化通量足迹模型之一，并且引入了大气边界层高度的概念，认为在当通量足迹贡献率大于80%时即具有空间代表性。KIM等[14]基于KLJUN模型分析了美国哈佛森林多年通量足迹形态的气候态特征，认为森林景观变化直接影响了通量源区的形态和面积，且认为高大植被较多的森林生态系统CO2通量观测源区对观测高度(zm)即通量塔视场的依赖性较大。其他通量足迹模型例如FSAM模型、KM模型和HSIEH模型等也是应用至森林生态系统较多的模型。金莹等[30]基于FSAM模型分析了杨树Populus林的通量源区，认为通量源区的大小主要受到大气稳定度(zm/L，L为莫宁-奥布霍夫长度)影响，当大气处于稳定状态时通量源区面积大于当大气处于不稳定状态下的通量源区面积，且生长季的通量源区面积大于非生长季的通量源区面积。龚笑飞等[31]基于FSAM模型对安吉毛竹Phyllostachys edilus林通量源区的分析与金莹等[30]的分析结果类似。NEFTEL等[18]基于KM01模型所开发的ART footprint tool首次增加了通量源区的区域选取功能，即可完成对多个感兴趣区域通量足迹贡献率的计算，首次应用至农田和草地生态系统，并且认为该通量足迹计算工具可应用至下垫面较均一的环境中的其他生态系统(如森林生态系统)的通量足迹检验和CO2通量数据质量检验。OGUNJEMIYO等[32]在对黑云杉Picea mariana和短叶松Pinus banksiana混交林进行通量足迹分析时结合卫星和航空器观测数据对下垫面不同景观结构的CO2通量源区与通量观测值的预测。张慧[33]基于FSAM模型、KM01模型和HSIEH模型共3个通量足迹模型分析了千烟洲中亚热带红壤丘陵区人工林的CO2通量足迹，认为参数化方式和建模方式的不同导致了3个通量足迹模型计算结果的差异，且同样认为通量源区的大小对观测高度的依赖性较大。这一结论与KIM等[14]的研究结果一致。

通量源区/足迹分析是检验CO2通量数据来源的重要方法，是CO2通量数据预处理的必要过程之一[8]。通常在搭建通量塔前需要对局域尺度的上风向风区和通量源区进行估算，以及下垫面冠层高度的计算后再确定通量塔位置和传感器架设高度[8]。通量数据处理软件EddyPro也集成了相关的通量足迹计算功能，包括10%～90%通量足迹最远点在上风向的距离，通量足迹贡献峰值在上风向的距离等[16]。研究人员可依据下垫面实际情况使用相关的软件或者平台，对局域尺度的CO2通量源区进行估算，亦可配合CO2通量观测值分析下垫面中不同森林景观类型对大气CO2的贡献强度。

3.2 森林生态系统CO2通量时间动态变化特征

由于森林生态系统受到光照、气温、降水等环境控制因素的影响，且由于自身光合作用和呼吸作用，一般局域尺度森林生态系统CO2通量的日变化特征呈现“U”形特征。即一般在当地时间的正午时分(12:00-14:00)达到碳汇的最大强度，夜间由于生态系统呼吸作用，CO2通量观测值达到最大值。白天由于植物光合作用吸收大气CO2，导致涡动相关系统CO2通量30 min观测值(Re)为负值。夜间整个生态系统进入呼吸作用，导致涡动相关系统CO2通量30 min观测值为正值，且夜间CO2通量变化较平稳[3]。

局域尺度的森林生态系统CO2通量观测一般依托于通量观测塔(tower-based)[34-39]，进行自下而上的森林下垫面与大气间CO2交换观测。单个森林通量观测站点会搭配1套涡动相关系统、1套微气象梯度观测系统以及1个数据采集器。部分通量观测站点会在通量塔不同高度搭载多套红外气体分析仪和三维超声风速温度测量仪器，以满足在不同大气条件下(稳定、中性和不稳定)对感兴趣区域进行均一性观测[25]。

森林植被高度决定通量观测塔架设高度。森林生态系统的通量观测塔架设高度一般高于城市、湿地和草地等陆地生态系统[35]。目前，已知较高的森林通量塔安装在美国威斯康星州，观测高度为距离地面447.0 m[35]。位于北京市八达岭林场的通量塔架设高度较低(11.7 m)，其原因是该地区下垫面树种的平均高度较低(约4.0 m)，其较低的下垫面冠层高度决定了通量塔高度[34]。不同森林生态系统由于地理位置和环境差异，导致其CO2通量特征有一定空间分异性，其中主要表现为热带地区森林生态系统的碳汇峰值最强(-1.32 mg·m-2·s-1)，且热带地区的CO2通量最大值也高于其他温带森林生态系统CO2通量峰值(0.66 mg·m-2·s-1)[38]。森林通量站点中CO2通量极值一般出现于当地的生长季，其原因可能是下垫面中的绿色植物在生长季产生了较强的光合作用，但在夜晚植物自身生长也增加了生态系统呼吸作用。目前，对森林生态系统CO2通量动态特征的相关研究也已经比较成熟，一般认为，在日尺度下森林生态系统在白天由于光合作用形成碳汇，夜间由于呼吸作用形成碳源，在更长的时间尺度例如多年局域尺度，森林生态系统CO2通量估算中热带雨林表现为较稳定的碳汇，但是在生长季/雨季中会表现为碳源。如张一平等[38]研究的森林生态系统，由于阴雨天气较多导致光合有效辐射降低，地上凋落物增加也增强了土壤呼吸，使局域森林生态系统在生长季/雨季表现为碳源，但在全年尺度上依然表现为微弱的碳汇。在温带森林生态系统中，唐祥等[34]认为：成熟林因为较高的叶面积指数比新造林具有更强的碳汇潜力，且其研究的中国北京市八达岭林场CO2通量特征表现为7月的日CO2吸收值最大。在亚热带森林生态系统中，纪小芳等[27]研究的浙江凤阳山混交林生态系统在全年角度表现为碳汇，在7月净CO2交换水平达到最低，碳汇强度达到最大。

开展基于涡动相关系统在森林生态系统CO2通量观测，可以为长时间序列森林生态系统碳交易(carbon budget)、碳平衡(carbon balance)进行评估，森林通量观测站点长时间序列CO2通量数据积累可以帮助我们更好地理解森林生态系统对大气CO2的调节，即碳源和碳汇状态转换的时空特征。BRACHO等[25]基于涡动相关系统观测的美国东南部种植的湿地松Pinus elliottii林1998-2008年的CO2通量数据，分析了多年尺度湿地松林的碳平衡状态，结果表明：相对年轻的湿地松林分在1998-2001年表现为碳汇，年碳交易分别为-12.68、-8.85、-5.28、-2.37 Mg·hm-2·a-1，而相对老龄的林分表现为碳源。研究结果显示：年平均叶面积指数为1～2 m2·m-2时，年轻的林分表现为较强的碳汇。随着年平均叶面积指数的上升，林分渐渐转换为碳源。相对老龄林分的年叶面积指数为6～8 m2·m-2，且在年尺度上均表现为碳源，年碳交易为4.91～8.18 Mg·hm-2·a-1，该局域尺度生长季降水是影响湿地松林生态系统年碳交易的主要因素。PITA等[23]分析了4类不同树种林分CO2通量的差异和季节分异，认为在长时间序列下桉树Eucalyptus有效的叶片气孔控制可以应对地中海气候下饱和水汽压差的变化，适宜在地中海地区种植。杨树Populus在干旱期间的主要表现为叶面积指数下降。由于温带地区的苏格兰松Pinus sylvestris林有地下水位的补给导致受到水分胁迫的压力较小，且长时间序列的总初级生产力受气候变化的影响较小，但是4类林分中月均尺度总初级生产力最大值出现在5月的洛克里斯蒂的苏格兰松林，为7.9 g·m-2·d-1，表现为较强吸收CO2的能力。王春林等[28]所研究的鼎湖山常绿针阔叶混交林碳交易表明：在2003和2004年该局域尺度森林生态系统表现为碳汇，年碳交易分别为-563.0，-441.2 g·m-2·a-1，这与欧洲通量网的大部分森林站点的结果一致，且认为在夜间大气稳定度高的状态下土壤和森林冠层呼吸作用所产生的CO2不能被传感器有效观测到，可能会对CO2通量观测值产生影响。由于涡动相关系统对夜间生态系统呼吸的低估可能会导致高估该生态系统年尺度下的碳汇能力，其他中国通量网的注册站点，例如：LIU等[39]基于涡动相关系统研究千烟洲常绿针叶林碳交易，在2003年同样表现为碳汇；张一平等[38]研究的西双版纳热带季雨林在2003和2004年也表现为微弱碳汇。

3.3 森林生态系统CO2通量影响因子

森林生态系统碳循环受到环境以及气候变化的影响较大。基于涡动相关系统CO2通量观测可结合梯度微气象观测系统的气象观测值如气温、降水、土壤温度、水汽压差和光量子通量密度(DPPF)等环境因子与CO2通量进行关系分析，研究环境控制因子中对森林生态系统碳循环影响和森林下垫面对气候变化的响应[8]。目前，研究环境因子与CO2通量关系一般通过CO2通量30 min观测值(μmol·m-2·s-1，g·m-2·s-1)，CO2通量日累积值(μmol·m-2·d-1，g·m-2·d-1)和CO2通量月累积值(μmol·m-2·月-1，g·m-2·月-1)与环境控制因子进行回归分析建立线性或者非线性的关系模型。这里的环境因子包括不同时间尺度的空气温平均值(°)、降水量累积值(mm)、光合有效辐射平均值(μmol·m-2·s-1)等[3]。

不同森林生态系统对以上环境因子变化的响应不同。目前，一般认为气温是森林生态系统碳循环的主要影响因子[3]，降水会影响光合有效辐射后对森林生态系统与大气间CO2交换造成影响[38]，其他环境因子对森林生态系统碳循环的影响包括生态系统呼吸、净CO2交换量、总初级生产力、净初级生产力等[3]。

叶面积指数、增强植被指数、归一化植被指数等[40]植被指数产品(https://modis.gsfc.nasa.gov/data/dataprod/)也开始被应用于森林生态系统碳循环的辅助测量，配合地面通量塔的观测数据，可研究植被生长状况对局域尺度生态系统CO2通量的影响，以及研究植被生理活动对气候变化的响应。

4 森林生态系统CO2通量模型研究

涡动相关系统CO2通量观测值除了用于分析CO2通量动态变化的现象特征外，亦可将观测数据用于建立生态过程模型(统计模型/机理模型)，以及将CO2通量观测值作为地面参考数据(ground reference data)拆分出总初级生产力、净初级生产力等通量数据与遥感观测数据，进行交叉验证，推动全球森林生态系统碳循环的整合分析，完成从局域尺度CO2通量测量到全球生态系统碳循环研究的尺度上推。

4.1 基于统计模型的森林生态系统CO2通量研究

基于涡动相关系统CO2通量观测使用统计模型，即非线性/线性逻辑模型等方法进行生态过程建模，主要包括植物群落光合作用模型、威布尔累积分布函数、参数化逻辑函数和人工神经网络模型等[24]，其函数功能主要包括下垫面植物生理/物候过程的提取、CO2通量值预测、CO2通量观测值与遥感数据的交叉验证等。RICHARDSON等[41]使用国际通量网多站点森林CO2通量(GEP/NEP)，基于物候模型分析了温带森林生态系统CO2通量对全球气候变化的响应，认为全球增温会延长温带森林生态系统生长季长度。GU等[24]开发了植物群落光合作用物候模型(phenology model)。该模型主要基于单日总初级生产力30 min观测最大值，提取其动态变化的关键时间点，来划分局域尺度生态系统的物候阶段。其物候阶段主要分为：准备阶段(pre-phase)、恢复阶段(recovery phase)、稳定阶段(stable phase)、衰落阶段(senescence phase)、终止阶段(termination phase)，以及计算生长季长度、生长季中点、生长季开始和结束点等物候特征参数。NIU等[42]使用国际通量网CO2通量数据，基于物候模型，分析了北半球生态系统CO2通量对年均气温变化的响应，认为所分析的森林通量站点的物候特征均对年均气温变化的敏感性较高。GONSAMO等[43]使用国际通量网数据，基于一个非线性逻辑函数(double logistic function)分析了北美温带森林生态系统物候特征，并且基于遥感数据(NDVI/NDII)开发了一个新的物候指数(PI)，在将该物候指数与通量数据的交叉验证中，遥感数据与涡动通量数据显示了较好的一致性，均较好地体现了下垫面森林景观的物候过程。基于逻辑算法的下垫面物候特征参数提取还有很多其他的函数可用，例如：非对称高斯函数、D-L拟合和S-G滤波法等[44]。在人工神经网络模型的应用方面，MANCUSO等[45]讨论了人工神经网络在橄榄树Olea europaea物候过程建模的潜力，通过训练反向传播神经网络并对其进行了测试，预测了下垫面中橄榄树生物事件的时间节点，并且结合了微气象数据预测橄榄树物候对未来气候变化的响应。PAPALE等[46]使用EUROFLUX等通量数据来训练神经网络模拟器(neural network simulator)，讨论了欧洲范围内森林碳通量的空间(1 km×1 km)和时间(周)格局，认为在欧洲南部的生长季更长(约32周)。HE等[47]使用人工神经网络技术，基于中国通量网中3个不同生态系统(森林、草地和农田)的CO2通量和微气象数据，来训练人工神经网络并预测CO2通量，结果表明：该技术可以成功地预测生态系统CO2通量，相关系数为0.75～0.86。上述模型和方法均可应用于森林生态系统物候特征参数提取的研究，一般统计模型多基于1个或者多个非线性逻辑函数对CO2通量观测数据进行拟合分析，并且依据拟合后的CO2通量变化斜率的计算划分不同的物候阶段。这些方法对温带、亚热带(暖季)植被的拟合效果较好。由于热带地区常绿植物较多，物候特征不明显，上述物候模型在热带地区的应用效果差于温带和亚热带地区。GU等[24]开发的物候模型集成了CO2通量增长率计算的功能，并且针对植物群落生物事件划分了多个关键物候阶段，是具有代表性的统计模型之一。近年来，随着遥感技术的发展，MODIS遥感产品为科研工作者提供了总初级生产力和净初级生产力等遥感植被产品，其中总初级生产力的产品周期为8 d，净初级生产力的产品周期为1 a[48-49]。基于遥感数据与地面观测数据的交叉验证也较多，且较多集中于总初级生产力产品的扩展应用[48]。杜启勇等[49]同样基于遥感和地面通量数据进行了验证分析，认为在森林生态系统总初级生产力估算值和实测值在年碳交易尺度上误差较小，但低估了其他生态系统的总初级生产力。刘啸添等[50]基于遥感数据和地面通量观测数据提取了局域尺度的温带针阔混交林的物候特征，认为叶绿素荧光与地面观测总初级生产力值的物候特征一致性较高，较真实地反应了地面植被的生长/生理过程。遥感数据与地面观测数据的交叉验证，可以为提高遥感数据观测精度提供参考，也方便了局域尺度CO2通量观测的尺度上推[51-52]。

4.2 基于机理模型的森林生态系统CO2通量研究

基于植物生理/机理的CO2通量模型研究多集中于净CO2交换数据拆分、CO2通量预测、CO2通量数据插补方法等研究。目前应用较多的植物生理模型主要包括CANVEG模型、CANOAK模型、CHANGE模型和BKPF模型等[53-55]。唐欢等[48]基于FLUXNET和中国通量网的涡动相关通量观测数据和遥感总初级生产力产品进行了交叉验证，认为在生长季中期出现了遥感产品与地面数据不一致的情况，主要原因可能是模型输入参数对环境变量的敏感度差异。LAI等[53]基于CANVEG模型对美国北卡罗来纳州的杜克森林的CO2通量分布进行预测和模拟。MEDLYN等[54]基于涡动相关系统CO2通量观测数据评估了一套生态系统碳交换模型，认为其中的总初级生产力模型(sun-shade GPP model)表现较好，2套生态系统呼吸模型(RE model 1/2)即基于土壤温度的呼吸模型和基质回收模型，由于受到局域尺度云杉Picea asperata林的环境差异限制不能完全表现真实的生态系统呼吸情况。BALDOCCHI等[55]基于生物物理模型CANOAK，对一个温带落叶林生态系统和大气之间的CO2交换进行了不同时间尺度的模拟，认为生长季长度会影响森林生态系统与大气间碳交易估算，主要表现为生长季节长度每增加1 d会使生态系统的净CO2交换量减少5.9 g·m-2。PARK等[56]使用生物地球化学模型(CHANGE)模拟了1998-2006年期间西伯利亚东部的落叶松Larix cajanderi林生态系统与大气间CO2交换，以了解该生态系统对气候变化的响应及其控制因素，研究结果表明：生态系统净CO2交换变化与净初级生产力密切相关，且土壤含水量是影响该落叶松林生态系统中CO2通量的决定因素。XIE等[57]使用通量塔观测数据结合CLM模型，讨论了干旱对亚热带森林生态系统的影响，研究结果表明：由于长时间的干旱，生态系统总初级生产总值下降了76%，水分条件影响了该生态系统的碳固存能力。CHEN等[58]使用BKPF模型预测了中国东北寒带和温带森林交界处的森林生态系统与大气间CO2交换，讨论了在大气CO2浓度上升的背景下，不同林分叶面积指数和生物量的差异。SHI等[59]使用多层气孔导度耦合模型模拟了温带混交林生态系统CO2通量，研究结果表明：在2003年5月至2007年9月，通量模拟值与CO2通量观测值有较好的相关性(R2=0.734)。李雪建等[60]基于2个竹林生态系统通量观测站点2014-2015年的通量观测数据和遥感数据，使用BEPS模型预测2种竹林生态系统总初级生产力、净CO2交换和生态系统呼吸等碳循环数据，结果表明：使用BEPS模型的模拟结果较好地体现了下垫面碳循环过程。陈晨等[61]使用2011年帽儿山生态站CO2通量数据，基于卡尔曼滤波的顺序同化技术，对BEPS模型的关键参数进行优化，结果表明：该方法明显改善了模型模拟碳水通量的能力。杨延征等[62]利用集成生物圈IBIS模型模拟分析了1960-2006年中国陆地生态系统碳收支时空变异特征和趋势，认为大兴安岭、小兴安岭和长白山等地区的森林生态系统是较强的碳汇。王萍[63]使用IBIS模型对2004-2005年大小兴安岭的植被净初级生产力进行了定量估算，结果表明：大小兴安岭森林植被的年均净初级生产力为494.7 g·m-2·a-1，热量条件是其主要影响因子。王培娟等[64]使用BEPS模型分析了在长白山自然保护区森林植被净初级生产力，认为叶面积指数是影响森林净初级生产力的重要影响因子。以上关于CO2通量的建模在局域尺度的运用多有限制，不同生理模型对叶面积指数、气温和降水等环境因素的敏感性不同，且影响机制复杂。当需要对局域尺度CO2通量进行预测时，应该准确评估当地的实际环境情况，考虑水文、生理、生化和生态因素与气象条件之间的相互作用[56]，选择合适的CO2通量模型，亦可选择统计模型和机制模型交叉验证的方式进行森林生态系统CO2通量预测。

5 结论与展望

5.1 结论

涡动相关通量观测系统是目前应用最广泛的森林生态系统CO2通量观测技术，该系统可以直接观测森林冠层与大气间的碳交换，但是在应用该系统观测森林生态系统CO2通量时，需要对局域尺度的通量源区和上风风区延伸度进行估算。此后对CO2通量观测结果需进行必要的修正和通量数据质量控制，对不合格或者缺失的数据需依据实际情况进行插补。

森林生态系统CO2通量源区受观测高度的影响较大，下垫面森林景观的更迭会直接影响通量源区的形态和长度。CO2通量动态特征具有较明显的日、月和季节分异，具体表现在日尺度上30 min的CO2通量观测值多呈现“U”形特征，且森林生态系统在白天由于光合作用多表现为碳汇，夜间由于呼吸作用变为碳源。亚热带和温带森林系统的CO2通量月累积值约在当地时间夏季达到最低，即碳汇能力最强。热带森林生态系统较之其他气候带的森林生态系统表现为更强的碳汇。森林生态系统长时间序列碳交易受到叶面积指数、树龄、降水等多种因素的影响，年轻林分是更强的碳汇。

森林生态系统CO2通量数据的建模主要集中于CO2通量数值预测和下垫面物候特征的提取，植被生长曲线和通量值模拟主要通过统计模型和机制模型2种模式进行。基于CO2通量通过建立统计模型提取其物候特征是CO2通量建模的重要方向，统计模型所提取的物候特征参数主要包括：生长季长度、生长季中点、生长季开始和结束日期、CO2通量增长率等。基于机理模型的CO2通量数值预测/模拟受气候背景的影响较大，不同局域尺度的CO2通量预测需准确评估当地微气象环境情况和植被生长状况，其中空气温度、土壤含水量是主要环境控制因子。

基于遥感技术的辅助测量是森林生态系统碳循环研究的重要手段，针对遥感数据和地面观测数据的建模和验证主要集中在遥感总初级生产力和地面观测数据的线性相关性分析，其他例如叶面积指数、增强植被指数和归一化植被指数等遥感植被产品对辅助研究森林生长状态对其CO2通量的影响也是重要的碳循环研究内容。

5.2 研究展望

基于涡动相关技术的森林生态系统CO2通量观测研究为森林生态系统碳循环研究提供了观测手法，是目前使用最广泛且科学的温室气体观测的技术手段。基于该技术的CO2通量观测和研究分析是现在和未来全球变化生态学的研究热点。但由于受观测设备故障、恶劣天气事件和人为活动等因素的影响，CO2通量数据的质量检验和数据插补依然是未来CO2通量研究的前提之一。

森林生态系统CO2通量动态特征的现象讨论趋于成熟，未来局域尺度森林生态系统CO2通量研究应该集中于提高CO2通量足迹模型计算精度和讨论局域森林生态系统中不同林分对大气CO2的贡献强度，结合通量和微气象观测数据建立植物生理过程参数化模型和设计控制实验，提取其物候特征并且预测气候和景观格局变化对森林碳交换的影响。由于国际通量网和其他区域通量共享系统提供了数据平台，区域-全球尺度森林生态系统CO2通量研究未来将关注多站点通量，气象数据长时间序列的整合分析，主要讨论CO2通量气候态特征与碳源/汇的空间格局。为了解未来气候变化背景下全球植被生理活动过程和其温室气体排放动态，使用CO2通量数据与区域大气模型、生态生理模型等地球系统模型结合，讨论气象/气候系统、能量传输、生物圈对森林生态系统净CO2交换影响，亦可使用该过程对生态系统净CO2交换进行预测是未来重要的研究方向。

在中国生态文明建设的背景下，未来森林生态系统CO2通量研究可结合社会经济、土地利用等数据，讨论城市化过程对森林生态系统碳循环的影响，完善森林生态系统服务功能体系，评估局域自然-社会系统碳平衡，为合理森林布局、保护生态环境和社会可持续发展提供理论依据和参考。