<p>An important influence that distinguishes the dynamic response of tunnels embedded in saturated frozen soil from that of unfrozen soil is the freezing of liquid-phase water in the pores, and the formation of ice alters the propagation of loads through the frozen soil. The dynamic response of the tunnel embedded in saturated frozen half-space when loads are applied to the lining invert is addressed by the wave function method. The saturated frozen soil is regarded as a poroelastic medium, the lining is regarded as a hollow cylinder, and the total wave field in the surrounding medium consists of outgoing cylindrical waves and down-going plane waves. The plane wave function and cylindrical wave function could be converted to obtain the boundary conditions that are convenient for solving the ground surface and soil-lining interface in the rectangular and cylindrical coordinate systems. The parametric analysis demonstrates that the dynamic response of a tunnel embedded in a saturated frozen half-space and the surrounding medium could be influenced by the permeability, porosity, ice content of frozen soil and the change of the tunnel depth. This method provides a novel approach to the safety and stability of tunnels in cold regions and serves as a reference for predicting the vibration of tunnels embedded in saturated frozen half-spaces.</p>