Participants
Twenty-one professional soccer players aged 28.3±3.9 yrs from the Iranian Persian Gulf Pro League volunteered to participate in this study for a full season. Goalkeepers were excluded from study participation. The athletic/fitness coaches of the club, after obtaining permission from the relevant authorities and the head coach of the club, designed and programmed the soccer training. Before commencing the study, it also received the approval of the research ethics committee from the University of Isfahan (IR.UI.REC.1399.064). All players were informed of the purpose of the study before completing the informed consent. All stages of this study were carried out based on the ethical principles in the Helsinki Declaration.
Study design
The study was performed longitudinally in a full season for 48 weeks in the Persian Gulf Pro League and knockout tournament. GPS (GPSPORTS systems Pty Ltd, Model: SPI High-Performance Unit (HPU), Canberra, Australia) was used for monitoring the external load at each training and match sessions during the whole season. All non-contact injuries were recorded during the season. During the whole season, 7-weeks congested (i.e., two or more matches within 7-days), 30-weeks non-congested; 44 matches, 200 training sessions, and 14126.8 minutes of time played and sessions were held. All workload parameters including AW, CW, ACWR, and Δ-Acute were calculated. Afterwards, each variable was divided into two levels, high and low, and subsequently, the relationship between the variables was measured.
Procedures
External load
GPS receiver specifications. GPSPORTS systems Pty Ltd, recorded all players' activities in training sessions and matches. The GPS-based tracking systems for professional athletes, model SPI HPU features included: 15 Hz position GPS, distance, and speed measurement; accelerometer: 100 Hz, 16 G Tri-Axial-Track impacts, accelerations, and decelerations as well as data source BL; Mag: 50 Hz, Tri-Axial; dimensions: 74 mm × 42 mm × 16 mm; SPI HPU based on Mining/Industrial Strength Electronics design; water resistance and data transmission: infra-red and weighs 56 g. Previous studies have shown that the GPS unit was tested for having a very high accuracy, demonstrated validity and inter-unit reliability, also, the intraclass correlation coefficients were high (>0.95) [20]. Data collecting during training sessions and matches was performed in favorable weather and GPS satellite status.
Data collection. Data collection was completed as in previous studies [21, 22]. At pre-session, we placed upright tracking units in the pouch of the manufacturer supplied belt, then the green light (GPS tracking) flashes were checked. At post-session, after verifying each unit was working properly, tracking units that collected the players' data were placed on the docking station. After 10 minutes, units turned off automatically, and data that were downloaded into docking memory were deleted from the units to prepare for the next session. The GPS units were tuned to the default SPI IQ Absolutes in this study. Duration in minutes and BL was calculated by accelerometer data, and it is designed to reflect both volume and intensity events of the accelerometer (acceleration). BL had replaced the original GPSports BL variable; it is an integrated loading variable used as a training load marker (BL) and work rate marker (BL/min), which we applied as a criterion for the training load in the current study. As to BL calculation, the following steps were repeated for each acceleration level: initialize the BL count to 0; magnitude of the acceleration vector (V) was calculated for the current acceleration (V = ax2 + ay2 + az2); normalize the magnitude vector (NV) by subtracting a national 1G (NV = V - 1.0 G); afterwards unscaled BL (USBL) was calculated through the formula USBLC = NV + [NV3]; then, the scaled BL (SBLC) was calculated taking into account the accelerometer logging rate (100 HZ) and exercise factor (EF) (SBLC = USBLC/100/EF); ultimately, final BL was calculated (BL = BL + SBLC).
Calculation of workload parameters
Workload calculated. In this study, BL was used as the criterion for the training workload [21, 22]. Weekly training workload was used to calculate other workload parameters.
Acute workload (AW) calculation. We recorded the mean weekly AW of the team over the 48 weeks of the season.
Chronic workload (CW) calculation. We recorded the average weekly CW of the team between weeks 4 and 48 (45 weeks). The CW of each player was calculated with the following formula [23, 24]: where n = week number
ACWR calculation. The mean ratio between the AW and CW of the team was calculated and recorded with this variable. The ratio between the individual AW and CW was calculated with the following formula [23, 25]: where n = week number
Δ-Acute load calculation. The AW variation between weeks was calculated through the formula: where n = week number
AW level division. Difference between "high load" and "low load" weeks according to the average weekly AW of the team. "high load" was defined as AW ≥ 571 and "low load" was defined as AW < 571. The cut-off point was established as follows: all the weeks of the season were ordered from highest to lowest AW load, the upper 1.3 was taken as high load and the lower 2.3 as low load.
CW level division. The difference between "high load" and "low load" weeks according to the average weekly CW of the team was calculated. "high load" was defined as CW ≥ 541 and "low load" was defined as CW < 541. The cut-off point was established as follows: all the weeks of the season, with values of CW (45 last weeks) were ordered from highest to lowest CW load, the upper 1.3 was taken as high load and the lower 2.3 as low load.
ACWR level division. The difference between "high load" and "low load" weeks according to the average weekly ACWR of the team was calculated. “high load” was defined as ACWR ≥ 1.18 and “low load” was defined as ACWR < 1.18. The cut-off point was established as follows: all the weeks of the season, with values of ACWR were ordered from highest to lowest ACWR load, the upper 1.3 was taken as high load and the lower 2.3 as low load.
Δ-Acute level division. The difference between "high load" and "low load" weeks according to the average weekly Δ-Acute of the team. “high load” was defined as Δ-Acute ≥ 1.19 and “low load” was defined as Δ-Acute < 1.19. The cut-off point was established as follows: all the weeks of the season, with values of Δ-Acute, were ordered from highest to lowest Δ-Acute load, the upper 1.3 was taken as high load and the lower 2.3 as low load.
Recording and calculating injury
Information on injuries was updated daily by the team's specialized medical staff. Based on a previous study, all injuries were recorded by type, location of the injury, and timing of the injury [26]. The information used for the injuries is as follows:
The number of registered injuries. The total number of non-contact injuries per week for the team, over the 48 weeks of the season, was recorded with this variable.
Weekly injury. The existence or not of a non-contact injury in each of the 48 weeks of the season was recorded.
Statistical analyses
The statistical software IBM Statistics 25 and R Studio 3.6.2. were used for statistical analyses. Accordingly, data were presented as means and standard deviations. A descriptive statistical analysis was realized, indicating the median values and interquartile range of the levels "high load" and "low load" for the variables "AW", "CW", "ACWR" and "Δ-AW," as well as the total values. Non-parametric Mann-Whitney U tests were realized to compare the median of the load levels of the previous variables, checking the existence of statistically significant differences between them. A normality test, Kolmogorov-Smirnov, was performed, determining that the variables "number of injuries" and level of “AW", “CW", “ACWR" and "Δ-AW" did not follow a normal distribution. Additionally, a descriptive analysis of the number of injuries produced in the weeks of high and low load of each one of the variables was completed, as well as the calculation of the median of each one of them, both for the two levels of load as well as for the total. In the purpose of detecting statistically significant inter-group differences between the median of injuries of the "high load" and "low load" levels of the variables "high load" and "low load" for the variables "AW", "CW", "ACWR" and "Δ-AW", non-parametric tests were carried out, taking into account, as factors, the load levels of each variable. A contrast of proportions was made to check the existence of significant differences between the levels of "high load" and "low load" of each variable and the weeks with the injury. To estimate the injury risk associated with high or low load level, the OR and RR were calculated. Finally, the variable "number of injuries" followed a Poisson distribution, so the Poisson test was performed to obtain lambda values (average number of injuries per week for each level of load), and the expected time until a new injury occurs. For checking possible significant differences between load levels, in addition to calculating the rate ratio, their confidence intervals 95% (CI 95%) were stated.