Emergent radiation from the atmosphere bounded by the two half Lambert surfaces composed of different albedo,II

For the evaluation of the effect of the non-uniform surface albedo on the emergent radiation from the atmosphere, the emergent radiation from the atmosphere bounded by the two half Lambert surfaces composed of different albedo is computed. This paper is the improved version of the previous paper (Takashima and Masuda, 1991). The atmosphere is assumed to be homogeneous, which is composed of aerosol, molecules, and absorbent gases. Their optical thicknesses are (1) 0.25, 0.23, and 0.02, and (2) 0.75, 0.23, and 0.02, respectively. The model aerosol is of the oceanic and water soluble types.In the computational procedure, the emergent radiation is calculated approximately by the contributions due to the multiple scattering in the atmosphere, and due to the diffusely or directly transmitted radiation through the atmosphere which is reflected by the surfaces once (4 interactive radiative modes between atmosphere and surface). Furthermore, to perform the hemispherical integration processing the radiative interaction, the transmission function based on the single scattering in the atmosphere is introduced and then the transmission function is averaged over the hemisphere with weighting function. The numerical simulation exhibits the extraordinary effect near the two half surface boundary of different albedoes. The effect decreases exponentially with the distance from the boundary. The effect depends on the atmospheric aerosol type, optical thickness, and surface albedo. The present version enables us to quantitatively discuss the radiative transfer trend near the boundary of two half surfaces. The upward radiance would simply be evaluated using the present scattering approximation method if the surface albedo is less than 0.3. The present method is thought of as a first step extending the one-dimensional radiative transfer model to two-dimensional using the doubling-adding method.     

Atmospheric effect on the upwelling radiation at the top of the atmosphere over a stream

A new version is adopted for the evaluation of the upwelling radiation from atmosphere bounded by the surface, where the surface is composed of two half semi-infinite Lambert surfaces and a stream is inserted between them. The contrast of the stream is discussed with respect to the atmospheric effect. The width of the stream is considered to be 0.5, 1, and 3km; The solar and observational direction is located in the normal plane to the stream. The observational site is located at altitude 30km. The horizontal distance of observational site to the stream is fixed to 6.28 . The atmosphere is assumed to be homogeneous, which is composed of aerosol and molecules, where the model aerosol is of the oceanic type.In the computational procedure, a probability of radiation interacting with respective half surfaces and the stream are calculated based on the assumption of single scattering in the atmosphere, where isotropic scattering is undertaken. By use of this probability, the emergent radiation at the top of the atmosphere is calculated approximately by considering the radiative interactions between atmosphere and surfaces up to twice. The numerical simulation exhibits the extraordinary effect near the stream. The contrast of the stream depends upon the albedo of the surrounding surfaces. It increases with the increase of the stream width and decreases with the optical thickness.     

Simulation of atmospheric effect on the emergent radiation over a circular lake

The upwelling radiation at the top of the atmosphere is computed over a circular lake which is located in the uniform Lambert surface, using a modified version of the doubling-adding method. The radiance over the lake is discussed with respect to the atmospheric effect. The radius of the lake is assumed to be 0.5, 1, and 3 km. The observational site is located at altitude 30 km. The zenith of the observational site is located in the plane which is determined by the zenith of the center of the lake and incident solar direction. The zenith angle of the observational site to the center of the lake is fixed to 6.28°. The atmosphere is assumed to be homogeneous, which is composed of aerosol and molecule, where the model aerosol is of the oceanic or the water soluble types.Numerical simulation exhibits an extraordinary effect near the lake. The radiance of the lake against the surrounding depends upon the albedo of the surrounding surface. It increases with the increase of the size of the lake and decreases with the optical thickness. At large optical depth, the radiance depends upon the aerosol characteristics. It shows little dependence on the solar zenith angle if less than 60°.     

Simulation of atmospheric effects on the emergent radiation over a checkerboard type of terrain

To evalute the effect of the non-uniform surface on the radiation field, the upwelling radiation at the top of the atmosphere bounded by the checkerboard type of terrain is computed using the modified doubling method. The terrain is composed of the square Lambert surfaces with two different albedoes. The dimension of the each square is assumed to be 0.5–6 km. The radiance of the terrain is discussed with respect to the atmospheric effect. The observational site is located at altitude 30 km. The corresponding projected point on the ground is located at the center of a square. The solar and observational direction is located in the plane parallel to the checkerboard squares. The atmosphere is assumed to be homogeneous, which is composed of aerosol and molecules, where the model aerosol is of the oceanic or the water soluble types.Numerical simulation exhibits the extraordinary effect near the edge of each squares. The radiance of the terrain depends upon the difference of albedoes and size of squares. It increases with the increase of the dimension of the square. It decreases with the optical thickness. At large optical thickness, the variation of radiation with zenith direction depends upon the aerosol characteristics. It shows little dependence on the solar zenith angle if less than 20°.     

A new approach to evaluation of the effect of the two half-Lambert surfaces composed of different albedoes on the emergent radiation at the top of the atmosphere

A new way is adopted for the evaluation of the upwelling radiation from atmosphere bounded by two half-Lambert surfaces. The atmosphere is assumed to be homogeneous, and is composed of aerosol, molecules, and absorbent gases, where the model aerosol is of the oceanic and water soluble types.In the computational procedure, an iterative doubling-adding equation is expanded into a series of the radiative interaction modes between atmosphere and surface. Next, a probability of radiation interacting with respective half surfaces is calculated based on the assumption of single-scattering in the atmosphere. On the basis of this probability, the emergent radiation at the top of the atmosphere is approximately calculated by considering the radiative intractions to be twice as large. The effect of the multiple-scattering is fully taken into account. A numerical simulation exhibits the extraordinary effect near the two half-surface boundary of different albedoes. The effect of the other half-surface on the radiance decreases monotonically with the distance from the boundary. The present new version enable us to quantitatively discuss radiative transfer near the boundary of two half-surfaces even if the optical thickness is large and (or) surface albedo is great.     

A method of evaluation of the atmospheric effect on the emergent radiation over the infinite-cliff

To evaluate the effect of the cliff on the radiation field, the upwelling radiation at the top of the atmosphere is computed over the cliff using the reflection and transmission functions derived from the doubling-adding method. The model is defined by the plane-parallel homogeneous atmosphere, which is composed of aerosol and molecules, and is bounded by the top level surface, cliff and low level surface. These surfaces may be assumed to be the Lambertian.In the computational procedure, the equation for the emergent radiation is expanded into a series of radiative interaction modes among atmosphere, surfaces and the cliff. In respective modes, probabilities of respective interactions are firstly evaluated. With the aid of these probabilities, the emergent radiation is calculated using the doubling-adding method for the model atmosphere bounded by the surfaces and cliff, where the above radiative interactions are considered upto twice as large to obtain the enough accuracy of simulation. The multiple scattering is considered.     

A spherical model for computing polarized radiation in Titan's atmosphere

The Huygens descent through Titan's atmosphere in January 2005 will provide invaluable information about Titan's atmospheric composition and aerosol properties. The Descent Imager/Spectral Radiometer (DISR) will perform upward and downward looking radiation observations at various spectral ranges and spatial resolutions. To prepare the DISR data interpretation we have developed a new model for radiation transfer in Titan's atmosphere. The model solves for the full three-dimensional polarized radiation field in spherical geometry. However, the atmosphere itself is assumed to be spherically symmetric. The model is initialized with a fast-to-compute plane–parallel solution based on the doubling and adding algorithm that incorporates a spherical correction for the incoming direct solar beam. The full three-dimensional problem is then solved using the characteristics method combined with the Picard iterative approximation as described in Rozanov et al. (J. Quant. Spectrosc. Radiat. Transfer 69 (2001) 491). Aerosol scattering properties are calculated with a new microphysical model. In this formulation, aerosols are assumed to be fractal aggregates and include methane gas absorption embedded into the extinction coefficient. The resulting radiance of the model atmosphere's internal field is presented for two prescribed DISR wavelengths.     

Evaluation of the internal radiation field in a composite planetary atmosphere with the lambert surface calculated by matrix method

If the atmosphere is simulated by a number of homogeneous sublayers, it was shown (Takashima, 1973a) that the intensity and polarization parameters emerging from any boundary of internal sublayer's field can be determined, provided that the diffuse reflection and transmission matrices of radiation of sublayers are known. Furthermore, if the surface (ground) is assumed to reflect light in accordance with the Lambert law, it is shown in this paper that the internal radiation field at boundary of any sublayer can be also computed in terms of the diffuse radiation matrices of sublayers rather than in terms of that of the entire atmosphere (Sekera, unpublished). The effect of polarization is included.     

Reflection and transmission by a vertically inhomogeneous planetary atmosphere

By appealing to the reciprocity principle simple expressions are derived for the plane albedo and the transmissivity of a vertically inhomogeneous, plane parallel atmosphere. The plane albedo is shown to equal the angular distribution of the reflected intensity for isotropie Illumination of unit intensity incident at the top of the atmosphere, while the transmissivity equals the angular distribution of the transmitted intensity for isotropie illumination of unit Intensity incident at the bottom of the atmosphere. Chandrasekhar's solution of the planetary problem (including a Lambert reflecting lower boundary) in terms of the solution to the standard problem (no reflecting ground) is extended to apply to an inhomogeneous atmosphere resting on a surface that reflects radiation anisotropically but with no dependence on the direction of incidence (anisotropic Lambert reflector). The computational aspects are discussed and a procedure for computing the planetary albedo and transmissivity Is outlined for a vertically inhomogeneous, anisotropically scattering atmosphere overlying a partially reflecting surface. Numerical verification and illustration are also provided and it is shown that the assumed vertical variation of the single scattering albedo strongly affects the plane albedo but only weakly the transmissivity.     

Characteristics of the radiation emerging from the top of a Rayleigh atmosphere—II Total upward flux and albedo

The solution of the equation of radiative transfer in a medium exhibiting Rayleigh scattering, as developed by S. Chandrasekhar, has been used for an extensive series of computations(³) of the characteristics of the scattered and diffusely reflected radiation emerging from the top of an atmospheric model which corresponds in many respects to the sunlit portion of the earth's atmosphere. The first part of this two-part discussion dealt with the intensity, degree of polarization, plane of polarization and the neutral points of the emergent light as functions of sun elevation, direction in the downward hemisphere, optical thickness of the model atmosphere and reflectivity of the underlying surface. This second part is concerned with the upward flux obtained by an integration of the intensity over the entire hemisphere, for the incident radiation (a) being independent of wavelength or (b) having the spectral distribution of the extra-terrestrial solar radiation. Integration with respect to wavelength in the latter case, together with an approximation for the sphericity of the atmosphere, yields a value of 7.6 per cent for the earth's planetary albedo due to scattering by the clear atmosphere. An approximation for ozone absorption decreases the computed albedo to 6.9 per cent.     

Multiple scattering in the atmosphere with a rough surface

In the present paper, the intensity of radiation emergent from the atmosphere bounded by a rough surface is discussed with the aid of the superposition method derived by Mukai (1973). The merit of this method is to express the laws of diffuse reflection and transmission for the planetary problem with a rough surface in terms of a scattering and a transmission function for the standard problem.Here the bottom surface is assumed to reflect light in accordance with the slope distribution given by Cox and Munk (1954a, b). The results are discussed in terms of the optical properties and roughness of the bottom surface.     

Diffusion of radiation in planetary atmospheres

The problem of interaction of the solar radiation with the turbid Earth atmosphere, containing complicated polydispersive aerosol systems, is discussed in this paper. Equations for computing the angular functions ofn-th order scattering are derived. On the basis of these functions the spectral radiance, radiation flows and radiation balance of the atmosphere in the short-wave spectral range are calculated. The relations obtained can be used to calculate the complex index of refraction, distribution function and other characteristics of the submicron aerosol fraction, by solving the inverse problems.     

Inverse radiation modeling of Titan's atmosphere to assimilate solar aureole imager data of the Huygens probe

During the descent of the Huygens probe through Titan's atmosphere in January 2005, the Descent Imager/Spectral Radiometer (DISR) will perform upward and downward looking measurements at various spectral ranges and spatial resolutions. This internal radiation density could be estimated by radiative transfer calculations for Titan's atmosphere. However, to do this, the optical properties—i.e. volume extinction coefficient, single scattering albedo and scattering phase function—have to be prescribed at every altitude, and these are apriori not known. Herein, an inverse approach is investigated, which retrieves the single scattering albedo and the phase function of the aerosols from DISR observations. The method uses data from a DISR subinstrument, the Solar Aureole imager (SA), to estimate the optical properties of the atmospheric layer between two successive observation altitudes. A unique solution for one layer can in principle be calculated directly from a linear system of equations, but due to the sparseness of the data and the unavoidable noise in the measurements, the inverse problem is ill-posed. The problem is stabilized by the regularization method requiring smoothness of the resultant solution. A consistent set of solutions for all layers is obtained by iterating several times downward and upward through the layers. The method is tested in a simulated radiation density scenario for Titan, which is based on a microphysical aerosol model for the haze layer. Within this scenario, the expected coverage of SA data allows a reconstruction of the angular dependence of the scattering phase function with an explained variance better than 90%.     

Polarization of resonance lines in the case of a Voigt absorption profile

Multiple resonance scattering of radiation in a spectral line is considered in the case of a Voigt absorption profile. The scattering is assumed to take place in a nonmagnetic semi-infinite atmosphere with uniformly distributed sources of unpolarized radiation. Polarization characteristics have been obtained for the emergent radiation by numerically solving the Ambartsumian-Chandrasekhar matrix integral equation.     

Rayleigh-Cabannes scattering in planetary atmospheres II. Constant internal sources in nonconservative atmospheres

The vector equation of radiative transfer is solved for non-conservative homogeneous plane-parallel atmosphere using the method of discrete ordinates. The scattering processes in the atmosphere bounded by a Lambert bottom are described by the Rayleigh-Cabannes phase matrix. The primary radiation field is generated by constant internal sources. A package of FORTRAN subroutines is compiled to find the axial radiation field for such an atmosphere at arbitrary optical depth.     

A new method for exact solution of transfer equations in finite media

In this paper we develop a new exact method combined with finite Laplace transform and theory of linear singular operators to obtain a solution of transport equation in finite plane-parallel steady-state scattering atmosphere both for angular distribution of radiation from the bounding faces of the atmosphere and for intensity of radiation at any depth of the atmosphere. The emergent intensity of radiation from the bounding faces are determined from simultaneous linear integral equations of the emergent intensity of radiation in terms ofX andY equations of Chandrasekhar. The intensity of radiation at any optical depth for a positive and negative direction parameter is derived by inversion of the Laplace transform in terms of intergrals of the emergent intensity of radiation. A new expression of theX andY equation is also derived for easy numerical computation. This is a new and exact method applicable to all problems in finite plane parallel steady scattering atmosphere.     

A discrete spherical harmonics method for radiative transfer analysis in inhomogeneous polarized planar atmosphere

A discrete spherical harmonics method is developed for the radiative transfer problem in inhomogeneous polarized planar atmosphere illuminated at the top by a collimated sunlight while the bottom reflects the radiation. The method expands both the Stokes vector and the phase matrix in a finite series of generalized spherical functions and the resulting vector radiative transfer equation is expressed in a set of polar directions. Hence, the polarized characteristics of the radiance within the atmosphere at any polar direction and azimuthal angle can be determined without linearization and/or interpolations. The spatial dependent of the problem is solved using the spectral Chebyshev method. The emergent and transmitted radiative intensity and the degree of polarization are predicted for both Rayleigh and Mie scattering. The discrete spherical harmonics method predictions for optical thin atmosphere using 36 streams are found in good agreement with benchmark literature results. The maximum deviation between the proposed method and literature results and for polar directions \(\vert \mu \vert \geq0.1 \) is less than 0.5% and 0.9% for the Rayleigh and Mie scattering, respectively. These deviations for directions close to zero are about 3% and 10% for Rayleigh and Mie scattering, respectively.     

Polarization effect on radiative transfer in planetary composite atmospheres with interacting interface

A procedure of computing the radiance and the polarization parameters of radiation diffusely reflected and transmitted by an inhomogeneous, plane-parallel terrestrial atmosphere bounded by a ruffled ocean surface is discussed with the aid of the adding method. If the atmosphere and the ocean are simulated by a number of homogeneous sublayers, the matrices of radiation reflected and transmitted diffusely by the atmosphere-ocean system can be expressed in terms of these matrices of sublayers by using only a couple of iterative equations in which the polarity effect of radiation is included. Furthermore, the upwelling radiance and the polarization degree of radiation at the top of the atmosphere can be calculated by using a single iterative equation without requiring the equation for the diffuse transmission matrix of radiation. The ruffled ocean surface can be treated as an interacting interface, where the transmitted radiation from beneath the ocean surface into the atmosphere is also taken into account into the derivation of equations. Finally, sample computations of the upwelling radiance and the polarization degree of radiation from the top of the atmosphere are carried out at the wavelength of 0.60 micron.     

Radiation transfer through atmospheric aerosol media with anisotropic scattering

Radiation transfer in atmospheric aerosol media with general boundary conditions has been studied for anisotropic scattering. The considered aerosol medium assumed to have specular and diffused reflecting boundary surfaces and in the presence of internal source. The radiation transfer scattering parameters as single scattering albedo, asymmetry factor, scattering, absorption, extinction efficiencies and anisotropic scattering coefficient have been calculated using the Mie theory. The problem with general boundary conditions is solved in terms of the solution of source-free problem with simply boundary conditions. Pomraning-Eddington approximation is used to solve the source-free problem. For the sake of comparison, a weight function is introduced and used in two special forms. The calculated partial heat fluxes with the two methods are compared and showed good agreement. Some of our results are found in a good agreement with published data.     

Average photon path-length of radiation emerging from finite atmospheres

The determination of the average path-length of photons emerging from a finite planeparallel atmosphere with molecular scattering is discussed. We examine the effects of polarisation on the average path-length of the emergent radiation by comparing the results with those obtained for the atmosphere where the scattering obeys the scalar Rayleigh function. Only the axial radiation field is considered for both cases.To solve this problem we have used the integro-differential equations of Chandrasekhar for the diffuse scattering and transmission functions (or matrices). By differentiation of these equations with respect to the albedo of single scattering we obtain new equations the solution of which gives us the derivatives of the intensities of the emergent radiation at the boundaries.As in the case of scalar transfer the principles of invariance by Chandrasekhar may be used to find an adding scheme to obtain both the scattering and transmission matrices and their derivatives with respect to the albedo of single scattering. These derivatives are crucial in determining the average path length.The numerical experiments have shown that the impact of the polarisation on the average pathlength of the emergent radiation is the largest in the atmospheres with optical thickness less than, or equal to, three, reaching 6.9% in the reflected radiation.