2008
Известия ТИНРО
Том 153
УДК 628.394(26)
Т.С. Лишавская, А.В. Мощенко
Дальневосточный научно-исследовательский гидрометеорологический институт, Институт биологии моря ДВО РАН, г. Владивосток [email protected], [email protected]
НЕКОТОРЫЕ МЕТАЛЛЫ В ДОННЫХ ОТЛОЖЕНИЯХ ПРИБРЕЖНЫХ АКВАТОРИЙ ОСТРОВА САХАЛИН
В июне—июле 2001 г. с целью определения фонового состояния морской среды на участках будущей установки нефтедобывающих платформ в районах трасс морского трубопровода и других объектов инфраструктуры по фазе 2 развития проекта "Сахалин-2" были выполнены экологические съемки на акваториях ряда районов сахалинского шельфа. В результате определены фоновые концентрации кислоторастворимых форм Al, As, Ba, Cd, Cr, Cu, Fe, Pb, Hg, Zn в донных отложениях северо-восточного шельфа о. Сахалин, портов Кайгон (зал. Набиль), Поронайск (зал. Терпения), Холмск и зал. Анива. Антропогенное обогащение донных отложений металлами в 2001 г. было минимальным на исследованных акваториях северо-восточного шельфа о. Сахалин и максимальным — в порту Холмск; в порту Поронайск степень техногенного воздействия была несколько выше, чем в зал. Анива. При помощи факторного и корреляционного анализа показано, что степень соответствия друг другу пространственных распределений кислоторастворимых форм металлов в донных отложениях возрастает при увеличении уровня загрязнения среды; увеличивается и сила связи их концентраций и некоторых физико-химических характеристик грунтов.
Lishavskaya T.S., Moshchenko A.V. Metals in bottom sediments at the coast of Sakhalin Island // Izv. TINRO. — 2008. — Vol. 153. — P. 295-311.
Background concentrations of acid-soluble forms of Al, As, Ba, Cd, Cr, Cu, Fe, Pb, Hg, and Zn were measured in bottom sediments on the northeastern shelf of Sakhalin Island and in the Nabil Bay (port Kajgon), the Terpeniya (Patience) Bay (port Poronajsk), the Aniva Bay, and in the port Kholmsk in June—July, 2001 to characterize marine environments at the sites of projected construction of oil-producing platforms, along the possible rout of underwater pipeline, and at the sites of other infrastructure objects of the project "Sakhalin-2", phase 2. Anthropogenic enrichment of bottom sediments by metals was minimal at the northeastern shelf and maximal in the port Kholmsk; in the port Poronajsk the degree of technogenic influence was higher than in the Aniva Bay. There is determined by means of factor and correlation analysis that degree of matching the spatial distributions of acid-soluble forms of different metals in bottom sediments increases with pollution; moreover, correlation between their concentrations and various physical and chemical parameters of sediments increases, as well. This relation allows to receive integral parameters for characteristic of anthropogenic influence on marine environments.
Введение
Комплексные экологические исследования северо-восточного шельфа о. Сахалин были развернуты в конце 80-х гг. прошлого века. В 1990-е гг. специалисты Дальневосточного научно-исследовательского гидрометеорологического института (ДВНИГМИ) выполнили несколько экспедиций в районы Пильтун-Астохского и Луньского месторождений нефти. Были изучены фоновые характеристики планктонных и бентосных сообществ, а также содержание и распределение различных классов загрязняющих веществ в морской воде и донных отложениях (Tkalin, 1993; Tkalin, Belan, 1993; Ткалин и др., 1999).
В 1992, 1995 и 1996 гг. компания "Сахалин Энерджи Инвестмент Компании", оператор проекта "Сахалин-2", организовала фоновые экологические исследования на Пильтун-Астохской нефтегазоносной площади, результаты которых изложены в отчетах ДВНИГМИ (CSA, 1996, 1997). Летом 1998 г. на Пиль-тун-Астохской площади была установлена нефтедобывающая платформа "Мо-ликпак" и пробурена разведочная скважина. Перед установкой платформы, в июне 1998 г., была выполнена комплексная съемка по программе экологического мониторинга. В дальнейшем изменения в морской среде отслеживали в сентябре—октябре 1998-2006 гг. (Олейник, 1999; Ткалин и др., 1999; Белан, Олейник, 2000; Лишавская и др., 2005; Мощенко и др., 2005; Lishavskaya et al., 2005).
В водных экосистемах соединения тяжелых металлов входят в число приоритетных загрязняющих веществ и могут оказывать негативное воздействие на гидробионтов. Кроме того, изменение концентрации металлов в различных компонентах экосистем указывает на характер и интенсивность процессов миграции вещества и может быть использовано для оценки общего уровня антропогенной нагрузки (Шулькин, 2004).
Цель данной работы — определение фоновых концентраций и оценка распределения кислоторастворимых форм металлов в грунтах на участках установки нефтедобывающих платформ, в районах инсталляции морского трубопровода и других объектов инфраструктуры проекта "Сахалин-2" (фаза 2), а также выяснение зависимости их концентраций и распределения от основных физико-химических характеристик донных отложений.
Материалы и методы
Работы выполняли в течение рейса НИС "Павел Гордиенко" в июне—июле 2001 г. на двух полигонах в районах Пильтун-Астохского (ПА) и Луньского (ЛУ) месторождений, вдоль меридионального и широтных разрезов, соответствующих трассе трубопровода (полигон МШР), в портах Кайгон (зал. Набиль), Поронайск (зал. Терпения), Холмск и зал. Анива (рис. 1). Донные отложения отбирали дночерпателем Ван Вина с площадью захвата 0,2 м2. Пробы для анализа металлов (Al, As, Ba, Cd, Cr, Cu, Fe, Pb, Hg, Zn), нефтяных углеводородов (НУ), органического углерода (C ) и определения гранулометрического состава брали из ненарушенного поверхностного слоя грунта (2 см) с помощью тефлоновых пробоотборников и замораживали. Для измерения активной реакции среды (рН) и окислительно-восстановительного потенциала (ЕН) из дночерпателя отбирали пробы объемом не менее 1 л, чтобы посредством центрифугирования получить достаточное количество поровых вод (10 мл).
Содержание алюминия, железа и цинка измеряли атомно-абсорбцион-ным методом в пламенном режиме (в дальнейшем — FLAA) на спектрофотометре фирмы "Perkin-Elmer", модель В-3030, с дейтериевой коррекцией фона (UNEP, 1995). Концентрации мышьяка, бария, кадмия, меди, хрома и свинца определяли этим же методом, но в режиме электротермической атомизации (GFAA) на спектрофотометре той же фирмы, модель Z-3030, с графитовой печью HGA-600 и коррекцией фона на основе эффекта Зеемана. Содержание ртути измеряли методом "холодного пара" (CVAA ) анализатором MAS-50.
296
гон Луньского месторождения; 4, 5, 6, 7 — соответственно широтные разрезы 1, 2, 3 и 4; 8, 9 — районы предполагаемой установки буровых платформ; 10 — место расположения платформы "Моликпак", 11 — основной полигон в зал. Анива напротив пос. Пригородное; 12 — участок выемки грунта в бухте Лососей, 13 — участок свалки грунта (5 станций по крестообразной схеме)
Fig. 1. Schematic map of the region studied at the shelf of Sakhalin Island: 1 — main polygon at Piltun-Astokhskoe deposit; 2 — meridian section; 3 — main polygon at Lunskoe deposit; 4, 5, 6, 7 — latitudinal sections 1, 2, 3, and 4, correspondingly; 8, 9 — reputed regions of the installation of oil-extracting platforms; 10 — site of Moliqpak platform location, 11 — main polygon in Aniva Bay opposite Prigorodnoe township; 12 — escavation site in Lososey Bay, 13 — damping site (5 station along cruciform scheme)
Металлы извлекали из осадков по стандартным методикам (UNEP, 1995). При подготовке образца для анализа методами GFAA и CVAA 1-2 г гомогенизированного грунта естественной влажности взвешивали с точностью до 0,01 г, переносили в стеклянный стакан, приливали 10 мл разбавленной (1 : 1) HNO3, перемешивали, закрывали пробу стеклом и нагревали на водяной бане при температуре 95 °С в течение 10-15 мин. Пробу охлаждали, добавляли 5 мл концентрированной HNO3, закрывали стеклом и нагревали 30 мин без кипячения. Вновь
охлаждали, приливали еще 5 мл концентрированной HNO3, упаривали на водяной бане при температуре 95 °С до 5 мл, добавляли 2 мл бидистиллированной воды, 3 мл 30 %-ной H2O2 и нагревали до прекращения вспенивания. Затем прибавляли по 1 мл H2O2 (не более 10 мл) с последующим нагреванием до минимального вспенивания. Образец разбавляли, фильтровали и доводили объем до 50 мл. При подготовке образца для анализа методом FLAA после обработки H2O2 к пробе прибавляли 5 мл концентрированной HCI, 10 мл воды, нагревали 15 мин на водяной бане, фильтровали и также доводили раствор до 50 мл.
Влияние матрицы анализируемых образцов на правильность анализа оценивали по извлечению известного количества внесенных металлов. Контроль качества включал оценку возможного загрязнения проб в ходе подготовки и анализа, оценку правильности проводимых измерений и оценку воспроизводимости результатов. Пробы грунтов анализировали сериями, включавшими по 11 образцов каждая. Для оценки возможного загрязнения проб в ходе подготовки и измерений одновременно с полевыми образцами в каждой серии проводили анализ холостой пробы. Воспроизводимость результатов подтверждали повторным анализом одного полевого образца в каждой серии анализов. Все измерения содержания металлов выполняли в НПО "Тайфун" (г. Обнинск). Методы анализа, пределы обнаружения металлов в донных отложениях и результаты контроля качества суммированы в табл. 1.
Нефтяные углеводороды экстрагировали гексаном, экстракт очищали с помощью колоночной хроматографии на окиси алюминия и газохроматографи-чески определяли суммарную концентрацию НУ пламенно-ионизационным детектором. Интервал измерения — 0,5-500,0 мкг/г массы сухого образца. Для контроля качества использовали суррогатные и внутренние стандарты, вносимые в пробы на разных стадиях анализа (UNEP, 1992; Методика ..., 2000).
Органический углерод определяли методом высокотемпературного каталитического сжигания в токе кислорода на приборе Т0С-5000А фирмы " Shimadzu", снабженном специальным модулем SSM-5000A для сожжения твердых образцов и суспензий (900 оС). В результате сжигания весь Сорг окисляется до С02, концентрация которого фиксировалась инфракрасным газовым анализатором. Пределы измеряемых концентраций С — 0,1-30,0 мг; стандартное отклонение — 1-8 %. °рг
Гранулометрический состав донных отложений анализировали ситовым и ареометрическим методами. В результате анализа определяли процентное содержание фракций > 10,0, 10,0-5,0, 5,0-2,0, 2,0-1,0, 1,0-0,5, 0,50-0,25, 0,250,10, 0,10-0,05, 0,05-0,01, 0,010-0,005 и менее 0,005 мм (UNEP, 1995).
Активную реакцию среды определяли рН-метром МР230 фирмы "METTLER TOLEDO" с комбинированным электродом InLab-413 в режиме автоматической температурной компенсации. Относительная точность измерения ±0,004 ед. рН. Окислительно-восстановительный потенциал измеряли мультипроцессорным рН-метром HI-9321 с комбинированным платиновым игольчатым электродом HI-3230 (электродная система Pt/PtO) фирмы "HANNA instruments" с точностью ±1 мВ. Стабильность работы платинового электрода проверяли по стандартному раствору HI-7020.
При статистической обработке результатов измерений применяли стандартные процедуры и тесты, предлагаемые пакетом прикладных программ STATISTICA 6.0 (Боровиков, Боровиков, 1998). Группирование полигонов по концентрациям металлов выполняли с помощью кластерного анализа. Структуру взаимосвязей металлов исследовали факторным анализом. Факторы выделяли методом главных компонент. Для упрощения их структуры и получения легко интерпретируемых решений использовали различные способы вращения осей — варимакс, квартимакс и т.п. (Афифи, Эйзен, 1982; Ким и др., 1989; Носов, 1990; StatSoft, 2001). Количество факторов определяли с соблюдением принципа факторной
298
экономии на основе критериев Кайзера и "каменной осыпи". В соответствии с принципом экономии из двух конкурирующих решений выбирают модели с минимальным числом общих факторов. По критерию Кайзера принимают факторы с собственными числами > 1, т.е. когда изменения фактора объясняют > 10 % дисперсии измеряемых переменных. Тест "каменной осыпи" предполагает выбор числа факторов, приходящихся на область наибольшего уплощения кривой на графике зависимости величины собственного числа от своего номера. Ортогональность искомых факторов проверяли иерархическим анализом косоугольных факторов. Построение карт производили посредством ПП Surfer 6.0 (Surfer Mapping System, 1993-2000).
Результаты и их обсуждение
Исследованные акватории заметно различаются по диапазонам глубин опробования. Некоторые из них простираются от литорали до средней и даже нижней зон шельфа, до глубины 49,5-75,0 м (полигон МШР и зал. Анива), другие приурочены лишь или к его верхней (порты Кайгон, Поро-найск, Холмск), или к средней и нижней частям (ПА и ЛУ месторождения) (табл. 2). Все рассматриваемые прибрежные участки о. Сахалин характеризуются значительной изменчивостью физико-химических свойств донных отложений. Так, грунты его северо-восточного шельфа варьируют от гравийных до алевропелитовых, от хорошо сортированных до несортированных. По содержанию Сорг осадки разных участков этого района отличаются друг
см
га <u ^
^ (Л
-н
к к ср
<и о
га
CQ
н о
<и у
га *
к ч
о р
н ж
о X
3 н га н
Л
со
<и р
н о о
га р
га Ж о =к
га р
х к К
ж
<и *
о ч н о
X 3 ж ж о к
CQ
о -
га н <и
к
S
ж
<и ^
р
га Ж
о о
3 -
<и К
<и р
с
3 К о н <и
<v
х
*
0
<и у
S н
S -
га
1
<
-а
с
га
-о ш -3 з
-о
<и
о .о
.о га
о <v
си -а
-о
о
га
С <
■cf
СО
см СО
О О
а
<и со
га р
О о
=s 3 ж ж
га
CQ
о р
S
а к
S
н р
<и О
=s
3 не
S о
S т
н ^
О С
^ 2
к S
О Ч
о ж
<и р
<и
о ж га са
о р
S S
S
н р
<и О
к га Ж
ж
<и у
g о ~ о
CQ g
т ^
га *
га о^
<и S
ж
<и
о О 'Р
CQ О
н
га
К °
га & СЬ 4
см
<V ^OI
£
4 о са
г! § „
С о р
с
к к
5 г^ с 2
ж
о о
н ж
<и S <и
(П
CM LO
I I
о от
ОО
СО ^
^ S2 i 1
^ оО
CD
I I I I
CM Q 3 CO CD
. io>
- CD CD —i
I I
-o
CM
LO CO LO
S <=> °q t-<q от ^ oo ю ^
CDCMCOO^OoO~^CO
OnfflTtlMCOONlDOO
^oo'o <D~cO~£ ^^^ ^соотсоотюр,—
ooooooorooo
0CM0CD01--0—'^O О CD~ LO CO CO C^ 1—1 CM t-T ЮП^МОО Ю ЮСО CO '—I '—I ^ —I
ю от
■cf «Э Lf^ «э Lf^ ОО от
ci СЭ ОТ ОТ ОО сг от ОО С5 ОТ от С5 от
Ю5 С5 - СО
С5 ao tC ОТ ОТ ОО Li^ О ооою ОТ ОО ОО
С5 ao С5 1Я О со ci
СО с^ Li^ о с^ ai
о
CM Ю
CO N о о о
LO C5 C5 C5 1 C5 ^^ О vvvvoovvov
опмммооююо
C5 C5 C5 C5 C5 О
VV о О V
V V
CM Ю
ONNOWWOOWO сгсгсгсгсгсгсгсгсгсо lii C5 C5 <D C5 C5 ci C5 dS О
ьоооооьооь
~ to a "a з ^ e
CQ
к о н aj
S <a a
5 ю
g о
s о к
«з 2
p X
о к
0 га \о oa га о
1 ср
5 а
S к
к
н р
<v о
н о
^J га н
% гБ
2 и
S О
ф с
^ m
га Ж га К
I s
§ ®
з га
u \о
cq о
о „ •
=s 3 ж ж о к а
о р
О
О -
« с^ О
S
И s
° I
га 1
13
S: s
^ к
. . а
cq га
о со
со s
к S
S °
га ^ Ж га га
=s
о X о <и У К
К
<и
СЬ р
Н
О р
н X
3 S
a
^ cu
а л
к
*
<и р
Некоторые физико-химические характеристики донных отложений исследованных акваторий о. Сахалин
Table 2
Some physical and chemical characteristics of bottom sediments of the studied water areas of Sakhalin Island
Район т Н г PG PS АР А1 Ре С орг рН ЕН РНС
Пильтун-Астохское месторождение, 0,474 0,449 7,6 89,9 2,4 2,01 0,41 0,53 7,411 191,4 0,56
п = 57, глубина 26,0-55,5 м 0,139 0,223 0,0 3,3 0,1 0,10 0,00 0,06 7,000 61,0 <0,50
4,683 0,808 72,2 84,9 18,4 18,10 4,90 2,80 7,670 258,0 0,94
Полигон МШР, n = 111, 0,321 0,396 3,4 89,9 6,7 5,69 1,02 1,26 7,443 163,1 1,34
глубина 0-49,5 м 0,033 0,122 0,0 0,4 0,1 0,10 0,00 0,04 6,940 29,0 <0,50
1,357 0,851 54,8 91,5 92,9 92,70 14,90 6,96 7,840 233,0 6,96
Луньское месторождение, п = 43, 0,222 0,459 2,5 85,2 12,4 9,26 3,09 2,64 7,292 156,4 0,79
глубина 36,5-61,5 м 0,060 0,234 0,0 3,9 0,3 0,32 0,00 0,32 6,930 93,0 <0,50
1,327 0,816 30,7 81,7 88,7 83,20 8,40 11,00 7,440 205,0 5,96
Порт Кайгон, п = 10, глубина 1,8-9,0 м 0,906 0,525 16,2 83,0 0,7 5,20 1,23 0,29 7,443 139,9 1,01
0,216 0,313 0,0 7,0 0,1 0,10 0,00 0,06 5,880 89,0 <0,50
3,518 0,762 57,3 66,3 1,9 92,70 14,90 1,31 7,760 197,0 2,52
Северо-восточный шельф в целом, 0,342 0,423 4,3 88,9 6,7 0,53 0,20 1,35 7,404 169,4 0,84
п = 221, глубина 0-61,5 м 0,033 0,122 0,0 0,4 0,1 0,10 0,00 0,04 6,930 29,0 <0,50
4,683 0,851 72,2 91,5 92,9 1,40 0,50 11,00 7,840 258,0 6,96
Порт Поронайск, п = 10, 0,130 0,402 0,4 67,4 32,2 28,93 3,27 4,76 6,837 166,3 2,29
глубина 0-3,5 м 0,024 0,091 0,0 13,5 0,4 0,20 0,20 0,54 6,050 66,0 <0,50
0,268 0,635 3,9 96,2 96,5 80,70 15,80 17,44 7,270 224,0 5,18
Зал. Анива, п = 26, глубина 0-75 м 0,699 0,641 18,0 36,8 45,3 34,47 10,84 13,78 7,440 143,1 4,02
0,011 0,277 0,0 3,3 0,7 0,40 0,00 0,68 6,650 38,0 <0,50
6,039 0,934 88,3 82,2 95,6 81,90 42,00 47,65 7,770 182,0 26,40
Порт Холмск, п = 7, глубина 4,0-9,5 м 0,208 0,644 3,8 48,9 47,3 37,61 9,71 40,78 7,381 138,1 247,60
0,033 0,410 0,0 14,2 0,2 0,20 0,00 0,72 6,800 114,0 <0,50
0,903 0,786 17,2 67,3 75,7 64,30 18,70 94,30 7,680 161,0 527,00
Примечание, т — средний размер зерен, мм; Нг — энтропия гранулометрического распределения; Рй, РБ, АР, А1, Ре — доли гальки и гравия, песка, алевропелитов, алевритов и пелитов, %; Сорг — содержание органического углерода, %; ЕН — окислительно-восстановительный потенциал, мВ; РИС — нефтяные углеводороды, мкг/г; первая строка — среднее, вторая — минимальное, третья — максимальное значение.
от друга более чем в 275 раз. Окислительно-восстановительные условия изменяются от промежуточных между окислительными и восстановительными до восстановительных, а реакция среды поровых вод — от слабокислой (почти нейтральной) до слабощелочной. Концентрации НУ, измеренные в донных отложениях северо-восточного шельфа, а также в портах Кайгон и Поронайск, невелики и сравнимы с их содержанием в осадках приполярных морей, но существенно, в десятки раз, возрастают в грунтах акватории порта Холмск.
Фоновые концентрации металлов и экологическое состояние морской среды. По содержанию металлов наиболее похожими и наименее обогащенными оказались осадки северо-восточного шельфа о. Сахалин — ПА и ЛУ месторождений, полигона МШР и порта Кайгон (рис. 2, табл. 3). В этом отношении порт Поронайск более похож на зал. Анива, чем на районы, изученные на северо-востоке, а наиболее оригинальным и максимально загрязненным оказался порт Холмск. Концентрации кислоторастворимых форм металлов, приведенные в табл. 3, следует использовать как фоновые для характеристики и мониторинга среды в процессе выполнения второй фазы работ по проекту Сахалин-2.
пам ' ' ' ' ' ' ' '
МШР
лм ■_
Кайгон
Поронайск
зал. Анива-
Холмск-
0 10 20 30 40 50 60 70 80 90
Квадрат Эвклидова расстояния
Рис. 2. Дендрограмма сходства исследованных акваторий по содержанию кислото-растворимых форм металлов в донных отложениях, полученная методом Уорда: ПАМ — Пильтун-Астохское и ЛМ — Луньское месторождения, МШР — полигон меридионального и широтных разрезов
Fig. 2. Similarity dendrogramm of the water areas studied obtain by using the concentrations of metal acid-soluble forms in bottom sediments: ПАМ — Piltun-Astokhskoe and ЛМ — Lunskoe deposits, МШР — polygon of meridian and latitudinal sections
На основе анализа экспериментальных данных о частоте биологических реакций при действии на организмы загрязненных осадков были определены два уровня концентраций различных поллютантов — ERL (effects range-low) и ERM (effects range-medium) (Long et al., 1995). При концентрации ниже ERL неблагоприятное воздействие наблюдается редко (в среднем от 1,9 до 27,3 % случаев). В диапазоне значений, равных и выше ERL, но ниже ERM, отрицательное воздействие отмечено в 11,1-75,0 % случаев. Концентрации на уровне ERM и выше оказывают вредное влияние в 16,9-100,0 % случаев.
Количество большинства кислоторастворимых форм металлов на исследованных участках северо-восточного шельфа о. Сахалин, а также в порту Поронайск заметно ниже соответствующих ERL-концентраций, что свидетельствует о
Статистические характеристики содержания кислоторастворимых форм металлов в донных отложениях некоторых прибрежных акваторий о. Сахалин, мкг/г сухого веса грунта
Table 3
Statistical characteristics of the contents of acid-soluble forms of metals in bottom sediments of some inshore areas of Sakhalin Island
Район
Пильтун- Полигон Луньское Порт Северо- Порт Зал. Анива, Порт
г-. Астохское МШР, месторож- Кайгон, восточный Поронайск, п=26 Холмск, рпг* рол л*
С'ЛсМсНТ 111 1/"ч 1 1/л т LLKL, LLKIVI
месторож- n= 111 дение, п=10 шельф п=10 п=7
дение, п=57 п=43 в целом,
п=221
Алюминий 845/379 228-2420 1351/1036 169-5690 1925/692 170-3910 822/223 571-1350 1308/899 169-5690 4529/1509 2990-7280 7632/2279 4600-15020 7567/1940 3490-8860 - -
Мышьяк 2.99/1.83 1,25-12,90 2.35/1.04 0,60-5,96 3.39/0.63 2,25-5,49 1.81/0.54 1,29-2,93 2.69/1.30 0,60-12,90 3.53/1.38 1,82-6,70 10.50/5.59 4,65-28,10 7.44/2.35 4,34-11,40 8,2 70,0
Барий 4.86/4.50 0,94-29,60 7.37/24.30 0,97-258,00 11.93/13.10 4,32-91,90 2.76/1.80 0,96-7,11 7.40/18.43 0,94-258,00 83.91 /77.61 25,50-246,00 32.13/31.00 4,78-171,00 146/137 22,50-439,00 - -
Кадмий 0.013/0.007 0,002-0,029 0.017/0.029 0,001-0,288 0.030/0.014 0,012-0,057 0.013/0.007 0,007-0,025 0.018/0.023 0,001-0,288 0.032/0.032 0,006-0,110 0.106/0.083 0,021-0,263 0.374/0.341 0,017-0,930 1,2 9,6
Хром 3.96/3.90 0,15-29,20 5.42/5.76 0,78-40,60 4.70/1.81 1,88-7,93 1.45/0.67 0,39-2,69 4.72/4.69 0,15-40,60 7.16/3.29 3,40-12,60 20.51/8.22 7,72-39,50 22.31/10.20 7,67-40,90 81,0 370,0
Медь 0.57/0.31 0,11-1,53 1.00/0.88 0,08-4,97 1.25/0.38 0,67-2,31 0.64/0.37 0,27-1,61 0.92/0.71 0,08-4,97 7.15/4.51 2,79-13,50 10.17/4.15 4,84-18,00 59.55/68.00 9,12-209,00 34,0 270,0
Железо 2341 /932 824-4160 2447/1107 741-5300 3379/393 2570-4170 1696/517 1100-2540 2567/1029 741-5300 8970/2181 6640-11970 6977/1939 1110-10600 11733/3910 5710-16970 - -
Свинец 1.75/0.60 0,76-3,48 1.90/0.80 0,69-5,22 2.80/0.43 2,04-3,79 1.75/0.35 1,35-2,38 2.03/0.78 0,69-5,22 5.49/2.89 1,73-10,50 7.35/2.48 2,49-12,40 71.88/128.00 7,60-361,00 46,7 218,0
Ртуть 0.005/0.001 0,005-0,009 0.005/0.001 0,005-0,009 0.005/0.001 0,005-0,008 0.005/0.000 0,005-0,006 0.005/0.001 0,005-0,009 0.044/0.022 0,022-0,089 0.029/0.038 0,005-0,124 0.175/0.146 0,039-0,456 0,15 0,71
Цинк 4.94/2.90 0,30-12,90 8.12/6.03 0,30-28,10 11.49/2.77 4,57-18,50 4.98/1.99 2,20-9,30 7.81/5.21 0,30-28,10 35.72/9.65 26,5-51,5 44.38/9.70 29,9-61,7 153.00/195.00 0,30-571,00 150,0 410,0
Примечание. Над чертой — среднее значение/среднеквадратическое отклонение, под чертой — пределы изменения; жирным шрифтом выделены концентрации, превышающие ERM, а жирным курсивом — £7?£-значения, прочерк — оценка отсутствует. *По: Long et al., 1995.
благоприятной экологической ситуации (табл. 3)*. На акватории ПА месторождения и зал. Анива величину ERL превышают лишь максимальные концентрации мышьяка, в порту Холмск — мышьяка и ртути. Однако на акватории порта Холмск наибольшие концентрации меди, свинца и цинка превышают ERM-значе-ния, что свидетельствует о катастрофическом состоянии морской среды, учитывая синергичность действия отдельных поллютантов.
Это соответствует выводам авторов, сделанным на основе анализа содержания органических загрязнителей — НУ, фенолов и детергентов — в грунтах рассматриваемых акваторий (Лишавская и др., 2005). Так, даже максимальные измеренные концентрации НУ на полигонах МШР, ПА и Луньского месторождений, а также в портах Кайгон и Поронайск были существенно ниже уровня ERL. В зал. Анива наибольшие концентрации НУ превышали ERL в 1,5-2,6 раза, а в порту Холмск уровень загрязнения осадков этими соединениями приближался к ERM-значению.
Структура взаимосвязей и распределение кислоторастворимых форм исследованных металлов. Концентрация металлов в грунтах является результатом действия целого ряда природных и антропогенных факторов и зависит от количества того или иного металла, попадающего в среду, расстояния от источника его поступления до точки измерения, системы течений, а также от физико-химических свойств самих металлов и донных отложений. Антропогенная нагрузка выражается в поступлении и накоплении на дне избыточного количества тяжелых металлов и других загрязняющих веществ (Шулькин, 2004). Поэтому химический состав грунтов является интегральным показателем уровня техногенного пресса и отражает наиболее значимые особенности распределения загрязнителей, в нашем случае металлов.
Для правильной оценки антропогенного влияния на концентрацию металлов в осадках необходимо учитывать природные факторы, основным из которых является изменчивость гранулометрического состава и содержания Сорг. Из-за тяготения большинства металлов к илистым фракциям, наблюдается увеличение их концентраций от песков к алевропелитам (Bradford, Luoma, 1980; Windom et al., 1989). Кроме того, содержание многих металлов положительно коррелирует с концентрацией Cope (Trefry, Presley, 1976; Daskalakis, O'Connor, 1995; Summers et al., 1996).
Учесть влияние механической дифференциации грунтов можно разными способами: анализом отдельных гранулометрических фракций, регрессионным анализом с использованием данных по фракционному составу, нормированием по содержанию тонкозернистых фракций и по концентрации элементов-индикаторов, т.е. элементов, распределение которых по размерным классам осадков достаточно известно и мало изменяется под влиянием антропогенных факторов (Шулькин, 2004). К таким геохимическим индикаторам относятся скандий, алюминий, железо.
Вариации (дисперсия) измеряемых переменных (в нашем случае — концентраций металлов в грунтах) могут быть как случайными, так и результатом воздействия неких скрытых (латентных, не измеряемых напрямую) факторов, которое проявляется в коррелированности переменных (Кулаичев, 1996; Боровиков, Боровиков, 1998). В контексте этой работы к таким факторам, скорее всего, следует отнести антропогенный привнос и механическую дифференциацию, а для оценки их влияния — применить процедуру факторного анализа.
Концентрации металлов можно представить в виде двух факторов, изменения которых объясняют в сумме более 85 % дисперсии измеряемых переменных (табл. 4). Справедливость двухфакторного решения подтверждается результатами теста
* Следует учитывать, что в полной мере сопоставление ERL и ERM-оценок, полученных Лонгом с соавторами (Long et al., 1995) методом полного разложения (валовые концентрации), справедливо для меди, кадмия, ртути и цинка, извлечение которых методами, используемыми в нашей работе, превышает 70 %.
"каменной осыпи" и критерием Кайзера (рис. 3). Вместе с тем первичные кластеры переменных заметно связаны друг с другом (коэффициент корреляции г = 0,583), что обусловлено в первую очередь заметными нагрузками первого фактора на концентрацию железа, а второго — на концентрацию кадмия. Первый фактор контролирует 64,990,6 % дисперсии концентраций Ба, Сй, Си, РЬ, Hg и ¿п. Второй — определяет 60,185,2 % изменчивости содержания Л1, As, Сг и Fe.
В составе первого фактора преобладают так называемые "металлы антропогенного происхождения", его вариации в наибольшей степени (> 63 % дисперсии) объясняются изменениями концентрации НУ, а максимальные значения приурочены к наиболее загрязненной акватории — порту Холмск (табл. 5, рис. 4). В пространственном распределении этого фактора на акваториях всех исследованных портов отчетливо прослеживается тенденция к снижению его величины от внутренних частей в сторону открытого моря (рис. 5). Второй фактор контролирует прежде всего концентрации алюминия и железа — элементов, по которым нормируют содержание других металлов при оценке уровня загрязнения, и соответственно его вариации определяются главным образом содержанием тонких фракций. Максимальные величины фактора отмечены в зал. Анива, грунты которого наиболее разнообразны среди всех исследованных акваторий, а диапазон глубин максимален. В пространственном распределении фактора 2 в портах отсутствует тенденция, наблюдаемая у фактора 1, и его значения могут быть повышены как во внутренних, так и во внешних частях гаваней. Таким образом, первый фактор, скорее всего, характеризует процессы антропогенного загрязнения среды, а второй — аккумуляцию металлов в результате действия комплекса естественных биогеохимических факторов седиментогене-за, и прежде всего механической дифференциации донных отложений.
В пользу такой интерпретации факторов свидетельствует положительная и статистически значимая корреляция фактора 2 с глубиной в большинстве исследованных районов — на полигонах ПА и МШР, в зал. Анива и порту Холмск (на других акваториях она статистически незначима) — и отсутствие таковой у фактора 1 (табл. 5). С глубиной гидродинамическая активность затухает, что способствует отложению более тонких фракций, обогащенных соответствующим набором элементов, прежде всего алюминием и железом. При этом в распределении фактора 2 на акваториях ПА и ЛУ месторождений прослеживается отчетливая тенденция к зональности с ростом значений с запада на восток, в сторону открытого моря (рис. 6). В пространственной изменчивости фактора 1 выражена в большей степени широтная зональность с увеличением значений с севера на юг на полигоне ПА и в центрально-южной области полигона ЛУ.
Таблица 4
Нагрузки факторов, выявленных методом главных компонент с последующим вращением
квартимакс при анализе концентраций кислоторастворимых форм металлов в донных отложениях исследованных акваторий (n = 264)
Table 4
Loadings of factors revealed by principal component method and quartimax axis rotation for the concentrations of metal acid-soluble forms in bottom sediments of the areas studied (n = 264)
Элемент Фактор ( 1 r1 2 = 0,581) 2
Al 0,249 0,923
As 0,076 0,822
Ba 0,846 0,286
Cd 0,806 0,460
Cr 0,307 0,811
Cu 0,951 0,262
Fe 0,476 0,775
Pb 0,952 0,084
Hg 0,904 0,253
Zn 0,920 0,335
Объясняемая дисперсия 5,227 3,331
Доля объясняемой дисперсии 0,523 0,333
Примечание. г1-2 — коэффициент корреляции первичных кластеров переменных; нагрузки > 0,7 выделены жирным шрифтом.
ю о
о ф
4 5 6 7
Число собственных значений
Рис. 3. Результаты теста "каменной осыпи" при факторном анализе концентраций кислото-растворимых форм металлов в донных отложениях исследованных акваторий шельфа о. Сахалин Fig. 3. Results of " scree-plot" test for factor analysis of the concentrations of metal acid-soluble forms in bottom sediments of the inshore water areas of Sakhalin Island
6 -
5 -
4 -
3 -
2 -
1 -
0 -
-1 --2
ПАМ МШР ЛМ Кайгон СВШ Поронайск зал. Анива Холмск
Полигон г(1-2)
ПАМ -0,681
МШР -0,018
ЛМ -0,105
Кайгон 0,089
Поронайск 0,446
зал. Анива -0,434
Холмск -0,917
□ Фактор 1 И Фактор 2
Д,
I
п
Рис. 4. Средние значения факторов, полученные при анализе концентраций кисло-торастворимых форм металлов в донных отложениях прибрежных акваторий о. Сахалин. Вертикальные линии — стандартная ошибка, r12 — коэффициент корреляции выделяемых факторов, СВШ — северо-восточный шельф о. Сахалин в целом
Fig. 4. Averaged scores of factors 1 and 2 revealed for the concentrations of metal acid-soluble forms in bottom sediments of the inshore water areas of Sakhalin Island. Whiskers — standard error, r12 — correlation coefficient for factors revealed
На эту зональность накладывается некоторая пятнистость, более заметная у фактора 2, что, очевидно, объясняется неравномерным накоплением тонких фракций по площадям указанных акваторий, которое, в свою очередь, может быть связано с особенностями рельефа дна ПА и ЛУ полигонов, характеризую-
Таблица 5
Коэффициенты корреляции Пирсона некоторых параметров среды и факторов, полученных при анализе концентраций металлов в донных отложениях района работ
Table 5
Pearson's correlation coefficients for some environmental parameters and factors revealed for the concentrations of metal acid-soluble forms in bottom sediments of the region studied
Район Фактор Глубина С орг НУ ЕН pH Алевриты (0,1-0,01 мм) Пелиты (<0,01 мм) Алевропелиты (<0,1 мм) Среднее значение модуля г ± SE
Все данные, п = 264 1 2 -0,151 0,090 0,527 0,536 0,786 0,044 -0,065 -0,253 -0,198 -0,071 0,303 0,598 0,267 0,630 0,313 0,643 0,326 + 0,082 0,358 + 0,089
Пильтун-Астохское месторождение, п = 57 1 2 -0,101 0,415 -0,217 0,412 0,087 -0,128 -0,096 0,122 -0,138 0,111 0,215 0,177 -0,040 0,062 0,183 0,172 0,135 + 0,023 0,200 + 0,048
Полигон МШР, п = 111 1 2 0,073 0,313 0,216 0,365 0,393 0,338 -0,075 -0,222 -0,095 -0,462 0,347 0,368 0,502 0,406 0,396 0,400 0,262 + 0,060 0,359 + 0,025
Луньское месторождение, п = 43 1 2 0,108 0,032 -0,171 0,516 0,103 0,219 0,150 -0,349 -0,040 0,016 0,081 -0,055 -0,092 0,173 0,059 -0,019 0,101 ±0,015 0,172 ±0,065
Кайгон, п = 10 1 2 -0,003 -0,304 0,529 -0,207 -0,149 -0,038 0,602 -0,174 -0,050 0,385 0,749 0,064 0,357 -0,334 0,671 -0,059 0,389 ±0,103 0,195 ±0,048
СВШ, п = 221 1 2 0,057 0,352 0,179 0,421 0,282 0,259 -0,117 -0,139 -0,098 -0,277 0,322 0,311 0,377 0,384 0,361 0,352 0,224 ±0,045 0,312 ±0,031
Поронайск, п = 10 1 2 0,399 -0,294 0,407 0,822 0,694 0,724 -0,291 -0,330 -0,552 -0,196 0,363 0,531 0,224 0,509 0,356 0,545 0,411 ±0,053 0,494 ± 0,076
Зал. Анива, п = 26 1 2 0,298 0,438 0,381 0,353 0,258 -0,160 0,148 -0,027 0,123 0,261 0,575 0,039 0,518 0,226 0,612 0,109 0,364 ± 0,067 0,202 ± 0,052
Холмск, п = 7 1 2 -0,670 0,830 0,367 -0,207 0,657 -0,495 -0,090 -0,136 -0,915 0,846 0,540 -0,272 0,625 -0,389 0,566 -0,301 0,554 ± 0,086 0,435 ± 0,096
Примечание. МШР — полигон меридионального и широтного разрезов, СВШ — северо-восточный шельф о. Сахалин в целом; г — коэффициент корреляции Пирсона, БЕ — стандартная ошибка; жирным шрифтом выделены статистически значимые г (р < 0,05); выделяемые факторы как в табл. 4.
-0.192 • -0.189 •
• -0.188
• -0.197
-0.181
• -0169 -0.168
a -0.167 • -0.107 • -0.122
1
^4.86
2.99 •
0£3 0^54 1.85 1*9° / i 1 \ \ 1^33 '
д !
073 Ц6
2.09
к1.03
1.68
-2.66
Рис. 5. Распределение значений факторов 1 (а, в, д) и 2 (б, г, е), полученных при анализе концентраций кислоторастворимых форм металлов в донных отложениях прибрежных акваторий о. Сахалин, в порту Кайгон (а, б), Поронайск (в, г) и Холмск (д, е) Fig. 5. Distribution of the scores of factor 1 (а, в, д) and 2 (б, г, е) revealed for the concentrations of metal acid-soluble forms in bottom sediments of the inshore water areas of Sakhalin Island in ports of Kaigon (а, б), Poronajsk (в, г) and Kholmsk (д, е)
щегося развитием холмов и ложбин между ними. Известно, что зерна тонких фракций скапливаются в ложбинах с несколько более затишными гидродинамическими условиями (Современное осадкообразование ..., 1997). Рост содержания тонких фракций из района ПА к ЛУ месторождению вызывает, очевидно, соответствующий тренд значений фактора 2 (рис. 7). Кроме концентрации НУ у фактора 1 и содержания тонких фракций у фактора 2, выделенные факторы мало зависят от других параметров донных отложений, за исключением Сорг, вариации которого объясняют примерно равное количество дисперсии обоих факторов (соответственно 27,8 и 28,8 %, табл. 5).
Рис. 6. Распределение значений факторов 1 (а, в) и 2 (б, г), полученных при анализе концентраций кислоторастворимых форм металлов в донных отложениях прибрежных акваторий о. Сахалин, в районах Пильтун-Астохского (а, б) и Луньского (в, г) месторождений
Fig. 6. Distribution of the scores of factor 1 (а, в) and 2 (б, г) revealed for the concentrations of metal acid-soluble forms in bottom sediments of the inshore water areas of Sakhalin Island at Piltun-Astokhskoe (а, б) and Lunskoe (в, г) deposits
Такая "стройная" картина корреляций выделенных факторов с параметрами среды наблюдается для полной выборки их значений и заметно нарушается на отдельных полигонах (табл. 5). Так, на полигоне ПА статистически значимые корреляции обнаружены только у одного из факторов — фактора 2 — с глубиной и Сорг; в зал. Анива и порту Холмск с долей тонких фракций положительно коррелирует первый, а не второй фактор (в Холмске эта связь статистически не значима). В порту Поронайск сила связи фактора 2 с содержанием НУ даже выше, чем у фактора 1; в порту Холмск отмечена сильная корреляция обоих факторов с рН, отрицательная у первого и положительная — у второго.
1,0-1 т 1,5
1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35
Номер станции
Рис. 7. Распределение значений факторов, полученных при анализе концентраций кислоторастворимых форм металлов в донных отложениях прибрежных акваторий о. Сахалин, вдоль меридионального разреза
Fig. 7. Distribution of the factor scores revealed for the concentrations of metal acid-soluble forms in bottom sediments of the inshore water areas of Sakhalin Island along the meridian section
По определению, в результате процедуры факторного анализа получаются новые переменные, которые не коррелируют между собой (ортогональны друг другу). Естественно, что и в нашем случае, несмотря на заметную связь первичных кластеров переменных, в полной выборке значений факторов отсутствует сколь бы то ни было существенная корреляция (r ~ 0,000). Однако для отдельных акваторий правило "ортогональности" явно нарушено.
В большинстве незагрязненных районов северо-восточного шельфа, за исключением ПА месторождения, связь факторов мала и статистически незначима (r = -0,105 ■ 0,089). Для грунтов порта Холмск, находящихся под сильным антропогенным прессом, отчетливо выражена их отрицательная связь (r = -0,917). В целом уровень их пропорциональности, прямой или обратной, возрастает по мере роста степени загрязнения среды. Кроме того, в среднем увеличивается и связь факторов с физико-химическими параметрами осадков (табл. 5).
На акватории порта Поронайск, где основное техногенное воздействие на морскую среду происходило в 80-е гг. прошлого века и к моменту отбора проб прошло почти десятилетие с тех пор, как оно существенно ослабело, обнаружена заметная положительная, хотя и статистически незначимая, вероятно из-за малого объема выборки, корреляция факторов (r = 0,446). При этом отмечено совпадение одного из самых высоких значений фактора 1 и максимальной величины фактора 2 на выходе из порта (см. рис. 5).
Очевидно, что в порту Поронайск ко времени исследований произошла переработка части "антропогенного материала" за счет механической дифференциации, что привело к его перераспределению и накоплению в типичных областях аккумуляции с относительно высоким содержанием в грунте частиц тонких фракций. О "застарелом" загрязнении порта говорит суммарное количество найденных здесь полиароматических углеводородов, составляющее 5-10 % от общего содержания НУ, что не характерно для свежих нефтяных выбросов (ДВНИГМИ ..., 2002). Такое повышение относительного содержания полиароматических углево-
дородов может наблюдаться на старых нефтяных разливах и при проливах тяжелых нефтяных фракций, например мазута, когда с течением времени происходит биодеградация наименее стойких алкановых и нафтеновых фракций нефти и накопление устойчивых полиядерных соединений — хризена, бензо(а)антрацена, бен-зо(е)пирена и т.п.
Следует обратить внимание на заметную и статистически значимую обратную связь выделяемых факторов на акватории ПА месторождения (r = -0,681). По аналогии с портом Холмск, где отрицательная корреляция еще выше, а загрязненность грунтов максимальна, возможно, что здесь начинают проявляться результаты интенсификации хозяйственной деятельности, которые еще не привели к какому-либо существенному обогащению грунтов поллютантами, но уже отражаются в структуре взаимосвязей отдельных элементов. Как и для всей выборки, измеренные в районе ПА концентрации наиболее логично представить в виде двух факторов, один из которых соответствует первому фактору для полной выборки с высокими нагрузками на медь, свинец и кадмий, а другой — второму (алюминий, железо, хром).
Заключение
Применение факторного анализа к концентрациям подвижных форм металлов, измеренным в осадках района работ, позволило получить интегральные характеристики как степени антропогенного воздействия на морскую среду, так и выраженности естественных процессов трансформации металлов (механической дифференциации). В 2001 г. антропогенное обогащение грунтов металлами было минимальным на исследованной акватории северо-восточного шельфа о. Сахалин и максимальным в порту Холмск; в порту Поронайск степень техногенного воздействия была несколько выше, чем в зал. Анива. Процессы механической дифференциации были выражены в наибольшей степени в зал. Анива, где был опробован наиболее широкий диапазон глубин, а состав грунтов был максимально разнообразным.
Список литературы
Афифи А., Эйзен К. Статистический анализ. Подход с использованием ЭВМ. — М.: Мир, 1982. — 488 с.
Белан Т.А., Олейник Е.В. Состав, распределение и современное состояние бентоса на Пильтун-Астохской нефтегазоносной площади // Тем. вып. ДВНИГМИ. № 2. — Владивосток: Дальнаука, 2000. — С. 166-177.
Боровиков В.П., Боровиков И.П. STATISTICA® — Статистический анализ и обработка данных в среде Windows®. — М.: Информационно-издательский дом "Фи-линъ", 1998. — 608 с.
ДВНИГМИ. Отчет о результатах фоновых морских экологических исследований летом 2001 г. в рамках проекта "Сахалин II". — Владивосток: ДВНИГМИ, 2002. — 256 с.
Ким Дж., Клекка У.Р., Мьюлер К.У. Факторный, дискриминантный и кластерный анализ. — М.: Финансы и статистика, 1989. — 312 с.
Кулаичев А.П. STADIA 6.0. Методы анализа и средства анализа данных в среде Windows. — М.: НПО "Информатика и компьютеры", 1996. — 256 с.
Лишавская Т.С., Мощенко А.В., Чернова А.С. Некоторые органические поллю-танты в донных отложениях прибрежных акваторий острова Сахалин // Изв. ТИНРО. — 2005. — Т. 142. — С. 296-309.
Методика определения содержания нефтяных углеводородов в объектах окружающей среды методом газовой хроматографии с пламенно-ионизационным детектором. — М.: НПО "Тайфун", 2000. — 45 с.
Мощенко А.В., Коновалова Т.В., Белан Т.А., Христофорова Н.К. Изменение биотопических условий и вариации распределения бентоса вблизи платформы Мо-ликпак (шельф северо-восточного Сахалина) // Изв. ТИНРО. — 2005. — T. 142. — С. 223-245.
Носов В.Н. Метод главных компонент и факторный анализ // Компьютерная биометрика. — М.: МГУ, 1990. — С. 162-187.
Олейник Е.В. Видовой состав и количественные характеристики зоопланктона и ихтиопланктона в летних и осенних сборах в районе северо-восточного прибрежья о. Сахалин // Тем. вып. ДВНИГМИ. № 2. — Владивосток: Дальнаука, 1999. — С. 161-166.
Современное осадкообразование в окраинных морях востока Азии (статистические модели). — Владивосток: Дальнаука, 1997. — 302 с.
Ткалин А.В., Ройл Д.Дж., Сергушева О.О. Некоторые результаты экологического мониторинга на Пильтун-Астохской площади в 1998 году // Тем. вып. ДВНИГМИ. № 2. — Владивосток: Дальнаука, 1999. — С. 145-160.
Шулькин В.М. Металлы в экосистемах морских мелководий. — Владивосток: Дальнаука, 2004. — 279 с.
Bradford W.L., Luoma S.N. Some perspectives on heavy metal concentrations in shellfish and sediment in San Francisco Bay, California // Contaminants and Sediments / Ed. R.A. Baker. — Ann Arbor: Ann Arbor Scientific Publishing, 1980. — P. 501-532.
CSA. Environmental survey report for the Piltun-Astokhskoye field offshore Sakhalin Island, Russia (the results of 1995 observations): Continental Shelf Associates. — 1996. — November. — 121 p.
CSA. Environmental survey report for the Piltun-Astokhskoye field 1996 offshore Sakhalin Island, Russia: Continental Shelf Associates. — 1997. — June. — 76 p.
Daskalakis K.D., O'Connor T.P. Distribution of chemical concentrations in US coastal and estuarine sediment // Mar. Environ. Res. — 1995. — Vol. 40. — P. 381-398.
Lishavskaya T.S., Moshchenko A.V., Konovalova T.V., Chernova A.S. Changes of the petroleum hydrocarbon content in marine environment near the operating oil platform MOLIKPAQ (Okhotsk Sea, northeast shelf of Sakhalin) // Pacific Oceanography. — 2005. — Vol. 3, № 1. — P. 49-54.
Long E.R., Macdonald D.D., Smith S.L., Calder F.D. Incidence of adverse biological effects within ranges of chemical concentrations in marine and estuarine sediments // Environment Management. — 1995. — Vol. 19. — P. 81-97.
StatSoft, Inc. Электронный учебник по статистике. — М.: StatSoft, 2001. http://www.statsoft.ru/home/ textbook/default.htm.
Summers J.K., Wade T.L., Engle V.D., Malaeb Z.A. Normalization of metal concentrations in estuarine sediments from the Gulf of Mexico // Estuaries. — 1996. — Vol. 19. — P. 581-594.
Surfer Mapping System. — Copyright © 1993-2000. http:// www.goldensoftware.
com.
Tkalin A.V. Background pollution characteristics of the NE Sakhalin Island shelf // Mar. Pollut. Bull. — 1993. — Vol. 26. — P. 704-705.
Tkalin A.V., Belan T.A. Background ecological condition of the NE Sakhalin Island shelf // Ocean Res. — 1993. — Vol. 15. — P. 169-176.
Trefry J., Presley B.J. Heavy metals in sediments from San Antonio Bay and the northern Gulf of Mexico // Environ. Geol. — 1976. — Vol. 1. — P. 283-294.
UNEP: Determination of petroleum hydrocarbons in sediments. Reference Methods for Marine Pollution Studies. — 1992. — № 20. — 75 p.
UNEP: Manual for the geochemical analyses of marine sediments and suspended particulate matter. Reference Methods for Marine Pollution Studies. — 1995. — № 63. — 74 p.
Windom H.L., Schropp S.J., Calder F.D. et al. Natural trace metal concentrations in estuarine and coastal marine sediments of the southeastern United States // Environ. Sci. Technol. — 1989. — Vol. 23. — P. 314-320.
Поступила в редакцию 11.12.07 г.