Next Article in Journal
Trees Diversity and Species with High Ecological Importance for a Resilient Urban Area: Evidence from Cotonou City (West Africa)
Next Article in Special Issue
Analysis of the Composite Risk Grade for Multi Extreme Climate Events in China in Recent 60 Years
Previous Article in Journal
Evaluation of Gridded Meteorological Data for Crop Sensitivity Assessment to Temperature Changes: An Application with CERES-Wheat in the Mediterranean Basin
Previous Article in Special Issue
Extreme Recreational Conditions in the Black Sea Resorts Associated with the North Atlantic Climate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Sensitivity Analysis of Heat Stress Indices

Laboratory of Innovative Technologies, University of Picardie Jules Verne, 80000 Amiens, France
*
Author to whom correspondence should be addressed.
Climate 2023, 11(9), 181; https://doi.org/10.3390/cli11090181
Submission received: 19 June 2023 / Revised: 19 August 2023 / Accepted: 29 August 2023 / Published: 30 August 2023
(This article belongs to the Special Issue Climate and Weather Extremes: Volume II)

Abstract

:
More than 40 heat indices are being used across the world to quantify outdoor thermal comfort. The selection of an Outdoor Heat Stress Index (OHSI) depends on several parameters, including clothing, age, awareness, local environment, food consumption, human activities, and resources. This study investigates various indicators of heat stress, including (i) OHSIs officially used to quantify heat stress worldwide, (ii) the estimation methods of these indices, and (iii) the sensitivity analysis of indices, namely, Corrected Effective Temperature (CET), Heat Index (HI), Wet Bulb Globe Temperature (WBGT), Universal Thermal Climate Index (UTCI), Discomfort Index (DI), Summer Simmer Index (SSI), and Predicted Mean Vote (PMV). The results indicate the degree of sensitivity of indices, with the HI being the most sensitive for estimating heat stress. Additionally, the WBGT, HI, and CET are recommended indices that can be directly measured using sensors instead of relying on calculated indices that are based on estimation techniques and some ideal physical assumptions.

1. Introduction

Heat stress is a condition that can result from prolonged exposure to high temperatures and humidity. It can range in severity from mild discomfort to a life-threatening medical emergency, depending on the individual’s susceptibility and the intensity and duration of the heat exposure [1]. It usually occurs when the body is unable to maintain a healthy temperature in response to a hot environment [2]. Heat stress can negatively affect public health and quality of life in both urban and rural areas [3,4]. It can lead to heavy sweating, muscle cramps, headaches, and decreased productivity [5,6]. Moreover, exposure to extreme heat can harm the body and lungs, exacerbating respiratory conditions such as asthma, particularly when breathing hot and humid air [7].
In different countries, heat stress is quantified with the help of specific indices. The thermal stress index is a quantitative measure that integrates one or more of the physical, thermal, and personal factors effecting heat transfer between the environment and a person [8]. Many heat stress indices have been developed and classified based on thermal comfort assessment, physiological strain, physical factors of the environment, and “rational” heat balance equation [9]. The thermal comfort of an individual depends mainly on the individual’s activities and other factors including behavioral activities, age group, health condition, and local climate and environment [10]. The collective factors are presented in Table 1. Elderly individuals and young children are particularly sensitive to heat stress due to their weaker tolerance for extreme heat events and greater need for comfort [11].
In relation to heat stress, rational heat measures refer to methods for assessing the thermal environment and predicting its effects on human health and comfort [12]. The factors listed in Table 1 can be used to evaluate the vulnerability of residents to heat stress. To address heat stress issues, different countries have adopted various techniques for assessing thermal comfort, such as Wet Bulb Global Temperature (WBGT), Heat Stress Index (HSI), and Thermal Work Limit (TWL). These indices are commonly used by climate, environmental, and national weather agencies [13]. Most indices can be directly measured and calibrated using equations, while others require the use of different models and methods.
Urban areas are often subject to higher temperatures due to factors such as the absorption and re-emission of heat by buildings and pavement, reduced vegetation cover, air pollution, and increased energy use. However, it should be noted that rising temperatures can also affect rural areas, particularly in the context of climate change, with some rural areas warming up faster than nearby urban areas [14,15].
This study is justified by the lack of arguments on the use of so many indices existing throughout the world to quantify the heat stress indices in different regions. Our hypothesis is that some parameters have a significant impact on different heat stress indices, while others may have a negligible effect.
Our research questions were as follows: What is the inventory of different External Heat Stress Indices (OHSI) around the world?; What are their current estimation methods?; How do CET, HI, SSI, PMV, DI, WBGT, and UTCI react to changes in different parameters?; What parameter variations induce changes from comfortable zone to warmer zones in the available thermal comfort scales? Our goal was to analyze the most commonly used ones by understanding their variation trend and exhibiting the most sensitive heat index parameters.
Therefore, the objectives of this study are to conduct a survey for OHSI officially used around the world, to analyze estimation methods of heat stress and currently available modified thermal comfort scales, and to perform a sensitivity analysis of CET, HI, SSI, PMV, DI, WBGT, and UTCI.
The novelty of this study lies in its comprehensive analysis of various heat stress indices used around the world, providing a unique perspective on their sensitivity to different relevant parameters. This sensitivity analysis is carried out for indices with available formula in two ways: 1. By using the partial derivatives with respect to their explanatory variables; and 2. By exploring the effect on the comfort zones of the realistic variations of the variables involved in the heat stress index considered. The first is theoretically valid for infinitesimal variations, while the second is obtained by simulating the value of the thermal stress index for chosen feasible values of its variables.
The structure of this paper is as follows: Section 2 presents the research methodology, Section 3 provides direct formulas for estimating different heat stress indices, Section 4 presents a survey map of the OHSI used across the world by official agencies, Section 5 and Section 6 discusses the results of sensitivity and thermal comfort scale variations, and Section 7 provides the conclusions and future perspectives.

2. Materials and Methods

The adopted methodology in our current study is described in Figure 1. The data were collected through OHSI surveys on internet sources, listed in Appendix A Table A1. Following this, sensitivity analysis of OHSI was performed based on updated estimation models.
The sensitivity of heat stress indices is paramount in understanding their real-world implications. A highly sensitive index might offer precise readings in specific conditions but may falter in others. On the other hand, a less sensitive index might provide a broader overview but miss nuances. This study delves into the sensitivity of various indices, highlighting their strengths and weaknesses, and offering recommendations for their optimal use.
The mathematical models of heat stress indices were computed to quantify the sensitivity of indices and evolving variables, such as air temperature and humidity. In addition, the process of identifying the range of values for each variable that corresponds to thermal comfort was conducted using minimum and maximum functions.

3. Estimation of Outdoor Thermal Indices

Heat stress can be quantified using more than 40 indices worldwide [16]. These indices are classified into three groups: “direct,” “empirical,” and “rational.” While direct indices are simpler and can be monitored using instruments like Kestrel [17], KIMO [18], and WBGT meter [19] in both indoor and outdoor environments, rational and empirical indices are more complex. These indices take into account meteorological, physiological, and behavioral factors, making it difficult to estimate thermal stress for an individual. To differentiate between empirical and rational indices, empirical indices rely on statistical analysis of observed data, while rational indices use physical models to estimate heat stress. Table 2 provides calibration methods and formulas for the 16 most commonly used OHSIs depicted from the relevant references in this paper. Other possible OHSIs may be studied by the reader by following our proposed methodology.

4. Outdoor Heat Indices in Different Regions

Different regions in the world, characterized by unique climatic conditions such as sub-tropical, tropical, and Mediterranean climates, influence human thermal conditions and physiology in distinct ways. Many countries monitor heat stress using a variety of parameters and indicators to inform locals about heat event warnings. Parameters like air temperature, relative humidity, global radiation, and wind speed are commonly used in the measurement of heat stress indices.
The WBGT is the most widely used index in various countries including Australia, the U.S.A, Europe, Japan, and Colombia. This is primarily because it can be measured by sensors and estimated by a simple mathematical equation. The choice of the heat index depends on several factors specific to an area such as weather patterns (hot and dry, semi-humid, cold winter regions, etc.), the age group of the population (elderly, children, youth), and activities.
To provide a comprehensive view of the indices used worldwide, we conducted a survey study. The results of this survey, including the agencies officially involved in the assessment of heat stress, are presented in Appendix A and visually represented in Figure 2.

5. Sensitivity Analysis

5.1. Sensitivity Analysis of Heat Indices versus Variables

In this section, we aim to investigate the sensitivity of heat stress indices to small variations around a reference point. To achieve this, we employ a multi-variable partial differential equation to analyze the sensitivities of the following operational indices: HI, WBGT, PMV, SSI, CET, DI, and UTCI. Simulations are started by defining a set of values ( S 0 ) for the input parameters that correspond to the comfort zone.
Table 3 gives the values for S0, where the input values for each index have been selected based on a comfortable situation. For example, with Tgo = 45 °C and Two = 16.27 °C, CET indicates comfortable conditions, but small variations could lead to slight discomfort.
Subsequently, the sensitivity of each index is computed with respect to small variations around the reference values.
Consider a heat index y as a function of explicative variables x 1 , , x n , that is,
y = f x 1 , , x n  
where y (thermal index) is the function of x (input variables e.g., T g ,   T a ,   T m r t ,   R H ,   v ,   T w ,   m e t ,   m w ). Then, a small variation of y can be expressed as
Δ y = f x 1 0 Δ x 1 + f x 2 0 Δ x 2 + + f x n 0 Δ x n
This study is relevant to know the trend of the variation and to exhibit the most sensitive parameters in the heat index:
  • A positive sign means that the index increased with the parameter, whereas a negative sign means that the index increases (resp. decreases) when the parameter decreases (resp. increases);
  • A high value of sensitivity means a high influence of the corresponding parameter on the given index;
  • We can see that relative humidity has a negligible impact on DI and UTCI.

5.2. Variations in Thermal Comfort Zones versus Heat Indices

In this section, we analyze the influence of variables for widely used heat indices. More precisely, simulations have been performed for the parameter ranges for each thermal comfort zone defined according to the referenced Table 4.
The formulas presented in Table 2 from references [20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38] were used to analyze the sensitivity of discomfort zones. The input mentioned in Table 3 was used to simulate each index, and it was observed that the HI, WBGT, and UTCI results indicated a different comfort zone, as shown in Table 5 with color codes. It was also noticed that, due to the significant variation coefficient of HI, a small increase in its variables could lead to an increase in the level of discomfort, while the results for other indices remained in the same comfort zone. Further, x1, … xn were simulated by using the min–max function, where xmin is the lowest possible (S0) and xmax is maximum possible value under realistic situations in the summer season ( T g = 50 ,   T a = 50   ° C ,   T m r t = 50   ° C ,   R H = 100 % ,   v = 10 m / s ,   T w = 25   ° C ,   m e t = 4 ,   m w = 0.9 ,   α = 3 ). The simulations were performed by changing only one variable at each time step. The obtained results are plotted and shown in Figure 3, which demonstrates the sensitivity of heat indices showing different thermal discomfort zones within the same range. The use of sensors and instruments for measuring the variables can significantly reduce errors and improve the accuracy of these simulations. This approach provides a more comprehensive understanding of how different variables can influence thermal discomfort, which is crucial for designing appropriate intervention strategies to mitigate the negative impacts of heat stress on human health and productivity. Therefore, it is recommended to use sensors and measure the variables accurately and continuously to obtain reliable data. Simultaneously, this can avoid over- and underestimation of heat stress.

6. Discussion

The findings of this research, as discussed in Table 6, shed light on the effectiveness of simulated indices, their advantages, associated challenges, and recommendations for their application. The innovation of this paper lies in its comprehensive analysis of heat stress indices, provides understanding that can guide policymakers and researchers.

7. Conclusions

This study concludes that every region across the globe chooses heat indices according to some specific parameters. The factors that influence the selection are awareness among people (especially vulnerable ones), physical interviews (knowing about the age, sex, sensations, clothing styles, and activities), the correlation between immune systems, and the number of heart-health events. While numerous countries prioritize other environmental concerns overheat stress indices, it is essential to recognize that each country’s choice of index is influenced by its unique climatic conditions and challenges. The comprehensive analysis of seven OHSIs (CET, WBGT, UTCI, HI, DI, SSI, PMV) was programmed and simulated to analyze the variation coefficient of their evolving parameters. The study provides a foundation for understanding of indices effectiveness and advantages. While each index has its unique strengths and limitations, their collective evaluation offers a holistic perspective on heat stress assessment. The results show that air temperature is the most sensitive parameter, especially for estimating HI, SSI, DI, and UTCI, whereas relative humidity is negligibly sensitive except in the case of SSI and HI. It is crucial for countries to consider mean radiant temperature in their indices, given its profound impact on perceived temperature and comfort.
In conclusion, we recommend the use of sensor-measurable indices such as WBGT, HI, and CET. Sensors and instruments are a more reliable method for measuring heat stress compared to mathematical models. It can avoid error by directly measuring parameters such as temperature, humidity, and air velocity, while models rely on assumptions that may not accurately reflect the actual conditions. Real-time measurements from sensors and instruments are critical in rapidly changing environments, such as industrial settings or during athletic events, and they provide more detailed information about the distribution of heat stress. While mathematical models have their advantages, they are often based on simplifications and assumptions that may not be as reliable as direct measurements. Overall, sensors and instruments are the preferred choice for accurately measuring heat stress in real-world settings.

Author Contributions

The paper was a collaborative effort between the authors A.R. and A.M.Q. contributed collectively to developing the methodology of this study, survey, sensitivity analysis, and the manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This paper has been produced within the COOL-TOWNS (Spatial Adaptation for Heat Resilience in Small and Medium Sized Cities in the 2 Seas Region) project which receives funding from the Interreg 2 Seas programme 2014-2020 co-funded by the European Regional Development Fund under subsidy contract N° 2S05-040.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. List of Outdoor thermal Stress Indices are officially used in Countries.
Table A1. List of Outdoor thermal Stress Indices are officially used in Countries.
Region/Country NameTSIEvolve Parameters Sources Agency Name
ColumbiaWBGT T w ,   T a ,   T g ,   R H ,   v https://columbiaweather.com/products/weather-stations/wet-bulb-globe-temperature/ (accessed on 4 November 2021). Columbia weather systems
Dallas WBGThttps://perryweather.com/weather-station/ (accessed 4 November on 2021).Perry weather
Tulsa https://www.weather.gov/arx/wbgt (accessed on 4 November 2021).National weather services
CaliforniaSSI (Pacific Ocean and Mediterranean weather) R H ,   T a in Fahrenheithttps://www.weather.gov/ (accessed on 8 October 2021).
Maryland Heat index (elder population)https://www.weather.gov/safety/heat-index (accessed on 11 November 2021).NOAA National Weather Services (Weather prediction center)
MiamiOxford Index T w ,   T a https://www.miamioh.edu/cas/academics/centers/erc/weather-station/index.html (accessed on 11 November 2021).Ecology Research Center
Eastern and western regionsHSI T a ,   R H ,   M ,   convective   heat   exchange ,   radiant   heat   exchange ,   v https://www.ncdc.noaa.gov/societal-impacts/heat-stress/climatology (accessed on 11 November 2021).National Centers for Environmental Information
AustraliaWBGT (for civilians) T w ,   T a ,   T g ,   R H ,   v http://www.bom.gov.au/ (accessed on 11 November 2021).Australian Government Bureau of meteorology
AT (for workers) T a ,   R H ,   v
TWL Physiological   data   like   height ,   age ,   sweat   rate ,   T a ,   T w ,   T g .
Japan (Tokyo)WBGT T w ,   T a ,   T g ,  RH, WShttps://mainichi.jp/english/articles/20180719/p2a/00m/0na/004000c (accessed on 18 November 2021).Ministry of environment
Europe UTCI (semi humid hot summer and cold winter regions) T a , Average radiation temperature, WS, RHhttps://climate-adapt.eea.europa.eu/knowledge/european-climate-data-explorer/health/thermal-comfort-indices-universal-thermal-climate-index#details (accessed on 15 December 2022). European Environmental agency
Required SweatingRH, Sweat rate, Mhttps://www.weather.gov/oun/safety-summer-heathumidity (accessed on 11 November 2021). National Weather service (W.S)
South Africa PMV (overrated summertime sensation and underestimated in winter) I c l ,   M ,   R H , T m r t .https://customweather.com/
(accessed on 11 November 2021).
Pietermaritzburg, South Africa
ChinaSET (subtropical regions) R H ,   v ,   M ,   radiation   temperature ,   T w ,   T a http://en.weather.com.cn/ (accessed on 15 December 2022). Weather China
Sweden (Stockholm) HSI (elderly population) T a , RH, M, convective heat exchange, radiant heat exchange, WShttps://www.smhi.se/en/q/Stockholm/2673730 (accessed on 15 December 2022).SMHI
UAE (Abu Dhabi) TWL (occupational) for workers Physiological   data   like   height ,   age ,   sweat   rate ,   T a ,   T w ,   T g https://weather.com/weather/today/l/Abu+Dhabi+Abu+Dhabi+Emirate+United+Arab+Emirates?placeId=0755f9b1a0f85388ca0d9510010eed3e6274c95ec9ecc1a8353af4782d304238 (accessed on 15 December 2022).The weather channels
Bangladesh DI T a , RHhttp://live.bmd.gov.bd/ (accessed on 15 December 2022).Bangladesh Meteorological department
Egypt DI (hot dry climate)https://www.weather-forecast.com/maps/Egypt (accessed on 15 December 2022).Weather Forecast
SudanCET T g ,   T w ,   R H ,   v https://worldweather.wmo.int/en/country.html?countryCode=203 (accessed on 15 December 2022).World Meteorological organization (M.D)
NigeriaCEThttps://www.nimet.gov.ng/ (accessed on 15 December 2022).NiMet (M.D)
Iran PET (hot dry climate) T a ,   v https://worldweather.wmo.int/en/country.html?countryCode=114 (accessed on 15 December 2022).World Meteorological Organization
CanadaHU (humid weather patterns) T a ,  dew point temperature, RH factor, molecular weight of water, latent heat and gas constanthttp://ec.gc.ca/meteo-weather/default.asp?lang=En (accessed on 11 June 2023). Weather and Meteorology. Retrieved 19 May 2016

References

  1. Sylla, M.B.; Faye, A.; Giorgi, F.; Diedhiou, A.; Kunstmann, H. Projected heat stress under 1.5 C and 2 C global warming scenarios creates unprecedented discomfort for humans in West Africa. Earth Future 2018, 6, 1029–1044. [Google Scholar]
  2. Cheung, S.S.; McLellan, T.M.; Tenaglia, S. The thermophysiology of uncompensable heat stress. Sports Med. 2000, 29, 329–359. [Google Scholar]
  3. Kovats, R.S.; Hajat, S. Heat stress and public health: A critical review. Annu. Rev. Public Health 2008, 29, 41–55. [Google Scholar]
  4. Smith, D.J.; Pius, L.M.; Plantinga, L.C.; Thompson, L.M.; Mac, V.; Hertzberg, V.S. Heat stress and kidney function in farmworkers in the US: A scoping review. J. Agromed. 2022, 27, 183–192. [Google Scholar]
  5. Krishnamurthy, M.; Ramalingam, P.; Perumal, K.; Kamalakannan, L.P.; Chinnadurai, J.; Shanmugam, R.; Srinivasan, K.; Venugopal, V. Occupational heat stress impacts on health and productivity in a steel industry in Southern India. Saf. Health Work. 2017, 8, 99–104. [Google Scholar] [PubMed]
  6. Yeo, T.P. Heat stroke: A comprehensive review. AACN Adv. Crit. Care 2004, 15, 280–293. [Google Scholar]
  7. Romaszko-Wojtowicz, A.; Cymes, I.; Dragańska, E.; Doboszyńska, A.; Romaszko, J.; Glińska-Lewczuk, K. Relationship between biometeorological factors and the number of hospitalizations due to asthma. Sci. Rep. 2020, 10, 9593. [Google Scholar] [PubMed]
  8. Epstein, Y.; Moran, D.S. Thermal comfort and the heat stress indices. Ind. Health 2006, 44, 388–398. [Google Scholar]
  9. Prek, M.; Butala, V. Comparison between Fanger’s thermal comfort model and human exergy loss. Energy 2017, 138, 228–237. [Google Scholar]
  10. Kim, J.; Schiavon, S.; Brager, G. Personal comfort models–A new paradigm in thermal comfort for occupant-centric environmental control. Build. Environ. 2018, 132, 114–124. [Google Scholar]
  11. Aubrecht, C.; Özceylan, D. Identification of heat risk patterns in the US National Capital Region by integrating heat stress and related vulnerability. Environ. Int. 2013, 56, 65–77. [Google Scholar] [CrossRef]
  12. Beshir, M.Y.; Ramsey, J.D. Heat stress indices: A review paper. Int. J. Ind. Ergon. 1988, 3, 89–102. [Google Scholar]
  13. Miller, V.S.; Bates, G.P. The thermal work limit is a simple reliable heat index for the protection of workers in thermally stressful environments. Ann. Occup. Hyg. 2007, 51, 553–561. [Google Scholar] [PubMed]
  14. Change, I.C. Impacts, adaptation and vulnerability. In Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2014; p. 1132. [Google Scholar]
  15. Meng, X.; Xue, S.; An, K.; Cao, Y. Physiological Indices and Subjective Thermal Perception of Heat Stress-Exposed Workers in an Industrial Plant. Sustainability 2022, 14, 5019. [Google Scholar] [CrossRef]
  16. Burton, I.; Ebi, K.L.; McGregor, G. (Eds.) Biometeorology for Adaptation to Climate Variability and Change. In Biometeorology for Adaptation to Climate Variability and Change; Springer: Dordrecht, The Netherlands, 2009; Volume 1, pp. 1–5. [Google Scholar]
  17. Kestrel 5400 User Manual PDF Download | ManualsLib. Available online: https://www.manualslib.com/manual/1399319/Kestrel-5400.html (accessed on 2 October 2022).
  18. Instruments de Mesure | Sauermann Group. Available online: https://sauermanngroup.com/fr-FR/instruments-de-mesure (accessed on 2 October 2022).
  19. LSI Heat Shield Heat Stress Monitors | Sensidyne. Available online: https://www.sensidyne.com/heatshield-heat-stress-monitor/ (accessed on 2 October 2022).
  20. 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]
  21. UTCI—Universal Thermal Climate Index. Available online: http://www.utci.org/utci_doku.php (accessed on 2 October 2022).
  22. 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] [PubMed]
  23. Błażejczyk, K. Mapping of UTCI in local scale (the case of Warsaw). Pr. I Stud. Geogr. WGSR UW 2011, 47, 275–283. [Google Scholar]
  24. Belding, H.S.; Hatch, T.F. Index for evaluating heat stress in terms of resulting physiological strains. ASHRAE Trans. 1955, 27, 129–136. [Google Scholar]
  25. An Outdoor Thermal Comfort Index (Out_SET*)–Part 1–The Model and Its Assumptions. Available online: https://www.researchgate.net/publication/268983057_An_outdoor_thermal_comfort_index_Out_SET-Part_1-The_model_and_its_assumptions (accessed on 3 October 2022).
  26. Moran, D.S.; Laor, A.; Heled, Y.; Pandolf, K.B.; Matthew, W.T.; Gonzalez, R.R. Evaluation and refinement of the environmental stress index for different climatic conditions. J. Basic Clin. Physiol. Pharmacol. 2003, 14, 1–16. [Google Scholar] [CrossRef] [PubMed]
  27. ISO 7730:2005; Ergonomics of the Thermal Environment—Analytical Determination and Interpretation of Thermal Comfort Using Calculation of the PMV and PPD Indices and Local Thermal Comfort Criteria. ISO: London, UK, 2005. Available online: https://www.iso.org/standard/39155.html (accessed on 3 October 2022).
  28. Fanger, P.O. Thermal Comfort. Analysis and Applications in Environmental Engineering; Danish Technical Press: Copenhagen, Denmark, 1970. [Google Scholar]
  29. Brake, R.; Bates, G. A valid method for comparing rational and empirical heat stress indices. Ann. Occup. Hyg. 2002, 46, 165–174. [Google Scholar]
  30. Giles, B.D.; Balafoutis, C.; Maheras, P. Too hot for comfort: The heatwaves in Greece in 1987 and 1988. Int. J. Biometeorol. 1990, 34, 98–104. [Google Scholar] [CrossRef]
  31. Sharma, M.R.; Ali, S. Tropical summer index—A study of thermal comfort of Indian subjects. Build. Environ. 1986, 21, 11–24. [Google Scholar] [CrossRef]
  32. Webb, C.G. Thermal Discomfort in an Equatorial Climate. A Nomogram for the Equatorial Comfort Index. J. Inst. Heat. Vent. Eng. 1960, 28, 297–304. [Google Scholar]
  33. Humidex; A Method of Quantifying Human Discomfort Due to Excessive Heat and Humidity | WorldCat.org. Available online: https://www.worldcat.org/title/humidex-a-method-of-quantifying-human-discomfort-due-to-excessive-heat-and-humidity/oclc/6952677 (accessed on 3 October 2022).
  34. Pepi, J.W. The Summer Simmer Index. Weatherwise 2010, 40, 143–145. [Google Scholar] [CrossRef]
  35. Lind, A.R.; Hellon, R.F. Assessment of physiological severity of hot climates. J. Appl. Physiol. 1957, 11, 35–40. [Google Scholar] [CrossRef] [PubMed]
  36. Nikolopoulou, M. Designing Open Spaces in the Urban Environment: A Bioclimatic Approach; Centre for Renewable Energy Sources: Pikermi Attiki, Greece, 2004. [Google Scholar]
  37. Australia’s Official Weather Forecasts & Weather Radar—Bureau of Meteorology. Available online: http://www.bom.gov.au/ (accessed on 3 October 2022).
  38. Heat Index—Wikipedia. Available online: https://en.wikipedia.org/wiki/Heat_index (accessed on 3 October 2022).
  39. Minard, D.; Belding, H.S.; Kingston, J.R. Prevention of heat casualties. J. Ame. Med. Assoc. 1957, 165, 1813–1818. [Google Scholar] [CrossRef]
  40. Rothfusz, L.P. The Heat Index Equation (Or, More than You Ever Wanted to Know About Heat Index); National Oceanic and Atmospheric Administration, National Weather Service, Office of Meteorology: Fort Worth, TX, USA, 1990; p. 640. [Google Scholar]
Figure 1. Schematic diagram of the adopted methodology: This figure presents a visual representation of the workflow used in this study for the sensitivity analysis of heat stress indices. It outlines the steps taken from data collection to analysis, providing a clear overview of the research process.
Figure 1. Schematic diagram of the adopted methodology: This figure presents a visual representation of the workflow used in this study for the sensitivity analysis of heat stress indices. It outlines the steps taken from data collection to analysis, providing a clear overview of the research process.
Climate 11 00181 g001
Figure 2. Global distribution of Outdoor Heat Stress Indices: This map illustrates the different indices officially used around the world. Each country is color-coded based on the specific index it uses, providing a global perspective on the prevalence and distribution of these indices.
Figure 2. Global distribution of Outdoor Heat Stress Indices: This map illustrates the different indices officially used around the world. Each country is color-coded based on the specific index it uses, providing a global perspective on the prevalence and distribution of these indices.
Climate 11 00181 g002
Figure 3. Sensitivity analysis of Simulated Heat Stress Indices: This visualization illustrates the responsiveness of each heat stress index to variations in environmental parameters. By comparing the indices, readers can discern the relative sensitivity of each, understanding which indices might be more robust or susceptible to specific changes in conditions. Such insights are crucial for policymakers and researchers aiming to select or design indices that are both accurate and resilient to varying climatic scenarios.
Figure 3. Sensitivity analysis of Simulated Heat Stress Indices: This visualization illustrates the responsiveness of each heat stress index to variations in environmental parameters. By comparing the indices, readers can discern the relative sensitivity of each, understanding which indices might be more robust or susceptible to specific changes in conditions. Such insights are crucial for policymakers and researchers aiming to select or design indices that are both accurate and resilient to varying climatic scenarios.
Climate 11 00181 g003
Table 1. Important factors involved in heat stress.
Table 1. Important factors involved in heat stress.
EnvironmentalPhysiologicalBehavioral
Air temperatureAgeActivity
Wind velocityGenderClothing
Relative humidityHealth status (body temperature, rate of heartbeat, blood pressure)Insulation
Mean radiant temperatureSweat ratePosture
Globe temperature
Location
Table 2. Outdoor heat stress indices and their methods for estimation.
Table 2. Outdoor heat stress indices and their methods for estimation.
Heat Stress IndicesFormulaNomenclature
WBGT: Wet bulb globe temperature [20]WBGT = 0.7 T w + 0.2 T g   + 0.1 T a   C r e s = heat exchange by convection in breathing
E c = heat exchange by evaporation on the skin
E r e s = evaporative heat exchange in breathing
E q = required evaporative heat
E m a x = maximum evaporative heat
H = sensitive heat losses
m e t = metabolic rate
T m r t = mean radiant temperature
m w = effective mechanical power
R H = relative humidity
S R = solar radiation
T w = wet-bulb temperature
T a = air temperature
T g = globe temperature
v = wind speed/air velocity
v p = vapor pressure
UTCI: Universal Thermal Climate Index [21,22,23]Calculated by 6th-order polynomial
U T C I = f ( T a   ; T m r t ; R H ; v ) U T C I = 3.21 + 0.872 T a   + 0.2459 T m r t 2.5078 v 0.0176 R H
HSI: Heat stress index [24] H S I = 100 E q E m a x
Out_SET: Outdoor standard effective temperature [25] O U T S E T = ( W B G T 11.76 ) / 0.405
ESI: Environmental Stress Index [26] E S I = 0.62 T a   0.007 R H + 0.002 . S R + 0.0043 ( T a   R H ) 0.078 0.1 + S R 1
PMV: Predicted Mean Vote [27,28] P M V = 0.303 e 2.1 × m e t + 0.028 × m e t m w H E c C r e s E r e s α
CET: Corrected Effective Temperature [29]CET= 1.21 T g 0.21 T w / 1 + 0.029 T g T w )
DI: Discomfort Index [30] D I = T a 0.55 ( 1 0.01 R H ) T a 14.5
TSI: Tropical summer index [31] T S I = 1 / 3 T w + 3 / 4 T g 2 v 1 / 2
ECI: Equatorial comfort Index [32] E C I = 0.574 T a + 0.271 v p 1.81 v 1 / 2 + 4.2
HU: Humidex [33] H U = T a   + ( 5 / 9 ) ( e 10 ) e = 6.11 × e 5417.7530 × 1 273.16 1   Dewpoint   in   Kelvin  
SSI: Summer simmer index [34] S S I = 1.98 T a ( 0.55 0.0055 . R H ) T a 58 56.83
OI: Oxford index [35]OI = 0.85 T w + 0.15 T a
ASV: Actual Sensation Vote value for Europe [36] A S V = 0.049 T a + 0.001 . S R 0.051 v + 0.014 . R H 2.079
AT: Australian apparent temperature [37] A T = T a + 0.33 v p 0.70 v 4.00
HI: Heat index [38] H I = 42.379 + 2.05 . T a + 10.14 . R H 0.22 . T a . R H 6.84 × 10 3 . T a 2 5.48 × 10 2 . R H 2 + 1.23 × 10 3 . T a 2 . R H + 8.5 × 10 4 . T a . R H 2 1.99 × 10 6 . T a 2 . R H ²
Important Remark: Each heat index has different assumptions/calibrations (such as body size, physical fitness, etc.) that consider temperature and humidity differently. A high-heat-stress event indicated by one index does not necessarily transfer to another index. For example, the original equation for WBGT was derived and calibrated using US Marine Corps Marines during basic training [39] and in good physical condition. HI was calibrated for an “average” American male [40]. We bring the attention of the reader only to the indices which can be measured using sensors, i.e., CET, HSI, HI, WBGT, and Thermal Work Limit (TWL). All other indices presented in Table 2 are also estimated using their respective formulas or physical models that depend on assumptions e.g., age group, activities, clothing, area, etc., but cannot be directly measured.
Table 3. Sensitivity analysis of heat stress indices.
Table 3. Sensitivity analysis of heat stress indices.
Heat Stress IndicesS0Sensitivity around S0Partial Differential
CET C E T o = 27.83 °C
T g o = 45 °C
T ω o = 16.27 °C
C E T T g 0 = 0.19
C E T T w 0 = −0.08
C E T T g 0 = 1.21 1 + 0.029 T g 0 T ω 0 1.21 × 0.029 × T g 0 1 + 0.029 T g 0 T ω 0 2 C E T T w 0 = 0.21 1 + 0.029 T g 0 T ω 0 + 0.21 × T ω 0 × 0.029 1 + 0.029 T g 0 T ω 0 2
PMV P M V o = 0.379
m e t 0 = 1.68
m w 0 = 0.1
α = −0.33
P M V m e t 0 = 3.6 × 10 3 P M V m w 0 = −3.6 × 10 2

P M V α 0 = −0.33
P M V m 0 = β   w h e r e   β = 2.1 × 0.303 × e 2.1 m 0 × m 0 + 0.303 × e 2.1 m 0 + 0.028 2.1 × 0.303 × e 2.1 m 0 × ω 0 P M V ω 0 = 0.028 0.303 × e 2.1 m 0 P M V α 0 = 1
HI H I o = 78.74 °F
T a o = 81 °F
R H 0 = 30%
H I T a o 0 = 1.71
H I R H 0 = 0.32
H I T 0 = 2.049 0.22 R H 0 2 6.84 × 10 3 T a o + 2 1.22 × 10 3 T a o ( R H 0 ) + 8.5 × 10 4 ( R H 0 ) ² 1.99 × 10 6 T a o ( R H 0 ) ² H I R H 0 = 10.143 0.22 ( T a o ) 2 ( 5.48 × 10 2 ) ( R H 0 ) + 1.22 × 10 3 T a o ² + 2 ( 8.5 × 10 4 ) T a o ( R H 0 )   2 ( 1.99 × 10 6 ) T a o ² ( R H 0 )
SSI S S I o = 78.52 °F
T a o = 81 °F
R H 0 = 30%
S S I T a o 0 = 0.89
S S I R H 0 = 0.25
S S I T 0 = 1.98 0.45 + 0.0055 R H 0 S S I R H 0 = 1.98 [ 0.0055 T a o 0.0055 58 ]
DI D I o = 20.24 °C
T a o = 27.22 °C
R H 0 = 30%
D I T a o 0 = 0.452
D I R H 0 = 0.07
D I T 0 = 0.45 + 0.0055 R H 0 D I R H 0 = 0.0055 T a o 0.07975
WBGT W B G T o = 23.11 °C
T g o = 45 °C
T ω o = 16.27 °C
T a o = 27.22 °C
W B G T T g 0 = 0.2
W B G T T w 0 = 0.7
W B G T T a o 0 = 0.1
WBGT = 0.7 T w + 0.2 T g   + 0.1 T a  
UTCI U T C I 0 = 28.52 °C
T a o = 27.22 °C
T m r t = 27 °C
v o = 2 m/s
R H 0 = 30%
U T C I T a o 0 = 0.87
U T C I R H 0 = −1.76 × 10 2 U T C I T m r t 0 = 0.24
U T C I v 0 = −2.5
U T C I = 3.21 + 0.872 T a   + 0.2459 T m r t 2.5078 v 0.0176 R H
Table 4. Reference comfort scale for widely used Thermal Stress Indices.
Table 4. Reference comfort scale for widely used Thermal Stress Indices.
Comfort Zone (Z)WBGT
°C
CET
°C
UTCI
°C
PMVDI
°C
SSI
°F
HI
°F
extreme cold stress (Z-5). <−40−3–−2.5
very strong cold stress (Z-4) −35–−25
strong cold stress (Z-3) −25–−13−2.5–−1.5
moderate cold (Z-2)<17<17−13–0−1.5–−0.5
slightly cold (Z-1) 0–8−0.5–0
comfortable (Z0) 17–308–250–0.5 77–83
slightly warm (Z1)17–2330–3425–310.5–1.522–2483–91
moderate warm (Z2)23–2834–4531–371.5–2.524–2891–10080–91
strong heat stress (Z3)28–3045–4937–462.5–3>28100–11290–105
very hot (Z4)>30 46 to above3 112–125105–130
extremely dangerous (Z5) >49 >3 125–150>130
Table 5. Sensitivity of comfort level—this table illustrates the variations in comfort levels across different heat indices, highlighting the differential responses to environmental parameters.
Table 5. Sensitivity of comfort level—this table illustrates the variations in comfort levels across different heat indices, highlighting the differential responses to environmental parameters.
Index y Δ y
HI 78.7480.78
SSI 78.5279.67
PMV 0.370.29
WBGT 23.1123.11
UTCI 28.5228.52
CET 27.8327.94
DI 20.2420.76
Table 6. Comparative analysis of simulated indices.
Table 6. Comparative analysis of simulated indices.
Heat Stress IndicesEffectiveness AdvantageLimitationRecommendation
WBGTEffective with direct solar radiation and high humidity.It is a comprehensive measure for heat stress that can be calibrated directly by using sensors.It does not account for wind speed and clothing insulation, which can affect the perceived temperature.It is recommended for outdoor activities.
CETEffective in humid regions, making it suitable for tropical climates.Incorporates modern adjustments to older indices, improving calibration accuracy.Underestimates outdoor heat stress with varying conditions.Best suited for indoor environments.
UTCIEffective in a variety of climates and conditions from cold to hot and dry to humid.It integrates air temperature, wind speed, humidity, and solar radiation.Requires a complex calculation.Suitable for urban planning and public health advisories.
PMVPMV is effective in controlled indoor environments.It provides a complete scale to indicate the level of human comfort.Requires detailed input data for accurate prediction, making calibration complex.Ideal for building design and HVAC system optimization.
DIEffective in various climates.Simple to calculate and can be measured directly by using sensors.Primarily considers temperature and humidity, which might not be sufficient in extreme conditions.
Limited scale for indicating thermal discomfort.
Suitable for general weather forecasts and public advisories.
SSIEffective in hot and high-humidity conditions.Defined thermal comfort scale for hot conditions.Limited to summer season.Useful for public health advisories during summer months.
HIAdaptable to a range of hot conditions, especially in high humidity.Simple calculations and can be measured by sensors.Sensitive and might not be effective in dry or extremely cold conditions.Useful for high-humidity areas.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rachid, A.; Qureshi, A.M. Sensitivity Analysis of Heat Stress Indices. Climate 2023, 11, 181. https://doi.org/10.3390/cli11090181

AMA Style

Rachid A, Qureshi AM. Sensitivity Analysis of Heat Stress Indices. Climate. 2023; 11(9):181. https://doi.org/10.3390/cli11090181

Chicago/Turabian Style

Rachid, Ahmed, and Aiman Mazhar Qureshi. 2023. "Sensitivity Analysis of Heat Stress Indices" Climate 11, no. 9: 181. https://doi.org/10.3390/cli11090181

APA Style

Rachid, A., & Qureshi, A. M. (2023). Sensitivity Analysis of Heat Stress Indices. Climate, 11(9), 181. https://doi.org/10.3390/cli11090181

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