On the kinetics of anaerobic power.  
Jump to Full Text  
MedLine Citation:

PMID: 22830586 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

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. 
Authors:

John F Moxnes; Kjell Hausken; Øyvind Sandbakk 
Related Documents
:

7882126  Adherence to pelvic floor muscle exercise and longterm effect on stress urinary incont... 2938306  Noninvasive screening procedure for evaluation of bladder emptying function. 10491346  Nonsteroidal antiinflammatory drugs and renal response to exercise: a comparison of i... 16217326  Factors that distinguish continent from incontinent older adults with detrusor overacti... 20930016  Comparison of 6min walk test distance and cardiopulmonary exercise test performance in... 21318926  Changing exercise through targeting affective or cognitive attitudes. 
Publication Detail:

Type: Journal Article Date: 20120725 
Journal Detail:

Title: Theoretical biology & medical modelling Volume: 9 ISSN: 17424682 ISO Abbreviation: Theor Biol Med Model Publication Date: 2012 
Date Detail:

Created Date: 20130121 Completed Date: 20130624 Revised Date: 20130712 
Medline Journal Info:

Nlm Unique ID: 101224383 Medline TA: Theor Biol Med Model Country: England 
Other Details:

Languages: eng Pagination: 29 Citation Subset: IM 
Affiliation:

Department of Human Movement Science, Norwegian University of Science and Technology, 7491, Trondheim, Norway. 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Anaerobiosis Biomechanics* Computer Simulation Exercise* Humans Kinetics Male Models, Biological* Muscle, Skeletal / physiology* 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): Theor Biol Med Model Journal ID (isoabbrev): Theor Biol Med Model ISSN: 17424682 Publisher: BioMed Central 
Article Information Download PDF Copyright ©2012 Moxnes et al.; licensee BioMed Central Ltd. openaccess: Received Day: 13 Month: 3 Year: 2012 Accepted Day: 20 Month: 6 Year: 2012 collection publication date: Year: 2012 Electronic publication date: Day: 25 Month: 7 Year: 2012 Volume: 9First Page: 29 Last Page: 29 PubMed Id: 22830586 ID: 3548746 Publisher Id: 17424682929 DOI: 10.1186/17424682929 
On the kinetics of anaerobic power  
John F Moxnes1  Email: johnf.moxnes@ffi.no 
Kjell Hausken2  Email: kjell.hausken@uis.no 
Øyvind Sandbakk3  Email: oyvind.sandbakk@svt.ntnu.no 
1Department for Protection, Norwegian Defence Research Establishment, P.O. Box 25, 2007, Kjeller, Norway 

2Faculty of Social Sciences, University of Stavanger, 4036, Stavanger, Norway 

3Department of Human Movement Science, Norwegian University of Science and Technology, 7491, Trondheim, Norway 
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.
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).
ATP is produced by three different sources:
(1a)
Qtmetabolicpower−Qrrestingpower=modelQawtaerobicpower+QGtpowerfromG+QCKt︸powerfromCK,a 
(1b)
ItrateofATPconsumed−Irrestingrate=modelIawtrateofATPconsumedwithATPproducedaerobically+IGtrateofATPconsumedwithATPproducedfromG+ICKt︸rateofATPconsumedwithATPproducedfromCK 
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. QG is ATP production by glycolysis/glycogenolysis and CK by phosphocreatine break down to creatine. Qaw is aerobic power due to internal and external work. The aerobic power is Qa=Qr+Qaw.
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:
(2a)
PaWorkpower=defηηaQaw,HawHeatpower=Qaw−Pa=Qaw1−ηηa 
(2b)
PGWorkpower=defηηGQG,HGHeatpower=QG−PG=QG1−ηηG 
(2c)
PCKWorkpower=defηηCKQCK,HCKHeatpower=QCK−PCK=QCK1−ηηCK 
(2d)
PTTotalworkpower=defPa+PG+PCK=ηηaQaw+ηηGQG+ηηCKQCK, 
(2e)
HHeatpowerduetowork=defHaw+HG+HCK 
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 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 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 Q¯vir which is a straight line that crosses Q_{max}. Thus we set:
(3a)
Q¯vir−Qr=modelQun+cP,Q¯a=MinQmax,Q¯vir, 
(3b)
P¯vir=defQ¯vir−Qrηaη=Qun+cPηaη,ηr=def1/cηaη 
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. P¯vir includes both internal and external work power. Q¯vir is input to the further modeling. We set that Q¯vir is a linear interpolation function through our data points. For values above or below the maximum and minimum measured Q¯vir we used linear extrapolation. During steady state, Qaw=Q¯a−Qr=MinQmax−Qr,Q¯vir−Qr.
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, Qat=modelMinQmax,Qvirt. To account for the delay mathematically, we used a first order differential equation of the virtual power, to read:
(4a)
Q·virt=modelQ¯vir−Qvirtτa,τa=30s,Qvirt0=Qr 
(4b)
Qa=modelMinQmax,Qvirt 
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 Qa=MinQmax,Qvir=Qvir. This gives from equation (4a and b) that:
(5)
Q·at=Q¯vir−Qatτa 
For steady work rate equation (5) gives the analytical solution:
(6)
Qat=Q¯vir−Q¯vir−Qat0Exp−t−t0/τa 
If the work rate is sufficiently low, Q¯a=Q¯vir. Notable is that solution (6) will only apply for a restricted time period unless Q¯a=Q¯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:
(7)
QGt=modelm×λ×C·t,λ=3×20:J/kgmmol/L,λ=3ml/kgO2,1mlO2=20J 
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:
(8)
EGt0,t=def∫t0tQGt′dt′=mλ∫t0tC·t′dt′=mλCt−Ct0 
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:
(9)
Qaw+QG+QCK=1ηPa/ηa+PG/ηa+PCK/ηa+PG/ηG−PG/ηa+PCK/ηCK−PCK/ηa=PTηηa+PGη1/ηG−1/ηa︸QGt1−ηa/ηa+PCKη1/ηG−1/ηa︸QCKt1−ηCK/ηa⇒QCKt=ηaηCKPTtηηa−Qawt−ηGηaQGt 
The metabolic power is then given by:
(10)
Qt−Qr=Qaw+QG+QCK=Qaw+QG+ηaηCKPTηηa−Qaw−ηGηaQG=Qawt1−ηaηCK+QGt1−ηGηCK+ηaηCKPTtηηa 
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 C.=QG=Q.aw=0. Thus Qaw is a constant and Q−Qr=Qaw 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 PT=ηηaQaw. However, during steady state (below LT) Qaw=Qun+cP. Thus, below the lactate threshold PT=Qun+cPηηa=P¯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:
(11)
PTt=modelP¯vir=Q¯vir−Qrηaη=Qun+cPηaη 
Inserting (11) into (9)(10) gives the alactic power as:
(12a)
QCK1t=ηaηCKQ¯vir−Qr−Qat−Qr−ηGηaQGt=ηaηCKQ¯vir−Qat−ηGηaQGt 
and the metabolic rate as:
(12b)
Q1t−Qr=Qat−Qr1−ηaηCK+QGt1−ηGηCK+ηaηCKQ¯vir−Qr 
To calculate QCK1t for a work rate P(t) in (12a) we inserted Q¯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:
(13)
ECK1(t0,t)=def∫t0tQCK1(t')dt'=∫t0tηaηCKQ¯vir−Qa(t')−ηGηaQG(t')dt'=ηaηCKEDF(t0,t)−ηGηCKEG(t0,t)EDF(t0,t)=def∫t0tQ¯vir(t')−Qa(t')dt' 
The anerobic power and the anaerobic energy used from time t0 to time t for Model 1 are:
(14a)
Qan1(t)=defQCK1(t)+QG(t)=ηaηCKQ¯vir−Qa(t)−ηGηaQG(t)+QG(t)=ηaηCKQ¯vir−Qa(t)+1−ηGηaQG(t) 
and:
(14b)
Ean1(t0,t)=def∫t0tQan1(t')dt'=ηaηCKEDF(t0,t)+1−ηGηaEG(t0,t) 
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 EDFt0,t=Eant0,t only if ηG=ηCK=η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:
(15)
QCK2(t)=modelθQ.a(t) 
where θ is a constant of proportionality parameter that we determine below. Equation (15) gives that:
(16)
ECK2t0,t=def∫t0tQCK2t′dt′=θ∫t0tQ·at′dt′=θQat−Qat0 
If the virtual metabolic power is below Q_{max} and Q¯a=Q¯vir, we achieved from equations (5) and (16) that:
(17)
ECK2t0,t=θ∫t0tQ·at′dt′=θτa∫t0tQ¯a−Qat′dt′︸Area=θQat−Qat0 
The alactic energy can be found by calculating the area between 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 θ/τa=ηa/ηCK.
The mechanical efficiency was assumed to be η ≈ 0.5. For the chemical efficiency related to aerobic or lactic power ηG=ηaarebetween0:3and0:7. We used ηG=ηa≈0.6. This gave η×ηG=η×ηa=0.3. Since c was around 6.6 we found that ηr=1/c×η×ηG=1/6.6×0.5×0.6=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 θ=ηa/ηCKτa=20s.
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 Qat2=Qat0 and C(t_{0}) = C(t_{2}). This means that EGt0,t2=0 and EGt0,t1=−EGt1,t2. We defined −E_{DF}(t_{1}, t_{2}) as the energy depth E_{DB}(t_{1}, t_{2}), to read:
(18)
EDBt1,t2=def−EDFt1,t2=−∫t1t2Qat′−Q¯virdt′,EPOCt1,t2=def−∫t1t2Qat′−Q¯adt′ 
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 Q¯vir=Q¯a. However, the latter is usually the case during recovery.
For Model 1 in equations (13) we achieved:
(19)
ECK1t0,t2=ηaηCKEDFt0,t2−ηGηCKEGt0,t2︸=0=ηaηCKEDFt0,t1+ηaηCKEDFt1,t2=0⇒EDFt0,t1=−EDFt1,t2=EDBt1,t2 
In this case the alactic energy used from time t_{0} to time t_{1} for Model 1 was therefore given by:
(20)
ECK1t0,t1=ηaηCKEDBt1,t2+ηGηCKEGt1,t2=θτaEDBt1,t2+ηGηCKEGt1,t2 
For Model 2 we have:
(21)
ECK2t0,t2=∫t0t2QCK2t′dt′=θ∫t0t2Q·at′dt′=θ∫t0t1Q·at′dt′︸ECK2t0,t1+θ∫t1t2Q·at′dt′︸ECK2t1,t2=0⇒ECK2t0,t1=−ECK2t1,t2 
We defined EPOC_{Alt} as an alternative variant of the excess postexercise oxygen consumption:
(22)
EPOCAltt1,t2=def−τa∫t1t2Q·at′dt′=τaQat1−Qat2 
Thus, finally:
(23)
ECK2t0,t1=∫t0t1QCK2t′dt′=θ∫t0t1Q·at′dt′=θτa×τaQat1−Qat2=θτa×EPOCAltt1,t2 
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 Q¯vir=Q¯a . When using equation (5) this gave:
(24)
IfQvir=QaandQ¯vir=Q¯a,then EPOCAltt1,t2=−τa∫t1t2Q·at′dt′=−∫t1t2Qat′−Q¯virdt′=EDBt1,t2=−∫t1t2Qat′−Q¯adt′=EPOCt1,t2 
From equations (20)(23) this gave:
(25a)
ECK1t0,t1=ηaηCKEDBt0,t1+ηGηCKEGt1,t2=θτaEPOCt1,t2︸firstterm+ηGηCKEGt1,t2︸secondterm 
(25b)
ECK2t0,t1=θτaEPOCt1,t2 
These two methods gave different results due to fact that the second term for ECK1t0,t1 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:
(26)
Ean1t0,t1=∫t0t1Qan1t′dt′=ECK1t0,t1+EG1t0,t1=ηaηCKEDFt0,t1+1−ηGηaEGt0,t1=θτaEPOCt1,t2+ηGηa−1EGt1,t2 
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).
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.
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χCt. No saturation was achieved when χ→0, since Limχ→0TanhχCt/χ=Ct. The lactate threshold Q_{LT} was achieved when 1−χp0d0Qmax−Qa=0⇒Qa=QLT=Qmax−χp0d0,QLTQmax=1−χp0d0Qmax. 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 Cs≈1+p0d0Qmax1−Qa/Qmax, 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 p0d0Qmax1−Qa/Qmax. 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:
Interval:Q¯vir=Qmax/2+Qmax/2×Sin2πt/T,T=360sContinuous:Q¯vir={Qmax/2+Q5−Qmax/2×Sin2πt/Twhen0≤t≤1/4TQmax/2+Q5−Qmax/2when1/4T<t≤TQ5=0.87Qmax 
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.
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.
In equation (25) of Moxnes and Hausken [^{10}] we find that:
(A1)
C·t=p0Qatp1Qat−d0QatQmax−Qatd1Qat1−Exp−αCtα 
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 TanhχCtχ instead of 1−Exp−αCtα to describe saturation. To read:
(A2)
C·t=p0DQat−d0×TanhχCtχ×DQat×Qmax−QatDQat=defQatp1Qat=d1Qat 
We set that:
(A3)
DQa=modQa−α0β0p0Tanhβ0p0Qa 
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:
(A4)
Cs=12χLn1+p0χd0Qmax1−Qa/Qmax1−p0χd0Qmax1−Qa/Qmax,Limχ→0Cs=p0d0Qmax1−Qa/Qmax 
We use 1/β0p0=0.6Qmax, where 0.6Q_{max} is around peak fat metabolism. All together we set:
(A5)
Dx=x−α00.6QmaxTanhx0.6Qmax,1/β0p0=0.6Qmax,α0=0.9,Ct0=0.045kg/m3=0.5mmol/L,χ=2.54m3/kg,p0=10−5kg/(m3s)/(J/s),d0=1.1×10−7/(J/s)2/s,Qmax=1886J/s 
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:
(B1)
P·virt=modelmodP¯vir−Pvirtτa,Pa=modelMinηηaQmax−Qb,Pvirt,τa=30s,Pvirt0=0,Qa=Pa/ηηa 
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):
(B2)
P·virt/ηηa=P¯vir/ηηa−Pvirt/ηηaτa⇒Q·virt=P¯vir/ηηa︷Q¯vir−Qvirtτa=Q¯vir−Qvirtτa,Qa=modelMinQmax,Qvirt,b,τa=30s,Qvirt0=Qb 
The authors declare that they have no competing interests.
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.
References
Cabrera ME,Saidel GM,Kalhan SC,Lactate metabolism during exercise: analysis by an integrative systems modelAm J PhysiolYear: 277R1522R1536  
Stirling JR,Zakynthinaki MS,Saltin B,A model of oxygen uptake kinetics in response to exercise: Including a means of calculating oxygen demand/deficit/debtBull Math BiolYear: 200567989101510.1016/j.bulm.2004.12.00515998492  
Murias JM,Spencer MD,Kowalchuk JM,Paterson DH,Influence of phase I duration on phase II VO2 kinetics parameters estimates in older and young adultsAm J Physiol Regul Integr Comp PhysiolYear: 2011301R218R22410.1152/ajpregu.00060.201121490368  
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 humansJ Sports Med Phys FitnessYear: 20004018519411125760  
Gaesser GA,Poole DC,The slow component of oxygen uptake in humansExerc Sport Sci RevYear: 19962435708744246  
Brooks GA,Lactate production under fully aerobic conditions: the lactate shuttle during rest and exerciseFed ProcYear: 198545292429293536591  
Brooks GA,Current concepts in lactate exchangeMed Sci Sports ExercYear: 1991238959061956262  
Brooks GA,Lactate shuttles in natureBiochem Soc TransYear: 20023025826412023861  
Brooks GA,Link between glycolytic and oxidative metabolismSports MedYear: 20073734134310.2165/000072562007370400001717465603  
Moxnes JF,Hausken K,A mathematical model for the training impulse and lactate influx and outflux during exerciseInt J Mod Phys CYear: 200920114717710.1142/S0129183109013522  
Carlson FD,Siger A,The mechanochemistry of muscular contraction I. The isimetric twitchJ Gen PhysiolYear: 196044336010.1085/jgp.44.1.3313690828  
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 humansJ PhysiolYear: 199951892193210.1111/j.14697793.1999.0921p.x10421675  
Margaria R,Biomechanics and energetic of muscular exerciseYear: 1976Oxford: Oxford University Press  
Morton RH,A three component model of human bioenergeticsJ Math BiolYear: 19862445146610.1007/BF012368923805905  
Crowther GJ,Kemper WF,Carey MF,Conley KE,Control of glycolysis in contracting skeletal muscle. II. Turning it offAm J Physiol Endocrinol MetabYear: 20022827479  
Jubrias SA,Esselman PC,Price LB,Cress ME,Conley KE,Large energetic adaptations of elderly muscle to resistance and endurance trainingJ Appl PhysiolYear: 2001901663167010.1063/1.137445211299253  
Lanza IR,Wigmore DM,Befroy DE,KentBraun JA,In vivo ATP production during freeflow and ischaemic muscle contractions in humansJ PhysiolYear: 200657735336710.1113/jphysiol.2006.11424916945975  
Medbø JI,Mohn AC,Tabata I,Bahr R,Vaage O,Sejersted OM,Anaerobic capacity determined by maximum accumulated O2 deficitJ Appl PhysiolYear: 19886450603356666  
Bangsbo J,Quantification of anaerobic energy production during intense exerciseMed Sci Sports ExercYear: 19983047529475643  
Gastin PB,Energy system interaction and relative contribution during maximal exerciseSports MedYear: 20013172574110.2165/000072562001311000000311547894  
Noordhof DA,de Kroning JJ,Foster C,The maximal accumulated oxygen deficit method: a valid and reliable measure of anaerobic capacity?Sports MedYear: 20104028530210.2165/115303900000000000000020364874  
Moxnes JF,Sandbakk K,The kinetics of lactate production and removal during wholebody exerciseTheor Biol Med ModelYear: 20129710.1186/174246829722413898  
Jeneson JAL,Schmitz JPJ,van den Broek NMA,van Riel AAW,Hilbers PAJ,Nicolay K,Prompers JJ,Magnitude and control of mitocondrial sensitivity to ADPAm J Physiol Endocrinol MetabolYear: 2009297E774E78410.1152/ajpendo.00370.2009  
Mahler M,First order kinetics of muscle oxygen consumption, and an equivalent proportionality between Qo2 and Phosphorylcreatine levelJ Gen PhysiolYear: 19858613516510.1085/jgp.86.1.1354031824  
Meyer RA,A linear model of muscle respiration explains monoexponential phosphocreatine changesJ Appl PhysiolYear: 198875648656  
GonzalezAlonso J,Quistorff B,Krustrup P,Bangsbo J,Saltin B,Heat production in human skeletal muscle at the onset of intense dynamic exerciseJ PhysiolYear: 2000524260361510.1111/j.14697793.2000.00603.x10766936  
di Prampero PE,Ferretti G,The energetics of anaerobic muscle metabolism: a reappraisal of older and recent conceptsRespir PhysiolYear: 19991182–310311510647856  
di Prampero PE,Factors limiting maximal performance in humansEur J PhysiolYear: 20039042042910.1007/s004210030926z  
Losnegard T,Myklebust H,Hallen J,Energy system contribution as a determinant of performance in elite skiersYear: 201110.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 exerciseAm J Physiol Endocrinol MetabolYear: 2001280E956E964  
Bertuzzi RC,Franchini E,Ugrinowitsch C,Kokubun E,LimaSilva AE,Pires FO,Nakamura FY,Kiss MA,Predicting MAOD using only a supramaximal exhaustive testInt J Sports MedYear: 201031747748110.1055/s0030125337520432195  
Cerretelli P,Pendergast DR,Paganelli WC,Rennie DW,Effects of specific muscle training on VO2on response and early blood lactateJ Appl PhysiolYear: 197947761769511683  
Binzoni T,Ferretti G,Sechenker K,Cerretelli P,Phosphocreatine hydrolysis in 31PNMR at the onset of constantload exercise in humansJ Appl PhysiolYear: 199273164416491447115  
Hill VH,Lupton H,Muscular exercise, lactic acid, and the supply and utilization of oxygenQ J MedYear: 19231613517110.1093/qjmed/os16.62.135  
Hill AV,Long CNH,Lupton H,Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts IIIIProc Roy SocYear: 19249643847510.1098/rspb.1924.0037  
Hill AV,Long CNH,Lupton H,Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts IVVIProc Roy SocYear: 1924978413810.1098/rspb.1924.0045  
Hill AV,Long CNH,Lupton H,Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts VIIVIIIProc Roy SocYear: 19249715517610.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 contractionAm J PhysiolYear: 1933106689715  
Gaesser GA,Brooks GA,Metabolic bases of excess postexercise oxygen consumption: a reviewMed Sci Sports ExercYear: 19841629436369064  
Laursen PB,Jenkins DG,The scientific basis for highintensity interval training: Optimising training programmes and maximizing performance in highly trained endurance athletesSports MedYear: 200232537310.2165/000072562002320100000311772161 
Figures
[Figure ID: F1] 
Figure 1
The steady state aerobic powerQ¯aand the virtual steady stateQ˜¯in J/s as functions of work rate P in J/s when running with poles at the angle of inclination of 0.105 radians for 5min steady state work rates. The lower horizontal line is the lactate threshold = 1654 J/s, and the upper horizontal line the maximal aerobic power Qmax=1886J/s. ■: Experimental data, v1=2.08m/s,Q¯vir=0.59Qmax, v2=2.33m/s,Q¯vir=0.66Qmax, v3=2.64m/s,Q¯vir=0.74Qmax, v4=2.83m/s,Q¯vir=0.80Qmax, v5=3.08m/s,Q¯vir=0.87Qmax, v6=3.19m/s,Q¯vir=0.90Qmax. 
[Figure ID: F2] 
Figure 2
Steady state lactate concentration C in mmol/L as a function of the fraction of the maximum aerobic power Q_{max}when running with poles at the angle of inclination of 0.105 for 5min steady state work rates.v1=2.08m/s,Q¯vir=0.59Qmax, v2=2.33m/s,Q¯vir=0.66Qmax, v3=2.64m/s,Q¯vir=0.74Qmax, v4=2.83m/s,Q¯vir=0.80Qmax, v5=3.08m/s,Q¯vir=0.87Qmax, C=Cs=12Ln1+p0χd0Qmax1−Qa/Qmax1−p0χd0Qmax1−Qa/Qmax/χ, χ=2.54m3/kg,p0=10−5kg/m3s/J/s,d0=1.110−7/J/s2/s,Qmax=1886J/s,QLT/Qmax=0.88. 
[Figure ID: F3] 
Figure 3
The theoretical and experimental blood lactate concentration C in mmol/L as a function of time t in seconds when running with poles at the angle of inclination of 0.105 radians. The eight symbols from below are experimental values from the different tests v_{18}. The straight curves represent results from a onecompartment model by Moxnes and Hausken [^{2}], whereas the three stapled curves represent results from the two compartment model by Moxnes and Sandbakk [^{25}] for v_{5}, v_{6} and v_{8}. Lines from below: black =v1=2.08m/s,Q¯vir=0.59Qmax, orange =v2=2.33m/s,Q¯vir=0.66Qmax, cyan =v3=2.64m/s,Q¯vir=0.74Qmax, red =v4=2.83m/s,Q¯vir=0.80Qmax, blue =v5=3.08m/s,Q¯vir=0.87Qmax until a recovery phase starts at t = 800s with v = 1.67m/s, Q¯vir=0.47Qmax, green =v6=3.19m/s,Q¯vir=0.90Qmax, yellow =v7=3.33m/s,Q¯vir=0.94Qmax, pink =v8=3.88m/s,Q¯vir=1.09Qmax.LT=0.88Qmax. 
[Figure ID: F4] 
Figure 4
Powers Q in J/s as function of time t in seconds when running with poles at the angle of inclination of 0.105 radians.v = 3.08 m/s, Q¯vir=0.87Qmax when t < 800 s. v = 1.67 m/s, Q¯vir=0.47Qmax when t > 800 s. ____ ___: Aerobic power Qa=Qb+Qaw, ________: Lactic anaerobic power (Q_{G}), : Alactic anaerobic power (Q_{CK}) for Model 1, .........: Alactic anaerobic power (Q_{CK}) for Model 2. 
[Figure ID: F5] 
Figure 5
The lactic and alactic anerobic energy utilization E(t)/E(0) as a function of time t in seconds for different work rates when running with poles at the angle of inclination of 0.105. ______: Lactic energy utilization EGt/EG0,EG0=870J/kg×m,     : Alactic energy utilization ECKt/ECK0,ECK0=830J/kg×m using Model 1, ......: Alactic energy utilization ECKt/ECK0,ECK0=830J/kg×m using Model 2. From below in different colors: pink =v8=3.88m/s,Q¯vir=1.09Qmaxyellow =v7=3.33m/s,Q¯vir=0.94Qmax, green =v6=3.19m/s,Q¯vir=0.90Qmax, blue =v5=3.08m/s,Q¯vir=0.87Qmax, black =v5=3.08m/s,Q¯vir=0.87Qmax until a recovery phase starts at t = 800s with v = 1.67m/s, Q¯vir=0.47Qmax, red =v4=2.83m/s,Q¯vir=0.80Qmax, The exercises are terminated for v_{8}, v_{7} and v_{6}. 
[Figure ID: F6] 
Figure 6
The metabolic power as a fraction of the maximum aerobic power Q/Qmax as a function of time t in seconds for an interval and a continuous exercise. The horizontal lines represent the lactate threshold. : Continuous. _____: Interval. 
[Figure ID: F7] 
Figure 7
The lactate concentration C in mmol/L as a function of time t in seconds for interval and continuous exercise during the protocol showed in Figure6. ....: Continuous. : Interval. 
[Figure ID: F8] 
Figure 8
The anaerobic lactic and alactic powers (Q) in J/s as a function of time t in seconds for interval (black lines) and continuous (red lines) exercise during the protocol showed in Figure6. ______: Lactic anaerobic power (Q_{G}). : Alactic anaerobic power (Q_{CK}), Model 1. ......: Alactic anaerobic power (Q_{CK}), Model 2. 
[Figure ID: F9] 
Figure 9
Calculations of the energy utilization E(t)/E(0) as a function of time t in seconds for (black lines) and continuous (red lines) exercises during the protocol showed in Figure6. _______: Lactic anaerobic energy utilization EGt/EG0,EG(0)=870J/kg×m.      : Alactic anaerobic energy utilization ECKt/ECK0,ECK0=830J/kg×m Model 1. ......: Alactic anaerobic energy utilization ECKt/ECK0,ECK0=830J/kg×m Model 2. 
Article Categories:
Keywords: Blood lactate, Crosscountry skier, Mathematical model, Oxygen uptake, Phosphocreatine. 
Previous Document: SolventDependent Gating Motions of an Extremophilic Lipase from Pseudomonas aeruginosa.
Next Document: Breaking Bad News in a Southeast European Population: A Survey among Cancer Patients in Albania.