Описана модернизированная численная модель распространения набора спектральных гармоник внутренних гравитационных волн (ВГВ) в неоднородной атмосфере от земной поверхности до нижней термосферы. Учтена диссипация ВГВ из-за турбулентной вязкости и теплопроводности, радиационного теплообмена и ионного трения в нижней ионосфере. Выполнено численное моделирование распространения спектра ВГВ в фоновых полях ветра и температуры, соответствующих координатам станций Звенигород (56° с.ш., 37° в.д.), Маймага (63° с.ш., 130° в.д.) и Торы (52° с.ш., 103° в.д.), на которых проводятся систематические наблюдения ночной эмиссии гидроксила. Рассчитаны сезонные вариации стандартных отклонений возмущений компонент горизонтальной скорости, создаваемых модельным спектром ВГВ. Обнаружено общее сходство сезонных вариаций модельных амплитуд ВГВ и наблюдаемой дисперсии мезомасштабных возмущений вращательной температуры гидроксила в Звенигороде и Торах. Это является свидетельством того, что интенсивность мезомасштабных возмущений температуры вблизи мезопаузы может зависеть от интенсивности ВГВ, распространяющихся из нижней атмосферы, и от профилей фоновых характеристик средней атмосферы на пути распространения волновых пакетов в различные сезоны и в разных географических пунктах.