Spinal lordosis optimizes the requirements for a stable erect posture
© Wagner et al; licensee BioMed Central Ltd. 2012
Received: 9 February 2012
Accepted: 16 April 2012
Published: 16 April 2012
Lordosis is the bending of the lumbar spine that gives the vertebral column of humans its characteristic ventrally convex curvature. Infants develop lordosis around the time when they acquire bipedal locomotion. Even macaques develop a lordosis when they are trained to walk bipedally. The aim of this study was to investigate why humans and some animals develop a lumbar lordosis while learning to walk bipedally.
We developed a musculoskeletal model of the lumbar spine, that includes an asymmetric, dorsally shifted location of the spinal column in the body, realistic moment arms, and physiological cross-sectional areas (PCSA) of the muscles as well as realistic force-length and force-velocity relationships. The model was used to analyze the stability of an upright body posture. According to our results, lordosis reduces the local joint torques necessary for an equilibrium of the vertebral column during an erect posture. At the same time lordosis increases the demands on the global muscles to provide stability.
We conclude that the development of a spinal lordosis is a compromise between the stability requirements of an erect posture and the necessity of torque equilibria at each spinal segment.
Keywordsmuscle physiology lordosis evolution spine stability biomechanics motor control Submitted to: Theoretical Biology and Medical Modelling
The lumbar region of the back is supported by short deep muscles, that connect the vertebrae, and long superficial muscles, that connect the thoracic cage and the pelvis. The first are usually referred to as local muscles. According to one view, the local muscles provide the stability of the vertebral column, whereas the superficial ones, the global muscles, would be the mobilizers [5–8]. This view has been challenged because it has been shown that global muscles also contribute to spinal stability [9–11]. Moreover, the local muscles are, in contrast to the global muscles, characterized by small lever arms and small cross-sectional areas, so that these muscles cannot generate large torques. This would be an undesirable property if the stability of the spine would depend only on the local muscles.
It has been shown that the spinal column of a standing human stores elastic energy [12, 13], but this elasticity cannot explain the efficiency of walking . Also, it has been suggested that lordosis in the lower back region might minimize the external moment of the centre of mass of the upper body, while retaining a stable hip joint position .
In the present work we start from the premise that stability control is central to the lumbar spine . Since the human lumbar spine is a loaded chain of joints (the intervertebral discs between the vertebrae), controlling its stability is inherently complex. Whereas the global muscles with their large moment arms are powerful enough, they can only control the chain, but not the individual joints. The local muscles might potentially stabilize each joint if they had the strength.
We thus hypothesize that the motor system selects a configuration in which the required local stabilizing torques on each of the lumbar joints is minimal.
We chose the concept of self-stability  to test this hypothesis. Self-stability is the stable performance of a musculoskeletal system without neuronal feedback. The reasoning underlying this approach is that neuronal control is time-delayed, and thus, if the musculoskeletal system is mechanically stable already, this enormously reduces the problem of stable control. The concept of self-stability relies fundamentally on the non-linear mechanical properties of muscles, as explained in the methods. For the model, we assume that every degree of freedom of each of the joints must be self-stable at any time, in order to maintain reliable physiological functioning of the spine.
The musculoskeletal model
The spinal column had a horizontal eccentricity E [m] (Figure 2). The eccentricity E is defined as the distance of the spinal column to the symmetrical axis. For example, in the sagittal plane, the eccentricity is equivalent to the dorsal location of the spine in the body. The spinal curvature was implemented as a cubic spline through P 1, P 3, and P 2. P 3 was located halfway between P 1 and P 2, whereas P 1 was located 176 mm above P 2 . As displayed in Figure 2, the spinal curvature parameter Λ [m] was defined as the horizontal position of P 3, with respect to the midline between P 1 and P 2. Both parameters were varied from 0 to 8 cm in 100 equidistant steps.
Three antagonistic pairs of global muscles were included, i.e. straight muscles acting in parallel to the spine and two oblique arrangements (Figure 2). Each global muscle consisted of five parallel muscle fibers, to simulate physiologically realistic surfaces of attachment of the muscles.
with optimal muscle length l opt = 1.2l0 and l0 muscle length at equilibrium. The specific constants k1 = 0.96, k2 = 0.35 were chosen such that the muscles were acting on the ascending limb of the force-length relationship [9, 19], which improves the self-stability .
The Hill-type muscle properties a [N], b [m/s] and c [W] can be derived from the physiological cross-sectional area PCS A [cm2], the proportion of fast-twitch-fibers FT , and the optimum muscle length l opt . The isometric force was estimated as f iso = 25 N/ cm2 ... PCS A cm2 and the maximum contraction velocity as v max = (6 + 10 ... FT) ... l opt , with FT = 0.5. From this the Hill-type muscle properties were calculated as: a = f iso / 4, b = v max / 4, c = b(f iso + a), A = f ecc ... f iso , B = (A − f iso )b2/c, C = (A − f iso )b, where the eccentric force enhancement f ecc = 1.5 . Finally, the torque generated by the i 'th muscle fibre was calculated as M i = h i ... f li ... f Hi . Where h i [m] is the lever arm of muscle i, with respect to the CoR.
with g gravitational acceleration and the point mass m with respect to the vertical axis, r the distance between the CoR and the point mass m (Figure 2B).
To guarantee an equilibrium at the CoR the counter torques CT local were calculated, representing the effect of local muscles, to counteract the sum of the torques generated by global muscles and the external torque due to the point mass m. The dynamic properties of the local muscles were neglected, because the torques are small compared with those of the global muscles.
Stability analysis and calculation of counter torques
The term . is always negative due to the negative slope of the force-velocity relationship . Hence, for a stable situation it is necessary that is negative, too. This was done for every CoR, resulting in five Jacobian matrices, each with two eigenvalues. The minimum PCS A needed to sustain stability was calculated numerically as the sum of the PCS A s of the muscle fibers. In order to find this minimum PCS A, the PCS A s of the three antagonistic muscles were varied iteratively.
Furthermore, it was necessary to derive the local counter torques CT k for each CoR to guarantee for an equilibrium. Finally, the Euclidian norm ||CT|| of the vector CT = (CT1, CT2, CT3, CT4, CT5) was calculated. The eccentricity E and the spinal curvature Λ were varied between 0 cm up to 8 cm in equidistant steps of 1 mm, respectively. From this a contour plot was drawn to analyze the sensitivity of the minimum PCS A, as well as the minimum counter torque in dependence of the eccentricity and spinal curvature.
The human spine is densely packed in an elastic, ligamentous sheath. We assume the spinal joints as linear elastic springs. In resting position the spring constant has been measured to be about K = (180π/π) ... 40° ... 10 = 15 N m/rad in the lumbar region of the spine .
for the eigenvalues λ1,2.
Note, that this threefold reduction in minimum PCS A with an increasing eccentricity is remarkable. The force-velocity Hill properties of the global muscles make one eigenvalue of the Jacobian always negative. Therefore, the other eigenvalue of the Jacobian causes this reduction in PCS A. Since the derivatives of the muscular torques M i with respect to α and ω are highly nonlinear, the decrease of the minimum PCS A is nonlinear too.
The analysis has so far ignored the contribution of elastic elements to the spinal system. Since the ligamentous sheath tends to erect the spine, adding elasticity tends to stabilize the spinal system, as can be seen from the negative sign of K in Equation 2. As a result, the required minimal PCSA of the global musculature is reduced. Apart from this quantitative reduction, notice that the addition of elastic components does not affect the relationship for the required PCSA qualitatively (see Additional file 1).
The same holds true for the local counter torques. The addition of elastic components does not change the range of optimal Λ and E.
The present study explored the optimality landscape of a spinal curvature by means of a simple model for mechanical stability. In the landscape where the lumbar spinal system is self-stable there is a range where the requirements to the local muscles are within the physiological constraints. The physiologically feasible spinal curvature depends on the eccentricity of the spine in the body. The optimal lumbar lordosis corresponds with literature findings, as we discuss below.
There are two different fundamental objectives for building a model. One objective is to build a model as detailed, complete and realistic as possible. It is evident that such a model would not only be highly complex and numerically expensive, it also would have not much explanatory power, for the underlying causality of the phenomena would still remain inaccessible. The second and more promising objective, therefore, is to build a model that is maximally simple, i.e. a model that contains as few conceptual elements as necessary to still explain those phenomena that are in need of explanation . As matters are, a model requires a particularly delicate balance between physiological detail and conceptual abstraction, hence, it is hard to draw direct conclusions to a living person. We deliberately chose to design the model as simple as possible, in order to acquire maximal insight into the general principles. The most dramatic simplification may be the exclusion of dynamics. We employed a stability analysis of an upright body posture based on the eigenvalue theory of Lyapunov . This approach treats a special case from the infinite range of possible dynamical states. The stability analysis thus estimates a property that is valid for dynamical systems in general. This special case therefore represents a minimal requirement.
The results indicate that the sagittal and the coronal plane each represent a solution on a continuous range (i.e., the presence of lumbar lordosis and the absence of scoliosis). Therefore, the expansion to a three-dimensional model would not change the general conclusion as movements in intermediate planes simply represent the intermediate results of the planar model.
As mentioned above, the muscular arrangement acting around the lumbar spine has been classified into global mobilizing muscles, and local stabilizing ones [5, 6]. The approach of our model is different. Instead of assigning different supposed functions to different muscle groups, we derived the model on the presumption that stability is a central requirement for a proper functioning of the spinal system. Following this view, it has been found in earlier studies that self-stability of the spine requires an increased coactivation of antagonistic muscles . Our simulations strongly suggest that self-stability is indeed critical in the spinal human system. The self-stable area in Figure 4 is within the physiological domain (Position C), but it is small. Therefore, the mechanical stability of the spine may be an important driving force in the ontogeny and evolution of humans.
The simulations of the present model require an eccentricity of E ≈ 3 cm and a lumbar lordosis of Λ ≈ 4 cm. Both values are within the physiological range (see Figure 4, position C). For a symmetrical arrangement of the musculoskeletal system, such as the coronal plane of the lumbar trunk (Figure 4, position A), no local counter torques are necessary to generate an equilibrium. For the sagittal plane the musculoskeletal arrangement becomes asymmetrical, i.e., the spinal column is dorsally shifted (Figure 4, position B). Here, to guarantee an equilibrium, local counter torques are necessary which exceed ±70 Nm. Therefore, for a dorsally shifted, but straight lumbar spinal column, extremely strong local muscles would be required.
A realistic estimate of the maximal counter torques of the local muscles is about ±30 Nm. This limit is drawn as a contour in Figure 4. Assuming as a conservative estimate, a lever arm of 7.0 cm, this limit is equivalent to a PCS A of 17 cm2. For a muscle with a circular cross section, the corresponding muscular diameter is about 4.5 cm. The lever arm is the sum of the radius of the local muscles and the radius of the vertebral bodies, which will certainly not be more than 7 cm. Thus, it is clear that all arrangements outside the defined tolerance area are unrealistic for physiological reasons.
The model reveals a remarkable phenomenon. The range of E and Λ where the norm of the local torques is physiological, coincides with the range where the highest PCS A is required (cf. Figures 3 and 4). This is not a problem, because the required PCS A is well within the physiological range of the global muscles. The reason why the PCS A is high with small eccentricities is that the global muscles on both sides of the spine need to co-contract to obtain self-stability . For eccentric configurations, the global muscles with the small lever arm are already contracted in order to maintain the body in an upright position. Therefore, a co-contraction is not necessary for self-stability and thus, the required PCS A is much reduced.
That the required local torques are minimal in a similar range where the global PCS A shows the highest values is caused by a different mechanism. The geometrical range where small local torques are required is where the joints between the spinal segments are close to the gravity vector passing from the body mass m (cf. the dashed line in Figure 2A).
It should be mentioned here that additional stabilizing structures, such as the ligamentous sheath, ligaments and the intervertebral discs, will have a positive influence on the stability. This was confirmed by the simulations where we introduced a stiffness parameter K (eqn. 2). As a result, the required PCSA is halved (see Additional file 1).
The local muscles in the lumbar region are located only on the dorsal side of the spinal column. This asymmetric configuration does not contradict with the model because these muscles are in part antagonized by the elastic elements.
Reflexes may potentially also stabilize the back. However, reflexes are delayed by up to 100 ms following a perturbation, and due to the electromechanical processes the force generation of the muscles takes another 50-100 ms. Therefore, maintaining stability with time-delayed reflexive control is challenging, especially under loaded conditions . Reflex loops will increase the stability margin and thus reduce the required co-contraction of the local muscles. If the spinal column is already a self-stable system, the reflex delays are much less problematic for stable control. On the contrary, the combination self-stability and the reflex systems should provide stable spinal control in the limited dynamic range.
An important finding is that the required local torques are almost independent of the eccentricity parameter, E, as long as an optimal lordosis is adopted. Thus, the medial position of the spine in the coronal plane provides body symmetry. On the other hand, the eccentric position of the spine in the sagittal plane allows a dynamic adaptation of the lordosis to the current weight distribution. The cost of this configuration is that the compressive load on the spine is much higher than it would be in a symmetrical configuration.
A low tonus or weakness of the local muscles will move the spinal system away from the self-stable range. As a result, the lower back will develop a hypo- or hyperlordosis in the sagittal plane, or a scoliosis in the coronal plane. The latter is shown in Figure 4 Position D. Since the spinal system is no longer self-stable, the model predicts that a scoliosis will develop rapidly and is difficult to reverse, as is indeed the case.
Bipedalism and lumbar lordosis are often regarded as a cause of low back pain, through the increased spinal shearing forces and increased risk of spondylolisthesis [4, 8, 16, 29, 30]. In this respect it is of utmost importance to understand the function of lumbar lordosis. The human anatomy has a dorsally shifted spinal column. This anatomical arrangement causes varying local torques at the spinal segments. These local torques have to be counteracted by local muscles. Because the stability requirement requires a considerable coactivation of the muscles, the local torques for an equilibrium increase. The introduction of a lumbar lordosis reduces these local torques dramatically (Figure 4).
A convex curvature, analogous to lumbar lordosis can be found in the necks of many animals, especially long-necked birds and dinosaurs. In the light of the model, these spinal curvatures might also be optimal adaptations to the eccentric location of the head on the neck and the eccentric location of the cervical vertebrae in the neck.
As a conclusion it was possible to self-stabilize the spine in every single CoR with a single activation pattern of global muscles. Our simulations support the hypothesis, that lumbar lordosis is a mechanical consequence of the dorsally shifted spine in the sagittal plane to minimize the local counter torques at the spinal segments while maintaining self-stability. On the other hand, due to the symmetrical arrangement of the muscles in the coronal plane, a scoliosis represents a pathological disorder (see Figure 4, position D).
More research is needed to explore the proposed relationship in living subjects.
This work was supported by the Federal Ministry of Education and Research BMBF [grant 01EC1003A].
This work was supported by the Federal Ministry of Education and Research BMBF [grant 01EC1003A]. We thank the anonymous reviewers for their constructive and helpful comments.
- Preuschoft H, Hayama S, Gunther MM: Curvature of the lumbar spine as a consequence of mechanical necessities in Japanese macaques trained for bipedalism. Folia Primatol. 1988, 50 (1-2): 42-58. 10.1159/000156333.View ArticlePubMedGoogle Scholar
- Robinson JT: Early hominid posture and locomotion. 1973, Chicago: University of Chicago PressGoogle Scholar
- Schilling N, Arnold D, Wagner H, Fischer MS: Evolutionary aspects and muscular properties of the trunk-Implications for human low back pain. Pathophys. 2005, 12 (4): 233-242. 10.1016/j.pathophys.2005.09.005.View ArticleGoogle Scholar
- Whitcome KK, Shapiro LJ, Lieberman DE: Fetal load and the evolution of lumbar lordosis in bipedal hominins. Nature. 2007, 450 (7172): 1075-1078. 10.1038/nature06342.View ArticlePubMedGoogle Scholar
- Bergmark A: Stability of the lumbar spine. Acta Orthop Scandinav. 1989, 230S: 1-54.View ArticleGoogle Scholar
- Comerford MJ, Mottram SL: Movement and stability dysfunction-Contemporary developments. Man Ther. 2001, 6: 15-26. 10.1054/math.2000.0388.View ArticlePubMedGoogle Scholar
- Gibbons SGT, Comerford MJ: Strength versus stability. Part 1: Concepts and terms. Orthop Division Rev. 2001, 2: 21-27.Google Scholar
- Panjabi MM: A hypothesis of chronic back pain: Ligament subfailure injuries lead to muscle control dysfunction. Eur Spine J. 2006, 15 (5): 668-676. 10.1007/s00586-005-0925-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Cholewicki J, McGill SM: Mechanical stability of the in vivo lumbar spine: Implications for injury and chronic low back pain. Clin Biomech. 1996, 11 (2): 1-15.View ArticleGoogle Scholar
- Kavcic N, Grenier S, McGill SM: Determining the stabilizing role of individual torso muscles during rehabilitation exercises. Spine. 2004, 29 (11): 1254-1265. 10.1097/00007632-200406010-00016.View ArticlePubMedGoogle Scholar
- Wagner H, Anders C, Puta C, Petrovitch A, Morl F, Schilling N, Witte H, Blickhan R: Musculoskeletal support of lumbar spine stability. Pathophys. 2005, 12 (4): 257-265. 10.1016/j.pathophys.2005.09.007.View ArticleGoogle Scholar
- Preuschoft H: Mechanisms for the acquisition of habitual bipedality: Are there biomechanical reasons for the acquisition of upright bipedal posture?. J Anat. 2004, 204 (5): 363-384. 10.1111/j.0021-8782.2004.00303.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Witte H, Recknagel S, Rao JG, Wüthrich M, Lesch C: Is elastic energy storage of quantitative relevance for the functional morphology of the human locomotor apparatus?. Acta Anat. 1997, 158 (2): 106-111. 10.1159/000147919.View ArticlePubMedGoogle Scholar
- Alexander RM: Human energetics. Making headway in Africa. Nature. 1986, 319 (6055): 623-624.View ArticleGoogle Scholar
- Wagner H, Blickhan R: Stabilizing function of skeletal muscles: An analytical investigation. J Theor Biol. 1999, 199 (2): 163-179. 10.1006/jtbi.1999.0949.View ArticlePubMedGoogle Scholar
- Bogduk N: Clinical anatomy of the lumbar spine and sacrum. 1997, Edinburgh: Churchill Livingstone, 3Google Scholar
- Hatze H: Estimation of myodynamic parameter values from observations on isometrically contracting muscle groups. Eur J Appl Physiol. 1981, 46 (4): 325-338. 10.1007/BF00422120.View ArticleGoogle Scholar
- Otten E: Optimal design of vertebrate and insect sarcomeres. J Morphol. 1987, 191: 49-62. 10.1002/jmor.1051910106.View ArticlePubMedGoogle Scholar
- Cholewicki J, McGill SM: EMG assisted optimization: A hybrid approach for estimating muscle forces in an indeterminate biomechanical model. J Biomech. 1994, 27 (10): 1287-1289. 10.1016/0021-9290(94)90282-8.View ArticlePubMedGoogle Scholar
- Wagner H, Blickhan R: Stabilizing function of antagonistic neuromusculoskeletal systems: An analytical investigation. Biol Cybern. 2003, 89: 71-79.PubMedGoogle Scholar
- Hill AV: The heat of shortening and the dynamic constants of muscle. Proc R Soc Lond B. 1938, 126: 136-195. 10.1098/rspb.1938.0050.View ArticleGoogle Scholar
- Siebert T, Sust M, Thaller S, Tilp M, Wagner H: An improved method to determine neuromuscular properties using force laws-From single muscle to applications in human movements. Hum Mov Sci. 2007, 26 (2): 320-341. 10.1016/j.humov.2007.01.006.View ArticlePubMedGoogle Scholar
- Thaller S, Wagner H: The relation between Hill's equation and individual muscle properties. J Theor Biol. 2004, 231 (3): 319-332. 10.1016/j.jtbi.2004.06.027.View ArticlePubMedGoogle Scholar
- Till O, Siebert T, Rode C, Blickhan R: Characterization of isovelocity extension of activated muscle: A Hill-type model for eccentric contractions and a method for parameter determination. J Theor Biol. 2008, 255 (2): 176-187. 10.1016/j.jtbi.2008.08.009.View ArticlePubMedGoogle Scholar
- Lyapunov AM: Problèeme général´ de la stabilité du mouvement. (reprint 1949). Ann Fac Sci Toulouse. 1892, 2 (9): 203-474.Google Scholar
- Arjmand N, Shirazi-Adl A: Model and in vivo studies on human trunk load partitioning and stability in isometric forward flexions. J Biomech. 2006, 39 (3): 510-521. 10.1016/j.jbiomech.2004.11.030.View ArticlePubMedGoogle Scholar
- Full RJ, Koditschek DE: Templates and anchors: neuromechanical hypotheses of legged locomotion on land. J Exp Biol. 1999, 202 (Pt 23): 3325-32.PubMedGoogle Scholar
- Franklin TC, Granata KP: Role of reflex gain and reflex delay in spinal stability-A dynamic simulation. J Biomech. 2007, 40 (8): 1762-1767. 10.1016/j.jbiomech.2006.08.007.View ArticlePubMedGoogle Scholar
- Östgaard H, Andersson G, Schultz A, Miller J: Influence of some biomechanical factors on low-back pain in pregnancy. Spine. 1993, 18: 61-65. 10.1097/00007632-199301000-00010.View ArticlePubMedGoogle Scholar
- Liebetrau A, Puta C, Schinowski D, Wulf T, Wagner H: Is there a correlation between back pain and stability of the lumbar spine in pregnancy? A model-based hypothesis. Der Schmerz. 2012, 26: 36-45. 10.1007/s00482-011-1125-1.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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.