Next Article in Journal
3D Chromatin Architecture Re-Wiring at the CDH3/CDH1 Loci Contributes to E-Cadherin to P-Cadherin Expression Switch in Gastric Cancer
Previous Article in Journal
Reduced Antibiotic Resistance in the Rhizosphere of Lupinus albus in Mercury-Contaminated Soil Mediated by the Addition of PGPB
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 12
Issue 6
10.3390/biology12060802
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Temperature–Humidity-Dependent Wind Effects on Physiological Heat Strain of Moderately Exercising Individuals Reproduced by the Universal Thermal Climate Index (UTCI)

by
Peter Bröde
1,* and
Bernhard Kampmann
2
1
Leibniz Research Centre for Working Environment and Human Factors at TU Dortmund (IfADo), Ardeystrasse 67, D-44139 Dortmund, Germany
2
Department of Occupational Health Science, School of Mechanical Engineering and Safety Engineering, University of Wuppertal, Gaußstraße 20, D-42119 Wuppertal, Germany
*
Author to whom correspondence should be addressed.
Biology 2023, 12(6), 802; https://doi.org/10.3390/biology12060802
Submission received: 28 April 2023 / Revised: 29 May 2023 / Accepted: 31 May 2023 / Published: 31 May 2023
(This article belongs to the Section Physiology)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

With climate change exacerbating heat extremes, there is a growing need for sustainable measures that can reduce physiological heat strain. However, public health policies have issued warnings against the use of electric fans and ventilators at ambient temperatures exceeding typical skin temperatures of 35 °C. Aiming to extend previous findings for sedentary persons to moderately exercising individuals, we analyzed the heart rates, sweat rates, and core and skin temperatures recorded in 198 climatic chamber experiments with 3h exposures to treadmill work under widely varying heat stress conditions. Our findings suggest mitigating wind effects at temperatures above 35 °C when humidity is elevated. Moreover, the close agreement of the observed effects with the corresponding wind effects predicted by the Universal Thermal Climate Index (UTCI) demonstrates the potential of the UTCI for evaluating sustainable strategies of heat stress mitigation for moderately exercising individuals.

Abstract

Increasing wind speed alleviates physiological heat strain; however, health policies have advised against using ventilators or fans under heat wave conditions with air temperatures above the typical skin temperature of 35 °C. Recent research, mostly with sedentary participants, suggests mitigating the effects of wind at even higher temperatures, depending on the humidity level. Our study aimed at exploring and quantifying whether such results are transferable to moderate exercise levels, and whether the Universal Thermal Climate Index (UTCI) reproduces those effects. We measured heart rates, core and skin temperatures, and sweat rates in 198 laboratory experiments completed by five young, semi-nude, heat-acclimated, moderately exercising males walking the treadmill at 4 km/h on the level for three hours under widely varying temperature–humidity combinations and two wind conditions. We quantified the cooling effect of increasing the wind speed from 0.3 to 2 m/s by fitting generalized additive models predicting the physiological heat stress responses depending on ambient temperature, humidity, and wind speed. We then compared the observed wind effects to the assessment performed by the UTCI. Increasing the wind speed lowered the physiological heat strain for air temperatures below 35 °C, but also for higher temperatures with humidity levels above 2 kPa water vapor pressure concerning heart rate and core temperature, and 3 kPa concerning skin temperature and sweat rate, respectively. The UTCI assessment of wind effects correlated positively with the observed changes in physiological responses, showing the closest agreement (r = 0.9) for skin temperature and sweat rate, where wind is known for elevating the relevant convective and evaporative heat transfer. These results demonstrate the potential of the UTCI for adequately assessing sustainable strategies for heat stress mitigation involving fans or ventilators, depending on temperature and humidity, for moderately exercising individuals.

Graphical Abstract

1. Introduction

1.1. Sustainable Heat Stress Mitigation by Wind

Climate change accelerating with fossil energy use will exacerbate heat extremes, which will provoke physiological heat strain in terms of increased heart rates, core and skin temperatures, and sweat rates, imposing health risks [1,2], and potentially affecting human performance and productivity [3,4]. Elevating air movements by outdoor wind or ventilators is a cost-effective and sustainable measure for mitigating the physiological heat strain that is predominantly considered indoors and at stationary workplaces [5,6,7]. However, policies and public health authorities have advised against electric fan use under heat wave conditions indoors with air temperatures above the typical skin temperature of 35 °C, because convective cooling will then turn to convective heating of the body [8,9,10]. On the other hand, the well-known enhancement of sweat evaporation with increased wind speed [11,12,13] challenged these postulations and advocated for an evidence-based approach [14]. Recent studies, mostly with sedentary participants [15,16], have demonstrated the cooling effects of ventilation due to enhanced sweat evaporation at even higher temperatures at elevated humidity levels.
Comparative studies involving experimental heat exposures are usually restricted to a few well-defined climatic conditions [15]. Thus, for separating beneficial cooling from detrimental heating wind effects over a grid of temperature–humidity conditions in heat wave scenarios, several studies applied biophysical modeling based on human heat balance calculation [17,18,19,20]. Interestingly, qualitatively similar threshold curves based on modeling the naturally aspirated wet-bulb temperature were developed one century ago [21]. These curves indicated wind cooling at elevated humidity levels but additional heat loads in hot–arid climates, as flagged by the term ‘poison wind’ (simoom) for potentially fatal outdoor conditions in Southwest Asian desert regions [21]. Notably, the simulations of heat wave scenarios focused on indoor settings and only considered resting conditions concerning metabolic rates of about 1 MET (1 MET = 58.2 W/m2), whereas physical loads at workplaces and many home and outdoor activities are associated with metabolic rates higher than 2 MET [22,23].
Concerning physiological heat strain under moderate exercise levels, the recorded heart rates, core and skin temperatures, and sweat rates in Figure 1 from an extensive database of climate chamber experiments [24,25] exemplify the aforementioned findings for an acclimated young male walking at 4 km/h on a treadmill, corresponding to a 2.3 MET metabolic rate. While increasing the air velocity from 0.3 to 2 m/s lowered the physiological heat strain in a hot–humid climate with air temperature exceeding skin temperature (Figure 1a), increased wind amplified the heat strain under hot–dry conditions (Figure 1b). Figure A1 in Appendix A includes a supporting example from another participant. In both cases, the transient overshooting of the sweat rate with low wind (va = 0.3 m/s) under warm–humid conditions (Figure 1a) may be explained by hidromeiosis [26], which did not occur when increasing the wind speed to va = 2.0 m/s enhanced sweat evaporation [27]. On the other hand, under hot–dry conditions (Figure 1b), the already high level of evaporative efficiency could not be further elevated by increased airflow and thus could not compensate for the higher convective heat gain aggravating physiological strain, in particular rising sweat rates [28].
A recent study with exercising participants [28] examined how wind speed at various temperature–humidity combinations affected physical work capacity, defined by the workload related to metabolic rate, which the participants could tolerate with their heart rate clamped at 130 beats per minute (bpm). Coupling their experimental data with the heat balance calculations, they showed beneficial wind effects for air temperatures up to 44 °C with relative humidity exceeding 50% but detrimental wind effects for higher temperatures. Similar findings have been reported concerning the shift of the critical, i.e., tolerable, humidity level by wind in hot environments [29]. Another study [30] showed beneficial wind effects for sports activities at specific conditions at 37 °C air temperature and 50% relative humidity.
In spite of the aforementioned research, and although Figure 1 and Figure A1 exhibit motivating examples, more systematic studies are required for quantification and evidence-based conclusions. Specifically, with a focus on moderately active persons, there is a lack of data from controlled laboratory experiments specifying wind effects on physiological heat strain in terms of heart rates, core and skin temperatures, and sweat rates while covering a wide grid of air temperature–humidity combinations, which determine the indoor thermal environment.
In addition to data requirements, there seems to be a need for improving the thermo-physiological modeling beyond the heat balance approach, as recent findings suggest that the simple models will underestimate physiological heat strain under electrical fan use, e.g., for the vulnerable elderly population [31].
When looking for advanced, but easy-to-use modeling approaches, the Universal Thermal Climate Index (UTCI) [32], developed for moderately active persons (2.3 MET), constitutes a noteworthy alternative. The index was constructed from the dynamic physiological thermal strain simulated by the advanced UTCI-Fiala model of human thermoregulation [33] coupled with an adaptive clothing model [34]. The UTCI allows for the assessment of the thermal environment, covering the range from extreme cold to extreme heat stress conditions. The model was extensively validated against laboratory data and field observations [35] and showed good agreement with thermal environment standards and experimental data in occupational settings [36,37,38] concerning the moderating influence of humidity, wind, and radiant heat load in addition to air temperature on both cold and heat stress. Though the UTCI relies on sophisticated models, the operational procedure [39] provides algorithms, making it easily applicable, e.g., for assessing the wind-cooling potential of actual and future urban scenarios in relation to wind direction [40]. Concerning the effects on physical work capacity, the UTCI outperformed any other considered thermal index in capturing not only the effects of wind over a wide grid of temperature–humidity conditions [28], but also in combination with heat radiation [41].

1.2. Study Objectives

Based on a comprehensive physiological heat strain database [24], this study aimed to explore and quantify whether reports on beneficial wind effects with air temperature exceeding skin temperature under humid conditions for sedentary individuals are transferable to moderate exercise levels, and whether the Universal Thermal Climate Index (UTCI) reproduces these effects.
Our approach was to determine the transition from cooling to heating wind effects on heart rates, sweat rates, and core and skin temperatures measured in experiments with exercising participants over a comprehensive matrix of temperature–humidity combinations and to compare the outcomes to the UTCI assessment.

2. Materials and Methods

2.1. Experimental Data

This experimental study is based on a database compiled from about 2400 climate chamber experiments with 33 young male adults, conducted previously at IfADo [24,42], under widely varying climatic conditions, clothing, and activity levels. The experiments were carried out according to the ethical principles of the Declaration of Helsinki after approval by IfADo’s local Ethics Committee.
We searched our database for series of heat stress experiments with participants exercising on a treadmill under both reference (RefWind, with air velocity va = 0.3 m/s) and high wind (HiWind, va = 2 m/s) conditions. Here, 2 m/s constituted the maximum value of controllable wind speed in the climatic chamber, corresponding to average airflow conditions observable in (semi-)outdoor workplaces and in coal mines [24]. Inclusion criteria were a minimum number of 15 experiments per series with comparable workloads and clothing. We retrieved 198 trials organized in 10 series, which originated from five acclimated, semi-nude (basic clothing insulation Icl = 0.1 clo, 1 clo = 0.155 K·m2·W−1), young, and fit males under RefWind and HiWind conditions, respectively. The number of experiments in each series varied inter-individually and depended on wind conditions between 16 and 25 experiments, with total numbers of 97 trials for RefWind and 101 for HiWind. The average personal characteristics (mean ± SD, with range in brackets) of the five participants were 20.1 ± 0.9 (19–22) years of age, 1.87 ± 0.02 (1.84–1.88) m of body height, 70.5 ± 2.1 (68–73) kg of body weight, 1.94 ± 0.02 (1.92–1.97) m2 of body surface area, and 47.9 ± 6.4 (43.2–57.4) mL/min/kg of maximal oxygen consumption. Before exposure, the participants had undergone a heat acclimation protocol lasting at least three weeks [43].
We will briefly summarize the procedures as detailed descriptions are available elsewhere [24]. Each trial consisted of treadmill work with constant workload of walking 4 km/h on the level, which corresponded to the activity level of 2.3 MET, as assumed for UTCI [39]. Each trial lasted for at least three hours and consisted of 30 min work periods interrupted by 3 min breaks for determining body weight loss, from which sweat rate (SR) was calculated correcting for ad libitum fluid ingestion. Under both HiWind and RefWind, the participants were exposed to varying levels of heat stress, with conditions characterized by different combinations of air temperature (Ta; range 25–55 °C) and humidity, expressed as water vapor pressure (pa; 0.5–5.3 kPa). Mean radiant temperature (Tmrt) was equal to Ta. Trials were stopped prematurely if rectal temperature exceeded 38.5 °C or on the participant’s demand.
Rectal temperatures (Tre) were recorded continuously using a thermistor probe (YSI 401, YSI Inc., Yellow Springs, OH, USA) inserted 10 cm past the anal sphincter. Local skin temperatures were measured with thermistors (YSI 427, YSI Inc., Yellow Springs, OH, USA) at the forehead (Tsk,head), chest (Tsk,chest), back (Tsk,back), upper arm (Tsk,arm), thigh (Tsk,thigh), and lower leg (Tsk,leg), and were used to calculate mean skin temperature (Tsk) as weighted average of the local skin temperatures according to [24], shown in Equation (1):
Tsk = 0.05 × Tsk,head + 0.2 × Tsk,chest + 0.15 × Tsk,back + 0.2 × Tsk,arm + 0.25 × Tsk,thigh + 0.15 × Tsk,leg
Heart rates (HR) were obtained using ECG electrodes and stored in one-minute intervals, as were Tre and Tsk. As illustrated by the yellow shaded areas in Figure 1, the averages of Tre, HR, Tsk, and SR over the third hour of exposure, representing steady-state [25], were submitted to the following data analysis.

2.2. Data Analysis and Statistics

Penalized regression splines fitted by generalized additive models (GAMs) provide a flexible modeling framework for irregularly distributed data [44], allowing for non-linear [45] and random effects [46].
E[Yij] = µ + s(IDi) + te(Ta,ij, pa,ij) + teΔv(Ta,ij, pa,ij) + εij
Using the notation from Equation (2) [44], we predicted HR, Tre, Tsk, and SR, respectively, as responses from experiment j of participant i (Yij) by separate GAMs. These models included an overall intercept µ, and a bivariate penalized tensor regression spline te(Ta,ij, pa,ij) for the effects of air temperature and humidity depending on the (Ta, pa) combinations, accounting for the repeated measurements by subject-specific intercepts as random coefficients s(IDi) [46]. Adding another bivariate spline teΔv(Ta,ij, pa,ij) as so-called factor smooth interaction [44] of the regression splines with the wind condition (RefWind vs. HiWind), we obtained estimates for the wind effect Δv, i.e., the response difference under HiWind compared to RefWind over the (Ta, pa) grid supplemented by p-values [47]. These were then visualized in a difference plot [48].
More specifically, s(IDi) refers to a cubic regression spline with basis dimension (or rank) set to k = 5 (i.e., with k − 1 = 4 as upper limit of the associated degrees of freedom). Similarly, te(Ta,ij, pa,ij) and teΔv(Ta,ij, pa,ij) represent tensor products of two cubic regression spline bases for the Ta and pa dimension, respectively, each of rank k = 9, and are hence of total rank k = 81. Maximum likelihood parameter estimates with standard errors and p-values were obtained assuming Gaussian error ϵij [47]. Calculations were performed with the R software version 4.2.1 [49] using the package mgcv [44] together with mgcViz [48] and tidymv [50].

2.3. UTCI Calculations

We calculated UTCI values using the regression polynomial, as described in the UTCI operational procedure [39] for given combinations of air temperature (Ta), ambient water vapor pressure (pa), mean radiant temperature (Tmrt), and air velocity 10 m above ground (va,10m).

2.3.1. UTCI Sensitivity to Wind

Following meteorological conventions, UTCI calculations rely on air velocity 10 m above ground (va,10m). For conversion to any other measurement height, e.g., 1 m for person level, the operational procedure [39] provides a logarithmic formula, shown in Equation (3), indicating that air velocity at person level (va,1m) is computed as the 10 m value (va,10m) divided by 1.5, in accordance with international standards [51].
va,1m = va,10m × log(1/0.01)/log(10/0.01) = va,10m/1.5
As the UTCI person is assumed to move on the level at 4 km/h, corresponding to a walking speed of vw = 1.1 m/s, this is taken into account by calculating the resulting or relative air velocity at person level (var,1m), according to Equation (4). Here, α denotes the angle between the directions of walking and wind assigned to zero for indicating the same direction. As UTCI does not assume a specific angle, var,1m is calculated by integrating Equation (4) over all α between zero and 2π [34,52].
v a r , 1 m = v w − v a , 1 m × c o s α 2 + v a , 1 m × s i n α 2
Figure 2 presents the resulting var,1m used by UTCI for calculating the convective and evaporative heat loss [33]. Notably, reducing wind speed below va,10m = 0.5 m/s (va,1m = 0.3 m/s) will hardly impact var,1m, which is limited by vw, while var,1m will approach va,1m for va,10m above 3 m/s.

2.3.2. Wind Effect on UTCI

We computed the effect of actual wind speed va,10m on UTCI (ΔvUTCI) compared to the reference wind speed (va,10m,ref = 0.5 m/s) while keeping the other parameters constant, as shown in Equation (5):
∆vUTCI = UTCI(Ta;pa;va,10m;Tmrt) − UTCI(Ta;pa;va,10m,ref;Tmrt)
Matching the range of experimental climatic conditions (Section 2.1), but respecting the upper limits of UTCI validity [39,53], we computed UTCI and ΔvUTCI for combinations of Ta between 25 and 50 °C and of pa from 0.1 to 5 kPa with relative humidity of rH ≤ 100% while setting Tmrt = Ta, as in the experiments. We performed the calculations for the reference wind speed (va,10m,ref = 0.5 m/s), matching RefWind with va,1m = 0.3 m/s (Figure 2), and increased wind speeds with va,10m = 3 m/s, matching HiWind with va,1m = 2 m/s.
Supplemental calculations were performed for assessing the modifying effect of radiant heat load by increasing ΔTmrt = Tmrt−Ta from 0 to 30 K in steps of 10 K. In addition, we considered conditions with higher wind speeds of va,10m = 4 and 6 m/s, respectively. Here, va,10m = 4 m/s corresponds to an increase of 1.7 m/s in relative air velocity at body level (var,1m), according to Figure 2, matching the settings of the climate chamber experiments, whereas va,10m = 6 m/s (va,1m = 4 m/s) parallels the wind speed applied in several simulations and experimental studies [15,17,18,19,20,28].

2.3.3. UTCI Prediction of Physiological Wind Effects

We compared the UTCI assessment of the wind effects ΔvUTCI with the wind effects obtained for the physiological variables by scatterplots, including Spearman correlation coefficients (rs). In addition, we quantified the predictive performance of UTCI in the binary classification of the occasion of heating wind effects (Δv > 0) on physiological variables by calculating the sensitivity, specificity, and overall accuracy.

3. Results

In order to illustrate the structure of our data analysis approach, and deviating from the order of topics in the methods section with the experimental study preceding UTCI calculations, we first present in Section 3.1 the wind effects (HiWind vs. RefWind) evaluated by the UTCI. Section 3.2 contains the corresponding analyses for the experimental heat strain data, while Section 3.3 relates the physiological wind effects to the UTCI assessment. Eventually, Section 3.4 shows the results for the supplemental UTCI analyses considering higher wind speeds and the influence of thermal radiation.

3.1. Wind Effects Assessed by UTCI

The psychrometric charts in Figure 3 display the UTCI contours for selected index values, representing the upper limits of the UTCI heat stress categories in relation to the air temperature and humidity for the different wind speeds (Figure 3a). The contours were bended upward to the left, indicating that the heat stress increased with temperature and humidity, as expected [36,39]. The dashed lines (HiWind) were above the solid lines (RefWind), indicating wind cooling at temperatures below 34–35 °C and in hot–humid conditions. On the other hand, dashed lines below solid lines indicating heating by wind were found in hot–dry conditions. The dot-dashed line connecting the intersections of the RefWind and HiWind contours represents the zero wind effect ΔvUTCI = 0.
Figure 3b quantifies the magnitude of the wind effect ΔvUTCI depending on the temperature and humidity. Wind cooling (blueish) increased with a decreasing temperature and increasing humidity, while heating due to wind (reddish) increased with an increasing temperature and decreasing humidity.

3.2. Wind Effects on Physiological Heat Strain

Table 1 summarizes the results of the GAMs fitted separately to the physiological heat strain variables. While the data for Tre were complete, there were a few missing values for HR, Tsk, and SR. The estimated mean values of 103 bpm (HR), 37.6 °C (Tre), 35.5 °C (Tsk), and 744 g/h (SR), respectively, were typical for light to moderate activities under heat stress [24,38]. The goodness-of-fit was good to excellent, with more than three quarters of the variance explained (R2 > 75%) and small residual standard errors for HR, Tre, and Tsk. SR showed almost a 90% explained variance but a slightly increased standard error, which, however, was still below 125 g/h, considered the relevant accuracy limit acceptable in occupational or military settings [54]. Moreover, graphical model checking, as displayed in Appendix A in Figure A2, did not reveal any problematic issues concerning the underlying model assumptions.
For all heat strain variables, the splines te(Ta, pa) modeling the influence of temperature and humidity were statistically highly significant (p < 0.0001). Similarly, the moderating effects of wind speed, as considered by the bivariate interaction splines teΔv(Ta, pa), were statistically significant for HR, Tre, Tsk, and SR.
Inter-individual variability, as indicated by the random term s(ID), played a considerable role (p < 0.0001) for Tre, Tsk, and SR, respectively, whereas individual influences on HR appeared to be less pronounced (p = 0.20).
Analogous to the UTCI results from Figure 3a, Figure 4 displays the effect of temperature and humidity separately for RefWind and HiWind, as contour lines for the physiological heat strain responses of HR, Tre, Tsk, and SR. Although the shape of the lines varied between the four variables and additionally depended on the strain level, a general pattern emerged, with contours bended upward to the left indicating strain levels increasing with temperature and humidity, as had also occurred for UTCI (Figure 3a). Again, similar to the UTCI, the dashed lines (HiWind) were above the solid lines (RefWind) at low temperatures and under hot–humid conditions, demonstrating reduced strain due to wind cooling, whereas in hot–dry climates, the opposite was found, with solid lines above dashed lines indicating higher strain levels with increased wind speed.
The dot-dashed lines represent the zero wind effect Δv = 0, thus demonstrating that the threshold separating cooling from heating effects was bended upward to the right, resemblingthe UTCI (Figure 3b), but with a greater curvature. They indicated reduced heat strain with HiWind at low temperatures and in hot–humid conditions, with limiting vapor pressure above 2 kPa for HR and Tre and above 3 kPa for Tsk and SR, respectively.
The difference plots in Figure 5 show the intensity of the cooling (blueish) and heating (reddish) wind effects over the temperature–humidity grid estimated by GAM, where white areas indicate indifferent temperature–humidity combinations with statistically non-significant differences from Δv = 0.
Whereas no significant core temperature increase due to wind was observed, detrimental wind effects with increased HR, Tsk, and SR were found under hot–arid conditions with a temperature above 40 °C and vapor pressure below 1 kPa. Temperatures below 35 °C yielded no detrimental wind effects.
Though significant decreases in the core temperature and heart rates by wind under high temperatures started at vapor pressures above 2 kPa, pronounced effects with reductions of more than 10 bpm HR or 0.15 °C Tre only appeared for very humid conditions with vapor pressures approaching 4 kPa.

3.3. UTCI Assessment Related to Physiological Wind Effects

The scatterplots in Figure 6a comparing the wind effects (Δv) on the UTCI with the physiological effects show positive correlations with the closest agreement (r ≈ 0.9) for skin temperature and sweat rate, also concerning the transition from cooling (Δv < 0) to heating (Δv > 0) effects. Slightly lower correlations were observed with the effects on heart rates and core temperature.
In addition to correlation analyses, the plots give an indication of the capacity of the UTCI for correctly classifying conditions with physiological wind cooling or heating [56]. The dashed reference lines in Figure 6a divide the plot into quadrants, with the lower left quadrant showing the physiological cooling events correctly assessed by the UTCI and the upper right containing the correctly assessed heating events, which dominate concerning Tsk and SR.
Whereas hardly any data points were present in the upper left quadrant with wind cooling falsely predicted by the UTCI, the conditions collected in the lower right quadrants (falsely predicted heating) concur with a more conservative UTCI assessment of the transition from cooling to heating wind concerning HR and Tre.
When quantifying the UTCI performance for predicting the occasion of heating wind effects (Δv > 0) as sensitivity, specificity, and accuracy in Figure 6b, these falsely predicted heating events resulted in lowered specificity and overall accuracy regarding HR and Tre, whereas well-balanced figures emerged for Tsk and SR.

3.4. UTCI Assessment with Higher Wind Speeds and Thermal Radiation

Figure A3 in Appendix A visualizes the wind effects (Δv) on the UTCI over the temperature–humidity grid for the different wind speeds and radiant heat loads, expressed as ΔTmrt = Tmrt – Ta, with the upper left panel with va,10m = 3 m/s and ΔTmrt = 0 K repeating the data from Figure 3b. Figure 7 summarizes the wind effects ΔvUTCI depending on the wind speed and radiant heat load as boxplots calculated over the temperature–humidity grid.
As expected, increasing the wind speed amplified both the cooling (ΔvUTCI < 0) and heating effects (ΔvUTCI > 0) at the extremes, as indicated by the minima and maxima, but also the first and third quartiles in Figure 7. In addition, the threshold line for ΔvUTCI = 0 slightly shifted toward higher temperatures for higher wind speeds (Figure A3), thus contributing to the decreasing trend with wind speed observed for the median effect in Figure 7.
Concerning the influence of heat radiation on ΔvUTCI, all percentiles showed a decreasing trend with an increasing radiant heat load ΔTmrt (Figure 7). The reduced reddish surface areas with an increased radiant heat load in Figure A3 suggest that this was attributable to a massive shift in the threshold line toward higher temperatures. Thus, wind-cooling effects occurred for a higher portion of climatic conditions over the temperature–humidity grid, shifting the ΔvUTCI distribution toward more negative values (Figure 7).
The wind effect ΔvUTCI at va,10m = 4 m/s (with ΔTmrt = 0 K) corresponded to the effect of increasing relative air velocity (var,1m) by 1.7 m/s (Figure 2), conforming to the change in the (relative) air velocity realized in the experiments. As the relevant UTCI calculations of convective and evaporative heat losses depend on the relative air velocity [33], we repeated the correlational analyses from Figure 6 for this condition, as shown in Figure A4. The outcome was very similar as before in Figure 6, though we observed slightly lowered correlations of the physiological wind effects with ΔvUTCI.

4. Discussion

This study achieved its primary aim to provide heat stress thresholds over a wide matrix of temperature–humidity combinations separating cooling from heating wind effects based on experimental evidence for moderately exercising individuals.

4.1. Temperature-Humdity-Dependent Wind Effect Thresholds

With a view on sustainable heat stress mitigation for moderate activity levels occurring during domestic or industrial work, we assessed the wind effects on physiological heat strain parameters connected with the use of, e.g., electric fans or ventilators. Linking a comprehensive database originating from climate chamber experiments [24] with flexible spline regression models [44], we quantified the physiological heat strain responses of moderately exercising individuals to increasing wind speed over a wide range of temperature–humidity combinations. Thus, we could derive the temperature–humidity-dependent thresholds marking the transition from cooling (Δv < 0) to heating (Δv > 0) wind effects on physiological heat strain based on the experimental data. This is a distinguishing feature and particular strength of our study because corresponding experiments with sedentary participants have been limited to a few selected climatic conditions [15,16], so mapping the wind effect thresholds in relation to air temperature and humidity for low-active persons had to rely on simulation studies [17,18,19,20]. Recently, a study with exercising participants [28] aimed at the interaction of wind with combinations of temperature and humidity over a likewise extensive grid but focused on the physical work capacity with cardiac strain clamped at 130 bpm HR.
Our experimental findings suggest beneficial, i.e., strain-reducing wind effects for temperatures below 35 °C, conforming to common public health recommendations [8,9,10]. However, beneficial wind effects were also observed at higher temperatures with humidity levels above 2 kPa water vapor pressure concerning heart rate and core temperature, and 3 kPa concerning skin temperature and sweat rate, respectively (Figure 5). This agreed qualitatively with the corresponding charts produced in previous studies [19,28], though there were discrepancies concerning the humidity level separating beneficial wind cooling from detrimental heating effects, which was specified at 50% relative humidity in the former studies, while our study suggested thresholds in terms of absolute humidity (vapor pressure). Some of these discrepancies may be explained by the 3 h exposures in our study, whereas shorter exposures (1 h) had been used in previous studies [28], with longer exposure times usually resulting in higher strain levels [25,57]. The lower humidity limits for skin temperature and sweat rates compared to core temperature and heart rates in our study may reflect the direct influence of wind speed on the heat transfer from the skin surface by convection and evaporation and on sweating efficiency [11,12,13,58].
For the UTCI, no limiting humidity in terms of vapor pressure was identified, but wind-cooling effects were evident for any temperature condition with a relative humidity above 50% (Figure 3b). The UTCI’s assessment of wind effects showed high positive correlations with the observed changes in physiological responses (Figure 6a), exhibiting the closest agreement (r = 0.9) for skin temperature and sweat rate, where wind is known for elevating the relevant convective and evaporative heat transfer as discussed before. The UTCI also showed good agreement concerning the transition from cooling (Δv < 0) to heating (Δv > 0) wind effects on skin temperatures and sweat rates. The UTCI’s assessment was more conservative, i.e., overpredictive, concerning the potential hazards due to additional heating by wind with respect to core temperature and heart rates (Figure 6b).
The expanded UTCI calculations suggested that further increasing wind speed would amplify the cooling and heating effects depending on temperature and humidity (Figure 7 and Figure A3). In addition, these calculations pointed to a relevant interaction with the radiant heat load, e.g., arising from outdoor exposure to solar radiation. The elevated radiant heat load will increase the UTCI by approximately 3 K per 10 K increase in ΔTmrt [36]. On the other hand, the limit of beneficial wind cooling was shifted toward higher temperatures, counteracting, at least partly, the heat gained from the increased radiant heat load. Similar interaction effects with wind decreasing the heat gain from radiant heat have been reported from thermal manikin measurements with work clothes [59,60]. However, as our results concerning the effects of heat radiation are purely modeled outcomes using UTCI, further supporting experimental studies with human participants concerning the interaction of wind with radiant heat load are required.

4.2. Limitations and Outlook

The selected study population of young, healthy, and fit males certainly limits the generalizability of our findings, though recent research suggests only minimal influences of individual factors on the heat stress assessment under moderate activity levels [61], as applied in our study. As another limitation, the responses of semi-nude participants (Icl = 0.1 clo) might differ from the effects with clothed humans, e.g., as assumed by the UTCI with clothing insulation not falling below Icl = 0.3 clo [34]. On the other hand, prior studies with sedentary participants [15,16] applied similar minimal clothing settings. In addition, limited differences were found concerning the wind effects on the physical work capacity with minimal clothing compared to light and covering clothing [28], and the corresponding temperature–humidity charts in that study were produced for a person wearing Icl = 0.28 clo.
Given these limitations, future experimental and modeling studies should consider the moderating influence of individual factors such as age, sex, or body constitution and of solar radiation and/or clothing with relevance, e.g., for desert conditions (simoom) [21], but also for (semi-)outdoor workplaces.

5. Conclusions

By widening the previous findings from resting persons to moderately exercising individuals, our results confirm the wind-cooling effects for air temperatures below 35 °C, and corroborate and quantify the potential wind-cooling effects at air temperatures above typical skin temperature values in humid climates.
In addition, our findings suggest that the UTCI has the capability to adequately assess the heat stress mitigation by wind depending on temperature and humidity under such conditions. Thus, the UTCI constitutes an easy-to-use tool that may be helpful for establishing evidence-based policies aiming at sustainable heat mitigation strategies under climate change scenarios.

Author Contributions

Conceptualization, P.B. and B.K.; methodology, P.B.; software, P.B.; validation, P.B. and B.K.; formal analysis, P.B.; investigation, B.K.; resources, P.B.; data curation, P.B. and B.K.; writing—original draft preparation, P.B.; writing—review and editing, P.B. and B.K.; visualization, P.B.; supervision, B.K.; project administration, P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable, since this paper makes use of data from prior studies [24,44], which were conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Ethics Committee of IfADo.

Informed Consent Statement

All participants of the prior studies [24,42] provided informed consent.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to data privacy issues.

Acknowledgments

We like to thank the members of the former Environmental Physiology department at IfADo for skillfully running the experiments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

This appendix contains the following supplemental figures:
  • An additional example of recordings of physiological heat strain variables depending on temperature, humidity, and air velocity;
  • Goodness-of-fit plots for the GAMs fitted to the physiological heat strain variables;
  • Influence of radiant heat and wind speed on wind effects assessed by UTCI (ΔvUTCI);
  • Correlations between the wind effects of the physiological heat strain variables and UTCI calculated for va,10m = 4 m/s.
Biology 12 00802 g0a1 550
Figure A1. Time course of heart rate (HR), rectal (Tre) and mean skin temperature (Tsk), and sweat rate (SR) of an acclimated participant during experiments in a hot−humid climate (a) or in a hot−dry climate (b) under both reference (va = 0.3 m/s, red lines) and high air velocity (va = 2 m/s, blue lines) conditions. The yellow shaded area marks the time interval with averaged values used for analyses. Note that the reference wind trial in (a) was prematurely aborted in the 8th working period.
Figure A1. Time course of heart rate (HR), rectal (Tre) and mean skin temperature (Tsk), and sweat rate (SR) of an acclimated participant during experiments in a hot−humid climate (a) or in a hot−dry climate (b) under both reference (va = 0.3 m/s, red lines) and high air velocity (va = 2 m/s, blue lines) conditions. The yellow shaded area marks the time interval with averaged values used for analyses. Note that the reference wind trial in (a) was prematurely aborted in the 8th working period.
Biology 12 00802 g0a1
Biology 12 00802 g0a2 550
Figure A2. Diagnostic goodness−of−fit plots for the generalized additive models (GAMs) fitted to (A) heart rate (bpm), (B) rectal temperature (°C), (C) mean skin temperature (°C), and (D) sweat rate (g/h), respectively. The subplots show Q−Q plots of standardized residuals vs. the theoretical Gaussian quantiles (upper left panels), histograms of residual distribution (upper right), residuals vs. predicted values (lower left), and correlation of observed to fitted values (lower right).
Figure A2. Diagnostic goodness−of−fit plots for the generalized additive models (GAMs) fitted to (A) heart rate (bpm), (B) rectal temperature (°C), (C) mean skin temperature (°C), and (D) sweat rate (g/h), respectively. The subplots show Q−Q plots of standardized residuals vs. the theoretical Gaussian quantiles (upper left panels), histograms of residual distribution (upper right), residuals vs. predicted values (lower left), and correlation of observed to fitted values (lower right).
Biology 12 00802 g0a2
Biology 12 00802 g0a3 550
Figure A3. Effect of radiant heat and wind speed on wind effects assessed by UTCI. Change in UTCI at increased wind speeds (ΔvUTCI) compared to the reference wind condition (va,10m = 0.5 m/s) depending on air temperature and vapor pressure considering the influence of radiant heat load (ΔTmrt = Tmrt − Ta). Note that the upper left plot replicates the data of Figure 3b.
Figure A3. Effect of radiant heat and wind speed on wind effects assessed by UTCI. Change in UTCI at increased wind speeds (ΔvUTCI) compared to the reference wind condition (va,10m = 0.5 m/s) depending on air temperature and vapor pressure considering the influence of radiant heat load (ΔTmrt = Tmrt − Ta). Note that the upper left plot replicates the data of Figure 3b.
Biology 12 00802 g0a3
Biology 12 00802 g0a4 550
Figure A4. Correlation of wind effects from Figure 5 for heart rate (ΔvHR), rectal temperature (ΔvTre), mean skin temperature (ΔvTsk), and sweat rate (ΔvSR), respectively, with the wind effect on UTCI (ΔvUTCI) calculated for an elevated wind speed (4 m/s).
Figure A4. Correlation of wind effects from Figure 5 for heart rate (ΔvHR), rectal temperature (ΔvTre), mean skin temperature (ΔvTsk), and sweat rate (ΔvSR), respectively, with the wind effect on UTCI (ΔvUTCI) calculated for an elevated wind speed (4 m/s).
Biology 12 00802 g0a4

References

  1. Ebi, K.L.; Capon, A.; Berry, P.; Broderick, C.; de Dear, R.; Havenith, G.; Honda, Y.; Kovats, R.S.; Ma, W.; Malik, A.; et al. Hot weather and heat extremes: Health risks. Lancet 2021, 398, 698–708. [Google Scholar] [CrossRef] [PubMed]
  2. Romanello, M.; Di Napoli, C.; Drummond, P.; Green, C.; Kennard, H.; Lampard, P.; Scamman, D.; Arnell, N.; Ayeb-Karlsson, S.; Ford, L.B.; et al. The 2022 report of the Lancet Countdown on health and climate change: Health at the mercy of fossil fuels. Lancet 2022, 400, 1619–1654. [Google Scholar] [CrossRef] [PubMed]
  3. Foster, J.; Smallcombe, J.W.; Hodder, S.; Jay, O.; Flouris, A.D.; Nybo, L.; Havenith, G. An advanced empirical model for quantifying the impact of heat and climate change on human physical work capacity. Int. J. Biometeorol. 2021, 65, 1215–1229. [Google Scholar] [CrossRef] [PubMed]
  4. Bröde, P.; Fiala, D.; Lemke, B.; Kjellstrom, T. Estimated work ability in warm outdoor environments depends on the chosen heat stress assessment metric. Int. J. Biometeorol. 2018, 62, 331–345. [Google Scholar] [CrossRef] [PubMed]
  5. Jay, O.; Capon, A.; Berry, P.; Broderick, C.; de Dear, R.; Havenith, G.; Honda, Y.; Kovats, R.S.; Ma, W.; Malik, A.; et al. Reducing the health effects of hot weather and heat extremes: From personal cooling strategies to green cities. Lancet 2021, 398, 709–724. [Google Scholar] [CrossRef]
  6. Malik, A.; Bongers, C.; McBain, B.; Rey-Lescure, O.; Dear, R.d.; Capon, A.; Lenzen, M.; Jay, O. The potential for indoor fans to change air conditioning use while maintaining human thermal comfort during hot weather: An analysis of energy demand and associated greenhouse gas emissions. Lancet Planet. Health 2022, 6, e301–e309. [Google Scholar] [CrossRef]
  7. Jay, O.; Hoelzl, R.; Weets, J.; Morris, N.; English, T.; Nybo, L.; Niu, J.; de Dear, R.; Capon, A. Fanning as an alternative to air conditioning—A sustainable solution for reducing indoor occupational heat stress. Energy Build. 2019, 193, 92–98. [Google Scholar] [CrossRef]
  8. CDC. Frequently Asked Questions (FAQ) About Extreme Heat; Centers for Disease Control and Prevention: Atlanta, GA, USA, 2012. Available online: https://www.cdc.gov/disasters/extremeheat/faq.html (accessed on 15 January 2020).
  9. EPA. Excessive Heat Events Guidebook; EPA 430-B-16-001; United States Environmental Protection Agency: Washington, DC, USA, 2016. Available online: https://www.epa.gov/sites/production/files/2016-03/documents/eheguide_final.pdf (accessed on 15 January 2020).
  10. McGregor, G.R.; Bessemoulin, P.; Ebi, K.; Menne, B. Heatwaves and Health: Guidance on Warning-System Development; WMO-No. 1142; World Meteorological Organization and World Health Organization: Geneva, Switzerland, 2015; Available online: https://library.wmo.int/doc_num.php?explnum_id=3371 (accessed on 15 January 2020).
  11. Candas, V.; Libert, J.P.; Vogt, J.J. Influence of air velocity and heat acclimation on human skin wettedness and sweating efficiency. J. Appl. Physiol. 1979, 47, 1194–1200. [Google Scholar] [CrossRef]
  12. Adams, W.C.; Mack, G.W.; Langhans, G.W.; Nadel, E.R. Effects of varied air velocity on sweating and evaporative rates during exercise. J. Appl. Physiol. 1992, 73, 2668–2674. [Google Scholar] [CrossRef]
  13. Nadel, E.R.; Stolwijk, J.A. Effect of skin wettedness on sweat gland response. J. Appl. Physiol. 1973, 35, 689–694. [Google Scholar] [CrossRef]
  14. Jay, O.; Capon, A. Use of physiological evidence for heatwave public policy. Lancet Planet. Health 2018, 2, e10. [Google Scholar] [CrossRef]
  15. Morris, N.B.; English, T.; Hospers, L.; Capon, A.; Jay, O. The Effects of Electric Fan Use Under Differing Resting Heat Index Conditions: A Clinical Trial. Ann. Intern. Med. 2019, 171, 675–677. [Google Scholar] [CrossRef]
  16. Ravanelli, N.M.; Hodder, S.G.; Havenith, G.; Jay, O. Heart Rate and Body Temperature Responses to Extreme Heat and Humidity with and Without Electric Fans. JAMA 2015, 313, 724–725. [Google Scholar] [CrossRef]
  17. Morris, N.B.; Chaseling, G.K.; English, T.; Gruss, F.; Maideen, M.F.B.; Capon, A.; Jay, O. Electric fan use for cooling during hot weather: A biophysical modelling study. Lancet Planet. Health 2021, 5, e368–e377. [Google Scholar] [CrossRef]
  18. Tartarini, F.; Schiavon, S.; Jay, O.; Arens, E.; Huizenga, C. Application of Gagge’s energy balance model to determine humidity-dependent temperature thresholds for healthy adults using electric fans during heatwaves. Build. Environ. 2021, 207, 108437. [Google Scholar] [CrossRef]
  19. Hospers, L.; Smallcombe, J.W.; Morris, N.B.; Capon, A.; Jay, O. Electric fans: A potential stay-at-home cooling strategy during the COVID-19 pandemic this summer? Sci. Total Environ. 2020, 747, 141180. [Google Scholar] [CrossRef]
  20. Jay, O.; Cramer, M.N.; Ravanelli, N.M.; Hodder, S.G. Should electric fans be used during a heat wave? Appl. Ergon. 2015, 46, 137–143. [Google Scholar] [CrossRef]
  21. Normand, C.W.B. The effect of high temperature, humidity, and wind on the human body. Q. J. R. Meteorol. Soc. 1920, 46, 1–14. [Google Scholar] [CrossRef]
  22. Ainsworth, B.E.; Haskell, W.L.; Herrmann, S.D.; Meckes, N.; Bassett, D.R.J.; Tudor-Locke, C.; Greer, J.L.; Vezina, J.; Whitt-Glover, M.C.; Leon, A.S. 2011 Compendium of Physical Activities: A Second Update of Codes and MET Values. Med. Sci. Sport. Exerc. 2011, 43, 1575–1581. [Google Scholar] [CrossRef]
  23. ISO 8996; Ergonomics of the Thermal Environment—Determination of Metabolic Rate. International Organisation for Standardisation: Geneva, Switzerland, 2021.
  24. Kampmann, B. Zur Physiologie der Arbeit in warmem Klima. Ergebnisse aus Laboruntersuchungen und aus Feldstudien im Steinkohlenbergbau. Habilitation Thesis, Bergische Universität Wuppertal, Wuppertal, Germany, 2000. [Google Scholar]
  25. Kampmann, B.; Bröde, P. Do one-hour exposures provide a valid assessment of physiological heat strain? Z. Für Arb. 2022, 76, 105–117. [Google Scholar] [CrossRef]
  26. Brown, W.K.; Sargent, F. Hidromeiosis. Arch. Environ. Health Int. J. 1965, 11, 442–453. [Google Scholar] [CrossRef]
  27. Candas, V.; Libert, J.; Vogt, J. Effect of hidromeiosis on sweat drippage during acclimation to humid heat. Eur. J. Appl. Physiol. Occup. Physiol. 1980, 44, 123–133. [Google Scholar] [CrossRef]
  28. Foster, J.; Smallcombe, J.W.; Hodder, S.; Jay, O.; Flouris, A.D.; Havenith, G. Quantifying the impact of heat on human physical work capacity; part II: The observed interaction of air velocity with temperature, humidity, sweat rate, and clothing is not captured by most heat stress indices. Int. J. Biometeorol. 2021, 66, 507–520. [Google Scholar] [CrossRef]
  29. Kamon, E.; Avellini, B. Wind speed limits to work under hot environments for clothed men. J. Appl. Physiol. 1979, 46, 340–345. [Google Scholar] [CrossRef]
  30. Graham, C.; Lynch, G.P.; English, T.; Hospers, L.; Jay, O. Optimal break structures and cooling strategies to mitigate heat stress during a Rugby League match simulation. J. Sci. Med. Sport 2021, 24, 793–799. [Google Scholar] [CrossRef]
  31. Wang, F.; Deng, Q.; Lei, T.-H.; Wang, X.; Wang, A.R. Biophysical modelling predicts unreliable core temperature responses on healthy older adults using electric fans at residential homes during heatwaves. Build. Environ. 2023, 228, 109888. [Google Scholar] [CrossRef]
  32. Jendritzky, G.; de Dear, R.; Havenith, G. UTCI—Why another thermal index? Int. J. Biometeorol. 2012, 56, 421–428. [Google Scholar] [CrossRef]
  33. 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]
  34. Havenith, G.; Fiala, D.; Blazejczyk, K.; Richards, M.; Bröde, P.; Holmér, I.; Rintamaki, H.; Ben Shabat, Y.; Jendritzky, G. The UTCI-clothing model. Int. J. Biometeorol. 2012, 56, 461–470. [Google Scholar] [CrossRef]
  35. Psikuta, A.; Fiala, D.; Laschewski, G.; Jendritzky, G.; Richards, M.; Blazejczyk, K.; Mekjavic, I.B.; Rintamäki, H.; de Dear, R.; Havenith, G. Validation of the Fiala multi-node thermophysiological model for UTCI application. Int. J. Biometeorol. 2012, 56, 443–460. [Google Scholar] [CrossRef]
  36. Bröde, P.; Blazejczyk, K.; Fiala, D.; Havenith, G.; Holmér, I.; Jendritzky, G.; Kuklane, K.; Kampmann, B. The Universal Thermal Climate Index UTCI Compared to Ergonomics Standards for Assessing the Thermal Environment. Ind. Health 2013, 51, 16–24. [Google Scholar] [CrossRef]
  37. Kampmann, B.; Bröde, P.; Fiala, D. Physiological responses to temperature and humidity compared to the assessment by UTCI, WGBT and PHS. Int. J. Biometeorol. 2012, 56, 505–513. [Google Scholar] [CrossRef]
  38. Ioannou, L.G.; Tsoutsoubi, L.; Mantzios, K.; Vliora, M.; Nintou, E.; Piil, J.F.; Notley, S.R.; Dinas, P.C.; Gourzoulidis, G.A.; Havenith, G.; et al. Indicators to assess physiological heat strain—Part 3: Multi-country field evaluation and consensus recommendations. Temperature 2022, 9, 274–291. [Google Scholar] [CrossRef]
  39. Bröde, P.; Fiala, D.; Blazejczyk, K.; Holmér, I.; Jendritzky, G.; Kampmann, B.; Tinz, B.; Havenith, G. Deriving the operational procedure for the Universal Thermal Climate Index (UTCI). Int. J. Biometeorol. 2012, 56, 481–494. [Google Scholar] [CrossRef]
  40. Sadeghi, M.; de Dear, R.; Wood, G.; Samali, B. Development of a bioclimatic wind rose tool for assessment of comfort wind resources in Sydney, Australia for 2013 and 2030. Int. J. Biometeorol. 2018, 62, 1963–1972. [Google Scholar] [CrossRef]
  41. Foster, J.; Smallcombe, J.W.; Hodder, S.; Jay, O.; Flouris, A.D.; Nybo, L.; Havenith, G. Quantifying the impact of heat on human physical work capacity; part III: The impact of solar radiation varies with air temperature, humidity, and clothing coverage. Int. J. Biometeorol. 2022, 66, 175–188. [Google Scholar] [CrossRef]
  42. Wenzel, H.G.; Mehnert, C.; Schwarzenau, P. Evaluation of tolerance limits for humans under heat stress and the problems involved. Scand. J. Work. Environ. Health 1989, 15, 7–14. [Google Scholar]
  43. Kampmann, B.; Bröde, P. Heat Acclimation Does Not Modify Q10 and Thermal Cardiac Reactivity. Front. Physiol. 2019, 10, 1524. [Google Scholar] [CrossRef]
  44. Wood, S.N. Generalized Additive Models: An Introduction with R, 2nd ed.; Chapman & Hall/CRC: Boca Raton, FL, USA, 2017. [Google Scholar]
  45. Zuur, A.F.; Ieno, E.N.; Walker, N.J.; Saveliev, A.A.; Smith, G.M. Things are not Always Linear; Additive Modelling. In Mixed Effects Models and Extensions in Ecology with R; Zuur, A.F., Ieno, E.N., Walker, N.J., Saveliev, A.A., Smith, G.M., Eds.; Springer: New York, NY, USA, 2009; pp. 35–69. [Google Scholar] [CrossRef]
  46. Wood, S.N. A simple test for random effects in regression models. Biometrika 2013, 100, 1005–1010. [Google Scholar] [CrossRef]
  47. Wood, S.N. On p-values for smooth components of an extended generalized additive model. Biometrika 2013, 100, 221–228. [Google Scholar] [CrossRef]
  48. Fasiolo, M.; Nedellec, R.; Goude, Y.; Wood, S.N. Scalable visualisation methods for modern Generalized Additive Models. J. Comput. Graph. Stat. 2020, 29, 78–86. [Google Scholar] [CrossRef]
  49. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022. [Google Scholar]
  50. Coretta, S. Tidymv: Tidy Model Visualisation for Generalised Additive Models. R Package Version 3.3.1. 2022. Available online: https://CRAN.R-project.org/package=tidymv (accessed on 9 February 2023).
  51. ISO 11079; Ergonomics of the Thermal Environment—Determination and Interpretation of Cold Stress When Using Required Clothing Insulation (IREQ) and Local Cooling Effects. International Organisation for Standardisation: Geneva, Switzerland, 2007.
  52. ISO 9920; Ergonomics of the Thermal Environment—Estimation of Thermal insulation and Water Vapour Resistance of a Clothing Ensemble. International Organisation for Standardisation: Geneva, Switzerland, 2007.
  53. Bröde, P. Issues in UTCI Calculation from a Decade’s Experience. In Applications of the Universal Thermal Climate Index UTCI in Biometeorology: Latest Developments and Case Studies; Krüger, E.L., Ed.; Springer International Publishing: Cham, Switzerland, 2021; pp. 13–21. [Google Scholar] [CrossRef]
  54. Cheuvront, S.; Montain, S.; Goodman, D.; Blanchard, L.; Sawka, M. Evaluation of the limits to accurate sweat loss prediction during prolonged exercise. Eur. J. Appl. Physiol. 2007, 101, 215–224. [Google Scholar] [CrossRef]
  55. Wood, S.N. Ref.df in Mgcv Gam. Available online: https://stat.ethz.ch/pipermail/r-help/2019-March/462135.html (accessed on 9 February 2023).
  56. Richmond, V.L.; Davey, S.; Griggs, K.; Havenith, G. Prediction of Core Body Temperature from Multiple Variables. Ann. Occup. Hyg. 2015, 59, 1168–1178. [Google Scholar] [CrossRef]
  57. Smallcombe, J.W.; Foster, J.; Hodder, S.G.; Jay, O.; Flouris, A.D.; Havenith, G. Quantifying the impact of heat on human physical work capacity; part IV: Interactions between work duration and heat stress severity. Int. J. Biometeorol. 2022, 66, 2463–2476. [Google Scholar] [CrossRef]
  58. Nadel, E.R.; Bullard, R.W.; Stolwijk, J.A. Importance of skin temperature in the regulation of sweating. J. Appl. Physiol. 1971, 31, 80–87. [Google Scholar] [CrossRef]
  59. Bröde, P.; Kuklane, K.; Candas, V.; den Hartog, E.A.; Griefahn, B.; Holmér, I.; Meinander, H.; Nocker, W.; Richards, M.; Havenith, G. Heat Gain from Thermal Radiation Through Protective Clothing with Different Insulation, Reflectivity and Vapour Permeability. Int. J. Occup. Saf. Ergon. 2010, 16, 231–244. [Google Scholar] [CrossRef]
  60. Den Hartog, E.A.; Havenith, G. Analytical study of the heat loss attenuation by clothing on thermal manikins under radiative heat loads. Int. J. Occup. Saf. Ergon. 2010, 16, 245–261. [Google Scholar] [CrossRef]
  61. Wolf, S.T.; Havenith, G.; Kenney, W.L. Relatively minor influence of individual characteristics on critical wet-bulb globe temperature (WBGT) limits during light activity in young adults (PSU HEAT Project). J. Appl. Physiol. 2023, 134, 1216–1223. [Google Scholar] [CrossRef]
Biology 12 00802 g001 550
Figure 1. Time course of heart rate (HR), rectal (Tre) and mean skin temperature (Tsk), and sweat rate (SR) of an acclimated participant during experiments in a hot–humid climate (a) or in a hot–dry climate (b) under both reference (va = 0.3 m/s, red lines) and high air velocity (va = 2 m/s, blue lines) conditions. The yellow shaded area marks the time interval with averaged values used for analyses. Note that the reference wind trial in (a) was prematurely aborted in the 8th working period.
Figure 1. Time course of heart rate (HR), rectal (Tre) and mean skin temperature (Tsk), and sweat rate (SR) of an acclimated participant during experiments in a hot–humid climate (a) or in a hot–dry climate (b) under both reference (va = 0.3 m/s, red lines) and high air velocity (va = 2 m/s, blue lines) conditions. The yellow shaded area marks the time interval with averaged values used for analyses. Note that the reference wind trial in (a) was prematurely aborted in the 8th working period.
Biology 12 00802 g001
Biology 12 00802 g002 550
Figure 2. Relative wind speed at person level (var,1m) when moving with the UTCI reference walking speed (vw) of 4 km/h, corresponding to 1.1 m/s, in relation to wind speed measured at person level (va,1m) and 10 m above ground (va,10m). The insert provides a detailed view of the gray-shaded region for va,10m ≤ 4 m/s. Figures include solid lines of identity, and horizontal dashed lines indicating walking speed (vw).
Figure 2. Relative wind speed at person level (var,1m) when moving with the UTCI reference walking speed (vw) of 4 km/h, corresponding to 1.1 m/s, in relation to wind speed measured at person level (va,1m) and 10 m above ground (va,10m). The insert provides a detailed view of the gray-shaded region for va,10m ≤ 4 m/s. Figures include solid lines of identity, and horizontal dashed lines indicating walking speed (vw).
Biology 12 00802 g002
Biology 12 00802 g003 550
Figure 3. Psychrometric charts showing (a) UTCI contours related to air temperature and humidity (vapor pressure) representing the upper limits of ‘no thermal stress’ (UTCI = 26 °C), and ‘moderate’ (32 °C), ‘strong’ (38 °C), and ‘very strong’ (46 °C) heat stress, as well as an illustrative value for ‘extreme heat stress’ (55 °C), respectively. The influence of increased wind speed (HiWind with va,10m = 3 m/s, corresponding to va,1m = 2 m/s: dashed lines) is compared to the reference RefWind (va,10m = 0.5 m/s, corresponding to va,1m = 0.3 m/s: solid lines). The dot−dashed line connecting the intersections of RefWind with HiWind contours represents zero wind effects (Δv = 0). This line is redrawn in (b) as white contour together with the colored regions, indicating temperature–humidity combinations with cooling (blueish) or heating (reddish) wind effects on UTCI (ΔvUTCI), according to Equation (5), respectively, where color intensity indicates magnitude. Gray contours mark relative humidity (rH) levels. All values were calculated without additional radiant heat load (ΔTmrt = Tmrt − Ta = 0 K).
Figure 3. Psychrometric charts showing (a) UTCI contours related to air temperature and humidity (vapor pressure) representing the upper limits of ‘no thermal stress’ (UTCI = 26 °C), and ‘moderate’ (32 °C), ‘strong’ (38 °C), and ‘very strong’ (46 °C) heat stress, as well as an illustrative value for ‘extreme heat stress’ (55 °C), respectively. The influence of increased wind speed (HiWind with va,10m = 3 m/s, corresponding to va,1m = 2 m/s: dashed lines) is compared to the reference RefWind (va,10m = 0.5 m/s, corresponding to va,1m = 0.3 m/s: solid lines). The dot−dashed line connecting the intersections of RefWind with HiWind contours represents zero wind effects (Δv = 0). This line is redrawn in (b) as white contour together with the colored regions, indicating temperature–humidity combinations with cooling (blueish) or heating (reddish) wind effects on UTCI (ΔvUTCI), according to Equation (5), respectively, where color intensity indicates magnitude. Gray contours mark relative humidity (rH) levels. All values were calculated without additional radiant heat load (ΔTmrt = Tmrt − Ta = 0 K).
Biology 12 00802 g003
Biology 12 00802 g004 550
Figure 4. Psychrometric charts including contour lines with colors indicating separate levels of predicted values fitted by GAMs to (a) heart rates (HR), (b) rectal temperatures (Tre), (c) mean skin temperatures (Tsk), and (d) sweat rates (SR) in relation to air temperature and vapor pressure under RefWind (solid lines) and HiWind conditions (dashed lines), respectively. Dot-dashed lines connecting the intersections of RefWind with HiWind contours indicate zero wind effects (Δv = 0).
Figure 4. Psychrometric charts including contour lines with colors indicating separate levels of predicted values fitted by GAMs to (a) heart rates (HR), (b) rectal temperatures (Tre), (c) mean skin temperatures (Tsk), and (d) sweat rates (SR) in relation to air temperature and vapor pressure under RefWind (solid lines) and HiWind conditions (dashed lines), respectively. Dot-dashed lines connecting the intersections of RefWind with HiWind contours indicate zero wind effects (Δv = 0).
Biology 12 00802 g004
Biology 12 00802 g005 550
Figure 5. Difference plots with contours of the wind effects (ΔV) estimated by GAMs in relation to air temperature and vapor pressure for (a) heart rates (HR), (b) rectal temperatures (Tre), (c) mean skin temperatures (Tsk), and (d) sweat rates (SR). White areas indicate temperature−humidity regions with non−significant wind effects (p > 0.05 testing for Δv = 0 vs. Δv ≠ 0), the blueish areas mark significant reductions in physiological strain by wind cooling, while reddish areas depict significantly increased strain levels due to heating wind effects, respectively. Dashed lines show relative humidity levels, and the gray shaded areas indicate combinations not supported by data. Experimental conditions are marked by ‘+’ (RefWind) and ‘×’ (HiWind), respectively.
Figure 5. Difference plots with contours of the wind effects (ΔV) estimated by GAMs in relation to air temperature and vapor pressure for (a) heart rates (HR), (b) rectal temperatures (Tre), (c) mean skin temperatures (Tsk), and (d) sweat rates (SR). White areas indicate temperature−humidity regions with non−significant wind effects (p > 0.05 testing for Δv = 0 vs. Δv ≠ 0), the blueish areas mark significant reductions in physiological strain by wind cooling, while reddish areas depict significantly increased strain levels due to heating wind effects, respectively. Dashed lines show relative humidity levels, and the gray shaded areas indicate combinations not supported by data. Experimental conditions are marked by ‘+’ (RefWind) and ‘×’ (HiWind), respectively.
Biology 12 00802 g005
Biology 12 00802 g006 550
Figure 6. (a) Correlation of wind effects from Figure 5 for heart rate (ΔvHR), rectal temperature (ΔvTre), mean skin temperature (ΔvTsk), and sweat rate (ΔvSR), respectively, with the wind effect on UTCI (ΔvUTCI), shown in Figure 3b, calculated for the temperature–humidity combinations of the experiments (cf. Figure 5). Spearman correlations (rs) and p−values are shown with smoothing splines and 95% confidence bands. Dashed vertical and horizontal lines indicate zero wind effects (Δv = 0) for UTCI and the physiological responses, respectively, where negative values (Δv < 0) indicate wind cooling and positive values (Δv > 0) heating effects. (b) UTCI performance in predicting heating wind effects (Δv > 0) on HR, Tre, Tsk, and SR, respectively, expressed by sensitivity, specificity and accuracy.
Figure 6. (a) Correlation of wind effects from Figure 5 for heart rate (ΔvHR), rectal temperature (ΔvTre), mean skin temperature (ΔvTsk), and sweat rate (ΔvSR), respectively, with the wind effect on UTCI (ΔvUTCI), shown in Figure 3b, calculated for the temperature–humidity combinations of the experiments (cf. Figure 5). Spearman correlations (rs) and p−values are shown with smoothing splines and 95% confidence bands. Dashed vertical and horizontal lines indicate zero wind effects (Δv = 0) for UTCI and the physiological responses, respectively, where negative values (Δv < 0) indicate wind cooling and positive values (Δv > 0) heating effects. (b) UTCI performance in predicting heating wind effects (Δv > 0) on HR, Tre, Tsk, and SR, respectively, expressed by sensitivity, specificity and accuracy.
Biology 12 00802 g006
Biology 12 00802 g007 550
Figure 7. Boxplots including the percentiles P0 (minimum), P25 (1st quartile), P50 (median), P75 (3rd quartile), and P100 (maximum) summarizing the effect of three levels of increased wind speed (va,10m) on UTCI (ΔvUTCI) over air temperatures between 25 and 50 °C and vapor pressure of 0.1−5 kPa in relation to radiant heat load (ΔTmrt = Tmrt − Ta), cf. Figure A3.
Figure 7. Boxplots including the percentiles P0 (minimum), P25 (1st quartile), P50 (median), P75 (3rd quartile), and P100 (maximum) summarizing the effect of three levels of increased wind speed (va,10m) on UTCI (ΔvUTCI) over air temperatures between 25 and 50 °C and vapor pressure of 0.1−5 kPa in relation to radiant heat load (ΔTmrt = Tmrt − Ta), cf. Figure A3.
Biology 12 00802 g007
Table
Table 1. Fitted generalized additive models (GAMs) predicting experimental data on heart rate (HR), rectal temperature (Tre), skin temperature (Tsk), and sweat rate (SR), respectively, depending on combinations of air temperature (Ta) and vapor pressure (pa), considering wind condition (RefWind, HiWind) as modifier and including a random subject effect. Abbreviations for model parameters, as in Equation (2).
Table 1. Fitted generalized additive models (GAMs) predicting experimental data on heart rate (HR), rectal temperature (Tre), skin temperature (Tsk), and sweat rate (SR), respectively, depending on combinations of air temperature (Ta) and vapor pressure (pa), considering wind condition (RefWind, HiWind) as modifier and including a random subject effect. Abbreviations for model parameters, as in Equation (2).
HR (bpm)Tre (°C)Tsk (°C)SR (g/h)
Observations (#missing values)189 (9)198 (0)189 (9)186 (13)
Goodness-of-fit
Adjusted R2 (%)78.876.288.689.9
Residual standard error7.60.20.4120.1
Intercept µ
Mean estimate102.537.635.5744.4
SE0.80.10.227.1
p-value<0.0001<0.0001<0.0001<0.0001
s(ID)
edf0.53.83.83.4
Ref.df4.04.04.04.0
F-value0.227.124.36.1
p-value0.2018<0.0001<0.0001<0.0001
te(Ta, pa)
edf14.66.510.310.0
Ref.df19.08.214.013.4
F-value18.934.731.551.3
p-value<0.0001<0.0001<0.0001<0.0001
teΔv(Ta, pa)
edf4.04.05.54.5
Ref.df4.14.06.24.8
F-value10.87.04.88.3
p-value<0.0001<0.00010.0001<0.0001
Notes: SE: standard error; edf: effective degrees-of-freedom; Ref.df: reference degrees-of-freedom [55]; s(ID): random subject effect spline; te(Ta, pa): temperature–humidity effect (bivariate tensor product spline); teΔv(Ta, pa): temperature–humidity-dependent wind effect of HiWind compared to RefWind (bivariate interaction spline).
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

Bröde, P.; Kampmann, B. Temperature–Humidity-Dependent Wind Effects on Physiological Heat Strain of Moderately Exercising Individuals Reproduced by the Universal Thermal Climate Index (UTCI). Biology 2023, 12, 802. https://doi.org/10.3390/biology12060802

AMA Style

Bröde P, Kampmann B. Temperature–Humidity-Dependent Wind Effects on Physiological Heat Strain of Moderately Exercising Individuals Reproduced by the Universal Thermal Climate Index (UTCI). Biology. 2023; 12(6):802. https://doi.org/10.3390/biology12060802

Chicago/Turabian Style

Bröde, Peter, and Bernhard Kampmann. 2023. "Temperature–Humidity-Dependent Wind Effects on Physiological Heat Strain of Moderately Exercising Individuals Reproduced by the Universal Thermal Climate Index (UTCI)" Biology 12, no. 6: 802. https://doi.org/10.3390/biology12060802

APA Style

Bröde, P., & Kampmann, B. (2023). Temperature–Humidity-Dependent Wind Effects on Physiological Heat Strain of Moderately Exercising Individuals Reproduced by the Universal Thermal Climate Index (UTCI). Biology, 12(6), 802. https://doi.org/10.3390/biology12060802

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