Skip to main content
Advertisement
  • Loading metrics

Association between temperature variability and daily hospital admissions for cause-specific cardiovascular disease in urban China: A national time-series study

  • Yaohua Tian,

    Roles Conceptualization, Formal analysis, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Hui Liu,

    Roles Methodology, Visualization, Writing – review & editing

    Affiliations Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China, Medical Informatics Center, Peking University, Beijing, China

  • Yaqin Si,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliations Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China, Beijing HealthCom Data Technology Co. Ltd, Beijing, China

  • Yaying Cao,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Jing Song,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Man Li,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Yao Wu,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Xiaowen Wang,

    Roles Software, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Xiao Xiang,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Juan Juan,

    Roles Formal analysis, Writing – original draft

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Libo Chen,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Beijing HealthCom Data Technology Co. Ltd, Beijing, China

  • Chen Wei,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Beijing HealthCom Data Technology Co. Ltd, Beijing, China

  • Pei Gao ,

    Contributed equally to this work with: Pei Gao, Yonghua Hu

    Roles Conceptualization, Data curation, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    yhhu@bjmu.edu.cn (YH); peigao@bjmu.edu.cn (PG)

    Affiliations Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China, Key Laboratory of Molecular Cardiovascular (Peking University), Ministry of Education, Beijing, China

  • Yonghua Hu

    Contributed equally to this work with: Pei Gao, Yonghua Hu

    Roles Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

    yhhu@bjmu.edu.cn (YH); peigao@bjmu.edu.cn (PG)

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

Abstract

Background

Epidemiological studies have provided compelling evidence of associations between ambient temperature and cardiovascular disease. However, evidence of effects of daily temperature variability on cardiovascular disease is scarce and mixed. We aimed to examine short-term associations between temperature variability and hospital admissions for cause-specific cardiovascular disease in urban China.

Methods and findings

We conducted a national time-series analysis in 184 cities in China between 2014 and 2017. Data on daily hospital admissions for ischemic heart disease, heart failure, heart rhythm disturbances, and ischemic stroke were obtained from the database of Urban Employee Basic Medical Insurance (UEBMI) including 0.28 billion enrollees. Temperature data were acquired from the China Meteorological Data Sharing Service Center. Temperature variability was calculated from the standard deviation (SD) of daily minimum and maximum temperatures over exposure days. City-specific associations between temperature variability and cardiovascular disease were examined with overdispersed Poisson models controlling for calendar time, day of the week, public holiday, and daily mean temperature and relative humidity. Random-effects meta-analyses were performed to obtain national and regional average associations. We also plotted exposure-response relationship curve using a natural cubic spline of temperature variability. There were 8.0 million hospital admissions for cardiovascular disease during the study period. At the national-average level, a 1-°C increase in temperature variability at 0–1 days (TV0–1) was associated with a 0.44% (0.32%–0.55%), 0.31% (0.20%–0.43%), 0.48% (0.01%–0.96%), 0.34% (0.01%–0.67%), and 0.82% (0.59%–1.05%) increase in hospital admissions for cardiovascular disease, ischemic heart disease, heart failure, heart rhythm disturbances, and ischemic stroke, respectively. The estimates decreased but remained significant when controlling for ambient fine particulate matter (PM2.5), NO2, and SO2 pollution. The main limitation of the present study was the unavailability of data on individual exposure to temperature variability.

Conclusions

Our findings suggested that short-term temperature variability exposure could increase the risk of cardiovascular disease, which may provide new insights into the health effects of climate change.

Author summary

Why was this study done?

  • Epidemiological studies have reported associations between daily temperature and mortality or morbidity from cardiovascular disease.
  • Evidence of effects of temperature variability on cardiovascular disease is scarce and mixed.

What did the researchers do and find?

  • We conducted a national time-series study using data on 8.0 million hospital admissions for cause-specific cardiovascular disease in 184 cities in China from 2014 to 2017.
  • City-specific associations between temperature variability and hospital admissions for cardiovascular disease were estimated with overdispersed Poisson models. Random-effects meta-analyses were applied to obtain national and regional average associations.
  • Our results showed that temperature variability was associated with increased hospital admissions for ischemic heart disease, heart failure, heart rhythm disturbances, and ischemic stroke.

What do these findings mean?

  • To our knowledge, this is the first national study in China to examine the associations between temperature variability and daily hospital admissions for cause-specific cardiovascular disease.
  • Our findings suggested that general practitioners, cardiologists, and public health officials should consider the detrimental effects of temperature variability in their mitigation strategies.

Introduction

Cardiovascular disease is a major cause of death and disability worldwide [1,2]. Climate change has been considered potentially the greatest threat to human health of the 21st century [35]. In the past few decades, increasing epidemiological studies have reported associations of mortality and morbidity from cardiovascular disease with ambient temperature [6,7]. However, daily mean temperature failed to reflect intra- and interday variations of temperature. Temperature variability is an important meteorological indicator reflecting climate changes, such as rapid temperature fluctuations within a certain period (e.g., intra- and interday changes in temperature) [8,9]. It may also pose a major danger to human health. To date, however, scientific evidence is still inadequate regarding the association between daily temperature variation and cardiovascular disease morbidity, especially in developing countries. Quantifying the effect of temperature variability on cardiovascular disease has major public health implications, particularly in the context of increasing climate change and unstable weather patterns [10,11].

In China, climate change has emerged as a significant public health problem [12]. Recently, several studies have examined the association between temperature variability and daily mortality [1315]. However, a common limitation of these studies was their focus on a single city or several sites in a local area. Previous studies have demonstrated that climate conditions, sociodemographic status, and locations of study sites could modify the health effects of temperature variability [8,16]. China has a vast territory and varied terrain, as well as a great diversity of climates, including tropical in southern regions to subarctic in the extreme north. Therefore, the associations between short-term temperature variability and cause-specific cardiovascular disease morbidity remain unclear at the national level. A national analysis is important to comprehensively assess the health effects of temperature variation in China.

In this study, we designed a national time-series study to examine the short-term associations between temperature variability and hospital admissions of cause-specific cardiovascular disease for 2014–2017 in urban China as well as to investigate individual-level and city-level potential effect modifiers.

Methods

Study sites

A total of 184 cities were included in this analysis. The inclusion of these cities was on the basis of the availability of both weather and health data. Cities with records of less than 1 year were excluded due to the feasibility of model fit. Individuals’ detailed information of the disease diagnosis was required to identify the cardiovascular admissions. Cities with no information on International Classification of Diseases (ICD) code or those whose text of disease diagnosis cannot be classified as categories of cardiovascular disease were also excluded. S1 Fig presents the locations of the 184 cities, representing a geographic distribution across China.

Data collection

There are three main health insurance programs in China: the Urban Employee Basic Medical Insurance (UEBMI) for urban working and retired employees, the Urban Residence Basic Medical Insurance for urban residents without formal employment, and the New Rural Cooperative Medical Scheme for rural residents. These three programs covered more than 92% of the population by 2011 [17].

Daily counts of hospital admissions within each city were obtained from the database of UEBMI. In 2016, there were 0.28 billion beneficiaries in 31 provincial administrative regions, representing approximately one-fifth of the total population in Mainland China. The size of this population allows for investigation of cause-specific cardiovascular diseases that have been associated with temperature variability [18]. The number of city residents, people enrolled in the database, and coverage rate of the population by UEBMI are presented in S1 Table. The number of city’s residents is according to the 2010 census. Among the 184 cities, Jiayuguan is the smallest city with a population of 23.19 thousand, and Chongqing is the largest city with a population of 28.8 million. The medical information recorded on the database includes age, sex, date of medical service, diagnosis, and financial cost. Hospital admissions for each health condition were identified based on the primary diagnoses (text of disease diagnosis or ICD code). Cardiovascular admissions included ischemic heart disease (ICD-10 codes I20–I25), heart failure (I50), heart rhythm disturbances (I47–I49), and ischemic stroke (I63) [19]. Admissions aged <18 years were too few and therefore were not included. Because the data used for this study were collected for administrative purpose without any individual identifiers, this study was exempted from Institutional Review Board approval by the Ethics Committee of Peking University Health Science Center, Beijing, China. The need for informed consent was also waived by the Institutional Review Board.

Data on temperature and relative humidity for each city were acquired from the China Meteorological Data Sharing Service Center. There are one to three monitoring stations in each city. We calculated daily maximum, minimum, and mean temperatures and mean relative humidity by averaging all valid monitoring measurements in each city. During the study period, the missing rate of temperature and relative humidity was only 0.25%. To control for the potential confounding effects of ambient air pollution, we obtained data on fine particulate matter (PM2.5), NO2, and SO2 concentrations for each city from the National Air Pollution Monitoring System. There are 1 to 17 monitors in each city. We calculated daily mean levels of PM2.5, NO2, and SO2 by averaging all valid monitoring measurements in each city. During the study period, the missing rates were 0.66% for PM2.5, 0.15% for NO2, and 0.15% for SO2, respectively. Days with missing information on meteorological variables or the air pollution data were excluded from analysis.

Exposure definition

In this study, we employed a composite indicator of intraday and interday temperature variability to assess the effects of temperature variation on daily hospital admissions. The definition of the temperature variability indicator has been described previously [8]. Temperature variability index was calculated from the standard deviation (SD) of daily minimum and maximum temperatures over the exposure days. For example, temperature variability for the preceding 2 days’ exposure (TV0–1) was calculated as SD (MinTemplag0, MaxTemplag0, MinTemplag1, MaxTemplag1) such that MinTemplagi and MaxTemplagi represent the minimum and the maximum temperature for the preceding day i; temperature variability for the preceding 3 days’ exposure was calculated by TV0–2 = SD (MinTemplag0, MaxTemplag0, MinTemplag1, MaxTemplag1, MinTemplag2, MaxTemplag2), and so on. This index could account for both intraday and interday temperature variability, as well as the delayed effects of temperature variation [8].

Statistical analysis

The associations between temperature variability and daily hospital admissions for cardiovascular disease were estimated by a commonly used two-stage approach [8,9]. The method and the model used in this study were designed before the analyses were conducted, and the prespecified statistical analysis plan was present in Supporting Information (S1 Appendix). In the first stage, quasi-Poisson regression models were applied for estimating city-specific estimates of hospital admissions associated with temperature variability exposure. Previous studies have observed that mortality risk is almost linearly elevated with increases in temperature variability [8,9,20]. We then designed a linear term for temperature variability to be included in the model. Furthermore, consistent with previous studies [8,9,18], several confounders were included in the model, as follows:

  1. A natural cubic spline of calendar time with 7 degrees of freedom (df) per year to adjust for unmeasured time trends longer than 2 months in hospital admissions
  2. A natural cubic spline of 3-day moving average relative humidity with 3 df
  3. Indicator variables for day of the week and public holiday
  4. To adjust for the confounding effect of daily mean temperature, a distributed lag nonlinear model (DLNM) for daily mean temperature was also included in the model [8,9,21].

Specifically, we used two natural cubic spline functions with 4 df, respectively, for daily mean temperature and lag over time up to 21 days to accommodate the nonlinear and lagged effects of ambient temperature. The selection of 21 days was based on the evidence that the effect of cold temperature was delayed and lasted for several weeks, while the effect of hot temperature was acute and generally presented within 1 week [22,23]. Three internal knots were placed at equally spaced temperature percentiles (25th, 50th, and 75th). Consequently, the model was as shown below: such that E(Yt) is the expected count of admissions on day t, β represents the log-relative risk of admission associated with a unit increase of temperature variability, “Temperature” indicates a two-dimensional cross-basis matrix produced by DLNM, and ns() indicates natural cubic spline function. Public holidays and day of the week were included in the model as indicator variables, and “relative humidity” indicates 3-day moving average relative humidity. In the second stage, the national and regional average estimates of the associations between temperature variability and hospital admissions were obtained by combining city-specific estimates using a random-effects meta-analysis. I2 statistics were calculated to examine the heterogeneity of city-specific estimates. We conducted separate analyses for temperature variability during different exposure days from TV0–1 to TV0–7 [8,18].

To quantify the uncertainty in heterogeneity estimates, we also calculated 95% confidence interval (CI) for I2 [24]. To assess the shape association between temperature variability and the outcome, we used the same approach as the “Air Pollution and Health: A European Approach” (APHEA) project and a recent multicity analysis on the association between temperature variability and mortality [9,25,26]. A natural cubic spline was applied. We also investigated potential effect modifiers of the association between temperature variability and daily hospital admissions for cardiovascular disease. We conducted stratified analyses by age, sex, and geographical region. Due to the regional differences in weather conditions, geography, and culture, the 184 cities were grouped into south and north regions [27,28]. The P values for difference in estimates between the strata were obtained via Z-test [29].

In addition, we explored the confounding influence of air pollutants or potential effect modifications by cities’ characteristics. To assess the confounding effect of air pollutants, the model was further adjusted for the linear terms of daily PM2.5, NO2, and SO2 concentrations. To evaluate the potential effect modifications by cities’ characteristics, city-specific relative risk (and CIs) as the outcome were meta-regressed on each continuous variable of city characteristics. City-level characteristics included annual-average temperature variability levels, temperature and relative humidity, gross domestic product (GDP) per capita, and coverage of population by UEBMI. We performed random-effects meta-regression model using city-level estimations from the first stage as: yi = α + βxi + μi + ei such that yi is the city-specific estimation, xi is the city-level characteristic variable, μi is a normal error term with known SDs within-city standard error that may vary across units, and ei is a normal error with variance to be estimated, assumed equal across units [30]. The model was implemented by the “metareg” function in STATA.

Sensitivity analysis

Sensitivity analyses were conducted by changing the maximum lag of temperature from 21 to 28 days, df for daily mean temperature (3–6), and df for time trend (6–8 per year), respectively. We also assessed the robustness of the results after controlling for daily minimum or maximum temperature instead of daily mean temperature. In addition, we used penalized spline functions for time trend and meteorological variables. Finally, we reported estimated effects using interquartile range (IQR) increase of temperature variability.

Statistical analyses were performed in R software version 3.2.2 (R Foundation for Statistical Computing, Vienna, Austria) and STATA software version 12 (StataCorp, College Station, TX). The results were presented as percentage change and 95% CI in hospital admissions per 1-°C increase in temperature variability. Percentage change equals relative risk minus 1 and then multiplied by 100.

Results

Table 1 presents the demographic characteristics of people enrolled in the UEBMI in the 184 Chinese cities in 2017. Of the 184 cities, 94 cities were located in south region and 90 cities were in the north. Overall, there were 54.4% male patients and 4.9% patients aged ≥75 years. The sex and age distributions of participants in the south versus north regions were similar. For the 184 cities, there were 8.0 million hospital admissions for cardiovascular disease from January 1, 2014, through December 31, 2017. S2 and S3 Tables present city-specific summary statistics on annual-average hospital admissions for cause-specific cardiovascular disease, weather variables, temperature variability, and air pollutants. Summary statistics of daily hospital admissions for cardiovascular disease, weather conditions, temperature variability, and air pollutants are presented in Table 2. On average, we recorded 44 cardiovascular admissions per day across the cities, with a range from 1 to 302. The annual-average daily mean temperature was 14.4°C, with a range from −0.2°C to 24.4°C. The annual-average TV0–1 was 5.7°C, ranging from 2.9°C to 9.3°C. The distributions of temperature variability were similar at different exposure days (TV0–1 to TV0–7). The distribution of daily temperature variability at different exposure days is presented in S4 Table.

thumbnail
Table 1. Demographic characteristics of people enrolled in the UEBMI in the 184 Chinese cities in 2017.

https://doi.org/10.1371/journal.pmed.1002738.t001

thumbnail
Table 2. Summary statistics of daily hospital admissions for cardiovascular disease, weather conditions, air pollutants, and temperature variability in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1002738.t002

Fig 1 presents the national-average exposure-response curve between TV0–1 and hospital admissions for cardiovascular disease. We noted a broadly linear association. The curve was relatively flat at low levels (<3°C) and increased rapidly at levels greater than 3°C.

thumbnail
Fig 1. National-average exposure-response association curve between TV0–1 and daily hospital admissions for cardiovascular disease in 184 cities in China, 2014–2017.

TV0–1, temperature variability at 0–1 days.

https://doi.org/10.1371/journal.pmed.1002738.g001

Fig 2 summarizes the national-average percentage changes in hospital admissions for cardiovascular diseases associated with 1-°C increase in temperature variability at different exposure days. We observed weak or moderate between-city heterogeneity for the associations between temperature variability and hospital admissions for cardiovascular disease (I2 = 47.8%; 95% CI 38.0%–56.0%), ischemic heart disease (I2 = 36.5%; 24.1%–47.0%), heart failure (I2 = 3.8%; 0%–15.1%), heart rhythm disturbances (I2 = 21.2%; 5.1%–34.5%), and ischemic stroke (I2 = 51.4%; 42.1%–59.3%). The largest estimates were observed at exposure 0–1 days for all health outcomes except heart failure. For a 1-°C increase in TV0–1, we observed significant increases of 0.44% (0.32%–0.55%) in hospital admissions for cardiovascular disease, 0.31% (0.20%–0.43%) for ischemic heart disease, 0.48% (0.01%–0.96%) for heart failure, 0.34% (0.01%–0.67%) for heart rhythm disturbances, and 0.83% (0.60%–1.06%) for ischemic stroke.

thumbnail
Fig 2. National-average PC with 95% CI in daily hospital admissions for cause-specific cardiovascular disease per 1-°C increase in temperature variability at different exposure days in 184 Chinese cities, 2014–2017.

CI, confidence interval; PC, percentage change; TV0–1, temperature variability at 0–1 days.

https://doi.org/10.1371/journal.pmed.1002738.g002

Fig 3 shows the results of the stratified analyses. The associations between TV0–1 and cardiovascular disease varied by age. For cardiovascular disease, the estimates were larger in people aged 65–74 years and ≥75 years than in people aged 18–64 years (0.81 [0.59–1.03] in aged ≥75 versus 0.55 [0.34–0.75] in aged 65–74 versus 0.19 [0.03–0.34] in aged 18–64 group [P < 0.05]), similar to in ischemic heart disease and ischemic stroke (0.75 [0.52–0.98] in aged ≥75 versus 0.08 [−0.07 to 0.23] in aged 18–64 group [P < 0.001] for ischemic heart disease and 0.90 [0.60–1.21] in aged ≥75 versus 0.32 [0.08–0.56] in aged 18–64 group [P = 0.003] for ischemic stroke). The estimates were broadly similar between males and females (Fig 3, all P > 0.05), between northern and southern regions (Fig 3, all P > 0.05), and among different climate temperature groups (S5 Table). We divided the cities into four groups based on the quantiles of their annual-average temperatures: cold area (≤25th), moderate cold area (25th–50th), moderate hot area (50th–75th), and hot area (>75th). We observed consistent significant associations in different climate areas. A 1-°C increase in TV0–1 corresponded to a 0.38% (0.15%–0.62%), 0.60% (0.34%–0.86%), 0.43% (0.13%–0.73%), and 0.55% (0.18%–0.93%) higher cardiovascular disease admissions in cold, moderate cold, moderate hot, and hot areas, respectively. The estimates for all health conditions were higher in warm season than in cool season, although the difference was only significant for ischemic heart disease (P = 0.027).

thumbnail
Fig 3. National-average PC with 95% CI in daily hospital admissions for cause-specific cardiovascular disease per 1-°C increase in TV0–1, stratified by sex, age, and geographical region.

CI, confidence interval. CI, confidence interval; PC, percentage change; TV0–1, temperature variability at 0–1 days.

https://doi.org/10.1371/journal.pmed.1002738.g003

The associations between TV0–1 and cardiovascular disease after further adjustment of air pollutants are shown in Table 3. The estimates decreased but remained significant after adjusting for PM2.5 (0.39% change; 0.22%–0.56%), NO2 (0.19% change; 0%–0.38%), and SO2 (0.30% change; 0.13%–0.47%). Table 4 presents the results of effect modification by city-level characteristics. Annual-average temperature variability, temperature and relative humidity levels, GDP per capita, and coverage rate of population by UEBMI did not significantly modify the associations between TV0–1 and hospital admissions for cardiovascular disease (P > 0.05).

thumbnail
Table 3. Results on the associations between TV0–1 and hospital admissions for cardiovascular disease after controlling for the air pollutants.

https://doi.org/10.1371/journal.pmed.1002738.t003

thumbnail
Table 4. Meta-regression results of the modification effects of city-level characteristics on the associations between TV0–1 and hospital admissions for cardiovascular disease in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1002738.t004

Table 5 shows the results of sensitivity analyses. The association between TV0–1 and cardiovascular disease was insensitive to changing maximum lag of temperature (28 days) and df for time (6–8 per year) and temperature (3–6). The estimate slightly changed when using penalized spline functions in the model (0.48% change; 0.31%–0.65%). The association remained after adjustment of daily minimum and maximum temperature. We also applied a threshold approach to estimate the association and found that a 1-°C increase in TV0–1 was associated with a 0.48% (0.33%–0.64%) increase in hospital admissions for cardiovascular disease if temperature variability > 3°C. The estimates for associations between temperature variability and cause-specific cardiovascular diseases are presented in S6 Table. The estimated effects in IQR increase of TV0–1 (3.1°C) are presented in S7 Table.

thumbnail
Table 5. Results of sensitivity analyses on the associations between TV0–1 and hospital admissions for cardiovascular disease in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1002738.t005

Discussion

In this national study, we found that short-term temperature variability exposure was associated with increased hospital admissions for cardiovascular disease, after controlling for daily mean temperature. The associations were more evident in the elderly but did not vary substantially by climates (annual-average temperature variability, temperature and relative humidity levels), demographic characteristics (sex and GDP per capita), and geographical region. These associations were robust to the adjustment of air pollutants (PM2.5, NO2, and SO2).

A number of epidemiological studies have estimated the cardiovascular disease risk in association with intraday or interday temperature variation separately [20,31,32]. However, only capturing the impact of intraday or interday temperature variation might underestimate the true effect of temperature variability [8,18]. Few studies to date have considered both intra- and interday temperature variability in an attempt to comprehensively evaluate the acute effect of temperature variability on cardiovascular disease. A pooled estimate of 12 counties across Hubei Province in China reported a 1.72% (0.69%–2.76%) increase in cardiovascular disease mortality in relation to an increment of 1°C in TV0–7 [15]. A recent study with 31 major Chinese cities demonstrated that a 1-°C increase in TV0–7 was associated with a 0.60% (0.25%–0.94%) increase in cardiovascular disease mortality [18]. However, air pollution, an important environmental risk factor for cardiovascular disease, was not controlled for in these two studies. Therefore, the two studies might have overestimated the effect of temperature variability on cardiovascular disease. In this study, the estimates substantially decreased when controlling for PM2.5, NO2, and SO2, which was consistent with a recent study that reported a weakened association of temperature variability with cardiovascular disease mortality after adjusting for air pollution [14]. Specifically, to further explore the modification effect of NO2 on the association between TV0–1 and cardiovascular disease, we divided the cities into three groups based on the tertiles of their annual average NO2 levels. The association was stronger in cities with higher annual average NO2 levels (S8 Table). The association between temperature variability and cardiovascular admissions was attenuated but remained significant after adjustment of air pollutants, suggesting that air pollution is a confounder of the association.

Overall, we found significant associations between temperature variability and hospital admissions for cause-specific cardiovascular disease in China, and these associations were strongest at lag 0–1 days. Several studies have reported stronger effects of daily temperature variation on cardiovascular disease mortality at longer exposure days [14,15]. Differences in these findings might be attributable to different health outcomes. In this study, we utilized hospital admission data rather death data. Cardiovascular disease admissions can differ markedly from death events from cardiovascular disease by volume, severity, and demographics. In addition, cardiovascular disease admissions are more sensitive to short-term exposure to environmental factors and therefore could better explore the temporal pattern between temperature variability and clinical presentation of cardiovascular disease [33]. A recent multicountry study reported that the highest estimates of temperature-variability–mortality appeared at different exposure days in different countries [8]. Therefore, it was plausible to hypothesize that study sites’ characteristics, such as climates, geographical region, and population susceptibility, might also modify the lag pattern between temperature variability exposure and cardiovascular disease.

We observed a broadly linear exposure-response curve between temperature variability and cardiovascular disease, which is consistent with previous findings [8,9,20]. We also noted that the curve was relatively flat at low levels, in line with the curve for the association between temperature variability and mortality observed in Melbourne, Australia [9]. Human core body temperature is normally maintained within a narrow range via heat exchange with the surrounding environments. People may be able to adapt well to changes in temperature when they are relatively minor. However, people may have difficulty with thermoregulation to sudden, dramatic temperature changes not only physiologically but also in terms of behavioral patterns [3436]. Previous studies have demonstrated that unstable temperatures could alter heart rate, blood pressure, blood cholesterol level, inflammatory reaction, and immune function [3740]. These physiological changes may trigger cardiovascular disease [41].

In line with previous studies [15,18,20,42], we found that the elderly have a higher risk of cardiovascular disease during days with great temperature variation compared with the younger group. The decline of thermoregulatory function occurs naturally in aging adults [43], which might partly explain the increased vulnerability to temperature variability in the elderly. In addition, elderly people generally have a higher prevalence of chronic health conditions, which might also contribute to the increased vulnerability. China is experiencing rapid population aging, therefore the cardiovascular disease burden attributed to temperature variability is expected to increase. Moreover, unlike other risk factors, temperature variability exposure is ubiquitous. As a result, temperature-variability–related burden of cardiovascular disease should be an increasingly serious environmental health issue in China.

We observed weak or moderate between-city heterogeneity for the associations between temperature variability and hospital admissions for heart failure, heart rhythm disturbances, and ischemic heart disease, indicating that the effects of temperature variability on these health outcomes might not vary substantially by geographical and demographic characteristics in China. This hypothesis is further supported by the absence of evidence for effects modification by city’s characteristics, as shown in the meta-regression analyses. However, we observed relatively stronger heterogeneity of city-specific estimates for ischemic stroke. This might be attributable to the substantial regional variation in baseline incidence of ischemic stroke. A recent nationwide population-based study in China demonstrated significant regional variations in annual incidence rate of stroke [44]. City-specific characteristics, such as socio-demography, chronic diseases burden, or lifestyles, might also be responsible for the heterogeneity. Further investigation with more detailed information on individuals’ lifestyle risk factors is required to confirm this finding.

This study has several strengths. We analyzed multicity data using a uniform statistical method, thus enhancing internal comparison of results across these cities, as well as external comparison with other studies. Our study should provide more representative estimates of the associations between temperature variability and hospital admissions for cause-specific cardiovascular disease. In addition, we included 184 cities involving great diversity in temperature variability levels, meteorological conditions, and geography, enabling us to evaluate the potential effect modification by these characteristics. Finally, a range of analyses were performed to test the robustness of the findings, including additional adjustment of air pollution.

Several limitations should also be acknowledged. First, we used temperature data from fixed monitors as a surrogate of personal exposure, which can cause exposure measurement errors. This kind of nondifferential error is likely to bias the estimates downward [45]. Second, misclassification bias caused by diagnostic error should be considered when interpreting the results. However, this error is likely to be random and typically reduce the precision of the estimates and leads to underestimation of estimates [46]. Finally, children and adolescents (age group 0–18) were not included in this assessment because UEBMI is for the urban working and retired employees, though young children may have more difficulty in regulating their temperature, similarly to the elderly.

In conclusion, this national study in 184 Chinese cities, which is, to our knowledge, the largest study worldwide to date, provides evidence of a relationship between daily temperature variation and increased hospital admissions for cardiovascular disease. The findings have significant implications for general practitioners, cardiologists, and public health officials.

Supporting information

S1 Table. City-specific number of people enrolled in the UEBMI in 2017, the number of residents, and the coverage rate based on the UEBMI in 184 Chinese cities.

UEBMI, Urban Employee Basic Medical Insurance.

https://doi.org/10.1371/journal.pmed.1002738.s002

(DOCX)

S2 Table. Summary statistics on annual-average hospital admissions for cause-specific cardiovascular disease in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1002738.s003

(DOCX)

S3 Table. Summary statistics on annual-average TV0–1, weather conditions, and air pollutants in 184 Chinese cities, 2014–2017.

TV0–1, temperature variability at 0–1 days.

https://doi.org/10.1371/journal.pmed.1002738.s004

(DOCX)

S4 Table. The distribution of daily temperature variability at different exposure days in 184 cities in China, 2014–2017.

https://doi.org/10.1371/journal.pmed.1002738.s005

(DOCX)

S5 Table. National-average PC with 95% CI in daily hospital admissions for cardiovascular disease associated with a 1-°C increase in TV0–1 in 184 Chinese cities by climate region, 2014–2017.

CI, confidence interval; PC, percentage change; TV0–1, temperature variability at 0–1 days.

https://doi.org/10.1371/journal.pmed.1002738.s006

(DOCX)

S6 Table. National-average PC with 95% CI in daily hospital admissions for cause-specific cardiovascular disease per 1-°C increase in TV0–1 if temperature variability > 3°C in 184 Chinese cities, 2014–2017.

CI, confidence interval; PC, percentage change; TV0–1, temperature variability at 0–1 days.

https://doi.org/10.1371/journal.pmed.1002738.s007

(DOCX)

S7 Table. National-average PC with 95% CI in daily hospital admissions for cause-specific cardiovascular disease per IQR increase in TV0–1 (3.1°C) in 184 Chinese cities, 2014–2017.

CI, confidence interval; IQR, interquartile range; PC, percentage change; TV0–1, temperature variability at 0–1 days.

https://doi.org/10.1371/journal.pmed.1002738.s008

(DOCX)

S8 Table. National-average PC with 95% CI in daily hospital admissions for cardiovascular disease associated with a 1-°C increase in TV0–1 in 184 Chinese cities, 2014–2017, classified by city-specific annual average NO2 levels.

CI, confidence interval; PC, percentage change; TV0–1, temperature variability at 0–1 days.

https://doi.org/10.1371/journal.pmed.1002738.s009

(DOCX)

S1 Fig. The locations of 184 cities analyzed in this study.

The most populous city in each province was marked by name.

https://doi.org/10.1371/journal.pmed.1002738.s010

(TIF)

References

  1. 1. GBD 2016 Causes of Death Collaborators. Global, regional, and national age-sex specific mortality for 264 causes of death, 1980–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet. 2017;390(10100):1151–210. pmid:28919116
  2. 2. GBD 2016 DALYs and HALE Collaborators. Global, regional, and national disability-adjusted life-years (DALYs) for 333 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet. 2017;390(10100):1260–344. pmid:28919118
  3. 3. Haines A, Kovats RS, Campbell-Lendrum D, Corvalan C. Climate change and human health: impacts, vulnerability and public health. Public Health. 2006;120(7):585–96. pmid:16542689
  4. 4. Costello A, Abbas M, Allen A, Ball S, Bell S, Bellamy R, et al. Managing the health effects of climate change: Lancet and University College London Institute for Global Health Commission. Lancet. 2009;373(9676):1693–733. pmid:19447250
  5. 5. Patz JA, Campbell-Lendrum D, Holloway T, Foley JA. Impact of regional climate change on human health. Nature. 2005;438(7066):310–7. pmid:16292302
  6. 6. Phung D, Thai PK, Guo Y, Morawska L, Rutherford S, Chu C. Ambient temperature and risk of cardiovascular hospitalization: An updated systematic review and meta-analysis. Sci Total Environ. 2016;550:1084–102. pmid:26871555
  7. 7. Moghadamnia MT, Ardalan A, Mesdaghinia A, Keshtkar A, Naddafi K, Yekaninejad MS. Ambient temperature and cardiovascular mortality: a systematic review and meta-analysis. PeerJ. 2017;5:e3574. pmid:28791197
  8. 8. Guo Y, Gasparrini A, Armstrong BG, Tawatsupa B, Tobias A, Lavigne E, et al. Temperature Variability and Mortality: A Multi-Country Study. Environ Health Perspect. 2016;124(10):1554–9. pmid:27258598
  9. 9. Cheng J, Xu Z, Bambrick H, Su H, Tong S, Hu W. The mortality burden of hourly temperature variability in five capital cities, Australia: Time-series and meta-regression analysis. Environ Int. 2017;109:10–9. pmid:28923460
  10. 10. Epstein PR. Climate change and human health. N Engl J Med. 2005;353(14):1433–6. pmid:16207843
  11. 11. Hansen J, Sato M, Ruedy R. Perception of climate change. Proc Natl Acad Sci U S A. 2012;109(37):E2415–23. pmid:22869707
  12. 12. Kan H, Chen R, Tong S. Ambient air pollution, climate change, and population health in China. Environ Int. 2012;42:10–9. pmid:21440303
  13. 13. Zhang Y, Yu C, Yang J, Zhang L, Cui F. Diurnal Temperature Range in Relation to Daily Mortality and Years of Life Lost in Wuhan, China. Int J Environ Res Public Health. 2017;14(8):891.
  14. 14. Luo K, Li R, Wang Z, Zhang R, Xu Q. Effect modification of the association between temperature variability and daily cardiovascular mortality by air pollutants in three Chinese cities. Environ Pollut. 2017;230:989–99. pmid:28763936
  15. 15. Zhang Y, Yu C, Bao J, Li X. Impact of temperature variation on mortality: An observational study from 12 counties across Hubei Province in China. Sci Total Environ. 2017;588:196–203.
  16. 16. Lee W, Bell ML, Gasparrini A, Armstrong BG, Sera F, Hwang S, et al. Mortality burden of diurnal temperature range and its temporal changes: A multi-country study. Environ Int. 2018;110:123–30. pmid:29089167
  17. 17. Yip WC, Hsiao WC, Chen W, Hu S, Ma J, Maynard A. Early appraisal of China's huge and complex health-care reforms. Lancet. 2012;379(9818):833–42. pmid:22386036
  18. 18. Yang J, Zhou M, Li M, Liu X, Yin P, Sun Q, et al. Vulnerability to the impact of temperature variability on mortality in 31 major Chinese cities. Environ Pollut. 2018;239:631–7. pmid:29709834
  19. 19. Dominici F, Peng RD, Bell ML, Pham L, McDermott A, Zeger SL, et al. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA. 2006;295(10):1127–34. pmid:16522832
  20. 20. Zhan Z, Zhao Y, Pang S, Zhong X, Wu C, Ding Z. Temperature change between neighboring days and mortality in United States: A nationwide study. Sci Total Environ. 2017;585:1152–61.
  21. 21. Gasparrini A, Armstrong B, Kenward MG. Distributed lag non-linear models. Stat Med. 2010;29(21):2224–34. pmid:20812303
  22. 22. Guo Y, Gasparrini A, Armstrong B, Li S, Tawatsupa B, Tobias A, et al. Global variation in the effects of ambient temperature on mortality: a systematic evaluation. Epidemiology. 2014;25(6):781–9. pmid:25166878
  23. 23. Gasparrini A, Guo Y, Hashizume M, Lavigne E, Zanobetti A, Schwartz J, et al. Mortality risk attributable to high and low ambient temperature: a multicountry observational study. Lancet. 2015;386(9991):369–75. pmid:26003380
  24. 24. Ioannidis JP, Patsopoulos NA, Evangelou E. Uncertainty in heterogeneity estimates in meta-analyses. BMJ. 2007;335(7626):914–6. pmid:17974687
  25. 25. Samoli E, Analitis A, Touloumi G, Schwartz J, Anderson HR, Sunyer J, et al. Estimating the exposure-response relationships between particulate matter and mortality within the APHEA multicity project. Environ Health Perspect. 2005;113(1):88–95. pmid:15626653
  26. 26. Chen R, Yin P, Meng X, Liu C, Wang L, Xu X, et al. Fine Particulate Air Pollution and Daily Mortality. A Nationwide Analysis in 272 Chinese Cities. Am J Respir Crit Care Med. 2017;196(1):73–81. pmid:28248546
  27. 27. Almond D, Chen Y, Greenstone M, Hongbin L. Winter heating or clean air? Unintended impacts of China's Huai river policy. Am Econ Rev. 2009;99(2):184–90.
  28. 28. Ebenstein A, Fan M, Greenstone M, He G, Zhou M. New evidence on the impact of sustained exposure to air pollution on life expectancy from China's Huai River Policy. Proc Natl Acad Sci U S A. 2017;114(39):10384–9. pmid:28893980
  29. 29. Altman DG, Bland JM. Interaction revisited: the difference between two estimates. BMJ. 2003;326(7382):219. pmid:12543843
  30. 30. Knapp G, Hartung J. Improved tests for a random effects meta-regression with a single covariate. Stat Med. 2003;22(17):2693–710. pmid:12939780
  31. 31. Zhou X, Zhao A, Meng X, Chen R, Kuang X, Duan X, et al. Acute effects of diurnal temperature range on mortality in 8 Chinese cities. Sci Total Environ. 2014;493:92–7. pmid:24937494
  32. 32. Kim J, Shin J, Lim YH, Honda Y, Hashizume M, Guo YL, et al. Comprehensive approach to understand the association between diurnal temperature range and mortality in East Asia. Sci Total Environ. 2016;539:313–21. pmid:26363726
  33. 33. Villeneuve PJ, Chen L, Stieb D, Rowe BH. Associations between outdoor air pollution and emergency department visits for stroke in Edmonton, Canada. Eur J Epidemiol. 2006;21(9):689–700. pmid:17048082
  34. 34. Havenith G. Temperature regulation, heat balance and climatic stress. In: Kirch W, Menne B, Bertollini R, editors. Extreme weather events and public health responses. Berlin: Springer; 2005. pp. 69–80.
  35. 35. Martinez-Nicolas A, Meyer M, Hunkler S, Madrid JA, Rol MA, Meyer AH, et al. Daytime variation in ambient temperature affects skin temperatures and blood pressure: Ambulatory winter/summer comparison in healthy young women. Physiol Behav. 2015;149:203–11. pmid:26072176
  36. 36. Garrett AT, Goosens NG, Rehrer NJ, Patterson MJ, Cotter JD. Induction and decay of short-term heat acclimation. Eur J Appl Physiol. 2009;107(6):659–70. pmid:19727796
  37. 37. Keatinge WR, Coleshaw SR, Cotter F, Mattock M, Murphy M, Chelliah R. Increases in platelet and red cell counts, blood viscosity, and arterial pressure during mild surface cooling: factors in mortality from coronary and cerebral thrombosis in winter. BMJ. 1984;289(6456):1405–8. pmid:6437575
  38. 38. Keatinge WR, Coleshaw SR, Easton JC, Cotter F, Mattock MB, Chelliah R. Increased platelet and red cell counts, blood viscosity, and plasma cholesterol levels during heat stress, and mortality from coronary and cerebral thrombosis. Am J Med. 1986;81(5):795–800. pmid:3776986
  39. 39. Halonen JI, Zanobetti A, Sparrow D, Vokonas PS, Schwartz J. Associations between outdoor temperature and markers of inflammation: a cohort study. Environ Health. 2010;9(42):9–42.
  40. 40. Halonen JI, Zanobetti A, Sparrow D, Vokonas PS, Schwartz J. Outdoor temperature is associated with serum HDL and LDL. Environ Res. 2011;111(2):281–7. pmid:21172696
  41. 41. Feigin VL, Roth GA, Naghavi M, Parmar P, Krishnamurthi R, Chugh S, et al. Global burden of stroke and risk factors in 188 countries, during 1990–2013: a systematic analysis for the Global Burden of Disease Study 2013. Lancet Neurol. 2016;15(9):913–24. pmid:27291521
  42. 42. Ding Z, Li L, Xin L, Pi F, Dong W, Wen Y, et al. High diurnal temperature range and mortality: Effect modification by individual characteristics and mortality causes in a case-only analysis. Sci Total Environ. 2016;544:627–34. pmid:26674692
  43. 43. Li M, Gu S, Bi P, Yang J, Liu Q. Heat waves and morbidity: current knowledge and further direction-a comprehensive literature review. Int J Environ Res Public Health. 2015;12(5):5256–83. pmid:25993103
  44. 44. Wang W, Jiang B, Sun H, Ru X, Sun D, Wang L, et al. Prevalence, Incidence, and Mortality of Stroke in China: Results from a Nationwide Population-Based Survey of 480 687 Adults. Circulation. 2017;135(8):759–71. pmid:28052979
  45. 45. Guo Y, Barnett AG, Tong S. Spatiotemporal model or time series model for assessing city-wide temperature effects on mortality? Environ Res. 2013;120:55–62. pmid:23026801
  46. 46. Wellenius GA, Schwartz J, Mittleman MA. Air pollution and hospital admissions for ischemic and hemorrhagic stroke among medicare beneficiaries. Stroke. 2005;36(12):2549–53. pmid:16254223