CC BY
23
НАУКИ О ЗЕМЛЕ
-I-
УДК [550.46+556.51]:004.94
Ю.Б. Кирста, А.В. Пузанов Yu.B. Kirsta, A.V. Puzanov
МОДЕЛИРОВАНИЕ ГИДРОХИМИЧЕСКОГО СТОКА ГОРНЫХ РЕК: 1. МОДЕЛЬ СТОКА ФОСФАТОВ
MODELING OF HYDROCHEMICAL COMPOSITION OF MOUNTAIN RIVER RUNOFF:
1. MODEL FOR PHOSPHATE RUNOFF
Ключевые слова: системно-аналитическое моделирование, гидрохимический сток, фосфаты, горные реки, Алтай, Саяны.
На примере 34 средних и малых рек Алтае-Саянской горной страны с помощью системно-аналитического моделирования разработана универсальная математическая модель сезонной и многолетней динамики стока фосфатов (анионов PO43"). Входными факторами и переменными модели являются: пространственно обобщенные по Алтае-Саянской горной стране нормированные месячные осадки и среднемесячные температуры воздуха; водные стоки, рассчитываемые для отдельных ландшафтов в речных бассейнах по ранее разработанной модели водного стока горных рек; картографическая информация о площади и средней высоте бассейнов; высоте замыкающих их речных створов; длине речных русел; площади пашни. Для учета ландшафтной структуры речных бассейнов выделены 13 типологических групп геосистем (ландшафтов). Путем решения обратной задачи в ходе системно-аналитического моделирования для четырех гидрологических сезонов года найдены постоянные среднесезонные концентрации фосфатов в водных стоках, поступающих в реки с каждого из выделенных ландшафтов. Для последних во всех 34 речных бассейнах рассчитана сезонная и многолетняя (1951-2003 гг.) динамики стока PO43". Пахотные земли, при их значительных площадях, дают на протяжении всего года заметный вклад в сток фосфатов, что может объясняться использованием на сельскохозяйственных полях минеральных удобрений. Объ-
единенные в один комплекс модели пространственного обобщения и нормировки среднемесячных температур и месячных осадков, водного и гидрохимического стоков позволяют рассчитать сезонную и многолетнюю динамику стока фосфатов для любых речных бассейнов Алтае-Саянской горной страны. Все модели применимы при ограниченной экспериментальной гидрометеорологической и гидрохимической информации по речным бассейнам.
Keywords: system-analytical modeling, hydro-chemical runoff, phosphates, mountain rivers, the Altai, the Sayan.
A universal mathematical model for seasonal and long-term dynamics of phosphates (anions PO43—) runoff was developed based on the system-analytical modeling and the example of 34 mid-size and small rivers of the Altai-Sayan mountain country. The input factors and variables of the model are as following: spatially generalized for the Altai-Sayan mountain country normalized monthly precipitation and mean monthly air temperature; water runoff estimated for individual landscapes in river basins by the previously developed model for water runoff of mountain rivers; the cartographic information on the area and average altitude of the basins; the height of the outlet; the length of river channels; the area of arable land. A total of 13 typological geosystem groups (landscapes) were selected to account for a landscape structure of river basins. The solution of the inverse problem during the system-analytical modeling for the four hydrological seasons of the year allowed
finding the permanent seasonal average concentrations of phosphates in the water runoff entering the rivers from each of the specified landscapes. The seasonal and long-term (1951-2003) dynamics of PO43— runoff was evaluated for the landscapes in all 34 river basins. Large areas of arable lands contribute significantly to the runoff of phosphates throughout the year that is associated with the use of fertilizers. The combined models of spatial generalization
and normalization of mean monthly temperature and monthly precipitation, water and hydrochemical runoff allow us to estimate seasonal and long-term dynamics of phosphates runoff for all river basins of the Altai-Sayan mountain country. All the models are applicable even in case of limited experimental hy-drometeorological and hydrochemical data on river basins.
Кирста Юрий Богданович, д.б.н., проф., гл. н.с. Институт водных и экологических проблем СО РАН, г. Барнаул. Тел.: (3852) 666-457. E-mail: kirsta@iwep.ru.
Пузанов Александр Васильевич, д.б.н., проф., ВрИО директора, Институт водных и экологических проблем СО РАН, г. Барнаул. Тел.: (3852) 666-460. E-mail: iwep@iwep.ru.
Kirsta Yuriy Bogdanovich, Dr. Bio. Sci., Prof., Chief Staff Scientist, Institute for Water and Environmental Problems, Sib. Branch of Rus. Acad. of Sci., Barnaul. Ph.: (3852) 666-457. E-mail: kirsta@iwep.ru. Puzanov Aleksandr Vasilyevich, Dr. Bio. Sci., Prof., Acting Director, Institute for Water and Environmental Problems, Sib. Branch of Rus. Acad. of Sci., Barnaul. Ph.: (3852) 666-460. E-mail: iwep@iwep.ru.
Введение
Решение проблем рационального использования водных ресурсов, в том числе для орошения, требует разработки адекватных математических моделей, количественно характеризующих качество этих ресурсов. Гидрохимический сток рек является одной из важнейших характеристик качества их вод. Анализ закономерностей формирования гидрохимического стока особенно трудоемок для горных территорий из-за их сложного гидрогеологического строения, разнообразного почвенно-растительного покрова, сложной пространственно-временной структуры климатических полей, трудоемкости и больших материальных затрат на подобные исследования. Крайне затруднена количественная оценка воздействия на сток разнообразных факторов среды.
Цель исследования — разработка модели гидрохимического стока на примере растворимой формы фосфатов (анионов Р043—). Нами используется системный подход и хорошо зарекомендовавший себя метод системно-аналитического моделирования (САМ) с привлечением ГИС-технологий [1, 2].
Сток Р04 зависит от идущих в речных бассейнах геохимических процессов, внесения на сельскохозяйственные поля фосфорных удобрений и возможного «сухого» выпадения фосфатов из атмосферы на поверхность водосбора [3]. Важной особенностью поведения фосфатов в почвах являются их слабая растворимость и иммобилизация (переход в нерастворимую форму) [4].
Объекты и методы
Рассматриваемая в данном исследовании территория Алтае-Саянской горной страны представляет собой часть мирового водораздела между гумидной областью Северного Ледовитого океана и аридной бессточной об-
ластью Центральной Азии. Климат резко кон-тинентален, с холодной зимой и холодным летом. У рек наибольший водный сток наблюдается в теплое время года и составляет до 80-90% годового.
Для выполнения САМ гидрохимического (анионов Р043—) стока средних и малых рек страны взяты 34 бассейна (рис. 1), для которых ранее была разработана имитационная модель водного стока [5]. Динамика водного и гидрохимического стоков рассчитывалась для периода 1951-2003 гг. по 4 характерным гидрологическим периодам/сезонам: первый (зимняя межень, ХН-Ш месяцы), второй (весенне-летнее половодье, третий (летняя межень, VII-VШ), четвертый (осенняя межень с возможными паводками при сильных дождях, 1Х-Х1). Учитывалось, что в каждом сезоне должно иметься одно или более экспериментальных измерений концентраций Р04. После их усреднения были получены по четыре среднесезонных концентрации для каждого года наблюдений. В среднем по всем 34 бассейнам концентрации Р04 составили 0,019; 0,029; 0,018; 0,017 мг/л для 1-, 2-, 3-, 4-го сезонов соответственно. Эти значения были использованы для нормировки наблюдаемых значений концентраций и перехода, тем самым к безразмерным единицам измерений.
Температуры воздуха и осадки характеризовались нами через доли/проценты от их среднемноголетних месячных значений. Для таких нормированных температур и осадков с помощью разработанного метода обобщения получена единая для Алтае-Саянской горной страны их многолетняя помесячная динамика [6]. Показано, что эта динамика является адекватной площадной характеристикой изменения метеорологических полей всей страны.
84*0,0"Е БВ'О'О'Е
84°0'0"Е 88°0ХГ¥
Рис. 1. Карта-схема расположения 34 модельных речньх бассейнов Алтае-Саянской горной страны
Таблица
Найденные в ходе САМ среднесезонные концентрации фосфатов (анионов РО/~) в водных стоках, поступающих с ландшафтов Алтае-Саянской горной страны
Группы геосистем (ландшафтов) Концентрации PO4 по 4 гидрологическим сезонам года (XII-III, IV-VI, VII-VIII, IX-XI мес.), мг/л
1 2 3 4
1. Гляциально-нивальные высокогорья (характеризуются параметром a( в (2)) 0,019 0 0,014 0
2. Гольцово-альпинотипные высокогорья и среднегорья, псевдогольцовые низкогорья (a2) 0 0 0 0
3. Тундрово-степные и крио-фитно-степные высокогорья (a3) 0,022 0 0,048 0
4. Лесные высокогорья, среднегорья и низкогорья (a4) 0 0,008 0,003 0,008
5. Экспозиционно-лесостепные и степные высокогорья и среднегорья (a5) 0,024 0,011 0,002 0
6. Лесостепные, степные низкогорья и предгорья (a6) 0 0,137 0,012 0,015
7. Межгорные котловины с различными вариантами степей и лесостепей (a7) 0 0,147 0,074 0,078
8. Степные и лесостепные подгорные и возвышенные аккумулятивные равнины (ag) 0,095 0,089 0,037 0,023
9. Недренируемые слабо-проточные интразональные и ин-трапоясные ландшафты (a9) 0 0 0 0
10. Долины горных рек (a10) 0 0,074 0,066 0,020
11. Долины равнинных рек (an) 0 0,111 0,009 0,008
12. Лесные возвышенные и подгорные равнины (a12) 0 0 0 0
13. Аквальные ландшафты (a13) 0,006 0 0 0
Подземный сток (b в (2)) 0,018 0,016 0,006 0,007
Пашня (d в (2)) 0,038 0,059 0,036 0,035
Для учета ландшафтной структуры речных бассейнов горной страны выделено 12 типологических групп геосистем и отдельно 13-я для аквальных ландшафтов, имеющих незначительную площадь (табл.) [7]. Для всех групп в каждом речном бассейне по модели водного стока горных рек были рассчитаны формируемые ими вклады в этот сток [5].
В целом, база данных для разработки модели стока фосфатов включала:
• 1240 среднесезонных концентраций фосфатов (анионов Р043—) в разные годы в водных стоках 34 речных бассейнов, то есть 1240/4=310 величин для каждого сезона;
• расчетный вклад в водный сток от каждой группы геосистем в каждом речном бассейне по каждому году гидрохимических наблюдений (3120, 3480, 3180, 3350 значений водных стоков для 1-, 2-, 3-, 4-го гидрологических сезонов соответственно);
• обобщенные по территории Алтае-Саянской горной страны нормированные месячные осадки и среднемесячные температуры воздуха за период 1951-2003 гг. (по 636 значений);
• площадь и средняя высота каждого из 34 речных бассейнов, высота его замыкающего створа, длина речного русла (от истока до створа);
• доля площади пашни (при ее наличии) в каждом речном бассейне.
Результаты и их обсуждение
Основу САМ составляет определение структуры и параметров имитационных моделей сложноорганизованных природных систем путем решения оптимизационными методами обратной задачи [1, 2, 5]. Путем последовательной проверки различных физически и гидрохимически непротиворечивых зависимостей, описывающих формирование стока фосфатов (Р04) под воздействием факторов среды, определяется система уравнений, дающая наименьшую квадратичную невязку (сумму квадратов разностей) между рассчитанными и наблюдаемыми концентрациями вещества в воде. Построение модели стока Р04 выполняется в известной среде программирования МА^АВ.
Для описания нетривиальных зависимостей процессов от факторов среды в САМ применяется универсальная функция Н, задаваемая выражением:
Н (X1, X 2, Yl, Y 2, Zl, Z 2, X)
71 + 21- (X - XI), 7 2 - 71
-(X - XI) + 71,
X 2 - XI
если X < XI
XI < X < X 2,
если
XI * X 2
72 + 22 - (X - X2), если X > X2 (^
где Х1, Х2, У1, У2, Z1, Z2 — подбираемые параметры;
X — какая-либо переменная модели.
Функция Н является непрерывной кусочно-линейной функцией из трех произвольных линейных фрагментов и позволяет аппроксимировать широкий спектр различных зависимостей между переменными и факторами среды путем изменения значений своих параметров (рис. 2).
Y=H
(Х2^2)
(Х1,Л)
X
Рис. 2. Непрерывная кусочно-линейная функция Н(Х1,Х2,У1ГУ2Г21Г22ГХ) из трех
линейных фрагментов с произвольно меняемыми параметрами (уравнение (1)
В результате САМ стока фосфатов получены следующие универсальные уравнения, описывающие этот сток:
для первого гидрологического сезона —
Е ^Н (С1, ¿1,1,1, с2, съ, Р) сток = к , (2а)
Н (с4, с4,1,1, с5, с6, К') + Ьц' для второго, третьего и четвертого гидрологических сезонов —
Е°к&Н (с,, ¿1,1,1, с2, Сз, Р)Н
сток = к , (2б)
(с4, с4,1,1, с5, сб, К ) + Ьд'+4>[&& где Р — обобщенные по территории Алтае-Саян нормированные осадки [6] за 1Х-Х1 месяцы предшествующего года для 1-го сезона или за VII-VШ, 1Х-Х1 месяцы для 2-,
3-, 4-го сезонов соответственно;
ак — параметры, отвечающие постоянной среднесезонной концентрации вещества
(Р04) в расчетном водном стоке Qk , формируемом к-й группой геосистем за счет осадков Р, к=Ж3;
Ь — параметр, сопоставляемый с постоянной среднесезонной концентрацией вещества в расчетном приходящем (или уходящем) среднесезонном подземном водном
стоке ц' , который формируется в бассейне i почвенно-грунтовыми водами и водами зон трещиноватых пород;
К — средний поперечный уклон бассейна i, рассчитываемый по картографическим
данным как тангенс угла наклона склонов относительно горизонтали (рис. 3) [8]; Н — кусочно-линейная функция (1); с1ч6 — параметры, отражающие влияние
на сток вещества осадков Р и уклона К ;
^ — параметр, характеризующий увеличение концентрации вещества от каждого
процента площади Б1 в расчетном водном
стоке О1 ;
Б' — относительная площадь пахотных земель (в долях/процентах от площади бассейна /).
Вместе с моделью пространственного обобщения и нормировки среднемесячных температур и месячных осадков [6], моделью водного стока [5] уравнения (2а), (2б) составляют полную имитационную модель климата, водного и гидрохимического (Р04) стоков горных рек. Все параметры уравнений (2) определены в ходе САМ через решение обратной задачи по ежегодно наблюдаемым среднесезонным стокам вещества, найденным как О'С'. Величина О' характеризует среднесезонный водный сток для замыкающего створа бассейна с номером /=1^34 в текущем году, который рассчитывается по модели водного стока и нормируется на свое среднемноголетнее наблюдаемое значение в
бассейне / [5]. С' представляет собой наблюдаемые концентрации вещества в речном стоке для бассейна /, нормированные на их среднюю по 34 бассейнам величину за конкретный сезон. Отметим, что в САМ одновременно выполняются идентификация и верификации разрабатываемой модели (подробнее см. [2, 5]).
Значение поперечного уклона К в (2) рассчитывалось как тангенс угла между поверхностью бассейна и горизонталью (рис. 3) по формуле:
К' =■
h
уклона К речных бассейнов. В (2) отражена также роль пашни. Вклад Ь - ц1 учитывает приток или отток вещества при положительном или отрицательном ц' , рассчитываемом
по модели водного стока. Значение ц' отражает вклад осадков предыдущего и текущего сезонов, а также обмен влагой с поч-венно-грунтовыми водами и водами зон трещиноватых пород.
1/2 L
(средняя высота бассейна ') - (высота створа) 1/2(площадь бассейна ')/(длина речного русла)
Такой расчет К формально исключал влияние геометрического уклона речного русла, который не должен влиять на гидрохимический сток с ландшафтов.
В правой части (2) суммируются вклады в сезонный сток Р04 от каждой группы геосистем, обеспечиваемые поверхностным, внут-рипочвенным и подземным водными стоками. Через Н (с1, с1, 1,1, с2, с3, Р) учитывается
влияние на вклад к-й группы геосистем осадков Р текущего гидрологического сезона, а
через Н(с4, с4,1,1, с5, сб, К ) — поперечного
I
Рис. 3. Поперечный разрез речного бассейна и схема определения его
среднего поперечного уклона К' по средним высоте (^ и ширине !)
В ходе САМ через решение обратной задачи по стоку Р04 для каждого гидрологического сезона найдены все параметры а, Ь, с, ^ в (2) (табл.). Значения параметров ак , где к=1^13, характеризуют постоянные среднесезонные концентрации Р04 в водных стоках от каждой из 13 групп геосистем Алтае-Саянской горной страны. Отметим, что какие-либо экспериментальные измерения концентраций Р04 для отдельных ландшафтов вообще отсутствовали, и ак характеризуют их теоретические значения, найденные с помощью САМ.
Из данных таблицы следует, что наиболее значительные концентрации фосфатов в ландшафтных стоках отвечают лесостепным, степным низкогорьям и предгорьям (а6), межгорным котловинам (а7), равнинам (а8) и долинам рек (а10, ап). Пахотные земли при значительных площадях дают заметный вклад в сток Р04 на протяжении всего года вследствие использования минеральных удобрений при возделывании сельскохозяйственных культур.
На рисунке 4 приведена типичная зависимость стока фосфатов (г/сек.) от осадков и поперечного уклона речных бассейнов для четырех гидрологических сезонов года. Отметим, что для зимней межени (рис. 4, А) приведены осадки за осенний сезон, так как зимой осадки остаются на поверхности в виде снега и не участвуют в формировании водного и гидрохимического стоков зимней межени. Мы видим, что сток
Ь
PO4 в целом падает при уменьшении осадков Р. Эта зависимость очевидна, поскольку меньшие осадки вымоют из почв и меньшее количество растворимых фосфатов (PO4). Обращает на себя внимание заметное падение стока PO4 с уклоном К в весенне-летнее половодье (рис. 4, Б) и такое же падение при больших осадках в другие сезоны года. Такая зависимость может объясняться преобладанием в таких случаях поверхностного/подповерхностного водного стока над грунтовым и подземным. Первый сток вымы-
вает фосфаты из верхнего 10 мм слоя почвы, где они преимущественно находятся [9,
10]. С уменьшением уклона К этот сток будет сокращаться, так как все большая доля поступившей в почву влаги будет просачиваться вглубь и уносить с собой часть растворенных фосфатов PO4 с их последующей иммобилизацией в нижележащих слоях. Поэтому уменьшение К вызовет и сокращение общего стока PO4.
Рис. 4. Зависимость стока фосфатов (Р04, г/с) от гипотетически разного поперечного уклона речного бассейна и нормированных на свое среднемноголетнее значение осадков для верховья р. Катунь (створ с. Тюнгур): А — зимняя межень (осадки за 1Х-Х1 месяцы предшествующего года); Б — весенне-летнее половодье (осадки за /У-У/месяцы текущего года); В — летняя межень (осадки за УН-УН/месяцы); Г — осенняя межень (осадки за 1Х-Х1 месяцы)
Выводы
На примере рек Алтае-Саянской горной страны выполнен системный анализ формирования стока фосфатов (анионов PO4). Проанализировано влияние на сток месячных осадков, ландшафтной структуры речных бассейнов (площадей ландшафтов), среднего поперечного уклона поверхности бассейнов, площади пашни.
С помощью метода системно-аналитического моделирования разработана универсальная имитационная модель сезонной и многолетней динамики стока PO4. Установле-
ны количественные зависимости последнего от меняющихся факторов среды.
Объединенные вместе модели климата (пространственного обобщения и нормировки среднемесячных температур и месячных осадков) [6], водного стока [5] и гидрохимического стока (уравнения (2)) позволяют рассчитать сезонную и многолетнюю динамику стока PO4 для любых речных бассейнов Ал-тае-Саянской горной страны. Для этого необходима лишь картографическая информация о ландшафтной структуре исследуемого бассейна, данные о среднемесячных температу-
рах воздуха и месячных осадках, а также 1-2-летние наблюдения за водным стоком с целью расчета коэффициента перехода от его нормированных значений в модели водного стока к измерению в м3/с.
Библиографический список
1. Kirsta Yu.B. System-analytical modelling — Part I: General principles and theoretically best accuracies of ecological models. Soil-moisture exchange in agroecosystems // Ecol. Modelling. — 2006. — Vol. 191. — P. 315-330.
2. Кирста Ю.Б., Кирста Б.Ю. Информационно-физический закон построения эволюционных систем. Системно-аналитическое моделирование экосистем. — Барнаул: Изд-во Алт. гос. ун-та, 2014. — 283 с.
3. Савичев О.Г., Иванов А.О. Атмосферные выпадения в бассейне Средней Оби и их влияние на гидрохимический сток рек // Известия РАН. — Сер. геогр. — 2010. — № 1. — С. 63-70.
4. Fageria N. K. The use of nutrients in crop plants. — CRC Press, Taylor and Franзis Group, 2009. — 448 p.
5. Кирста Ю.Б., Пузанов А.В., Ловц-кая О.В., Лубенец Л.Ф., Кузняк Я.Э., Пахот-нова А.Ю. Имитационная математическая модель стока средних и малых рек для горных территорий // Известия Самарского научного центра РАН. — 2012. — Т. 14. — № 1 (9). — С. 2334-2342.
6. Кирста Ю.Б. Пространственное обобщение климатических характеристик для горных территорий // Мир науки, культуры, образования. — 2011. — № 3(28). — С. 330-337.
7. Кирста Ю.Б., Лубенец Л.Ф., Черных Д.В. Типизация ландшафтов для оценки речного стока в Алтае-Саянской горной стране // Устойчивое развитие горных территорий. — 2011. — № 2 (8). — С. 51-56.
8. Кирста Ю.Б., Пузанов А.В., Ловц-кая О.В., Лубенец Л.Ф. Универсальная математическая модель стока взвешенных веществ для бассейнов горных рек // Устойчивое развитие горных территорий. — 2012. — № 3-4 (13-14). — С. 46-53.
9. Neitsch S.L., Arnold J.G., Kiniry J.R., Williams J.R. Soil and Water Assessment Tool. Theoretical Documentation. Version 2009. — Texas Water Resources Institute, 2011. Available electronically from http://swat.tamu.edu/media/99192/swat200 9-theory.pdf
10. Бутько А.А., Родькин О.И. Прогнозирование поступления фосфора в водные объекты с территории агроландшафтов / / Экологический вестник: научно-практический журнал. — 2012. — № 2(20). — С. 76-83.
References
1. Kirsta Yu.B. System-analytical modelling — Part I: General principles and theoretically best accuracies of ecological models. Soil-moisture exchange in agroecosystems // Ecol. Modelling. — 2006. — Vol. 191. — P. 315-330.
2. Kirsta Yu.B., Kirsta B.Yu. Informatsionno-fizicheskii zakon postroeniya evolyutsionnykh sistem. Sistemno-analiticheskoe modelirovanie ekosistem. — Barnaul: Izd-vo Alt. gos. un-ta, 2014. — 283 s.
3. Savichev O.G., Ivanov A.O. Atmosfernye vypadeniya v basseine Srednei Obi i ikh vliyanie na gidrokhimicheskii stok rek // Izvestiya RAN. Ser. geogr. — 2010. — № 1. — S. 63-70.
4. Fageria N.K. The use of nutrients in crop plants. — CRC Press, Taylor and Francis Group, 2009. — 448 p.
5. Kirsta Yu.B., Puzanov A.V., Lovts-kaya O.V., Lubenets L.F., Kuznyak Ya.E., Pakhotnova A.Yu. Imitatsionnaya matematich-eskaya model' stoka srednikh i malykh rek dlya gornykh territorii // Izvestiya Samarskogo nauchnogo tsentra RAN. — 2012. — T. 14. — № 1 (9). — S. 2334-2342.
6. Kirsta Yu.B. Prostranstvennoe obob-shchenie klimaticheskikh kharakteristik dlya gornykh territorii // Mir nauki, kul'tury, obra-zovaniya. — 2011. — № 3 (28). — S. 330-337.
7. Kirsta Yu.B., Lubenets L.F., Chernykh D.V. Tipizatsiya landshaftov dlya otsenki rech-nogo stoka v Altae-Sayanskoi gornoi strane // Ustoichivoe razvitie gornykh territorii. — 2011. — № 2 (8). — S. 51-56.
8. Kirsta Yu.B., Puzanov A.V., Lovts-kaya O.V., Lubenets L.F. Universal'naya ma-tematicheskaya model' stoka vzveshennykh veshchestv dlya basseinov gornykh rek // Ustoichivoe razvitie gornykh territorii. — 2012. — № 3-4 (13-14). — S. 46-53.
9. Neitsch S.L., Arnold J.G., Kiniry J.R., Williams J.R. Soil and Water Assessment Tool. Theoretical Documentation. Version 2009. — Texas Water Resources Institute, 2011. (Available electronically from http://swat.tamu.edu/ media/99192/swat2009-theory.pdf.).
10. But'ko A.A., Rod'kin O.I. Prognoziro-vanie postupleniya fosfora v vodnye ob"ekty s territorii agrolandshaftov // Ekologicheskii vestnik: nauchno-prakticheskii zhurnal. — 2012. — № 2 (20). — S. 76-83.
+ + +
