Yes / A set of shallow-water equations (SWEs) based on a k-epsilon Reynold stress model is established to simulate the turbulent flows over a complex roughness bed. The fundamental equations are discretized by the second-order finite-difference method (FDM), in which spatial and temporal discretization are conducted by staggered-grid and leap-frog schemes, respectively. The turbulent model in this study stems from the standard k-epsilon model, but is enhanced by replacing the conventional vertical production with a more rigorous and precise generation derived from the energy spectrum and turbulence scales. To verify its effectiveness, the model is applied to compute the turbulence in complex flow surroundings (including a rough bed) in an abrupt bend and in a natural waterway. The comparison of the model results against experimental data and other numerical results shows the robustness and accuracy of the present model in describing hydrodynamic characteristics, especially turbulence features on the complex roughness bottom. / National Key Research and Development Program of China (Grant No: 2016YFE0122500, 2013CB036401 and 2013CB036402), China Postdoctoral Science Foundation (Grant No: 2016M591184) and Programme of Introducing Talents of Discipline to Universities (Grant No: BC2018038) / Research Development Fund Publication Prize Award winner, June 2019.