Monte Carlo simulations of light scattering by composite particles in a planetary surface

Composite particles containing internal scatterers have been proposed as an explanation for the fact that most photometric studies of planetary surfac es based on Hapke's bidirectional reflectance model have found the planetar y particles to exhibit moderately backscattering phase functions. However, an implicit assumption made in this explanation is that the scattering by c omposite particles containing multiple internal inclusions in a planetary s urface can still be adequately computed using standard radiative transfer t heory assuming the composite particles to be the fundamental scattering uni t even though the particles are necessarily in close proximity to each othe r. This assumption was explored by J, K, Hillier (1997, Icarus 130, 328-335 ) using a Monte Carlo routine. However, in this initial study several simpl ifying assumptions were made. The internal scatterers were assumed to be is otropic and scattering off of the surface and absorption within the composi te particle were ignored. While these assumptions are not very realistic, i t was believed that the study could still provide insight into the light sc attering by such surfaces, Here we relax these assumptions in order to exam ine the light scattering by more realistic particles. Almost all of the con clusions reached in the earlier paper remain valid, As before, we find that classical radiative transfer (assuming a random distribution of scattering particles) coupled with the assumption that the composite particle is the fundamental scatterer provides a good approximation in the high porosity li mit, However, even for porosities as high as 90% the effects of close packi ng are clearly seen with the radiative transfer calculation underestimating the scattering by similar to 10% at high phase angles. In contrast to the earlier study we find that the radiative transfer calculation tends to over estimate, not underestimate, the scattering at high emission but moderate p hase angles. As the porosity is lowered further, the discrepancy becomes mo re severe and can reach 100% or greater. In particular, our main conclusion remains intact: the parameters derived using the classical radiative trans fer theory will yield results intermediate between those of the composite a s a whole and those of the internal scatterers and thus one should exercise caution in interpreting the results of models based on classical radiative transfer theory in terms of the physical properties of the surface particl es. (C) 2001 Academic Press.