This paper investigates the integral properties of the mean temperature dissipation a(00/0y)2 and the temperature variance production P-theta = R-v theta(partial derivative Theta/partial derivative y) in fully developed turbulent channel flow, and in particular the effects of Prandtl number Pr = v/alpha and Reynolds number Re-tau = delta u(tau)/v (and/or Peclet number Pe(tau) = Re tau Pr. Here Theta is the mean temperature and R-v theta is the wall-normal turbulent heat flux. y is the wall normal direction and delta is the channel half-height. v and alpha are the fluid viscosity and thermal diffusivity, respectively. ur is the friction velocity. The main findings of the present work include: (1) Identities are derived from the mean thermal energy equation for the global integral of the mean temperature dissipation and temperature variance production for four heating conditions: constant heat source (CHS), constant wall heat flux (CHF), constant temperature difference (CTD), and constant wall temperature (CTD). Identities for the CHS and CHF conditions are consistent with those derived by Pirozzoli et al. (2016) and Abe and Antonia (2017). The identities are verified using publicly available direct numerical simulation (DNS) data under three heating conditions CHS, CHF, and CTD. (2) The global integral of the mean temperature dissipation is found to depend strongly on the fluid Prandtl number, but it is largely independent of the Reynolds number of the flow. In inner-scaling, this global integral can be approximated as(1/Pr)integral(delta+)(0) (partial derivative Theta(+)/partial derivative y(+))(2)dy(+) approximate to 9.13Pr(1-n) where Theta(+) = Theta/theta(tau) is the inner-scaled mean temperature and theta(tau) is the friction temperature. y(+) = yu(tau)/v is the inner-scaled wall-normal distance. n approximate to 1/2 for Pr <= 1 and n approximate to 1/3 for Pr >= 1. The strong dependence of (1/Pr)integral(delta+)(0) (partial derivative Theta(+)/partial derivative y(+))(2)dy(+) on the Prandtl number is directly rooted in the molecular thermal diffusion sub-layer. (3) The global integral of temperature variance production is found to strongly depend on both Prandtl number and Reynolds number. (3.1) At sufficiently high Reynolds number and Peclet number, the integral of temperature variance production is found to exhibit a logarithmic-like layer similar to that of the mean temperature profile. (3.2) The effect of Prandtl number is mainly caused by the molecular thermal diffusion sub-layer and buffer layer. (3.3) Between the peak wall-normal turbulent heat flux R-v theta location, y(m tau), and the channel centerline, delta, the integral of temperature variance production increases with Peclet number in a logarithmic-like fashion. (3.4) At low Prandtl number Pr <= 1, the integral of temperature variance production is found to be equally partitioned around y(mi). (C) 2018 Elsevier Ltd. All rights reserved.