A three-dimensional finite-difference time-domain method that employs pure scattered-field formulation and perfectly matched layers (PML) as the absorbing boundary condition is developed for solving subsurface-scattering. ...