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