 Research
 Open Access
On the kinetics of anaerobic power
 John F Moxnes^{1},
 Kjell Hausken^{2} and
 Øyvind Sandbakk^{3}Email author
https://doi.org/10.1186/17424682929
© Moxnes et al.; licensee BioMed Central Ltd. 2012
 Received: 13 March 2012
 Accepted: 20 June 2012
 Published: 25 July 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.
Keywords
 Blood lactate
 Crosscountry skier
 Mathematical model
 Oxygen uptake
 Phosphocreatine
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
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}$.
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.
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 “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 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}}$.
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].
Experimentally, E_{ G }(t_{0}, t) can be found by measuring changes in the lactate concentration before and after exercise.
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 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.
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.
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.
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.
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.
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.
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.

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
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 }.
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
Appendix B: Model for parameters varying with the lactate level
Declarations
Authors’ Affiliations
References
 Cabrera ME, Saidel GM, Kalhan SC: Lactate metabolism during exercise: analysis by an integrative systems model. Am J Physiol. , 277: R1522R1536.PubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 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.PubMedGoogle Scholar
 Gaesser GA, Poole DC: The slow component of oxygen uptake in humans. Exerc Sport Sci Rev. 1996, 24: 3570.View ArticlePubMedGoogle Scholar
 Brooks GA: Lactate production under fully aerobic conditions: the lactate shuttle during rest and exercise. Fed Proc. 1985, 45: 29242929.Google Scholar
 Brooks GA: Current concepts in lactate exchange. Med Sci Sports Exerc. 1991, 23: 895906.View ArticlePubMedGoogle Scholar
 Brooks GA: Lactate shuttles in nature. Biochem Soc Trans. 2002, 30: 258264.View ArticlePubMedGoogle Scholar
 Brooks GA: Link between glycolytic and oxidative metabolism. Sports Med. 2007, 37: 341343. 10.2165/0000725620073704000017.View ArticlePubMedGoogle Scholar
 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.View ArticleGoogle Scholar
 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.PubMed CentralView ArticlePubMedGoogle Scholar
 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.PubMed CentralView ArticlePubMedGoogle Scholar
 Margaria R: Biomechanics and energetic of muscular exercise. 1976, Oxford: Oxford University PressGoogle Scholar
 Morton RH: A three component model of human bioenergetics. J Math Biol. 1986, 24: 451466. 10.1007/BF01236892.View ArticlePubMedGoogle Scholar
 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.Google Scholar
 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.View ArticlePubMedGoogle Scholar
 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.PubMed CentralView ArticlePubMedGoogle Scholar
 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.PubMedGoogle Scholar
 Bangsbo J: Quantification of anaerobic energy production during intense exercise. Med Sci Sports Exerc. 1998, 30: 4752.View ArticlePubMedGoogle Scholar
 Gastin PB: Energy system interaction and relative contribution during maximal exercise. Sports Med. 2001, 31: 725741. 10.2165/0000725620013110000003.View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 Moxnes JF, Sandbakk K: The kinetics of lactate production and removal during wholebody exercise. Theor Biol Med Model. 2012, 9: 710.1186/1742468297.PubMed CentralView ArticlePubMedGoogle Scholar
 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.View ArticleGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 Meyer RA: A linear model of muscle respiration explains monoexponential phosphocreatine changes. J Appl Physiol. 1988, 75: 648656.Google Scholar
 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.PubMed CentralView ArticlePubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 di Prampero PE: Factors limiting maximal performance in humans. Eur J Physiol. 2003, 90: 420429. 10.1007/s004210030926z.View ArticleGoogle Scholar
 Losnegard T, Myklebust H, Hallen J: Energy system contribution as a determinant of performance in elite skiers. 2011, 10.1249/MSS.0b013e3182388684.Google Scholar
 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.Google Scholar
 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.View ArticlePubMedGoogle Scholar
 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.PubMedGoogle Scholar
 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.PubMedGoogle Scholar
 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.View ArticleGoogle Scholar
 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.View ArticleGoogle Scholar
 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.View ArticleGoogle Scholar
 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.View ArticleGoogle Scholar
 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.Google Scholar
 Gaesser GA, Brooks GA: Metabolic bases of excess postexercise oxygen consumption: a review. Med Sci Sports Exerc. 1984, 16: 2943.PubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.