Study design
This study was based on the correlational model, which is a quantitative research method, to evaluate the relationship between DF-ROM and the SSJL task. The study was conducted using a single-blind cross-sectional study design. The Quality Output Checklist and Content Assessment (QuOCCA) checklist was used to enhance the methodological quality of the study [13], and it is presented in Appendix 2. Additionally, The study protocol was pre-registered on Open Science Framework (https://doi.org/10.17605/OSF.IO/KP7XA accessed date: 11.12.2023), and details of the study files are presented on the website (https://osf.io/mtney/ accessed date: 11.12.2023).
Participants
102 female soccer players (age: 21.30 ± 2.14 years, height: 165.10 ± 6.05 cm, weight: 56.56 ± 3.97 kg, BMI: 20.91 ± 2.25 kg∙m2) from Iran's Shahin and Khatoon soccer clubs were included in the study. The participants consisted of 55 amateur and 47 elite-level female soccer players. The participants were selected based on the following criteria in the study: (i) being over 18 years of age, (ii) having at least three years of experience in the field for amateur players, (iii) having competed in the premier, first, or second leagues for professional players (iv) participating in at least two training sessions per week, (v) not having suffered any severe lower limb injuries in the past year (such as ankle or knee sprains, or ACL tears), (vi) being out of the menstrual cycle.
The study sample was determined using the simple random method, and the sample size was determined based on the results of a previous study [9]. A priori power analysis was performed using G*Power software (version 3.1, University of Dusseldorf, Germany) based on the following reference values: α = 0.05, β = 0.80, r = 0.35, two-tailed, and bivariate normal model test. The results showed that the minimum sample size should be 61 participants, while the study included 102 female soccer players. This study was approved by the Vasile Alecsandri University Ethics Committee (decision no: 42/2, 14.12.2023) and was conducted by the Declaration of Helsinki. All subjects provided signed informed consent before being included.
Experimental Procedure
While this study was conducted in Kerman, Iran, measurements were performed on the natural grass fields at the Kargar and Fajr stadiums. A track and field coach with seven years of experience in sports carried out the study protocol. The study protocol was explained to the participants before the measurements to prevent any potential bias caused by learning. Participants were asked to refrain from eating for two hours before taking measurements. Additionally, they were allowed to drink water 30 minutes before measurements. Since the study protocol included jumping and landing tasks that required a high biomechanical load, a standardized 15-minute warm-up protocol was implemented to prevent injuries to the athletes. This warm-up protocol consisted of a 7-minute jog, double-leg squats (2 sets x 10 reps), and dynamic stretches for the hamstrings, quadriceps, and calves (2 sets x 10 reps). Before commencing the measurements, the device calibrations were verified, and all participants underwent the tests while wearing identical attire. Then non-dominant limb test, weight-bearing lunge test, vertical jump test, and soccer-specific jump landing test were performed, respectively. The best test results of the soccer players were used for statistical analysis.
Measurements
Non-dominant limb test.
The non-dominant leg of female soccer players was assessed to measure the DF-ROM level. The ball strike test was used to determine the dominant foot. This test was carried out according to the instructions specified by the researchers, and female soccer players kicked the ball into a target located four meters away. The leg used when hitting the target was considered the dominant leg. Researchers reported 100% agreement between the dominant leg reported by participants and the leg used during the ball-striking test (van Melick et al., 2017).
Weight − bearing lunge Test.
The weight − bearing lunge test was used in the study to determine the DF-ROM level of female soccer players. Researchers stated that the test has high validity (r = 0.71) and reliability (ICC = 0.80 to 0.99) for evaluating DF-ROM [14, 15]. This test was performed using the Inclinometer & Bubble Level (InclineSense Technologies, Iran) mobile app, which has high validity and reliability [16]. Measurements were performed according to instructions in previous study protocol [17]. Before measurement, each device was placed with its long axis on the ground and calibrated to 0° next to the participant. The closed kinetic chain test, which simulates ankle and knee joint loading during the landing mission and provides complete DF-ROM measurement, was chosen to perform the measurement. Participants were asked to place their hands shoulder-width apart on the wall. The non-dominant leg was placed 10 cm from the wall with the knee in line with the second toe. The dominant was placed behind the non-dominant leg with the knee extended. A mobile phone was placed on the flat surface behind the Achilles tendon, approximately one centimeter above the posterior calcaneal tuberosity and perpendicular to the tibia. An inclinometer application called Inclinometer & Bubble Level (or BublePro App) was installed on the Samsung Galaxy A72 to measure ankle DF-ROM. The subjects continued the test by flexing the knee until the heel lifted off the ground. DF-ROM measurement was completed before the heel was lifted off the ground. The test was repeated three times for the non-dominant leg, and the maximum value was used for data analysis. The details of the test protocol are presented in Fig. 1.
Vertical jump test.
The vertical jump height of female soccer players was determined with a traditional vertical jump test to perform the soccer-specific jump test. The vertical jump was measured by having each participant stand beside a wall with their non-dominant side. With feet firmly planted on the ground and hips positioned as close as possible to the wall, the participant extended the arm nearest the wall and reached as high as possible. This point of maximum reach was then marked on the wall. Throughout the exercise, it was essential to ensure that no part of the feet lost contact with the ground (Burkett et al., 2005). After marking the wall, the players jumped as high as possible. The wall was marked again at the highest point reached by the participants. The difference between these two measurements was recorded as the participant's maximum vertical jump. Participants were allowed to perform a maximum of three jumps, and the height of the ball was determined by taking the average of the two highest measurements.
Soccer-specific jump landing test (SSJL)
The SSJL test, an adapted version of the Landing Error Scoring System (LESS), was utilized to assess the landing biomechanics of soccer players. The reliability of the test was confirmed in previous studies [18]. The measurement protocol suggested by the literature was considered to increase the validity and reliability of the test [19, 20]. The test was recorded with an external video, and all female football players continued regular football training that activated the neuromuscular system. Since the age factor affects the LESS score, female soccer players of similar ages were included in this study [19, 20].
To perform the test, participants had to jump over a 7.5 − centimeter cone with the landing point positioned at a distance equal to half their height. Participants were instructed to jump vertically after landing on both feet and head a soccer suspended at a point equivalent to half of their maximum vertical jump height, landing approximately at the same location [9]. Additionally, the test was performed without shoes since footwear can impact the jumping performance of the participants. To record the tests, two digital video cameras (GoPro HERO 9, USA) were positioned 3 meters in front and to the right of the participants [9]. The test was performed three times, and the best result was used for data analysis. Each participant's landing technique was evaluated using Kinovea (version 0.7.10, USA), a valid and reliable software [21]. The evaluation was based on the criteria outlined in Appendix 1, adapted from the Landing Error Scoring System (LESS) introduced by researchers [22]. The LESS test was administered to evaluate each participant's landing technique. It was assessed 17 items on a scale ranging from 0 to 19. A higher score indicating more landing errors suggested potentially high − risk movement patterns. All participants were evaluated by an assessor trained in LESS. Movements of the lower extremities and trunk were analysed at two specific moments: (i) between the first contact of the foot with the ground and (ii) between the first contact with the land and maximum knee flexion. Because several studies suggest that the non-dominant extremity is more susceptible to injury than the dominant extremity during landing strategies in female athletes [23, 24], the evaluation was conducted solely on the non − dominant extremity to streamline the process. The mean score of three jump landing attempts was calculated to obtain an overall assessment of the landing technique.
Statistical Analysis
The study evaluated the relationship between SSJL and DF-ROM of elite and amateur female soccer players. Data were analysed by a researcher blinded to participant characteristics and not involved in the data collection process. The Kolmogorov-Smirnov test was used to check the assumption of normality. The relationship between variables was assessed using Pearson's product-moment correlation analysis. The correlation coefficient is interpreted based on the following references [25]: insignificant (< 0.10), small (0.10 − 0.29), moderate (0.30 − 0.49), strong (0.50 − 0.69), very strong (0.70 − 0.89), or excellent (> 0.90).
In addition to correlation analysis, regression analyses were also performed in the study. While potential outliers in the study data were analysed using Cook's distance, heteroskedasticity was assessed using the Breusch-Pagan test [26]. While outliers were identified in the study data (Fig. 1), heteroscedasticity was also observed (p < 0.05). Therefore, the Akaike Information Criterion (AIC) test was performed to determine the most appropriate method for the regression analysis in the study. A lower AIC value indicates a more appropriate regression analysis method [27]. Linear regression analysis was compared with quantile regression (QR) analysis to determine whether DF-ROM is a predictor of SSJL. Details regarding the comparison results are presented in Table 1.
Table 1
Comparison results of linear and quantile regression analyses based on Akaike's Information Criterion.
𝜏 | | Amateur Female Soccer Players | | Elite Female Soccer Players |
| LR | | QR | | LR | | QR |
| df | AIC | | df | AIC | | df | AIC | | df | AIC |
0.10 | | 3 | 200.95* | | 2 | 215.62 | | 3 | 162.65* | | 2 | 165.78 |
0.25 | | 3 | 200.95* | | 2 | 219.73 | | 3 | 162.65 | | 2 | 160.03* |
0.50 | | 3 | 200.95* | | 2 | 202.75 | | 3 | 162.65 | | 2 | 161.57* |
0.75 | | 3 | 200.95 | | 2 | 194.17* | | 3 | 162.65* | | 2 | 171.44 |
0.90 | | 3 | 200.95 | | 2 | 196.71* | | 3 | 162.65* | | 2 | 196.15 |
Legend. df: Freedom of degrees; AIC: Akaike’s Information Criterion; *: It is represents a more appropriate method to analyse data; 𝜏: Percentile of data (i.e., 0.90 represents the 10% that contains the highest landing scores in the dataset.); LR: Lineer regression; QR: Quantile regression. |
The study data included outliers and revealed heteroscedasticity. Additionally, the AIC results showed conflicting outcomes for different percentiles. Therefore, both linear regression and QR analyses were performed in the study. While evaluating the relationship between two variables in a data set, QR analysis is based on the median values instead of the mean values [28]. In the study, QR analysis was calculated using the following equation:
$$\text{Y} = {\beta }₀\left(\text{q}\right) + {\beta }₁\left(\text{q}\right) \text{*} \text{X}₁ + {\beta }₂\left(\text{q}\right) \text{*} \text{X}₂ + ... + {\beta }\text{ᵣ}\left(\text{q}\right) \text{*} \text{X}\text{ᵣ} + {\epsilon }$$
.
Y represents the dependent variable in this equation, while β defines the estimation coefficients. The expression x represents the independent variable, and q represents the percentiles. Finally, the symbol ε represents the error term in the equation. QR analysis was performed based on previous studies, and five percentiles (10%, 25%, 50%, 75%, 90%) were selected [29]. While the Q10 and Q25 percentiles included female soccer players with lower SSJL landing quality scores, the Q75 and Q90 percentiles included female soccer players with higher SSJL landing quality scores.
R (ver. 4.2.2., Core Team, Vienna, Austuria) software was used for the statistical analysis and data visualization of this study. These operations were performed using the {quantreg}, {olsrr}, {ggplot2}, {gridExtra}, {patchwork}, {gtsummary}, and {performance} packages. Statistical significance was set at p < 0.05 for all analyses. The R codes of the study are presented in Appendix 3.