Direct numerical simulations (DNS) of turbulent thermal convection in a Pr = 0.7 fluid up to Ra = 10 12 are used to study the statistics of thermal plumes. At various vertical locations in a cylindrical set-up with aspect ratio Γ = width/height = 1/3, plumes are identified and their properties extracted. It is found that plumes are much less likely to be emitted from plate regions with large wind shear. Close to the plates, the plumes have a unimodal log-normal distribution, whereas at more central locations the distribution becomes weakly bimodal, which can be traced back to clustering of the plumes and influence of the large-scale circulation. The number of hot plumes decreases with height. The width of the plumes scales with Ra approximately as Nu −1 , indicating that it is determined by the thermal boundary layer thickness.