Introduction

Soccer necessitates high-intensity efforts interspersed with low-intensity activities (Gil et al., 2007). The escalating demands, intensity and speed of the game require soccer players to have highly developed aerobic and anaerobic mechanisms (Chena et al., 2022; Stølen et al., 2005). This demand for a high level of aerobic capacity is a key factor differentiating among elite soccer players (Slimani et al., 2019).

Elevated aerobic capacity enables the maintenance of peak performance throughout the season and expedites recovery following intense efforts. In developing athletes, it permits tolerance of numerous training units (Deprez et al., 2015), facilitating greater high-intensity sprint and run frequencies as well as increased ball contact (Buchheit et al., 2010; Clemente et al., 2020).

Furthermore, superior aerobic capacity correlates with enhanced cognitive function, including internal inhibition, attention, and memory—crucial for skill acquisition (Scudder et al., 2014). As indicated in studies on the subject, better-performing young athletes exhibit improved reaction times, selective attention, executive control and accuracy (Pontifex et al., 2011; Voss et al., 2011). Therefore, monitoring aerobic capacity, even in early training stages, is crucial. This emphasis on physiological development, alongside technical skill development, is paramount in optimising training efficacy among youth soccer players.

Comprehensive assessment of aerobic capacity, in addition to evaluating the maximum values of a given indicator, should include evaluation of variables at the anaerobic threshold level. Maximal aerobic speed (MAS) should also be considered, a variable often overlooked in the limited research on this developmental period (Clemente et al., 2020; Mendez-Villanueva et al., 2010). This is necessary to better understand the relationships among ontogenetic development, physiological indices and regular sports participation in childhood.

For such assessments, MAS is a key variable, as it determines the threshold running speed beyond which uncompensated metabolic acidosis begins, making it a good indicator of the fitness level and conditioning (Baker and Heaney, 2015). However, there are few studies on the effects of training loads among children in the first stages of training, or the effects of training on changes in individual aerobic capacity indices.

Therefore, the aim of this study was to determine developmental changes in aerobic capacity of boys aged 7–12 participating in soccer training. An additional objective was to assess the influence of applied training loads on changes in aerobic capacity.

Methods

Participants

The study involved 54 prepubertal boys from a Polish soccer academy, selected using identical criteria, including fitness tests and coach evaluations. The boys were divided into two groups: Group 1 included 29 boys aged 7, born in the same calendar year, beginning their training, while Group 2 comprised 25 boys aged 10, with three years of training experience. Complete data were collected for 42 boys (n = 21 per group). All participants had high physical activity levels, as assessed by parent interviews and questionnaires. Training regimens differed between groups and evolved with age, reflecting standard practices within the academy. The 7-year-olds trained for 12 weeks (September to December) with two weekly 60-min sessions, but did not play league matches. Boys aged 8–12 trained for 35–36 weeks (January to December), playing matches in the local junior league I, according to their age category. The 8-year-olds had two weekly 60-min sessions and 20 matches annually, 9-year-olds had two weekly 75-min sessions and 22 matches annually and one 7-day summer camp, while 10-year-olds had three weekly 75-min sessions and 22–26 matches (60 min each), alongside a 10-day summer camp. The 11- and 12-year-olds had three weekly 75-min sessions and 24–29 matches annually (70 min each), alongside a 12-day summer camp.

Inclusion criteria were as follows: physician clearance for exercise testing (including medical history, electrocardiogram, and blood and urine analysis), consistent training and match participation, and written parental consent. Exclusion criteria encompassed illness, injury, or training interruptions exceeding one week. Data collection spanned three years for Group 1 (ages 7–10) and two years for Group 2 (ages 10–12). Exclusion primarily resulted from illness, injury, or physician-mandated withdrawal. Ethical approval was obtained from the Bioethics Committee of the Regional Medical Chamber in Cracow, Cracow, Poland (protocol code: 16/KBL/OIL/2006; approval date: 15 March 2006), and all testing procedures were conducted under medical supervision.

The average annual training volume varied due to attendance, which decreased with age because of injuries or other activities.

Measures

Anthropometric measurements, including body mass and tissue composition, were obtained using a TANITA Body Composition Analyzer (Tanita Company, model TBF-300, Tokyo, Japan). Body height was estimated via a standard anthropometer Seca 213 (Seca gmbh and co. kg., model 213, Hamburg, Germany). Testing was performed on a mechanical treadmill (Cardionics, type 2113, Jönköping, Sweden). Respiratory variables were measured using a computerized ergospirometer (Medikro Oy, type 919, Kuopio, Finland). The "Polar-Elektro S-610i" type device (Polar Electro Oy, type S-610i, Helsinki, Finland) was applied to measure the heart rate during exercise testing.

Measurements were conducted annually in November/December (transition season) and included aerobic capacity indices as well as somatic characteristics. The study comprised three phases: familiarisation, during which participants were acquainted with the exercise test procedures; general measurements, conducted to record somatic indices; and aerobic capacity testing. All measurements were conducted by two experienced laboratory technicians in a controlled laboratory environment (temperature: 22°C, atmospheric pressure: 999 mBar). The study involved a three-year observation of Group 1 (ages 7–10) and a two-year observation of Group 2 (ages 10–12).

Design and Procedures

Aerobic capacity indices were measured under laboratory conditions using a graded exercise test performed on a mechanical treadmill, to determine maximum oxygen uptake (VO2max) and the anaerobic threshold (AT). Participants began with a four-minute warm-up at a speed of 6 km/h. The load was then increased every 2 min by elevating the treadmill speed by 1 km/h. The treadmill inclination angle remained constant at 0° throughout the test. Participants continued the test until volitional exhaustion.

Achievement of VO2max was confirmed by the following criteria:

  • attainment of the maximum heart rate (HRmax),

  • a plateau or a decrease in oxygen uptake (VO2) at the end of exercise despite continuing to run and a further increase in load,

  • a respiratory exchange ratio above 1.0.

The anaerobic threshold (AT) was determined via change dynamics in ventilatory variables, specifically: maximum value of the carbon dioxide percentage in exhaled air (tCO2), a significant increase in minute ventilation (VE), and the minimum value of ventilatory equivalent for carbon dioxide (VE/VCO2).

Each training session was documented on a training load sheet, which allowed to comprehensively assess work time across five metabolic intensity zones: aerobic-sustained, aerobic, mixed, anaerobic-glycolytic and anaerobic-phosphagen.

In addition, to increase the accuracy of training load evaluation and ascertain the actual physiological responses to the training loads, the heart rate was measured using the Polar Electro S-610i sport-tester. Training load analysis, including the determination of work time in particular metabolic zones, was performed using the Polar Precision Performance SW 4.0 program (Polar Electro Oy, Polar Precision Performance SW 4.0, Helsinki, Finland).

Statistical Analysis

All data were statistically processed using Statistica 13.1 software (Tibco Software Inc., Statistica 13.1, Santa Clara, California, USA). Descriptive statistics were calculated for all variables, including arithmetic means (x̄) and standard deviations. Repeated measures analysis of variance with Greenhouse-Geisser correction was applied to determine the significance of differences between groups, with the Sidak correction applied for multiple comparisons within groups. Normality of distribution for all variables was confirmed using the Shapiro-Wilk test. An alpha level of p < 0.05 was considered statistically significant for all analyses. The Student's t-test for independent samples was used to determine the significance of differences observed between Groups 1 and 2. Pearson's r correlation analysis was implemented to assess the relationship between measured variables and training loads within the group. The strength of correlations was interpreted using the following criteria: <0.3 (or 0 to −0.3) biologically negligible; 0.31–0.5 (or −0.31 to −0.5) weak correlation; 0.51–0.7 (or −0.51 and −0.7) moderate correlation; 0.71–0.9 (or −0.71 to 0.9) strong correlation and >0.9 (or < −0.9) very strong correlation (Miot, 2018).

A priori power analysis was performed utilizing G*Power 3.1 software (Dusseldorf, Germany) to determine the required sample size for the study (Faul et al., 2007). Considering the experiment's structure of two groups and eight measurement points, alongside a predefined alpha of 0.05 and a test power of 0.95, it was calculated that a total of 36 subjects (thirteen per group) would be required to adequately power the study.

Results

The athletes exhibited a characteristic increase in body height and mass across the studied age period. These differences were statistically significant in all age groups. A similar trend was observed for lean body mass, with significant gains differentiating among all age groups. The percentage of body fat (%F) decreased with age, although these differences were not statistically significant. Group-related differences in fat mass were significant between 7- and 9-year-olds, 7- and 10-year-olds, as well as 8- and 9-year-olds. Significant differences were also noted in Group 2 between 10- and 11-year-olds, as well as 10- and 12-year-olds.

Analysis of birthdates revealed over-representation of athletes born in the first two quarters of the year in both groups (Table 1).

Table 1

Characteristics of subjects and aerobic capacity indices achieved in the graded test.

Distribution according to the players’ quarter of the birth date in Groups 1 and 2
IIIIIIIV
Gr. 171040
Gr. 27842
Anthropometric variables
Age (years)BH (cm)BM (kg)F (%)FM (kg)LBM (kg)
Gr. 17128.3 ± 5.125.1 ± 2.113.7 ± 2.13.5 ± 0.621.6 ± 1.9
8134.1 ± 5.027.9 ± 2.613.6 ± 2.23.8 ± 0.824.0 ± 2.1
9139.0 ± 6.130.9 ± 3.413.5 ± 3.14.3 ± 1.326.6 ± 2.5
10145.3 ± 5.634.2 ± 4.012.9 ± 3.34.5 ± 1.629.7 ± 2.8
Gr. 210144.1 ± 4.433.6 ± 4.211.4 ± 4.54.0 ± 2.129.6 ± 2.6
11149.0 ± 4.837.6 ± 4.111.6 ± 3.84.5 ± 1.933.1 ± 2.9
12156.0 ± 5.742.9 ± 5.010.7 ± 3.34.7 ± 1.938.3 ± 3.9
Maximum values of aerobic capacity indices
Age (years)t (min)V (km/h)VO2max (ml/kg/min)VO2max (l/min)VEmax (l/min)HRmax (beat/min)
Gr. 1714.1 ± 2.311.6 ± 1.246.2 ± 4.01.19 ± 0.1452.0 ± 8.1203 ± 7.0
816.3 ± 1.912.5 ± 0.959.1 ± 4.51.68 ± 0.1559.8 ± 8.1206 ± 7.8
916.5 ± 1.912.8 ± 0.963.4 ± 6.81.98 ± 0.2863.0 ± 11.5205 ± 5.3
1017.5 ± 1.713.1 ± 0.859.5 ± 4.22.07 ± 0.1770.1 ± 7.2204 ± 7.0
Gr. 21016.8 ± 2.612.8 ± 1.353.1 ± 8.51.77 ± 0.2573.1 ± 13.7206 ± 10.2
1118.8 ± 2.113.8 ± 1.061.5 ± 5.12.29 ± 0.2379.4 ± 9.2207 ± 8.9
1218.8 ± 2.013.8 ± 1.062.7 ± 5.72.65 ± 0.3589.3 ± 12.7206 ± 8.4
Values obtained during exercise intensity corresponding to anaerobic threshold (AT)
Age (years)V/AT (km/h)VO2/AT (ml/kg/min)VO2/AT (l/min)VE/AT (l/min)HR/AT (beat/min)%HRmax%VO2max
Gr. 178.7 ± 1.237.4 ± 3.70.96 ± 0.1134.4 ± 4.9180 ± 8.288.5 ± 3.781.0 ± 5.3
89.1 ± 1.143.3 ± 3.81.23 ± 0.1237.3 ± 5.0178 ± 7.386.3 ± 3.373.4 ± 6.0
99.1 ± 1.244.8 ± 7.01.38 ± 0.2636.3 ± 8.2176 ± 5.885.8 ± 3.270.6 ± 7.9
109.5 ± 1.142.9 ± 4.51.48 ± 0.1941.3 ± 6.1176 ± 6.686.4 ± 2.972.2 ± 5.7
Gr. 2108.5 ± 1.538.9 ± 6.61.30 ± 0.2242.8 ± 9.5177 ± 11.786.2 ± 3.973.4 ± 7.5
119.0 ± 1.241.5 ± 5.31.56 ± 0.2442.3 ± 7.5175 ± 11.984.6 ± 4.567.6 ± 7.0
129.2 ± 0.942.8 ± 4.61.83 ± 0.3148.4 ± 8.2177 ± 9.386.0 ± 4.268.4 ± 6.9

[i] BH: body height, BM: body mass, F: percentage of body fat, FM: fat mass, LBM: lean body mass, Gr. 1, Gr. 2: Group 1, Group 2, t: duration of exercise, V: running speed at the end of the graded test, VO2max: maximum oxygen uptake, VEmax: maximum minute ventilation, HRmax: maximum heart rate, V/AT: running speed at the anaerobic threshold, VO2/AT: oxygen uptake at the anaerobic threshold, VE/AT: minute ventilation at the anaerobic threshold, HR/AT: heart rate at the anaerobic threshold, %HRmax: percentage of the maximum heart rate at the anaerobic threshold, %VO2max: percentage of maximum oxygen uptake at the anaerobic threshold

Aerobic capacity indices measured during maximal exercise significantly differed among the age groups in most cases (Tables 2 and 3). The most pronounced differences were observed in Group 1 within the initial two years of training.

Table 2

Differences between measurements of exercise duration (t), running speed (V) and maximum oxygen uptake (VO2max: in values relative to subjects’ body mass) obtained during the graded test.

Age (years)t (min)V (km/h)VO2max (ml/kg/min)
Analysis of variance for repeated measures with Greenhouse-Geisser correction
dfB; dfWFpdfB; dfWFpdfB; dfWFp
7–103; 6019.32<0.0013; 6014.89<0.0013; 6063.67<0.001
10–122; 4024.88<0.0012; 4017.47<0.0011.34; 26.7521.56<0.001
Differences between measurements (multiple comparisons with Sidak correction)
Standard errorpStandard errorpStandard errorp
7–80.490.0010.240.0051.06<0.001
7–90.580.0030.310.0071.46<0.001
7–100.50<0.0010.27<0.0011.48<0.001
8–90.380.9900.190.7981.260.015
8–100.430.0670.200.0611.13<0.001
9–100.360.0870.190.4311.520.980
10–110.25<0.0010.18<0.0011.830.001
10–120.36<0.0010.21<0.0011.88<0.001
11–120.340.9990.190.9920.870.450
Differences between age groups: 1 (7–10) and 2 (10–12 years), cross-sectional equations and the Student's t-test
tptptp
7–11−−6.91<0.001−6.29<0.001−10.8<0.001
7–12−7.00<0.001−6.46<0.001−10.91<0.001
8–11−4.08<0.0014.17<0.001−1.630.111
8–12−4.11<0.001−4.36<0.001−2.300.027
9–11−3.670.001−3.340.0021.030.309
9–12−3.700.001−3.530.0010.360.718

[i] t: duration of exercise, V: running speed at the end of the graded test, VO2max: maximum oxygen uptake (in values relative to subjects’ body mass), p level: significance of differences when p < 0.05

Table 3

Differences between measurements of maximum oxygen uptake (VO2max: in global values), maximum minute ventilation (VEmax) and the maximum heart rate (HRmax) obtained during the graded test.

Age (years)VO2max (l/min)VEmax (l/min)HRmax (beat/min)
Analysis of variance for repeated measures with Greenhouse-Geisser correction
dfB; dfWFpdfB; dfWFpdfB; dfWFp
7–103; 60177.41< 0.0013; 6063.67<0.0012.26; 45.291.470.233
10–122; 40101.94< 0.0011.34; 26.7521.56<0.0011.51; 30.191.260.289
Differences between measurements (multiple comparisons with Sidak correction)
Standard errorpStandard errorpStandard errorp
7–80.03<0.0011.06<0.001--
7–90.05<0.0011.46<0.001--
7–100.05<0.0011.48<0.001--
8–90.04<0.0011.260.015--
8–100.03<0.0011.13<0.001--
9–100.050.2791.520.98--
10–110.05<0.0011.830.001--
10–120.07<0.0011.88<0.001--
11–120.06<0.0010.870.450--
Differences between age groups: 1 (7–10) and 2 (10–12 years), cross-sectional equations, Student's t-test
tptptp
7–11−18.92<0.001−10.8<0.001−1.740.090
7–12−17.83<0.001−10.91<0.001−1.090.281
8–11−10.23<0.001−1.630.111−0.480.634
8–12−11.74<0.001−2.300.0270.170.865
9–11−3.96<0.0011.030.309−0.830.414
9–12−6.95<0.0010.360.718−0.090.930

[i] VO2max: maximum oxygen uptake (in global values), VEmax: maximum minute ventilation, HRmax: maximum heart rate, p level: significance of differences when p < 0.05

The results obtained during exercise at anaerobic threshold intensities showed significant differences in oxygen uptake (VO2/AT) across all age groups except for participants between the age of 9 and 10 (Table 4). Significant differences, although not universally observed, were also found in minute ventilation (VE/AT) and the percentage of maximum oxygen uptake (%VO2max) (Table 5).

Table 4

Differences between measurements of exercise duration (t/AT), running speed (V/AT) and oxygen uptake (VO2/AT) relative to body mass and global values obtained during the graded test at exercise intensity corresponding to the anaerobic threshold (AT).

Age (years)t/AT (min)V/AT (km/h)VO2/AT (ml/kg/min)VO2/AT (l/min)
Analysis of variance for repeated measures with Greenhouse-Geisser correction
dfB; dfWFpdfB; dfWFpdfB; dfWFpdfB; dfWFp
7–103.00; 002.500.0683; 603.700.0162.38; 47.5412.45<0.0012.30; 45.9360.78<0.001
10–122; 402.680.0812; 402.450.1002; 403.330.0462; 4038.99<0.001
Differences between measurements (multiple comparisons with Sidak correction)
Stand. errorpStand. errorpStand. errorpStand. errorp
7–8--0.220.2511.03<0.0010.03<0.001
7–9--0.250.4941.42<0.0010.05<0.001
7–10--0.310.0891.450.0060.04<0.001
8–9--0.230.9991.480.8950.050.027
8–10--0.250.7450.930.9990.03<0.001
9–10--0.180.2321.390.7230.040.173
10–11----1.560.2810.05<0.001
10–12----1.710.0970.07<0.001
11–12----1.350.7460.060.001
Differences between age groups: 1 (7–10) and 2 (10–12 years), cross-sectional equations, Student's t-test
tptptptp
7–11−0.790.434−0.800.429−2.930.006−10.45<0.001
7–12−1.230.226−1.620.113−4.14<0.001−11.92<0.001
8–11−0.070.9440.550.5821.240.222−5.67<0.001
8–12−0.470.638−0.150.8780.400.694−8.16<0.001
9–11−0.100.9180.400.6941.720.093−2.300.027
9–12−0.500.623−0.290.7731.120.271−5.03<0.001

[i] t/AT: duration of exercise at the anaerobic threshold, V/AT: running speed at the anaerobic threshold, VO2/AT: oxygen uptake at the anaerobic threshold: values relative to subject's body mass (ml/kg/min) and global values (l/min), Stand. error: standard error, p level: significance of differences when p < 0.05

Table 5

Differences between measurements of minute ventilation (VE/AT), the heart rate (HR/AT), the percentage of the maximum heart rate (%HRmax) and the percentage of maximum oxygen uptake (%VO2max) obtained during the graded test at exercise intensity corresponding to the anaerobic threshold (AT).

Age (years)VE/AT (l/min)HR/AT (beat/min)%HRmax%VO2max
Analysis of variance for repeated measures with Greenhouse-Geisser correction
dfB; dfWFpdfB; dfWFpdfB; dfWFpdfB; dfWFp
7–103; 6011.41<0.0013; 601.390.2543; 602.950.0403; 6014.20<0.001
10–122; 407.820.0012; 400.190.8312; 401.190.3142; 404.630.015
Differences between measurements (multiple comparisons with Sidak correction)
Stand. errorpStand. errorpStand. errorpStand. errorp
7–80.850.018--1.140.3661.36<0.001
7–91.250.609--1.080.1141.96<0.001
7–101.08<0.001--1.070.3051.77<0.001
8–91.560.990--0.880.9871.860.622
8–101.220.021--0.88<0.0011.600.976
9–101.250.004--0.860.9831.770.944
10–111.830.990----1.980.027
10–121.950.026----2.330.129
11–121.390.001----1.790.968
Differences between age groups: 1 (7–10) and 2 (10–12 years); cross-sectional equations and the Student's t-test
tptptptp
7–11−4.06<0.0011.330.1903.060.0046.96<0.001
7–12−7.07<0.0011.040.3032.100.0426.62<0.001
8–11−2.610.0130.800.4291.420.1632.870.007
8–12−5.61<0.0010.410.6860.340.7382.520.016
9–11−2.520.0160.250.8060.950.3491.290.204
9–12−4.99<0.001−0.280.782−0.170.8670.970.336

[i] VE/AT: minute ventilation at the anaerobic threshold, HR/AT: heart rate at the anaerobic threshold, %HRmax: percentage of the maximum heart rate at the anaerobic threshold, %VO2max: percentage of maximum oxygen uptake at the anaerobic threshold, Stand. error: standard error, p level: significance of differences when p < 0.05

Total work time (per metabolic zone) demonstrated a steady increase with age, except for the 10- to 11-year-olds in Group 2. The "intake" group began training in September, with less than three months prior to the first assessment point, meaning the training period lasted only four months. Thus, the lowest training loads were observed in the 7-year-olds. From 9 years onward, annual training loads were influenced not only by regular training sessions but also by sports camps. Pearson's r correlation analysis did not reveal a significant effect of the applied training loads on the level of aerobic capacity indices, with the exception of %HRmax and %VO2max in 8-year-old boys (Table 6).

Table 6

Annual volume of training measures, their structure in particular metabolic zones and analysis of Pearson's r correlation analysis between training loads and selected dependent variables.

Age (years)Operating time per intensity zone (min)Annual training volume (min)
IIIIIIIVVVI
72324496473322-1383
815132060217411793-5957
9172627503038241126-7881
10166829713120263226-8248
10157331694097223154-9216
11173331503831267159-9140
12158137303917169249-9646
Pearson's r correlation analysis between training loads and selected dependent variables
Age (years)Dependent variabler
7V/AT−0.16−0.17−0.21−0.19−0.21-−0.25
VO2max (ml.kg.min-1)0.160.31−0.180−0.07-0.05
%VO2max−0.03−0.01−0.050.05−0.25-−0.04
%HRmax0.0100−0.280.23-0
8V/AT−0.13−0.10−0.21−0.24−0.14-−0.15
VO2max (ml.kg.min-1)0.300.330.210.250.31-0.28
%VO2max−0.52*−0.52*−0.58**−0.62**−0.55**-−0.55**
%HRmax−0.55**−0.55**−0.57**−0.66**−0.56**-−0.57**
9V/AT−0.13−0.09−0.15−0.27−0.12-−0.13
VO2max (ml.kg.min-1)0.130.190.130.030.13-0.15
%VO2max−0.21−0.19−0.22−0.24−0.20-−0.21
%HRmax−0.32−0.29−0.35−0.44*−0.34-−0.32
10V/AT−0.14−0.10−0.09−0.14−0.12-−0.11
VO2max (ml.kg.min-1)−0.0300.01−0.05−0.02-−0.01
%VO2max0.160.180.190.140.17-0.18
%HRmax0.050.010.0200-0.02
10V/AT0.050.120.150.120.09-0.13
VO2max (ml.kg.min-1)0.190.270.270.210.23-0.26
%VO2max0.010.00−0.05−0.02−0.06-−0.02
%HRmax0.01−0.03−0.11−0.08−0.12-−0.06
11V/AT−0.04−0.02−0.05−0.03−0.04-−0.04
VO2max (ml.kg.min-1)0.190.220.170.190.18-0.19
%VO2max−0.10−0.13−0.09−0.10−0.11-−0.11
%HRmax−0.15−0.16−0.13−0.14−0.14-−0.14
12V/AT−0.12−0.09−0.11−0.13−0.10-−0.10
VO2max (ml.kg.min-1)0.260.230.220.220.23-0.23
%VO2max−0.32−0.30−0.33−0.30−0.32-−0.32
%HRmax−0.14−0.12−0.14−0.11−0.14-−0.13

[i] V/AT: running speed at the anaerobic threshold, VO2max: maximum oxygen uptake, %VO2max: percentage of maximum oxygen uptake at the anaerobic threshold, %HRmax: percentage of the maximum heart rate at the anaerobic threshold, r: r-Pearson correlation coefficient <0.3 (or 0 to −0.3) biologically negligible; 0.31–0.5 (or −0.31 to −0.5) weak correlation; 0.51–0.7 (or −0.51 and −0.7) moderate correlation; 0.71–0.9 (or −0.71 to 0.9) strong correlation and >0.9 (or <−0.9) very strong correlation, * p < 0.05, ** p < 0.01

Discussion

In this longitudinal study, dynamics of aerobic capacity indices were examined, relative to training loads. Body composition was also assessed in 7- to 12-year-old prepubertal boys participating in soccer training. The longitudinal design, with repeated measurements performed for the same individuals, allowed to note long-term training effects. Importantly, the inclusion of a group beginning training at age 7 provided a unique opportunity to observe the initial impact of training on physiological adaptations. This is particularly relevant given the limited research on early training adaptations in children and the distinct physiological responses to exercise in this population compared to adults (Clemente et al., 2020).

The findings of the present study allow to demonstrate a gradual and systematic increase in absolute VO2max (l/min) among the participants. This increase, driven by muscle tissue growth and training experience, reflects developmental functional transformations occurring during childhood. Furthermore, relative VO2max (ml/kg/min) also increased progressively from early school years, likely due to the combined effects of regular training and age-related decreases in the body fat percentage. The most pronounced increases in VO2max were observed between the age of 7 and 10, with a 74% increase in absolute VO2max and a 29% increase in relative VO2max. This period highlights the synergistic influence of training and biological maturation on aerobic capacity development. Between 10 and 12 years of age, the increases were less pronounced (50% and 18% for absolute and relative VO2max, respectively).

While the findings regarding increases in relative VO2max align with those from previous reports (Baquet et al., 2003), the timing of the largest annual increments differs from some studies. For instance, it has been suggested that the most significant increases occur between 11 and 13 years of age, coinciding with the pubertal growth spurt (Lemura et al., 1999). In contrast, we observed the largest inter-annual changes between 7 and 8 years of age and the 10–11-year-old age group. This discrepancy underscores the potential for individual variation in the timing of aerobic capacity development and suggests that factors beyond pubertal maturation may influence these adaptations in youth athletes. Further research is needed to elucidate the complex interplay of biological and training-related factors contributing to aerobic capacity development in young athletes.

It is widely accepted that high VO2max levels are crucial for success in soccer (Mohr et al., 2003; Ortiz et al., 2024; Teplan et al., 2012). Elite soccer players in the Spanish Primera Division, for example, exhibit average VO2max of 66.4 ± 7.6 ml/kg/min (Casajús, 2001). The findings of this study allow to suggest that the studied boys possessed an adequate physiological predisposition for soccer. Their VO2max levels exceeded those noted for untrained peers and were comparable to those of trained peers from other studies (Strøyer et al., 2004).

The observed high VO2max level is undoubtedly the cumulative effect of exceptional genetic predispositions regarding aerobic capacity and early training adaptation. This aligns with research confirming the possibility of VO2max exceeding 55 ml/kg/min in children with high levels of physical activity, particularly in endurance and team sports (Aandstad et al., 2006). Our observations are also consistent with the findings of a recent study on the impact of intensive game-based training which is the model used in leading soccer academies on the fitness levels of young soccer players (Jastrzębski et al., 2025).

It is quite interesting that despite a higher frequency of training sessions (Bunc et al., 2002), the 10-year-old soccer players in our study did not demonstrate a strong influence of training on aerobic capacity. This observation is further supported by the lack of correlation between training loads and VO2max levels (both absolute and relative), except between the ages of 7 and 8. These findings suggest that differences in VO2max between groups may be primarily attributed to biological development and innate predisposition, rather than training or pubertal status.

The impact of training on VO2max in pre-adolescent children remains a subject of debate. While in some studies long-term benefits are reported of childhood training on adult aerobic capacity (Baquet et al., 2003; Janz et al., 2002), in others, improvement is demonstrated in VO2max following short-term vigorous training interventions (Baquet et al., 2002; Clemente et al., 2020). However, these findings are not universally supported (Williams et al., 2000).

Given the inconclusive nature of training effects on VO2max, a broader perspective on performance assessment is warranted. Therefore, in our study, we also incorporated the assessment of threshold values, such as MAS. These are known to be highly sensitive to training adaptations and commonly used in training programmes for both adult and junior soccer players. However, research on MAS in children and adolescents remains limited (Baker and Heaney, 2015; Clemente et al., 2020).

Although our study was not focused on the relative age effect (RAE), we observed over-representation of players born in the first two quarters of the year in both groups. This suggests that RAE should be considered during players’ recruitment, selection and evaluation, especially given its prevalence in soccer at various training levels (Helsen et al., 2005; Yagüe et al., 2020).

Monitoring and analysing changes in aerobic capacity from initial training stages is of significant predictive value for talent identification, performance prediction and training optimisation (Coppola and Raiola, 2019; Raiola and D'isanto, 2016). While using it solely for selection has limitations due to the multifaceted nature of athletic development, understanding a young player's aerobic potential, especially VO2max with its strong genetic influence, can be a valuable tool for identifying promising talent (Calandro et al., 2020; Ceruso et al., 2019). This is particularly relevant considering the demands of professional soccer and the importance of appropriately managing training loads (Mohr et al., 2003; Sargienko and Starosta, 2003).

One limitation of our study is the small sample size, which may affect the generalizability of our findings. Moreover, there are no precise methods for measuring training loads in young children. The child's body reactions to the same training stimuli may vary significantly. This requires further research and analysis.

Conclusions

The findings of the present study allow to indicate that between the age of 8 and 12, physical training is not the primary factor influencing aerobic capacity. This is evidenced by the lack of correlations between applied training loads and VO2max levels, as well as no significant impact on the anaerobic threshold. However, monitoring initial VO2max levels remains valuable for identifying individuals with high aerobic potential, often linked to success at elite levels of soccer.

A comprehensive approach is essential, including factors such as MAS, which showed sensitivity to training in this study, highlighting its potential in optimizing training for this age group. Coaches should avoid focusing on increasing aerobic exercise volume and instead prioritize speed, agility, and other high-intensity activities. Developing technical skills through high-intensity exercises is more beneficial during this stage of development.