A 2-D subcrustal velocity model for the central Indian continental lithosphere has been derived by travel time and relative amplitude modeling of a digitally normalized analog seismic record section of the Hirapur-Mandla DSS profile, using a ray-tracing technique. Some prominent wave groups with apparent velocities slightly higher than the Moho reflection phase (P M P) are identified on the normalized record sections assembled with a reduction velocity of 6 km s )1 . We interpret these phases as the wide-angle reflections from subcrustal lithospheric boundaries. Comparison of synthetic seismograms with the observed record section shows that the observed phases cannot be explained either by multiples or by the P-to-S converted phase (P M S) from the Moho. Subcrustal velocity models either with a velocity increase or with a single low velocity layer (LVL) also do not provide a satisfactory fit. We infer that a subcrustal velocity model with two alternate LVLs (velocity 7.2 km s )1 ), separated by a 6-km thick high velocity layer (velocity 8.1 km s )1 ), can satisfy both the observed travel times and amplitudes. The prominent reflection phases are modeled at depths of 49, 51, 57 and 60 km. It is inferred that the subcrustal lithosphere in the central Indian region has a lamellar structure with varying structural and mechanical properties. The alternating LVLs, occurring at relatively shallow depths below Moho, may be associated with the zones of weakness and lower viscosity suggesting continued mobility, with a possible thermal source in the upper mantle. This explains the source of observed high heat flow values in the central Indian region.