Johannes Vorwerk

A guideline for head volume conductor modeling in EEG and MEG

For accurate EEG/MEG source analysis it is necessary to model the head volume conductor as realistic as possible. This includes the distinction of the different conductive compartments in the human head. In this study, we investigated the influence of modeling/not modeling the conductive compartments skull spongiosa, skull compacta, cerebrospinal fluid (CSF), gray matter, and white matter and of the inclusion of white matter anisotropy on the EEG/MEG forward solution. Therefore, we created a highly realistic 6-compartment head model with white matter anisotropy and used a state-of-the-art finite element approach. Starting from a 3-compartment scenario (skin, skull, and brain), we subsequently refined our head model by distinguishing one further of the above-mentioned compartments. For each of the generated five head models, we measured the effect on the signal topography and signal magnitude both in relation to a highly resolved reference model and to the model generated in the previous refinement step. We evaluated the results of these simulations using a variety of visualization methods, allowing us to gain a general overview of effect strength, of the most important source parameters triggering these effects, and of the most affected brain regions. Thereby, starting from the 3-compartment approach, we identified the most important additional refinement steps in head volume conductor modeling. We were able to show that the inclusion of the highly conductive CSF compartment, whose conductivity value is well known, has the strongest influence on both signal topography and magnitude in both modalities. We found the effect of gray/white matter distinction to be nearly as big as that of the CSF inclusion, and for both of these steps we identified a clear pattern in the spatial distribution of effects. In comparison to these two steps, the introduction of white matter anisotropy led to a clearly weaker, but still strong, effect. Finally, the distinction between skull spongiosa and compacta caused the weakest effects in both modalities when using an optimized conductivity value for the homogenized compartment. We conclude that it is highly recommendable to include the CSF and distinguish between gray and white matter in head volume conductor modeling. Especially for the MEG, the modeling of skull spongiosa and compacta might be neglected due to the weak effects; the simplification of not modeling white matter anisotropy is admissible considering the complexity and current limitations of the underlying modeling approach.

Investigation of tDCS volume conduction effects in a highly realistic head model.

OBJECTIVE We investigate volume conduction effects in transcranial direct current stimulation (tDCS) and present a guideline for efficient and yet accurate volume conductor modeling in tDCS using our newly-developed finite element (FE) approach. APPROACH We developed a new, accurate and fast isoparametric FE approach for high-resolution geometry-adapted hexahedral meshes and tissue anisotropy. To attain a deeper insight into tDCS, we performed computer simulations, starting with a homogenized three-compartment head model and extending this step by step to a six-compartment anisotropic model. MAIN RESULTS We are able to demonstrate important tDCS effects. First, we find channeling effects of the skin, the skull spongiosa and the cerebrospinal fluid compartments. Second, current vectors tend to be oriented towards the closest higher conducting region. Third, anisotropic WM conductivity causes current flow in directions more parallel to the WM fiber tracts. Fourth, the highest cortical current magnitudes are not only found close to the stimulation sites. Fifth, the median brain current density decreases with increasing distance from the electrodes. SIGNIFICANCE Our results allow us to formulate a guideline for volume conductor modeling in tDCS. We recommend to accurately model the major tissues between the stimulating electrodes and the target areas, while for efficient yet accurate modeling, an exact representation of other tissues is less important. Because for the low-frequency regime in electrophysiology the quasi-static approach is justified, our results should also be valid for at least low-frequency (e.g., below 100 Hz) transcranial alternating current stimulation.

Combining EEG and MEG for the Reconstruction of Epileptic Activity Using a Calibrated Realistic Volume Conductor Model

To increase the reliability for the non-invasive determination of the irritative zone in presurgical epilepsy diagnosis, we introduce here a new experimental and methodological source analysis pipeline that combines the complementary information in EEG and MEG, and apply it to data from a patient, suffering from refractory focal epilepsy. Skull conductivity parameters in a six compartment finite element head model with brain anisotropy, constructed from individual MRI data, are estimated in a calibration procedure using somatosensory evoked potential (SEP) and field (SEF) data. These data are measured in a single run before acquisition of further runs of spontaneous epileptic activity. Our results show that even for single interictal spikes, volume conduction effects dominate over noise and need to be taken into account for accurate source analysis. While cerebrospinal fluid and brain anisotropy influence both modalities, only EEG is sensitive to skull conductivity and conductivity calibration significantly reduces the difference in especially depth localization of both modalities, emphasizing its importance for combining EEG and MEG source analysis. On the other hand, localization differences which are due to the distinct sensitivity profiles of EEG and MEG persist. In case of a moderate error in skull conductivity, combined source analysis results can still profit from the different sensitivity profiles of EEG and MEG to accurately determine location, orientation and strength of the underlying sources. On the other side, significant errors in skull modeling are reflected in EEG reconstruction errors and could reduce the goodness of fit to combined datasets. For combined EEG and MEG source analysis, we therefore recommend calibrating skull conductivity using additionally acquired SEP/SEF data.

Influence of the head model on EEG and MEG source connectivity analyses

The results of brain connectivity analysis using reconstructed source time courses derived from EEG and MEG data depend on a number of algorithmic choices. While previous studies have investigated the influence of the choice of source estimation method or connectivity measure, the effects of the head modeling errors or simplifications have not been studied sufficiently. In the present simulation study, we investigated the influence of particular properties of the head model on the reconstructed source time courses as well as on source connectivity analysis in EEG and MEG. Therefore, we constructed a realistic head model and applied the finite element method to solve the EEG and MEG forward problems. We considered the distinction between white and gray matter, the distinction between compact and spongy bone, the inclusion of a cerebrospinal fluid (CSF) compartment, and the reduction to a simple 3-layer model comprising only the skin, skull, and brain. Source time courses were reconstructed using a beamforming approach and the source connectivity was estimated by the imaginary coherence (ICoh) and the generalized partial directed coherence (GPDC). Our results show that in both EEG and MEG, neglecting the white and gray matter distinction or the CSF causes considerable errors in reconstructed source time courses and connectivity analysis, while the distinction between spongy and compact bone is just of minor relevance, provided that an adequate skull conductivity value is used. Large inverse and connectivity errors are found in the same regions that show large topography errors in the forward solution. Moreover, we demonstrate that the very conservative ICoh is relatively safe from the crosstalk effects caused by imperfect head models, as opposed to the GPDC.

Combined EEG/MEG Can Outperform Single Modality EEG or MEG Source Reconstruction in Presurgical Epilepsy Diagnosis

We investigated two important means for improving source reconstruction in presurgical epilepsy diagnosis. The first investigation is about the optimal choice of the number of epileptic spikes in averaging to (1) sufficiently reduce the noise bias for an accurate determination of the center of gravity of the epileptic activity and (2) still get an estimation of the extent of the irritative zone. The second study focuses on the differences in single modality EEG (80-electrodes) or MEG (275-gradiometers) and especially on the benefits of combined EEG/MEG (EMEG) source analysis. Both investigations were validated with simultaneous stereo-EEG (sEEG) (167-contacts) and low-density EEG (ldEEG) (21-electrodes). To account for the different sensitivity profiles of EEG and MEG, we constructed a six-compartment finite element head model with anisotropic white matter conductivity, and calibrated the skull conductivity via somatosensory evoked responses. Our results show that, unlike single modality EEG or MEG, combined EMEG uses the complementary information of both modalities and thereby allows accurate source reconstructions also at early instants in time (epileptic spike onset), i.e., time points with low SNR, which are not yet subject to propagation and thus supposed to be closer to the origin of the epileptic activity. EMEG is furthermore able to reveal the propagation pathway at later time points in agreement with sEEG, while EEG or MEG alone reconstructed only parts of it. Subaveraging provides important and accurate information about both the center of gravity and the extent of the epileptogenic tissue that neither single nor grand-averaged spike localizations can supply.

Comparison of Boundary Element and Finite Element Approaches to the EEG Forward Problem

The accurate simulation of the electric fields evoked by neural activity is crucial for solving the inverse problem of EEG. Nowadays, boundary element methods (BEM) are frequently applied to achieve this goal, usually relying on the simpli-fication of approximating the human head by three nested compartments with isotropic conductivities (skin, skull, brain). Here, including the highly-conducting cerebrospinal fluid (CSF) is a difficult task due to the complex geomet-rical structure of the CSF, demanding a high number of additional nodes for an accurate modeling and thus a strongly increased computational effort. Though, CSF conductivity is well-known and nearly not varying inter-individually and its significant influence on EEG forward simulation has been shown. The CSF can be included at negligible computational costs when applying finite element (FE) forward approaches. In this study we compare the accuracy and performance of state-of-the-art BE and FE approaches in both artificial and realistic three layer head models, showing that all approaches lead to high numerical accuracies. Furthermore, we give an impression of the significant influence of modeling the CSF compartment as disregarding this compartment leads to model errors that lie clearly above the observed numerical errors.

The role of blood vessels in high-resolution volume conductor head modeling of EEG

Reconstruction of the electrical sources of human EEG activity at high spatiotemporal accuracy is an important aim in neuroscience and neurological diagnostics. Over the last decades, numerous studies have demonstrated that realistic modeling of head anatomy improves the accuracy of source reconstruction of EEG signals. For example, including a cerebrospinal fluid compartment and the anisotropy of white matter electrical conductivity were both shown to significantly reduce modeling errors. Here, we for the first time quantify the role of detailed reconstructions of the cerebral blood vessels in volume conductor head modeling for EEG. To study the role of the highly arborized cerebral blood vessels, we created a submillimeter head model based on ultra-high-field-strength (7 T) structural MRI datasets. Blood vessels (arteries and emissary/intraosseous veins) were segmented using Frangi multi-scale vesselness filtering. The final head model consisted of a geometry-adapted cubic mesh with over 17 × 106 nodes. We solved the forward model using a finite-element-method (FEM) transfer matrix approach, which allowed reducing computation times substantially and quantified the importance of the blood vessel compartment by computing forward and inverse errors resulting from ignoring the blood vessels. Our results show that ignoring emissary veins piercing the skull leads to focal localization errors of approx. 5 to 15 mm. Large errors (>2 cm) were observed due to the carotid arteries and the dense arterial vasculature in areas such as in the insula or in the medial temporal lobe. Thus, in such predisposed areas, errors caused by neglecting blood vessels can reach similar magnitudes as those previously reported for neglecting white matter anisotropy, the CSF or the dura — structures which are generally considered important components of realistic EEG head models. Our findings thus imply that including a realistic blood vessel compartment in EEG head models will be helpful to improve the accuracy of EEG source analyses particularly when high accuracies in brain areas with dense vasculature are required.

Using reciprocity for relating the simulation of transcranial current stimulation to the EEG forward problem.

To explore the relationship between transcranial current stimulation (tCS) and the electroencephalography (EEG) forward problem, we investigate and compare accuracy and efficiency of a reciprocal and a direct EEG forward approach for dipolar primary current sources both based on the finite element method (FEM), namely the adjoint approach (AA) and the partial integration approach in conjunction with a transfer matrix concept (PI). By analyzing numerical results, comparing to analytically derived EEG forward potentials and estimating computational complexity in spherical shell models, AA turns out to be essentially identical to PI. It is then proven that AA and PI are also algebraically identical even for general head models. This relation offers a direct link between the EEG forward problem and tCS. We then demonstrate how the quasi-analytical EEG forward solutions in sphere models can be used to validate the numerical accuracies of FEM-based tCS simulation approaches. These approaches differ with respect to the ease with which they can be employed for realistic head modeling based on MRI-derived segmentations. We show that while the accuracy of the most easy to realize approach based on regular hexahedral elements is already quite high, it can be significantly improved if a geometry-adaptation of the elements is employed in conjunction with an isoparametric FEM approach. While the latter approach does not involve any additional difficulties for the user, it reaches the high accuracies of surface-segmentation based tetrahedral FEM, which is considerably more difficult to implement and topologically less flexible in practice. Finally, in a highly realistic head volume conductor model and when compared to the regular alternative, the geometry-adapted hexahedral FEM is shown to result in significant changes in tCS current flow orientation and magnitude up to 45° and a factor of 1.66, respectively.

Comparison Study for Whitney (Raviart–Thomas)-Type Source Models in Finite-Element-Method-Based EEG Forward Modeling

This study concentrates on finite-element-method (FEM)-based electroencephalography (EEG) forward simulation in which the electric potential evoked by neural activity in the brain is to be calculated at the surface of the head. The main advantage of the FEM is that it allows realistic modeling of tissue conductivity inhomogeneity. However, it is not straightforward to apply the classical model of a dipolar source with the FEM, due to its strong singularity and the resulting irregularity. The focus of this study is on comparing different methods to cope with this problem. In particular, we evaluate the accuracy of Whitney (Raviart-Thomas)-type dipole-like source currents compared to two reference dipole modeling methods: the St. Venant and partial integration approach. Common to all these methods is that they enable direct approximation of the potential field utilizing linear basis functions. In the present context, Whitney elements are particularly interesting, as they provide a simple means to model a divergence-conforming primary current vector field satisfying the square integrability condition. Our results show that a Whitney-type source model can provide simulation accuracy comparable to the present reference methods. It can lead to superior accuracy under optimized conditions with respect to both source location and orientation in a tetrahedral mesh. For random source orientations, the St. Venant approach turns out to be the method of choice over the interpolated version of the Whitney model. The overall moderate differences obtained suggest that practical aspects, such as the focality, should be prioritized when choosing a source model.

A Mixed Finite Element Method to Solve the EEG Forward Problem

Finite element methods have been shown to achieve high accuracies in numerically solving the EEG forward problem and they enable the realistic modeling of complex geometries and important conductive features such as anisotropic conductivities. To date, most of the presented approaches rely on the same underlying formulation, the continuous Galerkin (CG)-FEM. In this article, a novel approach to solve the EEG forward problem based on a mixed finite element method (Mixed-FEM) is introduced. To obtain the Mixed-FEM formulation, the electric current is introduced as an additional unknown besides the electric potential. As a consequence of this derivation, the Mixed-FEM is, by construction, current preserving, in contrast to the CG-FEM. Consequently, a higher simulation accuracy can be achieved in certain scenarios, e.g., when the diameter of thin insulating structures, such as the skull, is in the range of the mesh resolution. A theoretical derivation of the Mixed-FEM approach for EEG forward simulations is presented, and the algorithms implemented for solving the resulting equation systems are described. Subsequently, first evaluations in both sphere and realistic head models are presented, and the results are compared to previously introduced CG-FEM approaches. Additional visualizations are shown to illustrate the current preserving property of the Mixed-FEM. Based on these results, it is concluded that the newly presented Mixed-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches, which motivates a further evaluation of the Mixed-FEM for applications in bioelectromagnetism.