Neural ODEs fuse physics-based models with deep learning: neural networks substitute terms in differential equations that represent the mechanistic structure of the system. The approach combines the flexibility of machine learning with physical constraints for inter- and extrapolation. We demonstrate that neural ODE models achieve state-of-the-art predictive performance while keeping full interpretability of model states and processes in hydrologic modelling over multiple catchments.