Bioelectric source analysis in the human brain from scalp electroencephalography (EEG) signals is sensitive to the conductivities of different head tissues. The conductivity of tissues is subject dependent, so non-invasive methods for conductivity estimation are necessary to fine tune EEG models. To do so, the EEG forward problem solution (so-called lead field matrix) must be computed for a large number of conductivity configurations.Computing a lead field requires a matrix inversion which is computationally intensive for realistic head models. Thus, the required time for computing a large number of lead fields can become impractical. In this work, we propose to approximate the lead field matrix for a set of conductivity configurations, using the exact solution only for a small set of support points in the conductivity space. Our approach accelerates the computation time, while controlling the approximation error.Our method is tested on simulated and measured EEG data for brain and skull conductivity estimation. This test demonstrates that the approximation does not introduce any bias and runs significantly faster than if exact lead field were to be computed.Index Terms-EEG forward problem, EEG inverse problem, conductivity estimation, lead field matrix approximation.K. Maksymenko is with Athena,