Advective-dispersive equation with spatial fractional derivatives evaluated with tracer transport data

The classical model used to describe solute transport in soil is based on the advective-dispersive equation (ADE) in which an analog of Fick's law is used to model dispersion. The fractional ADE (FADE) has been proposed to address discrepancies between experimental solute concentrations and those predicted with the ADE. The order of the fractional derivative or Levy exponent, alpha , characterizes the deviation of the FADE solutions from the ADE, which is a specific case of the FADE for alpha = 2. The objective of this work was to test the hypothesis that using the FADE with values of alpha other than 2 can provide more accurate simulations of solute transport in soils. We fitted the FADE to 47 published data sets on tracer breakthrough in disturbed and undisturbed soil columns. The FADE was solved numerically with a mass-conserving boundary condition. While 19 breakthrough curves were best fitted with the ADE, the rest were fitted better using the FADE with alpha <2. The RMSEs of the FADE and ADE were close when the FADE alpha was >1.5. In contrast, the FADE RMSEs were, on average, 1.5 times smaller than the ADE RMSE when the FADE was fitted with alpha < 1.5. The value of alpha apparently reflected the structure of the void space available for flow and transport. Considered as a generalization of the classical ADE, the FADE can be a useful model if tails of the breakthrough curve are of special interest.