Next Article in Journal
A Brief Up-Date of the Use of Sodium Oxybate for the Treatment of Alcohol Use Disorder
Next Article in Special Issue
Mortality during a Large-Scale Heat Wave by Place, Demographic Group, Internal and External Causes of Death, and Building Climate Zone
Previous Article in Journal
Accumulation of Heavy Metals in Vegetable Species Planted in Contaminated Soils and the Health Risk Assessment
Previous Article in Special Issue
Evaluating the Appropriateness of Downscaled Climate Information for Projecting Risks of Salmonella
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 13
Issue 3
10.3390/ijerph13030284
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Patterns of Heat-Related Cardiovascular Mortality in the Czech Republic

by
Aleš Urban
1,2,*,
Katrin Burkart
3,
Jan Kyselý
1,4,5,
Christian Schuster
6,
Eva Plavcová
1,
Hana Hanzlíková
1,7,
Petr Štěpánek
5,8 and
Tobia Lakes
6
1
Institute of Atmospheric Physics, Czech Academy of Sciences, Boční II 1401, 14131 Prague 4, Czech Republic
2
Faculty of Science, Charles University, Albertov 6, 12843 Prague 2, Czech Republic
3
Department of Environmental Health Science, Mailman School of Public Health, Columbia University, 722 W 168th Street, New York, NY 10032, USA
4
Faculty of Environmental Sciences, Czech University of Life Sciences, Kamýcká 129, 16521 Prague 6, Czech Republic
5
Global Change Research Centre, Czech Academy of Sciences, Bělidla 986, 60300 Brno, Czech Republic
6
Department of Geography, Geoinformation Science Lab, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany
7
Institute of Geophysics, Czech Academy of Sciences, Boční II 1401, 14131 Prague 4, Czech Republic
8
Czech Hydrometeorological Institute, Regional Office Brno, Kroftova 2578, 61667 Brno, Czech Republic
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2016, 13(3), 284; https://doi.org/10.3390/ijerph13030284
Submission received: 30 November 2015 / Revised: 31 January 2016 / Accepted: 23 February 2016 / Published: 4 March 2016
(This article belongs to the Special Issue Climate Change and Human Health)
Browse Figures
Versions Notes

Abstract

:
The study examines spatial patterns of effects of high temperature extremes on cardiovascular mortality in the Czech Republic at a district level during 1994–2009. Daily baseline mortality for each district was determined using a single location-stratified generalized additive model. Mean relative deviations of mortality from the baseline were calculated on days exceeding the 90th percentile of mean daily temperature in summer, and they were correlated with selected demographic, socioeconomic, and physical-environmental variables for the districts. Groups of districts with similar characteristics were identified according to socioeconomic status and urbanization level in order to provide a more general picture than possible on the district level. We evaluated lagged patterns of excess mortality after hot spell occurrences in: (i) urban areas vs. predominantly rural areas; and (ii) regions with different overall socioeconomic level. Our findings suggest that climatic conditions, altitude, and urbanization generally affect the spatial distribution of districts with the highest excess cardiovascular mortality, while socioeconomic status did not show a significant effect in the analysis across the Czech Republic as a whole. Only within deprived populations, socioeconomic status played a relevant role as well. After taking into account lagged effects of temperature on excess mortality, we found that the effect of hot spells was significant in highly urbanized regions, while most excess deaths in rural districts may be attributed to harvesting effects.

1. Introduction

Harmful effects of high atmospheric temperature on human health and lives have been documented in mid-latitudes [1,2] as well as in tropical climates [3,4]. Due to the rise in temperature with climate change, increased frequency and magnitude of extreme heat events are projected for the 21st century in comparison to recent climatic conditions [5,6,7]. Therefore, increased heat-related mortality is expected in future [2,8]. On the other hand, studies have documented declining impact of heat waves in developed countries due to behavioural and technological adaptation [9,10,11], including more effective public warning and response [12]. Therefore, uncertainty in the projection of temperature-related mortality is still large [2,8,13,14].
Young children, the elderly and persons in poor physical and medical condition [15,16], and especially those with cardiovascular and respiratory diseases, are most affected by heat stress [1,17,18]. In the Czech Republic, excess all-cause and cardiovascular mortality due to hot conditions have been observed in several studies that considered the population as a whole [19,20]. The elderly, women and people with chronic cardiovascular diseases (CVDs) were found to be the populations groups with the largest heat-related mortality [21,22,23].
Many studies have investigated the effects of heat waves on urban populations, which are those most affected by heat-related mortality [24,25,26]. a higher percentage of elderly, socially deprived and isolated populations (people living alone, immigrants, ethnic minorities) have been identified as factors increasing heat-related mortality [27,28,29]. In addition to socioeconomic factors, such land cover characteristics as the amount of impervious surface, green space, and unvegetated areas in the neighbourhood have been associated with increased heat vulnerability [27,30,31,32] and heat-related mortality [29,33,34]. Heat stress has nevertheless been found to be a significant risk also for the rural population [27,35,36,37]. Analyses in the Czech Republic have revealed comparably short-term response of cardiovascular mortality in both urban and rural population regardless of what temperature measure was used [38,39]. Spatial patterns in temperature-health relationships taking socioeconomic and environmental factors into account have not yet been investigated in the Czech Republic, however.
During heat waves, excess mortality often occurs for several days which are then followed by a decrease of mortality rates to below expected levels (so-called “mortality displacement” or “harvesting”; [40,41,42]). Although magnitudes of the short-term mortality displacement have been observed to vary considerably among populations of different cities [1,43,44,45,46], studies assessing spatial variability in heat-related mortality according to urban vs. rural and socioeconomic differences in a specific geographic area do not usually address the harvesting issue (including [38,39]).
In this study, we investigate spatial as well as temporal patterns for the effects of high temperature extremes on cardiovascular mortality in the Czech Republic at a district level. Effects of demographic, socioeconomic, and physical–environmental factors on excess cardiovascular mortality in 76 administrative districts are identified. Groups of districts with similar characteristics are compared in order to provide a more general picture and to evaluate different temporal patterns of excess mortality after hot spell occurrence in (i) urban areas vs. predominantly rural areas, and (ii) regions with different overall socioeconomic level.

2. Data and Methods

2.1. Study Area

The Czech Republic is a Central European country with a temperate climate and relatively varied landscape (Figure 1). Its approximately 10.2 million total inhabitants (as of 2001, with only minor changes during 1994–2009) are distributed into 77 districts (Nomenclature of Units for Territorial Statistics level 4 (NUTS 4) according to EUROSTAT [47]), ranging from 42,000 inhabitants in the Jeseník District to about 1.2 million inhabitants in the capital city of Prague. While our mortality data start in 1994, the Jeseník District (CZ0711—Czech Statistical Office coding [48]) was established in 1996 by its separation from the Šumperk District. In order to retain homogeneity in mortality data, we considered Jeseník and Šumperk as a single district in our study (CZ0715). Thus, we analysed 76 districts.

2.2. Meteorological Data

We considered average daily temperature in each district as a proxy variable for ambient thermal conditions. For its calculation, we used a high-resolution regular temperature grid covering the whole of the Czech Republic (see Figure 1). The GriSt data set with resolution 25 × 25 km is based on interpolated mean daily temperature data from irregularly spaced meteorological stations operated by the Czech Hydrometeorological Institute as described in Kyselý and Plavcová [49]. Such data were not available for other meteorological variables needed as inputs for calculating biometeorological indices (see Błażejczyk et al. [50]). Every grid box provides information about mean daily temperature from 1 January 1994 to 31 December 2009 and average altitude. Total population in each grid box in 2001 was calculated from the Gridded Population of the World (v3) [51]. Finally, the average daily temperature in each district was calculated from grids falling into the district area weighted by population within these grids (Figure 2). Average altitude of grids falling into each district was calculated as well. Zonal statistics tools in QGIS Desktop 2.4.0 were used for transformation of the grid data to districts.

2.3. Mortality Data

We analysed daily mortality from cardiovascular diseases (CVDs; defined as codes I00–I99 according to the International Statistical Classification of Diseases, 10th Revision [ICD-10]) in each district during 1994–2009. Mortality from CVDs comprises more than 50% of total mortality in the Czech Republic with a total record of 930,659 deaths occurring during the study period [52]. As the population count and structure vary distinctly between districts, daily mortality rates were calculated by the direct standardization procedure, using the mid-year population of each district in the Czech Republic and the standard WHO European population as the standard [53]. Daily baseline mortality for individual districts was determined using a single location-stratified generalized additive model (“mgcv” package in R (version 2.15.2)) [54]. A spline function with seven degrees of freedom (df) per year (total df = 112) was used according to Bhaskaran et al. [55], in order to take account for long-term trend and seasonality. Additional components taking account for the day of week and the district were defined by categorical variables. Mean relative deviations of mortality from the baseline (excess mortality) on days with average temperature above the 90th percentile (hot days) in summer seasons (June–August) 1994–2009 were calculated for all districts (DevCVD). The 90th percentile instead of the 95th percentile was chosen because it yields larger sample size, which is beneficial for the analysis at district level. Statistical significance of mean relative deviations was evaluated by comparison with the 90% confidence bounds around the zero line, estimated from the 5% and 95% quantiles of a distribution calculated by the Monte Carlo method (cf. Plavcová and Kyselý [56]). Spatial autocorrelation of mortality deviations on all summer days as well as on hot days was tested by Global Moran’s I criteria [57] in ArcGIS 10 software.

2.4. Analysis of Spatial Patterns of Heat-Related Mortality

We used demographic and socioeconomic characteristics from the Census 2001 database (provided by the Czech Statistical Office), land cover characteristics from the CORINE land cover 2000 database (the Czech Environmental Information Agency), and mean summer temperature (summer T) and mean altitude calculated from the GriSt data set in order to identify effects of demographic, socioeconomic, and physical–environmental factors on spatial differences in heat-related mortality. The selection of demographic and socioeconomic variables was based on a literature review [30,58,59,60,61] and adapted with respect to the limitations of the Census database. We applied a modified population density variable according to the OECD’s international terminology [62] (hereafter termed OECD criterion) which defines a level of urbanization by percentage of inhabitants living in municipalities with population density less than 150 inhabitants per 1 km2. An index of socioeconomic status (SES) for each district was calculated as a sum of z-scores for 3 factors of social deprivation related to elevated mortality in the Czech Republic [63]. Percentages of unemployed population (% unemployed), people without secondary school diploma (% low education), and single-person households (% singles) were included into the index. The final index is a dimensionless variable with a range between 2.66 (highest SES) and −5.88 (lowest SES) across the districts in the Czech Republic.
Land cover structure of each district was calculated using zonal statistics tools in QGIS Desktop 2.4.0. Percentage of impervious surface (% impervious) representing residential, industrial, commercial and transport areas was chosen similarly to other studies as a land cover variable enhancing heat stress [27,28,29,36,37,61,64]. The specific list of all variables used is presented in Table 1. Spearman’s correlation coefficients (and their p-values) of relationships between mean relative mortality deviations on hot days and independent variables were calculated using the “rcorr” function in the R package “Hmisc” (version 3.14-3) [65]. To check for multicollinearity between explanatory variables, we used variance inflation factors (VIFs) calculated in the ordinary least squares (OLS) regression in ArcGIS. If two or more explanatory variables had a large VIF (>7.5), the variable leading to the OLS regression model with lower percentage of variance explained was excluded from further analysis. After reduction of variables, we performed a bidirectional stepwise regression procedure testing all remaining explanatory variables, in order to identify the most significant variables affecting excess CVD mortality due to heat. The “stepAIC” function in the “MASS” package in R (version 2.15.2) adds or drops variables to or from the model repeatedly, fits those models, and computes a table of the changes in fit according to the Akaike information criterion (AIC) and amount of variance explained [66], i.e., the model with the lowest AIC and all significant parameters (as determined by an F test, p < 0.05) was selected. Spatial autocorrelation of the model residuals was evaluated with Global Moran’s I criteria. Ultimately, linear relationships of CVD mortality deviations with independent variables representing socioeconomic status (SES), population density (OECD criterion), and environmental effects (summer T) were examined individually for districts with low (<−0.50 StdDev), intermediate, and high (>0.50 StdDev) SES.

2.5. Analysis of Temporal Patterns of Heat-Related Mortality in Groups of Districts

Groups of districts with similar characteristics were defined in order to obtain greater population sample sizes and provide a more general picture of results on the district level. The districts with high SES and low SES were split into urban (OECD < 25%) and rural (OECD > 37.5%) districts according to the OECD criterion of the population density calculated for the Czech Republic by Blatecká [67]. Mortality rates and temperature data of all districts in each of the four groups (high SES-Urban, high SES-Rural, low SES-Urban, and low SES-Rural) were averaged and baseline mortality for each group using the same GAM formula as for individual districts was again calculated. The four groups of districts are displayed in Figure 3 and their characteristics are presented in Table 2.
Hot spells were defined as at least two consecutive days with average temperature above the 90th percentile of its distribution in summer seasons (June–August) 1994–2009 in a given group of districts. A percentile definition was chosen because it leads to similar numbers of days in hot spells in each group of districts. Relative deviations of CVD mortality from the baseline on days D − 2 (2 days before the beginning of a hot spell) up to D + 14 were averaged over all hot spells in order to assess lagged patterns of mortality deviations. Statistical significance of mean relative deviations was evaluated by comparison with the 90% and 95% confidence bounds around the zero line, estimated from the 2.5%, 5%, 95% and 97.5% quantiles of a distribution calculated by the Monte Carlo method (cf. Plavcová and Kyselý [56]). As the effect of a heat event on significantly elevated mortality usually occurs up to 3 days after the event [1,19,42], cumulative excess CVD mortality on days D + 0 to D + 3 and D + 4 to D + 14 after a hot spell’s onset was compared in order to estimate magnitude of the mortality displacement effect.

3. Results

3.1. Heat-Related Mortality in the Czech Republic at District Level

About 148 hot days were defined in each district during 1994–2009. The highest and statistically significant (p < 0.1) heat-related mortality (DevCVD) was generally experienced in the northwest and southeast Czech Republic and in large municipalities (Prague-code CZ0110, Pilsen-CZ0323, Brno-CZ0622, Ostrava-CZ0816; Figure 4). These areas are mostly located in the lowest and warmest parts of the Czech Republic (cf. Figure 1 and Figure 2). However, the test for spatial autocorrelation (Global Moran’s I) did not show any significant spatial clustering of mean mortality deviations on all summer days as well as on hot days.

3.2. Explaining Spatial Patterns of Heat-Related Mortality

Spearman’s correlation analysis revealed generally very weak and insignificant correlation of relative CVD mortality deviations (DevCVD) with most of the socioeconomic factors (Table 3). On the contrary, variables representing effects of physical environment and urbanization level (altitude, summer T, and % impervious) correlated significantly (p < 0.05) with mortality anomalies on hot days (r = 0.37, −0.37, and 0.34 respectively), although correlation coefficients for the significant variables were still rather weak. Variables correlating significantly with mortality deviations correlated also significantly with one another. In particular, altitude and summer T showed strong negative correlation (r = −0.92). These two variables were also associated with large VIFs (>7.5) which indicates their strong collinearity. Since summer T led to higher percentage of variance explained, altitude was excluded from the regression analysis.
The stepwise regression procedure identified district’s summer T as the most significant factor positively linked to excess mortality (Table 4). The increased population density (decreased OECD criterion) was chosen by the model as the second variable significantly associated with increased mortality. Global Moran’s I criteria close to 0 indicated no spatial autocorrelation of residuals from the stepwise regression model.
In the analyses for all districts considered together, regression procedure did not reveal any significant effect of socioeconomic status (SES) on heat-related mortality (DevCVD). However, regression analysis between DevCVD and selected explanatory variables showed significant changes in the relationships (Table 5, Figure 5) when performed individually for districts with low, intermediate and high SES. In high SES districts, the decreased OECD criterion (i.e., increased percentage of urban population) was the only variable significantly (p < 0.1) linked to increase in DevCVD. In the largest group of districts with the intermediate SES, only increased summer T was significantly (p < 0.05) associated with increased DevCVD. But in districts with low SES, increased DevCVD was significantly (p < 0.05) linked to increased summer T, decreased OECD criterion as well as decreased SES (considered as explanatory variable). Decreased SES was also significantly related to a decrease of the OECD criterion in the same category.
Results obtained from the analysis at the district level gave us general information about spatial distribution of heat-related excess mortality due to CVDs. Although in most districts the mean mortality deviations were positive and significantly different from zero (Figure 4), numbers of cases especially in rural districts were rather low and baseline estimates in such districts were affected by overdispersion. Therefore, the districts were divided into groups of districts according to the SES index and the OECD criterion in order to compare effects of high temperature in larger population samples.

3.3. Explaining Lagged Patterns of Heat-Related Mortality in Groups of Districts

About 35 hot spells were defined in each group of districts during 1994–2009. In all groups, significant excess mortality occurred between days D + 0 and D + 3 after a hot spell’s onset followed by a decline in mortality deviations (Figure 6). While in rural districts the lagged mortality deviations (D + 4 and later) were predominantly negative, in urban districts they were rather positive during the whole two-week period after a hot spell’s onset. When cumulative mortality deviations on days D + 0 to D + 3 and D + 4 to D + 14 after a hot spell’s onset were compared (Figure 7), the immediate cumulative effect (ƩD + 0…D + 3) was largest in the most densely populated and the warmest group of districts with high socioeconomic status (high SES-Urban, 58.2% of daily mortality), while it was lowest in the least populated group with low socioeconomic status (low SES-Rural, 20.4%). When taking into account the lagged effects (ƩD + 4…D + 14), however, we observed negative relative mortality deviations in rural groups of districts which results in total heat-effects during the two-week period close to zero (8.9% and −3.8% in high SES-Rural and low SES-Rural, respectively). On the contrary, the cumulative lagged effects were positive in urban districts. In the urban district with low SES (low SES-Urban), the lagged effect was comparable with the immediate one (26.3%), while in the high SES-Urban group the additional increase in excess mortality was rather small (8.4%).

4. Discussion

We analysed spatial and temporal patterns of links between daily CVD mortality and mean daily temperature during summer (June–August), 1994–2009 in the Czech Republic. We further evaluated these patterns with respect to demographic, socioeconomic, and physical-environmental differences. The largest relative mortality deviations on days that exceeded the 90th percentile of summer temperature distribution (hot days) were generally observed in the warmest and most densely populated regions of the Czech Republic.

4.1. Spatial Patterns of Heat-Related Mortality

When all districts in the Czech Republic were considered together, characteristics of physical environment (mean summer temperature, altitude) and urbanization level (population density, percentage of impervious surface) had a significant influence on heat-related CVD mortality, while socioeconomic status did not show any significant effect. On the other hand, a comparison of districts with low and high index of socioeconomic status (SES) revealed that below a certain threshold, SES has a relevant influence on excess mortality.
Although previous studies have demonstrated differences in the heat–health relationship between different cities and different population groups, there are only few studies including also rural areas and considering socioeconomic and land cover factors affecting spatial differences in heat-related mortality. In the U.S., Sheridan and Dolney [35] found a comparable effect of heat on mortality in urban and rural counties in Ohio. In Massachusetts [60] and Georgia [36], demographic (percentage of elderly, percentage of African-Americans) and socioeconomic (social isolation, poverty) factors were more important than was the level of urbanization (as represented by population density and proportion of the impervious surface). Recently, Kovach et al. [37] in North Carolina found greater rates of heat-related illnesses in rural than urban locations, and a greater proportion of workers in agriculture with less access to air conditioning has been mentioned among the reasons of high heat-vulnerability in rural populations [35,37]. On the other hand, decreased vegetation, living in poverty, and low education were associated with increases in heat-related illnesses in urban locations [37]. However, excess morbidity and mortality due to heat in U.S. cities is reduced by the protective effect of air conditioning [35,68]. Air conditioning is much less widespread in Europe than in American households [69], and population in the most densely populated areas is usually most affected by heat in Europe [24,26]. Moreover, in contrast to the U.S. and many western-European countries, social differences between urban, suburban and rural populations in the Czech Republic are much smaller, as the country ranks as having among the lowest income inequality of OECD countries [70], and the lowest estimated percentage of population severely deprived with respect to health, education and/or standard of living in the European Union (the Human Poverty Index, HPI; [71]). In addition, none of the aforementioned studies considered the effects of physical–geographic conditions such as local climate and topography. Strong collinearity between districts’ average temperature and altitude found in our results suggests that highest heat-related mortality is generally observed in districts with low altitude, relatively hotter local climate and higher level of urbanization. Nevertheless, the significant relationship between decreased socioeconomic status and increased heat-related mortality in the most deprived district group highlights the role of social deprivation in vulnerability to heat across the population of the Czech Republic.

4.2. Lagged Patterns of Heat-Related Mortality

Another limitation of many studies assessing spatial variability in heat-related mortality within specific geographic areas is that they do not take into account the mortality displacement (harvesting effect) issue. Analysis of the four groups of districts with larger sample sizes allowed us to study lagged patterns of mortality deviations with more statistical power than for individual districts. While in urban districts the lagged sums of mortality deviations were still positive, we observed large decreases in cumulative relative mortality deviations in rural districts. A decline in mortality a few weeks after a heat wave due to mortality displacement has been documented in many studies [42,43,44,45,46]. Although estimates of the mortality displacement magnitude vary based on the methodology used and between different populations [1,44], our findings suggest that, at least in the Czech Republic, the effect of mortality displacement after hot spells is small and insignificant in urban districts, while most of the excess deaths in rural districts may be attributed to harvesting.
The highest short-term effect of hot spells in the high SES-Urban group suggests that high socioeconomic status was a risk factor of increased excess CVD mortality in the stratified analysis. However, the regression analysis at the district level showed SES as a relevant factor only in the low SES districts. Therefore, highest population count and density, highest average summer temperature due to low altitude, and generally worst thermal conditions in comparison to other groups of districts due to high proportion of artificial surface are probably the main reasons for the largest short-term effect of hot spells on excess CVD mortality in the high SES-Urban group. Many studies have documented strongest heat effects on mortality in the most densely populated areas [24,25,26] due to urban heat island effects [72,73], although this may be altered by large-scale factors such as altitude and topography [74] as well as such fine-scale factors as proportions of impervious, green or water surfaces [34,75,76,77].
Short-term heat effect was reduced in the urban group of districts with low socioeconomic status (low SES-Urban) in comparison to high SES-Urban. However, when the cumulative relative mortality deviations during the two-week period after a hot spell’s onset were considered, excess CVD mortality in low SES-Urban was comparable with that in the high SES-Urban group (62% and 67%, respectively). This may be caused not only by relatively high levels of urbanization but also by social deprivation of populations in low SES-Urban districts, inasmuch as regression analysis at the district level had shown a significant relationship between socioeconomic status and heat-related CVD mortality for districts with low socioeconomic status. Dzúrová [78] had demonstrated significant relationships between environmental deprivation, social deprivation, and poor health condition within the population of the Czech Republic. Long-term exposure to high concentrations of particulate matter (PM) has been associated with development of chronic diseases of the veins and heart [79,80,81]. People with chronic CVDs have been most affected by excess mortality during hot spells in the Czech Republic [22,38]. Highest concentrations of PM10 and O3 have been associated with oppressive weather in Prague, although they were not identified as significant factors linked to excess CVD mortality by a stepwise regression model within the oppressive air masses [82]. Since heavy industry and coal mining are typical economic activities in districts of the low SES-Urban group, poor health condition due to long-term exposure to high air pollution concentrations [83] and in general highest environmental deprivation within the Czech Republic in these districts may have a significant effect on lagged excess mortality due to CVDs during hot spells. However, inasmuch as air pollution characteristics are not spatially representative because they vary considerably across space, and causal relationships among temperature, air pollution, and human health are associated with large uncertainty [84], we did not consider the effect of air pollution in this study.

4.3. Limitations of the Study

A usual issue of studies assessing spatial variability of weather-mortality relationships is to find a trade-off between geographic specificity and sample size when selecting a unit of analysis [60]. The length of time period (number of days) and number of events (e.g., deaths) per day are factors determining statistical power and precision of time series modelling [55]. For example, Hondula and Barnett [85] employed a hierarchical Bayesian model treating each unit as not completely independent from every other. We applied a single generalized additive model stratifying the whole dataset to individual districts by a categorical variable, which allowed us to model baseline mortality consistently in all districts and helped to deal with small population counts. Mortality deviations in individual districts are still subject to great uncertainty, as numbers of deaths in districts of the Czech Republic show high variability and only a small percentage of total variance can be explained by meteorological factors [82]. Therefore larger regions should be considered for heat-related mortality assessment and biometeorological forecast (e.g., [86]) in the Czech Republic rather than individual districts.
Our findings also highlight the importance of relating heat risk factors to appropriate health variables when evaluating vulnerable populations and regions. This is an issue in particular for studies mapping heat vulnerability without validation on health outcomes [30,58,64]. As demonstrated by Reid et al. [31], a heat vulnerability index calculated from socioeconomic and environmental risk factors can be consistently associated with poor health on both normal and hot days. In our study, low socioeconomic status (SES) was strongly correlated with high mortality (not shown), while the association with excess mortality on hot days was weak. Excess deaths during a heat wave can be considered as an extreme outcome of heat stress. Previous studies have revealed no significant effect of high temperatures on hospitalizations due to CVDs [23,38,82,87]. Nevertheless, health outcomes such as emergency calls [27,88] and hospital admissions for respiratory and heat-related diseases [31,89] might be significantly related to high temperatures and more affected by socioeconomic differences than cardiovascular mortality.
Another limitation of this study is that the associations observed for data aggregated at district level may not necessarily represent associations existing at an individual level—the so-called ecological fallacy issue [90]. Similarly, the associations between mortality and explanatory variables observed for the whole group of districts may vary spatially, as demonstrated by the stratified analysis for districts with low, intermediate and high SES. However, as we did not find any significant spatial autocorrelation in the residuals from baseline mortality nor in residuals from the step-wise regression procedure, we did not use any spatially varying coefficient model.
Finally, the study may potentially be subject of the modifiable areal unit problem as all variables were aggregated to district spatial units. Repeating the study at different analysis level may produce variations in the results [91]. However, the results for the individual districts and for aggregated regions supported each other.
Another issue of this study is the use of gridded temperature data. As the spatial distribution of the grid is 25 × 25 km, it does not fully represent fine-scale thermal and environmental conditions in individual municipalities. This is particularly true in the low SES-Urban districts which are located in areas with large differences in topography and thermal conditions ranging from relatively large cities situated in a warm climatic region to surrounding mountain areas with cooler climate [92]. In order to take this issue into account, temperature data was weighted by population counts in every grid box when mean daily temperature in districts was calculated and percentiles were used instead of absolute temperature thresholds for the definition of hot days. Nevertheless, the effect of average summer temperature and average altitude on mortality deviations on hot days may be partly underestimated in the correlation and regression analyses. On the other hand, the gridded temperature data, representing spatial means rather than point values [49], allow capturing temperature conditions in individual districts better than raw station data or their averages.

5. Conclusions

In this study we analysed spatial patterns of links between high temperature extremes and daily cardiovascular mortality at a district level. We considered differences between districts’ demographics, socioeconomic status, as well as physical-environmental conditions.
The largest relative mortality deviations on hot days were generally observed in the warmest regions of the Czech Republic that involved most urbanized areas. We revealed significant effects of physical environment (altitude, mean summer temperature, artificial surface) and urbanization (population density) on excess cardiovascular mortality on hot days when all districts in the Czech Republic were considered together. The level of socioeconomic status played a relevant role only in the most deprived populations.
Analysis of mortality deviation patterns for larger groups of districts during two weeks after a hot spell’s onset revealed no mortality displacement in large municipalities, while most of the excess deaths in rural districts may be attributed to harvesting. Highest population count and density, highest average temperature due to low altitude, and generally worst thermal conditions due to high proportion of artificial surface were the factors associated with largest excess cardiovascular mortality due to hot spells. In addition, long-term exposure to highest environmental and socioeconomic deprivation within the Czech Republic in urban districts with low socioeconomic status may have a significant effect on lagged excess mortality associated with hot spells due to the higher percentage of people with chronic cardiovascular disease. Inasmuch as coal mining and associated heavy industry (energetics, metallurgy) are typical economic activities in these districts, their populations are potentially at risk of increasing social deprivation in future due to a decline of mining and industrial transformation. The significant relationship between decreased socioeconomic status and increased heat-related mortality within the most deprived district group highlights the role of socioeconomic status in adaptability to heat across the population.
We identified districts and regions in the Czech Republic whose populations are most at risk of excess mortality due to heat stress and might be at increasing risk in the future. Our results are potentially useful for better targeting biometeorological forecasts and warnings, although larger regions rather than individual districts need to be considered to achieve a reasonable statistical power. Such other health outcomes as emergency calls and hospital admissions for respiratory and heat-related diseases should be considered in follow-up research regarding the effects of heat stress on public health in the Czech Republic.

Acknowledgments

The study was supported by the Czech Science Foundation under project P209/11/1985. Thanks are due to Bohumil Kříž and Jan Kynčl, National Institute of Public Health, for preparing epidemiological data sets as well as for useful discussions concerning their interpretation, and Ladislav Kážmér, Institute of Sociology CAS, for useful discussion concerning the Census database and its interpretation. Special thanks go to the Erasmus+ programme that enabled Aleš Urban a traineeship at Humboldt University in Berlin.

Author Contributions

Jan Kyselý defined the initial focus of this study. Aleš Urban together with Katrin Burkart, Jan Kyselý, Christian Schuster, and Tobia Lakes carried out the final research design. Hana Hanzlíková and Aleš Urban were responsible for gathering, cleaning and standardization of the mortality data. Petr Štěpánek gathered the temperature data and prepared the gridded dataset. Aleš Urban was responsible for gathering, cleaning and analysis of the socioeconomic and environmental data and their transformation to districts. Aleš Urban, Katrin Burkart and Jan Kyselý contributed to defining the time series analysis approach. Eva Plavcová performed nonparametric tests of mortality deviations. Aleš Urban performed most analyses, prepared initial interpretation of results and drafted the first manuscript. Katrin Burkart, Jan Kyselý, Christian Schuster, and Eva Plavcová provided critical revisions of the manuscript. All authors contributed to the final interpretation of results, and read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Basu, R. High ambient temperature and mortality: A review of epidemiologic studies from 2001 to 2008. Environ. Health 2009. [Google Scholar] [CrossRef] [PubMed]
  2. Gosling, S.N.; Lowe, J.A.; McGregor, G.R.; Pelling, M.; Malamud, B.D. Associations between elevated atmospheric temperature and human mortality: A critical review of the literature. Clim. Chang. 2009, 92, 299–341. [Google Scholar] [CrossRef]
  3. Hajat, S.; Kosatky, T. Heat-related mortality: A review and exploration of heterogeneity. J. Epidemiol. Community Health 2010, 64, 753–760. [Google Scholar] [CrossRef] [PubMed]
  4. Burkart, K.; Khan, M.M.H.; Schneider, A.; Breitner, S.; Langner, M.; Krämer, A.; Endlicher, W. The effects of season and meteorology on human mortality in tropical climates: A systematic review. Trans. R. Soc. Trop. Med. Hyg. 2014, 108, 393–401. [Google Scholar] [CrossRef] [PubMed]
  5. Gosling, S.N.; McGregor, G.R.; Lowe, J.A. Climate change and heat-related mortality in six cities Part 2: Climate model evaluation and projected impacts from changes in the mean and variability of temperature with climate change. Int. J. Biometeorol. 2009, 53, 31–51. [Google Scholar] [CrossRef] [PubMed]
  6. Ballester, J.; Rodó, X.; Giorgi, F. Future changes in Central Europe heat waves expected to mostly follow summer mean warming. Clim. Dyn. 2009, 35, 1191–1205. [Google Scholar] [CrossRef]
  7. IPCC. Summary for Policymakers. In Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Field, C.B., Barros, V.R., Dokken, D.J., Mach, K.J., Mastrandrea, M.D., Bilir, T.E., Chatterjee, M., Ebi, K.L., Estrada, Y.O., Genova, R.C., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014; pp. 1–32. [Google Scholar]
  8. Huang, C.; Barnett, A.G.; Wang, X.; Vaneckova, P.; Fitzgerald, G.; Tong, S. Projecting future heat-related mortality under climate change scenarios: A systematic review. Environ. Health Perspect. 2011, 119, 1681–1690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Matzarakis, A.; Muthers, S.; Koch, E. Human biometeorological evaluation of heat-related mortality in Vienna. Theor. Appl. Climatol. 2010, 105, 1–10. [Google Scholar] [CrossRef]
  10. Kyselý, J.; Plavcová, E. Declining impacts of hot spells on mortality in the Czech Republic, 1986–2009: Adaptation to climate change? Clim. Chang. 2011, 113, 437–453. [Google Scholar] [CrossRef]
  11. Bobb, J.F.; Peng, R.D.; Bell, M.L.; Dominici, F. Heat-related mortality and adaptation to heat in the United States. Environ. Health Perspect. 2014, 122, 811–816. [Google Scholar] [CrossRef] [PubMed]
  12. Toloo, G.S.; Fitzgerald, G.; Aitken, P.; Verrall, K.; Tong, S. Are heat warning systems effective? Environ. Health 2013. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Boeckmann, M.; Rohn, I. Is planned adaptation to heat reducing heat-related mortality and illness? A systematic review. BMC Public Health 2014. [Google Scholar] [CrossRef] [PubMed]
  14. Gosling, S.N.; McGregor, G.R.; Lowe, J.A. The benefits of quantifying climate model uncertainty in climate change impacts assessment: An example with heat-related mortality change estimates. Clim. Chang. 2012, 112, 217–231. [Google Scholar] [CrossRef]
  15. Mercer, J.B. Cold—An underrated risk factor for health. Environ. Res. 2002, 12, 8–13. [Google Scholar] [CrossRef]
  16. Kenney, W.L.; Munce, T.A. Invited review: Aging and human temperature regulation. J. Appl. Physiol. 2003, 95, 2598–2603. [Google Scholar] [CrossRef] [PubMed]
  17. McMichael, A.J.; Woodruff, R.E.; Hales, S. Climate change and human health: Present and future risks. Lancet 2006, 367, 859–869. [Google Scholar] [CrossRef]
  18. Cheng, X.; Su, H. Effects of climatic temperature stress on cardiovascular diseases. Eur. J. Intern. Med. 2010, 21, 164–167. [Google Scholar] [CrossRef] [PubMed]
  19. Kyselý, J.; Huth, R. Heat-related mortality in the Czech Republic examined through synoptic and “traditional” approaches. Clim. Res. 2004, 25, 265–274. [Google Scholar] [CrossRef]
  20. Kyselý, J.; Kříž, B. Decreased impacts of the 2003 heat waves on mortality in the Czech Republic: An improved response? Int. J. Biometeorol. 2008, 52, 733–745. [Google Scholar] [CrossRef] [PubMed]
  21. Kyselý, J.; Plavcová, E.; Davídkovová, H.; Kynčl, J. Comparison of hot and cold spell effects on cardiovascular mortality in individual population groups in the Czech Republic. Clim. Res. 2011, 49, 113–129. [Google Scholar] [CrossRef]
  22. Davídkovová, H.; Plavcová, E.; Kynčl, J.; Kyselý, J. Impacts of hot and cold spells differ for acute and chronic ischaemic heart diseases. BMC Public Health 2014. [Google Scholar] [CrossRef] [PubMed]
  23. Hanzlíková, H.; Plavcová, E.; Kynčl, J.; Kříž, B.; Kyselý, J. Contrasting patterns of hot spell effects on morbidity and mortality for cardiovascular diseases in the Czech Republic, 1994–2009. Int. J. Biometeorol. 2015. [Google Scholar] [CrossRef] [PubMed]
  24. Hajat, S.; Kovats, R.S.; Lachowycz, K. Heat-related and cold-related deaths in England and Wales: Who is at risk? Occup. Environ. Med. 2007, 64, 93–100. [Google Scholar] [CrossRef] [PubMed]
  25. Tan, J.; Zheng, Y.; Tang, X.; Guo, C.; Li, L.; Song, G.; Zhen, X.; Yuan, D.; Kalkstein, A.J.; Li, F. The urban heat island and its impact on heat waves and human health in Shanghai. Int. J. Biometeorol. 2010, 54, 75–84. [Google Scholar] [CrossRef] [PubMed]
  26. Gabriel, K.M.A.; Endlicher, W.R. Urban and rural mortality rates during heat waves in Berlin and Brandenburg, Germany. Environ. Pollut. 2011, 159, 2044–2050. [Google Scholar] [CrossRef] [PubMed]
  27. Uejio, C.K.; Wilhelmi, O.V.; Golden, J.S.; Mills, D.M.; Gulino, S.P.; Samenow, J.P. Intra-urban societal vulnerability to extreme heat: The role of heat exposure and the built environment, socioeconomics, and neighborhood stability. Health Place 2011, 17, 498–507. [Google Scholar] [CrossRef] [PubMed]
  28. Hondula, D.M.; Davis, R.E.; Leisten, M.J.; Saha, M.V.; Veazey, L.M.; Wegner, C.R. Fine-scale spatial variability of heat-related mortality in Philadelphia County, USA, from 1983 to 2008: A case-series analysis. Environ. Health 2012. [Google Scholar] [CrossRef] [PubMed]
  29. Klein Rosenthal, J.; Kinney, P.L.; Metzger, K.B. Intra-urban vulnerability to heat-related mortality in New York City, 1997–2006. Health Place 2014, 30, 45–60. [Google Scholar] [CrossRef] [PubMed]
  30. Reid, C.E.; O’Neill, M.S.; Gronlund, C.J.; Brines, S.J.; Brown, D.G.; Diez-Roux, A.V.; Schwartz, J. Mapping community determinants of heat vulnerability. Environ. Health Perspect. 2009, 117, 1730–1736. [Google Scholar] [PubMed]
  31. Reid, C.E.; Mann, J.K.; Alfasso, R.; English, P.B.; King, G.C.; Lincoln, R.A.; Margolis, H.G.; Rubado, D.J.; Sabato, J.E.; West, N.L.; et al. Evaluation of a heat vulnerability index on abnormally hot days: An environmental public health tracking study. Environ. Health Perspect. 2012, 120, 715–720. [Google Scholar] [CrossRef] [PubMed]
  32. Bao, J.; Li, X.; Yu, C. The construction and validation of the heat vulnerability index, a review. Int. J. Environ. Res. Public Health 2015, 12, 7220–7234. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, Y.; Dadvand, P.; Barrera-Gomez, J.; Sartini, C.; Mari-Dell’Olmo, M.; Borrell, C.; Medina-Ramon, M.; Sunyer, J.; Basagana, X. Differences on the effect of heat waves on mortality by sociodemographic and urban landscape characteristics. J. Epidemiol. Community Health 2013, 67, 519–525. [Google Scholar] [CrossRef] [PubMed]
  34. Burkart, K.; Meier, F.; Schneider, A.; Breitner, S.; Canário, P.; Joao Alcoforado, M.; Scherer, D.; Endlicher, W. Modification of heat-related mortality in an elderly urban population by vegetation (Urban Green) and proximity to water (Urban Blue): Evidence from Lisbon, Portugal. Environ. Health Perspect. 2015. [Google Scholar] [CrossRef]
  35. Sheridan, S.C.; Dolney, T.J. Heat, mortality, and level of urbanization: Measuring vulnerability across Ohio, USA. Clim. Res. 2003, 24, 255–265. [Google Scholar] [CrossRef]
  36. Maier, G.; Grundstein, A.; Jang, W.; Li, C.; Naeher, L.P.; Shepherd, M. Assessing the performance of a vulnerability index during oppressive heat across Georgia, United States. Weather Clim. Soc. 2014, 6, 253–263. [Google Scholar] [CrossRef]
  37. Kovach, M.M.; Konrad, C.E., II; Fuhrmann, C.M. Area-level risk factors for heat-related illness in rural and urban locations across North Carolina, USA. Appl. Geogr. 2015, 60, 175–183. [Google Scholar] [CrossRef]
  38. Urban, A.; Davídkovová, H.; Kyselý, J. Heat- and cold-stress effects on cardiovascular mortality and morbidity among urban and rural populations in the Czech Republic. Int. J. Biometeorol. 2014, 56, 1057–1068. [Google Scholar] [CrossRef] [PubMed]
  39. Urban, A.; Kyselý, J. Comparison of UTCI with other thermal indices in the assessment of heat and cold effects on cardiovascular mortality in the Czech Republic. Int. J. Environ. Res. Public Health 2014, 11, 952–967. [Google Scholar] [CrossRef] [PubMed]
  40. Huynen, M.M.; Martens, P.; Schram, D.; Weijenberg, M.P.; Kunst, A.E. The impact of heat waves and cold spells on mortality rates in the Dutch population. Environ. Health Perspect. 2001, 109, 463–470. [Google Scholar] [CrossRef] [PubMed]
  41. Laschewski, G.; Jendritzky, G. Effects of the thermal environment on human health: An investigation of 30 years of daily mortality data from SW Germany. Clim. Res. 2002, 21, 91–103. [Google Scholar] [CrossRef]
  42. Kyselý, J. Mortality and displaced mortality during heat waves in the Czech Republic. Int. J. Biometeorol. 2004, 49, 91–97. [Google Scholar] [CrossRef] [PubMed]
  43. Gosling, S.N.; McGregor, G.R.; Páldy, A. Climate change and heat-related mortality in six cities part 1: Model construction and validation. Int. J. Biometeorol. 2007, 51, 525–540. [Google Scholar] [CrossRef] [PubMed]
  44. Baccini, M.; Kosatsky, T.; Biggeri, A. Impact of summer heat on urban population mortality in Europe during the 1990s: An evaluation of years of life lost adjusted for harvesting. PLoS ONE 2013. [Google Scholar] [CrossRef] [PubMed]
  45. Saha, M.V.; Davis, R.E.; Hondula, D.M. Mortality displacement as a function of heat event strength in 7 US cities. Am. J. Epidemiol. 2014, 179, 467–474. [Google Scholar] [CrossRef] [PubMed]
  46. Zaninović, K.; Matzarakis, A. Impact of heat waves on mortality in Croatia. Int. J. Biometeorol. 2014, 58, 1135–1145. [Google Scholar] [CrossRef] [PubMed]
  47. EUROSTAT. Local Administrative Units. Available online: http://ec.europa.eu/eurostat/web/nuts/local-administrative-units (accessed on 20 October 2015).
  48. Czech Statistical Office—Nomenclature of Territorial Units for Statistics (CZ-NUTS)—2004. (in Czech Language). Available online: https://www.czso.cz/csu/czso/klasifikace_uzemnich_statistickych_jednotek_-cz_nuts- (accessed on 20 October 2015).
  49. Kyselý, J.; Plavcová, E. A critical remark on the applicability of E-OBS European gridded temperature data set for validating control climate simulations. J. Geophys. Res. 2010, 115, 1–14. [Google Scholar] [CrossRef]
  50. Błażejczyk, K.; Epstein, Y.; Jendritzky, G.; Staiger, H.; Tinz, B. Comparison of UTCI to selected thermal indices. Int. J. Biometeorol. 2012, 56, 515–535. [Google Scholar] [CrossRef] [PubMed]
  51. Gridded Population of the World, Version 3 (GPWv3) Data Collection. Center for International Earth Science Information Network (CIESIN), Colombia Univesity. 2005. Available online: http://sedac.ciesin.columbia.edu/data/collection/gpw-v3 (accessed on 20 September 2015).
  52. Davídkovová, H.; Kyselý, J.; Kříž, B.; Vojtíšek, P.; Bobák, M. Trends in cardiovascular mortality and hospitalisations, and potential contribution of in hospital case-fatality rates to changes in national mortality in the Czech Republic 1994–2009. Heart 2013, 99, 409–416. [Google Scholar] [CrossRef] [PubMed]
  53. Ahmad, O.; Boschi-Pinto, C.; Lopez, A.D.; Murray, C.J.L.; Lozano, R.; Inoue, M. Age Standardization of Rates: A New WHO Standard; World Health Organization: Geneva, Switzerland, 2001. [Google Scholar]
  54. Wood, S. Generalized Additive Models: An Introduction with R; Chapman & Hall/CRC: London, UK, 2006. [Google Scholar]
  55. Bhaskaran, K.; Gasparrini, A.; Hajat, S.; Smeeth, L.; Armstrong, B. Time series regression studies in environmental epidemiology. Int. J. Epidemiol. 2013, 42, 1187–1195. [Google Scholar] [CrossRef] [PubMed]
  56. Plavcová, E.; Kyselý, J. Relationships between sudden weather changes in summer and mortality in the Czech Republic, 1986–2005. Int. J. Biometeorol. 2010, 54, 539–551. [Google Scholar] [CrossRef] [PubMed]
  57. Getis, A.; Ord, J.K. The analysis of spatial association by use of distance statistics. Geogr. Anal. 1992, 24, 189–206. [Google Scholar] [CrossRef]
  58. Vescovi, L.; Rebetez, M.; Rong, F. Assessing public health risk due to extremely high temperature events: Climate and social parameters. Clim. Res. 2005, 30, 71–78. [Google Scholar] [CrossRef]
  59. Wichmann, J.; Andersen, Z.J.; Ketzel, M.; Ellermann, T.; Loft, S. Apparent temperature and cause-specific mortality in Copenhagen, Denmark: A case-crossover analysis. Int. J. Environ. Res. Public Health 2011, 8, 3712–3727. [Google Scholar] [CrossRef] [PubMed]
  60. Hattis, D.; Ogneva-Himmelberger, Y.; Ratick, S. The spatial variability of heat-related mortality in Massachusetts. Appl. Geogr. 2012, 33, 45–52. [Google Scholar] [CrossRef]
  61. Harlan, S.L.; Declet-barreto, J.H.; Stefanov, W.L.; Petitti, D.B. Neighborhood effects on heat deaths: Social and environmental predictors of vulnerability in Maricopa County, Arizona. Environ. Health Perspect. 2013, 121, 197–204. [Google Scholar] [PubMed]
  62. Spiezia, V. Measuring Regional Economies. Statistics Brief OECD, No. 6. 2003. Available online: http://www.oecd.org/dataoecd/2/15/15918996.pdf (accessed on 14 August 2013).
  63. Šplíchalová, A.; Šlachtová, H.; Tomášková, H.; Polaufová, P. Analysis of mortality in the Czech Republic using the newly constructed socio-economic deprivation index. Epidemiology 2007. [Google Scholar] [CrossRef]
  64. Johnson, D.P.; Stanforth, A.; Lulla, V.; Luber, G. Developing an applied extreme heat vulnerability index utilizing socioeconomic and environmental data. Appl. Geogr. 2012, 35, 23–31. [Google Scholar] [CrossRef]
  65. Harrell, F.E.; Hmisc, S. Function Library. 2004. Available online: http://biostat.mc.vanderbilt.edu/s/Hmisc (accessed on 8 October 2015).
  66. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
  67. Blatecká, K. The Basic Characteristics of Czech Rural Areas. Ph.D Thesis, Masaryk University, Brno, Czech Republic, 2006. Available online: http://is.muni.cz/th/137827/esf_m (accessed on 14 August 2013) . [Google Scholar]
  68. Medina-Ramón, M.; Schwartz, J. Temperature, Temperature Extremes, and Mortality: A study of acclimatization and effect modification in 50 United States cities. Occup. Environ. Med. 2007, 64, 827–834. [Google Scholar] [CrossRef] [PubMed]
  69. Henderson, G. Home air conditioning in Europe—How much energy would we use if we became more like American households? In Energy Savings: What Works & Who Delivers; Attali, S., Tillerson, K., Eds.; European Council for an Energy Efficient Economy: Stockholm, Sweden, 2005; pp. 541–550. [Google Scholar]
  70. OECD. Income Distribution Database (IDD): Gini, Poverty, Income, Methods and Concepts. Available online: http://www.oecd.org/social/income-distribution-database.htm (accessed on 20 September 2015).
  71. Bubbico, R.L.; Dijkstra, L. Regional Focus. Available online: http://ec.europa.eu/regional_policy/sources/docgener/focus/2011_02_hdev_hpov_indices.pdf (accessed on 20 September 2015).
  72. Arnfield, A.J. Two decades of urban climate research: A review of turbulence, exchanges of energy and water, and the urban heat island. Int. J. Climatol. 2003, 23, 1–26. [Google Scholar] [CrossRef]
  73. Fischer, E.M.; Oleson, K.W.; Lawrence, D.M. Contrasting urban and rural heat stress responses to climate change. Geophys. Res. Lett. 2012, 39, 1–8. [Google Scholar] [CrossRef]
  74. Alcoforado, M.-J.; Andrade, H. Nocturnal urban heat island in Lisbon (Portugal): Main features and modelling attempts. Theor. Appl. Climatol. 2006, 84, 151–159. [Google Scholar] [CrossRef]
  75. Unger, J. Intra-urban relationship between surface geometry and urban heat island: Review and new approach. Clim. Res. 2004, 27, 253–264. [Google Scholar] [CrossRef]
  76. Dütemeyer, D.; Barlag, A.B.; Kuttler, W.; Axt-Kittner, U. Measures against heat stress in the city of Gelsenkirchen, Germany. DIE ERDE 2013, 144, 181–201. [Google Scholar]
  77. Dugord, P.-A.; Lauf, S.; Schuster, C.; Kleinschmit, B. Land use patterns, temperature distribution, and potential heat stress risk—The case study Berlin, Germany. Comput. Environ. Urban Syst. 2014, 48, 86–98. [Google Scholar] [CrossRef]
  78. Dzúrová, D. Zdravotní Stav Obyvatel a Kvalita životního Prostředí: Regionální Souvislosti. Ph.D. Thesis, Faculty of Science, Charles University, Prague, Czech Republic, 1993; p. 136. [Google Scholar]
  79. Brook, R.D.; Franklin, B.; Cascio, W.; Hong, Y.; Howard, G.; Lipsett, M.; Luepker, R.; Mittleman, M.; Samet, J.; Sidney, C.; et al. Air pollution and cardiovascular disease: A statement for healthcare professionals from the expert panel on population and prevention science of the American Heart Association. Circulation 2004, 109, 2655–2671. [Google Scholar] [CrossRef] [PubMed]
  80. Chen, H.; Goldberg, M.S. The effects of outdoor air pollution on chronic illnesses. McGill J. Med. 2009, 12, 58–64. [Google Scholar] [PubMed]
  81. Kaufman, J.D. Does Air pollution accelerate progression of atherosclerosis? JAC 2010, 56, 1809–1811. [Google Scholar] [CrossRef] [PubMed]
  82. Urban, A.; Kyselý, J. Application of spatial synoptic classification in evaluating links between heat stress and cardiovascular mortality and morbidity in Prague, Czech Republic. Int. J. Biometeorol. 2015. [Google Scholar] [CrossRef] [PubMed]
  83. Czech Hydrometeorological Institute. Air Pollution in the Czech Republic. Available online: http://portal.chmi.cz/files/portal/docs/uoco/isko/tab_roc/tab_roc_EN.html (accessed on 20 October 2015).
  84. Buckley, J.P.; Samet, J.M.; Richardson, D.B. Commentary: Does air pollution confound studies of temperature? Epidemiology 2014, 25, 242–245. [Google Scholar] [CrossRef] [PubMed]
  85. Hondula, D.M.; Barnett, A.G. Heat-related morbidity in Brisbane, Australia: Spatial variation and area-level predictors. Environ. Health Perspect. 2014, 122, 831–836. [Google Scholar] [CrossRef] [PubMed]
  86. Novák, M. Use of the UTCI in the Czech Republic. Geogr. Pol. 2013, 86, 21–28. [Google Scholar] [CrossRef]
  87. Zacharias, S.; Koppe, C.; Mücke, H.-G. Influence of heat waves on ischemic heart diseases in Germany. Climate 2014, 2, 133–152. [Google Scholar] [CrossRef]
  88. Wolf, T.; McGregor, G.; Analitis, A. Performance assessment of a heat wave vulnerability index for greater London, United Kingdom. Weather Clim. Soc. 2014, 6, 32–46. [Google Scholar] [CrossRef]
  89. Scherber, K.; Langner, M.; Endlicher, W. Spatial analysis of hospital admissions for respiratory diseases during summer months in Berlin taking bioclimatic and socio-economic aspects into account. Geogr. Soc. Berl. 2014, 144, 217–237. [Google Scholar]
  90. Wakefield, J.; Salway, R. A statistical framework for ecological and statistical studies. J. R. Stat. Soc. 2001, 164, 119–137. [Google Scholar] [CrossRef]
  91. Openshaw, S.; Taylor, P.J. The modifiable areal unit problem. In Quantitative Geography: A British View; Wrigley, N., Bennett, R.J., Eds.; Routledge: Boston, MA, USA, 1981; pp. 60–69. [Google Scholar]
  92. Tolasz, R.; Tolasz, R.; Míková, T.; Valeriánová, A.; Voženílek, V. Climate Atlas of Czechia; Czech Hydrometeorological Institute: Prague, Czech Republic, 2007; p. 254. [Google Scholar]
Ijerph 13 00284 g001 1024
Figure 1. Topography of the Czech Republic (Digital Elevation Model provided by ARCDATA PRAHA) and distribution of regular temperature grid points (GriSt data set as described in [49]).
Figure 1. Topography of the Czech Republic (Digital Elevation Model provided by ARCDATA PRAHA) and distribution of regular temperature grid points (GriSt data set as described in [49]).
Ijerph 13 00284 g001
Ijerph 13 00284 g002 1024
Figure 2. Average summer temperature (June–August 1994–2009) in districts of the Czech Republic calculated from the GriSt data set. Districts are labelled according to the Czech Statistical Office coding [48].
Figure 2. Average summer temperature (June–August 1994–2009) in districts of the Czech Republic calculated from the GriSt data set. Districts are labelled according to the Czech Statistical Office coding [48].
Ijerph 13 00284 g002
Ijerph 13 00284 g003 1024
Figure 3. Four groups of districts identified according to high (>0.50 StdDev) and low (<−0.50 StdDev) index of socioeconomic status (SES) and urban (OECD criterion <25%) and rural (OECD criterion > 37.5%) population.
Figure 3. Four groups of districts identified according to high (>0.50 StdDev) and low (<−0.50 StdDev) index of socioeconomic status (SES) and urban (OECD criterion <25%) and rural (OECD criterion > 37.5%) population.
Ijerph 13 00284 g003
Ijerph 13 00284 g004 1024
Figure 4. Spatial distribution of mean relative mortality deviations due to CVDs on hot days (DevCVD). Districts are labelled according to the Czech Statistical Office coding [48].
Figure 4. Spatial distribution of mean relative mortality deviations due to CVDs on hot days (DevCVD). Districts are labelled according to the Czech Statistical Office coding [48].
Ijerph 13 00284 g004
Ijerph 13 00284 g005 1024
Figure 5. Scatterplots comparing mean relative excess mortality due to CVDs on hot days (DevCVD) with index of socioeconomic status (SES), percentage of rural population (OECD criterion), and mean summer temperature (Summer t) for districts with low (<−0.50 StdDev), intermediate (-), and high (>0.50 StdDev) socioeconomic status (SES). Coefficients of regression lines are reported in Table 5. See Table 1 for description of variables examined.
Figure 5. Scatterplots comparing mean relative excess mortality due to CVDs on hot days (DevCVD) with index of socioeconomic status (SES), percentage of rural population (OECD criterion), and mean summer temperature (Summer t) for districts with low (<−0.50 StdDev), intermediate (-), and high (>0.50 StdDev) socioeconomic status (SES). Coefficients of regression lines are reported in Table 5. See Table 1 for description of variables examined.
Ijerph 13 00284 g005
Ijerph 13 00284 g006 1024
Figure 6. Mean relative excess mortality due to CVDs on days D − 2 to D + 14 around a hot spell’s onset in the groups of districts. Confidence bounds around the zero line are indicated by dashed (90%) and solid (95%) lines, respectively.
Figure 6. Mean relative excess mortality due to CVDs on days D − 2 to D + 14 around a hot spell’s onset in the groups of districts. Confidence bounds around the zero line are indicated by dashed (90%) and solid (95%) lines, respectively.
Ijerph 13 00284 g006
Ijerph 13 00284 g007 1024
Figure 7. Cumulative excess mortality (relative to the daily baseline mortality) due to CVDs on days D + 0 to D + 3 (ƩD + 0…D + 3) and D + 4 to D + 14 (ƩD + 4…D + 14) after a hot spell’s onset in the groups of districts. * denotes mean cumulative excess mortality above the 95% quantile of a distribution calculated by the Monte Carlo method.
Figure 7. Cumulative excess mortality (relative to the daily baseline mortality) due to CVDs on days D + 0 to D + 3 (ƩD + 0…D + 3) and D + 4 to D + 14 (ƩD + 4…D + 14) after a hot spell’s onset in the groups of districts. * denotes mean cumulative excess mortality above the 95% quantile of a distribution calculated by the Monte Carlo method.
Ijerph 13 00284 g007
Table
Table 1. District characteristics tested in the analysis.
Table 1. District characteristics tested in the analysis.
VariableDescription
DevCVDMean relative CVD mortality deviations on hot days with average temperature above the 90th percentile in June–August 1994–2009
SESIndex of socioeconomic status; sum of z-scores for % low education, % unemployed and % singles
% elderly% of population older than 65 years
% low education% of population without secondary school diploma
% unemployed% of unemployed population
% singles% of single-person households
OECD% of inhabitants in municipalities with population density less than 150 inhabitants/km2
Summer T (°C)Mean summer (1 June–31 August 1994–2009) temperature in °C
Altitude (m a.s.l.)Average altitude in metres above sea level
% impervious% of impervious surface area (categories 1.1 and 1.2 in the CORINE land cover classification)
Table
Table 2. Total population and average demographic, socioeconomic, and physical–environmental characteristics (in 2001) of the four groups of districts presented in Figure 3.
Table 2. Total population and average demographic, socioeconomic, and physical–environmental characteristics (in 2001) of the four groups of districts presented in Figure 3.
CharacteristicHigh SES-UrbanHigh SES-RuralLow SES-RuralLow SES-Urban
Population1,871,0951,197,212694,1151,170,971
OECD5.9348.9550.8810.51
SES2.061.56−1.85−3.62
% elderly15.5613.9112.2011.94
% low education49.6765.0868.5965.73
% unemployed3.813.156.648.16
% singles33.8227.0129.2333.54
Summer T (°C)18.417.417.617.5
Altitude (m a.s.l)327459392415
% impervious29.824.454.4013.42
Table
Table 3. Spearman’s correlation matrix comparing the mean relative excess mortality due to CVDs on hot days (DevCVD) with independent variables used in the analysis (see Table 1 for description). Significant coefficients with p-value < 0.05 are in bold.
Table 3. Spearman’s correlation matrix comparing the mean relative excess mortality due to CVDs on hot days (DevCVD) with independent variables used in the analysis (see Table 1 for description). Significant coefficients with p-value < 0.05 are in bold.
Independent VariableDevCVDSES% Elderly% low Education% Unemployed% SinglesOECDSummer T (°C)Altitude (m a.s.l.)
SES0.011
% elderly0.140.611
% low education−0.20−0.52−0.471
% unemployed0.12−0.75−0.530.331
% singles0.13−0.38−0.04−0.240.091
OECD−0.220.110.070.47−0.24−0.561
Summer T (°C)0.370.180.38−0.190.10−0.08−0.191
Altitude (m a.s.l.)−0.37−0.05−0.310.17−0.13−0.090.26−0.921
% impervious 0.340.050.12−0.410.240.23−0.740.64−0.66
Table
Table 4. Independent variables (see Table 1) chosen by the stepwise regression model as significantly related to the mean relative excess mortality due to CVDs on hot days (DevCVD) in individual districts. Regression coefficients are reported, along with their associated p-values in parentheses.
Table 4. Independent variables (see Table 1) chosen by the stepwise regression model as significantly related to the mean relative excess mortality due to CVDs on hot days (DevCVD) in individual districts. Regression coefficients are reported, along with their associated p-values in parentheses.
Independent VariableDevCVD
SES----
% elderly----
% low education----
% unemployed----
% singles----
OECD−0.096 (0.043)
Summer T (°C)2.922 (0.002)
% impervious ----
R20.191
Table
Table 5. Linear regression models for mean relative excess mortality due to CVDs on hot days (DevCVD) against index of socioeconomic status (SES), percentage of rural population (OECD criterion), and mean summer temperature (Summer T) for districts with low, intermediate, and high socioeconomic status (SES class). Regression coefficients are reported, along with their associated p-values in parentheses. See Table 1 for description of variables examined. Significant regression coefficients with p-value < 0.05 are in bold.
Table 5. Linear regression models for mean relative excess mortality due to CVDs on hot days (DevCVD) against index of socioeconomic status (SES), percentage of rural population (OECD criterion), and mean summer temperature (Summer T) for districts with low, intermediate, and high socioeconomic status (SES class). Regression coefficients are reported, along with their associated p-values in parentheses. See Table 1 for description of variables examined. Significant regression coefficients with p-value < 0.05 are in bold.
SES ClassDevCVD vs. SESDevCVD vs. Summer T (°C)DevCVD vs. OECDSES vs. OECD
Low−3.099 (0.034)4.307 (0.025)−0.193 (0.044)0.041 (0.000)
Intermediate−3.399 (0.490)3.537 (0.008)−0.073 (0.407)0.008 (0.293)
High1.474 (0.842)1.234 (0.181)−0.164 (0.078)−0.008 (0.347)
R20.0860.1600.1070.878

Share and Cite

MDPI and ACS Style

Urban, A.; Burkart, K.; Kyselý, J.; Schuster, C.; Plavcová, E.; Hanzlíková, H.; Štěpánek, P.; Lakes, T. Spatial Patterns of Heat-Related Cardiovascular Mortality in the Czech Republic. Int. J. Environ. Res. Public Health 2016, 13, 284. https://doi.org/10.3390/ijerph13030284

AMA Style

Urban A, Burkart K, Kyselý J, Schuster C, Plavcová E, Hanzlíková H, Štěpánek P, Lakes T. Spatial Patterns of Heat-Related Cardiovascular Mortality in the Czech Republic. International Journal of Environmental Research and Public Health. 2016; 13(3):284. https://doi.org/10.3390/ijerph13030284

Chicago/Turabian Style

Urban, Aleš, Katrin Burkart, Jan Kyselý, Christian Schuster, Eva Plavcová, Hana Hanzlíková, Petr Štěpánek, and Tobia Lakes. 2016. "Spatial Patterns of Heat-Related Cardiovascular Mortality in the Czech Republic" International Journal of Environmental Research and Public Health 13, no. 3: 284. https://doi.org/10.3390/ijerph13030284

APA Style

Urban, A., Burkart, K., Kyselý, J., Schuster, C., Plavcová, E., Hanzlíková, H., Štěpánek, P., & Lakes, T. (2016). Spatial Patterns of Heat-Related Cardiovascular Mortality in the Czech Republic. International Journal of Environmental Research and Public Health, 13(3), 284. https://doi.org/10.3390/ijerph13030284

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

Article Metrics

Back to TopTop