Abstract
Introduction
The aim of this study was to analyze changes in heart rate variability (HRV) during exercise in hot environments and recovery to baseline values depending on relative humidity.
Methods
Ten recreational runners participated in this study. Each participant performed 2 trials consisting of 30 min of continuous running under hot and dry (HD) (38°C and 28% relative humidity) and hot and humid (HH) conditions (38°C and 64% relative humidity) at their common 10 km race-running rhythm. HRV and body mass were assessed pre- and post-trial; the rating of perceived exertion and HRV were assessed during the trial; and HRV measurements were repeated 2, 4, 8, and 24 h postexercise. Primary HRV outcomes were root mean square of the successive differences, high frequency power, stress score, and sympathetic/parasympathetic ratio. One-way analysis of variance testing was used to analyze differences.
Results
No significant difference in body mass occurred across the different conditions or distances covered (P>0.05). Rating of perceived exertion presented the highest correlation values with stress score (r=0.729 for HD; r=0.568 for HH) and sympathetic/parasympathetic ratio (r=0.621 for HD; r=0.519 for HH) during exercise. HRV recovered to baseline values more quickly after exercising under dry conditions (4 h) than under humid conditions (between 8 and 24 h).
Conclusions
Stress score and sympathetic/parasympathetic ratio seem to be the best HRV predictors of internal load. Although there are no differences in HRV during recovery at the same time points in both conditions, the recovery is slower after exercise in HH than in HD.
Introduction
Heart rate (HR) variability (HRV) has been commonly used to assess autonomic nervous system (ANS) modulation. Variables of the time domain calculated from R-R intervals (RR), such as root mean square of successive differences of RR (RMSSD) or consecutive intervals that differ by more than 50 ms (pNN50) are related to parasympathetic activity, and increased values of these variables are related to physical training adaptations.1,2 Additionally, frequency domain variables, including very low frequency, low frequency (LF), and high frequency (HF) power, may be predictors of physical fitness. 3 Poincare plot analysis of HRV is also a common practice for assessment of the ANS because the standard deviation of instantaneous beat-to-beat RR interval variability (SD1) is identical to RMSSD. 4 Likewise, Poincare plot variables have been used to propose new indexes for the assessment of ANS balance. 5
HRV has become an interesting tool for training adaptation control. 6 –10 In regard to endurance exercise, there is a significant decrease in HRV during exercise for athletes with higher values of maximal oxygen consumption (VO2max) than for athletes with lower VO2max levels.11,12 Previous research has used HRV to assess ventilatory thresholds, 13 demonstrating the relationship between HRV and endurance exercise intensity. During exercise, both time and frequency domain variables decrease compared to resting values.14,15 HRV values start to recover immediately after exercise, but baseline levels are not reached 10 min 15 nor 60 min after exercise. 16 The intensity of the exercise has shown to be a determinant in the time needed to reach full recovery. 17 This means that the higher the exercise intensity, the more time it will take to return to baseline values.18,19 Nevertheless, exercise duration is not directly related to the recovery of HRV because different durations at the same intensity do not result in different recovery times.18,20
When performing resting measurements, heat has an important effect on HRV, with lower values in RMSSD and in the standard deviation of RR intervals, 21 as well as in pNN50 22 for hot conditions. Additionally, hot environments have a direct effect on frequency domain variables, such as LF and HF power, 23 –26 total power (TP), and very low frequency power, 22 –24 showing that higher temperatures provoke a decrease in cardiac variability. Higher relative humidity (RH) shows increments in body, skin, and rectal temperature and HR during exercise under hot and neutral conditions27,28 but not in the rate of perceived exertion. 28 The effect of RH on HRV has not been analyzed before; thus, it might be important to describe the effect RH may have on HRV.
Therefore, the aims of this study were first to analyze changes in HRV during exercise in hot environments and second to analyze recovery to baseline values depending on RH.
Methods
PARTICIPANTS
The present study was approved by the San Antonio University of Murcia ethics committee. Ten healthy recreational runners (28±5 y old, 172±9 cm tall, 65±10 kg, and running 39±5 km·wk-1) voluntarily participated in this study. All participants signed an informed consent form before the start of the measurements. Those under pharmacological or medical treatment were excluded from the study.
PROTOCOL
Each participant performed 2 trials on 2 different days separated by 1 week. Each trial was performed at the same time of the day (1400) in 2 different environmental conditions: a hot and dry environment (HD) with a temperature of 38°C and 28% RH and a hot and humid (HH) environment at 38°C and 64% RH. The orders of the trials were randomly assigned for every subject. Subjects were asked not to exercise or to consume any alcohol or caffeine 24 h before and after both tests. Before the first trial, anthropometric measurements were assessed; body mass and stretch stature were also assessed using a balance with stadiometer (Seca 720; Hamburg, Germany). Ambient temperature and RH were assessed using an environmental meter Kestrel 5000 (Minneapolis, MN).
Before trials, a 10-min warm-up was carried out, consisting of lower body dynamic stretching exercises. The trials consisted of 30 min of running in both environments on a running track. For both of the trials, the athletes were asked to run at the speed they would use during a 10 km race and were allowed to drink water ad libitum; water intake was recorded. No specific intensity was set in order to extrapolate the results of recreational running practitioners; average speed was calculated at the end of the trial. BM was assessed before and after the test but not after posterior recovery. Every 5 min, the participants were asked for a rating of perceived exertion (RPE) 29 score from 6 to 20. HRV was assessed before each trial (Pre); during the trial; immediately after the test (Post); and 2, 4, 8, and 24 h after the end of the trial. Participants were asked not to practice intense exercise until the 24 h measurement was conducted. Pre- and post-trial and recovery measurements were carried out at 25°C. HRV measurements were performed with the participants lying down, except during the exercise. Subjects wore a HR sensor Polar H7 (Kempele, Finland) to assess beat-to-beat recordings during the measurements. Data were analyzed during 5-min intervals according to the procedures determined by the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 30 for a total 6 HRV assessments during exercise.
HRV ANALYSIS
Further examinations of HRV variables were done using the software Kubios HRV version 3.0. This software was also used to apply medium threshold filters for artifact removal if required. Time domain, frequency domain, and Poincare plot were retrieved. The variables analyzed were RR, RMSSD expressed in natural logarithm transformed values (RMSSDln), and pNN50 for time domain variables. Fast Fourier transform was used to calculate the spectral components of the frequency domain. The HF power (0.15–1.0 Hz), LF power (0.04–0.15 Hz), and TP (0.04–1.0 Hz) components were calculated as integrals of the respective power spectral density curve. These variables were also expressed in natural logarithm transformed values (HFln, LFln, and TPln, respectively). Finally, Poincare plot variables such as SD1 and the standard deviation of continuous long-term R-R interval variability (SD2) were calculated. Stress score was analyzed by the formula 1000 × 1/SD2 5 and was presented as a natural logarithm transformed value. Finally, sympathetic/parasympathetic ratio (S/PS) was calculated. 5
STATISTICS
Statistical analyses were conducted with IBM statics SPSS v20.0 software (SPSS Inc, Chicago IL). The Shapiro-Wilk test was used to investigate normality of the data, and analysis of variance was used for paired samples to determine significant differences among measurements. The relationship between variables was analyzed by the Pearson correlation coefficient. Values of P<0.05 were considered significant differences. Results are expressed as mean values±SD.
Results
Body mass markedly decreased after exercise in both HD and HH, from 64.9±10.7 kg to 64.1±10.4 kg (P<0.05) and from 65.0±10.8 kg to 64.2±10.6 kg (P<0.05), respectively.
Table 1 shows the differences between HD and HH for distance covered and percentage decrease in BM. No significant differences between HD and HH (P>0.05) were detected for these variables. Similarly, there were no significant differences in water intake between HD and HH (P>0.05).
Difference between HD and HH for distance covered and body mass change
HD, hot and dry environment; HH, hot and humid environment.
Figure 1 shows an overview of RR during exercise and recovery, and Figure 2 shows how RMSSDln, HFln, Sscoreln, and S/PS ratio changed during exercise and postexercise. It may be observed as well how these variables change during recovery after the test. Note that no differences appeared among any of the variables when comparing HD and HH at the same time point (P>0.05). Furthermore, Figure 3 shows the data for pNN50 for pre- and posttest and recovery. Data for LFln and TPln during the recovery period are presented in Table 2.

Changes in RR during exercise and recovery postexercise. HD, hot and dry environment; HH: hot and humid environment. Values are presented as mean±SD. X-axis is a discontinuous scale. *P<0.05 from PRE for HD. †P<0.05 from PRE for HH.

Changes in RMSSDln, HFln, Sscoreln, and S/PS during exercise and recovery postexercise. Values are presented as mean±SD. X-axis is a discontinuous scale. HD, hot and dry environment; HH, hot and humid environment. *P<0.05 from PRE for HD. †P<0.05 from PRE for HH.

Changes in pNN50 pre- and postexercise and during recovery. pNN50: consecutive intervals that differ by more than 50 ms. HD, hot and dry environment; HH, hot and humid environment. *P<0.05 for HD. †P<0.05 for HH.
Recovery period after exercise for LFln and TPln
LFln, low frequency natural logarithm; TPln, total power natural logarithm; HD, hot and dry environment; PRE, pretrial; POST, post-trial;HH, hot and humid environment.
P<0.05 from PRE.
Table 3 shows the RPE values during exercise under both conditions at all-time points. For both conditions, RPE values in the first 5 min were significantly lower compared to 20, 25, and 30 min of exercise (P<0.05). There were no significant differences between conditions at the same time point. In Table 4, correlations between HRV variables and RPE during exercise are presented separately for HD and HH.
RPE values during exercise under HD and HH
RPE, rate of perceived exertion; HD, hot and dry environment; HH, hot and humid environment.
Correlation indexes between HRV and RPE
HRV, heart rate variability; RPE, rate of perceived exertion; HD, hot and dry environment; HH, hot and humid environment; RR, R-R intervals; SDNN, standard deviation of N-N intervals; HFln, high frequency natural logarithm; LFln, low frequency natural logarithm; TPln, total power natural logarithm.
Discussion
The main finding of this study was that the recovery of HRV to baseline values after exercise under hot conditions does not appear immediately. The effects of high temperatures on HRV are more remarkable in higher RH conditions than in drier environments, where the recovery is faster. The impact of RH has no effect on HRV when exercising at the same intensity under hot conditions.
Parasympathetic-related variables such as RMSSD and HF started to fall dramatically from the beginning of the exercise and then remained constant during the activity from the first 5 min in the frequency domain (HF) and from 10 min in the time domain (RMSSD). In both HD and HH, all participants were asked to run at the same intensity. We did not observe any difference between trials during exercise in HRV and identified no changes in HRV during exercise under different RH conditions in hot conditions. In agreement with our results, previous research states that RMSSD and HF tended to zero during moderate intensity. 17 The recovery periods for both variables were equal; values similar to baseline were identified 4 h after exercise in the HD environment. Conversely, under the HH environment, more than 4 h of recovery was needed. These results may imply that parasympathetic recovery takes place more quickly when the activity is done under HD conditions than under HH conditions. Analyzing more time domain variables, the pNN50 recovery period was even longer after exercise under the HH environment, with significantly lower values in pNN50 after 8 h of recovery compared to baseline values. Other authors did not observe significant recovery of HF after 30 min 31 or of RMSSD and HF after 60 min of recovery. 16 On the contrary, recovery to baseline values was identified in prior investigations after continuous running at anaerobic threshold intensity after 30 min for RMSSD and 60 min for pNN50. 18
Similar results have been observed for RMSSD and HF after interval training in the first 2 h of recovery, whereas different results were determined after moderate intensity endurance training. 2 These data indicate that higher intensity exercise can result in a longer time needed to return to baseline values of HRV; therefore, parasympathetic activity is more susceptible to change because of the high intensity. Furthermore, higher values of RH during exercise in hot conditions seem to be related to slower recovery after exercise. TP and LF power values after continuous exercise do not return to similar baseline levels after 30 min of recovery regardless of exercise intensity. 31 These results are in agreement with those presented here for HH but not for HD. After exercising in HD, TP values did not differ from initial values, with significantly lower values up to 2 h after exercise. Results in LFln values were similar in both HD and HH to those reported in previous research. 31 The TP changes during exercise did not differ among different intensities 32 ; however, our results showed a greater effect on TP during recovery after exercise under HH than for HD conditions.
Previous research 7 ,18,33,34 identified that exercise duration had no effect on time and frequency domain changes. Time and frequency domain variables remained constant during exercise. On the contrary, S/PS did not remain stable during exercise; it increased during the whole trial and did not reach a specific limit. Future research might study how long S/PS can continue increasing in relation to exercise intensity and duration. S/PS returned to baseline values faster in HD than for HH conditions, showing that higher ambient RH produces a greater impact on sympathetic activity during high temperatures. This effect was higher than in time and frequency domains; after 8 h of recovery after HH conditions, S/PS baseline values were not reached, whereas baseline was reached in less than 4 h after exercise in HD conditions. Unfortunately, acute effects on Sscore and S/PS after exercise have not been deeply analyzed in previous studies. No differences among HD and HH conditions were observed during recovery in Sscore. Results were not different from baseline immediately after the exercise; thus, Sscore did not yield any information on differences in the recovery period after exercise in HD and HH conditions. Sscore during exercise markedly increased in the first 10 min, and at minute 15 it remained constant until the end of the exercise, with a decrement in the last 5 min of exercise. Even though significant differences appear in the HRV recovery under humid environments, these differences do not exist at the same time point of recovery. Significance does not always indicate meaningfulness, however, so future research should study this effect of RH under hot conditions with regard to intensity of exercise and controlling for a longer acute effect. This way, more consistent results regarding the effect of RH could be retrieved.
Although we have observed that exercising under HH conditions leads to a faster recovery than HD for ANS, RPE did not differ between conditions during exercise. Additionally, RPE showed an identical increment tendency, with the same differences intratrial in both tests. The relationship between RPE and HRV has been previously described. High correlations were observed between RPE and HF, LF, and TP,7,9 revealing results similar to those reported in HD; this was true only for TP values in HH conditions. Higher correlations with RPE were observed for RR; these results are predictable because RR is proportional to HR (HR) (RR = 60 / HR), and HR presents high correlations with RPE.35,36 In addition, Sscore presented a high correlation value with RPE, as seen in Table 4, which may suggest that Sscore is the best predictor of internal load during exercise if we compare its values to values obtained in the time and frequency domain variables in the correlation with RPE. Sscore presented the highest correlation value with RPE for both HD and HH, but a higher value was obtained for HD. We suppose the fact that Sscore and S/PS ratio presented higher values in the correlation with RPE might be due to the tendency to zero time and frequency domain variables during exercise whereas Sscore and S/PS continued to increase with behavior similar to RPE.
Body mass decreased significantly in both runs, but there were no differences between them. Prior researchers obtained similar results when comparing exercise at 23% RH to 43%, 52%, and 61% RH. 27 Conversely, these authors observed significantly higher values in tympanic temperature when exercising at 71% RH compared to 23%. 27 Thus, higher values of RH would be needed to observe a higher increment in tympanic temperature between trials.
LIMITATIONS OF THE STUDY
In this study, exercise intensity was not controlled. As seen, intensity has a direct effect on HRV; thus, further research should establish a certain intensity. Furthermore, there was a lack of core and skin temperature data, which would have been useful to interpret why higher humidity provokes this response on HRV.
Conclusions
HRV does not behave differently during exercise in hot conditions depending on RH. The recovery of initial values of HRV does not happen immediately after exercise; the time needed to return to baseline levels may be up to 8 h. Higher RH may lead to a slower recovery of ANS modulation after continuous exercise under high-temperature conditions. The recovery to baseline conditions for practitioners took longer after endurance exercise HH conditions than in HD conditions. On the other hand, Sscore appeared to be significant in the assessment of internal load based on the relationship with RPE.
Footnotes
Author Contributions
Study concept and design (OAA); acquisition of the data (OAA, DLP); analysis of data (OAA, FA, PM, JNO); manuscript writing (OAA, DLP); critical revision of the manuscript (PM, FA, JNO); reading and approval of the final manuscript (OAA, DLP, FA, JNO, PM).
Financial/Material Support
None.
Disclosures
None.
