Automatic Detection of Aerobic Threshold through Recurrence Quantification Analysis of Heart Rate Time Series

During exercise with increasing intensity, the human body transforms energy with mechanisms dependent upon actual requirements. Three phases of the body’s energy utilization are recognized, characterized by different metabolic processes, and separated by two threshold points, called aerobic (AerT) and anaerobic threshold (AnT). These thresholds occur at determined values of exercise intensity(workload) and can change among individuals. They are considered indicators of exercise capacities and are useful in the personalization of physical activity plans. They are usually detected by ventilatory or metabolic variables and require expensive equipment and invasive measurements. Recently, particular attention has focused on AerT, which is a parameter especially useful in the overweight and obese population to determine the best amount of exercise intensity for weight loss and increasing physical fitness. The aim of study is to propose a new procedure to automatically identify AerT using the analysis of recurrences (RQA) relying only on Heart rate time series, acquired from a cohort of young athletes during a sub-maximal incremental exercise test (Cardiopulmonary Exercise Test, CPET) on a cycle ergometer. We found that the minima of determinism, an RQA feature calculated from the Recurrence Quantification by Epochs (RQE) approach, identify the time points where generic metabolic transitions occur. Among these transitions, a criterion based on the maximum convexity of the determinism minima allows to detect the first metabolic threshold. The ordinary least products regression analysis shows that values of the oxygen consumption VO(2), heart rate (HR), and Workload correspondent to the AerT estimated by RQA are strongly correlated with the one estimated by CPET (r > 0.64). Mean percentage differences are <2% for both HR and VO(2) and <11% for Workload. The Technical Error for HR at AerT is <8%; intraclass correlation coefficients values are moderate (≥0.66) for all variables at AerT. This system thus represents a useful method to detect AerT relying only on heart rate time series, and once validated for different activities, in future, can be easily implemented in applications acquiring data from portable heart rate monitors.