“…Because the spatial variations in hydraulic parameters such as hydraulic conductivity, in turn, cause the spatial variation in the pore water velocity, the implication is that the first‐order approximation should be applicable for solute transport in mildly heterogeneous media ( ; e.g., Gelhar, ; Rubin, ). In addition, the use of the first‐order approximation (Equation ) to evaluate the displacement covariance or travel time statistics has been validated by Monte Carlo simulations for the cases of nonstationary flow field (e.g., Butera & Tanda, , Wu & Hu, ) or radial flow field (e.g., Lu & Zhang, ; Stauffer, Attinger, Zimmermann, & Kinzelbach, ) in heterogeneous media. The results of Monte Carlo simulations in above‐mentioned references indicate that reasonable agreement with those obtained using the first‐order approximation theory for the case of the variance of log‐hydraulic conductivity <1.…”