<p><strong>Abstract.</strong> The role of thermally-driven local downslope or katabatic flows in the dynamics and turbulent features of the stable boundary layer (SBL) is investigated using observations. Measurements are carried out in a relatively flat area 2-km away from the steep slopes of the Guadarrama Mountain Range (Spain). Forty katabatic events are selected from an observational database spanning the 2017-summer period, by using an objective and systematic algorithm that is able to account for local and synoptic forcings. We subsequently classify the katabatic events into weak, moderate and intense according to the observed maximum wind speed. This classification enables us to contrast the main differences in dynamics and thermal structure. We find that the stronger katabatic events are associated with an earlier onset time of these flows. We relate it to very low soil-moisture values (<&#8201;0.07&#8201;m<sup>3</sup>&#8201;m<sup>&#8722;3</sup>, i.e. smaller than the median during the analysed period) and a weak synoptic wind (<i>V</i><sub>850</sub>&#8201;<&#8201;6&#8201;m&#8201;s<sup>&#8722;1</sup>) having the same direction as the katabatic. The relative flatness of the area favours the formation of very stable boundary layers characterized by a longwave radiative cooling of around 60&#8211;70&#8201;W&#8201;m<sup>&#8722;2</sup> and very weak turbulence (friction velocity (<i>u</I><sub>*</sub>)&#8201;<&#8201;0.1&#8201;m&#8201;s<sup>&#8722;1</sup>). They occur when katabatics are weak, and are occasionally associated with the formation of skin flows, that are manifested as weak jets (<i>U</i>&#8201;<&#8201;1&#8201;m&#8201;s<sup>&#8722;1</sup>) at 3&#8201;m. Intense katabatics, instead, are characterised by a strong and increasing bulk shear (the maximum <i>u</i><sub>*</sub> is close to 1&#8201;m&#8201;s<sup>&#8722;1</sup>) that avoids the development of the surface-based thermal inversion, giving rise to the so-called weakly stable boundary layer. We identify the transition between the two regimes for a threshold katabatic wind speed of around 1.5&#8201;m&#8201;s<sup>&#8722;1</sup>, in agreement with the hockey-stick transition hypothesis. Our analysis is extended by calculating non-dimensional numbers to characterize the transition: the shear capacity, the bulk Richardson number and <i>z/L</i>. On the other hand, by inspecting individual weak and intense events, we further explore the interaction between katabatic flows and turbulence, and the impact on CO<sub>2</sub> concentration. By relating the dynamics of the two regimes with the CO<sub>2</sub> budget, we are able to estimate the contribution of the different terms. For the intense event, indeed, we infer a horizontal transport of 67&#8201;ppm in 3&#8201;h driven by the katabatic advection.</p>