Seasonal Changes and Relationships in Training Loads, Neuromuscular Performance, and Recovery and Stress State in Competitive Female Soccer Players

Background: The purpose of this study was to examine seasonal changes in training load (TL), neuromuscular performance, subjective recovery, and stress state, and to investigate the relationships between acute and chronic TL and neuromuscular performance in competitive female soccer players.

Methods: Nine competitive female soccer players (20.0 ± 1.7 years; 60.3 ± 6.3 kg; 164.0 ± 5.8 cm) completed the Short Recovery and Stress Scale and the countermovement jump (CMJ) with polyvinyl chloride pipe (CMJ0) and 20 kg barbell (CMJ20) at 2–3 h before 1st match (NC₁), 6th match (NC₂), 9th match (C₁), and 15th match (C₂) of the competitive season. TL included total distance, high-speed running, and PlayerLoad. Acute and chronic TL was calculated by using the average of 2 days (D₂), 7 days (D₇), and 21 days (D₂₁) prior to four different match play.

Results: Significant decreases were found from NC₁ to C₁ in D₇ total distance [p = 0.03, Cohen's effect size (d_z) = 1.40]. D₇ total distance and PlayerLoad significantly decreased from NC to C₁ and C₂ (p = 0.001–0.01, d_z = 1.40–1.72). Significant increases were observed from NC₁ to C₁ in CMJ0 jump height (p = 0.03, d_z = 1.40), (p = 0.021, d_z = 1.44), and peak power (p = 0.03, d_z = 1.32). Significant negative correlations were observed for D₇ total distance and CMJ0 jump height (p = 0.02, r = 0.79) and peak power (p = 0.03, r = 0.71) at C₂, while significant positive correlations were observed at C₁ for D₇ PlayerLoad and CMJ0 jump height (p = 0.02, r = 0.80).

Conclusion: Polyvinyl chloride pipe (CMJ0) jump height and peak power may increase from preseason to the midcompetitive season. Seasonal variations may affect the relationships between D₇ TL and CMJ0 performance.

Introduction

Soccer is an intermittent sport consisting of walking, running, sprinting, and changing of direction, kicking, and heading for 90 min (Stølen et al., 2005). Due to the physical demands of the match, the practitioners (e.g., sports coaches, strength and conditioning coaches, sports scientists) require long-term strategic training plans to gradually accumulate training load (TL) and to develop physical capacity of the players during a pre-season (Stølen et al., 2005). However, it may be difficult to improve the physical capacity of the players during this period for the National Collegiate Athletic Association (NCAA) soccer. The NCAA soccer schedule allows only 2 to 3 weeks of pre-season training followed by 12 to 16 weeks of the competitive season (Sams et al., 2018, 2020; Walker et al., 2019; McFadden et al., 2020). As a result, the NCAA coaching staff often plans higher TLs during pre-season compared to the competitive season (McFadden et al., 2020), which increases the risk of injuries (Agel and Schisel, 2013), psychological stress, and physiological damage (Walker et al., 2019). In addition, the NCAA restricts players from reporting their training to the coaching staff during the summer break (May–July) and evidence indicates that players accumulate limited TL during this period (Sams et al., 2018, 2020; Walker et al., 2019; McFadden et al., 2020). Therefore, the lack of training in the summer often results in a short, high volume pre-season preparation at the NCAA soccer, increasing injuries, and compromising the physical and psychological preparation of the player for the competitive season (Eckard et al., 2018).

Integration of an athlete monitoring program would be useful to develop physical and psychological preparation of the player in NCAA soccer. Athlete monitoring programs primarily aim to maximize sports performance and improve program efficacy while minimize fatigue during the competitive season (Ishida et al., 2021b). Specifically, one aspect of an athlete monitoring program is aimed at quantifying dose-response relationships in order to manage fatigue and to improve training efficacy during the competitive season (Halson, 2014; Gabbett et al., 2017). The monitoring data can inform the practitioners about a current physical and psychological state of the player and assist their decision-making for TL management. Common measures of an athlete monitoring program, particularly for fatigue management, include the Global Navigation Satellite System (GNSS), subjective recovery and stress state, and the countermovement jump (CMJ) (Sams et al., 2018; Travis et al., 2018, 2020; Walker et al., 2019; McFadden et al., 2020; Draper et al., 2021; Ishida et al., 2021c). The GNSS devices have been used to quantify external TL, including the TL of soccer match play, which have been associated with acute muscle damage and the alterations in neuromuscular performance (de Hoyo et al., 2016; Russell et al., 2016; Coppalle et al., 2019; Wiig et al., 2019; Ishida et al., 2021a). Ishida et al. (2021b) reported the strong negative correlations between the changes in loaded CMJ, peak power (PP), and total running distance (p = 0.02, r = 0.65) at 12 h post-match in the NCAA female soccer players. Therefore, the assessment and manipulation of TL using GNSS can be beneficial for maximizing neuromuscular performance.

Measures of subjective recovery and stress state and CMJ are common monitoring tools to quantify the response to TL. The combined use of these measures with GNSS allows for quantifying a dose–response relationship (Halson, 2014; Impellizzeri et al., 2019; Draper et al., 2021; Ishida et al., 2021b). The Short Recovery and Stress Scale (SRSS) is a reliable questionnaire (Kellmann and Kölling, 2019) and consists of eight subscales including physical performance capability (PPC), mental performance capability (MPC), emotional balance (EB), overall recovery (OR), muscular stress (MS), lack of activation (LA), negative emotional state (NES), and overall stress (OS). Current literature (Wiewelhove et al., 2016; Hitzschke et al., 2017; Pelka et al., 2018; Travis et al., 2020) has shown that the SRSS scales are reflective of acute high TLs. The CMJ is also an easy and reliable measure to assess the neuromuscular performance in female soccer players (Andersson et al., 2008; Nedelec et al., 2014; de Hoyo et al., 2016; Claudino et al., 2017). Current evidence also indicates that CMJ performance alternations can be affected by TL and are associated with neuromuscular and physiological damage (Andersson et al., 2008; de Hoyo et al., 2016; Silva et al., 2018; Hader et al., 2019). For example, Andersson et al. (2008) found that statistically significant decreases were observed in CMJ jump height with substantially increased serum creatine kinase (CK) at 21 h after a soccer match play in the elite female soccer players. Therefore, the SRSS and CMJ measure monitoring subjective and objective responses to TL that may provide a better understanding of the response to the manipulations and variations of TL.

Among the soccer players, seasonal variability of TL, physical capacity, performance, and subjective recovery and stress state should be considered. Several studies (Malone et al., 2015; Clemente et al., 2019, 2020; Nobari et al., 2021b) indicate that TL can be substantially higher during a preseason in soccer than the typical loads during the competitive season. For example, Clemente et al. (2020) reported that the weekly total distance covered was considerably higher during pre-season than the end of the competitive season in the professional male soccer players. As TL tends to decrease, physical performance tends to increase as the competitive season progresses (Dragijsky et al., 2017; Sams et al., 2018; Emmonds et al., 2020). However, seasonal variability could be problematic, particularly for quantifying the relationship between TL and the alterations of neuromuscular performance. Although the meta-analyses (Silva et al., 2018; Hader et al., 2019) showed TL is inversely related to physical performance, no investigations have been performed to quantify the relationship interacting with seasonal variability. This information would provide the practitioners with an understanding of how monitoring measures can be useful for associating neuromuscular performance and subjective recovery status with TL variations.

With respect to the NCAA female soccer season, athlete monitoring measures may be incorporated to appropriately guide physical preparation of the player from the pre-season to the end of the competitive season (Ishida et al., 2021b). Although the seasonal variations in TL and neuromuscular performance in the soccer players would occur (Dragijsky et al., 2017; Sams et al., 2018; Walker et al., 2019; Emmonds et al., 2020), few investigations have been performed as to how measures of neuromuscular performance, subjective recovery, and stress state will change across the competitive season. Therefore, the purpose of this study was to investigate as follows: (1) seasonal changes in TL, neuromuscular performance, and subjective recovery and stress state and (2) to examine the relationship between TL and neuromuscular changes in the division I collegiate female soccer players. It was hypothesized that: (1) TL, neuromuscular performance, and subjective recovery and stress state would vary and (2) negative relationships would be observed between TL and neuromuscular performance across the competitive season in the Division I collegiate female soccer players.

Materials and Methods

Design

Data collection occurred throughout the 2019 NCAA soccer season consisting of 2 weeks of pre-season and 11 weeks of the competitive season, respectively. Data were collected by using a GNSS for TL. The average of 2 days (D₂), 7 days (D₇), and 21 days (D₂₁) GNSS TL was calculated prior to four different match play: 1st match of the competitive season (NC₁; 1st non-conference match play), 6th match-play of the competitive season (NC₂; 6th non-conference match play), 9th match of the competitive season (C₁; 1st conference match play), and 15th match of the competitive season (C₂; 6th conference match play). Neuromuscular performance and subjective recovery state were evaluated via the CMJ and SRSS assessments at 3 h prior to NC₁, NC₂, C₁, and C₂. Data collection at NC₁, NC₂, C₁, and C₂ were performed at the first match-play of each week.

Subjects

In this study, the NCAA division I nine female soccer players were included [age 20.0 ± 1.7 years (age range: 18–22 years); body mass 60.3 ± 6.3 kg; height 164.0 ± 5.8 cm; resistance training experience 1–4 years]. Demographic information was collected on the first day of the pre-season. The inclusion criteria for this study were as follows: (a) players were outfield players (defender, midfielder, or forward) and (b) must have completed all testing sessions. This study included six starters (played 74 to 100% of match time across all the NCAA matches) and three non-starters (played 28 to 44% of match time across all the NCAA matches). Six players were excluded from this study because they could not complete CMJ tests. Depending upon the travel schedule, the players completed 0 to 1 strength maintenance weight training sessions each week during pre-season and non-conference play (day 1–51). Weight training was performed two times per week during the conference play (day 52–82). During this study, all the participants were informed of the risks and benefits and testing procedures before participation. The participants signed informed consent and this study was approved by the University Institutional Review Board.

Training Load Measure

D₂, D₇, and D₂₁ TLs were measured by using a GNSS (10 Hz) and accelerometry (triaxial; 100 Hz) units (Catapult OptimEye S5, Catapult Innovations, Team Sport 5.0, Melbourne, Australia). A team of sports scientists powered all the units at least 10 min before the on-field warm-up. Players wore the GNSS units in a vest, which positioned the unit between the shoulder blades. Training and match-derived TL data included the warm-up until the end of the session. Variables of interest related to training volume included total distance (m), high-speed running distance (HSR; m), and PlayerLoad (au). HSR was defined as running above 15 km h⁻¹. PlayerLoad was calculated as the square root of the sum of the squared differences of acceleration in all the three axes divided by the device sampling frequency of 100 Hz (Nicolella et al., 2018). According to the previous literature (Scott et al., 2016; Nicolella et al., 2018; Nikolaidis et al., 2018), a 10-Hz GNSS unit demonstrates good-to-moderate reliability for total distance [coefficient of variation (CV) = 1.9%] and running involving accelerations (CV = 1.9–4.3%) and PlayerLoad (CV = 0.0–3.0% in anterior-posterior, medial-lateral, and vertical axes). The average TL of 2, 7, and 21 days prior to an average match play was calculated (Carey et al., 2017; Sams et al., 2018).

Short Recovery and Stress Scale

The player recovery and stress state were used to assess the subjective recovery and stress state across the 10 week competitive season via the SRSS. The players completed the SRSS from an online-based application (Google Forms, Google, California, United States). All the players were fully familiarized with the procedures using a pre-season match. Prior to the SRSS, hydration status was assessed using a refractometer (ATAGO Corporation Limited, Tokyo, Japan) before the SRSS measurement. The players were considered as hydrated if specific gravity of urine was <1.020. The SRSS is rated by using a seven-point Likert scale from 0 (does not fully apply) to 6 (fully applies) and consists of eight subscales including PPC, MPC, EB, OR, MS, LA, NES, and OS. The Recovery Scale (RS) includes PPC, MPC, EB, and OR, while the Stress Scale (SS) includes MS, LA, NES, and OS. The SRSS has shown acceptable internal reliability (α = 0.74 and α = 0.78) (Kellmann and Kölling, 2019).

Countermovement Jump

The players completed a standardized dynamic warm-up followed by submaximal CMJs at 75% and 100% of perceived maximal efforts, respectively. The players then performed three maximal CMJ trials with a polyvinyl chloride pipe (CMJ0) and with 20 kg barbell (CMJ20) on dual portable force plates (PASPORT Force Platform, PASCO, California, United States of America) by using a sampling frequency of 1,000 Hz. Each load was held across the back on the shoulders. For the CMJ tests, the players stood still on the force plates for at least 1 s and then vertically jumped after flexing the hip, knee, and ankle joints on the command of “3, 2, 1, jump!”. Approximately, a 1-min interval was provided between the CMJ0 and CMJ20 trials. After CMJ testing, the raw data were converted into a comma-separated values file and then analyzed by using a Microsoft Excel sheet (Microsoft Excel, Microsoft, Washington, United States of America) (Chavda et al., 2018). The mean of two trials with the best jump heights (JHs) was used for analysis. Body mass (BM; kg), JH from impulse (cm), modified reactive index (RSI; m s⁻¹), peak force (PF; N), relative PF (RPF; N kg⁻¹), eccentric impulse (EI; N s⁻¹), concentric impulse (CI; N s⁻¹), (PP; W), relative PP (RPP; W kg⁻¹), eccentric peak power (EPP; W), and concentric peak power (CPP; W) were included as variables of interest. The test-retest reliability of the variables was acceptable in CMJ0 [CV = 2.1–5.9%; intraclass correlation coefficient (ICC) = 0.86–0.97] and CMJ20 (CV = 2.7–6.2%; ICC = 0.76–0.92).

Statistical Analysis

All the statistical procedures were performed by using the statistical software RStudio (version 1.1.463) and the packages dplyr (0.8.5), rstatix (0.4.0), and stats (3.5.3). One-way repeated analysis of variance was conducted to examine the difference in the CMJ0 and CMJ20 kinetic variables and the TLs across the four different periods (NC₁, NC₂, C₁, and C₂). For the SRSS and the CMJ variables (JH in CMJ0 and JH and PF in CMJ20) that did not meet the assumption of normality, the Friedman test was performed to identify the differences between the periods. When necessary, post-hoc testing with the Bonferroni correction was performed. A Cohen's d_z effect sizes (d_z) were also calculated by using standardized mean difference and were classified as follows; d_z < 0.2 = trivial, 0.2 ≤ d_z < 0.6 = small, 0.6 ≤ d_z < 1.2 = moderate, 1.2 ≤ d_z < 2.0 = large, and d_z ≥ 2.0 = very large (Hopkins et al., 2009). The Pearson coefficient correlation tests were also conducted to examine the relationship between TL and the changes from NC₁ to NC₂, C₁, and C₂ in the selected CMJ0 and CMJ20 kinetic variables (JH, RSI, PF, and PP). The correlation coefficient magnitudes were determined by using Hopkin's classification (Hopkins et al., 2009) and were classified as follows: r < 0.10 = trivial, 0.10 ≤ r < 0.30 = small, 0.30 ≤ r < 0.50 = moderate, 0.5 ≤ r < 0.7 = large, 0.70 ≤ r < 0.90 = very large, 0.90 ≤ r < 1.00 = nearly perfect, and 1.0 = perfect. All the data were expressed as mean ± SD. Statistical significance was set at p ≤ 0.05.

Results

Seasonal Changes in Training Loads, Countermovement Jump, Short Recovery, and Stress State

Statistically significant differences were observed between NC₁ and C₁ in D₂ total distance (p = 0.02, d_z = 1.53) and PlayerLoad (p = 0.05, d_z = 1.24). D₂ total distance and PlayerLoad were also statistically decreased from NC₂ to C₁ (total distance, p < 0.001, d_z = 2.68; PlayerLoad, p = 0.04, d_z = 1.28) and C₂ (total distance, p = 0.05, d_z = 1.86; PlayerLoad, p = 0.03, d_z = 1.35). Additionally, statistically significant differences were found from NC₁ to C₁ in D₇ total distance (p = 0.03, d_z = 1.40). D₇ total distance and PlayerLoad were statistically decreased from NC₂ to C₁ (total distance, p = 0.03, ES = 1.40; PlayerLoad, p = 0.02, d_z = 1.72) and C₂ (total distance, p = 0.03, d_z = 1.40; PlayerLoad, p = 0.02, d_z = 1.49) (Table 1). However, no statistical differences were noted for D₂₁ total distance, HSR, and PlayerLoad (p > 0.05).

TABLE 1

Table 1. Seasonal changes in training loads.

Statistically significant differences were observed from NC₁ to C₁ in CMJ0 JH (p = 0.03, d_z = 1.40), RSI (p = 0.02, d_z = 1.44), PP (p = 0.034, d_z = 1.32), and CAP (p = 0.01, d_z = 1.74) (Figure 1). Statistically significant differences were also noted from NC₁ to C₂ for CMJ0 JH (p = 0.015, d_z = 1.53), PP (p = 0.01, d_z =1.57), and RPP (p = 0.03, d_z =1.40). Additionally, CMJ20 JH, PP, and RPP showed statistically significant differences from NC₁ to C₂ (JH, p = 0.019, d_z = 1.47; PP, p = 0.03, d_z = 1.37; RPP, p = 0.02, d_z = 1.49). However, no statistically significant changes were seen in BM and CMJ0 and CMJ20 PF, RPF, and EAP (p > 0.05) (Table 2). No significant changes were observed in any of the SRSS items across time (p > 0.05).

FIGURE 1

Figure 1. Seasonal changes in unloaded countermovement jump height (A), modified strength index (RSI) (B), and peak power (C). NC₁ = 1st non-conference match-play. NC₂ = 7th non-conference match-play. C₁ = 1st conference match-play. C₂ = 6th conference match-play. *Denotes p ≤ 0.05 from NC₁. ^§Denotes p ≤ 0.05 from NC. ^¶Denotes p ≤ 0.05 from C₁.

TABLE 2

Table 2. Seasonal changes in unloaded and loaded countermovement jump kinetics.

Seasonal Relationships Between Training Loads and the Changes in Countermovement Jump Performance

Very large positive correlations were observed at C₁ between D₂ total PlayerLoad and the changes from NC₁ to C₁ in CMJ0 JH (p = 0.05, r = 0.67); D₇ total PlayerLoad and CMJ0 JH (p = 0.02, r = 0.80), RSI (p = 0.03, r = 0.80), and PF (p = 0.02, r = 0.77). In CMJ20, very large positive correlations were also found between D₂ total PlayerLoad and the changes from NC₁ to C₁ in CMJ0 JH (p < 0.001, r = 0.88), RSI (p = 0.001, r = 0.86), and PF (p = 0.001, r = 0.93); D₇ total PlayerLoad and JH (p = 0.01, r = 0.87) and PF (p < 0.001, r = 0.90). At C₂, very large negative correlations were observed between D₇ total distance and the changes from NC₁ to C₂ CMJ0 JH (p = 0.02, r = 0.79) and PP (p = 0.03, r = 0.71); D₂₁ total distance and CMJ0 JH (p = 0.04, r = 0.70) and PP (p = 0.04, r = 0.69) (Figure 2).

FIGURE 2

Figure 2. Correlations between percentage changes from 1st non-conference match-play to 6th conference match-play inunloaded countermovement jump (CMJO) jump height and peak power and total distance. (A) Between jump height and D₇ total distance, (B) between jump height and D₂₁ total distance, (C) peak power and D₇ total distance, and (D) peak power and D₂₁ total distance. D₇ = average of 7 days. D₂₁ = average of 21 days.

Discussion

The purposes of this study were to investigate: (1) seasonal changes in TL, neuromuscular performance, subjective recovery, and stress state and (2) the relationships between TL and neuromuscular changes in division I collegiate female soccer players. The main findings of this study were: (1) D₇ total distance and PlayerLoad were lower at the conference period than the pre-season and non-conference periods, (2) CMJ0 JH, PP, and RSI were statistically higher at C₁ than NC₁, (3) positive correlations were found between D₇ PlayerLoad and unloaded CMJ JH, PP, and PF at C₁ while negative correlations were observed between D₂ and D₇ total distance and CMJ0 JH and PP at C₂. This study showed that neuromuscular performance gradually increased from the pre-season to the conference play and there was a negative relationship between D₇ TL and CMJ at C₂.

In professional and NCAA soccer, a pre-season is considered as a physical preparatory phase to develop physical capabilities for 2–4 weeks prior to the competitive season. However, practitioners are aware that the length of the pre-season would not be sufficient to adequately improve characteristics of fitness (Dragijsky et al., 2017; Sams et al., 2018; Emmonds et al., 2020; McFadden et al., 2020). The short pre-season may also not provide sufficient recovery time between training sessions, potentially resulting in greater physiological and psychological damage (McFadden et al., 2020; Nobari et al., 2021a,b). Our data indicate that D₇ total distance and PlayerLoad at C₁ were statistically lower than NC₁ and NC₂. This agrees with current evidence that soccer players will typically accumulate higher TLs during the pre-season and the early competitive season (Malone et al., 2015; Clemente et al., 2019, 2020). For example, Malone et al. (2015) found that total distance was substantially higher during the early competitive season (weeks 7–12) compared to the end of the competitive season (weeks 37–42) in elite English male soccer players (p < 0.05, ES = 0.84). Therefore, similar to previous investigations (Malone et al., 2015; Clemente et al., 2019, 2020), TL was highest during the pre-season and the early competitive season in the Division I NCAA female soccer players. However, the NCAA female soccer players may not accumulate sufficient TL after the summer break (May–July). Therefore, practitioners may need to manipulate and progress from the pre to the early competitive season to minimize the risk of injuries.

CMJ0 JH, RSI, and PP statistically increased from NC₁ to C₁, although weight training frequency was inconsistent and limited from NC₁ to C₁. Similar to our findings, several studies (Dragijsky et al., 2017; Sams et al., 2018; Emmonds et al., 2020) have shown that neuromuscular performance improves among soccer players from pre-season to the mid or end of the competitive season. Sams et al. (2018) reported that squat JH showed a moderate increase from baseline to the 8th match of the competitive season (p = 0.039, ES = 1.01) in the NCAA Division I female soccer players. Improved neuromuscular performance at the mid and end of the competitive season may be explained by the training status of the players prior to the NCAA pre-season. The NCAA prohibits division I soccer teams from starting pre-season until 2–3 weeks prior to the first competition (National Collegiate Athletic National Collegiate Athlete Association, 2021). Nonetheless, the NCAA also restricts strength and conditioning coaches from having mandatory physical training sessions and monitoring their TL during the summer period (12–14 weeks from May to early August). The NCAA restrictions could be detrimental for reasonable TL maintenance or accumulation during the summer and can result in a sudden increase in TL at pre-season, increasing the risk of injuries, muscle damage, and autonomic nervous system fatigue (Agel and Schisel, 2013; Walker et al., 2019; Sekiguchi et al., 2021). The 2–3 weeks of pre-season, as a result of the NCAA regulation, may leave players physically underprepared during the non-conference play and the 12–16 weeks of the competitive season.

The correlations between D₇ TL and CMJ0 changes were positive at C₁, while very large negative correlations were observed between D₇ total distance and CMJ0 JH and PP at C₂. These findings may indicate that the relationship between TL and neuromuscular performance may be altered across the competitive season. Our finding at C₂ agrees with the numerous previous findings (Andersson et al., 2008; Silva et al., 2018; Hader et al., 2019) indicating that jump performance can be negatively affected both acutely and chronically by training. For example, a meta-analysis by Hader et al. (2019) reported that a large negative correlation was found between high-speed running distance (>19.8 km/h) and the CMJ0 PP (r = 0.52, 95% CI 0.64–0.40) at 24 h post-match. However, the correlations between D₇ TL and CMJ0 changes were positive at C₁. The disagreement between our finding and previous literature (Andersson et al., 2008; Silva et al., 2018; Hader et al., 2019) may be explained by the training status of the player. The NCAA soccer players may be in physically undertrained status due to the lack of insufficient training volume prior to the pre-season. When prescribing intense/high volume during 2–3 weeks of a pre-season, the training prescription may increase the neuromuscular performance of the player from the pre-season to the mid and late competitive season (Dragijsky et al., 2017; Sams et al., 2018; Emmonds et al., 2020) resulting in the positive correlations between D₇ TL and CMJ0 changes from NC₁ to C₁. Based on our findings and previous literature (Dragijsky et al., 2017; Sams et al., 2018; Emmonds et al., 2020), the effects of seasonal variations may affect the assessment of neuromuscular performance in relation to TL in NCAA soccer.

Limitations

There are three main limitations of this study. First, the sample size of this study was limited. The data were collected as a part of ongoing athlete monitoring, so this study could not maintain large sample sizes due to injuries affecting jump testing. Second, there were no measures of additional physical performance abilities such as maximum strength, sprinting, change of directions, and intermittent endurance performance over the period. Third, no internal TL measures were included in this study. Future investigation should include other performance tests and examine the effects of weight training sessions with larger sample sizes.

Practical Application

In the NCAA Division I female soccer players, CMJ0 and CMJ20 may increase from the first match play to the midcompetitive and late season. Practitioners (e.g., sports coaches, strength and conditioning coaches, sports scientists) should be aware that longer-term strategic training plans may be required to develop and maximize neuromuscular performance at the end of the competitive season. The practitioners should also carefully impart to collegiate athletes the importance of the quantification of the summer and pre-season TL for maximizing neuromuscular performance at the early competitive season. In addition, care should also be taken when analyzing and interpreting the relationship between acute TL and CMJ0 performance due to the seasonal variations associated with physical preparedness of the player. Our findings demonstrate that CMJs may be a worthwhile test to quantify neuromuscular alternations at the mid and end of the competitive season.

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 studies involving human participants were reviewed and approved by East Tennesse State University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

AI, CB, AS, MS, and JG contributed to study design and implementation. AI and JG carried out all the data collection and analysis. All authors contributed to data analysis, interpretation, discussion of the results, editing and reviewing of the article, and read and agreed to the final version of the submitted manuscript.

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.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspor.2021.757253/full#supplementary-material

References

Agel, J., and Schisel, J. (2013). Practice injury rates in collegiate sports. Clin. J. Sport. Med. 23, 33–38. doi: 10.1097/JSM.0b013e3182717983

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, H. A., Raastad, T., Nilsson, J. E., Paulsen, G., Garthe, I., and Kadi, F. (2008). Neuromuscular fatigue and recovery in elite female soccer: effects of active recovery. Med. Sci. Sports Exerc. 40, 372–380. doi: 10.1249/mss.0b013e31815b8497

PubMed Abstract | CrossRef Full Text | Google Scholar

Carey, D. L., Blanch, P., Ong, K.-L., Crossley, K. M., Crow, J., and Morris, M. E. (2017). Training loads and injury risk in australian football-differing acute: chronic workload ratios influence match injury risk. Br. J. Sports Med. 51, 1215–1220. doi: 10.1136/bjsports-2016-096309

PubMed Abstract | CrossRef Full Text | Google Scholar

Chavda, S., Bromley, T., Jarvis, P., Williams, S., Bishop, C., Turner, A. N., et al. (2018). Force-time characteristics of the countermovement jump: analyzing the curve in excel. Strength Cond. J. 40, 67–77. doi: 10.1519/SSC.0000000000000353

CrossRef Full Text | Google Scholar

Claudino, J. G., Cronin, J., Mezêncio, B., McMaster, D. T., McGuigan, M., Tricoli, V., et al. (2017). The countermovement jump to monitor neuromuscular status: a meta-analysis. J. Sci. Med. Sport. 20, 397–402. doi: 10.1016/j.jsams.2016.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Clemente, F. M., Owen, A., Serra-Olivares, J., Nikolaidis, P. T., van der Linden, C. M. I., and Mendes, B. (2019). Characterization of the weekly external load profile of professional soccer teams from portugal and the netherlands. J. Hum. Kinet. 66, 155–164. doi: 10.2478/hukin-2018-0054

PubMed Abstract | CrossRef Full Text | Google Scholar

Clemente, F. M., Silva, R., Castillo, D., Los Arcos, A., Mendes, B., and Afonso, J. (2020). Weekly load variations of distance-based variables in professional soccer players: a full-season study. Int. J. Environ. Res. Public. Health. 17:3300. doi: 10.3390/ijerph17093300

PubMed Abstract | CrossRef Full Text | Google Scholar

Coppalle, S., Rave, G., Ben Abderrahman, A., Ali, A., Salhi, I., Zouita, S., et al. (2019). Relationship of pre-season training load with in-season biochemical markers, injuries and performance in professional soccer players. Front. Physiol. 10:e00409. doi: 10.3389/fphys.2019.00409

PubMed Abstract | CrossRef Full Text | Google Scholar

de Hoyo, M., Cohen, D. D., Sañudo, B., Carrasco, L., Álvarez-Mesa, A., Del Ojo, J. J., et al. (2016). Influence of football match time-motion parameters on recovery time course of muscle damage and jump ability. J. Sports. Sci. 34, 1363–1370. doi: 10.1080/02640414.2016.1150603

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Draper, G., Wright, M., Chesterton, P., and Atkinson, G. (2021). The tracking of internal and external training loads with next-day player-reported fatigue at different times of the season in elite soccer players. Int. J. Sports Sci. Coach. 16, 793–803. doi: 10.1177/1747954121988960

CrossRef Full Text | Google Scholar

Eckard, T. G., Padua, D. A., Hearn, D. W., Pexa, B. S., and Frank, B. S. (2018). The relationship between training load and injury in athletes: a systematic review. Sports Med. 48, 1929–1961. doi: 10.1007/s40279-018-0951-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Emmonds, S., Sawczuk, T., Scantlebury, S., Till, K., and Jones, B. (2020). Seasonal changes in the physical performance of elite youth female soccer players. J. Strength Cond. Res. 34, 2636–2643. doi: 10.1519/JSC.0000000000002943

PubMed Abstract | CrossRef Full Text | Google Scholar

Gabbett, T. J., Nassis, G. P., Oetter, E., Pretorius, J., Johnston, N., Medina, D., et al. (2017). The athlete monitoring cycle: a practical guide to interpreting and applying training monitoring data. Br. J. Sports Med. 51, 1451–1452. doi: 10.1136/bjsports-2016-097298

PubMed Abstract | CrossRef Full Text | Google Scholar

Hader, K., Rumpf, M. C., Hertzog, M., Kilduff, L. P., Girard, O., and Silva, J. R. (2019). Monitoring the athlete match response: can external load variables predict post-match acute and residual fatigue in soccer? A systematic review with meta-analysis. Sports Med. Open 5, 1–19. doi: 10.1186/s40798-019-0219-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Halson, S. L. (2014). Monitoring training load to understand fatigue in athletes. Sports Med. 44, 139–147. doi: 10.1007/s40279-014-0253-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Hitzschke, B., Wiewelhove, T., Raeder, C., Ferrauti, A., Meyer, T., Pfeiffer, M., et al. (2017). Evaluation of psychological measures for the assessment of recovery and stress during a shock-microcycle in strength and high-intensity interval training. Perform. Enhanc. Health 5, 147–157. doi: 10.1016/j.peh.2017.08.001

CrossRef Full Text | Google Scholar

Hopkins, W. G., Marshall, S. W., Batterham, A. M., and Hanin, J. (2009). Progressive statistics for studies in sports medicine and exercise science. Med. Sci. Sports Exerc. 41, 3–13. doi: 10.1249/MSS.0b013e31818cb278

PubMed Abstract | CrossRef Full Text | Google Scholar

Impellizzeri, F. M., Marcora, S. M., and Coutts, A. J. (2019). Internal and external training load: 15 years on. Int. J. Sports Physiol. Perform. 14, 270–273. doi: 10.1123/ijspp.2018-0935

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishida, A., Bazyler, C. D., Sayers, A. L., Mizuguchi, S., and Gentles, J. A. (2021a). Acute effects of match-play on neuromuscular and subjective recovery and stress state in division I collegiate female soccer players. J. Strength Cond. Res. 35, 976–982. doi: 10.1519/JSC.0000000000003981

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishida, A., Bazyler, C. D., Sayers, A. L., Stone, M. H., and Gentles, J. A. (2021b). Evidence and application of athlete monitoring programs in national collegiate athletic association women's soccer: a narrative review. Strength Cond. J. ahead of print, 1–13. doi: 10.1519/SSC.0000000000000670

CrossRef Full Text | Google Scholar

Ishida, A., Travis, S. K., and Stone, M. H. (2021c). Associations of body composition, maximum strength, power characteristics with sprinting, Jumping, and intermittent endurance performance in male intercollegiate soccer players. J. Funct. Morphol. Kinesiol. 6:7. doi: 10.3390/jfmk6010007

PubMed Abstract | CrossRef Full Text | Google Scholar

Kellmann, M., and Kölling, S. (2019). Recovery and Stress in Sport: A Manual for Testing and Assessment. Abingdon: Routledge.

Google Scholar

Malone, J. J., Di Michele, R., Morgans, R., Burgess, D., Morton, J. P., and Drust, B. (2015). Seasonal training-load quantification in elite english premier league soccer players. Int. J. Sports Physiol. Perform. 10, 489–497. doi: 10.1123/ijspp.2014-0352

PubMed Abstract | CrossRef Full Text | Google Scholar

McFadden, B. A., Walker, A. J., Bozzini, B. N., Sanders, D. J., and Arent, S. M. (2020). Comparison of internal and external training loads in male and female collegiate soccer players during practices vs. games. J. Strength Cond. Res. 34, 969–974. doi: 10.1519/JSC.0000000000003485

PubMed Abstract | CrossRef Full Text | Google Scholar

National Collegiate Athlete Association (2021). Rules and Bylaws: NCAA Division I Manual. Indianapolis, IN. Available online at: https://web3.ncaa.org/lsdbi/reports/getReport/90008 (accessed September 30, 2021).

Nedelec, M., McCall, A., Carling, C., Legall, F., Berthoin, S., and Dupont, G. (2014). The influence of soccer playing actions on the recovery kinetics after a soccer match. J. Strength Cond. Res. 28, 1517–1523. doi: 10.1519/JSC.0000000000000293

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicolella, D. P., Torres-Ronda, L., Saylor, K. J., and Schelling, X. (2018). Validity and reliability of an accelerometer-based player tracking device. PLoS ONE 13:e0191823. doi: 10.1371/journal.pone.0191823

PubMed Abstract | CrossRef Full Text | Google Scholar

Nikolaidis, P. T., Clemente, F. M., van der Linden, C. M. I., Rosemann, T., and Knechtle, B. (2018). Validity and reliability of 10-Hz global positioning system to assess in-line movement and change of direction. Front. Physiol. 9:e00228. doi: 10.3389/fphys.2018.00228

PubMed Abstract | CrossRef Full Text | Google Scholar

Nobari, H., Alves, A. R., Haghighi, H., Clemente, F. M., Carlos-Vivas, J., Pérez-Gómez, J., et al. (2021a). Association between training load and well-being measures in young soccer players during a season. Int. J. Environ. Res. Public Health. 18:4451. doi: 10.3390/ijerph18094451

PubMed Abstract | CrossRef Full Text | Google Scholar

Nobari, H., Mainer-Pardos, E., Adsuar, J. C., Franco-García, J. M., Rojo-Ramos, J., Cossio-Bolaños, M. A., et al. (2021b). Association between endocrine markers, accumulated Workload, and fitness parameters during a season in elite young soccer players. Front. Psychol. 12:3592. doi: 10.3389/fpsyg.2021.702454

PubMed Abstract | CrossRef Full Text | Google Scholar

Pelka, M., Schneider, P., and Kellmann, M. (2018). Development of pre- and post-match morning recovery-stress states during in-season weeks in elite youth football. Sci. Med. Football. 2, 127–132. doi: 10.1080/24733938.2017.1384560

CrossRef Full Text | Google Scholar

Russell, M., Sparkes, W., Northeast, J., Cook, C. J., Bracken, R. M., and Kilduff, L. P. (2016). Relationships between match activities and peak power output and creatine kinase responses to professional reserve team soccer match-play. Hum. Mov. Sci. 45, 96–101. doi: 10.1016/j.humov.2015.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Sams, M. L., Sato, K., DeWeese, B. H., Sayers, A. L., and Stone, M. H. (2018). Quantifying changes in squat jump height across a season of men's collegiate soccer. J. Strength Cond. Res. 32, 2324–2330. doi: 10.1519/JSC.0000000000002118

PubMed Abstract | CrossRef Full Text | Google Scholar

Sams, M. L., Wagle, J. P., Sato, K., DeWeese, B. H., Sayers, A. L., and Stone, M. H. (2020). Using the session rating of perceived exertion to quantify training load in a men's college soccer team. J. Strength Cond. Res. 34, 2793–2799. doi: 10.1519/JSC.0000000000003793

PubMed Abstract | CrossRef Full Text | Google Scholar

Scott, M. T. U., Scott, T. J., and Kelly, V. G. (2016). The validity and reliability of global positioning systems in team sport: a brief review. J. Strength Cond. Res. 30, 1470–1490. doi: 10.1519/JSC.0000000000001221

PubMed Abstract | CrossRef Full Text | Google Scholar

Sekiguchi, Y., Huggins, R. A., Curtis, R. M., Benjamin, C. L., Adams, W. M., Looney, D. P., et al. (2021). Relationship between heart rate variability and acute: chronic load ratio throughout a season in ncaa d1men's soccer players. J. Strength Cond. Res. 35, 1103–1109. doi: 10.1519/JSC.0000000000002853

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva, J. R., Rumpf, M. C., Hertzog, M., Castagna, C., Farooq, A., Girard, O., et al. (2018). Acute and residual soccer match-related fatigue: a systematic review and meta-analysis. Sports Med. 48, 539–583. doi: 10.1007/s40279-017-0798-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Stølen, T., Chamari, K., Castagna, C., and Wisløff, U. (2005). Physiology of soccer: an update. Sports Med. 35, 501–536. doi: 10.2165/00007256-200535060-00004

PubMed Abstract | CrossRef Full Text | Google Scholar

Travis, S. K., Goodin, J. R., Beckham, G. K., and Bazyler, C. D. (2018). Identifying a test to monitor weightlifting performance in competitive male and female weightlifters. Sports (Basel) 6, 1–12. doi: 10.3390/sports6020046

PubMed Abstract | CrossRef Full Text | Google Scholar

Travis, S. K., Mizuguchi, S., Stone, M. H., Sands, W. A., and Bazyler, C. D. (2020). Preparing for a national weightlifting championship: a case series. J. Strength Cond. Res. 34, 1842–1850. doi: 10.1519/JSC.0000000000003312

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, A. J., McFadden, B. A., Sanders, D. J., Rabideau, M. M., Hofacker, M. L., and Arent, S. M. (2019). Biomarker response to a competitive season in division I female soccer players. J. Strength Cond. Res. 33, 2622–2628. doi: 10.1519/JSC.0000000000003264

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiewelhove, T., Raeder, C., Meyer, T., Kellmann, M., Pfeiffer, M., and Ferrauti, A. (2016). Effect of repeated active recovery during a high-intensity interval-training shock microcycle on markers of fatigue. Int. J. Sports. Physiol. Perform. 11, 1060–1066. doi: 10.1123/ijspp.2015-0494

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiig, H., Raastad, T., Luteberget, L. S., Ims, I., and Spencer, M. (2019). External load variables affect recovery markers up to 72 h after semiprofessional football matches. Front. Physiol. 10:689. doi: 10.3389/fphys.2019.00689

PubMed Abstract | CrossRef Full Text | Google Scholar