Abstract. In geotechnical earthquake engineering, wave propagation plays a fundamental role in engineering applications related to the dynamic response of geotechnical structures and to site response analysis. However, current engineering practice is primarily concentrated on the investigation of shear wave propagation and the corresponding site response only to the horizontal components of the ground motion.Due to the repeated recent observations of strong vertical ground motions and compressional damage of engineering structures, there is an increasing need to carry out a comprehensive investigation of vertical site response and the associated compressional wave propagation, particularly when performing the seismic design for critical structures (e.g. nuclear power plants and high dams). Therefore, in this paper, the compressional wave propagation mechanism in saturated soils is investigated by employing hydro-mechanically (HM) coupled analytical and numerical methods. A HM analytical solution for compressional wave propagation is first studied based on Biot's theory, which shows the existence of two types of compressional waves (fast and slow waves) and indicates that their characteristics (i.e. wave dispersion and attenuation) are highly dependent on some key geotechnical and seismic parameters (i.e. the permeability, soil stiffness and loading frequency). The subsequent HM Finite Element (FE) study reproduces the duality of compressional waves and identifies the dominant permeability ranges for the existence of the two waves. In particular the existence of the slow compression wave is observed for a range of permeability and loading frequency that is relevant for geotechnical earthquake engineering applications. In order to account for the effects of soil permeability on compressional dynamic soil behaviour and soil properties (i.e. P-wave velocities and damping ratios), the coupled consolidation analysis is therefore recommended as the only tool capable of accurately simulating the dynamic response of geotechnical structures to vertical ground motion at intermediate transient states between undrained and drained conditions. 2