杨树人工林生态系统通量贡献区分析

采取通量源区模型FSAM(Flux Source Area Model)利用2009年北京大兴区杨树人工林生态系统碳水通量涡度相关观测资料,分析了不同大气条件下生态系统的通量贡献区分布特征。研究结果表明:(1)该站通量贡献区随大气稳定条件增强而增大。除326.25°—3.75°方向生长季不稳定条件下外,生长季的通量贡献区范围普遍大于非生长季的贡献区范围;(2)通量贡献区与观测高度、冠层高度、地表粗糙度、风向以及大气稳定度有关,当风速风频较大,大气不稳定时,湍流扩散作用强烈,贡献区范围较小;(3)该观测场在2009年以不稳定大气条件为主,通量信息主要来源于距离观测塔50—400 m范围,且69.3%的信息来源于通量塔偏北风与偏南风方向,其中42.56%的信息来自于偏南风方向;(4)随着大气稳定程度加强,通量来源最少区从塔偏西方转为偏东方,在大气稳定度条件和风向的共同作用下,生长季时主要通量贡献区在塔偏南方向,而非生长季时主要通量贡献区在塔偏北方向。     

华北平原冬小麦农田生态系统通量贡献区

利用2013—2014年涡度相关系统观测的华北平原冬小麦农田生态系统通量数据,结合通量贡献区模型FSAM,分析华北平原冬小麦农田生态系统通量贡献区的时空分布特点,对比研究不同大气稳定层结条件和生长期内通量贡献区的分布差异.结果表明: 在主风风向上,冬小麦整个生育期内大气稳定条件下的通量贡献区范围大于不稳定条件下的贡献区范围.在0°～90°主风风向上,生长初期稳定条件下通量贡献区范围比不稳定条件下大17.8 m左右,生长末期稳定条件下的通量贡献区范围比不稳定条件下大11 m左右.生长初期的通量贡献最大值点位置比生长末期距观测点位置远15 m(大气稳定条件)和12.4 m(大气不稳定条件);通量贡献最大值点在稳定条件下比不稳定条件下距观测点位置远5 m(生长初期)和2.4 m(生长末期).在非主风风向上,当风向为90°～180°时,生长初期和生长末期不同大气条件下的最大通量值分别位于距观测点的67.8、53.4和47.0、30.8 m.当风向为270°～360°时,生长初期和生长末期不同大气条件下的最大通量值位于距观测点的58.8、42和41.1、33.1 m.在整个生育期尺度上,观测塔的通量信息主要来自东北、西南和东南方向,其所占比例分别为35.4%、32.5%和19.4%.冬小麦整个生育期内通量贡献区的主要变化发生在观测点东北方向16.0～173.8 m和西南方向14.7～209 m,通量信息全部来源于农田生态系统.两个典型日期的通量贡献区日变化特征明显,通量贡献区范围随大气稳定条件和风向改变而发生变化.夜晚通量信息全部来源于农田生态系统,白天少部分通量信息来源于居民区和果园.本文的定量化结果可为农田生态系统通量贡献区的研究提供依据.     

大气热层结条件对林冠显热的影响

应用欧拉二阶闭合模型研究了大气热层结条件下森林冠层显热通量源汇分布和通量特征.结果表明:白天,冠层上的不稳定层结和冠层下的稳定层结是森林冠层大气层结的一种特有现象;温度廓线的变化表明林冠高度2/3处存在较强的热源;冠层内大气处于弱稳定状态时,热量继续向上输送,呈现出热通量的反梯度输送.显热通量日变化的模拟值与实测值吻合,其R2=0.9035(P0.01).在显热收支方程中添加浮力项,可提高反演模型在实际大气中的模拟精度,从而改善模型对热通量收支的模拟能力.     

长白山阔叶红松林二氧化碳湍流交换特征

采用开路式涡度相关技术,研究了长白山阔叶红松林森林-大气界面的CO₂湍流交换特征.结果表明,在近中性大气层结条件下,冠层上方垂直风速和CO₂浓度功率谱在惯性子区基本符合-2/3定律,垂直方向主导湍涡尺度约为40 m.湍流通量贡献区主要在0.01～2 Hz频率范围内,冠层上方低频传输的湍涡贡献了更多的CO₂通量.这说明开路式涡度相关仪器系统可以满足冠层上方湍流通量观测的基本要求.但通过涡度相关法实测获得的森林-大气CO₂通量仍存在夜间低估现象,非湍流过程的增加是涡度相关技术应用的主要制约因素.因此,需要对弱湍流条件下的CO₂通量做相应的修订.     

中亚热带森林区通量观测的源面积探讨

利用 2002~2003年典型中亚热带森林地区 (千烟洲站 ) 通量观测资料, 运行FSAM (Flux_source_area_mod el) 解析模型初步确定该站的通量源面积。结果表明 :通量源面积随测量高度、测量时间和大气稳定度变化 ;风向的季节变化, 使得通量源面积分布方向也变化 ;该测站有足够大的风距, 测量结果能代表仪器所在地的现场特征, 即森林下垫面特征。另外, 相同时间的风向和通量数据绘图表明 :在Z/L满足模型运行要求时, 通量数据几乎集中在一个方向上 (随风向变化 ) 。模型结果在水平均匀湍流下得到, 满足模型基本要求。在该测站中, FSAM模型只能模拟昼夜层结稳定, 湍流不活跃时的通量源面积。今后重点是进一步查明森林内、外湍流特点, 改进模型, 使之对该测站的适用性更广。     

农林复合生态系统与低层大气间的通量研究

根据近年在黄淮海平原大面积农林复合区的观测资料,采用简单的一维模型,讨论了此类复合系统与低层大气之间的热量、动量和水汽通量,比较了不同的结构林网形成的通量差异。结果表明,由于林冠层的摩擦作用,农林复合生态系统上空风速随高度增大较快,动量输送总是向下的。在层结稳定的夜间,其数值较大,而在湍流比较活跃的正午前后,由于垂直风切变减弱,动量输送也相应减小。夜间农林作物冠层的辐射冷却加强,近地层逆温梯度增大,加之作物冠层内部因辐射冷却而造成水汽凝结,冠层上部水汽大于冠层内部,此时热通量和水汽通量均向下。这一现象在白天日出后开始逆转,日出后,地面升温加快,并逐渐形成超绝热梯度,致使热量通量向上。此时农林作物蒸腾加剧,冠层高度范围内水汽含量增大,造成水汽向上输送。文中还讨论了一维模式的缺点,指出了进一步研究生态系统与大气相互作用的重要性.     

凤阳山针阔混交林通量观测源区分布及特征

基于涡度相关技术的通量观测要求下垫面均匀一致,然而在实际观测中,地形往往非常复杂。尤其是丘陵山区森林植被覆盖的区域,不同植被类型的空间分布对该区域的碳通量观测影响很大。本研究利用Kljun模型和ART Footprint Tool,对浙江凤阳山针阔混交林森林生态系统2017年全年的观测数据进行分析,探究该区域涡度相关系统在不同季节、不同大气稳定度条件下的通量源区变化情况。结果表明：受地形的影响,研究区内全年盛行东北风（0-90°）和西南风（180-270°）;因此,通量的源区分布也主要在东北和西南方向;当通量贡献率达到90%时,源区长度最大不超过7000m,当通量贡献率在80%时,源区长度不超过3000m;在大气稳定条件下,其源区分布范围总是比大气不稳定条件下的源区范围广。在该源区分布范围内,主要以针阔混交林分为主,也分布有少量杉木、柳杉和毛竹等纯林,源区贡献从大到小依次为针阔混交林、阔叶林、杉木林、毛竹林、柳杉林、黄山松林。     

科尔沁草地不同大气稳定度下湍流特征谱分析

利用科尔沁半干旱风沙草原区涡度相关观测数据,计算并分析了不同大气稳定度条件下湍流的功率谱和协谱.结果表明:经向风速与纬向风速的功率谱在大气不稳定层结、近中性层结、稳定层结条件下谱形均相似,并遵循-2/3定律;近中性层结条件下垂直风速的功率谱没有惯性子区,此时仪器对脉动信号的捕捉不够完全,稳定层结条件下垂直风速功率谱表明观测层产生了泄流和平流;CO_2和H_2O浓度的功率谱分析表明,红外气体分析仪对高频信号的响应能力能够满足实际观测的需要,其功率谱的谱峰随大气稳定度的增强而向高频转移;垂直风速与3种标量的协谱谱形相似,不同大气稳定度下协谱惯性子区长度随大气稳定度的增强而变短,且均符合-4/3定律,说明三维超声风速仪与红外气体分析仪的空间间隔不会对通量观测结果造成显著影响.     

涡度相关技术及其在陆地生态系统通量研究中的应用

通量观测是定量描述土壤-植被-大气间物质循环和能量交换过程的基础。涡度相关技术作为直接测量植被冠层与大气间能量与物质交换通量的技术手段, 已经逐步发展成为国际通用的通量观测标准方法。随着涡度相关技术在全球碳水循环研究中的广泛应用, 长期连续的通量观测正在为准确评价生态系统碳固持能力、水分和能量平衡状况、生态系统对全球气候变化的反馈作用、区域和全球尺度模型的优化与验证、极端事件对生态系统结构与功能影响等方面的研究提供重要数据支撑和机制理解途径。通过站点尺度通量长期动态观测, 明确了不同气候区和植被类型生态系统碳水通量强度基线及其季节与年际变异特征。通过多站点联网观测, 在区域和全球尺度研究生态系统碳通量空间变异特征, 揭示了区域尺度上温度和降水对生态系统碳通量空间格局的生物地理学控制机制。该文概括地介绍了涡度相关技术的基本原理、假设与系统构成, 总结了涡度通量长期联网观测在陆地生态系统碳水通量研究中的主要应用, 并对通量研究发展前景进行了展望。     

涡度相关技术及其在陆地生态系统通量研究中的应用

通量观测是定量描述土壤-植被-大气间物质循环和能量交换过程的基础。涡度相关技术作为直接测量植被冠层与大气间能量与物质交换通量的技术手段,已经逐步发展成为国际通用的通量观测标准方法。随着涡度相关技术在全球碳水循环研究中的广泛应用,长期连续的通量观测正在为准确评价生态系统碳固持能力、水分和能量平衡状况、生态系统对全球气候变化的反馈作用、区域和全球尺度模型的优化与验证、极端事件对生态系统结构与功能影响等方面的研究提供重要数据支撑和机制理解途径。通过站点尺度通量长期动态观测,明确了不同气候区和植被类型生态系统碳水通量强度基线及其季节与年际变异特征。通过多站点联网观测,在区域和全球尺度研究生态系统碳通量空间变异特征,揭示了区域尺度上温度和降水对生态系统碳通量空间格局的生物地理学控制机制。该文概括地介绍了涡度相关技术的基本原理、假设与系统构成,总结了涡度通量长期联网观测在陆地生态系统碳水通量研究中的主要应用,并对通量研究发展前景进行了展望。     

区域尺度无人机涡动相关通量观测系统的应用研究

陆地生态系统的水热循环与碳循环是陆地表层系统中物质能量循环的核心,其中区域尺度地表水、热、碳通量的直接观测是当下陆地生态系统通量观测与模拟研究中的热点与难点。机载涡动相关方法能够直接观测区域尺度生态系统通量,基于无人机平台的涡动相关通量观测技术同时兼具了区域覆盖性与经济灵活性等优点,是机载通量观测技术的最新发展方向。在介绍机载涡动相关通量观测方法的主要技术原理、观测特点以及无人机通量观测系统组成的基础上,通过在相对均匀的区域开展无人机与地面通量观测对比试验,采用谱分析、观测结果对比以及源区分析等方式对无人机通量观测系统的性能进行了初步评价。结果表明：无人机通量观测系统能够实现对大气高频湍流信号的有效采样;无人机与地面观测的湍流通量具有较好的一致性,但是感热和CO₂通量出现了低估、潜热和摩擦风速出现了高估;观测平台与仪器的差异、垂直通量辐散、大气边界层条件、不同的地面源区及地表异质性的影响是造成二者差异的潜在主要因素。最后对未来研究目标进行了展望,以进一步推动该技术在相关领域中的应用。     

Towards physiologically meaningful water‐use efficiency estimates from eddy covariance data

Intrinsic water‐use efficiency (iWUE) characterizes the physiological control on the simultaneous exchange of water and carbon dioxide in terrestrial ecosystems. Knowledge of iWUE is commonly gained from leaf‐level gas exchange measurements, which are inevitably restricted in their spatial and temporal coverage. Flux measurements based on the eddy covariance (EC) technique can overcome these limitations, as they provide continuous and long‐term records of carbon and water fluxes at the ecosystem scale. However, vegetation gas exchange parameters derived from EC data are subject to scale‐dependent and method‐specific uncertainties that compromise their ecophysiological interpretation as well as their comparability among ecosystems and across spatial scales. Here, we use estimates of canopy conductance and gross primary productivity (GPP) derived from EC data to calculate a measure of iWUE (G₁, “stomatal slope”) at the ecosystem level at six sites comprising tropical, Mediterranean, temperate, and boreal forests. We assess the following six mechanisms potentially causing discrepancies between leaf and ecosystem‐level estimates of G₁: (i) non‐transpirational water fluxes; (ii) aerodynamic conductance; (iii) meteorological deviations between measurement height and canopy surface; (iv) energy balance non‐closure; (v) uncertainties in net ecosystem exchange partitioning; and (vi) physiological within‐canopy gradients. Our results demonstrate that an unclosed energy balance caused the largest uncertainties, in particular if it was associated with erroneous latent heat flux estimates. The effect of aerodynamic conductance on G₁ was sufficiently captured with a simple representation. G₁ was found to be less sensitive to meteorological deviations between canopy surface and measurement height and, given that data are appropriately filtered, to non‐transpirational water fluxes. Uncertainties in the derived GPP and physiological within‐canopy gradients and their implications for parameter estimates at leaf and ecosystem level are discussed. Our results highlight the importance of adequately considering the sources of uncertainty outlined here when EC‐derived water‐use efficiency is interpreted in an ecophysiological context.     

中国陆地生态系统通量观测站点空间代表性

涡度相关技术是测定大气与陆地生态系统之间CO₂交换、水分和能量通量最直接的方法,可用于研究土壤、植被与大气间的CO₂交换及其调控机制。收集了11个影响净碳交换量的主要变量信息,包括气象因素、土壤因素和地形因素的非生物因子、实际植被状态以及植被生产力,采用多元地理变量空间聚类分析方法,绘制出不同聚类数（25、50、75、85、100、150和200类）的通量生态区。结合中国现有通量观测站点的空间分布格局,与新生成的通量生态区和已有的自然地理区划进行对比分析,发现由于中国地形复杂,生态系统类型多样,现有85个涡度相关通量观测站点仅能刻画部分中国生态系统类型的净碳交换量时空特征,通量生态区划分为100-150类比较合适。考虑到涡度相关通量观测运行成本,通量站点可增加至150个,从而使得优化后的通量观测网络能够代表中国主要类型的生态系统,并且有利于通量观测数据与遥感资料的有效结合,提高碳水通量观测从站点扩展到区域尺度的精度,从而更好地检验过程机理模型的模拟结果。     

Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: past, present and future

The eddy covariance technique ascertains the exchange rate of CO₂ across the interface between the atmosphere and a plant canopy by measuring the covariance between fluctuations in vertical wind velocity and CO₂ mixing ratio. Two decades ago, the method was employed to study CO₂ exchange of agricultural crops under ideal conditions during short field campaigns. During the past decade the eddy covariance method has emerged as an important tool for evaluating fluxes of carbon dioxide between terrestrial ecosystems and the atmosphere over the course of a year, and more. At present, the method is being applied in a nearly continuous mode to study carbon dioxide and water vapor exchange at over a hundred and eighty field sites, worldwide. The objective of this review is to assess the eddy covariance method as it is being applied by the global change community on increasingly longer time scales and over less than ideal surfaces. The eddy covariance method is most accurate when the atmospheric conditions (wind, temperature, humidity, CO₂) are steady, the underlying vegetation is homogeneous and it is situated on flat terrain for an extended distance upwind. When the eddy covariance method is applied over natural and complex landscapes or during atmospheric conditions that vary with time, the quantification of CO₂ exchange between the biosphere and atmosphere must include measurements of atmospheric storage, flux divergence and advection. Averaging CO₂ flux measurements over long periods (days to year) reduces random sampling error to relatively small values. Unfortunately, data gaps are inevitable when constructing long data records. Data gaps are generally filled with values produced from statistical and empirical models to produce daily and annual sums of CO₂ exchange. Filling data gaps with empirical estimates do not introduce significant bias errors because the empirical algorithms are derived from large statistical populations. On the other hand, flux measurement errors can be biased at night when winds are light and intermittent. Nighttime bias errors tend to produce an underestimate in the measurement of ecosystem respiration. Despite the sources of errors associated with long‐term eddy flux measurements, many investigators are producing defensible estimates of annual carbon exchange. When measurements come from nearly ideal sites the error bound on the net annual exchange of CO₂ is less than ±50 g C m^?2 yr^?1. Additional confidence in long‐term measurements is growing because investigators are producing values of net ecosystem productivity that are converging with independent values produced by measuring changes in biomass and soil carbon, as long as the biomass inventory studies are conducted over multiple years.     

A new model of net ecosystem carbon exchange for the deciduous-dominated forest by integrating MODIS and flux data

The eddy covariance technique provides continuous measurements of plot-level net ecosystem carbon exchange (NEE) across a wide range of vegetation types. However, these NEE estimates only represent fluxes at the tower footprint scale. To quantify the NEE over regions or continents, flux tower measurements need to be up-scaled to large areas. In the present study, we propose a new NEE model solely based on Moderate Resolution Imaging Spectroradiometer data, including enhanced vegetation index (EVI), land surface water index (LWSI), land surface temperature (LST), and Terra nighttime LST′. Site-specific data from the deciduous-dominated Harvard Forest flux site were used. Analysis covered six years (2001-2006) of CO₂ flux data. The data of the first four years were used for model building and the rest as validation set. Compared with the model based solely on EVI, we also introduced LST and LSWI into the new model. The results showed that this method could further improve the precision (R² and RMSE reached 0.857 and 1.273, respectively) and generally capture the expected seasonal patterns of NEE.     

FLUXNET and modelling the global carbon cycle

Measurements of the net CO₂ flux between terrestrial ecosystems and the atmosphere using the eddy covariance technique have the potential to underpin our interpretation of regional CO₂ source–sink patterns, CO₂ flux responses to forcings, and predictions of the future terrestrial C balance. Information contained in FLUXNET eddy covariance data has multiple uses for the development and application of global carbon models, including evaluation/validation, calibration, process parameterization, and data assimilation. This paper reviews examples of these uses, compares global estimates of the dynamics of the global carbon cycle, and suggests ways of improving the utility of such data for global carbon modelling. Net ecosystem exchange of CO₂ (NEE) predicted by different terrestrial biosphere models compares favourably with FLUXNET observations at diurnal and seasonal timescales. However, complete model validation, particularly over the full annual cycle, requires information on the balance between assimilation and decomposition processes, information not readily available for most FLUXNET sites. Site history, when known, can greatly help constrain the model‐data comparison. Flux measurements made over four vegetation types were used to calibrate the land‐surface scheme of the Goddard Institute for Space Studies global climate model, significantly improving simulated climate and demonstrating the utility of diurnal FLUXNET data for climate modelling. Land‐surface temperatures in many regions cool due to higher canopy conductances and latent heat fluxes, and the spatial distribution of CO₂ uptake provides a significant additional constraint on the realism of simulated surface fluxes. FLUXNET data are used to calibrate a global production efficiency model (PEM). This model is forced by satellite‐measured absorbed radiation and suggests that global net primary production (NPP) increased 6.2% over 1982–1999. Good agreement is found between global trends in NPP estimated by the PEM and a dynamic global vegetation model (DGVM), and between the DGVM and estimates of global NEE derived from a global inversion of atmospheric CO₂ measurements. Combining the PEM, DGVM, and inversion results suggests that CO₂ fertilization is playing a major role in current increases in NPP, with lesser impacts from increasing N deposition and growing season length. Both the PEM and the inversion identify the Amazon basin as a key region for the current net terrestrial CO₂ uptake (i.e. 33% of global NEE), as well as its interannual variability. The inversion's global NEE estimate of −1.2 Pg [C] yr⁻¹ for 1982–1995 is compatible with the PEM‐ and DGVM‐predicted trends in NPP. There is, thus, a convergence in understanding derived from process‐based models, remote‐sensing‐based observations, and inversion of atmospheric data. Future advances in field measurement techniques, including eddy covariance (particularly concerning the problem of night‐time fluxes in dense canopies and of advection or flow distortion over complex terrain), will result in improved constraints on land‐atmosphere CO₂ fluxes and the rigorous attribution of mechanisms to the current terrestrial net CO₂ uptake and its spatial and temporal heterogeneity. Global ecosystem models play a fundamental role in linking information derived from FLUXNET measurements to atmospheric CO₂ variability. A number of recommendations concerning FLUXNET data are made, including a request for more comprehensive site data (particularly historical information), more measurements in undisturbed ecosystems, and the systematic provision of error estimates. The greatest value of current FLUXNET data for global carbon cycle modelling is in evaluating process representations, rather than in providing an unbiased estimate of net CO₂ exchange.     

Aircraft‐derived regional scale CO2 fluxes from vegetated drained thaw‐lake basins and interstitial tundra on the Arctic Coastal Plain of Alaska

The landscape surface of the Barrow Peninsula of Alaska is a mosaic of small ponds, thaw lakes, different aged vegetated drained thaw‐lake basins (VDTLBs), and interstitial tundra which have been dynamically formed by both short‐ and long‐term processes. We used a combination of tower‐ and aircraft‐based eddy covariance measurements to characterize the spatial and temporal patterns of CO₂, latent, and sensible heat fluxes along with MODIS NDVI, and were able to scale the aircraft‐based CO₂ fluxes to the 1802 km² Barrow Peninsula region. During typical 2006 summer conditions, the midday hourly CO₂ flux over the region was ?2.04 × 10⁵ kg CO₂ h^?1. The CO₂ fluxes among the interstitial tundra, Ancient, and Old VDTLBs, as well as between the Medium and Young VDTLBs were not significantly different. Combined, the interstitial tundra and Old and Ancient VDTLBs represent～67% of the Barrow Peninsula surface area, accounting for ～59% of the regional flux signal. Although the Medium and Young VDTLBs represent ～11% of the surface area, they account for a large portion, ～35%, of the total regional flux. The remaining ～22% of the surface area are lakes and contributed the remaining ～6% of the total regional flux. Previous studies treated vegetated areas of the region as a single surface type with measurements from a few study sites; doing so could underestimate the regional flux by ～22%. Here, we demonstrate that aircraft‐based systems have the ability to cover large spatial scales while measuring the turbulent fluxes across a number of surfaces and combined with ground‐ and satellite‐based measurements provide a valuable tool for both scaling and validation of regional‐scale fluxes.     

Identification of vegetation and soil carbon pools out of equilibrium in a process model via eddy covariance and biometric constraints

Assumptions of steady‐state conditions in biogeochemical modelling are often invoked because knowledge on the development status of the modelling domain is generally unavailable. Here, we investigate the role of vegetation pool sizes on nonequilibrium conditions through model‐data integration approaches for a set of sites using eddy covariance CO₂ flux data. The study is based on the Carnegie–Ames–Stanford Approach (CASA) model, modified (CASA_G) in order to evaluate the sensitivity of simulated net ecosystem production (NEP) fluxes to vegetation pool sizes. The experimental design is based on the inverse model optimization of different parameter vectors performed at the measurement site level. Each parameter vector prescribes different simulation dynamics that embody different model structural assumptions concerning (non)steady‐state conditions in vegetation and soil carbon pools. We further explore the potential of assimilating biometric constraints through the cost function for sites where in situ information on aboveground biomass or wood pools is available. The integration of biometric data yields marked improvements in the simulation of vegetation C pools compared to single constraints with eddy flux data. Overall, it is necessary to relax both vegetation and soil carbon pools for consistency with the observed data streams. Multiple constraints approaches also leads to variable model performance among the different experimental setups and model structures. We identify and assess the limitations of various model structures and the role of multiple constraints approaches for tackling issues of equifinality. These studies emphasize the need for establishing consistent data sets of fluxes and biometric data for successful model‐data fusion.