Ice cover in an open channel can influence the flow structure, such as the flow velocity, Reynolds stress and turbulence intensity. This study analyzes the vertical distributions of velocity, Reynolds stress and turbulence intensity in fully and partially ice-covered channels by theoretical methods and laboratory experiments. According to the experimental data, the vertical profile of longitudinal velocities follows an approximately symmetry form. Different from the open channel flow, the maximum value of longitudinal velocity occurs near the middle of the water depth, which is close to the channel bed with a smoother boundary roughness compared to the ice cover. The measured Reynolds stress has a linear distribution along the vertical axis, and the vertical distribution of measured turbulence intensity follows an exponential law. Theoretically, a two-power-law function is presented to obtain the analytical formula of the longitudinal velocity. In addition, the vertical profile of Reynolds stress is obtained by the simplified momentum equation and the vertical profile of turbulence intensity is investigated by an improved exponential model. The predicted data from the analytical models agree well with the experimental ones, thereby confirming that the analytical models are feasible to predict the vertical distribution of velocity, Reynolds stress and turbulence intensity in ice-covered channels. The proposed models can offer an important theoretical reference for future study about the sediment transport and contaminant dispersion in ice-covered channels.