Numerical Simulation of the Effects of Surface Roughness on Light Scattering by Hexagonal Ice Plates

Cirrus clouds have an extensive global coverage and their high altitude means they play a critical role in the atmospheric radiation balance. Hexagonal ice plates and columns are two of the most abundant species present in cirrus and yet there remains a poor understanding of how surface roughness affects the scattering of light from these particles. To advance current understanding, the scattering properties of hexagonal ice plates with varying surface roughness properties are simulated using the discrete dipole approximation and the parent beam tracer physical–optics method. The ice plates are chosen to have a volume-equivalent size parameter of 2πr/λ=60, where r is the radius of the volume-equivalent sphere, and a refractive index n=1.31+0i at a wavelength λ=0.532 µm. The surface roughness is varied in terms of a characteristic length scale and an amplitude. The particles are rotated into 96 orientations to obtain the quasi-orientation averaged scattering Mueller matrix and integrated single-scattering parameters. The study finds that the scattering is largely invariant with respect to the roughness length scale, meaning it can be characterised solely by the roughness amplitude. Increasing the amplitude is found to lead to a decrease in the asymmetry parameter. It is also shown that roughness with an amplitude much smaller than the wavelength has almost no effect on the scattering when compared with a pristine smooth plate. The parent beam tracer method shows good agreement with the discrete dipole approximation when the characteristic length scale of the roughness is several times larger than the wavelength, with a computation time reduced by a factor of approximately 500.