Abstract. The amount of solar radiation reaching the Earth surface (surface solar irradiance, SSI) is critical for a variety of applications, ranging from surface–atmosphere interactions to solar energy. SSI is characterized by a large spatiotemporal variability, in particular in the presence of cumulus clouds. This results in complex spatial patterns of shadows and sunlight directly related to clouds' geometry and physical properties. Although key in many respects, the instantaneous spatial distribution of SSI remains largely unexplored. Here, we use unique observations from a dense network of pyranometers deployed during the HOPE field campaign to investigate the SSI spatial distribution. For cumulus scenes, bimodal distributions are found, with one mode corresponding to cloud shadows and the other to sunlit areas with enhanced SSI exceeding clear-sky values. Combining large-eddy simulations of cumulus clouds with Monte Carlo ray tracing, we demonstrate the capability of advanced numerical tools to reproduce the observed distributions and quantify the impact of cloud geometrical and physical properties on both modes. In particular, cloud cover strongly modulates their amplitudes, in addition to their position and width, which are also sensitive to cloud height, geometrical depth, and liquid water content. Combining observations and simulations, we also explore sampling strategies to estimate the SSI spatial distribution with a limited number of sensors, suggesting that 10 pyranometers integrated over 10 min can capture most details of the full distribution. Such a strategy could be used for future campaigns to further investigate SSI distributions and their impact on land–atmosphere exchanges or photovoltaic farm management.