Time and time again, transport behaviors that cannot be described by traditional simple advection dispersion equations (ADE) are observed in hydrologic systems, both in surface and subsurface environments. Such departures from classical predictions have been termed anomalous. In all cases, the presence of processes and heterogeneities spanning a broad range of temporal and spatial scales leads to a breakdown in assumptions inherent to classical ADEs and complicates modeling efforts. Developing approaches that account for these complexities naturally leads to the use of upscaling and stochastic modeling frameworks. We briefly introduce these ideas from both Lagrangian and Eulerian perspectives. We start with classical ADE models, highlighting many of the assumptions they are built on. By relaxing these assumptions, models capable of reproducing observed anomalous behaviors have been developed. Here, we briefly review three of the most popular models: continuous time random walks (CTRW), fractional advection dispersion equations (fADE), and multirate mass transfer (MRMT). While all three have enjoyed success, they too have their limitations. We conclude this article with a discussion of these shortcomings and present three very recent models that aim to address them as well as a brief discussion on future challenges.