“…There are solutions based on approximations such as the reflection-coefficient method [11,12], which is computationally faster, but as it assumes that the waves incident on the ground are plane waves, it presents losses of accuracy when the approximation is not valid [13]. Meanwhile, these problems can be solved by the finite-difference time-domain (FDTD) method [14][15][16], time-domain integral equation (TDIE) method [17][18][19], and discontinuous Galerkin time-domain (DGTD) method [20,21]. Most of these studies are carried out in free space or isotropic half-space because the reflected field of multilayered and even-layered anisotropic half-space is difficult to obtain.…”