The regulation and spatial differences of urban runoffs are of great concern in contemporary hydrological research. However, owing to a shortage of basic data sources and restrictions on urban hydrological simulation functions, simulating and investigating the regulation mechanism behind rainfall-runoff processes remain significantly challenging. In this study, the Time Variant Gain Model (TVGM), a hydrological nonlinear system model, was extrapolated to the hydrodynamic model of an urban drainage network system by integrating it with the widely used Stormwater Management Model (SWMM) to adequately simulate urban runoff events while considering various underlying surfaces and runoff routing modes, such as surface, drainage network and river runoff, in urban regions (i.e., TVGM-SWMM). Moreover, runoff events were characterized using the following four runoff regime metrics: runoff coefficient, capture ratio of annual runoff volume, standardized flood timescale, and the ratio of occurrence time differences between flow and rainfall peak to event duration (peak flow delay time). The characteristics and spatial differences of urban runoff regulations were investigated, and the key impact factors and their relative contributions were identified using multivariate statistical analyses. Four communities were selected as our study areas, consisting of communities from Beijing, Shenzhen, Wuhan, and Chongqing. Our results showed that the TVGM-SWMM performed considerably better than SWMM alone. The comprehensive simulation accuracy of 60% of the events (12/20) improved by 486%, with the bias improving the most, followed by the efficiency coefficient. Barring the runoff coefficient, significant spatial differences were identified at the patch scale for the runoff regime metrics, with differences of 0.43, 0.22, and 0.16 ( p <0.05). The key impact factors were the pipe length ( r =0.51) in the drainage network system and the forest area ratios ( r =0.56), sponge measures ( r =0.52), grassland ( r =0.48), and impervious surface ( r =0.46) in the underlying surfaces. The contributions of the drainage network system and the underlying surfaces were 4.27% and 37.83%, respectively. Regulation in the Beijing community, dominated by grassland regulation, delayed and reduced the peak flow and total runoff volume. In the Shenzhen community, sharp and thin runoff events were mainly generated by impervious surfaces and were not adequately regulated. Forest regulation was the dominant regulation type in the Wuhan community, which reduced the total runoff volume and delayed the peak flow. Waterbody regulation was the primary regulation type in the Chongqing community, which reduced the total runoff volume and peak flow. This study aims to introduce a comprehensive theoretical and technical assessment of the hydrological effects of urbanization and the performance of sponge city construction and provide a reference for urban hydrological model improvements in China.