Articles | Volume 21, issue 6
https://doi.org/10.5194/acp-21-5063-2021
https://doi.org/10.5194/acp-21-5063-2021
Research article
 | 
31 Mar 2021
Research article |  | 31 Mar 2021

Development and intercity transferability of land-use regression models for predicting ambient PM10, PM2.5, NO2 and O3 concentrations in northern Taiwan

Zhiyuan Li, Kin-Fai Ho, Hsiao-Chi Chuang, and Steve Hung Lam Yim
Abstract

To provide long-term air pollutant exposure estimates for epidemiological studies, it is essential to test the feasibility of developing land-use regression (LUR) models using only routine air quality measurement data and to evaluate the transferability of LUR models between nearby cities. In this study, we developed and evaluated the intercity transferability of annual-average LUR models for ambient respirable suspended particulates (PM10), fine suspended particulates (PM2.5), nitrogen dioxide (NO2) and ozone (O3) in the Taipei–Keelung metropolitan area of northern Taiwan in 2019. Ambient PM10, PM2.5, NO2 and O3 measurements at 30 fixed-site stations were used as the dependent variables, and a total of 156 potential predictor variables in six categories (i.e., population density, road network, land-use type, normalized difference vegetation index, meteorology and elevation) were extracted using buffer spatial analysis. The LUR models were developed using the supervised forward linear regression approach. The LUR models for ambient PM10, PM2.5, NO2 and O3 achieved relatively high prediction performance, with R2 values of > 0.72 and leave-one-out cross-validation (LOOCV) R2 values of > 0.53. The intercity transferability of LUR models varied among the air pollutants, with transfer-predictive R2 values of > 0.62 for NO2 and < 0.56 for the other three pollutants. The LUR-model-based 500 m × 500 m spatial-distribution maps of these air pollutants illustrated pollution hot spots and the heterogeneity of population exposure, which provide valuable information for policymakers in designing effective air pollution control strategies. The LUR-model-based air pollution exposure estimates captured the spatial variability in exposure for participants in a cohort study. This study highlights that LUR models can be reasonably established upon a routine monitoring network, but there exist uncertainties when transferring LUR models between nearby cities. To the best of our knowledge, this study is the first to evaluate the intercity transferability of LUR models in Asia.

Dates
1 Introduction

Air pollution has been reported to be positively associated with a variety of health effect endpoints, such as lung function and respiratory and cardiovascular diseases (Çapraz et al., 2017; Sun et al., 2010; Yin et al., 2020; Zhou et al., 2020). Exposure assessment of air pollution is a critical component of epidemiological studies (Cai et al., 2020; Hoek et al., 2008; Li et al., 2017). Cohort studies focusing on the long-term effect on specific diseases of exposure to air pollution require accurate exposure estimates for a large group of participants (e.g., thousands or more) over a defined time period (Brokamp et al., 2019; Morley and Gulliver, 2018; Zhou et al., 2020). Different air quality prediction methods, such as air dispersion models, atmospheric chemical transport models, satellite remote sensing and various statistical methods, have been developed and applied to estimate air pollution (Yim et al., 2019a, b; Tong et al., 2018a, b; Luo et al., 2018; Shi et al., 2020a) and population exposure (Gu and Yim 2016; Gu et al., 2018; Hao et al., 2016; Li et al., 2020; Hou et al., 2019; Michanowicz et al., 2016; Wang et al., 2019, 2020; Yim et al., 2019c). Among these exposure assessment methods, land-use regression (LUR) is a widely used modeling approach to characterize long-term-average air pollutant concentrations at a fine spatial scale, which provides high-spatial-resolution estimates of exposure for use in epidemiological studies (Bertazzon et al., 2015; Eeftens et al., 2016; Jones et al., 2020).

The LUR method is based on the principle that ambient air pollutant concentrations at fixed-site measurement stations are linearly associated with different environmental features (e.g., land use, population density, road network and meteorological conditions) surrounding the stations (Anand and Monks, 2017; Lu et al., 2020; Naughton et al., 2018; Wu et al., 2017). In a city or even at a smaller-spatial-scale area, the LUR method is comparable to or sometimes even better than the approaches of satellite-remote-sensing-based air quality retrievals and air dispersion models in characterizing spatiotemporal variation in air pollution (Marshall et al., 2008; Shi et al., 2020b). Following feasible procedures of data processing and analysis, established air pollution LUR models can be applied to predict concentrations of air pollutants at locations without measurements at multiple spatial scales or at residential locations of participants in epidemiological studies (Li et al., 2021; Liu et al., 2016; Shi et al., 2020b).

In recent years, a large number of air pollution LUR studies have been conducted in different areas around the world (Jones et al., 2020; Lee et al., 2017; Liu et al., 2016, 2019; Lu et al., 2020; Miri et al., 2019; Ross et al., 2007; Wu et al., 2017). However, the development and application of LUR models in the Taiwan region were limited (Hsu et al., 2019). In addition, most previous Taiwan LUR studies used data from purpose-designed monitoring networks or combined purpose-designed and routine monitoring networks (Ho et al., 2015; Lee et al., 2014, 2015). For example, Lee et al. (2015) established LUR models for ambient particles of aerodynamic diameter less than or equal to 2.5 µm (PM2.5) using a purpose-designed monitoring network of 20 sites in the Taipei metropolis. The purpose-designed monitoring campaign has the advantage of capturing short-term air pollution exposure profiles (Jones et al., 2020), but it typically requires extra human labor and resources (e.g., experimental materials) (Hoek et al., 2008). Moreover, it is almost impossible to conduct long-term measurement (e.g., over years) using purpose-designed monitoring networks (Ho et al., 2015; Lee et al., 2017). As a result, a general limitation of LUR models upon purpose-designed monitoring networks is that the established models may only reflect the situation during the measurement period (Hoek et al., 2008; Shi et al., 2020b). Therefore, the development of long-term average LUR models for specific air pollutants using only routine monitoring networks should be explored, which is especially critical for epidemiological studies.

The application of established LUR models to areas outside the study area can reduce extra efforts to develop new models (Poplawski et al., 2009). To date, a few studies have evaluated the transferability of air pollution LUR models within a city and between cities or countries (Allen et al., 2011; Patton et al., 2015; Vienneau et al., 2010; Yang et al., 2020). Direct transferability refers to predictor variables and coefficients of LUR models both being transferred (Allen et al., 2011), whereas transferability with calibration means that model coefficients are calibrated using air pollutant measurements from the target areas (Yang et al., 2020). Direct transferability is more meaningful because it can be applied in areas without air quality measurements (Allen et al., 2011; Yang et al., 2020). Previous studies on the transferability of LUR models concluded that the predictive performances of LUR models from one area to another were not consistent, ranging from poor (Marcon et al., 2015) to relatively acceptable predictive accuracy (Poplawski et al., 2009; Wang et al., 2014). Therefore, more studies should be conducted to assess the transferability of air pollution LUR models.

In this study, annual-average LUR models and spatial-distribution maps were developed for ambient particles of aerodynamic diameter less than or equal to 10 µm (PM10), PM2.5, nitrogen dioxide (NO2) and ozone (O3) in northern Taiwan in 2019. In addition, the transferability of LUR models between cities in the study area was evaluated. The remainder of this paper is organized as follows: the “Materials and methods” section describes the study area, data collection and processing, LUR model establishment and validation, and prediction of the air pollution exposure surface. The “Results and discussion” section presents an overview of measurement data; established LUR models and their comparison with previous LUR models in Taiwan; the transferability of LUR models; the spatial-distribution maps of ambient PM10, PM2.5, NO2 and O3 concentrations; and PM2.5 exposure estimates for a cohort study. The “Conclusions” section summarizes the main results and demonstrates the implications of the present study.

2 Materials and methods

2.1 Study area

The Taipei–Keelung metropolitan area (TKMA), located in northern Taiwan, includes Taipei City, New Taipei City and Keelung City. The TKMA is the political, cultural and socioeconomic center of Taiwan. It covers an area of approximately 2457 km2 and has 48 administrative districts (Chiu et al., 2019; Wang et al., 2018). The TKMA had a population of about 7.03 million in 2019 (TWMOI, 2020), accounting for approximately 30 % of the total population of Taiwan (Fig. 1a). The population densities of Taipei City, New Taipei City and Keelung City were 10 175, 2021 and 2826 people km−2, respectively, in 2019 (TWMOI, 2020). The numbers of registered motor vehicles were 1.76 million, 3.21 million and 0.28 million in Taipei City, New Taipei City and Keelung City, respectively, by the end of 2018 (TWMOTC, 2020).

https://acp.copernicus.org/articles/21/5063/2021/acp-21-5063-2021-f01

Figure 1The characteristics of the study area. (a) Population density and the location of air-quality-monitoring stations. A total of 30 air-quality-monitoring stations were included in this study. (b) Digital elevation. (c) Land-use types. (d) The road network.

The TKMA is situated in the subtropical region and on the downwind side of mainland China. The built-up area of the TKMA is located in the central part of the Tamsui River basin, surrounded by mountains, agricultural land and forests (Fig. 1b and c). The characteristics of the basin terrain can constrain the diffusion of polluted air masses and thus favor the accumulation of air pollution in urban areas (Yu and Wang, 2010). Local emission sources of air pollutants in the TKMA include vehicular exhaust, industrial emissions and various sources related to residential activities (e.g., cooking) (Chen et al., 2020; Ho et al., 2018; Wu et al., 2017). In winter time, the long-distance transport of dust and polluted air masses under the northeast monsoon from the Asian continent results in a significant increase in concentrations of air pollutants (Chi et al., 2017; Chou et al., 2010).

2.2 Data collection and processing

The Taiwan Environmental Protection Administration (TWEPA) operates 20 central air-quality-monitoring stations in the TKMA, of which 12 stations are in New Taipei City, 7 are in Taipei City, and 1 station is in Keelung City (https://airtw.epa.gov.tw/ENG/default.aspx, last access: 9 July 2020). In addition, the Taipei Environmental Protection Agency (TPEPA) operates 10 local air-quality-monitoring stations (https://www.tldep.gov.taipei/EIACEP_EN/Air_NormalStation.aspx, last access: 9 July 2020). In total, these stations include 21 general stations, 6 traffic stations, 2 background stations and 1 country park station (Fig. 1a). Detailed descriptions of sampling stations, measurement instruments, and quality assurance and control procedures are available in TWEPA (2020). Hourly measurements of ambient PM10, PM2.5, NO2 and O3 concentrations and the meteorological variables of temperature, wind speed and relative humidity at the central stations from 1 January to 31 December 2019 were collected from the Environment Resource database of TWEPA (https://erdb.epa.gov.tw/DataRepository/EnvMonitor/AirQualityMonitorDayData.aspx, last access: 9 July 2020). In addition, hourly concentrations of ambient PM10, PM2.5, NO2 and O3 at the local stations from 1 January to 31 December 2019 were downloaded from the TPEPA website (https://www.tldep.gov.taipei/Public/DownLoad/AirAutoHour.aspx, last access: 9 July 2020). We calculated daily average values of air pollutant concentrations and meteorological variables from hourly data and calculated the annual-average values from daily averaged data for the development of LUR models. Daily and annual-average estimates for the air pollutants require at least 75 % data completeness (Cai et al., 2020); otherwise there was no value estimate for that day or year.

As presented in Table S1 in the Supplement and Fig. 1, the potential predictor variables of the road network, land-use data, normalized difference vegetation index (NDVI), population density and digital elevation data, which were frequently used in previous LUR studies, were collected. Land-use information was taken from the Land Use Investigation of Taiwan conducted by the National Land Surveying and Mapping Center (https://www2.nlsc.gov.tw/LUI/Home/Content_Home.aspx, last access: 9 July 2020). The Taiwan land-use status is classified into nine main categories, 41 subcategories and 103 detailed items. As shown in Fig. 1c, the nine main land-use categories are agriculture, forest, transportation, water bodies, built-up areas, public utilities, recreation, mining or salt production, and others (Chen et al., 2020). The road network from the Taiwan Ministry of Transportation and Communications includes three types of road: local roads, major roads and expressways (Fig. 1d). The NDVI and elevation data were extracted from the database of the Resources and Environmental Sciences Data Center, Chinese Academy of Sciences (http://www.resdc.cn, last access: 9 July 2020).

The values of potential predictor variables in buffer sizes of 50, 100, 300, 500, 700, 1000, 2000, 3000, 4000 and 5000 m surrounding the sampling stations were summarized for use in LUR model development. To ensure the consistency of results between model training and cross-validation, we included only the potential predictor variables with at least seven stations (i.e., around 25 % of all stations) exhibiting different values and where the minimum or maximum values lay within three times the 10th to the 90th percentile range below or above the 10th and the 90th percentile (Wolf et al., 2017).

2.3 Model development and validation

The LUR models of ambient PM10, PM2.5, NO2 and O3 for the entire study area (the area-specific LUR models) were established using all 30 air-quality-monitoring stations. In addition, city-specific LUR models for New Taipei City and Keelung City were developed using the 13 air-quality-monitoring stations located in these two cities, and the established models were directly transferred to Taipei City. Similarly, city-specific LUR models for Taipei City were developed using the 17 quality-monitoring stations located in this city, and the established models were directly transferred to New Taipei City and Keelung City. In this study, we did not consider the calibration of model coefficients because we planned to evaluate the direct transferability of city-specific LUR models to another nearby city area when there were no routine air quality measurements.

There is no standard modeling method for developing LUR models (Hoek et al., 2008). In this study, the supervised forward linear regression method (Cai et al., 2020; Eeftens et al., 2016; Xu et al., 2019) was used to develop the LUR models. This modeling method can ensure that only predictor variables following the plausible direction of effect are included, and meanwhile the predictive accuracy of the established model is maximized. In brief, all potential predictor variables were included as candidate independent variables, and a prior direction was assigned for each category of variable based on the atmospheric mechanism. The model construction started by including the predictor variable with the highest adjusted explained variance (R2). The remaining predictor variables were entered into the model if they met all of the following criteria: (1) the gain of the adjusted R2 was no less than 1 %; (2) the direction of effect of the predictor variable was pre-defined; (3) variables were added into the model when the probability of F was less than 0.05 and removed when the probability of F was greater than 0.10; (4) variables already included in the model retained the same direction of effect; and (5) following previous studies (Chen et al., 2020; Marcon et al., 2015; Wang et al., 2014), the predictor variables with variance inflation factor (VIF) values larger than 3 were dropped to make a tradeoff between model interpretation and the predictive accuracy (Eeftens et al., 2016). Multiple buffer sizes of a specific variable (e.g., the length of local roads) could be selected in the final model as long as they followed the selection criteria (Henderson et al., 2007).

Standard diagnostic tests were applied to ensure that the LUR models were reasonably established (Li, 2020; Wolf et al., 2017). The Cook's distance value was calculated to detect the outliers of data points (i.e., stations) (Jones et al., 2020). Air pollutant observations with a Cook's distance value greater than 1 would be excluded, and the LUR model for this air pollutant would be re-established (Weissert et al., 2018; Wolf et al., 2017). In addition, Moran's I values on the concentration residuals of the final LUR models were calculated using ArcGIS software to evaluate the spatial autocorrelation (Bertazzon et al., 2015; Lee et al., 2017; Liu et al., 2016). The R2 and root mean square error (RMSE) were estimated to evaluate the performance of the models (Li et al., 2021). Furthermore, leave-one-out cross-validation (LOOCV) was employed to evaluate the predictive capacity of the LUR models (Liu et al., 2019; Shi et al., 2020b; Yang et al., 2020).

Spatial analysis and calculations were performed using ArcGIS software, version 10.6 (ESRI Inc., Redlands, CA, USA). The statistical analysis was performed using R software, version 3.5.2 (R Core Team, 2018).

2.4 Air pollution surface prediction

The entire study area of the TKMA was divided into 9839 500 m × 500 m grid cells. The air pollutant concentrations at the centroids of the grid cells were estimated using the established area-specific LUR models. When the LUR models estimated negative concentration values, the concentration values of the grid cells were set to 0; when air pollutant concentration estimates exceeded the maximum observed concentrations by more than 20 %, the concentrations of grid cells were set to 120 % of the maximum observed concentrations (Henderson et al., 2007). The area-specific LUR-model-based negative and high concentration estimates accounted for only 0 %, 4 %, 2 % and 0 % of PM10, PM2.5, NO2 and O3 estimates, respectively. Then the spatial-distribution maps of ambient PM10, PM2.5, NO2 and O3 concentrations were created using the kriging interpolation method (Cai et al., 2020).

3 Results and discussion

3.1 Descriptive statistics of the air quality data

In general, the included air-quality-monitoring stations were situated at different types of land use across the TKMA (Table 1 and Fig. 1c), which suggests that the collected data set has relatively good representativeness. The annual-average PM10 concentration of 39.3 µg m−3 at background stations was the highest, followed in descending order by traffic stations with 33.6 µg m−3, general stations with 28.5 µg m−3 and the country park station with 15.7 µg m−3. The traffic stations and country park station had the highest and lowest annual-average PM2.5 concentrations, respectively. The annual-average PM2.5 concentrations at general stations of 13.7 µg m−3 and background stations of 13.2 µg m−3 were comparable. Except for the country park station, the annual-average PM10 and PM2.5 concentrations at other types of stations were higher than the air quality guidelines (AQGs) for PM10 and PM2.5 of 20.0 and 10.0 µg m−3, respectively, proposed by the World Health Organization (WHO) (WHO, 2006). The annual-average NO2 concentration of 24.6 ppb at the traffic stations was the highest, followed by general stations with 14.3 ppb. The annual-average NO2 concentrations at background stations (3.81 ppb) and the country park station (1.89 ppb) were significantly lower than those of general and traffic stations because they were farther away from traffic emissions. The annual-average NO2 concentration at traffic stations (24.6 ppb) was slightly higher than the WHO NO2 AQG of 40.0 µg m−3 (about 21.3 ppb) (WHO, 2006), while other types of stations had annual-average NO2 concentrations lower than this AQG. In contrast to NO2, the background stations (41.7 ppb) and the country park station (39.8 ppb) had higher annual-average O3 concentrations than those of traffic stations (21.6 ppb) or general stations (29.4 ppb) (Table 1).

Table 1Statistical description of measured air pollutants by different types of stations.

N means the number of stations for this type; SD means the standard deviation; min and max refer to the minimum and maximum values of the air pollutant concentrations, respectively. * There is only one country park station; therefore there are no estimates of SD, min and max values.

Download Print Version | Download XLSX

3.2 The area-specific LUR models

Figure S1 in the Supplement shows that Cook's distance values were below 1 for all the stations of the area-specific LUR models, suggesting that there were no station outliers in developing these LUR models. For PM10 and PM2.5 LUR models, Cook's distance values ranged from almost 0.00 to around 0.72. The Cook's distance values of the NO2 LUR model were between almost 0.00 and 0.28, whereas the Cook's distance values of the O3 LUR model were between almost 0.00 and 0.38 (Fig. S1). The final area-specific LUR models and their corresponding predictive accuracy are summarized in Table 2 and Fig. 2. The model R2 values ranged from 0.72 for PM2.5 to 0.91 for NO2, indicating a good fit for all air pollutants. PM10, NO2 and O3 LUR models performed well, with LOOCV R2 values being < 0.10 lower than the model R2 values. For PM2.5, the model was not as robust as those of other air pollutants, with the LOOCV R2 value being 0.19 lower than the model R2 value (Fig. 2). The reason for this is that the PM2.5 concentrations among the stations were not as discrete as those of other air pollutants (Table 1 and Fig. 2). The significance of the predictor variables (p value) and VIF values all met the requirements for LUR model development. Moran's I values were 0.0047, 0.072, 0.023 and 0.055 for the LUR models of ambient PM10, PM2.5, NO2 and O3. In addition, z-score values were 0.83, 0.79, 1.2 and 0.34 for ambient PM10, PM2.5, NO2 and O3 LUR models, respectively, indicating that the spatial patterns of concentration residuals of the LUR models do not appear to be significantly different from random (Fig. S2).

Table 2Description of the 2019 annual-average LUR models for ambient PM10, PM2.5, NO2 and O3 in the TKMA.

LU2_300, LU2_700: the area of forest in buffer sizes of 300 m and 700 m. LU5_50: the area of built-up land in a buffer size of 50 m. LU7_1000 and LU7_2000: the area of recreational land in buffer sizes of 1000 and 2000 m. LU3_50, LU3_100, LU3_300, LU3_700 and LU3_3000: the area of transportation land in buffer sizes of 50, 100, 300, 700 and 3000 m. LU4_500: the area of water body in a buffer size of 500 m. R2_100: the area of major roads in a buffer size of 100 m. Dis_Major: the distance to the nearest major roads. RL3_1000: the length of local roads in a buffer size of 1000 m. LU6_300: the area of public utilization land in a buffer size of 300 m. VIF: the variance inflation factor. LOOCV: leave-one-out cross-validation. RMSE: root mean square error. NA: not available.

Download Print Version | Download XLSX

https://acp.copernicus.org/articles/21/5063/2021/acp-21-5063-2021-f02

Figure 2A comparison of LUR-predicted concentrations and observed concentrations of the studied air pollutants and the LOOCV-predicted concentrations and observed concentrations of the studied air pollutants. (a) PM10, (b) PM2.5, (c) NO2 and (d) O3. N is the sample size, and the solid line is the 1:1 line.

Download

The final area-specific LUR models consisted of three (for O3), four (for NO2) and five predictor variables (for PM10 and PM2.5) (Table 2). Consistent with the previous LUR studies of De Hoogh et al. (2018), Eeftens et al. (2016), Jones et al. (2020), Weissert et al. (2018) and Wolf et al. (2017), the established LUR models contained at least one traffic-related predictor variable in buffer sizes ranging from 50 to 3000 m. Traffic emission is a major source of air pollution in urban areas of the TKMA (Lee et al., 2014; Wu et al., 2017). For instance, it was reported that gasoline and diesel vehicle emissions contributed approximately half of PM2.5 concentrations in Taipei City based on source apportionment analysis (Ho et al., 2018). Several previous LUR studies selected the population density variable as the final explanatory variable in their PM2.5 and NO2 LUR models (Ji et al., 2019; Meng et al., 2015; Rahman et al., 2017). However, it was not included in our final LUR models. A possible explanation is that the population density variable is moderately or highly correlated with the variables (e.g., the area of recreational land) included in our final LUR models.

As shown in Table 2, PM10 and PM2.5 LUR models included predictor variables of both small and large buffer sizes. The LUR model for PM10 included the area of forest land in a buffer size of 300 m, the area of built-up land in a buffer size of 50 m, the area of recreational land in a buffer size of 2000 m, the area of transportation land in a buffer size of 100 m and the area of water body land in a buffer size of 500 m. The PM2.5 LUR model included the area of transportation land within a 300 m buffer, the area of major roads within a 100 m buffer, the area of forest land within a 700 m buffer, the area of recreational land within a 2000 m buffer and the distance to the nearest major roads. For PM10 and PM2.5 LUR models, the direction of effect for transportation land and traffic roads was positive, while the direction of effect of other predictor variables was negative. Forest and urban green space land (i.e., recreational land) were included in both PM10 and PM2.5 LUR models (Table 2). Ji et al. (2019), Jones et al. (2020) and Miri et al. (2020) included forest land or urban green space as the predictor variables in their final city-scale PM LUR models, demonstrating the mitigation effect of these land-use types on PM concentrations. Chen et al. (2019) and Jeanjean et al. (2016) reported the effectiveness of urban green space in mitigating PM pollution. The water body type of land use reduced PM10 concentrations, as evidenced by the negative regression coefficient (Table 2). The water bodies can make PM10 absorb moisture and increase sedimentation. In addition, large areas of water provide good conditions for the dispersion of air pollutants (Zhu and Zhou, 2019).

For the NO2 LUR model, the four predictor variables included were the area of transportation land in buffer sizes of 3000 and 50 m, the area of recreational land in a 1000 m buffer, and the sum of the length of local roads in a 1000 m buffer. The direction of effect for the recreational land was negative, while other predictor variables showed a positive effect (Table 2). The O3 LUR model included predictor variables with relatively small buffer sizes of less than 700 m. The three predictor variables were the area of transportation land in buffer sizes of 700 and 50 m and the area of public utilization land within a 300 m buffer. The directions of effect for these three variables were all negative (Table 2). The traffic-related predictor variables were important variables in predicting NO2 and O3 concentrations but in different directions of effect. Consistent with previous studies by De Hoogh et al. (2016), Eeftens et al. (2016), Lee et al. (2014) and Liu et al. (2019), the established NO2 LUR model also revealed the mitigation effect of urban green space (i.e., recreational land) on NO2 concentration.

A comparison of this study with previous LUR studies in Taiwan is presented in Table S2. The predictive performance of the LUR model for ambient PM10 in this study was slightly worse than that of Lee et al. (2015), with an R2 value of 0.87. In addition, the R2 and LOOCV R2 values (0.72 and 0.53, respectively) of the PM2.5 LUR model in this study were lower than those of Ho et al. (2015) (an R2 value of 0.75 and an LOOCV R2 value of 0.62), Lee et al. (2015) (an R2 value of 0.95 and an LOOCV R2 value of 0.91) and Wu et al. (2017) (an R2 value of 0.90 and an LOOCV R2 value of 0.83) but higher than that of Wu et al. (2018), with an R2 value of 0.66. The NO2 LUR model performed better than that of Lee et al. (2014) and was comparable to that of Chen et al. (2020). Hsu et al. (2019) developed an O3 LUR model for the whole Taiwan region, with an R2 value of 0.74 (Hsu et al., 2019). Our study established a reasonable LUR model for ambient O3 in the TKMA with an R2 value of 0.80 and an LOOCV R2 value of 0.70, which is a relatively high predictive performance. Compared with PM10, PM2.5 and NO2, the establishment of O3 LUR models has been limited in these previous Taiwan LUR studies (Table S2) or in most of the LUR studies in other areas, but it is essential to establish O3 LUR models given that O3 is a toxic photochemical pollutant threatening human health and the ecosystem (Ning et al., 2020; Yim et al., 2019b).

https://acp.copernicus.org/articles/21/5063/2021/acp-21-5063-2021-f03

Figure 3The changes in R2 values for direct transfer of ambient PM10, PM2.5, NO2 and O3 LUR models between Taipei City and New Taipei City and Keelung City.

Download

3.3 Transferability of the city-specific LUR models

The city-specific LUR models for ambient PM10, PM2.5, NO2 and O3 in Taipei City, New Taipei City and Keelung City are shown in Tables S3 and S4, respectively. The model R2 values of the Taipei City PM10, PM2.5, NO2 and O3 LUR models were 0.91, 0.64, 0.89 and 0.76, respectively (Table S3), while the New Taipei City and Keelung City PM10, PM2.5, NO2 and O3 LUR models had R2 values of 0.63, 0.65, 0.95 and 0.93, respectively (Table S4). In general, for each specific air pollutant, the predictive performance of these city-specific LUR models can be slightly higher or lower than those of the area-specific LUR models. Figure 3 shows the transferability of LUR models between Taipei City and New Taipei City and Keelung City. The city-specific LUR models performed worse in another city area than in the city where these models were established. For instance, the transfer-predictive R2 values of the Taipei LUR models were 0.31, 0.04, 0.62 and 0.56 for predicting ambient PM10, PM2.5, NO2 and O3 in New Taipei City and Keelung City, respectively (Fig. 3). These values were substantially lower than the corresponding R2 values of the Taipei LUR models. The NO2 LUR models showed good transferability between the two city areas, with transfer-predictive R2 values higher than 0.62. However, the PM10, PM2.5 and O3 LUR models performed poorly when they were transferred between the two city areas, with transfer-predictive R2 values of < 0.31, < 0.37 and < 0.56, respectively (Fig. 3). Similar to the previous studies of Marcon et al. (2015) and Yang et al. (2020), these results suggested that there may be large uncertainties in transferring LUR models between cities and even between nearby cities with similar geographic and urban-design characteristics. The use of novel cost-effective methods (e.g., low-cost air quality sensors or a satellite remote-sensing approach) is therefore recommended to assess air pollution and associated population exposure in cities with limited fixed-site measurement stations.

3.4 Spatial maps

LUR-model-derived air-pollution-spatial-distribution maps provide valuable and useful air pollutant concentration surfaces in the TKMA. In general, there was a good agreement between LUR-model-based concentration estimates and observations for PM10, PM2.5, NO2 and O3 (Fig. 4). For PM10 and PM2.5, there were certain differences between LUR-model-based concentration estimates and observations at the country park station (Fig. 4). A possible reason for this difference may be that the kriging interpolation method removed low-concentration estimates at this small area when the concentration estimates at nearby areas were higher.

https://acp.copernicus.org/articles/21/5063/2021/acp-21-5063-2021-f04

Figure 4The spatial distribution of ambient air pollutant concentrations derived from established LUR models. (a) PM10, (b) PM2.5, (c) NO2 and (d) O3. The colored circles represent the observations from stations.

Table 3Pearson correlation coefficients (PCCs) among the estimated concentrations of ambient PM10, PM2.5, NO2 and O3.

** Correlation is significant at the 0.01 level (two-tailed).

Download Print Version | Download XLSX

High concentrations of ambient PM10, PM2.5 and NO2 were predicted in the urban areas of Taipei City, New Taipei City and Keelung City and along the road network. The estimated PM10 and PM2.5 concentrations in urban areas were around 35.0 to 40.9 µg m−3 and around 12.0 to 17.0 µg m−3, respectively, whereas the urban areas had NO2 concentrations of around 12.0 to 31.7 ppb (Fig. 4). This spatial-distribution pattern is understandable given that the traffic-related predictor variables were included in the final PM10, PM2.5 and NO2 LUR models. A similar spatial pattern of PM2.5 concentrations was reported by Wu et al. (2017), which documented that high PM2.5 concentrations were distributed mainly in the urban areas of the TKMA, and there were also scattered points of high PM2.5 concentrations in its outer ring. However, the estimated 2019 annual-average PM2.5 concentrations in this study were significantly lower than those for 2006–2012 estimated by Wu et al. (2017). There was a clear decreasing trend in PM2.5 concentrations in the whole of Taiwan over the past decade (Ho et al., 2020; Jung et al., 2018). For example, Jung et al. (2018) reported that the estimated PM2.5 concentrations declined by 1.7 and 1.6 µg m−3 in the morning and afternoon, respectively, per year over the whole of Taiwan during the period 2005–2015. O3 showed a generally opposite spatial-variability pattern compared with the other three air pollutants, with lower concentrations (< about 32.0 ppb) in urban areas than in rural areas (Fig. 4). A possible explanation for this finding is that high concentrations of NO and NO2 in urban areas react with O3, resulting in a decrease in O3 concentration (Hsu et al., 2019; Vardoulakis et al., 2011).

Correlations of estimated concentrations of PM10, PM2.5, NO2 and O3 in the TKMA are shown in Table 3. Consistent with previous studies by Hoek et al. (2008), Lu et al. (2020), Vardoulakis et al. (2011) and Wolf et al. (2017), the spatial-distribution maps revealed high spatial correlations among the four air pollutants. PM10 concentrations had strong positive correlations with PM2.5 and NO2, suggesting common sources of these three air pollutants. In contrast to this, PM10 concentrations were negatively correlated with O3 concentrations, with a Pearson correlation coefficient (PCC) value of 0.730. Similarly, PM2.5 concentrations had a strong positive correlation with NO2 concentrations but showed a significant negative correlation with O3 concentrations. The concentrations of NO2 and O3 were negatively correlated because of the NOx titration effect in urban areas, with a PCC value of 0.920. Similar findings were reported by De Hoogh et al. (2018) and Lu et al. (2020).

3.5 Air pollutant exposure estimates for a cohort study

Air pollutant concentrations measured at nearby fixed-site stations are often used to represent exposures in epidemiological studies (Lin et al., 2016; Shi et al., 2020c), but the spatial resolution of these estimates is relatively coarse due to the limited number of sampling stations (Bertazzon et al., 2015). In recent years, LUR modeling has become a more widely applied method to estimate air pollution exposures at a fine spatial scale (Lee et al., 2014; Wolf et al., 2017). Figure S3 shows that there are differences between LUR-model-based air pollution exposure estimates and nearby-station measurements at residential locations of participants in a cohort study conducted in the TKMA. The average values of the LUR-estimated PM10, PM2.5, NO2 and O3 exposure concentrations were 36.0 µg m−3, 14.2 µg m−3, 18.0 ppb and 29.2 ppb, respectively, whereas the corresponding nearby-station measurements were 27.7 µg m−3, 13.8 µg m−3, 16.3 ppb and 28.6 ppb, respectively (Table S5). Compared with LUR-model-based estimates, the nearby-station measurements underestimated PM10, PM2.5, NO2 and O3 exposures of cohort participants by 8.23 µg m−3, 0.41 µg m−3, 1.73 ppb and 0.60 ppb, respectively (Table S5). In addition, the concentration ranges of LUR-estimated annual-average PM10 (13.0–45.2 µg m−3), PM2.5 (6.96–19.9 µg m−3), NO2 (0.70–32.2 ppb) and O3 (17.5–44.0 ppb) exposure concentrations were larger than those of nearby-station measurements for PM10 (22.3–40.3 µg m−3), PM2.5 (10.6–21.3 µg m−3), NO2 (2.90–32.2 ppb) and O3 (15.2–42.2 ppb) (Table S5 and Fig. 5). This indicates that the LUR-model-based exposure estimates can capture the large spatial variability in air pollutant exposure among the cohort participants. Similar findings have been reported in studies of Lee et al. (2014) and Marshall et al. (2008). Furthermore, the LUR-model-based PM10, PM2.5, NO2 and O3 exposure estimates and nearby-station measurements were weakly correlated, with linear regression R2 values ranging from 0.05 for PM10 to 0.19 for NO2 (Fig. 6). A possible explanation is that LUR-model-based exposure estimates generally accounted for neighborhood-scale variations in air pollutant concentrations, while the nearby-station measurements usually only revealed the urban-scale variability in air pollution (e.g., urban area versus suburban area versus rural area) (Marshall et al., 2008). The LUR-model-based exposure estimates and nearby-station measurements should be further validated if the air quality measurement data at residential locations of cohort participants (if not all, at least some of the participants) are available.

https://acp.copernicus.org/articles/21/5063/2021/acp-21-5063-2021-f05

Figure 5Box plots of nearby-station air pollutant measurements and LUR-model-based estimates of air pollutant concentration. (a) PM10, (b) PM2.5, (c) NO2 and (d) O3. The triangle symbol in each box is the mean value, the solid line is the median value, the box extends from the 25th to the 75th percentile, the whiskers (error bars) below and above the box are the 10th and 90th percentiles, and the lower and upper cycle symbols are outliers.

Download

https://acp.copernicus.org/articles/21/5063/2021/acp-21-5063-2021-f06

Figure 6The linear regression of nearby-station air pollutant measurements and LUR-model-based air pollutant concentration estimates. (a) PM10, (b) PM2.5, (c) NO2 and (d) O3.

Download

3.6 Limitations

This study is subject to several limitations. First, apart from the variables used in this study, more predictor variables (e.g., localized emission data and urban building morphology data) should be included and tested to develop LUR models. For example, Wu et al. (2017) and Chen et al. (2020) assessed the roles of two culturally specific emission sources, Chinese restaurants and temples, on the development of ambient PM2.5 and NO2 LUR models in Taiwan. More studies should be conducted to test the influence of different potential predictor variables on the development of LUR models (Hoek et al., 2008). Second, like most linear regression techniques, the supervised forward linear regression method is not proficient in modeling extreme values (Jones et al., 2020). In addition, there may be complex and non-linear relationships between the explanatory variables and air pollutant concentrations (Wang et al., 2020). Other types of linear regression methods (Hoek et al., 2008; Shi et al., 2020b) and the novel machine learning algorithms (Wang et al., 2020) can be tested in estimating surface-level air pollutant concentrations in the further study. Third, the kriging interpolation method tends to remove air pollutant peak concentrations, resulting in an underestimation of air pollution exposure at pollution hot spots. Other spatial-mapping methods should be considered in further studies. It is recommended that air pollutant concentrations at residential locations of participants should be estimated directly for cohort studies. Fourth, there may be uncertainty in spatial estimations of air pollutant concentrations with a limited number of sampling stations. Further studies are warranted to evaluate the influence of the number of sampling stations and their spatial distributions on the development of LUR models and the air pollution spatial maps.

4 Conclusions

Following standard development procedures, the annual-average LUR models of ambient PM10, PM2.5, NO2 and O3 were established in the TKMA of northern Taiwan using only data from the routine monitoring network. These LUR models were reasonable, based on the evaluation metrics of Cook's distance, VIF, Moran's I and p values. The R2 values of the LUR models for ambient PM10, PM2.5, NO2 and O3 were 0.80, 0.72, 0.91 and 0.80, respectively. The traffic-related predictor variables were the major explanatory factors in the LUR models for all the studied air pollutants.

The predictive performance varied greatly among air pollutants in examining the transferability of city-specific LUR models between New Taipei City and Keelung City and Taipei City, with relatively high transfer-predictive R2 values for NO2. Therefore, this study highlights that the established LUR models in a city area can result in a large estimation bias when applied to another nearby city area with similar geographic and urbanization conditions. The transferability may even be uncertain in a city with complex terrain (Yim et al., 2007, 2014). It is necessary to conduct more studies to evaluate and improve the intercity transferability of LUR models.

The spatial-distribution maps of the four air pollutants showed that the developed LUR models are reasonable in modeling the spatial variabilities in air pollution. Ambient PM10, PM2.5 and NO2 shared similar spatial variations, with relatively high concentrations in urban areas and along the road network. Ambient O3 presented a generally opposite spatial variability compared with PM10, PM2.5 and NO2. These estimated air pollution concentration surfaces provide information for the management of air pollution and exposure estimates for epidemiological studies. Compared with nearby-station measurements, the LUR-model-based concentration estimates captured a wider range of exposure to PM10, PM2.5, NO2 and O3 for participants in a cohort study in the TKMA. Further studies should pay more attention to utilizing other data sources (e.g., satellite remote-sensing data) with comprehensive spatiotemporal coverage to validate the LUR-model-based estimations of air pollutant concentrations.

Data availability

The model data presented in this article are available from the authors upon request (yimsteve@gmail.com).

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/acp-21-5063-2021-supplement.

Author contributions

SHLY planned, supervised and sought funding for this study. ZL performed the data analysis and prepared the paper with contributions from all co-authors.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Air Quality Research at Street-Level (ACP/GMD inter-journal SI)”. It is not associated with a conference.

Acknowledgements

We would like to thank the Taiwan Environmental Protection Administration for providing air quality and meteorological data.

Financial support

This work is funded by the Vice-Chancellor's Discretionary Fund of The Chinese University of Hong Kong (grant no. 4930744), the Dr. Stanley Ho Medicine Development Foundation (grant no. 8305509), and the project from the ENvironmental SUstainability and REsilience (ENSURE) partnership between the CUHK and UoE.

Review statement

This paper was edited by Sunling Gong and reviewed by two anonymous referees.

References

Allen, R. W., Amram, O., Wheeler, A. J., and Brauer, M.: The transferability of NO and NO2 land use regression models between cities and pollutants, Atmos. Environ., 45, 369–378, 2011. 

Anand, J. S. and Monks, P. S.: Estimating daily surface NO2 concentrations from satellite data – a case study over Hong Kong using land use regression models, Atmos. Chem. Phys., 17, 8211–8230, https://doi.org/10.5194/acp-17-8211-2017, 2017. 

Bertazzon, S., Johnson, M., Eccles, K., and Kaplan, G. G.: Accounting for spatial effects in land use regression for urban air pollution modeling, Spatial and Spatiotemporal Epidemiology, 14, 9–21, 2015. 

Brokamp, C., Brandt, E. B., and Ryan, P. H.: Assessing exposure to outdoor air pollution for epidemiological studies: Model-based and personal sampling strategies, J. Allergy Clin. Immun., 143, 2002–2006, 2019. 

Cai, J., Ge, Y., Li, H., Yang, C., Liu, C., Meng, X., Wang, W., Niu, C., Kan, L., Schikowski, T., and Yan, B.: Application of land use regression to assess exposure and identify potential sources in PM2.5, BC, NO2 concentrations, Atmos. Environ., 223, 117267, https://doi.org/10.1016/j.atmosenv.2020.117267, 2020. 

Çapraz, Ö., Deniz, A., and Doğan, N.: Effects of air pollution on respiratory hospital admissions in İstanbul, Turkey, 2013 to 2015, Chemosphere, 181, 544–550, 2017. 

Chen, M., Dai, F., Yang, B., and Zhu, S.: Effects of neighborhood green space on PM2.5 mitigation: Evidence from five megacities in China, Build. Environ., 156, 33–45, 2019. 

Chen, T. H., Hsu, Y. C., Zeng, Y. T., Lung, S. C. C., Su, H. J., Chao, H. J., and Wu, C. D.: A hybrid kriging/land-use regression model with Asian culture-specific sources to assess NO2 spatial-temporal variations, Environ. Pollut., 259, 113875,https://doi.org/10.1016/j.envpol.2019.113875, 2020. 

Chi, K. H., Li, Y. N., and Hung, N. T.: Spatial and temporal variation of PM2.5 and atmospheric PCDD/FS in Northern Taiwan during winter monsoon and local pollution episodes, Aerosol Air Qual. Res., 17, 3151–3165, 2017. 

Chiu, H. W., Lee, Y. C., Huang, S. L., and Hsieh, Y. C.: How does periurbanization teleconnect remote areas? An emergy approach, Ecol. Model., 403, 57–69, 2019. 

Chou, C. C.-K., Lee, C. T., Cheng, M. T., Yuan, C. S., Chen, S. J., Wu, Y. L., Hsu, W. C., Lung, S. C., Hsu, S. C., Lin, C. Y., and Liu, S. C.: Seasonal variation and spatial distribution of carbonaceous aerosols in Taiwan, Atmos. Chem. Phys., 10, 9563–9578, https://doi.org/10.5194/acp-10-9563-2010, 2010. 

De Hoogh, K., Gulliver, J., van Donkelaar, A., Martin, R. V., Marshall, J. D., Bechle, M. J., Cesaroni, G., Pradas, M. C., Dedele, A., Eeftens, M., and Forsberg, B.: Development of West-European PM2.5 and NO2 land use regression models incorporating satellite-derived and chemical transport modelling data, Environ. Res., 151, 1–10, 2016. 

De Hoogh, K., Chen, J., Gulliver, J., Hoffmann, B., Hertel, O., Ketzel, M., Bauwelinck, M., van Donkelaar, A., Hvidtfeldt, U. A., Katsouyanni, K., and Klompmaker, J.: Spatial PM2.5, NO2, O3 and BC models for Western Europe–Evaluation of spatiotemporal stability, Environ. Int., 120, 81–92, 2018. 

Eeftens, M., Meier, R., Schindler, C., Aguilera, I., Phuleria, H., Ineichen, A., Davey, M., Ducret-Stich, R., Keidel, D., Probst-Hensch, N., and Künzli, N.: Development of land use regression models for nitrogen dioxide, ultrafine particles, lung deposited surface area, and four other markers of particulate matter pollution in the Swiss SAPALDIA regions, Environ. Health, 15, 53, https://doi.org/10.1186/s12940-016-0137-9, 2016. 

Gu, Y. and Yim, S. H. L.: The air quality and health impacts of domestic trans-boundary pollution in various regions of China, Environ. Int., 97, 117–124, 2016. 

Gu, Y., Wong, T. W., Law, C. K., Dong, G. H., Ho, K. F., Yang, Y., and Yim, S. H. L.: Impacts of sectoral emissions in China and the implications: air quality, public health, crop production, and economic costs, Environ. Res. Lett., 13, 084008, https://doi.org/10.1088/1748-9326/aad138, 2018. 

Hao, H., Chang, H. H., Holmes, H. A., Mulholland, J. A., Klein, M., Darrow, L. A., and Strickland, M. J.: Air pollution and preterm birth in the US State of Georgia (2002–2006): associations with concentrations of 11 ambient air pollutants estimated by combining Community Multiscale Air Quality Model (CMAQ) simulations with stationary monitor measurements, Environ. Health Persp., 124, 875–880, 2016. 

Henderson, S. B., Beckerman, B., Jerrett, M., and Brauer, M.: Application of land use regression to estimate long-term concentrations of traffic-related nitrogen oxides and fine particulate matter, Environ. Sci. Technol., 41, 2422–2428, 2007. 

Ho, C. C., Chan, C. C., Cho, C. W., Lin, H. I., Lee, J. H., and Wu, C. F.: Land use regression modeling with vertical distribution measurements for fine particulate matter and elements in an urban area, Atmos. Environ., 104, 256–263, 2015. 

Ho, C. C., Chen, L. J., and Hwang, J. S.: Estimating ground-level PM2.5 levels in Taiwan using data from air quality monitoring stations and high coverage of microsensors, Environ. Pollut., 264, 114810, https://doi.org/10.1016/j.envpol.2020.114810, 2020. 

Ho, W. Y., Tseng, K. H., Liou, M. L., Chan, C. C., and Wang, C. H.: Application of positive matrix factorization in the identification of the sources of PM2.5 in Taipei City, Int. J. Env. Res. Pub. He., 15, 1305, https://doi.org/10.3390/ijerph15071305, 2018. 

Hoek, G., Beelen, R., De Hoogh, K., Vienneau, D., Gulliver, J., Fischer, P., and Briggs, D.: A review of land-use regression models to assess spatial variation of outdoor air pollution, Atmos. Environ., 42, 7561–7578, 2008. 

Hou, X., Chan, C. K., Dong, G. H., and Yim, S. H. L.: Impacts of transboundary air pollution and local emissions on PM2.5 pollution in the Pearl River Delta region of China and the public health, and the policy implications, Environ. Res. Lett., 14, 034005, https://doi.org/10.1088/1748-9326/aaf493, 2019. 

Hsu, C. Y., Wu, J. Y., Chen, Y. C., Chen, N. T., Chen, M. J., Pan, W. C., Lung, S. C. C., Guo, Y. L., and Wu, C. D.: Asian culturally specific predictors in a large-scale land use regression model to predict spatial-temporal variability of ozone concentration, Int. J. Env. Res. Pub. He., 16, 1300, https://doi.org/10.3390/ijerph16071300, 2019. 

Jeanjean, A. P. R., Monks, P. S., and Leigh, R. J.: Modelling the effectiveness of urban trees and grass on PM2.5 reduction via dispersion and deposition at a city scale, Atmos. Environ., 147, 1–10, 2016. 

Ji, W., Wang, Y., and Zhuang, D.: Spatial distribution differences in PM2.5 concentration between heating and non-heating seasons in Beijing, China, Environ. Pollut., 248, 574–583, 2019. 

Jones, R. R., Hoek, G., Fisher, J. A., Hasheminassab, S., Wang, D., Ward, M. H., Sioutas, C., Vermeulen, R., and Silverman, D. T.: Land use regression models for ultrafine particles, fine particles, and black carbon in southern California, Sci. Total Environ., 699, 134234, https://doi.org/10.1016/j.scitotenv.2019.134234, 2020. 

Jung, C. R., Hwang, B. F., and Chen, W. T.: Incorporating long-term satellite-based aerosol optical depth, localized land use data, and meteorological variables to estimate ground-level PM2.5 concentrations in Taiwan from 2005 to 2015, Environ. Pollut., 237, 1000–1010, 2018. 

Lee, J. H., Wu, C. F., Hoek, G., de Hoogh, K., Beelen, R., Brunekreef, B., and Chan, C. C.: Land use regression models for estimating individual NOx and NO2 exposures in a metropolis with a high density of traffic roads and population, Sci. Total Environ., 472, 1163–1171, 2014. 

Lee, J. H., Wu, C. F., Hoek, G., de Hoogh, K., Beelen, R., Brunekreef, B., and Chan, C. C.: LUR models for particulate matters in the Taipei metropolis with high densities of roads and strong activities of industry, commerce and construction, Sci. Total Environ., 514, 178–184, 2015. 

Lee, M., Brauer, M., Wong, P., Tang, R., Tsui, T. H., Choi, C., Cheng, W., Lai, P. C., Tian, L., Thach, T. Q., and Allen, R.: Land use regression modelling of air pollution in high density high rise cities: A case study in Hong Kong, Sci. Total Environ., 592, 306–315, 2017. 

Li, Q. X.: Statistical modelling experiment of land precipitation variations since the start of the 20th century with external forcing factors, China Sci. Bull., 65, 2266–2278, 2020 (in Chinese). 

Li, Z., Che, W., Frey, H. C., Lau, A. K., and Lin, C.: Characterization of PM2.5 exposure concentration in transport microenvironments using portable monitors, Environ. Pollut., 228, 433–442, 2017. 

Li, Z., Yim, S. H. L., and Ho, K. F.: High temporal resolution prediction of street-level PM2.5 and NOx concentrations using machine learning approach, J. Clean. Prod., 268, 121975, https://doi.org/10.1016/j.jclepro.2020.121975, 2020. 

Li, Z., Tong, X., Ho, J. M. W., Kwok, T. C., Dong, G., Ho, K. F., and Yim, S. H. L.: A practical framework for predicting residential indoor PM2.5 concentration using land-use regression and machine learning methods, Chemosphere, 265, 129140, https://doi.org/10.1016/j.chemosphere.2020.129140, 2021. 

Lin, H., Liu, T., Xiao, J., Zeng, W., Li, X., Guo, L., Zhang, Y., Xu, Y., Tao, J., Xian, H., and Syberg, K. M.: Mortality burden of ambient fine particulate air pollution in six Chinese cities: results from the Pearl River Delta study, Environ. Int., 96, 91–97, 2016. 

Liu, C., Henderson, B. H., Wang, D., Yang, X., and Peng, Z. R.: A land use regression application into assessing spatial variation of intra-urban fine particulate matter (PM2.5) and nitrogen dioxide (NO2) concentrations in City of Shanghai, China, Sci. Total Environ., 565, 607–615, 2016. 

Liu, Z., Guan, Q., Luo, H., Wang, N., Pan, N., Yang, L., Xiao, S., and Lin, J.: Development of land use regression model and health risk assessment for NO2 in different functional areas: A case study of Xi'an, China, Atmos. Environ., 213, 515–525, 2019. 

Lu, M., Soenario, I., Helbich, M., Schmitz, O., Hoek, G., van der Molen, M., and Karssenberg, D.: Land use regression models revealing spatiotemporal co-variation in NO2, NO, and O3 in the Netherlands, Atmos. Environ., 223, 117238, https://doi.org/10.1016/j.atmosenv.2019.117238, 2020. 

Luo, M., Hou, X., Gu, Y., Lau, N. C., and Yim, S. H. L.: Trans-boundary air pollution in a city under various atmospheric conditions, Sci. Total Environ., 618, 132–141, 2018. 

Marcon, A., de Hoogh, K., Gulliver, J., Beelen, R., and Hansell, A. L.: Development and transferability of a nitrogen dioxide land use regression model within the Veneto region of Italy, Atmos. Environ., 122, 696–704, 2015. 

Marshall, J. D., Nethery, E., and Brauer, M.: Within-urban variability in ambient air pollution: comparison of estimation methods, Atmos. Environ., 42, 1359–1369, 2008. 

Meng, X., Chen, L., Cai, J., Zou, B., Wu, C. F., Fu, Q., Zhang, Y., Liu, Y., and Kan, H.: A land use regression model for estimating the NO2 concentration in Shanghai, China, Environ. Res., 137, 308–315, 2015. 

Michanowicz, D. R., Shmool, J. L., Tunno, B. J., Tripathy, S., Gillooly, S., Kinnee, E., and Clougherty, J. E.: A hybrid land use regression/AERMOD model for predicting intra-urban variation in PM2.5, Atmos. Environ., 131, 307–315, 2016. 

Miri, M., Ghassoun, Y., Dovlatabadi, A., Ebrahimnejad, A., and Löwner, M. O.: Estimate annual and seasonal PM1, PM2.5 and PM10 concentrations using land use regression model, Ecotox. Environ. Safe., 174, 137–145, 2019. 

Morley, D. W. and Gulliver, J.: A land use regression variable generation, modelling and prediction tool for air pollution exposure assessment, Environ. Modell. Softw., 105, 17–23, 2018. 

Naughton, O., Donnelly, A., Nolan, P., Pilla, F., Misstear, B. D., and Broderick, B.: A land use regression model for explaining spatial variation in air pollution levels using a wind sector based approach, Sci. Total Environ., 630, 1324–1334, 2018. 

Ning, G., Yim, S. H. L., Yang, Y., Gu, Y., and Dong, G.: Modulations of synoptic and climatic changes on ozone pollution and its health risks in mountain-basin areas, Atmos. Environ., 240, 117808, https://doi.org/10.1016/j.atmosenv.2020.117808, 2020. 

Patton, A. P., Zamore, W., Naumova, E. N., Levy, J. I., Brugge, D., and Durant, J. L.: Transferability and generalizability of regression models of ultrafine particles in urban neighborhoods in the Boston area, Environ. Sci. Technol, 49, 6051–6060, 2015. 

Poplawski, K., Gould, T., Setton, E., Allen, R., Su, J., Larson, T., Henderson, S., Brauer, M., Hystad, P., Lightowlers, C., and Keller, P.: Intercity transferability of land use regression models for estimating ambient concentrations of nitrogen dioxide, J. Expo. Sci. Env. Epid., 19, 107–117, 2009. 

R Core Team: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, available at: https://www.Rproject.org/ (last access: 9 July 2020), 2018. 

Rahman, M. M., Yeganeh, B., Clifford, S., Knibbs, L. D., and Morawska, L.: Development of a land use regression model for daily NO2 and NOx concentrations in the Brisbane metropolitan area, Australia, Environ. Modell. Softw., 95, 168–179, 2017. 

Ross, Z., Jerrett, M., Ito, K., Tempalski, B., and Thurston, G. D.: A land use regression for predicting fine particulate matter concentrations in the New York City region, Atmos. Environ., 41, 2255–2269, 2007. 

Shi, C., Nduka, I. C., Yang, Y., Huang, Y., Yao, R., Zhang, H., He, B., Xie, C., Wang, Z., and Yim, S. H. L.: Characteristics and meteorological mechanisms of transboundary air pollution in a persistent heavy PM2.5 pollution episode in Central-East China, Atmos. Environ., 223, 117239, https://doi.org/10.1016/j.atmosenv.2019.117239, 2020a. 

Shi, T., Dirienzo, N., Requia, W. J., Hatzopoulou, M., and Adams, M. D.: Neighbourhood scale nitrogen dioxide land use regression modelling with regression kriging in an urban transportation corridor, Atmos. Environ., 223, 117218, https://doi.org/10.1016/j.atmosenv.2019.117218, 2020b. 

Shi, W., Sun, Q., Du, P., Tang, S., Chen, C., Sun, Z., Wang, J., Li, T., and Shi, X.: Modification Effects of Temperature on the Ozone–Mortality Relationship: A Nationwide Multicounty Study in China, Environ. Sci. Technol., 54, 2859–2868, 2020c. 

Sun, Q., Hong, X., and Wold, L. E.: Cardiovascular effects of ambient particulate air pollution exposure, Circulation, 121, 2755–2765, 2010. 

Tong, C. H. M., Yim, S. H. L., Rothenberg, D., Wang, C., Lin, C. Y., Chen, Y. D., and Lau, N. C.: Assessing the impacts of seasonal and vertical atmospheric conditions on air auality over the Pearl River Delta Region, Atmos. Environ., 180, 69–78, 2018a. 

Tong, C. H. M., Yim, S. H. L., Rothenberg, D., Wang, C., Lin, C. Y., Chen, Y. D., and Lau, N. C.: Projecting the impacts of atmospheric conditions under climate change on air quality over Pearl River Delta region, Atmos. Environ., 193, 79–87, 2018b. 

TWEPA (Taiwan Environmental Protection Administration): Air Quality Annual Report of R. O. C. (Taiwan), available at: https://airtw.epa.gov.tw/CHT/Themes/LinkOut.aspx, last access: 9 July 2020. 

TWMOI (Taiwan Ministry of the Interior): Statistical Yearbook of Interior, available at: https://www.moi.gov.tw/files/site_stuff/321/2/year/year_en.html#2Population, last access: 9 July 2020. 

TWMOTC (Taiwan Ministry of Transportation and Communications): Annual Transportation Report, available at: https://www.motc.gov.tw/en/home.jsp?id=610&parentpath=0,154, last access: 9 July 2020. 

Vardoulakis, S., Solazzo, E., and Lumbreras, J.: Intra-urban and street scale variability of BTEX, NO2 and O3 in Birmingham, UK: Implications for exposure assessment, Atmos. Environ., 45, 5069–5078, 2011. 

Vienneau, D., De Hoogh, K., Beelen, R., Fischer, P., Hoek, G., and Briggs, D.: Comparison of land-use regression models between Great Britain and the Netherlands, Atmos. Environ., 44, 688–696, 2010. 

Wang, J., Cohan, D.S., and Xu, H.: Spatiotemporal ozone pollution LUR models: Suitable statistical algorithms and time scales for a megacity scale, Atmos. Environ., 237, 117671, 2020. 

Wang, M., Beelen, R., Bellander, T., Birk, M., Cesaroni, G., Cirach, M., Cyrys, J., de Hoogh, K., Declercq, C., Dimakopoulou, K., and Eeftens, M.: Performance of multi-city land use regression models for nitrogen dioxide and fine particles, Environ. Health Persp., 122, 843–849, 2014. 

Wang, M. Y., Yim, S. H., Wong, D. C., and Ho, K. F.: Source contributions of surface ozone in China using an adjoint sensitivity analysis, Sci. Total Environ., 662, 385–392, 2019. 

Wang, M. Y., Yim, S. H., Dong, G. H., Ho, K. F., and Wong, D. C.: Mapping ozone source-receptor relationship and apportioning the health impact in the Pearl River Delta region using adjoint sensitivity analysis, Atmos. Environ., 222, 117026, https://doi.org/10.1016/j.atmosenv.2019.117026, 2020. 

Wang, S. H., Huang, S. L., and Huang, P. J.: Can spatial planning really mitigate carbon dioxide emissions in urban areas? A case study in Taipei, Taiwan, Landscape Urban Plan., 169, 22–36, 2018. 

Weissert, L. F., Salmond, J. A., Miskell, G., Alavi-Shoshtari, M., and Williams, D. E.: Development of a microscale land use regression model for predicting NO2 concentrations at a heavy trafficked suburban area in Auckland, New Zealand, Sci. Total Environ., 619, 112–119, 2018. 

WHO (World Health Organization): WHO Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide – global update 2005, WHO/SDE/PHE/OEH/06.02, available at: https://www.who.int/airpollution/publications/aqg2005/en/ (last access: 9 July 2020), 2006. 

Wolf, K., Cyrys, J., Harciníková, T., Gu, J., Kusch, T., Hampel, R., Schneider, A., and Peters, A.: Land use regression modeling of ultrafine particles, ozone, nitrogen oxides and markers of particulate matter pollution in Augsburg, Germany, Sci. Total Environ., 579, 1531–1540, 2017. 

Wu, C. D., Chen, Y. C., Pan, W. C., Zeng, Y. T., Chen, M. J., Guo, Y. L., and Lung, S. C. C.: Land-use regression with long-term satellite-based greenness index and culture-specific sources to model PM2.5 spatial-temporal variability, Environ. Pollut., 224, 148–157, 2017. 

Wu, C. D., Zeng, Y. T., and Lung, S. C. C.: A hybrid kriging/land-use regression model to assess PM2.5 spatial-temporal variability, Sci. Total Environ., 645, 1456–1464, 2018. 

Xu, M., Sbihi, H., Pan, X., and Brauer, M.: Local variation of PM2.5 and NO2 concentrations within metropolitan Beijing, Atmos. Environ., 200, 254–263, 2019. 

Yang, Z., Freni-Sterrantino, A., Fuller, G. W., and Gulliver, J.: Development and transferability of ultrafine particle land use regression models in London, Sci. Total Environ., 740, 140059, https://doi.org/10.1016/j.scitotenv.2020.140059, 2020. 

Yim, S. H. L., Fung, J. C. H., Lau, A. K. H., and Kot, S. C.: Developing a high-resolution wind map for a complex terrain with a coupled MM5/CALMET system, J. Geophys. Res., 112, D05106, https://doi.org/10.1029/2006JD007752, 2007. 

Yim, S. H. L., Fung, J. C. H., and Ng, E. Y. Y.: An assessment indicator for air ventilation and pollutant dispersion potential in an urban canopy with complex natural terrain and significant wind variations, Atmos. Environ., 94, 297–306, 2014. 

Yim, S. H. L., Gu, Y., Shapiro, M. A., and Stephens, B.: Air quality and acid deposition impacts of local emissions and transboundary air pollution in Japan and South Korea, Atmos. Chem. Phys., 19, 13309–13323, https://doi.org/10.5194/acp-19-13309-2019, 2019a. 

Yim, S. H. L., Hou, X., Guo, J., and Yang, Y.: Contribution of local emissions and transboundary air pollution to air quality in Hong Kong during El Niño-Southern Oscillation and heatwaves, Atmos. Res., 218, 50–58, 2019b. 

Yim, S. H. L., Wang, M., Gu, Y., Yang, Y., Dong, G., and Li, Q.: Effect of urbanization on ozone and resultant health effects in the Pearl River Delta region of China, J. Geophys. Res.-Atmos., 124, 11568–11579, 2019c. 

Yin, P., Guo, J., Wang, L., Fan, W., Lu, F., Guo, M., Moreno, S. B., Wang, Y., Wang, H., Zhou, M., and Dong, Z.: Higher risk of cardiovascular disease associated with smaller size-fractioned particulate matter, Environ. Sci. Tech. Let., 7, 95–101, 2020. 

Yu, H. L. and Wang, C. H.: Retrospective prediction of intraurban spatiotemporal distribution of PM2.5 in Taipei, Atmos. Environ., 44, 3053–3065, 2010. 

Zhou, Y., Ma, J., Wang, B., Liu, Y., Xiao, L., Ye, Z., Fan, L., Wang, D., Mu, G., and Chen, W.: Long-term effect of personal PM2.5 exposure on lung function: A panel study in China, J. Hazard. Mater., 393, 122457, https://doi.org/10.1016/j.jhazmat.2020.122457, 2020.  

Zhu, D. and Zhou, X.: Effect of urban water bodies on distribution characteristics of particulate matters and NO2, Sustain. Cities Soc., 50, 101679, https://doi.org/10.1016/j.scs.2019.101679, 2019. 

Download
Short summary
This study established land-use regression (LUR) models using only routine air quality measurement data to support long-term health studies in an Asian metropolitan area. The established LUR models captured the spatial variability in exposure to air pollution with remarkable predictive accuracy. This is the first Asian study to evaluate intercity transferability of LUR models, and it highlights that there exist uncertainties when transferring LUR models between nearby cities.
Altmetrics
Final-revised paper
Preprint