This study includes the results of a set of numerical simulations carried out for sands containing plastic/non-plastic fines, and silts with relative densities of approximately 30−40% under different surcharges on the shallow foundation using FLAC 2D. Each model was subjected to three ground motion events, obtained by scaling the amplitude of the El Centro (1940), Kobe (1995) and Kocaeli (1999) Q12earthquakes. Dynamic behaviour of loose deposits underlying shallow foundations is evaluated through fully coupled nonlinear effective stress dynamic analyses. Effects of nonlinear soil structure interaction (SSI) were also considered by using interface elements. This parametric study evaluates the effects of soil type, structure weight, liquefiable soil layer thickness, event parameters (e.g., moment magnitude of earthquake (Mw), peak ground acceleration PGA, PGV/PGA ratio and the duration of strong motion (D5−95) and their interactions on the seismic responses. Investigation on the effects of these parameters and their complex interactions can be a valuable tool to gain new insights for improved seismic design and construction.