Представлены результаты математического моделирования и прогнозирования динамики полей термокарстовых озер на территории многолетней мерзлоты Западной Сибири до конца нынешнего столетия. Описана геоимитационно-регрессионная модель динамики термокарстово-озерных полей, свойства которой определены на основе статистического распределения координат местоположения и площадей озер, полученных экспериментальным путем с использованием разновременных космических снимков Landsat за период 1973-2009 гг. В модели учтена взаимосвязь термокарстово-озерных и климатических изменений, установленная на основе многомерного регрессионного анализа экспериментальных данных. Прогнозирование динамики полей термокарстовых озер в зоне мерзлоты Западной Сибири осуществлено с помощью компьютерных экспериментов на модели с учетом известных прогнозов изменения среднегодовой температуры для исследуемой территории. Показано, что повышение температуры будет сопровождаться сокращением площади термокарстовых озер в зоне мерзлоты Западной Сибири.