Global Estimates of Ambient Fine Particulate Matter Concentrations from Satellite-Based Aerosol Optical Depth: Development and Application
Aaron van Donkelaar1, Randall V. Martin1,2, Michael Brauer3, Ralph Kahn4, Robert Levy4, Carolyn Verduzco1, Paul J. Villeneuve5,6
1Department of Physics and Atmospheric Science, Dalhousie University, Halifax, Nova Scotia, Canada; 2Harvard-Smithsonian Center for Astrophysics, Cambridge, Massachusetts, USA; 3School of Environmental Health, University of British Columbia, British Columbia, Canada; 4NASA Goddard Space Flight Center, Greenbelt, Maryland, USA; 5Dalla Lana School of Public Health, University of Toronto, Toronto, Ontario, Canada; 6Population Studies Division, Health Canada, Ottawa, Ontario, Canada
Environ Health Perspect 118:847-855 (2010). http://dx.doi.org/10.1289/ehp.0901623 [online 16 March 2010]
Background: Epidemiologic and health impact studies of fine particulate matter with diameter < 2.5 µm (PM2.5) are limited by the lack of monitoring data, especially in developing countries. Satellite observations offer valuable global information about PM2.5 concentrations.
Objective: In this study, we developed a technique for estimating surface PM2.5 concentrations from satellite observations.
Methods: We mapped global ground-level PM2.5 concentrations using total column aerosol optical depth (AOD) from the MODIS (Moderate Resolution Imaging Spectroradiometer) and MISR (Multiangle Imaging Spectroradiometer) satellite instruments and coincident aerosol vertical profiles from the GEOS-Chem global chemical transport model.
Results: We determined that global estimates of long-term average (1 January 2001 to 31 December 2006) PM2.5 concentrations at approximately 10 km × 10 km resolution indicate a global population-weighted geometric mean PM2.5 concentration of 20 µg/m3. The World Health Organization Air Quality PM2.5 Interim Target-1 (35 µg/m3 annual average) is exceeded over central and eastern Asia for 38% and for 50% of the population, respectively. Annual mean PM2.5 concentrations exceed 80 µg/m3 over eastern China. Our evaluation of the satellite-derived estimate with ground-based in situ measurements indicates significant spatial agreement with North American measurements (r = 0.77; slope = 1.07; n = 1057) and with noncoincident measurements elsewhere (r = 0.83; slope = 0.86; n = 244). The 1 SD of uncertainty in the satellite-derived PM2.5 is 25%, which is inferred from the AOD retrieval and from aerosol vertical profile errors and sampling. The global population-weighted mean uncertainty is 6.7 µg/m3.
Conclusions: Satellite-derived total-column AOD, when combined with a chemical transport model, provides estimates of global long-term average PM2.5 concentrations.
Key words: aerosol, aerosol optical depth, AOD, particulate matter, PM2.5
Address correspondence to A. van Donkelaar, Department of Physics and Atmospheric Science, Dalhousie University, 6300 Coburg Rd., Halifax, NS, Canada B3H 3J5. Telephone: (902) 494-1820. Fax: (902) 494-5191. E-mail: Aaron.van.Donkelaar@dal.ca
This study was funded by contracts from Health Canada (4500171909 and 4500220294). A.v.D. was supported by graduate fellowships from the Natural Sciences and Engineering Research of Canada and the Killam Trust.
The authors declare they have no actual or potential competing financial interests.
Received 26 October 2009; accepted 25 February 2010; online 16 March 2010.
Chronic exposure to airborne fine particulate matter with diameter < 2.5 µm (PM2.5) is associated with adverse human health impacts including morbidity and mortality (e.g., Dockery et al. 1993; McDonnell et al. 2000; Pope et al. 2009). Several national environmental agencies in North America and Europe monitor PM2.5 concentrations at numerous sites throughout their jurisdictions, but even these relatively dense networks have limited geographic coverage. Few long-term measurement sites exist elsewhere in the world, particularly in rapidly developing countries where concentrations and estimated health impacts are greatest (Cohen et al. 2004). Point measurements collected at monitoring sites are not necessarily representative of regional concentration, and regional variability is difficult to assess from point measurements alone. In recent years, application of satellite observation to surface air quality has advanced considerably (Hoff and Christopher 2009; Martin 2008). In fact, global aerosol observations from satellite remote sensing could substantially improve estimates of population exposure to PM2.5.
Since the mid 2000s, the MODIS (Moderate Resolution Imaging Spectroradiometer) and MISR (Multiangle Imaging Spectroradiometer) instruments onboard the National Aeronautics and Space Administration’s (NASA) Terra satellite has provided global observations of aerosol optical depth (AOD), a measure of light extinction by aerosol in the atmospheric column above the earth’s surface. Terra’s sun-synchronous orbit encircles the earth approximately 15 times each day, with each pass crossing the equator at approximately 1030 hours local solar time. Observations of AOD from Terra provide daily insight into the global distribution of column-integrated aerosol. However, the applicability of AOD to surface air quality depends on several factors, including the vertical structure, composition, size distribution, and water content of atmospheric aerosol.
Many studies have investigated the relationship between total-column AOD and surface PM2.5 measurements. Most have developed simple empirical relationships between these two variables (e.g., Engel-Cox et al. 2004a; Wang and Christopher 2003); more recent investigations often have used local meteorological information to better relate AOD and PM2.5 (e.g., Koelemeijer et al. 2006; Liu et al. 2005) or to filter the AOD (e.g., Gupta et al. 2006). Some studies have employed light detection and ranging (LIDAR) instruments to capture the vertical aerosol distribution at specific locations (e.g., Engel-Cox et al. 2006; Schaap et al. 2009). Schaap et al. (2009) noted that locally derived AOD–PM2.5 relationships cannot be extended easily to other regions because of variation in meteorology and aerosol composition. Unique, local, time-dependent AOD–PM2.5 relationships are necessary to infer global estimates of PM2.5. Ground-based measurements of aerosol vertical profiles and properties have insufficient coverage to estimate global AOD–PM2.5 relationships.
Global chemical transport models (CTMs) resolve atmospheric composition at a resolution of hundreds of kilometers horizontally by hundreds of meters vertically, with a temporal frequency of tens of minutes. Liu et al. (2004) first estimated surface-level PM2.5 from MISR observations by using CTM output to represent local AOD–PM2.5 conversion factors over the contiguous United States. van Donkelaar et al. (2006) extended the approach used by Liu et al. (2004) to estimate PM2.5 from both MODIS and MISR observations and investigated the factors affecting the agreement between AOD and surface-level PM2.5. Statistical models have also been used to relate AOD to PM2.5. For example, Liu et al. (2007) used MISR-retrieved spherical versus nonspherical particle fraction, in addition to model-derived vertical distribution, to separate mineral dust from other aerosol species. More recently, Paciorek and Liu (2009) probed the limitations of using AOD without accounting for vertical distribution or speciation and concluded that agreement with ground-based monitors based on this approach might depend on factors other than satellite observations.
We developed a global satellite-based estimate of surface PM2.5 at a spatial resolution of 0.1° × 0.1°, or approximately 10 km × 10 km at midlatitudes. We developed an approach for combining MODIS and MISR AOD into a single improved estimate of AOD. Using this methodology, we calculated AOD–PM2.5 conversion factors with a global CTM and produced and applied these factors to the AOD. We present a global estimate of PM2.5 concentrations and validate it with ground-based (in situ) observations. We estimate global exposure to outdoor ambient PM2.5 using our satellite-derived product to demonstrate potential application for global health studies. We then examined sources of error.
Materials and Methods
Satellite observations. The MODIS instrument measures a wide range of spatial and spectral information from its orbit aboard the Terra satellite. The near-daily global coverage from the MODIS AOD retrieval (Levy et al. 2007) is advantageous due to frequent measurements. The MISR instrument (Diner et al., 1998), which is also on board Terra, offers smaller spatial and spectral ranges, but views each scene on the earth from nine different angles. This additional angular information allows the MISR AOD retrieval (Diner et al. 2005; Martonchik et al. 2009) to reduce algorithmic assumptions and retrieval bias, as well as obtain information about microphysical properties and plume heights in aerosol source regions (Kahn et al. 2007). Neither instrument can retrieve AOD in cloudy conditions.
We used the MODIS BRDF/Albedo product (MOD43, Collection 5; Schaaf et al. 2002) to distinguish surface types, in conjunction with ground-based retrievals of AOD, and to identify regions of high bias in both MODIS and MISR AOD. We defined these surface types for each month according to the ratio of surface albedo for different wavelengths, similar to assumptions inherent in the MODIS AOD retrieval. We removed AOD that was retrieved from either instrument with an anticipated bias greater than the larger of ± (0.1 or 20%), based on comparison with the Aerosol Robotic Network (AERONET; Holben et al. 1998) sun photometer measurements of AOD. Remaining MODIS and MISR AODs were averaged to produce a single value at a given grid cell. The Supplemental Material (doi:10.1289/ehp.0901623) describes in detail the satellite retrievals and this bias filtration. We restricted our subsequent analysis to locations with at least 50 successful satellite retrievals for 2001–2006 to yield a nearly complete (95%) global geographic coverage.
Estimating PM2.5 from AOD. Estimating ground-level concentrations of dry 24-hr PM2.5 (micrograms per cubic meter) from satellite observations of total-column AOD (unitless) requires a conversion factor that accounts for their spatially and temporally varying relationship:
PM2.5 = η × AOD 
η is a function of the factors that relates 24-hr dry aerosol mass to satellite observations of ambient AOD: aerosol size, aerosol type, diurnal variation, relative humidity, and the vertical structure of aerosol extinction (van Donkelaar et al. 2006). Following the methods of Liu et al. (2004, 2007) and van Donkelaar et al. (2006), we used a global 3-D CTM [GEOS-Chem; geos-chem.org; see Supplemental Material (doi:10.1289/ehp.0901623)] to calculate the daily global distribution of η.
The GEOS-Chem model solves for the temporal and spatial evolution of aerosol (sulfate, nitrate, ammonium, carbonaceous, mineral dust, and sea salt) and gaseous compounds using meteorological data sets, emission inventories, and equations that represent the physics and chemistry of atmospheric constituents. The model calculates the global 3-D distribution of aerosol mass and AOD with a transport time step of 15 min. We applied the modeled relationship between aerosol mass and relative humidity for each aerosol type to calculate PM2.5 for relative humidity values that correspond to surface measurement standards [European Committee for Standardization (CEN) 1998; U.S. Environmental Protection Agency 1997] (35% for theUnited States and Canada; 50% for Europe). We calculated daily values of η as the ratio of 24-hr ground-level PM2.5 for a relative humidity of 35% (U.S. and Canadian surface measurement gravimetric analysis standard)and of 50% (European surface measurement standard) to total-column AOD at ambient relative humidity. We averaged the AOD between 1000 hours and 1200 hours local solar time, which corresponded to the Terra overpass period. We interpolated values of η from 2° × 2.5°, the resolution of the GEOS-Chem simulation, to 0.1° × 0.1° for application to satellite AOD values.
We compared the original MODIS and MISR total-column AOD with coincident ground-based measurements of daily mean PM2.5. Canadian sites are part of the National Air Pollution Surveillance Network (NAPS) and are maintained by Environment Canada (http://www.etc.cte.ec.gc.ca/NAPS/index_e.html). The U.S. data were from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network (http://vista.cira.colostate.edu/improve/Data/data.htm) and from the U.S. Environmental Protection Agency Air Quality System Federal Reference Method sites (http://www.epa.gov/air/data/index.html). Validation of global satellite-derived PM2.5 estimates was hindered by the lack of available surface-measurement networks in many parts of the world. To supplement this lack of available surface measurements, we collected 244 annually representative, ground-based PM2.5 data from both published and unpublished field measurements outside the United States and Canada [see Supplemental Material (doi:10.1289/ehp.0901623)].
In Figure 1A and 1B, we show the mean AOD for retrievals from 1 January 2001 to 31 December 2006 over North America from MODIS and MISR. Both data sets exhibit similar AOD values of 0.15–0.25 over the eastern United States, which reflect a combination of anthropogenic and biogenic sources. Several individual cities can be clearly identified in mean MODIS AOD for the Great Lakes region. A large AOD enhancement over the southwestern United States appears in the MODIS retrievals but is absent from the MISR retrievals. Figure 1C presents the mean combined MODIS and MISR AODs over North America. Our filtration of these two AOD products removes the biased AOD observed by MODIS over the western United States. The combined product is dominated by MODIS in the east because of finer temporal sampling. MISR dominates in the west because of its accuracy.
In Table 1, we provide statistics that compare the spatial variation in 6-year mean AOD retrievals with measurements of daily 24-hr average PM2.5 sampled on the same days as successful satellite observations. Both the MODIS and MISR instruments indicate some relationship between retrieved total-column AOD and in situ PM2.5, both with spatial correlation coefficients of 0.39. A simple average of the daily AOD from both instruments yields a correlation of 0.44. Combining retrievals from these instruments as described in the Methods section increases the correlation to 0.61. Additional information is required to quantitatively estimate PM2.5 concentrations from AOD, as presented below.
Figure 2 shows the annual mean distribution of daily η values used to relate satellite- observed total-column AOD to PM2.5 at 35% relative humidity. Average values of η are typically 20–130 µg/m3. High values of η over regions with large dust concentrations (Prospero et al. 2002) reflect, in part, the low hygroscopicity of dust. Values of η are lower for hygroscopic aerosols, as their dry volume is significantly smaller than under ambient conditions. Ground-level aerosol sources in industrial regions lead to vertical profiles that peak near ground and to moderate values of η. Western North America is characterized by low η, which provides additional insight into the poor AOD–PM2.5 correlations (Engel-Cox et al. 2004b; Hu 2009) associated with this region and in agreement with Liu et al. (2007), who found that transported dust aloft affects the western North America AOD–PM2.5 relationship. η is related to land types only insofar as these are typified by particular aerosol types, meteorology, and vertical structures. Temporal variation in η is considerable.
Figure 3A shows the 6-year mean of 24-hr average satellite-derived surface PM2.5 over North America as calculated from Equation 1 at a daily time scale. A large-scale PM2.5 enhancement is apparent over the eastern United States. The western and northern parts of the continent are generally characterized by low concentrations, with a few exceptions. Geographic mean PM2.5 concentrations over eastern and western North America are 6.9 µg/m3 and 6.2 µg/m3, respectively. Application of η (Figure 2) increased the spatial contrast relative to Figure 1, which reflects ground-level aerosol sources in the east and aerosols aloft in the north and west.
Satellite-derived PM2.5 and comparison with surface measurements. (A) Mean satellite-derived PM2.5 between 2001 and 2006; white space denotes water or < 50 AOD measurements. (B) Average coincident values of both measured and satellite-estimated PM2.5. The solid black line denotes unity; thin dotted lines show uncertainty of ± (1 µg/m3 + 15%); and the dashed line represents the best fit (Hirsh and Gilroy 1984). (C) Positions and mean values of coincidently measured surface sites.
We evaluated the satellite-derived PM2.5 with surface monitors. Figure 3C shows the annual mean of 24-hr PM2.5 concentrations measured with the surface monitors and sampled on the same days as the satellite-derived PM2.5. Ground-level measurements show features similar to our satellite-derived product. Figure 3B quantitatively compares satellite-derived and ground-level measured PM2.5. We found significant cross-sectional correlation between average coincidently sampled satellite-derived and ground-based PM2.5 across North America (r = 0.77; slope = 1.07; bias = –1.75 µg/m3). Many factors contribute to the scatter of points, including differences between what satellite and in situ measurements represent, that do not necessarily indicate errors in either measurement.
Global estimates of PM2.5 concentrations. In Figure 4, we present the 6-year mean of our global satellite-derived PM2.5. This figure, and all subsequent figures, are at 50% relative humidity, which is in agreement with European ground-based measurements. We rejected points created with < 50 values, enabling 95% global geographic coverage. The satellite-derived PM2.5 include an adjustment for discontinuous sampling, as described in the error analysis. The annual mean PM2.5 concentrations vary spatially by more than an order of magnitude. Values are < 10 µg/m3 for large regions of the earth. In contrast, PM2.5 concentrations of 60–90 µg/m3 are found over eastern China, with values > 100 µg/m3 for its major industrial regions. The Indo-Gangetic plain, from New Delhi eastward contains the highest PM2.5 concentrations in India, with values of 80–100 µg/m3, especially in winter (e.g., Di Girolamo et al. 2004). Concentrations elsewhere in northern India are 15–60 µg/m3. The effects of biomass burning on PM2.5 levels are visible in central South America and central Africa, where we estimated concentrations of 10–17 µg/m3. Dust transport in the fine mode is substantial (Jones and Christopher 2007) and contributes to large-scale PM2.5 of approximately 20–50 µg/m3 in the Middle East.
Global satellite-derived PM2.5 averaged over 2001–2006. White space indicates water or locations containing < 50 measurements. Circles correspond to values and locations of comparison sites outside Canada and the United States; the black box outlines European sites.
Figure 4 also shows locations of ground-based measurements and values outside North America that were used for comparison. Despite increased uncertainty because of temporal sampling differences, significant overall agreement exists (r = 0.83; slope = 0.86; intercept = 1.15 µg/m3; n = 244). Similar agreement is obtained when all sites except Europe and North America are considered (r = 0.83; slope = 0.91; intercept = –2.64 µg/m3; n = 84).
Figure 5 overlays contours of population density and surface elevation onto satellite-derived PM2.5 for regions of major anthropogenic sources: eastern North America, western Europe, and eastern Asia. Some relationships are apparent between PM2.5, topography and population. Heavily populated and highly polluted, low-lying regions of eastern China and the Po Valley of northern Italy contrast sharply with neighboring higher altitude regions. The Appalachian Mountains in eastern North America emerged as a relatively clean region. Many PM2.5 enhancements were associated with urban or industrial areas, but these relationships are complex.
Regional satellite-derived PM2.5 concentrations. Columns show mean satellite-derived PM2.5 for 2001–2006 at locations that contain at least 50 measurements. Contours denote population density (left) and surface elevation (right). Crosses indicate city centers. Note the different color scales for each region. Altitude data are from the U.S. Geological Survey (1996).
Error analysis. The dominant sources of error in satellite-derived PM2.5 arose from uncertainties in both AOD retrieval and aerosol vertical structure (van Donkelaar et al. 2006). The residual AOD bias after data filtering is within the larger of ± (0.1 or 20%), as evaluated with ground-based AERONET measurements. We evaluated the GEOS-Chem simulation of the aerosol vertical profile using observations from the Cloud-Aerosol LIDAR and Infrared Pathfinder Satellite Observation (CALIPSO) satellite (Winkler et al. 2007). The GEOS-Chem simulation generally captures to within 5% the fraction of AOD within the boundary layer [see Supplemental Material (doi:10.1289/ehp.0901623)]. We estimate the error in satellite-derived PM2.5 as the change in PM2.5 that occurs when η and AOD are adjusted by their uncertainty, approximated as the GEOS-Chem vertical profile bias and residual satellite AOD bias, respectively.
Figure 6 shows the error distribution of coincidently sampled satellite-derived PM2.5. Arid regions are typically overpredicted and populated regions of East Asia underpredicted. We found that 1 SD of the global error distribution is within ± 15% of the satellite-derived value. We tested this uncertainty estimate by comparing coincident PM2.5 observations for North America (Figure 2) and find that 1 SD of the data lies within ± (1 µg/m3 + 15%). The necessary inclusion of a small absolute term suggests that our uncertainty estimate may be underestimated at low PM2.5 values and supports the presence of a small negative bias (Figure 3B).
Nonuniform and incomplete sampling by satellites have the potential to create bias in long-term mean observations (Levy et al. 2009; Paciorek and Liu 2009). Here we investigate how nonrandom sampling of AOD by satellite observations affects the representation of annual mean PM2.5. The total number of successful satellite retrievals are shown in Figure 7A and are summarized regionally as population-weighted mean in Table 2. Lower sampling was fortuitously collocated with lower population. The global population-weighted mean of observations per 0.1° × 0.1° box was 297. The percent difference between a GEOS-Chem simulation of PM2.5 sampled coincidently with daily satellite-derived PM2.5 versus a complete annual mean of the simulated values is presented in Figure 7B. Most regions exhibited a sampling-induced uncertainty (1 SD) within ± 20% of simulated PM2.5. Regions of low sampling did not necessarily demonstrate enhanced uncertainty and vice versa. Sampling error of satellite-derived PM2.5 is larger in regions influenced by biomass burning, mineral dust, or persistent cloud because of a combination of large seasonal variability and nonrepresentative sampling. We applied the ratio of complete to coincident mean simulated PM2.5 to reduce uncertainty from sampling variability.
Validation of this ratio is inhibited by the lack of in situ measurements in the regions most significantly affected by intermittent sampling. Statistical comparison over the United States and Canada of noncoincident satellite-derived and in situ PM2.5 decreases the agreement relative to a coincident comparison (noncoincident: slope = 1.13; r = 0.70 vs. coincident: slope = 1.07; r = 0.77). This finding supports the need for sampling error correction. Uncertainties derived from both the PM2.5 estimate and sampling can vary substantially on the regional scale. Testing the combined uncertainty of ± 25% from both sources reveals that approximately 1 SD of the North American data falls within this overall error envelope. Globally, the population-weighted mean uncertainty in satellite-derived PM2.5 is 6.7 µg/m3.
Global ambient PM2.5: application to population exposure. Pope et al. (2009) estimated that a decrease of 10 µg/m3 in long-term PM2.5 exposure increases life expectancy by 0.61 ± 0.30 years for persons in the United States. We estimated global long-term exposure to ambient PM2.5 at a spatial resolution of 0.1° using our satellite-derived values for 2001–2006 and the Gridded Population of the World (GPW; Tobler et al. 1997) data for 2005 from the Socioeconomic Data and Applications Center (GPW version 3; http://sedac.ciesin.columbia.edu/). Figure 8 shows the global and regional distributions of long-term ambient PM2.5 exposure; these results are summarized in Table 2. All regions exhibited nonlinear relationships between population and PM2.5 concentrations. Eastern and central Asia have the highest levels of PM2.5 concentrations, with 38–50% of the regional population exceeding the World Health Organization (WHO) Air Quality Interim Target-1 (WHO 2005) of 35 µg/m3. According to the WHO Guidelines, concentrations at this level and higher are associated with an approximately 15% increased risk of mortality, relative to the Air Quality Guideline (AQG) of 10 µg/m3. Globally, 80% of the population live in regions that exceed the AQG. These PM2.5 estimates should be of considerable value for assessing the chronic health impacts of air pollution, especially in regions with sparse ground-based monitoring.
Cumulative distribution of regional, annual mean PM2.5 estimated from satellite-derived PM2.5 at a resolution of 0.1° × 0.1° for 2001–2006. The top axis identifies WHO AQG and Interim Target (IT) associated with each concentration level. Regions are outlined in Figure 6.
A major challenge for global epidemiologic studies and assessments of air pollution health impacts is the lack of representative exposure estimates (Cohen et al. 2004). Extensive ground-based monitoring networks exist in some parts of the world, but major portions of the globe are not covered. The situation is especially acute in developing countries with large populations and high pollution levels and where monitoring with traditional ground-based sampling techniques is limited. Although measurements from ground monitors are currently the gold standard for epidemiologic studies, these are not only sparse, but may represent only a small spatial extent in heterogeneous regions (Chen et al. 2006). Satellite observations offer area-integrated values with global coverage, providing valuable additional information for global health studies.
In our study, we produced a satellite-derived climatology of PM2.5 concentrations. These estimates should facilitate studies of chronic exposure to particulate matter, similar to those already conducted in Europe and North America (e.g., Beelen et al. 2008; Dockery et al. 1993; Pope et al. 2002, 2009), in regions of the world currently without extensive ground-based monitoring networks. Although a growing number of studies are assessing the impacts of short-term exposure to particulate matter in previously underrepresented regions of the world (e.g., Romieu et al. 2009; Wong et al. 2008), studies of long-term exposure also incorporate impacts related to chronic disease and therefore provide a more comprehensive estimate of overall health effects (Kunzli et al. 2001). Our estimates suggest that a concentration of 20 µg/m3 represents the global population-weighted geometric mean of PM2.5 and that 80% of the global population resides in locations where ambient concentrations exceed the WHO AQG of 10 µg/m3. By applying the satellite-derived PM2.5 data set, we also identified global regions and areas of specific concern; half (50%) of the eastern Asian population lives in regions that exceed the WHO Air Quality Interim Target-3 of 35 µg/m3 and are therefore at increased risk from air pollution-related health impacts. These results highlight the potential use of satellite aerosol observations to contribute to studies on the chronic effects of air pollution at regional and global scales.
Several notable developments over previous work were included in our estimates. We combined AODs from two satellite instruments (MODIS and MISR) to improve the correlation of AOD versus ground-based PM2.5 measurements. Extending the satellite data over 6 years (2001–2006) reduced sampling errors. The unprecedented global spatial resolution of 0.1° × 0.1° retains variation relevant to population distribution. A CTM (GEOS-Chem) was applied to account for aerosol vertical distribution, a key factor affecting the relationship between satellite-retrieved, total column AOD and near-surface PM2.5. We found significant spatial agreement between mean coincident satellite-derived and ground-based PM2.5 for North America (slope = 1.07; r = 0.77; n = 1057), as well as evidence of global agreement with noncoincident measurements from published and unpublished data (slope = 0.86; r = 0.83; n = 244). Notably, this level of agreement with ground-based PM2.5 is significantly better than that obtained using a global CTM (GEOS-Chem) without satellite data (Supplemental Material, available online at doi:10.1289/ehp.0901623). Detailed spatial structure in the satellite-derived PM2.5 concentrations reflect multiple influences.
We assessed the uncertainty in the satellite-derived product through comparison with independent observations and error propagation. We estimated our coincident satellite-derived PM2.5 to be accurate at the 1-SD level to within ± 15% of the satellite-derived value using the relative AOD vertical profile measured by the CALIPSO satellite and the total column AOD from ground-based measurements (AERONET). We found evidence that the effect of nonuniform satellite sampling typically biases annual mean satellite-derived PM2.5 by < ± 20% of the satellite-derived value. Larger effects are expected over regions influenced by substantial seasonal variation, by persistent cloud, or for individual, severe pollution events. The overall combined PM2.5 uncertainty of ± 25% indicates a mean global, population-weighted uncertainty in PM2.5 concentration of 6.7 µg/m3.
Additional developments could continue to reduce error in the satellite-derived PM2.5 estimates presented here. Increased satellite coverage would reduce sampling concerns and might allow for satellite-derived PM2.5 to be applied to studies of temporal or spatiotemporal variation. Further improvements to the AOD retrieval (e.g., Drury et al. 2008) would improve accuracy and reduce sampling bias by reducing data rejection. Simulating the AOD–PM2.5 conversion factors at finer spatial resolution would better capture their variability, which is especially important in regions of sharp topographic or emissions gradients. Further development of aerosol speciation capability (e.g., Liu et al. 2007) and satellite-based estimates of additional species, such as NO2 (Lamsal et al. 2008), would be valuable to more specifically estimate pollutant concentrations.
- Beelen R, Hoek G, van den Brandt PA, Goldbohm RA, Fischer P, Schouten LJ, et al. 2008. Long-term effects of traffic-related air pollution on mortality in a Dutch cohort (NLCS-Air Study). Environ Health Perspect 116:196–202.
- Chen C-H, Liu W-L, Chen C-H. 2006. Development of a multiple objective planning theory and system for sustainable air quality monitoring networks. Sci Total Environ 354:1–19.
- Cohen AJ, Anderson HR, Ostro B, Pandey KD, Krzyzanowski M, Kuenzli N, et al. 2004. Mortality impacts of urban air pollution. In: Comparative Quantification of Health Risks: Global and Regional Burden of Disease Due to Selected Major Risk Factors, Vol. 2 (Ezzati M, Lopez AD, Rodgers A, Murray CJL, eds.). Geneva:World Health Organization.
- Di Girolamo L, Bond TC, Bramer D, Diner DJ, Fettinger F, Kahn RA, et al. 2004. Analysis of Multi-Angle Imaging Spectroradiometer (MISR) aerosol optical depths over greater India during winter 2001–2004. Geophys Res Lett 31:L23115; doi: 10.1029/2004GL021273.
- Diner D, Beckert J, Reilly T, Bruegge C, Conel J, Kahn R, et al. 1998. Multiangle Imaging SpectroRadiometer (MISR) instrument description and experiment overview. IEEE Trans Geosci Rem Sens 36:1072–1087.
- Diner DJ, Braswell BH, Davies R, Gobron N, Hu J, Jin Y, et al. 2005. The value of multiangle measurements for retrieving structurally and radiatively consistent properties of clouds, aerosols, and surfaces. Remote Sens Environ 97:495–518.
- Dockery DW, Pope CA, Xu X, Spengler JD, Ware JH, Fay ME, et al. 1993. An association between air-pollution and mortality in 6 U.S. cities. N Engl J Med 329(24):1753–1759.
- Drury E, Jacob DJ, Wang J, Spurr RJD, Chance K. 2008. Improved algorithm for MODIS satellite retrievals of aerosol optical depths over western North America. J Geophys Res 113:D16204.; doi: 10.1029/2007JD009573
- Engel-Cox JA, Hoff RM, Haymet ADJ. 2004. Recommendations on the use of satellite remote-sensing data for urban air quality. J Air Waste Manage Assoc 54(11):1360–1371.
- Engel-Cox JA, Hoff RM, Rogers R, Dimmick F, Rush AC, Szykman JJ, et al. 2006. Integrating lidar and satellite optical depth with ambient monitoring for 3-dimensional particulate characterization. Atmos Environ 40 (8056); doi:10.1016/ j.atmosenv.2006.02.039 [Online 18 October 2006].
- Engel-Cox JA, Holloman CH, Coutant BW, Hoff RM. 2004. Qualitative and quantitative evaluation of MODIS satellite sensor data for regional and urban scale air quality. Atmos Environ 38(16):2495–2509.
- European Committee for Standardization (CEN). 1998. Air quality determination of the PM10 fraction of suspended particulate matter reference method and field test procedure to demonstrate reference equivalence of measurement methods (European Standard EN 12341).
- Gupta P, Christopher SA, Wang J, Gehrig R, Lee Y, Kumar N. 2006. Satellite remote sensing of particulate matter and air quality assessment over global cities. Atmos Environ 40:5880; doi: 10.1016/j.atmosenv.2006.03.016.
- Hirsh RM, Gilroy EJ. 1984. Methods of fitting a straight line to data: examples in water resources. J Am Water Resour Assoc 20(5):705–711.
- Hoff RM, Christopher SA. 2009. Remote sensing of particulate pollution from space: have we reached the promised land? J Air Waste Manag Assoc 59:645–675.
- Holben BN, Eck TF, Slutsker I, Tanre D, Buis JP, Setzer A, et al. 1998. AERONET – a federated instrument network and data archive for aerosol characterization. Remote Sens Environ 66(1):1–16.
- Hu Z. 2009. Spatial analysis of MODIS aerosol optical depth, PM2.5 and chronic coronary heart disease. Int J Health Geogr 8:27; doi: 10.1186/1476-072X-8-27.
- Jones AJ, Christopher SA. 2007. MODIS derived fine mode fraction characteristics of marine, dust, and anthropogenic aerosols over the ocean, constrained by GOCART, MOPITT, and TOMS. J Geophys Res 112 (D22204); doi:10. 1029/2007JD008974 [Online 29 November 2007].
- Kahn RA, Li WH, Moroney C, Diner DJ, Martonchik JV, Fishbein E. 2007. Aerosol source plume physical characteristics from space-based multiangle imaging. J Geophys Res 112:D11205; doi: 10.1029/2006JD007647.
- Koelemeijer RBA, Homan CD, Matthijsen J. 2006. Comparison of spatial and temporal variations of aerosol optical thickness and particulate matter over Europe. Atmos Environ 40:5304; doi: 10.1016/j.atmosenv.2006.04.044.
- Künzli N, Medina S, Kaiser R, Quénel P, Horak FJ, Studnicka M. 2001. Assessment of deaths attributable to air pollution: should we use risk estimates based on time series or on cohort studies? Am J Epidemiol 153(11):1050–1055.
- Lamsal LN, Martin RV, van Donkelaar A, Steinbacher M, Celarier EA, Bucsela E, et al. 2008. Ground-level nitrogen dioxide concentrations inferred from the satellite-borne ozone monitoring instrument. J Geophys Res 113:16308; doi: 10.1029/2007JD009235.
- Levy R, Leptoukh G, Kahn R, Zubko V, Gopalan A, Remer L. 2009. A critical look at deriving monthly aerosol optical depth from satellite data. IEEE Trans Geosci Remote Sens 47:2942–295; doi: 10.1109/TGRS.2009.2013842.
- Levy RC, Remer LA, Mattoo S, Vermote EF, Kaufman YJ. 2007. Second-generation operational algorithm: retrieval of aerosol properties over land from inversion of moderate resolution imaging spectroradiometer spectral reflectance. J Geophys Res 112:13211; doi: 10.1029/2006JD007811.
- Liu Y, Koutrakis P, Kahn R. 2007. Estimating fine particulate matter component concentrations and size distribution using satellite-retrieved fractional aerosol optical depth: Part 1 – method development. J Air Waste Manage Assoc 57(11):1351–1359.
- Liu Y, Park RJ, Jacob DJ, Li QB, Kilaru V, Sarnat JA. 2004. Mapping annual mean ground-level PM2.5 concentrations using multiangle imaging spectroradiometer aerosol optical thickness over the contiguous United States. J Geophys Res 109:22206; doi: 10.1029/2004JD005025.
- Liu Y, Sarnat JA, Kilaru A, Jacob DJ, Koutrakis P. 2005. Estimating ground-level PM2.5 in the eastern United States using satellite remote sensing. Environ Sci Technol 39(9):3269–3278.
- Martin RV. 2008. Satellite remote sensing of surface air quality. Atmos Environ 42:7823–7843.
- Martonchik JV, Kahn RA, Diner DJ. 2009. Retrieval of aerosol properties over land using MISR observations. In: Satellite Aerosol Remote Sensing Over Land (Kokhanovsky AA, de Leeuw G, eds). Berlin:Springer, 267–291.
- McDonnell WF, Nishino-Ishikawa N, Petersen FF, Chen LH, Abbey DE. 2000. Relationships of mortality with the fine and coarse fractions of long-term ambient PM10 concentrations in nonsmokers. J Expos Anal Environ Epidemiol 10(5):427–436.
- Paciorek CJ, Liu Y. 2009. Limitations of remotely-sensed aerosol as a spatial proxy for fine particulate matter. Environ Health Perspect 117:904–909.
- Pope CA III, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, et al. 2002. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. JAMA 287(9):1132–1141.
- Pope CA III, Ezzati M, Dockery DW. 2009. Fine-particulate air pollution and life expectancy in the United States. N Engl J Med 360(4):376–386.
- Prospero JM, Ginoux P, Torres O, Nicholson SE, Gill TE. 2002. Environmental characterization of global sources of atmospheric soil dust identified with the nimbus 7 Total Ozone Mapping Spectrometer (TOMS) absorbing aerosol product. Rev Geophys 40:1002; doi: 10.1029/2000RG000095.
- Romieu I, Gouveia N, Cifuentes L, Ponce de Leon A, Junger W, Miranda V, et al. 2009. In: Health Effects Institute Annual Conference 2009 Program and Abstracts. Boston:Health Effects Institute. Mortality effects of air pollution in Latin American cities: results from the Escala study [Abstract]. Available: http://www.healtheffects.org/Pubs/AnnualConferenceProgram2009.pdf [accessed 13 May 2010].
- Schaaf CB, Gao F, Strahler AH, Lucht W, Li X, Tsang T, et al. 2002. First operation BRDF, albedo nadir relectance products from MODIS. Remote Sens Environ 83:135–148.
- Schaap M, Apituley A, Timmermans RMA, Koelemeijer RBA, de Leeuw G. 2009. Exploring the relation between aerosol optical depth and PM2.5 at Cabauw, the Netherlands. Atmos Chem Phys 9:909–925.
- Tobler W, Deichmann U, Gottsegen J, Maloy K. 1997. World population in a grid of spherical quadrilaterals. Int J Popul Geogr 3:203–225.
- U.S. Environmental Protection Agency. 1997. Reference Method for the Determination of Fine Particulate Matter as PM2.5 in the Atmosphere. 40CFR50, Appendix L.
- U.S. Geological Survey. 1996. Earth Resources Observation and Science (EROS) Center. GTOPO30. Available: http://eros.usgs.gov/#/Find_Data/Products_and_Data_Available/gtopo30_info) [accessed 11 May 2010].
- van Donkelaar A, Martin RV, Park RJ. 2006. Estimating ground-level PM2.5 using aerosol optical depth determined from satellite remote sensing. J Geophys Res 111:21201; doi: 10.1029/2005JD006996.
- Wang J, Christopher SA. 2003. Intercomparison between satellite- derived aerosol optical thickness and PM2.5 mass: implications for air quality studies. Geophys Res Lett 30:2095; doi: 10.1029/2003GL018174.
- WHO (World Health Organization). 2005. Air Quality Guidelines Global Update 2005. Copenhagen, Denmark:WHO/Europe.
- Winkler DM, Hunt WH, McGill MJ. 2007. Initial performance assessment of CALIOP, Geophys Res Lett 34:L19803; doi: 10.1029/2007GL030135.
- Wong C-M, Vichit-Vadakan N, Kan H, Qian Z, Teams PP. 2008. Public Health and Air Pollution in Asia (PAPA): a multicity study of short-term effects of air pollution on mortality. Environ Health Perspect 116:1195–1202.
Come visit with EHP’s Science Editor, Jane Schroeder, at the 2016 ISEE Meeting in Rome, Italy, from 1–4 September 2016. This is a great opportunity to meet a member of our team, and to learn more about the journal. A number of our Associate Editors will be in attendance as well. Jane is also co-hosting a half-day writing and publishing workshop on 4 September, following the conference.
EHP is pleased to present the abstracts from the 28th annual meeting of the International Society for Environmental Epidemiology (ISEE), held in Rome, Italy, 1–4 September 2016, and hosted by the Department of Epidemiology Lazio Regional Health Service, ASL Roma 1, and the Italian Epidemiological Association. The focus of this year’s conference is current and future challenges in exposure assessment, study design, and data analyses.
Featured Children’s Health
Birgit Claus Henn, Adrienne S. Ettinger, Marianne R. Hopkins, Rebecca Jim, Chitra Amarasiriwardena, David C. Christiani, Brent A. Coull, David C. Bellinger, and Robert O. Wright
Diane Gilbert-Diamond, Jennifer A. Emond, Emily R. Baker, Susan A. Korrick, and Margaret R. Karagas
Benjamin B. Green, Margaret R. Karagas, Tracy Punshon, Brian P. Jackson, David J. Robbins, E. Andres Houseman, and Carmen J. Marsit