Can a Time Fractional-Derivative Model Capture Scale-Dependent Dispersion in Saturated Soils?

Time nonlocal transport models such as the time fractional advection-dispersion equation (t-fADE) were proposed to capture well-documented non-Fickian dynamics for conservative solutes transport in heterogeneous media, with the underlying assumption that the time nonlocality (which means that the current concentration change is affected by previous concentration load) embedded in the physical models can release the effective dispersion coefficient from scale dependency. This assumption, however, has never been systematically examined using real data. This study fills this historical knowledge gap by capturing non-Fickian transport (likely due to solute retention) documented in the literature (Huang et al. 1995) and observed in our laboratory from small to intermediate spatial scale using the promising, tempered t-fADE model. Fitting exercises show that the effective dispersion coefficient in the t-fADE, although differing subtly from the dispersion coefficient in the standard advection-dispersion equation, increases nonlinearly with the travel distance (varying from 0.5 to 12 m) for both heterogeneous and macroscopically homogeneous sand columns. Further analysis reveals that, while solute retention in relatively immobile zones can be efficiently captured by the time nonlocal parameters in the t-fADE, the motion-independent solute movement in the mobile zone is affected by the spatial evolution of local velocities in the host medium, resulting in a scale-dependent dispersion coefficient. The same result may be found for the other standard time nonlocal transport models that separate solute retention and jumps (i.e., displacement). Therefore, the t-fADE with a constant dispersion coefficient cannot capture scale-dependent dispersion in saturated porous media, challenging the application for stochastic hydrogeology methods in quantifying real-world, preasymptotic transport. Hence improvements on time nonlocal models using, for example, the novel subordination approach are necessary to incorporate the spatial evolution of local velocities without adding cumbersome parameters.