TH-F3-4 SCIS&ISIS2006 @ Tokyo, Japan (September 20-24, 2006)

EEG Changes After Bhramari Pranayama Rajkishore Prasad1, Fumitoshi Matsuno1, Hovagim Bakardjian 2, Francois Vialatte2 and Andrzej Cichocki2 1. Man Machine Interface Lab., University of Electro Communication, Japan 2. Science Institute, RIKEN, Japan E-mails :1.{ kishore,matsuno}@hi.mce.uec.ac.jp, 2. {hova, fvialatte, cichocki}@brain.riken.jp

Abstract—This paper presents changes in the EEG pattern observed after Bhramari Pranayama (BP) or like breathing exercise. BP has been found effective in healing many mental problems such as tension, stress, hypertension, etc and it brings relaxation in the practitioners of it. It has also been said that BP is very quick in producing results. The analysis of EEG data recorded before and after BP shows similar and positive changes in different brainwave patterns, related to relaxed state, both for an experienced and a new subjects. Index Terms—yoga, EEG, relaxation vibro-acoustic, Bhramari pranayama. Fig.1 One of the easiest body posture for doing Bhramari Pranayama. Subject plugs ear by finger and makes humming sound with exhalation through nasal airways. I. INTRODUCTION hramari Pranayama (BP) is one of the important types of eye brows, middle fingers are used to press root of the nose B yoga related with breath control. Today yoga is well near the corners of the eyes and other fingers are placed along introduced and is well accepted health caring practice in nose such that little fingers lie below. Different variants of BP everyday's life. Regular practice of Yoga is preventive, curative differ mainly in body posture. In on of its difficult form and recuperative. It helps healing process to accelerate process humming sound is made both during exhaling and inhaling. of recovery [1],[2]. The most ancient use of the word yoga can While doing BP concentration is kept at the third eye point be found in the oldest ritual books, known as VEDA, of Indian located between eyebrows as shown in Fig.1. The volume of civilization, however, it has been assumed that Yoga was humming should be loud enough for vibrations in the brain. originated in India nearly 6000 years ago as a part of Indian The humming sound should be steady and not wavering like a traditional treatment system known as Aryuveda. All the motorcycle starting up. It should be effortless and mellow. ancient records available on the effects and benefits of yoga Healing effects of BP are numerous, however, it has been have been results of subjective observations of different sages, mainly recommended for reducing mental problems such as rishis (ancient Indian researchers) and yoga masters, for tension, cerebral tension, hypertension, anxiety and controlling thousands of years. However, in the modern age, studies have blood pressure [6],[7],[8]. It has also been said that it brings been started on different yogic acts using modern scientific bliss and improves psychic sensitivity and for subtle methods and many validatory results have been reported too sound vibration. However, very few scientific studies on [3],[4],[5]. effects of this pranayama have been done. In [7], authors made Pranayama is also a type of yoga popularly known as Hath studies on effects and benefits of BP on the pregnant women yoga. Practice of Hath yoga in some sense is common with and have concluded with great usability of it for the preparation Raja Yoga or Power Yoga. According to Hath Yoga Pradipika, for labor. In one of our very recent studies [9], first two authors an ancient Indian text on Hath yoga, Hath yoga leads to Raja have shown that humming sound produced during BP and that Yoga [6]. The Sanskrit word Pranayama contains two of bumble bee bears similarity in spectral contents and have meaningful segments namely (this means vital force) proposed that humming sound plays key role in the healing and (this means control). So pranayama literally means a process of the BP. In the present study we provide some of our yogic act performed for controlling flow of the vital energy that observations on changes in EEG brainwaves after performing governs all the physiological process of our body. There are BP. To best of our knowledge no EEG study has been made on many types of pranayama among which BP is an important BP in past, however, many EEG studies have been made on item [6]. The term Bhramari is also a Sanskrit word which in different types of pranayama. For example, in [10],[11] EEG this context means like a bumble bee. Thus in BP subject has to studies on alternate nostrils breathing pranayama has been produce sound like bumble bee while exhaling strictly through made and it has been shown that it improves balances between nasal airways, keeping oral cavity closed at the lips, ears closed right and left hemispheres of the brain. In other studies on the by fingers and eyelids are also shut down. There are many yogic breathing it has been reported that pranayamas are variants of BP. One of the easiest body posture during BP is effective in improving cognitive performance, bringing shown in Fig.1.In another variants ear is closed by thumb by relaxation, improving spatial memory performances [12]. In pressing ear flaps and putting index finger on forehead along [13], EEG study on Sudarshan Kriya, a type of pranayama, has been done. Based on such studies we also decided to make EEG

- 390 -

studies on the BP as it has been found useful in reducing mental Separation (BSS) algorithm which estimates independent problems and most importantly it has been said that it is very components of the signal and arrange them in the decreasing quick in bringing effects. order of their variance and linear predictability. The EEG Rest of the paper is organized as follows. In the Section II study components corresponding to artifacts will have lower values methods have been dealt. In Section III EEG data analysis of variance and linear predictability which can be eliminated methods and results have been presented. Section IV presents and rest components can be used to reconstruct the sensor conclusions and future research work which are followed by signals known as deflated signals [15]. Similar approaches references. have been used here too. Fig.3 shows raw EEG data at 21 sensors and Fig.4 shows EEG data obtained after AMUSE II. EEG RECORDING AND METHODS filtering where it can be seen how components have been In this preliminary study two subjects were chosen one with 4 arranged in order of increasing complexities. months regular experience of BP and other with no prior It is an established fact that level of is related to experience of BP, however, second subject was trained for how different brainwaves signals known as delta (δ ) wave with to do BP before the EEG recording. Subjects were chosen from different cultural and ethnic groups. Subjects with such contrast were chosen to see if the BP produces quick effect or not and if it can works in similar ways on different people. In this study our aims have been to objectify effect of BP, so experiment was designed to record EEG before BP and after BP. The BP was done for 15-20 rounds by each subject. EEG was recorded with 128 channel BIOSEMI system at RIKEN, Japan. The sensor layout used in the EEG recording is shown in Fig. 2. For the purpose of analysis, we selected the sensor positions located in the part of , and from both the left and right sides of the head. The reason for this is based on the fact that in BP concentration is kept at the third eye point which has biological connection with the pineal gland. Also, the humming sound is produced in BP and its interaction will be related with the temporal lobe of the brain. In this way total 21 sensor signals were chosen covering all these regions. The selected EEG signals were filtered using AMUSE algorithm [14] . AMUSE algorithm is a Blind Signal

Fig.3 Raw EEG signals at for 21 sensors (x-axis time in sec, scale is arbitrary).

Fig.2 Sensor layout of 128 Channel biosemi system. Radial Fig.4 AMUSE filtered EEG signals arranged in order of positions of the sensor locations were a little modified to forbid increasing complexities (x-axis time in sec, scale is partial overlaps in figure. arbitrary).

- 391 -

frequency range 0.2 Hz to 4 Hz, theta (θ ) wave with frequency state. This is related with waking and responsible for range 4Hz to 8 Hz, alpha (α ) wave with frequency range 8Hz flashing vivid imageries before the minds eye. It has been to 14 Hz, beta ( β ) wave with frequency range 14 Hz to 30 Hz identified as gateway of learning. Its dominance increases and Gamma ( γ ) wave with frequency range 30 Hz to 80 Hz. creativity, learning, reduces stress, awakens institution and other extraordinary skills. There have been The EEG signal obtained from AMUSE algorithm were passed consistent reporting in many previous research papers that through zero phase shift band pass filters with pass bands during relax state theta activity increases [16],[17],[18]. As it corresponding to that of δ ,θ,α, β,and γ brainwaves to filter can be imbued from Fig.5, Fig 6, Fig.7 and Fig.8 that there are out these signals. The temporal variation of Relative Spectral considerable changes in RSP of theta activity of the both Power (RSP) defined as the ratio of Band Spectral Power (BSP) subjects. However, amount of change in theta power in the first and Total Spectral Power (TSP) as follows subject is more. From this it can be explained why BP has been described as the yogic practice for improving intuition. The

BSPi increase in theta power has also been associated with decreased RSPi = , i ∈{δ ,θ,α, β,γ}, (1) activity in the brain arousal system and onset process TSP [18],[19],[20] . Thus BP might also be improving sleep due increase in theta power. Alpha waves are present during closed were calculated for each brainwave patterns. The TSP was eye state, light . It gives deep relaxation. But past calculated from spectral components falling within frequency reports in this regard are contradictory. Some reports say that range of 0.5 Hz to 80 Hz. For the spectral analysis by the Fast alpha power increase and some research reports say that alpha Fourier Transform method, Hanning window of 1s with 50% power decreases in relaxation [21], [22]. In the case of BP it can overlap was used. The EEG data were sampled at 256 Hz. be seen in figures displaying RSP for both subjects that the In order to measure interplay and functional interactions alpha power is increased. The beta waves are dominating in between different parts of the brain, coherence function normal waking state, alertness, concentration and visual between different sensor signals are measured. The coherence acuteness.The temporal variation of RSP inour case for beta function is estimated from cross power spectrum and is known as Magnitude Squared Coherence (MSCH). For two signals x and y, the MSCH at frequency f can be given by

P Before BP

S 1 R

2 After BP

d e

S f z 1 ( ) i 2 xy (2) l a MSCH ( f ) = m r

o 1 Sxx ( f ) * S yy ( f ) N 0 10 20 30 40 50 60 Time(s) Before BP

P After BP Sxy(f)| Sxx(f)| Syy(f)| S 1 where | is cross-power spectrum, | and | are R

d

e 0.8

auto-power spectra. The value of MSCH lies between 0 and 1. z i l a 0.6 A higher value of coherence within this range for two different m r

o EEG signals indicates that the corresponding two brain areas N 0 10 20 30 40 50 60 Time(s) are working together in a holistic and synchronous fashion. Before BP P After BP

Thus the coherence values quantify the linear dependencies S 1 R

between two different brain locations. d e

z 0.5 i l

a m r

o 0 N 0 10 20 30 40 50 60 III. RESULTS AND DISCUSSION Before BP Time(s)

P After BP

The estimated temporal variation of RSP, as per Eq.(1), from S 1 R

d

the deflated EEG signals, corresponding to left side of the head, e

z 0.5 i l for the first subject is shown in Fig.5. The RSP for the right side a m r of the head is shown in the Fig.6. The temporal variations of o 0 N 0 10 20 30 40 50 60 RSP for the second subject are shown in Fig.7 and Fig.8. These Time(s) Before BP Gamma wave figures revel the effect of BP on the pattern of different P After BP

S 1 R brainwaves which are representatives of subject's mental state. d e

z 0.5 i The changes in delta and beta cortical activities are highly l a m r effected by artifacts so changes in alpha, theta and gamma o 0 N 0 10 20 30 40 50 60 brainwaves are important in this report. These plots show that Time(s) the BP produces very obvious and noticeable changes in all the brainwave patterns. Changes in brainwave patterns and their interpretations have been studied by different researchers since long. Brain activity of each frequency range has characteristic Fig.5 Showing temporal variation of RSP in delta, theta, blueprints and produces distinctive state of consciousness. alpha, beta and gamma brainwave pattern before and after Delta brain activity is related with sleep and deep sleep. Theta BP for the first subject.(Data used from the left side of the brainwaves also occur in sleep and dominate in the meditation head )

- 392 -

diagonal of the coherence matrix, to keep clarity of the figure. The shown coherence values are averaged over frequency bins. Delta wave P Before BP

S 1 In both subjects there is change in coherency of EEG of R

After BP d

e different scalp areas. Both subjects show increase in coherency

z 1 i l

a among electrodes placed in frontal areas. However, at present m r

o 1 we have no explanation for pattern of changes taking place in N 0 10 20 30 40 50 60 Time(s) Before BP coherence values. Theta wave P After BP

S R 0.9 d e

z 0.8 i l IV. DISCUSSION AND RESEARCH ISSUES a 0.7 m r

o

N 0 10 20 30 40 50 60 The results presented in this paper on BP are from two subjects Time(s) Before BP so statistical reliability may be in questions. But this is our Alpha wave P After BP precursor paper on the explored topic and trends found in S 1 R

d results show some interesting things for further deep e

z 0.5 i l exploration of the matter. The effects of BP on theta and a m r

o 0 gamma brainwave are important. It is also important to note

N 0 10 20 30 40 50 60 Time(s) that the second subject of our experiment does not has any Before BP Beta wave

P practice of BP before EEG recording, however, observed

S 1 After BP R

effects on both subjects show some similar trends. This shows d e

z 0.5 i l a m r

o 0 N 0 10 20 30 40 50 60 Delta wave

Time(s) P 1 Before BP S R

Gamma wave P d

S After BP 1 e Before BP z R 0.9 i

l d a After BP e m z 0.5 i r l o a 0.8 N m

r 0 10 20 30 40 50 60

o 0

N 0 10 20 30 40 50 60 Time(s) Time(s) Theta wave

P 1 Before BP S R

After BP d e

z 0.5 i l a

Fig 6 Showing temporal variation of RSP in delta, theta, m r o 0 alpha, beta and gamma brainwave pattern before and after N 0 10 20 30 40 50 60 BP for the first subject.(Data used from the right side of Time(s)

the head ) Alpha wave

P 1 Before BP S

R After BP

d e

z 0.5 i l wave is different in both subjects. In the first subject it is a m r decreasing (in lower range) after BP while in the second subject o 0 N 0 10 20 30 40 50 60 it is increasing after BP. So it is difficult to make preliminary Time(s) comments on such trends with observations on only two Beta wave subjects. The gamma wave patterns are high frequency P 1 Before BP S

R After BP brainwaves and are found in the every part of the brain. It has d e

z 0.5 i been hypothesized that the gamma brainwaves are related with l a m the performance of higher mental activities and perceptual r o 0 tasks. Such tasks need functional integration and N 0 10 20 30 40 50 60 synchronization of different parts of brain which is established Time(s) by gamma brainwaves. Both subjects show improved power in Gamma wave P 1 Before BP S

the gamma brainwave patterns after BP. Difference is clearer in R After BP

d e the second subject. Significant change in gamma wave pattern z 0.5 i l a

by BP is very important for future research. There has also been m r o 0 reporting in the past on interaction of auditory stimuli and N 0 10 20 30 40 50 60 gamma brainwaves in many subjects [23], [24]. Time(s) As per Eq. (2), the coherence function for each sensor pair was calculated. The coherence matrices for the first subject for EEG data at each of 21 sensors are shown in Fig. 9 for the both (before BP and after BP) conditions. For the second subject Fig.7 Showing temporal variation of RSP in delta, theta, alpha, coherence matrices for the both conditions are shown in Fig.10 beta and gamma brainwave pattern before and after BP for the and Fig.11, respectively, in which the self coherence values for second subject.(Data used from the right side of the head ) all the sensors were set to zero, which were really 1 lying on

- 393 -

Delta wave P

S 1 R

d e

z 0.9 i l

a Before BP m r After BP

o 0.8

N 0 10 20 30 40 50 60 Time(s) Theta wave P Before BP

S 1 R

After BP d e

z 0.5 i l a m r

o 0

N 0 10 20 30 40 50 60 Time(s) Alpha wave P Before BP

S 1 R

After BP d e

z 0.5 i l a

m r

o 0

N 0 10 20 30 40 50 60 Fig.10 Coherence matrix for the second subject before BP. Time(s) Beta wave

P Before BP

S 1

R After BP

d e

z 0.5 i l a m r

o 0

N 0 10 20 30 40 50 60 Time(s) Gamma wave P Before BP

S 1 R

After BP d e

z 0.5 i l a m r

o 0

N 0 10 20 30 40 50 60 Time(s)

Fig.8 Showing temporal variation of RSP in delta, theta, alpha, beta and gamma brainwave pattern before and after BP for the second subject.(Data used from the left side of the head )

Fig.11. Coherence matrix for the second subject after BP

BP is capable in bringing effect immediately. One of the most active components of BP is humming sound which we suppose playing key role in producing results. The brainwave phenomena and sound wave are different things but their interaction is very common experience in daily life, e.g. a child sleeps on hearing lullaby by his/her mother. The interaction of sound wave and brainwave has been discovered in past [25] and recently it has gained clinical importance, leading to business, for bringing relaxation by artificially generated external sounds. However, in BP sound is internal and subject's own sound so it's important to investigate interaction of own sound with change in brainwave patterns. Since the results presented in this paper are only from two subjects it's essential to do experiments with large number of subjects to get statistically acceptable conclusions and we left it for future studies. The dynamics of brainwave changes during the BP will help in Fig.9 Coherence matrix for the first subject before and after understanding mechanism of BP by which it works to bring BP. positive effects. Such explorations will be useful in developing

- 394 -

some biomimetic system for the challenged people. It will be interesting to carry out analysis in this direction.

REFERENCES [1] I. Whicher, D. Carpenter (Eds.)," Yoga the Indian tradition", Rutledge Curzon, hardback, 2003. [2] P. Sinha, Yogic cure for common diseases, Orient Paperback, New Delhi 1996. [3] B.R.Cahn, J. Polich, "Meditation States and traits: EEG, ERP and nuroimaging studies", Psychological Bulletin, Vol.132, No.2, 180-211, 2006. [4] C. Wood, "Mood change and of vitality: a comparison of effects of relaxation, visualization and yoga", Jr. Royal Society of medicine, vol.86, 1993. [5] Bhala, Balmakund, “Yoga in modern medicine”, International Journal of Alternative and Complementary Medicine, April, 1996. [6] Svatmarama, "The Pradipika" Translated into English by Brian Dana Akers. Woodstock: YogaVidya.com, 2002. (http://www.conjurer.org/library/yoga/hyp.pdf) [7] J.Rodenbeck, "Bhramari the bumble bee', Bindu, no.10,(http://www.scand-yoga.org/english/ bindu/ bindu 10/index.html) [8] M.Grady, "Sound of bee', Available at http://www. himalayaninstitute.org/ArticleRead.aspx?code=138 [9] R.K.Prasad, F. Matsuno,' Bhramari pranayama: An excellent yoga to use your own voice in healing", Accepted for ICPTJ, Tokyo, 2006. [10] Jella SA, Shannahoff-Khalsa DS," The effects of unilateral nostril breathing on cognitive performance. Int J Neurosci ; 73(1-2): 61-8, Nov.,1993 [11] Wood C. Mood change and perceptions of vitality: a comparison of the effects of relaxation, visualisation and yoga. J R Soc Med 1993 May; 86(5): 254-8. [12] Telles S, Nagarathna R, Nagendra HR. Breathing through a particular nostril can alter metabolism and autonomic activities. Indian J Physiol Pharmacol (India) 1994 Apr; 38(2): 133-7. [13] Bhatia, M., Kumar, A., Kumar, N., Pandey, R.M., and Kochupillai, V. Electrophysiologic evaluation of Sudarshan Kriya: an EEG, BAER, and study. Indian J. Physiol. Pharmacol. 47, 157-163, 2003. [14] A. Cichocki, S. Amari, Adaptive blind signal and image processing, New York:, Wiley, 2003. [15] A. Cichocki, S. Amari, et.al., "ICALAB toolboxes", available at http://www.bsp.brain.riken.jp/ICALAB/ [16] Banquet, J. P. (1972). EEG and meditation. and Clinical Neurophysiology, 33, 454. [17] Banquet, J. P. (1973). Spectral analysis of the EEG in meditation. Electroencephalography and Clinical Neurophysiology,35, 143–151. [18] G.D.Jacobe, R.Friedman, "EEG spectral analysis of relaxation techniques' Applied Psychophysiology and Biofeedback, Vol. 29, No. 4, December 2004 (_C 2004) [19] Canteros, J. L., Atienza, M., Stickgold, R., & Hobson, J. A. Nightcap: A reliable system for determining sleep-onset latency. Sleep, 25, 238–245, 2002. [20] West, M. A. Meditation and the EEG. Psychological Medicine, 10, 369–375, 1980. [21] Jacobs, G. D., & Lubar, J. F. (1989). Spectral analysis of the central nervous system effects of the relaxation response elicited by autogenic training. Behavioral Medicine, 15, 125–132. [22] Plotkin, W B., & Rice, K. M. (1981). Biofeedback as a placebo: Anxiety reduction facilitated by training in either suppression or enhancement of alpha brainwaves. Journal of Consulting and Clinical Psychology, 49, 590-96. [23] H. Titenen, et.al., "Selective enhances the auditory 40-Hz transient response in humans", Nature, Vol.364, July,1993. [24] C. Pantev, et.al., " Human auditory evoked gamma band magnetic fields", Proc. Natl. Acad. Sci, USA, Vol.88, pp.8996-9000, Octb,1991. [25] R. A. Monroe," Hemi-Sync method of inducing mental, emotional and physical states of consciousness, including specific mental activity, in human beings", US Patent Number: 5213562; Issue Year: 1993;

- 395 -