Spatiotemporal Exposure Assessment of PM 2.5 Concentration Using a Sensor-Based Air Monitoring System

: Sensor-based air monitoring instruments (SAMIs) can provide high-resolution air quality data by offering a detailed mapping of areas that air quality monitoring stations (AQMSs) cannot reach. This enhances the precision of estimating PM 2.5 concentration levels for areas that have not been directly measured, thereby enabling an accurate assessment of exposure. The study period was from 30 September to 2 October 2019 in the Guro-gu district, Seoul, Republic of Korea. Four models were applied to assess the suitability of the SAMIs and visualize the temporal and spatial distribution of PM 2.5 . Assuming that the PM 2.5 concentrations measured at a SAMI located in the center of the Guro-gu district represent the true values, the PM 2.5 concentrations estimated using QGIS spatial interpolation techniques were compared. The SAMIs were used at seven points (S 1 –S 7 ) according to the distance. Models 3 and 4 accurately estimated the unmeasured points with higher coefficients of determination (R 2 ) than the other models. As the distance from the AQMS increased from S 1 to S 7 , the R 2 between the observed and estimated values decreased from 0.89 to 0.29, respectively. The auxiliary installation of SAMIs could resolve regional concentration imbalances, allowing for the accurate estimation of pollutant concentrations and improved risk assessment for the population.


Introduction
Air pollution can have significant adverse effects on human health.The World Health Organization (WHO) has reported that nine out of 10 people worldwide breathe air containing high levels of pollutants, and approximately 7 million people die prematurely every year due to air pollution [1].Among these, fine particulate matter (PM 2.5 ) is the main risk factor [2].The International Agency for Research on Cancer has designated PM 2.5 as a Group 1 human carcinogen, which can cause cancer [3].PM 2.5 adversely affects health, causing respiratory diseases, such as asthma and chronic obstructive pulmonary disease, cardiovascular diseases, skin diseases, and allergies [4].Therefore, it is important to recognize the adverse effects of atmospheric PM 2.5 and develop strategies to manage exposure in the general population.
The Korean Ministry of Environment has operated an air quality monitoring station (AQMS) throughout the country since 1995 to determine whether national air pollution and air quality standards have been met.In addition, Air Korea (www.airkorea.or.kr (accessed on 1 April 2024)) provides exposure information to air pollutants by providing real-time air quality information measured by the AQMS to the public in the form of the air quality index through the Internet and smartphones.In previous epidemiological and risk assessment studies, the PM 2.5 concentration produced by AQMSs was used to compare mortality or human effects in the area or analyze correlations between population and prevalence of asthma, cardiovascular disease, and other diseases [5][6][7].These monitoring stations may provide accurate and reliable data on air quality in a specific area [8].However, the AQMSs cannot accurately assess personal or population exposure owing to spatial constraints, resulting in lower accuracy for air pollutants far from the monitoring station.Therefore, even though AQMS provides accurate data, it cannot explain important local variations crucial for public health protection in detail [9,10].In addition, there is a limitation to increasing the spatial density owing to the high cost of additional installations.
Currently, the trend to increase air quality data collection using a variety of methodologies other than AQMS has been growing globally.Sensor platforms for monitoring a variety of air pollutants are accessible, and new devices in this field are constantly being released [11].In recent years, many studies using sensor-based air monitoring instruments (SAMIs) have been conducted to address the spatial limitations of existing AQMSs [12][13][14].SAMIs can provide high-resolution air quality data by providing a more detailed mapping of areas that AQMSs cannot reach [15].
Owing to the nature of air pollutants, it is important to increase their spatiotemporal resolution because the influence of pollutants varies from region to region, and their distribution might not be uniform.This improvement in spatiotemporal resolution can increase the accuracy of concentration estimations for unmeasured points in modeling methods, such as Geographic Information Systems (GIS).
The objectives of this study were to enhance spatiotemporal resolution by complementing the spatial constraints of AQMSs through the auxiliary installation of SAMIs and to utilize this setup to improve the accuracy of PM 2.5 concentration estimation for unmeasured points.This enhances the precision of estimating PM 2.5 concentration levels for areas that have not been directly measured, thus enabling accurate assessments of exposure.This manuscript was organized to estimate the spatiotemporal concentrations using SAMIs and QGIS and to provide an optimal installation distance through a comparison of applied models.

Target Area and Locations of SAMIs and AQMSs
The study area was the Guro-gu district, Seoul, Republic of Korea (Figure 1).It had a population of 444,146 in January 2020, accounting for 4.2% of the population of Seoul, with an area of 20.12 km 2 .Four models were applied using AQMSs and SAMIs for the temporal and spatial visualization of

PM2.5 Concentration Measurement
Outdoor PM2.5 concentrations, measured at one AQMS in the Guro district and 11 nearby AQMSs, were collected from the data provided by Air Korea (www.airkorea.or.kr (accessed on 30 September 2019)) of the Ministry of Environment, as shown in Figure 2. The PM2.5 concentrations measured by the SAMIs were collected every minute and transmitted to the Government Cloud (G-Cloud), which is a cloud computing service developed for the Korean government's public institution and launched in 2012 [16].The collection period for all data was real-time concentration.
The collection period of all data was a real-time concentration in 1 h units from 30 September to 2 October 2019.The principle of PM2.5 concentration measurement by SAMI

PM 2.5 Concentration Measurement
Outdoor PM 2.5 concentrations, measured at one AQMS in the Guro district and 11 nearby AQMSs, were collected from the data provided by Air Korea (www.airkorea.or.kr (accessed on 30 September 2019)) of the Ministry of Environment, as shown in Figure 2. The PM 2.5 concentrations measured by the SAMIs were collected every minute and transmitted to the Government Cloud (G-Cloud), which is a cloud computing service developed for the Korean government's public institution and launched in 2012 [16].The collection period for all data was real-time concentration.
The collection period of all data was a real-time concentration in 1 h units from 30 September to 2 October 2019.The principle of PM 2.5 concentration measurement by SAMI was the light scattering spectrometer, and the specification was a sensor model of SD-16C (Dust Mon, Sentry Co., Ltd., Seoul, Republic of Korea), a flow rate of 0.5 L/min, and a resolution of 0.1 µg/m 3 .The comprehensive specifications of the SAMI are shown in Figure 3.
This instrument in SAMIs was designed to maintain a temperature of 20 • C and a relative humidity of less than 70% with a pretreatment control device and an integrated system because of measurement errors that may occur due to various environmental factors, as shown in Figure 3.In our previous study, it was found that for the accuracy of the PM 2.5 concentration, the coefficient of determination (R 2 ) was 0.96, obtained through a co-location test comparing the values monitored by installing a sensor-based measuring instrument in the same location as the air pollution monitoring station [17].Therefore, among the SAMIs installed in Guro-gu, we considered S 0 , located at the central point, to be equivalent to the AQMS, as shown in Figure 2. Compared with AQMSs, SAMIs provided slightly lower observed values, and the data were applied without applying a separate correction factor.For all PM 2.5 concentration data, outliers were removed using the interquartile range to increase the reliability of the values [18].
strument in the same location as the air pollution monitoring station [17].Therefore, among the SAMIs installed in Guro-gu, we considered S0, located at the central point, to be equivalent to the AQMS, as shown in Figure 2. Compared with AQMSs, SAMIs provided slightly lower observed values, and the data were applied without applying a separate correction factor.For all PM2.5 concentration data, outliers were removed using the interquartile range to increase the reliability of the values [18].

Spatial Interpolation
An interpolation methodology using GIS was used to estimate the PM2.5 concentration at unmeasured points.GIS provides information, such as maps, pictures, and diagrams, through the location data of objects with geographical locations and attribute data, which include information on the characteristics of those locations.For spatial modeling, the outdoor daily average PM2.5 concentration was calculated and used as attribute data, and the satellite coordinate system (WGS-84) of the location where the SAMI was installed in each region was used as the location data.Three spatial interpolation techniques were used: Inverse Distance Weighted (IDW), Ordinary Kriging (OK), and Universal Kriging (UK).The IDW technique was interpolated and visualized using the open-source program QGIS (version 3.22), and the kriging was interpolated for 1 h using the open-source program SAGA GIS (version 7.8.2) and visualized using QGIS.

Model Evaluation via Spatial Analysis
One station (S0) was selected from among the measurement stations located at the center of the Guro-gu district, as shown in Figure 2. Based on the AQMS points located in the Guro-gu district, seven stations (S1-S7) were selected based on the distance to the SAMI.The distance from station 1 to station 7 ranges from 1.4 km to 6.3 km, as shown in Figure 4. Verification was compared to the estimated value of S0 using the four models mentioned above.Three statistical analyses were performed to confirm the accuracy:

Spatial Interpolation
An interpolation methodology using GIS was used to estimate the PM 2.5 concentration at unmeasured points.GIS provides information, such as maps, pictures, and diagrams, through the location data of objects with geographical locations and attribute data, which include information on the characteristics of those locations.For spatial modeling, the outdoor daily average PM 2.5 concentration was calculated and used as attribute data, and the satellite coordinate system (WGS-84) of the location where the SAMI was installed in each region was used as the location data.Three spatial interpolation techniques were used: Inverse Distance Weighted (IDW), Ordinary Kriging (OK), and Universal Kriging (UK).The IDW technique was interpolated and visualized using the open-source program QGIS (version 3.22), and the kriging was interpolated for 1 h using the open-source program SAGA GIS (version 7.8.2) and visualized using QGIS.

Model Evaluation via Spatial Analysis
One station (S 0 ) was selected from among the measurement stations located at the center of the Guro-gu district, as shown in Figure 2. Based on the AQMS points located in the Guro-gu district, seven stations (S 1 -S 7 ) were selected based on the distance to the SAMI.The distance from station 1 to station 7 ranges from 1.4 km to 6.3 km, as shown in Figure 4. Verification was compared to the estimated value of S 0 using the four models mentioned above.Three statistical analyses were performed to confirm the accuracy: coefficient of determination (R 2 ), root mean squared error (RMSE), and mean absolute error (MAE).

PM2.5 Concentrations according to the Models
The distribution of PM2.5 concentrations, measured according to each model, is shown in Figure 5.As a result of analyzing the concentration of PM2.5 observed on September 30, Model 1 showed the lowest value among the four models at 29.6 ± 10.9 µg/m 3 .Average PM2.5 concentrations in Models 3, 4, and 2 were 62.6 ± 33.6 µg/m 3 , 55.3 ± 27.7

PM 2.5 Concentrations According to the Models
The distribution of PM 2.5 concentrations, measured according to each model, is shown in Figure 5.As a result of analyzing the concentration of PM 2.5 observed on 30 September, Model 1 showed the lowest value among the four models at 29.6 ± 10.9 µg/m 3 .Average PM 2.5 concentrations in Models 3, 4, and 2 were 62.6 ± 33.6 µg/m 3 , 55.3 ± 27.7 µg/m 3 , and 48.9

Verification between Estimation and Observation
The accuracy of each model was evaluated via statistical verification by comparing the measured values obtained from the central monitoring station in the Guro district (S 0 ) with the values estimated via interpolation in each model.The interpolation results for each estimation method at 13:00 on 30 September are shown in Figure 6.Model 1 was interpolated at a single point, resulting in R 2 , RMSE, and MAE values of 0.66, 45.0 µg/m 3 , and 34.0 µg/m 3 , respectively.These values remained consistent across all the interpolation methods.The R 2 increased or remained similar as the number of measurement points increased using the IDW technique, whereas the RMSE and MAE decreased, indicating an improvement in the accuracy of the model.Similar to the IDW technique, the accuracy of the model improved as the number of measurement points increased, even when the OK technique was used for statistical validation.The R 2 of Model 3 was 0.95, which was similar to that of Model 4, which included the monitoring stations in the vicinity.Model 3 was the most accurate in predicting S 0 , with an RMSE of 13.6 µg/m 3 and an MAE of 8.9 µg/m 3 .The UK technique showed that the R 2 values for Model 1 and Model 2, considering AQMSs, were 0.66 and 0.63, respectively.The corresponding RMSE values were 45.0 and 34.9.Models 3 and 4, considering the SAMIs, had the same results as OK, as shown in Table 1.ilar to that of Model 4, which included the monitoring stations in the vicinity.Model 3 was the most accurate in predicting S0, with an RMSE of 13.6 µg/m 3 and an MAE of 8.9 µg/m 3 .The UK technique showed that the R 2 values for Model 1 and Model 2, considering AQMSs, were 0.66 and 0.63, respectively.The corresponding RMSE values were 45.0 and 34.9.Models 3 and 4, considering the SAMIs, had the same results as OK, as shown in Table 1.

Comparison by SAMIs According to Distance from the AQMSs
The correlations between SAMIs (S 1 -S 7 ) according to the distance from the AQMSs in the Guro-gu district through each spatial interpolation technique were compared using OK (Figure 7).S 1 and S 2 present the scatter plot of R 2 derived using the IDW and UK method, respectively.

Discussion
This study was conducted to confirm whether the auxiliary installation of a SAMI overcomes the limitations of the existing AQMSs and improves the accuracy of calculating local residents' exposure to PM2.5.After classifying the four models using AQMSs and SAMIs, the PM2.5 concentrations measured were collected at the measurement stations of each model.Spatial interpolation techniques, such as IDW and kriging, were performed through the QGIS program to visualize the PM2.5 concentration using the collected PM2.5 concentration, and the concentration for unmeasured points was estimated.The estimated values obtained through interpolation were statistically verified by the comparison with the observed values measured using SAMIs (S0).
The AQMS is limited in its ability to measure air quality beyond a certain distance, which can reduce the accuracy of exposure assessments.In addition, owing to cost constraints, the spatial resolution cannot be increased.Therefore, installing a SAMI as an auxiliary installation could overcome spatial limitations at a relatively low cost and address the limitations of an AQMS by collecting more data on PM2.5 concentrations.In this study, cross validation was performed for each interpolation technique to compare and evaluate the accuracy of the models.
In previous epidemiological studies, PM2.5 concentrations measured at the AQMSs were estimated through spatial interpolation using GIS and subsequently utilized for exposure assessment [19].Spatial interpolation is a technique that estimates the values of unobserved points using the values of observed points.It can be used to predict the concentration of air pollutants in a range that cannot be measured using an AQMS [20].The accuracy of this process depends on the precision, quantity, and distribution of known points [21].In a study conducted in Delhi, India, pollutants generated from vehicles were visualized using the kriging technique in GIS.Subsequently, respiratory and Among the AQMSs and SAMIs, the R 2 of the monitoring station (S 1 ) located at 1.4 km was the highest at 0.81-0.91,and Model 3 showed the highest R 2 at 0.91.Monitoring stations (S 2 -S 4 ) within a 2 km radius in Model 3, which includes a SAMI with a relatively large number of monitoring stations, recorded relatively high R 2 values of 0.80, 0.72, and 0.86.Although there was a difference in the R 2 values according to the classification of the monitoring station, the trend of the graph decreased in all four models as the distance from the AQMSs increased.

Discussion
This study was conducted to confirm whether the auxiliary installation of a SAMI overcomes the limitations of the existing AQMSs and improves the accuracy of calculating local residents' exposure to PM 2.5 .After classifying the four models using AQMSs and SAMIs, the PM 2.5 concentrations measured were collected at the measurement stations of each model.Spatial interpolation techniques, such as IDW and kriging, were performed through the QGIS program to visualize the PM 2.5 concentration using the collected PM 2.5 concentration, and the concentration for unmeasured points was estimated.The estimated values obtained through interpolation were statistically verified by the comparison with the observed values measured using SAMIs (S 0 ).
The AQMS is limited in its ability to measure air quality beyond a certain distance, which can reduce the accuracy of exposure assessments.In addition, owing to cost constraints, the spatial resolution cannot be increased.Therefore, installing a SAMI as an auxiliary installation could overcome spatial limitations at a relatively low cost and address the limitations of an AQMS by collecting more data on PM 2.5 concentrations.In this study, cross validation was performed for each interpolation technique to compare and evaluate the accuracy of the models.
In previous epidemiological studies, PM 2.5 concentrations measured at the AQMSs were estimated through spatial interpolation using GIS and subsequently utilized for exposure assessment [19].Spatial interpolation is a technique that estimates the values of unobserved points using the values of observed points.It can be used to predict the concentration of air pollutants in a range that cannot be measured using an AQMS [20].The accuracy of this process depends on the precision, quantity, and distribution of known points [21].In a study conducted in Delhi, India, pollutants generated from vehicles were visualized using the kriging technique in GIS.Subsequently, respiratory and cardiovascular deaths were estimated based on exposure to air pollutants [22].Spatial interpolation was used to evaluate the efficacy of monitoring sites located in Seoul, Republic of Korea, and to validate the accuracy of atmospheric PM 2.5 concentration measurements.These results indicate that additional monitoring stations are required for high-resolution monitoring [23].In a study utilizing wind-field data, the spatial interpolation results were significantly enhanced compared with conventional methods, demonstrating the effectiveness of incorporating meteorological factors into PM 2.5 concentration estimations [24].Similarly, research on regional PM 2.5 concentration estimation found that using optimized spatial models with radial basis function interpolation greatly improved prediction accuracy [25].These findings suggest the potential for advancing interpolation techniques through various innovative approaches.
Similar to previous studies, this study also demonstrates the potential of overcoming the limitations of the measurement imbalance between AQMS regions by predicting the concentration in unmeasured regions through spatial interpolation.The auxiliary installation of SAMIs has the potential to improve the accuracy of forecasting ambient PM 2.5 levels and reduce costs by overcoming the limitations of the AQMSs.Several studies have conducted comparative evaluations to assess the accuracy of both the AQMSs and SAMIs [26].
In a previous study conducted in December 2018 to ensure the accuracy of SAMI, we conducted a co-location test with AQMS at the site S 0 over approximately 334 h across 15 days, resulting in a correlation coefficient (R 2 ) of 0.96 [17].Although the SAMIs were only placed outdoors at public buildings by the Guro-gu district, the monitoring locations were carefully selected.The selection criteria included the ability to provide representative data, quality assurance, site security, power access, and specific building consideration.The SAMI network adhered to a neighborhood spatial scale, as specified by the U.S. EPA's five categories, with sensors spaced 0.5 to 4.0 km apart [27].To assess the accuracy of 12 measuring instruments equipped with low-cost sensors, the study compared the PM 2.5 concentrations measured under identical conditions to those obtained from an ambient air monitoring station operated using the Federally Equivalent Methods in Riverside, California, USA [28].The data collected from the six low-cost sensors showed a relatively high correlation, with R 2 values ranging from 0.7 to 0.75.Furthermore, during a field test conducted in Australia, a correlation analysis was performed by computing the daily average concentration of ambient PM 2.5 measured by an air quality monitoring station and a low-cost IoT sensor [29].The results indicated a strong correlation between R 2 and the range of 0.63-0.87,which suggests a relatively high level of performance.Similar to the results by Feenstra et al. and Tagle et al. [28,29], the R 2 value of the predictive model in this study, which included SAMIs, showed a strong correlation ranging from 0.88 to 0.95.
The accuracy of estimating PM 2.5 concentration improved when SAMIs were used in combination with the AQMS compared with the results obtained through single-point interpolation when using only the AQMS.By estimating the AQMS using each spatial interpolation technique with Model 1 as a single point, the R 2 value was found to be 0.66.The accuracy improved significantly when SAMIs were installed as an auxiliary, increasing from approximately 23% to 55%.Moreover, the RMSE and MAE values were relatively low, indicating an increase in the accuracy of the model estimation of PM 2.5 .One simulation study demonstrated that the accuracy of PM 2.5 concentration predictions decreased as the number of monitoring stations decreased from 22 to 12 [30].In a study that evaluated the performance of a model estimating PM 2.5 , the accuracy improved in cities with a higher number of monitoring stations compared with those with fewer monitoring stations.This suggests that the accuracy of concentration prediction improves as the number of monitoring stations increases, owing to the increase in spatial resolution density [31].Furthermore, SAMIs can overcome spatial limitations and gather data from many locations [32].This can improve the accuracy of predicting PM 2.5 concentrations and contribute to more precise assessments of exposure.
Seven stations (S 1 -S 7 ) of the SAMIs were selected based on their distances from the AQMS.The estimated values were obtained through spatial interpolation and then compared to the observed values using statistical verification techniques.All four models (Models 1, 2, 3, and 4) showed that the R 2 value between the observation and estimation decreased with the distance from the installed AQMS.According to the installation and operation guidelines for AQMSs, the distance between air quality monitoring stations should be 2 km to avoid redundancy [33].This suggests that the AQMS represents the PM 2.5 concentration for an area within a radius of approximately 2 km and is not indicative of the concentration beyond this distance.SAMIs are considered to have spatial representativeness when the R 2 value is above 0.80 within a radius of approximately 1.5 km.Therefore, a clear estimation of unmeasured points was possible by installing SAMI at approximately 3 km intervals.The installation requirements for monitoring stations vary depending on the presence or absence of air pollutant sources.However, for large cities, such as Seoul, the accuracy of PM 2.5 concentration estimation could improve if measurement stations are installed at 3 km intervals.By comparing the statistical values of Models 3 and 4, it was observed that the R 2 value of Model 4 decreased.These results are contrary to the research direction, which stated that the accuracy of the PM 2.5 concentration estimation would increase with an increase in the number of monitoring stations.Owing to the nature of the interpolation technique, it is sensitive to points in close proximity.As a result, AQMSs located outside the Guro-gu district did not have a significant impact on the improvement in PM 2.5 concentration estimation.As the distance from the AQMS increased, the R 2 value between the observed and estimated values of the selected SAMIs decreased.Combining the AQMSs and SAMIs in a region could enhance the precision of PM 2.5 concentration estimations for unmeasured locations.
This study identified the limitations of AQMS in accurately measuring local PM 2.5 concentrations due to their restricted spatial coverage.By integrating SAMI, we were able to effectively address these shortcomings, specifically through the auxiliary installation of SAMIs, which significantly enhanced the spatial resolution and accuracy of PM 2.5 concentration estimations.Furthermore, our research established the optimal deployment range for SAMIs, ensuring effective spatial representativeness.However, the accuracy of the statistical verification for SAMI may be compromised because of the short duration of only 3 days.In addition, the measured outdoor PM 2.5 results may be biased due to missing values in the concentration data, and the study area was limited to the Guro-gu district, Seoul.
Despite these limitations, this study provides significant insights into improving air quality monitoring and demonstrates the potential of using SAMIs to enhance the accuracy and spatial resolution of PM 2.5 measurements.In future studies, extending the data collection period of SAMI would be essential to ensure sufficient data availability.Additionally, it could be necessary to collect and compare data from both urban and rural areas, as well as to gather data during significant air quality events, like yellow dust storms.Although the sensors used have a 2-year warranty, research should be carried out to detect sensor malfunctions and resulting data errors during the measurement period.

Conclusions
This study classified four models based on the increase in the number of monitoring stations.Using QGIS, interpolation techniques, such as IDW and kriging, were used to visualize the concentration, and statistical verification between the estimated concentration value and the observed value was performed.Based on the IDW, OK, and UK techniques in Model 3, the R 2 values were 0.88, 0.95, and 0.95, respectively.These values indicate a higher correlation than the results of the interpolation obtained using Model 1.In addition, compared with Model 1, the RMSE and MAE showed low values, indicating that the PM 2.5 concentration estimation model had high accuracy.Following the assessment of the ambient PM 2.5 concentration estimation accuracy using spatial interpolation techniques, it was observed that the kriging technique demonstrated a higher accuracy than the IDW technique in each model.When statistical verification was performed on the monitoring stations according to distance, R 2 decreased as the distance from the AQMS increased.This suggests that the auxiliary installation of a SAMI is necessary.Consequently, the spatial limitations and cost aspects of the AQMS can be supplemented through the SAMIs.This increase in spatial resolution would improve the accuracy of PM 2.5 concentration estimation for unmeasured points, resolve the concentration distribution imbalance in a specific space, and efficiently manage air quality.It is possible to calculate the exact amount of PM 2.5 exposure due to the increase in the measurement density.Therefore, it can be used for accurate exposure and risk assessments of the entire population and subgroups.

13 Figure 1 .
Figure 1.Location of the study area.

Figure 1 .
Figure 1.Location of the study area.

Figure 1 .
Figure 1.Location of the study area.

Figure 2 .
Figure 2. Measurement points of AQMSs and SAMIs according to the four kinds of models.

Figure 2 .
Figure 2. Measurement points of AQMSs and SAMIs according to the four kinds of models.

Figure 3 .
Figure 3. Specification of the SAMI (sensor-based air monitoring instrument) used in this study.

Figure 3 .
Figure 3. Specification of the SAMI (sensor-based air monitoring instrument) used in this study.

Figure 4 .
Figure 4. Selection of SAMI (S 1 -S 7 ) to evaluate PM 2.5 concentration estimation accuracy according to the distance from AQMS.

Figure 5 .
Figure 5. Distribution of PM 2.5 concentrations according to each model.

Figure 6 .
Figure 6.Visualization of PM2.5 concentration estimated according to the four kinds of models.Figure 6. Visualization of PM 2.5 concentration estimated according to the four kinds of models.

Figure 6 .
Figure 6.Visualization of PM2.5 concentration estimated according to the four kinds of models.Figure 6. Visualization of PM 2.5 concentration estimated according to the four kinds of models.

Figure 7 .
Figure 7. Scatter plot of the coefficient of determination derived using the Ordinary Kriging (OK) method for S1-S7 according to distance from the AQMS.

Figure 7 .
Figure 7. Scatter plot of the coefficient of determination derived using the Ordinary Kriging (OK) method for S 1 -S 7 according to distance from the AQMS.

Table 1 .
Comparison of statistical verification (R 2 , RMSE, and MAE) for each interpolation technique for S 0 .