<p>The local contribution of clouds to the surface energy balance and temperature variability is an important topic in order to apprehend how this intake affects local climate variability and extreme events, how this contribution varies from one place to another, and how it evolves in a warming climate. The scope of this study is to understand how clouds impact temperature variability, to quantify their contribution, and to compare their effects to other surface processes. To do so, we develop a method to estimate the different terms that control temperature variability at the surface (&#8706;T<sub>2m</sub> /&#8706;t) by using this equation: <strong>&#8706;T<sub>2m</sub> /&#8706;t=R+HA+HG+Adv</strong> where R is the radiation that is separated into the cloud term (R<sub>cloud</sub>) and the clear sky one (R<sub>CS</sub>), HA the atmospheric heat exchange, HG the ground heat exchange, and Adv the advection. These terms are estimated hourly, almost only using direct measurements from SIRTA-ReOBS dataset (an hourly long-term multi-variables dataset retrieved from SIRTA, an observatory located in a semi-urban area 20-km South-West of Paris; Chiriaco et al., 2019) for a five-years period. The method gives good results for the hourly temperature variability, with a 0.8 correlation coefficient and a weak residual term between left part (directly measured) and right part of the equation.</p><p>A bagged decision trees analysis of this equation shows that R<sub>CS</sub> dominates temperature variability during daytime and is mainly modulated by cloud radiative effect (R<sub>cloud</sub>). During nighttime, the bagged decision trees analysis determines that R<sub>cloud</sub> is the term controlling temperature changes. When a diurnal cycle analysis (split into seasons) is performed for each term, HA becomes an important negative modulator in the late afternoon, chiefly in spring and summer, when evaporation and thermal conduction are increased. In contrast, HG and Adv terms do not play an essential role on temperature variability at this temporal scale and their contribution is barely considerable in the one-hour variability, but still they remain necessary in order to obtain the best coefficient estimator between the directly measured observations and the method estimated. All terms except advection have a marked monthly-hourly cycle.</p><p>Next steps consist in characterize the types of clouds and study their physical properties corresponding to the cases where R<sub>cloud</sub> is significant, using the Lidar profiles also available in the SIRTA-ReOBS dataset.</p>