“…The solution is obtained in an iterative fashion directly in the time domain because FDTD is effectively a wave propagator. GprMax3D is part of a suite of GPR modelling algorithms under the name GprMax which has been used for a wide range of GPR simulations by a number of researchers (Giannopoulos and Diamanti, 2008;Jeannin et al, 2006;Persico and Soldovieri, 2008;Saintenoy et al, 2008;Tsoflias and Becker, 2008;Wilson et al, 2009). Some of its features included are: Perfectly Matched Layer (PML) boundary conditions to efficiently truncate the computational domain (Gedney, 1996), modelling of frequency dependent materials, interface roughness, heterogeneous media properties and fine features (Diamanti and Giannopoulos, 2009;Giannopoulos, 2005Giannopoulos, , 2008Giannopoulos and Diamanti, 2003, 2005.…”