The conventional transfer matrix method, adopted formerly in a semi analytical solution for layered elastic or viscoelastic asphalt pavement system, has inherent deficiencies, such as illcondition matrix, numerical overflow and error accumulation, due to exponential items. This phenomenon is more evident in multilayered dynamic analysis with imperfect interfaces. Moreover, several studies revealed that pavement materials exhibit transverse isotropy in service. Consequently, a novel semi analytical solution methodology, wave propagation approach, was proposed herein to calculate the dynamic responses of asphalt pavement under an impact load considering the transversely isotropic material model and the imperfect interfaces. First, the transfer matrix was established based on matrix theory and wave propagation approach, while the relation between the state vector and wave vector in the transformed domain was constructed simultaneously. Then, combined with the boundary conditions and interface contact conditions, the solution of the wave