Radiation transport plays an important role in stellar atmospheres, but the effects of turbulence are being obscured by other effects such as stratification. Using radiative hydrodynamic simulations of forced turbulence, we determine the decay rates of sinusoidal large-scale temperature perturbations of different wavenumbers in the optically thick and thin regimes. Increasing the wavenumber increases the rate of decay in both regimes, but this effect is much weaker than for the usual turbulent diffusion of passive scalars, where the increase is quadratic for small wavenumbers. The turbulent decay is well described by an enhanced Newtonian cooling process in the optically thin limit, which is found to show a weak increase proportional to the square root of the wavenumber. In the optically thick limit, the increase in turbulent decay is somewhat steeper for wavenumbers below the energy-carrying wavenumber of the turbulence, but levels off toward larger wavenumbers. In the presence of turbulence, the typical cooling time is comparable to the turbulent turnover time. We observe that the temperature takes a long time to reach equilibrium in both the optically thin and thick cases, but in the former, the temperature retains smaller scale structures for longer.