Turbulence exhibits a striking duality: it drives concentrated substances apart, enhancing mixing and transport, while simultaneously drawing particles and bubbles into collisions. Little experimental data exist to clarify the latter process due to challenges in techniques for resolving bubble pairs from afar to coalescence via turbulent entrainment, film drainage and rupture. In this work, we tracked pairs of bubbles across nearly four orders of magnitude in spatial resolution, capturing the entire dynamics of collision and coalescence. The resulting statistics show that critical variables exhibit scalings with bubble size in ways that are different from some classical models, which were developed based on assumptions that bubble collision and coalescence only mirror the key scales of the surrounding turbulence. Furthermore, contrary to classical models which suggest that coalescence favours slow collision velocity, we find a ‘Goldilocks zone’ of relative velocities for bubble coalescence, where there is an optimal coalescence velocity that is neither too high nor too low. This zone arises from the competition between bubble–bubble and bubble–eddy interactions. Incorporating this zone into the new model yields excellent agreement with experimental results, laying a foundation for better predictions for many multiphase flow systems.