Land-use regression with long-term satellite-based greenness index and culture-specific sources to model PM_2.5 spatial-temporal variability

This study utilized a long-term satellite-based vegetation index, and considered culture-specific emission sources (temples and Chinese restaurants) with Land-use Regression (LUR) modelling to estimate the spatial-temporal variability of PM_2.5 using data from Taipei metropolis, which exhibits typical Asian city characteristics. Annual average PM_2.5 concentrations from 2006 to 2012 of 17 air quality monitoring stations established by Environmental Protection Administration of Taiwan were used for model development. PM_2.5 measurements from 2013 were used for external data verification. Monthly Normalized Difference Vegetation Index (NDVI) images coupled with buffer analysis were used to assess the spatial-temporal variations of greenness surrounding the monitoring sites. The distribution of temples and Chinese restaurants were included to represent the emission contributions from incense and joss money burning, and gas cooking, respectively. Spearman correlation coefficient and stepwise regression were used for LUR model development, and 10-fold cross-validation and external data verification were applied to verify the model reliability. The results showed a strongly negative correlation (r: −0.71 to −0.77) between NDVI and PM_2.5 while temples (r: 0.52 to 0.66) and Chinese restaurants (r: 0.31 to 0.44) were positively correlated to PM_2.5 concentrations. With the adjusted model R² of 0.89, a cross-validated adj-R² of 0.90, and external validated R² of 0.83, the high explanatory power of the resultant model was confirmed. Moreover, the averaged NDVI within a 1750 m circular buffer (p < 0.01), the number of Chinese restaurants within a 1750 m buffer (p < 0.01), and the number of temples within a 750 m buffer (p = 0.06) were selected as important predictors during the stepwise selection procedures. According to the partial R², NDVI explained 66% of PM_2.5 variation and was the dominant variable in the developed model. We suggest future studies consider these three factors when establishing LUR models for estimating PM_2.5 in other Asian cities.