“…Hence, they can modulate the memory of the unconfined aquifer flow, which, in turn, can modulate the memory of the watershed baseflow. Meanwhile, the Caputo derivative, as defined in its special form D β x i 0 f (x i ) in space in this study, was shown by Kavvas and Ercan (2016) and Ercan and Kavvas (2017) to be a nonlocal quantity where the effect of the boundary conditions on the groundwater flow within the flow domain can have long spatial memories with the decrease in the powers of the spatial fractional derivatives from unity. Similarly, it was shown by Kavvas et al (2017a) that the Caputo derivative in time, D α 0 f (t), as defined in this study, is nonlocal in time and can carry the effect of initial conditions on the groundwater flow for long time periods as the power in the time fractional derivative decreases from 1.…”