i.e., HadCM3Q0, Q3, Q10, Q13, and Q15). The ensemble simulations are then synthesized through a Bayesian hierarchical model to develop probabilistic projections of future temperature outcomes with consideration of some uncertain parameters involved in the regional climate modeling process. The results suggest that there would be a consistent increasing trend in the near-surface air temperature with time periods from 2030s to 2080s. The most likely mean temperature in next few decades (i.e., 2030s) would be [−2, 2] °C in northern Ontario, [2, 6] °C in the middle, and [6, 12] °C in the south, afterwards the mean temperature is likely to keep rising by ~ 2 °C per 30-years period. The continuous warming across the Province would drive the lowest mean temperature up to 2 °C in the north and the highest mean temperature up to 16 °C in the south. In addition, the spread of the most likely ranges of future outcomes shows a consistent widening trend from 2030s to 2080s, implying that long-term climate change is more difficult to predict than near-term change because many more uncertain or unknown factors may continue to emerge during the long-term simulation.