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 on the parameters and, consequently, better simulate the hydrogeological processes. In this chapter a fractional governing partial differential equation on unconfned groundwater (fractional Boussinesq equation [FBE]) is derived using the fractional mass conservation law. The FBE is a generalization of the Boussinesq equation (BE) that can be used in both homogeneous and heterogeneous unconfned aquifers. Compared to the BE, the FBE includes an additional parameter which represents the heterogeneity degree of the porous medium. This parameter changes within the range of 0 to 1 in the non-linear form of the FBE. The smaller the value of the heterogeneity degree, the more heterogeneous the aquifer is, and vice versa. To investigate the applicability of the FBE to real problems in groundwater flow, a fractional Glover–Dumm equation (FGDE) was obtained using an analytical solution of the linear form of the FBE for onedimensional unsteady flow towards parallel subsurface drains. The FGDE was ftted to water table profles observed at laboratory and feld scales, and its performance was compared to that of the Glover–Dumm equation (GDE). The parameters of the FGDE and the GDE were estimated using the inverse problem method. The results indicate that one can recognize the heterogeneity degree of porous media examined according to the obtained values for the indicator of the heterogeneity degree. The FGDE and the GDE showed similar performances in homogeneous soil, while the performance of the FGDE was signifcantly better than that of the GDE in heterogeneous soil. In summary, the FBE can be used as a highly general differential equation governing groundwater flow in unconfned aquifers.
