The breaststroke technique is considered one of the least economic of the four swimming styles (Barbosa et al., 2006), which can lead to early fatigue while swimming. Neuromuscular fatigue can be defined as: (i) the failure to maintain the required or expected force, accompanied by changes in muscular activity (Dimitrova et al., 2003); and (ii) the inability of skeletal muscle to generate high levels of muscular strength or maintain these levels over time (Enoka and Stuart, 1992). Additionally, the manifestations of fatigue have been associated with (Allen et al., 1995; Pagala et al., 1994): (i) a decline in muscle tension produced during and after submaximal and maximal exercise; (ii) an inability to maintain a given exercise intensity over time, reducing the speed of contraction and increased muscle relaxation time, as well as; (iii) the variation of intra and extracellular concentrations of certain metabolites and ions.

Neuromuscular fatigue can be divided into central and peripheral fatigue. Central fatigue comprises of decreases in the voluntary activation of the muscle, which is due to decreases in the number of recruited motor units and their discharge rate (González et al., 2012). Peripheral fatigue involves changes in neuromuscular transmission, muscle action potential propagation and decreases in the contractile tension of the muscle fibers (Boyas and Guevel, 2011). Peripheral fatigue during exercise is considered an impairment of the peripheral mechanisms from excitation to muscle contraction. Peripheral regulation is, therefore, related to a perturbation of calcium ion movements, an accumulation of phosphate, and/or a decrease of adenosine triphosphate stores (Boyas and Guevel, 2011).

Some studies attempted to relate the onset of fatigue with the execution of various sport techniques (Camata et al., 2011; Giangarra et al., 1993; Pink et al., 1993; Robineau et al., 2012). While swimming, neuromuscular mechanisms related to fatigue remain unclear, and the few studies which carried out research, focused mainly on the front-crawl technique (Caty et al., 2006; Figueiredo et al., 2011; Ikuta et al., 2012; Stirn et al., 2011).

Caty et al., 2006 observed a decrease in the instantaneous mean frequency in the extensor carpi ulnaris muscles (11. 4% and 8.5%, respectively) after a 4x50m high intensity front crawl. Stirn et al., 2011 reported that at the end of a 100m front crawl at a maximal effort, the mean power frequency decreased by 20-25%. Ikuta et al., 2012 suggested that the decrease in swimming velocity was related to a decrease in the activity of several muscles coordinated with each other, and that a compensating strategy was involved between the pectoralis major and other muscles during the last lap of a 4x50m front crawl test.

Important muscles activated in breaststroke swimming seem to be the biceps brachii, triceps brachii, subscapularis, latissimus dorsi, pectoralis major, supraspinatus, infraspinatus, serratus anterior, and deltoid anterior, teres minor and trapezio ) for the upper limbs (Conceição et al., 2010; Nuber et al., 1986; Ruwe et al., 1994; Yoshizawa et al., 1976). While gluteus maximus, vastus medialis, rectus femoris, biceps femoris, abductor magnus, quadriceps, gastrocnemius, tibialis anterior, abductor hallucius, abductor digiti minimi, flexor digitorum brevis are for the lower limbs (McLeod, 2010; Yoshizawa et al., 1976).

It is clear that EMG can be useful in tracking muscle fatigue for many reasons (the relationships found between sEMG features and muscle fatigue, and the possibility of recording them in almost any type of situation). Therefore, some studies have examined the relationship between muscle fatigue and EMG variables and, consequently, the possibility of using EMG models to accurately track muscle fatigue (González-Izal et al., 2012).

The linear techniques that are used to estimate muscle fatigue are based on linear regression, which relates changes in EMG parameters to changes in power loss (as a direct measurement of muscle fatigue). The non-linear techniques that are used to estimate muscle fatigue are based on neural networks to relate EMG parameters to muscle fatigue. Both linear and non-linear techniques provided good results for mapping changes in power or force loss during dynamic exercises based on sets of EMG parameters (González-Izal et al., 2012).

Apart from this, however, we have no idea of the appearance of fatigue of the upper limbs in breaststroke. Since there is a lack of studies to quantify the neuromuscular fatigue in the upper limbs, using spectral parameters, we believe that this research would have the potential to make a contribution to the limited body of knowledge in this field. The main aim for swimming coaches and researchers relies on the identification of the factors that might predict performance with higher validity and accuracy.

The aims of this study were: i) to analyze the activation patterns (duration of active and non-active phase) of the upper limbs’ four muscles during each lap for 200m breaststroke, ii) quantify neuromuscular fatigue, through kinematics and physiologic assessment. It was hypothesized that an increase in signal amplitude and a decrease in spectral parameters due to repetitive sub-maximal contractions, characterized a non-linear fatigue process. Also considered was that the fatigue process occurs differently for each of the muscles studied.

Nine male swimmers (age: 22.3 ± 2.9 years; height: 1.81 ± 0.05 m; body mass: 73.60 ± 3.82 kg; mean ± SD) volunteered to participate in this study and provided written consent. They were all swimmers competing at the national level with an average personal best result for 200-m breaststroke (149.44 ± 6.59 s, corresponding respectively to 643.75 ± 53.77 FINA ranking points). All the procedures were approved by the institutional Ethics Committee and carried out according to the Helsinki Declaration.

The experiments were performed in a 50m indoor swimming pool at a water temperature of 27.5 °C and 75% humidity. Subjects performed a standard warm-up of 800m front crawl, and a specific warm-up of 200m breaststroke at a medium level of effort. After a twenty minute rest subjects performed a maximal 200m breaststroke trial with a push off start.

Surface EMG signals from the biceps brachii (BB), deltoid anterior (DA), pectoralis major (PM) and triceps brachii (TB) on the right side of the body were collected. These muscles were selected according to their importance in breaststroke (Conceição et al., 2010; Nuber et al., 1986; Ruwe et al., 1994; Yoshizawa et al., 1976).

Bipolar surface electrodes were used (10-mm diameter discs, Plux, Lisbon, Portugal) with an inter-electrode distance of 20mm. Electrodes on the upper part of the PM were placed in the middle of the line that connects the acromion process and the manubrium (sternum) two fingers below the clavicle (Stirn et al., 2011). The electrodes on the long head of the BB, DA and TB were placed in accordance with SENIAM recommendations (Herrmens and Freriks, 1999).

The skin was shaved, rubbed with sandpaper and cleaned with alcohol so that the inter-electrode resistance did not exceed 5 KOhm (Basmajian and De Luca, 1985). The ground electrode was positioned over the cervical vertebrae. Transparent dressings with labels (Hydrofilm^®, 10 cm x 12,5 cm, USA) were used to cover the electrodes to isolate them from the water (Hohmann et al., 2006). All cables were fixed to the skin by adhesive tape in several places to minimize their movement and any signal interference. Swimmers wore a full body swimming suit to further immobilize the cables (Fastskin Speedo^®, Speedo Aqualab, USA).

The wireless EMG device (BioPLUX.research, Lisbon, Portugal) with 8 analog channels (12-bit), sampling rate at 1 kHz, was put in a waterproof bag (84x53x18mm) and placed inside the swimmer´s cap. Data was transmitted to the PC in real time via Bluetooth.

The EMG and video data were recorded simultaneously and were not synchronized with video recording, because the aim of this study was to understand the appearance of fatigue in relation to race parameters rather than detailed within cycle swimming kinematics.

All EMG analysis was conducted with a MATLAB routine (Mathworks, Inc.,Natick MA,USA). The process of determining the muscle activity boundaries consists of finding the neighborhood points, where the energy was 30% of muscle activation maximum peak within a stroke (Stirn et al., 2011). These were calculated by segmenting the muscle input signal energy according to the same criteria described in Stirn et al., 2011. Starting from the raw signal, DC components were removed and thereafter filtered with a fifth-order Butterworth band pass filter where the lower and upper cut-off frequencies were set to 10 and 500Hz respectively. The signal energy was then determined with a 250 ms sliding window (Stirn et al., 2011) and according to: (1) Even though the high frequencies of the input signal were filtered with a Butterworth filter, muscle energy is very noisy and presents several local maximum peaks that didn´t correspond to the muscle active window center, as shown in Figure 2. To overcome this difficulty, a strategy to determine the muscle’s “true” maximum energy peaks was devised. Each stroke taken by a swimmer produces patterns in the signal, these patterns are mainly translated by a periodicity in EMG energy, see Figure 1.

By defining the signal mean period, one can use this information to determine the maximum peak candidates with the highest and with minimal differences between two maximum candidates and the expected period.

Once the maximum candidates have been determined, the muscle activity boundaries were then selected by determining the neighborhood points where the energy is 30% of the determined maximum peaks. For each muscle activation, the active phase was defined as the part of the EMG signal for which the energy was at least 30% of the local maximum energy value, for each particular muscle activation. The raw EMG segments belonging to the active phases were extracted and used in the calculation of the active phase duration and the amplitude frequency analysis. The non-active phase was defined as the time interval between the two consecutives active phases as shown in Figure 2.

The temporal evolution of the active and non-active phases average durations during stroke were calculated for each muscle for the entire 200m. Linear regression curves were fitted to the data and the durations of the fitted curves at the beginning and completion of the swim were compared. The comparisons of the durations of the active and non-active phases of each muscle during the stroke were made to detect changes that might occur due to fatigue.

The average amplitude of EMG of each active phase was estimated using the average rectified value (ARV) of the EMG. ARV was calculated in accordance with SENIAM recommendations (Herrmens and Freriks, 1999) and plotted as a function of time. A linear regression curve was fitted to the data and the ARV values of the curve at the time of commencement and the end of the swimming bout were compared (Figure 3).

For the frequency analysis, each extracted segment was zero-padded to the total length of 1s (i.e., 1,000 samples). In this way, a uniform frequency resolution (the frequency sampling interval) was used for all signal segments. The power spectral density (PSD) for each segment was estimated using the periodogram method (Proakis and Manolakis, 1996). The periodogram of a continuous signal segment x(t) of length T is defined as: (2) Even though the periodogram is a non-consistent estimate of PSD (its variance is large and does not become zero with increasing length of the signal), it was demonstrated that using more sophisticated methods for PSD estimation, does not significantly improve the estimation of power spectrum central frequency measures (the mean or the median frequency) (Farina and Merletti, 2000).

Since the differences were insignificant, we decided to use the periodogram estimate in our study. As a measure of the central tendency of PSD we used the mean frequency of the PSD (MNF). MNF is defined as the first moment of the PSD and was developed by Kwatny et al., 1970. For the continuous spectrum spanning the frequencies between 0 and fMAX is defined as: (3) The value of MNF was calculated for each segment and used as a frequency parameter of muscle fatigue.

The values of MNF belonging to each muscle were plotted as a function of time. Since all data exhibited a trend toward linear decrease with time, a linear model was fitted to all MNF data sets to the initial MNF value (the value of MNF at the time of the first stroke) and the final MNF value (the value of MNF at the time of the last stroke) and labeled as MNF_beg and MNF_end, respectively. To normalize results between subjects, the final MNF values were expressed as a percentage of the initial values and labeled MNF_n as: (4) The slope of the regression line was also calculated as an estimate of the rate of change over time.

Blood samples were collected from the earlobe at 3, 5, 7 minutes after the bout. Post exercise blood lactate concentrations were measured using a hand-held portable apparatus (Lactate Pro Analyser™, Arkay, Kyoto, Japan). According Pyne et al., 2001 this equipment exhibits a high accuracy compared with other lactate analyzers and high degree of versatility under a variety of testing conditions.

The swims were videotaped simultaneously on sagittal plane with a pair of cameras providing a dual projection from both underwater (Sony Mini Dv DCR-HC42E, EUA) and above (Sony Mini Dv DCR-HC42E, EUA) the water surface. The cameras were placed stationary at 25-m of the headwall, on a lateral wall, perpendicular to the swimmer’s trajectory and 10-m away. The study comprised the kinematic analysis of stroke cycles (Ariel Performance Analysis System, Ariel Dynamics Inc., USA) at a sample rate of 50Hz. Zatsiorsky’s model with an adaptation by De Leva, 1996 was used with 8 anatomical landmarks and a trunk division in two articulated parts (Barbosa et al., 2010; Hirata and Duarte, 2007; Lafond et al., 2004; Pavol et al., 2002). To create a single image of dual projection as described previously (Barbosa et al., 2006; 2010), the independent digitalization from both cameras was reconstructed with the support of a calibration volume (16 point) and a 2D-DLT algorithm (Abdel-Aziz and Karara, 1971).

The stroke cycle was measured between the 18^th m and 22^th m, and identified at the end of the leg´s recovery (i.e., when the knees were in full flexion and feet on eversion).

Stroke mechanics were measured by the stroke cycle period (P, s), the stroke rate (SR = 1/P, Hz), the stroke length (SL, m) and the mean swimming speed of the full stroke (v, m·s^-1). Finally, the swimming efficiency was estimated using the stroke index (SI = v SL, m²·c^-1·s^-1) as suggested by Costill et al., 1985. This index assumes that, at a given v, the swimmer with greater SL has the most efficient swimming technique.

The assumptions of normality of data (Kolmogorov-Smirnov test) and homogeneity of variance (Levene’s test) were both confirmed for all parameters. All results are reported as the mean values along with standard deviation values. The repeated measures ANOVA, with subsequent Tukey’s test as post hoc analysis where applicable, were used for multiple comparisons between the laps (p < 0.05).

Figure 4 presents the mean and SD of the kinematics parameters (Swimming Speed, SL, SR, SI) in the 4^th lap of the 200m breaststroke. Swimming speed decreased from 1.38 ± 0.09 m·s^-1 in the 1^st lap to 1. 14 ± 0.08 m·s^-1 in the 4^th lap. Significant variations were observed on swimming speed between the 1^st lap and 2^nd lap, and between the 1^st lap and 3^rd lap [F = 21.27, p < 0.05].

SL decreased from 2.23 ± 0.18 m in the 1^st lap to 1.92 ± 0.15 m in 4^th lap. Significant variations were observed on the SL between the 1^st lap and 3^rd lap, and between the 1^st lap and 4^th lap [F = 4.41, p < 0.05].

SR decreased from 37.58± 4.90 stroke min^-1 in the 1^st lap to 34. 80±2.83 stroke min^-1 in 3^rd lap and then increased in the 4^th lap (35.91 ± 2.99 stroke·min^-1). No significant variations were observed between laps [F = 0.92, p > 0. 05]. Lastly, the SI was higher in the 1^st lap, and decreasing from the 1^st (3.07 ± 0.25 m²·s^-1) to the last lap (2.19 ± 0.29 m2·s^-1). Significant variations were observed on the SI between the 1^st lap and 2^nd lap, between the 1^st lap and 3^rd lap and between the 1^st lap and 4^th lap, [F = 16.94, p < 0.05].

The means and standard deviations in blood lactate concentration (La) measured after the 200 m breaststroke are shown in the Figure 5. The highest value (13.02 ± 1.72 mmol·l^-1) was observed at the 3^rd minute after the trial. Significant variations were observed in La between the 1^st minute and 3^rd minute, and between the 1^st minute and the 5th minute [F = 2. 47, p < 0.05].

The duration of the active phase (Table 1) was, for the BB muscle, higher in the beginning phase compared to the end phase of the swim for all four laps. The DA muscle presented a higher value in the beginning phase in the 1^st lap and 3^rd lap, but exhibited a higher value in the end phase for the 2^nd and 4^th lap. The PM muscle demonstrated a higher duration of active phase at the end for all laps with the exception of the 3^rd lap (PM_beg: 0.69 ± 0.62 s and PM_end: 0.55 ± 0.33 s) where the active phase was higher at the beginning than the end of the session. The TB muscle showed a higher active phase at the beginning in the 1^st and 3^rd lap, and higher values in the end in the 2^nd and 4^th lap.

The long duration of active phase was shown in BB at the beginning of the swim in the 1^st lap (0.99 ± 0.78 s) and the shorter for the TB at the end of swimming the 3^rd lap (0.52 ± 0. 15 s). Moreover, for all the muscles, the active phase showed the same trend in the 3^rd lap. The duration was higher in the beginning compared to the end of the swim. The differences between beginning and the end for each muscle were not statistically significant for p < 0.05.

The non-active phase in Table 2, showed an increase in the end of the test for all the muscles in the 3^rd lap and a decrease in the 4^th lap. The muscles BB, DA and PM demonstrated the same behavior in the relative duration of the non-active phase for each lap of the 200m breaststroke, that is, in the 1^st and 3^rd lap they had a longer duration at the end compared to the beginning of the swim, unlike in the 2^nd and 4^th lap where the higher duration was observed in the beginning of the swim.

The TB showed a different behavior compared to the other muscles studied, which presented a higher duration in the non-active phase in the beginning compared to the end of the swimming the 1^st and 4^th lap, and in the end for the 2^nd and 3^rd laps, moreover for all the muscles under study the non-active phase showed the same behavior in the 3^rd lap, where the duration was higher in the beginning which also happened in the active phase. The differences between the beginning and the end for each muscle were not statistically significant for p < 0.05.

The relative duration of activation (RAF) in Figure 6, decreases at the beginning phase in the 1^st and 3^rd lap for the DA and PM, whereas for the 2^nd and 4^th lap it was higher in the end for both the DA and PM. The BB demonstrated a RAF higher in the end phase in the 1^st and 2^nd lap, and the reverse in the 3^rd and 4^th lap. The TB decreased at the beginning phase over the first three laps, in contrast to the 4^th lap where the value was higher in the end. The differences between the beginning and the end for each muscle were not statistically significant for p < 0.05.

Figure 7 demonstrated that the average rectified value increases at the end phase compared to the beginning phase for all muscles comparing the 1^st lap with the 4^th lap, except for the DA (DA_beg:0.42 ± 0.13; DA_end: 0.42 ± 0.10) that demonstrated higher values in the beginning phase. The BB and TB muscles have similarities throughout the laps, that is, in the first three laps (1^st, 2^nd and 3^rd laps) the ARV was higher in the end compared to the beginning and in the 4^th lap the reverse was observed. The PM showed similarities in the 1^st and 3^rd lap which had higher values in the beginning, and in the 2^nd and 4^th lap had higher values at the end.

In Figure 8, the mean MNF value for the BB increased in all the laps at the end except the 3^rd lap (BB_end: 82.47 ± 16.20 Hz BB_beg: 86.15 ± 11.72 Hz). DA decreases in all the laps at the beginning except in the 3^rd lap (DA_end: 77.38 ± 13.68 Hz, DA_beg: 74.54 ± 15. 52 Hz). The PM decreases at the beginning in the 1^st lap (PM_beg: 89.65 ± 23.61 Hz; PM_end: 89.53 ± 17.74 Hz) and 3^rd lap (PM_beg: 79.37 ± 10.58 Hz; PM_end: 76.16 ± 16.92 Hz), and increases at the beginning in the 2^nd lap (PM_beg: 82.13 ± 15.33 Hz; PM_end: 86.11 ± 18.12 Hz) and 4^th lap (PM_beg: 74. 33 ± 14.61 Hz; PM_end: 75.18 ± 13.92 Hz). TB decreases at the beginning in the 1^st lap (TB_beg: 78.17 ± 10.83 Hz; TB_end: 74.79 ± 16.37 Hz) and in the 4^th lap (TB_beg: 72.45 ± 11.83 Hz; TB_end: 65.72 ± 12. 94 Hz), whereas in the 2^nd (TBend: 73.72 ± 13.68 Hz; TB_beg: 72.41 ± 12.65 Hz) and 3^rd lap (TB_end: 71.32 ± 11.26 Hz; TB_beg: 68.78 ± 11.03 Hz) we can observe an increase at the end phase compared to the beginning. Moreover, the MNF decreases at the end of swimming the 4^th lap relative to the 1^st lap for all muscles under observation for the 200m breaststroke. Note that the differences between the beginning and the end for each muscle were not statistically significant for p < 0.05.

The aims of this study were: i) to analyze four muscle activation patterns (duration of active and non-active phases) of the upper limbs during each lap for the 200m breaststroke, ii) quantify neuromuscular fatigue, through kinematics and physiologic assessment. The main findings were that some indicators of neuromuscular fatigue begin to appear shortly into the 2^nd lap, but are more evident in the 3^rd lap (i.e amplitudes increase and spectral variables decrease) in a non-linear fashion.

The average decrease in the swimming speed from the 1^st lap to the 4^th lap shows a much greater decline in speed compared to the results of similar studies (Takagi et al., 2004; Thompson et al., 2000; 2004). In breaststroke, an increase in swimming speed is associated with an increase in SR, but has a decline more in SL relative to the other swimming strokes (Thompson et al., 2004), supporting the changes in the swimming technique under the 200m breaststroke (Takagi et al., 2004). In this study SL suffered a decrease in the 1^st lap to the 4^th lap, whereas the SR decreased from the 1^st lap to the 3^rd lap and then increased in the 4^th lap.

As swimming speed decreased, different SR and SL combinations were observed, giving the best performance in the face of the fatigue task constraints, and decreasing the swimmer’s ability to maintain a constant SL. Increasing SR in the last 50m lap is due to the inability to generate sufficient extra effort to overcome high drag. The decrease in the swimming speed SL, SI and the increase of SR in the last lap suggested a trend for the inability to generate sufficient extra effort needed to overcome high drag throughout the test (Craig and Pendergast, 1979; Chatard et al., 2003; Thompson et al., 2004). This change in the stroke kinematics due to a decrease in the power output can be a consequence of peripherical fatigue phenomena (Ikuta et al., 2012; Stirn et al., 2011).

The blood lactate concentration levels after the test reached values that would be expected after a 200m breaststroke race and are in agreement with the results in other studies in which similar tests were performed with swimmers of the same or similar competitive level (e.g., Capelli et al., 1998; Lomax and Castle, 2011). These results suggest that the swimmers performed the swimming test very close to their maximum effort level. The average blood lactate concentration measured after 3 minutes was 13.2 ± 1.72 mmol·l^-1, representing the highest values. Increases of lactate concentration are responsible for fatigue by changes in the intracellular pH. As a result, muscle fiber conduction velocity decreases and changes the shape of the motor unit action potential waveform (Cifrek et al., 2009).

The amplitude of EMG signals is influenced by the number of active motor units (Moritani et al., 1986), their discharge rates, and the shape and propagation velocity of the intracellular action potentials (Dimitrova and Dimitrov, 2002).

The amplitude signal of EMG provided by the ARV demonstrated an increase and maintenance at the end phase compared to the beginning for BB and PM, comparing the 1^st lap with the 4^th lap, and a small decrease for DA and TB that demonstrated higher values in the beginning phase compared to the end phase. The swimmers, during the end of the swim, increase the period in the recovery phase, decreasing the amplitude of the EMG signal for DA and TB. These muscles are responsible for the movement of antepulsion of the arms up to the point where it is stretched in front.

EMG amplitude has been observed to increase during submaximal dynamic exercise (Tesch et al., 1990) and to decrease during exercises at maximal levels of voluntary contraction (Komi and Tesch, 1979). An increase of the EMG amplitude during swimming has already been observed in previous studies (Rouard and Clarys, 1995; Rouard et al., 1997; Stirn et al., 2011; Wakayoshi et al., 1994). Ikuta et al., 2012 noted that the mean amplitude value of the PM was significantly higher for the 4^th lap than for the 2^nd and 3^rd laps in 4x50 m front crawl. The authors noted that the increase activity of the PM may have been a compensatory strategy to maintain swimming speed during the test swim. Moreover, Nuber et al., 1986 verified that the PM continued to activate through recovery in the breaststroke.

So, it seems that a decrease in the active phase duration and the increase in the non-active phase on the 3^rd lap, with all the muscles monitored, reveal that in this lap muscles start to display the onset of neuromuscular fatigue indicating that this is a milestone of the 200m breaststroke race.

The MNF decreases at the 4^th lap in comparison to the 1^st lap for all muscles. Dynamic fatiguing tasks have been found to show a decrease in the mean power frequency (Tesch et al., 1990). A number of authors reported similar data in several other human muscles (Basmajian and De Luca, 1985; De Luca, 1984). The same is also apparent for swimming, or at least for front-crawl (Figueiredo et al., 2011; Stirn et al., 2011). The reduction of muscle fiber conduction velocity is one of the causes of signal power spectrum shift toward lower frequencies, and also of the increase in the EMG signal amplitude because of a spatial low-pass filtering effect of tissue as a volume conductor (De Luca, 1984).

However, during a fatiguing exercise, two oppositing effects might occur (Petrofsky and Lind, 1980): (i) a decrease in the mean frequency as discussed above and; (ii) increase in intramuscular temperature due to the exercise, which causes an increase in the mean frequency. So it is possible that during some types of exercise, the two effects compensate each other, and the decrements found in the mean frequency can be small or non-significant.

A decrease in swimming speed, SL, SI and an increase in the SR in the last lap, can be associated to the higher La after the test, to a non-significant increase of the amplitude and a non-significant decrease in the spectral variables. These suggest that the presence of fatigue in the upper limbs occurs in a different way for the different muscles studied, by a non-linear process. However, more investigation in order to better understand this phenomenon is required in the near future.

As blood flow is restricted, at a certain level of contraction, by intramuscular pressure, muscle becomes ischemic. Myoelectric manifestations of muscle fatigue might be affected by this event (Merletti et al., 1984). However, these are evidences for land-based exercises, and to the best of our knowledge none of this data has been collected during swimming.

The shape of the ARV during the four laps represent the increased activation patterns in the 2^nd and 3th lap in accordance with the normal adaptation and the tactic to use the remaining capacity in the final lap to achieve the best result. The higher activation by DA and TB showed that these muscles play a significant role in this event (Yoshizawa et al., 1976; Conceição et al., 2010).

With this research we can provide quick practical information for swimmers and coaches to adjust training methods.

This study showed a trend for a neuromuscular fatigue state in a non-linear fashion, with the 2^nd and 3^rd lap being the critical moments of a 200m breaststroke race. Decreases in swimming speed are related to SL, SI, and increases of the SR in the last lap, increases in the blood lactate concentration after the test, and a non-significant increase of the amplitude and non-significant decrease in the spectral parameters.

Both the technical and physiologic training in breaststroke swimmers should take into account the onset of fatigue in a non-linear fashion, including specific training to prevent sharp falls in the 2^nd and 3^rd laps of competitive racing.