Polarons play a major role in the description of optical, electrical and dielectrical properties of several ferroelectric oxides. The motion of those particles occurs by elementary hops among the material lattice sites. In order to compute macroscopic transport parameters such as charge mobility, normal (i.e. Fickian) diffusion laws are generally assumed. In this paper we show that when defect states able to trap the polarons for long times are considered, significant deviations from the normal diffusion behaviour arise. As an example of this behavior, we consider here the case of lithium niobate (LN). This can be considered as a prototypical system, having a rich landscape of interacting polaron types and for which a significant wealth of information is available in literature. Our analysis considers the case of a stoichiometric, defect-free lithium niobate containing a certain concentration of small electron polarons hopping on regular Nb sites, and compares it to the material in congruent composition, which is generally found in real-life applications and which is characterized by a large concentration of antisite NbLi defects. While in the first case the charge carriers are free polarons hopping on a regular Nb sublattice, in the second case a fraction of polarons is trapped on antisite defects. Thus, in the congruent material, a range of different hopping possibilities arises, depending on the type of starting and destination sites. We develop a formalism encompassing all these microscopic processes in the framework of a switching diffusion model which can be well approximated by a mobile–immobile transport model providing explicit expressions for the polaron mobility. Finally, starting from the Marcus–Holstein’s model for the polaron hopping frequency we verify by means of a Monte Carlo approach the diffusion/mobility of the different polarons species showing that, while free polarons obey the laws for normal diffusion as expected, bound polarons follow an anomalous diffusion behaviour and that in the case of the congruent crystal where mixed free and bound polaron transport is involved, our expressions indeed provide a satisfactory description.