Skip to main content

Effects of two neuromuscular training programs on running biomechanics with load carriage: a study protocol for a randomised controlled trial



In recent years, athletes have ventured into ultra-endurance and adventure racing events, which tests their ability to race, navigate, and survive. These events often require race participants to carry some form of load, to bear equipment for navigation and survival purposes. Previous studies have reported specific alterations in biomechanics when running with load which potentially influence running performance and injury risk. We hypothesize that a biomechanically informed neuromuscular training program would optimize running mechanics during load carriage to a greater extent than a generic strength training program.


This will be a two group, parallel randomized controlled trial design, with single assessor blinding. Thirty healthy runners will be recruited to participate in a six weeks neuromuscular training program. Participants will be randomized into either a generic training group, or a biomechanically informed training group. Primary outcomes include self-determined running velocity with a 20 % body weight load, jump power, hopping leg stiffness, knee extensor and triceps-surae strength. Secondary outcomes include running kinetics and kinematics. Assessments will occur at baseline and post-training.


To our knowledge, no training programs are available that specifically targets a runner’s ability to carry load while running. This will provide sport scientists and coaches with a foundation to base their exercise prescription on.

Trial registration

ANZCTR (ACTRN12616000023459) (14 Jan 2016)

Peer Review reports


Ever since the “Battle of Marathon” between Greece and Persia was recorded [1], the accomplishment of running a 42 km marathon is seen as the ultimate achievement for a distance runner. However, the last two decades has seen both recreational and elite level runners striving for distances well beyond a standard marathon. Interest and participation in ultra-endurance races [2], multi-stage racing events [3], and off-terrain trail and adventure races have risen [4] as individuals seek new ways to challenge human limits. These races not only test a runner’s speed and endurance, but also their ability to navigate and survive over undulating terrains and harsh environments [3]. Navigation and survival requires routine access to specialized equipment and sustenance. This requirement necessitates athletes to compete with externally carried loads [5]. Few studies have considered the role of load carriage on the potential impact on a runner’s health and performance [6]. In addition, no studies have considered if runners can be trained to adapt to external loads in running.

Load carriage in running poses two fundamental problems to athletes and occupational personnel: 1) an increased injury rate, 2) and increased metabolic energy expenditure that may reduce performance [7]. In the adult population, the increased overuse injury rate associated with load carriage has largely been investigated in the military setting, where load carriage biomechanics have been predominantly investigated while walking [7]. A previous study reported that 8 % of the 5000 injuries reported in the Australian Defence Force from January 2009 to December 2010, were related to heavy load carriage [8]. Of these injuries, 56 % affected the lower limb and were classified as muscular stress related [8]. Although no causative studies have been performed, it is likely that load carriage while running may exacerbate the already high incidence of running related injuries [9]. In addition, when an individual runs with load, the energy demand involved in maintaining constant running speed is increased [10]. Minimising the reduction in running speed associated with load carriage is important for the survivability of military personnel, the performance of athletes, and the overall efficiency of movement in recreational runners [11, 12].

The risk of injury and reduced performance associated with load carriage in running, points to the need for a preconditioning program for these athletes. There is convincing evidence that resistance based neuromuscular training programs are effective at reducing running related injuries (RRI) during body weight (BW) running (i.e. running with no external load) [13], and improving BW running performance [14]. However, current training programs have been developed using BW running research [15, 16], rather than loaded running research. The only studies that have attempted to define best training practices for load carriage gait has been performed in the military setting [17]. A limitation in existing training studies has been that exercise prescription has not been explicitly informed from biomechanical studies of load carriage gait. Rather, training was of a generalised nature, targeting the large muscle groups of the lower limb [17]. The type of exercises and mode of contractions used for preconditioning programs should be specific to the gait pattern required for athletes, and be based on prior knowledge of biomechanical adaptations during load carriage.

Potential adaptive and mal-adaptive biomechanics

Studies using computed muscle control and induced acceleration analysis have identified the integrated roles of lower limb muscles in BW running. Collectively, the functions of these muscles are to provide a vertical force to accelerate/decelerate body weight, and horizontal forces to accelerate/decelerate inertial mass [18, 19]. When additional load is imposed on a runner, greater vertical and horizontal forces are needed to accelerate and decelerate an increased total weight and total inertial mass, respectively. Biomechanical changes to running with load are classified as adaptive if they enable an increase in baseline motor function (Table 1). For example, an increase in ankle power absorption in mid stance with load may be adaptive as it transfers power away from proximal segments to the foot [6]. This may aid in increased elastic-energy recovery at the ankle plantar flexor muscle-tendon unit, which may be essential to sustain faster running velocities during load carriage.

Table 1 Biomechanical adaptations of load carriage to potentially optimize metabolic cost and minimise injury risk

On the contrary not all biomechanical changes with load may be adaptive. Some mechanical changes are likely to be mal-adaptive as they may contribute to a greater risk of incurring RRI or represent an inefficient running style. Poor hip control of non-sagittal plane rotations has been documented to increase the risk of developing RRIs [20]. At the kinematic level, load carriage has been associated with increased hip adduction at terminal stance [6] (Table 1). This increase in non-sagittal plane movements may represent suboptimal muscle capacity and motor control [21]. In addition, poor proximal trunk-pelvis control in running may result in energy being wasted in maintaining postural balance and inter-segmental alignment. Remediating mal-adaptive mechanical changes whilst enhancing adaptive changes could improve biomechanical indices of running performance and injury risk during load carriage.


Load carriage in running is increasingly common in running related sports. The ability to positively and predictably adapt to the imposed load when running necessitates an evidence-based training program. Existing training studies for load carriage performance in the military setting cannot be immediately applied to load carriage running, as most studies investigated performance in walking. This is because running and walking involve different movement dynamics, making the extrapolation of results from walking studies problematic when applied to running. For example, the hip contributes approximately 20 % of total positive power in the stance phase of BW walking, but less than 10 % of total positive power in the same phase of BW running [22]. Second, studies that have investigated ways to improve load carriage performance have adopted a non-randomized design [17]. Reported effect sizes of benefit in intervention studies were larger in trials without a randomized design compared to one with a randomized design [23]. Lastly, studies on load carriage do not appear to specifically target the known neuromuscular demands of load carriage gait patterns [17]. Therefore, the purpose of this investigation is to compare the effects of a biomechanically informed neuromuscular training program to a generic standard best-practice resistance training program on changes in the biomechanics of running with load.


To compare changes in (1) self-determined running velocity with and without load carriage, (2) lower limb running kinematics and kinetics, (3) jumping power and hopping stiffness, (4) and isokinetic knee and ankle extensor strength in healthy adult runners participating in a biomechanically informed training program compared to a generic resistance training program. This generic resistance training program may be seen as the current “gold-standard” program based on current best evidence [17].


Research design

The study is a single blinded, parallel-grouped randomized controlled trial which will be designed and reported according to the Consolidated Standards of Reporting Trials (CONSORT) statement (Fig. 1) [24]. This study is registered with the Australian New Zealand Clinical Trials Registry (ACTRN: ACTRN12616000023459). The Curtin University Human Research Ethics Committee (RD-41-14) has approved this study protocol. All participants will provide written informed consent prior to study inclusion.

Fig. 1
figure 1

CONSORT Flow Diagram

Participants, setting and recruitment

Runners with a variety of training experience residing in Western Australia will be invited to participate. All assessments and intervention will be conducted within Curtin University, Perth, Australia. Participants between 18 and 60 years old who are in good general health, and have been running or participating in running-related sports with a minimum cumulated total of 4 km/week or 45 min/week over the past 12 months, will be recruited. Exclusion criteria include: the presence of any disorders that could affect their gait and load carrying ability; medical conditions that preclude heavy resistance training and strenuous running; presence of a training-loss running related injury within the last three months [25]; current running related pain (except blisters or muscle soreness) [25]; lower limb surgery within the past 12 months; and females who are pregnant.

Sample size calculation

This study was powered on the effects of a core stability program on changes to hopping leg stiffness [26]. Sample size was planned based on a two way, repeated measures ANOVA, using the Hotelling-Lawley Trace to test for an intervention by time interaction [27]. Previous studies on leg stiffness reported a standard deviation of 3600 N/m [26], and a correlation between repeated measures of 0.80 [28]. For a desired power of 0.80, and a Type 1 error rate of 0.05, 24 participants are needed to detect a between group mean difference of 3000 N/m. In order to account for a 20 % dropout over the six week intervention period, 30 participants will be tested.

Randomization, allocation and blinding

Prior to randomization, participants will be stratified into two groups based on their gender. Previous studies have identified gender differences in BW running mechanics and different associative relationships between running mechanics and economy [16, 29]. Permuted block randomization will be performed within each stratum, using two different block sizes (two blocks of four and three blocks of eight) [30]. Each block consists of either four or eight group assignments (half of the assignments to one of two groups), to ensure that at the end of each block, the number of participants allocated to either intervention group will be balanced. The sequence within each block, and the order of all five blocks per stratum are randomized. The randomization sequence will be generated with an online random sequence generator used in previous a study [31]. When a participant has provided signed-informed consent, an external allocator not involved in the experiment, will sequentially draw an envelope (lowest numbered to highest) from either of two containers, depending on the stratum. The allocator will write the participant’s name, identifier number, date, and allocator’s signature on the envelope, which will be ink-printed onto the treatment allocation card via carbon paper. The envelope’s seal will be broken and the participant and the trainer will be informed about the allocated intervention group. Participant and trainer blinding will not be feasible due to the nature of prescribed intervention. Outcome assessor (for biomechanical and strength assessments) will be blinded to the allocation of participants to intervention groups.

Subjective assessment

The following data will be collected at baseline: 1) participant’s demographics, 2) self-reported running training, load carriage, and strength training history; 3) self-reported medical status (including the Physical Activity Readiness Questionnaire); 4) self-reported running overuse injury history [25, 32].

Three dimensional motion capture - load carriage running protocol

Participants will be carrying a backpack (CAMELBAK, H.A.W.G.® NV,14 l) fitted with chest strap and hip belt secured snuggly. Participants will be wearing their personal running attire and running shoes for all assessments. Participants will run over ground, in a straight line at two velocities: 1) self-determined velocity, and 2) 3.5 m/s, over two external load conditions (0 and 20 % BW). A lead up distance of at least 20 m to the edge of the first force plate and tail off distance of 10 m after the edge of the last force plate will be given to enable sufficient distance for acceleration and deceleration. A velocity of 3.5 m/s was predetermined as it closely represented the running velocity of the fastest men in 24 h ultra-marathons [33]. Sand bags will be filled to 20 % BW and secured within the backpack. A load of 20 % BW will be used as this represents a common relative load borne by tactical athletes during periods of running [11]. Indirect evidence from recommended backpack volume during ultra-endurance and adventure racers, suggest that a 20 % BW load is a reasonable approximate to actual load magnitude carried [5]. The order of load-velocity testing will be completely randomized using an online random sequence generator [6]. Timing gates (SMARTSPEED Pro, Fusion Sport Pty Ltd, Australia) will be positioned on both sides of the force plates (AMTI, Watertown, MA) (2000 Hz) 5 m apart, to monitor running velocity. Familiarisation trials will be given to practice the required running velocity. Participants will be required to complete a minimum of five successful running trials for each condition. Each trial would be interspersed with a 30 s rest break. A successful trial is defined when the right limb contacts the middle of the force platform without alteration of running pattern, within ±10 % of the prescribed speed. A minimum of five minutes rest will be given after each condition.

Three dimensional motion capture – jumping and hopping protocol

Squat jumps (SJ), countermovement jumps (CMJ), and single leg vertical hoping (SL hop) will be performed on the force plates. For both the SJ and CMJ, the assessor will demonstrate and instruct technique, and participants will have to practise these jump techniques until performance is stable. For all tests, participants will be required to fix their arms at 90° abduction (“T” pose), to limit the influence of upper body on jump performance [34]. For the SJ, participants will perform a maximal concentric vertical jump from an initial squat depth of approximately 90° knee angle (visual estimation). For the CMJ, participants will descend into a squat depth of approximately 90° knee angle (visual estimation), followed without pause by a rapid maximal vertical jump and landing in a comfortable squat position [35]. When the assessor is satisfied with the practice performances, three maximal SJ trials and three maximal CMJ trials will be performed and recorded. Both SJ and CMJ will be performed with and without an external 20 % BW load. Each trial will be interspersed with a 30s rest and each test interspersed with a one minute rest to avoid fatigue. SL hopping will be performed for both legs separately over four conditions: 1) self-paced hopping frequency (BW), 2) self-paced hopping frequency with 20 % BW carriage, 3) hopping at a frequency of 3.0 Hz (BW), and 4) hopping at a frequency of 3.0 Hz with 20 % BW. Hopping frequency will be set using a handheld digital metronome. Each hop condition will last 10 s, with a one minute break interspersed [36]. Trials will be repeated if hopping frequency is not maintained. The order of testing will be randomized using an online random sequence generator (

Data measurement and processing

An 18 camera motion capture system (Vicon T-series, Oxford Metrics, UK) (250Hz), synchronised to three consecutive in-ground force plates (three metre long in direction of progression) will be used to collect motion and force data for running, hopping and maximal jump tasks. Data will be captured and stored using manufacturer supplied software (Vicon Nexus, v2.3, Oxford Metrics, UK). Data processing will be performed in Vicon Nexus and Visual 3D (C-motion, Germantown, MD). Marker trajectories and ground reaction forces will be filtered at identical cut-off values, for use in inverse-dynamics calculations [37]. For joint angles, raw marker trajectories will be filtered at a separate cut-off frequency. The choice of cut-off frequencies for motion and force data will be based on past research [37]. A Cardan XYZ rotation sequence will be used to calculate 3D joint angles [38], and both kinematics and kinetics will be expressed in an orthogonal frame in the proximal segment [39]. Data will be computed during the stance and swing period of running. A threshold of 20 N in ground reaction force will used to determine initial contact and toe-off. Kinetic variables will be normalized using base factors of gravitational constant g (9.81 m/s2), leg length, L (m), and body mass, M (kg) [40].

Biomechanical modelling

Individual retro-reflective markers will be attached to anatomical bony landmarks, and marker cluster-shells to limb segments of the thorax, pelvis, bilateral thigh, shank, and foot segments [18]. An eight segment, 27 degrees of freedom (DOFs) model will be constructed from a static standing trial [18]. The position and orientation of each segment will be calculated using an inverse kinematic (IK) algorithm in Visual 3D [41].

Strength assessment (Isokinetic dynamometry)

An isokinetic dynamometer (HUMAC NORM, Computer Sports Medicine Inc., Stoughton, MA) will be used for strength testing, according to manufacturer’s guidelines. For the ankle plantar flexor assessment, participants will lay supine (inclination at 30° above horizontal) with the hip and knee positioned at 60° and 80° of flexion, respectively [42]. Based on the manufacturer’s manual, the dynamometer’s axis will be positioned in line with the lateral malleolus of the test limb. Next, the input arm penetration depth and foot plate penetration depth will be adjusted to approximate the dynamometer’s axis to the ankle axis visually. The testing ankle range will be set from 10° dorsiflexion to 30° plantarflexion. For knee extension testing, participants will be seated in the machine with the hip flexed to 85°. The axis of rotation will be aligned to the femoral condyles with the knee flexed at 90°. Testing range will be set from 0° (complete knee extension) to 90° knee flexion. For all testing, appropriate stabilization of segments will be applied using Velcro straps according to manufacturer’s testing guidelines, and gravity correction mode will be used [43]. Participants will first perform a standardised warm up protocol, consisting of 10 knee extension/flexion repetitions and ankle plantar flexion/dorsiflexion repetitions (90°/s), at a submaximal effort [44]. Concentric knee extensor and ankle plantar flexor torque will be assessed at an angular velocity of 60°/s. For each muscle assessment, two sets of six maximal concentric contractions will be performed, with a between set rest period of 1 min. Between test and between side rest periods of 3 min will be provided.


Familiarization phase (all participants) (two weeks)

All participants will first be enrolled into a two week preparatory training phase prior to baseline testing and randomization. During this phase, all participants will perform the same set of exercises (see Additional file 1 ‘Familiarization phase’). This preparatory phase will control for the effect of motor learning on improvements in performance on the assessments [4547]. Participants will be encouraged not to alter any of their personal training regimes throughout the entire program. Self-reported training for the period of the intervention will be recorded by the participants in a supplemented training diary log book. Variables to report for self-resistance training include external mass magnitude, sets, and repetitions. Variables to report for cardiovascular training include duration, and type of exercise undertaken.

Training phase (six weeks) - standardized warm up (both groups)

Participants from both groups will begin each exercise session with a 15 min warm up consisting of four active, dynamic stretches consisting of 1) lunge, 2) ‘Good Morning’ hamstrings, 3) squats, and 4) calf raises off a step (see Additional file 1 ‘warm up’). Each dynamic stretch will be performed using only the BW as resistance. Each dynamic stretch will involve two sets of ten repetitions [48].

Generic neuromuscular resistance training group (GT group)

The principle governing this training program is that participants perform progressively heavier isoinertial (constant external mass) resistance training, at intensities from approximately 80 % progressing to approximately 88 % of one repetition maximum (1 RM). This program will involve three training sessions per week for six weeks (total 18 sessions). Inter-set rest duration of up to three minutes will be provided [49]. Exercises will include isoinertial bilateral leg press (Cybex® Plate Loaded Squat Press, Cybex International, Inc.), unilateral calf raises (Cybex® Plate Loaded Seated Calf, Cybex International, Inc.), and lunge (Cybex® Plate Loaded Smith Press, Cybex International, Inc.). For the leg press, foot placement will be shoulder width, and the depth of foot placement on the plate will be such that at 90° of knee flexion, the tip of the toes are in line with the knee and shoulders. For the calf raises, the foot will be positioned at the level of the 1st metatarsophalangeal head (i.e. “ball” of the foot). For the lunge, the length of foot placements will be determined as a position that would enable approximately 90° knee flexion of the lead leg at the lowest position of the lunge. The foot of the trailing leg will be positioned such that the trailing knee is slightly posterior to the hip at the lowest position of the lunge. These exercises were selected as they represented generic lower limb exercises used in current load carriage training studies [17]. The intensity, repetitions, sets, rest duration, and description of each exercise will be gradually built up over the six weeks, and is described in the Additional file 1 (see ‘Exercises progression table (General training group)’).

Biomechanically informed neuromuscular resistance training group (BIT group)

The principle of this training program is that key neuromuscular requirements of load carriage running are targeted by specific neuromuscular exercises. This program will involve three training sessions per week for six weeks (total 18 sessions), and will involve SL hopping, CMJ, and hip flexor pull (Cybex® Bravo Pull, Cybex International, Inc.). For the plyometric component (hop and CMJ), intensity will be varied using a weighted vest. A maximal mass of 20 % BM will be added, as previous studies have demonstrated a reduction in peak power with heavier loads [50]. A previously published review indicated that sessions incorporating more than 50 foot contacts per session resulted in the most benefit for jump performance [51]. In order to maintain peak power application for the CMJ and hip flexor pull, a cluster set method (multiple sets of two to three repetitions with 10 s inter-set breaks) will be used [52]. The hip flexor pull will involve a range of 5° hip extension to 90° hip flexion (visual estimation), and a single hand-hold support will be used to maintain balance. For the CMJ, a depth of 80° to 90° knee flexion will be visually estimated and used for all participants. For SL hopping, participants will be encouraged to generate hopping power from the ankle joint, with the knee kept in a relatively ‘isometric’, slightly flexed posture. SL hopping will not involve a cluster-set method as the exercise intrinsically involves continuous, repetitive cycles of fast stretch-shortening cycle. The intensity, repetitions, sets, rest duration, and description of each exercise will be gradually built up over the six weeks, and is listed in the Additional file 1 (see ‘Exercises progression table (Biomechanically informed training group)’).

Augmented feedback (both groups)

Augmented feedback (AF) during all exercises for both groups will be provided to participants, using the principles of motor skill learning [53] (see Additional file). This is to enhance the learning and retention of optimal exercise performance in both groups, with the intention that sub-optimal lower limb kinematics with load may be corrected post-intervention. First, AF that directs an individual’s attention to the consequence of a movement (i.e. external focus of attention) has been shown to result in better motor learning and retention. Second, AF will be provided before (demonstration and instruction), during (mirror feedback and physical/verbal guidance), and after (knowledge of performance) each set of exercise in the initial stages, progressing to feedback delivered only after each exercise set. Feedback based on knowledge of performance (KP) will be provided in a prescriptive sense (i.e. what you should perform) at the initial stages, progressing to descriptive sense (i.e. what was performed) in the later stages of training. This is to allow participants to self-formulate correctional motor strategies in the later stages. The frequency of AF will be reduced from occurring at every set in the early stage of learning, to the last set of an exercise in the later stages. Previous research has shown that introducing a time delay from motor task completion to feedback delivery, especially when participants self-evaluate their performance during this time lag, improves motor skill learning [53].

Determining initial training loads and progression

During the familiarization phase, for all isoinertial exercises a 10 RM will be utilised. A 1 RM load will then be derived from a 10RM load using a regression table for novice strength trainers (Table 2) [54]. Load intensity will progress from 80 % of estimated 1 RM (equivalent to a 10RM load) in the first two weeks, to 84 % of estimated 1RM in the next two weeks (equivalent to an 8RM load), to 88 % of estimated one RM (equivalent to a 6RM load) in the final two weeks. The number of repetitions performed per set will be two repetitions less than the repetition maximum [55]. Estimated 1 RM load for each exercise will be adjusted by a weekly increase of approximately 2.5 % to account for progression in strength. The rate of progression was based on a previous study on time course for strength gains, which demonstrated approximately 20 % increase in measured 1 RM in six weeks [56]. For the hopping and CMJ, a weekly increase in load carried (approximately 5 % BW per week) will be used, until a limit of 20 % BW is reached.

Table 2 Guide for determining one repetition maximum in novice weight trainers [54]

Between group differences in training volume

Training volume as quantified by the number of sets, repetitions, and load magnitude will not be exactly matched between groups. This could mediate any potential between group intervention effects [57]. However, the aim of this study is to test two ecologically realistic model of training on load carriage running mechanical outcomes.

Intervention adherence

Participants attending ≥ 70 % of all training sessions (≥13/18 sessions) will be classified as high adherence, whereas those attending < 70 % will be classified as low adherence. Adherence will be calculated from attendance records in each participant’s exercise training records. Adherence to prescribed neuromuscular training has been previously reported to be an important effect modifier in these programs [58, 59]. Efforts to increase participant adherence include, weekly mobile text (short messaging service) reminders and an exercise diary.

Dependent variables and statistics

For the SJ and CMJ, peak power using inverse dynamics and the force plate approach will be derived [60]. For SL hopping, leg stiffness at each condition will be derived. For the self-paced running tasks with and without load carriage, average self-paced running velocity will be derived over a complete stride. Discrete variables of individual joint positive and negative work, total and net joint work for stance and swing phase of all running trials will be derived. Spatio-temporal variables of stance and swing duration, stride length, and cadence will be derived for all running trials. Time series of the three dimensional joint angles, moments, and powers of all three joints will be extracted for all running trials. Three-dimensional leg stiffness in running will be calculated in a three-step process from an adapted method in a previous study [61]. First, a three dimensional leg length will be defined as the vector from the hip joint centre to the centre of pressure. Second, the component of the resultant three dimensional ground reaction force (GRF) projected onto the leg vector will be calculated (taking the dot product of the GRF vector with the unit vector of the leg). Lastly, leg stiffness will be derived using the ratio of projected GRF (at the time of peak resultant GRF) to the change in leg length (between initial contact to peak GRF). For strength analysis, average peak concentric torque and power, and absolute peak concentric torque and power, of the knee extensors and ankle plantar flexors will be extracted.

Descriptive statistics (mean and standard deviation) will be calculated for baseline demographics of participants. Between groups difference in baseline demographics will be calculated using t-test or non-parametric test where appropriate. Analysis will be based on an intention-to-treat (ITT) using the multiple imputation method [62]. A repeated measures linear mixed model with time, group, and their interaction as fixed effects, and participants clustered within groups as random effects will be used to analyse our discrete dependent variables [63]. For the linear mixed model, significance will be set at α = 0.05. Descriptive statistics and linear mixed modelling will be performed in R software within RStudio (Version 0.98.1062, RStudio, Inc.). Statistical testing between group and within group mechanical wave form data (kinematics, kinetics) will be analysed using Statistical Parametric Mapping (SPM). Statistical significance will be inferred using Random Field Theory (RFT), with appropriate Bonferroni correction applied to retain a family-wise error rate of α = 0.05. SPM will be performed using the latest version of spm1d package (, installed in Python 2.7, and implemented in Enthought Canopy 1.5.4 (Enthought Inc., Austin, USA).


Carrying some form of external load is becoming increasingly ubiquitous in running related sports, such as adventure racing and ultra-endurance events. Running mechanics alterations with load carriage could represent adaptive or mal-adaptive mechanics. Mechanical changes like increased joint power may represent attempts at maintaining constant running velocity, support an increased weight, maintain postural control and/or attenuate excessive impact shocks. In addition, some mechanical changes are likely to represent a failed capacity of lower limb muscles to cope with the additional load, that result in a reduction in running performance and an increased risk of future injuries. The long term sequela not only has an effect at an individual level, but could affect long term sporting participation and health care costs. In addition, runners may have to compromise running economy and running velocity when load carriage is involved if lower limb muscles are not tuned to the specific neuromuscular demands. Velocity decrements as a result of load carriage would result in compromised survivability in combat soldiers, and reduced performance in competing running athletes. This study will provide preliminary evidence of the potential efficacy of a targeted neuromuscular training program or a best-practice strength training program on improvements in strength, stiffness, running velocity and biomechanics during load carriage running.



Augmented feedback


Body weight


Countermovement jumps


Consolidated Standards of reporting trials


Ground reaction force


Knowledge of performance


Repetition maximum


Running related injuries


Squat jumps

SL hop:

Single leg vertical hoping


  1. Grogan R. Run, pheidippides, Run! the story of the battle of marathon. Br J Sports Med. 1981;15(3):186–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Cejka N, Rust CA, Lepers R, Onywera V, Rosemann T, Knechtle B. Participation and performance trends in 100-km ultra-marathons worldwide. J Sports Sci. 2014;32(4):354–66.

    Article  PubMed  Google Scholar 

  3. Joslin J, Worthing R, Black T, Grant WD, Kotlyar T, Wojcik SM. Analysis of weight change and Borg rating of perceived exertion as measurements of runner health and safety during a 6-day, multistage, remote ultramarathon. Clin J Sport Med. 2014;24(3):245–50.

    Article  PubMed  Google Scholar 

  4. Knechtle B, Rosemann T, Zingg MA, Rust CA. Increase in participation but decrease in performance in age group mountain marathoners in the ‘Jungfrau Marathon’: a Swiss phenomenon? Springer Plus. 2015;4:523.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Marais J, de Speville L. Adventure racing. 1st ed. Champaign: Human Kinetics; 2004.

    Google Scholar 

  6. Liew XWB, Morris S, Netto K. Joint power and kinematics coordination in load carriage running: implications for performance and injury. Gait Posture. 2016;47:74–9.

    Article  PubMed  Google Scholar 

  7. Knapik JJ, Reynolds KL, Harman E. Soldier load carriage: historical, physiological, biomechanical, and medical aspects. Mil Med. 2004;169(1):45–56.

    Article  PubMed  Google Scholar 

  8. Orr RM, Johnston V, Coyle J, Pope R. Reported Load Carriage Injuries of the Australian Army Soldier. J Occup Rehabil. 2015;25(2):316–22.

    Article  PubMed  Google Scholar 

  9. Lopes AD, Hespanhol Junior LC, Yeung SS, Costa LO. What are the main running-related musculoskeletal injuries? a systematic review. Sports Med. 2012;42(10):891–905.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Teunissen LP, Grabowski A, Kram R. Effects of independently altering body weight and body mass on the metabolic cost of running. J Exp Biol. 2007;210(Pt 24):4418–27.

    Article  PubMed  Google Scholar 

  11. Carlton SD, Orr RM. The impact of occupational load carriage on carrier mobility: a critical review of the literature. Int J Occup Saf Ergon. 2014;20(1):33–41.

    Article  PubMed  Google Scholar 

  12. Solomonson AA, Dicks ND, Kerr WJ, Pettitt RW. Influence of Load Carriage on High-Intensity Running Performance Estimation. J Strength Cond Res. 2016;30(5):1391–6.

    Article  PubMed  Google Scholar 

  13. Lauersen JB, Bertelsen DM, Andersen LB. The effectiveness of exercise interventions to prevent sports injuries: a systematic review and meta-analysis of randomised controlled trials. Br J Sports Med. 2014;48(11):871–7.

    Article  PubMed  Google Scholar 

  14. Barnes KR, Kilding AE. Strategies to improve running economy. Sports Med. 2015;45(1):37–56.

    Article  PubMed  Google Scholar 

  15. Clansey AC, Hanlon M, Wallace ES, Nevill A, Lake MJ. Influence of tibial shock feedback training on impact loading and running economy. Med Sci Sports Exerc. 2014;46(5):973–81.

    Article  PubMed  Google Scholar 

  16. Barnes KR, McGuigan MR, Kilding AE. Lower-body determinants of running economy in male and female distance runners. J Strength Cond Res. 2014;28(5):1289–97.

    Article  PubMed  Google Scholar 

  17. Knapik JJ, Harman EA, Steelman RA, Graham BS. A systematic review of the effects of physical training on load carriage performance. J Strength Cond Res. 2012;26(2):585–97.

    Article  PubMed  Google Scholar 

  18. Hamner SR, Delp SL. Muscle contributions to fore-aft and vertical body mass center accelerations over a range of running speeds. J Biomech. 2013;46(4):780–7.

    Article  PubMed  Google Scholar 

  19. Hamner SR, Seth A, Delp SL. Muscle contributions to propulsion and support during running. J Biomech. 2010;43(14):2709–16.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Lankhorst NE, Bierma-Zeinstra SM, van Middelkoop M. Factors associated with patellofemoral pain syndrome: a systematic review. Br J Sports Med. 2013;47(4):193–206.

    Article  PubMed  Google Scholar 

  21. Esculier J-F, Roy J-S, Bouyer LJ. Lower limb control and strength in runners with and without patellofemoral pain syndrome. Gait Posture. 2015;41(3):813–9.

    Article  PubMed  Google Scholar 

  22. Schache AG, Brown NA, Pandy MG. Modulation of work and power by the human lower-limb joints with increasing steady-state locomotion speed. J Exp Biol. 2015;218(Pt 15):2472–81.

    Article  PubMed  Google Scholar 

  23. Savovic J, Jones H, Altman D, Harris R, Juni P, Pildal J, Als-Nielsen B, Balk E, Gluud C, Gluud L, et al. Influence of reported study design characteristics on intervention effect estimates from randomised controlled trials: combined analysis of meta-epidemiological studies. Health Technol Assess. 2012;16(35):1–82.

    Article  CAS  PubMed  Google Scholar 

  24. Moher D, Schulz KF, Altman DG. The CONSORT statement: revised recommendations for improving the quality of reports of parallel-group randomised trials. Lancet. 2001;357(9263):1191–4.

    Article  CAS  PubMed  Google Scholar 

  25. Kluitenberg B, van der Worp H, Huisstede BMA, Hartgens F, Diercks R, Verhagen E, van Middelkoop M. The NLstart2run study: Training-related factors associated with running-related injuries in novice runners. J Sci Med Sport. 2016;9(8):642–6.

    Article  Google Scholar 

  26. Dupeyron A, Hertzog M, Micallef JP, Perrey S. Does an abdominal strengthening program influence leg stiffness during hopping tasks? J Strength Cond Res. 2013;27(8):2129–33.

    Article  PubMed  Google Scholar 

  27. Guo Y, Logan HL, Glueck DH, Muller KE. Selecting a sample size for studies with repeated measures. BMC Med Res Methodol. 2013;13:100.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Pruyn EC, Watsford ML, Murphy AJ, Pine MJ, Spurrs RW, Cameron ML, Johnston RJ. Seasonal variation of leg stiffness in professional Australian rules footballers. J Strength Cond Res. 2013;27(7):1775–9.

    Article  PubMed  Google Scholar 

  29. Phinyomark A, Hettinga BA, Osis ST, Ferber R. Gender and age-related differences in bilateral lower extremity mechanics during treadmill running. PLoS ONE. 2014;9(8):e105246.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Doig GS, Simpson F. Randomization and allocation concealment: a practical guide for researchers. J Crit Care. 2005;20(2):187–91. discussion 191–183.

    Article  PubMed  Google Scholar 

  31. Baltich J, Emery C, Stefanyshyn D, Nigg B. The effects of isolated ankle strengthening and functional balance training on strength, running mechanics, postural control and injury prevention in novice runners: design of a randomized controlled trial. BMC Musculoskelet Disord. 2014;15(1):407.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Clarsen B, Myklebust G, Bahr R. Development and validation of a new method for the registration of overuse injuries in sports injury epidemiology: the Oslo sports trauma research centre (OSTRC) overuse injury questionnaire. Br J Sports Med. 2013;47(8):495–502.

    Article  PubMed  Google Scholar 

  33. Zingg M, Rüst CA, Lepers R, Rosemann T, Knechtle B. Master runners dominate 24-h ultramarathons worldwide—a retrospective data analysis from 1998 to 2011. Extrem Physiol Med. 2013;2:21.

    Article  PubMed  PubMed Central  Google Scholar 

  34. McErlain-Naylor S, King M, Pain MT. Determinants of countermovement jump performance: a kinetic and kinematic analysis. J Sports Sci. 2014;32(19):1805–12.

    Article  PubMed  Google Scholar 

  35. Gheller RG, Dal Pupo J, Ache-Dias J, Detanico D, Padulo J, dos Santos SG. Effect of different knee starting angles on intersegmental coordination and performance in vertical jumps. Hum Mov Sci. 2015;42:71–80.

    Article  PubMed  Google Scholar 

  36. Hobara H, Inoue K, Kanosue K. Effect of hopping frequency on bilateral differences in leg stiffness. J Appl Biomech. 2013;29(1):55–60.

    Article  PubMed  Google Scholar 

  37. Van den Bogert A, de Koning J. On optimal filtering for inverse dynamics analysis. In: IXth biennial conference of the Canadian society for biomechanics, August 21–24 1996. British Columbia: Vancouver; 1996.

    Google Scholar 

  38. Cole GK, Nigg BM, Ronsky JL, Yeadon MR. Application of the joint coordinate system to three-dimensional joint attitude and movement representation: a standardization proposal. J Biomech Eng. 1993;115(4A):344–9.

    Article  CAS  PubMed  Google Scholar 

  39. Schache AG, Baker R. On the expression of joint moments during gait. Gait Posture. 2007;25(3):440–52.

    Article  PubMed  Google Scholar 

  40. Huang TWP, Kuo AD. Mechanics and energetics of load carriage during human walking. J Exp Biol. 2014;217(4):605–13.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Lu TW, O’Connor JJ. Bone position estimation from skin marker co-ordinates using global optimisation with joint constraints. J Biomech. 1999;32(2):129–34.

    Article  CAS  PubMed  Google Scholar 

  42. Webber SC, Porter MM. Reliability of ankle isometric, isotonic, and isokinetic strength and power testing in older women. Phys Ther. 2010;90(8):1165–75.

    Article  PubMed  Google Scholar 

  43. de Araujo Ribeiro Alvares JB, Rodrigues R, de Azevedo Franke R, da Silva BGC, Pinto RS, Vaz MA, Baroni BM. Inter-machine reliability of the biodex and cybex isokinetic dynamometers for knee flexor/extensor isometric, concentric and eccentric tests. Phys Ther Sport. 2015;16(1):59–65.

    Article  PubMed  Google Scholar 

  44. Baroni BM, Rodrigues R, Franke RA, Geremia JM, Rassier DE, Vaz MA. Time course of neuromuscular adaptations to knee extensor eccentric training. Int J Sports Med. 2013;34(10):904–11.

    Article  CAS  PubMed  Google Scholar 

  45. Feltner ME, MacRae PG. Time course of changes in novice jumpers’ countermovement vertical jump performance. Percept Mot Skills. 2011;112(1):228–42.

    Article  PubMed  Google Scholar 

  46. Amarante do Nascimento M, Januario RS, Gerage AM, Mayhew JL, Cheche Pina FL, Cyrino ES. Familiarization and reliability of one repetition maximum strength testing in older women. J Strength Cond Res. 2013;27(6):1636–42.

    Article  PubMed  Google Scholar 

  47. Glass SC. Effect of a learning trial on self-selected resistance training load. J Strength Cond Res. 2008;22(3):1025–9.

    Article  PubMed  Google Scholar 

  48. Yessis M. Explosive plyometrics: ultimate athlete concepts. 2009.

    Google Scholar 

  49. de Salles BF, Simao R, Miranda F, Novaes Jda S, Lemos A, Willardson JM. Rest interval between sets in strength training. Sports Med. 2009;39(9):765–77.

    Article  PubMed  Google Scholar 

  50. Suzovic D, Markovic G, Pasic M, Jaric S. Optimum load in various vertical jumps support the maximum dynamic output hypothesis. Int J Sports Med. 2013;34(11):1007–14.

    Article  CAS  PubMed  Google Scholar 

  51. de Villarreal ES, Kellis E, Kraemer WJ, Izquierdo M. Determining variables of plyometric training for improving vertical jump height performance: a meta-analysis. J Strength Cond Res. 2009;23(2):495–506.

    Article  PubMed  Google Scholar 

  52. Moreno SD, Brown LE, Coburn JW, Judelson DA. Effect of cluster sets on plyometric jump power. J Strength Cond Res. 2014;28(9):2424–8.

    Article  PubMed  Google Scholar 

  53. Schmidt RA, Lee TD editors. Motor learning and performance : a situation-based learning approach 5edn. Champaign, IL: Human Kinetics; 2014.

  54. Dan Baker Strength and Power Training. [].

  55. Sampson JA, Groeller H. Is repetition failure critical for the development of muscle hypertrophy and strength? Scand J Med Sci Sports. 2016;26(4):375–83.

    Article  CAS  PubMed  Google Scholar 

  56. Abe T, DeHoyos DV, Pollock ML, Garzarella L. Time course for strength and muscle thickness changes following upper and lower body resistance training in men and women. Eur J Appl Physiol. 2000;81(3):174–80.

    Article  CAS  PubMed  Google Scholar 

  57. Krieger JW. Single versus multiple sets of resistance exercise: a meta-regression. J Strength Cond Res. 2009;23(6):1890–901.

    Article  PubMed  Google Scholar 

  58. Sugimoto D, Myer GD, Bush HM, Hewett TE. Effects of compliance on trunk and hip integrative neuromuscular training on hip abductor strength in female athletes. J Strength Cond Res. 2014;28(5):1187–94.

    PubMed  PubMed Central  Google Scholar 

  59. Sugimoto D, Myer GD, Bush HM, Klugman MF, Medina McKeon JM, Hewett TE. Compliance with neuromuscular training and anterior cruciate ligament injury risk reduction in female athletes: a meta-analysis. J Athl Train. 2012;47(6):714–23.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Jandacka D, Uchytil J, Farana R, Zahradnik D, Hamill J. Lower extremity power during the squat jump with various barbell loads. Sports Biomech. 2014;13(1):75–86.

    Article  PubMed  Google Scholar 

  61. Coleman DR, Cannavan D, Horne S, Blazevich AJ. Leg stiffness in human running: comparison of estimates derived from previously published models to direct kinematic-kinetic measures. J Biomech. 2012;45(11):1987–91.

    Article  PubMed  Google Scholar 

  62. Powney M, Williamson P, Kirkham J, Kolamunnage-Dona R. A review of the handling of missing longitudinal outcome data in clinical trials. Trials. 2014;15:237.

    Article  PubMed  PubMed Central  Google Scholar 

  63. West B, Welch K, Galecki A, editors. Linear mixed models: a practical guide using statistical software. 2nd ed. Boca Raton: Chapman and Hall/CRC; 2014.

    Google Scholar 

  64. Siegel KL, Kepple TM, Stanhope SJ. Joint moment control of mechanical energy flow during normal gait. Gait Posture. 2004;19(1):69–75.

    Article  PubMed  Google Scholar 

  65. Lai A, Schache AG, Lin YC, Pandy MG. Tendon elastic strain energy in the human ankle plantar-flexors and its role with increased running speed. J Exp Biol. 2014;217(Pt 17):3159–68. Accessed 1 Feb 2016.

    Article  PubMed  Google Scholar 

  66. Dorn TW, Schache AG, Pandy MG. Muscular strategy shift in human running: dependence of running speed on hip and ankle muscle performance. J Exp Biol. 2012;215(Pt 11):1944–56.

    Article  PubMed  Google Scholar 

  67. Brughelli M, Cronin J. A review of research on the mechanical stiffness in running and jumping: methodology and implications. Scand J Med Sci Sports. 2008;18(4):417–26.

    Article  CAS  PubMed  Google Scholar 

  68. Caron RR, Wagenaar RC, Lewis CL, Saltzman E, Holt KG. Center of mass trajectory and orientation to ankle and knee in sagittal plane is maintained with forward lean when backpack load changes during treadmill walking. J Biomech. 2013;46(1):70–6.

    Article  PubMed  Google Scholar 

  69. Moore IS, Jones AM, Dixon SJ. Reduced oxygen cost of running is related to alignment of the resultant GRF and leg axis vector: A pilot study. Scand J Med Sci Sports. 2016;26(7):809–15.

    Article  CAS  PubMed  Google Scholar 

  70. Silder A, Besier T, Delp SL. Running with a load increases leg stiffness. J Biomech. 2015;48(6):1003–8.

    Article  PubMed  Google Scholar 

  71. Liew XWB, Morris S, Netto K. The effects of load carriage on joint work at different running velocities. J Biomech. 2016; doi:10.1016/j.jbiomech.2016.08.012. [Epub ahead of print].

  72. Dierks TA, Manal KT, Hamill J, Davis IS. Proximal and distal influences on hip and knee kinematics in runners with patellofemoral pain during a prolonged run. J Orthop Sports Phys Ther. 2008;38(8):448–56.

    Article  PubMed  Google Scholar 

  73. Bonacci J, Saunders PU, Hicks A, Rantalainen T, Vicenzino BG, Spratford W. Running in a minimalist and lightweight shoe is not the same as running barefoot: a biomechanical study. Br J Sports Med. 2013;47(6):387–92.

    Article  PubMed  Google Scholar 

  74. Rooney BD, Derrick TR. Joint contact loading in forefoot and rearfoot strike patterns during running. J Biomech. 2013;46(13):2201–6.

    Article  PubMed  Google Scholar 

Download references


Not applicable.


No funds were received in support of this work. No benefits in any form have been or will be received from a commercial party related directly or indirectly to the subject of this manuscript. Mr Bernard Liew is currently supported by an institutional doctoral scholarship.

Availability of data and materials

Not applicable.

Authors’ contributions

BL, SM, and KN conceived of the study and participated in its design and coordination. JK and BA provided in depth contributions to the design of the resistance training program and the assessment used in this study. BL is responsible for the coordination of the study, and is involved with the assessment of all participants. BL was responsible for training the trainers in conducting the resistance training program. All five authors were responsible for providing statistical supervision for this study. All five authors were involved in drafting the protocol, and contributed to and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Written informed consent was provided by the participant for publication of their image in the Additional file of this manuscript.

Ethics approval and consent to participate

This study protocol design follows the Declaration of Helsinki and is approved by the Curtin University Human Research Ethics Committee (RD-41-14). All participants will provide written informed consent prior to study inclusion.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Bernard X. W. Liew.

Additional file

Additional file 1:

Effects of two neuromuscular training programs on running biomechanics with load carriage: a randomized controlled trial – a study protocol. (DOCX 312 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liew, B.X.W., Morris, S., Keogh, J.W.L. et al. Effects of two neuromuscular training programs on running biomechanics with load carriage: a study protocol for a randomised controlled trial. BMC Musculoskelet Disord 17, 445 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Running biomechanics
  • Load carriage
  • Randomized clinical trial
  • Neuromuscular training
  • Resistance training