Dengue is an important emerging vector-borne disease in Bhutan. This study aimed to quantify the spatial and temporal patterns of dengue and their relationship to environmental factors in dengue-affected areas at the sub-district level. A multivariate zero-inflated Poisson regression model was developed using a Bayesian framework with spatial and spatiotemporal random effects modelled using a conditional autoregressive prior structure. The posterior parameters were estimated using Bayesian Markov Chain Monte Carlo simulation with Gibbs sampling. A total of 708 dengue cases were notified through national surveillance between January 2016 and June 2019.