Abstract. The Shallow Subsurface Radar (SHARAD), a sounding radar equipped on the Mars Reconnaissance Orbiter (MRO), has produced highly valuable information about the Martian subsurface. In particular, the complicated substructures of Mars such as polar deposit, pedestal crater and the other geomorphic features involving possible subsurface ice body has been successfully investigated by SHARAD. In this study, we established a 3D subsurface mapping strategy employing the multiple SHARAD profiles. A number of interpretation components of SHARAD signals were integrated into a subsurface mapping scheme using radargram information and topographic data, then applied over a few mid latitude Lobate Debris Aprons (LDAs). From the identified subsurface layers of LDA, and the GIS data base incorporating the other interpretation outcomes, we are expecting to trace the origin of LDAs. Also, the subsurface mapping scheme developed in this study will be further applied to other interesting Martian geological features such as inter crater structures, aeolian deposits and fluvial sediments. To achieve higher precision sub-surface mapping, the clutter simulation employing the high resolution topographic data and the upgraded clustering algorithms assuming multiple sub-surface layers will be also developed.