Comparison of Spatial Modelling Approaches on PM(10) and NO(2) Concentration Variations: A Case Study in Surabaya City, Indonesia

Because of fast-paced industrialization, urbanization, and population growth in Indonesia, there are serious health issues in the country resulting from air pollution. This study uses geospatial modelling technologies, namely land-use regression (LUR), geographically weighted regression (GWR), and geographic and temporal weighted regression (GTWR) models, to assess variations in particulate matter (PM(10)) and nitrogen dioxide (NO(2)) concentrations in Surabaya City, Indonesia. This is the first study to implement spatiotemporal variability of air pollution concentrations in Surabaya City, Indonesia. To develop the prediction models, air pollution data collected from seven monitoring stations from 2010 to 2018 were used as dependent variables, while land-use/land cover allocations within a 250 m to 5000 m circular buffer range surrounding the monitoring stations were collected as independent variables. A supervised stepwise variable selection procedure was applied to identify the important predictor variables for developing the LUR, GWR, and GTWR models. The developed models of LUR, GWR, and GTWR accounted for 49%, 50%, and 51% of PM(10) variations and 46%, 47%, and 48% of NO(2) variations, respectively. The GTWR model performed better (R(2) = 0.51 for PM(10) and 0.48 for NO(2)) than the other two models (R(2) = 0.49–0.50 for PM(10) and 0.46–0.47 for NO(2)), LUR and GWR. In the PM(10) model four predictor variables, public facility, industry and warehousing, paddy field, and normalized difference vegetation index (NDVI), were selected during the variable selection procedure. Meanwhile, paddy field, residential area, rainfall, and temperature played important roles in explaining NO(2) variations. Because of biomass burning issues in South Asia, the paddy field, which has a positive correlation with PM(10) and NO(2), was selected as a predictor. By using long-term monitoring data to establish prediction models, this model may better depict PM(10) and NO(2) concentration variations within areas across Asia.