New secondary-scattering correction in DISORT with increased efficiency for forward scattering

Abstract We present an alternative method to calculate the directional distribution after secondary scattering of light in an atmosphere, and apply it to the correction developed by Nakajima and Tanaka (1988) [1] as implemented in the DISORT radiative transfer solver. This method employs the scattering phase functions directly, instead of expanding over their Legendre moments as in the original formulation, and hence is not compromised in cases where a prohibitive number of moments is required to maintain accuracy. The new approach is designed to be particularly efficient in the strongly forward-scattering case, which arises for example in problems involving cloud-ice or dust particles. We have implemented this in a newly rewritten C-code version of DISORT that provides additional computational efficiencies via dynamic and cache-aware memory allocation. The new version uses less memory and runs considerably faster than the original, while producing results with equal or greater accuracy.

[1]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[2]  Claudia Emde,et al.  A polarized discrete ordinate scattering model for simulations of limb and nadir long‐wave measurements in 1‐D/3‐D spherical atmospheres , 2004 .

[3]  B. Mayer.  Radiative transfer in the cloudy atmosphere , 2009 .

[4]  Bernhard Mayer,et al.  Atmospheric Chemistry and Physics Technical Note: the Libradtran Software Package for Radiative Transfer Calculations – Description and Examples of Use , 2022 .

[5]  A. Kokhanovsky,et al.  SCIATRAN 2.0 – A new radiative transfer model for geophysical applications in the 175–2400 nm spectral region , 2004 .

[6]  K. Evans.  The Spherical Harmonics Discrete Ordinate Method for Three-Dimensional Atmospheric Radiative Transfer , 1998 .

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

[8]  B. Mayer,et al.  ALIS: An efficient method to compute high spectral resolution polarized solar radiances using the Monte Carlo approach , 2011, 1901.01842.

[9]  Knut Stamnes,et al.  On the computation of angular distributions of radiation in planetary atmospheres , 1982 .

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

[11]  Vladimir V. Rozanov,et al.  The solution of the vector radiative transfer equation using the discrete ordinates technique : Selected applications , 2006 .

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

[13]  P. Koepke,et al.  Optical Properties of Aerosols and Clouds: The Software Package OPAC , 1998 .

[14]  Tatsuya Yokota,et al.  Matrix formulations of radiative transfer including the polarization effect in a coupled atmosphere–ocean system , 2010 .

[15]  Bernhard Mayer,et al.  Solar radiation during a total eclipse: A challenge for radiative transfer , 2007 .

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

[17]  Bryan A. Baum,et al.  Bulk Scattering Properties for the Remote Sensing of Ice Clouds. Part I: Microphysical Data and Models. , 2005 .

[18]  J. Lenoble.  Radiative transfer in scattering and absorbing atmospheres: Standard computational procedures , 1985 .

[19]  Knut Stamnes,et al.  General Purpose Fortran Program for Discrete-Ordinate-Method Radiative Transfer in Scattering and Emitting Layered Media: An Update of DISORT , 2000 .

[20]  W. Wiscombe.  Improved Mie scattering algorithms. , 1980, Applied optics.

[21]  Bernhard Mayer,et al.  Efficient unbiased variance reduction techniques for Monte Carlo simulations of radiative transfer in cloudy atmospheres: The solution , 2011 .

[22]  C. Emde,et al.  Simulation of solar radiation during a total solar eclipse: a challenge for radiative transfer , 2007 .

[23]  G. Rybicki.  Radiative transfer , 2019, Climate Change and Terrestrial Ecosystem Modeling.