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%.