Fast and accurate radiance calculations using truncation approximation for anisotropic scattering phase functions

The phase function for solar light scattering by large particles such as cloud droplets is strongly anisotropic due to very strong peaking in the forward direction. This creates numerical difficulties when attempting to calculate accurate reflected and transmitted radiances, which are important for remote sensing of atmospheric and surface properties. A popular approach uses the delta function to approximate the forward-scattering peak in a fraction of energy and a limited number of polynomial terms or a geometrically truncated function for the remaining fraction (so-called truncation approximations). This article compares and discusses several methods for fast and accurate calculations using truncation approximations. When using a single truncation approximation for all scattering orders, large biases appear in directions near the solar and anti-solar points. As shown here, high accuracy can be obtained using different truncation approximations depending on the order of scattering. Of particular importance is the use of phase functions close to the exact phase functions for the first few orders of scattering. Applying the method in combination with the Monte Carlo (MC) method, in which the truncation fraction for a scattering order depends on the scattering angle at the previous scattering event, obtains accurate radiance calculations under almost all geometrical and optical conditions, including in directions near the solar point. Because the method also reduces computational noise due to the MC sampling of radiance, it is useful for fast and accurate radiance calculations for cloudy atmospheres.

[1]  Anthony B. Davis,et al.  3D Radiative Transfer in Cloudy Atmospheres , 2005 .

[2]  K. Stamnes,et al.  Radiative Transfer in the Atmosphere and Ocean , 1999 .

[3]  V. Budak,et al.  On the solution of a vectorial radiative transfer equation in an arbitrary three-dimensional turbid medium with anisotropic scattering , 2008 .

[4]  M. Modest.  Radiative heat transfer , 1993 .

[5]  Alexander A. Kokhanovsky,et al.  Light Scattering Media Optics: Problems and Solutions , 2010 .

[6]  Rhys Goldstein,et al.  Monte Carlo Simulation of Solar Reflectances for Cloudy Atmospheres , 2003 .

[7]  Estimates of radiation over clouds and dust aerosols: Optimized number of terms in phase function expansion , 2009 .

[8]  Hironobu Iwabuchi,et al.  Efficient Monte Carlo Methods for Radiative Transfer Modeling , 2006 .

[9]  Boris A. Kargin,et al.  The Monte Carlo Methods in Atmospheric Optics , 1980 .

[10]  Teruyuki Nakajima,et al.  Algorithms for radiative intensity calculations in moderately thick atmospheres using a truncation approximation , 1988 .

[11]  W. Wiscombe.  The Delta–M Method: Rapid Yet Accurate Radiative Flux Calculations for Strongly Asymmetric Phase Functions , 1977 .

[12]  T. Nakajima,et al.  Scaling algorithms for the calculation of solar radiative fluxes , 2007 .

[13]  Gary G. Gibson,et al.  δ-Fit: A fast and accurate treatment of particle scattering phase functions with weighted singular-value decomposition least-squares fitting , 2000 .

[14]  V. Antyufeev,et al.  Solution of the generalized transport equation with a peak-shaped indicatrix by the Monte Carlo method , 1996 .