We present an extension of the Polarizable Continuum Model (PCM) to the calculation of solvent effects on indirect spin-spin coupling constants for Hartree-Fock wave functions and Density Functional Theory. This is achieved by implementing the PCM model for singlet and triplet linear response functions. The new code is used for calculating the solvent effects on the indirect spin-spin coupling constants of benzene. For the (1)J((HC)-C-13) coupling constants, our calculated solvent shifts are in good agreement with experimental observations when geometry relaxation is taken into account. However, our results do not support the extrapolated gas-phase value for this coupling constant. A new experimentally derived (1)J((HC)-C-13) for a vibrating benzene molecule at 300 K is proposed.