Skip to main content

Spatial analysis of MODIS aerosol optical depth, PM2.5, and chronic coronary heart disease



Numerous studies have found adverse health effects of acute and chronic exposure to fine particulate matter (PM2.5). Air pollution epidemiological studies relying on ground measurements provided by monitoring networks are often limited by sparse and unbalanced spatial distribution of the monitors. Studies have found correlations between satellite aerosol optical depth (AOD) and PM2.5 in some land regions. Satellite aerosol data may be used to extend the spatial coverage of PM2.5 exposure assessment. This study was to investigate correlation between PM2.5 and AOD in the conterminous USA, to derive a spatially complete PM2.5 surface by merging satellite AOD data and ground measurements based on the potential correlation, and to examine if there is an association of coronary heart disease with PM2.5.


Years 2003 and 2004 daily MODIS (Moderate Resolution Imaging Spectrometer) Level 2 AOD images were collated with US EPA PM2.5 data covering the conterminous USA. Pearson's correlation analysis and geographically weighted regression (GWR) found that the relationship between PM2.5 and AOD is not spatially consistent across the conterminous states. The average correlation is 0.67 in the east and 0.22 in the west. GWR predicts well in the east and poorly in the west. The GWR model was used to derive a PM2.5 grid surface using the mean AOD raster calculated using the daily AOD data (RMSE = 1.67 μg/m3). Fitting of a Bayesian hierarchical model linking PM2.5 with age-race standardized mortality rates (SMRs) of chronic coronary heart disease found that areas with higher values of PM2.5 also show high rates of CCHD mortality: = 0.802, posterior 95% Bayesian credible interval (CI) = (0.386, 1.225).


There is a spatial variation of the relationship between PM2.5 and AOD in the conterminous USA. In the eastern USA where AOD correlates well with PM2.5, AOD can be merged with ground PM2.5 data to derive a PM2.5 surface for epidemiological study. The study found that chronic coronary heart disease mortality rate increases with exposure to PM2.5.


Numerous epidemiological studies have indicated that exposure to fine particulate matter (particles smaller than 2.5 micrometers, PM2.5) is associated with asthma, respiratory infections, lung cancer, cardiovascular problems, and premature death [15]. A few have examined coronary heart disease, finding evidence for acute effects on mortality and hospital admissions [68]. Recently, attention has focused on whether there is an association between chronic exposure to air pollution and coronary heart disease [9]. An ecological study at the census enumeration district level found an association between nitrogen oxides, and to a lesser extent particulate matter (PM10) and carbon monoxide, and coronary heart disease mortality in Sheffield, UK [9]. The objective of this study was to examine if there is an association of coronary heart disease with chronic exposure to PM2.5. The study adopted an ecological method using aggregate disease data at the county level and average PM2.5 concentration value for each county.

Air pollution epidemiological studies often rely on ambient observations to provide metrics of exposure, as in studies of PM and cardiovascular diseases [2, 4, 912]. Methods of exposure assessment in those studies include averaging multiple monitors within each enumeration unit or study site [4, 10, 11], assigning the exposure value of the nearest monitor to each case/control [2, 12] and spatial interpolation/modelling method [9]. Ground monitoring data often lacks spatially complete coverage. Ground monitors are rare in rural areas. Assessment of the exposure to air pollution using in situ observations is hampered by the sparse and unbalanced spatial distribution of the monitors.

The repetitive and broad-area coverage of satellites may allow atmospheric remote sensing to offer a unique opportunity to monitor air quality at continental, national and regional scales. Recent studies have established quantitative relationships between satellite derived aerosol optical depth (AOD), which describes the mass of aerosols in an atmospheric column, and PM2.5 using linear regression models [1323]. Except long-range dust or pollution transport events, AOD is dominated by near-surface emissions sources [24]. AOD retrieved at visible wavelengths is most sensitive to particles between 0.1 and 2 micrometers [25]. Several studies have merged AOD with ground PM2.5 measurement to derive PM2.5 surfaces [2628]. A study in a region centered in Massachusetts [26] examined the benefits of using AOD retrieved by the Geostationary Operational Environmental Satellite (GOES) in conjunction with land use and meteorological information to estimate ground-level PM2.5 concentrations. Another project [27] combined MODIS (Moderate Resolution Imaging Spectrometer) AOD data with US EPA PM2.5 data to estimate a PM2.5 surface in Atlanta, Georgia. Existing studies estimating PM2.5 surfaces using AOD data use uniform linear relationships between AOD and PM2.5. However, studies have found that correlation between PM2.5 and AOD is not spatially consistent [29] due to variation in terrain, land cover, selection of aerosol model in the AOD retrieval algorithm and meteorological factors.

This paper quantitatively examined the relationship between PM2.5 ground measurements and MODIS AOD data in the conterminous USA using Pearson's correlation analysis and geographically weighted regression (GWR). For the region with high correlations, the GWR model was used to calculate a PM2.5 surface based on the AOD data and the spatially varying relationships between PM2.5 and AOD. A Bayesian hierarchical simulation model was used to link chronic coronary heart disease mortality rate with PM2.5.


MODIS data

The MODIS sensor flies on polar-orbiting and sun-synchronous Terra and Aqua satellites. MODIS performs measurements in the visible to thermal infrared spectrum region. The MODIS sensor was expected to be the key for monitoring global aerosol properties. Not only have MODIS aerosol products been used to answer scientific questions about radiation and climate, they are being used for applications not previously intended, including monitoring surface air quality for health [15, 16, 3032].

One of the fundamental aerosol products from MODIS is spectral AOD. The MODIS level 2 files are produced everyday at the spatial resolution of a 10 km × 10 km (at nadir) and represent the first level of MODIS aerosol retrieval. The latest version of the MODIS aerosol retrieval algorithm is Collection 5 (C005) [33]. The aerosol retrieval makes use of seven wavelength bands (channels 1 to 7), and a number of other bands to help with cloud rejection and other screening procedures. The aerosol algorithm relies on calibrated, geolocated reflectance Level 1B data. In addition, the MODIS algorithm uses cloud mask product [34], atmospheric profile product and ancillary data from NCEP (National Center for Environmental Prediction) analyses, including hourly meteorological analysis and daily ozone analysis. Different dynamic aerosol models (biomass burning, dust aerosol, and aerosol from industrial/urban origin) are used to determine the aerosol optical properties used in the algorithm for different parts of the land. Over the land of the conterminous US, the aerosol retrieval algorithm uses the urban/industrial aerosol model for the east and the biomass burning/dust model for the west. The splitting line is at approximately -100° longitude.

Daily level 2 MODIS data (2003–2004) were obtained from the NASA Level 1 and Atmosphere Archive and Distribution System (LAADS Web) [35]. A two-year average AOD raster data layer (10 km by 10 km grid) was calculated. Data from both Terra and Aqua satellites were used. MODIS AOD data are not available every day due to cloud cover. Data for cold seasons (October to March) were not used in the two-year average calculation and correlation analysis. Cloud cover, snow reflectivity, and diminished vertical mixing all reduce the accuracy of ground-level pollutant levels measured in winter. During warm seasons, vertical columns in the atmosphere are more integrated. AOD measures correlate best with ground-based monitoring in warm months, likely because of stronger boundary layer mixing during the warmer months [36].

Pearson's correlation and geographically weighted regression

PM2.5 ground data was obtained from US EPA Air Quality System (AQS) online Data Mart [37]. There were 877 monitoring sites. PM2.5 data was collated with AOD both temporally and spatially. For each MODIS AOD image scene and each monitor within the scene, PM2.5 measurement within one hour of the imaging time was assigned to the pixel containing the monitor. Pearson's correlation value was calculated for each monitoring site. The Pearson's correlation analysis examined the temporal relationship between PM2.5 and AOD for each site. The relationship shows how AOD changes as PM2.5 changes over time at a sampling site. To visualize the spatial variation of the correlations, inverse distance weighted interpolation was used to create a correlation surface. It should be noted that the correlation surface was used for visualization purpose only. It was not intended for the use of quantitative analysis, since the interpolation was based on the unbalanced distribution of the monitors and the accuracy was not spatially consistent.

A geographically weighted regression (GWR) model was also fitted to examine the relationship between PM2.5 (dependent variable) and AOD (independent variable) using the two-year average PM2.5 and AOD values. The mean PM2.5 value was geo- collated with the mean AOD value for each of the 877 sites. GWR is a local form of linear regression used to model spatially varying relationships [38]. GWR is one of several spatial regression techniques, increasingly used in geography and other disciplines. GWR provides a local model of the variable or process we are trying to understand/predict by fitting a regression equation to every feature in the dataset. GWR constructs these separate equations by incorporating the dependent and explanatory variables of features falling within a kernel of each target feature. GWR accounts for the effect of local multicolliearity which occurs when the values for a particular explanatory variable cluster spatially. GWR creates coefficient raster surfaces for the model intercept and each explanatory variable. GWR used in the study analyzed the spatial varying relationship between two-year mean AOD and mean PM2.5. Like the correlation surface created from Pearson's correlation analysis, GWR revealed how the relationship between AOD and PM2.5 changes across the space. Pearson's correlation analysis used daily time-series of data, and GWR used mean data. While both analyses revealed the spatial variation of the relationship between AOD and PM2.5, GWR also calculated a mean PM2.5 surface that was needed for disease modeling since disease data was a two-year aggregate. In this study, a fixed kernel was used to fit the regression. The bandwidth value was chosen by using the corrected Akaike Information Criterion (AICc). This option tries to identify an optimal fixed distance. For the region with high correlations, the coefficient raster surfaces were applied to the two-year mean AOD raster to calculate a PM2.5 surface. The accuracy of the estimated PM2.5 surface was assessed by comparison to the two-year mean PM2.5 ground measurements.

Disease data

Modelling of disease and PM2.5 was focused on the region where higher positive correlation between AOD and PM2.5 was found. Chronic coronary heart disease (CCHD) (International Classification of Disease version 10 codes: I25.0-I25.6, I25.8, I25.9) data at the county level were extracted for the period from 2003 to 2004 from the National Center for Health Statistics Compressed Mortality File 1999–2005 in the CDC WONDER online database [39]. CCHD mortality count and population at risk was retrieved by county, race (White, Black or African American, Other race) and age (5 years interval from 1–4 to 85+). Aggregated CCHD mortality count and population at risk were also retrieved by race and age groups for the focused study area to be used as standard population in calculating CCHD mortality rate adjusting for race and age effects.

Race and age adjusted rate was calculated using indirect standardization [40] for each county. Rate adjustment is a technique for removing the effects of race and age from crude rates, so as to allow meaningful comparisons across populations with different underlying race and age structures. In the US, researchers normally use the year 2000 US decennial census data as a standard population to standardize rates. In the present study, the use of an internal standard population consisting of race-age-specific rates from a super-population containing the regions to compare would be better than an external standard based on individuals from an entirely separate population [41]. The indirect standardization first calculated expected number of CCHD deaths for each county. The calculated death count is the number of cases that would be expected in the study population if people in the study population contracted the disease at the same rate as people in the standard population. Standardized mortality rates (SMRs) were calculated by dividing the observed count by the expected value. According to Curtin and Klein [42], one of the problems with rate adjustment is that rates based on small numbers of deaths will exhibit a large amount of random variation. In the aggregate CCHD mortality count and population data set, counts of twenty or less were flagged as statistically "unreliable". CCHD counts were "suppressed" when the data meets the criteria for confidentiality constraints. Counts for counties with census year populations of less than 100,000 were replaced with "suppressed" if the number of cases is five or less and the count is based on only one or two years of data. All unreliable and suppressed data were not used in calculating standardized rates and spatial analysis and modelling thereafter.

Bayesian hierarchical modelling of CCHD and PM2.5

A GIS zonal statistical function was used to calculate the mean PM2.5 value for each county by averaging PM2.5 values of all cells within the county. A Bayesian hierarchical model was used to explore the association between CCHD mortality and PM2.5. Simulation-based algorithms for Bayesian inference allow us to fit very complicated hierarchical models, including those with spatially correlated random effects. In this ecological study, there could exist spatial autocorrelation within values of the CCHD mortality and PM2.5. The following model was fitted allowing a convolution prior for the random effects:

where i is the index for a county, O is observed CCHD death count, E is expected death count reflecting race-age-standardized values. For model specification, an improper (flat) prior for the intercept parameters β0 and a uniform prior distribution for the fixed-effect (β1) were assumed. By fixed effect we mean it applies equally to all the counties. Two sets of county-specific random effects were included in the model. The first set b i is spatially structured random effects assigned an intrinsic Gaussian conditional auto-regression (CAR) prior distribution [43]. The second set of random effects h i is assigned an exchangeable (non-spatial) normal prior. The random effect for each county is thus the sum of a spatially structured component b i and an unstructured component h i . This is termed a convolution prior [43, 44]. The model is more flexible than assuming only CAR random effects, since it allows the data to decide how much of the residual disease risk is due to spatially structured variation, and how much is unstructured over-dispersion. To complete the model specification, conjugate inverse-gamma prior distributions were assigned to the variance parameters associated with the exchangeable and/or CAR priors.

The Markov chain Monte Carlo (MCMC) simulation computation technique and Gibbs sampling algorithm were used to fit the Bayesian model. Having specified the model as a full joint distribution on all quantities, whether parameters or observables, we wish to sample values of the unknown parameters from their conditional (posterior) distribution given those stochastic nodes that have been observed. MCMC methods perform Monte Carlo simulations generating parameter values from Markov chains having stationary distributions identical to the joint posterior distribution of interest. After these Markov chains converge to a stationary distribution, the simulated parameter values represent a correlated sample of observations from the posterior distribution. The basic idea behind the Gibbs sampling algorithm is to successively sample from the conditional distribution of each node given all the others. Under broad conditions this process eventually provides samples from the joint posterior distribution of the unknown quantities. Summaries of the post-convergence MCMC samples provide posterior inference for model parameters. The result of such analysis is the posterior distribution of an intensity function with covariate effects.

The model was fitted using the WinBUGS software – an interactive Windows version of the BUGS (Bayesian inference Using Gibbs Sampling) program for Bayesian analysis of complex statistical models using MCMC techniques [45]. A queen's case spatial adjacency matrix (w ij = 1 when county i and j share a boundary or a vertice, w ij = 0 otherwise) that is required as input for the conditional autoregressive distribution was created using the Adjacency for WinBUGS Tool developed by Upper Midwest Environmental Sciences Center of the US Geological Survey (USGS) [46]. Typical disease mapping applications consider adjacency-based weights. Other weighting options (e.g., distance-decay weights) also appear in the literature (e.g., [47]) but are much less widely applied. Two Markov chains were simulated in the present study. The MCMC samplers were given initial values for each stochastic node. A total of 80,000 iterations with 60,000 burn-in were run. "Burn-in" denotes iterations that were discarded due to non-convergence of the Markov chains at the early stages of the simulation. A visual inspection of the Markov chain trace plots of the sample values versus iteration has indicated that convergence has been reached after 60,000 iterations. Inference was based on iterations 60,001 to 80,000.


Pearson's correlation

Figure 1 shows the 20 monitoring sites with top significant correlations (r > 0.8, p < 0.05). Eighteen of the twenty sites are located east of the -100° longitude line. The average correlation is 0.67 in the east and 0.22 in the west. The highest correlation (r = 0.96) is at Nashville, Tennessee. The interpolated correlation surface is shown in Figure 2. The surface indicates that the eastern part of the USA has higher correlation between PM2.5 and AOD.

Figure 1
figure 1

Monitoring sites with top 20 correlations of PM 2.5 with AOD.

Figure 2
figure 2

Surface of Pearson's correlation between PM 2.5 and AOD.


Results from the GWR show that all the monitoring sites have a condition number less than 30. The condition number evaluates local collinearity. In the presence of strong local collinearity, results become unstable. Results associated with condition numbers larger than 30 may be unreliable. Figure 3 shows a map of local R square. R2 values range between 0.0 and 1.0 and indicate how well the local regression model fits observed PM2.5 values. Very low values indicate the local model is performing poorly. It can be seen that GWR predicts well in the eastern USA and poorly in the west. Figure 4 shows the coefficient raster surface for the explanatory variable AOD. The map exhibits regional variation in the explanatory variable. The relationship between PM2.5 and AOD is not spatially consistent (stationary) across the conterminous states. Like Pearson's correlation surface, eastern USA shows higher AOD coefficient values, while values in the west are lower. Even negative values are found in part of the western region. The Pearson's correlation and GWR analyses revealed that it is appropriate to estimate a PM2.5 surface using AOD for disease modelling for the region to the east of the -100° longitude line. Application of the GWR model to the two-year mean AOD raster resulted in a continuous PM2.5 surface with a root mean square error (RMSE) of 1.67 μg/m3. Compared to the average PM2.5 value of 10.18 μg/m3 for all the monitors, the PM2.5 estimation achieved an accuracy of 84%.

Figure 3
figure 3

Local R square of geographically weighted regression.

Figure 4
figure 4

Coefficient raster surface for AOD from geographically weighted regression.

Bayesian model

The number of data points (counties) in the Bayesian modelling is 2,306. The actual number of counties to the east of the -100 degree longitude is 2,506. Some counties with unreliable and suppressed disease data were omitted from the analysis. Figure 5 displays a trace of the 80,000 updates of the fixed effects in each of model parameters in β = (β0, β1)'. The horizontal axis represents the simulation iteration number. The vertical axis represents simulated parameter values. Positive parameter β1 values indicate positive relationship between PM2.5 and CCHD mortality. The red trace is for one Markov chain, and the blue for the other. For each of the two parameters (the intercept and the PM2.5 effect), the simulation moves away from the initial seed values. After 60,000 updates, each simulation generates values fluctuating around within a consistent range of values representing the posterior distribution of each model parameter, indicating that convergence has been reached. Inferences were made about the parameters of the model using the simulated values on iterations 60,001 to 80,000. Table 1 provides the estimated posterior mean, median, and associated 95% credible set for each of the fixed effects. A 95% credible set defines an interval having a 0.95 posterior probability of containing the parameter of interest which is assumed to be a random variable in Bayesian statistics. The 2.5%, 50% (median) and 97.5% quantiles and posterior mean were calculated via an approximate algorithm [48]. Summaries of the postconvergence MCMC samples provide posterior inference for model parameters. For example, the sample mean of the postconvergence sampled values for a particular model parameter provides an estimate of the marginal posterior mean and a point estimate of the parameter itself. The interval defined by the 2.5th and 97.5th quantiles of the postconvergence sampled values for a model parameter provides a 95% interval estimate of the parameter. Standard deviations and Monte Carlo (MC) errors were calculated to assess the accuracy of the simulation. The MC error is an estimate of the difference between the mean of the sampled values and the true posterior mean. As a rule of thumb, the simulation should be run until the Monte Carlo error for each parameter of interest is less than about 5% of the sample standard deviation. The MC errors calculated from the iterations 60,001 to 80,000 for both the parameters are less than 5% of the corresponding standard deviations, suggesting an accurate posterior estimate for each parameter. In addition, Figure 6 provides kernel estimates of the corresponding posterior densities. The horizontal axis represents simulated parameter value. The vertical axis represents the density of each parameter value. The results in Figures 5 and 6 and Table show an association between PM2.5 pollution on CCHD. The 95% credible set covers positive values: = 0.802, CI = (0.386, 1.225). The positive boundary values of the CI indicate that areas with higher values of PM2.5 also show high rates of CCHD mortality.

Figure 5
figure 5

Trace plots of the 80,000 Markov Chain Monte Carlo (MCMC) updates. Simulation trace plots for the intercept and the effect of PM2.5 on chronic coronary heart disease from the Bayesian hierarchical model with a convolution prior. Horizontal axis represents iteration number and vertical axis represents simulated parameter value. The red trace is for one Markov chain, and the blue for the other.

Figure 6
figure 6

Kernel estimates of the posterior densities of the fixed effects in the Bayesian hierarchical model. Horizontal axis represents simulated parameter values and vertical axis represents the density of each value.

Table 1 Results of Bayesian hierarchical modeling


The Pearson's correlation and GWR analyses found that the relationship and between PM2.5 and MODIS AOD is not spatially consistent across the conterminous states. In the west AOD poorly correlates with PM2.5. The east region exhibits high positive correlations between PM2.5 and AOD where a PM2.5 surface can be estimated using AOD data for assessment of PM2.5 effect on disease. The correlation variation is likely due to the difference in terrain, AOD retrieval algorithm [29] and meteorological conditions. The algorithm to calculate AOD is based on dark surface pixels, contrast with surface, and assumptions in the model on the types of pollutants and the terrain [33]. The lower correlations in general can be found in the more arid parts as in the west of the USA where the higher surface reflectance reduces contrasts and where the model assumes more dust and smoke than urban/industrial model. The MODIS aerosol retrieval algorithm uses the urban/industrial aerosol model for eastern USA and other densely populated regions, such as Western Europe and eastern China [49]. Urban-industrial aerosols are mainly from fossil fuel combustion in populated industrial regions and dominated by fine particles. The positive relationship in the east could exist in other regions of the world where the aerosol retrieval algorithm uses the urban/industrial aerosol model. AOD measurements calculated from the sulphate-rich urban/industrial aerosol model may be more useful in predicting ground-based particulate matter levels. The inherent differences in the datasets may also explain some of the correlation variation [29]. MODIS AOD measures aerosol scattering in a total column from ground to satellite while EPA PM monitors measure a size-based subset of PM in the boundary layer in a specific location.

The Bayesian model found a spatial association between age-race-standardized mortality rate of chronic coronary heart disease and PM2.5. There is an excess risk of CCHD mortality in areas with high PM2.5 levels. Fine aerosol particles tend to penetrate into the gas-exchange regions of the lung, and ultrafine particles are able to pass through the lungs to enter the blood circulation and affect other organs [50]. A study indicates that PM2.5 leads to high plaque deposits in arteries, causing vascular inflammation and atherosclerosis – a hardening of the arteries that reduces elasticity, which can lead to heart attacks and other cardiovascular problems [51, 52].

This study has several strengths. First, the present study calculated PM2.5 surface by merging satellite derived AOD with ground measurements using GWR. To our knowledge, no applications of satellite remote sensing data in spatial modeling of regional PM2.5 concentrations have considered spatially varying relationships. Existing studies often use a uniform linear relationship to estimate PM2.5 using AOD. While the uniform relationship could apply to a smaller region (such as a metropolitan area), it might not work for a larger region (such as the eastern USA including 37 states in the study). Second, ecological studies are particularly useful when disease data at individual level is not available and individual level of exposure is either difficult or impossible to obtain, or can be only measured imprecisely [47]. Ecological studies are more useful for generating and testing hypothesis [53]. The association between CCHD mortality rate and PM2.5 revealed by the Bayesian model can be taken as indicative (though not necessarily causative) of a potential air pollution effect. Third, the MCMC computation techniques and Bayesian analysis allow us to use more elaborate models without conventional model restrictions. The model also allows us to treat the CCHD death counts as outcomes without the awkward transformations of the outcome and covariates of interest required to meet standard assumptions for linear regression models. Fourth, the MODIS aerosol data lends itself to population-based exposure assessment and the ecological approach. Moderate resolution satellite image pixels share the same feature dimension with the disease data enumeration districts in that both pixels and county polygons are area features. Derivation of a continuous PM2.5 surface using AOD alleviates the non- straightforwardness and noncompliance of comparing discrete point monitoring data with area-aggregated disease data. A continuous PM2.5 surface calculated by merging AOD and ground data accounting for varying correlations better represents aggregate exposure than averaging ground measurements over a limited number of monitoring sites.

There are also several limitations to this study. First, the correlation and GWR analyses did not account for meteorological parameters (such as mixing height, relative humidity, air temperature, and wind speed), aerosol vertical distribution, and aerosol properties. Although AOD is a highly significant predictor of PM2.5, those parameters could also modify the relationship between AOD and PM2.5 concentrations. Second, the study estimated an ambient PM2.5 surface to assess the health effect of fine particulate air pollution. However, ambient concentrations do not necessarily represent actual exposures, which can be influenced by the infiltration of ambient concentrations into indoor facilities (such as automobiles, homes, schools, and work places) and the activity of individuals (such as outdoor exercise, walking, commuting, etc.). Third, the use of aggregated data and therefore inferences based on the analysis cannot be directly transferred to the individual level. An inherent limitation of an ecological study is that it uses aggregate data and does not have the ability to incorporate individual information. The population-based ecological analysis did not consider population dynamics, for example, population migration. Fourth, chronic heart disease is a long term medical condition that has been shown to progress at varying rates over a lifetime with many other environmental, biomedical and social factors in effect. The study examined chronic heart disease over a short timeframe. Last, it must be noted that there are limitations of AOD data. AOD is a quantitative measure of total column aerosol, which is the mass of aerosols within a measured column extending from Earth to the satellite sensor. AOD does not correlate well with PM2.5 during cold seasons. This study used warm season AOD data but the disease data covered whole seasons. AOD does not differentiate between fine-mode and course-mode particles. Ideally, fine-mode optical depth should be used for the health effect assessment because fine mode particles which dominate urban/industrial pollution cause the most severe health problem. However, land-based measures of fine-mode fraction of AOD can only be used as a qualitative indicator of whether AOD values are dominated by natural or anthropogenic emissions [33]. Additionally, AOD does not specify the location of aerosols within a column and the AOD values do not necessarily represent ground conditions. This lack of vertical information emphasizes the importance of combining the satellite image with vertical profiles. Emerging Light Detection and Ranging (LIDAR) systems could provide vertical resolution for AOD and quantify pollutant levels on the ground.


There is a spatial variation of the relationship between PM2.5 and AOD in the conterminous USA. In the eastern USA where AOD correlates well with PM2.5, AOD can be merged with ground PM2.5 data to derive a PM2.5 surface for epidemiological study. The study found that chronic coronary heart disease mortality rate increases with exposure to PM2.5. High risk of CCHD mortality was found in areas with elevated levels of fine particulate air pollution. The association between CCHD mortality and PM2.5 would justify the need of further toxicological approach for investigating the biological mechanism of the adverse effect of fine particulate air pollution on heart. The evidence of raised CCHD mortality risk in high pollution areas would support targeting of policy interventions on such areas to reduce pollution levels. Aerosol remote sensing could help fill pervasive data gaps that impede efforts to study air pollution and protect public health.


  1. Pope CA: Epidemiology of fine particulate air pollution and human health: biologic mechanisms and who's at risk?. Environ Health Perspect. 2000, 108 (Supl 4): 713-723. 10.2307/3454408.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  2. Peters A, Dockery DW, Muller JE, Mittleman MA: Increased particulate air pollution and the triggering of myocardial infarction. Circulation. 2001, 103: 2810-2815.

    Article  PubMed  CAS  Google Scholar 

  3. Anderson HR, Atkinson RW, Bremner SA, Marston L: Particulate air pollution and hospital admissions for cardiorespiratory diseases: are the elderly at greater risk?. Eur Respir J Suppl. 2003, 40 (Suppl): 39s-46s. 10.1183/09031936.03.00402203.

    Article  PubMed  CAS  Google Scholar 

  4. Zanobetti A, Schwartz J: The effect of particulate air pollution on emergency admissions for myocardial infarction: a multicity case-crossover analysis. Environ Health Perspect. 113 (8): 978-982.

  5. Hu Z, Liebens J, Rao KR: Linking stroke mortality with air pollution, income, and greenness in northwest Florida: an ecological geographical study. Int J Health Geogr. 2008, 7: 20-

    Article  PubMed  PubMed Central  Google Scholar 

  6. Poloniecki JD, Atkinson RW, de Leon AP, Anderson HR: Daily time series for cardiovascular hospital admissions and previous day's air pollution in London, UK. Occup Environ Med. 1997, 54: 535-540. 10.1136/oem.54.8.535.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Schwartz J: Harvesting and long-term exposure effects in the relationship between air pollution and mortality. Am J Epidemiol. 2000, 151: 440-448.

    Article  PubMed  CAS  Google Scholar 

  8. Le Tertre A, Medina S, Samoli E, Forsberg B, Michelozzi P, Boumghar A, Vonk JM, Bellini A, Atkinson R, Ayres JG, Sunyer J, Schwartz J, Katsouyanni K: Short-term effects of particulate air pollution on cardiovascular diseases in eight European cities. J Epidemiol Community Health. 2002, 56 (10): 773-779. 10.1136/jech.56.10.773.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  9. Maheswaran R, Haining RP, Brindley P, Law J, Pearson T, Fryers PR, Wise S, Campbell MJ: Outdoor air pollution, mortality, and hospital admissions from coronary heart disease in Sheffield, UK: a small-area level ecological study. Eur Heart J. 2005, 26: 2543-2549. 10.1093/eurheartj/ehi457.

    Article  PubMed  CAS  Google Scholar 

  10. von Klot S, Peters A, Aalto PR, Hörmann A, Kulmala M, Lanki T, Löwel H, Pekkanen J, Picciotto S, Sunyer J, Forastiere F: Ambient air pollution is associated with increased risk of hospital cardiac readmissions of myocardial infarction survivors in five European cities. Circulation. 2005, 112: 3073-3079. 10.1161/CIRCULATIONAHA.105.548743.

    Article  PubMed  CAS  Google Scholar 

  11. Murakami Y, Ono M: Myocardial infarction deaths after high level exposure to particulate matter. J Epidemiol Community Health. 2006, 60: 262-266. 10.1136/jech.2005.038315.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Mille KA, Siscovick DS, Sheppard L, Shepher K, Sullivan JH, Anderson GL, Kaufman JD: Long-Term Exposure to Air Pollution and Incidence of Cardiovascular Events in Women. N Engl J Med. 2007, 356 (5): 447-458. 10.1056/NEJMoa054409.

    Article  Google Scholar 

  13. Chu DA, Kaufman YJ, Ichoku C, Remer LA, Tanre' D, Holben BN: Validation of the MODIS aerosol optical depth retrieval over land. Geo Res Lett. 2002, 29-

    Google Scholar 

  14. Wang J, Christopher SA: Intercomparison between satellite-derived aerosol optical thickness and PM2.5 mass: implications for air quality studies. Geo Res Lett. 2003, 30: 2095-10.1029/2003GL018174.

    Article  Google Scholar 

  15. Hutchison KD: Applications of MODIS satellite data and products for monitoring air quality in the state of Texas. Atmos Environ. 2003, 7: 32403-2412.

    Google Scholar 

  16. Hutchison KD, Smith S, Faruqui SJ: Correlating MODIS aerosol optical thickness data with ground-based PM2.5 observations across Texas for use in a real-time air quality prediction system. Atmos Environ. 2005, 39: 7190-7203. 10.1016/j.atmosenv.2005.08.036.

    Article  CAS  Google Scholar 

  17. Hu Z, Liebens J, Rao KR: Assessing health effects of fine aerosol particles using MODIS atmospheric data. Proceedings of the 31st International Geographical Congress. 2008, Tunis, Tunisia

    Google Scholar 

  18. Kumar N, Chu A, Foster A: An empirical relationship between PM2.5 and aerosol optical depth in Delhi Metropolitan. Atmos Environ. 2007, 41: 4492-4503. 10.1016/j.atmosenv.2007.01.046.

    Article  CAS  PubMed Central  Google Scholar 

  19. Schafer K, Harbusch A, Emeis S, Koepke P, Wiegner M: Correlation of aerosol mass near the ground with aerosol optical depth during two seasons in Munich. Atmos Environ. 2008, 42: 4036-4046. 10.1016/j.atmosenv.2008.01.060.

    Article  Google Scholar 

  20. Vidot J, Santer R, Ramon D: Atmospheric particulate matter (PM) estimation from SeaWiFS imagery. Remote Sens Environ. 2007, 111: 1-10. 10.1016/j.rse.2007.03.009.

    Article  Google Scholar 

  21. Liu Y, Sarnat J, Kilaru V, Jacob D, Koutrakis P: Estimating ground-level PM2.5 in the eastern United States using satellite remote sensing. Environ Sci Technol. 2005, 39 (9): 3269-3278. 10.1021/es049352m.

    Article  PubMed  CAS  Google Scholar 

  22. Liu Y, Franklin M, Kahn R, Koutrakis P: Using aerosol optical thickness to predict ground-level PM2.5 concentrations in the St. Louis area: a comparison between MISR and MODIS. Remote Sens Environ. 2007, 107 (1–2): 33-44. 10.1016/j.rse.2006.05.022.

    Article  Google Scholar 

  23. van Donkelaar A, Martin RV, Park RJ: Estimating ground-level PM2.5 using aerosol optical depth determined from satellite remote sensing. J Geophys Res. 2006, 111 (D21): 10.1029/2005JD006996.

  24. Seinfeld JH, Pandis SN: Atmospheric chemistry and physics: from air pollution to global change. 1998, New York: John Wiley & Sons

    Google Scholar 

  25. Kahn R, Banerjee P, McDonald D, Diner D: Sensitivity of multiangle imaging to aerosol optical depth and to pure-particle size distribution and composition over ocean. J Geophys Res. 1998, 103 (D24): 32195-32213. 10.1029/98JD01752.

    Article  Google Scholar 

  26. Liu Y, Paciorek CJ, Koutrakis P: Estimating regional spatial and temporal variability of PM2.5 concentrations using satellite data, meteorology, and land use information. Environ Health Perspec.

  27. HELIX-Atlanta (Health & Environment Linked for information exchange, Atlanta).

  28. Apituley MS, Koelemeijer R, Timmermans R, Schoemaker R, de Leeuw G: Construction of satellite derived PM2.5 maps using the relationship between AOD and PM2.5 at the Cabauw experimental site for atmospheric research (CESAR) – the Netherlands. IEEE International Geoscience & Remote Sensing Symposium: 6–11 July 2008. 2008, Boston, Massachusetts, USA

    Google Scholar 

  29. Engel-Cox JA, Holloman CH, Coutant BW, Hoff RM: Qualitative and quantitative evaluation of MODIS satellite sensor data for regional and urban scale air quality. Atmos Environ. 2004, 38: 2495-2509. 10.1016/j.atmosenv.2004.01.039.

    Article  CAS  Google Scholar 

  30. Chu DA, Kaufman YJ, Zibordi G, Chern JD, Mao J, Li C, Holben BN: Global monitoring of air pollution over land from EOS-Terra MODIS. J Geophys Res. 2003, 108 (D21): 4661-10.1029/2002JD003179.

    Article  Google Scholar 

  31. Al-Saadi J, Szykman J, Pierce RB, Kittaka C, Neil D, Chu DA, Remer L, Gumley L, Prins E, Weinstock L, MacDonald C, Wayland R, Dimmick F, Fishman J: Improving national air quality forecasts with satellite aerosol observations. Bull Amer Met Soc. 2005, 86: 1249-1261. 10.1175/BAMS-86-9-1249.

    Article  Google Scholar 

  32. Gupta P, Christopher SA, Wang J, Gehrig R, Lee Y, Kumar N: Satellite remote sensing of particulate matter and air quality assessment over global cities. Atmos Environ. 2006, 40: 5880-5892. 10.1016/j.atmosenv.2006.03.016.

    Article  CAS  Google Scholar 

  33. Algorithm for remote sensing of tropospheric aerosol from MODIS: Collection 005.

  34. Ackerman SA, Strabala KI, Menzel WP, Frey RA, Moeller CC, Bumley LE: Discriminating clear sky from clouds with MODIS. J Geophys Res. 1998, 103: 139-140. 10.1029/1998JD200032.

    Google Scholar 

  35. NASA Level 1 and Atmosphere Archive and Distribution System (LAADS Web).

  36. Integrated Earth Observations: Application to Air Quality and Human Health.

  37. US EPA Air Quality System (AQS) online Data Mart.

  38. Fotheringham SA, Chris B, Martin C: Geographically Weighted Regression: the Analysis of Spatially Varying Relationships. 2002, John Wiley & Sons

    Google Scholar 

  39. Centers for Disease Control and Prevention, National Center for Health Statistics.

  40. Mausner RB, Kramer S: Epidemiology: An Introductory Text. 1985, Philadelphia: W. B. Saunders

    Google Scholar 

  41. Waller LA, Gotway CA: Applied Spatial Statistics for Public Health Data. 2004, John Wiley & Sons, Inc

    Book  Google Scholar 

  42. Curtin LR, Klein RJ: Direct Standardization (Age-Adjusted Death Rates). Statistical Notes; No. 6. 1995, Hyattsville, Maryland: National Center for Health Statistics

    Google Scholar 

  43. Besag J, York Y, Mollie A: Bayesian image restoration with two applications in spatial statistics (with discussion). Ann I Stat Math. 1991, 43: 1-59. 10.1007/BF00116466.

    Article  Google Scholar 

  44. Mollie A: Bayesian mapping of disease. Markov Chain Monte Carlo in Practice. Edited by: Gilks WR, Richardson S, Spiegelhalter DJ. 1986, New York: Chapman and Hall, 359-379.

    Google Scholar 

  45. Lunn DJ, Thomas A, Best N, Spiegelhalter D: WinBUGS – a Bayesian modelling framework: concepts, structure, and extensibility. Stat Comput. 2000, 10: 325-337. 10.1023/A:1008929526011.

    Article  Google Scholar 

  46. Adjacency For WinBUGS Tool.

  47. Best N: Bayesian ecological modeling. Disease Mapping and Risk Assessment for Public Health. Edited by: Lawson A, Biggeri A, Böhning D, Lesaffre E, Viel JF, Bertollini R. 1999, West Sussex, England: John Wiley & Sons Ltd, 193-201.

    Google Scholar 

  48. Tierney L: A space-efficient recursive procedure for estimating a quantile of an unknown distribution. SIAM J Sci Stat Comp. 1993, 4: 706-711.

    Article  Google Scholar 

  49. Remer LA, Kaufman YJ, Tanre D, Mattoo S, Chu DA, Martins JV, Li RR, Ichoku C, Levy RC, Kleidman RG, Eck TF, Vermote E, Holben BN: The MODIS aerosol algorithm, products, and validation. J Atmos Sci. 2005, 62: 947-972. 10.1175/JAS3385.1.

    Article  Google Scholar 

  50. Nemmar A, Hoet PH, Vanquickenborne B, Dinsdale D, Thomeer M, Hoylaerts MF, Vanbilloen H, Mortelmans L, Nemery B: Passage of inhaled particles into the blood circulation in humans. Circulation. 2002, 105: 411-414. 10.1161/hc0402.104118.

    Article  PubMed  CAS  Google Scholar 

  51. Pope CA, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, Thurston GD: Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. J Amer Med Assoc. 2002, 287: 1132-1141. 10.1001/jama.287.9.1132.

    Article  CAS  Google Scholar 

  52. Suwa T, Hogg JC, Quinlan KB, Ohgami A, Voncent R, van Eeden SF: Particulate air pollution induces progression of atherosclerosis. J AM Coll Cardiol. 2002, 39: 935-942. 10.1016/S0735-1097(02)01715-1.

    Article  PubMed  CAS  Google Scholar 

  53. Rytkönen M, Moltchanova E, Ranta J, Taskinen O, Tuomilehto J, Karvonen M: The incidence of type 1 diabetes among children in Finland – rural-urban difference. Health Place. 2003, 9: 315-325. 10.1016/S1353-8292(02)00064-3.

    Article  PubMed  Google Scholar 

Download references


This study is a component of the "Assessment of Environmental Pollution and Community Health in Northwest Florida" supported by U.S. EPA Cooperative Agreement Award X-9745002 to the University of West Florida. The content of this report are solely the responsibility of the author and do not necessarily represent the official views of the U.S. EPA. The author thanks Tianping Bi, a visiting scholar of the University of West Florida from Northeast University of China, for help with satellite image batch processing.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Zhiyong Hu.

Additional information

Competing interests

The author declares that they have no competing interests.

Authors' contributions

ZH conceived of the study, conducted the research and wrote the manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Hu, Z. Spatial analysis of MODIS aerosol optical depth, PM2.5, and chronic coronary heart disease. Int J Health Geogr 8, 27 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: