Prestimulus amplitudes modulate P1 latencies and evoked traveling alpha waves

Abstract

Traveling waves have been well documented in the ongoing, and more recently also in the evoked EEG. In the present study we investigate what kind of physiological process might be responsible for inducing an evoked traveling wave. We used a semantic judgment task which already proved useful to study evoked traveling alpha waves that coincide with the appearance of the P1 component. We found that the P1 latency of the leading electrode is significantly correlated with prestimulus amplitude size and that this event is associated with a transient change in alpha frequency. We assume that cortical background excitability, as reflected by an increase in prestimulus amplitude, is responsible for the observed change in alpha frequency and the initiation of an evoked traveling trajectory.

Introduction

Brain oscillations as measured by the EEG are manifested by rhythmic fluctuations in amplitude size that show a relation to membrane currents of masses of neurons (for a review, cf. Buzsáki et al., 2012). They reflect phases of low and high excitability which are associated with phases of decreased and increased firing rate. This basic timing mechanism of oscillations has been well described for a variety of different rhythms (Buzsáki, 2006; Steriade, 2006). The general hypothesis is that temporal aspects of neuronal firing, which are crucial for information transmission, are largely and causally organized by oscillations (Buzsáki and Draguhn, 2004; Romei et al., 2010). Because cortical neurons form large networks (Hagmann et al., 2008), information transfer between neurons implies a spatiotemporal distribution of oscillations which is manifested by traveling waves (Wu et al., 2008).

Traveling waves have been observed very early in electrophysiological research (Adrian and Matthews, 1934; Petsche and Marko, 1955; Hughes, 1995) and have been documented on all levels of measurements, including the synaptic level (Williams, 1992; Wu et al., 2008; Nadasdy, 2009), the local field potential recorded from microelectrodes in animals and humans (Rubino et al., 2006; Takahashi et al., 2011; Nauhaus et al., 2012), the ECoG, EEG, and MEG (Freeman and Barrie, 2000; Alexander et al., 2013) in a large variety of studies focussing, e.g., on sensory encoding processes (Ermentrout and Kleinfeld, 2001; Dong and Olson, 2008), information transfer between cortical regions (Rubino et al., 2006) and methodological aspects (Alexander et al., 2006; Hallatschek, 2010). They have been observed in the ongoing EEG (Bahramisharif et al., 2013), in evoked potentials (Alexander et al., 2006, 2013; Klimesch et al., 2007a), in the flicker induced steady state evoked potential (Silberstein et al., 1990; Burkitt et al., 2000), and in the sleep EEG (Massimini et al., 2004). The physiological process underlying traveling waves represents a crucial feature for theories on encoding (Nadasdy, 2010) and for brain theories describing the functional interplay between different neuronal networks (Nunez, 2000; Nunez and Srinivasan, 2014). The most general finding is that traveling waves are ubiquitous in the ongoing EEG, but occur intermittently in episodes (Patten et al., 2012). A special and important phenomenon are evoked traveling waves (Alexander et al., 2006, 2009; Klimesch et al., 2007a) because they link the ongoing EEG with properties of the event-related potential (ERP). In a study, focusing on the within trial spatial phase distribution, Alexander et al. (2013) have demonstrated that the trial averages of traveling waves reproduced the topography of ERPs. Particularly interesting is the observation that topographical latency differences of the P1 component can be described in terms of an evoked traveling alpha wave with a posterior lateral to medial traveling movement in lexical and semantic decision tasks (Klimesch et al., 2007a; Fellinger et al., 2012; Zauner et al., 2014). This means that the phase topography of evoked alpha, but not that of other EEG frequencies, (Gruber et al., 2005; Klimesch et al., 2007a) coincides with the topography of the P1 components. Since the pioneering work of Berger (1929), alpha waves are known as an important phenomenon of the ongoing EEG. An interesting question, arising from these observations is which processes underlie the transition from an ongoing to an evoked traveling wave. This is the question, we ask in the present study.

A simple method for detecting an evoked traveling wave is based on two requirements. First, an evoked component must exhibit a spatiotemporal distribution that allows for determining a leading and trailing electrode site. Second, the evoked component must have a dominant frequency characteristic within a typical frequency range (such as e.g., alpha). The P1 component fulfills these requirements. Several studies have shown that the P1 has a typical spatiotemporal distribution and a characteristic frequency in the alpha frequency range (Gruber et al., 2005; Klimesch et al., 2007b; Fellinger et al., 2012). The leading site is that site with the shortest P1 latency, whereas the trailing site is that with the longest P1 latency. An evoked traveling wave has – by definition – a clear onset of traveling movement peri- or poststimulus.

What kind of physiological process would be capable of inducing a topographical specific frequency change at the leading site? To our knowledge, this question has not yet been addressed in an empirical study. In an attempt to answer this question, we proceed from the idea that a traveling trajectory can be induced by a transient frequency change. When we have exactly the same oscillatory frequency at two recording sites, the spatiotemporal phase relationship will remain constant. If, however, the oscillatory frequency exhibits a transient change at one of the two sites, a clear traveling trajectory may emerge as illustrated in Figure 1A. This transient change in frequency, which is manifested by a change in the period of one or more oscillatory cycles, can be understood as a special type of phase ‘reset’ (cf. Cases 2a,b in Figure 1A). We assume that this frequency change occurs at the leading site, because it is that site from where the traveling wave spreads to other sites.

FIGURE 1

What could be the physiological mechanism underlying a transient frequency change? We proceed from the idea that a change in ‘background’ excitation and/or inhibition might be related to a change in oscillatory frequency (Nunez and Srinivasan, 2014). For slow EEG fluctuations (with a frequency of about 1 Hz and slower) there is good evidence that they reflect cyclical variations in the excitability of neuronal ensembles (e.g., Bishop, 1933; Steriade et al., 1993; Contreras et al., 1996; Sanchez-Vives and McCormick, 2000). This research has shown that action potentials are generated during an excitatory up state but not during an inhibitory down state (Amzica and Steriade, 1997). The scalp surface polarity of such an up and down state is difficult to predict, because it depends (beside other factors) on the exact location of the source. Nonetheless, based on these findings one may speculate, that the magnitude of a slow deflection (in the sub-delta or delta frequency range) reflects a change in background excitation and is correlated with P1 latency in a sense, that a large deflection is related to a large frequency change and a large change in P1 latency. There is yet another, additional, possibility that must be considered. Research on rat hippocampus gamma oscillations has shown that instantaneous oscillatory frequency (as measured in terms of the duration of the gamma period in single trials) depends on the extent of excitation and inhibition (Whittington et al., 1995; Traub et al., 1996; Atallah and Scanziani, 2009). In an interesting study by Atallah and Scanziani (2009) clear evidence was found that fluctuations in (instantaneous) gamma amplitude reflect changes in synaptic excitation and are associated with fluctuations in (instantaneous) gamma period. The basic finding was that an increase in amplitude is closely associated with a lengthening in the immediately following period, and – vice versa – a decrease in amplitude is associated with a shortening in the immediately following period. We refer to this finding as ‘cycle to cycle fluctuations in amplitude and period.’ It is manifested by a significant positive correlation between amplitude and period on a cycle per cycle basis.

These two mechanisms, slow waves (associated with up and down states) and cycle to cycle fluctuations in oscillatory amplitude size and period length must not be considered mutually exclusive interpretations. It might well be possible that slow wave components are associated with increased oscillatory cycle to cycle fluctuations in different frequency ranges.

The measurement of cycle to cycle fluctuations of a frequency of interest – in our case the alpha band - requires a specific analyzing procedure. First, for single trials, the time points of peaks and troughs are identified for the band pass filtered data. Then, and most importantly, these time points are used to determine the amplitudes of peaks and troughs in the raw data. The reasons for this procedure is, that alpha band pass filtering abolishes the influence of (and a possible interaction with) slow components and in addition tends to reduce asymmetric and short lasting amplitude fluctuations. This method, which determines the amplitudes of peaks and troughs in the raw data, was suggested by researchers who found evidence for asymmetric alpha amplitude fluctuations (e.g., Nikulin et al., 2007; Mazaheri and Jensen, 2008) will be applied here in addition to traditional ERP analyses.

In the present study, we used a semantic (living vs. non-living) judgment task that already proved useful to investigate evoked alpha waves (Zauner et al., 2014). One of the basic findings was that words with many semantic features tend to reduce the latency difference between leading and trailing sites but for yes response trials (to living objects) only. Here, we want to test the above described hypotheses regarding the predicted association between the amplitude of a slow prestimulus ERP deflection on the one hand and a change in P1 latency at the leading site on the other hand. In addition, we want to test, whether a transient frequency change in alpha can be observed at the subject level (i.e., for the averaged EEG data as reflected by the ERP) but also on the single trial level (i.e., for the not averaged data). This latter method aims to detect the single trial P1 and to analyze alpha amplitude to period fluctuations for several cycles immediately preceding the P1.

Materials and Methods

Subjects

A sample of 34 subjects (all different from those, used in Zauner et al., 2014), students of the University of Salzburg, participated in the present study after giving informed consent. The sample consisted of 14 male (mean age = 25.1 ± 2.5 years) and 20 female subjects (mean age = 23.2 ± 3.4 years). These subjects had normal or corrected-to-normal vision, did not report neurological diseases and were not on psychotropic medication. All subjects were compensated by course credits. The experiment was conducted according to the code of ethics (World Medical Association, 1996) and was approved by the Ethics Committee of the University of Salzburg.

Stimulus Material and Task

The stimulus material and task was identical with a previously published study by Zauner et al. (2014). Here we give just a brief outline of the most important features of the stimulus material and task design. Subjects performed a semantic (living vs. non-living) judgment task. A total of 280 German words were presented in a randomized sequence. For each word, subjects had to decide, whether it represents a living or non-living object. The living and non-living category consisted of 140 words each, which were subdivided into 70 words with a large number of visual-semantic features (NOF+) and 70 words with a small NOF-. An example for a NOF+ word is eagle, an example for a NOF- word is hamster. The words were taken from McRae et al.’s (2005) “semantic feature production norms” and were translated into German.

Subjects were seated in a comfortable chair in front of a computer monitor (75 Hz refresh rate) at a distance of about 130 cm. The fixation cross was replaced by a word written in upper-case-letters either belonging to the living or non-living category. The word was presented for 1000 ms in dark gray (horizontal angle: 2,8∘–4,3∘; vertical angle = 0,66∘) within a bright gray box (9,7∘ × 2,6∘) to ensure comfortable reading and to hold visual surface features constant between trials. The interval between the onset of the fixation cross and the onset of the word varied between 400 and 600 ms in 50 ms intervals in order to reduce onset expectations. The intertrial interval varied between 1900 and 2100 ms. Subjects indicated by button press on the keyboard with their left index finger when the word denoted a living (yes-response) and the right index finger when the word denoted a non-living object (no-response).

EEG Recordings

We used a 64-channel BrainAmp amplifier (BrainProducts, Inc., Gilching, Germany) for EEG recording. Signals were online referenced against the nose and subsequently (off-line) re-referenced to digitally averaged [(A1 + A2)/2] ear lobes. Recording bandwidth was set from 0.016 to 100 Hz and a notch filter at 50 Hz. Signals were digitized at a sampling rate of 1000 Hz 60 Ag–AgCl-electrodes were mounted using an EasyCap. For the present study we used only 17 posterior electrodes (P7, P5, P3, P1, Pz, P2, P4, P6, P8, PO7, PO3, POz, PO4, PO8, O1, Oz, O2) for data analysis. Impedances were kept below 8 kΩ. To control for vertical and horizontal eye movements two bipolar EOG-channels were mounted. After re-referencing, epochs containing eye artifacts were corrected by the Gratton et al. (1983) and muscle artifacts were rejected. BrainVisionAnalyzer (BrainProducts, Inc.) was used for data analyses. Epochs consisted of EEG segments ranging from -600 to 1000 ms relative to the stimulus.

Data Analysis: ERPs, Prestimulus Power, and Transient Changes in Evoked Alpha Frequency

For the calculation of ERPs we used two different filters, a broad band filter between 0.5 and 70 Hz. For alpha bandpass filtering we used an equiripple filter (7.5–12 Hz passband, 24 db at cutoff frequencies of 6 and 14 Hz). Prestimulus power was determined as the average of rectified voltage within a time window of – 100 to 200 ms relative to stimulus onset (0 ms) for the broadly filtered data (i.e., the raw data), a subdelta band (0.5–2.5 Hz), and the alpha band.

The alpha filtered data were used to calculate P1 peak latencies and amplitude. This was done because the P1 peak is easier to detect in the filtered data. In the ERPs of the broadly filtered data, evoked beta activity can be observed that leads to double peaks in the P1 time window and makes the detection of the P1 ambiguous. For the calculation of evoked frequency changes we used a narrow alpha band pass filter (with a width of 2 Hz) that was adjusted to individual alpha peak frequency (IAF) as center frequency (IAF ± 1 Hz). IAF was determined as the spectral component within an extended alpha frequency range of 7–14 Hz which showed the largest power during a pre-task eyes closed resting period of 5 min. The obtained IAF’s were averaged over four electrodes (P3, P4, Po3, and Po4). The segmented data were averaged and peaks were detected semi-automatically.

The typical latency of the visual P1 is between 80 and 120 ms post-stimulus (cf. Figure 1B for an illustration of the P1 latency differences). The semi-automatic algorithm detects positive peaks within a defined time window. We used a window of 70–185 ms because we expected a large variation in P1 latency due to a traveling movement. The algorithm defines any positive going deflection as positive peak even if the absolute amplitude has a negative voltage. By visual inspection we accepted only those peaks which had a positive amplitude larger than 0.5 μV. This procedure was done separately for each of the four conditions (living NOF+, NOF- and non-living NOF+, NOF-). In cases of missing values we used the respective values of the remaining condition(s). For each subject, response type (yes and no response to words denoting a living or non-living object) and semantic feature condition (NOF+ and NOF-), we determined that electrode with the shortest P1 latency and that with the longest P1 latency. These electrodes are referred to as leading and trailing sites in the following. As depicted in Figure 1C, PO8 was the most frequent leading site and Pz the most frequent trailing site. We then calculated the latency differences between the leading and trailing sites as a gross measure of the traveling movement of the evoked alpha wave for each subject, response type and feature condition.

We also used a peak detection procedure to determine transient changes in alpha period (as measure of ‘instantaneous’ alpha frequency) around stimulus onset. For this purpose we used the averaged alpha filtered data for the leading electrode and determined the interpeak latencies between six consecutive positive peaks. As starting point, we used the second positive alpha peak that occurred poststimulus. From this peak we calculated the interpeak latencies between six consecutive positive peaks (i.e., comprising the second positive peak poststimulus plus five consecutive peaks) going backward in time into the prestimulus interval. This procedure yielded five different period estimates for each subject, averaged over all conditions. In the following, the first period value is referred to as poststimulus period. The next following value was termed peristimulus period. The last three values were averaged and termed prestimulus period. This was done to reduce differential influences of the filter stemming from evoked components.

Single Trial Data Analysis: Prestimulus Amplitude and Transient Changes in Alpha Frequency

To test, whether fluctuations in prestimulus amplitude in single trials are associated with fluctuations in P1 latency, we used a recently suggested method for the detection of single trial P1 components (Gruber et al., 2014). This method is based on the hypothesis that ongoing alpha develops more or less seamlessly (but with some jitter) into the P1 component. It selects those trials, where a positive alpha peak can be observed within the time window of the P1 component of the ERP. The selection algorithm comprises four steps. First, we define a selection window as an interval of ±32 ms around the individual ERP P1 latency at the leading electrode PO8. This window – referred to as window 1 in the following – was chosen to be long enough to cover a half cycle of slow alpha of about 8 Hz. Second, within this window, we search for a positive voltage alpha peak (the single trial P1) in the alpha filtered data. The search criterion is a positive peak with positive amplitude in the individually determined 64 ms time window. Third, those trials for which the peak detection yielded a positive result constitute the subset of selected trials. It comprises those trials, in which ongoing alpha develops more or less seamlessly (within the defined selection window) into the P1. For these trials – which are simply termed selected trials in the following – the single trial P1 latencies (the latency of the selected positive alpha peak) were determined. The average number of selected trials per subject was 172 out of 206. This means that 83% of all trials were alpha phase aligned to the individual P1 component. This subset of trials constitutes the selected trials. The remaining 17% of trials are those that are not phase aligned to the P1. They are termed rejected trials. For the investigation of fluctuations in alpha frequency in single trials, we used the selected trials, obtained in step three as described above. In these data, we used the single trial P1 as starting point for determining the interpeak latencies between five consecutive positive peaks going backward in time into the prestimulus interval as illustrated in Figure 2A. This procedure allowed us to measure the length of five consecutive alpha periods by going backward in time from the single trial P1 to the prestimulus period. These are termed p+(1) … p+(5) for the periods between positive alpha peaks and p-(1) … p-(5) for the respective periods between negative alpha peaks starting from the first negative amplitude preceding the single trial P1 (cf. Figure 2A).

FIGURE 2

To study fluctuations in instantaneous EEG, amplitudes that are phase locked to the positive alpha peaks of each of the five cycles (preceding the single trial P1), we determined for each trial the time points of these peaks. At each of these five time points we determined the amplitude of the raw EEG data and the alpha filtered data. For the raw (r) data, the five amplitude estimates are termed ar+(1) … ar+(5), denoting the first … fifth raw amplitude (coinciding with the positive alpha peak) preceding the single trial P1. For the alpha filtered data (a), the five amplitude estimates are termed aa+(1) … aa+(5), denoting the first … fifth positive alpha amplitude preceding the single trial P1. The same procedure was applied for the amplitudes at the negative alpha peaks and the respective values are termed ar-(1) … ar-(5) and aa-(1) … aa-(5), respectively.

Single Trial Data Analysis: Control Analysis

The above described procedure aims at identifying those trials that exhibit a more or less seamless development of alpha phase with respect to the P1 component of the ERP. As a consequence, trials not meeting this criterion (which depends on the width of the selection window) are rejected. In order to test, whether the expected cycle to cycle fluctuations in amplitude and period depend on the width of the selection window, we ran a control analysis with an extended selection window (of ±50 ms), which was chosen to be long enough to cover more than a half cycle of slow alpha of about 8 Hz. This selection window of the control analysis is referred to as window 2 in the following. Because with this broad window more than one positive peak can appear, we selected that positive peak as single trial P1 which is closest to the individual P1 of the ERP. In contrast to the selection procedure with time window 1, which yielded 17% rejected trials, the percentage of rejected trials with time window 2 was only 2%.

Statistical Analysis

We first wanted to determine, whether the present study replicates the findings from Zauner et al. (2014). For this purpose, we used a two way ANOVA with the factors response type (yes and no responses) and semantic condition (NOF+/-). The dependent measure was latency difference between leading and trailing sites. Across subjects, correlations were calculated between prestimulus power and mean P1 latency (averaged over all conditions) for the most frequent leading electrode PO8. To test for transient changes in instantaneous alpha frequency, we calculated a priori t-tests to test the hypothesis that the predicted transient change in alpha period occurs selectively for the peristimulus period. In addition, we performed a two way repeated measures ANOVA with electrode (Pz, PO8) and time (pre-, peri-, and poststimulus period) as factors. The three levels of factor time represent (i) the first poststimulus alpha period (measured from the second positive evoked peak that appears after stimulus onset to the next positive peak, by going backward in time), (ii) the next period (backward in time) constitutes the peristimulus period. (iii) The prestimulus period is the average of the three consecutive periods that follow the peristimulus period (again determined by going backward in time). The dependent measure was the duration of the alpha period determined for IAF.

For the single trial analysis, we calculated correlations between the single trial P1 latencies, the five preceding amplitude estimates at the positive and negative alpha peaks, and between the amplitude estimates and the immediately preceding periods. Because our hypothesis predicts a positive association between amplitude size and period, we used one-sided significance values at the 5% level. We also calculated one way ANOVA’s with factor TIME (with five levels, denoting the first … fifth cycle preceding the P1) and amplitude/period as dependent measures.

Results

Behavioral Data and Replication Analysis

The overall mean RT was 688.94 ms and the mean percentage of correct responses was 96%. The RT’s in the four conditions, yes responses NOF+/-, and no responses NOF+/-, were in that order: 648.9 (SD = 58.4), 677.5 (SD = 63.9), 704.3 (SD = 79.0), 725.1 (SD = 86.4). IAF varied between 8.1 and 12.9 Hz around a mean of 10.3 Hz (SD = 1.01).

The two way ANOVA with semantic condition (NOF+/-) and response type (yes vs. no response) as factors and latency difference (between leading and trailing sites) as dependent measure yielded a significant interaction [F_(1,33) = 4.987, p < 0.05). The main effects did not reach significance. Inspection of the respective means indicates that a large number of features (NOF+) tends to increase latency differences for no responses, but tends to decrease latency difference for yes responses. These findings replicate the basic aspects of the Zauner et al. (2014) study.

Prestimulus Power, P1 Latency, and Peristimulus Changes in Alpha Period

Visual inspection of the ERPs at the most frequent leading and trailing sites reveals large differences in prestimulus power (for the raw data) and P1 latency as depicted in Figure 1B. The average P1 latencies are 116.0 ms (SD = 16.3 ms) for PO8 and 135.8 ms (SD = 10.4 ms) for Pz. The calculation of correlations between the respective variables showed a highly significant association for the leading electrode (PO8; r = -0.56, p < 0.01) but not for the trailing electrode (Pz; r = -0.12). As shown in Figure 1D, at PO8, large prestimulus power is significantly associated with short P1 latencies. The correlations between PO8 P1 latency and the subdelta and alpha filtered prestimulus power yielded a significant result for the subdelta power only (r = -0.53; testing for symmetrical distributions of both variables, using the Kolmogorov–Lillieforce procedure, yielded no significant deviation). This finding suggests that the critical component for the correlation is a slow wave in the subdelta frequency range. Furthermore, dependent t-tests revealed a larger P1-amplitude for Po8 electrode compared to Pz electrode [t₍₃₃₎ = 3.198, p = 0.003).

Changes in alpha period are depicted in Figure 1E. A priori t-tests showed that only the peristimulus alpha period was longer for Po8 as compared to Pz [t₍₃₃₎ = 2.274, p = 0.030]. The pre- and poststimulus alpha- period length did not differ significantly. The two-way repeated measures ANOVA yielded a significant main effects for electrode site [leading vs. trailing, F_(1,33) = 9.178, p = 0.005] and time [F_(2,66) = 3.864, p = 0.027]. The interaction, however, did not reach significance. The significant main effect suggests slightly different frequencies at the two electrodes. However, the calculation of a t-test between IAF at Pz and PO8 (as measured during the resting period) did not show significant differences.

Single Trial Analysis of Peri-/Prestimulus Amplitude

The single trial analysis focused on two different, but closely interrelated questions. One question refers to the association between single trial P1 latency and the magnitude of peri-/prestimulus amplitude, another to single trial fluctuations in amplitude and period.

Peri-/Prestimulus Amplitude and Single Trial P1 Latency

The correlations between the single trial P1 latencies and the amplitude estimates of the raw EEG [ar+(1) … ar+(5); ar-(1) … ar-(5)] at the time points of the five preceding positive and negative alpha peaks are summarized in Table 1A. The results for positive peaks show that 22 of the 34 subjects exhibit significant negative correlations between latency and ar+(1). These correlations indicate that an increase in ar+ is associated with a decrease in the single trial P1 latency. The respective correlations with the amplitudes at the negative alpha peaks are somewhat weaker, but exhibit an analogous relationship. The larger the magnitude of the negative amplitude [ar-(1)], the shorter is P1 latency.

Single Trial Fluctuations in Amplitude and Period

When we focus on the cycle per cycle correlations between amplitude estimates [(ra+(1) … ra+(5) and ra-(1) … ra-(5)] and the respective, corresponding single trial periods [p+(1) … p+(5) and p-(1) … p-(5)] we see primarily positive correlations for positive amplitudes and negative correlations for negative amplitudes. These findings provide evidence for single trial fluctuations in amplitude and period and show that an increase in the magnitude of the raw amplitude tends to increase the period for the immediately following cycle.

It must be noted that the positive association between amplitude and period as depicted in Table 1B does not contradict the negative association between amplitude and single trial P1 latency as shown in Table 1A. The reason lies in the different ways P1 latency and period are measured. Latency is measured in relation to stimulus onset, but period is measured from the single trial P1 to the preceding amplitudes. Thus, as illustrated in Figure 2B, a long period preceded by a large amplitude is associated with a comparatively short latency, whereas a short period preceded by a small amplitude is associated with a long P1 latency.

The correlations for the alpha filtered single trial amplitudes exhibit a very similar pattern of results. However, the number of significant correlations is much smaller than for the raw data.

Single Trial Analysis: ANOVA Results

The results of the two ANOVA’s with TIME as factor and the subject averages for p+ and ar+ as dependent measures yielded significant effects [F_(2.10,69.44) = 12.63, p < 0.001 and F_(2.39,78.77) = 35.48, p < 0.001; the Greenhouse Geisser corrected df-values are depicted]. The findings show two different aspects. First, fluctuations in amplitude and period are highly significant between cycles. Secondly, the fluctuations in both variables, as depicted in Figure 3, co-vary in the same way, as was statistically documented by the correlational analysis as summarized in Table 1B. Amplitude size and the length of the immediately following period are perfectly associated: a large amplitude is followed by a long period and a small amplitude is followed by a short period. Most interestingly, the ANOVA with the subject averages for p- did not show significant differences. This suggests that the association between amplitude size and period is stronger for positive cycles.

FIGURE 3

Control Analysis

The control analysis showed very similar results, with a tendency to exhibit somewhat more significant correlations as a comparison of the findings for window 1 and window 2 in Tables 1A,1B reveals. Also, the respective ANOVA’s with TIME as factor and the subject averages for p+ and ar+ as dependent measures yielded very similar effects [F_(2.04,67.33) = 8.96, p < 0.001 and F_(2.29,75.67) = 33.84, p < 0.001; the Greenhouse Geisser corrected df-values are depicted]. Again, the ANOVA with the subject averages for p– did not show significant differences.

Discussion

The most important findings are the significant correlations between prestimulus power and P1-latencies and the peristimulus increase in the alpha period at the leading electrode PO8 (cf. Figures 1B,D,E). In the ERP analysis, the latter effect is weak, because only a priori t-tests showed significant results. In the single trial analysis, however, highly significant differences in the magnitude of amplitudes and the length of alpha periods preceding the P1 were observed at the leading electrode as depicted in Figure 3. Most importantly, the fluctuations in alpha period are associated with amplitude (in the raw EEG, measured at the time point of alpha peaks) in a way that, a large amplitude is associated with a prolonged period in the immediately following cycle (cf. Figure 2B for a paradigmatic illustration of this finding). As an example, let us consider the first positive amplitude preceding the P1 and the first period (i.e., interpeak latency between the P1 and that first positive amplitude) as depicted in Figure 3. These two values are reflected by the first data points of the red and blue lines at time point 1 (i.e., at the time of the first positive alpha peak in the first cycle). As shown in Figure 3, amplitude size and period length co-vary over time (over cycles). It should be noted that the association between amplitude and alpha period is much weaker, if amplitudes are measured at the time points of the negative alpha peaks and interpeak latencies are calculated between the negative peaks (cf. Figure 2A for an illustration for the measurement of positive and negative cycles). We consider this an important finding, suggesting that amplitudes at positive and negative alpha peaks may be modulated differentially (cf. Nikulin et al., 2007; Mazaheri and Jensen, 2008) and that the modulation of the positive peaks is the critical factor that is associated with changes in P1 latency. We also have to emphasize that the results obtained for alpha amplitude size (instead of the amplitude of the raw EEG measured at the alpha peaks) are similar but much weaker as compared to the raw data. This latter finding most likely is due to the influence of the alpha band pass filter which abolishes slow drifts and tends to reduce the magnitude of short lasting, transient fluctuations in positive or negative peak amplitudes.

Taken together, these findings suggest that, in general, an increase in prestimulus amplitudes (of the raw EEG) is associated with a shortening of the P1 latency. The single trial analysis shows that amplitudes decrease from cycle 5 to cycle 2 (cf. Figure 3) in a similar way as power decreases (from about – 400 ms prestimulus) in the grand average ERP (cf. Figure 3). Because cycles (with a period around 100 ms) are counted relative to the P1, cycle 5 represents – in terms of time prestimulus – a time window around 400 ms. Thus, the time course of prestimulus amplitudes as measured at the positive alpha peaks in the single trial analysis and the time course of the slow component in the ERP match each other. In addition, the single trial analysis has revealed that an increase in the amplitude in cycle 1 (immediately preceding the P1) is associated with a lengthening in the alpha period (in that cycle) and a decrease in P1 latency. This latter finding resembles closely Case 2b (illustrated in Figure 1A) which describes a situation, where a peristimulus lengthening of the alpha period leads to a shortening in the P1 latency. Thus, as illustrated by Case 2a in Figure 1, a transient frequency change – elicited by a transient increase in the alpha period – may be the key factor for triggering a traveling wave with a well defined trajectory and that a topographical specific amplitude increase triggers that change.

It is also important to note that the results of the single trial analysis contradict the evoked model of ERP generation but also the phase reset model (at least in its most radical formulation). The reason is that 83% of all trials were found to be already aligned in phase in a way that voltage positive alpha peaks develop more or less seamlessly into the P1 (with a maximal jitter of ±32 ms). These findings replicate those of a visual target detection study (Gruber et al., 2014), which in addition has shown that phase aligned trials (i.e., selected trials) are associated with shorter detection times. We cannot rule out that the appearance of the fixation point, preceding the stimulus word, may play a role for the observed phase alignment. Because the interval was jittered (randomly between 400 and 600 ms) a direct influence of the fixation point seems unlikely. Seamless alpha is one argument against the evoked and phase reset model. An additional argument is that the rejected trials exhibit a single trial P1 that is counter phase to the positive peak in the P1 time window. This fact is in part due to the width of the selection window (±32 ms) which is close to a half period of slow alpha (with a frequency of about 8 Hz). If the P1 would be generated by fixed latency fixed polarity components, opposite polarities in single trials should not emerge within the time window of the P1. A reformulation of the phase reset model in the sense that phase is not abruptly reset but ‘reorganized’ in most trials (for a review cf. Klimesch et al., 2007c) would be very well in line with the results of the present study and those reported by Gruber et al. (2014).

Although the evoked model is not in line with the results of the single trial analysis of the leading site, the evoked model is capable of explaining latency differences between PO8 and Pz which we have interpreted in terms of a traveling wave. The theoretical basis for the evoked traveling wave model and that of the traditional evoked model are radically different. The former assumes that an event related phase reorganization at the leading site triggers an evoked traveling wave of different evoked components that move to the trailing site. In contrast, the traditional ERP model- based on the fixed polarity, fixed latency concept – assumes that different evoked components are generated at different sites within typical latency windows. As an example, according to the evoked traveling wave model, the C1 component (the negative peaks preceding the P1 at PO8 and Pz; cf. Figure 1) would be considered the negative peak of an evoked alpha wave that develops into the P1, and topographical latency differences are the result of the traveling movement. According to the evoked model topographical latency differences of one component may be due to a superposition with another component. With respect to our example it may be argued that the long latency of the P1 at Pz may be due to the influence of a pronounced C1 which delays the appearance of the P1. Because the C1 is usually larger at midline as compared to lateral sites (e.g., Clark et al., 1995), the delaying influence of the C1 on the P1 would be larger at midline as compared to lateral topographies. On the basis of our findings, we cannot rule out this alternative interpretation of the observed P1 latency differences between PO8 and Pz, although in an earlier study, we could show that the traveling movement also comprises the negative polarity peak preceding the P1 (Klimesch et al., 2007a). An empirical evaluation of these two conflicting interpretations would be possible in an experiment that aims to vary the C1 component according to the cruciform model which considers the different polarities stemming from activations of the lower versus upper banks of the calcarine fissure (e.g., Clark et al., 1995; Di Russo et al., 2002). For lower field presentations with a positive polarity C1, the evoked traveling model would predict a negative component within a half cycle of alpha. The evoked model would predict a superposition with a positive component within the typical latency window of the P1. A single trial analysis of the C1 component – in an analogous way as we have performed for the P1 – could clarify this question if lower and upper field presentations are compared.

Our findings are in perfect agreement with predictions of the global wave model, suggested by Nunez and Srinivasan (2014). One of the central assumptions is that the modulation density of action potentials is a function of cortical background excitability and inhibitory feedback strength. Quantification of the model predicts that an increase in parameter β (reflecting the degree of cortical background excitability) is associated with an increase in oscillatory amplitude but a decrease in frequency. This is exactly the pattern of results, we found for the peri- and poststimulus interval at the leading electrode. Our findings are to our knowledge the first for the human scalp EEG, showing cycle per cycle fluctuations of amplitude and period in alpha, resembling closely findings on rat hippocampal gamma (Whittington et al., 1995; Traub et al., 1996; Atallah and Scanziani, 2009).

For the physiological interpretation of traveling waves, different conceptualisations have been used. As an example, spatial phase shifts may originate from a phase lag between neighboring, coupled neuronal oscillators (Ermentrout and Kleinfeld, 2001; Wu et al., 2008). Here, however, we refer to the water wave analogy, Nunez (2000) uses to illustrate the dynamics of cortical traveling waves, spreading via myelinated fibers that connect neighboring cortical regions. This model has the advantage that it has a strong emphasis on large scale white matter connectivity, which presumably plays also an important role for cognitive processes and the scalp EEG as well (Hipp et al., 2012). It should be noted, however, that this model, which assumes fast intercortical propagation has been challenged recently. Hindriks et al. (2014) point out that local field potential recordings in animal studies show slow propagation velocities that suggest intra- instead of intercortical propagation. Using a combination of EEG analysis and biophysical modeling, the authors demonstrate that the fast scalp velocities can be accounted for by slow traveling oscillations. The central idea is that a slow traveling movement in a small cortical region may ‘project’ via volume conduction to a large scalp region, thereby mimicking a fast traveling movement. For our results, this would mean that the evoked traveling alpha movement would be restricted to a comparably small cortical region.

With respect to the functional meaning of the evoked traveling alpha wave, we have suggested that it may reflect a top down controlled process that is associated with an early categorization of the presented stimulus (Klimesch, 2011). Because for words, topographical P1 latencies are sensitive to lexical and (possibly already) semantic features (Zauner et al., 2014), it was suggested that the evoked traveling alpha wave reflects early stages of access to lexical and possibly also semantic memory. The posterior topography may be associated with the processing of visual and graphemic features that enable access to lexical and semantic memory.

Our findings suggest a close association between alpha and the P1. This, however, does not mean that other frequencies do not play an important role. Higher frequencies, particularly in the beta (e.g., Gruber et al., 2005) and gamma range (e.g., Porcaro et al., 2011) may also transiently phase align with alpha within the P1 latency window.

Conclusion

The important conclusion is that single trial fluctuations in amplitude size are associated with fluctuations in the length of alpha periods in the peri- and prestimulus period. These fluctuations modulate alpha frequency and are capable of inducing a traveling trajectory that may be interpreted as a spreading activation process within a neural network that is associated with access to memory. We assume that the topography of prestimulus amplitude increases reflects a top down process that controls the poststimulus flow of spreading activation, as manifested by evoked traveling alpha waves.

References

• 1

  AdrianE. D.MatthewsB. H. (1934). The interpretation of potential waves in the cortex.J. Physiol.81440–471. 10.1113/jphysiol.1934.sp003147

  • CrossRef
  • Google Scholar

• 2

  AlexanderD. M.FlynnG. J.WongW.WhitfordT. J.HarrisA. W.GalletlyC. A.et al (2009). Spatio-temporal EEG waves in first episode schizophrenia.Clin. Neurophysiol.1201667–1682. 10.1016/j.clinph.2009.06.020

  • CrossRef
  • Google Scholar

• 3

  AlexanderD. M.JuricaP.TrengoveC.NikolaevA. R.GepshteinS.ZvyagintsevM.et al (2013). Traveling waves and trial averaging: the nature of single-trial and averaged brain responses in large-scale cortical signals.Neuroimage7395–112. 10.1016/j.neuroimage.2013.01.016

  • CrossRef
  • Google Scholar

• 4

  AlexanderD. M.TrengoveC.WrightJ.BoordP.GordonE. (2006). Measurement of phase gradients in the EEG.J. Neurosci. Methods156111–128. 10.1016/j.jneumeth.2006.02.016

  • CrossRef
  • Google Scholar

• 5

  AmzicaF.SteriadeM. (1997). The K-complex: Its low (<1 Hz) rhythmicity and relation to delta waves.Neurology49952–959. 10.1212/WNL.49.4.952

  • CrossRef
  • Google Scholar

• 6

  AtallahB.ScanzianiM. (2009). Instantaneous modulation of gamma oscillation frequency by balancing excitation with inhibition.Neuron62566–577. 10.1016/j.neuron.2009.04.027

  • CrossRef
  • Google Scholar

• 7

  BahramisharifA.van GervenM. A.AarnoutseE. J.MercierM. R.SchwartzT. H.FoxeJ. J.et al (2013). Propagating neocortical gamma bursts are coordinated by traveling alpha waves.J. Neurosci.3318849–18854. 10.1523/JNEUROSCI.2455-13.2013

  • CrossRef
  • Google Scholar

• 8

  BergerH. (1929). Über das Elektroenkephalogramm des Menschen.Arch. Psychiat. Nervenkr.87527–570. 10.1007/BF01797193

  • CrossRef
  • Google Scholar

• 9

  BishopG. (1933). Cyclic changes in excitability of the optic pathway of the rabbit.Am. J. Physiol.103213–224.

  • Google Scholar

• 10

  BurkittG. R.SilbersteinR. B.CaduschP. J.WoodA. W. (2000). Steady-state visual evoked potentials and traveling waves.Clin. Neurophysiol.111246–258. 10.1016/S1388-2457(99)00194-7

  • CrossRef
  • Google Scholar

• 11

  BuzsákiG. (2006). Rhythms of the Brain.New York: Oxford University Press.10.1093/acprof:oso/9780195301069.001.0001

  • CrossRef
  • Google Scholar

• 12

  BuzsákiG.AnastassiouC.KochC. (2012). The origin of extracellular fields and currents – EEG, ECoG, LFP and spikes.Nat. Rev.13407–420. 10.1038/nrn324

  • CrossRef
  • Google Scholar

• 13

  BuzsákiG.DraguhnA. (2004). Neuronal oscillations in cortical networks.Science3041926–1929. 10.1126/science.1099745

  • CrossRef
  • Google Scholar

• 14

  ClarkV.FanS.HillyardS. (1995). Identification of early visually evoked potential generators by retinotopic and topographic analysis.Hum. Brain Mapp.2170–187. 10.1002/hbm.460020306

  • CrossRef
  • Google Scholar

• 15

  ContrerasD.TimofeevI.SteriadeM. (1996). Mechanisms of long-lasting hyperpolarizations underlying slow sleep oscillations in cat corticothalamic networks.J. Physiol.494251–264. 10.1113/jphysiol.1996.sp021488

  • CrossRef
  • Google Scholar

• 16

  Di RussoF.MartinezA.SerenoM.PitzalisS.HillyardS. (2002). Cortical sources of the early components of the visual evoked potential.Hum. Brain Mapp.1595–111. 10.1002/hbm.10010

  • CrossRef
  • Google Scholar

• 17

  DongW.OlsonE. S. (2008). Supporting evidence for reverse cochlear traveling waves.J. Acoust. Soc. Am.123222–240. 10.1121/1.2816566

  • CrossRef
  • Google Scholar

• 18

  ErmentroutG. B.KleinfeldD. (2001). Traveling electrical waves in cortex: insights from phase dynamics and speculation on a computational role.Neuron2933–44. 10.1016/S0896-6273(01)00178-7

  • CrossRef
  • Google Scholar

• 19

  FellingerR.GruberW.ZaunerA.FreunbergerR.KlimeschW. (2012). Evoked traveling alpha waves predict visual-semantic categorization-speed.Neuroimage593379–3388. 10.1016/j.neuroimage.2011.11.010

  • CrossRef
  • Google Scholar

• 20

  FreemanW. J.BarrieJ. M. (2000). Analysis of spatial patterns of phase in neocortical gamma EEGs in rabbit.J. Neurophysiol.841266–1278.

  • Google Scholar

• 21

  GrattonG.ColesM. G.DonchinE. (1983). A new method for off-line removal of ocular artifact.Electroencephalogr. Clin. Neurophysiol.55468–484. 10.1016/0013-4694(83)90135-9

  • CrossRef
  • Google Scholar

• 22

  GruberW. R.KlimeschW.SausengP.DoppelmayrM. (2005). Alpha phase synchronization predicts P1 and N1 latency and amplitude size.Cereb. Cortex15371–377. 10.1093/cercor/bhh139

  • CrossRef
  • Google Scholar

• 23

  GruberW.ZaunerA.LechingerJ.SchabusM.KutilR.KlimeschW. (2014). Alpha phase, temporal attention, and the generation of early event related potentials.Neuroimage103119–129. 10.1016/j.neuroimage.2014.08.055

  • CrossRef
  • Google Scholar

• 24

  HagmannP.CammounL.GigandetX.MeuliR.HoneyC. J.WedeenV. J.et al (2008). Mapping the structural core of human cerebral cortex.PLoS Biol.6:e159. 10.1371/journal.pbio.0060159

  • CrossRef
  • Google Scholar

• 25

  HallatschekO. (2010). The noisy edge of traveling waves.Proc. Natl. Acad. Sci. U.S.A.1081783–1787. 10.1073/pnas.1013529108

  • CrossRef
  • Google Scholar

• 26

  HindriksR.van PuttenM. J. A. M.DecoG. (2014). Intra-cortical propagation of EEG alpha oscillations.Neuroimage103444–453. 10.1016/j.neuroimage.2014.08.027

  • CrossRef
  • Google Scholar

• 27

  HippJ. F.HawellekD. J.CorbettaM.SiegelM.EngelA. K. (2012). Large-scale cortical correlation structure of spontaneous oscillatory activity.Nat. Neurosci.15884–890. 10.1038/nn.3101

  • CrossRef
  • Google Scholar

• 28

  HughesJ. R. (1995). The phenomenon of traveling waves: a review.Clin. Electroencephalogr.261–6. 10.1177/155005949502600103

  • CrossRef
  • Google Scholar

• 29

  KlimeschW. (2011). Evoked alpha and early access to the knowledge system: the P1 inhibition timing hypothesis.Brain Res.140852–71. 10.1016/j.brainres.2011.06.003

  • CrossRef
  • Google Scholar

• 30

  KlimeschW.HanslmayrS.SausengP.GruberW.DoppelmayrM. (2007a). The P1 and traveling alpha waves: evidence for evoked oscillations.J. Neurophysiol.971311–1318. 10.1152/jn.00876.2006

  • CrossRef
  • Google Scholar

• 31

  KlimeschW.SausengP.AnslmayrS. (2007b). EEG alpha oscillations: the inhibition- timing hypothesis.Brain Res. Rev.5363–88. 10.1016/j.brainresrev.2006.06.003

  • CrossRef
  • Google Scholar

• 32

  KlimeschW.SausengP.HanslmayrS.GruberW.FreunbergerR. (2007c). Event- related phase reorganization may explain evoked neural dynamics.Neurosci. Biobehav. Rev.311003–1016. 10.1016/j.neubiorev.2007.03.005

  • CrossRef
  • Google Scholar

• 33

  MassiminiM.HuberR.FerrarelliF.HillS.TononiG. (2004). The sleep slow oscillation as a traveling wave.J. Neurosci.246862–6870. 10.1523/JNEUROSCI.1318-04.2004

  • CrossRef
  • Google Scholar

• 34

  MazaheriA.JensenO. (2008). Asymmetric amplitude modulations of brain oscillations generate slow evoked responses.J. Neurosci.287781–7787. 10.1523/JNEUROSCI.1631-08.2008

  • CrossRef
  • Google Scholar

• 35

  McRaeK.CreeG. S.SeidenbergM. S.McNorganC. (2005). Semantic feature production norms for a large set of living and nonliving things.Behav. Res. Methods37547–559. 10.3758/BF03192726

  • CrossRef
  • Google Scholar

• 36

  NadasdyZ. (2009). Information encoding and reconstruction from the phase of action potentials.Front. Syst. Neurosci.3:6. 10.3389/neuro.06.006.2009

  • CrossRef
  • Google Scholar

• 37

  NadasdyZ. (2010). Binding by asynchrony: the neuronal phase code.Front. Neurosci.4:e00051. 10.3389/fnins.2010.00051

  • CrossRef
  • Google Scholar

• 38

  NauhausI.BusseL.RingachD. L.CarandiniM. (2012). Robustness of traveling waves in ongoing activity of visual cortex.J. Neurosci.323088–3094. 10.1523/JNEUROSCI.5827-11.2012

  • CrossRef
  • Google Scholar

• 39

  NikulinV. V.Linkenkaer-HansenK.NolteG.LemmS.MüllerK. R.IlmoniemiR. J.et al (2007). A novel mechanism for evoked responses in the human brain.Eur. J. Neurosci.253146–3154. 10.1111/j.1460-9568.2007.05553.x

  • CrossRef
  • Google Scholar

• 40

  NunezP. L. (2000). Toward a quantitative description of large-scale neocortical dynamic function and EEG.Behav. Brain Sci.23371–437. 10.1017/S0140525X00003253

  • CrossRef
  • Google Scholar

• 41

  NunezP. L.SrinivasanR. (2014). Neocortical dynamics due to axon propagation delays in cortico-cortical fibers: EEG traveling and standing waves with implications for top-down influences on local networks and white matter disease.Brain Res.1542138–166. 10.1016/j.brainres.2013.10.036

  • CrossRef
  • Google Scholar

• 42

  PattenT. M.RennieC. J.RobinsonP. A.GongP. (2012). Human cortical traveling waves: dynamical properties and correlations with responses.PLoS ONE7:e38392. 10.1371/journal.pone.0038392

  • CrossRef
  • Google Scholar

• 43

  PetscheH.MarkoA. (1955). Toposkopische Untersuchungen zur Ausbreitung des Alpharhythmus.Wien. Z. Nervenheilk.1287–100.

  • Google Scholar

• 44

  PorcaroC.OstwaldD.HadjipapasA.BarnesG.BagshawA. (2011). The relationship between the visual evoked potential and the gamma band investigated by blind and semi-blind methods.Neuroimage561059–1071. 10.1016/j.neuroimage.2011.03.008

  • CrossRef
  • Google Scholar

• 45

  RomeiV.GrossJ.ThutG. (2010). On the role of prestimulus alpha rhythms over occipito-parietal areas in visual input regulation: correlation or causation?J. Neurosci.308692–8697. 10.1523/JNEUROSCI.0160-10.2010

  • CrossRef
  • Google Scholar

• 46

  RubinoD.RobbinsK. A.HatsopoulosN. G. (2006). Propagating waves mediate information transfer in the motor cortex.Nat. Neurosci.91549–1557. 10.1038/nn1802

  • CrossRef
  • Google Scholar

• 47

  Sanchez-VivesM.McCormickD. (2000). Cellular and network mechanisms of rhythmic recurrent activity in neocortex.Nat. Neurosci.31027–1034. 10.1038/79848

  • CrossRef
  • Google Scholar

• 48

  SilbersteinR. B.SchierM. A.PipingasA.CiorciariJ.WoodS. R.SimpsonD. G. (1990). Steady-state visually evoked potential topography associated with a visual vigilance task.Brain Topogr.3337–347. 10.1007/BF01135443

  • CrossRef
  • Google Scholar

• 49

  SteriadeM. (2006). Grouping of brain rhythms in corticothalamic systems.Neuroscience1371087–1106. 10.1016/j.neuroscience.2005.10.029

  • CrossRef
  • Google Scholar

• 50

  SteriadeM.NunezA.AmzicaF. (1993). A novel slow (<1 Hz) oscillation of neocortical neurons in vivo: depolarizing and hyperpolarizing components.J. Neurosci.133252–3265.

  • Google Scholar

• 51

  TakahashiK.SalehM.PennR. D.HatsopoulosN. (2011). Propagating waves in human motor cortex.Front. Hum. Neurosci.5:40. 10.3389/fnhum.2011.00040

  • CrossRef
  • Google Scholar

• 52

  TraubR.WhittingtonM.CollingS.BuzsakiG.JefferysJ. (1996). Analysis of gamma rhythms in the rat hippocampus in vitro and in vivo.J. Physiol.493471–484. 10.1113/jphysiol.1996.sp021397

  • CrossRef
  • Google Scholar

• 53

  WhittingtonM.TraubR.JefferysJ. G. (1995). Synchronized oscillations in interneuron networks driven by metabotropic glutamate receptor activation.Nature373612–615. 10.1038/373612a0

  • CrossRef
  • Google Scholar

• 54

  WilliamsT. L. (1992). Phase coupling by synaptic spread in chains of coupled neuronal oscillators.Science258662–665. 10.1126/science.1411575

  • CrossRef
  • Google Scholar

• 55

  World Medical Association. (1996). Declaration of Helsinki.Br. Med. J.3131445–1449.

  • Google Scholar

• 56

  WuJ. Y.HuangX.ZhangC. (2008). Propagating waves of activity in the neocortex: what they are, what they do.Neuroscientist14487–502. 10.1177/1073858408317066

  • CrossRef
  • Google Scholar

• 57

  ZaunerA.GruberW.HimmelstoßN.LechingerJ.KlimeschW. (2014). Lexical access and evoked traveling alpha waves.Neuroimage91252–261. 10.1016/j.neuroimage.2014.01.041

  • CrossRef
  • Google Scholar