Subsurface infrastructure at shallow depth is becoming more common due to increase in urban settlements and space constraints. Among these infrastructural elements, train tunnels are essential for urban transportation networks. The stability of such tunnels is even more challenging in seismically active regions due to the potential impact of seismic forces. Past studies have shown that circular tunnels built in the soils that are prone to liquefaction experience uplift and deformations due to increase in the generation of excess pore water pressures and ground settlement. However, few studies have focused on understanding the behaviour of tunnels in liquefiable saturated sands during seismic loading. In this study, the advanced constitutive model PM4Sand is adopted to study the seismic response of shallow circular tunnel embedded in saturated sand, using numerical modelling in PLAXIS 2D. The PM4Sand constitutive model is able to capture the generation of excess pore water pressure during liquefaction and post-liquefaction phenomena, as well as soil dilatancy behaviour. The numerical model is first validated with the results of centrifuge model tests available in existing literature. The study aims to evaluate the parameters for the PM4Sand model through a method of back calculation by utilizing the data presented in the study of Chian et al., (2014) and subsequently validate these obtained values through a cross-validation procedure. The study highlights the complex interaction between the tunnel structures, seismic forces, and soil response in terms of acceleration response, magnitude of excess pore pressure generation and uplift displacement.