AbstractA numerical method is given for calculating resistivity and induced polarization anomalies produced by a surface polarization model. Surface polarization is generated when a purely electronic conductor is located in an electrolyte environment. The system that develops on the boundary between the conductor and the electrolyte is described macroscopically by a net surface charge distribution and an electric double layer. An integral equation is derived for the potential by assuming that the electronic conductor forms an equipotential system and that the polarization impedance across the boundary is linear. The integral equation is solved by means of the method of subsections. As an application some numerical modeling results are presented. The surface impedance values used in calculations are based on laboratory measurements that are briefly described. Implications of the results for scale modeling are discussed.