This article describes a plane-wave time-domain (PWTD) algorithm that facil
itates the fast evaluation of transient wave fields produced by surface sca
ttering. The algorithm presented relies on a Whittaker-type expansion of tr
ansient fields in terms of propagating plane waves. The incorporation of th
e PWTD scheme into existing marching-on-in-time- (MOT-) based integral-equa
tion solvers is elucidated. It is shown that the computational cost of perf
orming a surface-scattering analysis, using two-level and multi-level PWTD-
enhanced MOT schemes, scales as O(NtNs1.5 logN(s)) and O(NtNs log(2) N-s),r
espectively, when the surface source density is represented by N-s spatial
and N-t temporal samples. Hence, the computational cost of the proposed alg
orithms scales much more favorably than that of classical MOT schemes, whic
h scale as O(NtNs2). Therefore, PWTD-enhanced MOT schemes make possible the
analysis of broadband scattering from structures of unprecedented dimensio
ns.