Abstract: In this paper, a parallel time domain electric field integral equation is applied to numerical simulation of two-dimensional electromagnetic scattering from time-evolving ocean surfaces. The proper choice of temporal basis functions and precise evaluations of matrix elements are used to overcome the late-time instability. Quadratic B-spline function is employed as the temporal basis function, while triangle function is employed as the spatial basis function. The matrix elements are precisely evaluated with analytical formulae for time integrals terms and singular terms. The tapered modulated Gaussian pulse is used to reduce the edge effects resulted from cutting down the infinite sea surfaces. A scheme of message passing interface and compressed storage of sparse matrix are adopted to improve the computational efficiency and reduce the computational storage requirement. Numerical results demonstrate that the method proposed is accurate and efficient for analyzing transient scattering from time-evolving ocean surfaces.