CC BY
10
УДК [574.5:551.463.5] (262.5.04): [520.626.05+681.785]-021.272
DOI: 10.22449/0233-7584-2021-2-233-246
Биооптические характеристики прибрежных вод Черного моря вблизи Севастополя: оценка точности спутниковых продуктов, восстановленных по данным MODIS и VIIRS
Е. Ю. Скороход 1 Т. Я. Чурилова 1, Т. В. Ефимова 1, Н. А. Моисеева 1, В. В. Суслин 2
1 Институт биологии южных морей имени А. О. Ковалевского РАН, Севастополь, Россия 2 Морской гидрофизический институт РАН, Севастополь, Россия н elenaskorokhod@ibss-ras
Поступила в редакцию 16.04.2020 г., после доработки - 26.10.2020 г.
Цель. Цель работы - оценить точность спутниковых продуктов для прибрежных вод Севастополя, восстановленных стандартными алгоритмами по данным спектрорадиометров MODIS, установленных на искусственных спутниках Земли Aqua и Terra, и VIIRS на спутнике Suomi NPP. Методы и результаты. Отбор проб in situ проводился на станции, расположенной на расстоянии двух миль от бухты Севастопольской в точке с координатами 44° 37' 26" с. ш. и 33° 26' 05" в. д. Для определения концентрации хлорофилла а использовали спектрофотометрический метод. Спектральные показатели поглощения света оптически активными компонентами измеряли в соответствии с современным протоколом NASA. Использовались данные MODIS и VIIRS 2-го уровня с пространственным разрешением 1 км вокруг станции. Обработку спутниковых данных проводили при помощи программного обеспечения SeaDAS 7.5.3, разработанного NASA. Исследование показало, что содержание оптически активных компонентов в прибрежных водах Севастополя по данным MODIS и VIIRS при использовании стандартных алгоритмов определяется некорректно: в сравнении с данными измерений in situ значения концентрации хлорофилла а в среднем весной меньше в 1,6 раза, а летом - больше в 1,4 раза; вклад пигментов фитопланктона в общее поглощение света на длине волны 443 нм в среднем меньше в 8,7 раза; поглощение света окрашенным растворенным органическим веществом в сумме с неживой взвесью в среднем больше в 2,2 раза.
Выводы. Стандартные алгоритмы NASA не применимы для расчета биооптических показателей вод (концентрация хлорофилла a, показатель поглощения света пигментами фитопланктона и показатель поглощения света окрашенным растворенным органическим веществом в сумме с неживой взвесью) в прибрежном районе Черного моря вблизи Севастополя. Чтобы использовать спутниковые данные для экологического мониторинга, необходимо развивать региональный алгоритм, учитывающий оптические особенности вод в этом районе, в частности в сложной прибрежной зоне.
Ключевые слова: хлорофилл a, фитопланктон, окрашенное растворенное органическое вещество, неживое взвешенное вещество, дистанционное зондирование, MODIS, VIIRS, Черное море
Благодарности: работа выполнена по теме ГЗ № АААА-А19-119061190081-9 и при финансовой поддержке проекта РФФИ № 18-45-920070. Авторы выражают глубокую признательность м. н. с. Е. Г. Сахоню и вед. инженеру Е. А. Землянской (ФИЦ ИнБЮМ, г. Севастополь) за помощь в отборе проб при регулярном биооптическом мониторинге в бухте Севастопольской. Также авторы благодарят Goddard Space Flight Center и Ocean Biology Processing Group (NASA) за предоставленные данные и их предварительную обработку.
Для цитирования: Биооптические характеристики прибрежных вод Черного моря вблизи Севастополя: оценка точности спутниковых продуктов, восстановленных по данным MODIS и VIIRS / Е. Ю. Скороход [и др.] // Морской гидрофизический журнал. 2021. Т. 37, № 2. С. 233246. doi:10.22449/0233-7584-2021-2-233-246
© Скороход Е. Ю., Чурилова Т. Я., Ефимова Т. В., Моисеева Н. А., Суслин В. В., 2021
Bio-Optical Characteristics of the Black Sea Coastal Waters near Sevastopol: Assessment of the MODIS and VIIRS Products Accuracy
E. Yu. Skorokhod 1 T. Ya. Churilova \ T. V. Efimova \ N. A. Moiseeva 1, V. V. Suslin 2
1 A. O. Kovalevsky Institute of Biology of the Southern Seas, Russian Academy of Science,
Sevastopol, Russia
2 Marine Hydrophysical Institute, Russian Academy of Science, Sevastopol, Russia H elenaskorokhod@ibss-ras
Purpose. The purpose of the work is to evaluate accuracy of the satellite products for the coastal waters near Sevastopol, reconstructed by the standard algorithms based on the MODIS and VIIRS (installed at the artificial Earth satellites Aqua and Terra, and at Suomi NPP, respectively) data. Methods and Results. In situ sampling was carried out at the station (44°37'26" N and 33°26'05" E) located at a distance of two miles from the Sevastopol Bay. The chlorophyll a concentration was measured by the spectrophotometric method. The spectral light absorption coefficients by optically active components were measured in accordance with the current NASA protocol. The spectroradiometers MODIS and VIIRS Level 2 data with spatial resolution 1 km in nadir around the in situ station (44°37'26"±0°00'32" N and 33°26'05"±0°00'54" E) were used. The satellite products were processed by the SeaDAS 7.5.3 software developed in NASA. The research showed that the standard NASA algorithms being applied to the MODIS and VIIRS data, yielded incorrect values of the optically active components' content in the Black Sea coastal waters near Sevastopol as compared to the data of in situ measurements in the same region: the satellite-derived "chlorophyll a concentration" was on average 1.6 times lower in spring, and 1.4 times higher in summer; the contribution of phytoplankton pigments to total light absorption at 443 nm was underestimated in 8.7 times; the light absorption by colored detrital organic matter was overestimated in 2.2 times.
Conclusions. The NASA standard algorithms are inapplicable to calculating bio-optical indices in the coastal waters of the Black Sea near Sevastopol since they provide incorrect values of the satellite products (Ca-s, aph-s(443) and acDM-s(443)). Operative ecological monitoring based on satellite data requires development of a regional algorithm taking into account the seawater optical features in the region and in the coastal zone, in particular.
Keywords: chlorophyll a, phytoplankton, colored dissolved organic matter, non-algal particulate, remote sensing, MODIS, VIIRS, Black Sea
Acknowledgments: The investigation was carried out within the framework of the state task on theme No. АААА-А19-119061190081-9 and at the RFBR financial support of project No. 18-45-920070. The authors are very thankful to the junior researcher Sakhon E. G. and the leading engineer Zemlianskaia E. A. (IBSS, Sevastopol) for their assistance in sampling during regular bio-optical monitoring in the Sevastopol Bay. The authors are grateful to Goddard Space Flight Center and Ocean Biology Processing Group (NASA) for the data provided and their preliminary processing.
For citation: Skorokhod, E.Yu., Churilova, T.Ya., Efimova, T.V., Moiseeva, N.A. and Suslin, V.V., 2021. Bio-Optical Characteristics of the Black Sea Coastal Waters near Sevastopol: Assessment of MODIS and VIIRS Products Accuracy. Physical Oceanography, [e-journal] 28(2), pp. 215-227. doi:10.22449/1573-160X-2021-2-215-227
Введение
Дистанционное зондирование активно используется для решения различных научных задач, в том числе оперативного мониторинга состояния водных экосистем. Большими преимуществами дистанционного зондирования являются возможность проведения долговременных серий наблюдений с большим пространственным покрытием.
В водных экосистемах отдельного внимания заслуживают прибрежные воды, так как они находятся под сильным антропогенным воздействием. С речным/береговым стоком в прибрежные воды поступает взвешенное и окрашенное растворенное органическое вещество вместе с биогенными элементами, что приводит к увеличению биомассы фитопланктона в прибрежных водах [1]. Фитопланктон, неживое взвешенное вещество и окрашенное растворенное органическое вещество являются оптически активными компонентами среды (ОАК). Содержание ОАК в воде и их спектральные оптические свойства (поглощение и рассеяние света) [2-4] обусловливают восходящее из водной толщи излучение, которое регистрируется спутниковыми инструментами [5]. Вариабельность вклада ОАК в общее поглощение света является одной из причин неточности моделирования спутниковых продуктов при использовании стандартных алгоритмов.
Концентрация хлорофилла а (Ca-с) является одним из стандартных спутниковых продуктов, характеризующих продуктивность вод. Для использования спутниковых продуктов необходимо убедиться, что они корректно характеризуют состояние исследуемой акватории. В настоящее время не известны работы, в которых были бы представлены результаты сравнения спутниковых продуктов MODIS и VIIRS с близкими в пространстве и во времени данными измерений in situ в прибрежных водах Черного моря в разные сезоны года.
Многолетний массив данных регулярного биооптического мониторинга, осуществляемого в прибрежных водах Черного моря в районе бухты Севастопольской начиная с 2009 г., представляет собой уникальную базу для валидации спутниковых данных.
Цель работы состоит в оценке точности восстановленных спутниковых продуктов стандартными алгоритмами по результатам спектральных измерений спектрорадиометрами MODIS, установленными на искусственных спутниках Земли Aqua (MA) 1 и Terra (MT) 2, и VIIRS 3 на Suomi NPP (V). С этой целью выполнено сравнение концентрации хлорофилла a и показателей поглощения света ОАК, рассчитанных по спутниковым и по натурным данным в прибрежных водах Севастополя в отдельные сезоны.
Предварительные результаты были представлены на всероссийской научной конференции [6].
Материал и методы
Данные in situ. В рамках регулярного биооптического мониторинга отобраны пробы в поверхностном слое Черного моря на станции, расположенной напротив бухты Севастопольской на расстоянии двух миль от берега в точке с координатами 44° 37' 26" с. ш. и 33° 26' 05" в. д. Глубина на станции составляет ~ 62 м.
Для определения концентрации хлорофилла a (Ca-i) использовали спектрофо-тометрический метод [7]. Показатели поглощения света взвешенным веществом
1 URL: https://oceancolor.gsfc.nasa.gov/data/aqua/ (дата обращения: 21.10.2020).
2 URL: https://oceancolor.gsfc.nasa.gov/data/terra/ (дата обращения: 21.10.2020).
3 URL: https://oceancolor.gsfc.nasa.gov/data/viirs-snpp/ (дата обращения: 21.10.2020).
(ap-i(X)) и окрашенным растворенным органическим веществом (acDOM-i(^)) измеряли в соответствии с протоколом NASA [8, 9]. Пробы морской воды осаждали на стекловолоконных фильтрах типа GF/F (Whatman) диаметром 25 мм с порами 0,7 мкм при слабом вакууме (не более 0,2 атм). Для получения пробы окрашенного растворенного органического вещества морскую воду пропускали через мембранные фильтры диаметром 47 мм с порами 0,2 мкм.
Спектрофотометрические измерения проводились сразу после отбора проб на двулучевых спектрофотометрах: SPECORD M40 (Carl Zeiss Jena) в 20092014 гг. и LAMBDA 35 (PerkinElmer) с интегрирующей сферой в 2015-2019 гг.
Поглощение света взвешенным веществом (ap_i(X)) разделяли на поглощение света пигментами фитопланктона (aph-i(X)) и поглощение света неживой взвесью (aNAP-i(^)) путем обесцвечивания пигментов взвеси на фильтре в соответствии с методами [10, 11].
Спутниковые данные. Использовались данные 2-го уровня 4 с пространственным разрешением 1 км.
При анализе были использованы следующие спутниковые продукты:
- концентрация хлорофилла a (Ca-c), мг-м-3, рассчитанная на основе двух стандартных алгоритмов NASA: Ocean chlorophyll algorithm (OCx) [12] и Hu color index (CI) [13];
- показатель поглощения света пигментами фитопланктона на длине волны 443 нм (aph-o(443)), м-1, рассчитанный по модели NASA Generalized Inherent Optical Property (GIOP) [14, 15];
- показатель поглощения света acDOM-o(^) в сумме с aNAP-o(^) на длине волны 443 нм, (acDM-o(443)), м-1, рассчитанный по модели NASA GIOP.
Обработка спутниковых данных проводилась при помощи программного обеспечения SeaDAS 7.5.3, разработанного NASA 5. Спутниковые данные выбирали из области вокруг станции. Область имеет форму квадрата с центром в точке с координатами (44° 37' 26" ± 0° 00' 32") с. ш. и (33° 26' 05" ± 0° 00' 54") в. д.
При выборе спутниковых данных предпочтение отдавали более географически близким к точке натурных измерений. При отсутствии спутниковых данных в выбранной области стандартные продукты получали путем аппроксимации значений соседних пикселей (минимум два прилегающих пикселя с противоположных сторон).
Результаты и обсуждение
Сравнение спутниковых продуктов с данными измерений in situ проводили в течение всего годового цикла развития фитопланктона. С этой целью был создан массив данных, содержащий результаты натурных и спутниковых измерений.
Концентрация хлорофилла а (Ca). Ошибка расчета Ca-c в зимний период в среднем составила ±38 %. Отклонения от результатов натурных измерений наблюдались как в сторону более высоких, так и в сторону более низких значений (рис. 1, a).
4 URL: https://oceancolor.gsfc.nasa.gov/docs/format/l2nc/ (дата обращения: 21.10.2020).
5 URL: https://seadas.gsfc.nasa.gov/ (дата обращения: 21.10.2020).
Р и с. 1. Сравнение спутниковых продуктов Ca-c (a) и aph-с (443) (b) с результатами натурных биооптических измерений (Ca-i и aph-i(443)). Цветные линии - линии линейных регрессий; красная линия - линейная регрессия для всех спектрорадиометров
F i g. 1. Comparison of the satellite-derived products Ca-s (a) and aph-s(443) (b) with the results of field bio-optical measurements (Ca-i and aph-i(443)). Color lines are the linear regression lines; red line is the linear regression for all the spectroradiometers
Существенное отклонение значений Ca-c от результатов натурных измерений Ca-i отмечено в весенний период. Явно выражен характер ошибки: при высоких значениях Ca-i (Ca-i > 1,1 мг-м-3) Ca-c была ниже Ca-i в 1,0-5,4 раза. Погрешность определения Ca-c весной в среднем составила ±35 %.
Летом характер ошибки менялся: преимущественно наблюдалось превышение значений Ca-c по сравнению с Ca-i в 0,36-2,7 раза, в среднем в 1,4 раза. Значения Ca небольшие, но ошибка расчета Ca-c была весьма существенной -в среднем составила ±55 %.
Осенью, так же как и весной, Ca-с < Ca-i (при Ca-i > 1,1 мг-м-3), преимущественно по данным MT. Ошибка определения Ca-c по данным всех спектрора-диометров в тот сезон в среднем составила ±36 %.
То есть на протяжении всего года спутниковый продукт Ca-c определяется с большой относительной ошибкой (в среднем ±41 %) по данным всех спектрора-диометров, однако все же можно выделить V (ошибка в среднем ±43 % при выборке в 29 значений) как самый неточный из рассматриваемых спектрорадиомет-ров иMA как наилучший (ошибка в среднем ±38 % при выборке в 38 значений).
Спутниковый продукт Ca-c генерируется двумя алгоритмами: CI и OCx. При Ca-с ниже 0,15 мг-м-3 используется алгоритм CI, при Ca-с выше 0,20 мг-м-3 -алгоритм OCx. При значениях Ca-с в диапазоне от 0,15 до 0,20 мг-м-3 алгоритмы CI и OCx сочетаются при помощи взвешенного подхода. Входными данными при генерации Ca-с алгоритмами CI и OCx являются значения спектрального коэффициента яркости Rrs(^). Кроме Rrs(^), алгоритм OCx включает в себя постоянные табличные коэффициенты a,.
Использование одних и тех же коэффициентов без учета типа вод может являться одной из причин возникновения ошибки в определении Ca-с на региональном уровне.
Показатель поглощения света пигментами фитопланктона на длине волны 443 нм (flph(443)). Значения этого параметра, полученные по данным спектрорадиометров (oph-c(443)), на протяжении года почти всегда ниже результатов наблюдений in situ (aph-i(443)). Диапазоны изменения afph-c(443) в среднем в два раза уже диапазона изменения aph-i(443) (рис. 1, b). Таким образом, поглощение света пигментами фитопланктона в прибрежных водах Черного моря, сгенерированное стандартными алгоритмами по спутниковым данным, существенно недооценивается.
Связь между Ca и aph(443). Ранее для большого массива океанических данных была установлена зависимость между Ca и Oph(^), которая описывается степенной функцией [16, 17]. Коэффициенты этой степенной зависимости различны как для разных акваторий [18, 19], так и для одной акватории в разные сезоны [20, 21].
Связь между Ca и aph(X) позволяет достаточно точно восстанавливать Ca по показателю поглощения света с учетом сезонных особенностей влияния условий среды обитания фитопланктона на содержание пигментов и светопогло-щающую способность клеток. Этот подход, особенность которого состоит в оценке Ca на основе Oph(^), используется в активно развивающемся в настоящее время трехканальном алгоритме [22].
Если Ca и aph(^) определены верно, то должна прослеживаться их взаимосвязь. На рис. 2 видно, что степенная зависимость между данными натурных 238 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 37 № 2 2021
наблюдений (aph-i(443) и Ca_i) отмечается во все сезоны, тогда как связь между спутниковыми продуктами (aph-o(443) и Ca-с) почти отсутствует, что свидетельствует о некорректности оценок этих параметров по стандартным алгоритмам NASA.
1
0.1
§ 0,01
з, ?
0.001 0.0001
0,1 1 10
С,, mi м !
С. мгм ;
1 0.1
£
5 0,01 ■j
£
0.001
Q 0001 I_1........*_I........I
0.1 1 10 С,, мг '.г -
Зима
In situ У= 0,04? R- 0,77, и 21 ДЗ: У 0,023 R- = 0,34, п = 6
Весна
insifw: У= 0,058Л2 = 0,51, и = 32 ДЗ: 7= 0,015 А1",«2 = 0,048, л = 15
Лето
Insitu У= 0,048-J?-51, R2 — 0,32, и = 1S ДЗ: 0 OlO A"-45, R-= 0.25, it = 15
Осень
Insitu-. Y= О ■)'■(> Х- R- = 0.55, п = 2(1 ДЗ: Y= R- = 0.26, я = 5
• In sitll
• Дистанционное зондирование (ДЗ)
Р и с. 2. Связь между aph(443) и Ca. Цветные линии - линии тренда F i g. 2. Relation between aph(443) and Ca. Color lines denote the regression lines
Показатель поглощения света окрашенным растворенным органическим веществом в сумме с показателем поглощения света неживой взвесью на длине волны 443 нм (acdm(443)). Диапазон изменчивости aCDM-o(443) (от 0,010 до 1,02 м-1) на протяжении года в пять раз превышает диапазон изменчивости aCDM-i(443) (от 0,060 до 0,23 м-1) (рис. 3, a). Погрешность определения aCDM-o(443) в среднем за год составила ±160 %. При этом наименьшие ошибки наблюдались весной: средняя относительная погрешность определения aCDM-o(443) составила ±82 %.
Спутниковые продукты, восстановленные по данным квазисинхронных измерений разными спектрорадиометрами, различались, но в целом точность определения спутниковых продуктов на протяжении года была сопоставимой: относительная ошибка 8(Ca.c) составила ±38, ±41, ±43 % для MA, MT, V соответственно, а 5(aph-o(443)) ±64, ±54, ±75 %, за исключением 5(aCDM-o(443)), равной ±158, ±207, ±117 %.
Показатель поглощения света всеми ОАК на длине волны 443 нм (fltot(443)). Ошибка определения Ca.c, aph-o(443) и aCDM-o(443), возможно, является следствием некорректного разделения общего поглощения света взвешенным и окрашенным растворенным органическим веществом на aph-c(443) и aCDM-o(443).
Для проверки данной гипотезы проведено сравнение значений atot(443) по данным спектрорадиометров и измерений in situ. Поглощение света чистой водой на длине волны 443 нм (aw(443)) не учитывали, так как оно постоянно.
На протяжении года разница между максимальным и минимальным значениями по спутниковым данным (0,014-0,21 м-1) и по данным измерений in situ (0,051-0,25 м-1) была одинаковой и составила 0,20 м-1, при этом диапазон спутниковых значений был несколько смещен (на 0,037 м-1) относительно данных in situ в сторону меньших величин (рис. 3, b).
Р и с. 3. Сравнение значений acDM-c(443) (a), atot-c(443) (b) с результатами натурных измерений (acDM-i(443) и atot-i(443)). Обозначения см. на рис. 1
F i g. 3. Comparison of the acDM-s (443) (a) and atot-s (443) (b) values with the results of field measurements (acDM-i(443) and atot-i(443)). Designations as in Fig. 1
Однако следует отметить, что сравнение данных в парах квазисинхронных измерений показало, что atot-c(443) ниже atot-i(443) в 0,92-4,1 раза. В одном случае было установлено существенное отличие спутниковых (V) данных (ол-с(443) = 0,014 м-1) от результатов наблюдений in situ (atot-i(443) = 0,21 м-1). Сравнение всех пар
данных выявило, что значения atot-c(443) и atot-i(443) не совпадают. При этом значения atot-c(443) по данным разных спектрорадиометров, полученным в течение одних суток, существенно различаются (в 0,57-7,2 раза) между собой. Из этого следует, что проблема заключается не в точности разделения atot-o(^) на показатели поглощения света ОАК (пигментами фитопланктона и окрашенным растворенным органическим веществом в сумме с неживой взвесью), а в его изначально неточном определении. Сравнение показало, что оценка atot-o(443) ниже atot-i(443) в большинстве случаев (в 0,92-15 раз, в среднемв 2,9 раза).
Вклад ОАК в общее поглощение света по данным измерений in situ (atot-i(^)). Относительный вклад отдельных ОАК в общее поглощение света меняется в зависимости от длины волны, времени года и характеризует оптический тип вод [23-26]. Вклад ОАК в общее поглощение света (atot-i(^)) на длинах волн 440 и 488 нм в отдельные сезоны в поверхностном слое моря представлен на рис. 4.
Р и с. 4. Вклад пигментов фитопланктона (aph-i(X)/atot-i(X)), окрашенного растворенного органического вещества (acDOM-i(X)/atot-i(X)) и неживого взвешенного вещества (aNAP-i (X)/atot-i(X)) в общее поглощение света (atot-i(X)) на длинах волн 440 и 488 нм
F i g. 4.Contribution of the phytoplankton pigments (aph-i(X)/atot-i(X)), the colored dissolved organic matter (acDOM-i(X)/atot-i(X)), and the non-algal particulate (aNAP-i (X)/atot-i(X)) to the light absorption budget (atot-i(X)) at 440 nm and 488 nm
Для исследованных длин волн получено, что поглощение света окрашенным растворенным органическим веществом превалирует над поглощением света пигментами фитопланктона и поглощением света взвешенным веществом во все сезоны. На протяжении года наибольший вклад окрашенного растворенного органического вещества отмечен в осенний период. В среднем его значение составляло 59 и 47 % на длинах волн 440 и 488 нм соответственно. На длине волны 488 нм вклад пигментов фитопланктона в общее поглощение возрастает в сравнении с данными, полученными на длине волны 440 нм. Анализ сезонной динамики показал, что наибольший вклад пигментов фитопланктона в общее поглощение приходится на весну и в среднем составляет 35 и 44 % на длинах волн 440 и 488 нм соответственно. Отмеченная сезонная динамика значений показателей поглощения света отдельных ОАК среды, а так-
же их относительного вклада в общее поглощение света на длинах волн в диапазоне излучения, в котором яркость моря наиболее чувствительна к оптическим свойствам среды, вероятно, является причиной грубых ошибок при оценке ОАК с помощью стандартных алгоритмов.
Исследование показало, что рассмотренные спутниковые продукты некорректно соотносятся с биооптическими показателями вод исследуемой акватории, а также с их пропорциями (таблица). Сезонный ход концентрации хлорофилла а сглаживается: значения Са-с изменяются от 0,18-1,8 мг-м-3 летом до 0,39-1,4 мг-м-3 весной в отличие от отмеченного в природе увеличения Са-i от минимальных значений (0,22-1,5 мг-м-3) летом до максимальных (0,38-2,3 мг-м-3) весной. Ошибка Са-с в зависимости от сезона варьировалась от ±34 до ±55 %, при этом она не являлась однотипной: в весенний период значения Са-с в среднем были ниже Са-i в 1,6 раза, а в летний -выше в 1,4 раза.
Отношение поглощения света пигментами фитопланктона к поглощению окрашенным растворенным органическим веществом в сумме с неживой взвесью (ярь(443)/астм(443) Ratio of light absorption by the phytoplankton pigments to absorption by the colored dissolved organic matter in total with non-algal particles (aph(443)/acDM(443))
Значение / Value in situ MA MT V
Зима / Winter
Min 12/88 3,0/97 - 31/69
Max 51/49 80/20 - 84/16
Mean 32/68 38/62 50/50* 58/42
Весна / Spring
Min 20/80 1,0/99 7,0/93 2/98
Max 58/42 15/85 81/19 60/40
Mean 38/62 9,0/91 36/64 21/79
Лето / Summer
Min 6,0/94 1/99 0,40/100 3,0/97
Max 37/63 14/86 16/84 55/45
Mean 26/71 7,0/93 8,0/92 18/82
Осень / Autumn
Min 24/76 5,0/95 30/70 -
Max 39/61 88/12 88/12 -
Mean 33/67 45/55 51/49 34/66*
* Одна пара данных / One data pair.
Неточность определения концентрации хлорофилла a ведет к ошибкам в определении многих других показателей, расчеты которых базируются на этом параметре. В частности, это приводит к некорректной оценке биомассы фитопланктона, первичной продукции и скорости роста планктонных микроводорослей [27, 28]. Значения acDм-с(443) выше значений acDм-l(443) (в среднем в 2,2 раза), а aph-с(443) ниже aph-l(443) (в среднем в 9 раз) на протяжении всего года, что ведет к некорректному результату расчета относительного вклада фитопланктона и неживых компонентов оптически активных веществ в среде 242 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 37 № 2 2021
(CDOM и NAP), а следовательно, и к ошибке определения оптических свойств вод, в частности прозрачности, спектральных свойств квантовой облученности, проникающей в водную толщу [4, 24, 29] и фотосинтетического потенциала микроводорослей [30-32].
Результаты исследования свидетельствует о том, что применяемые стандартные алгоритмы, основанные на биооптических показателях вод 1 -го типа [33], не подходят для расчета спутниковых продуктов для прибрежных вод Черного моря (в частности, в районе бухты Севастопольской).
Чтобы использовать данные дистанционного зондирования, требуется развивать региональные алгоритмы, учитывающие сезонные особенности спектральных биооптических показателей вод и их соотношение с концентрацией основного фотосинтетически активного пигмента - хлорофилла а в исследуемом регионе. В частности, можно адаптировать трехканальный алгоритм [22], который разделяет общее поглощение света на поглощение растворенным окрашенным органическим веществом в сумме с неживой взвесью и пигментами фитопланктона. В расчете по этому алгоритму используются коэффициенты яркости моря в диапазоне длин волн 480-560 нм. По этим коэффициентам восстанавливается концентрация хлорофилла а с учетом связи между показателями поглощения света пигментами и концентрации хлорофилла а [20, 34].
Заключение
На основе сформированного многолетнего массива натурных биооптических данных проведена оценка точности спутниковых продуктов, рассчитанных по данным спектрорадиометров MA, MT и V. Получены следующие результаты:
1. В среднем относительная ошибка стандартного спутникового продукта «концентрация хлорофилла а» в зависимости от сезона составила от ±34 до ±55 %, при этом ошибка не являлась однотипной: в весенний период значения Са-с в среднем были ниже Ca-i в 1,6 раза, а в летний - выше в 1,4 раза.
2. В процессе исследования выявлено, что при расчете спутниковых продуктов с применением стандартных алгоритмов NASA вклад пигментов фитопланктона в общее поглощение света остается недооцененным: в среднем значения показателя поглощения света пигментами фитопланктона на длине волны 443 нм (aph-o(443)) были в девять раз ниже результатов натурных измерений. В это же время показатель поглощения света растворенным окрашенным органическим веществом в сумме с показателем поглощения света неживой взвесью на длине волны 443 нм (acDM-c(443)) в среднем был переоценен в 2,2 раза.
3. Существенные погрешности в расчете спутниковых продуктов и разница между показателями поглощения света всеми ОАК воды (atot(^)) ставят под сомнение применимость стандартных алгоритмов для исследуемого района моря.
4. Распределение вклада ОАК в общее поглощение света меняется на протяжении года: наибольший вклад CDOM на длинах волн 443 и 488 нм (59 и 47 %) приходится на осень, а наибольший вклад aph-i(X) (35 и 43 %) - на весну.
СПИСОК ЛИТЕРАТУРЫ
1. State of the Environment of the Black Sea (2001 - 2006/7) / T. Oguz (ed.). Istanbul, Turkey, 2008. 448 p. (Publications of the Commission on the Protection of the Black Sea Against Pollution (BSC) ; 2008-3).
2. Surface chlorophyll in the Black Sea over 1978-1986 derived from satellite and in situ data / О. V. Kopelevich [et al.] // Journal of Marine Systems. 2002. Vol. 36, iss. 3-4. P. 145-160. doi:10.1016/S0924-7963(02)00184-7
3. ЕрловН. Г. Оптика моря. Л. : Гидрометиздат, 1980. 248 с.
4. Kirk J. T. O. Light and photosynthesis in aquatic ecosystems. 3d edition. Cambridge : Cambridge University Press, 2011. 662 p.
5. Gordon H. R., Wang M. Retrieval of water-leaving radiance and aerosol optical thickness over the oceans with SeaWiFS: a preliminary algorithm // Applied Optics. 1994. Vol. 33, iss. 3. P. 443-452. doi:10.1364/ao.33.000443
6. Сравнение стандартных продуктов со спектрорадиометров MODIS Aqua/Terra и VIIRS с результатами биооптических измерений в прибрежных водах Севастополя / Е. Ю. Скороход [и др.] // Понт Эвксинский - 2019 : материалы XI Всероссийской научно-практической конференции молодых ученых по проблемам водных экосистем, посвященной памяти д. б. н., проф. С. Б. Гулина, г. Севастополь, 23-27 сентября 2019 г. Севастополь, 2019. С. 134-135.
7. Jeffrey S. W., Humphrey G. F. New spectrophotometry equations for determining chlorophylls a, b, c1 and c2 in higher plants, algae and natural phytoplankton // Biochemie und Physiology der Pflanzen. 1975. Vol. 167, iss. 2. P. 191-194.
8. Roesler C. S. In Situ Bio-Optical Observations on NERACOOS Buoy A01 (2005-2017): multichannel calibrated chlorophyll fluorescence, turbidity, and multispectral incident irradiance and upwelling radiance. Boston : Massachusetts Water Resources Authority, 2018. 19 p.
9. Twardowski M., Rottgers R., Stramski D. The absorption Coefficient, An Overview // IOCCG Ocean Optics and Biogeochemistry Protocols for Satellite Ocean Colour Sensor Validation. Dartmouth, NS, Canada, 2018. Vol. 1.0 : Inherent Optical Property Measurements and Protocols: Absorption Coefficient. Chapter 1. P. 1-17. http://dx.doi.org/10.25607/OBP-119
10. Estimation of the spectral absorption coefficients of phytoplankton in the sea / M. Kishino [et al.] // Bulletin of Marine Science. 1985. Vol. 37, no. 2. P. 634-642.
11. Tassan S., Ferrari G. M. An alternative approach to absorption measurements of aquatic particles retained on filters // Limnology and Oceanography. 1995. Vol. 40, iss. 8. P. 1358-1368. https://doi.org/10.4319/lo.1995.40.8.1358
12. Ocean color chlorophyll algorithms for SeaWiFS / J. E. O'Reilly [et al.] // Journal of Geophysical Research: Oceans. 1998. Vol. 103, iss. C11. P. 24937-24953. https://doi.org/10.1029/98JC02160
13. Hu C., Lee Z., Franz B. Chlorophyll a algorithms for oligotrophic oceans: A novel approach based on the three-band reflectance difference // Journal of Geophysical Research: Oceans. 2012. Vol. 117, iss. C1. C01011. doi:10.1029/2011JC007395
14. Retrieving marine inherent optical properties from satellites using temperature and salinity-dependent backscattering by seawater / P. J. Werdell [et al.] // Optics Express. 2013. Vol. 21, iss. 26. P. 32611-32622. doi:10.1364/OE.21.032611
15. Generalized ocean color inversion model for retrieving marine inherent optical properties / P. J. Werdell [et al.] // Applied Optics. 2013. Vol. 52, iss. 10. P. 2019-2037. http://doi.org/10.1364/ao.52.002019
16. Variability in the chlorophyll-specific absorption coefficients of natural phytoplankton: Analysis and parameterization / A. Bricaud [et al.] // Journal of Geophysical Research: Oceans. 1995. Vol. 100, iss. C7. P. 13321-13332. doi:10.1029/95JC00463
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
Variations of light absorption by suspended particles with chlorophyll a concentration in oceanic (case 1) waters: Analysis and implications for bio-optical models / A. Bricaud [et al.] // Journal of Geophysical Research: Oceans. 1998. Vol. 103, iss. C13. P. 31033-31044.
Cleveland J. S. Regional models for phytoplankton absorption as a function of chlorophyll a concentration // Journal of Geophysical Research: Oceans. 1995. Vol. 100, iss. C7. P. 1333313344. https://doi.org/10.1029/95JC00532
Light absorption properties of southeastern Bering Sea waters: Analysis, parameterization and implications for remote sensing / P. Naik [et al.] // Remote Sensing of Environment. 2013. Vol. 134. P. 120-134. doi:10.1016/j.rse.2013.03.004
Light Absorption by Phytoplankton in the Upper Mixed Layer of the Black Sea: Seasonality and Parametrization / T. Churilova [et al.] // Frontiers in Marine Science. 2017. Vol. 4. 90. doi:10.3389/fmars.2017.00090
Dynamics in pigment concentration and light absorption by phytoplankton, non-algal particles and colored dissolved organic matter in the Black Sea coastal waters (near Sevastopol) / T. Efimova [et al.] // Proceedings of SPIE. 2018. Vol. 10833 : 24th International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics. 108336C. doi: 10.1117/12.2504657
Suslin V., Churilova T. A regional algorithm for separating light absorption by chlorophyll-a and coloured detrital matter in the Black Sea, using 480-560 nm bands from ocean colour scanners // International Journal of Remote Sensing. 2016. Vol. 37, iss. 18. P. 4380-4400. doi: 10.1080/01431161.2016.1211350
Phytoplankton absorption predicts patterns in primary productivity in Australian coastal shelf waters / C. M. Robinson [et al.] // Estuarine, Coastal and Shelf Science. 2017. Vol. 192. P. 1-16. doi:10.1016/j.ecss.2017.04.12
Annual variability in light absorption by particles and colored dissolved organic matter in the Crimean coastal waters (the Black Sea) / T. Churilova [et al.] // Proceedings of SPIE. 2017. Vol. 10466 : 23rd International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics. 104664B. doi:10.1117/12.2288339
Optical characteristics of two contrasting Case 2 waters and their influence on remote sensing algorithms / M. Darecki [et al.] // Continental Shelf Research. 2003. Vol. 23, iss. 3-4. P. 237250. doi:10.1016/s0278-4343(02)00222-4
Light absorption properties and absorption budget of Southeast Pacific waters / A. Bricaud [et al.] // Journal of Geophysical Research: Oceans. 2010. Vol. 115, iss. C8. C08009. doi: 10.1029/2009JC005517
Comparison of algorithms for estimating ocean primary production from surface chlorophyll, temperature, and irradiance / J. Campbell [et al.] // Global Biogeochemical Cycles. 2002. Vol. 16, iss. 3. P. 9-1-9-15. doi:10.1029/2001GB001444
Реанализ долговременных рядов изменения биомассы фитопланктона в открытой части Черного моря по результатам натурных и спутниковых наблюдений / О. В. Кривенко [и др.] // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. Севастополь : МГИ, 2012. Вып. 26, т. 2. С. 185-194.
The Effect of Optical Properties on Secchi Depth and Implications for Eutrophication Management / E. T. Harvey [et al.] // Frontiers in Marine Science. 2019. Vol. 5. 496. doi:10.3389/fmars.2018.00496
Markager S., Vincent W. F. Light absorption by phytoplankton: development of a matching parameter for algal photosynthesis under different spectral regimes // Journal of Plankton Research. 2001. Vol. 23, iss. 12. P. 1373-1384. https://doi.org/10.1093/plankt/23.12.1373
Bracher A. U., Tilzer M. M. Underwater light field and phytoplankton absorbance in different surface water masses of the Atlantic sector of the Southern Ocean // Polar Biology. 2001. Vol. 24, iss. 9. P. 687-696. https://doi.org/10.1007/s003000100269
32. «Цветение» фитопланктона и фотосинтетически активная радиация в прибрежных водах / Т. Я. Чурилова [и др.] // Журнал прикладной спектроскопии. 2019. Т. 86, № 6. С. 976-985.
33. Morel A., Prieur L. Analysis of variations in ocean color // Limnology and Oceanography. 1977. Vol. 22, iss. 4. P. 709-722. https://doi.org/10.4319/lo.1977.22.4.0709
34. Phytoplankton light absorption in the deep chlorophyll maximum layer of the Black Sea / T. Churilova [et al.] // European Journal of Remote Sensing. 2019. Vol. 52, iss. sup. 1. P. 123-136. doi:10.1080/22797254.2018.1533389
Об авторах:
Скороход Елена Юрьевна, младший научный сотрудник, Научно-исследовательский центр геоматики, ФИЦ ИнБЮМ (Россия, 299011, г. Севастополь, пр. Нахимова, д. 2), SPIN-код: 3314-5775, ORCID ID: 0000-0002-3057-3964, elenaskorokhod@ibss-ras
Чурилова Татьяна Яковлевна, кандидат биологических наук, руководитель Научно-исследовательского центра геоматики, ФИЦ ИнБЮМ (Россия, 299011, г. Севастополь, пр. Нахимова, д. 2), SPIN-код: 2238-9533, ORCID ID: 0000-0002-0045-7284, tanya.churilova@ibss-ras.ru
Ефимова Татьяна Владимировна, младший научный сотрудник, Научно-исследовательский центр геоматики, ФИЦ ИнБЮМ (Россия, 299011, г. Севастополь, пр. Нахимова, д. 2), SPIN-код: 1668-0742, ORCID ID: 0000-0003-3908-4160, tefimova@ibss-ras.ru
Моисеева Наталия Александровна, младший научный сотрудник, Научно-исследовательский центр геоматики, ФИЦ ИнБЮМ (Россия, 299011, г. Севастополь, пр. Нахимова, д. 2), SPIN-код: 8946-3315, ORCID ID: 0000-0003-1356-7981, nataliya-moiseeva@yandex.ru
Суслин Вячеслав Владимирович, кандидат физико-математических наук, ведущий научный сотрудник, руководитель отдела динамики океанических процессов, ФГБУН ФИЦ МГИ (Россия, 299011, г. Севастополь, ул. Капитанская, д. 2), SPIN-код: 1681-7926, ORCID ID: 0000-0002-8627-7603, slava.suslin@mhi-ras.ru
