A study of long-term changes in Black Sea ecosystems based on data assimilation of remote measurements in a numerical model

Full Text

ВВЕДЕНИЕ

Черное море – практически замкнутый водоем; оно соединяется со Средиземным морем системой узких проливов. Уникальная в силу своей замкнутости экосистема Черного моря за последние десятилетия подверглась значительной трансформации. В работах [9, 19–21] на основе качественного анализа биогеохимических измерений описаны эти изменения, связанные как с климатическими условиями региона, так и с антропогенной нагрузкой, вызванной стоками рек. С 1960-х гг. экосистема Черного моря претерпела существенную перестройку по биомассе, таксономическому составу и структуре групп планктона. В верхнем слое моря произошли также изменения сезонных циклов основных компонентов экосистемы. В частности, естественный годовой цикл фитопланктона с весенним и осенним цветением сменился более сложным, характеризующимся несколькими максимумами. Эти трансформации структуры пищевой цепи сопровождались изменениями вертикальной гидрохимической структуры вод моря.

Наряду с анализом доступных данных биологических измерений достаточно давно начались работы по математическому моделированию динамики экосистемы Черного моря. Изменчивость морской экосистемы и ее реорганизация в течение разных периодов изучалась с помощью одномерных междисциплинарных моделей [11, 13, 15, 16, 18]. Одномерные модели позволили достичь общего качественного соответствия между модельными вертикальными распределениями биогеохимических параметров и вертикальными распределениями, полученными по данным наблюдений. Кроме того, они дали возможность достоверно описывать сезонную изменчивость параметров экосистемы. Однако особенности динамики Черного моря и главным образом интенсивная синоптическая изменчивость существенно влияют на биогеохимические процессы, один из них – перенос течениями богатых биогенными элементами вод северо-западного шельфа (СЗШ) моря. На изображениях Черного моря, полученных с помощью спутниковых цветовых сканеров, видно, что перенос продуктивных вод вдоль западного берега моря поддерживает повышенные концентрации хлорофилла вблизи Анатолийского побережья. Неоднородности распределения цвета моря, связанные с изменениями продуктивности вод бассейна, обусловлены синоптическими вихрями, меандрами основного Черноморского течения и филаментами, формирующимися в процессе трансфронтального переноса. Учесть все эти явления и их влияние на морскую экосистему можно только с помощью трехмерных моделей. Такие модели созданы и использовались для изучения поведения экосистемы Черного моря, в том числе ее длительной эволюции [2, 5, 10, 12, 14].

Появление регулярных измерений спутниковыми цветовыми сканерами внесло новый значительный вклад в изучение морских экосистем, позволив получать поверхностные распределения концентрации пигмента хлорофилла-а, что особенно важно для задач мониторинга окружающей среды. Вместе с тем существенный недостаток этих данных – большое количество их пропусков из-за наличия облачности, особенно в зимнее время. Один из перспективных способов решить эту проблему – ассимиляция данных измерений в математические модели морской среды. При этом использование трехмерных моделей позволяет, помимо заполнения пробелов в данных, обусловленных облачностью, экстраполировать эти данные с поверхности в глубину моря. Данная работа посвящена изучению изменчивости экосистемы Черного моря на основе ассимиляции спутниковых измерений в трехмерной математической модели.

МЕТОД ИССЛЕДОВАНИЯ

Трехмерная модель экосистемы Черного моря представляет собой междисциплинарную модель, в которую входят блок, описывающий циркуляцию вод Черного моря, и биогеохимический блок, описывающий взаимодействие между различными компонентами экосистемы на нижнем уровне пищевой цепи. В качестве модели циркуляции в Черном море использовалась модель, разработанная в Морском гидрофизическом институте [1], – z-уровенная модель, основанная на аппроксимации системы примитивных уравнений динамики океана и включающая в себя уравнения для горизонтальных компонент скорости течения, переноса–диффузии тепла и соли, неразрывности и гидростатики. Позже модель была переведена на сетку с шагом 4.8 км, что позволило помимо крупномасштабной циркуляции достаточно аккуратно описывать синоптические процессы (радиус деформации Кибеля–Россби для первой бароклинной моды в глубоководной части Черного моря равен ~25 км) [3, 14]. Используемый в работе вариант модели содержит 35 расчетных уровней, сгущающихся к морской поверхности, где дискретность равна 5 м. Вертикальная турбулентная диффузия и вязкость параметризуются с помощью коэффициентов, зависящих от глубины и времени.

Биогеохимический блок модели экосистемы Черного моря включает в себя 15 переменных состояния и описывает процессы в верхнем 200-метровом слое моря. Шаг по пространству и расчетные горизонты соответствуют модели циркуляции. В горизонтальной плоскости расчетная сетка имеет шаг 4.8 км, по вертикали модель экосистемы имеет 18 расчетных уровней в z-координатах, сгущающихся к поверхности моря. Биологическую часть представляют две группы фитопланктона (диатомовые и флагелляты), две размерные группы зоопланктона: микрозоопланктон (<0.2 мм) и мезозоопланктон (0.2–3.0 мм), медуза Aurelia aurita и гребневик Mnemiopsis leidyi. В пищевую цепь входят, кроме того, нефотосинтезирующий бактериопланктон, растворенная и взвешенная органика, а также всеядная динофлагеллята Noctiluca scintillans. Единственным биогенным элементом, лимитирующим рост фитопланктона, в данной модели рассматривается азот. Азотный цикл включает в себя три неорганических соединения: аммоний, нитраты и нитриты. Все перечисленные компоненты экосистемы в модели измеряются в мМ N/м³. В качестве отдельных переменных состояния в модель также входят растворенные кислород и сероводород. Временная эволюция компонентов экосистемы описывается уравнениями переноса–диффузии (1), включающими в правую часть слагаемые типа источник–сток:

$\frac{\partial F}{\partial t}+\frac{\partial (uF )}{\partial x}+\frac{\partial (vF )}{\partial y}+\frac{\partial \left(\right(w+wS)F )}{\partial z}=Kh{\nabla }^{2}F+\frac{\partial }{\partial z}\left(Kv\frac{\partial F}{\partial z}\right)+ R(F )$ (1)

R(F) описывает биогеохимические взаимодействия между переменными состояния F; w_s – скорость седиментации диатомовых водорослей и взвешенной органики (для остальных компонентов она равна нулю); u, v, w – компоненты скорости течений; K_h, K_v – коэффициенты горизонтальной и вертикальной турбулентной диффузии соответственно.

Члены правой части уравнения (1) типа источник–сток R(F) для биологической части модели аналогичны тем, которые использовались в одномерной модели [17].

R(P_d) = (1–χ)σ_dФf_p(T)P_d – G_S(Pd )H_S – G_l(P_d)H_l – G_n(P_d)Z_n – λ_dP_d; (2)

R(P_f) = (1–χ)σ_dФf_p(T )P_f – G_S(P_f )H_S – G_l (P_f)H_l – G_n(P_f )Z_n – λ_f P_f ; (3)

R(H_S) = γ_S[G_S (P_f )+G_S (P_d )+G_S (P_n)+G_S (B)]H_S – G_l(H_S)H_l – G_a(H_S)Z_a – G_m(H_S)Z_m – G_n(H_S )Z_n – μ_SH_S – λ_S H_S²; (4)

R(H_l ) = γ_l [G_l (P_f )+G_l (P_d )+G_l (P_n ) + G_l (H_S )+G_l (Z_n )]H_l – G_a (H_l )Z_a – G_m (H_l )Z_m – μ_l H_l – λ_l H_l² ; (5)

R(Z_n ) = γ_n [G_n (P_f )+G_n (P_d )+G_n (H_S )+G_n (P_n )]Z_n – G_i(Z_n )H_l – μ_nZ_n – λ_n Z_n ; (6)

R(Z_a ) = γ_a [G_a (H_S )+G_a (H_l )]Z_a – μ_a Z_a – λ_a Z_a ; (7)

R(Z_m ) = γ_m [G_m (H_S )+G_m (H_l )]Z_m – μ_m Z_m – λ_m Z_m ; (8)

R(B) = [G_b(D_n)+G_b(P_n)]B – G_S(B)H_S– μ_bB , (9)

где P_d и P_f – концентрации диатомовых и флагеллят групп фитопланктона; H_s, H_l – концентрации микро- и мезозоопланктона соответственно; Z_n – Noctiluca scintillans; Z_a – Aurelia aurita; Z_m – Mnemiopsis leidyi; В – концентрация бактериопланктона; G_i (G_i) – члены уравнений, описывающие питание i-м потребителем j-й компоненты F_j; f_p (T) – температурная лимитирующая функция; Ф – минимум от световой и азотной лимитирующих функций; χ – коэффициент экссудации; γ_i – эффективность ассимиляции; λ_i – коэффициент смертности; µ_i – коэффициент экскреции.

Соответствующие члены уравнения, выражающие концентрации взвешенного P_n и растворенного D_n органического вещества, нитрата [NO₃], нитрита [NO₂], аммония [NH₄], растворенного кислорода [O₂] и сероводорода [H₂ SH₂S] выбраны в соответствии с работой [18]:

R(P_n) = (1–γ_S [G_S (P_d )+G_S (P_f )+G_S (B)]H_S + (1 – γ_l [G_l (P_f )+G_S (P_d )+G_l (H_S )+G_l (Z_n )]H_l + (1 – γ_n [G_n (P_f )+G_n (P_d )+G_n (H_S )Z_n + (1 – γ_a [G_a (H_S )+G_a (H_l )]Z_a + (1 – γ_m [G_m (H_S )+G_m (H_l )]Z_m + [λ_d + λ_f + λ_SH_S² + λ_nZ_n + λ_aZ_a + λ_mZ_m ] – [γ_S G_S (P_n )H_S + γ_l G_l (P_n )H_l + γ_nG_n (P_n )Z_n + G_b (P_n )B] – ε[ f_n (O₂) + f_n (O₂)]P_n ; (10)

R(D_n ) = χФf_p (T) (σ_f P_f + σ_dP_d) + (1–κ)ε[ f_n (O₂) + f_n (O₂)]P_n ; (11)

R(NO₃ ) = ₂ f_n (O₂) [NO₂] – ₃ f_d (O₃) [NO₃] – β_n/β_tФf_p (T )(σ_f P_f + σ_d P_d ) – α₉ P₉ [H₂S][NO₃] ; (12)

R(NO₂) = ₁ f_n (O₂)[NH₄] – ₂ f_n (O₃)[NO₂] + ₃ f_d (O₂)[NO₃] – ₄ f_d (O₂)[NO₂] ; (13)

R(NH₄) = ε[ f_n (O₂) + f_n (O₂)]P_n – β_n/β_tФf_p (T )(σ_f P_f + σ_d P_d ) + μ_S H_S + μ_l H_l + μ_n Z_n + μ_a Z_a+ μ_m Z_m + μ_bB + ₁ f_n (O₂)[NH₄] ; (14)

R(O₂) = α₀Фf_p (T )(σ_f P_f + σ_d P_d ) – α_bεf_n (O₂)P_n – α₁₁ f_n (O₂)[NH₄] – α₂₂ f_n (O₂)[NO₂] – α₅₅ [O₂][H₂S] ; (15)

R(H₂S) = ₅[O₂][H₂S] – ₅[H₂S][NO₃] , (16)

где α_i – стехиометрический коэффициент,

_i – скорости реакций,  – коэффициент деминерализации детрита в аммоний, ε – скорость разложения детрита, β_α и β_n – лимитирующие функции для аммония и нитратов, β_τ = β_α + β_n, f_n O₂ и f_d O₂ – лимитирующие функции, контролирующие аэробное и анаэробное разложение детрита.

Поля скоростей течений, температура, соленость и коэффициенты турбулентного перемешивания, входящие в качестве коэффициентов в уравнение (1), обеспечиваются моделью циркуляции. В данном исследовании связь модели циркуляции с биогеохимической моделью осуществлялась в режиме off line, т.е. сначала рассчитывались необходимые гидрофизические поля, затем они использовались для моделирования эволюции параметров морской экосистемы. При расчете циркуляции Черного моря ассимилировались спутниковые данные по температуре воды на поверхности моря и аномалиям уровня моря, что повысило точность полученных гидрофизических полей [6].

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

Важная особенность настоящей работы – ассимиляция в модели экосистемы Черного моря данных измерений, полученных со спутниковых цветовых сканеров. На основе продуктов SeaWiFS, MODIS и MERIS по разработанному специально для Черного моря алгоритму подготовлены двухнедельные поля поверхностных концентраций хлорофилла-а [7, 8]. Далее значения концентраций хлорофилла-а были пересчитаны в концентрации фитопланктона, выраженные в мМ N/м³, по известным соотношениям между концентрацией хлорофилла-а и содержанием органического азота в морском фитопланктоне. Процедура ассимиляции основана на принципе релаксации полей, полученных по модели, к данным, подготовленным по изображениям сканеров. В правую часть уравнения (1) для двух компонент фитопланктона добавляется слагаемое типа источник–сток, зависящее от разности измеренной концентрации и полученной по модели:

1/τ (P_S – ϑ_i P_i(O)) , (17)

где P_s – данные спутниковых измерений, P_i (0) — концентрация одной из двух компонент фитопланктона на поверхности моря, ϑ_i — вес соответствующей компоненты в суммарной концентрации, τ – коэффициент релаксации. Процедура обновления слагаемого (17) проводится раз в сутки; соответственно, значение коэффициента релаксации равно суткам. Для того чтобы устранить появление или убыль в системе дополнительного азота, вызванного источником (17), аналогичный источник противоположного знака вводится в уравнение для нитратов. В отличие от (17), он работает на всех горизонтах, его амплитуда пропорциональна величине концентрации на данной глубине, а интегральное действие равно по величине суммарному действию источников в уравнениях для фитопланктона. Подробно алгоритм ассимиляции описан в работе [4].

РЕЗУЛЬТАТЫ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ

Расчет по модели экосистемы Черного моря проводился на срок 15 лет — с 1998 по 2012 г. Такой выбор объясняется тем, что данные, полученные со световых сканеров, стали доступны лишь с сентября 1997 г. Изменчивость средней за год биомассы фито- и зоопланктона в верхнем 50-метровом слое моря приведена на рис. 1. В целом биомасса фитопланктона на рассматриваемом промежутке времени уменьшается, соответственно падает биомасса зоопланктона, который питается фитопланктоном. В Черном море основное питательное вещество, обеспечивающее рост биомассы фитопланктона, — нитраты. Они поступают в поверхностный слой моря со стоком рек в основном на северо-западном шельфе и в глубоководной части моря из нижележащих слоев, имеющих максимальную концентрацию нитратов. Второй механизм насыщения фотической зоны моря нитратами действует в основном в зимнее время в период интенсивного ветрового и конвективного перемешивания. На рис. 2 показана зависимость от времени средней по площади бассейна концентрации нитратов на поверхности моря и в верхнем 50-метровом слое, который примерно соответствует толщине фотической зоны Черного моря.

 

Рис. 1. Изменчивость среднегодовых значений биомассы фитопланктона (а) и зоопланктона (б) в верхнем 50-метровом слое воды в Черном море, вычисленная по результатам моделирования за период 1998–2012 гг.

 

Рис. 2. Графики изменения средней поверхностной концентрации нитратов (а) и средней концентрации в верхнем 50-метровом слое воды в Черном море, полученные по результатам моделирования (б).

 

Графики изменения поверхностной концентрации нитратов имеют ярко выраженные пики в зимнее время, обусловленные поступлением нитратов из слоя максимума к поверхности в глубоководной части моря. Значительные величины поверхностной концентрации в летнее время обеспечиваются выносом рек. Ход средней концентрации в 50-метровом слое отличается от описанного выше. Средняя концентрация нитратов имеет минимум в начале года, что вызвано зимне-весенним цветением фитопланктона, во время которого потребляется максимальное количество нитратов. В целом за рассматриваемый 15-летний период средняя концентрация нитратов имеет отрицательный линейный тренд и, как следствие, уменьшаются биомассы фитопланктона и зоопланктона (рис. 1).

Внутригодовую изменчивость параметров экосистемы в верхнем слое Черного моря можно рассмотреть на основе климатических данных, полученных путем осреднения по 15 годам рассчитанных полей. Более подробно остановимся на распределении нитратов и суммарного фитопланктона. Эти два параметра экосистемы Черного моря выбраны из-за того, что первый – основное питательное вещество в системе, а фитопланктон – продуцент. На рис. 3 представлен годовой ход средней по площади бассейна концентрации нитратов и суммарного фитопланктона в верхнем слое Черного моря, полученный по результатам моделирования.

 

Рис. 3. Сезонный ход концентрации нитратов (а) и фитопланктона в верхнем слое моря (б), полученный по результатам моделирования. I–XII – месяцы года.

 

Максимум концентрации нитратов наблюдается на глубине ~80 м. Выше и ниже расположены слои верхнего и нижнего нитроклина. В зимнее время благодаря интенсивному перемешиванию максимум концентрации заметно падает за счет переноса большого количества нитратов к поверхности моря; соответственно, растет поверхностная концентрация. В остальные сезоны поверхностная концентрация нитратов в глубоководной части моря мала. В это время значения средней концентрации нитратов на поверхности (см. рис. 3) заметные за счет влияния северо-западного шельфа моря, где количество нитратов велико в течение всего года из-за их поступлений с реками. На рис. 4 представлены средние профили нитратов для всех сезонов (черная сплошная линия соответствует зиме, черная пунктирная – весне, серая сплошная – лету, серая пунктирная – осени) для глубоководной части Черного моря и СЗШ. В глубоководной части моря максимум сохраняется примерно на одной и той же глубине для всех сезонов, а концентрация нитратов на поверхности близка к нулю, за исключением зимы. На шельфе поверхностная концентрация сохраняет большие значения в течение всего года, причем в зимнее время профиль распределения нитратов имеет максимум на поверхности.

 

Рис. 4. Климатические профили распределения средних по площади концентраций нитратов (а, б) и фитопланктона (в, г). Профили (а) и (в) соответствуют глубоководной части Черного моря, (б) и (г) – северо-западному шельфу.

 

Из-за высокой концентрации питательных веществ, в зимнее время происходит интенсивное цветение фитопланктона на поверхности моря, которое в весеннее время с увеличением интенсивности солнечного света опускается в более глубокие слои. Затем интенсивность поступления нитратов в поверхностный слой моря из нижележащих слоев ослабевает и максимум концентрации фитопланктона опускается на глубину ~30 м (летний подповерхностный максимум). В конце осени – начале зимы максимальные концентрации фитопланктона вновь наблюдаются на поверхности моря. Причем, согласно рис. 4, если в глубоководной части профили распределения концентраций фитопланк-тона соответствуют описанной выше схеме, то на шельфе максимум концентрации наблюдается в приповерхностном слое для всех сезонов. Это объясняется достаточно высоким содержанием нитратов на шельфе и более высокой освещенностью поверхностного слоя по сравнению с нижележащими слоями.

Распределение нитратов и фитопланктона на поверхности Черного моря для четырех сезонов, полученное путем осреднения за 15 лет, показано на рис. 5 и 6.

 

Рис. 5. Климатические распределения поверхностной концентрации нитратов для четырех сезонов.

 

Рис. 6. Климатические распределения поверхностной концентрации фитопланктона для четырех сезонов.

 

Максимальные концентрации нитратов наблюдаются в зимний сезон, когда они значительны и в глубоководной части бассейна. В летнее время большие концентрации заметны только на СЗШ возле устьев рек, впадающих в Черное море. В глубоководной части моря поверхностная концентрация нитратов близка к нулю. В осенний сезон повышенное содержание нитратов наблюдается вдоль всего западного берега. Это происходит за счет переноса богатых биогенами вод северо-западного шельфа Основным черноморским течением.

Зимне-весеннее цветение фитопланктона определяет его максимальную поверхностную концентрацию в зимний период. В этот сезон высокие значения отслеживаются по всей акватории Черного моря. Поверхностная концентрация фитопланктона в глубоководной части бассейна ненамного ниже, чем на шельфе. Похожая картина наблюдается и во время осеннего цветения. Весной и особенно летом поверхностная концентрация фитопланктона в глубоководной части Черного моря резко отличается от его концентрации на СЗШ, где она существенно больше. Повышенные концентрации в эти сезоны прослеживаются также вдоль западного берега, куда биопродукция выносится циклоническим течением с северо-западного шельфа.

На рис. 7 и 8 представлены зональные распределения концентраций нитратов и фитопланктона на основе полученных по результатам моделирования климатических полей.

 

Рис. 7. Зональные распределения концентрации нитратов для четырех сезонов.

 

Рис. 8. Зональные распределения концентрации фитопланктона для четырех сезонов.

 

Основной элемент вертикальной структуры распределения нитратов (рис. 7) – максимум концентрации на глубине ~80 м, который прослеживается во все сезоны и вдоль всего разреза. Летом этот максимум достигает наибольшей величины. В зимний сезон из-за интенсивного перемешивания содержание нитратов в приповерхностном слое выравнивается с содержанием нитратов в более глубоких слоях. Следует отметить также повышение концентрации в верхнем слое моря возле западного берега, вызванное стоком рек на СЗШ, главным образом Дуная. Другая особенность распределения нитратов – опускание изолиний нижнего нитроклина возле берегов, что вызвано поднятием вод в центре бассейна, что в свою очередь вызвано крупномасштабной циклонической циркуляцией в Черном море.

Распределение фитопланктона на разрезе в зимний и осенний сезоны имеет монотонный характер: концентрация повышается к поверхности моря до максимальных значений. Весной и летом максимум концентрации фитопланктона опускается на глубину ~30 м, что связано с увеличением световой радиации на глубине и уменьшением концентрации нитратов в поверхностном слое глубоководной части моря. В эти сезоны также наблюдается увеличение концентрации фитопланктона в поверхностном слое моря возле западного берега (рис. 8).

СРАВНЕНИЕ РЕЗУЛЬТАТОВ МОДЕЛИРОВАНИЯ С ДАННЫМИ НАБЛЮДЕНИЙ

Для оценки качества полученных по результатам моделирования биогеохимических полей необходимо сравнить их с соответствующими данными измерений. Биологических данных in situ в Черном море за рассматриваемый период было приведено немного. Авторы статьи использовали для этих целей междисциплинарную океанографическую базу данных для Черного моря, составленную в рамках проекта NATO “Science for Stability TU-Black Sea project”. Для сравнения с результатами моделирования использовались наиболее представленные данные in situ измерений кислорода, нитратов и хлорофилла-а. На рис. 9 приведена карта с обозначенными положениями станций, на которых проводились измерения этих параметров и которые далее были использованы для сравнения (часть данных отбракована). Результаты моделирования интерполировались по пространству и по времени в те точки пространства и времени, где имелись данные измерений. Больше всего имеется данных измерений по кислороду, меньше всего – по хлорофиллу-а.

 

Рис. 9. Места выполнения измерений концентрации кислорода (серые кружки), нитратов (треугольники) и хлорофилла-а (квадратики). Линия внутри бассейна Черного моря – изобата 200 м.

 

На рис. 10 представлены профили концентрации кислорода, полученные по данным измерений, и соответствующие им профили, полученные по результатам моделирования для глубоководной части моря и СЗШ. Также представлены профили среднего и среднеквадратичного отклонения модельных значений от измеренных.

 

Рис. 10. Профили концентрации кислорода, полученные по данным измерений (сплошная линия), и соответствующие им профили, полученные по результатам моделирования (пунктир), для глубоководной части моря (а) и СЗШ (в); б и г – соответствующие профили среднего отклонения (сплошная линия) и среднеквадратичного отклонения (пунктир).

 

Для глубоководной части моря профили концентрации кислорода, полученные по результатам измерений и моделирования, достаточно близки между собой до глубины ~70 м, затем среднее отклонение модельного профиля от измеренного резко возрастает. Это происходит из-за того, что оксиклин, по результатам моделирования, шире в этом месте, чем по данным измерений. На шельфе средние профили практически совпадают до глубины ~30 м, затем расходятся.

На рис. 11 представлены профили распределения нитратов и отклонений данных измерений от результатов моделирования для глубоководной и шельфовой частей Черного моря. Для глубоководной части средний профиль концентрации нитратов, полученный по результатам моделирования, достаточно близок к профилю, полученному по данным измерений. Основное отличие в том, что верхний нитроклин по модельным данным расположен несколько выше, чем по данным измерений. Это приводит к систематической ошибке, как видно на рис. 11б. Среднеквадратичное отклонение имеет два максимума на горизонтах 40 и 110 м. Эти горизонты соответствуют средним положениям верхнего и нижнего нитроклинов. Хотя в средних профилях они достаточно близки, а нижние нитроклины по результатам моделирования и по данным измерений практически совпадают, большие градиенты концентрации нитратов даже при незначительных ошибках в положении нитроклина приводят к значительным величинам среднеквадратичного отклонения. На СЗШ концентрации нитратов систематически превышают величины, полученные по данным контактных измерений. Это, по-видимому, объясняется неточным заданием количества поступающих со стоком рек нитратов, которое известно только оценочно.

 

Рис. 11. Профили концентрации нитратов, полученные по данным измерений (сплошная линия), и соответствующие им профили, полученные по результатам моделирования (пунктир), для глубоководной части моря (а) и СЗШ (в); б и г – соответствующие профили среднего отклонения (сплошная линия) и среднеквадратичного отклонения (пунктир).

 

Следующий параметр, по которому проводилось сравнение результатов моделирования с данными измерений, – хлорофилл-а. Модельные концентрации хлорофилла пересчитывались из суммарной концентрации фитопланктона. На рис. 12 приведены средние профили и отклонения модельных значений хлорофилла-а от данных измерений для глубоководной части Черного моря. Данные для шельфа практически отсутствуют, поэтому они не анализируются. Профиль концентрации хлорофилла-а, полученный по результатам моделирования, довольно близок к профилю, полученному по данным измерений. Наибольшие среднеквадратичные отклонения наблюдаются у поверхности и на горизонте 20 м. На поверхности эта ошибка вызвана, по-видимому, тем, что ассимилируемые спутниковые данные – это композиция за две недели, поэтому они сглажены. Данные контактных измерений, с которыми сравниваются результаты моделирования, могут содержать элементы мелкомасштабной изменчивости, которые пропадают в используемых спутниковых данных. Максимум на глубине 20 м связан с тем, что в весенний и летний сезоны (когда было выполнено основное количество измерений) на этой глубине наблюдаются максимальные градиенты концентрации фитопланктона (рис. 8).

 

Рис. 12. Профили концентрации хлорофилла, мг/м3, полученные по данным измерений (сплошная линия), и соответствующие им профили, полученные по результатам моделирования (пунктир), для глубоководной части моря (а); соответствующие профили среднего отклонения (сплошная линия) и среднеквадратичного отклонения (пунктир) (б).

 

ЗАКЛЮЧЕНИЕ

Представлены результаты изучения динамики экосистемы Черного моря за 15-летний период — с 1998 по 2012 г. Исследуемые биогеохимические поля на регулярной сетке получены с использованием трехмерной модели нижнего пищевого уровня морской экосистемы с ассимиляцией данных спутниковых цветовых сканеров.

Результаты анализа выявили межгодовые изменения параметров экосистемы в верхнем слое Черного моря. Биомассы фитопланктона и зоопланктона имеют отрицательный линейный тренд, который сопровождается уменьшением содержания нитратов в поверхностном слое. Анализ сезонной изменчивости параметров экосистемы проведен на основе климатических полей концентраций нитратов и фитопланктона, полученных осреднением по всему периоду моделирования. Он показал, что результаты численного моделирования качественно правильно описывают сезонный ход и пространственное распределение основных параметров экосистемы Черного моря.

Проведено сравнение полученных полей с доступными данными измерений концентраций кислорода, нитратов и хлорофилла-а. Это сравнение показало, что полученные в результате моделирования параметры экосистемы Черного моря не только качественно, но и количественно согласуются с имеющимися данными измерений. Наибольшее несоответствие между результатами моделирования и данными измерений наблюдается на шельфе. Это связано в первую очередь с недостаточной точностью величин задаваемых потоков нитратов в местах впадения рек в Черное море. Другой причиной может быть неточность воспроизведения гидродинамических полей на шельфе моделью циркуляции. В глубоководной части моря наибольшие ошибки возникают для горизонтов, на которых наблюдаются максимумы градиентов соответствующих параметров.

About the authors

V. L. Dorofeev

Marine Hydrophysical Institute, Russian Academy of Sciences

Author for correspondence.
Email: l.sukhikh@gmail.com
Russian Federation, 2A, Kapitanskaya street, Sevastopol, 299011

L. I. Sukhikh

Marine Hydrophysical Institute, Russian Academy of Sciences

Email: l.sukhikh@gmail.com
Russian Federation, 2A, Kapitanskaya street, Sevastopol, 299011

References

Supplementary files