Realistic hydrological response is sensitive to the spatial variability of landscape properties. For a grid-based distributed rainfall-runoff model with a hypothesis of a uniform grid, the high-frequency information within a grid cell will be gradually lost as the resolution of the digital elevation model(DEM) grows coarser. Therefore, the performance of a hydrological model is usually scale-dependent. This study used the GridXinanjiang(GXAJ) model as an example to investigate the effects of subgrid variability on hydrological response at different scales. With the aim of producing a more reasonable hydrological response and spatial description of the landscape properties, a new distributed rainfall-runoff model integrating the subgrid variability(the GXAJSV model) was developed. In this model, the topographic index is used as an auxiliary variable correlated with the soil storage capacity. The incomplete beta distribution is suggested for simulating the probability distribution of the soil storage capacity within the raster grid. The Yaogu Basin in China was selected for model calibration and validation at different spatial scales.Results demonstrated that the proposed model can effectively eliminate the scale dependence of the GXAJ model and produce a more reasonable hydrological response. For a grid-based distributed rainfall-runoff model with a hypothesis of a uniform grid, the high-frequency information within a grid cell will be gradually lost as the resolution of the digital Therefore, the performance of a hydrological model is usually scale-dependent. This study used the GridXinanjiang (GXAJ) model as an example to investigate the effects of subgrid variability on hydrological response at different scales. With the aim of producing a more reasonable hydrological response and spatial description of the landscape properties, a new distributed rainfall-runoff model integrating the sub-variability (the GXAJSV model) was developed. In this model, the topographic index is used as an auxiliary variable correlated with the soil storage capacity. The incomplete beta distribution is suggested for simulating the probability distribution of the The storage of within the raster grid. The Yaogu Basin in China was selected for model calibration and validation at different spatial scales. Results demonstrated that the proposed model can effectively eliminate the scale dependence of the GXAJ model and produce a more reasonable hydrological response.