Abstract
Introduction:
Hot-dry (HD) and warm-humid (WH) environments with equivalent wet-bulb globe temperature (WBGT) may be expected to produce similar physiological strain; however, previous studies have reported divergent thermoregulatory and cardiovascular responses. This study compared heat storage, thermal, and cardiovascular responses between HD and WH conditions with matched WBGT during cycling.
Methods:
Ten healthy, active, non-heat-acclimatized males completed two randomized crossover trials consisting of 60 minutes of cycling at an individualized workload corresponding to 55% of peak power output in HD (39.2°C, 30% relative humidity [RH]) and WH (34.2°C, 55% RH) conditions with equivalent WBGT (~30°C). Heat storage, dry heat loss (DHL), and evaporative heat loss (EHL) were estimated using partitional calorimetry. Rectal temperature (Tre), mean skin temperature (Tsk), and cardiovascular responses were measured throughout the trials.
Results:
No condition differences were observed for cumulative heat storage, Tre, heart rate, blood pressure, or forearm skin blood flux (all p > .050). Tsk (HD: 36.28 ± 0.67°C, WH: 35.89 ± 0.70°C; p = .030, ηp² = .465) and EHL (HD: 450 ± 53 W, WH: 319 ± 30 W; p < .001, ηp² = .932) were greater in HD, whereas DHL was greater in WH (HD: −80 ± 25 W; WH: 48 ± 25 W; p < .001, ηp² = .950). Conclusion: Despite similar heat storage and cardiovascular strain, matched-WBGT environments altered the mechanisms of heat exchange. Thus, ambient temperature, humidity, and workload should be considered alongside WBGT when monitoring heat strain and selecting heat-mitigation strategies.
Conclusion:
Despite similar heat storage and cardiovascular strain, matched-WBGT environments altered the mechanisms of heat exchange. Thus, ambient temperature, humidity, and workload should be considered alongside WBGT when monitoring heat strain and selecting heat-mitigation strategies.
Introduction
Global heatwaves, which have increased in frequency, duration, and intensity, have led to numerous health consequences over the last several decades (). An estimate suggests that around a third of globe’s population is exposed to deadly heatwaves and this percentage will increase to half by 2100 (). With a warming climate, the potential for more frequent and severe heat events could significantly worsen heat-related morbidity and mortality (), especially for settings such as athletic (), occupational (), and military operations () that occur during extreme heat stress with moderate to vigorous intensity physical activity. Notably, researchers found that heat disproportionately kills young people aged 18 to 34 years old (). Even though young adults are usually more physiologically robust to heat, social, economic, and behavioral factors can lead to higher heat-related mortality in this age group ().
Heatwaves have constantly threatened human health and safety (). For instance, Europe sustained a severe heatwave that resulted in over 70,000 heat-related deaths in 2003 (). Two notable heatwave types, namely, hot-dry (HD) and warm-humid (WH), endanger occupational activity, military operations, and athletic events. It is unclear, however, whether there is a difference in physiological stress and subsequent risk of heat-related illness between HD and WH heatwaves (). It is imperative to differentiate the impact on physiological responses for these differing types of heat stress. Especially when a single metric (e.g., wet-bulb globe temperature [WBGT]) is used to quantify heat stress, regardless of type. In addition, a previous study reported that WBGT underestimates heat strain in high relative humidity (RH) environments (). Others revealed that HD and WH conditions with matched WBGT induced divergent thermal and cardiovascular responses (; ; ; ). It is necessary to identify heat strain between HD and WH conditions with identical heat stress values and subsequently help develop a targeted policy and safety strategy to mitigate future heat-related morbidity and mortality.
Furthermore, many previous studies have investigated the effect of heat stress on thermoregulatory and cardiovascular responses, mostly focused on minimal to low intensity physical activity due to the safety considerations (; ; ; ; ; ; ). found heat storage and gastrointestinal temperature were no different between HD and WH conditions for minimal and low intensity physical activity. also indicated no differences in rectal temperature and heart rate between HD and WH conditions during low intensity walking. In contrast, suggested that physiological stress (heart rate, core temperature, and whole body sweat rate) was greater in HD conditions than in WH conditions during low intensity cycling. In consideration of the demand on adjusted safety requirements at varying exercise intensity, it is urgent to evaluate thermoregulatory and cardiovascular differences during more intense exercise with extreme heat stress.
Despite the general population reducing outdoor physical activity to avoid extreme heat stress, occupational groups (e.g., construction and utility workers), athletic events, and military operations occur in thermally stressful environments with moderate or vigorous activity. However, most studies comparing matched WBGT during exercise in the heat focus on minimal to low intensity (; ; ; ; ). Given that it is inherent that moderate to vigorous activity occurs daily through recreation, industry and the military, it is imperative we investigate this work rate in terms of varying heat stress. This area of research is essential to develop, or confirm, safety for individuals in hot environments.
Therefore, the primary aim of this investigation was to compare thermoregulatory and cardiovascular responses between HD and WH conditions with matched WBGT during 60 minutes of moderate intensity cycling among healthy male adults. We hypothesized that heat storage, rectal temperature, and heart rate would not differ between HD and WH conditions. Still, there would be more dry heat gain and evaporative heat loss (EHL) in HD compared to WH conditions. Further, we assumed dry heat loss (DHL) would be greater in WH than in HD conditions.
Materials and methods
Participants
Fifteen healthy, active [physically active at least two days/week; peak oxygen consumption (VO2peak) ≥ 45 ml/kg/min] and non-heat-acclimatized males volunteered and completed this study. Data from five participants were excluded prior to analysis because of factors that could substantially influence physiological responses to exercise in the heat. Three participants subsequently disclosed illness that occurred during the experimental period after atypical physiological responses prompted additional questioning, and two participants did not comply with pre-trial preparation instructions (one did not refrain from exercise within 24 h of testing and one reported altered sleep prior to a trial). These exclusions were made to preserve participant safety and data integrity and were not based on the study outcomes. Demographics for our final N = 10 participants appear in Table 1.
Table 1
| Characteristic | Mean (SD) | Range |
|---|---|---|
| Age (y) | 29 (7) | 19–39 |
| Height (m) | 1.79 (0.11) | 1.62–1.92 |
| Body Mass (kg) | 77.4 (9.3) | 59.9–87.3 |
| Body Fat (%) | 18.2 (5.6) | 7.5–24.3 |
| BMI (kg/m2) | 23.9 (1.4) | 21.2–25.7 |
| AD (m2) | 1.96 (0.18) | 1.64–2.15 |
| AD/wt (m2/kg) | 0.025 (0.001) | 0.025–0.027 |
| VO2peak (mL/kg/min) | 51.9 (5.3) | 44.9–62.0 |
| VO2peak (L/min) | 4.00 (0.53) | 3.14–4.79 |
Participant characteristics (n = 10).
BMI, body mass index; AD, body surface area; AD/wt, body surface area-to-mass-ratio; VO2peak, peak oxygen consumption.
All experimental procedures were approved by the University Institutional Review Board prior to data collection (IRB approved protocol number: 2309490663). Written informed consent was obtained from all participants before participation. All participants completed a medical history questionnaire and verified no exclusion criteria. Participants were originally excluded if they self-reported previous heat illness within the past three years, had a current musculoskeletal injury, or were taking medications known to affect thermoregulation. We excluded participants who use tobacco or electronic cigarettes. Participants were not allowed to participate if they reported cardiovascular, metabolic, renal, gastrointestinal or liver disease on screening medical history forms. Only male participants were included because fluctuations in ovarian hormones across the menstrual cycle can influence thermoregulatory responses, including sweating, skin blood flux, and core temperature regulation (). Restricting the sample to males reduced physiological variability and allowed for a more controlled examination of environmental effects on thermoregulatory and cardiovascular responses. All participants were asked to refrain from alcohol for 24 hours leading into trials and caffeine ingestion was matched between trials. Strenuous activity was restricted for 24 hours before both preliminary and experimental trials. Participants wore a T-shirt or cycling jersey, shorts, socks, and cycling shoes that were matched for both trials.
Preliminary trial
During a preliminary visit, participants were familiarized with the experimental protocol and a body composition was assessed via dual energy x-ray absorptiometry prior to the experimental trial (Lunar Prodigy, General Electric, Madison, WI, USA). VO2peak was assessed via a cycle ergometer (Racer Mate, Seattle, WA) graded exercise protocol. Participants warmed up at 100W for five minutes at a self-selected cadence. Following the warm-up, resistance was increased by 25-50W every two minutes until the participant could not maintain 60 revolutions per minute (RPM) on the cycle ergometer. During the VO2peak test, heart rate (Tango II, Suntech, Medical Inc, Morrisville, NC, USA), rating of perceived exertion (), respiratory exchange ratio (RER), and relative VO2 were recorded every 15s. The VO2peak was determined when at least two of the following physiological criteria were met: a) RER > 1.10, b) A failure of heart rate to increase with an increase in workload and/or c) RPE >17 at the final stage. VO2peak was considered the mean of the two highest 15-second values achieved during testing.
Experimental trial
Two experimental trials followed a repeated-measures, randomized crossover design. All exercise trials were conducted in a climate-controlled environmental chamber (University of Arkansas, Fayetteville, AR, USA). Photographs illustrating the climate-controlled environmental chamber, exercise equipment, and placement of the Kestrel 4400 are provided in Supplementary Figure 1. Participants completed matched exercise in two separate environmental conditions with target setpoints of HD (39.2 °C and 30% RH) and WH (34.2 °C and 55% RH), both designed to achieve an equivalent WBGT of approximately 30 °C. Environmental conditions were established and stabilized before participants entered the chamber and were maintained throughout the exercise protocol. Trials were separated by at least six days and occurred at the same time of day to account for circadian rhythm. On arrival, participants provided a 24-hour urine sample to ensure euhydration, defined as urine specific gravity ≤ 1.020 (USG, Master-SUR, Atago Co Ltd, Tokyo, Japan) (). One of our participants’ USG was >1.020, so we provided 500 mL of bottled water before the trial commenced. Nude body mass was assessed prior to trial instrumentation. Body mass with instrumentation was obtained when the participant mounted the cycle ergometer on a scale (VS-2501, WeighSouth.Inc, Asheville, NC, USA) before starting experimental trials. Participants then transitioned to the environmental chamber for a10 minutes seated baseline period. Participants then began a warm up for 2-5 minutes at a self-selected load, matched between trials. Participants then completed 60 minutes of cycling on a cycle ergometer (RacerMate, Seattle, WA, USA) at an individualized workload corresponding to 55% of the peak power output achieved during the graded VO2peak test. The same individualized absolute workload was maintained during the HD and WH trials, and power output and cycling cadence were monitored throughout each trial to ensure adherence to the prescribed external workload. During 60 minutes of cycling, rectal temperature (Tre), mean skin temperature (Tsk), and body mass with instrumentation were measured every 10 minutes. Blood pressure, heart rate, forearm skin blood flux (FSBF), and local skin temperature were measured every 12 minutes. WBGT, ambient temperature, RH, and wind speed were measured every 15 minutes. Watts and cadence were measured by VeloTron CS 2008 (Racer Mate, Seattle, WA). Participants were provided water (warmed to 39 °C) every 15 minutes to replace approximately 50% of estimated body mass losses. Partial fluid replacement was used to attenuate dehydration while avoiding complete fluid replacement, thereby minimizing the confounding effects of excessive dehydration or overhydration ().
After the 60 minutes of cycling, participants completed a 5 minutes cooldown on the cycle ergometer at a self-selected cadence and workload. Following the cooldown, participants got off the cycle ergometer and sat in a chair. Maximum skin blood flux was determined by heating forearm skin to 42°C for 30 minutes during recovery (; ; ). During a 30 minutes recovery period, participants were seated to have maximum skin blood flux measured via laser Doppler probe heating. After this, researchers removed instrumentation and measured nude body mass. Participants were provided refreshments, and future trials were scheduled, if necessary. A schematic overview of the experimental protocol, including the timeline, environmental conditions, exercise workload, physiological measurements, and principal instruments, is presented in Figure 1.
Figure 1
Instrumentation
Total body heat loss and heat gain parameters were calculated using partitional calorimetry as described by . Briefly, heat storage was calculated as metabolic heat production minus DHL, respiratory heat loss, and EHL. DHL was calculated as the sum of convective and radiative heat exchange, whereas conductive heat exchange was considered negligible. EHL was estimated from whole-body sweat rate, latent heat of vaporization, and sweating efficiency. Full equations (Equations 1–), assumptions, and constants used for these calculations are provided in the Appendix. Tre was determined by a rectal thermocouple (RET-1, Physitemp Instrumental Inc, Clifton, NJ, USA) inserted ~15 cm past the anal sphincter before participants started experimental trials. Tsk was continuously measured using wireless temperature sensors (iButton, Maxim Integrated, San Jose, CA, USA) attached with porous zinc oxide tape to the chest, triceps, thigh, and calf on the right side of the body. Four-site weighted mean skin temperature was calculated using the equation: Tsk = 0.3 (Tchest + Ttriceps) + 0.2 (Tthigh + Tcalf). This four-site method was selected because it is commonly used in exercise and environmental physiology research and provides a practical estimate of whole-body mean skin temperature while minimizing participant instrumentation burden during cycling exercise. Although forehead and neck sites may contribute meaningfully to regional heat exchange, these sites are not included in the original Ramanathan four-site equation; therefore, the current site selection was used to remain consistent with the established equation and facilitate comparison with prior studies. BP and heart rate were measured via an automated blood pressure cuff (Tango II, Suntech, Medical Inc, Morrisville, NC, USA) with a 3-lead ECG and microphone. Red blood cell flux was assessed by laser Doppler flowmetry on the left dorsal forearm with a probe (laser-Doppler perfusion monitor and Probe 2b, Moor Instruments, Wilmington, DE, USA) held in place by a local heater (Periflux System 5000, Perimed, Ardmore, PA, USA) (). Tre, Tsk, FSBF, and local skin temperature data were continuously recorded and stored using computer software for further analysis (LabChart 7; AD Instruments, Colorado Springs, CO, USA). Whole body sweat rate (WBSR) was calculated using the equation: WBSR = [(pre-body mass − post-body mass) − urine produced + fluid ingested]/exercise duration (). Ambient temperature, RH, WBGT, and wind speed were measured every 15 minutes throughout each experimental trial using a portable weather monitor (Kestrel 4400, Nielsen-Kellerman, Boothwyn, PA, USA). According to the manufacturer, the Kestrel 4400 has an accuracy of ±1.0 °C for ambient temperature, ± 3% for relative humidity, and ±3% of reading for wind speed measurements. The instrument was positioned adjacent to the cycle ergometer at approximately torso height and away from direct obstruction by the participant to verify the stability of environmental chamber conditions. Environmental conditions remained stable throughout the trials, as evidenced by the minimal variability in ambient temperature and relative humidity measurements presented in Table 2. The measurement ranges and manufacturer-reported accuracies of the principal instruments used during the experimental trials are summarized in Supplementary Table 1.
Table 2
| Variable | HD | WH |
|---|---|---|
| Temperature (°C) | 39.07 (0.21) | 34.46 (0.81) ** |
| RH (%) | 33.25 (0.81) | 58.79 (2.19) ** |
| WBGT | 30.1 (0.22) | 30.0 (0.27) |
| Wind Speed (m/s) Fluid Intake (L) | 2.44 (0.01) 0.82 (0.56) | 2.42 (0.05) 0.65 (0.37) |
| WBSR (L/h) | 1.68 (0.60) | 1.39 (0.36) |
| USG | 1.011 (0.005) | 1.012 (0.004) |
Measured environmental conditions used to verify chamber stability, including ambient temperature, relative humidity (RH), wet-bulb globe temperature (WBGT), and wind speed, as well as fluid intake, whole-body sweat rate (WBSR), and urine specific gravity (USG) during the experimental trials [n = 10, mean (SD)].
p < 0.05 between experimental conditions, **p < 0.01 between experimental conditions.
Statistical analysis
Data were analyzed with GraphPad Prism software (v.10.6.1). Descriptive statistics were calculated for all variables and are presented as mean ± standard deviation (SD). Linear mixed-effects models were used to assess differences in biophysical, thermal, and cardiovascular responses between HD and WH conditions over time, with condition and time treated as fixed effects and participant included as a random effect. Normality of residuals was assessed using the Shapiro–Wilk test. For variables that violated normality assumptions, non-parametric analyses were performed using the Friedman test to evaluate changes over time within each condition. When significant main effects were observed, post hoc pairwise comparisons were conducted using the Wilcoxon signed-rank test with Bonferroni adjustment. For parametric analyses, post hoc pairwise comparisons were performed with Bonferroni-adjusted alpha levels to identify significant differences between time points and conditions. Binary variables were analyzed using McNemar’s test. Statistical significance was set a priori at p ≤ 0.05.
Results
Table 2 presents the measured environmental conditions used to verify chamber stability, as well as WBSR, fluid intake, and USG for each experimental condition. As intended, ambient temperature (p < 0.001) and relative humidity (p < 0.001) differed significantly between HD and WH conditions. However, no differences were observed in WBGT (p = 0.566) or wind speed (p = 0.832) between conditions. Both HD and WH conditions were classified as moderate heat stress (Category 3; orange) based on WBGT criteria (). No significant differences were observed between conditions for WBSR (p = 0.176), fluid intake (p = 0.127), or USG (p = 0.756).
Biophysical, thermal, and cardiovascular responses for the HD and WH conditions, together with the corresponding condition-effect statistical results and effect size are summarized in Table 3.
Table 3
| Variable | n | HD, Mean (SD) | WH, Mean (SD) | Condition effect (p) | Effect size |
|---|---|---|---|---|---|
| DHL (W) | 10 | −80 (25) | 48 (25) | < 0.001 | ηp² = 0.950 |
| EHL (W) | 10 | 450 (53) | 319 (30) | < 0.001 | ηp² = 0.932 |
| Heat storage (W) | 10 | 285 (115) | 245 (81) | 0.323 | ηp² = 0.347 |
| Cumulative heat storage (kJ) | 10 | 1058 (350) | 977 (288) | 0.364 | d = 0.300 |
| Tre (°C) | 10 | 38.31 (0.79) | 38.29 (0.71) | 0.885 | ηp² = 0.009 |
| Tsk (°C) | 8 | 36.28 (0.67) | 35.89 (0.70) | 0.030 | ηp² = 0.465 |
| Heart rate (bpm) | 10 | 146 (40) | 143 (40) | 0.690 | ηp² = 0.147 |
| SBP (mmHg) | 9 | 166 (39) | 164 (38) | 0.876 | ηp² < 0.001 |
| DBP (mmHg) | 9 | 62 (17) | 61 (16) | 0.737 | ηp² = 0.143 |
| MAP (mmHg) | 9 | 93 (17) | 96 (13) | 0.926 | ηp² = 0.121 |
| FSBF (AU) | 10 | 113 (59) | 95 (44) | 0.324 | ηp² = 0.001 |
| %CVCmax (AU/mmHg) | 10 | 0.65 (0.35) | 0.61 (0.27) | 0.646 | ηp² < 0.001 |
Comparison of biophysical, thermal, and cardiovascular responses between hot-dry (HD) and warm-humid (WH) conditions.
Hot-dry; WH, warm-humid; DHL, dry heat loss; EHL, evaporative heat loss; Tre, rectal temperature; Tsk, mean skin temperature; SBP, systolic blood pressure; DBP, diastolic blood pressure; MAP, mean arterial pressure; FSBF, forearm skin blood flux; CVCmax, maximal cutaneous vascular conductance; AU, arbitrary units; ηp², partial eta squared; d, Cohen’s d. Effect sizes are reported as partial eta squared (ηp²) for outcomes analyzed using linear mixed-effects models and Cohen’s d for cumulative heat storage. Statistical significance was set at p ≤ 0.05.
Biophysical responses
DHL was significantly greater in WH than HD (p < 0.001), with no main effect of time (p = 0.438) or condition × time interaction (p = 0.130). DHL remained negative throughout the exercise protocol in HD, indicating net dry heat gain, whereas positive values in WH indicated net dry heat loss (Figure 2A; Table 3). Conversely, EHL was significantly greater in HD than WH (p < 0.001), with no main effect of time (p = 0.170) or condition × time interaction (p = 0.523). Heat storage did not differ over time (p = 0.116), between conditions (p = 0.323), or for the condition × time interaction (p = 0.430; Figure 2B; Table 3). Similarly, cumulative heat storage did not differ between conditions (p = 0.364; Figure 2C; Table 3).
Figure 2
Thermal responses
Tsk data were unavailable for two participants because the iButton sensors detached during exercise or the recorded temperature data could not be retrieved. Therefore, analyses involving Tsk included eight participants. Tre increased significantly over time (p < 0.001), with no main effect of condition (p = 0.885) or condition × time interaction (p = 0.782; Figure 3A; Table 3). Tsk demonstrated significant main effects of time (p = 0.033) and condition (p = 0.030), with greater values in HD than WH, but no condition × time interaction (p = 0.299; Figure 3B; Table 3).
Figure 3
Cardiovascular responses
Heart rate increased significantly over time (p < 0.001), with no main effect of condition (p = 0.244) or condition × time interaction (p = 0.514; Figure 4A; Table 3). SBP and DBP also increased over time (both p < 0.001), whereas MAP demonstrated a significant time effect (p = 0.040). No condition effects or condition × time interactions were observed for SBP, DBP, or MAP (all p > 0.05; Figure 4B; Table 3). FSBF and %CVCmax increased significantly over time (both p < 0.001), with no condition effects or condition × time interactions (all p > 0.05; Figures 4C, D; Table 3).
Figure 4
Discussion
The current study was designed to compare heat storage and physiological responses to different environmental conditions (HD versus WH) but similar WBGT during moderate intensity cycling. Our primary findings support similar thermoregulation responses to previous literature based on minimal or light exercise. With moderate intensity exercise, our data indicate no differences in heat storage, Tre, heart rate, SBP, DBP, MAP, FSBF, and %CVCmax between HD and WH conditions. However, EHL and Tsk in HD were significantly greater than in WH, while DHL in WH was significantly greater than HD. Hence, the results of the present study demonstrate that moderate heat stress of varying conditions, but similar WBGT, stimulated matched cardiovascular responses, yet altered thermoregulation to accommodate moderate intensity 60 minutes exercise in the heat. Likely, convective heat losses were increased during WH conditions to allow thermoregulatory similarities, despite our uncompensable heat stress. Our findings suggest that WBGT alone may not fully characterize the physiological mechanisms of heat exchange across different environmental conditions, particularly during moderate-intensity exercise. Occupational, military, or athletic groups should be cautious when following any guidelines, recommendations, standards, or protocols for heat stress based on the WBGT assessment, if work rate is not taken into consideration for activity modifications.
The primary mechanism of heat dissipation in the current study can be explained by differences in DHL and EHL between HD and WH conditions. DHL in HD was negative, which revealed heat gain compared to our WH condition. In the present study of both HD and WH environmental conditions, heat balance was mainly determined by convection and evaporation. Under our HD condition, the ambient temperature was greater than Tsk (Table 2), and heat transferred from the environment to the body. In addition, fan use increases ambient air movement, which further increases heat gain via convection when there is a gradient from the environment to the skin (). According to previous research on ambient temperature thresholds where fan use is beneficial, using a fan is an effective strategy for cooling the human body when air temperature is < 39°C (). In contrast, convective heat loss is enhanced when the air temperature is lower than Tsk; as seen in our WH condition. In addition, fan use in our WH condition (34.5°C) boosts convection due to ambient temperature < 39°C. Regarding EHL in our WH condition, high relative humidity decreases evaporative efficacy (). A recent study also demonstrated that sweat can be fully evaporated at low relative humidity (< 35%), but evaporation of sweat is restricted when relative humidity is > 55% (). Under our WH condition, EHL was limited due to high relative humidity (58.79%). In contrast, sweat can be fully evaporated in our HD condition because of lower relative humidity (33.25%). In conclusion, our data suggest that heat balance in HD is limited by convection rather than evaporation of sweat, whereas heat balance in WH was limited by the evaporation of sweat.
No difference in heat storage and cumulative heat storage between HD and WH conditions indicated that heat loss from evaporation in HD was partly offset by the heat gain from convection and radiation, consistent with previous studies with WBGT 27.8–30°C (; ). However, found heat storage in HD (ambient temperature: 40°C, RH: 30%, WBGT: 32°C) was greater compared to WH (ambient temperature: 31°C, RH: 80%, WBGT: 32°C) conditions during 60 minutes of cycling at 40% V O2peak. The divergent results from previous studies indicated that there may be a threshold near 32°C WBGT before the difference in heat storage is identified.
Consistent with our heat storage findings, an equivalent WBGT ~30°C resulted in a similar Tre between HD and WH conditions during 60 minutes of moderate-intensity cycling. Previous studies agree with our findings with matched WBGT (27.8–30°C) of combining varying ambient temperature (20–46°C) and relative humidity (10–84%) (; ; ; ). For example, one of the studies used similar WBGT (30°C) and exercise type (cycling) to the present investigation (). Previous data suggest Tre was no different between HD and WH conditions during 60 minutes of cycling at a fixed work rate of 350 watts. In contrast, reported that core temperature in HD (ambient temperature: 40°C, RH: 30%) was significantly greater than WH (ambient temperature: 31°C, RH: 80%) at equivalent WBGT during 60 minutes of cycling at 40% V O2peak. More extreme heat stress in that study may explain the differences in core temperature between HD and WH conditions (). Future research should investigate whether there is a threshold of WBGT before differences may be observed between HD and WH conditions.
Tsk reflects temperature changes in the peripheral shell of the human body. In the current study, Tsk was greater in HD than WH, consistent with previous studies comparing HD and WH conditions with equivalent WBGT values ranging from 27.8 to 32 °C (; ; ; ). The greater Tsk observed in HD was likely attributable, in part, to the higher ambient temperature in HD (~39 °C) than WH (~34.5 °C). In addition, greater dry heat gain occurred in HD than WH. This transfer of dry heat from the environment to the skin may have contributed to the higher Tsk observed in HD. Greater interindividual variability in Tsk was also evident in the HD condition, including at the initial measurement, and remained apparent throughout exercise. Because this variability was present at the beginning of the trial, it may not be attributable solely to subsequent exercise-induced thermal responses and may partly reflect individual differences in regional skin temperature or local factors associated with skin-temperature measurement. Furthermore, Tsk analyses included eight participants because data from two participants were unavailable due to iButton sensor detachment or unreadable recordings. Therefore, the condition effect for Tsk should be interpreted cautiously given the small sample size and observed interindividual variability.
Integrative cardiovascular responses to heat stress during prolonged exercise depends on thermal stress and exercise intensity (). In the current study, we found no cardiovascular differences between conditions. Heart rate is widely used to assess physiological strain in combination with core temperature during exercise in the heat (). During exercise in the heat, more blood is distributed to the skin for heat dissipation. Heart rate increases to maintain cardiac output and meet the demands for maintaining blood pressure, heat dissipation, and exercising muscle perfusion (). Our data confirmed that heart rate increased throughout 60 minutes of cycling in both HD and WH conditions. In contrast, there was no difference in HR between HD and WH, which is consistent with findings from previous studies (; ; ). Heart rate in the current study reached over 86% age-predicted maximum heart rate after 30 minutes of cycling (Figure 3), which corresponds to vigorous intensity exercise (). The heart rate data demonstrated that the exercise intensity ranged from moderate to vigorous during our trials. In the current study, SBP and MAP increased and reached a maximum around 30 minutes of cycling and were maintained throughout the protocol. In contrast, DBP was maintained during the 60 minutes of cycling.
Skin blood flow increases substantially in response to thermal stress via vasodilation (). In the present study, FSBF increased throughout the 60 minutes of cycling. Although Tsk in HD condition was higher than WH, FSBF was no different between HD and WH conditions. The amplitude of change in the skin blood flux during exercise in the heat is mediated by the core-to-skin temperature gradient (; ). The rate of change in skin blood flux is a priority driven by central (i.e., visceral and brain) thermoreceptors at any temperature gradient compared to peripheral (i.e., skin) thermoreceptors (). Our data revealed that Tre was not significantly different between HD and WH conditions, which likely accounted for similar responses in skin blood flux between HD and WH conditions, despite Tsk differences. Altogether, under exercise-heat stress, the cardiovascular system dynamically adjusted and balanced the competition between skin and muscle for large fractions of cardiac output ().
Examining the effectiveness of heat stress indices in differentiating the severity of environmental heat stress is important for developing guidelines and strategies to mitigate heat-related health consequences. Both HD and WH conditions in the present study were classified as moderate heat stress (Category 3; orange) according to WBGT criteria (), for which recommendations permit up to 120 minutes of continuous activity under these environmental conditions, depending on work intensity and other modifying factors. However, these recommendations may not fully account for differences in metabolic demands and environmental characteristics. Although substantial increases in rectal temperature occurred during the 60 minutes exercise protocol, the present data do not permit accurate prediction of physiological responses beyond the experimental duration. Therefore, caution should be exercised when extrapolating these findings to longer periods of exercise or work in the heat. Nevertheless, our findings suggest that both environmental characteristics and metabolic demands should be considered when applying WBGT-based heat stress recommendations during moderate-intensity exercise.
Our findings should be taken with caution based on limitations. First, it would be rare to continue steady-state physical activity or work in living situations when people experience prolonged heat stress. Our experimental design did not allow for heat stress mitigating strategies, such as reducing workload, rest breaks or slowing down. Second, the present study was limited to 60 minutes of exercise and therefore cannot characterize the longer-term trajectory of core temperature responses in either environmental condition. Future studies employing longer exercise durations are needed to determine whether thermoregulatory responses eventually diverge between hot-dry and warm-humid environments. Third, mean skin temperature was estimated using the Ramanathan four-site equation. Although this method is widely used, it may not fully capture regional skin temperature differences at sites not measured in the present study, such as the forehead or neck, which may contribute to heat exchange during exercise in the heat. Fourth, although the randomized crossover design reduced between-subject variability by allowing each participant to serve as his own control, the exclusion of five participants resulted in a final sample of ten participants, which may have reduced statistical power to detect smaller differences in secondary outcomes and limits the generalizability of the findings. Furthermore, an a priori power analysis was not performed, and therefore the possibility of a type II error for some outcomes cannot be excluded. Fifth, environmental conditions were monitored using a portable weather monitor rather than a fixed laboratory-grade environmental monitoring system. Although the Kestrel 4400 may provide lower precision than laboratory-grade instruments, the small variability in measured environmental conditions suggests that this limitation likely had little influence on the study findings. Finally, our data were collected in young, physically fit men in whom thermoregulatory capacity is likely preserved. The relatively homogeneous aerobic fitness of the participants (all VO2peak ≥ 45 mL/kg/min), combined with the small sample size, limited the suitability of exploratory analyses examining whether VO2peak or body fat percentage influenced the observed responses. Future studies with larger samples and broader ranges of aerobic fitness and body composition are needed to determine whether these characteristics modify responses under HD and WH conditions. Consequently, these findings may not be generalizable to women, older adults, children, or other vulnerable populations. Females may demonstrate different thermoregulatory responses due to differences in body composition, sweating characteristics, skin blood flow, and hormonal fluctuations throughout the menstrual cycle. Future studies should investigate whether the present findings are consistent across these populations.
Conclusion
The findings of our study provide important insight regarding heat storage, thermoregulatory and cardiovascular responses to heat stress of varying environmental conditions with equivalent WBGT during moderate-intensity exercise. Our data reveal hot environments of equivalent WBGT alter thermoregulation to accommodate for differences in humidity, which confirmed previous findings and extend data to include more extreme environments and greater exercise intensity. It is noteworthy that our data support WBGT as a reliable heat stress indicator for assessing thermal and cardiovascular strain in varying environmental conditions during moderate intensity exercise. However, these findings suggest that WBGT alone may not fully capture physiological strain under all combinations of environmental conditions and metabolic demands, and caution should be taken when occupational, military, or athletic groups use guidelines, recommendations, standards, or protocols for heat stress based on the WBGT assessment, unless work rate is considered. Future investigations should examine similar heat stress in the aged, women, children, or other vulnerable populations.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by University Institutional Review Board at University of Arkansas. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
XZ: Formal analysis, Data curation, Conceptualization, Writing – original draft, Methodology. BM: Writing – review & editing, Methodology, Data curation, Supervision, Conceptualization. CB: Data curation, Methodology, Supervision, Writing – review & editing.
Funding
The author(s) declared that financial support was not received for this work and/or its publication.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2026.1856143/full#supplementary-material
References
1
Alber-WallerströmB.HolmérI. (1985). Efficiency of sweat evaporation in unacclimatized man working in a hot humid environment. Eur. J. Appl. Physiol. Occup. Physiol.54, 480–487. doi: 10.1007/BF00422956
2
ArmstrongL. E. (2007). Assessing hydration status: the elusive gold standard. J. Am. Coll. Nutr.26, 575S–584S. doi: 10.1080/07315724.2007.10719661
3
BartmanN. E.VargasN. T.CavuotoL. A.HostlerD.PryorR. R. (2024). Heat strain differences walking in hot-dry and warm-wet environments of equivalent wet bulb globe temperature. Temp. (Austin).11, 333–349. doi: 10.1080/23328940.2024.2384185
4
BeigtanM.GonçalvesM.WeonB. M. (2024). Heat transfer by sweat droplet evaporation. Environ. Sci. Technol.58, 6532–6539. doi: 10.1021/acs.est.4c00850
5
BorgG. (1970). Perceived exertion as an indicator of somatic stress. Scand. J. Rehabil. Med.2, 92–98. doi: 10.2340/1650197719702239298
6
BuddG. M. (2008). Wet-bulb globe temperature (WBGT)--its history and its limitations. J. Sci. Med. Sport.11, 20–32. doi: 10.1016/j.jsams.2007.07.003
7
CharkoudianN. (2003). Skin blood flow in adult human thermoregulation: how it works, when it does not, and why. Mayo. Clin. Proc.78, 603–612. doi: 10.4065/78.5.603
8
CharkoudianN.StachenfeldN. (2016). Sex hormone effects on autonomic mechanisms of thermoregulation in humans. Auton. Neurosci.196, 75–80. doi: 10.1016/j.autneu.2015.11.004
9
ChouT.-H.CoyleE. F. (2023). Cardiovascular responses to hot skin at rest and during exercise. Temp. (Austin).10, 326–357. doi: 10.1080/23328940.2022.2109931
10
ConsolazioC. F.JohnsonE. R.FecoraJ. L. (1963). Physiological Measurements of Metabolic Functions in Man (New York: Blakiston Division, McGraw-Hill).
11
CottleR. M.VecellioD. J.WolfS. T.KenneyW. L. (2022). Core temperature responses to compensable versus uncompensable heat stress in young adults (PSU HEAT Project). J. Appl. Physiol. (1985)133, 1011–1018. doi: 10.1152/japplphysiol.00388.2022
12
CramerM. N.JayO. (2019). Partitional calorimetry. J. Appl. Physiol. (1985)126, 267–277. doi: 10.1152/japplphysiol.00191.2018
13
Du BoisD.Du BoisE. F. (1916). Clinical calorimetry: tenth paper a formula to estimate the approximate surface area if height and weight be known. Arch. Intern. Med. (Chic).XVII, 863–871. doi: 10.1001/archinte.1916.00080130010002
14
EbiK. L.CaponA.BerryP.BroderickC.de DearR.HavenithG.et al. (2021). Hot weather and heat extremes: health risks. Lancet398, 698–708. doi: 10.1016/S0140-6736(21)01208-3
15
FeitoY.MagalM. (Eds.) (2022). “ American college of sports medicine,” in Acsm’s Fitness Assessment Manual, Sixth edition ( Wolters Kluwer, Philadelphia Baltimore New York London Buenos Aires Hong Kong Sydney Tokyo).
16
FosterJ.SmallcombeJ. W.HodderS.JayO.FlourisA. D.NyboL.et al. (2022). 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.66, 175–188. doi: 10.1007/s00484-021-02205-x
17
GibsonO. R.JamesC. A.MeeJ. A.WillmottA. G. B.TurnerG.HayesM.et al. (2020). Heat alleviation strategies for athletic performance: a review and practitioner guidelines. Temp. (Austin).7, 3–36. doi: 10.1080/23328940.2019.1666624
18
GrundsteinA.WilliamsC.PhanM.CooperE. (2015). Regional heat safety thresholds for athletics in the contiguous United States. Appl. Geogr.56, 55–60. doi: 10.1016/j.apgeog.2014.10.014
19
GubernotD. M.AndersonG. B.HuntingK. L. (2014). The epidemiology of occupational heat exposure in the United States: a review of the literature and assessment of research needs in a changing climate. Int. J. Biometeorol.58, 1779–1788. doi: 10.1007/s00484-013-0752-x
20
HayesM.CastleP. C.RossE. Z.MaxwellN. S. (2014). The influence of hot humid and hot dry environments on intermittent-sprint exercise performance. Int. J. Sports. Physiol. Perform.9, 387–396. doi: 10.1123/ijspp.2012-0247
21
International Organization for Standardization (2007). Ergonomics of the Thermal Environment — Estimation of Thermal Insulation and Water Vapour Resistance of a Clothing Ensemble (Geneva: International Organization for Standardization). Available online at: https://www.iso.org/standard/39257.html (Accessed March 20, 2026).
22
JohnsonJ. M.O’LearyD. S.TaylorW. F.KosibaW. (1986). Effect of local warming on forearm reactive hyperaemia. Clin. Physiol.6, 337–346. doi: 10.1111/j.1475-097x.1986.tb00239.x
23
JohnsonJ. M.ParkM. K. (1979). Reflex control of skin blood flow by skin temperature: role of core temperature. J. Appl. Physiol. Respir. Environ. Exerc. Physiol.47, 1188–1193. doi: 10.1152/jappl.1979.47.6.1188
24
KamonE.AvelliniB. (1976). Physiologic limits to work in the heat and evaporative coefficient for women. J. Appl. Physiol.41, 71–76. doi: 10.1152/jappl.1976.41.1.71
25
KamonE.AvelliniB.KrajewskiJ. (1978). Physiological and biophysical limits to work in the heat for clothed men and women. J. Appl. Physiol. Respir. Environ. Exerc. Physiol.44, 918–925. doi: 10.1152/jappl.1978.44.6.918
26
KeatisuwanW.OhnakaT.TochiharaY. (1996). Physiological responses of men and women during exercise in hot environments with equivalent WBGT. Appl. Hum. Sci.15, 249–258. doi: 10.2114/jpa.15.249
27
KelloggD. L.LiuY.KosibaI. F.O’DonnellD. (1999). Role of nitric oxide in the vascular effects of local warming of the skin in humans. J. Appl. Physiol. (1985)86, 1185–1190. doi: 10.1152/jappl.1999.86.4.1185
28
KenneyW. L.ZemanM. J. (2002). Psychrometric limits and critical evaporative coefficients for unacclimated men and women. J. Appl. Physiol. (1985)92, 2256–2263. doi: 10.1152/japplphysiol.01040.2001
29
LotensW. A.HavenithG. (1991). Calculation of clothing insulation and vapour resistance. Ergonomics34, 233–254. doi: 10.1080/00140139108967309
30
McculloughE. A.JonesB.HuckJ. (1985). A Comprehensive Data Base for Estimatng Clothing Insulation ( Ashrae Transactions). Available online at: https://www.semanticscholar.org/paper/A-comprehensive-data-base-for-estimatng-clothing-Mccullough-Jones/3d1c0ff407191580efbc9d031fa21435e3626fee (Accessed June 24, 2026).
31
McDermottB. P.AndersonS. A.ArmstrongL. E.CasaD. J.CheuvrontS. N.CooperL.et al. (2017). National athletic trainers’ Association position statement: fluid replacement for the physically active. J. Athl. Train.52, 877–895. doi: 10.4085/1062-6050-52.9.02
32
MoraC.DoussetB.CaldwellI. R.PowellF. E.GeronimoR. C.BieleckiC. R.et al. (2017). Global risk of deadly heat. Nat. Clim. Change7, 501–506. doi: 10.1038/nclimate3322
33
MoranD. S.DeGrootD. W.PotterA. W.CharkoudianN. (2023). Beating the heat: military training and operations in the era of global warming. J. Appl. Physiol. (1985)135, 60–67. doi: 10.1152/japplphysiol.00229.2023
34
MorrisN. B.ChaselingG. K.EnglishT.GrussF.MaideenM. F. B.CaponA.et al. (2021). Electric fan use for cooling during hot weather: a biophysical modelling study. Lancet Planet. Health5, e368–e377. doi: 10.1016/S2542-5196(21)00136-4
35
NilssonG. E.TenlandT.ObergP. A. (1980). Evaluation of a laser Doppler flowmeter for measurement of tissue blood flow. IEEE Trans. Biomed. Eng.27, 597–604. doi: 10.1109/TBME.1980.326582
36
PérgolaP. E.JohnsonJ. M.KelloggD. L.KosibaW. A. (1996). Control of skin blood flow by whole body and local skin cooling in exercising humans. Am. J. Physiol.270, H208–H215. doi: 10.1152/ajpheart.1996.270.1.H208
37
PérgolaP. E.KelloggD. L.JohnsonJ. M.KosibaW. A.SolomonD. E. (1993). Role of sympathetic nerves in the vascular effects of local temperature in human forearm skin. Am. J. Physiol.265, H785–H792. doi: 10.1152/ajpheart.1993.265.3.H785
38
PériardJ. D.EijsvogelsT. M. H.DaanenH. A. M. (2021). Exercise under heat stress: thermoregulation, hydration, performance implications, and mitigation strategies. Physiol. Rev.101, 1873–1979. doi: 10.1152/physrev.00038.2020
39
Perkins-KirkpatrickS. E.LewisS. C. (2020). Increasing trends in regional heatwaves. Nat. Commun.11, 3357. doi: 10.1038/s41467-020-16970-7
40
RamanathanN. L. (1964). A New Weighting System for Mean Surface Temperature of the Human Body. J Appl Physiol19, 531–533. doi: 10.1152/jappl.1964.19.3.531
41
RavanelliN.CoombsG. B.ImbeaultP.JayO. (2018). Maximum skin wettedness after aerobic training with and without heat acclimation. Med. Sci. Sports. Exercise50, 299–307. doi: 10.1249/MSS.0000000000001439
42
RobineJ.-M.CheungS. L. K.Le RoyS.Van OyenH.GriffithsC.MichelJ.-P.et al. (2008). Death toll exceeded 70,000 in Europe during the summer of 2003. C.R. Biol.331, 171–178. doi: 10.1016/j.crvi.2007.12.001
43
RowellL. B. (1993). Human cardiovascular control.New York, NY: Oxford University Press. pp. 204–254.
44
SeoY.PowellJ.StrauchA.RobergeR.KennyG. P.KimJ.-H. (2019). Heat stress assessment during intermittent work under different environmental conditions and clothing combinations of effective wet bulb globe temperature (WBGT). J. Occup. Environ. Hyg.16, 467–476. doi: 10.1080/15459624.2019.1612523
45
SmolanderJ.IlmarinenR.KorhonenO.PyykköI. (1987). Circulatory and thermal responses of men with different training status to prolonged physical work in dry and humid heat. Scand. J. Work. Environ. Health13, 37–46. doi: 10.5271/sjweh.2083
46
StapletonJ. M.WrightH. E.HardcastleS. G.KennyG. P. (2012). Body heat storage during intermittent work in hot-dry and warm-wet environments. Appl. Physiol. Nutr. Metab.37, 840–849. doi: 10.1139/h2012-053
47
TanabeS.NaritaC.OzekiY.KonishiM. (2000). Effective radiation area of human body calculated by a numerical simulation. Energy Build.32, 205–215. doi: 10.1016/S0378-7788(00)00045-1
48
VargasN. T.ChapmanC. L.JohnsonB. D.GathercoleR.CramerM. N.SchladerZ. J. (2019). Thermal behavior alleviates thermal discomfort during steady-state exercise without affecting whole body heat loss. J. Appl. Physiol.127, 984–994. doi: 10.1152/japplphysiol.00379.2019
49
WinslowC.-E. A.HerringtonL. P.GaggeA. P. (1936). The determination of radiation and convection exchanges by partitional calorimetry. Am. J. Physiology-Legacy. Content.116, 669–684. doi: 10.1152/ajplegacy.1936.116.3.669
50
WolfS. T.FolkertsM. A.CottleR. M.DaanenH. A. M.KenneyW. L. (2021). Metabolism- and sex-dependent critical WBGT limits at rest and during exercise in the heat. Am. J. Physiol. Regul. Integr. Comp. Physiol.321, R295–R302. doi: 10.1152/ajpregu.00101.2021
Appendix
Metabolic heat production calculations
Metabolic heat production (Hprod) was estimated with the equation ():
where M is metabolic energy expenditure and Wk is external work rate. M was calculated using the equation ():
VO2 and CO2 were measured from expired gas for each subject (see above for procedure), and RER was calculated using VO2 and VCO2:
Wk was calculated using the 55% of wattage reached during the last stage of the V• O2peak test.
Rate of heat storage calculations
Heat storage (S) was calculated at 30 and 60 min throughout the protocol, calculated as:
where Hdry represents the sum of dry heat losses from radiation (R) and convection (C), Hres represents heat loss from respiration, and EHL represents heat loss from evaporation. The rate of heat loss from conduction was considered negligible and thus was eliminated from the equation ().
Total Hdry was calculated as the sum of C+R where C Equation 5 and R Equation 6 were calculated as:
where To is the operational temperature, Tsk is mean skin temperature, and hc and hr are the convective heat transfer coefficient and radiative heat transfer coefficient, respectively (). Since radiant heat was considered negligible, ambient temperature (Ta) was used instead of operational temperature (To): To=Ta.
Hc value for cycle ergometry was calculated using the equation ():
Where fped is pedalling frequency (revolutions/min), Vair is air velocity (m/s).
Hr was calculated as ():
where ϵ is the nondimensional emissivity of the body surface, assumed 0.95 for clothed skin, σ is the Stefan-Boltzmann constant (5.67 × 10−8 W/m2/K4), Ar/AD was assumed 0.70 for a sitting individual (; ).
BSA was calculated as:
Whole body sweat rate (WBSR) was calculated as:
where Mpre is the initial mass determination, Mpost is the final mass determination, and 60 is the total minutes of the protocol.
Hres was calculated using M as:
where Pa is the partial pressure of water vapor in ambient air as kPa calculated as ():
where Pa,s is the saturated water vapor pressure of the air:
EHL was estimated using measurements of WBSR Equation (1) and sweating efficiency (Seff) with the equation:
Where λ is the latent heat vaporization of sweat in (2,426 J/g), and 60 represents the conversion from minutes to seconds for calculating evaporative heat loss in watts.
Sweating efficiency was calculated as:
Skin wettedness required to achieve heat balance (ω) was calculated as:
where ωmax was assumed 0.85 for an aerobically trained, unacclimated subject (). Psk,s is water vapor pressure at the skin when saturated, Re,cl is the evaporative heat transfer resistance of clothing ensemble assumed 0.005, similar to measured Re,cl for an athletic ensemble of shorts, t shirt, socks, and sneakers (), he is the evaporative heat transfer coefficient of air adjacent to the clothing layer, fcl is the clothing area factor (1.031 ND), and AD is body surface area calculated from the standard DuBois equation ().
he was calculated as:
where Pb is barometric pressure in mmHg.
fcl was calculated as:
where 0.1 clo is the insulation factor (Icl) of subject attire ().
Cumulative S was calculated as the sum of total heat storage in kJ determined by the S (in W) at each time-point (30 and 60 min):
Where s is the interval (in second) between the previous time-point. The summation of all two time-points was cumulative S.
Summary
Keywords
biophysical, cardiovascular response, exercise, heat balance, thermal strain
Citation
Zhao X, McDermott BP and Butts CL (2026) Thermoregulation in hot-dry and warm-humid heat stress during cycling among males. Front. Physiol. 17:1856143. doi: 10.3389/fphys.2026.1856143
Received
15 April 2026
Revised
10 July 2026
Accepted
21 July 2026
Published
05 August 2026
Volume
17 - 2026
Edited by
Zbigniew Waśkiewicz, Jerzy Kukuczka Academy of Physical Education in Katowice, Poland
Reviewed by
Shu Sun, University of Wisconsin-Madison, United States
Barlo Hillen, Johannes Gutenberg University Mainz, Germany
Updates
Copyright
© 2026 Zhao, McDermott and Butts.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Brendon P. McDermott, brendonm@uark.edu
†ORCID: Xiujing Zhao, orcid.org/0009-0001-3064-8460; Brendon P. McDermott, orcid.org/0000-0001-7238-8145; Cory L. Butts, orcid.org/0000-0002-8720-5206
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.