A three-dimensional thermal lattice Boltzmann model (TLBM) using multirelaxation time method was used to simulate stratified atmospheric flows over a ridge. The main objective was to study the efficacy of this method for turbulent flows in the atmospheric boundary layer, complex terrain flows in particular. The simulation results were compared with results obtained using a traditional finite difference method based on the Navier-Stokes equations and with previous laboratory results on stably stratified flows over an isolated ridge. The initial density profile is neutral stratification in the boundary layer, topped with a stable cap and stable stratification aloft. The TLBM simulations produced waves, rotors, and hydraulic jumps in the lee side of the ridge for stably stratified flows, depending on the governing stability parameters. The Smagorinsky turbulence parameterization produced typical turbulence spectra for the velocity components at the lee side of the ridge, and the turbulent flow characteristics of varied stratifications were also analyzed. The comparison of TLBM simulations with other numerical simulations and laboratory studies indicated that TLBM is a viable method for numerical modeling of stratified atmospheric flows. To our knowledge, this is the first TLBM simulation of stratified atmospheric flow over a ridge. The details of the TLBM, its implementation of complex boundaries and the subgrid turbulence parameterizations used in this study are also described in this article.