We derive a phenomenological continuum saltation model for aeolian sand transport that can serve as an efficient tool for geomorphological applications. The coupled differential equations for the average density and velocity of sand in the saltation layer reproduce both the known equilibrium relations for the sand flux and the time evolution of the sand flux as predicted by microscopic saltation models. The three phenomenological parameters of the model are a reference height for the grain-air interaction, an effective restitution coefficient for the grain-bed interaction, and a multiplication factor characterizing the chain reaction caused by the impacts leading to a typical time or length scale of the saturation transients. We determine the values of these parameters by comparing our model with wind tunnel measurements. Our main interest are out of equilibrium situations where saturation transients are important, for instance at phase boundaries (ground/sand) or under unsteady wind conditions. We point out that saturation transients are indispensable for a proper description of sand flux over structured terrain, by applying the model to the windward side of an isolated dune, thereby resolving recently reported discrepancies between field measurements and theoretical predictions.