Modeling the Expenditure and Recovery of Anaerobic Work Capacity in Cycling †

The objective of this research is to model the expenditure and recovery of Anaerobic Work Capacity (AWC) as related to Critical Power (CP) during cycling. CP is a theoretical value at which a human can operate indefinitely and AWC is the energy that can be expended above CP. There are several models to predict AWC-depletion, however, only a few to model AWC recovery. A cycling study was conducted with nine recreationally active subjects. CP and AWC were determined by a 3-min all-out test. The subjects performed interval tests at three recovery intervals (15 s, 30 s, or 60 s) and three recovery powers (0.50CP, 0.75CP, and CP). It was determined that the rate of expenditure exceeds recovery and the amount of AWC recovered is influenced more by recovery power level than recovery duration. Moreover, recovery rate varies by individual and thus, a robust mathematical model for expenditure and recovery of AWC is needed.


Introduction
In this research, experiments are developed to model the expenditure and recovery of Anaerobic Work Capacity (AWC) for power levels above and below Critical Power (CP).CP is defined as a power level at which a human can operate indefinitely and AWC is the amount of anaerobic energy that can be expended at power levels above CP [1].The advent of miniature, real-time sensors has provided opportunities for health monitoring [2] as well as exercise-fatigue modeling.These devices include heartrate monitors, Near-Infrared Spectroscopy (NIRS) devices, GPS-enabled wrist watches, and power meters.The motivation of this research is to allow for prolonged physical exertion and to understand how energy and power are used within human work expenditure.It should be noted that the study performed in this research is focused on cycling.However, the goal is to apply the learnings to other exercise forms like running and sprint-work.Understanding the limitations of the human body can lead to optimized exercise and performance.This can be achieved by optimally reaching fatigue in training and thereby avoiding exhaustion in races.

A Review of Anaerobic Work Capacity and Critical Power
Hill [3] first suggested the hyperbolic relationship between average speed and time while studying athletic records in different sports.Monod and Scherrer [1] then quantified AWC and CP, which were determined through multiple cycling tests to exhaustion across several days.Monod and Sherrer [1] defined Critical Power (CP) as the power a muscle can "keep up for a very long time without fatigue".Anaerobic Work Capacity was defined as Limit Work and was equal to the product of power and time to exhaustion given by, () where, P is power in Watts, tLim is time-to-exhaustion in seconds, and WLim is Limit Work in Joules.This WLim has since been referred to as Anaerobic Work Capacity.To determine CP, a subject is required to exercise till exhaustion at a constant power (P) on at least three different lab visits.These tests yield a hyperbolic plot of power versus time with the asymptote being the CP.The hyperbolic relationship has been transformed into a linear relationship (P vs. 1/tLim) by many researchers, [1,[4][5][6] to name a few.Where, CP is determined by the y-intercept of the plot and AWC is determined by the slope.Morton in [7] developed a ramp exercise to exhaustion protocol to determine CP and AWC, which has been validated in [8].To avoid multiple lab visits to determine CP and AWC, Vanhatalo and colleagues [9] proposed a 3 min all-out test (3mAOT) protocol.The test involves a subject pedaling at all-out intensity for 3 min.CP is given by the average power output of the last 30 s of the test and the area under the power output curve above CP gives the subjects' AWC [9,10].The test has since been validated in [11,12].Morton and Billat [13] suggested that recovery of AWC varies linearly with expenditure of AWC.Ferguson and colleagues [14] investigated charge and discharge of AWC in cycling and suggested that the recovery of AWC is non-linear.Skiba and colleagues [15,16] also investigated recovery of AWC and proposed an exponential model for the "reconstitution" (or recovery as it will be here) of AWC of an athlete given by, where W' is AWC, W'bal is the amount of AWC remaining in Joules, W'exp is the amount of AWC that has been burned in Joules, (t-u) is time spent below CP in seconds, and τW' is the time constant of reconstitution of AWC in seconds.The time constant follows a non-linear trend with DCP, which is the difference between the average recovery power for the entire exercise and CP, given by: While Skiba's work provides a baseline equation to model expenditure and recovery, there are several limitations in the model.First, the recovery rate is dependent upon a DCP calculated for the entire exercise duration.Second, the mathematical model allows for a negative balance of AWC [16].A negative balance of AWC could be due to the standard deviations associated with CP and AWC.However, it would mean that the athlete is drawing more from their AWC stores than they did during the baseline CP testing.Additionally, the integral term of Equation ( 2) results in Joules-second and this does not match with the units of AWC in Joules.Finally, there are two versions of Equation ( 2) [15,16].

Hypothesized Model of Recovery and Expenditure
The model presented in this research is based on the relationship shown in Equation (1) for the expenditure of AWC.Equation ( 1) is modified to reflect both expenditure and recovery of AWC as,

 
() where, ΔWAn is a change in energy or an amount of work expended or recovered from a single duration of exertion or recovery, P is the power level operated at for an expenditure interval, CP is the Critical Power of the subject, and t is the time of the expenditure interval.
To model a specific power level, β(t) is defined as a percent of Critical Power: With this relationship, if β is less than 1, the subject is recovering; if β is equal to 1, the subject is at their CP; and if β is greater than 1, the subject is depleting their energy stores.The expenditure/recharge factor, Φ is represented in a piece-wise form as a function of β, Combining Equations ( 4)-( 6) results in the hypothesized model of expenditure and recovery as, ,

Participants
Nine recreationally active subjects (8 male, 1 female aged 23-44) voluntarily participated in the study.An informed consent form and a health survey were signed and submitted by each participant before testing began.The subjects were instructed to avoid strenuous exercise 24 h before each test and to arrive well-hydrated and well-rested.It was suggested that they have a carbohydrate-rich snack (e.g., a gel-pack or a banana) about an hour before testing and to avoid caffeine for at least 3 h before testing.All tests were performed under the compliance of Clemson University IRB2016-169.

Equipment and Data Collection
The 3-min all out test and 8 interval tests were conducted on a system comprising of a CompuTrainer cycling trainer with Racermate software, a BSX Insight muscle oxygenation (SmO2) sensor, and a Garmin Heart Rate monitor.The CompuTrainer system includes a cadence sensor which reads to the control panel.Data was analyzed using Microsoft Excel and Minitab.

Experimental Design for Determining Phi, Φ
An experimental protocol (see Figure 1) was set up so that a known amount of AWC was burned, followed by a recovery interval, and then a constant power interval where time to exhaustion was measured.In Figure 1, CP4 refers to the work-rate at which a subject will reach exhaustion in 4 min.From Equation (7), Φ is determined as, ) Three values of β (50% CP, 75% CP, and 100% CP) and three values of tn (15 s, 30 s, and 60 s) were selected.The combination of β = 50% CP and t2 = 60 s was eliminated because pilot testing indicated complete recovery of AWC.On the first day, the subjects performed a modified version of Vanhatalo's 3-min all out test.For the following tests, on separate days, the subjects followed the protocol shown in Figure 1.The point of complete depletion of AWC was determined to be the time at which a subject's preferred cadence fell by >10 rpm for >10 s, despite continued encouragement.The subjects then cooled down until their heart rate and SmO2 values returned to those of resting.
The subjects performed two tests each day and the testing days were separated by at least 24 h.Positive encouragement was provided throughout all the tests.

Determining Critical Power and Anaerobic Work Capacity
The validity of the 3mAOT using CompuTrainer has been documented by Clark and colleagues in [17].To ensure that the 3mAOT on CompuTrainer was yielding comparable results, the average and standard deviation from Vanhatalo's 3-min all out test [9] were studied.The subjects were grouped into two categories based on the standard deviation of their CP (see Table ).If a subject were to pace themselves during the 3mAOT, their CP and AWC would be miscalculated, which could lead to issues in determining statistical significance in the behavior of Φ. Results from Group 1 are in agreement with the results seen in [9].Results similar to that of Group 2 can also be seen in 3mAOT performed by inexperienced cyclists as seen in Johnson and colleagues investigation in [12].

Results and Discussions
As discussed earlier, isolating a recovery interval is the key to determining Φ for each specific β and recovery time, t2.Using Equation (10), the amount of Anaerobic Work expended in both intervals was calculated as the sum of the Anaerobic Work expended each second: where Δt = 1 s and Φ =1 (expenditure intervals).The value for β was calculated to be the average power across the recovery interval.On some days, the subjects could not hold their CP4 for the 4 min minimally necessary to calculate Φ.This resulted in a negative Φ as it was assumed that the subject started the test with 100% AWC balance.It is possible that these subjects had not fueled in the same way the day before the testing and were not testing at the same initial anaerobic work balance as they were for the all-out test.These data points were not considered for statistical analysis.
An Anderson-Darling test was performed with α = 0.05 (level of significance) to calculate the normality of the data points.The p-value was calculated to be 0.375 and the null-hypothesis that the data follows a normal distribution could not be rejected.A t-test was then carried out to determine the probability of the mean value of Φ being equal to 1 with α = 0.05.For this test, only valid data points from the first testing session of each day were used, to follow the practices found in literature of giving subjects at least 24 h to recover before performing another test.The p-value was determined to be 0.029 (<0.05), indicating that the subjects did not recover AWC at the same rate as they expended their AWC stores.This aligns with the work done by Ferguson and colleagues [14], Skiba and colleagues [15,16], and Chidnok and colleagues [18].Parallels cannot be drawn to Chidnok's study as the power levels are based on CP in the study presented here.From the data collected, a regression equation for Φ in terms of β and Δt was determined to be: 8.07 14.80 0.0180 t Equation ( 12) does not represent a universal recovery rate for all subjects and cannot be used to determine recovery rates for athletes.An Analysis of Variance (ANOVA) was performed on the regression to determine the probability of the coefficients of β and Δt being equal to zero with α = 0.05 and the p-value was determined to be 0.001.Therefore, it can be said that the regression equation provides a better model than the constant term alone.The ANOVA showed that the coefficient of β was significant (p-value = 0), whereas coefficient of Δt was not significant alone (p-value = 0.524).The recovery factor Φ followed different trends for different subjects.It might be prudent to investigate the recovery of AWC among similar athletes as the participants in this study were recreationally active people pursuing different sports.The results suggest that it might be worthwhile to look at subjects individually rather than as a group.These findings suggest that an optimal control system can be designed with objectives of minimizing race time and complete depletion of AWC by the end of the race.The system would need to be implemented into a humanin-the-loop feedback control system, as there are factors affecting the race day performance that cannot always be accounted for in training.It would be possible to achieve this human-in-the-loop system using a real-time multiple sensor network comprising of a power meter, cadence, GPS, and bike computer.For example, optimal power can be determined based on pre-race planning as well as in-race modifications through feedback control approaches.A cyclist would be prompted to adjust power by accounting for factors such as wind speed, ground conditions, and deviation from the optimal power.This approach can be adapted to optimize distance, efficiency, or time.The findings from this research provide encouraging signs to model the recovery of AWC.The authors are currently conducting tests with a larger participant pool to develop models of individual specific expenditure and recovery and a more generalized model of the population.

Conclusions and Future Work
The objective of the work presented in this paper was to understand the underpinnings of recovery of AWC and to model the same in terms of CP.In current literature, the lack of a robust understanding of recovery of AWC was identified.A model is proposed, and experiments were conducted to understand expenditure and recovery.Through the experiments, it was determined that that sub-CP recovery is not equal to above-CP expenditure.A regression equation for rate of recovery of AWC was determined in terms of recovery duration and recovery power.From the regression, it was determined that the amount of AWC recovered is more dependent upon the recovery power level than the time spent in recovery.It was observed that the rate of recovery varied significantly between subjects and hence it is important to look at athlete-specific recovery rates.The findings suggest a human in the loop optimal control system can be designed to optimize performance.Furthermore, AWC balance could be linked to overtraining and injury not only in cycling and running, but other sports involving multiple sprints like football and soccer.Another future research direction would be to determine the effect of power (as %V ̇O2max ) on the value of Φ.This would help in gaining more insights into the relationships between power output and the physiological mechanisms occurring during fatigue.

Figure 1 .
Figure 1.Visual Representation of Experimental Protocol.

Table 1 .
Average and standard deviation of final 30 s of 3-min test.