J. Fluid Mech. (2013), vol. 721, pp. 199–224. c Cambridge University Press 2013 199 doi:10.1017/jfm.2013.73 Thermal boundary layer structure in turbulent Rayleigh–Bénard convection in a rectangular cell Quan Zhou 1,2 and Ke-Qing Xia 2, 1 Shanghai Institute of Applied Mathematics and Mechanics, Shanghai University, Shanghai 200072, China 2 Department of Physics, The Chinese University of Hong Kong, Shatin, Hong Kong, China (Received 3 October 2012; revised 24 January 2013; accepted 29 January 2013) We report high-spatial-resolution measurements of the thermal boundary layer (BL) properties in turbulent thermal convection. The experiment was made near the lower conducting plate of a water-filled rectangular convection cell of height 0.76 m, with a Prandtl number Pr = 4.3 and over the Rayleigh-number range 2 × 10 10 < Ra < 7 × 10 11 . Time series of the local temperature at various vertical distance z from the plate were measured. Statistical properties of the profiles of the temperature, i.e. the mean temperature T , fluctuating temperature root mean square (r.m.s.) σ T , temperature skewness S T , and flatness F T , and those of the temperature time derivative, i.e. the r.m.s. σ T , skewness S T and flatness F T of the derivative, are studied. It is found that most of these quantities exhibit some degree of invariability with Ra, especially for the regime inside the thermal BL. When comparing with the mean temperature profiles, the profiles of the second moment of temperature seem to possess a higher level of universality. It is shown that the distance δ σ from the plate to the maximal temperature r.m.s. position provides a natural length scale for the characterization of the thermal BL, as the statistical properties of the temperature field, such as its r.m.s., skewness and flatness, are all sharply different below and above this length scale, i.e. below δ σ , σ T increases linearly with the vertical distance z from the plate and S T is close to zero and F T is close to three and both quantities remains nearly constant, whereas above δ σ the decay of σ T obeys a logarithmic behaviour and S T and F T both exhibit a hill-like structure. It is also found that near the plate T , σ T and σ T all increase linearly with z. Our observations further reveal that such linear dependence occurs within a self-similar region of the thermal BL, where the temperature probability density functions can be scaled onto a single distribution that differs slightly from the Gaussian distribution. The Ra-dependencies of various thermal BL properties are also studied and our results yield δ th /H = (6.85 ± 0.70)Ra 0.33±0.03 , δ σ /H = (2.86 ± 0.30)Ra 0.31±0.03 and δ σ /H = (25 ± 3)Ra 0.38±0.05 , where H is the height of the cell, δ th and δ σ are the BL thicknesses determined respectively from the profiles of T and σ T . Key words: enard convection, boundary layer structure, turbulent convection † Email address for correspondence: kxia@phy.cuhk.edu.hk