JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117, D19202, doi:10.1029/2012JD017799, 2012
Asian dust height and infrared optical depth retrievals over land from hyperspectral longwave infrared radiances Zhigang Yao,1 Jun Li,1 Hyo-Jin Han,2 Allen Huang,1 B. J. Sohn,2 and Peng Zhang3 Received 17 March 2012; revised 14 August 2012; accepted 18 August 2012; published 4 October 2012.
[1] The dust top height and infrared optical depth over land are retrieved from the Atmospheric Infrared Sounder (AIRS) longwave infrared measurements by using a one-dimensional-variation retrieval algorithm for different Asian dust storms. By combining particle size measurements from a 10-year ground observation data set from the Dunhuang Skynet station located to the east of the Taklimakan Desert in China and the Optical Properties of Aerosols and Clouds data set of optical properties, the mineral dust scattering and absorption coefficients are obtained and then used to compute brightness temperatures with RTTOV 9.3. The retrieved dust thermal infrared optical depths are compared with the Ozone Monitoring Instrument and Moderate Resolution Imaging Spectroradiometer (MODIS) products. The retrieved dust top heights are compared against the extinction backscatter profiles obtained from the Cloud-Aerosol Lidar with Orthogonal Polarization lidar. Infrared optical depths from AIRS correlate favorably with visible optical depths from MODIS, and dust top heights agree reasonably with lidar observations for the single-layered dust storms over the Taklimakan Desert. Citation: Yao, Z., J. Li, H.-J. Han, A. Huang, B. J. Sohn, and P. Zhang (2012), Asian dust height and infrared optical depth retrievals over land from hyperspectral longwave infrared radiances, J. Geophys. Res., 117, D19202, doi:10.1029/2012JD017799.
1. Introduction [2] Dust storms occur more frequently now due to climatic variability and land use change [Jickells et al., 2005; Lee and Sohn, 2011]. Asian dust generated mostly in the springtime in northwest China, one of the major sources of large-scale dust storms, is transported frequently toward the Pacific Ocean. Each year, thousands of tons of dust from the Taklimakan Desert are blown into the Pacific Ocean. In recent years, researchers have found that Asian dust storms could affect the radiation budget through the regulation of cirrus cloud formation and the marine ecosystem through nutrient supply [Eguchi et al., 2009; Uno et al., 2009]. Furthermore, the springtime aerosol load in the United States Pacific Northwest is related to Asian dust emissions [Fischer et al., 2009]. [3] Accurate dust storm information is critically important for climate change studies and numerical weather prediction 1 Cooperative Institute for Meteorological Satellite Studies, University of Wisconsin-Madison, Madison, Wisconsin, USA. 2 School of Earth and Environmental Sciences, Seoul National University, Seoul, South Korea. 3 National Satellite Meteorological Center, China Meteorological Administration, Beijing, China.
Corresponding author: Z. Yao, Cooperative Institute for Meteorological Satellite Studies, University of Wisconsin-Madison, 1225 West Dayton St., Madison, WI 53706, USA. (
[email protected]) ©2012. American Geophysical Union. All Rights Reserved. 0148-0227/12/2012JD017799
[Hansen et al., 1997; Lubin et al., 2002; Seinfeld et al., 2004; Sun et al., 2001; Intergovernmental Panel on Climate Change, 2007]. The study conducted by De Souza-Machado et al. [2006] indicates that the dust contamination has significantly negative impacts on atmospheric parameter retrievals and degrades the use of the data in weather forecasting. Furthermore, the retrieval of surface properties from space is affected by dust contamination [Díaz et al., 2000]. Additionally, current understanding of the radiative forcing of dust aerosols is limited, especially over dust source regions where ground observations are few [Kim et al., 2005] and satellite retrievals at visible to nearinfrared wavelengths are often difficult to carry out because of the high surface albedo [Kaufman et al., 2002]. [4] High spectral resolution remote sensing provides the capability to detect dust storms and retrieve their heights and optical depths [Sokolik, 2002; Pierangelo et al., 2004; De Souza-Machado et al., 2006]. Sokolik [2002] reported that the radiative signatures of dust are significantly different from those of clouds and the greenhouse gases in the atmosphere. De Souza-Machado et al. [2006] investigated infrared spectral signatures of dust over ocean using the Atmospheric Infrared Sounder (AIRS) infrared measurements. Pierangelo et al. [2004] developed an approach based on precomputed lookup tables to retrieve dust altitude and infrared optical depth over the Atlantic Ocean. The dust optical depth and the dust height in the Mediterranean are also retrieved by using the Thermal Infrared (TIR) observations from the AIRS [De Souza-Machado et al., 2010].
D19202
1 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
[5] Although great efforts have been devoted to extracting dust properties from high spectral resolution infrared measurements, certainly more work is needed on the dust retrieval from these data, especially over land [Peyridieu et al., 2010; De Souza-Machado et al., 2010]. Dust parameter retrievals over land are difficult mainly due to poor knowledge of the surface emissivity and inaccurate surface temperature as well as less accurate temperature of the underlying atmosphere. The retrievals may be complicated further by the fact that the single scattering properties of dust originating from different areas differ primarily due to different dust composition [Su and Toon, 2010]. [6] This study focuses on retrieving the altitude and thermal optical depth of Asian dust over land by using the TIR hyperspectral sounding data from AIRS. The mineral dust optical properties derived from a 10-year ground size distribution observation data set from the Skynet station [Kim et al., 2004] in Dunhuang (40.146 N 90.799 E), which is located to the east of the Taklimakan, are used in the radiative transfer model [Han et al., 2012]. A one-dimensional variation retrieval (1DVAR) algorithm is developed to extract dust parameters by using RTTOV 9.3 (Radiative Transfer for (A)TOVS) forward model. The retrieved dust TIR optical depths are compared with the Ozone Monitoring Instrument (OMI) aerosol index [Torres et al., 2007] and the Moderate Resolution Imaging Spectroradiometer (MODIS) Deep Blue product [Hsu et al., 2004]. The retrieved dust top heights are compared against the extinction backscatter profiles obtained from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) lidar [Winker et al., 2007]. [7] This paper is organized as follows: section 2 introduces the data used in this study, section 3 outlines the dust height and infrared optical depth retrieval method, section 4 investigates the accuracy of the dust property retrieval in different assumed conditions, and section 5 presents the detailed analysis of dust top heights and optical depths retrieved from the AIRS measurements by comparison with other A-Train instruments for four different Asian dust storms. Finally, conclusions are presented in section 6.
2. Data [8] The A-Train constellation of satellites provides multi sensor observations to study cloud and aerosol properties, as well as weather and atmospheric chemistry [Stephens et al., 2002]. The measurements are made over the same geographic locations and separated by less than 15 min. The instruments include MODIS and AIRS on Aqua, OMI on AURA, and CALIOP on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO). Since the satellites in this constellation fly closely to each other, the dust retrievals from AIRS can be easily collocated and compared to the products from the OMI, MODIS, and CALIOP in the study. 2.1. AIRS [9] The AIRS instrument [Chahine et al., 2006], onboard the Earth Observing System (EOS) Aqua satellite, was launched on 4 May 2002 and is the first of a new generation of high spectral resolution IR sounders. It has 2378 channels measuring outgoing radiances in three bands: 8.8–15.4, 6.2–
D19202
8.2, 3.75–4.58 mm at a spatial resolution of 13.5 km at nadir. AIRS was primary designed to retrieve atmospheric temperature, water vapor, trace gases, and surface skin temperature with a 2000 km swath for weather prediction and climate studies [Chahine et al., 2006]. Its measurements also contain information about dust properties [Ackerman et al., 2004; Pierangelo et al., 2004; De Souza-Machado et al., 2006]. Additionally, it can be used to retrieve cloud microphysical properties by taking full advantage of high spectral resolution AIRS longwave cloud-sensitive radiance measurements [Li et al., 2005; Weisz et al., 2007]. 2.2. OMI [10] The OMI on the EOS Aura platform is a hyperspectral push-broom sensor measuring solar reflected and backscattered light in two spectral bands: the UV channel (270– 350 nm) and the VIS channel (350–500 nm) [Levelt et al., 2006]. This instrument has a 2600 km wide viewing swath such that it provides daily global coverage at a spatial resolution 13 24 km at nadir. [11] Arid and semiarid regions are problematic in retrievals using visible and near-IR measurements, as these surfaces have high reflectance in these spectral bands. However, low surface albedo in the UV spectral region for all land and water surfaces means that OMI can be used to retrieve aerosol properties [Torres et al., 2007]. There are two aerosol inversion algorithms applied to the OMI measurements: the OMI near-UV (OMAERUV) algorithm and the multiwavelength algorithm (OMAERO). The two algorithms make use of the sensitivity of UV measurements to aerosol absorption [Torres et al., 1998, 2007]. The OMAERUV is used to derive aerosol extinction and absorption optical depth with two UV wavelengths. This algorithm produces the Aerosol Index, and aerosol optical depth at 388 nm. OMI aerosol retrievals have been validated against AERONET observations [Torres et al., 2007; Curier et al., 2008; Brinksma et al., 2008]. The UV-absorbing aerosols, such as desert dust carbonaceous and volcanic ash, yield positive aerosol index values. The OMI aerosol index is also sensitive to the aerosol height. [12] In this study, the OMI aerosol index-UV is used to evaluate the dust AOD (aerosol optical depth) retrievals from AIRS over land. The OMI aerosol index is calculated as a measure of the difference between the measured backscattered UV radiance from an atmosphere containing aerosols and that of a pure Rayleigh scattering atmosphere [Torres et al., 2007]. 2.3. MODIS [13] The MODIS instrument on board the Aqua satellite in the A-Train has high spatial resolution and a near-daily global coverage. It provides observations of the Earth’s surface and atmosphere using 36 channels centered in the visible (VIS), near-infrared (NIR), and infrared (IR) regions of the spectrum from 0.4 to 14.5 mm. The MODIS retrieval of AOD over land employs primarily three spectral channels (0.47, 0.66, and 2.1 mm) [Kaufman et al., 1997a]. However, the standard MODIS AOD product cannot provide information about aerosols over bright surfaces such as deserts because it is retrieved using NIR wavelengths (2.1 and 3.8 mm) [Kaufman et al., 1997b]. In this paper, the MODIS Deep Blue (DB) product (10 km resolution) is used to evaluate the dust AOD from AIRS. The DB algorithm has
2 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
larger sensitivity to aerosols over bright surfaces because it primarily employs two blue channels (0.412 and 0.470 mm), for which surface reflectances are much darker, to infer aerosol properties [Hsu et al., 2004]. It is reported that the uncertainties of the MODIS DB product are around 20–30% by comparison with observations from AERONET sites in China and Mongolia [Hsu et al., 2006]. Li et al. [2012] use the ground-based remote sensing of AOD from Sun photometers at four sites in Xinjiang during the years 2002–2003 to validate aerosol products, including DB of MODIS. The results show that the DB algorithm is able to retrieve AOD at Tazhong in the Taklimakan desert, where the correlation coefficient between the DB product and ground-based Sun photometer measurements reaches 0.94. 2.4. CALIOP [14] Operating at optical wavelengths, CALIOP is a spacebased lidar on board the CALIPSO satellite platform, and provides unique measurements of the global vertical distributions of clouds and aerosols, including dust [Winker et al., 2007]. CALIOP transmits and receives backscattered light at laser wavelengths of 532 nm and 1064 nm. It has a typical horizontal averaging interval of 5 km for aerosols. The observations from CALIOP about 45 s behind AIRS will be used to validate the AIRS dust height retrievals. More information on CALIPSO instruments and products is available at http://www-calipso.larc.nasa.gov/.
3. Physical Basis and Method 3.1. Physical Basis [15] The TIR observations could be used to detect dust over desert both day and night [Wald et al., 1998; Sokolik, 2002]. The impacts of the aerosol optical depth, aerosol layer altitude, surface emissivity, and size distribution on the AIRS observations were also investigated by Pierangelo et al. [2004]. Their study demonstrated that the longwave IR channels are very sensitive to the altitude and aerosol optical depth. As a result, it is possible to simultaneously retrieve the AOD and altitude of a dust layer using AIRS measurements simultaneously. 3.2. Forward Model [16] In the physical retrieval, a forward model is used to model satellite observations as accurately as possible. The forward model F(x) relates the atmospheric state vector x to the brightness temperature measurement vector Y. F, therefore, has the same dimension as Y. In this study, the dust optical properties used are calculated with the Optical Properties of Aerosols and Clouds (OPAC) mineral refractive index (http://www.lrz.de/uh234an/www/radaer/opac. html). The particle size distributions are obtained by taking an average of all available dust observations at Dunhuang. Then we combine three modes for mineral dust components (nucleus, accumulation, and coarse modes) from OPAC database, and the number mixing ratios for the three modes are 0.862, 0.136, and 0.002, respectively. Considering the number size distribution, the mixing ratio of nucleus mode is higher than the other modes. However, for the volumetric size distribution considered, the coarse mode is dominant. It is a typical pattern of dust size distributions. For example,
D19202
desert aerosol type of OPAC has a very low number mixing ratio of the coarse mode. [17] Given the atmospheric parameter profiles, surface temperature, surface emissivity, and dust parameters, the AIRS observation vector has the form Y ¼ FðxÞ;
ð1Þ
where the individual elements are calculated by RTTOV 9.3, which is a development of the fast radiative transfer model originally developed at European Centre for Medium-Range Weather Forecasts (ECMWF) in the early 1990s [Eyre, 1991] for TOVS. The model allows rapid simulations of radiances for satellite infrared or microwave nadir scanning radiometers given atmospheric parameter profiles. This model, which includes a multiple scattering code, can be used to simulate radiances at infrared wavelengths in the presence of aerosols. Han et al. [2012] has used this model to compare the simulated measurements to the AIRS observations in cases of dusty skies. 3.3. Surface Emissivity [18] Dust retrieval over continental surfaces, such as the Taklimakan Desert, is possible because of the availability of surface emissivity data sets from MODIS, AIRS, and Infrared Atmospheric Sounding Interferometer (IASI) observations. These monthly mean data sets are derived from the clear sky observations, although they may be somewhat contaminated by dust and cloudy pixels which are misidentified as clear. This information paves the way for retrieving dust properties over land with infrared measurements. [19] Seemann et al. [2008] derived a high spatial and moderate spectral resolution global database of Land Surface Emissivity (LSE) by using a procedure called the baseline fit (BLF) method. On the basis of an analysis of MODIS/UCSB and ASTER/JPL libraries, they extracted the most representative spectra (about 300) and identified 10 wavelengths (“hinge points”) situated within strategic regions of the emissivity spectrum. According to a conceptual model of LSE, the method adjusts a baseline emissivity spectrum based on the operational MODIS product MYD11 from Aqua or MOD11 from Terra. This approach captures the most important features of the emissivity spectrum. In their study, the impact of LSEs on atmospheric total precipitable water (TPW) retrievals with MODIS radiance measurements was evaluated over land. Substantial improvement was shown over traditional retrievals conducted with the typical assumption of constant LSEs. The surprisingly better results by using the MODIS products in the 11–12 mm range are reported by Li et al. [2010], although the method could not provide the LSE quantitative accuracy due to the lack of the true emissivity. The original data set provides monthly averaged global IR LSE spectra from 3.7 to 14.3 mm at a 0.05 0.05 spatial resolution based on the analysis of MODIS data. To make this data set closer to AIRS observations in their spatial resolution, the data are degraded to a 0.25 0.25 degree latitude-longitude grid scale in the following study. 3.4. Channel Selection [20] It is well known that the physical retrievals strongly depend on the accuracy of the forward model for selected
3 of 15
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
D19202
Table 1. AIRS Selected Channels Number, Wave Number, Wavelength, Surface Transmittance Using the Tropical Atmospheric Parameters in RTTOV 9.3 Channel Number
Wave Number (cm1)
Wavelength (mm)
Surface Transmittance
671 672 673 676 914 924 929 933
870.792 871.114 871.436 872.404 965.215 969.222 971.237 972.856
11.484 11.450 11.475 11.463 10.360 10.318 10.296 10.279
0.29 0.17 0.33 0.54 0.65 0.63 0.47 0.55
channels. Significant discrepancies in the brightness temperatures could be caused by the differences in dust composition, which could strongly depend on geographic source regions [Claquin et al., 1999; Sokolik, 2002; Su and Toon, 2010]. As a result, the channels should be carefully selected after detailed radiative transfer simulations using the available dust composition in a particular area. [21] Han et al. [2012] calculated the aerosol scattering and absorption coefficients by using the mineral aerosol size distribution data measured at the Dunhuang dust observation station located at 94.8 N and 40.1 E from October 1998 to January 2007 and the OPAC data set of optical properties. Then, these coefficients and the CALIPSO dust observations in 2008 are used to simulate the AIRS observations with RTTOV 9.3. Their results show that the discrepancies between the simulated and observed brightness temperatures are less than 0.5 K in the 10.2–11.5 mm band when the dust effect is taken into account. Consequently, only eight thermal channels in this band as shown in Table 1 are used to retrieve the dust top height and infrared optical depth. The sensitivity of this band to the dust has been extensively investigated by other studies [Pierangelo et al., 2004; De Souza-Machado et al., 2006; Peyridieu et al., 2010]. Channels 672 and 914 have been used by Peyridieu et al. [2010]. Although the shortwave infrared channels could provide more information on the dust optical depths, they could be contaminated by the effects of surface reflectance of downwelling solar radiance, degrading the retrievals during the daytime. 3.5. Retrieval Algorithm [22] The general idea of the 1DVAR theory is to minimize a penalty function J(x), which measures the degree of fit of the dust measurements to the background information and possibly to other physical constraints. The 1DVAR or Optimal Estimation Method (OEM) tries to minimize the following penalty function [Rodgers, 1976]: x xa Þ þ ½y Fð^xÞT S1 J ð^xÞ ¼ ð^x xa ÞT S1 xÞ; ð2Þ a ð^ ɛ ½y Fð^
where xa is the background information vector of x, Sa is the assumed background error covariance matrix that constrains the solution, and Sɛ is the observation error covariance matrix, which includes instrument noise plus the assumed forward model error. To numerically solve the equation, a quasi-Newtonian iteration is used to obtain [Eyre, 1989]: 1 d^xnþ1 ¼ KTn S1 ɛ ⋅ K n þ Sa
1
⋅ KTn S1 xn Þ; ð3Þ ɛ ⋅ ðdYn þ K n ⋅ d^
D19202
where d^xn ¼ x^n xa , dYn ¼ Ym Yð^xn Þ, and K ¼ ∂F=∂^ x represents the linear or tangent model of the forward model F, which is used to calculate the AIRS observations from the inputs including the atmospheric parameter profiles, surface parameters, and dust properties. In the actual retrieval process, the iteration is stopped when the dYn is less than 0.5 K or the dYn difference from two successive iterations is less than 0.005 K. [23] The variance Sa approaching zero implies that the constraint is a very good estimate of the true parameters and the solution approaches the constraint xa. Conversely, as Sa is approaching infinity, the constraint is a very poor estimate of the true structure and the solution tends toward that specified by the exact inversion of the forward model. [24] Although the observations are sensitive to atmospheric and surface state, dust height, and dust amount as well as its composition, the dust retrieval will be unstable if all parameters are varied simultaneously. Consequently, in this study, the vector x only contains the dust top height and density number (corresponding to the optical depth) that need to be solved, and other inputs for the forward model calculation are assumed to be known. The vector Y m is the observed AIRS brightness temperatures used in the retrieval process. The vector Y(xn) is the calculated AIRS brightness temperatures from a dust state xn by using a forward transfer model. The superscripts n and n + 1 indicate the iteration number and ^x indicates an estimate or retrieved value of x. [25] De Souza-Machado et al. [2010] illustrated that the dust optical depth retrieval is less sensitive to the geometrical thickness of the dust layer than to the dust layer height. To reduce the retrieval uncertainty, the dust geometrical depth is set to five layers, which is about 1 km in the lower troposphere. The initial dust top is set to log (617 hPa) with an uncertainty of 0.1. The initial density number is set to 2000, which result in an optical depth of 1 in the longwave infrared band, with an uncertainty of 2000. The brightness temperature uncertainty from the model and the observation is assumed to be 0.5 K at each selected channel. The matrix Sɛ is a fixed diagonal matrix in which each diagonal element is the square of the assumed uncertainty of the forward model and AIRS observation, while the off-diagonal elements are set to zero. Similarly, in the matrix Sa, each diagonal element is the square of the assumed uncertainty of the initial guess and all off-diagonal elements are set to zero. [26] The optical depth at each channel could be calculated in the RTTOV model from the retrieved density number. As with Peyridieu et al. [2010] and De Souza-Machado et al. [2010], the 10.4 mm wavelength, which roughly corresponds to the maximum of the IR emission from Earth and the absorption of the dust, is chosen as the reference wavelength for infrared optical depth.
4. Simulation Analysis 4.1. Retrieval Accuracy [27] For a simple and robust evaluation of AIRS dust retrievals, a series of simulation experiments are conducted. The simulation analysis could take advantage of the fact that the true values are known. The simulation experiments are conducted as follows: (1) specify true dust parameters xtrue; (2) simulate observations ysim using the forward model; (3) specify the constraint xa, its error covariance Sa and
4 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
D19202
Figure 1. Simulation analysis of dust (a) optical depth and (b) top height retrieval errors from AIRS as a function of the dust optical depth and top height. observation error covariance Sɛ; (4) retrieve xret; (5) calculate the bias xret xtrue. The final retrieval is a weighted average of two quantities: an a priori estimate of parameters and those for which the predicted brightness temperatures from a forward model exactly match the observations. The weights are obtained from the covariance of the error constraint on the one hand and the effective accuracy of the observations on the other hand. In addition to the accuracy of the radiometer itself, effective accuracy depends on the sensitivity of the brightness temperatures to changes in the dust parameters. [28] Figure 1 shows the results of the simulation experiments for different dust top heights and optical depths. For a true vector of dust parameters identical to the constraint values, the retrievals produce identical results with zero bias. As the true vector is allowed to vary above or below the constraint, the retrieved parameters are always biased somewhat toward the constraint because the result is a compromise between the constraint and the observations. Our findings show that better results could be obtained in the case of optical depths greater than 0.5 and dust heights greater than 2.5 km. The worst retrievals occur when the optical depth is close to zero. In this case, the dust signature in the observations could be comparable with the forward and observation errors. When the true dust top height is between 3 and 5 km, the bias in the retrieved dust infrared AOD is less than 25%. 4.2. Retrieval Sensitivity to Other Parameter Errors [29] The retrieval uncertainty not only depends on the sensitivity of TIR observations to dust parameters to be retrieved but also on the error associated with each of these measurements from the instrument itself and from the assumptions of the surface parameters and the atmospheric state. Instrument error primarily results from calibration issues and is about 0.2 K [Chahine et al., 2006]. The assumed error from the forward model used to simulate radiances, however, is generally larger. The remainder of this section is intended to test the algorithm’s ability to estimate dust parameters in the presence of additional
sources of error. In this study, “reference” values (AOD = 1, dust top = 3.93 km (at 617 hPa), surface emissivity = 0.98, atmospheric profiles = tropical climate) are used in the following simulation analysis. [30] The surface emissivity uncertainty is considered as the primary error source for the dust property retrievals, especially for semitransparent dust. Pierangelo et al. [2004] shows that if the dust optical depth is 1.5 at an altitude of 2.4 km, an increase in the surface emissivity form 0.8 to 0.9 changes the observations by 2 K or less. In this study, when the emissivity bias is set to 0.1, the dust top height retrieval error changes by +0.24 km, and the dust optical thickness error changes by 18%. Similarly, when the emissivity bias is 0.02, the retrieval errors for dust height and optical thickness change by +0.05 km and 4%, respectively. Furthermore, when a bias of 2 K is added to the surface temperature, the retrieval errors for dust height and optical thickness change by +0.15 km and 12%, respectively, in the reference case. [31] The errors in temperature profiles from ECMWF reanalysis data should be far less than 2 K, which was reported by Eyre et al. [1993]. The impacts of atmospheric temperature errors on the brightness temperature calculation at these window channels are significantly less than those at sounding channels. Given that the Numerical Weather Prediction (NWP) error for the dust storm is generally worse, a bias of 2 K is added to the temperature profile to investigate the dependence of dust retrievals on the temperature profile errors. The results show that the dust height error changes by 0.65 km and the optical thickness error changes by 3%. It suggests that the temperature profile accuracy has significant impact on the dust top height retrieval. [32] More than 20% of the globe is covered by cirrus clouds [Liou, 1986]. Consequently, dust observations may be contaminated by high-level cirrus. Hong et al. [2006] simulated high spectral resolution infrared observations under cirrus only, dust only, and cirrus and dust combined. Their study shows that the simulated observations are sensitive to cirrus in the case of the low-level dust simultaneously covered by cirrus. In theory, the bias in the dust top
5 of 15
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
D19202
D19202
Figure 2. (a) MODIS composite image, (b) AIRS AOD at 10.4 mm, (c) OMI index, and (d) MODIS DB AOD for the 03/13/2006 dust storm. height decreases as the cirrus approaches the low-level dust. Further uncertainties associated with multilayered dust and the variation in dust composition, although they are certainly important, are beyond the scope of this paper and will not be discussed.
5. Application Results [33] One of the world’s largest shifting-sand deserts and a main source of atmospheric dust in northern China [Xuan and Sokolik, 2002], the Taklimakan Desert lies in the Tarim Basin, between the Tien Shan Mountains in the north and Kunlun Mountains in the south. The basin’s lowest point is roughly 154 m below sea level, and salt collects in the basin due to a lack of drainage (http://www.eoearth.org/ article/Taklimakan_desert). Because of its aridity and abundant sand, this desert is a regular source of dust storms in Asia. To evaluate the retrievals simply and robustly, we have excluded cases for which multiple aerosol layers are detected and measured by the lidar. Three dust storms in the Taklimakan Desert are analyzed in detail by using the 1 DVAR retrieval algorithm in the following study. Additionally, one storm originating from Mongolia is also selected to evaluate the dependence of the retrievals on the dust sources. In the retrieval process, the maximum number of iterations is limited to 10. In the iteration, the log of the dust top pressure must be 0.1 greater than that of the surface pressure, which is about 1 km in height above the surface level. [34] Before retrieving the dust parameters, dust detection is performed by the algorithm described by De SouzaMachado et al. [2010], which uses a summed score from brightness temperature differences from five AIRS channels
(822.4, 900.3, 961.1, 1129.0, and 1231.3 cm1) in the 10 mm atmospheric window region in order to identify dust pixels. In this study, comparisons against visible images show that this algorithm produces few false dust pixels for the selected cases. To further reduce false detection, the pixels with brightness temperatures at 822.4 cm1 greater than their brightness temperatures at 900.3, 961.1, and 1129.0 cm1 are excluded. The results show that these additional criteria could produce more consistent results with the MODIS visible composite images. In the retrieval, the atmospheric profiles and surface temperature are assumed to be known and taken from the ECMWF reanalysis data. The surface infrared emissivity spectrum is taken from the University of Wisconsin-Madison (UW) baseline fit (BLF) emissivity data sets [Seemann et al., 2008]. 5.1. Case 1: Dust Storm of 13 March 2006 [35] On 13 March 2006, a dust and clouds combination lingered over the Taklimakan Desert. In Figure 2a, an Aqua’s MODIS visible composite image shows that the dust storm sweeps over most of the desert. From the northeast toward the center, the desert remains relatively dust free and its floor is clearly shown. White clouds fringe the peripheral edge of the desert. [36] Because OMI, MODIS, and AIRS view the same area within 15 min, the AOD distribution should be very similar. Figures 2b and 2c also show the AIRS retrieved infrared AOD and OMI aerosol index for the dust storm. The heavy dust areas detected by AIRS are consistent with those shown in Figure 2. The aerosol distribution pattern from AIRS agrees well with that from the OMI aerosol index. Additionally, the MODIS DB visible AOD is also shown in
6 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
D19202
Figure 3. Comparison of the AIRS retrieved AOD with (a) the OMI index and (b) the MODIS DB AOD for the 03/13/2006 dust storm. Figure 2d and its distribution pattern also agrees very well with the AIRS AOD retrievals. [37] Figure 3a illustrates the relationship between the infrared AOD from AIRS and the aerosol index from OMI. In the collocation, the distance between the AIRS and OMI pixel centers is less than 7 km, and the AIRS pixels with retrieved AOD differences greater than 0.5 from the adjacent ones are excluded. As expected, a favorable correlation of 0.91 between the AIRS dust thermal optical depth and OMI aerosol index is obtained for this particular dust storm because the two instruments could be used to extract the similar features of the dust storm. It is noted that there are some pixels below the regression line. In the case of low OMI index values, the AIRS underestimation could be caused by the surface parameter errors due to high atmospheric transparence as mentioned above. [38] A clear correlation is also illustrated between the AIRS infrared AOD and the MODIS visible AOD in Figure 3b. The slope of the linear regression between the MODIS visible AOD and the AIRS infrared AOD is about 1.35. It is noted that the ratio between the MODIS AOD and the AIRS AOD increases with the decrease of the AIRS AOD. It means that the relationship between the visible and infrared AOD values varies with the dust optical thickness and the MODIS is more sensitive to the light dust than the AIRS. 5.2. Case 2: Dust Storm of 10 May 2007 [39] A springtime dust storm covering most of the Taklimakan Desert in western China was captured by Aqua’s MODIS at 07:30 UTC on 10 May 2007, which is shown in Figure 4a. In this image, the dust appears thickest in the west and even more concentrated at the western edge of the desert. North of the desert, snow cover caps the Tian Shan Mountains. [40] The vertical cross section along the CALIPSO track is also shown in Figure 4c. It presents the CALIOP backscatter values (km1 sr1) from 0 to 12 km, with the color bar showing the backscatter times 1000. The surface level is represented by the solid black line. For the sake of clarity, CALIOP backscatter data less than 0.75 103 km1 sr1
are masked. The dashed line indicates the AIRS retrieved dust top heights. It is noted that a single dust layer is visible on the right half of the Taklimakan Desert. A thick layer of dust is present from 1 to 2 km, with a tenuous layer between 2 and 4 km. Its top height is about 3 km in the south and then gradually increases to 4 km in the north. It is encouraging that the AIRS retrievals could accurately capture the dust top height variation between 39 and 42 latitude. Also, it is noted that the AIRS retrieved top is significantly underestimated near 38 latitude. This should be caused by the relatively thin aerosol shown in Figure 4d, which is consistent with the simulation results in the case of the AOD less than 0.5 and the aerosol top greater than 5 km in Figure 1b. Between 41 and 42 latitude, the CALIOP dust top is about 1 km greater than the AIRS retrieval because the aerosol near the top is very thin as shown in Figure 4c. [41] For this dust storm the optical depth generally decreases as the dust top height increases. From the horizontal distribution of the dust top height shown in Figure 4b, it is noted that some retrievals are near 5 km, which are significantly larger than neighboring values, thus these pixels could be contaminated by high-level cirrus. The optical depth as a function of latitude along the CALIPSO track is also shown in Figure 4d. There are some discrepancies between the CALIOP retrieved optical depths and those retrieved from the passive instruments, which is similar to the case presented by De SouzaMachado et al. [2010]. However, the variation trend between the AIRS retrieved infrared AOD and the MODIS visible AOD along the track between 38.5 and 41 latitude is consistent. [42] From Figures 4e, 4f, and 4g, it can be seen that the horizontal distribution patterns of the retrieved AIRS AOD and the OMI index as well as the MODIS AOD are very similar. Owing to the relatively higher resolution, the MODIS AOD image could capture detailed variation as shown in the visible composite image, especially in the west of the dust storm. Figure 5b shows the comparisons of the infrared AOD retrievals with the optical products from OMI and MODIS. The correlation coefficients are very similar to those for the first case. It is noted that the offset of the
7 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
Figure 4. (a) MODIS composite image, (b) AIRS retrieved dust top height (km), (c) comparison of AIRS dust top height with the CALIPSO measurements, (d) optical parameters from different measurements (AIRS AOD multiplied by 4 shown in the optical properties along CALIOP track), (e) AIRS AOD at 10.4 mm, (f) OMI aerosol index, and (g) MODIS DB AOD for the 05/10/2007 dust storm. 8 of 15
D19202
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
D19202
Figure 5. Comparison of the AIRS retrieved AOD with (a) the OMI index and (b) MODIS DB AOD for the 05/10/2007 dust storm. regression between MODIS visible and AIRS infrared AOD values is about 0.81. This positive offset means that the visible AOD from MODIS may be overestimated and/or the infrared AOD from AIRS may be underestimated, especially for the light dust areas. The underestimation of the AIRS retrieval may result from the underestimation of the ECMWF surface temperature and/or the surface emissivity used in this study, especially in the case of semitransparent dust, as presented in the previous simulation studies. In general, AOD retrievals are mainly influenced by two factors: surface albedo and aerosol type. At low aerosol loading, the surface albedo is the primary factor causing errors; at high aerosol loading, aerosol model including the radiative transfer model and the OPAC database would affect the AOD retrieval. The first case is a heavy dust storm, while this one is a relatively light dust storm. [43] Figure 6 demonstrates the brightness temperature differences (BTD) between the observed and the simulated measurements before and after the iterations in the retrieval. Before the iteration, the BTDs are dramatically larger in areas with more transparent dust in the middle of the Taklimakan Desert as shown in the MODIS composite image, which means that the initial AOD is too large for these areas.
However, it clearly shows that the iterative process significantly decreases the brightness temperature discrepancies, which could contribute to the final retrieval results. 5.3. Case 3: Dust Storm of 8 April 2007 [44] On 8 April 2007, dust plumes blew out of the Taklimakan Desert toward the east. They are captured by the MODIS image shown in Figure 7a. In this image, the dust seems more concentrated in the east of the desert. White clouds exist along the outside of the desert edge. The CALIOP cross section is also shown in Figure 7. Compared to the two previous storms, this one is relatively more transparent. Along the CALIPSO track, the high cirrus clouds and low level dust are all captured by the lidar, which is shown in Figure 7c. [45] According to the lidar, a dust layer with a top height between 5 km and 4 km and dust base between 2 km and 3 km is present between 38 to 42 latitude. A cirrus layer with geometrical depth ranging from 500 m to 1 km is also present between 39 to 41.5 latitude, separated from the dust layer below by about 6 km. The comparison of the AIRS dust top and the CALIOP backscatters demonstrates that the dust top is retrieved accurately between 39 and 41 latitude, although the high level cirrus, which has a visible
Figure 6. Differences between the observed and the simulated brightness temperatures at 10.4 mm (left) before and (right) after iteration retrievals for the 05/10/2007 dust storm. 9 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
Figure 7. (a) MODIS composite image, (b) AIRS retrieved dust top height (km), (c) comparison of AIRS dust top height with the CALIPSO measurements, (d) cloud and dust optical depths along CALIOP track, (e) AIRS AOD at 10.4 mm, (f) OMI aerosol index, and (g) MODIS DB AOD for the 04/08/2007 dust storm. 10 of 15
D19202
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
D19202
Figure 8. Comparison of the AIRS retrieved AOD with (a) the OMI index and (b) MODIS DB AOD for the 04/08/2007 dust storm. optical depth of less than 0.5 according to Figure 7d, are located above the dust layer with a visible optical depth of greater than 0.5. However, the AIRS dust top near 41.5 latitude is overestimated due to cirrus with a visible optical depth of greater than 0.5 overlapping the dust with a visible optical depth near 0. Between 37 and 38 latitude, the retrieved height seems dominated by the lower level clouds covered by optically thin dust. Additionally, the underestimation near 38.5 latitude may be caused by inaccurate surface and atmospheric parameters used in the retrieval. From Figures 7a and 7b, it is illustrated that the dust detection method could not mask the pixels in the presence of this optically thin cirrus. Some areas covered by thin cirrus are identified as dusty pixels and the retrieved dust top heights are significantly greater than those over the nearby areas. Over this desert, the dust could be lifted up above the surrounding mountains by the strong surface wind accompanying a cold front. As a result, it is difficult to use the negative slope between BT820 and BT960 as well as low BT820 values to indicate thin high cloud over the dust due to the fact that the high level dust could also generate a similar spectrum. Consequently, it is nontrivial to differentiate the high level dust from the high cirrus clouds overlapping a dust layer. [46] From Figures 7e and 7g, it can be seen that the coverage of the dust from AIRS is similar to that from MODIS and the dust variation patterns are also very similar. Figure 7f shows that the heavy dust area from the OMI index agrees well with the results from AIRS. Figure 8 shows the comparisons of the infrared AOD retrievals with the optical products from OMI and MODIS. Although the AIRS infrared AOD correlates well with the MODIS visible AOD and slightly worse with the OMI index, the AIRS AOD and dust top height should be overestimated due to the presence of higher level cirrus, which could be easily seen from the MODIS image, over the northern part of the desert according to the analysis in section 4. The values greater than 3 in the MODIS visible AOD may also be contaminated by the cirrus. The correlation between the results from different instruments will be more complicated to explain due to their different sensitivities to the cirrus overlapping the dust. For
this case, the majority of the AIRS infrared AOD values are less than 1, and the ratio of the MODIS visible AOD and the AIRS infrared AOD is greater than those for the two previous cases. This further confirms that ratio is greater for the light dust than for the heavy dust over this area. 5.4. Case 4: Dust Storm of 27 May 2008 [47] On 27 May 2008, a springtime dust storm spreading from Mongolia through northern China was captured by Aqua’s MODIS image shown in Figure 9a. In this image, dust and clouds form a giant swirl at the borders of Mongolia, Russia, and China. Although clouds mask much of the dust and the arid landscape below, isolated tan streaks mixing with the clouds hover over eastern Mongolia. Dust plumes blow in a counterclockwise direction. Clear skies at southwest of the clouds allow a view of the plumes’ origin. A large amorphous mass of uniform tan formed by the dust exists in the southeast. [48] As the CALIOP backscatter shows in Figure 9c, a dust layer with a top height between 4 km and 3 km and dust base between 1 km and 3 km is present from 41 to 45.5 latitude. The geometrical thickness of the dust layer near 42 latitude reaches 3 km. From this location to 44.5 latitude, the geometrical thickness gradually decreases to 0.5 km. The AIRS retrieved dust top heights agree well with the lidar observations between 41.5 and 44.5 latitude. It seems that the retrieved dust top is insensitive to the bottom of the dust layer. The retrieved results near 45 latitude are about 1 km higher than the lower dust layer top. This difference could be explained by the higher thin dust or cirrus over the dust layer near the height of 3 km. It demonstrates that the algorithm works reasonably for the dust top height retrieval for this particular dust storm originating from Mongolia. [49] For the AOD retrieval, the results reveal that the correlation between the AIRS retrieval and the OMI and MODIS products is remarkably worse than those for the three previous cases. Although the distribution pattern of the AIRS AOD appears to agree well with that of the OMI aerosol index shown in Figure 9e, the poor correlation between the AIRS infrared AOD and the OMI aerosol index
11 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
D19202
Figure 9. (a) MODIS composite image, (b) AIRS retrieved dust top height (km), (c) comparison of AIRS dust top height with the CALIPSO measurements, (d) AIRS AOD at 10.4 mm, (e) OMI aerosol index, and (f) MODIS DB AOD for the 05/27/2008 dust storm. partly resulted from the unreasonable variations between the adjacent observation lines along the satellite track. The significant variations of the horizontal distribution in dust details in the distribution may have been caused by the relatively poor horizontal resolution of the OMI. Additionally, although the heavy dust areas and the detailed variation of the dust storm are clearly shown in the MODIS visible composite image, they are not illustrated well by the MODIS visible AOD image as shown in Figure 9f, in which most of the values is significantly lower than the AIRS retrievals. It appears that the MODIS product is underestimated for this case. [50] Also, it can be seen that the AIRS infrared AOD results are more consistent with the MODIS composite
image than OMI and MODIS aerosol products, especially for the three heavy dust areas. As a result, it is difficult to conclude that the AIRS AOD retrieval is poor for this particular case. Strictly speaking, it is also hard to say that the optical properties from the dust originating from the Taklimakan Desert could be applicable for the dust originating from other areas. 5.5. Summary of Results for All Cases [51] The relation between the AIRS retrieved AOD and the OMI aerosol index is summarized in Table 2. For the three single-layered dust storms over the Taklimakan Desert. The offset and slope of the regression between the AIRS infrared AOD and OMI aerosol index are stable. It is also
12 of 15
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
D19202
D19202
Table 2. Correlation Coefficients (R), Regression Offsets, and Slopes Between AIRS Infrared AOD and OMI Aerosol Index for Different Asian Dust Storms AIRS Infrared AOD Versus OMI Index
AIRS Infrared AOD Versus MODIS DB AOD
Case
Source
R
Offset
Slope
R
Offset
Slope
1 2 3 4
Taklimakan Taklimakan Taklimakan Mongolian
0.91 0.84 0.77 0.54
2.63 2.64 2.57 2.92
1.67 1.62 1.67 0.72
0.95 0.87 0.79 0.35
0.71 0.81 0.81 0.17
1.35 1.45 1.72 0.21
noted that the UV AI is greater than 2 when the infrared AOD is near zero. Actually, the ultraviolet could detect small absorbing aerosol particles, which could not be sensed by the longwave infrared channels. Consequently, some dust pixels with an OMI aerosol index less than 2 could not be correctly identified by using the AIRS longwave infrared band. In other words, this algorithm works well for the optically thick dust storms evaluated in this work but not for tenuous dust. [52] The relation between the AIRS infrared and MODIS DB AOD values is also shown in Table 2. For the three single-layered dust storms, the AIRS infrared AOD is correlated very well to the MODIS DB AOD. The slopes are comparable to the ratio of 1.5–2.5 reported by Peyridieu et al. [2010] for three regions of the tropical North Atlantic and one region of the northwestern Indian Ocean. They also find a slight increase of this ratio when the distance to the sources increases, which may be caused by a small loss in coarse size particles due to gravitational settling. It is also noted that the ratio for the third case is greater than those for the two former cases. From Figures 3b, 5b and 8b, it is shown that the ratio is about 2 when the AIRS AOD is less than 1. It means that the relationship between the AIRS and MODIS DB AODs for the light dust is stable over the Taklimakan desert. With the increase of the AIRS AOD, the ratio will decrease because of the increase of the coarse size particles. [53] Furthermore, cases of lower level dust overlapped by cirrus and double-layered dust are difficult to address due to their complexities. When multiple layers coexist in a vertical column and a single-dust layer assumption is imposed on a retrieval of the dust properties, the resulting dust top tends to be incorrect. The followup effort of the present study will focus on differentiating dust from the cirrus-contaminated observations over the typical dust regions. [54] Although the method could be applied to dust storms occurring over other areas, the quality of the retrievals could not be guaranteed due to the discrepancies in dust composition, which could lead to differences in the absorption and scattering properties. The correlation between the AIRS retrievals and OMI and MODIS optical products significantly degraded for a dust storm originating from Mongolia, which may be partly caused by the different dust composition although the offset is still comparable with those from other cases. 5.6. Impact of Surface Temperature Uncertainty Assumption on Retrievals [55] Simulation analysis shows that surface temperature error has a strong impact on the AIRS AOD retrieval. In
Figure 10. Comparison of the AIRS AOD values from different surface temperature uncertainty assumptions ((a) 2 K, (b) 5 K, and (c) 10 K) for the 04/08/2007 dust storm.
13 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
addition, The ECMWF surface temperature may not be trusted during the day since solar warming of the surface is so variable. Consequently, one may have to simultaneously adjust the surface temperature in the retrieval. To investigate this impact, the experiments are also conducted by adjusting the surface temperature, in which the surface temperature as a variable is added to the vector x described in section 3.5 and other corresponding vectors are also changed easily. In theory, the sensitivity of the AOD retrieval to the surface temperature error increases with the decrease of the dust optical thickness. So, this analysis will mainly focus on the case 3, which is a relatively light dust storm and in which the majority of the visible AOD value are less than 1.5. [56] Comparisons of the AIRS AOD values from different surface temperature uncertainty assumptions against the MODIS DB visible AOD products are shown in Figure 10. It shows that the correlation is gradually improved with the increase of the surface temperature uncertainty. It also clearly reveals that the offset of the linear regression decreases as the surface temperature uncertainty increases. Further analysis shows that the improvements mainly result from the AOD increase in the light dusty pixels, which should be caused by the increase of the surface temperature in the retrieval. The results indicate that the inclusion of the surface temperature as a variable could make the dust AOD retrieval slightly better for this dust storm. For the second case, which is a relatively heavy dust storm, the correlation coefficients of the AOD values do not change for different assumptions as the surface is mainly blacked out under conditions of heavy dust loading.
6. Conclusions and Discussion [57] It is important to detect dust and estimate its properties to evaluate the impact of dust on regional weather and climate systems as well as improve atmospheric parameter retrievals. However, the quality of satellite dust retrievals over land critically depends upon the modeling accuracy of surface properties and dust composition. Also, comparing aerosol data sets from different satellite-based instruments is challenging, and large discrepancies may occur because each retrieval algorithm has its own limitations and shortcomings [Liu and Mishchenko, 2008]. Good correlations between the AIRS infrared AOD and MODIS DB for different one-layered dust storms over the Taklimakan Desert reveal that hyperspectral infrared radiances in the longwave infrared band could be used to retrieve the dust storm optical property over land. The work reported here could help develop an operational dust retrieval algorithm for future AIRS processing. Particularly, it could provide the dust height in the day and night that other passive instruments cannot. The algorithm could also be used to study day/night differences of the optical depth and height of dust close to its sources. Furthermore, the synergetic use of the AIRS and MODIS for the improvement of the dust properties will be conducted in the near future by using the stable relation between the AIRS AOD and MODIS DB. Additionally, the retrievals could be used to validate the dust forecast results. Also, assimilation of the results are expected to help dust forecast. [58] Although the AIRS observations contain abundant dust information and the simulation retrievals by using more channels show that the dust properties could be accurately
D19202
extracted, the retrieval capability from AIRS in shortwave IR channels is limited by the reflected solar radiation. In particular, the simulated observations are generally over 10 K less than the observations in the shortwave band in the daytime. Further experiments also confirm that the iteration could not converge due to very large brightness temperature values in the shortwave band. Consequently, dust properties could not be effectively extracted by using short-wave infrared channels in the day. Additionally, to extract dust information globally, the dust composition and their optical properties over different areas should be investigated by in situ measurements, which will significantly benefit the retrievals. [59] Acknowledgments. The authors would like to thank Timothy J. Schmit for his suggestions on improving the manuscript. Elisabeth Weisz is also thanked for her help on using AIRS data. This work is partly supported by basic research program 2010CB950802, the National Oceanic and Atmospheric Administration (NOAA) GOES-R Algorithm Working Group (AWG) and GOES-R Risk Reduction programs NA10NES4400010 and the National Natural Science Foundation of China (NSFC) grant 41075021. Hyo-jin Han and B. J. Sohn were supported by the National Space Lab program funded by the National Research Foundation of Korea (S10801000184-08A0100-18410) and the Korea Meteorological Administration Research and Development Program under grant CATER 20122060. The authors would like to specifically thank the AIRS, MODIS, OMI, and CALIPSO science teams for providing satellite data and ECMWF for sharing the model output data. We appreciate the anonymous reviewers for their valuable comments and suggestions to improve the quality of this paper.
References Ackerman, S., H.-L. Huang, P. Antonelli, R. Holz, H. Revercomb, D. Tobin, K. Baggett, and J. E. Davies (2004), Detection of clouds and aerosols using infrared hyperspectral observations, paper presented at 20th International Conference on Interactive Information and Processing Systems for Meteorology, Oceanography, and Hydrology, Am. Meteorol. Soc., Seattle, Wash. Brinksma, E. J., et al. (2008), The 2005 and 2006 DANDELIONS NO2 and aerosol intercomparison campaigns, J. Geophys. Res., 113, D16S46, doi:10.1029/2007JD008808. Chahine, M. T., et al. (2006), AIRS: Improving weather forecasting and providing new data on greenhouse gases, Bull. Am. Meteorol. Soc., 87, 911–926, doi:10.1175/BAMS-87-7-911. Claquin, T., M. Schulz, and Y. J. Balkanski (1999), Modeling the mineralogy of atmospheric dust sources, J. Geophys. Res., 104, 22,243–22,256, doi:10.1029/1999JD900416. Curier, L., J. P. Veefkind, R. Braak, O. Torres, B. Veihelmann, and G. de Leeuw (2008), Retrieval of aerosols optical properties from OMI radiances using a multi-wavelength algorithm: Application and validation to western Europe, J. Geophys. Res., 113, D17S90, doi:10.1029/ 2007JD008738. De Souza‐Machado, S., L. L. Strow, H. Motteler, and S. Hannon (2006), Infrared dust spectral signatures from AIRS, Geophys. Res. Lett., 33, L03801, doi:10.1029/2005GL024364. De Souza-Machado, S. G., et al. (2010), Infrared retrievals of dust using AIRS: Comparisons of optical depths and heights derived for a North African dust storm to other collocated EOS A-Train and surface observations, J. Geophys. Res., 115, D15201, doi:10.1029/2009JD012842. Díaz, J., F. Expósito, C. Torres, V. Carreo, and A. Redondas (2000), Simulation of mineral dust effects on UV radiation levels, J. Geophys. Res., 105, 4979–4991, doi:10.1029/1999JD901058. Eguchi, K., I. Uno, K. Yumimoto, T. Takemura, A. Shimizu, N. Sugimoto, and Z. Liu (2009), Trans-Pacific dust transport: Integrated analysis of NASA/CALIPSO and a global aerosol transport model, Atmos. Chem. Phys., 9, 3137–3145, doi:10.5194/acp-9-3137-2009. Eyre, J. R. (1989), Inversion of cloudy satellite sounding radiances by nonlinear optimal estimation. I: Theory and simulation, Q. J. R. Meteorol. Soc., 115, 1001–1026, doi:10.1002/qj.49711548902. Eyre, J. R. (1991), A fast radiative transfer model for satellite sounding systems, Res. Dep. Tech. Memo. 176, Eur. Cent. for Medium-Range Weather Forecast., Reading, U. K. Eyre, J., G. Kelly, A. McNally, E. Anderson, and A. Persson (1993), Assimilation of TOVS radiance information through one-dimensional variational analysis, Q. J. R. Meteorol. Soc., 119, 1427–1463, doi:10.1002/ qj.49711951411.
14 of 15
D19202
YAO ET AL.: DUST RETRIEVALS OVER LAND FROM AIRS
Fischer, E. V., N. C. Hsu, D. A. Jaffe, M.-J. Jeong, and S. L. Gong (2009), A decade of dust: Asian dust and springtime aerosol load in the U.S. Pacific Northwest, Geophys. Res. Lett., 36, L03821, doi:10.1029/ 2008GL036467. Han, H.-J., B. J. Sohn, A. A. Huang, E. Weisz, R. Saunders, and T. Takamura (2012), An improved radiance simulation for hyperspectral infrared remote sensing of Asian dust, J. Geophys. Res., 117, D09211, doi:10.1029/ 2012JD017466. Hansen, J., M. Sato, and R. Ruedy (1997), Radiative forcing and climate response, J. Geophys. Res., 102, 6831–6864, doi:10.1029/96JD03436. Hong, G., P. Yang, H.-L. Huang, S. A. Ackerman, and I. N. Sokolik (2006), Simulation of high-spectral-resolution infrared signature of overlapping cirrus clouds and mineral dust, Geophys. Res. Lett., 33, L04805, doi:10.1029/2005GL024381. Hsu, N. C., S. C. Tsay, M. D. King, and J. R. Herman (2004), Aerosol properties over bright-reflecting source regions, IEEE Trans. Geosci. Remote Sens., 42, 557–569, doi:10.1109/TGRS.2004.824067. Hsu, N. C., S. Tsay, M. D. King, and J. R. Herman (2006), Deep blue retrievals of Asian aerosol properties during ACE-Asia, IEEE Trans. Geosci. Remote Sens., 44(11), 3180–3195, doi:10.1109/TGRS.2006.879540. Intergovernmental Panel on Climate Change (2007), Climate Change 2007—The Fourth Assessment Report of the IPCC, edited by S. Solomon et al., Cambridge Univ. Press, Cambridge, U. K. Jickells, T., et al. (2005), Global iron connections between desert dust, ocean biogeochemistry and climate, Science, 308, 67–71, doi:10.1126/ science.1105959. Kaufman, Y. J., D. Tanré, L. A. Remer, E. F. Vermote, A. Chu, and B. N. Holben (1997a), Operational remote sensing of tropospheric aerosol over land from EOS moderate resolution imaging spectroradiometer, J. Geophys. Res., 102, 17–51. Kaufman, Y. J., A. E. Wald, L. A. Remer, B.-C. Gao, R.-R. Li, and L. Flynn (1997b), The MODIS 2.1-mm channel-correlation with visible reflectance for use in remote sensing of aerosol, IEEE Trans. Geosci. Remote Sens., 35, 1286–1298, doi:10.1109/36.628795. Kaufman, Y. J., D. Tanre, and O. Boucher (2002), A satellite view of aerosols in climate systems, Nature, 419, 215–223, doi:10.1038/nature01091. Kim, D.-H., B.-J. Sohn, T. Nakajima, T. Takamura, T. Takemura, B.-C. Choi, and S.-C. Yoon (2004), Aerosol optical properties over East Asia determined from ground-based sky radiation measurements, J. Geophys. Res., 109, D02209, doi:10.1029/2003JD003387. Kim, D.-H., B.-J. Sohn, T. Nakajima, and T. Takamura (2005), Aerosol radiative forcing over East Asia determined from ground-based solar radiation measurements, J. Geophys. Res., 110, D10S22, doi:10.1029/ 2004JD004678. Lee, E.-H., and B. J. Sohn (2011), Recent increase of dust trend over Mongolia and Inner Mongolia regions and its association with changed surface climate conditions, Atmos. Environ., 45(27), 4611–4616. Levelt, P. F., E. Hilsenrath, G. W. Leppelmeier, G. H. J. van den Ooord, P. K. Bhartia, J. Taminnen, J. F. de Haan, and J. P. Veefkind (2006), Science objectives of the Ozone Monitoring Instrument, IEEE Trans. Geosci. Remote Sens., 44(5), 1093–1101, doi:10.1109/TGRS.2006.872333. Li, J., et al. (2005), Retrieval of cloud microphysical properties from MODIS and AIRS, J. Appl. Meteorol., 44, 1526–1543, doi:10.1175/ JAM2281.1. Li, X., et al. (2012), Validation of MODIS and Deep Blue aerosol optical depth retrievals in an arid/semi-arid region of northwest China, Particuology, 10, 132–139, doi:10.1016/j.partic.2011.08.002. Li, Z., J. Li, X. Jin, T. J. Schmit, E. E. Borbas, and M. D. Goldberg (2010), An objective methodology for infrared land surface emissivity evaluation, J. Geophys. Res., 115, D22308, doi:10.1029/2010JD014249. Liou, K.-N. (1986), Influence of cirrus clouds on weather and climate processes: A global perspective, Mon. Weather Rev., 114, 1167–1199, doi:10.1175/1520-0493(1986)1142.0.CO;2. Liu, L., and M. Mishchenko (2008), Toward unified satellite climatology of aerosol properties: Direct comparisons of advanced level 2 aerosol
D19202
products, J. Quant. Spectrosc. Radiat. Transfer, 109, 2376–2385, doi:10.1016/j.jqsrt.2008.05.003. Lubin, D., S. Satheesh, G. McFarquar, and A. Heymsfield (2002), Longwave radiative forcing of Indian Ocean tropospheric aerosol, J. Geophys. Res., 107(D19), 8004, doi:10.1029/2001JD001183. Peyridieu, S., A. Chedin, D. Tanre, V. Capelle, C. Pierangelo, N. Lamquin, and R. Armante (2010), Saharan dust infrared optical depth and altitude retrieved from AIRS: A focus over North Atlantic—Comparison to MODIS and CALIPSO, Atmos. Chem. Phys., 10, 1953–1967, doi:10.5194/acp-10-1953-2010. Pierangelo, C., A. Chedin, S. Heilliette, N. Jacquinet-Husson, and R. Armante (2004), Dust altitude and infrared optical depth from AIRS, Atmos. Chem. Phys., 4, 1813–1822, doi:10.5194/acp-4-1813-2004. Rodgers, C. D. (1976), Retrieval of atmospheric temperature and composition from remote measurements of thermal radiation, Rev. Geophys., 14, 609–624, doi:10.1029/RG014i004p00609. Seemann, S. W., E. E. Borbas, R. O. Knuteson, G. R. Stephenson, and H.-L. Huang (2008), Development of a global infrared land surface emissivity database for application to clear sky sounding retrievals from multispectral satellite radiance measurements, J. Appl. Meteorol. Climatol., 47, 108–123, doi:10.1175/2007JAMC1590.1. Seinfeld, H. J., et al. (2004), Regional climatic and atmospheric chemical effects of Asian dust and pollution, Bull. Am. Meteorol. Soc., 85(3), 367–380, doi:10.1175/BAMS-85-3-367. Sokolik, I. (2002), The spectral radiative signature of wind-blown mineral dust: Implications for remote sensing in the thermal IR region, Geophys. Res. Lett., 29(24), 2154, doi:10.1029/2002GL015910. Stephens, G. L., et al. (2002), The CloudSat mission and the A-Train: A new dimension to space-based observations of clouds and precipitation, Bull. Am. Meteorol. Soc., 83, 1771–1790, doi:10.1175/BAMS-83-121771. Su, L., and O. B. Toon (2010), Saharan and Asian dust: Similarities and differences determined by CALIPSO, AERONET and a coupled climateaerosol microphysical model, Atmos. Chem. Phys. Discuss., 10, 29,513–29,567, doi:10.5194/acpd-10-29513-2010. Sun, J., M. Zhang, and T. Liu (2001), Spatial and temporal characteristics of dust storms in China and its surrounding regions, 1960–1999: Relations to source area and climate, J. Geophys. Res., 106, 10,325–10,333. Torres, O., P. K. Bhartia, J. R. Herman, and Z. Ahmad (1998), Derivation of aerosol properties from satellite measurements of backscattered ultraviolet radiation: Theoretical basis, J. Geophys. Res., 103, 17,099–17,110, doi:10.1029/98JD00900. Torres, O., A. Tanskanen, B. Veihelmann, C. Ahn, R. Braak, P. Bhartia, P. Veefkind, and P. Levelt (2007), Aerosols and surface UV products from Ozone Monitoring Instrument observations: An overview, J. Geophys. Res., 112, D24S47, doi:10.1029/2007JD008809. Uno, I., K. Eguchi, K. Yumimoto, T. Takemura, A. Shimizu, M. Uematsu, Z. Liu, Z. Wang, Y. Hara, and N. Sugimoto (2009), Asian dust transported one full circuit around the globe, Nat. Geosci., 2, 557–560, doi:10.1038/ngeo583. Wald, A., Y. Kaufman, D. Tanré, and B.-C. Gao (1998), Daytime and nighttime detection of mineral dust over desert using infrared spectral contrast, J. Geophys. Res., 103, 32,307–32,313, doi:10.1029/98JD01454. Weisz, E., J. Li, J. Li, D. K. Zhou, H.-L. Huang, M. D. Goldberg, and P. Yang (2007), Cloudy sounding and cloud-top height retrieval from AIRS alone single field-of-view radiance measurements, Geophys. Res. Lett., 34, L12802, doi:10.1029/2007GL030219. Winker, D. M., B. H. Hunt, and M. J. McGill (2007), Initial performance assessment of CALIOP, Geophys. Res. Lett., 34, L19803, doi:10.1029/ 2007GL030135. Xuan, J., and I. N. Sokolik (2002), Sources of wind-blown dust in Northern China: Towards developing an Asian Dust Databank, in Proceedings of ICAR5/GCTE-SEN Joint Conference, USDA Pub. 2–02, edited by J. A. Lee and T. M. Zobeck, pp. 357–361, Int. Cent. for Arid and Semiarid Lands Stud., Tex. Tech Univ., Lubbock.
15 of 15