Abstract. Soil moisture is the main driver of temporal changes in values of the soil thermal conductivity. The latter is a key variable in land surface models (LSMs) used in hydrometeorology, for the simulation of the vertical profile of soil temperature in relation to soil moisture. Shortcomings in soil thermal conductivity models tend to limit the impact of improving the simulation of soil moisture in LSMs. Models of the thermal conductivity of soils are affected by uncertainties, especially in the representation of the impact of soil properties such as the volumetric fraction of quartz (q), soil organic matter, and gravels. As soil organic matter and gravels are often neglected in LSMs, the soil thermal conductivity models used in most LSMs represent the mineral fine earth, only. Moreover, there is no map of q and it is often assumed that this quantity is equal to the volumetric fraction of sand. In this study, q values are derived by reverse modelling from the continuous soil moisture and soil temperature sub-hourly observations of the Soil Moisture Observing System – Meteorological Automatic Network Integrated Application (SMOSMANIA) network at 21 grassland sites in southern France, from 2008 to 2015. The soil temperature observations are used to retrieve the soil thermal diffusivity (Dh) at a depth of 0.10 m in unfrozen conditions, solving the thermal diffusion equation. The soil moisture and Dh values are then used together with the measured soil properties to retrieve soil thermal conductivity (λ) values. For ten sites, the obtained λ value at saturation (λsat) cannot be retrieved or is lower than the value corresponding to a null value of q, probably in relation to a high density of grass roots at these sites or to the presence of stones. For the remaining eleven sites, q is negatively correlated with the volumetric fraction of solids other than sand. The impact of neglecting gravels and organic matter on λsat is assessed. It is shown that these factors have a major impact on λsat.