Parameter estimation in surface exchange models using nonlinear inversion: how many parameters can we estimate and which measurements are most useful?

Models of mass and energy exchanges between the biosphere and the atmosphere generally contain a nonlinear dependence between fluxes and model parameters, and thus estimation of these parameters from measurements in a heterogeneous landscape depends on the scale of the observations. The scale‐dependence of a typical surface‐exchange model (the CSIRO Biospheric Model, CBM) is examined using the diurnal variation of hourly fluxes of CO₂, latent heat, sensible heat and soil heat. The fluxes were measured using micrometeorological techniques over six sites in a grazing/pasture system in SE Australia during a period of three weeks in 1995. Nonlinear parameter inversion was used to determine model parameters. Analysis of the covariance of the estimates of the parameters and the unexplained residuals of the model showed that a maximum of three or four parameters could be determined independently from the observations for all six sites. Estimates of a key model parameter, j_max, the mean of maximum potential electron transport rate of all leaves within the canopy, was best determined by the measurements of net CO₂ flux at all sites examined. Measurements of ground heat flux provide little information about any of the model parameters in CBM. Because of nonlinearities in the surface exchange model, calculated fluxes will be in error if parameters for the component vegetation types are simply averaged in proportion to their areal fraction. The magnitude of these errors was examined for CBM using a hypothetical land surface consisting of two surface types, each with different parameter values. Predictions of net CO₂, latent heat and ground heat fluxes using a linear combination of model parameters for the two surface types were quite similar with those found using optimal estimates of the parameters for the landscape, but were significantly poorer for sensible heat fluxes.     

季节性干旱对中亚热带人工林显热和潜热通量日变化的影响

植被与大气间的显热和潜热通量的日变化是大气过程和植被生理过程的显著标志。本研究利用ChinaFLUX千烟洲站典型的夏季雨热不同季的季节性干旱的试验条件,探讨了2003年季节性干旱对该生态系统显热和潜热通量日变化变异幅度和峰值时间的影响。研究表明：显热通量的日变化变异幅度年平均值为176 W/m2。潜热通量的日变化变异幅度年平均值为171 W/m2。显热通量到达日变化峰值的时间平均为11:57。全年潜热通量的日变化都在午后达到峰值,平均值为12:33。季节性干旱造成显热通量的日变异幅度明显增大,从144W m-2增加到321 W m-2。而潜热通量的日变异幅度明显降低,从324 W/m2减小到198 W/m2。,显热和潜热通量日变异幅度的相对变化明显增大,从-165 W/m2增加到76 W/m2,气温和饱和水汽压差是影响显热和显热日变异幅度及其相对变化的主要控制因素。干旱胁迫期,深层水对显热通量日变化变异幅度及其与潜热通量日变化变异幅度的相对变化的作用更显著,而潜热通量日变化变异幅度与气象要素关系不显著。季节性干旱造成显热通量日变化的峰值时间和显热和潜热通量日变化峰值时间的相对变化明显向下午偏移,显热通量日变化的峰值从上午11:31到中午12:17,相对变化从1小时到1小时20分钟。季节性干旱对潜热通量日变化峰值时间没有显著的影响。非干旱胁迫期,显热通量日变化峰值时间和显热及潜热通量日变化峰值时间的相对变化均与气温负相关,而干旱胁迫期,则与气温正相关。潜热通量日变化峰值时间与气象要素关系均不显著。该生态系统显热和潜热通量日变化峰值的相对变化主要受降水量的季节分配控制,在干旱胁迫期降水的作用更加明显。潜热和显热通量日变化峰值时间的相对变化总体上都受植被与大气间的耦合程度控制。     

玉米农田水热通量动态与能量闭合分析

 基于锦州农田生态系统野外观测站玉米农田涡度相关系统近2年的水热通量观测数据，分析了玉米农田水热通量的日际、年际变化特征及其能量 平衡状况。结果表明： 1）玉米农田水热通量日变化与年变化均呈单峰型二次曲线，峰值出现在12∶00～13∶00左右，与净辐 射的日变化、年 变化同步，潜热通量最大可达到655 w•m^－2（出现在2004年7月8日1 3∶00），显热通量最大值大约为369 w•m^－2（出现在2004年5月31日13∶ 00）。2）玉米农田水热通量强度与局地的环境条件密切相关：显热通量与大气压的年变化呈负相关，潜热通量与气温年变化呈正相关。水热通 量受降水的影响较大，对降水的反应较敏感。其中，潜热通量(LE)不仅与降水的强度有关，而且随着降水的季节分布的不同而出现不同的响应 ，即使同样量级的降水在夜间与白天对LE的影响也是不同的。3）玉米农田通量观测呈现能量不闭合现象，主要原因可能是未包含0～5 cm土壤 热储量与冠层热储量，造成大约15.5%的能量损失。     

开放式空气C02浓度增高对水稻冠层能量平衡的影响

大气CO2浓度升高对植物冠层能量平衡的影响是导致植物生长发育和水分利用率发生变化的环境物理原因．利用位于江苏省无锡市安镇的农田自由开放式空气CO2浓度增高(FACE)系统平台,进行水稻冠层微气候和土壤热通量的连续观测,并结合能量平衡分析,研究了FACE对水稻冠层能量平衡的影响．结果表明,水稻冠层显热和潜热通量FACE与对照的差异日最大值出现在14：00左右,与空气相对湿度日最低值出现时间一致;潜热通量FACE与对照的差异日最大值变化在—15-—65J·m^-2·s^-1之间,显热通量FACE与对照的差异最低值变化在12—55J·m^-2·s^-1之间;显热和潜热通量FACE与对照的差异日最大值随冠层上方辐射平衡增加而增大．水稻冠层白天总显热通量FACE均高于对照,而总潜热通量FACE均低于对照．白天总显热和潜热通量FACE与对照的差异在同一生育期内随冠层上方净辐射增强而增大,在不同生育期随生育期推进而减少．开花期至蜡熟期,水稻冠层白天总潜热通量FACE比对照平均低6．7％．FACE使水稻冠层白天总显热通量及其占冠层上方辐射平衡的比例减少,而使总潜热通量及其占冠层上方辐射平衡的比例增大,但对土壤热通量及夜间显热和潜热通量的影响不大．开花期至蜡熟期水稻冠层白天总显热、潜热通量占冠层上方净辐射总量的比例FACE与对照之差平均为5．5％．     

Coupling heat flux dynamics with meteorological conditions in the green roof ecosystem

Green roofs can notably modify the thermal properties of the building envelope and adjacent air to bring environmental benefits. This study investigates the heat flux dynamics of the tropical green roof ecosystem to provide a scientific basis for design and management. Green roof experimental plots were established to monitor the total solar radiation, net radiation, and micrometeorological parameters. The data permit calculation of sensible and latent heat fluxes using the Bowen ratio energy balance (BREB) method. The results demonstrated the life cycle characteristics of heat flux components. The dynamic changes of sensible (H), latent (λE) and soil (G) heat fluxes were denoted by single-peak quadratic curves. Net radiation (R_n) was largely determined by quantity and variation trends of λE, reaching at 1300 h a maximum λE of 655 W m^?2 and maximum H of 369 W m^?2. Temporal heat-flux fluctuations were strongly correlated with meteorological variables. Extreme values of H and λE correlated well with precipitation and temperature (R² = 0.78). Dynamics of heat-flux magnitude and partitioning demonstrated notable differences by daily and season periods. They displayed considerable variations in flux partitioning, with Bowen ratios strongly correlated with weather conditions and vegetation types. The energy budget of the green roof ecosystem is unbalanced with a heat loss of about 15.5% caused by soil and canopy heat reserve. The passive indoor cooling effect under the green roof is attributed to the unbalanced energy closure.     

Interannual and seasonal variations of energy and water vapor fluxes above a tropical seasonal rain forest in Xishuangbanna, SW China

The values and variation characteristics of energy components, their relationship with net radiation and the characteristics of water balance in the forest were analyzed, based on the observation data of energy fluxes, meteorological parameters and biomass in a tropical seasonal rain forest in Xishuangbanna from January 2003 to December 2004. The results show that annual net radiation was 3516.4 MJ/(m² · a) and 3516.6 MJ/(m² · a) in 2003 and 2004, respectively, of which 46% and 44% were used in latent heat flux, and 12% and 11% were lost as sensible heat flux. Annual mean canopy surface conductance was 10.3 mm/s and 10.0 mm/s in 2003 and 2004, respectively. Moreover, canopy surface conductance was lower in dry-hot seasons than in fog-cool and rainy seasons. Canopy surface conductance correlated significantly and positively with leaf area index, but negatively with water vapor pressure deficit. In general, canopy surface conductance was not affected directly by soil water content, but highly depended on soil moisture status when soil water content was below 0.15 m³/m³. Annual total evapotranspiration of this forest ecosystem in dry seasons was lower than that in rainy seasons, which was considered as one of the most important reasons that tropical seasonal rain forest could survive and flourish in Xishuangbanna at limit of water and heat.     

Measuring and modelling the radiation balance of a heterogeneous shrubland

The accurate quantification of the energy available for sensible and latent heat transfer from plant canopies is essential for the prediction of impacts of climate on vegetation water use and growth. Unlike agricultural fields and extensive forests of more humid zones, vegetation growing in semi‐arid climates is usually sparse creating a heterogeneous surface of shrubs, annuals and bare soil. Under these conditions many of the assumptions of the basic equations used in microclimatology, which assume a uniform vegetated surface, may be violated. It is proposed here that heterogeneous canopies may require a formulation of their energy balance that includes a measure of the canopy complexity in order to both interpret field measurements and be used in predictive models. This paper explores the need for a more complex formulation of the vegetation energy balance through a series of experiments on a sparse clumped shrubland of Retama sphaerocarpa in the Tabernas Desert, Almería, south‐east Spain. These experiments investigated the importance of the radiative properties of each surface on energy balance of soil, annuals and shrubs individually, and the surface as a whole. The study evaluated the use of the fractional vegetative cover (f) and the radiative characteristics of each surface (reflection coefficients and emissivities) for calculating net radiation partitioning between shrubs and bare soil. Results indicated that partitioning of net radiation between components could be accurately calculated from values of fractional vegetative cover, reflection coefficients and emissivities for both bare soil and plant surfaces. A sensitivity analysis showed the importance of specific radiation properties of each surface. Measurements of horizontal long‐wave fluxes between components showed that the location of a plant with respect to other plants made little difference to its long‐wave energy balance. The results also emphasized the importance of soil water content on the energy balance, through its influence on albedo and soil heat storage. This was particularly true when measurements of soil heat flux were used to measure the available energy for soil under shrubs and bare soil because of strong hysteresis cycles. These cycles were larger in the bare soil than in the substrate under shrubs.     

Structural complexity and land‐surface energy exchange along a gradient from arctic tundra to boreal forest

Question: Current climate changes in the Alaskan Arctic, which are characterized by increases in temperature and length of growing season, could alter vegetation structure, especially through increases in shrub cover or the movement of treeline. These changes in vegetation structure have consequences for the climate system. What is the relationship between structural complexity and partitioning of surface energy along a gradient from tundra through shrub tundra to closed canopy forest? Location: Arctic tundra‐boreal forest transition in the Alaskan Arctic. Methods: Along this gradient of increasing canopy complexity, we measured key vegetation characteristics, including community composition, biomass, cover, height, leaf area index and stem area index. We relate these vegetation characteristics to albedo and the partitioning of net radiation into ground, latent, and sensible heating fluxes. Results: Canopy complexity increased along the sequence from tundra to forest due to the addition of new plant functional types. This led to non‐linear changes in biomass, cover, and height in the understory. The increased canopy complexity resulted in reduced ground heat fluxes, relatively conserved latent heat fluxes and increased sensible heat fluxes. The localized warming associated with increased sensible heating over more complex canopies may amplify regional warming, causing further vegetation change in the Alaskan Arctic.     

Evapotranspiration under advective conditions

Arid and semi-arid regions are heterogeneous landscapes in which irrigated fields are surrounded by arid areas. The advection of sensible heat flux from dry surfaces is a significant source of energy that has to be taken into consideration when evaluating the evaporation from crops growing in these areas. The basic requirement of most of the common methods for estimating evapotranspiration [Bowen ratio, aerodynamic and Penman–Monteith (PM) equation] is that the horizontal fluxes of sensible and latent heat are negligible when compared to the corresponding vertical fluxes. We carried out measurements above an irrigated tomato field in a desert area. Latent and sensible heat fluxes were measured using a four-level Bowen machine with aspirated psychrometers. Our results indicate that under advective conditions only measurements carried out in the lowest layer are satisfactory for the estimation of latent heat fluxes and that the use of the PM equation with an appropriately parameterized canopy resistance may be preferable.This work was carried out by the author during her sojourn at the Wyler Department for Drayland Agriculture, Jacob Blaustein Institute for Desert Research, Ben Gurion University of the Negev, Israel.Manufacturers names are provided for the benefit of the reader and do not imply endorsement by the authors.     

极端干旱区大气边界层厚度时间演变及其与地表能量平衡的关系

主要采用ECMWF的地表和大气产品分析了中国西北极端干旱区大气边界层厚度与地表能量通量的时间变化特征,同时,结合探空加强观测分析了大气边界层演变的可能因素.得出:西北极端干旱区大气边界层厚度呈现出季节性的年际和年代际变化,夏季大气边界层厚度呈下降趋势,春、秋季节呈现出先增加后降低的趋势,冬季以阶段性降低趋势为主,20世纪80年代是大气边界层厚度的转折时期;感热通量是极端干旱区大气边界层发展的主要热力因素;由于夏季净辐射量、地气温差、粗糙度以及风速等因子随时间演变而呈降低趋势,潜热通量呈增加趋势,导致了边界层高度形成的热力作用减弱,边界层厚度降低;同时,粗糙度和风速也是大气边界层发展的主要动力因素,由于边界层粗糙度和风速降低,促使垂直风切变减小,湍流动力作用减弱,也会导致边界层厚度降低.     

Remote sensing of forest gas exchange: Considerations derived from a tomographic perspective

The global exchange of gas (CO₂, H₂O) and energy (sensible and latent heat) between forest ecosystems and the atmosphere is often assessed using remote sensing (RS) products. Although these products are essential in quantifying the spatial variability of forest–atmosphere exchanges, large uncertainties remain from a measurement bias towards top of canopy fluxes since optical RS data are not sensitive for the vertically integrated forest canopy. We hypothesize that a tomographic perspective opens new pathways to advance upscaling gas exchange processes from leaf to forest stands and larger scales. We suggest a 3D modelling environment comprising principles of ecohydrology and radiative transfer modelling with measurements of micrometeorological variables, leaf optical properties and forest structure, and assess 3D fields of net CO₂ assimilation (A_n) and transpiration (T) in a Swiss temperate forest canopy. 3D simulations were used to quantify uncertainties in gas exchange estimates inherent to RS approaches and model assumptions (i.e. a big‐leaf approximation in modelling approaches). Our results reveal substantial 3D heterogeneity of forest gas exchange with top of canopy A_n and T being reduced by up to 98% at the bottom of the canopy. We show that a simplified use of RS causes uncertainties in estimated vertical gas exchange of up to 300% and that the spatial variation of gas exchange in the footprint of flux towers can exceed diurnal dynamics. We also demonstrate that big‐leaf assumptions can cause uncertainties up to a factor of 10 for estimates of A_n and T. Concluding, we acknowledge the large potential of 3D assessments of gas exchange to unravelling the role of vertical variability and canopy structure in regulating forest–atmosphere gas and energy exchange. Such information allows to systematically link canopy with global scale controls on forest functioning and eventually enables advanced understanding of forest responses to environmental change.     

基于改进SW模型的千烟洲人工林蒸散组分拆分及其特征

蒸散组分拆分是准确评估陆地生态系统生产力以及估算水分利用效率的重要基础。利用改进后的Shuttleworth-Wallace模型,将蒸散拆分为植被蒸腾、土壤蒸发和冠层截留蒸发,并采用Monte Carlo随机参数化方案对模型参数进行优化。将模型与千烟洲亚热带人工针叶林站点的2011年涡度相关及小气候观测资料结合,对千烟洲人工林蒸散及其组分进行模拟。研究结果表明:半小时尺度上蒸散量模拟值与实测值的一致性在晴天和雨天都较高。半小时尺度上全年蒸散模拟值与实测值的决定系数、均方根误差和平均偏差为0.73、1.55 mmol m~(-2)s~(-1)和0.21 mmol m~(-2)s~(-1)。蒸散是该生态系统水分输出的最主要贡献项,占全年降水的80%。在蒸散中,植被蒸腾约占总蒸散量的85%,可推测2011年千烟洲人工林生态系统有较高的水分利用效率。该生态系统的蒸腾量季节变化明显,主要受饱和水汽压差和气温两种环境因素以及植被的叶面积指数影响且与三者均呈正相关;土壤蒸发约占总蒸散量的5%,季节变化平缓;模拟的冠层截留蒸发量约占总蒸散量的10%,季节变化大,与降水量呈正相关,与暴雨频次呈负相关,说明冠层无法有效截留强降水。该模型参数较少、时间分辨率高且可以有效模拟蒸散及其组分特征,是陆地生态系统水分循环过程研究有力的模型工具。     

土壤-植物-大气连续体水热动态模拟的研究

本文从能量平衡原理和质量守恒定律出发,描述了土壤-植物-大气连续体中的热量转换和水分输送,模拟了系统中水分和热量的动态变化过程,并用所建立的模拟模型计算了冬小麦群落的冠层温度、叶水势及系统的潜热与显热变化关系,结果表明该模型有一定的可靠性。     

无水分胁迫下行作物蒸发散与双涌源能量分配和交换关系

以内蒙古浑善达克沙地人工草地种植的行作物——青贮玉米为研究对象,将FAO-56的双作物系数法与双涌源能量平衡模型相结合,计算了太阳入射能量按叶面积指数（LAI）分配到两个涌源（冠层、土壤表面）的有效能量A_c和A_s、潜热通量λE_c和λE_s以及显热通量H_c和H_s.分析两个涌源在有效能量驱使下的潜热和显热通量相互作用.结果表明,1）在无水分胁迫条件下,冠层H_c与λE_c相互作用使冠层吸收微热平流,强化蒸腾作用,加大蒸腾量.蒸腾（以潜热通量表示）超过冠层有效能量的增量(λE_cⁱ-A_cⁱ).最大值出现在生长发育阶段LAI为0.6的7月15日到LAI为2.4 的8月9日之间,其平均值为4.32 MJ·m^-2·d^-1.2）无水分胁迫情况下,λE_s和H_s相互作用,除强湿润过程后的1~2 d外,其他各天土壤均以低于土壤表面有效能量的速率蒸发.蒸发强度取决于土壤表面有效能量消散为土壤潜热通量的百分比,这个百分比最小值出现在生长中期阶段,其平均值为11.5%;最大值出现在生长初始阶段,其平均值为51.9%.3）两个涌源潜热通量是蒸散过程中能量交换的主要成分,在生长发育、中期、后期阶段转换为两个涌源潜热通量的有效能量均占总能量的83%以上.     

我国亚热带毛竹林生长季能量通量过程及闭合度分析

为探讨我国亚热带毛竹林(Phyllostachys edulis)生长季的能量平衡关系,利用开路涡度相关法,对2011年毛竹林生长季的能量通量的变化特征进行了研究,并应用能量平衡比率法和线性回归2种方法,分析了能量闭合的特点。结果表明,我国亚热带毛竹林生长季的净辐射总量为1738.2 MJ m–2,显热通量为354.3 MJ m–2,潜热通量为1146.0 MJ m–2,土壤热通量为58.9 MJ m–2,土壤为热汇,显热通量占净辐射的20.4%,潜热通量占65.9%,土壤热通量占3.4%。毛竹林生长季的能量闭合度为0.89,月平均闭合度为0.91,但仍有11%的能量不闭合。可见,毛竹林生长季以潜热能量散失形式为主,各能量分量均以净辐射变化为基础,且日变化基本呈单峰型曲线。     

Evaluation of actual evapotranspiration of a Quercus ilex L. stand by the Bowen Ratio-Energy Budget method

Actual evapotranspiration from a closed-canopy Quercus ilex L. stand has been estimated by applying the Bowen Ratio-Energy Budget method. Daily water loss was 3.5 mm day^–1, with a peak rate near 0.6 mm hour^–1. The phenomenon of thermal inversion, quite common in mediterranean climates, seemed to play a significant role in reducing evapotranspiration, by promoting dew formation and delaying the establishment of fluxes of latent and sensible heat away from the canopy. Dew, which may form over many hours in the night, appears to be a major sink of available energy in the early morning and may represent a useful water source for stressed foliage. The alternating processes of condensation and evaporation may have a beneficial effect on the closed stand micro-environment.     

LAS测算森林冠层上方温度结构参数的可行性

利用位于河南省济源市的华北低丘山地30年生栓皮栎-侧柏-刺槐人工混交林2010年4月至8月每月连续7d,LAS直接测算的森林冠层上方湍流结构参数,与经过湍流谱方法计算处理的三维超声风速/温度仪的观测数据比较,分析LAS测算低丘山地森林冠层温度湍流结构的可行性。结果表明:水平风速和温度湍流谱都有明显的惯性区出现(斜率-2/3);LAS直接测算的湍流温度结构参数与利用该惯性区的数据计算的结果具有较好的一致性,说明在起伏非均匀下垫面上,采用LAS观测湍流结构的变化情况具有较好的可行性。     

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

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

Energy exchange within individual layers of a meadow

Summary Measurements of the radiation extinction in a meadow at Baumkirchen (Tyrol) show that the decrease in the photosynthetically active radiation (PhAR: 400–700 nm) is different to that of the total net-radiation in a characteristic way. The photosynthetically active radiation is distributed evenly to all vegetation layers, the active surface comprises practically the entire 90 cm high canopy. The total radiation energy (net-radiation) is absorbed and 45% is converted into sensible and latent heat only in a 25 cm wide layer, i.e., between 30 and 55 cm within the canopy. A second active surface lies in the lowermost 10 cm of the meadow and at the soil surface where additional 28% of the radiation energy is transformed into heat.     

内蒙古典型草原小叶锦鸡儿灌丛化对水分再分配和利用的影响

灌丛化是全球草原地区存在的主要环境问题。通过对内蒙古典型草原区小叶锦鸡儿灌丛和草地斑块冠层降雨再分配、地表径流、土壤含水量的对比观测,研究了小叶锦鸡儿灌丛化对该区水分再分配和利用的影响。结果表明,灌丛和草地斑块的冠层截留量分别占降雨量的20.86%和7.88%,灌丛和草地斑块的平均地表径流系数分别为5.95%和17.19%。土壤含水量观测结果显示,0—60 cm土层中,降雨事件过程中,灌丛斑块较草地斑块能捕获更多水分,灌丛斑块植被冠层下方土壤含水量高于草地斑块;而在雨后无有效降水补充土壤水分的前提下,0—60 cm土层中,灌丛斑块土壤水分蒸散发量高于草地斑块,其中0—10cm土层中灌丛斑块土壤水分蒸散发速率低于草地斑块,10—60 cm土层中灌丛斑块土壤水分蒸散发速率高于草地斑块。研究认为,在水分为关键性限制因子的干旱半干旱区,小叶锦鸡儿灌丛化过程增加草原生态系统中水分分布的空间异质性,灌丛斑块能捕获、利用更多水分以维持更多的生物量。