 Research
 Open Access
 Published:
On the kinetics of anaerobic power
Theoretical Biology and Medical Modelling volume 9, Article number: 29 (2012)
Abstract
Background
This study investigated two different mathematical models for the kinetics of anaerobic power. Model 1 assumes that the work power is linear with the work rate, while Model 2 assumes a linear relationship between the alactic anaerobic power and the rate of change of the aerobic power. In order to test these models, a cross country skier ran with poles on a treadmill at different exercise intensities. The aerobic power, based on the measured oxygen uptake, was used as input to the models, whereas the simulated blood lactate concentration was compared with experimental results. Thereafter, the metabolic rate from phosphocreatine break down was calculated theoretically. Finally, the models were used to compare phosphocreatine break down during continuous and interval exercises.
Results
Good similarity was found between experimental and simulated blood lactate concentration during steady state exercise intensities. The measured blood lactate concentrations were lower than simulated for intensities above the lactate threshold, but higher than simulated during recovery after high intensity exercise when the simulated lactate concentration was averaged over the whole lactate space. This fit was improved when the simulated lactate concentration was separated into two compartments; muscles + internal organs and blood. Model 2 gave a better behavior of alactic energy than Model 1 when compared against invasive measurements presented in the literature. During continuous exercise, Model 2 showed that the alactic energy storage decreased with time, whereas Model 1 showed a minimum value when steady state aerobic conditions were achieved. During interval exercise the two models showed similar patterns of alactic energy.
Conclusions
The current study provides useful insight on the kinetics of anaerobic power. Overall, our data indicate that blood lactate levels can be accurately modeled during steady state, and suggests a linear relationship between the alactic anaerobic power and the rate of change of the aerobic power.
Background
During human exercise, it is well known that intracellular adenosine triphosphate (ATP) can be produced aerobically in the mitochondria by oxidative phosphorylation, anaerobically due to glycolysis or glycogenolysis, or by the breakdown of phosphocreatine (PCr) into Creatine (Cr) in the Creatine Kinease (CK) reaction. With the aerobic energy as input, the current paper develops mathematical models that simulate the kinetics of anaerobic ATP production and power.
The rate of oxygen (O_{2}) consumption can be set to the sum of a constant rate (resting rate of O_{2} consumption), a rate due to unloaded body movements and at a rate that is used to perform work [1]. For moderate constant work rates, the O_{2} consumption increases to a steady state condition. However, at the onset of exercise or due to a change in work rate, there is a rate of change of VO_{2} (or aerobic power) that is named as the VO_{2} kinetics. Pulmonary rate of O_{2} consumption (VO_{2p}) has been used as a proxy for VO_{2}. For moderate intensity exercise at constant work rate below the lactate threshold three distinct phases have been observed for VO_{2p}. Phase I is the cardio dynamic phase, which represents the circulatory transit delay from muscles to lungs. Phase II is the mono exponential increase of VO_{2p}. This phase reflects the adjustment of VO_{2} due to the use of active skeletal muscles. Phase III is the steady state phase of VO_{2p} and VO_{2} during moderate exercise intensities [2, 3]. For work rates associated with sustained acidosis, the monoexponential component is slowed compared with lower intensities below the lactate threshold. In addition, a delayed slow component is superimposed. The slow component begins around 100200 s into the exercise and can result in either a delayed submaximal steady state or a steady state equal to the maximal oxygen uptake (VO_{2} max) [4]. However, the mechanism of this slow component has not been resolved [5]. In this paper, we model phase II and III because only these two phases are considered relevant for the anaerobic alactic power model.
When exercise intensity increases and the rate of ATP production by oxidative sources becomes insufficient, anaerobic ATP production is required. When ATP is produced by glycolysis or glycogenolysis the endpoint is pyruvate, which can be reduced to lactate or oxidized to CO_{2} or H_{2}O. The blood lactate concentration in the lactate pool is the result of the appearance of lactate from working muscles and various tissues and the disappearance of lactate in the skeletal muscles, the heart, the liver and the kidney cortex [6–9]. During steady state, lactate production (influx) is equal to lactate removal (outflux). As a result, the lactate concentration in the lactate pool stays constant, and the rate of oxygen consumption is the measure of the whole body energy expenditure regardless of the magnitude of lactate production and removal or the absolute blood lactate concentration [10]. The concept of a maximal lactate steady state can be defined as the highest level of intensity where a steady state condition of lactate can be obtained, which is also referred to as the lactate threshold. At exercise intensities above the lactate threshold the rise in the lactate concentration could be attributed to an increase in the rate of lactate appearance or the result of a decrease in the rate of lactate disappearance [10].
The maximum anaerobic energy that can be utilized is proportional to the sum of Cr and lactate that can accumulate in the body. PCr is an energy buffer that supports the transient failure of other metabolic pathways to support ATP production. The equilibrium constant of the CK reaction is around 20 and the slightest drop in ATP allows the reaction to proceed to ATP production [11]. Thus, the ATP concentration stays nearly constant until almost all the PCr is utilized. Rossiter et al. [12] found that the PCr levels follow an exponential time course after changes in work rate before approaching a steady state condition at moderate exercise intensities. In such cases, a strong similarity has been reported for the time constants of the VO_{2} kinetics and the PCr consumption [12]. Margaria [13] was the first to propose a hydraulic model for examining the whole body energy process during exercise. Despite this breakthrough, Margaria’s model was not quantified. Morton [14] presented a generalization of this model that was solved mathematically and compared with experimental data. Here Morton [14] modeled aerobic power and alactic anaerobic power with an exponential time devolvement without any anaerobic lactic power during steady state exercises below the lactate threshold. However, for exercise intensities that are above the lactate threshold, the anaerobic glycolytic energy supply is significant. The association between PCr and VO_{2} rate constants for exercises at such intensities has not yet been systematically reported. Furthermore, during recovery after high intensity exercise, the level of PCr must be restored, the pH must be reestablished and ADP removed. While the PCr recovery is mainly due to oxidative ATP synthesis, the PCr stores may be rebuilt by anaerobic glycolysis [15–17]. Altogether, these insights provide a theoretical background for the models developed in this paper.
The O_{2}deficit formula of Medbø et al. [18] is an alternative model for the anaerobic power that accounts for lactic and alactic anaerobic power in which the chemical coupling efficiencies are assumed to be similar. Medbø et al. [18] suggested that the O_{2}deficit can be calculated by assuming that the metabolic power at intensities above VO_{2}max can be estimated by extrapolating the steady state linear relationship between work rate and VO_{2} at submaximal intensities. The validity of the O_{2}deficit method has been widely debated [e.g., [19–21]. However, a rationale for the O_{2}deficit model is the assumption that the chemical coupling efficiencies of the three sources of ATP synthesis are similar.
The current study investigated two different mathematical models for the kinetics of anaerobic power during whole body exercise at different intensities. In order to test the models during exercise, oxygen uptake and blood lactate concentration were measured while a cross country skier ran with poles on a treadmill. Aerobic power, based on oxygen uptake measurements, was used as input. The lactic anaerobic power was calculated with a model presented by Moxnes and Hausken [10] using lactate concentration averaged over the whole lactate space, and with a model of Moxnes and Sandbakk [22] where the lactate concentration was separated into two compartments: muscles + internal organs and blood. The current study has input from these two previous studies, and focuses on alactic power as a novel contribution. Therefore, the power due to PCr break down was calculated theoretically and compared against the results of Jeneson et al. [23]. Finally, the models were used to compare PCr break down during continuous and interval exercise.
Methods
Overall design
The current study simulated the kinetics of anaerobic powers during whole body exercise at different exercise intensities through the use of Mathematica 8 (Wolfram Research Inc., Champaign, IL, USA). In a first model (Model 1), we suggested that for all exercise intensities, the work power would be linear to the work rate. As an alternative model (Model 2) where the lactic power and energy were the same as for Model 1, we hypothesized that the alactic anaerobic power would be proportional to the rate of change of the aerobic power [24, 25]. For both models we used a first order differential equation for the time development of aerobic power as a function of work rate. As most parameters in the models needed to be fitted to each individual athlete, we fitted all parameters to a male Norwegian national level crosscountry skier, with a body mass of m = 77.5 kg, body height of 181 cm and 600 h of endurance training per year. Thereafter, the simulations were compared with experimental data where this skier was running with poles on a treadmill (see details below).
The aerobic and anaerobic powers
ATP is produced by three different sources:
In these equations, “model” means model assumption. Q_{ r } is the resting metabolic power, set to 80 J/s based on oxygen uptake measurements of this skier during rest in our laboratory. {Q}_{G} is ATP production by glycolysis/glycogenolysis and CK by phosphocreatine break down to creatine. {Q}_{a}^{w} is aerobic power due to internal and external work. The aerobic power is {Q}_{a}={Q}_{r}+{Q}_{a}^{w}.
For each ATP used (produced and consumed) heat and work are outputs. Heat when producing ATP differs between aerobic and anaerobic sources. However, heat when consuming ATP is the same for the aerobic and anaerobic sources. We defined work power and heat power due to oxidative phosphorylation (a), glycolysis/glycogenolysis (G) and phosphocreatine break down to creatine (CK) as:
where “def” means definition. η_{ a }, η_{ G } and η_{ CK } are the chemical coupling efficiencies when producing ATP aerobically, anaerobically by glycolysis/glycogenolysis (G), and anaerobically by CK, respectively. The two chemical coupling efficiencies η_{ G } and η_{ a } are similar, whereas η_{ CK } is larger [26]. η is the mechanical efficiency, which means the work per unit use of energy of 1 ATP. P_{ T } includes both internal and external work.
At exercise intensities that exceed the maximal aerobic power Q_{ max }, the virtual steady state is the steady state that would be attained if it was possible to carry out the exercise under purely aerobic conditions [27]. This virtual steady state is never reached as the aerobic power ends when Q_{ max } is achieved. Figure 1 shows the steady state aerobic power {\overline{Q}}_{a} for an elite skier running with poles at the angle of inclination of 0.105 radians on a large treadmill. A linear relationship applies for {\overline{Q}}_{a} versus work rate up to approximately 95% of Q_{ max }. When repeating the experiment with the same skier, the difference did not exceed 25 W for this skier which is considered an upper bound error margin.
We defined a steady state virtual power {\overline{Q}}_{\mathit{vir}} which is a straight line that crosses Q_{ max }. Thus we set:
where c is a parameter. P is the work rate. Q_{ un } is the metabolic power due to unloaded running with poles. In general, Q_{ un } depends on the angle of inclination and the cycle frequency of the activity. For a given incline the frequency depends of the work rate. Thus Q_{ un } becomes a function of the work rate even for a fixed incline. {\overline{P}}_{\mathit{vir}} includes both internal and external work power. {\overline{Q}}_{\mathit{vir}} is input to the further modeling. We set that {\overline{Q}}_{\mathit{vir}} is a linear interpolation function through our data points. For values above or below the maximum and minimum measured {\overline{Q}}_{\mathit{vir}} we used linear extrapolation. During steady state, {Q}_{a}^{w}={\overline{Q}}_{a}{Q}_{r}=Min\left({Q}_{\mathit{max}}{Q}_{r},{\overline{Q}}_{\mathit{vir}}{Q}_{r}\right).
The aerobic power is delayed with a time lag during steady state work rate. Q_{ a }(t) is the aerobic power and Q_{ vir } is the virtual aerobic power, {Q}_{a}\left(t\right)\stackrel{model}{=}Min\left({Q}_{\mathit{max}},{Q}_{\mathit{vir}}\left(t\right)\right). To account for the delay mathematically, we used a first order differential equation of the virtual power, to read:
The “dot” above a variable means time derivative. τ_{ a } is the socalled efolding time, which is the time interval when an exponentially growing quantity increases by a factor of e. Thus τ_{ a } is a time parameter that scales from the initiation of the activity until the aerobic power reaches a steady state asymptotically when assuming a constant work rate. In practice, reaching steady state means having less than 1% change in aerobic power per second. The cardio dynamic phase is not necessary to account for since only the true aerobic power is used as input for our anaerobic models. Typical values for τ_{ a } are shown to be 1036 s for moderate intensity exercise [22, 27–33]. For example, di Prampero [27] suggested 1024 s, Ceretelli et al. [32] found that τ_{ a } increases linearly with the concentration of lactate up to 36 s, and Bizoni et al. [33] found τ_{ a } = 23 s for all work rates. We set that τ_{ a } = 30 s as a compromise. Equation (4b) ensures that the aerobic power is less than the maximal aerobic power.
If the virtual metabolic power is equal to or below Q_{ max } we have that {Q}_{a}=Min\left({Q}_{\mathit{max}},{Q}_{\mathit{vir}}\right)={Q}_{\mathit{vir}}. This gives from equation (4a and b) that:
For steady work rate equation (5) gives the analytical solution:
If the work rate is sufficiently low, {\overline{Q}}_{a}={\overline{Q}}_{\mathit{vir}}. Notable is that solution (6) will only apply for a restricted time period unless {\overline{Q}}_{a}={\overline{Q}}_{\mathit{vir}}.
The anaerobic power due to anaerobic glycolysis or glycogenolysis can be calculated from the total lactate concentration by using the equation from di Prampero and Ferretti [27], to read:
C is the lactate concentration defined as the amount of lactate per unit volume of lactate space (including muscles and blood). m is the mass of the body. λ converts the lactate per unit body mass to oxygen consumption. Note that during dynamic situations the blood lactate concentration differs significantly from the muscle lactate concentration. In this article, we measured the blood lactate concentration, which is different from the muscle lactate concentration. However, earlier research has demonstrated that the lactate concentration can be calculated as summarized in Appendix A [10].
The lactic energy used from time t_{0} to time t is:
Experimentally, E_{ G }(t_{0}, t) can be found by measuring changes in the lactate concentration before and after exercise.
From (2) we have that:
The metabolic power is then given by:
Steady state is achieved at exercise intensities below the lactate threshold. We defined steady state with the lactate concentration regarded as steady (i.e., no lactic power) and the aerobic power was steady, to read \stackrel{.}{C}={Q}_{G}={\stackrel{.}{Q}}_{a}^{w}=0. Thus {Q}_{a}^{w} is a constant and Q{Q}_{r}={Q}_{a}^{w} during steady state. Hence, it follows that Q_{ CK } = 0. Thus the alactic power is zero. To achieve Q_{ CK } = 0 during steady state in (9) we must have {P}_{T}=\eta {\eta}_{a}{Q}_{a}^{w}. However, during steady state (below LT) {Q}_{a}^{w}={Q}_{\mathit{un}}+cP. Thus, below the lactate threshold {P}_{T}=\left({Q}_{\mathit{un}}+cP\right)\eta {\eta}_{a}={\overline{P}}_{\mathit{vir}} during steady state.
To close the equations in (9) and (10), which produce an indeterminate solution by themselves unless we have steady state, we needed a model for the work power P_{ T } in general. As Model 1 we forecasted that the work power P_{ T } was similar to the steady state virtual aerobic metabolic power. Thus, Model 1:
Inserting (11) into (9)(10) gives the alactic power as:
and the metabolic rate as:
To calculate {Q}_{\mathit{CK}}^{1}\left(t\right) for a work rate P(t) in (12a) we inserted {\overline{Q}}_{\mathit{vir}} from (3), Q_{ a }(t) from (4) and Q_{ G }(t) from (7). It should be noted that the mechanical efficiency η is absent in (12a and b). We conceived that the work power P_{ T } was proportional to the rate of consumption of ATP. Thus, an equivalent model to (11) would be that the rate of ATP consumption is linear with the work rate.
The alactic energy is for Model 1:
The anerobic power and the anaerobic energy used from time {t}_{0} to time t for Model 1 are:
and:
where E_{ DF }(t_{0}, t) is in the literature named the oxygen deficit of the exercise. Medbø et al. [18] and Losnegard et al. [29] calculated the anaerobic energy during exercise as E_{ DF }(t_{0}, t) (assuming that the rate of change of O_{2} consumption is proportional to the aerobic power). Equation (14b) shows that {E}_{\mathit{DF}}\left({t}_{0},t\right)={E}_{\mathit{an}}\left({t}_{0},t\right) only if {\eta}_{G}={\eta}_{\mathit{CK}}={\eta}_{a}. Thus, we deduced that the O_{2}deficit model for the anaerobic energy was equal to (14b) if the chemical coupling efficiencies were alike.
Model 1 is a general version of the O_{2}deficit model. As the O_{2}deficit model has been widely debated, we proposed an alternative in Model 2 where the lactic power and energy were the same as for Model 1. However, we hypothesized that the alactic power was proportional to the rate of change of the aerobic power [24, 25]. Thus, for Model 2:
where θ is a constant of proportionality parameter that we determine below. Equation (15) gives that:
If the virtual metabolic power is below Q_{ max } and {\overline{Q}}_{a}={\overline{Q}}_{\mathit{vir}}, we achieved from equations (5) and (16) that:
The alactic energy can be found by calculating the area between {\overline{Q}}_{a} and Q_{ a } in a power time diagram. This area is multiplied with θ/τ_{ a } to find the alactic energy used. θ/τ_{ a } can be considered as the effectiveness of alactic ATP production relative to the aerobic ATP production. Thus we set \theta /{\tau}_{a}={\eta}_{a}/{\eta}_{\mathit{CK}}.
The mechanical efficiency was assumed to be η ≈ 0.5. For the chemical efficiency related to aerobic or lactic power {\eta}_{G}^{}={\eta}_{a}^{}\text{are}\text{between}\mathrm{0:3}\text{and}\mathrm{0:7}. We used {\eta}_{G}={\eta}_{a}\approx 0.6. This gave \eta \times {\eta}_{G}=\eta \times {\eta}_{a}=0.3. Since c was around 6.6 we found that {\eta}_{r}=1/\left(c\times \eta \times {\eta}_{G}\right)=1/\left(6.6\times 0.5\times 0.6\right)=0.51. GonzalesAlonso et al. [26] concluded from experimental data that the heat per use of ATP was around two times larger for oxidative phosphorylation and anaerobic glycolysis compared to ATP from CK. For the alactic power we therefore forecasted that η_{ CK } = 0.95. This gave that \theta =\left({\eta}_{a}/{\eta}_{\mathit{CK}}\right){\tau}_{a}=20\phantom{\rule{0.5em}{0ex}}s.
It has been shown that during intensities above the lactate threshold, ATP utilization increases and mechanical efficiency decreases at constant work rates. This may be explained by a change in fiber type recruitment, an elevated temperature, lowered pH or increased Pi levels [30]. We forecasted that τ_{ a }, η, η_{ a }, η_{ G } and η_{ CK } mainly depend on the lactate concentration and change in muscle pH, and assumed that a lower mechanical efficiency of muscle contractions and a lower P/O_{2} ratio could explain the drift in the steady state ATP consumption for a given work rate. To account for time varying parameters see the model in Appendix B.
Increased fitness after exercise is developed during the recovery period after the exercise. We defined an exercise which starts at time t_{0} and ends at time t_{1}, and a subsequent recovery period which starts at time t_{1} and ends at time t_{2}. Time t_{2} is defined as {Q}_{a}\left({t}_{2}\right)={Q}_{a}\left({t}_{0}\right) and C(t_{0}) = C(t_{2}). This means that {E}_{G}\left({t}_{0},{t}_{2}\right)=0 and {E}_{G}\left({t}_{0},{t}_{1}\right)={E}_{G}\left({t}_{1},{t}_{2}\right). We defined −E_{ DF }(t_{1}, t_{2}) as the energy depth E_{ DB }(t_{1}, t_{2}), to read:
We also defined the socalled excess post–exercise oxygen consumption (EPOC). In general E_{ DB }(t_{1}, t_{2}) is different from EPOC(t_{1}, t_{2}) unless {\overline{Q}}_{\mathit{vir}}={\overline{Q}}_{a}. However, the latter is usually the case during recovery.
For Model 1 in equations (13) we achieved:
In this case the alactic energy used from time t_{0} to time t_{1} for Model 1 was therefore given by:
For Model 2 we have:
We defined EPOC_{ Alt } as an alternative variant of the excess postexercise oxygen consumption:
Thus, finally:
Here, the alactic energy used from time t_{0} to time t_{1}, could be found by measuring the aerobic power at time t_{1} and t_{2}. The respective aerobic powers could be subtracted and then multiplied with θ to find the alactic energy.
It was also of interest to compare Model 1 and Model 2 more directly when Q_{ vir } = Q_{ a } and {\overline{Q}}_{\mathit{vir}}={\overline{Q}}_{a} . When using equation (5) this gave:
From equations (20)(23) this gave:
These two methods gave different results due to fact that the second term for {E}_{\mathit{CK}}^{1}\left({t}_{0},{t}_{1}\right) lacks in Model 2.
Hill et al. [34–37] have explained that the oxidative removal of lactate increases aerobic power. Using equation (25a) we achieved that:
If η_{ G } = η_{ a }, EPOC(t_{1}, t_{2}) corresponds to the total anaerobic energy, i.e. lactic and alactic anaerobic energy. Indeed, Margaria et al. [38] later on modified the concept of Hill et al. [34–37], and suggested that the increased aerobic power consisted of the fast alactic component and the slower lactic component. Finally, Gaesser and Brooks [39] introduced the term “excess postexercise oxygen consumption”, which also included the much more prolonged increase in aerobic power that is observed for hours after exercise. Model 1 did not account for this very slow component that lasts for hours. Model 2 in equation (25b) said that EPOC corresponded to the alactic component. However, the development of the energy in equation (25b) depended of the time dynamics for the aerobic power modeled according to equation (4).
Experimental tests
The derived mathematical simulations were compared with experimental data from an elite crosscountry skier while running with poles on a treadmill. The mass of the skier was m = 77.5 kg in all tests. All treadmill tests were performed on a 6 × 3 m motordriven treadmill (Bonte Technology, Zwolle, The Netherlands). Inclination and speed were calibrated using the Qualisys Pro Reflex system and the Qualisys Track Manager software (Qualisys AB, Gothenburg, Sweden). The treadmill belt consisted of a nonslip rubber surface that allowed the skier to use his own poles (pole length: 80% of body height) with special carbide tips. Gas exchange values were measured by opencircuit indirect calorimetry using an Oxycon Pro apparatus with a mixing chamber (Jaeger GmbH, Hoechberg, Germany). Before each measurement, the VO_{2} and VCO_{2} gas analyzers were calibrated using highprecision gases (16.00 ± 0.04% O_{2} and 5.00 ± 0.1% CO_{2}, RiessnerGase GmbH & co, Lichtenfels, Germany). The inspiratory flow meter was calibrated with a 3 L volume syringe (Hans Rudolph Inc., Kansas City, MO). Heart rate (HR) was measured with a heart rate monitor (Polar S610, Polar Electro OY, Kempele, Finland), using a 5 s interval for data storage. Blood lactate concentration (BLa) was measured on 5 μL samples taken from the fingertip by a Lactate Pro LT1710 t (ArkRay Inc, Kyoto, Japan).
Eight experimental protocols were performed, each on separate days with a minimum of 48 h between. The order of tests was performed as presented below. In order to investigate whole body exercise, running with poles was used in all tests. Before each testing session a standardized, testspecific 20min warmup was performed. Training on the days before testing was standardized, and the subject drank a standard fluid with sugar and electrolytes during all breaks while testing.
On the first test day, the skier performed six 5 min bouts with constant work rates at 0.105 inclines in radians. Five speeds at 0.25 m/s intervals below the lactate threshold were chosen, starting at 7.5 km/h = 2.08 m/s, followed by 2.33 m/s, 2.58 m/s, 2.83 m/s and 3.08 m/s. The sixth speed was increased by 0.125 m/s, giving 3.19 m/s, which is slightly above the lactate threshold. The starting speed was chosen based on experience from earlier tests of this athlete. 5min bouts were used to obtain steady state conditions. A 2min break with lowintensity walking was mandatory between each of the exercise bouts. Gas exchange values and heart rates were determined by the average of the last minute during each stage and blood lactate was measured directly after finishing each stage. The reason for using the last minute to assess respiratory variables was that the athletes are able to keep a more steady technique and physiological responses after 34 min.
In followup tests, we aimed to investigate whether our simulations fitted experimental data when varying the exercise time and during recovery after highintensity exercise. Also, velocities slightly above or below the lactate threshold were tested in order to test the sensitivity of our model for lactate. Thus, the skier performed the following six tests on separate days with at least 2 days in between each test, at 0.105 radians incline with blood lactate measured directly after each stage:

2.87 m/s was performed for 800 s.

3.08 m/s was performed over 800 s, followed by a 5min recovery phase at 1.67 m/s with blood lactate measured after 2 and 5 min. Additionally, a 3.05 m/s stage was performed to exhaustion.

3.08 m/s was performed over 2000 s.

3.19 m/s stages were performed for 400 s and to exhaustion.

A 3.33 m/s stage was performed for 400 s.

3.88 m/s stages were performed for 150 s, 200 s and to exhaustion.
Finally, maximal metabolic power was tested on a separate day at an incline of 0.105 radians with a starting speed of 3 m/s. The speed was increased by 0.3 m/s every minute until exhaustion. VO_{2} was measured continuously, and the average of the three highest 10 s consecutive measurements determined VO_{2}max and was used to calculate the maximal metabolic power. The test was considered to be a maximal effort if the following three criteria were met: 1) a plateau in VO_{2} was obtained with increasing exercise intensity, 2) respiratory exchange ratio above 1.10, and 3) blood lactate concentration exceeding 8 mmol/L.
Results
Figure 2 illustrates that the experimental steady state blood lactate values showed good similarity with the simulated values based on Moxnes and Hausken’s [10] one compartment model for the average concentration in the total lactate pool. The blood lactate concentration equals the average concentration during steady state. At high lactate concentrations, the drain of lactate is saturated, which can be modeled by the function Tanh\left(\chi \phantom{\rule{0.5em}{0ex}}C\left(t\right)\right). No saturation was achieved when χ→0, since \underset{\chi \to 0}{Lim}\phantom{\rule{0.5em}{0ex}}Tanh\left(\chi \phantom{\rule{0.5em}{0ex}}C\left(t\right)\right)/\chi =C\left(t\right). The lactate threshold Q_{ LT } was achieved when 1\frac{\chi {p}_{0}}{{d}_{0}\left({Q}_{\mathit{max}}{Q}_{a}\right)}=0\Rightarrow {Q}_{a}={Q}_{\mathit{LT}}={Q}_{\mathit{max}}\frac{\chi {p}_{0}}{{d}_{0}},\frac{{Q}_{\mathit{LT}}}{{Q}_{\mathit{max}}}=1\frac{\chi {p}_{0}}{{d}_{0}{Q}_{\mathit{max}}}. Q_{ max } was known and we used the numerical value of p_{0} from Moxnes and Hausken [10]. Thus, we only needed to fit the two parameters χ/d_{0} and χ to the steady state measured values.
In the first step, we applied visual curve fitting, which means that we chose plausible values for the two parameters. We then plotted and compared the solution visually with the experimental data. The values of the two parameters were changed repeatedly until a good visual fit was obtained, while ensuring that the parameters had physiologically trustworthy numerical values. In the second step, a least square fit to the data was performed to produce ex post best fit estimates of these two parameters using the visual estimates as starting guess points and choosing a range around each starting point of the parameters. The method was performed separately for each parameter keeping the other parameter fixed. Steps 1 and 2 were repeated sufficiently until we were certain that we had obtained the optimal values for each of the two parameters. For robustness in the calculations, we applied the steady state solution. If we applied that {C}_{s}\approx 1+\frac{{p}_{0}}{{d}_{0}{Q}_{\mathit{max}}\left(1{Q}_{a}/{Q}_{\mathit{max}}\right)}, the parameter χ got removed. The lactate threshold became Q_{ a } = Q_{ max }. χ determined the lactate threshold by the exact relation. The rest of the lactate curve was dependent on \frac{{p}_{0}}{{d}_{0}{Q}_{\mathit{max}}\left(1{Q}_{a}/{Q}_{\mathit{max}}\right)}. This gave a linear dependency of the lactate curve on p_{0}/d_{0} and an inverse dependency of the lactate curve on Q_{ max } − Q_{ a }.
Figure 3 demonstrates the measured and simulated lactate concentration during steady state work rate. In this simulation, we only applied visual fitting to determine the d_{0} parameter, and found that the simulated lactate concentration was somewhat higher than the measured (blood lactate) concentration for intensities above the lactate threshold. During recovery, the measured blood lactate concentration was higher than the simulated lactate concentration. Thus, we also calculated values according to Moxnes and Sandbakk [22] where the lactate concentration was separated into two compartments: muscles + internal organs and blood. In that case, Figure 3 shows that the simulation results are more similar to the experimental results. However, some discrepancies still appear for the higher work rates.
Figure 4 shows the different anaerobic alactic powers for intensities below the lactate threshold when using Models 1 and 2. The calculated alactic power for Model 1 showed a somewhat complicated behavior since it also depended on the model for the lactic power. A special feature was the overshoot that applied during recovery because the lactic power becomes negative. For Model 2 the lactic power was only different from zero when the rate of change of the aerobic power was different from zero.
Figure 5 illustrates the lactic and alactic storage of energy. We found different behaviors for Models 1 and 2. For Model 1, the alactic energy storage started to increase from a local minimum at the time the aerobic power reached a steady state. The aerobic power reached steady state before the lactic power did for v_{4} and v_{5}. When the exercise was terminated for v_{6}, v_{7} and v_{8}, the alactic energy was fully restored. The alactic energy resources increased due to use of lactic power. Although the details of the time history of alactic energy depend on the chemical coupling coefficients, the overall scenario was much the same for different numerical values of the coefficients. It is important to note that the alactic power depended on the model for the lactic power. When O_{2} was excluded, the PCr store was rebuilt by anaerobic glycogenolysis or glycolysis [15–17]. In Model 2, the alactic energy showed a simpler behavior than in Model 1 since the alactic power only depended on the rate of change of the aerobic power. When comparing the solution in Figure 5 in the current paper with the invasive measurements shown in Figure 2 in the paper Jeneson et al. [23] it appeared that Model 2 gave a better fit than Model 1.
To study the kinetics of anaerobic powers more in detail we compared interval exercise with continuous exercise by using Models 1 and 2. These were simulated exercises matched for energy expended. Interval exercise changed in intensity with durations above and below the lactate threshold as described below. The continuous exercise maintained steady state intensity. The two exercises were given by:
T is the time variable, 0<=t<=T, and T = 360 s is a parameter for the time duration, i.e. how long the exercise lasted. Figure 6 shows the metabolic powers for interval and continuous exercises as functions of time, whereas Figure 7 shows the simulated lactate concentration for the two different exercises. The lactate values were reduced to the initial levels between intervals, whereas lactate continued to rise during continuous exercise. Figure 8 shows the lactic power and the alactic powers calculated by Models 1 and 2. Overall, the alactic energy stores were rebuilt during recovery with interval exercise, which allowed for the alactic energy to be used at the beginning of each interval. Figure 9 shows the lactic and alactic energy utilization and shows that during exercise the alactic energy storage decreased with time for Model 2. For Model 1 a local minimum appeared around the time when the aerobic power reached steady state. During interval exercise, the alactic energy storage as a function of time did not differ much between the two models.
Conclusions
The current study investigated two different mathematical models for the kinetics of anaerobic power during whole body exercise at different exercise intensities. The results indicate that blood lactate levels can be accurately modeled during steady state, and suggest a linear relationship between the alactic anaerobic power and the rate of change of the aerobic power. Overall, we propose that the current simulation models provide useful insight into how the anaerobic powers during practical training and testing should be interpreted.
When comparing the experimental blood lactate values with Moxnes and Hausken’s [10] model for lactic power, we found similar results during steady state exercise intensities. Uncertainties of less than ± 0.5 mmol/L using these mathematical models are indicated. For intensities above the lactate threshold, the measured blood lactate values were significantly lower than those calculated by the Moxnes and Hausken [10] model, whereas during recovery after high intensity exercise, the blood lactate concentration was higher than calculated. This could be explained by a possible delay in blood lactate concentration compared to the average concentration in the total lactate pool that applies theoretically. Therefore, we compared the experimental results with a twocompartment model by Moxnes and Sandbakk [22], where lactate concentration was separated into muscle and blood compartments. The employment of this model gave a better fit between simulated and experimental results, especially for the highest work rates and during recovery. However, some discrepancies still appeared for the higher work rates, which is a topic for further research.
The calculation of alactic storage of energy in Model 1, which is a generalization of the O_{2}deficit model, showed that the alactic energy storage increased due to use of lactic power from the time the aerobic power reached a steady state. However, the various inferences about the alactic storage of energy during steady state work rate have some shortcomings. For example, the results may have been influenced by the fact that the O_{2}deficit is composed of both lactic and alactic metabolisms. Thus, no definite conclusion could be drawn for Model 1 for alactic power. For Model 2, the alactic energy showed a simpler behavior since it only depended on the rate of change of the aerobic power. When comparing Models 1 and 2 with the invasive measurements by Jeneson et al. [23], we concluded that Model 2 gave the best behavior. Thus, a linear relationship between the alactic anaerobic power and the rate of change of the aerobic power was suggested. In future research, it remains a challenge to confirm the contribution of lactic and alactic anaerobic metabolism and power.
The simulations of interval exercises showed a drop in lactate levels between the intervals, and that the alactic energy stores were rebuilt during recovery. Practically, this allows athletes to work at high exercise intensities for longer time durations without lowering the pH values in the muscles compared to continuous exercise at similar intensities. This might be one reason for the effectiveness of interval exercise, as verified both by the scientific literature [40] and by elite athletes’ practical training programs where interval exercises are an important component. When the alactic anaerobic powers during continuous exercises from Models 1 and 2 were compared, Model 1 showed a local minimum appearing around the time when the aerobic power reached steady state for continuous exercise, whereas for Model 2 the alactic energy storage decreased with time. During interval exercise, the alactic energy storage as a function of time was similar for the two models.
Modeling in human biology is a challenge since one is confronted with conceiving a simple, but realistic representation of complex phenomena. The parameter values used in the current study will be dependent on the fitness level of the individual being tested and on the concentration of glycogen in the muscles. To construct valid simulation models during dynamic exercise such models need to be developed for each individual.
Appendix A: The lactate model
In equation (25) of Moxnes and Hausken [10] we find that:
where p_{1}() and d_{1}() are dimensionless functions and α is a parameter to describe the saturation of the concentration of lactate C(t). We use \frac{Tanh\left(\chi \phantom{\rule{0.12em}{0ex}}C\left(t\right)\right)}{\chi} instead of \frac{\left(1Exp\left(\alpha \phantom{\rule{0.12em}{0ex}}C\left(t\right)\right)\right)}{\alpha} to describe saturation. To read:
We set that:
and β_{0} is a parameter that must be fitted to the experimental data, as earlier done by Moxnes and Sandbakk [22] when studying a two compartment model of lactate. The steady state solution is given by:
We use 1/\left({\beta}_{0}\phantom{\rule{0.12em}{0ex}}{p}_{0}\right)=0.6\phantom{\rule{0.12em}{0ex}}{Q}_{\mathit{max}}, where 0.6Q_{ max } is around peak fat metabolism. All together we set:
Appendix B: Model for parameters varying with the lactate level
To account for situations with varying τ_{ a }, η_{ a }, η_{ G } and η_{ CK } we can use a first order differential equation of the virtual aerobic metabolic power, to read:
Equation (B1) allows the calculation of the aerobic power when the efficiencies and the time parameter τ_{ a } are time dependent. This might appear due to changes in lactate levels. We assume that η × η_{ a } is constant through time. This gives from (B1):
References
Cabrera ME, Saidel GM, Kalhan SC: Lactate metabolism during exercise: analysis by an integrative systems model. Am J Physiol. , 277: R1522R1536.
Stirling JR, Zakynthinaki MS, Saltin B: A model of oxygen uptake kinetics in response to exercise: Including a means of calculating oxygen demand/deficit/debt. Bull Math Biol. 2005, 67: 9891015. 10.1016/j.bulm.2004.12.005.
Murias JM, Spencer MD, Kowalchuk JM, Paterson DH: Influence of phase I duration on phase II VO2 kinetics parameters estimates in older and young adults. Am J Physiol Regul Integr Comp Physiol. 2011, 301: R218R224. 10.1152/ajpregu.00060.2011.
Billat VL, Bocquet V, Slawinski J, Laffitte L, Demarle A, Chassaing P, Koralsztein JP: Effects of prior intermitted runs av v VOmax on oxygen kinetics during an allout severe run in humans. J Sports Med Phys Fitness. 2000, 40: 185194.
Gaesser GA, Poole DC: The slow component of oxygen uptake in humans. Exerc Sport Sci Rev. 1996, 24: 3570.
Brooks GA: Lactate production under fully aerobic conditions: the lactate shuttle during rest and exercise. Fed Proc. 1985, 45: 29242929.
Brooks GA: Current concepts in lactate exchange. Med Sci Sports Exerc. 1991, 23: 895906.
Brooks GA: Lactate shuttles in nature. Biochem Soc Trans. 2002, 30: 258264.
Brooks GA: Link between glycolytic and oxidative metabolism. Sports Med. 2007, 37: 341343. 10.2165/0000725620073704000017.
Moxnes JF, Hausken K: A mathematical model for the training impulse and lactate influx and outflux during exercise. Int J Mod Phys C. 2009, 20 (1): 147177. 10.1142/S0129183109013522.
Carlson FD, Siger A: The mechanochemistry of muscular contraction I. The isimetric twitch. J Gen Physiol. 1960, 44: 3360. 10.1085/jgp.44.1.33.
Rossiter HB, Ward SA, Doyle VL, Howe FA, Griffiths JR, Whipp BJ: Interference from pulmonary O2 uptake with respect to intramuscular (phosphcreatine) kinetcs during moderate exercise in humans. J Physiol. 1999, 518: 921932. 10.1111/j.14697793.1999.0921p.x.
Margaria R: Biomechanics and energetic of muscular exercise. 1976, Oxford: Oxford University Press
Morton RH: A three component model of human bioenergetics. J Math Biol. 1986, 24: 451466. 10.1007/BF01236892.
Crowther GJ, Kemper WF, Carey MF, Conley KE: Control of glycolysis in contracting skeletal muscle. II. Turning it off. Am J Physiol Endocrinol Metab. 2002, 282: 7479.
Jubrias SA, Esselman PC, Price LB, Cress ME, Conley KE: Large energetic adaptations of elderly muscle to resistance and endurance training. J Appl Physiol. 2001, 90: 16631670. 10.1063/1.1374452.
Lanza IR, Wigmore DM, Befroy DE, KentBraun JA: In vivo ATP production during freeflow and ischaemic muscle contractions in humans. J Physiol. 2006, 577: 353367. 10.1113/jphysiol.2006.114249.
Medbø JI, Mohn AC, Tabata I, Bahr R, Vaage O, Sejersted OM: Anaerobic capacity determined by maximum accumulated O2 deficit. J Appl Physiol. 1988, 64: 5060.
Bangsbo J: Quantification of anaerobic energy production during intense exercise. Med Sci Sports Exerc. 1998, 30: 4752.
Gastin PB: Energy system interaction and relative contribution during maximal exercise. Sports Med. 2001, 31: 725741. 10.2165/0000725620013110000003.
Noordhof DA, de Kroning JJ, Foster C: The maximal accumulated oxygen deficit method: a valid and reliable measure of anaerobic capacity?. Sports Med. 2010, 40: 285302. 10.2165/1153039000000000000000.
Moxnes JF, Sandbakk K: The kinetics of lactate production and removal during wholebody exercise. Theor Biol Med Model. 2012, 9: 710.1186/1742468297.
Jeneson JAL, Schmitz JPJ, van den Broek NMA, van Riel AAW, Hilbers PAJ, Nicolay K, Prompers JJ: Magnitude and control of mitocondrial sensitivity to ADP. Am J Physiol Endocrinol Metabol. 2009, 297: E774E784. 10.1152/ajpendo.00370.2009.
Mahler M: First order kinetics of muscle oxygen consumption, and an equivalent proportionality between Qo2 and Phosphorylcreatine level. J Gen Physiol. 1985, 86: 135165. 10.1085/jgp.86.1.135.
Meyer RA: A linear model of muscle respiration explains monoexponential phosphocreatine changes. J Appl Physiol. 1988, 75: 648656.
GonzalezAlonso J, Quistorff B, Krustrup P, Bangsbo J, Saltin B: Heat production in human skeletal muscle at the onset of intense dynamic exercise. J Physiol. 2000, 524 (2): 603615. 10.1111/j.14697793.2000.00603.x.
di Prampero PE, Ferretti G: The energetics of anaerobic muscle metabolism: a reappraisal of older and recent concepts. Respir Physiol. 1999, 118 (2–3): 103115.
di Prampero PE: Factors limiting maximal performance in humans. Eur J Physiol. 2003, 90: 420429. 10.1007/s004210030926z.
Losnegard T, Myklebust H, Hallen J: Energy system contribution as a determinant of performance in elite skiers. 2011, 10.1249/MSS.0b013e3182388684.
Bangsbo J, Krustrup P, GonzalezAlonso J, Saltin B: ATP production and efficiency of human skeletal muscle during intense exercise: effects of previous exercise. Am J Physiol Endocrinol Metabol. 2001, 280: E956E964.
Bertuzzi RC, Franchini E, Ugrinowitsch C, Kokubun E, LimaSilva AE, Pires FO, Nakamura FY, Kiss MA: Predicting MAOD using only a supramaximal exhaustive test. Int J Sports Med. 2010, 31 (7): 477481. 10.1055/s00301253375.
Cerretelli P, Pendergast DR, Paganelli WC, Rennie DW: Effects of specific muscle training on VO2on response and early blood lactate. J Appl Physiol. 1979, 47: 761769.
Binzoni T, Ferretti G, Sechenker K, Cerretelli P: Phosphocreatine hydrolysis in 31PNMR at the onset of constantload exercise in humans. J Appl Physiol. 1992, 73: 16441649.
Hill VH, Lupton H: Muscular exercise, lactic acid, and the supply and utilization of oxygen. Q J Med. 1923, 16: 135171. 10.1093/qjmed/os16.62.135.
Hill AV, Long CNH, Lupton H: Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts IIII. Proc Roy Soc. 1924, 96: 438475. 10.1098/rspb.1924.0037.
Hill AV, Long CNH, Lupton H: Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts IVVI. Proc Roy Soc. 1924, 97: 84138. 10.1098/rspb.1924.0045.
Hill AV, Long CNH, Lupton H: Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts VIIVIII. Proc Roy Soc. 1924, 97: 155176. 10.1098/rspb.1924.0048.
Margaria R, Edwards HT, Dill OB: The possible mechanisms of contracting and paying the oxygen dept and the role of lactic acid in muscular contraction. Am J Physiol. 1933, 106: 689715.
Gaesser GA, Brooks GA: Metabolic bases of excess postexercise oxygen consumption: a review. Med Sci Sports Exerc. 1984, 16: 2943.
Laursen PB, Jenkins DG: The scientific basis for highintensity interval training: Optimising training programmes and maximizing performance in highly trained endurance athletes. Sports Med. 2002, 32: 5373. 10.2165/0000725620023201000003.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
JM performed the mathematical simulations, ØS performed all laboratory testing and KH supported in finishing the paper in its final form. All authors contributed with important intellectual content in all parts of the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Moxnes, J.F., Hausken, K. & Sandbakk, Ø. On the kinetics of anaerobic power. Theor Biol Med Model 9, 29 (2012). https://doi.org/10.1186/17424682929
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/17424682929
Keywords
 Blood lactate
 Crosscountry skier
 Mathematical model
 Oxygen uptake
 Phosphocreatine