In this work we perform global fits of microscopic decoherence models of neutrinos to all available current data, including LSND and KamLAND spectral distortion results. In previous works on related issues the models used were supposed to explain LSND results by means of quantum gravity induced decoherence. However those models were purely phenomenological without any underlying microscopic basis. It is one of the main purposes of this article to use detailed microscopic decoherence models with complete positivity, to fit the data.The decoherence in these models has contributions not only from stochastic quantum gravity vacua operating as a medium, but also from conventional uncertainties in the energy of the (anti)neutrino beam. All these contributions lead to oscillation-length independent damping factors modulating the oscillatory terms from which one obtains an excellent fit to all available neutrino data, including LSND and Kamland spectral distortion.