Barry Shender

ABSTRACT Military pilots are subjected to high magnitude inertial loads applied to the head-neck complex during high-G maneuvers. Cervical spinal soft-tissue injuries have occurred in this population [1–3]. Acute injury rates were reported between 54 and 89%, most commonly resulting in muscle or neck pain. Early cervical spine degenerative changes were also identified for fighter pilots [4]. Because the neck muscles are responsible for maintaining head-neck stability, one study hypothesized that cervical injuries in aviators may result from insufficient neck muscle strength to support the head-neck complex during high-G maneuvers [5]. This hypothesis is supported by the finding that pilots participating in pre-injury neck strengthening exercises demonstrated fewer injuries [1]. Although clinical data on the subject are limited, female pilots may be more susceptible to neck injury due to more slender necks and cervical columns that may be less resistant to bending [6, 7]. Differences in neck muscle geometry, in terms of cross-sectional area and positioning, may also lead to differing injury rates. Previous investigations of neck muscle geometry using contemporary medical imaging modalities were conducted with subjects in supine position [8–11], which removes the axial loads of the head and superior cervical structures due to gravity and likely changes neck muscle geometry. To date, no study has outlined gender-dependent neck muscle geometry determined using MRI of subjects in upright, sitting posture. The present hypothesis was that significant gender differences exist in neck muscle geometry.

A uniaxial tensile loading study of 13 lumbar porcine ligaments under varying environmental temperature conditions. To investigate a possible temperature dependence of the material behavior of porcine lumbar anterior longitudinal ligaments. Temperature dependence of the mechanical material properties of ligament has not been conclusively established. The anterior longitudinal ligaments (ALLs) from domestic pigs (n = 5) were loaded in tension to 20% strain using a protocol that included fast ramp/hold and sinusoidal tests. These ligaments were tested at temperatures of 37.8 degrees C, 29.4 degrees C, 21.1 degrees C, 12.8 degrees C, and 4.4 degrees C. The temperatures were controlled to within 0.6 degrees C, and ligament hydration was maintained with a humidifier inside the test chamber and by spraying 0.9% saline onto the ligament. A viscoelastic model was used to characterize the force response of the ligaments. The testing indicated that the ALL has strong temperature dependence. As temperature decreased, the peak forces increased for similar input peak strains and strain rates. The relaxation of the ligaments was similar at each temperature and showed only weak temperature dependence. Predicted behavior using the viscoelastic model compared well with the actual data (R2 values ranging from 0.89 to 0.99). A regression analysis performed on the viscoelastic model coefficients confirmed that relaxation coefficients were only weakly temperature dependent while the instantaneous elastic function coefficients were strongly temperature dependent. The experiment demonstrated that the viscoelastic mechanical response of the porcine ligament is dependent on the temperature at which it is tested; the force response of the ligament increased as the temperature decreased. This conclusion also applies to human ligaments owing to material and structural similarity. This result settles a controversy on the temperature dependence of ligament in the available literature. The ligament viscoelastic model shows a significant temperature dependence on the material properties; instantaneous elastic force was clearly temperature dependent while the relaxation response was only weakly temperature dependent. This result suggests that temperature dependence should be considered when testing ligaments and developing material models for in vivo force response, and further suggests that previously published material property values derived from room temperature testing may not adequately represent in vivo response. These findings have clinical relevance in the increased susceptibility of ligamentous injury in the cold and in assessing the mechanical behavior of cold extremities and extremities with limited vascular perfusion such as those of the elderly.

Experimental testing incorporating lumbar columns and isolated components is essential to advance the understanding of injury tolerance and for the development of safety enhancements. This study incorporated a whole column axial acceleration model and an isolated vertebral body model to quantify compression rates during realistic loading and compressive tolerance of vertebrae. Eight lumbar columns and 53 vertebral bodies from 23 PMHS were used. Three-factor ANOVA was used to determine significant differences (p<0.05) in physiologic and failure biomechanics based on compression rate, spinal level, and gender. Results demonstrated a significant increase in ultimate force (i.e., fracture) from lower to higher compression rates. Ultimate stress also increased with compression rate. Displacement and strain to failure were consistent at both compression rates. Differences in ultimate mechanics between vertebral bodies obtained from males and females demonstrated non-significant trends, with female vertebral bodies having lower ultimate force that would be associated with decreased injury tolerance. This was likely a result of smaller vertebrae in that population. Combined with existing literature, results presented in this manuscript contribute to the understanding of lumbar spine tolerance during axial loading events that occur in both military and civilian environments with regard to effects of compression rate and gender.

Clinical studies have indicated that thoracolumbar trauma occurs in the civilian population at its junction. In contrast, injury patterns in military populations indicate a shift to the inferior vertebral levels of the lumbar spine. Controlled studies offering an explanation for such migrations and the associated clinical biomechanics are sparse in literature. The goals of this study were to investigate the potential roles of acceleration loading on the production of injuries and their stability characteristics using a human cadaver model for applications to high-speed aircraft ejection and helicopter crashes. Biomechanical laboratory study using unembalmed human cadaver lumbar spinal columns. Thoracolumbar columns from post-mortem human surrogates were procured, x-rays taken, intervertebral joints and bony components evaluated for degeneration, and fixed using polymethylmethacrylate. The inferior end was attached to a platform via a load cell and uniaxial accelerometer. The superior end was attached to the upper metal platform via a semi-circular cylinder. The pre-flexed specimen was preloaded to simulate torso mass. The ends of the platform were connected to the vertical post of a custom-designed drop tower. The specimen was dropped inducing acceleration loading to the column. Axial force and acceleration data were gathered at high sampling rates, filtered, and peak accelerations and inertia-compensated axial forces were obtained during the loading phase. Computed tomography images were used to identify and classify injuries using the three-column concept (stable vs. unstable trauma). The mean age, total body mass, and stature of the five healthy degeneration-free specimens were 42 years, 73 kg, and 167 cm. The first two specimens subjected to peak accelerations of approximately 200 m/sec(2) were classified as belonging to high-speed aircraft ejection-type and the other three specimens subjected to greater amplitudes (347-549 m/sec(2)) were classified as belonging to helicopter crash-type loadings. Peak axial forces for all specimens ranged from 4.8 to 7.2 kN. Ejection-type loaded specimens sustained single-level injuries to the L1 vertebra; one injury was stable and the other was unstable. Helicopter crash-type loaded specimens sustained injuries at inferior levels, including bilateral facet dislocation at L4-L5 and L2-L4 compression fractures, and all specimens were considered unstable at least at one spinal level. These findings suggest that the severity of spinal injuries increase with increasing acceleration levels and, more importantly, injuries shift inferiorly from the thoracolumbar junction to lower lumbar levels. Acknowledging that the geometry and load carrying capacity of vertebral bodies increase in the lower lumbar spine, involvement of inferior levels in trauma sparing the superior segments at greater acceleration inputs agree with military literature of caudal shift in injured levels. The present study offers an experimental explanation for the clinically observed caudal migration of spinal trauma in military populations as applied to high-speed aircraft ejection catapult and helicopter crashes.

This study determined bone mineral density (BMD) of cervical, thoracic, and lumbar vertebrae in healthy asymptomatic human subjects. To test the hypothesis that BMD of neck vertebrae (C2-C7) is equivalent to BMD of lumbar vertebrae (L2-L4). BMD of lumbar vertebrae is correlated to their strength. Although numerous studies exist quantifying BMD of the human lumbar spine, such information for the cervical spine is extremely limited. In addition, BMD correlations are not established between the two regions of the spinal column. Adult healthy human female volunteers with ages ranging from 18 to 40 years underwent quantitative computed tomography (CT) scanning of the neck and back. All BMD data were statistically analyzed using paired nonrepeating measures ANOVA techniques. Significance was assigned at a P < 0.05. Linear regression analyses were used to compare BMD as a function of level and region; +/-95% confidence intervals were determined. When data were grouped by cervical (C2-C7), thoracic (T1), and lumbar (L2-L4) spines, mean BMD was 260.8 +/- 42.5, 206.9 +/- 33.5, and 179.7 +/- 23.4 mg/mL. Average BMD of cervical vertebrae was higher than (P < 0.0001) thoracic and lumbar spines. Correlations between BMD and level indicated the lowest r value for T1 (0.42); in general, the association was the strongest in the lumbar spine (r = 0.89-0.95). The cervical spine also responded with good correlations among cervical vertebrae (r ranging from 0.66 to 0.87). The present study failed to support the hypothesis that BMD of lumbar spine vertebrae is equivalent to its cranial counterparts. The lack of differences in BMD among the three lumbar vertebral bodies confirms the appropriateness of using L2, L3, or L4 in clinical or biomechanical situations. However, significant differences were found among different regions of the vertebral column, with the cervical spine demonstrating higher trabecular densities than the thoracic and lumbar spines. In addition, the present study found statistically significant variations in densities even among neck vertebrae.

In contrast to clinical studies wherein loading magnitudes are indeterminate, experiments permit controlled and quantifiable moment applications, record kinematics in multiple planes, and allow derivation of moment-angulation corridors. Axial and coronal moment-angulation corridors were determined at every level of the subaxial cervical spine, expressed as logarithmic functions, and level-specificity of range of motion and neutral zones were evaluated. segmental primary axial and coupled coronal motions do not vary by level. Although it is known that cervical spine responses are coupled, segment-specific corridors of axial and coronal kinematics under axial twisting moments from healthy normal spines are not reported. Ten human cadaver columns (23-44 years, mean: 34 +/- 6.8) were fixed at the ends and targets were inserted to each vertebra to obtain kinematics in axial and coronal planes. The columns were subjected to pure axial twisting moments. Range of motion and neutral zone for primary-axial and coupled-coronal rotation components were determined at each spinal level. Data were analyzed using factorial analysis of variance. Moment-rotation angulations were expressed using logarithmic functions, and mean +/-1 standard deviation corridors were derived at each level for both components. Moment-angulations responses were nonlinear. Each segmental curve for both components was well represented by a logarithmic function (r2 > 0.95). Factorial analysis of variance indicated that the biomechanical metrics are spinal level-specific (P < 0.05). Axial and coronal angulations of cervical spinal columns show statistically different level-specific responses. The presentation of moment-angulation corridors for both metrics forms a dataset for the normal population. These segment-specific nonlinear corridors may help clinicians assess dysfunction or instability. These data will assist mathematical models of the spine in improved validation and lead to efficacious design of stabilizing systems.

Survival analyses of a large cohort of published lumbar spine compression fatigue tests. To produce the first large-scale evaluation of human lumbar spine tolerance to repetitive compressive loading and to evaluate and improve guidelines for human exposure to whole-body vibration and repeated mechanical shock environments. Several studies have examined the effects of compressive cyclic loading on the lumbar spine. However, no previous effort has coalesced these studies and produced an injury risk analysis with an expanded sample size. Guidelines have been developed for exposure limits to repetitive loading (e.g., ISO 2631-5), but there has been no large-scale verification of the standard against experimental data. Survival analyses were performed using the results of 77 male and 28 female cadaveric spinal segment fatigue tests from 6 previously published studies. Segments were fixed at each end and exposed to axial cyclic compression. The effects of the number of cycles, load amplitude, sex, and age were examined through the use of survival analyses. Number of cycles, load amplitude, sex, and age all are significant factors in the likelihood of bony failure in the spinal column. Using a modification of the risk prediction parameter from ISO 2631-5, an injury risk model was developed, which relates risk of vertebral failure to repeated compressive loading. The model predicts lifetime risks less than 7% for industrial machinery exposure from axial compression alone. There was a 38% risk for a high-speed planing craft operator, consistent with epidemiological evidence. A spinal fatigue model which predicts the risk of in vitro lumbar spinal failure within a narrow confidence interval has been developed. Age and sex were found to have significant effects on fatigue strength, with sex differences extending beyond those accounted for by endplate area disparities.

Clinical literature consistently identifies women as more susceptible to trauma-related neck pain, commonly resulting from soft tissue cervical spine injury. Structural gender differences may explain altered response to dynamic loading in women leading to increased soft tissue distortion and greater injury susceptibility. Identify anatomic gender differences in cervical spinal geometry that contribute to decreased column stability in women. Previous studies investigating male and female vertebral and vertebral body geometry demonstrated female vertebral dimensions were smaller. However, populations were not size matched and parameters related to biomechanical stability were not reported. Computed tomography scans of the cervical spine were obtained from size-matched young healthy volunteers. Geometrical dimensions were obtained at the C4 level and analysis of variance determined significant gender differences. Two volunteer subsets were size matched based on sitting height and head circumference. All geometrical measures were greater in men for both subsets. Vertebral width and disc-facet depth were significantly greater in men. Additionally, segmental support area, combining interfacet width and disc-facet depth, was greater in men, indicating more stable intervertebral coupling. Present results of decreased linear and areal cervical dimensions leading to decreased column stability may partially explain increased traumatic injury rates in women.

The failure responses of the anterior longitudinal ligament, posterior longitudinal ligament, and ligamentum flavum were examined in vitro under large strain-rate mechanical loading. To quantify the failure properties for 3 cervical spinal ligaments at strain rates associated with traumatic events. There exists little experimentation literature for fast-rate loading of the cervical spine ligaments. The small amount of available information is framed only in extensive experimental coordinates, and not in the context of strains. Bone-ligament-bone complexes were strained at fast rates, in an incrementally increasing loading protocol using a servohydraulic mechanical test frame. Failure loads and displacements were converted to engineering and true stress and strain values, and compared for the different ligaments (anterior longitudinal ligament, posterior longitudinal ligament, and ligamentum flavum), spinal levels (C3-C4, C5-C6, and C7-T1), and for male versus female specimens. There were no significant differences in force or true stress for gender or spinal level. There was a significant difference in force and true stress for ligament type. A difference was found between the posterior longitudinal ligament and ligamentum flavum for failure force, and between the ligamentum flavum and both the anterior and posterior longitudinal ligaments for failure true stress. No significant differences were found in true strain for ligament, gender, or spinal level. The mean ligament failure true strain was 0.81. Failure true strains were approximately 57% of the failure engineering strains. Once the injury mechanisms of the cervical spine are fully understood, computational models can be employed to understand the potentially traumatic effects of clinical procedures, and mitigate injury in impact, falls, and other high-rate scenarios. The soft tissue failure properties in this study can be used to develop failure tolerances in fast-rate loading scenarios. Failure properties of the anterior and posterior longitudinal ligaments were similar, and the same properties can be used to model both ligaments.

Aging, trauma, or degeneration can affect intervertebral kinematics. While in vivo studies can determine motions, moments are not easily quantified. Previous in vitro studies on the cervical spine have largely used specimens from older individuals with varying levels of degeneration and have shown that moment-rotation responses under lateral bending do not vary significantly by spinal level. The objective of the present in vitro biomechanical study was, therefore, to determine the coronal and axial moment-rotation responses of degeneration-free, normal, intact human cadaveric cervicothoracic spinal columns under the lateral bending mode. Nine human cadaveric cervical columns from C2 to T1 were fixed at both ends. The donors had ranged from twenty-three to forty-four years old (mean, thirty-four years) at the time of death. Retroreflective targets were inserted into each vertebra to obtain rotational kinematics in the coronal and axial planes. The specimens were subjected to pure lateral bending moment with use of established techniques. The range-of-motion and neutral zone metrics for the coronal and axial rotation components were determined at each level of the spinal column and were evaluated statistically. Statistical analysis indicated that the two metrics were level-dependent (p < 0.05). Coronal motions were significantly greater (p < 0.05) than axial motions. Moment-rotation responses were nonlinear for both coronal and axial rotation components under lateral bending moments. Each segmental curve for both rotation components was well represented by a logarithmic function (R(2) > 0.95). Range-of-motion metrics compared favorably with those of in vivo investigations. Coronal and axial motions of degeneration-free cervical spinal columns under lateral bending showed substantially different level-dependent responses. The presentation of moment-rotation corridors for both metrics forms a normative dataset for the degeneration-free cervical spines.

Experimental studies indicate age and degeneration affect spinal biomechanics. In vitro biomechanical experimentation is used to validate finite element cervical spine models. A high percentage of experimental studies have utilized older specimens. Computer models based on these experimental studies may not accurately represent the normal population. Younger full-column and C5-C6 motion segments were tested under pure sagittal plane moments. A review of literature was conducted, and results from previous studies were compared to present data to determine whether age was an influencing factor in spinal biomechanics. Findings indicate younger specimens under equivalent pure moment loading magnitudes underwent greater ranges of motion between 0.5 and 2.5 Nm. Based on these preliminary findings, validation of finite element modeling to ensure biofidelity should consider age as a factor that may affect biomechanics.

The objective of this study was to determine the bone mineral density (BMD) of cervical vertebrae and correlate with the lumbar spine. Fifty-seven young adult healthy male volunteers, ranging from 18 to 41 years of age, underwent quantitative computed tomography (QCT) scanning of C2-T1 and L2-L4 vertebrae. To account for correlations, repeated measures techniques were used to compare data as a function of spinal level and region. Linear regression methods were used (+/-95% CI) to compare data as a function of spinal level and region. The mean age and body height were 25.0 +/- 5.8 years and 181.0 +/- 7.6 cm. BMD decreased from the rostral to caudal direction along the spinal column. Grouped data indicated that the neck is the densest followed by the first thoracic vertebra and low back with mean BMD of 256.0 +/- 48.1, 194.3 +/- 44.2, and 172.2 +/- 28.4 mg/cm(3), respectively; differences were statistically significant. While BMD did not vary significantly between the three lumbar bodies, neck vertebrae demonstrated significant trends. The matrix of correlation coefficients between BMD and spinal level indicated that the relationship is strong in the lumbar (r = 0.92-0.96) and cervical (r = 0.73-0.92) spines. Data from the present study show that the trabecular bony architecture of the neck is significantly different from the low back. These quantitative BMD data from a controlled young adult healthy human male volunteer population may be valuable in establishing normative data specifically for the neck. From a trabecular bone density perspective, these results indicate that lumbar vertebrae cannot act as the best surrogates for neck vertebrae. Significant variations in densities among neck vertebrae, unlike the low back counterpart, may underscore the need to treat these bones as different structures.

Military aviators are susceptible to spinal injuries during high-speed ejection scenarios. These injuries commonly arise as a result of strains induced by extreme flexion or compression of the spinal column. This study characterizes the vertebral motion of two postmortem human surrogates (PMHS) during a simulated catapult phase of ejection on a horizontal decelerator sled. During testing, the PMHS were restrained supinely to a mock ejection seat and subjected to a horizontal deceleration profile directed along the local z-axis. Two midsized males (175.3 cm, 77.1 kg; 185.4 cm, 72.6 kg) were tested. High-rate motion capture equipment was used to measure the three-dimensional displacement of the head, vertebrae, and pelvis during the ejection event. The two PMHS showed generally similar kinematic motion. Head injury criterion (HIC) results were well below injury threshold levels for both specimens. The specimens both showed compression of the spine, with a reduction in length of 23.9 mm and 45.7 mm. Post-test autopsies revealed fractures in the C5, T1, and L1 vertebrae. This paper provides an analysis of spinal motion during an aircraft ejection.The injuries observed in the test subjects were consistent with those seen in epidemiological studies. Future studies should examine the effects of gender, muscle tensing, out-of-position (of head from neutral position) occupants, and external forces (e.g., windblast) on spinal kinematics during aircraft ejection.

This study presents the results of seven aerospace manikin and three post mortem human surrogate (PMHS) horizontal deceleration sled tests. The objective of this study was to establish a body of baseline data that examines the ability of small (fifth percentile) manikins to predict whole-body kinematics associated with aircraft ejection, and whether currently available head and neck injury criteria are applicable in these situations. Subjects were exposed to a short-duration local z-axis sled pulse while horizontally seated and restrained in an ejection seat. Test subjects included instrumented fifth percentile female and male manikins, and two small (163.8 cm, 48.3 kg; 143.5 cm, 48.6 kg) female and one small (166.2 cm, 54.3 kg) male PMHS. The anterior (local x-axis) translations of the PMHS heads were less than those observed in the manikin tests, but the local z-axis translations of the PMHS heads were greater than those of the manikins. Z-axis translations of the manikins' T1 were generally similar to those of the PMHS T1, but the anterior x-axis translations of T1 were greater in the PMHS. The neck injury criterion (Nij) tended to under-predict observed injury (primarily ruptures of the posterior ligaments at C4-5, T2-3), and the Beam Criterion (BC) tended to over-predict observed injury for small occupants. The USN/USAF neck injury criteria (NIC) performed best in predicting the observed injuries. Present manikin designs do not predict the kinematics of PMHS in ejection tests. Further refinement of existing injury criteria is required to accurately predict location and severity of ejection-induced injuries.

Vertebral fracture is the most common severe injury during high-speed pilot ejection. However, the loading paradigm experienced by pilots may also lead to soft-tissue spinal injuries that are more difficult to quantify and can lead to long-term deficits. This manuscript describes a new experimental protocol to simulate the effects of pilot ejection on the tissues of the head-neck complex. The model permits precise control of head-neck complex initial positioning, detailed analysis of head and spinal kinematics and upper and lower neck loads, and the ability to thoroughly investigate and identify soft-tissue injuries through upper and lower neck injury criteria, radiography, manual palpation, and cryomicrotomy. For the current test, peak acceleration of +14.8 Gz was similar to actual ejection events and duration of the acceleration pulse was approximately 100 ms. The specimen was oriented in flexion prior to initiation of inferior-to-superiorly directed acceleration. Subfailure ligamentum flavum injuries were sustained at the C4-C5 and C5-C6 cervical spinal levels and identified by increased segmental motions during the simulated ejection, increased laxity following testing, and cryomicrotomy. Upper and lower neck injury criteria did not predict these soft-tissue injuries. This experimental model can be used for detailed analysis of the effects of gender, head-neck orientation, helmet instrumentation, and acceleration pulse characteristics on cervical spine injury potential during pilot ejection events.

Neck muscles are important in the static and dynamic stability of the head-neck complex. Deep neck muscles act to maintain upright posture and superficial muscles are responsible for gross movements. Previous studies have quantified neck muscle geometry using traditional supine magnetic resonance imaging (MRI). However, supine orientation removes the vertical load on the cervical spine from the head-neck complex and changes the relative orientation of the spine and neck muscles. Therefore, the purpose of this study was to demonstrate the feasibility of upright MRI to obtain neck muscle morphometric data on a spinal level-by-level basis for subjects in upright seated positions. Upright MRI scans were obtained of the neck region for six younger male volunteers in neutral and flexed positions. Planar images were oriented parallel to the intervertebral disc space at each level. Cross-sectional area (CSA) and orientation of neck muscles were quantified at four spinal levels. Area and position of all four muscles were significantly dependent upon spinal level. Average CSA of the sternocleidomastoid, longus colli, levator scapulae, and trapezius muscles in neutral position were 512, 113, 281, and 174 mm2. Head-neck position significantly affected area and position of the sternocleidomastoid and position of posterior neck muscles. Comparison of neck muscle areas from the present study to a previous study incorporating supine MRI demonstrated differing trends between anterior and posterior neck muscles that may be attributable to upright orientation of volunteers and planar image orientation in the present study. Differences between supine and upright MRI identified in the present study may warrant incorporation of this technique in future spinal imaging studies.