This paper aims at investigating the effects of variations in thrust hydrodynamic bearing (HDB) parameters such as axial stiffness and damping coefficients on the axial vibration of disk-spindle systems in hard disk drives. For a parametric study, a closed-form axial frequency response function (FRF) of HDB spindle systems is derived as a function of the axial stiffness and damping coefficients of thrust HDBs. It is known that the axial vibration of the disk-spindle system is composed of two main parts: the vibration of the rigid hub in the axial direction and the disk deflection in the transverse direction. The results from this research clearly show that the vibration amplitudes at low frequency range is dominated by the axial vibration of the hub, and the amplitude of the unbalanced (0,0) mode is dominated by the disk deflection. The parametric study reveals that at low frequency range an increase in the bearing stiffness significantly reduces the hub axial vibration, and hence the axial vibration of the disk-spindle system. Surprisingly, a too much increase in the damping results in a higher amplitude of the unbalanced (0,0) mode. This is because a heavy damping constrains the hub vibration to nearly no motion, resulting in a direct transmission of vibration from the base to disk. To confirm the parametric study, a vibration test was performed on two HDB spindle motors with identical design but different fluid viscosity. The higher viscosity represents the higher axial stiffness and damping in the thrust bearing. The test result indicates that the spindle motor with higher viscosity has a larger unbalanced (0,0) mode amplitude when subjected to an axial base excitation.