This is a proposal to probe local part of the interplanetary dust (IPD) cloud complex and retrieve mean volume emissivity of the local IPDs at mid-infrared wavelengths. This will be done by monitoring, with Infrared Camera (IRC) aboard the ASTRO-F, the annual modulation of the zodiacal emission. In pointing mode of the ASTRO-F mission the spacecraft can make attitude maneuvering over approximately ±1 ̊ range centered at solar elongation 90 ̊ in the ecliptic plane. The attitude maneuvering combined with high sensitivity of the IRC will provide us with a unique opportunity observationally to take derivatives of the zodiacal emission brightness with respect to the solar elongation. From the resulting differential of the brightness over the ±1̊ range, one can directly determine the mean volume emissivity of the local IPDs with a sufficient accuracy to de-modulate the annual emissivity variations due to the Earth's elliptical motion and the dis-alignment of the maximum IPD density plane with respect to the ecliptic. The non-zero eccentricity (e⊕= 0.0167) of the Earth's orbit combined with the sensitive temperature dependence of the Planck function would bring modulations of amplitude at least 3.34% to the zodiacal emission brightness at mid-infrared wavelengths, with which one may determine the IPD temperature T(r) and mean number density n(r) as functions of heliocentric distance r. This will in turn fix the power-law exponent δ in the relation T(r) = T_o(r/r_o)-δ for the dust temperature and v in n(r) = n_o(r/r_o)-v for the density. We discuss how one may de-couple the notorious degeneracy of cross-section, density, reference temperature To and exponent δ.

We have written a code called QDM_sca, which numerically solves the problem of radiative transfer in an anisotropically scattering, spherical atmosphere. First we formulate the problem as a second order differential equation of a quasi-diffusion type. We then apply a three-point finite differencing to the resulting differential equation and transform it to a tri-diagonal system of simultaneous linear equations. After boundary conditions are implemented in the tri-diagonal system, the QDM_sca radiative code fixes the field of specific intensity at every point in the atmosphere. As an application example, we used the code to calculate the brightness of atmospheric diffuse light(ADL) as a function of zenith distance, which plays a pivotal role in reducing the zodiacal light brightness from night sky observations. On the basis of this ADL calculation, frequent uses of effective extinction optical depth have been fully justified in correcting the atmospheric extinction for such extended sources as zodiacal light, integrated starlight and diffuse galactic light. The code will be available on request.