Abstract. Process-based geomorphic transport laws enable to assess the impact of rainfall variability on bedrock river incision over geological timescales. However, isolating the role of rainfall variability on erosion remains difficult in natural environments in part because the variability of rock strength and its resistance to incision are poorly constrained. Here, we explore spatial differences in the rate of bedrock river incision in the Tropical Andes. The Ecuadorian Andes are characterized by strong rainfall gradients due to orographic precipitation sourced in the Amazon basin. In addition, the tectonic configuration has generated a profound lithological heterogeneity. The relative role of either these controls in modulating river incision on millennial time scales, however, remains unclear. Using ¹⁰Be catchment-wide erosion rates, meteorological and hydrological data, as well as data on bedrock erodibility, we provide quantitative constraints on the importance of rainfall variability and lithological variations. Explicit incorporation of rock erodibility in river incision models predicated on the stream power equation enables us to identify a first order control of lithology on river incision rates. Rainfall variability based on a spatially and temporally explicit hydrological dataset and a stochastic-threshold river incision model explain regional differences in river incision that cannot be attributed to topographical and/or lithological variability.