Propagation of a surface discharge on water is evaluated using simulation results of the voltage/current waveforms and discharge contact area. Voltage/current waveforms are calculated using an exponential function which assumes the resistance of water decreases with increasing discharge contact area. A conductive disk having a given potential is used to model the discharge contact area and the temporal variation of the radius of the disk is discussed as the discharge propagation. The calculation in the current field is replaced by a calculation of the electrostatic field using the similarity between current and electrostatic fields. The calculation of the electrostatic field is conducted by a charge simulation method. The discharge phenomena are classified into two stages which are the breakdown of needle-to-water gap and the surface discharge propagation on water. The electrostatic field calculation at needle-to-water gap is performed to determine the initial discharge contact area. Expansion of the contact area in the creepage direction, is evaluated by the electric field calculation at the edge of the conductive disk. The expansion terminates when the field becomes lower than 26 kV/cm. The velocity of the expansion increases with the applied voltage. The maximum contact radius increases with decreasing conductivity of water under the same applied voltage. The tendency of simulation results of the expansion of the contact area are consistent with the previous observation results of the discharge propagation.