Cancellation of EEG and MEG Signals Generated by Extended and Distributed Sources

doi:10.1002/hbm.20851

Hum Brain Mapp. Author manuscript; available in PMC 2011 January 1.

Published in final edited form as:

Hum Brain Mapp. 2010 January; 31(1): 140–149.

doi: 10.1002/hbm.20851

PMCID: PMC2797557

NIHMSID: NIHMS142166

Cancellation of EEG and MEG Signals Generated by Extended and Distributed Sources

Seppo P. Ahlfors,^1,² Jooman Han,¹ Fa-Hsuan Lin,^1,³ Thomas Witzel,^1,² John W. Belliveau,^1,² Matti S. Hämäläinen,^1,² and Eric Halgren⁴

¹ Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital/Harvard Medical School, Charlestown, MA 02129, USA

² Harvard-MIT Division of Health Sciences and Technology, Cambridge, MA 02135, USA

³ Institute of Biomedical Engineering, National Taiwan University, Taipei, Taiwan

⁴ Department of Radiology, University of California, San Diego, San Diego, CA 92037, USA

Corresponding author: Dr. Seppo P. Ahlfors, Athinoula A. Martinos Center, Massachusetts General Hospital, 149 13^th Street, Rm 2301, Charlestown, MA 02129, Tel: (617) 726-0663, Fax: (617) 726-7422, Email: seppo/at/nmr.mgh.harvard.edu

The publisher's final edited version of this article is available at Hum Brain Mapp

See other articles in PMC that cite the published article.

Abstract

Extracranial patterns of scalp potentials and magnetic fields, as measured with electro- and magnetoencephalography (EEG, MEG), are spatially widespread even when the underlying source in the brain is focal. Therefore, loss in signal magnitude due to cancellation is expected when multiple brain regions are simultaneously active. We characterized these cancellation effects in EEG and MEG using a forward model with sources constrained on an anatomically accurate reconstruction of the cortical surface. Prominent cancellation was found for both EEG and MEG in the case of multiple randomly distributed foci of simulated activity, even when the number of simultaneous dipoles was small. Substantial cancellation occurred also for locally extended patches of simulated activity, when the patches extended to opposite walls of sulci and gyri. For large patches, a difference between EEG and MEG was seen, presumably due to selective cancellation of tangentially vs. radially oriented sources. Cancellation effects can be of importance when electrophysiological data are related to hemodynamic measures. Furthermore, selective cancellation may be used to explain some observed differences between EEG and MEG in terms of focal vs. widespread cortical activity.

Keywords: Magnetoencephalography, Electroencephalography, Forward model, Inverse problem, Cerebral cortex, Cortical patch analysis

INTRODUCTION

Electro- and magnetoencephalography (EEG, MEG) are methods to study electrical brain activity by recording scalp potentials and extra-cranial magnetic fields [Cohen and Halgren, 2003; Gevins et al., 1995; Hämäläinen et al., 1993]. EEG and MEG provide excellent time resolution for noninvasive detection of human brain activity, and estimates of the locations of the underlying sources can be constructed on the basis of the spatial distribution of the EEG and MEG [Baillet et al., 2001]. In part, the spatial resolution in EEG and MEG source estimates is limited by the spatial spread of the signal patterns generated even by focal sources. The spread depends on the distance between the sources and the sensors; the minimum distance is limited by the thickness of the skull and scalp, and, in the case of MEG, the thickness of the Dewar vessel containing the superconducting SQUID sensors. When multiple sources are simultaneously active, their overlapping signal patterns are superimposed, and cancellation of signals of opposite polarity may occur. Cancellation reduces the overall signal-to-noise ratio, thereby contributing to the difficulty of source estimation. In general, when interpreting EEG and MEG data, it is important to take into consideration the overall pattern of simultaneous activation across the brain. This is different from such techniques as functional magnetic resonance imaging (fMRI) in which the signals from individual voxels can be observed largely independently. Here, we attempt to characterize cancellation effects in EEG and MEG signal patterns.

The sources of EEG and MEG signals are described in terms of primary currents [Plonsey, 1969; Tripp, 1983]. The active primary currents give rise to the electric potential distribution, which can measured on the scalp using EEG. In addition, in a conductive medium, there will be associated passive volume currents at locations where the gradient of the electric potential differs from zero. However, the primary currents, together with the conductivity geometry of the head, determine the electric potentials and the volume currents; therefore, it is convenient to express the signals conceptually as functions of the primary currents. The main contribution to the primary currents detectable with EEG and MEG comes from ionic currents in the apical dendrites of cortical pyramidal cells, which are mostly aligned perpendicular to the local cortical surface [Hämäläinen et al., 1993; Lopes da Silva, 1991; Okada et al., 1999]. These currents result from both post-synaptic and active trans-membrane currents [Murakami and Okada, 2006]. For symmetry reasons, the signals from other types of active currents at the cellular level are expected to cancel out macroscopically or generate electric potentials and magnetic fields which fall off rapidly with distance. These invisible currents include those from trans-membrane currents themselves, post-synaptic currents in stellate cells or in the basal dendrites of pyramidal cells, as well as the quadrupole-like current patterns associated with action potentials [Humphrey, 1968a; Humphrey, 1968b; Lorente de No, 1947]. Therefore, it is reasonable to assume that the primary currents underlying the extracranial EEG and MEG signals are located in the cortical gray matter, oriented perpendicular to the cortical surface.

In EEG and MEG source analysis, the basic element for the primary current distribution can be taken to be a current dipole, which summarizes the net effect of the microscopic currents within a volume of a few cubic millimeters, or a few square mm of cortical surface. At the intermediate spatial scale of about 1 cm and above, the folding of the cerebral cortex adds a macroscopic geometric factor, which can be examined using the detailed anatomical information provided by structural MRI. For example, the spatial frequency spectrum of cortical gyrification has been related to the pattern of spatial coherence in scalp EEG [Freeman et al., 2003]. Cancellation of the EEG and MEG signals occurs when the cortical activity extends over a region where the surface normal, and therefore, also the orientation of the source elements, changes. In particular, cancellation occurs when the activation involves opposing walls of a sulcus or a gyrus. Figure 1A depicts simulated EEG and MEG signals from two current dipoles of opposite direction in the primary visual cortex in the calcarine region, which mostly cancel out when active simultaneously. Even at the spatial scale of several centimeters, corresponding to activity distributed over multiple regions across the whole cortex, large amount of cancellation of spatially overlapping EEG or MEG signal patterns can occur. Partial cancellation is illustrated in Fig. 1B with sources in the inferior occipitotemporal region of the left and right hemispheres.

Figure 1

Examples of cancellation of EEG and MEG signals. (A) Simulated data generated by two current dipoles in the primary visual cortex of the left and right hemispheres. The maps depict spatial patterns of signals over the occipital region when each dipole (more ...)

We refer to local patches of activity as extended sources, whereas distributed source patterns may consist of multiple local foci or patches of activity in separate regions of the cortex. This terminology is useful also in the context of source estimation: on one hand, a single equivalent dipole is usually an adequate model, not only for very small foci, but also for moderately extended source patches, such as those corresponding to various types of sensory evoked responses. On the other hand, distributed source models are particularly useful for widely spread activation that consists of multiple simultaneously active functional areas, as is often the case, e.g., with language-related activity. In the present study we explore the cancellation effects in EEG and MEG signals using a realistic reconstruction of the cortical surface [Dale et al., 1999; Dale and Sereno, 1993; Fischl et al., 1999]. For quantitative analysis, we calculated a measure of cancellation for multiple randomly distributed foci of simultaneous activity across the cortex as well as for locally extended patches of activity.

MATERIALS AND METHODS

Forward Model

The forward model can be used to calculate the signals in a set of sensors, due to source elements of unit magnitude. The relationship between the EEG and MEG signals (signal vector v, with the dimension equaling the number of sensors) and the source amplitudes (source vector q, with the dimension equaling the number of dipole elements) is linear: v = A q, where A is the forward matrix often also called the lead-field matrix. The column vectors a_j of A correspond to the signal patterns generated by a single source element of unit amplitude whereas each row of A contains samples of the lead-field or sensitivity pattern of a given sensor. To calculate A, one needs to specify the location and orientation of each source element, the electromagnetic properties of the tissues within the head (in particular, the conductivity geometry), and the location, orientation, and configuration of the sensors. Here we restricted the sources to be located on a realistic reconstruction of the cortical surface, oriented perpendicular to the surface, used the boundary element model (BEM) for the conductivity geometry of the head, and chose commonly used whole-head EEG and MEG sensor array configurations.

The forward model was based on data from a 42-year old male volunteer subject. The procedures were approved by the Massachusetts General Hospital Institutional Review Board. High-resolution structural T1-weighted magnetic resonance images (MRIs) were acquired on a 1.5T Siemens scanner (TR = 2530 ms, TE = 3.25 ms, flip angle = 7°, 128 sagittal slices, slice thickness = 1.3 mm, voxel size = 1.3×1.0×1.3 mm³). A surface representation for the cerebral cortex was constructed from the MRI using the FreeSurfer software [Dale et al., 1999; Fischl et al., 1999]. A mesh of about 300 000 vertices was created, with an average distance of about 1 mm between neighboring source points. Each vertex on the surface corresponding to the gray-white matter border was selected to represent the location of an MEG and EEG source element (current dipole). Selecting the inner surface of the cerebral cortex, rather than, e.g., the middle of the gray matter, helped to avoid numerical instabilities that can occur if the dipoles are very close to the brain-skull interface, where the conductivity change is large. In our forward model, sources were included also in the non-cortical medial regions corresponding to corpus callosum and thalamus that complete the closed surface representation of each hemisphere. The orientation of each source element was assumed to match the local surface normal.

The conductivity geometry of the head was modeled using a boundary element model with three compartments: the brain, the skull, and the scalp. The boundaries were determined from the MRI by identifying the inner and outer surfaces of the skull and the outer surface of the skin [Hamalainen and Sarvas, 1989; Oostendorp and van Oosterom, 1989: Liu, 2002 #15]. Each surface consisted of 5120 vertices; the linear collocation method was used in the BEM calculations [Mosher et al., 1999]. This provides a reasonable accuracy for most cortical source locations [Yvert et al., 1996]. The accuracy of the forward model for the most superficial sources close to the inner surface of the skull may suffer from the sparseness of the BEM vertices [Fuchs et al., 2001]. However, these inaccuracies are expected to have only small effect on the cancellation explored in the present study, except if very large values for the elements of the forward matrix are obtained because of numerical instabilities in the BEM calculations. To exclude this possibility, we checked that there were no outliers (very large values) among the elements of the forward matrix. The conductivity ratios of 1:0.0125:1 were assumed for brain:skull:scalp [Geddes and Baker, 1967].

We used existing MEG and EEG sensor configurations and realistic position based on actual experimental setup, with a 306-channel MEG system (VectorView, Elekta Neuromag) and an EEG cap with 60 scalp electrodes. Of the MEG sensors, 102 were magnetometers, measuring the magnetic field component approximately radial to the surface of the head, and 204 were planar first-order gradiometers, measuring the tangential derivative of the radial component of the magnetic field. The location and orientation of the MEG sensor array with respect to the head was determined with the help of small marker coils approximating magnetic dipoles attached to the head. The position of the marker coils and the EEG electrodes with respect to anatomical landmarks (nasion, pre-auricular points, and several additional points on the scalp and face) were determined with a 3D digitizer (Polhemus, Inc.).

Cancellation index

To quantify cancellation effects, we compared the length of the signal vector generated by n simultaneous sources

(1)

with the sum of the signal vector lengths for the same sources when active individually

(2)

Here a_ij are the elements of the forward matrix, M is the number of sensors. For cortical patch analysis, the forward solution a_j for each source element in Eqs. (1) and (2) was scaled by the area δ_j associated with the source, i.e., a_j was replaced by δ_j a_j. The cancellation index I_C was defined as:

(3)

which ranges from 0 to 1. If the net signal pattern for simultaneous sources is zero, then α = 0 and I_C = 1, corresponding to full (100%) cancellation. If there is no overlap between the signal patterns of the individual sources, then α = β and I_C = 0, indicating absence of cancellation.

Simulated Source Patterns

We calculated the cancellation index for two types of source configurations: (1) individual source elements (current dipoles) distributed randomly across the cortical surface, and (2) patches of uniform activity within the cortex.

For the distributed dipole configuration, the number n of simultaneous source elements ranged from a single dipole up to the total number of vertices in the reconstructed cortex representation. The amplitudes for the source elements were derived from a random zero-mean Gaussian distribution. The cancellation index I_C(n) was estimated by averaging the values obtained for 1000 random selections of n dipoles.

The cortical patches consisted of neighboring source elements within a given distance (“radius”) from the center of the patch. The distances were computed along the cortical surface using the Dijkstra algorithm [Lin et al., 2006]. The radius of a patch was varied from 1 to 50 mm (the largest patches corresponding to n ~ 15,000 dipole elements). Two hundred locations for the patches were chosen randomly on the cortical surface.

To relate the cancellation index for patches with the local curvature properties of the cortex, we computed an orientation nonuniformity index I_O, defined as:

(4)

where

(5)

is the averaged dipole moment vector within a patch, and

(6)

is the area of the patch; δ_j is the area associated with a single vertex and n_j is the corresponding surface normal vector. If all dipole elements within a patch have the same orientation, then I_O = 0. If the dipole moment vectors are oriented such that their sum is zero, then I_O = 1. This is analogous to a “closed source” configuration, in which the net dipole moment cancels out, resulting in minimal measurable signal when distance between the sensors and the patch is large compared with the size of the patch.

RESULTS

Randomly Distributed Dipoles

As shown in Fig. 2, the mean cancellation index I_C for randomly distributed dipoles on the cortical surface increased monotonically as a function of the number of source elements. Notable cancellation occurred even with a small number of sources: I_C reached the value of about 0.5 with only 5 simultaneous dipoles. A small difference between EEG and MEG was seen for n < 100, with I_C being larger for EEG than for MEG. For large values of n, I_C approached unity; e.g., with 10,000 dipoles I_C was ~0.99, suggesting that substantial cancellation of EEG and MEG signals takes place when many sources are simultaneously active.

Figure 2

Cancellation index I_C for randomly distributed sources. The index was computed for an array of 60 EEG electrodes, 102 MEG magnetometers, and 204 MEG gradiometers. The mean and standard deviation over 1000 random selections of n dipoles are shown. For (more ...)

Extended Patches of Activity

An example of EEG and MEG signal patterns generated by uniformly activated cortical patches of two different sizes are shown in Fig. 3. The signal pattern due to a small patch is likely to be similar to that of a single dipole (Fig. 3A), whereas substantial cancellation is expected to occur when the patch is large enough to extend to opposite walls of a sulcus or a gyrus (Fig. 3B and C). To emphasize the cancellation effect per source element, the maps in Fig. 3 were scaled by the patch area. Note, however, that the absolute peak signal magnitudes in this example were larger for the large patch than for the small patch, though only by a factor of about 7 for EEG and 2 for MEG, i.e., much less than 25, which was the ratio of the areas.

Figure 3

Examples of EEG and MEG signals for extended patches of uniform cortical activity. A patch with a 4-mm (A) and 20-mm (B) radius in the left parietal cortex is shown in red on a cortical surface representation reconstructed from anatomical MRI. The dots (more ...)

The mean value of the cancellation index I_C for patches of different sizes is depicted in Fig. 4. For very small patches, I_C was close to zero, consistent with neighboring dipole elements having a similar orientation and consequently the signal patterns summating with little cancellation in either EEG or MEG. For intermediate patches (radii between 2 and 10 mm), I_C increased rapidly as a function of the patch size, reaching I_C ~ 0.5 for a radius of 10 mm; this suggests that the local curvature affected the orientation of the dipole elements within a patch, resulting in cancellation of the signals patterns. For large patches, the cancellation index depended only weakly on the size of the patch, reaching I_C ~ 0.9 for MEG and ~0.8 for EEG for patches with a radius of 50 mm. Thus, in contrast to the case of randomly distributed dipoles (cf. Fig. 2), there appeared to be a prominent difference between EEG and MEG in the cancellation of signals from large patches of uniform cortical activity. Note that it is essential here that all the dipole elements within a patch have the same polarity, corresponding to currents flowing either outward or inward with respect to the cortical surface.

Figure 4

Cancellation index I_C for cortical patches. The mean and standard deviation over 200 randomly located simulated patches of uniform activity for each of the different radii shown.

The relationship between the cancellation index I_C for cortical patches and the corresponding orientation non-uniformity index I_O is shown in Fig. 5. I_C is a measurement of the actual cancellation of the EEG and MEG signals, whereas I_O only takes into account the mean dipole moment vector across the source elements within a patch, ignoring the location of each element. Figure 5 indicates that I_C and I_O were highly correlated. In the case of EEG, the two indices were close to being equal for most of the patches, suggesting that in the case of EEG the cancellation effects can be well approximated by the orientation nonuniformity, i.e., the cancellation was proportional to the net dipole moment obtained by summing the individual dipole moment vectors. For MEG, however, there were many patches for which I_C > I_O, i.e., points that were above the diagonal in Fig. 5. For smaller patches, for which I_C ~ 0.2–0.6, the occasional large values of I_C for MEG can be understood in terms of the orientation of the net dipole moment being radial, resulting in very little MEG signal, i.e., α ~ 0 in Eq. (1). For large values of the indices, the points for MEG clustered around I_C ~ 0.9 and I_O ~ 0.8 instead of the diagonal, suggesting that there was more cancellation than would be predicted simply by the net dipole moment vector. Inspection of the distribution as a function of the patch size indicated that these points correspond to the largest patches. These off-diagonal points found for MEG but not for EEG are presumably related to selective cancellation of tangential as compared with radial source components.

Figure 5

Comparison of the cancellation index I_C and the orientation non-uniformity measure I_O for cortical patches. Data from simulated patches with different radii up to 50 mm are shown for EEG (A), MEG magnetometers (B), and MEG gradiometers (C).

DISCUSSION

We quantified cancellation effects in EEG and MEG signal patterns due to cortical sources using an anatomically accurate model of the cortex. Substantial cancellation was found both for randomly distributed sources and for extended patches of activity. Below we discuss implications of signal cancellation that depends on the distribution and extent of cortical activation to the interpretation of EEG and MEG data, particularly in terms of the integration of multimodal brain imaging data and different characteristics of EEG and MEG.

Relevance of signal cancellation in EEG and MEG studies

Several sources of data suggest that both spontaneous and event-related cortical activity is likely to be spatially extended [Halgren, 2004]. The cortex is highly interconnected, with each neuron receiving and sending ~10,000–100,000 synapses with other neurons. A variety of data, including direct intracranial EEG recordings obtained in subjects undergoing diagnostic procedures prior to surgical therapy for epilepsy [Halgren et al., 1998], indicates that cortical activity can spread widely in less than 30 ms. This suggests that with the possible exception of early sensory evoked responses, patterns of activation are often sufficiently widespread that the effects of cancellation are relevant to interpretation of EEG and MEG data.

Prominent cancellation is expected when the two walls of a sulcus are simultaneously active. Functional areas often cover both sides of a given sulcus. A well-known example is the primary visual cortex, which is organized so that visual stimuli extending across the horizontal meridian will produce sources that cancel across the calcarine fissure. Cancellation may also occur where the opposite sides of the sulcus belong to distinct neurocognitive systems. For example, in the Sylvian fissure, primary and association auditory cortices in the ventral bank are contrasted with somatosensory and other cortices dorsally [Hari et al., 1993]. Also, bilateral activation of regions facing the interhemispheric fissure may produce cancellation (cf. Fig. 1A); such regions include the visual medial occipital, cingulate, and supplementary motor areas [Lang et al., 1991].

An important goal of many MEG and EEG studies is to determine the spatiotemporal patterns of source activity underlying the measured signals. Characterization of cancellation effects, combined with sensitivity analyses [de Jongh et al., 2005; Goldenholz et al., 2009; Hillebrand and Barnes, 2002], can provide valuable insights into what can and cannot be detected with EEG and MEG. The present study was based on the properties of the forward matrix only; thus, the results are not specific to any particular source estimation method. In the calculation of the electric potentials and magnetic fields we assumed a piecewise homogeneous conductor of arbitrary shape and implemented the computations with a boundary-element method, commonly used for EEG and MEG source estimation. More detailed head models that take into account local variations in conductivity and its anisotropy are not expected to change the major patterns of the forward solutions [Haueisen et al., 2002], and thus not to change much the cancellation measures.

The cancellation index for distributed sources is expected to depend somewhat on the design of the sensor array. For example, the overlap of signals, and therefore also the cancellation index I_C, for sources located far apart is likely to be smaller for differential sensors, such as planar MEG gradiometers compared with magnetometers, or bipolar or Laplacian EEG derivation compared with a referential derivation. In general, meaningful comparisons between different sensor array designs require that the signal-to-noise ratio of the measurements is taken into account [Ahonen et al., 1993], which was beyond the scope of the present study. The similarity of the graphs for MEG magnetometers and gradiometers in Fig. 2, however, suggest that the cancellation index is not overly sensitive to the specific sensor configuration.

Also, the exact number of sensors in the array does not have a major effect on the cancellation analysis. For example, the cancellation index in Figs. 2 and and44 were almost identical for MEG gradiometers and MEG magnetometers, even though the number of gradiometers was twice the number of magnetometers. It appears that the arrays used here adequately captured those spatial frequencies that dominate the signal patterns and, therefore, also the cancellation effects. Increasing the sensor density would result in signals that are similar in neighboring sensors; this is expected to affect the cancellation index minimally, as its definition was such that, e.g., doubling the number of sensors by including each channel twice does not alter the value of the index. It is unlikely that the differences seen between EEG and MEG were due to the differences in the number of sensors. As will be discussed below, the main differences in the cancellation effects found between EEG and MEG are consistent with selective cancellation of radial vs. tangential source components, which is independent of the specific sensor arrays used.

Implications for multimodal spatiotemporal brain imaging

Cancellation effects are of concern when electrophysiological EEG and MEG data are compared with hemodynamic measures (fMRI, positron emission tomography, optical imaging) [Dale and Halgren, 2001] or with invasive intracranial recordings [Halgren et al., 1998; Lu and Williamson, 1991]. At the spatial scale of up to a few square millimeters of cortical surface, the net signal from a complex pattern of microscopic neural currents was represented by a single current dipole element in our forward model. This is adequate in practice, as only the dipolar term of the multipole expansion for a complex pattern of microscopic source currents contributes significantly to the recorded potentials and fields when the distance from the source to the sensors is much larger than the extent of the source [Nunez, 1995]. Details of the local current patterns and microscopic cancellation effects are beyond the spatial resolution of presently available non-invasive imaging techniques. Thus, detailed quantitative comparisons with electrophysiological and hemodynamic measures have to rely on information provided by invasive intracranial recordings [Schroeder et al., 1998; Ulbert et al., 2001; Vaughan, 1982] and biophysical modeling [Jones et al., 2007; Murakami and Okada, 2006].

At the intermediate spatial scale of a few square centimeters of cortex, represented here by patches of extended activity, anatomical information about the local folding of the cortical surface can be used to estimate cancellation effects. In addition, information provided by fMRI about the spatial extent of local regions of activity could be incorporated in the cancellation analysis. Assumptions about activity within cortical patches can also be incorporated in inverse modeling [Fuchs et al., 1999; Han et al., 2007; Kincses et al., 1999; Lapalme et al., 2006; Lin et al., 2006].

At the whole-cortex scale, cancellation of EEG and MEG can be large even for randomly distributed simultaneous sources. The reduction in the signal-to-noise ratio per source element due to cancellation may result in some active areas being detected with fMRI but not with EEG or MEG. This appears not to be a serious problem, however, for approaches in which fMRI is used to suggest likely locations for the EEG and MEG sources: simulation studies have indicated that the inclusion of such non-contributing regions as fMRI priors have only little effect on the EEG and MEG source estimates [Ahlfors and Simpson, 2004; Im et al., 2005; Liu et al., 1998].

Differences between EEG and MEG in cancellation of signals from extended sources

The cancellation analysis with extended patches of activity indicated a systematic difference between EEG and MEG that may have interesting implications to the interpretation of combined recordings of EEG and MEG. If a patch of uniform activity extends over multiple sulci and gyri, selective cancellation of signals from sources of particular orientation can occur. The dipole moments of source elements in the two walls of a sulcus point to opposite directions, resulting in strong cancellation of extracranial signals, whereas the sources on the crown of a gyrus and the bottom of a sulcus have the same direction and, thus, produce similar signal patterns that summate with minimal cancellation. In the lateral, dorsal, frontal, and posterior regions of the cerebral cortex, the cancelling sources in sulcal walls are mainly tangentially oriented, whereas the non-cancelling ones are radial. Since MEG is insensitive to radially oriented sources [Grynszpan and Geselowitz, 1973], strong cancellation is expected for MEG, whereas the radial sources will still generate prominent EEG signals [Eulitz et al., 1997; Freeman et al., 2009].

The roles of tangential and radial sources are interchanged in the ventral and medial parts of the cortex. For these regions, the main orientation of the non-opposing gyral source elements is tangential, whereas the cancelling sources in the sulcal walls have a predominantly radial orientation. Therefore, signals from extended sources in these regions should be present in both EEG and MEG. However, as mentioned above, bilateral activity of functional areas in the medial wall of the cerebral hemispheres is likely to result in strong cancellation [Lang et al., 1991].

Selective cancellation of radial vs. tangential sources could underlie some of the differences seen between EEG and MEG in cognitive event-related responses [Eulitz et al., 1997]. For example, the amplitude of the P300 component is typically large in central scalp EEG electrodes, whereas in MEG the P300 is often small and scattered across sensors [Siedenberg et al., 1996]. This is consistent with signal patterns generated by large patches of activity, for which the EEG is expected to mostly reflect the superficial radial sources on the top of the gyri, but the corresponding MEG only reflects the residual signals remaining from incomplete cancellation of tangential sources. A widely distributed pattern of cortical activity related to the P300 has also been observed in intracortical recordings [Halgren et al., 1998] and in fMRI [Kiehl and Liddle, 2003]. The difference in the selective cancellation of EEG and MEG signals could also be helpful in identifying the presence of extended regions of synchronized oscillatory activity [Freeman et al., 2009].

Selective cancellation of spontaneous background activity is also likely to contribute to differences in the SNR of event-related activity in EEG and MEG. For example, inter-ictal epileptogenic discharges are sometimes detected either in EEG or in MEG but not necessarily in both [Iwasaki et al., 2005; Knake et al., 2006; Lin et al., 2003]; explanation for this may lie in the combination of differential sensitivity of EEG vs. MEG to the source of interest [de Jongh et al., 2005; Goldenholz et al., 2009; Hillebrand and Barnes, 2002] as well as in the selective cancellation of widespread background brain activity [de Munck et al., 1992]. These considerations support the view that EEG and MEG provide complementary information of human brain activity [Cohen and Cuffin, 1983; Hari and Ilmoniemi, 1986; Wood et al., 1985] and that these two closely related measures should be combined for optimal results. Both experimental [Sharon et al., 2007] and theoretical [Molins et al., 2008] results have suggested that the accuracy of source estimates can improve beyond the capabilities of EEG or MEG alone by combining data acquired with the two methods.

CONCLUSIONS

Substantial amount of cancellation can occur in extracranial EEG and MEG signals generated by simultaneously active sources. For extended patches of activity, selective cancellation of signals due to sulcal vs. gyral sources can result in a systematic difference between EEG and MEG. These effects are independent of the source estimation method employed and can be important in explaining differences in experimental data obtained with EEG and MEG.

Acknowledgments

This work was supported by Whitaker foundation grant RG-01-0294 and by NIH grants NS57500, NS44623, and NS18741. This work was supported in part by The National Center for Research Resources (P41RR14075) and the MIND Research Network (MRN). We thank David Cohen, Anders Dale, Ksenija Marinkovic, Maria Mody, and Suresh Narayanan for useful discussions and comments.

References

Ahlfors SP, Simpson GV. Geometrical interpretation of fMRI-guided MEG/EEG inverse estimates. Neuroimage. 2004;22:323–32. [PubMed]
Ahonen AI, Hamalainen MS, Ilmoniemi RJ, Kajola MJ, Knuutila JE, Simola JT, Vilkman VA. Sampling theory for neuromagnetic detector arrays. IEEE Trans Biomed Eng. 1993;40:859–69. [PubMed]
Baillet S, Mosher JC, Leahy RM. Electromagnetic brain mapping. IEEE Signal Proc Mag. 2001;18:14–30.
Cohen D, Cuffin BN. Demonstration of useful differences between magnetoencephalogram and electroencephalogram. Electroencephalogr Clin Neurophysiol. 1983;56:38–51. [PubMed]
Cohen D, Halgren E. Magnetoencephalograpy (Neuromagnetism) In: Adelman G, Smith BH, editors. Encyclopedia of Neuroscience. 3. Amsterdam: Elsevier; 2003. in press.
Dale AM, Fischl B, Sereno MI. Cortical surface-based analysis. I Segmentation and surface reconstruction. NeuroImage. 1999;9:179–94. [PubMed]
Dale AM, Halgren E. Spatiotemporal mapping of brain activity by integration of multiple imaging modalities. Curr Opin Neurobiol. 2001;11:202–8. [PubMed]
Dale AM, Sereno MI. Improved localization of cortical activity by combining EEG and MEG with MRI cortical surface reconstruction: A linear approach. J Cogn Neurosci. 1993;5:162–176.
de Jongh A, de Munck JC, Goncalves SI, Ossenblok P. Differences in MEG/EEG epileptic spike yields explained by regional differences in signal-to-noise ratios. J Clin Neurophysiol. 2005;22:153–8. [PubMed]
de Munck JC, Vijn PC, Lopes da Silva FH. A random dipole model for spontaneous brain activity. IEEE Trans Biomed Eng. 1992;39:791–804. [PubMed]
Eulitz C, Eulitz H, Elbert T. Differential outcomes from magneto- and electroencephalography for the analysis of human cognition. Neurosci Lett. 1997;227:185–8. [PubMed]
Fischl B, Sereno MI, Dale AM. Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system. NeuroImage. 1999;9:195–207. [PubMed]
Freeman WJ, Ahlfors SP, Menon V. Combining fMRI with EEG and MEG in order to relate patterns of brain activity to cognition. Int J Psychophysiol 2009 [PMC free article] [PubMed]
Freeman WJ, Holmes MD, Burke BC, Vanhatalo S. Spatial spectra of scalp EEG and EMG from awake humans. Clin Neurophysiol. 2003;114:1053–68. [PubMed]
Fuchs M, Wagner M, Kastner J. Boundary element method volume conductor models for EEG source reconstruction. Clin Neurophysiol. 2001;112:1400–7. [PubMed]
Fuchs M, Wagner M, Kohler T, Wischmann HA. Linear and nonlinear current density reconstructions. J Clin Neurophysiol. 1999;16:267–95. [PubMed]
Geddes LA, Baker LE. The specific resistance of biological materials - a compendium of data for the biomedical engineer and physiologists. Med Biol Eng. 1967;5:271–293. [PubMed]
Gevins A, Leong H, Smith ME, Le J, Du R. Mapping cognitive brain function with modern high-resolution electroencephalography. Trends Neurosci. 1995;18:429–36. [PubMed]
Goldenholz DM, Ahlfors SP, Hamalainen MS, Sharon D, Ishitobi M, Vaina LM, Stufflebeam SM. Mapping the signal-to-noise-ratios of cortical sources in magnetoencephalography and electroencephalography. Hum Brain Mapp. 2009;30:1077–86. [PMC free article] [PubMed]
Grynszpan F, Geselowitz DB. Model studies of the magnetocardiogram. Biophys J. 1973;13:911–25. [PubMed]
Halgren E. How can intracranial recordings assist MEG source localization? Neurol Clin Neurophysiol. 2004;2004:86. [PubMed]
Halgren E, Marinkovic K, Chauvel P. Generators of the late cognitive potentials in auditory and visual oddball tasks. Electroencephalogr Clin Neurophysiol. 1998;106:156–64. [PubMed]
Hämäläinen MS, Hari R, Ilmoniemi RJ, Knuutila J, Lounasmaa OV. Magnetoencephalography - theory, instrumentation, and applications to noninvasive studies of the working human brain. Reviews of Modern Physics. 1993;65:413–497.
Hamalainen MS, Sarvas J. Realistic conductivity geometry model of the human head for interpretation of neuromagnetic data. IEEE Trans Biomed Eng. 1989;36:165–171. [PubMed]
Han J, Kim JS, Chung CK, Park KS. Evaluation of smoothing in an iterative lp-norm minimization algorithm for surface-based source localization of MEG. Phys Med Biol. 2007;52:4791–803. [PubMed]
Hari R, Ilmoniemi RJ. Cerebral magnetic fields. Crit Rev Biomed Eng. 1986;14:93–126. [PubMed]
Hari R, Karhu J, Hamalainen M, Knuutila J, Salonen O, Sams M, Vilkman V. Functional organization of the human first and second somatosensory cortices: a neuromagnetic study. Eur J Neurosci. 1993;5:724–34. [PubMed]
Haueisen J, Tuch DS, Ramon C, Schimpf PH, Wedeen VJ, George JS, Belliveau JW. The influence of brain tissue anisotropy on human EEG and MEG. Neuroimage. 2002;15:159–66. [PubMed]
Hillebrand A, Barnes GR. A quantitative assessment of the sensitivity of whole-head MEG to activity in the adult human cortex. Neuroimage. 2002;16:638–50. [PubMed]
Humphrey DR. Re-analysis of the antidromic cortical response. I. Potentials evoked by stimulation of the isolated pyramidal tract. Electroencephalogr Clin Neurophysiol. 1968a;24:116–29. [PubMed]
Humphrey DR. Re-analysis of the antidromic cortical response. II. On the contribution of cell discharge and PSPs to the evoked potentials. Electroencephalogr Clin Neurophysiol. 1968b;25:421–42. [PubMed]
Im CH, Jung HK, Fujimaki N. fMRI-constrained MEG source imaging and consideration of fMRI invisible sources. Hum Brain Mapp. 2005;26:110–8. [PubMed]
Iwasaki M, Pestana E, Burgess RC, Luders HO, Shamoto H, Nakasato N. Detection of epileptiform activity by human interpreters: blinded comparison between electroencephalography and magnetoencephalography. Epilepsia. 2005;46:59–68. [PubMed]
Jones SR, Pritchett DL, Stufflebeam SM, Hamalainen M, Moore CI. Neural correlates of tactile detection: a combined magnetoencephalography and biophysically based computational modeling study. J Neurosci. 2007;27:10751–64. [PMC free article] [PubMed]
Kiehl KA, Liddle PF. Reproducibility of the hemodynamic response to auditory oddball stimuli: a six-week test-retest study. Hum Brain Mapp. 2003;18:42–52. [PubMed]
Kincses WE, Braun C, Kaiser S, Elbert T. Modeling extended sources of event-related potentials using anatomical and physiological constraints. Hum Brain Mapp. 1999;8:182–93. [PubMed]
Knake S, Halgren E, Shiraishi H, Hara K, Hamer HM, Grant PE, Carr VA, Foxe D, Camposano S, Busa E, Witzel T, Hamalainen MS, Ahlfors SP, Bromfield EB, Black PM, Bourgeois BF, Cole AJ, Cosgrove GR, Dworetzky BA, Madsen JR, Larsson PG, Schomer DL, Thiele EA, Dale AM, Rosen BR, Stufflebeam SM. The value of multichannel MEG and EEG in the presurgical evaluation of 70 epilepsy patients. Epilepsy Res. 2006;69:80–6. [PubMed]
Lang W, Cheyne D, Kristeva R, Beisteiner R, Lindinger G, Deecke L. Three-dimensional localization of SMA activity preceding voluntary movement. A study of electric and magnetic fields in a patient with infarction of the right supplementary motor area. Exp Brain Res. 1991;87:688–95. [PubMed]
Lapalme E, Lina JM, Mattout J. Data-driven parceling and entropic inference in MEG. Neuroimage. 2006;30:160–71. [PubMed]
Lin FH, Belliveau JW, Dale AM, Hamalainen MS. Distributed current estimates using cortical orientation constraints. Hum Brain Mapp. 2006;27:1–13. [PubMed]
Lin YY, Shih YH, Hsieh JC, Yu HY, Yiu CH, Wong TT, Yeh TC, Kwan SY, Ho LT, Yen DJ, Wu ZA, Chang MS. Magnetoencephalographic yield of interictal spikes in temporal lobe epilepsy. Comparison with scalp EEG recordings. Neuroimage. 2003;19:1115–26. [PubMed]
Liu AK, Belliveau JW, Dale AM. Spatiotemporal imaging of human brain activity using fMRI constrained MEG data: Monte Carlo simulations. Proceedings of the National Academy of Sciences of the United States of America. 1998;95:8945–8950. [PubMed]
Lopes da Silva F. Neural mechanisms underlying brain waves: From neural membranes to networks. Electroencephalography and Clinical Neurophysiology. 1991;79:81–93. [PubMed]
Lorente de No R. A study of nerve physiology. Studies of the Rockefeller Institute. 1947;132 ch 16.
Lu ZL, Williamson SJ. Spatial extent of coherent sensory-evoked cortical activity. Exp Brain Res. 1991;84:411–6. [PubMed]
Molins A, Stufflebeam SM, Brown EN, Hamalainen MS. Quantification of the benefit from integrating MEG and EEG data in minimum l2-norm estimation. Neuroimage. 2008;42:1069–77. [PubMed]
Mosher JC, Leahy RM, Lewis PS. EEG and MEG: forward solutions for inverse methods. IEEE Trans Biomed Eng. 1999;46:245–59. [PubMed]
Murakami S, Okada Y. Contributions of principal neocortical neurons to magnetoencephalography and electroencephalography signals. J Physiol. 2006;575:925–36. [PubMed]
Nunez P. Neocortical Dynamics and Human EEG Rhythms. New York: Oxford University Press; 1995.
Okada Y, Lahteenmaki A, Xu C. Comparison of MEG and EEG on the basis of somatic evoked responses elicited by stimulation of the snout in the juvenile swine. Clin Neurophysiol. 1999;110:214–29. [PubMed]
Oostendorp TF, van Oosterom A. Source parameter estimation in inhomogeneous volume conductors of arbitrary shape. IEEE Trans Biomed Engr. 1989;36:382–391. [PubMed]
Plonsey R. Bioelectric Phenomena. New York: McGraw-Hill; 1969.
Schroeder CE, Mehta AD, Givre SJ. A spatiotemporal profile of visual system activation revealed by current source density analysis in the awake macaque. Cereb Cortex. 1998;8:575–92. [PubMed]
Sharon D, Hamalainen MS, Tootell RB, Halgren E, Belliveau JW. The advantage of combining MEG and EEG: comparison to fMRI in focally stimulated visual cortex. Neuroimage. 2007;36:1225–35. [PMC free article] [PubMed]
Siedenberg R, Goodin DS, Aminoff MJ, Rowley HA, Roberts TP. Comparison of late components in simultaneously recorded event-related electrical potentials and event-related magnetic fields. Electroencephalogr Clin Neurophysiol. 1996;99:191–7. [PubMed]
Tripp JH. Physical concepts and mathematical models. In: Williamson SJ, Romani G-L, Kaufman L, Modena I, editors. Biomagnetism: An Interdisciplinary Approach. New York: Plenum Press; 1983. pp. 101–139.
Ulbert I, Karmos G, Heit G, Halgren E. Early discrimination of coherent versus incoherent motion by multiunit and synaptic activity in human putative MT+ Hum Brain Mapp. 2001;13:226–38. [PubMed]
Vaughan HG., Jr The neural origins of human event-related potentials. Ann N Y Acad Sci. 1982;388:125–38. [PubMed]
Wood CC, Cohen D, Cuffin BN, Yarita M, Allison T. Electrical sources in human somatosensory cortex: identification by combined magnetic and potential recordings. Science. 1985;227:1051–3. [PubMed]
Yvert B, Bertrand O, Echallier JF, Pernier J. Improved dipole localization using local mesh refinement of realistic head geometries: an EEG simulation study. Electroencephalogr Clin Neurophysiol. 1996;99:79–89. [PubMed]