Next Article in Journal
Feasibility and Reliability of a Physical Fitness Test Battery in Individuals with Down Syndrome
Next Article in Special Issue
An Occupational Heat–Health Warning System for Europe: The HEAT-SHIELD Platform
Previous Article in Journal
Cumulative Impact of Environmental Pollution and Population Vulnerability on Pediatric Asthma Hospitalizations: A Multilevel Analysis of CalEnviroScreen
Previous Article in Special Issue
Overview of Existing Heat-Health Warning Systems in Europe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heat Warnings in Switzerland: Reassessing the Choice of the Current Heat Stress Index

1
Federal Office of Meteorology and Climatology MeteoSwiss, 8058 Zurich Airport, Switzerland
2
Meteorology Group, Department of Applied Mathematics and Computer Sciences, University of Cantabria, 39005 Santander, Spain
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(15), 2684; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152684
Submission received: 17 June 2019 / Revised: 20 July 2019 / Accepted: 25 July 2019 / Published: 27 July 2019

Abstract

:
High temperatures lead to heat-related human stress and an increased mortality risk. To quantify heat discomfort and the relevant dangers, heat stress indices combine different meteorological variables such as temperature, relative humidity, radiation and wind speed. In this paper, a set of widely-used heat stress indices is analyzed and compared to the heat index currently used to issue official heat warnings in Switzerland, considering 28 Swiss weather stations for the years 1981–2017. We investigate how well warnings based on the heat index match warning days and warning periods that are calculated from alternative heat stress indices. The latter might allow for more flexibility in terms of specific warning demands and impact-based warnings. It is shown that the percentage of alternative warnings that match the official warnings varies among indices. Considering the heat index as reference, the simplified wet bulb globe temperature performs well and has some further advantages such as no lower bound and allowing for the calculation of climatological values. Yet, other indices (e.g., with higher dependencies on humidity) can have some added value, too. Thus, regardless of the performance in terms of matches, the optimal index to use strongly depends on the purpose of the warning.

1. Introduction

Scientific evidence from the analysis of long-term climate records and future climate projections suggests an increase in the frequency and duration of extreme temperature events during this century [1,2,3]. The IPCC Special Report on Extremes ranks heat waves among the most severe risks associated with climate change [4]. It is further widely acknowledged that higher temperatures lead to heat-related human stress and discomfort [5,6], which is likewise expected to increase under a warmer climate (e.g., [7]), also in Switzerland [8]. Robust and effective heat warnings are essential to increasing people’s awareness of the risk of extreme heat situations [9]. With multiple societal sectors being in need of reliable warnings, an effective warning system is supposed to be flexible in order to account for variable demographical (children and the elderly as especially vulnerable groups) and occupational demands (indoor and outdoor activities with different levels of physical labor, clothing etc.) [10]. Additionally, regional-specific climatic conditions and phenomena associated with different levels of acclimatization among individuals, need to be taken into consideration [10,11]. In cities, people are potentially more affected by high temperatures than those in the surrounding areas, especially during night-time, as urban fabrics absorb and store more energy than rural areas, i.e., the Urban Heat Island—UHI—effect (e.g., [12,13]). Diversifying heat warnings is thus of major relevance for being able to issue user-tailored, impact-based warnings, that can guarantee the highest effectiveness. Within some meteorological projects or services, these requirements have been identified and systematically implemented already. The European project HEAT-SHIELD, for instance, includes the development of a personalized heat warning system particularly for occupational settings [14]. The underlying heat stress index wet bulb globe temperature (wbgt), which is solely based on meteorological parameters, allows for ISO-conformal guidance (i.e., according to the International Organization for Standardization, ISO) to reduce health impacts on the work force and to prevent associated productivity losses for the considered industries [15,16]. The diversified warning system of the German Meteorological Service (Deutscher Wetterdienst) represents another example. It combines perceived temperature (“Gefühlte Temperatur”) with daily minimum temperature inside buildings [17] to account for the typically higher heat stress in urban settings (UHI effect). User-tailored warnings are issued for all German municipalities exceeding 100,000 inhabitants [18]. Conscious of the special vulnerability of elderly people, the warning system is additionally adapted for high-risk groups [18].
Not all heat stress indices are based on meteorological data exclusively (so-called direct thermal indices; [19,20]) such as the heat index (hi) used in the heat warning systems in Switzerland and the United States (see Section 2; [21,22]). Some of them are based on objective and subjective strain (so-called empirical thermal indices; [19,20]) such as the physiological strain index (PSI; [23]), and still others further include physiological parameters and are thus able to account for the interrelation between metabolic activities, clothing and environmental parameters (so-called rational thermal indices; [19,20,24]). The universal thermal climate index (utci), for instance, is based on the heat balance of the human body. More precisely, it uses a multi-node human heat balance model, to include the meteorologically-induced heat stress a human body experiences when trying to maintain a thermal equilibrium with the outdoor surrounding [25,26,27,28]. The human heat balance model is augmented with a sophisticated clothing model, defining the effective clothing insulation and vapor resistance values for each of the thermo-physiological model’s body segments over various climatic conditions [29]. Although rational thermal indices are often based on highly advanced physiological models of the body’s response, simpler, direct thermal indices have certain advantages: they are especially suitable for long-term studies such as historical analyses of the climate [30] or climate projections. Other advantages are their less intense computational costs and higher flexibility. For instance, purely meteorologically-based indices only have to be computed once for a given meteorological setting, and can in a second step be interpreted or modified in terms of the specific individual setting (clothing, activity, grade of acclimatization etc.), which allows for user-friendly applications.
Direct thermal indices are often the preferred variables used to issue heat warnings in European countries, including Switzerland [31]. We present a case study on the comparison of the current heat warning system in Switzerland to a (virtual) new system that makes use of different and potentially more appropriate heat stress indices, which allows for a more flexible use and diversification in terms of specific warning demands. This study can be considered as a first step towards the general desire for expanding the current warning system to a more impact-oriented and diverse system, which, by definition, cannot rely on one single threshold or index, but needs to consider a broader set of user needs. The paper is structured as follows: First, we introduce the heat warning system in Switzerland. Second, the considered heat stress indices and data are described and we explain the employed methods. Third, we compare heat warnings based on different indices with the official heat warnings and analyze their potential matches. The paper ends with the conclusions.

2. Heat Warning System in Switzerland

The current heat warning system in Switzerland is based on the hi, which is purely meteorological, i.e., a direct thermal index, and combines temperature and relative humidity. The hi equation, derived from the apparent temperature [21], is a refinement of a result obtained by multiple regression analyses [22]. If daily maxima hi exceeds the Swiss-wide threshold of 90°F for at least three consecutive days, a heat warning of danger level 3 is issued. If the threshold of daily maxima hi of 93°F is surpassed for at least five consecutive days a warning of danger level 4 is delivered [32]. These are the warning levels currently in use in Switzerland. In terms of the criteria used for issuing the different levels of risk, the Swiss heat warning system does not take into account any epidemiological analysis of the relationship between thermal stress indicator and health impact (mortality). Instead, it relates to the US heat warning system, also based on hi, [33], which was adopted after summer 2003. The National Oceanic and Atmospheric Administration (NOAA; [33]) issues heat warnings in four categories, scaled in terms of possible heat-related health issues: when hi is expected to exceed a threshold of 80–90 °F (category “caution”), when it surpasses 90–105 °F (category “extreme caution”), when hi reaches 105–110 °F for two consecutive days (category “danger”) and when it surpasses 130 °F (category “extreme danger”). Note that these thresholds used in the US depend on the local climate and vary among different regions. As high danger levels (category “danger” and “extreme danger”) used by NOAA do not apply to Switzerland, MeteoSwiss relates to the lower risk level of 90 °F only and implemented an additional risk level of 93 °F. Compared to higher thresholds used in the US system, these better reflect extreme heat conditions prevailing at regional level in Switzerland.
hi is generally considered to perform well and offers robust results in the context of weather forecasting [8]. However, the index is rather inflexible and limited in its usage for further applications. Its lower bound at the level of 80 °F, for instance, does not allow for the calculation of continuous climatologies. Moreover, no associated ISO standards exist. It is exactly here where our analyses try to add important value. In the present case study, we compare the current heat warning system to a (virtual) new system that makes use of different and potentially more appropriate heat stress indices, which allows for a more flexible use and diversification in terms of specific warning demands.

3. Materials and Methods

3.1. Indices and Data

Eight widely-used heat stress indices are tested against the currently used heat index (hi) to issue official heat warnings in Switzerland, namely: wet bulb temperature (wbt; [34]), wet bulb globe temperature in the shade (wbgt.shade; [35,36]), wet bulb globe temperature in the sun (wgbt.sun; [36,37]), simplified wet bulb globe temperature (swbgt; [25,38]), apparent temperature (apparentTemp; [21,38,39]; corresponding to the shade apparent temperature following Morabito et al. [40]), effective temperature (effectiveTemp; [24,41,42]), humidity index (humidex; [38,43]) and discomfort index (discomInd; [24,44]). Depending on the considered index, different meteorological input parameters are required (Table 1): 2 m temperature (temperature), solar shortwave downwelling radiation (radiation), 10 m wind speed (wind) and/or relative humidity (relative humidity). A fifth input variable, the dew point temperature (dew point), is derived from relative humidity and temperature data. All indices (except hi) are provided in degrees Celsius and allow for the calculation of continuous climatologies. The hi, in contrast, is set to zero by definition below a value of 80. The wbt and wbgt (in the shade and in the sun) have been thoroughly tested and approved already in related contexts. The first has been used for heat stress analyses in the recently published CH2018 Swiss Climate Scenarios [45]. The wbt represents the temperature registered by a sensor that is covered with a wetted cloth. At 100% relative humidity the wbt equals the (dry bulb) air temperature. The wbt can be calculated by means of thermodynamic equations [46] or by using the empirical formula used in [34], as done for the CH2018 scenarios [45]. The wbgt, in contrast, is a widely used index in indoor (wbgt.shade) and outdoor (wbgt.sun) occupational environments. The wbgt.sun is a combination of the natural wet bulb temperature (measured with a wetted thermometer exposed to the wind and heat radiation at the site, it represents the cooling of the body via sweat evaporation), the black globe temperature (measured inside a 150 mm diameter black globe, it represents the heat absorption from radiation), and the air temperature (measured with a “normal” thermometer shaded from direct heat radiation). The wbgt.shade constitutes an approximation of the former, assuming no strong source of radiation and a constant value for wind speed of 1 m/s [25]. Both indices provide values that can be modified to account for clothing and offer standardized reference values for acclimatized and unacclimatized people with recommended rest/work cycles at different metabolic rates (ISO 7243: 2017) [16,47]. Detailed information and implementation on the remaining indices employed is provided in the referenced literature. All the heat stress indices have been calculated with the R package HeatStress (https://github.com/anacv/HeatStress, doi:10.5281/zenodo.3264929.).
The observational data used for the index calculation are provided by a network of 28 automated stations in Switzerland (Figure 1), operated by the Swiss national weather service MeteoSwiss (SwissMetNet stations) and in accordance with World Meteorological Organization (WMO) standards. Due to these standards, MeteoSwiss does not operate weather stations in urban areas. All sites are consistent with the official heat warning-relevant elevation of <600 m in Switzerland [48]. We use long-term (1981–2017) observational data to compute the indices at hourly resolution in a retrospective approach. To compare with the official warnings based on daily maxima hi, we aggregate the data and focus on maximum daily heat stress values exclusively. By investigating the daily maximum heat stress, there are no major differences between urban and rural sites as there are for daily minimum temperatures, due to the UHI effect (e.g., [12,13,49]), which justifies the use of high-qualitative rural station data.

3.2. Threshold Calculation

For testing the current heat warning system against a virtual new system that makes use of a different heat stress index, we investigate how well warnings that are based on hi match warning days (days above the threshold) and warning periods (consecutive days above the threshold) that are calculated from the alternative heat stress indices under consideration. In order to identify these warnings with the alternative indices, index-specific thresholds need to be determined. We test two different approaches: frequency-based approach (site-specific thresholds, i.e., different for every station) and averaging-approach (Swiss-wide thresholds, i.e., identical for all stations), to compare to the official Swiss-wide threshold of hi exceeding 90 °F [32]. Both threshold approaches can be likewise applied to the danger level 4 warnings by using the hi threshold of 93 °F [32]. In this paper, however, we focus on the more common danger level 3 warnings exclusively.

3.2.1. Site-Specific Thresholds

Individual thresholds for each index and for each station are derived by using a frequency-based approach, meaning that the total number of official warning days at each station stays the same regardless of the index used. In more detail, we first identify the site-specific number of official warning days by counting the values exceeding the hi threshold of 90 °F. In a next step, site-specific values of the alternative indices are each sorted in ascending order. The one value that is followed by the exact same number of exceeding values as the number of official warning days is considered as the threshold value for the respective station and index.

3.2.2. Swiss-Wide Thresholds

Official heat warnings in Switzerland are based on a hi threshold valid for whole Switzerland [32]. To ensure comparability, Swiss-wide and site-independent thresholds for each index are derived by using a straightforward averaging approach. To reduce the influence of extremes and to gain more robust results, all station values of the respective index for the period between 1981–2017 are pooled together and compared to hi values exceeding 90 °F (Figure 2). The index values that correspond to the positions of 88 °F < hi < 92 °F are averaged and the obtained value is considered the Swiss-wide threshold for the employed index (Figure 2). In contrast to the site-specific thresholds, the assumption of having the same total number of official warning days regardless of the index used is no longer valid for the Swiss-wide threshold. Instead, the total number of warning days at a given station based on an alternative index typically differs from the total number of official warning days.

3.3. Matches of Warnings for Single and Consecutive Warning Days

3.3.1. Matches of Single Warning Days

We use a straightforward matching approach to identify the percentage of corresponding single warning days based on hi and the potential heat index alternatives. Official warning days are checked against days that indicate an index value exceeding the threshold (see Section 3.2) in order to examine the number of temporal matches. The sum of matching days relative to the total number of official warning days in the considered period 1981–2017 gives the overall percentage of matches (Figure 3). In order to obtain comprehensive site- and index-specific information on the percentage of matching single warning days, this approach is applied to each index and each station separately. The single-day matches calculation is performed based on the site-specific thresholds for the alternative indices (Section 3.2.1) to ensure the same frequency of warning days as for the hi.

3.3.2. Matches of Consecutive Warning Days

Official level 3 heat warnings in Switzerland are based on at least three consecutive days above the hi threshold of 90 °F [32]. In line with this regulation, we identify for each index and for each station warning periods of at least three days in a row. To better account for the different temporal durations of the official warning periods, we distinguish between three-day events, four-day events and five-(or more) day events. By comparing the warning periods based on hi and the potential heat index alternatives (both using site-specific and Swiss-wide thresholds), we can identify and quantify potential temporal overlaps, which we refer to as matches (Figure 4).
Concerning these warning matches, two kinds of information are relevant for us (Figure 4): (i) The “hit rate” provides information on the number of official warning periods captured by the alternative index relative to the total number of official warning periods within the hi respective category (three-day events, four-day events and five- (or more) day events), regardless of the correct duration of the warning event (100% for four-day events in the example, since one out of one event is captured). Besides the number of captured warnings, the “hit rate” identifies potentially “missed” official warning periods within the respective category; that is, no warnings are issued by the alternative index, despite an official warning period. (ii) Whenever there is a hit (overlap of two periods for at least one day), the percentage amount of overlapping warning days compared to the duration of the respective official warning period of a category is calculated (“quality of hit(s)”, 25% in the example below, as one out of four days are captured). On the basis of these categories (three-day events, four-day events and five- (or more) day events), we then average the percentages within one category for each index and each station (only values exceeding zero are considered) to receive the mean station- and index-specific “quality of hit(s)”.

4. Results

4.1. Threshold Calculation

4.1.1. Site-Specific Thresholds

Figure 5 shows the site-specific threshold anomalies, i.e., the local differences to the spatial mean, for all alternative indices considered as well as the spatial averages of the site-specific thresholds and the corresponding standard deviations.
The values obtained as reference thresholds change among indices as a result of the different input variables and their different combinations into a single number. Note, however, the partly pronounced spatial differences of thresholds within one index. This becomes obvious when considering the respective spatial standard deviations of the thresholds of wbt (standard deviation 0.49) or apparentTemp (standard deviation 0.60). In contrast, the site-specific thresholds of swbgt (standard deviation 0.05) or discomInd (standard deviation 0.09) vary markedly less and are thus more homogenous in space (Figure 5 and Table 2). The large spatial standard deviations of the site-specific indices are somewhat surprising since all index thresholds are derived by a frequency-based approach: by definition the number of days above the defined threshold of any alternative index and the number of days above the hi threshold are the same. The threshold of 90°F for hi has been applied to all stations. Thus, the threshold value of the respective index at “warmer” stations with higher hi values (i.e., more days above the threshold 90°F) is supposed to be similar to the threshold value of “colder” stations with lower hi values (i.e., less days above the threshold 90 °F). Hence, significant differences in absolute threshold levels within one index are not to be expected in the first place. However, since these considerations do not match the actual results of the majority of indices, we have proceeded with further analyses to better understand spatial threshold differences within a given index.
Large differences are noticeable, for instance, when comparing the absolute thresholds of wbt in Chur (21.5 °C) and Lugano (23.5 °C). We assume differing two-dimensional distributions of the input parameters temperature and relative humidity to be responsible for the differing thresholds of the two sites. A kernel density analysis allows for investigating the respective distributions for each station for the considered time period 1981–2017. The representation of the two-dimensional densities shows the probability of having different combinations of temperature and relative humidity (Figure 6). In Chur (Figure 6a), the histograms of both input parameters show unimodal-like patterns, whereas in Lugano (Figure 6b) the distributions (humidity in particular) feature two distinct maxima, resulting in different two-dimensional distribution patterns with two maxima in Lugano and only one maximum in Chur. The distribution of humidity in Chur, skewed to the right, is associated with a reduced probability of having high values of both variables at the same time, whereas in Lugano, relative humidity is multimodal with more weight on values of higher relative humidity, shifting the combined density towards warmer and moister conditions. These site-specific differences in distributional patterns, in addition to the differing mathematical combination of temperature and humidity in hi and the alternative index, are the reason for the spatial variation of the indices. The resulting distribution of wbt values in Chur (Figure 6c) and Lugano (Figure 6d) resembles very well the different distributional patterns of its input variables. In particular, temperature seems to be the dominant factor in determining the wbt distribution, pushing wbt in Lugano to higher density values in the upper tail. Local differences in thresholds result from differing distributions of input variables, i.e., localized meteorological characteristics. A country-wide threshold for such indices might overlook local characteristics related to different health risks.

4.1.2. Swiss-Wide Thresholds

Official heat warnings in Switzerland are based on a single hi threshold valid for the whole of Switzerland. Table 2 contains the corresponding Swiss-wide thresholds for the alternative indices (1) obtained as the spatial average of the site-specific thresholds (and the spatial standard deviation, as in Figure 5) and (2) derived by the straightforward averaging approach described in Section 3.2.2. Values range between 22.0 (wbt) and 36.6 °C (humidex). Comparing the two approaches, the Swiss-wide thresholds show very similar results (Table 2, first and third column).
Both the employed frequency-based approach to calculate site-specific thresholds and the averaging-approach to calculate Swiss-wide thresholds compare well to previous works. Results for the mean site-specific and the Swiss-wide thresholds of wbt (22.4 °C and 22.0 °C) are supported by the Swiss climate scenarios [45] that use a wbt threshold of 22.0 °C. This threshold has been derived as a high percentile of the present-day distribution of daily values at Swiss measurement sites. Also, apparentTemp thresholds (site-specific threshold of 32.6 °C and Swiss-wide threshold of 32.2 °C) resemble very well the threshold as calculated in [5], that is 32.0 °C (although their implementation does not consider the effect of wind). The value as used in [5] corresponds to the 98th percentile of the pooled maximum daily apparentTemp distribution across the eight considered meteorological stations in Switzerland.

4.2. Matches of Warnings for Single and Consecutive Warning Days

4.2.1. Matches of Single Warning Days (for Site-Specific Thresholds)

For all alternative indices, Figure 7 shows the percentage of matches of single warning days based on site-specific thresholds (Section 4.1.1) compared to official warning days based on hi values exceeding 90 °F.
Interestingly, among the indices, the percentages of matches differ considerably, with the spatial mean values ranging between 53.8 (wbt) and 95.5% (swbgt). To better explain the varying matching percentages among the indices and to identify the indices’ sensitivities towards their input variables, we conducted a response surface analysis (RSA, Figure 8). RSA is a statistical approach that examines the relationship between explanatory input variables and one (or more) response variable(s). By altering one input variable at a time, while the others remain constant, it is possible to identify the influence of each variable individually on the response variable. For an exemplary analysis we chose the indices with the lowest (wbt) and the highest (swbgt) matches, both based on the same input variables temperature and relative humidity. The results reveal similar sensitivities (slopes) towards temperature (Figure 8a,c), but in terms of humidity, wbt shows a stronger dependency than swbgt does (larger slope in Figure 8b compared to Figure 8d). We assume these differences in sensitivity to be responsible for the diverging percentages of matches. With swbgt yielding the highest matches, we further expect this index to best resemble the sensitivity of hi towards its input parameters. It is, however, important to mention that lower matches (e.g., wbt) should not be classified as poor results: Different sensitivities to the respective input variables might offer additional benefits for effective heat warnings, depending on the purpose of the warning.
The different sensitivities of the heat stress indices towards their input variables becomes apparent when considering the temporal evolution of the indices. As Figure 9 indicates (example of Basel, summer 2003), hi (Figure 9a) and swbgt (Figure 9c) seem to combine their input variables relative humidity and temperature in a very similar way, resulting in the exact same days above the respective thresholds (hi threshold of 90 °F and site-specific thresholds for alternative indices). This becomes clearer, when focusing on the first half of August 2003 when hi and swbgt show very high values for almost the whole period. In the case of wbt (Figure 9b), however, the high index values in the beginning of the period decrease much earlier (at the end of the first week of August 2003) than the hi and swbgt values do, leading to fewer values exceeding the wbt threshold during this period and thus to lower percentages of matches. Note that, by definition, the total number of warning days in the whole considered period of 1981–2017 stays the same regardless of the index used (see Section 3.2.1), meaning in total 119 warning days in Basel. Yet, zooming in to a shorter period results in differing numbers of warning days for the respective indices. In summer 2003, there are 17 official warning days, 5 wbt-based warning days and 17 swbgt-based warning days. In the case of wbt, warning days are less concentrated in summer 2003 and more widely spread among the considered period 1981–2017 than the swbgt-based warning days, that have the same temporal distribution as hi in this specific summer.

4.2.2. Matches of Consecutive Warning Days (for Site-Specific and Swiss-Wide Thresholds)

We choose the stations Basel (BAS), Chur (CHU), Geneva (GVE), Lugano (LUG) and Zurich Fluntern (SMA) (see Figure 1) for an exemplary analysis of matching consecutive warning days based on site-specific (Section 4.1.1) and Swiss-wide (Section 4.1.2) thresholds. Picking these sites is motivated by our effort to include representative stations of different climate regions in Switzerland. We focus again on swbgt and wbt.
Figure 10 presents matches of overlapping warning periods based on wbt (site-specific thresholds) compared to hi. The alternative index wbt is not able to capture all official warning events (low “hit rate”). Yet, when there is a hit, the “quality of hits” is larger than 70% in most of the cases.
As expected, warnings based on swbgt (site-specific thresholds) show higher “hit rates” than warnings based on wbt (Figure 11), since for swbgt the matches of single warning days already show clear similarities with the official warnings (Figure 7). In general, all warning periods based on swbgt show an overlap with the official warning periods (“hit rate” of 100%). Also, the quality of overlapping warnings is very high: The majority of warnings based on swbgt is consistent with the official warnings with mean “qualities of hits” of 100%. Still, there is an interesting exception for Chur (CHU). For the five- (or more) day warning events, even two swbgt-based warning periods overlap with the official warning period (“hit rate” of 200%). Nevertheless, the mean “quality of hit(s)” is low (35%). Information on the matches of overlapping warning periods based on the remaining indices compared to hi (Figure S1–S6) is provided as supplementary material.
Figure 12 illustrates the matches of overlapping warning periods in accordance with index-specific, but Swiss-wide thresholds for wbt (Table 2, third column). Bear in mind that, under this approach, the frequency of warning days is not necessarily the same for the hi and the alternative indices. Depending on the station, these matches (“hit rate” and “quality of hit(s)”) are higher (lower) compared to using site-specific thresholds, as a larger (smaller) number of days exceeds the warning level, which increases (reduces) the probability of overlaps with official warning periods. By focusing on selected stations (Figure 1) the pronounced station dependency of using Swiss-wide instead of site-specific thresholds becomes clear.
As seen in Figure 10, considering the category of four-day warning events in Lugano (LUG), for instance, one out of three official warning periods is captured when using wbt with its site-specific threshold of 23.5 °C (“hit rate”: 33%; mean “quality of hit”: 75%). Issuing warnings based on the wbt Swiss-wide threshold of 22.0 °C, as expected, raises the number of warnings in Lugano (from 111 to 509 days), due to a lowered threshold: According to Figure 12, three out of three official warning periods are then captured with a mean “quality of hits” of 100%. In Chur (CHU), in contrast, the individual threshold (21.5 °C) is lower than the Swiss-wide threshold (22.0 °C), resulting in a decreased number of warnings (from 70 to 33 days) and a reduced probability of overlaps. Consistent with our assumptions, the mean “quality of hit(s)” in the category of four-day events drops from 100% to 25%, with still capturing the official warning period. Yet, one needs to keep in mind that the change in the total number of warning days does not allow for a fair comparison in terms of the influence of using a site-specific or a Swiss-wide threshold, as the case of Lugano confirms: the higher “hit rates” and “quality of hits” stem from the higher total number of warning days when using a Swiss-wide threshold compared to a site-specific threshold, as it increases the probability of overlaps. Moreover, a large number of total warning days might also derive into many warning episodes which are not registered by the hi, increasing the so-called false alarm rate in the field of verification of weather forecasts.
This difference in the number of warning days is especially pronounced for indices with highly inhomogeneous site-specific warning levels like wbt (mean site-specific threshold 22.4 °C, standard deviation 0.49, Swiss-wide threshold 22.0 °C; Table 2) and less pronounced for indices with rather homogeneous site-specific thresholds like swbgt (mean site-specific threshold 25.5 °C, standard deviation 0.05, Swiss-wide threshold 25.4 °C; Table 2). For the latter, the total number of single warning days shows only minor changes between site-specific (WDsite) and Swiss-wide thresholds (WDCH): Basel (WDsite #119 vs. WDCH #131), Chur (WDsite #70 vs. WDCH #71), Geneva (WDsite #96 vs. WDCH #101), Lugano (WDsite #111 vs. WDCH #136) and Zurich (WDsite #31 vs. WDCH #38). The analysis of swbgt-based warning periods using the Swiss-wide threshold thus presents the exact same results for the considered five stations as the analysis based on site-specific thresholds does (Figure 11).
Still, the considered performance metrics (“hit rate” and “quality of hits”) are most informative when interpreted in a setting with site-specific thresholds in which, by definition, the total number of warning days for each index stays the same compared to the official number of warning days.

5. Summary and Conclusions

The current heat warning system in Switzerland is based on the heat index, hi, which performs well and offers robust results [8]. However, the index is rather inflexible and limited in its use for further applications: it presents a lower bound and does not allow for the calculation of continuous climatologies. Furthermore, no associated ISO standards exist. The desire towards a more flexible and impact-oriented heat warning system is related to the use of several thresholds or even several indices to be able to meet the large diversity of user demands.
In this paper, we compare the current heat warning system to a (virtual) new system that makes use of a different and potentially more appropriate heat stress index. For our analysis we consider eight widely used heat stress indices, all of them provided in degrees Celsius and subject to continuous distributions with no systematic lower bound. The respective single warning days and warning periods derived by the computed thresholds of each index are then analyzed in terms of potential matches or overlaps with official warnings.
For single warning days using site-specific thresholds, our results indicate that the percentage of matches varies considerably among indices. The different sensitivities towards the indices’ input variables might be the reason for the differing matches among the analyzed indices. We assume the index that best resembles the dependencies of hi towards its input variables offers highest matches. Among the employed indices, swbgt shows the highest (95.5%) spatially averaged matches, followed by discomInd (90.3%) and effectiveTemp (86.2%). wbt offers the lowest spatially averaged matches (53.8%).
In terms of consecutive warning days, two types of information are relevant for matches of warning periods: the “hit rate” (the relation of the absolute number of official warning periods that can be captured—at least one day of overlap—by the potential heat index alternative and the total number of official warning periods) and the “quality of hit(s)” (the mean percentage amount of overlapping warnings within one period relative to the duration of the respective official warning period).
When using wbt (site-specific thresholds), not all the official warnings are captured (low “hit rate”). Yet, when there is a hit, the “quality of hits” is larger than 70% in most of the cases. swbgt (site-specific thresholds) offers very high “hit rates” and “quality of hit(s)” among all warning categories, as already indicated by high matches of single warning days.
The influence of using a Swiss-wide threshold rather than a site-specific threshold, however, cannot be estimated in terms of the proposed performance metrics, as the total number of warning days at a given site may change substantially, when changing between the two employed thresholds (e.g., for wbt).
To summarize, swbgt resembles official heat warnings very well, both for individual and consecutive warning days. It even has further advantages: With having no lower boundary, the index allows for the calculation of continuous climatologies. However, other indices with higher dependencies on humidity like wbt, could give some additional information as well, despite their lower matches to the official warnings. Buzan et al. [38] confirm the different dependencies of heat stress metrics on their input variables: some may have a stronger dependency on regionally extreme temperature, while others have a stronger dependency on regionally extreme humidity. Through its major role in affecting the body’s hydration state humidity has a large effect on the heat-related human stress level: once the air temperature is higher than the skin temperature, the body can no longer cool via direct heat exchange, but needs to start evaporating, which is most efficient when a pronounced humidity gradient is present between the skin surface and the surrounding air. Otherwise, the body core temperature rises, causing severe health effects like heat exhaustion, heat cramps or even strokes [50]. An index placing weight not only on temperature but also strongly on humidity, could therefore be an advantage especially in heat-related human stress research areas.
Thus, regardless of the performance in terms of matches, the optimal index to use strongly depends on the purpose of the warning. According to Coccolo et al. [24], users need to choose the best thermal stress index according to their research’s needs, as a perfect solution to quantify the outdoor human comfort does not exist. Aiming at an effective, integrated and impact-oriented warning system, the combination of several indices, or at least different user-tailored thresholds depending on the specific warning demand, needs to be considered.
In most European countries the respective heat-health warning system (HHWS) is based on several thresholds at two to four risk levels, using single- or multiple-parameter methods, i.e., they use a single metric of temperature or multiple variables combined in one index [31]. The cited work [31] also found out that, in most cases, the warning levels are determined from epidemiological studies, linking heat to mortality (or other health indicators). Those relationships are established through statistical models, such as the distributed lag nonlinear model [51]. As the Swiss heat warning system does not take into account any epidemiological analysis, a way to improve the warning system could be to use a thermal stress indicator or heat stress indicator and to perform epidemiological analyses to identify location-specific thresholds to be then implemented operationally, building upon previous works such as [5]. This is important because, as described in a recent report of WHO-WMO on HHWS [9], usually the thresholds are response-specific, that means, the threshold values are determined at a level related to negative human response as suggested by the long-term relationship between some heat stress indicator and mortality. The threshold for action could then be the index value at which mortality (or any other health outcome) starts to rise substantially. Although this work is devoted to the heat warning system in Switzerland, it could serve as an example of application to provide support in the revision process of any heat warning system and paves the way to further reassessments of heat warnings.

Supplementary Materials

The following supporting information is available at https://0-www-mdpi-com.brum.beds.ac.uk/1660-4601/16/15/2684/s1, as part of the online article: Figure S1: As Figure 10 but for wbgt.shade. Figure S2: As Figure 10 but for wbgt.sun. Figure S3: As Figure 10 but for apparentTemp. Figure S4: As Figure 10 but for effectiveTemp. Figure S5: As Figure 10 but for humidex. Figure S6: As Figure 10 but for discomInd.

Author Contributions

All authors conceptualized the research and developed the methodology. A.B. carried out the analysis together with the co-authors A.C. and S.K. and conducted the majority of the manuscript writing. A.C., S.K. and C.S. reviewed and further edited the article. C.S. provided supervision during the entire research process.

Funding

Ana Casanueva was funded by the HEAT-SHIELD Project (HORIZON 2020, research and innovation program under the grant agreement 668786).

Acknowledgments

We thank Regula Gehrig (MeteoSwiss) for her valuable assistance in the context of heat stress indices. The authors also thank three anonymous referees for their useful comments that helped to improve the original manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fischer, E.M.; Schär, C. Consistent geographical patterns of changes in high-impact European heatwaves. Nat. Geosci. 2010, 3, 398–403. [Google Scholar] [CrossRef]
  2. Stocker, T.F.; Qin, D.; Plattner, G.-K.; Tignor, M.; Allen, S.K.; Boschung, J.; Nauels, A.; Xia, Y.; Bex, V.; Midgley, P.M. (Eds.) Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC); Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; pp. 1–1535. [Google Scholar]
  3. Morabito, M.; Crisci, A.; Messeri, A.; Messeri, G.; Betti, G.; Orlandini, S.; Raschi, A.; Maracchi, G. Increasing heatwave hazards in the Southeastern European Union capitals. Atmosphere 2017, 8, 115. [Google Scholar] [CrossRef]
  4. Cardona, O.D.; van Aalst, M.K.; Birkmann, J.; Fordham, M.; McGregor, G.; Perez, R.; Pulwarty, R.S.; Schipper, E.L.F.; Sinh, B.T. Determinants of risk: Exposure and vulnerability. In Managing the Risks of Extreme Events and Disasters to Advance. A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change (IPCC); Field, C.B., Barros, V., Stocker, T.F., Qin, D., Dokken, D.J., Ebi, K.L., Mastrandrea, M.D., Mach, K.J., Plattner, G.-K., Allen, S.K., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2012; pp. 65–108. [Google Scholar]
  5. Ragettli, M.S.; Vicedo-Cabrera, A.M.; Schindler, C.; Röösli, M. Exploring the association between heat and mortality in Switzerland between 1995 and 2013. Environ. Res. 2017, 158, 703–709. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Toloo, G.S.; FitzGerald, G.; Aitken, P.; Verrall, K.; Tong, S. Are heat warning systems effective? Environ. Health 2013, 12, 27. [Google Scholar] [CrossRef] [PubMed]
  7. Kjellstrom, T.; Weaver, H.J. Climate Change and health: Impacts, vulnerability, adaptation and mitigation. N. S. W. Public Health Bull. 2009, 20, 5–9. [Google Scholar] [CrossRef] [PubMed]
  8. Federal Office of Meteorology and Climatology (MeteoSwiss). Hitzebericht. Fachber. MeteoSchweiz, submitted.
  9. World Meteorological Organization (WMO). Heatwaves and Health: Guidance on Warning-System Development; McGregor, G.R., Ed.; WMO-No. 1142; World Meteorological Organization: Geneva, Switzerland, 2015; pp. 1–114. [Google Scholar]
  10. Pogačar, T.; Casanueva, A.; Kozjek, K.; Ciuha, U.; Mekjavić, L.B.; Kajfež Bogataj, L.; Črepinšek, Z. The effect of hot days on occupational heat stress in the manufacturing industry: Implications for workers’ well-being and productivity. Int. J. Biometeorol. 2018, 62, 1251–1264. [Google Scholar] [CrossRef] [PubMed]
  11. Gao, C.; Kuklane, K.; Östergren, P.-O.; Kjellstrom, T. Occupational heat stress assessment and protective strategies in the context of climate change. Int. J. Biometeorol. 2018, 62, 359–371. [Google Scholar] [CrossRef] [PubMed]
  12. Oke, T.R. The energetic basis of urban heat island. Q. J. R. Meteorol. Soc. 1982, 108, 1–24. [Google Scholar] [CrossRef]
  13. Parlow, E.; Vogt, R.; Feigenwinter, C. The urban heat island of Basel—Seen from different perspectives. DIE ERDE 2014, 145, 96–110. [Google Scholar]
  14. Morabito, M.; Messeri, A.; Noti, P.; Casanueva, A.; Crisci, A.; Kotlarski, S.; Orlandini, S.; Schwierz, C.; Spirig, C.; Kingma, B.; et al. An occupational heat warning system for Europe: The HEAT-SHIELD platform. Int. J. Environ. Res. Public Health 2019. (under review). [Google Scholar]
  15. Heat-Shield. About Heat-Shield. Available online: https://www.heat-shield.eu/about (accessed on 29 April 2019).
  16. International Organization for Standardization (ISO). Hot Environments—Ergonomics of the Thermal Environment—Assessment of Heat Stress Using the WBGT (Wet Bulb Globe Temperature) Index, ISO Standard 7243; International Standards Organization: Geneva, Switzerland, 2017; pp. 1–18. [Google Scholar]
  17. Pfafferott, J.; Becker, P. Erweiterung des Hitzewarnsystems um die Vorhersage der Wärmebelastung in Innenräumen. Bauphysik 2008, 30, 237–243. [Google Scholar] [CrossRef]
  18. German Meteorological Service (DWD). Hitzewarnung. Available online: https://www.dwd.de/DE/leistungen/hitzewarnung/hitzewarnung.html (accessed on 29 April 2019).
  19. Blazejczyk, 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]
  20. Parsons, K.C. Human Thermal Environments: The Effects of Hot, Moderate, and Cold Environments on Human Health, Comfort and Performance; Taylor and Francis: London, UK; New York, NY, USA, 2003; p. 527. [Google Scholar]
  21. Steadman, R.G. The assessment of sultriness. Part I: A temperature-humidity index based on human physiology and clothing science. J. Appl. Meteorol. 1979, 18, 861–873. [Google Scholar] [CrossRef]
  22. Rothfusz, L.P. The Heat Index Equation; NWS Southern Region Technical Attachment, SR/SSD 90–23; National Oceanic and Atmospheric Administration, National Weather Service, Office of Meteorology: Fort Worth, TX, USA, 1990. [Google Scholar]
  23. Moran, D.S.; Shitzer, A.; Pandolf, K.B. A physiological strain index to evaluate heat stress. Am. J. Physiol. Regul. Integr. Comp. Physiol. 1998, 275, 129–134. [Google Scholar] [CrossRef] [PubMed]
  24. Coccolo, S.; Kämpf, J.; Scartezzini, J.-L.; Pearlmutter, D. Outdoor human comfort and thermal stress: A comprehensive review on models and standards. Urban Clim. 2016, 18, 33–57. [Google Scholar] [CrossRef]
  25. Lemke, B.; Kjellstrom, T. Calculating Workplace WBGT from Meteorological Data: A Tool for Climate Change Assessment. Ind. Health 2012, 50, 267–278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Di Napoli, C.; Pappenberger, F.; Cloke, H. Assessing heat-related health risk in Europe via the Universal Thermal Climate Index (UTCI). Int. J. Biometeorol. 2018, 62, 1155–1165. [Google Scholar] [CrossRef] [Green Version]
  27. Fiala, D.; Havenith, G.; Bröde, P.; Kampmann, B.; Jendritzky, G. UTCI-Fiala multi-node model of human heat transfer and temperature regulation. Int. J. Biometeorol. 2012, 56, 429–441. [Google Scholar] [CrossRef]
  28. Jendritzky, G.; de Dear, R.; Havenith, G. UTCI–Why another thermal index? Int. J. Biometeorol. 2012, 56, 421–428. [Google Scholar] [CrossRef]
  29. Havenith, G.; Fiala, D.; Blazejczyk, K.; Richards, M.; Bröde, P.; Holmer, I.; Rintamaki, H.; Benshabat, Y.; Jendritzky, G. The UTCI-clothing model. Int. J. Biometeorol. 2012, 56, 461–470. [Google Scholar] [CrossRef]
  30. Emmanuel, R. Thermal comfort implications of urbanization in a warm-humid city: The Colombo Metropolitan Region (CMR), Sri Lanka. Build. Environ. 2005, 40, 1591–1601. [Google Scholar] [CrossRef]
  31. Casanueva, A.; Burgstall, A.; Kotlarski, S.; Messeri, A.; Morabito, M.; Flouris, A.; Nybo, L.; Spirig, C.; Schwierz, C. Overview of existing heat-health warning systems in Europe. Int. J. Environ. Sci. Public Health 2019. submitted. [Google Scholar] [CrossRef]
  32. Federal Office of Meteorology and Climatology (MeteoSwiss). Explanation of the Danger Levels: Heat Wave. Available online: https://www.meteoswiss.admin.ch/home/weather/gefahren/ explanation-of-the-danger-levels/heat-wave.html (accessed on 29 April 2019).
  33. National Oceanic and Atmospheric Administration (NOAA). Heat Index. Available online: https://web.archive.org/web/20110629041320/http:/www.crh.noaa.gov/pub/heat.php (accessed on 17 July 2019).
  34. Stull, R. Wet-Bulb Temperature from Relative Humidity and Air Temperature. J. Appl. Meteorol. Climatol. 2011, 50, 2267–2269. [Google Scholar] [CrossRef] [Green Version]
  35. Bernard, T.E.; Pourmoghani, M. Prediction of Workplace Wet Bulb Global Temperature. Appl. Occup. Environ. Hyg. 1999, 14, 126–134. [Google Scholar] [CrossRef] [PubMed]
  36. Yaglou, C.P.; Minard, D. Control of heat casualties at military training centers. Arch. Indust. Health 1957, 16, 302–316. [Google Scholar]
  37. Liljegren, J.C.; Carhart, R.; Lawday, P.; Tschopp, S.; Sharp, R. Modeling wet bulb globe temperature using standard meteorological measurements. J. Occup. Environ. Hyg. 2008, 5, 645–655. [Google Scholar] [CrossRef] [PubMed]
  38. Buzan, J.; Oleson, K.; Huber, M. Implementation and comparison of a suite of heat stress metrics within the Community Land Model version 4.5. Geosci. Model Dev. 2015, 8, 151–170. [Google Scholar] [CrossRef] [Green Version]
  39. Steadman, R.G. A Universal Scale of Apparent Temperature. J. Clim. Appl. Meteorol. 1984, 23, 1674–1687. [Google Scholar] [CrossRef]
  40. Morabito, M.; Crisci, A.; Messeri, A.; Capecchi, V.; Modesti, P.; Gensini, G.; Orlandini, S. Environmental Temperature and Thermal Indices: What Is the Most Effective Predictor of Heat-Related Mortality in Different Geographical Contexts? Sci. World J. 2014, 2014, 961750. [Google Scholar] [CrossRef]
  41. Houghton, F.C.; Yaglou, C.P. Determining equal comfort lines. J. Am. Soc. Heat. Vent. Eng. 1923, 29, 165–176. [Google Scholar]
  42. Missenard, F.A. Température Effective d’une Atmosphere Généralisation Température Résultante d’un Milieu. Encyclopédie Industrielle et Commerciale, Etude Physiologique et Technique de la Ventilation; Librerie de l’Enseignement Technique: Paris, France, 1933; pp. 131–185. [Google Scholar]
  43. Masterson, J.; Richardson, F.A. Humidex, a Method of Quantifying Human Discomfort Due to Excessive Heat and Humidity; Environment Canada: Downsview, ON, Canada, 1979. [Google Scholar]
  44. Thom, E.C. The discomfort index. Weatherwise 1959, 12, 57–61. [Google Scholar] [CrossRef]
  45. National Centre for Climate Services. CH2018—Climate Scenarios for Switzerland. In Tech. Rep; National Centre for Climate Services: Zurich, Switzerland, 2018; pp. 1–271. [Google Scholar]
  46. Davies-Jones, R. An Efficient and Accurate Method for Computing the Wet-Bulb Temperature along Pseudoadiabats. Mon. Weather Rev. 2008, 136, 2764–2785. [Google Scholar] [CrossRef]
  47. Lundgren Kownacki, K.; Gao, C.; Kuklane, K.; Wierzbicka, A. Heat Stress in Indoor Environments of Scandinavian Urban Areas: A Literature Review. Int. J. Environ. Res. Public Health 2019, 16, 560. [Google Scholar] [CrossRef] [PubMed]
  48. Federal Office of Meteorology and Climatology (MeteoSwiss). Der Hitzesommer 2015 in der Schweiz. Fachber. MeteoSchweiz 2016, 272, 1–47. [Google Scholar]
  49. Federal Office of Meteorology and Climatology (MeteoSwiss). Städtische Wärmeinseln in der Schweiz—klimatologische Studie mit Messdaten in fünf Städten. Fachber. MeteoSchweiz 2018, 273, 1–62. [Google Scholar]
  50. Davis, R.E.; McGregor, G.R.; Enfield, K.B. Humidity: A review and primer on atmospheric moisture and human health. Environ. Res. 2016, 144, 106–116. [Google Scholar] [CrossRef] [PubMed]
  51. Gasparrini, A.; Armstrong, B.; Kenward, M.G. Distributed lag non-linear models. Stat. Med. 2010, 29, 2224–2234. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location and elevation of the 28 considered automated stations of the Swiss national weather service MeteoSwiss (SwissMetNet stations). The selected sites are consistent with the heat warning-relevant elevation of <600 m [48]. The stations Basel (BAS), Chur (CHU), Geneva (GVE), Lugano (LUG) and Zurich Fluntern (SMA) serve as exemplary sites for later analyses.
Figure 1. Location and elevation of the 28 considered automated stations of the Swiss national weather service MeteoSwiss (SwissMetNet stations). The selected sites are consistent with the heat warning-relevant elevation of <600 m [48]. The stations Basel (BAS), Chur (CHU), Geneva (GVE), Lugano (LUG) and Zurich Fluntern (SMA) serve as exemplary sites for later analyses.
Ijerph 16 02684 g001
Figure 2. Exemplary Swiss-wide threshold calculation for wet bulb temperature (wbt) on the basis of 88 °F < hi < 92 °F. The markers represent virtual wbt and heat index (hi) pairs for all stations together.
Figure 2. Exemplary Swiss-wide threshold calculation for wet bulb temperature (wbt) on the basis of 88 °F < hi < 92 °F. The markers represent virtual wbt and heat index (hi) pairs for all stations together.
Ijerph 16 02684 g002
Figure 3. Exemplary matching of warning days (n = 10) based on hi (left) and wbt (right) for a station x. The black boxes mark matching warnings. Two out of ten wbt-based warning days match the official hi-based warning days, resulting in an overall 20% match using wbt in station x.
Figure 3. Exemplary matching of warning days (n = 10) based on hi (left) and wbt (right) for a station x. The black boxes mark matching warnings. Two out of ten wbt-based warning days match the official hi-based warning days, resulting in an overall 20% match using wbt in station x.
Ijerph 16 02684 g003
Figure 4. Exemplary warning periods for a station x, based on hi (left) and wbt (right). The warning periods are separated into three-, four- and five- (or more) day events according to their lengths. The black boxes mark warning periods that offer at least one day of overlapping warnings (a). In this example, the official four-day warning event is captured by one overlap (5th June 2017) when using wbt, which results in a “hit rate” of 100% in the 4-day event category and a “quality of hit(s)” of 25% (b).
Figure 4. Exemplary warning periods for a station x, based on hi (left) and wbt (right). The warning periods are separated into three-, four- and five- (or more) day events according to their lengths. The black boxes mark warning periods that offer at least one day of overlapping warnings (a). In this example, the official four-day warning event is captured by one overlap (5th June 2017) when using wbt, which results in a “hit rate” of 100% in the 4-day event category and a “quality of hit(s)” of 25% (b).
Ijerph 16 02684 g004
Figure 5. Spatial anomalies of site-specific thresholds (color-coded) and spatial averages and standard deviations of the site-specific thresholds (on top of each map) for the considered indices wet bulb temperature (wbt), wet bulb globe temperature in the shade (wbgt.shade), wet bulb globe temperature in the sun (wgbt.sun), simplified wet bulb globe temperature (swbgt), apparent temperature (apparentTemp); effective temperature (effectiveTemp), humidity index (humidex) and discomfort index (discomInd). Anomalies are computed wrt. the respective spatially averaged site-specific threshold. Blue colors apply to negative threshold anomalies, red colors to positive threshold anomalies.
Figure 5. Spatial anomalies of site-specific thresholds (color-coded) and spatial averages and standard deviations of the site-specific thresholds (on top of each map) for the considered indices wet bulb temperature (wbt), wet bulb globe temperature in the shade (wbgt.shade), wet bulb globe temperature in the sun (wgbt.sun), simplified wet bulb globe temperature (swbgt), apparent temperature (apparentTemp); effective temperature (effectiveTemp), humidity index (humidex) and discomfort index (discomInd). Anomalies are computed wrt. the respective spatially averaged site-specific threshold. Blue colors apply to negative threshold anomalies, red colors to positive threshold anomalies.
Ijerph 16 02684 g005
Figure 6. The 2-dimensional Kernel density plots of temperature and relative humidity (a,b) and the resulting distribution of wbt (c,d) for Chur and Lugano for the period 1981–2017. In Figure 6a,b, the top and right axes show the histograms of temperature (red) and relative humidity (blue). The density plots combine the distributions of both variables with red colors referring to higher probabilities and blue colors to unlikely conditions.
Figure 6. The 2-dimensional Kernel density plots of temperature and relative humidity (a,b) and the resulting distribution of wbt (c,d) for Chur and Lugano for the period 1981–2017. In Figure 6a,b, the top and right axes show the histograms of temperature (red) and relative humidity (blue). The density plots combine the distributions of both variables with red colors referring to higher probabilities and blue colors to unlikely conditions.
Ijerph 16 02684 g006
Figure 7. Percentage of matches of single warning days based on site-specific thresholds compared to official warning days based on the hi exceeding 90 °F, analyzed for the considered indices wbt, wbgt.shade, wbgt.sun, swbgt, apparentTemp, effectiveTemp, humidex, discomInd. Blue colors apply to lower matches, red colors to higher matches.
Figure 7. Percentage of matches of single warning days based on site-specific thresholds compared to official warning days based on the hi exceeding 90 °F, analyzed for the considered indices wbt, wbgt.shade, wbgt.sun, swbgt, apparentTemp, effectiveTemp, humidex, discomInd. Blue colors apply to lower matches, red colors to higher matches.
Ijerph 16 02684 g007
Figure 8. Response surface analysis of wbt (a,b) and swbgt (c,d).
Figure 8. Response surface analysis of wbt (a,b) and swbgt (c,d).
Ijerph 16 02684 g008
Figure 9. The daily evolution of hi (a), wbt (b) and swbgt (c) in Basel in summer 2003 and the respective warnings days: Green vertical lines (ac) and green dots (a) show the official warning days based on hi, blue dots indicate warning days based on wbt (b) and red dots show warning days based on swbgt (c). Dashed horizontal lines mark the site-specific thresholds of each index.
Figure 9. The daily evolution of hi (a), wbt (b) and swbgt (c) in Basel in summer 2003 and the respective warnings days: Green vertical lines (ac) and green dots (a) show the official warning days based on hi, blue dots indicate warning days based on wbt (b) and red dots show warning days based on swbgt (c). Dashed horizontal lines mark the site-specific thresholds of each index.
Ijerph 16 02684 g009
Figure 10. Matches of overlapping warning periods based on wbt (site-specific thresholds) compared to hi. The warning periods result from the use of site-specific thresholds (i.e., fixed total number of warning days for hi and every alternative index, WDsite, see Section 3.2.1): BAS (WDsite #119), CHU (WDsite #70), GVE (WDsite #96), LUG (WDsite #111) and SMA (WDsite #31). The figure presents matches for three- (yellow), four- (orange) and five- (or more, dark orange) day official warning events. The shaded bars give the number of official warning periods of the respective duration. The dashed bars show the number of warning periods captured by wbt (indicating the “hit rate” when compared to the length of the entire bar). The percentage number on top of each bar gives the mean “quality of hit(s)” (Section 3.3.2). In Basel (BAS) and for the case of three-day warnings, for instance, four warning periods based on wbt show an overlap with (in total) eight official three-day warning events (“hit rate” of 50%). In the mean, these overlapping warnings show a 92% match (“quality of hits”). The asterisk applies to stations without any official warning periods in the respective category.
Figure 10. Matches of overlapping warning periods based on wbt (site-specific thresholds) compared to hi. The warning periods result from the use of site-specific thresholds (i.e., fixed total number of warning days for hi and every alternative index, WDsite, see Section 3.2.1): BAS (WDsite #119), CHU (WDsite #70), GVE (WDsite #96), LUG (WDsite #111) and SMA (WDsite #31). The figure presents matches for three- (yellow), four- (orange) and five- (or more, dark orange) day official warning events. The shaded bars give the number of official warning periods of the respective duration. The dashed bars show the number of warning periods captured by wbt (indicating the “hit rate” when compared to the length of the entire bar). The percentage number on top of each bar gives the mean “quality of hit(s)” (Section 3.3.2). In Basel (BAS) and for the case of three-day warnings, for instance, four warning periods based on wbt show an overlap with (in total) eight official three-day warning events (“hit rate” of 50%). In the mean, these overlapping warnings show a 92% match (“quality of hits”). The asterisk applies to stations without any official warning periods in the respective category.
Ijerph 16 02684 g010
Figure 11. As Figure 10 but for swbgt.
Figure 11. As Figure 10 but for swbgt.
Ijerph 16 02684 g011
Figure 12. Matches of overlapping warning periods based on wbt (Swiss-wide threshold) compared to hi. The warning periods result from the use of Swiss-wide thresholds (i.e., the number of warning days is different for the hi and the alternative index, see Section 3.2.2). Compared to the number of warning days calculated with the site-specific threshold (WDsite), the number of warning days changes individually for each station when using a Swiss-wide threshold (WDCH): BAS (WDsite #119 vs. WDCH #154), CHU (WDsite #70 vs. WDCH #33), GVE (WDsite #96 vs. WDCH #53), LUG (WDsite #111 vs. WDCH #509) and SMA (WDsite #31 vs. WDCH #38). The figure presents matches for three- (yellow), four- (orange) and five- (or more, dark orange) day official warning events. The shaded bars give the number of official warning periods for the respective duration. The dashed bars show the number of warning periods captured by wbt (indicating the “hit rate” when compared to the length of the entire bar). The percentage on each bar gives the mean “quality of hit(s)”. In Basel (BAS) and for the case of three-day warnings, for instance, six warning periods based on wbt show an overlap with (in total) eight official three-day warning events (“hit rate” of 75%). In the mean, these overlapping warnings show a 94% match (“quality of hits”). The asterisk applies to stations without any official warning periods in the respective category.
Figure 12. Matches of overlapping warning periods based on wbt (Swiss-wide threshold) compared to hi. The warning periods result from the use of Swiss-wide thresholds (i.e., the number of warning days is different for the hi and the alternative index, see Section 3.2.2). Compared to the number of warning days calculated with the site-specific threshold (WDsite), the number of warning days changes individually for each station when using a Swiss-wide threshold (WDCH): BAS (WDsite #119 vs. WDCH #154), CHU (WDsite #70 vs. WDCH #33), GVE (WDsite #96 vs. WDCH #53), LUG (WDsite #111 vs. WDCH #509) and SMA (WDsite #31 vs. WDCH #38). The figure presents matches for three- (yellow), four- (orange) and five- (or more, dark orange) day official warning events. The shaded bars give the number of official warning periods for the respective duration. The dashed bars show the number of warning periods captured by wbt (indicating the “hit rate” when compared to the length of the entire bar). The percentage on each bar gives the mean “quality of hit(s)”. In Basel (BAS) and for the case of three-day warnings, for instance, six warning periods based on wbt show an overlap with (in total) eight official three-day warning events (“hit rate” of 75%). In the mean, these overlapping warnings show a 94% match (“quality of hits”). The asterisk applies to stations without any official warning periods in the respective category.
Ijerph 16 02684 g012
Table 1. Employed heat stress indices with respective input parameters.
Table 1. Employed heat stress indices with respective input parameters.
Heat Stress IndexTemperatureRadiationWindRelative HumidityDew Point
hix x
wbtx x
wbgt.shadex x
wbgt.sunxxx x
swbgtx x
apparentTempx xx
effectiveTempx xx
humidexx x
discomIndx x
Table 2. Considered heat stress indices with respective mean site-specific (including the spatial standard deviation) and Swiss-wide threshold.
Table 2. Considered heat stress indices with respective mean site-specific (including the spatial standard deviation) and Swiss-wide threshold.
Heat Stress IndexSite-Specific Threshold (°C)Swiss-Wide Threshold (°C)
MeanStandard Deviation
hi--32.2 *
wbt22.40.4922.0
wbgt.shade25.00.2124.8
wbgt.sun27.90.3827.6
swbgt25.50.0525.4
apparentTemp32.60.6032.2
effectiveTemp25.70.3025.4
humidex36.90.4036.6
discomInd26.10.0926.0
* or 90 °F.

Share and Cite

MDPI and ACS Style

Burgstall, A.; Casanueva, A.; Kotlarski, S.; Schwierz, C. Heat Warnings in Switzerland: Reassessing the Choice of the Current Heat Stress Index. Int. J. Environ. Res. Public Health 2019, 16, 2684. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152684

AMA Style

Burgstall A, Casanueva A, Kotlarski S, Schwierz C. Heat Warnings in Switzerland: Reassessing the Choice of the Current Heat Stress Index. International Journal of Environmental Research and Public Health. 2019; 16(15):2684. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152684

Chicago/Turabian Style

Burgstall, Annkatrin, Ana Casanueva, Sven Kotlarski, and Cornelia Schwierz. 2019. "Heat Warnings in Switzerland: Reassessing the Choice of the Current Heat Stress Index" International Journal of Environmental Research and Public Health 16, no. 15: 2684. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152684

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