A Path‐Tracing Monte Carlo Library for 3‐D Radiative Transfer in Highly Resolved Cloudy Atmospheres

Interactions between clouds and radiation are at the root of many difficulties in numerically predicting future weather and climate and in retrieving the state of the atmosphere from remote sensing observations. The large range of issues related to these interactions, and in particular to three-dimensional interactions, motivated the development of accurate radiative tools able to compute all types of radiative metrics, from monochromatic, local and directional observables, to integrated energetic quantities. In the continuity of this community effort, we propose here an open-source library for general use in Monte Carlo algorithms. This library is devoted to the acceleration of path-tracing in complex data, typically high-resolution large-domain grounds and clouds. The main algorithmic advances embedded in the library are those related to the construction and traversal of hierarchical grids accelerating the tracing of paths through heterogeneous fields in null-collision (maximum cross-section) algorithms. We show that with these hierarchical grids, the computing time is only weakly sensitivive to the refinement of the volumetric data. The library is tested with a rendering algorithm that produces synthetic images of cloud radiances. Two other examples are given as illustrations, that are respectively used to analyse the transmission of solar radiation under a cloud together with its sensitivity to an optical parameter, and to assess a parametrization of 3D radiative effects of clouds.

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

[2]  Roger Davies,et al.  Transport of infrared radiation in cuboidal clouds , 1981 .

[3]  J. Gautrais,et al.  Integral formulation of null-collision Monte Carlo algorithms , 2013 .

[4]  Donald E. Knuth,et al.  Literate Programming , 1984, Comput. J..

[5]  Christophe Spiesser,et al.  Radiative, conductive and convective heat-transfers in a single Monte Carlo algorithm , 2016 .

[6]  Jacques Gautrais,et al.  Addressing nonlinearities in Monte Carlo , 2016, Scientific Reports.

[7]  E. Mlawer,et al.  Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave , 1997 .

[8]  B. Barkstrom,et al.  Cloud-Radiative Forcing and Climate: Results from the Earth Radiation Budget Experiment , 1989, Science.

[9]  Anthony B. Davis,et al.  Nonlocal independent pixel approximation: direct and inverse problems , 1998, IEEE Trans. Geosci. Remote. Sens..

[10]  B. Mayer,et al.  Effects of 3-D thermal radiation on the development of a shallow cumulus cloud field , 2017 .

[11]  T. Mckee,et al.  Effects of Cloud Shape on Scattered Solar Radiation , 1978 .

[12]  B. Mayer,et al.  Representing 3‐D cloud radiation effects in two‐stream schemes: 1. Longwave considerations and effective cloud edge length , 2016 .

[13]  R. Hogan,et al.  Entrapment: An Important Mechanism to Explain the Shortwave 3D Radiative Effect of Clouds , 2019, Journal of the Atmospheric Sciences.

[14]  E. Gelbard,et al.  Monte Carlo Principles and Neutron Transport Problems , 2008 .

[15]  Hideki Kobayashi,et al.  A coupled 1-D atmosphere and 3-D canopy radiative transfer model for canopy reflectance, light environment, and photosynthesis simulation in a heterogeneous landscape , 2008 .

[16]  Ares Lagae,et al.  Out‐of‐Core Construction of Sparse Voxel Octrees , 2013, HPG '13.

[17]  Anthony B. Davis,et al.  Radiative smoothing in fractal clouds , 1995 .

[18]  A. de Latailladea,et al.  Monte Carlo method and sensitivity estimations , 2002 .

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

[20]  Kei Iwasaki,et al.  Toward Optimal Space Partitioning for Unbiased, Adaptive Free Path Sampling of Inhomogeneous Participating Media , 2011, Comput. Graph. Forum.

[21]  Rintaro Okamura,et al.  Multispectral Monte Carlo radiative transfer simulation by the maximum cross-section method , 2017 .

[22]  H. Barker,et al.  Estimation of Errors in Two-Stream Approximations of the Solar Radiative Transfer Equation for Cloudy-Sky Conditions , 2015 .

[23]  B. Mayer,et al.  A three-dimensional parallel radiative transfer model for atmospheric heating rates for use in cloud resolving models—The TenStream solver , 2015 .

[24]  Leonidas J. Guibas,et al.  Bidirectional Estimators for Light Transport , 1995 .

[25]  P. Mascart,et al.  Giga-LES of Hector the Convector and Its Two Tallest Updrafts up to the Stratosphere , 2016 .

[26]  G. A. Mikhaĭlov,et al.  Optimization of Weighted Monte Carlo Methods , 1992 .

[27]  S. Blanco,et al.  The practice of recent radiative transfer Monte Carlo advances and its contribution to the field of microorganisms cultivation in photobioreactors , 2013 .

[28]  Frédéric Parol,et al.  Cloud heterogeneity on cloud and aerosol above cloud properties retrieved from simulated total and polarized reflectances , 2017, Atmospheric Measurement Techniques.

[29]  Leonidas J. Guibas,et al.  Robust Monte Carlo methods for light transport simulation , 1997 .

[30]  B. Mayer,et al.  Representing 3‐D cloud radiation effects in two‐stream schemes: 2. Matrix formulation and broadband evaluation , 2016 .

[31]  K. Evans,et al.  Three-dimensional solar radiative transfer in small tropical cumulus fields derived from high-resolution imagery , 2001 .

[32]  W. Collins,et al.  Radiative forcing by long‐lived greenhouse gases: Calculations with the AER radiative transfer models , 2008 .

[33]  D. Glouchkov,et al.  Monte Carlo simulation of photon transport for optically thick, differentially moving plasmas: II. Escape factors for differentially moving spheres , 2003 .

[34]  Brent Burley,et al.  The Design and Evolution of Disney’s Hyperion Renderer , 2018, ACM Trans. Graph..

[35]  A. Marshak,et al.  A method for analyzing how various parts of clouds influence each other's brightness , 2003 .

[36]  A. Marshak,et al.  Solar radiation transport in the cloudy atmosphere: a 3D perspective on observations and climate impacts , 2010 .

[37]  Sandrine Bony,et al.  An Assessment of the Primary Sources of Spread of Global Warming Estimates from Coupled Atmosphere–Ocean Models , 2008 .

[38]  Ken Museth,et al.  VDB: High-resolution sparse volumes with dynamic topology , 2013, TOGS.

[39]  Robert F. Cahalan,et al.  The I3RC - Bringing Together the Most Advanced Radiative Transfer Tools for Cloudy Atmospheres , 2005 .

[40]  B. Mayer,et al.  The role of 1-D and 3-D radiative heating in the organization of shallow cumulus convection and the formation of cloud streets , 2017 .

[41]  R. Davies,et al.  The Effect of Finite Geometry on the Three-Dimensional Transfer of Solar Irradiance in Clouds , 1978 .

[42]  Andrew A. Lacis,et al.  Scattering, Absorption, and Emission of Light by Small Particles , 2002 .

[43]  Klaus Pfeilsticker,et al.  The Monte Carlo atmospheric radiative transfer model McArtim: Introduction and validation of Jacobians and 3D features , 2011 .

[44]  J. Potter The Delta Function Approximation in Radiative Transfer Theory , 1970 .

[45]  V. Eymet,et al.  Null-collision meshless Monte-Carlo-Application to the validation of fast radiative transfer solvers embedded in combustion simulators , 2013 .

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

[47]  B. Mayer,et al.  The Neighboring Column Approximation (NCA) - A fast approach for the calculation of 3D thermal heating rates in cloud resolving models , 2016 .

[48]  Robin J. Hogan,et al.  A Flexible and Efficient Radiation Scheme for the ECMWF Model , 2018, Journal of Advances in Modeling Earth Systems.

[49]  S. Bony,et al.  Marine boundary layer clouds at the heart of tropical cloud feedback uncertainties in climate models , 2005 .

[50]  Ingo Wald,et al.  Embree: a kernel framework for efficient CPU ray tracing , 2014, ACM Trans. Graph..

[51]  B. A. Kargin,et al.  Elements of Radiative-Transfer Theory Used in the Monte Carlo Methods , 1980 .

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

[53]  J. Redelsperger,et al.  A turbulence scheme allowing for mesoscale and large‐eddy simulations , 2000 .

[54]  Evgueni I. Kassianov,et al.  Spectral dependence of radiative horizontal transport in stratocumulus clouds and its effect on near-IR absorption , 2002 .

[55]  M. Wendisch,et al.  IPRT polarized radiative transfer model intercomparison project – Phase A , 2015, 1901.01813.

[56]  J. Deardorff Stratocumulus-capped mixed layers derived from a three-dimensional model , 1980 .

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

[58]  Roger Davies,et al.  Effects of Cloud Heterogeneities on Shortwave Radiation: Comparison of Cloud-Top Variability and Internal Heterogeneity , 1999 .

[59]  Jan Novák,et al.  Residual ratio tracking for estimating attenuation in participating media , 2014, ACM Trans. Graph..

[60]  Véronique Ducrocq,et al.  The Meso-NH Atmospheric Simulation System. Part I: adiabatic formulation and control simulations , 1997 .

[61]  Frédéric André,et al.  A symbolic approach for the identification of radiative properties , 2017 .

[62]  B. Mayer,et al.  A Fast Three-Dimensional Approximation for the Calculation of Surface Irradiance in Large-Eddy Simulation Models , 2008 .

[63]  Michael I. Mishchenko,et al.  Bidirectional Reflectance of Flat, Optically Thick Particulate Layers: An Efficient Radiative Transfer Solution and Applications to Snow and Soil Surfaces , 1999 .

[64]  James T. Kajiya,et al.  The rendering equation , 1986, SIGGRAPH.

[65]  Meng Gao,et al.  The impacts of 3D radiative transfer effects on cloud radiative property simulations and retrievals , 2013 .

[66]  Paul W. Stackhouse,et al.  The Effect of Cumulus Cloud Field Anisotropy on Domain-Averaged Solar Fluxes and Atmospheric Heating Rates , 2007 .

[67]  Robin J. Hogan,et al.  Tripleclouds: An Efficient Method for Representing Horizontal Cloud Inhomogeneity in 1D Radiation Schemes by Using Three Regions at Each Height , 2008 .

[68]  K Weise,et al.  Uncertainty treatment in Monte Carlo simulation , 1997 .

[69]  Leonidas J. Guibas,et al.  Metropolis light transport , 1997, SIGGRAPH.

[70]  Turner Whitted,et al.  An improved illumination model for shaded display , 1979, CACM.

[71]  W. Rossow,et al.  The International Satellite Cloud Climatology Project (ISCCP) Web Site An Online Resource for Research , 2004 .

[72]  G. Caniaux,et al.  A Numerical Study of the Stratiform Region of a Fast-Moving Squall Line. Part II: Relationship between Mass, Pressure, and Momentum Fields , 1995 .

[73]  W. A. Coleman Mathematical Verification of a Certain Monte Carlo Sampling Technique and Applications of the Technique to Radiation Transport Problems , 1968 .

[74]  Markus Wagner,et al.  Interactive Rendering with Coherent Ray Tracing , 2001, Comput. Graph. Forum.

[75]  Hironobu Iwabuchi,et al.  Fast and accurate radiance calculations using truncation approximation for anisotropic scattering phase functions , 2009 .

[76]  S. Blanco,et al.  Monte carlo estimates of domain-deformation sensitivities. , 2005, Physical review letters.

[77]  A. Macke,et al.  Monte Carlo radiative transfer calculations for inhomogeneous mixed phase clouds , 1999 .

[78]  C. Lac,et al.  Evaluation of turbulence parametrizations in convective clouds and their environment based on a large‐eddy simulation , 2019, Quarterly Journal of the Royal Meteorological Society.

[79]  Robert Pincus,et al.  The Accuracy of Determining Three-Dimensional Radiative Transfer Effects in Cumulus Clouds Using Ground-Based Profiling Instruments , 2005 .

[80]  R. Ellingson,et al.  Scattering effects on longwave fluxes in broken cloud fields , 1996 .

[81]  J. Joseph,et al.  The delta-Eddington approximation for radiative flux transfer , 1976 .

[82]  C. Gautier,et al.  A Three-Dimensional Radiative Transfer Model to Investigate the Solar Radiation within a Cloudy Atmosphere. Part II: Spectral Effects , 1998 .

[83]  Inna,et al.  Radiation in numerical weather prediction , 2017 .

[84]  G. Caniaux,et al.  A Numerical Study of the Stratiform Region of a Fast-Moving Squall Line. Part I: General Description and Water and Heat Budgets , 1994 .

[85]  Q. Fu,et al.  On the correlated k-distribution method for radiative transfer in nonhomogeneous atmospheres , 1992 .

[86]  Katsuhisa Koura,et al.  Null‐collision technique in the direct‐simulation Monte Carlo method , 1986 .

[87]  Robert L. Cook,et al.  Distributed ray tracing , 1984, SIGGRAPH.

[88]  L. Girolamo,et al.  Design and Verification of a New Monochromatic Thermal Emission Component for the I3RC Community Monte Carlo Model , 2018 .

[89]  M. Wendisch,et al.  IPRT polarized radiative transfer model intercomparison project - Three-dimensional test cases (phase B) , 2019, 1901.01828.

[90]  Céline Cornet,et al.  Effects of cirrus heterogeneity on lidar CALIOP/CALIPSO data , 2017 .

[91]  Martin Köhler,et al.  Observations and numerical simulations of the diurnal cycle of the EUROCS stratocumulus case , 2004 .

[92]  E. Clothiaux,et al.  Assessing 1D atmospheric solar radiative transfer models: Interpretation and handling of unresolved clouds , 2003 .

[93]  Alexander Marshak,et al.  Solar zenith and viewing geometry‐dependent errors in satellite retrieved cloud optical thickness: Marine stratocumulus case , 2009 .

[94]  H R Skullerud,et al.  The stochastic computer simulation of ion motion in a gas subjected to a constant electric field , 1968 .

[95]  Jean-François Cornet,et al.  Monte Carlo advances and concentrated solar applications , 2014 .

[96]  S. Blanco,et al.  Radiative transfer and spectroscopic databases: A line-sampling Monte Carlo approach , 2016 .

[97]  Johannes Hanika,et al.  Monte Carlo Methods for Volumetric Light Transport Simulation , 2018, Comput. Graph. Forum.

[98]  A. P. Siebesma,et al.  A Large Eddy Simulation Intercomparison Study of Shallow Cumulus Convection , 2003 .

[99]  Jean-Christophe Golaz,et al.  Large‐eddy simulation of the diurnal cycle of shallow cumulus convection over land , 2002 .

[100]  Pierre Aumond,et al.  Overview of the Meso-NH model version 5.4 and its applications , 2018, Geoscientific Model Development.

[101]  R. Hogan,et al.  Incorporating the Effects of 3D Radiative Transfer in the Presence of Clouds into Two-Stream Multilayer Radiation Schemes , 2013 .

[102]  Hartwig Deneke,et al.  Large‐eddy simulations over Germany using ICON: a comprehensive evaluation , 2017 .

[103]  Anthony B. Davis,et al.  Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated , 1994 .

[104]  T. Smith,et al.  The C.I.E. colorimetric standards and their use , 1931 .

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

[106]  Alexander Keller,et al.  Unbiased Global Illumination with Participating Media , 2008 .

[107]  Z. X. Li,et al.  Interpretation of Cloud-Climate Feedback as Produced by 14 Atmospheric General Circulation Models , 1989, Science.

[108]  G. Feingold,et al.  Cloud droplet growth in shallow cumulus clouds considering 1-D and 3-D thermal radiative effects , 2019, Atmospheric Chemistry and Physics.

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

[110]  B. Mayer,et al.  Validating the MYSTIC three-dimensional radiative transfer model with observations from the complex topography of Arizona's Meteor Crater , 2010 .

[111]  J. Dufresne,et al.  A boundary-based net-exchange Monte Carlo method for absorbing and scattering thick media , 2005, 1902.09966.

[112]  S. Schafer What is the global impact of 3D cloud-radiation interactions? , 2017 .

[113]  Robert Pincus,et al.  Computational Cost and Accuracy in Calculating Three-Dimensional Radiative Transfer: Results for New Implementations of Monte Carlo and SHDOM , 2009 .

[114]  Arthur Appel,et al.  Some techniques for shading machine renderings of solids , 1968, AFIPS Spring Joint Computing Conference.

[115]  Jan Novák,et al.  Spectral and decomposition tracking for rendering heterogeneous volumes , 2017, ACM Trans. Graph..

[116]  Céline Cornet,et al.  Three-dimensional polarized Monte Carlo atmospheric radiative transfer model (3DMCPOL): 3D effects on polarized visible reflectances of a cirrus cloud , 2010 .