This　paper　examines　non-homogeneous　geo-material　halfspaces　under　the　axisymmetric　indentation　of　a　circular　rigid　plate.　A　multilayered　model　is　proposed　to　approximate　arbitrary　variations　of　the　elastic　parameters　in　geo-materials　with　depth.　The　non-homogeneous　geo-material　halfspaces　can　have　variable　shear　modulus　and　Poisson＇s　ratio　with　depth.　Numerical　verification　is　performed　for　the　contact　problem　of　a　non-homogeneous　halfspace　with　variable　shear　modulus　and　constant　Poisson＇s　ratio.　Results　illustrate　the　effects　of　material　inhomogeneity　on　the　elastic　fields　of　a　non-homogeneous　halfspace.　©　2024　Elsevier　Ltd