# Experimental assessment of cervical ranges of motion and compensatory strategies

Céline Niewiadomski, Rohan-Jean Bianco, Sanae Afquir, Morgane Evin,
Pierre-Jean Arnoux

► To cite this version:

Céline Niewiadomski, Rohan-Jean Bianco, Sanae Afquir, Morgane Evin, Pierre-Jean Arnoux. Experimental assessment of cervical ranges of motion and compensatory strategies. Chiropractic & Manual Therapies, 2019, 27 (9), 9p. 10.1186/s12998-018-0223-x. hal-02937761

HAL Id: hal-02937761

https://hal.science/hal-02937761

Submitted on 16 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

CrossMark

# Experimental assessment of cervical ranges of motion and compensatory strategies

Céline Niewiadomski¹,<fnref n="3" />†, Rohan-Jean Bianco<fnref n="2" />,<fnref n="3" />†, Sanae Afquir<fnref n="2" />,<fnref n="3" />, Morgane Evin<fnref n="2" />,<fnref n="3" />* and Pierre-Jean Arnoux<fnref n="2" />,<fnref n="3" />

## Abstract

**Background:** Literature is still limited regarding reports of non-invasive assessment of the cervical range of motion in normal subjects. Investigations into compensatory motions, defined as the contribution of an additional direction to the required motion, are also limited.

The objectives of this work were to develop and assess a reliable method for measuring the cervical range of motion in order to investigate motion and compensatory strategies.

**Methods and data collection:** Ninety-seven no neck-related pain subjects (no severe cervical pathology, 57 women, age: 28.3 ± 7.5y. old, BMI: 22.5 ± 3.2 kg/m<fnref n="2" />) underwent a non-invasive cervical range of motion assessment protocol. In-vivo head's motion relative to the thorax was assessed through the measurement of the main angular amplitudes in the 3 directions (flexion/extension, axial rotations and lateral inclinations) and associated compensatory motions using an opto-electronic acquisition system.

**Results:** The principal motion reproducibility resulted in intra-class correlation coefficients ranging from 0.81 to 0.86. The following maximum ranges of motion were found: 127.4 ± 15.1° of flexion/extension, 89.3 ± 12° of lateral inclinations and 146.4 ± 13° of axial rotations after 6 outlier exclusions. Compensatory motions highly depend on the associated principal motion: for flexion/extension: (3.5 ± 7.6;-2.1 ± 7.8°), for rotation: (25.7 ± 17.9°;0.4 ± 4.7)°, for inclination: (22.9 ± 34.7°;-0.04 ± 8.7°). Age, BMI and weight significantly correlated with flexions (p < 0.032). Motion patterns were identified through clustering.

**Conclusions:** This kinematic analysis has been proven to be a reliable diagnostic tool for the cervical range of motion. The non-unicity and variability of motion patterns through the clustering of motion strategy identification have been shown. Compensatory motions contributed to such motion pattern definition despite presenting significant intra-individual variability.

**Keywords:** Rotation, Cervical motion, Compensatory motion, Motion pattern, Motion strategy, Experimental, Non-invasive motion measurements

## Introduction

Cervicalgia, perceived as incapacitating for patients, is a medical, social and economic burden on society [1]. This type of pathology has been reported to be associated with advanced age and smoking or exposure to tobacco in childhood [1]. The number of cervicalgia diagnoses is

expected to increase as a majority of people are likely to "experience some degree of neck pain in their lifetime" [2]. Neck pain diagnosis is required to accurately manage and address cervicalgia treatment. Indeed, cervical pain in cervicalgia could lead to a modification of the range of motion (ROM), as pain will prevent the subject from fully performing a specific task. The ROM of the cervical spine could then be assessed quantitatively and non-invasively and has been shown to be restricted in patients suffering from cervical radiculopathy and cervicogenic headaches [3].

The cervical segment is a complex articular part of the human body; it can be divided into two parts: the lower

<figure><img src="image_3.png" /></figure>

BMC

© The Author(s). 2019 **Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

<footnote>* Correspondence: morgane.evin@ifsttar.fr</footnote>

<footnote>†Céline Niewiadomski and Rohan-Jean Bianco contributed equally to this work.</footnote>

<footnote n="2">²iLab Spine - Laboratoire de Biomécanique Appliquée, UMRT24 IFSTTAR - Aix Marseille Université, Faculté de Médecine secteur-Nord, 51 Bd. P. Dramard, 13015 Marseille, France</footnote>

<footnote n="3">³iLab-Spine - Laboratoire international en imagerie et biomécanique du Rachis, Marseille, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
cervical part (C3-C7) and the upper one (C0-C2) [4, 5]. According to Watier et al. [6], the spine ROM is defined in each direction of the space and the ROM can be assessed for each functional unit. A wide range of protocols is available to assess the cervical spine ROM. These protocols assess the ROM of the head using: motion acquisition systems [7, 8], inclinometers [9, 10], goniometers [11] or radiography imaging [4]. Main ROM values were provided either by dividing the cervical units as in Frobin et al. [12], or by studying the global ROM. In the no neck-related pain population, the ranges of cervical motion were defined using different protocols. For flexion/extension, a maximum ROM of 87.2° to 145° was found; lateral inclination ranged from 59° to 186°, and rotation from 75° to 175° [6]. In 25 healthy volunteers, flexion was found to equal 68 ± 6.2°, extension 68.3 ± 7.3°, right flexion 49.8 ± 7.5°, left flexion 52.6 ± 7.6°, right rotation 78 ± 6.4° and left rotation 77 ± 7.7° [13]. In 13 healthy volunteers, similar results were found: 79.4 ± 11.7° in extension, 62.2 ± 11.1° in flexion, 46.6 ± 6° in left bending, 48.6 ± 6.9° in right bending, 69.8 ± 7.1° and 71.2 ± 6.4° in left and right rotations [14].

The resolution of the acquisition system and the error associated with data acquisition are one of the issues encountered in the available protocols for cervical ROM measurement. Thus, regarding the measurement of the angular error, very few angular error measurements were reported when others reported 0.1° [15] to 3.92° [16]. Another issue concerns the definition of the reference positioning of the subject before and during the acquisition. Walmsley et al. [17] showed that depending on the subject's head positioning in neutral position, the ROMs were modified. Such results were confirmed by Wang et al. [18] for flexion extension while contradicted by Sato et al. [19] in a small population.

The complex structure of the spine functional units results in the coupling between large ROMs resulting from combined uni-axial motions. Motions of the adjacent vertebrae are induced by articular surface geometry (different orientation along the spine level) and enable the coupling between axial homo-lateral rotation and lateral inclination motions [20]. At the lower cervical spine level and resulting from the zygapophyseal articulation orientation and from cervical lordosis, lateral inclination could be associated with a homo-lateral rotation and light extension [21]. The intra-individual measurements have shown combined ROM in simple motion, which highlights the variability of motion patterns for a similar instruction or final position [21]. Additionally, compensatory motions are defined as the contribution of an additional direction to the required motion while an alteration of the amplitude of the required motion could be noticed [22, 23]. At the level of the cervical spine structure, several ROM strategies

could be used to reach one given objective. Quantitative descriptions of compensatory motions, defined as patterns established unintentionally by the subject in case of cervical pains, are highly limited [24–26] depending on the ROM acquisition technique and the recruitment criteria.

In this context, this work aims 1) to develop a quantitative and reproducible method for the measurement of cervical spine mobility, 2) to measure main and compensatory cervical motions in normal subjects and 3) to characterize motion patterns in this population according to subject characteristics.

# Methods

## Data collection

This study was performed on 101 volunteers according to the inclusion criteria: older than 18 and younger than 65 years old, ability to agree to survey participation, presenting no severe neck pains. Exclusion criteria defined as reported in medical records included severe cervical pathologies, bone fractures or cervical pains, cervicobrachial neuralgia, cardiac or vascular pathologies, pulmonary, cephalic and visceral pathologies, pregnancy or heavy medical treatment. This experimental setting received the approval of the Agence Régionale de Santé for experimentation on healthy volunteers (N°2017–5).

## Experimental layout

Subjects were seated on a specific experimental chair (Fig. 1a). Hip flexion was set to a 90° angular positioning in using the adjustable foot support and a level to ensure horizontality of the thighs. Subjects were positioned on the chair with straps across their pelvis and thorax to avoid improper compensatory motions of the pelvic and scapular waists. The thoraco-lumbar spine was fixed with an adjustable lumbar support to promote physiological lordosis (pelvis set with a slightly retroversion of iliac muscles) and thus create a compensatory physiological dorsal kyphosis that enabled free cervical positioning. The backrest of the chair was in contact with the subject only at lumbar level in order to set the thoracic spine free. Subjects were asked to find a neutral position (horizontal eye trajectory set with a reference on the wall), hands on thighs, palms up to avoid support of the superior limbs and improper movements of the scapular waist. Method repeatability was assessed on 50 subjects selected randomly. These subjects underwent the experimental protocol twice within a fifteen-minute interval.

Subjects were asked to perform several types of motion, starting from a neutral position and returning to this neutral position between each instruction. They were asked to follow their own pace in order to ensure that they stay within their comfort range of motion. Each motion had to be performed with maximum
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Experimental setting. Positioning of a subject on the experimental chair with thorax ROM reduction straps during a complex motion instruction (a). Codamotion (Charnwood Dynamics Ltd.) acquisition system for opto-electronic motion. Aluminium support with 3 video cameras and location sensors (b). Headband and opto-electronic sensors (c)</figcaption></figure>

amplitude while staying within the comfort zone (non-forced motion and no pain) and avoiding thoracic motions. The first sequence of motions was composed of 5 flexions followed by 5 extensions and 30 s of rest in the neutral position. The second sequence included 5 axial rotations to the right followed by 5 axial rotations to the left and 30 s of rest in the neutral position. Third sequence was composed of 5 inclinations to the right followed by 5 inclinations to the left and 30 s of rest in the neutral position.

The motion acquisition system used to capture motions of the subjects was a Codamotion (Codamotion CX1, Charnwood Dynamics Ltd., UK) system. It is composed of 3 pre-calibrated video cameras in an aluminium support and of active opto-electronic markers (4 sensors per sensor box) located in front of the subject's seat and enabling maximum visibility of the markers (Fig. 2-ab). Marker displacements were tracked in the 3 directions of the space, using specific infrared frequencies. Two sensors were placed on the shoulders of the subject at the coracoid apophysis level. A rigid adjustable headband instrumented with 4 sensors (fixed in a T shape on the anterior part of the head band) was placed on the subject's head (Fig. 1-c). Acquisition was performed at a frequency of 100 Hz and spatial resolutions of 0.1 mm. Measurements were performed with the ODIN software (Charnwood Dynamics Ltd.). The acquisition system was

aligned to define the origin at the centre of the chair: Z axis defined as the caudo-cranial direction, X axis as the medio-lateral direction on the subject's left and Y axis in the antero-posterior direction. The relative locations of the 4 sensors on the headband were fixed during the experiment. Any missing data caused by the obstruction of sensor visualisation by the cameras was computed with a specific algorithm using the 3 known sensor locations to deduce the fourth one.

## Kinematic data analysis

Exported data was post processed using custom-made Matlab (The MathWorks, Inc.) scripts for this study.

For each motion, neutral initial positioning was used as a reference. The angular variations of the 4 sensor positions were measured in each plane projection (XY, XZ, YZ) and in the neutral position. The angular variations of the sensors for the flexion-extension motions (sagittal plan), axial rotation (transverse plan) and inclination (frontal), were measured for each acquired time step. For each motion, maximal amplitudes as well as standard deviations on 5 similar motions were computed to avoid false motions or effects related to subject weariness (Fig. 2).

Neutral positions were quantified for each motion sequence (flexion/extension, inclination and rotation) to identify the variability of the initial position between two
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Data analysis representation. Right/left rotations with planar projections (transvers. Sagittal and frontal) and 3D representation. Subject SE-S3: 24 years old, male. Compensatory motion. Positioning of a subject on experimental chair with thorax ROM reduction straps during a complex motion instruction (a). Right/left rotations with planar projections (transvers, sagittal and frontal) and (b) 3D representation. Subject SE-S3: 24 years ols, male. (c) Compensatory motions</figcaption></figure>

consecutive motions. Compensatory motions were defined as maximum amplitude on the two other planes at the same time as the associated principal motions: for flexion/extension motion, transverse and frontal compensatory motions were computed, for right and left rotations, sagittal and frontal compensatory motions and for right and left inclination, sagittal and transverse compensatory motions.

## Statistical methods

First, univariate statistical analysis was performed to insure data homogeneity. Principal component analysis (PCA) was performed on the principal motions analysis and used to identify outliers, taking into consideration the first five axes (presenting 80% of inertia). Method repeatability was tested on 50 subjects on account of time constraints and sufficient in terms of statistical power for such quantification. It was assessed using intra-class correlation coefficient (ICC – consistency ICC for one-way random effects model) for each principal motion. Clusters of the 97 subject population studied (without outliers) were identified according to the coordinates of the first five axes of the PCA of principal and reliable compensatory motions (ICC > 0.75) using agglomerative hierarchical clustering (AHC- Warel method). Additionally, ANOVA for the assessment of cluster group influence and linear regressions between morphological characteristics

(gender, age, height and weight) and global ROM measurements (flexion + extension, right/left rotation and right/left inclination) were tested. Statistical analyses were performed using R software. Significant results were considered with a *p* value < 0.05.

## Results

### Method reliability

The neutral position as selected as the minimum angular value between two consecutive motions was found to equal 3.0 ± 2.5°, with an intra-individual variability within the same motion of 6.8 ± 3.8° on average in flexion/extension, respectively 1.6 ± 2.5° and 5.3 ± 2.7° in inclination and 2.2 ± 2.2° and 4.2 ± 2.1° in rotation (Table 1).

Reproducibility between the test and retest on 50 subjects of the method was assessed using an intra-class coefficient. This resulted in an ICC of 0.81 for flexions/extensions, of 0.86 for right/left inclinations and of 0.83 for right/left rotations. Regarding compensatory motions, the ICC averaged 0.41 and 0.26 for flexions/extensions in transverse and frontal planes respectively, 0.81 and 0.75 for sagittal and transverse planes in inclinations and 0.89 and 0.38 for sagittal and frontal planes in rotations.

The missing data, which had to be computed by the dedicated algorithm using the 3 known sensor locations, was 9.7% on average.
<table><caption>**Table 1** Maximum amplitude of ROM principal motions. Intra-individual variability of the neutral position. Repeatability of the ROM measurement in 50 no neck-related pain volunteers. Annotation: values are presented as mean (standard deviation)</caption><thead><tr><th></th><th>Total</th><th>Neutral position</th><th>ICC (95%)</th></tr></thead><tbody><tr><td>Gender (Female)</td><td>101 (57)</td><td></td><td></td></tr><tr><td>Age (y.)</td><td>28.3 (7.5)</td><td></td><td></td></tr><tr><td>Weight (kg)</td><td>66.2 (12.4)</td><td></td><td></td></tr><tr><td>Height (cm)</td><td>171 (8.8)</td><td></td><td></td></tr><tr><td>BMI (kg/m2)</td><td>22.5 (3.2)</td><td></td><td></td></tr><tr><td>Flexion (°)</td><td>65.3 (12.2)</td><td>3 (2.5)</td><td>0.81 (0.68-0.88)</td></tr><tr><td>Extension (°)</td><td>-60.9 (9.2)</td><td>6.8 (3.8)</td><td>0.81 (0.69-0.89)</td></tr><tr><td>Right Rot. (°)</td><td>71.2 (8.1)</td><td>2.2 (2.2)</td><td>0.82 (0.7-0.89)</td></tr><tr><td>Left Rot. (°)</td><td>-74.5 (8.1)</td><td>4.2 (2.1)</td><td>0.83 (0.72-0.9)</td></tr><tr><td>Right Incl. (°)</td><td>43.3 (12.2)</td><td>5.3 (2.7)</td><td>0.92 (0.86-0.95)</td></tr><tr><td>Left Incl. (°)</td><td>-45.3 (9.2)</td><td>1.6 (2.5)</td><td>0.80 (0.67-0.88)</td></tr></tbody></table>

## Ranges of motion

The cohort included 42 males and 55 females (age: 28.3 years range 21–53, median: 25 years; weight: 66.4 kg range 41 to 106 kg; height: 171.3 cm, range 152 to 196 cm). Ranges of principal motions are presented in Table 1 for the entire population without outliers. Eighty subjects practiced a physical activity on a regular basis (no professional practice) and 49 had suffered cervical pain (related to car accidents or sports) (simple whiplash, physiological curve inversion, cervical sprain, associated cervicalgia). Maximum flexion/extension angular motions were found of 126.8 (15.7)°, maximum right and left inclinations equalled 88.8(12.4)° and axial rotations equalled 146.1(13.3)° (Table 3). Intra-individual standard

deviation was found 2.4(1.7)° for flexion/extension, 1.7(1.3)° for rotations and 1.4(1.7)° for inclinations.

## Cluster analysis and motion patterns

Analysis of the ascending hierarchy classification of the population enabled the definition of 3 clusters. The first cluster included 30 subjects more advanced in age, with a similar gender ratio to the entire population, a higher BMI, a lower average ROM compared to the entire population, and the highest reduction in flexion/extension (–12.3% on average for both) and inclinations (–11.2%, rotations: –10%) (Table 2).

In the second cluster, presenting an equal number of male subjects, slightly heavier and with a similar BMI, the ranges of motion were found to be higher than in the entire population (+6% for flexion/extension, +13.7% for inclinations and 5.3% for rotations). The third cluster of 35 subjects presented a similar average ROM to the entire population, more women, younger and with a lower BMI. Significant differences between clusters were found for all range of motion measurements (AOV: *p* < 0.003) except for inclinations.

Averaged times required to perform each task are reported in Table 2. Times for flexion/extension motions were found to have required slightly longer (383 ± 81 ms) than rotations (303 ± 70 ms) and inclinations (292 ± 69 ms). Average times did not differ significantly between clusters.

## Compensatory motions analysis

Compensatory analysis presented high inter-individual variability and values depending on the associated

<table><caption>**Table 2** Ranges of principal motion of painless ROM within the 3D space in 97 subjects. Reference values for cluster populations. Annotation: values are presented as mean (standard deviation)</caption><thead><tr><th></th><th>Cluster 1</th><th>Cluster 2</th><th>Cluster 3</th><th>Total</th><th>p</th></tr></thead><tbody><tr><td>Population (Female)</td><td>30 (16)</td><td>32 (15)</td><td>35 (24)</td><td>97 (55)</td><td></td></tr><tr><td>Age (y.)</td><td>31.8 ± 9.7</td><td>27.7 ± 6.8</td><td>26 ± 4.8</td><td>28.3 ± 7.5</td><td>0.002</td></tr><tr><td>Weight (kg)</td><td>69.9 ± 15.2</td><td>66.8 ± 11.1</td><td>63 ± 10.5</td><td>66.4 ± 12.5</td><td>0.2</td></tr><tr><td>Height (cm)</td><td>171.2 ± 9.5</td><td>172.8 ± 8</td><td>169.9 ± 8.7</td><td>171.3 ± 8.7</td><td>ns</td></tr><tr><td>BMI (kg/m2)</td><td>23.7 ± 4.1</td><td>22.3 ± 2.7</td><td>21.7 ± 2.6</td><td>22.5 ± 3.2</td><td>0.1</td></tr><tr><td>Flexion (°)</td><td>55.3 ± 9</td><td>71.2 ± 8.6</td><td>69.9 ± 10.1</td><td>65.8 ± 11.6</td><td>&lt; 0.001</td></tr><tr><td>Extension (°)</td><td>-55.7 ± 9.1</td><td>-63.3 ± 8.8</td><td>-63.3 ± 7.5</td><td>-61 ± 9.1</td><td>&lt; 0.001</td></tr><tr><td>Right Rot. (°)</td><td>64.6 ± 7.3</td><td>75.6 ± 5.2</td><td>74.1 ± 6</td><td>71.7 ± 7.8</td><td>&lt; 0.002</td></tr><tr><td>Left Rot. (°)</td><td>-66.9 ± 6.7</td><td>-78.3 ± 6.4</td><td>-77.6 ± 4.4</td><td>-74.5 ± 7.8</td><td>&lt; 0.003</td></tr><tr><td>Right Incl. (°)</td><td>38.4 ± 6.2</td><td>49.6 ± 6.4</td><td>42 ± 5.2</td><td>43.4 ± 7.5</td><td>ns</td></tr><tr><td>Left Incl. (°)</td><td>-40.4 ± 6.2</td><td>-51.4 ± 5.2</td><td>-44.2 ± 4.5</td><td>-45.4 ± 6.9</td><td>ns</td></tr><tr><td>Mean Time Flex/Ext</td><td>362 ± 89.2</td><td>391.1 ± 68.1</td><td>393.6 ± 83.4</td><td>383 ± 81</td><td>ns</td></tr><tr><td>Mean Time Rot</td><td>287.4 ± 64.4</td><td>314.4 ± 65.4</td><td>305.9 ± 76.5</td><td>303 ± 69.5</td><td>ns</td></tr><tr><td>Mean Time Incl.</td><td>282.9 ± 65.4</td><td>295.9 ± 62.1</td><td>296.7 ± 77.5</td><td>292.2 ± 68.6</td><td>ns</td></tr></tbody></table>
principal motions (Table 3). For flexion/extension, the amplitudes of compensatory motions were similar between transverse and frontal planes (maximum 11°, -0.4 ± 4.7 and 3.4 ± 7.6 in the entire population), while main compensatory motions were found in sagittal planes for both rotation (maximum 33.7, 22.7 ± 33.9 in the entire population vs. frontal: vs. maximum 1.8, -2 ± 8 in the entire population) and inclination (transverse: maximum 36.8, 25.7 ± 20.6 in the entire population vs. sagittal: maximum -1.9, 0.4 ± 9.1 in the entire population).

## Correlations with subject characteristics

Correlation analysis between subject characteristics (gender, weight, height, BMI and age) and principal ROMs showed significant a correlation between weight and flexion (r = -0.19, p = 0.03). BMI significantly correlated with flexion (r = -0.29, p = 0.003). Age correlates with all principal motions (p > 0.15, Table 4).

## Discussion

### Reliability of the methods

Results on the reliability of the principal motion measurement method, assessed during a test and retest at a fifteen-minute interval, showed good reproducibility (ICC ≥80). The reproducibility of a subject's response to an instruction will be influenced by the reproducibility of the neutral position which has been assessed here. It was shown to be lower than 6.8° maximum in extension (ranging between 1.6 and 5.3 for others motions). Additionally, intra-individual standard deviation was limited (found maximum in flexion/extension at 2.4° on average) and far under the neutral position showing higher reproducibility of maximum amplitude value. The reproducibility of compensatory motion analysis was irregular depending on the associated principal motion. Compensatory motions associated with inclinations and those within the sagittal plane associated with the rotation motions could be considered as reproducible (ICC > 0.75) while others should be

used carefully. Such results could be explained by motion mechanisms (inclination is more likely to induce a higher compensatory motion than flexion/extension; similarly, rotation is more likely to induce a compensatory motion in the sagittal plane).

### ROM results

The angular values of the main motions in this cohort were similar to the ones reported in the literature by several authors [6]: flexion and extension: 87.2° to 145° (our result on average: 126.8°), lateral inclination on average: 59° to 186° (our result: 88.8°) and rotation 75° to 175° (our result on average: 146.1°). Depending on the range of age, weight, height and gender of the participants, results vary from one study to another. When similar age ranges to our cohort were considered, similar results were found [27].

The average time required to perform the task was not found to differ between clusters in our result, however the differences in terms of age and BMI, while significant, remained limited. Such differences may not be substantial enough to show changes in terms of average time previously reported in the literature [28].

### Compensatory motions analysis

As previously highlighted in the literature [20, 21], coupling between cervical spine motions has been shown in the lateral inclination motion. Coupling between homo-lateral rotations was highlighted in the lower cervical vertebrae (C3-C7), taking into account sub-occipital and cranial rotations in opposite direction to the inclination [29]. Such compensatory motions were reported to show a functional outcome of the head and spine and were explained by the anatomy of the inclination and articular facet. Quantifications of such compensatory motions were limited [30] and our results show the feasibility of such quantifications and their inputs (cluster analysis) in the cervical motion analysis despite the high

<table><caption>**Table 3** Compensatory motion amplitudes and principal motions in 97 subjects. Reference values for cluster populations. Annotation: values are presented as mean (standard deviation)</caption><thead><tr><th></th><th>Cluster 1</th><th>Cluster 2</th><th>Cluster 3</th><th>Total</th><th>p</th></tr></thead><tbody><tr><td>Population (Female)</td><td>30 (16)</td><td>32 (15)</td><td>35 (24)</td><td>97 (55)</td><td></td></tr><tr><td>Flexion/extension (°)</td><td>111 ± 13.6</td><td>134.5 ± 8.4</td><td>133.2 ± 12.6</td><td>126.8 ± 15.7</td><td>&lt; 0.001</td></tr><tr><td>Rotation (°)</td><td>131.4 ± 10.9</td><td>153.8 ± 8.1</td><td>151.7 ± 7.9</td><td>146.1 ± 13.3</td><td>&lt; 0.001</td></tr><tr><td>Inclination (°)</td><td>78.9 ± 9.3</td><td>101 ± 9.3</td><td>86.2 ± 7</td><td>88.8 ± 12.4</td><td>0.04</td></tr><tr><td>Trans Flexion/Extension (°)</td><td>-0.7 ± 4.7</td><td>0.8 ± 5.2</td><td>0.9 ± 4.3</td><td>0.4 ± 4.7</td><td>ns</td></tr><tr><td>Front Flexion/Extension (°)</td><td>5 ± 8.1</td><td>1.5 ± 8.2</td><td>3.9 ± 6.2</td><td>3.4 ± 7.6</td><td>ns</td></tr><tr><td>Sag. Right/Left Rot. (°)</td><td>13.8 ± 15.4</td><td>28.6 ± 38.6</td><td>24.9 ± 39.8</td><td>22.7 ± 33.9</td><td>ns</td></tr><tr><td>Front. Right/Left Rot. (°)</td><td>2.6 ± 6</td><td>-0.9 ± 7.4</td><td>-6.9 ± 7.4</td><td>-2 ± 8</td><td>&lt; 0.001</td></tr><tr><td>Sag. Right/Left Incl. (°)</td><td>0 ± 8.6</td><td>5.4 ± 9.1</td><td>-3.7 ± 7.4</td><td>0.4 ± 9.1</td><td>0.07</td></tr><tr><td>Trans. Right/Left Incl. (°)</td><td>17.9 ± 15.5</td><td>26.2 ± 26</td><td>31.9 ± 16.7</td><td>25.7 ± 20.6</td><td>0.006</td></tr></tbody></table>
<table><caption>**Table 4** Linear regression between subject characteristics and principal motions: adjusted R value (p value): *** depicted *p* < 0.001. ** depicted *p* < 0.01 and * depicted *p* > 0.05</caption><thead><tr><th>r(p)</th><th>Gender (M/F)</th><th>Weight (kg)</th><th>Height (cm)</th><th>BMI (kg/m2)</th><th>Age (years)</th></tr></thead><tbody><tr><td>Flexion/extension</td><td>0.07 (0.24)</td><td>-0.19 (0.032)*</td><td>0.1 (0.81)</td><td>-0.29 (0.003)**</td><td>-0.35(&gt; 0.001)***</td></tr><tr><td>Rotation</td><td>0.09 (0.45)</td><td>-0.07 (0.23)</td><td>0.1 (0.96)</td><td>-0.13 (0.102)</td><td>-0.26 (0.006)**</td></tr><tr><td>Inclination</td><td>-0.05 (0.14)</td><td>0.09 (0.18)</td><td>0.03 (0.35)</td><td>0.06 (0.247)</td><td>-0.23 (0.015)*</td></tr></tbody></table>

inter-individual variability. Such inter-individual variability could then be explained by both lifestyle and low significance dysfunction.

High values were found in the sagittal planes for rotations and inclinations. For right and left axial rotations, all subjects presented large variations in the sagittal plane. For inclinations, compensatory motions in the sagittal plane were also found to be necessary in order to perform the required motion. Flexion/extension presented the more balanced measurement in terms of compensatory motions.

## Cluster analysis

To perform the required instruction, the subject will use different cervical mobility patterns (principal motions and compensatory motions). This cluster analysis revealed significant differences in the compensatory ROM, while principal motion differences were found to differ significantly between clusters. This highlights the non-unique motion pattern to reach a similar displacement. Thus, frontal plane contributions to principal rotation and transverse plane contributions to inclination motions were highlighted. Several results could explain differing motion patterns within the population: anatomical variability (skeletal), muscle strength and activation during the task, as well as physical activity and lifestyle. While the time required to perform a task has been proven to influence motion patterns [31], no significant differences in average times were reported in our results.

## Correlation with subject characteristics

The correlation analysis showed that all range of motion will reduce with age, which was already highlighted previously in the literature [30, 32]. The effects of aging tested here are limited as far as the correlation coefficient between motions and age are concerned, probably due to the limited number of available volunteers. However, as highlighted in the first section of the discussion, our results were found to be similar to the literature whenever the age of the cohort was similar [13, 14]. This could suggest that our cohort age would require to be widened in order to validate the correlation between age and principal motion. Such influence has been explained in the literature through the physiological degenerative mechanism of the spine increasing disk and posterior articulation stiffness, decreasing ROM on the transverse

and sagittal planes [33, 34]. Such interpretation highlights the main role of intervertebral disks in axial rotation and flexion/extension as well as the significance of their alterations.

Additionally, the slight correlation between weight and flexion could be explained by the increase in the muscular and adipose mass of the visceral area and of the neck. Similar results were found in Malmström et al. [30] regarding weight, height and age while showing that older, overweight subjects had decreased overall ROM. In contrast, slightly younger females with lower weight compared to the group average had higher or similar ROM (cluster analysis: cluster 1 vs. cluster 3).

## Limitations

Experimental acquisition systems as well as protocol and data analysis have limitations, which partially explain the variability of the measurements found in the literature. In this study, error related to protocol was intended to be limited by initial subject positioning (head and thorax postures) as well as by the codaMotion acquisition method, which enabled initialisation of the subject posture to a reference similar in all subjects while remaining specific to each subject. Despite these efforts, by increasing the number of sensors on the cranium, face, sternum and clavicles, we would improve the reliability and precision of the results and would be adding data to characterize ROM. Inter-days and inter-operator reproducibility for marker positioning should also have been tested. Similarly, a second acquisition system could help in measuring front and back motions. Additional data analysis, including helicoidally defined reference system analysis, has been developed but could however not be tested with the present set of data [35]. Medical questionnaires or indexes could have been used to further describe the selected population as in [25]. Symmetrical analysis could also be performed, requiring the laterality of the subjects. Finally, despite the limited number of volunteers we were able to recruit, results show a correlation between ranges of motion, BMI and weight. A larger cohort could however be considered to further confirm these results and enable a closer investigation into the correlation with age. Further analysis is needed to fully assess the benefits of compensatory motion analysis regarding cervicalgia treatment.
## Conclusion

This work has provided description and testing of a robust and reliable protocol in order to measure main cervical ROM through head to torso motion detection. Compensatory motion measurements as well as the reliability of the method have been quantified and the principal motion amplitude could thus be defined as a normal reference value for cervical ROM in the specific age range of our cohort. Cluster analysis highlighted the need to assess compensatory motions and revealed the non-unicity of the motion pattern when following a similar instruction in no neck-related pain volunteers. Additionally, age, BMI and weight were shown to correlate with mostly flexion-related ranges of motion. Further work is needed to understand how such quantitative methods could objectively contribute to diagnostic and treatment.

## Acknowledgements

We would like to thank Max Py for his technical assistance.

## Funding

There is no funding to report.

## Availability of data and materials

Data is available on request to the writer.

## Authors' contributions

CN and RB acquired the data and set the protocol. RB and ME post-processed the data. SA and ME statistically analysed the results. PJA designed the study and all authors contributed to the writing of the manuscript.

## Ethics approval and consent to participate

All volunteers who took part in this study gave their full written consent on the day of the experiment and have received the information on the protocol. This experimental setting received the approval of the Agence Régionale de Santé for experiments on healthy volunteers (N°2017-5).

## Consent for publication

All the authors consent to the writing of this manuscript. Written informed consent was obtained from the participants for publication of their individual details and accompanying images in this manuscript. The consent form is held by the authors and is available for review by the Editor-in-chief.

## Competing interests

The authors declare that they have no competing interest.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

### Author details

¹Collège d'Ostéopathie de Provence, Aix en Provence, France. ²iLab Spine - Laboratoire de Biomécanique Appliquée, UMRT24 IFSTTAR - Aix Marseille Université, Faculté de Médecine secteur-Nord, 51 Bd. P. Dramard, 13015 Marseille, France. ³iLab-Spine - Laboratoire international en imagerie et biomécanique du Rachis, Marseille, France.

Received: 14 July 2018 Accepted: 13 December 2018

Published online: 22 January 2019

## References

1. Hogg-Johnson S, van der Velde G, Carroll LJ, Holm LW, Cassidy JD, Guzman J, et al. The burden and determinants of neck pain in the general population: results of the bone and joint decade 2000-2010 task force on neck pain and its associated disorders. Spine. 2008;33:S39-51.

2. Haldeman S, Carroll LJ, Cassidy JD. The empowerment of people with neck pain: introduction: the bone and joint decade 2000-2010 task force on neck pain and its associated disorders. Spine. 2008;33:S38.

3. Snodgrass SJ, Cleland JA, Haskins R, Rivett DA. The clinical utility of cervical range of motion in diagnosis, prognosis, and evaluating the effects of manipulation: a systematic review. Physiotherapy. 2014;100:290-304.

4. Rousseau M-A, Laporte S, Chavary-Bernier E, Lazennec J-Y, Skalli W. Reproducibility of measuring the shape and three-dimensional position of cervical vertebrae in upright position using the EOS stereoradiography system. Spine. 2007;32:2569-72.

5. Chicheportiche V. Anatomie et imagerie du rachis cervical supérieur et de la charnière cranio-vertébrale. Rev Rhum Monogr. 2013;80:26-31.

6. Watier B. Comportement mécanique du rachis cervical: une revue de littérature. ITBM-RBM. 2006;27:92-106.

7. Amiri M, Jull G, Bullock-Saxton J. Measuring range of active cervical rotation in a position of full head flexion using the 3D Fastrak measurement system: an intra-tester reliability study. Man Ther. 2003;8:176-9.

8. Bechara BP, Bell KM, Hartman RA, Lee JY, Kang JD, Donaldson WF. In vivo analysis of cervical range of motion after 4- and 5-level subaxial cervical spine fusion. Spine. 2012;37:E23-9.

9. Audette I, Dumas J-P, Côté JN, De Serres SJ. Validity and between-day reliability of the cervical range of motion (CROM) device. J Orthop Sports Phys Ther. 2010;40:318-23.

10. Yoo W-G, Park S-Y, Lee M-R. Relationship between active cervical range of motion and flexion-relaxation ratio in asymptomatic computer workers. J Physiol Anthropol. 2011;30:203-7.

11. Hostler D, Colburn D, Seitz SR. A comparison of three cervical immobilization devices. Prehospital Emerg Care. 2009;13:256-60.

12. Frobin W, Leivseth G, Biggemann M, Brinckmann P. Sagittal plane segmental motion of the cervical spine. A new precision measurement protocol and normal motion data of healthy adults. Clin Biomech Bristol Avon. 2002;17:21-31.

13. Chiu TTW, Sing KL. Evaluation of cervical range of motion and isometric neck muscle strength: reliability and validity. Clin Rehabil. 2002;16:851-8.

14. Guo L-Y, Lee S-Y, Lin C-F, Yang C-H, Hou Y-Y, Wu W-L, et al. Three-dimensional characteristics of neck movements in subjects with mechanical neck disorder. J Back Musculoskelet Rehabil. 2012;25:47-53.

15. Feipel V, Rondelet B, Le Pallec J, Rooze M. Normal global motion of the cervical spine: an electrogoniometric study. Clin Biomech Bristol Avon. 1999;14:462-70.

16. Petersen CM, Johnson RD, Schuit D. Reliability of cervical range of motion using the OSI CA 6000 spine motion analyser on asymptomatic and symptomatic subjects. Man Ther. 2000;5:82-8.

17. Walmsley RP, Kimber P, Culham E. The effect of initial head position on active cervical axial rotation range of motion in two age populations. Spine. 1996;21:2435-42.

18. Wang X, Lindstroem R, Carstens NPB, Graven-Nielsen T. Cervical spine reposition errors after cervical flexion and extension. BMC Musculoskelet Disord. 2017;18:102.

19. Sato T, Masui K. Effect of initial head position on range of active cervical flexion. J Phys Ther Sci. 2012;24:1141-4.

20. Boussion L, Chêze L. Étude cinématique tridimensionnelle du rachis cervical. Datareves177901230011109-11036 [Internet]. 2011 [cited 2017 Sep 5]; Available from: http://www.em-consulte.com/en/article/280320

21. Harrison DE, Harrison DD, Cailliet R, Janik TJ, Troyanovich SJ. Cervical coupling during lateral head translations creates an S-configuration. Clin Biomech Bristol Avon. 2000;15:436-40.

22. Lysell E. Motion in the cervical spine. An experimental study on autopsy specimens. Acta Orthop Scand. 1969;Suppl 123:1+.

23. Berthoz A, Graf W, Vidal PP. The head-neck sensory motor system: Oxford University Press; 1992. http://www.oxfordscholarship.com/view/10.1093/ acprof:oso/9780195068207.001.0001/acprof-9780195068207.

24. Ernst MJ, Crawford RJ, Schelldorfer S, Rausch-Osthoff A-K, Barbero M, Kool J, et al. Extension and flexion in the upper cervical spine in neck pain patients. Man Ther. 2015;20:547-52.

25. Rudolfsson T, Björklund M, Djupsjöbacka M. Range of motion in the upper and lower cervical spine in people with chronic neck pain. Man Ther. 2012;17:53-9.
26. Stenneberg MS, Busstra H, Eskes M, van Trijffel E, Cattrysse E, Scholten-Peeters GGM, et al. Concurrent validity and interrater reliability of a new smartphone application to assess 3D active cervical range of motion in patients with neck pain. Musculoskelet Sci Pract. 2018;34:59–65.

27. Ferrario VF, Sforza C, Serrao G, Grassi G, Mossi E. Active range of motion of the head and cervical spine: a three-dimensional investigation in healthy young adults. J Orthop Res. 2002;20:122–9.

28. Zawadka M, Skublewska-Paszkowska M, Gawda P, Lukasik E, Smolka J, Jablonski M. What factors can affect lumbopelvic flexion-extension motion in the sagittal plane?: a literature review. Hum Mov Sci. 2018;58:205–18.

29. Ishii T, Mukai Y, Hosono N, Sakaura H, Nakajima Y, Sato Y, et al. Kinematics of the upper cervical spine in rotation: in vivo three-dimensional analysis. Spine. 2004;29:E139–44.

30. Malmström E-M, Karlberg M, Fransson PA, Melander A, Magnusson M. Primary and coupled cervical movements: the effect of age, gender, and body mass index. A 3-dimensional movement analysis of a population without symptoms of neck disorders. Spine. 2006;31:E44–50.

31. Srinivasan D, Mathiassen SE, Samani A, Madeleine P. The combined influence of task accuracy and pace on motor variability in a standardised repetitive precision task. Ergonomics. 2015;58:1388–97.

32. Dvorak J, Antinnes JA, Panjabi M, Loustalot D, Bonomo M. Age and gender related normal motion of the cervical spine. Spine. 1992;17:S393–8.

33. Trott PH, Pearcy MJ, Ruston SA, Fulton I, Brien C. Three-dimensional analysis of active cervical motion: the effect of age and gender. Clin Biomech Bristol Avon. 1996;11:201–6.

34. Park MS, Moon S-H, Lee H-M, Kim T-H, Oh JK, Nam JH, et al. Age-related changes in cervical sagittal range of motion and alignment. Glob Spine J. 2014;4:151–6.

35. Ellingson AM, Yelisetti V, Schulz CA, Bronfort G, Downing J, Keefe DF, et al. Instantaneous helical Axis methodology to identify aberrant neck motion. Clin Biomech Bristol Avon. 2013;28:731–5.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

BMC