Первый член правой части уравнения описывает миграцию сульфатов в почве в результате действия процессов диффузии. К — коэффициент конвективной диффузии в пористой среде и равен.
Синтез целостной модели, представленной системой дифференциальных уравнений, решение которой будет имитировать динамику процессов транспорта серы в почве лесной экосистемы, осуществлялся на основе закона сохранения вещества: скорость изменения количества серы в резервуаре равна сумме входящих потоков минус сумма исходящих потоков. Таким балансовым методом была получена следующая система дифференциальных уравнений динамики переменных состояния:
Для получения решения этой системы дифференциальных уравнений был применен метод конечно-разностной аппроксимации, который позволяет свести задачу к решению системы алгебраических уравнений, что может быть эффективно сделано с помощью ЭВМ. Для решения всей системы уравнений указанным методом была написана программа на алгоритмическом языке ФОРТРАН, которая была реализована на ЕС-ЭВМ. Шаг по времени 10 сут, прогнозирование велось на 10 лет. В результате работы по программе рассчитывается сезонная динамика переменных состояния модели и интенсивности процессов, слагающих БГХЦ серы в почве лесной экосистемы для двух случаев: «чистого» участка с поступлением серы с атмосферными осадками 12 кг/га в год и загрязненного с поступлением серы из атмосферы, равным 42 кг/га в год, и с растительным опадом.