To understand past Martian climates, it is important to know the distribution and nature of water ice on Mars. Impact craters are widely used ubiquitous indicators for the presence of subsurface water or ice on Mars. Remote sensing observations and numerical simulation are powerful tools for investigating morphological and topographic features on planetary surfaces, and we can use the morphology of layered ejecta craters and hydrocode modeling to constrain possible layering and impact environments. The approach of this work consists of three stages：Firstly, the morphological characteristics of the Martian layered ejecta craters are performed based on Martian images and DEM data. Secondly, numerical modeling layered ejecta are performed through the hydrocode iSALE (impact-SALE). We present hydrocode modeling of impacts onto targets with a single icy layer within an otherwise uniform basalt crust to quantify the effects of subsurface H₂O on observable layered ejecta morphologies. The model setup is based on a layered target made up of a regolithic layer (described by the basalt ANEOS), on top an ice layer (described by ANEOS equation of H₂O ice), in turn on top of an underlying basaltic crust. The bolide is a 0.8 km diameter basaltic asteroid hitting the Martian surface vertically at a velocity of 12.8 km/s. Finally, the numerical results are compared with the MOLA DEM profile in order to analyze the formation mechanism of Martian layered ejecta craters. Our simulations suggest that the presence of an icy layer significantly modifies the cratering mechanics, and many of the unusual features of SLE craters may be explained by the presence of icy layers. Impact cratering on icy satellites is significantly affected by the presence of subsurface H₂O.