Microwave Brightness of Polar Firn as Measured by Nimbus 5 and 6 ESMR

Abstract

The microwave emission from a half-space medium characterized by coordinate dependent scattering and absorbing centers has been calculated by numerically solving the radiative transfer equation by the method of invariant imbedding. A Mie scattering phase function and surface polarization have been included in the calculation. Also included are the physical temperature profile and the temperature variation of the index of refraction for ice. Using published values of grain-size and temperature-profile data of polar firn, the brightness temperature has been calculated for the 1.55 cm and 0.8 cm wavelengths. For selected regions in Greenland and Antarctica, the results of our calculations are in reasonable agreement with the observed Nimbus-5 and Nimbus-6 ESMR data.

Résumé

L’émission en micro-ondes issue d’un milieu hémispatial caractérisé par des centres coordonnés dépendants de dispersion et d’absorption a été calculée par résolution numérique de l’équation de transfert de chalcur avec la méthode de l’invariant inclus. On a pris en compte dans le calcul une fonction pour la phase de dispersion Mie et une polarisation en surface. On a également tenu compte du profil physique des températures et de la variation avec la température de l’indice de réfraction de la glace. En utilisant les valeurs publiées pour la tailie des grains et les données des profils de température du névé polaire, la température de brillance a été calculée pour les longueurs d’onde de 1,55 cm et de 0,8 cm. Pour des régions sélectionnées du Groenland et de l’Antarctique, les résultats de nos calculs sont en accord raisonnable avec les données observées par ESMR de Nimbus-5 et Nimbus-6.

Zusammenfassung

Durch numerische Lösung der Gleichung für die Strahlungsübertragung mit der Methode der invarianten Einbettung wurde die Mikrowellen-Emission im Halbraum über einem Medium, das durch Zentren mit ortsabhängiger Streuung und Absorption charakterisiert ist, berechnet. In die Rechnung wurde eine Mie’sche Phasenfunktion der Streuung und eine Polarisation an der Oberfläche einbezogen. Desgleichen wurde das physikalische Temperaturprofil und die temperaturabhängige Schwankung des Refraktionsindex für Eis berüchsichtigt. Auf der Basis veröffentlichter Werte für die Korngrösse und das Temperaturprofil in polarem Firn wurde die Farbtemperatur für Wellenlängen von 1,55 cm und 0,8 cm berechnet. Für ausgewählte Gebiete in Grönland und Antarktika stehen die Rechenergebnisse in guter Übereinstimmung zu beobachteten Werten der Satelliten Nimbus-5 und Nimbus-6 ESMR.

Introduction

Brightness temperatures obtained from the electrically scanning microwave radiometers (ESMR) on board the Nimbus-5 and Nimbus-6 spacecrafts (ESMR-5, operating at a wavelength of 1.55 cm and ESMR-6, at 0.8 cm, respectively) over Greenland and Antarctica (Gloersen and others, 1974, and this paper) have shown a lack of correlation with the physical surface temperature. Model studies of microwave emission from a half-space with scattering and absorbing centers were carried out by England (1974) and by Chang and others (1976). They calculated the brightness temperature for a model snow field consisting of randomly-spaced ice spheres of different radii and dielectric properties. The scattering and extinction cross-sections were calculated using the Mie scattering theory and the brightness temperature values were then obtained by numerically solving the radiative transfer equation. Results of the calculation showed that the emerging microwave radiation originates deep within the medium. These results not only provided an explanation for the lack of correlation of the observed brightness temperatures and the physical surface temperatures but also opened up the possibility of remotely sensing such parameters as the snow accumulation rate and snow temperature profile.

In the model calculations performed by Chang and others (1976) and by England (1974 and 1975), the scattering and absorption properties have been assumed to be independent of the snow depth. This is in contrast to the actual situation where the crystal size does vary with the snow depth (Gow, 1969, 1971). Also in the former work the physical temperature of snow has been taken to be independent of depth and in the latter calculation a linear temperature variation has been used. The accuracy of the numerical results obtained by England (1975) is difficult to assess because of the convergence problem associated with his method of solution (England, 1974). By using approximate, integrable analytic expressions for the radiative transfer equation, Zwally (1977) has obtained microwave brightness temperatures for polar ice sheet, taking into account the grain size and temperature profile variations in season and location. However, in order to obtain agreement between theory and experiment, he found it necessary to lower the calculated value of scattering coefficient by an order of magnitude. Here we have succeeded in providing an explanation for the radio brightness temperatures observed over Greenland and Antarctica with the ESMR-5 and ESMR-6 without the need for such an adjustment to the scattering coefficient. The expression for the variation of snow grain size with depth and location used in the calculation are those compiled by Zwally (1977) by fitting the crystal size data measured by Gow (1969, 1971) at different snow depths. The radiative transfer equation (Equation (1)), has been solved numerically using these snow parameters by the method of invariant imbedding (Chandrasekhar, 1950, p. 161–82; Redheffer, 1962; Bellman and others, 1963, p. 5–35; Preisendorfer, 1965; Grant and Hunt, 1969).

Numerical solution for the radiative transfer equation

The radiative transfer equation for an axially symmetric inhomogeneous medium in which all interactions are linear can be written in the form of an integro-differential equation (Grant and Hunt, 1969)

(1)

where the radiation intensity I(x, µ) is at depth x traveling in the direction making an angle whose cosine is µ with the normal toward the direction of increasing x (Fig. 1).

Fig 1 Radiation intensity of 1(x, μ) is at depth × traveling in a direction making an angle whose cosine is μ with the normal toward the direction of increasing x.

The functions σ(x), ω(x), B(x), and ρ(x, µ, µ′) are prescribed functions of their arguments. They are referred to as the extinction per unit length, the single scattering albedo, the source, and the phase function, respectively. For a non-uniform medium these functions are generally piecewise continuous functions of depth subject to the conditions

(2)

In the present work, the following normalization for the phase function will be used

(3)

for all values of x. Instead of working with depth x, one generally works with a dimensionless depth variable called optical depth τ, defined in differential form as

(4)

In terms of optical depth, Equation (1) reduces to

(5)

This equation of radiative transfer was solved numerically by the invariant imbedding technique (Grant and Hunt, 1969). This technique is based upon the principle that the radiation emerging from a semi-infinite plane parallel medium is invariant with respect to addition (or subtraction) of layers of arbitrary thickness to (or from) the medium. Expressed numerically, this technique solves the radiative transfer equation via matrix recurrence relations. Then by using the proper scattering phase function and the boundary conditions (England, 1974, 1975), the brightness temperature emerging from the snow field can be calculated.

Brightness temperature of the polar firn

The calculation of the brightness temperature of the polar firn requires the specification of a single scattering albedo ω(τ), physical temperature within the medium B(τ), the scattering phase matrix ρ(τ, µ, µ′), and the reflection coefficient at the medium boundary. In the single-particle scattering approach, the albedo can be calculated from the knowledge of the scattering and the extinction coefficients of the particle. These coefficients depend upon the radius of the particle and its refractive index. For ice particles, with which we are concerned, the specification of size is sufficient to determine the scattering coefficient. The extinction coefficient of the particle on the other hand depends on the imaginary part of the index of refraction which is less sensitive to the crystal size. Calculations of Chang and others (1976) show that for ice at 273 K the extinction coefficient is about an order of magnitude higher than for ice at 253 K. Although this sensitivity can be used to distinguish temperature profiles, it also introduces some uncertainty in quantitative comparison with remotely-measured brightness temperature. It is also to be noted that the imaginary part of the index of refraction changes considerably with the presence of impurities (Hoekstra and Cappillino, 1971). For continental ice sheets the effect of impurities may be negligible.

Zwally (1977) has performed a regression study of the crystal-size data for different snow depths obtained by Gow (1969, 1971) at various locations in Greenland and Antarctica. In Table I we show the crystal-size profiles, for which we have calculated the brightness temperature.

Table I. Location and mean annual surface temperature and mean crystal-size profile (Zwally, 1977)

We have mentioned that absorption per unit length is quite insensitive to the crystal size and depends largely on the imaginary part of the refractive index n″. Although the choice of n″ is not completely arbitrary, its exact value is difficult to obtain. In this study we have developed an interpolation/extrapolation algorithm to determine the index of refraction of ice from experimental data (Evans, 1965, fig. 7). Table II shows the experimental data for n″ and the interpolated and extrapolated values for 1.55 cm and 0.8 cm wavelengths. The interpolated values for 3.2 cm differ from the measured values reported by Cumming (1952). These differences are probably caused by the structure of the different ice types tested.

Table II. The imaginary part index of refraction for ice at different wavelengths and temperatures

The single-scattering albedo is defined as

(6)

Therefore, the depth variation of ω depends on the depth variation of γ _s and γ _a. Since γ _a varies with the physical temperature (Cumming, 1952) of the snow, an analytic expression for depth and seasonal variation of polar snow temperature has been used (Zwally, 1977).

(7)

The maximum surface temperature occurs at time t = 0 and t = 365 d and minimum temperature is for time t = (365/2) d. The mean temperature of the surface and the asymptotic temperature is 250 K and the peak-to-peak variation at the surface is 30 K. This mean temperature is quite acceptable for the locations given in Table I.

The correct value of n″ is crucial because the absorption coefficient γ _a depends on this value. To illustrate the effect of different values of n″, a single calculation has been performed. The Camp Century, Greenland, site was chosen for the comparison. Figure 2 shows the variation of vertically-polarized brightness temperature by fixing the n″ in each calculation, while the other parameters were remained the same as in the other calculations.

Fig. 2 Dependence of 0.8 cm vertically polarized brightness temperature on n″ the imaginary part of index of refraction.

Calculated brightness-temperature values for several locations and times are given in Table III. The agreement between the calculated and the observed ESMR-5 values are considered quite good. The r.m.s. deviation between the calculated brightness temperature and the observed brightness temperature for ESMR-5 is 8.4 K. The agreement between the calculated brightness temperature and the observed ESMR-6 values is not as good as the ESMR-5, with r.m.s. deviation of 14.5 K and 20.4 K for ESMR-6 vertical and horizontal polarization respectively. The differences in r.m.s. deviation for the calculated ESMR-5 and ESMR-6 values are attributable at least in part to the following factors: (1) the observed seasonal variations in the brightness temperature larger than those calculated indicate that the microwave emission for 0.8 cm wavelength largely emanates from a thin surface layer and thus the values assumed for n′ at this wavelength were too low. As a result, the calculated extinction coefficient is too small, resulting in too large a penetration depth, (2) sensitivity to near-surface moisture due to solar heating of snow, in the presence of below-freezing air temperature, and (3) the inherent difficulties in the calibration of ESMR-6 due to the spacecraft thermal conditions.

Table III. Comparison of observed and calculated brightness temperatures of locations in Table I

Summary and conclusions

A microscopic single-particle scattering model has been used along with the measured variation of crystal size and temperature with depth to provide a quantitative explanation of observed brightness temperatures of the South Polar and North Polar regions. This extends the calculations performed by Chang and others (1976) by taking into account the variations of grain size and temperature with depth. The calculated brightness temperatures for 1.55 cm wavelength are in good agreement with the observations, but one of the calculated brightness temperatures for 0.8 cm wavelength differs by 32 K with the observations. The probable causes of these discrepancies are: (1) uncertainty in the calibration of the ESMR-6 instrument, and (2) uncertainty in the values of n′. In some areas of Greenland, percolation is known to occur, which is accompanied by a large change in n′. Since we felt that sufficiently accurate values of n′ in this region are not available, we did not consider these areas in our comparisons.

Based on the single-particle scattering model, the calculated brightness temperatures correspond to the data obtained from the ESMR 5 and ESMR 6 over selected sites in Greenland and Antarctica. Using essentially the same temperature and grain-size profile that were used by Zwally (1977), we find that our calculated brightness temperatures for 1.55 cm wavelength are significantly higher than his. This discrepancy is probably due to the approximate analytic solution of the radiative transfer equation used in his calculations of the brightness-temperature values. In view of these findings, our calculations would provide a more realistic brightness-temperature estimate for the polar firn.

Apart from the independent-particle scattering model discussed above, there are alternate explanations for the source of microwave scattering and absorption for radio brightness temperatures of Antarctica and continental glaciers which have been studied by Gurvich and others (1973), Stogryn (1974), and Tsang and Kong (1976). In these studies, the source of microwave scattering is the fluctuation of the dielectric constant of the media (Tatarskiy, 1959, English translation, p. 59–70). A special case of such fluctuations occurring where freeze-thaw cycles have taken place, accompanied by percolation of melt water and formation of sub-surface ice lenses, has been identified as important particularly at the longer wavelengths where single-particle scattering is of lesser importance (private communication from P. Gudmundsen in 1978). The calculated brightness temperature based on these models, depends upon the variance and the correlation length of the fluctuation of the dielectric media. Physically this scattering mechanism is as plausible as scattering by independent ice grains, and it is reasonable to believe that these two scattering mechanisms co-exist in the medium and should be considered together in the quantitative calculation.