315 Extrapolation of Simulated Scattering Matrix for Large Size Parameter Values

Wednesday, 9 July 2014
Teemu Mäkinen, Finnish Meteorological Institute, Helsinki, Finland; and T. Nousiainen

Handout (454.8 kB)

In order to compute the radiative effects due to atmospheric aerosol particles, their single-scattering properties are needed. For non-spherical, inhomogeneous particles such as mineral dust, solving these properties accurately can be extremely challenging. For particles clearly larger than the wavelength, the computation of these properties require considerable simplifications to the characteristics of the model particles, the physics in the computational method, or both.

In a recent paper[1] light scattering by single, inhomogeneous dust particles was simulated using the discrete-dipole approximation, based on shapes and compositions derived directly from measurements of real particles. The components of the scattering matrix S as a function of the scattering angle θ were calculated for size parameters x = 2πaeqλ from 0.5 to 16. Computational cost makes it infeasible to use the approach directly to size parameters much larger than that. Instead of resorting to simplifications to deal with larger particle size parameters, we propose an alternative method based on extrapolation, which may prove to be preferable at least in some cases. A simple extrapolation scheme is described, providing a reasonable solution to extending the range of available values. Within the scope of this abstract, we restrict ourselves to discussing the scheme for component S11 of the matrix.

A stable extrapolation can only be obtained if the asymptotic behavior of a function is known. Ignoring the extinction effect, the relation dθS11(x,θ) x2 holds for large x. However, while the behavior of S11 for side scattering angles can to a reasonable degree be separated as S11(x,θ) = f(x)g(θ), the prominent back- and forward-scattering peaks only become asymptotically separated in a different set of coordinates, as S11(x,θ) = f(xcos(θ))g(xsin(θ)), necessitating the separate modeling of the peaks.

Here, the extrapolating size parameter functions f(x) for all features were defined as parametrizable hyperbolic curves with the property f(x) 2 as x . The asymptotic shape of g(θ) and the parameters of f(x) for the side-scattering model were represented by spline curves. The back- and forward-scattering peaks were modeled with a series of generalized sinc functions, as ncnsincn2(qxsin(θ)) where q is a scale parameter. The parameters of this model were determined by fitting it to simulated values for size parameters between 10 and 16. This allows the evaluation of S11 for arbitrary values of x and θ as shown in Fig. 1. Since the approach is stable, an error estimate can be derived by comparing the model directly to simulated values.

A similar approach of identifying the features and modeling them separately applies to other components of the scattering matrix, as well. Furthermore, it is possible to take into account theoretical values for the x limit, essentially turning the extrapolation into interpolation.


Figure 1: Extrapolated behavior of the S11 component of the scattering matrix of a Dolomite particle. The values up to size parameter 16 (dashed line) are provided by a simulation; those values are extrapolated up to size parameter 32 with the described method.


[1]   H. Lindqvist, O. Jokinen, K. Kandler, D. Scheuvens, and T. Nousiainen. Single scattering by realistic, inhomogeneous mineral dust particles with stereogrammetric shapes. Atmos. Chem. Phys., 14:143–157, 2014.

, 14:143–157, 2014.

- Indicates paper has been withdrawn from meeting
- Indicates an Award Winner