Langmuir circulations (LCs) generated by the interaction between wind-driven currents and surface waves can engulf the whole water column in neutrally stratified shallow water and interact with the turbulence in the bottom boundary layer. In this study, we perform a mechanistic study using wall-resolved large-eddy simulations (LES) based on the Craik-Leibovich equations to investigate the effects of LCs on turbulence statistics in the bottom half of shallow water. The highest. Reynolds number considered in this paper, Re-tau = 1000, is larger than the values considered in wall-resolved LES studies of shallow-water Langmuir turbulence reported in literature. The logarithmic layer is diagnosed based on a plateau region in the profile of a diagnostic function. It is found that the logarithmic layer disrupted at Re-tau = 395 reappears at Re-tau = 1000, but the von Karman constant is slightly different from the traditional value 0.41. To study the effects of LCs on turbulence statistics, LCs are extracted using streamwise averaging. The velocity fluctuations u'(i) are decomposed into a LC-coherent part u(i)(L) and a residual turbulence part u(i)(T). It is found that the profiles of LC-coherent. Reynolds shear stress -< u(L)v(L)> obtained at various Reynolds numbers are close to each other in the water-column coordinate y/h, with h being the half-water depth. As the Reynolds number (or, by definition, the ratio between the outer and inner length scales) increases, the influence of LCs on the near-bottom momentum transfer is reduced, which is responsible for the reappearance of the logarithmic layer. At all of the Reynolds numbers under investigation, the peaks of < u(L)u(L)> are collocated in the water-column coordinate y/h, while those of < u(T)u(T)> are collocated in the inner-scale coordinate y/(v/u(tau)). Due to the increase in the distance between the peaks of < u(L)u(L)> and < u(T)u(T)> with the Reynolds number, the profile of < u' u'> forms a bimodal shape at Re-tau 700 and 1000.