We present a detailed simulational study of the vertical propagation of weak and strong impulses in deep gravitationally compacted granular columns [see R. S. Sinkovits and S. Sen, Phys. Rev. Lett. 74, 2686 (1995)]. The intergrain potential is assumed to be [Formula Presented], [Formula Presented], where [Formula Presented] is the overlap between the grains. Due to gravitational compaction, the magnitude of the overlap between the grains increases progressively with increasing depth. Therefore the sound velocity increases as an impulse travels vertically downward into a granular column. For weak impulses, our large scale simulational studies show that the sound velocity [Formula Presented], where [Formula Presented] is the depth at which [Formula Presented] is measured. This result, which has been obtained from particle dynamical studies, is in perfect agreement with the predictions based upon elasticity theory. We then extend our analysis to show that (i) for columns with small void fractions, [Formula Presented], [Formula Presented] and (ii) for large amplitude impulses, the velocity of the perturbation [Formula Presented] behaves very differently compared to [Formula Presented] at shallow depths with [Formula Presented] as [Formula Presented]. We also present a detailed numerical study of the velocity power spectra of the individual grains as a function of depth [Formula Presented]. We close with a discussion of the effects of both light and heavy impurities on the vertical sound and shock propagation. © 1996 The American Physical Society.