Radiative transfer, interception and scattering in coniferous forests : models and applications for production ecology and remote sensing

This work develops methods to account for shoot structure in models of coniferous canopy radiative transfer. Shoot structure, as it varies along the light gradient inside canopy, affects the efficiency of light interception per unit needle area, foliage biomass, or foliage nitrogen. The clumping of needles in the shoot volume also causes a notable amount of multiple scattering of light within coniferous shoots. The effect of shoot structure on light interception is treated in the context of canopy level photosynthesis and resource use models, and the phenomenon of within-shoot multiple scattering in the context of physical canopy reflectance models for remote sensing purposes. Light interception. A method for estimating the amount of PAR (Photosynthetically Active Radiation) intercepted by a conifer shoot is presented. The method combines modelling of the directional distribution of radiation above canopy, fish-eye photographs taken at shoot locations to measure canopy gap fraction, and geometrical measurements of shoot orientation and structure. Data on light availability, shoot and needle structure and nitrogen content has been collected from canopies of Pacific silver fir (Abies amabilis (Dougl.) Forbes) and Norway spruce (Picea abies (L.) Karst.). Shoot structure acclimated to light gradient inside canopy so that more shaded shoots have better light interception efficiency. Light interception efficiency of shoots varied about two-fold per needle area, about four-fold per needle dry mass, and about five-fold per nitrogen content. Comparison of fertilized and control stands of Norway spruce indicated that light interception efficiency is not greatly affected by fertilization. Light scattering. Structure of coniferous shoots gives rise to multiple scattering of light between the needles of the shoot. Using geometric models of shoots, multiple scattering was studied by photon tracing simulations. Based on simulation results, the dependence of the scattering coefficient of shoot from the scattering coefficient of needles is shown to follow a simple one-parameter model. The single parameter, termed the recollision probability, describes the level of clumping of the needles in the shoot, is wavelength independent, and can be connected to previously used clumping indices. By using the recollision probability to correct for the within-shoot multiple scattering, canopy radiative transfer models which have used leaves as basic elements can use shoots as basic elements, and thus be applied for coniferous forests. Preliminary testing of this approach seems to explain, at least partially, why coniferous forests appear darker than broadleaved forests in satellite data.

Comparative evaluation of inversion approaches of the radiative transfer model for estimation of crop biophysical parameters

The inversion of canopy reflectance models is widely used for the retrieval of vegetation properties from remote sensing. This study evaluates the retrieval of soybean biophysical variables of leaf area index, leaf chlorophyll content, canopy chlorophyll content, and equivalent leaf water thickness from proximal reflectance data integrated broadbands corresponding to moderate resolution imaging spectroradiometer, thematic mapper, and linear imaging self scanning sensors through inversion of the canopy radiative transfer model, PROSAIL. Three different inversion approaches namely the look-up table, genetic algorithm, and artificial neural network were used and performances were evaluated. Application of the genetic algorithm for crop parameter retrieval is a new attempt among the variety of optimization problems in remote sensing which have been successfully demonstrated in the present study. Its performance was as good as that of the look-up table approach and the artificial neural network was a poor performer. The general order of estimation accuracy for para-meters irrespective of inversion approaches was leaf area index > canopy chlorophyll content > leaf chlorophyll content > equivalent leaf water thickness. Performance of inversion was comparable for broadband reflectances of all three sensors in the optical region with insignificant differences in estimation accuracy among them.

Sun glint estimation in marine satellite images: a comparison of results from calculation and radiative transfer modeling

The intensity and location of Sun glint in two Medium Resolution Imaging Spectrometer (MERIS) images was modeled using a radiative transfer model that includes elevation features as well as the slope of the sea surface. The results are compared to estimates made using glint flagging and correction approaches used within standard atmospheric correction processing code. The model estimate gives a glint pattern with a similar width but lower peak level than any current method, or than that estimated by a radiative transfer model with surfaces that include slope but not height. The MERIS third reprocessing recently adopted a new slope statistics model for Sun glint correction; the results show that this model is an outlier with respect to both the elevation model and other slope statistics models and we recommend that its adoption should be reviewed.

Enhancement of optically thick to thin line intensities in solar and stellar coronal plasmas through radiative transfer effects: an angularly resolved study.

An analysis of radiative transfer effects present in the Fe XV ion stage of solar and stellar coronal plasmas provides a general explanation of line radiation intensity enhancement above the optically thin limit. Full linearization radiation transfer is compared with the escape factor method and found to be in good agreement at the lower column densities. An angular study of the enhancement shows that symmetry factors are of great importance. This gives a possible reason for the indeterminate status of opacity in relation to coronal lines of distant stellar sources, where only emission integrated across the whole surface is detected.

Multidimensional simulations of radiative transfer in Type Ia supernovae

A three-dimensional Monte Carlo code for modelling radiation transport in Type Ia supernovae is described. In addition to tracking Monte Carlo quanta to follow the emission, scattering and deposition of radiative energy, a scheme involving volume-based Monte Carlo estimators is used to allow properties of the emergent radiation field to be extracted for specific viewing angles in a multidimensional structure. This eliminates the need to compute spectra or light curves by angular binning of emergent quanta. The code is applied to two test problems to illustrate consequences of multidimensional structure on the modelling of light curves. First, elliptical models are used to quantify how large-scale asphericity can introduce angular dependence to light curves. Secondly, a model which incorporates complex structural inhomogeneity, as predicted by modern explosion models, is used to investigate how such structure may affect light-curve properties. © 2006 RAS.

Line-driven Disk Winds in Activ Galactic Nuclei: The Critical Importance of Ionization and Radiative Transfer

Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGNs) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga & Kallman. To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations.

A fast radiative transfer model for the assimilation of water vapor radiances from the Kalpana very high resolution radiometer

A fast radiative transfer model (RTM) to compute emitted infrared radiances for a very high resolution radiometer (VHRR), onboard the operational Indian geostationary satellite Kalpana has been developed and verified. This work is a step towards the assimilation of Kalpana water vapor (WV) radiances into numerical weather prediction models. The fast RTM uses a regression‐based approach to parameterize channel‐specific convolved level to space transmittances. A comparison between the fast RTM and the line‐by‐line RTM demonstrated that the fast RTM can simulate line‐by‐line radiances for the Kalpana WV channel to an accuracy better than the instrument noise, while offering more rapid radiance calculations. A comparison of clear sky radiances of the Kalpana WV channel with the ECMWF model first guess radiances is also presented, aiming to demonstrate the fast RTM performance with the real observations. In order to assimilate the radiances from Kalpana, a simple scheme for bias correction has been suggested.

The full-spectrum correlated-k method for longwave atmospheric radiative transfer using an effective Planck function

The correlated k-distribution (CKD) method is widely used in the radiative transfer schemes of atmospheric models and involves dividing the spectrum into a number of bands and then reordering the gaseous absorption coefficients within each one. The fluxes and heating rates for each band may then be computed by discretizing the reordered spectrum into of order 10 quadrature points per major gas and performing a monochromatic radiation calculation for each point. In this presentation it is shown that for clear-sky longwave calculations, sufficient accuracy for most applications can be achieved without the need for bands: reordering may be performed on the entire longwave spectrum. The resulting full-spectrum correlated k (FSCK) method requires significantly fewer monochromatic calculations than standard CKD to achieve a given accuracy. The concept is first demonstrated by comparing with line-by-line calculations for an atmosphere containing only water vapor, in which it is shown that the accuracy of heating-rate calculations improves approximately in proportion to the square of the number of quadrature points. For more than around 20 points, the root-mean-squared error flattens out at around 0.015 K/day due to the imperfect rank correlation of absorption spectra at different pressures in the profile. The spectral overlap of m different gases is treated by considering an m-dimensional hypercube where each axis corresponds to the reordered spectrum of one of the gases. This hypercube is then divided up into a number of volumes, each approximated by a single quadrature point, such that the total number of quadrature points is slightly fewer than the sum of the number that would be required to treat each of the gases separately. The gaseous absorptions for each quadrature point are optimized such that they minimize a cost function expressing the deviation of the heating rates and fluxes calculated by the FSCK method from line-by-line calculations for a number of training profiles. This approach is validated for atmospheres containing water vapor, carbon dioxide, and ozone, in which it is found that in the troposphere and most of the stratosphere, heating-rate errors of less than 0.2 K/day can be achieved using a total of 23 quadrature points, decreasing to less than 0.1 K/day for 32 quadrature points. It would be relatively straightforward to extend the method to include other gases.

The full-spectrum correlated k method for longwave atmospheric radiative transfer: treatment of gaseous overlap

The correlated k-distribution (CKD) method is widely used in the radiative transfer schemes of atmospheric models, and involves dividing the spectrum into a number of bands and then reordering the gaseous absorption coefficients within each one. The fluxes and heating rates for each band may then be computed by discretizing the reordered spectrum into of order 10 quadrature points per major gas, and performing a pseudo-monochromatic radiation calculation for each point. In this paper it is first argued that for clear-sky longwave calculations, sufficient accuracy for most applications can be achieved without the need for bands: reordering may be performed on the entire longwave spectrum. The resulting full-spectrum correlated k (FSCK) method requires significantly fewer pseudo-monochromatic calculations than standard CKD to achieve a given accuracy. The concept is first demonstrated by comparing with line-by-line calculations for an atmosphere containing only water vapor, in which it is shown that the accuracy of heating-rate calculations improves approximately in proportion to the square of the number of quadrature points. For more than around 20 points, the root-mean-squared error flattens out at around 0.015 K d−1 due to the imperfect rank correlation of absorption spectra at different pressures in the profile. The spectral overlap of m different gases is treated by considering an m-dimensional hypercube where each axis corresponds to the reordered spectrum of one of the gases. This hypercube is then divided up into a number of volumes, each approximated by a single quadrature point, such that the total number of quadrature points is slightly fewer than the sum of the number that would be required to treat each of the gases separately. The gaseous absorptions for each quadrature point are optimized such they minimize a cost function expressing the deviation of the heating rates and fluxes calculated by the FSCK method from line-by-line calculations for a number of training profiles. This approach is validated for atmospheres containing water vapor, carbon dioxide and ozone, in which it is found that in the troposphere and most of the stratosphere, heating-rate errors of less than 0.2 K d−1 can be achieved using a total of 23 quadrature points, decreasing to less than 0.1 K d−1 for 32 quadrature points. It would be relatively straightforward to extend the method to include other gases.