Abstract
ideo-based biomechanical running gait analysis is widely used to optimise technique, guide footwear selection, and identify orthopaedic risk factors. Despite the increasing availability of such assessments, it is often assumed—without strong empirical support—that key kinematic parameters of running gait remain stable over short periods of time. This study aimed to examine the short-term stability of key joint angles during running using a standard 2D video-based kinematic analysis. Specifically, it was investi- gated whether these angles change within the first 4 min of treadmill running under three defined conditions: barefoot at 12 km h −1 , shoed at 12 km h −1 , and shoed at 14 km h −1 , in a homogeneous sample of twelve young, trained, male recreational soccer players.Methods: Participants completed three four-minute runs. Joint angles were quantified manually from 2D video recordings. Temporal variation was analysed using repeated-measures statistics, intraclass correlation coefficients (ICCs), and minimal detectable changes (MDCs).Results: Six out of nine joint angles showed statistically significant temporal changes, mainly in hip extension, knee flexion, the Duchenne angle, the Trendelenburg angle, the leg axis angle, and heel-bottom angle. Lower leg angle and Achilles tendon angle remained stable. ICCs showed moderate to excellent agreement, indicating high within-session consistency across all angles.Discussion: Under the applied study protocol, significant short-term
(MDCs).Results: Six out of nine joint angles showed statistically significant temporal changes, mainly in hip extension, knee flexion, the Duchenne angle, the Trendelenburg angle, the leg axis angle, and heel-bottom angle. Lower leg angle and Achilles tendon angle remained stable. ICCs showed moderate to excellent agreement, indicating high within-session consistency across all angles.Discussion: Under the applied study protocol, significant short-term variations were observable in several joint angles during the first four minutes of running. These findings highlight the importance of analysing multiple strides and considering measurement uncertainty when interpreting short-duration running kinematics. Keywords:kinematics; 2D motion analysis; treadmill; footwear; running speed 1. Introduction Endurance running is one of the most popular recreational sport activities in Central Europe and North America [1,2]. It has been demonstrated to exert positive effects on the cardiovascular system, reducing resting pulse and blood pressure, and strengthening the immune system [2,3]. However, running is also associated with a high prevalence of overuse injuries, particularly in the lower extremities [1,4,5]. Biomechanical running gait analyses play a pivotal role in preventing such overuse injuries by helping to identify imbalances and improving load distribution and efficiency of movements [4,6]. In addition, Appl. Sci.2026,16, 133 https://doi.org/10.3390/app16010133
Appl. Sci.2026,16, 133 2 of 13 running gait analyses may help to find adequate footwear and optimise performance on an individual level [2]. There are numerous approaches to conducting a running gait assessment, with optical video analysis motion capture traditionally considered the gold standard [2]. As technology advances, three-dimensional (3D) optical motion analysis systems have been increasingly utilised to provide more detailed and comprehensive insights into running kinematics than two-dimensional (2D) optical systems. Conventional 3D motion capture systems are most often based on infrared reflective optics and are as such typically elaborate, costly, not portable and only applicable to indoor environments, hindering their in-field use outside the laboratory [2]. In the quest to overcome these general limitations of traditional optical motion capture, advances in artificial intelligence (AI) applied to easily accessible smartphone-based video analysis have recently gained traction as a more cost-effective and accessible alternative [7,8]. Regardless of the measurement system used, the basis of any running gait analysis is the evaluation of spatio-temporal biomechanical parameters at specific phases of the running gait cycle. However, defining the precise conditions under which these parameters are to be assessed proves challenging, given that runners continually adapt their running gait in response to internal and external factors, such as intensity and environmental condi- tions [9,10]. Consequently, it is imperative to establish explicit criteria to ensure that the findings of a running gait analysis are translatable into practical exercise recommenda- tions. The literature identifies three main influencing factors of running gait: duration and distance [11–13], intensity [9], and the surface of an overground running session [10,14,15]. As for distance, it has been observed, for instance, that hip extension during the swing phase decreases over the course of a marathon, while knee flexion increases during the swing phase and decreases during stance over that distance [13]. Step protocols until voluntary exhaustion have also been shown to result in a reduced knee flexion angle at initial contact and during the stance phase. However, these changes typically did not exceed the threshold of 5% [11]. Similarly, alterations in kinematic parameters have also been observed at shorter distances.
the swing phase and decreases during stance over that distance [13]. Step protocols until voluntary exhaustion have also been shown to result in a reduced knee flexion angle at initial contact and during the stance phase. However, these changes typically did not exceed the threshold of 5% [11]. Similarly, alterations in kinematic parameters have also been observed at shorter distances. As reported by Williams et al. [16], a significant increase in knee flexion angle may be seen during the initial and final phases of 1500 to 5000 m treadmill and track runs. Additionally, Derrick et al. [12] reported a time-dependent change in knee flexion angle at initial contact (IC) and peak knee angle during the stance phase after a 3200 m treadmill or track run. Running intensity is usually expressed in terms of running speed or pace and may be influenced by external conditions, including incline, climate, and altitude. Running speed has been found to affect, among other parameters, leg stiffness as well as ankle and knee forces [17]. For joint angles, however, conflicting results have been reported: While some studies found no evidence for an effect of running speed on sagittal joint angles [18], others did report changes in ankle, knee, and hip flexion angles as a result of altered running speed [9,19]. The running surface is widely known to influence running biomechanics [14,15]. For instance, grass or woodchip tracks induce significantly lower peak tibial accelerations than synthetic tracks or concrete roads [5,14,20]. Moreover, treadmill running cannot fully simulate actual overground running: While overground running results in greater hip flexion at foot strike and greater ankle dorsiflexion at stance, it shows reduced ankle eversion [21]. In addition, several studies have demonstrated biomechanical differences between barefoot running and shoed running. Barefoot running typically results in reduced peak ground reaction forces (GRF), greater foot and ankle plantarflexion, and increased knee flexion at initial contact [18,22,23]. https://doi.org/10.3390/app16010133
Appl. Sci.2026,16, 133 3 of 13 In light of these findings, it is clear that a valid implementation of a running gait analysis should include different durations, intensities and surfaces. However, in real- world testing scenarios, this scientifically based requirement is often not compatible with practical economic constraints. For example, in typical running gear stores, running gait analysis is often reduced to a short session of treadmill running at a single low to moderate pace. Even if equipped with advanced AI-based instant video analysis, modern retail stores rarely spend time evaluating more than 2–5 consecutive strides to obtain the key spatio-temporal metrics that form the basis of biomechanical running gait analysis. Interestingly, even in several recent studies that were conducted in well-equipped modern kinesiology laboratories, only 12–15 consecutive steps were assessed [9,24], despite the fact that technological capacity would have allowed for a deeper data base and neither time nor cost constraints were likely to be relevant in this context. In essence, there often seems to be the implicit notion that for a well-defined running condition at submaximal intensity, the running gait is expected to be stable, suggesting that it does change significantly over several minutes. Despite the apparent importance of the validity of running gait analyses, there is no explicit scientific evidence for this assumption in the extant literature. Therefore, the aim of the present study was to examine the short-term stability of key joint angles during treadmill running using standard 2D video-based kinematic analysis. Based on initial findings on early-phase adaptations and footwear-related effects, it was hypothesised that joint angles might show small but detectable short-term changes already during the first minutes of running. We further hypothesised that these temporal adjustments might differ between barefoot and shoed running as well as between moderate and faster running speeds. Nonetheless, given the short duration of each running bout, we expected overall intra-individual consistency to remain relatively high. Specifically, we investigated whether these angles change within the first 4 min under three defined running conditions: barefoot at 12 km h −1 , shoed at 12 km h −1 , and shoed
as well as between moderate and faster running speeds. Nonetheless, given the short duration of each running bout, we expected overall intra-individual consistency to remain relatively high. Specifically, we investigated whether these angles change within the first 4 min under three defined running conditions: barefoot at 12 km h −1 , shoed at 12 km h −1 , and shoed at 14 km h −1 . Based on our initial assumptions, we expected small but detectable short-term changes in key joint angles during the first minutes of running. Specifically, we hypothesised that (1) in the sagittal plane, knee flexion at mid-stance would show a gradual increase over time whereas hip extension at toe-off would exhibit a gradual decrease, (2) frontal-plane angles at mid-stance would moderately increase with time, and (3) these changes would be most pronounced in the faster shoed condition (14 km h −1 ). 2. Materials and Methods 2.1. Participants Twelve male amateur soccer players aged 22 to 36 years (body mass 86.9±13.2 kg; height: 180.1±6.0 cm) participated in this study. All participants were active soccer players of the regional team competing in the state league of Saxony Anhalt (7th league in Germany) and reported to be free of injuries and illnesses. This cohort was chosen to reflect a typical target group of healthy, non-specialist recreational athletes, who may occasionally seek counselling on optimal footwear choices in sport equipment stores or request basic biomechanical assessment at orthopaedic facilities. The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics committee of the Department of Engineering and Industrial Design of the Magdeburg-Stendal University of Applied Sciences (Approval number: EKIWID-2024-11-001). All subjects gave their written informed consent prior to participation and were free to withdraw from the study at any time. https://doi.org/10.3390/app16010133
Appl. Sci.2026,16, 133 4 of 13 2.2. Study Design The employed video-based biomechanical assessment was based on a traditional 2D video analysis in the frontal and sagittal planes. The running movement was recorded by two cameras (Canon XH G1S; Canon XH G1, Tokyo, Japan, 50 Hz, 1920×1080 pixels, FHD) placed according to Marquardt [25], one laterally on the runner’s left side and one behind the treadmill, both positioned perpendicular to the runner’s centre of mass at a distance of 2 m. The treadmill used in this study was a motorised h/p/cosmos Saturn ® (250/100) model (h/p/cosmos sports & medical GmbH, Nußdorf-Traunstein, Germany). For each joint angle, clearly defined body or limb positions were analysed as illustrated in Figure. In the sagittal plane, angles were measured at initial contact, at mid-stance, at toe-off, and during the swing phase. In the frontal plane, all angles were evaluated during mid-stance as viewed from behind. These definitions follow the standard 4-point method by Marquardt [25]. Participants were instructed to wear tight shorts and their preferred running shoes during the shoed trials. To facilitate video-based joint angle monitoring, 17 optical mark- ers were placed on the participants’ body and shorts according to Marquardt’s 4-point method [25]. After preparation, participants completed a standardised warm-up consisting of a 4-min shoed treadmill run at a constant speed of 10 km h −1 . The test protocol consisted of three 4-min bouts of treadmill running, separated by 4-min periods of rest. The warm-up duration of 4 min was selected to align with common practice in clinical and commercial treadmill gait assessments, where only brief (and often even shorter) familiarisation periods are typically used. This duration allowed participants to reach a comfortable steady-state running pattern while minimising pre-fatigue prior to data collection. All bouts were performed in a fixed order, beginning with the running condition (1) barefoot at a speed of 12 km h −1 (Barefoot12), followed by (2) shoed at 12 km h −1 (Shoed12), and (3) shoed at 14 km h −1 (Shoed14), reflecting standard running analysis protocols used for amateur athletes. Perceived exertion was assessed at the
data collection. All bouts were performed in a fixed order, beginning with the running condition (1) barefoot at a speed of 12 km h −1 (Barefoot12), followed by (2) shoed at 12 km h −1 (Shoed12), and (3) shoed at 14 km h −1 (Shoed14), reflecting standard running analysis protocols used for amateur athletes. Perceived exertion was assessed at the start of each bout using the Borg scale [26]. Following the common 4-point method of Marquardt [25], the videos were analysed by the same researcher using a commercial video analysis software (utilius fairplay 5, ccc software, Leipzig, Germany) to evaluate four joint angles in the sagittal plane and five angles in the frontal plane, all measured in degrees ( ◦ ). The sagittal joint angles included the lower leg angle at initial contact (LLIC), knee flexion angle in mid-stance (KFMS), hip extension angle at toe-off (HETO), and lower leg angle during swing phase (LLSP) (Figurea–d). The frontal plane angles analysed were the Duchenne angle (DA), Trendelenburg angle (TB), leg axis angle (LA), Achilles tendon angle (AT) and heel-bottom angle (HB) (Figuree–i). 2.3. Statistical Analysis To determine the nine joint angles, eight contiguous segments [ti,ti+1], i = 1, 2, . . . , 8, of each 4-min run were evaluated, starting att1= 30 s and progressing in steps of 20 s up to t8= 190 s, i.e., [30 s, 50 s], [50 s, 70 s], [70 s, 90 s] etc. For each segment, the first five strides were evaluated and averaged for subsequent analysis. Data curation was performed using MATLAB (R2023a, MathWorks Inc., Natick, MA, USA), while all statistical analyses were performed using JASP (version 0.17.1, JASP Team (2023), Amsterdam, The Netherlands). For all statistical tests, the significance level was set top< 0.05. Two-way repeated measures ANOVA with Greenhouse-Geisser correction was used to analyse whether the running condition (factor 1) and/or time in terms of segment number i (factor 2) influenced the outcome of individual joint angles. If a significant main effect was observed, effect sizes were assessed using partial eta squared (ηp 2 ), and post-hoc tests were performed. Results
top< 0.05. Two-way repeated measures ANOVA with Greenhouse-Geisser correction was used to analyse whether the running condition (factor 1) and/or time in terms of segment number i (factor 2) influenced the outcome of individual joint angles. If a significant main effect was observed, effect sizes were assessed using partial eta squared (ηp 2 ), and post-hoc tests were performed. Results https://doi.org/10.3390/app16010133
Appl. Sci.2026,16, 133 5 of 13 are presented as mean±standard deviation (SD) unless otherwise stated. All variables were analysed using an intraclass correlation coefficient (ICC(1, 3)), indicating sufficiently high reliability, with all ICC values exceeding 0.75. For further evaluation of the results, the minimal detectable change (MDC), indicating the smallest measurable difference that can be attributed to the effect of time rather than measurement error, was calculated according to MDC=1.96 SE× √ 2, with SE denoting the standard error [27]. Sagittal plane (a) ( b) ( c) ( d) Lower leg angle at initial contact (LLIC) Knee flexion angle in mid-stance (KFMS) Hip extension angle at toe-off (HETO) Lower leg angle during swing phase (LLSP) Frontal plane (e) ( f) ( g) ( h) ( i) Duchenne angle (DA) Trendelenburg angle (TB) Leg axis angle (LA) Achilles tendon angle (AT) Heel-bottom angle (HB) Figure 1.Definition of the analysed joint angles using a model representative of the study pop- ulation. Joint angles are visualised by a blue reference line and an orange line representing the measured segment orientation used to determine the respective angle. Subfigures (a–d) depict sagittal-plane angles: (a) lower leg angle at initial contact (LLIC), (b) knee flexion angle in mid-stance (KFMS), (c) hip extension angle at toe-off (HETO), and (d) lower leg angle during the swing phase (LLSP). Subfigures (e–i) show frontal-plane angles recorded from behind: (e) Duchenne angle (DA), (f) Trendelenburg angle (TB), (g) leg axis angle (LA), (h) Achilles tendon angle (AT), and (i) heel- bottom angle (HB). The camera positions may have been adjusted for visual clarity. 3. Results Table the time points analysed, together with the number of soccer players included in the respective angle measurement. While eight players could be included for the analysis of the https://doi.org/10.3390/app16010133
Appl. Sci.2026,16, 133 6 of 13 Duchenne angle, 11 players were included for the Achilles tendon and heel-bottom angles. It is noteworthy that the minimum and maximum mean values for approximately 85% of the angles measured occur either at the first two (t1ort2) or last two time points (t7ort8). These minima are marked in blue colour in Table, while the maxima are highlighted in green colour. In contrast to the other angles, the heel-bottom angle shows its maximum at the beginning of the measurements (t1) and its minimum at the end of the run (t7). Looking at the median values, the tendency of an increase of angles over time is less obvious. Minima are mostly observed during the first four time points (t1,t2,t3,t4), while maxima are present in the last four time points (t5,t6,t7,t8). Figures during the first 4 min of each running condition. Qualitative assessment shows that some angles exhibit slight increasing or decreasing trends over time, whereas others remain relatively stable. The corresponding quantitative examination of stability over time confirms significant temporal changes in several joint angles (Table). Specifically, significant time effects were observed for KFMS, HETO, DA, TB and LA (allp< 0.05), whereas LLIC, LLSP, AT and HB showed no significant temporal change. In contrast, the running condition had a significant effect on almost all angles except for HB, with substantial effects observed for LLIC, KFMS, HETO, LLSP, DA, TB, LA, and AT. Importantly, no interaction effects of time×running condition were detected for any joint angle (all p> 0.14), indicating that temporal changes did not differ between the three conditions. Post-hoc analyses showed that temporal increases were most pronounced for HETO and TB, while KFMS and DA exhibited smaller but still detectable changes over time. Figure 2.Experimentally observed joint angles in the sagittal plane as a function of the three running conditions: Barefoot at 12 km h −1 (green), shoed at 12 km h −1 (blue) and shoed at 14 km h −1 (red). The angles shown are the lower leg angle at initial contact (LLIC), the knee flexion angle in mid-stance (KFMS), hip extension angle at
observed joint angles in the sagittal plane as a function of the three running conditions: Barefoot at 12 km h −1 (green), shoed at 12 km h −1 (blue) and shoed at 14 km h −1 (red). The angles shown are the lower leg angle at initial contact (LLIC), the knee flexion angle in mid-stance (KFMS), hip extension angle at toe-off (HETO) and the lower leg angle during swing phase (LLSP). Red crosses mark outliers (outside the whiskers representing 1.5 times the interquartile range). https://doi.org/10.3390/app16010133
Description
The study investigates joint angle stability in running under different conditions.