A Comparison of Two Canopy Radiative Models in Land Surface Processes

This paper compares the predictions by two radiative transfer models-the two-stream approximation model and the generalized layered model （developed by the authors） in land surface processes -for different canopies under direct or diffuse radiation conditions. The comparison indicates that there are significant differences between the two models, especially in the near infrared （NIR） band. Results of canopy reflectance from the two-stream model are larger than those from the generalized model. However, results of canopy absorptance from the two-stream model are larger in some cases and smaller in others compared to those from the generalized model, depending on the cases involved. In the visible （VIS） band, canopy reflectance is smaller and canopy absorptance larger from the two-stream model compared to the generalized model when the Leaf Area Index （LAI） is low and soil reflectance is high. In cases of canopies with vertical leaf angles, the differences of reflectance and absorptance in the VIS and NIR bands between the two models are especially large. Two commonly occurring cases, with which the two-stream model cannot deal accurately, are also investigated. One is for a canopy with different adaxial and abaxial leaf optical properties; and the other is for incident sky diffuse radiation with a non-uniform distribution. Comparison of the generalized model within the same canopy for both uniform and non-uniform incident diffuse radiation inputs shows smaller differences in general. However, there is a measurable difference between these radiation inputs for a canopy with high leaf angle. This indicates that the application of the two-stream model to a canopy with different adaxial and abaxial leaf optical properties will introduce non-negligible errors.     

A Comparison of Two Canopy Radiative Models in Land Surface Processes

This paper compares the predictions by two radiative transfer models-the two-stream approximation model and the generalized layered model (developed by the authors) in land surface processes-for different canopies under direct or diffuse radiation conditions. The comparison indicates that there are significant differences between the two models, especially in the near infrared (NIR) band. Results of canopy reflectance from the two-stream model are larger than those from the generalized model. However, results of canopy absorptance from the two-stream model are larger in some cases and smaller in others compared to those from the generalized model, depending on the cases involved. In the visible (VIS) band, canopy reflectance is smaller and canopy absorptance larger from the two-stream model compared to the generalized model when the Leaf Area Index (LAI) is low and soil reflectance is high. In cases of canopies with vertical leaf angles, the differences of reflectance and absorptance in the VIS and NIR bands between the two models are especially large.Two commonly occurring cases, with which the two-stream model cannot deal accurately, are also investigated. One is for a canopy with different adaxial and abaxial leaf optical properties; and the other is for incident sky diffuse radiation with a non-uniform distribution. Comparison of the generalized model within the same canopy for both uniform and non-uniform incident diffuse radiation inputs shows smaller differences in general. However, there is a measurable difference between these radiation inputs for a canopy with high leaf angle. This indicates that the application of the two-stream model to a canopy with different adaxial and abaxial leaf optical properties will introduce non-negligible errors.     

A simplified scheme of the generalized layered radiative transfer model

In this paper, firstly, a simplified version （SGRTM） of the generalized layered radiative transfer model （GRTM） within the canopy, developed by us, is presented. It reduces the information requirement of inputted sky diffuse radiation, as well as of canopy morphology, and in turn saves computer resources. Results from the SGRTM agree perfectly with those of the GRTM. Secondly, by applying the linear superposition principle of the optics and by using the basic solutions of the GRTM for radiative transfer within the canopy under the condition of assumed zero soil reflectance, two sets of explicit analytical solutions of radiative transfer within the canopy with any soil reflectance magnitude are derived： one for incident diffuse, and the other for direct beam radiation. The explicit analytical solutions need two sets of basic solutions of canopy reflectance and transmittance under zero soil reflectance, run by the model for both diffuse and direct beam radiation. One set of basic solutions is the canopy reflectance αf （written as α1 for direct beam radiation） and transmittance βf （written as β1 for direction beam radiation） with zero soil reflectance for the downward radiation from above the canopy （i.e. sky）, and the other set is the canopy reflectance （αb） and transmittance βb for the upward radiation from below the canopy （i.e., ground）. Under the condition of the same plant architecture in the vertical layers, and the same leaf adaxial and abaxial optical properties in the canopies for the uniform diffuse radiation, the explicit solutions need only one set of basic solutions, because under this condition the two basic solutions are equal, i.e., αf = αb and βf = βb. Using the explicit analytical solutions, the fractions of any kind of incident solar radiation reflected from （defined as surface albedo, or canopy reflectance）, transmitted through （defined as canopy transmittance）, and absorbed by （defined as canopy absorptance） the canopy and other properties per     

Global radiation within corn

Global and net solar radiation profiles were measured by traversing sensors at four heights in a square-sown plot and a row plot of field corn on four relatively cloudless days in August 1972. The fluxes and their vertical distributions are discussed. A numerical model of the short-wave radiation fluxes in a canopy is presented. Using leaf area index and fixed leaf radiative properties, calculated values of radiation are within 10% of measurements in most instances. The performance under predominantly beam or diffuse radiation is similar and model values of crop albedo compare well with values calculated from radiation measurements above the crops.     

Optical properties of a vertically inhomogeneous mid-latitude mid-level mixed-phase altocumulus in the infrared region

A vertically inhomogeneous mid-latitude mixed-phase altocumulus cloud was observed around 17:26 UTC on Oct. 14, 2001 during the 9th Cloud Layer Experiment (CLEX9). In this study the microphysical and optical properties of this cloud are investigated on the basis of in-situ observed vertical profiles of particle size and habit distributions. Two cloud models, assuming that the cloud properties were vertically homogeneous and inhomogeneous, are adopted to derive the bulk optical and radiative properties of this cloud. The observed microphysical properties are combined with the theoretical solutions to the scattering and absorption properties of individual cloud particles to determine the bulk optical properties at various heights within the cloud layer. The single-scattering properties of spherical liquid water droplets and nonspherical ice crystals are obtained from the Lorenz–Mie theory and an existing database, respectively. The bulk microphysical and optical properties associated with the inhomogeneous model depend strongly on the height above the cloud-base whereas the dependence is smoothed out in the case of the homogeneous model. Furthermore, the transfer of infrared radiation is simulated in conjunction with the two cloud models. It is shown that the brightness temperatures at the top of the atmosphere in the case of the homogeneous model can be 1.5% (3.8 K) higher than their counterpart associated with the inhomogeneous cloud model. This result demonstrates that the effect of the vertical inhomogeneity of a mixed-phase cloud on its radiative properties is not negligible.     

陆面过程中冠层四流辐射传输模式的模拟性能检验

利用实际观测资料检验新发展的植被冠层四流辐射传输模式的模拟性能.将浙江大学田间观测数据用于模式试验,其结果表明:在可见光波段,因叶片散射率值较小,四流模式和二流模式模拟结果差别不大,模拟的冠层反照率都接近观测值.这表明在可见光波段,二流模式已经能够较好模拟芯层反照率,四流模式能够提高的精度范围有限;在近红外波段,因叶片散射率值较大,两个模式模拟结果差别较明显,四流模式模拟的冠层反照率相对二流模式的模拟结果更接近观测值.利用第2次雪模式比较计划SNOWMIP资料和改进的10层陆面模式BATS进行耦合试验,将四流模式以及二流模式均耦合到该陆面模式中,其目的是为考察四流模式对陆面模式模拟地表反射太刚短波辐射通量的影响.耦合后得到的结果与SNOWMIP中加拿大BERMS草地站和森林站的观测资料进行了对比.对比结果表明,采用网流模式、二流模式、BATS原辐射传输模式后,耦合四流模式的陆面模式模拟地表反射太阳短波辐射通量最接近观测值.说明采用四流模式能够改善陆面模式对地表反射太阳短波辐射通啦的模拟.     

阿尔卑斯山杉林冠层影响辐射传输的个例分析

利用瑞士Alptal观测站杉树林冠层上方、下方的辐射观测资料,分析了冠层对短波辐射的减弱及对长波辐射的增幅作用及其季节变化。结果表明,对比较密集的常绿针叶林,冠层对入射短波辐射的透过率随着太阳高度的降低而减小,春季以后趋于稳定;冠层对长波辐射的增幅作用随天气状况而变化,这种增幅作用在晴空条件下最显著,可达1.5倍。在冬季,因为太阳辐射较弱,冠层对长波辐射的增幅作用超过对短波辐射的减弱从而增加地面净辐射。在其它季节,太阳辐射比较强,冠层对短波辐射的减弱超过对长波辐射的增幅作用而减少地面净辐射。地面净辐射与冠层上方气温的变化趋势虽然在有些时段一致,但在伴随降雪过程的降温时段,地面净辐射与气温的变化趋势近乎反相,在积雪融化时段,地面净辐射的增加比气温升高更显著,尤其是在白天。     

Dense canopy albedo as a function of illumination direction: Dependence on structure and leaf transmittance

Summary The dependence of the albedo on illumination direction is analyzed by constructing a canopy model in which the individual leaves are planar, Lambertian reflectors. Leaf transmission is treated separately, and is assumed to be proportional to the cosine of the zenith angle of the leaf normal. Effects of shading and obscuration are formulated assuming random distribution of the leaves with respect to the viewing and illumination directions. Variants of the model with different azimuthal leaf distributions are created and discussed as explicit expressions of the viewing/illumination geometry and the canopy characteristics.The canopy spectral single-scattering albedo factors, which are the conditional probabilities of a photon escaping from the canopy after a first reflection or a first transmission, are obtained by numerical integration over a hemisphere of the bidirectional reflectance and transmittance factors. Our analysis identifies the ratio of the projection of leaf area on the vertical plane perpendicular to the principal plane to that on a horizontal plane as the parameter that controls the dependence of the albedo on the solar zenith angle. The albedo factor due to the leaf transmittance generally increases with the zenith angle of illuminating beam more sharply than that due to the leaf reflectance.Model variants with various azimuthal distributions are compared with measured albedo of soybeans. Second and higher order scatterings are accounted for in a simplified way. The degree to which a model variant fits the measured albedo and its change with the solar zenith angle depends both on the leaf inclination angleand the azimuthal distribution of the leaf area.With 6 Figures     

一个简单的陆面过程模式

本模式为针对大气环流模式所发展的一个简单的陆面过程模式，它包含:（1）地表温度计算，（2）冠层叶面贮水量和土壤湿度计算，（3）陆面与大气之间的水分和能量交换。对于表面温度和含水量的计算，采用的是联立求解计算方案，即耦合计算。植被冠层叶面的辐射特性和冠层形态对冠层中的辐射交换的影响得到有效和尽可能简单的模拟。另外，植被的气孔阻抗、表面与大气之间的水热交换通量和土壤中的水热输导作了较为细致的描写。利用此模式开展了对两个不同覆盖类型的陆面过程的模拟，模拟和观测的表面通量、温度和湿度较为相近。     

A Multi-layer Radiation Model for Urban Neighbourhoods with Trees

A neighbourhood-scale multi-layer urban canopy model of shortwave and longwave radiation exchange that explicitly includes the radiative effects of tall vegetation (trees) is presented. Tree foliage is permitted both between and above buildings, and mutual shading, emission and reflection between buildings and trees are included. The basic geometry is a two-dimensional canyon with leaf area density profiles and probabilistic variation of building height. Furthermore, the model accounts for three-dimensional path lengths through the foliage. Ray tracing determines the receipt of direct shortwave irradiance by building and foliage elements. View factors for longwave and shortwave diffuse radiation exchange are computed once at the start of the simulation using a Monte Carlo ray tracing approach; for subsequent model timesteps, matrix inversion rapidly solves infinite reflections and interception of emitted longwave between all elements. The model is designed to simulate any combination of shortwave and longwave radiation frequency bands, and to be portable to any neighbourhood-scale urban canopy geometry based on the urban canyon. Additionally, the model is sufficiently flexible to represent forest and forest-clearing scenarios. Model sensitivity tests demonstrate the model is robust and computationally feasible, and highlight the importance of vertical resolution to the performance of urban canopy radiation models. Full model evaluation is limited by the paucity of within-canyon radiation measurements in urban neighbourhoods with trees. Where appropriate model components are tested against analytic relations and results from an independent urban radiation transfer model. Furthermore, system response tests demonstrate the ability of the model to realistically distribute shortwave radiation among urban elements as a function of built form, solar angle and tree foliage height, density and clumping. Separate modelling of photosynthetically-active and near-infrared shortwave bands is shown to be important in some cases. Increased canyon height-to-width ratio and/or tree cover diminishes the net longwave radiation loss of individual canyon elements (e.g., floor, walls), but, notably, has little effect on the net longwave loss of the whole urban canopy. When combined with parametrizations for the impacts of trees on airflow and hydrological processes in the urban surface layer, the new radiation model extends the applicability of urban canopy models and permits more robust assessment of trees as tools to manage urban climate, air quality, human comfort and building energy loads.     

影响植被内部辐射状况的冠层结构特征研究

植被冠层结构特征是影响植被内部辐射状况的重要因子之一,对于植被群体的光合特性来说,冠层结构特性影响了光合有效辐射的截获、吸收和透射。本文着重分析了冠层结构的数量特征（如叶面积指数）和冠层几何结构特征,综合考虑叶倾角、叶方位角和Ｇ函数等群体几何特征,以便直接用于遥感光合植被生物量机理估测模型的建立。     

Penetration of sunlight into a canopy: Explicit models based on vertical and horizontal leaf projections

The projections of leaf areas onto a horizontal plane and onto a vertical plane are examined for their utility in characterizing canopies for sunlight penetration (direct beam only) models. These projections exactly specify the penetration if the projections on the principal plane of the normals to the top surfaces of the leaves are in the same quadrant as the Sun. Inferring the total leaf area from these projections (and therefore the penetration as a function of the total leaf area) is possible only with a large uncertainty (up to ±32%) because the projections are a specific measure of the total leaf area only if the leaf angle distribution is known. It is expected that this uncertainty could be reduced to more acceptable levels by making an approximate assessment of whether the zenith angle distribution is that of an extremophile canopy. An extremophile canopy would have the maximum leaf area possible for given set of projections. Simple leaf projection measurements would then become a practical substitute for detailed measurements of the leaf angle distribution. This is not true if for a fraction of the canopy, the leaf normal projections fall in the non-solar quadrant. In this case, accurate and detailed information about the leaf orientation is required for assessing the penetration; the horizontal and vertical projections are inadequate for this purpose.     

Modelling daily global radiation at the forest floor

A simple model is developed which estimates daily global radiation at the floor of a non-homogeneous Eucalyptus forest. Model input parameters are easily derived from field measurements and consist of individual tree location, tree height, maximum canopy width and its corresponding height, height of the lowest branch and trunk thickness. In addition, the model requires values for global and diffuse irradiance in the open. The tree canopy is represented as a series of spheres containing leaves which are homogeneously spaced but are oriented in the vertical plane. This configuration closely approaches that of actual eucalyptus trees. A Monte-Carlo approach is used to estimate the albedo of the unit sphere as a function of solar zenith angle. At a given combination of solar zenith and azimuth angle, the model estimates the solar irradiance at a specific forest floor location.The model, when tested against pyranometer measurements, predicted daily solar irradiance with a correlation of 0.98 and a standard error of 0.98 MJ m^-2 day^-1. This good performance is attributed to the spatial averaging of the radiation fluxes over the entire day, and the relatively low sensitivity of the calculated solar irradiance to sphere albedo.     

冰云短波辐射特性参数化

该文对卷云和高层云分别构造了15种冰晶尺度分布，3种云高和云厚，总共30种冰水含量、90种冰水程长模式。对尺度小于与大于30 μm的冰晶粒子的单次散射特性分别用表面积等效冰球Mie理论和射线光学理论进行计算。采用delta-Eddington法计算了在4种不同地面反射率下冰云的多重散射辐射传输特性，提出了冰云辐射特性的参数化公式，指出参数化公式可以用于气候模式或大气环流模式。     

THEORETICAL CALCULATION MODEL OF SHORTWAVE RADIATION FOR THE EARTH-ATMOSPHERE SYSTEM AND ITS TESTING RESULTS

On the basis of radiation transfer theory,adopting improved two-stream algorithm incorporated with addingalgorithm,we build up a theoretical calculation model of shortwave radiation for the earth-atmosphere system whichcan be applied with satellite data.The model can calculate direct solar radiation,scattering solar radiation,heating rateand other physical quantities of radiation field at every layer of the atmosphere and on the earth's surface,if the under-ground reflectance or the planetary albedo obtained from satellite can be known.The model can be used in clear orcloudy atmosphere,and its calculating speed is fairly fast.We think that the model can be incorporated into large-scaleand mesoscale climatic models for the consideration of radiation calculation,and also it is useful for the utilization of so-lar energy.     

Atmospheric radiative equilibria in a simple column model

 An analytic radiative-equilibrium model is formulated where both short- and longwave radiation are treated as two-stream (down- and upward) fluxes. An equilibrium state is defined in the model by the vertical temperature profile. The sensitivity of any such state to the model atmosphere’s optical properties is formulated analytically. As an example, this general formulation is applied to a single-column 11-layer model, and the model’s optical parameters are obtained from a detailed radiative parametrization of a general circulation model. The resulting simple column model is then used to study changes in the Earth-atmosphere system’s radiative equilibrium and, in particular, to infer the role of greenhouse trace gases, water vapor and aerosols in modifying the vertical temperature profile. Multiple equilibria appear when a positive surface-albedo feedback is introduced, and their stability is studied. The vertical structure of the radiative fluxes (both short- and longwave) is substantially modified as the temperature profile changes from one equilibrium to another. These equilibria and their stability are compared to those that appear in energy-balance models, which heretofore have ignored the details of the vertical temperature and radiation profiles. Received: 22 December 1995/Accepted: 17 January 1997     

盘锦湿地芦苇群落冠层内辐射分布特征

基于盘锦湿地芦苇群落冠层内不同层次的微气象要素与生物学特性观测,探讨芦苇群落冠层内总辐射分布与叶面积指数(LAI)的关系。结果表明:芦苇群落内太阳辐射的垂直分布在生长前期(5月)和后期(9月)呈S型曲线,生长盛期(6月和7月)呈指数曲线变化。冠层内不同层次的太阳辐射透射率随叶面积指数的增加而减少,消光系数(k)存在明显的季节变化,表现为从萌芽期到成熟期逐渐减小趋势(最小值k=0.436)。     

双叶模型在冬小麦田冠层CO2通量多层模拟中的应用

综合考虑农田生态系统中水、热、CO2输送所涉及的大气、水文、生物等生物物理过程,以Farquhar等提出的叶片尺度光合作用生物化学过程机理模型为理论基础,对其进行空间尺度扩展,并改进冠层分层方法,建立了均匀农田与大气之间物质输送和能量交换的多层模式,在模式中运用双叶模型,同时考虑叶片氮素水平垂直差异,对2008年4—5月华北平原冬小麦生长旺季农田生态系统中冠层CO2通量进行了模拟研究,并利用涡度相关观测的通量数据对模型的有效性加以验证,结果表明:在冠层多层空间,小麦拔节至孕穗期和开花至乳熟期叶片氮含量随冠层高度的衰减系数分别为0.793(R2=0.698)和1.374(R2=0.728),冠层内叶片氮含量的空间分布可以用以相对累积叶面积指数为自变量的函数来描述;模型分别计算各层阴、阳叶的光截取、气孔传导、光合作用等,最终计算冠层上方CO2通量,冬小麦农田净生态系统生产力模拟值与实测值相关显著(R2=0.78),模拟的CO2通量日变化特征晴天昼间比阴雨天和夜间的效果好;在考虑丛聚影响的叶片非随机分布的密集农田中,阴叶对总初始生产力的贡献率在35.7%左右,对生产力贡献很重要。分层统计显示,作物最终产量的形成主要...     

Measurements of chemically and biologically effective radiation reaching the ground

Solar ultraviolet radiation at the surface has been measured at Potsdam on cloudless days by spectrometer OL 752/10. The measurements are compared with broad-band filter measurements and with model calculations using a modified version of Green's model, which is independent of the measurements. Input data to the model such as atmospheric ozone and aerosol optical thicknesses were measured by a Dobson and Brewer spectrophotometer as well as a Linke Feussner pyrheliometer, respectively. Differences between the model and the measurements are discussed in terms of uncertainties in the calibration and errors of instruments as well as uncertainties in the model calculations including the errors of input data. It is demonstrated that different chemically and biologically effective radiances can be determined from only one set of measured spectral irradiance components, i.e. global radiation and diffuse downward and upward directed radiation. Examples of diurnal variations of the photochemical production of ozone and hydroxyl radicals as determined from spectral irradiance measurements and measured concentrations of relevant trace gases are given.It is shown from the measured irradiance that relations between different effects of radiation to the biosphere depend on solar zenith angle, and to a certain extent also on atmospheric ozone. This has to be taken into account when adverse effects of changing UV radiation are evaluated. Radiation Amplification Factors derived from measurements correspond to those determined from model calculations.     

A Vegetated Urban Canopy Model for Meteorological and Environmental Modelling

An urban canopy model is developed for use in mesoscale meteorological and environmental modelling. The urban geometry is composed of simple homogeneous buildings characterized by the canyon aspect ratio (h/w) as well as the canyon vegetation characterized by the leaf aspect ratio (σ _l ) and leaf area density profile. Five energy exchanging surfaces (roof, wall, road, leaf, soil) are considered in the model, and energy conservation relations are applied to each component. In addition, the temperature and specific humidity of canopy air are predicted without the assumption of thermal equilibrium. For radiative transfer within the canyon, multiple reflections for shortwave radiation and one reflection for longwave radiation are considered, while the shadowing and absorption of radiation due to the canyon vegetation are computed by using the transmissivity and the leaf area density profile function. The model is evaluated using field measurements in Vancouver, British Columbia and Marseille, France. Results show that the model quite well simulates the observations of surface temperatures, canopy air temperature and specific humidity, momentum flux, net radiation, and energy partitioning into turbulent fluxes and storage heat flux. Sensitivity tests show that the canyon vegetation has a large influence not only on surface temperatures but also on the partitioning of sensible and latent heat fluxes. In addition, the surface energy balance can be affected by soil moisture content and leaf area index as well as the fraction of vegetation. These results suggest that a proper parameterization of the canyon vegetation is prerequisite for urban modelling.