Finite-difference time domain (FDTD) techniques are widely used to model the propagation of viscoelastic waves through complex and heterogeneous structures. However, in the specific case of media mixing liquid and solid, attempts to model continuous media onto a Cartesian grid produces errors when the liquid-solid interface between different media do not align precisely with the Cartesian grid. The increase in spatial resolution required to eliminate this grid staircasing effect can be computationally prohibitive. Here, a modification to the Virieux staggered-grid FDTD scheme called the superposition method is presented. This method is intended to reduce this staircasing effect while keeping a manageable computational time. The method was validated by comparing low-spatial-resolution simulations against simulations with sufficiently high resolution to provide reasonably accurate results at any incident angle. The comparison of the root-mean-square of the stress amplitude maps showed that the amplitude of artifactual waves could be reduced by several orders of magnitude when compared to the Virieux staggered-grid FDTD method and that the superposition method helped to significantly reduce the staircasing effect in FDTD simulations.