AVHRR surface temperature and narrow-band albedo comparison with ground measurements for the Greenland ice sheet

Abstract

An ice-surface temperature retrieval algorithm for the Greenland ice sheet was developed using NOAA 11 thermal radiances from channels 4 and 5. Temperature, pressure and humidity profiles, cloud observations and skin temperatures from the Swiss Federal Institute of Technology (ETH) camp, located at the equilibrium line altitude at 49°17′W, 69°34′N, were used in the LOWTRAN 7 model. Through a statistical analysis of daily clear sky profiles, the coefficients that correct for the atmospheric effects were determined for the ETH-Camp field season (May to August). Surface temperatures retrieved by this method were then compared against the in situ observations with a maximum difference of 0.6 K. The NOAA 11 narrow-band planetary albedo values for channels 1 and 2 were calculated using pre-launch calibration coefficients. Scattering and absorption by the atmosphere were modelled with LOWTRAN 7. Then, narrow-band albedo values for the AVHRR visible and near infrared channels were compared with in situ high resolution spectral reflectance measurements. In the visible band (580–680 nm), AVHRR-derived narrow-band albedo and the in situ measurements corrected with radiative transfer model LOWTRAN 7 showed a difference of less than 2%. For the near infrared channel (725–1100 nm) the difference between the measured and modelled narrow-band albedo was 14%. These discrepancies could be either the result of inaccurate aerosol scattering modelling (lack of the in situ observation), or the result of sensor drift due to degradation.

Introduction

Past and ongoing research indicated that several polar surface properties can be mapped through a combination of multispectral satellite data. These properties include spectrally integrated surface albedo, surface temperature and short- and long-wave radiation balance, which in turn is an important measure for the discussion of the long-term climate stability. The Greenland ice sheet provides an ideal case for the calibration of satellite-derived parameters due to its large size and différent homogenous snow and ice surfaces over large areas. Further, the study of seasonal and interannual variations of the above mentioned parameters provide an ideal data set to determine the “noise level” for future climate trend studies. Here the technique and some preliminary results are presented to estimate ice surface temperatures and surface albedo values from different spectral channels of the NOAA AVHRR.

Data and Method

AVHRR calibration and navigation

The advanced very high resolution radiometer (AVHRR) data used in this study were obtained from National Océanographie and Atmospheric Research Laboratories (NOARL). This is the local area coverage (LAC) data in level lb format from the NOAA-11 satellite, and includes 5 channels: 2 visible and 3 infrared. The level lb format includes information for the navigation and calibration of the images (Lauritson and others, 1979), but this information was re-done on site to gain more control over the processing of the data. The navigation of the AVHRR data can be used to better depict the actual motion of the satellite (Baldwin and Emery, 1993). A small registration error is present in the resulting images that are removed by linearly nudging the images onto a map grid. Images from each of the days are are then compared with each other. Registration errors can then be reduced to ±½ pixel. The navigation of the data was mapped to a polar Stereographic projection with a 1 km per pixel resolution. The resolution of the scan spot from the AVHRR instrument varies from 1.1 km at nadir to about 6.8 km at a scan angle of 55°. This is not a linear function; a scan angle of 30° still results in a pixel resolution of about 1.5 km.

Calibration of the AVHRR imagery is channel dependent. Calibration of the reflected solar energy “visible” channels, 1 and 2, is done by using pre-launch values particular to each satellite. These values (Kidwell, 1991) convert the instrument counts into per cent albedo. Solar zenith angles are computed during the navigation procedure for each pixel in the image and corrections made to the per cent albedo. Channel 3 is a combination of reflected and emitted energy during daylight conditions and was not used in this paper. The emitted infrared channels, 4 and 5, require several steps for calibration. First, continuous on-board sensing of space and black-body targets every few scan lines provides information for computing slope and intercept values for the infrared channels (Lauritson and others, 1979). Energy values are calculated, and then using Planck’s equation with pre-launch wavenumbers, brightness temperatures can be found. Non-linear corrections that are dependent upon the scene temperature are then applied. These non-linear corrections were on the order of −1.0° for the observed area.

Ice-surface Temperature

A common approach for estimating surface temperatures is to relate satellite data to surface temperature observations with a regression model. The radiative transfer equation can be applied to model satellite radiances and brightness temperatures. A common approach to retrieve surface temperatures is to use the split-window infrared AVHRR channels 4 and 5 (10.5–11.5 μm, 11.5–12.5 μm) as was applied for open-ocean areas (McClain and others, 1985; Minnett, 1990) as well as sea-ice areas (Key and Haefliger, 1992). To correct for atmospheric attenuation of the satellite measured brightness temperatures over the Greenland ice sheet the surface temperatures are regressed against modelled AVHRR brightness temperatures as described by Key and Haefliger (1992). Atmospheric attenuation was derived with the radiative transfer model LOWTRAN 7 (Kneizys and others, 1988) using radiosonde profile measurements from the ETH-Camp (Federal Institute of Technology (ETH), Department of Geography (Fig. 1)) as input data. A total of 78 radiosondes were launched from 8 June through 29 August 1990, and 194 from 13 May through 21 August 1991. The profiles include pressure, temperature and relative humidity for every second up to a height of 20–25 km (Fig. 2). Continous incoming and outgoing longwave radiation measurements and daily cloud observations (7 times each day) were also available from the ETH Camp. The longwave radiation at the ETH Camp was calculated as the difference between the all-wave radiation, measured by pyrradiometers, and the short-wave radiation, measured by pyranometers, and compensated for by the emission-loss from the instrument. The pyrradiometers were calibrated over a melting ice surface. Only situations with less than ⅜ cloud cover were used in the analysis, which are referred to as “clear sky” conditions. A total of 90 profiles remained with clear sky conditions.

Fig. 1. Overview of Disco Bay in West Greenland (560 × 560km). The rectangular area around the ETH Camp (152 × 227 km) is the test area to which the coefficients are applied (dotted lines = coast, solid line = ice margin).

Fig. 2. Radiosonde profile measurements of temperature, relative humidity, and wind speed from the ETH Camp, located at the equilibrium line altitude on the Greenland ice sheet (49°17′ w, 69°34′ N).

The surface temperature at the camp is calculated from the outgoing longwave irradiance with the Stefan-Boltzmann law: R = σT ⁴, assuming the snow surface is a black-body. In general the longwave irradiance measurements have an error of ± 10Wm⁻² (±3%) or even better in case of the Greenland measurements (personal communication from T. Konzelmann). To every measured outgoing longwave irradiance 3% are added and subtracted, respectively. This way three new profiles are created out of each of the 90 original clear sky ones, where the surface temperatures are the minimum

and the maximum

calculated temperatures, respectively. This results in a sample size of 270 clear sky profiles.

To eliminate the influence of the atmosphere on the signal reaching the satellite, the measurements of at least two spectral bands are required. A good overview of the theoretical algorithms for satellite-derived SSTs is presented in Barton and others (1989). The method used in this study is the same as described in Key and Haefliger (1992) and reads as following:

(1)

T ₄ and T ₅ are the satellite-measured brightness temperatures in Kelvin and θ is the sensor scan angle (0–55°). The coefficients are determined through a least-squares regression procedure. Surface temperatures are regressed against modelled brightness temperatures. LOWTRAN 7 is used to simulate radiances in the AVHRR thermal channels 4 and 5, which are then converted to brightness temperature.

The coefficients for the estimation of ice surface temperature (1ST) according to Equation (1) are shown in Table 1. The root mean square (RMS) error of 0.3 K is well within the error of the outgoing longwave measurements which is ±1 K for the Greenland measurements. Table 2 shows the calculated 1ST for the AVHRR pixel representing the location of the ETH Camp. Due to the uncertainty of the exact location (AVHRR geolocation error), the temperature of 3 × 3 pixels are calculated and compared to the in situ measurement of the camp. The comparison of AVHRR-derived ice surface temperatures with the in situ measurements show differences of up to 0.6 K. The regional variation of ice surface temperatures around the ETH Camp are shown in Figure 3 (see also Figure 1 for overview of area).

Fig. 3. Ice-surface temperatures around the ETH camp (see Fig. 2). The variation from day to day is shown. The image of 26 June contains an error in a few scan lines, visible in the lower right corner.

Table I. Coefficients and RMS error based on NOAA 11 AVHRR channels 4 and 5

Table 2. Measured and calculated AVHRR ice-surface temperature (1ST) at the ETH Camp

Planetary and surface narrow-band albedo

In the presented case study, the planetary narrow-band albedo for AVHRR channels 1 and 2, measured on 23 May 1991, will be compared against ground-based values measured with a portable spectrometer. The Geophysical Research Corporation (GER) portable field spectrometer measures reflectance values for 875 different spectral bands between 300 and 2500 nm (Table 3). The spectrometer was regulated with a laptop computer (SNAP 1 + 1) which was also used for data acquisition. The field-of-view of the optical head is 7° in horizontal and 1.6° in vertical direction (rectangular area of sensitivity given by the photodetector array). An integration sphere was mounted in front of the optical head for hemispheric measurements which could be directed towards the zenith or perpendicular to the ground which was used to measure the spectral reflectance of the snow surface. The spectral reflectance of the snow surface is shown in Figure 4. The snow was fresh, less than one day old, with grain sizes between 0.1–0.2 mm in diameter. The prevailing air temperature was well below freezing. In addition to the spectrometer measurements broadband global radiation, and diffuse and reflected shortwave radiation were measured with a set of pyranometers as part of the ongoing radiation balance experiment at the ETH Camp. The reflected shortwave radiation was measured at 1.5 m and at 27 m above the surface. The broadband albedo values derived from the pyranometers were 84.3% for the lower instrument, and 82.2% for the instrument at 27 m. The slightly lower albedo value for the later can be explained by the larger field-of-view which included part of the field camp as well as surface roughness in the scale order of meters and larger. The spectrally-integrated reflectance derived from the spectrometer is 82.9% which shows the level of accuracy (δ = 0.7%) which can be expected using two different kind of instruments.

Fig. 4. Spectral reflectance of new snow (grain size 0.1–0.2 mm), pre-mell period.

Table 3. specification for the GER-spectrometer

The narrow-band albedo measured by the AVHRR sensors (i = 1,2) is computed as a linear function of the input value C and the pre-launch calibration coefficients as follows (Kidwell, 1991):

(2)

with S(1) = 0.095, S(2) = 0.1061,I(1) = −3.8,I(2) = −3.6 for NOAA 11.

This calibration takes the spectral response of the channels (i = 1,2) into account by multiplying the spectral irradiance viewed at the top of the atmosphere with the spectral response function. NOAA 11 satellite data from the ETH Camp area were collected on 23 May 1991, 15.11 GMT with a satellite zenith angle relative to the camp of 27.37° and a solar zenith angle of 49.16°. The narrow-band planetary albedo values for AVHRR channels 1 and 2 are listed in Table 4. The accuracy of the AVHRR-derived narrow-band albedo values are discussed in the following.

Table 4. Narrow-band albedo values for AVHRR channels 1 and 2 as measured by A VHRR at the top of the atmosphere (a_p — planetary narrow-band albedo), and calculated using ground-based spectrometer values a_g and radiative transfer modelling

The atmospheric absorption and scattering was modelled using the radiative transfer code LOWTRAN 7. Two model runs were performed: (a) using the in situ radiosonde profile measurements of temperature and humidity as input data, and (b) using the standard subarctic atmosphere. Figure 5 depicts the solar irradiance at the top of the atmosphere and the irradiance as modelled with LOWTRAN 7 at the surface using radiosonde profile measurements. The narrow-band albedo at the top of the atmosphere dp was then calculated by integrating over the AVHRR spectral band-width (λ₁, λ₂), as following:

(3)

where S(λ)is the solar irradiance at the top of the atmosphere, ∅_i(λ)is the spectral response function of the AVHRR channel (i = 1,2), a_g (λ) is the narrow-band reflectance at the ground, and τ(λ) is the calculated spectral transmission for two path lengths: (a) between the top of the atmosphere and the surface, and (b) between the surface and the satellite. The integration was done for both AVHRR channels with (a) radiosonde profile data and (b) standard sub-Arctic atmosphere (Table 4). The difference between measured and calculated narrow-band planetary albedo for channel 1 is 0.5% using radiosonde profile data in the LOWTRAN 7 model. For the standard sub-Arctic atmosphere, the difference increased to 2.9%. However, for AVHRR channel 2 the calculated and measured narrow-band albedo values difler by approximately 9%. Surface narrow-band albedo a_g can also be derived from planetary narrow-band albedo a_p values as following:

Fig. 5. Solar irradiance outside the atmosphere, and at the surface as modelled by LOWTRAN 7 radiative transfer model with the in situ radiosonde data for 23 May, 15.11 GMT.

Table 5. Narrow-band surface albedo values a_g compared against calculated surface narrow-band albedo from AVHRR channels 1 and 2 using transmission coefficients τ₁ (sun-surface), and τ₂ (surface-satellite)

(4)

The transmission values were calculated with the radiative transfer model which takes the intrinsic atmospheric reflectance due to the Rayleigh and aerosol scattering of the standard sub-Arctic atmosphere into account. The surface narrow-band albedo values derived from the AVHRR planetary albedo are given in Table 5. The differences have the same order of magnitude for the calculated planetary narrow-band albedo which is not surprising as both computations use the LOWTRAN 7 model with identical input parameters.

So far, we have neglected the bidirectional reflectance (BR) of a snow surface as a possible error due to large scanner angle and large solar zenith angle (Steffen, 1987). It can be assumed that the snow has an isotropic reflectance for the small satellite zenith angle encountered during this study. For satellite zenith angles of 50° or larger, the BR of snow has to be taken into account, especially for wet snow surface. The in situ measurements of aerosol distribution with height were not available. This might be an additional error source. Aerosol content can vary with time and space significantly which can reduce the reflectance by several per cents. If we assume that the radiative transfer model is correct, we can speculate that a major sensor drift of AVHRR channel 2 has occurred since launch of NOAA 11 on 8 November 1988. For the visible channel of NOAA 9 degrading rates of up to 6% per year were found (Staylor, 1990). This case study should be judged as a preliminary analysis, as more cases have to be compared to get a statistically significant error analysis. Further, we have to find a means to validate the LOWTRAN 7 derived spectral transmission values, in order to separate the modelling error from the error caused by the sensor drift.

Summary and Conclusion

The relationship between clear sky NOAA 11 thermal radiances and the measured surface temperature of the Greenland ice sheet is examined through forward calculations of radiative transfer equation. Pressure, temperature and relative humidity profiles and cloud observation from the ETH Camp are used. Using the split window channels 4 and 5 and scan angle, the RMS error in the estimated 1ST is 0.3 K.

Narrow-band surface albedo values can be esümated from AVHRR channel 1 and 2 measurements with an accuracy of 2% and 14% respectively, using LOWTRAN 7 radiative transfer modelling and radiosonde profile information. More case studies are needed to discuss this technique in more depth and to derive statistically significant values.

Acknowledgements

This research was funded by NASA Ocean Science Branch (Grant NAG 5–882), and by the Swiss National Science Foundation (Grants 20–2571.88 and 21–27449.89). Thanks are due to the ETH expedition team for collection of the radiation and radiosonde data, and to A. Ohmura and T. Konzelmann, Swiss Federal Institute of Technology, and to J. Key, University of Colorado, for useful discussions.