In ice gouging analyses, the seabed is often simplified as a uniform material domain, over-looking the potential complexities inherent in layered seabed formations, which are wide-spread in numerous Arctic regions. This study investigates the pipeline response into the distinct seabed configurations, including layered soft over stiff clay, layered stiff over soft clay, and uniformly soft and uniformly stiff compositions, to ice gouging. This was achieved through comprehensive large deformation finite element (LDFE) analysis, employing a Coupled Eulerian Lagrangian (CEL) algorithm. Incorporation of the strain-rate dependency and strain-softening effects involved the development of a user-defined subroutine and incremental update of the undrained shear strength within the Abaqus software. The pipeline is explicitly incorporated into the model, allowing for a detailed investigation of its response to ice gouging phenomena within these layered seabed scenarios. This research demonstrates that simplifying a layered seabed into a uniform can be misleading. Such an approach can lead to significant errors, potentially overestimating or underestimating pipeline response to ice gouging. By incorporating these complexities, engineers can develop more optimized designs with enhanced safety margins in Arctic environments prone to ice gouging phenomena.