PM2.5 Pollutant in Asia-A Comparison of Metropolis Cities in Indonesia and Taiwan.

Air pollution has emerged as a significant health, environmental, economic, and social problem all over the world. In this study, geospatial technologies coupled with a LUR (Land Use Regression) approach were applied to assess the spatial-temporal distribution of fine particulate (PM2.5). In-situ observations of air pollutants from ground monitoring stations from 2016-2018 were used as dependent variables, while the land-use/land cover, a NDVI (Normalized Difference Vegetation Index) from a MODIS sensors, and meteorology data allocations surrounding the monitoring stations from 0.25-5 km buffer ranges were collected as spatial predictors from GIS and remote sensing databases. A linear regression method was developed for the LUR model and 10-fold cross-validation was used to assess the model robustness. The R2 model obtained was 56% for DKI Jakarta, Indonesia, and 83% for Taipei Metropolis, Taiwan. According to the results of the PM2.5 model, the essential predictors for DKI Jakarta were influenced by temperature, NDVI, humidity, and residential area, while those for the Taipei Metropolis region were influenced by PM10, NO2, SO2, UV, rainfall, spring, main road, railroad, airport, proximity to airports, mining areas, and NDVI. The validation of the results of the estimated PM2.5 distribution use 10-cross validation with indicated R2 values of 0.62 for DKI Jakarta and 0.84 for Taipei Metropolis. The results of cross-validation show the strength of the model.

1. Introduction

Ambient air pollution has been related to increased levels of mortality and morbidity in megacities [1]. It is closely related to people’s daily lives. The relatively high growth of socio-economic, industrial, and urbanization activities in urban and suburban areas has great potential for increasing energy consumption, which is one of the factors of air pollution. There are several types of outdoor air pollution, such as Cox, NOx, Sox, Ox, and one of the pollutants that is quite dangerous for human health is PM2.5. PM2.5, usually called fine particulate matter, is an air pollutant consisting of a mixture of solid and liquid particles. PM2.5 has a diameter of less than 2.5 , and thus, it is often referred to as smooth PM, which also consists of ultrafine particles with a diameter of less than 0.1 [2]. Epidemiological research associated with the study of PM2.5 exposure to the incidence of lung cancer, respiratory problems, and early death has shown a positive association [3,4]. The results show that it harms human health.

Several studies have shown that the concentrations of PM (particulate matter) can be affected by geographical features such as land use information, meteorology, and satellite data, including temporal variation parameters [5]. Methods for evaluating local urban variability to fine particulate matter were needed for these studies. Various models have been established around the world to explore the statistical correlations between ground monitoring stations of PM2.5 and the variables derived from geographical data information due to develops in the performance of GIS (Geographical Information System) technology.

Research is being conducted on the impact of geographical parameters on improvement air quality has been carried out. However, it is a big challenge to access PM2.5 data, particularly in developing countries. During the last period, some methods have been established to tackle challenges related to air pollution, such as interpolation using kriging and IDW (Inverse Distance Weighting) [6,7], and LUR (Land Use Regression) [8,9,10]. The interpolation of pollutant concentrations is based on the monitoring sites with densely clustered stations, whereas it is difficult to monitor locations with few stations. LUR models have proven to be relevant for these approaches in recent years. The method is to develop statistical regression models based on GIS platforms. These can be used to estimate air pollutant levels in a particular site by establishing a statistical correlation between pollutant observations and potential prediction variables [11,12,13].

In major cities, the types of land use can affect the level of PM through urban area development [14]. The land changes of forest, grasslands, agriculture in residential areas, industrial sites, and commercial centers frequently lead to increased emission levels. Air pollution concentration is related to changes in meteorological conditions. The development of prediction method starts by analyzing periods of serious atmospheric pollution, which are correlated with meteorological conditions monitored during those periods. These elements are considered as predictors [15]. Meteorological conditions such as wind speed, wind direction, and rainfall have an impact on particulate matter exposure [16]. Furthermore, satellite data using remote sensing, such as identification of greenness information, are efficient for measuring a wide range with multitemporal variations and easily available. The Normalized Difference Vegetation Index (NDVI) from Moderate Resolution Imaging Spectroradiometer (MODIS) sensors have been commonly used in greenness inventory management, which is a satellite-based method that has established a strong dynamic range and responsiveness for recording and calculating spatial-temporal variations in vegetation density [17,18]. NDVI has been used to determine green space area, which has a negative correlation with particulate matter [19].

Indonesia, with an area 1,922,570 km² (land only), has five ground stations monitoring PM2.5 concentration, one of the station locations is in DKI Jakarta [20]. In contrast, Taiwan, with an area of 36,000 km², has 78 ground stations to monitor air quality concentrations. For controlling and preventing air pollution from exceeding air quality standards, monitoring mass concentrations and multi-element identification of PM2.5 can be used for particle characterization and estimation of pollutant sources. In general, air quality monitoring, especially for PM2.5 is performed in various major cities in each country.

This study compares the air quality DKI Jakarta, in Indonesia and the Taipei Metropolis, in Taiwan using LUR model development for land use, meteorology, and greenness related to PM2.5. All the parameters have multi-temporal variability in order to achieve better performance in the estimation of PM2.5. The result acquired from this research will be particularly useful when developing LUR models in each country, in epidemiological studies and environmental health research in each country.

2. Materials and Methods

2.1. Study Area

The research area was in DKI Jakarta, in Indonesia and the Taipei Metropolis, in Taiwan, which are both capital cities. Taipei Metropolis includes 41 districts, with a total area of 2324 km2. The population densities of Taipei City and New Taipei City were 9818 and 1947 persons/km2, respectively [21]. According to the statistical data from the Land Use Investigation of Taiwan, Taipei Metropolis is covered by 68.43% forest, 8.44% buildings, 7.05% agricultural land, and 5.07% transportation infrastructure [22]. DKI Jakarta covers an area 662.33 km2, with 44 districts. However, only the major cities in DKI Jakarta were taken as the study area for Indonesia. The distribution of population density in DKI Jakarta was 15,763 persons/km2 for South Jakarta City, 15,385 persons/km2 for East Jakarta City, 19,143 persons/km2 for Central Jakarta City, 19,516 persons/km2 for West Jakarta City, and 12,146 persons/km2 for North Jakarta City [23]. Based on the classification from satellite data around 69.91% is built up area, 16.48% is mixed tree vegetation, 7.42% is grass land, and 1.63% is paddy field [24].

2.2. PM2.5 Concentration Data

PM2.5 concentration data of PM2.5 in Jakarta, Indonesia, were obtained from the two monitoring stations through AirNow DOS, which collects PM2.5 concentration data from the U.S. Embassy in DKI Jakarta. In Taiwan, ground monitoring measurements of PM2.5 mass concentration in Taipei Metropolis were obtained from 17 automatic monitoring stations established by the Environmental Protection Administration Executive Yuan R.O.C (Taiwan) through the Taiwan Air Quality Monitoring Network. The stations distributed within the study area, including 14 general stations, two traffic stations, and one national park station, use beta ray analyzers. Daily concentration observations from 2016 to 2018 were aggregated into monthly averages. The illustrations of the study area and the monitoring stations are shown in Figure 1.

Figure 1

Distribution of PM2.5 monitoring stations, (a) in DKI Jakarta, Indonesia and (b) in Taipei Metropolis, Taiwan.

2.3. Geographic Information Datasets

Topography data for DKI Jakarta were provided by the Geospatial Information Agency. Daily meteorological data in DKI Jakarta were collected by Meteorological, Climatological, and Geophysical Agency (BMKG), which is a government agency. The data from monitoring stations on Java Island in Indonesia, including temperature, wind direction, wind speed, relative humidity, solar radiation, and rainfall, were obtained from the BMKG database center (Data Online—BMKG Database Center: http://dataonline.bmkg.go.id/home) from 1 January 2016 to 31 December 2018. Vegetation Indices (MOD13Q1) version 6 data are produced every 16 days at 250 × 250 m spatial resolution as a level 3 product by a Terra Moderate Resolution Imaging Spectroradiometer (MODIS). The algorithm chosen is based on the best pixel value from all acquisitions for each 16-day period, which include various criteria such as low clouds, low view angle, and the highest NDVI/EVI value [25]. Topography data for Taiwan can be accessed in open data provided by the government (https://data.gov.tw/en), the GIS-T Transportation Network Geographic Information Warehousing System (https://gist.motc.gov.tw/gist_web/GistMapGeneral/MapTopic), and the Geographic Information Map Cloud Service Platform (https://www.tgos.tw). Figure 2a describes the land use information in the big cities of DKI Jakarta, which is dominated by built up areas, such as residential areas, buildings, and industrial and trade areas. Figure 2b illustrates the land use information in Taipei Metropolis, which is dominated by vegetation, such as forest.

Figure 2

Land use map of (a) DKI Jakarta and (b) Taipei Metropolis.

The estimated values of the predictor variables of the monitoring site coordinates were calculated through the GIS platform. The geospatial variables for land use and NDVI data were compiled in the range of a circular radius from 250 to 5000 m at each PM2.5 monitoring station to describe the area around it.

Data collection needed in this research included land use data in each region related to the land use data format with the shapefile format, meteorological data, and NDVI data from the MODIS sensor. NDVI data from the MODIS sensor can be downloaded at (https://ladsweb.modaps.eosdis.nasa.gov/search/order/1); the data downloaded only includes the area related to this study. The data processing step was then carried out after the data was collected, to create a database of land use data, meteorological data, and NDVI data. Land use data was formed by raster calculation of the focal statistics radius of 250–5000 m using ArcGIS and Python, for each type of land use. Meteorological data from each station point were obtained using the Inverse Distance Weighting method. All the raster maps (50 × 50 m2) were developed for the PM2.5 model to create each element based on the predictors, such as land use types, meteorology conditions, and NDVI. The focal statistics function in ArcMap was used to describe the predictor data for each type of land use and NDVI. Circular buffers with radiuses of 250–5000 m around each PM2.5 monitoring stations were used to generate maps.

2.4. LUR Modelling and Validation

LUR models were established for PM2.5 based on ground station monitoring data from the monitoring network and predictor variables. In addition, R x64 3.5.2 software (RStudio, Auckland, New Zealand) was used for statistical analysis, model development, and validating the models. ArcGIS 10.3 and Python were used to identify geographical variables (e.g., land use, meteorology, and NDVI) and create the map for prediction of PM2.5. The PM2.5 concentration models for DKI Jakarta and Taipei Metropolis were developed based on monthly variance from 1 January 2016 to 31 December 2018.

Stepwise linear regression is a combination of forward and backward selection techniques. This type of regression is used to calculate the percentage explanation variability in order to optimize the LUR models. Explanation of R2 was used for optimizing the percentage variability and assessing LUR models. The univariate regression analysis was carried out with the R2 performance and the coefficient correlations were listed for all predictor variables were listed in order to observe the direction and significance of the association. The LUR model with the highest R2 and the correct direction according to the predetermined criteria, was considered as the preliminary model. The high rank variable was selected and then the model was run to find the next significant variable. The Variance Inflation Factor (VIF) were examined to identify multicollinearity in the model development. When performing stepwise linear regression, the given conditions had values of < 0.1 and VIF < 3 [19,26].

The PM2.5 prediction maps were generated for each cell with a resolution of 50 × 50 m using the regression Equation (1). The formula from the final result of the stepwise regression guidelines is shown as, follows: where, is the PM2.5 concentrations; is the constant intercept of regression equation; to are regression coefficients; and to are the potential predictors of the central point.

Cross-validation was used to test the strength of the model. This study conducted 10-fold cross-validation. The 10-fold cross-validation was used to test the performance of the LUR model. The validation method determined by 90% of the data was used as training data, and the remaining 10% was used for validation [19]. The models were chosen with high R2, adj R2, and low RMSE for monthly concentrations.

3. Results

3.1. PM2.5 Concentrations between Indonesia and Taiwan

The results of the temporal trend are shown in Figure 3. According to the figure, the DKI Jakarta region has a higher level of pollution during the dry season than the rainy season, whereas for the Taipei City area has a high level of pollution in the spring season and a low level of pollution in the summer. Air quality guidelines (AQGs) with the base target level of PM2.5 is specified by the World Health Organization. The AQGs of PM2.5 for short-term (24-h average) and long-term (annual average) exposure are 25 and 10 . Based on identified health effects, PM2.5 pollution index standard is used for long-term exposure. These are the lowest levels with a total increase in cardiopulmonary and lung cancer mortality of 95% to PM2.5 [27]. The guidelines on this matter have been established in Indonesia. The threshold of PM2.5 is 65 for the 24-h average and 15 for the annual average [28].

Figure 3

Comparison of time series trend of PM2.5 concentrations between DKI Jakarta and Taipei Metropolis from 2016 to 2018.

3.2. LUR Development

This study proposes two models of PM2.5 estimation, and the resulting model includes spatial-temporal PM2.5 estimation for DKI Jakarta, Indonesia, and Taipei Metropolis, Taiwan. The results of Spearman correlation were selected to look for variables that had high potential/correlation and had the same direction of correlation, for example, having a negative or positive correlation.

3.2.1. PM2.5 Model in Indonesia

According to the results of correlation analysis, all variables selected have an intuitive relationship with PM2.5, such as variables related to residential areas that have a positive correlation and have a negative correlation with NDVI/Greenness. This study used a total of 496 variables. The criteria are shown in Table 1, along with the estimated efficiency and strength of the developed model. The indicated R2 value is 56% for variations of PM2.5 in DKI Jakarta, Indonesia. The results show that six variables were statistically significant predictors. Five variables, including temperature, NDVI with a radius of 1500 m, NDVI with a radius of 1750 m, humidity, and residential areas have values of < 0.01, and NDVI with a radius of 4750 m has a value of < 0.05. VIF was applied to assess the collinearity of the predictors in the developed model. The proposed the LUR model has a VIF value of less than 3. It shows that there were no cases of multicollinearity between predictor variables.

Table 1

Coefficient estimates of the developed LUR (Land Use Regression) model for DKI Jakarta.

Variable	β	ρ	VIF	Partial R2	Model Performance
Intercept	−60.566	0.360			R2 = 0.56
Temperature	5.797	<0.001	1.57	0.325	ADJ R2 = 0.52
NDVI1500 m	−77.419	<0.001	1.07	0.098	RMSE = 8.19
NDVI4750 m	−54.467	<0.05	1.04	0.044
NDVI1750 m	−50.729	<0.01	1.20	0.27
Relative humidity	−0.779	<0.005	1.82	0.36
Residential Area4250 m	0.039	<0.05	1.39	0.26

Descriptive statistical results are shown in Table 2 for the DKI Jakarta. The results were used to find information from the selected variables considered in this study.

Table 2

Descriptive statistics of the selected predictors.

Variable (Unit)	Min	Max	Standard Deviation	Variance	Medium (25–75% Percentile)
a NDVI1500 m	1.43 × 10−1	3.36 × 10−1	5.27 × 10−2	2.79 × 10−3	0.256 (2.00 × 10−1–2.86 × 10−1)
b NDVI1750 m	1.08 × 10−1	3.43 × 10−1	5.99 × 10−2	3.59 × 10−3	0.272 (2.29 × 10−1–2.93 × 10−1)
c NDVI4750 m	1.35 × 10−1	2.97 × 10−1	4.22 × 10−2	1.79 × 10−3	0.236 (1.91 × 10−1–2.63 × 10−1)
Relative humidity	67.603	94.194	4.985	24.858	78.409 (75.741–80.981)
Temperature	30.357	33.872	7.43 × 10−1	0.5522	32.614 (32.085–32.897)
d Residential area4250 m	337.099	451.852	57.778	3,338.41	394.476 (337.099–451.852)

a Average NDVI within a radius of 1500 m; b Average NDVI within a radius of 1750 m; c Average NDVI within a radius of 4750 m; d Residential area within a radius of 4250 m.

3.2.2. PM2.5 model in Taiwan

Based on the result of Spearman correlation analysis, all selected variable have powerful correlations, such as PM10, NO2, SO2, fall and spring season, UV, rainfall, major roads with a radius 250 m, railways with radius of 4000 and 5000 m, airports with nearest distance and radius 2500 and 5000 m, quarrying sites with radius 5000 m, and NDVI with a radius 4000 m, have powerful correlation. The UV, rainfall, NDVI, fall season, and airport with nearest distance variables have a negative correlation to determine the quality of PM2.5. This study used a total of 69 variables. The criteria are shown in Table 3, along with the estimated efficiency and strength of the developed model. The indicated R2 value is 84% for PM2.5 variations in Taipei Metropolis. The results show 14 variables are statistically significant predictors, 13 variables including PM10, NO2, SO2, UV, spring, main road, railroad, airport, airport closest distance, mine area, NDVI had a value of < 0.01, and the level of rainfall has a value of = 0.07. The VIF value of the LUR model has a value of less than 3. This shows that there are no cases of multicollinearity between predictor variables.

Table 3

Coefficient estimates of the developed LUR model for Taipei Metropolis.

Variable	β	ρ	VIF	Partial R2	Model Performance
Intercept	1.978	0.16			R2 = 0.84
PM10	0.302	<0.001	2.37	0.56	ADJ R2 = 0.84
NO2	0.266	<0.001	2.98	0.09	RMSE = 2.57
SO2	2.101	<0.001	2.37	0.01
UV	−0.326	<0.001	2.39	0.01
Rainfall	−0.074	0.07	1.16	0.001
Fall	−0.501	0.05	1.21	0.002
Spring	1.828	<0.001	1.94	0.001
Major Road250 m	6.42 × 10−3	<0.001	1.38	0.01
Railway4000 m	1.151	<0.001	1.51	0.02
Railway5000 m	0.615	<0.001	1.47	0.02
Airport2500 m	0.092	<0.001	1.75	0.11
Airport5000 m	0.043	<0.001	2.94	0.01
Airportnearest distance	−1.03 × 10−4	<0.01	2.81	0.004
Quarryingsite5000 m	0.346	<0.001	1.29	0.003
NDVI4000 m	−4.46 × 10−4	<0.001	1.56	0.01

Descriptive statistical results are shown in Table 4 for Taipei Metropolis. The results were used to find information from the selected variables, which was used to develop the LUR model.

Table 4

Descriptive statistics of the selected predictors of Taiwan.

Variable (Unit)	Min	Max	Standard Deviation	Variance	Medium (25–75% Percentile)
PM10	7.1	58.3	8.79	77.2	33.6 (30–41.9)
NO2	2.8	34.1	4.6	21.2	19.2 (17.1–22.1)
SO2	1.73	4.11	0.457	0.21	2.75 (2.45–3.13)
UV	1.83	10.62	2.35	5.51	5.9 (3.76–8.21)
Rainfall	0.04	14.39	2.68	7.2	2.28 (1.38–3.89)
Fall	0	1	0.434	0.189	0 (0–1)
Spring	0	1	0.432	0.187	0 (0–0)
a Major Road250 m	0	586.4	144.8	20,989	30.9 (0–123.4)
b Railway4000 m	0	2.863	0.709	0.502	0.124 (0–0.622)
c Railway5000 m	0	6.29	1.97	3.9	0.318 (0–0.875)
d Airport2500 m	0	77.8	17.9	320.5	0 (0–0)
e Airport5000 m	0	54.7	19.7	387	0 (0–8.12)
f Airportnearest distance	1433	17,289	4431.3	19,636,747	7356 (4425.5–11,242)
g Quarryingsite5000 m	0	5.73	1.33	1.77	0 (0–0.318)
h NDVI4000 m	585	9985	1335.3	1,783,200	8826 (7985–NA)

a Major road within a radius of 250 m; b Railway within a radius of 4000 m; c Railway within a radius of 5000 m; d Airport within a radius of 2500 m; e Airport within a radius of 5000 m; f Airport within the nearest distance; g Quarry site within a radius of 5000 m; h Average NDVI within a radius of 4000 m.

3.3. Model Performance

In this study, PM2.5 monthly estimation results of the DKI Jakarta area shown in Figure 4 have a high average PM2.5 value. PM2.5 in the region fluctuated, decreased in 2017, and showed high values in 2016 and 2018. The energy consumption of urban residents and the volume of vehicles continuously have an impact on urban air quality. The Normalized Difference Vegetation Index shows that it can reduce the level of PM2.5 in urban areas through a negative correlation.

Figure 4

Prediction maps of spatial-temporal variability of PM2.5 concentration using the developed LUR model in DKI Jakarta: (a) 2016–2018, (b) 2016, (c) 2017, (d) 2018.

Estimation of PM2.5 in the Taipei Metropolis area, shown in Figure 5, has a high average PM2.5 value. PM2.5 in the region has decreased from year to year, from 2016–2018. Urban population activity, airport emission, and the volume of vehicles continuously have an impact on urban air quality.

Figure 5

Prediction maps of spatial-temporal variability of PM2.5 concentration using the developed LUR model in Taipei Metropolis: (a) 2016–2018, (b) 2016, (c) 2017, (d) 2018.

The range value (green to red) of PM2.5 was from 0–55.726 μg/m3. The value showed that density of PM2.5 particles was 55.726 μg/3, indicating the worst air quality. The resolution of the map shows the monthly averages from 2016–2018.

3.4. Model Validation

Figure 6 shows the comparison of R2, adj R2, and RMSE values. The results of R2 were 62% and 83% for DKI Jakarta and Taipei Metropolis, respectively. The 10-cross validation confirms the robustness of the PM2.5 model. Figure 6a shows that the average accuracy from 2016–2018 has an R2 value of 0.617 for DKI Jakarta. The validation shows that the model developed is entirely accurate. Figure 6b shows the comparison between model predictions and observations at each point in 2018. Figure 6 shows the comparison of predictions between DKI Jakarta and New Taipei City for PM2.5 in 2018.

Figure 6

Result of 10-fold cross-validation of (a) DKI Jakarta and (b) Taipei Metropolis.

4. Discussion

This study implemented land use, meteorology and greenness data and long-term PM2.5 monitoring data, from DKI Jakarta in Indonesia and Taipei Metropolis in Taiwan. The study covered the capital cities of the two countries. The variables were concerned with multi-temporal variability data. The LUR model has been suitable method for estimating the concentration of pollutants in several areas, especially for PM2.5 [9,11,12]. The model shows good predictive ability (R2 = 56% for DKI Jakarta, and 84% for Taipei Metropolis), with spatial resolution of 50 × 50 meters. LUR can be used to determine the relations between PM2.5 and various variables. The results were derived from the proposed two model development approaches. The model can be validated using 10-fold cross-validation which shows the model robustness. It can also be shown by the R2 value of cross validation (R2 = 56% for DKI Jakarta, and 84% for Taipei Metropolis). The results spatial-temporal maps can be obtained using the LUR model in each country.

A limitation of our study was the lack of traffic-related and satellite data (such as, AOD (Aerosol Optical Depth) from MODIS sensor) to determine the level of PM2.5 concentrations. Land use information in the two countries is based on different information standards, which could affect the number of variables used in selecting correlations. However, AOD (MCD19A2 products) from the MODIS sensor did not show pixel values, due to the cloud density/null value in the study area. It was difficult to obtain more specific GIS data for these sources.

However, this study showed that DKI Jakarta has a lower R2 value rather than Taiwan. This could be because to the total numbers of monitoring stations are not balanced between these regions. The small number of monitoring stations in DKI Jakarta could affect LUR model development due to the number of predictors factored, which might have influenced the performance of LUR model establishment. The varying distribution of station monitoring networks across countries could contribute to variations in PM2.5 level.

The proportion of residential areas, major roads, railways, airports, and quarrying sites has important effects across land use characteristics in the model of PM2.5. However, according to the land use data, residential areas showed a significant positive correlation with PM2.5 concentration in DKI Jakarta. Studies of PM, the correlation have demonstrated high correlation between PM2.5 level and residential areas, especially in high-density residential areas [29,30].

Several studies have demonstrated that meteorological factor (such as humidity, temperature, wind speed and wind direction) highly correlate with air pollution, espacially PM2.5 concentrations [31,32,33]. Statistical modelling of PM2.5 in previous studies consistently showed higher concentrations in spring periods compared to fall periods [34]. Coefficient estimation from the LUR model showed that PM2.5 has a positive correlation with spring periods and a negative correlation with fall periods.

NDVI was obtained from satellite data using MODIS sensors to calculate greenness in the research area. However, NDVI was not differentiated among greenness types (such as planted area, green spaces, public parks, etc.). The difference of greenness type might be relevant to air pollutants. Previous studies have shown that greenness can affect air pollution as well [35].The annual averages of PM2.5 in Taipei Metropolis and DKI Jakarta exceed the national policy threshold standard in Indonesia and Taiwan, which is 15 . Meanwhile, the concentration of PM2.5 both countries also exceed the threshold standard of the WHO. Accurate and frequent ground monitoring stations of air quality concentrations are necessary for the capital cities in each country. They can be used to assess air quality, identify the most relevant potential sources, strengthen management to control air quality, and provide advice to policymakers, especially to the governments. By a adopting remote sensing and GIS techniques, a comprehensive strategy can be used to develop an air quality monitoring network, which allows for controlling air quality for the protection of human health [36].

5. Conclusions

This study analyzed recent trends of PM2.5 concentrations from 2016–2018 in two countries in Asia. LUR models were used to predict of PM2.5 level. The assessment of PM2.5 based on spatial-temporal was strongly influenced by land use, meteorological conditions, and MODIS NDVI. PM2.5 pollution in the Taipei Metropolis region, has a positive correlation with PM10, SO2, NO2, spring conditions, main roads, railways, airports, mining areas, and has a negative correlation with UV, rainfall, airports within close distances. PM2.5 pollution in DKI Jakarta, was strongly influenced by humidity, NDVI, temperature and residential areas. PM2.5 pollution in DKI Jakarta, has a positive correlation with residential areas, temperatures and has a negative correlation with NDVI, humidity. The R2 values of the resulting model were 0.84 and 0.56 for Taipei Metropolis and DKI Jakarta, respectively. Meanwhile, the 10-cross validation result shows R2 values of 0.83 and 0.61 for Taipei Metropolis and DKI Jakarta, respectively.