Dynamic navicular motion measured using a stretch sensor is different between walking and running, and between over-ground and treadmill conditions
© Barton et al.; licensee BioMed Central. 2015
Received: 30 July 2014
Accepted: 6 February 2015
Published: 25 February 2015
Non-invasive evaluation of in-shoe foot motion has traditionally been difficult. Recently a novel ‘stretch-sensor’ was proposed as an easy and reliable method to measure dynamic foot (navicular) motion. Further validation of this method is needed to determine how different gait analysis protocols affect dynamic navicular motion.
Potential differences in magnitude and peak velocity of navicular motion using the ‘stretch sensor’ between (i) barefoot and shod conditions; (ii) overground and treadmill gait; and/or (iii) running and walking were evaluated in 26 healthy participants. Comparisons were made using paired t-tests.
Magnitude and velocity of navicular motion was not different between barefoot and shod walking on the treadmill. Compared to walking, velocity of navicular motion during running was 59% and 210% higher over-ground (p < 0.0001) and on a treadmill (p < 0.0001) respectively, and magnitude of navicular motion was 23% higher during over-ground running compared to over-ground walking (p = 0.02). Compared to over-ground, magnitude of navicular motion on a treadmill was 21% and 16% greater during walking (p = 0.0004) and running (p = 0003) respectively. Additionally, maximal velocity of navicular motion during treadmill walking was 48% less than walking over-ground (p < 0.0001).
The presence of footwear has minimal impact on navicular motion during walking. Differences in navicular motion between walking and running, and treadmill and over-ground gait highlight the importance of task specificity during gait analysis. Task specificity should be considered during design of future research trials and in clinical practice when measuring navicular motion.
KeywordsNavicular drop Navicular motion Foot kinematics In-shoe Stretch-sensor
Excessive or abnormal foot pronation mechanics are thought to result in lower extremity mal-alignment and pathology due to joint coupling with the tibia and femur [1,2]. Based on this theory, several clinical measures of foot pronation have been proposed, attempting to identify those at greatest risk of injury and to guide treatment decisions (e.g. foot orthoses prescription). The reason for multiple proposed measures of foot pronation may be due to the controversy related to which segments and planes of movement are associated. Importantly, there is also large variation in foot motion between individuals in relation to how specific joints move and interact with each other .
One common clinical method to measure foot pronation proposed by Brody , static navicular drop, involves calculation of the difference in vertical navicular height between sub-talar joint neutral and relaxed standing positions . Menz  proposed that the addition of navicular drift or “medial bulging” should be combined with navicular drop measurement. Additionally, he proposed these measures may offer significant benefits over more traditional measures of foot pronation such as rearfoot angle, providing an indication of talonavicular (midfoot) motion . Although, there is a paucity of research evaluating the link between navicular drift and lower limb injury, current evidence indicates greater navicular drop may be a risk factor for both medial tibial stress syndrome (MTSS)  and patellofemoral pain (PFP) . However, associated effect sizes are small [6,7], possibly owing to the multifactorial nature of these conditions. Additionally, static navicular drop is a poor predictor of dynamic foot function , taking into account only sagittal plane motion without consideration to other planes of movement .
Non-invasive measurement of in-shoe foot motion has been traditionally difficult, due to the need for skin-marker placement over the foot. One alternative is to place markers over the shoe, however, previous research indicates this method may poorly reflect true foot motion. Previous studies have reported both over- [9,10] and under-estimation  of actual foot motion when measured by markers on the shoe. Shoe windows have also been cut, allowing marker placement, but this affects the shoes structural integrity , and is not feasible in a clinical setting. More recently, Christensen et al.  proposed the use of a ‘stretch-sensor’ as an easy and efficient method to measure and provide insight into shod and barefoot tri-planar navicular motion. This may prove to be a valuable tool in identifying injury risk, guiding treatment decisions, and measuring the effectiveness of interventions designed to control navicular motion. However, further validation is needed to determine how different gait analysis protocols affect dynamic navicular motion.
Evaluation of running and walking biomechanics is frequently completed on a treadmill. Previously, most research has focused on sagittal plane differences and has measured predominantly proximal motion. Studies comparing treadmill to over-ground gait indicate decreased peak and range of knee flexion during both walking and running on a treadmill [14-19]; inconsistent differences for hip flexion during running with both increased [20,21] and decreased  peaks on a treadmill; decreased ankle dorsiflexion range of motion, velocity and peak [19,22] when running on a treadmill; and greater rearfoot/ankle eversion during running (1.5 – 6.30) [18,19,23] on a treadmill, although one study did report no significant difference in rearfoot motion . To the author’s knowledge, there does not appear to be any comparison of midfoot or navicular motion during treadmill and over-ground walking or running. Additionally, no previous research has compared navicular motion between barefoot and shod conditions, or between walking and running.
This study investigates whether differences in magnitude and peak velocity of navicular motion exist between (i) running and walking; (ii) overground and treadmill gait; and/or (iii) barefoot and shod conditions. It was hypothesised that magnitude and velocity of navicular motion would be greater during running compared to walking.
This cross-sectional study compared navicular motion during five different test conditions in 26 (11 female and 15 male) participants. Mean (standard deviation) for age, height, weight, and foot length of participants was 27 (7) years, 178 (9) cm, 69 (11) kg and 25.3 (1.9) cm respectively. The 26 participants (age 19-57) were conveniently sampled from the School of Engineering, Aarhus University, and surrounding community. Prospective participants with present injury or pain in the lower extremities or back preventing them from walking or running; or with medical or neurological conditions were excluded. The study was approved by Aarhus University and conducted in accordance with the Helsinki Declaration  and all participants were given written and verbal information about the project and signed an informed consent before participating.
Measure of dynamic navicular motion (using the stretch-sensor)
Data acquisition device
The data acquisition device was custom-made hardware that enabled recording of the elongation of a stretch sensor with a sample rate of 200Hz. The device contained a three axis gyroscope and a three axis accelerometer that was sampled simultaneously with the measurements stored on a SD card together with the data from the stretch sensor. The device was mounted just above the malleoli to minimise the movement of the device caused by contractions of the muscles in the lower leg.
All participants followed the same order of test conditions. Acclimatisation was determined by both visual observation by the investigator and participants answering ‘yes’ when asked “do you feel that you are walking or running as you normally would?” Test conditions in order were as follows:
Participant’s own standard athletic trainers and socks
Self- selected over-ground walking pace
Self-selected over-ground running pace
Treadmill walking (same pace as during over-ground walking)
Treadmill running (same pace as during over-ground running).
Barefoot (with socks)
Treadmill walking (same pace as during over-ground walking).
Primary outcomes included the magnitude and peak velocity of navicular motion measured by the stretch sensor. Velocity of navicular motion (i.e. peak velocity) was estimated by differentiating the elongation of the stretch-sensor, specifically convoluting the elongation of the stretch-sensor with the filter mask [-0.5 0 0.5]. The intra and intertester reliability of the sensor has previously been reported as acceptable with ICC(2.1) > 0.76 for barefoot measurements and ICC(2.1) 0.65 for shod measurements . Additionally, conservative estimates for 95% limits of agreement for intra- and inter-rater reliability have been established as between -2.4 and 2.6 mm . Testing was completed by the same tester, who had used the stretch-sensor approximately 100 times prior to study commencement.
The sample size was based on pilot data and aimed at detecting a 20% difference in navicular motion between shod overground walking and running. Using a standard deviation of 30%, power of 80% and alpha at 0.05, at least 20 participants were needed.
Separate one way analysis of variance (ANOVA) were completed for magnitude and velocity of navicular motion, with paired samples t-tests used to compare differences between test conditions. Bonferroni adjustment was not made for pairwise comparisons to ensure potentially clinically meaningful findings were not missed due to stringent statistical correction . Effect size (ES) were calculated using the mean difference between conditons divided by the standard deviation of the first condition . Additionally the percentage of differences were calculated for each comparison made. Following criteria proposed by Hume et al. , ES were allocated as small (<0.6), medium (0.61 – 1.19), or large (≥1.2). Pearson’s product moment correlation evaluated the association between the magnitude and velocity of navicular motion between test conditions.
Data of one participant was excluded from all trials due to poor data quality. Additionally, another participant was excluded from treadmill running, and four other participants were excluded from barefoot treadmill walking due to poor data quality. ANOVA results indicated a significant effect of condition on magnitude of navicular motion (F = 6.9, p = 0.002) and on velocity of navicular motion (F = 52.6, p = <0.001).
Barefoot versus shod treadmill walking
Navicular motion and velocity of navicular motion during over ground walking and running and during treadmill walking and running
TrW Vs TrBF
TrR Vs TrW
OGR Vs OGW
TrW Vs OGW
TrR Vs OGR
Difference (95% CI)
Difference (95% CI)
Difference (95% CI)
Difference (95% CI)
Difference (95% CI)
(n = 25, 72-134 steps)
(n = 24, 108-156 steps)
(n = 21, 80-156 steps)
(n = 25, 53-72 steps)
(n = 25, 64-117 steps)
Navicular Motion – mm (95%CI)
(-0.5 to 1.5)
(-0.1 to 2.4)
(0.3 to 2.9)
(0.8 to 2.3)
(0.5 to 2.3)
p = 0.34
p = 0.07
p = 0.02
p = 0.0004
p = 0.003
Velocity* - mm/s (95%CI)
(-8 to 44)
(216 to 330)
(117 to 179)
(86 to 152)
(-22 to 26)
p = 0.16
p < 0.0001
p < 0.0001
p < 0.0001
p = 0.85
Running versus walking
Compared to walking, velocity of navicular motion during running was 59% (ES = 1.03) and 210% (ES = 1.63) higher over-ground and on a treadmill respectively (see Table 1). Magnitude of navicular motion was 23% (ES = 0.64) higher during over-ground running compared to over-ground walking. Additionally, there was a statistical trend (p = 0.07) towards a 14% (ES = 0.31) increase in the magnitude of navicular motion during treadmill running compared to treadmill walking.
Over-ground versus treadmill walking and running
Compared to over-ground, magnitude of navicular motion on a treadmill was 21% (ES = 0.39) and 16% (ES = 0.36) greater during walking and running respectively (see Table 1). Compared to over-ground, maximal velocity of navicular motion was 48% (ES = 1.45) less during walking on a treadmill, but there were no differences in maximal velocity during running between the two conditions (see Table 1).
There were significant associations between magnitude (r = 0.88, p < 0.0001) and velocity of navicular motion (r = 0.66, p = 0.003) during over-ground and treadmill walking. Magnitude (r = 0.88, p < 0.0001) and velocity (r = 0.95, p < 0.0001) of navicular motion during over-ground running was strongly associated with treadmill running.
Traditionally, the need to place markers directly over the foot has made measuring in-shoe foot motion difficult. In this study, a viable method of measuring in-shoe navicular motion was used to identify possible differences between barefoot and shod walking, running and walking, and between treadmill and over-ground gait. Overall, these findings highlight the need for task specificity during gait analysis of navicular motion, with differences between running and walking, and treadmill and over-ground gait indicated.
Magnitude and velocity of navicular motion is greater during running compared to walking, regardless of whether comparison is made on a treadmill or over-ground. Specifically, small to moderate increases in magnitude, and moderate to large increases in velocity were found, indicating running is likely to have a greater effect on velocity than magnitude of navicular motion. Greater magnitude and velocity of navicular motion is likely the result of greater vertical ground reaction forces occurring during running [32,33]. Importantly, greater navicular motion during running highlights the need for task specificity when evaluating foot function in a research or clinical setting.
Evaluation of running and walking biomechanics is frequently completed on a treadmill. Previous research has indicated likely differences between treadmill and over-ground gait at the hip [19-21], knee [14-19], and ankle/rearfoot [19,22]. To the author’s knowledge, this is the first study to compare midfoot motion between treadmill and over-ground walking or running. Although a large reduction of navicular motion velocity occurred during treadmill walking, moderate to large increases in magnitude of navicular motion resulted during treadmill walking and running compared to over-ground. Importantly, similar to sagittal plane motion at the hip, knee and ankle, this indicates navicular motion on a treadmill may not accurately reflect over-ground conditions, particularly during walking. This should be considered during gait assessment in both research and clinical practice settings.
Greater navicular motion during treadmill gait in this study is consistent with increases in ankle eversion (1.5 – 6.30) previously reported during treadmill running [18,19,23]. This consistency between different studies is logical, since navicular and rearfoot motion may both be considered proxies for foot pronation , and strong correlations between the two have been previously reported . Reduced velocity of navicular motion during walking on a treadmill compared to over-ground is also consistent with previous reports of reduced internal inversion moments during treadmill walking , which is a measure of forces required to control foot pronation. Further research concurrently collecting stretch sensor data with ground reaction forces, muscle function and more proximal kinematics is needed to understand potential mechanisms for the differences in navicular motion between treadmill and overground gait.
Foot motion during gait is commonly measured barefoot due to issues with marker placement on the shoe . This method is often criticised as barefoot may not accurately reflect shod kinematics. Findings in this study indicate the magnitude and velocity of navicular motion was not significantly different between barefoot and shod treadmill walking. Additionally, there were strong correlations between the two conditions for both magnitude and velocity of navicular motion, indicating footwear may have minimal impact on navicular motion during walking. Whilst walking findings from this study align with previous bone-pin research indicating minimal difference in calcaneal and tibial movement patterns between barefoot and shod running , other research indicates barefoot running may be associated with reduced and earlier rearfoot eversion . Additionally, increased step rate, reduced ankle dorsiflexion, and lower vertical ground reaction forces have also been reported during barefoot running . Further research is needed to explore if differences in navicular motion exist between barefoot and shod running, and the potential implications of this.
Limitations and future research
Poor data quality for some conditions led to the exclusion of some participant data. Specifically, data for treadmill running was excluded in one participant, and data from four other participants were excluded from barefoot treadmill walking. Poor data quality in these instances is attributed to accumulation of sweat during testing, and subsequent loosening of the stretch sensor attachments, with the data lost coming from the final two conditions tested. Poor data quality was determined by very large navicular drop, e.g. 20 or 30mm, a velocity of more than 2000 mm/s, or a loss of the normal rhythm (see Figure 2). To address this in future research, the use of skin adhesive, evaluating fewer conditions, and/or reapplying and recalibrating the ‘stretch-sensor’ during testing is recommended.
Participants performed testing conditions in a set order to allow matching of outdoor walking and running speed. This may have led to systematic differences between conditions as result of fatigue , potentially limiting generalisability of findings. However, Dierks et al.  has previously reported greater excursion and peak velocity of rearfoot eversion during running following an exhaustive running protocol, which averaged 45 minutes and all 20 runners reached a rated perceived exertion of 15. Increases were 1.20 (15%) and 12.80 per second (11%) respectively. Considering the protocol used in the current study was far less fatiguing, and differences identified were of larger magnitudes (16 – 210%), it is unlikely that fatigue substantially contributed to current findings. Additionally, unpublished data indicates that navicular motion measured by the stretch sensor remains consistent until at least 30 minutes of fast barefoot walking. Nonetheless, further research is warranted to confirm findings from this study, and to understand the impact of fatigue on navicular motion measured by the stretch sensor during various conditions. Aclimatisation to treadmill running varied between two and five minutes, and not six minutes recommended by Matsas . However, all participants reported feeling comfortable in each testing condition prior to commencement of data collection. Further research is needed to determine an adequate aclimitisation period to ensure consistent navicular motion during treadmill gait.
Due to the anatomical placement of the stretch sensor over the navicular, it is thought to provide insight into tri-planar talonavicular (midfoot) motion during gait . However, further research comparing motion measured using the stretch sensor to three dimensional motion of the midfoot and its individual components including the talonavicular joint, as well as rearfoot and forefoot segments is needed to understand the possible biomechanical information it can provide.
Measurement of navicular motion using the ‘stretch-sensor’ is a reliable method of measuring in-shoe navicular motion during gait, and can detect differences between varied gait analysis protocols. Magnitude and velocity of navicular motion is higher during running compared to walking, magnitude of navicular motion is higher during walking and running on a treadmill compared to over-ground, and velocity of navicular motion is lower during walking on a treadmill compared to overground. These differences highlight the importance of task specificity during gait analysis, and should be considered during design of future research trials and during gait assessment in clinical practice.
- Powers CM. The influence of altered lower-extremity kinematics on patellofemoral joint dysfunction: a theoretical perspective. J Orthop Sports Phys Ther. 2003;33(11):639–46.View ArticlePubMedGoogle Scholar
- Tiberio D. The effect of excessive subtalar joint pronation on patellofemoral mechanics: a theoretical model. J Orthop Sports Phys Ther. 1987;9(4):160–5.View ArticlePubMedGoogle Scholar
- Nester CJ. Lessons from dynamic cadaver and invasive bone pin studies: do we know how the foot really moves during gait? J Foot Ankle Res. 2009;2:18.View ArticlePubMedPubMed CentralGoogle Scholar
- Brody DM. Techniques in the evaluation and treatment of the injured runner. Orthop Clin North Am. 1982;13(3):541–58.PubMedGoogle Scholar
- Menz HB. Alternative techniques for the clinical assessment of foot pronation. J Am Podiatr Med Assoc. 1998;88(3):119–29.View ArticlePubMedGoogle Scholar
- Newman P, Witchalls J, Waddington G, Adams R. Risk factors associated with medial tibial stress syndrome in runners: a systematic review and meta-analysis. Open Access J Sports Med. 2013;4:229–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Boling MC, Padua DA, Marshall SW, Guskiewicz K, Pyne S, Beutler A. A prospective investigation of biomechanical risk factors for patellofemoral pain syndrome: the Joint Undertaking to Monitor and Prevent ACL Injury (JUMP-ACL) cohort. Am J Sports Med. 2009;37(11):2108–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Rathleff MS, Nielsen RG, Kersting UG. Navicula drop test ad modum Brody: does it show how the foot moves under dynamic conditions? J Am Podiatr Med Assoc. 2012;102(1):34–8.View ArticlePubMedGoogle Scholar
- Reinschmidt C, van Den Bogert AJ, Murphy N, Lundberg A, Nigg BM. Tibiocalcaneal motion during running, measured with external and bone markers. Clin Biomech (Bristol, Avon). 1997;12(1):8–16.View ArticleGoogle Scholar
- Reinschmidt C, Stacoff A, Stüssi E. Heel movement within a court shoe. Med Sci Sports Exerc. 1992;24(12):1390–5.View ArticlePubMedGoogle Scholar
- Sinclair J, Greenhalgh A, Taylor P, Edmundson C, Brooks D, Hobbs S. Differences in tibiocalcaneal kinematics measured with skin- and shoe-mounted markers. Hum Mov. 2013;14(1):64–9.Google Scholar
- Arnold JB, Bishop C. Quantifying foot kinematics inside athletic footwear: a review. Footwear Sci. 2013;5(1):55–62.View ArticleGoogle Scholar
- Christensen BH, Andersen KS, Pedersen KS, Bengtsen BS, Simonsen O, Kappel SL, et al. Reliability and concurrent validity of a novel method allowing for in-shoe measurement of navicular drop. J Foot Ankle Res. 2014;7(1):12.View ArticlePubMedPubMed CentralGoogle Scholar
- Riley PO, Paolini G, Della Croce U, Paylo KW, Kerrigan DC. A kinematic and kinetic comparison of overground and treadmill walking in healthy subjects. Gait Posture. 2007;26(1):17–24.View ArticlePubMedGoogle Scholar
- Strathy GM, Chao EY, Laughman RK. Changes in knee function associated with treadmill ambulation. J Biomech. 1983;16(7):517–22.View ArticlePubMedGoogle Scholar
- Murray MP, Spurr GB, Sepic SB, Gardner GM, Mollinger LA. Treadmill vs. floor walking: kinematics, electromyogram, and heart rate. J Appl Physiol. 1985;59(1):87–91.PubMedGoogle Scholar
- Matsas A, Taylor N, McBurney H. Knee joint kinematics from familiarised treadmill walking can be generalised to overground walking in young unimpaired subjects. Gait Posture. 2000;11(1):46–53.View ArticlePubMedGoogle Scholar
- Riley PO, Dicharry J, Franz J, Della Croce U, Wilder RP, Kerrigan DC. A kinematics and kinetic comparison of overground and treadmill running. Med Sci Sports Exerc. 2008;40(6):1093–100.View ArticlePubMedGoogle Scholar
- Sinclair J, Richards J, Taylor P, Edmundson C, Brooks D, Hobbs S. Three-dimensional kinematic comparison of treadmill and overground running. Sports Biomech. 2013;12(3):272–82.View ArticlePubMedGoogle Scholar
- Alton F, Baldey L, Caplan S, Morrissey MC. A kinematic comparison of overground and treadmill walking. Clin Biomech (Bristol, Avon). 1998;13(6):434–40.View ArticleGoogle Scholar
- Wall JC, Charteris J. A kinematic study of long-term habituation to treadmill walking. Ergonomics. 1981;24(7):531–42.View ArticlePubMedGoogle Scholar
- Fellin RE, Manal K, Davis IS. Comparison of lower extremity kinematic curves during overground and treadmill running. J Appl Biomech. 2010;26(4):407–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Nigg BM, De Boer RW, Fisher V. A kinematic comparison of overground and treadmill running. Med Sci Sports Exerc. 1995;27(1):98–105.PubMedGoogle Scholar
- Association WM. Ethical principles for medical research involving human subjects. Eur J Emerg Med. 2001;8(3):221–3.View ArticleGoogle Scholar
- Kappel SL, Rathleff MS, Hermann D, Simonsen O, Karstoft H, Ahrendt P. A novel method for measuring in-shoe navicular drop during gait. Sensors. 2012;12(9):11697–711.View ArticlePubMedPubMed CentralGoogle Scholar
- Wolf P, Stacoff A, Liu A, Nester C, Arndt A, Lundberg A, et al. Functional units of the human foot. Gait Posture. 2008;28(3):434–41.View ArticlePubMedGoogle Scholar
- Rathleff MS, Nielsen RG, Simonsen O, Olesen CG, Kersting UG. Perspectives for clinical measures of dynamic foot function-reference data and methodological considerations. Gait Posture. 2010;31(2):191–6.View ArticlePubMedGoogle Scholar
- Kavanagh JJ, Menz HB. Accelerometry: a technique for quantifying movement patterns during walking. Gait Posture. 2008;28(1):1–15.View ArticlePubMedGoogle Scholar
- Perneger T. What's wrong with Bonferroni adjustments? Br Med J. 1998;316:1236–8.View ArticleGoogle Scholar
- Morris S. Estimating effect sizes from pretest-posttest-control group designs. Organ Res Methods. 2008;11(2):364–86.View ArticleGoogle Scholar
- Hume P, Nigg BM, Reinschmidt C, van den Bogert AJ, Lundberg A. Effectiveness of foot orthoses for treatment and prevention of lower limb injuries - A review. Sports Med. 2008;38(9):759–79.View ArticlePubMedGoogle Scholar
- Dugan SA, Bhat KP. Biomechanics and analysis of running gait. Phys Med Rehabil Clin N Am. 2005;16(3):603–21.View ArticlePubMedGoogle Scholar
- Nilsson J, Thorstensson A. Ground reaction forces at different speeds of human walking and running. Acta Physiol Scand. 1989;136(2):217–27.View ArticlePubMedGoogle Scholar
- Cornwall MW, McPoil TG. Relative movement of the navicular bone during normal walking. Foot Ankle Int. 1999;20(8):507–12.View ArticlePubMedGoogle Scholar
- Stacoff A, Nigg BM, Reinschmidt C, van den Bogert AJ, Lundberg A. Tibiocalcaneal kinematics of barefoot versus shod running. J Biomech. 2000;33(11):1387–95.View ArticlePubMedGoogle Scholar
- Morley JB, Decker LM, Dierks T, Blanke D, French JA, Stergiou N. Effects of varying amounts of pronation on the mediolateral ground reaction forces during barefoot versus shod running. J Appl Biomech. 2010;26(2):205–14.View ArticlePubMedGoogle Scholar
- Hall JPL, Barton C, Jones PR, Morrissey D. The biomechanical differences between barefoot and shod distance running: a systematic review and preliminary meta-analysis. Sports Med. 2013;43(12):1335–53.View ArticlePubMedGoogle Scholar
- Morin JB, Samozino P, Millet GY. Changes in running kinematics, kinetics, and spring-mass behavior over a 24-h run. Med Sci Sports Exerc. 2011;43(5):829–36.View ArticlePubMedGoogle Scholar
- Dierks TA, Manal KT, Hamill J, Davis I. Lower extremity kinematics in runners with patellofemoral pain during a prolonged run. Med Sci Sports Exerc. 2011;43(4):693–700.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.