A thin cirrus cloud thermal infrared radiative transfer model has been developed for application to cloudy satellite data assimilation. This radiation model was constructed by combining the Optical Path Transmittance (OPTRAN) model, developed for the speedy calculation of transmittances in clear atmospheres, and a thin cirrus cloud parameterization using a number of observed ice crystal size and shape distributions. Numerical simulations show that cirrus cloudy radiances in the 800–1130-cm−1 thermal infrared window are sufficiently sensitive to variations in cirrus optical depth and ice crystal size as well as in ice crystal shape if appropriate habit distribution models are selected a priori for analysis. The parameterization model has been applied to the Atmospheric Infrared Sounder (AIRS) on board the Aqua satellite to interpret clear and thin cirrus spectra observed in the thermal infrared window. Five clear and 29 thin cirrus cases at nighttime over and near the Atmospheric Radiation Measurement program (ARM) tropical western Pacific (TWP) Manus Island and Nauru Island sites have been chosen for this study. A χ2-minimization program was employed to infer the cirrus optical depth and ice crystal size and shape from the observed AIRS spectra. Independent validation shows that the AIRS-inferred cloud parameters are consistent with those determined from collocated ground-based millimeter-wave cloud radar measurements. The coupled thin cirrus radiative transfer parameterization and OPTRAN, if combined with a reliable thin cirrus detection scheme, can be effectively used to enhance the AIRS data volume for data assimilation in numerical weather prediction models.
Satellite data assimilation in numerical weather prediction models requires an efficient and accurate radiative transfer model for the computation of radiances and Jacobians. Present thermal infrared radiative transfer models for satellite data assimilation have been developed primarily for clear conditions (i.e., pure absorbing atmospheres). However, many studies have found that a great majority of satellite observations is “contaminated” by clouds. For example, Saunders (2000) showed that only 6.5% of the High Resolution Infrared Radiometer Sounder (HIRS) footprints were clear based on the analysis of the collocated HIRS 17-km and the Advanced Very High Resolution Radiometer Global Area Coverage data. Using Atmospheric Infrared Sounder (AIRS; Aumann et al. 2003) simulated data, it was reported that less than 1% of the AIRS-observed fields of view were found to be cloud free at the instrumental noise level (Chahine et al. 2006). Of the large number of cloud-contaminated satellite observations, Wylie and Menzel (1999) reported that 34% of HIRS observations between 65°S and 65°N contained high-level clouds (both semitransparent and opaque) above 6 km. In addition, Wang et al. (1996) illustrated that 50%–70% of the Limb-sounding Stratospheric Aerosol and Gas Experiment data were contaminated by cirrus clouds.
Satellite-derived cloud climatology shows that cirrus clouds regularly cover about 20%–30% of the globe (Liou 1986). In particular, Wylie et al. (2005) reported that optically thin cirrus clouds with visible optical depths less than 1.4 were found in 20% of the HIRS data from 1979 to 2001. The effect of cirrus clouds on the energy balance of the earth–atmosphere system is a topic of critical importance because on the one hand, they affect solar radiation, referred to as the albedo effect, and on the other hand, they trap a significant amount of thermal infrared radiation emitted from the atmosphere below and the surface, referred to as the greenhouse effect. The competition between the two effects determines the gain and loss of radiative energy and hence the warming and cooling of the atmosphere. It is clear that thin cirrus clouds play an important role in global climate study. Moreover, in view of the fact that many advanced infrared and microwave sensors have been and will be deployed in space in the near future and that there appears to be a need for the incorporation of thin cirrus cloudy radiance data in numerical weather prediction models, it is important to develop a radiative transfer model for infrared satellite data simulation that can be effectively applied to thin cirrus cloudy atmosphere where upwelling radiances only marginally deviates from clear radiances. In our study, we are particularly interested in thin cirrus clouds with optical depths less than about 1.0.
We have developed an efficient radiative transfer model by combining the Optical Path Transmittance (OPTRAN) program (McMillin et al. 1995b; Kleespies et al. 2004) for the speedy computation of clear-column transmittances in numerical weather prediction models, and a thin cirrus radiative transfer parameterization based on prescribed size and habit distribution models and the associated scattering and absorption properties. We use this model to simulate infrared spectral radiances, and apply it to the AIRS data to generate clear and thin cirrus cloud brightness temperature (BT) spectra. Cirrus optical depth and cloud microphysical properties, including size and habit factors, have been inferred by a χ2-minimization program involving differences between the modeled and measured AIRS BT spectra. Moreover, we compare these inferred values to the independently derived cirrus optical depth and ice crystal size from collocated and coincident ground-based Department of Energy’s (DOE’s) Atmospheric Radiation Measurement (ARM) program millimeter-wave cloud radar (MMCR) retrievals. We have also carried out numerical experiments to study the sensitivity of the BT spectra to cirrus optical depth and ice crystal size and habit.
This paper is organized as follows. Section 2 presents the infrared radiative transfer model along with a thin cirrus parameterization. Also included are discussions of the database for ice crystal size and habit distributions, the single-scattering properties for nonspherical ice crystals, and the methodology for computing the bulk properties of cirrus clouds. Section 3 describes the calculation of clear and thin cirrus cloudy radiances, and the χ2-minimization method for inferring optical depth and ice crystal size and habit. In section 4 we show results of the sensitivity study, apply our methodology to three AIRS data granules over the ARM tropical western Pacific (TWP) sites, and compare the inferred cloud properties with the ground-based MMCR retrievals. Finally, a summary is given in section 5.
2. A thin cirrus radiative transfer parameterization
a. Radiative transfer model
Based on the principle of radiative transfer, the upwelling radiance Iν at the top of the atmosphere (TOA) in cirrus cloudy atmospheres is contributed by the radiance transmitted through the cloud from the atmosphere and surface below and the radiance emitted from the cloud. Thus, we have the following linear equation (Ou et al. 1993):
where the two terms on the right denote the transmission and emission contributions, respectively, and I0 denotes the upwelling thermal infrared radiances reaching the cloud base. Since tropical cirrus clouds locate in the upper troposphere above which temperature is cold and the water vapor amount is low, the contribution of emissions from water vapor above cirrus is small compared to the transmission through the cloud and the emission from the cloud itself, and can be neglected. Thus, I0 is simply referred to as the clear radiance in this study. Here Bν(Tc) is the cloud Planck radiance at a wavenumber ν, and the cloud temperature is denoted by Tc. The cloud emittance is denoted by ɛν, which under local thermodynamic equilibrium can be defined for remote sensing application (Minnis et al. 1998) as follows:
where μ is the cosine of the viewing angle, τIR is the infrared optical depth of a given cloud layer, ων is the single-scattering albedo, and (1 − ων)τIR is defined as the spectrally dependent cloud absorption optical depth. Equation (2) is applicable to cirrus clouds with optical depths approaching 0.
To evaluate τIR in Eq. (2), a relationship between visible and infrared optical depths is needed. Because the typical ice particle size is much larger than the visible wavelength, the mean extinction efficiency of an ensemble of ice crystals at a visible wavelength is approximately given by the geometric-optics asymptotic value of 2. The infrared optical depth can then be approximated by (Minnis et al. 1998)
where 〈Qext,IR〉 is the mean extinction efficiency at an infrared wavelength to be defined in section 2c, and τ is the cirrus cloud optical depth at a visible wavelength (e.g., 0.65 μm).
b. Ice crystal size distribution and habit distribution models
To build size distribution (SD) and habit distribution (HD) models for computing radiances in various thin cirrus conditions, we have acquired nine size distributions based on the microphysics measurements from the Central Equatorial Pacific Experiment (CEPEX) that was conducted during March and April 1993, as given in Rolland et al. (2000) and McFarquhar and Heymsfield (1996). The microphysical measurements were made by the Aeromet Learjet between 20°S and 2°N latitude and between 165°E and 170°W longitude, with operations centered at Nadi, Fiji.
We used the microphysical observations corresponding to three anvils that were carried out on 4 April, 1 April, and 17 March 1993. On 4 April, the Learjet flew through a cumulonimbus complex northeast of Fiji between 2150 and 0207 UTC. On 1 April, the airborne probes sampled a well-defined cirrus anvil associated with a dissipating nocturnal thunderstorm north of Nadi. The cirrus sampled on 17 March was the outflow of a much deeper convective disturbance. The microphysical data from these flights were expected to be representative of the ice crystal size distributions typical of tropical cirrus anvil and convective clouds. Airborne in situ instruments included the 2D cloud probe (2D-C), which sampled ice crystals with maximum dimension from 30 to 1800 μm with a resolution of 30 and 40 μm in the y and x directions, respectively, and the video ice particle sampler (VIPS), which measured smaller ice crystals with a low detection limit of approximately 5–10 μm.
To characterize the ice crystal SD, a mean effective size De is defined as 1.5 times the ratio of ensemble particle volume and the particle projected area as follows (Foot 1988; Francis et al. 1994; Wyser and Yang 1998):
where l and k are indices for HD and SD, ranging from 1 to 11 and 1 to 9, respectively. Here Vi(L) and Ai(L) are the volume and projected area of an ice particle with a maximum dimension L for the ith habit, respectively; fi,l(L) is the fraction of the ith habit in the lth HD; nk(L)dL is the number density of ice crystals in the size interval (L, L+dL) for the kth SD; Lmin k and Lmax k are the minimum and maximum L for the kth SD, respectively; and Vl,k and Al,k are the total volume and projected area for the set of the lth HD and the kth SD, respectively. There is only one unique De l,k value for a given set of the lth HD and the kth SD (hereafter De l,k will be simply referred to as De). As shown in Fu and Liou (1993) for parameterization of the broadband radiative properties of ice clouds, and in Rolland et al. (2000) for the satellite remote sensing of cirrus clouds, the De is the key parameter representing the scattering and absorption properties of a spectrum of sizes and shapes.
Figure 1a shows the four composite ice crystal SDs (CEPEX01–04) based on an averaging procedure such that the associated ice water content (IWC) is defined in the intervals of four bins: 10−4 to 10−3, 10−3 to 10−2, 10−2 to 10−1, and 10−1 to 0.32 g m−3. Their mean effective sizes are 30, 42, 80, and 120 μm, respectively, following Eq. (4) and assuming 100% solid hexagonal columns. These composite SDs are obtained from thousands of 2D-C-observed data by averaging all spectra within IWC intervals. Since the 2D-C data are skewed toward ice crystals larger than about 90–100 μm, especially at aircraft flight speeds (∼200 m s−1), the SDs for ice crystals smaller than 90 μm were obtained from a parameterization scheme developed by McFarquhar and Heymsfield (1997). A slightly bimodal structure is noted for all the SDs. The “fine” mode, consisting of small particles, is located near the low end of the size spectrum, while the “coarse” mode is around 200–600 μm. This bimodal structure in SD is characteristic based on analyses of the CEPEX data, particularly when a significant number of ice crystals were present (McFarquhar and Heymsfield 1996). Both the numbers of large and small ice crystals increase as IWC increases. However, more large ice crystals were commonly found at higher IWCs, consistent with the fact that large ice crystals contribute more significantly to the amount of IWC (McFarquhar and Heymsfield 1997).
Figure 1b shows the other five SDs (CEPEX05–09), which were obtained at 2250, 2341, 2342, 2347, and 2348 UTC on 4 April 1993. The mean effective sizes of the five SDs are 24, 26, 38, 81, 92 μm, respectively, assuming 100% solid hexagonal columns. Data from the VIPS measurement were included to account for ice crystals with maximum dimensions less than about 120 μm. For ice particles larger than 120 μm, the SDs were obtained from 2D-C measurements. This particular date was chosen because the measurements from VIPS for the largest sizes and those from 2D-C for smallest sizes agreed better than the 1 April case. Note that VIPS data were not available for the 17 March case (McFarquhar and Heymsfield 1996). Except for the two SDs with the largest mean effective sizes of 81 and 92 μm, the coarse mode is not as obvious as that appears in the averaged SDs. The SD with the smallest mean effective size of 24 μm has the largest portion of ice crystals smaller than 40 μm but a limited number of large crystals. Although small ice crystals are less important in determining IWC, they contribute much more significantly to the ice crystal scattering properties. For example, McFarquhar and Heymsfield (1996) found that ice crystals smaller than 40 μm contribute up to 40% to cross-sectional area based on analysis of the 1 April and 4 April data, whereas their contribution to IWC is less than 20% in these two cases.
On the basis of a number of observations depicted in Table 1, we have established 11 HDs for tropical cirrus clouds in association with our study. The first six (HD1–6) were complied based on the analysis of in situ airborne observations during CEPEX (Rolland et al. 2000). In their study the existence of small ice crystals near cloud top was recognized. One measured CEPEX HD (HD7) was taken from the 2D-C data presented in McFarquhar et al. (1999) for the 4 April case. They found that the most predominant crystal habit on this particular date was polycrystals, which, in our study, are considered as aggregates. Yang et al. (2000) and Zhang et al. (2004) incorporated the droxtal shape in their HD model to represent small ice crystals observed at temperature below −35°C. Macke et al. (1996) showed that light scattering from hollow columns is representative of the scattering from most columnlike particles, including hollow columns, solid bullets, and hollow bullets. In our study, the mixture of droxtals and hollow columns is denoted as HD9, while the case involving 100% hollow columns is denoted as HD10. The case HD11 contains 100% solid columns. In a comprehensive study of ice crystal SDs for remote sensing application, Baum et al. (2005) found that the BT computed from HD11 tends to agree more closely with the observed tropical BT data, especially when small ice particles are predominant, a common feature in thin cirrus clouds.
The HD and SD models in the present study were determined based on tropical data. Because midlatitude cirrus clouds are formed by different mechanisms and under different environmental conditions, their microphysical properties could be quite different from those in the Tropics. Thus, it would appear that separate HDs and SDs should be constructed for midlatitude cirrus clouds for future remote sensing study.
c. Single-scattering properties for cirrus clouds
Following Yang et al. (2000, 2003, 2005) and Zhang et al. (2004), we have considered seven ice crystal habits including droxtals, spheroids, plates, solid and hollow columns, bullet rosettes, and aggregates. For each habit, a separate library of the microphysical properties was developed, including maximum dimension, volume, projected area, and the single-scattering properties (single-scattering albedo, extinction efficiency, and asymmetry factor) for the spectral range between 3 and 100 μm. Moreover, 45 ice crystal size bins with a maximum dimension varying from 2 to 10000 μm were used to account for SDs. For the purpose of computational efficiency, we used a fourth-order polynomial interpolation method to obtain the single-scattering properties of each ice crystal habit and size at the AIRS channel wavelengths. The bulk single-scattering properties of cirrus clouds used in radiative transfer simulations were subsequently obtained by integrating the single-scattering properties of individual ice crystals over a given set of SD and HD. Thus, the wavenumber-dependent bulk single-scattering property is determined as follows:
where ζ can be single-scattering albedo ω, asymmetry factor g, or extinction efficiency Qext; ζi(L, ν) is defined as a function of maximum dimension L and wavenumber ν for the ith ice crystal habit. Here 〈ζl,k(ν)〉 denotes the bulk single-scattering property at wavenumber ν for a combination of the lth HD and kth SD. Subscripts l and k will be omitted for simplicity. Only 〈ω(ν)〉 and 〈Qext(ν)〉 were used in the present thermal infrared radiative transfer model for thin cirrus.
3. Interpretation of clear and cloudy spectra
a. Computation of clear radiance using OPTRAN
The computationally efficient, parameterized transmittance code referred to as OPTRAN was developed by the scientists at the National Oceanic and Atmospheric Administration (NOAA; McMillin et al. 1995a) for use in satellite data assimilation in association with NWP modeling. The OPTRAN model incorporates a number of methodologies described in a series of papers since the 1970s. These advances include the rapid computation of atmospheric transmittance for constant-mixing-ratio gases at nadir view (McMillin and Fleming 1976), for different zenith angles (Fleming and McMillin 1977), and for gases with variable mixing ratios (McMillin et al. 1979). More recently, Kleespies et al. (2004) reported improvements for the OPTRAN code and an implementation that has been achieved for operational use. The most distinct approach in OPTRAN is the reversal of the usual roles of the pressure and absorber amount. The advantages of this performance include that 1) pressure is available for use as a predictor together with water vapor, ozone, and other gaseous components; 2) an arbitrary pressure profile is permitted in the calculation; and 3) the zenith angle is implicitly included in the absorber profile and explicit treatment is not necessary. In this approach, transmittance is computed through a recurrence procedure in which only simple calculations are involved.
The OPTRAN method has been applied to a number of instruments by Kleespies et al. (2004) who have tested the OPTRAN’s accuracy by comparing with the calculation from a line-by-line radiative transfer model (Clough et al. 1992). Differences between the two in terms of brightness temperature were less than 0.2 K over the infrared spectrum. OPTRAN has also been utilized for the calculation of transmittances of the clear atmosphere for ∼600 AIRS channels that are dominated by water vapor absorption (Strow et al. 2003).
b. Computation of thin cirrus radiance
The atmospheric profiles used for radiance calculations are adopted from the AIRS level 2 Support Retrieval product datasets, including temperature, water vapor, and ozone, which are reported on 100 pressure levels from 1100 to 0.016 hPa. The surface skin temperature, surface pressure, and the index of the last physically meaningful profile level above the mean surface are utilized to determine the surface emission and the number of actual pressure levels above the surface. It is necessary to check this index in the calculation and extrapolate in the log (pressure) domain to obtain the surface value of water vapor and ozone (Olsen 2003). The satellite viewing angle is taken from the AIRS viewing geometry given in the AIRS full swath data fields. The spectral IR surface emissivity and bidirectional surface reflectivity are adopted from the AIRS level 2 standard retrieval product and given in order of increasing frequency. Manus Island cases were analyzed using version 3 (Olsen 2003) of the AIRS retrievals, while Nauru Island cases used version 4 (Olsen 2005). The selection of these cases will be discussed in section 4b. Employing all the aforementioned AIRS variables, clear radiances can be computed from OPTRAN.
The bulk single-scattering properties for clouds were calculated according to Eq. (5), and the cloud-top temperature was extracted from the AIRS level 2 standard product in which it was given within a pixel. Using the relationship between τIR and τ [Eq. (3)], the cloud emittance ɛν is computed based on Eq. (2). The cloudy radiance is then evaluated following Eq. (1), where the clear radiance is obtained from OPTRAN. The simulated BT spectrum, which is converted from the cloudy radiance spectrum, is a function of τ, De, and HD models. Using the aforementioned radiative transfer parameterization, we have constructed look-up tables, which contain the simulated BT spectra in terms of these three parameters. The look-up tables are applied to a minimization method to search for the optimal combination of cirrus parameters that are associated with the minimum residual.
Many efforts have been devoted to the validation of AIRS products (Fetzer 2006). By comparing the observations of AIRS and the airborne Scanning High-Resolution Interferometer throughout the infrared spectrum, Tobin et al. (2006) showed that the absolute values of mean BT differences were found to be less than 0.2 K for most channels. The AIRS sounding accuracy, in terms of RMS differences between AIRS retrieval and radiosonde, is: 1 K per 1 km layers for temperature and 15% per 2 km layers for water vapor (Divakarla et al. 2006). Kahn et al. (2006, manuscript submitted to J. Geophys. Res.) found that AIRS-measured cloud-top heights statistically agree with MMCR and Micropulse lidar (MPL) measurements at ARM TWP sites and for cases associated with thin cirrus. Uncertainty in the cloud-top height is from −1.1 to 2.1 km on average.
c. Determination of cirrus optical depth and particle size using a χ2-minimization method
For each observed spectrum, a total of 9900 modeled spectra are calculated. These correspond to combinations of 9 CEPEX size distributions, 11 tropical cirrus habit distributions based on tropical measurements, and 100 visible cirrus optical depths varying from 0.01 to 1.0 inputted to the cloudy radiative transfer model. Within the 800–1130 cm−1 thermal window, we select 14 “clean” subbands for their minimal water vapor absorption effect. They cover 220 AIRS channels and their starting and ending wavenumbers are given in Table 2. These subbands have also been used in the Atmospheric Emitted Radiance Interferometer (AERI) analysis (DeSlover et al. 1999). Not all AERI subbands were used in our simulation. Only those covering more than 10 AIRS channels were adopted for obtaining statistically meaningful average BTs. Then a mean BT is calculated for each subband and is denoted as BTj(Model), where the index j ranges from 1 to 14, corresponding to the 14 subbands selected. As a result, a total of 9900 subband BTj(Model) spectra are computed. For each observed AIRS BT spectrum, we also obtain one set of 14 mean BTs, BTj(AIRS). Then we define a parameter χ2 as the sum of squared differences between BTj(Model) and BTj(AIRS):
where χ2 varies with different SD, HD, and τ. Note that the χ2 defined above is not related to the standard statistical χ2 test. We search for the minimum χ2 to determine De, HD, and τ corresponding to the best-fit spectrum.
The search method may fail if the curve of χ2 with respect to one of the three parameters is flat, indicating that there is no constraint on this parameter. For this reason, we carefully examined the χ2 functions in reference to the three parameters. Figure 2 shows an example of the χ2 patterns with respect to SD, HD, and τ. We found that none of the curves displays flatness, illustrating that the minimization method has constraint on all three parameters. In particular, χ2−SD and χ2−τ shown in Figs. 2a,c reveal a more tightly constraint of χ2 on SD and τ. The χ2−HD curve in Fig. 2b does not show a systematic pattern because of the arbitrary order of HDs in Table 1.
4. Results and discussions
a. Sensitivity study
We used the standard tropical atmospheric profile (McClatchey et al. 1972) to investigate the sensitivity of thermal IR window BT spectra to optical depth and mean effective size that are determined by both SD and HD functions. Surface and cloud-top temperature are prescribed to be 300 and 200 K, respectively, while surface infrared emissivity and reflectivity are set to be 1.0 and 0.0, respectively, to represent the thermal infrared radiation features over the tropical ocean. Two experiments have been carried out. The same habit distribution (HD11) and different size distributions (De = 42 μm for CEPEX02 and De = 80 μm for CEPEX03) were used in the first experiment, while in the second experiment, the size distribution was fixed (CEPEX05) and different habit distributions (De = 15 μm for HD7 and De = 28 μm for HD9) were used.
Differences of the BT spectra are due to the different radiative properties of ice clouds associated with various combinations of HDs and SDs. Figure 3 shows the bulk extinction efficiency and single-scattering albedo for the four HD–SD combinations that were used in the two experiments within the spectral range 750–1150 cm−1. The bulk extinction efficiency and single-scattering albedo strongly depend on wavenumber in both experiments; however they are less dependent on the mean effective ice crystal size in the first experiment than those in the second.
Figure 4 shows the subband BT spectra simulated in the first experiment for ice clouds with optical depths of 0.1 (left panel), 0.5 (middle panel), and 1.0 (right panel). Clear BT spectra are also shown in each panel by the solid curves with squares for comparison with thin cirrus cloudy spectra. The spectra are displayed in terms of mean BTs for the 14 subbands introduced in section 3c. Between 800 and 960 cm−1, the slope of the mean subband BTs, defined as ΔBT/Δν, increases with decreasing mean effective size. Between 1050 and 1130 cm−1, BTs decrease by more than 10 K as optical depth increases from 0.1 to 1.0. Compared to the BT spectra between 1050 and 1130 cm−1, the spectra between 800 and 960 cm−1 show a stronger sensitivity to wavenumber. Also, the sensitivity of the BT slope to the mean effective size is larger. These results are in agreement with the findings reported by Huang et al. (2004) and Wei et al. (2004). The left panel shows that even for the cirrus optical depth as small as 0.1, the BT difference between cirrus and clear can be as large as 3 K. This sensitivity forms the basis for the retrieval of thin cirrus optical depth using the thermal infrared spectra.
In the second experiment, we used the same SD (CEPEX05) but varying HD. The results are shown in Fig. 5. Although different mean effective sizes (De = 15 μm for HD7 and De = 28 μm for HD9) are produced, BT slopes over 800–960 cm−1 show a small variation despite the fact that different mean effective sizes are used. We also find that colder BTs are associated with larger mean effective sizes. The BT differences between the two spectra using different HDs can be 1, 5, and more than 10 K corresponding to the cirrus optical depths of 0.1, 0.5, and 1.0, respectively. It appears evident that both HD and SD are important in determining the ice cloud bulk radiative properties and in defining the shape of the BT spectra.
b. Interpretation of the AIRS clear and cloudy spectra over the ARM TWP site
The AIRS is a grating spectrometer on board the NASA’s Aqua satellite, launched on 4 May 2002. It covers a spectral region from 3.7 to 15.4 μm and has 2378 channels with a nominal spectral resolution of λ/Δλ = 1200. The nadir spatial resolution is 13.5 km from the nominal 705-km orbital height. It views the ground through a cross-track rotary scan mirror that provides a ±49.5° ground coverage every 2.67-s scan cycle. A total of 90 ground footprints containing all 2378 spectral samples are observed in each cycle. According to the prelaunch testing of AIRS, the absolute radiometric accuracy is 0.2 K (Pagano et al. 2003).
We applied the methodology developed in this paper to three nighttime AIRS granules; one is clear and the other two contain thin cirrus. An AIRS data granule contains the data gathered in 6 min for an area of about 2400 km by 2400 km. We selected granule 158 (28 May 2003, clear), granule 159 (20 June 2003, thin cirrus), and granule 147 (1 March 2005, thin cirrus). All the granules are over the ARM TWP sites and adjacent oceans, among which the first two are over the site at the Manus Island (2.006°S, 147.425°E) and the last one is over the Nauru Island (0.521°S, 166.916°E), ideal for comparison with independent ARM measurements. In the cloud-screening method developed by Aumann et al. (2004), a spectral cloud filter was combined with a spatial coherence filter to yield cloud-free spectra with an accuracy of 0.2 K, the same as the AIRS prelaunch radiometric accuracy. Only 1%–2% of the AIRS data passed this test. The clear granule was selected based on this approach, which indicated that clear pixels appeared in this granule. Further, for identification of thin cirrus we follow the detection scheme developed by Kahn et al. (2005), in which the 3.8- and 10.4-μm BT differences were used along with the total column precipitable water measurements. The collocated and coincident browse images obtained by the Moderate Resolution Imaging Spectroradiometer (MODIS) were used to provide a supplementary reference to the cloudy conditions in the selected granule.
Because of the large variation of the cloud properties in an AIRS granule, we selected a domain inside the granule based on the coincident MODIS browse images. For the 28 May 2003 case, the domain covers an area over 1.5°S–2.5°S and 147°E–148°E. Within this area, five AIRS pixels are considered “clear” according to the detection schemes described above. Figure 6 shows the average residuals between the AIRS-observed and model-simulated BT spectra for the five clear pixels for the individual AIRS channels as well as averaged over the subbands. The zero points indicate the AIRS observations were missing. For the 9.6-μm ozone band, the comparison is less satisfactory. A similar pattern is displayed in the residual analysis of the five AIRS pixels (bottom panel). The averaged (AIRS-observed BTs minus model calculation) differences (averaged over five pixels) indicate that, except for the 9.6-μm ozone band, differences between computed results and observed data are generally less than 1 K, with the AIRS BTs being colder. The colder bias in the AIRS observation may indicate the potential contamination by thin cirrus clouds in the selected scene. For channels in the subbands, the BT differences between the AIRS observation and model results are within −0.5 K, except for the band at 850 cm−1, which has a value of about −1.0 K. The reason for discrepancy in the ozone band could be related to uncertainty of the ozone profile in the AIRS data product.
To simulate BTs corresponding to the thin cirrus AIRS granule on 20 June 2003, we have selected an area covering the region 1.0°S–3.0°S and 147.0°E–149.0°E that contains 15 pixels; the results are shown in Fig. 7. In the first pixel (upper left), the cirrus has a visible optical depth of 0.16 and a mean effective size of 81 μm, a result of a combination of the CEPEX08 SD and the HD with 100% hexagonal columns (HD11). Curves in upper right demonstrate another pixel for a cloud with a visible optical depth of 0.20 and a mean effective size of 78 μm having a different HD (HD9). The average residuals (AIRS minus model) for the 15 thin cirrus cloudy pixels are shown in the lower panel. Similar to the clear case, residuals are within 0.5 K in the subbands between 800 and 960 cm−1 region, while the largest residuals occur in the ozone band. Differences in the 1050–1130 cm−1 region are larger but still within −1 K for the subband channels. The AIRS cirrus BTs are warmer than the model results. It is possible that the neglect of multiple scattering may contribute to the bias in this model calculation.
In granule 159 (20 June 2003) the selected 15 pixels were all thin cirrus cases with the AIRS-inferred optical depths less than 0.3, while granule 147 (1 March 2005) contains optically thicker cirrus clouds. The domain selected around the Nauru Island site covers the region 0°–1.5°S and 165°–167.5°E and contains 14 AIRS pixels, 7 of which have the AIRS-inferred optical depths smaller than 0.3 and their simulation results (not shown) are similar to those for granule 159 in Fig. 7. However, there are three thicker cirrus pixels with residuals larger than 1 K in the subbands that have optical depths of 0.69, 0.4, and 0.32, respectively. This suggests that multiple-scattering contributions become important as cirrus optical depth becomes larger than 0.3, which is the upper limit for application of the current radiative transfer model.
Based on the AIRS validation studies discussed in section 3b, we have carried out different numerical experiments for the case on 20 June 2003 to study the effects of input parameter uncertainty and instrument noise on the retrieval accuracy. The AIRS BT spectra, temperature profile, water vapor profile, cloud-top temperature, ozone profiles, and surface temperature have been perturbed by ±0.2 K, ±1.0 K, ±15%, ±6.5 K, ±20%, and ±1.0 K, respectively. These are bias values and not standard deviations. Results show that uncertainties in the water vapor profile are the major error source, which causes an error in optical depth less than 0.05. Errors in the effective particle diameter on average are less than 8 μm.
c. Validation using ARM TWP data
To have an independent check on the model results, we compare the optimal mean effective ice crystal size and optical depth to those determined from the MMCR retrieval developed by Mace et al. (2002). The extinction coefficient can be calculated by using the following equations (Fu and Liou 1993):
where IWC and De denote the IWC and mean effective size retrieved from the MMCR data, respectively; βext is the IR extinction coefficient; Δz is the cloud thickness obtained from the MMCR retrieval; and c and d are constants. For columns, c ≈ −6.656 × 10−3 and d ≈ 3.686 (Fu and Liou 1993). Based on these parameterization equations, optical depth can be computed as a function of time and altitude.
To validate the calculations for the two cirrus granules, we used the mean wind of the cloud layer on the basis of the nearest wind-sounding profile in the collocation method, and converted the MMCR retrieval results from time domain to spatial domain to compare with satellite measurements. This assumes that the wind speed and direction remained the same during the time period between the AIRS overpass of ARM sites and the collocation times of different pixels. Among the 29 thin cirrus cloudy pixels reported in this paper, only those along the wind direction could be collocated with the ground-based measurements. Cirrus parameters were assumed to remain unchanged during the collocation time period. If a pixel is too far away from ARM sites, then the collocation time of this pixel would be far apart from the AIRS overpass, and collocation and comparison of this AIRS pixel and the MMCR time series would not be considered “true collocation” for validation. Since AIRS have a nadir spatial resolution of 13.5 km (Pagano et al. 2003), and the wind speed at the cirrus layers was very small (5 m s−1 on 20 July 2003 and 10 m s−1 on 1 March 2005), it took about 1 hr for clouds in one pixel to drift to the adjacent pixel. A 1-hr average was taken for the ground-based measurements.
Given the criteria of the collocation method, only four pixels were considered as “true collocation” for this comparison in the granule 159 from 20 June 2003 (pixels 14, 10, 8, and 4). As shown in the Table 3, the AIRS-inferred optical depth is slightly larger than the MMCR retrieval value, while the AIRS-inferred mean effective size is smaller than the MMCR retrieval result. This discrepancy could be explained by the fact that MMCR generally misses small particles (Comstock et al. 2002). We note that for pixel 14, MMCR completely missed the thin cirrus layer, with τ and De approximately equal to zero. However, the AIRS-inferred result shows that a thin cloud was present having a De of 19.6 μm and a τ of 0.01. This may be due to MMCR’s missing small ice crystals, but may also be due in part to our analysis and the possible spatial and temporal variations in cirrus.
Applying the aforementioned criteria of collocation to the 11 pixels (granule 147, 1 March 2005) with the condition that the simulated BT spectra closely agree with AIRS-measured spectra, four pixels have been selected for comparison (6, 9, 10, and 13), as shown by Table 4. The AIRS-inferred and MMCR-retrieved optical depths for these four pixels are from 0.05 to 0.37 and from 0 to 0.18, respectively. The effective particle sizes range from 20 to 91.5 μm and 0 to 72 μm based on the AIRS and MMCR retrievals, respectively. Pixel 6 collocated with a thin cirrus layer between 17 and 18 km that occurred 4 hr earlier than the AIRS overpass at the Nauru Island on that day. This thin cirrus layer was missed by the MMCR detection and only captured by the simultaneous MPL measurements. For this reason, the results of MMCR retrievals in this pixel are ∼0 for both optical depth and particle size. As expected, for pixels 9 and 13, the AIRS-inferred optical depths (mean effective sizes) are larger (smaller) than the MMCR retrieval values. However, exception was found for the mean effective size of pixel 10, in which the AIRS-inferred optical depth is 0.37. Although good agreement between the calculated and AIRS-measured BTs has also been obtained for this pixel, the cirrus clouds with optical depths larger than the 0.3 threshold may be in error because of the neglect of multiple scattering in retrieval analysis.
Although the normalized MPL backscatter profiles with all known instrument artifacts removed was not made available for this case study, we have made a preliminary comparison with the MPL retrieval based on raw data. Because pixels 10 and 13 contain thicker cirrus, MPL did not have backscattering signals at the collocated times. Therefore, comparisons with MPL were made only for pixels 6 and 9 where the MPL-retrieved optical depths are 0.02 and 0.15, respectively (J. M. Comstock 2005, personal communication).
In this study, we have developed a fast radiative transfer model for the interpretation of AIRS data by combining the OPTRAN model with a thin cirrus parameterization. This effort includes the construction of size and habit distribution models for cirrus clouds and the calculation of their bulk optical properties. As a demonstration, we analyze 5 clear pixels and 29 thin cirrus cloudy pixels in three AIRS data granules coincident with the ARM program sites at Manus Island and Nauru Island in the tropical western Pacific. In the clear cases, differences between the AIRS-observed spectral BTs and the modeled results are within −0.5 to −1 K for the subband average at 800–1130 cm−1. In the thin cirrus case, we find that differences in the subbands channels are within +0.5 to −1 K. Relatively larger differences on the order of 4–5 K occur in some absorption lines, including those in the ozone 9.6-μm band region for both clear and cloudy cases. These differences may be related to uncertainties in the atmospheric state, which is taken from the AIRS level 2 operational support product.
We have also developed a χ2-minimization method to infer the optical depth and mean effective ice crystal size in thin cirrus conditions. We show that the AIRS-inferred optical depth and mean effective size are consistent with the results independently determined from the collocated and coincident MMCR measurements available at the ARM TWP Manus Island and Nauru Island sites for the cases selected.
Finally, we have investigated the sensitivity of the spectral BT to cirrus optical depth and ice crystal size and habit. Cirrus cloudy radiances in the 800–1130-cm−1 window are shown to be sensitive to all three parameters. The degree of their sensitivity varies spectrally, suggesting that it appears feasible to use BT spectra to infer them simultaneously. Sufficient sensitivity of the ice crystal habit in the thermal infrared spectrum illustrated in our analysis appears to be a new finding. Additional data analysis is required to determine the different roles of ice crystal size and habit in shaping the AIRS spectra. In the present theoretical modeling for fast radiance calculations in thin cirrus conditions, the effect of multiple scattering of thermal emission has not been accounted for. For cirrus clouds with optical depths larger than about 0.3, multiple scattering may become significant, an area requiring further investigation. The errors in the computation caused by uncertainties in AIRS-measured BT spectra, atmospheric profiles, and surface and cloud temperature have been studied and they are relatively small compared with the sensitivity of thermal infrared spectrum to optical depth. Last, updated SD and HD models that can realistically represent the tropical and midlatitude cirrus cloud conditions should be constructed in future remote sensing studies.
We thank Annmarie Eldering for suggesting further examination of applying the χ2-minimization method to AIRS radiances, and Jennifer M. Comstock for providing MPL raw data retrieval. This research has been supported by DOE Grant DE-FG03-00ER62904 and NOAA Grant NA03NES4400011.
Corresponding author address: Ms. Qing Yue, Department of Atmospheric and Oceanic Sciences, University of California, Los Angeles, 405 Hilgard Ave., Los Angeles, CA 90095–1565. Email: email@example.com
This article included in the Assimilation of Satellite Cloud and Precipitation Observations special collection.