All natural surfaces exhibit nanoscale roughness (NR) and chemical heterogeneity (CH) to some extent. Expressions were developed to determine the mean interaction energy between a colloid and a solid water interface, as well as for colloid colloid interactions, when both surfaces contain binary NR and CH. The influence of heterogeneity type, roughness parameters, solution ionic strength (IS), mean zeta potential, and colloid size on predicted interaction energy profiles was then investigated. The role of CH was enhanced on smooth surfaces with larger amounts of CH, especially for smaller colloids and higher IS. However, predicted interaction energy profiles were mainly dominated by NR, which tended to lower the energy barrier height and the magnitudes of both the secondary and primary minima, especially when the roughness fraction was small. This dramatically increased the relative importance of primary to secondary minima interactions on net electrostatically unfavorable surfaces, especially when roughness occurred on both surfaces and for conditions that produced small energy barriers (e.g., higher IS, lower pH, lower magnitudes in the zeta potential, and for smaller colloid sizes) on smooth surfaces. The combined influence of roughness and Born repulsion frequently produced a shallow primary minimum that was susceptible to diffusive removal by random variations in kinetic energy, even under electrostatically favorable conditions. Calculations using measured zeta potentials and hypothetical roughness properties demonstrated that roughness provided a viable alternative explanation for many experimental deviations that have previously been attributed to electrosteric repulsion (e.g., a decrease in colloid retention with an increase in solution IS; reversible colloid retention under favorable conditions; and diminished colloid retention and enhanced colloid stability due to adsorbed surfactants, polymers, and/or humic materials).