Real-Time Estimation of Aerobic Threshold and Exercise Intensity Distribution Using Fractal Correlation Properties of Heart Rate Variability: A Single-Case Field Application in a Former Olympic Triathlete

Heart rate dynamics and asymmetry during sympathetic activity stimulation and post-stimulation recovery in ski mountaineers-a pilot exploratory study

There is a lack of studies on non-linear heart rate (HR) variability in athletes. We aimed to assess the usefulness of short-term HR dynamics and asymmetry parameters to evaluate the neural modulation of cardiac activity based on non-stationary RR interval series by studying their changes during sympathetic nervous system activity stimulation (isometric handgrip test) and post-stimulation recovery in professional ski mountaineers. The correlation between the changes in the parameters and the respiratory rate (RespRate) and also the duration of the career was analyzed. Short-term (5 min) and ultra-short-term (1 min) rates of patterns with no variations (0V), number of acceleration runs of length 1 (AR1), and short-term Porta's Index were greater, whereas Guzik's Index (GI) was smaller during sympathetic stimulation compared to rest. GI increased and the number of AR1 decreased during recovery. Greater increases in GI and RMSSD were associated with greater decreases in RespRate during recovery. Greater increases in RespRate from rest to short-term sympathetic stimulation were associated with greater increases in 0V (Max-min method) and AR1 but also with greater decreases in decelerations of short-term variance and accelerations and decelerations of long-term variance. Greater increases in 0V (Max-min method) and number of AR1 during sympathetic stimulation were associated with a shorter career duration. Greater decreases in these parameters during recovery were associated with a longer career duration. Changes in measures of HR dynamics and asymmetry, calculated based on short-term non-stationary RRi time series induced by sympathetic stimulation and post-stimulation recovery, reflected sympathovagal shift and were associated with condition-related alterations in RespRate and career duration in athletes who practice ski mountaineering.

Correlation properties of heart rate variability to assess the first ventilatory threshold and fatigue in runners

PURPOSE

The short-term scaling exponent alpha1 of detrended fluctuation analysis (DFA-a1) of heart rate variability (HRV) has shown potential to delineate the first ventilatory threshold (VT1). The aims of this study were to investigate the accuracy of this method for VT1 determination in runners using a consumer grade chest belt and to explore the effects of acute fatigue.

METHODS

We compared oxygen uptake (V̇O2) and heart rate (HR) at gas exchange VT1 to V̇O2 and HR at a DFA-a1 value of 0.75. Gas exchange and HRV data were obtained from 14 individuals during a treadmill run involving two incremental ramps. Agreement was assessed using Bland-Altman analysis and linear regression.

RESULTS

Bland-Altman analysis between gas exchange and HRV V̇O2 and HR at VT1 during the first ramp showed a mean (95% limits of agreement) bias of -0.5 (-6.8 to 5.8) ml∙kg-1∙min-1, and -0.9 (-12.2 to 10.5) beats∙min-1, with R2 of 0.83 and 0.56, respectively. During the second ramp, the differences were -7.3 (-18.1 to 3.5) ml∙kg-1∙min-1 and -12.3 (-30.4 to 5.9) beats∙min-1, with R2 of 0.62 and 0.43, respectively.

CONCLUSION

A chest-belt derived DFA-a1 of 0.75 is closely related to gas exchange VT1, with the variability in accuracy at an individual level being similar to gas exchange methods. This suggests this to be a useful method for exercise intensity demarcation. The altered relationship during the second ramp indicates that DFA-a1 is only able to accurately demarcate exercise intensity thresholds in a non-fatigued state, but also opens opportunities for fatigue-based training prescription.

Heart Rate Variability-Derived Thresholds for Exercise Intensity Prescription in Endurance Sports: A Systematic Review of Interrelations and Agreement with Different Ventilatory and Blood Lactate Thresholds

BACKGROUND

Exercise intensities are prescribed using specific intensity zones (moderate, heavy, and severe) determined by a 'lower' and a 'higher' threshold. Typically, ventilatory (VT) or blood lactate thresholds (LT), and critical power/speed concepts (CP/CS) are used. Various heart rate variability-derived thresholds (HRVTs) using different HRV indices may constitute applicable alternatives, but a systematic review of the proximity of HRVTs to established threshold concepts is lacking.

OBJECTIVE

This systematic review aims to provide an overview of studies that determined HRVTs during endurance exercise in healthy adults in comparison with a reference VT and/or LT concept.

METHODS

A systematic literature search for studies determining HRVTs in healthy individuals during endurance exercise and comparing them with VTs or LTs was conducted in Scopus, PubMed and Web of Science (until January 2022). Studies claiming to describe similar physiological boundaries to delineate moderate from heavy (HRVTlow vs. VTlow and/or LTlow), and heavy from severe intensity zone (HRVThigh vs. VThigh and/or LThigh) were grouped and their results synthesized.

RESULTS

Twenty-seven included studies (461 participants) showed a mean difference in relative HR between HRVTlow and VTlow of - 0.6%bpm in weighted means and 0.02%bpm between HRVTlow and LTlow. Bias between HR at HRVTlow and VTlow was 1 bpm (limits of agreement (LoA): - 10.9 to 12.8 bpm) and 2.7 bpm (LoA: - 20.4 to 25.8 bpm) between HRVTlow and LTlow. Mean difference in HR between HRVThigh and VThigh was 0.3%bpm in weighted means and 2.9%bpm between HRVThigh and LThigh while bias between HR at HRVThigh and VThigh was - 4 bpm (LoA: - 17.9 to 9.9 bpm) and 2.5 bpm (LoA: - 12.1 to 17.1 bpm) between HRVThigh and LThigh.

CONCLUSION

HRVTlow seems to be a promising approach for the determination of a 'lower' threshold comparable to VTlow and potentially for HRVThigh compared to VThigh, although the latter needs further empirical evaluation. LoA for both intensity zone boundaries indicates bias of HRVTs on an individual level. Taken together, HRVTs can be a promising alternative for prescribing exercise intensity in healthy, male athletes undertaking endurance activities but due to the heterogeneity of study design, threshold concepts, standardization, and lack of female participants, further research is necessary to draw more robust and nuanced conclusions.

Validation of a non-linear index of heart rate variability to determine aerobic and anaerobic thresholds during incremental cycling exercise in women

Studies highlight the usage of non-linear time series analysis of heart rate variability (HRV) using the short-term scaling exponent alpha1 of Detrended Fluctuation Analysis (DFA-alpha1) during exercise to determine aerobic and anaerobic thresholds. The present study aims to further verify this approach in women. Gas exchange and HRV data were collected from 26 female participants with different activity levels. Oxygen uptake (VO₂) and heart rate (HR) at first (VT1) and second ventilatory thresholds (VT2) were compared with DFA-alpha1-based thresholds 0.75 (HRVT1) and 0.50 (HRVT2). Results: VO₂ at VT1 and VT2 were 25.2 ml/kg/min (± 2.8) and 31.5 ml/kg/min (± 3.6) compared with 26.5 ml/kg/min (± 4.0) and 31.9 ml/kg/min (± 4.5) for HRVT1 and HRVT2, respectively (ICC_3,1 = 0.77, 0.84; r = 0.81, 0.86, p < 0.001). The mean HR at VT1 was 147 bpm (± 15.6) and 167 bpm (± 12.7) for VT2, compared with 152 bpm (± 15.5) and 166 bpm (± 13.2) for HRVT1 and HRVT2, respectively (ICC_3,1 = 0.87, 0.90; r = 0.87, 0.90, p < 0.001). Bland-Altman analysis for VT1 vs. HRVT1 showed a mean difference of - 1.3 ml/kg/min (± 2.4; LoA: 3.3, - 6.0 ml/kg/min) for VO₂ and of - 4.7 bpm (± 7.8; LoA: 10.6, - 20.0 bpm) for HR. VT2 vs. HRVT2 showed a mean difference of - 0.4 ml/kg/min (± 2.3; LoA: 4.1, - 4.9 ml/kg/min) for VO₂ and 0.5 bpm (± 5.7; LoA: 11.8, - 10.8 bpm) for HR. DFA-alpha1-based thresholds showed good agreement with traditionally used thresholds and could be used as an alternative approach for marking organismic transition zones for intensity distribution in women.

Fractal correlation properties of HRV as a noninvasive biomarker to assess the physiological status of triathletes during simulated warm-up sessions at low exercise intensity: a pilot study

BACKGROUND

The non-linear index alpha 1 of Detrended Fluctuation Analysis (DFA a1) of heart rate variability, has been shown to be a marker of fatigue during endurance exercise. This report aims to explore its ability to assess the physiological status as a surrogate metric for "readiness to train" while performing simulated warm-up sessions the day after two different exercise sessions.

METHODS

11 triathletes were recruited to determine the first ventilatory threshold (VT1) during a baseline assessment and to perform 10-min of cycling at 90% of VT1 (simulating a warm-up bout) before (PRE) and within 36 h after (POST) light and heavy running exercise. RR intervals were recorded for DFA a1 analysis along with neuromuscular testing to verify the effects of the performed exercise sessions. In addition to common statistical methods, magnitude-based inferences (MBI) were applied to assess the changes in true score and thus also the practical relevance of the magnitude.

RESULTS

Rating of perceived exertion for the heavy exercise session showed a significant higher rating as opposed to the light exercise session (p < 0.001, d = 0.89). In regard of MBIs, PRE versus POST comparisons revealed a significant reduced DFA a1 with large effect size after the heavy exercise session (p = 0.001, d = - 1.44) and a 99% chance that this negative change was clinically relevant.

CONCLUSIONS

Despite inter-individual differences, DFA a1 offers potential to assess physiological status and guide athletes in their training as an easy-to-apply monitoring procedure during a standardized warm-up. A regular assessment including individual data history and statistical references for identification of response is recommended. Further data are necessary to confirm the results in a larger and more homogeneous population.

Fractal Correlation Properties of Heart Rate Variability as a Biomarker for Intensity Distribution and Training Prescription in Endurance Exercise: An Update

While established methods for determining physiologic exercise thresholds and intensity distribution such as gas exchange or lactate testing are appropriate for the laboratory setting, they are not easily obtainable for most participants. Data over the past two years has indicated that the short-term scaling exponent alpha1 of Detrended Fluctuation Analysis (DFA a1), a heart rate variability (HRV) index representing the degree of fractal correlation properties of the cardiac beat sequence, shows promise as an alternative for exercise load assessment. Unlike conventional HRV indexes, it possesses a dynamic range throughout all intensity zones and does not require prior calibration with an incremental exercise test. A DFA a1 value of 0.75, reflecting values midway between well correlated fractal patterns and uncorrelated behavior, has been shown to be associated with the aerobic threshold in elite, recreational and cardiac disease populations and termed the heart rate variability threshold (HRVT). Further loss of fractal correlation properties indicative of random beat patterns, signifying an autonomic state of unsustainability (DFA a1 of 0.5), may be associated with that of the anaerobic threshold. There is minimal bias in DFA a1 induced by common artifact correction methods at levels below 3% and negligible change in HRVT even at levels of 6%. DFA a1 has also shown value for exercise load management in situations where standard intensity targets can be skewed such as eccentric cycling. Currently, several web sites and smartphone apps have been developed to track DFA a1 in retrospect or in real-time, making field assessment of physiologic exercise thresholds and internal load assessment practical. Although of value when viewed in isolation, DFA a1 tracking in combination with non-autonomic markers such as power/pace, open intriguing possibilities regarding athlete durability, identification of endurance exercise fatigue and optimization of daily training guidance.

Validity of Detrended Fluctuation Analysis of Heart Rate Variability to determine intensity thresholds in professional cyclists

BACKGROUND

The evaluation of performance in endurance athletes and the subsequent individualization of training is based on the determination of individual physiological thresholds during incremental tests. Gas exchange or blood lactate analysis are usually implemented for this purpose, but these methodologies are expensive and invasive. The short-term scaling exponent alpha 1 of detrended Fluctuation Analysis (DFA-α1) of the Heart Rate Variability (HRV) has been proposed as a non-invasive methodology to detect intensity thresholds.

PURPOSE

The aim of this study is to analyse the validity of DFA-α1 HRV analysis to determine the individual training thresholds in elite cyclists and to compare them against the lactate thresholds.

METHODOLOGY

38 male elite cyclists performed a graded exercise test to determine their individual thresholds. HRV and blood lactate were monitored during the test. The first (LT1 and DFA-α1-0.75, for lactate and HRV, respectively) and second (LT2 and DFA-α1-0.5, for lactate and HRV, respectively) training intensity thresholds were calculated. Then, these points were matched to their respective power output (PO) and heart rate (HR).

RESULTS

There were no significant differences (p > 0.05) between the DFA-α1-0.75 and LT1 with significant positive correlations in PO (r = 0.85) and HR (r = 0.66). The DFA-α1-0.5 was different against LT2 in PO (p = 0.04) and HR (p = 0.02), but it showed significant positive correlation in PO (r = 0.93) and HR (r = 0.71).

CONCLUSIONS

The DFA1-a-0.75 can be used to estimate LT1 non-invasively in elite cyclists. Further research should explore the validity of DFA-α1-0.5.

An Index of Non-Linear HRV as a Proxy of the Aerobic Threshold Based on Blood Lactate Concentration in Elite Triathletes

A non-linear index of heart rate (HR) variability (HRV) known as alpha1 of Detrended Fluctuation Analysis (DFA a1) has been shown to change with increasing exercise intensity, crossing a value of 0.75 at the aerobic threshold (AT) in recreational runners defining a HRV threshold (HRVT). Since large volumes of low-intensity training below the AT is recommended for many elite endurance athletes, confirmation of this relationship in this specific group would be advantageous for the purposes of training intensity distribution monitoring. Nine elite triathletes (7 male, 2 female) attended a training camp for diagnostic purposes. Lactate testing was performed with an incremental cycling ramp test to exhaustion for the determination of the first lactate threshold based on the log–log calculation method (LT1). Concurrent measurements of cardiac beta-to-beat intervals were performed to determine the HRVT. Mean LT1 HR of all 9 participants was 155.8 bpm (±7.0) vs. HRVT HR of 153.7 bpm (±10.1) (p = 0.52). Mean LT1 cycling power was 252.3 W (±48.1) vs. HRVT power of 247.0 W (±53.6) (p = 0.17). Bland–Altman analysis showed mean differences of −1.7 bpm and −5.3 W with limits of agreement (LOA) 13.3 to −16.7 bpm and 15.1 to −25.6 W for HR and cycling power, respectively. The DFA a1-based HRVT closely agreed with the LT1 in a group of elite triathletes. Since large volumes of low-intensity exercise are recommended for successful endurance performance, the fractal correlation properties of HRV show promise as a low-cost, non-invasive option to that of lactate testing for identification of AT-related training boundaries.

Analysis of Fractal Correlation Properties of Heart Rate Variability during an Initial Session of Eccentric Cycling

Eccentric cycling (ECC) has attracted attention as a method to improve muscle strength and aerobic fitness in populations unable to tolerate conventional methods. However, agreement on exercise prescription targets have been problematic. The current report is an initial exploration of a potentially useful tool, a nonlinear heart rate (HR) variability (HRV) index based on the short-term scaling exponent alpha1 of detrended fluctuation analysis (DFA a1), which has been previously shown to correspond to exercise intensity. Eleven male volunteers performed 45 min of concentric (CON) cycling and ECC separated by 1 month. Work rates were matched for HR (~50% of the maximal HR) during the first 5 min and remained stable thereafter. HRV, HR, oxygen consumption (VO2), and cycling power were monitored and evaluated at elapsed times of 10 (T10) and 45 (T45) minutes duration. HR significantly increased between ECC T10 and ECC T45 (p = 0.003, d = 1.485), while DFA a1 significantly decreased (p = 0.004, d = 1.087). During CON, HR significantly increased (p < 0.001 d = 1.570) without significant DFA a1 change (p = 0.48, d = 0.22). Significantly higher HR was observed at T45 in ECC than in CON (p = 0.047, d = 1.059). A session of unaccustomed ECC lead to decreased values of DFA a1 at T45 in comparison to that seen with CON at similar VO2. ECC lead to altered autonomic nervous system balance as reflected by the loss of correlation properties compared to CON.

Aerobic Threshold Identification in a Cardiac Disease Population Based on Correlation Properties of Heart Rate Variability

An index of heart rate (HR) variability correlation properties, the short-term scaling exponent alpha1 of detrended fluctuation analysis (DFA a1) has shown potential to delineate the first ventilatory threshold (VT1). This study aims to extend this concept to a group of participants with cardiac disease. Sixteen volunteers with stable coronary disease or heart failure performed an incremental cycling ramp to exhaustion PRE and POST a 3-week training intervention. Oxygen uptake (VO₂) and HR at VT1 were obtained from a metabolic cart. An ECG was processed for DFA a1 and HR. The HR variability threshold (HRVT) was defined as the VO₂, HR or power where DFA a1 reached a value of 0.75. Mean VT1 was reached at 16.82 ± 5.72 mL/kg/min, HR of 91.3 ± 11.9 bpm and power of 67.8 ± 17.9 watts compared to HRVT at 18.02 ± 7.74 mL/kg/min, HR of 94.7 ± 14.2 bpm and power of 73.2 ± 25.0 watts. Linear relationships were seen between modalities, with Pearson’s r of 0.95 (VO₂), 0.86 (HR) and 0.87 (power). Bland–Altman assessment showed mean differences of 1.20 mL/kg/min, 3.4 bpm and 5.4 watts. Mean peak VO₂ and VT1 did not change after training intervention. However, the correlation between PRE to POST change in VO₂ at VT1 with the change in VO₂ at HRVT was significant (r = 0.84, p < 0.001). Reaching a DFA a1 of 0.75 was associated with the VT1 in a population with cardiac disease. VT1 change after training intervention followed that of the HRVT, confirming the relationship between these parameters.

Fractal correlation properties of heart rate variability as a biomarker of endurance exercise fatigue in ultramarathon runners

Although heart rate variability (HRV) indexes have been helpful for monitoring the fatigued state while resting, little data indicate that there is comparable potential during exercise. Since an index of HRV based on fractal correlation properties, alpha 1 of detrended fluctuation analysis (DFA a1) displays overall organismic demands, alteration during exertion may provide insight into physiologic changes accompanying fatigue. Two weeks after collecting baseline demographic and gas exchange data, 11 experienced ultramarathon runners were divided into two groups. Seven runners performed a simulated ultramarathon for 6 h (Fatigue group, FG) and four runners performed daily activity over a similar period (Control group, CG). Before (Pre) and after (Post) the ultramarathon or daily activity, DFA a1, heart rate (HR), running economy (RE) and countermovement jumps (CMJ) were measured while running on a treadmill at 3 m/s. In Pre versus Post comparisons, data showed a decline with large effect size in DFA a1 post intervention only for FG (Pre: 0.71, Post: 0.32; d = 1.34), with minor differences and small effect sizes in HR (d = 0.02) and RE (d = 0.21). CG showed only minor differences with small effect sizes in DFA a1 (d = 0.19), HR (d = 0.15), and RE (d = 0.31). CMJ vertical peak force showed fatigue-induced decreases with large effect size in FG (d = 0.82) compared to CG (d = 0.02). At the completion of an ultramarathon, DFA a1 decreased with large effect size while running at low intensity compared to pre-race values. DFA a1 may offer an opportunity for real-time tracking of physiologic status in terms of monitoring for fatigue and possibly as an early warning signal of systemic perturbation.