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.


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 PM 2.5 . PM 2.5, usually called fine particulate matter, is an air pollutant consisting of a mixture of solid and liquid particles. PM 2.5 has a diameter of less than 2.5 µm, and thus, it is often referred to as smooth PM, which also consists of ultrafine particles with a diameter of less than 0.1 µm [2]. Epidemiological research associated with the study of PM 2.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 PM 2.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 2 (land only), has five ground stations monitoring PM 2.5 concentration, one of the station locations is in DKI Jakarta [20]. In contrast, Taiwan, with an area of 36,000 km 2 , 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 PM 2.5 can be used for particle characterization and estimation of pollutant sources. In general, air quality monitoring, especially for PM 2.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 PM 2.5 . All the parameters have multi-temporal variability in order to achieve better performance in the estimation of PM 2.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.

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 km 2 . The population densities of Taipei City and New Taipei City were 9818 and 1947 persons/km 2 , 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 km 2 , 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/km 2 for South Jakarta City, 15,385 persons/km 2 for East Jakarta City, 19,143 persons/km 2 for Central Jakarta City, 19,516 persons/km 2 for West Jakarta City, and 12,146 persons/km 2 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].

PM 2.5 Concentration Data
PM 2.5 concentration data of PM 2.5 in Jakarta, Indonesia, were obtained from the two monitoring stations through AirNow DOS, which collects PM 2.5 concentration data from the U.S. Embassy in DKI Jakarta. In Taiwan transportation infrastructure [22]. DKI Jakarta covers an area 662.33 km 2 , 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/km 2 for South Jakarta City, 15,385 persons/km 2 for East Jakarta City, 19,143 persons/km 2 for Central Jakarta City, 19,516 persons/km 2 for West Jakarta City, and 12,146 persons/km 2 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].

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

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

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 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. 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 m 2 ) 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. 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 PM 2.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 m 2 ) were developed for the PM 2.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 PM 2.5 monitoring stations were used to generate maps.

LUR Modelling and Validation
LUR models were established for PM 2.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 PM 2.5 . The PM 2.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 R 2 was used for optimizing the percentage variability and assessing LUR models. The univariate regression analysis was carried out with the R 2 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 R 2 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 PM 2.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, Y is the PM 2.5 concentrations; β 0 is the constant intercept of regression equation; β 1 to β n are regression coefficients; and X 1 to X n 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 R 2 , adj R 2 , and low RMSE for monthly concentrations.

PM 2.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 PM 2.5 is specified by the World Health Organization. The AQGs of PM 2.5 for short-term (24-h average) and long-term (annual average) exposure are 25 µg/m 3 and 10 µg/m 3 . Based on identified health effects, PM 2.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 PM 2.5 [27]. The guidelines on this matter have been established in Indonesia. The threshold of PM 2.5 is 65 µg/Nm 3 for the 24-h average and 15 µg/Nm 3 for the annual average [28].

LUR Development
This study proposes two models of PM2.5 estimation, and the resulting model includes spatialtemporal 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.

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 R 2 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. 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.

LUR Development
This study proposes two models of PM 2.5 estimation, and the resulting model includes spatial-temporal PM 2.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.

PM 2.5 Model in Indonesia
According to the results of correlation analysis, all variables selected have an intuitive relationship with PM 2.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 R 2 value is 56% for variations of PM 2.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. 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.

PM 2.5 model in Taiwan
Based on the result of Spearman correlation analysis, all selected variable have powerful correlations, such as PM 10 , NO 2 , SO 2 , 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 PM 2.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 R 2 value is 84% for PM 2.5 variations in Taipei Metropolis. The results show 14 variables are statistically significant predictors, 13 variables including PM 10 , NO 2 , SO 2 , 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. 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.

Model Performance
In this study, PM 2.5 monthly estimation results of the DKI Jakarta area shown in Figure 4

Model Performance
In this study, PM2.5 monthly estimation results of the DKI Jakarta area shown in Figure 4    Estimation of PM 2.5 in the Taipei Metropolis area, shown in Figure 5, has a high average PM 2.5 value. PM 2.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. The range value (green to red) of PM2.5 was from 0-55.726 μg/m 3 . 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. Figure 6 shows the comparison of R 2 , adj R 2 , and RMSE values. The results of R 2 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 R 2 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.

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 The range value (green to red) of PM 2.5 was from 0-55.726 µg/m 3 . The value showed that density of PM 2.5 particles was 55.726 µg/ 3 , indicating the worst air quality. The resolution of the map shows the monthly averages from 2016-2018. Figure 6 shows the comparison of R 2 , adj R 2 , and RMSE values. The results of R 2 were 62% and 83% for DKI Jakarta and Taipei Metropolis, respectively. The 10-cross validation confirms the robustness of the PM 2.5 model. Figure 6a shows that the average accuracy from 2016-2018 has an R 2 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 PM 2.5 in 2018. The range value (green to red) of PM2.5 was from 0-55.726 μg/m 3 . 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. Figure 6 shows the comparison of R 2 , adj R 2 , and RMSE values. The results of R 2 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 R 2 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.

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

Discussion
This study implemented land use, meteorology and greenness data and long-term PM 2.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 PM 2.5 [9,11,12]. The model shows good predictive ability (R 2 = 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 PM 2.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 R 2 value of cross validation (R 2 = 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 PM 2.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 R 2 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 PM 2.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 PM 2.5 concentration in DKI Jakarta. Studies of PM, the correlation have demonstrated high correlation between PM 2.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 PM 2.5 concentrations [31][32][33]. Statistical modelling of PM 2.5 in previous studies consistently showed higher concentrations in spring periods compared to fall periods [34]. Coefficient estimation from the LUR model showed that PM 2.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 PM 2.5 in Taipei Metropolis and DKI Jakarta exceed the national policy threshold standard in Indonesia and Taiwan, which is 15 µg/m 3 . Meanwhile, the concentration of PM 2.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].

Conclusions
This study analyzed recent trends of PM 2.5 concentrations from 2016-2018 in two countries in Asia. LUR models were used to predict of PM 2.5 level. The assessment of PM 2.5 based on spatial-temporal was strongly influenced by land use, meteorological conditions, and MODIS NDVI. PM 2.5 pollution in the Taipei Metropolis region, has a positive correlation with PM 10 , SO 2 , NO 2 , spring conditions, main roads, railways, airports, mining areas, and has a negative correlation with UV, rainfall, airports within close distances. PM 2.5 pollution in DKI Jakarta, was strongly influenced by humidity, NDVI, temperature and residential areas. PM 2.5 pollution in DKI Jakarta, has a positive correlation with residential areas, temperatures and has a negative correlation with NDVI, humidity. The R 2 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 R 2 values of 0.83 and 0.61 for Taipei Metropolis and DKI Jakarta, respectively.