“…Chang and Yeh () developed relationships between the stochastic structures of ln K and those of induced excess pressure head based on the assumption that a load decreasing linearly with time is applied vertically on the upper part of a finite‐thickness aquifer. Their results governing the groundwater flow field will be applied to investigate the behaviour of variability of solute travel time, which are given by (Eqs 17 and 31 of Chang & Yeh, ) where Φ represents the mean of excess pressure head, and ϕ represents small pressure head fluctuations about the mean. In addition, Λ denotes the uniform decreasing rate of vertical loading, η = K G t /( SL 2 ), S is the storage coefficient defined as S = ρ f g ( β p + nβ f ), ρ f is the density of fluid, β p and β f are the coefficients of soil compressibility and fluid compressibility, respectively, K G = exp[ F ], F is the mean of log hydraulic conductivity, n is the porosity, L is the vertical length of the flow domain, ξ = Z / L , and …”