The present work proposed to study the oxygen uptake slow component (VO2 SC) of breaststroke swimmers at four different intensities of submaximal exercise, via mathematical modeling of a multi-exponential function. The slow component (SC) was also assessed with two different fixed interval methods and the three methods were compared. Twelve male swimmers performed a test comprising four submaximal 300 m bouts at different intensities where all expired gases were collected breath by breath. Multi-exponential modeling showed values above 450 ml·min−1 of the SC in the two last bouts of exercise (those with intensities above the lactate threshold). A significant effect of the method that was used to calculate the VO2 SC was revealed. Higher mean values were observed when using mathematical modeling compared with the fixed interval 3rd min method (F=7.111; p=0.012; η2=0.587); furthermore, differences were detected among the two fixed interval methods. No significant relationship was found between the SC determined by any method and the blood lactate measured at each of the four exercise intensities. In addition, no significant association between the SC and peak oxygen uptake was found. It was concluded that in trained breaststroke swimmers, the presence of the VO2 SC may be observed at intensities above that corresponding to the 3.5 mM-1 threshold. Moreover, mathematical modeling of the oxygen uptake on-kinetics tended to show a higher slow component as compared to fixed interval methods.
During the step transition from rest to constant-load exercise, oxygen uptake (VO2) has shown varied responses to different workloads, being directly dependent on the type of exercise (Billat et al., 1998) and the level of endurance (Carter et al., 2000; Casaburi et al., 1987). In 1982, Whipp et al. (1982) theorized a model of VO2 onkinetics describing three phases. The first one, the cardio-dynamic component, corresponds to a fast increase in alveolar O2, enabling a 15-20 s transient plateau after on-transition (Demarle et al., 2005). The second phase, called the primary component, is characterized by an exponential increase of VO2 due to muscle oxygen demand, lasting about 2-3 min in healthy subjects (Barstow and Mole, 1991). The third phase is described by an additional slow rise in VO2, superimposing the primary component initiated at exercise onset. This slow phase, commonly called the slow component (SC), continues until either a delayed steady state is attained or maximal VO2 is reached (Whipp, 1994). The magnitude of the SC is considered to have a physiological meaning only when it is higher than 200 ml·min−1 (Billat et al., 2000).
The type of muscle contraction (Pringle et al., 2002), the pattern of motor unit recruitment (Carter et al., 2000) and total muscle mass involved in exercise (Ryschon et al., 1997) have direct influence on the SC. Indeed, the SC has been reported to be sensitive to different types of exercise, being larger for running when compared to cycling exercise (Billat et al., 1998). The SC, which was first observed by Astrand and Saltin (1961), has been assessed through a number of methods including: the difference between VO2 at the 6th and the 3rd min of exercise (Koppo and Bouckaert, 2001), the difference between VO2 at the last and the 3rd min of exercise (Patterson and Whipp, 1991), the difference between VO2 at the last and the 2nd min of exercise (Koppo and Bouckaert, 2002), a single exponential function (Wasserman and Whipp, 1983), a linear function superimposing the exponential component (Casaburi et al., 1987) and the amplitude of a multi exponential function (Santos and Gianella-Neto, 2006). Green and Dawson (1996) observed that the SC was more pronounced between the 3rd and 6th min of exercise, which may explain the popularity of this method among researchers. Notwithstanding, Bearden and Moffat (2001) suggested that this method tended to underestimate the SC when compared with mathematical modeling.
Studies about the on-kinetics of VO2 during swimming are scarce. The SC has already been reported in front crawl swimming (Demarie et al., 2001; Fernandes et al., 2003, 2008; Machado et al., 2006; Reis et al., 2010), but only one study has attempted to evaluate the SC during breaststroke swimming (Reis et al., 2010). As they are, alternated and simultaneously swimming techniques, respectively, their kinematics (Thompson et al., 2004) and energetics (Barbosa et al., 2008) differ. Reis et al. (2010) concluded, by observing the difference between VO2 at the end and at the 3rd min of exercise, that a SC seemed to be present only at swimming intensities above that corresponding to the lactic threshold, as occurs in land-based activities (Billat et al., 1998; Machado et al., 2006). However, to date, no study has modeled with a multi exponential function the kinetics of VO2 in breaststroke swimming.
Therefore, the aim of the present study was to analyze the VO2 slow component of trained breaststroke swimmers during submaximal incremental exercise, using a multi exponential function. It was hypothesized that the SC would only appear at swimming intensities higher than the 3.5 mM-1 lactate threshold. In addition, a comparison of the multi exponential function with the method of fixed intervals was carried out. Potential relationships between the SC and blood lactate concentration (LA-) as well as peak VO2 were also investigated.
Material and Methods
Twelve male high-level swimmers volunteered to participate in the study. Mean (±SD) values for age, body mass and height, the arm spam as well as the estimated body fat percentage were 16.9 ±2.8 yrs, 69.7 ±9.0 kg, 178 ±5.0 cm, 186 ±6.0 cm and 11.6 ±2.3%, respectively. All subjects were involved in systematic training (from 9 to 12 training sessions a week) and participated regularly in swimming breaststroke events. The swimmers’ mean personal best performance at the time of the measurements in the 100 m and 200 m breaststroke short course were 70.20 ±5.02 and 153.67 ±10.40 s, respectively. All participants signed a written informed consent form and those under 18 years of age had their parent’s written consent to participate in the study. All ethical standards for research with human subjects were in accordance with the Declaration of Helsinki and were approved by the Review Board of the University of Trás-os-Montes & Alto Douro.
Measurement of pulmonary gas exchange
Pulmonary gas exchange variables, particularly VO2, were calculated from the gas concentration and flow using a telemetric metabolic cart (K4b2, Cosmed, Italy). Subjects breathed through a low-resistance attached snorkel valve (AquaTrainer, Cosmed, Italy). Before each test, a reference air calibration of the device was performed using a gas sample with 16% O2 and 5% CO2 concentrations, according to the manufactured instructions, as described by Reis et al. (2010). The data was telemetrically displayed and analyzed breath-by-breath.
VO2 on-kinetics model
Prior to VO2 kinetic modeling, a linear interpolation was applied to all breath-by-breath collected data to fit the time response to 1 s intervals. Thereon, an average filter with 11 samples (data ± 5 values) was used to smooth the data. According to the short-term fluctuations observed in the data, the number of samples used before and after, the central value was chosen to eliminate the high-frequency oscillations and simultaneously to preserve the information from the original data.
The kinetics of VO2 was modeled by the following exponential function:
A non-linear regression analysis was applied to fit the time responses of VO2, in which coefficients were obtained by the least squares method, being the oxygen uptake slow component determined by mathematical modeling using the magnitude of the second exponential amplitude (MMSC).
The calculation of the VO2 SC through the fixed interval methods consisted of subtracting average VO2 observed in the last min of each bout by:
average VO2 observed in the 3rd min of exercise (∆SC3min);
average VO2 observed in the 2nd min of exercise (∆SC2min).
The second exponential rise was considered to represent a SC when it surpassed 200 ml·min−1 (Billat et al., 2000).
All mathematical procedures and modeling were done using MATLAB R2010b (Mathworks, USA) for Windows®.
Measurement of blood lactate
Immediately after the completion of each bout, capillary blood samples were collected from the fingertip to determine LA-values (Accusport, Boehringer, Germany). Before each test a calibration of the device was performed with several YSI 1530 Standard Lactate Solutions (2, 4, 8 and 16 mM−1). The lactate threshold was accepted to be the value corresponding to 3.5 mM1 of LA- (Fernandes et al., 2003).
Experimental testing was conducted in an indoor 25 m swimming pool, with water temperature of 28.5o C and air humidity of 55%. During the five days prior to the experiment subjects were not involved in high-intensity training and limited their training program to a single daily low-intensity swimming session. The experimental procedures were conducted in the morning between 10 and 12 am.
All subjects performed a graded individualized swimming test, with four 300 m exercise bouts, each one performed at a controlled constant speed (ensuring that all subjects performed, at least, 5 min of exercise in each bout); a recovery period of 15 min was introduced between successive bouts. Initial speed corresponded to 50% of the swimmer´s estimated best performance at the 300 m breaststroke event and subsequent increments were equal to 10% of the 300 m best performance speed, controlled by an electronic pacer (TAR. 1.1, GBK-electronics, Aveiro, Portugal), with successive flashing lights placed at the bottom of the pool. In-water starts and open turns were used, and only a dry land warm-up was performed before the test. The performance in the 300 m breaststroke was estimated taking into account the performance at the 100 m and 200 m breaststroke (Reis et al., 2010):
where v100 and v200 represent the velocity of the 100 m and 200 m breaststroke event and v300 and v400 represent the estimation of the velocity of the 300 m and 400 m breaststroke event.
Normal distribution of data was confirmed by the Shapiro-Wilk test and sphericity by the Mauchly’s test. To compare the SC according to the different intensities and calculation methods, a two factor repeated measure analysis of variance (ANOVA) with a post-hoc Bonferroni test was performed. With the aim of identifying potential relationships between the SC, LA– and peak oxygen uptake (VO2peak), Pearson or Spearman correlations were applied, whenever appropriate. The level of significance was set at 5% (p ≤ 0.05). Results are presented as means and standard deviations. All data analysis was performed using Statistical Package for Social Sciences (SPSS 17.0 - Science, Chicago, USA) for Windows®.
In Figure 1, the first bout of a representative subject is shown where, as expected, no slow component was observed in VO2 on-kinetics.
In Figure 2, VO2 kinetics of another swimmer is shown during the fourth bout of the experimental protocol, being evident the appearance of the slow component.
Blood lactate in the four bouts of exercise was 2.33 ±0.53, 2.28 ±0.52, 4.03 ±2.45 and 5.72 ±1.73 mM−1. Mean peak VO2 value obtained for the total sample during the incremental protocol was 4.4 ±0.90 l·min−1. Table 1 presents the variables extracted from mathematical modeling of the VO2 on-kinetics and Table 2 compares the magnitude of the slow component reached by the three methods. Table 3 presents the descriptive data of absolute VO2, relative VO2 and the number of subjects that surpassed the 200 ml·min−1 as to the slow component in the four bouts of exercise.
Mean values (±SD) of the parameters extracted from the bi-exponential modeling of oxygen uptake.
|1st stage||2nd stage||3rd stage||4th stage|
|A0, ml·min−1||278.71 ±67.21||270.56 ±59.39||270.79 ±57.74||269.48 ±45.74|
|A1, ml·min−1||708.69 ±285.96||1035.71 ± 386.18||1571.68 ±739.26||1988.16 ±565.06|
|A2, ml·min−1||113.08 ±162.51||227.00 ±269.93||497.76 ±276.07||473.36 ±404.22|
|tau1, s||48.07 ±31.69||60.83 ±27.23||71.09 ±22.36||72.77 ±19.86|
|tau2, s||13.03 ±16.14||14.30 ±13.45||17.83 ±10.04||22.39 ±13.73|
|TD, s||94.82 ±109.14||131.45 ±90.61||147.40 ±34.42||148.78±56.17|
Mean values (±SD) of the slow component calculated by the three methods in each of the four bouts (with indication of mean blood lactate values at each stage)
|1st stage 2.33 ±0.53||2nd stage 2.28 ±0.52||3rd stage 4.03 ±2.45||4th stage 5.72 ±1.73|
|MMSC, ml·min−1||113.08 ±162.51||227.00 ±269.93||497.76 ±276.07||473.36 ±404.22|
|∆SC3min, ml·min−1||107.9 ±54.9||88.3 ±73.6||311.0 ±153.4||325.1 ±160.5|
|∆SC2min, ml·min−1||121.6 ±66.2||109.1 ±90.8||373.2 ±185.7||384.8 ±197.9|
Mean values (±SD) of absolute VO2, relative VO2 and the number of subjects with an evident slow component in each of the four bouts (with indication of mean blood lactate values at each stage)
|1st stage 2.33 ±0.53||2nd stage 2.28 ±0.52||3rd stage 4.03 ±2.45||4th stage 5.72 ±1.73|
|VO2, ml·min−1||1958.05 ±191.74||1970.38 ±251.79||2547.74 ±383.90||4167.85 ±428.67|
|VO2, ml·min−1·kg−1||27.91 ±2.73||28.09 ±3.89||36.32 ±5.47||59.41 ±6.11|
|n > 200 ml·min−1||3||4||9||8|
The repeated measures analysis of variance revealed a significant effect of the method that was used to calculate the VO2 SC. Higher mean values were observed when using mathematical modeling compared with the ∆SC3min method (F=7.111; p=0.012; η2=0.587); furthermore, differences were detected among the two interval methods. The effect of exercise intensity was also evident (F=6.680; p=0.011; η2=0.690), and significant differences were observed between the second and the third bout, and between the second and the fourth bout of the exercise. Both the method and intensity effects presented a strong statistical power (0.830 and 0.865, respectively). However, there was no significant effect of the interaction method vs. intensity on the VO2 SC (F=2.014; p=0.012; η2=0.668).
No significant relationship was found between the SC determined by whichever method and blood lactate measured at each of the four exercise intensities. In addition, no significant relationship between the SC and peak VO2 was found.
The aim of the present study was to analyze the SC of trained breaststroke swimmers during submaximal exercise using a multi exponential function. Testing procedures were conducted in ecological swimming conditions, using a modern telemetric procedure for collecting and measuring the breath-by-breath expired gases, and a suitable modified snorkel and valve system (Keskinen et al., 2003) specific for breath-by-breath data analysis. It was hypothesized that the SC would only appear during swimming at intensities above 3.5 mM-1 of the LA- threshold, i.e., at the heavy intensity domain.
The main finding of the present study was that the SC determined by mathematical modeling was evident in most subjects at intensities higher than that corresponding to 3.5 mM-1. The results of both the interval methods indicated the presence of the SC in the last two bouts of exercise, both with exercise intensity above the 3.5 mM-1 lactate threshold.
The decision to invite highly-trained individuals to participate in the present study relied on the fact that top level swimmers were expected to be characterized by maximal physical fitness, very close to human genetic limits (Fernandes et al., 2008). Another very important reason was the fact that athletes were the main target population in terms of the possible applicability of this type of research.
The protocol was designed to investigate VO2 kinetics in different intensity domains. The SC is commonly described in studies where the intensity corresponds to the 3.5 mM-1 lactate threshold (Fernandes et al., 2003, 2008; Reis et al., 2010, 2012). In the protocol of Demarie et al. (2001) constant-velocity front crawl swimming in a flume above critical velocity was used. Fernandes et al. (2003) used a test protocol where subjects swam at a pace corresponding to their VO2max velocity until volitional exhaustion. This was repeated in 2006 (Machado et al., 2006) and 2008 (Fernandes et al., 2008) in tests where swimmers swam until exhaustion at the previously determined VO2max velocity. However, to the best of our knowledge, the SC has not been investigated yet using mathematical modeling of the moderate-domain in breaststroke swimming.
In the present study, the SC was assessed through mathematical modeling (Santos and Gianella-Neto, 2006), a more precise and accurate method (Barstow and Mole, 1991) than rather simplistic calculation of the increase in VO2 between the second or third minute and the last minute of exercise (Bearden and Moffatt, 2001; Machado et al., 2006). Moreover, the mathematical model used in the present study allows the discrimination of different components of VO2 kinetics. The obtained results are markedly different from the findings of all other studies published about the SC in swimming. As an example we may indicate the single work found in literature with the same methodological characteristics (Fernandes et al., 2008), however, performed with front crawl swimmers, where the SC values were positive in all subjects and the mean of the sample was equal to 356.27 ±168.16 ml·min−1.
To our knowledge, the study by Reis et al. (2010) is the only report in the literature with data for the VO2SC in breaststroke swimming. However, their SC was calculated by a fixed interval method (difference between mean VO2 at the last and third minute of exercise). A SC above 350 ml·min−1 was verified at all four exercise intensities of testing, which comprehended different domains of intensity. In fact, according to the literature, we expected to observe the SC in the last two stages of effort, stages related with the heavy and maximal domain of intensity (Denadai et al., 2000). On the other hand, the moderate domain, where the first two bouts of the present study are included, is typically reported in the literature as of an exercise intensity which does not enable the presence of the SC (Gaesser and Poole, 1996). However, our results do not confirm such an assumption for breaststroke swimming, as the 2nd bout of exercise showed a mean value higher than 200 ml·min−1, a limit often referred as being the lower threshold for an evident SC (Billat et al., 2000).
The absence of a relationships between post-exercise LA- and the SC, irrespective of the method used for its assessment, is in agreement with some authors (Mole and Hoffmann, 1999; Ogita and Tabata, 1992) who stated it was unlikely that blood lactate per se could be responsible for the SC phenomenon, but rather the concurrent acidosis. Hence, our SC observations could alternatively be explained by larger rates of type II fiber recruitment and additional energy cost of contraction (Whipp, 1994) at higher exercise intensities. However, Casaburi et al. (1987) and Billat et al. (1998) found a significant correlation between the SC and LAobtained after high intensity exercise.
In summary, our results indicated that mathematical modeling of VO2 on-kinetics during submaximal breaststroke swimming revealed, in most cases, a higher SC as compared to fixed interval methods. The SC values observed with the aid of mathematical modeling confirmed the hypothesis of a presence of the SC at intensities above the 3.5 mM-1 lactate threshold.
Astrand PO, Saltin B. Maximal oxygen uptake and heart rate in various types of muscular activity. J Appl Physiol, 1961; 16: 977-81
Barbosa TM, Fernandes RJ, Keskinen KL, Vilas-Boas JP. The influence of stroke mechanics into energy cost of elite swimmers. Eur J Appl Physiol, 2008; 103(2): 139-49
Barstow TJ, Mole PA. Linear and nonlinear characteristics of oxygen uptake kinetics during heavy exercise. J Appl Physiol, 1991; 71(6): 2099-106
Bearden SE, Moffatt RJ. VO2 slow component: to model or not to model? Med Sci Sports Exerc, 2001; 33(4): 677-80
Billat VL. VO2 slow component and performance in endurance sports. Br J Sports Med, 2000; 34(2): 83-85
Billat VL, Richard R, Binsse VM, Koralsztein JP, Haouzi P. The VO2 slow component for severe exercise depends on type of exercise and is not correlated with time to fatigue. J Appl Physiol, 1998; 85(6): 2118-24
Carter H, Jones AM, Barstow TJ, Burnley M, Williams CA, Doust JH. Oxygen uptake kinetics in treadmill running and cycle ergometry: a comparison. J Appl Physiol, 2000; 89(3): 899-907
Casaburi R, Storer TW, Ben-Dov I, Wasserman K. Effect of endurance training on possible determinants of VO2 during heavy exercise. J Appl Physiol, 1987; 62(1): 199-207
Demarie S, Sardella F, Billat V, Magini W, Faina M. The VO2 slow component in swimming. Eur J Appl Physiol, 2001; 84(1-2): 95-9
Demarle AP, Slawinski JJ, Laffite LP, Bocquet VG, Koralsztein JP, Billat VL. Decrease of O2 deficit is a potential factor in increased time to exhaustion after specific endurance training. J Appl Physiol, 2001; 90(3): 947-53
Denadai BS, Greco CC, Teixeira M. Blood lactate response and critical speed in swimmers aged 10-12 years of different standards. J Sports Sci, 2000; 18(10): 779-84
Fernandes RJ, Keskinen KL, Colaco P, Querido AJ, Machado LJ, Morais PA, Marinho DA, Vilas-Boas JP. Time limit at VO2max velocity in elite crawl swimmers. Int J Sports Med, 2008; 29(2): 145-50
Fernandes RJ, Cardoso CS, Soares SM, Ascensao A, Colaco PJ, Vilas-Boas JP. Time limit and VO2 slow component at intensities corresponding to VO2max in swimmers. Int J Sports Med, 2003; 24(8): 576-81
Gaesser GA, Poole DC. The slow component of oxygen uptake kinetics in humans. Exerc Sport Sci Rev, 1996; 24: 35-71
Green S, Dawson BT. Methodological effects on the VO2-power regression and the accumulated O2 deficit. Med Sci Sports Exerc, 1996; 28(3): 392-7
Keskinen KL, Rodriguez FA, Keskinen OP. Respiratory snorkel and valve system for breath-by-breath gas analysis in swimming. Scand J Med Sci Sports, 2003; 13(5): 322-9
Koppo K, Bouckaert J. The effect of prior high-intensity cycling exercise on the VO2 kinetics during high-intensity cycling exercise is situated at the additional slow component. Int J Sports Med, 2001; 22(1): 21-6
Koppo K, Bouckaert J. The decrease in VO2 slow component induced by prior exercise does not affect the time to exhaustion. Int J Sports Med, 2002; 23(4): 262-7
Machado FA, Guglielmo L, Greco CC, Denadai BS. Slow component of O2 in children during running exercise performed at heavy intensity domain: analysis with different mathematical models. Rev Bras Med Esporte, 2006; 12(6): 308-13
Machado L, Querido AJ, Keskinen KL, Fernandes RJ, Vilas Boas JP. Mathematical modelling of the slow component of oxygen kinetics in front crawl swimming. Rev Port Cien Desp, 2006; 6: 142-4
Mole PA, Hoffmann JJ. VO2 kinetics of mild exercise are altered by RER. J Appl Physiol, 1999; 87(6): 2097-106
Ogita F, Tabata I. Oxygen uptake during swimming in a hypobaric hypoxic environment. Eur J Appl Physiol Occup Physiol, 1992; 65(2): 192-6
Patterson D, Whipp BJ. Asymmetries of oxygen uptake transients at the on- and offset of heavy exercise in humans. J Physiol, 1991; 443: 575-86
Pringle JS, Carter H, Doust JH, Jones AM. Oxygen uptake kinetics during horizontal and uphill treadmill running in humans. Eur J Appl Physiol, 2002; 88(1-2): 163-9
Reis JF, Alves FB, Bruno PM, Vleck V, Millet GP. Effects of aerobic fitness on oxygen uptake kinetics in heavy intensity swimming. Eur J Appl Physiol, 2012; 112(5): 1689-97
Reis VM, Marinho DA, Barbosa FP, Reis AM, Guidetti L, Silva AJ. Examining the accumulated oxygen deficit method in breaststroke swimming. Eur J Appl Physiol, 2010; 109(6): 1129-35
Reis VM, Marinho DA, Policarpo FB, Carneiro AL, Baldari C, Silva AJ. Examining the accumulated oxygen deficit method in front crawl swimming. Int J Sports Med, 2010; 31(6): 421-7
Ryschon TW, Fowler MD, Wysong RE, Anthony A, Balaban RS. Efficiency of human skeletal muscle in vivo: comparison of isometric, concentric, and eccentric muscle action. J Appl Physiol, 1997; 83(3): 867-74
Santos EL, Gianella-Neto A. Estimation of the VO2 peak, ventilatory threshold and the respiratory compensation point based on the gas exchange kinetics during the transition to constant-load exercise. Rev Bras Eng Bio, 2006; 22(2): 107-18
Thompson KG, MacLaren DP, Lees A, Atkinson G. The effects of changing pace on metabolism and stroke characteristics during high-speed breaststroke swimming. J Sports Sci, 2004; 22(2): 149-57
Wasserman DH, Whipp BJ. Coupling of ventilation to pulmonary gas exchange during nonsteady-state work in men. J Appl Physiol, 1983; 54(2): 587-93
Whipp BJ. The slow component of O2 uptake kinetics during heavy exercise. Med Sci Sports Exerc, 1994; 26(11): 1319-26
Whipp BJ, Ward SA, Lamarra N, Davis JA, Wasserman K. Parameters of ventilatory and gas exchange dynamics during exercise. J Appl Physiol, 1982; 52(6): 1506-13