Seasonal impact of diurnal temperature range on intracerebral hemorrhage in middle-aged and elderly people in central China
Article information
Abstract
OBJECTIVES
This study investigated the seasonal impact of diurnal temperature range (DTR) on hospitalization rates for intracerebral hemorrhage (ICH) in middle-aged and elderly adults.
METHODS
We collected data on the DTR and hospitalization records of ≥45-year-old patients with ICH in 2019 in Hunan Province, central China. Time-series analyses were performed using a distributed lag non-linear model.
RESULTS
Overall, 54,690 hospitalizations for ICH were recorded. DTR showed a non-linear relationship with ICH hospitalization in both middle-aged and elderly populations (45-59 and ≥60 years, respectively). During spring, a low DTR coupled with persistently low temperatures increased ICH risk in both age groups, while a high DTR was associated with an increased risk in the middle-aged group only (relative risk [RR], 1.24; 95% confidence interval [CI], 1.21 to 1.27). In the summer, a low DTR combined with persistently high temperatures was linked to a higher risk exclusively in the middle-aged group. A high DTR in the autumn was correlated with increased risk in both age groups. In winter, either a low DTR with a continuously low temperature or a high DTR elevated the risk solely in the elderly population (RR, 1.37; 95% CI, 1.00 to 1.69). In the elderly group, the impact of DTR on hospitalization risk manifested within a 5-day period.
CONCLUSIONS
The impact of DTR on ICH hospitalization risk differed significantly across seasons and between age groups. Elderly individuals demonstrated greater sensitivity to the impact of DTR. Weather forecasting services should emphasize DTR values, and interventions targeting sensitive populations are needed.
Key Message
In a large-scale, study of 54,690 hospitalizations for intracerebral hemorrhage (ICH) from China, we identified the seasonal impact of diurnal temperature range (DTR) on hospitalization rates for intracerebral hemorrhage (ICH) in middle aged and elderly adults. Specifically, the impact of DTR on ICH hospitalization risk differed significantly across seasons and between age groups. Elderly individuals demonstrated greater sensitivity to the impact of DTR. The findings of this study could be instrumental in the prevention and control of ICH.
INTRODUCTION
Intracerebral hemorrhage (ICH) is a clinical syndrome characterized by the non-traumatic rupture of blood vessels within the brain, resulting in the accumulation of blood in the brain parenchyma. This can lead to cognitive impairment, mental and emotional disorders, and other symptoms [1,2]. ICH is known for its rapid onset and potentially serious consequences, with the highest mortality rate among all stroke types [3,4]. According to 2019 data on the global burden of disease, the incidence of ICH in China is between 49.4 and 71.9 per 100,000 people. While this rate is the highest in the world, it represents a decrease of 2.1% compared to 1990. However, the mortality rate has increased by 14.4%, making ICH the stroke subtype with the greatest mortality in China [5,6]. ICH is particularly prevalent among middle-aged and elderly individuals [7,8]. In 2018, the total direct medical costs for ICH in China amounted to 20.512 billion yuan per year, which accounted for 3.91% of the nation’s total medical expenses—an enormous financial burden [9-11]. Therefore, the early prevention of ICH is crucial.
The diurnal temperature range (DTR) reflects the variation in temperature within a single day and is defined as the difference between the daily maximum and minimum temperatures [12]. It represents a key meteorological indicator of global climate change. As with temperature, exposure to DTR is unavoidable, and its link to the development of ICH is an important public health issue [13-18]. Studies have examined the association between ICH mortality and DTR, suggesting that a higher DTR may elevate the risk of mortality in patients with ICH [19-22]. Accordingly, most research has concentrated on mortality outcomes. However, the growing burden and costs associated with stroke care underscore the pressing need for effective prevention strategies [23-28]. Nonetheless, the evidence pertaining to the relationship between DTR and the onset of or hospitalization for ICH remains ambiguous, largely due to inadequate control for confounding factors such as ambient temperature, limited geographic scope of studies, and insufficient consideration of varying temperature exposures and physical adaptations among subgroups [29-33]. The impact of DTR on hospitalizations for ICH warrants further study [34,35], as do the health implications of DTR across different seasons and subpopulations [36-38].
Hunan Province, a representative region of central China, spans from 108°47´ to 114°15´ east longitude and from 24°38´ to 30°08´ north latitude. In Changsha, the capital city of Hunan, the incidence of ICH is 131.0 per 100,000 persons, which is substantially higher than the global average [39]. The region experiences 4 distinct seasons, with sufficient heat, concentrated rainfall, variable spring temperatures, dry summers and autumns, a brief but intense cold period, and an extended period of summer heat. These climatic conditions provide Hunan with a geographic advantage for studying the effects of DTR on ICH hospitalization rates [40].
We employed a distributed lag non-linear model to examine the influence of DTR on ICH. Our analysis centered on DTR as the primary exposure, incorporating its interaction with temperature in different seasons as well as additional variables, such as social factors, within a population of middle-aged and elderly individuals in central China.
MATERIALS AND METHODS
Data collection
The study included a total of 54,690 patients aged 45 years or older who were admitted to secondary (554 in total) or tertiary (103 in total) hospitals with a primary diagnosis of ICH (I61) between January 1, 2019 and December 31, 2019. These data were obtained from the direct statistical reporting system of the Hunan Health Commission. The geographic distribution of the secondary and tertiary hospitals is depicted in Supplementary Material 1. Daily meteorological data were obtained from the China Meteorological Data Network (https://data.cma.cn/) and categorized into 4 seasons: spring (March to May), summer (June to August), autumn (September to November), and winter (December to February) [12]. Daily air quality data were sourced from the National Air Pollution Monitoring System (http://envi.ckcest.cn/environment/data_Integration/data_Integration.jsp). Information on the population, economy, and health resources for each district/county in Hunan Province for the year 2019 was extracted from the Hunan Statistical Yearbook 2020 (http://tjj.hunan.gov.cn/).
For admitted patients, the inclusion criteria encompassed: (1) a hospitalization date between January 1, 2019 and December 31, 2019; (2) a primary discharge diagnosis of ICH, under the International Classification of Diseases, 10th revision (code I61); (3) an age of 45 years or older (patients 45 to 59 years old were categorized as middle-aged, while those 60 years old or older were considered elderly) [41]; and (4) a minimum hospital stay of 1 day. The exclusion criteria were: (1) an address that was either unclear or not located in Hunan Province and (2) ICH resulting from trauma. The data collection process is illustrated in Figure 1.
Exposure assessment
Hunan Province is equipped with 29 environmental/meteorological monitoring stations and 71 air quality monitoring stations. To obtain daily average meteorological factors and air pollution levels for each county or district, spatiotemporal Kriging interpolation analysis was employed. This method was used to supplement the daily meteorological and air quality data for each district and county of Hunan Province for the year 2019.
Data checking and screening included several processes. (1) Logical verification involved removing admissions with clear logical errors and duplications. For example, individuals with hospitalization costs of 0 were excluded, and multiple admissions of a patient on the same day were counted as a single admission. Otherwise, admissions were treated as separate, independent events. (2) Regarding missing data, admissions lacking essential information (such as present residence, household registration address, or discharge diagnosis information) were omitted from the analysis. (3) The Q-test was employed to identify univariate outliers. Any outlier data identified were then verified and either removed or adjusted accordingly.
Statistical analysis
In this study, percentiles were used to describe meteorological factors and air pollutant levels (Supplementary Material 2). We tested the correlation between air pollutant concentrations and meteorological factors (Supplementary Material 3). Any pair of variables with a correlation coefficient greater than 0.8 was considered strongly correlated, and a single variable was selected for inclusion in the model according to the clinical significance and goodness-of-fit. Additionally, we conducted an over-dispersion test; this yielded a p-value of less than 0.05, indicating that the data were over-dispersed. Consequently, a quasi-Poisson regression was employed to fit the time-series model.
The selection of independent variables was informed by a literature review and a priori knowledge from experts about factors known to be associated with time, as well as their presence in the dataset. The controlled potential covariates included [42-45]: (1) environmental factors: diurnal mean temperature (denoted in the model as “temp”), relative humidity (“relahumidity”), particulate matter with diameters of 2.5 micrometers and smaller (PM2.5), particulate matter with diameters of 10 micrometers and smaller (PM10), nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), and ozone (O3); and (2) social factors: gross regional product (“GDP”), average population (“popu”), the number of medical institutions in each district/county in 2019 (“hospitals”), day of the week (“DOW”), and public holidays (“holiday”).
Based on the data presented, we applied 2 sets of time-series quasi-Poisson models to estimate the associations between DTR and ICH hospitalization, utilizing the distributed lag non-linear model framework. In the first time-series quasi-Poisson model analysis—model (1)—we adjusted for seasonality to estimate the overall impact of DTR on hospitalization for ICH among middle-aged and elderly populations. Correlation analysis indicated that PM2.5 was strongly correlated with PM10 (rs= 0.83, p<0.05). Given the practical significance and goodness-of-fit, PM10 was excluded from this analysis. The initial model was formulated as follows:
Here, α represents the intercept, β symbolizes the log-relative risk of hospitalization associated with a unit increase in DTR, t denotes the day of the year (from 1 to 365), i is the code number for the district/county, and l represents the lag in days. E(Yt,i) indicates the expected count of hospitalizations in area i on day t. DTRt,i,l is a cross-basis matrix that assesses the non-linear and lag effects of DTR on ICH [46]. Time denotes long-term effects; accordingly, we set the degrees of freedom (df) to 12 [47]. Consistent with previous research, our primary model incorporated a natural cubic spline of calendar days with 12 df per year to adjust for seasonality and long-term trends. We tested different df for time, ranging from 4 to 16, and found that 12 yielded the lowest Akaike information criterion (AIC). The definitions of the other variables are available in the previous text. The nodes of the exposure-response relationship were established at the 25th, 50th, and 75th percentiles. Based on the AIC, we determined the maximum lag days to be 14. The environmental variables, including relative humidity, diurnal mean temperature, PM2.5, PM10, NO2, SO2, CO, and O3, were modeled using natural cubic spline functions (ns), with the df set at 3.
In the second time-series quasi-Poisson model, labeled model (2), we conducted separate analyses for each age group in the different seasons. Specifically, we calculated the relative risk (RR) by using the DTR associated with the lowest risk as a reference to estimate the effect of DTR during each season. We constructed models individually for the 2 age groups in the 4 seasons. The methodology was as follows: (1) the median DTR of each season was used as the initial reference for model fitting; (2) the DTR corresponding to the lowest RR was identified for each season and served as the updated reference; and (3) the model was refitted using this updated reference DTR to estimate the RR.
The model for each season was expressed as follows. The meanings of the variables are consistent with model (1).
Unless otherwise specified, all p-values represent 2-sided probabilities, with the significance level set at α=0.05. To ensure the robustness of our main results, we conducted a sensitivity analysis by varying the df (from 3 to 5) and the maximum lag days (from 14 to 21) for DTR. Additionally, sensitivity analyses were performed by altering the number of times of admission (1, 2, > 2). All analyses were conducted using the “dlnm” package (version 2.4.7) in R version 4.2.1 (R Foundation for Statistical Computing, Vienna, Austria). The names, identification numbers, and telephone numbers of patients were desensitized during the data collection phase. We adhered to data confidentiality principles throughout the study.
Ethics statement
This study is an ecological analysis of hospitalization records. Patient names, identification numbers, and telephone numbers were desensitized during the data collection phase. Data confidentiality principles were followed throughout the study. Approval for this study was granted by the Xiangya School of Public Health, Central South University (approval No. XYGW-2019-040).
RESULTS
Table 1 presents the daily hospital admissions for middle-aged and elderly individuals with ICH in Hunan Province during the year 2019. The total number of admissions was 54,690, corresponding to a hospitalization rate of 0.24%. Of these admissions, 16,824 were middle-aged patients, representing 44.43% of the total and exhibiting a hospitalization rate of 0.13%. The remaining 37,866 patients were elderly, comprising 55.57% of the admissions and representing a hospitalization rate of 0.40%.
The ICH hospitalization rate and the number of admissions were higher in the elderly than in the middle-aged population. Supplementary Material 4 presents the average length of ICH hospitalization in Hunan in 2019. Male and elderly individuals experienced longer hospital stays, with the duration of ICH hospitalizations being longest during the winter for all groups.
Table 2 and Supplementary Material 5 indicate that the highest number of ICH hospitalizations among the middle-aged and elderly populations occurred in January, with March and November also displaying high admission numbers. Hospital admissions were more frequent in spring and winter and less frequent in summer. Additionally, the monthly variation in admission person-times was greater in the elderly group compared to the middle-aged group, with a range of 2,154 compared to 875, respectively.
Supplementary Material 2 presents the daily meteorological and air quality data for 122 districts/counties in Hunan Province for 2019. Supplementary Material 3 displays the correlation results between these variables. Correlation analysis indicated that DTR was positively correlated with daily average temperature (rs=0.34, p < 0.05). Additionally, PM2.5 was strongly correlated with PM10 (rs=0.83, p<0.05). Due to considerations of practical significance and goodness-of-fit, PM10 was excluded from this analysis.
Supplementary Material 6 presents a line graph depicting the monthly average temperature and DTR in Hunan Province for 2019. The DTR trend exhibited a bimodal pattern, with peaks in April and October.
Figure 2 illustrates the overall impact of DTR on ICH among the middle-aged and elderly populations, suggesting that a higher DTR may elevate the risk of hospitalization for ICH in both age groups. Notably, in the elderly group, the risk of ICH decreased as DTR increased when the range was within 4°C.
Table 3 and Figure 3 illustrate the cumulative impact of DTR on ICH hospitalization among middle-aged and elderly populations in different seasons. Table 3 reveals that during spring, a higher DTR was associated with a greater risk of hospitalization for both middle-aged and elderly people. When the DTR was low (less than 3°C for middle-aged and less than 4°C for elderly individuals), both age groups also faced an elevated risk of hospitalization. In the summer, a low DTR was associated with a higher risk of ICH hospitalization for middle-aged individuals. Conversely, in autumn, a high DTR was linked to a greater risk of ICH hospitalization for both middle-aged and elderly populations. During winter, the risk of ICH hospitalization in the elderly rose regardless of whether the temperature difference was high (greater than 3°C) or low (less than 3°C).
Table 4 presents the cumulative lag effects of extreme DTR on hospitalization for ICH in middle-aged individuals, categorized by season. The P2.5 and P97.5 values represent extremely low and extremely high DTR values, respectively. In the spring, only the extremely low DTR exhibited a lag effect in this population. The initial effect persisted until lag0-4 days and subsided; the effect then reemerged at lag0-11 (RR, 1.21; 95% CI, 1.02 to 1.48), with a peak at lag0-12 (RR, 1.32; 95% CI, 1.11 to 1.42). During the summer, both extremely high and extremely low DTR began to influence hospitalizations starting at lag0-14. No lag effects of extreme DTR were observed in autumn or winter. Overall, middle-aged individuals should be aware of the lag effect occurring approximately 14 days after extreme DTR in the summer and about 11 days following extremely low DTR in the spring.
Table 4 presents the cumulative lag effect of extreme DTR on ICH hospitalization among the elderly, with P2.5 and P97.5 again representing extremely low and extremely high DTR, respectively. In spring, the impact of extremely low DTR began to appear at lag0-5 days (RR, 1.20; 95% CI, 1.01 to 1.36), and continued to increase until lag0-14. In comparison, the influence of extremely high DTR diminished as the number of lag days increased, disappearing by lag0-9 (RR, 1.23; 95% CI, 0.99 to 1.38). Regarding summer, the lag effects of both extremely low and extremely high DTR were statistically insignificant. In autumn, the effect of an extremely high DTR was observed at lag0-14. In winter, the impact of an extremely high DTR began at lag0-1 (RR, 1.21; 95% CI, 1.04 to 1.42) and grew until lag0-10 (RR, 1.39; 95% CI, 1.04 to 1.64), exhibiting a peak before declining. The effect of extremely low DTR started at lag0-4 and persisted until lag0-14. Overall, among the elderly, the heightened risk of ICH associated with DTR in the spring and winter manifested rapidly and lasted for approximately 1 week to 2 weeks.
The main results of the sensitivity analyses remained robust when we adjusted the df for DTR from 3 to 5 (Supplementary Material 7). Consistent findings were also observed when we varied the maximum lag for DTR from 14 days to 21 days (Supplementary Material 8) and when we modified the times of admission (Supplementary Material 9).
DISCUSSION
This study presents a comprehensive analysis of the effects of daily fluctuations in DTR on hospitalizations for ICH. Our research is unique in that it centers on DTR as the primary exposure, considering its interplay with seasonal temperatures and variables such as social factors among middle-aged and elderly populations in central China. The findings of this study could be instrumental in the prevention and control of ICH.
We observed that the most prominent difference in hospitalizations occurred between December 31 and January 1. This discrepancy may stem from the annual settlement of medical insurance and the impact of holidays at the end of the year, which could influence reporting practices. Consequently, hospitalizations on January 1 were markedly higher than those on December 31. Our findings also indicate that the effects of DTR on ICH hospitalizations in the middle-aged demographic are predominantly delayed, manifesting after an extended period (> 10 days), except for low DTR during spring (< 2°C). The lag in the impact of DTR on ICH onset may be attributed to the time required for certain high-risk factors to develop or change [35]. The risk incrementally increases with rising DTR in spring and autumn, suggesting that the body’s thermoregulatory system may struggle to swiftly adapt to these environmental shifts [33]. In spring, a low DTR (< 3°C) coupled with sustained low temperatures provoked a rapid increase in the risk of ICH hospitalization among middle-aged individuals. Similarly, a persistent high-temperature environment with low DTR during summer elevated the risk of hospitalization for ICH in this age group. Potential underlying mechanisms include elevated blood pressure, peripheral vasoconstriction and dehydration, changes in salt intake, and heightened blood viscosity due to continuous exposure to extreme temperatures, whether low or high [34,38]. This implies that comprehensive health education on ICH, along with timely weather forecasts that highlight DTR, could effectively mitigate the risk of brain hemorrhage [34]. Additionally, the establishment of a high-temperature health alert system, greater utilization of air conditioning, avoidance of prolonged outdoor labor, and vigilant monitoring of health status may reduce the likelihood of hospitalization for ICH in middle-aged individuals [44].
The impact of DTR on ICH hospitalizations in the elderly resembled that in the middle-aged group in spring, but the effect manifested more quickly after exposure. The diminished thermoregulatory function in elderly people makes them particularly susceptible to temperature fluctuations, and the higher prevalence of chronic diseases in this population may further increase their sensitivity. Consequently, the risks associated with DTR manifested more rapidly in the elderly than in the middle-aged group [27,38]. Early warning systems and medical services must offer convenient access and provide prompt support for the elderly. This population should also be vigilant about staying warm to prevent exposure to sustained low temperatures in spring and winter. Additionally, caregivers should monitor weather forecasts and promptly assist elderly individuals to mitigate their risk [14]. In autumn, a high DTR may increase the risk of ICH hospitalization among the elderly. Therefore, among older adults, it is crucial to wear proper clothing, limit time spent outdoors, and maintain a healthy living environment with appropriate temperatures to reduce the risk of ICH [2].
This study has several limitations. First, the exact time of ICH onset prior to hospitalization was unknown, which means that using hospital admission as the primary outcome could result in an inaccurate assessment of the impact of DTR. However, patients with ICH are usually admitted to the hospital immediately for treatment. The median time from onset to admission has been reported to be 6.25 hours, with an interquartile range of 2.5 hours to 24.0 hours [47]. Additionally, individuals with mild or asymptomatic ICH may not always seek hospitalization, potentially leading to an underestimation of the true incidence of ICH. Despite this, ICH is responsible for a high proportion of stroke-related mortality and morbidity, and few demonstrably effective treatments are available for its acute management or prevention. The proportion of severe admissions is much higher than that of mild or asymptomatic admissions. Therefore, hospitalization is a reasonable proxy for estimating incidence based on the available data [49,50]. Second, the environmental exposure measurement used in this study is an average estimate for the local population and does not account for individual variations in exposure. This limitation hampers the ability to control for individual-level confounding factors, indicating a need for further research. Finally, the ecological nature of this study means that its findings cannot be used to establish a causal relationship.
After adjusting for other confounding factors, we observed a distinct non-linear relationship between DTR and hospitalization for ICH among middle-aged and elderly populations. The influence of DTR on ICH hospitalization risk varied significantly across seasons and between age groups, with the elderly participants demonstrating greater sensitivity to DTR. We suggest that weather forecasting services should report and highlight DTR values, and targeted interventions should be developed for populations that are particularly vulnerable to these changes.
Supplementary materials
Supplementary materials are available at https://doi.org/10.4178/epih.e2024053.
Notes
Data availability
The datasets used and/or analyzed during this study are available from the corresponding author on reasonable request.
Conflict of interest
The authors have no conflicts of interest to declare for this study.
Funding
This research was funded as part of the project of the Natural Science Foundation of Hunan Province, China (grant Nos. 2020 JJ4772 and 2021JJ30893) and the National Science & Technology Foundational Resource Investigation Program of China (grant Nos. 2018FY100900 and 2018FY100901).
Author contributions
Conceptualization: Wang S, Shi J, Zhou H. Data curation: Wang S, Zhou H, Shi Q, Yang F. Formal analysis: Wang S, Zhou H. Funding acquisition: Shi J, Deng J. Project administration: Shi J. Writing – original draft: Wang S. Writing – review & editing: Wang S, Yin J, Zhou H, Lai J, Xiao G, Tong Z, Deng J, Yang F, Shi Q, Shi J.
Acknowledgements
We thank the Information Statistics Center of Health Commission of Hunan Province for their support with data access and all agencies who provided the meteorological data.