Applying single and two layer canopy models to derive conductances of a Scots pine plantation from micrometeorological measurements

Summary A single layer (Penman-Monteith) and a two layer (modified Shuttleworth-Wallace) evapotranspiration (ET) model are used alternatively to derive conductances related to the dominant fluxes of water vapor from a semi-closed Scots pine plantation. The derivations are based on micrometeorological measurements of above canopy energy flux densities and a simple resistance network. For a period of consecutive fine weather days, below canopy net radiation and below canopy ET were about 20 percent of the corresponding above canopy values. Resulting conductances for latent heat flux agreed well with porometric measurements of pines and understory scaled to canopy level. The shift from single to two layer modelling reduced the canopy conductance to pine conductance by the fraction of understory ET.However, characteristics of porometer results and micrometeorologically derived conductances were quite different: The porometer estimates of conductance were highly variable due to stomatal response to local environmental conditions or natural variability within the tree canopy and vegetation patches which characterized the forest understory. Micrometeorologically derived conductances integrate spatially resulting in relatively smooth and repetitive daily patterns that lack the information of small scale variability. This is seen as a favorable feature of micrometeorological derived conductances when used for the parameterization of atmospheric models for climate research as long as small scale bio-diversity is irrelevant.With 5 Figures     

基于碳水通量耦合原理改进Penman-Monteith蒸散发模型

陆地蒸散(ET)涵括地表和潮湿叶片的蒸发和植物的蒸散发,是陆地水循环的重要组成部分。Penman-Monteith方程是估算陆地蒸散的重要方法,方程中的叶片或冠层气孔导度是提高估算精度的关键因子。根据碳水循环的耦合原理,植物光合作用模型可用于估算叶片或冠层气孔导度。植物光合作用模型可分为三类:1)使用总冠层导度的大叶模型(BL),2)区别阴、阳叶冠层导度的双大叶模型(TBL),3)区别阴、阳叶叶片导度的双叶模型(TL)。与这三类光合作用模型相对应,衍生出基于不同导度计算方法的三种蒸散估算模型。三种蒸散模型之间的主要区别在于是否进行从叶片尺度到冠层尺度的气孔导度集成。这三种模型中,双叶模型使用叶片尺度的气孔导度,集成度最低。反之,大叶模型使用冠层尺度的气孔导度,集成度最高。由于在Penman-Monteith中,蒸腾和气孔导度之间的关系是非线性的,气孔导度的集合会导致负偏差。因此,与通量测量相比,大叶蒸散模型的估算偏差最大,而双叶蒸散模型的估算偏差最小。     

A comparison of two methods of estimating transpiration rates from a Sitka spruce plantation

Estimates of hourly transpiration from a 16–17 yr old Sitka spruce forest were calculated from the Penman-Monteith combination equation and compared with estimates from an eddy correlation/energy balance method.Canopy conductances were estimated from stomatal conductances measured using null balance diffusion porometers and took account of canopy variations of stomatal conductance and needle area index.Vertical heat fluxes were measured by the eddy correlation method; transpiration fluxes were then estimated from an energy balance of the forest.There was not a 1:1 relationship between the estimates of transpiration from the two methods. The major sources of error were concluded to be (i) difficulties of estimating the variation in stomatal conductance and leaf area through the canopy, (ii) errors in the value of total leaf area index, and (iii) errors in stomatal conductance measurements.The eddy correlation method was suggested as the more useful for future studies of the variation of forest transpiration in time or space, because the Penman-Monteith equation requires extensive biological measurements.     

盘锦湿地芦苇叶片气孔导度的模拟

基于2005年5～9月盘锦湿地芦苇叶片气体交换观测数据,针对芦苇叶片气孔导度与光合速率以及光合速率与光合有效辐射之间的关系进行分析。结果表明:芦苇叶片气孔导度与光合速率的关系可应用Ball-Berry模型描述,光合速率与光合有效辐射的关系可应用非直角双曲线光合模型描述,联合Ball-Berry模型与非直角双曲线光合模型可通过环境变量求解叶片气孔导度。模型考虑了气孔导度与光合之间的相互作用。利用实验数据对气孔导度模型验证表明,叶片气孔导度模拟值和观测值回归方程的斜率为0.95,方程决定系数R=0.82(P<0.05)。     

Introduction of a sub-grid hydrology in the ISBA land surface model

In atmospheric models, the partitioning of precipitation between infiltration and runoff has a major influence on the terrestrial water budget, and thereby on the simulated weather or climate. River routing models are now available to convert the simulated runoff into river discharge, offering a good opportunity to validate land surface models at the regional scale. However, given the low resolution of global atmospheric models, the quality of the hydrological simulations is much dependent on various processes occurring on unresolved spatial scales. This paper focuses on the parameterization of sub-grid hydrological processes within the ISBA land surface model. Five off-line simulations are performed over the French Rhône river basin, including various sets of parameterizations related to the sub-grid variability of topography, precipitation, maximum infiltration capacity and land surface properties. Parallel experiments are conducted at a high (8 km by 8 km) and low (1° by 1°) resolution, in order to test the robustness of the simulated water budget. Additional simulations are performed using the whole package of sub-grid parameterizations plus an exponential profile with depth of saturated hydraulic conductivity, in order to investigate the interaction between the vertical soil physics and the horizontal heterogeneities. All simulations are validated against a dense network of gauging measurements, after the simulated runoff is converted into discharge using the MODCOU river routing model. Generally speaking, the new version of ISBA, with both the sub-grid hydrology and the modified hydraulic conductivity, shows a better simulation of river discharge, as well as a weaker sensitivity to model resolution. The positive impact of each individual sub-grid parameterization on the simulated discharges is more obvious at the low resolution, whereas the high-resolution simulations are more sensitive to the exponential profile with depth of saturated hydraulic conductivity.     

The conductance of a maize crop and the underlying soil to ozone under various environmental conditions

Flux measurements of ozone and water vapour employing the eddy correlation technique were used to determine the surface conductance and canopy conductance to ozone. In the surface conductance to ozone, all surfaces at which ozone is destroyed and the transport process to these surfaces are included. The canopy conductance to ozone represents the ozone uptake of transpiring plant parts. The surface conductance to ozone of the maize crop and the underlying soil was generally larger than the canopy conductance to ozone. This means that beside the uptake by stomata, there was another important ozone sink. Under wet soil surface conditions, the surface conductance and the canopy conductance to ozone coincided. This indicates that the resistance of wet soil and the remaining plant parts (cuticle) to ozone was much larger than the stomatal or soil resistance. On the other hand, under dry soil conditions the conductances differ, largely caused by a variation in the transport process to the soil. The transport of ozone to soil increased with increasing friction velocity (u _*) and decreased with increasing atmospheric stability, leaf area index (LAI) or crop height (h). These effects for midday (unstable) conditions were parameterized with an in-crop aerodynamic resistance,r _inc in a very straightforward way;r _inc=13.9 LAIh/u _*+67 (cc.=0.77). If the ozone flux in air pollution models is described with a simple resistance model (Big Leaf model), the extra destruction at the soil should be modelled using an in-crop aerodynamic resistance. For these measurements the ozone flux to the soil was 0–65% of the total ozone flux measured above the crop. Under wet soil conditions, this was less than 20%; under dry soil conditions, this was 30–65%.     

Analysis of the variability of canopy resistance over a desert steppe site in Inner Mongolia,China

ABSTRACT Canopy resistance substantially affects the partitioning of available energy over vegetated surfaces. This study analyzed the variability of canopy resistance and associated driving environmental factors over a desert steppe site in Inner Mongolia, China, through the use of eddy-flux and meteorological data collected from 2008 to 2010. Distinct seasonal and interannual variabilities in canopy resistance were identified within those three years, and these variabilities were controlled primarily by precipitation. Strong interannual variability was found in vapor pressure deficit （VPD）, similar to that of air temperature. Based on the principal component regression method, the analysis of the relative contribution of five major environmental factors [soil-water content （SWC）, leaf-area index （LAI）, photosynthetically active radiation （Kp）, VPD, and air temperature] to canopy resistance showed that the canopy-resistance variation was most responsive to SWC （with 〉 35% contribution）, followed by LAI, especially for water-stressed soil conditions （〉 20% influence）, and VPD （consistently with an influence of approximately 20%）. Canopy-resistance variations did not respond to Kp due to the small interannual variability in Kp during the three years. These analyses were used to develop a new exponential function of water stress for the widely used Jarvis scheme, which substantially improved the calculation of canopy resistance and latent heat fluxes, especially for moist and wet soils, and effectively reduced the high bias in evaporation estimated by the original Jarvis scheme. This study highlighted the important control of canopy resistance on plant evaporation and growth for the investigated desert steppe site with a relatively low LA1.     

Implementation of a Stable PBL Turbulence Parameterization for the Mesoscale Model MM5: Nocturnal Flow in Complex Terrain

The difficulties associated with the parameterization of turbulence in the stable nocturnal planetary boundary layer (PBL) have been a great challenge for the nighttime predictions from mesoscale meteorological models such as MM5. As such, there is a general consensus on the need for better stable boundary-layer parameterizations. To this end, two new turbulence parameterizations based on the measurements of the Vertical Transport and Mixing (VTMX) field campaign were implemented and evaluated in MM5. A unique aspect of this parameterization is the use of a stability-dependent turbulent Prandtl number that allows momentum to be transported by the internal waves, while heat diffusion is impeded by the stratification. This improvement alleviates the problem of over-prediction of heat diffusion under stable conditions, which is a characteristic of conventional atmospheric boundary-layer schemes, such as the Medium Range Forecast (MRF) and Blackadar schemes employed in MM5. The predictions made with the new PBL scheme for the complex terrain airshed of Salt Lake City were compared with those made with a default scheme of MM5, and with observations made during the VTMX campaign. The new schemes showed an improvement in predictions, particularly for the nocturnal near-surface temperature. Surface wind predictions also improved slightly, but not to the extent of temperature predictions. The default MRF scheme showed a significantly higher surface temperature than observed, which could be attributed to the enhanced vertical heat exchange brought about by its turbulence parameterization. The modified parameterizations reduced the surface sensible heat flux, thus enhancing the strength of the near-surface inversion and lowering the temperature towards the observed values.     

Estimation of Regional Evapotranspiration by Combining Aircraftand Ground-Based Measurements

Two methods are examined for combining measurements from instrumented aircraftand towers to estimate regional scale evapotranspiration. Aircraft data provided spatially averaged values of properties of the surface, the evaporative fraction and maximum stomatal conductance. These quantities are less sensitive to meteorological conditions than the turbulent fluxes of heat and water vapour themselves. The methods allowed aircraft data collected over several days to be averaged and thus to reduce the random error associated with the temporal under-sampling inherent in aircraft measurements. Evaporative fraction is estimated directly from the aircraft data, while maximum stomatal conductance is estimated by coupling the Penman–Monteith equation to a simple model relating surface conductance to the incoming shortwave radiation and specific humidity saturation deficit. The spatial averages of evaporative fraction and maximum stomatal conductance can then be used with routine tower data to estimate the regional scale evapotranspiration. Data from aircraft flights and six ground based sites during the OASIS field campaign in south–east New South Wales in 1995 have been used to check the methods. Both the evaporative fraction and the maximum stomatal conductance derived from the aircraft data give information on the spatial variability of the surface energy budget at scales from 10 to 100 km. Daily averaged latent heat fluxes estimated using these methods for the OASIS study region agree with the available observations in quasi-stationary conditions or in weakly non-stationary conditions when the data from several aircraft flights are averaged to reduce the impact of short term imbalances in the surface energy budget.     

Use of atmospheric radiation measurement program data from Barrow, Alaska, for evaluation and development of snow-albedo parameterizations

Summary Snow albedo is determined from the ratio of out-going to incoming solar radiation using three years of broadband shortwave radiometer data obtained from the Barrow, Alaska, Atmospheric Radiation Measurement (ARM) site. These data are used for the evaluation of various types of snow-albedo parameterizations applied in numerical weather prediction or climate models. These snow-albedo parameterizations are based on environmental conditions (e.g., air or snow temperature), snow related characteristics (e.g., snow depth, snow age), or combinations of both. The ARM data proved to be well suited for snow-albedo evaluation purposes for a low-precipitation tundra environment. The evaluation confirms that snow-age dependent parameterizations of snow albedo work well during snowmelt, while parameterizations considering meteorological conditions often perform better during snow accumulation. Current difficulties in parameterizing snow albedo occur for long episodes of snow-event free conditions and episodes with a high frequency of snow events or strong snowfall. In a further step, the first two years of the ARM albedo dataset is used to develop a snow-albedo parameterization, and the third year’s data serves for its evaluation. This parameterization considers snow depth, wind speed, and air temperature which are found to be significant parameters for snow-albedo modeling under various conditions. Comparison of all evaluated snow-albedo parameterizations with this new parameterization shows improved snow-albedo prediction. Correspondence: Nicole M?lders, Geophysical Institute and College of Natural Science and Mathematics, University of Alaska Fairbanks, 903 Koyukuk Drive, P.O. Box 757320, Fairbanks, AK 99775-7320, USA     

A history of mesoscale model development

The development of atmospheric mesoscale models from their early origins in the 1970’s until the present day is described. Evolution has occurred in dynamical and physics representations in these models. The dynamics has had to change from hydrostatic to fully nonhydrostatic equations to handle the finer scales that have become possible in the last few decades with advancing computer power, which has enabled real-time forecasting to go to finer grid sizes. Meanwhile the physics has also become more sophisticated than the initial representations of the major processes associated with the surface, boundary layer, radiation, clouds and convection. As resolutions have become finer, mesoscale models have had to change paradigms associated with assumptions related to what is considered sub-grid scale needing parameterization, and what is resolved well enough to be explicitly handled by the dynamics. This first occurred with cumulus parameterization as real-time forecast models became able to represent individual updrafts, and is now starting to occur in the boundary layer as future forecast models may be able resolve individual thermals. Beyond that, scientific research has provided a greater understanding of detailed microphysical and land-surface processes that are important to aspects of weather prediction, and these parameterizations have been developing complexity at a steady rate. This paper can just give a perspective of these developments in the broad field of research associated with mesoscale atmospheric model development.     

Factors affecting the canopy resistance of a Douglas-fir forest

The physiological nature of canopy resistance was studied by comparing the stomatal and canopy resistance of a 10-m high Douglas-fir forest. Stomatal resistance of the needles was measured using porometry, while the canopy resistance was calculated using energy balance/Bowen ratio measurements of evapotranspiration. A typical steady increase in the forest canopy resistance during daytime hours, even at high soil water potentials, was observed. A similar trend in the stomatal resistance indicated that increasing canopy resistance during the daytime was caused by gradually closing stomata. During a dry period when soil water potentials declined from 0 to –10.5 bars, the mean daytime value of canopy resistance increased in proportion to the mean daytime value of the stomatal resistance. Values of canopy resistance calculated from stomatal resistance and leaf area index measurements agreed well with those calculated from energy balance measurements. The dependences of stomatal resistance on light, vapour pressure deficit, twig and soil water potentials art summarized.     

A Dynamic Urban Canopy Parameterization for Mesoscale Models Based on Computational Fluid Dynamics Reynolds-Averaged Navier–Stokes Microscale Simulations

Urban canopy parameterizations (UCPs) are necessary in mesoscale modelling to take into account the effects of buildings on wind and turbulent structures. This study is focused on the dynamical part of UCPs. The main objective is twofold: first, computing important UCP input parameters (turbulent length scales and the sectional drag coefficient) by means of Reynolds-averaged Navier–Stokes (RANS) simulations of turbulent flow; and second, comparing UCP variables with spatially-averaged properties obtained from RANS simulations for the same configurations. The results show the importance of using a suitable parameterization of the drag force for different packing densities. An urban canopy parameterization that is a compromise between simplicity and accuracy is proposed. This scheme accounts for the variation of drag coefficients with packing densities, and has a parameterization of turbulent length scales. The technique adopted ensures that, at least for the simple configurations studied, the urban canopy parameterization gives values of spatially-averaged variables similar to those computed from a more complex simulation, such as RANS that resolves explicitly the flow around buildings.     

Energy balance in urban Mexico City: observation and parameterization during the MILAGRO/MCMA-2006 field campaign

The parameterization of the energy balance from a residential and commercial neighborhood of Mexico City was investigated using direct measurements of radiative and heat fluxes carried out during the MILAGRO/MCMA-2006 field campaign as a reference. The measured fluxes were used to evaluate different models of the energy balance based on parameterizations that require standard meteorological observations: ambient temperature, relative humidity, atmospheric pressure and cloudiness. It was found that these models reproduce with reasonable accuracy the diurnal features of the radiative and heat fluxes. The largest differences between modeled and observed fluxes correspond to the incoming longwave radiation, mainly due to errors in the cloudiness data. This paper contributes to the understanding of the energy partitioning in (sub)tropical urban environments, particularly in the developing world, where energy balance models have not been evaluated.     

利用遥感信息研究区域冬小麦气孔导度的时空分布

气孔导度是影响作物蒸散和作物的光合速率进而影响作物产量的重要因子。文中通过利用NOAA－AVHRR数据首次对华北平原典型区冬小麦气孔导度分布进行了研究,给出了华北平原典型区冬小麦不同生长季节的气孔导度空间分布状况,为进一步研究田间水分和作物蒸散对产量影响以及建立遥感作物水分胁迫生物量模型和监测不同生育期的农田缺水等提供依据。     

Intra-field variability of scalar flux densities across a transition between a desert and an irrigated potato field

This paper reports on measurements of sensible and latent heat and CO₂ fluxes made over an irrigated potato field, growing next to a patch of desert. The study was conducted using two eddy correlation systems. One measurement system was located within the equilibrium boundary layer 800 m downwind from the edge of the potato field. The other measurement system was mobile and was placed at various downwind positions to probe the horizontal transition of vertical scalar fluxes. Latent (LE) and sensible (H) heat fluxes, measured at 4 m above the surface, exhibited marked variations with downwind distance over the field. Only after the fetch to height ratio exceeded 75 to 1 didLE andH become invariant with downwind distance. When latent and sensible heat fluxes were measured upwind of this threshold, significant advection of humidity-deficit occurred, causing a vertical flux divergence ofH andLE.The measured fluxes of momentum, heat, and moisture were compared with predictions from a second-order closure two-dimensional atmospheric boundary layer model. There is good agreement between measurements and model predictions. A soil-plant-atmosphere model was used to examine nonlinear feedbacks between humidity-deficits, stomatal conductance and evaporation. Data interpretation with this model revealed that the advection of hot dry air did not enhance surface evaporation rates near the upwind edge of the potato field, because of negative feedbacks among stomatal conductance, humidity-deficits, andLE. This finding is consistent with results from several recent studies.     

Regional climate modelling of the 2006 West African monsoon: sensitivity to convection and planetary boundary layer parameterisation using WRF

Regional climate model (RCM) is a valuable scientific tool to address the mechanisms of regional atmospheric systems such as the West African monsoon (WAM). This study aims to improve our understanding of the impact of some physical schemes of RCM on the WAM representation. The weather research and forecasting model has been used by performing six simulations of the 2006 summer WAM season. These simulations use all combinations of three convective parameterization schemes (CPSs) and two planetary boundary layer schemes (PBLSs). By comparing the simulations to a large set of observations and analysis products, we have evaluated the ability of these RCM parameterizations to reproduce different aspects of the regional atmospheric circulation of the WAM. This study focuses in particular on the WAM onset and the rainfall variability simulated over this domain. According to the different parameterizations tested, the PBLSs seem to have the strongest effect on temperature, humidity vertical distribution and rainfall amount. On the other hand, dynamics and precipitation variability are strongly influenced by CPSs. In particular, the Mellor?CYamada?CJanjic PBLS attributes more realistic values of humidity and temperature. Combined with the Kain?CFritsch CPS, the WAM onset is well represented. The different schemes combination tested also reveal the role of different regional climate features on WAM dynamics, namely the low level circulation, the land?Catmosphere interactions and the meridional temperature gradient between the Guinean coast and the Sahel.     

城市冠层过程的研究与进展

介绍了城市冠层过程研究的概况和历史,概括了城市冠层的特征及其作用于局地气候和中尺度系统的主要方式,总结了城市冠层参数化方案或模式的研究进展.比较了几种参数化方法和模式的特点,提出目前存在的一些问题.认为：现有的参数化方案和模式尚未能完整体现城市下垫面特征、准确反映人类活动影响,二者的发展有赖于对城市下垫面建筑特征更合理、细致地描述,对人为热量、水汽影响更准确的估计与刻划.与局地因素密切相关的城市冠层模式还需要更多地在模式中引入局地差异影响.同时指出,要提高耦合城市过程的中尺度模式预报水平,需要相应改进包括云、降水、次网格地形及边界层动力学等的参数化.     

MM5模式中城市冠层参数化方案的设计及其数值试验

文中在综合国外一些较先进的中尺度模式城市作用参数化方案的基础上 ,从城市下垫面结构对城市边界层大气作用的物理机制及实际应用两方面出发 ,对城市下垫面结构和人为活动等因素对边界层结构的影响及中尺度模式中城市化作用的合理体现等问题进行了较全面的考虑 ,改进和设计出能够较全面、细致地描述城市结构对大气边界层动力、热力结构的影响 ,且适合中尺度模式结构特点的城市冠层参数化方案 (UCP) ,并实现了其与MM5模式的耦合。进行了耦合后的UCP方案及采用原城市作用方案的MM5模式对BECAPEX试验期间北京地区气象条件多重嵌套细尺度进行了模拟试验 ,并与观测结果对比 ,结果表明 :相比于MM 5模式中原有表示城市作用的参数化方案来讲 ,设计的UCP方案在很大程度上提高了MM 5模式对城市边界层热力和动力结构的模拟能力。     

Global and regional coupled climate sensitivity to the parameterization of rainfall interception

A coupled land?Catmosphere model is used to explore the impact of seven commonly used canopy rainfall interception schemes on the simulated climate. Multiple 30-year simulations are conducted for each of the seven methods and results are analyzed in terms of the mean climatology and the probability density functions (PDFs) of key variables based on daily data. Results show that the method used for canopy interception strongly affects how rainfall is partitioned between canopy evaporation and throughfall. However, the impact on total evaporation is much smaller, and the impact on rainfall and air temperature is negligible. Similarly, the PDFs of canopy evaporation and transpiration for six selected regions are strongly affected by the method used for canopy interception, but the impact on total evaporation, temperature and precipitation is negligible. Our results show that the parameterization of rainfall interception is important to the surface hydrometeorology, but the seven interception parameterizations examined here do not cause a statistically significant impact on the climate of the coupled model. We suggest that broad scale climatological differences between coupled climate models are not likely the result of how interception is parameterized. This conclusion is inconsistent with inferences derived from earlier uncoupled simulations, or simulations using very simplified climate models.