Next Article in Journal
Development of Machine Learning Models for Prediction of Smoking Cessation Outcome
Next Article in Special Issue
Load Measures in Training/Match Monitoring in Soccer: A Systematic Review
Previous Article in Journal
A Cross-Country Examination on the Fear of COVID-19 and the Sense of Loneliness during the First Wave of COVID-19 Outbreak
Previous Article in Special Issue
Effects of Ballroom Dance on Physical Fitness and Reaction Time in Experienced Middle-Aged Adults of Both Genders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinematic Analysis of Water Polo Player in the Vertical Thrust Performance to Determine the Force-Velocity and Power-Velocity Relationships in Water: A Preliminary Study

1
Department of Medicine Systems, “Tor Vergata” University of Rome, via Montpellier 1, 00133 Rome, Italy
2
Sport Engineering Lab, Department of Industrial Engineering, “Tor Vergata” University of Rome, via del Politecnico 1, 00133 Rome, Italy
3
Department for Life Quality Studies, University of Bologna, 47900 Rimini, Italy
4
Robotics and Artificial Intelligence Lab, ENEA “Casaccia” Research Centre, via Anguillarese, 00301 Rome, Italy
5
Italian Swimming Federation, Stadio Olimpico Curva Nord, 00135 Rome, Italy
6
Department of Human Science and Promotion of Quality of Life, San Raffaele Open University of Rome, via di val Cannuta 247, 00166 Rome, Italy
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(5), 2587; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18052587
Submission received: 7 February 2021 / Revised: 24 February 2021 / Accepted: 27 February 2021 / Published: 5 March 2021
(This article belongs to the Collection Sports Medicine and Physical Fitness)

Abstract

:
Background: To date, studies on muscle force and power-velocity (F-v and P-v) relationships performed in water are absent. Aim: The goal of this study is to derive the F-v and P-v regression models of water polo players in water vertical thrust performance at increasing load. Methods: After use of a control object for direct linear transformation, displacement over the water and elapsed time was measured, by using a high-speed 2D-videoanalysis system, on 14 players involved in the study. Results: Intra-operator and player’s performance interclass correlation coefficient (ICC) reliability showed an excellent level of reproducibility for all kinematic and dynamic measurements considered in this study with a coefficient of variation (CV) of less than 4.5%. Results of this study have shown that an exponential force-velocity relationship seems to explain better the propulsive force exerted in the water in lifting increasing loads compared to the linear one, while the power and velocity have been shown to follow a second-order polynomial regression model. Conclusion: Given the accuracy of the video analysis, the high reliability and the specificity of the results, it is pointed out that video analysis can be a valid method to determine force-velocity and power-velocity curves in a specific environment to evaluate the neuromuscular profile of each water polo player.

1. Introduction

Water polo is characterized by a complex number of movements: swim with speed changes, faster counterattack actions, frequents changes from horizontal to vertical positions, shots, blocks and fight to gain or maintain the position in water [1]. Most of these actions (handlings, shots, fight) performed at high intensity require a vertical position in water [2]. There are two actions in movement of lower limbs of the water polo player that can be identified: the eggbeater kick (cyclic movement) [3,4] and the breaststroke kick (ballistic movement) [5]. The latter skill, involving maximal lower limbs muscle power, is usually adopted in trunk vertical thrust over the water level to complete the pass, in overall shots and in goalkeeper save actions. Indeed, some studies have found in elite female water polo players a significant correlation between the shot speed and the vertical thrust over the water level performed with breaststroke kick [6,7]. Net of sex differences, it is plausible to consider the vertical thrust one of the main skills also for men’s water polo.
From a biomechanical point of view, the maximal vertical thrust is obtained through the breaststroke kick techniques performed with quick movements on horizontal foot plan in extensive abduction, hip in flexion position and fast flexion-extension of knee [8]. Relative to muscular power and strength of lower limbs, it is common practice to test and condition the water polo players directly in the gravitational environment [1,9] without taking into account the specificity principle of neuromuscular and biomechanical performance that has to be transferred directly to the vertical thrust performed in water [10]. Relative to exercises performed on dry land, some authors showed a poor relationship between ground vertical jump and vertical thrust in water [11,12]. Recently, some authors, using different strength and power training methods performed on dry land or combined (dry-land and in-water) or in water only, showed a positive effect on some of the water polo skills performance with different results related to the method used [12,13,14]. Nevertheless, to date, studies on muscle force and power-velocity (F-v and P-v) relationships performed in water are absent [15]. In fact, the relationship between force and muscular contraction velocity has been determined in athletes to evaluate the dynamic neuromuscular characteristics in isotonic or ballistic conditions [15,16,17]. For this reason, individual power load-based training is difficult be carry out in the water taking into account that this is not specific if performed in a gravitational field.
Usually, in gravitational field, the most used devices for their practical applications in determining the mentioned above curves are linear encoders that, through a derivation process of measured space-time values, are able to calculate force and power parameters in relation to displaced mass [16,18]. Therefore, also taking into account the logistical difficulties in applying whatever isoinertial dynamometer in an aquatic environment, it remains mandatory to find a reliable and non-invasive assessment system. The practical goal of this study has been to verify an easy and reliable method, through a 2-D motion analysis approach, whose validity on kinematic measurements has already been shown [19], to assess the vertical displacement reached over the water level—net of the submerged breaststroke kick technique—and the related derivate kinematic as well as dynamic parameters. Furthermore, this needs to determine the accuracy of the measurement system together with the intra-rater and neuromuscular performance reliability of the assessment method used. In order to obtain in the aquatic environment F-v and P-v relationships like those obtained in gravitational field, it a test protocol was used at increasing loads performing the vertical thrust with a breaststroke technique.

2. Materials and Methods

2.1. Subjects

Fourteen male sub-elite level water polo players, (age 22.7 ± 5; Body Weight, 72.9 ± 8.2 kg; height 178.9 ± 5.2 cm, Body Mass Index 22.8 ± 2.2 kg/m2) participating in the regional championships (Serie C level) organized by the Italian Swimming Federation participated in this study. The body mass and height of the subjects were measured to the nearest 0.5 kg and 0.5 cm, respectively (Seca Beam Balance-Stadiometer, Germany). The players with physical problems (pain or injuries) or with low compliance training were excluded from the study. Written informed consent was obtained from participants (n = 14) before being tested. The research was approved by the Internal Research Board of “Tor Vergata” University of Rome. All procedures were carried out in accordance with the Declaration of Helsinki.

2.2. Experimental Design

This study requires the subjects to perform in the water vertical thrust tests at increasing load. This has been applied to the subjects by using the Water polo Overload Test/Training (WOT) [1] equipment, shown in Figure 1, that consists of a harness made of belts that are worn by the player and a load which can be fastened to its lower extremity, and does not interfere in any way with the legs’ movements.
The subjects, once they reached the position between the posts, spent a few seconds floating with the eggbeater kick technique to achieve and to maintain the optimal start position keeping the acromion at the same water level as before to perform, by using the breaststroke kick technique, an explosive boost to raise vertically the body as high as possible. In addition, to avoid any coordinative influence, the subjects held their upper limbs to their shoulders during the performance. Then, wearing the WOT system, they started to perform the increasing load test starting from free load condition which represents the reference trial for the test specificity. The player performed vertical thrusts increasing the load by 5 kg at each step (5, 10, 15, 20, 25 kg) where the 25 kg was the maximum load vertically raised at the limit of the buoyancy. The best trial of three measurements in terms of displacement and verticality at each increased load performance was selected for statistical analysis. Each subject completed raised load test with almost three-minute rest time between the trials enough to recovery from single boost performance.
In order, to determine day-to-day reliability, the subjects underwent the same protocol after two rest days. The measurements were performed in the same swimming pool where the subjects usually train with the water temperature of 29 °C, pH of 7.2–7.6, and an environmental temperature that ranges between 24–26 °C at a humidity of 75%. These parameters remained unchanged in both the test days. One week before the test administration, the subjects performed some simulations for familiarization with the equipment. On the first test day, the anthropometric data were recorded and the subjects performed, after a warm-up, an incremental loads protocol test with the WOT. The warm-up exercises were completed in 15 min and consisted general to specific skills performed with a progressive increase intensity.
The subject, to maintain the assigned start position must produce, properly moving his limbs (eggbeater kick), an upward floating force (Fek) equal to the difference between the sum of his body weight (WB) plus the weight of the eventual additional load (WL) and the sum of the respective buoyant forces (FbB and FbL) [20]. Equation (1) takes into account a further force (Fdw) related to the power that is wasted to accelerate the water downwards and that does not contribute to the thrust.
F e k F d w = ( W B F b B ) + ( W L F b L )
Considering that the body weight and the buoyant force, respectively (VB [m3] is the volume of the body, ρs its density [kg/m3], ρw the water density (995.96 kg/m3 at 29 °C), g [m/s2] the acceleration due to the gravity and f the fraction of the submerged body), are W B = g V B ρ B and F b B = g f V B ρ w . The Equation (1) can be written as follow
F e k = W B ( 1 f 0 ρ w ρ B ) + ( W L F b L ) + F d w
where f0 is the fraction of the submerged body at starting position (i.e., the volume of the whole body without the head and neck).
Conversely, when the subject has to perform the vertical thrust, he has to provide, by moving his legs with a breaststroke kick, an upward force (Fbk) greater than the sum of the weight force of both load WL and body WB, the friction forces Ffr, the buoyant forces Fb and the losses Fdw as reported in the follow expression:
F b k > ( W B F b B + F f r B ) + ( W L F b L + F f r L ) + F d w
where FbB and FbL represent, respectively, the buoyant force of the subject and the load, while FfrB and FfrL are the respective friction forces [20].
The buoyant and friction forces of the load have been evaluated starting from manufacturing features (material, shape and dimensions). Moreover, because the load remains entirely immersed during the whole test, the relative buoyant force always presents the same value. Furthermore, since the friction depends on the displacement velocity, it will be possible, due to the features of the WOT, to use the same value of the vertical thrust velocity computed for the subject.
A different approach is required for the calculation of the same forces for the human body. Indeed, the buoyant force depends on the fraction of the volume of the immersed body that tends to vary during each test because the height of the vertical thrust changes at different loads. Therefore, the accuracy on the computation of such a term depends on a proper estimation of both the volumes of the different parts of the body and its density. In this study, this value of density ρB has been simply estimated, for each subject, starting from his weight and height by using the procedure suggested in [21,22]. Moreover, in order to identify the right fraction of the submerged body volume, we used the mean relative percentage values of the volumes of the different segments of the human body [23,24]. Finally, we chose to evaluate the upward force performed with the breaststroke kick considering the buoyant force (FbB0) at the start position only (i.e., when it presents its maximum level) and where the estimation of the immersed body volume shows the minimum error. Consequently, the calculation of the upward force will be underestimated, in the same way, for all the subjects.
In this context, to avoid the difficulties in evaluating the body volumes at different vertical thrust height, we consider this force minus the relative body weight (i.e., F b k = F b k W B ). Therefore, neglecting the skin friction drag of the body FfrB and the losses Fdw, the breaststroke force equation to lift the loads becomes
F b k = m L ( a + g ) F b B 0 F b L + F f r L
where a is value of the acceleration in the thrust and mL the mass of the load. Thus, the corresponding mechanical power P b k relative to the force exerted with the breaststroke kick during the vertical thrust can be computed as:
P b k = F b k · v

2.3. Experimental Procedure

Each trial was recorded at 240 fps (time resolution ~4 ms) with a high-speed camera (Casio Exilim EX-ZR 3700—Japan) that was positioned at a distance of 2.30 m perpendicular at the sagittal plan of the subject in water. To verify the verticality of upper body displacement over the water level, a second camera was placed orthogonally (and at the same distance) to the first one so that the subject lay in the center of view angles of both cameras. No subject that performed a jump too far from his vertical was considered in this study.
The video analysis procedure allows, by processing the acquired videos, the value of the displacement Δd of the vertical jump to be obtained and the time Δt required by the subject to reach the maximum elevation. In detail, the duration of the rising phase of each thrust was obtained by multiplying the frame time by the number of frames between the start of the movement (i.e., the frame where is observed the starting vertical movement) from the buoyance position and the point where the subject reaches the higher position. The starting position was identified where the subject stands stably with the acromion at the water level.
Moreover, the height of each thrust has been evaluated by measuring (in number of pixels) the distance between the position, in the two different frames of start and top position, of the marker placed on the center of the subject’s headgear with respect to the level of the water (Figure 1). The values of mean velocity, force and power were calculated starting from these values while the muscular force and the relative power produced by the subject were computed starting from the maximum displacement reached in the jump by using Equations (4) and (5) respectively. A single operator provided the acquisition of these values, by using specific tools available inside the video analysis software BioMovie ERGO© (by Infolabmedia, Italy).

2.4. Video Analysis System Accuracy

The size of the images obtained by the camera was 432 × 320 pixels. The calibration factor KC [pix/cm] has been evaluated by using a 2D-DLT (2D-direct linear transformation) [25] with vertical (post) and horizontal (crossbar) reference objects in the picture placed at the same distance of the subject (i.e., the subject and reference object are in the same calibration plane). Considering as negligible the horizontal displacement of the athletes during vertical jump performance, the post height (86 cm) only was considered for the calculation of the factor Kc that has been evaluated as 0.717 cm/pix.
Moreover, the relative errors (in percentages) of measured displacement ( ε d % ) and of measured time ( ε t % ) can be evaluated as:
ε d % =   ε u   K C d 0 kg · 100  
where d0kg is the average of the height reached by the subjects during the trials at free load and ε u is the uncertainty error, due to the motion blur [26] in the estimation of the maximum reference point in the vertical displacement.
The time absolute error of the camera can be assumed as equal to a frame time ( ε t = 4 ms) with a negligible jitter considering that the inaccuracy of the internal clock oscillator of the camera can be estimated at less than 0.1 μs. The percentage errors for the forces as well as the power were evaluated according to the usual methodologies for the error propagation [27].
The estimation of the error for the buoyancy force relative to the different subjects was computed to take into account a value equal to 3.5% for the percentage error in the measure of the body volumes ( ε V B % ) as reported by [22].

2.5. Statistical Analysis

Data in text, tables and figures are expressed as mean ± standard deviation (SD). The Kolmogorov-Smirnov test was used to validate the assumption of normality. Since no significant deviations from normality were detected, the coefficient of variation (CV), interclass correlation coefficients (ICC), standard error of measurement (SEM) and 95% confidence interval (95% CI) were calculated to determine the day-to-day reliability for displacement, time, velocity, acceleration, force and power. Moreover, the ICC was used as assessment test of consistency, repeatability of quantitative measurements made by same operator and to evaluate the athlete’s performance in two different days. Paired t-tests with Bonferroni adjustment and the Pearson correlation coefficient (r) were used for between-group comparisons, for test-re-test measurements repeatability and to determine the level of specificity among selected variables of the test. In addition, the effect sizes (ES) were also calculated using Cohen’s d between the pre-test and post-test means [28], where small effect was 0.1, moderate 0.3 and large was 0.5 [29]. The level of statistical significance was set at p < 0.05. The IBM-SPSS 20.0 (SPSS, Inc., Chicago, IL, USA) was used for statistical analysis.

3. Results

3.1. System Accuracy

In order to evaluate the displacement relative error we apply Equation (2) where d(0kg) is equal to 68 cm and the uncertainty error ε u set to 3 pixels, the ε d % is equal to 3.11% while the percentage errors for the velocity, acceleration, force and power can be estimated as ε v % = 3.37 % , ε a % = 4.25 % , ε F b k % = 6.89 % , ε P b k % = 7.58 % respectively.

3.2. Reliability

Test-retest values of Mean, SD, SEM, ICC, Pearson correlation coefficient (r) and the CV relative to the displacement, velocity, acceleration, force and power performed in the same day and day-to-day are reported in Table 1. The average displacement decreases from 0.69 m (without load) to 0.15 m (load 25 kg), with r ranging from 0.87 at 5 kg to 0.99 at 0, 10 and 25 kg respectively. The thrust performance time (s) decreases at increasing loads with r ranging from 0.86 at 15 kg to 0.99 at 25 kg. Also, the vertical velocity (m/s) decreases as the loads increase with r ranging from 0.95 at 5 kg to 0.99 at 25 kg. Also, the acceleration (m/s2) decreases at increase load with r ranging from 0.93 at 20 kg to 0.99 at 5 kg with high correlation values. By contrast with the kinematic parameters, the force increases proportionally to the load ranging from 20.31 N at 5 kg to 304.35 N at 25 kg with high r values ranging from 0.95 at 20 kg to 0.99 at 5 kg. The power increases progressively from 5 kg (r = 0.99) to reach its maximal value at 20 kg (442.70 W with r = 0.94) and then decreases at 25 kg (313.80 W). The ICC of all parameters, expressed in detail in Table 1, showed an excellent level of reproducibility for all measurements. The CV, while remaining low in the kinematic parameters, tends to increase in the dynamic ones reaching its maximum value of 4.32 at the Pbk 10 kg. In addition, the SEM values observed in day-to-day trials are very low in all kinematic parameters considered for each load. The effect size (ES) calculated between pre-test and post-test means, showed a magnitude ranging from small to moderate in all kinematic and dynamic observed parameters (Table 1). The level of statistical significance was set at p < 0.05. An IBM-SPSS 20.0 (SPSS, Inc., Chicago, IL, USA) was used for statistical analysis.

3.3. Specificity

For the specificity of the method analyzed in this study, the vertical thrust without overloads was considered as a specific water polo skill and, therefore, correlated with the same skill performed at increasing loads. The analysis of correlation between displacement and the force, power and velocity at each load showed a strong correlation with low load (until 20 kg). As the loads increase, these correlations tend to decrease until it becomes not significant at 20 and 25 kg for velocity while for the force and power became non-significant at 25 kg only (Figure 2).

3.4. Force-Velocity and Power-Velocity Relationships

Taking in account the means and SD values of force, power and velocity obtained by the measurements showed in Table 1, it was possible to determine a linear relationship between force and velocity and a quadratic curve between power and velocity (Figure 3).
It is worth noting that force and velocity values presented an inverse trend at increasing loads while the power reached the minimum value at 5 kg condition, reached a higher value at 20 kg, and then decreased again at 25 kg load. Both curves, depicted in Figure 2, show the linear and quadratic equation with a high correlation value (r = 0.92 and 0.99 for F-v and P-v curves respectively). Moreover, with a more accurate analysis of the F-v curve, it is interesting to highlight that the values recorded up to 20 kg maintain a linear relationship while at 25 kg the curve tends to assume an exponential like shape (Figure 4).
Therefore, the following exponential equation (Equation (7)) seems to fit better the behavior of the F-v relationship of the incremental loads test (r = 0.99; p < 0.001) than the previous linear one (r = 0.92):
F b k ( v ) = F 0   e   1 2 ( v a ) b
where F0 is the maximum value of the force recorded at the lowest value of velocity (constant a) of the vertical thrust performed in the test, v is the velocity value recorded at each load while the constant b allows us to model the growth in the exponential rate.

4. Discussion

The results of this study confirm the accuracy of the kinematic parameters measured with the video analysis system. Displacement, time and the calculated parameters as velocity and acceleration showed error values contained below 4.5% in any ballistic performance (breaststroke kick) load conditions, while the dynamic derivate as force and power showed the maximum error below 8%. It needs to be underlined that, for each parameter (measured or calculated), the relative error was less than the mean differences observed among athletes in each load condition performance.
The level of reproducibility of all parameters assessed in this study was very high between the two trials performed in two different days (Table 1) in terms of correlation (r from 0.86 to 0.99) and CV (<4.5%). Thus, the intra-rater reliability on the video analysis system used in this study and the water polo player’s performance provide a consistent result, with an excellent level of ICC, satisfying the basic requirement of any assessment method [16].
Usually, the methods used to determine the F-v and P-v curves of leg extensor muscles are the half-squat weightlifting or jumping test performed at increasing load in gravitational environment. Considering that the specificity represents the most important discriminant criterion of a test [30], it should be emphasized there is scant specificity from biomechanical and neuromuscular points of view between dry half squat, vertical jump and vertical thrust on the water performance [11]. The strong significant correlation showed in this study between the free load vertical displacement and the other kinematic and dynamic parameters obtained at increasing loads (Figure 2), gave to this method a high level of specificity from biomechanical and neuromuscular points of view. Biomechanically, lifting the upper body over the water level means apply a lift force able to counteract the drag force. Indeed, by using the breaststroke kick technique, Sanders [8] showed that the lift forces in the water polo boots are developed through the synergic action of feet where their velocity action is obtained using the anteroposterior and mediolateral directions, followed by the knee extension and trunk straightening from their start angle with respect to the horizontal plane. In this context, squat weightlifting or dry-land jump involves the neuromuscular system in a different biomechanics condition [31]. In addition, Platanou [11] observed no correlation between the vertical thrust on water and the vertical jump on dry land (r = 0.25). From a neuromuscular point of view, in this study the relationship between vertical thrust tends to decrease at increasing load just to become minimal in correspondence of the maximal strength (25 kg) (Figure 3). Furthermore, in the water the muscle contraction does not use the same strategies related to the stretching-shortening cycle and the performance is not characterized neither by the use of elastic energy nor by stretch reflex, typical features of natural gravity movements on developing the ground reaction force. Currently, the methods used to assess the power and strength of the leg extensor and arm muscles are performed in a gravity environment [32]. This study represents the first tentative, in aquatic environment, able to determine the linear F-v and parabolic P-v relationship of lower limb muscles during a vertical thrust performance directly on the water. Both curves maintain the same characteristics of the F-v and P-v relationship observed on an athlete’s performance made in a gravitational environment using a leg or arm extensors isotonic [16,18] or isokinetic devices [33] or ballistic movement [34]. According to Jaric [32], the linear relationship of F-v and consequent parabolic P-v relationship performed in a multi-joint performance showed a strong correlation revealing a high reliability of all the parameters considered in this study as reported in Table 1. Moreover, the second-order polynomial regression of P-v has shown a Pmax in correspondence with 20 kg that represents the optimal load averagely expressed by the analyzed subjects. In this context, also the high values of specificity observed with low loads tend to decrease becoming not significant after the Pmax load, probably due to a different neuromuscular pattern. In fact, according to the motor unit size recruitment principle of Henneman [35], by using an increasing loads protocol, the water polo players exhibited a decreasing heights and muscle contraction velocity on vertical thrust in relation to increasing muscular strength (increased loads) (Figure 3). Although the force and velocity values satisfy the linearity of the relationship, it is worth noting that, as shown in Figure 4, these values recorded at 25 kg tend to lose this linearity. Indeed, it can be presumed that the force exerted by the lower limbs in holding and lifting very heavy loads reached a saturation level (plateau) without ever reaching the maximum isometric force, which is impossible to obtain in a water environment, as instead it is observed in a gravitational field. In this condition, it seems conceivable to consider that the F-v curve could switch its linearity in an exponential like relationship with heavy loads (i.e., when the vertical displacement will become negligible to the further increasing load without muscular force increase). In this case, the eggbeater and breaststroke kick are performed alternatively to maintain buoyancy, as happens in games during the hard attacks and blocks between centre forward and defender. Indeed, it is feasible to assume that the limiting factor of the maximal force exerted by a breaststroke kick is based upon the maximal buoyant force sustained with eggbeater kick performance [36]. In this context, the exponential model represented by Equation (7), with the strongest relationship (r = 0.99) compared with the previous linear curve (r = 0.92; p < 0.001), seems to better explain the development of the propulsive force exerted in the water by water polo players with the breaststroke kick technique [8,11]. Furthermore, the P-v curve (Figure 3) is not influenced by the linear or exponential F-v curve maintaining the same parabolic trend. The scant correlation observed between the maximal strength and boots performance shows that Pmax load could be considered as a reference load to plan strength or velocity conditioning training in water polo players.
Then, in accordance with the incremental load method for eggbeater kick used by Melchiorri [1] the method for breaststroke kick used in this study, seems to be able to overcome the specificity limits of all strength and power monitoring and training method performed on the land for water polo players.
As a limitation, the use of a more performing camera or a new sensor system able to detect the space-time variations of players on the water, should allow the improvement of such a measure with a consequent reduction of the error. In addition, the indirect assessment of the volumes of the different sections of the human body could lead to a less accurate estimation than the real buoyant force of the body players at the different heights reached during the vertical thrust.

5. Conclusions

Considering the accuracy and the reliability recorded between two consecutive trials and two days’ video analysis measurements and the high specificity of the breaststroke kick performed at increasing loads, it is reasonable to consider the validity of this method. Thus, the kinematic assessment of the water polo player performing specific neuromuscular and biomechanical patterns in specific environments could reduce the bias in the assessment and training. In line with these considerations, this easy and practical method could provide coaches and trainers specific indications of the individual linear (especially at light loads) or exponential F-v and quadratic P-v relationships of each water polo player, useful for strength and power monitoring and conditioning to perform directly on the water without time spent in a non-specific regimen. Future studies should be required to verify these preliminary results calculating more accurately the human volumes and densities. Moreover, the same method should be also verified on female water polo players.

Author Contributions

G.A. and V.B. contributed to Conceptualization and Methodology. C.R. and V.V. did the data collection. A.Z., C.R. and G.M. analyzed the data; G.A. and E.P. helped in the literature review. V.B. and C.R. involved in manuscript revision. G.A., V.B., C.R. and E.P. wrote the manuscript. G.A. and V.B. supervised all the phases of the study. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the internal Research Board of the “Tor Vergata” University of Rome.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All study data are included in the present manuscript.

Acknowledgments

We want to acknowledge the contribution of Cristiano Maria Verrelli, Niloofar Lamouchideli and Fabrizio Tufi.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Melchiorri, G.; Viero, V.; Triossi, T.; Tancredi, V.; Galvani, C.; Bonifazi, M. Testing and Training of the Eggbeater Kick Movement in Water Polo: Applicability of a New Method. J. Strength Cond. Res. 2015, 29, 2758–2764. [Google Scholar] [CrossRef] [PubMed]
  2. Melchiorri, G.; Castagna, C.; Sorge, R.; Bonifazi, M. Game Activity and Blood Lactate in Men’s Elite Water-Polo Players. J. Strength Cond. Res. 2010, 24, 2647–2651. [Google Scholar] [CrossRef]
  3. Homma, M.; Homma, M. Swimming: Coaching Points for the Technique of the Eggbeater Kick in Synchronized Swimming Based on Three-Dimensional Motion Analysis. Sports Biomech. 2005, 4, 73–87. [Google Scholar] [CrossRef]
  4. Sanders, R.H. Analysis of the Eggbeater Kick Used to Maintain Height in Water Polo. J. Appl. Biomech. 1999, 15, 284–291. [Google Scholar] [CrossRef]
  5. Tsunokawa, T.; Nakashima, M.; Takagi, H. Use of Pressure Distribution Analysis to Estimate Fluid Forces around a Foot during Breaststroke Kicking. Sports Eng. 2015, 18, 149–156. [Google Scholar] [CrossRef]
  6. McCluskey, L.; Lynskey, S.; Leung, C.K.; Woodhouse, D.; Briffa, K.; Hopper, D. Throwing Velocity and Jump Height in Female Water Polo Players: Performance Predictors. J. Sci. Med. Sport 2010, 13, 236–240. [Google Scholar] [CrossRef] [PubMed]
  7. Platanou, T.; Varamenti, E. Relationships between Anthropometric and Physiological Characteristics with Throwing Velocity and on Water Jump of Female Water Polo Players. J. Sports Med. Phys. Fit. 2011, 51, 185. [Google Scholar]
  8. Sanders, R.H. A Model of Kinematic Variables Determining Height Achieved in Water Polo Boosts. J. Appl. Biomech. 1999, 15, 270–283. [Google Scholar] [CrossRef]
  9. Veliz, R.R.; Requena, B.; Suarez-Arrones, L.; Newton, R.U.; de Villarreal, E.S. Effects of 18-Week in-Season Heavy-Resistance and Power Training on Throwing Velocity, Strength, Jumping, and Maximal Sprint Swim Performance of Elite Male Water Polo Players. J. Strength Cond. Res. 2014, 28, 1007–1014. [Google Scholar] [CrossRef]
  10. Gobbi, M.; D’ercole, C.; D’ercole, A.; Gobbi, F. The Components of the Jumps in Expert and Intermediate Water Polo Players. J. Strength Cond. Res. 2013, 27, 2685–2689. [Google Scholar] [CrossRef]
  11. Platanou, T. On-Water and Dryland Vertical Jump in Water Polo Players. J. Sports Med. Phys. Fit. 2005, 45, 26–31. [Google Scholar]
  12. Martin, M.S.; Blanco, F.P.; Villarreal, E.S.D. Effects of Different In-Season Strength Training Methods on Strength Gains and Water Polo Performance. Int. J. Sports Physiol. Perform. 2021, 1, 1–10. [Google Scholar] [CrossRef]
  13. De Villarreal, E.S.; Suarez-Arrones, L.; Requena, B.; Haff, G.G.; Ramos-Veliz, R. Effects of Dry-Land Vs. In-Water Specific Strength Training on Professional Male Water Polo Players’ Performance. J. Strength Cond. Res. 2014, 28, 3179–3187. [Google Scholar] [CrossRef]
  14. Sáez de Villarreal, E.; Suarez-Arrones, L.; Requena, B.; Haff, G.G.; Ramos Veliz, R. Enhancing Performance in Professional Water Polo Players: Dryland Training, in-Water Training, and Combined Training. J. Strength Cond. Res. 2015, 29, 1089–1097. [Google Scholar] [CrossRef]
  15. Kaneko, M. The Relation between Force, Velocity and Mechanical Power in Human Muscle. Taiikugaku Kenkyu 1970, 14, 143–147. [Google Scholar] [CrossRef] [Green Version]
  16. Bosco, C.; Belli, A.; Astrua, M.; Tihanyi, J.; Pozzo, R.; Kellis, S.; Tsarpela, O.; Foti, C.; Manno, R.; Tranquilli, C. A Dynamometer for Evaluation of Dynamic Muscle Work. Eur. J. Appl. Physiol. Occup. Physiol. 1995, 70, 379–386. [Google Scholar] [CrossRef] [PubMed]
  17. Cavagna, G.A.; Dusman, B.; Margaria, R. Positive Work Done by a Previously Stretched Muscle. J. Appl. Physiol. 1968, 24, 21–32. [Google Scholar] [CrossRef] [PubMed]
  18. Sreckovic, S.; Cuk, I.; Djuric, S.; Nedeljkovic, A.; Mirkov, D.; Jaric, S. Evaluation of Force-Velocity and Power-Velocity Relationship of Arm Muscles. Eur. J. Appl. Physiol. 2015, 115, 1779–1787. [Google Scholar] [CrossRef] [PubMed]
  19. Zult, T.; Allsop, J.; Tabernero, J.; Pardhan, S. A Low-Cost 2-D Video System Can Accurately and Reliably Assess Adaptive Gait Kinematics in Healthy and Low Vision Subjects. Sci. Rep. 2019, 9, 1–11. [Google Scholar] [CrossRef] [Green Version]
  20. Davidovits, P. Physics in Biology and Medicine; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  21. Haycock, G.B.; Schwartz, G.J.; Wisotsky, D.H. Geometric Method for Measuring Body Surface Area: A Height-Weight Formula Validated in Infants, Children, and Adults. J. Pediatrics 1978, 93, 62–66. [Google Scholar] [CrossRef]
  22. Sendroy Jr, J.; Collison, H.A. Determination of Human Body Volume from Height and Weight. J. Appl. Physiol. 1966, 21, 167–172. [Google Scholar] [CrossRef] [Green Version]
  23. Daniell, N.; Olds, T.; Tomkinson, G. Volumetric Differences in Body Shape among Adults with Differing Body Mass Index Values: An Analysis Using Three-Dimensional Body Scans. Am. J. Hum. Biol. 2014, 26, 156–163. [Google Scholar] [CrossRef]
  24. Var, S.M.; Marangoz, I. Leg Volume and Mass Scales of Elite Male and Female Athletes in Some Olympic Sports. World J. Educ. 2018, 8, 54–58. [Google Scholar]
  25. Brewin, M.A.; Kerwin, D.G. Accuracy of Scaling and DLT Reconstruction Techniques for Planar Motion Analyses. J. Appl. Biomech. 2003, 19, 79–88. [Google Scholar] [CrossRef]
  26. Tiwari, S.; Shukla, V.; Singh, A.; Biradar, S. Review of Motion Blur Estimation Techniques. J. Image Graph. 2013, 1, 176–184. [Google Scholar] [CrossRef]
  27. Rabinovich, S.G. Measurement Errors and Uncertainties: Theory and Practice; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006; ISBN 978-0-387-29143-7. [Google Scholar]
  28. Cohen, J. Statistical Power Analysis for the Behavioral Sciences; Academic Press: Cambridge, MA, USA, 2013; ISBN 978-1-4832-7648-9. [Google Scholar]
  29. Cooper, H.; Hedges, L.V. The Handbook of Research Synthesis; Russell Sage Foundation: New York, NY, USA, 1993; ISBN 978-1-61044-137-7. [Google Scholar]
  30. Sale, D.; MacDougall, D. Specificity in Strength Training: A Review for the Coach and Athlete. Can. J. Appl. Sport Sci. 1981, 6, 87–92. [Google Scholar] [PubMed]
  31. Jacobs, R.; Bobbert, M.F.; van Ingen Schenau, G.J. Mechanical Output from Individual Muscles during Explosive Leg Extensions: The Role of Biarticular Muscles. J. Biomech. 1996, 29, 513–523. [Google Scholar] [CrossRef]
  32. Jaric, S. Force-Velocity Relationship of Muscles Performing Multi-Joint Maximum Performance Tasks. Int. J. Sports Med. 2015, 36, 699–704. [Google Scholar] [CrossRef]
  33. Perrine, J.J.; Edgerton, V.R. Muscle Force-Velocity and Power-Velocity Relationships under Isokinetic Loading. Med. Sci. Sports 1978, 10, 159–166. [Google Scholar] [CrossRef]
  34. Bosco, C.; Mognoni, P.; Luhtanen, P. Relationship between Isokinetic Performance and Ballistic Movement. Eur. J. Appl. Physiol. Occup. Physiol. 1983, 51, 357–364. [Google Scholar] [CrossRef] [PubMed]
  35. Henneman, E.; Somjen, G.; Carpenter, D.O. Functional Significance of Cell Size in Spinal Motoneurons. J. Neurophysiol. 1965, 28, 560–580. [Google Scholar] [CrossRef] [PubMed]
  36. Kawai, E.; Tsunokawa, T.; Sakaue, H.; Takagi, H. Propulsive Forces on Water Polo Players’ Feet from Eggbeater Kicking Estimated by Pressure Distribution Analysis. Sports Biomech. 2020, 1–15. [Google Scholar] [CrossRef]
Figure 1. Frontal view of water polo overload test (WOT) conditions and acting forces. Buoyancy force of the subject (FbB), Body weight of the subject (WB), eggbeater kick force (Fek), buoyancy force of the load (FbL), weight of the load (WL), force relative to power that is wasted to accelerate water downwards (Fdw).
Figure 1. Frontal view of water polo overload test (WOT) conditions and acting forces. Buoyancy force of the subject (FbB), Body weight of the subject (WB), eggbeater kick force (Fek), buoyancy force of the load (FbL), weight of the load (WL), force relative to power that is wasted to accelerate water downwards (Fdw).
Ijerph 18 02587 g001
Figure 2. Median correlation coefficients and their ranges obtained comparing the vertical thrust’s height free load with individual Fbk, Pbk and v at 5–25 kg.
Figure 2. Median correlation coefficients and their ranges obtained comparing the vertical thrust’s height free load with individual Fbk, Pbk and v at 5–25 kg.
Ijerph 18 02587 g002
Figure 3. Linear F-v (grey line and squares) and second-order polynomial regression P-v (black line and dots) with the relative regression equations building on vertical thrust performed at incremental loads (from 5 to 25 kg). Both curves are depicted according the average and SD of velocity, force and power at different load as shown in Table 1.
Figure 3. Linear F-v (grey line and squares) and second-order polynomial regression P-v (black line and dots) with the relative regression equations building on vertical thrust performed at incremental loads (from 5 to 25 kg). Both curves are depicted according the average and SD of velocity, force and power at different load as shown in Table 1.
Ijerph 18 02587 g003
Figure 4. Exponential F-v with the relative regression equations building on vertical thrust performed at incremental loads (from 5 to 25 kg). The curve is depicted according to the average and SD of velocity and force with the different loads as shown in Table 1.
Figure 4. Exponential F-v with the relative regression equations building on vertical thrust performed at incremental loads (from 5 to 25 kg). The curve is depicted according to the average and SD of velocity and force with the different loads as shown in Table 1.
Ijerph 18 02587 g004
Table 1. Day-to-day repeatability of average displacement (m), time (s), velocity (m/s) and acceleration (m/s2), force and power ± Standard Deviation (SD) of the vertical thrust exercise (in water) performed by 14 water polo players with increasing loads, r Pearson correlation coefficient; CV, Coefficient of Variation for repeated measurements; ICC, Interclass Correlation Coefficient; 95% Confidence Interval (CI); SEM, Standard Error of Measurement; and ES, Effect Size, for each load.
Table 1. Day-to-day repeatability of average displacement (m), time (s), velocity (m/s) and acceleration (m/s2), force and power ± Standard Deviation (SD) of the vertical thrust exercise (in water) performed by 14 water polo players with increasing loads, r Pearson correlation coefficient; CV, Coefficient of Variation for repeated measurements; ICC, Interclass Correlation Coefficient; 95% Confidence Interval (CI); SEM, Standard Error of Measurement; and ES, Effect Size, for each load.
Different DayDay 1Day 2rCVICC95% CISEMES
ParametersMean ± SDMean ± SD
Displacement (m)
D0 kg0.69 ± 0.020.69 ± 0.010.990.310.990.941 to 0.9980.00130.004
D5 kg0.55 ± 0.040.57 ± 0.050.871.610.850.315 to 0.9780.01030.147
D10 kg0.49 ± 0.050.50 ± 0.050.991.090.980.910 to 0.9980.0031−0.093
D15 kg0.44 ± 0.040.43 ± 0.040.942.430.940.718 to 0.9920.00610.091
D20 kg0.35 ± 0.030.34 ± 0.040.971.720.950.595 to 0.9930.00340.219
D25 kg0.15 ± 0.050.15 ± 0.050.990.940.990.983 to 0.9990.00130.004
Time (s)
T0 kg0.258 ± 0.0090.257 ± 0.0100.950.890.940.713 to 0.9920.00130.134
T5 kg0.244 ± 0.0080.245 ± 0.0100.911.280.940.648 to 0.9910.0018−0.173
T10 kg0.235 ± 0.0090.237 ± 0.0130.901.820.930.561 to 0.9900.0025−0.191
T15 kg0.225 ± 0.0100.225 ± 0.0110.861.780.930.494 to 0.9910.00230.017
T20 kg0.205 ± 0.0080.206 ± 0.0110.961.310.930.627 to 0.9910.0015−0.015
T25 kg0.141 ± 0.0230.143 ± 0.0230.990.730.990.878 to 0.9990.0008−0.075
Velocity (m/s)
V0 kg2.68 ± 0.172.69 ± 0.170.980.900.980.891 to 0.9970.0013−0.062
V5 kg2.30 ± 0.202.30 ± 0.230.952.130.950.733 to 0.9940.0283−0.023
V10 kg2.11 ± 0.232.12 ± 0.250.981.410.980.916 to 0.9980.0173−0.022
V15 kg1.95 ± 0.181.94 ± 0.210.971.980.960.808 to 0.9950.02220.066
V20 kg1.71 ± 0.141.65 ± 0.120.961.630.930.403 to 0.9900.01590.279
V25 kg1.03 ± 0.171.02 ± 0.170.991.380.980.903 to 0.9990.01220.073
Acceleration (m/s2)
Acc0 kg10.41 ± 1.0410.51 ± 1.110.971.720.970.843 to 0.9960.1041−0.094
Acc5 kg9.46 ± 1.099.43 ± 1.230.991.530.980.912 to 0.9980.08390.025
Acc10 kg8.99 ± 1.218.96 ± 1.400.963.110.960.752 to 0.9940.16160.024
Acc15 kg8.69 ± 0.958.66 ± 0.950.972.940.950.712 to 0.9930.14720.031
Acc20 kg8.34 ± 0.788.15 ± 0.790.932.370.920.552 to 0.9880.11320.241
Acc25 kg7.33 ± 0.467.15 ± 0.590.941.820.890.234 to 0.9920.10610.323
Force Bk (N)
F5 kg20.31 ± 11.3919.63 ± 11.740.992.370.990.976 to 0.9990.30830.058
F10 kg103.55 ± 17.01103.39 ± 14.000.982.760.970.809 to 0.9972.03930.009
F15 kg183.20 ± 19.31182.55 ± 18.170.982.160.970.750 to 0.9962.80660.030
F20 kg257.43 ± 18.51253.36 ± 17.950.951.530.980.787 to 0.9962.10820.223
F25 kg304.35 ± 14.59300.41 ± 17.150.961.020.930.431 to 0.9952.50360.248
Power Bk (W)
P5 kg48.27 ± 29.3547.01 ± 30.060.993.470.990.978 to 0.9990.95580.042
P10 kg221.99 ± 58.20222.95 ± 64.920.984.320.970.862 to 0.9975.5571−0.015
P15 kg361.15 ± 70.34358.17 ± 84.300.974.300.960.784 to 0.9958.94330.038
P20 kg442.70 ± 62.11421.90 ± 55.860.943.300.890.239 to 0.9858.25650.353
P25 kg313.80 ± 55.93306.06 ± 55.890.972.470.970.733 to 0.9986.50530.138
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Annino, G.; Romagnoli, C.; Zanela, A.; Melchiorri, G.; Viero, V.; Padua, E.; Bonaiuto, V. Kinematic Analysis of Water Polo Player in the Vertical Thrust Performance to Determine the Force-Velocity and Power-Velocity Relationships in Water: A Preliminary Study. Int. J. Environ. Res. Public Health 2021, 18, 2587. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18052587

AMA Style

Annino G, Romagnoli C, Zanela A, Melchiorri G, Viero V, Padua E, Bonaiuto V. Kinematic Analysis of Water Polo Player in the Vertical Thrust Performance to Determine the Force-Velocity and Power-Velocity Relationships in Water: A Preliminary Study. International Journal of Environmental Research and Public Health. 2021; 18(5):2587. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18052587

Chicago/Turabian Style

Annino, Giuseppe, Cristian Romagnoli, Andrea Zanela, Giovanni Melchiorri, Valerio Viero, Elvira Padua, and Vincenzo Bonaiuto. 2021. "Kinematic Analysis of Water Polo Player in the Vertical Thrust Performance to Determine the Force-Velocity and Power-Velocity Relationships in Water: A Preliminary Study" International Journal of Environmental Research and Public Health 18, no. 5: 2587. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18052587

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop