Оценка выноса биогенных веществ с малых водосборных бассейнов Калининградского/Вислинского залива посредством численного моделирования

УДК 556.53

ОЦЕНКА ВЫНОСА БИОГЕННЫХ ВЕЩЕСТВ С МАЛЫХ ВОДОСБОРНЫХ

БАССЕЙНОВ КАЛИНИНГРАДСКОГО/ВИСЛИНСКОГО ЗАЛИВА ПОСРЕДСТВОМ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ

Д. А. Домнин, Б. В. Чубаренко

NUTRIENT REMOVAL ASSESSMENT FROM SMALL CATCHMENT AREAS OF THE VISTULA LAGOON USING NUMERICAL MODELING

D. A. Domnin, B. V. Chubarenko

Калининградский/Вислинский залив Балтийского моря преобразует биогенную нагрузку, поступающую с водосбора площадью 23,9 тыс. км2. Проанализированы два частных водосбора: р. Приморской (120 км , небольшие поселения и сельскохозяйственные угодья) и трансграничный между Россией и Польшей р. Банувки-Мамоновки (350 км , относительно крупные населённые пункты, предприятия пищевой промышленности и сельского хозяйства).

С помощью численного моделирования произведены оценки выноса с водосбора общего азота и общего фосфора и удержания этих элементов в его пределах. В качестве исходных данных использованы морфометрические характеристики водосбора, информация о структуре землепользования и точечных источниках сброса биогенных веществ, временные ряды данных о расходе воды, температуре воздуха и атмосферных осадках. Гидрологические параметры были получены с использованием модели HYPE, учитывающей испарение в атмосферу и инфильтрацию влаги в почву. Удержание и транспортировка биогенных веществ рассчитаны с помощью модели FyrisNP.

Вынос биогенных элементов с водосборов рек Банувки-Мамоновки (20082011 гг.) и Приморской (2010-2014 гг.) составляет 400 и 43 т общего азота и 25 и 6 т общего фосфора в год, а удержание на этих территориях оценивается в 20 и 8 % для азота и 31 и 2 % для фосфора соответственно. Определены вклады различных источников в итоговый вынос биогенных веществ. Наибольшее количество общего азота и фосфора в замыкающих створах приходит от пахотных земель (50-80 %), точечные источники составляют меньшую часть (5-30 %).

Показано, что оценка нагрузки в пропорции к площади земель различного типа является слишком неопределенной, и только путем моделирования можно учесть все аспекты взаимодействий.

математическое моделирование, водосборный бассейн, биогенная нагрузка, Калининградский залив, Вислинский залив

The Vistula Lagoon of the Baltic Sea converts the nutrient load coming from the

о

catchment area of 23,9 thousand km . Two water catchment areas have been analyzed: the catchment of the Primorskaya River (120 km2, small settlements and agricultural lands) and the cross-border catchement of the Banuvki-Mamonovka River between

2

Russia and Poland (350 km , relatively large populated areas, food and agriculture enterprises).

Using numerical modeling, estimates have been made of the removal from the water catchment areas of total nitrogen and total phosphorus and their retention within the catchment areas. Morphometric characteristics of the catchment area, information on the structure of land use and point sources of nutrient discharge, time series of data on water consumption, air temperature and precipitation have been used as the initial data. Hydrological parameters have been obtained using the HYPE model, taking into account evaporation and infiltration into the soil. The retention and transportation of nutrients are calculated using the FyrisNP model.

Emission of nutrients from catchments of the Banowka-Mamonovka River (20082011) and the Primorskya River (2010-2014) is 400 and 43 tons of total nitrogen and 25 and 6 tons of total phosphorus per year, the retention in these catchment areas are estimated as 20 and 8 % for nitrogen and 31 and 2 % for phosphorus, respectively.

Source apportionment has been made for the nutrient load discharging from both catchments to the Vistula Lagoon. The greatest amount of nutrients in final discharge is coming from the arable land (50-80 %), point sources constitute a smaller proportion (530 %).

It has been shown that estimations of emission based on the proportion of land use types brings too big uncertainty, and only modeling may consider all nonlinear aspects.

mathematical modeling, catchment area, nutrient load, Vistula Lagoon

ВВЕДЕНИЕ

Эвтрофикация, вызванная переизбытком биогенных элементов, является одной из ключевых проблем качества вод Балтийского моря [1]. В 2007 г. был принят План действий ХЕЛКОМ [2] по защите окружающей среды Балтийского моря, определяющий снижение биогенной нагрузки со стороны окружающих его стран. Согласно ему в Калининградской области необходимо уменьшение сбросов на 25 % по азоту и на 50 % по фосфору от существующего уровня [3]. Учитывая разрозненность малых поселений и объектов сельского хозяйства, недостаточно распределять эту квоту по уменьшению нагрузки просто пропорционально площади водосбора, а необходима оценка (аналогично [4]) конкретных условий и роли отдельных водосборов в формировании биогенного стока, поступающего в море. В качестве тестовых были выбраны водосборы малых рек Банувки-Мамоновки и Приморской, на которых имеются точки мониторинга. Эти водосборы расположены в бассейне Калининградского/Вислинского залива1, их общая площадь составляет 3% от площади водосбора залива.

1 Залив по-разному именуется в различных источниках. Немецкое название всего залива (838 км2) Frisches Haff, в польской литературе залив называется Wislany Zalew. В англоязычных изданиях чаще всего употребляется этноним «the Vistula Lagoon». В российских официальных изданиях

северная (российская) часть залива (56.2% площади акватории) называется Калининградским заливом, южная (польская) часть - Вислинским заливом. Здесь и далее, следуя [5], используется название «Калининградский/Вислинский залив», хотя авторы ранее [6] придерживались терминологии [7] и именовали весь залив Вислинским.

ОБЪЕКТЫ, МАТЕРИАЛЫ И МЕТОДЫ Водосборный бассейн р. Банувки-Мамоновки является трансграничным (рис. 1). Река начинается в Польше под названием Банувка, впадает в Калинин-градский/Вислинский залив Балтийского моря на территории Калининградской

области России под названием Мамоновка [7]. Площадь её водосбора составляет 2 2 и ц 350 км , из них 140 км находятся в юго-западной части Калининградской области, а 210 км2 - в северной части Варминско-Мазурского воеводства Польши [8].

Рис. 1. Расположение водосборных бассейнов рек Банувки-Мамоновки (1) и Приморской (2), их замыкающих створов (3, 4, соответственно) в пределах

водосбора Калининградского/Вислинского залива (по материалам [8]) Fig. 1. Locations of catchment basins of the Banowka-Mamonovka and Pri-morskaya Rivers within the Vistula Lagoon catchment

Водосбор (120 км ) р. Приморской (длина без притоков 30 км) расположен на Самбийском п-ове в западной части Калининградской области. Река впадает в Приморскую бухту Калининградского/Вислинского залива [8].

Для моделирования речного стока использовался гидрологический модуль модели HYPE [9] Шведского метеорологического и гидрологического института (SMHI), а для оценки стока биогенных веществ и определения вклада различных источников - модель FyrisNP [10] Шведского университета сельскохозяйственных наук (SLU).

Модель HYPE воспроизводит временную зависимость с шагом в одни сутки расхода воды в замыкающих створах речных водосборных бассейнов с учётом испарения и инфильтрации в почву по данным о морфометрии рельефа, структуре землепользования и почв, временным рядам атмосферных осадков и температуры воздуха [9].

Модель FyrisNP по данным о поступлении валового азота и фосфора оценивает их вынос с водосборного бассейна с учётом удержания. Временной шаг

модели составляет один месяц. Удержание биогенных веществ за счёт осаждения в реках и озёрах, усвоения растениями и денитрификации рассчитывается в зависимости от температуры воды, концентрации питательных веществ, слоя стока воды, площади водных объектов [10].

Для модельных инсталляций HYPE использовались данные по температуре воздуха и атмосферным осадкам со станции Мамоново (водосбор р. Банувки-Мамоновки) и со станции Калининград (водосбор р. Приморской) [11, 12]). Ка-либрация гидрологической модели HYPE проводилась по ежедневным данным для речного стока, а модели FyrisNP - по данным концентрации общего азота и общего фосфора в воде в ключевые гидрологические сезоны. Для р. Банувки-Мамоновки использованы данные станции Мамоново (7 км выше устья), для р. Приморской - станции Приморск (1,5 км выше устья). Период моделирования для р. Банувки-Мамоновки - 2007-2011 гг., для р. Приморской - 2010-2014 гг.

В качестве точечных источников загрязнения принимались только населённые пункты, загрязнения от предприятий-водопользователей не учитывались в силу конфиденциальности информации. Сточные воды от населенных пунктов попадают в поверхностные водные объекты слабо очищенными [13]. Для оценки нагрузки использовались данные региональной статистики о количестве населения в пределах муниципальных образований и информация [8]. Предполагалось, что нагрузка от каждого жителя составляет 5,5 г азота и 1,2 г фосфора в день, а в соответствии с оценкой [14] в конечном итоге сточные воды дают около 40 % общего азота и 50 % общего фосфора от этого количества. Было принято, что 20 % населения не подключено к канализации. Для этой группы степень удержания через инфильтрацию в почву составляла 76 % для азота и 88 % для фосфора [14].

Поскольку для периодов моделирования производительность сельскохозяйственного сектора в российских частях водосборов была низкой, большинство земель в моделях описаны как временно неиспользуемые земли («открытые»). Для польской части водосбора р. Банувки-Мамоновки структура землепользования определена по методике [15] на 2005 г. Для расчётов биогенной нагрузки, поступающей с различных типов земель, использованы данные об удельных концентрациях азота и фосфора в водах, стекающих с разных типов земель, принятые для юго-восточной Швеции [14].

МОДЕЛЬНЫЕ ИНСТАЛЛЯЦИИ

Водосбор р. Банувки-Мамоновки был разделён на десять частных водосборов - по пять в польской и российской частях (наименьшая площадь - 9,

2 u

наибольшая - 78 км ), а р. Приморской - на 15 частных водосборов. Площадь наименьшего составила 1,5, а наибольшего - 19 км2.

Исследование чувствительности гидрологической модели показало, что решение максимально зависит от параметров, влияющих на испарение с водосборной площади.

Для р. Банувки-Мамоновки коэффициент корреляции между измеренными и рассчитанными значениями расхода в контрольном створе для калибрационного

о

периода (2008-2010 гг.) составил 0,93, средний измеренный расход - 4,3 м /с, а модельный - 4,4 м /с. (табл. 1). Для р. Приморской коэффициент корреляции равнялся 0,88 для калибрационного периода 2010-2013 гг., средний расход воды по

данным контактных измерений составил 1,3 м3/с, по данным модельного расчёта -1,2 м3/с (табл. 2).

Модель FyrisNP имеет два эмпирических калибрационных параметра (со и кув), которые определяют, соответственно, удержание и смыв биогенных элементов в зависимости от температуры (в диапазоне 0-20°С) и расхода воды. При расчётах выноса биогенных веществ с территории исследуемых водосборов критерием верности вычислений было совпадение средних значений концентрации за моделируемый период (табл. 3, рис. 2, 3). Такой подход обусловлен тем, что измеренные данные концентрации азота и фосфора являются разовыми, а при моделировании используется осреднённый месячный шаг.

Тот факт, что не удалось получить единого набора калибровочных параметров для обоих водосборных бассейнов, означает, что, несмотря на их сходство, неопределенности, связанные с пространственными данными, осадками и данными о биогенной нагрузке в водосборе, всё-таки высоки. А значит, модели выноса биогенных веществ должны быть специально настроены для каждого водосбора.

Таблица 1. Средние, максимальные и минимальные значения расхода воды р. Банувки-Мамоновки, полученные путём измерений и моделирования для калибрационного (2008-2010 гг.) и верификационного (2011 г.) периодов

Table 1. Mean, maximum and minimum water discharge for the Banowka-Mamonovka river, obtained by measurements and simulation for calibration (2008-2010) and verification (2011) periods_

Временные серии Сред- Среднее, м3/с Макси- Мини-

нее, маль- маль-

км3/год ное, м3/с ное, м3/с

Результаты калибрации модели

Данные измерений 0,14 4,3 12,7 1,2

Результат моделирования, 2008-2010 гг. 0,14 4,4 14,0 0,8

Результаты верификации модели

Данные измерений 0,14 4,3 6,2 2,3

Результат моделирования, 2011 г. 0,13 4,0 8,2 1,2

Таблица 2. Средние, максимальные и минимальные значения расхода воды р. Приморской, полученные путём измерений и моделирования для калибрационного (2010-2013 гг.) и верификационного (2014 г.) периодов

Table 2. Mean, maximum and minimum water discharge for the Primorskaya river, obtained by measurements and simulation for calibration (2010-2013) and verification (2014) periods_

Временные серии Среднее, км3/год Среднее, м3/с Максимальное, м3/с Минимальное, м3/с

1 2 3 4 5

Результаты калибрации модели

Данные измерений 0,04 1,3 4,1 0,1

Результат моделирования, 20102013 гг. 0,04 1,2 8,5 0,2

Окончание табл. 2

1 2 3 4 5

Результаты верификации модели

Данные измерений 0,6 0,02 2,1 0,1

Результат моделирования, 2014 г. 0,6 0,02 4,0 0,1

Таблица 3. Средние концентрации общего азота и фосфора в воде рек Банувки-Мамоновки и Приморской, полученные путём измерений и моделирования Table 3. Mean concentration of total nitrogen and phosphorus in the Banowka-Mamonovka

and Primorskaya rivers according to measurements and simulation

Река (период моделирования) Общий азот, мг/л Общий фосфор, мг/л

Изме- По Изме- По

ренный модели ренный модели

Банувка-Мамоновка (2008-2011гг.) 2,22 2,33 0,14 0,14

Приморская (2010-2014 гг.) 1,24 1,15 0,26 0,22

a) б)

Рис. 2. Вариации концентраций общего азота (а) и общего фосфора (б) в р. Банувке-Мамоновке для периода 2008-2011 гг. по данным измерений и моделирования Fig. 2. Concentration variations of total nitrogen (a) and total phosphorus (b) in the Banowka-Mamonovka river for the period of 2008-2011 according to measurements

and modeling

a) б)

Рис. 3. Вариации концентраций общего азота (а) и общего фосфора (б) в р. Приморской для периода 2010-2014 гг. по данным измерений и моделирования

Fig. 3. Concentration variations of total nitrogen (a) and total phosphorus (b) in the Primorskaya river for the period of 2010-2014 according to measurements and modeling

ОБСУЖДЕНИЕ

Среднегодовой вынос общего азота с водосборного бассейна р. Банувки-Мамоновки в Калининградский/Вислинский залив составил 400 т/год за период 2008-2011 гг., а общего фосфора - 25 т/год. При этом удержание в водосборном бассейне общего азота - 20, а общего фосфора - 31 %. Нагрузка от пахотных земель является наиболее значительной: 82 % - для общего азота и 61 % - для общего фосфора (рис. 4).

В среднем со всего водосборного бассейна р. Приморской удерживается 8 % общего азота и 2 % общего фосфора. Средняя годовая нагрузка в Калининградский/Вислинский залив за период 2010-2014 гг. составляла 43 т/год общего азота и 6,2 т/год общего фосфора. Наибольший вклад вносили пахотные земли, точечные источники и фоновая нагрузка (рис. 5).

Для р. Приморской, водосбор которой полностью находится в пределах Калининградской области, наиболее значимой является нагрузка от пахотных земель - 58 % для азота и 49 % для фосфора (рис. 6). Открытые территории занимают почти 60 % её водосбора, но дают только 23 % азота и 18 % фосфора. Также меньшую долю биогенов, по сравнению с занимаемой ими площадью, привносят леса (8 % азота, 1 % фосфора) и заболоченные территории (2 % азота, 0,1 % фосфора). При этом точечные источники выступают одним из наиболее значительных биогенных загрязнителей в водосборном бассейне. Количество поступающего от них азота составляет 9 % от общего количества, а фосфора - 32 %. Тем самым проверена несостоятельность подхода, когда вклад в вынос биогенов от пространственных источников оценивается пропорционально занимаемым ими площадям.

Источники поступления Роб

1

1-Леса

A. 7 (13%>  2 - Заболоченные

территории

МШк^ 3 - Пашни

4 - Пастбища

3

5 - Открытые территории

5(7%)/

6 - Основные точечные

источники

4(1%)  7 - Второстепенные

точечные источники

a) б)

Рис. 4. Распределение вклада различных источников в вынос с водосборного бассейна р. Банувки-Мамоновки общего азота (a) и общего фосфора (б) по результатам моделирования для периода 2008-2011 гг. Fig. 4. Sources apportionment of nitrogen (a) and phosphorus (b) in the Banowka-Mamonovka catchment area by the modeling results for the period 2008-2011

a) б)

Рис. 5. Распределение вклада различных источников в вынос с водосборного бассейна р. Приморской общего азота (a) и общего фосфора (б) по результатам

моделирования для периода 2010-2014 гг. Fig. 5. Sources apportionment of nitrogen (a) and phosphorus (b) in the Primorskaya catchment area by the modeling results for the period 2010-2014

Леса 60

40

Точечные источники 20 у Заболоченные территории

Открытые А fiV 'Ч - О-.Ч Площадь, % N общ, %

территории Пашни Р общ, %

Рис. 6. Соотношение долей вклада различных источников в вынос биогенов с водосбора р. Приморской и занимаемые ими площади (в %) Fig. 6. Fractions (in %) of nutrient sources in the nitrogen and phosphorus emission from the Primorskaya River catchment and shares of these sources in total catchment area

ВЫВОДЫ

Посредством численной модели, откалиброванной по данным концентраций биогенных веществ в реках Банувке-Мамоновке (2008-2011 гг.) и Приморской (2010-2014 гг.), определено, что вынос биогенных элементов с их водосборов составляет 400 и 43 т общего азота и 25 и 6 т общего фосфора в год, а удержание на этих территориях оценивается 20 и 8 % для азота и 31 и 2 % для фосфора соответственно.

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

Полученные результаты показали, что оценка нагрузки в соответствии с учётом пропорций земель различного типа использования является слишком неопределенной. Не существует альтернативы применению численных моделей (даже самых простых), поскольку при оценке выноса биогенных веществ с водосборных территорий должны быть приняты во внимание многие процессы (в том числе и нелинейные).

БЛАГОДАРНОСТИ Расчёты для водосбора р. Банувки-Мамоновки выполнены при поддержке проекта BaltHazAR «Создание потенциала в рамках экологического мониторинга для получения данных загрязнения из различных источников, например, для HELCOM PLCs» (2012 г.), для водосбора р. Приморской - гранта РФФИ-Бонус 14-05-91730 «Уменьшение биогенной нагрузки с сельскохозяйственных территорий в Балтийское море (Soils2Sea)» (2014-2016 гг.). Подготовка модели и сбор сопутствующих данных проводился в рамках темы № 0149-2014-0017 государственного задания АО ИО РАН.

СПИСОК ИСПОЛЬЗОВАННЫХ ЛИТЕРАТУРНЫХ ИСТОЧНИКОВ

1. State and Evolution of the Baltic Sea, 1952-2005: A Detailed 50-Year Survey of Meteorology and Climate, Physics, Chemistry, Biology, and Marine Environment / Ed. R. Feistel, G.Nausch, N. Wasmund. 2008. 704 p. ISBN: 9780471979685

2. HELCOM Baltic Sea Action Plan. Krakow, 2007, 101 p.

3. HELCOM Copenhagen Ministerial Declaration. Copenhagen, 2013, 19 p.

4. Hessea, C. Assessment of climate change impacts on water quantity and quality of the multi-river Vistula Lagoon catchment / C. Hessea, V. Krysanova, A. Stefa-nova, M. Bieleckab, D. Domnin // Hydrological Sciences Journal. - 2015. - vol. 60. -Issue 5. - P. 890-911.

5. Регион Калининградского/Вислинского залива: современное состояние и сценарий развития / под ред. В. Кушевски, Г. М. Федорова, Б. В. Чубаренко,

B.А. Гриценко. - Калининград: Изд-во БФУ им. И.Канта, 2014. - 216 с.

6. Домнин, Д. А. Трансграничные водосборы Юго-Восточной Балтики / Д. А. Домнин, Б. В. Чубаренко // География и природные ресурсы. - 2012. - № 3. -

C. 69-76.

7. Гидрометеорологические условия шельфовой зоны морей СССР [Hy-drometeorological conditions of the shelf zone of the seas of the USSR]. - Т. 1: Балтийское море. Вып. 3: Куршский и Вислинский заливы [The Curonian and the Vistula Lagoons] / под ред. Ф. С. Терзиева. - Ленинград: Гидрометеоиздат, 1985. - 72 с.

8. Domnin, D. Vistula Lagoon Catchment: Atlas of water use. / D. Domnin, B. Chubarenko, A. Lewandowski - Moscow: Exlibris Press, 2015. - 106 p. ISBN 978-59900699-4-7.

9. Donnelly, C., High-resolution, large-scale hydrological modelling tools for Europe. / C. Donnelly, J. Dahne, J. Rosberg , J. Stromqvist, W. Yang, B. Arheimer // Proc. of the Sixth World FRIEND Conference, Fez, Morocco, October 2010. IAHS Publ. 340, 2010.

10. Hansson, K. The FyrisNP model Version 3.1 - A tool for catchment-scale modelling of source apportioned gross and net transport of nitrogen and phosphorus in rivers / K. Hansson, M. Wallin, F. Djodjic, C. Orback. - Uppsala: Dept. of Aquatic Science and Assessment 2008. - 28 p.

11. RP5.RU: Архив погоды в Мамоново [Weather records in Mamonovo] [Электронный ресурс]. URL: http://rp5.ru/5193/ru

12. RP5.RU: Архив погоды в Калининграде [Электронный ресурс]. URL: http://rp5.ru/ Архив_погоды_в_Калининграде

13. Великанов, Н. Л. Калининградская область: особенности использования водных ресурсов / Н. Л. Великанов, Е. Д. Проскурнин. - Калининград: ФГУИПП «Янтарный сказ», 2003. - 128 c.

14. Swedish Environmental Protection Agency. What does domestic wastewater contain? // Swedish EPA Report. 1995. 4425.

15. Corine Land Cover technical guide - Addendum 2000 / Prepared by M. Bossard, J. Feranec, J. Otahel - Published European Environment Agency. - 2000. -104 p.

REFERENCES

1. State and Evolution of the Baltic Sea, 1952-2005: A Detailed 50-Year Survey of Meteorology and Climate, Physics, Chemistry, Biology, and Marine Environment. Ed. R. Feistel, G.Nausch, N. Wasmund. 2008. ISBN: 9780471979685.

2. HELCOM Baltic Sea Action Plan. Krakow, 2007, 101 p.

3. HELCOM Copenhagen Ministerial Declaration. Copenhagen, 2013, 19 p.

4. Hessea C., Krysanova V., Stefanova A., Bieleckab M., Domnin D. Assessment of climate change impacts on water quantity and quality of the multi-river Vistula Lagoon catchment. Hydrological Sciences Journal. 2015, vol. 60, iss. 5, pp. 890-911.

5. Region Kaliningraskogo/Vislinskogo zaliva: sovremennoe sostoyanie i scenrii razvitia [The Vistula Lagoon area: current state and development scenario]. Pod. red. Kushevski V., Fedorova G. M., Chubarenko B. V., Gricenko V. A. Kaaliningrad, BFU im. I. Kanta, 2014, 216 p.

6. Domnin D. A., Chubarenko B. V. Transgranichnye vodosbory Jugo-Vostochnoj Baltiki [Trans-border catchment areas of the south-eastern Baltics]. Geo-grafija iprirodnye resursy, 2012, no. 3, pp. 69-76.

7. Gidrometeorologicheskie uslovia shelfovoi zony morei SSSR [Hydrometeoro-logical conditions of the shelf zone of the seas of the USSR]. Leningrad, Hyddromete-oizddat, vol. 1, iss. 3, 1985, 72 p.

8. Domnin D., Chubarenko B., Lewandowski A. Vistula Lagoon Catchment: Atlas of water use. Moscow, Exlibris Press, 2015, 106 p.

9. Donnelly C., Dahne J., Rosberg J., Stromqvist J., Yang W., Arheimer B. High-resolution, large-scale hydrological modelling tools for Europe. Proc. of the Sixth World FRIEND Conference, Fez, Morocco, October 2010. IAHS Publ. 340, 2010.

10. Hansson K., Wallin M., Djodjic F., Orback C. The FyrisNP model Version 3.1 - A tool for catchment-scale modelling of source apportioned gross and net transport of nitrogen and phosphorus in rivers. Uppsala: Dept. of Aquatic Science and Assessment 2008, 28 p.

11. RP5. RU: Arhiv pogody v Mamonovo, available at: http://rp5.ru/5193/ru (Accessed 10.03.2017).

12. RP5.RU: Arhiv pogody v Kaliningrade [Weather records in Kaliningrad]. Avai-lable at: http://rp5.ru/ (Accessed 10 March 2017).

13. Velikanov N. L., Proskurnin E. D. Kaliningradskaja oblast': osobennosti ispol'zovanija vodnyh resursov [The Kaliningrad region: peculiarities of using water resources]. Kaliningrad, FGUIPP «Jantarnyj skaz», 2003, 128 p.

14. Swedish Environmental Protection Agency. What does domestic wastewater contain? Swedish EPA Report. 1995. 4425.

15. Corine Land Cover technical guide - Addendum 2000. Prepared by M. Bossard, J. Feranec, J. Otahel. Published European Environment Agency. 2000, 104 p.

ИНФОРМАЦИЯ ОБ АВТОРАХ

Домнин Дмитрий Александрович - Атлантическое отделение Института океанологии им. П. П. Ширшова РАН; научный сотрудник;

E-mail: dimanisha@rambler.ru

Domnin Dmitry Alexandrovich - Atlantic Branch of P. P. Shirshov's Institute of Oceanology of RAS; scientific researcher;

E-mail: dimanisha@rambler.ru

Чубаренко Борис Валентинович - Атлантическое отделение

Института океанологии им. П. П. Ширшова РАН; кандидат физико-математических наук; зав. лабораторией;

E-mail: chuboris@mail.ru

Chubarenko Boris Valentinovich - Atlantic Branch of P. P. Shirshov's Institute of Oceanology of RAS; PhD; head of the laboratory;

E-mail: chuboris@mail.ru