Dynamic head-neck stabilization and modulation with perturbation bandwidth investigated using a multisegment neuromuscular model

The human head-neck system requires continuous stabilization in the presence of gravity and trunk motion. We investigated contributions of the vestibulocollic reflex (VCR), the cervicocollic reflex (CCR), and neck muscle co-contraction to head-in-space and head-on-trunk stabilization, and investigated modulation of the stabilization strategy with the frequency content of trunk perturbations and the presence of visual feedback.We developed a multisegment cervical spine model where reflex gains (VCR and CCR) and neck muscle co-contraction were estimated by fitting the model to the response of young healthy subjects, seated and exposed to anterior-posterior trunk motion, with frequency content from 0.3 up to 1, 2, 4 and 8 Hz, with and without visual feedback.The VCR contributed to head-in-space stabilization with a strong reduction of head rotation (<8 Hz) and a moderate reduction of head translation (>1 Hz). The CCR contributed to head-on-trunk stabilization with a reduction of head rotation and head translation relative to the trunk (<2 Hz). The CCR also proved essential to stabilize the individual intervertebral joints and prevent neck buckling. Co-contraction was estimated to be of minor relevance. Control strategies employed during low bandwidth perturbations most effectively reduced head rotation and head relative displacement up to 3 Hz while control strategies employed during high bandwidth perturbations reduced head global translation between 1 and 4 Hz. This indicates a shift from minimizing head-on-trunk rotation and translation during low bandwidth perturbations to minimizing head-in-space translation during high bandwidth perturbations. Presence of visual feedback had limited effects suggesting increased usage of vestibular feedback.     

Analysis of control strategies for VIVA OpenHBM with active reflexive neck muscles

Modeling muscle activity in the neck muscles of a finite element (FE) human body model can be based on two biological reflex systems. One approach is to approximate the Vestibulocollic reflex (VCR) function, which maintains the head orientation relative to a fixed reference in space. The second system tries to maintain the head posture relative to the torso, similar to the Cervicocolic reflex (CCR). Strategies to combine these two neck muscle controller approaches in a single head-neck FE model were tested, optimized, and compared to rear-impact volunteer data. The first approach, Combined-Control, assumed that both controllers simultaneously controlled all neck muscle activations. In the second approach, Distributed-Control, one controller was used to regulate activation of the superficial muscles while a different controller acted on deep neck muscles. The results showed that any muscle controller that combined the two approaches was less effective than only using one of VCR- or CCR-based systems on its own. A passive model had the best objective rating for cervical spine kinematics, but the addition of a single active controller provided the best response for both head and cervical spine kinematics. The present study demonstrates the difficulty in completely capturing representative head and cervical spine responses to rear-impact loading and identified a controller capturing the VCR reflex as the best candidate to investigate whiplash injury mechanisms through FE modeling.

     

Self-organization of trajectory formation

 Most studies examining the stability and change of patterns in biological coordination have focused on identifying generic bifurcation mechanisms in an already active set of components (see Kelso 1994). A less well understood phenomenon is the process by which previously quiescent degrees of freedom (df ) are spontaneously recruited and active df suppressed. To examine such behavior, in part I we study a single limb system composed of three joints (wrist, elbow, and shoulder) performing the kinematically redundant task of tracing a sequence of two-dimensional arcs of monotonically varying curvature, κ. Arcs were displayed on a computer screen in a decreasing and increasing κ sequence, and subjects rhythmically traced the arcs with the right hand in the sagittal plane at a fixed frequency (1.0 Hz), with motion restricted to flexion-extension of the wrist, elbow, and shoulder. Only a few coordinative patterns among the three joints were stably produced, e.g., in-phase (flexion-extension of one joint coordinated with flexion-extension of another joint) and antiphase (flexion-extension coordinated with extension-flexion). As κ was systematically increased and decreased, switching between relative phase patterns was observed around critical curvature values, κ_c. A serendipitous finding was a strong 2:1 frequency ratio between the shoulder and elbow that occurred across all curvature values for some subjects, regardless of the wrist-elbow relative phase pattern. Transitions from 1:1 to 2:1 frequency entrainment and vice versa were also observed. The results indicate that both amplitude modulation and relative phase change are utilized to stabilize the end-effector trajectory. In part II, a theoretical model is derived from three coupled nonlinear oscillators, in which the relative phases (φ) between the components and the relative joint amplitudes (ρ) are treated as collective variables with arc curvature as a control parameter. Received: 2 February 1996/Accepted in revised form: 13 December 1996     

Spike discharges of skeletomotor neurons during random noise modulated transmembrane current stimulation and muscle stretch

 Spike discharges of skeletomotor neurons innervating triceps surae muscles elicited by white noise modulated transmembrane current stimulation and muscle stretch were studied in decerebrated cats. The white noise modulated current intensity ranged from 4.3 to 63.2 nA peak-to-peak, while muscle stretches ranged from 100 μm to 4.26 mm peak-to-peak. The neuronal responses were studied by averaging the muscle length records centered at the skeletomotor action potentials (peri-spike average, PSA) and by Wiener analysis. Skeletomotor spikes appeared after a sharp peak in PSA of the injected current, preceded by a longer-lasting smaller wavelet of either depolarizing or hyperpolarizing direction. The PSA amplitude was not related to the injected current amplitude nor showed any differences related to the motor unit type. The PSA amplitudes were virtually independent of the stretching amplitude σ, after an initial increase with stretching amplitudes in the range of 15–40 μm (S.D.), or 100–270 μm peak-to-peak.Analyses of cross-spectra indicated a small or absent increase in gain with frequency in response to injected current, but about 20 dB/decade in the range 10–100 Hz in response to muscle stretch. The peaks of both Wiener kernels in response to current injection appear to decrease with the amplitude of injected current, but this decrease was not statistically significant. The narrow first-order kernels suggest that the transfer function between the current input and spike discharge is lowpass with a wide passband, i.e. there is very little change in dynamics. The values of the second-order kernels appear to be nonzero only along the main diagonal. This is characteristic of a simple Hammerstein type cascade, i.e. a zero memory nonlinearity followed by a linear system. Small values of second-order kernels away from the origin and narrow first-order kernels suggest that the linear cascade contributes very little to the overall dynamic response.In contrast to Wiener kernels found in response to current injection, the Wiener kernels in response to stretch showed a decreasing trend with stretch amplitude. The size of the second-order kernels decreased to a somewhat larger extent with input amplitude than that of the first-order kernels, indicating an amplitude-dependent nonlinearity. Overall, the transformation between length and spike output was described as an LNNL cascade with second-order nonlinearities. Received: 1 April 1993/Accepted in revised form: 24 March 1994     

Degradation of anthracene by laccase of Trametes versicolor in the presence of different mediator compounds

 Laccase of Trametes versicolor was generally able to oxidize anthracene in vitro. After 72 h incubation about 35% of the anthracene was transformed stoichiometrically to 9,10-anthraquinone. Transformation of anthracene increased rapidly in the presence of different mediators that readily generate stable radicals: 2,2′-azino-bis-(3-ethylbenzthiazoline-6-sulfonic acid) (ABTS) and 1-hydroxybenzotriazole. For the reaction, the presence of both the laccase and the mediator was necessary. In the presence of 0.005 mM 1-hydroxybenzotriazole this conversion had removed 47% of the anthracene after 72 h; 75% of the substrate was oxidized during this period when ABTS (1 mM) was used as mediator. In contrast to reactions without or with only low concentrations of a mediator, there was a discrepancy between the disappearance of anthracene and the formation of 9,10-anthraquinone in mediator-forced reactions. Coupling-products of mediators with anthracene degradation products were found. Anthracene disappeared nearly completely after incubation for 72 h with laccase in a 0.1 mM solution of 1-hydroxybenzotriazole and was transformed to 9,10-anthraquinone in about 80% yield; 90% of the substrate was transformed in the presence of ABTS (2.0 mM) resulting again in 80% quinone. Phenothiazine was not effective in this system. Received: 9 April 1996/Received last revision: 7 June 1996/Accepted: 10 June 1996     

Frequency dependence of the action-perception cycle for postural control in a moving visual environment: relative phase dynamics

When standing human subjects are exposed to a moving visual environment, the induced postural sway displays varying degrees of coherence with the visual information. In our experiment we varied the frequency of an oscillatory visual display and analysed the temporal relationship between visual motion and sway. We found that subjects maintain sizeable sway amplitudes even as temporal coherence with the display is lost. Postural sway tended to phase lead (for frequencies below 0.2 Hz) or phase lag (above 0.3 Hz). However, we also observed at a fixed frequency, highly variable phase relationships in which a preferred range of phase lags is prevalent, but phase jumps occur that return the system into the preferred range after phase has begun drifting out of the preferred regime. By comparing the results quantitatively with a dynamical model (the sine-circle map), we show that this effect can be understood as a form of relative coordination and arises through an instability of the dynamics of the action-perception cycle. Because such instabilities cannot arise in passively driven systems, we conclude that postural sway in this situation is actively generated as rhythmic movement which is coupled dynamically to the visual motion. Received: 7 September 1993/Accepted in revised form: 2 May 1994     

Tests of a linear model of visual-vestibular interaction using the technique of parameter estimation

The goal of this study was to test whether a superposition model of smooth-pursuit and vestibulo-ocular reflex (VOR) eye movements could account for the stability of gaze that subjects show as they view a stationary target, during head rotations at frequencies that correspond to natural movements. Horizontal smooth-pursuit and the VOR were tested using sinusoidal stimuli with frequencies in the range 1.0–3.5 Hz. During head rotation, subjects viewed a stationary target either directly or through an optical device that required eye movements to be approximately twice the amplitude of head movements in order to maintain foveal vision of the target. The gain of compensatory eye movements during viewing through the optical device was generally greater than during direct viewing or during attempted fixation of the remembered target location in darkness. This suggests that visual factors influence the response, even at high frequencies of head rotation. During viewing through the optical device, the gain of compensatory eye movements declined as a function of the frequency of head rotation (P < 0.001) but, at any particular frequency, there was no correlation with peak head velocity (P > 0.23), peak head acceleration (P > 0.22) or retinal slip speed (P > 0.22). The optimal values of parameters of smooth-pursuit and VOR components of a simple superposition model were estimated in the frequency domain, using the measured responses during head rotation, as each subject viewed the stationary target through the optical device. We then compared the model's prediction of smooth-pursuit gain and phase, at each frequency, with values obtained experimentally. Each subject's pursuit showed lower gain and greater phase lag than the model predicted. Smooth-pursuit performance did not improve significantly if the moving target was a 10 deg × 10 deg Amsler grid, or if sinusoidal oscillation of the target was superimposed on ramp motion. Further, subjects were still able to modulate the gain of compensatory eye movements during pseudo-random head perturbations, making improved predictor performance during visual-vestibular interactions unlikely. We conclude that the increase in gain of eye movements that compensate for head rotations when subjects view, rather than imagine, a stationary target cannot be adequately explained by superposition of VOR and smooth-pursuit signals. Instead, vision may affect VOR performance by determining the context of the behavior. Received: 16 June 1997 / Accepted: 5 December 1997     

Leaf movement of bush bean: a biometeorological perspective

 Leaf movements of bush bean plants were studied at the relatively low photon flux density of 0.2 mmol/m² per s, and air temperatures of 25° and 35° C in a growth chamber. A beta-ray gauge system was used to monitor continuously pulvinus water status and bending. Leaf angles were below the horizontal and were linearly related to the soil water content (R≥−0.91 at 25° C and R≥−0.93 at 35° C). The beta-ray transmission maxima coincided with the stem temperature minima in darkness and vice versa when brightness prevailed as the growth chamber temperature varied with the photoperiod. Leaf angle increased linearly with increased beta-ray transmission. The Q₁₀ temperature coefficient, a measure of the metabolic energy requirement for leaf movement between 25° and 35° C was estimated at 1.8, and the corresponding mean Arrhenius constant at 423 kJ/mol for bush bean. Received: 19 July 1996 / Accepted: 9 September 1996     

Heat stress in Greece

 For 12 selected synoptic stations of the Greek Weather Service, the daily 12 UTC values of the thermal index Predicted Mean Vote (PMV) were calculated for the years 1980 to 1989. The locally varied occurrence of diverse thermal sensation and particularly of strong heat stress were analysed in relation to the human-biometeorological significance. With the help of a statistical model, PMV values of individual stations were transformed into a high-resolution bioclimatic map. The map presents the average annual number of days with at least strong heat stress (PMV>3.0). Received: 13 June 1996 / Revised: 10 February 1997 / Accepted: 18 February 1997     

Investigating the effects of opioid drugs on electrocortical activity using wavelet transform

 Fetal electrocortical activity (ECoG) is characterized by two distinct patterns: HVSA (high voltage, slow activity) and LVFA (low voltage, fast activity). Using the wavelet transform (WT), we recently reported that the frequency characteristics of these two ECoG patterns undergo significant maturational changes prior to birth (Akay et al. 1994a). We now report that fetal ECoG can also be significantly affected by pharmacological agents. In this paper, we compared the effects of two opioid drugs (morphine and [D-Pen², D-Pen⁵]enkephalin, DPDPE) on fetal ECoG, using the chronically instrumented fetal lamb model. Morphine was infused intravenously (i.v.) at 2.5 mg/h, while DPDPE was infused into the lateral cerebroventricle (i.c.v.) at 30 μg/h. The ECoG was analyzed using WT. We performed multiresolution decomposition for four sets of parameters D _{2 j} where −1<j<−4. The four series WTs represent the detail signal bandwidths: (1) 16–32 Hz, (2) 8–16 Hz, (3) 4–8 Hz, (4) 2–4 Hz. The data were subjected to statistical analysis using the Kolmogorov–Smirnov (KS) test. Both morphine and DPDPE resulted in a significant increase in power in the first wavelet band, while power was reduced in the second, third and fourth wavelet bands. In addition, both drugs resulted in a disruption of the normal cyclic pattern between the two ECoG patterns. There was a difference in the time course of action between morphine and DPDPE. This is the first occasion in which continuous ECoG has been subjected to rigorous statistical analysis. The results suggest that the WT-KS method is most suitable for quantitating changes in the ECoG induced by pharmacological agents. Received: 21 January 1994/Accepted in revised form: 15 September 1994     

Surrogate data analysis of sleep electroencephalograms reveals evidence for nonlinearity

 We tested the hypothesis of whether sleep electroencephalographic (EEG) signals of different time windows (164 s, 82 s, 41 s and 20.5 s) are in accordance with linear stochastic models. For this purpose we analyzed the all-night sleep electroencephalogram of a healthy subject and corresponding Gaussian-rescaled phase randomized surrogates with a battery of five nonlinear measures. The following nonlinear measures were implemented: largest Lyapunov exponent L1, correlation dimension D2, and the Green-Savit measures δ2, δ4 and δ6. The hypothesis of linear stochastic data was rejected with high statistical significance. L1 and D2 yielded the most pronounced effects, while the Green-Savit measures were only partially successful in differentiating EEG epochs from the phase randomized surrogates. For L1 and D2 the efficiency of distinguishing EEG signals from linear stochastic data decreased with shortening of the time window. Altogether, our results indicate that EEG signals exhibit nonlinear elements and cannot completely be described by linear stochastic models. Received: 21 December 1995/Accepted in revised form: 19 March 1996     

Banana waste as substrate for α-amylase production by Bacillus subtilis (CBTK 106) under solid-state fermentation

 Bacillus subtilis CBTK 106, isolated from banana wastes, produced high titres of α-amylase when banana fruit stalk was used as substrate in a solid-state fermentation system. The effects of initial moisture content, particle size, cooking time and temperature, pH, incubation temperature, additional nutrients, inoculum size and incubation period on the production of α-amylase were characterised. A maximum yield of 5 345 000 U mg^-1 min^-1 was recorded when pretreated banana fruit stalk (autoclaved at 121 °C for 60 min) was used as substrate with 70% initial moisture content, 400 μm particle size, an initial pH of 7.0, a temperature of 35 °C, and additional nutrients (ammonium sulphate/sodium nitrate at 1.0%, beef extract/peptone at 0.5%, glucose/sucrose/starch/maltose at 0.1% and potassium chloride/sodium chloride at 1.0%) in the medium, with an inoculum-to-substrate ratio of 10% (v/w) for 24 h. The enzyme yield was 2.65-fold higher with banana fruit stalk medium compared to wheat bran. Received: 18 April 1995/Received last revison: 6 May 1996/Accepted: 9 May 1996     

Disposable sensor for biochemical oxygen demand

 Disposable-type microbial sensors were prepared for the determination of biochemical oxygen demand (BOD). The yeast, Trichosporon cutaneum, was directly immobilized on the surface of miniature oxygen electrodes using an ultraviolet crosslinking resin (ENT-3400). The oxygen electrodes (15 mm× 2 mm×0.4 mm) were made on silicon substrates using micromachining techniques. They were Clark-type two-electrode systems with−1021 mV applied to the working electrode. Typical response times of the BOD sensors were in the range of 7–20 min. At 20°C, the sensors’ dynamic range was from 0 to 18 mg/l BOD when a glucose/glutamate BOD standard solution was used. The lower limit of detection was 0.2 mg/l BOD. This value was one order of magnitude lower than that of sensors previously reported. The sensors’ operational lifetime of 3 days was satisfactory for a disposable type. The sensors’ responses were reproducible to within 8% relative standard deviation. The BOD sensors’ were applied to untreated and treated waste waters from industrial effluents and municipal sewage. BOD values determined using these sensors correlated well with those determined by the conventional 5-day BOD determination method. Received: 22 December 1995/Received revision: 19 February 1996/Accepted: 17 March 1996     

Growth of Corynebacterium glutamicum in ammonium- and potassium-limited continuous cultures under high osmotic pressure

 In order to determine the possible effect of nutrient limitations on the response of Corynebacterium glutamicum to a saline osmotic up-shock, the bacteria were grown in continuous cultures, at osmotic pressures of 0.4 osmol/kg and 1.2 osmol/kg, under ammonia and potassium limitation. At the low osmolality of 0.4 osmol/kg, the glutamate and proline levels of 15 mg/g and 5 mg/g dry weight respectively were lower than previously reported in glucose-limited continuous cultures (50 mg/g and 10 mg/g dry weight respectively). On the other hand, the internal trehalose pool was much higher at 40 mg/g dry weight. When the medium osmolality was increased to 1.2 osmol/kg by NaCl addition, under ammonia limitation, the proline content rose from 5 mg/g to 20 mg/g dry weight and the trehalose content from 40 mg/g to 70 mg/g dry weight, whereas the intracellular pool of glutamate remained essentially constant. An increase in the internal sodium content was also observed. Similar results were found for the internal pool of glutamate, proline and trehalose when C. glutamicum was grown under potassium limitations at an osmolality of 1.2 osmol/kg. There were also higher levels of sodium ions, glutamine and alanine. According to the present results, whereas proline was previously reported to be the dominantly accumulated osmoprotectant in C. glutamicum grown under glucose limitations, under ammonia and potassium limitations trehalose represented the dominantly synthesized metabolite. Received: 19 December 1995/Received revision: 9 April 1996/Accepted: 15 April 1996     

黑蝉(C.atrata Fabricius.)鸣声的方向性和第三气门的功能

黑蝉鸣声的波形结构无明显的方向性.单音节的重复周期和调幅脉冲列的间隔(I_1和I_2)分别为9.787±0.813ms、2.286±0.093ms和1.874±0.063ms.幅值特性有明显的方向性.主峰频率(MPF=5.47±0.11kHz)的幅值,头向和背向分别比尾向下降5.9dB和3.9dB,侧向和腹向分别增高1.1dB和2.3dB.两侧第三气门受阻后鸣声的波形结构和音色都产生明显变化.I_1和I_2分别为0.912±0.156ms和1.099±0.113ms,约为正常值的40—59%.有三个谱带,MPF为5775Hz,两侧谱带的峰值频率为4575Hz和7025Hz,分别下降1.5dB和3.4dB.     

Auditory role of lateral trunk channels in cobitid fishes

In cobitid fishes the anterior part of the swimbladder is encapsulated by bone to varying extent. This might diminish the auditory sensitivity of these otophysine fishes by reducing the vibrations of the swimbladder wall in the sound field. However, according to prior studies the auditory thresholds of the cobitid Botia modesta is similar to that of other otophysine fishes. According to anatomical investigation B. modesta has a cranial encapsulation of the anterior part of the swimbladder (camera aerea Weberiana) as expected and in addition special channels stretching laterally from the swimbladder to the outer body wall. These lateral trunk channels are filled with fat and lymph. They form a muscle-free acoustic window beneath the skin, which could be demonstrated by measuring the auditory brainstem response at 400 Hz, 800 Hz, 1500 Hz, and 3000 Hz. Filling the lateral trunk channels with wettex (cotton/rayon staple) resulted in an increase of the auditory thresholds by 13.6–17.6 dB, indicating mechanical damping of the swimbladder. Our experiments demonstrate that the intact lateral trunk channels enhance the hearing sensitivity of cobitid fishes. Accepted: 15 December 1999     

Dynamics of the head-neck complex in response to the trunk horizontal vibration: modeling and identification

Although many studies exist concerning the influence of seat vibration on the head in the seated human body, the dynamic response of the head-neck complex (HNC) to the trunk vibration has not been well investigated. Little quantitative knowledge exists about viscoelastic parameters of the neck. In this study, the dynamics of the HNC is identified when it is exposed to the trunk horizontal (fore-and-aft) vibration. The frequency response functions between the HNC angular velocity and the trunk horizontal acceleration, corresponding to four volunteers, are obtained in the frequency range of 0.5 Hz to 10 Hz. A fourth-order mathematical model, derived by considering a double-inverted-pendulum model for the HNC, is designed to simulate the dynamic response of the HNC to the trunk horizontal vibration. The frequency domain identification method is used to determine the coefficients of the mathematical model of the HNC. Good agreement has been obtained between experimental and simulation results. This indicates that the system, similar to the designed fourth-order model, has mainly two resonance frequencies. The viscoelastic parameters of the neck, including the spring and damping coefficients, are then obtained by use of the optimization method.     

Use of an industrial effluent as a carbon source for denitrification of a high-strength wastewater

Denitrification of a high-strength synthetic wastewater (150 g NO^- ₃ l^-1) was carried out using a wine distillery effluent as an example of an industrial carbon source (22.7 g chemical oxygen demand l^-1). Two configurations were tested: one consisted of an acidogenesis reactor followed by a denitrifying reactor and the other was a single reactor directly fed with the raw effluents. In both cases, denitrification was achieved at a nitrate load of 9.54 g NO^- ₃ l^-1 day^-1 (2.19 g N as NO^- ₃ l^-1 day^-1) with good specific reduction rates: 32.6 mg and 35.2 mg N as NO _x  g volatile suspended solids h^-1, calculated on a single day, for the two-step and the one-step process respectively. Dissimilatory nitrate reduction to ammonium did not occur, even in the one-step process. Received: 26 October 1995/Received revision: 15 February 1996/Accepted: 20 February 1996     

Maximum diving depths of northern rockhopper penguins (Eudyptes chrysocome moseleyi) at Amsterdam Island

 The mean maximum dive depth from 49 foraging bouts by northern rockhopper penguins, measured using capillary-tube depth gauges, was 66±4 m (12–168 m). There were no differences in the maximum dive depths between male and female penguins. Northern rockhopper penguins dived deeper in early than in late creche stages (83±7 vs 57±4 m), and this was associated with probable dietary changes, squid dominating the diet by mass (44%) in November, and fish (64%) in December 1994 at Amsterdam Island. Received: 10 January 1996/Accepted: 31 March 1996     

Selection of an adhesive phenotype of Streptococcus salivarius subsp. thermophilus for use in fixed-bed reactors

 Streptococcus salivarius subsp. thermophilus was cultivated in a chemostat in order to obtain an adhesive phenotype of this strain. When the system was operated at low dilution rates (D<0.2 h^-1) for about 4 weeks, the strain formed a visible film on the surface of the culture vessel. The biofilm cells were not washed out even when dilution rates were increased (D=6.9 h^-1), and this resulted in a high biomass productivity (P=4.1 g l^-1h^-1). On the other hand, when the culture was grown at dilution rates faster than 0.2 h^-1, only the free suspended cells were present in the culture broth, and were washed out at velocities of about 1.0 h^-1. The biomass productivity was consequently lower (P=1.33 g l^-1h^-1) than in the previous case. The selected adhesive phenotype was grown on different glass beads and the possibility of lactate fermentation in a continuous and semicontinuous mode was demonstrated. Received: 16 August 1995/Received revision: 18 March 1996/Accepted: 25 March 1996