UBC Faculty Research and Publications

Global Estimates of Ambient Fine Particulate Matter Concentrations from Satellite-Based Aerosol Optical… van Donkelaar, Aaron; Martin, Randall V.; Brauer, Michael; Kahn, Ralph; Levy, Robert; Verduzco, Carolyn; Villeneuve, Paul J. Jun 30, 2010

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


52383-van_Donkelaar_A_et_al_Global_estimates_ambient_fine_particulate.pdf [ 12.67MB ]
JSON: 52383-1.0103477.json
JSON-LD: 52383-1.0103477-ld.json
RDF/XML (Pretty): 52383-1.0103477-rdf.xml
RDF/JSON: 52383-1.0103477-rdf.json
Turtle: 52383-1.0103477-turtle.txt
N-Triples: 52383-1.0103477-rdf-ntriples.txt
Original Record: 52383-1.0103477-source.json
Full Text

Full Text

Environmental Health Perspectives • volume 118 | number 6 | June  2010 847ResearchChronic 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 environ-mental 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 meas-urement 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 con-centration, 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 satel-lite remote sensing could substantially improve estimates of population exposure to PM2.5. Since the mid 2000s, the MODIS (Moder-ate Resolution Imaging Spectro radiometer) 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 atmo-spheric aerosol.Many studies have investigated the rela-tionship 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. 2008). Schaap et al. (2008) 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 esti-mates 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 hor-izontally by hundreds of meters vertically, with a temporal frequency of tens of min-utes. 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 fac-tors 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 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.caSupplemental Material is available online (doi:10.1289/ehp.0901623 via http://dx.doi.org/).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.Global Estimates of Ambient Fine Particulate Matter Concentrations from Satellite-Based Aerosol Optical Depth: Development and ApplicationAaron van Donkelaar,1 Randall V. Martin,1,2 Michael Brauer,3 Ralph Kahn,4 Robert Levy,4 Carolyn Verduzco,1 and Paul J. Villeneuve5,61Department 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, CanadaBackground: 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 opti-cal depth (AOD) from the MODIS (Moderate Resolution Imaging Spectroradiometer) and MISR (Multiangle Imaging Spectroradiometer) satellite instruments and coincident aerosol vertical pro-files 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 esti-mate 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 sam-pling. 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. Environ Health Perspect 118:847–855 (2010). doi:10.1289/ehp.0901623 [Online 16 March 2010]van Donkelaar et al.848 volume 118 | number 6 | June 2010 • Environmental Health Perspectivesand 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 esti-mate 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 demon-strate potential application for global health studies. We then examined sources of error.Materials and MethodsSatellite observations. The MODIS instru-ment measures a wide range of spatial and spectral information from its orbit aboard the Terra satellite. The near-daily global cover-age from the MODIS AOD retrieval (Levy et al. 2007) is advantageous due to frequent measure ments. 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 differ-ent angles. This additional angular information allows the MISR AOD retrieval (Diner et al. 2005; Martonchik et al. 2009) to reduce algo-rithmic 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 accord-ing 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 instru-ment with an anticipated bias greater than the larger of ± (0.1 or 20%), based on com-parison with the Aerosol Robotic Network (AERONET; Holben et al. 1998) sun pho-tometer 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 [1]η is a function of the factors that relates 24-hr dry aerosol mass to satellite observa-tions of ambient AOD: aerosol size, aero-sol type, diurnal variation, relative humidity, and the vertical structure of aerosol extinc-tion (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 rep-resent the physics and chemistry of atmo-spheric 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 Figure 1. Mean AOD for 2001–2006 from the MODIS (A) and MISR (B) satellite instruments. (C) Data from the combined product developed here. White space denotes water or < 50 measurements.MODISMISRCombined0. (unitless)ABCSatellite-based global estimates of PM2.5Environmental Health Perspectives • volume 118 | number 6 | June 2010 849between 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 sur-face measurement standard) to total-column AOD at ambient relative humidity. We aver-aged 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 coinci-dent 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 avail-able 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)].ResultsIn 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 combi-nation of anthropogenic and biogenic sources. Several individual cities can be clearly identi-fied 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 tem-poral sampling. MISR dominates in the west because of its accuracy.In Table 1, we provide statistics that com-pare 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 cor-relation coefficients of 0.39. A simple aver-age of the daily AOD from both instruments yields a correlation of 0.44. Combining retriev-als 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 distribu-tion 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 Table 1. Comparison of coincidently sampled 6-year mean measurementsa of daily 24-hr average PM2.5 with AOD and satellite-derived PM2.5. Retrievals Slopeb Intercept r nMODIS AOD 0.020 0.10 0.39 1,218MISR AOD 0.010 0.11 0.39 353Average AOD 0.015 0.06 0.44 1,236Combined AOD 0.017 0.10 0.61 1,057Satellite-derived PM2.5 1.066 –1.75 0.77 1,057aA minimum of 50 measurements is required for each point. bCalculated with reduced major-axis linear regression (Hirsh and Gilroy 1984).Figure 2. Annual mean η (ratio of PM2.5 to AOD) for 35% relative humidity. White space indicates water.20 40 60 80η (µg/m3)100 120 140van Donkelaar et al.850 volume 118 | number 6 | June 2010 • Environmental Health PerspectivesFigure 3. 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 dot-ted 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 meas ured surface sites.181614121086420Satellite-derived PM2.5 (µg/m3 )Measured in situMeasured in situSatellite-derivedSatellite-derived PM2.5 (µg/m3 )3025201510500 10y = 1.07x – 1.75; r = 0.77n = 1,05720 30AB CFigure 4. 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.20Satellite-derived PM2.5 (µg/m3)1050 15 8050Satellite-based global estimates of PM2.5Environmental Health Perspectives • volume 118 | number 6 | June 2010 851low η, 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 rela-tionship. η 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 charac terized by low concentrations, with a few excep-tions. Geographic mean PM2.5 concentra-tions over eastern and western North America are 6.9 µg/m3 and 6.2 µg/m3, respectively. Application of η (Figure 2) increased the spa-tial contrast relative to Figure 1, which reflects ground-level aerosol sources in the east and aerosols aloft in the north and west.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 sam-pled on the same days as the satellite-derived PM2.5. Ground-level measurements show fea-tures similar to our satellite-derived product. Figure 3B quantitatively compares satellite-derived and ground-level measured PM2.5. We found significant cross-sectional correla-tion 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 measure-ments represent, that do not necessarily indi-cate 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 fig-ure, 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 adjust-ment 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 east-ward contains the highest PM2.5 concentra-tions 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 contrib-utes to large-scale PM2.5 of approximately 20–50 µg/m3 in the Middle East.Figure 4 also shows locations of ground-based measurements and values outside North America that were used for comparison. Despite increased uncertainty because of tem-poral sampling differences, significant over-all agreement exists (r = 0.83; slope = 0.86; intercept = 1.15 µg/m3; n = 244). Similar agree-ment 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 anthropo-genic 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 con-trast sharply with neighboring higher alti-tude regions. The Appalachian Mountains in Figure 5. 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).Eastern North AmericaEastern AsiaWestern EuropePopulation density Surface elavation1816141210864200055101015154570100202535Satellite-derived PM2.5 (µg/m3 )Satellite-derived PM2.5 (µg/m3 )Satellite-derived PM2.5 (µg/m3 )50/km250/km250 m100 m300 m500 m50 m100 m100 m300 m500 m500 m500/km2500/km2100/km2100/km21,000/km2500/km2van Donkelaar et al.852 volume 118 | number 6 | June 2010 • Environmental Health Perspectiveseastern North America emerged as a relatively clean region. Many PM2.5 enhancements were associated with urban or industrial areas, but these relationships are complex.Error analysis. The dominant sources of error in satellite-derived PM2.5 arose from uncertainties in both AOD retrieval and aero-sol vertical structure (van Donkelaar et al. 2006). The residual AOD bias after data fil-tering is within the larger of ± (0.1 or 20%), as evaluated with ground-based AERONET meas urements. We evaluated the GEOS-Chem Figure 6. Estimate of the satellite-derived PM2.5 bias, defined as (satellite-derived PM2.5 – truth) ÷ truth. Boxed areas outline the regions used in Figure 8; see Supplemental Material (doi:10.1289/ehp.0901623) for the two subregions. PM2.5 bias estimate (%)100402510–10–40–100 –25Figure 7. Satellite-derived PM2.5 sampling and its estimated induced uncertainty. (A) Total number of values used from satellite per 0.1° grid box. (B) Percentage change in average coincidently sampled simulated PM2.5 concentrations relative to a full-year average.1,000800550300200100–60–35–101060350ObservationsSampling-induced uncertainty (%)ABSatellite-based global estimates of PM2.5Environmental Health Perspectives • volume 118 | number 6 | June 2010 853simulation 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 distri-bution 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 inves-tigate how nonrandom sampling of AOD by satellite observations affects the representa-tion 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 sam-pled 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 simu-lated 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 persis-tent 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 rela-tive 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 approxi-mately 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) esti-mated 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 expo-sure 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 Figure 8. 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.1009080706050403020100Population (%)AQG5 10 15 25PM2.5 exposure (µg/m3)35 50 100IT-3 IT-2 IT-1WHO air quality guideline and interim targetsWorldNorth AfricaSouth AfricaWestern EuropeEastern EuropeEastern North AmericaSouth AmericaCentral AsiaEastern AsiaSouth Asia/AustraliaWestern North AmericaTable 2. Regional PM2.5 statistics, number of observations, and population in excess of WHO AQG and Interim Targetsa. Population-weighted total observationsPopulation-weighted statistics [µg/m3] Population [million people (%)]Region Mean SD GM GSD Total AQG (10 µg/m3) IT-3 (15 µg/m3) IT-2 (25 µg/m3) IT-1 (35 µg/m3)World 27 23 20 2.3 297 6,400 5,100 (81) 4,000 (63) 2,300 (37) 1,600 (25)Eastern Asia 44 29 34 2.2 270 1,900 1,800 (93) 1,600 (83) 1,200 (65) 940 (50)Central Asia 31 16 27 1.8 230 1,500 1,400 (94) 1,300 (86) 780 (54) 560 (38)North Africa 26 10 24 1.5 195 540 530 (98) 450 (85) 250 (47) 8.7 (17)South Africa 11 6 9 1.9 368 540 240 (47) 110 (21) 7.7 (1.5) 0.1 (0.0)Eastern Europe 15 5 14 1.4 437 450 400 (88) 210 (47) 14 (3.0) 1.8 (0.4)South America 7 4 5 2.0 361 400 52 (13) 18 (4.5) 5.3 (1.3) 0.0 (0.0)Eastern North America 13 5 12 1.5 476 350 250 (72) 110 (32) 17 (4.7) 0.0 (0.0)South Asia/Australia 12 6 10 1.8 304 310 180 (58) 72 (24) 13 (4.4) 1.3 (0.4)Western Europe 17 5 16 1.4 311 260 250 (94) 170 (63) 15 (5.9) 3.7 (1.4)Western North America 11 5 10 1.6 366 120 56 (46) 22 (18) 0.3 (0.2) 0.0 (0.0)Abbreviations: GM, geometric mean; GSD, geometric standard deviation; IT, interim target.aData from WHO (2005). Regions are outlined in Figure 6.van Donkelaar et al.854 volume 118 | number 6 | June 2010 • Environmental Health Perspectivessummarized 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 concen-trations, with 38–50% of the regional popula-tion 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 approxi-mately 15% increased risk of mortality, rela-tive 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.DiscussionA major challenge for global epidemiologic studies and assessments of air pollution health impacts is the lack of representative expo-sure 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 cur-rently the gold standard for epidemiologic studies, these are not only sparse, but may rep-resent only a small spatial extent in hetero-geneous 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, simi-lar 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 under-represented 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 concentra-tions 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 previ-ous work were included in our estimates. We combined AODs from two satellite instru-ments (MODIS and MISR) to improve the correlation of AOD versus ground-based PM2.5 measurements. Extending the satel-lite data over 6 years (2001–2006) reduced sampling errors. The unprecedented global spatial resolution of 0.1° × 0.1° retains varia-tion 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 agree-ment 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 agree-ment with ground-based PM2.5 is signifi-cantly 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 con-centrations reflect multiple influences.We assessed the uncertainty in the satel-lite-derived product through comparison with independent observations and error propaga-tion. 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 meas-ured by the CALIPSO satellite and the total column AOD from ground-based measure-ments (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 improve-ments 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 fac-tors 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 addi-tional species, such as NO2 (Lamsal et al. 2008), would be valuable to more specifically estimate pollutant concentrations.RefeRencesBeelen 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 pol-lution. 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 [Online 15 December 2004].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 mor-tality 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 [Online 21 August 2008].Engel-Cox JA, Hoff RM, Haymet ADJ. 2004a. 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. 2004b. 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 demon-strate 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 [Online 20 July 2006].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.Satellite-based global estimates of PM2.5Environmental Health Perspectives • volume 118 | number 6 | June 2010 855Hu 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 [Online 12 May 2009].Jones AJ, Christopher SA. 2007. MODIS derived fine mode fraction characteristics of marine, dust, and anthropo-genic 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 characteris tics from space-based multiangle imaging. J Geophys Res 112:D11205; doi:10.1029/2006JD007647 [Online 7 June 2007].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 [Online 21 June 2006].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 [Online 28 August 2008].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–2956; doi:10.1109/TGRS.2009.2013842 [Online 31 March 2009].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 [Online 13 July 2007].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 concentra-tions using multiangle imaging spectroradiometer aero-sol optical thickness over the contiguous United States. J Geophys Res 109:22206; doi:10.1029/2004JD005025 [Online 24 November 2004].Liu Y, Sarnat JA, Kilaru A, Jacob DJ, Koutrakis P. 2005. Estimating ground-level PM2.5 in the eastern United States using satel-lite 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 non-smokers. J Expos Anal Environ Epidemiol 10(5):427–436.Paciorek CJ, Liu Y. 2009. Limitations of remotely-sensed aero-sol 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 [Online 4 September 2002].Romieu I, Gouveia N, Cifuentes L, Ponce de Leon A, Junger W, Miranda V, et al. 2010. Mortality effects of air pollution in Latin American cities: results from the Escala study [Abstract]. In: Health Effects Institute Annual Conference 2009 Program and Abstracts. Boston:Health Effects Institute, 46. 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 prod-ucts from MODIS. Remote Sens Environ 83:135–148.Schaap M, Apituley A, Timmermans RMA, Koelemeijer RBA, de Leeuw G. 2009. Exploring the relation between aero-sol 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 [Online 2 November 2006].Wang J, Christopher SA. 2003. Intercomparison between satellite- derived aerosol optical thickness and PM2.5 mass: impli-cations for air quality studies. Geophys Res Lett 30:2095; doi:10.1029/2003GL018174 [Online 6 November 2003].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 [Online 3 October 2007].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.


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items