Heat-Attributable Deaths between 1992 and 2009 in Seoul, South Korea

Background Climate change may significantly affect human health. The possible effects of high ambient temperature must be better understood, particularly in terms of certain diseases’ sensitivity to heat (as reflected in relative risks [RR]) and the consequent disease burden (number or fraction of cases attributable to high temperatures), in order to manage the threat. Purpose This study investigated the number of deaths attributable to abnormally high ambient temperatures in Seoul, South Korea, for a wide range of diseases. Method The relationship between mortality and daily maximum temperature using a generalized linear model was analyzed. The threshold temperature was defined as the 90th percentile of maximum daily temperatures. Deaths were classified according to ICD-10 codes, and for each disease, the RR and attributable fractions were determined. Using these fractions, the total number of deaths attributable to daily maximum temperatures above the threshold value, from 1992 to 2009, was calculated. Data analyses were conducted in 2012–2013. Results Heat-attributable deaths accounted for 3,177 of the 271,633 deaths from all causes. Neurological (RR 1.07; 95% CI, 1.04–1.11) and mental and behavioral disorders (RR 1.04; 95% CI, 1.01–1.07) had relatively high increases in the RR of mortality. The most heat-sensitive diseases (those with the highest RRs) were not the diseases that caused the largest number of deaths attributable to high temperatures. Conclusion This study estimated RRs and deaths attributable to high ambient temperature for a wide variety of diseases. Prevention-related policies must account for both particular vulnerabilities (heat-sensitive diseases with high RRs) and the major causes of the heat mortality burden (common conditions less sensitive to high temperatures).


Introduction
The progression of climate change has led to a concomitant worldwide increase in the incidence of extreme heat spells. Following several particularly dramatic examples, such as the heat waves in Chicago (1995) and Western Europe (2003), many studies have been published on high temperatures and human health. [1][2][3][4] It has been well documented that exposure to elevated ambient temperatures significantly increases the risk of mortality from cardiovascular problems. [5][6][7][8] The normal human body temperature ranges from 36.1°C to 37.8°C and is maintained by the hypothalamus in the brain. When the degree of environmental heat exceeds the regulatory capacity of the hypothalamus, the body's core temperature rises and exerts substantial stress on the organs, particularly the cardiovascular system. [6,[8][9][10] Most studies on the relationship between high temperature and mortality have focused on overall death and examined temperature-related cardiovascular and respiratory conditions, since these account for the majority of fatalities. [6,7,[11][12][13][14][15] While some evidence suggests a substantial effect of high temperatures on mortality rate, [13,14,[16][17][18][19][20] the relationship between heat and specific causes of death has not yet been examined in a metropolitan city in Asia.
In South Korea, an extreme heat watch/warning system was launched in 2007. These heat warning systems mainly concentrate on servicing the elderly, nursing homes, and childcare professionals. Present results show that heat affects a broad range of diseases and conditions, suggesting that there needs to be further development of inclusive heat-associated mortality prevention policies and programs for specific diseases.
Prior studies on temperature and health employed a variety of statistical designs and modeling, such as case-control, [1,2,9] time series, [6,12,[21][22][23] and case-crossover [18,21,22,24] methods. There has been extensive research on estimating threshold temperatures and applying relative risk (RR) derived from time series analysis. [11,15,23,[25][26][27] Of the many previous epidemiological analyses on high temperatures, few have used both relative risk and attributable death to evaluate the effect of high temperature on mortality according to disease state. [18,19,[28][29][30][31]"Attributable death" refers to the number of deaths that would be averted if a given risk factor was absent. [32] Estimating the number of deaths attributable to high temperatures for specific diseases will enable policymakers to anticipate the health gains that can be achieved from the implementation of preventative measures. [4,33,34] Unlike previous studies that focused intensively on identifying the RRs of highly prevalent diseases, the present study investigates not just the RRs, but a number of heat-attributable deaths that occurred in conjunction with cause-specific disease in Seoul, South Korea. Quantifying both RRs and attributable deaths will provide an understanding of the high temperature effect and the burden of mortality attributable to heat.

Study Area
Seoul is the largest city in South Korea, with a reported population of 9,794,304 in 2010 (S1 Fig.). The total land area of the South Korea is 99,720 square kilometers. The territories of Seoul are about 605.17 square kilometers. It is the same size as the city and county of San Francisco. The population density of Seoul is 17,473 person per square kilometers in 2010. It is more than 2.5 times the population density of San Francisco (6,664 person per square kilometer). [35] Seoul has a humid, subtropical climate and four distinct seasons. The East Asian monsoon affects the weather throughout the year, [36] approaching from the southeast region of the Korean Peninsula during the summer (June to August) and bringing hot and humid weather with abundant precipitation (accounting for more than 60% of annual rainfall). [36,37] The "urban heat island effect" is particularly marked in Seoul because the main commercial and industrial sectors are located there. [38] Over the past century, the average annual temperature has risen by 1.7°C in South Korea, greater than the increase in global average temperatures (0.74°C). [39] In Korea, heat wave advisories have been issuing warnings from summer 2007. The warning is issued if the temperature is expected to remain above 35 degrees Celsius for two or more consecutive days. The Korean climate change assessment reported that consecutive days of very high temperatures are increasing in magnitude and duration. [40] Mortality Data This study examined the influence of temperature on mortality by disease category. The daily mortality data from 1992 to 2009 were provided by Statistics Korea [35] (Table 1). Deaths were classified in accordance with International Classification of Disease, 10 th Revision (ICD-10) codes (Table 1).

Environmental Data
Measurements of relative humidity, air pressure, and ambient temperature from 1992 to 2009 were obtained from the Korea Meteorological Administration (KMA) [41]. Three-hour maximum temperature, mean relative humidity, and air pressure were measured from the representative synoptic surface observation station in Seoul (S1 Fig.). To control for potential confounding effects, 24-hour mean concentration of particulate matter (PM10) and 8-hour maximum concentration of ozone (O3) were obtained from the Korean National Institute of Environmental Research. Concentrations of air pollutants were measured every 15 min at 27 monitoring stations in Seoul.

Statistical Analysis
Graphs that depict the relationship between temperature and mortality generally take the shape of a U, V, or J, and show an increase in mortality rate above a specific threshold temperature. [15,23,42,43] The modeling approaches were adopted from well-established studies on air pollution [44,45] and high ambient temperature. [15,28,43] An overdispersed Poisson generalized linear model (GLM) was used to analyze the relationship between mortality and daily maximum temperature from 1992 to 2009. [45,46] PM10 (24-hour average), O3 level (daily maximum), relative humidity, day of the week, air pressure, and long-term time trends were selected as confounding variables. The threshold temperature was set at 90% of the maximum daily temperature during the study period, which was consistent with the findings from a previous study in Seoul.
[11] S1 Table shows the number of days over the study period that had daily maximum temperatures higher than the threshold value. The analysis was based on the measurements taken from May to September of each year.
Daily attributable fractions and RRs were estimated using the threshold temperature (29.5°C) as the reference point, and then the total number of deaths attributable to high temperatures from 1992 to 2009 in Seoul were calculated. [19,28] The attributable fraction of each disease on days higher than the threshold temperature was computed as follows: Attributable fraction (AF) of deaths AF = (RR − 1)/RR, [19,47] The number of attributable deaths relative to the total number of deaths was also added to the result.

Sensitivity Analyses
Sensitivity analyses were performed to explore the possible influences of PM10 and O3 on attributable deaths; however, this was only investigated for the 2001-2009 period (the only time period for which data were available; S1 Table). In addition, the RRs were compared using different lag days: single-day lag 1, lag 2, and moving average lag 0-2. The threshold temperature was also varied, being set at 93%, 95%, and 99% of the maximum daily temperature. All of these analyses were performed using SAS version 9.3 (SAS Institute Inc., USA) and R version 2.14 for Windows (http://cran.r-project.org/). All analyses were conducted in 2012-2013.

Ethics Statement
This study qualified for exemption by the Institutional Review Board of the Graduate School of Public Health at Seoul National University. (IRB No. 65-2013-12-05)

Results
During the study period, 271,633 (40%) of the 676,509 deaths occurred from May to September. Table 1 shows the descriptive statistics for cause-specific mortality in Seoul, based on ICD-10 codes. The main cause of death was cardiovascular disease (CVD). Table 2 summarizes the meteorological data. The average daily maximum temperature from May to September was 26.5°C, and 16.9°C for the full year. The pattern of daily all-cause mortality and temperature shows seasonal trends. All-cause mortality has a J-shaped relationship with temperature (S2 Fig.). The 1994 summer was the hottest, and there was some variation in each year. However, the average temperature remained stable (S3 Fig.). In 1994, temperature may have contributed to a higher number of deaths. Although there was temperature variation with regard to time and season, our statistical model controlled for the unmeasured time-varying potential confounder and showed the short-term temperature effect on mortality. Fig. 1 displays the relations between daily maximum temperature and mortality. Above a threshold temperature of 29.5°C, a rise in temperature of 1°C resulted in an increase in death from all causes (RR 1.03; 95% confidence interval (CI) 1.02-1.03). This included cardiovascular conditions, (RR 1.04; 95% CI, 1.03-1.04); respiratory conditions (RR 1.02; 95% CI, 1-1.04);   The RR and attributable deaths were compared, revealing that most heat-sensitive conditions (those with the highest RR) tended to be relatively uncommon causes of death, and the numbers of attributable deaths were less than for the major disease groups such as CVD (Fig. 2). The two LOWESS (locally weighted scatterplot smoothing) lines (span fraction = 1.0) of main-category and subcategory diseases were added. Both lines indicated a negative relationship between RRs and attributable deaths. Non-significant causes were excluded in this Fig.

Sensitivity Analysis
Several additional analyses were performed to explore the sensitivity of the estimated RR to model specification. Specifically, sensitivity to 1) the inclusion of air pollutants, 2) exposure lag, and 3) threshold temperature were examined.
In addition to the adverse health effects of high temperatures, measured air pollution variables should be considered. In this study, adjusting air pollution in the model yielded nearly the same estimates as did the non-adjusted model (S2 Table). Consequently, the models used in this study are believed to have adequately captured the main effects of temperature on mortality.
The lag in the time series models was also changed to single lag 1, 2, and moving average lag 0-2 days, which yielded similar results (S3 Table). For all-cause deaths, the estimated attributable deaths of lag 1 was the largest. The lag 1 estimates for each disease showed a pattern similar to that of lag 0, but were slightly larger. Estimates for endocrine, nervous, and genitourinary system conditions were slightly larger for lag 0.
In addition, four different temperatures, which corresponded to the 90th, 93th, 95th, and 99th percentiles for daily maximum temperature during the 19 years, were compared. The number of days over the study period where the maximum temperature was higher than the threshold value is shown in S1 Table. When the threshold temperature increases, the RR estimates increase (S4 Table). Although the RR increases, the corresponding days in which the temperature is higher than the threshold decreases, so the number of attributable days appears to decrease when the threshold increases.

Discussion
Several previous epidemiological studies have reported a correlation between temperature changes and disease-related mortality. In this study, the RRs above the threshold temperature and the number of heat-attributable deaths for each type of disease from 1992 to 2009 were calculated, particularly including sub-categories of the major disease groups. Results showed that almost all diseases were associated with a greater mortality risk during periods of unusual heat. The proportion of deaths attributable to days with a maximum temperature greater than 29.5 degrees ranged from 0.83% to 3.27%.
The present investigation builds on previous work on mortality in Seoul, and may be the first study conducted in an East Asian metropolitan city that examined the effects of heat on disease-specific mortality. Several studies that were conducted in cities in Europe and the US have identified specific populations vulnerable to mortality during high ambient temperatures. [18,19,24,48] The present study found that high temperature was linked with increases in mortality for most primary causes of death, including cardiovascular, respiratory, endocrine, nutritional, metabolic, neurological, genitourinary, digestive, mental, and behavioral disorders.
The observations herein are similar to those reported previously on the relationship between total mortality and ambient temperature. In this study, the results indicated that the overall mortality risk increased by 2.7% for every 1°C increase above threshold temperature. In a study conducted in Wales and 10 regions of London, total mortality risk increased by 2.1% (95% CI 1.6-2.6) for every 1°C increase in threshold temperature. [19] A study conducted by Chung et al. in Seoul reported a 2.7% increase (95% CI, 2.2-3.1) in daily non-accidental mortality, using a threshold temperature of 30.1-33.5°. [11] Note that this earlier study in Seoul used a different temperature metric (apparent temperature), included all seasons of the year, and was based on a shorter period of observation (1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006); however, the present results are still similar.
Also similar to previous reports, the present results show that the mortality rate from cardiovascular conditions increases with high ambient temperatures (3.5%; 95% CI, 2.6-4.5). Specific cardiovascular causes of death that are significantly affected by high temperatures include sudden death (4.4%), cerebrovascular diseases (4.3%), hypertensive diseases (3.7%), and ischemic heart disease (2.1%).
Data on the relationship between the mortality rate from respiratory conditions and high ambient temperatures have been inconsistent. Several studies reported no significant change in the number of respiratory deaths as a function of high temperature during the summer or in hot cities. [6,11,48,49] The present results do show that total respiratory mortality increases when the temperature rises, but an effect of high temperature on mortality due to pneumonia and chronic obstructive pulmonary disease (COPD) was not found in the present study. However, it is clear that deaths from asthma were sensitive to high temperatures.
The results of this study reveal a relatively high increase in the RR of mortality from neurological (7.4%) and mental and behavioral disorders (4.1%). This increase in mortality may occur when extremely high temperatures induce neurological dysfunction and multi-organ failure. [8,9,[50][51][52][53] However, other risk factors cannot be overlooked. It is well-known that patients taking antipsychotic or psychotropic medications may be vulnerable to high temperatures, and the health effects observed in this study might be partially drug-related. [3,9,[53][54][55] Additionally, a significant heat-associated increase in the RR of mortality from genitourinary conditions was observed (5%). Thus, the present results show that patients with a wide range of pre-existing chronic conditions may be more susceptible to high ambient temperatures.
Unlike previous studies that focused solely on identifying the RRs of highly prevalent diseases, the present study determined the RRs and attributable deaths for ICD-10 subcategory diseases. Some of these conditions, especially asthma and psychoactive-substance-use-related mental disorders, were found to be particularly heat sensitive. Furthermore, there was a negative relationship between RR and attributable death (log scale) in this study. The most heatsensitive diseases (those with the highest RRs) were not the diseases that caused the largest numbers of deaths attributable to high temperatures.
There were several limitations of the present study. First, the maximum daily temperature (lag 0) was used in the final model. Although sensitivity analyses using different lag days (single lag 1, single lag 2, and moving average lag 0-2) were performed, the results did not substantially change (S3 Table). High ambient temperatures generally had an immediate effect, but it is possible that lag patterns differ between diseases. [26,56] Second, threshold temperature used herein corresponded to the 90th percentile of the maximum daily temperature, which has been used previously. [11,19,57] In Japan, researchers examined the optimum temperature at which mortality is the lowest, observing that this point was 80-85% of the daily maximum temperature. [57] However, there are alternate methods of estimating the threshold temperature. [27,58,59] One of the most commonly used methods is grid search. It searches for the point that minimizes deviation among figures of the simulated model. Some studies compare complex models that assume a non-linear relationship between mortality and temperature. [23,28] Depending on the method used, the observed heat effect can vary. Moreover, the threshold temperature may differ for each disease because the underlying mechanisms by which temperature influences mortality is likely vary by disease. Despite this, using the percentage distribution method has been reported as effective and having advantages compared with alternate methods and their outcomes. [28,56,57,60] The sensitivity analysis results on different percentiles of the maximum daily temperature were included (S4 Table). There were no substantial change in the results between the 93 rd and 95 th percentiles.
Third, the effect of air pollutants for the entire study period could not be analyzed because air pollution data were not available until 2001. To estimate how the missing information on air pollutants may have affected the results another sensitivity analysis was conducted, after adjusting for this information-the results of the adjusted model and unadjusted model were robust (S2 Table). Thus, air pollution did not appear to be a concern.
Finally, despite the consistency in the results, there may be some contextual variability in the effects of temperature on mortality across locations within Seoul and among different populations, levels of social development, and individual sensitivities and adaptive capacities. Due to the urban heat island effect, there is also temperature variability within a metropolitan region that is not measured in this study because temperature was recorded at only one weather station.

Conclusion
This study estimated the RRs and attributable deaths related to high ambient temperature and included a wide variety of diseases. The present results indicate that high ambient temperature increases mortality for many conditions and that the effect varies with each disease. Some relatively uncommon conditions-that is, those with a low number of attributable deaths-appear to be highly sensitive to heat, warranting close attention by interventions that are designed to manage the risk of high temperature. In this, as in other areas of public health, policymakers must consider both vulnerable sub-groups and the predominant causes of disease burden.

Author Contributions
Conceived and designed the experiments: CK YL AW HK. Performed the experiments: CK YL AW HK. Analyzed the data: CK YL HK. Contributed reagents/materials/analysis tools: CK HK. Wrote the paper: CK YL AW HK.