We propose a simple and robust approach for investigating uncertainty in the results of inversion in geophysics. We apply this approach to inversion of Surface Nuclear Magnetic Resonance (SNMR) data, which is also known as Magnetic Resonance Sounding (MRS). Solution of this inverse problem is known to be non-unique. We inverse MRS data using the well-known Tikhonov regularization method, which provides an optimal solution as a trade-off between the stability and accuracy. Then, we perturb this model by random values and compute the fitting error for the perturbed models. The magnitude of these perturbations is limited by the uncertainty estimated with the singular value decomposition (SVD) and taking into account experimental errors. We use 106 perturbed models and show that the large majority of these models, which have all the water content within the variations given by the SVD estimate, do not fit data with an acceptable accuracy. Thus, we may limit the solution space by only the equivalent inverse models that fit data with the accuracy close to that of the initial inverse model. For representing inversion results, we use three equivalent solutions instead of the only one: the “best” solution given by the regularization or other inversion technic and the extreme variations of this solution corresponding to the equivalent models with the minimum and the maximum volume of water. For demonstrating our approach, we use synthetic data sets and experimental data acquired in the framework of investigation of a hard rock aquifer in the Northern Ireland (County Donegal).