Next Article in Journal
Changes in Screening Test Volume in the National Breast and Cervical Cancer Early Detection Program during the COVID-19 Pandemic, 2020–2022
Next Article in Special Issue
Association between Heat Vulnerability Index and Stroke Severity
Previous Article in Journal
Taxonomy, Bio-Ecology and Insecticide Resistance of Anopheline Vectors of Malaria in Sri Lanka
Previous Article in Special Issue
Hourly Associations between Heat Index and Heat-Related Emergency Medical Service (EMS) Calls in Austin-Travis County, Texas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elucidating Uncertainty in Heat Vulnerability Mapping: Perspectives on Impact Variables and Modeling Approaches

1
Chungnam Economy Promotion Agency, Support Center for Community Revitalization, Yesan-gun 32416, Republic of Korea
2
Department of Geography, Kongju National University, Kongju-si 32588, Republic of Korea
3
Korea Adaptation Center for Climate Change, Korea Environment Institute, Sejong-si 30116, Republic of Korea
4
Research Center for Atmospheric Environment, Global Campus, Hankuk University of Foreign Studies, Yongin-si 17035, Republic of Korea
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2024, 21(7), 815; https://doi.org/10.3390/ijerph21070815
Submission received: 15 April 2024 / Revised: 27 May 2024 / Accepted: 28 May 2024 / Published: 21 June 2024

Abstract

:
Heat vulnerability maps are vital for identifying at-risk areas and guiding interventions, yet their relationship with health outcomes is underexplored. This study investigates the uncertainty in heat vulnerability maps generated using health outcomes and various statistical models. We constructed vulnerability maps for 167 municipalities in Korea, focusing on the mild and severe health impacts of heat waves on morbidity and mortality. The outcomes included incidence rates of heat-related outpatient visits (morbidity) and attributable mortality rates (mortality) among individuals aged 65 years and older. To construct these maps, we utilized 11 socioeconomic variables related to population, climate, and economic factors. Both linear and nonlinear statistical models were employed to assign these socioeconomic variables to heat vulnerability. We observed variations in the crucial socioeconomic variables affecting morbidity and mortality in the vulnerability maps. Notably, nonlinear models depicted the spatial patterns of health outcomes more accurately than linear models, considering the relationship between health outcomes and socioeconomic variables. Our findings emphasize the differences in the spatial distribution of heat vulnerability based on health outcomes and the choice of statistical models. These insights underscore the importance of selecting appropriate models to enhance the reliability of heat vulnerability maps and their relevance for policy-making.

1. Introduction

Climate change has induced a surge in extreme weather events worldwide, including heatwaves, cold waves, typhoons, floods, and droughts [1]. Heatwaves have become increasingly frequent in Europe, Australia, and South Korea, primarily due to the compounding effects of global warming and urbanization [2,3]. Numerous regions of South Korea have experienced unprecedented temperature spikes since 2018, marking a significant departure from historical meteorological records [4].
The adverse health effects of heat waves have been extensively documented in various countries. For instance, in 2003, approximately 72,210 deaths were attributed to heatwaves across 15 European nations [5]. Similarly, the UK reported 892 excess deaths during a heatwave in 2019, while [6], Portugal and Spain experienced 1,063 and 679 excess deaths, respectively, during heatwaves in June and July 2022 [7]. South Korea also faced a grim toll, with an estimated 800–900 excess deaths attributed to heatwaves in 2018 [8,9].
Future projections indicate that climate change and urbanization will exacerbate the intensity, frequency, and duration of heatwaves [1]. Given the profound implications of heatwaves on human health, understanding the regions susceptible to heatwaves is of paramount importance for public health and safety. Effective response strategies and protective measures hinge on the accurate identification of vulnerable areas.
Heat vulnerability maps serve as invaluable tools for identifying regions vulnerable to heat waves. By visually highlighting areas with heightened vulnerability, these maps can assist policymakers in devising targeted intervention measures. Additionally, they raise awareness of health-related issues and aid public health organizations in formulating appropriate response strategies. Moreover, heat vulnerability maps can be used to anticipate spatial pattern changes in vulnerability, rendering them indispensable for urban planning and sustainable development.
Previous studies on heat vulnerability maps have predominantly focused on analyzing regional vulnerability using socioeconomic variables. Various methodologies, such as principal component (PC) analysis and equal weighting assumptions, have been employed for variable selection and weighting [10,11,12,13,14,15]. Research on the uncertainty of heat vulnerability maps is of great interest to researchers. Some studies have shown that caution is required when relying solely on heat vulnerability maps constructed using a single method compared to those constructed using different methods [16,17].
However, there remains a dearth of studies exploring the uncertainty surrounding the explanatory power of health outcomes in heat vulnerability maps. Health outcomes characterize heat-related outpatient visits as a mild impact, and emergency room admissions and mortality as a severe impact [18]. Several studies have compared the explanatory power of health outcomes and heat vulnerability maps but only used mortality as an indicator [17,19]. Most other studies have not attempted to compare health outcomes and heat vulnerability maps. Therefore, little is known about the relationship between heat vulnerability maps and health outcomes. For instance, it is difficult to determine which heat impact indicator’s spatial distribution is better explained by heat vulnerability maps and whether the same variables can be used to construct vulnerability maps for mild and severe health outcomes. Understanding the relationship between health outcomes and heat vulnerability maps is essential for identifying vulnerability factors and developing policies. Differences in socioeconomic vulnerability factors, depending on the impact indicator, could provide critical information for policy development.
Previous studies have assumed a linear relationship between independent variables and health outcomes and constructed heat vulnerability maps through a combination of independent variables. However, the possibility of a nonlinear expression of the influence of specific variables due to interactions between independent variables and the existence of thresholds in the influence of independent variables on dependent variables has been overlooked.
Additionally, the impact of heatwaves varies depending on socioeconomic conditions such as age, occupation, and income, with older adults being well known as a group vulnerable to heatwaves [20,21,22]. Limited mobility and chronic illnesses (e.g., cardiovascular diseases) of the elderly may affect their ability to adequately protect themselves and respond promptly to extreme weather [23]. Furthermore, their acquisition and utilization of information related to heatwaves are limited [24,25].
The increasing number of elderly people vulnerable to heatwaves could exacerbate the societal and economic risks of heatwaves, leading to significant social problems. South Korea is one of the countries experiencing rapid aging globally, with the proportion of the elderly population increasing from 7.2% in 2000 to 16.6% in 2021 [26]. Sim et al. calculated the exposure risk of the population to heatwaves based on future scenarios of heatwave days and national population estimates, forecasting that rapid aging would be a significant factor in increasing the exposure risk to heatwaves [27]. An increase in the elderly population in the future is expected to escalate the scale of heatwave damage alongside an increase in extreme weather phenomena [28]. Therefore, despite several studies, continued interest and research on the elderly are needed.
The HVI (Heat Vulnerability Index) mapping can be used as a proxy measure for potential health impacts when comprehensive health data is not available. Furthermore, it provides various insights that are difficult to identify from simple health impact indicators, offering advantages in policy development and future risk prediction. HVI mapping takes a multifaceted approach by integrating demographic, socioeconomic, and climate variables instead of relying on a single indicator or variable. This allows for a more comprehensive assessment of heat vulnerability in specific regions, enabling the analysis of key factors contributing to regional vulnerability. This approach offers a more holistic understanding compared to analyzing individual health outcomes alone. It provides policymakers with the information necessary for developing tailored policies and effectively allocating resources to address socioeconomic weaknesses in specific regions. Additionally, while direct health impact indicators show the current state or past impacts, HVI can help predict potential future risks. This plays a crucial role in minimizing health impacts through preventive measures. Therefore, the approach to improving HVI mapping proposed in this study will play a vital role in predicting future risks, developing tailored policies, and enhancing preventive measures.
The objective of this study was to enhance the understanding of the uncertainties that need to be considered when constructing heat vulnerability maps. First, we aimed to examine the relationship between heat vulnerability outcomes and health outcomes, proposing a method for combining heat vulnerability outcomes based on linear and nonlinear relationships. Furthermore, this study aimed to construct heat vulnerability maps based on the nonlinear relationship between heat vulnerability and health outcomes. Finally, we aimed to investigate differences in socioeconomic vulnerability factors associated with various health outcomes of heatwaves, thereby providing critical insights for policy development.

2. Materials and Methods

2.1. Data

The spatial scope of this study encompassed 167 municipalities in South Korea (Figure 1). The analysis covered the period from 2007 to 2020, spanning 14 years. Two negative health outcomes of for heatwaves were considered: morbidity, representing a mild impact, and mortality, representing a severe impact. Morbidity was measured using the incidence rate of heat-related illnesses during the summer months (June–August), whereas mortality was assessed using the proportion of excess deaths attributable to heatwaves.
Data on outpatient visits for heat-related illnesses were obtained from the National Health Insurance Service of South Korea, utilizing healthcare big data. Heat-related illnesses (HRI) were defined using the T67 code from the 10th Revision of the International Classification of Diseases (ICD-10). Mortality data were acquired from Statistics Korea’s cause-of-death statistics, including all deaths from external and other causes, and extracted using ICD-10 codes A-R. In this study, older adults were defined as individuals aged 65 and above. The number of outpatient visits for heat-related illnesses among older adults during the summer months was 70,685, and the total number of older adult deaths over the entire period was 4,006,289.
A total of 11 socioeconomic variables were used to create the heat vulnerability map. A list of variables and a map for each variable can be seen in Table 1 and Figure 1, respectively. Population variables include factors related to the elderly population, occupation, and household structure. The variables included the proportion of the population aged 65 or older (pop_65), the proportion of the population in rural areas (agri_65), and the proportion of the population living alone (single_65).
Economic variables include factors that reflect the income level of the elderly population and social infrastructure. The income level of the elderly population (income_65) was calculated as the ratio of low-income households among the population aged 65 or older. For low-income households, data from the lower insurance premium class (bottom 20%) were used. Financial independence refers to the ratio of self-revenue, which is the basis of a local government’s independent financial capacity, to its budget size. It is an indicator showing the degree of self-reliance based on each local government’s internal revenue, and the higher the ratio, the better the self-reliant basis of the local government.
For climate variables, temperature and humidity data from each region from 2007 to 2020 (14 years) were used to explain exposure to heat waves. The 90th percentile temperature (t90) of each region represents the level of high temperature throughout the year. The number of days exceeding the 95th percentile (tmax95) indicates the number of days when the daily maximum temperature in summer was above the 95th percentile, reflecting the intensity and frequency of heat waves. The daily average temperature (tag) in summer for each region reflects the average summer temperature.
Summers in Korea tend to have high humidity, which is closely related to the occurrence of heat-related illnesses. Therefore, we included humidity data in our analysis. The number of days (rhav80) with relative humidity above the 80th percentile in summer for each region indicates the intensity and frequency of high humidity. Additionally, the average summer relative humidity (rhav_avg) of each region was included as a variable.
We adopted the daily maximum temperature for our study. Although some existing studies have used average temperature [29], those analyses included the influence of winter. In contrast, this study focuses solely on the summer period. Several studies in Korea have demonstrated the health effects of the daily maximum temperature, and the heatwave warning system has long been based on this metric [18,22]. Therefore, our analysis was performed using the daily maximum temperature.
Furthermore, the ratio of forest to total area was used as a proxy variable for assessing the impact of green space on alleviating the thermal environment.
Temperature and humidity data, consisting of climate variables, were obtained from the Automated Synoptic Observation System of the Korea Meteorological Administration. Area average data for each region were assigned using the Thiessen polygon method. Commonly used interpolation methods for meteorological applications include nearest-station assignment, inverse-distance weighting, inverse-distance-square weighting, Thiessen polygon method, orthogonal-polynomial approximation, Lagrange method, interpolation by splines, kriging, and interpolation by empirical orthogonal functions. Each method has its merits and is applicable based on the temporal scale, spatial scale, stationarity, and variability of the field under consideration [30].
The Thiessen polygon approach determines the average weather by considering only the representative area of a station, which is determined by the bisector between each pair of stations. Although this method has limitations, such as not accounting for topographic effects, it is useful for producing long-term regional average weather data [31].

2.2. Ethical Approval

The study protocol underwent thorough review and approval by the Institutional Review Board (IRB) of Kongju National University (Confirmation No. 2022-46). Given the observational nature of this study and the use of anonymized statistical data, the IRB committee waived the requirement for informed consent. All methodologies adhered to the guidelines stipulated by the Korean government for health and medical data usage.

2.3. Methodology for Impact Indicator Calculation

The attributable mortality rate of excess deaths was calculated based on the statistical relationship between the daily maximum temperature and mortality. Initially, the relationship between the daily maximum temperature and mortality was analyzed regionally. The Distributed Lag Non-Linear model (DLNM) was adopted for this purpose. The DLNM, known for its ability to account for the lag effects of independent variables, has been widely utilized in recent epidemiological research [22]. The model was formulated as shown in Equation (1), and analysis was conducted using R version 4.3.0 with the dlnm package.
I n y i j = α + β T m a x i j , j l + N S s n i j + N S r h a v g i j + N S d o y i j + γ w e e k d a y i j + ε
where Yij represents the number of deaths recorded in region j on day i. Tmax denotes the daily maximum temperature, and l represents the lag effect. Considering the tendency for the impact of heatwaves to be strongest on the same day and rapidly decrease within a short period [32,33,34], we set the lag period to three days. Here, Rhavg refers to the average daily humidity, doy represents the day of the year in terms of Julian days, and sn denotes the serial number. The degree of freedom for this variable was set to 6 and multiplied by 14 years (the duration of the study data). Weekdays indicate the days of the week. Natural cubic splines (NS) were used to examine the non-linear relationship between the dependent and independent variables. Based on the relative risk of death calculated from temperature using the aforementioned method, excess deaths for each region were calculated using Equation (2) [8].
Y i j = D i j × e β × Δ T i j 1 e β × Δ T i j × H i j
where i and j represent the variables for the date and region, respectively, Yij denotes the number of excess deaths that occurred in region j on day i. D represents the daily deaths, and β indicates the slope of the relative risk above the threshold temperature. The threshold temperature for region j was determined by observing the non-linear relationship between the temperature and deaths mentioned above. ΔT represents the difference between the threshold and the temperature for each day. If the temperature in region j on day i exceeded the threshold temperature, H was 1; otherwise, it was 0.
The contribution rate of excess deaths for each region was calculated by averaging the percentage of excess deaths compared to the total number of deaths per year. The heat-related illness incidence rate was calculated by averaging the number of heat-related illness cases per ten thousand elderly population per year.
Throughout all stages of the analysis, including the calculation of excess death contribution rates, heat-related illness incidence rates, and socioeconomic variables, data for seven metro cities (Seoul, Busan, Daegu, Incheon, Gwangju, Daejeon, and Ulsan) were exclusively used for each city. However, for other cities and rural areas, data were integrated from adjacent regions because the sample size for these areas was sometimes inadequate due to low variability and limited data availability to establish statistical models. Therefore, to ensure an adequate sample size for the analysis, data from adjacent regions were integrated. Consequently, the interpretation of the results considered these areas from a regional perspective, rather than focusing on individual cities outside metro areas.

2.4. Methodology for Analysis of Relationships between Variables

The analysis of the relationships between health outcomes and socioeconomic variables proceeded as follows: Firstly, the socioeconomic variables were transformed into components without multicollinearity issues using a PC analysis. Subsequently, these components were utilized as independent variables, while the relative risks of morbidity and mortality were employed as dependent variables to analyze their relationships. Both linear and nonlinear models, represented by Equations (3) and (4) respectively, were used in this analysis.
Y = β0 + β1X1 + β2X2 + … + βnXn + ϵ
Y = β0 + f1(X1) + f1(X2) + … + fn(Xn) + ϵ
where Y represents the dependent variable, which can be either the incidence rate of heat-related illnesses or the attributable mortality rate. β0 denotes the intercept of the model, while X1–Xn represents the independent variables. β1–βn are the coefficients for each independent variable. f1–fn represents the smoothing functions for each independent variable to account for nonlinearity. For linear and nonlinear models, we used R’s linear model (LM) and the generalized additive model (GAM), respectively.

2.5. Methodology for Generating Heat Vulnerability Maps

Heat vulnerability maps were constructed using two approaches. The first was based on the results of a linear model. Utilizing the linear relationship analysis mentioned earlier, vulnerability maps were created using PCs that exhibited statistically significant relationships (p < 0.05) with health outcomes (incidence rate of heat-related illnesses and attributable mortality rate) at a 95% confidence level. The PC values were assigned to nine risk levels, with the direction of risk allocation determined based on the relationship (negative or positive) between the PCs and the health outcomes. The final heat vulnerability map was constructed by averaging the maps assigned to the PCs. In this study, we refer to this method as the Heat Vulnerability Index based on the Linear Model (HVI-LM).
The second approach was based on the results of a nonlinear model (NLM). This map was created based on previously analyzed nonlinear relationships. Firstly, we identified the thresholds in the nonlinear relationship between the health outcomes and PCs. The changes observed above these thresholds were used to assign the values of the PCs to the nine risk levels. We refer to this method as the Heat Vulnerability Index based on the Nonlinear Model (HVI-NLM).

3. Results

3.1. Health Outcoms of Heatwave

The incidence rate of heat-related illnesses and the attributable mortality rate from heat waves exhibited spatial variation. The incidence of heat-related illnesses was generally higher in the southwestern region of Korea and lower in the northeastern region (Figure 2a). In contrast, the areas with higher attributable mortality rates were dispersed across various regions (Figure 2b). Due to the complexity of the spatial distribution of mortality compared to that of diseases, it may be difficult to statistically explain the spatial distribution of mortality. This observation reinforces the findings of the analysis of the relationship between health outcomes and PCs.

3.2. Principal Components

PC analysis of the socioeconomic variables revealed four PCS, named PC1, PC2, PC3, and PC4. Table 2 summarizes the results of the PC analysis conducted using the 11 variables. Four PCs showed initial eigenvalues exceeding 1, and the cumulative explained variance for the first four components was 86.08%.
Table 3 lists the weights of the variables for the four PCs. Based on this table, we identified the key variables of each PC and defined the nature of each component. PC1 represented the demographic characteristics of the elderly population in the area and socioeconomic conditions. This component consisted of variables such as the proportion of the elderly population (pop_65), the proportion of the elderly engaged in agriculture (agri_65), the proportion of the elderly living alone (single_65), and the financial autonomy of the local government (financial). PC2 represented the climatic component of the area, comprising temperature (t90 and tmax95) and humidity variables (rhav80 and rhav_avg). PC3 consisted of the cover ratio of the forest (per_forest), and average temperature (tavg), defining a combination of temperature and a variable that can moderate its effects. Lastly, PC4 was composed of the low-income elderly population ratio (lowincome_65) and temperature variable, defining a component of the economic conditions of the elderly and the possibility of exposure to heatwaves.

3.3. Relationship between PCs and Impacts

We derived three key findings regarding the relationship between the PCs and health outcomes. First, the NLM showed higher statistical explanatory power for the spatial distribution of health outcomes than the linear model. Table 4 summarizes the p-values of the PCs in the analysis of their relationships with the health outcomes. The R2 values of the LM were 0.381 and 0.154 for morbidity and mortality, respectively. In contrast, the R2 values for the NLM were 0.518 and 0.486 for morbidity and mortality, respectively. The NLM performed better than the LM in terms of both morbidity and mortality. Secondly, the second PC explained morbidity better than mortality. For LM, the R2 value for morbidity was 0.381, which was higher than that for the mortality rate of 0.154. Similarly, in the NLM, the R2 for morbidity was 0.518, compared to 0.486 for mortality. In the NLM, the explained deviations were 55.7% for morbidity and 52.5% for mortality. Thirdly, statistically significant PCs varied between morbidity and mortality. PC1 was significantly associated with morbidity, whereas PC4 was significantly associated with mortality. PC3 was a statistically significant variable in all analyses.
In the LM model, PC2 had a negative relationship, while PC3 had a positive relationships with the health outcomes. The linear regression slopes of PC2 for morbidity and mortality were −1.1151 and −0.1113, respectively. On the other hand, the linear regression slopes for PC3 in terms of morbidity and mortality were 3.1630 and 0.1010, respectively. In the case of morbidity, the slope of PC1 was 2.0553, while in the case of mortality, the slope of PC4 was −0.1414.
Similar relationships between health outcomes and PCs were observed in the NLM (Figure 3). PC1 exhibited a nonlinear relationship with morbidity, showing an increasing trend in morbidity when PC1 was not significant in mortality. In morbidity, PC2 and PC3 had statistically significant negative and positive correlations, respectively. PC3 did not show a statistically significant relationship with morbidity. Mortality showed a nonlinear relationship with PC3 and PC4. PC1 and PC2 did not show a statistically significant relationship with mortality. PC3 exhibited a positive relationship with mortality. PC4 served as a threshold point at approximately −1, with no statistically significant change in mortality above −1. However, mortality increased distinctly below −1.

3.4. Comparison between HVI-LM and HVI-NLM

Based on the analysis results, the values of the PCs were assigned to a heat vulnerability index (HVI) in nine intervals ranging from −4 to 4 (Table 5). Each interval was considered using probabilities of 1%, 5%, 10%, 25%, 50%, 75%, 90%, 95%, and 99% after standardizing the PCs using means and standard deviations. When assigning the PCs based on the results of HVI-LM, if the linear regression coefficient was positive, the values of each PC were sequentially assigned from −4 to 4. Conversely, if the linear regression coefficient was negative, the values of each PC ranged from 4 to −4. For the HVI-NLM, the HVI was assigned only to the intervals of PCs where a relationship between the dependent and independent variables was observed.
For instance, in a nonlinear relationship, when PC1 exceeded approximately 1, the incidence rate of heat-related illnesses showed a dramatic increase. When PC1 was near 0, no statistically significant change was recorded in the incidence rate of heat-related illnesses. Therefore, the HVI gradually increased from the negative extreme value of PC1 to the positive values. However, the HVI was not assigned around 0 for PC1. Unlike linear relationships, assigning an HVI of 0 for intervals where statistically significant relationships were not observed in nonlinear relationships could result in excessive smoothing of the final results because the HVI was averaged for all PCs.
After assigning the HVI based on the range of each PC, a heat vulnerability map was created by averaging the HVIs of the PCs (Figure 4). The results showed significant differences in the spatial distribution of heat vulnerability predictions depending on the impact of heat waves (morbidity and mortality) and vulnerability assignment methods (linear and nonlinear). In the heat vulnerability maps for morbidity and mortality, areas with high vulnerability to morbidity were concentrated in the southwest region, whereas vulnerability to mortality was generally observed in coastal areas such as the western and eastern coastal regions.
Furthermore, even when the same PC was used, distinct differences were observed between (a) HVI-LM-morbidity and (c) HVI-NLM-morbidity, as well as between (b) HVI-LM-mortality and (d) HVI-NLM-mortality. In (a) HVI-LM-morbidity, vulnerable areas were concentrated in the southwest region, whereas in (c) HVI-NLM-morbidity, vulnerable areas included broader inland areas in the west and southeast compared with (a). In (b) HVI-LM-mortality, vulnerability in the northeast region was classified as low, whereas in (d) HVI-NLM-mortality, vulnerability was classified as high. Differences were also noted in the vulnerability between the two maps in the southeastern region.

4. Discussion

This study demonstrated significant differences in heat vulnerability maps based on health outcomes and statistical models. The results revealed notable variations in socioeconomic variables that were statistically significant for morbidity and mortality risks. For morbidity, PC1 (size of the elderly population and socioeconomic conditions) and PC3 (temperature and the effect of green space) showed statistically significant explanatory power. In contrast, PC2 (climate conditions) and PC4 (income of the elderly and temperature) demonstrated statistically significant explanatory power for mortality. These differences contributed to variations in the heat vulnerability maps for morbidity and mortality.
Furthermore, this study identified nonlinear relationships between socioeconomic variables and health outcomes. Several existing studies that have constructed heat vulnerability maps [13,14,15,16,17] have assumed a linear relationship between socioeconomic variables and vulnerability, assigning vulnerability based on changes in socioeconomic variables. However, our study indicates that NLMs better explain the relationship between socioeconomic variables and health outcomes than LMs. Heat vulnerability maps based on nonlinear relationships matched the spatial patterns of impacts better than those using linear models.
The implications of these findings underscore the importance of considering health outcomes and statistical modeling factors when constructing heat vulnerability maps. Given the utilization of heat vulnerability maps as a tool for formulating heat response policies, differences in health outcomes have significant implications. A heat vulnerability map that does not specifically identify heat-related health risks may have limitations when utilized in policymaking. Regions with high morbidity risks may require education and preventive measures to raise awareness about heat-related illnesses, whereas regions with high mortality risks may require strengthening of healthcare facilities, emergency services, and support for vulnerable populations, such as low-income individuals and the disabled. Therefore, ensuring clarity in health outcomes when developing heat vulnerability maps can enhance policy utility.
Previous studies have not addressed the various types of damages that can be identified using heat vulnerability maps. According to prior research, because the temperature thresholds for morbidity and mortality differ, different response policies are required in terms of timing and risk communication. Since morbidity tends to precede mortality, studies have indicated challenges in providing adequate warnings for a large number of morbidity cases when heat warning systems are primarily designed to address mortality [18,35]. Similarly, the distinct distribution of heat vulnerability for morbidity and mortality, along with differences in the PCs associated with mortality and morbidity, implies the need for tailored response policies for each region. Therefore, constructing separate heat vulnerability maps for morbidity and mortality could assist in formulating customized policies for each region.
However, this study has certain limitations that highlight the need for further investigation. First, additional research is required to enhance the explanatory power of the spatial distribution of heat-related damage. The explanatory power for morbidity and mortality in this study was 55.7% and 52.5%, respectively, indicating that there is room for improvement. This could involve exploring dynamic data that reflect societal activities. In an existing case targeting a single city, the coefficient of determination between the heat vulnerability index and heat-related disease deaths was 0.58 using an improved method and 0.32 using the existing method [16]. This suggests that it is challenging to accurately reproduce actual impacts with the heat vulnerability index. The selection of different variables can lead to significantly different results on heat vulnerability maps. Therefore, efforts to develop variables and find appropriate combinations are crucial [36].
For example, there is a need to develop variables for social ties [37], effective heat wave warning systems [24,38,39] and the quality and accessibility of medical services [40]. Additionally, variables representing the operational and activity levels of social and cultural facilities should be considered. Cultural facilities not only serve as shelters from heat but also provide gathering spaces and activities. In Seoul, for example, cultural facilities extend their operating hours during the summer and offer various programs. The operations of cultural facilities encompass structural (installation, development, and maintenance), nonstructural (investment, support, and monitoring), and social (collaboration, education, promotion, and guidelines) measures. Quantifying regional differences in such dynamic variables can be incorporated into the construction of heat vulnerability maps.
Additionally, in terms of modeling, it is necessary to apply the NLMs identified in this study while determining the extent to which variables should be minimized or maximized to enhance the model’s explanatory power. Finding the appropriate balance is crucial. Although the integration of factor analysis and DLNM used in this study was essential for minimizing variables and identifying statistically significant causal relationships, the DLNM model, which reduces the variance of model estimates in a nonlinear fashion, poses challenges in interpreting the results. Future research is required to enhance the overall reliability of DLNM and select significant factors that can differentiate between urban and non-urban areas.
Secondly, when using NLMs, it is necessary to establish a more systematic process for assigning heat vulnerabilities. The finding that NLMs yield better results than LMs is significant. However, in NLMs, the assignment of vulnerability values based on researcher observations introduces the potential for subjective judgment, which can hinder the construction of heat vulnerability maps based on NLMs. To address this issue, it is necessary to accumulate various research cases based on NLMs and establish a systematic assignment process.
Despite these challenges, this study identified the impact variables associated with heat-related morbidity and mortality and revealed the spatial patterns of heat vulnerability for morbidity and mortality at the municipal level. Addressing the limitations and conducting further research will contribute to the refinement and effectiveness of heat vulnerability mapping.

5. Conclusions

This study highlights the differences in heat vulnerability maps based on health outcomes and statistical models, underscoring the significant impact of socioeconomic variables on morbidity and mortality risks. Our findings indicate that the variables significantly associated with morbidity and mortality differed. Additionally, non-linear models offered a more accurate representation of the relationship between socioeconomic variables and health outcomes compared to linear models. These results emphasize the importance of selecting appropriate statistical models to enhance the accuracy of heat vulnerability maps and the need to configure variables differently depending on the specific health outcomes.
Several limitations were identified in the heat vulnerability mapping process. Firstly, there is the issue of low explanatory power. Addressing this challenge requires efforts to develop various variables beyond those proposed in this study. Variables such as social cohesion, effective heat warning systems, and the quality and accessibility of healthcare may be crucial. Furthermore, research is needed to minimize subjective judgment when using nonlinear models. To enhance model reliability, establishing a systematic process for assigning the results of nonlinear models to heat vulnerability is essential.
In conclusion, our study underscores the importance of considering both health outcomes and statistical modeling factors in constructing heat vulnerability maps. These maps serve as crucial tools to inform heat response policies and facilitate targeted interventions based on specific regional vulnerabilities. Insights from this study can assist national and local governments in prioritizing areas for heat response and management. Additionally, our findings offer valuable guidance for researchers developing heat vulnerability maps and advocate for the integration of advanced statistical modeling techniques with health outcomes.

Author Contributions

Conceptualization, S.J. and C.Y.; methodology, S.J.; software, Y.K. and Y.L; formal analysis, Y.K. and Y.L.; writing—original draft preparation, S.J.; writing—review and editing, C.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a grant (2020-MOIS35-002) of the Policy-linked Technology Development Program on Natural Disaster Prevention and Mitigation funded by Ministry of Interior and Safety (MOIS, Republic of Korea).

Institutional Review Board Statement

The study protocol was thoroughly reviewed and approved by the Institutional Review Board (IRB) of Kongju National University (Confirmation No. 2022-46).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from Korea’s National Health Insurance Service; however, restrictions apply to the availability of these data, which were used under license for the current study and are not publicly available. However, the final analysis data are available from the authors upon reasonable request and with permission from the National Health Insurance Service. Please contact the author (Yeon-yeop Lim, [email protected]) if someone wants to request data from this study and submit the revised manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. IPCC. Climate Change 2022: Impacts, Adaptation, and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Pörtner, H.-O., Roberts, D.C., Tignor, M., Poloczanska, E.S., Mintenbeck, K., Alegría, A., Craig, M., Langsdorf, S., Löschke, S., Möller, V., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2022; p. 3056. [Google Scholar]
  2. Climate Change Service. Record-Breaking Temperatures for June. 2019. Available online: https://climate.copernicus.eu/record-breaking-temperatures-june (accessed on 22 March 2024).
  3. Bureau of Meteorology, AUS. Climate Change—Trends and Extremes. 2020. Available online: https://www.bom.gov.au/climate/change (accessed on 1 June 2020).
  4. Korea Meteorological Administration. Korean Climate Change Assessment Report; Korea Meteorological Administration: Daejeon, Republic of Korea, 2020. [Google Scholar]
  5. UNDRR (UN Office for Disaster Risk Reduction). The Human Cost of Disaster: An Overview of the Last 20 Years (2000–2019); UNDRR: Geneva, Switzerland, 2020. [Google Scholar]
  6. PHE (Public Health England). PHE Heatwave Mortality Monitoring; PHE: London, UK, 2019; pp. 1–9.
  7. Korea Meteorological Administration. Abnormal Climate Report 2022; Korea Meteorological Administration: Daejeon, Republic of Korea, 2023. [Google Scholar]
  8. Lim, Y.-H.; Lee, H.; Hong, Y.-C. Attributable all-cause mortality during heat waves in South Korea, 2006–2018. Public Health Wkly. Rep. 2019, 12, 1435–1441. [Google Scholar]
  9. Park, J.; Chae, Y. Analysis of heat-related illness and excess mortality by heat waves in South Korea in 2018. J. Korean Geogr. Soc. 2020, 55, 391–408. [Google Scholar] [CrossRef]
  10. Bai, L.; Woodward, A.; Cirendunzhu, C.; Liu, Q. County-level heat vulnerability of urban and rural residents in Tibet, China. Environ. Health 2016, 15, 3. [Google Scholar] [CrossRef]
  11. Harlan, S.L.; Declet-Barreto, J.H.; Stefanov, W.L.; Petitti, D.B. Neighborhood effects on heat deaths: Social and environmental predictors of vulnerability in Maricopa County, Arizona. Environ. Health Perspect. 2013, 121, 197–204. [Google Scholar] [CrossRef]
  12. Inostroza, L.; Palme, M.; de la Barrera, F. Heat vulnerability index: Spatial patterns of exposure, sensitivity and adaptive capacity for Santiago, Chile. PLoS ONE 2016, 11, e0162464. [Google Scholar] [CrossRef]
  13. Harlan, S.L.; Brazel, A.J.; Prashad, L.; Stafanov, W.; Larsen, L. Neighborhood microclimates and vulnerability to heat stress. Soc. Sci. Med. 2006, 63, 2847–2863. [Google Scholar] [CrossRef]
  14. Aubrecht, C.; Özceylan, D. Identification of heat risk patterns in the U.S. National Capital Region by integrating heat stress and related vulnerabilities. Environ. Int. 2013, 56, 65–77. [Google Scholar] [CrossRef]
  15. Zhu, Q.; Liu, T.; Lin, H.; Xiao, J.; Luo, Y.; Zeng, W.; Zeng, S.; Wei, Y.; Chu, C.; Baum, S.; et al. Spatial distribution of health vulnerability to heat waves in Guangdong Province, China. Glob. Health Action 2014, 7, 25051. [Google Scholar] [CrossRef]
  16. Guo, X.; Huang, G.; Jia, P.; Wu, J. Estimating fine-scale heat vulnerability in Beijing through two approaches: Spatial patterns, similarities, and divergence. Remote Sens. 2019, 11, 2358. [Google Scholar] [CrossRef]
  17. Liu, X.; Yue, W.; Yang, X.; Hu, K.; Zhang, W.; Huang, M. Mapping urban heat vulnerability of extreme heat in Hangzhou via comparing two approaches. Complexity 2020, 2020, 9717658. [Google Scholar] [CrossRef]
  18. Chae, Y.; Park, J. Analysis of effectiveness of impact based heatwave warning considering severity and likelihood of health impact in Seoul, Korea. Int. J. Environ. Res. Public Health 2021, 18, 2380. [Google Scholar] [CrossRef]
  19. Kim, J.; Jung, T.Y. Assessing health sector climate vulnerability in 226 local entities of South Korea based on principal component analysis. Urban Clim. 2023, 49, 101521. [Google Scholar] [CrossRef]
  20. Yardley, J.; Sigal, R.J.; Kenny, G.P. Heat health planning: The importance of social and community factors. Glob. Environ. Chang. 2011, 21, 670–679. [Google Scholar] [CrossRef]
  21. Park, J.; Chae, Y. Analysis of time-dependent effects of ambient temperatures on health by vulnerable groups in Korea in 1999–2018. Sci. Rep. 2023, 13, 922. [Google Scholar] [CrossRef]
  22. Kang, Y.; Baek, I.; Park, J. Assessing heatwave effects on disabled persons in South Korea. Sci. Rep. 2024, 14, 3459. [Google Scholar] [CrossRef]
  23. Flynn, A.; McGreevy, C.; Mulkerrin, E.C. Why do older patients die in a heatwave? QJM 2005, 98, 227–229. [Google Scholar] [CrossRef]
  24. Takahashi, N.; Nakao, R.; Ueda, K.; Ono, M.; Kondo, M.; Honda, Y.; Hashizume, M. Community trial on heat-related illness prevention behaviors and knowledge for the elderly. Int. J. Environ. Res. Public Health 2015, 12, 3188–3214. [Google Scholar] [CrossRef]
  25. 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]
  26. Lee, Y. Is an aging society a social problem? Aging from a complex perspective. J. Korean Urban Geogr. Soc. 2023, 26, 1–16. [Google Scholar]
  27. Shim, C.; Kim, O.; Han, J.; Song, S.; Na, G.; Kim, K. A Projected Population Exposed to Future Hot Weather Events Considering the Korean Society with Low Birth Rate and Rapid Aging: For Supporting Regional Policy Measures; Korea Environment Institute: Sejong, Republic of Korea, 2019; pp. 1–192. [Google Scholar]
  28. Lim, Y.; Kim, O.S.; Park, J. Projection of mortality risk of elderly to high ambient temperature under climate change scenarios in Korea during the 2050s and the 2090s. Landsc. Geogr. 2023, 33, 51–64. [Google Scholar]
  29. Gasparrini, A.; Guo, Y.; Hashizume, M.; Lavigne, E.; Zanobetti, A.; Schwartz, J.; Tobias, A.; Tong, S.; Rocklöv, J.; Forsberg, B.; et al. Mortality risk attributable to high and low ambient temperature: A multicountry observational study. Lancet 2015, 386, 369–375. [Google Scholar] [CrossRef] [PubMed]
  30. Park, J.; Jang, D.-H. Application of MK-PRISM for interpolation of wind speed and comparison with co-kriging in South Korea. GIScience Remote Sens. 2016, 53, 421–443. [Google Scholar] [CrossRef]
  31. Shen, S.S.; Dzikowski, P.; Li, G.; Griffith, D. Interpolation of 1961–97 daily temperature and precipitation data onto Alberta polygons of ecodistrict and soil landscapes of Canada. J. Appl. Meteorol. Climatol. 2001, 40, 2162–2177. [Google Scholar] [CrossRef]
  32. Hess, J.J.; LM, S.; Knowlton, K.; Saha, S.; Dutta, P.; Ganguly, P.; Tiwari, A.; Jaiswal, A.; Sheffield, P.; Sarkar, J.; et al. Building resilience to climate change: Pilot evaluation of the impact of India’s First Heat Action Plan on All-Cause Mortality. J. Environ. Public Health 2018, 1, 7973519. [Google Scholar] [CrossRef]
  33. Heo, S.; Nori-Sarma, A.; Lee, K.; Benmarhnia, T.; Dominici, F.; Bell, M.L. The use of a quasi-experimental study on the mortality effect of a heat wave warning system in Korea. Int. J. Environ. Res. Public Health 2019, 16, 2245. [Google Scholar] [CrossRef]
  34. Royé, D.; Codesido, R.; Tobías, A.; Taracido, M. Heat wave intensity and daily mortality in four of the largest cities of Spain. Environ. Res. 2020, 182, 109027. [Google Scholar] [CrossRef]
  35. Williams, L.; Erens, B.; Ettelt, S.; Hajat, S.; Manacorda, T.; Mays, N. Evaluation of the Heatwave Plan for the UK: Final Report; Policy Innovation Evaluation Research Unit, London School of Hygiene and Tropical Medicine: London, UK, 2019. [Google Scholar]
  36. Conlon, K.C.; Mallen, E.; Gronlund, C.J.; Berrocal, V.J.; Larsen, L.; O’Neill, M.S. Mapping human vulnerability to extreme heat: A critical assessment of heat vulnerability indices created using principal components analysis. Environ. Health Perspect. 2020, 128, 097001. [Google Scholar] [CrossRef] [PubMed]
  37. Klinenberg, E. Heat Wave: A Social Autopsy of Disaster in Chicago; University of Chicago Press: Chicago, IL, USA, 2015. [Google Scholar]
  38. Ebi, K.L.; Teisberg, T.J.; Kalkstein, L.S.; Robinson, L.; Weiher, R.F. Heat watch/warning systems save lives: Estimated costs and benefits for Philadelphia 1995–98. Bull. Am. Meteorol. Soc. 2004, 85, 1067–1074. [Google Scholar] [CrossRef]
  39. Fouillet, A.; Rey, G.; Wagner, V.; Laaidi, K.; Empereur-Bissonnet, P.; Le Tertre, A.; Frayssinet, P.; Bessemoulin, P.; Laurent, F.; De Crouy-Chanel, P.; et al. Has the impact of heat waves on mortality changed in France since the European heat wave of summer 2003? A study of the 2006 heat wave. Int. J. Epidemiol. 2008, 37, 309–317. [Google Scholar] [CrossRef]
  40. Yang, J.; Yin, P.; Sun, J.; Wang, B.; Zhou, M.; Li, M.; Tong, S.; Meng, B.; Guo, Y.; Liu, Q. Heatwave and mortality in 31 major Chinese cities: Definition, vulnerability and implications. Sci. Total Environ. 2019, 649, 695–702. [Google Scholar] [CrossRef]
Figure 1. Socioeconomic maps used in the study.
Figure 1. Socioeconomic maps used in the study.
Ijerph 21 00815 g001
Figure 2. Spatial pattern of health outcomes in Korea (2007–2020).
Figure 2. Spatial pattern of health outcomes in Korea (2007–2020).
Ijerph 21 00815 g002
Figure 3. Nonlinear relationship between Z-scores of PCs and morbidity and mortality (Solid blue lines represent the morbidity and mortality contribution rates. The dashed blue lines represent the 95% confidence level).
Figure 3. Nonlinear relationship between Z-scores of PCs and morbidity and mortality (Solid blue lines represent the morbidity and mortality contribution rates. The dashed blue lines represent the 95% confidence level).
Ijerph 21 00815 g003
Figure 4. Comparison of heat vulnerability indices.
Figure 4. Comparison of heat vulnerability indices.
Ijerph 21 00815 g004
Table 1. Socioeconomic variables considered for constructing heat vulnerability maps in the study.
Table 1. Socioeconomic variables considered for constructing heat vulnerability maps in the study.
TypeVariableDescriptionMeanRangePeriod
Populationpop_65Population aged 65 and higher (%)19.67.9–33.52007–2020
agri_65Percentage of agricultural workers among the population aged 65 and higher (%)41.90.5–65.42007–2020
single_65Percentage of single-person households among the population aged 65 and higher (%)22.814.7–29.52010–2020
Economiclowincome_65Percentage of low-income households among the population aged 65 and higher (%)7.74.5–10.62010–2020
financialMunicipal fiscal self-reliance (%)25.110.4–57.92007–2020
Climate& Environmentt9090th percentile temperature (°C) for each region (2007–2020)30.326.6–32.42007–2020
tmax95Number of days with temperatures above the 95th percentile (2007–2020) compared to those with 90th percentile temperature (°C) in summer (JJA) from 1999 to 20205.83.6–6.92007–2020
tavgAverage temperature (°C) in summer24.221.9–25.62007–2020
rhav80Number of days with relative humidity exceeding 80% in summer (JJA)38.716.4–74.82007–2020
rhav_avgAverage relative humidity (%) in summer (JJA)77.068.7–86.02007–2020
per_forestForest occupies the total area (%)60.025.8–86.02007–2020
Table 2. Summary of principal component analysis results.
Table 2. Summary of principal component analysis results.
ComponentInitial Eigenvalues
Total% of VarianceCumulative
14.3239.2939.29
22.8525.9465.23
31.6314.8180.04
41.968.6888.72
50.635.6994.41
60.252.2496.65
70.151.3798.02
80.110.9999.01
90.050.5099.51
100.040.3899.89
110.010.11100.00
Table 3. Matrix of the first four principal components.
Table 3. Matrix of the first four principal components.
VariablePC1PC2PC3PC4
pop_650.9060.2380.160−0.191
agri_650.8440.2360.227−0.314
single_650.8630.3030.2360.015
lowincome_650.4080.200−0.2090.804
financial−0.927−0.185−0.025−0.077
t90−0.2720.8320.395−0.064
tmax950.0360.6650.1480.321
tavg−0.5160.3310.7390.059
rhav800.452−0.7990.3300.118
rhav_avg0.500−0.7430.3670.109
per_forest0.4710.374−0.697-0.171
The colored areas are where high impact occurred.
Table 4. p-values of principal components in the regression models.
Table 4. p-values of principal components in the regression models.
ModelOutcomePC1PC2PC3PC4Adjusted R2
LMMorbidity<0.001<0.001<0.0010.1550.381
Mortality0.0780.0060.0150.0010.154
NLMMorbidity<0.0010.001<0.0010.1840.518
Mortality0.1500.264<0.001<0.0010.486
Statistically significant parts are marked in bold.
Table 5. Assigning values of the heat vulnerability index.
Table 5. Assigning values of the heat vulnerability index.
MethodImpactRange of ComponentsAssigned HVI
PC1PC2PC3PC4
HVI-LMMorbidity<−2.326−44−4
−2.326–−1.645−33−3
−1.645–−1.282−22−2
−1.282–−0.674−11−1
−0.674–0.674000
0.674–1.2821−11
1.282–1.6452−22
1.645–2.3263−33
2.326<4−44
Mortality<−1.960 4−44
−1.960–−1.645 3−33
−1.645–−1.150 2−22
−1.150–−0.675 1−11
−0.675–0.675 000
0.675–1.150 −11−1
1.150–1.645 −22−2
1.645–1.960 −33−3
1.960 < −44−4
HVI-NLMMorbidity<−1.960−34−4
−1.960–−1.645−23−3
−1.645–−1.150−12−2
−1.150–−0.67501−1
−0.675–0.675 00
0.675–1.1500−11
1.150–1.6452−22
1.645–1.9603−33
1.960 <4−44
Mortality<−1.960 −44
−1.960–−1.645 −33
−1.645–−1.150 −22
−1.150–−0.675 11
−0.675–0.675 00
0.675–1.150 1
1.150–1.645
1.645–1.960
1.960<
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

Jeong, S.; Lim, Y.; Kang, Y.; Yi, C. Elucidating Uncertainty in Heat Vulnerability Mapping: Perspectives on Impact Variables and Modeling Approaches. Int. J. Environ. Res. Public Health 2024, 21, 815. https://doi.org/10.3390/ijerph21070815

AMA Style

Jeong S, Lim Y, Kang Y, Yi C. Elucidating Uncertainty in Heat Vulnerability Mapping: Perspectives on Impact Variables and Modeling Approaches. International Journal of Environmental Research and Public Health. 2024; 21(7):815. https://doi.org/10.3390/ijerph21070815

Chicago/Turabian Style

Jeong, Sockho, Yeonyeop Lim, Yeji Kang, and Chaeyeon Yi. 2024. "Elucidating Uncertainty in Heat Vulnerability Mapping: Perspectives on Impact Variables and Modeling Approaches" International Journal of Environmental Research and Public Health 21, no. 7: 815. https://doi.org/10.3390/ijerph21070815

APA Style

Jeong, S., Lim, Y., Kang, Y., & Yi, C. (2024). Elucidating Uncertainty in Heat Vulnerability Mapping: Perspectives on Impact Variables and Modeling Approaches. International Journal of Environmental Research and Public Health, 21(7), 815. https://doi.org/10.3390/ijerph21070815

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