Comment on “Modeling groundwater flow with random hydraulic conductivity using radial basis function partition of unity method” by Shile et al. (2024)

The goal of this article is to explain some striking discrepancies between the Monte Carlo inferences for flow in heterogeneous aquifers presented in (Shile et al., 2024) and reliable and verifiable results previously published. Comparisons with statistical inferences done within the same numerical setup and using the same benchmark codes demonstrate that quantifiable aspects, such as a too small departure from the linear theory and the decrease of the standard deviations with the increase of the aquifer heterogeneity, cannot be produced with the tools used by the authors. • Monte Carlo results of Shile et al. (2024) are compared to previously published ones. • Comparisons are also made with other methods for the same groundwater-flow benchmark. • Striking discrepancy with reliable and verifiable results indicate data fabrication.