**НЕПРЕРЫВНОЕ ИЗМЕРЕНИЕ**

**ТЕМПЕРАТУРЫ, СОЛЕНОСТИ И КОНЦЕНТРАЦИИ ХЛОРОФИЛЛА**

**В ПОВЕРХНОСТНОМ СЛОЕ БАРЕНЦЕВА МОРЯ (ИЮНЬ, 1993)**

**ШАВЫКИН А.А.**

**Мурманский морской биологический институт КНЦ РАН,**

**ул. Владимирская 17, Мурманск, 183010, Россия.**

1.  ВВЕДЕНИЕ

(Особенности непрерывного измерения концентрации хлорофилла

фитопланктона в воде)

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

Метод измерения кон­центрации хлорофилла фитопланктона по его in vivo
флуоресценции существует с середины шестидесятых годов (Lorenzen, 1966);
с его помощью выявлены важ­ные закономерности пространственного
распределения хлорофилла фитопланк­тона (Карабашев, 1987; и там же см.
ссылки на литературу). Однако в силу ряда методических трудностей
(влияния внешнего освещения, видового состава фито­планктона,
растворенного органического вещества на измерение отношения
флуоресценция/хлорофилл ) этот метод использовался и используется
большей частью как качественный или полуколичественный метод.
Использование in vivo интенсивности флуоресценции хлорофилла
фитопланктона (ИФХ) для измерения концентрации хлорофилла фитопланктона
(КХ) требует решения ряда методических вопросов, что связано с
указанными особенностями измеряемого объекта.

> На основании накопленного опыта и опубликованных данных можно
> утверждать, что существует несколько факторов, влияющих на связь между
> ИФХ и КХ, которые приводят к тому, что погрешность определения КХ по
> ИФХ может быть порядка нескольких сот и даже тысяч процентов, если эти
> факторы при измерениях ИФХ не учитываются. К факторам, учет которых
> обязателен и возможен уже в настоящее время, относятся (Шавыкин,
> 1997): а) фон водной системы - интенсивность флуоресценции
> растворенного органического вещества (ИФР); б) внешняя освещенность
> флуориметрируемого объема (клеток фитопланктона) в момент измерения;
> в) дисперс­ность флуоресцирующих частиц; г) видовой состав
> фитопланктона.
>
> Мы использовали один из вариантов реализации флуориметрического метода
> измерения КХ, когда учитывается и устраняется влияние ИФР и внешней
> освещенности, а относительно частый отбор проб дает возможность снять
> влияние видового состава, путем использования различных градуировочных
> уравнений (Шавыкин, 1995; Шавыкин, 1997).

2.  АППАРАТУРА

(Краткое описание флуориметрического комплекса, использовавшегося для
измерения концентрации хлорофилла, температуры и солености)

Измерения проводились на флуориметрическом комплексе (рис. 1), в состав
которого входят:

-   двухканальный флуориметр \"Квант-7\" с проточными кюветами;

-   безнасосная проточная система подачи воды к флуориметру;

-   измеритель облученности поверхности моря HIQ-2;

-   измеритель температуры и электропроводности воды ЭПТ-65;

-   спутниковая навигационная система GPS RAYSTAR-900;

-   персональная ЭВМ IBM PC 386.

    Флуориметр \"Квант-7\" - серийный флуориметр, разработанный
    Барнаульским ОКБА НПО \"Химавтоматика\". Проточные кюветы
    флуориметра имеют внутреннее поперечное сечение 10х10 мм. Источник
    света - лампа накаливания, частота модуляции светового сигнала - 5
    Гц. Процессор флуориметра производит усреднение измеряемых значений
    за 3 с. Первый канал флуориметра был настроен на измерение
    интенсивности in vivo флуоресценции хлорофилла фитопланктона.
    Возбуждение флуоресценции осуществлялось через стеклянный
    светофильтр СЗС-22 в полосе 400-500 нм (Карабашев, 1987).
    Интенсивность флуоресценции регистрировалась ФЭУ-100 через вторичный
    светофильтр КС-17, пропускающий излучение с длинами волн большими
    660 нм. По второму каналу флуориметра проводилось измерение
    интенсивности флуоресценции РОВ; первичный светофильтр: УФС-6
    (полоса пропускания: 320-380 нм), вторичный светофильтр состоит из
    двух светофильтров: СЗС-22 и ЖС-17, полоса пропускания: 480-500 нм
    (Карабашев, 1987).

    Безнасосная система подачи анализируемой воды к флуориметру.
    Ис­поль­зовалась навесная гидродинамическая (безнасосная) система
    подачи воды к флуориметру. Анализируемая вода отбиралась с горизонта
    2.5 м с помощью водозаборника, разработанного в Атлантическом
    отделении Института океанологии АН СССР (Лила, Ханаев, 1983).
    Водозаборник укреплялся на штанге (рис.1), установленной на
    расстоянии примерно равном 0.3 длины судна от его кормы. Расход воды
    в такой системе обеспечивается за счет динамического подпора на
    входном конце шланга и меняется с изменением скорости движения
    судна. Использовался прозрачный силиконовый шланг (длина 35 м при
    длине от входного отверстия насадки до штуцера проточной кюветы
    флуориметра 15.5 м, внутренний диаметр шланга 8 мм).

    Пример записи аппаратной функции измерительной системы при разных
    скоростях представлен на рис. 2. График зависимости задержки tзад
    гидродинами­ческой системы прокачки от скорости судна V~с~ ,
    определенной по спутниковой навигационной системе GPS, приведённой
    на рис. 3. Подробное описание процедуры определения разрешающей
    способности системы и полученные при этом выводы приведены в работе
    (Шавыкин, 1995).

    Измеритель освещенности поверхности моря. Для учета влияния внешней
    облученности на отношение F~хл~ / C~хл~ и расчета градуировочных
    уравнений C~хл~ = f(F~хл~), как уже отмечалось, необходимо прово­дить
    соответствующие измерения освещенности поверхности моря солнечным
    светом. Мы использовали для этих целей гидроквантометр HIQ-2,
    измеряющий фотосинтетически активную квантовую радиацию. Датчик
    устанавливался на верхней палубе судна. Спектральный диапазон
    измерений ра­вен 360 - 710 нм. Диапазон измерений облученности:
    1-3000 мкЭйнштейн⋅м^-2^⋅с^-1^ с тремя линейными диапазонами 1-30,
    1-300, 1-3000. Погрешность прибора (accuracy) составляет ±3 % от
    полной шкалы. Гидроквантометр имеет аналоговый выход. Измеряемые
    значения облученности вводились в ЭВМ через плату ана­лого-цифрового
    преобразователя (АЦП) DT-2811.

    Измеритель температуры и электропроводности воды. Для измерения
    гидрологических характеристик использовался датчик температуры и
    солености ЭПТ-65, установленный на глубине 2.5 м в месте от­бора воды
    в систему прокачки (рис.1). Аналоговый сигнал датчика вводился в ЭВМ
    через плату АЦП DT-2811.

    *Канал измерения электропроводности:*

-   диапазон преобразования равен 0.5 - 5 См·м^-1^ ;

-   по паспортным данным коэффициент преобразования составляет 1
    В/См·м^-1^ ; погрешность преобразования равна ± 1 %;

-   коэффициент преобразования, определенный по сопоставлению с дан­ными,
    полученными на солемере ГМ-65 (семь совместных измерений датчиком in
    situ и на пробах, отобранных батометром), равен 0.991 В/См·м^-1^ и
    находится в пределах погрешности;

-   погрешность измерения солености, определенная по сопоставлению с
    данными, полученными на солемере - не превышает 0.05 ^0^/~00~ ;

    Для пересчета электропроводности в единицы солености использовалась
    Практическая шкала солености 1978 г. (О введении\..., 1982).

    *Канал измерения температуры:*

-   диапазон преобразования равен 0-17 °С;

-   по паспортным данным коэффициент преобразования равен 0.3 В/ °С;
    по­греш­ность преобразования составляла ±0.1 °С;

-   коэффициент преобразования, определенный по сопоставлению с дан­ными,
    полученными на погружных опрокидывающихся термометрах (десять
    со­вместных измерений in situ), и после введения постоянной поправки
    0.130 °С, ра­вен 0.30 В/ °С;

-   погрешность измерения температуры, определенная по сопоставлению с
    дан­ными, полученными на батометрической серии, не превышает 0.05 °С.

Для проверки градуировки датчика в ходе всего рейса, примерно 1 раз в
су­тки, выполнялись совместные измерения температуры и электропроводности
из­мерителем ЭПТ-65 и стандартным батометрическим методом.

Измерение концентрации хлорофиллов экстрактным спектрофо­тометрическим
методом. Чтобы рассчитать градуировочное уравнение C~хл~ = f(Fхл), для
ряда проб измеряли F~хл~, при этом анализируемая флуориметром вода
отбиралась в стеклян­ную емкость, и фитопланктон из отобранной воды
осаждался на мембранный фильтр (ВЛАДИПОР № 5, диаметр - 60 мм, размер
пор - 0.5 мкм). Концентрация хлорофиллов определялась по методике
описанной в ГОСТе 17.1.04.02-90, BaCO~3~ не применялся. Фильтры с
осадком хранились в морозильнике и обрабатывались на берегу. Объем воды,
из которой отфильтровывались клетки фитопланктона, - 1.5 - 4 литра.
Объем экстракта - 5-10 мл. Измерения оптической плотности выпол­нялись на
спектрофотометре СФ-46. Градуировочные уравнения C~хл~ = f(F~хл~)
рас­считывали, используя метод наименьших квадратов. Под C~хл~ далее
подразумева­ется концентрация хлорофилла [а]{.underline}.

3\. МЕТОДИКА НЕПРЕРЫВНОГО ИЗМЕРЕНИЯ ХЛОРОФИЛЛА

Результаты представленные здесь получены на материале собранном в июне
1993г. в 72-м рейсе НИС «Дальние Зеленцы» в Баренцевом море.

Модель измерения. В общий сигнал, измеряемый на флуориметре по каналу
измерения хлорофилла (F~пр~, вклад вносит и рассеяние света -
молекулярное и комбинационное. С учетом этого, можно исходить из
следующей модели измерений (Шавыкин, 1995, Шавыкин, 1997):

F~пр~ = F~хл~ + S~дв~ + F~ров~ , (\*)

где:

  -------- --------------------------------------------------------------
  F~пр~    \- интенсивность флуоресценции пробы, измеренная по каналу,
           настроенному на измерение хлорофилла;

  F~хл~    \- интенсивность in vivo флуоресценции хлорофилла
           фитопланктона, иначе в тексте мы обозначаем ее как ИФХ;

  S~дв~    \- интенсивность света, рассеянного дистиллированной водой по
           каналу измерения хлорофилла;

  F~ров~   \- интенсивность флуоресценции растворенного органического
           вещества соответственно по каналу измерения хлорофилла.
  -------- --------------------------------------------------------------

При работе в Баренцевом море в июне 1993 года в ходе изме­рений с
фильтратом и дистиллированной водой было установлено, что в данном
случае влияние РОВ крайне мало (в пределах погрешности измерения
значения интен­сивности флуоресценции фильтрата морской воды и
дистиллированной воды в районе работ совпадали между собой), таким
образом, членом F~ров~ в нашей модели измерения можно пренебречь.
Соответственно, результаты изме­рений РОВ по второму каналу не
учитываются (предполагается, что составляющую F~ров~ в уравнении (\*)
можно учитывать, используя измерения РОВ по второму каналу - подробно
см. (Шавыкин, 1995, Шавыкин, 1997)). Было определено, что S~дв~ = 0.350.
Далее мы принимаем разность F~пр~ - S~дв~ за интенсивность флуоресценции
хлорофилла F~хл~ .

Оценка вида градуировочного уравнения. Из совместных измерений F~хл~ и
C~хл~ , выполненных в Баренцевом море, видно, что почти для одного и
того же значения F~хл~ значения C~хл~ могут различаться в несколько раз
и даже почти на по­рядок (рис. 4). Поэтому вряд ли целесообразно
использовать градуировочное уравнение, рассчитанное непосредственно
только по измерениям C~хл~ и F~хл~. Вместе с тем, точки в координатах
{E; C~хл~ / F~хл~ } можно разделить на три группы, для которых
характерна линейная зависимость отноше­ния C~хл~ / F~хл~ от E: группы А,
В и С (различные обозначения на рис. 5).

Указанное разделение на три района независимо обосновывается
соотно­шениями между пигментами фитопланктона: 1) доли хл [а]{.underline}
в сумме молекул хл [а]{.underline}, [в]{.underline}, [с]{.underline} и
феофитина [а]{.underline} и 2) соотношением между хл [а]{.underline} и
хл [в]{.underline} (рис. 6). Видно, что все точки (пункты отбора проб)
можно отнести к трем различным районам (А, В, С), характеризующимися
своими соотношениями между пигментами. При этом учи­тывалось и
пространственное расположение (\"связность\") точек на маршруте (рис.
7). Районы считаются различными, если при выбранном уровне значимо­сти
(0.10) средние значения отношений пигментов для этих районов значимо
различаются между собой, что имело место и в нашем случае (см.
доверительные интервалы на рис. 6). Подобное разделение на три района
примерно отражает и деление на рай­оны с различными видовыми комплексами
фитопланктона, полученное гидро­биологами при анализе количественного и
качественного состава фитопланктона - рис. 7 (Шавыкин, 1997). Всего нами
была отобрана 21 проба для определения концентрации пигментов
фитопланктона.

Градуировочные уравнения для отдельных районов (j = A, B, C)
представлены в табл. 1.

Таблица 1

Градуировочные уравнения для пересчета F~хл~ в C~хл~ для различных

районов (C~хл~ - нмоль/л; F~хл~ - отн.ед.; E - мкЭ·м^-2^·с^-1^ )

+------------+---------------------------+------+--------+----+------+
| Район,     | Модель                    | П    |        |    |      |
|            |                           | арам |        |    |      |
| о          |                           | етры |        |    |      |
| бозначение |                           | мо   |        |    |      |
|            |                           | дели |        |    |      |
+------------+---------------------------+------+--------+----+------+
| и число    | Y = a + b · x             | S~a~ | S~b~   | S~ | R2   |
|            |                           |      |        | y. |      |
| точек      |                           |      |        | x~ |      |
+------------+---------------------------+------+--------+----+------+
|            |                           |      |        |    |      |
+------------+---------------------------+------+--------+----+------+
|            | Уравнение C~хл~/ F~хл~ =  |      |        |    |      |
|            | f (E):                    |      |        |    |      |
+------------+---------------------------+------+--------+----+------+
| A ()       | C~хл~ / F~хл~ =0.40 +     | 0.14 | 0      | 0. | 0.57 |
|            | 0.00081 ⋅ E               |      | .00035 | 21 |      |
| N = 6      |                           |      |        |    |      |
+------------+---------------------------+------+--------+----+------+
| B ()       | C~хл~ / F~хл~ =0.76 +     | 0.15 | 0      | 0. | 0.88 |
|            | 0.0033 ⋅ E                |      | .00047 | 24 |      |
| N = 9      |                           |      |        |    |      |
+------------+---------------------------+------+--------+----+------+
| C ()       | C~хл~ / F~хл~ =0.98 +     | 0.17 | 0      | 0. | 0.75 |
|            | 0.00094 ⋅ E               |      | .00031 | 18 |      |
| N = 5      |                           |      |        |    |      |
+------------+---------------------------+------+--------+----+------+

Примечание:

  ----------------- --------------------------------------------------------------
  Для района B (при 
  N=9) отсутствует  
  точка,            
  обозначенная "х"  
  (рис. 4-7) -      
  удалена на основе 
  статистического   
  критерия.         

  S~a~ , S~b~       \- стандартные отклонения для коэффициентов a и b,

  S~y.x~            \- стандартное отклонение значений Y при заданных значениях X,

  R^2^              \- квадрат коэффициента корреляции.
  ----------------- --------------------------------------------------------------

Использование освещенности поверхности моря. Далее предполагается, что
значения F~хл~, используемые для вычисления C~хл~ сдвинуты на время
задержки, определяемое скоростью судна. Это одна из первых операций при
обработке данных. Потому далее все измеряемые значения (время,
навигаци­онные параметры, температура, электропроводность, освещенность,
интенсивность флуоресценции по двум каналам) соответствуют одной и той
же анализируемой воде.

При расчетах параметров градуировочных уравнений (табл. 1)
использова­лась средняя освещенность поверхности моря за интервал отбора
проб для спек­трофотометрических измерений C~хл~ (две-три минуты).
Оказалось, что такая мо­дель, где значения всех величин соответствуют
одному конкретному моменту времени не вполне корректна: сильные
флуктуации освещенности E ( E(0) на рис. 8) никак не отражаются на
регистрируемой интенсивности флуоресценции F~хл~. Это несколько
неожиданный результат, если учитывать, во-первых, опублико­ванные данные
(Карабашев, 1987; Abbot et al., 1982; Loftus, Seliger, 1975; Stram­sky,
Dickey, 1992), в соответствии с которыми резкие изменения Е обычно сразу
же приводят к резким изменениям F~хл~, и, во-вторых, выявленной в данной
работе зависимостью F~хл~ / C~хл~ от Е (рис. 5, табл. 1).
Соответственно, использование полученных градуировочных уравнений
приводит к тому, что рассчитанные значения C~хл~ также \"приобретают\"
пики и провалы как и освещенность Е. Это может быть устранено, если
использовать не мгновенное значение Е, а ее среднее значе­ние за время
∆t~уср~ предшествующее моменту измерения F~хл~. В нашем случае время
усреднения ∆t~уср~ должно быть порядка 10-20 минут. Мы брали при
расчетах ∆t~уср~ = 20 мин. Подробное обоснование см. в (Шавыкин, 1995,
Шавыкин, 1997).

4\. РЕЗУЛЬТАТЫ НЕПРЕРЫВНЫХ ИЗМЕРЕНИЙ И ИХ ОБРАБОТКА

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

Горизонтальное распределение температуры, солености, плотности морской
воды и кон­центрации хло­ро­фил­ла см. на рис. 10a-10g. Для обработанных
фрагментов галсов координаты начальных и конечных точек галса, время
выполнения измерений на галсе, средние и некоторые другие значения
величин, которые измерялись на галсе, представлены в таблице 2.

1.  Основные процедуры обработки

Одна из первых процедур обработки массивов данных - визуальный просмотр
графиков получаемых значений изме­ряемых величин и отбраковка резко
выделяющихся значений. Для проведения дальнейших расчетов (вычисления
спектральных и корреляцион­ных характеристик и др.), в которых
предполагается, что значения измеренных величин равномерно
дискретизированы в пространстве, использовалась проце­дура равномерной
дискретизации см. ( Шавыкин, 1997).

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

Устранение задержки, обусловленной движением воды в системе про­качки.
Интенсивность флуоресценции F~пр~ измеряется по первому каналу
флуори­метра после прохождения анализируемой водой подающего шланга.
Соответственно, эти значения сдвинуты во времени по отношению к
значениям величин температуры, электропроводности, освещенности,
навигационным параметрам. Мы всегда использовали режим, при котором
T~зад~ для каждого момента времени определяется по таблице значений
\[V~с~ , T~зад~\] введенной в калибровочный файл. Соответственно
проводился сдвиг значений F~хл~ , чтобы привести все значения к одному
моменту времени. С учетом сдвига значений КХ, ряд значений КХ оказался
на 4-6 значений короче, чем ряд значений T и S.

Усреднение значений освещенности поверхности моря. С учетом методики
выполнения измерений, усреднение значений освещенности поверхности моря
производится по формулам, аналогичным для усреднения скорости судна, но
значения освещенности усредняются за интервал ∆T~E~ только до заданного
момента времени. Длина файла этой процедурой не изменяется. ∆T~E~ = 20
мин.

4.2. Спектральные плотности измеренных величин

Для анализа изменчивости полей T, S, Chl в указанном масштабе (100 м -
100 км) нами использовался математический аппарат цифрового
спект­раль­но­го анализа. Формулы расчета спектров соот­ветст­вуют таковым в
(Дженкинс, Ваттс, 1971, 1972). Наши данные соответствуют значениям
сигнала, отсчитанного через интервал в 40 м.

Графики спектральной плотности температуры и солености воды,
кон­центрации хлорофилла в воде для семи анализируемых галсов
предс­тав­ле­ны на рис. 11 - 13. На этих же рисунках показаны
аппроксимирующие прямые
![](media/image1.emf){width="0.9583333333333334in"
height="0.2361111111111111in"} для различ­ных значений *n*. Введены
следующие обозначения спектров галсов относящихся к различным районам:
для района A - штриховая линия, района B - тонкая сплошная и для района
С - толстая сплошная линия. На приведенных рисунках видно сходство
соответствующих графиков между собой.

Для детального анализа спектральных плотностей температуры, солености и
концентрации хлорофилла на отдельных рисунках были построены
нормированные спектральные плотности этих величин для каждого из галсов
в отдельности (см. примеры спектров для трех галсов на рис. 14, 15, 16).
Под рисунками приведены уравнения аппроксимирующих прямых вида
![](media/image2.emf){width="0.7222222222222222in"
height="0.19444444444444445in"}, полученные методом наименьших квадратов
(МНК) по выделенному интервалу волновых чисел.

Подробный анализ полученных спектров представлен в работе (Шавыкин,
Бойцов, 1997). Основные выводы этого анализа сводятся к следующему.

Спектры измеренных параметров (температуры, солености воды, концентрации
хлорофилла в воде) в основном соответствуют известным теоретическим
моделям. На спектрах всех величин выделяются участки "закона пяти
третей" и кубической зависимости. Для спектров солености вполне
отчетливо можно выделить горизонтальный участок (n≈0), который
располагается между участками с n=3 и n=5/3. Значение волнового числа k
для точки перегиба (перехода зависимости с n=3 в зависимость с n=5/3 )
для температуры примерно постоянно для всех галсов и составляет 2-3
км^-1^ , для хлорофилла она меняется в диапазоне 0.3 - 2.0 км^-1^ в
зависимости от скорости поверхностного течения. Все спектры в
определенной степени согласуются друг с другом: 1) имеют общий участок с
наклоном n≈3 и n≈5/3; 2) значения волнового числа k для точек перегиба
спектров концентрации хлорофилла и температуры близки к друг другу и их
несовпадение определяется рядом специфических для данного
океанологического параметра особенностей; 3) кроме того, для одного и
того же галса точки перегиба кривых спектров температуры и, в
определенной степени и концентрации хлорофилла, отвечают значению
волнового числа k, лежащему примерно в середине диапазона,
соответствующего горизонтальному участку спектра солености.

ЛИТЕРАТУРА

**ГОСТ 17.1.04.02.-90.** Вода. Методика спектрофотометрического
определения хлорофилла [а]{.underline}. М.: Изд-во стандартов, 1990. -
16 с.

**Дженкинс Г., Ваттс Д***.* Спектральный анализ и его приложения.
Вып. 1. - М.: Мир, 1971. - 320 с.

**Дженкинс Г., Ваттс Д.** Спектральный анализ и его приложения. Вып.
2. - М.: Мир, 1972. - 288 с.

**Карабашев Г.С.** Флюоресценция в океане - Л.: Гидрометеоиздат, 1987. -
200 с.

**Лила А.Н., Ханаев С.А.** Двухканальный проточный флуориметр ДПФ-81 для
морских исследований // Океанология. - 1983. - Т.23, вып. 5. - С.
897-899.

**О введении шкалы практической солености,** 1978 и нового
международного уравнения состояния морской воды, 1980 // Океанология. -
1982. - Т. 22, № 2. - С. 337-343.

**Шавыкин А.А.** Флуориметрический метод непрерывного измерения
концентрации хлорофилла фитопланктона в воде в проточном режиме.:
Припер. - Апатиты, 1995. 54 с.

**Шавыкин А.А.** Контактный флуориметрический метод непрерывного
измерения концентрации хлорофилла фитопланктона в воде // Планктон морей
Западной Арктики. - Апатиты, 1997. - С. 185- 223.

**Шавыкин А.А., Бойцов В.Д.** Характеристики пространственной
изменчивости хлорофилла фитопланктона, температуры, солености
поверхностного слоя южной части Баренцева моря // Планктон морей
Западной Арктики. - Апатиты, 1997. - С. 223 - 265.

**Abbott M.R., Richerson P.J., Powell T.M.** In situ response of
hytoplankton fluorescence to rapid variations in light // Limnol.
Oceanogr. - 1982. - V. 27. N 2. - P. 218-225.

**Loftus M.E., Seliger H.H.** Some limitation of the in vivo
fluorescence technique // Chesapeake Sci. - 1975. - V. 16, N 2. - P.
79-92.

**Lorenzen C.J.** A method for the continuous measurement of in vivo
chlorophyll concentration // Deep - Sea Res. - 1966. - V. 13. - P. 223 -
227.

**Stramsky R.W., Dickey T.D.** Shot-term variations of the biooptical
properties of the ocean in response to cloud-induced irradiance
fluctuations // J. Geophys. Res. - 1992. - V.97, N C4. - P. 5713- 5721.