# Associations Between Variations in Accumulated Workload and Physiological Variables in Young Male Soccer Players Over the Course of a Season

Hadi Nobari, Ana Ruivo Alves, Filipe Manuel Clemente, Jorge Perez-Gomez, Cain C. T. Clark, Urs Granacher, Hassane Zouhal

► To cite this version:

Hadi Nobari, Ana Ruivo Alves, Filipe Manuel Clemente, Jorge Perez-Gomez, Cain C. T. Clark, et al. Associations Between Variations in Accumulated Workload and Physiological Variables in Young Male Soccer Players Over the Course of a Season. Frontiers in Physiology, 2021, 12, 10.3389/fphys.2021.638180. [https://doi.org/10.3389/fphys.2021.638180](https://doi.org/10.3389/fphys.2021.638180). [hal-03217864](https://HAL.inistre.fr/hal-03217864v1)

HAL Id: hal-03217864

https://hal.science/hal-03217864v1

Submitted on 5 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Associations Between Variations in Accumulated Workload and Physiological Variables in Young Male Soccer Players Over the Course of a Season

Hadi Nobari<fnref n="1" />,²,³*, Ana Ruivo Alves<fnref n="4" />,⁵, Filipe Manuel Clemente⁶, Jorge Pérez-Gómez², Cain C. T. Clark⁷, Urs Granacher⁸*† and Hassane Zouhal<fnref n="9" />*†

## OPEN ACCESS

**Edited by:**
Emiliano Cé,
University of Milan, Italy

**Reviewed by:**
Luiz Guilherme Antonacci
Guglielmo,
Federal University of Santa Catarina,
Brazil
João Pedro Duarte,
University of Coimbra, Portugal

**\*Correspondence:**
Hadi Nobari
hadi.nobari1@gmail.com
Urs Granacher
urs.granacher@uni-potsdam.de
orcid.org/0000-0002-7095-813X
Hassane Zouhal
hassane.zouhal@univ-rennes2.fr

† These authors share last authorship

**Specialty section:**
This article was submitted to
*Exercise Physiology,*
a section of the journal
*Frontiers in Physiology*

**Received:** 05 December 2020

**Accepted:** 09 February 2021

**Published:** 18 March 2021

**Citation:**
Nobari H, Alves AR, Clemente FM, Pérez-Gómez J, Clark CCT, Granacher U and Zouhal H (2021) Associations Between Variations in Accumulated Workload and Physiological Variables in Young Male Soccer Players Over the Course of a Season. *Front. Physiol.* 12:638180.
doi: 10.3389/fphys.2021.638180

This study sought to analyze the relationship between in-season training workload with changes in aerobic power (VO₂max), maximum and resting heart rate (HRmax and HRrest), linear sprint medium (LSM), and short test (LSS), in soccer players younger than 16 years (under-16 soccer players). We additionally aimed to explain changes in fitness levels during the in-season through regression models, considering accumulated load, baseline levels, and peak height velocity (PHV) as predictors. Twenty-three male sub-elite soccer players aged 15.5 ± 0.2 years (PHV: 13.6 ± 0.4 years; body height: 172.7 ± 4.2 cm; body mass: 61.3 ± 5.6 kg; body fat: 13.7% ± 3.9%; VO₂max: 48.4 ± 2.6 mL·kg⁻<fnref n="1" />·min⁻<fnref n="1" />), were tested three times across the season (i.e., early-season (EaS), mid-season (MiS), and end-season (EnS) for VO₂max, HRmax, LSM, and LSS. Aerobic and speed variables gradually improved over the season and had a strong association with PHV. Moreover, the HRmax demonstrated improvements from EaS to EnS; however, this was more evident in the intermediate period (from EaS to MiS) and had a strong association with VO₂max. Regression analysis showed significant predictions for VO₂max [F(2, 20) = 8.18, *p* ≤ 0.001] with an *R*² of 0.45. In conclusion, the meaningful variation of youth players' fitness levels can be observed across the season, and such changes can be partially explained by the load imposed.

**Keywords:** internal load, heart rate, linear sprint, aerobic power, football

## INTRODUCTION

Soccer, an intermittent-activity profile sport, integrates explosive activities that require high power output of the lower limb muscles (e.g., sprints, changes of direction, jumps), being interspersed by low-intensity activities with brief recovery intervals (Mohr et al., 2005; Bradley et al., 2014). Young soccer players cover distances of 10–13 km throughout matches and achieve around

<footnote n="1">¹ Department of Physical Education and Sports, University of Granada, Granada, Spain, ² HEME Research Group, Faculty of Sport Sciences, University of Extremadura, Cáceres, Spain, ³ Sports Scientist, Sepahan Football Club, Isfahan, Iran,</footnote>

<footnote n="4">⁴ Department of Arts, Humanities and Sport, Polytechnic Institute of Beja, Beja, Portugal, ⁵ Research Center in Sport Sciences, Health Sciences and Human Development, Vila Real, Portugal, ⁶ Escola Superior Desporto e Lazer, Instituto Politécnico de Viana do Castelo, Viana do Castelo, Portugal, ⁷ Centre for Intelligent Healthcare, Coventry University, Coventry, United Kingdom, ⁸ Division of Training and Movement Sciences, University of Potsdam, Potsdam, Germany,</footnote>

<footnote n="9">⁹ Movement, Sport, Health and Sciences Laboratory (M2S), University of Rennes 2, Rennes, France</footnote>
1,300 individual activities (every 4–6 s), such as accelerations/decelerations (>2 m/s²) and changes of direction, with more than 50 turns per match (Mohr et al., 2005; Dalen et al., 2016). Based on antecedence in the literature, the capacity to perform fast and powerful movements in team sports, specifically in soccer, represents one of the most important skills to develop during training in order to improve performance (Zamparo et al., 2015; Christopher et al., 2016). Indeed, the implementation of neuromuscular training (i.e., strength, plyometric, sprint training) in youth is effective in promoting the athletic development and lowering the risk of injuries, due to the increased plasticity of the neuromuscular system, before, during, and after the period of peak height velocity (PHV) (Rumpf et al., 2013).

In addition to physical performance, soccer technical performance requires the progressive development of neuromuscular abilities in match-related power activities, such as sprinting and jumping (Barnes et al., 2014). However, over periods of high exposure to endurance training (e.g., preseason), improving these abilities presents a challenge, largely because training-induced power and speed adaptations might be impaired by concurrent training (Loturco et al., 2015; Noon et al., 2015). In fact, evidence regarding the interference effect in professional soccer players and elite youth players is well-reported (Loturco et al., 2015; Noon et al., 2015). However, the findings in young players are not well-described. Accordingly, it is necessary to better understand potential interference effects of concurrently performed power and endurance training to optimize players’ physiological adaptations over the soccer season using individualized training load data. Indeed, this knowledge could be helpful to short- and long-term coaching decision-making (Barker and Armstrong, 2010; Álvarez-Kurogi et al., 2019).

As recognized in the literature, there are various key points related to training effectiveness, such as the period of the season (e.g., greater fitness variations are described in preseason compared to in-season), the players' competitive level (e.g., amateurs report superior adaptations than elite players succeeding soccer actions), and the baseline level or the training load imposed (Caldwell and Peters, 2009; Kalapotharakos et al., 2011; Coratella et al., 2016; Clemente et al., 2019). Interestingly, studies based on the analysis of performance during the season have predominantly focused on professional and semiprofessional adult players (Haritonidis et al., 2004; Caldwell and Peters, 2009), yet there is a distinct dearth of studies involving young players (Buchheit et al., 2010; Dragijsky et al., 2017). In fact, managing physical fitness parameters in younger ages is fundamental to the training process, because it is an imperative period of physical development (Lloyd and Oliver, 2012). Thus, it is pertinent that researchers and coaches understand the factors that might affect performance to improve the training process. This will provide useful information on young players' strengths and weaknesses over the season.

In sports practice, the physical fitness level is specific for each individual, and the accurate assessment is imperative to achieve training progress by enabling physiological adaptations (McArdle et al., 2011). Accordingly, heart rate frequency, maximum oxygen

uptake (VO₂max), and the individual anaerobic threshold are useful indicators of cardiorespiratory fitness (McArdle et al., 2011). Previous studies have examined the effects of periodized training on physical fitness (i.e., aerobic power, speed, repeated sprint ability) in soccer players (Owen et al., 2012; Fitzpatrick et al., 2018). However, there is a lack of knowledge regarding the physiological responses to a given training load in young soccer players (Fitzpatrick et al., 2018).

In terms of training load imposed, it is proposed that the load should be sufficient to stimulate and increase players' physical fitness (Impellizzeri et al., 2005). Indeed, given that variegated training loads may yield different adaptations in players, it is important to monitor load, in order to identify possible relationships between the dose and the effect. Although changes in physical fitness variables, preceded by specific training periods, are evident in soccer players, the association between such changes and the accumulated load needs to be further investigated (Clemente et al., 2019). There is only one study available (Clemente et al., 2020a) that has examined the relationship between cardiorespiratory fitness and the accumulated training load parameters over a 4-months in-season period. Results from Clemente et al. (2020a) indicated that training effort is strongly associated with changes in aerobic power. Indeed, no study, to the authors' knowledge, has examined the relationship between training load, maturity status, and different physiological variables (e.g., VO₂max, maximum and resting heart rate [HRmax and HRrest], and speed variables) in elite young male soccer players. Moreover, the relationship between training load during the competitive season and performance is very important to facilitate improvements in strength and conditioning practices. For this reason, the aim of this study was to analyze the relationships between maturity status and training load with variations in physiological variables in soccer players younger than 16 years (under-16 [U16]). Based on the relevant literature (Clemente et al., 2020b; Nobari et al., 2020b), we hypothesized that the accumulated training load might partially explain variation of youth players' fitness levels across the soccer season.

# MATERIALS AND METHODS

## Participants

The study sample consisted of 23 young male sub-elite soccer players who compete at a national U16 level. Players' mean chronological ages were 15.5 ± 0.2 years, VO₂max, 48.5 ± 2.6 mL·kg⁻¹·min⁻¹, and all players were post-PHV (+ 1.9 ± 0.3 years) (Table 1). Among the 23 participants, there were nine defenders, six midfielders, four wingers, and four forward. Inclusion criteria were players (i) who participated in at least 90% of training sessions throughout the season, (ii) did not use dietary supplements during the study, (iii) who were uninjured over the course of the study, and (iv) who did not participate in another training program along with this study. Each player who did not participate in a match during the week participated in a separate training session as a replacement. Exclusion criteria comprised (i) goalkeepers,
<table><caption>TABLE 1 | Descriptive characteristics of 23 soccer players under 16.</caption><thead><tr><th>Variables</th><th>Mean ± SD</th><th>Confidence interval 95%</th></tr></thead><tbody><tr><td>Body height (cm)</td><td>172.7 ± 4.2</td><td>(171–174.4)</td></tr><tr><td>Body mass (kg)</td><td>61.3 ± 5.6</td><td>(59–63.6)</td></tr><tr><td>Sitting height (cm)</td><td>96.6 ± 2.1</td><td>(91.8–93.4)</td></tr><tr><td>Age at PHV (years)</td><td>13.6 ± 0.4</td><td>(13.5–13.7)</td></tr><tr><td>Maturity offset (years)</td><td>1.9 ± 0.3</td><td>(1.7–2.0)</td></tr><tr><td>Chronological age (years)</td><td>15.5 ± 0.2</td><td>(15.3–15.5)</td></tr><tr><td>Training experience (years)</td><td>6.2 ± 1.6</td><td>(5.6–6.9)</td></tr><tr><td>VO<sub>2max</sub> (mL kg<sup>-1</sup>·min<sup>-1</sup>)</td><td>48.4 ± 2.6</td><td>(47.3–49.4)</td></tr><tr><td>Body fat (%)</td><td>13.7 ± 3.9</td><td>(12.1–15.3)</td></tr></tbody></table>

PHV, peak height velocity; VO<sub>2max</sub>, maximal oxygen consumption.

due to the different variations in the physical demands with outfield players. Prior to the start of the study, explanations on the goals and procedures of the study were provided to all participants and their parents/legal representatives. All players were informed about potential risks and benefits of participating in the study. The study was conducted in accordance with the latest version of the Declaration of Helsinki and approved by the Ethics Committee of the University of Isfahan with the ethical code number IR/UI.REC.1397.181. Prior to study commencement, all players and their parents/guardians signed a written informed consent form.

## Experimental Approach to the Problem

The present study was performed during the competitive season, lasting 20 weeks, and included three stages of evaluation (before, mid, and after the end of the season). The season was divided into three periods: (i) early-season (EaS) = week (W) 1 to W7; (ii) mid-season (MiS) = W8–W13; and (iii) end-season (EnS) = W14–W20 according to **Figure 1**. The first assessment was in the week before the start of the season, the second test point was in the 11th week of the mid-season, and the third assessment was in the week after the season ended. During each test, players were measured on consecutive days. Individual indicators of each player were assessed on the first test day (i.e., anthropometric and body composition variables) and maturity status to determine age at PHV. On the next test, the linear sprint short test (LSS), the linear sprint medium test (LSM), and the intermittent fitness test 30–15 (30–15IFT) were performed. Evaluations of each stage were performed at the same time, temperature, and humidity (Medicine ACoS, 2013). Maximal heart rate (HRmax) was also recorded during the 30–15IFT test. Resting heart rate (HRrest) was measured in lying position 1 day after resting when waking up in the morning. A familiarization session was scheduled 1 week prior to the assessment. For this cohort study, 30 min after each training session, all players reported load of training, and then each “training load” was calculated, alongside training time, to discern accumulated workload for any period. Each week, the training days included the following: (i) technical and skills with moderate-intensity aerobic training; (ii) tactical and full body-weight training; (iii) high-intensity training or small side games; (iv) tactical with agility and coordination training; (v) and finally,

an official or friendly match. Athletes were kindly asked not to change their diet throughout the study period.

## Anthropometrics and Body Composition

Stature was measured with a stadiometer [Seca model 213, Germany (accuracy of ± 5 mm)]. Participants stood barefoot and touched their lower back, head, shoulder blades, buttocks, and heel to the stadiometer, with their feet placed together. After adopting the correct position, height was measured vertically by the rod from above the head. To measure the sitting height, participants were required to sit on a 50-cm-height table (forward face) with their lower back as close as possible to the stadiometer. They kept their upper body straight and put their hands on their feet. Then, the height measured between the highest point of the head (vertex) and the (ischial spines) plane where they were sitting. Body mass was measured using weighing scales [Seca model 813, United Kingdom (accuracy ± 0.1 per kg)]. Each subject, with light clothes, stood, and the scale recorded the subjects' mass. The technical error of anthropometric measurements, interobserver and intraobserver, with the same observer was less than 3% (Nobari et al., 2020a,b). A practical method for the determination of maturity offset was applied using three somatic dimensions (height, sitting height, leg length) and chronological age as well as their interactions: $[\text{maturity offset} = -9.236 + 0.0002708 (\text{leg length} \times \text{sitting height}) - 0.001663 (\text{chronological age} \times \text{leg length}) + 0.007216 (\text{chronological age} \times \text{sitting height}) + 0.02292 (\text{mass by height ratio})]$. This equation has previously been validated by Mirwald et al. (2002).

The subcutaneous fat thickness of the triceps and subscapular points of the body were measured. Percentage body fat was estimated using the equation of Slaughter et al. (1988). Skin thickness was obtained by using calibrated Lafayette skinfold calipers (United States, with an accuracy of ± 0.1 mm). All measurements were performed twice on the right side of the body; the final score was recorded as the mean of two measurements. If technical error was higher than 5%, the measurements were executed again, and the median value was recorded (Nobari et al., 2020a,b). All anthropometric and body composition measurements were taken by a skilled practitioner, with a 5-year history in this area, and all measurements were assessed in the morning (Arazi et al., 2015; Rahmat et al., 2016).

## Monitoring Internal Training Loads

Players were monitored daily for their rating perceived of exertion (RPE) using the Category-Ratio-10 Borg's scale, which is a valid and reliable scale to estimate the intensity of a session (Foster et al., 2001). To the question, “How intense was your session?” players answered in the interval range of 1 (minimum effort) and 10 (maximal effort). Players answered 30 min after the end of each training session. Additionally, the duration of the training sessions, in minutes, was recorded. As a measure of internal load, the session RPE was calculated by multiplying the score in Category-Ratio-10 scale (11 numbers) by the duration of the session in minutes. Players were previously familiarized with the scale, having used it in the previous 2 years in the club. In this study, the accumulated load (for training and
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Research outline of the monitoring on training, match, and assessed sessions during the competition season. EaS, early season = before league for assessments and W1–W7 for the cohort study; MiS, midseason = mid league for assessments and W8–W13 for the cohort study; EnS, end-season = after league for assessments and W14–W20 for the cohort study; W, week; TS, training sessions.</figcaption></figure>

competition) was used for 20 weeks. These weeks of the full-competition season were divided into three periods: early-season (EaS), W1–W7; mid-season (MiS), W8–W13; and end-season (EnS), W14–W20, for analysis used the mean of the accumulated load of periods (MALP).

## Intermittent Fitness Test 30–15 (30–15IFT)

The 30–15IFT was used to calculate the VO₂max of the players. The test started at an average speed of 8 km·h⁻¹ between a 40 m shuttle. Players started running 30 s after hearing the first beep (from line A to B and from B to C and then back). Each stage consisted of running for 30 s and then resting for 15 s. Thereafter, every stage had a 0.5 km·h⁻¹ increment until volitional exhaustion (Buchheit et al., 2010). Each person's score was the last stage of speed that he could not continue, or three times not making it to the 2 m lines. At this time, the final velocity of the IFT (VIFT) was recorded. All these steps were performed after a standardized warm-up of 10 min. The following equation was used to determine the VO₂max (mL·kg⁻¹·min⁻¹) of each player (Buchheit et al., 2010) = 28.3 – (2.15 × 1) – (0.741 × age) – (0.0357 × mass) + (0.0586 × age × VIFT) + (1.03 × VIFT). The test–retest reliability was calculated, using intraclass coefficient (ICC), to be 0.86.

## Heart Rate Measurements

The HR variables were measured using Mi-Band 3, Xiaomi Company, made in China. The HRmax was recorded during the IFT30–15 test, to measure each player's HRrest at each assessment stage. In the morning, when waking up, they performed three assessments over 3 days. Information was received by the researcher every day, and the average of 3 days was considered as a criterion for analysis.

## Linear Sprint Test

To measure linear sprint performance, the LSS 10 m and the LSM 20 m tests were used. These two variables were performed with the Newtest Powertimer 300-series device made in Finland, and times were recorded. Each distance was repeated three times with 3 min of rest. The best time in each distance was considered as a criterion for analysis. The ICCs in these tests was LSS = 0.91

and LSM = 0.94. The photocells were adjusted from the athlete's pelvis height. All players followed a standardized warm-up of 10–15 min, including jogging, running drills (e.g., high knee, butt kicks, carioca), and two to three linear sprints. After warming up, the players stood at a distance of 70 cm from the start line, and then, on the command of the instructor, the player started covering the specified distance at full speed (Mirkov et al., 2008).

## Statistical Analyses

Based on a previous work that highlighted large to very large correlations between physiological variables in soccer players (Clemente et al., 2019), we computed the sample size necessary to achieve a power of, at least, 0.90. Accordingly, using a two-tailed hypothesis, expected large to very large correlations, and an α error probability of 0.05, 21 participants would be required to achieve the aforementioned power. The normality of the data was checked by Shapiro–Wilk test. Data are presented as means and standard deviations. In cases where data were normally distributed, Pearson correlation coefficients were calculated for measures of physical fitness and training load used for normal data (i.e., HRmax and MALP of between periods), and if normality of data was violated, Spearman correlation coefficients were computed. The following threshold was used to determine the amplitude of correlation levels (Batterham and Hopkins, 2006): < 0.1 = trivial, 0.1–0.3 = small, 0.3–0.5 = moderate, 0.5–0.7 = large, 0.7–0.9 = very large, and > 0.9 = nearly perfect. Multiple linear (least squares) regression analysis was calculated to predict the percentage of change in fitness levels (i.e., VO₂max, HR, and speed variables) based on accumulated load, baseline levels, and PHV. In accordance with previous recommendations (Jenkins and Quintana-Ascencio, 2020), the Akaike information criterion (AIC) for each model's regression was additionally calculated, to support inferences about the model's suitability. We analyzed the ICCs for the reliability evaluated for within-session outcome of all measures but used the 30–15IFT test–retest reliability. The ICC > 0.7 was suitable (Baumgartner and Chung, 2001). Statistical significance was accepted, a priori, at *p* < 0.05. All statistical analyses were performed using GraphPad Prism 8.0.1 software, except for multiple linear regression and AIC, which were calculated using R software (version 1.2.5019).
# RESULTS

**Figure 2** shows the correlation coefficient between PHV with mean differences between assessments and the accumulated load of between periods at the 95% confidence interval (EaS to MiS, MiS to EnS, and EaS to EnS). There were significant correlations between VO₂max at EaS to MiS (r = 0.49; 95% CI, 0.09–0.76; p = 0.02), LSS at EaS to MiS (r = −0.46; 95% CI, −0.74 to −0.05; p = 0.02), and LSS at MiS to EnS (r = 0.51; 95% CI, 0.11–0.76; p ≤ 0.001) with PHV. Likewise, VO₂max at EaS to EnS (r = 0.67 large; 95% CI, 0.35–0.85; p ≤ 0.001) and HRmax at EaS to MiS (r = 0.46; 95% CI, 0.052–0.74; p = 0.03) were associated with VO₂max at MiS to EnS, respectively. In addition, HRmax at EaS to MiS (r = 0.48; 95% CI, 0.07–0.75; p = 0.02) was related to VO₂max at EaS to EnS. For HRmax, there were strong associations between EaS to EnS and MiS to EnS (r = 0.94; 95% CI, 0.86–0.95; p ≤ 0.001). Further, HRrest at EaS to EnS was related to HRrest at EaS to MiS (r = 0.41; 95% CI, −0.01 to 0.72; p = 0.05), and MiS to EnS (r = 0.75; 95% CI, 0.48–0.89; p ≤ 0.0001), respectively.

Overall, the relations MALP at EaS to MiS with HRmax at EaS to EnS (r = −0.42; 95% CI, 0.01–0.71]; p = 0.047) and MALP at MiS to EnS (r = −0.69; 95% CI, −0.86 to −0.38]; p ≤ 0.001) were moderate and large (Table 2).

Multiple linear regression analyses were performed to predict the percentage of change in fitness levels (i.e., VO₂max, HR, and speed variables) based on accumulated load, baseline levels, and PHV. The regression for VO₂max was significant [F(8, 14) = 3.72, p = 0.02], with an R² of 0.68. Participants predicted VO₂max (Y) was equal to β0 + β1 (accumulated load) + β2

(PHV) + β3 (Maturity offset) + β4 (EaS VO₂max) + β5 (baseline HRmax) + β6 (EaS HRrest) + β7 (EaS LSS) + β8 (EaS LSM). Where accumulated load was coded or measured as arbitrary units, maturity status coded or measured as years, fitness status coded or measured as mL·kg⁻¹·min⁻¹, bpm, and s. Participants predicted VO₂max levels during the season decreased −0.84 mL·kg⁻¹·min⁻¹ for each mL·kg⁻¹·min⁻¹ of VO₂max, as well as decreasing −0.30 mL·kg⁻¹·min⁻¹ for each bpm of HRₜ, both baseline VO₂max and HRₜ were significant predictors of aerobic power levels (Table 3). An additional predictive equation was made for these two variables (VO₂max and HRₜ), which was significant [F(2, 20) = 8.18, p ≤ 0.001], with an R² of 0.45. Participants' predicted VO₂max (Y) was equal to β0 + β1 (VO₂max) + β2 (HRₜ). Participants' VO₂max levels during the season decreased −0.89 mL·kg⁻¹·min⁻¹ for each mL·kg⁻¹·min⁻¹ of VO₂max, as well as decreases −0.17 mL·kg⁻¹·min⁻¹ for each bpm of HRₜ (Table 3). Residual plots are depicted in Figure 3.

For LSS [F(8, 14) = 0.38, p = 0.92, R² = 0.18] and LSM [F(8, 14) = 0.55, p = 0.80, R² = 0.24], multiple linear regression did not yield any significant findings (Table 4).

For HRmax [F(8, 14) = 1.75, p = 0.17, R² = 0.50] and HRrest [F(8, 14) = 1.30, p = 0.32, R² = 0.43], multiple linear regression did not yield any significant findings (Table 5).

Post hoc power analysis was conducted to detect the achieved power for each multiple linear regression analysis, based on a two-tailed hypothesis, a sample size of 23, nine predictors, and the achieved effect size. For percentage of change in VO₂max, with accumulated load, baseline fitness levels, and PHV, the achieved

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Correlation coefficients (95% CI) the peak height velocity soccer player with the VO₂max (ml·kg⁻¹·min⁻¹), maximal oxygen consumption; HRmax (bpm), maximal heart rate; HRₜ (bpm), resting heart rate; LSS (s), linear sprint short test; LSM (s), linear sprint medium test; MALP (A.U.), the mean of the accumulated load of period; A.U., arbitrary unit.</figcaption></figure>
<table><caption>TABLE 2 | Pearson and Spearman correlation analysis between the mean differences between assessments and the average accumulated load of between periods.</caption><thead><tr><th>Variable</th><th>β0</th><th>β1</th><th>β2</th><th>β3</th><th>β4</th><th>β5</th><th>β6</th><th>β7</th><th>β8</th><th>β9</th><th>β10</th><th>β11</th><th>β12</th><th>β13</th><th>β14</th><th>β15</th><th>β16</th><th>β17</th><th>β18</th><th>β19</th><th>β20</th><th>β21</th></tr></thead><tbody><tr><td>PHV (β0)</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VO<sub>2max</sub> 1 (β1)</td><td><strong>0.49</strong></td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VO<sub>2max</sub> 2 (β2)</td><td>-0.18</td><td>-0.36</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VO<sub>2max</sub> 3 (β3)</td><td>0.15</td><td>0.31</td><td><strong>0.67</strong></td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HR<sub>max</sub> 1 (β4)</td><td>0.02</td><td>0.15</td><td><strong>0.46</strong></td><td><strong>0.48</strong></td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HR<sub>max</sub> 2 (β5)</td><td>0.21</td><td>-0.05</td><td>-0.07</td><td>-0.09</td><td>-0.38</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HR<sub>max</sub> 3 (β6)</td><td>0.20</td><td>0.01</td><td>0.01</td><td>0.01</td><td>-0.04</td><td><strong>0.94</strong></td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HR<sub>rest</sub> 1 (β7)</td><td>0.06</td><td>0.05</td><td>-0.08</td><td>-0.10</td><td>0.17</td><td>0.10</td><td>0.14</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HR<sub>rest</sub> 2 (β8)</td><td>0.09</td><td>-0.12</td><td>0.05</td><td>0.02</td><td>0.16</td><td>-0.06</td><td>0.00</td><td>-0.25</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HR<sub>rest</sub> 3 (β9)</td><td>0.16</td><td>0.01</td><td>-0.03</td><td>0.02</td><td>0.26</td><td>0.07</td><td>0.16</td><td><strong>0.41</strong></td><td><strong>0.75</strong></td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>LSS 1 (β10)</td><td><strong>-0.46</strong></td><td><strong>-0.57</strong></td><td>0.07</td><td>-0.23</td><td>-0.27</td><td>0.19</td><td>0.13</td><td>-0.19</td><td>-0.23</td><td>-0.39</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>LSS 2 (β11)</td><td><strong>0.51</strong></td><td>0.38</td><td>0.01</td><td>0.21</td><td>-0.13</td><td>0.17</td><td>0.15</td><td>-0.25</td><td>0.26</td><td>0.08</td><td>-0.17</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>LSS 3 (β12)</td><td>0.13</td><td>-0.01</td><td>-0.03</td><td>-0.03</td><td>-0.29</td><td>0.22</td><td>0.17</td><td>-0.33</td><td>0.14</td><td>-0.10</td><td><strong>0.54</strong></td><td><strong>0.69</strong></td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>LSM 1 (β13)</td><td>-0.05</td><td>-0.01</td><td>0.33</td><td>0.23</td><td>0.25</td><td>-0.13</td><td>-0.10</td><td>-0.17</td><td>0.00</td><td>-0.14</td><td>-0.15</td><td>0.12</td><td>-0.07</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>LSM 2 (β14)</td><td>0.08</td><td>0.15</td><td>-0.16</td><td>-0.23</td><td>-0.20</td><td>-0.01</td><td>-0.06</td><td>0.29</td><td>-0.20</td><td>0.06</td><td>-0.25</td><td>0.22</td><td>-0.01</td><td>-0.17</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>LSM 3 (β15)</td><td>-0.05</td><td>0.11</td><td>-0.01</td><td>-0.09</td><td>-0.12</td><td>-0.06</td><td>-0.11</td><td>0.15</td><td>-0.25</td><td>-0.08</td><td>-0.20</td><td>0.23</td><td>-0.02</td><td><strong>0.44</strong></td><td><strong>0.72</strong></td><td>1.00</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MALP 1 (β16)</td><td>0.10</td><td>0.03</td><td>0.35</td><td>0.31</td><td>0.11</td><td>0.35</td><td><strong>0.42</strong></td><td>-0.08</td><td>0.14</td><td>0.11</td><td>-0.03</td><td>0.24</td><td>0.08</td><td>0.33</td><td>0.20</td><td>0.36</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MALP 2 (β17)</td><td>-0.01</td><td>0.13</td><td>-0.16</td><td>-0.04</td><td>-0.07</td><td>-0.31</td><td>-0.36</td><td>-0.13</td><td>-0.19</td><td>-0.24</td><td>0.09</td><td>-0.15</td><td>0.01</td><td>-0.21</td><td>-0.24</td><td>-0.32</td><td><strong>-0.69</strong></td><td>1.00</td><td></td><td></td><td></td><td></td></tr><tr><td>MALP 3 (β18)</td><td>0.11</td><td>0.20</td><td>0.22</td><td>0.32</td><td>0.03</td><td>0.01</td><td>0.02</td><td>-0.27</td><td>-0.09</td><td>-0.18</td><td>0.09</td><td>0.09</td><td>0.11</td><td>0.11</td><td>-0.07</td><td>0.01</td><td>0.30</td><td><strong>0.48</strong></td><td>1.00</td><td></td><td></td><td></td></tr><tr><td>VIFT 1 (β19)</td><td><strong>-0.42</strong></td><td>-1.00</td><td>0.10</td><td><strong>-0.59</strong></td><td>-0.19</td><td>0.06</td><td>-0.01</td><td>-0.04</td><td>-0.11</td><td>-0.13</td><td>0.34</td><td>-0.25</td><td>0.19</td><td>-0.02</td><td>0.04</td><td>0.02</td><td>-0.03</td><td>-0.13</td><td>-0.21</td><td>1.00</td><td></td></tr><tr><td>VIFT 2 (β20)</td><td>0.20</td><td>0.10</td><td>-1.00</td><td><strong>-0.75</strong></td><td>-0.22</td><td>0.07</td><td>-0.01</td><td>0.09</td><td>-0.09</td><td>-0.06</td><td>0.04</td><td>-0.02</td><td>0.02</td><td>-0.18</td><td>0.01</td><td>-0.09</td><td>-0.35</td><td>0.16</td><td>-0.22</td><td>-0.09</td><td>1.00</td></tr><tr><td>VIFT 3 (β21)</td><td>-0.12</td><td><strong>-0.59</strong></td><td><strong>-0.74</strong></td><td><strong>-1.00</strong></td><td>-0.30</td><td>0.10</td><td>-0.01</td><td>0.05</td><td>-0.15</td><td>-0.13</td><td>0.26</td><td>-0.18</td><td>0.15</td><td>-0.15</td><td>0.04</td><td>-0.06</td><td>-0.30</td><td>0.04</td><td>-0.32</td><td><strong>0.60</strong></td><td><strong>0.74</strong></td><td>1.00</td></tr></tbody></table>

Significant differences (p ≤ 0.05) are highlighted in bold font.

PHV, peak height velocity; VO<sub>2max</sub>, maximal oxygen consumption; HR<sub>max</sub>, maximal heart rate; HR<sub>rest</sub>, resting heart rate; LSS, linear sprint short test; LSM, linear sprint medium test; MALP, the mean of the accumulated load of period; VIFT, final velocity of the 30–15<sub>IFT</sub>; 1, 2, and 3 = the mean differences between assessments (EaS to MiS, and MiS to EnS and EaS to EnS), respectively.
<table><caption>TABLE 3 | Multiple linear regression analysis: percentage of change in VO₂max with accumulated load, baseline fitness levels, and PHV.</caption><thead><tr><th>Variable</th><th>β</th><th>Estimate</th><th>|t|</th><th>P-value</th><th>95% CI for estimated</th><th>Total predict</th></tr></thead><tbody><tr><td>VO<sub>2max</sub> (%)</td><td>β0</td><td>37.13</td><td>0.67</td><td>0.51</td><td>-81.9 to 156.1</td><td><b>R<sup>2</sup> = 0.68</b></td></tr><tr><td>ACL (A.U.)</td><td>β1</td><td>0.00</td><td>0.91</td><td>0.38</td><td>-0.002 to 0.001</td><td><b>Adjusted</b></td></tr><tr><td>Maturity offset (years)</td><td>β2</td><td>2.91</td><td>1.23</td><td>0.24</td><td>-2.17 to 7.98</td><td><b>R<sup>2</sup> = 0.50</b></td></tr><tr><td>PHV (years)</td><td>β3</td><td>3.96</td><td>1.30</td><td>0.21</td><td>-2.56 to 10.48</td><td><b>P = 0.02</b></td></tr><tr><td>VO<sub>2max</sub> (mL·kg<sup>-1</sup>·min<sup>-1</sup>)</td><td>β4</td><td>-0.84</td><td>3.78</td><td><b>≤0001</b></td><td>-1.32 to -0.36</td><td><b>AIC = 112.1</b></td></tr><tr><td>HR<sub>max</sub> (bpm)</td><td>β5</td><td>-0.08</td><td>0.64</td><td>0.53</td><td>-0.36 to 0.2</td><td></td></tr><tr><td>HR<sub>rest</sub> (bpm)</td><td>β6</td><td>-0.30</td><td>2.08</td><td><b>0.06</b></td><td>-0.61 to 0.01</td><td></td></tr><tr><td>LSS (s)</td><td>β7</td><td>5.60</td><td>1.74</td><td>0.10</td><td>-1.31 to 12.50</td><td></td></tr><tr><td>LSM (s)</td><td>β8</td><td>-0.89</td><td>0.45</td><td>0.66</td><td>-5.15 to 3.36</td><td></td></tr></tbody></table>

Significant differences ($p \le 0.05$) are highlighted in bold font.

PHV, peak height velocity; VO₂max, maximal oxygen consumption; HRmax, maximal heart rate; HRrest, resting heart rate; LSS, linear sprint short test; LSM, linear sprint medium test; ACL, all of accumulated load of 20 weeks; %, percentage of change in between assessments from EaS to EnS; CI, confidence interval; AIK, Akaike information criterion; A.U., arbitrary unit.

<table><caption>TABLE 4 | Multiple linear regression analysis: percentage of change in speed variables with accumulated load, baseline fitness levels, and PHV.</caption><thead><tr><th>Variable</th><th>β</th><th>Estimate</th><th>|t|</th><th>P-value</th><th>95% CI for estimated</th><th>Total predict</th></tr></thead><tbody><tr><td>LSS (%)</td><td>β0</td><td>-94.42</td><td>0.65</td><td>0.529</td><td>-408.2 to 219.3</td><td><b>R<sup>2</sup> = 0.18</b></td></tr><tr><td>ACL (A.U.)</td><td>β1</td><td>0.001</td><td>0.48</td><td>0.639</td><td>-0.003 to 0.004</td><td><b>Adjusted</b></td></tr><tr><td>Maturity offset (yrs)</td><td>β2</td><td>5.20</td><td>0.84</td><td>0.414</td><td>-8.18 to 18.58</td><td><b>R<sup>2</sup> = -0.29</b></td></tr><tr><td>PHV (yrs)</td><td>β3</td><td>5.58</td><td>0.7</td><td>0.498</td><td>-11.60 to 22.76</td><td><b>P = 0.92</b></td></tr><tr><td>VO<sub>2max</sub> (mL·kg<sup>-1</sup>·min<sup>-1</sup>)</td><td>β4</td><td>-0.89</td><td>1.51</td><td>0.153</td><td>-2.15 to 0.37</td><td><b>AIC = 156.71</b></td></tr><tr><td>HR<sub>max</sub> (bpm)</td><td>β5</td><td>0.31</td><td>0.93</td><td>0.371</td><td>-0.41 to 1.03</td><td></td></tr><tr><td>HR<sub>rest</sub> (bpm)</td><td>β6</td><td>-0.22</td><td>0.59</td><td>0.566</td><td>-1.04 to 0.59</td><td></td></tr><tr><td>LSS (s)</td><td>β7</td><td>0.11</td><td>0.01</td><td>0.989</td><td>-18.11 to 18.33</td><td></td></tr><tr><td>LSM (s)</td><td>β8</td><td>-3.15</td><td>0.60</td><td>0.557</td><td>-14.37 to 8.08</td><td></td></tr></tbody></table><table><thead><tr><th>Variable</th><th>β</th><th>Estimate</th><th>|t|</th><th>P-value</th><th>95% CI for estimated</th><th>Total predict</th></tr></thead><tbody><tr><td>LSM (%)</td><td>β0</td><td>-104.2</td><td>1.08</td><td>0.298</td><td>-310.8 to 102.4</td><td><b>R<sup>2</sup> = 0.24</b></td></tr><tr><td>ACL (A.U.)</td><td>β1</td><td>=0001</td><td>0.95</td><td>0.359</td><td>-0.001 to 0.003</td><td><b>Adjusted</b></td></tr><tr><td>Maturity offset (yrs)</td><td>β2</td><td>2.56</td><td>0.62</td><td>0.544</td><td>-6.25 to 11.37</td><td><b>R<sup>2</sup> = -0.20</b></td></tr><tr><td>PHV (yrs)</td><td>β3</td><td>7.65</td><td>1.45</td><td>0.169</td><td>-3.67 to 18.96</td><td><b>P = 0.80</b></td></tr><tr><td>VO<sub>2max</sub> (mL·kg<sup>-1</sup>·min<sup>-1</sup>)</td><td>β4</td><td>0.28</td><td>0.73</td><td>0.478</td><td>-0.5484 to 1.11</td><td><b>AIC = 137.48</b></td></tr><tr><td>HR<sub>max</sub> (bpm)</td><td>β5</td><td>0.09</td><td>0.41</td><td>0.686</td><td>-0.38 to 0.57</td><td></td></tr><tr><td>HR<sub>rest</sub> (bpm)</td><td>β6</td><td>0.05</td><td>0.21</td><td>0.838</td><td>-0.49 to 0.59</td><td></td></tr><tr><td>LSS (s)</td><td>β7</td><td>2.06</td><td>0.37</td><td>0.718</td><td>-9.93 to 14.06</td><td></td></tr><tr><td>LSM (s)</td><td>β8</td><td>-3.61</td><td>1.05</td><td>0.313</td><td>-10.99 to 3.78</td><td></td></tr></tbody></table>

Significant differences ($p = 0.05$) are highlighted in bold.

PHV, Peak height velocity; VO₂max, Maximal oxygen consumption; HRmax, Maximal heart rate; HRrest, resting heart rate; LSS, Linear sprint short test; LSM, linear sprint medium test; ACL, All of accumulated load of 20 weeks; %, percentage of change in between assessments from EaS to EnS; CI, Confidence interval; AIK, Akaike information criterion; A.U., Arbitrary unit.

power was 0.95. For percentage of change in speed variables (LSS and LSM), with accumulated load, baseline fitness levels, and PHV, the achieved power was 0.68, due to the trivial effect sizes. Finally, for percentage of change in HR measurements (HRmax and HRrest), with accumulated load, baseline fitness levels, and PHV, the achieved power was 0.88.

## DISCUSSION

This study sought to discern the relationship between in-season training workload with changes in VO₂max, heart rate,

and speed variables in male U16 soccer players. Accordingly, large and moderate correlations were reported between aerobic power, HR, speed variables, and MALP between the different periods of the soccer season. We additionally aimed to explain changes in fitness levels during the in-season through regression models considering accumulated load, baseline levels, and PHV as predictors. Significant predictors (i.e., VO₂max and HRrest) of aerobic power levels during the competition season were found, but not for the speed variables (i.e., LSS and LSM), whereas for HR, only during the competition season, HRmax was reported as a significant predictor of percentage of change in HRmax levels.
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

**FIGURE 3** | Multiple linear regression analysis was calculated to predict the percentage of change in fitness levels includes: VO₂max, (A1); LSS (B1); LSM (C1); HRmax (D1); and HRrest (E1) based on accumulated load, baseline levels, and PHV with residual plot for each analysis VO₂max, (A2); LSS (B2); LSM (C2); HRmax (D2); and HRrest (E2), the difference between the actual value of the dependent variable and the value predicted by the residual provided. VO₂max, maximal oxygen consumption.
<table><caption>TABLE 5 | Multiple linear regression analysis: percentage of change in HR measurements with accumulated load, baseline fitness levels, and PHV.</caption><thead><tr><th>Variable</th><th>β</th><th>Estimate</th><th>|t|</th><th>P-value</th><th>95% CI for estimated</th><th>Total predict</th></tr></thead><tbody><tr><td>HR<sub>max</sub> (%)</td><td>β0</td><td>110.3</td><td>1.57</td><td>0.138</td><td>-40.24 to 260.8</td><td><b>R<sup>2</sup> = 0.50</b></td></tr><tr><td>ACL (A.U.)</td><td>β1</td><td>≤-0.001</td><td>0.85</td><td>0.409</td><td>-0.002 to ≤ 0.001</td><td><b>Adjusted</b></td></tr><tr><td>Maturity offset (years)</td><td>β2</td><td>-1.46</td><td>0.49</td><td>0.634</td><td>-7.88 to 4.96</td><td><b>R<sup>2</sup> = 0.21</b></td></tr><tr><td>PHV (years)</td><td>β3</td><td>-5.44</td><td>1.42</td><td>0.179</td><td>-13.68 to 2.81</td><td><b>P = 0.17</b></td></tr><tr><td>VO<sub>2max</sub> (mL·kg<sup>-1</sup>·min<sup>-1</sup>)</td><td>β4</td><td>0.27</td><td>0.95</td><td>0.357</td><td>-0.34 to 0.88</td><td><b>AIC = 122.93</b></td></tr><tr><td>HR<sub>max</sub> (bpm)</td><td>β5</td><td>-0.47</td><td>2.93</td><td><b>0.011</b></td><td>-0.82 to -0.13</td><td></td></tr><tr><td>HR<sub>rest</sub> (bpm)</td><td>β6</td><td>0.36</td><td>1.97</td><td>0.069</td><td>-0.03 to 0.75</td><td></td></tr><tr><td>LSS (s)</td><td>β7</td><td>-3.12</td><td>0.77</td><td>0.456</td><td>-11.86 to 5.62</td><td></td></tr><tr><td>LSM (s)</td><td>β8</td><td>0.79</td><td>0.31</td><td>0.759</td><td>-4.6 to 6.17</td><td></td></tr></tbody></table><table><thead><tr><th>Variable</th><th>β</th><th>Estimate</th><th>|t|</th><th>P value</th><th>95% CI for estimated</th><th>Total predict</th></tr></thead><tbody><tr><td>HR<sub>rest</sub> (%)</td><td>β0</td><td>60.38</td><td>0.74</td><td>0.468</td><td>-113.1 to 233.9</td><td><b>R<sup>2</sup> = 0.43</b></td></tr><tr><td>ACL (A.U.)</td><td>β1</td><td>≤ 0.001</td><td>0.93</td><td>0.368</td><td>≤-0.001 to 0.003</td><td><b>Adjusted</b></td></tr><tr><td>Maturity offset (years)</td><td>β2</td><td>-4.42</td><td>1.28</td><td>0.221</td><td>-11.82 to 3</td><td><b>R<sup>2</sup> = 0.01</b></td></tr><tr><td>PHV (years)</td><td>β3</td><td>-6.33</td><td>1.43</td><td>0.175</td><td>-15.8 to 3.18</td><td><b>P = 0.32</b></td></tr><tr><td>VO<sub>2max</sub> (mL·kg<sup>-1</sup>·min<sup>-1</sup>)</td><td>β4</td><td>-0.1</td><td>0.31</td><td>0.765</td><td>-0.8 to 0.6</td><td><b>AIC = 129.46</b></td></tr><tr><td>HR<sub>max</sub> (bpm)</td><td>β5</td><td>-0.02</td><td>0.13</td><td>0.901</td><td>-0.42 to 0.38</td><td></td></tr><tr><td>HR<sub>rest</sub> (bpm)</td><td>β6</td><td>-0.13</td><td>0.62</td><td>0.546</td><td>-0.58 to 0.32</td><td></td></tr><tr><td>LSS (s)</td><td>β7</td><td>-5.1</td><td>1.1</td><td>0.296</td><td>-15.17 to 4.98</td><td></td></tr><tr><td>LSM (s)</td><td>β8</td><td>3.406</td><td>1.18</td><td>0.259</td><td>-2.79 to 9.61</td><td></td></tr></tbody></table>

Significant differences ($p \le 0.05$) are highlighted in bold font. PHV, peak height velocity; $VO_{2max}$, maximal oxygen consumption; $HR_{max}$, maximal heart rate; $HR_{rest}$, resting heart rate; LSS, linear sprint short test; LSM, linear sprint medium test; ACL, accumulated load of 20 weeks; %, percentage of change in between assessments from EaS to EnS; CI, confidence interval; AIK, Akaike information criterion; A.U., arbitrary unit.

The capacity of soccer players to develop and preserve high physical fitness levels over the season is paramount to success (Krustrup et al., 2005; Bloomfield et al., 2007). The current study indicated that aerobic capacity improved from the EaS to the EnS, with more evident increases from the MiS to the EnS. Moreover, $VO_{2max}$ between EaS to MiS was significantly related to PHV. Following the present study, Clemente et al. (2020a) also reported significant improvements in maximal aerobic speed over the season in young soccer players, with the greatest improvements at the EnS (4 months after baseline test). However, Dragijsky et al. (2017) showed an improvement in endurance, with the highest levels reported in the EaS, followed by a partial fall in the MiS, as well as at the EnS. Findings from this study can most likely be explained by the possible changes in anthropometric traits, heart and lung mass, and nervous system maturation associated with normal growth and maturation leading to improvement in aerobic capacity (Malina et al., 2015; Clemente et al., 2020a).

A significant improvement in LSS from EaS to MiS was reported, but greater levels were evident from MiS to EnS and were significantly related to PHV. A differential timing of adolescent spurts body dimensions, functional capacities, and sports skills during youth has been reported (Malina et al., 2015). Indeed, for example, data pertaining to speed suggest peak gains occur before PHV in boys (Beunen et al., 1988), whereas tests of strength and power achieved peak gains after PHV (Beunen and Malina, 1988; Malina et al., 2004), and peak velocity in $VO_{2max}$ was reached concomitant with PHV in boys and girls (Mirwald and Bailey, 1987; Geithern et al., 2004). Thus,

considering the strong relationship between PHV and speed variables (e.g., maximal sprints) (Rommers et al., 2019), it is crucial to monitor PHV in order to explore sensitive periods of trainability that allow individualized training (Balyi and Hamilton, 2004).

Heart rate has been considered a useful predictor of aerobic performance changes over a soccer season (Krustrup et al., 2003; Impellizzeri et al., 2005). The $HR_{max}$ (EaS-EnS period) highlighted improvements mainly from EaS to MiS and was related to $VO_{2max}$. In accordance with our results, significant effects ($p < 0.05$) on $HR_{max}$ percentage in elite soccer players, after the EaS, were found by Kalapotharakos et al. (2011), whereas with reference to $HR_{rest}$, it was reported to develop over the season with higher relation between EaS-EnS and MiS-EnS. The present study appears to confirm the close inverse relationship between physical fitness and HR reported in previous studies (Kingwell and Jennings, 1993; Rabbia et al., 2002). Previous research has also reported significant positive changes in HR after 4 weeks EaS, albeit in futsal players (Soares-Caldeira et al., 2014).

In the analysis of speed variables, a gradual improvement in LSS performance over the season was found and was highly related to the $VO_{2max}$ (EaS-MIS), as well as to the LSS (EaS-MiS and MiS-EnS). These results were also concordant with Dragijsky et al. (2017), who found significant improvements in running speed performance at the end of the season compared to the early season in young elite soccer players. The observed improvements in the present study may be attributed to the shuttle run and small-sided soccer games (anaerobic speed training with required change
direction and LSS) implemented in the EaS. With regard to the accumulated load, a strong association between MALP (MiS-EnS) and $HR_{max}$ (EaS-EnS) was found, indicating that the intermediate period of the season may be more positive for optimizing $HR_{max}$. To the best of our knowledge, only one study (Clemente et al., 2020a) has analyzed the relationship between cardiorespiratory fitness and the accumulated training load in young soccer players. These authors found a large correlation between training load and improvements in aerobic power over time. The strong association obtained in the present study may be attributed to the natural morphological changes occurring during growth spurt, leading to the development of players' motor skills. The results of the present study highlight the effects of training load and physiological changes over the season on performance, which should be considered by practitioners and researchers in their future work.

Predictive models are now ubiquitous in sports sciences research; nevertheless, most are developed to justify acute performance responses based on biomechanical and anthropometric measures (Barbosa et al., 2010; Morais et al., 2016) or even to predict aerobic capacity or maximal strength (Hoffman and Kang, 2003). Only one study, to the authors' knowledge, provided indirect predictive models of the responses to different training program interventions in prepubertal children (Alves et al., 2020). Interestingly, the present study sought to understand and predict responses on fitness levels considering accumulated load, baseline levels, and PHV. First, analysis of $VO_{2max}$ showed significant results, where it was possible to predict the aerobic power levels of young soccer players during the season. Additionally, it was also found that both baseline $VO_{2max}$ and $HR_{rest}$ were significant predictors of aerobic power levels. Contrarily, in the speed and physiological variables, respectively, there were no significant results. In fact, in LSS and LSM, poor predictors of speed were reported. It is conceivable that predictability would be more attainable if additional measurements were included, such as motor control, psychological aspects, or behavior genetics (Vandendriessche et al., 2011; Schutte et al., 2016). Nevertheless, despite the novelty of the present study, there are some limitations that should be addressed: (i) the number of included players was quite low, where adding more teams may have obtained more representative sample of the soccer young population; however, a priori power calculations indicated sufficient numbers were recruited, whereas post hoc power calculations indicated the achieved power; and (ii) the absence of psychological component evaluation, which could influence the motivation and performance.

## CONCLUSION

Aerobic and speed performance improved over the season and had a strong association with PHV. Moreover, the physiological

variables demonstrated improvements from the EaS to EnS; however, improvements were more evident in the intermediate period, from EaS to MiS, and significantly associated with $VO_{2max}$. Thus, it is advisable that practitioners and researchers consider such changes in future work. Additionally, through multiple linear regression analysis, we were able to discern significant predictors of aerobic power levels of young soccer players based on their $VO_{2max}$ baseline and $HR_{rest}$, whereas for speed and physiological variables, no significant results were found. Therefore, the results of the present study highlight the effects of training load and physiological changes over the season on performance, which should be considered by practitioners and researchers in their future work.

## DATA AVAILABILITY STATEMENT

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

## ETHICS STATEMENT

The study was conducted in accordance with the latest version of the Declaration of Helsinki and approved by the Ethics Committee of the University of Isfahan with the ethical code number IR.UI.REC.1397.181. Prior to study commencement, all players and their parents/guardians signed a written informed consent form. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

## AUTHOR CONTRIBUTIONS

HN, FC, and JP-G designed the study and drafted the manuscript. HN performed the experiments. HN, AA, CC, UG, and HZ participated in the data analysis and drafted the manuscript. HN, FC, UG, HZ, and JP-G critically revised the manuscript. All authors read and approved the final version of the manuscript.

## FUNDING

The authors acknowledge the support of the Deutsche Forschungsgemeinschaft (DFG) and Open Access Publishing Fund of the University of Potsdam, Germany.

## ACKNOWLEDGMENTS

We would like to thank all the players and the staff from Sepahan F.C.

## REFERENCES

* Álvarez-Kurogi, L., Onetti, W., Fernández-García, J. C., and Castillo-Rodríguez, A. (2019). Does the psychological profile influence the position of promising

young futsal players? *PLoS One* 14:e0224326. doi: 10.1371/journal.pone.0224326

* Alves, A. R., Marta, C., Neiva, H. P., Nunes, C., Izquierdo, M., and Marques, M. C. (2020). Modeling fitness variable responses to training in
prepubescent children. *J. Strength Cond. Res.* 34, 2352–2359. doi: 10.1519/JSC. 0000000000002201

Arazi, H., Mirzaei, B., and Nobari, H. (2015). Anthropometric profile, body composition and somatotyping of national Iranian cross-country runners. *Turk. J. Sports Exerc.* 17:35. doi: 10.15314/tjse.49873

Balyi, I. and Hamilton, A. (2004). *Long-Term Athlete Development: Trainability in Childhood and Adolescence—Windows of Opportunity, Optimal Trainability*. National Coaching Institute British Columbia and Advanced Training and Performance Ltd.

Barbosa, T. M., Costa, M., Marinho, D. A., Coelho, J., Moreira, M., and Silva, A. J. (2010). Modeling the links between young swimmers' performance: energetic and biomechanic profiles. *Pediatr. Exerc. Sci.* 22, 379–391.

Barker, A. R., and Armstrong, N. (2010). "Exercise testing elite young athletes," in *The Elite Young Athlete*, eds N. Armstrong and A. M. McManus (Basel: KARGER), 106–125. doi: 10.1159/000320642

Barnes, C., Archer, D., Hogg, B., Bush, M., and Bradley, P. (2014). The evolution of physical and technical performance parameters in the english premier league. *Int. J. Sports Med.* 35, 1095–1100. doi: 10.1055/s-0034-1375695

Batterham, A. M., and Hopkins, W. G. (2006). Making meaningful inferences about magnitudes. *Int. J. Sports Physiol. Perform.* 1, 50–57. doi: 10.1123/ijspp.1.1.50

Baumgartner, T. A., and Chung, H. (2001). Confidence limits for intraclass reliability coefficients. *Meas. Phys. Educ. Exerc. Sci.* 5, 179–188. doi: 10.1207/S15327841MPEE0503_4

Beunen, G., and Malina, R. M. (1988). Growth and physical performance relative to the timing of the adolescent spurt. *Exerc. Sports Sci. Rev.* 16, 503–540.

Beunen, G. P., Malina, R. M., Van 't Hof, M. A., Simons, J., Ostyn, M., Renson, R., et al. (1988). *Adolescent Growth and Motor Performance: a Longitudinal Study of Belgian Boys*. Champaign, IL: HUman Kinetics.

Bloomfield, J., Polman, R., and O'Donoghue, P. (2007). Physical demands of different positions in FA premier league soccer. *J. Sports Sci. Med.* 6, 63–70.

Bradley, P. S., Dellal, A., Mohr, M., Castellano, J., and Wilkie, A. (2014). Gender differences in match performance characteristics of soccer players competing in the UEFA champions league. *Hum. Mov. Sci.* 33, 159–171. doi: 10.1016/j.humov.2013.07.024

Buchheit, M., Mendez-Villanueva, A., Delhomel, G., Brughelli, M., and Ahmaidi, S. (2010). Improving repeated sprint ability in young elite soccer players: repeated shuttle sprints vs. explosive strength training. *J. Strength Cond. Res.* 24, 2715–2722. doi: 10.1519/JSC.0b013e3181bf0223

Caldwell, B. P., and Peters, D. M. (2009). Seasonal variation in physiological fitness of a semiprofessional soccer team. *J. Strength Cond. Res.* 23, 1370–1377. doi: 10.1519/JSC.0b013e3181a4e82f

Christopher, J. Beato, M., and Hulton, A. T. (2016). Manipulation of exercise to rest ratio within set duration on physical and technical outcomes during small-sided games in elite youth soccer players. *Hum. Mov. Sci.* 48, 1–6. doi: 10.1016/j.humov.2016.03.013

Clemente, F. M., Nikolaidis, P. T., Rosemann, T., and Knechtle, B. (2019). Dose-response relationship between external load variables, body composition, and fitness variables in professional soccer players. *Front. Physiol.* 10:443. doi: 10.3389/fphys.2019.00443

Clemente, F. M., Silva, A. F., Alves, A. R., Nikolaidis, P. T., Ramirez-Campillo, R., Lima, R., et al. (2020a). Variations of estimated maximal aerobic speed in children soccer players and its associations with the accumulated training load: comparisons between non, low and high responders. *Physiol. Behav.* 224:113030. doi: 10.1016/j.physbeh.2020.113030

Clemente, F. M., Silva, A. F., Clark, C. C. T., Conte, D., Ribeiro, J., Mendes, B., et al. (2020b). Analyzing the seasonal changes and relationships in training load and wellness in elite volleyball players. *Int. J. Sports Physiol. Perform.* 15, 731–740. doi: 10.1123/ijspp.2019-0251

Coratella, G., Beato, M., and Schena, F. (2016). The specificity of the Loughborough intermittent shuttle test for recreational soccer players is independent of their intermittent running ability. *Res. Sports Med.* 24, 363–374. doi: 10.1080/15438627.2016.1222279

Dalen, T., Jørgen, I., Gertjan, E., Geir Havard, H., and Ulrik, W. (2016). Player load, acceleration, and deceleration during forty-five competitive matches of elite soccer. *J. Strength Cond. Res.* 30, 351–359. doi: 10.1519/JSC.0000000000001063

Dragijsky, M., Aly, T., Zahalka, F., Kunzmann, E., and Hank, M. (2017). Seasonal variation of agility, speed and endurance performance in young elite soccer players. *Sports* 5:12. doi: 10.3390/sports5010012

Fitzpatrick, J. F., Hicks, K. M., and Hayes, P. R. (2018). Dose-response relationship between training load and changes in aerobic fitness in professional youth soccer players. *Int. J. Sports Physiol. Perform.* 13, 1365–1370. doi: 10.1123/ijspp.2017-0843

Foster, C., Florhaug, J. A., Franklin, J., Gottschall, L., Hrovatin, L. A., Parker, S., et al. (2001). A new approach to monitoring exercise training. *J. Strength Cond. Res.* 15, 109–115.

Geithern, C. A., Thomis, M. A., Eynde Vanden, B., Maes, H. H. M., Loos, R. J. F., Peeters, M., et al. (2004). Growth in peak aerobic power during adolescence. *Med. Sci. Sports Exerc.* 36, 1616–1624. doi: 10.1249/01.MSS.0000139807.72229.41

Haritonidis, K., Koutlianos, N., Kouidi, E., Haritonidou, M., and Deligiannis, A. (2004). Seasonal variation of aerobic capacity in elite soccer, basketball and volleyball players. *J. Hum. Mov. Stud.* 46, 289–302.

Hoffman, J. R., and Kang, J. (2003). Strength changes during an in-season resistance-training program for football. *J. Strength Cond. Res.* 17:109. doi: 10.1519/1533-42872003017<0109:SCDAIS>2.0.CO;2

Impellizzeri, F. M., Mognoni, P., Sassi, A., and Rampinini, E. (2005). "Validity of a submaximal running test to evaluate aerobic fitness changes in soccer players," in *The Proceedings of the Science and Football V. the Proceedings of the Fifth World Congress on Science and Football* (New York, NY: Routledge).

Jenkins, D. G., and Quintana-Ascencio, P. F. (2020). A solution to minimum sample size for regressions. *PLoS One* 15:e0229345. doi: 10.1371/journal.pone.0229345

Kalapotharakos, V. I., Ziogas, G., and Tokmakidis, S. P. (2011). Seasonal aerobic performance variations in elite soccer players. *J. Strength Cond. Res.* 25, 1502–1507. doi: 10.1519/JSC.0b013e3181da85a9

Kingwell, B. A., and Jennings, G. L. (1993). Effects of walking and other exercise programs upon blood pressure in normal subjects. *Med. J. Aust.* 158, 234–238.

Krustrup, P., Mohr, M., Amstrup, T., Rysgaard, T., Johansen, J., Steensberg, A., et al. (2003). The yo-yo intermittent recovery test: physiological response. *Reliability, and Validity. Med. Sci. Sports Exerc.* 35, 697–705. doi: 10.1249/01.MSS.0000058441.94520.32

Krustrup, P., Mohr, M., Ellingsgaard, H., and Bangsbo, J. (2005). Physical demands during an elite female soccer game: importance of training status. *Med. Sci. Sports Exerc.* 37, 1242–1248. doi: 10.1249/01.mss.0000170062.73981.94

Lloyd, R. S., and Oliver, J. L. (2012). The youth physical development model. *Strength Cond. J.* 34, 61–72. doi: 10.1519/SSC.0b013e31825760ea

Loturco, I., Pereira, L. A., Kobal, R., Zanetti, V., Gil, S., Kitamura, K., et al. (2015). Half-squat or jump squat training under optimum power load conditions to counteract power and speed decrements in Brazilian elite soccer players during the preseason. *J. Sports Sci.* 33, 1283–1292. doi: 10.1080/02640414.2015.1022574

Malina, R. M., Bouchard, C., and Oded, B. O. (2004). *Growth, Maturation, and Physical Activity*. Champaign, IL: Human Kinetics.

Malina, R. M., Rogol, A. D., Cumming, S. P., Coelho e Silva, M. J., and Figueiredo, A. J. (2015). Biological maturation of youth athletes: assessment and implications. *Br. J. Sports Med.* 49, 852–859. doi: 10.1136/bjsports-2015-094623

McArdle, W., Katch, F. I., and Katch, V. L., (2011). *Essentials of Exercise Physiology*. Philadelphia: Lippincott Williams and Wilkins.

Medicine ACoS. (2013). *ACSM's Guidelines for Exercise Testing and Prescription*, 10th Edn. Philadelphia, PA: Lippincott Williams & Wilkins.

Mirkov, D., Nedeljkovic, A., Kukolj, M., Ugarkovic, D., and Jaric, S. (2008). Evaluation of the reliability of soccer-specific field tests. *J. Strength Cond. Res.* 22, 1046–1050. doi: 10.1519/JSC.0b013e31816eb4af

Mirwald, R. L., and Bailey, D. A. (1987). Maximal aerobic power—a longitudinal analysis. *Br. J. Sports Med.* 21:140.

Mirwald, R. L., Baxter-Jones, A. D. G., Bailey, D. A., and Beunen, G. P. (2002). An assessment of maturity from anthropometric measurements. *Med. Sci. Sports Exerc.* 34, 689–694. doi: 10.1097/00005768-200204000-00020

Mohr, M., Krustrup, P., and Bangsbo, J. (2005). Fatigue in soccer: a brief review. *J. Sports Sci.* 23, 593–599. doi: 10.1080/0264041000021286

Morais, J. E., Silva, A. J., Marinho, D. A., Marques, M. C., Batalha, N., and Barbosa, T. M. (2016). Modelling the relationship between biomechanics and performance of young sprinting swimmers. *Eur. J. Sports Sci.* 16, 661–668. doi: 10.1080/17461391.2016.1149227

Nobari, H., Aquino, R., Clemente, F. M., Khalafi, M., Adsuar, J. C., and Pérez-Gómez, J. (2020a). Description of acute and chronic load, training monotony
and strain over a season and its relationships with well-being status: a study in elite under-16 soccer players. *Physiol. Behav.* 225:113117. doi: 10.1016/j.physbeh.2020.113117

Nobari, H., Tubagi Polito, L. F., Clemente, F. M., Pérez-Gómez, J., Ahmadi, M., Garcia-Gordillo, M. Á, et al. (2020b). Relationships between training workload parameters with variations in anaerobic power and change of direction status in elite youth soccer players. *Int. J. Environ. Res. Public Health* 17:7934. doi: 10.3390/ijerph17217934

Noon, M. R., James, R. S., Clarke, N. D., Akubat, I., and Thake, C. D. (2015). Perceptions of well-being and physical performance in English elite youth footballers across a season. *J. Sports Sci.* 33, 2106–2115. doi: 10.1080/02640414.2015.1081393

Owen, A. L., Wong, D. P., Paul, D., and Dellal, A. (2012). Effects of a periodized small-sided game training intervention on physical performance in elite professional soccer. *J. Strength Cond. Res.* 26, 2748–2754. doi: 10.1519/JSC.0b013e318242d2d1

Rabbia, F., Grosso, T., Cat Genova, G., Conterno, A., De Vito, B., Mulatero, P., et al. (2002). Assessing resting heart rate in adolescents: determinants and correlates. *J. Hum. Hypertens.* 16, 327–332. doi: 10.1038/sj.jhh.1001398

Rahmat, A. J., Arsalan, D., Bahman, M., and Hadi, N. (2016). Anthropometrical profile and bio-motor abilities of young elite wrestlers. *Phys. Educ. Students* 20, 63–69. doi: 10.15561/20755279.2016.0608

Rommers, N., Mostaert, M., Goossens, L., Vaeyens, R., Witvrouw, E., Lenoir, M., et al. (2019). Age and maturity related differences in motor coordination among male elite youth soccer players. *J. Sports Sci.* 37, 196–203. doi: 10.1080/02640414.2018.1488454

Rumpf, M. C., Cronin, J. B., Oliver, J. L., and Hughes, M. G. (2013). Vertical and leg stiffness and stretch-shortening cycle changes across maturation during maximal sprint running. *Hum. Mov. Sci.* 32, 668–676. doi: 10.1016/j.humov.2013.01.006

Schutte, N. M., Nederend, I., Hudziak, J. J., de Geus, E. J. C., and Bartels, M. (2016). Differences in adolescent physical fitness: a multivariate approach and meta-analysis. *Behav. Genet.* 46, 217–227. doi: 10.1007/s10519-015-9754-2

Slaughter, M. H., Lohman, T. G., Boileau, R., Horswill, C. A., Stillman, R. J., Van Loan, M. D., et al. (1988). Skinfold equations for estimation of body fatness in children and youth. *Hum. Biol.* 60, 709–723.

Soares-Caldeira, L. F., de Souza, E. A., de Freitas, V. H., de Moraes, S. M. F., Leicht, A. S., and Nakamura, F. Y. (2014). Effects of additional repeated sprint training during preseason on performance, heart rate variability, and stress symptoms in futsal players. *J. Strength Cond. Res.* 28, 2815–2826. doi: 10.1519/JSC.000000000000461

Vandendriessche, J. B., Vandorpe, B., Coelho-e-Silva, M. J., Vaeyens, R., Lenoir, M., Lefevre, J., et al. (2011). Multivariate association among morphology, fitness, and motor coordination characteristics in boys age 7 to 11. *Pediatr. Exerc. Sci.* 23, 504–520. doi: 10.1123/pes.23.4.504

Zamparo, P., Bolomini, F., Nardello, F., and Beato, M. (2015). Energetics (and kinematics) of short shuttle runs. *Eur. J. Appl. Physiol.* 115, 1985–1994. doi: 10.1007/s00421-015-3180-2

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Nobari, Alves, Clemente, Pérez-Gómez, Clark, Granacher and Zouhal. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.