Nuclear magnetic resonance (NMR) logging provides a relatively new approach for estimating the hydraulic conductivity κ of unconsolidated aquifers. We have evaluated results from model validation and uncertainty quantification using direct push measurements of NMR mean relaxation times and κ in sands and gravels at three field sites.We have tested four models that have been proposed for predicting κ from NMR data, including the Schlumberger-Doll research, Seevers, and sum-of-echoes equations, all of which use empirically determined constants, as well as the Kozeny-Godefroy model, which predicts κ from several physical parameters. We have applied four methods of analysis to reanalyze NMR and κ data from the three field sites to quantify how well each model predicted κ from the mean log NMR relaxation time T2ML given the uncertainties in the data. Our results show that NMR-estimated porosity does not improve prediction of κ in our data set for any model and that all of the models can predict κ to within an order of magnitude using the calibrated constants we have found. We have shown the value of rigorous uncertainty quantification using the methods we used for analyzing κ-NMR data sets, and we have found that incorporating uncertainty estimates in our analysis gives a more complete understanding of the relationship between NMR-derived parameters and hydraulic conductivity than can be obtained through simple least-squares fitting. There is little variability in our data set in the calibrated constants we find, given the uncertainty present in the data, and therefore we suggest that the constants we find could be used to obtain first-order estimates of hydraulic conductivity in unconsolidated sands and gravels at new sites with NMR data available.