В области криолитозоны на фоне изменения климата и возрастающей техногенной нагрузки с каждым годом возрастает актуальность изучения условий формирования ресурсов и режима над- и межмерзлотных вод, которые наиболее чувствительны к внешнему воздействию. Одним из методов изучения динамики подземных вод при невозможности проведения систематических круглогодичных наблюдений является численное моделирование. Цель выполненных исследований – создать математическую модель фильтрации надмерзлотно-межмерзлотных подземных вод, широко распространенных на правобережье р. Лена в Центральной Якутии, и оценить межгодовую изменчивость их ресурсов. Трехмерная стандартная геофильтрационная модель была создана в программе MODFLOW-USG. Исходными данными для ее построения послужили результаты многолетних полевых гидрогеологических и геокриологических исследований, полученные на водосборной площади одного из круглогодично действующих наледеобразующих источников на бестяхской террасе р. Лена. Для калибровки численной модели использованы данные натурных наблюдений за уровнем воды в скважине, расположенной в зоне транзита межмерзлотных вод. Корректировка фильтрационных параметров и граничных условий модели проводилась до тех пор, пока расхождение между расчетными и фактическими значениями уровня межмерзлотных вод не составило менее 0,1 м для более чем 90 % измерений. Результаты моделирования позволили оценить динамику водообильности межмерзлотного водоносного талика за последние 15 лет, уточнить влияние различных метеорологических факторов и процессов сезонного промерзания-протаивания пород на формирование режима подземных вод, количественно оценить внутригодовую изменчивость дебита источника. Полученные результаты могут послужить основой для постановки задач по прогнозированию изменения мерзлотно-гидрогеологических условий под действием колебаний климата, а также рассматриваться как эффективный метод количественной оценки динамики подземных вод в области распространения многолетнемерзлых пород.