Next Article in Journal
Board Diversity and Corporate Sustainability Performance: Do CEO Power and Firm Environmental Sensitivity Matter?
Next Article in Special Issue
Investigating Land Cover Changes and Their Impact on Land Surface Temperature in Khyber Pakhtunkhwa, Pakistan
Previous Article in Journal
Hybrid Integrated Computing Algorithm for Sustainable Tourism
Previous Article in Special Issue
Geochemical and Isotopic Evidence for Investigating the Impacts of Landfills on Groundwater: A Case Study in the Campania Region (Southern Italy)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relationships between Thermal Environment and Air Pollution of Seoul’s 25 Districts Using Vector Autoregressive Granger Causality

1
Korea Adaptation Center for Climate Change, Korea Environment Institute, Sejong 30116, Republic of Korea
2
Department of Urban Design and Planning, Hongik University, Seoul 04066, Republic of Korea
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(23), 16140; https://doi.org/10.3390/su152316140
Submission received: 30 September 2023 / Revised: 6 November 2023 / Accepted: 14 November 2023 / Published: 21 November 2023

Abstract

:
Rising temperatures and heightened air pollution are widespread across many parts of the world today. Despite some initial attempts for analysis, the intricate interconnection between the two still requires further investigation. This study focuses on Seoul, South Korea, by adopting vector-autoregressive-based Granger causality tests to unravel the nuances of these relationships at the district level. While bidirectional Granger causality links between temperature and urban heat island intensity, as well as between PM10 concentration and urban pollution island intensity, are found in many cases, our findings reveal diverse causal relationships that are evident in the districts. These findings underscore the necessity for evidence-based strategies to guide planners and policymakers in addressing the challenges of rising temperatures and air pollution in urban areas.

1. Introduction

Of the many environmental challenges South Korea faces today, rising temperatures and increasing air pollution levels are widely perceived by locals as one of the most critical [1,2]. The country’s overall temperature rose by up to 2.6 times faster than the global average over the past century [3], and air pollution levels of most of the country’s cities, usually represented among the locals by the concentration of particulate matter 10 μm or less in diameter (PM10), frequently exceed the World Health Organization (WHO)’s guidelines and remain higher than those of major global cities [4,5,6].
The consequential impacts of these challenges are tangibly evident. The country’s rising temperatures, especially during hot seasons, inflate the demand and consumption of energy [7,8], increase heat-related illnesses and deaths [9,10,11], and put its ecosystem and natural resources at higher risk. The high concentrations of PM10 exacerbate respiratory, cardiovascular, and skin diseases [12] and increase premature deaths [13].
To generate policy measures to tackle rising temperatures and air pollution levels, research has been identifying each of their causes. First, warming is caused by excessive emission of greenhouse gases, produced naturally or by anthropogenic activities, into the Earth’s atmosphere. They encompass reliance on fossil fuels for transportation, heating, cooling, and power generation, alongside urban land transformations catering to burgeoning infrastructural demands, particularly within densely populated cities [14,15,16,17]. Second, rising air pollution levels, especially the PM10 concentrations, highly depend on anthropogenic sources, including automobiles, industrial facilities, power plants, and buildings, in addition to natural sources like volcano eruptions, wildfires, and the decomposition of organic substances [18,19,20].
Notably, a larger body of literature finds that temperatures and concentrations of PM10 operate as a cause for each other in cities using statistical approaches in which usually independent and dependent variables are fixed. They identify temperature as a statistically significant or highly important explanatory variable when predicting PM10 concentrations [21,22,23,24,25,26,27], as heat increases energy use, thereby emitting more air pollutants and generating secondary aerosols in the atmosphere [28,29,30]. Another substantial body presents evidence that particulates induce discernible temperature fluctuations, often facilitated through radiative forcing mechanisms [31,32,33,34].
To address thermal environment and air pollution levels at the same time, some researchers use correlations to explain the relationship [35,36,37]. One German study, using an attribution method, finds that a highly complex interaction exists between temperatures and PM10 concentrations, as the former intensify turbulence mixing and the latter affect radiation transfer processes at the same time [38]. A study from China adopts Granger causality tests to reveal the intricate interactions between fine particles and urban heat islands at seasonal, monthly, and daily scales [39]. Another study from northeastern United States develops a vulnerability index to address their interaction [40]. However, it is difficult to find studies that attempt to explain the intricate causal relationships among multiple time series variables.
These studies commonly suggest that the relationship between temperatures and PM10 concentrations manifests in diverse ways. Yet, they also clearly demonstrate that it is difficult to converge to a firm agreement on the causal relationship between the two. Studies remark that existence, directions, and strengths of the relationships may significantly differ by spatial and temporal conditions set for each study and strongly call for further investigation across diverse contexts [35,39,41].
We look into Seoul, South Korea, which stands as the country’s capital city and also as one of the fastest warming and air-polluted cities in Far East Asia. We unravel the intricate causal relationships between the thermal environment and air pollution in Seoul’s 25 administrative districts, which form the bedrock of the city’s local policymaking. More specifically, we investigate temperatures and PM10 concentrations, as well as urban heat island intensity (UHII) and urban pollution island intensity (UPII), to yield deeper insights. We use data collected at each of the 25 districts’ weather and air monitoring stations from 1 January 2015 to 31 December 2019, the most recent five-year period prior to the outbreak of COVID-19, and employ a vector autoregressive (VAR)-based Granger causality test. Our findings may help to provide a strategic framework for crafting effective and targeted interventions tailored for each district. They can also serve as a valuable reference for urban regions worldwide that are grappling with similar challenges.

2. Materials and Methods

2.1. Case Context

As previously mentioned, Seoul (37°33′36″ N 126°59′24″ E) has long been experiencing rapid warming and a high concentration of PM10. Figure 1 reports that its annual average temperatures remained generally below 12 °C in the 1960s but rose sharply since then and now nears 14 °C. This corresponds to an increase rate of 0.3 °C per decade, which is more than twice the global average of 0.18 °C per decade since 1981 [42]. Figure 2 shows Seoul’s monthly PM10 concentration trends since 2015. Although there has been a gradual decline over time with seasonal fluctuations and a significant drop since early 2020, the concentrations still remain well above the WHO’s guidelines, 15 μg/m3. While the concentrations of most other major air pollutants like sulfur dioxide, nitrogen dioxide, carbon monoxide, and ozone remain low, PM10 stands out as a significant concern for local researchers and policymakers and is regarded as one of the few pollutants that still requires substantial improvement [43,44].

2.2. Variables and Data

Seoul, for a better accommodation of its 9.4 million residents and 5.7 million jobs on approximately 605 square kilometers of land, subdivides itself into 25 districts mostly based on population, as Figure 3 illustrates. Since the mid-1900s, Seoul’s neighborhoods have typically expanded in a way that broadly adheres to district borders. Today, districts ensure a relatively even provision of public service throughout Seoul and actively play a central role in developing significant planning and policy decisions. For this reason, we carry out analysis at the district level.
Our variables primarily fall into two categories: thermal environment and air pollution. For the former, we employ daily average temperatures and UHII, which is measured by the difference between the urban air temperature and the background rural temperature [45,46]. For the latter, we utilize daily average PM10 concentrations and UPII, which is determined by the difference in the concentration between urban areas and rural areas [41,47].
Data for the two thermal environment variables are sourced from the 29 weather stations situated in Seoul, as presented in Figure 4. There is at least one weather station in each district that makes it possible to assign relevant data to each district. Of the 29, 2 stations were excluded since they sit on high mountains or along rivers. For the two districts with two weather stations, we use averages of the two readings. We downloaded relevant data from the Open Met Data Portal (https://data.kma.go.kr/, accessed on 29 September 2023) of the Korea Meteorological Administration. Likewise, data for the two air pollution variables originate from the city’s 40 air quality monitoring stations, also illustrated in Figure 4, ensuring representation of at least 1 station in each district. We use averages for districts with two stations. Relevant data are downloaded from the Seoul Metropolitan Government’ Seoul Atmospheric Environment Information (https://cleanair.seoul.go.kr/, accessed on 29 September 2023).
While constructing the time series for daily average temperatures and PM10 concentrations, we encountered missing data points, constituting about 3% of the dataset. With this relatively minor percentage, we opted for linear interpolation as our preferred method, predicting values using neighboring data points.
By drawing from prior research [48,49,50], we compute the UHII and UPII by quantifying the difference between the average measurements taken within the confines of Seoul and those obtained at a reference location situated outside the city limits. We use the following formulae to calculate UHII and UPII:
U H I I i = T ¯ i T ¯ r
U P I I i = P ¯ i P ¯ r
where U H I I i is the UHII of location i , T ¯ i is the average temperature of location i , T ¯ r is the average temperature of the reference point r , U P I I i is the UPII of location i , P ¯ i is the average PM10 concentration of location i , and P ¯ r is the average PM10 concentration of the reference point.
From the selection of three potential reference points commonly employed by local researchers—Yangpyeong, Icheon, and Dongducheon—depicted in Figure 4, our preference aligns with Yangpyeong. It receives particular favor due to its distinctly rural characteristics, coupled with its alignment in terms of altitude and latitude with Seoul [51,52,53].

2.3. Analysis

Given the statistical properties of the variables used in this study, we adopt a VAR model approach. It is structured by relating each endogenous variable to the past values of all endogenous variables. This expansion of the univariate autoregressive model results in a vector autoregressive model encompassing multiple time series variables. Originally developed in the 1980s [54], the VAR model offers a straightforward methodology for analyzing and predicting multiple interconnected indicators where all the variables are treated as endogenous [55,56,57]. The VAR model for p -order lagged stationary time series data can be expressed as follows:
V A R p : y t = v + A 1 y t 1 + + A p y t p + u t
where v is the constant vector, A i is the regression parameter matrix of independent variables, and u t is the random error term.
In the context of time series data analysis, the underlying assumption is that the data exhibit stationarity, implying consistent statistical characteristics over time. Nonetheless, numerous real-world datasets deviate from this assumption and display non-stationarity. Such non-stationarity manifests as fluctuations in the mean or variance of the data across time, indicating the presence of a unit root, which can lead to spurious regressions. Therefore, it is necessary to verify whether the time series represents stationarity, which can be confirmed through a unit root test. We use the augmented Dickey–Fuller (ADF) test [58] as it is one of the most favored methods among a few. If a unit root exists, differencing is applied to data in order to remove any sort of stochasticity.
Once it is established that the dataset either possesses stationarity or has been transformed through differencing, the next step involves selecting a suitable lag for analysis. This selection is informed by criteria such as the Akaike information criterion (AIC), the Schwarz Bayesian information criterion (SBIC), and the Hannan–Quinn information criterion (HQIC). In this analysis, we employ all three of these criteria.
Granger causality test [59] is a statistical hypothesis that determines whether one time series is useful for predicting another time series. When the prediction of variable y is improved by including the past value of variable x , rather than solely relying on predicting y ’s current value using its own past values, x is said to Granger-cause y . We construct the two following regressions:
y t = i = 1 m α i x t i + i = 1 n β i y t i + u 1 t
x t = i = 1 p γ i x t i + i = 1 q δ i y t i + u 2 t
where α , β , γ , and δ are regression coefficients, and u t is the error terms.
As Figure 5 illustrates, the four variables construct six pairs of variables for six instances for conducting Granger causality tests: (1) daily average temperature and UHII, (2) daily average temperature and daily average PM10 concentration, (3) daily average temperature and UPII, (4) daily average PM10 concentration and UHII, (5) daily average PM10 concentration and UPII, and (6) PM10 concentration and UPII.
In VAR analysis, employing an impulse response function (IRF) is a typical procedure used to identify the influence of modifications in one endogenous variable on another endogenous variable usually in a graphical manner and to observe how the system dynamically reacts to shocks. To account for the immediate impacts of shocks, we employ the orthogonalized impulse response function (IRF), achieved through the utilization of Cholesky decomposition to transform the residual terms to become independent [60,61,62].
Forecast error variance decomposition (FEVD) is a post-analysis step in VAR models. It reveals how much variation in an endogenous variable can be attributed to specific shocks or impulses, separating their distinct effects. Like IRF, it entails orthogonalizing residual terms to determine each variable’s contribution to forecast error variance [63,64].

3. Results

Table 1 displays the descriptive statistics for the variables under analysis. Among the two thermal environment variables, the daily average temperature exhibits a mean of 13.5 and a standard deviation of 10.5. The mean and standard deviation for UHII are 0.9 and 1.4, respectively. Shifting to the air pollution variables, the mean of daily average PM10 concentration is 45.4, accompanied by a standard deviation of 27.7. The mean UPII is 1.8, with a standard deviation of 11.8.

3.1. Unit Root Test

Table 2 shows the ADF test results. Of the four variables in their original forms, we found that the average daily temperature exhibits unit roots in all 25 districts, while the remaining three showed no evidence of unit roots. To eliminate the seasonality present in the temperature data, seasonal differencing was applied to all variables to ensure normality and establish a stable time series for further analysis. As a result, seasonally differenced data successfully pass unit root tests in all cases, as shown in the table.

3.2. Lag Length Selection

We employed AIC, SBIC, and HQIC statistics to ascertain the optimal lag length for each of the 25 districts, taking into consideration that the causal effects may occur gradually and manifest in changes later in time in VAR models. As illustrated in Table 3, five districts, namely Dobong, Gangbuk, Jongro, Jungrang, and Songpa, exhibited a lag length of 5, while that of the remaining twenty was determined to be 6.

3.3. Granger Causality Tests

Using stationary time series data and optimal lag lengths, Granger causality tests were performed for all 25 districts. These tests aimed to reveal causal relationships among four variables: daily average temperature, UHII, daily average PM10 concentration, and UPII. The results for each district are available in Appendix A due to their volume. For brevity, Table 4 summarizes twelve distinct Granger-causal relationships.
Type 1 comprises the districts of Geumcheon, Mapo, Yangcheon, and Yeongdeungpo, each exhibiting the following characteristics: (1) daily average temperature Granger-causes UHII, daily average PM10 concentration, and UPII; (2) UHII Granger-causes daily average temperature, daily average PM10 concentration, and UPII; (3) daily average PM10 concentration Granger-causes UPII; and (4) UPII Granger-causes daily average temperature and daily average PM10 concentration. Type 2 includes Gangdong, Jungrang, Songpa, and Yongsan districts and yields the following causalities: (1) daily average temperature Granger-causes UHII and daily average PM10 concentration; (2) UHII Granger-causes daily average temperature and UPII; (3) daily average PM10 concentration Granger-causes UPII; and (4) UPII Granger-causes daily average PM10 concentration. Dongdaemun, Gangbuk, Jung, and Seongbuk districts constitute Type 3, displaying the following relationships: (1) daily average temperature Granger-causes UHII and daily average PM10 concentration; (2) UHII Granger-causes daily average temperature, daily average PM10 concentration, and UPII; (3) daily average PM10 concentration Granger-causes UPII; and (4) UPII Granger-causes daily average PM10 concentration and UHII. Type 4 is represented by the districts of Gangnam, Gwangjin, and Seongdong, all of which share the following: (1) daily average temperature Granger-causes UHII, daily average PM10 concentration, and UPII; (2) UHII Granger-causes daily average temperature, daily average PM10 concentration, and UPII; (3) daily average PM10 concentration Granger-causes UPII; and (4) UPII Granger-causes daily average PM10 concentration.
Type 5 is composed of Dobong and Gangseo districts, which present the following: (1) daily average temperature Granger-causes UHII and daily average PM10 concentration; (2) UHII Granger-causes daily average temperature, daily average PM10 concentration, and UPII; (3) daily average PM10 concentration Granger-causes daily average temperature and UPII; and (4) UPII Granger-causes daily average PM10 concentration. Type 6 consists of Seocho and Seodaemun districts, which share the following: (1) daily average temperature Granger-causes UHII and daily average PM10 concentration; (2) UHII Granger-causes daily average temperature, daily average PM10 concentration, and UPII; (3) daily average PM10 concentration Granger-causes UPII; and (4) UPII Granger-causes daily average temperature and daily average PM10 concentration. Types 7 through 12 are represented by the districts of Dongjak, Eunpyeong, Guro, Gwanak, Jongro, and Nowon, respectively. Like the previous six, they also yield distinctive mixes of Granger causality pairs.

3.4. IRFs and FEVD Results

Table 5 displays IRFs for the Granger-causal relationships within this analysis. It highlights a selected subset of exemplary functions, each exemplifying distinctive patterns, from a total of 400 functions, with 16 functions dedicated to each district. Notably, the table excludes the scenario in which the impulse variable is pm10 and the response variable is uhii, as no causal relationship was identified in that case.
To elaborate on a few, firstly, when the impulse variable is temp and the response variable is uhii, as shown in the first row of the table, the effect of a one-standard-deviation shock in daily average temperature becomes apparent. Initially, UHII exhibits a positive response, starting at 0 in period 1. A subsequent negative reaction occurs between periods 3 and 4, followed by a gradual convergence toward 0. Predominantly positive trends persist until later periods, with a significant positive reaction around the sixth or seventh period. However, despite this, the impact diminishes after 4 days, as indicated by the confidence interval encompassing 0. This pattern in Gangdong district aligns with similar trends observed in other districts. Second, when the impulse variable is uhii and the response variable is upii, the effect of a one-standard-deviation shock in UHII yields impacts on UPUII in three distinctive patterns. In the Gangdong case, an initial positive response in the first period was swiftly countered by a negative reaction in the second period. Stability marked periods 3 to 5, with short fluctuations in periods 6 and 7. Other periods showed no notable response. Seodaemun exhibited a similar but less significant pattern up to the fifth period, with significant reactions in the sixth and seventh stages. Nowon displayed an initial positive response, followed by alternating negative and positive reactions. Statistically significant responses were observed in the first, third, and seventh stages, with no substantial reactions in other periods.
Overall, the analysis revealed that all the IRFs exhibited immediate shocks in the early stages, typically up to the seventh period, but then gradually converged to 0 or reached statistically insignificant levels of response.
As indicated in Table 6, which presents the FEVD results for all relationships, during the 10th period, a significant portion of the variance in each of the four variables can be attributed to their own past values. Notably, daily average temperature demonstrates the highest degree of self-explanation, ranging from 93.79% to 97.73%. Following closely is daily average PM10 concentration, accounting for 86.05% to 92.68% of its own variance. UPII explains 67.45% to 87.95% of its variance, while UHII shows the lowest self-explanation, ranging from 58.96% to 85.12%. This suggests that the two intensity variables, UHII and UPII, are relatively more influenced by their own past values than by other variables. UHII and UPII are also somewhat affected by daily average temperature and daily average PM10 concentration, respectively. For a more compact presentation, the FEVD results for each district by each causal relationship are presented in Appendix B.

4. Discussion

What has been presented thus far warrants in-depth discussion. First, in terms of the thermal environment, our analysis has revealed a bidirectional Granger causality between daily average temperature and UHII in most cases. The IRFs showed that, with the exception of the Nowon district, a substantial positive response was exhibited when subjected to shocks in daily average temperature and UHII. The FEVD results emphasize that beyond their inherent influence, daily average temperature was significantly influenced by UHII, and UHII was predominantly shaped by daily average temperature across all districts. Second, concerning the air pollution variables, daily average PM10 concentration and UPII, Granger causality was observed in both directions throughout Seoul, with an initial positive response in each case. Third, in districts categorized as types 1 and 2, our analysis indicates that the thermal environment Granger-caused air pollution in one direction. The IRFs revealed that the air pollution variables initially exhibited a notable positive response to shocks in the thermal environment, signifying the rapid deterioration of air pollution due to changes in the thermal environment. The FEVD results indicate that the overall influence of temperature on air pollution variables exceeded that of UHII. Fourth, in the majority of districts in Seoul, specifically types 3 to 12, a complex interaction between the thermal environment and air pollution was observed. The IRF analysis reveals an initial positive response, indicating that air pollution worsens due to changes in the thermal environment and vice versa. However, the response of the thermal environment to UPII shocks generally displayed an initial negative reaction. The FEVD results indicate that, in general, the thermal environment has a greater influence on air pollution than the reverse.
Figure 6 presents a comprehensive summary. (1) Daily average temperature and UHII display two-way Granger-causal relationships in 24 districts, with only one district showing daily average temperature as the Granger cause of UHII. (2) Daily average PM10 concentration and UPII reveal two-way Granger-causal relationships across all 25 districts. (3) Daily average temperature and daily average PM10 concentration exhibit two-way Granger-causal relationships in 6 districts. Meanwhile, daily average temperature Granger-causes daily average PM10 concentration in 19 districts. (4) UHII and UPII establish two-way Granger-causal relationships in 5 districts, with UHII acting as the Granger cause of UPII in 20 districts. (5) Daily average temperature and UPII demonstrate two-way Granger-causal relationships in five districts. In four districts, daily average temperature Granger-causes UPII, while UPII Granger-causes daily average temperature in five districts. (6) Daily average PM10 concentration and UHII do not exhibit any two-way Granger-causal relationships. However, daily average PM10 concentration Granger-causes UHII in 4 districts, while UHII Granger-causes daily average PM10 concentration in 16 districts.
Our findings are generally in line with what previous studies set forth. The one-way causal relationships, where either thermal environment is the independent variable and air pollution is the dependent [21,22,23,24,25,26,27,28,29,30] or the other way around [31,32,33,34], are evident in all districts of Seoul. At the same time, each district yields at least one two-way causal relationship, as another body of literature suggests [35,36,37,38,39,40].
Our contribution lies in moving beyond surface-level observations. Unlike prior studies, we reveal the coexistence of both one-way and two-way causal relationships within Seoul. Furthermore, our findings highlight the complexity of these relationships across the city’s 25 districts, indicating diverse combinations of causal interactions. This nuanced understanding underscores the variability of causal relationships within a city.

5. Conclusions

This study investigated the relationship between the thermal environment and air pollution in Seoul’s 25 districts, conducting a VAR Granger causality test to figure out causal relationships and IRF and FEVD to discover dynamic relationships. Major findings of this study can be summarized as follows. There exists a bidirectional Granger causality relationship between two thermal environment variables, namely the daily average temperature and UHII, in nearly all districts examined. Similarly, a bidirectional relationship is evident between the two air pollution variables, daily average PM10 concentration and UPII, across all 25 districts. Simultaneously, it is noteworthy that a greater number of districts in general exhibit a Granger causality relationship from thermal environment variables to air pollution variables.
The conclusion of this study yields several policy recommendations. First is the concurrent management of temperature and UHIs. Our findings support the idea of simultaneously addressing temperature and UHII. While many climate action plans focus on reducing temperature through various adaptation strategies, they often overlook temperature variations across urban regions. Second, the simultaneous management of PM10 concentration and UPII is crucial. Current air pollution policies have primarily concentrated on reducing air pollutants, neglecting the concept of urban pollution islands. Therefore, alongside existing air quality improvement goals, it is imperative to establish air quality standards that account for differences in urban and suburban air quality levels. Third, especially in most districts where one-way Granger-causal relationships exist, with thermal environment variables influencing air pollution variables but not vice versa, more emphasis on the former is crucial. Lastly, this holistic approach to managing the thermal environment and air pollution should be guided by the region-wide real-time monitoring of temperature and air pollution levels so as to provide an evidence-based roadmap for policymakers.
This study acknowledges several limitations. First, although we have uncovered complex links between the thermal environment and air pollution in Seoul’s 25 districts, the use of a VAR Granger causality test fails to clearly present specific factors that contribute to the diverse causal relationships. Second, although we are using five years of time series data, results from the statistical models do not capture relationships that may change over time. Third, the study’s reliance on pre-COVID-19 data means that it does not account for potential shifts in these relationships post-COVID-19. Lastly, our data are confined within the spatial boundaries of Seoul and exploring data beyond the city may yield varied results, potentially leading to different insights and recommendations.
Despite its limitations, this study represents one of the pioneering endeavors to comprehensively unveil the intricate interplay between the thermal environment and air pollution. The insights derived from Seoul have the potential to significantly aid other regions worldwide that are grappling with comparable environmental issues. These findings can serve as valuable guidance for communities striving to create healthier and more sustainable urban environments globally.

Author Contributions

Conceptualization, J.Y. and H.K.; methodology, J.Y., H.K. and J.L.; software, H.K. and J.L.; formal analysis, J.Y., H.K. and J.L.; investigation, J.Y. and H.K.; resources, J.Y.; data curation, J.Y.; writing—original draft preparation, J.Y. and H.K.; writing—review and editing, H.K. and J.L.; visualization, J.Y. and H.K.; supervision, H.K.; funding acquisition, H.K. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Korea Agency for Infrastructure Technology Advancement (KAIA) grant funded by the Ministry of Land, Infrastructure and Transport (Grant RS-2022-00143404). This work is supported by the Korea Agency for Infrastructure Technology Advancement (KAIA) grant funded by the Ministry of Land, Infrastructure and Transport (Grant RS-2021-KA162410).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors thank the anonymous reviewers for their excellent comments.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Table A1. Granger causality test results between temp and uhii.
Table A1. Granger causality test results between temp and uhii.
District χ 2
H0: temp Does Not Granger-Cause uhiiH0: uhii Does Not Granger-Cause temp
Dobong96.915 ***68.582 ***
Dongdaemun94.680 ***121.470 ***
Dongjak110.530 ***45.903 ***
Eunpyeong35.205 ***10.805 ***
Gangbuk106.890 ***115.870 ***
Gangdong99.279 ***105.600 ***
Gangnam88.008 ***107.100 ***
Gangseo88.260 ***17.048 **
Geumcheon81.065 ***73.020 ***
Guro87.597 ***58.479 ***
Gwanak82.853 ***33.236 ***
Gwangjin94.715 ***104.050 ***
Jongro96.569 ***102.390 ***
Jung87.518 ***91.097 ***
Jungrang78.183 ***114.230 ***
Mapo91.687 ***73.281 ***
Nowon41.358 ***25.353 ***
Seocho101.490 ***86.658 ***
Seodaemun100.610 ***69.943 ***
Seongbuk85.126 ***97.699 ***
Seongdong111.600 ***110.940 ***
Songpa106.770 ***121.750 ***
Yangcheon94.990 ***81.147 ***
Yeongdeungpo83.670 ***85.782 ***
Yongsan79.787 ***83.881 ***
** p < 0.01, *** p < 0.001.
Table A2. Granger causality test results between pm10 and upii.
Table A2. Granger causality test results between pm10 and upii.
District χ 2
H0: pm10 Does Not Granger-Cause upiiH0: upii Does Not Granger-Cause pm10
Dobong46.778 ***54.926 ***
Dongdaemun53.082 ***80.225 ***
Dongjak48.955 ***58.427 ***
Eunpyeong79.051 ***62.192 ***
Gangbuk68.752 ***61.114 ***
Gangdong43.302 ***58.150 ***
Gangnam37.563 ***67.334 ***
Gangseo25.130 ***51.894 ***
Geumcheon96.954 ***59.293 ***
Guro42.355 ***60.106 ***
Gwanak55.410 ***66.740 ***
Gwangjin50.616 ***50.857 ***
Jongro86.934 ***94.985 ***
Jung81.111 ***65.221 ***
Jungrang75.393 ***77.423 ***
Mapo40.888 ***61.009 ***
Nowon30.332 ***52.086 ***
Seocho44.565 ***32.921 ***
Seodaemun65.699 ***77.948 ***
Seongbuk48.804 ***56.387 ***
Seongdong59.608 ***71.646 ***
Songpa95.204 ***62.781 ***
Yangcheon45.908 ***55.287 ***
Yeongdeungpo29.768 ***59.574 ***
Yongsan67.606 ***64.298 ***
*** p < 0.001.
Table A3. Granger causality test results between temp and pm10.
Table A3. Granger causality test results between temp and pm10.
District χ 2
H0: temp Does Not Granger-Cause pm10H0: pm10 Does Not Granger-Cause temp
Dobong12.341 *15.085 **
Dongdaemun24.603 ***7.665
Dongjak29.509 ***12.884 *
Eunpyeong36.295 ***7.158
Gangbuk15.717 **9.262
Gangdong23.579 **11.097
Gangnam21.831 **7.690
Gangseo27.268 ***16.201 **
Geumcheon28.610 ***9.398
Guro31.106 ***12.635 *
Gwanak29.126 ***13.249 *
Gwangjin27.984 ***9.163
Jongro13.036 *8.252
Jung21.159 **5.125
Jungrang17.742 **6.693
Mapo27.330 ***11.495
Nowon33.256 ***21.092 **
Seocho24.037 **9.521
Seodaemun23.466 **9.901
Seongbuk22.736 **11.084
Seongdong28.668 ***8.392
Songpa13.313 *7.478
Yangcheon30.191 ***7.817
Yeongdeungpo28.869 ***10.831
Yongsan22.492 **8.159
* p < 0.05, ** p < 0.01, *** p < 0.001.
Table A4. Granger causality test results between temp and upii.
Table A4. Granger causality test results between temp and upii.
District χ 2
H0: temp Does Not Granger-Cause upiiH0: upii Does Not Granger-Cause temp
Dobong2.7516.929
Dongdaemun10.4386.050
Dongjak14.393 *6.309
Eunpyeong9.18920.958 **
Gangbuk4.0715.890
Gangdong10.6375.436
Gangnam13.573 *10.603
Gangseo10.79311.522
Geumcheon14.250 *13.487 *
Guro16.554 *14.845 *
Gwanak11.2535.342
Gwangjin16.267 *7.908
Jongro6.01415.156 *
Jung9.7639.022
Jungrang8.5477.445
Mapo19.389 **14.407 *
Nowon12.52113.704 *
Seocho10.06312.712 *
Seodaemun10.42418.683 **
Seongbuk7.03511.483
Seongdong17.917 **6.8461
Songpa8.1976.428
Yangcheon14.606 *15.362 *
Yeongdeungpo22.464 **12.795 *
Yongsan12.07711.077
* p < 0.05, ** p < 0.01.
Table A5. Granger causality test results between pm10 and uhii.
Table A5. Granger causality test results between pm10 and uhii.
District χ 2
H0: pm10 Does Not Granger-Cause uhiiH0: uhii Does Not Granger-Cause pm10
Dobong5.61916.130 **
Dongdaemun6.90738.797 ***
Dongjak9.8259.300
Eunpyeong6.4036.532
Gangbuk4.97032.701 ***
Gangdong6.39534.032 ***
Gangnam7.97435.918 ***
Gangseo7.45413.849 *
Geumcheon7.49813.037 *
Guro6.4269.684
Gwanak7.9847.869
Gwangjin8.16036.288 ***
Jongro5.02519.412 **
Jung4.02131.833 ***
Jungrang5.20822.893 ***
Mapo9.14420.616 **
Nowon2.62310.725
Seocho11.00818.450 **
Seodaemun8.58415.009 *
Seongbuk9.61218.942 **
Seongdong9.63024.130 ***
Songpa5.38715.492 **
Yangcheon8.92420.153 **
Yeongdeungpo10.02418.994 **
Yongsan6.89920.797 **
* p < 0.05, ** p < 0.01, *** p < 0.001.
Table A6. Granger causality test results between uhii and upii.
Table A6. Granger causality test results between uhii and upii.
District χ 2
H0: uhii does Not Granger-Cause upiiH0: upii Does Not Granger-Cause uhii
Dobong12.639 *6.067
Dongdaemun44.503 ***12.788 ***
Dongjak17.386 **5.861
Eunpyeong14.831 *7.908
Gangbuk35.170 ***12.682 *
Gangdong38.899 ***7.907
Gangnam38.273 ***7.328
Gangseo17.048 **5.201
Geumcheon20.021 **7.723
Guro16.247 **8.526
Gwanak13.784 **2.980
Gwangjin32.207 ***12.273
Jongro29.905 ***21.792 **
Jung25.207 ***12.987 *
Jungrang21.357 **9.580
Mapo31.171 ***6.564
Nowon16.341 *6.151
Seocho27.478 ***7.369
Seodaemun14.445 *8.091
Seongbuk30.854 ***17.183 **
Seongdong34.680 ***4.935
Songpa18.350 **2.557
Yangcheon29.573 ***12.446
Yeongdeungpo33.116 ***12.533
Yongsan31.187 ***11.566
* p < 0.05, ** p < 0.01, *** p < 0.001.

Appendix B

Table A7. FEVD results when the response variable is temp.
Table A7. FEVD results when the response variable is temp.
Districttemptemppm10tempuhiitempupiitemp
Dobong0.95960.00740.02870.0044
Dongdaemun0.93790.00520.05240.0045
Dongjak0.96940.00740.01950.0037
Eunpyeong0.97730.00450.00680.0115
Gangbuk0.94110.00610.04810.0047
Gangdong0.94040.00790.04870.0030
Gangnam0.94310.00480.04570.0065
Gangseo0.96170.00660.02710.0046
Geumcheon0.95510.00510.03190.0079
Guro0.96000.00570.02610.0082
Gwanak0.97390.00780.01530.0030
Gwangjin0.94400.00550.04560.0049
Jongro0.94480.00470.04030.0103
Jung0.95090.00390.03890.0063
Jungrang0.94000.00520.04980.0049
Mapo0.95520.00570.03120.0079
Nowon0.96740.01110.01550.0060
Seocho0.94950.00590.03750.0071
Seodaemun0.95200.00590.03190.0102
Seongbuk0.94600.00510.04230.0065
Seongdong0.94250.00610.04750.0039
Songpa0.94110.00660.04900.0033
Yangcheon0.95200.00350.03520.0094
Yeongdeungpo0.95280.00410.03560.0075
Yongsan0.95300.00440.03590.0068
Range0.9379~0.97730.0035~0.01110.0068~0.05240.0030~0.0115
Table A8. FEVD results when the response variable is uhii.
Table A8. FEVD results when the response variable is uhii.
Districttempuhiipm10uhiiuhiiuhiiupiiuhii
Dobong0.22800.00410.76490.003
Dongdaemun0.20910.00770.77630.0069
Dongjak0.29740.00420.69620.0022
Eunpyeong0.36820.01050.61610.0052
Gangbuk0.25530.00360.73390.0072
Gangdong0.16800.01160.81640.0040
Gangnam0.26320.00440.72790.0044
Gangseo0.28330.00610.70900.0016
Geumcheon0.33980.00460.65170.0040
Guro0.30090.00570.69090.0026
Gwanak0.33550.00440.65940.0007
Gwangjin0.20620.00650.78050.0069
Jongro0.32940.00240.65710.0111
Jung0.40040.00400.58960.006
Jungrang0.22110.00400.76850.0064
Mapo0.39780.00450.59490.0028
Nowon0.14080.00380.85120.0042
Seocho0.32050.00400.67260.0029
Seodaemun0.39000.00410.60250.0035
Seongbuk0.29300.00480.69400.0081
Seongdong0.23510.00520.75700.0028
Songpa0.25020.00480.74360.0014
Yangcheon0.31090.00430.67830.0065
Yeongdeungpo0.31440.00310.67620.0063
Yongsan0.29950.00230.69290.0053
Range0.1408~0.40040.0023~0.01160.5896~0.85120.0007~0.0111
Table A9. FEVD results when the response variable is pm10.
Table A9. FEVD results when the response variable is pm10.
Districttemppm10pm10pm10uhiipm10upiipm10
Dobong0.04530.91630.00880.0296
Dongdaemun0.08160.86050.01380.0441
Dongjak0.07050.89440.00370.0314
Eunpyeong0.03660.92680.00350.0331
Gangbuk0.05420.90110.01300.0316
Gangdong0.06520.88800.01460.0322
Gangnam0.06490.88400.01270.0383
Gangseo0.07030.89460.00560.0295
Geumcheon0.05490.90770.00520.0322
Guro0.06380.89950.00410.0326
Gwanak0.07380.88820.00350.0345
Gwangjin0.05980.89760.01280.0299
Jongro0.05810.88270.00770.0516
Jung0.07360.87770.01100.0376
Jungrang0.04940.90050.00810.0420
Mapo0.06550.89420.00700.0333
Nowon0.06120.90250.00850.0279
Seocho0.06880.90490.00790.0184
Seodaemun0.04690.90650.00510.0415
Seongbuk0.05730.90170.00840.0326
Seongdong0.07200.88000.00970.0384
Songpa0.04930.91050.00850.0316
Yangcheon0.07400.88830.00780.0299
Yeongdeungpo0.07470.88470.00690.0338
Yongsan0.06650.88870.00760.0372
Range0.0366~0.08160.8605~0.92680.0035~0.01460.0184~0.0516
Table A10. FEVD results when the response variable is upii.
Table A10. FEVD results when the response variable is upii.
Districttempupiipm10upiiuhiiupiiupiiupii
Dobong0.02280.23880.00780.7306
Dongdaemun0.03190.06770.02480.8756
Dongjak0.03550.15410.00910.8013
Eunpyeong0.0110.24180.00790.7393
Gangbuk0.0190.1850.01870.7773
Gangdong0.03010.07160.02360.8748
Gangnam0.02660.1880.02160.7639
Gangseo0.04860.19540.0110.7451
Geumcheon0.01220.13730.01170.8388
Guro0.02960.19470.00940.7663
Gwanak0.03670.14120.0080.8141
Gwangjin0.02340.20310.01820.7553
Jongro0.02340.17140.01230.7929
Jung0.02240.08890.01410.8747
Jungrang0.01130.22390.01260.7522
Mapo0.03370.23030.01550.7205
Nowon0.04890.0590.01260.8795
Seocho0.03750.20840.01450.7395
Seodaemun0.01020.25070.00880.7303
Seongbuk0.02620.19660.01680.7604
Seongdong0.0330.15510.01860.7933
Songpa0.0170.15850.01210.8124
Yangcheon0.02360.14240.01680.8172
Yeongdeungpo0.04830.26170.01550.6745
Yongsan0.02990.17810.01530.7768
Range0.0102~0.04890.0590~0.26170.0078~0.02480.6745~0.8795

References

  1. Son, G.; Hwang, W. Consumers’ Views on ESG and Environmental Friendly Consumption Behavior; KB Financial Group: Seoul, Republic of Korea, 2021. [Google Scholar]
  2. Hwang, I.; Baek, J. Market Mechanisms for Addressing Climate Change and Air Pollution in Seoul; Seoul Institute: Seoul, Republic of Korea, 2020. [Google Scholar]
  3. Korea Meteorological Administration. Korean Climate Change Assessment Report 2020: Scientific Basis for Climate Change (Summary for Policymakers); Korea Meteorological Administration: Seoul, Republic of Korea, 2020.
  4. OECD. How’s Life? 2020: Measuring Well-Being. Available online: https://www.oecd.org/wise/how-s-life-23089679.htm (accessed on 22 August 2023).
  5. Kim, H. Seasonal Impacts of Particulate Matter Levels on Bike Sharing in Seoul, South Korea. Int. J. Environ. Res. Public Health 2020, 17, 3999. [Google Scholar] [CrossRef] [PubMed]
  6. Kim, H. Land Use Impacts on Particulate Matter Levels in Seoul, South Korea: Comparing High and Low Seasons. Land 2020, 9, 142. [Google Scholar] [CrossRef]
  7. Shin, J.; Yang, H.; Kim, C. The Relationship between Climate and Energy Consumption: The Case of South Korea. Energy Sources Part A Recovery Util. Environ. Eff. 2023, 45, 6456–6471. [Google Scholar] [CrossRef]
  8. Lee, K.; Baek, H.-J.; Cho, C. The Estimation of Base Temperature for Heating and Cooling Degree-Days for South Korea. J. Appl. Meteorol. Climatol. 2014, 53, 300–309. [Google Scholar] [CrossRef]
  9. Lim, Y.-H.; Lee, K.-S.; Bae, H.-J.; Kim, D.; Yoo, H.; Park, S.; Hong, Y.-C. Estimation of Heat-Related Deaths during Heat Wave Episodes in South Korea (2006–2017). Int. J. Biometeorol. 2019, 63, 1621–1629. [Google Scholar] [CrossRef] [PubMed]
  10. Park, J.; Chae, Y.; Choi, S.H. Analysis of Mortality Change Rate from Temperature in Summer by Age, Occupation, Household Type, and Chronic Diseases in 229 Korean Municipalities from 2007–2016. Int. J. Environ. Res. Public Health 2019, 16, 1561. [Google Scholar] [CrossRef] [PubMed]
  11. Kim, J.; Song, K.J.; Hong, K.J.; Ro, Y.S. Trend of Outbreak of Thermal Illness Patients Based on Temperature 2002–2013 in Korea. Climate 2017, 5, 94. [Google Scholar] [CrossRef]
  12. Kim, O.-J.; Lee, S.H.; Kang, S.-H.; Kim, S.-Y. Incident Cardiovascular Disease and Particulate Matter Air Pollution in South Korea Using a Population-Based and Nationwide Cohort of 0.2 Million Adults. Environ. Health 2020, 19, 113. [Google Scholar] [CrossRef]
  13. Kim, J.-H.; Oh, I.-H.; Park, J.-H.; Cheong, H.-K. Premature Deaths Attributable to Long-Term Exposure to Ambient Fine Particulate Matter in the Republic of Korea. J. Korean Med. Sci. 2018, 33, e251. [Google Scholar] [CrossRef]
  14. Wheeler, S.M. Planning for Sustainability: Creating Livable, Equitable and Ecological Communities; Routledge: New York, NY, USA, 2013; ISBN 978-0-415-80989-4. [Google Scholar]
  15. Newman, P.; Beatley, T.; Boyer, H. Resilient Cities: Responding to Peak Oil and Climate Change; Island Press: Washington, DC, USA, 2009; ISBN 978-1-59726-498-3. [Google Scholar]
  16. Newman, P.; Kenworthy, J. Sustainability and Cities: Overcoming Automobile Dependence; Island Press: Washington, DC, USA, 1999; ISBN 978-1-55963-660-5. [Google Scholar]
  17. Calthorpe, P. Urbanism in the Age of Climate Change, 2nd ed.; Island Press: Washington, DC, USA, 2013; ISBN 978-1-59726-721-2. [Google Scholar]
  18. Moussiopoulos, N. Air Quality in Cities; Springer: Berlin/Heidelberg, Germany, 2003; ISBN 978-3-642-05646-8. [Google Scholar]
  19. Zhang, D.; Zhou, C.; He, B.-J. Spatial and Temporal Heterogeneity of Urban Land Area and PM2.5 Concentration in China. Urban Clim. 2022, 45, 101268. [Google Scholar] [CrossRef]
  20. Ross, Z.; Jerrett, M.; Ito, K.; Tempalski, B.; Thurston, G.D. A Land Use Regression for Predicting Fine Particulate Matter Concentrations in the New York City Region. Atmos. Environ. 2007, 41, 2255–2269. [Google Scholar] [CrossRef]
  21. Stafoggia, M.; Bellander, T.; Bucci, S.; Davoli, M.; de Hoogh, K.; de’ Donato, F.; Gariazzo, C.; Lyapustin, A.; Michelozzi, P.; Renzi, M.; et al. Estimation of Daily PM10 and PM2.5 Concentrations in Italy, 2013–2015, Using a Spatiotemporal Land-Use Random-Forest Model. Environ. Int. 2019, 124, 170–179. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, Z.; Wang, J.; Hart, J.E.; Laden, F.; Zhao, C.; Li, T.; Zheng, P.; Li, D.; Ye, Z.; Chen, K. National Scale Spatiotemporal Land-Use Regression Model for PM2.5, PM10 and NO2 Concentration in China. Atmos. Environ. 2018, 192, 48–54. [Google Scholar] [CrossRef]
  23. Kim, H.; Hong, S. Relationship between Land-Use Type and Daily Concentration and Variability of PM10 in Metropolitan Cities: Evidence from South Korea. Land 2022, 11, 23. [Google Scholar] [CrossRef]
  24. Ahn, H.; Lee, J.; Hong, A. Urban Form and Air Pollution: Clustering Patterns of Urban Form Factors Related to Particulate Matter in Seoul, Korea. Sustain. Cities Soc. 2022, 81, 103859. [Google Scholar] [CrossRef]
  25. Park, Y.; Shin, J.; Lee, J.Y. Spatial Association of Urban Form and Particulate Matter. Int. J. Environ. Res. Public Health 2021, 18, 9428. [Google Scholar] [CrossRef]
  26. McCarty, J.; Kaza, N. Urban Form and Air Quality in the United States. Landsc. Urban Plan. 2015, 139, 168–179. [Google Scholar] [CrossRef]
  27. Sarrat, C.; Lemonsu, A.; Masson, V.; Guedalia, D. Impact of Urban Heat Island on Regional Atmospheric Pollution. Atmos. Environ. 2006, 40, 1743–1758. [Google Scholar] [CrossRef]
  28. Akbari, H.; Pomerantz, M.; Taha, H. Cool Surfaces and Shade Trees to Reduce Energy Use and Improve Air Quality in Urban Areas. Sol. Energy 2001, 70, 295–310. [Google Scholar] [CrossRef]
  29. Akbari, H. Shade Trees Reduce Building Energy Use and CO2 Emissions from Power Plants. Environ. Pollut. 2002, 116, S119–S126. [Google Scholar] [CrossRef]
  30. Stone, B. Urban Heat and Air Pollution: An Emerging Role for Planners in the Climate Change Debate. J. Am. Plan. Assoc. 2005, 71, 13–25. [Google Scholar] [CrossRef]
  31. Jin, M.; Shepherd, J.M.; Zheng, W. Urban Surface Temperature Reduction via the Urban Aerosol Direct Effect: A Remote Sensing and WRF Model Sensitivity Study. Adv. Meteorol. 2011, 2010, e681587. [Google Scholar] [CrossRef]
  32. Jin, M.S.; Kessomkiat, W.; Pereira, G. Satellite-Observed Urbanization Characters in Shanghai, China: Aerosols, Urban Heat Island Effect, and Land–Atmosphere Interactions. Remote Sens. 2011, 3, 83–99. [Google Scholar] [CrossRef]
  33. Li, H.; Sodoudi, S.; Liu, J.; Tao, W. Temporal Variation of Urban Aerosol Pollution Island and Its Relationship with Urban Heat Island. Atmos. Res. 2020, 241, 104957. [Google Scholar] [CrossRef]
  34. Cao, C.; Lee, X.; Liu, S.; Schultz, N.; Xiao, W.; Zhang, M.; Zhao, L. Urban Heat Islands in China Enhanced by Haze Pollution. Nat. Commun. 2016, 7, 12509. [Google Scholar] [CrossRef] [PubMed]
  35. Ngarambe, J.; Joen, S.J.; Han, C.-H.; Yun, G.Y. Exploring the Relationship between Particulate Matter, CO, SO2, NO2, O3 and Urban Heat Island in Seoul, Korea. J. Hazard. Mater. 2021, 403, 123615. [Google Scholar] [CrossRef] [PubMed]
  36. Lai, L.-W.; Cheng, W.-L. Air Quality Influenced by Urban Heat Island Coupled with Synoptic Weather Patterns. Sci. Total Environ. 2009, 407, 2724–2733. [Google Scholar] [CrossRef]
  37. Xu, L.Y.; Xie, X.D.; Li, S. Correlation Analysis of the Urban Heat Island Effect and the Spatial and Temporal Distribution of Atmospheric Particulates Using TM Images in Beijing. Environ. Pollut. 2013, 178, 102–114. [Google Scholar] [CrossRef]
  38. Li, H.; Meier, F.; Lee, X.; Chakraborty, T.; Liu, J.; Schaap, M.; Sodoudi, S. Interaction between Urban Heat Island and Urban Pollution Island during Summer in Berlin. Sci. Total Environ. 2018, 636, 818–828. [Google Scholar] [CrossRef]
  39. Zheng, Z.; Ren, G.; Wang, H.; Dou, J.; Gao, Z.; Duan, C.; Li, Y.; Ngarukiyimana, J.P.; Zhao, C.; Cao, C.; et al. Relationship Between Fine-Particle Pollution and the Urban Heat Island in Beijing, China: Observational Evidence. Bound.-Layer Meteorol. 2018, 169, 93–113. [Google Scholar] [CrossRef]
  40. Sabrin, S.; Karimi, M.; Nazari, R. Developing Vulnerability Index to Quantify Urban Heat Islands Effects Coupled with Air Pollution: A Case Study of Camden, NJ. ISPRS Int. J. Geo-Inf. 2020, 9, 349. [Google Scholar] [CrossRef]
  41. Ulpiani, G. On the Linkage between Urban Heat Island and Urban Pollution Island: Three-Decade Literature Review towards a Conceptual Framework. Sci. Total Environ. 2021, 751, 141727. [Google Scholar] [CrossRef] [PubMed]
  42. Lindsey, R.; Dahlman, L. Climate Change: Global Temperature. Available online: http://www.climate.gov/news-features/understanding-climate/climate-change-global-temperature (accessed on 23 August 2023).
  43. Han, B.-S.; Park, K.; Kwak, K.-H.; Park, S.-B.; Jin, H.-G.; Moon, S.; Kim, J.-W.; Baik, J.-J. Air Quality Change in Seoul, South Korea under COVID-19 Social Distancing: Focusing on PM2.5. Int. J. Environ. Res. Public Health 2020, 17, 6208. [Google Scholar] [CrossRef] [PubMed]
  44. Kim, S.-U.; Kim, K.-Y. Physical and Chemical Mechanisms of the Daily-to-Seasonal Variation of PM10 in Korea. Sci. Total Environ. 2020, 712, 136429. [Google Scholar] [CrossRef] [PubMed]
  45. Oke, T.R. Boundary Layer Climates; Routledge: New York, NY, USA, 1987; ISBN 978-0-415-04319-9. [Google Scholar]
  46. Oke, T.R. The Distinction between Canopy and Boundary-layer Urban Heat Islands. Atmosphere 1976, 14, 268–277. [Google Scholar] [CrossRef]
  47. Crutzen, P.J. New Directions: The Growing Urban Heat and Pollution “Island” Effect—Impact on Chemistry and Climate. Atmos. Environ. 2004, 38, 3539–3540. [Google Scholar] [CrossRef]
  48. Zhou, B.; Rybski, D.; Kropp, J.P. On the Statistics of Urban Heat Island Intensity. Geophys. Res. Lett. 2013, 40, 5486–5491. [Google Scholar] [CrossRef]
  49. Li, Y.; Schubert, S.; Kropp, J.P.; Rybski, D. On the Influence of Density and Morphology on the Urban Heat Island Intensity. Nat. Commun. 2020, 11, 2647. [Google Scholar] [CrossRef]
  50. Dewan, A.; Kiselev, G.; Botje, D.; Mahmud, G.I.; Bhuian, H.; Hassan, Q.K. Surface Urban Heat Island Intensity in Five Major Cities of Bangladesh: Patterns, Drivers and Trends. Sustain. Cities Soc. 2021, 71, 102926. [Google Scholar] [CrossRef]
  51. Kim, Y.-H.; Baik, J.-J. Maximum Urban Heat Island Intensity in Seoul. J. Appl. Meteorol. Climatol. 2002, 41, 651–659. [Google Scholar] [CrossRef]
  52. Park, J.; Kim, J.-H.; Sohn, W.; Lee, D.-K. Urban Cooling Factors: Do Small Greenspaces Outperform Building Shade in Mitigating Urban Heat Island Intensity? Urban For. Urban Green. 2021, 64, 127256. [Google Scholar] [CrossRef]
  53. Jeong, S.-J.; Park, H.; Ho, C.-H.; Kim, J. Impact of Urbanization on Spring and Autumn Phenology of Deciduous Trees in the Seoul Capital Area, South Korea. Int. J. Biometeorol. 2019, 63, 627–637. [Google Scholar] [CrossRef] [PubMed]
  54. Sims, C.A. Macroeconomics and Reality. Econometrica 1980, 48, 1–48. [Google Scholar] [CrossRef]
  55. Wang, N.; Guo, J.; Liu, X.; Fang, T. A Service Demand Forecasting Model for One-Way Electric Car-Sharing Systems Combining Long Short-Term Memory Networks with Granger Causality Test. J. Clean. Prod. 2020, 244, 118812. [Google Scholar] [CrossRef]
  56. Zhao, K.; Chen, D.; Zhang, X.; Zhang, X. How Do Urban Land Expansion, Land Finance, and Economic Growth Interact? Int. J. Environ. Res. Public Health 2022, 19, 5039. [Google Scholar] [CrossRef] [PubMed]
  57. Sunde, T. Energy Consumption and Economic Growth Modelling in SADC Countries: An Application of the VAR Granger Causality Analysis. Int. J. Energy Technol. Policy 2020, 16, 41–56. [Google Scholar] [CrossRef]
  58. Dickey, D.A.; Fuller, W.A. Distribution of the Estimators for Autoregressive Time Series with a Unit Root. J. Am. Stat. Assoc. 1979, 74, 427–431. [Google Scholar] [CrossRef]
  59. Granger, C.W.J. Investigating Causal Relations by Econometric Models and Cross-Spectral Methods. Econometrica 1969, 37, 424–438. [Google Scholar] [CrossRef]
  60. Liu, H.; Kim, H. Ecological Footprint, Foreign Direct Investment, and Gross Domestic Production: Evidence of Belt & Road Initiative Countries. Sustainability 2018, 10, 3527. [Google Scholar] [CrossRef]
  61. Pereira, A.M.; de Frutos, R.F. Public Capital Accumulation and Private Sector Performance. J. Urban Econ. 1999, 46, 300–322. [Google Scholar] [CrossRef]
  62. Chaiechi, T.; Nguyen, T.M.T. Measuring Urban Economic Resilience of Two Tropical Cities, Using Impulse Response Analysis. Bull. Appl. Econ. 2021, 8, 59–79. [Google Scholar] [CrossRef]
  63. Bildirici, M.; Ozaksoy, F. The Relationship between Woody Biomass Consumption and Economic Growth: Nonlinear ARDL and Causality. J. For. Econ. 2017, 27, 60–69. [Google Scholar] [CrossRef]
  64. Liao, W.-C.; Zhao, D.; Lim, L.P.; Wong, G.K.M. Foreign Liquidity to Real Estate Market: Ripple Effect and Housing Price Dynamics. Urban Stud. 2015, 52, 138–158. [Google Scholar] [CrossRef]
Figure 1. Annual average temperatures of Seoul from 1961 to 2021 (data source: Open MET Data Portal, Korea Meteorological Administration, https://data.kma.go.kr/cmmn/main.do, accessed on 29 September 2023).
Figure 1. Annual average temperatures of Seoul from 1961 to 2021 (data source: Open MET Data Portal, Korea Meteorological Administration, https://data.kma.go.kr/cmmn/main.do, accessed on 29 September 2023).
Sustainability 15 16140 g001
Figure 2. Monthly average PM10 concentration of Seoul from January 2015 to December 2021 (data source: Seoul Air Quality Information, https://cleanair.seoul.go.kr/, accessed on 29 September 2023).
Figure 2. Monthly average PM10 concentration of Seoul from January 2015 to December 2021 (data source: Seoul Air Quality Information, https://cleanair.seoul.go.kr/, accessed on 29 September 2023).
Sustainability 15 16140 g002
Figure 3. Location of Seoul (left) and its 25 districts (right) (base map source: Google Maps (https://www.google.com/maps, accessed on 29 September 2023)).
Figure 3. Location of Seoul (left) and its 25 districts (right) (base map source: Google Maps (https://www.google.com/maps, accessed on 29 September 2023)).
Sustainability 15 16140 g003
Figure 4. Locations of weather and air quality monitoring stations in Seoul and their candidate reference stations outside Seoul to calculate heat and pollution island intensities.
Figure 4. Locations of weather and air quality monitoring stations in Seoul and their candidate reference stations outside Seoul to calculate heat and pollution island intensities.
Sustainability 15 16140 g004
Figure 5. Six pairs of variables for Granger causality tests (variable names are presented in lowercase italic grey letters).
Figure 5. Six pairs of variables for Granger causality tests (variable names are presented in lowercase italic grey letters).
Sustainability 15 16140 g005
Figure 6. Granger-causal relationships for the 6 pairs of variables.
Figure 6. Granger-causal relationships for the 6 pairs of variables.
Sustainability 15 16140 g006
Table 1. Descriptive statistics.
Table 1. Descriptive statistics.
VariableUnitObs.MeanStd. Dev
NameDescription
tempDaily average temperature°C45,65013.510.5
uhiiUHII°C45,6500.91.4
pm10Daily average PM10 concentrationμg/m345,65045.427.7
upiiUPIIμg/m345,6501.811.8
Table 2. Unit root test results for original data and seasonally differenced data.
Table 2. Unit root test results for original data and seasonally differenced data.
DistrictOriginal DataSeasonally Differenced Data
tempuhiipm10upiitempuhiipm10upii
Dobong−2.413−11.627 ***−9.368 ***−11.627 ***−22.302 ***−20.416 ***−20.878 ***−21.293 ***
Dongdaemun−2.438−7.619 ***−14.201 ***−10.796 ***−22.272 ***−20.919 ***−20.966 ***−21.898 ***
Dongjak−2.337−8.463 ***−14.457 ***−10.289 ***−22.541 ***−20.815 ***−20.751 ***−21.619 ***
Eunpyeong−2.428−8.516 ***−14.090 ***−10.172 ***−22.584 ***−21.196 ***−20.400 ***−20.974 ***
Gangbuk−2.417−12.662 ***−14.403 ***−9.061 ***−22.407 ***−20.949 ***−20.710 ***−21.487 ***
Gangdong−2.370−8.046 ***−14.009 ***−9.761 ***−22.113 ***−20.370 ***−21.114 ***−22.228 ***
Gangnam−2.411−12.029 ***−14.508 ***−23.839 ***−22.470 ***−21.515 ***−20.846 ***−22.123 ***
Gangseo−2.219−9.118 ***−9.158 ***−9.643 ***−20.952 ***−20.973 ***−20.973 ***−20.623 ***
Geumcheon−2.476−7.666 ***−14.402 ***−10.253 ***−22.559 ***−21.356 ***−20.543 ***−21.628 ***
Guro−2.537−8.237 ***−9.180 ***−9.967 ***−22.511 ***−20.245 ***−20.592 ***−21.963 ***
Gwanak−2.448−8.413 ***−14.411 ***−8.704 ***−22.529 ***−21.348 ***−20.925 ***−21.952 ***
Gwangjin−2.372−10.675 ***−8.895 ***−9.328 ***−22.232 ***−20.675 ***−20.744 ***−22.044 ***
Jongro−2.427−9.588 ***−8.778 ***−9.547 ***−22.548 ***−21.437 ***−21.080 ***−22.444 ***
Jung−2.359−16.546 ***−14.963 ***−9.473 ***−20.883 ***−22.264 ***−20.986 ***−21.829 ***
Jungrang−2.413−15.223 ***−8.532 ***−8.867 ***−22.309 ***−21.27 ***−23.393 ***−21.217 ***
Mapo−2.464−15.688 ***−9.222 ***−9.873 ***−22.965 ***−21.318 ***−20.889 ***−21.553 ***
Nowon−2.301−24.224 ***−9.050 ***−24.444 ***−22.065 ***−20.589 ***−21.022 ***−21.574 ***
Seocho−2.391−8.709 ***−8.289 ***−10.111 ***−22.589 ***−21.375 ***−20.641 ***−21.977 ***
Seodaemun−2.493−9.008 ***−8.97 ***−10.494 ***−22.802 ***−21.484 ***−20.734 ***−21.074 ***
Seongbuk−2.408−8.203 ***−9.341 ***−9.939 ***−22.437 ***−21.13 ***−23.455 ***−21.763 ***
Seongdong−2.419−6.425 ***−8.231 ***−9.873 ***−22.47 ***−21.224 ***−20.685 ***−21.476 ***
Songpa−2.450−8.075 ***−9.627 ***−9.641 ***−22.439 ***−20.473 ***−20.473 ***−21.281 ***
Yangcheon−2.400−8.516 ***−14.090 ***−10.172 ***−22.584 ***−21.196 ***−20.400 ***−20.974 ***
Yeongdeungpo−2.428−8.549 ***−9.425 ***−23.188 ***−22.618 ***−21.386 ***−20.934 ***−21.514 ***
Yongsan−2.401−8.837 ***−9.428 ***−9.454 ***−22.547 ***−21.275 ***−20.705 ***−21.684 ***
*** p < 0.001.
Table 3. Selected lag lengths for each district.
Table 3. Selected lag lengths for each district.
DistrictSelected Lag Length (Days)
Dobong5
Dongdaemun6
Dongjak6
Eunpyeong6
Gangbuk5
Gangdong6
Gangnam6
Gangseo6
Geumcheon6
Guro6
Gwanak6
Gwangjin6
Jongro5
Jung6
Jungrang5
Mapo6
Nowon6
Seocho6
Seodaemun6
Seongbuk6
Seongdong6
Songpa5
Yangcheon6
Yeongdeungpo6
Yongsan6
Table 4. Types of Granger-causal relationships and their corresponding districts.
Table 4. Types of Granger-causal relationships and their corresponding districts.
TypeDistricts
Type 1Sustainability 15 16140 i001Geumcheon, Mapo,
Yangcheon, Yeongdeungpo
Type 2Sustainability 15 16140 i002Gangdong, Jungrang,
Songpa, Yongsan
Type 3Sustainability 15 16140 i003Dongdaemun, Gangbuk,
Jung, Seongbuk
Type 4Sustainability 15 16140 i004Gangnam, Gwangjun,
Seongdong
Type 5Sustainability 15 16140 i005Dobong, Gangseo
Type 6Sustainability 15 16140 i006Seocho, Seodaemun
Type 7Sustainability 15 16140 i007Dongjak
Type 8Sustainability 15 16140 i008Eunpyeong
Type 9Sustainability 15 16140 i009Guro
Type 10Sustainability 15 16140 i010Gwanak
Type 11Sustainability 15 16140 i011Jongro
Type 12Sustainability 15 16140 i012Nowon
Table 5. Exemplary IRFs for twelve pairs of impulse and response variables.
Table 5. Exemplary IRFs for twelve pairs of impulse and response variables.
Impulse VariableResponse VariableExemplary IRFs
tempuhiiSustainability 15 16140 i013
Gangdong
pm10Sustainability 15 16140 i014
Gwanak
Sustainability 15 16140 i015
Songpa
upiiSustainability 15 16140 i016
Gangnam
Sustainability 15 16140 i017
Guro
Sustainability 15 16140 i018
Geumcheon
uhiitempSustainability 15 16140 i019
Gangdong
Sustainability 15 16140 i020
Nowon
pm10Sustainability 15 16140 i021
Gangdong
Sustainability 15 16140 i022
Seodaemun
upiiSustainability 15 16140 i023
Gangdong
Sustainability 15 16140 i024
Seodaemun
Sustainability 15 16140 i025
Nowon
pm10tempSustainability 15 16140 i026
Gwanak
Sustainability 15 16140 i027
Nowon
upiiSustainability 15 16140 i028
Gangdong
upiipm10Sustainability 15 16140 i029
Gangdong
tempSustainability 15 16140 i030
Seodaemun
Sustainability 15 16140 i031
Nowon
uhiiSustainability 15 16140 i032
Gangbuk
Sustainability 15 16140 i033
Jongro
Note: no IRFs are generated in cases where the impulse variable is pm10 and the response variable is uhii due to the absence of identified Granger causalities.
Table 6. Summary of FEVD results.
Table 6. Summary of FEVD results.
Impulse VariableResponse Variable
tempuhiipm10upii
temp93.79%~97.73%14.08%~40.04%3.66%~8.16%1.02%~4.89%
uhii0.68%~5.24%58.96%~85.12%0.35%~1.46%0.78%~2.48%
pm100.35%~1.11%0.23%~1.16%86.05%~92.68%5.90%~26.17%
upii0.30%~1.15%0.07%~1.11%1.84%~5.16%67.45%~87.95%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Youn, J.; Kim, H.; Lee, J. Relationships between Thermal Environment and Air Pollution of Seoul’s 25 Districts Using Vector Autoregressive Granger Causality. Sustainability 2023, 15, 16140. https://doi.org/10.3390/su152316140

AMA Style

Youn J, Kim H, Lee J. Relationships between Thermal Environment and Air Pollution of Seoul’s 25 Districts Using Vector Autoregressive Granger Causality. Sustainability. 2023; 15(23):16140. https://doi.org/10.3390/su152316140

Chicago/Turabian Style

Youn, Jeemin, Hyungkyoo Kim, and Jaekyung Lee. 2023. "Relationships between Thermal Environment and Air Pollution of Seoul’s 25 Districts Using Vector Autoregressive Granger Causality" Sustainability 15, no. 23: 16140. https://doi.org/10.3390/su152316140

APA Style

Youn, J., Kim, H., & Lee, J. (2023). Relationships between Thermal Environment and Air Pollution of Seoul’s 25 Districts Using Vector Autoregressive Granger Causality. Sustainability, 15(23), 16140. https://doi.org/10.3390/su152316140

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop