Solute transport in bounded porous media characterized by generalized sub-Gaussian log-conductivity distributions

There are increasing evidences that probability distributions and associated statistical moments of a variety of hydrogeological and soil science variables and their spatial increments display distinctive scale-dependent features that are not captured by a typical Gaussian model. A Generalized Sub-Gaussian (GSG) model is able to capture key aspects of this pattern. We present the results of a suite of computational analyses set in a Monte Carlo framework and aimed at assessing the impact of a GSG structure of log hydraulic conductivity (Y) on transport of a conservative solute through a three-dimensional bounded porous medium under steady-state saturated Darcy flow. Our results indicate that the longitudinal spreading of a plume is on average significantly smaller for Sub-Gaussian than for Gaussian Y fields. Otherwise, the velocity field arising from a Sub-Gaussian Y field induces enhanced plume stretching with respect to what can be observed in a Gaussian Y setting, this aspect potentially influencing the strength of solute mixing within these two types of conductivity domains. We also find that, in some cases, it may be difficult to identify the nature of the underlying conductivity field relying solely on observations of solute concentrations migrating within the system. In this regard, we show that the action of local dispersion tends to mask the influence of Sub-Gaussianity on major transport metrics.