This paper compares six methods for data-driven rainfall–runoff simulation in terms of predictive accuracy, error structure, interpretability, and uncertainty. We demonstrate that autocorrelation in model errors can result in biased estimates of important values and show how certain model structures can be more easily interpreted to yield insights on physical watershed function. Finally, we explore how model structure can impact uncertainty in climate change sensitivity estimates.