Skip to main content

Spine system changes in soldiers after load carriage training in a plateau environment: a prediction model research



Low back pain is the most common spinal disorder among soldiers, and load carriage training (LCT) is considered the main cause. We aimed to investigate changes in the spine system of soldiers after LCT at high altitudes and the change trend of the lumbar spine and surrounding soft tissues under different load conditions.


Magnetic resonance imaging scans of the lumbar spines of nine soldiers from plateau troops were collected and processed. We used ImageJ and Surgimap software to analyze changes in the lumbar paraspinal muscles, intervertebral discs (IVDs), intervertebral foramina, and curvature. Furthermore, the multiple linear regression equation for spine injury owing to LCT at high altitudes was established as the mathematical prediction model using SPSS Statistics version 23.0 software.


In the paraspinal muscles, the cross-sectional area (CSA) increased significantly from 9126.4 ± 691.6 mm2 to 9862.7 ± 456.4 mm2, and the functional CSA (FCSA) increased significantly from 8089.6 ± 707.7 mm2 to 8747.9 ± 426.2 mm2 after LCT (P < 0.05); however, the FCSA/CSA was not significantly different. Regarding IVD, the total lumbar spine showed a decreasing trend after LCT with a significant difference (P < 0.05). Regarding the lumbar intervertebral foramen, the percentage of the effective intervertebral foraminal area of L3/4 significantly decreased from 91.6 ± 2.0 to 88.1% ± 2.9% (P < 0.05). For curvature, the lumbosacral angle after LCT (32.4° ± 6.8°) was significantly higher (P < 0.05) than that before LCT (26.6° ± 5.3°), while the lumbar lordosis angle increased significantly from 24.0° ± 7.1° to 30.6° ± 7.4° (P < 0.05). The linear regression equation of the change rate, FCSA% = − 0.718 + 23.085 × load weight, was successfully established as a prediction model of spinal injury after LCT at high altitudes.


The spinal system encountered increased muscle volume, muscle congestion, tissue edema, IVD compression, decreased effective intervertebral foramen area, and increased lumbar curvature after LCT, which revealed important pathophysiological mechanisms of lumbar spinal disorders in soldiers following short-term and high-load weight training. The injury prediction model of the spinal system confirmed that a load weight < 60% of soldiers’ weight cannot cause acute pathological injury after short-term LCT, providing a reference supporting the formulation of the load weight standard for LCT.


Chronic low back pain is the most common musculoskeletal disorder among soldiers. According to statistics, > 50% of soldiers have experienced low back pain during service [1], and the prevalence rate of low back pain in new recruits during basic training reaches 30.3% [2]. Chronic low back pain is also the most common illness in the United States (US) military that seriously affects the normal training of soldiers and reduces the service life and quality of active-duty soldiers [3]. At present, commonly considered main risk factors of low back pain in soldiers during military training include weight, training environment, and psychological and organizational management factors. Among them, excessive load carriage training (LCT) is considered the main cause of back pain and injuries among soldiers [4].

In LCT, a soldier carries a certain weight of marching combat equipment according to the requirements of the mission, battlefield environment, weather conditions, and other factors when performing training or combat tasks. Then, the soldiers carry different weights across various terrains to complete various operational tasks to a high standard. A previous study described that under general temperature conditions when the marching speed is 5 km/h, the maximum allowable load weight for soldiers is 25 kg, and the most acceptable weight is 20 kg [5]. However, during actual training or combat, soldiers always carry heavier loads than the current standards set by the military. For example, the US military standard requires 25 kg as the maximum load weight of field soldiers in marches; however, during the Iraq war, the actual load weight of US soldiers was approximately 35 kg, which could reach 42–75 kg when soldiers were approaching the enemy [6].

At present, “plateau adaptive training” is an important training method for endurance and perseverance. In a high-altitude environment (> 500 m above sea level), some factors such as low partial pressure of air and oxygen caused by high altitude and dry and cold climate reduce the alveolar oxygen partial pressure, blood oxygen saturation, energy consumption, and metabolic rate of soldiers during training. Thus, both the maximum and acceptable load weights are reduced accordingly. Yin et al. [7] suggested that at a constant marching speed of 4 km/h, when the altitude of the training environment increased from 3700 m to 4300 m, the appropriate load weight of soldiers should be reduced from 20.0 kg to 11.5 kg, and the reduction rate is close to 50%. To deal with the adverse effects of load weight reduction on military training in the plateau environment, at present, the military mainly increases the maximum load weight of soldiers to improve their aerobic endurance so as to strengthen the plateau LCT. However, this often causes an overload on the spine system, accelerates strain on the paraspinal muscles and intervertebral discs (IVDs), and even leads to changes in the spine curvature and spinal cord compression [8, 9]. Hrubec et al. [10] reported that 1095 soldiers had been treated for low back pain after World War II. Recent studies have found that soldiers have a higher risk of lumbar disc herniation than ordinary people. Moreover, some studies have reported that as the load weight of soldiers increases, the forward tilt angle of the lumbar spine also increases, which in turn affects the balance and increases the risk for a spinal disorder [11].

To the best of our knowledge, there is a lack of standards for plateau LCT. Existing military studies on plateau training have focused more on respiration, sleep, or psychological problems of soldiers [12], whereas little attention has been paid to changes in the spine system of soldiers who carry varying load weights during plateau training. Therefore, we aimed to investigate changes in the spine system of soldiers during plateau LCT and the trend of the lumbar spine change under different load conditions. By obtaining radiographic lumbar images of soldiers, we established a prediction model for lumbar spine injury after LCT at high altitude.


Experimental design

Nine soldiers from combat troops in Xining City, Qinghai Province (2443 m above sea level), volunteered for this study. All subjects were of Han nationality with the following characteristics: age, 18–31 years; weight, 60–80 kg; height, 168–182 cm; and service life, 1–10 years on the plateau (Table 1). The training scores and muscle strength of these soldiers were up to the excellent standard. Since the interventions in this study may lead to injury, we ensured that soldiers were trained safely. With the limited military training schedule and secrecy requirements, we could not directly obtain image data of the spines from an adequate number of soldiers. To simulate the LCT of plateau soldiers with different load weights and time, we used a uniform experimental design (Table 2). The load weight (percentage of soldiers’ weight) was Factor 1, which included 20, 40, and 60% of the soldiers’ weight, and the load time was Factor 2, which included 20 min, 40 min, and 60 min. The subjects were grouped by draw lots as a completely randomized grouping. In this experimental design, one participant was assigned to one test condition (i.e., subject 1 had a 40-min LCT with 20% load weight).

Table 1 Basic information of the subjects
Table 2 Uniform experimental design of plateau LCT

This study was approved by the ethics committee of the Qinghai Red Cross Hospital (KY-2020-15), and each subject provided written informed consent.

Research equipment

Lumbar spine magnetic resonance imaging (MRI) was performed before and after the LCT by using the GE Signa Explorer 1.5 T MRI system at the Qinghai Red Cross Hospital. As for the examination time, some studies have shown variations on the lumbar flexion range, lumbar IVD height, and other biomechanical characteristics at different times of the day, known as the diurnal effect [13, 14]. To avoid this effect, all soldiers were trained uniformly in the morning and MRI was completed within 30 min after LCT.

Because of the limited supine position during training, we could not directly obtain actual spine images of soldiers who were standing with different load weights after training. Therefore, a pressure device was designed in advance for this study to simulate the soldiers’ loaded backs when they were in the lying position (Fig. 1). Detailed steps were as follows: we used nylon wire through the Velcro fixing hole (7) to immobilize the legs, and a rubber belt was used to link the subject’s shoulder to the adjusting mechanism (4) through the mounting hole (5) (Fig. 1). The adjusting mechanism was then used to adjust the length of the belt to simulate the load weight on the back. This device is made of non-metallic materials and can be used during MRI to obtain images of the spine under different load weights. The effectiveness of this device was verified by simulation mechanics analysis using AnyBody Modeling System, version 4.0 software (AnyBody Technology, Aalborg, Denmark), that is, it can simulate the loaded state of the human body in the standing position (Fig. 2).

Fig. 1
figure 1

Pressure device for simulating the soldier back loaded in the supine position

Fig. 2
figure 2

Simulation mechanics analysis of the human body with load weight. a The simulation model with load weight by carrying a backpack. b The simulation model with load weight using the pressure device

According to the experimental design, to perform MRI before LCT, the subjects were lying supine on the device without support on the lumbar spine. Then, each subject carried a backpack weighting a percentage of their weight during LCT. After LCT, the subject was fixed to the device in the same posture, and the spring dynamometer was attached to one end of the rubber belt to apply a tension pressure, which is equal to the subject’s backpack weight, as the load weight. Finally, we affixed the rubber belt and performed MRI.

Research indicators

Lumbar paraspinal muscle

Changes in the paraspinal muscles can be detected based on the cross-sectional area (CSA), fat content, and proportions of fast-twitch fibers [15]. The CSA image was obtained from T2-weighted MRI scans. Because the boundary of the paraspinal muscle at the L3/4 vertebral level can be clearly identified, we used ImageJ version 1.51 software (National Institute of Mental Health, Bethesda, MD, USA) to manually mark the contour of the paraspinal muscle at this level and then summed the CSAs of the following paraspinal muscles: bilateral erector spinae, bilateral psoas major, and bilateral quadratus lumborum (Fig. 3). The CSA was used to quantify the level of hypertrophy or atrophy of the lumbar paraspinal muscles [16, 17].

Fig. 3
figure 3

Research indicators of the lumbar spine and measurement methods. a The location of the lumbar paraspinal muscles identified for the cross-sectional area (CSA) and functional cross-sectional area (FCSA) at the L3/4 level. b The measurement for intervertebral discs (IVD) height (anterior, middle, and posterior) from L1 to S1. c The measurement for the intervertebral foramen area (IF) and the effective intervertebral foraminal area (E-IF) in L3/4 segments. d The measurement for the lumbosacral angle (LSA) and the lumbar lordosis angle (LLA)

At the same time, Parkkola et al. [18] confirmed that when the paraspinal muscles undergo physiological changes, the CSA may not change correspondingly because some tissue fluids, fats, or other fibrous tissues fill the original locations of the muscle fibers. Ranson et al. [19] reported that the gray value threshold of the muscle fibers in T2-weighted MRI scans was approximately 0–120 pixels (30–240 pixels for the bones and 120–660 pixels for the fats). Based on this fact, we obtained the gray values of images of the paraspinal muscles before and after LCT. Then, to determine the functional CSA (FCSA) of the paraspinal muscle, the threshold analysis method was used to select the image area in the CSA within the gray value threshold of 0–120 pixels. FCSA was used to represent the total content of lean muscle fibers within the boundaries of the lumbar paraspinal muscle fascia to measure changes in the muscle fibers. A related control study by Fortin et al. [20] showed that this image measurement method has certain reliability and repeatability regardless of whether it is performed by one or more people, and the intra-group correlation coefficient is 0.99.

Lumbar IVD height

We measured the height of the lumbar IVD on T2-weighted MRI scans at the midsagittal plane with a slice thickness of 3.5 mm. According to the method of Dabbs, we measured the disc height at the anterior, middle, and posterior sections at L1/2, L2/3, L3/4, L4/5, and L5/S1 levels before and after LCT [21] (Fig. 3). The disc height was defined as the average of three height values at a given level. Finally, the difference in IVD height before and after LCT and the change in the total IVD height were compared.

Lumbar intervertebral foramen

We measured the lumbar intervertebral foramen on T2-weighted MRI scans at the sagittal plane. Clinical investigations showed that the intervertebral foramina in the L3/4 and L4/5 segments are located at the physiological bend of the lumbar spine. Due to stress concentration during weight loading, the annulus fibrosus is prone to damage under the action of external forces; thus, the degeneration in these segments is more severe than in other segments [22, 23]. Therefore, we used the measurement method similar to CSA: we artificially marked the contour of the L3/4 and L4/5 intervertebral foramina on both the left and right sides (the upper, lower, and outer borders of the intervertebral foramina are bone boundaries, and the inner boundary is the line connecting the upper and lower vertebral bodies, Fig. 3) and measured this area as the intervertebral foraminal (IF) area. Then, we kept the upper, lower, and outer borders unchanged, marked the inner border again as the IVD boundary (Fig. 3), and measured this area as the effective IF (E-IF) area. Finally, the ratio of these two areas (E-IF/IF) was defined as the percentage of the effective IF area (E-IF%), which represented the effective space rate for nerve root activity in the lumbar intervertebral foramen.

Lumbar curvature

We measured the lumbosacral angle (LSA) and lumbar lordosis in the midsagittal T2-weighted image using Surgimap version 2.2.15 software (Nemaris, NY, USA) to evaluate the lumbar curvature. The LSA was measured as the angle of intersection between the extension of the upper margin of the sacrum and the horizontal line [24]. Because the lumbosacral region is the joint point of stress when the body is upright or loaded, this angle can indicate the stability of the lumbar spine. The lumbar lordosis angle (LLA) was measured using Cobb’s method, which measures the angle between the L1 vertebral upper endplate and L5 vertebral lower endplate extension line (Fig. 3) [25].

Questionnaire survey

The subjective questionnaire survey was performed and analyzed on nine participants. The questionnaire included basic information, degree of low back pain, degree of fatigue, and other physical symptoms. The low back pain level was rated by using an 11-point numerical rating scale (NRS) (scored from 0 to 10, with 0 as painless, 1–3 as mild pain, 4–6 as moderate pain, 7–9 as severe pain, 10 as twinge), and the degree of fatigue was self-rated from 0 to 10 (0 as energetic, 10 as exhausted), which means that the higher the score, the more the physical exertion [26, 27]. All questionnaires were completed independently by the subjects before and after LCT.

Statistical analysis

All data of each research index were analyzed using SPSS Statistics version 23.0 (IBM Corp., Armonk, NY, USA). First, the mean value analysis was used to evaluate the characteristics of changes before and after LCT. Second, significance analysis (P < 0.05) and simulation regression analysis were performed for the change rates (Index%) of different indexes (LCT% = (post LCT value – pre LCT value)/pre LCT value). Detailed steps are as follows: 1) Shapiro-Wilk tests were used to check data normality (P > 0.05). 2) Paired sample t-tests were used to compare significant differences between each index before and after LCT (P < 0.05). 3) The main effect variance analysis method was used to analyze whether there were significant differences in the changes in various indicators between groups with different load weights and between those with different load times (P < 0.05). 4) Multiple comparison analysis was used to further compare differences in the index changes in the same load weight group or the same load time group (P < 0.05). 4) A specific index was selected to fit the multiple linear regression equation: the change rate before and after LCT = b1 + b2 × load weight + b3 × load time, which we used as a mathematical prediction model.


Lumbar paraspinal muscle indexes

The CSA increased from 9126.371 ± 691.614 mm2 to 9862.655 ± 456.416 mm2 after LCT, with an average increase of 8.07%. The paired t-test showed that the CSA was significantly different before and after LCT (t = − 6.368, P = 0.0002, Table 3). The average gray value of the paraspinal muscle was reduced from 80.396 ± 17.522 to 55.009 ± 10.641 after LCT, with an average decrease of 31.58%. The paired t-test showed significant differences before and after LCT (t = 5.968, P = 0.0003, Table 3). The FCSA before LCT was 8089.556 ± 707.739 mm2, and it increased to 8747.917 ± 426.245 mm2 after LCT, with an average increase of 8.14%. The paired t-test showed that the FCSA was significantly different before and after LCT (t = − 5.348, P = 0.0007, Table 3).

Table 3 Statistical analysis results of lumbar spine indexes

The FCSA/CSA changed from 0.886 ± 0.032 to 0.887 ± 0.024 after LCT, with an average change of 0.12%. No significant difference was found before and after LCT (t = −0.169, P = 0.8703). After the analysis of the main effect variance, no significant difference in the parameters of the lumbar paraspinal muscle between the two groups was found, regardless of whether they were grouped by different load weight or load time (Table 4).

Table 4 Results of the main effect variance analysis

Lumbar IVD height indexes

The disc height at each level of L1-S1 and the total height showed a consistent decreasing trend after LCT. The total disc height decreased from 87.243 ± 8.943 pixels to 84.759 ± 7.620 pixels, with an average decrease of 2.85%. The paired t-test showed a significant difference in the total disc height before and after LCT (t = 2.971, P = 0.0180, Table 3). Moreover, a significant difference was found in the analysis results of the disc height at L1/2 (t = 3.308, P = 0.0107) and L4/5 level (t = 3.509, P = 0.0080). A related study reported that the standard deviations of this measurement method between observers and within observers are 0.2 mm and 0.3 mm, respectively [28]. After the analysis of the main effect variance, no significant differences in the total height or each height at different levels of the lumbar IVD were found between the two groups, regardless of whether they were grouped by load weight or load time (Table 4).

Lumbar IF indexes

The E-IF% of L3/4 decreased from 91.6% ± 2.0 to 88.1% ± 2.9% after LCT, with an average reduction of 3.80%. The paired t-test indicated a significant difference in E-IF% of L3/4 (t = 8.214, P = 0.0004). The E-IF% of L4/5 decreased from 89.6% ± 1.9 to 88.2% ± 1.7% after LCT, with an average reduction of 1.66%. The paired t-test indicated no significant difference in this index (t = 2.553, P = 0.0511, Table 3). After the analysis of the main effect variance, no significant differences were noted in the E-IF% of L3/4 and L4/5 between the two groups, regardless of whether they were grouped by load weight or load time (Table 4).

Lumbar curvature index

The LSA increased from 26.622° ± 5.303° to 32.422° ± 6.816° after LCT, with an average increase of 21.79%. The paired t-test showed a significant difference in LSA (t = − 4.911, P = 0.0012, Table 3). The LLA increased from 24.011° ± 7.117° to 30.644° ± 7.415° after LCT, with an average increase of 27.63%. The paired t-test indicated a significant difference (t = − 4.851, P = 0.0013, Table 3).

The results of the main effect variance analysis showed no significant difference between the groups with different load times. However, a significant difference in the LSA was noted among different load weight groups (F = 8.228, P = 0.0191), while the LLA showed no significant difference (Table 4).

Questionnaire results

The NRS score results showed that eight soldiers had mild back pain (NRS score from 1 to 3) after LCT, while only one soldier suffered moderate back pain (NRS score 5). The fatigue degree score results indicated that all soldiers performed intense physical activities (fatigue degree score > 6) after LCT, and two of them felt totally exhausted (fatigue degree score 10) (Table 5).

Table 5 Results of the questionnaire survey (n)

Prediction model for spinal injury

To establish a prediction model of spinal injury for LCT, we selected the FCSA which had a significant difference and a certain clinical significance to fit the multiple linear regression equation. The result was as follows:

As regards the load weight, the results of the Pearson correlation coefficient analysis showed that the correlation coefficient between the FCSA% and load weight was 0.761, which was a high value (close to 0.8), indicating that the correlation is strong, and a linear relationship can be established. With regard to the load time, the correlation coefficients of FCSA% was − 0.176, and the relative value was low (< 0.6), suggesting the absence of a linear relationship.

To fit the multiple linear regression equation of FCSA%, the load weight value was introduced into the prediction model as a variable; the result of the analysis of variance of the multiple regression indicated that the F-value of this model was 9.612 with a P-value of 0.017 (< 0.1), indicating a linear relationship between FCSA% and load weight. In the regression coefficient analysis, the constant term was − 0.718, and the partial regression coefficient b1 was 23.085 after the t-test (P = 0.017 (< 0.1)). The regression coefficient of load time outside the equation of the model was re-tested, and its P-value was 0.515 (> 0.1) which could not be introduced. The maximum absolute value of the standardized residual of this model was 1.422 (< 3). The normal p-p plot of regression standardized residual showed that the scattered points were basically distributed on the diagonal, which indicated that the residuals could be judged to be normally distributed (Fig. 4). Finally, we obtained the multiple linear regression equation, FCSA% = − 0.718 + 23.085 × load weight, which can be used to describe the influence of load weight on the change rate of FCSA during LCT.

Fig. 4
figure 4

Normal p-p plot of regression standardized residual


Changes to the spine system

To compensate for plateau LCT, certain anatomical or physiological changes will occur in the lumbar spine system. The results showed that both CSA and FCSA had an increasing trend after LCT. However, no significant change was noted in the FCSA/CSA after LCT, proving that the proportion of the lean muscle mass remained basically unchanged after acute high-intensity loading and confirming that the mass and volume of the muscle fibers did not increase after short-term LCT according to the imaging study [29, 30]. That finding indicates that the volume of the paraspinal muscle tissue increased after loading, which may be clinically related to some physiological changes such as muscle congestion and tissue edema. Similarly, the gray value of the paraspinal muscles decreased significantly after LCT, indicating that the mass of the muscle fibers did not increase again. The lumbar paraspinal muscles function by resisting load, stretching or flexing the spine, maintaining the stability of the spine joints, and protecting the discs and ligaments from excessive strain [31]. Excessive LCT with long duration and short rest can lead to aseptic inflammatory reactions, such as muscle cell necrosis, muscle fiber rupture, tissue fluid exudation, and lactic acid accumulation, that result in training injuries of the paraspinal muscles [32].

In terms of lumbar IVDs, the results showed that the total height of the discs decreased after LCT, and the paired t-test results showed a significant difference in the change in the IVD height, which proved that the entire spinal system was compressed after loading and the IVD provided a certain cushioning effect. However, in view of clinical medicine, IVDs generally provide the secondary support to the spine during weight loading. The specific mechanism involves sensing the position of the vertebral body through nerve feedback to activate the contraction of the paraspinal muscles during movement and accomplishes stabilization and support of the articular structures [33, 34]. From an anatomical perspective, the bone, ligaments, and other tissues around the IVD have certain limitations on its deformation. Therefore, the height change of the IVD is limited under the loading state that the body can accept. Meanwhile, LCT significantly increases repetitive cycling pressure on the IVD, and the force perpendicular to the ground is also transmitted to the IVD. These mechanical factors eventually cause damage to the IVD [35].

Regarding the lumbar intervertebral foramen, the results showed that the E-IF% in the L3/4 and L4/5 levels decreased after LCT. From an anatomical view, when pathophysiological changes such as disc herniation, facet joint swelling, synovial swelling, or bony hyperplasia occur, foraminal stenosis may develop. When the E-IF is less than the diameter of the nerve root, nerve root compression occurs and manifested as pain in the lower back and buttocks and radiating pain or numbness in the limbs; this is an important mechanism for the development of lumbar spine diseases [36, 37]. Mayoux-Benhamou et al. [38] suggested that when the height of the IVD is reduced to 4 mm, that is, 35–50% of the normal height, the reduction in the foraminal area causes nerve compression. However, our results showed a 1-mm reduction in IVD after LCT, which cannot cause compression. Moreover, there is a known correlation between the IVD and the narrowest area of the neuroforamen (R2 = 0.56) [39]. Our study confirmed that for soldiers without previous spinal degenerative disease, posterior IVD herniation generated by LCT will occupy the available space of the intervertebral foramina and even cause irreversible lesions in severe cases, which is an important pathogenesis of training-related injury symptoms.

As for the lumbar curvature, the LLA increased by 27.63%, while the LSA increased by 21.79% after LCT in the midsagittal image, both showing an increasing trend and having significant differences after LCT. Anatomically, the LLA and LSA can indirectly reflect the stability of the lumbar spine [40, 41]. If the LLA or LSA is too large, it will affects the balance of the spinal system, increases the load on structures such as the paraspinal muscles and IVDs, and increases the probability of injury [42]. The increases in these two angles confirmed the compensatory effect of the lumbar curvature on weight loading. However, the variation in the lumbar curvature was confirmed that it was not the main compensatory mechanism for loading and had a certain upper limit. Moreover, the comparison results of the LSA between groups with different load weights showed significant differences, confirming that the change rate of the LSA has a certain correlation with the load weight that has the potential to become an evaluation indicator for LCT.

In accordance with the questionnaire results, although we found that all soldiers suffered varying degrees of low back pain after LCT, most of them just experienced mild pain and only one soldier experienced moderate pain [43]. These back pain symptoms related to LCT could be relieved by rest. Moreover, the fatigue degree score showed that all soldiers performed intense physical activities during LCT, which indicated a high intensity of the current training. In addition, none of soldiers had other physical symptoms after LCT. Beyond that, the radiologists reviewed all MRI scans after examination and did not find evidence of spinal injury. These results confirmed that the current LCT with different experimental conditions achieved desirable training effect without causing acute irreversible injury to the lumbar spine in this study.

Prediction model

The linear regression equation of the change rate, FCSA% = − 0.718 + 23.085 × load weight, was successfully established as a predictive model of soldiers’ spinal injury after plateau LCT. According to this equation, we believe that as the load weight increases, the FCSA% will increase gradually, and the FCSA value will present an increasing trend. Therefore, we can use this prediction model to estimate that when the load weight reaches 100% of the soldier’s weight, the FCSA is expected to increase by approximately 22.367% after LCT. Since FCSA represents the lean muscle volume within the paraspinal muscle fascia, the variation trend of the prediction model further confirmed that plateau LCT can induce a certain proportion of hypertrophy and proliferation of muscle fiber cells. However, according to the physiological structure of the paraspinal muscle, the FCSA cannot increase indefinitely, and both the CSA and FCSA will eventually reach a peak with the increase in training intensity. Therefore, as the load weight continues to increase, we speculate that the partial regression coefficient b1 (23.085) of this prediction model will gradually decrease and the constant term (− 0.718) will gradually increase; finally, the FCSA% will reach the upper limit, that is, the FCSA will reach the maximum threshold. However, this prediction model has not yet covered this range.

According to the results of questionnaire survey and radiologists review, it can be concluded that when the maximum load weight does not exceed 60% of their weight, soldiers will not experience paraspinal muscle injury related to short-term LCT, proving that a load weight < 60% is within the tolerable load range of the lumbar paraspinal muscles. As a consequence, this prediction model can be applied to the current conventional weight-bearing training for soldiers in the plateau (the proportion of weight-bearing is less than 60%), which can provide a reference for the formulation of the load weight standard for plateau soldiers under training or combat and improve the scientificity and safety of LCT in the plateau environment.


This study performed actual human experiments in plateau troops for LCT. Results should be interpreted in consideration of the limitations. First, due to military requirements and ethical principles, this experiment could not recruit a sufficient number of subjects. Although a uniform experimental design was used, it cannot completely offset the bias of test data and differences between subjects caused by the small sample size. We hope to recruit more subjects in subsequent studies to verify the present results and obtain further physiological changes in the spine system following plateau LCT. Second, due to ethical restrictions, we were unable to conduct LCT on soldiers using a load weight > 60% or a load time longer than 60 min at present; therefore, we could not directly determine the load weight and load time that may cause physiological damage. Thus, our prediction model also cannot cover the threshold range. We will further validate and optimize this model in subsequent experiments and estimate the training threshold through a more accurate mathematical model.


This study directly confirmed that the spinal system encountered the following changes after short-term and high-load weight training: increased muscle volume, muscle congestion, tissue edema, IVD compression, decreased E-IF area, increased lumbar curvature, and other physiological changes that reveal important pathophysiological mechanisms of lumbar spinal disorders in soldiers. Moreover, this study established an injury prediction model of the spinal system for plateau LCT through multiple regression analysis: FCSA% = − 0.718 + 23.085 × load weight. Combined with the results of the questionnaire survey and the radiologists review, this model confirmed that a load weight < 60% of soldiers’ weight cannot cause acute pathological injury in short-term LCT, providing a reference supporting the formulation of the load weight standard for the LCT and, thus, improving the scientificity and safety of military LCT in plateau environments. In addition, the simulated weight load device, image processing methods, and regression analysis method involved in this research have certain universality that can be applied to clinical research of spinal diseases and have considerable application prospects.

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available due to the individual privacy, but are available from the corresponding author on reasonable request.



Cross-sectional area


Percentage of effective intervertebral foramen area


Effective intervertebral foraminal


Functional cross-sectional area


Intervertebral foraminal


Intervertebral disc


Load carriage training


Lumbar lordosis angle


Lumbosacral angle


Magnetic resonance imaging


United States


  1. Larson MJ, Adams RS, Ritter GA, Linton A, Williams TV, Saadoun M, et al. Associations of early treatments for low-back pain with military readiness outcomes. J Altern Complement Med. 2018;24(7):666–76.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Moshe S, Zack O, Finestone AS, Mishal M, Segal N, Slodownik D, et al. The incidence and worsening of newly diagnosed low back pain in a population of young male military recruits. BMC Musculoskelet Disord. 2016;17:279.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Childs JD, Wu SS, Teyhen DS, Robinson ME, George SZ. Prevention of low back pain in the military cluster randomized trial: effects of brief psychosocial education on total and low back pain-related health care costs. Spine J. 2014;14(4):571–83.

    Article  PubMed  Google Scholar 

  4. Roy TC, Lopez HP, Piva SR. Loads worn by soldiers predict episodes of low back pain during deployment to Afghanistan. Spine (Phila Pa 1976). 2013;38(15):1310–7.

    Article  Google Scholar 

  5. Thelen MD, Koppenhaver SL, Gill NW, Shaffer SW. Effects of technique-focused training in conjunction with physical readiness training on Army physical fitness test performance. US Army Med Dep J. 2017;(3-17):79-83.

  6. Mala J, Szivak TK, Kraemer WJ. Improving performance of heavy load carriage during high-intensity combat-related tasks. Strength Cond J. 2015;37(4):43–52.

    Article  Google Scholar 

  7. Yin ZY, Zan JP, Sun XB, Nan WK, Wang HX, Wang CY, et al. A study on optimum load carriage of men at high altitudes. Space Med Med Eng. 1994;7(3):202–8.

    Google Scholar 

  8. Orr RM, Pope R, Johnston V, Coyle J. Soldier occupational load carriage: a narrative review of associated injuries. Int J Inj Control Saf Promot. 2014;21(4):388–96.

    Article  Google Scholar 

  9. Yin ZY, Zan JP, Ma Z, Nan WK, Sun XB, Wang HX, et al. Study of optimum load carriage of soldiers at high altitude of 3700 m. J Prev Med Chin Peopl Liber Army. 1994;12(6):422–7.

    Google Scholar 

  10. Hrubec Z, Nashold BS Jr. Epidemiology of lumbar disc lesions in the military in world war II. Am J Epidemiol. 1975;102(5):367–76.

    Article  CAS  PubMed  Google Scholar 

  11. Rodríguez-Soto AE, Berry DB, Palombo L, Valaik E, Kelly KR, Ward SR. The effect of load magnitude and distribution on lumbar spine posture in active-duty marines. Spine (Phila Pa 1976). 2017;42(5):345–51.

    Article  Google Scholar 

  12. Muza SR. Military applications of hypoxic training for high-altitude operations. Med Sci Sports Exerc. 2007;39(9):1625–31.

    Article  PubMed  Google Scholar 

  13. Adams MA, Dolan P, Hutton WC. Diurnal variations in the stresses on the lumbar spine. Spine (Phila Pa 1976). 1987;12(2):130–7.

    Article  CAS  Google Scholar 

  14. Kolf AK, Konieczny M, Hesper T, Hosalkar H, Schleich C, Antoch G, et al. T2* mapping of the adult intervertebral lumbar disc: normative data and analysis of diurnal effects. J Orthop Res. 2019;37(9):1956–62.

    Article  PubMed  Google Scholar 

  15. Mannion AF, Weber BR, Dvorak J, Grob D, Müntener M. Fibre type characteristics of the lumbar paraspinal muscles in normal healthy subjects and in patients with low back pain. J Orthop Res. 1997;15(6):881–7.

    Article  CAS  PubMed  Google Scholar 

  16. Fortin M, Macedo LG. Multifidus and paraspinal muscle group cross-sectional areas of patients with low back pain and control patients: a systematic review with a focus on blinding. Phys Ther. 2013;93(7):873–88.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Cooley JR, Walker BF, Ardakani EM, Kjaer P, Jensen TS, Hebert JJ. Relationships between paraspinal muscle morphology and neurocompressive conditions of the lumbar spine: a systematic review with meta-analysis. BMC Musculoskelet Disord. 2018;19(1):351.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Parkkola R, Kormano M. Lumbar disc and back muscle degeneration on MRI: correlation to age and body mass. J Spinal Disord. 1992;5(1):86–92.

    Article  CAS  PubMed  Google Scholar 

  19. Ranson CA, Burnett AF, Kerslake R, Batt ME, O'Sullivan PB. An investigation into the use of MR imaging to determine the functional cross sectional area of lumbar paraspinal muscles. Eur Spine J. 2006;15(6):764–73.

    Article  PubMed  Google Scholar 

  20. Fortin M, Battie MC. Quantitative paraspinal muscle measurements: inter-software reliability and agreement using OsiriX and ImageJ. Phys Ther. 2012;92(6):853–64.

    Article  PubMed  Google Scholar 

  21. Dabbs VM, Dabbs LG. Correlation between disc height narrowing and low-back pain. Spine (Phila Pa 1976). 1990;15(12):1366–9.

    Article  CAS  Google Scholar 

  22. Lao L, Daubs MD, Takahashi S, Lord EL, Cohen JR, Zhong G, et al. Kinetic magnetic resonance imaging analysis of lumbar segmental motion at levels adjacent to disc herniation. Eur Spine J. 2016;25(1):222–9.

    Article  PubMed  Google Scholar 

  23. Paholpak P, Nazareth A, Khan YA, Khan SU, Ansari F, Tamai K, et al. Evaluation of foraminal cross-sectional area in lumbar spondylolisthesis using kinematic MRI. Eur J Orthop Surg Traumatol. 2019;29(1):17–23.

    Article  PubMed  Google Scholar 

  24. Onyemaechi NO, Anyanwu GE, Obikili EN, Onwuasoigwe O, Nwankwo OE. Impact of overweight and obesity on the musculoskeletal system using lumbosacral angles. Patient Prefer Adher. 2016;10:291–6.

    Article  Google Scholar 

  25. Chen YL. Vertebral centroid measurement of lumbar lordosis compared with the cobb technique. Spine (Phila Pa 1976). 1999;24(17):1786–90.

    Article  CAS  Google Scholar 

  26. Thong ISK, Jensen MP, Miró J, Tan G. The validity of pain intensity measures: what do the NRS, VAS, VRS, and FPS-R measure? Scand J Pain. 2018;18(1):99–107.

    Article  PubMed  Google Scholar 

  27. Chiarotto A, Maxwell LJ, Ostelo RW, Boers M, Tugwell P, Terwee CB. Measurement properties of visual analogue scale, numeric rating scale and pain severity subscale of the brief pain inventory in patients with low back pain: a systematic review. J Pain. 2019;20(3):245–63.

    Article  PubMed  Google Scholar 

  28. Shymon S, Hargens AR, Minkoff LA, Chang DG. Body posture and backpack loading: an upright magnetic resonance imaging study of the adult lumbar spine. Eur Spine J. 2014;23(7):1407–13.

    Article  PubMed  Google Scholar 

  29. Davison MJ, Maly MR, Keir PJ, Hapuhennedige SM, Kron AT, Adachi JD, et al. Lean muscle volume of the thigh has a stronger relationship with muscle power than muscle strength in women with knee osteoarthritis. Clin Biomech (Bristol, Avon). 2017;41:92–7.

    Article  Google Scholar 

  30. Fortin M, Videman T, Gibbons LE, Battié MC. Paraspinal muscle morphology and composition: a 15-yr longitudinal magnetic resonance imaging study. Med Sci Sports Exerc. 2014;46(5):893–901.

    Article  CAS  PubMed  Google Scholar 

  31. Li SSW, Chow DHK. Effects of backpack load on critical changes of trunk muscle activation and lumbar spine loading during walking. Ergonomics. 2018;61(4):553–65.

    Article  PubMed  Google Scholar 

  32. Dupeyron A, Lecocq J, Vautravers P, Pélissier J, Perrey S. Muscle oxygenation and intramuscular pressure related to posture and load in back muscles. Spine J. 2009;9(9):754–9.

    Article  PubMed  Google Scholar 

  33. Leinonen V, Kankaanpää M, Luukkonen M, Hänninen O, Airaksinen O, Taimela S. Disc herniation-related back pain impairs feed-forward control of paraspinal muscles. Spine (Phila Pa 1976). 2001;26(16):E367–72.

    Article  CAS  Google Scholar 

  34. Chen WY, Wang K, Yuan WA, Zhan HS. Relationship between lumbosacral multifidus muscle and lumbar disc herniation. Chin J Orthopa Trauma. 2016;29(6):581–4 (Article in China).

    CAS  Google Scholar 

  35. Leinonen V, Kankaanpää M, Luukkonen M, Kansanen M, Hänninen O, Airaksinen O, et al. Lumbar paraspinal muscle function, perception of lumbar position, and postural control in disc herniation-related back pain. Spine (Phila Pa 1976). 2003;28(8):842–8.

    Google Scholar 

  36. Marawar SV, Madom IA, Palumbo M, Tallarico RA, Ordway NR, Metkar U, et al. Surgeon reliability for the assessment of lumbar spinal stenosis on MRI: the impact of surgeon experience. Int J Spine Surg. 2017;11(5):34.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Park HJ, Kim SS, Lee SY, Park NH, Chung EC, Rho MH, et al. A practical MRI grading system for cervical foraminal stenosis based on oblique sagittal images. Br J Radiol. 2013;86(1025):1–7.

    Google Scholar 

  38. Mayoux-Benhamou MA, Revel M, Aaron C, Chomette G, Amor B. A morphometric study of the lumbar foramen. Influence of flexion-extension movements and of isolated disc collapse. Surg Radiol Anat. 1989;11(2):97–102.

    Article  CAS  PubMed  Google Scholar 

  39. Chen D, Fay LA, Lok J, Yuan P, Edwards WT, Yuan HA. Increasing neuroforaminal volume by anterior interbody distraction in degenerative lumbar spine. Spine (Phila Pa 1976). 1995;20(1):74–9.

    Article  CAS  Google Scholar 

  40. Fatih KB, Ejder B. Comparison of lumbosacral alignment in geriatric and non-geriatric patients suffering low back pain. Pak J Med Sci. 2018;34(2):282–7.

    Google Scholar 

  41. Habibi Z, Maleki F, Meybodi AT, Mahdavi A, Saberi H. Lumbosacral sagittal alignment in association to intervertebral disc diseases. Asian Spine J. 2014;8(6):813–9.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Makirov SK, Yuz AA, Jahaf MT, Nikulina AA. Quantitative evaluation of the lumbosacral sagittal alignment in degenerative lumbar spinal stenosis. Int J Spine Surg. 2015;9:68.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Karcioglu O, Topacoglu H, Dikme O, Dikme O. A systematic review of the pain scales in adults: which to use? Am J Emerg Med. 2018;36(4):707–14.

    Article  PubMed  Google Scholar 

Download references


We would like to thank Editage ( for English language editing.


This work was supported by the National Key Research and Development Program of China (2018YFB1307603) and the National Natural Science Foundation of China (8174100706).

Author information

Authors and Affiliations



JTW, GL and SHL collected the basic and image data of soldiers. YZ and LD analyzed and interpreted the image data regarding the lumbar spine. HQ and LJY performed the statistical analysis. HQ was the major contributor in writing the manuscript. All authors read and approved the final manuscript.

Authors’ information

YZ is the professor, chief physician and doctoral supervisor of the Department of Orthopedics, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College. LD is the professor and doctoral supervisor of the school of Biological Science and Medical Engineering, Beihang University.

Corresponding authors

Correspondence to Li Ding or Yu Zhao.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the institutional review board of Qinghai Red Cross Hospital.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qu, H., Yu, LJ., Wu, JT. et al. Spine system changes in soldiers after load carriage training in a plateau environment: a prediction model research. Military Med Res 7, 63 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: