1

1 Medium term effects of physical conditioning on breath-hold diving performance

2

3 1. Introduction

4 In indoor diving competition, breath-hold divers compete to remain immersed

5 as long as possible in (STA), and to dive the longest distance in the

6 categories dynamic with fins (DYN) and dynamic no fins (DNF). Adverse effects, such

7 as sickness [Batle, 1999], narcosis or arterial embolism [Batle, 2002],

8 observed during deep events, are not going to occur in the pool practice where

9 submerges barely 2 meters; thus, hypoxic syncope is the main risk during apnea indoor

10 practice [Fitz-Clarke, 2006].

11“Recent improvements in equipment, nutrition, and training have led to increases in

12 apnea performance [Schagatay, 2009; Schagatay, 2010; Fernández, 2015]. Practice and

13 training is mainly responsible for performance in breath-hold divers; however, it is

14 difficult to establish a rigid training model that suits all divers with different levels of

15 ability, age, previous physical condition or anthropometry.

16 The effects of physical activity on apnea performance have not been clarified,

17 probably because the effects of physical training and apnea training were studied

18 independently. The effects of physical training on apnea performance are controversial.

19 While Bavagad showed in his study [Bagavad, 2014] a correlation between physical

20 conditioning with apnea performance, Schagatay [Schagatay, 2000] have not found did

21 not find effects on apnea performance; however, this latter research concluded that

22 physical training increased stamina, allowing prolongation of the struggle phase of

23 apnea. In addition, a few longitudinal research studies of apnea training have been

24 conducted [Schagatay, 2000; Bagavad, 2014; Fernández, 2015] .The study

1 2

25 aimed to analyze the effects of physical conditioning inclusion on apnea performance

26 after a 22-week structured apnea training program.

27

28 2. Materials and Methods

29 2.1. Participants

30 Twenty-nine male breath-hold divers (36 ± 5 years of age) with two years of

31 experience in breath-hold diving were divided in different training programs: (I) cross-

32 training in apnea and physical conditioning (CT; n=10); (II) apnea training (AT; n=10);

33 and (III) control group (CG; n=9). No training intervention was carried out on the CG;

34 however, regarding pre-post measurements, both the CG and the experimental groups

35 performed them without any distinction between them. The participants were informed

36 of the benefits and risks prior to signing the informed consent document to participate in

37 the research. Participants diagnosed with cardiac, metabolic or respiratory diseases were

38 excluded. The study was conducted in accordance with the Declaration of Helsinki

39 [Harris, 2011] and approved by the institutional Human Research Ethics Committee

40 (CSEULS-PI-114/2016).

41

42 2.2. Training procedures

43 All participants conducted 66 sessions of 60-minute of duration each. During the first

44 session, CT was performed for 15 minutes with a dry strength circuit of 10 calisthenic

45 exercises - squat, pull-up, push-up, squat, pull-up, push-up, squat, pull-up, push-up and

46 squat- with 50 repetitions and 10 seconds of recovery. This was followed by swimming

47 training in crawl, alternating each week between 45 min of swimming continuous

48 training (14 RPE) and high intensive interval training (HIIT) (20 RPE). Regarding

49 HIIT, volume and density were 3 series x 10 repetitions of 25 meters with 20 seconds of

2 3

50 recovery. During the second session, participants performed a specific hypoxic training

51 in dynamic, that consisted of 40 repetitions of 25 m dynamic to maximum

52 individual underwater swimming speed with 20 s of recovery each 25 m. In the third

53 session, divers conducted a specific hypoxic training in static. In this training divers

54 performed 3 maximal static apneas with 10 min of full recovery between each apnea.

55 The AT group, performed a specific hypoxic training in dynamic, on the first and

56 second session, and a specific hypoxic training in static, on the third session.

57 Glossopharyngeal insufflation (Seccombe, 2006) was not allowed during training for

58 any group.

59

60 2.3. Test procedures

61 Measurements were performed under similar environmental conditions in a sports

62 laboratory, a health center and a 25-meter pool with 2-meter depth. The test battery

63 sequence was performed chronologically in five visits in different days, as follows:

64 Visit 1– STA; Visit 2 – DNF; Visit 3– DYN; Visit 4 – body composition, count,

65 spirometry and an incremental test on treadmill; Visit 5 – resting metabolic rate, and

66 pulse-oximetry during a static apnea in dry conditions.

67

68 2.1.1. STA, DYN and DNF performance

69 Same protocol was performed, previous to apnea, to avoid possible warming effects.

70 Thus, before STA, a 15-minute warm-up consisted of 10 minutes of relaxation, two

71 minutes of static apnea and a three-minute countdown on the surface until maximal

72 attempt. Regarding the DYN and DNF, a 15-minute warm-up was performed

73 comprising 10 minutes of relaxation, 50 meters in the specific (fins or no fins) dynamic

74 discipline and a three-minute countdown on the surface. After warm-up, the divers

3 4

75 attempted to achieve the maximal individual time or distance. In all the disciplines,

76 during the last 30 seconds of the countdown, the participant placed the nose clip and

77 performed a maximal inspiration.

78

79 2.1.2. Physiological variables

80 Body fat and fat-free mass were measured by whole-body dual-energy X-ray

81 absorptiometry (GE Lunar Prodigy; GE Healthcare, Madison, Wisconsin).Regarding

82 hematological values, hemoglobin (Hb) and hematocrit were analyzed by a blood count

83 analysis at the health center and under fasting conditions; a portable spirometry test

84 (Spirostik, Geratherm) was performed as well to calculate vital capacity (VC). Before

85 starting the test, we explained the technique to the subject: sitting, with knees bent and

86 feet resting on the floor. In this way, the nasal clip was placed and the participant would

87 make a relaxed, but maximum inspiration, at the end of which the mouthpiece would be

88 placed; then, he was to perform slow expiration, maintaining a constant flow until

89 residual volume. This procedure was be repeated three times, with at least 30 seconds of

90 rest between each measurement and the highest value of the three measurements was

91 registered [Wanger, 2005].

92 The incremental test to measure maximal consumption(VO2max)and maximal

93 rate was performed on a treadmill (H/P/COSMOS 3P® 4.0, H / P / Cosmos Sports

94 & Medical, Nussdorf-Traunstein, Germany). The volume and composition of expired

95 gases were measured using a gas analyzer (Ultima CPX, Medical Graphics) and the

96 measured by electrocardiogram (WelchAllyn, CardioPerfect). After a five-

97 minute warm-up at 10 km•h-1, the speed was increased 1 km•h-1 every minute until

98 volitional exhaustion of the participant. Throughout the test the treadmill elevation was

99 maintained at 1%. VO2max was noted as the average of the two highest values of 15

4 5

100 seconds of VO2 consecutive reached toward the end of the test, as long as the values of

101VO 2 and HR have stabilized despite the increase in the intensity of the effort.

102 Resting metabolic rate was measured, breath-to-breath (Ultima CPX, Medical

103 Graphics), over 15 minutes by indirect calorimetry. According to a previous study

104 [Melanson, 2002], a minimum of 15 minutes of steady state, determined as <10%

105 fluctuation in oxygen consumption and <5% fluctuation in respiratory exchange ratio

106 between CO2 and O2, was considered criteria for valid resting metabolic rate(RMR).

107 Participants performed the test in fasting conditions after 12 hours without physical

108 effort. The environmental conditions established were as follows: a quiet laboratory

109 with the lights off, 550 meters of altitude, 22.5 ± 1°C, with 55 ± 3% relative humidity.

110 Heart rate and oxygen saturation were monitored during a maximal static apnea in

111 dry conditions. Before beginning the test, the diver rested for 10 minutes in the prone

112 position, with head and upper limbs leaning on a table placed in front of the stretcher.

113 During the last 30 seconds of the countdown, a nose clip was placed to avoid possible

114 air leakage. At that point, the participant performed a maximal followed by a

115 maximal inspiration; glossopharyngeal insufflation was not allowed. Throughout the

116 static apnea test, the average heart rate (HRavg) and the minimum heart rate (HRmin)

117 were recorded. To calculate the HRavg, only the heart rate data after first 30 seconds

118 were analyzed; i.e., once the heart rate was stabilized [Breskovic, 2012]. The oxygen

119 saturation was monitored by a pulse oximeter (CMS 50F) placed on the second finger of

120 the left hand. A minimal value of oxygen saturation was collected during the test

121 (SpO2min). Time during static apnea in dry conditions (STAdry) was recorded, also.

122

123 2.3. Statistical analyses

5 6

124 The effects of training on performance and physiological variables were analyzed

125 using repeated measures analysis of variance (ANOVA). In order to improve the power

126 of the test and to reduce Type I error, a repeated measures multiple analysis of

127 variance(MANOVA) was used in the case of a set of variables defined in a cohesive

128 theme [Tabachnick, 2007]. Moderate correlations(r=0.3 to 0.7) led to determine the

129 group of variables: STA, DYN and DNF in the case of performance variables, whereas

130 the groups VO2max, Hb and VC (related to the storage), percent of body fat (PBF)

131 and lean body mass (LBM) (measuring body composition) and RMR and HRavg

132 (corresponding to the individual metabolic rate) were established with the physiological

133 variables.

134 A post-hoc analysis was performed when significant differences were detected

135 between the groups. Regarding ANOVA, the specific contrast tests were used, when the

136 homogeneity of the variances was achieved, with Bonferroni correction. When

137 interaction occurs, the effect of one factor on a response depends on the level of

138 (an)other factor(s). Thus, we cannot consider the main effects of the factors separately

139 as the main effects and interaction need to be considered as a whole to describe the

140 relationship between input and output. However, discriminant analysis was used as a

141 post-hoc procedure of the MANOVA, with the objective of determining which

142 dependent variables discriminate between the training programs and analyzing their

143 contributions. The significance level considered in all the developed tests was 0.05.

144

145 3. Results

146 Tables 1 and 2 summarize the mean and standard deviations (SD) of the scores

147 obtained before and after the training interventions as well as the differences between

148 them. (Tables 1 and 2)

6 7

149 3.1. Performance variables

150 Total performance, referred as POINTS in this work, was used as a global performance

151 variable on a Apnea Indoor diving. It is constructed from the variables STA, DNF and

152 DYN. In particular, each second the athlete remains immersed in STA was multiplied

153 by 0.2, whereas each meter reached at DNF and DYN was multiplied by 0.5. AIDA

154 International – one of the official competition organizing bodies – uses these

155 scores during championships: The diver with the highest score wins the competition.

156 The analysis of variance of repeated measured data for the POINTS variable indicated

157 the effect of the training group (F = 6.270, p<0.01), time (F = 45.019, p<0.001) and the

158 interaction between both (F = 7.409, p<0.01) to be statistically significant. Profile

159 graphs (Figure 1) show weak discrepancies between CT and AT estimated marginal

160 means. Nevertheless, there are notable differences between control group and both CT

161 and AT estimated marginal means. In Table 1, significant differences between pairs of

162 measures of both time and training group factors were collected (Figure 1).

163 Group linear discriminant analysis identifies a set of linear combinations of

164 dependent variables (LDF) that discriminate groups. After a stepwise regression the

165 variables included in the model were STA2 and DNF2, leaving out DYN2. The number

166 of LDFs is the number of groups minus 1; thus two discriminant functions were

167 necessary to achieve significant discrimination: Function 1 = 0.575·ZSTA2 + 0.539·ZDNF2

168 and Function 2 = -1.126·ZSTA2 + 1.144·ZDNF2. , where ZX is the standardized score of the

169 variable X. Both LDFs were statistically significant (p<0.001 and p<0.01 respectively).

170 The coefficients of the LDFs are a measure of the importance of each variable for

171 distinguishing groups. Thus, the first function these coefficients indicates that STA2

172 was more important in distinguishing, while DNF2 had a larger contribution for the

173 second function. The graph in Figure 2 shows the subjects’ scores on the discriminant

7 8

174 functions and centroids of each group. We can see from the graph that the centroids are

175 positive in the first function for the groups CT and AT and negative for CG. Therefore,

176 the first function classifies as a control the individuals with the lowest score in the

177 variables STA2 and DNF2, whereas individuals with higher scores would be classified

178 in the CT group and, to a lesser extent, in the AT group. In the second function, we have

179 opposite sign coefficients and similar values for the discriminant variables. Because the

180 centroids that are the most distant are those that correspond to the groups CT and AT,

181 those of opposite sign, the function classifies the participants with greater scores in

182 STA2 and smaller scores in DNF2 within the AT group. In contrast, the participants

183 with a lower score in STA2 and a higher score in DNF2 are classified in the CT group

184 (Figure 2).

185 3.2. Physiological variables

186 The post-hoc analysis determined that the CT group was the only group in which the

187 difference of means was significant before and after training for the VC and VO2max

188 variables. We found a significant increase in VC and VO2max after training (see Table 2).

189 Mean of HRmin significantly decrease after training for CT and AT (see Table 2).

190 Analyzing this variable by pairs of training groups, significant differences were found

191 between the control group and the rest of the groups, being the mean of HRmin higher in

192 the CG than in the others (see Table 2).

193

194 4. Discussion

195 There is current interest in the role of physical training on apnea performance.

196 Previous research [Schagatay, 2000] analyzed the effects of physical training and apnea

197 training separately. In the current study, after the discriminant analysis, we concluded

198 that training groups can be distinguished in two areas: isolated apnea training appears to

8 9

199 be the most suitable for STA, whereas cross-training appears to be best for DNF.

200 Regarding DYN performance, although there were no significant differences between

201 the training groups (AT vs. CT), there was a trend (favorable to CT).

202 As hypothesis, the inclusion of physical training could have triggered an

203 improvement in the diver’s underwater swimming economy; however, this variable has

204 not been measured in our study. Another study [Nadel, 1974] proposed that trained

205 swimmers use less oxygen than those not trained to swim at the same speed. Also, a

206 previous review [Schagatay, 2010] concluded that a high proportion of ex-swimmers

207 represent the elite in the apnea dynamic discipline.

208 According to the correlation shown between physiological variables studied, groups

209 were defined into: lung storage, body composition, metabolic rate, and individual

210 tolerance to . The difference in means after training interventions was

211 significant only in the CT group for the VC and VO2max variables.

212 Maximal aerobic capacity is closely related to performance in endurance sports

213 (Bassett, 2000). Expected benefits, after physical training, would be an increased

214 storage of oxygen before apnea, a decreased heart rate during apnea and increased

215 motivation or stamina. According to the results in this study, the increased aerobic

216 capacity did not result in a lower heart rate during apnea, which is consistent with the

217 findings of Schagatay (Schagatay, 2000). On the other hand, the test for assessing

218 maximum aerobic capacity (incremental test on treadmill) was no specific (swimming);

219 thus, muscular fatigue could produce limitation on aerobic and anaerobic performance

220 during the incremental test on treadmill.

221 Physical training resulted in an increased VC; thus, the inclusion of physical activity

222 in the apnea training plan could be considered appropriate for increasing divers’ total

223 O2/CO2 lung storage. Glossopharyngeal insufflation was not allowed during training or

9 10

224 spirometry test; thus, there was no any improvement or adaptation produced by this

225 technique. Whittaker noted in his study [Whittaker, 2007] that after a specific

226 training consisting of thoracic stretching with glossopharyngeal insufflation divers

227 could increase total lung capacity up to 2 liters. However, the proposed training was

228 insufficient to cause long-term adaptations in hematological values. In addition to

229 training, it would be advisable to include hypoventilation training in low lung

230 volume instead of the traditional training in , mentioned in our study. A

231 study conducted in 2010 [Woorons, 2010] showed that hypoventilation training in low

232 lung volume produced both an increase in CO2 (hypercapnic effect) and

233 a decrease in O2 (hypoxic effect) in the blood and muscles.

234 No tendency was found regarding body composition after training programs,

235 suggesting that AT or CT was insufficient to cause significant long-term changes in

236 body composition. Previous study [Sarang, 2006] proposed a decrease greater than 30%

237 in O2consumption after the application of relaxation techniques; however, achieving the

238 that enables breath-hold divers to reduce their metabolic rate during apnea

239 requires many hours of deliberate practice and years of experience.

240 Decreased in HRmin after AT and CT, compared with the CG, is probably due to a

241 pronounced bradycardia induced by increased hypoxia during the static apnea in dry

242 conditions; i.e., breath-hold divers who followed periodized training achieved greater

243 tolerance to asphyxia and consequently an increase in apnea time. On the other hand,

244 there was no trend for the HRmin between experimental groups, suggesting a lack of

245 cardiac adaptation produced by physical exercise proposed. Physical and psychological

246 tolerance to effort have aspects in common during the struggle phase of apnea or

247 physical exertion; thus, a intensive physical training, as performed through HIIT, could

248 trigger an increased tolerance to the asphyxia produced during apnea. According to

10 11

249 results showed in our study, oxygen saturation reached during dry static apnea were far

250 from limits, i.e.; hypoxic syncope (50% SaO2 in amateurs and 30% in trained freedivers

251 [Andersson, 2004]). In this way, an increase in tolerance to asphyxia in our divers could

252 dilute the effect (on Apnea Indoor performance) of other factors such as: total storage of

253 O2/CO2, rate metabolic rate or underwater swimming economy. A

254 limitations of the present study was that O2 consumption and CO2 release rates were

255 not calculated during the breath-hold; moreover, there are no data on the partial

256 of O2 and CO2 in the alveolar air before the apnea and at the end of the apnea. The

257 availability of these data had would allowed would provide a more detailed analysis of

258 the identified correlations.

259

260 5. Conclusions

261 The inclusion of physical activity in apnea training increased VC and VO2max in breath

262 hold divers; in addition, divers who followed a mixed training, physical training and

263 hypoxic training, achieved increased DNF performance.

264

265 References

266 Andersson J. Cardiovascular and respiratory responses to apneas with and without face

267 immersion in exercising humans. J Appl Physiol. 2004; 96:1005-10

268 Bagavad G. Effect of physical training on breath holding time in Indian subjects. Indian

269 J Physiol Pharmacol. 2014; 58: 108-109.

270 Bassett DR. Limiting factors for maximum oxygen uptake and determinants of

271 endurance performance. Med Sci Sports Exerc. 2000; 32: 70-84.

11 12

272 Batle JM. caused by breath-hold diving hunting. A report of

273 about 30 cases. Proceedings of the 13th International Congress on Hyperbaric

274 Medicine, 1999.

275 Batle JM. How to avoid arterial gas embolism in free diving, the attitude a breath-hold

276 diver should take. A study of 35 cases, in a period 1995-2002. Proceedings of

277 the 14th International Congress on , 2002.

278 Breskovic T. Cardiovascular changes during underwater static and dynamic breath-hold

279 dives in trained divers. J Appl Physiol. 2012; 3:673-8.

280 Fernández FA. Periodization of Apnea Training. Doctoral thesis, University of Castilla-

281 la Mancha, Toledo, 2015.

282 Fitz-Clarke, J. Adverse events in competitive breath-hold diving. Undersea Hyperb

283 Med. 2006 Jan-Feb; 33(1): 55-62.

284 Harris D, Atkinson G. Update - Ethical standards in sport and exercise science research.

285 Int J Sports Med. 2011; 32: 819-21.

286 Melanson EL et al. Resistance and aerobic exercise have similar effects on 24-h nutrient

287 oxidation. Med Sci Sports Exerc. 2002; 34:1793-800.

288 Nadel ER. Energy exchanges of swimming man. J Appl Physiol. 1974; 36: 465-71.

289 Sarang P, Telles S. Oxygen Consumption and During and After two Yoga

290 Relaxation Techniques. Appl Psychophysiol Biofeedback. 2006; 31:143-53.

291 Schagatay E. Effects of physical and apnea training on apneic time and the diving

292 response in humans. Eur J Appl Physiol.2000; 82:161-9

293 Schagatay E. Predicting performance in competitive apnoea diving. Part I: static apnoea.

294 Diving Hyperb Med. 2009, 39: 88-99.

295 Schagatay E. Predicting performance in competitive apnea diving, part II: dynamic

296 apnoea. Diving Hyperb Med. 2010; 40: 11-22.

12 13

297 Seccombe, L., Features of glossopharyngeal breathing in breath-hold divers. J Appl

298 Physiol, 2006. 101: 799-801.

299 Tabachnick BG, FidellL S. Using multivariate statistics. Education Inc, eds. New York;

300 2007; 33-50.

301 Wanger J. Standardization of the measurement of . Eur Respir J. 2005; 26:

302 511-22.

303 Whittaker L. Going to extremes of lung volume. J Appl Physiol. 2007; 102: 831-833.

304 Woorons X. Exercise with hypoventilation induces lower muscle oxygenation and

305 higher blood lactate concentration: role of hypoxia and hypercapnia. Eur J Appl

306 Physiol. 2010; 110: 367-77.

307

308 Figure Legends

309 Figure 1a-d

310 Profile plots of the static apnea – STA (1a) -, no fins- DNF (1b) -,

311 dynamic apnea with fins – DYN (1c) – and the global performance indicator–POINTS

312 (1d)-. The estimated marginal means of each performance variable are the mean

313 responses of the performance variable for each combination of level pairs of both

314 factors: training group and time measure.

315

316 Figure 2

317 The graph depicts the subjects’ scores on both discriminant functions (Function 1 and

318 Function 2 over X- and Y-axis respectively) as well as centroids of each group. The

319 centroid are calculated by averaging all subjects’ scores belonging to the same training

320 group.

321

13 14

322 Table Legends

323 Table 1

324 Means and standard deviations (SD) of the scores obtained before and after the training

325 interventions as well as the differences between cross training (CT), apnea training (AT)

326 and control group (CG). POINTS variable, used as a performance indicator, is

327 constructed from the variables STA, DNF and DYN. Significant differences between

328 pre- post were marked by * (p<0.05), ** (p< 0.01) and NS (no significant differences).

329

330 Table 2

331 Means and standard deviations (SD) of the scores obtained before and after the training

332 interventions as well as the differences between cross training (CT), apnea training (AT)

333 and control group (CG). VO2max, maximum oxygen consumption; VC, vital capacity;

334 Hb, hemoglobin; BFP, body fat percentage; LBM, lean body mass; RMR, resting

335 metabolic rate; STAdry, static apnea time; HRavg, average heart rate during static apnea;

336HR min, minimum heart rate during static apnea. Significant differences between pre-

337 post were marked by * (p<0.05), ** (p< 0.01) and NS (no significant differences).

14 1. In order to increase dynamic apnea performance, breath-hold divers should introduce as much as possible physical conditioning associated with apnea training.

2. Physical conditioning, associated to apnea training, increases Vital Capacity and VO2max, significantly influencing apnea performance.

3. There is no significant evidence about the influence of apnea training on haematological values or body composition

1 1

1 Medium term effects of physical conditioning on breath-hold diving performance

2

3 Abstract

4 The current study aimed to analyze the effects of physical conditioning inclusion on

5 apnea performance after a 22-week structured apnea training program. Twenty-nine

6 male breath-hold divers participated and were allocated into: (1) cross-training in apnea

7 and physical activity (CT; n=10); (2) apnea training only (AT; n=10); and control group

8 (CG; n=9). Measures were static apnea (STA), dynamic with fins (DYN) and dynamic

9 no fins (DNF) performance, body composition, hemoglobin, vital capacity (VC),

10 maximal aerobic capacity (VO2max), resting metabolic rate, oxygen saturation, and

11 pulse during a static apnea in dry conditions at baseline and after the intervention. Total

12 performance, referred as POINTS (constructed from the variables STA, DNF and DYN)

13 was used as a global performance variable on apnea indoor diving. +30, +26 vs. +4

14 average POINTS of difference after-before training for CT, AT and CG respectively

15 were found. After a discriminant analysis, CT appears to be the most appropriate for

16 DNF performance..The post-hoc analysis determined that the CT was the only group in

17 which the difference of means was significant before and after training for the VC

18 (p<0.01) and VO2max (p<0,05) variables. Inclusion of physical activity in apnea training

19 increased VC and VO2max in breath hold divers; divers who followed a mixed training,

20 physical training and hypoxic training, achieved increased DNF performance.

21

22 Keywords: apnea, hypoxia, training, respiratory.

1 1

1 Medium term effects of physical conditioning on breath-hold diving performance

2

3 1. Introduction

4 In apnea indoor diving competition, breath-hold divers compete to remain immersed

5 as long as possible in static apnea (STA), and to dive the longest distance in the

6 categories dynamic with fins (DYN) and dynamic no fins (DNF). Adverse effects, such

7 as decompression sickness [Batle, 1999], narcosis or arterial embolism [Batle, 2002],

8 observed during deep events, are not going to occur in the pool practice where the diver

9 submerges barely 2 meters; thus, hypoxic syncope is the main risk during apnea indoor

10 practice [Fitz-Clarke, 2006].

11 Recent improvements in equipment, nutrition, and training have led to increases in

12 apnea performance [Schagatay, 2009; Schagatay, 2010; Fernández, 2015]. Practice and

13 training is mainly responsible for performance in breath-hold divers; however, it is

14 difficult to establish a rigid training model that suits all divers with different levels of

15 ability, age, previous physical condition or anthropometry.

16 The effects of physical activity on apnea performance have not been clarified,

17 probably because the effects of physical training and apnea training were studied

18 independently. The effects of physical training on apnea performance are controversial.

19 While Bavagad showed in his study [Bagavad, 2014] a correlation between physical

20 conditioning with apnea performance, Schagatay [Schagatay, 2000] did not find effects

21 on apnea performance; however, this latter research concluded that physical training

22 increased stamina, allowing prolongation of the struggle phase of apnea. In addition, a

23 few longitudinal research studies of apnea training have been conducted [Schagatay,

24 2000; Bagavad, 2014; Fernández, 2015] .The current study aimed to analyze the effects

1 2

25 of physical conditioning inclusion on apnea performance after a 22-week structured

26 apnea training program.

27

28 2. Materials and Methods

29 2.1. Participants

30 Twenty-nine male breath-hold divers (36 ± 5 years of age) with two years of

31 experience in breath-hold diving were divided in different training programs: (I) cross-

32 training in apnea and physical conditioning (CT; n=10); (II) apnea training (AT; n=10);

33 and (III) control group (CG; n=9). No training intervention was carried out on the CG;

34 however, regarding pre-post measurements, both the CG and the experimental groups

35 performed them without any distinction between them. The participants were informed

36 of the benefits and risks prior to signing the informed consent document to participate in

37 the research. Participants diagnosed with cardiac, metabolic or respiratory diseases were

38 excluded. The study was conducted in accordance with the Declaration of Helsinki

39 [Harris, 2011] and approved by the institutional Human Research Ethics Committee

40 (CSEULS-PI-114/2016).

41

42 2.2. Training procedures

43 All participants conducted 66 sessions of 60-minute of duration each. During the first

44 session, CT was performed for 15 minutes with a dry strength circuit of 10 calisthenic

45 exercises - squat, pull-up, push-up, squat, pull-up, push-up, squat, pull-up, push-up and

46 squat- with 50 repetitions and 10 seconds of recovery. This was followed by swimming

47 training in crawl, alternating each week between 45 min of swimming continuous

48 training (14 RPE) and high intensive interval training (HIIT) (20 RPE). Regarding

49 HIIT, volume and density were 3 series x 10 repetitions of 25 meters with 20 seconds of

2 3

50 recovery. During the second session, participants performed a specific hypoxic training

51 in dynamic, that consisted of 40 repetitions of 25 m dynamic apneas to maximum

52 individual underwater swimming speed with 20 s of recovery each 25 m. In the third

53 session, divers conducted a specific hypoxic training in static. In this training divers

54 performed 3 maximal static apneas with 10 min of full recovery between each apnea.

55 The AT group, performed a specific hypoxic training in dynamic, on the first and

56 second session, and a specific hypoxic training in static, on the third session.

57 Glossopharyngeal insufflation (Seccombe, 2006) was not allowed during training for

58 any group.

59

60 2.3. Test procedures

61 Measurements were performed under similar environmental conditions in a sports

62 laboratory, a health center and a 25-meter pool with 2-meter depth. The test battery

63 sequence was performed chronologically in five visits in different days, as follows:

64 Visit 1– STA; Visit 2 – DNF; Visit 3– DYN; Visit 4 – body composition, blood count,

65 spirometry and an incremental test on treadmill; Visit 5 – resting metabolic rate, and

66 pulse-oximetry during a static apnea in dry conditions.

67

68 2.1.1. STA, DYN and DNF performance

69 Same protocol was performed, previous to apnea, to avoid possible warming effects.

70 Thus, before STA, a 15-minute warm-up consisted of 10 minutes of relaxation, two

71 minutes of static apnea and a three-minute countdown on the surface until maximal

72 attempt. Regarding the DYN and DNF, a 15-minute warm-up was performed

73 comprising 10 minutes of relaxation, 50 meters in the specific (fins or no fins) dynamic

74 discipline and a three-minute countdown on the surface. After warm-up, the divers

3 4

75 attempted to achieve the maximal individual time or distance. In all the disciplines,

76 during the last 30 seconds of the countdown, the participant placed the nose clip and

77 performed a maximal inspiration.

78

79 2.1.2. Physiological variables

80 Body fat and fat-free mass were measured by whole-body dual-energy X-ray

81 absorptiometry (GE Lunar Prodigy; GE Healthcare, Madison, Wisconsin).Regarding

82 hematological values, hemoglobin (Hb) and hematocrit were analyzed by a blood count

83 analysis at the health center and under fasting conditions; a portable spirometry test

84 (Spirostik, Geratherm) was performed as well to calculate vital capacity (VC). Before

85 starting the test, we explained the technique to the subject: sitting, with knees bent and

86 feet resting on the floor. In this way, the nasal clip was placed and the participant would

87 make a relaxed, but maximum inspiration, at the end of which the mouthpiece would be

88 placed; then, he was to perform slow expiration, maintaining a constant flow until

89 residual volume. This procedure was be repeated three times, with at least 30 seconds of

90 rest between each measurement and the highest value of the three measurements was

91 registered [Wanger, 2005].

92 The incremental test to measure maximal oxygen consumption(VO2max)and maximal

93 heart rate was performed on a treadmill (H/P/COSMOS 3P® 4.0, H / P / Cosmos Sports

94 & Medical, Nussdorf-Traunstein, Germany). The volume and composition of expired

95 gases were measured using a gas analyzer (Ultima CPX, Medical Graphics) and the

96 heart rate measured by electrocardiogram (WelchAllyn, CardioPerfect). After a five-

97 minute warm-up at 10 km•h-1, the speed was increased 1 km•h-1 every minute until

98 volitional exhaustion of the participant. Throughout the test the treadmill elevation was

99 maintained at 1%. VO2max was noted as the average of the two highest values of 15

4 5

100 seconds of VO2 consecutive reached toward the end of the test, as long as the values of

101VO 2 and HR have stabilized despite the increase in the intensity of the effort.

102 Resting metabolic rate was measured, breath-to-breath (Ultima CPX, Medical

103 Graphics), over 15 minutes by indirect calorimetry. According to a previous study

104 [Melanson, 2002], a minimum of 15 minutes of steady state, determined as <10%

105 fluctuation in oxygen consumption and <5% fluctuation in respiratory exchange ratio

106 between CO2 and O2, was considered criteria for valid resting metabolic rate(RMR).

107 Participants performed the test in fasting conditions after 12 hours without physical

108 effort. The environmental conditions established were as follows: a quiet laboratory

109 with the lights off, 550 meters of altitude, 22.5 ± 1°C, with 55 ± 3% relative humidity.

110 Heart rate and oxygen saturation were monitored during a maximal static apnea in

111 dry conditions. Before beginning the test, the diver rested for 10 minutes in the prone

112 position, with head and upper limbs leaning on a table placed in front of the stretcher.

113 During the last 30 seconds of the countdown, a nose clip was placed to avoid possible

114 air leakage. At that point, the participant performed a maximal exhalation followed by a

115 maximal inspiration; glossopharyngeal insufflation was not allowed. Throughout the

116 static apnea test, the average heart rate (HRavg) and the minimum heart rate (HRmin)

117 were recorded. To calculate the HRavg, only the heart rate data after first 30 seconds

118 were analyzed; i.e., once the heart rate was stabilized [Breskovic, 2012]. The oxygen

119 saturation was monitored by a pulse oximeter (CMS 50F) placed on the second finger of

120 the left hand. A minimal value of oxygen saturation was collected during the test

121 (SpO2min). Time during static apnea in dry conditions (STAdry) was recorded, also.

122

123 2.3. Statistical analyses

5 6

124 The effects of training on performance and physiological variables were analyzed

125 using repeated measures analysis of variance (ANOVA). In order to improve the power

126 of the test and to reduce Type I error, a repeated measures multiple analysis of

127 variance(MANOVA) was used in the case of a set of variables defined in a cohesive

128 theme [Tabachnick, 2007]. Moderate correlations(r=0.3 to 0.7) led to determine the

129 group of variables: STA, DYN and DNF in the case of performance variables, whereas

130 the groups VO2max, Hb and VC (related to the lung storage), percent of body fat (PBF)

131 and lean body mass (LBM) (measuring body composition) and RMR and HRavg

132 (corresponding to the individual metabolic rate) were established with the physiological

133 variables.

134 A post-hoc analysis was performed when significant differences were detected

135 between the groups. Regarding ANOVA, the specific contrast tests were used, when the

136 homogeneity of the variances was achieved, with Bonferroni correction. When

137 interaction occurs, the effect of one factor on a response depends on the level of

138 (an)other factor(s). Thus, we cannot consider the main effects of the factors separately

139 as the main effects and interaction need to be considered as a whole to describe the

140 relationship between input and output. However, discriminant analysis was used as a

141 post-hoc procedure of the MANOVA, with the objective of determining which

142 dependent variables discriminate between the training programs and analyzing their

143 contributions. The significance level considered in all the developed tests was 0.05.

144

145 3. Results

146 Tables 1 and 2 summarize the mean and standard deviations (SD) of the scores

147 obtained before and after the training interventions as well as the differences between

148 them. (Tables 1 and 2)

6 7

149 3.1. Performance variables

150 Total performance, referred as POINTS in this work, was used as a global performance

151 variable on a Apnea Indoor diving. It is constructed from the variables STA, DNF and

152 DYN. In particular, each second the athlete remains immersed in STA was multiplied

153 by 0.2, whereas each meter reached at DNF and DYN was multiplied by 0.5. AIDA

154 International – one of the official freediving competition organizing bodies – uses these

155 scores during championships: The diver with the highest score wins the competition.

156 The analysis of variance of repeated measured data for the POINTS variable indicated

157 the effect of the training group (F = 6.270, p<0.01), time (F = 45.019, p<0.001) and the

158 interaction between both (F = 7.409, p<0.01) to be statistically significant. Profile

159 graphs (Figure 1) show weak discrepancies between CT and AT estimated marginal

160 means. Nevertheless, there are notable differences between control group and both CT

161 and AT estimated marginal means. In Table 1, significant differences between pairs of

162 measures of both time and training group factors were collected (Figure 1).

163 Group linear discriminant analysis identifies a set of linear combinations of

164 dependent variables (LDF) that discriminate groups. After a stepwise regression the

165 variables included in the model were STA2 and DNF2, leaving out DYN2. The number

166 of LDFs is the number of groups minus 1; thus two discriminant functions were

167 necessary to achieve significant discrimination: Function 1 = 0.575·ZSTA2 + 0.539·ZDNF2

168 and Function 2 = -1.126·ZSTA2 + 1.144·ZDNF2. , where ZX is the standardized score of the

169 variable X. Both LDFs were statistically significant (p<0.001 and p<0.01 respectively).

170 The coefficients of the LDFs are a measure of the importance of each variable for

171 distinguishing groups. Thus, the first function these coefficients indicates that STA2

172 was more important in distinguishing, while DNF2 had a larger contribution for the

173 second function. The graph in Figure 2 shows the subjects’ scores on the discriminant

7 8

174 functions and centroids of each group. We can see from the graph that the centroids are

175 positive in the first function for the groups CT and AT and negative for CG. Therefore,

176 the first function classifies as a control the individuals with the lowest score in the

177 variables STA2 and DNF2, whereas individuals with higher scores would be classified

178 in the CT group and, to a lesser extent, in the AT group. In the second function, we have

179 opposite sign coefficients and similar values for the discriminant variables. Because the

180 centroids that are the most distant are those that correspond to the groups CT and AT,

181 those of opposite sign, the function classifies the participants with greater scores in

182 STA2 and smaller scores in DNF2 within the AT group. In contrast, the participants

183 with a lower score in STA2 and a higher score in DNF2 are classified in the CT group

184 (Figure 2).

185 3.2. Physiological variables

186 The post-hoc analysis determined that the CT group was the only group in which the

187 difference of means was significant before and after training for the VC and VO2max

188 variables. We found a significant increase in VC and VO2max after training (see Table 2).

189 Mean of HRmin significantly decrease after training for CT and AT (see Table 2).

190 Analyzing this variable by pairs of training groups, significant differences were found

191 between the control group and the rest of the groups, being the mean of HRmin higher in

192 the CG than in the others (see Table 2).

193

194 4. Discussion

195 There is current interest in the role of physical training on apnea performance.

196 Previous research [Schagatay, 2000] analyzed the effects of physical training and apnea

197 training separately. In the current study, after the discriminant analysis, we concluded

198 that training groups can be distinguished in two areas: isolated apnea training appears to

8 9

199 be the most suitable for STA, whereas cross-training appears to be best for DNF.

200 Regarding DYN performance, although there were no significant differences between

201 the training groups (AT vs. CT), there was a trend (favorable to CT).

202 As hypothesis, the inclusion of physical training could have triggered an

203 improvement in the diver’s underwater swimming economy; however, this variable has

204 not been measured in our study. Another study [Nadel, 1974] proposed that trained

205 swimmers use less oxygen than those not trained to swim at the same speed. Also, a

206 previous review [Schagatay, 2010] concluded that a high proportion of ex-swimmers

207 represent the elite in the apnea dynamic discipline.

208 According to the correlation shown between physiological variables studied, groups

209 were defined into: lung storage, body composition, metabolic rate, and individual

210 tolerance to asphyxia. The difference in means after training interventions was

211 significant only in the CT group for the VC and VO2max variables.

212 Maximal aerobic capacity is closely related to performance in endurance sports

213 (Bassett, 2000). Expected benefits, after physical training, would be an increased

214 storage of oxygen before apnea, a decreased heart rate during apnea and increased

215 motivation or stamina. According to the results in this study, the increased aerobic

216 capacity did not result in a lower heart rate during apnea, which is consistent with the

217 findings of Schagatay (Schagatay, 2000). On the other hand, the test for assessing

218 maximum aerobic capacity (incremental test on treadmill) was no specific (swimming);

219 thus, muscular fatigue could produce limitation on aerobic and anaerobic performance

220 during the incremental test on treadmill.

221 Physical training resulted in an increased VC; thus, the inclusion of physical activity

222 in the apnea training plan could be considered appropriate for increasing divers’ total

223 O2/CO2 lung storage. Glossopharyngeal insufflation was not allowed during training or

9 10

224 spirometry test; thus, there was no any improvement or adaptation produced by this

225 breathing technique. Whittaker noted in his study [Whittaker, 2007] that after a specific

226 training consisting of thoracic stretching with glossopharyngeal insufflation divers

227 could increase total lung capacity up to 2 liters. However, the proposed training was

228 insufficient to cause long-term adaptations in hematological values. In addition to

229 hypoxia training, it would be advisable to include hypoventilation training in low lung

230 volume instead of the traditional training in hypercapnia, mentioned in our study. A

231 study conducted in 2010 [Woorons, 2010] showed that hypoventilation training in low

232 lung volume produced both an increase in CO2 concentrations (hypercapnic effect) and

233 a decrease in O2 (hypoxic effect) in the blood and muscles.

234 No tendency was found regarding body composition after training programs,

235 suggesting that AT or CT was insufficient to cause significant long-term changes in

236 body composition. Previous study [Sarang, 2006] proposed a decrease greater than 30%

237 in O2consumption after the application of relaxation techniques; however, achieving the

238 concentration that enables breath-hold divers to reduce their metabolic rate during apnea

239 requires many hours of deliberate practice and years of experience.

240 Decreased HR min after AT and CT, compared with the CG, is probably due to a

241 pronounced bradycardia induced by increased hypoxia during the static apnea in dry

242 conditions; i.e., breath-hold divers who followed periodized training achieved greater

243 tolerance to asphyxia and consequently an increase in apnea time. On the other hand,

244 there was no trend for the HRmin between experimental groups, suggesting a lack of

245 cardiac adaptation produced by physical exercise proposed. Physical and psychological

246 tolerance to effort have aspects in common during the struggle phase of apnea or

247 physical exertion; thus, a intensive physical training, as performed through HIIT, could

248 trigger an increased tolerance to the asphyxia produced during apnea. According to

10 11

249 results showed in our study, oxygen saturation reached during dry static apnea were far

250 from limits, i.e.; hypoxic syncope (50% SaO2 in amateurs and 30% in trained freedivers

251 [Andersson, 2004]). In this way, an increase in tolerance to asphyxia in our divers could

252 dilute the effect (on Apnea Indoor performance) of other factors such as: total storage of

253 O2/CO2, metabolic rate or underwater swimming economy. A limitation of the present

254 study was that O2 consumption and CO2 release rates were not calculated during the

255 breath-hold; moreover, there are no data on the of O2 and CO2 in the

256 alveolar air before the apnea and at the end of the apnea. The availability of these data

257 would provide a more detailed analysis of the identified correlations.

258

259 5. Conclusions

260 The inclusion of physical activity in apnea training increased VC and VO2max in breath

261 hold divers; in addition, divers who followed a mixed training, physical training and

262 hypoxic training, achieved increased DNF performance.

263

264 References

265 Andersson J. Cardiovascular and respiratory responses to apneas with and without face

266 immersion in exercising humans. J Appl Physiol. 2004; 96:1005-10

267 Bagavad G. Effect of physical training on breath holding time in Indian subjects. Indian

268 J Physiol Pharmacol. 2014; 58: 108-109.

269 Bassett DR. Limiting factors for maximum oxygen uptake and determinants of

270 endurance performance. Med Sci Sports Exerc. 2000; 32: 70-84.

271 Batle JM. Decompression sickness caused by breath-hold diving hunting. A report of

272 about 30 cases. Proceedings of the 13th International Congress on Hyperbaric

273 Medicine, 1999.

11 12

274 Batle JM. How to avoid arterial gas embolism in free diving, the attitude a breath-hold

275 diver should take. A study of 35 cases, in a period 1995-2002. Proceedings of

276 the 14th International Congress on Hyperbaric Medicine, 2002.

277 Breskovic T. Cardiovascular changes during underwater static and dynamic breath-hold

278 dives in trained divers. J Appl Physiol. 2012; 3:673-8.

279 Fernández FA. Periodization of Apnea Training. Doctoral thesis, University of Castilla-

280 la Mancha, Toledo, 2015.

281 Fitz-Clarke, J. Adverse events in competitive breath-hold diving. Undersea Hyperb

282 Med. 2006 Jan-Feb; 33(1): 55-62.

283 Harris D, Atkinson G. Update - Ethical standards in sport and exercise science research.

284 Int J Sports Med. 2011; 32: 819-21.

285 Melanson EL et al. Resistance and aerobic exercise have similar effects on 24-h nutrient

286 oxidation. Med Sci Sports Exerc. 2002; 34:1793-800.

287 Nadel ER. Energy exchanges of swimming man. J Appl Physiol. 1974; 36: 465-71.

288 Sarang P, Telles S. Oxygen Consumption and Respiration During and After two Yoga

289 Relaxation Techniques. Appl Psychophysiol Biofeedback. 2006; 31:143-53.

290 Schagatay E. Effects of physical and apnea training on apneic time and the diving

291 response in humans. Eur J Appl Physiol.2000; 82:161-9

292 Schagatay E. Predicting performance in competitive apnoea diving. Part I: static apnoea.

293 Diving Hyperb Med. 2009, 39: 88-99.

294 Schagatay E. Predicting performance in competitive apnea diving, part II: dynamic

295 apnoea. Diving Hyperb Med. 2010; 40: 11-22.

296 Seccombe, L., Features of glossopharyngeal breathing in breath-hold divers. J Appl

297 Physiol, 2006. 101: 799-801.

12 13

298 Tabachnick BG, FidellL S. Using multivariate statistics. Education Inc, eds. New York;

299 2007; 33-50.

300 Wanger J. Standardization of the measurement of lung volumes. Eur Respir J. 2005; 26:

301 511-22.

302 Whittaker L. Going to extremes of lung volume. J Appl Physiol. 2007; 102: 831-833.

303 Woorons X. Exercise with hypoventilation induces lower muscle oxygenation and

304 higher blood lactate concentration: role of hypoxia and hypercapnia. Eur J Appl

305 Physiol. 2010; 110: 367-77.

306

307 Figure Legends

308 Figure 1a-d

309 Profile plots of the static apnea – STA (1a) -, dynamic apnea no fins- DNF (1b) -,

310 dynamic apnea with fins – DYN (1c) – and the global performance indicator–POINTS

311 (1d)-. The estimated marginal means of each performance variable are the mean

312 responses of the performance variable for each combination of level pairs of both

313 factors: training group and time measure.

314

315 Figure 2

316 The graph depicts the subjects’ scores on both discriminant functions (Function 1 and

317 Function 2 over X- and Y-axis respectively) as well as centroids of each group. The

318 centroid are calculated by averaging all subjects’ scores belonging to the same training

319 group.

320

321 Table Legends

322 Table 1

13 14

323 Means and standard deviations (SD) of the scores obtained before and after the training

324 interventions as well as the differences between cross training (CT), apnea training (AT)

325 and control group (CG). POINTS variable, used as a performance indicator, is

326 constructed from the variables STA, DNF and DYN. Significant differences between

327 pre- post were marked by * (p<0.05), ** (p< 0.01) and NS (no significant differences).

328

329 Table 2

330 Means and standard deviations (SD) of the scores obtained before and after the training

331 interventions as well as the differences between cross training (CT), apnea training (AT)

332 and control group (CG). VO2max, maximum oxygen consumption; VC, vital capacity;

333 Hb, hemoglobin; BFP, body fat percentage; LBM, lean body mass; RMR, resting

334 metabolic rate; STAdry, static apnea time; HRavg, average heart rate during static apnea;

335HR min, minimum heart rate during static apnea. Significant differences between pre-

336 post were marked by * (p<0.05), ** (p< 0.01) and NS (no significant differences).

14

Fernández. Table.1

STA DNF DYN POINTS (s) (m) (m)

CT(n=10) PRE 185±46 65±14 71±18 100±24 POST 249±43 79±17 90±25 130±29 Dif 63±36** 14.50±14** 19±15** 30±15**

AT (n=9) PRE 182±46 55±12 65±13 96±19 POST 250±33 64±13 79±13 122±14 Dif 68±35** 9±15* 14±9** 26±17**

CG(n=10) PRE 152±46 45±10 59±9 82±15 POST 160±56 47±14 57±15 86±23 Dif 8±31NS 3±11NS -2±13NS 4±15NS Fernández. Table.2

VO2 max VC Hb BFP LBM RMR STA dry HRavg HRmin SpO2min (ml/kg/min) (L) (g/dl) (%) (%) (cal*min) (s) (%HRmax) (%HRmax) (%)

CT(n=10) PRE 47±4 5.33±0.72 15.9±1.2 26.2±6.3 70.8±5.9 1.1±0.2 184±45 30± 3 24±4 90±5 POST 51±6 5.88±1.13 15.9±1.1 25.4±6.1 71.5±5.7 1.1±0.3 206±34 29± 2 21±3 86±7 Dif 4±3* 0.55±0.50** -0.1±0.5NS -0.8±1.8NS 0.8±1.7NS 0.0±0.2NS 22±21** -1 ± 4NS -3±3** -4±8NS

AT (n=9) PRE 47±6 5.90±1.26 15.4±1.2 25.0±7.2 72.0±6.7 1.1±0.2 175±68 31 ± 3 22±3 85±11 POST 47±5 6.01±1.45 15.2±1.0 24.6±7.4 72.3±6.8 1.2±0.2 212±36 32 ± 3 20±1 84±11 Dif 0±4NS 0.11±0.45NS -0.3±0.5NS -0.3±0.8NS 0.3±0.7NS 0.1±0.2NS 37±26** 1± 2NS -2±2* -1±9NS

CG(n=10) PRE 53±10 5.20±1.04 15.5±0.8 22.4±7.0 74.4±6.6 1.2±0.2 167±55 31±5 27±3 90±8 POST 50±8 5.27±1.10 15.4±0.8 23.4±7.8 66.4±25.6 1.2±0.3 152±58 31±4 27±2 88±7 Dif -2±5NS 0.07±0.23NS -0.1±0.6NS 1.0±4.3NS -8.1±24.2NS 0.0±3.6NS -15±4NS 0±4NS 0±0NS -1±6NS