Improved methods are formulated for solution of the linearized Poisson-Boltzmann equation, to be used in conjunction with electronic structure calculation on a solute together with dielectric continuum representation of the salt-containing solvent. Volume polarization effects due to quantum mechanical penetration of solute charge density outside the cavity that excludes solvent are treated by exact and by approximate methods analogous to those previously developed for the salt-free case. With boundary element approaches, exact solutions lead to coupled equations for a pair of cavity surface distributions that mimic the polarization of the solvent dielectric and the ionic atmosphere. A novel means is found to effectively decouple these equations, yielding more efficient practical methods for their numerical solution. Detailed comparisons are given to related boundary element formulations previously reported in the literature, which neglect volume polarization, and analogous decoupling is also found for the pair of surface distributions invoked there. Illustrative results are provided for a simple spherical example. (C) 2004 American Institute of Physics.