- Original Paper
Fractal properties of postural sway during quiet stance with changed visual and proprioceptive inputs
The Journal of Physiological Sciences volume 61, pages 123–130 (2011)
We studied the fractal dynamics of postural sway during quiet stance with changed visual and/or proprioceptive information. Radius-vector length and angle of rotation of center of foot pressure (COP) displacements were used as parameters of postural sway. The experiments were performed on 24 healthy volunteers of both sexes aged 20–30 years. Using a battery of nonlinearity tests, differences in the degree of fractality of both parameters during stance were found. We found that the behavior of radius-vector was similar to fractional Brownian motion, while that of angle of rotation was similar to flicker noise. Quantitative parameters that can be used to characterize the changes in radius-vector length and angle of rotation include self-similarity intervals and fractal dimension. In healthy subjects, the process of postural control maintains its fractal structure independently of altered sensory information. We believe that this analysis provides information about new methods for evaluating postural sway behavior during quiet stance.
The maintenance of quiet upright stance involves complex mechanisms of integrating information from the visual, vestibular, and proprioceptive systems. The evaluation of postural sway measurements has been the subject of intense research in recent decades. Measurements of body balance are mostly based on displacements of the center-of-foot-pressure (COP) in the plane under the feet. Stabilograms describe the COP coordinates along the x-axis or medial-lateral (ML), and along the y-axis or anterior-posterior (AP), as time series: ML = f(t) and AP = f(t). The curve showing how the COP shifts in the Cartesian coordinate system [y = f(x)] is called a statokinesigram. The COP position in the statokinesigram is characterized by the radius-vector length (RD) and angle of rotation (β) between two consecutive COP points [1–4].
There are two main approaches for evaluating postural steadiness. The first one evaluates the COP displacements by standard statistical analysis of basic parameters such as mean amplitude, mean velocity, length of sway path, sway area, and frequency parameters of the power spectrum [5–7]. The second approach is based on recent nonlinear techniques for analysis of the postural control systems [8–12]. These studies have focused on the nonlinear properties of postural sway that are not evident in the classic basic descriptive parameters. The stabilogram diffusion analysis (SDA) and the detrended fluctuation analysis (DFA) are the most used nonlinear dynamic analyses of postural control [11, 12].
SDA is based on the theory of pure fractional Brownian motion [8, 9, 13]. Collins and De Luca  estimated the stochastic properties of the COP dynamics by the Hurst exponent. They hypothesized that there are two stochastic processes in the control of posture during quiet stance with different scaling exponents: a short-term process (t < 1 s) with an open-loop control mechanism and a long-term process with a closed-loop mechanism.
Delignieres et al.  argued that the SDA method was inappropriate for the pure fractional Brownian motion and proposed the DFA analysis to assess the fractal characteristics of COP. DFA is a method of integrating and detrending the data series; it allows for the scaling exponent to be calculated. Amoud et al. , however, reported differences between the postural sway of elderly and young persons using SDA and DFA and found that both analyses are relevant. The nonlinear properties of postural sways have also been analyzed using fractal dimensions calculated by means of Higuchi’s algorithm [15–18].
All these findings suggest that postural stability may be described by an exponential decreased law of 1/f α. Thus, it would be of great interest to investigate how this law will be modified when the visual and/or somatosensory inputs are altered and whether the components of decomposed postural sways follow the same law.
The aim of the present study was to determine the nonlinear properties of body sways measured in quiet stance under normal conditions and under conditions of sensory conflicted stance when the visual and proprioceptive inputs are changed. The postural sway was characterized by the radius-vector length and angle of rotation of COP displacements. We also estimated the separate contribution of these parameters in terms of the fractal dimension (D f) and structure functions (SF).
Experimental setup and data recording
Twenty-four healthy subjects, 10 males and 14 females, aged between 20 and 30 years (mean ± SD, 24.2 ± 3.6) without a history of vestibular, musculoskeletal, or neurological disorders were investigated. All participants were volunteers and gave their informed consent prior to taking part in the experiments, which were approved by the Bioethics Committee of the Institute of Neurobiology, BAS, Sofia.
During the experimental sessions the subjects stood freely upright on the platform with their heels separated by a distance of 2.5 cm, feet positioned with a 25° angle between their medial sides, and with arms freely hanging alongside the body. The quiet stance task was performed either on a stable support or on a foam pad with dimensions of 400 × 400 × 150 mm, density of 77.5 kg/m3, and elastic modulus of 74.9 N/m2.
Two visual input conditions were used: eyes open (gaze fixed on a point 2 m in front of the subject and at eye level) and eyes closed. Each trial lasted 30 s.
Postural sway was measured by a static posturographic system (Fig. 1) [19, 20]. A position-sensitive photodetector (PSD) fixed on the bottom of the platform was used to record the COP displacements in two orthogonal directions in the plane under the feet. A laser source of light was installed on top of the platform above the PSD so that it moved in all directions proportionally to the difference in pressure between the two feet. The action of the PSD, sensitive to visible light, was based on the lengthways photoeffect. The nonlinearity of the position curve of PSD in both directions was 1.3%. The photopotentials generated by PSD when the body was upright were proportional to the COP displacements in both X and Y directions. The sensitivity of the measurement of the platform was 0.1 mm/bit. The signals were digitized with sampling interval of 10 ms (sampling frequency F s = 100 Hz).
A digital Hamming low-pass filter with a cut-off frequency of 10 Hz was used to remove high frequency noise.
The radius-vector length and its angle were calculated according to the following equations:
where N is the number of samples, n the transitory value of postural sway, and AP and ML are the anterior-posterior and medial-lateral planes, respectively. Raw samples of statokinesigrams of one person are shown in Fig. 2.
The first step in the nonlinear analysis is to check the type of distribution as normal distribution and stationarity of data are required for conventional linear analyses. There are three possibilities in cases where the distribution is non-Gaussian: (1) the process is linear but non-Gaussian; (2) the process has linear dynamics, but the observations are obtained by nonlinear “static” transformation, and (3) the process has nonlinear dynamics .
In the present study we used the Kolmogorov-Smirnov test to check the amplitude distributions of RD and β. Three tests were applied to reveal the nonlinear dynamic structure of COP displacements during the quiet upright stance. The first one evaluated the self-similarity of the process of postural control by calculating the duration of self-similarity intervals and fractal dimension (D f). The power spectrum (PS) of postural sways exhibited an exponential decreasing shape (1/f α) when plotted on the double logarithmic scale. The self-similarity intervals included the linear part of the slope. The fractal dimension was calculated from the linear slope α of the log-log plot of the PS. According to Higuchi , there is a direct relation between α and D f. The fractal dimensions have to be calculated as follows: when 1 < α < 3, D f = (5 − α)/2; when α ≥ 3, D f = (7 − α)/2; and when α ≤ 1, D f = (3 − α)/2. In case of a fractional Brownian motion process with α = 2, D f is equal to 1.5.
The second test distinguishes the fractal from the chaotic nature of the signals through the inspection of the structure function (SF) shape presented in the log-log plot for the raw data epochs and its first derivative . SF was determined according to the formula:
where t = 1 to N − n and y(t) are the samples of the scalar signal.
According to Provenzale et al. , when the process is chaotic, the SF(n) increases for small n and saturates for large n for both the original signal and its first derivative. In contrast, the fractal noisy processes are characterized by increased SF(n) at the large scale for the real signal and saturated for the first derivative of the signal.
The third test is based on the surrogate data method for detection of possible chaotic dynamics of postural sways. The idea is to generate many realizations from a single real data recording by preserving its linear properties and by removing some nonlinear properties that may exist. The important characteristic of the chaotic dynamics is the correlation dimension D 2 computed by the algorithm of Grasberger-Procaccia .
In this study we used the amplitude-adjusted Fourier transform algorithm (AAFT)  to generate 10 surrogate data epochs. This algorithm removes the possible nonlinear properties of the real data epochs by randomizing the phases after the fast Fourier transformation (FFT), while preserving all of their linear properties. The null hypothesis was that the real time series was generated after the monotonic nonlinear transformation of a linear Gaussian process. D 2 was computed for values of embedding dimension m = 2, 4,…, 18. For reconstruction in the phase space, time delay τ = 2 was chosen according to the zero-crossing of the autocorrelation function. The correlation dimension of surrogate data (D 2 sur) was computed by the modified algorithm (re-embedding) procedure with m = 10 and m = 20 . The significance of the differences between D 2 for real data (D 2 real) and D 2 sur was statistically estimated.
All statistical analyses were performed using the statistical package Statistica 7.0 (Stat Soft, USA, 2004). The significances of the differences in mean values were assessed with Wilcoxon’s nonparametric rank test. Data are presented as mean value ± SD, and p ≤ 0.05 was considered significant.
By testing the null hypothesis using the maximum distance D between real and Gaussian cumulative distributions, we revealed that the amplitude distributions of RD and β for all experimental conditions were asymmetrical and non-Gaussian according to the Kolmogorov-Smirnov test. The null hypothesis should be rejected if D > D a, where D a is the critical value for a chosen level of significance. For all experimental conditions, D was higher than D a, and the null hypothesis was rejected with p < 0.05. This result indicates possible nonlinear dynamics of RD and β.
The power spectra of RD and β plotted on a double logarithmic scale consisted of two segments with different slopes (Fig. 3). The first segment started in the very low frequency area and was characterized by an almost flat curve with a zero slope, similar to the random white noise. The second segment had a linear decreasing curve, which was characteristic of a process with long-range correlation with fractal properties.
The slope (α) of the second segment was determined by fitting a linear polynomial regression function and was used to compute D f of the selected segment (Fig. 3). The slope is limited by two characteristic frequencies. Their reciprocal values determined the domain of the interval between trial numbers with a long-range correlation of RD and β, where the self-similarity was exhibited (Tables 1, 2).
For all conditions, the slopes of RD (αRD) were in the range of 1 < α < 3 (mean 2.59 ± 0.4) and D f was equal to 1.2 ± 0.19, indicating fractional Brownian processes within the selected frequency bands . Significant differences in the mean values of α RD across the four experimental conditions were not observed (Table 1). The initial frequency of intervals for RD self-similarity during the stance on a foam support with closed eyes of all subjects was higher (0.72 Hz) than measured in the other three experimental conditions (Table 1). The foam pad stance led to significant decreases in the duration of self-similarity intervals compared with stable surface stance with eyes closed (p < 0.01). The lack of vision led to a significant decrease in this parameter during the foam pad stance (p < 0.01). Mean duration of the self-similarity intervals for the other experimental conditions was about 5 s, and significant differences were not observed. There were also significant differences of the D f mean values between the stance on stable and foam support, both with open eyes (p < 0.05), and closed eyes (p < 0.01) (Table 1).
The slopes of angle of rotation of the COP displacements (αβ) were <1 (mean 0.8 ± 0.24) and D f = 1.09 ± 0.06, which is typical for a flicker-like noise process with blurred self-similarity (Table 2). The mean value of αβ during stance on foam support with closed eyes was significantly higher (p < 0.05) than in the other experimental conditions (Table 2). The initial frequencies of self-similarity intervals of β were shifted to higher frequencies of the spectrum when standing on the foam support with both open and closed eyes compared to when standing on the hard support. The durations of self-similarity intervals of the angles in subjects with both open and closed eyes when standing on the foam pad were significantly smaller (p < 0.05) than during the stable surface stance (Table 2). The β self-similarity interval was between 2 s for the stable support and 0.4 s for the foam support (Table 2). Significant differences in Df were observed between subjects on stable and foam supports (p < 0.001) with eyes closed and between subjects with open and closed eyes standing on the foam support (p < 0.001).
Figure 4 illustrates the double log plots of SF and their first derivatives for RD (Fig. 4a) and β (Fig. 4b). SF for RD increased for small n and saturated for large n, while the first derivatives of SF were almost constant (Fig. 4a). According to the criteria of Provenzale et al. , these changes in SF and its first derivative suggest a nonlinear stochastic process with fractal dynamics. The shapes of both SF and its first derivative for β were approximately identical with flat plateaus with zero slopes (Fig. 4b). These results suggested possible chaotic dynamics of β .
Some classes of the stationary random processes with 1/f α spectrum exhibit nonlinear fractal properties that could be erroneously classified as chaotic dynamics. The RD showed typical fractal dynamics, while β presented ambiguous behavior (fractal or chaotic). To clarify these findings, an additional test was applied to test the surrogate data.
In order to select the proper embedding dimension, the correlation integrals and their first derivatives (D 2 ) were computed for m = 2, 4,…, 18. A saturation of the plot D 2 for real data (D 2 real) was reached after m = 14 (Fig. 5). The null hypothesis Ho (D 2 real = D 2 sur) was not rejected (Wilcoxon’s test, p < 0.05). In other words, chaotic dynamics of β was not observed for all experimental conditions (Table 3). Rather, β had a fractal structure with a blurred self-similarity.
The results of this study support the hypothesis of fractal behavior of postural sway during quiet upright stance promoted by some authors, who have used various nonlinear methods and posturographic measurements for study [11–16]. We described a simple technique for estimating the nonlinear properties of two parameters (radius-vector length and angle of rotation) of COP displacements in support under the feet by an estimation of duration of the self-similarity intervals, fractal dimension, and structural function. Our results showed that in healthy subjects during quiet stance, RD and β have nonlinear properties that can still be described by the 1/f α.law, regardless of changes due to the altered visual and/or proprioceptive inputs. The length of the linear slope α in the log-log plot of the PS, the duration of the self-similarity intervals, and the fractal dimension were different for RD and β, which defined their different levels of the fractality. The RD structure shows typical fractional Brownian motion (α about 2 and D f between 1 and 1.5), while β showed a blurred self-similarity identical to the flicker noise (α < 1 and D f close to 1).
In order to maintain postural stability, the central nervous system integrates information from multiple sensory channels and generates complex motor responses. Two major mechanisms are responsible for postural adjustment: (1) compensatory postural adjustment (feedback mechanisms), activated by sensory events related to posture unsteadiness and (2) anticipatory postural adjustment (feedforward mechanisms), predicting the disturbances and producing preprogrammed responses for the maintenance of stability. Some of the compensatory postural adjustments are innate, while others are learned during growth and daily life. The minor role of the visual system in the fractal behavior of postural sway has been observed in healthy adults during stance on stable support [14–16]. In these studies, the examination of COP trajectories was based on the evaluation of medial-lateral and anterior-posterior postural sways. Furthermore, no significant differences between subjects with open and closed eyes standing on hard support were established using a variety of nonlinear analysis techniques including detrended fluctuation analysis and diffusion coefficients , fractal dimension by modified pixel dilation (mPD) method , and Higuchi’s algorithm . Our results are in line with the above suggestion as we also did not observe significant differences in our investigated parameters—duration of the self-similarity intervals and D f—between subjects with eyes open and closed during hard support stance. We assume that the central nervous system (CNS) has saved in memory a relatively stable model of postural control during quiet upright stance on firm support as this is the most widely encountered environment in daily life; the role of the visual system is less important here than proprioceptive and vestibular inputs [16, 26].
Collins and De Luca [8, 9] suggested that both closed-loop and open-loop mechanisms of postural control are present to control postural sway. Open-loop control and closed-loop control are closely related to the concepts of feedforward and feedback, respectively. All of these concepts are derived from systems and control engineering, and they provide descriptions of two different ways in which the nervous system behaves in controlling movement. Feedback control is often referred to as closed-loop control because the outgoing commands, the effectors, the feedback signals, and the control center all form a loop. Feedforward control is often referred to as open-loop control in order to emphasize that feedback sensory signals do not directly affect the timing of the response. In general, closed-loop control can supplement open-loop control in maintaining upright stance stability, since feedback can be utilized by the CNS to control postural sway when there are challenges to the central program. The long duration of self-similarity intervals of RD (about 5 s) and β (about 2 s) can be explained by this hypothesis [8, 9] for postural stability during stance on firm support.
The ability to assess pressure distribution and body orientation decreases when standing on a foam pad, where proprioception and tactile information are altered . In addition, the compliant visco-elastic surface reduces the effectiveness of the ankle torque required for the postural stabilization . These facts can explain why, during the upright stance on the foam pad, we observed significant differences in the mean duration of self-similarity intervals and mean values of fractal dimension D f under both conditions—open and closed eyes.
In case of a heavy sensory conflict, when postural control was relying on the vestibular system only (stance on foam support with closed eyes), the self-similarity intervals were shorter and the D f decreased (by about 1) for both the linear and angular parameters of the COP displacements. This indicates an enhanced contribution of the open-loop mechanism related to feedforward postural control. When the vestibular input was predominantly active, the CNS produced an adequate response to correct for postural instability. It is likely that the influence of vestibular information is essential to accomplish the necessary changes in postural control strategy from ankle to hip during stance on unstable support . The changed information from the visual and proprioceptive inputs in healthy subjects, however, impacts only the degree of fractality, while the total fractal structure of postural control process is preserved.
The visual system participates in the maintenance of postural stability, especially when other sensory systems are compromised [29, 30]. Our results during stance on the foam support with eyes open also support the suggestion that vision provides an effect of stabilization in the case of altered proprioceptive and tactile information.
In conclusion, the present results indicate that the COP displacements during quiet stance have complex dynamics and nonlinear characteristics. The self-similarity interval and fractal dimension of both radius-vector length and angle of rotation of COP are quantitative parameters related to the dynamic characteristics of the process-generating program of postural control and can be used for investigation of postural stability.
Bulbulian R, Hargan ML (2000) The effect of activity history and current activity on static and dynamic postural balance in older adults. Physiol Behav 70:319–325
Fujita T, Nakamura S, Ohue M, Fujii Y, Miyauchi A, Takagi Y, Tsugeno H (2005) Effect of age on body sway assessed by computerized posturography. J Bone Miner Metab 23:152–156
Karst GM, Venema DM, Roehrs TG, Tyler AE (2005) Center of pressure measures during standing tasks in minimally impaired persons with multiple sclerosis. J Neurol Phys Ther 29:170–180
Gerbino P, Griffin E, Zurakowski D (2007) Comparison of standing balance between female collegiate dancers and soccer players. Gait Posture 26:501–507
Dichgans J, Mauritz KH, Allum JHJ, Brandt T (1976) Postural sway in normals and atactic patients: analysis of the stabilizing and destabilizing effects of vision. Agressologie 17:15–24
Hufschmidt A, Dichgans J, Mauritz KH, Hufschmidt M (1980) Some methods and parameters of body sway quantification and their neurological application. Arch Psychiat Nervenkl 228:135–150
Prieto TR, Hoffmann R, Lovett E, Myklebust M (1996) Measures of postural steadiness: differences between healthy young and elderly adults. IEEE Trans Biomed Eng 43:956–966
Collins JJ, De Luca CJ (1993) Open-loop and closed-loop control of posture: a random-walk analysis of center-of-pressure trajectories. Exp Brain Res 95:308–318
Collins JJ, De Luca CJ (1995) Upright, correlated random walks: a statistical-biomechanics approach to the human postural control system. Chaos 5:57–63
Laughton CA, Slavin M, Katdare K, Nolan L, Bean JF, Kerrigan DC, Phillips E, Lipsitz LA, Collins JJ (2003) Aging, muscle activity, and balance control: physiologic changes associated with balance impairment. Gait Posture 18:101–108
Duarte M, Zatsiorsky VM (2000) On the fractal properties of natural human standing. Neurosci Lett 283:173–176
Delignieres D, Deschamps T, Legros A, Caillou N (2003) A methodological note on nonlinear time series analysis: is the open- and closed-loop model of Collins and De Luca (1993) a statistical artifact? J Mot Behav 35:86–97
Rougier R (1999) Influence of visual feedback on successive control mechanisms in upright quiet stance in humans assessed by fractional Brownian motion modelling. Neurosci Lett 266:157–160
Amoud H, Abadi M, Hewson D, Pellegrino VM, Doussot M, Duchene J (2007) Fractal time series analysis of postural stability in elderly and control subjects. J Neuroeng Rehabil 4:1–12
Manabe Y, Honda E, Shiro Y, Sakai K, Kohira I, Kashihara K, Shohmori T, Abe K (2001) Fractal dimension analysis of static stabilometry in Parkinson’s disease and spinocerebellar ataxia. Neurol Res 23:397–404
Doyle T, Dugan E, Humphrie B, Newton R (2004) Discriminating between elderly and young using a fractal dimension analysis of centre of pressure. Int J Med Sci 1:11–20
Higuchi T (1988) Approach to an irregular time series on the basis of the fractal theory. Phys D 31:277–283
Higuchi T (1990) Relationship between fractal dimension and the power low index for a time series. Phys D 46:254–264
Ratcheva T, Stambolieva K, Kostadinov K (1999) Posturograph with position-sensitive detector and method for its preparing. BG Patent #61749 http://v3.espacenet.com/publicationDetails/originalDocument?CC=BG&NR=61749B1&KC=B1&FT=D&date=19980529&DB=EPODOC&locale=en_EP
Stambolieva K, Racheva T, Kostadinov K (1998) Posturographic system with position-sensitive detector of registration. In: Proceedings of 7th international conference on electronics, Sofia, vol 3, pp 45–50
Popivanov D, Mineva A (1999) Testing procedures for non-stationarity and non-linearity in physiological signals. Math Biosci 157:303–320
Provenzale A, Smith LA, Vio R, Murante G (1992) Distinguishing between low-dimensional dynamics and randomness in measured time series. Phys D 58:31–49
Theiler J (1986) Spurious dimensions from correlation algorithms applied to limited time series data. Phys Rev A 34:2427–2433
Theiler J, Eubank S, Longin A, Galdrikian B, Farmer JD (1992) Testing for nonlinearity in time series: the method of surrogate data. Phys D 58:77–94
Fraedrich K, Wang R (1993) Estimating the correlation dimension of an attractor from noisy and small data sets based on re-embedding. Phys D 65:373–398
Tinetti ME, Williams TF, Mayewski R (1986) Fall risk index for elderly patients based on number of chronic disabilities. Am J Med 80:429–434
Perry S, McIllroy W, Maki B (2000) The role of plantar cutaneous mechanoreceptors in the control of compensatory stepping reactions evoked by unpredictable, multi-directional perturbation. Brain Res 877:401–406
MacLellan MJ, Patla AE (2006) Adaptations of walking pattern on a compliant surface to regulate dynamic stability. Exp Brain Res 173:521–530
Jeka J, Kiemel T, Creath R, Horak F, Peterka R (2004) Controlling human upright posture: velocity information is more accurate than position or acceleration. J Neurophysiol 92:2368–2379
Rosengren KS, Rajendran K, Contakos J, Chuang L, Peterson M, Doyle R, McAuley E (2007) Changing control strategies during standard assessment using computerized dynamic posturography with older women. Gait Posture 25:215–221
Shumway-Cook A, Horak F (1986) Assessing the influence of sensory interaction on balance. Phys Ther 66:1548–1559
The study was partly supported by the National Fund for Scientific Research, Grant TK_02/60. The author is grateful to Prof. Slavianka Moyanova, Prof. Lilia Christova, and Prof. Plamen Gatev for their fruitful comments.
Rights and permissions
About this article
Cite this article
Stambolieva, K. Fractal properties of postural sway during quiet stance with changed visual and proprioceptive inputs. J Physiol Sci 61, 123–130 (2011). https://doi.org/10.1007/s12576-010-0129-4
- Body sway
- Postural control
- Nonlinear analysis
- Fractal dynamics