Abstract. Land surface models such as JULES (the Joint UK Land Environment Simulator) are usually run on a regular, rectilinear grid, resulting in gridded outputs for variables such as soil moisture and water fluxes. Here we investigate a method of clustering grid cells with similar characteristics together in JULES. Clustering grid cells has the potential to reduce computational expense as well as providing an alternative to tiling approaches for capturing sub-grid heterogeneity. In this study, we cluster grid cells exclusively in the land surface part of modelling, i.e., separate from river routing. We compare gridded and clustered soil moisture outputs from JULES with measurements from the UK Centre for Ecology and Hydrology (UKCEH) COSMOS-UK network and show that the clustering approach can model soil moisture well while reducing computational expense. However, soil moisture results are dependent on the characteristics used to create the clusters. We investigate the effect of using clusters on predicted river flows, and compare routed JULES outputs with NRFA gauge data in the catchment. We show that less expensive JULES clustered outputs give similar river flow results to standard gridded outputs when routed at the grid resolution, and are able to match observed river flow better than gridded outputs when routed at higher resolution.