Study on the Influence of Aerosol Radiation Balance in One-Dimensional Atmospheric Medium Using P n-Approximation Method

A numerical study, used for the influence of aerosol on the atmospheric radiation transfer, is conducted in this paper. Based on the established atmospheric radiation transfer model, we calculated the effect of species, concentration, and height distribution of aerosols on atmosphere radiation. The calculation results show that the aerosol particles affect the atmospheric radiation balance greatly and thus are the important component of the radiation balance of the earth-troposphere system.


Introduction
Atmospheric aerosols are suspensions of small solid or liquid particles with the diameters between 0.001 and 10 m.Aerosol affects the Earth's radiation budget by direct and indirect effects and plays an important role in global climate systems.The direct effect is that aerosol particles scatter and absorb the solar radiation and surface infrared radiation [1]; the indirect effect is that aerosol particles act as cloud condensation nuclei and ice nuclei [2,3].When the climate system is in balance, its absorbing solar radiation energy is accurately equal to the infrared radiation energy that the earth and atmosphere emit into outer space.Any factor that disturbs this balance and then changes the climate is the radiative forcing factor.Its forcing on the earth-troposphere system is the radiative forcing [4,5].Radiative forcing (in W ⋅ m −2 ) is defined as the change in the net radiative flux at the tropopause due to some factors change.Aerosol is a very important correction factor in the atmospheric radiation transfer.So it is significant to study the influence of aerosols particles on the atmospheric transfer.
Many numerical methods have been developed to solve the problem of radiative heat transfer, such as the discrete ordinates method (DOM) [6,7], finite volume method (FVM) [8], finite element method (FEM), Monte Carlo method (MCM) [9], and spherical harmonics method (SHM).The spherical harmonics method, also known as the   -approximation method, was originally proposed by Jeans [10] for interstellar radiative transfer and further developed by Davison [11] and Kourganoff [12] for neutron transport.Khouaja et al. [13] developed the   -approximation method for radiative heat transfer.Recently, Larsen et al. [14] used a double-  approximation for the transfer equation in neutron transport, but only to low order in the expansion.McClarren et al. [15] analyzed numerical oscillations of the   -approximation method and also suggested improvements based on wavelet analysis.
We established atmospheric radiative transfer calculations model based on the spherical harmonics method in onedimensional media and analyzed the influence of aerosols on the atmospheric radiative transfer.The calculations from the species, number, and distribution height of aerosols were achieved and the results were given in Figure 1.

Numerical Model
The radiative transfer equation for one-dimensional absorbing, emitting, and scattering gray medium can be expressed as where  is the optical thickness given by  = (  +   ),  is the scattering albedo with  =   /(  +   ),  is the medium thickness, and   and   are the absorbing and scattering coefficients, respectively.The radiative intensity (, ) and scattering phase function Φ(  , ) can be expanded in terms of associated Legendre polynomials   ().
Using recurrence relations between these polynomials, (1) becomes where   is an asymmetry parameter.Multiplying (2) by   () and integrating over  ∈ [−1, 1], we obtain using the orthogonality properties of the Legendre polynomials  + 1 differential equations for coefficients   () that comprise the  + 1   -approximation equations.Consider the following: For one-dimensional opaque diffuse emitting and reflecting interfaces, odd terms in the Legendre polynomials form an orthogonal set.The boundary conditions are given by where  is the emissivity.The boundary conditions imposed on the spherical harmonics functions are given by ( Mathematical Problems in Engineering 3 The boundary conditions [16] are applicable in the onedimensional case: where  0 (  ⇀ s ) is the circumstance radiation,   0 and   are reflectivities of the two sides of the boundary surfaces, respectively, when considering total reflection, and the values of these depend on the medium refractive index .
Separating the large refractive index from small, the diffuse reflectance at the boundary surfaces is given by Siegel and Spuckler [17]: In contrast, separating small refractive index from large, the diffuse reflectance at the boundary surfaces is given by where  = ( ℎ /  ),  ℎ ≥   .In treating a semitransparent interface, we can use opaque interface treatment as reference.
Considering one-dimensional unsteady coupled radiative and conductive heat transfer, the general form of the energy equation of the internal radiation and the transient coupled conductive heat transfer is defined as Considering unsteady state, no internal heat source and coupled radiative and conductive heat transfer, the energy equation ( 9) can be discretized as follows: In order to make the problem have universal significance, the dimensionless form of the energy equation ( 10) will be evaluated.Introducing the following definitions: reference temperature   , the dimensionless temperature T  =    /  , unit heat capacity  = , the dimensionless time Δ t = (4 2  3  /)Δ, t = (4 2  3  /), the dimensionless length Δ  = Δ  /, the radiative and conductive parameter  = /4 2  3  , the dimensionless convective heat transfer coefficient  = ℎ/4 2  3  , the dimensionless radiative heat flux Φ = Φ/(4 2  4  ), (10) becomes According to [18], in considering one-dimensional nonsteady coupled radiative and conductive heat transfers, the thermal radiation source term Φ  () is found to be Here, the incident radiative power () = ∫ Ω (, Ω)Ω can be transformed to () = 2 ∫ 1 −1 (, ).According to the properties of Legendre polynomials and  0 () = 1, the above equation is rearranged into the form: Meanwhile, considering the physical properties of the inhomogeneous medium, we next develop a  approximation method for a multilayer medium.We redetermine the boundary conditions of each layer.The modified boundary conditions are Here,  ,0 and  ,  are the emission radiation intensities of two adjacent layers, and the initial values are obtained from the temperature field.The total emission radiation intensity is determined by an iterative calculation of intensities from successive layers, the last providing the required total.
The convergence condition of the multilayer model calculation is decided by the radiative intensity   .After the intensity calculation has converged, the iterative calculation of the temperature field is then performed.

Results and Discussion
According to statistics, 98.8% of the total atmospheric mass distributes in the atmosphere within 30 km above the ground, and with the height increasing to 50 km the proportion reaches 99.92%.Most of the atmospheric radiation (involving solar radiation, absorption, and scattering of the ground infrared radiation and the emission of the atmosphere itself) occurs within certain atmospheric thickness.Compared with the earth radius (6378 km), the atmospheric thickness can be regarded as one-dimensional radiative media.Therefore, we treated the atmosphere within 50 km above the ground as a multilayer one-dimensional media distributing along the radial of the earth.The upper surface of the medium layer is near to space and is considered as a semitransparent medium.The lower interface is close to ground and is considered as opaque interface.The emissivity and radiative temperature related to the ground are given.For the inner of the medium layer we mainly consider the gas absorption and emission of the radiation energy and the aerosol particles scattering and absorption of the radiation energy and consider the atmospheric coupled radiation and conduction in the vertical direction.
The atmospheric convection is very important for the energy transfer of the atmosphere.But we mainly discuss the influence of aerosols on the atmospheric radiative energy and investigate the atmospheric radiative balance and the variation in radiation characteristic, which is not affected by the atmospheric convection that acts only as the background.Therefore, we ignore the influence of the atmospheric convection.
The atmospheric radiation physical model is shown in Figure 1 and the atmospheric radiation transfer model is shown in Figure 2. The atmospheric internal heat source is set as  = 0.
The standard atmospheric condition is given in Tables 1-4.
Table 1 gives the main input parameters for the calculation of atmospheric radiation transfer.
Table 2 gives the parameters associated with layers divided along the vertical height.The density of each layer is determined by the calculation of the pressure.We consider (/ 0 ) = (/ 0 ) (the standard atmosphere pressure  is  1.225 kg/m 2 ) with the temperature fluctuating in small range.Therefore, the equation is as follows: Table 3 gives the wavelength range of each spectral band.The atmosphere optical thicknesses of the layers at different spectral bands are obtained from the HITRAN database under the standard atmospheric condition of America in 1976.The aerosol optical thicknesses are decided based on the Mie scattering theory using the typical aerosol particles size distribution assumption and combining the complex refractive index of different aerosol species.Table 3 lists the spectral optical thicknesses of the dust aerosol and the scattering albedo parameters in the height of 0-1 km.The other parameters are not given in this paper as there are too much data.
Table 4 gives the complex refractive index of different aerosols.From the results in the figure we can see that as the different scattering effect of the aerosol particles, their influence on the temperature is greatly different.Compared with dust aerosol, the marine and water-soluble aerosol particles have stronger scattering effect.They will scatter upwards more radiation energy and decrease the atmospheric temperature.The soot particles have stronger absorption and their existence will cause the atmospheric temperature increasing.

The Influence of the Scattering Phase Function Variation.
The aerosol particles scattering characteristics are considered as forward and back scatter (the scattering phase functions are Φ 1 (Θ) and Φ 2 (Θ), resp.; see ( 16) and ( 17)).Other parameters are not changed.The results are shown in Figure 4: It can be seen from Figure 4(a) that, compared with the anisotropic scattering, the variation in the scattering phase functions causes the atmospheric temperature to change.But the effect is not obvious.Figure 4(b) shows the difference value between the atmospheric temperature fields and the isotropic scattering with different scattering phase functions.From the figure we can see that the forward and back scatters have opposite effect on the atmospheric temperature fields but the effect degree is similar.

The Influence of the Aerosol Particles Concentration.
The aerosol particles transport with the tropospheric movement will cause the aerosol particles concentrations to vary.In Figure 5 we consider the influence of the variation in total particles number concentrations, and the particles number concentrations are increased or decreased by 50%.So the radiation characteristics of the aerosol particles are changed, the atmospheric energy balance is affected, and then the atmospheric temperature fields are influenced.From the results it can be seen that the variation in aerosol particles concentrations can affect the atmospheric temperature but not obviously compared with other factors.

The Influence of the Height Distribution of Aerosol Particles.
We investigate the effect of the medium layer height on the aerosol particles size distribution; that is, the atmospheric radiation temperature fields change with the aerosol particles  being in different height medium layers.The calculation conditions are as follows: the aerosol particles distribute in the height of 0∼1 km, 1∼3 km, and 3∼9 km, respectively.The calculation results are in Figure 6.
From the calculation results in Figure 6 it can be seen that, with the aerosol particles transporting to the high altitude, the atmospheric temperature is affected greatly if the radiation characteristics of the aerosol particles have no changes.With the distribution height of the aerosol particles increasing, the medium temperature above the medium layer in which the aerosol particles exist increases greatly but the below temperature decreases obviously.This phenomenon shows the scattering effect of the aerosol particles on the solar radiation energy, that is, as the scattering effect of the atmospheric aerosols, most of the solar radiation energy to the ground is scattered back to the high altitude and cannot reach lower atmosphere, which causes the lower atmospheric temperature to decrease to about 200 K.

Conclusion
We analyzed the influence of the aerosol particles characteristics on the atmospheric radiation transfer model.The results showed the different effect of aerosol types on the atmospheric radiation.The dust, marine, and water-soluble  aerosols had negative effects on the atmospheric absorption, but the soot aerosol could increase the atmospheric absorption.The scattering phase functions affected the ground temperature.The variation of the aerosol particles distribution height could determine the atmospheric temperature; that is, the ground temperature decreased greatly with the aerosol particles existing in the high altitude stable.

Figure 1 :
Figure 1: Radiative heat transfer model diagram of one-dimensional atmospheric layer.

Figure 2 :Figure 3 :
Figure 2: Physical model of one plane-parallel slab divided into  sublayers.

Figure 4 :
Figure 4: Influence of scattering phase function change on atmospheric temperature.

Figure 5 :
Figure 5: Influence of the aerosol concentration on atmospheric temperature.

Figure 6 :
Figure 6: Influence of temperature fields on different aerosol distribution altitude.

Table 1 :
The main input parameters of the atmospheric radiation calculation.

Table 3 :
Spectral bands for the atmospheric radiative transfer model and radiative characteristic parameters under the standard dust aerosol in the height of 0-1 km.