Gully erosion poses a severe threat to ecological security and economic sustainability globally. Soil erodibility, characterizing the susceptibility to gully erosion, is influenced by hydrothermal conditions that control pedogenic processes, especially for soils derived from the same parent rock. Soil erodibility characteristics in gully erosional areas have been widely investigated, particularly at a pedon scale. However, limited information is available on the regulation pathway of hydrothermal conditions on gully erosion. Therefore, gully erosion distributed in the granite weathering mantle was selected due to its increased density and sediment yield from north to south across subtropical China. To evaluate the impacts of granitic soil erodibility on the spatial distribution of gully erosion, a comprehensive soil erodibility index (CSEI) was constructed. This index integrated various soil erodibility indicators, including mean weight diameter (MWD), soil structural stability index (SSSI), penetration resistance (PR), soil erodibility factor (K-factor), and others. These soil erodibility indicators with inconsistent spatial variation across the study area showed the highest values at the sandy horizon (F=1.72 ∼ 6.53, p < 0.05), indicating the critical role of this horizon to gully erosion. While a consistent increasing trend with the increase of hydrothermal condition was observed, this spatial variation showed a less significant difference statistically. Moreover, the CSEI (0.12 ∼ 0.91) was mainly determined by cementation (R2 = 0.65 ∼ 0.83, p < 0.01) and hydrophysical properties (LL, PL, w) (R2 = 0.64 ∼ 0.69, p < 0.01), especially soil organic matter and free iron and aluminum oxides (R2 = 0.91). Furthermore, hydrothermal conditions indirectly determined the spatial variation of CSEI by altering soil cementation (SPC=0.51) and physical-related properties (SPC= − 0.15). Collectively, hydrothermal conditions determine gully erosion by altering granitic soil erodibility in subtropic China, and this alteration was triggered by an indirect way. The obtained results would facilitate a better understanding of gully erosion mechanisms and more effective gully control strategies.