Cold pools formed by precipitating convective clouds are an important source of mesoscale temperature variability. However, their sub-mesoscale (100 m to 10 km) structure has not been studied, which impedes validation of numerical models and understanding of their atmospheric and societal impacts. We quantify temperature variability in observed and simulated cold pools using variograms calculated from dense network observations collected during a field experiment and in high-resolution case study and idealized simulations. The temperature variance in cold pools is enhanced for spatial scales between ~5-15 km compared to pre-cold pool conditions, but the magnitude varies strongly with cold pool evolution and environment. Simulations capture the overall cold pool variogram shape well but underestimate the magnitude of the variability, irrespective of model resolution. Temperature variograms outside of cold pool periods are represented by the range of simulations evaluated here, suggesting that models misrepresent cold pool formation and/or dissipation processes.