The porous medium which groundwater flows in is heterogeneous at all scales. This complicates the simulation of groundwater flow. Fractional derivatives, because of their non-locality property, can reduce the scale effects ...