Todd C. Pataky

The effects of habitual footwear use: foot shape and function in native barefoot walkers†

The human foot was anatomically modern long before footwear was invented, and is adapted to barefoot walking on natural substrates. Understanding the biomechanics of habitually barefoot walkers can provide novel insights both for anthropologist and for applied scientists, yet the necessary data is virtually non-existent. To start assessing morphological and functional effects of the habitual use of footwear, we have studied a population of habitually barefoot walkers from India (n = 70), and compared them with a habitually shod Indian control group (n = 137) and a Western population (n = 48). We focused on foot metrics and on the analysis of plantar pressure data, which was performed using a novel, pixel based method (Pataky and Goulermas 2008, Journal of Biomechanics, 41, 2136). Habitually shod Indians wore less often, and less constricting shoes than Western people. Yet, we found significant differences with their habitually barefoot peers, both in foot shape and in pressure distribution. Barefoot walkers had wider feet and more equally distributed peak pressures, i.e. the entire load carrying surface was contributing more uniformly than in habitually shod subjects, where regions of very high or very low peak pressures were more apparent. Western subjects differed strongly from both Indian populations (and most from barefoot Indians), by having relatively short and, especially, slender feet, with more focal and higher peak pressures at the heel, metatarsals and hallux. The evolutionary history of humans shows that barefoot walking is the biologically natural situation. The use of footwear remains necessary, especially on unnatural substrates, in athletics, and in some pathologies, but current data suggests that footwear that fails to respect natural foot shape and function will ultimately alter the morphology and the biomechanical behaviour of the foot.

New insights into the plantar pressure correlates of walking speed using pedobarographic statistical parametric mapping (pSPM)

This study investigates the relation between walking speed and the distribution of peak plantar pressure and compares a traditional ten-region subsampling (10RS) technique with a new technique: pedobarographic statistical parametric mapping (pSPM). Adapted from cerebral fMRI methodology, pSPM is a digital image processing technique that registers foot pressure images such that homologous structures optimally overlap, thereby enabling statistical tests to be conducted at the pixel level. Following previous experimental protocols, we collected pedobarographic records from 10 subjects walking at three different speeds: slow, normal, and fast. Walking speed was recorded and correlated with the peak pressures extracted from the 10 regions, and subsequently with the peak pixel data extracted after pSPM preprocessing. Both methods revealed significant positive correlation between peak plantar pressure and walking speed over the rearfoot and distal forefoot after Bonferroni correction for multiple comparisons. The 10RS analysis found positive correlation in the midfoot and medial proximal forefoot, but the pixel data exhibited significant negative correlation throughout these regions (p<5x10(-5)). Comparing the statistical maps from the two approaches shows that subsampling may conflate pressure differences evident in pixel-level data, obscuring or even reversing statistical trends. The negative correlation observed in the midfoot implies reduced longitudinal arch collapse with higher walking speeds. We infer that this results from pre- or early-stance phase muscle activity and speculate that preferred walking speed reflects, in part, a balance between the energy required to tighten the longitudinal arch and the apparent propulsive benefits of the stiffened arch.

Human-like external function of the foot, and fully upright gait, confirmed in the 3.66 million year old Laetoli hominin footprints by topographic statistics, experimental footprint-formation and computer simulation

It is commonly held that the major functional features of the human foot (e.g. a functional longitudinal medial arch, lateral to medial force transfer and hallucal (big-toe) push-off) appear only in the last 2 Myr, but functional interpretations of footbones and footprints of early human ancestors (hominins) prior to 2 million years ago (Mya) remain contradictory. Pixel-wise topographical statistical analysis of Laetoli footprint morphology, compared with results from experimental studies of footprint formation; foot-pressure measurements in bipedalism of humans and non-human great apes; and computer simulation techniques, indicate that most of these functional features were already present, albeit less strongly expressed than in ourselves, in the maker of the Laetoli G-1 footprint trail, 3.66 Mya. This finding provides strong support to those previous studies which have interpreted the G-1 prints as generally modern in aspect.

Functional Evolution of the Feeding System in Rodents

The masticatory musculature of rodents has evolved to enable both gnawing at the incisors and chewing at the molars. In particular, the masseter muscle is highly specialised, having extended anteriorly to originate from the rostrum. All living rodents have achieved this masseteric expansion in one of three ways, known as the sciuromorph, hystricomorph and myomorph conditions. Here, we used finite element analysis (FEA) to investigate the biomechanical implications of these three morphologies, in a squirrel, guinea pig and rat. In particular, we wished to determine whether each of the three morphologies is better adapted for either gnawing or chewing. Results show that squirrels are more efficient at muscle-bite force transmission during incisor gnawing than guinea pigs, and that guinea pigs are more efficient at molar chewing than squirrels. This matches the known diet of nuts and seeds that squirrels gnaw, and of grasses that guinea pigs grind down with their molars. Surprisingly, results also indicate that rats are more efficient as well as more versatile feeders than both the squirrel and guinea pig. There seems to be no compromise in biting efficiency to accommodate the wider range of foodstuffs and the more general feeding behaviour adopted by rats. Our results show that the morphology of the skull and masticatory muscles have allowed squirrels to specialise as gnawers and guinea pigs as chewers, but that rats are high-performance generalists, which helps explain their overwhelming success as a group.

Generalized n-dimensional biomechanical field analysis using statistical parametric mapping

A variety of biomechanical data are sampled from smooth n-dimensional spatiotemporal fields. These data are usually analyzed discretely, by extracting summary metrics from particular points or regions in the continuum. It has been shown that, in certain situations, such schemes can compromise the spatiotemporal integrity of the original fields. An alternative methodology called statistical parametric mapping (SPM), designed specifically for continuous field analysis, constructs statistical images that lie in the original, biomechanically meaningful sampling space. The current paper demonstrates how SPM can be used to analyze both experimental and simulated biomechanical field data of arbitrary spatiotemporal dimensionality. Firstly, 0-, 1-, 2-, and 3-dimensional spatiotemporal datasets derived from a pedobarographic experiment were analyzed using a common linear model to emphasize that SPM procedures are (practically) identical irrespective of the datas physical dimensionality. Secondly two probabilistic finite element simulation studies were conducted, examining heel pad stress and femoral strain fields, respectively, to demonstrate how SPM can be used to probe the significance of field-wide simulation results in the presence of uncontrollable or induced modeling uncertainty. Results were biomechanically intuitive and suggest that SPM may be suitable for a wide variety of mechanical field applications. SPMs main theoretical advantage is that it avoids problems associated with a priori assumptions regarding the spatiotemporal foci of field signals. SPMs main practical advantage is that a unified framework, encapsulated by a single linear equation, affords comprehensive statistical analyses of smooth scalar fields in arbitrarily bounded n-dimensional spaces.

One-dimensional statistical parametric mapping in Python

Statistical parametric mapping (SPM) is a topological methodology for detecting field changes in smooth n-dimensional continua. Many classes of biomechanical data are smooth and contained within discrete bounds and as such are well suited to SPM analyses. The current paper accompanies release of ‘SPM1D’, a free and open-source Python package for conducting SPM analyses on a set of registered 1D curves. Three example applications are presented: (i) kinematics, (ii) ground reaction forces and (iii) contact pressure distribution in probabilistic finite element modelling. In addition to offering a high-level interface to a variety of common statistical tests like t tests, regression and ANOVA, SPM1D also emphasises fundamental concepts of SPM theory through stand-alone example scripts. Source code and documentation are available at: www.tpataky.net/spm1d/.

A dynamic model of the windlass mechanism of the foot: evidence for early stance phase preloading of the plantar aponeurosis

SUMMARY In the present study we have estimated the temporal elongation of the plantar aponeurosis (PA) during normal walking using a subject-specific multi-segment rigid-body model of the foot. As previous studies have suggested that muscular forces at the ankle can pre-load the PA prior to heel-strike, the main purpose of the current study was to test, through modelling, whether there is any tension present in the PA during early stance phase. Reflective markers were attached to bony landmarks to track the kinematics of the calcaneus, metatarsus and toes during barefoot walking. Ultrasonography measurements were performed on three subjects to determine both the location of the origin of the PA on the plantar aspect of the calcaneus, and the radii of the metatarsal heads. Starting with the foot in a neutral, unloaded position, inverse kinematics allowed calculation of the tension in the five slips of the PA during the whole duration of the stance phase. The results show that the PA experienced tension significantly above rest during early stance phase in all subjects (P<0.01), thus providing support for the PA-preloading hypothesis. The amount of preloading and the maximum elongation of the slips of the PA decreased from medial to lateral. The mean maximum tension exerted by the PA was 1.5 BW (body weight) over the three subjects.

Dynamics of longitudinal arch support in relation to walking speed: contribution of the plantar aponeurosis

The plantar aponeurosis (PA), in spanning the whole length of the plantar aspect of the foot, is clearly identified as one of the key structures that is likely to affect compliance and stability of the longitudinal arch. A recent study performed in our laboratory showed that tension/elongation in the PA can be predicted from the kinematics of the segments to which the PA is attached. In the present investigation, stereophotogrammetry and inverse kinematics were employed to shed light on the mechanics of the longitudinal arch and its main passive stabilizer, the PA, in relation to walking speed. When compared with a neutral unloaded position, the medial longitudinal arch underwent greater collapse during the weight‐acceptance phase of stance at higher walking speed (0.1°±1.9° in slow walking; 0.9°±2.6° in fast walking; P = 0.0368). During late stance the arch was higher (3.4°±3.1° in slow walking; 2.8°±2.7° in fast walking; P = 0.0227) and the metatarsophalangeal joints more dorsiflexed (e.g. at the first metatarsophalangeal joint, 52°±5° in slow walking; 64°±4° in fast walking; P < 0.001) during fast walking. Early‐stance tension in the PA increased with speed, whereas maximum tension during late stance did not seem to be significantly affected by walking speed. Although, on the one hand, these results give evidence for the existence of a pre‐heel‐strike, speed‐dependent, arch‐stiffening mechanism, on the other hand they suggest that augmentation of arch height in late stance is enhanced by higher forces exerted by the intrinsic muscles on the plantar aspect of the foot when walking at faster speeds.

The effect of running speed on knee mechanical loading in females during side cutting

BACKGROUND Side cutting involves mechanical loading of the knee which has been associated with anterior cruciate ligament injury risk. Despite a fast growing body of research, the relationship between loading mechanisms and running speed is still unclear. The aim of this study was to investigate how running speed determines a likely trade-off between task achievement and actual mechanical loading. METHODS Fourteen female participants (mean age=20.6±0.7yr, height=1.66±0.05m, mass=57.5±6.9kg) performed 45° side cutting manoeuvres at 2, 3, 4 and 5ms(-1) approach speeds. Three dimensional motion and ground reaction forces were recorded to calculate whole body centre of mass (CoM) velocity and lower limb kinematics and kinetics, focusing on knee flexion angle at touch-down and peak knee valgus loading during weight acceptance. One-way repeated measures ANOVA and one-dimensional statistical parametric mapping were used to identify significant speed effects on task achievement and mechanical loading. RESULTS Analysis of CoM velocities revealed that side cutting manoeuvres at higher running speeds matched the task requirements to a lesser extent. Despite a gradual increase of anterior-posterior deceleration and medio-lateral acceleration with running speed, knee loading mechanisms only reached meaningful levels from a 4ms(-1) running speed. CONCLUSION Our results confirmed a trade-off between task achievement and actual mechanical loading. This identified a need for standardisation of reporting running speeds. Taking into account also safety considerations, standardisation of a 4ms(-1) running speed is proposed for female athletes.

Registration of pedobarographic image data in the frequency domain

Image registration has been used to support pixel-level data analysis on pedobarographic image data sets. Some registration methods have focused on robustness and sacrificed speed, but a recent approach based on external contours offered both high computational processing speed and high accuracy. However, since contours can be influenced by local perturbations, we sought more global methods. Thus, we propose two new registration methods based on the Fourier transform, cross-correlation and phase correlation which offer high computational speed. We found out that both proposed methods revealed high accuracy for the similarity measures considered, using control geometric transformations. Additionally, both methods revealed high computational processing speed which, combined with their accuracy and robustness, allows their implementation in near-real-time applications. Furthermore, we found that the current methods were robust to moderate levels of noise, and consequently, do not require noise removal procedure like the contours method does.