Retrieval of High-Resolution Atmospheric Particulate Matter Concentrations from Satellite-Based Aerosol Optical Thickness over the Pearl River Delta Area , China

Satellite remote sensing offers an effective approach to estimate indicators of air quality on a large scale. It is critically significant for air quality monitoring in areas experiencing rapid urbanization and consequently severe air pollution, like the Pearl River Delta (PRD) in China. This paper starts with examining ground observations of particulate matter (PM) and the relationship between PM10 (particles smaller than 10 μm) and aerosol optical thickness (AOT) by analyzing observations on the sampling sites in the PRD. A linear regression (R2 = 0.51) is carried out using MODIS-derived 500 m-resolution AOT and PM10 concentration from monitoring stations. Data of atmospheric boundary layer (ABL) height and relative humidity are used to make vertical and humidity corrections on AOT. Results after correction show higher correlations (R2 = 0.55) between extinction coefficient and PM10. However, coarse spatial resolution of meteorological data affects the smoothness of retrieved maps, which suggests high-resolution and accurate meteorological data are critical to increase retrieval accuracy of PM. Finally, the model provides the spatial distribution maps of instantaneous and yearly average PM10 over the PRD. It is proved that observed PM10 is more relevant to yearly mean AOT than instantaneous values. OPEN ACCESS Remote Sens. 2015, 7 7915


Introduction
Particulate matters (PM), or aerosols, not only have direct and indirect radiative forcing effects on climate systems [1,2], but also reduce visibility and induce respiratory diseases thus affecting air quality and human health [3].The strong and consistent links between ambient PM exposure and adverse health effects have been shown in epidemiological and toxicological studies around the world [4][5][6][7].Therefore, monitoring air quality at high spatial-temporal resolution is of great significance, especially for rapidly growing megacities which are facing severe air pollution owing to industrial development and population expansion.
Current ground measurement networks, though well-calibrated, are spatially-limited and inadequate to evaluate time-space dynamics of air pollution and effects on human health.Advances in satellite remote sensing have provided a cost-effective approach to frequently monitor the spatial distribution and the regional transport of aerosol particles on a large scale [8][9][10][11][12][13][14], though the discontinuous satellite observations could hardly reflect the diurnal variations of particle matters.Aerosol optical thickness (AOT), derived from satellite observations, is a dimensionless measure of aerosol abundance and atmospheric extinction of radiance over the entire air column.Previous researches have revealed a strong positive relationship between satellite-based measurements of AOT and ground-based measurements of particulate matters smaller than 2.5 μm and 10 μm in aerodynamic diameters (PM2.5 and PM10, respectively) [15][16][17][18][19][20][21][22][23][24][25].A summary of linear AOT/PM relationships, intercepts and correlation coefficients in several previous studies is shown in Table 1.IDEA (Infusion satellite Data into Environmental air quality application), as a NASA-EPA-NOAA partnership, has produced the daily surface PM estimation across the United States from MODIS AOT products through predefined regression relation.Its time ranges from 1 January 2008 to the present, with a coarse spatial resolution of 10 km.
Located in south-central Guangdong, China (about 21.5°N-24°N and 112°E-115.5°E),the Pearl River Delta (PRD) (Figure 1) is an area that boasts the fastest economic growth and population increase in the world.It also has become an area of China which is facing the severest air pollution because of rapid industrial development, vegetation reduction, and heavy traffic pressure.Deaths and health effects resulting from air pollution in the PRD had been calculated from an economic perspective by scholars at Peking University.Economic losses caused by air pollution were estimated to be 292 billion RMB in 2006, accounting for 1.35% of the regional Gross Domestic Product (GDP) [26].As a consequence, the increase in hazy days and deteriorating air quality, greatly impairing the health of the population and the visibility over the region, has drawn widespread attention from the Chinese government and the public [27][28][29][30].However, there is still a lack of AOT-based estimations for PM concentrations in this area currently.By combining ground measurements and remote sensing observations, estimations of AOT-based PM concentrations, especially with high spatial resolution data, could offer spatially continuous mapping with high accuracy, which is useful for regional air quality monitoring.In this study, we developed correlative linear models between AOT and PM10 in the Pearl River Delta region, based upon ground-based observations from 2010 to 2011.We also explored the potential of retrieving the spatial distribution of particulate matters both daily and yearly from MODIS-derived AOT using these models, with an aim to provide more intuitive and reliable images of air quality in the area and reasonably illustrate the usefulness of satellite-derived AOT as a proxy for particulate matters.

Station Monitoring Data of PM10
Hourly air quality measurements from 15 monitoring stations Figure 1 from January 2008 to December 2008 were provided by the Department of Environmental Protection of Guangdong Province, China.Locations and the underlying surface types of the stations are shown in Table 2 and quoted from a report of monitoring results in 2008: Pearl River Delta regional air quality monitoring network.The measurements included concentrations of PM10 and major gaseous pollutants such as surface-level O3, SO2 and NO2.Also, the government has fully implemented the QA/QC program to ensure that the air quality data from the monitoring stations are highly accurate and reliable, and the control limit set for PM10 are ±10%.Nevertheless, the sparse distribution of stations cannot reliably reflect the spatial distribution of atmospheric particulate matter over the Pearl River Delta region.

Ground-Based Sampling Observations of AOT and PM
AOT ground observations were obtained using a MICROTOP II sun photometer Manufactured by U.S Solarlight Company.The MICROTOP II is a hand-held sun photometer for measuring aerosol optical thickness and direct solar irradiance at the following five discrete wavelengths: 380 nm, 500 nm, 870 nm, 936 nm and 1020 nm.The water vapor column can also be measured at the following three wavelengths: 870 nm, 936 nm and 1020 nm.Features of MICROTOP II include high accuracy, ease of use, portability and instantaneous results.
A DustTrak DRX aerosol particulate monitor was used to measure the particulate concentrations in a wide variety of environments.It provided reliable exposure assessment by simultaneously measuring size-segregated mass fraction concentrations corresponding to PM1, PM2.5, PM10 and Total PM size fractions.The one-minute resolution of DustTrak DRX offers us an insight into the occurrence of high concentrations of PM, and the clear duration of the dust plumes.The comparison among eight PM10 measurement instruments shows that multiple size ranges in the DustTrak DRX could identify the spatial non-uniformity for the sources [31].
For the limitation of monitoring devices and manpower, we chose several sampling sites, including urban areas, suburbs and their borders in the three cities of the Pearl River Delta, to avoid the uncertainty caused by observations on a single underlying surface Table 3.The total number of sample points is 305, with 162 points in the urban areas and 143 points in the suburbs.Rapid urbanization process makes the difference between the locations less apparent, thus the monitoring results reflect the air quality of urban actually.The sampling time was from 10:30 a.m. to 3:00 p.m. between October 2010 and November 2011, to be consistent with the pass time of MODIS, about 11 a.m. for Terra and 2 p.m. for Aqua.

Satellite Data
Data used in this study over the PRD from a Moderate-resolution imaging spectroradiometer (MODIS) (onboard the Terra and Aqua satellite) were downloaded from NASA LAADS (Level 1 and Atmosphere Archive and Distribution System).They are calibrated L1B data with 500 m spatial resolution including visible, near-infrared, and mid-infrared bands.The complete set of cloud-free remote sensing images during 2008 was acquired, together with monitoring data of air quality from 15 stations in 2008 provided by the Department of Environmental Protection of Guangdong Province, China.

Atmospheric Boundary Layer (ABL) Height and Relative Humidity (RH)
The National Centers for Environmental Prediction (NCEP) Climate Forest System Reanalysis (CFSR) provide global atmospheric analyzed products, such as boundary layer height and relative humidity, with a 0.5 degree horizontal resolution at 6-hourly intervals, initialized four times per day (0:00, 6:00, 12:00 and 18:00 UTC) [32] The products are proved to be accurate by using the 5-day forecast scores as a measure of the accuracy of initial states [33].We used data from 2, 3 and 5 January 2008 to perform height and humidity corrections on AOT.

Ground Observations of Particulate Matter
The concentrations of PM10, PM2.5 and PM1 are expressed in terms of milliequivalents of particles smaller than 10 μm, 2.5 μm and 1 μm in diameter per cubic meter of air in units of µg/m 3 .Ground-based observations of PM10, PM2.5 and PM1 at all the sampling sites from November 2010 to November 2011 were used to analyze the size distribution of particles.The percentage of PM1 from PM2.5 and PM10 calculated from all the sample points range from 90% to 100% and from 60% to 90%, respectively.Most of the atmospheric particles over the sampling sites are small-size respirable particles, consistent with research results by Wu [27,34,35].He pointed out that the ratio of PM2.5 to PM10 over PRD has obviously increased in the past two decades, reaching 60% to 80%, which is higher in dry seasons than in rainy seasons.Our results are slightly higher, perhaps caused by the difference between measuring instruments.The aerosol monitor used in our study is a continuous real-time 90° light-scattering laser photometer.Aerosol particles are absorbed into the optics chamber through the built-in air pump, and the concentration of the particles is measured by the amount of light scattering.The measurements of fine-mode particles are less accurate than those from β-ray decay method and tapered element oscillating microbalance (TEOM) method.

Station Monitoring Data of Particulate Matters in PRD
Measurements of PM10 from the 15 stations in the PRD at 11 a.m.(the passing time of Terra satellite) from 1 to 5 January are listed in Table 4.The PM10 concentration on most of the sites on 2 January increased especially at stations in Nansha (Guangzhou) and Foshan compared with those on 1 January.Dongguan, Foshan and Nansha became areas with peak value on 3 January and the value from Huijingcheng in Foshan was the highest, followed by Donghu in Jiangmen, Haogang School in Dongguan and Chenzhong in Zhaoqing on 4 January.In general, station measurements in Guangzhou, Conghua, Huizhou, Shenzhen and Hongkong did not change much during the period.Our field survey reveals that stations with small variations are mostly located in areas with high vegetation cover.This apparently enhances subsidence of atmospheric particles.
The hourly PM10 limit concentration in the Pearl River Delta is set to be 150 µg/m 3 , according to the ambient air quality standards issued by China's Ministry of Environmental Protection and China's General Administration of Quality Supervision in 2012.Inspection of Table 4 shows that on 1 and 2 January PM10 of all sites were below this standard.However, on days of severe air pollution such as 3, 4 and 5 January, measurements were far above the standard threshold 150 µg/m 3 at Jinjuzui and Huijincheng in Foshan, Haogang school in Dongguan, Donghu in Jiangmen and Wanqinsha in Guangzhou.Concentrations of PM10 from several stations, exceeding the standard (150 µg/m 3 ), are marked in bold format.

Correlation Test between Particulate Concentration and AOT
We developed linear relationships between PM10 and aerosol optical thickness (AOT) at 500 nm by using 305 datasets of ground-based observation from our instruments between 18 October and 29 November 2011 (Figure 2), with correlation coefficients of 0.67.In winter, the hazy weather occurred more frequently and visibility dramatically deteriorated over the Pearl River Delta region [35].The strong linear correlation between AOT and PM reveals that it is reasonable and feasible to retrieve the spatial distribution of PM from AOT. Factors affecting the accuracy of the models include: (a) AOT are column measurements whereas PM are recorded 5 feet above the ground; (b) the average particle sizes of aerosols increases with increasing relative humidity [36,37].Therefore, it is concluded that satellite-retrieved aerosol maps provide the possibility and capability of monitoring and characterizing spatio-temporal distribution of the surface PM with a reasonable degree of accuracy.
Figure 2. Linear regression between AOT at 500 nm and ground measurements PM10.

Retrieval of Atmospheric Particulate Concentrations
A unary linear regression is used to estimate the spatial distribution of particulate matter over the PRD, as expressed in Equation ( 1): where Y is the predicted variable, X are the independent variables of the model, a is a constant and b is the coefficient of the variables.In our model, X represents AOT values with high-resolution, and Y is predicted PM concentrations.As Figure 2   Retrieval of 500 m-resolution AOT in this study is based on an improved dark dense vegetation (DDV) method, proposed by our previous study [38].Aerosol products of satellite remote-sensing provide an effective way to study air pollution, but have a conflict between resolution and signal-to-noise ratio (SNR).Our method, on a basis of MOD04-C005 algorithm, hopes to obtain Aerosol Optical Thickness (AOT) with a higher spatial resolution in the case of moderate decreases in SNR [39].The theory of the method is that in the case of low surface reflectance, the 2.1-μm channel is transparent to most aerosol types so that its apparent reflectance can be considered to be equal to the surface reflectance [40,41].The relations among the surface reflectance in the blue, red and SWIR band are not a fixed ratio but depend on the scattering angle and NDVISWIR derived from reflectance at 1.24 and 2.12 μm bands, taking into account the angular variability and surface type, since the Earth's surface is not Lambertian and some surface types exhibit bidirectional reflectance functions (BRDF).Besides, before AOT retrieval, the MODIS cloud mask products (MOD35_L2) were applied to remove the cloud on MODIS L1B images.For the remaining potentially clear pixels, the selection of dark pixels by discarding brightest 50% and darkest 20% of reflectance at 0.66 μm, also reduce the cloud and surface contamination.Thus, the AOT can be retrieved from the surface reflectance due to the correlation between surface reflection in the blue (0.49 μm), red (0.66 μm), and 2.12 μm bands.Changing the movement pattern of the retrieval window, selecting a more suitable aerosol type, and storing the look-up table as a four-dimensional array, enhance the spatial resolution of AOT considerably relative to MODIS AOT products.The selection of effective dark pixels is conducted in a window of 20 × 20 pixels (10 km × 10 km), and the calculation of the surface reflectance at 0.47 and 0.66 μm are functions of the surface reflectance at 2.13 μm, NDVI and the scattering angle.We also adopted the aerosol type suitable for the study area proposed by Li [21], rather than using the urban or continental aerosols in the 6S model.All of these factors would reduce the errors of AOT retrieval on the bright surfaces, especially in urban areas.The AOT maps by our algorithm are much smoother and show more extensive value ranges with fewer null values.We also validated the retrieval results with ground measurements and 10 km-resolution MODIS AOT products, both showing high precision.We will further validate our retrieval results of aerosol optical thickness in specific seasons and regions once we get ground observations of AOT and PM on a long time scale.

Retrieval of the Instantaneous Particulate Distribution
Data involved in the retrieval include MODIS L1B data onboard the Terra satellite covering the whole Pearl River Delta, and PM10 measurements at 10:00 to 11:00 a.m. at the 15 stations on 2, 3 and 5 January 2008, together with atmospheric boundary height and relative humidity were downloaded from NCEP CFSR on the same days.The MODIS instrument passes the entire Earth twice per day, with Terra in the morning and Aqua in the afternoon, while the PM data at the stations were recorded hourly.To reduce the effects of the difference in temporal resolution on the relationship between PM10 and AOT, PM data recorded from 10:00 to 11:00 at the 15 point locations were used to match the Terra data in our study.After deriving the AOT from the MODIS images, the relationship between average AOT values within 500 m radius of the monitoring station and PM10 from all sites was established as shown in Figure 3.The accuracy of the relationship is affected by the imperfect match between the spatial-temporal resolutions of AOT-PM data, as the AOT was estimated at about 10:30 a.m. with a spatial resolution of 500 m, and the null values were excluded from the calculation.The PM data, on the other hand were recorded from 10:00 to 11:00 at the 15 point locations.In addition, the relative humidity and atmospheric boundary layer height have impacts on the correlation between column AOT and ground PM10 concentrations.The spatial distribution maps of PM10 on 2, 3 and 5 January using the linear model in Figure 3 are shown in Figure 4.The difference between the slope values of the models in Figures 2  and 3 is possibly caused by the measuring instruments, and the PM concentrations from ground observations are a little higher than that from the monitoring stations.Variations in local meteorological conditions and occurrence of multiple aerosol layers play important roles in the relationship between AOT and PM10.To estimate the PM concentrations directly from AOT may be somewhat tenuous, and the vertical distribution of aerosols and the relative humidity, related to atmospheric profiles, ambient condition, the size distributions and chemical compositions of aerosols, are especially important.As AOT reflects the aerosol optical properties of the total column whereas particle matter concentrations are usually recorded at 5 feet above the ground, the relation between them is greatly influenced by the vertical distribution of aerosols.As we know, AOT is the integral of the extinction coefficient ka at all altitudes along the vertical orientation, and is expressed as: where λ is the total atmospheric optical thickness, and λ, is the extinction coefficient at the altitude of z and the wavelength of λ.Additionally, the vertical distribution of λ, could be described as the negative exponent form as: ,0 ( , ) ( ) exp( / ) where , λ is the extinction coefficient at the wavelength of λ near the surface, and HA is the scale height of aerosol, approximately represented by the atmospheric boundary layer height (ABL) [42,43].
Combining Equations ( 2) and ( 3), we calculate λ as: On the other hand, the hygroscopic growth of particles has effects on the refraction and extinction index, as well as other optical properties of aerosols [44].The correlation between extinction coefficient ka and PM10 is also affected by the chemical components of particles and relative humidity of the air, since the PM10 concentrations measured by the instruments are almost the dry mass of the particles with aerodynamic diameter less than 10 μm.Thus, f (RH) is used to define the hygroscopic growing factor as [36,45]: where g is an empirical fie coefficient and set as 1 in our study.The dry , λ is obtained from Hence, we could calculate the aerosol extinction coefficient in the dry air near the surface , λ from Equation (1) to Equation ( 5), and express it as Wang et al. [24] once described the relationship between ka and PM concentration on the basis of the Mie theory as where Qext is the size distribution integrated extinction efficiency, reff is the effective radius, being approximately constant, and PMx is the mass concentration of PM.Therefore, it is possible to develop a linear correlation between , λ and PM10 concentrations, like Furthermore, the vertical and relative humidity corrections on AOT could enhance the relationship and increase the robustness of the estimate.Spatial distributions of ABL height and relative humidity on 2 and 3 January are displayed in Figure 5.The comparison of measurements at the weather stations and forecast of relative humidity is shown in Figure 6, with a correlation coefficient of 0.816.The deviation of the correlation is caused by the different temporal resolutions of these two datasets.Data from NCEP CFSR was at 6-hourly intervals while the station observations are hourly averaged.ABL height and relative humidity data at 6:00 UTC (14:00 pm local time), which are the closest datasets to the overpassing time of satellite, were applied to the vertical and humidity correction on MODIS Terra images to reduce the errors from temporal scale.The correlation between average corrected , λ and PM10 is shown in Figure 7.Comparison between Figures 3 and 7 indicates that the correlation (R 2 = 0.55) is slightly higher after this correction but not significant, mainly because of the low spatial resolution of meteorological data.The final retrieved PM10 images over the Pearl River Delta on 2, 3 and 5 January are displayed in Figure 8.The coarse spatial resolution of meteorological data leads to apparent grids in the images in Figure 8.We therefore averaged the values of ABL height and relative humidity from NECP CFSR over the PRD to smooth and refine the maps.To do so, we took HA and RH in Equation ( 2) to be single valued rather than maps with 0.5° resolution.The relationship between corrected AOT and PM10 is shown in Figure 9 with a coefficient of 0.531.It indicates that the correlation can be enhanced by vertical and humidity corrections, but the degree of improvement depends on the accuracy of the meteorological data.The retrieved PM10 images over the Pearl River Delta on 2, 3 and 5 January are displayed in Figure 10.

Retrieval of Yearly Average PM10 Distributions
We used all cloud-free MODIS data in 2008 to derive the corresponding AOT using the dark dense vegetation method.Averaging the results leads to an estimate of the yearly aerosol optical thickness over the PRD Figure 11.However, in summer the common occurrence of clouds makes AOT rarely observable, thus only a few MODIS images were studied.The average results would be a little lower since AOT values in summer are higher than these in winter [46].Additionally we estimated the PM concentrations in non-site areas from our regression model.According to the station measurements provided by Guangdong Environmental Protection Bureau, the days of PM10 exceeding the standard threshold 150 µg/m 3  in Liyuan, 35 in Zimaling Park, 2 in Tamen, 1 in Quanwan and 0 in Tangjia, respectively.Linear regression of these data is shown in Figure 12 with a correlation coefficient of 0.582.The slope coefficient in the linear regression model of yearly average AOT and PM10 is relatively low, probably because of the average calculation on AOT and fewer sample points.The correlation between yearly mean AOT and observed PM10 is slightly higher than that between the instantaneous values, mainly because averaging tends to eliminate outliers.The instantaneous AOT is influenced by many factors that show frequent fluctuations, such as temperature, monsoon, and precipitation, as well as human activities, such as fossil fuel burning.Retrieved AOT maps cannot reflect the true air quality when there appears to be a low or high value anomaly of these factors.Occurrence of clouds also makes aerosols sometimes rarely observable, accompanied by null values of AOT.Outliers and missing data may introduce bias or affect the representativeness of the results.The yearly average AOT can remove the occasional effects and eliminate the outliers of AOT by replacing them with the average values of the pixel.Thus, the larger the time scales, the higher the correlation between aerosol optical thickness and PM10 concentration.The yearly average PM10 is then retrieved through the regression equation (Figure 13).12, by putting in the yearly average AOT values at each monitoring station.The deviations between the actually-measured and predicted PM10 concentrations were calculated to indicate the stability of the model.The root mean square error (RMSE) between station measured PM10 and the predicted PM10 was calculated from the deviations to be ±10 µg/m 3 , indicative of the reliability of the retrieved results while the relative errors of the predicted results are from 1.75% to 27.12%, showing the limitations of the applicability of the model.Air quality at most sites is poor according to the standards of both yearly average PM10 (70 µg/m 3 ) laid down in the ambient air quality standards.Figure 13 show that areas with the worst air pollution are located in Foshan city, at the border of Shunde and Zhongshan, west of Gongguan, and northwest of Shenzhen.

Conclusions
Analysis on the size distribution of particles indicates that most of the atmospheric particulate matters in Pearl River Delta region are inhaled particles of small sizes.Simultaneous observations of AOT and PM10 are proved to be strongly correlated; R 2 being 0.670.It forms the basis for successful linear regression between AOT and PM.
A linear regression (R 2 = 0.51) is carried out using MODIS-derived 500 m-resolution aerosol optical thickness and PM10 concentration from monitoring stations on 2, 3 and 5 January 2008.However, the relationship between AOT and PM over PRD is not necessarily suitable for other regions because of regional variations in aerosols type, sources, and air pollution affecting the relationship between AOT and PM.In addition the relationship may vary over time, for seasonal differences in climate would affect the formation and movement of atmospheric particles.The model in our study successfully demonstrates the possibility and feasibility of retrieval particulate matters from satellite-derived AOT, especially high-resolution AOT maps, which makes air quality monitoring in smaller regions more reliable.The correlation between AOT and PM10 also increased a little after vertical and relative humidity correction.Our further work aims at using more accurate meteorological data to increase retrieval accuracy for air pollution assessment.
Retrieval of atmospheric particulate concentrations shows that yearly mean AOT and observed PM10 correlate better than instantaneous values, and the root mean square error (RMSE) of our predictive model was 0.01.The yearly average particle matter concentration in 2008 shows the four most polluted regions to be: Foshan city, the border between Shunde and Zhongshan, the area to the west of Dongguan, and the northwest of Shenzhen.
However, there still exist some uncertainties in the linear model.The relationship between AOT and PM varies with time and region, and the chemical reactions in the air are more complicated than the linear regression models reflect [47].Assimilating the satellite products and meteorological data into numerical models could be a good way to retrieve the spatial distributions of PM.Also, the relation between AOT and PM are influenced by the lack of accurate information on particle composition, size distribution and vertical profile.It is also affected by the accuracy and spatio-temporal resolution of the input data.The variations of water vapor could lead to the swelling or condensation of particles, and consequently the change of the microstructure and composition of the particles, which affects the relationship between AOT and PM.Vertical profile of aerosols could also introduce errors in the relationship, because of the significant contribution of aerosols above the boundary layer to extinction, especially when weak atmospheric convection happens [48].Therefore, information on the aerosol vertical profiles from ground-based lidars and space-borne lidars, such as the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), could be applied into the retrieval of PM from satellite data.The accuracy of the input data, like the satellite-derived AOT, the instruments for ground measurements and the meteorological data for corrections, would also affect the linear model between AOT and PM.Furthermore, the satellite images covering the local region are obtained twice a day, and have a relatively coarse spatial resolution, while the station data are measured on the ground all day.The imperfect match of spatio-temporal resolution between satellite data and surface observations brings errors to the model.The seasonal variation of the relationship is another issue that should be taken into account.The seasonal difference in the meteorological conditions, such as temperature, precipitation and relative humidity, and human activities could result in the change of the relationship.

Figure 1 .
Figure 1.Location of 11 cities and the 15 air quality monitoring stations in the Pearl River Delta region in Guangdong Province, China.The blue spots in the map represent the sampling sites selected to conduct our ground observations of PM and AOT.
shows, PM10 and AOT are highly linear-correlated, thus an observed set of values of station-measured PM10 concentrations and satellite-derived AOT values at the same place in 2008 could be used to develop the linear model in Figure 3.Then, through the established regression model, the unmeasured PM concentrations in non-station areas can be retrieved from 500 m-resolution aerosol optical thickness, for the satellite-derived AOT maps are spatially continuous, and show values both at station and non-station points.

Figure 3 .
Figure 3. Correlation diagram between average aerosol optical thickness within 500 m radius of the monitoring station and station-measured PM10 concentrations.

Figure 5 .Figure 6 .
Figure 5. Spatial distribution of atmospheric boundary layer (ABL) height (unit: meter) on the left and relative humidity (unit: %) on the right in PRD at 6:00 UTC on 2, 3 and 5 January 2008, provided by NCEP CFSR, with a resolution of 0.5 degree.

Figure 7 .
Figure 7. Correlation between extinction coefficient and surface-level PM10 after vertical and humidity correction using meteorologic data form NCEP CFSR.

Figure 9 .
Figure 9. Correlation between extinction coefficient corrected by average value of atmospheric boundary layer height and relative humidity, and the surface-level PM10.

Figure 11 .
Figure 11.A 500 m-resolution yearly average aerosol optical thickness in 2008 over the Pearl River Delta region retrieved from all cloud-free MODIS L1B data by using an improved Dark Dense Vegetation method.

Figure 12 .
Figure 12.Regression analysis between the yearly averaged retrieved AOT and ground-measured PM10 of the 15 monitoring stations in the PRD, 2008.

Table 1 .
Literature survey linear AOT-PM relationships, intercepts and correlation coefficients.

Table 2 .
Locations and types of the monitoring stations in the Pearl River Delta.

Table 3 .
Location of AOT and PM sampling sites in the Pearl River Delta.
Notes: N is the number of sampling dataset.

Table 5
lists the yearly average values of AOT, measured PM10 and regressive PM10 of the 15 monitoring stations in PRD, 2008.The predicted PM10 was derived from the model shown in Figure

Table 5 .
Yearly average values of AOT, measured PM10 and Predicted PM10 of the 15 monitoring stations in the PRD, 2008.