/home/claudia/yang_hong_emde

Script to calculate optical properties from single scattering data: calc_opt_prop.py

Synchronize with lx321:

rsync -rdltvu yang_hong_emde emde_cl@lx001.pa.op.dlr.de:/net/lx321/export/home/emde_cl/

Single scattering data calculated by Gang Hong using the codes by Ping Yang:

/home/claudia/yang_hong_phamat/88_bands

Format of data:

Wl(um) Size(um) Volume(um3) Area(um2) Extioneff Albedo GG F-delta Theta(0-180)

(separate files for each phase matrix element)

Data needs to be integrated over size distribution.

This method has been applied to derive HEY parameterisation for libRadtran.

Effective diameters obtained for alpha=1 and reff from 5 to 95. The effective diameter obtained by integrating over Yang/Hong data is similar to the assumed reff.

From SHDOM documentation:

<latex> n®=N_0 r^\alpha \exp(-\frac{\alpha+3}{r_{\rm eff}}) </latex>

Here r is the mass equivalent radius or effective size (Yang et al., 2005)

<latex> r=\frac{3V}{2A} </latex>

where V is the volume of the particle and A is the cross section.

Evans et al. 1998 study the radiative properties (submm range) of ice clouds depending on measured size distributions of the FIREII experiment. In the measurements alpha varies from -0.5 to 7 (larger values for smaller temperatures/ smaller particles). Setting alpha=1 was ok. to model radiative properties correctly.

We now keep alpha=1 and iteratively find the parameters lambda, that result in the desired effective radii (5,10,15,…,90). The resulting optical properties (found by integration over the gamma distributions) are similar to the ones found by integration over the 30 measured size distributions (see below). The main differences are found in the asymmetry parameter, which is systematically larger for the gamma distributions than for the measured distributions, indicating that the gamma distributions underestimate the amount of small particles (small asymmetry parameter). But in all cases the optical properties obtained from the gamma distributions are in the range of the measured distributions.

Optical properties at 1.85 microns. The markers show the optical properties obtained from the size distributions, the blue line shows the fit through these points (Key parameterization), and the red line shows results from the simple gamma distribution. To compute bext, omega, and g a gamma size distribution seems to be sufficient. Only some differences in the asymmetry parameter can occur, because the gamma distribution underestimates the amount of small particles.

Effective diameters, black line corresponds to Heymsfield fit and colored lines to different shapes and Yang/Hong data.

The Gamma-distributions by Heymsfield et al. (2002) that have been fitted to in-situ measurements of the TRMM field campaigns are defined similarly to Evans et al. (1998).

<m 12> N = N_0 D^\mu \exp(-\lambda D) </m>

T {\rm smaller} 27^\circ: \lambda=24.0 \exp(-0.049 T)

T \ge 27^\circ: \lambda=6.8 \exp(-0.096 T)\\[1ex]

\mu=0.13\lambda^{0.64}-2\\[1ex]

r_e=\frac{159 \Gamma(3.2+\mu)}{\Gamma(2.8+\mu)\lambda^{0.41}}

Here D is the maximum (!) dimension of the particles. More or less consistently with Evans et al. 1998, the parameter mu varies from about -1 to 5.

Applying these size distributions to the various shapes gives effective radii in the range from approx. 20 to 100 microns, for droxtals only effective radii larger than about 70 micron which shows that probably no droxtals have been measured. As stated, small particles were not fitted. Comparing the obtained effective radii to a fit of for effective radius given in Heymsfield 2002 (based on measured IWC and A) gives a plausible result if one assumes a particle mixture (droxtals excluded).

The 30 size distributions cover the following ranges of Reff:

It is possible to get a fit for the optical properties extinction, single scattering albedo, and the asymmetrie parameter as a function of effective radius:

Optical properties fit for solid columns at 450nm using the 30 size distributions (third order polynomial fit as in Key (2002).

But for the phase functions it is problematic. In the forward scattering region, the value of the phase function extends over several orders of magnitude. The fit returns negative values for the small Reff values. These negative values can be avoided by fixing the phase function of the smallest Reff to the value obtained by the size distribution corresponding to the smallest Reff. But this results in ugly discontinuities for some scattering angles and a strange behaviour in the forward scattering region. The maximum is not at 0 degrees scattering angle but slightly larger. Therefore it is not a good idea to fit each angle in the phase function separately.

Phase function fits for solid columns (450 nm). The first and the last point (smalles and larges Reff) are computed directly from a size distribution, the values in between are fitted using a polynominal of the 3rd order.

Phase functions obtained for different effective radii from the fit above, here for aggregates and at 1.35 microns wavelength. There is a strange behaviour in the forward peak which is not visible when the phase functions are averaged over a size sistribution, therefore this behaviour is an artefact from the fit.

Another possiblity is to use the measured size distribution and fit them to gamma distributions. It turned out that most of the distributions can only be fitted by bimodal gamma distribution. The fit of the gamma distribution parameters against reff is also not very nice (large scattering). Results in /yang_hong_emde/30_dist.

Parameters of gamma distributions fitted to the 30 size distributions. The effective radius corresoponds to the one calculated using solid columns.

Ivanova (2001) has generated a parameterization for ice cloud bimodal size-distribution. Here the parameters depend not only on Reff, but also on IWC, T, particle falling speed … My conclusion from the study with the 30 size distributions is that it is probably not possible to generate a realistic parameterization for the size distribution that depends only on the effective radius.