Recognition: unknown
Forward--Inverse Interplay in FEM-Based EEG Source Imaging: Distributional Signatures of Advanced Source Models and Inverse Solvers
Pith reviewed 2026-05-10 00:02 UTC · model grok-4.3
The pith
The success of an inverse solver in EEG source imaging depends on how closely the forward source model matches the solver's assumptions about the spatial distribution of brain activity.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Advanced source models implemented via the finite element method produce unique distributional signatures in inverse solutions, with a clear dependence such that point-like source models align successfully with inverse solvers presupposing single-point activity, as quantified by Earth Mover's Distance, depth bias scatter plots, amplitude distributions, and focality measures.
What carries the argument
Distributional quantitative measures such as Earth Mover's Distance and depth bias scatter plots that capture how different source model implementations affect the spatial characteristics of inverse reconstructions.
Load-bearing premise
The chosen distributional metrics and the specific implementations of the divergence-conforming and local subtraction models are representative enough to reveal the general forward-inverse interplay.
What would settle it
Repeating the experiments on real measured EEG data with independently known source locations to check if the dependence between source model and inverse success rate persists.
Figures
read the original abstract
Electroencephalography (EEG) source imaging aims to infer brain activity from electrical potentials measured on the scalp. This is a difficult problem because many different source patterns can explain the same measurements. The result depends strongly on two things: the forward model and the inverse method. In this work, we study how these two parts work together. We focus not only on where the activity is located, but also on how the reconstructed activity is distributed in space. We suggest that different source models create different signatures in the reconstructed activity. We use realistic head models and compute forward solutions with the finite element method using Zeffiro Interface and DUNEuro. We test different source models, including 2 implementations of a divergence-conforming model, and one implementation of Local subtraction approach. For inverse methods, we use advanced methods such as standardized hierarchical adaptive L1 regression (SHAL1R), standardized Kalman filtering (SKF), and classical dipole scanning. To understand the complex interplay between the forward and inverse approaches, we analyze the inverse source localization results using distributional quantitative measures, including Earth Mover's Distance and depth bias scatter plot, and qualitatively assess the amplitude distribution and focality. The results show that there is a strong dependence between the choice of source model and the success rate of a given inverse method: a source model that corresponds well with a single point-like source is a good match with an inverse method that presupposes such a source.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a simulation study of the forward-inverse interplay in FEM-based EEG source imaging. Using realistic head models computed via Zeffiro Interface and DUNEuro, it compares three source-model variants (two divergence-conforming implementations and one local-subtraction approach) against three inverse solvers (SHAL1R, SKF, and classical dipole scanning). Reconstruction quality is assessed with distributional metrics including Earth Mover's Distance, depth-bias scatter plots, amplitude distributions, and focality. The central claim is that a strong dependence exists between source-model choice and inverse-method performance, with point-like forward models pairing effectively with point-assuming inverses.
Significance. If the reported pairings hold under scrutiny, the work usefully demonstrates that distributional signatures can reveal assumption-matching effects between forward and inverse components, moving beyond scalar localization error. The empirical use of multiple advanced source-model implementations and quantitative metrics (EMD, depth bias) provides concrete evidence for the interplay in a controlled FEM setting. This could guide model selection in numerical EEG pipelines, though the simulation-only design limits direct translation to measured data.
major comments (2)
- [Methods] Methods section: the simulation protocol does not specify the number of Monte Carlo realizations, exact noise models, source-location sampling strategy, or quantitative differences between the two divergence-conforming implementations. Without these details it is impossible to judge whether the observed dependence is robust or sensitive to post-hoc choices.
- [Results] Results and Discussion: the claim of 'strong dependence' and differential 'success rates' rests on internal consistency within each source-model family; the manuscript should add a sensitivity test using source distributions that deviate from all tested models (e.g., extended patches or multi-dipole configurations) to show that the reported ranking is not an artifact of the closed simulation loop.
minor comments (2)
- [Abstract] Abstract: the phrasing '2 implementations of a divergence-conforming model' is repeated without clarifying their numerical distinctions (e.g., basis-function choice or boundary handling); a single sentence of differentiation would improve clarity.
- [Results] Figure captions and text: the depth-bias scatter plots and EMD values are presented without error bars or statistical comparison tests; adding these would strengthen the quantitative claims.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We address each major comment below and will revise the manuscript to incorporate the suggested improvements for clarity and robustness.
read point-by-point responses
-
Referee: [Methods] Methods section: the simulation protocol does not specify the number of Monte Carlo realizations, exact noise models, source-location sampling strategy, or quantitative differences between the two divergence-conforming implementations. Without these details it is impossible to judge whether the observed dependence is robust or sensitive to post-hoc choices.
Authors: We agree that these details are essential for reproducibility and for allowing readers to assess the robustness of the reported forward-inverse dependence. In the revised manuscript we will expand the Methods section to specify the number of Monte Carlo realizations performed, the precise noise model (including SNR and distribution), the source-location sampling strategy, and the quantitative differences between the two divergence-conforming implementations (e.g., basis-function order and stabilization techniques). revision: yes
-
Referee: [Results] Results and Discussion: the claim of 'strong dependence' and differential 'success rates' rests on internal consistency within each source-model family; the manuscript should add a sensitivity test using source distributions that deviate from all tested models (e.g., extended patches or multi-dipole configurations) to show that the reported ranking is not an artifact of the closed simulation loop.
Authors: We acknowledge that the current results are obtained within the tested source-model families and that an external sensitivity test would further strengthen the claim. In the revised manuscript we will add a dedicated sensitivity subsection that includes simulations with extended source patches and multi-dipole configurations, reporting the corresponding distributional metrics (EMD, depth-bias, focality) to demonstrate that the observed model-matching effects are not an artifact of the closed loop. revision: yes
Circularity Check
No circularity: empirical simulation study with independent benchmarks
full rationale
The paper reports results from FEM-based numerical simulations comparing source models (divergence-conforming variants, local subtraction) and inverse solvers (SHAL1R, SKF, dipole scanning) using metrics such as Earth Mover's Distance and depth-bias plots. The central observation of dependence between source-model choice and inverse success is an empirical finding from these controlled experiments, not a mathematical derivation. Ground-truth sources are generated within the simulation framework, but the study explicitly benchmarks against classical dipole scanning as an external reference method. No equations are presented that reduce predictions to fitted parameters defined from the same data, and no load-bearing self-citations or uniqueness theorems are invoked to force the conclusions. The work is self-contained as a comparative simulation study.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Finite element discretization of the EEG forward problem on realistic head models accurately represents volume conduction.
- domain assumption The distributional measures (Earth Mover's Distance, depth bias) are appropriate proxies for reconstruction quality and source model suitability.
Reference graph
Works this paper leans on
-
[1]
Standardized low-resolution brain electromagnetic tomography (sloreta),
R. D. Pascual-Marqui, “Standardized low-resolution brain electromagnetic tomography (sloreta),”Methods and Findings in Experimental and Clinical Pharmacology, vol. 24, pp. 5–12, 2002
2002
-
[2]
Standardized shrinking LORETA-FOCUSS (SSLOFO): A new algorithm for spatio-temporal EEG source reconstruction,
H. Liu, P. Schimpf, G. Dong, X. Gao, F. Yang, and S. Gao, “Standardized shrinking LORETA-FOCUSS (SSLOFO): A new algorithm for spatio-temporal EEG source reconstruction,” eng,IEEE Trans. Biomed. Eng., vol. 52, no. 10, pp. 1681–1691, 2005,ISSN: 0018-9294
2005
-
[3]
Standardized hierar- chical adaptive lp regression for noise robust focal epilepsy source reconstructions,
J. Lahtinen, A. Koulouri, S. Rampp, J. Wellmer, C. Wolters, and S. Pursiainen, “Standardized hierar- chical adaptive lp regression for noise robust focal epilepsy source reconstructions,”Clinical Neurophysiol- ogy, vol. 159, pp. 24–40, 2024
2024
-
[4]
Standardized kalman filtering for dynamical source localization of concurrent subcortical and cortical brain activity,
J. Lahtinen, P. Ronni, N. P. Subramaniyam, A. Koulouri, C. Wolters, and S. Pursiainen, “Standardized kalman filtering for dynamical source localization of concurrent subcortical and cortical brain activity,”Clinical Neuro- physiology, vol. 168, pp. 15–24, 2024
2024
-
[5]
Discrete, 3D distributed, linear imaging methods of electric neuronal activity. Part 1: Exact, zero error localization,
R. D. Pascual-Marqui, “Discrete, 3D distributed, linear imaging methods of electric neuronal activity. Part 1: Exact, zero error localization,” eng,arXiv.org, 2007, ISSN: 2331-8422
2007
-
[6]
J. Lahtinen, “On bias and its reduction via standardization in discretized electromagnetic source localization prob- lems,”Inverse Problems, vol. 40, no. 9, p. 095 002, Jul. 2024.DOI: 10.1088/1361-6420/ad5f53
-
[7]
The effects of peeling on finite element method- based eeg source reconstruction,
S. Söderholm, J. Lahtinen, C. H. Wolters, and S. Pursi- ainen, “The effects of peeling on finite element method- based eeg source reconstruction,”Biomedical Signal Pro- cessing and Control, vol. 89, p. 105 695, 2024
2024
-
[8]
Zeffiro interface for electromagnetic brain imaging,
Q. He, A. Rezaei, and S. Pursiainen, “Zeffiro interface for electromagnetic brain imaging,”Neuroinformatics, vol. 18, no. 2, pp. 237–250, 2020
2020
-
[9]
Duneuro: A software toolbox for forward modeling in bioelectromagnetism,
S. Schrader, A. Nüßing, C. Engwer, J. V orwerk, and C. H. Wolters, “Duneuro: A software toolbox for forward modeling in bioelectromagnetism,”PLOS ONE, vol. 16, no. 6, e0252431, 2021
2021
-
[10]
A realistic, accurate and fast source modeling approach for the eeg forward problem,
T. Miinalainen et al., “A realistic, accurate and fast source modeling approach for the eeg forward problem,” NeuroImage, vol. 184, pp. 56–67, 2019
2019
-
[11]
The local subtraction ap- proach for eeg and meg forward modeling,
M. B. Höltershinken et al., “The local subtraction ap- proach for eeg and meg forward modeling,”SIAM Journal on Scientific Computing, vol. 47, no. 1, B160–B189, 2025
2025
-
[12]
A. A. Vergani, “Hans berger (1873–1941): The german psychiatrist who recorded the first electrical brain sig- nal in humans 100 years ago,”Advances in Physiology Education, vol. 48, no. 4, pp. 878–881, 2024, PMID: 39236103.DOI: 10.1152/advan.00119.2024
-
[13]
F. Missey et al., “Non-invasive temporal interference stimulation of the hippocampus suppresses epileptic biomarkers in patients with epilepsy: Biophysical differ- ences between kilohertz and amplitude modulated stim- ulation,”Brain Stimulation, vol. 19, no. 1, p. 102 981, 2026,ISSN: 1935-861X.DOI: 10.1016/j.brs.2025.11.008
-
[14]
T. R. Knösche and J. Haueisen,EEG/MEG Source Recon- struction, Texbook for Electro- and Magnetoencephalog- raphy. Springer, 2022,ISBN: 978-3-030-74918-7
2022
-
[15]
R. Hari and A. Puce,MEG - EEG Primer. Oxford University Press, Sep. 2023,ISBN: 9780197542187.DOI: 10.1093/med/9780197542187.001.0001
-
[16]
A. Chaddad, Y . Wu, R. Kateb, and A. Bouridane, “Elec- troencephalography signal processing: A comprehensive review and analysis of methods and techniques.,”Sensors (Basel, Switzerland), vol. 23, no. 14, 2023,ISSN: 1424- 8220 (Electronic).DOI: 10.3390/s23146434
-
[17]
Icbm 152 nonlinear atlases,
NeuroImaging and Surgical Technologies Lab. “Icbm 152 nonlinear atlases,” Accessed: Apr. 16, 2026. [Online]. Available: https://nist.mni.mcgill.ca/icbm-152-nonlinear- atlases-2009/
2026
-
[18]
A mathematical evaluation of the core conductor model,
J. Clark and R. Plonsey, “A mathematical evaluation of the core conductor model,”Biophysical Journal, vol. 6, no. 1, pp. 95–112, 1966,ISSN: 0006-3495.DOI: 10.1016/ S0006-3495(66)86642-0
1966
-
[19]
A minimum entropy principle in the gas dynamics equations
M. Scherg and D. V on Cramon, “Two bilateral sources of the late aep as identified by a spatio-temporal dipole model,”Electroencephalography and Clinical Neuro- physiology/Evoked Potentials Section, vol. 62, no. 1, pp. 32–44, 1985,ISSN: 0168-5597.DOI: 10.1016/0168- 5597(85)90033-4
-
[20]
Magnetoencephalography— theory, instrumentation, and applications to noninvasive studies of the working human brain,
M. S. Hämäläinen, R. Hari, R. J. Ilmoniemi, J. Knuu- tila, and O. V . Lounasmaa, “Magnetoencephalography— theory, instrumentation, and applications to noninvasive studies of the working human brain,”Reviews of Modern Physics, vol. 65, no. 2, pp. 413–497, 1993
1993
-
[21]
Location of sources of evoked scalp potentials: Corrections for skull and scalp thicknesses,
J. P. Ary, S. A. Klein, and D. H. Fender, “Location of sources of evoked scalp potentials: Corrections for skull and scalp thicknesses,”IEEE Transactions on Biomedical Engineering, vol. BME-28, no. 6, pp. 447–452, 1981. DOI: 10.1109/TBME.1981.324817
-
[22]
Experimental tests of eeg source localization accuracy in spherical head models,
B. N. Cuffin, D. L. Schomer, J. R. Ives, and H. Blume, “Experimental tests of eeg source localization accuracy in spherical head models,”Clinical Neurophysiology, vol. 112, no. 1, pp. 46–51, 2001
2001
-
[23]
Experimental tests of eeg source localization accuracy in realistically shaped head models,
B. N. Cuffin, D. L. Schomer, J. R. Ives, and H. Blume, “Experimental tests of eeg source localization accuracy in realistically shaped head models,”Clinical Neurophys- iology, vol. 112, no. 12, pp. 2288–2292, 2001
2001
-
[24]
J. de Munck, “A linear discretization of the volume con- ductor boundary integral equation using analytically in- tegrated elements (electrophysiology application),”IEEE Transactions on Biomedical Engineering, vol. 39, no. 9, pp. 986–990, 1992.DOI: 10.1109/10.256433
-
[25]
Evaluation of boundary element methods for the eeg forward problem: Effect of linear interpolation,
H. Schlitt, L. Heller, R. Aaron, E. Best, and D. Ranken, “Evaluation of boundary element methods for the eeg forward problem: Effect of linear interpolation,”IEEE Transactions on Biomedical Engineering, vol. 42, no. 1, pp. 52–58, 1995.DOI: 10.1109/10.362919
-
[26]
C. H. Wolters, H. Köstler, C. Möller, J. Härdtlein, L. Grasedyck, and W. Hackbusch, “Numerical mathematics of the subtraction method for the modeling of a current dipole in eeg source reconstruction using finite element head models,”SIAM Journal on Scientific Computing, vol. 30, no. 1, pp. 24–45, 2008.DOI: 10.1137/060659053
-
[27]
Kaipio and E
J. Kaipio and E. Somersalo,Statistical and computational inverse problems. Springer Science & Business Media, 2006, vol. 160
2006
-
[28]
Unbiased nonlinear average age-appropriate brain templates from birth to adulthood,
V . Fonov, A. Evans, R. McKinstry, C. Almli, and D. Collins, “Unbiased nonlinear average age-appropriate brain templates from birth to adulthood,”NeuroImage, vol. 47, S102, 2009, Organization for Human Brain Mapping 2009 Annual Meeting,ISSN: 1053-8119.DOI: 10.1016/S1053-8119(09)70884-5
-
[29]
Justification of point electrode models in electrical impedance tomog- raphy,
M. Hanke, B. Harrach, and N. Hyvönen, “Justification of point electrode models in electrical impedance tomog- raphy,”Mathematical Models and Methods in Applied Sciences, vol. 21, no. 06, pp. 1395–1413, 2011.DOI: 10.1142/S0218202511005362
-
[30]
Advanced boundary electrode modeling for tes and parallel tes/eeg,
S. Pursiainen, B. Agsten, S. Wagner, and C. H. Wolters, “Advanced boundary electrode modeling for tes and parallel tes/eeg,”IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 26, no. 1, pp. 37–44, 2018.DOI: 10.1109/TNSRE.2017.2748930
-
[31]
Efficient computation of lead field bases and influence matrix for the fem-based eeg and meg inverse problem,
C. Wolters, L. Grasedyck, and W. Hackbusch, “Efficient computation of lead field bases and influence matrix for the fem-based eeg and meg inverse problem,”Inverse Problems, vol. 20, pp. 1099–1116, Aug. 2004.DOI: 10. 1088/0266-5611/20/4/007
2004
-
[32]
Weighted sparsity regular- ization for solving the inverse eeg problem: A case study,
O. L. Elvetun and N. Sudheer, “Weighted sparsity regular- ization for solving the inverse eeg problem: A case study,” Biomedical Signal Processing and Control, vol. 107, p. 107 673, 2025,ISSN: 1746-8094.DOI: 10.1016/j.bspc. 2025.107673
-
[33]
Selvaraju, Michael Cogswell, Ab- hishek Das, Ramakrishna Vedantam, Devi Parikh, and Dhruv Batra
Y . Rubner, C. Tomasi, and L. Guibas, “A metric for dis- tributions with applications to image databases,” inSixth International Conference on Computer Vision (IEEE Cat. No.98CH36271), 1998, pp. 59–66.DOI: 10.1109/ICCV. 1998.710701
-
[34]
Improving source reconstructions by combining bioelectric and biomagnetic data,
M. Fuchs et al., “Improving source reconstructions by combining bioelectric and biomagnetic data,”Electroen- cephalography and Clinical Neurophysiology, vol. 107, no. 2, pp. 93–111, 1998,ISSN: 0013-4694.DOI: 10.1016/ S0013-4694(98)00046-7
1998
-
[35]
F. Neugebauer et al., “Validating eeg, meg and combined meg and eeg beamforming for an estimation of the epileptogenic zone in focal cortical dysplasia,”Brain Sciences, vol. 12, no. 1, 2022,ISSN: 2076-3425.DOI: 10.3390/brainsci12010114
-
[36]
In silico study of the effects of cerebral circulation on source localization using a dynamical anatomical atlas of the human head,
J. Lahtinen, F. Moura, M. Samavaki, S. Siltanen, and S. Pursiainen, “In silico study of the effects of cerebral circulation on source localization using a dynamical anatomical atlas of the human head,” eng,Journal of neural engineering, vol. 20, no. 2, pp. 26 005–, 2023, ISSN: 1741-2560
2023
-
[37]
Malmivuo and R
J. Malmivuo and R. Plonsey,Bioelectromagnetism: Principles and Applications of Bioelectric and Bio- magnetic Fields. Oxford University Press, Oct. 1995, ISBN: 9780195058239.DOI: 10 . 1093 / acprof : oso / 9780195058239.001.0001
1995
-
[38]
Bioelectro- magnetism in human brain research: New applications, new questions.,
J. Gross, M. Junghöfer, and C. Wolters, “Bioelectro- magnetism in human brain research: New applications, new questions.,”The Neuroscientist : a review journal bringing neurobiology, neurology and psychiatry, vol. 29, no. 1, pp. 62–77, 2023,ISSN: 1089-4098 (Electronic). DOI: 10.1177/10738584211054742
-
[39]
Sparsity enables estimation of both subcortical and cortical activity from meg and eeg,
P. Krishnaswamy et al., “Sparsity enables estimation of both subcortical and cortical activity from meg and eeg,” eng,PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMER- ICA, vol. 114, no. 48, E10465–E10474, 2017,ISSN: 0027-8424
2017
-
[40]
Localization of deep brain activity with scalp and subdural eeg,
M. Fahimi Hnazaee et al., “Localization of deep brain activity with scalp and subdural eeg,” eng,NeuroImage (Orlando, Fla.), vol. 223, pp. 117 344–, 2020,ISSN: 1053-8119
2020
-
[41]
Reconstructing subcortical and corti- cal somatosensory activity via the ramus inverse source analysis technique using median nerve sep data,
A. Rezaei et al., “Reconstructing subcortical and corti- cal somatosensory activity via the ramus inverse source analysis technique using median nerve sep data,” eng, NeuroImage (Orlando, Fla.), vol. 245, pp. 118 726–, 2021,ISSN: 1053-8119
2021
-
[42]
Localization of realistic spatial patches of complex source activity in meg and eeg,
A. Giri, L. Hecker, J. C. Mosher, A. Adler, and D. Pan- tazis, “Localization of realistic spatial patches of complex source activity in meg and eeg,”IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 33, pp. 4241–4254, 2025.DOI: 10 . 1109 / TNSRE . 2025 . 3622587
2025
-
[43]
Forward- Inverse interplay in FEM-Based EEG source imaging: Distributional signatures of advanced source models and inverse solvers,
S. Söderholm, J. Lahtinen, and S. Pursiainen, “Forward- Inverse interplay in FEM-Based EEG source imaging: Distributional signatures of advanced source models and inverse solvers,” in2026 IEEE International Conference on Metrology for eXtended Reality, Artificial Intelligence and Neural Engineering (MetroXRAINE) (IEEE MetroX- RAINE 2026), Chemnitz, German...
2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.