Physiological and psychological assessment of sound

 The psycho-physiological effects of several sound stimulations were investigated to evaluate the relationship between a psychological parameter, such as subjective perception, and a physiological parameter, such as the heart rate variability (HRV). Eight female students aged 21–22 years old were tested. Electrocardiogram (ECG) and the movement of the chest-wall for estimating respiratory rate were recorded during three different sound stimulations; (1) music provided by a synthesizer (condition A); (2) birds twitters (condition B); and (3) mechanical sounds (condition C). The percentage power of the low-frequency (LF; 0.05≤0.15 Hz) and high-frequency (HF; 0.15≤0.40 Hz) components in the HRV (LF%, HF%) were assessed by a frequency analysis of time-series data for 5 min obtained from R-R intervals in the ECG. Quantitative assessment of subjective perception was also described by a visual analog scale (VAS). The HF% and VAS value for comfort in C were significantly lower than in either A and/or B. The respiratory rate and VAS value for awakening in C were significantly higher than in A and/or B. There was a significant correlation between the HF% and the value of the VAS, and between the respiratory rate and the value of the VAS. These results indicate that mechanical sounds similar to C inhibit the para-sympathetic nervous system and promote a feeling that is unpleasant but alert, also suggesting that the HRV reflects subjective perception. Received: 28 June 1996 / Revised: 6 January 1997 / Accepted: 10 January 1997     

Coarse-graining spectral analysis: new method for studying heart rate variability.

Heart rate variability (HRV) spectra are typically analyzed for the components related to low- (less than 0.15 Hz) and high- (greater than 0.15 Hz) frequency variations. However, there are very-low-frequency components with periods up to hours in HRV signals, which might smear short-term spectra. We developed a method of spectral analysis suitable for selectively extracting very-low-frequency components, leaving intact the low- and high-frequency components of interest in HRV spectral analysis. Computer simulations showed that those low-frequency components were well characterized by fractional Brownian motions (FBMs). If the scale invariant, or self-similar, property of FBMs is considered a new time series (x') was constructed by sampling only every other point (course graining) of the original time series (x). Evaluation of the cross-power spectra between these two (Sxx') showed that the power of the FBM components was preserved, whereas that of the harmonic components vanished. Subtraction of magnitude of Sxx from the autopower spectra of the original sequence emphasized only the harmonic components. Application of this method to HRV spectral analyses indicated that it might enable one to observe more clearly the low- and high-frequency components characteristic of autonomic control of heart rate.     

Fine structure of the low-frequency spectra of heart rate and blood pressure

Background  

The aim of this study was to explore the principal frequency components of the heart rate and blood pressure variability in the low frequency (LF) and very low frequency (VLF) band. The spectral composition of the R–R interval (RRI) and systolic arterial blood pressure (SAP) in the frequency range below 0.15 Hz were carefully analyzed using three different spectral methods: Fast Fourier transform (FFT), Wigner-Ville distribution (WVD), and autoregression (AR). All spectral methods were used to create time–frequency plots to uncover the principal spectral components that are least dependent on time. The accurate frequencies of these components were calculated from the pole decomposition of the AR spectral density after determining the optimal model order – the most crucial factor when using this method – with the help of FFT and WVD methods.     

Natural environmental associations in a 50-day human electrocardiogram

Meteorotropic associations of heart rate (HR) and HR variability (HRV) were investigated in a clinically healthy 48-year-old man in Kiev. His electrocardiogram (ECG) was determined over 50 days by fitting him with an ambulatorily wearable device; various natural physical environmental variables were also monitored. The mean inter-beat interval, the standard deviation of these intervals, the spectral power in several frequency ranges, the power ratio of the approx. 10.5-s/approx. 3.6-s spectral components and other aspects of HRV were computed over consecutive 14.4-min intervals. Together with ordinary meteorological variables, geomagnetic disturbances (GMD) and fluctuations of atmospheric pressure (FAP) in the range 0.01–0.10 Hz (10–100 s) were measured. The assessable infradian spectra (with frequencies lower than 1 cycle/28 hours) of all HRV parameters showed two major components with periods of about 3.5 and 10 days. Two environmental variables, FAP and wind speed (with which FAP is closely related), revealed both of these rhythms and showed the greatest cross-spectral coherence (0.70–0.98) with corresponding oscillations of HRV. Less specific but statistically significant product-moment correlations with major HRV indices were also found; most of these were with FAP, but correlations with air temperature, humidity, wind speed and geomagnetic disturbances were also found. Long-term ECG recording, essential in the detection of infradian rhythms, proved to be sensitive to physical environmental variables, notably meteorological ones. FAP, usually neglected since its role has not been considered in previous biometeorological studies, or some factor closely related to FAP is probably involved in synchronizing or influencing the approximately 3.5-day HR and HRV rhythms in humans. Received: 27 July 2000 / Revised: 13 November 2000 / Accepted: 28 November 2000     

Associations between Diurnal 24-Hour Rhythm in Ambulatory Heart Rate Variability and the Timing and Amount of Meals during the Day Shift in Rotating Shift Workers

It has not hitherto been clarified whether there is an association between dietary behavior and circadian variation in autonomic nervous system activity among shift workers. This study examines diurnal 24-h rhythm in heart rate variability (HRV) and dietary behavior among rotating shift workers, while taking into account the sleep-wake cycle and physical activity. The subjects were 11 female and 2 male nurses or caregivers working in a rotating 2-shift system at a health care facility. All the subjects were asked to undergo 24-h electrocardiograph and step count recordings, and to record the time of each meal and the amounts of each food and beverage consumed. Coarse graining spectral analysis was used for approximately 10-min segments of HRV to derive the total power (TOT: >0.04 Hz) of the periodic components and the integrated power of periodic components in the low-frequency (LF: 0.04–0.15 Hz) and high-frequency (HF: >0.15 Hz) ranges. Then the ratio of HF power to TOT (HF nu) and the ratio of LF power to HF power (LF/HF) were calculated to assess cardiac vagal tone and cardiac sympathovagal balance, respectively. Single cosinor analysis was used to obtain 24-h period variations in both variables of HRV. Acrophases of HF nu and LF/HF expressed in time since awakening were significantly (p<0.05) delayed for subjects having breakfast at a later time after awakening. Multivariable regression analysis indicated that the timing of breakfast, the ratio of energy intake at dinner to total energy intake, and total energy intake were correlated to the acrophases of HF nu and/or LF/HF. These results suggest that the phase angle between circadian variation in cardiac autonomic nervous system activity and the sleep-wake cycle may be associated with dietary behavior in shift workers.     

Intraindividual validation of heart rate variability indexes to measure vagal effects on hearts

Heart rate variability (HRV) has been widely used as a measure of vagal activation in physiological, psychological, and clinical examinations. We studied the within-subject quantitative relationship between HRV and vagal effects on the heart in different body postures during a gradually decreasing vagal blockade. Electrocardiogram and respiratory frequency were measured in subjects (8 endurance athletes and 10 participants of nonendurance sports) in supine, sitting, and standing postures before the blockade, under vagal blockade (atropine sulfate, 0.04 mg/kg), and four times during a 150-min recovery from the blockade. Fast Fourier transform was used to calculate low-frequency power (LFP, 0.04-0.15 Hz), high-frequency power (HFP, 0.15-0.40 Hz), and total power (TP, 0.04-0.40 Hz). A within-subject linear regression analysis of recovery time on each HRV index was conducted. Complete vagal blockade decreased all HRV significantly, particularly HFP (P < 0.001). A linear fit explained a large portion of the within-subject variance between recovery time and natural log-transformed (ln) HRV indexes in every posture, with coefficients of determination (R2) in the supine posture [means (SD)]: 98 (SD 2)% for mean R-R interval, 87 (SD 10)% for lnLFP, 87 (SD 13)% for lnHFP, and 91 (SD 10)% for lnTP. Neither body posture nor endurance-training background had an impact on R2 values. There was marked between-subject variation in the R2 values, slopes, and intercepts. In conclusion, all HRV, particularly HFP, is predominantly under vagal control. Within subjects, lnLFP, lnHFP, and lnTP increased linearly with the gradually decreasing vagal blockade in all postures.     

Inhibition of heart rate variability during sleep in humans by 6700 K pre-sleep light exposure

Two different spectral analyses of heart rate (HR) variability (HRV) were performed on seven young male subjects to evaluate the effects of different color temperatures of light exposure (6700 K, 5000 K, 3000 K) before sleep on cardiac vagal activity. In investigating HRV, we used an ordinary fast Fourier transform (FFT) and coarse graining spectral analysis (CGSA), which selectively extracts random fractal components from a given time series. The results showed that suppressions of HR during sleep after 6700 K light exposure were more inhibited than the other two lighting conditions. Increases in high-frequency (HF) components of HRV during sleep were also inhibited by 6700 K pre-sleep lighting. These results indicate that pre-sleep exposure to light of a higher color temperature may inhibit the enhancement of cardiac vagal activity during sleep. Moreover, significant HF alterations were shown in fractal-free HF (not in ordinary HF) components by CGSA. Because the HF component originates from respiratory sinus arrhythmia with periodical fluctuations, CGSA may be an appropriate approach for HRV evaluation during sleep.     

Frequency characteristics of long-term heart rate variability during constant-routine protocol

The effects of such behavioral factors as physical activity, food intake, and circadian rhythm on long-term heart rate variability (HRV) in humans remain poorly understood. We therefore studied their effects on HRV using a constant-routine protocol that included simultaneous core body temperature (CBT) correction. Seven healthy subjects completed the constant-routine and daily-routine protocols, during which HRV and CBT were continuously monitored. During the constant routine, subjects were kept awake for 27 h in a semirecumbent posture with minimal physical activity; small isocaloric meals were provided every 2 h. During the daily routine, subjects carried on their lives normally. Data were analyzed using generic spectral analysis based on a fast Fourier transform; coarse-graining spectral analysis was also used to eliminate periodicity due to the regular meals for raw HRV and for the CBT-corrected HRV without circadian and/or low-frequency ultradian components. The results showed that 1) the power spectra of HRV in the constant routine and daily routine had similar power-law scalings at frequencies above approximately 10(-3.5) Hz, while 2) below that crossover frequency, HRV was smaller in the constant routine than in the daily routine, with the difference becoming significant (P < 0.05) at <10(-4) Hz, 3) coarse-graining spectral analysis eliminated diet-induced peaks in generic spectral analysis-based HRV spectra during the constant routine and emphasized the crossover at approximately 10(-3.5) Hz, and 4) CBT correction did not alter the results. Below a frequency of approximately 10(-3.5) Hz (a period >1 h), HRV is strongly influenced by behavioral factors; above that crossover frequency, HRV is behavior independent, possibly reflecting an intrinsic regulatory system.     

A measure of heart rate variability is sensitive to orthostatic challenge in women with chronic fatigue syndrome

The use of symptoms generated by head up tilt (HUT) is not a useful tool in identifying chronic fatigue syndrome (CFS). We investigated whether heart rate variability (HRV) assessed early during HUT might be useful. A sample of 46 female subjects (24 with CFS and 22 sedentary, age-matched healthy controls; CON) who had exhibited no difference in time to syncope during tilt was examined for HRV responses to 10 min of 70 degrees HUT after 5 min of baseline in the supine position. HRV data were analyzed by the method of coarse graining spectral analysis. Variables compared between groups included mean and standard deviation (SD(RRI)) of RR intervals (RRI), amplitudes of low- (A(LF); 0.04-0.15 Hz) and high-frequency (A(HF); >0.15 Hz) harmonic as well as aperiodic, fractal (A(FR); 1/f(beta)) spectral components, the spectral exponent beta, and the difference in these values between baseline and HUT for each subject. In the supine baseline, only mean RRI was significantly (P < 0.01) lower in CFS than in CON. During HUT, however, mean RRI (P < 0.01), SD(RRI) (P < 0.01), A(HF) (P < 0.05), and A(FR) (P < 0.01) were significantly lower in CFS than in CON. When the difference in values between baseline and HUT for each subject was examined, only the difference for A(FR) (deltaA(FR)) was significantly (P < 0.01) lower in CFS than in CON, suggesting that A(FR)is a disease-specific response of HRV to HUT. When a cut-off level was set to deltaA(FR) = -2.7 msec, the sensitivity and the specificity in differentiating CFS from controls were 90% and 72%, respectively. The data suggest that a decrease in aperiodic fractal component of HRV in response to HUT can be used to differentiate patients with CFS from CON.     

Effects of immersion in water containing high concentrations of CO2 (CO2-water) at thermoneutral on thermoregulation and heart rate variability in humans

Immersion in high concentrations of CO₂ dissolved in freshwater (CO₂-water) might induce peripheral vasodilatation in humans. In this study, we investigated whether such immersion could affect the autonomic nervous system in humans using spectral analysis of heart rate variability. Ten healthy men participated in this study. Tympanic temperature, cutaneous blood flow and electrocardiogram (ECG) were measured continuously during 20 min of immersion in CO₂-water. The ECG was analyzed by spectral analysis of R–R intervals using the maximal entropy method. The decrease in tympanic temperature was significantly greater in CO₂-water immersion than in freshwater immersion. Cutaneous blood flow at the immersed site was significantly increased with CO₂-water immersion compared to freshwater. The high frequency component (HF: 0.15–0.40 Hz) was significantly higher in CO₂-water immersion than in freshwater immersion, but the low frequency (LF: 0.04–0.15 Hz) /high frequency ratio (LF/HF ratio) was significantly lower in CO₂-water immersion than in freshwater immersion. The present study contributes evidence supporting the hypothesis that CO₂-water immersion activates parasympathetic nerve activity in humans.     

The effect of sinus node depression on heart rate variability in humans using zatebradine, a selective bradycardic agent

Zatebradine is a bradycardic agent with a selective effect on the pacemaker current in the sinus node. The effect of such drugs on heart rate variability is not known. Thirty-six patients without structural heart disease were randomly assigned to receive 10 mg of zatebradine i.v. (n = 24) or isotonic saline (n = 12). Heart rate variability (HRV) was recorded as power in the very low frequency (VLF, 0.003-0.040 Hz), low frequency (LF, 0.040-0.150 Hz), and high frequency (HF, 0.150-0.400 Hz) spectral bands as well as total power (TP, 0.003-0.400 Hz) during 5-min ECG acquisitions at baseline, 30, and 60 min following the start of the infusion. No change in heart rate variability was detected in the control group. Zatebradine significantly reduced heart rate variability at 60 min in all frequency bands: VLF (-12+/-4%, p<0.001), LF (-19+/-4%, p<0.001), and HF (-26+/-5%, p<0.001). The reduction in HRV following zatebradine is due to depression of sinus node response to all external stimuli and underscores the need for documentation of normal sinus node function in HRV research.     

High resolution electrocardiography in healthy dogs: time domain parameters and comparison of the non-stationary (Wigner distribution) versus standard stationary frequency domain analysis methods

Fractionated heart activation can be detected as late potentials from surface recordings of signal-averaged electrocardiograms (SA ECG) which are considered as a marker of sustained ventricular tachycardia. For animal studies, reference values in time and frequency domain analyses are essentially missing. In the present study, we have established reference values in SA ECG time domain analysis and time-frequency representation of heart activation in healthy dogs. A group of 25 healthy mongrel dogs (body weight 12-15 kg) was investigated. Wigner distribution and our modification of Fast Fourier transform (FFT), gliding window FFT, was applied in SA ECG frequency domain analysis. Reference values in time domain SA ECG were established. Time and voltage criteria were adapted to short duration of heart cycle and fast voltage decrement of the QRS complex in dogs. Wigner distribution and gliding window FFT were applied in order to describe mean heart activation in the frequency domain. Contribution of higher frequencies (30-80 Hz) was detected by both frequency analysis methods in the second third of ventricular activation in healthy animals. Presented results could offer a basis for further experimental arrhythmologic studies.     

Fractal and periodic heart rate dynamics in fetal sheep: comparison of conventional and new measures based on fractal analysis

The physiological significance of spectral and fractal components of spontaneous heart rate (HR) variability in the fetus remains unclear. To examine the relationship between circadian rhythms in different measures of HR variability, R-R interval time series obtained by fetal ECGs were recorded continuously over 24 h in five pregnant sheep at 116-125 days gestation. Conventional measures of short-term (STV) and long-term variability (LTV), low-frequency (LF; 0.025-0.15 cycles/beat) and high-frequency (HF; 0.2-0.5 cycles/beat) spectral powers, the LF-to-HF ratio, and fractal dimension values were calculated from 24-h ECG recordings and quantified every 60 min. STV, LTV, and LF and HF spectral powers were minimal during the day but increased significantly to their highest values at night. We found a significant positive correlation between these measures, whereas the cosinor method showed significant similarity between their circadian rhythm patterns. Fetal R-R intervals also exhibited fractal structures. Fetal HR variability had a fractal structure, which was similar between day and night. These results suggested that the circadian rhythms exhibited by STV and LTV during the day were mainly due to changes in frequency components rather than to fractal components of fetal HR fluctuation.     

Temporal contribution of body movement to very long-term heart rate variability in humans

A newly developed, very long-term ( approximately 7 days) ambulatory monitoring system for assessing beat-to-beat heart rate variability (HRV) and body movements (BM) was used to study the mechanism(s) responsible for the long-period oscillation in human HRV. Data continuously collected from five healthy subjects were analyzed by 1) standard auto- and cross-spectral techniques, 2) a cross-Wigner distribution (WD; a time-frequency analysis) between BM and HRV for 10-s averaged data, and 3) coarse-graining spectral analysis for 600 successive cardiac cycles. The results showed 1) a clear circadian rhythm in HRV and BM, 2) a 1/f (beta)-type spectrum in HRV and BM at ultradian frequencies, and 3) coherent relationships between BM and HRV only at specific ultradian as well as circadian frequencies, indicated by significant (P < 0.05) levels of the squared coherence and temporal localizations of the covariance between BM and HRV in the cross-WD. In a single subject, an instance in which the behavioral (mean BM) and autonomic [HRV power >0.15 Hz and mean heart rate (HR)] rhythmicities were dissociated occurred when the individual had an irregular daily life. It was concluded that the long-term HRV in normal humans contained persistent oscillations synchronized with those of BM at ultradian frequencies but could not be explained exclusively by activity levels of the subjects.     

Preeclamptic pregnancy is associated with increased sympathetic and decreased parasympathetic control of HR

Previous work from our laboratory using heart rate variability (HRV) has demonstrated that women before menopause have a more dominant parasympathetic and less effective sympathetic regulations of heart rate compared with men. Because it is still not clear whether normal or preeclamptic pregnancy coincides with alternations in the autonomic functions, we evaluated the changes of HRV in 17 nonpregnant, 17 normotensive pregnant, and 11 preeclamptic women who were clinically diagnosed without history of diabetic neuropathy, cardiac arrhythmia, and other cardiovascular diseases. Frequency-domain analysis of short-term, stationary R-R intervals was performed to evaluate the total variance, low-frequency power (LF; 0.04-0.15 Hz), high-frequency power (HF; 0.15-0.40 Hz), ratio of LF to HF (LF/HF), and LF in normalized units (LF%). Natural logarithm transformation was applied to variance, LF, HF, and LF/HF for the adjustment of the skewness of distribution. We found that the normal pregnant group had a lower R-R value and HF but had a higher LF/HF and LF% compared with the nonpregnant group. The preeclamptic group had lower HF but higher LF/HF compared with either the normal pregnant or nonpregnant group. Our results suggest that normal pregnancy is associated with a facilitation of sympathetic regulation and an attenuation of parasympathetic influence of heart rate, and such alterations are enhanced in preeclamptic pregnancy.     

Spectral analysis of photoplethysmographic signals: The importance of preprocessing

Heart rate variability (HRV) is an important and useful index to assess the responses of the autonomic nervous system (ANS). HRV analysis is performed using electrocardiography (ECG) or photoplethysmography (PPG) signals which are typically subject to noise and trends. Therefore, the elimination of these undesired conditions is very important to achieve reliable ANS activation results. The purpose of this study was to analyze and compare the effects of preprocessing on the spectral analysis of HRV signals obtained from PPG waveform. Preprocessing consists of two stages: filtering and detrending. The performance of linear Butterworth filter is compared with nonlinear weighted Myriad filter. After filtering, two different approaches, one based on least squares fitting and another on smoothness priors, were used to remove trends from the HRV signal. The results of two filtering and detrending methods were compared for spectral analysis accomplished using periodogram, Welch's periodogram and Burg's method. The performance of these methods is presented graphically and the importance of preprocessing clarified by comparing the results. Although both filters have almost the same performance in the results, the smoothness prior detrending approach was found more successful in removing trends that usually appear in the low frequency bands of PPG signals. In conclusion, the results showed that trends in PPG signals are altered during spectral analysis and must be removed prior to HRV analysis.     

The effect of dim light at night on cerebral hemodynamic oscillations during sleep: A near-infrared spectroscopy study

Recent studies have reported that dim light at night (dLAN) is associated with risks of cardiovascular complications, such as hypertension and carotid atherosclerosis; however, little is known about the underlying mechanism. Here, we evaluated the effect of dLAN on the cerebrovascular system by analyzing cerebral hemodynamic oscillations using near-infrared spectroscopy (NIRS). Fourteen healthy male subjects underwent polysomnography coupled with cerebral NIRS. The data collected during sleep with dim light (10 lux) were compared with those collected during sleep under the control dark conditions for the sleep structure, cerebral hemodynamic oscillations, heart rate variability (HRV), and their electroencephalographic (EEG) power spectrum. Power spectral analysis was applied to oxy-hemoglobin concentrations calculated from the NIRS signal. Spectral densities over endothelial very-low-frequency oscillations (VLFOs) (0.003–0.02 Hz), neurogenic VLFOs (0.02–0.04 Hz), myogenic low-frequency oscillations (LFOs) (0.04–0.15 Hz), and total LFOs (0.003–0.15 Hz) were obtained for each sleep stage. The polysomnographic data revealed an increase in the N2 stage under the dLAN conditions. The spectral analysis of cerebral hemodynamics showed that the total LFOs increased significantly during slow-wave sleep (SWS) and decreased during rapid eye movement (REM) sleep. Specifically, endothelial (median of normalized value, 0.46 vs. 0.72, p = 0.019) and neurogenic (median, 0.58 vs. 0.84, p = 0.019) VLFOs were enhanced during SWS, whereas endothelial VLFOs (median, 1.93 vs. 1.47, p = 0.030) were attenuated during REM sleep. HRV analysis exhibited altered spectral densities during SWS induced by dLAN, including an increase in very-low-frequency and decreases in low-frequency and high-frequency ranges. In the EEG power spectral analysis, no significant difference was detected between the control and dLAN conditions. In conclusion, dLAN can disturb cerebral hemodynamics via the endothelial and autonomic systems without cortical involvement, predominantly during SWS, which might represent an underlying mechanism of the increased cerebrovascular risk associated with light exposure during sleep.     

Detection of low- and high-frequency rhythms in the variability of skin sympathetic nerve activity

Spectral analysis of skin blood flow has demonstrated low-frequency (LF, 0.03-0.15 Hz) and high-frequency (HF, 0.15-0.40 Hz) oscillations, similar to oscillations in R-R interval, systolic pressure, and muscle sympathetic nerve activity (MSNA). It is not known whether the oscillatory profile of skin blood flow is secondary to oscillations in arterial pressure or to oscillations in skin sympathetic nerve activity (SSNA). MSNA and SSNA differ markedly with regard to control mechanisms and morphology. MSNA contains vasoconstrictor fibers directed to muscle vasculature, closely regulated by baroreceptors. SSNA contains both vasomotor and sudomotor fibers, differentially responding to arousals and thermal stimuli. Nevertheless, MSNA and SSNA share certain common characteristics. We tested the hypothesis that LF and HF oscillatory components are evident in SSNA, similar to the oscillatory components present in MSNA. We studied 18 healthy normal subjects and obtained sequential measurements of MSNA and SSNA from the peroneal nerve during supine rest. Measurements were also obtained of the electrocardiogram, beat-by-beat blood pressure (Finapres), and respiration. Spectral analysis showed LF and HF oscillations in MSNA, coherent with similar oscillations in both R-R interval and systolic pressure. The HF oscillation of MSNA was coherent with respiration. Similarly, LF and HF spectral components were evident in SSNA variability, coherent with corresponding variability components of R-R interval and systolic pressure. HF oscillations of SSNA were coherent with respiration. Thus our data suggest that these oscillations may be fundamental characteristics shared by MSNA and SSNA, possibly reflecting common central mechanisms regulating sympathetic outflows subserving different regions and functions.