В последние годы все больший интерес вызывают гибридные методы интерполяции, которые дают лучшие результаты по сравнению с одиночными подходами. Один из таких гибридных методов известен как регрессионный кригинг, который сочетает в себе регрессионную зависимость отклика от вспомогательных переменных с кригингом остатков от регрессии. В данной статье рассмотрена методология цифрового почвенного картографирования с использованием регрессионного кригинга на примере почвенного покрова Лесной опытной дачи (ЛОД) РГАУ-МСХА имени К.А. Тимирязева. Для пространственной интерполяции использовали 87 образцов почв, в которых были определены содержание гумуса, pH водной и солевой вытяжек, а также набор переменных-предикторов (цифровая модель рельефа и производные от нее). В результате построены карты исследуемых свойств почв и изучено влияние характеристик рельефа на пространственную дифференциацию свойств почв. В результате 10-кратной перекрестной проверки показано, что полученная модель для содержания гумуса объясняет 88% дисперсии (RMSE 0.34), модель для величины pH водной вытяжки объясняет 89% дисперсии (RMSE 0.25), модель для величины pH солевой вытяжки - также 89% дисперсии (RMSE 0.31). Моделирование пространственной дифференциации исследуемых свойств почв ЛОД показало, что рельеф территории оказывает значительное влияние на формирование почвенного покрова. Наиболее значимыми переменными в дифференциации содержания гумуса в почвах ЛОД оказались плановая и кросс-секционная кривизны поверхности, значения абсолютных высот, крутизна склонов, относительные превышения в окрестности 100 м, топографический индекс влажности.