Retrieval of Aerosol Optical Depth over Arid Areas from MODIS Data

Moderate Resolution Imaging Spectroradiometer (MODIS) data have been widely applied for the remote sensing of aerosol optical depth (AOD) because the MODIS sensor features a short revisit period and a moderate spatial resolution. The Dense Dark Vegetation (DDV) method is the most popular retrieval method. However, the DDV method can only be used to retrieve the AOD with high precision when the surface reflectance in the visible spectrum is low, such as over dense vegetation or water. To obtain precise AOD values in areas with higher reflectance, such as arid areas, Land Surface Reflectance (LSR) must be estimated accurately. This paper proposes a method of estimating LSR for AOD retrieval over arid areas from long-term series of MODIS images. According to the atmospheric parameters (AOD and water vapor), the clearest image without clouds was selected from the long-term series of continuous MODIS images. Atmospheric correction was conducted based on similar ground-measured atmospheric parameters and was used to estimate the LSR and retrieve the AOD at adjacent times. To validate this method, aerosol inversion experiments were performed in northern Xinjiang, in which the inverted AOD was compared to ground-measured AOD and MODIS aerosol products (MOD04). The AOD retrieved using the new algorithm was highly consistent with the AOD derived from ground-based measurements, with a correlation coefficient of 0.84. Additionally, 82.22% of the points fell within the expected error defined by NASA. The precision of the retrieved AOD data was better than that of MOD04 AOD products over arid areas.


Introduction
Because of their excellent abilities to absorb and scatter short-wavelength solar radiation, aerosols are important in the radiant energy balance of the earth-atmosphere system and greatly influence regional weather variations and global climate change. In addition, because aerosols originate from diverse sources, have complex compositions, and are affected by human and natural activities, determination of their optical properties and spatial-temporal variations is difficult.
The use of inversion methods in aerosol remote sensing dates back to the 1970s, and several inversion algorithms have been developed. The first algorithms mainly focused on areas with low surface reflectance. For example, Griggs [1] and Mekler et al. [2] determined the distribution of aerosols above terrestrial water areas using Earth Resources Technology Satellite (ERTS) data; Griggs [3] studied aerosols over the Mediterranean Sea and Atlantic Ocean using National Oceanic and Atmospheric Administration (NOAA) 5 and Geostationary Operational Environmental Satellites (GOES) data; and Nagaraja et al. [4] applied an aerosol operational method for areas over oceans using NOAA's Advanced Very High Resolution Radiometer (AVHRR) data. Based on previous research, Durkee et al. [5] proposed a multi-channel aerosol retrieval algorithm with improved precision.

Fundamentals
Aerosol remote sensing over land is usually performed using the visible waveband, which has a shorter wavelength and is sensitive to aerosol scattering. In the visible waveband, the top of the atmosphere (TOA) radiance L measured by a satellite sensor is a function of the solar and view zenith and azimuth angles, and it can be estimated using Equation (1): where τ a is the atmospheric optical depth resulting from gas molecules and aerosol particles; L 0 is the path radiance that is produced by scattering of the particles and gas molecules in the atmosphere; ρ s is the surface reflectance; S(τ a ) is an atmospheric backscattering ratio to account for the multiple reflection between the surface and atmosphere; T(τ a , µ s ) is the transmittance from the ground to the TOA; T(τ a , µ v ) is the transmittance from the top to the ground; µ s and µ v are the cosine of the solar zenith angle and the cosine of the view zenith angle, respectively; and φ is the relative azimuth angle between the solar and view angles. The energy sources of the apparent reflectance or radiance values acquired by the satellite sensor can be divided into two parts: the purely atmospheric contribution, which is addressed in the first part of Equation (1), and a combined atmosphere and land surface contribution, which is addressed in the second part of Equation (1). The first part of Equation (1) corresponds to small surface reflectance values because aerosols mainly influence the electromagnetic waves in the visible waveband; thus, the reflectance values are sensitive to aerosol variations. This relationship also explains why the identification of aerosols from remote sensing data is more accurate over oceans and why aerosol data collection from the red and blue bands improves with increasing vegetation coverage over land. Figure 1 shows the spectral curves obtained from the visible to near-infrared (NIR) bands of several typical land features (vegetation, sand, soil, water, and rock) extracted from the USGS spectral library. In the visible waveband, water and vegetation have lower reflectance values, with a reflectance <0.08 for water in all the experimental bands. In addition, the reflectance of healthy vegetation is <0.06 in the red and blue bands, while soil and sand have relatively higher reflectance values of 0.10 and 0. 17, respectively. grows slowly in arid areas; and (2) the directional reflectance has little impact because arid areas have low vegetation coverage.

Fundamentals
Aerosol remote sensing over land is usually performed using the visible waveband, which has a shorter wavelength and is sensitive to aerosol scattering. In the visible waveband, the top of the atmosphere (TOA) radiance L measured by a satellite sensor is a function of the solar and view zenith and azimuth angles, and it can be estimated using Equation (1): where a  is the atmospheric optical depth resulting from gas molecules and aerosol particles; 0 L is the path radiance that is produced by scattering of the particles and gas molecules in the atmosphere; s  is the surface reflectance; ( ) a S  is an atmospheric backscattering ratio to account for the multiple reflection between the surface and atmosphere; ( , ) a s T   is the transmittance from the ground to the TOA; ( , ) a v T   is the transmittance from the top to the ground; s  and v  are the cosine of the solar zenith angle and the cosine of the view zenith angle, respectively; and  is the relative azimuth angle between the solar and view angles. The energy sources of the apparent reflectance or radiance values acquired by the satellite sensor can be divided into two parts: the purely atmospheric contribution, which is addressed in the first part of Equation (1), and a combined atmosphere and land surface contribution, which is addressed in the second part of Equation (1). The first part of Equation (1) corresponds to small surface reflectance values because aerosols mainly influence the electromagnetic waves in the visible waveband; thus, the reflectance values are sensitive to aerosol variations. This relationship also explains why the identification of aerosols from remote sensing data is more accurate over oceans and why aerosol data collection from the red and blue bands improves with increasing vegetation coverage over land. Figure 1 shows the spectral curves obtained from the visible to near-infrared (NIR) bands of several typical land features (vegetation, sand, soil, water, and rock) extracted from the USGS spectral library. In the visible waveband, water and vegetation have lower reflectance values, with a reflectance < 0.08 for water in all the experimental bands. In addition, the reflectance of healthy vegetation is < 0.06 in the red and blue bands, while soil and sand have relatively higher reflectance values of 0.10 and 0.17, respectively.  The 6S (Second Simulation of the Satellite Signal in the Solar Spectrum) basic radiative transfer model is used to simulate satellite observations that consider elevated targets, molecular and aerosol scattering and gaseous absorption (including H 2 O, O 3 , CO 2 , etc.) [21]. Figure 2 shows the simulation The surface reflectance range of the red and blue bands is 0.001-0.20, and this range includes most typical land features, excluding ice and snow; and (4) the 6S model defines the commonly used aerosol models, including urban/industrial, continental, desert, and biomass-burning aerosol models. These parameters were chosen to simulate the sensitivity of the AOD, and the results are shown in Figure 2.
As shown in Figure 2, the urban/industrial aerosol model exhibited a significant damping effect because the Ångström index was generally larger than 1.1 [22], and the aerosol particles influence the propagation of electromagnetic waves [23]. Because this study considers the arid area of Northwest China, the main type of land cover is desert, with a small amount of farmland. The analysis results from desert, continental and biomass-burning aerosol models are shown in Figure 2. When the surface reflectance of the blue band is high (~0.20), the slope of the simulated curve becomes smaller (Figure 2A,C); thus, the satellite apparent reflectance is not sensitive to changes in the AOD. A similar situation occurs when the surface reflectance in the red band is greater than 0.15 ( Figure 2B,D). In other words, when the surface reflectance of the blue or red band is smaller than a certain threshold value, the aerosol path radiance and the radiance of the land surface can be identified in the satellite signal, but their degrees are lower than the low surface reflectance. To guarantee the accuracy of AOD retrieval, and based on the surface reflectance distribution in Figure 1, only the blue band reflectance values less than 0.12 and the red band reflectance values less than 0.15 were considered. The 6S (Second Simulation of the Satellite Signal in the Solar Spectrum) basic radiative transfer model is used to simulate satellite observations that consider elevated targets, molecular and aerosol scattering and gaseous absorption (including H2O, O3, CO2, etc.) [21]. Figure 2 shows the simulation results of the sensitivity of the apparent reflectance and AOD under different conditions: (1) the spectral parameters are the MODIS blue and red bands; (2) the geometric parameters are 15 , 0 , 100 The surface reflectance range of the red and blue bands is 0.001-0.20, and this range includes most typical land features, excluding ice and snow; and (4) the 6S model defines the commonly used aerosol models, including urban/industrial, continental, desert, and biomass-burning aerosol models. These parameters were chosen to simulate the sensitivity of the AOD, and the results are shown in Figure 2. As shown in Figure 2, the urban/industrial aerosol model exhibited a significant damping effect because the Ångström index was generally larger than 1.1 [22], and the aerosol particles influence the propagation of electromagnetic waves [23]. Because this study considers the arid area of Northwest China, the main type of land cover is desert, with a small amount of farmland. The analysis results from desert, continental and biomass-burning aerosol models are shown in Figure 2 When the surface reflectance of the blue band is high (~0.20), the slope of the simulated curve becomes smaller (Figure 2 A,C); thus, the satellite apparent reflectance is not sensitive to changes in the AOD. A similar situation occurs when the surface reflectance in the red band is greater than 0.15 (Figure 2 B,D). In other words, when the surface reflectance of the blue or red band is smaller than a certain threshold value, the aerosol path radiance and the radiance of the land surface can be identified in the satellite signal, but their degrees are lower than the low surface reflectance. To guarantee the accuracy of AOD retrieval, and based on the surface reflectance distribution in Figure 1, only the blue band reflectance values less than 0.12 and the red band reflectance values less than 0.15 were considered. The MODIS DDV algorithm over land primarily uses three spectral channels centered at wavelengths of 0.47, 0.66, and 2.13 μm. The surface reflectance of the red and blue bands is obtained by two empirical relationships [9]  , respectively. In addition, areas beyond this range cannot be retrieved with high precision because the ability to identify aerosols from the radiance (or apparent reflectance) decreases, making it difficult to determine the surface reflectance [14].

Inversion Algorithm
A clear image was selected from a series of MODIS data taken at adjacent times to provide surface reflectance information. Clear images were selected based on ground measurements with small AODs, few clouds and stable atmospheric conditions. The observation sites were reasonably distributed throughout the research region, and the AOD data were obtained from the images via interpolation. Then, after atmospheric correction, the surface reflectance was acquired and applied to AOD remote sensing using other images after spatial matching. The surface reflectance was calculated using Equation (2): After acquiring the surface reflectance using Equation (2) and substituting it into Equation (1), all atmospheric contributions, including the atmospheric aerosols and molecules within the visible The MODIS DDV algorithm over land primarily uses three spectral channels centered at wavelengths of 0.47, 0.66, and 2.13 µm. The surface reflectance of the red and blue bands is obtained by two empirical relationships [9] (ρ s 0.66 = ρ 2.13 /2 and ρ s 0.47 = ρ 2.13 /4) developed for use over vegetated surfaces. In addition, to minimize the error, ρ 2.13 was limited to less than 0.25. Therefore, the surface reflectance ranges in the red and blue bands were limited to ρ s 0.66 ≤ 0.125 and ρ s 0.47 ≤ 0.0625, respectively. In addition, areas beyond this range cannot be retrieved with high precision because the ability to identify aerosols from the radiance (or apparent reflectance) decreases, making it difficult to determine the surface reflectance [14].

Inversion Algorithm
A clear image was selected from a series of MODIS data taken at adjacent times to provide surface reflectance information. Clear images were selected based on ground measurements with small AODs, few clouds and stable atmospheric conditions. The observation sites were reasonably distributed throughout the research region, and the AOD data were obtained from the images via interpolation. Then, after atmospheric correction, the surface reflectance was acquired and applied to AOD remote sensing using other images after spatial matching. The surface reflectance was calculated using Equation (2): After acquiring the surface reflectance using Equation (2) and substituting it into Equation (1), all atmospheric contributions, including the atmospheric aerosols and molecules within the visible wavelength range, were obtained, where the Rayleigh scattering of the atmospheric molecules was relatively important at longer wavelengths [24]. To extract the aerosol information, the atmospheric molecular Rayleigh scattering effect was removed using Equation (3), which results in the AOD (τ Aero ): At sea level, the Rayleigh optical thickness (ROT) over the visible light wavelength can be estimated as a function of the wavelength λ [25]: At a given height above sea level, the ROT can be determined as a function of atmospheric pressure or elevation [26]: where z is the elevation of the ground above sea level in kilometers (km) and 8.5 is the exponential scale height of the atmosphere. This paper utilizes high-temporal resolution MODIS data. Based on the atmospheric parameters from ground-measured data, after selecting a clear image from the MODIS images, the surface reflectance can be obtained using the atmospheric corrections described above. Then, the aerosol optical parameters (the single scattering albedo is 0.827, and the asymmetry factor is 0.910) in the study area are retrieved using sun photometric data, and the 6S model is used to construct a Look-Up table (LUT). To eliminate the influence of Rayleigh scattering on atmospheric molecules, a digital elevation model (DEM) based on the MODIS geolocation product data (MOD03) was used to estimate the height (z) and calculate τ Ray for each pixel. A Rayleigh correction was applied only to land pixels, and water pixels were masked based on the land mask data from the MOD03 product. The flow chart in Figure 3 shows the procedure for retrieving AOD from MODIS data.
Atmosphere 2016, 7, 134 6 of 16 wavelength range, were obtained, where the Rayleigh scattering of the atmospheric molecules was relatively important at longer wavelengths [24]. To extract the aerosol information, the atmospheric molecular Rayleigh scattering effect was removed using Equation (3), which results in the AOD ( Aero At sea level, the Rayleigh optical thickness (ROT) over the visible light wavelength can be estimated as a function of the wavelength  [25]: At a given height above sea level, the ROT can be determined as a function of atmospheric pressure or elevation [26]: where z is the elevation of the ground above sea level in kilometers (km) and 8.5 is the exponential scale height of the atmosphere. This paper utilizes high-temporal resolution MODIS data. Based on the atmospheric parameters from ground-measured data, after selecting a clear image from the MODIS images, the surface reflectance can be obtained using the atmospheric corrections described above. Then, the aerosol optical parameters (the single scattering albedo is 0.827, and the asymmetry factor is 0.910) in the study area are retrieved using sun photometric data, and the 6S model is used to construct a Look-Up table (LUT). To eliminate the influence of Rayleigh scattering on atmospheric molecules, a digital elevation model (DEM) based on the MODIS geolocation product data (MOD03) was used to estimate the height (z) and calculate Ray  for each pixel. A Rayleigh correction was applied only to land pixels, and water pixels were masked based on the land mask data from the MOD03 product. The flow chart in Figure 3 shows the procedure for retrieving AOD from MODIS data.

Look-Up Tables
To quickly and conveniently calculate AOD using the radiative transfer equation, an LUT was constructed using the 6S model. All of these methods rely on LUTs and require several parameters, such as observation geometry, atmospheric conditions, aerosol parameters, and surface reflectance.
(1) Observation geometry. During satellite imaging, the geometric relationship between the sun, target, and satellite can be described by the solar zenith angle, solar azimuth, viewing zenith angle, and viewing azimuth. The sun and viewing azimuths can be described by their relative azimuths. In addition, based on a simulation of the analytical sensitivity between the angle and satellite apparent reflectance provided by the 6S model, three angle parameters were set: the solar zenith angle and viewing zenith angle were set 0 • -80 • with a step size of 5 • and 17 parameters, and the relative azimuth angle was set from 0 • to 180 • with a step size of 10 • and 19 parameters.

Data Source and Experimental Area
MODIS sensors are located on the sun-synchronous polar satellites Terra and Aqua, which pass over the study area at 10:30 a.m. and 1:30 p.m. local time, respectively. Daily coverage is achieved at high latitudes (above 30 • ). MODIS has 36 spectral channels that obtain images in the visible and thermal infrared bands (0.4-14 µm) with spatial resolutions of 1 km, 500 m and 250 m. MODIS Level 1B (L1B) data were applied in this inversion experiment.
The autonomous territory of Xinjiang province is located in the hinterland of the Eurasian continent, far from the ocean, and it is dominated by a typical continental arid climate. Because it has low precipitation, sparse vegetation, poor soil, and fragile ecosystems, this region is highly sensitive to various natural and human activities. Ground-measured AOD data were used for atmospheric correction to obtain clear images and validate the inversion method. Continuous observation experiments were performed for two months (July and August 2014) at the Dahuangshan and Wucaiwan sites (the specific locations are shown in Figure 4).
The measurement instruments used at the Wucaiwan and Dahuangshan sites included a CE318 sun photometer produced by the CIMEL Co. and a Microtops II hand-held sun photometer produced by the Solar Light Company. The CE318 sun photometer automatically tracks and measures direct solar and sky radiation, which are used to calculate atmospheric transmissivity, AOD, and ozone content. The CE318 has 8 channels and provides cloud-free AOD observations every 15 min, with an uncertainty of~0.01-0.02 [27]. The Microtops II sun photometer has five measurement channels and is easy to operate and transport, with an uncertainty of~0.015-0.02 [28]. The band information for the above two instruments is shown in Table 1. Before the experiment, we conducted a comparison study at the Dahuangshan site from 7 July to 10 July 2014, with the two instruments and observed a maximum absolute error of 0.005 when the AOD at 500 nm varied from 0.14-1.63. The observation time in the study area spanned from 10 July to 24 August 2014, during which observations of several atmospheric parameters were made, including AOD, direct solar radiation, and atmospheric moisture content.

Surface Reflectance
Based on the ground measurement data (shown in Figure 5), the AODs measured by the CE318 and Microtops II sun photometers were closest on 30 July 2014, with an absolute difference of 0.0146. The imaging time of the MODIS sensor on this date was 13:35 within one hour before and after imaging, and the maximum and minimum absolute deviations of the AOD values at 500 nm obtained from the two instruments were 0.0234 and 0.0011, respectively. The average water vapor measured by CE318 during this time period was 1.4804 cm. No cloud cover occurred over the study area during this time. Thus, the atmospheric conditions were relatively stable over the study area. Therefore, we selected the 30 July 2014, MODIS image as a clear reference image from the long-term series of MODIS data. Using the measured atmospheric parameters (AOD = 0.2435 and water vapor = 1.4804 cm), atmospheric correction was conducted using the radiative transfer model to obtain the surface reflectance of the red and blue bands. The results are shown in Figure 6. The surface reflectance in the blue band was generally smaller than that in the red band. The local desert   1  340  2  1  380  4  2  380  2  2  500  10  3  440  10  3  870  10  4  500  10  4  936  10  5  670  10  5  1020  10  6  870  10  7  1020  10  8 1640 60

Surface Reflectance
Based on the ground measurement data (shown in Figure 5), the AODs measured by the CE318 and Microtops II sun photometers were closest on 30 July 2014, with an absolute difference of 0.0146. The imaging time of the MODIS sensor on this date was 13:35 within one hour before and after imaging, and the maximum and minimum absolute deviations of the AOD values at 500 nm obtained from the two instruments were 0.0234 and 0.0011, respectively. The average water vapor measured by CE318 during this time period was 1.4804 cm. No cloud cover occurred over the study area during this time. Thus, the atmospheric conditions were relatively stable over the study area. Therefore, we selected the 30 July 2014, MODIS image as a clear reference image from the long-term series of MODIS data. Using the measured atmospheric parameters (AOD = 0.2435 and water vapor = 1.4804 cm), atmospheric correction was conducted using the radiative transfer model to obtain the surface Atmosphere 2016, 7, 134 9 of 16 reflectance of the red and blue bands. The results are shown in Figure 6. The surface reflectance in the blue band was generally smaller than that in the red band. The local desert area and the ice-covered area on top of the Tianshan Mountains had reflectance values greater than 0.2. Figure 7 shows a cumulative percent distribution diagram of the surface reflectance according to the atmospheric correction results, the blue and red lines represent the surface reflectance distribution of blue and red bands, respectively. As shown in Figure 7, approximately 99% of the surface reflectance values were below 0.12 in the blue band, and almost 58% of the surface reflectance values were below 0.15 in the red band. As illustrated in Figure 2, when the surface reflectance in the blue and red bands is less than 0.12 or 0.15, respectively, changes in AOD can still cause changes in apparent reflectance. In addition, to minimize the error, the surface reflectance values of the blue and red bands were limited to less than 0.12 and 0.15, respectively.  Figure 7 shows a cumulative percent distribution diagram of the surface reflectance according to the atmospheric correction results, the blue and red lines represent the surface reflectance distribution of blue and red bands, respectively. As shown in Figure 7, approximately 99% of the surface reflectance values were below 0.12 in the blue band, and almost 58% of the surface reflectance values were below 0.15 in the red band. As illustrated in Figure 2, when the surface reflectance in the blue and red bands is less than 0.12 or 0.15, respectively, changes in AOD can still cause changes in apparent reflectance. In addition, to minimize the error, the surface reflectance values of the blue and red bands were limited to less than 0.12 and 0.15, respectively.   Figure 7 shows a cumulative percent distribution diagram of the surface reflectance according to the atmospheric correction results, the blue and red lines represent the surface reflectance distribution of blue and red bands, respectively. As shown in Figure 7, approximately 99% of the surface reflectance values were below 0.12 in the blue band, and almost 58% of the surface reflectance values were below 0.15 in the red band. As illustrated in Figure 2, when the surface reflectance in the blue and red bands is less than 0.12 or 0.15, respectively, changes in AOD can still cause changes in apparent reflectance. In addition, to minimize the error, the surface reflectance values of the blue and red bands were limited to less than 0.12 and 0.15, respectively.

Inversion Experiments
During the study period, the MODIS reflectance data covering the study area at a 500 m resolution (MOD02HKM) in July and August 2014 were selected using ground-based remote sensing measurements. In accordance with the selected MODIS data, inversion experiments were performed based on the obtained surface reflectance data. Figure 8 shows the spatial distribution of AOD over the study area retrieved from MODIS images at a 500 m resolution. These results demonstrate that the new algorithm can achieve a continuous AOD distribution, even in desert areas with high reflectance. The new algorithm describes the aerosol distribution and variability in detail and at a higher spatial resolution (500 m) than the current operational MODIS AOD product (MOD04) (10 or 3 km). As shown in Figure 8, although the retrieved overall AOD value is small, some large values appear at the edges of the clouds. The white areas in the map indicate that no data were retrieved in these areas, generally because they were covered by clouds or snow.

Inversion Experiments
During the study period, the MODIS reflectance data covering the study area at a 500 m resolution (MOD02HKM) in July and August 2014 were selected using ground-based remote sensing measurements. In accordance with the selected MODIS data, inversion experiments were performed based on the obtained surface reflectance data. Figure 8 shows the spatial distribution of AOD over the study area retrieved from MODIS images at a 500 m resolution. These results demonstrate that the new algorithm can achieve a continuous AOD distribution, even in desert areas with high reflectance. The new algorithm describes the aerosol distribution and variability in detail and at a higher spatial resolution (500 m) than the current operational MODIS AOD product (MOD04) (10 or 3 km). As shown in Figure 8, although the retrieved overall AOD value is small, some large values appear at the edges of the clouds. The white areas in the map indicate that no data were retrieved in these areas, generally because they were covered by clouds or snow.

Inversion Experiments
During the study period, the MODIS reflectance data covering the study area at a 500 m resolution (MOD02HKM) in July and August 2014 were selected using ground-based remote sensing measurements. In accordance with the selected MODIS data, inversion experiments were performed based on the obtained surface reflectance data. Figure 8 shows the spatial distribution of AOD over the study area retrieved from MODIS images at a 500 m resolution. These results demonstrate that the new algorithm can achieve a continuous AOD distribution, even in desert areas with high reflectance. The new algorithm describes the aerosol distribution and variability in detail and at a higher spatial resolution (500 m) than the current operational MODIS AOD product (MOD04) (10 or 3 km). As shown in Figure 8, although the retrieved overall AOD value is small, some large values appear at the edges of the clouds. The white areas in the map indicate that no data were retrieved in these areas, generally because they were covered by clouds or snow.

Validation
The AODs measured at different wavelengths were used to calculate the inversion results at 550 nm. The relationship between the AOD and wavelength can be expressed by the Ångström equation as follows [29]: where τ a (λ) is the AOD at λ, λ is the wavelength, α is the wavelength exponent closely related with the aerosols size, and β is the turbidity factor, which are written as: Because the solar spectrometer lacked a 550 nm channel, we chose the two channels closest to 550 nm and then calculated the AOD at 550 nm using Equation (6). Based on the available channels for the CE318 and Microtops II sun photometers, and the influence of water vapor on the radiative transfer [30]. As a result of previous studies [31][32][33][34], we selected the 500 nm and 670 nm channels and 500 nm and 870 nm channels, respectively.
To decrease the influences of atmospheric instability, such as pixel pollution due to the presence of thin clouds, ragged clouds, and cloud edges, the average inversion results for a 3 × 3 pixel area were selected as the final verification result. Forty-five pairs of valid data were obtained at the Dahuangshan and Wucaiwan sites. Figure 9 shows a comparison of the time series between the measured AODs and the retrieved results. The overall trends of these results exhibited good agreement, although the retrieved values were slightly smaller than the measured values. Figure 9A shows the Microtops II results at the Dahuangshan site. The maximum absolute difference between the inversion values was 0.12. Figure 9B shows the results of the CE318 at the Wucaiwan site, where the maximum absolute difference was 0.27. The largest differences occurred on 9 August, and the measured value was largest on that day. Based on the true color composite image, this error may be due to the impacts of clouds.

Validation
The AODs measured at different wavelengths were used to calculate the inversion results at 550 nm. The relationship between the AOD and wavelength can be expressed by the Ångström equation as follows [29]: where ( ) a   is the AOD at  ,  is the wavelength, α is the wavelength exponent closely related with the aerosols size, and β is the turbidity factor, which are written as: Because the solar spectrometer lacked a 550 nm channel, we chose the two channels closest to 550 nm and then calculated the AOD at 550 nm using Equation (6). Based on the available channels for the CE318 and Microtops II sun photometers, and the influence of water vapor on the radiative transfer [30]. As a result of previous studies [31][32][33][34], we selected the 500 nm and 670 nm channels and 500 nm and 870 nm channels, respectively.
To decrease the influences of atmospheric instability, such as pixel pollution due to the presence of thin clouds, ragged clouds, and cloud edges, the average inversion results for a 3 × 3 pixel area were selected as the final verification result. Forty-five pairs of valid data were obtained at the Dahuangshan and Wucaiwan sites. Figure 9 shows a comparison of the time series between the measured AODs and the retrieved results. The overall trends of these results exhibited good agreement, although the retrieved values were slightly smaller than the measured values. Figure 9A shows the Microtops II results at the Dahuangshan site. The maximum absolute difference between the inversion values was 0.12. Figure 9B shows the results of the CE318 at the Wucaiwan site, where the maximum absolute difference was 0.27. The largest differences occurred on 9 August, and the measured value was largest on that day. Based on the true color composite image, this error may be due to the impacts of clouds. The MODIS aerosol retrievals, except in coastal zones, lie within the retrieval error range of (0.05 0.2 ) a a        , which is the expected error (EE) defined by NASA [35]. In this paper, we use the NASA standard to verify the accuracy of the inversion results. Table 2 and Figure 10 show the statistics of the AOD values retrieved from MODIS images compared to the AOD measurements at the two sites. The validation results show that the AOD retrieved from the MODIS images by the algorithm outlined in this paper exhibits high consistency with the AODs measured by sun photometers (the R 2 is 0.84), and a total of 45 contrast points were obtained. For all points, The MODIS aerosol retrievals, except in coastal zones, lie within the retrieval error range of ∆τ a = ±(0.05 ± 0.2 × τ a ), which is the expected error (EE) defined by NASA [35]. In this paper, we use the NASA standard to verify the accuracy of the inversion results. Table 2 and Figure 10 show the statistics of the AOD values retrieved from MODIS images compared to the AOD measurements at the two sites. The validation results show that the AOD retrieved from the MODIS images by the algorithm outlined in this paper exhibits high consistency with the AODs measured by sun photometers (the R 2 is 0.84), and a total of 45 contrast points were obtained. For all points, the average absolute error was 0.60, the relative error was 21.85%, and 82.22% of the points were within the NASA EE range. Because the Dahuangshan site is located near the highly vegetated Tianshan Mountains, the surface reflectance is low at the site. Additionally, the surface reflectance can be relatively accurately determined, and the inversion result is satisfactory. the average absolute error was 0.60, the relative error was 21.85%, and 82.22% of the points were within the NASA EE range. Because the Dahuangshan site is located near the highly vegetated Tianshan Mountains, the surface reflectance is low at the site. Additionally, the surface reflectance can be relatively accurately determined, and the inversion result is satisfactory.  The global aerosol data provided by MOD04 products achieve a certain level of accuracy. Most validation studies have concluded that, in general, the MODIS algorithm used to retrieve aerosol products provides results comparable to sun photometer-based products and that an expected error envelope could be defined that contained at least a 67% (approximately one standard deviation) agreement over vegetated areas [16]. The MODIS algorithm has been widely applied for climate modeling, air quality testing, etc. [36]; however, a similar result was found in northern desert areas, with the ratio of MODIS to observations ranging from 15% to 55%. Approximately 7% to 39% of MODIS products fall within NASA's EE range, and most of the inversion results are higher than the ground-measured values [37]. Figure 11 shows the distribution of comparable MODIS MOD04 products during the study campaign, illustrating that this product was suitable only for areas covered by dense vegetation in Tianshan, and it had a substantially altered spatial distribution and lower resolution (10 km). Figure 12 shows that the scatter plots of the AOD retrieved using the proposed method are validated by the ground-based sun photometer AOD measurements. In Figure 12, the black dashed lines represent NASA's EE lines, the black solid lines are the 1:1 lines, and the red solid lines represent the regression lines. Figure 12A,B show the results retrieved using our proposed algorithm and the MOD04 AOD results, which have R 2 values of 0.84 and 0.10, respectively. The AOD results obtained using the new algorithm better fit the ground-measured values, with more points lying within NASA's EE range. The global aerosol data provided by MOD04 products achieve a certain level of accuracy. Most validation studies have concluded that, in general, the MODIS algorithm used to retrieve aerosol products provides results comparable to sun photometer-based products and that an expected error envelope could be defined that contained at least a 67% (approximately one standard deviation) agreement over vegetated areas [16]. The MODIS algorithm has been widely applied for climate modeling, air quality testing, etc. [36]; however, a similar result was found in northern desert areas, with the ratio of MODIS to observations ranging from 15% to 55%. Approximately 7% to 39% of MODIS products fall within NASA's EE range, and most of the inversion results are higher than the ground-measured values [37]. Figure 11 shows the distribution of comparable MODIS MOD04 products during the study campaign, illustrating that this product was suitable only for areas covered by dense vegetation in Tianshan, and it had a substantially altered spatial distribution and lower resolution (10 km). Figure 12 shows that the scatter plots of the AOD retrieved using the proposed method are validated by the ground-based sun photometer AOD measurements. In Figure 12, the black dashed lines represent NASA's EE lines, the black solid lines are the 1:1 lines, and the red solid lines represent the regression lines. Figure 12A,B show the results retrieved using our proposed algorithm and the MOD04 AOD results, which have R 2 values of 0.84 and 0.10, respectively. The AOD results obtained using the new algorithm better fit the ground-measured values, with more points lying within NASA's EE range.

Conclusions and Discussion
Based on the land surface characteristics in arid areas, this paper proposes a method for estimating land surface reflectance (LSR) and retrieving aerosol information over arid areas using long-term series of MODIS images. The clearest image was selected from a series of continuous MODIS data with matching spatial characteristics and at adjacent times. Atmospheric correction was

Conclusions and Discussion
Based on the land surface characteristics in arid areas, this paper proposes a method for estimating land surface reflectance (LSR) and retrieving aerosol information over arid areas using long-term series of MODIS images. The clearest image was selected from a series of continuous MODIS data with matching spatial characteristics and at adjacent times. Atmospheric correction was

Conclusions and Discussion
Based on the land surface characteristics in arid areas, this paper proposes a method for estimating land surface reflectance (LSR) and retrieving aerosol information over arid areas using long-term series of MODIS images. The clearest image was selected from a series of continuous MODIS data with matching spatial characteristics and at adjacent times. Atmospheric correction was performed using ground-measured atmospheric parameters and then used to estimate the LSR for AOD retrieval from data collected at adjacent times. Validation was performed by comparing the results from two different sun photometers (CE318 and Microtops II) located at two different sites (Dahuangshan and Wucaiwan) in arid regions; both yielded comparable results. A validation analysis using ground-measured data and the MOD04 aerosol product allowed us to draw the following conclusions: (1) Based on the 6S model simulation analysis, the sensitivity of the apparent reflectance and the AOD under different conditions, the path radiance of aerosols, and the radiance of the land surface based on the satellite signal can be identified when the surface reflectance values of blue and red wavelengths are less than 0.12 and 0.15, respectively, except when using the urban/industrial aerosol model. To obtain AOD, surface reflectance must be accurately determined; (2) The AOD results retrieved using the new algorithm represent a more accurate estimation of aerosol loading over the arid study region and can describe the aerosol distribution and variability in greater detail at a high spatial resolution (500 m); (3) This validation shows that the retrieved AODs are highly correlated with the AODs from the ground stations. The R 2 reached 0.84, while the relative error was 21.85%, and 82.2% of the retrieved AODs fell within the EE range defined by NASA.
Over arid coverage areas, the proposed inversion algorithm based on sequences of images with adjacent times can compensate for some of the shortcomings of the DDV method. However, the following three limitations remain: (1) surface directional reflectance still induces some errors, although the errors are generally less than 0.035; (2) the algorithm is not suitable for areas with extremely large reflectance, such as ice-covered areas; and (3) the requirements of ground-measured data and the large area of automatic retrieval are restricting.