Разработана математическая модель прогноза температурного режима породного целика и развала многолетнеморозных горных пород после взрывной отбойки, учитывающая особенности криолитозоны и динамику горных работ. На примере Кангаласского буроугольного месторождения Якутии проведены расчеты прогноза температурного режима массива многолетнемерзлых горных пород до и после взрывной отбойки, а также в процессе экскавации взорванных пород вскрыши. Установлено, что на формирование температурного поля во взорванном массиве горных пород оказывает влияние период проведения взрывных работ. Результаты математического моделирования показывают существенную неоднородность распределения температурного поля в развале взорванных горных пород, что должно учитываться при выборе схемы его отработки.