This paper proposes a Bayesian framework for localization of multiple sources in the event of accidental hazardous contaminant release. The framework assimilates sensor measurements of the contaminant concentration with an integrated multizone computational fluid dynamics (multizone-CFD) based contaminant fate and transport model. To ensure online tractability, the framework uses deep Gaussian process (DGP) based emulator of the multizone-CFD model. To effectively represent the transient response of the multizone-CFD model, the DGP emulator is reformulated using a matrix-variate Gaussian process prior. The resultant deep matrix-variate Gaussian process emulator (DMGPE) is used to define the likelihood of the Bayesian framework, while Markov Chain Monte Carlo approach is used to sample from the posterior distribution. The proposed method is evaluated for single and multiple contaminant sources localization tasks modeled by CONTAM simulator in a single-story building of 30 zones, demonstrating that proposed approach accurately perform inference on locations of contaminant sources. Moreover, the DMGP emulator outperforms both GP and DGP emulator with fewer number of hyperparameters.arXiv:1806.08069v1 [stat.AP] 21 Jun 2018 approaches do not provide full uncertainty analysis of contaminant source location and characteristics. Primarily due to their ability to fully quantify the uncertainties, the Bayesian framework for contaminant source localization and characterization is gaining prominence in various studies [16][17][18][19]. However, barring few notable exceptions of conjugacy like linear models with Gaussian priors, estimation of the Bayesian posterior distribution is analytically intractable. Thus, the Bayesian framework is often implemented by approximating the posterior distribution using sampling methods like Markov Chain Monte Carlo (MCMC) [20].Implementation of the MCMC based Bayesian inference framework for contaminant source localization and characterization is challenging due to: 1) for acceptable accuracy, MCMC requires 10 4 − 10 7 samples [16,19,21,22]; 2) transient nature of the contaminant dispersion phenomena [19,23]. Note that each MCMC sample requires simulating the contaminant fate and transport model, rendering the Bayesian framework computationally intractable for real-time contaminant source localization and characterization [16]. To resolve the first challenge, many state of the art approaches relies on simplification of the contaminant fate and transport model [18]. However, the computational cost of these simplified models is still prohibitive for online applications. Moreover, these models use model order reduction which results in a loss of spatial information.Alternatively, recent studies have explored various machine learning algorithms to develop computationally efficient emulators to infer the indoor events [24]. A similar approach is adopted in this paper, where the contaminant fate and transport model is replaced by a statistical emulator in the Bayesian framework. Notwithstanding t...