## Servicios Personalizados

## Revista

## Articulo

## Indicadores

- Citado por SciELO
- Accesos

## Links relacionados

- Citado por Google
- Similares en SciELO
- Similares en Google

## Compartir

## Gayana (Concepción)

##
*versión impresa* ISSN 0717-652X*versión On-line* ISSN 0717-6538

### Gayana (Concepc.) v.68 n.2 supl.TIIProc Concepción 2004

#### http://dx.doi.org/10.4067/S0717-65382004000300044

Gayana 68(2): 558-564, 2004
Keldysh Institute of Applied Mathematics, Russian Academy of Sciences, 4, Miusskaya Ploschadj, Moscow, 125047 Russia E-mail: tamaras@keldysh.ru
Boundary-value problem of radiation transfer in optically thick layers is considered while describing ocean, cloudiness, aerosol particle outbreaks, dust traces and other specific effects resulted from large fires (forest, peat-lands, in steppe regions, anthropogenic). Space and angular distributions of radiation inside the relevant layer of these media as well as reflected and passed through the layer radiation are formed as a result of multiple scattering and absorption. A new approach has been proposed to radiation transfer modeling in thick layers based on an influence function method. An internal border is inferred that separates the relevant layer on two regions. The influence function is calculated in the first region as a complete solution of the planar problem and as an asymptotic solution in the second region considering both cases of this last solution as that of a azimuth symmetric problem or as a complete solution of kinetic equation. The total radiation in the layer and outside it is found by a matrix functional known as the optical transfer operator the nuclei of which are given by the influence functions for each region. Forward and back hemispheres of the phase function serve as reflection and extinction coefficients of the internal border.
The different models and approximations of the optical transfer operator are employed in multidimentional problems of radiative correction for remote sensing of ocean, cloudiness, in theories of light and image propagation in turbin media, and in theoretical and computational fundamentals of optical-electronic remote sensing systems. Boundary-value problem of radiation transfer in optically thick finite planar layers is considered which are represented, in particular, by: · cloudiness; Spatial and angular distributions of radiation inside such layer as well as reflected and transmitted radiation are formed for these media as a result of multiple scattering of radiation and its absorption [1, 2]. One-dimensional models of solar radiation transfer in optically thick and not so thick layers of the two-media systems of the types like · "the atmosphere-ocean",· "the atmosphere-cloudiness", · "the atmosphere-hydrometeores", · "the atmosphere-gas/aerosol outbreaks cloud", · "the atmosphere-dust plume", i.e. in the atmosphere with admixtures resulted from vast fires influence (forest, peat, steppe, anthropogenic), for instance, which can result in an extended curtain, are elaborated by us using the following two approaches. The related calculations are conducted in the first model by the iterative characteristics method (ICM) [3] with procedures of speeding up the iteration convergence for the planar system with two layers on altitude: the lower layer cloudiness or the layer with admixtures, the upper layer free atmosphere. Radiation of the system in the second model is obtained by the functional that is called as the optical transfer operator (OTO) via the influence function (IF) of the atmosphere (IFA) and the layer with admixtures (IFLA) [3 5].
The first methodological investigations, algorithms and calculation programs were made by T.A. Sushkevich in 1965. Methods of speeding up the iteration convergence of successive approximations had been actively developed in the 1960-th and 1970-th (T.A. Germogenova and T.A. Sushkevich, V.Ya. Goldin Miln's problem and asymptotical methods had also been developed that time (M.V. Maslennikov and M.G. Kuzmina, T.A. Germogenova and N.V. Konovalov, G.A. Mikhailov A natural aspiration was to use the known approximate analytical & asymptotical solutions: - for modeling radiation transfer in two-media systems with the first layer being optically thin and the second optically thick; The asymptotical or diffuse radiation regime is known to be reached in any optically thick layer on a sufficiently large optical depth due to multiple scattering. Properties of the diffuse regime are maintained from physical considerations to be in the following: - the role of direct radiation (passed without scattering) is negligibly small comparing to the role of scattering radiation; Analytical solutions are obtained for any scattering layer with large optical depth which are expressed by asymptotic formulas of radiation transfer theory. Authors of such solutions often directly indicate, but sometimes keep silence that the errors of these analytical solutions are higher than 10-20% and even more in the layer nearest to the upper boundary and are increasing as the radiation absorption is enhancing. The errors due to not taking the azimuthal dependence into account have not yet been investigated. Theoretically and practically incorrect computation task is known to take the azimuthal dependence into account through short (3-5 terms) Fourier series. The authors of paper [2] have succeeded in the realization of the azimuthal harmonics calculation to determine the angular structure of reflected radiation from cloudiness while having in mind that this part of radiation is anisotropic and its description through albedo is rather inaccurate.
In fact, our approach (the iterative characteristics method/ICM with procedures of speeding up the iteration convergence and with taking large anisotropy of scattering into account) to the boundary problem solution of radiation transfer theory in optically thick homogeneous and inhomogeneous layers with three variables (altitude, zenith angle and azimuth) without any preliminary expansion in azimuthal harmonics series is the only method that enables to calculate on updated computers with large resources of operative and external memory and study angular and spatial distributions not only reflected and transmitted radiation, but that on any depth inside the layer. Speeding up procedures in the ICM: - relaxation speeding up (convergence speeding up not more than in two time); Radiation fields have been calculated by an iterative method of characteristics with speeding up procedures inside and outside of the layer. An external mono-directional flux is incident for the layer and the calculations have been conducted for a great number of optical models with different thickness of the layers, scattering phase functions, single scattering albedo all are corresponding to conservative scattering without absorption and weak, medium and heavy absorption as well. Additional calculations have been also conducted for the same models using asymptotic solutions. An adjustment of asymptotic regime through the layer thickness for outgoing and incident fluxes together with angular distributions and symmetry on azimuth has been investigated using updated graphical and visualization techniques. Azimuth harmonics have been calculated in accordance with the angular distributions. Applicability and accuracy problems for asymptotic solutions of the transfer equation have been studied. Numerical research has been conducted of how the asymptotical regime is established and what is the area of the asymptotical approach applicability. Radiation fields have been calculated by the iterative characteristics method inside and outside the layer for a large set of optical models (more than 200 calculation models) while an external mono-directional flux is incident with different: - directions of the solar flux (cosines are from 1.0 to 0.1), Calculations have been carried out for a number of the models based on the asymptotical solutions (using formulas of I.N. Melnikova, I.N. Minin, E.G. Yanovitsky, N.V. Konovalov).
An adjustment of the asymptotical regime along the layer thickness for outgoing and down-welling fluxes as well as angular distribution and azimuthal symmetry has been studied by updated graphical and visualization means. Azimuthal harmonics have been calculated using angular distributions inside the system and on its boundaries which are gained by the iterative method on the finite difference grid of large dimension (100 zenith angles, 91 azimuth, 50-100 layers). Applicability and accuracy problems of the asymptotical solutions of radiation transfer equation have been studied. Some illustrations of the conducted research are given at fig. (a h). Model of layer: optical thickness t Some important conclusions have been obtained based on the numerical studies: 1. about spatial and angular distributions of radiation intensity inside and on the boundaries of an optically thick layer;
1. About spatial and angular distributions of radiation intensity inside and on the boundaries of an optically thick layer: - dependence of reflected radiation from the thickness layer, phase function, degree of absorption, direction of solar flux; 2. About efficiency of the expansion procedure on azimuthal harmonics series: - practically all one-dimensional problems are solved by numerical and approximate analytical methods for single azimuthal harmonics with two variables (thickness and zenith angle); 3. About applicability and accuracy of asymptotic approximations for homogeneous and inhomogeneous optically thick layers; - all known approaches to solution of the problems for inhomogeneous layers utilize an adding procedure for the layers with reflection R and transmission T coefficients for separate layers; these coefficients are calculated by approximate asymptotical formulas without azimuthal dependence and the approaches have a low accuracy; 4. A new approach to modeling the radiation fields in optically thick layers based on the radiation transfer model in two media systems with optical transfer operator and with influence function method. This new approach is proposed to modeling the radiation transfer in optically thick finite and semi-infinite layers by the influence function method [3-5]. - The overall problem is in dividing a thick layer on sub-layers known methods can be used to and after that the problem is in obtaining the summation solution. This new original approach formulated within a universal linear-system approach is particular perspective for studying processes of solar radiation transformation in any two-media system, of short-wave radiation exchange between the atmosphere and any layer with admixtures as well as of light scattering regime in a layer with admixtures, in a free atmosphere, on the boundary "free atmosphere-layer with admixtures". This approach is particular effective if two media are such different that for modeling the radiation transfer in them modelers have to use various methods and different approximations to solve kinetic equation. The same approach is effective for modeling the angular and spatial radiation distributions in optically thick layers (clouds, dusts, etc.) which can be represented as two sub-layers with different radiation regimes.
About opportunities to use a new approach based on optical transfer operator and influence functions method: - Optically thick layer can be homogeneous or inhomogeneous on altitude; Our approach is completely different with the known adding layers method. In the adding layers method: - Only planar models; In IF and OTO method: - Research and analysis, but not only simple calculation of the multi-layer system; The influence functions method possesses an important priority: having available the IFA and IFLA set, it is possible: - To make different structures of radiation filed in any system; IFA and IFLA can be calculated besides that by various numerical and analytical methods with different degree of accuracy and approximation: - The influence function is calculated in the first area as a total solution of the planar problem, This approach enables to model radiation field in a wide range of variations of optical characteristics for the layers, to analyze mechanisms of radiation characteristics formation inside and outside the layers as well as to estimate the influence of each area. Reflected radiation and distribution of radiation characteristics inside any layer near to the boundary with an incident flux are calculated in the proposed approach with much higher accuracy. This statement is important to verify and improve remote sensing methods for clouds and outbreaks of different admixtures as well as methods of operative express-analysis of space data, in particular, in regions with natural environmental disasters including consequences of wide-scale fires.
This work has been supported by The Russian Foundation for Basic Research (projects 02-01-00135, 03-01-00132).
King M.D., Radke L.F., Hobbs P.V., 1990: Determination of the spectral absorption of solar radiation by marine stratocumulus clouds from airborne measurements within clouds // Journal of the Atmospheric Sciences. Vol. 47, 7, 894-907. [1] Melnikova I.N., Dlugach Z.M., Nakajima T., Kawamoto K., 2000: On the reflection function calculation in case of optically thick scattering layer, for Henyey-Greenstein phase function // Applied Optics. Vol. 39, 24, 4195-4204. [2] Sushkevich T.A., Strelkov S.A., Ioltukhovsky A.A., 1990: Characteristic Method to the Atmospheric Optics Problems. (in Russian) Moscow, Nauka. (Book review in Transport theory and statistical physics, 1993. Vol. 22, 4, 587-591.). [3] Sushkevich T.A., 1996: A solution to the boundary-value problem of the transport theory for a planar layer with a horizontally inhomogeneous boundary of an interface between two media // Doklady Mathematics. Vol. 54, 2, 775-779. Translated from Doklady Akademii Nauk. 1996. Vol. 350, 4, 460-464. [4] Sushkevich T.A., 2000: Linear-system approach and the theory of optical transfer operator // Atmospheric Oceanic Optics. V. 13, 8, 692-700. [5] |