The maximum velocity profile in turbulent wave boundary layer flows has been experimentally investigated under both regular and irregular wave conditions. Four types of seabed models are adopted, i.e., smooth, sand-covered, uniform-sphere-covered and nonuniform-stone covered. The results show that the maximum overshoot increases with the decreasing A/ks (A is the semi-excursion of fluid particles in the free stream and ks is the bottom roughness), but it is not notably influenced by the irregularity of the waves nor the nonuniformity of bottom roughness elements. Two-dimensional numerical simulations are carried out to reveal the physics behind the large overshoot behavior. It is found that the gap flow accelerates around roughness elements and subsequently contributes to the overshoot. A three-parameter defect function is proposed that well describes the maximum velocity profile, and all three model parameters are correlated with the governing flow parameter A/ks.