Time and space nonlocalities underlying fractional-derivative models: Distinction and literature review of field applications

We investigate the spatiotemporal nonlocality underlying fractional-derivative models as a possible explanation for regional-scale anomalous dispersion with heavy tails. Properties of four fractional-order advection–dispersion equation (fADE) models were analyzed and compared systematically, including the space fADEs with either maximally positive or negative skewness, the time fADE with a temporal fractional-derivative 0<γ<10<γ<1, and the extension of the time fADE with 1<γ<21<γ<2. Space fADEs describe the dependence of local concentration change on a wide range of spatial zones (i.e., the space nonlocality), while time fADEs describe dynamic mass exchange between mobile and multiple immobile phases and therefore record the temporal history of concentration “loading” (i.e., the time-nonlocality). We then applied the fADEs as models of anomalous dispersion to four extensively-studied, regional-scale, natural systems, including a hillslope composed of fractured soils, a river with simultaneous active flow zones and various dead-zones, a relatively homogeneous glaciofluvial aquifer dominated by stratified sand and gravel, and a highly heterogeneous alluvial aquifer containing both preferential flowpaths and abundant aquitards. We find that the anomalous dispersion observed at each site might not be characterized reasonably or sufficiently by previous studies. In particular, the use of the space fADE with less than maximally positive skewness implies a spatial dependence on downstream concentrations that may not be physically realistic for solute transport in watershed catchments and rivers (where the influence of dead-zones on solute transport can be described by a temporal, not spatial, fractional model). Field-scale transport studies show that large ranges of solute displacement can be described by a space nonlocal, fractional-derivative model, and long waiting times can be described efficiently by a time-nonlocal, fractional model. The unknown quantitative relationship between the nonlocal parameters and the heterogeneity, and the similarity in concentration profiles that are solutions to the different nonlocal transport models, all demonstrate the importance of distinguishing the representative nonlocality (time and/or space) for any given regional-scale anomalous dispersion process.