Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 8249.
Published online 2017 August 15. doi:  10.1038/s41598-017-08421-z
PMCID: PMC5557761

Fine Structure of Posterior Alpha Rhythm in Human EEG: Frequency Components, Their Cortical Sources, and Temporal Behavior


Heterogeneity of the posterior alpha rhythm (AR) is a widely assumed but rarely tested phenomenon. We decomposed the posterior AR in the cortical source space with a 3-way PARAFAC technique, taking into account the spatial, frequency, and temporal aspects of mid-density EEG. We found a multicomponent AR structure in 90% of a group of 29 healthy adults. The typical resting-state structure consisted of a high-frequency occipito-parietal component of the AR (ARC1) and a low-frequency occipito-temporal component (ARC2), characterized by individual dynamics in time. In a few cases, we found a 3-component structure, with two ARC1s and one ARC2. The AR structures were stable in their frequency and spatial features over weeks to months, thus representing individual EEG alpha phenotypes. Cortical topography, individual stability, and similarity to the primate AR organization link ARC1 to the dorsal visual stream and ARC2 to the ventral one. Understanding how many and what kind of posterior AR components contribute to the EEG is essential for clinical neuroscience as an objective basis for AR segmentation and for interpreting AR dynamics under various conditions, both normal and pathological, which can selectively affect individual components.


The most prominent activity recorded with electroencephalography (EEG) from the scalp of a healthy human at rest is the near-sinusoidal oscillations within the 7–13 Hz range called the alpha rhythm (AR). Multi-rhythmicity of this frequency band has been known for a long time1, 2. Empirical evidence shows that the waking AR consists of several oscillatory components, which differ in scalp topography and reactivity to sensory stimulation or to motor activity. The family of AR includes the central mu rhythm, the tau rhythm, recorded over the temporal cortex, the anterior AR, and the widespread and ample posterior AR2, 3. The latter, also termed the occipital or visual AR, is either considered a unitary phenomenon46 or divided into 2 or 3 sub-bands711.

By and large, the authors of EEG studies assume the heterogeneity of the posterior AR as a well-established fact, based on which they demonstrate individual properties of preselected AR sub-bands under various physiological and pathological conditions. A non-exhaustive list of supporting evidence includes: different responses of the low- and high-frequency ARs to visual stimulation1113; frequency-specific AR correlates to memory and attention14, 15; the selective efficiency of the high-frequency AR in neurofeedback training16; and separate alterations of the low- and high-frequency ARs due to neurodegenerative processes in old age8, 17.

Although the above findings suggest that the surface-recorded posterior AR represents a combination of rhythms, their number as well as their cortical sources require a systematic investigation, which would involve an analytical determination of AR components (ARCs) based on EEG spectra obtained with a high frequency resolution. Along these lines, Chiang and colleagues18 studied the AR structure in 1500 healthy subjects by means of cluster analysis. Of all participants, 44% showed a two-cluster AR, while in the rest of the population, multiple peaks remained unresolved. Another study19 based on a similar approach identified multiple ARCs only in 1.6% of 1215 subjects. The inability to separate ARCs in many EEGs can be attributed to the methods that took into account the visibility of separate spectral peaks, but underestimated the spatial and temporal dynamics of the AR because of crude spatial sampling (19 electrodes) and short EEG recordings (2 min).

Alternatively, factor analysis was applied to the EEG power of healthy subjects20 and patients21. In both studies, 2 independent ARCs, which, in general, corresponded to empirically determined low- and high-frequency ARs, were reported. However, being a bilinear decomposition technique, factor analysis cannot fully account for the multi-way nature of EEG data. Specifically, the spatial aspect of AR was ignored since only a few preselected EEG derivations were used. Moreover, these attempts of AR decomposition were based on short recordings of low-density referential EEG, thus leaving beyond the scope of analysis the sources of ARCs, as well as the extent of their independence, stability, and dynamics at various time scales.

This motivated us to test the hypothesis that the posterior AR in healthy young adults has a structure characterized by separate components with individual spatial, frequency, and temporal features with advanced techniques involving a multi-linear decomposition of mid-density EEG. Considering the assumed role of the AR in the mechanisms of cognition and its dramatic changes during human development and aging, detailed knowledge of AR component structure can help us better understand the role of oscillatory activity in information processing and develop efficient biomarkers of the brain functional state.



Thirty-three healthy subjects between 20 and 45 years old were recruited for this study. Here we report the data from 29 of the subjects, since EEG data from 4 subjects were excluded from the analysis owing to the excessive muscle artifacts (1 EEG), no visible AR (2 EEGs), and the interrupted recording session (1 EEG). The group consisted of 17 men and 12 women (mean age 31.6). Prospective participants underwent a brief interview, which included the MoCA (Montreal Cognitive Assessment) test used to control for cognitive status22. Only individuals with MoCA score ≥26 (mean group score was 28.4) and without cognitive complaints, psychoactive drug use, head traumas or other conditions that interfere with cognition, were enrolled in the project.

All methods and procedures in this study conform to the Declaration of Helsinki (1964) of the World Medical Association concerning human experimentation. It was approved by the local Ethics Committee (Commission cantonale d′éthique de la recherche sur l′être humain (CER-VD)). To obtain written informed consent, we provided the essential information about the research to each participant.

EEG recording and pre-processing

The experimental protocol started with interviewing a subject (see Subjects), followed by placing and fitting the EEG cap on the subject’s head. During the recording session (1.5–2.0 hours), subjects were seated in a comfortable chair and instructed to stay awake with eyes closed (Resting state with Eyes Closed, REC). The four REC segments (REC1, REC2, REC3, REC4) alternated with the resting state with eyes open or else with the hand movement task and 2 haptic cognitive tasks. To keep the subjects in a stable awake state over the entire recording session, we alternated the task and REC conditions such that an uninterrupted REC duration could not exceed 3 minutes. Here we report only REC data, which consist of 4 concatenated periods about 8–10 minutes each.

The EEG data were collected by means of mobile electroencephalograph eegosports using a 64-channel waveguardTM original cap (ANT Neuro b.v., Enschede, The Netherlands). The recordings were made with a CPz reference and a high-cutoff filter set to 100 Hz. The electrode impedances were kept under 30 kΩ (10 kΩ in most cases), which was well below the recommended maximum of 50 kΩ for high-impedance eego amplifiers.

The signals were digitized at a rate of 500 samples/s with a 24-bit analog-to-digital converter. They were band-pass filtered from 1 to 45 Hz using a zero-phase filter and edited offline. The artifactual segments were identified by a thorough visual inspection and removed. The signals from mastoid channels were removed from all recordings because of excessive noise in many EEGs. For further analysis, the EEG signals from the remaining 62 EEG channels were re-referenced against the common-average reference.

General design of EEG analysis

A straightforward approach to the analysis of ARCs would be a decomposition of the source-space EEG spectrum by a 3-way parallel factor analysis (PARAFAC), which would produce information about frequency, spatial, and temporal features of ARCs. The main steps would be: calculating the sensor-space EEG spectrum, estimating the source spectrum and, finally, applying PARAFAC to the source spectrum. However, since PARAFAC is more robust for processing data with a higher signal-to-noise ratio and a smaller number of parameters to estimate, we first applied the method to the sensor-space EEG spectrum (PARAFAC1). That allowed us to obtain reliable estimates of the frequency and temporal features of ARCs and to use these estimates as fixed parameters in the application of PARAFAC to the amplitude spectrum of the signals in the source space (PARAFAC2). Figure 1 shows the actual pipeline of our analysis: (1) calculating the EEG spectrum, (2) applying PARAFAC in the sensor space, (3) estimating the source EEG spectrum, and (4) applying PARAFAC in the source space. The following paragraphs provide the necessary details of all the analysis steps.

Figure 1
Extraction of alpha rhythm components (ARCs) from EEG signal. The top row shows 4 EEG segments from corresponding REC periods in a single subject. To extract the time-varying EEG power spectrum matrix, Short-Time Fourier Transformation (STFT) was applied ...

Spectral analysis of EEG

We applied a time-frequency analysis to the concatenated EEG recordings from 4 REC periods (Fig. 1) using discrete Short-Time Fourier Transform (STFT) of the time-windowed data. If we consider EEGs as a matrix D ∈ RNm×Nτ, where R refers to real numbers, N m is the number of channels, and N τ is the number of time points over all available REC periods, then the STFT of EEG at a channel m is calculated as


where W(τ  t) is a windowing function of length L centered at time t, N f is the number of discrete frequencies, N t is the number of time windows, and C denotes complex numbers. The choice of 5 second Hann windows W with 50% overlap resulted in N t = 781 ± 100 (mean ± standard deviation) per subject, a 0.2 Hz frequency resolution of the spectral data (df), and N f = 51 for the 5–15 Hz extended alpha range that was considered in our analysis.

The cross-spectrum density matrices of scalp EEGs were calculated as


where ϕf,t ∈ CNm×1 is the STFT of all sensors at time t and frequency f, and 2(W(τ)2¯) is the cross-spectrum normalization factor, with W(τ)2¯ being the average of window value squared (for the Hann window it is 0.375).

The EEG amplitude spectrum (amplitude spectral density, ASD) at frequency f and time window t was calculated as the square root of the diagonal of CSD[var phi] f,t. The 3-way EEG amplitude spectrum matrix ASDϕ ∈ RNm×Nf×Nt, with frequency, channel, and time dimensions served as the input to the first PARAFAC analysis.

Source analysis of EEG

To demonstrate the source-space distribution of ARCs, we used the LOw-Resolution Electromagnetic Tomography (LORETA)23. The lead field matrix was computed using the 3-shell Locally Spherical Model with Anatomical Constraints (LSMAC) method for 3005 uniformly distributed sources located in the gray matter of the MNI (Montreal Neurological Institute) average brain. The resulting inverse matrix T ∈ RNs×Nm (where N s is the number of sources) was down-sampled to 387 uniformly distributed sources in order to reduce the computation time. Using the inverse matrix and sensor-space cross-spectrum density CSD[var phi] f,t the source cross-spectra were computed as

CSDjf,tT . CSDϕf,t . TTCSDjf,t ∈ CNs×Ns.

The source amplitude spectrum ASDj ∈ RNs×Nf×Nt was obtained as the square root of the diagonal of CSDj.

Parallel factor analysis

In order to extract the ARCs with their specific frequency, spatial, and temporal features, we employed the PARAFAC method. It is a method for decomposing multi-way data into a small number of hidden variables called components or factors, which explain most of the variance of the original data while reducing its overall dimensionality. Among the multi-way decomposition methods, PARAFAC represents the simplest and the most restricted generalization of bilinear Principle Component Analysis (PCA) to higher order data2426. In contrast to PCA, PARAFAC provides a unique solution, provided the number of components is correctly chosen and the signal-to-noise ratio is appropriate. The method was validated by Bro et al.24, 27 using simulations and empirical data with known underlying structure. It was applied to EEG for the first time to decompose the topographic components of event related potentials28, 29. In more recent applications, PARAFAC was used for the identification of seizure foci in ictal EEG of epileptic patients30, 31, the “atomic” decomposition of normal EEG32, and for separating the EEG sources on the basis of their phase coupling patterns that vary in space, frequency, and time33.

The latter applications of PARAFAC in the EEG area were essentially aimed at introducing this technique as an instrument of EEG decomposition32 or proposing an extension of the method33 rather than addressing a particular neuroscience problem. The validation of PARAFAC with simulations and its comparisons with other decomposition methods have been reviewed there in detail, and we refer the reader to these reports for an in-depth analysis of these issues31, 32, 34. Here we use PARAFAC to show the fine structure of the AR as an individual EEG phenotype. To this end, PARAFAC application is supplemented by a combination of the experimental design features and analysis techniques tuned to the requirements of testing the narrow-band EEG activity, including a high frequency resolution, EEG source reconstruction, long-term and longitudinal EEG recordings.

PARAFAC decomposes multi-way data into a sum of components, which result from a multi-linear production of loading and score vectors. In this study we applied a 3-way (frequency, time, space) PARAFAC decomposition of the EEG spectrum with non-negativity constraints on the loadings to keep the results interpretable. Since in a 3-way PARAFAC analysis loading and score vectors are treated equally, they are not distinguished. To formalize, PARAFAC decomposes the 3-way data matrix ASD[var phi] into K components as follows:


Here, the loading matrices A, B, and C with elements a fk, b mk, and c tk, correspond to the frequency, space (channels or sources), and time modes of the original data. Intuitively, each matrix contributes a “signature” of its corresponding mode to each of the K components.

A graphical illustration of this model is presented in Fig. 1. In the PARAFAC analysis, the solution is obtained by the alternating least squares method. This iterative technique estimates the parameters of each mode while assuming that the parameters of two other modes are known. It iterates until convergence. The second and third loading vectors are normalized so that the variance of the data is kept in the first mode. The obtained loading matrices have dimensionless values. To implement PARAFAC, we used the N-way toolbox35 in Matlab.

As a result of the first PARAFAC analysis (PARAFAC1 in Fig. 1), for each subject the ASD[var phi] in the sensor space was decomposed into a combination of ARCs with specific frequency, channel topography, and time signatures. As an input to the second PARAFAC analysis, we used the 3-way source amplitude spectrum density matrix ASDj ∈ RNf×Ns×Nt with the frequency, source, and time dimensions. The two loading matrices corresponding to the frequency and temporal features had the fixed values calculated in the sensor-space PARAFAC1. Therefore, the source-space model estimated the locations of ARC sources using their predefined frequency and temporal signatures. The source locations were then parcellated according to the Brodmann atlas using the Talairach Daemon software36.

PARAFAC Model verification

To develop a reliable model, we checked the residuals and averages, applied degeneracy tests for the stability and reliability of the model, and performed Core Consistency Diagnostic (CORCONDIA)27, 29 using the N-way toolbox. Selection of the appropriate number of components is a critical issue in developing a valid PARAFAC model. Selecting too few components results in a model that is unable to identify the true underlying variables. In contrast, too many components would excessively model noise and/or the true underlying variables would be modeled by correlated components. A number of approaches have been proposed for the selection of a proper number of components including cross-validation, judging residuals, and CORCONDIA24, 27. In the preliminary analysis, we used both cross-validation and CORCONDIA, which produced similar results. Ultimately, we decided on CORCONDIA, which has been shown to provide appropriate results in simulations and real datasets with known structure27 by testing the model complexity.

The approach comes down to decomposing the data with increasing number of components and calculating the core consistency value, which varies between the maximum of 100% in case of perfect model and 0 or below when the model is invalid. A model with the largest number of components and its core consistency exceeding some threshold is considered as valid. Here, following a proposal by Bro27, we set the threshold at 90%.

To remove the noise components, we focused on their frequency and spatial characteristics. The components with flat spectra and non-posterior topography (e.g., those with highest loadings on frontal and anterior temporal regions sensitive to eye- and muscle-related artifacts) were removed from the model. Although this report is focused on the posterior AR, the analysis of whole-head EEG suggested a possibility to identify other ARCs, for instance, the central mu rhythm. For those who are interested in the latter, we provide relevant results in the Supplementary materials.

Group-level ARC statistics

For each of the 26 participants whose AR was successfully decomposed into at least two ARCs, the highest frequency component was labeled ARC1 and the lowest frequency one ARC2. To compare the source distributions of these two components, we estimated the significant differences of the spatial loadings in the contrasts ARC1 > ARC2 and ARC2 > ARC1 by means of the one-sided cluster-based permutation test37. This non-parametric method extracts the largest significant cluster of sources using permutation test and cluster-based statistics to correct for multiple comparisons.

For each subject, the loading vectors of the source-space PARAFAC2 components were normalized so that their maximum values were set to 1, and, for each source, we computed the one-sided t-statistics by comparing the normalized values of the two components. Then, based on spatial adjacency, we selected the clusters of sources with t-values above the 99th percentile of the t-distribution, and calculated cluster-level statistics as the sums of the t-statistics for all sources within each cluster. The ARC labels were permuted 5000 times, and, for each permutation, the largest cluster-level statistic was selected, thus providing a random distribution of the largest cluster-level statistics. Finally, the P-value for the largest cluster in the non-permuted data was computed as the number of cluster-level statistics in the random distribution greater than the non-permuted cluster-level statistic divided by the number of permutations.

Analysis of ARC replicability

To evaluate the temporal stability of the AR structure over a single recording session (1.5–2 hours), for the REC1, REC2, REC3, and REC4 periods of all the subjects with 2 or 3 ARCs, we computed the mean temporal loadings of ARC1 and ARC2 normalized so that the maximum loading was set to 1. To check for possible task effects on the structure of the resting-state AR, we applied a group-level 1-way repeated-measures ANOVA with task as a factor (4 levels) and temporal loadings of an ARC averaged over each of the REC periods as a dependent variable.

To estimate the source and frequency stability of the AR structure at a time scale of 2 to 13 months, we compared the concatenated data from four REC segments between two separate sessions in 5 subjects available for follow-up. The similarity between loading vectors of ARCs obtained from two different EEG sessions was calculated using Tucker Congruence Coefficients (TCC) as


where x k and y k are the elements of loading vectors of a component extracted from the two PARAFAC models. TCC for the non-negative loadings varies between 0 and 1. The TCC value of >0.85 shows highly similar loadings and the TCC value of >0.95 points to nearly identical loadings. We compared loading vectors from spatial and frequency modes.


Decomposition of individual AR with PARAFAC

The individual amplitude spectra of the common-average-referenced resting-state EEGs suggested that in many subjects more than one rhythm in the alpha range could be recorded from the posterior sensors (Figs 2A, ,3A3A and and4A).4A). Indeed, the decomposition with the PARAFAC method showed that 26 out of 29 subjects had more than one ARC. The typical structure of the AR band, observed in 24 participants (83%) involved 2 ARCs (Fig. 2B). In two thirds of these subjects, the source distributions of the two ARCs were markedly different (Fig. 2C, subject 16; Fig. 4C, subject 4), while in the remaining third the two ARCs had significant overlap of their spectral signatures in the alpha band and similar source distributions (Fig. 2C, subject 22).

Figure 2
Individual examples of two-component structure of alpha rhythm. (A) The amplitude spectrum densities (ASD) of an EEG from posterior electrodes (O1, P1, P5, and TP7 according to the Extended 10/20 System) are presented in the conventional alpha frequency ...
Figure 3
Individual examples of three- and one-component structures of alpha rhythm. All designations are as in Fig. 2. See Results for details.
Figure 4
Follow-up illustrates individually stable multi-component structure of alpha rhythm. The representative examples show the component structure of the alpha rhythm in subject 4 with two ARCs and subject 49 with three ARCs, based on two EEG recordings separated ...

The sources of the higher-frequency ARC (ARC1), with its peak frequency around 10 Hz, dominated in the occipito-parietal regions of the neocortex, while the sources of the lower-frequency (around 9 Hz) ARC2 were dominant in the occipito-temporal regions. The maximum loading values of ARC1 were located in the dorsal BA18/19 close to the midline or in the neighboring BA20 in all but 2 subjects, in whom they were localized to BA37. The maxima of ARC2 were located in BA37 (fusiform gyrus) or in the neighboring sites of BA19 or BA20 in all but one subject, in whom the maximum was found in BA18. This 2-component structure was present in both sexes and across the whole age spectrum of our sample (20–45 years).

In 2 male subjects, we found a more complex AR structure represented by 3 ARCs (Fig. 3, Subject 52; Fig. 4, Subject 49). These included the low-frequency occipito-temporal ARC (similar to ARC2 in the 2-component AR) and the two high-frequency ARCs (peak frequency ≥10 Hz) with predominantly occipito-parietal topography, closely resembling ARC1. These 3 ARCs manifested co-existence and individual behavior over the studied temporal intervals (Fig. 3D). The ARs of the remaining 3 subjects could not be decomposed by our analysis (Fig. 3, Subject 42). In these cases, the single ARC had a peak frequency of at least 10 Hz and originated from posterior sources that varied among the subjects, making it difficult to identify the component as any of the previously described ARCs.

Individual ARC replicability

Throughout a single recording session, the AR structure was relatively stable in all the subjects with 2 or 3 ARCs (Figs 2D, ,3D3D and and4D).4D). Specifically, in this group, the mean temporal loadings for ARC1 in 4 REC periods were 0.35–0.38 and for ARC2, 0.37–0.40. The repeated-measures ANOVA showed no significant modulation of the temporal loadings of ARC1 or ARC2 by the tasks.

To estimate the stability of the AR structure on a large time scale, we compared spatial and frequency loadings from baseline and follow-up EEG recordings of 5 subjects. In all these subjects, ARC1 showed stability of its spatial and frequency features over the two recordings (TCC > 0.85). The spatial loadings of ARC2 also remained stable in all the subjects. In contrast, the frequency loadings of ARC2 were very similar only in 3 subjects, whereas 2 others showed less similarity (0.75 < TCC < 0.80). In subject 49 (Fig. 4) with a 3-component AR structure, the spatial and frequency loadings of the 3rd component were practically equal in the two recordings (TCC > 0.94).

Group-level analysis of AR components

Since the posterior AR reveals at least 2 distinct components in 90% of our subjects, one can give a group-level description of the 2-component structure of AR. Specifically, ARC1 has a peak frequency of 10.4 ± 0.6 Hz (mean ± standard deviation) and shows higher loading values (ARC1 > ARC2, P = 0.006) in a large (~16 cm3) bilateral cluster (Fig. 5) that includes primary and association visual, parietal, and posterior cingulate cortices (BA7, 17, 18, 19, 23). In contrast, ARC2 with a peak frequency of 9.4 ± 0.7 Hz has significantly higher loading values (ARC2 > ARC1, P = 0.041) for the left-hemisphere cluster (~7 cm3) mainly distributed in the inferior temporal and fusiform gyri (BA20, 36, 37), whereas the right-hemisphere cluster did not survive corrections for multiple comparisons with the cluster-based permutation test.

Figure 5
Group level statistics shows different sources of alpha rhythm components. Top row: Normalized group-averaged spatial loadings of ARC1 rendered on the average MNI brain are shown in four views. Second row: Normalized group-averaged spatial loadings of ...


As outlined in the Introduction, heterogeneity of the posterior AR has been assumed but rarely explicitly tested, and only with two-dimensional decomposition techniques, which left the ARCs unresolved in one of the three EEG dimensions20, 21.

Here, with the multidimensional PARAFAC technique, which takes into account the spatial, frequency, and temporal properties of EEG, we demonstrate that a typical structure of the resting-state posterior AR consists of 2 components, originating from the occipito-parietal (ARC1) and occipito-temporal (ARC2) cortices.

Previous EEG/MEG studies consistently localized the equivalent current dipoles (ECDs) of the posterior AR along the brain midline in the parieto-occipital and, to a lesser extent, in the calcarine sulci3843. The EEG/MEG signals were band-pass filtered and an ECD was computed every several milliseconds, resulting in cortical clusters of ECDs. Band-pass filtering of raw EEG/MEG signals within a frequency range spanning the entire alpha band makes problematic the resolution of sources co-existing at neighboring frequencies. This explains why the dominant occipito-parietal component was overemphasized in the literature, while the weaker occipito-temporal component remained under-reported.

The rare reports that addressed the component structure of the posterior AR either did not localize the ARC sources20, 21, 44 or localized one of the components to the temporal rather than the occipito-temporal region45. Srinivasan and colleagues showed that in MEG both occipital and temporal sensors detect strong ARs, while in the EEG of the same subjects temporal sensors record a much weaker AR than occipital ones. Besides, the temporal AR in MEG had lower peak frequencies than the occipital AR, while this difference was not obvious in EEG. These findings share important features with our results: the occipito-temporal ARC2 reported here is usually weaker and has a lower peak frequency than the occipito-parietal ARC1, and their separation without a decomposition technique is challenging (Figs 2A, ,3A3A and and4A).4A). However, in Srinivasan’s study, the sources of the temporal AR were only broadly defined by the sensor locations covering the entire convexity of the temporal lobe, making it unclear to what extent the EEG-based ARC2 that we describe is homologous to the low-frequency AR recorded by MEG.

The recently reported decomposition of MEG activity by Independent Component Analysis in the cortical source space revealed occipitо-parietal 10 Hz components in all the subjects, and temporal 8–10 Hz components in some of them45. The sources of the latter were localized to the superior temporal cortex, suggesting that these components have a different origin from the EEG-based ARC2, which we localized to the fusiform gyrus and neighboring sites. Considering that the temporal AR was mainly recorded with MEG, localized to the midtemporal region, and never reported in the majority of the tested subjects, we conclude that, in previous EEG studies, the occipito-temporal ARC2 was overlooked or not accurately differentiated from the temporal auditory AR.

A prerequisite for AR “visibility” in EEG are layers of radial dipoles occupying cortical territories from several square millimeters to several square centimeters that generate coordinated rhythmic activity1. Classical intracranial human studies showed that a strong continual AR occurs over the convexity of the parietal, posterior temporal, and occipital cortices of an exposed brain46, 47, leaving the issue of its possible heterogeneity unresolved. This question was examined in a recent electrocorticographic study of epileptic patients, which found two widespread posterior clusters of spontaneous oscillations with dissimilar spatial and frequency features48. A relatively narrow cluster peaking at 10 Hz was limited to the occipital and parietal regions, whereas a broad cluster peaking at 7 Hz was mainly located in the occipital and inferior temporal cortices. Considering the EEG slowing in epilepsy49, 50, the broad occipito-temporal cluster can be tentatively matched to ARC2.

Spontaneous alpha-range oscillations with their peak frequencies and functional properties distinct in the ocipital and inferotemporal cortices were also recorded in awake macaques51, 52. In the occipital areas, pyramids of layer 5 were driving the oscillations, whereas in the inferotemporal cortex, pyramids of layers 2/3 served as a pacemaker (ibid). Since human BA37/19 (lateral occipital complex), where we localized the low-frequency ARC2, is a homolog of the macaque inferotemporal cortex53, it is conceivable that, in humans, the laminar mechanisms of the occipito-temporal ARC may also differ from occipito-parietal ones. The fact that the occipito-temporal ARC2 is less powerful than the occipito-parieatal ARC1 favors this suggestion: based on the animal model, a smaller contribution of ARC2 to the surface EEG can be expected due to the tendency of radially distributed dendrites, inherent in the temporal cortices, to form a closed field that falls off faster than a radial dipole configuration typical for the occipital cortex54, 55.

In the two subjects with a three-component AR structure, the low-frequency ARC2 was found in the occipito-temporal cortices, whereas the two high-frequency ARCs were localized to the occipito-parietal regions. This finding closely resembles observations by40, 43, 56. In43, the rhythms peaking around 10 Hz formed two clusters, with the stronger one located in the parietal region and the weaker one in the occipital areas. Here, with the extremely limited number of subjects exhibiting this AR structure, the results of the source reconstruction did not allow us to separate the sources of these high-frequency ARCs. The uncommonness of this configuration in a resting state does not rule out its higher incidence under other conditions. Indeed, as many as three posterior clusters of ARCs with distinct spectral peaks were identified in the subjects involved in a visuo-spatial task57.

Repeated observations of the AR show that its frequency and spatial characteristics are subject-specific and manifest the highest heritability among conventional EEG frequency bands46, 58. These studies considered the AR as a unitary rhythm and, therefore, were biased in favor of the dominant high-frequency occipito-parietal component. To our best knowledge, our study is the first report on the replicability of separate ARCs over weeks or months. Importantly, we estimated not only frequency features of the ARCs, but also their cortical topography. The source distributions of both components as well as the frequency of ARC1 remained highly stable, whereas the ARC2 frequency showed somewhat higher variability, suggesting state-dependent effects on this component.

While the long-term stability of the AR structure remains to be confirmed in adults, the AR “signature”, i.e., its unique component structure, has a good chance to be an individual characteristic, as suggested by an early longitudinal study of healthy children59. These authors examined EEG spectra with a high frequency resolution of 0.1 Hz that allowed the differentiation of the posterior AR into 2 components with distinct topographies, peak frequencies, and reactivities. These individual AR structures, found in the EEGs of 7-year-old children, proved to be well-preserved in the follow-up recordings 4 years later.

Implications for human neuroscience

According to our analysis, in healthy adults the posterior resting-state AR typically consists of two distinct, occipito-parietal and occipito-temporal, components. The cortical topography, individual stability, and similarity to the primate alpha-band organization of this 2-component AR structure link it to the functional modules of the posterior brain, specifically, to the dorsal (ARC1) and ventral (ARC2) visual streams. The next step would be an investigation of functional properties of the two ARCs in paradigms involving tasks that differentially engage the dorsal and ventral visual pathways.

Understanding how many rhythms in the alpha range contribute to the posterior surface EEG is essential for clinical neuroscience as an objective basis for the AR segmentation in the frequency domain. Specifically, the uncommonness of mono- and three-rhythmic structure in the resting-state EEG limits segmentation options in cross-sectional studies. The mechanical (not data-driven) subdivision of the AR into a priori defined 3 sub-bands, often encountered in the literature, might provide misleading effects, while treating the AR as a single component (also frequently encountered) results in a mixture of effects originating from the components with different properties. It should be emphasized that unmixing the AR components with the PARAFAC technique per se does not require a source analysis, i.e., can be used as a decomposition technique in various clinically motivated analyses of surface EEG.

Finally, the interpretation of AR dynamics under various psychiatric and neurodegenerative conditions requires revision in view of the 2-component model, since apparently similar changes can result from different pathological processes selectively affecting AR components. For instance, the posterior AR often “slows down” in patients with obsessive-compulsive disorder, depression, and Alzheimer’s disease6062. Possible scenarios leading to such changes include a reduction in high-frequency ARC power, an increase in low-frequency ARC power, a decrease in frequency of both components, etc. To learn what really occurs, the study of differential dynamics of the alpha rhythm components is essential.

Data availability statement

The datasets analyzed during the current study are available from the corresponding author on reasonable request.

Electronic supplementary material


This work was supported by Swiss National Science Foundation and Swiss Research Program Nano-Tera grant to Maria G. Knyazeva (BodyPoweredSenSE project). The authors would like to thank Prof. J.-F. Demonet, Director, Leenaards Memory Centre, Hospitalier Universitaire Vaudois, for permanent support, helpful comments, and valuable discussion regarding the research and the manuscript. The authors are grateful to all the participants for their willingness to freely contribute time and effort to this study.

Author Contributions

Author Contributions

The authors meet the criteria for authorship recommended by the International Committee of Medical Journal Editors and take full responsibility for all contents of the manuscript. All authors (E.B., V.V., and M.K.) contributed to the study design. E.B. and M.K. implemented experimental work. All authors (E.B., V.V., and M.K.) were involved in the data analysis and interpretation. All the authors participated in the manuscript drafting, reviewing, and editing. All the authors approved the final manuscript. E.B. and M.K. are the principal guarantors of this work; they have full access to all the data and take responsibility for the integrity of the data and accuracy of their analysis.


Competing Interests

The authors declare that they have no competing interests.


Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-08421-z

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


1. Lopes da Silva FH. Neural mechanisms underlying brain waves: from neural membranes to networks. Electroencephalography and clinical neurophysiology. 1991;79:81–93. doi: 10.1016/0013-4694(91)90044-5. [PubMed] [Cross Ref]
2. Shaw, J. C. The brain’s alpha rhythms and the mind. BV Elsevier Science (2003).
3. Niedermeyer E. Alpha rhythms as physiological and abnormal phenomena. International Journal of Psychophysiology. 1997;26:31–49. doi: 10.1016/S0167-8760(97)00754-X. [PubMed] [Cross Ref]
4. Smit CM, Wright MJ, Hansell NK, Geffen GM, Martin NG. Genetic variation of individual alpha frequency (IAF) and alpha power in a large adolescent twin sample. International Journal of Psychophysiology. 2006;61:235–243. doi: 10.1016/j.ijpsycho.2005.10.004. [PubMed] [Cross Ref]
5. Anokhin AP. Genetic psychophysiology: Advances, problems, and future directions. International Journal of Psychophysiology. 2014;93:173–197. doi: 10.1016/j.ijpsycho.2014.04.003. [PMC free article] [PubMed] [Cross Ref]
6. Salmela E, et al. Evidence for genetic regulation of the human parieto-occipital 10-Hz rhythmic activity. European Journal of Neuroscience. 2016;44:1963–1971. [PMC free article] [PubMed]
7. Cantero JL, Atienza M, Salas RM. Human alpha oscillations in wakefulness, drowsiness period, and REM sleep: different electroencephalographic phenomena within the alpha band. Neurophysiologie Clinique/Clinical Neurophysiology. 2002;32:54–71. doi: 10.1016/S0987-7053(01)00289-1. [PubMed] [Cross Ref]
8. Moretti DV, et al. Individual analysis of EEG frequency and band power in mild Alzheimer’s disease. Clinical Neurophysiology. 2004;115:299–308. doi: 10.1016/S1388-2457(03)00345-6. [PubMed] [Cross Ref]
9. Klimesch W. Alpha-band oscillations, attention, and controlled access to stored information. Trends in cognitive sciences. 2012;16:606–617. doi: 10.1016/j.tics.2012.10.007. [PMC free article] [PubMed] [Cross Ref]
10. Reichert JL, Kober SE, Witte M, Neuper C, Wood G. Age-related effects on verbal and visuospatial memory are mediated by theta and alpha II rhythms. International Journal of Psychophysiology. 2016;99:67–78. doi: 10.1016/j.ijpsycho.2015.11.004. [PubMed] [Cross Ref]
11. Sigala R, Haufe S, Roy D, Dinse HR, Ritter P. The role of alpha-rhythm states in perceptual learning: insights from experiments and computational models. Frontiers in computational neuroscience. 2014;8:36. doi: 10.3389/fncom.2014.00036. [PMC free article] [PubMed] [Cross Ref]
12. Hanslmayr S, et al. Alpha phase reset contributes to the generation of ERPs. Cerebral Cortex. 2007;17:1–8. doi: 10.1093/cercor/bhj129. [PubMed] [Cross Ref]
13. Mazaheri, A. & Jensen, O. Rhythmic pulsing: linking ongoing brain activity with evoked responses. Front Hum Neurosci4 (2010). [PMC free article] [PubMed]
14. Klimesch W, Doppelmayr M, Russegger H, Pachinger T, Schwaiger J. Induced alpha band power changes in the human EEG and attention. Neuroscience letters. 1998;244:73–76. doi: 10.1016/S0304-3940(98)00122-0. [PubMed] [Cross Ref]
15. Palva S, Palva JM. Functional roles of alpha-band phase synchronization in local and large-scale cortical networks. Frontiers in psychology. 2011;2:204. doi: 10.3389/fpsyg.2011.00204. [PMC free article] [PubMed] [Cross Ref]
16. Zoefel B, Huster RJ, Herrmann CS. Neurofeedback training of the upper alpha frequency band in EEG improves cognitive performance. NeuroImage. 2011;54:1427–1431. doi: 10.1016/j.neuroimage.2010.08.078. [PubMed] [Cross Ref]
17. Babiloni C, et al. Occipital sources of resting-state alpha rhythms are related to local gray matter density in subjects with amnesic mild cognitive impairment and Alzheimer’s disease. Neurobiology of aging. 2015;36:556–570. doi: 10.1016/j.neurobiolaging.2014.09.011. [PMC free article] [PubMed] [Cross Ref]
18. Chiang A, Rennie C, Robinson P, Van Albada S, Kerr C. Age trends and sex differences of alpha rhythms including split alpha peaks. Clinical Neurophysiology. 2011;122:1505–1517. doi: 10.1016/j.clinph.2011.01.040. [PubMed] [Cross Ref]
19. Nolte G, et al. Identifying true brain interaction from EEG data using the imaginary part of coherency. Clinical Neurophysiology. 2004;115:2292–2307. doi: 10.1016/j.clinph.2004.04.029. [PubMed] [Cross Ref]
20. Kubicki S, Herrmann W, Fichte K, Freund G. Reflections on the topics: EEG frequency bands and regulation of vigilance. Pharmacopsychiatry. 1979;12:237–245. doi: 10.1055/s-0028-1094615. [PubMed] [Cross Ref]
21. Lopes da Silva, F. H. In Niedermeyer’s Electroencephalography: basic principles, clinical applications, and related fields 1203–1226 (Lippincott Williams & Wilkins, 2011).
22. Nasreddine ZS, et al. The Montreal Cognitive Assessment, MoCA: a brief screening tool for mild cognitive impairment. Journal of the American Geriatrics Society. 2005;53:695–699. doi: 10.1111/j.1532-5415.2005.53221.x. [PubMed] [Cross Ref]
23. Pascual-Marqui RD, Michel CM, Lehmann D. Low resolution electromagnetic tomography: a new method for localizing electrical activity in the brain. International Journal of psychophysiology. 1994;18:49–65. doi: 10.1016/0167-8760(84)90014-X. [PubMed] [Cross Ref]
24. Bro R. PARAFAC. Tutorial and applications. Chemometrics and intelligent laboratory systems. 1997;38:149–171. doi: 10.1016/S0169-7439(97)00032-4. [Cross Ref]
25. Carroll JD, Chang J-J. Analysis of individual differences in multidimensional scaling via an N-way generalization of “Eckart-Young” decomposition. Psychometrika. 1970;35:283–319. doi: 10.1007/BF02310791. [Cross Ref]
26. Harshman, R. A. Foundations of the PARAFAC procedure: Models and conditions for an “explanatory” multi-modal factor analysis. (1970).
27. Bro R, Kiers HA. A new efficient method for determining the number of components in PARAFAC models. Journal of chemometrics. 2003;17:274–286. doi: 10.1002/cem.801. [Cross Ref]
28. Möcks J. Decomposing event-related potentials: a new topographic components model. Biological Psychology. 1988;26:199–215. doi: 10.1016/0301-0511(88)90020-8. [PubMed] [Cross Ref]
29. Field AS, Graupe D. Topographic component (parallel factor) analysis of multichannel evoked potentials: practical issues in trilinear spatiotemporal decomposition. Brain topography. 1991;3:407–423. doi: 10.1007/BF01129000. [PubMed] [Cross Ref]
30. Acar E, Aykut-Bingol C, Bingol H, Bro R, Yener B. Multiway analysis of epilepsy tensors. Bioinformatics. 2007;23:i10–i18. doi: 10.1093/bioinformatics/btm210. [PubMed] [Cross Ref]
31. De Vos M, et al. Canonical decomposition of ictal scalp EEG reliably detects the seizure onset zone. NeuroImage. 2007;37:844–854. doi: 10.1016/j.neuroimage.2007.04.041. [PubMed] [Cross Ref]
32. Miwakeichi F, et al. Decomposing EEG data into space–time–frequency components using parallel factor analysis. NeuroImage. 2004;22:1035–1045. doi: 10.1016/j.neuroimage.2004.03.039. [PubMed] [Cross Ref]
33. van der Meij R, van Ede F, Maris E. Rhythmic Components in Extracranial Brain Signals Reveal Multifaceted Task Modulation of Overlapping Neuronal Activity. PloS one. 2016;11:e0154881. doi: 10.1371/journal.pone.0154881. [PMC free article] [PubMed] [Cross Ref]
34. Mengelkamp, J., Weis, M. & Husar, P. In Signal Processing Conference (EUSIPCO), 2013 Proceedings of the 21st European. 1–5 (IEEE).
35. Andersson CA, Bro R. The N-way toolbox for MATLAB. Chemometrics and intelligent laboratory systems. 2000;52:1–4. doi: 10.1016/S0169-7439(00)00071-X. [Cross Ref]
36. Lancaster JL, et al. Automated Talairach atlas labels for functional brain mapping. Human brain mapping. 2000;10:120–131. doi: 10.1002/1097-0193(200007)10:3<120::AID-HBM30>3.0.CO;2-8. [PubMed] [Cross Ref]
37. Maris E, Oostenveld R. Nonparametric statistical testing of EEG-and MEG-data. Journal of neuroscience methods. 2007;164:177–190. doi: 10.1016/j.jneumeth.2007.03.024. [PubMed] [Cross Ref]
38. Ciulla C, Takeda T, Endo H. MEG characterization of spontaneous alpha rhythm in the human brain. Brain topography. 1999;11:211–222. doi: 10.1023/A:1022233828999. [PubMed] [Cross Ref]
39. De Jongh A, et al. The localization of spontaneous brain activity: first results in patients with cerebral tumors. Clinical Neurophysiology. 2001;112:378–385. doi: 10.1016/S1388-2457(00)00526-5. [PubMed] [Cross Ref]
40. Hari R, Salmelin R. Human cortical oscillations: a neuromagnetic view through the skull. Trends in neurosciences. 1997;20:44–49. doi: 10.1016/S0166-2236(96)10065-5. [PubMed] [Cross Ref]
41. Lopes da Silva FH. Functional localization of brain sources using EEG and/or MEG data: volume conductor and source models. Magnetic resonance imaging. 2004;22:1533–1538. doi: 10.1016/j.mri.2004.10.010. [PubMed] [Cross Ref]
42. Liljeström M, Kujala J, Jensen O, Salmelin R. Neuromagnetic localization of rhythmic activity in the human brain: a comparison of three methods. NeuroImage. 2005;25:734–745. doi: 10.1016/j.neuroimage.2004.11.034. [PubMed] [Cross Ref]
43. Salmelin R, Hari R. Characterization of spontaneous MEG rhythms in healthy adults. Electroencephalography and clinical neurophysiology. 1994;91:237–248. doi: 10.1016/0013-4694(94)90187-2. [PubMed] [Cross Ref]
44. Katada A, Ozaki H, Suzuki H, Suhara K. Developmental characteristics of normal and mentally retarded children’s EEGs. Electroencephalography and clinical neurophysiology. 1981;52:192–201. doi: 10.1016/0013-4694(81)90166-8. [PubMed] [Cross Ref]
45. Srinivasan R, Winter WR, Nunez PL. Source analysis of EEG oscillations using high-resolution EEG and MEG. Progress in brain research. 2006;159:29–42. doi: 10.1016/S0079-6123(06)59003-X. [PMC free article] [PubMed] [Cross Ref]
46. Jasper H, Penfield W. Electrocorticograms in man: effect of voluntary movement upon the electrical activity of the precentral gyrus. Archiv für Psychiatrie und Nervenkrankheiten. 1949;183:163–174. doi: 10.1007/BF01062488. [Cross Ref]
47. Perez-Borja C, Chatrian GE, Tyce FA, Rivers MH. Electrographic patterns of the occipital lobe in man: a topographic study based on use of implanted electrodes. Electroencephalography and clinical neurophysiology. 1962;14:171–182. doi: 10.1016/0013-4694(62)90026-3. [PubMed] [Cross Ref]
48. Groppe DM, et al. Dominant frequencies of resting human brain activity as measured by the electrocorticogram. NeuroImage. 2013;79:223–233. doi: 10.1016/j.neuroimage.2013.04.044. [PMC free article] [PubMed] [Cross Ref]
49. Dodrill CB, Wilkus RJ, Tassinari C, Scheffner D. Neuropsychological correlates of the electroencephalogram in epileptics: II. The waking posterior rhythm and its interaction with epileptiform activity. Epilepsia. 1976;17:101–109. doi: 10.1111/j.1528-1157.1976.tb03388.x. [PubMed] [Cross Ref]
50. Stoller A. Slowing of the alpha-rhythm of the electro-encephalogram and its association with mental deterioration and epilepsy. The British Journal of Psychiatry. 1949;95:972–984. doi: 10.1192/bjp.95.401.972. [PubMed] [Cross Ref]
51. Bollimunta A, Chen Y, Schroeder CE, Ding M. Neuronal mechanisms of cortical alpha oscillations in awake-behaving macaques. The Journal of Neuroscience. 2008;28:9976–9988. doi: 10.1523/JNEUROSCI.2699-08.2008. [PMC free article] [PubMed] [Cross Ref]
52. Mo J, Schroeder CE, Ding M. Attentional modulation of alpha oscillations in macaque inferotemporal cortex. The Journal of Neuroscience. 2011;31:878–882. doi: 10.1523/JNEUROSCI.5295-10.2011. [PubMed] [Cross Ref]
53. Orban GA, Van Essen D, Vanduffel W. Comparative mapping of higher visual areas in monkeys and humans. Trends in cognitive sciences. 2004;8:315–324. doi: 10.1016/j.tics.2004.05.009. [PubMed] [Cross Ref]
54. Mitzdorf, U. Current source-density method and application in cat cerebral cortex: investigation of evoked potentials and EEG phenomena. (American Physiological Society, 1985). [PubMed]
55. Lopes da Silva, F. H. & Rotterdam, A. In Electroencephalography: basic principles, clinical applications, and related fields 93–109 (Lippincott Williams & Wilkins, 2005).
56. Manshanden I, De Munck JC, Simon NR, Lopes da Silva FH. Source localization of MEG sleep spindles and the relation to sources of alpha band rhythms. Clinical Neurophysiology. 2002;113:1937–1947. doi: 10.1016/S1388-2457(02)00304-8. [PubMed] [Cross Ref]
57. Makeig S, et al. Electroencephalographic brain dynamics following manually responded visual targets. PLoS Biol. 2004;2:e176. doi: 10.1371/journal.pbio.0020176. [PMC free article] [PubMed] [Cross Ref]
58. Van Beijsterveldt C, Van Baal G. Twin and family studies of the human electroencephalogram: a review and a meta-analysis. Biological Psychology. 2002;61:111–138. doi: 10.1016/S0301-0511(02)00055-8. [PubMed] [Cross Ref]
59. Farber D, Vildavsky VY. Heterogeneity and age dynamics of the alpha-rhythm of the electroencephalogram. Human Physiology. 1996;22:517–524.
60. Schulman JJ, et al. Imaging of thalamocortical dysrhythmia in neuropsychiatry. Frontiers in human neuroscience. 2011;5:69. doi: 10.3389/fnhum.2011.00069. [PMC free article] [PubMed] [Cross Ref]
61. Knyazeva MG, et al. Evolution of source EEG synchronization in early Alzheimer’s disease. Neurobiology of aging. 2013;34:694–705. doi: 10.1016/j.neurobiolaging.2012.07.012. [PubMed] [Cross Ref]
62. Barzegaran E, van Damme B, Meuli R, Knyazeva MG. Perception-related EEG is more sensitive to Alzheimer’s disease effects than resting EEG. Neurobiology of aging. 2016;43:129–139. doi: 10.1016/j.neurobiolaging.2016.03.032. [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group