- Research article
- Open Access
The effect of various weight-bearing activities on the motion of lumbar facet joints in vivo
Journal of Orthopaedic Surgery and Research volume 17, Article number: 114 (2022)
Lumbar facet joints (LFJs) are usually related to the pathogenesis of the spine. The purpose of this paper is to study the effects of lifting load on the motion of lower lumbar facet joints in vivo.
Ten healthy volunteers aged 25 ≤ 39 years, 5 males and 5 females, were recruited. Using a dual fluoroscopy imaging system (DFIS) combined with CT, firstly, the L3-S1 segment image scanned by CT was converted into a three-dimensional model. Then, the lumbar motion images of L3-S1 vertebrae taken by the DFIS under different loads (0 kg, 5 kg, 10 kg) and different body postures (maximum flexion and extension, maximum left and right bending, and maximum left and right torsion) were captured. Finally, in the Rhino software, the instantaneous motion state of the lumbar spine is reproduced by translation and rotation according to the anatomical structure of the lumbar spine and the previous images. With the help of computer software, a Cartesian coordinate system was placed in the center of each articular surface to measure the kinematics of the articular process and to obtain 6DOF data under different loads (0 kg, 5 kg, 10 kg) in the lumbar facet joints.
In the flexion and extension of the trunk, weight bearing reduced the translational range in the mid-lateral direction. In the L3/4 segment, the lateral translational range of the left and right facet joints gradually decreased with increasing load, and the translational range at 0 kg was significantly greater than that at 10 kg (left side: 0 kg, 0.86° ± 0.57°, 10 kg, 0.24° ± 0.26°, p = 0.01; right side: 0 kg, 0.86° ± 0.59°, 10 kg, 0.26° ± 0.27°, p = 0.01). In the L5/S1 segment, the translation range of the LFJ at 0 kg was significantly greater than that at 10 kg (p = 0.02). Other bending and rotation movements were not found to cause differential changes in the 6DOF of the LFJ. In bending, the rotation range was the largest in the L3/4 segment (p < 0.05) and gradually decreased from top to bottom. At the same level, there were significant differences in the translation range of the left and right facets in the anterior posterior and craniocaudal directions (p < 0.05).
Increasing the load has a significant impact on the coupled translational movement of lumbar facet joints. The asymmetric translational movement of the left and right facet joints may be a factor that accelerates the degeneration of facet joints.
Facet joints play an important role in the stability and movement of the entire spine . In addition to restricting the movement of the vertebrae, the main function of joints can also bear different types of loads, such as compression, stretching, shearing, and torsion [2, 3]. Lumbar facet joint osteoarthritis is also considered the main cause of chronic low back pain in 15–45% of patients [4,5,6]. Heavy lifting (over 25 kg) and trunk flexion and rotation are considered to be moderate risk factors for low back pain [7, 8].
From a clinical point of view, an accurate description of the kinematics of normal lumbar facet joints can help us to improve our understanding, diagnosis and treatment of degenerative spine diseases. A literature review shows that data reports on the motion patterns of lumbar facet joints are mainly based on cadaver models , animal models [9, 10], computed tomography (CT) [11,12,13] and magnetic resonance imaging (MRI)  had been gradually applied to describe the morphology of the lumbar spine and the classification of osteoarthritis. The dual fluoroscopy imaging system DFIS is one of the latest technologies used to study the range and direction of motion of small joints in vivo.
This technology had been widely used in various bone, spine and joint studies. There are no restrictions on the subjects' activities, which overcomes the limitations of MRI/CT in space, time and body position, and can obtain the in vivo kinematics data of different weight-bearing and activity states under the real physiological conditions of the human body. In terms of LFJs, Li and his colleagues  were the first to use this technique to quantify the in vivo kinematics of lumbar facets in nonfunctional weight-bearing positions. Yin  also studied the effect of degenerative facets on the range of motion of the lumbar spine through this system. Although these studies have greatly analyzed the kinematics of LFJs in vivo under functional loading, their research direction is limited to either non-weight bearing or flexion and extension. Because the movement of the spine is complex, most daily life and work require complex movements such as flexion, lateral bending and rotation under appropriate weight bearing. Thus, it is very important to understand the in vivo kinematics of LFJs under different postures and various loads. These data are of great significance to study the mechanism of delaying the disease progression of the lower lumbar spine.
Therefore, the purpose of this paper is to study the effects of lifting load (0 kg, 5 kg, 10 kg) on the flexion and extension, left–right bending and left–right rotation of lower lumbar facet joints. We used the 3D modeling method of a dual fluoroscopy imaging system (DFIS) combined with CT to measure the facet motion of the lumbar L3-S1 vertebral body. This technique has been validated by our team for accuracy and repeatability and can be used to noninvasively measure spinal motion during weight-bearing functional activities . We hypothesize that the in vivo kinematics of lumbar facets are affected by load materials and methods.
In this study, 10 healthy asymptomatic volunteers, 5 males and 5 females aged between 25 and 39 years, were recruited. The inclusion criteria were (1) normal spine development and lumbar mobility within the normal range and (2) no previous symptoms of low back pain and history of lumbar trauma and sprain. (3) BMI was within the normal range, 18.5 < BMI < 25, and (4) bone mineral density was normal. The exclusion criteria were (1) previous symptoms of low back pain, history of low back trauma and sprain, (2) spinal deformities, (3) currently pregnant, (4) BMI > 25 or BMI < 18.5, and (5) Pfirrmann disk degeneration classification > grade II. This study was approved by the medical ethics committee of Tianjin Hospital.
Original CT image acquisition of lumbar spine
All subjects were scanned in a supine and relaxed position using computed tomography (CT) (Sensation 16 Siemens, Germany). The subjects were instructed to avoid any aggressive activities within 2 h before the experiment and rest on the examination table for approximately 20 min before CT scanning. The high-resolution axial plane image was scanned to obtain a parallel digital image with a thickness of 0.625 mm and a resolution of 512 * 512 pixels, and the CT image was saved in medical digital imaging and communication format (DICOM). The DICOM images were imported into solid modeling software (mimics version 19.0), and the 3D anatomical vertebral body model of L3-S1 was constructed using the established and validated protocol (Fig. 1a, b) . After outputting the 3D model in binary stereolithography (STL) format, it was imported into solid modeling software (rhinoceros® Robert McNeel & Associates, Seattle, Washington). The vertebral bodies were organized into a group to establish a 3D model of the L3-S1 vertebral body (Fig. 1c).
Dual fluoroscopy imaging system
The dual fluoroscopy imaging system consists of two "C" arm x-ray machines of the same model. Two fixed "C" arm x-ray machines were placed in the vertical position (Fig. 2a. F1, F2) to capture lumbar images in different postures from two orthogonal directions at the same time. Under the guidance of two spinal surgeons, the subject's lumbar spine was observed to ensure that it is always within the cross-projection range of the two "C" arms during movement and that the lumbar L3-S1 segment is basically within the field of vision of two fluoroscopes. Each subject was required to perform maximum flexion to maximum extension, maximum left to right bending, and maximum left to right rotation at loads of 0 kg, 5 kg (carrying 5 kg sandbags) and 10 kg (carrying 5 kg sandbags at the front and rear). The sandbag is located on the chest and back to avoid restricting waist movement (Fig. 2b, c). Each image acquisition position was maintained for approximately 2 s to obtain an instantaneous x-ray double oblique perspective image. During this period, the researchers assisted in restraining the hip and knee joints to reduce errors caused by other joint activities, and the subjects wore custom lead clothes to protect their thyroids and gonads (Fig. 2a).
Three-dimensional reconstruction of lumbar facet joint motion
The in vivo position of the vertebral body in different weight-bearing positions can be reproduced in the modeling software with DFIS . The three-dimensional model obtained in previously was imported into the modeling software (Rhino), then we outlined the anatomical characteristics of the lumbar spine in each posture, which must to be outlined the significant structures such as vertebral body, superior and inferior articular process and spinous process in detail (Fig. 3, F1a, F2a). In Rhino, through vertebral fluoroscopy the translation and rotation of the three-dimensional vertebral model, the spatial position of the lumbar spine was adjusted so that the three-dimensional lumbar model could match two mutually perpendicular double-oblique x-ray fluoroscopy images at the same time to realize the matching of two-dimensional and three-dimensional images (Fig. 3, F1b, F2b). The above can ensure the accuracy of matching results to a great extent to reduce the penetration of the model. Another important way to reduce model penetration and ensure data scientificity is that this series of test processes require system calibration. According to the C-arm used in the research, we had specially made the calibration system instrument for DFIS and compiled the corresponding correction processing software program . The repeatability of this technique in reproducing the 6DOF kinematics of the human spine in vivo is less than 0.3 mm in translation and less than 0.7 degrees in rotation [17, 18]. This technique can reproduce the position of the vertebral body in different postures. By matching the L3-S1 standard model, we reproduced the in vivo motion of lumbar facet joints under different loads (0 kg, 5 kg, 10 kg) and different postures.
Calculation of motion range of lumbar facet joint
In the L3-S1 model obtained in the supine position, a right-hand Cartesian coordinate system was established in the center of the vertebral body (Fig. 4 a). The articular surface of the articular process of the target vertebral body were divided into nine parts. The center was designated as the volume center of the facet capsule (Fig. 4 b1). The facet capsule coordinate system can be established by moving the coordinate system to the center point (Fig. 4 b2), which could calculate the relative activity data of two adjacent joints (Fig. 4c). The X axis (red) was defined as perpendicular to the sagittal plane of the cone and pointing to the left, the Y axis (green) was defined as parallel to the sagittal plane of the cone and pointing to the back, and the Z axis (blue) was defined as perpendicular to the coronal plane of the cone and pointing to the left cephalic side.
The ROM size of LFJ was divided into rotation angle and space vector size. The rotation angles along the X, Y and Z axes were α, β and γ, which represents the rotation range of the vertebral body in the movement process. The X (red), y (green) and Z (blue) axes in the spatial coordinates represent the direction and the vector size in space, respectively (Fig. 4a). The angle and translation of lumbar facet joint movement were calculated by the range of motion difference method. The final difference is taken as the absolute value to observe the range of motion in the facet joint.
Two-way repeated-measures ANOVA was used to compare the range of motion of the lumbar facet joints at the L3/4, L4/5 and L5/S1 vertebral levels under weight-bearing conditions. Kinematics was the dependent variable, and weight bearing and vertebral body level were the independent variables. When a statistically significant difference was detected, a post hoc Newman–Keuls test was performed, and the level of significance was again chosen at p < 0.05. A paired t test was used to compare the difference in bilateral facet joint displacement. If p < 0.05, it was statistically significant. All analyses were performed using SPSS software (SPSS 19.0; IBM, Armonk, NY, USA), with continuous variables expressed as X ± S.
In the flexion and extension of the trunk, the lumbar facet joint mainly rotated around the medial and lateral axes. The average rotation ranges of L3/4, L4/5 and L5/S1 levels under load were 2.72° ± 2.45° to 3.58° ± 2.73°, 2.58° ± 1.81° to 3.59° ± 2.68°, and 2.16° ± 1.63° to 2.72° ± 2.45°, respectively. The load condition had no significant effect on the rotation range of the main rotating shaft and the coupling shaft (p > 0.05), and the data showed that there was no significant difference in the rotation range of the main rotating shaft in different horizontal segments (Fig. 5a). The translation direction was mainly craniocaudal, ranging from 1–2.5 mm. The load had no effect on it. Interestingly, for the coupled translation, we found that the load reduced the translation in the mid-lateral direction. In the L3/4 segment, the translation range of the left and right facets under 0 kg load was greater than under 5 kg and 10 kg and significantly greater than that under a 10 kg load (left: 0 kg, 0.86 ± 0.57°, 10 kg, 0.24 ± 0.26°, p = 0.01; right: 0.86° ± 0.59°, 0.26° ± 0.27°, p = 0.01). In the L5/S1 segment, the 0 kg load translation range of the left facet was significantly greater than the 10 kg load translation range (p = 0.02). In addition, when the left and right facets at the same level were loaded with 0 kg, the translation range of the L4/5 segment in the anteroposterior direction was greater on the left than on the right (p < 0.05). When loading 10 kg, the translation range of the L3/4 segment in the craniocaudal direction was greater on the left than on the right (p < 0.05) (Table 1).
In the lateral bending movement, the facet joint mainly rotated around the anterior and posterior axes. The average rotation ranges of L3/4, L4/5 and L5/S1 levels under load were 3.2° ± 2.03° to 4.93° ± 4.35°, 1.75° ± 1.24° to 2.6° ± 1.72°, and 0.9° ± 0.85° to 1.75° ± 0.75°, respectively. The data showed that when loading, the rotation range of the L3/4 segment was the largest (p < 0.05), and the rotation range gradually decreased from top to bottom (Fig. 5b). The main direction of translation was the craniocaudal direction, which was less than 2.5 mm. We did not find that weight bearing led to significant changes in the kinematics of lumbar facet joints (p > 0.05). In addition, for the left and right facets at the same level, the translation range of L3/4 in the craniocaudal direction and anteroposterior direction at 5 kg was greater on the left than on the right (p < 0.05). At a 10 kg load, the L5/S1 level was significantly lower on the left than on the right in anteroposterior direction (p < 0.05) (Table 1).
In terms of rotational motion, the lumbar facet joint has a more complex motion dominated by coupling motion. The average rotation ranges of L3/4, L4/5 and L5/S1 around the primary rotation axis craniocaudal axis were 1.05° ± 0.8° to 1.86° ± 2.28°, 1.03° ± 0.74° to 1.62° ± 2.32°, and 0.83° ± 0.67° to 1.17° ± 1.11°, respectively. We did not find that weight bearing led to significant changes in the kinematics of lumbar facets (p > 0.05), and the data showed that there was no significant difference in the rotation range of the main rotation axis at different horizontal segments (Fig. 5c). In addition, when the left and right facets at the same level were loaded with 5 kg, the left side of the translation range in the front and rear directions of L3/4 was less than that on the right side (left: 0.48 mm ± 0.3, right: 0.81 mm ± 0.45, p < 0.05) (Table 1).
In this study, the effects of lifting load (0 kg, 5 kg, 10 kg) on the flexion and extension, left–right bending and left–right rotation of lower lumbar facet joints were measured. In general, the lifting load can affect the coupling displacement of small joints. The data showed that during the flexion and extension exercises, it was found in the mid-lateral direction that the translation range of the L3/4 and L5/S1 segments at 0 kg weight was significantly greater than that at 10 kg (p < 0.05).
Overall, the results seem to be the same as those of previous studies. In Chowdhury's  study, it was found that the lifting load had a significant impact on LFJ movement, which was reflected in increasing the rotation range and up-down translation range of the main rotating shaft during flexion and extension movement, especially when the lifting load was 13.5 kg. In our study, it was found that the rotation range of the L3/4 and L4/5 levels was also the largest when the load was increased to the maximum 10 kg, but there was no significant difference with other loads. Interestingly, the data in this paper showed that lifting the load could significantly reduce the translation of the external axis. This suggests that physiological load may limit joint activity and discourage excessive activity, which may lead to interfacet impact, thus increasing the risk of facet degeneration. The absence of differences in other kinematics does not mean that weight bearing has no effect on joint kinematics, which may be because weight bearing is not common enough to affect the daily activities of normal adults. However, the literature shows that if the load is too large, such as more than 25 kg, it may become a medium risk factor for low back pain .
We found that the horizontal motion patterns of the upper and lower lumbar spines around the main rotation axis are inconsistent in lateral bending. Under 0 kg load and 10 kg load, the rotation range showed a large range at the L3/4 level (p < 0.05), and the rotation range gradually decreased from top to bottom. When we tested the subjects' in vivo movements, the pelvis and knee joints of each volunteer were fixed. Therefore, during lateral bending movement, the body center of gravity may be located in the upper segment, and the craniolumbar level is more active. At the same time, we observed that the rotating joint showed strong coupling motion, and the maximum rotation range was not even the main rotation axis but the coupling rotation in the front and rear directions. The possible reason for this result is related to the morphological transformation of facets. The direction of the articular surface is transformed from the sagittal plane and vertical direction to the coronal plane and horizontal direction from top to bottom, which leads to complex coupling motion during rotation. At the same time, this can also confirm that the direction of displacement is mainly in the craniocaudal direction.
From the perspective of LFJ biomechanics, paired facets and intervertebral disks jointly support and stabilize the spine and reduce spinal injury by limiting the movement of all motion planes. This plays a very important role in the load-bearing movement and stability of the whole spine . Under normal human load, LFJs can carry 3–25% of the human load, and the rest is carried by intervertebral disks and ligaments. However, if osteoarthritis occurs, the joint load of abnormal joint force can be as high as 47% , which also increases the risk of low back pain. Yin et al.  found that with the increase in facet degeneration, the initial increase in lumbar range of motion decreased, which would change the stability of the lumbar spine. Therefore, it is very important to delay facet degeneration. The literature shows that lumbar facet joint asymmetry may be one of the risk factors for intervertebral disk degeneration . In our study, asymmetric displacement of the left and right articular processes was observed in all movements (p < 0.05). Any asymmetric displacement of the facet joint will increase the risk of intercartilage impact, leading to osteophyte formation, cartilage erosion, fiber formation or shedding of the facet joint, narrowing of the joint space and nerve foramen, and the formation of synovial cysts, resulting in disorder of the movement mode of the facet joint, resulting in lower back pain.
The findings of this study may have important clinical significance. There is sufficient literature to show that the overweight movement of the lumbar spine will lead to the loss of articular cartilage water, which cannot effectively reduce the impact load and eventually lead to the initial destruction of subchondral bone and osteoarthritis . Therefore, this study can provide accurate data on the internal motion state in the process of joint weight bearing to understand the mechanism of spinal disease and guide us to carry out correct motion postures and reasonable weight-bearing ranges in daily life. At the same time, these data can also provide an important reference for the selection of surgical methods and the design of new spinal implants in lumbar surgery to reduce the risk of adjacent vertebral disease after lumbar surgery.
Our study also has some limitations. First, the samples in this study were healthy adults. In future studies, relevant pathological control groups should be added to explore the direct relationship between disease and weight bearing. Second, due to the limited visual field, we only studied the range of motion of the L3-S1 horizontal segment. Third, the relationship between the intervertebral disk and facet joint should be discussed. In future, the common changes between them after weight bearing should be further studied by the finite element method. Despite the above limitations, our study still provides accurate data on the in vivo movement of facet joints in different postures under weight-bearing conditions.
This in vivo study provides us with accurate data on the weight-bearing of facet joints in various body positions. It can guide us how to avoid waist injury due to incorrect weightbearing in our daily life, and provide a theoretical basis for studying the pathological changes of spine diseases, and provide reference for guiding clinical surgical procedures and the selection of new spinal implants.
Availability of data and materials
Please contact the corresponding author for data requests.
Lumbar facet joints
Dual fluoroscopy imaging system
Ranges of motion
6 Degrees of freedom
Digital Imaging and Communications in Medicine
Haher TR, O’Brien M, Dryer JW, Nucci R, Zipnick R, Leone DJ. The role of the lumbar facet joints in spinal stability Identification of alternative paths of loading. Spine. 1994;19(23):2667–70.
Shirazi-Adl A, Drouin G. Load-bearing role of facets in a lumbar segment under sagittal plane loadings. J Biomech. 1987;20(6):601–13.
Adams MA, Hutton WC. The mechanical function of the lumbar apophyseal joints. Spine. 1983;8(3):327–30.
Cohen SP, Raja SN. Pathogenesis, diagnosis, and treatment of lumbar zygapophysial (facet) joint pain. J Am Soc Anesthesiol. 2007;106(3):591–614.
Friedrich KM, Nemec S, Peloschek P, Pinker K, Weber M, Trattnig S. The prevalence of lumbar facet joint edema in patients with low back pain. Skeletal Radiol. 2007;36(8):755–60.
Manchikanti L, Hirsch JA, Pampati V. Chronic low back pain of facet (zygapophysial) joint origin: is there a difference based on involvement of single or multiple spinal regions? Pain Physician. 2003;6(4):399–405.
Hoogendoorn WE, Bongers PM, De Vet HCW, Douwes M, Koes BW, Miedema MC, et al. Flexion and rotation of the trunk and lifting at work are risk factors for low back pain: results of a prospective cohort study. Spine (Philadelphia, Pa 1976). 2000;25(23):3087–92.
Bigos SJ, Spengler DM, Martin NA, Zeh J, Fisher L, Nachemson A, et al. Back injuries in industry: a retrospective study. II. Injury factors. Spine. 1986;11(3):246–51.
Wood KB, Schendel MJ, Pashman RS, Buttermann GR, Lewis JL, Ogilvie JW, et al. In vivo analysis of canine intervertebral and facet motion. Spine. 1992;17(10):1180–6.
Buttermann GR, Schendel MJ, Kahmann RD, Lewis JL, Bradford DS. In vivo facet joint loading of the canine lumbar spine. Spine. 1992;17(1):81–92.
Ochia RS, Inoue N, Takatori R, Andersson GBJ, An HS. In vivo measurements of lumbar segmental motion during axial rotation in asymptomatic and chronic low back pain male subjects; 2007. p. 1394–9.
Simon P, Orías AAE, Andersson GB, An HS, Inoue N. In vivo topographic analysis of lumbar facet joint space width distribution in healthy and symptomatic subjects. Spine. 2012;37(12):1058.
Svedmark P, Tullberg T, Noz ME, Maguire G, Zeleznik MP, Weidenhielm L, et al. Three-dimensional movements of the lumbar spine facet joints and segmental movements: in vivo examinations of normal subjects with a new non-invasive method. Eur Spine J. 2012;21(4):599–605.
Fujiwara A, Tamai K, Yamato M, An HS, Yoshida H, Saotome K, et al. The relationship between facet joint osteoarthritis and disc degeneration of the lumbar spine: an MRI study. Eur Spine J. 1999;8(5):396–401.
Kozanek M, Wang S, Passias PG, Xia Q, Li G, Bono CM, et al. Range of motion and orientation of the lumbar facet joints in vivo. Spine (Philadelphia, Pa 1976). 2009;34(19):E689–96.
Yin J, Liu Z, Li C, Luo S, Lai Q, Wang S, et al. Effect of facet-joint degeneration on the in vivo motion of the lower lumbar spine. J Orthop Surg Res. 2020;15(1):340.
Bai J, Hu Y, Du L, He J, Kai L, Liu Z, et al. Assessing validation of dual fluoroscopic image matching method for measurement of in vivo spine kinematics. Chin Med J-Peking. 2011;124(11):1689–94.
Wang S, Passias P, Li G, Li G, Wood K. Measurement of vertebral kinematics using noninvasive image matching method-validation and application.; 2008. p. E355–61.
Chowdhury SK, Byrne RM, Zhou Y, Zhang X. Lumbar facet joint kinematics and load effects during dynamic lifting. Hum Factors. 2018;60(8):1130–45.
Claeson AA, Barocas VH. Computer simulation of lumbar flexion shows shear of the facet capsular ligament. Spine J. 2017;17(1):109–19.
Yang KH, King A. Mechanism of facet load transmission as a hypothesis for low-back pain. Spine. 1984;9(6):557–65.
Karacan I, Aydin T, Sahin Z, Cidem M, Koyuncu H, Aktas I, et al. Facet angles in lumbar disc herniation: their relation to anthropometric features. Spine. 2004;29(10):1132–6.
Workman J, Thambyah A, Broom N. The influence of early degenerative changes on the vulnerability of articular cartilage to impact-induced injury. Clin Biomech. 2017;43:40–9.
Supported by partially National Natural Science Foundation of China (81472140, 81572199).
Ethics approval and consent to participate
The ethical committee approved the study of the participating hospitals. All subjects signed informed consent by each patient. All clinical investigations had been conducted according to the principles expressed in the Declaration of Helsinki.
Consent for publication
Written informed consent for publication was obtained from all participants.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wen, W., Xu, H., Zhang, Z. et al. The effect of various weight-bearing activities on the motion of lumbar facet joints in vivo. J Orthop Surg Res 17, 114 (2022). https://doi.org/10.1186/s13018-022-03016-5
- Lumbar facet joints
- Range of motion
- Body position