In this work, we developed a novel qGP-DRT framework that enables the simultaneous deconvolution of DRTs from multiple EIS datasets. This new model inherits key features from the existing fGP-DRT model, including the optimization of its hyperparameters through evidence maximization, robustness against experimental noise, non-negativity of the recovered DRTs, and utilization of the complete impedance spectra. However, unlike fGP-DRT, the qGP-DRT model